/FAIL/LAD_DAMA

Ключевое слово формата блока Описывает модель разрушения Ладевезе при расслоении (межламинарном переломе).

Эта модель разрушения доступна для ортотропных твердых тел и толстых оболочек. Его также можно было бы использовать

с Plyxfem в свойстве оболочки /PROP/TYPE17 как интерактивная модель разрушения материала. Эта модель отказа совместима с /MAT/LAW12 (3D_COMP), /MAT/LAW14 (COMPSO) и /MAT/LAW25 (COMPSH) и /MAT/LAW1 (ELAST) (только при использовании с Поликсфем).

Формат

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

/FAIL/LAD_DAMA/mat_ID/unit_ID

K1

K1

K2

K2

K3

K3

\(\gamma_{1}\)

\(\gamma_{1}\)

\(\gamma_{2}\)

\(\gamma_{2}\)

\(Y_{0}\)

\(Y_{0}\)

\(Y_{c}\)

\(Y_{c}\)

k

k

a

a

\(\tau_{max}\)

\(\tau_{max}\)

Ifail_sh

Ifail_so

Дополнительная линия .. 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 цифр)

K1

Межламинарная жесткость в направление 1. По умолчанию = 1030 (реальное)

\([\frac{Pa}{m}]\)

K2

Межламинарная жесткость в направление 2. По умолчанию = 1030 (реальное)

\([\frac{Pa}{m}]\)

K3

Межламинарная жесткость в направление 3. По умолчанию = 1030 (реальное)

\([\frac{Pa}{m}]\)

\(\gamma_{1}\)

Коэффициент связи между режим расслоения I и режим II. По умолчанию = 0 (реальный)

\(\gamma_{2}\)

Коэффициент связи между режим расслоения I и режим III. По умолчанию = 0 (реальный)

\(Y_{0}\)

Нанести энергетический урон за начало расслоения.По умолчанию = 1030 (Реальное)

\(Y_{c}\)

Критический энергетический урон параметр для полного расслоения. По умолчанию = \(2\sqrt{Y_{0}}\) (Реал)

k

Скорость распространения трещины постоянная времени. По умолчанию = 0 (реальная)

\([\frac{m}{s}]\)

a

Скорость распространения трещины множитель.По умолчанию = 1030 (Реальный)

\(\tau_{max}\)

Динамическое время релаксации. 3По умолчанию = 1030 (Реал)

\([s]\)

Ifail_sh

Флаг отказа оболочки. = 1 (по умолчанию) Снаряд удаляется, если достигнут критерий повреждения за один слой. = 2 Оболочка удаляется, если критерий повреждения достигнут для всех слои оболочки. (Целое число)

Ifail_so

Твердый флаг отказа. = 1 (по умолчанию) Тело удаляется, если критерий повреждения достигнут за один точка интеграции. = 3 Внеплоскостное напряжение устанавливается равным нулю, если достигнуто повреждение. за одну точку интегрирования твердого тела ( \(\sigma_{33}=\sigma_{23}=\sigma_{13}=0\) ). (Целое число)

fail_ID

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

Пример (составной)

#RADIOSS STARTER

/UNIT/1

unit for mat and failure

#              MUNIT               LUNIT               TUNIT

                   g                  mm                  ms

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

/MAT/COMPSH/1/1

composite example

#              RHO_I

             .001506

#                E11                 E22                NU12     Iform                           E33

              144000               10000                 .25         0                         20000

#                G12                 G23                 G31              EPS_f1              EPS_f2

                4200                4200                4200                   0                   0

#             EPS_t1              EPS_m1              EPS_t2              EPS_m2                dmax

                   0                   0                   0                   0                   0

#              Wpmax               Wpref      Ioff                         ratio

             1000000                   0         0                             0

#                  b                   n                fmax

                   0                   0             1000000

#            sig_1yt             sig_2yt             sig_1yc             sig_2yc               alpha

               10100               10100               10100               10100                   0

#           sig_12yc            sig_12yt                c_12          Eps_rate_0       ICC

               10068               10068                   0                   0         0

#          GAMMA_ini           GAMMA_max               d3max

                   0                   0                   0

#  Fsmooth                Fcut

         0                   0

/FAIL/LAD_DAMA/1/1

#                 K1                  K2                  K3              Gamma1              Gamma2

                2000                2000                2000               1E-20               1E-20

#                 Y0                  YC                   K                   A             Tau_max

                  40                 160              100000                   1                 .01

# Ifail_sh  Ifail_so

         1         3

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

#enddata

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

Комментарии

  1. Ущерб от аварии в Ладевезе

    модель для расслоения:

