Научная статья на тему 'МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПНЕВМАТИЧЕСКОГО АМОРТИЗАТОРА С ВОЗДУХОПРОНИЦАЕМОЙ ТКАНЕВОЙ ОБОЛОЧКОЙ'

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПНЕВМАТИЧЕСКОГО АМОРТИЗАТОРА С ВОЗДУХОПРОНИЦАЕМОЙ ТКАНЕВОЙ ОБОЛОЧКОЙ Текст научной статьи по специальности «Механика и машиностроение»

CC BY
41
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Труды МАИ
ВАК
Ключевые слова
ПНЕВМАТИЧЕСКИЙ АМОРТИЗАТОР / ПНЕВМОАМОРТИЗАТОР / ТКАНЕВАЯ ОБОЛОЧКА / СИСТЕМА МЯГКОЙ ПОСАДКИ / ВОЗДУХОПРОНИЦАЕМОСТЬ / СТАТИСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ПРИЗЕМЛЕНИЯ

Аннотация научной статьи по механике и машиностроению, автор научной работы — Аверьянов Игорь Олегович, Васильченко Андрей Геннадьевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Аверьянов Игорь Олегович, Васильченко Андрей Геннадьевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MATHEMATICAL MODEL OF AIR DAMPER WITH PERMEABLE FABRIC SHELL

One of the phase of landing process for cargo dropped system with air dampers (System) is when the cargo moves on the damper during its compression. It may lead to cargo bounce and sometimes to its overturn. The System trajectory on this phase depends on experiment conditions: landing velocity, mass that comes to the damper, environment conditions and regime of gas expelling from air damper camera that is controlled with valves. Generally, materials of dampers shells are not hermetic - air permeability may have an influence on System landing dynamics. In this article we describe the mathematical model of the air damper with air permeability fabric shell for its further using in the tasks of System landing parameters prediction with variable experiment conditions. Considering dynamic model consists of cargo (mass) and cylindrical air damper with a valve. The law of the mass moving describes by one-dimensional equilibrium equation, where reaction from the damper is defined with equation of state for ideal gas. Forth integration method is used to solutions on each step. The model of gas expelling through the shell has been developed. Analysis of existing permeability theories shows that Rysev’s law is the the most appropriate for the considering task: this theory is based on Rahmatulin’s non-linear law of pressure gradient to gas velocity, but transformed with similarity theory to nondimensional parameters, related only to fabric characteristics - it makes possible not to consider experiment conditions if we knew them. Based on assumption that gas conditions inside the camera keep constant during the integration step of dynamic mass-to-damper model, we find the value of flow rate through the fabric shell on this step (with formulas of Rysev’s law) and add it to the value of valve flow rate. There are two experiment-calculated tasks are presented to evaluate the influence of permeability on System dynamics: 1 - we put the mass on air damper and find fabric characteristics; 2 - we make a drop test with the same mass and damper. Analysis shows well correlation of experimental results and calculated results received from the developed model. Comparison of the solutions with and without permeable shells of the air damper shows the essential difference on the phase of mass bounce. Thus, we can conclude that the influence of gas permeability through the damper shell has to be taking into account to define the dynamics of mass-to-air damper system, especially in the tasks of bounce and overturn prediction.

Текст научной работы на тему «МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПНЕВМАТИЧЕСКОГО АМОРТИЗАТОРА С ВОЗДУХОПРОНИЦАЕМОЙ ТКАНЕВОЙ ОБОЛОЧКОЙ»

Труды МАИ. 2022. № 125 Trudy MAI, 2022, no. 125

Научная статья УДК 539.39

DOI: 10.34759/^-2022-125-11

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПНЕВМАТИЧЕСКОГО АМОРТИЗАТОРА С ВОЗДУХОПРОНИЦАЕМОЙ ТКАНЕВОЙ

ОБОЛОЧКОЙ

Игорь Олегович Аверьянов1Н, Андрей Геннадьевич Васильченко2

1АО «МКПК «Универсал» им. А.И. Привалова», Москва, Россия 2 АО «НИИ Парашютостроения», Москва, Россия

1 i. averyanov@mail.

[email protected]

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

Ключевые слова: пневматический амортизатор, пневмоамортизатор, тканевая оболочка, система мягкой посадки, воздухопроницаемость, статистическое моделирование процесса приземления.

Для цитирования: Аверьянов И.О., Васильченко А.Г. Математическая модель пневматического амортизатора с воздухопроницаемой тканевой оболочкой // Труды МАИ. 2022. № 125. БОТ: 10.34759/Ы-2022-125-11

