Пиротехническая химия
Главная Начинающим пиротехникам Статьи Добавить статью Добавить материалы на сайт Поиск по сайту Карта книг Карта сайта
Книги в помощь
Военная история Изготовление и применение ВВ Пиротехника в военном деле Разное по пиротехнике Физика в пиротехнике Химия ВВ и составов
Новые книги
Яковлев Г.П. "122 мм самоходная пушка образца 1944 г." (Военное дело)

Суворов С. "Бронированная машина пехоты БМП -3 часть 1" (Военное дело)

Суарес Г. "Тактическое преимущество " (Военное дело)

Стодеревский И.Ю. "Автобиография записки офицера спецназа ГРУ " (Военное дело)

Соколов А.Н. "Альтернатива. Непостроенные корабли Российского императорского флота" (Военное дело)
Физика взрыва - Орленко Л.П.
Орленко Л.П. Физика взрыва. Под редакцией Орленко Л.П. — M.: ФИЗМАТЛИТ, 2002. — 832 c.
ISBN 5-9221-0219-2
Скачать (прямая ссылка): orlfizvzr2002.djvu
Предыдущая << 1 .. 233 234 235 236 237 238 < 239 > 240 241 242 243 244 245 .. 394 >> Следующая


Интегрирование системы (12.10)-(12.12) проводилось с использованием явной полностью консервативной разностной схемы с порядком аппроксимации 0(h2 + г) [12.21]. Однородность схемы достигается сквозным расчетом зоны контактного разрыва. При этом, для сохранения локальной аппроксимации при переходе через контактный разрыв в области ПД, использовалась логарифмическая сетка с размером ячеек, уменьшающимся к контактному разрыву. Расчет фронта УВ проводился с использованием характеристик системы (12.10 совместно с соотношениями динамической совместности [12.22], что обеспечивает выделение сильного разрыва и повышает точность расчета параметров на фронте. Для «размазывания» внутренних скачков использовалась искусственная вязкость, вводимая на основе анализа дифференциальных свойств решения [12.23]:

12.2. Взрыв зарядов конденсированного BB

485

где

О

Si

S2,

при

при

если

если

O2u

д2т д2и

д2т

dr — S3

dr-63

ди дг

ди

дг

^ О,

>0;

с — скорость звука; Si,S2 — коэффициенты при линейной и квадратичной составляющих псевдовязкого давления; q, S1, 62, S3 — коэффициенты дифференциального анализатора.

Наилучшие результаты получены при значениях констант Si = 0,25; Si = 1; S1 = S2 = 2; 03 = 3. Точность решения контролировалась проверкой интегрального закона сохранения энергии, который соблюдался с погрешностью не более 0,5 % в любой момент времени вплоть до завершения расчетов.

Расчеты проведены для

Таблица 12.2

Детонационные характеристики зарядов тэна при различной плотности

Рвв, кг/м3
1600
1200
800
400

D, м/с
7750
6040
4450
3020

QlO-6, Дж/кг
5,85
5,65
5,45
5,25

к
2,953
2,80
2,51
1,94

Го/Ге
0,0137
0,0153
0,0177
0,0226

воздуха с начальной температурой T0 = 15°С и зарядов тэна с плотностью рвв = 1600; 1200; 800; 400 кг/м3, вычисленные детонационные характеристики которых приведены в табл. 12.2 (Q — удельная теплота взрыва, А; — показатель адиабаты ПД на фронте ДВ).

На рис. 12.6-12.9 представлены результаты расчетов в

виде зависимостей от расстояния избыточного давления на фронте Арт, импульса

T Т„

избыточного давления і = J Ap dt, импульса скоростного потока j = J pu2dt

о о

и длительности фазы сжатия т в воздушной УВ для зарядов с плотностью

рвв = 1600кг/м3(сплошные линии) и рвв = 400кг/м3 (штриховые линии). Экспериментальные значения для Арт [12.11] нанесены на рис. 12.6 точечными линиями. Все результаты представлены в безразмерном виде относительно масштабов: рм = 101325 Па; рм = 1,2249KrZM3; им = (рм/рм)1^2 = 287,6 м/с; re = (mQ/рм)1/3; tM = ге/им', ім = Рм^м, где т — масса заряда. Значения безразмерного радиуса заряда приведены в последней строке табл. 12.2.

По избыточному давлению на фронте волны (рис. 12.6) результаты расчетов достаточно точно совпадают с экспериментальными данными, однако на расстояниях г/ге > 0,2.. .0,3 отчетливо просматривается тенденция к повышению Арт с уменьшением плотности заряда. Импульс скоростного потока (рис. 12.8) начинает возрастать при понижении плотности заряда на расстояниях г/ге > 0,09... 0,11, а импульс избыточного давления и длительность фазы сжатия (рис. 12.7, 12.9) возрастают как в ближней (г/ге < 0,1... 0,2), так и в дальней (г/ге > 0,3... 0,4) зонах взрыва. Немонотонный характер изменения і и т в диапазоне 0,05 < г/ге < 0,2 отражает сложный волновой характер течения среды, связанный с формированием и распространением вторичной УВ. Увеличение значения г во вторичной УВ выделено на рис. 12.7 заштрихованной областью. В дальней зоне взрыва при переходе от плотности рвв = 1600кг/м3 к рвв = 1400кг/м3 значения Apm,j,i увеличиваются соответственно на 10; 12 и 15 %. В табл. 12.3 для рассматриваемых вариантов приведены значения максимального радиуса пузыря продуктов взрыва

486

12. Взрыв в воздухе

Pu 1000

100

10

о



- 1
........3










0,01

0,1

rlr. 1

Рис. 12.6. Избыточное давление на фронте УВ при сферическом взрыве зарядов тэна: 1 — PBB = 1600 кг/м3; 2 — рвв = 400кг/м3; 3 — экспериментальные данные.

HL

0,3

0,1

0,03


-1
----2

Vl











0,01

0,1

rlr. 1

Рис. 12.7. Импульс избыточного давления в УВ при сферическом взрыве зарядов тэна: обозначения соответствуют рис. 12.6

HL

0,1

0,01

0,01


-1
----2








\\
\

0,1

rlr. 1

Рис. 12.8. Импульс скоростного потока в УВ при сферическом взрыве зарядов тэна: обозначения соответствуют рис. 12.6

0,1

0,01

-1
----2


— ___ N
А /і Ii
Ii
//
ж
\ //
\ її


«/ /

0,01

0,1

rlr, 1

Рис. 12.9. Длительность фазы сжатия в УВ при сферическом взрыве зарядов тэна: обозначения соответствуют рис. 12.6

rm , КПД взрыва г), под которым подразумевается доля выделившейся энергии, излучаемая в воздух при первой пульсации ПД, и температура воздуха на границе с ПД при распаде разрыва Ti и в момент достижения пузырем максимального радиуса T2.
Предыдущая << 1 .. 233 234 235 236 237 238 < 239 > 240 241 242 243 244 245 .. 394 >> Следующая
Реклама
 
 
Авторские права © 2010 PiroChem. Все права защищены.