跳到主要内容

比力

解决运载体相对地球运动加速度计所敏感出的比力表达式

表达式
f=v˙+ωep×v+2ωie×vg\boldsymbol{f}=\dot{\boldsymbol{v}}+\omega_{e p} \times \boldsymbol{v}+2 \omega_{i e} \times \boldsymbol{v}-\boldsymbol{g}

加速度计力学模型:

image-20230612152409110

注意到敏感质量本身也会受到地球、月球、太阳等天体引力的作用,这个值会叠加到测量值上

x=mk(aG)\boldsymbol{x}=\frac{m}{k}(\boldsymbol{a}-\boldsymbol{G})

这个位移通过位移传感器就能转化为电信号,把加速度计的输入量(aG)(\boldsymbol{a}-\boldsymbol{G})称为“比力”,即:

f=aGf=a-G

其中引力加速度G包括了地球、月球、太阳等天体的影响

G=Ge+Gm+Gs+i=1n3Gi\boldsymbol{G}=\boldsymbol{G}_e+\boldsymbol{G}_m+\boldsymbol{G}_s+\sum_{i=1}^{n-3} \boldsymbol{G}_i
提示

比力代表了作用在单位质量上的弹簧力。因为比力的大小与弹簧变形量成正比,而加速度计输出电压的大小正是与弹簧变形量成正比,所以加速度计实际感测的量并非是运载体的加速度,而是比力。

另外,比力具有与加速度相同的量纲

在哥氏加速度一节中介绍了地球表面附近的运载体所受绝对加速度为

d2Rdt2i=d2R0dt2i+d2rdt2e+2ωie×drdte+ωie×(ωie×r)\left.\frac{d^2 \boldsymbol{R}}{d t^2}\right|_i=\left.\frac{d^2 \boldsymbol{R}_0}{d t^2}\right|_i+\left.\frac{d^2 \boldsymbol{r}}{d t^2}\right|_e+2 \omega_{i e} \times\left.\frac{d \boldsymbol{r}}{d t}\right|_e+\boldsymbol{\omega}_{i e} \times\left(\boldsymbol{\omega}_{i e} \times \boldsymbol{r}\right)

将这个绝对加速度带入比力的表达式f=aGf=a-G得:

f=d2R0dt2i+d2rdt2e+2ωie×drdte+ωie×(ωie×r)(Ge+Gm+Gs+i=1n3Gi)\boldsymbol{f}=\left.\frac{d^2 \boldsymbol{R}_0}{d t^2}\right|_i+\left.\frac{d^2 \boldsymbol{r}}{d t^2}\right|_e+2 \boldsymbol{\omega}_{i e} \times\left.\frac{d \boldsymbol{r}}{d t}\right|_e+\boldsymbol{\omega}_{i e} \times\left(\boldsymbol{\omega}_{i e} \times \boldsymbol{r}\right)-\left(\boldsymbol{G}_e+\boldsymbol{G}_m+\boldsymbol{G}_s+\sum_{i=1}^{n-3} \boldsymbol{G}_i\right)

下面对上式进行简化:

  • 一般来说,地球公转引起的向心加速度 d2R0dt2i\left.\frac{d^2 \boldsymbol{R}_0}{d t^2} \right|_i 与太阳引力加速度 Gs\boldsymbol{G}_s 的量值大致相等,可以约掉
  • 月球的引力加速度量值在106Ge10^{-6}G_e量级
  • 金星和木星的引力加速度量值都在108Ge10^{-8}G_e量级

最终加速度计测到的比力为:

f=d2rdt2e+2ωie×drdte+ωie×(ωie×r)Ge\boldsymbol{f}=\left.\frac{d^2 \boldsymbol{r}}{d t^2}\right|_e+2 \omega_{i e} \times\left.\frac{d \boldsymbol{r}}{d t}\right|_e+\omega_{i e} \times\left(\omega_{i e} \times \boldsymbol{r}\right)-\boldsymbol{G}_e

假设安装加速度计测量坐标系pp 系,它相对地球坐标系的转动角速度为 ωep\omega_{e p} ,则有

dvdte=dvdtp+ωep×v\left.\frac{d \boldsymbol{v}}{d t}\right|_e=\left.\frac{d \boldsymbol{v}}{d t}\right|_p+\boldsymbol{\omega}_{e p} \times \boldsymbol{v}

由此可以得到运载体相对地球运动时加速度计所敏感的比力表达式:

f=dvdte+2ωie×vgf=dvdtp+ωep×v+2ωie×vgf=v˙+ωep×v+2ωie×vg\begin{gathered} \boldsymbol{f}=\left.\frac{d \boldsymbol{v}}{d t}\right|_e+2 \boldsymbol{\omega}_{i e} \times \boldsymbol{v}-\boldsymbol{g} \\ \boldsymbol{f}=\left.\frac{d \boldsymbol{v}}{d t}\right|_p+\boldsymbol{\omega}_{e p} \times \boldsymbol{v}+2 \boldsymbol{\omega}_{i e} \times \boldsymbol{v}-\boldsymbol{g} \\ \boldsymbol{f}=\dot{\boldsymbol{v}}+\boldsymbol{\omega}_{e p} \times \boldsymbol{v}+2 \boldsymbol{\omega}_{i e} \times \boldsymbol{v}-\boldsymbol{g} \end{gathered}
f=v˙+ωep×v+2ωie×vg\boldsymbol{f}=\dot{\boldsymbol{v}}+\omega_{e p} \times \boldsymbol{v}+2 \omega_{i e} \times \boldsymbol{v}-\boldsymbol{g}

式中各物理量的含义:

  • v˙\dot{\boldsymbol{v}}:运载体相对于地球的速度在测量坐标系中的变化率
  • ωep×v\boldsymbol{\omega}_{e p} \times \boldsymbol{v}:测量坐标系相对地球转动引起的向心加速度
  • 2ωie×v2 \omega_{i e} \times \boldsymbol{v}:运载体相对地球速度与地球自转加速度相互影响所形成的哥氏加速度
  • g\boldsymbol{g}地球重力加速度

aB=ωep×v+2ωie×vg\boldsymbol{a}_B=\boldsymbol{\omega}_{e p} \times \boldsymbol{v}+2 \boldsymbol{\omega}_{i e} \times \boldsymbol{v}-\boldsymbol{g}

则有

faB=v˙\boldsymbol{f}-\boldsymbol{a}_B=\dot{\boldsymbol{v}}

导航计算中需要的是运载体相对地球的加速度 v˙\dot{v} ,因此必须从加速度计测得的比力中补偿有害加速度 aBa_B 的影响,才能得到运载体相对地球 的加速度,经过数学运算进而获得运载体相对地球的速度、位置等导航参数。

注意

此时可以知道,加速计测量出的加速度还有“向心加速度”和“哥氏加速度”组成的有害加速度的影响

根据安装加速度计的测量坐标系的不同,可以得到运载体相对地球运动时加速度计所敏感的不同比力表达式,通常称为比力方程

由于比力方程表明了加速度计所敏感的比力与运载体相对地球的加速度之间的关系,所以它是惯性系统的一个基本方程。不论惯性系统的具体方案和结构如何,该方程都是适用的。