======================= /MAT/LAW121 (PLAS_RATE) ======================= Ключевое слово в формате блока Упругопластический материал, зависящий от скорости деформации с изотропным критерием текучести фон Мизеса. Этот материальный закон доступен как для твердых тел, так и для ракушки. Формат ------ .. 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/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/mat_ID/unit_ID", "/MAT/LAW121/mat_ID/unit_ID or /MAT/PLAS_RATE/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}`", "", "", "", "", "", "", "", "" "E", "E", ":math:`\nu`", ":math:`\nu`", "Ирес", "Ивиск", "Fcut", "Fcut", "DTMIN", "DTMIN" "Fct_SIG0", "", "Xscale_SIG0", "Xscale_SIG0", "Yscale_SIG0", "Yscale_SIG0", "", "", "", "" "Fct_YOUN", "", "Xscale_YOUN", "Xscale_YOUN", "Yscale_YOUN", "Yscale_YOUN", "", "", "", "" "Fct_TANG", "", "Xscale_TANG", "Xscale_TANG", "TANG", "TANG", "", "", "", "" "Fct_FAIL", "Если не удалось", "Xscale_FAIL", "Xscale_FAIL", "Yscale_FAIL", "Yscale_FAIL", "", "", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "mat_ID", "Идентификатор материала.(Целое число, максимум 10 цифр)", "" "unit_ID", "(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)", "" "mat_title", "Название материала.(Персонаж, максимум 100 символов)", "" ":math:`\rho_{i}`", "Начальный плотность.(Реальная)", ":math:`[\frac{kg}{m^{3}}]`" "E", "Янг модуль.(Реальный)", ":math:`[Pa]`" ":math:`\nu`", "Пуассона соотношение.(Реальное)", "" "Ирес", "Метод разрешения пластичности. = 0 Установите на 2. = 1 NICE (ошибка следующего приращения) явный метод. = 2 (по умолчанию) Итерационный полунеявный метод Ньютона (Cutting самолет). (Целое число)", "" "Ивиск", "Формулировка зависимости скорости деформации. = 0 (по умолчанию) Масштабированный предел текучести. = 1 Вязко-пластическая рецептура. (Целое число)", "" "Если не удалось", "Переменная критерия отказа (используется только Fct_FAIL > 0). = 0 (по умолчанию) Эквивалентный стресс. =1 Эффективный пластический штамм. = 2 Максимальное главное напряжение или абсолютное значение минимальное главное напряжение. = 3 Максимальное главное напряжение. (Целое число)", "" "Fcut", "Частота среза для скорости деформации фильтрация (использовал только Ivisc= 0).По умолчанию = 10000 Гц (реальная)", ":math:`[Hz]`" "DTMIN", "Минимальный размер шага по времени для автоматическое удаление элемента. По умолчанию = 0,0 (Реальное)", ":math:`[s]`" "Fct_SIG0", "Предел текучести по сравнению с эффективным идентификатор функции скорости деформации. (Целое число)", "" "Xscale_SIG0", "Масштабный коэффициент для абсцисс (деформация ставка) в Fct_SIG0.По умолчанию = 1,0. (Реал)", ":math:`[\frac{1}{s}]`" "Yscale_SIG0", "Масштабный коэффициент для ординаты (напряжения) в Fct_SIG0.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" "Fct_YOUN", "Модуль Юнга в сравнении с эффективным идентификатор функции скорости деформации (доступен только Ivisc = 0).По умолчанию = 1.0 (Реал)", "" "Xscale_YOUN", "Масштабный коэффициент для абсцисс (деформация ставка) в Fct_YOUN.По умолчанию = 1,0. (Реал)", ":math:`[\frac{1}{s}]`" "Yscale_YOUN", "Масштабный коэффициент для ординаты (напряжения) в Fct_YOUN.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" "Fct_TANG", "Касательный модуль пластичности в зависимости от эффективная функция скорости деформации идентификатор.(Целое число)", "" "Xscale_TANG", "Масштабный коэффициент для абсцисс (деформация ставка) в Fct_TANG.По умолчанию = 1,0. (Реал)", ":math:`[\frac{1}{s}]`" "TANG", "Если Fct_TANG > 0 Масштабный коэффициент для ординаты (напряжения) в Fct_TANG. По умолчанию = 1,0 (реальное) Если Fct_TANG = 0 Пластическая касательная TANG. По умолчанию = 0,0 (реальное)", ":math:`[Pa]`" "Fct_FAIL", "Переменная критерия отказа по сравнению с эффективная функция скорости деформации идентификатор.(Целое число)", "" "Xscale_FAIL", "Масштабный коэффициент для абсцисс (деформация ставка) в Fct_FAIL.По умолчанию = 1,0. (Реал)", ":math:`[\frac{1}{s}]`" "Yscale_FAIL", "Масштабный коэффициент для ординаты (напряжения) в Fct_FAIL.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" Пример (Сталь) -------------- .. code-block:: #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /UNIT/1 Unit system Mg mm s #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /MAT/LAW121/1/1 Steel # RHO_I 7.85E-9 # E Nu Ires Ivisc Fcut DTMIN 194200.0 0.3 1 0 1000 0 # Fct_SIG0 Xscale_SIG0 Yscale_SIG0 5 0 0 # Fct_YOUN Xscale_YOUN Yscale_YOUN 6 0 0 # Fct_TANG Xscale_TANG TANG 7 0 0 # Fct_FAIL Ifail Xscale_FAIL Yscale_FAIL 8 1 0 0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/5 Yield stress versus effective strain rate 0.0 272.2 1000.0 572.2 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/6 Young’s modulus versus effective strain rate 0.0 192400.0 1000.0 272400.0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/7 Plastic tangent modulus versus effective strain rate 0.0 4500.0 1000.0 6500.0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/8 Failure criterion variable versus effective strain 0.0 0.15 1000.0 0.25 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #ENDDATA /END #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| Комментарии ----------- 1. Этот материальный закон учитывает изотропную и линейную упругость. 2. Пластическое поведение описывается классической формулой фон Мизеса. функция текучести пластичности, описываемая: :math:`f=\sigma_{VM}−\sigma_{Y}` С: :math:`\sigma_{VM}=\sqrt{\frac{3}{2}s:s}` Где, :math:`s` это девиатор тензор напряжений. 3. Предел текучести зависит от скорости деформации и описывается как: :math:`\sigma_{Y}=\sigma_{0}(\dot{\epsilon})+\frac{E(\dot{\epsilon})E_{t}(\dot{\epsilon})}{E(\dot{\epsilon})−E_{t}(\dot{\epsilon})}\epsilon_{p}` Где, :math:`\sigma_{0}` Зависящий от скорости деформации начальный предел текучести, определяемый формулой функция Fct_SIG0 . :math:`E(\dot{\epsilon})` Модуль Юнга, который может опционально зависеть от деформации ставка (если Fct_YOUN определяется). :math:`E_{t}(\dot{\epsilon})` Касательный модуль пластичности, который может опционально зависеть от скорость деформации (если Fct_TANG есть определено). :math:`\epsilon_{p}` Кумулятивная пластическая деформация. .. note:: Неравенство :math:`E(\dot{\epsilon})>E_{t}(\dot{\epsilon})` всегда должно выполняться. Если нет, решатель автоматически ограничить значение касательного модуля :math:`E_{t}(\dot{\epsilon})=0.99\timesE(\dot{\epsilon})` . 4. Доступны две формулы зависимости скорости деформации в зависимости от стоимость флага Ивиск . - Ивиск = 0 : масштабированный предел текучести формулировка. В этом случае скорость деформации соответствует полной эффективная скорость деформации рассчитывается как: :math:`\dot{\epsilon}=\sqrt{\frac{2}{3}\dot{\epsilon}':\dot{\epsilon}'}` Где, :math:`\dot{\epsilon}'` - девиаторная полная деформация тензор скорости. Для этой формулировки все скорости деформации зависимые переменные обновляются в начале временного шага и являются постоянными во время алгоритма отображения возврата. Все штаммы Использование этой формулировки позволяет использовать переменные, зависящие от скорости. скорость деформации в этом случае фильтруется по умолчанию как: :math:`\dot{\epsilon}_{_{f}}^{n}=\alpha\dot{\epsilon}^{n}+(1−\alpha)\dot{\epsilon}_{f}^{n−1}` Где, :math:`\alpha=2\piF_{cut}\Delta t` . Затем вы сможете адаптировать фильтрация путем указания значения частоты среза Fcut. Этот состав легче и быстрее, чем вязкопластик. один описан ниже. - Ивиск = 1 : вязкопластик формулировка. В этом случае скорость деформации соответствует Скорость пластической деформации. :math:`\dot{\epsilon}=\dot{\epsilon}_{p}` В этом случае скорость деформации изначально равна null, вычисляется и обновляется во время сопоставления возврата процедура. Все изменения переменных, зависящих от скорости деформации, приниматься во внимание, что приводит к увеличению стоимости вычислений. Однако фильтрация скорости деформации не требуется, и вы можете ускорить запустить симуляцию, используя явное отображение возврата NICE (Ирес = 1). Эта формулировка Юнга Модуль не может изменяться в зависимости от скорости деформации. 5. Вы можете выбрать алгоритм отображения возврата, чтобы найти компромисс. между точностью и стоимостью. Доступны два метода сопоставления возврата. - Ires = 1: NICE явный метод (ошибка следующего приращения). Эта процедура требует только 1 итерация и использует способность самокорректировки, чтобы по-прежнему сохранять хорошая точность. Этот метод позволяет значительно ускорить моделирование. - Ires = 2: секущая плоскость (Итерационный метод Ньютона). Эта процедура требует нескольких итерации (обычно от 3 до 5) для решения поведения материала нелинейные уравнения. Это дороже, чем явное NICE метод, но обеспечивает очень хорошую точность. Этот метод установлен по умолчанию. 6. Вы можете добавить сбой в поведение материала. Для этого функция Fct_FAIL должны быть определены для описания эволюция переменной критерия разрушения со скоростью деформации. Природа переменной критерия отказа выбирается с помощью флага Если не удалось : - Ifail = 0: критерий выполняется, когда напряжение по Мизесу достигает критического значения определено в функции - Ifail = 1: критерий выполняется, когда совокупная пластическая деформация достигает значения определено в функции - Ifail = 2: критерий выполняется, когда максимальное главное напряжение или абсолютное значение минимальное главное напряжение достигает критического значения, определенного в функция - Ifail = 3: критерий выполняется, когда максимальное главное напряжение достигает критического значение, определенное в функции 7. Автоматическое удаление элемента также можно настроить, если временной шаг элемента уменьшается до достижения критического значения, определяемого DTMIN.