CC BY
2Электронный журнал «Труды МАИ». Выпуск № 75
www.mai.ru/science/trudy/
УДК 539.3
Верификация численной модели взаимодействия прямоугольной
пластины с поверхностью воды
Крупенин А.М., Мартиросов М.И.
Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, A-80, ГСП-3, 125993, Россия
e-mail: zeus-russ@yandex.ru
Аннотация
В статье изучаются вопросы численного моделирования эксперимента по сбросу плоских тел в жидкость. Моделируется вертикальный сброс на воду прямоугольной однослойной пластиныиз однородного изотропного материала. Вода моделируется идеальной сжимаемой жидкостью. Учитывается наличие воздушной прослойки между пластиной и водой. Изучается начальный этап взаимодействия, когда гидродинамические силы и давления достигают максимальных значений. Учитывается влияние на динамику пластины гравитационных сил. Задача решается в связной плоскосимметричной постановке.
Ключевые слова: удар, прямоугольная пластина,взаимодействие с жидкостью, численное моделирование.
Введение
В современныхавиационных, ракетных и космических системах применяются транспортные средства и аппараты, вступающие в ударное взаимодействие с жидкостью в процессе эксплуатации или на аварийных режимах работы (экранопланы, спускаемые капсулы и платформы с грузами, гидросамолеты).
Для эффективного проектирования подобных конструкций необходимо учитывать различные эффекты взаимодействия их с жидкостью и применять современные решения, способствующие сохранению прочности и надежности при уменьшении массы(слоистые конструкции, композиционные материалы). Сложный характер движения жидкости в сочетании с нелинейным поведением материалов конструкции исключает возможность аналитического решения поставленной проблемы.
В общем случае задачи взаимодействия конструкций с жидкостью крайне сложны. В прошлом исследователи пользовались упрощенными математическими моделями. Главное упрощение состояло в разделении задачи на две самостоятельных (несвязная постановка): взаимодействие абсолютно твердого тела с жидкостью и исследование напряженно-деформированного состояния конструкции под действием гидродинамических нагрузок. Решения, полученные подобным образом, не могут описать всех нюансов процесса ударного взаимодействия, но позволяют получить важные для экспериментального исследования и практического применения результаты.
Применительно к рассматриваемому вопросу можно привести следующие, полученные ранее, результаты. В работе [1] М. В. Келдыш исследовал задачу об ударе жесткой пластины шириной 2ао несжимаемую жидкость конечной глубины И. Им было показано, что при к > 5авлияние дна уже незначительно.
При ударе затупленных тел о поверхность жидкости граница контакта будет расширяться со сверхзвуковой скоростью. В таких условиях необходимо учитывать сжимаемость жидкости [2].
Сегодня, благодаря бурному развитию вычислительной техники,возможно рассматривать сложные задачи в полной постановке и численно исследовать их, учитывая множество различных факторов. Однако, в связи с наличием разнообразных моделей описания поведения материалов и сложностью отклика материалов на импульсное воздействие, необходимо проверять достоверность получаемых численно результатов. Для этого проводят верификацию численных моделей, сравнивая полученные численно результаты с экспериментальными данными. Вопросу верификации численного эксперимента и посвящена данная работа.
В статье, в прикладном программном пакете ANSYSAUTODYN, численно моделируется эксперимент из работы [3]. Полученные результаты обрабатывались в программном комплексе Mathcad 14.В численном эксперименте рассматривается два случая: с учетом воздушной прослойки между падающей однослойной упругой пластиной и поверхностью воды и без него. В обоих случаях учитывается влияние на падающую в воду пластину гравитационных сил. Задача решается в плоскосимметричной связной постановке. Полученные результаты сравниваются с приведенными в [3]. Проводится анализ максимальных гидродинамических давлений в зависимости от скорости соударения.
Описание эксперимента
В статье [3] приводятся результаты и описание экспериментального исследования плоского удара прямоугольных пластин о воду. На поверхности пластин установлены электромагниты, общая мощность которых достаточна для подъема пластины. С помощью крана и электромагнитов пластина поднималась на заданную высоту, после чего электромагниты выключались, и она падала на поверхность воды. В средней части пластины на ее наружной поверхности вмонтированы мембранные датчики давлений. В процессе испытаний менялась высота падения пластины. Пластины изготовлены из пенопласта ПС-1 с кажущейся плотностью 200 кг/м (другие характеристики материала в статье не указаны).
В данной статье рассмотрим моделирование падения только одного типа пластин размером 1.65м х0.55 м (высота пластин в [3] не указана). Вес пластины40 кг. Недостающий геометрический параметр пластины (высота) легко найти, зная массу, плотность, длину и ширину.
где т - масса пластины, р - плотность материала, 5 - площадь пластины, И - высота пластины.
Недостающие для моделирования характеристики материала найти сложнее. В работе [4] приведены модуль упругости и модуль сдвига для пенопласта ПС-1 с кажущейся плотностью 100 кг/м . В работе [5] — модуль упругости для интересующего нас пенопласта. Чтобы найти модуль сдвига, сделаем допущение, что оба пенопласта обладают одинаковым коэффициентом Пуассона.
Е100 = 5.88 10 Па - модуль упругости на сжатие для пенопласта с кажущейся
-5
плотностью в 100 кг/м .
7
Е200 = 9.8 10 Па - модуль упругости на сжатие для пенопласта с кажущейся
"5
плотностью в 200 кг/м3.
С100 = 2.45 10 Па - модуль сдвига пенопласта с кажущейся плотностью в 100
кг/м3.
Тогда коэффициент Пуассонау равен:
Модуль сдвига для пенопласта с кажущейся плотностью в 200 кг/м3 равен:
Для моделирования эксперимента в ANSYSAUTODYNнеобходим мгновенный модуль объемной деформации К.
В результате испытаний в [3] полученымаксимальные гидродинамические
л
давленияв центре пластины от высоты падения в кг/см .Переведенные в Паскали, они сведены в таблицу 1.
Таблица 1
Номер эксперимента Высота падения, м Ударные давления, Па
1 0.5 1.382-107
2 0.9 1.509-107
3 1.0 1.813107
4 1.2 2.234 107
5 1.3 1.813107
6 2.3 2.46107
В [3] делается вывод о том, что зависимость между ударными давлениями и скоростью удара близка к линейной.
Однако, в статье [3] получено по одному результату для каждого случая, поэтому ни о какой статистической выборке не может идти и речи. Полученные результаты представляют исключительно научную ценность, как позволяющие ориентировочно оценить уровень нагрузок, действующих на пластины со стороны воды.
На рисунке 1 изображены результаты из таблицы 1 (за исключение результата за номером 5).
Рис.1 Максимальные давления в центре пластины от высоты падения
Описание моделирования
Численное моделирование рассматриваемой задачи проводилось в программном продукте АКБУЗАиТОВУК, результаты обрабатывались в пакете МаШеаё 14.
В программном комплексе АКБУБАиТОВУК используется явный метод интегрирования уравнений. Разрешающие соотношения в лагранжевыхдекартовых координатах представляют собой:
Уравнение сохранения массы:
где введены следующие обозначения:
..:■ - 'г- - начальная плотность и начальный объем,
т, V - текущая масса и текущий объем.
Уравнения сохранения количества движения:
где сг,; - тензор напряжений,
- компоненты объемных сил,
л "■" I - компоненты перемещений в соответствующих направлениях.
Уравнение сохранения энергии:
где е - удельная энергия,
- тензор деформаций.
Точкой традиционно обозначается частная производная по времени. В общем случае материалы имеют сложный отклик на динамическую нагрузку. В ANSYSAUTODYN модель материала разбивается на три части:
уравнение состояния, которое выражает изменение объема и описывается шаровым
тензором, модель прочности, которая выражает изменение формы и описывается девиатором, и модель разрушения. Модель материала:
где р - гидростатическое давление,
- главные напряжения,
. - девиатор тензора напряжений.
В качестве уравнения состояния в работе принято соотношение:
где р - гидростатическое давление, А" - мгновенный модуль объемной деформации,
- сжатие,
,.:■ - плотность материала,
..:■ - - начальная плотность.
В качестве модели прочности принято соотношение:
где С - модуль сдвига,
- объем.
Вместе они являются эквивалентом закона Гука.
Для моделирования воды и воздуха использовались стандартные модели, заложенные в ANSYSAUTODYN.
Полное моделирование падение пластин ресурсоемко, поэтому моделировалось падение пластин с высоты в 0.1 м с заданной начальной скоростью, которая вычислялась в зависимости от высоты падения. В дальнейшем пластина под действием гравитационного ускорения продолжает падение до встречи с жидкостью. Так же моделировался случай, когда пластина лежит на поверхности воды и в начальный момент времени приобретает скорость.
На рисунке 2 представлен фрагмент конечно элементной модели эксперимента с учетом воздушной прослойки. Бассейн и воздух моделировались областью шириной 4 м и высотой 7 м, заполнялись элементами Эйлера 0.01 м*0.01 м (400 элементов по ширине и 700 элементов по высоте). Воздух занимает 1 мверхней части моделируемой области (в численном эксперименте с учетом воздушной прослойки воздух занимает 1.1 м). Количество элементов Эйлера, моделирующих жидкость и воздух, равно280000. Пластина моделировалась элементами Лагранжа размерами 0.028 м*0.028 м (8 элементов по высоте и 20 элементов по ширине). Количество элементов, моделирующих пластину, равно 160.
Цифрой 1 на рисунке 2отмечена контрольная точка (геометрический центр нижней поверхности пластины), в которой снимались параметры гидродинамических давлений.
/
/
Рис. 2 Фрагмент конечно элементной модели с учетом воздушной прослойки На рисунке 3 представлен фрагмент конечно элементной модели без учета воздушной прослойки. На рисунке 4 представлен вид бассейна (вся область моделирования).
Рис. 3 Фрагмент конечно элементной модели без учета воздушной прослойки
Рис. 4Вид бассейна (вся область моделирования)
Вычисление скорости
В [3] автором высказано предположение, что для определения скорости соударения пластины о жидкость по начальной высоте падения, формулой, полученной для тела свободно падающего под действием гравитационных сил, пользоваться нельзя(в силу наличия сопротивления воздуха). Для определения скорости падения пластины рассмотрим уравнение свободного падения тела с учетом аэродинамического сопротивления воздуха(точкой традиционно обозначается производная по времени):
где m - масса, падающей пластины, g -ускорение свободного падения, Сх - коэффициент сопротивления, рв - плотность воздуха, V - скорость, падающей пластины, 5 - площадь пластины в плане.
Второе слагаемое в правой части уравнения есть сила аэродинамического сопротивления воздуха.
Разделим уравнение на т и введем обозначение:
2 9 а = —
Тогда уравнение примет вид:
где t - время.
Разделяя переменные, получим:
Интегрируя, получим:
где С\- постоянная интегрирования.
где у0- начальная скорость.
V
Выражениедля V будет иметь вид: (О = ■ - ЫУ),
Интегрируяпо времени, получим:
где С2- постоянная интегрирования, - перемещение.
где х0- начальное перемещение.
Найдем скорость в зависимости от пройденного перемещения. Для этого выразим время из уравнения для перемещения и подставим его в уравнение для скорости. Получаем следующие выражения:
V
(У) = л/—а2 • 1д(_—агссоБ[ е дг-)).
Напомним, что скорость в зависимости от перемещения в случае свободного падения тела без учета аэродинамического сопротивления воздуха, имеет вид:
На рисунке 5 приведен график зависимости скорости соударения от начальной высоты падения для узкой пластины без груза.
"5 Л
Расчетные характеристики: ш = 40 кг, Сх = 1, рв = 1.225 кг/м , 5 = 0.907 м .
£
6
О
й 4
сх ^
и о
2 О
0 0.5 1 1.5 2
Начальная высота падения, м
-С учетом аэродинамического сопротивления
• • • • Без учета аэродинамического сопротивления
Рис. 5Скорость соударения в зависимости от начальной высоты падения Из рисунка видно, что на рассматриваемых высотах отклонения незначительные. А при сбросе грузов с парашютом скорость падения вообще постоянная. Однако, в приведенных выше формулах не учитывается деформирование водной поверхности воздушной подушкой.
По полученным выше формулам рассчитаны начальные характеристики для численного моделирования падающей пластины с высоты 0.1 м. В численном
эксперименте значения начальных скоростей соударения рассмотрено для 18 случаев. Вычисленные характеристики занесены в таблицу 2.За начальную скорость в численном эксперименте без учета воздушной прослойки принималась скорость соударения.
Таблица 2
№ Высота падения, м Начальная скорость, м/с Скорость соударения, м/с Время полного падения, с Время падения с 0.1 м, с
1 0.5 2.792 3.12 0.32 0.034
2 0.9 3.938 4.174 0.429 0.025
3 1.0 4.174 4.397 0.453 0.023
4 1.2 4.608 4.81 0.496 0.021
5 1.3 4.81 5.003 0.517 0.02
6 1.4 5.003 5.188 0.536 0.02
7 2.3 6.467 6.608 0.689 0.015
8 2.7 7.012 7.14 0.747 0.014
9 3.0 7.39 7.511 0.788 0.013
10 3.5 7.974 8.085 0.852 0.012
11 4.0 8.511 8.614 0.912 0.012
12 4.5 9.01 9.105 0.968 0.011
13 5.0 9.476 9.565 1.022 0.011
14 - 20 21.773 - 0.005
15 - 30 29.202 - 0.0033
16 - 40 37.187 - 0.0025
17 - 50 45.437 - 0.002
18 - 100 88.004 - 0.001
Результаты
Результаты моделирования сведены в таблицу 3. Из таблицы видно, что давления, полученные в численном эксперименте с учетом воздушной прослойки, близки с давлениями полученными в [3]. Давления же без учета воздушной прослойки значительно выше.
На рисунке 6 изображены результаты из таблицы 3. Из численного эксперимента без учета воздушной прослойки видно, что гидродинамические давления зависят от скорости начального взаимодействия линейно. Наличие воздушной прослойки приводит к колебанию пластины и поверхности воды до удара, т.е. к их неравномерному взаимодействию (номера 3,4,5,6,7 в таблице 3), и тем больше, чем меньше скорость, падающей пластины. Также из рисунка 6 видно, что при увеличении скорости взаимодействия влияние воздушной прослойки уменьшается.
Скорость удара, определяемая по полученным выше формулам, для случая с учетом воздушной прослойки несколько завышена. Таким образом, если скорость удара принять меньшей, то и график зависимость гидродинамических давлений от скорости, будет стремиться к линейному, что совпадает с результатами в [3].
О 20 40 60 80
-с учетом воздушной прослойки
• • • • без учетаа воздушной прослойки
Рис. 6 Зависимость гидродинамических давлений от скорости соударения
Таблица 3
№ Скорость соударения, м/с Давления в [3], Па Расчетные давления с учетом воздушной прослойки, Па Разница расчетных давлений с учетом воздушной прослойки и экспериментальным, % Расчетные давления без учета воздушной прослойки, Па
1 3.12 1.382 105 1.534 105 9.9 6.33105
2 4.174 1.509 105 1.801105 16.2 7.588 105
3 4.397 1.813105 2.152105 15.75 7.704 105
4 4.81 2.234 105 2.125105 -5.13 8.584 105
5 5.003 1.813105 2.279 105 20.5 8.67105
6 5.188 - 2.127105 - 8.75105
7 6.608 2.46105 3.089 105 20.4 10.877 105
8 7.14 - 2.711105 - 11.736 105
9 7.511 - 3.2105 - 12.344 105
10 8.085 - 3.756 105 - 13.29105
11 8.614 - 4.12105 - 14.161105
12 9.105 - 4.466 105 - 14.958 105
13 9.565 - 4.604 105 - 15.717105
14 21.773 - 18.173105 - 36.208 105
15 29.202 - 29.819105 - 45.3105
16 37.187 - 44.778 105 - 56.523 105
17 45.437 - 58.525 105 - 69.234 105
18 88.004 - 136105 - 140.28 105
Заключение
В работе проведено численное моделирование эксперимента по вертикальному сбросу однослойных упругих пластин прямоугольной формы на воду. Рассмотрено влияние на гидродинамические давленияскорости начального взаимодействия пластины и жидкости, а также влияние воздушной прослойки на величину максимальных гидродинамических давлений.
В результате получено, что в программном комплексе ANSYSAUTODYN возможно адекватно моделировать сбросы грузов на воду. Численные значения гидродинамических давлений близки к экспериментальным. Гидродинамические давления зависят от скорости взаимодействия по линейному закону, что тоже хорошо согласуется с экспериментом [3].
Влияние воздушной прослойки существенно сказывается на величине максимальных давлений только на малых скоростях взаимодействия. При увеличении скоростей взаимодействия влияние воздушной прослойки уменьшается.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 11-01-00540 а, 12-01-00566_а).
Библиографический список
1. Келдыш М. В. Удар пластины о воду, имеющую конечную глубину// Труды ЦАГИ, 1935. Вып. 152. С. 13-20.
2. Горшков А. Г., Григолюк Э. И. Взаимодействие упругих конструкций с жидкость (удар и погружение). Л.:Судостроение. 1976. 200с.
3. Бельгова М. А. Экспериментальное исследование плоского удара// Труды ЦНИИМФ, 1969. Вып. 117. С. 15-28.
4. Пособие по физико-механических характеристикам пенопластов и сотопластов. М.: Стройиздат. 1977. 79с.
5. Справочник по пластическим массам// под редакцией Гарбара М. И., М.: Химия. 1969. 520с.
