/MAT/LAW121 (PLAS_RATE)

Ключевое слово в формате блока Упругопластический материал, зависящий от скорости деформации: с изотропным критерием текучести фон Мизеса. Этот материальный закон доступен как для твердых тел, так и для ракушки.

Формат


/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID	/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/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	E	\(\nu\)	\(\nu\)	Ирес	Ивиск	Fcut	Fcut	DTMIN	DTMIN
Fct_SIG0		Xscale_SIG0	Xscale_SIG0	Yscale_SIG0	Yscale_SIG0
Fct_YOUN		Xscale_YOUN	Xscale_YOUN	Yscale_YOUN	Yscale_YOUN
Fct_TANG		Xscale_TANG	Xscale_TANG	TANG	TANG
Fct_FAIL	Если не удалось	Xscale_FAIL	Xscale_FAIL	Yscale_FAIL	Yscale_FAIL

Определение

Поле	Содержание	Пример единицы СИ
mat_ID	Идентификатор материала.(Целое число, максимум 10 цифр)
unit_ID	(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)
mat_title	Название материала.(Персонаж, максимум 100 символов)
\(\rho_{i}\)	Начальный плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
E	Янг модуль.(Реальный)	\([Pa]\)
\(\nu\)	Пуассона соотношение.(Реальное)
Ирес	Метод разрешения пластичности. = 0 Установите на 2. = 1 NICE (ошибка следующего приращения) явный метод. = 2 (по умолчанию) Итерационный полунеявный метод Ньютона (Cutting самолет). (Целое число)
Ивиск	Формулировка зависимости скорости деформации. = 0 (по умолчанию) Масштабированный предел текучести. = 1 Вязко-пластическая рецептура. (Целое число)
Если не удалось	Переменная критерия отказа (используется только Fct_FAIL > 0). = 0 (по умолчанию) Эквивалентный стресс. =1 Эффективный пластический штамм. = 2 Максимальное главное напряжение или абсолютное значение минимальное главное напряжение. = 3 Максимальное главное напряжение. (Целое число)
Fcut	Частота среза для скорости деформации фильтрация (использовал только Ivisc= 0).По умолчанию = 10000 Гц (реальная)	\([Hz]\)
DTMIN	Минимальный размер шага по времени для автоматическое удаление элемента. По умолчанию = 0,0 (Реальное)	\([s]\)
Fct_SIG0	Предел текучести по сравнению с эффективным идентификатор функции скорости деформации. (Целое число)
Xscale_SIG0	Масштабный коэффициент для абсцисс (деформация ставка) в Fct_SIG0.По умолчанию = 1,0. (Реал)	\([\frac{1}{s}]\)
Yscale_SIG0	Масштабный коэффициент для ординаты (напряжения) в Fct_SIG0.По умолчанию = 1,0. (Реал)	\([Pa]\)
Fct_YOUN	Модуль Юнга в сравнении с эффективным идентификатор функции скорости деформации (доступен только Ivisc = 0).По умолчанию = 1.0 (Реал)
Xscale_YOUN	Масштабный коэффициент для абсцисс (деформация ставка) в Fct_YOUN.По умолчанию = 1,0. (Реал)	\([\frac{1}{s}]\)
Yscale_YOUN	Масштабный коэффициент для ординаты (напряжения) в Fct_YOUN.По умолчанию = 1,0. (Реал)	\([Pa]\)
Fct_TANG	Касательный модуль пластичности в зависимости от эффективная функция скорости деформации идентификатор.(Целое число)
Xscale_TANG	Масштабный коэффициент для абсцисс (деформация ставка) в Fct_TANG.По умолчанию = 1,0. (Реал)	\([\frac{1}{s}]\)
TANG	Если Fct_TANG > 0 Масштабный коэффициент для ординаты (напряжения) в Fct_TANG. По умолчанию = 1,0 (реальное) Если Fct_TANG = 0 Пластическая касательная TANG. По умолчанию = 0,0 (реальное)	\([Pa]\)
Fct_FAIL	Переменная критерия отказа по сравнению с эффективная функция скорости деформации идентификатор.(Целое число)
Xscale_FAIL	Масштабный коэффициент для абсцисс (деформация ставка) в Fct_FAIL.По умолчанию = 1,0. (Реал)	\([\frac{1}{s}]\)
Yscale_FAIL	Масштабный коэффициент для ординаты (напряжения) в Fct_FAIL.По умолчанию = 1,0. (Реал)	\([Pa]\)

Пример (Сталь)

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

/UNIT/1

Unit system

                  Mg                  mm                   s

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

/MAT/LAW121/1/1

Steel

#              RHO_I

             7.85E-9

#                  E                  Nu      Ires     Ivisc                Fcut               DTMIN

            194200.0                 0.3         1         0                1000                   0

# Fct_SIG0                   Xscale_SIG0         Yscale_SIG0

         5                             0                   0

# Fct_YOUN                   Xscale_YOUN         Yscale_YOUN

         6                             0                   0

# Fct_TANG                   Xscale_TANG                TANG

         7                             0                   0

