/FAIL/HC_DSSE

Ключевое слово формата блока Модель вязкого разрушения на основе деформации: Хосфорда-Кулона

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

Деформация разрушения описывается функцией Хосфорда-Кулона (для твердых тел см. /FAIL/EMC). Работает только с законами упругопластического материала.

> 28. Этот критерий неудачи был разработан Кеунхваном Паком (Массачусетский институт Technology MIT) и Дирк Мор (Швейцарский федеральный технологический институт ETH Zurich). 1

Формат

/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 — Накопление урона

параметры

a

a

b

b

c

c

d

d

\(n_{f}\)

\(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", "", "", "", "", "", "", "", "", ""

Определение

Поле

Содержание

Пример единицы СИ

mat_ID

Идентификатор материала.(Целое число, максимум 10 цифр)

unit_ID

Идентификатор объекта.(Целое число, максимум 10 цифры)

Ifail_sh

Флаг отказа оболочки. = 1 (по умолчанию) Оболочка удаляется, если критерий поврежденности удовлетворяется за одно интегрирование. точка или слой. = 2 Оболочка удаляется, если критерий повреждения полностью удовлетворен или определен соотношение в P_thickfail, оболочка слои. Для однослойной оболочки тензор напряжений равен нулю. = 3 Ни один элемент не удаляется, ни тензор напряжений не устанавливается на ноль. (Целое число)

P_thickfail

Коэффициент интеграции по толщине точки, которые должны выйти из строя, прежде чем элемент будет удален. (только оболочки). По умолчанию = 1,0 (Реал)

I-флаг

Флаг типа ввода. = 0 (по умолчанию) Ввод параметров с помощью a, b, в, д, \(n_{f}\) . = 1 Ввод теста отказа с использованием c1, c2, с3, с4 1. (Целое число)

a

Параметр модели отказа а.(Реал)

b

Параметр модели отказа б.(Реал)

c

Параметр модели отказа в.(Реал)

d

Параметр модели отказа д.(Реал)

c2

Деформация разрушения при чистом сдвиге (трехосность \(\eta=0\) ).(Реал)

c3

Деформация разрушения при одноосном растяжении (трехосность \(\eta=\frac{1}{3}\) ).(Реал)

c4

Деформация разрушения при плоскостном растяжении (трехосность = \(\eta=\frac{1}{\sqrt{3}}\) ).(Реал)

Inst_str

Деформация нестабильности при плоской деформации напряжение для локализованного образования шейки. (Реальное)

\(n_{f}\)

Параметр \(n_{f}\) .По умолчанию = 0,1 (Реальное)

fail_ID

(Необязательно) Критерии отказа идентификатор.(Целое число, максимум 10 цифр)

Пример

В этом примере используется ввод параметра HC-DSSE (a, b,

параметры c, d с I-флаг=0).

\(I−Flag=0a=1.742b=0.7c=0.029d=1.6\)

\(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.

Определение отказа в этих двух примерах эквивалентно.

Пример (ввод параметра HC-DSSE)

#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----|

Пример (физический ввод)

#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. Этот критерий отказа определен в

    пространство эквивалентной пластической деформации при разрушении в зависимости от трехосности напряжений (состояние

    напряжение) для плосконапряженного состояния. Это позволяет определить деформацию разрушения, которая

    варьируется в зависимости от различных условий нагрузки. Он состоит из двух кривых. Один из них – перелом

    локус (то есть ниже красной кривой), а другой — локализованный локус перешейка (синий

    кривая).

images/fail_hc_dsse_starter_r_hc_dsse_graph.png*(Рисунок2.)*

Первый принимает следующую форму: \(\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}}}\) \(f_{1}\)

,

\(f_{2}\) и \(f_{3}\) являются функциями угла Лоде \(\theta\)

:

\(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))\) \(\theta=1−\frac{2}{\pi}arccos(−\frac{27}{2}\eta(\eta^{2}−\frac{1}{3}))\) \(dD_{HC}=\frac{d\bar{\epsilon}_{p}}{\bar{\epsilon}_{HC}^{pr}(\eta)}\) Последнее определяется только между одноосными и

условия экви-двухосного растяжения для ( \(\frac{1}{3}<\eta<\frac{2}{3}\)

), при котором может возникнуть локализованное образование шейки.

\(\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}}\) \(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}}\) \(dD_{DSSE}=\frac{d\bar{\epsilon}_{p}}{\bar{\epsilon}_{DSSE}^{pr}(\eta)}\) - Универсальный показатель

\(p\) =0,01 используется.

  • Локализованное сужение прогнозируется нелокальным способом, когда все точки интеграции

    или слои имеют

    \(dD_{DSSE}>1\) .

  • Параметры кривой разрушения a , b , c , и d , должен быть откалиброван на основе испытаний

    результаты.

    \(\eta\)

это трехосность.
  • Альтернативно, переключив I-Flag на 1,

позволяет Radioss выполнять внутреннюю настройку параметров.
  • Разлома оболочки в зоне сжатия нет ( \(\eta<−\frac{1}{3}\) ).

  • Ущерб рассчитывается как: \(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 Пак, Кынхван и Дирк Мор.

Комбинированная модель образования шейки и разрушения для прогнозирования пластического разрушения с конечной оболочкой

элементы