Original article

MATHEMATICAL MODEL OF AIR DAMPER WITH PERMEABLE

FABRIC SHELL

Igor O. Averyanov1H, Andrey G. Wasilchenko2

1MKPK «Universal» by A.I. Privalov, Moscow, Russia 2JSC "NII Parachutostroyeniya", Moscow, Russia 1 i. averyanov@mail. ruB [email protected]

Abstract. One of the phase of landing process for cargo dropped system with air dampers (System) is when the cargo moves on the damper during its compression. It may lead to cargo bounce and sometimes to its overturn. The System trajectory on this phase depends on experiment conditions: landing velocity, mass that comes to the damper, environment conditions and regime of gas expelling from air damper camera that is controlled with valves. Generally, materials of dampers shells are not hermetic - air permeability may have an influence on System landing dynamics. In this article we describe the mathematical model of the air damper with air permeability fabric shell for its further using in the tasks of System landing parameters prediction with variable experiment conditions.

Considering dynamic model consists of cargo (mass) and cylindrical air damper with a valve. The law of the mass moving describes by one-dimensional equilibrium equation, where reaction from the damper is defined with equation of state for ideal gas. Forth integration method is used to solutions on each step.

The model of gas expelling through the shell has been developed. Analysis of existing permeability theories shows that Rysev's law is the the most appropriate for the considering task: this theory is based on Rahmatulin's non-linear law of pressure gradient to gas velocity, but transformed with similarity theory to nondimensional parameters, related only to fabric characteristics - it makes possible not to consider experiment conditions if we knew them. Based on assumption that gas conditions inside the camera keep constant during the integration step of dynamic mass-to-damper model, we find the value of flow rate through the fabric shell on this step (with formulas of Rysev's law) and add it to the value of valve flow rate.

There are two experiment-calculated tasks are presented to evaluate the influence of permeability on System dynamics: 1 - we put the mass on air damper and find fabric characteristics; 2 - we make a drop test with the same mass and damper. Analysis shows well correlation of experimental results and calculated results received from the developed model. Comparison of the solutions with and without permeable shells of the air damper shows the essential difference on the phase of mass bounce.

Thus, we can conclude that the influence of gas permeability through the damper shell has to be taking into account to define the dynamics of mass-to-air damper system, especially

in the tasks of bounce and overturn prediction.

Keywords: soft landing system, air damper, dropped cargo landing, fabric permeability, Rysev's model

For citation: Averyanov I.O., Wasilchenko A.V. Mathematical model of air damper with permeable fabric shell. Trudy MAI, 2022, no. 125. DOI: 10.34759/trd-2022-125-11

Введение

Одним из этапов процесса посадки груза с закрепленным на нём пневматическим амортизационным устройством (далее всё вместе - объект) [1-2] является движение груза на амортизационном устройстве (АУ, амортизатор), сопровождающееся возможным отскоком объекта, в некоторых случаях приводящим к его опрокидыванию [3]. Траектория движения приземляющегося объекта на данном этапе во многом зависит от условий эксперимента: скорости посадки, приходящей на амортизатор массы груза, параметров среды, а также режима стравливания газа из камеры амортизатора, регулирование которого производится с помощью клапанов. Материалы оболочек камер амортизаторов в общем случае не являются герметичными - их воздухопроницаемость может оказывать влияние на динамику приземления объекта. В данной работе описывается разработанная математическая модель пневматического амортизатора с воздухопроницаемой тканевой оболочкой для использования в задачах прогнозирования параметров приземления объекта, характеризующихся вариативностью условий проводимых экспериментов [4].

Обзор теорий воздухопроницаемости

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

V = к ■Ар (1)

Работы многих исследователей показали, что при возрастании скорости протекания среды зависимость (1) перестаёт быть линейной [5]. Рахматуллиным, занимавшимся исследованиями воздухопроницаемости тканевых материалов, была предложена следующая форма квадратичной зависимости перепада давлений от скорости протекания среды, в которой коэффициенты а и Ь подбирались под условия проводимых экспериментов

Ар = а ■ V + Ь ■ V2 (2)

Таким образом, связь перепада давления и скорости проницания становилась всё более сложной; уравнения (1) и (2) - это последовательные разложения зависимости Ар (V.) в ряд всё с большим количеством членов.

Установленные линейная (1) и нелинейная (2) зависимости в том или ином виде используются в различных отраслях народного хозяйства, в том числе в текстильной промышленности при определении воздухопроницаемости

технических тканей. Так современные методы моделирования процессов взаимодействия воздухонаполняемых оболочек с твёрдыми телами (например процесс посадки объекта с системой воздушной амортизации [6-12], процесс раскрытия автомобильных подушек безопасности [13]) основаны на использовании численных методов, в которых проницаемость тканевых материалов задаётся линейным законом: в частности, согласно [14] условие проницаемости записывается как

В уравнении (3): CDARCY - коэффициент Дарси; CDF - дополнительный коэффициент, применяемый для учёта вязкости воздушного потока на больших скоростях протекания; PORE - коэффициент пористости материала; V - скорость протекания среды в направлении i; PERMi - коэффициент проницания; дРа/дх -градиент давления по направлению i; f - некоторая функция накладываемая на градиент давления.

Анализируя формулу (3), можно заключить следующее: формула представляет собой модифицированную запись зависимости (1) при коэффициенте CDF = 0 и использовании выражения для градиента давления без наложения какой-либо функции f; в общем виде формула (3) выражает нелинейную зависимость перепада давления от скорости протекания среды и при задании соответствующих параметров может являться аналогом зависимости (2).

дР

+ CDF ■ |V I) ■ PORE ■ V = PERMi • f ^, i = 1,2,3 1 Л} 1 1 дхг

