/MAT/LAW116

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

неудача.

Этот материал совместим только с твердыми шестигранными элементами ( /BRICK ) и когезионное твердое свойство ( /PROP/TYPE43 (CONNECT) ). .. note:

Не совместим ни с одним
 модель неудачи. Все повреждения и отказы определяются внутри материала.
 напрямую.

Формат

/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID	/MAT/LAW116/mat_ID/unit_ID
mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title	mat_title
\(\rho_{i}\)	\(\rho_{i}\)
\(E_{I}\)	\(E_{I}\)	\(E_{II}\)	\(E_{II}\)	Толстый	Толстый	Имасс	Идель	Икрит
\(GC_{I_ini}\)	\(GC_{I_ini}\)	\(GC_{I_inf}\)	\(GC_{I_inf}\)	\(\dot{\epsilon}_{G_{I}}\)	\(\dot{\epsilon}_{G_{I}}\)	\(fG_{I}\)	\(fG_{I}\)
\(GC_{II_ini}\)	\(GC_{II_ini}\)	\(GC_{II_inf}\)	\(GC_{II_inf}\)	\(\dot{\epsilon}_{G_{II}}\)	\(\dot{\epsilon}_{G_{II}}\)	\(fG_{II}\)	\(fG_{II}\)
\(\sigma_{A_I}\)	\(\sigma_{A_I}\)	\(\sigma_{B_I}\)	\(\sigma_{B_I}\)	\(\dot{\epsilon}_{I}\)	\(\dot{\epsilon}_{I}\)	Iorder_I	Ifail_I
\(\sigma_{A_II}\)	\(\sigma_{A_II}\)	\(\sigma_{B_II}\)	\(\sigma_{B_II}\)	\(\dot{\epsilon}_{II}\)	\(\dot{\epsilon}_{II}\)	Iorder_II	Ifail_II

Определение

Поле	Содержание	Пример единицы СИ
mat_ID	Идентификатор материала.(Целое число, максимум 10 цифр)
unit_ID	(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)
mat_title	Название материала.(Персонаж, максимум 100 символов)
\(\rho_{i}\)	Начальный плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
\(E_{I}\)	Модуль Юнга (жесткости) в нормальное направление на единицу длины.(Реальное)	\([\frac{Pa}{m}]\)
\(E_{II}\)	Модуль сдвига (жесткости) в направление касательной на единицу длины. По умолчанию = \(E_{II}=E_{I}\) (Реал)	\([\frac{Pa}{m}]\)
Толстый	Ссылка связная Толщина.(Реальная)	\([m]\)
Имасс	Флаг расчета массы. = 1 (по умолчанию) Масса элемента рассчитывается с использованием плотности и среднего значения. площадь. = 2 Масса элемента рассчитывается с использованием плотности и объем. (Целое число)
Идель	Флаг отказа с указанием номера точек интеграции для удаления элемента (между 1 и 4).По умолчанию = 1 (Целое число)
Икрит	Флаг инициации выхода и повреждения. = 1 (по умолчанию) На основе квадратичного номинального напряжения. = 2 На основе максимального номинального напряжения. (Целое число)
\(GC_{I_ini}\)	Начальное критическое выделение энергии скорость для режима I (нормальное направление).(Реальная)	\([J]\)
\(GC_{I_inf}\)	Верхняя граница критической энергии скорость выпуска. Указывает зависимость скорости деформации \(GC_{I}\) .По умолчанию = 0,0 (Реал)	\([J]\)
\(\dot{\epsilon}_{G_{I}}\)	Эталонная (нижняя граница) скорость деформации для зависимости скорости деформации GC. Должно быть определено, если \(GC_{I_inf}>0\) .(Реал)	\([Hz]\)
\(fG_{I}\)	Фактор формы для выделения энергии скорость до отказа в режиме I.(Реальная)
\(GC_{II_ini}\)	Начальное критическое выделение энергии скорость для режима II (сдвиг).(Реальная)	\([J]\)
\(GC_{II_inf}\)	Верхняя граница критической энергии скорость выпуска. Указывает зависимость скорости деформации \(GC_{II}\) .По умолчанию = 0,0 (Реал)	\([J]\)
\(\dot{\epsilon}_{G_{II}}\)	Эталонная (нижняя граница) скорость деформации для зависимости скорости деформации GC. Должно быть определено, если \(GC_{II_inf}>0\) .(Реал)	\([Hz]\)
\(fG_{II}\)	Фактор формы для выделения энергии скорость до отказа в режиме II.(Реальная)
\(\sigma_{A_I}\)	Статический предел текучести в режиме Я.(Реал)	\([Pa]\)
\(\sigma_{B_I}\)	Предел текучести, зависящий от скорости деформации срок в режиме I.(Реальный)	\([Pa]\)
\(\dot{\epsilon}_{I}\)	Эталонная (нижняя граница) скорость деформации значение зависимости скорости предела текучести в режиме I. Должно быть определенный \(\sigma_{B_I}>0\) .(Реал)	\([Hz]\)
Iorder_I	Порядок зависимости предела текучести от скорость деформации в режиме I. = 1 (по умолчанию) Линейно-логарифмическая зависимость скорости деформации. = 2 Квадратичная логарифмическая зависимость деформации ставка. (Целое число)
Ifail_I	Критерии отказа, определенные \(fG_{I}\) : = 1 (по умолчанию) Коэффициент энергии разрушения. = 2 Соотношение смещений переломов. (Целое число)
\(\sigma_{A_II}\)	Статический предел текучести в режиме II.(Реал)	\([Pa]\)
\(\sigma_{B_II}\)	Предел текучести, зависящий от скорости деформации срок в режиме II.(Реальный)	\([Pa]\)
\(\dot{\epsilon}_{II}\)	Эталонная (нижняя граница) скорость деформации значение зависимости предела текучести от скорости текучести в режиме II. Должно быть определено, если \(\sigma_{B_II}>0\) .(Реал)	\([Hz]\)
Iorder_II	Порядок зависимости предела текучести скорость деформации в режиме II. = 1 (по умолчанию) Линейно-логарифмическая зависимость скорости деформации. = 2 Квадратичная логарифмическая зависимость деформации ставка. (Целое число)
Ifail_II	Критерии отказа, определенные \(fG_{II}\) : = 1 (по умолчанию) Коэффициент энергии разрушения. = 2 Соотношение смещений переломов. (Целое число)

