========== /FAIL/PUCK ========== Ключевое слово формата блока Описывает модель отказа шайбы. Формат ------ .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID", "/FAIL/PUCK/mat_ID/unit_ID" ":math:`\sigma_{1}^{t}`", ":math:`\sigma_{1}^{t}`", ":math:`\sigma_{2}^{t}`", ":math:`\sigma_{2}^{t}`", ":math:`\bar{\sigma}_{12}`", ":math:`\bar{\sigma}_{12}`", ":math:`\sigma_{1}^{c}`", ":math:`\sigma_{1}^{c}`", ":math:`\sigma_{2}^{c}`", ":math:`\sigma_{2}^{c}`" ":math:`p_{12}^{+}`", ":math:`p_{12}^{+}`", ":math:`p_{12}^{−}`", ":math:`p_{12}^{−}`", ":math:`p_{22}^{−}`", ":math:`p_{22}^{−}`", ":math:`\tau_{max}`", ":math:`\tau_{max}`", "Ifail_sh", "Ifail_so" "Fcut", "Fcut", "", "", "", "", "", "", "", "" Дополнительная линия .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "fail_ID", "", "", "", "", "", "", "", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "mat_ID", "Идентификатор материала.(Целое число, максимум 10 цифр)", "" "unit_ID", "Идентификатор объекта.(Целое число, максимум 10 цифр)", "" ":math:`\sigma_{1}^{t}`", "Продольная прочность на разрыв. По умолчанию = 1030 (Реальный)", ":math:`[Pa]`" ":math:`\sigma_{2}^{t}`", "Поперечная прочность на разрыв. По умолчанию = 1030 (Реал)", ":math:`[Pa]`" ":math:`\bar{\sigma}_{12}`", "Прочность на сдвиг. По умолчанию = 1030 (Реал)", ":math:`[Pa]`" ":math:`\sigma_{1}^{c}`", "Продольное сжатие сила.По умолчанию = 1030 (Реальный)", ":math:`[Pa]`" ":math:`\sigma_{2}^{c}`", "Поперечное сжатие сила.По умолчанию = 1030 (Реальный)", ":math:`[Pa]`" ":math:`p_{12}^{+}`", "Коэффициент конверта отказа 12 (+).По умолчанию = 0 (Реальный)", "" ":math:`p_{12}^{−}`", "Коэффициент конверта отказа 12 (-).По умолчанию = 0 (Реальное)", "" ":math:`p_{22}^{−}`", "Коэффициент конверта отказа 22 (-).По умолчанию = 0 (Реальное)", "" ":math:`\tau_{max}`", "Динамическое время релаксации. 5По умолчанию = 1030 (Реальное)", ":math:`[s]`" "Ifail_sh", "Флаг модели отказа оболочки. = 1 (по умолчанию) Оболочка удаляется, если повреждение достигнуто за один слой. = 2 Оболочка удаляется, если повреждения достигнуты для всех слоев оболочки. (Целое число)", "" "Ifail_so", "Флаг модели твердого отказа. = 1 (по умолчанию) Тело удаляется, если повреждение достигнуто для одной точки интеграции твердый. (Целое число)", "" "Fcut", "Тензорная фильтрация напряжений частота.По умолчанию = 0,0 (Реальный)", ":math:`[\frac{1}{s}]`" "fail_ID", "Критерии отказа идентификатор. 4(Целое число, максимум 10 цифр)", "" Пример (составной) ------------------ Сила ( :math:`\sigma_{1}^{t}` , :math:`\sigma_{2}^{t}` , :math:`\sigma_{1}^{c}` , :math:`\sigma_{2}^{c}` , :math:`\bar{\sigma}_{12}` ) взяты из следующих тестов. m1 есть направление волокон. .. image:: images/fail_puck_starter_r_fail_puck_example.png :alt: fail_puck_example .. code-block:: #RADIOSS STARTER #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /UNIT/1 unit for mat and failure # MUNIT LUNIT TUNIT g mm ms #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #- 1. MATERIALS: #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /MAT/COMPSH/1/1 composite example # RHO_I .0015 # E11 E22 NU12 Iform E33 114000 9650 .025 0 0 # G12 G23 G31 EPS_f1 EPS_f2 6000 6000 6000 0 0 # EPS_t1 EPS_m1 EPS_t2 EPS_m2 dmax 0 0 0 0 0 # Wpmax Wpref Ioff ratio 0 0 4 0 # b n fmax 0 0 0 # sig_1yt sig_2yt sig_1yc sig_2yc alpha 1E30 1E30 1E30 1E30 0 # sig_12yc sig_12yt c_12 Eps_rate_0 ICC 1E30 1E30 0 0 0 # GAMMA_ini GAMMA_max d3max 0 0 0 # Fsmooth Fcut 0 0 /FAIL/PUCK/1/1 # Sigma1_T Sigma2_T Sigma_12 Sigma1_C Sigma2_C 1720 55.2 103 765 503 # P+12 P-12 P-22 Tau_max Ifail_sh Ifail_so 0 0 0 .005 1 0 # Fcut 0.1 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #enddata #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| Комментарии ----------- 1. Эта модель отказа доступна для Shell. и Твердый. 2. Критерии режима отказа записаны как: - Недостаточность фракции клетчатки: Режим разрушения волокна при растяжении: :math:`\sigma_{11}>0` :math:`e_{f}=\frac{\sigma_{11}}{\sigma_{1}^{t}}` Режим разрушения сжимающего волокна: :math:`\sigma_{11}<0` :math:`e_{f}=\frac{|\sigma_{11}|}{\sigma_{1}^{c}}` - Межволоконный отказ: Режим А, если :math:`\sigma_{22}>0` : .. image:: images/fail_puck_starter_r_fail_puck_modea.png :alt: fail_puck_modeA *(Рисунок 1.)* :math:`e_{f}=\frac{1}{\bar{\sigma}_{12}}[\sqrt{(\frac{\bar{\sigma}_{12}}{\sigma_{2}^{t}}−p_{12}^{+})^{2}\sigma_{22}^{2}+\sigma_{12}^{2}}+p_{12}^{+}\sigma_{22}]` Режим C, если :math:`\sigma_{22}<0` : .. image:: images/fail_puck_starter_r_fail_puck_modec.png :alt: fail_puck_modeC *(Рисунок 2.)* :math:`e_{f}=[(\frac{\sigma_{12}}{2(1+p_{22}^{−})\bar{\sigma}_{12}})^{2}+(\frac{\sigma_{22}}{\sigma_{2}^{c}})^{2}](\frac{\sigma_{2}^{c}}{−\sigma_{22}})` Режим Б .. image:: images/fail_puck_starter_r_fail_puck_modeb.png :alt: fail_puck_modeB *(Рисунок 3.)* :math:`e_{f}=\frac{1}{\bar{\sigma}_{12}}(\sqrt{\sigma_{12}^{2}+(p_{12}^{−}\sigma_{_{22}})^{2}}+p_{12}^{−}\sigma_{_{22}})` Если параметр повреждения :math:`e_{f}\ge1.0` , напряжения уменьшаются с помощью экспоненциальной зависимости функция, позволяющая избежать числовых нестабильностей. Используется техника релаксации путем уменьшения стресс постепенно: :math:`\sigma(t)=f(t)⋅\sigma_{d}(t_{r})` С, :math:`f(t)=exp(−\frac{t−t_{r}}{\tau_{max}})` и :math:`t\get_{r}` Где, :math:`t` Время :math:`t_{r}` Время начала релаксации, когда принимаются критерии повреждения :math:`\tau_{max}` Время динамического расслабления :math:`\sigma_{d}(t_{r})` Стресс в начале повреждения 3. Величина ущерба D равна :math:`0\leD\le1` . Статус перелома: - Бесплатно, если :math:`0\leD<1` - Неудача, если :math:`D=1` С :math:`D=Max(e_{f}(tensile),e_{f}(compression),e_{f}(ModeA),e_{f}(ModeB),e_{f}(ModeC))` . Это значение урона отображается с помощью /ANIM/BRICK/DAMA, /ANIM/SHELL/DAMA, H3D/SOLID/DAMA или /H3D/SHELL/DAMA. 4. fail_ID используется с /STATE/BRICK/FAIL и /INIBRI/FAIL. Нет значения по умолчанию ценность. Если строка пуста, для переменных модели отказа в /INIBRI/FAIL (записано в файле .sta с /STATE/BRICK/FAIL вариант). 5. После достижения критерия отказа тот :math:`\tau_{max}` Значение определяет период времени, в течение которого напряжение в вышедший из строя элемент постепенно сводится к нулю. Когда напряжение достигает 1% от значения напряжения при начало сбоя, элемент удаляется. Это необходимо, чтобы избежать нестабильности. происходит из-за внезапного удаления элемента и сбоя» цепная реакция » в соседние элементы. Даже если критерий отказа достигнут, значение по умолчанию :math:`\tau_{max}=1.0E30` не приводит к удалению элемента. Поэтому рекомендуется определить :math:`\tau_{max}` В 10 раз больше, чем шаг моделирования. 6. Чтобы избежать « цепная реакция » при удалении элементов вы также можете определить частоту фильтрации тензора напряжений F резать . Таким образом, тензор напряжений, используемый для расчета критерия PUCK, сначала подвергается фильтрации. по формуле ниже: :math:`\sigma_{n+1}^{filt}=\alpha\sigma_{n+1}+1−\alpha\sigma_{n}^{filt}` С :math:`\alpha=\frac{2\pi⋅F_{cut}⋅\Delta t}{2\pi⋅F_{cut}⋅\Delta t+1}` Где, :math:`\Delta t` текущий временной шаг. Если нет частоты фильтрации задано (Fcut = 0,0), эффект фильтрации деактивирован. 7. Различные виды отказов можно отобразить с помощью опции /H3D/ELEM/FAILURE/ID = FAILID с опцией MODE (= I or ВСЕ ): - РЕЖИМ1 : Разрушение волокна при растяжении. :math:`e_{f}(tensile)` - РЕЖИМ2 : Сбой компрессионного волокна :math:`e_{f}(compression)` - РЕЖИМ3 : Межволоконный отказ A :math:`e_{f}ModeA` - РЕЖИМ4 : Межволоконный отказ B :math:`e_{f}ModeB` - РЕЖИМ5 : Межволоконный отказ C :math:`e_{f}ModeC` . Глобальный индекс отказов (который соответствует максимуму между всеми режимами) может быть построено с использованием /H3D/ELEM/FAILURE/ID=FAILID без опция MODE.