/MAT/LAW38 (VISC_TAB)

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

Формат


/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID	/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_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}\)
E0	E0	\(\nu_{t}\)	\(\nu_{t}\)	\(\nu_{c}\)	\(\nu_{c}\)	\(R_{\nu}\)	\(R_{\nu}\)	Ифлаг	Итота
\(\beta\)	\(\beta\)	H	H	RD	RD	KR	KD	\(\theta\)	\(\theta\)
Каир	fct_IDp	FscaleP	FscaleP
P0	P0	RP	RP	Pмакс.	Pмакс.	\(\Phi\)	\(\Phi\)
fct_IDul		Fscaleвыгрузить	Fscaleвыгрузить	\(\dot{\epsilon}_{unload}\)	\(\dot{\epsilon}_{unload}\)	a	a	b	b
Функция		Отрезать	Отрезать	Инста
Эфинал	Эфинал	\(\epsilon_{final}\)	\(\epsilon_{final}\)	\(\lambda\)	\(\lambda\)	Виск	Виск	Тол	Тол
Fмасштаб1	Fмасштаб1	Fмасштаб2	Fмасштаб2	Fшкала3	Fшкала3	Fшкала4	Fшкала4	Fшкала5	Fшкала5
\(\dot{\epsilon}_{1}\)	\(\dot{\epsilon}_{1}\)	\(\dot{\epsilon}_{2}\)	\(\dot{\epsilon}_{2}\)	\(\dot{\epsilon}_{3}\)	\(\dot{\epsilon}_{3}\)	\(\dot{\epsilon}_{4}\)	\(\dot{\epsilon}_{4}\)	\(\dot{\epsilon}_{5}\)	\(\dot{\epsilon}_{5}\)
fct_ID1L	fct_ID2L	fct_ID3L	fct_ID4L	fct_ID5L
fct_ID1ul	fct_ID2ul	fct_ID3ul	fct_ID4ul	fct_ID5ul

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