Пример

#RADIOSS STARTER

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

/UNIT/1

unit for mat

                  Mg                  mm                   s

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

/MAT/LAW116/3/1

MAT_COHESIVE_MIXED_MODE_ELASTOPLASTIC_RATE

#              RHO_I

              1.2E-9

#                 E1                  E2               Thick     Imass      Idel    Icrit

                3000                1000               0.200         2         1        0

#            GC1_INI             GC1_INF              SRATG1                 FG1

               2.000               3.000               1.500                 0.7

#            GC2_INI             GC2_INF              SRATG2                 FG2

                9.00                   0                   0                 0.4

#              SIGA1               SIGB1              SRATE1   Iorder1    Ifail1

               33.00               1.500          2.50000E-5         1         2

#              SIGA2               SIGB2              SRATE2   Iorder2    Ifail2

               26.00               1.300          1.00000E-5         1         2

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

#enddata

Комментарии

1. Упругая жесткость – это

   определяется с помощью:

Где, GP Пластическая энергия в условиях постоянного стресса GC Общая энергия \(i\) ={I,II} Режим I (нормальный) и режим II (сдвиг) Форма закона разделения тяги определяется следующим образом: - Критерии отказа

определяется соотношением энергии разрушения

(

I fail_i =1) \(0\lefG_{i}=\frac{GC_{i}(\dot{\epsilon}_{eq})}{GC_{i}(\dot{\epsilon}_{eq})}<1−\frac{\sigma(\dot{\epsilon}_{eq})^{2}}{2GC_{i}(\dot{\epsilon}_{eq})E_{i}}<1\)

• Критерии отказа

  определяется соотношением смещений переломов

  (

  I fail_i =2) \(0\lefG_{i}=\frac{\delta_{i2}−\delta_{i1}}{\delta_{if}−\delta_{i1}}<1\)

2. Предел текучести определяется

   как:

• Когда I order_i =1: \(\sigma(\dot{\epsilon}_{eq})=\sigma_{A_i}+\sigma_{B_i}.[max(0,ln(\frac{\dot{\epsilon}_{eq}}{\dot{\epsilon}_{i}}))]\)

• Когда I order_i =2: \(\sigma(\dot{\epsilon}_{eq})=\sigma_{A_i}+\sigma_{B_i}.[max(0,ln(\frac{\dot{\epsilon}_{eq}}{\dot{\epsilon}_{i}}))]^{2}\) Где,

\(i\)

={I,II}, режим I и режим

3. Эквивалентная скорость деформации

   определяется с помощью:

\(\dot{\epsilon}_{eq}=\frac{\sqrt{\dot{\Delta }^{2}_{I}+\dot{\Delta }^{2}_{II}}}{Thick}\) Где, \(\dot{\Delta }_{I}\) Нормальная скорость. \(\dot{\Delta }_{II}\) Скорость сдвига.

4. Зависимый от скорости перелом

   энергии определяются с помощью:

\(GC_{i}(\dot{\epsilon}_{eq})=GC_{i_ini}+(GC_{i_inf}−GC_{i_ini}).exp(−\frac{\dot{\epsilon}_{G_{i}}}{\dot{\epsilon}_{eq}})\) Где,

\(i\)

={I,II}, режим I и режим II. 5. Закон о пределе текучести и повреждении

схема:

6. Для урожайности и урона на основе

   по квадратичному номинальному напряжению

   (

Икрит = 1 ): - Начало урожайности в смешанном режиме

