/MAT/LAW74
- Ключевое слово формата блока Этот закон описывает ортотропную трехмерную модель Thermal Hill.
материал и применим только к твердым элементам. Предел текучести может зависеть от скорости деформации, или как от скорости деформации, так и от температуры.
Формат
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/mat_ID/unit_ID |
/MAT/LAW74/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 |
\(\upsilon\) |
\(\upsilon\) |
\(\epsilon_{p}^{max}\) |
\(\epsilon_{p}^{max}\) |
\(\epsilon_{t}\) |
\(\epsilon_{t}\) |
\(\epsilon_{m}\) |
\(\epsilon_{m}\) |
fct_IDE |
Эйнф |
Эйнф |
CE |
CE |
|||||
Фгладкий |
Мангольд |
Мангольд |
Fcut |
Fcut |
|||||
\(\sigma_{11}^{y}\) |
\(\sigma_{11}^{y}\) |
\(\sigma_{22}^{y}\) |
\(\sigma_{22}^{y}\) |
\(\sigma_{33}^{y}\) |
\(\sigma_{33}^{y}\) |
||||
\(\sigma_{12}^{y}\) |
\(\sigma_{12}^{y}\) |
\(\sigma_{23}^{y}\) |
\(\sigma_{23}^{y}\) |
\(\sigma_{31}^{y}\) |
\(\sigma_{31}^{y}\) |
||||
Tab_ID |
\(\sigma_{scale}\) |
\(\sigma_{scale}\) |
\(\dot{\epsilon}_{scale}\) |
\(\dot{\epsilon}_{scale}\) |
|||||
Ti |
Ti |
\(\rho_{0}C_{p}\) |
\(\rho_{0}C_{p}\) |
Определение
Поле |
Содержание |
Пример единицы СИ |
|---|---|---|
mat_ID |
Идентификатор материала.(Целое число, максимум 10 цифр) |
|
unit_ID |
Идентификатор объекта.(Целое число, максимум 10 цифр) |
|
mat_title |
Название материала.(Символ, максимум 100 символов) |
|
\(\rho_{i}\) |
Начальная плотность.(Реальная) |
\([\frac{kg}{m^{3}}]\) |
E |
Начальный Янг модуль.(Реальный) |
\([Pa]\) |
\(\upsilon\) |
Коэффициент Пуассона.(Реальный) |
|
fct_IDE |
Идентификатор функции для масштабного коэффициента модуля Юнга, когда модуль Юнга является функцией пластической деформации. = 0 (по умолчанию) В этом случае эволюция Юнга зависит от Einf и CE. (Целое число) |
|
Эйнф |
Насыщенный модуль Юнга для инфинитива пластическая деформация.(Настоящая) |
\([Pa]\) |
CE |
Параметр модуля Юнга эволюция.(Реальная) |
|
\(\epsilon_{p}^{max}\) |
Пластическая деформация отказа. По умолчанию = 1030 (Реал) |
|
\(\epsilon_{t}\) |
Деформация разрушения при растяжении, при которой напряжение начинает уменьшаться. По умолчанию = 1,0 × 1030 (Реальное) |
|
\(\epsilon_{m}\) |
Максимальная деформация разрушения при растяжении, при которой напряжение в элементе установлено на ноль. По умолчанию = 2,0 × 1030. (Реал) |
|
Фгладкий |
Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число) |
|
Мангольд |
Коэффициент закалки. = 0 Закалка – полностью изотропная модель. = 1 Для закалки используется кинематическая модель Прагера-Циглера. = между 0 и 1 Упрочнение интерполируется между двумя моделями. (Настоящий) |
|
Fcut |
Частота среза для скорости деформации фильтрация.По умолчанию = 1,0 × 1030 (Реальное) |
\([Hz]\) |
\(\sigma_{11}^{y}\) |
Выход в направлении 1.(Реальный) |
\([Pa]\) |
\(\sigma_{22}^{y}\) |
Выход в направлении 2.(Реальный) |
\([Pa]\) |
\(\sigma_{33}^{y}\) |
Выход в направлении 3.(Реал) |
\([Pa]\) |
\(\sigma_{12}^{y}\) |
Предел текучести в направлении сдвига 12.(Реал) |
\([Pa]\) |
\(\sigma_{23}^{y}\) |
Предел текучести в направлении сдвига 23.(Реал) |
\([Pa]\) |
\(\sigma_{31}^{y}\) |
Предел текучести в направлении сдвига 31.(Реал) |
\([Pa]\) |
Tab_ID |
Идентификатор таблицы предела текучести определение. 6 (целое число) |
|
\(\sigma_{scale}\) |
Масштабный коэффициент напряжения текучести. Набор по умолчанию. до 1,0 (реальный) |
\([Pa]\) |
\(\dot{\epsilon}_{scale}\) |
Масштабный коэффициент скорости деформации. Установлен по умолчанию. до 1,0 (реальный) |
\([\frac{1}{s}]\) |
Ti |
Начальная температура. По умолчанию установлено значение 293 К (Реал) |
\([K]\) |
\(\rho_{0}C_{p}\) |
Удельная теплоемкость на объем единица.(Реальная) |
\([\frac{J}{m^{3}⋅K}]\) |
Пример (алюминий)
#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/LAW74/1/1
Aluminum
# RHO_I
.0027
# E NU EPSILON_P_MAX EPSILON_T EPSILON_M
60400 .33 0 0 0
# fct_ID EINF CE
0 0 0
# FSMOOTH C_HARD FCUT
1 0 10
# SIGMA11Y SIGMA22Y SIGMA33Y
1 1 1
# SIGMA12Y SIGMA23Y SIGMA31Y
1 1 1
# TABLE SIGMA_SCALE EPSPT_SCALE
10 0 0
# TI RHO0_CP
0 0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#- 3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/TABLE/1/10
table
#DIMENSION
3
# fct_ID X Z
38 0 293
38 10 293
39 11 293
40 20 293
38 0 400
38 10 400
39 11 400
40 20 400
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/38
function_38
# X Y
0 90
2.5E-4 100
.001 104.5
.009 121
.01 136
.02 143.5
.04 163
.07 169.5
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/39
function_39
# X Y
0 108
2.5E-4 120
.001 125.4
.009 145.2
.01 163.2
.02 172.2
.04 195.6
.07 203.4
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/40
function_40
# X Y
0 126
2.5E-4 140
.001 146.3
.009 169.4
.01 190.4
.02 200.9
.04 228.2
.07 237.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
Комментарии
Этот материальный закон следует использовать с
набор свойств /PROP/TYPE6 (SOL_ORTH), /PROP/TYPE14 (SOLID), /PROP/TYPE20 (TSHELL) или /PROP/TYPE21 (TSH_ORTH).
Предел текучести определяется пользователем.
функция, а предел текучести сравнивается с эквивалентным напряжением:
\(\sigma_{eq}=\sqrt{F(\sigma_{2}−\sigma_{3})^{2}+G(\sigma_{3}−\sigma_{1})^{2}+H(\sigma_{1}−\sigma_{2})^{2}+2L\sigma_{23}^{2}+2M\sigma_{31}^{2}+2N\sigma_{12}^{2}}\) Где, HILL
- параметры:
\(F=\frac{1}{2}(\frac{1}{\sigma_{22}^{2}}+\frac{1}{\sigma_{33}^{2}}−\frac{1}{\sigma_{11}^{2}}),G=\frac{1}{2}(\frac{1}{\sigma_{11}^{2}}+\frac{1}{\sigma_{33}^{2}}−\frac{1}{\sigma_{22}^{2}}),H=\frac{1}{2}(\frac{1}{\sigma_{11}^{2}}+\frac{1}{\sigma_{22}^{2}}−\frac{1}{\sigma_{33}^{2}})\) \(L=\frac{1}{2\sigma_{23}^{2}},M=\frac{1}{2\sigma_{31}^{2}},N=\frac{1}{2\sigma_{12}^{2}}\) Где,
\(\sigma_{11},\sigma_{22},\sigma_{33},\sigma_{12},\sigma_{23}\) и \(\sigma_{31}\) представляют компоненты напряжений либо в ортотропной кадр, если используется ортотропное свойство, или в ортогональной изопараметрической рама.
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}\) , напряжение снижается до 0 (но элемент не
удалил).
Таблицу для определения предела текучести можно
быть двухмерным или трехмерным.
Если таблица двумерная, предполагается, что ее параметры представляют соответственно
пластическая деформация и скорость деформации
\((\epsilon^{p},\dot{\epsilon})\) . Тогда, если
\(\epsilon_{m−1}^{p}\le\epsilon^{p}\le\epsilon_{m}^{p}\) и \(\dot{\epsilon}_{n−1}\le\dot{\epsilon}\le\dot{\epsilon}_{n}\) доходность линейно интерполируется между четырьмя значениями таблицы, соответствующей \((\epsilon_{i}^{p},\dot{\epsilon}_{j}),i=m−1,m;j=n−1,n\)
- .
Если таблица трехмерная, предполагается, что ее параметры представляют соответственно
пластическая деформация, скорость деформации и температура
\((\epsilon^{p},\dot{\epsilon},T)\) . Тогда, если
\(\epsilon_{m−1}^{p}\le\epsilon^{p}\le\epsilon_{m}^{p}\) и \(\dot{\epsilon}_{n−1}\le\dot{\epsilon}\le\dot{\epsilon}_{n}\) и \(T_{q−1}\leT\leT_{q}\) доходность линейно интерполируется между восемью значения таблицы, соответствующие \((\epsilon_{i}^{p},\dot{\epsilon}_{j},T_{k}),i=m−1,m;j=n−1,n;k=q−1,q\)
- .
Если
\((\epsilon^{p},\dot{\epsilon})\) или \((\epsilon^{p},\dot{\epsilon},T)\) выходит за пределы таблицы, предел текучести составляет получено методом линейной экстраполяции. Таким образом, необходимо внести в таблицу статические кривые, соответствующие нулевой скорости деформации (запись \(\dot{\epsilon}=0\) должно принадлежать определению таблицы).
Ценности
таблицы — значения предела текучести.
Если предел текучести зависит еще и от температуры, то
таблица трехмерная:
Если опция /HEAT/MAT не связана с
идентификатор материала, предполагаются адиабатические условия и вычисляется температура как:
\(Τ=T_{i}+\frac{E_{int}}{\rhoC_{p}(Volume)}\) Где, E интервал Внутренняя энергия, рассчитанная по \(\rho\) , Объем Плотность и объем тока C p Теплоемкость на единицу массы В противном случае необходимо задать вопрос о формулировке метода конечных элементов для теплопередачи.
for (Iform =1 в опции /HEAT/MAT); начальный тогда температура и удельная тепловая мощность в опции /HEAT/MAT будут быть использован.
Эволюция Янга
модуль:
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 и Einf — соответственно начальный и
асимптотическое значение модуля Юнга, \(\bar{\epsilon}_{p}\) накопленная эквивалентная пластическая деформация.
Note
- Если fct_IDE = 0 и CE = 0, модуль Юнга
E остается постоянным.