(3)

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

Герасимато, взяв за основу зависимость (2) и используя методы теории размерности и подобия [15], предложил обратиться к физической сути явления проницания и ввести безразмерные переменные, характеризующие закон протекания среды через ткань. С целью получения безразмерных критериев самой ткани и придания им физического смысла Рысевым разработана видоизменённая запись структуры закона протекания среды через ткань - введены безразмерный коэффициент сопротивления ткани с^ и эффективный диаметр нити ёэф. Получив

своё развитие в работах [16-18], данный подход к оценке воздухопроницаемости в настоящее время используется при проектирования парашютных систем [18].

Согласно [16-18], закон протекания газообразной среды через ткань запишется следующим образом

с 1

°¥ = 1 + (4)

С¥0

Ар

С¥

р- V 2/2'

КеГ = V-р-Лэф/Л

В формулах (4): сг - безразмерный коэффициент сопротивления ткани в конкретных условиях эксперимента; Яег - число Рейнольдса из теории подобия, характеризующее процесс протекания среды через ткань; С - предельный безразмерный коэффициент сопротивления ткани при Яег ^да, в дальнейшем названный «сопротивляемостью»; V - скорость протекания воздуха через ткань; Ар - избыточное давления; р - плотность среды; / - коэффициент кинематической вязкости. Таким образом, коэффициенты а и Ь из (2) (в которых смешаны характеристики ткани и среды) выражены через два параметра с^ и Лэф,

являющихся характеристиками конкретной ткани.

В случае, когда условия эксперимента известны, оболочки воздушных

амортизаторов изготовлены из обладающих малой пористостью тканей и их

нагруженность избыточным внутренним давлением невысока для оценки влияния

воздухопроницаемости оболочек амортизаторов на динамику приземления объекта

возможно использование линейной зависимости в виде (1). Однако в общем случае,

при решении задач динамики работы тканевых пневматических амортизаторов,

следует использовать нелинейную теорию протекания среды через ткань -

зависимости (2)-(3) и (4). При этом, использование закона протекания, выраженного

зависимостями (4), позволяет по известным характеристикам ткани проводить

расчётную оценку параметров воздухопроницаемости тканевой оболочки

амортизатора при различных условиях эксперимента [4]. Следует также отметить,

8

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

Таким образом, задачей данного исследования является разработка математической модели (ММ) пневматического амортизатора с воздухопроницаемой тканевой оболочкой с использованием нелинейной теории протекания среды через ткань (4) с целью уточнения динамики приземления объекта на этапе обжатия его амортизационной системы.

ММ пневматического амортизатора с газонепроницаемой оболочкой

Предварительно описывается ММ тканевого пневматического амортизатора (амортизатора, АУ) с газонепроницаемой оболочкой, аналогичная [19-21]: стравливание газа из амортизатора осуществляется только с помощью регулируемых клапанов; закон изменения площади стравливающих клапанов известен и является функцией избыточного давления ^ = Г (Ар); заполняющий

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

Рисунок 1. Схема работы амортизатора

Тогда, с учётом принятых допущений, сопротивление амортизатора обжатию является функцией скорости его обжатия Р = f(V), определяемой скоростью вертикального поступательного движения груза массы т. Уравнение движения находящегося в поле сил тяжести g груза для случая вертикальной посадки запишется следующем образом:

Р + mg = та (5)