Поле	Содержание	Пример единицы СИ
mat_ID	Идентификатор материала.(Целое число, максимум 10 цифр)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
mat_title	Название материала.(Символ, максимум 100 символов)
\(\rho_{i}\)	Начальная плотность.(Реальная)	\([\frac{kg}{m^{3}}]\)
E0	Минимальный модуль растяжения, используемый для интерфейс и расчет временного шага. (Реальный)	\([Pa]\)
\(\nu_{t}\)	Максимальный коэффициент Пуассона в напряжение.По умолчанию = 10-30 (Реальное)
\(\nu_{c}\)	Максимальный коэффициент Пуассона в сжатие.(Реальное)
\(R_{\nu}\)	Экспонента коэффициента Пуассона расчет.(Реальный)
Ифлаг	Флаг типа формулировки анализа. 4 = 0 Вязкоупругость рассчитывается в каждом направлении главного напряжения. = 1 Поведение линейное как при растяжении, так и при сжатии. (Целое число)
Итота	Флаг инкрементной формулировки. = 0 (по умолчанию) Поведение при напряжении линейное. =1 Поведение при растяжении считывается по кривым напряжения. (Целое число)
\(\beta\)	Скорость релаксации для разгрузки. По умолчанию. = 10-30 (Реал)
H	Коэффициент гистерезиса для выгрузка.По умолчанию = 1,0 (Реальный)
RD	Коэффициент демпфирования скорости деформации. По умолчанию. = 0,5 (Реальный)
KR	Флаг модели восстановления для выгрузки (петля гистерезиса). = 0 (по умолчанию) Отсутствие восстановления напряжения при разгрузке (кривая разгрузки = кривая нагрузки) = 1 Восстановление напряжения при разгрузке – это: \(\sigma=\sigma⋅H⋅min(1,1−e^{−\beta\epsilon(t)})\) = 2 Восстановление напряжения при разгрузке – это: \(\sigma=\sigma⋅{1−H⋅[1−(\frac{E^{int}}{E_{max}^{int}})^{\beta}]}\) Где, \(E^{int}\) и \(E_{max}^{int}\) — текущая внутренняя энергия и максимальная внутренняя энергии соответственно. 6 (целое число)
KD	Флаг модели распада, тип гистерезиса. = 0 (по умолчанию) Распад активен во время загрузки и разгрузки. = 1 Распад активен только во время загрузки. = 2 Распад активен во время разгрузки. (Целое число)
\(\theta\)	Коэффициент интегрирования для мгновенное обновление модуля. По умолчанию = 0,67 (Реальное)
Каир	Флаг расчета содержания воздуха. 7 = 0 (по умолчанию) Нет замкнутого воздуха = 1 Вычисление содержания замкнутого воздуха активно = 2 Считайте гидростатическую кривую (идентификатор функции определяется fct_IDp). Разница между учитываются чистое сжатие и гидростатика. (Целое число)
fct_IDp	Идентификатор кривой давления (давление относительно относительного объема).(Целое число)
FscaleP	Шкала кривой давления фактор.(Реальный)	\([Pa]\)
P0	Атмосферный давление.(Реальное)	\([Pa]\)
RP	Скорость релаксации давления. По умолчанию = 10-30 (Реал)
Pмакс.	Максимальное давление воздуха. По умолчанию = 1030 (Реал)	\([Pa]\)
\(\Phi\)	Пористость (плотность пены/плотность полимер).(Настоящий)
fct_IDul	Идентификатор функции выгрузки. > 0 Когда скорость деформации разгрузки равна статической, для разгрузки будет использоваться только функция fct_IDul. (Целое число)
Fscaleвыгрузить	Масштаб функции разгрузки коэффициент.По умолчанию = 1,0 (Реальный)	\([Pa]\)
\(\dot{\epsilon}_{unload}\)	Скорость деформации разгрузки (должна быть больше чем \(\dot{\epsilon}_{1}\) ).(Реал)	\([\frac{1}{s}]\)
a	Экспонента стресса интерполяция.По умолчанию = 1,0 (Реальное)
b	Экспонента стресса интерполяция.По умолчанию = 1,0 (Реальное)
Функция	Количество функций, определяющих скорость зависимость (пять или меньше).(Целое число)
Отрезать	Напряжение отсечки напряжения.Элемент удаляется, когда одна точка интеграции элемента превышает напряжение отсечки напряжения значение.По умолчанию = 1030 (Реальное)	\([Pa]\)
Инста	Флаг контроля нестабильности материала. = 0 (по умолчанию) Отсутствие контроля нестабильности материала. = 1 Контроль нестабильности материала. (Целое число)
Эфинал	Максимальный модуль растяжения. По умолчанию = E0 (Реал)	\([Pa]\)
\(\epsilon_{final}\)	Абсолютное значение конечной деформации модуль.По умолчанию = 1,0 (Реальный)
\(\lambda\)	Интерполяция модуля коэффициент.По умолчанию = 1,0 (Реальный)
Виск	Максимальная вязкость. 10По умолчанию = 1030 (Реал)	\([Pa⋅s]\)
Тол	Допуск по главному направлению update.Default = 1,0 (реальный)
Фскалей	Масштабный коэффициент для кривой я.(Реал)	\([Pa]\)
\(\dot{\epsilon}_{i}\)	Инженерная скорость деформации для кривой я.(Реал)	\([\frac{1}{s}]\)
fct_IDiL	Загрузка идентификатора функции для кривой я.(Целое число)
fct_IDiul	Идентификатор функции разгрузки для кривой я.(Целое число)

Пример (пена)

#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/VISC_TAB/1/1

Foam

#              RHO_I

               2E-10

#                E_0                nu_t                nu_c                 R_V     Iflag     Itota

                 200                   0                   0                   0         0         0

#               Beta                   H                 R_D       K_R       K_D                Teta

                   0                   0                   0         0         0                   0

#    K_air  fct_ID_p            Fscale_P

         0         0                   1

#                 P0                  Rp                Pmax                 Phi

                   0                   0                   0                   0

#funID_unl                 Fscale_unload        Eps_._unload                   a                   b

         0                             0                   0                   0                   0

#  N_funct                       CUT_off   I_insta

         1                             0         0

#            E_final           Eps_final              Lambda                Visc                 Tol

                   0                   0                   0                   0                   0

#      Fscale_i

                   1

#      Eps_._i

                   0

#    func_ID_iload

         4

#    func_ID_iunload

         0

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

#-  3. FUNCTIONS:

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

/FUNCT/4

function_4

#                  X                   Y

                  -1                -200

                   1                 200

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

#ENDDATA

/END

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

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

Инженерный стресс по сравнению с

инженерная деформация может быть введена как функция для различных скоростей деформации. Стресс и Деформация положительна при сжатии и отрицательна при растяжении. По умолчанию (Итота=0), поведение растяжения линейно-упругое. используя модуль Юнга E0. Если Itota=1, инженерное поведение при напряжении и деформации следует вводить с помощью функций, fct_IDiL, с кривой напряжения-деформации определяется как при сжатии, так и при растяжении

Когда функции напряжения-деформации

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

функции:

