/MAT/LAW76 (SAMP)

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

Формат


/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID	/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/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\)
tab_IDt	tab_IDc	tab_IDs
Fшкалет	Fшкалет	Фскалек	Фскалек	Fмасштабы	Fмасштабы			XFAC	XFAC
\(\nu_{p}\)	\(\nu_{p}\)	fct_IDpr	Fscalepr	Fscalepr	Фгладкий	Fcut	Fcut
\(\epsilon_{p}^{f}\)	\(\epsilon_{p}^{f}\)	\(\epsilon_{p}^{r}\)	\(\epsilon_{p}^{r}\)
fct_ID1			Fмасштаб1	Fмасштаб1
Яформа	IQUAD	ICONV

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

Поле	Содержание	Пример единицы СИ
mat_ID	Материал идентификатор.(Целое число, максимум 10 цифр)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
mat_title	Материал заголовок.(Символ, максимум 100 символов)
\(\rho_{i}\)	Начальный плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
E	Начальный Янг модуль.(Реальный)	\([Pa]\)
\(\nu\)	Пуассона соотношение.(Реальное)
tab_IDt	Таблица пределов текучести при растяжении идентификатор (напряжение в сравнении с пластической деформацией растяжения с возможность зависимости скорости деформации).(Целое число)
tab_IDc	Предел текучести при сжатии идентификатор таблицы (напряжение в сравнении с пластической деформацией сжатия с возможность зависимости скорости деформации).(Целое число)
tab_IDs	Таблица пределов текучести при сдвиге идентификатор (напряжение в сравнении с пластической деформацией сдвига с возможностью зависимости скорости деформации).(Целое число)
Fшкалет	Масштабный коэффициент для ординаты (стресс) для tab_IDt.По умолчанию = 1,0 (Реальный)	\([Pa]\)
Фскалек	Масштабный коэффициент для ординаты (стресс) для tab_IDc.По умолчанию = 1,0 (Реальный)	\([Pa]\)
Fмасштабы	Масштабный коэффициент для ординаты (стресс) для tab_IDs.По умолчанию = 1,0 (Реальный)	\([Pa]\)
XFAC	Масштабный коэффициент для вторая запись (скорость деформации) из трех таблиц (tab_IDt, tab_IDc и tab_IDs). 6По умолчанию = 1,0 (реальное)
\(\nu_{p}\)	Пластиковый Пуассон соотношение.(Реальное)
fct_IDpr	Пластиковый коэффициент Пуассона идентификатор функции ( \(\nu_{p}\) против пластика штамм).(Реальный)
Fscalepr	Масштабный коэффициент для ординаты ( \(\nu_{p}\) ) в fct_IDpr.По умолчанию = 1,0 (Реальный)
Фгладкий	Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)
Fcut	Частота среза для фильтрация скорости деформации. По умолчанию = 1030 (Реал)	\([Hz]\)
\(\epsilon_{p}^{f}\)	Неудачная пластическая деформация (начало элементарного урона). По умолчанию = 2e30 (Реальное)
\(\epsilon_{p}^{r}\)	Максимальная пластическая деформация (элемент удален). По умолчанию = 2e30 (реальный)
fct_ID1	Идентификатор функции повреждения (повреждение в зависимости от пластической деформации). 2 (целое число)
Fмасштаб1	Масштабный коэффициент для ординаты для fct_ID1. 2По умолчанию = 1,0 (Реальное)
Яформа	Флаг формулировки. 4 = 0 (по умолчанию) Нет соответствующей формулировки. = 1 Связанная с ним формулировка фон Мизеса. (Целое число)
IQUAD	Флаг поверхности текучести. 3 = 0 (по умолчанию) Поверхность текучести линейна по фон Мизесу. = 1 Поверхность текучести квадратична по Мизесу. (рекомендуется). (Целое число)
ICONV	Флаг условия выпуклости. =0 (по умолчанию) Никакой обработки для обеспечения стабильности материала. = 1 Выпуклость предела текучести (стабильность материала) равна уверен. (Целое число)

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

#RADIOSS STARTER

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

/UNIT/1

unit for mat

                  kg                  mm                  ms

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

/MAT/LAW76/1/1

LAW76_Material

#              RHO_I

                1E-6

#                  E                  nu

               100.0                  .3

#  TAB_IDt   TAB_IDc   TAB_IDs

      1000      1001      1003

#           Fscale_t            Fscale_c            Fscale_s                                    XFAC

               1.000               1.000               1.000                                   1.000

#               Nu_p  fct_IDpr           Fscale_pr   Fsmooth      Fcut

                 0.5         0                   0         1      1e30

#            EPS_f_p             EPS_r_p

                   0                   0

#funct_ID1                                Fscale_1

         0

#    IFORM     IQUAD     ICONV

         0         0         1

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

/TABLE/1/1000

curve_list TENSION strain rates

         2

     10010                        1.0e-4

     10020                           1.0

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

/FUNCT/10010

eps_vs_sigma funct dt=1.0e-4

              0.0000             .100000

              1.0000             .200000

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

/FUNCT/10020

