/МАТ/ЗАКОН5 (JWL)

Radioss 2025.1

Ключевое слово в блочном формате

Этот закон описывает уравнение состояния Джонса-Уилкинса-Ли (JWL) для продуктов детонации взрывчатых веществ. Доступно опциональное моделирование дожигания.

Формат

/MAT/LAW5/mat_ID/unit_ID
или
/MAT/JWL/mat_ID/unit_ID
mat_title
ρ_i ρ_0 A B R1 R2 ω D P_CJ E_0 E_add I_BFRAC Q_OPT P_0 P_sh B_unreacted

Вставка, если E_add > 0 и Q_OPT = 0,1,2 (временное управление дожиганием):

T_start T_stop

Вставка, если E_add и Q_OPT = 3 (расширение Миллера):

a m n

Определение

Поле | Содержание ———————–|———————————– mat_ID | Идентификатор материала. (Целое число, максимум 10 цифр) unit_ID | Идентификатор единицы. (Целое число, максимум 10 цифр) mat_title | Название материала. (Символы, максимум 100 символов) ρ_i | Начальная плотность. (Реальное число) [кг/м³] ρ_0 | Эталонная плотность, используемая в уравнении состояния. По умолчанию ρ_0 = ρ_i. (Реальное число) [кг/м³] A | Параметр A уравнения состояния. (Реальное число) [Па] B | Параметр B уравнения состояния. (Реальное число) [Па] R1 | Параметр R1 уравнения состояния. (Реальное число) R2 | Параметр R2 уравнения состояния. (Реальное число) ω | Параметр ω уравнения состояния. (Реальное число) D | Скорость детонации. (Реальное число) [м/с] P_CJ | Давление Чепмена-Жуге. (Реальное число) [Па] E_0 | Энергия детонации на единицу объема. (Реальное число) [Дж/м³] E_add | Дополнительная энергия на единицу объема. = 0 — параметры дожигания не действуют. (Реальное число) [Дж/м³] I_BFRAC | Флаг вычисления доли сгорания. 0 = Объемное сжатие + Время горения. 1 = Только объемное сжатие. 2 = Только время горения. (Целое число) Q_OPT | Опциональная модель дожигания (если E_add > 0). 0 = Мгновенный выброс в T_start. 1 = Постоянная скорость от T_start до T_stop. 2 = Линейная скорость от T_start до T_stop. 3 = Расширение Миллера. (Целое число) P_0 | Начальное давление. (Реальное число) [Па] P_sh | Сдвиг давления. (Реальное число) [Па] B_unreacted | Неотреагировавший модуль сдвига взрывчатого вещества. (Реальное число) [Па] T_start | Время начала для дополнительной энергии (Q_OPT = 0, 1, 2). (Реальное число) [с] T_stop | Время окончания для дополнительной энергии (Q_OPT = 0, 1, 2). (Реальное число) [с] a | Опциональный параметр Миллера, если Q_OPT = 3. (Реальное число) [с⁻¹Па⁻ⁿ] m | Опциональный параметр Миллера, если Q_OPT = 3. (Реальное число) n | Опциональный параметр Миллера, если Q_OPT = 3. (Реальное число)

Пример (TNT)

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/JWL/2/123
TNT - данные из примера 46 с единицей: (г-см-мкс) - Стандартный JWL, Без дожигания
#              RHO_I
                1.63
#                  A                   B                  R1                  R2               OMEGA
              3.7121               .0323                4.15                 .95                  .3
#                  D                P_CJ                  E0                Eadd   I_BFRAC     Q_OPT
                .693                 .21                 .07                   0         0         0
#                 P0                 Psh          Bunreacted
                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/123
Система единиц расширения Миллера
                   г                  см                 мкс

Комментарии

Давление JWL выглядит как:

\[P_{jwl} = A(1 - \omega R_1 V)e^{-R_1 V} + B(1 - \omega R_2 V)e^{-R_2 V} + \omega(E + Q)V\]

Radioss выводит:

\[P = B_{frac} \cdot P_{jwl} + (1 - B_{frac})(P_0 + B_{unreacted} \cdot \mu) - P_{sh}\]

Где \(V = \frac{V}{V_0}\) — относительный объем, \(E = E_{int}/V_0\) — внутренняя энергия на единицу начального объема, \(\omega = \gamma - 1\), с \(\gamma = C_p/C_V\) (адвокатика) и обычное \(\mu = \rho/\rho_0 - 1\).

Для более подробной информации, пожалуйста, обратитесь к документации по материалу и моделям отказа в Reference Guide.