/FAIL/WIERZBICKI
================

Radioss - 2025.1

Описание
--------

Ключевое слово блока описывает модель отказа BAO-XUE-Wierzbicki.

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

.. code-block::

    /FAIL/WIERZBICKI/mat_ID/unit_ID C1 C2 C3 C4 m n I_fail_sh I_fail_so I_moy
    fail_ID

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

+----------+-----------------------------------------------------------+
| Поле     | Содержание                                                |
+==========+===========================================================+
| mat_ID   | Идентификатор материала. (Целое число, максимум 10 цифр) |
+----------+-----------------------------------------------------------+
| unit_ID  | Идентификатор единицы. (Целое число, максимум 10 цифр)   |
+----------+-----------------------------------------------------------+
| C1       | Первый параметр. (Вещественное число)                    |
+----------+-----------------------------------------------------------+
| C2       | Второй параметр. (Вещественное число)                    |
+----------+-----------------------------------------------------------+
| C3       | Третий параметр. (Вещественное число)                    |
+----------+-----------------------------------------------------------+
| C4       | Четвертый параметр. (Вещественное число)                 |
+----------+-----------------------------------------------------------+
| m        | Пятый параметр. (Вещественное число)                     |
+----------+-----------------------------------------------------------+
| n        | Показатель упрочнения. (Вещественное число)              |
+----------+-----------------------------------------------------------+
| I_fail_sh| Флаг модели отказа оболочки.                             |
|          | 1 - Оболочка удаляется, если D >= 1 для любой точки      |
|          | интеграции или слоя.                                     |
|          | 2 - Для каждой точки интеграции тензор напряжений        |
|          | обнулен, если D >= 1, и оболочка удаляется, если D >= 1  |
|          | для всех точек интеграции или слоя.                      |
|          | (Целое число)                                            |
+----------+-----------------------------------------------------------+
| I_fail_so| Флаг модели отказа твердого элемента.                    |
|          | 1 - Твердый элемент удаляется, если D >= 1 для одной     |
|          | точки интеграции.                                        |
|          | 2 - Для каждой точки интеграции девиаторный тензор       |
|          | напряжений исчезает, если D >= 1.                        |
|          | (Целое число)                                            |
+----------+-----------------------------------------------------------+
| I_moy    | Флаг модели отказа в 3D (твердые элементы).              |
|          | 0 - Использует следующие вычисления:                     |
|          | :math:`\bar{\eta} = \frac{\sigma_m}{\sigma_{VM}}`        |
|          | :math:`\bar{\xi} = \frac{27J_3}{2\sigma_{VM}^3}`         |
|          | 1 - Использует интегральные расчетные формулы.           |
|          | (Целое число)                                            |
+----------+-----------------------------------------------------------+
| fail_ID  | Идентификатор критерия отказа.                           |
|          | (Целое число, максимум 10 цифр)                          |
+----------+-----------------------------------------------------------+

Пример (Металл)
---------------

.. code-block::

    #RADIOSS STARTER
    #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
    /UNIT/1
    unit for mat
    #              MUNIT               LUNIT               TUNIT
                      Mg                  mm                   s
    #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
    #-  1. MATERIALS:
    #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
    /MAT/PLAS_JOHNS/1/1
    Metal
    #              RHO_I
                  2.8E-9                   0
    #                  E                  Nu
                   60000                  .3
    #                  a                   b                   n           EPS_p_max            SIG_max0
                     370                 190                  .4                   0                   0
    #                  c           EPS_DOT_0       ICC   Fsmooth               F_cut               Chard
                       0                   0         0         0                   0                   0
    #                  m              T_melt              rhoC_p                 T_r
                       0                   0                   0                   0
    /FAIL/WIERZBICKI/1/1
    #                 C1                  C2                  C3                  C4                   m
                    0.87                1.77                0.21                0.56                   2
    #                  n  Ifail_sh  Ifail_so      Imoy
                     0.2         0         1         0
    #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
    #enddata
    #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

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

Для :math:`\bar{\varepsilon}_f` формула используется при флаге отказа:

.. figure:: ./media/8e176e72_fail_wierzbicki.png
   :align: center

Где,

.. math::

    \bar{\varepsilon}_\text{max} = C_1 e^{-C_2 \bar{\eta}}
    \bar{\varepsilon}_\text{min} = C_3 e^{-C_4 \bar{\eta}}

Для *Твердых элементов*:

Где,

- :math:`\sigma_{m}` - гидростатическое напряжение
- :math:`\sigma_{VM}` - напряжение по Мизесу
- :math:`J_3` - третий инвариант девиаторного напряжения