============= /FAIL/HC_DSSE ============= Ключевое слово формата блока Модель вязкого разрушения на основе деформации: Хосфорда-Кулона с областью эквивалентности оболочки и твердого тела. Критерии разрушения оболочек на основе нелинейных деформаций с линейным накоплением повреждений. Деформация разрушения описывается функцией Хосфорда-Кулона (для твердых тел см. /FAIL/EMC). Работает только с законами упругопластического материала. > 28. Этот критерий неудачи был разработан Кеунхваном Паком (Массачусетский институт Technology MIT) и Дирк Мор (Швейцарский федеральный технологический институт ETH Zurich). 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 "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID", "/FAIL/HC_DSSE/mat_ID/unit_ID" "Ifail_sh", "P_thickfail", "P_thickfail", "I-флаг", "", "", "", "", "", "" If I-флаг = 0 ; Карта 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 "a", "a", "b", "b", "c", "c", "d", "d", ":math:`n_{f}`", ":math:`n_{f}`" If I-флаг = 1 ; Карта 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 "c2", "c2", "c3", "c3", "c4", "c4", "Inst_str", "Inst_str", ":math:`n_{f}`", ":math:`n_{f}`" Дополнительная линия .. 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 цифры)", "" "Ifail_sh", "Флаг отказа оболочки. = 1 (по умолчанию) Оболочка удаляется, если критерий поврежденности удовлетворяется за одно интегрирование. точка или слой. = 2 Оболочка удаляется, если критерий повреждения полностью удовлетворен или определен соотношение в P_thickfail, оболочка слои. Для однослойной оболочки тензор напряжений равен нулю. = 3 Ни один элемент не удаляется, ни тензор напряжений не устанавливается на ноль. (Целое число)", "" "P_thickfail", "Коэффициент интеграции по толщине точки, которые должны выйти из строя, прежде чем элемент будет удален. (только оболочки). По умолчанию = 1,0 (Реал)", "" "I-флаг", "Флаг типа ввода. = 0 (по умолчанию) Ввод параметров с помощью a, b, в, д, :math:`n_{f}` . = 1 Ввод теста отказа с использованием c1, c2, с3, с4 1. (Целое число)", "" "a", "Параметр модели отказа а.(Реал)", "" "b", "Параметр модели отказа б.(Реал)", "" "c", "Параметр модели отказа в.(Реал)", "" "d", "Параметр модели отказа д.(Реал)", "" "c2", "Деформация разрушения при чистом сдвиге (трехосность :math:`\eta=0` ).(Реал)", "" "c3", "Деформация разрушения при одноосном растяжении (трехосность :math:`\eta=\frac{1}{3}` ).(Реал)", "" "c4", "Деформация разрушения при плоскостном растяжении (трехосность = :math:`\eta=\frac{1}{\sqrt{3}}` ).(Реал)", "" "Inst_str", "Деформация нестабильности при плоской деформации напряжение для локализованного образования шейки. (Реальное)", "" ":math:`n_{f}`", "Параметр :math:`n_{f}` .По умолчанию = 0,1 (Реальное)", "" "fail_ID", "(Необязательно) Критерии отказа идентификатор.(Целое число, максимум 10 цифр)", "" Пример ------ В этом примере используется ввод параметра HC-DSSE (a, b, параметры c, d с I-флаг=0). :math:`I−Flag=0a=1.742b=0.7c=0.029d=1.6` ⇔ :math:`I−Flag=1c2=0.8(Pure−Shear)c3=0.7(Uniaxial−Tension)c4=0.58(Plane−Strain−Tension)` и один физический вход использования (вход напряжения при отказе с I-Флаг=1). HC-DSSE параметры а, b, c, d будут рассчитываться с использованием кривая, подобранная с помощью Radioss Starter. Определение отказа в этих двух примерах эквивалентно. .. image:: images/fail_hc_dsse_starter_r_hc_dsse_example1.png *(Рисунок 1.)* Пример (ввод параметра HC-DSSE) ------------------------------- .. code-block:: #RADIOSS STARTER #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /UNIT/1 unit for mat # MUNIT LUNIT TUNIT kg mm ms #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #- 1. MATERIALS: /MAT/PLAS_TAB/2/1 DP600 from SSAB Homepage # RHO_I 7.8E-6 0 # E Nu Eps_p_max Eps_t Eps_m 210 .3 0 0 0 # N_funct F_smooth C_hard F_cut Eps_f 1 0 0 0 0 # fct_IDp Fscale Fct_IDE EInf CE 0 0 0 0 0 # func_ID1 func_ID2 func_ID3 func_ID4 func_ID5 14 # Fscale_1 Fscale_2 Fscale_3 Fscale_4 Fscale_5 1 # Eps_dot_1 Eps_dot_2 Eps_dot_3 Eps_dot_4 Eps_dot_5 0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FAIL/HC_DSSE/2/1 # Ishell P_thickfail I_Flag 1 .5 0 # a b c d n_f 1.742 0.7 0.029 1.6 0.1 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/14 Mat_Curev Quasi-static DOCOL DP 600 (Material from SSAB Homepage 2010) # X Y 0 .306 .00112 .415 .00218 .445 .003 .461 .00404 .474 .00517 .489 .00613 .498 .0071 .505 .00806 .512 .00901 .522 .0102 .53 .0121 .543 .013 .55 .014 .555 .015 .561 .0159 .567 .0171 .572 .0181 .577 .0204 .592 .0303 .632 .0405 .663 .0502 .687 .06 .706 .0702 .722 .0807 .737 .09 .749 .0997 .758 .101 .759 .11 .768 .15000001 .805 .2 .84 .30000001 .9 .5 1 1 1.21 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #enddata /END #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| Пример (физический ввод) ------------------------ .. code-block:: #RADIOSS STARTER #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /UNIT/1 unit for mat # MUNIT LUNIT TUNIT kg mm ms #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #- 1. MATERIALS: /MAT/PLAS_TAB/2/1 DP600 from SSAB Homepage # RHO_I 7.8E-6 0 # E Nu Eps_p_max Eps_t Eps_m 210 .3 0 0 0 # N_funct F_smooth C_hard F_cut Eps_f 1 0 0 0 0 # fct_IDp Fscale Fct_IDE EInf CE 0 0 0 0 0 # func_ID1 func_ID2 func_ID3 func_ID4 func_ID5 14 # Fscale_1 Fscale_2 Fscale_3 Fscale_4 Fscale_5 1 # Eps_dot_1 Eps_dot_2 Eps_dot_3 Eps_dot_4 Eps_dot_5 0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FAIL/HC_DSSE/2/1 # Ishell P_thickfail I_Flag 1 .5 1 # c2 c3 c4 Inst_str n_f 0.8 0.7 0.58 0.1 0.1 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/14 Mat_Curev Quasi-static DOCOL DP 600 (Material from SSAB Homepage 2010) # X Y 0 .306 .00112 .415 .00218 .445 .003 .461 .00404 .474 .00517 .489 .00613 .498 .0071 .505 .00806 .512 .00901 .522 .0102 .53 .0121 .543 .013 .55 .014 .555 .015 .561 .0159 .567 .0171 .572 .0181 .577 .0204 .592 .0303 .632 .0405 .663 .0502 .687 .06 .706 .0702 .722 .0807 .737 .09 .749 .0997 .758 .101 .759 .11 .768 .15000001 .805 .2 .84 .30000001 .9 .5 1 1 1.21 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #enddata /END #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| Комментарии ----------- 1. Этот критерий отказа определен в пространство эквивалентной пластической деформации при разрушении в зависимости от трехосности напряжений (состояние напряжение) для плосконапряженного состояния. Это позволяет определить деформацию разрушения, которая варьируется в зависимости от различных условий нагрузки. Он состоит из двух кривых. Один из них – перелом локус (то есть ниже красной кривой), а другой — локализованный локус перешейка (синий кривая). .. image:: images/fail_hc_dsse_starter_r_hc_dsse_graph.png *(Рисунок 2.)* Первый принимает следующую форму: :math:`\bar{\epsilon}_{HC}^{pr}(\eta,\theta)=b(1+c)^{\frac{1}{n_{f}}}{[\frac{1}{2}((f_{1}−f_{2})^{a}+(f_{2}−f_{3})^{a}+(f_{1}−f_{3})^{a})]^{\frac{1}{a}}+c(2\eta+f_{1}+f_{3})}^{−\frac{1}{n_{f}}}` :math:`f_{1}` , :math:`f_{2}` и :math:`f_{3}` являются функциями угла Лоде :math:`\theta` : :math:`with{f_{1}(\theta)=\frac{2}{3}cos(\frac{\pi}{6}(1−\theta))f_{2}(\theta)=\frac{2}{3}cos(\frac{\pi}{6}(3+\theta))f_{3}(\theta)=−\frac{2}{3}cos(\frac{\pi}{6}(1+\theta))` :math:`\theta=1−\frac{2}{\pi}arccos(−\frac{27}{2}\eta(\eta^{2}−\frac{1}{3}))` :math:`dD_{HC}=\frac{d\bar{\epsilon}_{p}}{\bar{\epsilon}_{HC}^{pr}(\eta)}` Последнее определяется только между одноосными и условия экви-двухосного растяжения для ( :math:`\frac{1}{3}<\eta<\frac{2}{3}` ), при котором может возникнуть локализованное образование шейки. :math:`\bar{\epsilon}_{DSSE}^{pr}(\eta)=b{[\frac{1}{2}((g_{1}−g_{2})^{d}+g_{1}^{d}+g_{2}^{d})]^{\frac{1}{d}}}^{−\frac{1}{p}}` :math:`with{g_{1}(\eta)=\frac{2}{3}\eta+\sqrt{\frac{1}{3}−\frac{3}{4}\eta^{2}}g_{2}(\eta)=\frac{2}{3}\eta−\sqrt{\frac{1}{3}−\frac{3}{4}\eta^{2}}` :math:`dD_{DSSE}=\frac{d\bar{\epsilon}_{p}}{\bar{\epsilon}_{DSSE}^{pr}(\eta)}` - Универсальный показатель :math:`p` =0,01 используется. - Локализованное сужение прогнозируется нелокальным способом, когда все точки интеграции или слои имеют :math:`dD_{DSSE}>1` . - Параметры кривой разрушения a , b , c , и d , должен быть откалиброван на основе испытаний результаты. :math:`\eta` это трехосность. - Альтернативно, переключив I-Flag на 1, позволяет Radioss выполнять внутреннюю настройку параметров. - Разлома оболочки в зоне сжатия нет ( :math:`\eta<−\frac{1}{3}` ). - Ущерб рассчитывается как: :math:`D=\sum\frac{\Delta \bar{\epsilon}_{p}}{\bar{\epsilon}_{HC}^{pr}(\eta)}` - Контур зоны можно построить в H3D файлы с использованием /H3D/ELEM/HC_DSSE_Z : - Зона 1: Ниже кривых HC и DSSE - Зона 2: выше кривой DSSE и ниже кривой HC. - Зона 3: выше кривых HC и DSSE - Фактор отказа можно отобразить в файлах H3D, используя /H3D/ELEM/HC_DSSE_F. 1 Пак, Кынхван и Дирк Мор. " Комбинированная модель образования шейки и разрушения для прогнозирования пластического разрушения с конечной оболочкой элементы .. Инженерная механика разрушения 182 (2017): 32-51.