CC BY
26
УДК 621:531.3
Теоретическое исследование виброакустических характеристик при динамических испытаниях на циклическую прочность
С. А. Шамшура, И. В. Богуславский, А. Н. Чукарин
(Донской государственный технический университет)
Приведены результаты теоретических исследований возбуждения вибраций и излучения шума оборудования для динамических испытаний лонжеронов вертолётов. Получены зависимости для определения спектров вибраций в рабочей зоне операторов.
Ключевые слова: виброакустические характеристики, динамические испытания.
Введение. Вопросы вибро- и шумозащиты всегда были актуальны в машиностроительном производстве. В настоящее время это связано с усилением контроля соблюдения санитарных норм и правил [1] на рабочих местах, что в конечном итоге повышает производительность труда.
Измерения и анализ уровней шума и вибраций согласно [2], проведённые сотрудниками университета на участке динамических испытаний лонжеронов лопастей вертолётов на циклическую прочность ОАО «Роствертол», показали, что проведение мероприятий по шумо- и виброзащите позволит существенно улучшить условия труда. В ходе экспериментов было установлено, что стенд для динамических испытаний лонжеронов лопастей вертолёта выступает и как источник шума, генерируемого колебаниями протяжённого лонжерона, вибровозбудителя и тросов натяжения, и как источник вибраций, возникающих главным образом в результате колебаний опор стенда и передаваемых в цеховое пространство через их основание.
Результаты исследований. Рассмотренная в данной работе математическая модель стенда динамических испытаний лонжеронов лопастей вертолёта является базовой для дальнейших работ по расчёту средств виброзащиты в лаборатории динамических испытаний (ЛДИ) ОАО «Роствертол».
Стенд для динамических испытаний лонжеронов лопастей вертолёта представляет собой сложную конструкцию, состоящую из следующих подсистем (рис. 1):
• подсистема опоры со стороны лонжерона 1;
• подсистема опоры со стороны системы натяжения 2;
• подсистема тросов натяжения 3)
• подсистема лонжерона 4;
• подсистема вибровозбудителя 5.
Рис. 1. Структурная схема стенда динамических испытаний лопасти вертолёта
Рассмотрим более подробно структуру стенда. Опоры со стороны лонжерона і и со стороны системы натяжения 2, массами М1 и М2 соответственно, соединены с полом цеха упруго-
1857
диссипативными связями. Жёсткость с9 и диссипация /7. связей (/ =1,2, ] =1,2) в общем случае нелинейно зависят от смещений х1 и скоростей х1.
Подсистема натяжения 3 представляет собой систему стальных тросов, с помощью которой осуществляется предварительный натяг лонжерона в стенде. С точки зрения излучения звуковых колебаний эта подсистема представляет собой набор струн. Однако для задачи исследования передачи вибрации в опоры стенда её можно заменить упруго-диссипативной связью, как это сделано на структурной схеме.
Возбудитель колебаний представляет собой два эксцентрика, закреплённых на правом конце лонжерона, и для уменьшения колебаний в продольном направлении вращающихся навстречу друг другу. Вращение производится от электродвигателя постоянного тока (на схеме не показан) с приводом через упругую муфту, практически исключающую передачу колебаний на двигатель.
Эксцентрики имеют форму полуцилиндров, сумма несбалансированных масс эксцентриков равна т, расстояние от несбалансированной массы до оси вращения — г, циклическая частота
вращения со может регулироваться. Общая масса эксцентриков с системой закрепления равна М3.
Лонжерон лопасти вертолёта 4 имеет каплевидное сечение (рис. 2). Конструкция стенда позволяет рассматривать лонжерон как балку, шарнирно закреплённую с одного конца. Известно, что такая балка будет иметь п собственных частот и её математическая модель может быть аппроксимирована в виде системы с п сосредоточенными массами и соответствующими частотами и коэффициентами затухания.
Учтены также особенности закрепления опор, которые позволяют представить их как массивные призмы, шарнирно закреплённые одним углом и связанные с основанием упругодиссипативной связью.
Возможность моделирования подсистемы лонжерона 4 (рис. 3) в виде гармонического осциллятора [3] дополнительно можно обосновать резонансным режимом работы стенда динамических испытаний, обеспечиваемым настройкой силы натяжения тросов согласно [4] и подстройкой частоты вращения вибровозбудителя для получения заданной амплитуды колебаний.
Рис. 2. Сечение лонжерона лопасти вертолёта
1858
Воспользуемся формализмом Лагранжа 2-го рода [5] для получения уравнений динамики. Кинетическая энергия Т, потенциальная энергия Р, диссипативная функция Я и обобщённая сила , действующая на систему, имеют следующий вид:
Г = 2
РЛ
2
+ Мз (*з + Уз)2
с, (о,/, )2 + с2 (а,/, )2 +с3( ^ + *3 )2 + у2 - ^ - а2/21 + с4 Г ^ - *3 )2 + у2 -/,-а^
(1)
Я = і 2
(ОіА )2 + К (<Уі )2 + (>1*!+у! - <У2) + Ь (№+у! - оА)
С2Уз = /г(а))зіп(а)ґ).
Амплитуда возбуждающей силы зависит от частоты вращения эксцентриков следующим образом [6]:
^ (со) = пю2г. (2)
Вычисляя частные производные по координатам и импульсам [7] и подставляя их в уравнение Лагранжа, получаем систему уравнений динамики моделируемой системы:
¿А + ГіАЇ ~ (л/*з+Уз - ¿Л) /2 + - с4 ^ - х3 )2 +/32 - /, - ах/2) /2 = 0;
72а2 + /72а2/2 - /73 (7^? + /з2 - «Уг) К + с2а2/2 - с3 (7(^+>гз)2+/з2 - ¡-2 ~ °2>2) >2 = 0;
Мз*з Мз
Г • , ^ 1 2 2 <У2 ^
+ /74Л-3 1 + С3 (*3 + ¿2 )
1 №+УІ; V №+УІ;
1-
¿2 + а2/2
л/(^ +Хз)2 +/з
с4{х3-^)
1-
/1+о1/2
^ -^з)2 +/з2
= 0;
(3)
/
Мз/з лз/з
Г . , А
^ _ 2 2
V а/*з + Уз ,
(
(
+ /74/з
аЛ
л
V л/^з + Уз ,
+ сз/з
1-
¿2 + а2/2
+ С4/з
Л
1-
/1+о1/2
^-*з)2 + Уз
= тш2г віпші.
Полученная модель содержит нелинейные составляющие и может быть линеаризована [8] в окрестности нулевых значений л-3, у2, а! и а2. Элементы матрицы частных производных
>4 = ) имеют вид (для простоты вычислений не учтён множитель -1):
5/;
5/;
5а^
^ = -1/2 С^("2/1+2хз)/:
5Л",
— = с I2 + с I2' ___— = 0 1
11 4 2 ' ^ '
5а,
5/;
5^ 5^ = М2+л4/22; О II
С4/3/2 Л4Х3/2
/74/3/2
А2 -2^з +^з2 +/з2 ' ^ - 21^Х3 +Х2+у2 ' ЗЛ-3 + /з2 # + /з2 #
— = 0;
5а,
-2- = с2/2 + с3/2 ; = 0;
5а, 2 1 3 2 5а,
^ = /72/2+/7з/22;
5а,
1859
5л-,
^3 ( ^2 + Х3 ) ^2 .
д?з
5а,
+ х3) +у.
сЛхз~к)1:
Ж
4^ - 2/1Х3 +^з2 +/з2
5/3
5а,
СзУз4_______ . ^2 _
+ ^з)2+/з2 5Хз С3 {1~2 + Хз) І2 . <5/д
/73X3/2
л/*з + /з2
5/,
' 5/3
2 5^
3
/74 Х3 /2
^з + /з2
5/3
5а,
^зУ 3^2 .
4*1+/з2 #
Ь3х312
4*1+/з2
5/3
5л-,
= С,
1-
¿2 + а2/2
^з{^-2+Хз) {1-2 + ^2)
^{¡-2+Хз)2 +УІ ) ((¿2+Хз)2+/з2)
Сз(^+^з)2 (^1+аЛ).
+ СЛ
1-
/1+а1/2
дМ +^з)2 +/:
V 3/2
((^ + ^з )2 + Хз )
5/3 _ с3 (¿2 + а2/2)у3 с4(/1+Хз)(/1+а1/2)/з
5/3 ((¿2+^з)2+/з2)3/2
5/3
5л-,
= /7,
1-
а2/2
л/хз + /з
^з*з®2^2
И +/з2)
3/2 + Л4
((*1 +^з)2 +/з2)
^ \ <у2
1-
4*1 + /з2
/?4 -^~3 01 /2
/ * 2 * 2 (*з +Уз)
5/3
5/.
^зОАУ
3л3и2'2/3 + ^4*З^АУЗ
3 И+/з2)3/2 (*з+Уз)312
5/4 ^ лУ 3^2 5/^ ^ С3У3І2 . 5/4 _ ^уА . 5/4 _ ЬзУзК
-2^X3 +ХІ + УІ ' д02 ^ +Х3)2 +/з? А 4*3 + /з2 д62 4*3 + /з"
5/4 _ с3 Уз{Ц. + ^2^){Ц. + Хз) ! С3 М^і+аЛК*!-* з)
дХз ((¿2 + ^3 )2 + Хз )3 2 ^+хз)2 +УҐ)'2
\ ( \
5/4 _г 1 ¿2 + а2/2 Сз/3 (¿2 + °2^2 ) , 1 + ^2 СзУз^+аЛ) .
^3 о/з ^(¿2 +^3)2 +/З2у / 2 \3/2 4 ((¿2+*з) + Хз2) ^-^з)2 +/з у ({к - хз)2 + У2з)3'2
_ ^3X3 а2/2у3 | ЬцХ^АУз
дх.
3 {*з+Уз)312 (*1+У2з)
ч 3/2
5/4
5/3
= /7,
1-
а2/2
7Х3 + /з
^зУз®2І;
3/2 + Л4
1-
аЛ
7хз + /з
3/2
(*3 +/з)
Н+/з2)
После подстановки значений координат в точке равновесия X*
Х*=|а* а* а* а* л-3* /* л-3* /3*) = {0 0 0 0 0 0 0 0}
получаем матрицу частных производных, соответствующую собственной матрице линеаризованной системы в явном виде:
(4)
1860
Е.
5х
С,/,
0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
1 +сЛ Л 0 М2+М2 Л 0 с4/2 Л 0 Л4/2 Л Л4/2 Л
0 С2/1 + с3/2 0 0 /73/2 /73/2
Л Л Л Л Л
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1
с4/і ^зА Л4/2 /73/2 С3 +С4 0 /73 +/74 0
М3 Мз Мз Мз Мз Я73
0 0 Л4/2 /73/2 0 0 0 /73+/74
Мз Мз
(5)
(6)
и система уравнений динамики может быть записана в виде:
\Х = АХ + Ви)
V = СХ.
Экспериментальные измерения, проведённые на участке динамических испытаний, показали, что амплитуды вибраций во всех частотных диапазонах на опорах стенда динамических испытаний на 15...20 дБ меньше амплитуд колебаний на конце лонжерона, закреплённого на вибровозбудителе. То есть амплитуды колебаний опор на несколько порядков меньше амплитуды колебаний вибровозбудителя. Исходя из этого, можно осуществить дальнейшее упрощение модели (3) стенда динамических испытаний, рассматривая его с точки зрения анализа формы возникающих колебаний как упругую подсистему, прикреплённую к неподвижным опорам (рис. 4, а), а с точки зрения передачи вибрации через опоры — как упругие подсистемы опор, на которые действует внешняя сила, характеристики которой определяются из предыдущей подсистемы (рис. 4, б).
«а, £3 ,—3—,
Уз
а)
Х3
£
•2
/
/
Рис. 4. Декомпозиция динамической модели рис. 3 на подсистемы: а — вибровозбудителя; б — опор
Таблица 1
Режимы динамических испытаний лонжеронов
Испытываемый участок Статическая нагрузка Переменная нагрузка аисп, МПа Частота колебаний образца, Гц
Средний участок с двумя секциями каркаса и в плоскости тяги 0р =88,2 МПа {(?= 1176 Н) 54 68,6 88 41 42 42
Концевой участок с двумя секциями каркаса в плоскости тяги ар = 58,8 МПа {(?= 1176 Н) 35 54 68,6 88 38 38 40 42
Концевой участок с двумя секциями каркаса в плоскости тяги ар = 58,8 МПа 54 68,6 88 38 40 42
1861
(7)
(Ю)
Соответствующие системы уравнений, описывающие динамические свойства полученных подсистем, представлены ниже. Для подсистемы вибровозбудителя:
Í М3х3 + (с3 + с4 )л"3 + (h3 + h4}x3 + £1у3 — 0,
|М3у3 + (h3 + Л4) Уз = ты2 г sin сot.
Для подсистем опор:
J& + h1l2á1 +CJ& =/2/r1(f),. (8)
_72q2 + /íj/j o2 +c2¡i a2 — l2F2 (t). (9)
Соответственно силы F1 (t) и F2 (t) получаются в линейном приближении из стационарных решений системы уравнений (7) по формулам:
И(0 = сЛ(0 + /7Л(0'
lF2(f) = c3*3(f) + h3x3(t).
Численный расчёт удобно проводить, представив уравнения (7) — (9) в форме Коши:
Z1=Z2'
z2 =(-(сз +C4)zl -(/73 +/74 )z2 +ElZ3)/M3;
Ч =z4;
z4 =(-(/73 +h4)z4 + ГГКх)2 Г Sin (OÍ )/М3 1 где z,=x3, z2=x3, z3=y3, Z4=y3)
í¿5 =Zer
U = (-МЧ - crfzs + (0) / Л;
í¿7 = z6;
[z8 = (~h2l2ze - c2I2z7 + l2F2 (t))/j2, где f z7 -a2, z8-á2.
(11)
(12)
(13)
Таблица 2
Параметры модели стенда динамических испытаний
Обозначение Значение Размерность Пояснения
Сі 5-10b Н/м Жёсткость опоры со стороны лонжерона
Сі 10у H/m Жёсткость конструкции со стороны тросов натяжения
Сз 2,2-104 H/m Приведённая жёсткость системы тросов натяжения
й 3,5-10b H/m Приведённая жёсткость лонжерона
Лі 1000 H-c/m Диссипация в опоре со стороны лонжерона
1000 H-c/m Диссипация в опоре со стороны тросов натяжения
fb 10 H-c/m Приведённая диссипация в системе тросов натяжения
Л, 10 H-c/m Приведённая диссипация в лонжероне
/і 0,5 m Длина образующей опоры
h 0,4 M Высота образующей опоры
Li 2 M Длина лонжерона
Li 1,5 M Длина системы тросов
Ji 1 H-m¿ Момент инерции опоры со стороны лонжерона
Ji 1 H-m¿ Момент инерции опоры со стороны тросов натяжения
Mi 10 КГ Масса подсистемы возбудителя
m 0,4 КГ Масса груза вибровозбудителя
r 0,05 M Радиус эксцентриситета груза вибровозбудителя
1862
а)
6)
в)
Рис. 5. Результаты моделирования работы стенда динамических испытаний лонжеронов лопастей вертолётов: а — смещение в горизонтальном направлении; б — смещение в вертикальном направлении; в — переменная составляющая растягивающей силы, действующей на лонжерон
1863
Таким образом, определение коэффициента передачи вибрации через опоры стенда динамических испытаний производится в два этапа: на первом этапе определяется установившееся решение системы x(f) x(f) (7), а на втором этапе полученные установившиеся решения подставляются в выражение для внешних сил уравнений (8) и (9), для которых в свою очередь так же находятся установившиеся решения (t), (t), о2 (t), á2 (t).
Моделирование динамики стенда циклических испытаний лонжерона лопасти вертолёта. В соответствии с программой испытаний лопасти несущего винта определены условия и режимы испытания лонжеронов лопастей несущего винта вертолётов. Для рассматриваемого в данной работе стенда динамических испытаний, на котором испытываются средние и концевые участки лонжерона, определены следующие условия испытаний (табл. 1).
Результаты предварительных измерений и анализа конструкции стенда динамических испытаний позволили оценить значения постоянных с i модели стенда. Они приведены в табл. 2.
Результаты моделирования, выполненные на основе [9], приведены на рис. 5. Как показало моделирование, при быстром выводе оборотов двигателя вибровозбудителя в зону резонанса возможны колебания амплитуды упругих сил, действующих на лонжерон, что может привести к его повреждению. По этой причине при запуске стенда динамических испытаний применяется плавный выход в зону резонанса, с относительно медленным повышением частоты вращения вибровозбудителя (рис. 6).
Зависимость возбуждающей силы в этом случае может быть описана формулой
F(t,со) = /n(co(l-ef/r)) rsin(co(l-ef/r)f), (14)
учитывающей постепенное возрастание как частоты, так и амплитуды колебаний возбуждающей силы. Зависимость возбуждающей силы от времени приведена на рис. 7, а, а соответствующие графики изменения координат состояния моделируемой системы — на рис. 7, б, в.
Заключение. Приведённые в работе результаты математического моделирования позволяют сделать вывод о работоспособности разработанных моделей стенда динамических испытаний лонжеронов лопастей вертолётов и их применимости для расчёта характеристик источников шума и вибрации. Результаты моделирования показали соответствие экспериментальным измерениям как по качественному поведению системы, так и по количественным показателям — амплитудам колебаний и величинам сил.
1864
0 1 2 3 4 5
Время, с
а)
^ -I-----I-------I-------I-------I------I-------I-------I-------Г
0 1 2 3 4 5
Время,с
6)
Рис. 7. Результаты моделирования работы стенда динамических испытаний лонжеронов лопастей вертолётов при плавном увеличении частоты вращения: а — смещение в горизонтальном направлении; б — смещение в вертикальном направлении; в — переменная составляющая растягивающей силы, действующей на лонжерон
1865
Библиографический список
1. СН 2.2.4/2.1.8562-96. Шум на рабочих местах, в помещениях жилых, общественных зданий и на территории жилой застройки. — Москва, 1996.
2. ГОСТ ССБТ 12.1.043-90. Вибрация. Методы измерений на рабочих местах в производственных помещениях. — Москва: Издательство стандартов, 2001.
3. Рабинович, М. И. Введение в теорию колебаний и волн / М. И. Рабинович, Д. И. Трубец-ков. — Ижевск: НИЦ «Регулярная и хаотическая динамика», 2000. — 560 с.
4. Программа испытаний лопасти несущего винта. — Москва: МВЗ им. Миля, 2002. — 25 с.
5. Лич, Дж. К. Классическая механика / Дж. К. Лич. — Москва: Изд-во иностранной литературы, 1961. — 173 с.
6. Яворский, Б. М. Справочник по физике / Б. М. Яворский, А. А. Детлаф. — Москва: Наука, 1973. - 942 с.
7. Тер-Крикоров, А. М. Курс математического анализа: уч. пособие для вузов / А. М. Тер-Крикоров, М. И. Шабунин. — Москва: Наука, 1988. — 816 с.
8. Арнольд, В. А. Математические методы классической механики / В. А. Арнольд. — Москва: Наука, 1974. — 432 с.
9. Потёмкин, В. Г. Инструментальные средства MATLAB 5.x / В. Г. Потёмкин. — Москва: Диалог-МИФИ, 2000. — 336 с.
Материал поступил в редакцию 01.12.2011.
References
1. SN 2.2.4/2.1.8562-96. Shum па rabochix mestax, v pomeshheniyax zhily'x, obshhestvenny'x zdanij i na territorii zhiloj zastrojki. — Moskva, 1996. — In Russian.
2. GOST SSBT 12.1.043-90. Vibraciya. Metody' izmerenij na rabochix mestax v proizvodstvenny'x pomeshheniyax. — Moskva: Izdatel'stvo standartov, 2001. — In Russian.
3. Rabinovich, М. I. Vvedenie v teoriyu kolebanij i voln / М. I. Rabinovich, D. I. Trubeczkov. — Izhevsk: NICz «Regulyarnaya i xaoticheskaya dinamika», 2000. — 560 s. — In Russian.
4. Programma ispy'tanij lopasti nesushhego vinta. — Moskva: MVZ im. Milya, 2002. — 25 s. — In Russian.
5. Lich, Dzh. K. Klassicheskaya mexanika / Dzh. K. Lich. — Moskva: Izd-vo inostrannoj literatury', 1961. — 173 s. — In Russian.
6. Yavorskij, В. M. Spravochnik po fizike / В. M. Yavorskij, A. A. Detlaf. — Moskva: Nauka, 1973. — 942 s. — In Russian.
7. Ter-Krikorov, A. M. Kurs matematicheskogo analiza: uch. posobie dlya vuzov / A. M. Ter-Krikorov, М. I. Shabunin. — Moskva: Nauka, 1988. — 816 s. — In Russian.
8. Arnol'd, V. A. Matematicheskie metody' klassicheskoj mexaniki / V. A. Arnol'd. — Moskva: Nauka, 1974. — 432 s. — In Russian.
9. Potyomkin, V. G. Instrumental'ny'e sredstva MATLAB 5.x / V. G. Potyomkin. — Moskva: Di-alog-MIFI, 2000. — 336 s. — In Russian.
THEORETICAL STUDY ON VIBROACOUSTICS UNDER DYNAMIC TEST FOR CYCLICAL STRENGTH S. A. Shamshura, I. V. Boguslavskiy, A. N. Chukarin
(Don State Technical University)
The theoretical study on the excitation of the equipment vibrations and noise emission under the dynamic testing of helicopter longerons is resulted. The dependences for determining vibration spectra in the operating space are obtained.
Keywords: vibroacoustics, dynamic test.
1866
