===================== /MAT/LAW43 (HILL_TAB) ===================== Ключевое слово формата блока Этот закон описывает ортотропный материал Хилла и применимо только к элементам оболочки. Этот закон отличается от LAW32 (HILL) только тем, что ввод предела текучести (здесь он определяется пользовательской функцией). Формат ------ .. 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/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_TAB/mat_ID/unit_ID", "/MAT/LAW43/mat_ID/unit_ID or /MAT/HILL_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" ":math:`\rho_{i}`", ":math:`\rho_{i}`", "", "", "", "", "", "", "", "" "E", "E", ":math:`\nu`", ":math:`\nu`", "", "", "", "", "", "" "fct_IDE", "", "Эйнф", "Эйнф", "CE", "CE", "", "", "", "" "r00", "r00", "r45", "r45", "r90", "r90", "Мангольд", "Мангольд", "Доходность0", "" ":math:`\epsilon_{p}^{max}`", ":math:`\epsilon_{p}^{max}`", ":math:`\epsilon_{t}`", ":math:`\epsilon_{t}`", ":math:`\epsilon_{m}`", ":math:`\epsilon_{m}`", "Fcut", "Fcut", "Фгладкий", "" Каждая функция пластичности в строке .. 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_IDi", "", "Фскалей", "Фскалей", ":math:`\dot{\epsilon}_{i}`", ":math:`\dot{\epsilon}_{i}`", "", "", "", "" Определение ----------- .. 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`", "Коэффициент Пуассона.(Реальный)", "" "fct_IDE", "Идентификатор функции для масштабного коэффициента модуля Юнга, когда модуль Юнга является функцией пластической деформации. 12По умолчанию = 0: в этом случае эволюция Модуль Юнга зависит от Einf и CE. (Целое число)", "" "Эйнф", "Насыщенный модуль Юнга для инфинитива пластическая деформация.(Настоящая)", ":math:`[Pa]`" "CE", "Параметр эволюции модуля Юнга. 12(Реал)", "" "r00", "Параметр Ланкфорда 0 градусов. 3По умолчанию = 1,0 (реальное)", "" "r45", "Параметр Ланкфорда 45 градусов. По умолчанию = 1,0 (Реальный)", "" "r90", "Параметр Ланкфорда 90 градусов.По умолчанию = 1,0 (Реальный)", "" "Мангольд", "Коэффициент закалки. = 0 Закалка – полностью изотропная модель. = 1 Для закалки используется кинематическая модель Прагера-Циглера. = между 0 и 1 Упрочнение интерполируется между двумя моделями. (Настоящий)", "" "Доходность0", "Флаг стресса доходности. = 0 Средний входной предел текучести. = 1 Предел текучести в ортотропном направлении 1. (Целое число)", "" ":math:`\epsilon_{p}^{max}`", "Пластическая деформация разрушения. По умолчанию = 1,0. × 1030 (Реал)", "" ":math:`\epsilon_{t}`", "Деформация разрушения при растяжении, при которой напряжение начинает уменьшаться. По умолчанию = 1,0 × 1030 (Реальное)", "" ":math:`\epsilon_{m}`", "Максимальная деформация разрушения при растяжении, при которой напряжение в элементе установлено на ноль. По умолчанию = 2,0 × 1030. (Реал)", "" "Fcut", "Частота среза для скорости деформации фильтрация.(Реальная)", ":math:`[Hz]`" "Фгладкий", "Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)", "" "fct_IDi", "Кривые пластичности .ith идентификатор функции.(Целое число)", "" "Фскалей", "Масштабный коэффициент для i-я функция. По умолчанию установлено значение 1,0. (Реал)", "" ":math:`\dot{\epsilon}_{i}`", "Скорость деформации для i-я функция.(Реальная)", ":math:`[\frac{1}{s}]`" Пример (Металл) --------------- .. code-block:: #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/HILL_TAB/1/1 metal # RHO_I 7.8E-09 # E NU 206000 .3 #FUNCT_IDE EINF CE 0 0 0 # r00 r45 r90 C_hard Iyield0 1.73 1.34 2.24 0 0 # EPSP_max EPS_t1 EPS_m Fcut Fsmooth 0 0 0 15000 1 # func_IDi Fscale_i EPS_dot_i 5 0 0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #- 3. FUNCTIONS: #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/5 metal # X Y 0 260 .002 270 .005 280 .01 297 .02 322 .05 370 .1 422 .15 457 .2 485 .3 528 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #ENDDATA /END #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| Комментарии ----------- 1. Этот материальный закон необходимо использовать с набором свойств /PROP/TYPE9 (SH_ORTH) или /PROP/TYPE10 (SH_COMP). 2. Предел текучести определяется пользовательскую функцию, а предел текучести сравнивается с эквивалентным напряжением: :math:`\sigma_{eq}=\sqrt{A_{1}\sigma_{1}^{2}+A_{2}\sigma_{2}^{2}−A_{3}\sigma_{1}\sigma_{2}+A_{12}\sigma_{12}^{2}}` 3. Углы для параметров Ланкфорда равны определяется относительно ортотропного направления 1. :math:`R=\frac{r_{00}+2r_{45}+r_{90}}{4}H=\frac{R}{1+R}A_{1}=H(1+\frac{1}{r_{00}})A_{2}=H(1+\frac{1}{r_{90}})A_{3}=2HA_{12}=2H(r_{45}+0.5)(\frac{1}{r_{00}}+\frac{1}{r_{90}})r_{00}=\frac{A_{3}}{2A_{1}−A_{3}}r_{45}=\frac{1}{2}(\frac{A_{12}}{A_{1}+A_{2}−A_{3}}−1)r_{90}=\frac{A_{3}}{2A_{2}−A_{3}}` Параметры Ланкфорда ra – коэффициент пластической деформации в плоская и пластическая деформация по толщине :math:`\epsilon_{33}` . :math:`r_{\alpha}=\frac{d\epsilon_{\alpha+\pi/2}}{d\epsilon_{33}}` Где α — угол к ортотропному направлению 1. Это Параметры Лэнкфорда ra могут быть определяется путем простого испытания на растяжение под углом α. Более высокое значение R означает лучшую формуемость. 4. Если последняя точка первой (статическая) функция равна 0 в стрессе, значение по умолчанию :math:`\epsilon_{p}^{max}` устанавливается на соответствующее значение :math:`\epsilon_{p}` . 5. Удаление элемента: - Один раз :math:`\epsilon_{p}` (пластическая деформация) достигает :math:`\epsilon_{p}^{max}` , в одной точке интеграции элемент удаляется. - If :math:`\epsilon_{1}` достигает :math:`\epsilon_{t}` , напряжение уменьшается по следующему соотношению: :math:`\sigma=\sigma(\frac{\epsilon_{m}−\epsilon_{1}}{\epsilon_{m}−\epsilon_{t}})` - If :math:`\epsilon_{1}` (крупнейший основной штамм) достигает :math:`\epsilon_{m}` ( :math:`\epsilon_{1}>\epsilon_{m}` ), напряжение в элементе снижается до 0 (но элемент не удаляется). - Один раз :math:`\epsilon_{1}` (крупнейший основной штамм) достигает :math:`\epsilon_{f}` (максимальная деформация разрушения при растяжении), элемент удален. 6. Максимальное количество кривых, которые можно ввести 10. 7. If :math:`\dot{\epsilon}\le\dot{\epsilon}_{n}` , доходность интерполируется между :math:`f_{n}` и :math:`f_{n−1}` . 8. If :math:`\dot{\epsilon}\le\dot{\epsilon}_{1}` , функция :math:`f_{1}` используется. 9. Выше :math:`\dot{\epsilon}_{max}` , доходность экстраполируется. .. image:: images/mat_law43_hill_tab_starter_r_mat_law43_yield.png :alt: mat_law43_yield *(Рисунок 1.)* 10. Радиальный возврат недоступен. (только итерационная пластичность). 11. Если предел текучести был полученное в ортотропном направлении 1, определим Доходность0 =1; иначе Доходность0 =0. 12. Эволюция Янга модуль: - If fct_ID E > 0 , кривая определяет масштабный коэффициент для модуля Юнга эволюция с эквивалентной пластической деформацией :math:`\bar{\epsilon}_{p}` , что означает, что модуль Юнга масштабируется по функция :math:`f(\bar{\epsilon}_{p})` : :math:`E(t)=f(\bar{\epsilon}_{p})E` Начальное значение масштабного коэффициента должно быть равно 1 и оно уменьшается. - If fct_ID E = 0 , модуль Юнга рассчитывается как: :math:`E(t)=E−(E−E_{inf})(1−exp(−C_{E}\bar{\epsilon}_{p}))` Где, :math:`E` и :math:`E_{inf}` являются соответственно начальным и асимптотическим значением модуль Юнга и :math:`\bar{\epsilon}_{p}` – накопленная эквивалентная пластическая деформация. .. note:: Если fct_IDE = 0 и CE = 0, то Юнга модуль E остается постоянным. 13. Параметры F гладкий и F резать позволяют включить фильтрацию по скорости деформации. Можно установить три случая: - Если Fsmooth = 0 + Fcut = 0,0, фильтрация скорости деформации отключается. - Если Fsmooth = 1 + Fcut = 0,0, фильтрация скорости деформации использует частоту среза по умолчанию 10 кГц. - Если Fcut ≠ 0, Fsmooth автоматически устанавливается на 1, а фильтрация скорости деформации использует частоту среза. предоставленный вами.