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

УДК 22.251

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

© В.М. Дубровин, Т. А. Бутина, Н.С. Полякова МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

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

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

Введение. Моделированию волновых процессов взаимодействия ударных волн с конструкциями посвящены многие исследования [1-15], в том числе при воздействии волн, вызванных электромагнитными источниками [11], а также для конструкций из неклассических, гибких материалов [4-6], обладающих высокими демпфирующими свойствами.

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

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

В общем случае при нормальном падении ударной волны на жесткую преграду давление на ее внешней поверхности за счет отражения повышается. Значение этого давления может быть найдено по формуле [1-3]

АР = К АР

отр -"-отр^" ф 5

где Котр = 2 н----; ДРф — избыточное давление во фронте

р 1 + ак/ДРф ф

ударной волны; ак = 7,2 Па.

Последнее соотношение показывает, что для слабых ударных волн типа акустических (ДРф ~ 0) Котр = 2, для сильных ударных

волн (ДРф ^го) Котр ^ 8. График зависимости от избыточного давления во фронте ударной волны представлен на рис. 1.

Процесс взаимодействия ударной волны с цилиндрической оболочкой условно можно разделить на два характерных периода: 1) период погружения — от момента встречи оболочки с ударной волной до момента полного погружения в ударную волну; 2) период обтекания — от момента полного погружения оболочки в ударную волну до конца фазы движения.

Определение нагрузок на цилиндрическую оболочку при ее взаимодействии с ударной волной сводится к определению давления в каждой точке ее поверхности для различных моментов времени. Определение нагрузок на цилиндрическую оболочку в период погружения является задачей в достаточной мере сложной. Ее решение в соответствии с имеющимися теоретическими исследованиями [4-7] может быть проведено следующим образом. Рассматривается волна заданной формы импульса, падающая на оболочку перпендикулярно ее образующей. Схема взаимодействия волны с оболочкой в плоской постановке представлена на рис. 2. Результирующее давление Р

