2. 实验数据获取与预处理

材料参数标定的精度上限，几乎完全由实验数据的质量决定。无论后续用最小二乘、优化算法还是贝叶斯反演，输入曲线中的噪声、漂移、单位错误或坐标错位都会被参数拟合过程放大，最终表现为“仿真能对上一条曲线、换一条工况就崩”。本节围绕拉伸、压缩、剪切和动态力学分析（DMA）四类常见实验，说明数据采集的关键控制点，以及去噪、平滑、滤波、单位统一和曲线对齐的具体做法。

一、四类实验的采集要点

单轴拉伸是获取弹性模量、屈服强度和塑性硬化段的主力实验。采集时需注意三点：一是引伸计标距必须与后续仿真中用于计算应变的标距一致，若实验用 25 mm 引伸计而仿真按 50 mm 标距换算，弹性段斜率会直接偏差；二是采样率要覆盖断裂瞬间，金属拉伸建议不低于 50 Hz，高分子材料因断裂伸长率大、软化明显，建议 100 Hz 以上；三是必须记录载荷—位移原始数据，而非只保存工程应力—应变，因为真实应力—应变需要由原始截面积和瞬时长度重新计算。

单轴压缩常用于泡沫、橡胶和混凝土类材料。压缩实验最大的干扰来自压板与试样的摩擦约束，它会造成“鼓肚”并虚高流动应力。采集时应记录压板位移与载荷，同时标注是否使用了润滑（如聚四氟乙烯薄膜+润滑脂）。对于橡胶类超弹性材料，还需注意加载—卸载循环中的 Mullins 效应，通常取第 2 或第 3 次循环的数据用于标定，以排除首次软化的历史影响。

剪切实验（含纯剪切、轨道剪切、薄壁圆筒扭转）用于标定剪切模量和剪切屈服。其采集难点在于剪切应变往往由夹具位移换算，包含夹具柔度。建议先做一次空载或刚性试样标定，测出系统柔度 C_m，再从总位移中扣除：γ = (δ_总 − F·C_m)/h，其中 h 为试样有效剪切高度。

DMA用于获取黏弹性参数（储能模量 E′、损耗模量 E″、损耗因子 tanδ）。采集要点是：频率扫描、温度扫描和应变扫描三种模式目的不同，频率扫描用于构建主曲线，温度扫描用于时温等效移位；必须确认处于线性黏弹区（应变扫描中 E′ 基本不随应变变化），否则大应变下的非线性会污染主曲线。DMA 输出的原始量是动态刚度与相位角，需按试样几何换算为模量。

二、去噪与平滑：先判断噪声来源

噪声不是单一概念。常见有三类：高频随机噪声（电子干扰、振动），表现为曲线毛刺；脉冲异常点（信号丢包、夹具打滑），表现为孤立尖峰；低频漂移（温漂、零点漂移），表现为整条曲线缓慢偏移。处理顺序应是先去脉冲、再去高频、最后处理漂移，否则平滑会把尖峰“抹开”成一段虚假的缓变段。

脉冲点可用 3σ 准则或中值滤波识别：对每个点计算其与局部中值的偏差，超过 3 倍局部标准差即标记为异常并剔除。高频噪声推荐 Savitzky-Golay 平滑，它在拟合多项式的同时保留峰值高度和位置，比移动平均更适合屈服点、tanδ 峰这类特征。窗口长度取数据点数的 1%～3%，多项式阶数取 2～3 阶；窗口过大将压低屈服峰，导致屈服强度被低估。

低频漂移可用多项式拟合基线后相减，或对弹性段做线性回归求零点偏移。例如拉伸曲线在应力为零时应变读数不为零，说明引伸计存在初始偏置，应整体平移使弹性段起点归零。

三、滤波与单位统一

滤波主要用于 DMA 和动态信号。DMA 的相位角对噪声敏感，建议对动态刚度和相位角分别做低通滤波，截止频率取测试频率的 3～5 倍，既能保留响应又不引入相位滞后。若使用带通滤波，中心频率应锁定测试频率，带宽不宜过窄，否则会削掉非线性谐波成分。

单位统一是最容易被忽视却最致命的环节。Abaqus 没有内置单位制，必须自洽。推荐统一采用 mm-N-MPa-s 体系：长度 mm、力 N、应力 MPa、密度 t/mm³、时间 s。常见错误是把实验的 MPa 与仿真的 Pa 混用，或把密度写成 kg/m³（应为 7.8e-9 t/mm³ 量级）。建议在数据表头显式标注单位，并在导入前做一次量纲检查：弹性模量应在 10²～10⁵ MPa 量级，屈服强度应低于模量两个数量级左右，超出即需复查。

四、曲线对齐与多工况归一

曲线对齐包含三个层面：零点对齐（消除初始偏置）、标距对齐（引伸计标距与仿真标距一致）、工况对齐（不同应变率、温度下的曲线按同一参考状态归一）。对于率相关材料，常以参考应变率下的曲线为基准，其余曲线通过移位因子平移，这与 DMA 的时温等效主曲线构建思路一致。

多工况数据在标定前应统一到同一应变定义。工程应变与真实应变的换算关系为 ε_true = ln(1+ε_eng)，σ_true = σ_eng(1+ε_eng)，该式在颈缩前有效；颈缩后需用 Bridgman 修正或直接截断至颈缩点，避免用错误的外推段标定硬化参数。