1. 课程导入与工程问题界定
在材料与结构仿真中，几何模型决定“形状”，网格决定“离散精度”，边界条件决定“外部约束”，而本构模型决定材料如何把变形转化为应力。这句话看起来抽象，但它直接决定了仿真结果是否可信。一个悬臂梁的有限元模型即使网格再细、载荷再准，如果本构模型选错，位移和应力分布仍可能偏离真实情况数倍。因此，本构模型不是仿真流程中的“可选项”，而是连接力学输入与工程判断的核心环节。
从工程视角看，本构模型至少承担三项职责：第一，描述材料在不同应力状态下的响应规律，例如线弹性、塑性屈服、硬化、蠕变、损伤与失效；第二，为求解器提供应力—应变增量关系，使每个积分点能够根据当前应变增量更新应力；第三，为工程判断提供可解释的参数体系，例如弹性模量、泊松比、屈服强度、硬化指数、蠕变参数等。换句话说，本构模型既是数学表达，也是工程假设的集合。
选型失误的工程后果通常不会以“报错”形式出现，而是以“结果看似合理但实际错误”的形式潜伏。常见后果包括：用线弹性模型分析金属塑性成形，导致应力被严重高估、回弹预测失真；用理想弹塑性模型分析循环载荷，忽略随动硬化，导致棘轮效应或安定性判断错误；用室温参数分析高温蠕变，忽略时间相关变形，导致长期位移和松弛载荷被低估；用各向同性损伤模型分析纤维增强复合材料，忽略方向相关失效，导致破坏位置和承载能力判断偏离。更隐蔽的是，错误模型可能给出“漂亮”的云图，让工程师误以为问题已经解决。
下面这张图把本构模型在仿真链路中的位置、输入输出以及选型失误的传导路径放在一起，便于建立整体认识。
从图中可以看到，本构模型位于“工程输入”和“求解器”之间，是应力更新的核心。选型失误不会在输入阶段暴露，而是沿着求解链路传导到结果，最终表现为工程判断失效。因此，本课程的任务不是简单罗列模型公式，而是建立一套可操作的选型与验证方法。
为了把“选型”从经验判断变成可复现的工程决策，本课程建议采用如下任务地图：第一步，界定问题——明确材料类型、载荷类型、温度、时间尺度、失效模式和精度要求；第二步，建立模型谱系——理解线弹性、弹塑性、蠕变、损伤、超弹性、粘弹性等模型的适用边界；第三步，参数标定——从试验数据中提取模型参数，并检查参数是否落在合理范围；第四步，选型决策——根据问题特征选择候选模型，并用简单算例做敏感性对比；第五步，验证与复盘——用独立试验或解析解验证，记录选型依据和误差来源。
下面给出一个可直接用于工程项目的本构模型选型检查表示例。字段和判断依据均围绕材料与结构仿真场景设计，可根据具体项目增删。
| 检查项 | 需要确认的内容 | 示例判断 |
|---|---|---|
| 材料类型 | 金属、聚合物、复合材料、混凝土、岩土等 | 低碳钢，室温，准静态 |
| 主导变形机制 | 弹性、塑性、蠕变、粘弹性、损伤、断裂 | 塑性主导，需考虑硬化 |
| 载荷特征 | 单调/循环、速率相关、温度相关、多轴 | 单调拉伸，应变率约 1e-3 /s |
| 时间尺度 | 短时、长时、松弛、蠕变 | 短时，可忽略蠕变 |
| 失效模式 | 屈服、颈缩、断裂、疲劳、失稳 | 屈服后塑性流动，关注颈缩前响应 |
| 参数来源 | 标准试验、文献、反演、经验值 | 单轴拉伸试验 + 文献硬化参数 |
| 验证方式 | 解析解、独立试验、网格敏感性、参数敏感性 | 与拉伸曲线对比，误差控制在 5% 以内 |
以低碳钢单向拉伸为例，若应变率约 1e-3 /s、温度约 20 ℃、关注屈服后到颈缩前的应力—应变响应，则线弹性模型只能描述初始阶段，无法给出屈服后的塑性流动；理想弹塑性模型可以给出屈服平台，但无法描述硬化；各向同性硬化弹塑性模型（如 Voce 或 Swift 硬化）能够较好拟合屈服后的上升段。此时选型的关键不是“哪个模型更复杂”，而是“硬化参数是否可从试验曲线稳定提取”。如果试验数据只覆盖 2% 应变，却外推到 20% 应变，参数外推风险就会成为主要误差来源。
本课程后续章节将围绕上述任务地图展开：先建立模型谱系，再逐类深入弹性、塑性、蠕变、损伤、超弹性、粘弹性等模型，随后进入参数标定、选型决策、验证与复盘。学习路径不是线性背诵公式，而是“问题—模型—参数—验证”的循环迭代。每学完一类模型，都应回到本节检查表，确认它适合解决哪类工程问题、需要哪些输入、结果如何判断。
最后需要强调：本构模型选型没有唯一正确答案，只有与工程问题匹配的合理选择。一个在金属成形中表现良好的模型，放到橡胶密封件分析中可能完全失效；一个在室温下足够精确的模型，放到高温蠕变场景中可能严重低估变形。因此，本课程的目标不是让你记住所有模型，而是让你在面对新问题时，能够快速界定问题、筛选候选模型、评估参数风险，并给出可验证的选型结论。