CC BY
43
1546
Механика деформируемого твердого тела Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (4), с.1546-1548
УДК 539.3
МОДЕЛИРОВАНИЕ ВЗАИМОДЕЙСТВИЯ УДАРНЫХ ВОЛН С ПОРИСТЫМИ ГАЗОПРОНИЦАЕМЫМИ ПАКЕТАМИ МЕТАЛЛИЧЕСКИХ СЕТОК
© 2011 г. А.В. Кочетков1, В.И. Романов2, И.А. Ходыкина1
1НИИ механики Нижегородского госуниверситета им. Н.И. Лобачевского 2Российский федеральный ядерный центр - Всероссийский НИИ экспериментальной физики,
г. Саров
kochetkov@mech.unn.ru
Поступила в редакцию 15.06.2011
Формулируются начально-краевые задачи, описывающие взаимодействие ударных волн с газопроницаемыми пакетами металлических плетеных сеток. Рассматриваются среды из двух взаимопроникающих кон -тинуумов - деформируемый пакет сеток, как скелет двухфазной среды, и поровый газ, движущийся внутри скелета. Численное решение задачи проводится на основе схемы Годунова в эйлерово-лагранжевых подвижных разно стных сетках. Приводятся результаты экспериментальных исследований влияния внутренних плоских и цилиндрических пакетов плетеных сеток на деформации силовых оболочек камер при внутреннем взрывном нагружении. Численные результаты соответствуют экспериментальным данным.
Ключевые слова: взрыв, плетеные сетки, двухфазная среда, деформация, моделирование, эксперимент.
Одним из эффективных способов ослабления газодинамических нагрузок от ударных волн является использование пакетов мелкоячеистых проволочных тканевых сеток, размещаемых перед преградой. В ряде исследований даются оценки эффективности подобных многослойных экранов с точки зрения снижения газодинамических параметров проходящих ударных волн. За исключением [1], пакет сетки предполагался недефор-мируемым и неподвижным в пространстве.
Динамическое поведение пакета сетки с содержащимся поровым газом описывается на основе уравнений динамики двух взаимопроникающих континуумов, каждый из которых имеет свои скорости, напряжения (давления) и температуры. При формулировке уравнений использованы обычные для подобных смесей предположения [2]. Наряду с этими допущениями полагаем, что проволочки металлической сетки ориентированы по двум ортогональным направлениям в плоскости слоев и эти направления являются осями ортотропии пакета сеток, как скелета анизотропной сплошной среды. Вдоль этих направлений происходит упругопластическое растяжение скелета среды. Третью ось ортотропии образует направление, нормальное к плоскостям сеток, вдоль которого происходит упругопластичес-кое сжатие скелета.
В декартовой (V = 0) или цилиндрической системе координат (V = 1, 02 — ось вращения) урав-
нения динамики порового газа в форме законов сохранения массы, импульса и энергии имеют вид:
д (P1V1) + д о
= —V
dt dryr^r P°ßrM1 д^У^) + д
+ ^ (P1U ßr«i) (P?ß z"1) =
+ ^ (ßrP +ßr P0<) +
дг
r дt
+ f (ßz P0«1V1) = —V-^ — nfr
дг r
+f- (ßrP0V1)+ + <ß +ß zp0v?) ^vM^ — nfz,
(1)
d_ дt
д а
+ д- ßr
дг
a1p1
e1 +
r
2 2 V\ Wf + vf
w
P°W1
+
yj
22 u2 + v,2 e1 + —-L
+ P1u1
+
+
дг
= —V
P0v1
e +
u{ + Vf 2
+P1V1
ßrV1
r — nq
P1
e1 +
u{ + v{ 2
+ P1
Р1 = (к — 1)р°е1, Т1 = е1/с(7]). Здесь и далее нижний индекс 1 относится к газу, 2 — к твердой
2
с)
ß
z
Моделирование взаимодействия ударных волн с пористыми газопроницаемыми пакетами
1547
компоненте; ^ — время; г, г — координаты;р — давление; р0 — истинная плотность газа; и, V — компоненты скорости; е — внутренняя энергия; Т — температура; cv — удельная теплоемкость; к — показатель адиабаты; /г, д — межфазные силы и тепловой поток, действующие на «элементарный узел» (рис. 1, штриховая линия) твердой компоненты со стороны порового газа, п — количество таких узлов в единице объема смеси. Из системы (1) при а1 = вг = вг = 1 и п = 0 следуют уравнения газовой динамики для однородной газовой среды. Тогда уравнения динамики скелета среды примут вид:
dp2 , d
+ -
d
p2U2
■ "(P2u2) + 3- (P2U2> = -V , dt dr dz r
d(p2u2) d . 2 . д . 0 .
- + (P2U22 -Crr ) + T- (P2U2V2 - Srz) = dr dz
dt
V 2
= - (2Srr + Szz - P2U2) + nfr r
d(P2V2) , d
'2V2 dt
2
+ 3- (P2U2V2 - Szr ) + dr
V
+ "dz(P2V2 -azz) = T(Srz -p2u2v2) + nfz
de'
dt
■ + Ur,
du2
dz
de'zz
— e
dt ^dv
- + u2
2
r de'rr dr
dv2Л
dr
dgZz dr du2
+v
de'
dz
1
---e v
+ v2
du2
"dr 3 de'
dz
v
dr dz
dv2
1
---ev
dz 3
/ _ , / / ч
e ее =-(err +ezzЛ
+v
de rz
2 "dT
du dz
— >
de rz dt du
- + u2
+
2
2 + dv2 dr
V
dz
dev dt
de rz dr dv2
dr
/
3p2 dt
Здесь 51.., £'/ — девиаторные компоненты тензо-
У ^
ров напряжений и деформаций; полные деформации рассчитываются обычным образом:
еу = 4 +6. ^ (''. = г, г, 9), 5. — символ Кронекера.
Законы связи между компонентами напряжений и деформаций строятся в направлениях осей ортотропии.
По нормали к слоям сетки считаем напряжения зависящими от плотности о пп = /п (р 2, р 2),
по оси ортотропии вдоль проволочек сетки на плоскости rz отт = fT(eTT, e*). При V = 1 по окружному направлению Оее = fT(еее, еее). Сдвиговые компоненты тензора напряжений в системе координат, связанной с этими осями, полагаем равными нулю. Параметры, отмеченные верхним индексом *, представляют собой максимальные значения, достигнутые частицей при нагр ужении в соответствующем направлении и необходимы для описания разгрузки частиц среды, испытывающей необратимые деформации. При взаимодействии твердой и газовой компонентов учитываем в качестве межфазных сил: силы сопротивления частичек твердой фазы при их обтекании поро-вым газом; силы Стокса вязкого трения. Также учитывается конвективный теплообмен через межфазную поверхность.
Численное решение уравнений производится по схеме С. К. Годунова [3]. С целью ее адаптации к уравнениям динамики деформирующихся взаимопроникающих континуумов разработаны алгоритмы расчета задачи распада произвольного разрыва для высокопористой упругопластичес-кой среды, а также алгоритмы решения уравнений динамики порового газа в эйлерово-лагран-жевых сетках с учетом изменяющихся коэффициентов проходных сечений и объемной концентрации пор. Для расчета параметров на контактных границах «двухфазная среда - газ» используется алгоритм, реализующий последовательно две задачи о распаде произвольного разрыва - «скелет - газ», «поровый газ - газ». Потоки импульса и энергии от обеих фаз усредняются с учетом значений коэффициентов проходных сечений и используются как граничные потоки для интегрирования в области однородного газа. На рис. 1 показана схема экспериментов со взрывной камерой (1 - базальтопластиковый слой, 2 - сталь, 3 -пакет сеток), а на рис. 2 - сравнение расчетных и экспериментальных зависимостей окружных деформаций базальтопластикового слоя от времени на внешнем силовом корпусе при z = 0.
Наблюдается хорошее соответствие численных и экспериментальных результатов в течение первого периода нестационарных колебаний.
Таким образом, предложенная математическая и численная модель деформирования газопроницаемых пакетов плетеных сеток может успешно использоваться для оценки их защитных свойств при взрывных воздействиях.
Авторы выражают благодарность М. А. Сыру-нину за предоставленные экспериментальные материалы и С.В. Крылову за помощь в расчетах.
e
2
1548
А.В. Кочетков, В.И. Романов, И.А. Ходыкина
Работа выполнена при частичном финансировании РФФИ (грант Ж»09-08-00711а).
Список литературы
1. Абакумов А.И. и др. Численная модель деформирования противоосколочной сетки при взрывном
нагружении // Тр. ВНИИЭФ. Математическое моделирование физических процессов. 2006. №10. С. 16-30.
2. Нигматулин Р.И. Динамика многофазный сред. Ч.1. М.: Наука, 1987.
3. Численное решение многомерных задач газовой динамики / Под ред. С.К. Годунова. М.: Наука, 1976.
MODELLING THE INTERACTION OF SHOCK WAVES WITH POROUS GAS-PERMEABLE METAL GRID PACKAGES
A. V. Kochetkov, V.I. Romanov, I.A.Hodykina
Initial boundary-value problems describing the interaction of shock waves with gas-permeable packages of woven metal grids are formulated. Media consisting of two mutually permeating continua- a deformable grid package as the skeleton of the two-phase medium and a pore gas flowing inside the skeleton - are considered. The problem is numerically analyzed using Godunov's scheme with Euler-Lagrange movable difference meshes. The results of experimentally studying the effect of internal plane and cylindrical packages of woven grids on the deformation of the supporting shells of the chambers under internal explosive loading are presented. The numerical results agree with the experimental data.
Keywords: explosion, woven grids, two-phase medium, deformation, modeling, experiment.
