/MAT/LAW60 (PLAS_T3)

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

Формат


/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/mat_ID/unit_ID	/MAT/LAW60/mat_ID/unit_ID or /MAT/PLAS_T3/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\)	\(\epsilon_{p}^{max}\)	\(\epsilon_{p}^{max}\)	\(\epsilon_{t}\)	\(\epsilon_{t}\)	\(\epsilon_{m}\)	\(\epsilon_{m}\)
Функция	Фгладкий	Мангольд	Мангольд	Fcut	Fcut
fct_IDp	Fшкала	Fшкала	fct_IDE	Эйнф	Эйнф	CE	CE
fct_ID1	fct_ID2	fct_ID3	fct_ID4	fct_ID5

Читать только если 6 ≤ N функция ≤ 10 .. csv-table:

:header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
:widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

"fct_ID6", "fct_ID7", "fct_ID8", "fct_ID9", "fct_ID10", "", "", "", "", ""

Всегда читайте .. csv-table:

:header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
:widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

"Fмасштаб1", "Fмасштаб1", "Fмасштаб2", "Fмасштаб2", "Fшкала3", "Fшкала3", "Fшкала4", "Fшкала4", "Fшкала5", "Fшкала5"

Читать только если 6 < N функция < 10 .. csv-table:

:header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
:widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

"Fмасштаб6", "Fмасштаб6", "Fscale7", "Fscale7", "Fшкала8", "Fшкала8", "Fшкала9", "Fшкала9", "Fшкала10", "Fшкала10"

Всегда читайте .. csv-table:

:header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
:widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

":math:`\dot{\epsilon}_{1}`", ":math:`\dot{\epsilon}_{1}`", ":math:`\dot{\epsilon}_{2}`", ":math:`\dot{\epsilon}_{2}`", ":math:`\dot{\epsilon}_{3}`", ":math:`\dot{\epsilon}_{3}`", ":math:`\dot{\epsilon}_{4}`", ":math:`\dot{\epsilon}_{4}`", ":math:`\dot{\epsilon}_{5}`", ":math:`\dot{\epsilon}_{5}`"

Читать только если 6 ≤ N функция ≤ 10 .. csv-table:

:header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
:widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

":math:`\dot{\epsilon}_{6}`", ":math:`\dot{\epsilon}_{6}`", ":math:`\dot{\epsilon}_{7}`", ":math:`\dot{\epsilon}_{7}`", ":math:`\dot{\epsilon}_{8}`", ":math:`\dot{\epsilon}_{8}`", ":math:`\dot{\epsilon}_{9}`", ":math:`\dot{\epsilon}_{9}`", ":math:`\dot{\epsilon}_{10}`", ":math:`\dot{\epsilon}_{10}`"

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

