================= /MAT/LAW76 (SAMP) ================= Ключевое слово формата блока Этот закон описывает полуаналитический упругопластический материал с использованием определяемых пользователем функций для наклепанной части для растяжение, сжатие и сдвиг (напряжение как функция деформации). Формат ------ .. 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/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/mat_ID/unit_ID", "/MAT/LAW76/mat_ID/unit_ID or /MAT/SAMP/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`", "", "", "", "", "", "" "tab_IDt", "tab_IDc", "tab_IDs", "", "", "", "", "", "", "" "Fшкалет", "Fшкалет", "Фскалек", "Фскалек", "Fмасштабы", "Fмасштабы", "", "", "XFAC", "XFAC" ":math:`\nu_{p}`", ":math:`\nu_{p}`", "fct_IDpr", "Fscalepr", "Fscalepr", "Фгладкий", "Fcut", "Fcut", "", "" ":math:`\epsilon_{p}^{f}`", ":math:`\epsilon_{p}^{f}`", ":math:`\epsilon_{p}^{r}`", ":math:`\epsilon_{p}^{r}`", "", "", "", "", "", "" "fct_ID1", "", "", "Fмасштаб1", "Fмасштаб1", "", "", "", "", "" "Яформа", "IQUAD", "ICONV", "", "", "", "", "", "", "" Определение ----------- .. 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`", "Пуассона соотношение.(Реальное)", "" "tab_IDt", "Таблица пределов текучести при растяжении идентификатор (напряжение в сравнении с пластической деформацией растяжения с возможность зависимости скорости деформации).(Целое число)", "" "tab_IDc", "Предел текучести при сжатии идентификатор таблицы (напряжение в сравнении с пластической деформацией сжатия с возможность зависимости скорости деформации).(Целое число)", "" "tab_IDs", "Таблица пределов текучести при сдвиге идентификатор (напряжение в сравнении с пластической деформацией сдвига с возможностью зависимости скорости деформации).(Целое число)", "" "Fшкалет", "Масштабный коэффициент для ординаты (стресс) для tab_IDt.По умолчанию = 1,0 (Реальный)", ":math:`[Pa]`" "Фскалек", "Масштабный коэффициент для ординаты (стресс) для tab_IDc.По умолчанию = 1,0 (Реальный)", ":math:`[Pa]`" "Fмасштабы", "Масштабный коэффициент для ординаты (стресс) для tab_IDs.По умолчанию = 1,0 (Реальный)", ":math:`[Pa]`" "XFAC", "Масштабный коэффициент для вторая запись (скорость деформации) из трех таблиц (tab_IDt, tab_IDc и tab_IDs). 6По умолчанию = 1,0 (реальное)", "" ":math:`\nu_{p}`", "Пластиковый Пуассон соотношение.(Реальное)", "" "fct_IDpr", "Пластиковый коэффициент Пуассона идентификатор функции ( :math:`\nu_{p}` против пластика штамм).(Реальный)", "" "Fscalepr", "Масштабный коэффициент для ординаты ( :math:`\nu_{p}` ) в fct_IDpr.По умолчанию = 1,0 (Реальный)", "" "Фгладкий", "Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)", "" "Fcut", "Частота среза для фильтрация скорости деформации. По умолчанию = 1030 (Реал)", ":math:`[Hz]`" ":math:`\epsilon_{p}^{f}`", "Неудачная пластическая деформация (начало элементарного урона). По умолчанию = 2e30 (Реальное)", "" ":math:`\epsilon_{p}^{r}`", "Максимальная пластическая деформация (элемент удален). По умолчанию = 2e30 (реальный)", "" "fct_ID1", "Идентификатор функции повреждения (повреждение в зависимости от пластической деформации). 2 (целое число)", "" "Fмасштаб1", "Масштабный коэффициент для ординаты для fct_ID1. 2По умолчанию = 1,0 (Реальное)", "" "Яформа", "Флаг формулировки. 4 = 0 (по умолчанию) Нет соответствующей формулировки. = 1 Связанная с ним формулировка фон Мизеса. (Целое число)", "" "IQUAD", "Флаг поверхности текучести. 3 = 0 (по умолчанию) Поверхность текучести линейна по фон Мизесу. = 1 Поверхность текучести квадратична по Мизесу. (рекомендуется). (Целое число)", "" "ICONV", "Флаг условия выпуклости. =0 (по умолчанию) Никакой обработки для обеспечения стабильности материала. = 1 Выпуклость предела текучести (стабильность материала) равна уверен. (Целое число)", "" Пример (Материал) ----------------- .. 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/LAW76/1/1 LAW76_Material # RHO_I 1E-6 # E nu 100.0 .3 # TAB_IDt TAB_IDc TAB_IDs 1000 1001 1003 # Fscale_t Fscale_c Fscale_s XFAC 1.000 1.000 1.000 1.000 # Nu_p fct_IDpr Fscale_pr Fsmooth Fcut 0.5 0 0 1 1e30 # EPS_f_p EPS_r_p 0 0 #funct_ID1 Fscale_1 0 # IFORM IQUAD ICONV 0 0 1 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /TABLE/1/1000 curve_list TENSION strain rates 2 10010 1.0e-4 10020 1.0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/10010 eps_vs_sigma funct dt=1.0e-4 0.0000 .100000 1.0000 .200000 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/10020 eps_vs_sigma funct dt=1.0e-4 0.0000 .100000 1.0000 .200000 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /TABLE/1/1001 curve_list COMPRESSION strain rates 2 10030 1.0e-4 10040 1.0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/10030 eps_vs_sigma funct dt=1.0e-4 0.0000 .200000 1.0000 .400000 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/10040 eps_vs_sigma funct dt=1.0e-4 0.0000 .200000 1.0000 .400000 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /TABLE/1/1003 curve_list SHEAR strain rates 2 10050 1.0e-4 10060 1.0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/10050 eps_vs_sigma funct dt=1.0e-4 0.0000 .050000 0.5000 .060000 1.0000 .065000 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/10060 eps_vs_sigma funct dt=1.0e-4 0.0000 .050000 0.5000 .060000 1.0000 .065000 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #ENDDATA /END Комментарии ----------- 1. Это материал совместим с панцирем, толстой оболочкой и твердыми элементами. 2. Материальный ущерб может быть смоделировал двумя способами: - :math:`\epsilon_{p}^{f}` (начало повреждать) и :math:`\epsilon_{p}^{r}` (удаление элемента): :math:`D=\frac{\epsilon_{p}-\epsilon_{p}^{f}}{\epsilon_{p}^{r}-\epsilon_{p}^{f}}` Где, :math:`\epsilon_{p}` это текущий накопленный пластик напряжение. - Функция повреждения fct_ID 1 .. image:: images/mat_law76_samp_starter_r_fail_tab1_damage.png *(Рисунок 1.)* Если используется функция повреждения fct_ID1, то :math:`\epsilon_{p}^{f}` и :math:`\epsilon_{p}^{r}` будет проигнорирован. Ущерб переменная затем влияет на жесткость материала и, таким образом, вычисление тензора напряжений как: :math:`\sigma=\sigma_{eff}1-D` Где, :math:`\sigma_{eff}` Неповрежденное (эффективное) тензорное напряжение :math:`\sigma` Поврежденный (текущий) тензор напряжений :math:`D` Переменная урона 3. Выбор поверхности текучести: :math:`f={\sigma_{VM}−A_{0}−A_{1}​P−A_{2}P^{2}​IQUAD=0\sigma_{VM}^{2}−A_{0}−A​​​​​1P−A_{2}P^{2}IQUAD=1` Где, :math:`P=−\frac{\sigma_{xx}+\sigma_{yy}+\sigma_{zz}}{3}` :math:`\sigma_{VM}=\sqrt{\frac{3}{2}[(\sigma_{xx}+P)^{2}+(\sigma_{yy}+P)^{2}+(\sigma_{zz}+P)^{2}+2\sigma_{xy}^{2}+2\sigma_{yz}^{2}+2\sigma_{xz}^{2}]}` :math:`A_{0}` , :math:`A_{1}` и :math:`A_{2}` коэффициенты рассчитываются по упрочнению кривая растяжения, сжатия и сдвига. Для фон Мизеса и Поверхность текучести Друкера-Прагера, IQUAD=0 можно использовать. Однако в некоторых ситуациях Radioss может быть сложно адаптироваться к :math:`A_{0}` , :math:`A_{1}` и :math:`A_{2}` коэффициенты при использовании IQUAD=0, и более простая подгонка получено с использованием IQUAD=1. 4. Выбор пластика формулировка: - Для формулировки, не связанной с пластичностью, I форма =0: Функция правила пластического потока, :math:`g` , используется для описания пластической деформации приращение :math:`d\epsilon_{p}=d\lambda\frac{\partialg}{\partial\sigma}` . В этом случае :math:`\frac{\partialg}{\partial\sigma}` не является нормальным к поверхности текучести :math:`f` и :math:`g` не связан с поверхностью текучести :math:`f` . .. image:: images/mat_law76_samp_starter_r_mat_law76_plasticity_formulation.png *(Рисунок 2.)* Правило пластического течения, :math:`g` , определяется: :math:`g=\sqrt{\sigma_{VM}^{2}+\alphaP^{2}}` Для таких материалов, как почва или камень, обычно используется знак «нет». связанная формула пластичности, Iform=0. - Для связанной пластичности: I форма = 1, :math:`g=f` В этом случае скорость пластической деформации равна функция вектора нормали поверхности текучести :math:`f` . В таких материалах, как металл, обычно используется соответствующая формула пластичности. :math:`d\epsilon_{p}=d\lambda\frac{\partialf}{\partial\sigma}=d\lambda\frac{\partialg}{\partial\sigma}` .. image:: images/mat_law76_samp_starter_r_mat_law76_iform1.png *(Рисунок 3.)* 5. Выпуклость флаг условия ICONV=1 используется для обеспечения стабильность материального закона за счет выпуклости поверхности текучести. Доходность Поверхность может быть гиперболической при низких значениях текучести при сдвиге при растяжении и сжатии. В в этом случае нет однозначного решения и Radioss будет обновить (увеличить) предел текучести при сдвиге, чтобы обеспечить выпуклость поверхности текучести. Следовательно, предел текучести при сдвиге может отличаться от входной кривой. 6. В таблицах должно быть максимум размерность равна 2. Первая запись — это пластическая деформация, вторая запись — скорость деформации. 7. Пользовательские переменные USR2, USR3, USR4 используются для вывода компонентов пластической деформации при растяжении, сжатии и сдвиге. вывод доступен как для оболочек, так и для твердых тел в истории времени и в анимации. файл.