参考资料

地球的力学扁率可以通过弹性动力学理论来计算。根据地球作为一个弹性椭球体的假设，力学扁率与地球的平均密度（ρ）、引力加速度（g）、自转角速度（ω）、平均半径（R）、弹性模量（E）和泊松比（v）相关。具体计算公式为：

[ f = \frac{(1 - v)(R \omega^2)^2}{2gE} ]

其中：

• ( f ) 是地球的扁率；
• ( v ) 是泊松比；
• ( R ) 是地球的平均半径；
• ( \omega ) 是地球的自转角速度；
• ( g ) 是地球表面的重力加速度；
• ( E ) 是地球的弹性模量。

通过将地质历史时期的地球半径、质量和角速度变化值代入上述公式，可以计算出各地质历史时期地球扁率值，作为地球扁率变化的下限值。例如，根据新星云假说下地质历史时期半径、质量和角速度变化值代入公式，可以计算出各地质历史时期地球扁率值。

此外，地球扁率对天体引潮力的影响也需要考虑。在高精度天体引潮力计算中，需要考虑地球扁率的影响，通过DE421历表计算得到太阳系各个天体在GCRS中的坐标，并基于天球参考系经典的岁差章动转换方法、IERS 2010规范推荐的转换参数以及IERS发布的EOP 08 C04数据模型，综合考虑影响天体坐标转换的各种因素进行计算。