Моделирование аэрогидроупругости несущих поверхностей в авиатике мягких летательных аппаратов Текст научной статьи по специальности «Механика и машиностроение»

УДК 519.6, 519.95

Ю.М. Давыдов Институт машиностроительной механики, г. Москва

МОДЕЛИРОВАНИЕ АЭРОГИДРОУПРУГОСТИ НЕСУЩИХ ПОВЕРХНОСТЕЙ В АВИАТИКЕ МЯГКИХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ

Abstract

The aviatics of soft flight apparatus is quite science-capacious field. Calculations by using Large-Particles Method permit to reduce the weight and price ofparachutes. All types of modern and perspective parachutes (cosmic atmospheric re-entry, hypersonic, space shuttle drogue, round, conical, gliding-type, ribbon, with constructive holes etc.) are investigated by this method. The algorithm for solving of problem on stress-strain state in canopy is proposed. Canopy may be considered as pliable shell or measure. All practically important problems of aeroelasticity of parachute system (canopy opening; parachute braking with form changing; interaction of shoch waves and parachute, etc.) are investigated by this algorithm of Large-Particles Method..

В настоящее время активно развивается перспективная область машиноведения - авиатика мягких летательных аппаратов. К аппаратам данного класса относится огромный класс летающих устройств: парашюты, парапланы, баллюты, дельтапланы, аэростаты и др. Автором исследованы все эти летающие устройства, в том числе рассмотрены все типы парашютов: планирующие (однооболочковые и двухоболочковые), ленточные, тормозные, спасательные и др. Исследованы многокупольные системы. Особая сложность исследования парашютных систем заключается в том, что форма купола (и строп) не является неизменной, как в случае самолёта, дирижабля или спускаемого космического аппарата, а определяется в результате аэроупругого взаимодействия парашюта с окружающей средой.

Автором исследованы практически все задачи авиатики мягких летательных аппаратов: раскрытие купола парашютных систем, торможение парашюта с

изменением его формы, взаимодействие ударных волн с парашютом и др. В аэродинамических расчётах им учтены упругость купола и строп (в том числе рассмотрены соответствующие разрушению предельные прочностные характеристики ткани купола и лент строп), вес парашюта и строп, реальный вес спасаемого объекта, внутреннее трение материала строп и купола.

С начала XX века центр проведения всех основных работ по теории и практике конструирования парашютов находится в нашей стране. Г.Е. Котельников в 1911 году изобрёл ранцевый парашют. В течение многих десятилетий основополагающими книгами по парашютостроению являлись монографии Н.А. Лобанова [1,2 и др.]. Значительный вклад в элементарную теорию парашютов внесли работы Х.А. Рахматулина [3,4 и др.]. Важные разделы, посвященные парашютостроению, содержатся в фундаментальных монографиях академика А.М. Матвеенко [5,6 и др.] и Генерального конструктора, академика Г.И. Северина [6 и др.].

Натурные и аэрофизические эксперименты дают небольшое количество информации, к тому же они экономически весьма дороги. Поэтому основным способом

исследования парашютных систем в настоящее время является численное моделирование [7-18 и др.]. Наиболее мощным методом численного эксперимента, с помощью которого решены всевозможные задачи парашютостроения на всех режимах как дозвуковых, так и сверхзвуковых, является метод крупных частиц [19 и др.]. Единообразие вычислений этим методом позволяет впервые корректно поставить и решить многие комплексные практически важные задачи и достичь высокой точности получаемых результатов [20].

Опишем алгоритм метода крупных частиц для решения задачи о напряженно-деформированном состоянии купола парашюта. Купол рассматривается как мягкая оболочка (мембрана). В данной работе мы будем использовать традиционный подход теории упругости.

Для того чтобы показать, каким образом получаются соотношения метода крупных частиц применительно к расчету напряженно-деформированного состояния мягких оболочек, рассмотрим произвольную мягкую оболочку (рис. 1,а), с которой свяжем сопутствующую криволинейную систему координат а1, а2. Оболочка движется относительно внешней декартовой системы отсчета х1, х2, х3.