К(

6

4

2'

О

ДРф, кг/см

Рис. 1. Зависимость коэффициента Рис. 2. Схема взаимодействия отражения ударной волны от избы- ударн°й в°лны с °болочк°й:

точного давления во фронте волны 1— поперечное сечение оболочки; 2 — падающая волна; 3 — отраженная волна

в каждой точке цилиндрической поверхности складывается из давления в падающей волне Р1 и давления в отраженной волне Р2. Для определения этих давлений справедливы соотношения

Р = -

АфР Эф

Я Эт

Р2 =-

АфР Эу

г=1

Я Эт

г=1

где Аф — скорость распространения фронта ударной волны; р(т) — плотность среды в ударной волне; Я — радиус оболочки; ф(г, 9, т) , у (г, 9, т) — потенциалы скоростей в падающей и отраженной

волнах; 9 — угловая координата; т = —---безразмерный параметр

Я

времени; г — безразмерный радиус в долях Я.

В общем случае потенциал ф( г, 9, т) есть функция, определяемая видом падающей волны и контуром обтекаемого тела, т. е. является известной функцией. Потенциал у(г, 9, т) должен удовлетворять волновому уравнению

Э 2у 1 Эу 1 Э 2у Э 2у Э2г г Эг г2 Э92 " Эт2 !

граничным условиям

Эф Эу

Эг Эг

= 0,

г=1

условиям излучения на бесконечности и нулевым начальным услови-Эу

ям у =-= 0

Эт

т=0

Для определения потенциала отраженной волны у(г, 9, т) воспользуемся его представлением на многомерной римановой поверхности с углом 9, непрерывно меняющимся в пределах -да <9<да. Выполнив преобразования Фурье по 9 вида

¥ (и) = { / (9)е-2™9а9

и преобразования Лапласа по т вида

да

¥ (8 ) = | / (т)е-а т,

получим уравнения (2) и граничные условия (3) следующего вида:

д у * 1 ду * дт 2

V

т дт дФ * 1 ду:

5 + — 1у* = 0;

дт2

т дт

= 0,

т=1

где 5 = а + /ю — комплексная переменная; V — вещественный параметр; Ф * (г, 0,5), у * (г, 0,5) — изображение по Лапласу потенциалов ф(г, 0, т), у(г, 0, т).

Знак «*» означает, что для этих функций выполнено преобразование Фурье.

Решение уравнения (4) с граничными условиями (5) и условиями изучения на бесконечности можно представить в виде

1 ЭФ *

у* =

sKV дт

т=1

= к J_v(s)- Jv (s) ;

, JV

где KV (s) — функция Макдональда; KV (s) =

2 sin vk

функция Бесселя мнимого аргумента; J-V — функция Бесселя, полученная из JV заменой v на -V.

Соответствующее изображение для давления отражения на цилиндрической поверхности

P* =

рДф Kv (s) дФ*

R KV(s) дт

т=1

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

Kvífl = -1 + I s-1

KV (s) 2

'V 2

Л

-2 I 2 J I -3 s -I V — Is

0

где 0( 4 ) — остаточный член разложения.

В рассматриваемом случае изображение для давления отражения примет вид

РАФ

Pi=-

R

дФ 1 дФ -+-

3 дФ 3 дФ

-+-Т--+

V2 дФ*

дт 2s дт 8s2 дт 8s3 дт 2s2 дт

V дФ

s3 дт

т=1

Обратное преобразование Фурье от первых четырех слагаемых в последнем выражении может быть выполнено непосредственно, а для обращения последних двух членов необходимо иметь выражение для падающей волны. Если рассматривать ударную волну, в которой давление изменяется по закону

Ap(t) = АР

( t Л" 1 - —

V т+У

где т+ — период фазы сжатия в ударной волне; " = 1,АРф — показатель степени бинома, зависящий от избыточного давления (для малых т), потенциал скоростей можно записать в виде

АРфR -а(гcos6-1+т) TJ, a i , ч

Ф = —— e v H (r cos 0-1 + т),

рАфа

где а =- — постоянная, определяющая скорость падения давле-

тэф

ния за фронтом ударной волны; тэф = 2 т+ — эффективное время,

п +1

определяемое из условия равенства импульсов, действительной и приведенной эпюр давления; Н — функция Хевисайда.

В рассматриваемом случае преобразования Лапласа и Фурье дают

Ф, = ^да е„.со>^ = ^(^г

рАф 8 + а-"да рАф а 8 + а

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

Р =

раф

R

до i до з до з до --+----^-+ —^-

дг 2s дг 8s2 дг 8s3 дг

+

г=1

АРф e-s 7 ^ . „0 , АРф e-s 7

f v2 J(s)eiv0dv + ^-^- f v2J'v(s)eiv0dv.

2

2а 8 + а -да а 82

Интегралы в последнем выражении могут быть вычислены по методу перевала с использованием асимптотического разложения функций Бесселя при 8 ^да в виде [8, 9]

Jv(s) =

а-

V 2 + s

2 I 2 2 v Vv + s - varcsh—

В рассматриваемом случае седловая точка v0 = is sin 9 . В связи с тем что функция, стоящая под знаком интеграла, имеет точки ветвления V = ±is, результат, получаемый по методу перевала, применим

к

для углов 9< —. После преобразования выражение для давления отражения на поверхности оболочки может быть записано в виде

P2 =

Р^х

R

ЭФ 1 ЭФ

3 ЭФ 3 ЭФ

Эг 2s Эг 8s2 Эг 8s3 Эг 2s2

3 2 1 2

-Ф sin2 9 cos 9-—Ф sin2 9 cos 9

г=1

Пользуясь теоремой о свертке, оригинал этого выражения можно представить в виде

P = -

рАф

R

Г^dTi +1 Г£Ф(T-Ti)dv

Эг 2 0 Эг 1 2 0 Эг 1

— f^i (Т-Ц)2 dT  1 Эф 9 cos 9 i sin2

16 J Эг

i +---- sin 9 cos 9-ф sin 9 cos 9

1 2 Эг

г=1

Применение асимптотической формы для больших интервалов времени, соизмеримых с периодом обтекания ударной волной оболочки, не представляется возможным, так как задание отношения функций Макдональда в виде (6) не дает предельного перехода при 5 ^ 0. Полученная формула (7) дает точную картину распределения давления на жестком цилиндре при т ^ 0. Кроме того, она позволяет оценить точность приближенных решений для давления при т ^ 0. В качестве такого приближенного решения можно рассматривать уравнение

Эу Эу у = Z Sn , Эг Эт г

где gn — параметр, учитывающий влияние течения за фронтом ударной волны.

Изображение для потенциала у в этом случае с учетом гранич-

1 ЭФ

ных и начальных условий имеет вид у =

s + Sn Эг

г=1

1

s

s

На основании исходного уравнения (8) и теоремы о свертке функций можно записать

Р2 =-

рДф

Я

5ф Еп

дг

дг

г=1

В этом случае при т^- 0 решение совпадает с асимптотическим, так как у удовлетворяет соответствующему предельному переходу при 8 ^ 0 в поле изображений.

Результаты численных расчетов по разработанной модели. На рис. 3 приведены графики зависимости давления во фронтальной точке (9 = 0) оболочки от параметра т, полученные по асимптотической формуле (7) и приближенной формуле (9). Здесь же для сравнения приведена кривая, соответствующая известному решению Блейха [10]. Как видно из представленных графиков для т < 2,0, оба решения практически совпадают, после чего начинается их расхождение. Для т > 2,0 можно пользоваться только формулой (7).

(Рх+Р2)1АР 2

Рис. 3. Зависимость давления во фронтальной точке оболочки от безразмерного момента времени:

1 — точное решение; 2 — решение при сохранении двух членов разложения; 3 — приближенное решение; 2 % 4 — приближенное решение Блейха

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

(Р,+Р2)/ЛР

1

0 1 2

Рис. 4. Зависимость давления в различных точках оболочки от безразмерного параметра времени:

1 — 9 = 0°; 2 — 9 = 60°; 3 — 9 = 80°; 4 — 9 = 90°; 5 — 9 = 100°; 6 — 9 = 120°; 7 — 9 = 140°; 8 — 9 = 160°; 9 — 9 = 180°

Полученные результаты справедливы для слабых ударных волн,

ЛРф

когда —— < 0,3, где Р0 — давление в невозмущенной среде перед

Р0

ЛРф

фронтом ударной волны. При -> 0,3 полученное давление следу-

Р0

ет умножать на коэффициент Котр.

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

поверхности оболочки

давление в ней мгновенно возраста-

ет до максимального значения P = P1 + P2, а затем падает до давления обтекания по мере продвижения фронта ударной волны в течение времени до полного погружения оболочки в ударную волну. Для точек тыльной поверхности характерно наличие времени нарастания давления до максимальной величины, вполне определенного для каждой точки этой поверхности оболочки. Время нарастания давления до максимума может быть определено по формуле

t = ^^cos9. На рис 5, 6 представлены характерные зависимости

АФ

Рис. 5. Характер изменения давления в точках передней поверхности оболочки при взаимодействии с ударной волной

изменения давления в произвольной точке передней и тыльной поверхностей оболочки в течение времени ее обтекания ударной волной. На рис. 7 представлен характерный график изменения максимального избыточного давления по периметру оболочки.

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

АР

где Po — давление в невозмущенной среде; AP(t) — местное давление в ударной волне; Pe — коэффициент давления; U(t) — скорость

При стационарном обтекании оболочки потоком несжимаемой воздушной массы изменение Pe может быть представлено в виде

pe = 1 - 4 sin2 e.

Как показано в [Ii], коэффициенты давления, вычисляемые по формуле (10) и полученные экспериментальным путем с помощью продувок в аэродинамических трубах для передней поверхности оболочки, достаточно хорошо согласуются. Для тыльной части поверхности наблюдается некоторое расхождение результатов, объясняемое отрывом струи воздушного потока от поверхности оболочки вблизи миделевого сечения и образованием вихрей. Тем не менее это расхождение, согласно [i2], не ведет к существенным погрешностям в определении коэффициента давления. Кроме того, для давления APe < 0,8 кг/см2 сжимаемость воздуха можно не учитывать при определении нагрузок на оболочку. При больших давлениях сжимаемость воздуха может быть учтена введением поправочного коэффициента в. Тогда скоростной напор с учетом сжимаемости воздуха

qi = q(1 + в).

Зависимость в от числа Маха приведена на рис. 8.

воздуха в ударной волне.

Рис. 7. Изменение максимального давления по периметру оболочки

При взаимодействии с ударной волной оболочка вследствие упругих деформаций начинает излучать волну. Обозначив через у и Р3 соответственно потенциал скоростей и давление в излученной волне, по аналогии с (1) и (2) будем иметь

Рз = -

рАф 5ф1

R Эх

r=1

2 М

Рис. 8. Зависимость поправочного коэффициента от числа Маха

Э 2Ф1

Эг 2

1 +_L .ЭФ!

r Эг r 2 эе2

э V

Эх2

Для отыскания у (г, 0, т) имеем граничные условия на поверхности оболочки вида

1 Э^1

Аф Эг

r=1

Э^

"Э?

где н,(г, 0, т) — прогиб оболочки, положительной внутрь.

Для определения Р3 используем асимптотическое представление для .

Проведя преобразования по Лапласу и Фурье, аналогичные преобразованиям для потенциала у(г, 0, т), получим изображение для давления излучения

Рз =

рАф

2 г

R

* 1 * 3 *

-sw + —w--w +---

2 8s 8s2

V * V *

-w +--w —-w I.

2s s2 J

Для последнего выражения легко может быть выполнено обратное преобразование Лапласа, однако обратное преобразование Фурье для двух последних слагаемых в общем случае при неизвестной функции w(e, х) выполнить нельзя. Формула (12) совместно с формулой (6) может быть использована в таком виде, если для решения уравнений движения оболочки преобразование Лапласа и Фурье известны.

Если прогиб w(e, х) определяется в виде разложения по периодическим функциям типа sin е и cos е, то может быть использован метод разделения переменных. Для форм колебаний, близких к осесим-

метричным, согласно [7], в формуле (12) можно пренебречь последними двумя слагаемыми. Тогда давление излучения можно записать в виде

Р* =

„ г\2 ( 1 х ~ х Л

рДф' 1 - -

Я

1 3 3

'(0, х) + 2* (0, х)- -1W (0, Х1)Х1 + -1* (0, х1)ёХ1

V

.

Если предположить, что каждый кольцевой элемент оболочки излучает цилиндрическую волну, то для ^(г, 0, х) можно приближенно записать уравнение типа (8). Как было отмечено ранее, это уравнение при gn = 0,5 есть асимптотическая форма волнового уравнения для осесимметричного движения. Решение этого уравнения при граничных условиях (11) может быть представлено в виде

Р3 = -М

3 Я

1 - 1 г /л \-2(х-х1) * — * + —I *(0;х1 )е 2 ах1

2 4 о

Последнее выражение можно представить в форме Р3 | * (в, х, )е "1 (х-х1 У х,

о

Если в асимптотическом разложении (6) использовать только один член, то для Р3 получим формулу

рАф*

Р3 = .

3 Я

Сравнение, приведенное для давления излучения, показывает, что формулу (15) можно использовать для малых х (х < 0,1). При больших х эта формула может привести к неверным результатам. То же самое верно для формул (12) и (13). Формула (14) справедлива для форм прогиба оболочки, близких к осесимметричным. Однако для больших интервалов взаимодействия оболочки с ударной волной осесимметричная форма колебаний играет основную роль, поэтому целесообразно использовать эту формулу для оценки давления излучения на поверхности оболочки. Для малых х функция прогибов сильно изменяется в окружном направлении и основную роль играют высшие гармоники колебаний, в связи с чем целесообразно использовать асимптотическую формулу (12).

Представляет интерес величина давления излучения по сравнению с давлением в ударной волне. На рис. 9 приведены графики изменения давления Р3, выраженного в процентах по отношению ДРф, для различных оболочек. Как следует из полученных графиков, отклонение давления излучения составляет на основном участке примерно 5 % от ДРф

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

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

• неравномерно по контуру поперечного сечения оболочки;

• мгновенно нарастает при соприкосновении с ударной волной для всех точек передней поверхности оболочки;

• имеет конечную скорость нарастания при погружении оболочки в ударную волну для всех точек тыльной поверхности оболочки.

Для процесса взаимодействия ударной волны с цилиндрической оболочкой характерны два периода: 1) период погружения — от момента встречи оболочки с ударной волной до момента полного погружения в нее; 2) период обтекания — от момента полного погружения оболочки в ударную волну до конца фазы движения.

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

ЛИТЕРАТУРА

[1] Карлов Н.В., Кириченко Н.А. Колебания, волны, структуры. Москва, Физматгиз, 2003, 496 с.

[2] Регирер Л. Ударные волны. Москва, Наука, 1979, 542 с.

[3] Алешков Ю.З. Теория взаимодействия волн с преградами. Ленинград, Изд-во ЛГУ, 1990, 371 с.

[4] Димитриенко Ю.И., Беленовская Ю.В., Анискович В.А. Численное моделирование ударно-волнового деформирования гибких броневых композитных материалов. Наука и образование: электронное научно-техническое издание, 2013, № 12. doi: 10.7463/1213.0665297

Рис. 9. Зависимость давления излучения от избыточного давления во фронте ударной волны

[5] Димитриенко Ю.И., Дзагания А.Ю., Беленовская Ю.В., Воронцова М.А. Численное моделирование проникания ударников в анизотропные упру-гопластические преграды. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2008, № 4, с. 100-117.

[6] Димитриенко Ю.И., Димитриенко И.Д. Моделирование динамических процессов деформирования гибких тканевых композиционных материалов. Инженерный журнал: наука и инновации, 2014, вып. 5(29). URL: http://engjournal.ru/search/author/40/page1.html

[7] Пиковский А., Розенблюм М., Куртс Ю.М., Синхронизация: фундаментальное нелинейное явление. Москва, Техносфера, 2003, 496 с.

[8] Кубенко В.Д. Проникновение упругих оболочек в сжимаемую жидкость. Киев, Наукова думка, 1984, 158 с.

[9] Димитриенко Ю.И. Механика сплошной среды. Т. 2. Универсальные законы механики и электродинамики сплошной среды. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2011, 560 с.

[10] Димитриенко Ю.И. Механика сплошной среды. Т. 4. Основы механики твердого тела. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2013, 624 с.

[11] Березин А.В., Жуков Д.А., Жуковский М.Е., Конюков В.В., Крайнюков В.И., Марков М.Б., Помазан Ю.В., Потапенко А.И. Моделирование электромагнитных эффектов в сложных конструкциях при воздействии импульсных излучений. Математическое моделирование и численные методы, 2015, № 2 (6), с. 58-72.

[12] Дубровин В.М., Бутина Т.А. Моделирование динамической устойчивости цилиндрической оболочки при действии осевой сжимающей нагрузки. Математическое моделирование и численные методы, 2015, № 2 (6), с. 46-57.

[13] Дубровин В.М., Бутина Т.А. Моделирование процесса взаимодействия ударной волны с цилиндрической оболочкой. Инженерный журнал: наука и инновации, 2014, вып. 4 (28). URL: http:// engjournal.ru/ catalog/ mathmodel/ solid/ 1233.html

[14] Бутина Т.А., Дубровин В.М. Устойчивость цилиндрической оболочки при комбинированном нагружении. Инженерный журнал: наука и инновации, 2012, вып. 2. URL: http://engjournal.ru/articles/44/44.pdf

[15] Дубровин В.М., Бутина Т.А. Моделирование динамической устойчивости цилиндрической оболочки при действии внешнего избыточного давления. Инженерный журнал: наука и инновации, 2014, вып. 6 (30). URL: http:// engjournal.ru/catalog/mathmodel/solid/1237.html

Статья поступила в редакцию 15.03.2015

Ссылку на эту статью просим оформлять следующим образом:

Дубровин В.М., Бутина Т.А. Моделирование процесса взаимодействия ударной волны с цилиндрической оболочкой с учетом волны излучения. Математическое моделирование и численные методы, 2015, № 4 (8), с. 38-52.

Дубровин Виктор Митрофанович родился в 1935 г., окончил Саратовский государственный университет в 1958 г. Канд. техн. наук, доцент МГТУ им. Н.Э. Баумана. Специалист в области прочности и устойчивости деформируемых систем. Автор пяти изобретений. e-mail: dubrovinvm1934@yandex.ru

Бутина Татьяна Александровна родилась в 1950 г., окончила Московский физико-технический институт в 1974 г. Канд. физ.-мат. наук, доцент МГТУ им. Н.Э. Баумана. Специалист в области прочности и устойчивости деформируемых систем. e-mail: butina_ta@mail.ru

Полякова Надежда Сергеевна родилась в 1950 г., окончила МГУ им. М.В. Ломоносова в 1972 г. Канд. физ.-мат. наук, доцент МГТУ им. Н.Э. Баумана. Специалист в области функционального анализа. e-mail: polyakova.nadiya@mail.ru

Modeling of the process of interaction of the shock wave with cylindrical shell considering wave radiation effect

© V.M. Dubrovin, T.A. Butina, N.S. Polyakova Bauman Moscow State Technical University, Moscow, 105005, Russia

The objective of this research is to examine the shock wave with cylindrical shell and to propose a method for calculating its dynamic stability under axial compressive time-varying load. For weak shock waves we conducted comparative analysis of the exact solution and the existing approximate solutions. We evaluated the wave radiation effect after the shell deformation. The case of linearly varying load was considered as an example.

Keywords: incident shock wave, reflected shock wave, radiated shock wave, overpressure, cylindrical shell, integral transformation, asymptotic representation.

REFERENCES

[1] Karlov N.V., Kirichenko N.A. Kolebaniya, volny, struktury [Vibrations, waves, structures]. Moscow, Fizmatlit Publ., 2011, 497 p.

[2] Regirer L.S. Udarnye volny [Shock waves]. Moscow, Nauka Publ., 1979, 542 p.

[3] Aleshkov Yu.Z. Teoriya vzaimodeystviya voln s pregradami [The theory of the interaction of waves and obstacles]. Leningrad, LGU Publ., 1990, 371 p.

[4] DimitrienkoYu.I., Belenovskaya Yu.V., Aniskovich V.A. Computational Modeling of shock-wave deforming of flexible armour composite materials. Nauka i obrazovanie: elektronnoe nauchno-tekhnicheskoe izdanie — Science and Education: Electronic Scientific and Technical Edition, 2013, no. 12. doi: 10.7463/1213.0665297

[5] Dimitrienko Yu.I., Dzaganiya A.Yu., Belenovskaya Yu.V., Vorontsova M.A. Vest-nik MGTU im. N.E. Baumana. Ser. Estestvennye nauki — Herald of Bauman Moscow State Technical University. Ser. Natural Science, 2008, no. 4, pp. 100-117.

[6] Dimitrienko Yu.I., Dimitrienko I.D. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2014, iss. 5 (29). Available at: http://engjournal.ru/search/author/40/page 1 .html

[7] Pikovsky A., Rozenblyum M., Kurts Yu.M. Sinkhronizatsiya: Fundamental-noe nelineynoe yavlenie [Synchronization: The fundamental nonlinear phenomenon]. Moscow, Tekhnosfera Publ., 2003, 496 p.

[8] Kubenko V.D. Proniknovenie uprugikh obolochek v szhimaemuyu zhidkost [Penetration of elastic shells in compressible fluid]. Kiev, Naukova dumka Publ., 1984, 158 p.

[9] DimitrienkoYu.I. Mekhanika sploshnoi sredy.Tom 2. Universalnye zakony mekhaniki i elektrodinamiki sploshnykh sred [Continuum mechanics. Vol. 2. General lows of the continuum mechanics and electrodynamics]. Moscow, BMSTU Publ., 2011, 590 p.

[10] DimitrienkoYu.I. Mekhanika sploshnoi sredy. Tom 4. Osnovy mekhaniki tver-dogo tela [Continuum mechanics. Vol. 4. Fundamentals of solid mechanics]. Moscow, BMSTU Publ., 2013, 624 p.

[11] Berezin A.V., Zhukov D.A., Zhukovskiy M.E., Konukov V.V., Krainukov V.I., Markov M.B., Pomazan Yu.V., Potapenko A.I. Matematicheskoe modeliro-vanie i chislennye metody — Mathematical Modeling and Computational Methods, 2015, no. 2 (6), pp. 58-72.

В.М. flyópoem, T.A. Eymma, H.C. nomKoea

[12] Dubrovin V.M., Butina T.A. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2015, no. 2 (6), pp. 46-57.

[13] Dubrovin V.M, Butina T.A. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2014, iss. 4 (28). Available at: http://engjournal.ru/catalog/mathmodel/solid/1233.html

[14] Dubrovin V.M., Butina T.A. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2012, iss. 2. Available at: http://engjournal.ru/articles/44/44.pdf

[15] Dubrovin V.M., Butina T.M. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2014, iss. 6 (30). Available at: http:// engj ournal.ru/catalog/mathmodel/solid/1237.html

Dubrovin V.M. (b. 1935) graduated from the Faculty of Mathematics and Mechanics, Saratov State University in 1958. Cand. Sci. (Eng.), Associated Professor of the Computational Mathematics and Mathematical Physics and of the Higher Mathematics Departments at Bauman Moscow State Technical University. Specialist in the field of strength, stability of deformable systems. Research interests: dynamics of strength and stability of deformable systems; creep of structural materials. e-mail: vmdubrovinvm1934@mail.ru

Butina T.A. (b. 1950) graduated from the Faculty of Management and Applied Mathematics of Moscow Institute of Physics and Technology in 1974. Cand. Sci. (Phys. & Math.), Associated Professor of the Computational Mathematics and Mathematical Physics Department at Bauman Moscow State Technical University. A specialist in the field of strength and stability of deformable systems. e-mail: butinata@mail.ru

Polyakova N.S. (b. 1950) graduated from Lomonosov Moscow State University in 1972. Cand. Sci. (Phys.-Math.), Assoc. Professor of the Computational Mathematics and Mathematical Physics Department at Bauman Moscow State Technical University. A specialist in the field of functional analysis. e-mail: polyakova.nadiya@mail.ru