Динамическое деформирование цилиндрических стеклопластиковых трубопроводов, заполненных жидкостью, при внутреннем импульсном нагружении

Динамика конструкций и сооружений

ДИНАМИЧЕСКОЕ ДЕФОРМИРОВАНИЕ ЦИЛИНДРИЧЕСКИХ

СТЕКЛОПЛАСТИКОВЫХ ТРУБОПРОВОДОВ, ЗАПОЛНЕННЫХ ЖИДКОСТЬЮ, ПРИ ВНУТРЕННЕМ ИМПУЛЬСНОМ НАГРУЖЕНИИ

А.Г. ШЕВЛЯКОВ, канд. техн. наук

Московский государственный университет природообустройства

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

КЛЮЧЕВЫЕ СЛОВА: тонкостенный цилиндрический трубопровод,

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

Начало исследованиям распространения волны давления в жидкости в цилиндрических трубопроводах заложено в классических работах [1, 2], в которых впервые на базе простейших математических моделей осмыслено явление гидравлического удара.

В монографиях [3, 4] исследованы движения потоков реальных жидкостей в трубах, моделируемых уравнениями тонкостенных оболочек. Общая теория деформирования оболочек, выполненных из композитных материалов, представлена в работах [5, 6]. Частные линейные задачи о распространении волновых процессов в тонкостенных оболочках с заполнителем (и без него) решены в статьях [7, 8, 9]. В книге [10] разработаны надежные методы исследования статики и динамики композитных тонкостенных элементов, склонных к потере устойчивости и испытывающих значительные прогибы в условиях статического и динамического нагружения. Предлагаемые математические модели и методы их исследования приложены к широким классам практически важных задач. Прочность цилиндрических оболочек из композитных материалов без заполнителя и с упругим заполнителем исследовалась авторами работ [11, 12]. Динамика тонкостенной изотропной цилиндрической оболочки большой протяженности, нагруженной давлением слабой воздушной ударной волны, движущейся по её наружной поверхности, в линейной постановке рассмотрена авторами [13]. В работе [14] предложена линейная одномерная модель динамического деформирования гидроупругих цилиндрических и сферических систем, находящихся под действием внутренних источников, расположенных на оси симметрии или в центре сферической оболочки. Внутренний объем композитной конструкции полностью заполнен идеальной сжимаемой жидкостью. В монографии [15] разработаны важные для практики линейные модели колебаний замкнутых тонко-

стенных оболочек, частично заполненных идеальной несжимаемой жидкостью. Исследованы случаи свободных и вынужденных колебаний гидроупругих систем. В рамках различных теорий тонкостенных цилиндрических оболочек исследованы процессы их деформирования, обусловленные внешним кратковременным давлением [16, 19]. Учтены анизотропные и вязкоупругие свойства материалов цилиндрических оболочек. Прикладные теории и расчеты трубопроводов и трубопроводных систем исследованы в работах [20-23]. Статьи [25-28] посвящены экспериментальным исследованиям параметров воздушных ударных волн и волновым процессам, протекающим при действии подземных ядерных взрывов. Рассмотрены обобщенные оценки действия взрывных нагрузок. Описаны способы создания взрывобезопасных объектов.

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

На первом этапе построено решение нелинейной сферически симметричной задачи о распространении одномерной волны в окрестности заряда.

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

Принимается, что источник динамического давления имеет первоначально сферическую форму радиуса г0 и помещен на оси гидроупругой системы (рис.1, а), причем Г) << R ^ - радиус трубопровода).

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

Пусть в начальный момент времени t = 0 происходит детонация химического ВВ. Следовательно, при t > 0 внутри сферы радиуса г0 происходит фор-

мирование газовой полости, радиус которой r = a(t) возрастает с течением времени. Вначале пренебрегаем сжимаемостью жидкости. В сферических координатах r, р, а это условие записывается в форме известного равенства div(v ) = 0 или в развернутом виде

Ч^-Ь^р s.np) + -^ = 0

r2 or r sin р др r sin р да

(r > a(t), 0 < р < л, 0 < а < 2л).

Здесь vr, vр, va - проекции вектора скорости жидких частиц на оси сферической системы координат.

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

