1. 课程导入与工程问题界定
在工程设计与失效分析中，我们真正关心的往往不是“材料叫什么名字”，而是“在给定载荷、温度与加载路径下，它会产生多大的应力、应变与位移”。材料力学试验的任务，是用标准试样和受控加载，把这种力学响应测量出来；本构反演的任务，则是把离散的试验数据转化为可嵌入有限元或结构分析程序的本构模型参数。两者合在一起，构成从“物理现象”到“数学描述”的闭环。
举个具体例子。某汽车悬架控制臂采用 6061-T6 铝合金，设计时需要评估其在路面冲击下的塑性变形风险。试验室对圆棒试样做单轴拉伸，得到名义应力—应变曲线：弹性段斜率约 68.9 GPa，0.2% 屈服强度约 276 MPa，抗拉强度约 310 MPa，断后伸长率约 12%。但有限元软件需要的不是这条曲线本身，而是弹性模量、泊松比、屈服面参数、硬化曲线等模型输入。若直接把工程应力—应变曲线输入，忽略颈缩后的真实应力修正，就会在高应变区低估承载能力。因此，本构反演的第一步，是明确“试验测什么、模型要什么、两者如何对应”。
1. 试验设计——确定加载方式、试样几何、测量点与采样频率，使数据能覆盖目标工况；
2. 参数识别——从试验数据中提取弹性常数、屈服准则参数、硬化参数、黏弹性或损伤参数；
3. 模型验证——用独立试验（如不同应变率、不同应力三轴度）检验反演参数的外推能力。
典型工程场景远不止单轴拉伸。例如：
- 金属成形：需要大应变下的真实应力—应变曲线和各向异性屈服准则（如 Hill48、Barlat Yld2000），试验包括单轴拉伸、平面应变拉伸和液压胀形。
- 聚合物与橡胶：呈现超弹性、黏弹性与 Mullins 效应，需做单轴、等双轴、平面拉伸及循环加载，反演 Mooney-Rivlin、Ogden 或 Prony 级数参数。
- 焊接接头与增材制造：材料不均匀，需用微区压痕、数字图像相关（DIC）和反演方法识别局部本构参数。
- 高温结构：蠕变与应力松弛试验，反演 Norton-Bailey 或 Garofalo 模型参数，用于汽轮机、管道寿命评估。
这些场景的共同点是：试验数据有限、噪声不可避免、模型参数之间存在相关性。因此，本构反演不是“把曲线拟合一下”那么简单，而是一个带约束的优化问题。以单轴拉伸为例，若采用 Voce 硬化模型 σ = σy + Q(1 − e−bεp)，待定参数为初始屈服 σy、饱和增量 Q 和硬化速率 b。反演时通常以试验应力与模型应力的残差平方和为目标函数，用最小二乘或遗传算法求解。但若试验只到 5% 应变，外推到 20% 应变时误差可能超过 15%，所以必须用不同应力三轴度的试样（如缺口拉伸、剪切试样）交叉验证。
从学习路径看，本课程按“试验基础—本构理论—反演算法—工程应用”四段推进。前几章先建立应力、应变、弹性常数和试验标准的基本概念，随后进入塑性、超弹性、黏弹性等本构模型的形式与适用条件，再学习最小二乘、正则化、贝叶斯反演等参数识别方法，最后通过金属成形、橡胶密封、焊接接头和高温蠕变四个综合案例完成闭环训练。每一章都会给出可复现的示例数据、参数表和判断依据，而不是停留在概念介绍。
需要特别强调的是，本构反演的结果不是一组“真值”，而是一组在给定试验条件和模型假设下“足够好用”的参数。工程判断的关键在于：明确参数的适用范围、不确定度和外推风险。例如，用室温、准静态拉伸数据反演出的 Johnson-Cook 参数，直接用于高速冲击仿真时，必须补充应变率相关试验，否则模型给出的流动应力可能严重偏低。因此，本课程始终把“试验设计—参数识别—验证确认”作为一个整体来训练，而不是孤立地讲某一个模型或某一种算法。