\(\sigma=f(\epsilon,\dot{\epsilon})\) для данного \(\dot{\epsilon}\) , прочитайте два значения функции по адресу \(\epsilon\) для двух непосредственно более низких и более высоких скоростей деформации. .. image:: images/mat_law38_visc_tab_starter_r_starter_mat_plas_tab2_zoom70.png

alt

starter_mat_plas_tab2

(Рисунок 1. Две кривые скорости деформации (можно ввести до пяти))

с \(\sigma=\sigma_{2}+(\sigma_{1}−\sigma_{2})[1−(\frac{\dot{\epsilon}−\dot{\epsilon}_{1}}{\dot{\epsilon}_{2}−\dot{\epsilon}_{1}})^{a}]^{b}\) Параметры

\(a\) и \(b\) определить форму интерполяционной функции внутри каждого интервал. Если \(a\) = \(b\) = 1, интерполяция линейный.

«Связанный» набор принципалов

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

\(\nu_{ij}=\nu_{c}+(\nu_{t}−\nu_{c})(1−exp(−R_{v}|\epsilon_{ij}|))\)

в напряжении ( \(\epsilon_{ij}\ge0\)

)

\(\nu_{ij}=\nu_{c}\)

в сжатии.: Где, \(\epsilon_{ij}=\frac{(\epsilon_{i}+\epsilon_{j})}{2}\) \(\epsilon_{ij}\ge0\)

Тип формулировки анализа

Ифлаг . Iflag=0: соответствует табличному материалу из вязкоупругого пенопласта.

(вязкоупругость рассчитывается в каждом направлении главного напряжения).
Еслиlag=1: поведение будет линейным как при растяжении, так и при сжатии.

следуя за отношениями Крюка.
Для сжатия модуль Юнга

\(E_{0}\) и коэффициент Пуассона \(\nu_{c}\) используются.

В напряжении мгновенный Янг

коэффициент модуля \(E_{t}\) используется.

Мгновенный модуль Юнга равен

обновлено с использованием:
\(E_{t}=E_{final}+(E_{0}−E_{final})[1−e^{−\lambda(V_{R}−1+\epsilon_{final})}]\) с \(E_{0}<E<E_{f}_{inal}\) Где, \(E_{0}\) Минимальный модуль напряжения \(E_{final}\) Максимальный модуль растяжения \(V_{R}\) Относительный объем, рассчитанный в Радиосс \(\epsilon_{final}\) Абсолютное значение деформации, соответствующее максимальному сжатию

модуль.

Мгновенный модуль используется только для растяжения.

Для стабильности,

\(\dot{\epsilon}\) фильтруется с использованием: \(\dot{\epsilon}_{filt}^{n}=\dot{\epsilon}_{filt}^{n−1}+R_{D}(\dot{\epsilon}^{n}−\dot{\epsilon}_{filt}^{n−1})\)

Гистерезис применяется в линейном

случай натяжения. Если KR=1, Гистерезис применяется только при сжатии. Если KR=2, гистерезис применяется как при сжатии, так и при растяжении.

Для давления воздуха

\(P_{air}\) (когда K воздух = 1 ) Если fct_IDp≠0: \(P_{air}=Fscale_{p}.f(\frac{V}{V_{0}})\) Где,

\(f\) относится к номеру функции fct_IDp.

Если fct_IDp=0: \(P_{air}=P_{0}\frac{(1−\frac{V}{V_{0}})}{(\frac{V}{V_{0}}−\Phi)}\) Релаксация применяется как: \(P_{air}=min(P_{air},P_{max})exp(−R_{p}t)\) Где,

\(R_{p}\) - скорость релаксации давления и \(t\) это время.

Во время разгрузки без

кривая разгрузки определена

fct_ID июль = fct_ID ul = 0 , \(\sigma\) рассчитывается по первой кривой нагрузки, fct_ID 1L . .. image:: images/mat_law38_visc_tab_starter_r_starter_mat_visc_tab_zoom76.png

alt

starter_mat_visc_tab

(Рисунок 2.)

Если кривая разгрузки определена,

\(\sigma\) интерполируется между первой кривой нагрузки fct_ID1L и заданная выгрузка кривая fct_IDul или fct_IDiul. В этом случае fct_ID1L должен соответствовать квазистатическое состояние.

Функции разгрузки

fct_IDiul (строка 12) используются только в том случае, если кривая разгрузки fct_IDul нет определен.

If

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

временной шаг:

\(\sigma\le\sigma_{1}+Visc(\dot{\epsilon}−\dot{\epsilon}_{1})\)

поведение не зависит от скорости деформации, когда

\(\dot{\epsilon}\le\dot{\epsilon}_{1}\) .

/VISC/PRONY can be used with this material law to include

viscous effects.