В этом частном случае уравнение несжимаемости (1) значительно упрощается и записывается так

~

r2 dr

Это уравнение допускает решение в замкнутом виде

v

где C (t ) - произвольная функция времени (индекс r у скорости v опущен).

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

p

Здесь у - показатель адиабаты газа (постоянная величина); Р0 - начальное значение взрывного давления, r0- радиус сферического заряда в начальный момент времени t = 0 .

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

v

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

C

Поле скоростей частиц в сферическом случае выражается конечной зависимостью

v

Поле давлений в сферически симметричном случае описывается нелинейным дифференциальным уравнением первого порядка

Ov Ov 1 dp

— + v— = , r > a

dt dr p dr

После подстановки сюда выражения (7) и последующего интегрирования по r имеем следующее нелинейное дифференциальное уравнение второго порядка

-2

г г р

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

аа + 3 а2 = 1 р0 Г^У.

2 р \а)

Произвольная функция времени, появляющаяся в процессе интегрирования, полагается равной нулю. После введения обозначения с2 = р0/р это уравнение можно записать более компактно

аа +3а2 = с2^ У

Для концентрированных продуктов взрыва показатель адиабаты газа у чаще всего выбирают равным 4/3, поэтому предыдущее уравнение записывается так

а 2 = с 2 (Г014

аа + -а = с* .

Это уравнение (12) при начальных условиях а(?)| = г0, а(?)| í=0 = у0 проинтегрировано численно. Следовательно, известно давление р(г, ^) во всех точках сферы, радиус которой равен радиусу R цилиндрического трубопровода. Максимальное значение этого давления может достигать величин (30 - 40) кПа.

С этого момента начинается второй этап решения задачи. Он характеризуется тем, что во внутренней области трубопровода формируется гидродинамическая волна сжатия, которая распространяется вдоль оси конструкции.

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

Отнесем движение гидроупругой системы к цилиндрическим координатам Г[, г и в с началом в центре торца (г = 0). Для описания осесимметричных деформаций цилиндрического трубопровода используется линейная система динамических уравнений ортотропной оболочки [12] в стандартных обозначениях

„ , 82и h 8w , 82и

В^—- + --= р0 h—— - q1,

1 8г2 12 R 82 Ио 82 1

„ h 8и „ h „ h3 84 w , 82 w ---+ В2 —- w + В,--— = -р0 h—-— q2.

12 2 2 1 4 0 2 2

R 8г 2 R2 1 12 8г4 0 812 Для ортотропного материала независимые параметры анизотропии В1, В2 и ^12 выражаются через технические постоянные по следующим формулам

В: = ~Е~, В2 =, 612 =

1 -Уу2 1 -Уу2 1 -Уу2

Здесь у2- модуль нормальной упругости при растяжении в направлении основы; Е2 - в направлении утка; у1 - коэффициент Пуассона, характеризующий сокращение в направлении основы при растяжении в направлении утка; у2 - коэффициент Пуассона, характеризующий сокращение в направлении утка

при растяжении в направлении основы; р 0 - плотность материала трубопровода; h - толщина стенки трубопровода.

Радиальный прогиб w( г, t) считается положительным, если он направлен к центру кривизны; и( г, t)- осевое перемещение точек трубопровода. Для металлического трубопровода соотношения (14) принимают вид

В = В 2 , G =

1 -V 1 -V

Через ^ и обозначены поверхностные гидродинамические силы, действующие на внутреннюю поверхность конструкции. Заполняющую трубопровод жидкость считаем идеальной (невязкой), следовательно ^ = 0. Гидродинамическое давление в линейном приближении выражается зависимостью

д2 =-рж дF/Ыг= к,

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

д2 F 1 дF д2 F 1 д2 F 1 д2 Р

+ -—+ + = --Г,

дт^ т[ дг1 дг2 г1 дв2 сЖ дг2 1

Здесь сж - скорость гидродинамической волны. В случае осесимметрично-го движения последнее слагаемое слева в уравнении (17) обращается в нуль тождественно. Функция Р(т[,г,в,г) позволяет выразить вектор скорости частиц сжимаемой жидкости по известной зависимости

у^, Уг, vв) = gтadF,

из которой проекции вектора скорости на оси цилиндрической системы координат выражаются по следующим формулам

дР 1 дР дР

Уп =-т-, Ув = , Уг .

1 дт1 Т1 дв дг

В случае осесимметричной задачи компонента у в тождественно обращается в нуль (у в = 0) . Радиальная компонента скорости у^ в точках оси симметрии (т = 0, 0 < г < ж) равна нулю (ут (0, г, г) = 0), а продольная Уг зависит только от переменных т1 и г.

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

I дР ^

Ут\ =— =—, 0 < г <ж, г > 0.

Т1 !т1 =R дт „ дг

1 т1 =R

Если жидкость, заполняющая внутренний объем трубопровода, предполагается вязкой, то на внутренней поверхности (т1 = Я, 0 < г < ж) следует дополнительно использовать граничное условие

дР

= —, 0 < г <ж,

дг

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

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

В граничных точках торца (z = 0, 0 < r1 < R) действует давление p(t), вычисленное на первом этапе данной задачи

8F

- Рж — = P

8t

В точках контура торца (z = 0, r1 = R, 0 < в < 2л) выполняются граничные условия симметрии задачи, в соответствии с которыми при z = 0 в нуль обращаются 8w/ 8z и 83 w / 8z3 (производная от прогиба и поперечная сила).

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

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

В качестве примера рассмотрен трубопровод, выполненный из стеклопластика СВАМ-ИММ со следующими упругими характеристиками в главных направлениях упругости [5]

Ej = 34 • 106 кПа, E2 = 8,86 • 106 кПа, G12 = 2,25 • 106 кПа, v1 = 0,3 , v2 = 0,078 .

Трубопровод, изготовленный из композитного материала, заполнен водой с рж = 1000 кг/м3, сж = 1500 м/с. Толщина стенки трубопровода h = 0,01 м; радиус срединной поверхности R = 0,3 м; радиус заряда r0 равен 0,02 м; показатель адиабаты продуктов взрыва у = 1,4; их плотность р = 1167 кг/м3; плотность стеклопластика р0 = 2200 кг/м3; начальное давление продуктов взрыва p0 = 4-103 кПа. Начальная скорость при взрывном процессе составляет 200 м/с. Выбраны

следующие параметры расчетной сетки: Ar1 = Az = 10 3 м, At = 7,5 • 10 6 c .

Кривыми 1, 2 и 3 на рис.2 показано изменение во времени кратковременного импульса избыточного гидродинамического давления на внутреннюю поверхность трубопровода в трех точках наблюдения z = 3 см; z = 6 см и z = 12 см, отстоящих от торца z = 0. Анализ численных результатов показывает, что пиковое значение давления на стенке достигается в момент подхода в заданную точку поверхности фронта гидродинамической волны сжатия. Максимум давления на стенку трубопровода по мере роста расстояния точки z = const от торца z = 0 падает, а время достижения пикового значения давления растет. Импульс давления растягивается во времени. Зависимость давления на внутреннюю поверхность трубопровода от времени в точке наблюдения z = 20 см показана кривой на рисунке 3. Давление в этой точке носит осциллирующий характер; по-прежнему, ведущим фактором, определяющим давление на внутреннюю стенку трубопровода, является распространяющаяся вдоль трубы гидродинамическая волна. Вместе с тем, на этом графике уже видно заметное влияние упругих поперечных колебаний стенки цилиндрического трубопровода. Небольшие амплитуды отрицательных значений избыточного давления могут вызвать отрыв гидравлического потока от внутренней поверхности трубопровода.

Зависимость давления в точке г = 15 см на оси г = 0 от времени (рис.4) носит пульсирующий характер. Амплитуда первой пульсации связана с приходом гидродинамической волны в точку; последующие пульсации давления обусловлены упругими поперечными колебаниями стенки трубопровода.

Изменение во времени кольцевой деформации ед = конструкции в сечении г = 15см графически представлено кривой на рис.5. Ход этого графика показывает, что деформации растяжения в этом сечении чередуются с деформациями сжатия.

Из анализа числовых результатов можно сделать общий вывод:

1) главный вклад в динамическое нагружение внутренней поверхности трубопровода вносит гидродинамическая волна сжатия, распространяющаяся вдоль его оси;

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

3) возможен отрыв (отход) потока жидкости от внутренней поверхности конструкции.

1,5

1,0

0,5

р • 10 ~3, кПа

^2

^-3

)) ч___ t, мкс

