4. 加速度计设计与误差建模
加速度计是惯性导航系统中最基础的敏感元件之一,它测量的是载体相对惯性空间的比力,而非直接测量运动加速度。理解这一点是设计加速度计与建立误差模型的前提:在静止状态下,加速度计输出的是当地重力加速度的反作用比力;在自由落体状态下,输出为零。惯导解算中,比力经姿态矩阵投影、扣除重力后积分,才能得到速度与位置。因此,加速度计的标度因数、零偏、安装误差和非线性误差会直接进入导航结果,且随时间累积。
典型加速度计的工作原理
工程上常见的加速度计包括石英挠性摆式、MEMS电容式和振梁式等。以应用最广的摆式加速度计为例,其核心是一个由挠性支承悬挂的检测质量(摆组件)。当载体存在沿输入轴的比力时,检测质量因惯性产生相对位移,位移被电容或光电传感器转换为电信号,经伺服放大器驱动力矩器产生平衡力矩,使质量回到零位。此时力矩器电流正比于输入比力,通过精密采样电阻转换为电压输出。这种“闭环力平衡”结构使标度因数稳定、线性度好,是中等精度惯导的常用方案。
MEMS电容式加速度计则采用梳齿结构,检测质量在比力作用下改变差分电容,通过电容-电压转换和闭环或开环解调得到输出。其体积小、成本低,但零偏稳定性和标度因数重复性通常弱于石英器件。振梁式加速度计利用谐振梁频率随轴向力变化的特性,输出为频率信号,抗干扰能力强,常用于高精度场景。
标度因数与零偏
标度因数 K₁ 定义为输出变化量与输入比力变化量之比,单位通常为 V/(m/s²) 或 mA/g。理想情况下 K₁ 为常数,但实际中它随温度、时间和输入量程变化。标度因数误差通常用相对误差 δK/K 表示,例如 100 ppm 表示在 1 g 输入下输出偏差为满量程的万分之一。零偏 K₀ 是输入为零时的输出,单位与输出相同,折算到输入轴后常用 g 或 m/s² 表示。零偏是惯导位置误差的主要来源之一:1 mg 的零偏在 100 s 内可造成约 0.5 m 的位置误差,在 1000 s 内可达数十米。
标度因数与零偏均需通过标定确定。以六位置标定为例,将加速度计分别置于 ±X、±Y、±Z 六个方向,每个位置采集静态输出。设第 i 个位置输入比力为 f_i,输出为 U_i,则可通过最小二乘拟合 U_i = K₁·f_i + K₀,得到 K₁ 和 K₀ 的估计值。若需同时辨识安装误差和非线性项,则需增加位置数量并采用多元回归。
安装误差与非线性误差
安装误差描述加速度计敏感轴与载体坐标系之间的失准角。对于三轴加速度计组合,理想情况下三个敏感轴正交且与载体轴重合,但实际存在小角度偏差。设加速度计 j 的敏感轴在载体坐标系中的方向余弦为 (cos α_j, cos β_j, cos γ_j),则输出可写为 U_j = K₁_j · (f_x cos α_j + f_y cos β_j + f_z cos γ_j) + K₀_j。当失准角较小时,cos α ≈ 1,sin α ≈ α,可线性化为 U_j = K₁_j · (f_j + Σ_{i≠j} ε_ji · f_i) + K₀_j,其中 ε_ji 为安装误差角。安装误差会导致轴间耦合,例如沿 X 轴的加速度在 Y 轴输出中产生虚假信号。
非线性误差表现为输出与输入不成正比,通常用二次项 K₂ 和三次项 K₃ 描述:U = K₁·f + K₀ + K₂·f² + K₃·f³。在 ±1 g 量程内,K₂ 和 K₃ 通常很小,但在大过载或高精度应用中不可忽略。非线性误差的标定需要多组不同幅值的输入,通过多项式拟合分离各阶系数。
误差建模与补偿方法
完整的加速度计误差模型可写为矩阵形式:U = K · f + B + N(f),其中 U 为三轴输出向量,K 为 3×3 标度因数与安装误差矩阵,f 为三轴比力向量,B 为零偏向量,N(f) 为非线性项向量。补偿过程即求逆:f̂ = K⁻¹ · (U − B − N(f))。由于 N(f) 与 f 相关,实际中常采用迭代或直接忽略高阶项。
标定流程通常包括:温度循环标定、多位置静态标定和动态标定。温度标定建立 K₁、K₀ 随温度变化的模型,常用多项式拟合;多位置静态标定确定 K 矩阵和 B 向量;动态标定则利用离心机或振动台施加已知正弦或阶跃输入,检验幅频特性和非线性。补偿后的残差应满足指标要求,例如零偏残差小于 50 μg,标度因数残差小于 100 ppm。
示例:某石英挠性加速度计标定参数
以下为一组教学示例参数(非实测数据,仅用于说明标定结果的形式):
| 参数 | 符号 | 标定值 | 单位 |
|---|---|---|---|
| 标度因数 | K₁ | 1.250 | V/g |
| 零偏 | K₀ | 0.0032 | V |
| 二次非线性 | K₂ | 0.00015 | V/g² |
| 三次非线性 | K₃ | −0.00002 | V/g³ |
| 安装误差角 | ε_xy | 0.0008 | rad |
| 零偏稳定性 | — | 15 | μg |
若输入比力为 0.5 g,则理论输出为 U = 1.250×0.5 + 0.0032 + 0.00015×0.25 − 0.00002×0.125 ≈ 0.6283 V。补偿时先减去零偏,再除以标度因数,并迭代扣除非线性项,即可得到比力估计值。实际工程中,这些系数会随温度变化,因此需在温度范围内建立系数-温度查找表或拟合多项式,在导航解算中实时补偿。
加速度计误差建模的核心在于:明确每一项误差的物理来源,通过标定获得系数,在解算中实时补偿,并用残差检验验证补偿效果。只有将标度因数、零偏、安装误差和非线性误差统一纳入模型,才能保证惯导系统在长时间工作下的精度和可靠性。