Интегрируя уравнения движения груза как твёрдого тела по времени с заданным шагом Лt с использованием метода конечных разностей [22], процесс обжатия амортизатора сводится к рассмотрению его г квазистационарных

состояний, в каждом из которых оболочка амортизатора находится в равновесном состоянии.

Режим истечения газа из АУ со стравливающим клапаном и газонепроницаемой стенкой на I - м шаге интегрирования по времени определяется таким параметром, как действительный мгновенный секундный расход

\

(6)

я-т.

Здесь: / - коэффициент расхода воздуха; ^ - коэффициент, зависящий от природы газа; ^ - площадь стравливающего клапана; р - давление в АУ; Я - газовая постоянная воздуха (287,05 Дж/кг-град); Т - абсолютная температура газа.

Указанные в правой части формулы (6) переменные, характеризующие состояние газа в камере амортизатора ^, р и Т, а также передаваемое от амортизатора на груз усилие ^, определяются по известным методикам, изложенным, например в [19-21]. В то же время понятно, что в [19-21] описан предельный случай воздухопроницаемой ткани для пневматического амортизатора, которой в природе нет - любая ткань имеет хоть малую, но проницаемость. Учтём реальные характеристики проницаемости ткани по зависимости (4).

1

ММ пневматического амортизатора с воздухопроницаемой оболочкой

Принимаются следующие допущения:

1) оболочка амортизатора, представляющая собой тканевый пакет (внутренняя и внешняя оболочки с армирующими лентами), может быть представлена по аналогии с [16-18] в качестве эквивалентной однородной ткани с характеристиками и с^ ;

2) характеристики воздухопроницаемости рассматриваемой эквивалентной ткани не зависят от её деформации;

3) динамический процесс истечения воздуха через оболочку может быть рассмотрен как сумма квазистационарных состояний на каждом шаге интегрирования по времени Лt процесса обжатия амортизатора со стравливающим клапаном.

Рассматриваемый процесс протекания газообразной среды через ткань представляет собой процесс преодоления этой среды некоего барьера в виде тканевой стенки и пристеночных - областей [18]. Схематичный график перепада давлений на стенке оболочки представлен на рисунке 2 - а). На нём тканевая стенка, разделяющая пространство на 2 области, с левой стороны обдувается воздушным потоком, имеющим скорость V. Создаваемая сила, воздействующая на стенку, вызвана перепадом давлений Лр = рх - р2, где р1 и р2 давления до ткани и за ней. За пределом 8 -области справа скорость воздушного потока по-прежнему равна V. При этом сама ткань характеризуется такими параметрами как эффективный диаметр нити йэф и сопротивляемость с^ .

Слой

Рисунок 2. Схема протекания воздуха через тканевую оболочку: а) модель воздухопроницаемости ткани [18]; б) предложенная модель воздухопроницаемости

Предлагаемая схема протекания воздуха через оболочку представлена на рисунке 2 - б). На ней тканевая оболочка представляется в виде границы, разделяющей замкнутый объём и атмосферу. Воздух внутри оболочки (слева от границы) представляется в виде прилегающих друг к другу слоёв. Равномерность расположения слоёв воздуха обуславливается сделанным ранее допущением о том, что внутри АУ распределение давления и температуры равномерное. Перепад давлений Ар = р1 - ра осуществляется скачком на границе самой оболочки.

За некоторое время Аt под воздействием перепада давлений Ар слой воздуха

проникает через оболочку и оказывается вне объёма АУ. При этом этот слой

13

воздуха в процессе проникновения через оболочку приобретает некоторую скорость, среднее значение которой равно VВП. Положив, что шаг интегрирования по времени процесса обжатия АУ Лt и время проникновения через границу одного слоя Лt связаны соотношением Лt = п • Лt, а также учитывая допущение № 3 о квазистационарности процесса на интервале времени Лt (в смысле постоянства значений перепада давления Лр, плотности и температуры среды внутри АУ на рассматриваемом шаге интегрирования по времени уравнений, описывающих обжатие АУ), можно заключить, что

1 п

Vвп =1 ^ (7)

п м

Здесь: Viввп - средняя скорость потока на -м шаге интегрирования процесса обжатия АУ; VJВП - средняя скорость слоя газа, проходящего через границу на у -м шаге интегрирования процесса протекания.

Таким образом, объём воздуха, вышедшего на г -м шаге интегрирования из АУ через оболочку площадью ^ за время Лt

Щвп = X - Л • 3) = Vвп •Л/ • 3 (8)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

у=1

Пусть известны характеристики ткани йэф и с^ . Тогда для г -го шага

