LS-DYNA混凝土材料参数

12th International LS-DYNA® Users Conference
Constitutive Modeling(1)
Performance of LS-DYNA® Concrete Constitutive Models
Youcai Wu, John E. Crawford, Joseph M. Magallanes
1 Introduction
Concrete is a very common material in modern civil constructions, such as highways, bridges, skyscrapers, etc. The safety of these structures under blast, impact, and ballistic loads has been one of the primary concerns of designers in recent years. However, full-scale structural response test data is costly and difficult to obtain for these types of loads. As a consequence, responses predicted by physics-based numerical analyses have been important resources for both academics and structural engineers to determine the behaviors of reinforced concrete (RC) structures under these loads. Consequently, practitioners are continuously looking for verified and validated numerical models, so that consistent sets of virtual response data can be generated, and, from which one can study structural behaviors, formulate simplified engineering models, and develop new structural and material designs to have improved resistance to these loads. To realistically predict the behaviors of RC structures under various types of loads, the concrete constitutive model needs to be shown to simulate known behaviors at smaller material specimen levels up to the full-scale structural level. Although it is very difficult to predict the behavior of concrete exactly, the constitutive model should capture the most basic behaviors of concrete. Advances in finite element (FE) methods and material constitutive modelings have made it feasible to support engineers’ daily requirements in designs and assessments. Among the readily available FE softwares, LS-DYNA® is widely applied in analyzing structural responses to shock and impact loads, and it provides a variety of concrete constitutive models, such as Mat_Pseudo_Tensor (MAT016), Mat_Geologic_Cap_Model (MAT025), Mat_Concrete_Damage (MAT072), Mat_Soil_Concrete (MAT078), Mat_Winfrith_Concrete (MAT084), Mat_Brittle_Damage (MAT096), Mat_Johson_Holmquist_Concrete (MAT111), Mat_CSCM_Concrete (MAT159). Each of these models has its own advantages and disadvantages, therefore, FE analysts are actually facing a wide range of choices but they may

