========== /MAT/LAW66 ========== Ключевое слово формата блока Этот закон моделирует изотропное растяжение-сжатие. Закон об упругопластическом материале с использованием определяемых пользователем функций для упрочняющей части напряжение-деформация (пластическая деформация против напряжения). Этот закон можно определить для сжатия и напряжение. Формат ------ .. 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/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/mat_ID/unit_ID", "/MAT/LAW66/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:`\upsilon`", ":math:`\upsilon`", "Мангольд", "Мангольд", "Fcut", "Fcut", "Фгладкий", "Iyld_rate" "Pc", "Pc", "Pt", "Pt", "Ec", "Ec", "RPCT", "RPCT", "", "" Читать только если I yld_rate = 0, 1 or 2 .. 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_IDc", "fct_IDt", "Фскалек", "Фскалек", "Fшкалет", "Fшкалет", "", "", "", "" ":math:`\dot{\epsilon}_{0}`", ":math:`\dot{\epsilon}_{0}`", "c", "c", ":math:`\sigma_{y0}`", ":math:`\sigma_{y0}`", "VP", "", "", "" Читать только если I yld_rate = 3 .. 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_IDc", "fct_IDt", "Фскалек", "Фскалек", "Fшкалет", "Fшкалет", "", "", "", "" "Frate_IDc", "Frate_IDt", "Fscale_ratec", "Fscale_ratec", "Fscale_ratet", "Fscale_ratet", "", "", "", "" Читать только если I yld_rate = 4 .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "NFUNCC", "NFUNCT", "", "", "", "", "", "", "", "" Для каждого NFUNCC .. 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_IDc", "", ":math:`\dot{\epsilon}_{i}^{c}`", ":math:`\dot{\epsilon}_{i}^{c}`", "Фскалек", "Фскалек", "", "", "", "" Для каждого NFUNCT .. 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_IDt", "", ":math:`\dot{\epsilon}_{i}^{t}`", ":math:`\dot{\epsilon}_{i}^{t}`", "Fшкалет", "Fшкалет", "", "", "", "" Определение ----------- .. 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:`\upsilon`", "Коэффициент Пуассона.(Реальный)", "" "Мангольд", "Коэффициент закалки. = 0 Закалка – полностью изотропная модель. = 1 Для закалки используется кинематическая модель Прагера-Циглера. = между 0 и 1 Упрочнение интерполируется между двумя моделями. (Настоящий)", "" "Фгладкий", "Флаг опции плавной скорости деформации. = 0 (по умолчанию) Нет сглаживания скорости деформации. = 1 Сглаживание скорости деформации активно. (Целое число)", "" "Fcut", "Частота среза для скорости деформации фильтрация, Приложение: Filtering.Default = 1030 (Реал)", ":math:`[Hz]`" "Iyld_rate", "Влияние скорости на флаг предела текучести. = 1 (по умолчанию) Использование Каупер-Саймондса: :math:`1+(\frac{\dot{\epsilon}}{\dot{\epsilon}_{0}})^{\frac{1}{c}}` = 2 Используя: :math:`1+cLn(\frac{\dot{\epsilon}}{\dot{\epsilon}_{0}})` = 3 Используя две кривые нагрузки для масштабирования предела текучести. (fct_IDc) при сжатии и напряжение (fct_IDt). = 4 Используя разные функции сжатия и натяжения для разных значения скорости деформации. (Целое число)", "" "Pc", "Предельное давление при сжатии. По умолчанию = 0 (Реальный)", ":math:`[Pa]`" "Pt", "Предельное давление при растяжении. По умолчанию = 0 (Реал)", ":math:`[Pa]`" "Ec", "(Необязательно) Модуль Юнга сжатия. 2(Реал)", ":math:`[Pa]`" "RPCT", "Масштабный коэффициент, используемый на ПК и Пт. 2(Реал)", "" "fct_IDc", "Выход сжатия стресс.(Целое число)", "" "fct_IDt", "Выход напряжения стресс.(Целое число)", "" "Фскалек", "Масштабный коэффициент для ординаты (напряжения) в fct_IDc.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" "Fшкалет", "Масштабный коэффициент для ординаты (напряжения) в fct_IDt.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" "c", "Скорость деформации параметр.(Реальный)", "" ":math:`\dot{\epsilon}_{0}`", "Эталонная скорость деформации. По умолчанию = 1,0. (Реал)", ":math:`[\frac{1}{s}]`" ":math:`\sigma_{y0}`", "Начальный предел текучести. По умолчанию = 0. (Реал)", ":math:`[Pa]`" "VP", "Флаг выбора скорости деформации. = 0 Влияние скорости деформации на предел текучести зависит от общей деформации. ставка. = 1 Влияние скорости деформации на предел текучести зависит от пластической деформации. ставка. В этом случае фильтрация скорости деформации отсутствует, поэтому Fsmooth и Fcut не используются. Доступно только если Iyld_rate = 1 (Каупер Саймондс). 3 (Целое число)", "" "Frate_IDc", "Функция влияния скорости деформации сжатия идентификатор.(Целое число)", "" "Frate_IDt", "Функция эффекта скорости деформации растяжения идентификатор.(Целое число)", "" "Fscale_ratec", "Масштабный коэффициент для ординаты (напряжения) в Frate_IDc.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" "Fscale_ratet", "Масштабный коэффициент для ординаты (напряжения) в Frate_IDt.По умолчанию = 1,0. (Реал)", ":math:`[Pa]`" "NFUNCC", "Число сжатий функция.(Целое число)", "" "NFUNCT", "Количество напряжений функция.(Целое число)", "" ":math:`\dot{\epsilon}_{i}^{c}`", "это скорость деформации сжатия i =1,%%CP0000%.(Реальный)", ":math:`[\frac{1}{s}]`" ":math:`\dot{\epsilon}_{i}^{t}`", "это скорость деформации растяжения i=1,%%CP0000%.(Реальный)", ":math:`[\frac{1}{s}]`" Пример (алюминий) ----------------- .. code-block:: #RADIOSS STARTER /UNIT/1 unit for mat g mm ms #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #- 2. MATERIALS: #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /MAT/LAW66/1/1 Aluminum # RHO_I .0027 # E Nu C_hard F_cut F_smooth Iyld_rate 60400 .33 0 0 0 4 # P_c P_t 500 600 # NFUNCC NFUNCT 2 2 #funct_IDc Epsilon_c Fscalec 38 10 1 40 40 1.6 #funct_IDt Epsilon_t Fscalet 38 10 1 40 40 1.6 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #- 3. FUNCTIONS: #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/38 function_38 # X Y 0 90 .08 170 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/40 function_40 # X Y 0 90 .08 170 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| #ENDDATA /END #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| Пример (дополнительный модуль сжатия Юнга) ------------------------------------------ .. code-block:: #RADIOSS STARTER #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| ## Material Law #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /MAT/LAW66/1 20MAT124_20degree # RHO_I 1.25000000000000E-09 # E Nu C_hard F_cut F_smooth Iyld_rate 210000.0 0.33 0.0 0.0 0 1 # P_c P_t EC RPCT 0.0 0.0 70000. #funct_IDc funct_IDt Fscalec Fscalet 34 34 100.0 200.0 # Epsilon_0 c Sigma_Y0 VP 0.0 0.0 0.0 0 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| ##HWCOLOR curves 24 9 /FUNCT/24 20deg_TENSION # X Y 0.0 0.200535124 1.11148000000000E-04 0.23893938 2.47763000000000E-04 0.274602617 3.84248000000000E-04 0.308164207 5.81974000000000E-04 0.354639794 7.98180000000000E-04 0.395218545 1.03131900000000E-03 0.43242102 0.001382833 0.473801266 0.001747862 0.508073173 0.002139804 0.539268776 0.002704889 0.577068218 0.003299031 0.610692299 0.004078009 0.646563755 0.005715537 0.702190108 1.0 0.75 #---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----| /FUNCT/34 20deg_COMPRESSION # X Y 0.0 0.709520996 0.002556758 0.768542475 0.005112572 0.814793042 0.007700089 0.850991269 0.010231756 0.877769462 0.012846771 0.896395146 0.015380438 0.909432849 0.020500175 0.924867455 0.030979299 0.948100519 0.073545677 1.001804424 0.126630132 1.075808748 0.245056814 1.17707424 1.0 1.2 #enddata Комментарии ----------- 1. Это изотропный упругопластический закон. Предел текучести определяется с помощью сжатия и Предел текучести при растяжении по сравнению с эффективной пластической деформацией для обоих (сжатия и растяжения). При превышении оба давления P t и P c , определить, является ли предел текучести при растяжении или предел текучести при сжатии используется соответственно. Если давление находится между этими двумя значения, предел текучести определяется по формуле: Если :math:`−P_{t}\leP\leP_{c}` :math:`\sigma_{y}=\alpha\sigma_{y}^{t}(\epsilon_{p})+(1−\alpha)\sigma_{y}^{c}(\epsilon_{p})\alpha=\frac{P_{c}−P}{P_{e}+P_{t}}` Если :math:`P_{t}=P_{c}=0` , или давление выходит за пределы диапазона двух значений, выход стресс создается: :math:`\sigma_{y}=\sigma_{y}^{t}(\epsilon_{p})` если :math:`P\le0` :math:`\sigma_{y}=\sigma_{y}^{c}(\epsilon_{p})` если :math:`P>0` 2. If Ec определяется, Модуль Юнга определяется как: - Модуль Юнга равен E, если P > -RPCT * Пт - Модуль Юнга равен Ec, если P < RPCT * ПК - Линейная интерполяция выполняется между E и Ec, если -RPCT * Pt < P < RPCT * ПК 3. Предел текучести рассчитывается как: If VP= 1: :math:`\sigma_{y}(\epsilon_{p},\dot{\epsilon}_{p})=\sigma_{y}^{s}(\epsilon_{p})+\sigma_{y0}(\frac{\dot{\epsilon}_{p}}{\dot{\epsilon}_{0}})^{\frac{1}{c}}` если :math:`\sigma_{y0}>0` :math:`\sigma_{y}(\epsilon_{p},\dot{\epsilon}_{p})=\sigma_{y}^{s}(\epsilon_{p})[1+(\frac{\dot{\epsilon}_{p}}{\dot{\epsilon}_{0}})^{\frac{1}{c}}]` если :math:`\sigma_{y0}=0` If VP= 0: :math:`\sigma_{y}(\epsilon_{p},\dot{\epsilon}_{p})=\sigma_{y}^{s}(\epsilon_{p})` если :math:`\sigma_{y0}>0` :math:`\sigma_{y}(\epsilon_{p},\dot{\epsilon}_{p})=\sigma_{y}^{s}(\epsilon_{p})` если :math:`\sigma_{y0}=0` с :math:`\sigma_{y}^{t}(\epsilon_{p})` статический предел текучести и :math:`\sigma_{y0}` является начальным пределом текучести. 4. /VISC/PRONY can be used with this material law to include viscous effects.