/MAT/LAW43 (HILL_TAB)

Ключевое слово формата блока Этот закон описывает ортотропный материал Хилла и

применимо только к элементам оболочки. Этот закон отличается от LAW32 (HILL) только тем, что ввод предела текучести (здесь он определяется пользовательской функцией).

Формат

/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID	/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/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\)
fct_IDE		Эйнф	Эйнф	CE	CE
r00	r00	r45	r45	r90	r90	Мангольд	Мангольд	Доходность0
\(\epsilon_{p}^{max}\)	\(\epsilon_{p}^{max}\)	\(\epsilon_{t}\)	\(\epsilon_{t}\)	\(\epsilon_{m}\)	\(\epsilon_{m}\)	Fcut	Fcut	Фгладкий

Каждая функция пластичности в строке .. 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_IDi", "", "Фскалей", "Фскалей", ":math:`\dot{\epsilon}_{i}`", ":math:`\dot{\epsilon}_{i}`", "", "", "", ""

Определение

Поле	Содержание	Пример единицы СИ
mat_ID	Идентификатор материала.(Целое число, максимум 10 цифр)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
mat_title	Название материала.(Символ, максимум 100 символов)
\(\rho_{i}\)	Начальная плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
E	Модуль Юнга.(Реальный)	\([Pa]\)
\(\nu\)	Коэффициент Пуассона.(Реальный)
fct_IDE	Идентификатор функции для масштабного коэффициента модуля Юнга, когда модуль Юнга является функцией пластической деформации. 12По умолчанию = 0: в этом случае эволюция Модуль Юнга зависит от Einf и CE. (Целое число)
Эйнф	Насыщенный модуль Юнга для инфинитива пластическая деформация.(Настоящая)	\([Pa]\)
CE	Параметр эволюции модуля Юнга. 12(Реал)
r00	Параметр Ланкфорда 0 градусов. 3По умолчанию = 1,0 (реальное)
r45	Параметр Ланкфорда 45 градусов. По умолчанию = 1,0 (Реальный)
r90	Параметр Ланкфорда 90 градусов.По умолчанию = 1,0 (Реальный)
Мангольд	Коэффициент закалки. = 0 Закалка – полностью изотропная модель. = 1 Для закалки используется кинематическая модель Прагера-Циглера. = между 0 и 1 Упрочнение интерполируется между двумя моделями. (Настоящий)
Доходность0	Флаг стресса доходности. = 0 Средний входной предел текучести. = 1 Предел текучести в ортотропном направлении 1. (Целое число)
\(\epsilon_{p}^{max}\)	Пластическая деформация разрушения. По умолчанию = 1,0. × 1030 (Реал)
\(\epsilon_{t}\)	Деформация разрушения при растяжении, при которой напряжение начинает уменьшаться. По умолчанию = 1,0 × 1030 (Реальное)
\(\epsilon_{m}\)	Максимальная деформация разрушения при растяжении, при которой напряжение в элементе установлено на ноль. По умолчанию = 2,0 × 1030. (Реал)
Fcut	Частота среза для скорости деформации фильтрация.(Реальная)	\([Hz]\)
Фгладкий	Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)
fct_IDi	Кривые пластичности .ith идентификатор функции.(Целое число)
Фскалей	Масштабный коэффициент для i-я функция. По умолчанию установлено значение 1,0. (Реал)
\(\dot{\epsilon}_{i}\)	Скорость деформации для i-я функция.(Реальная)	\([\frac{1}{s}]\)

Пример (Металл)

#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----|

#-  2. MATERIALS:

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

/MAT/HILL_TAB/1/1

metal

#              RHO_I

             7.8E-09

#                  E                  NU

              206000                  .3

#FUNCT_IDE                          EINF                  CE

         0                             0                   0

#                r00                 r45                 r90              C_hard   Iyield0

                1.73                1.34                2.24                   0         0

#           EPSP_max              EPS_t1               EPS_m                Fcut   Fsmooth

                   0                   0                   0               15000         1

# func_IDi                      Fscale_i           EPS_dot_i

         5                             0                   0

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

#-  3. FUNCTIONS:

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

/FUNCT/5

metal

#                  X                   Y

                   0                 260

                .002                 270

                .005                 280

                 .01                 297

                 .02                 322

                 .05                 370

                  .1                 422

                 .15                 457

                  .2                 485

                  .3                 528

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

#ENDDATA

/END

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

Комментарии

1. Этот материальный закон необходимо использовать