Поле	Содержание	Пример единицы СИ
mat_ID	Материал идентификатор.(Целое число, максимум 10 цифр)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
mat_title	Материал заголовок.(Символ, максимум 100 символов)
\(\rho_{i}\)	Начальный плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
E	Янг модуль.(Реальный)	\([Pa]\)
\(\nu\)	Пуассона соотношение.(Реальное)
\(\epsilon_{p}^{max}\)	Неудачный пластик штамм.По умолчанию = 1,0 × 1030 (реальный)
\(\epsilon_{t}\)	Деформация разрушения при растяжении какое напряжение начинает уменьшаться. По умолчанию = 1,0 x 1030 (Реал)
\(\epsilon_{m}\)	Максимальное разрушение при растяжении напряжение, при котором элемент удаляется. По умолчанию = 2,0 x 1030 (Реал)
Функция	Количество функций. Это должно быть 4 < Nфункция < 10. По умолчанию ≤ 10 (целое число).
Фгладкий	Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)
Мангольд	Коэффициент закалки. = 0 Закалка представляет собой полную изотропную модель. = 1 При закалке используется кинематика Прагера-Циглера. модель. = значение от 0 до 1 Упрочнение интерполируется между двумя модели. (Настоящий)
Fcut	Частота среза для фильтрация скорости деформации. 7По умолчанию = 1,0 × 1030 (Реал)	\([Hz]\)
fct_IDp	Давление в зависимости от урожайности Факторная функция. 9По умолчанию = 0 (Целое число)
Fшкала	Масштабный коэффициент доходности коэффициент в fct_IDp.По умолчанию = 1,0 (реальный)	\([Pa]\)
fct_IDE	Идентификатор функции для масштабный коэффициент модуля Юнга, когда модуль Юнга равен Функция пластической деформации. 6По умолчанию = 0: в этом случае эволюция модуля Юнга зависит от Einf. и CE.(Целое число)
Эйнф	Насыщенный модуль Юнга для инфинитива пластического напряжения.(Настоящее)	\([Pa]\)
CE	Параметр для Янга эволюция модуля.(Реальный)
fct_ID1	Функция напряжения текучести идентификатор 1, соответствующий скорости деформации \(\dot{\epsilon}_{1}\) .(Целое число)
fct_ID2	Функция напряжения текучести идентификатор 2, соответствующий скорости деформации \(\dot{\epsilon}_{2}\) .(Целое число)
fct_ID3	Функция напряжения текучести идентификатор 3, соответствующий скорости деформации \(\dot{\epsilon}_{3}\) .(Целое число)
fct_ID4	Функция напряжения текучести идентификатор 4, соответствующий скорости деформации \(\dot{\epsilon}_{4}\) .(Целое число)
fct_ID5	Функция напряжения текучести идентификатор 5, соответствующий скорости деформации \(\dot{\epsilon}_{5}\) .(Целое число)
fct_ID6	Функция напряжения текучести идентификатор 6, соответствующий скорости деформации \(\dot{\epsilon}_{6}\) .(Целое число)
fct_ID7	Функция напряжения текучести идентификатор 7, соответствующий скорости деформации \(\dot{\epsilon}_{7}\) .(Целое число)
fct_ID8	Функция напряжения текучести идентификатор 8, соответствующий скорости деформации \(\dot{\epsilon}_{8}\) .(Целое число)
fct_ID9	Функция напряжения текучести идентификатор 9, соответствующий скорости деформации \(\dot{\epsilon}_{9}\) .(Целое число)
fct_ID10	Функция напряжения текучести идентификатор 10, соответствующий скорости деформации \(\dot{\epsilon}_{10}\) .(Целое число)
Fмасштаб1	Масштабный коэффициент для ординаты (напряжение) в fct_ID1.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fмасштаб2	Масштабный коэффициент для ординаты (напряжение) в fct_ID2.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fшкала3	Масштабный коэффициент для ординаты (напряжение) в fct_ID3.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fшкала4	Масштабный коэффициент для ординаты (напряжение) в fct_ID4.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fшкала5	Масштабный коэффициент для ординаты (напряжение) в fct_ID5.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fмасштаб6	Масштабный коэффициент для ординаты (напряжение) в fct_ID6.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fscale7	Масштабный коэффициент для ординаты (напряжение) в fct_ID7.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fшкала8	Масштабный коэффициент для ординаты (напряжение) в fct_ID8.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fшкала9	Масштабный коэффициент для ординаты (напряжение) в fct_ID9.По умолчанию = 1,0 (реальное)	\([Pa]\)
Fшкала10	Масштабный коэффициент для ординаты (напряжение) в fct_ID10.По умолчанию = 1,0 (реальное)	\([Pa]\)
\(\dot{\epsilon}_{1}\)	Скорость деформации 1.(Реальный)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{2}\)	Скорость деформации 2.(Реальный)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{3}\)	Скорость деформации 3.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{4}\)	Скорость деформации 4.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{5}\)	Скорость деформации 5.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{6}\)	Скорость деформации 6.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{7}\)	Скорость деформации 7.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{8}\)	Скорость деформации 8.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{9}\)	Скорость деформации 9.(Реал)	\([\frac{1}{s}]\)
\(\dot{\epsilon}_{10}\)	Скорость деформации 10.(Реал)	\([\frac{1}{s}]\)

Пример (алюминий)

#RADIOSS STARTER

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

/UNIT/1

unit for mat

                   g                  mm                  ms

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

#-  2. MATERIALS:

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

/MAT/LAW60/1/1

Aluminium_example

#              RHO_I

               .0027

#                  E                  Nu           Eps_p_max               Eps_t               Eps_m

               60400                 .33                   0                   0                   0

#  N_funct  F_smooth              C_hard               F_cut

         4         0                   0                   0

#  fct_IDp              Fscale   Fct_IDE                EInf                  CE

         0                   0         0                   0                   0

# Functions

         1         2         3         4

# Scale factors          Fscale_5

                   1                 1.2                 1.4                 1.6

# Strain rates          Eps_dot_5

                   0                  20                  30                  40

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

#-  3. FUNCTIONS:

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

/FUNCT/1

function_36

