/MAT/LAW77
- Ключевое слово в формате блока Этот закон о пеноматериалах с открытыми порами является обобщением
ЗАКОН70. Это объясняет течение невязкого сжимаемого идеального газа внутри пена и ее взаимодействие со структурой пены.
- ALE-моделирование потока газа и лагранжево моделирование деформации пены.
выполняются по той же системе элементов. Взаимодействие потока газа с конструкцией происходит посредством закона Дарси и прямого воздействия давления газа на конструкцию.
Формат
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/mat_ID/unit_ID |
/MAT/LAW77/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\) |
\(\nu\) |
Емакс |
Емакс |
\(\epsilon_{max}\) |
\(\epsilon_{max}\) |
FP0 |
FP0 |
Fcut |
Fcut |
Фгладкий |
NL |
НуЛ |
Ифлаг |
Форма |
Форма |
Хис |
Хис |
If N L ≠ 0 , каждая загрузка
функция на строку
fct_IDL |
\(\dot{\epsilon}_{L}\) |
\(\dot{\epsilon}_{L}\) |
FшкалаL |
FшкалаL |
If N uL ≠ 0 , каждая выгрузка
функция на строку
fct_IDuL |
\(\dot{\epsilon}_{uL}\) |
\(\dot{\epsilon}_{uL}\) |
FscaleuL |
FscaleuL |
Карты на газ и внешний газ .. 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:`\rho_{gas}`", ":math:`\rho_{gas}`", "P0", "P0", ":math:`\gamma`", ":math:`\gamma`", "", "", "R", "R"
":math:`\rho_{ext}`", ":math:`\rho_{ext}`", "Пекст", "Пекст", "Иклос", "Inc_gas", "", "", "", ""
":math:`\alpha`", ":math:`\alpha`", ":math:`\beta`", ":math:`\beta`", ":math:`\tau`", ":math:`\tau`", "K", "K", "", ""
"fct_IDK", "fct_IDR", "", "", "", "", "", "", "", ""
Определение
Поле |
Содержание |
Пример единицы СИ |
|---|---|---|
mat_ID |
Идентификатор материала (целое число, максимум 10 цифр) |
|
unit_ID |
Идентификатор объекта (целое число, максимум 10 цифр) |
|
mat_title |
Название материала(Символ, максимум 100 персонажи) |
|
\(\rho_{i}\) |
Начальная плотность (реальная) |
\([\frac{kg}{m^{3}}]\) |
E0 |
Начальный Янг модуль.(Реальный) |
\([Pa]\) |
\(\nu\) |
Коэффициент Пуассона.(Реальный) |
|
Емакс |
Максимальный модуль Юнга.(Реальный) |
\([Pa]\) |
\(\epsilon_{max}\) |
Эталонное значение деформации для максимального Использование модуля Юнга. По умолчанию = 1 (Реальное) |
|
FP0 |
Начальное давление пены. По умолчанию = 0 (Реал) |
\([Pa]\) |
Fcut |
Частота среза для скорости деформации фильтрация.По умолчанию = 1020 (Реальное) |
\([Hz]\) |
Фгладкий |
Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число) |
|
NL |
Количество функций загрузки.По умолчанию = 0 (целое число) |
|
НуЛ |
Количество функций разгрузки.По умолчанию = 0 (целое число) |
|
Ифлаг |
Флаг управления ответом на выгрузку. = 0 Поведение материала соответствует определенным кривым для нагрузки и разгрузка. = 1 (по умолчанию) Кривые нагрузки используются как для погрузки, так и для разгрузки. Для при разгрузке девиаторное напряжение повреждается при использовании квазистатической разгрузки кривая: \(\sigma=(1−D)(\sigma+P)−P\) Где D вычисляется с учетом квазистатическая кривая разгрузки, \(D=(\frac{\sigma_{unloading}}{\sigma_{quasi−static}})\) текущие напряжения, рассчитанные соответственно из кривые разгрузки и квазистатические кривые. = 2 Кривые нагрузки используются как для погрузки, так и для разгрузки. Для при разгрузке тензорное напряжение снижается за счет использования квазистатической разгрузки кривая: \(\sigma=(1−D)\sigma\) где D вычисляется с учетом квазистатического кривая разгрузки, \(D=(\frac{\sigma_{unloading}}{\sigma_{quasi−static}})\) — текущие напряжения, рассчитанные соответственно из разгрузки и квазистатические кривые. = 3 Кривые нагрузки используются как для погрузки, так и для разгрузки. Девиаторное напряжение разгрузки снижается за счет: \(\sigma=(1−D)(\sigma+P)−P\) \(D=(1−Hys)(1−(\frac{W_{cur}}{W_{max}})^{Shape})\) Где, Wcur и Wmax являются текущими и максимальная энергия. Если Iflag = 3, кривые разгрузки не используются. (Целое число) |
|
Форма |
Фактор формы описывает «выпуклость». кривой разгрузки. = 1 Соответствует линейной кривой разгрузки. < 1 Выпуклая кривая разгрузки. По умолчанию = 1,0 (реальное) |
|
Хис |
Коэффициент разгрузки гистерезиса. По умолчанию = 1.0 (Реал) |
|
fct_IDL |
Функция загрузки идентификатор.(Целое число) |
|
\(\dot{\epsilon}_{L}\) |
Скорость деформации для нагрузки функция.(Реальная) |
\([\frac{1}{s}]\) |
FшкалаL |
Масштабный коэффициент функции загрузки. По умолчанию = 1.0 (Реал) |
\([Pa]\) |
fct_IDuL |
Функция разгрузки идентификатор.(Целое число) |
|
\(\dot{\epsilon}_{uL}\) |
Скорость деформации для разгрузки функция.(Реальная) |
\([\frac{1}{s}]\) |
FscaleuL |
Масштабный коэффициент функции выгрузки. По умолчанию = 1,0 (Реальный) |
\([Pa]\) |
\(\rho_{gas}\) |
Плотность воздуха. По умолчанию = 0 (Реал) |
\([\frac{kg}{m^{3}}]\) |
P0 |
Начальное давление. По умолчанию = 0 (Реал) |
\([Pa]\) |
\(\gamma\) |
Гамма-константа для газа.По умолчанию (Реал) |
|
\(\rho_{ext}\) |
Внешняя плотность газа. По умолчанию: \(\rho_{gas}\) (Реал) |
\([\frac{kg}{m^{3}}]\) |
Пекст |
Внешнее давление. По умолчанию = P0 (Реал) |
\([Pa]\) |
Inc_gas |
Флаг обратного потока. = 0 (по умолчанию) Отсутствие обратного потока наружного воздуха в пенопласт. = 1 Допускается обратный поток. (Целое число) |
|
R |
Начальная пористость элемента объем, заполненный газом.(0 < R < 1) = 0,0 Нет газа в элементе, нет потока газа. = 1,0 Полный объем элемента используется для потока газа. (Настоящий) |
|
\(\alpha\) |
Линейный параметр для обобщенного Дарси Закон.По умолчанию = 0 (реальный) |
|
\(\beta\) |
Квадратичный параметр для обобщенного Закон Дарси.По умолчанию = 0 (Реальный) |
\([\frac{s}{m}]\) |
\(\tau\) |
Переходный параметр для обобщенных Закон Дарси.По умолчанию = 0 (Реальный) |
\([s]\) |
K |
Начальная проницаемость пены модуль.По умолчанию = 0 (Реальный) |
\([\frac{m^{2}}{Pa⋅s}]\) |
Иклос |
Открыть/закрыть свободную поверхность твердого тела флаг блока (поверхность не связана ни с одним твердотельным элементом). = 0 (по умолчанию) Свободный отток на свободной поверхности (открытый). = 1 Нет истечения на свободной поверхности (закрыть). = 2 Отсутствие утечки газа через внешние поверхности, находящиеся в контакте. Параметр Ibag= 1 должен быть активирован в соответствующий контакт. (Целое число) |
|
fct_IDK |
Функция масштабного коэффициента проницаемости (накипь в зависимости от относительной плотности пены). = 0 (по умолчанию) Модуль проницаемости постоянный. (Целое число) |
|
fct_IDR |
Функция масштабного коэффициента пористости (шкала относительно относительной плотности пены). = 0 (по умолчанию) Пористость постоянна и fct_IDR= fct_IDK. (Целое число) |
Пример (пористая пена)
#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/LAW77/1/1
Open cell foam
# RHO_I
4.5E-8
# E0 NU EMAX EPS_max FP0
.1 0 5 .99
# FCUT FSMOOTH NLOAD NUNLOAD IFLAG SHAPE HYS
.1 1 1 0 3 2 1E-20
# F_ID_ID SLOAD FSCALELOAD
1 0 .001
# RHO_AIR P0 GAMMA R
1.2E-9 1.0E-4 1.4 1.0
# RHO_EXT P_EXT ICLOSE INC_GAS
1.2E-9 1.0E-4 2 0
# ALPHA BETA T K
1 5 0 1e07
# F_ID_K F_ID_R
2 3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1
LoadCurve
# X Y
-.8 -.11
-.7 -.10
-.4 -.05
-.2 -.02
0 0
.2 .004
.4 .006
.6 .01
.7 .020
.8 .050
.99 134
/FUNCT/2
Funct_2
# X Y
0.0 0.2
0.555 0.2
0.909 1.0
1.0 1.0
/FUNCT/3
Funct_3
# X Y
0.0 0.25
0.555 0.25
0.909 0.55
1.0 0.55
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
Комментарии
Этот материальный закон можно использовать только с
твердые невырожденные гексаэлементы. Этот материал доступен только для следующих
параметры в свойстве твердого тела:
Изолид = 1 (Беличко)
Ismstr = 1 (малая деформация)
Iframe = 1 (не
совместное вращение)
Для напряжений выше последней функции нагрузки
поведение экстраполируется с использованием двух последних функций нагрузки. Во избежание огромных значения напряжения, рекомендуется повторить последнюю функцию нагрузки.
Невязкое течение идеального газа с использованием
Рассмотрен подход ALE.
Когда
\(\epsilon_{p}\) достигает \(\epsilon_{p}^{max}\) , в одной точке интегрирования, девиаторное напряжение
соответствующая целая точка постоянно установлена на 0; однако твердый элемент не
удален.
Для соединения пены с газом
использовать модифицированный закон Дарси (Дюпюи-Форхимера):
\(\alphaV+\betaV^{2}+\tau\frac{\partialV}{\partialt}=−K⋅Grad(P)\) Где, \(V\) Скорость жидкости \(P\) Давление жидкости \(K\) Модуль проницаемости
Чтобы сохранить структурное напряжение,
данных о деформации и расходе газа, в Engine необходимо использовать следующие параметры: файл:/STATE/BRICK/STRAIN/FULL/STATE/BRICK/STRESS/FULL/STATE/BRICK/AUX/FULL
Предварительное напряжение пенопласта и исходное состояние
воздух внутри пенопласта можно определить на основе моделирования предварительного напряжения и ввести в модель с картами /INIBRI/STRS_F, /INIBRI/STRA_F и /INIBRI/AUX.
Доступны следующие пользовательские переменные
для вывода данных о расходе газа в анимацию (
/ANIM/BRICK/TENS ) или история времени ( /TH/BRIC ) файл: - USR1: плотность газа - USR20: давление газа - USR21: значение R - USR22: значение K - Векторное поле скорости газа можно вывести в анимации с помощью
/ANIM/VECT/GVEL.