с набором свойств /PROP/TYPE9 (SH_ORTH) или /PROP/TYPE10 (SH_COMP).

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

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

\(\sigma_{eq}=\sqrt{A_{1}\sigma_{1}^{2}+A_{2}\sigma_{2}^{2}−A_{3}\sigma_{1}\sigma_{2}+A_{12}\sigma_{12}^{2}}\)

3. Углы для параметров Ланкфорда равны

   определяется относительно ортотропного направления 1.

\(R=\frac{r_{00}+2r_{45}+r_{90}}{4}H=\frac{R}{1+R}A_{1}=H(1+\frac{1}{r_{00}})A_{2}=H(1+\frac{1}{r_{90}})A_{3}=2HA_{12}=2H(r_{45}+0.5)(\frac{1}{r_{00}}+\frac{1}{r_{90}})r_{00}=\frac{A_{3}}{2A_{1}−A_{3}}r_{45}=\frac{1}{2}(\frac{A_{12}}{A_{1}+A_{2}−A_{3}}−1)r_{90}=\frac{A_{3}}{2A_{2}−A_{3}}\) Параметры Ланкфорда

ra – коэффициент пластической деформации в плоская и пластическая деформация по толщине \(\epsilon_{33}\)

.

\(r_{\alpha}=\frac{d\epsilon_{\alpha+\pi/2}}{d\epsilon_{33}}\) Где α — угол к ортотропному направлению 1. Это

Параметры Лэнкфорда ra могут быть определяется путем простого испытания на растяжение под углом α.

Более высокое значение

R означает лучшую формуемость.

4. Если последняя точка первой

   (статическая) функция равна 0 в стрессе, значение по умолчанию

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

5. Удаление элемента:

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

• If \(\epsilon_{1}\) достигает \(\epsilon_{t}\) , напряжение уменьшается по следующему соотношению: \(\sigma=\sigma(\frac{\epsilon_{m}−\epsilon_{1}}{\epsilon_{m}−\epsilon_{t}})\)

• If \(\epsilon_{1}\) (крупнейший основной штамм) достигает \(\epsilon_{m}\) ( \(\epsilon_{1}>\epsilon_{m}\) ), напряжение в элементе снижается до 0 (но элемент

  не удаляется).

• Один раз \(\epsilon_{1}\) (крупнейший основной штамм) достигает \(\epsilon_{f}\) (максимальная деформация разрушения при растяжении), элемент

  удален.

6. Максимальное количество кривых, которые

можно ввести 10.

7. If

\(\dot{\epsilon}\le\dot{\epsilon}_{n}\) , доходность интерполируется между \(f_{n}\) и \(f_{n−1}\) .

8. If

\(\dot{\epsilon}\le\dot{\epsilon}_{1}\) , функция \(f_{1}\) используется.

9. Выше

\(\dot{\epsilon}_{max}\) , доходность экстраполируется. .. image:: images/mat_law43_hill_tab_starter_r_mat_law43_yield.png

alt

mat_law43_yield

(Рисунок 1.)

10. Радиальный возврат недоступен.

(только итерационная пластичность).

11. Если предел текучести был

полученное в ортотропном направлении 1, определим Доходность0 =1; иначе Доходность0 =0.

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

    модуль:

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

  эволюция с эквивалентной пластической деформацией

  \(\bar{\epsilon}_{p}\) , что означает, что модуль Юнга масштабируется по

  функция

  \(f(\bar{\epsilon}_{p})\) : \(E(t)=f(\bar{\epsilon}_{p})E\) Начальное значение масштабного коэффициента должно быть равно 1 и оно

уменьшается.

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

\(E\) и \(E_{inf}\) являются соответственно начальным и асимптотическим значением модуль Юнга и \(\bar{\epsilon}_{p}\) – накопленная эквивалентная пластическая деформация.

Note

Если fct_IDE = 0 и CE = 0, то Юнга

модуль E остается постоянным.

13. Параметры

F гладкий и F резать позволяют включить фильтрацию по скорости деформации. Можно установить три случая: - Если Fsmooth = 0 + Fcut = 0,0, фильтрация скорости деформации отключается. - Если Fsmooth = 1 + Fcut = 0,0, фильтрация скорости деформации использует частоту среза по умолчанию 10 кГц. - Если Fcut ≠ 0, Fsmooth автоматически устанавливается на 1, а фильтрация скорости деформации использует частоту среза.

предоставленный вами.