========================== /MAT/LAW112 (PAPER or XIA) ========================== Ключевое слово формата блока Закон картона моделирует ортотропный и асимметричный упругопластический материал из предложено Ся, 2002 г. Основной принцип – полностью разъединить поведение в плоскости бумажного листа. и поведение вне плоскости. Предел текучести определяется для каждого направления нагрузка, растяжение и сжатие. Формат ------ .. 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/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/mat_ID/unit_ID", "/MAT/LAW112/mat_ID/unit_ID or /MAT/PAPER/mat_ID/unit_ID or /MAT/XIA/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}`", "", "", "", "", "", "", "", "" "E1", "E1", "E2", "E2", "E3", "E3", "Ирес", "Итаб", "Исглад", "" ":math:`\nu_{21}`", ":math:`\nu_{21}`", "G12", "G12", "G23", "G23", "G13", "G13", "", "" "K", "K", "E3C", "E3C", "CC", "CC", "", "", "", "" ":math:`\nu_{1p}`", ":math:`\nu_{1p}`", ":math:`\nu_{2p}`", ":math:`\nu_{2p}`", ":math:`\nu_{4p}`", ":math:`\nu_{4p}`", ":math:`\nu_{5p}`", ":math:`\nu_{5p}`", "", "" If Итаб = 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 "S01", "S01", "A01", "A01", "B01", "B01", "C01", "C01", "", "" "S02", "S02", "A02", "A02", "B02", "B02", "C02", "C02", "", "" "S03", "S03", "A03", "A03", "B03", "B03", "C03", "C03", "", "" "S04", "S04", "A04", "A04", "B04", "B04", "C04", "C04", "", "" "S05", "S05", "A05", "A05", "B05", "B05", "C05", "C05", "", "" "ASIG", "ASIG", "BSIG", "BSIG", "CSIG", "CSIG", "", "", "", "" "TAU0", "TAU0", "ATAU", "ATAU", "BTAU", "BTAU", "", "", "", "" If Итаб = 1, вставить табличные значения пределов текучести .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "TAB_YLD1", "TAB_YLD1", "MAT_Xscale1", "MAT_Xscale1", "MAT_Yscale1", "MAT_Yscale1", "", "", "", "" "TAB_YLD2", "TAB_YLD2", "MAT_Xscale2", "MAT_Xscale2", "MAT_Yscale2", "MAT_Yscale2", "", "", "", "" "TAB_YLD3", "TAB_YLD3", "MAT_Xscale3", "MAT_Xscale3", "MAT_Yscale3", "MAT_Yscale3", "", "", "", "" "TAB_YLD4", "TAB_YLD4", "MAT_Xscale4", "MAT_Xscale4", "MAT_Yscale4", "MAT_Yscale4", "", "", "", "" "TAB_YLD5", "TAB_YLD5", "MAT_Xscale5", "MAT_Xscale5", "MAT_Yscale5", "MAT_Yscale5", "", "", "", "" "TAB_YLDC", "TAB_YLDC", "MAT_XscaleC", "MAT_XscaleC", "MAT_YscaleC", "MAT_YscaleC", "", "", "", "" "TAB_YLDS", "TAB_YLDS", "MAT_XscaleS", "MAT_XscaleS", "MAT_YscaleS", "MAT_YscaleS", "", "", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "mat_ID", "Идентификатор материала.(Целое число, максимум 10 цифр)", "" "unit_ID", "(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)", "" "mat_title", "Название материала.(Персонаж, максимум 100 символов)", "" ":math:`\rho_{i}`", "Начальный плотность.(Реальная)", ":math:`[\frac{kg}{m^{3}}]`" "Ei", "Модуль Юнга в i-е ортотропное направление.(Реальное)", ":math:`[Pa]`" ":math:`\nu_{ij}`", "Коэффициент Пуассона, связанный с i-й и j-й ортотропный направление.(Реальное)", "" ":math:`G_{ij}`", "Модуль сдвига, связанный с i-й и j-й ортотропный направление.(Реальное)", ":math:`[Pa]`" "Ирес", "Метод разрешения пластичности. = 0 Установить на 2 = 1 NICE (ошибка следующего приращения) явный метод. = 2 (по умолчанию) Итерационный метод Ньютона – секущая плоскость. (Целое число)", "" "Итаб", "Расчет пределов текучести. = 0 Непрерывные напряжения текучести. = 1 Табличные значения пределов текучести. (Целое число)", "" "Исглад", "Тип интерполяции (в случае табличная функция доходности). = 1 (по умолчанию) Линейная интерполяция. = 2 Логарифмическая интерполяция по основанию 10. =3 Логарифмическая база интерполяции n. (Целое число)", "" "K", "Поверхность текучести в плоскости показатель.По умолчанию = 1,0 (Реальный)", "" "E3C", "Первое эластичное сжатие параметр.По умолчанию = E3 (Реальный)", ":math:`[Pa]`" "CC", "Второе эластичное сжатие параметр.По умолчанию = 1,0 (Реальное)", "" ":math:`\nu_{1p}`", "Растяжимый пластик Коэффициент Пуассона в направление 1.(Реальное)", "" ":math:`\nu_{2p}`", "Растяжимый пластик Коэффициент Пуассона в направление 2.(Реальное)", "" ":math:`\nu_{4p}`", "Коэффициент Пуассона сжимающей пластмассы в направлении 1.(Реальное)", "" ":math:`\nu_{5p}`", "Коэффициент Пуассона сжимающей пластмассы в направлении 2.(Реальное)", "" "S0i", "Начальный предел текучести в i-е направление загрузки. Каждое направление привязанные к данному направлению нагрузки в следующем порядке: я = 1 Натяжение в ортотропном направлении 1. я = 2 Натяжение в ортотропном направлении 2. я = 3 Плоский сдвиг. я = 4 Сжатие в ортотропном направлении 1. я = 5 Сжатие в ортотропном направлении 2. По умолчанию = 1.0e20 (реальное)", ":math:`[Pa]`" "A0i", "Первый параметр закалки в i-е направление нагрузки.(Реальное)", ":math:`[Pa]`" "B0i", "Второй параметр закалки в i-е направление нагрузки.(Реальное)", "" "C0i", "Третий параметр закалки в i-е направление нагрузки.(Реальное)", ":math:`[Pa]`" "ASIG", "Начальный предел текучести вне плоскости в сжатии. По умолчанию = 1.0e20 (Реальное)", ":math:`[Pa]`" "BSIG", "Первая внеплоскостная закалка параметр сжатия.(Реальный)", ":math:`[Pa]`" "CSIG", "Вторая внеплоскостная закалка параметр сжатия.(Реальный)", "" "TAU0", "Начальная текучесть при поперечном сдвиге стресс.По умолчанию = 1.0e20 (Реальный)", ":math:`[Pa]`" "ATAU", "Первая поперечная закалка при сдвиге параметр.(Реальный)", ":math:`[Pa]`" "BTAU", "Второе упрочнение поперечным сдвигом параметр.(Реальный)", "" "TAB_YLDi", "Табличный предел текучести – пластик штамм - идентификатор функции скорости деформации в i-м направление загрузки.(Целое число)", "" "MAT_Xscalei", "Масштабный коэффициент X табличного текучесть – пластическая деформация — функция скорости деформации в i-е направление загрузки. По умолчанию = 1,0 (Реал)", ":math:`[Hz]`" "MAT_Yscalei", "Масштабный коэффициент Y табличного текучесть – пластическая деформация — функция скорости деформации в i-е направление загрузки. По умолчанию = 1,0 (Реал)", ":math:`[Pa]`" Пример (бумага) --------------- .. code-block:: #RADIOSS STARTER /UNIT/1 unit for mat Mg mm s #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /MAT/LAW112/1/1 Xia # RHO_I 7.83E-10 # E1 E2 E3 Ires Itab Ismooth 4193 1554 1554 2 0 0 # nu21 G12 G23 G13 0.1011 988 76 76 # K E3C CC 2.0 47.2 24.46 # nu1p nu2p nu4p nu5p 0.555 0.1537 0.18 0.145 # S01 A01 B01 C01 12.0 19.0 260.0 800.0 # S02 A02 B02 C02 6.5 40.0 160.0 250.0 # S03 A03 B03 C03 6.0 11.0 100.0 125.0 # S04 A04 B04 C04 7.3 6.0 160.0 300.0 # S05 A05 B05 C05 6.3 9.0 310.0 225.0 # ASIG BSIG CSIG 16.55 16.55 3.16 # TAU0 ATAU BTAU 2.1 9.0 2.0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #ENDDATA Пример (в таблице) ------------------ .. code-block:: #RADIOSS STARTER /UNIT/1 unit for mat Mg mm s #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /MAT/LAW112/1/1 Xia_tab # RHO_I 7.83E-10 # E1 E2 E3 Ires Itab Ismooth 4193 1554 1554 1 1 1 # nu21 G12 G23 G13 0.1011 988 76 76 # K E3C CC 2.0 47.2 24.46 # nu1p nu2p nu4p nu5p 0.555 0.1537 0.18 0.145 # TAB_YLD1 MAT_Xscale1 MAT_Yscale1 25 1.0 1.0 # TAB_YLD2 MAT_Xscale2 MAT_Yscale2 25 1.0 0.35 # TAB_YLD3 MAT_Xscale3 MAT_Yscale3 25 1.0 0.75 # TAB_YLD4 MAT_Xscale4 MAT_Yscale4 25 1.0 0.6341 # TAB_YLD5 MAT_Xscale5 MAT_Yscale5 25 1.0 0.5 # TAB_YLDC MAT_XscaleC MAT_YscaleC 25 1.0 0.5 # TAB_YLDS MAT_XscaleS MAT_YscaleS 25 1.0 0.5 /FUNCT/46 ecoulement2 # plastic strain stress 0.0 12.00 0.012 32.979020979021 0.025 50.4615384615385 0.05 74.5 0.075 90.9473684210526 0.1 102.909090909091 0.125 112.00 0.15 119.142857142857 0.175 124.903225806452 0.2 129.647058823529 0.25 137.00 0.3 142.434782608696 0.4 149.931034482759 0.5 154.857142857143 1.0 165.846153846154 /TABLE/1/25 Yld Functions : plastic strain + strain rate dependency #DIMENSION 2 # FCT_ID strain rate Scale_y 46 0.0 1.00 46 1.0 1.10 46 5.0 1.15 46 10.0 1.20 46 100.0 1.25 46 100000.0 1.35 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #enddata Комментарии ----------- 1. Чтобы описать закон поведения материала картона, следующее ортотропное направление считается. .. image:: images/mat_law112_paper_starter_r_mat_law107_orthotropic_dir.png *(Рисунок 1.)* 2. Эластическое поведение этого материального закона ортотропен. Поведение в плоскости должно быть полностью несвязанный с поведением вне плоскости, вычисляемым как: :math:`{\sigma_{xx}=C_{11}\epsilon_{xx}+C_{12}\epsilon_{yy}\sigma_{yy}=C_{21}\epsilon_{xx}+C_{22}\epsilon_{yy}\sigma_{xy}=G_{12}\gamma_{xy}` С :math:`C=\frac{1}{1−\nu_{12}\nu_{21}}[E_{1}\nu_{12}E_{2}\nu_{21}E_{1}E_{2}]` Компоненты поперечного сдвига вычисляются как: :math:`{\sigma_{yz}=G_{23}\epsilon_{yz}\sigma_{zx}=G_{21}\epsilon_{zx}` Внеплоскостное упругое поведение (для твердых тел только) рассматривается как одноосная эквивалентная задача. Однако Расчет напряжения может различаться при растяжении и сжатии. Эластичность становится нелинейной при сжатии. нагрузки. :math:`\sigma_{zz}=E_{3}\epsilon_{zz}^{e}if\epsilon_{zz}^{e}\ge0\sigma_{zz}=E_{3C}(1−e^{−C_{c}\epsilon_{zz}^{e}})if\epsilon_{zz}^{e}<0` 3. В Ся 2002 г. формулировке критерий текучести в плоскости, обозначаемый как :math:`f` , определяется как: :math:`f=\sumI=16\chi_{I}(\frac{\sigma:N_{I}}{\sigma_{Y}^{I}})^{2k}−1` Где, :math:`\chi_{I}={1if\sigma:N_{I}>00otherwise` :math:`\chi_{I}` Переключение параметров. :math:`\sigma` Тензор напряжений Коши. :math:`N_{I}` Нормальное направление плоскостей текучести. :math:`\sigma_{Y}^{I}` Напряжения текучести. :math:`k` Положительное целое число. Каждое направление связано с определенным направлением нагрузки. в порядке, определенном ниже: 1 Натяжение в ортотропном направлении 1. 2 Натяжение в ортотропном направлении 2. 3 Положительный сдвиг в плоскости. 4 Сжатие в ортотропном направлении 1. 5 Сжатие в ортотропном направлении 2. 6 Отрицательный сдвиг в плоскости (тот же входной сигнал, что и положительный сдвиг в плоскости) срезать :math:`\sigma_{Y}^{6}=\sigma_{Y}^{3}` ). Вектор направления нормали к плоскостям текучести: :math:`N_{1}=[\frac{1}{\sqrt{1+\nu_{1p}^{2}}}−\frac{\nu_{1p}}{\sqrt{1+\nu_{1p}^{2}}}0000]N_{2}=[−\frac{\nu_{2p}}{\sqrt{1+\nu_{2p}^{2}}}\frac{1}{\sqrt{1+\nu_{2p}^{2}}}0000]N_{3}=[000100]N_{4}=[−\frac{1}{\sqrt{1+\nu_{4p}^{2}}}\frac{\nu_{4p}}{\sqrt{1+\nu_{4p}^{2}}}0000]N_{5}=[\frac{\nu_{5p}}{\sqrt{1+\nu_{5p}^{2}}}−\frac{1}{\sqrt{1+\nu_{5p}^{2}}}0000]N_{6}=[000−100]` Каждое направление :math:`I` затем связывается с конкретной доходностью стресс, выражением которого является: :math:`\sigma_{Y}^{I}=S_{I}^{0}+A_{I}^{0}tanh(B_{I}^{0}\epsilon_{p}^{f})+C_{I}^{0}\epsilon_{p}^{f}withI\in[1,6]` Где, :math:`\epsilon_{p}^{f}` это эквивалент пластика в плоскости деформация (связанная с функцией текучести :math:`f` ). Функция выхода вне плоскости обозначается как :math:`g` определяется как: :math:`g=−\sigma_{zz}−\sigma_{Y}^{C}` с :math:`\sigma_{Y}^{C}=A_{\sigma}+B_{\sigma}exp(C_{\sigma}\epsilon_{p}^{g})` Где, :math:`\epsilon_{p}^{g}` это внеплоскостной эквивалент пластика деформация (связанная с функцией текучести :math:`g` ). Выход поперечного сдвига функция: :math:`h=\frac{\sqrt{\sigma_{yz}^{2}+\sigma_{zx}^{2}}}{\sigma_{Y}^{S}}−1` Где, :math:`\sigma_{Y}^{S}=\tau_{0}+[A_{\tau}−min(0,\sigma_{zz})B_{\tau}]\epsilon_{p}^{h}` :math:`\epsilon_{p}^{h}` Эквивалентная пластическая деформация вне плоскости (связанная с функция доходности :math:`h` ). Если выбран вариант табличного предела текучести (Itab = 1), каждый предел текучести равен связанный с таблицей (TAB_YLDi) для определения напряжения эволюция с пластической деформацией, при нескольких скоростях пластической деформации. Два масштабные коэффициенты также могут быть определены в направлении X и Y для каждого стол. В этом случае параметры закалки :math:`S0i` , :math:`A0i` , :math:`B0i` , :math:`C0i` , :math:`A_{\sigma}` , :math:`B_{\sigma}` , :math:`C_{\sigma}` , :math:`\tau_{0}` , :math:`A_{\tau}` и :math:`B_{\tau}` игнорируются, а предел текучести становится: :math:`\sigma_{Y}^{I}=f_{Y}^{tab_YLDI}(\epsilon_{p}^{f},\dot{\epsilon}_{p}^{f})I\in[1,6]\sigma_{Y}^{C}=f_{Y}^{tab_YLDC}(\epsilon_{p}^{g},\dot{\epsilon}_{p}^{g})\sigma_{Y}^{S}=f_{Y}^{tab_YLDS}(\epsilon_{p}^{h},\dot{\epsilon}_{p}^{h})` Для поля вывода эквивалент «Глобальная» пластическая деформация рассчитывается как: :math:`\epsilon_{p}=\sqrt{(\epsilon_{p}^{f})^{2}+(\epsilon_{p}^{g})^{2}+(\epsilon_{p}^{h})^{2}}`