Дифференциальное уравнение движения мягкой оболочки в векторной форме запишем следующим образом:

Р^ = ^^1 +Г12е2 А ] +

Эt Эа1

+ ^ [(Т12 е1 +Т22е 2 )^1 ] + Ар£ , (1)

Эа 2

где р - масса единицы площади оболочки; верхний индекс "о" обозначает начальное (раскройное, недеформированное) состояние; g = д/§п§22—§п , ^ = д/^§17,

^2 g22 - соответственно дискриминант метрического тензора поверхности и

коэффициенты Ламе; glm - компоненты метрического тензора (1,т=1,2); и - вектор скорости оболочки; t - время; Т11, Т12, Т22 - компоненты тензора мембранных усилий

оболочки; е1 , е2 - единичные векторы локального базиса сопутствующей системы координат; Ар - вектор поверхностной нагрузки. При анализе деформированного состояния мягких оболочек воспользуемся степенями удлинения

1П = жп/ж°п = к!к , п=1,2, или величинами относительных удлинений

еп - 1

и углом сдвига

0 = Х° - С = Х° - агссо8( gl2І^11 g22 ).

Связь между мембранными усилиями и относительными удлинениями, и углом сдвига для изотропной оболочки можно записать в следующем виде:

Е Е

Т11 =----2 (е1 + ^2 ) , Т22 =-2 (£2 + У£1), Т12=С0,

1-V 1-V

где Е - жёсткость материала оболочки при растяжении; G = Е / [2(1+у)] - жёсткость материала оболочки при сдвиге; V - коэффициент Пуассона.

Вычисленные по относительным удлинениям мембраны усилия должны удовлетворять условиям существования мягкой оболочки [20]:

Т11 + Т12 + 2Т12 §'12> 0;

Тп Т12 - 2Т* > 0.

Предположим, что нашей задачей является определение напряжённо-деформированного состояния мягкой оболочки при существенном изменении во времени её формы, т.е. координат хх=хх(аьа2,0, Х=1, 2, 3, точек оболочки М(аьа2), под действием заданных нагрузок с учётом связей, наложенных на оболочку. Знание координат хx точек оболочки позволяет вычислить, как показано ниже, интересующие нас параметры - относительные удлинения и напряжения. Отметим, что задача расчёта мягкой оболочки при больших перемещениях является нелинейной, но применение при численных расчетах малых шагов по времени и пространству позволяет с приемлемой точностью линеаризовать исходную задачу. Рассмотрим способ получения разрешающих уравнений метода крупных частиц применительно к расчёту напряжённо-деформированного состояния мягкой каркасированной оболочки при больших перемещениях. Особенно просто получить соотношения метода крупных частиц для исследования динамики мягких оболочек, изготовленных из ткани. Координатные линии а1, а2 в этом случае целесообразно располагать вдоль основы и утка ткани, из которой изготовлена оболочка. Специфика работы таких оболочек проявляется в отсутствии сопротивления изгибу, сжатию и сдвигу. Исходное дифференциальное уравнение движения мягкой оболочки (1) без учета сдвиговых усилий в проекциях на оси координат хx запишем в следующем виде:

Эи х Э(Т1) х Э(Т2) х

Р0 ^7 = ^ +-Э^ +АРх. (2)

Эt Эа1 Эа 2

где Т1, Т2 - натяжения, действующие по направлениям координатных линий а1, а2; индекс X = 1, 2, 3 обозначает проекцию соответствующей величины на ось хх.

На оболочку наложим расчётную сетку (рис. 1,б). Узловым точкам присвоим номера г,]. Вспомогательной расчётной сеткой, сдвинутой относительно основной на полшага, разобьём оболочку на «крупные частицы». Сосредоточим массу каждой крупной частицы в соответствующей узловой точек г, ]. Величину массы будем вычислять по формуле

М] =+Р °] (Ь°+1/2,] + Ку2,])(Ь°]+1/2 + К]-1/2) / ^