интегрирования по времени, подставив второе выражение формулы (4) в первое и решив полученное квадратное уравнение относительно Яе, получим

Ке, = 1

1 2

л

8 - Ар ■ р ■ Л \

1 + ^ Р\ эф -1

С¥0 - Л

Выражая скорость из третьего выражения формулы (4), можно определить среднюю скорость потока воздуха через ткань

у ВП= Кег - л

Р,'Лэф

(10)

Определив объём вышедшего через оболочку воздуха ¥ВП по формуле (8) и

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

тВП = ШгВП-р, тВП

АОВП = 1 Аt

(11)

Действительный мгновенный секундный расход для АУ со стравливающим клапаном и воздухопроницаемой оболочкой получается совмещением формулы (6) и второго выражения формулы (11)

/гРГ

1 Л 1 Г 1

Я-Т

+ АО,

ВП

(12)

Таким образом, в ММ учтены реальные характеристики

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

1

Экспериментальное исследование влияния воздухопроницаемости тканевой оболочки амортизатора на динамику приземления системы груз-амортизатор

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

г = 1,00 с г = 2,00 с г = 5,00 с

Рисунок 3. Кинограмма обжатия амортизатора

На рисунке 4 сплошной тонкой линией представлены экспериментальные

данные по изменению перемещения верхнего основания амортизатора во времени,

полученные обработкой видеозаписи процесса обжатия (тонкими вертикальными пунктирными линиями обозначены моменты времени, соответствующие кадрам на рисунке 3), тонкой горизонтальной пунктирной линией показана остаточная высота амортизатора, сплошной жирной - результат расчёта с заложенными характеристиками тканевой оболочкой dэф = 0,0001 м, = 106.

Для демонстрации устойчивости алгоритма разработанной математической модели процесса динамического обжатия амортизатора на рисунке 4 штрих-пунктирной линией показан график расчётного перемещения верхнего основания, полученный для газонепроницаемой (герметичной) оболочки при А Ош = 0 в формуле (12).

§ ос

0.7

0.6

0.5

0.4

Ж 0.3 ю

ш =

| 0.2

о.

и Ш

0.1

\ / \./Л 7\ \ / V л. 'V \ л /■ \./л <\ У V V л Л У V

\;й/ Л 4 ^ Т расч.

\ 1" г

п \/ Т эксп, 4-->

1 _1_

И остаточн.

; ; ;

0.0 0.5 1.0 1.5 2.0 2.5

4,5 5.0 5,5

3.0 3.5 4.0 Время, с

■ Расчёт с учётом воздухопроницаемости оболочки АУ -Экспериментальные данные

Расчёт без учёта воздухопроницаемости оболочки АУ

6.0

6.5 7.С

Рисунок 4. Перемещение верхнего основания амортизатора при обжатии клетью, м

Анализ представленных на рисунке 4 графиков расчётных и экспериментальных значении перемещений верхнего основания амортизатора позволяет сделать заключение, о удовлетворительном совпадении траекторий движения при установлении режима обжатия (начиная с t = 1,0 с) и времени полного истечения воздуха из камеры амортизатора (~ 6,0 с). В то же время штрих-пунктирная линия расчётного перемещения верхнего основания, полученная для газонепроницаемой (герметичной) оболочки вообще не совпадает с данными эксперимента. Некоторые отличия (20%) расчётных и экспериментальных значений на первоначальном этапе обжатия АУ с воздухопроницаемой оболочкой объясняются несоответствием форм потери устойчивости оболочки амортизатора: на втором кадре кинограммы процесса обжатия (рисунок 3, t = 0,25 с) видно, что амортизатор приобрёл наклон и имеет место излом оболочки в месте контакта с клетью с правой стороны, в то время как в ММ данные эффекты не учитываются -деформации цилиндрической оболочки амортизатора симметричны относительно его центральной оси. После 2-й с от начала движения отличие расчёта и эксперимента снижается до 4%. При этом, не смотря на различия в амплитудах некоторых колебаний, в целом наблюдается соответствие периодов колебательных

пр°цессов: Трасч. « Тэксп. (см. рисунок 4).

Таким образом, по результатам данной расчётно-экспериментальной работы по плавному обжатию цилиндрического амортизатора продемонстрированы общая работоспособность ММ амортизатора с воздухопроницаемой оболочкой и

соответствие заложенных характеристик ткани ^ = 0,0001 м, сЖп = 106 физическим

свойствам тканевой оболочки рассматриваемого амортизатора.