eps_vs_sigma funct dt=1.0e-4

              0.0000             .100000

              1.0000             .200000

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

/TABLE/1/1001

curve_list COMPRESSION strain rates

         2

     10030                        1.0e-4

     10040                           1.0

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

/FUNCT/10030

eps_vs_sigma funct dt=1.0e-4

              0.0000             .200000

              1.0000             .400000

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

/FUNCT/10040

eps_vs_sigma funct dt=1.0e-4

              0.0000             .200000

              1.0000             .400000

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

/TABLE/1/1003

curve_list SHEAR strain rates

         2

     10050                        1.0e-4

     10060                           1.0

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

/FUNCT/10050

eps_vs_sigma funct dt=1.0e-4

              0.0000             .050000

              0.5000             .060000

              1.0000             .065000

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

/FUNCT/10060

eps_vs_sigma funct dt=1.0e-4

              0.0000             .050000

              0.5000             .060000

              1.0000             .065000

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

#ENDDATA

/END

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

Это

материал совместим с панцирем, толстой оболочкой и твердыми элементами.

Материальный ущерб может быть

смоделировал двумя способами:

\(\epsilon_{p}^{f}\) (начало повреждать) и \(\epsilon_{p}^{r}\) (удаление элемента): \(D=\frac{\epsilon_{p}-\epsilon_{p}^{f}}{\epsilon_{p}^{r}-\epsilon_{p}^{f}}\) Где,

\(\epsilon_{p}\) это текущий накопленный пластик напряжение.

Функция повреждения fct_ID 1 .. image:: images/mat_law76_samp_starter_r_fail_tab1_damage.png

(Рисунок 1.)

Если используется функция повреждения fct_ID1, то

\(\epsilon_{p}^{f}\) и \(\epsilon_{p}^{r}\) будет проигнорирован.

Ущерб

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

\(\sigma=\sigma_{eff}1-D\) Где, \(\sigma_{eff}\) Неповрежденное (эффективное) тензорное напряжение \(\sigma\) Поврежденный (текущий) тензор напряжений \(D\) Переменная урона

Выбор поверхности текучести:

\(f={\sigma_{VM}−A_{0}−A_{1}P−A_{2}P^{2}IQUAD=0\sigma_{VM}^{2}−A_{0}−A1P−A_{2}P^{2}IQUAD=1\) Где, \(P=−\frac{\sigma_{xx}+\sigma_{yy}+\sigma_{zz}}{3}\) \(\sigma_{VM}=\sqrt{\frac{3}{2}[(\sigma_{xx}+P)^{2}+(\sigma_{yy}+P)^{2}+(\sigma_{zz}+P)^{2}+2\sigma_{xy}^{2}+2\sigma_{yz}^{2}+2\sigma_{xz}^{2}]}\) \(A_{0}\)

,

\(A_{1}\) и \(A_{2}\) коэффициенты рассчитываются по упрочнению кривая растяжения, сжатия и сдвига.

Для фон Мизеса и

Поверхность текучести Друкера-Прагера, IQUAD=0 можно использовать. Однако в некоторых ситуациях Radioss может быть сложно адаптироваться к \(A_{0}\)

,: \(A_{1}\) и \(A_{2}\) коэффициенты при использовании IQUAD=0, и более простая подгонка получено с использованием IQUAD=1.

Выбор пластика

формулировка:

Для формулировки, не связанной с пластичностью, I форма =0: Функция правила пластического потока,

\(g\)

, используется для описания пластической деформации: приращение \(d\epsilon_{p}=d\lambda\frac{\partialg}{\partial\sigma}\)
. В этом случае: \(\frac{\partialg}{\partial\sigma}\) не является нормальным к поверхности текучести \(f\) и \(g\) не связан с поверхностью текучести \(f\)

.

Правило пластического течения,

\(g\)

, определяется:

\(g=\sqrt{\sigma_{VM}^{2}+\alphaP^{2}}\) Для таких материалов, как почва или камень, обычно используется знак «нет».

связанная формула пластичности, Iform=0.

Для связанной пластичности: I форма = 1, \(g=f\) В этом случае скорость пластической деформации равна

функция вектора нормали поверхности текучести \(f\)

. В таких материалах, как металл, обычно используется

соответствующая формула пластичности.: \(d\epsilon_{p}=d\lambda\frac{\partialf}{\partial\sigma}=d\lambda\frac{\partialg}{\partial\sigma}\) .. image:: images/mat_law76_samp_starter_r_mat_law76_iform1.png

(Рисунок 3.)

Выпуклость

флаг условия ICONV=1 используется для обеспечения стабильность материального закона за счет выпуклости поверхности текучести. Доходность Поверхность может быть гиперболической при низких значениях текучести при сдвиге при растяжении и сжатии. В в этом случае нет однозначного решения и Radioss будет обновить (увеличить) предел текучести при сдвиге, чтобы обеспечить выпуклость поверхности текучести. Следовательно, предел текучести при сдвиге может отличаться от входной кривой.

В таблицах должно быть максимум

размерность равна 2. Первая запись — это пластическая деформация, вторая запись — скорость деформации.

Пользовательские переменные USR2, USR3, USR4

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