УЧЕНЫЕ ЗАПИСКИ ЦАГИ
Том XXIII 1992 . № 3
УДК 629.735.33.015.4 : 533.69.048.5
МЕТОДИКА РАСЧЕТА ДИНАМИЧЕСКОГО НАГРУЖЕНИЯ И МАССЫ КРЫЛА САМОЛЕТА С УЧЕТОМ ВОЗДЕЙСТВИЯ АТМОСФЕРНОЙ ТУРБУЛЕНТНОСТИ
В. М. Литвинов, В. В. Политое
Разработаны аналитические методы расчета перегрузки, изгибающего и крутящего моментов, перерезывающей силы консоли крыла от воздействия непрерывной случайной атмосферной турбулентности.
На основании предложенных методов дан способ оценки массы силовых .элементов крыла по условиям прочности от воздействия атмосферной турбулентности. Проведен сравнительный анализ результатов расчета передаточных функций динамических нагрузок и массы крыла самолета по предложенному и численному методам.
Одной из основных характеристик, существенно влияющих на летнотехнические и экономические показатели летательных аппаратов (ЛА), является масса конструкции. Определяющее влияние на ее величину оказывают внешние динамические нагрузки, действующие на ЛА. Поэтому важно уже на начальных этапах проектирования уметь достаточно недежно определять эти нагрузки с целью прогнозирования массы силовых элементов, обеспечивающей необходимую прочность и жесткость конструкции.
Существующие методы определения динамических нагрузок [1] достаточно громоздки и требуют задания большого количества исходной информации, что приводит практически к невозможности их применения на стадии предварительного проектирования. Задача создания способа определения динамических нагрузок при наличии ограниченного количества исходных данных может быть решена путем разработки упрощенного алгоритма оценки динамического нагружения с помощью методов математической статистики и выделения характерных сечений. Такой подход был использован при решении задачи учета некоторых требований статической и динамической аэроупругости на весовые характеристики авиаконструкций [2, 3].
Предлагаемый аналитический способ расчета динамических нагрузок получен путем сведения системы дифференциальных уравнений, используемых в численных методах с учетом п тонов упругих колебаний, к упрощенной системе уравнений колебательных звеньев с учетом первых низших тонов изгиба и кручения консоли крыла.
1. Общая структура динамического нагружения. Для определения динамического нагружения упругой несущей консоли от воздействия непрерывной случайной атмосферной турбулентности используется статистический подход. Упругая несущая жонсоль рассматривается как стационарная линейная система, на вход которой действует случайная стационарная функция —
вертикальный порыв воздуха — с известной спектральной плотностью 5г((о), которая определяется согласно модели турбулентности Кармана.
Среднеквадратичное значение выходной величины (нагрузки) в линейной системе находится по соотношению
где Гн(/, (о) — частотная характеристика нагрузки (амплитуда и фаза) при действии периодической возмущающей силы; а» — угловая частота.
В качестве периодической возмущающей силы рассматривается воздействие на упругую несущую консоль вертикального синусоидального порыва с частотой со и единичной амплитудой № (/) = е‘“'.
По найденному значению а„ согласно [4, 5] можно определить эксплуатационное значение рассматриваемого вида нагрузки Рэи. В качестве Р3Я понимается перегрузка Дл’, изгибающий ДМ*ЗГ и крутящий ДМ’р моменты, перерезывающая сила Д(?э в каждом сечении крыла от воздействия непрерывной атмосферной турбулентности. Для определения массы конструкции необходимо знать расчетное значение Р* соответствующих динамических нагрузок, которое определяется следующим образом:
где / — коэффициент безопасности, принимаемый согласно требованиям [4].
Следовательно, рассматриваемая задача сводится к разработке алгоритма определения частотных характеристик нагрузок Гн(ш).
2. Алгоритм численного метода. Крыло ЛА схематизируется консольно-закрепленной упругой балкой (рис. 1), которая под действием аэродинамических сил от воздействия вертикального синусоидального порыва может совершать изгибные и крутильные упругие колебания. Принимается допущение, что порыв охватывает всю консоль одновременно. .•
Уравнение движения упругой системы в потоке воздуха запишем в обобщенных координатах с использованием метода Бубнова — Галеркина в следующем виде:
где <7, (/) — обобщенные координаты по времени; п — число тонов упругих колебаний, принятое в расчете.
(1)
(2)
Л
(3)
Рис. 1
Коэффициенты этого уравнения зависят от распределенных упругомассовых характеристик несущей консоли и аэродинамических сил, определяемых здесь по линейной теории скользящего крыла, основанной на гипотезе квазистационарности.
Угол атаки от действия порыва неспокойного воздуха определяем как
аш=ФЦ)/У.
Граничные условия для консольно-закрепленной балки
- К = Ф< = 0 при е = 0;
(4)
\"= = 0 при е= /.
Здесь и ф, — формы |-го тона собственных колебаний консоли.
Каждое уравнение (3) описывает совместные изгибно-крутильные колебания упругой консоли по одному тону.
3. Алгоритм упрощенного расчета. Из опыта расчетов по различным самолетам известно, что основная доля в динамическое нагружение упругих консолей вносится первыми низшими тонами собственных колебаний. Поэтому для получения упрощенного алгоритма ограничимся рассмотрением двухтрех первых изгибных и крутильных тонов. Более высокие тона могут в ряде случаев вносить заметный вклад в нагружение отдельных частей консоли (различных пилонов и узлов крепления подвешиваемых объектов). Однако при проведении оценочных расчетов весовых характеристик консолей этими тонами можно пренебречь из-за малого их влияния на общую массу консоли. Кроме того, на начальном этапе проектирования определить формы собственных колебаний высоких тонов практически невозможно. Любой тон колебаний (пусть это будет, к примеру, первый тон) включает в себя как изгибные у = !I (ф) <71 (0. так и крутильные Ф= ф1 (ф) (/) колебания.
Представляет интерес установившийся процесс колебаний от воздействия единичного порыва № (?) = е'“‘, решение которого можно записать следующим образом:
<7, (*) = («! +к/1)е'“', (5)
где и\, v\ — амплитудные значения и фазовый сдвиг по отношению к воздействующему порыву № (/) изгибных и крутильных колебаний по рассматриваемому (здесь первому) тону.
Таким образом, задача нахождения |7'н(1а))| сводится к определению и, и VI по каждому рассматриваемому тону.
Подставляя (51) в уравнения (3), можно получить выражения для нахождения Ы, и У;. При этом необходимо упростить выражения для коэффициентов уравнений (3) исходя из допущения, что уравнения (3) будут описывать изгибно-крутильные колебания только по одному тону. т. е. каждый тон рассматривается как совокупность изгибного (индекс г) и крутильного (индекс /) тонов колебаний по удовлетворяющим (4) заданным формам. Тогда решения уравнений (3) будем искать в следующем виде:
Я, — ир‘р‘, |
<7,- = ь |
Используя теорему о среднем, предварительно выделив характерные сечения на консоли и пренебрегая отдельными членами, мало влияющими на конечный результат, выражения для коэффициентов уравнений (3) запишем в следующем виде:
С И == т0.7^ср2 ; СЦ — С/7 = — т0,7а^0,7^срЗ '>
С// = т0,?(°2 "Ь х2^0,7^0,7^ср4 ’>
ь /ЛД2<3/ Й • Л —/ 1,875 \4 су Ь
н— \^Г/ 03 СР®’ “— I I ) 0,Злср5»
Ь/, = - с“кр6о,78тх*Ср7; ЬИ= с“ рй^7созхЛср8;
Ьц= - с“ рС^-х^д^^втхЛсрд; ьн = с“р(х0—^д)07^.7СО8ХЛср9;
Р^0,7^срП » ^(/ Р^0,7^“^ -^0^^^ср12»
<*Ц = - С“ • Р&0.7(*о - *д)0,Ар и ;
*11 = СУ. р(*о “ ^0.7 ^0.7 ^-4 Хо_ 16с«к(^-Хд)07)ЛсрИ ’
о — Га
сср15 *
2/ ^0,7$ (*о "^д)о,7 ^с(>16 ’
+ л,.
-Ь
.,(•
ЯС05%
16с^ (^О-^д)07 У
ср I
(7)
где &о,7= 60.7// — относительная хорда консоли в сечении £ = 0,7 (или среднее значение для сечений £ = 0,6; 0,7; 0,8); х0, ха — относительное положение оси жесткости и центра давления консоли в сечении £ = 0,7 (или среднее значение для сечений ё = 0,6; 0,7; 0,8); £/0,з О/0,з— жесткость консоли на изгиб и кручение в сечении £=0,3 (или среднее значение д^я сечейий ё = 0,2; 0,3; 0,4); /По,7 — значение погонной массы консоли в сечении £=0,7 (или среднее значение для сечений £=0,6; 0,7; 0,8); а, х1—относительное положение центра тяжести и относительный радиус инерции профиля консоли в сечении § = 0,7 (или среднее значение для сечений ё=0,6; 0,7; 0,8); Лср1 — статистические поправочные коэффициенты.
Значения поправочных коэффициентов кср1 зависят от характера распределения аэродинамической нагрузки по размаху консоли (циркуляции Г(е)), форм собственных колебаний /,(5) и ф/(е) и влияния осреднения массово-жесткостных и геометрических параметров по размаху. Величины Лср1 для предварительных оценок можно принять равными 1,0. При наличии прототипа консоли величины &ср, можно уточнить на основе данных эксперимента, расчетов по численным методам.
Подставляя (6) и (7) в (3), предварительно сделав упрощающие преобразования со сложными функциями согласно работам [2, 3), получим аналитические зависимости для искомых величин щ и V) в следующем виде:
и‘~ В к, - Ркс
С
С*6 —Я/Ц-О.ЗУ/!*,
1 -
и1
0,5АЬОТ/1-----------------—
V > + »>(*■.-лЛи/
1 + Ю(^ 16(В*2-/?*3-0,З^Л*4)
(8)
где
10£/„ ,
А = 0,5с“ рКСсовх — зіп*); В = —
0,5 в Г „
р___ Ро.з р._ 3 - ___лсоэх_____.
_ Ь1,712 ’ 4 Х° 16С?,(-*о--*д)о.7 ’
£, = -^у + 5(х0 —Хд)27; /г = 40т07(1—о)р2;
/? = 40^р2;
^0,7
т) — сужение консоли; /то7— значение погонного момента инерции консоли относительно оси жесткости в сечении 5 = 0,7 (или. среднее значение для сечений д = 0,6; 0,7; 0,8); А, — поправочные коэффициенты, учитывающие влияние возможного отличия реальных циркуляции и форм собственных колебаний от принятых их аппроксимаций; р — частота в герцах.
В зависимости от рассматриваемых тонов и заданных форм колебаний выражения для и, и и, могут быть уточнены с помощью поправочных коэффициентов кі.
На основании опыта расчетно-экспериментальных исследований частотных характеристик нагружения консолей можно рассматривать взаимодействие следующих тонов:
1. Первый тон изгиба плюс первый тон кручения (и, и и,):
•^^0,7 /. 1 \ ,п\
и! В _ 1 33Р I 1 _ 1 ЧЧ#?™ 1» ' /
1 5(В-1,33£Р) V
с-/г-о,зи(^-^ / 0.5АЬа1(\------‘ -Л
\ 1 + 10(дсо-хд) )
1 16(В-Л-0,3 УА)
2. Второй тон изгиба плюс первый тон кручения (и2 и у,);
'4^*0,7 ✓. 1 \ / * а\
“2 — 11,2735 —1,08 /=■ I ' 5(0,8356 — О.Ов/’О) Г ' ^
\ С-Л-0,ЗМ(хо-жа) /
3. Третий тон изгиба плюс второй тон кручения (ы3 и у2):
0,5АЕЬот , 1 ч , .
и3 о о пси____п и7с:« I * _пя7с;рл'1 1» V11/
и.О/і СО07 , \ ч
23,05В — 0,875/> I ( 5(23,05В —0,875^0) Г
\ С-Л-0,ЗКЛ(іо-ід) /
омь0,(1----------Д... Д
V 1 + ю(х0-ха) )
2 2(В-1,5Я-0,6УЛ)
В соотношениях (9) — (11) обозначения те же, что ив (8).
Соотношения (9) — (11) эквивалентны моделированию каждого тона колебаний консоли идеальным колебательным звеном. При этом на резонансной частоте ррез получается разрыв в графике амплитудно-частотной характеристики нагрузки. Для определения уровня амплитудно-частотной характеристики
нагрузки на резонансной частоте Т (ррез) можно использовать следующие зависимости:
/ ■? і ' \
Т&м,&(){Ррез) — ТьмАо(Рі) /?) <?ДЛ|Д<3] ’
Т (п \_______ Т <п\ і^+^-|)2<?1-]'/2.
Т \ipJP рез) *\у\Рі) ;2 ’
где
Рдм.ло —
,•2 _ .-2 /і+1 /і
(1 -/?)-(•-/?+.)
________Ц+і_________
(2/?7?+1-./?-/?+,)
1/2
1/2 •
Здесь
4(М, 4(? 1
/< Р‘/Ррез’ /і+1 Рі+і/Ррез>
(рі), Гд,(р,) — значения амплитудно-частотной характеристики нагрузок на частоте р, в районе резонансной частоты (/;, /,+| « 0,95-г-1,05); р,+і — частота, определяемая из условия Т (рі) — Т (р<+і).
С помощью (9) — (11) можно определить частотные характеристики нагрузок (перегрузок, изгибающих и крутящих моментов, перерезывающих сил) в различных сечениях консоли.
Амплитудно-частотную характеристику перегрузки в сечении д консоли в зависимости от частоты р действующего единичного синусоидального порыва можно определить по формуле
1 Д п
,(й| - ИИ=-£ (
ч*= 1 /— І I/
(12)
где g — ускорение свободного падения; о — расстояние по перпендикуляру от оси жесткости до заданной точки в сечении е консоли (знак «+» по направлению к задней кромке); п, т — число учитываемых тонов упругих колебаний (здесь п, т ^ 3 согласно (9) — (11)).
Амплитудно-частотную характеристику изгибающего момента в сечении д консоли можно определить так:
.<8>|-ИМЕМ-
(13)
(=1
Амплитудно-частотную характеристику перерезывающей силы в сечении д консоли можно найти с помощью следующей зависимости:
|7д<г(б)|‘
(14)
Используя статистические данные о законе изменения жесткости изгиба крыла вдоль его размаха £/(£) « £/(£= 0) (1 — е)2,5 [ 2], зависимость (14)
можно записать в виде
(т'доСе)! = £/(з£|/;ч|
5£/(е = 0) 21
1
(і5)
І=1
Амплитудно-частотную характеристику крутящего момента в сечении е консоли можно определить следующим образом:
г
|Т„.(!)| = 0/,(«) £ |ф>,| ■ (16)
У-1
Расчетные значения динамических приращений нагрузок в различных сечениях конструкции консоли находим согласно выражениям (1), (2) при принятых по [4] значениях интенсивности турбулентности Т и коэффициента безопасности / для рассматриваемого режима полета, используя результаты расчета |7'н(йо)| в соответствии с (12) — (16).
4. Алгоритм расчета добавочной массы конструкции крыла. Потребную массу крыла самолета тд, обеспечивающую прочность конструкции при действии динамических нагрузок от воздействия непрерывной случайной атмосферной турбулентности, определим как сумму
тА = "Ч.п + Ат,
где тг п — масса крыла, определенная с учетом требований статической прочности для исходного согласно [4] режима установившегося горизонтального полета; Ат— добавочная масса крыла, обеспечивающая прочность конструкции от действия динамических нагрузок при полете в неспокойном воздухе.
На начальном этапе проектирования добавочную массу крыла можно определить по следующей формуле:
Ат
где а—коэффициент, учитывающий кривизну профиля крыла (коэффициент использования высоты кессона) [2]; у — плотность материала; h (g) — толщина профиля консоли в сечении g; одоп, тДШ1 — допускаемые нормальные и касательные напряжения соответственно.
С целью дальнейшего уточнения потребной величины Ат можно воспользоваться зависимостью, полученной на основе гипотезы Мора:
О 4 '
, , °доп (l+c(s))AAJJ,r(s)+0,8<xA(6)6AQp -
где x=l-|--:ri-----------------------———-------------; c(g)— относительная толщина
доп 1,6оДЛ1£рут(£)
крыла в сечении g; 6 — коэффициент, учитывающий отношение толщину обшивки к средней толщине продольного набора крыла [2] .
Таким образом найденная для рассмотренного согласно [4] расчетного случая потребная масса тА должна быть сопоставлена с потребной массой тм для квазистатического нагружения при нормированной перегрузке и при весовом анализе использована большая из них, т. е. т — max (Атл Атм).
5. Пример расчета. В качестве объекта исследования взята консоль крыла большого удлинения. Рассмотрен режим: Н — О, V = 12,5 м/с. Исходные данные, необходимые для проведения численных расчетов, использовались и для расчетов по предлагаемому приближенному алгоритму, но уже в ограниченном объеме. На рис. 2 и 3 приведены графики спектральных плотностей и среднеквадратичных вертикальных ускорений в сечениях консоли g = 0,5 и 0,75. Сопоставление с результатами численного расчета показало, что предлагаемый алгоритм позволяет надежно оценить значения среднеквадратичных перегрузок (погрешность не превышала 15%).
Графики спектральных плотностей и среднеквадратичных изгибающих моментов для сечений £ = 0; 0,5; и 0,75 приведены на рис. 4—6. Полученное отличие аДЛ| , вычисленных по приближенному алгоритму и по результатам численных расчетов, не превышало 15%.
Предлагаемый алгоритм разработан для учета влияния упругости конструкции, динамического нагружения на весовые .характеристики консоли. Поэтому при расчете передаточных функций и среднеквадратичных нагрузок (изгибающих и крутящих моментов, перерезывающих сил) следует ограничиваться рассмотрением диапазона частот упругих колебаний.
Результаты расчетных исследований показывают, что предложенный алгоритм может быть использован на начальных этапах проектирования для определения добавочной массы продольного набора крыла, необходимой для обеспечения прочности от действия динамических нагрузок в условиях воздействия непрерывной атмосферной турбулентности. Погрешность определения указанной добавочной массы консоли крыла большого удлинения находится в пределах ± 5%.
ЛИТЕРАТУРА
1.Кузнецов О. А., П о л и т о в В. В. Определение нагрузок на упругий самолет с системой автоматического управления при полете в неспокойном воздухе // Труды ЦАГИ. — 1982. Вып. 2135.
2. Литвинов В. М. Распределение массы материала в крыле большого удлинения с учетом требований к его прочности и жесткости//Труды ЦАГИ. — 1984. Вып. 2234.
3. Литвинов В. М. Методика расчета массы воздушного винта самолета //Труды ЦАГИ.—1984. Вып. 2234.
4. Нормы летной годности гражданских самолетов СССР / Межведомственная комиссия по нормам летной годности гражданских самолетов и вертолетов СССР. — Изд. 3-е. 1984.
5. Тейлор Дж. Нагрузки, действующие на самолет. — М.: Машиностроение, 1971.
Рукопись поступила 25/VII 1989 г.