перемещение это:

\(\delta_{m1}=\delta_{I1}\delta_{II1}.\sqrt{\frac{1+\beta^{2}}{\delta_{II1}^{2}+(\beta.\delta_{I1})^{2}}}\) Где, \(\delta_{i1}=\frac{\sigma_{i}}{E_{i}}\) \(i\) ={I,II}, режим I и

режим II.

\(\beta=\frac{\Delta _{II}}{\Delta _{I}}\)

• Смешанный урон

  инициация – это:

  \(\delta_{m2}=\delta_{I2}\delta_{II2}.\sqrt{\frac{1+\beta^{2}}{\delta_{II2}^{2}+(\beta.\delta_{I2})^{2}}}\) Где, \(\delta_{i2}=\delta_{i1}+\frac{fG_{i}.GC_{i}}{\sigma_{i}}\) \(i\) ={I,II}, режим I и

  режим II.

7. Для урожайности и урона на основе

   по квадратичному номинальному напряжению

   (

Икрит = 2 ). - Начало урожайности в смешанном режиме

перемещение это:

Если

\(\beta\le\frac{\delta_{II1}}{\delta_{I1}}\)

:

\(\delta_{m1}=\delta_{I1}.\sqrt{1+\beta^{2}}\) Если

\(\beta>\frac{\delta_{II1}}{\delta_{I1}}\)

:

\(\delta_{m1}=\frac{\delta_{II1}}{\beta}.\sqrt{1+\beta^{2}}\) Где, \(\beta=\frac{\Delta _{II}}{\Delta _{I}}\) \(\Delta _{I}\) Смещение соответствует режиму I (нормальному). \(\Delta _{II}\) Смещение — режим II (сдвиг).

• Смешанный урон

  инициация – это:

  Если

\(\beta\le\frac{\delta_{II2}}{\delta_{I2}}\)

:

\(\delta_{m2}=\delta_{I2}.\sqrt{1+\beta^{2}}\) Если

\(\beta>\frac{\delta_{II2}}{\delta_{I2}}\)

:

\(\delta_{m2}=\frac{\delta_{II2}}{\beta}.\sqrt{1+\beta^{2}}\)

8. Окончательный урон в смешанном режиме

   есть

   (

Икрит = 1 , 2 ): \(\delta_{mf}=\frac{\delta_{m1}.(\delta_{m1}−\delta_{m2})E_{I}GC_{II}cos^{2}\gamma+GC_{I}.(2GC_{II}+\delta_{m1}.(\delta_{m1}−\delta_{m2})E_{II}sin^{2}\gamma)}{\delta_{m1}(E_{I}GC_{II}cos^{2}\gamma+E_{II}GC_{I}sin^{2}\gamma)}\) Где, \(\gamma=arccos(\frac{\Delta _{I}}{\Delta _{m}})\) \(\Delta _{m}\) Перемещение – смешанный режим.

9. Пластическая деформация – это

   определяется как:

• Режим I: \(\Delta p_{I}=max(\Delta p_{I}(t−1),\Delta p_{I}−\delta_{m1}cos\gamma,0)\) Где,

\((t−1)\) значение предыдущего времени шаг.

• Режим II: Если

\(\sqrt{(\Delta _{II−1}−\Delta p_{II−1}(t−1))^{2}+(\Delta _{II−2}−\Delta p_{II−2}(t−1))^{2}}>\delta_{m1}\)

Пластическая деформация – это

вычисляется для каждого направления 1 и 2 сдвига самолет.

\(\Delta p_{II−1}=\Delta p_{II−1}(t−1)+\Delta _{II−1}−\Delta _{II−1}(t−1)\) \(\Delta p_{II−2}=\Delta p_{II−2}(t−1)+\Delta _{II−2}−\Delta _{II−2}(t−1)\)

10. Значение стресса снижается

    линейно от начала повреждения до окончательного повреждения (

\(\Delta _{m}>\delta_{m2}\) ). \(D=max(\frac{\Delta _{m}−\delta_{m2}}{\delta_{mf}−\delta_{m2}},D(t−1),0)\) Снижение напряжения рассчитывается в нормальном режиме.

направление как:

Если

\(\Delta _{I}>\Delta p_{I}\)

, тогда

\(\sigma_{I}=E_{I}(\Delta _{I}−\Delta p_{I})\)

.

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

\(\sigma_{I}=E_{I}(1−D)(\Delta _{I}−\Delta p_{I})\)

.

Для каждого направления 1 и 2 в

плоскость сдвига.

\(\sigma_{II−1}=E_{II}(1−D)(\Delta _{II−1}−\Delta p_{II−1})\) \(\sigma_{II−2}=E_{II}(1−D)(\Delta _{II−2}−\Delta p_{II−2})\)

11. Соединительный элемент – это

    удалено, когда

\(\Delta _{m}>\delta_{mf}\) .