#                  X                   Y

                   0                  90

              2.5E-4                 100

                .001                 104

                .009                 121

                .017                 136

                .021                 143

                .036                 156

                .045                 162

                .055                 165

                .072                 170

                .075                 170

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

/FUNCT/2

function_36

#                  X                   Y

                   0                  90

              2.5E-4                 100

                .001                 104

                .009                 121

                .017                 136

                .021                 143

                .036                 156

                .045                 162

                .055                 165

                .072                 170

                .075                 170

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

/FUNCT/3

function_36

#                  X                   Y

                   0                  90

              2.5E-4                 100

                .001                 104

                .009                 121

                .017                 136

                .021                 143

                .036                 156

                .045                 162

                .055                 165

                .072                 170

                .075                 170

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

/FUNCT/4

function_36

#                  X                   Y

                   0                  90

              2.5E-4                 100

                .001                 104

                .009                 121

                .017                 136

                .021                 143

                .036                 156

                .045                 162

                .055                 165

                .072                 170

                .075                 170

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

#ENDDATA

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

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

Первая точка выхода

функции напряжений (пластическая деформация против напряжения) должны иметь пластическую деформацию

значение ноль. Если последняя точка первой (статической) функции равна 0 в

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

\(\epsilon_{p}^{max}\) устанавливается на соответствующее значение \(\epsilon_{p}\) .

If

\(\epsilon_{p}\) (пластическая деформация) достигает \(\epsilon_{p}^{max}\) , в одной точке интеграции элемент удаляется.

Если (наибольшая основная деформация)

\(\epsilon_{1}>\epsilon_{t}\) , стресс снижается с помощью: \(\sigma=\sigma(\frac{\epsilon_{m}−\epsilon_{1}}{\epsilon_{m}−\epsilon_{t}})\)

If

\(\epsilon_{1}>\epsilon_{m}\) , элемент удаляется.

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

в глобальной формулировке отсутствует (закалка полностью изотропна).

Для кинематической закалки и

зависимость от скорости деформации, предел текучести зависит от скорости деформации.

Входной сигнал фильтрации скорости деформации

(Fcut) доступен только для оболочек и твердотельных элементов.

Используется фильтрация по скорости деформации.

для сглаживания скорости деформации.

fct_ID

p используется для различения поведения при растяжении и сжатии определенных

материалов (т. е. выход, зависящий от давления). Это доступно для твердых

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

номинальный предел текучести на коэффициент текучести, соответствующий фактическому

давление.

If

\(\dot{\epsilon}_{n}\le\dot{\epsilon}\le\dot{\epsilon}_{n+1}\) , предел текучести представляет собой кубическую интерполяцию между

функции

f n-1 , f n , f n+1 и f n+2 .

If

\(\dot{\epsilon}\le\dot{\epsilon}_{1}\) , предел текучести экстраполируется между функциями \(f_{1}\) , \(f_{2}\) и \(f_{3}\) .

If

\(\dot{\epsilon}_{Nfunc−1}\le\dot{\epsilon}\le\dot{\epsilon}_{Nfunc}\) , доходность интерполируется между функциями f Нфунк-3 , f Нфунк-2 , f Нфунк-1 и f Нфунк .

If

\(\dot{\epsilon}>\dot{\epsilon}_{Nfunc}\) , доходность экстраполируется между функциями f Нфунк-2 , f Нфунк-1 и f Нфунк . .. image:: images/mat_law60_plas_t3_starter_r_mat_law60_yield.png

alt

mat_law60_yield

(Рисунок 2.)

Функции, описывающие деформацию

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

Значения скорости деформации должны быть

даны строго в порядке возрастания.

Эволюция Янга

модуль:

If fct_ID E > 0, кривая определяет масштабный коэффициент для эволюции модуля Юнга.

с эквивалентной пластической деформацией, что означает, что модуль Юнга равен

масштабируется функцией

\(f(\bar{\epsilon}_{p})\) : - \(E(t)=E⋅f(\bar{\epsilon}_{p})\)

Начальное значение шкалы

коэффициент должен быть равен 1 и он уменьшается.

If fct_ID E = 0, модуль Юнга рассчитывается как: \(E(t)=E−(E−E_{inf})[1−exp(−C_{E}\bar{\epsilon}_{p})]\)

Где, E и E инфа соответственно

начальное и асимптотическое значение модуля Юнга и

\(\bar{\epsilon}_{p}\) это накопленный эквивалент пластика

напряжение.

Note

Если fct_IDE = 0 и CE = 0, то Юнга
модуль E остается постоянным.