0

40

120

Рис. 2. Изменение во времени избыточного гидродинамического давления в трёх точках наблюдения ( г = 3 см; 6 см и 12 см) на внутренней поверхности трубопровода

Рис. 3. Характер изменения во времени избыточного гидродинамического давления в точке z = 20 см на внутренней поверхности трубопровода

вода в точке z = 15 см

Рис. 5. Изменение во времени кольцевой деформации трубопровода в сечении z = 15 см

Если внутри трубопровода или емкости взрывается мощное ВВ (гексоген, пластит и др.), то на внутреннюю поверхность действуют интенсивные кратковременные нагрузки. Они могут разрушить трубопровод (емкость), что нередко приводит к разгерметизации конструкции [25-28]. В статье [29] представлена зависимость G = £ркрак^ для расхода газа через отверстие площадью S, образовавшееся в цилиндрической трубе. Здесь ркр и акр - критические параметры вязкого газа, постоянный коэффициент (0<£ < 1), G - секундный расход газа, измеряемый в кг/с. В случае идеального газа £ = 1.

В качестве критерия прочности стенки трубопровода можно принять условие |стд|тах <[ст], где [ст]- допускаемое напряжение материала трубопровода на

растяжение-сжатие, <Уд - окружное напряжение. Кроме того, должно выполняться условие устойчивости трубопровода (как цилиндрической оболочки) |ст|кр < икр , где <Укр - критическое напряжение, причем &кр = 0,18Eh / R .

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

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

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

