2. 实验数据到仿真输入

材料模型的标定质量，首先取决于输入数据的质量。无论后续采用 *MAT_PLASTIC_KINEMATIC、*MAT_JOHNSON_COOK 还是 GISSMO 损伤模型，如果喂给参数识别算法的力—位移曲线本身含有零点漂移、高频振荡或不同通道之间的时间错位，那么拟合出的屈服强度、硬化模量乃至失效应变都会系统性偏离真实值。本节围绕拉伸、压缩、剪切、缺口拉伸和穿孔这几类常见试验，说明从原始采集数据到 LS-DYNA 可用的真应力—真应变（或等效塑性应变—损伤）曲线之间，需要完成哪些处理，以及每一步的判断依据。

一、原始数据的常见缺陷与清洗顺序

试验机采集到的通常是时间、载荷、位移（横梁位移或引伸计位移）以及可能的应变片、DIC 输出。原始数据进入标定流程前，建议按“去异常点 → 零点对齐 → 滤波 → 通道同步 → 曲线转换”的顺序处理，顺序颠倒会放大误差。例如先滤波再对齐零点，会把真实的初始非线性段一并抹掉。

去异常点主要针对采集丢帧、电磁干扰导致的单点跳变。判断方法不是简单看数值大小，而是看相邻点的斜率：若某点与前后点构成的一阶差分超过局部中位差分的 8～10 倍，可标记为可疑点并用线性插值替换。对于载荷通道，还需检查是否出现负值——拉伸试验中若载荷持续为负，说明夹持方向或传感器极性接反，属于必须回溯的硬件问题，不能靠算法修补。

二、零点对齐与刚度修正

零点对齐包含两个层面：载荷零点和位移零点。载荷零点应在试样正式受载前的稳定段取平均，而不是取第一个采样点。位移零点则更复杂：横梁位移包含试验机柔度，直接用于计算应变会低估刚度。工程上常用“刚度修正”处理：先做一组无试样或刚性试样的空载试验，得到机器柔度曲线 C_machine(F)，则试样真实位移 u_spec = u_crosshead − C_machine(F)·F。若没有空载试验，也可用初始线性段反推机器柔度，但精度较低，需在报告中注明。

三、滤波：截止频率与滤波器的选择

滤波的目的是抑制高频噪声，同时保留真实的材料响应特征。截止频率不能凭经验随意设定，应结合采样率和信号特征确定。一个实用做法是对载荷—时间信号做快速傅里叶变换，观察噪声主频与信号主频的分界。对于准静态试验（应变率约 10⁻³～10⁻¹ s⁻¹），截止频率通常取 50～200 Hz；对于中高应变率（10⁰～10² s⁻¹），截止频率需提高到 1～5 kHz，否则会削掉真实的应力波前缘。

滤波器类型推荐零相位 Butterworth 低通滤波（如 filtfilt 实现），避免普通 IIR 滤波引入的相位滞后。相位滞后会使载荷与位移通道之间产生虚假的时间差，在后续同步中表现为错误的初始斜率。滤波阶数一般取 2～4 阶，阶数过高会在截止频率附近产生振铃。

四、多通道同步

拉伸试验中载荷与引伸计、DIC 与力传感器往往由不同采集系统记录，采样率和触发时刻不同，必须做时间同步。同步的判据是：在弹性段，载荷—位移曲线应呈良好直线；若两通道存在时间偏移 Δt，曲线会在初始段出现“S”形或斜率异常。可通过互相关函数求取使弹性段线性度最优的 Δt，再对其中一个通道做时间平移。对于穿孔试验，冲头位移与载荷的同步尤为关键，Δt 达到 0.5 ms 就可能使峰值力对应的位移偏移数十微米。

五、曲线转换：工程量到真量

拉伸试验得到的是工程应力 σ_e = F/A₀ 和工程应变 ε_e = ΔL/L₀。进入 LS-DYNA 塑性材料卡前，需转换为真应力 σ_t = σ_e(1+ε_e) 和真应变 ε_t = ln(1+ε_e)。在颈缩之前，这一转换是严格的；颈缩之后，变形集中在颈缩区，均匀应变假设失效，需借助 Bridgman 修正或 DIC 局部应变测量。压缩试验因试样与压板摩擦导致鼓肚，需做摩擦修正或采用 Rastegaev 试样。剪切试验（如双缺口剪切）通常直接输出等效剪应力—剪应变，转换为等效应力—等效塑性应变时需乘以 √3 关系。

六、一个可复用的处理记录表

为了让标定过程可追溯，建议对每条曲线保留一份处理记录。以下字段可直接用于实验数据管理，示例数值为教学假设，不代表真实试验结果。