=========
/EOS/PUFF
=========

Ключевое слово формата блока описывает 
 линейное полиномиальное уравнение состояния 
 :math:`P(\rho,E)` 
.

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

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

   "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID", "/EOS/PUFF/mat_ID/unit_ID"
   "eos_title", "eos_title", "eos_title", "eos_title", "eos_title", "eos_title", "eos_title", "eos_title", "eos_title", "eos_title"
   "C1", "C1", "C2", "C2", "C3", "C3", ":math:`\gamma_{0}`", ":math:`\gamma_{0}`", "", ""
   "T1", "T1", "T2", "T2", "Es", "Es", "", "", "", ""
   "H", "H", "E0", "E0", "", "", "", "", "", ""

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

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

   "mat_ID", "Идентификатор материала.(Целое число, максимум   10 цифр)", ""
   "unit_ID", "Идентификатор объекта.(Целое число, максимум 10 цифр)", ""
   "eos_title", "Титул EOS.(Символ, максимум 100   персонажи)", ""
   "C1", "Коэффициент С1. 1(Реал)", ":math:`[Pa]`"
   "C2", "Коэффициент С2.(Реал)", ":math:`[Pa]`"
   "C3", "Коэффициент С3.(Реал)", ":math:`[Pa]`"
   ":math:`\gamma_{0}`", "Грюнейзен гамма.(Реал)", ""
   "T1", "Т1   коэффициент. 1По умолчанию = C1 (Реальный)", ":math:`[Pa]`"
   "T2", "Т2   коэффициент.(Реальный)", ":math:`[Pa]`"
   "Es", "Энергия сублимации на единицу исходной   объем.(Реальный)", ":math:`[\frac{J}{m^{3}}]`"
   "H", "Коэффициент Н. 1(Реал)", ""
   "E0", "Начальная энергия на единицу начальной   объем.(Реальный)", ":math:`[\frac{J}{m^{3}}]`"

Пример (медь)
-------------

.. code-block::

   #RADIOSS STARTER

   #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

   /UNIT/1

   unit for mat

                      g                  cm                 mus

   #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

   #-  2. MATERIALS:

   #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

   /MAT/HYD_JCOOK/1/1

   Copper

   #              RHO_I               RHO_0

                    8.9                   0

   #                  E                  nu

                    1.3                 .34

   #                  A                   B                   n              epsmax              sigmax

                   9E-4              .00292                 .31                   0                   0

   #               Pmin

                      0

   #                  C           EPS_DOT_0                   M               Tmelt                Tmax

                   .025                1E-6                1.09                1356                1E30

   #              RHOCP                                                         T_r

               3.432E-5                                                           0

   #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

   /EOS/PUFF/1/1

   Copper

   #                 C1                  C2                  C3                  G0

                 1.3597              2.7466              2.4479                1.97

   #                 T1                  T2                  ES

                      0                   0             0.12282

   #                  H                  E0               

                      1             1.25675                   

   #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

   #ENDDATA

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

1. Позволять
  :math:`\mu=\frac{\rho}{\rho_{0}}−1`
  ,
  :math:`\eta=1+\mu`
  ,
  :math:`x=1−\frac{\rho_{0}}{\rho}`
  и
  E
  , внутренняя энергия единицы начального объема,

          давление определяется:
  Регион 1: 
 :math:`\mu\ge0`
  :math:`P=(C_{1}\mu+C_{2}\mu^{2}+C_{3}\mu^{3})(1−\frac{\gamma\mu}{2})+\gamma(1+\mu)E`
  с 
 :math:`\gamma\rho=\gamma_{0}\rho_{0}`
  Регион 2: 
 :math:`\mu<0` 
 и 
 :math:`E\geE_{s}`
  :math:`P=(T_{1}\mu+T_{2}\mu^{2})(1−\frac{\gamma\mu}{2})+\gamma(1+\mu)E`
  Регион 3: 
 :math:`\mu<0` 
 и 
 :math:`E\geE_{s}`
  :math:`P=\eta[H+(\gamma_{0}−H)\sqrt{\eta}][E−E_{s}(1−exp(\frac{N(\eta−1)}{\eta^{2}}))]`
  с 
 :math:`N=\frac{C_{1}\eta}{\gamma_{0}E_{s}}`
2. Уравнения состояния используются
  Радиосс
  для расчета гидродинамического давления и являются

          совместимые с материальными законами:
  - /MAT/LAW2 (PLAS_JOHNS)
  - /MAT/LAW3 (HYDPLA)
  - /MAT/LAW4 (HYD_JCOOK)
  - /MAT/LAW6 (HYDRO or HYD_VISC)
  - /MAT/LAW10 (DPRAG1)
  - /MAT/LAW12 (3D_COMP)
  - /MAT/LAW36 (PLAS_TAB)
  - /MAT/LAW44 (COWPER)
  - /MAT/LAW49 (STEINB)
  - /MAT/LAW102 (DPRAG2)
  - /MAT/LAW103 (HENSEL-SPITTEL)
  - /MAT/LAW109
  - /MAT/LAW133 (GRANULAR)

1
«66 компьютерных программ PUFF»

        Броди, Хормут, Лаборатория вооружения ВВС, май 1966 г.