\(Y_{d_{3}}=\frac{\partialE_{D}}{\partiald_{3}}|_{\sigma=cst}=\frac{1}{2}\frac{〈\sigma_{33}〉^{2}}{K_{3}(1−d_{3})^{2}}ModeI\) \(Y_{d_{1}}=\frac{\partialE_{D}}{\partiald_{1}}|_{\sigma=cst}=\frac{1}{2}\frac{\sigma_{31}^{2}}{K_{1}(1−d_{1})^{2}}ModeII\) \(Y_{d_{2}}=\frac{\partialE_{D}}{\partiald_{2}}|_{\sigma=cst}=\frac{1}{2}\frac{\sigma_{32}^{2}}{K_{2}(1−d_{2})^{2}}ModeIII\) .. image:: images/fail_lad_dama_starter_r_fail_lad-dama_delam1a.png

alt

fail_lad-dama_delam1A

(Рисунок 1.)

2 Для элемента Quad 2D только режим II.

и режим III доступны.

Где,

\(d_{i}\) – параметры внутреннего повреждения, связанные с его режимом разрушения.

Закон развития повреждений контролируется

эквивалентная скорость выделения энергии повреждения.

\(Y=Y_{d_{3}}+\gamma_{1}Y_{d_{1}}+\gamma_{2}Y_{d_{2}}\)

с \(Y_{d_{i}}|_{t}=supY_{d_{i}}|_{\tau\let}\)

Эволюция параметров повреждений

сильно связан с коэффициентом связи \(\gamma_{1}\) и \(\gamma_{2}\)

. Эти два материальных параметра происходят из
испытания на расслаивание.

Для настоящей модели отказа учтите, что

\(d_{1}=d_{2}=d_{3}=d\)

.

Значение урона

\(d\) увеличивается при определенных скорость:

\(\dot{d}=\frac{k}{a}[1−exp(−a〈w(Y)−d〉)]\) если

\(d<1\)

.

В противном случае

\(d=1\)

Пока, \(a\)

является мерой пластичности разрушения, Чем ниже значение, тем более пластичным является разрушение.

\(\frac{a}{k}\)

– минимальная продолжительность отказа. продолжительность энергии между \(Y_{0}\) и \(Y_{c}\) должно быть как минимум равно \(\frac{a}{k}\)

.

\(〈x〉={xifx>00ifx<0\) Функция

\(w(Y)\) рассчитывается как:

\(w(Y)=\frac{〈\sqrt{Y}−\sqrt{Y_{0}}〉}{\sqrt{Y_{c}}−\sqrt{Y_{0}}}\) Если параметр повреждения

\(d\le1.0\)

, напряжения

\(\sigma_{33}\)

,

\(\sigma_{13}\) и \(\sigma_{23}\) уменьшаются в соответствии со следующим функция:

Применяется техника релаксации путем постепенного уменьшения

стресс:

\(\sigma(t)=f(t)⋅\sigma_{d}(t_{r})\) С, \(f(t)=exp(−\frac{t−t_{r}}{\tau_{max}})andt\get_{r}\) Где, \(t\) Время \(t_{r}\) Время начала релаксации, когда критерии повреждения

предполагалось

\(\tau_{max}\) Время динамического расслабления \(\sigma_{d}(t_{r})\) Стресс в начале повреждения

  1. fail_ID используется с /STATE/BRICK/FAIL и /INIBRI/FAIL. Нет

значение по умолчанию. Если строка пуста, значение для модели отказа выведено не будет. переменные в /INIBRI/FAIL (написаны на .sta файл с /STATE/BRICK/FAIL вариант).

  1. После критерия отказа

    достигается,

\(\tau_{max}\) Значение определяет период времени, в течение которого

напряжение в вышедшем из строя элементе постепенно сводится к нулю. Когда стресс

достигает 1% от значения напряжения в начале отказа, элемент удаляется. Это

необходимо, чтобы избежать нестабильности, возникающей из-за внезапного удаления элемента и

неудача»

цепная реакция »в соседних элементах.Даже если

критерий отказа достигнут, значение по умолчанию

\(\tau_{max}=1.0E30\) не приводит к удалению элемента. Следовательно, это

рекомендуется определить

\(\tau_{max}\) В 10 раз больше, чем время моделирования

шаг.

1 О. Алликс, П. Ладевез,

«Моделирование межламинарного интерфейса для прогнозирования расслоения», Композит

структура 22 (1992) 235-242

2 Л. Горнет, «Прогнозирование повреждения композитных материалов методом конечных элементов».

Структуры”