На рисунке 5 представлена кинограмма копрового сброса того же амортизатора, что был рассмотрен в предыдущем примере, но с открытым клапаном, закон стравливания которого известен. В начальный момент времени избыточное давление внутри камеры отсутствует.

г = 0,59 с г = 0,93 с г = 1,40 с

Рисунок 5. Кинограмма копрового сброса клети с амортизатором

Аналогично предыдущему случаю результат обработки видеозаписи копрового сброса представлен на рисунке 6 в виде графика изменения положения верхнего основания амортизатора (тонкая сплошная линия) во времени. За начальный момент времени принят момент касания нижнем основанием амортизатора пола (рисунок 5, г = 0,00 с).

0.00 ■ • '

0.0 0.1 0.2 0.3 ОД 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1Д 1.5

Время, с

Расчёте учётом воздухопроницаемости оболочки АУ -Экспериментальныеданные

--Расчёт без учёта воздухопроницаемости оболочки АУ

Рисунок 6. Вертикальное перемещения верхнего основания амортизатора при

копровом сбросе, м

Там же приведены расчётные значения соответствующих перемещений, полученные с использованием математической модели амортизатора, истечение газа из которого задавалось формулой (12): сплошная толстая линия соответствует случаю амортизатора с воздухопроницаемой стенкой, толстая пунктирная линия - с герметичной стенкой (в формуле (12) принималось ЫЭБП = 0).

Анализ полученных результатов позволяет сделать вывод об удовлетворительном совпадении результатов эксперимента с результатом расчёта по предлагаемой математической модели, учитывающей воздухопроницаемость тканевой оболочки амортизатора. Наблюдаемое расхождение этих результатов в конечной фазе приземления амортизатора (при t > 0,59 с), объясняется тем обстоятельством, что в эксперименте при складывании тканевая стенка амортизатора заслоняла стравливающий клапан, препятствуя выходу через него воздуха и тормозя дальнейший ход обжатия. Сопоставление результатов расчёта рассматриваемого процесса приземления объекта с учётом и без учёта воздухопроницаемости оболочки амортизатора, при одинаковых условиях эксперимента, показывает существенность влияния учёта эффекта воздухопроницаемости оболочек на траекторию движения приземляющегося объекта, прежде всего при его отскоке - рисунок 6, t = 0,36 с. Величина отскока по ММ герметичного амортизатора на 40% выше величины отскока по ММ воздухопроницаемого амортизатора при идентичности величины отскока по ММ воздухопроницаемого амортизатора данным эксперимента.

Заключение

На основании нелинейной теории протекания среды через ткань, записанной в форме (4), разработана математическая модель пневматического амортизатора со стравливающим клапаном и воздухопроницаемой тканевой оболочкой,

позволяющая по известным характеристикам тканевой оболочки йэф и с^ уточнить

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

Список источников

1. Трямкин А.В., Емельянов Ю.Н. Математическая модель процесса торможения объекта десантирования парашютно-реактивной системой // Труды МАИ. 2000. № 1. URL: http://trudymai.ru/published.php?ID=34731

2. Трямкин А.В., Скиданов С.Н. Исследование процесса наполнения парашютных систем // Труды МАИ. 2001. № 3. URL: https: //trudymai.ru/published.php?ID=34686

3. Аверьянов И.О., Зинин А.В. Математическая модель процесса приземления недеформируемого груза с амортизирующим устройством на жесткую площадку // Труды МАИ. 2022. № 124. URL: http://trudymai.ru/published.php?ID=167067. DOI: 10.34759/trd-2022-124-12

4. Аверьянов И.О. Применение метода Монте-Карло для прогнозирования надёжности процесса приземления систем мягкой посадки с пневмоамортизаторами

// Труды МАИ. 2020. № 115. URL: http://trudymai.ru/published.php?ID=119896. DOI: 10.34759/trd-2020-115-03

5. Куличенко А.В. Разработка моделей и экспериментальных методов изучения воздухопроницаемости текстильных материалов: дисс. ... д.т.н. - М.: МГУДТ, 2005. - 340 с.

6. Аверьянов И.О., Сулейманов Т.С., Тараканов П.В. Разработка обобщённой методики расчета систем мягкой посадки с пневмоамортизаторами // Труды МАИ. 2017. № 92. URL: http://trudymai.ru/published.php?ID=77448

7. Qu Pu, Yang Zhen, Shi Rui. Research on Airbags Landing System for Airborne Vehicle Airdrop // Journal of Information and Computational Science, 2015, vol. 12(5), pp. 2035-2042. DOI: 10.12733/jics20105798