LS-DYNA混凝土材料参数Material Name - CONC-35MPA(水泥混凝土) Equation of State P alphaReference density 2.75000E+00 (g/cm3 ) Porous density 2.31400E+00 (g/cm3 )Porous soundspeed 2.92000E-01 (cm/us ) Initial compaction pressure 2.33000E-04 (Mbar ) Solid compaction pressure 6.00000E-02 (Mbar ) Compaction exponent 3.00000E+00 (none ) Solid EOS PolynomialBulk Modulus A1 3.52700E-01 (Mbar ) Parameter A2 3.95800E-01 (Mbar ) Parameter A3 9.04000E-02 (Mbar ) Parameter B0 1.22000E+00 (none ) Parameter B1 1.22000E+00 (none ) Parameter T1 3.52700E-01 (Mbar )Parameter T2 0.00000E+00 (Mbar ) Reference Temperature 3.00000E+02 (K ) Specific Heat 6.54000E-06 (Terg/gK )Thermal Conductivity 0.00000E+00 (Terg/cmKus) Compaction Curve StandardStrength RHT ConcreteShear Modulus 1.67000E-01 (Mbar ) Compressive Strength (fc) 3.50000E-04 (Mbar ) Tensile Strength (ft/fc) 1.00000E-01 (none ) Shear Strength (fs/fc) 1.80000E-01 (none ) Intact Failure Surface1.60000E+00 (none )Constant AIntact Failure Surface6.10000E-01 (none )Exponent NTens./Comp. Meridian Ratio6.80500E-01 (none )(Q)Brittle to Ductile Transition 1.05000E-02 (none )G (elas.)/(elas.-plas.) 2.00000E+00 (none )Elastic Strength / ft 7.00000E-01 (none )Elastic Strength / fc 5.30000E-01 (none )Fractured Strength Constant B 1.60000E+00 (none )Fractured Strength Exponent6.10000E-01 (none )MCompressive Strain Rate Exp.3.20000E-02 (none )AlphaTensile Strain Rate Exp. Delta 3.60000E-02 (none )Max. Fracture Strength Ratio 1.00000E+20 (none )Use CAP on Elastic Surface? Y esDamage Constant, D1 4.00000E-02 (none ) Damage Constant, D2 1.00000E+00 (none ) Minimum Strain to Failure 1.00000E-02 (none ) Residual Shear ModulusFraction1.30000E-01 (none ) Tensile Failure Hydro (Pmin) Erosion Geometric Strain Erosion Strain2.00000E+00 (none ) Type of Geometric Strain Instantaneous Material Cutoffs-Maximum Expansion 1.00000E-01 (none ) Minimum Density Factor(Euler)1.00000E-04 (none ) Minimum Density Factor(SPH)2.00000E-01 (none ) Maximum Density Factor(SPH)3.00000E+00 (none ) Minimum Soundspeed 1.00000E-10 (cm/us ) Maximum Soundspeed 1.01000E+20 (cm/us ) Maximum Temperature 1.01000E+20 (K )Reference:Riedel W. "Beton unter dynamischen Lasten..." Phd Thesis. EMI. Freiburg. Germany. 2000Material Name - CONC140MPA（高强度很混凝土）Equation of State P alphaReference density 2.75000E+00 (g/cm3 ) Porous density 2.52000E+00 (g/cm3 ) Porous soundspeed 3.24200E-01 (cm/us ) Initial compaction pressure 9.33000E-04 (Mbar )Solid compaction pressure 6.00000E-02 (Mbar ) Compaction exponent 3.00000E+00 (none ) Solid EOS Polynomial Bulk Modulus A1 3.52700E-01 (Mbar ) Parameter A2 3.95800E-01 (Mbar ) Parameter A3 9.04000E-02 (Mbar ) Parameter B0 1.22000E+00 (none ) Parameter B1 1.22000E+00 (none ) Parameter T1 3.52700E-01 (Mbar ) Parameter T2 0.00000E+00 (Mbar ) Reference Temperature 3.00000E+02 (K ) Specific Heat 6.54000E-06 (Terg/gK ) Thermal Conductivity 0.00000E+00 (Terg/cmKus) Compaction Curve StandardShear Modulus 2.20600E-01 (Mbar )Compressive Strength (fc) 1.40000E-03 (Mbar )Tensile Strength (ft/fc) 1.00000E-01 (none )Shear Strength (fs/fc) 1.80000E-01 (none )Intact Failure Surface Constant1.60000E+00 (none )AIntact Failure Surface Exponent6.10000E-01 (none )NTens./Comp. Meridian Ratio6.80500E-01 (none )(Q)Brittle to Ductile Transition 1.05000E-02 (none )G (elas.)/(elas.-plas.) 2.00000E+00 (none )Elastic Strength / ft 7.00000E-01 (none )Elastic Strength / fc 5.30000E-01 (none ) Fractured Strength Constant B 1.60000E+00 (none ) Fractured Strength Exponent M 6.10000E-01 (none ) Compressive Strain Rate Exp.9.09000E-03 (none )AlphaTensile Strain Rate Exp. Delta 1.25000E-02 (none ) Max. Fracture Strength Ratio 1.00000E+20 (none ) Use CAP on Elastic Surface? Y esFailure RHT ConcreteDamage Constant, D1 4.00000E-02 (none ) Damage Constant, D2 1.00000E+00 (none ) Minimum Strain to Failure 1.00000E-02 (none ) Residual Shear Modulus1.30000E-01 (none )FractionTensile Failure Hydro (Pmin)Erosion Geometric StrainErosion Strain 2.00000E+00 (none )Type of Geometric Strain InstantaneousMaterial Cutoffs-Maximum Expansion 1.00000E-01 (none )Minimum Density Factor1.00000E-04 (none )(Euler)Minimum Density Factor2.00000E-01 (none )(SPH)Maximum Density Factor3.00000E+00 (none )(SPH)Minimum Soundspeed 1.00000E-10 (cm/us )Maximum Soundspeed 1.01000E+20 (cm/us )Maximum Temperature 1.01000E+20 (K )Reference:Riedel W. "Beton unter dynamischen Lasten..." Phd Thesis. EMI.Freiburg. Germany. 2000Material Name - CONCRETE-L贫混凝土Equation of State PorousReference density 2.44000E+00 (g/cm3 ) Solid Soundspeed 2.20000E-01 (cm/us ) Porous Soundspeed 2.20000E-01 (cm/us ) Density #1 2.34000E+00 (g/cm3 ) Density #2 2.35000E+00 (g/cm3 ) Density #3 2.40000E+00 (g/cm3 ) Density #4 2.46000E+00 (g/cm3 ) Density #5 2.50000E+00 (g/cm3 ) Density #6 0.00000E+00 (g/cm3 ) Density #7 0.00000E+00 (g/cm3 ) Density #8 0.00000E+00 (g/cm3 ) Density #9 0.00000E+00 (g/cm3 ) Density #10 0.00000E+00 (g/cm3 ) Pressure #1 0.00000E+00 (Mbar ) Pressure #2 2.50000E-04 (Mbar ) Pressure #3 7.00000E-04 (Mbar ) Pressure #4 1.30000E-03 (Mbar ) Pressure #5 2.50000E-03 (Mbar ) Pressure #6 0.00000E+00 (Mbar )Pressure #7 0.00000E+00 (Mbar ) Pressure #8 0.00000E+00 (Mbar ) Pressure #9 0.00000E+00 (Mbar ) Pressure #10 0.00000E+00 (Mbar ) Crush failure density 1.01000E+20 (g/cm3 ) Strength Drucker-PragerShear Modulus 7.88000E-02 (Mbar ) Pressure hardening type PiecewisePressure #1 0.00000E+00 (Mbar ) Pressure #2 8.00000E-04 (Mbar ) Pressure #3 1.10000E-03 (Mbar ) Pressure #4 2.00000E-03 (Mbar ) Pressure #5 0.00000E+00 (Mbar ) Pressure #6 0.00000E+00 (Mbar ) Pressure #7 0.00000E+00 (Mbar ) Pressure #8 0.00000E+00 (Mbar ) Pressure #9 0.00000E+00 (Mbar ) Pressure #10 0.00000E+00 (Mbar ) Yield Stress (zero plastic2.50000E-04 (Mbar ) strain)Yield Stress #2 1.10000E-03 (Mbar )Yield Stress #3 1.60000E-03 (Mbar ) Yield Stress #4 1.95000E-03 (Mbar ) Yield Stress #5 0.00000E+00 (Mbar ) Yield Stress #6 0.00000E+00 (Mbar ) Yield Stress #7 0.00000E+00 (Mbar ) Yield Stress #8 0.00000E+00 (Mbar ) Yield Stress #9 0.00000E+00 (Mbar ) Yield Stress #10 0.00000E+00 (Mbar ) Failure Hydro (Pmin) Hydro Tensile Limit -2.50000E-05 (Mbar ) Reheal Y esCrack Softening NoStochastic failure NoErosion Geometric Strain Erosion Strain 4.00000E+00 (none ) Type of Geometric Strain Instantaneous Material Cutoffs- Maximum Expansion 1.00000E-01 (none ) Minimum Density Factor(Euler)1.00000E-04 (none ) Minimum Density Factor(SPH)2.00000E-01 (none ) Maximum Density Factor(SPH)3.00000E+00 (none ) Minimum Soundspeed 1.00000E-10 (cm/us ) Maximum Soundspeed 1.01000E+20 (cm/us ) Maximum Temperature 1.01000E+20 (K )Reference:C.M. Wentzel et al "Concrete under Multi-Axial Dynamic Loading" - Tran. Load. & resp. Struct. 19 98。

