=========== /MAT/LAW116 =========== Ключевое слово формата блока Описывает смешанный режим, модель материала, зависящую от скорости деформации, с повреждениями и неудача. Этот материал совместим только с твердыми шестигранными элементами ( /BRICK ) и когезионное твердое свойство ( /PROP/TYPE43 (CONNECT) ). .. note:: Не совместим ни с одним модель неудачи. Все повреждения и отказы определяются внутри материала. напрямую. Формат ------ .. 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/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/mat_ID/unit_ID", "/MAT/LAW116/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}`", "", "", "", "", "", "", "", "" ":math:`E_{I}`", ":math:`E_{I}`", ":math:`E_{II}`", ":math:`E_{II}`", "Толстый", "Толстый", "Имасс", "Идель", "Икрит", "" ":math:`GC_{I_ini}`", ":math:`GC_{I_ini}`", ":math:`GC_{I_inf}`", ":math:`GC_{I_inf}`", ":math:`\dot{\epsilon}_{G_{I}}`", ":math:`\dot{\epsilon}_{G_{I}}`", ":math:`fG_{I}`", ":math:`fG_{I}`", "", "" ":math:`GC_{II_ini}`", ":math:`GC_{II_ini}`", ":math:`GC_{II_inf}`", ":math:`GC_{II_inf}`", ":math:`\dot{\epsilon}_{G_{II}}`", ":math:`\dot{\epsilon}_{G_{II}}`", ":math:`fG_{II}`", ":math:`fG_{II}`", "", "" ":math:`\sigma_{A_I}`", ":math:`\sigma_{A_I}`", ":math:`\sigma_{B_I}`", ":math:`\sigma_{B_I}`", ":math:`\dot{\epsilon}_{I}`", ":math:`\dot{\epsilon}_{I}`", "Iorder_I", "Ifail_I", "", "" ":math:`\sigma_{A_II}`", ":math:`\sigma_{A_II}`", ":math:`\sigma_{B_II}`", ":math:`\sigma_{B_II}`", ":math:`\dot{\epsilon}_{II}`", ":math:`\dot{\epsilon}_{II}`", "Iorder_II", "Ifail_II", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "mat_ID", "Идентификатор материала.(Целое число, максимум 10 цифр)", "" "unit_ID", "(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)", "" "mat_title", "Название материала.(Персонаж, максимум 100 символов)", "" ":math:`\rho_{i}`", "Начальный плотность.(Реальная)", ":math:`[\frac{kg}{m^{3}}]`" ":math:`E_{I}`", "Модуль Юнга (жесткости) в нормальное направление на единицу длины.(Реальное)", ":math:`[\frac{Pa}{m}]`" ":math:`E_{II}`", "Модуль сдвига (жесткости) в направление касательной на единицу длины. По умолчанию = :math:`E_{II}=E_{I}` (Реал)", ":math:`[\frac{Pa}{m}]`" "Толстый", "Ссылка связная Толщина.(Реальная)", ":math:`[m]`" "Имасс", "Флаг расчета массы. = 1 (по умолчанию) Масса элемента рассчитывается с использованием плотности и среднего значения. площадь. = 2 Масса элемента рассчитывается с использованием плотности и объем. (Целое число)", "" "Идель", "Флаг отказа с указанием номера точек интеграции для удаления элемента (между 1 и 4).По умолчанию = 1 (Целое число)", "" "Икрит", "Флаг инициации выхода и повреждения. = 1 (по умолчанию) На основе квадратичного номинального напряжения. = 2 На основе максимального номинального напряжения. (Целое число)", "" ":math:`GC_{I_ini}`", "Начальное критическое выделение энергии скорость для режима I (нормальное направление).(Реальная)", ":math:`[J]`" ":math:`GC_{I_inf}`", "Верхняя граница критической энергии скорость выпуска. Указывает зависимость скорости деформации :math:`GC_{I}` .По умолчанию = 0,0 (Реал)", ":math:`[J]`" ":math:`\dot{\epsilon}_{G_{I}}`", "Эталонная (нижняя граница) скорость деформации для зависимости скорости деформации GC. Должно быть определено, если :math:`GC_{I_inf}>0` .(Реал)", ":math:`[Hz]`" ":math:`fG_{I}`", "Фактор формы для выделения энергии скорость до отказа в режиме I.(Реальная)", "" ":math:`GC_{II_ini}`", "Начальное критическое выделение энергии скорость для режима II (сдвиг).(Реальная)", ":math:`[J]`" ":math:`GC_{II_inf}`", "Верхняя граница критической энергии скорость выпуска. Указывает зависимость скорости деформации :math:`GC_{II}` .По умолчанию = 0,0 (Реал)", ":math:`[J]`" ":math:`\dot{\epsilon}_{G_{II}}`", "Эталонная (нижняя граница) скорость деформации для зависимости скорости деформации GC. Должно быть определено, если :math:`GC_{II_inf}>0` .(Реал)", ":math:`[Hz]`" ":math:`fG_{II}`", "Фактор формы для выделения энергии скорость до отказа в режиме II.(Реальная)", "" ":math:`\sigma_{A_I}`", "Статический предел текучести в режиме Я.(Реал)", ":math:`[Pa]`" ":math:`\sigma_{B_I}`", "Предел текучести, зависящий от скорости деформации срок в режиме I.(Реальный)", ":math:`[Pa]`" ":math:`\dot{\epsilon}_{I}`", "Эталонная (нижняя граница) скорость деформации значение зависимости скорости предела текучести в режиме I. Должно быть определенный :math:`\sigma_{B_I}>0` .(Реал)", ":math:`[Hz]`" "Iorder_I", "Порядок зависимости предела текучести от скорость деформации в режиме I. = 1 (по умолчанию) Линейно-логарифмическая зависимость скорости деформации. = 2 Квадратичная логарифмическая зависимость деформации ставка. (Целое число)", "" "Ifail_I", "Критерии отказа, определенные :math:`fG_{I}` : = 1 (по умолчанию) Коэффициент энергии разрушения. = 2 Соотношение смещений переломов. (Целое число)", "" ":math:`\sigma_{A_II}`", "Статический предел текучести в режиме II.(Реал)", ":math:`[Pa]`" ":math:`\sigma_{B_II}`", "Предел текучести, зависящий от скорости деформации срок в режиме II.(Реальный)", ":math:`[Pa]`" ":math:`\dot{\epsilon}_{II}`", "Эталонная (нижняя граница) скорость деформации значение зависимости предела текучести от скорости текучести в режиме II. Должно быть определено, если :math:`\sigma_{B_II}>0` .(Реал)", ":math:`[Hz]`" "Iorder_II", "Порядок зависимости предела текучести скорость деформации в режиме II. = 1 (по умолчанию) Линейно-логарифмическая зависимость скорости деформации. = 2 Квадратичная логарифмическая зависимость деформации ставка. (Целое число)", "" "Ifail_II", "Критерии отказа, определенные :math:`fG_{II}` : = 1 (по умолчанию) Коэффициент энергии разрушения. = 2 Соотношение смещений переломов. (Целое число)", "" Пример ------ .. 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----| /MAT/LAW116/3/1 MAT_COHESIVE_MIXED_MODE_ELASTOPLASTIC_RATE # RHO_I 1.2E-9 # E1 E2 Thick Imass Idel Icrit 3000 1000 0.200 2 1 0 # GC1_INI GC1_INF SRATG1 FG1 2.000 3.000 1.500 0.7 # GC2_INI GC2_INF SRATG2 FG2 9.00 0 0 0.4 # SIGA1 SIGB1 SRATE1 Iorder1 Ifail1 33.00 1.500 2.50000E-5 1 2 # SIGA2 SIGB2 SRATE2 Iorder2 Ifail2 26.00 1.300 1.00000E-5 1 2 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #enddata Комментарии ----------- 1. Упругая жесткость – это определяется с помощью: .. image:: images/mat_law116_starter_r_mat_law116_elastic_stiffness.png *(Рисунок 1.)* Где, GP Пластическая энергия в условиях постоянного стресса GC Общая энергия :math:`i` ={I,II} Режим I (нормальный) и режим II (сдвиг) Форма закона разделения тяги определяется следующим образом: - Критерии отказа определяется соотношением энергии разрушения ( I fail_i =1) :math:`0\lefG_{i}=\frac{GC_{i}(\dot{\epsilon}_{eq})}{GC_{i}(\dot{\epsilon}_{eq})}<1−\frac{\sigma(\dot{\epsilon}_{eq})^{2}}{2GC_{i}(\dot{\epsilon}_{eq})E_{i}}<1` - Критерии отказа определяется соотношением смещений переломов ( I fail_i =2) :math:`0\lefG_{i}=\frac{\delta_{i2}−\delta_{i1}}{\delta_{if}−\delta_{i1}}<1` 2. Предел текучести определяется как: - Когда I order_i =1: :math:`\sigma(\dot{\epsilon}_{eq})=\sigma_{A_i}+\sigma_{B_i}.[max(0,ln(\frac{\dot{\epsilon}_{eq}}{\dot{\epsilon}_{i}}))]` - Когда I order_i =2: :math:`\sigma(\dot{\epsilon}_{eq})=\sigma_{A_i}+\sigma_{B_i}.[max(0,ln(\frac{\dot{\epsilon}_{eq}}{\dot{\epsilon}_{i}}))]^{2}` Где, :math:`i` ={I,II}, режим I и режим II. 3. Эквивалентная скорость деформации определяется с помощью: :math:`\dot{\epsilon}_{eq}=\frac{\sqrt{\dot{\Delta }^{2}_{I}+\dot{\Delta }^{2}_{II}}}{Thick}` Где, :math:`\dot{\Delta }_{I}` Нормальная скорость. :math:`\dot{\Delta }_{II}` Скорость сдвига. 4. Зависимый от скорости перелом энергии определяются с помощью: :math:`GC_{i}(\dot{\epsilon}_{eq})=GC_{i_ini}+(GC_{i_inf}−GC_{i_ini}).exp(−\frac{\dot{\epsilon}_{G_{i}}}{\dot{\epsilon}_{eq}})` Где, :math:`i` ={I,II}, режим I и режим II. 5. Закон о пределе текучести и повреждении схема: .. image:: images/mat_law116_starter_r_mat_law116_yield_stress_damage.png *(Рисунок 2.)* 6. Для урожайности и урона на основе по квадратичному номинальному напряжению ( Икрит = 1 ): - Начало урожайности в смешанном режиме перемещение это: :math:`\delta_{m1}=\delta_{I1}\delta_{II1}.\sqrt{\frac{1+\beta^{2}}{\delta_{II1}^{2}+(\beta.\delta_{I1})^{2}}}` Где, :math:`\delta_{i1}=\frac{\sigma_{i}}{E_{i}}` :math:`i` ={I,II}, режим I и режим II. :math:`\beta=\frac{\Delta _{II}}{\Delta _{I}}` - Смешанный урон инициация – это: :math:`\delta_{m2}=\delta_{I2}\delta_{II2}.\sqrt{\frac{1+\beta^{2}}{\delta_{II2}^{2}+(\beta.\delta_{I2})^{2}}}` Где, :math:`\delta_{i2}=\delta_{i1}+\frac{fG_{i}.GC_{i}}{\sigma_{i}}` :math:`i` ={I,II}, режим I и режим II. 7. Для урожайности и урона на основе по квадратичному номинальному напряжению ( Икрит = 2 ). - Начало урожайности в смешанном режиме перемещение это: Если :math:`\beta\le\frac{\delta_{II1}}{\delta_{I1}}` : :math:`\delta_{m1}=\delta_{I1}.\sqrt{1+\beta^{2}}` Если :math:`\beta>\frac{\delta_{II1}}{\delta_{I1}}` : :math:`\delta_{m1}=\frac{\delta_{II1}}{\beta}.\sqrt{1+\beta^{2}}` Где, :math:`\beta=\frac{\Delta _{II}}{\Delta _{I}}` :math:`\Delta _{I}` Смещение соответствует режиму I (нормальному). :math:`\Delta _{II}` Смещение — режим II (сдвиг). - Смешанный урон инициация – это: Если :math:`\beta\le\frac{\delta_{II2}}{\delta_{I2}}` : :math:`\delta_{m2}=\delta_{I2}.\sqrt{1+\beta^{2}}` Если :math:`\beta>\frac{\delta_{II2}}{\delta_{I2}}` : :math:`\delta_{m2}=\frac{\delta_{II2}}{\beta}.\sqrt{1+\beta^{2}}` 8. Окончательный урон в смешанном режиме есть ( Икрит = 1 , 2 ): :math:`\delta_{mf}=\frac{\delta_{m1}.(\delta_{m1}−\delta_{m2})E_{I}GC_{II}cos^{2}\gamma+GC_{I}.(2GC_{II}+\delta_{m1}.(\delta_{m1}−\delta_{m2})E_{II}sin^{2}\gamma)}{\delta_{m1}(E_{I}GC_{II}cos^{2}\gamma+E_{II}GC_{I}sin^{2}\gamma)}` Где, :math:`\gamma=arccos(\frac{\Delta _{I}}{\Delta _{m}})` :math:`\Delta _{m}` Перемещение – смешанный режим. 9. Пластическая деформация – это определяется как: - Режим I: :math:`\Delta p_{I}=max(\Delta p_{I}(t−1),\Delta p_{I}−\delta_{m1}cos\gamma,0)` Где, :math:`(t−1)` значение предыдущего времени шаг. - Режим II: Если :math:`\sqrt{(\Delta _{II−1}−\Delta p_{II−1}(t−1))^{2}+(\Delta _{II−2}−\Delta p_{II−2}(t−1))^{2}}>\delta_{m1}` Пластическая деформация – это вычисляется для каждого направления 1 и 2 сдвига самолет. :math:`\Delta p_{II−1}=\Delta p_{II−1}(t−1)+\Delta _{II−1}−\Delta _{II−1}(t−1)` :math:`\Delta p_{II−2}=\Delta p_{II−2}(t−1)+\Delta _{II−2}−\Delta _{II−2}(t−1)` 10. Значение стресса снижается линейно от начала повреждения до окончательного повреждения ( :math:`\Delta _{m}>\delta_{m2}` ). :math:`D=max(\frac{\Delta _{m}−\delta_{m2}}{\delta_{mf}−\delta_{m2}},D(t−1),0)` Снижение напряжения рассчитывается в нормальном режиме. направление как: Если :math:`\Delta _{I}>\Delta p_{I}` , тогда :math:`\sigma_{I}=E_{I}(\Delta _{I}−\Delta p_{I})` . в противном случае, :math:`\sigma_{I}=E_{I}(1−D)(\Delta _{I}−\Delta p_{I})` . Для каждого направления 1 и 2 в плоскость сдвига. :math:`\sigma_{II−1}=E_{II}(1−D)(\Delta _{II−1}−\Delta p_{II−1})` :math:`\sigma_{II−2}=E_{II}(1−D)(\Delta _{II−2}−\Delta p_{II−2})` 11. Соединительный элемент – это удалено, когда :math:`\Delta _{m}>\delta_{mf}` .