3. 试验机柔度标定与位移修正

为什么必须做柔度修正

单轴压缩试验中，引伸计或压板位移传感器测到的“位移”并不等于试样本身的轴向变形。加载链上还有试验机框架、上压头、下压头、球铰座等部件，它们同样在受力时发生弹性压缩。设传感器测得总位移为 ΔL_total，则：

ΔL_total = ΔL_specimen + ΔL_machine，其中 ΔL_machine = C_m · P，C_m 为机器柔度（mm/kN），P 为载荷。若不扣除 ΔL_machine，试样刚度会被系统性低估，弹性模量随之偏低。

标定的标准做法是空载压板试验：不放试样，让上下压头直接对压，记录载荷—位移曲线。此时试样变形为零，全部位移都是机器变形，曲线斜率即为机器柔度 C_m。

标定数据处理步骤

1. 截取有效段：空载曲线起始段存在压头贴合、间隙消除过程，需从曲线明显转为线性处开始取点；上限取到额定载荷的 60%~80%，避免压头局部屈服。
2. 线性拟合：对有效段做 ΔL = C_m · P + b 最小二乘拟合，斜率 C_m 即柔度，截距 b 反映初始间隙。
3. 逐点修正：对试样试验的每个采样点，ΔL_corr = ΔL_total − C_m · P。
4. 验证：比较修正前后弹性段斜率（刚度），修正后刚度应显著提高且更接近材料真实值。

关键实现（Python 示例）

下面代码完成三件事：对空载数据做线性拟合求 C_m，对试样数据逐点修正，并计算修正前后弹性段刚度。示例数据为教学构造，实际使用时替换为试验机导出文件。

import numpy as np

# ---------- 1. 空载压板试验数据（示例：P kN, dL mm） ----------
P_blank = np.array([0.0, 5.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0])
dL_blank = np.array([0.031, 0.094, 0.156, 0.219, 0.281, 0.344, 0.406, 0.469, 0.531])

# 有效段：取 P >= 10 kN 的线性段（避开贴合段）
mask = P_blank >= 10.0
Cm, b = np.polyfit(P_blank[mask], dL_blank[mask], 1)
print(f"机器柔度 Cm = {Cm:.5f} mm/kN, 截距 b = {b:.4f} mm")

# ---------- 2. 试样试验数据（示例） ----------
P_s = np.array([0, 5, 10, 15, 20, 25, 30, 35, 40], dtype=float)
dL_total = np.array([0.000, 0.062, 0.125, 0.187, 0.250, 0.312, 0.375, 0.437, 0.500])

# 逐点扣除机器变形
dL_machine = Cm * P_s
dL_corr = dL_total - dL_machine

# ---------- 3. 弹性段刚度对比（取 10~30 kN 段） ----------
seg = (P_s >= 10) & (P_s <= 30)
k_raw  = np.polyfit(dL_total[seg], P_s[seg], 1)[0]
k_corr = np.polyfit(dL_corr[seg],  P_s[seg], 1)[0]
print(f"修正前刚度 k_raw  = {k_raw:.2f} kN/mm")
print(f"修正后刚度 k_corr = {k_corr:.2f} kN/mm")
print(f"刚度提升比例 = {(k_corr/k_raw - 1)*100:.1f}%")

结果判读与验证要点