========
/BEM/DAA
========

Ключевое слово блочного формата Двойная асимптотическая аппроксимация для подводной среды 
 взрыв, где матрица масс жидкости рассчитывается методом граничных элементов.

Формат
------

.. csv-table::
   :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
   :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

   "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID", "/BEM/DAA/daa_ID/unit_ID"
   "daa_title", "daa_title", "daa_title", "daa_title", "daa_title", "daa_title", "daa_title", "daa_title", "daa_title", "daa_title"
   "surf_ID", "grav_ID", "", "", "", "", "", "", "", ""
   ":math:`\rho`", ":math:`\rho`", "C", "C", "Пмин", "Пмин", "", "", "", ""
   "Xs", "Xs", "Ys", "Ys", "Zs", "Zs", "", "", "", ""
   "Яформа", "Ипри", "Ипре", "", "", "Кформ", "Фрисерфинг", "Послеток", "Интегрировать", ""

Вставьте, если
Ипре
=
1
.. csv-table::
   :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
   :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

   "Pm", "Pm", ":math:`\theta`", ":math:`\theta`", "a", "a", ":math:`a_{\theta}`", ":math:`a_{\theta}`", "", ""

Вставьте, если
Ипре
=
2
.. csv-table::
   :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
   :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

   "fct_IDP", "", "FscaleP", "FscaleP", "", "", "", "", "", ""

.. csv-table::
   :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
   :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

   "Xc", "Xc", "Yc", "Yc", "Zc", "Zc", "", "", "", ""

Вставьте, если
grav_ID
>
0
or
Фрисерфинг
=
2
.. csv-table::
   :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
   :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

   "XA", "XA", "YA", "YA", "ZA", "ZA", "", "", "", ""
   "Дир-Х", "Дир-Х", "Дир-Y", "Дир-Y", "Дир-З", "Дир-З", "", "", "", ""

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

.. csv-table::
   :header: "Поле", "Содержание", "Пример единицы СИ"
   :widths: 33, 33, 33

   "daa_ID", "блок ДАА   идентификатор.(Целое число, максимум 10 цифр)", ""
   "unit_ID", "(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)", ""
   "daa_title", "блок ДАА   заголовок.(Символ, максимум 100 символов)", ""
   "surf_ID", "Идентификатор мокрой поверхности.   2 3 (целое число)", ""
   "grav_ID", "/GRAV   option identifier.(Integer)", ""
   ":math:`\rho`", "жидкость   плотность.(Реальная)", ":math:`[\frac{kg}{m^{3}}]`"
   "C", "Плавный звук   скорость.(Реальная)", ":math:`[\frac{m}{s}]`"
   "Пмин", "Отсечка давления (<   0).По умолчанию = -1030 (Реальный)", ":math:`[Pa]`"
   "Xs", "X-координата противостояния   точка. 3(Реал)", ":math:`[m]`"
   "Ys", "Y-координата противостояния   точка. 3(Реал)", ":math:`[m]`"
   "Zs", "Z-координата противостояния   точка. 3(Реал)", ":math:`[m]`"
   "Яформа", "Флаг БЭМ-решения. = 1 (по умолчанию) Интегрирование Гаусса. = 2 Аналитическая интеграция. (Целое число)", ""
   "Ипри", "Уровень флага распечатки. = 1 (по умолчанию) Уменьшенная распечатка. = 2 Полная распечатка. (Целое число)", ""
   "Ипре", "Флаг нагрузки под давлением.   6 = 1 Давление рассчитывается как   :math:`P_{i}(t)=P_{m}e^{−\frac{t}{\theta}}`  . = 2 Ввод по функции. (Целое число)", ""
   "Кформ", "Флаг анализа. = 1 (по умолчанию) Состав DAA1 = 2 Высокая частота (Целое число)", ""
   "Фрисерфинг", "Флаг свободной поверхности. 6 = 1 (по умолчанию) Нет = 2 Да (Целое число)", ""
   "Послеток", "Расчет послетока.   7 = 1 Нет = 2 (по умолчанию) Да (Целое число)", ""
   "Интегрировать", "Целочисленный флаг времени. = 1  Первый заказ = 2 (по умолчанию) Второй заказ (Целое число)", ""
   "Pm", "Максимальное давление. 5(Реал)", ":math:`[Pa]`"
   ":math:`\theta`", "Распад   время.(Реальное)", ":math:`[s]`"
   "a", "Постоянное максимальное давление.   5(Реал)", ""
   ":math:`a_{\theta}`", "Время спада давления   постоянный. 5(Реал)", ""
   "fct_IDP", "Функция давления падения   идентификатор.(Целое число)", ""
   "FscaleP", "Шкала ординат (давления)   фактор для   fct_IDP.(Реальный)", ":math:`[Pa]`"
   "XC", "X-координата взрывчатого вещества   заряд.(Реальный)", ":math:`[m]`"
   "YC", "Y-координата взрывчатого вещества   заряд.(Реальный)", ":math:`[m]`"
   "ZC", "Z-координата взрывчатого вещества   заряд.(Реальный)", ":math:`[m]`"
   "XA", "Координата X точки A на   свободная поверхность.(Реальная)", ":math:`[m]`"
   "YA", "Координата Y точки А на   свободная поверхность.(Реальная)", ":math:`[m]`"
   "ZA", "Координата Z точки A на   свободная поверхность.(Реальная)", ":math:`[m]`"
   "Дир-Х", "X-компонент нормы   к плоскости свободной поверхности.(Реальная)", ""
   "Дир-Y", "Y-компонент нормы   к плоскости свободной поверхности.(Реальная)", ""
   "Дир-З", "Z-компонент нормы   к плоскости свободной поверхности.(Реальная)", ""

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

