Рис. 2. Зависимости изменения максимального избыточного давления от
расстояния до центра взрыва
длительности положительной фазы сжатия, рассчитываемых в точках
на жесткой поверхности с шагом 0,5 м на удалении 1. . . 5 м и шагом
1 м на удалении 5. . . 25 м.
Результаты расчетов для варианта 1 приведены на рис. 2, где пред-
ставлены максимальные избыточные давления
Δ
p
m
, рассчитанные на
удалении от центра взрыва
r
(кривая
2
) и определенные по формуле
Садовского (кривая
1
для заряда удвоенной массы) [8]:
Δ
p
m
= 0
,
084
3
√
2
m
0
r
+ 0
,
27
3
√
2
m
0
r
2
+ 0
,
7
3
√
2
m
0
r
3
.
(1)
Сравнение их показывает, что имеет место систематическое зани-
жение результатов численного расчета по сравнению с экспериментом.
На ближней границе применимости формулы (1)
r
=
3
√
2
m
0
5
,
85
м
расхождение между результатами составило 25%, а на расстоянии
20 м от центра взрыва уменьшилось до 3,1%, что говорит об удовле-
творительном соответствии расчета и эксперимента.
Аналогичные результаты наблюдаются при сравнении длительно-
стей и импульсов положительной фазы сжатия [8], где максимальное
отличие в результатах не превышает 25%.
Таким образом, верификация показала, что результаты численного
моделирования взрыва полусферического заряда на жесткой поверх-
ности удовлетворительно совпадают с результатами экспериментов.
На рис. 3 приведены распределения давления по координате вдоль
жесткой поверхности, рассчитанные в близкие моменты времени во
всех вариантах расчета, когда фронт УВ в воздухе удалился примерно
на 8 м от места взрыва.
При взрыве полусферического заряда (рис. 3,
а
) распределение да-
вления за фронтом воздушной УВ аналогично распределению, полу-
ченному при решении одномерной сферической задачи в момент вре-
33
1,2 4,5,6,7