8. Per-Olof Marklund, Larsgunnar Nilsson. Simulation of airbag deployment a coupled fluid-structure approach // 7th International LS-Dyna Users Conference, 2002. URL: https://www.dynalook.com/conferences/international-conf-2002

9. Yves de Lassat de Pressigny, Vincent Lapoujade. Numerical simulation of ground impact after airdrop // 5th European LS-Dyna Users Conference, 2005. URL: https://www.dynalook.com/conferences/european-conf-2005

10. Yves de Lassat de Pressigny, Thierry Baylot. Simulation of the impact on ground of airdrop loads to define a standard worst case test // 6th European LS-Dyna Users Conference, 2007. URL: https://www.dynalook.com/conferences/european-conf-2007

11. Masoud Alizadeh, Ahmad Sedaghat, Ebrahim Kargar. Shape and orifice

optimization of airbag systems for UAV parachute landing // International Journal of

Aeronautical and Space Sciences, 2014, vol. 15(3), pp. 112-121. URL: https://doi.org/10.5139/IJASS.2014.15.3.335

12. Pan Cao, Xuyan Hou, Yongbin Wang, Meng Li, Xiaoshan Rao, Yuetian Shi. Flexible airbag cushioning for Martian landing based on discrete element method // Advances in Space Research, 2019, vol. 63, issue 8, pp. 2566-2583. URL: https://doi.org/10.1016/j.asr.2019.01.007

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

13. Huajian Zhou, Zhihua Zhong, Manjiang Hua. Design and Occupant-Protection Performance Analysis of a New Tubular Driver Airbag // Research Vehicle Engineering, 2018, vol. 4, issue 2, pp. 291-297. URL: https://doi.org/10.1016/j.eng.2018.03.015

14. LS-DYNA KEYWORD USER'S MANUAL. URL: www.dynasupport.com/manuals/

14. Седов Л.И. Методы подобия и размерности в механике. - М.: Наука, 1977. -440 с.

15. Рысев О.В., Вишняк А.А., Чуркин В.М., Юрцев Ю.Н. Динамика связанных тел в задачах движения парашютных систем. - М.: Машиностроение, 1992. - 288 с.

16. Рысев О.В., Пономарев А.Т., Васильев М.И., Вишняк А.А., Днепров И.В., Мосеев Ю.В. Парашютные системы. - М.: Наука.Физматлит, 1996. - 288 с.

17. Лялин В.В., Морозов В.И., Пономарев А.Т. Парашютные системы. Проблемы и методы их решения. - М.: Физматлит, 2009. - 576 с.

18. Емельянов Ю.Н., Павлов А.С., Титов В.А. Проектирование системы мягкой посадки приземляющегося объекта. - М.: Изд-во МАИ, 1988. - 62 с.

19. Пономарев П.А. Исследование и выбор рациональных параметров пневматического амортизатора для посадки дистанционно-пилотируемых летательных аппаратов. Дисс. ... к.т.н., Москва, МАИ, 2000, 145 с.

20. Пономарев П.А., Скиданов С.Н., Тимохин В.А. Расширение диапазона применения пневмоамортизаторов в системах мягкой посадки с использованием разрывных элементов // Труды МАИ. 2000. № 2. URL: http://trudymai.ru/published.php?ID=34708

21. Амосов А.А., Дубинский Ю.А., Копченова Н.В. Вычислительные методы для инженеров. - М.: Высшая школа, 1994. - 544 с.

References

1. Tryamkin A.V., Emel'yanov Yu.N. Trudy MAI, 2000, no. 1. URL: http: //trudymai .ru/eng/publ i shed.php?ID=34731

2. Tryamkin A.V., Skidanov S.N. Trudy MAI, 2001, no. 3. URL: https://trudymai .ru/eng/publ i shed.php?ID=34686

3. Aver'yanov I.O., Zinin A.V. Trudy MAI, 2022, no. 124. URL: http://trudymai.ru/eng/published.php?ID=167067. DOI: 10.34759/trd-2022-124-12

4. Aver'yanov I.O. Trudy MAI, 2020, no. 115. URL: http://trudymai.ru/eng/published.php?ID= 119896. DOI: 10.34759/trd-2020-115-03

5. Kulichenko A.V. Razrabotka modelei i eksperimental'nykh metodov izucheniya vozdukhopronitsaemosti tekstil'nykh materialov (Development of models and experimental

methods for research of fabric materials permeability), Doctor's thesis, Moscow, MSUDT, 2005, 340 p.

6. Aver'yanov I.O., Suleimanov T.S., Tarakanov P.V. Trudy MAI, 2017, no. 92. URL: http: //trudymai .ru/eng/publ i shed.php?ID=77448

