📚 课程目录第2章试验数据获取与预处理

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

材料参数标定的精度上限，几乎完全由试验数据的质量决定。无论后续在 Ansys Mechanical 中使用曲线拟合、参数反演还是优化算法，输入数据中的噪声、异常点和坐标系定义错误都会被放大成参数偏差。本节围绕“面向标定需求”这一核心，讨论如何设计拉伸、压缩、剪切、三点弯曲四类基础试验，以及拿到原始数据后如何完成去噪、平滑、滤波、异常点剔除和工程应力-应变到真实应力-应变的转换。

一、试验类型与标定目标的对应关系

不同本构模型需要不同应力状态的试验数据来约束参数。设计试验时，应先明确“这份数据要标定哪些参数”，再决定加载方式、试样尺寸和采集频率。

试验类型	典型试样	主要标定参数	关键采集量
单轴拉伸	哑铃形/狗骨形板材	弹性模量 E、泊松比 ν、屈服强度、硬化曲线	载荷 F、标距段伸长 ΔL、横向应变
单轴压缩	圆柱或立方体	压缩屈服、体积模量、Drucker-Prager/Cap 参数	载荷 F、轴向位移、端面摩擦状态
纯剪切/扭转	薄壁圆管或平板剪切	剪切模量 G、剪切屈服、各向异性参数	扭矩 T、扭转角 φ、剪切应变
三点弯曲	矩形梁	弯曲强度、层间/界面参数、损伤起始	跨中载荷 F、挠度 δ、跨距 L

以金属塑性标定为例，单轴拉伸提供等效应力-应变主曲线，压缩试验补充大应变段（拉伸在颈缩后难以继续提供均匀应变），剪切试验用于校核屈服面形状。若目标是标定超弹性或泡沫材料，压缩和三点弯曲的数据权重往往更高。试验设计阶段就应记录试样几何、加载速率、环境温度和夹持方式，因为这些信息直接影响后续参数的可移植性。

二、原始数据的去噪、平滑与滤波

试验机采集的载荷和位移信号通常包含三类干扰：高频电噪声、机械振动引起的周期性波动、以及夹持间隙或引伸计滑移造成的阶跃。处理顺序建议为：先剔除明显异常点，再做平滑或滤波，最后进行坐标转换。

异常点剔除可采用基于物理约束的判据，而不是单纯依赖统计。例如拉伸曲线中，若某点载荷突然下降超过前一点的 15%，且位移没有对应突增，则很可能是信号丢帧或夹持打滑。对于三点弯曲，若挠度增加而载荷不增反降，需检查是否已进入损伤软化段，不能一律当作异常点删除。

平滑方法的选择取决于噪声特征。移动平均适合低频漂移，但会削弱屈服拐点；Savitzky-Golay 滤波在保留峰值和拐点方面更好，常用窗口 11～21 点、多项式阶数 2～3。低通 Butterworth 滤波适合去除高频振动，截止频率可取采集频率的 1/10～1/5。以下是一个典型的 Python 预处理片段，用于演示流程，参数为示例值：

import numpy as np
from scipy.signal import savgol_filter, butter, filtfilt

# 假设已读入 raw_F, raw_dL，采样频率 fs = 100 Hz
# 1) 异常点剔除：载荷突降且位移未突增
mask = np.ones_like(raw_F, dtype=bool)
for i in range(1, len(raw_F)):
    dF = raw_F[i] - raw_F[i-1]
    dL = raw_dL[i] - raw_dL[i-1]
    if dF < -0.15 * raw_F[i-1] and dL < 0.02 * raw_dL[i-1]:
        mask[i] = False
F = raw_F[mask]; dL = raw_dL[mask]

# 2) Savitzky-Golay 平滑，窗口 15，二阶
F_s = savgol_filter(F, window_length=15, polyorder=2)

# 3) 低通滤波，截止频率 5 Hz
b, a = butter(2, 5/(0.5*100), btype='low')
F_f = filtfilt(b, a, F_s)

平滑后的曲线应回看是否出现“过平滑”：屈服平台被抹平、颈缩点被抬高、卸载段出现虚假回弹。若出现这些现象，应减小窗口或降低滤波强度。滤波只改变数据形态，不改变物理事实，任何处理都应在报告中保留原始数据与处理参数的对应关系。

三、工程应力-应变转真实应力-应变

Ansys Mechanical 中的塑性、超弹性和损伤模型大多以真实应力-真实应变或真实应力-塑性应变作为输入。工程应力 σ_eng = F/A₀，工程应变 ε_eng = ΔL/L₀。在均匀变形阶段，真实应力 σ_true = σ_eng(1+ε_eng)，真实应变 ε_true = ln(1+ε_eng)。进入颈缩后，上述公式不再成立，需要借助 Bridgman 修正或直接截断到颈缩点。

对于压缩试验，若端面摩擦较大，试样会出现鼓肚，实测应力偏高。此时应通过润滑或修正因子降低摩擦影响，或在标定中只使用小应变段。剪切试验中，若采用薄壁圆管扭转，剪切应力 τ = T·r/J，剪切应变 γ = φ·r/L，其中 r 为平均半径，J 为极惯性矩，L 为标距。三点弯曲则通常通过梁理论反算弯曲应力，但更推荐直接用载荷-挠度曲线参与反演标定，避免引入额外假设。

转换前务必统一单位：载荷 N、位移 mm、应力 MPa、应变无量纲。若试验机输出为 kN 和 μm，需先换算，否则后续参数会整体偏移 10³ 量级。

四、面向标定的数据交付检查

完成预处理后，交付给标定环节的数据应满足以下条件：曲线单调段无负斜率异常；屈服点附近保留足够数据密度；真实应力-应变曲线在颈缩点前连续；压缩和剪切数据已注明摩擦或夹持修正方式；每条曲线附带试样编号、试验类型、加载速率和预处理参数。建议将原始数据、处理脚本和最终曲线一并归档，便于后续复现和参数追溯。

需要强调的是，预处理不是“把曲线修好看”，而是让数据更接近材料真实响应。任何平滑、滤波和修正都应在报告中记录方法和参数，并在标定结果中评估其敏感性。若某组参数对预处理方式高度敏感，说明试验数据本身的信息量不足，应补充试验或改进采集方案，而不是继续调整滤波窗口。