lsdyna混凝土072材料对应的eos状态方程在LS-DYNA中，混凝土材料的EOS（Equation of State，状态方程）可以通过定义材料卡和定义参数卡来确定。

具体而言，若要为混凝土材料定义EOS状态方程，需要在材料卡中指定mat_eos_type、mat_eos_param和mat_eos_const参数。

mat_eos_type参数是一个整数值，用于指定EOS状态方程的类型。

对于混凝土材料，通常使用EOS模型的类型为15或72。

因此，在这种情况下，需要将mat_eos_type设置为72。

mat_eos_param参数用于指定EOS状态方程的参数。

对于混凝土材料，EOS参数包括密度和泊松比。

这些参数通常根据混凝土的实际特性进行实验测量或估算。

可以将这些参数设置为mat_eos_param(1)和mat_eos_param(2)。

mat_eos_const参数是一个整数值，用于指定额外的EOS状态方程参数。

对于混凝土材料，该参数通常保持为0即可。

因此，对于混凝土材料"072"，可以在材料卡中定义如下的EOS状态方程：*MAT_CONCRETE_DAMAGE_TITLE$# 应变速率效应模型$# 弹性材料模量参考应变率2100/s*MAT_ADD_ERATE_RATE_DEPENDENCE0 14.3*MATERIAL_CONCRETE_DAMAGE_TITLE12 1.9 0.30 1 65.0 2.4E+7*MAT_CONCRETE_DAMAGE8.00E8 4.80E8 00 0.0 0.00 0.0 0.0*EOS72 2250.0 0.20 0其中，*MAT_ADD_ERATE_RATE_DEPENDENCE用于定义应变速率效应模型（可选），*MATERIAL_CONCRETE_DAMAGE用于定义混凝土的损伤模型（可选）。

