4. 测试与仿真坐标系对齐

在应变测试与 FEA 对标中,最容易造成“数据对不上”的往往不是测量精度或网格质量,而是坐标系没有对齐。应变是张量,它的数值依赖于方向:同一个物理变形状态,在全局坐标系下可能表现为以 εx 为主,在局部坐标系下却可能表现为以 ε1 为主。如果试验和仿真各自使用不同的参考方向,即使物理状态完全一致,提取出来的应变分量也会不同,对标就失去了意义。

全局坐标系与局部坐标系的分工

全局坐标系是结构建模和装配的基准,通常记为 X-Y-Z,用于定义零件位置、约束和整体载荷方向。局部坐标系则附着在具体对象上,例如应变片粘贴位置、焊缝方向、复合材料铺层方向或某个关注截面。局部坐标系通常记为 1-2-3 或 x'-y'-z',其中 1 轴沿应变片敏感栅方向,2 轴在片基平面内与 1 轴垂直,3 轴为片基法向。

两者之间的转换由方向余弦矩阵完成。若局部坐标系三个轴在全局坐标系中的单位向量分别为 e1、e2、e3,则全局应变张量 ε 与局部应变张量 ε' 满足 ε' = Q ε QT,其中 Q 的行向量即为 e1、e2、e3。工程应变分量在转换时还需注意剪应变的因子 2 问题:张量剪应变 γ12 与工程剪应变 ε12 的关系为 γ12 = 2ε12,混用会导致剪应变对标出现成倍偏差。

应变片粘贴方向的定义与记录

应变片只能测量其敏感栅方向上的正应变。对于单轴片,粘贴方向就是 1 轴方向;对于 0°-45°-90° 三向片,三个敏感栅分别沿 0°、45°、90° 方向,可解出平面内的 εx、εy、γxy,再计算主应变 ε1、ε2 和主方向角 θp。主方向角通常以全局 X 轴为基准,逆时针为正,公式为 θp = 0.5·arctan(γxy / (εx − εy))。

实际操作中,粘贴方向必须可追溯。建议在测试记录中至少保留以下字段:

  • 测点编号与位置描述(如“左侧翼肋根部上表面,距前缘 120 mm”);
  • 应变片型号、栅长、敏感栅数量;
  • 参考基准(全局 X 轴或某条结构边线);
  • 粘贴角度 α(敏感栅 1 轴相对基准的夹角,逆时针为正);
  • 片基法向方向,用于判断是否处于平面应力状态;
  • 温度补偿方式与导线走向。

例如某测点记录为:基准为全局 X 轴,粘贴角 α = 30°,使用 0°-45°-90° 三向片。若实测得到 ε0 = 420 με、ε45 = 260 με、ε90 = −80 με,则可解出 εx = 420 με、εy = −80 με、γxy = 2ε45 − εx − εy = 2×260 − 420 + 80 = 180 με。主应变为 ε1,2 = (εx + εy)/2 ± 0.5·√((εx − εy)² + γxy²) = 170 ± 0.5·√(500² + 180²) ≈ 170 ± 265.7,即 ε1 ≈ 435.7 με、ε2 ≈ −95.7 με,主方向 θp = 0.5·arctan(180/500) ≈ 9.9°。这些数值是示例计算,用于说明方向转换过程,并非实测结果。

FEA 中的材料方向与结果坐标系

在 FEA 中,材料方向决定了各向异性材料(如复合材料、木材、增材制造件)的刚度分布。以复合材料层合板为例,铺层角是相对于材料 1 轴定义的,而材料 1 轴又通过局部坐标系或方向矢量指定。如果材料方向与试验粘贴方向不一致,仿真得到的应变场在物理上就对应了不同的铺层配置,对标必然失败。

结果坐标系则决定后处理中输出的应变分量。多数求解器默认输出全局坐标系下的分量,但也可以指定局部坐标系输出。对标时应统一到同一坐标系:要么把试验结果转换到全局,要么把仿真结果转换到应变片局部方向。推荐做法是在应变片局部坐标系下比较,因为应变片本身只敏感一个方向,局部比较能直接对应测量物理量。

常见错误:仿真结果直接读取全局 εx,而试验片粘贴角为 30°,两者方向不同却直接对比,导致误差被误判为模型误差。正确做法是先将仿真应变张量按粘贴角旋转,再与实测值比较。

一致性处理流程

为保证试验与仿真坐标系一致,建议按以下顺序处理:

  1. 确定全局坐标系,并在试验记录和 FEA 模型中明确标注其原点与轴向;
  2. 为每个测点建立局部坐标系,记录 1 轴方向(粘贴方向)和 3 轴方向(法向);
  3. 在 FEA 中为对应位置定义相同的局部坐标系,或记录材料方向与全局的夹角;
  4. 提取仿真应变张量,按局部坐标系方向余弦进行旋转,得到 ε11、ε22、γ12;
  5. 将旋转后的仿真值与试验解出的 εx、εy、γxy 或主应变进行对比;
  6. 若偏差较大,检查角度符号、剪应变因子 2、单位(με 与 ε)以及温度补偿是否一致。

下面用一张关系图概括试验侧与仿真侧在坐标系对齐中的关键对象和转换路径。

试验与仿真坐标系对齐关系 全局坐标系 X-Y-Z 结构基准 / 载荷方向 试验侧 应变片粘贴方向 α 敏感栅 1 轴 / 片基法向 3 轴 仿真侧 材料方向 / 铺层角 结果坐标系选择 局部坐标系 1-2-3 1 轴 = 粘贴方向 / 材料主方向 方向余弦矩阵 Q 坐标转换 ε' = Q · ε · Qᵀ 注意 γ = 2ε 的因子关系 同一局部坐标系下对比 ε11、ε22、γ12 或 ε1、ε2、θp 试验与仿真必须统一到同一局部坐标系后再比较应变分量

实操中的检查要点

在完成一次对标前,建议逐项核对以下内容:全局坐标系原点与轴向是否在试验报告和 FEA 模型中一致;每个测点的粘贴角是否以同一基准定义,逆时针为正还是顺时针为正;FEA 中材料 1 轴是否与试验粘贴方向重合;结果提取时是否已旋转到局部坐标系;剪应变使用的是张量剪应变还是工程剪应变;应变单位是否统一为 με;温度补偿是否在两侧都考虑。任何一项不一致,都可能让对标结果出现系统性偏差。

当试验与仿真在同一局部坐标系下比较时,若 ε11 偏差在 5% 以内、γ12 偏差在 10% 以内,通常可认为模型在弹性范围内具有较好的一致性。若偏差显著,应优先排查坐标系和方向定义,再检查边界条件、材料参数和网格收敛性。坐标系对齐不是后处理的小技巧,而是对标工作的前提条件。