2. 装配误差的数学描述与建模基础

在精密机构装配中,误差不是单一数值,而是具有方向、作用点与耦合关系的空间量。要把“某个零件偏了 5 μm”转化为可计算、可传递、可优化的模型,必须建立统一的数学描述。本节围绕标量、矢量、矩阵三种表示层次,重点讲解小位移旋量(Small Displacement Screw, SDS)与齐次变换矩阵,并给出误差旋量到雅可比矩阵的映射关系,最后明确建模假设与适用边界。

一、误差描述的三种层次

标量表示是最简单的形式,例如“孔心距偏差 +8 μm”“端面平行度 3 μm”。它只给出大小,不携带方向与作用点信息,适合单自由度、方向已固定的场合,如长度尺寸链的极值计算。但装配误差往往同时包含平动与转动,标量无法表达耦合。

矢量表示把误差写成三维向量,如位置误差 δ = [δx, δy, δz]ᵀ,方向信息被引入。然而刚体在空间中有 6 个自由度,纯矢量只能描述 3 个平动分量,转动误差仍需单独处理。

矩阵表示用齐次变换矩阵统一描述位姿,是装配误差建模的主流工具。一个 4×4 齐次变换矩阵 T 同时包含旋转矩阵 R(3×3)与平移向量 p(3×1):

T = [ R  p ]
    [ 0  1 ]

当误差为微小量时,旋转矩阵可线性化,从而引出小位移旋量理论。

二、小位移旋量理论

小位移旋量把刚体的微小运动表示为 6 维向量:

ξ = [ ω ; v ] = [ωx, ωy, ωz, vx, vy, vz]ᵀ

其中 ω 为微小转动矢量,v 为参考点处的微小线位移。其物理意义是:刚体绕某轴转动 ω,同时沿该轴平移 v。对于小角度(通常小于 1° 或 0.017 rad),旋转矩阵可近似为:

R ≈ I + [ω]×

这里 [ω]× 是反对称矩阵。该线性化使误差可以叠加,是误差传递分析能够使用线性代数的关键前提。例如某轴系装配中,轴承内圈相对外圈产生 ω = [0, 0, 0.0005]ᵀ rad、v = [2, 0, 0]ᵀ μm,则旋量直接给出该结合面的完整误差状态。

三、齐次变换矩阵与误差旋量的映射

设理想位姿为 T₀,实际位姿为 T = T₀ · ΔT,其中 ΔT 为误差变换。在小位移假设下:

ΔT ≈ I + [ξ]×

即误差旋量 ξ 与误差变换矩阵一一对应。当误差从局部坐标系传递到全局坐标系时,需用伴随变换(Adjoint)进行坐标变换:

ξ_global = Ad_T · ξ_local

这一步是装配误差传递的核心:局部结合面的微小偏差,经位姿变换后成为末端执行器或测量点的误差贡献。

四、误差旋量到雅可比矩阵的映射

雅可比矩阵 J 描述关节/结合面变量到末端位姿的微分关系。若各结合面误差旋量为 ξᵢ,则末端总误差为:

ξ_end = Σ Jᵢ · ξᵢ

其中 Jᵢ 是第 i 个误差源到末端的传递雅可比。它把“局部误差”映射为“全局误差”,并直接给出敏感度:Jᵢ 中某列数值越大,说明该方向误差对末端影响越显著。这正是后续敏感度分析与精度分配的计算基础。

装配误差数学描述与传递链路 标量表示 δ = 8 μm 仅大小,无方向 矢量表示 δ = [δx, δy, δz]ᵀ 含方向,缺转动 齐次变换矩阵 T = [R p; 0 1] 完整位姿描述 小位移旋量 ξ = [ω; v] 6维,线性可叠加 小角度线性化:R ≈ I + [ω]× ,ΔT ≈ I + [ξ]× 适用条件:ω < 0.017 rad(约1°) 局部误差旋量 ξ_local 伴随变换 ξ_global = Ad_T · ξ_local 雅可比映射 ξ_end = Σ Jᵢ · ξᵢ 末端总误差与敏感度 Jᵢ 列值越大 → 该误差源影响越显著 适用边界 小变形、小角度、线性叠加

五、建模假设与适用边界

上述模型建立在若干假设之上,超出边界会导致显著误差:

  • 小位移假设:转动角需小于约 1°。若结合面出现大角度倾斜(如螺栓未拧紧导致的翘曲),线性化失效,必须使用完整非线性变换。
  • 刚体假设:忽略零件自身弹性变形。对于薄壁件、悬臂结构,受力后的变形可能大于装配几何误差,需引入柔度矩阵。
  • 误差独立假设:各误差源互不相关。若同一工序同时产生位置与角度偏差且存在耦合,需用协方差矩阵描述相关性。
  • 静态假设:不考虑热变形、振动与磨损。精密机构在温度梯度下,热误差可能主导总误差,需单独建模。
工程判断:当末端误差中转动分量引起的位移(ω × L,L 为力臂)超过平动分量时,必须保留旋量中的转动项;反之可简化为纯平动模型,降低计算量。

六、建模流程与示例参数

以某三轴精密装配平台为例,说明建模步骤:

  1. 定义全局坐标系与各结合面局部坐标系,明确正方向约定(右手定则,Z 轴向上)。
  2. 测量或给定各结合面误差旋量。示例假设:结合面 1 的 ω₁ = [0, 0, 0.0003]ᵀ rad,v₁ = [3, 1, 0]ᵀ μm;结合面 2 的 ω₂ = [0.0002, 0, 0]ᵀ rad,v₂ = [0, 2, 1]ᵀ μm。
  3. 计算各结合面到末端的伴随变换矩阵 Ad_T,将局部旋量转换到全局坐标系。
  4. 构造雅可比矩阵 J,计算末端总误差旋量 ξ_end = J₁ξ₁ + J₂ξ₂。
  5. 检查 ξ_end 中转动分量与平动分量的量级,判断是否需要保留完整旋量模型。

示例计算结果(假设值,仅用于说明方法):若 J₁ 对 Z 向转动的放大系数为 1.0,J₂ 对 X 向转动的放大系数为 0.8,则末端 Z 向转角约为 0.0003 + 0.8×0.0002 = 0.00046 rad,对应 100 mm 力臂处产生约 46 μm 的位移。这一数值已远超单结合面的 3 μm 平动误差,说明转动误差在长力臂下不可忽略。

常见误区:直接把各结合面的标量误差相加作为总误差。这种做法忽略了方向与力臂放大效应,在精密装配中可能低估总误差数倍。正确做法是至少使用旋量模型,并通过雅可比矩阵完成传递。

综上,标量适合快速估算,矢量适合方向明确的平动问题,而小位移旋量加齐次变换矩阵与雅可比映射,构成了精密装配误差建模的完整数学框架。理解其假设与边界,才能在实际工程中判断模型是否可信、何时需要升级为非线性或柔性模型。