# Fct_FAIL     Ifail         Xscale_FAIL         Yscale_FAIL

         8         1                   0                   0

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

/FUNCT/5

Yield stress versus effective strain rate

                 0.0               272.2

              1000.0               572.2

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

/FUNCT/6

Young’s modulus versus effective strain rate

                 0.0            192400.0

              1000.0            272400.0

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

/FUNCT/7

Plastic tangent modulus versus effective strain rate

                 0.0              4500.0

              1000.0              6500.0

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

/FUNCT/8

Failure criterion variable versus effective strain

                 0.0                0.15

              1000.0                0.25

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

#ENDDATA

/END

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

Комментарии

Этот материальный закон учитывает изотропную и линейную упругость.
Пластическое поведение описывается классической формулой фон Мизеса.

функция текучести пластичности, описываемая:

\(f=\sigma_{VM}−\sigma_{Y}\) С: \(\sigma_{VM}=\sqrt{\frac{3}{2}s:s}\) Где,

\(s\) это девиатор тензор напряжений.

Предел текучести зависит от скорости деформации и описывается

как:

\(\sigma_{Y}=\sigma_{0}(\dot{\epsilon})+\frac{E(\dot{\epsilon})E_{t}(\dot{\epsilon})}{E(\dot{\epsilon})−E_{t}(\dot{\epsilon})}\epsilon_{p}\) Где, \(\sigma_{0}\) Зависящий от скорости деформации начальный предел текучести, определяемый формулой

функция

Fct_SIG0 . \(E(\dot{\epsilon})\) Модуль Юнга, который может опционально зависеть от деформации

ставка (если

Fct_YOUN определяется). \(E_{t}(\dot{\epsilon})\) Касательный модуль пластичности, который может опционально зависеть от

скорость деформации (если

Fct_TANG есть

определено).

\(\epsilon_{p}\) Кумулятивная пластическая деформация. .. note:
Неравенство :math:`E(\dot{\epsilon})>E_{t}(\dot{\epsilon})`
 всегда должно выполняться. Если нет, решатель автоматически
 ограничить значение касательного модуля
 :math:`E_{t}(\dot{\epsilon})=0.99\timesE(\dot{\epsilon})`
.

Доступны две формулы зависимости скорости деформации в зависимости от

стоимость флага

Ивиск . - Ивиск

= 0 : масштабированный предел текучести

формулировка. В этом случае скорость деформации соответствует полной

эффективная скорость деформации рассчитывается как:

\(\dot{\epsilon}=\sqrt{\frac{2}{3}\dot{\epsilon}':\dot{\epsilon}'}\) Где,

\(\dot{\epsilon}'\) - девиаторная полная деформация тензор скорости.

Для этой формулировки все скорости деформации

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

\(\dot{\epsilon}_{_{f}}^{n}=\alpha\dot{\epsilon}^{n}+(1−\alpha)\dot{\epsilon}_{f}^{n−1}\) Где,

\(\alpha=2\piF_{cut}\Delta t\)

.

Затем вы сможете адаптировать

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

Ивиск = 1 : вязкопластик

формулировка. В этом случае скорость деформации соответствует

Скорость пластической деформации.

\(\dot{\epsilon}=\dot{\epsilon}_{p}\) В этом случае скорость деформации изначально равна

null, вычисляется и обновляется во время сопоставления возврата процедура. Все изменения переменных, зависящих от скорости деформации, приниматься во внимание, что приводит к увеличению стоимости вычислений. Однако фильтрация скорости деформации не требуется, и вы можете ускорить запустить симуляцию, используя явное отображение возврата NICE (Ирес = 1). Эта формулировка Юнга Модуль не может изменяться в зависимости от скорости деформации.

Вы можете выбрать алгоритм отображения возврата, чтобы найти компромисс.

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

Ires = 1: NICE явный

метод (ошибка следующего приращения). Эта процедура требует только 1 итерация и использует способность самокорректировки, чтобы по-прежнему сохранять хорошая точность. Этот метод позволяет значительно ускорить моделирование.

Ires = 2: секущая плоскость

(Итерационный метод Ньютона). Эта процедура требует нескольких итерации (обычно от 3 до 5) для решения поведения материала нелинейные уравнения. Это дороже, чем явное NICE метод, но обеспечивает очень хорошую точность. Этот метод установлен по умолчанию.

Вы можете добавить сбой в поведение материала. Для этого

функция

Fct_FAIL должны быть определены для описания

эволюция переменной критерия разрушения со скоростью деформации. Природа

переменной критерия отказа выбирается с помощью флага

Если не удалось : - Ifail = 0: критерий

выполняется, когда напряжение по Мизесу достигает критического значения определено в функции

Ifail = 1: критерий

выполняется, когда совокупная пластическая деформация достигает значения определено в функции

Ifail = 2: критерий

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

Ifail = 3: критерий

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

Автоматическое удаление элемента также можно настроить, если временной шаг элемента

уменьшается до достижения критического значения, определяемого DTMIN.