ABAQUS强度仿真：从本构模型到失效准则的完整评估

ABAQUS在强度仿真中的优势在于其丰富的非线性本构模型库和强大的接触算法。与ANSYS偏重工程化不同，ABAQUS更深入材料的微观力学行为，支持用户自定义材料（UMAT/VUMAT）扩展，这使得它在复合材料、断裂力学和冲击动力学领域有独特地位。

强度仿真的层次框架

强度仿真不是单一分析，而是一个从宏观到微观的评估体系：

第一层：线弹性强度校核 — 计算von Mises应力或主应力，与材料许用应力比较。适用于小载荷、小变形的常规设计校核，大部分结构件的第一轮评估。

第二层：弹塑性极限分析 — 超过屈服后的塑性应变分布和极限载荷。需要定义完整的应力-应变曲线，判断结构是否发生塑性失稳。

第三层：失效与断裂分析 — 裂纹起裂、扩展路径、断裂韧性评估。需要定义损伤起始准则和损伤演化规律。

第四层：疲劳寿命预测 — 循环载荷下的裂纹萌生和扩展寿命。需要S-N曲线或Paris定律参数。

ABAQUS中的材料本构模型

金属弹塑性

ABAQUS的弹塑性材料定义以”屈服面+流动准则+硬化规律”三要素为核心：

屈服准则：

• von Mises（默认）：适用于各向同性金属
• Hill48：适用于各向异性金属（轧制板材、冲压件）
• Barlat89/2004：适用于铝合金的各向异性屈服

流动准则：

• 关联流动（默认）：塑性应变增量方向垂直于屈服面
• 非关联流动：塑性势与屈服函数不同，适用于岩土材料

硬化规律：

• 等向硬化（ISOTROPIC）：屈服面均匀膨胀，适用于单调加载
• 随动硬化（KINEMATIC）：屈服面平移不膨胀，适用于循环加载的Bauschinger效应
• 混合硬化（COMBINED）：膨胀+平移，最接近真实金属行为

ABAQUS输入格式（以等向硬化为例）：

*MATERIAL, NAME=STEEL
*ELASTIC
210000, 0.3
*PLASTIC, HARDENING=ISOTROPIC
345, 0.0
380, 0.02
420, 0.05
480, 0.1
550, 0.2

每行是屈服应力-等效塑性应变对。第一行必须是从零塑性应变开始（屈服点）。数据点需要覆盖到预期最大应变水平，超出范围ABAQUS会做线性外推，可能导致非物理行为。

损伤与失效

ABAQUS的损伤力学模型分两步定义：

损伤起始准则：

• ductile（延性）：基于应力三轴度和应变率的孔洞形核准则
• shear（剪切）：基于剪应力比和应变率的剪切带准则
• Johnson-Cook：温度和应变率相关的经验公式
• Hashin：复合材料分层准则

损伤演化规律：

• 位移控制：达到临界位移δ_f后单元完全失效
• 能量控制：达到断裂能G_c后单元失效
• 线性/指数软化：损伤变量D从0到1的演化方式

当D=1时，单元失去承载能力并被删除（删单元法）。这种方法模拟裂纹扩展不需要预定义裂纹路径，但结果对网格尺寸敏感——网格越细，裂纹路径越精确，但计算量越大。

接触算法与摩擦

ABAQUS的接触算法是其核心竞争力之一：

接触离散化

• 点-面接触（Node-to-Surface）：从面节点与主面段接触。传统方法，计算效率高，但可能穿透主面。
• 面-面接触（Surface-to-Surface）：考虑从面和主面的整体形状。精度更高，穿透更小，推荐用于大多数分析。

接触公式

• 罚函数法（Penalty）：在接触面间引入虚拟弹簧，穿透量与接触力成正比。计算稳定，是最常用的方法。弹簧刚度需要足够大以限制穿透，但不能太大导致收敛困难。ABAQUS自动选择刚度，但可以手动调整。
• 拉格朗日乘子法（Lagrange Multiplier）：精确满足非穿透条件，不产生穿透。但增加额外自由度，计算量大，且可能引起振荡。
• 增广拉格朗日法（Augmented Lagrange）：罚函数+拉格朗日修正，平衡精度和效率。

摩擦模型

库仑摩擦是最常用的模型：τ = μ*σ_n（μ为摩擦系数）。ABAQUS支持弹性滑移（弹性阶段小滑移）和有限滑移（大滑移）。摩擦系数可以定义为接触压力、滑移速度或温度的函数。

对于金属成形模拟，剪切摩擦模型（τ = m*σ_y，m为摩擦因子）在高压接触下比库仑模型更合理——库仑模型在高压下会预测超过材料剪切强度的摩擦应力，这是非物理的。

实际案例：螺栓连接强度分析

以M10螺栓连接两块钢板为例，展示ABAQUS强度仿真流程：

模型设置

• 螺栓：M10×1.5，10.9级（σ_y=940 MPa，σ_u=1040 MPa）
• 钢板：Q345B（σ_y=345 MPa），厚12mm
• 预紧力：40 kN（螺栓拉伸截面上施加）
• 工作载荷：板面切向力20 kN

分析步

1. Step 1：预紧 — 用*PRE-TENSION SECTION施加40 kN预紧力
2. Step 2：锁定预紧 — 将螺栓长度固定在预紧后的长度
3. Step 3：施加载荷 — 板面施加20 kN切向力

接触设置

• 螺栓头-板面：面-面接触，罚函数法，μ=0.15
• 螺纹-板孔：面-面接触，罚函数法，μ=0.15
• 两板之间：面-面接触，罚函数法，μ=0.30

网格

• 螺栓：C3D8R（六面体减缩积分），螺纹区域加密至0.5mm
• 钢板：C3D8R，孔边加密至1mm，远离孔区域3mm
• 总单元数：约15万

结果

• 预紧后螺栓最大应力：680 MPa（第一牙根部），低于屈服940 MPa
• 施加切向力后：接触面发生微量滑移（0.02mm），螺栓应力升至720 MPa
• 板孔最大挤压应力：380 MPa（孔壁），低于Q345B许用挤压强度（1.5×345=517.5 MPa）
• 安全裕度：螺栓1.31，板孔1.36

关键发现

螺栓第一牙承担约35%的预紧载荷（螺纹载荷分布不均），这是疲劳失效的起始点。如果需要更精确的螺纹应力分布，应建真实螺纹几何而非简化为圆柱面。但真实螺纹模型的网格数量会增加到50万以上，计算时间约10倍。

强度评估准则的选择

强度仿真不是追求”最精确”的结果，而是在合理的计算成本下给出足够安全的工程判断。对于关键承力件，建议同时用线弹性校核和弹塑性极限分析，前者保证工作载荷下的安全裕度，后者评估极限承载能力。