7. Qu Pu, Yang Zhen, Shi Rui. Research on Airbags Landing System for Airborne Vehicle Airdrop, Journal of Information and Computational Science, 2015, vol. 12(5), pp. 20352042. DOI: 10.12733/jics20105798

8. Per-Olof Marklund, Larsgunnar Nilsson. Simulation of airbag deployment a coupled fluid-structure approach, 7th International LS-Dyna Users Conference, 2002. URL: https://www.dynalook.com/conferences/international-conf-2002

9. Yves de Lassat de Pressigny, Vincent Lapoujade. Numerical simulation of ground impact after airdrop, 5th European LS-Dyna Users Conference, 2005. URL: https://www.dynalook.com/conferences/european-conf-2005

10. Yves de Lassat de Pressigny, Thierry Baylot. Simulation of the impact on ground of airdrop loads to define a standard worst case test, 6th European LS-Dyna Users Conference, 2007. URL: https://www.dynalook.com/conferences/european-conf-2007

11. Masoud Alizadeh, Ahmad Sedaghat, Ebrahim Kargar. Shape and orifice optimization of airbag systems for UAV parachute landing, International Journal of Aeronautical and Space Sciences, 2014, vol. 15(3), pp. 112-121. URL: https://doi.org/10.5139/IJASS.2014.15.3.335

12. Pan Cao, Xuyan Hou, Yongbin Wang, Meng Li, Xiaoshan Rao, Yuetian Shi. Flexible

airbag cushioning for Martian landing based on discrete element method, Advances in

Space Research, 2019, vol. 63, issue 8, pp. 2566-2583. URL: https://doi.org/10.10167i.asr.2019.01.007

13. Huajian Zhou, Zhihua Zhong, Manjiang Hua. Design and Occupant-Protection Performance Analysis of a New Tubular Driver Airbag, Research Vehicle Engineering, 2018, vol. 4, issue 2, pp. 291-297. URL: https://doi.org/10.10167i.eng.2018.03.015

14. LS-DYNA KEYWORD USER'S MANUAL. URL: www.dynasupport.com/manuals/

15. Sedov L.I. Metody podobiya i razmernosti v mekhanike (Methods of similarity and dimensionality in mechanics), Moscow, Nauka, 1977, 440 p.

16. Rysev O.V., Vishnyak A.A., Churkin V.M., Yurtsev Yu.N. Dinamika svyazannykh tel v zadachakh dvizheniya parashyutnykh system (Dynamics of connected bodies in the tasks of parachute systems movement), Moscow, Mashinostroenie, 1992, 288 p.

17. Rysev O.V., Ponomarev A.T., Vasil'ev M.I., Vishnyak A.A., Dneprov I.V., Moseev Yu.V. Parashyutnye sistemy (Parachute systems), Moscow, Nauka.Fizmatlit, 1996, 288 p.

18. Lyalin V.V., Morozov V.I., Ponomarev A.T. Parashyutnye sistemy. Problemy i metody ikh resheniya (Parachute systems. Problems and methods of their solving), Moscow, Fizmatlit, 2009, 576 p.

19. Emel'yanov Yu.N., Pavlov A.S., Titov V.A. Proektirovanie sistemy myagkoi posadki prizemlyayushchegosya ob"ekta (Developing of soft landing system), Moscow, Izd-vo MAI, 1988, 62 p.

20. Ponomarev P.A. Issledovanie i vybor ratsional'nykh parametrov pnevmaticheskogo amortizatora dlya posadki distantsionno-pilotiruemykh letatel'nykh apparatov (Analysis

and choice of rational parameters of pneumatic shock absorber for the landing aircrafts), Doctor's thesis, Moscow, MAI, 2000, 145 p.

21. Ponomarev P.A., Skidanov S.N., Timokhin V.A. Trudy MAI, 2000, no. 2. URL: http: //trudymai .ru/eng/published.php?ID=34708

22. Amosov A.A., Dubinskii Yu.A., Kopchenova N.V. Vychislitel'nye metody dlya inzhenerov (Numerical methods for engineers), Moscow, Vysshaya shkola, 1994, 544 p.

Статья поступила в редакцию 09.07.2022 Статья после доработки 11.07.2022 Одобрена после рецензирования 18.07.2022 Принята к публикации 25.08.2022

The article was submitted on 09.07.2022; approved after reviewing on 18.07.2022; accepted for publication on 25.08.2022

i Надоели баннеры? Вы всегда можете отключить рекламу.