где рг] - масса единицы площади крупной частицы; йу±1/2 - расстояние между смежными узловыми точками.

Расчёт движения крупных частиц ведётся шагами по времени, каждый из которых содержит три этапа. На первом (эйлеровом) этапе определяют скорости движения крупных частиц. Разностное уравнение, аппроксимирующее уравнение (2),

!+1 к а

разрешённое относительно скорости иг] в момент времени t =t +Аt, запишем следующим образом:

(и‘+1 )х ={(и‘ )5 + [(Т+1/2] -]2 -]2 ){ +

+ (Ар! )х (И+/2] + *М,'2] )(Л,']+1,2 + ]2) / 4](Аt/M,^ ))з, (3)

где Т±, ]] 2 - погонные усилия в оболочке, АР] - усредненный перепад давления, действующий на крупную частицу и в общем случае являющийся функцией координат, времени, скорости движения оболочки и т.д., Ь - коэффициент демпфирования, гасящий нефизичные осцилляции (в случае их возникновения в процессе решения), связанные со способом аппроксимации континуальной задачи дискретным аналогом. Коэффициент Ь подбирается опытным путем и зависит от условий конкретной задачи. С целью устранения возможной немонотонности решения можно также пользоваться корректировкой скорости [21].

Рис. 1. Мягкая оболочка (а) и расчетная сетка, наложенная на оболочку (б)

б

На втором (лагранжевом) этапе при использовании метода крупных частиц обычно вычисляются потоки массы через границы ячеек на временном слое /к+1. В нашем случае принимаем, что обмен массами между крупными частицами не происходит (АМіі = 0), а крупные частицы перемещаются в новое положение,

характеризуемое координатами х^+1. Используя вычисленные на первом этапе

скорости, находим эти координаты по формуле

XГ =Х)Г + (4+1)х А/.

(4)

На заключительном этапе по полученным координатам вычисляем новые расстояния между крупными частицами (узловыми точками):

к

к+1

Ч±1/2,} '■

1

3 2

І [(X 5 )& - (X х )Г‘ ]

х=1

2

Пн

£+1

У

Зная новые расстояния между узловыми точками ячеек, определяем степени удлинения, возникающие в оболочке:

где Еі±\щ, Еіл±т - приведенные жёсткости оболочки, зависящие от жёсткости материала оболочки и жёсткости каркаса, если последний проходит через данную крупную частицу. Приведенные жёсткости оболочки вычисляем по формуле:

где Е1, Е2 - погонные жёсткости оболочки вдоль координатных линий а1, а2, Ес -продольная жёсткость каркаса. Жёсткости оболочки и каркаса могут быть функциями степеней удлинения, скорости деформирования и зависят от вида напряжённо-деформированного состояния оболочки (одноосное, двухосное). В соответствии с особенностями работы мягких оболочек на вычисляемые натяжения накладываем условия:

Погонные условия в материале оболочки и натяжения в каркасе подсчитываем по формулам:

При этом вычисляемые натяжения подчиняются условиям типа (5). Так как уравнения (3) - (4) реализуют явную конечно-разностную схему, то для устойчивости счёта интервал времени выбирается из условия устойчивости Куранта:

Решение уравнений (3) - (4) осуществляется с использованием соответствующих начальных и граничных условий, выполняемых в узловых точках и зависящих от конкретной задачи. Стационарное решение получается в результате установления статического состояния.

По степеням удлинения вычисляем натяжения в оболочке:

(5)

&< ^шшл/Ршт/Е

тах •

В качестве тестового примера с помощью предложенного алгоритма метода крупных частиц была решена задача о раздувании эластичной прямоугольной мембраны, аналитическое решение которой для установившегося состояния известно. Рассматривалась первоначально плоская длинная прямоугольная мембрана,

закреплённая вдоль длинных сторон, нагруженная равномерным перепадом давления. Результаты расчёта в безразмерном виде и сопоставление их с аналитическим решением даны на рис. 2 - рис. 4.