Широким фронтом такие работы выполняются по Государственной научно-технической программе РФ "Безопасность населения и народно- хозяйственных объектов с учетом риска возникновения природных и техногенных катастроф". Настоящая работа выполнена в рамках Федеральной целевой научно-технической программы (ФЦНТП) "Исследования и разработки по приоритетным направлениям развития науки и техники", рассчитанной на 2002-2006 г.г.

Л и т е р а т у р а

1. Громека И.С. О скорости распространения волнообразного движения жидкости в упругих трубах. - М.: Изд-во АН СССР, Сборник сочинений, 1949. - 458 с.

2. Жуковский Н.Е. О гидравлическом ударе в водопроводных трубах. - М-Л.: Гос-техиздат, 1949. - 236 с.

3. Чарный И.А. Неустановившееся движение реальной жидкости в трубах. - М.: Недра, 1956. - 382 с.

4. Вольмир А. С. Оболочки в потоке жидкости и газа. Задачи аэроупругости. - М.: Наука, 1976. - 416 с.

5. Амбарцумян С.А. Общая теория анизотропных оболочек.-М.: Наука, 1974.-446 с.

6. Механика композитов в 12 томах. Под общей редакцией академика Гузя А.Н. Том 12, Прикладные исследования. - Киев: А.С.К., 2003 - 399 с.

7. Сабодаш П.Ф., Навал И.К. Осесимметричные волны в ортотропной цилиндрической оболочке// Механика композитных материалов, Рига, 1981, №5. - С.45-51.

8. Сабодаш П. Ф. Распространение волн в ортотропной оболочке вращения конечной длины// Известия АН СССР, Механика твердого тела, 1985, №2. - С.163-170.

9. Васильев В.В., Сабодаш П.Ф. Новые наукоемкие технологии в технике// Энциклопедия, том 5. Композиционные материалы. Расчет и конструирование. Российское космическое агентство. - М.: АО НИИ "ЭНЦИТЕХ", 1996. - 400 с.

10. Васильев В.В., Елпатьевский А.Н. Прочность цилиндрических оболочек из армированных материалов. - М.: Машиностроение, 1972. - 166 с.

11. Лихнов П.П. Динамика системы цилиндр-оболочка. - М.: Машиностроение, 1988. - 151 с.

12. Ляховенко И.А., Шенк Ю.В. Цилиндрическая оболочка при действии распространяющегося по её поверхности давления// Труды ЦАГИ, вып.2495, 1992.