请注意，以上仅为示例，具体的EOS参数需要根据实际情况进行定义和调整。

登录协同仿真时代
62 *MAT_VISCOUS_FOAM
粘性泡沫 : • 在碰撞模拟中吸能泡沫常用于模拟吸能材料（例 如假人） • 仅用于压缩载荷下的实体模型 • 由平行的非线性弹簧和粘性阻尼构成 • 需要有DENS, EX (初始杨氏模量E1) 和 NUXY • 弹性刚度E’ 定义为: E’ = E1 V-n1 其中 n1 是弹性刚度幂率 • 粘性系数V’ 定义为: V’ = V2 |1-V|n2 其中 V2 = 初始粘度系数 n2 = 粘度系数幂率
*mat_null+EOS(常见流体) *mat_high_explosive+EOS（炸药） *MAT_ELASTIC_PLASTIC_HYDRO+EOS（固体在高应变率下 的流体表现） *mat_3-parameter_barlat（三参数的BARLAT各向异性塑性金属 材料） *mat_transversely_elastic_plastic（横向各向异性金属材料） *mat_fabric（气囊、丝织物） *mat_composite_failure_options_model（实体或壳复合材料） *mat_seatbelt（安全带材料） *mat_thermal_options （材料热物理参数） *user_defined_material_models（用户自定义材料的参数输入）
登录协同仿真时代
泡沫
5 *MAT_SOIL_AND_FOAM 14 *MAT_SOIL_AND_FOAM_FAILURE 26 *MAT_HONEYCOMB 38 *MAT_BLATZ_FO_FOAM 53 *MAT_CLOSED_CELL_FOAM 57 *MAT_LOW_DENSITY_FOAM 62 *MAT_VISCOUS_FOAM 63 *MAT_CRUSHABLE_FOAM 75 *MAT_BILKHU/DUBOIS_FOAM 83 *MAT_FU_CHANG_FOAM 仅对体单元