Приняты обозначения: х/а - безразмерная координата, Я/а - радиус поперечного сечения деформированной мембраны; ж = Ара/Е - параметр, связывающий перепад давления Ар, ширину первоначально плоской мембраны а и жёсткость материала мембраны Е; //а - относительная стрелка прогиба мембраны; Т/Е - безразмерное натяжение мембраны. Численные значения исходных данных, принимавшиеся при расчёте мембраны, были следующими:

Рис. 2. Результаты решения задачи о раздувании эластичной прямоугольной мембраны

Рис. 3. Вычисленная величина относительной стрелки прогиба мембраны

Т/Е

Рис. 4. Вычисленная величина безразмерного натяжения мембраны

ширина мембраны а = 1,2 м, длина мембраны Ь = 6 м, перепад давления, действующий на мембрану, Ар = 3,93 кПа, варьируемая жёсткость материала мембраны Е =2,35 ... 78,5кН/м. В силу симметрии рассматривалась четвёртая часть мембраны, которая разбивалась на 36 крупных частиц. Связь между натяжениями и относительными удлинениями принималась подчиняющейся закону Гука. На рис. 2 представлены формы поперечного сечения мембраны в зависимости от параметра ж. Сплошной линией показаны результаты аналитического решения, треугольниками обозначены данные численного расчета.

Как видно из рис. 2, совпадение форм поперечного сечения мембраны, полученных аналитически и численно, наблюдается во всем диапазоне изменения параметра 0,06 < ж < 2. При этом относительные удлинения мембраны изменяются от 5 до 106%. На рис. 3 ив табл. 1 дано сопоставление безразмерной стрелки прогиба мембраны в зависимости от параметра ж при аналитическом и численном решениях. На рис. 4 и в табл. 2 приведено безразмерное натяжение мембраны в зависимости от того же параметра. В обоих случаях наблюдается вполне удовлетворительное совпадение аналитического и численного решений, что позволяет сделать вывод о хорошей сходимости численного расчета к точному решению данной нелинейной задачи.

Таблица 1.

Безразмерная стрелка прогиба мембраны при различных значениях ж

ж Численный расчёт Точное решение Погрешность, %

0,06 0,1467 0,1458 0,6173

0,96 0,4583 0,4544 0,8492

2,0 0,7250 0,7253 -0,0414

Таблица 2.

Безразмерное натяжение мембраны при различных значениях ж

ж Численный расчёт Точное решение Погрешность, %

0,06 0,05563 0,05581 0,3225

0,96 0,4860 0,4825 0,7254

2,0 1,058 1,070 1,075

Заметим, что реальная ткань является неизотропной оболочкой. Она неодинаково растягивается (и садится) в разных направлениях. Это связано с технологией

производства тканых материалов. В процессе тканья нити утка переплетаются с нитями основы. Нити утка проходят через ткань с меньшей волнистостью, чем нити основы. Нити основы получают большую волнистость, и поэтому по основе ткань тянется больше, чем по утку. В табл. 3 приведены полученные экспериментально показатели сопротивления ткани растяжению.

Таблица 3.

Показатели сопротивления ткани растяжению

Продолжительность Расстояние Расстояние Расстояние

испытаний (час) по основе (м) по утку (м) по диагонали (м)

0,5 9 7,5 24

6 9 7,5 24

24 10,5 7,85 24,25

12009 12 9,0 25

При кратковременных полётах парашютов указанные неравномерности

растяжения не имеют практического значения. В случае же долговременно

работающих тормозных устройств (типа плавучих якорей и др.) может существенно

исказиться форма первоначального раскроя.

Библиографический список

1. Н.А. Лобанов, Ф.С. Чуриков. Основы теории и расчёта парашютов. - М.-Л.: ОНТИ НКПТ СССР, Главная редакция авиационной литературы, 1936. - 160 с.

2. Н.А. Лобанов. Основы расчёта и конструирования парашюта. - М.:

Машиностроение, 1965. - 363 с.

3. Х.А. Рахматулин. Теория раскрытия парашюта // Техника воздушного флота. -1940. - № 8. - С. 79-89.