13. Ляховенко И.А., Меркурьев А.В. Гидроупругое нагружение емкостей из метал-ло-органо-пластика// Труды ЦАГИ, вып.2646, 2001. - С.36-49.

14. Рапопорт И.М. Колебания упругой оболочки, частично заполненной жидкостью. - М.: Машиностроение, 1967. - 358 с.

15. Бивин Ю.К., Найда А.А. Несущая способность цилиндрических оболочек при воздействии динамического внешнего давления// Прикладная механика, 1970, вып.6, №10. - С.28-34.

16. Буштырков А.А., Найда А.А. Динамическое поведение тонкостенных цилиндрических оболочек из анизотропных упругих материалов при кратковременном внешнем воздействии// Будiвельник, 1973, вып.21. - С. 37-44.

17. Галли М., Гьюни М., Рьетор Е. Исследование переходных процессов при колебаниях давления в вязкоупругих трубах// Теоретические основы инженерных расчетов (перевод с английского), 1979, том 101, №4. - С. 58-65.

18. Блинштейн Ю.М., Хубларян М.Г. Распространение волны возмущений в бесконечно длинном вязкоупругом трубопроводе// Изв. АН СССР, Механика жидкости и газа, 1975, №3. - С. 145-158.

19. Аксельрод Э.Л., Ильин В.П. Расчет трубопроводов. - Ленинград: Машиностроение, 1972. - 322 с.

20. Федоров А.В., Фомин В.М. Математическое моделирование работы магистральных трубопроводов. - Изд-во ИТПМ СО АН СССР, Препринт, №35, 1983. - 153 с.

21. Овчинников В.Ф., Прохорович В.А., Смирнов Л.В. Методика численного расчета свободных колебаний многопролетных разветвленных трубопроводов// Прикладные проблемы прочности и пластичности. Методы решения. Всесоюзный Межвузовский Сборник, Горьковский университет, 1988. - С. 104-121.

22. Исследование технологии изготовления и монтажа резервуаров и трубопроводов. - М.: Стройиздат, 1987. - 144 с.

23. Гвоздев А.А. К расчету конструкций на действие взрывной волны// Строительная промышленность, 1943, №1-2. - С. 34-44.

24. Содовский М.А., Похил П.Ф. Импульс взрыва и его зависимость от формы и размеров заряда, а также от свойств взрывчатого вещества (ВВ)// Физика взрыва. М., 1952, Изд-во АН СССР, Сборник №1. - С. 89-95.

25. Физика ядерного взрыва, том 2, Действие взрыва. - М., МО РФ, Центральный физико-технический институт, Наука, 1997. - 253 с.

26. Вейтман Я.И., Смолин Н.И. Ударные воздушные волны промышленных взрывов. - М.: Недра, 1981. - 150 с.

27. Взрывные явления. Оценка и последствия. Книга 2. Перевод с англ. под ред. акад. Зельдовича Я.Б. и д-ра физ.-мат. наук Гельфанда Б.Е. - М.: Мир, 1986. - 382 с.

28. Демченко В.Г., Демченко Г.В. Исследование процесса продольного разрушения газопроводов. - М.: Транспорт и подземное хранение газа. Стройиздат, 1994. - 156 с.

29. Епишкин Ю.А., Клепацный А.Р. и др. Истечение струи вязкого газа через отверстие малого диаметра или трещину// Космонавтика и ракетостроение, ЦНИИМаш, 2002, вып.23. - С. 5-12.

30. Сидоров В.И. Исследования по промышленной безопасности в рамках ГНТП "Безопасность"// Проблемы безопасности и чрезвычайных ситуаций. Информационный сборник. - М.: ВИНИТИ, №4, 2003. - С. 62-68.

31. Кедринский В.К. Волновые процессы и динамика структуры неоднородных сред при импульсивном нагружении// ПМТФ, 1997, №4. - С. 111-139.

DYNAMIC DEFORMING CYLINDRICAL GLASS-AND-PLASTIC PIPES, FILLED BY LIQUID, UNDER INNER IMPULSE LOADING

The manuscript is devoted to the working out of non-classical model of dynamical behavior of thin walled cylindrical pipe containing compressive liquid under action of blast pressure. A chemical explosive is situated on the axis of the pipe.

Shevlyakov A.G.