2. Abaqus/CAE建模实战

本章以一个可运行的非线性静力模型为主线，把 Abaqus/CAE 的建模链条完整走一遍。为了让过程可复现，我们选用一个经典教学对象：带中心圆孔的矩形薄板，尺寸 100 mm × 50 mm × 2 mm，中心圆孔直径 10 mm，材料为 Q235 低碳钢，左端固支，右端施加 2 mm 的位移载荷。这个模型几何简单、边界清晰，却能真实触发孔边应力集中、局部塑性以及几何非线性效应，非常适合作为非线性静力分析的入门载体。

需要先明确一个贯穿全章的约定：坐标系采用全局直角坐标系，X 轴沿板长方向（100 mm），Y 轴沿板宽方向（50 mm），Z 轴为厚度方向。位移载荷沿 +X 方向施加，因此下文所有“正向”“负向”均以此为准。切换参照系时必须重新说明，本章不再重复。

一、部件与装配：几何是分析的起点

在 Part 模块中创建 2D Planar、Deformable 部件，先绘制 100 × 50 的矩形，再在中心 (50, 25) 处绘制直径 10 的圆，用 Cut Extrude 切除。这里有一个容易被忽略的细节：圆孔边缘应保留足够的网格密度，因此建议在草图阶段就把孔边分成 16 段，为后续结构化网格打基础。部件创建完成后，进入 Assembly 模块，选择 Dependent (mesh on part) 实例化。对于单部件模型，装配看似多余，但它决定了后续载荷、边界条件和网格都作用在实例上，而不是部件本身，这是 Abaqus 数据层级的关键。

二、材料与截面：非线性从本构开始

非线性静力分析的“非线性”首先来自材料。在 Property 模块中创建材料 Q235，密度 7.85e-9 t/mm³（Abaqus 无单位制，需自行保持一致），弹性模量 210000 MPa，泊松比 0.3。塑性部分采用等向强化，输入真实应力—塑性应变数据。一个合理的教学示例如下：

Yield Stress (MPa)   Plastic Strain
235                  0.0
300                  0.02
360                  0.05
400                  0.10

注意：Abaqus 塑性数据要求的是真实应力与塑性应变，而非工程应力—应变。若只有工程数据，需先转换。截面类型选择 Homogeneous Solid、Section 为 Solid，厚度方向由几何决定。将截面指派给整个部件后，务必用 Assign Section 确认颜色变化，这是判断指派是否成功的直观依据。

三、分析步：把“非线性”显式打开

进入 Step 模块，创建 Static, General 分析步，时间总长设为 1，初始增量 0.05，最小 1e-5，最大 0.1。关键设置有三处：第一，Nlgeom 必须设为 ON，否则大变形效应被忽略；第二，增量步采用 Automatic，让求解器根据收敛情况自适应；第三，在 Incrementation 中打开 Adaptive Automatic Stabilization 作为收敛困难时的兜底。若模型出现局部失稳或接触，可进一步引入阻尼因子，但本章模型通常不需要。

四、载荷与边界条件：约束要“够用而不多余”

在 Load 模块中，对左端面施加 Encastre（U1=U2=U3=UR1=UR2=UR3=0）固支约束；对右端面施加位移边界条件 U1=2 mm，其余方向自由。这里要特别提醒：位移载荷比力载荷更容易收敛，尤其在接近塑性屈服时，力控制容易导致迭代发散。若必须用力控制，建议配合 Riks 弧长法，但本章不展开。

边界条件的施加对象是面还是边，会直接影响结果。对于 2D 平面应力模型，左端应约束整条边，而不是单个节点，否则会产生不真实的应力集中。

五、网格：质量决定收敛与精度

在 Mesh 模块中，全局种子设为 2 mm，孔边局部种子加密至 0.5 mm，单元类型选择 CPE4R（四节点双线性平面应力减缩积分）。减缩积分能缓解剪切自锁，但需注意沙漏控制，Abaqus 默认的 Enhanced 沙漏控制通常足够。网格划分后，用 Verify 检查单元质量，重点关注最小内角、长宽比和雅可比。孔边单元若出现严重畸变，应回到几何阶段调整分割策略。

六、作业提交与结果判读

在 Job 模块创建作业，提交前建议先做 Data Check，确认无未定义截面、无孤立节点、无过约束。提交后监控 .sta 文件，观察增量步是否顺利推进。若出现“Too many attempts made for this increment”，优先检查材料塑性数据是否合理、边界条件是否过约束、网格是否畸变。

计算完成后，进入 Visualization，查看 Mises 应力云图。对于本模型，孔边应力集中系数理论值约为 3，即孔边最大应力约为远场应力的 3 倍。若远场应力为 100 MPa，孔边应接近 300 MPa，已超过屈服强度 235 MPa，说明孔边进入塑性。此时应重点观察：塑性区是否局限于孔边、载荷—位移曲线是否出现刚度退化、最大塑性应变是否在可接受范围。若曲线在 2 mm 位移时仍未收敛，可减小增量步或引入稳定化阻尼。