微椭弹性地球运动的一个边界条件的推导.pdfVIP

第42卷第1期 天 文 vol 学 报 42，No．1 2001年2月 ACTAA卵RoNoMIcAsINIcA Feb．。2001 微椭弹性地球运动的一个边界条件的推导* 黄乘利1，2廖新浩1朱耀仲2，3金文敬1 (1中国科学院上海天文台上海200030) (2中国科学院动力大地测量学开放研究实验室武汉430077) (3中国科学院测量与地球物理研究所武汉430077) 摘要 在研究地球章动或镧汐理论时，常常需要对均匀自转、徽椭、弹性、自引力的地球 的运动方程组积分，并通过选取一组恰当的边界条件来定解在一阶扁率近似下，先将椭球形 参考边界上一个有关形变的连续量转化到等效球面上，然后作广义面球谤函数展开进行标量 化，并分解为球形与环形部分，截断后可导出3个标量常微分形式的边界条件． 关羹词 章动。地球自转，边界条件 中圈分类号：P1833 l 引言 由于地球自转角速度的存在、以及地球形状非球对称而微椭，使得对地球自由振荡、 固体地球对日月行星作用的响应等研究变得更为繁杂．在研究地球章动或潮汐理论时，常 常需要对均匀自转、微椭、弹性、自引力的地球运动方程组积分，并通过选取一组恰当的边 界条件来定解． Dahlen【l o基于动量平衡、牛顿引力定律和应力应变关系给出了线性化的弹性地球的 动力学场方程，即Lagrange运动方程和Poisson方程．smith_“在此基础上，对运动方程作 surface Ha彻onics 广义面球谐函数展开(Generalized expansi。n，GssH)，利用本 spherical 方程的标量形式，然后对标量运动方程进行截断，并缩写成简洁的大矩阵形式，从而推导 出均匀自转、微椭、弹性、自引力的地球的运动方程组及其截断的标量形式．目前绝大多数 有关非剐体地球的章动和潮汐理论研究【3。o都是基于文[2]给出的运动方程组． 文[2]给出的弹性地球运动方程是一组常微分方程常微分方程要定解，即对常微分 运动方程组积分时，必须给定一些初值或边界条件．这些边界条件也需相应地展开并截 断．对受本征模(尤其是自由核章动模)频率强共振的章动项，合适边界条件的选择显得尤 为关键Ls’61．文[2]在一阶扁率近似下同时也给出了一组相应的边界连续条件，但由于其 中过程过于繁杂而没有给出具体的推导细节 04 2000—01一17收到原稿，200018收到修改蔫 10)资助 天 文 学 报 为了将这些运动方程组和边界条件推广到二阶扁率精度，以便计算并给出一个更高 精度的非刚体地球的章动模型，本文在一阶扁率近似下，从详细推导微椭弹性地球运动的 一个边界条件人手，旨在将其数理理论理解清楚． 2基本假设与坐标系 假设地球由固体内核、流体外核、弹性地幔所构成，该球微椭、刚性均匀地自转，处于 流体静平衡态(记为HsE)且具有各向同性的弹性连续关系．由于地球自引力的影响，相 对于位移而言并非小量，故在章动及固体潮汐这类全球动力学问题中必须顾及自引力的 影响 章动坐标系V。是空间惯性系与体固系之间的一个中介坐标系，其原点与地球质心 重合，其z轴在惯性空间中指向一个固定不变的方向，即平均的地球形状轴，x和Y轴以 常角速度n绕z轴旋转，即yE是一个旋转的坐标系．v。与体固系三轴之间的Euler角 很小，可视为一阶小量，因此在一个较短的时间内，v。适合于考察岁差、章动． 假设未受扰动时(或受一小扰动发生形变前)地球处于HSE态，此时的地球正好与 VE周联在一起；当有外力扰动时，地球内部体元的应力为初始应力与外部扰动所致的附 加应力之和，附加应力及应变皆为小量；设地球为理想弹性体，应力、应变之间可用广义胡 克定律来描述．可以证明【7