1. 课程导入与非线性问题概览

在结构有限元分析中，线性分析建立在两个关键假设之上：材料本构关系是线性的（应力与应变成正比），且变形足够小以至于几何方程可以线性化。当结构中的应力水平超过材料屈服点、橡胶类材料经历大变形、聚合物在时间载荷下产生蠕变，或者裂纹在循环载荷下逐步扩展时，这些假设便不再成立。此时若仍用线性模型求解，得到的位移、应变和应力分布可能与真实物理行为产生显著偏差，甚至给出完全错误的承载能力判断。Abaqus 作为通用有限元平台，其核心优势之一正是对材料非线性的系统支持——它提供了从金属塑性、橡胶超弹性到粘弹性、损伤与断裂的完整材料库，并通过 Newton-Raphson 迭代与一致切线刚度矩阵保证非线性方程组的收敛求解。

材料非线性在结构仿真中的定位，可以从“输入—求解—输出”三个环节来理解。输入环节，用户需要选择恰当的本构模型并给出材料参数，例如金属塑性需要弹性模量、泊松比、屈服应力与硬化曲线；超弹性需要单轴、双轴或平面拉伸试验数据拟合得到的应变能势参数；粘弹性需要松弛模量或蠕变柔量。求解环节，Abaqus 在每个增量步内根据当前应变状态计算应力与切线刚度，若残差力不满足收敛容差则继续迭代，直到平衡或达到最大迭代次数。输出环节，除了常规的位移、应变、应力，还需要关注等效塑性应变、损伤变量、能量耗散等内变量，它们直接反映非线性效应的程度与位置。

常见的材料非线性类型可按物理机制分为几大类。塑性描述金属等材料在屈服后的不可逆变形，Abaqus 中常用 *PLASTIC 配合 *ELASTIC 定义，硬化方式包括各向同性硬化、随动硬化和混合硬化。超弹性用于橡胶、生物软组织等，Abaqus 提供 Neo-Hookean、Mooney-Rivlin、Ogden、Yeoh 等模型，参数通常由试验数据拟合获得。粘弹性描述时间或频率相关的耗散行为，可用 Prony 级数表示松弛模量，常用于聚合物、沥青和生物材料。损伤类模型则描述刚度退化与破坏，包括延性损伤、剪切损伤、Hashin 复合材料损伤以及基于断裂力学的 cohesive 单元与 XFEM。此外，还有多孔材料、土壤塑性、形状记忆合金等专用本构，它们共同构成了材料非线性的完整谱系。

学习路径上，建议先建立“弹性—塑性—超弹性—粘弹性—损伤”的递进认知。第一步，用简单单轴拉伸模型验证线弹性响应，确认单位制、边界条件和输出变量设置正确。第二步，引入金属塑性，在 Abaqus/Standard 中定义 *PLASTIC 硬化曲线，观察卸载后的残余应变与等效塑性应变分布。第三步，切换到超弹性模型，用单轴试验数据拟合 Mooney-Rivlin 参数，对比不同应变能势的预测能力。第四步，加入时间维度，用 *VISCOELASTIC 定义 Prony 级数，进行应力松弛或蠕变分析。第五步，引入损伤与断裂，从延性损伤起始准则到单元删除，再到 cohesive 界面，逐步理解刚度退化和能量耗散。每一步都应配合收敛性检查，因为材料非线性往往伴随迭代困难，需要调整增量步长、稳定化参数或求解器设置。

软件环境准备方面，Abaqus 提供 Standard 与 Explicit 两个主求解器。Standard 适合静态、准静态和低频动态问题，采用隐式积分与 Newton-Raphson 迭代，对材料非线性的处理成熟稳定；Explicit 适合高速冲击、接触复杂和高度非线性问题，采用显式时间积分，无需迭代但需要控制时间步长与质量缩放。材料非线性的大部分本构在两者中均可使用，但损伤演化、单元删除和接触行为在 Explicit 中往往更易收敛。建议安装 Abaqus 2023 或更新版本，并配置好 Python 脚本环境，便于参数化建模与批量后处理。此外，准备一个简单的单轴拉伸 INP 文件作为练习模板，包含材料定义、边界条件、加载幅值和历史输出请求，可以快速验证本构设置是否正确。

判断材料非线性分析是否可信，不能只看“计算是否完成”。需要检查三个层面：一是收敛性，增量步是否顺利推进，残差力是否在容差内；二是物理合理性，应力是否超过材料强度、等效塑性应变是否集中在预期位置、损伤变量是否在合理范围；三是能量平衡，塑性耗散、粘性耗散和损伤耗散是否与外部做功匹配。Abaqus 的 .sta、.msg 和 .dat 文件提供了收敛历史、警告与错误信息，后处理中的 ALLPD、ALLVD、ALLDMD 等能量变量可用于交叉验证。只有这三层都通过，材料非线性仿真结果才具备工程参考价值。