天球参考系转换、勒让德函数计算是精密引潮力计算的两个核心工作.为了确保计算结果的准确性,采用两种彼此相互独立的算法公式进行检核计算.对于天球参考系转换,分别采用基于春分点的岁差章动转换和基于CIO的无旋转原点转换两种方法,转换参数分别采用IERS 2003、2010规范推荐的4组参数模型;对于勒让德函数计算,分别采用勒让德函数的递推算法、完全规格化缔合勒让德函数的递推算法两种数学模型;计算结果验证了引潮力计算过程与数据成果的准确性.以德国BFO测站为例,计算得到1962年1月1日至2015年8月1日54年内,时间间隔为1小时的精密引潮力时间序列,统计结果表明引潮力的量值在149459.841×10~(~(-11))ms~(-2)(10~(~(-11))ms~(-2)=1 nGal)以内.同时对精密引潮力计算中的各类影响因素进行计算与分析,计算结果表明:岁差章动模型更新、以及两种天球参考系转换方法之间的差异对引潮力的影响可完全忽略不计;地球扁率、极移、参考框架偏差、历表更新对引潮力影响的量值分别在1.891×10~(~(-11))ms~(-2)、0.586×10~(~(-11))ms~(-2)、0.032×10~(-11)ms~(-2)、0.012×10~(-11)ms~(-2)以内,这些因素在高精度引潮力计算中不能忽略;而地球时与世界时转换、径向法向转换对引潮力的影响分别达到916×10~(-11)ms~(-2)、476×10~(-11)ms~(-2),忽略这两个因素将会导致计算结果错误.
The celestial reference system transformation and the Legendre functions calculation are the key to the precise Tidal Generating Force (TGF) calculation. In order to ensure the accuracy of the TGF, two independent algorithms are adopted to check the calculation results. The equinox-based precession- nutation transformation method and the CIO-based non-rotating origin transformation method are used respectively in the celestial reference system transformation. Their transformation parameters come from the four models recommended by IERS Conventions 2003 and 2010. As to Legendre functions calculation, the Legendre function recursive algorithm and the fully normalized associated Legendre function recursive algorithm are employed to verify the accuracy of TGF calculation process and results. Taken BFO station in Germany as an example, the precise TGF time series data with 1 hour interval are calculated form January 1, 1962 to August 1, 2015. The results indicate that the TGF is a data within 149459. 841 × 10-11 ms-2 ( 10-11 ms-2 = 1 nGal). Meanwhile, the influence factors during the precise TGF calculation are analyzed. The results show that, the effects of the precession- nutation models updates and the difference between two celestial reference system transformation methods on the TGF calculation can be completely ignored. The effects of the earth flattening, polar motion, frame bias and ephemeris updates on the TGF calculation are within the value of 1. 891 ×10 -11 ms -2,0. 586×10 -11 ms -2, 0. 032 × 10 - 11 ms -2,0. 012 ×10 - 11 ms -2 respectively, which can not be ignored in the precise TGF calculation. The effects of the difference between Terrestrial Time (TT) and Universal Time ( UT1 ), the radial and normal direction transformation on the TGF calculation are 916 ×10-11 ms-2 and 476 × 10-11 ms-2 respectively. These two factors can not be neglected, which may lead to the inaccuracy of calculation results.