4. Х.А. Рахматулин. Теория осесимметричного парашюта. Часть 1 //Обтекание парашютов и проницаемых тел потоком газа. Институт механики МГУ //Научные труды (№ 35). - М.: Изд-во МГУ, 1975. - С. 3-55.

5. М.Г. Акопов, А.М. Матвеенко, Г.И. Северин и др. Системы оборудования летательных аппаратов. - М.: Машиностроение, 1995. - 496 с.

6. С.М. Егер, А.М. Матвеенко, И. А. Шаталов. Основы авиационной техники. - М.: Машиностроение, 2003. - 720 с.

7. Ю.М. Давыдов. Решение задач парашютостроения методом крупных частиц. - В кн.: Сборник докладов научно-технической конференции. Вып. 1(11). - М.: НИИ парашютостроения, 1993. - С. 29-40.

8. Ю.М. Давыдов, В.А. Мозговой. Эффект аномального аэродинамического нагрева при спуске парашюта по траектории // Доклады академии наук. -1993. - Т. 330. -№ 1. - С. 48-51.

9. Ю.М. Давыдов, В.А. Тутурин. Взаимодействие ударной волны с аэроупругим парашютом // Доклады академии наук. -1996. - Т. 351. - № 3. - С. 329-331.

10. Ю.М. Давыдов. Устойчивость полёта парашютов // Доклады академии наук. - 1998.

- Т. 365. - № 5. - С. 626-631.

11. Yu.M. Davydov. Numerical Modelling of Aerohydrodynamic Problems of Parachute Systems. - Two-Days Workshop on Parachute Technologies and Space Recovery Systems. - Spain, Madrid, 1992. - Р. 24.

12. Yu.M. Davydov. Parachute and Other Wind Device Engineering. - East European Conference on Wind Engineering. - Poland, Warshaw, 1994. - Р. 31.

13. Yu.M. Davydov. Space Vehicle Re-Entry, Landing and Safety Parachute Systems. - 47th International Astronautical Congress. Report № IAF-96, - Vol.4.04. - China, Beijin, 1996. - Р. 9.

14. Yu.M. Davydov. Re-Entry and Landing Parachute Systems for Space Vehicles. - 48th International Astronautical Congress. Report № IAF-97. - Vol.4.09. - Italy, Turin, 1997.

- Р. 10.

15. Ю.М. Давыдов, О.В. Рысев. Гидродинамика парашютных систем. - М.: НИИ парашютостроения, 1991. - 176 с.

16. Ю.М. Давыдов. Моделирование нелинейных задач механики парашютов //Международный конгресс «Нелинейный анализ и его приложения». - М.: АНН, 1998. - С. 166.

17. Ю.М. Давыдов. Исследование аэродинамики, аэроупругости, нагрева, прочности и устойчивости полёта парашютных систем //Современные проблемы аэрогидродинамики. Т. 2. - М.: Изд - во ин-та прикладной механики РАН, 1999. -С. 35-54.

18. Ю.М. Давыдов. Влияние длины строп на аэродинамику парашюта //В кн.: III Международный конгресс «Актуальные проблемы механики сплошных и сыпучих сред», 2-я Международная конференция по парашютостроению и ветровой инженерии. - М.: Изд-во «Нефть и газ», 2000.

19. Ю.М. Давыдов. Крупных частиц метод //В кн.: Математическая энциклопедия. Т. 3.

- М.: Советская энциклопедия, 1982. - С. 125-189.

20. Ю.М. Давыдов. Аэродинамика, гидроупругость и устойчивость полета парашютных систем. - М.: НАПН РФ, 2000. - 256 с.; 2-е изд-е, дополненное, 2001.

- 306 с.

21. Б.В. Гулин, Р.И. Давыдов, В.В. Ридель. Численное исследование динамики мягкой оболочки в одноосном состоянии //Нелинейные проблемы аэроупругости: Труды семинара. Вып. 11. - Казань: Изд - во физ.-тех. ин-та Казанского филиала АН СССР, 1979. - С. 43-58.

Получено 12.08.2004.