========== /MAT/LAW77 ========== Ключевое слово в формате блока Этот закон о пеноматериалах с открытыми порами является обобщением ЗАКОН70. Это объясняет течение невязкого сжимаемого идеального газа внутри пена и ее взаимодействие со структурой пены. ALE-моделирование потока газа и лагранжево моделирование деформации пены. выполняются по той же системе элементов. Взаимодействие потока газа с конструкцией происходит посредством закона Дарси и прямого воздействия давления газа на конструкцию. Формат ------ .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "/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" ":math:`\rho_{i}`", ":math:`\rho_{i}`", "", "", "", "", "", "", "", "" "E0", "E0", ":math:`\nu`", ":math:`\nu`", "Емакс", "Емакс", ":math:`\epsilon_{max}`", ":math:`\epsilon_{max}`", "FP0", "FP0" "Fcut", "Fcut", "Фгладкий", "NL", "НуЛ", "Ифлаг", "Форма", "Форма", "Хис", "Хис" If N L ≠ 0 , каждая загрузка функция на строку .. 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_IDL", ":math:`\dot{\epsilon}_{L}`", ":math:`\dot{\epsilon}_{L}`", "FшкалаL", "FшкалаL", "", "", "", "", "" If N uL ≠ 0 , каждая выгрузка функция на строку .. 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_IDuL", ":math:`\dot{\epsilon}_{uL}`", ":math:`\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", "", "", "", "", "", "", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "mat_ID", "Идентификатор материала (целое число, максимум 10 цифр)", "" "unit_ID", "Идентификатор объекта (целое число, максимум 10 цифр)", "" "mat_title", "Название материала(Символ, максимум 100 персонажи)", "" ":math:`\rho_{i}`", "Начальная плотность (реальная)", ":math:`[\frac{kg}{m^{3}}]`" "E0", "Начальный Янг модуль.(Реальный)", ":math:`[Pa]`" ":math:`\nu`", "Коэффициент Пуассона.(Реальный)", "" "Емакс", "Максимальный модуль Юнга.(Реальный)", ":math:`[Pa]`" ":math:`\epsilon_{max}`", "Эталонное значение деформации для максимального Использование модуля Юнга. По умолчанию = 1 (Реальное)", "" "FP0", "Начальное давление пены. По умолчанию = 0 (Реал)", ":math:`[Pa]`" "Fcut", "Частота среза для скорости деформации фильтрация.По умолчанию = 1020 (Реальное)", ":math:`[Hz]`" "Фгладкий", "Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)", "" "NL", "Количество функций загрузки.По умолчанию = 0 (целое число)", "" "НуЛ", "Количество функций разгрузки.По умолчанию = 0 (целое число)", "" "Ифлаг", "Флаг управления ответом на выгрузку. = 0 Поведение материала соответствует определенным кривым для нагрузки и разгрузка. = 1 (по умолчанию) Кривые нагрузки используются как для погрузки, так и для разгрузки. Для при разгрузке девиаторное напряжение повреждается при использовании квазистатической разгрузки кривая: :math:`\sigma=(1−D)(\sigma+P)−P` Где D вычисляется с учетом квазистатическая кривая разгрузки, :math:`D=(\frac{\sigma_{unloading}}{\sigma_{quasi−static}})` текущие напряжения, рассчитанные соответственно из кривые разгрузки и квазистатические кривые. = 2 Кривые нагрузки используются как для погрузки, так и для разгрузки. Для при разгрузке тензорное напряжение снижается за счет использования квазистатической разгрузки кривая: :math:`\sigma=(1−D)\sigma` где D вычисляется с учетом квазистатического кривая разгрузки, :math:`D=(\frac{\sigma_{unloading}}{\sigma_{quasi−static}})` — текущие напряжения, рассчитанные соответственно из разгрузки и квазистатические кривые. = 3 Кривые нагрузки используются как для погрузки, так и для разгрузки. Девиаторное напряжение разгрузки снижается за счет: :math:`\sigma=(1−D)(\sigma+P)−P` :math:`D=(1−Hys)(1−(\frac{W_{cur}}{W_{max}})^{Shape})` Где, Wcur и Wmax являются текущими и максимальная энергия. Если Iflag = 3, кривые разгрузки не используются. (Целое число)", "" "Форма", "Фактор формы описывает «выпуклость». кривой разгрузки. = 1 Соответствует линейной кривой разгрузки. < 1 Выпуклая кривая разгрузки. По умолчанию = 1,0 (реальное)", "" "Хис", "Коэффициент разгрузки гистерезиса. По умолчанию = 1.0 (Реал)", "" "fct_IDL", "Функция загрузки идентификатор.(Целое число)", "" ":math:`\dot{\epsilon}_{L}`", "Скорость деформации для нагрузки функция.(Реальная)", ":math:`[\frac{1}{s}]`" "FшкалаL", "Масштабный коэффициент функции загрузки. По умолчанию = 1.0 (Реал)", ":math:`[Pa]`" "fct_IDuL", "Функция разгрузки идентификатор.(Целое число)", "" ":math:`\dot{\epsilon}_{uL}`", "Скорость деформации для разгрузки функция.(Реальная)", ":math:`[\frac{1}{s}]`" "FscaleuL", "Масштабный коэффициент функции выгрузки. По умолчанию = 1,0 (Реальный)", ":math:`[Pa]`" ":math:`\rho_{gas}`", "Плотность воздуха. По умолчанию = 0 (Реал)", ":math:`[\frac{kg}{m^{3}}]`" "P0", "Начальное давление. По умолчанию = 0 (Реал)", ":math:`[Pa]`" ":math:`\gamma`", "Гамма-константа для газа.По умолчанию (Реал)", "" ":math:`\rho_{ext}`", "Внешняя плотность газа. По умолчанию: :math:`\rho_{gas}` (Реал)", ":math:`[\frac{kg}{m^{3}}]`" "Пекст", "Внешнее давление. По умолчанию = P0 (Реал)", ":math:`[Pa]`" "Inc_gas", "Флаг обратного потока. = 0 (по умолчанию) Отсутствие обратного потока наружного воздуха в пенопласт. = 1 Допускается обратный поток. (Целое число)", "" "R", "Начальная пористость элемента объем, заполненный газом.(0 < R < 1) = 0,0 Нет газа в элементе, нет потока газа. = 1,0 Полный объем элемента используется для потока газа. (Настоящий)", "" ":math:`\alpha`", "Линейный параметр для обобщенного Дарси Закон.По умолчанию = 0 (реальный)", "" ":math:`\beta`", "Квадратичный параметр для обобщенного Закон Дарси.По умолчанию = 0 (Реальный)", ":math:`[\frac{s}{m}]`" ":math:`\tau`", "Переходный параметр для обобщенных Закон Дарси.По умолчанию = 0 (Реальный)", ":math:`[s]`" "K", "Начальная проницаемость пены модуль.По умолчанию = 0 (Реальный)", ":math:`[\frac{m^{2}}{Pa⋅s}]`" "Иклос", "Открыть/закрыть свободную поверхность твердого тела флаг блока (поверхность не связана ни с одним твердотельным элементом). = 0 (по умолчанию) Свободный отток на свободной поверхности (открытый). = 1 Нет истечения на свободной поверхности (закрыть). = 2 Отсутствие утечки газа через внешние поверхности, находящиеся в контакте. Параметр Ibag= 1 должен быть активирован в соответствующий контакт. (Целое число)", "" "fct_IDK", "Функция масштабного коэффициента проницаемости (накипь в зависимости от относительной плотности пены). = 0 (по умолчанию) Модуль проницаемости постоянный. (Целое число)", "" "fct_IDR", "Функция масштабного коэффициента пористости (шкала относительно относительной плотности пены). = 0 (по умолчанию) Пористость постоянна и fct_IDR= fct_IDK. (Целое число)", "" Пример (пористая пена) ---------------------- .. code-block:: #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. Этот материальный закон можно использовать только с твердые невырожденные гексаэлементы. Этот материал доступен только для следующих параметры в свойстве твердого тела: - Изолид = 1 (Беличко) - Ismstr = 1 (малая деформация) - Iframe = 1 (не совместное вращение) 2. Для напряжений выше последней функции нагрузки поведение экстраполируется с использованием двух последних функций нагрузки. Во избежание огромных значения напряжения, рекомендуется повторить последнюю функцию нагрузки. 3. Невязкое течение идеального газа с использованием Рассмотрен подход ALE. 4. Когда :math:`\epsilon_{p}` достигает :math:`\epsilon_{p}^{max}` , в одной точке интегрирования, девиаторное напряжение соответствующая целая точка постоянно установлена на 0; однако твердый элемент не удален. 5. Для соединения пены с газом использовать модифицированный закон Дарси (Дюпюи-Форхимера): :math:`\alphaV+\betaV^{2}+\tau\frac{\partialV}{\partialt}=−K⋅Grad(P)` Где, :math:`V` Скорость жидкости :math:`P` Давление жидкости :math:`K` Модуль проницаемости 6. Чтобы сохранить структурное напряжение, данных о деформации и расходе газа, в Engine необходимо использовать следующие параметры: файл:/STATE/BRICK/STRAIN/FULL/STATE/BRICK/STRESS/FULL/STATE/BRICK/AUX/FULL 7. Предварительное напряжение пенопласта и исходное состояние воздух внутри пенопласта можно определить на основе моделирования предварительного напряжения и ввести в модель с картами /INIBRI/STRS_F, /INIBRI/STRA_F и /INIBRI/AUX. 8. Доступны следующие пользовательские переменные для вывода данных о расходе газа в анимацию ( /ANIM/BRICK/TENS ) или история времени ( /TH/BRIC ) файл: - USR1: плотность газа - USR20: давление газа - USR21: значение R - USR22: значение K - Векторное поле скорости газа можно вывести в анимации с помощью /ANIM/VECT/GVEL.