1. Вся конструкция должна быть 
 смоделировано. Симметричный анализ не поддерживается.
2. Поверхность нормальная
  :math:`n`
  должен быть направлен в жидкость.
3. Точка отступления определяется с помощью

                        (
  X
  s
  ,
  Y
  s
  ,
  Z
  s
  ) — это место, где

                    падающая волна давления задана в момент времени t=0:
  .. image:: images/bem_daa_starter_r_bem_daa.png
     *(Рисунок 1.)*
4. Плоскую волну можно смоделировать 
 используя сферическую волну и помещая заряд взрывчатого вещества достаточно далеко.
5. Давление в точке зазора как

                    функция времени:
  :math:`P_{i}(t)=P_{m}e^{−\frac{t}{\theta}}`
  Где,
  :math:`P_{m}`
  Максимальное давление
  :math:`t`
  Время
  :math:`\theta`
  Время затухания
  Максимальное давление и время затухания могут быть рассчитаны

                            используя:
  :math:`P_{m}=K[\frac{W^{\frac{1}{3}}}{R}]^{a}`
  :math:`\theta=K_{\theta}W^{\frac{1}{3}}[\frac{W^{\frac{1}{3}}}{R}]^{a_{\theta}}`
  :math:`W`
  Взрывоопасная масса.
  :math:`R`
  Расстояние до взрыва.
  :math:`K`
  ,
  :math:`\alpha`
  ,
  :math:`K_{\theta}`
  и
  :math:`a_{\theta}`
  Константы в зависимости от взрывчатки.
  If
  :math:`W`
  in kg,
  :math:`R`
  в метре,
  :math:`P_{m}`
  в МПа и в мс.
  .. csv-table::
     :header: "", ":math:`K`", ":math:`\alpha`", ":math:`K_{\theta}`", ":math:`a_{\theta}`"
     :widths: 20, 20, 20, 20, 20
     "TNT", "52.12", "1.180", "0.0895", "-0.185"
     "PETN", "56.21", "1.194", "0.0860", "-0.257"
     "HBX", "53.51", "1.144", "0.0920", "-0.247"
6. Свободная поверхность – это плоскость, определенная 
 по точке и ее вектору нормали.
7. Нормальная скорость послетока равна

                    рассчитывается как:
  :math:`v_{afterflow}=\frac{cos\gamma}{\rhoR}\int_{0}^{t}P(t)dt`
  P
  Жидкостная точка.
  C
  Точка взрывчатого заряда.
  S
  Точка противостояния.
  .. image:: images/bem_daa_starter_r_bem_daa_afterflow.png
     *(Рисунок 2.)*

1
Литтлвуд, Дж. де Рунц Т.

                2004.»
Кодекс США
". Мастерская Mecalog, София Антиполис, Франция,

                2004 г.
2
Коул, Роберт Х.
Подводный взрыв
. Издательство Принстонского университета,

                1948 год