@inproceedings{oai:jaxa.repo.nii.ac.jp:00007265, author = {福島, 登志夫 and Fukushima, Toshio}, book = {アストロダイナミクスシンポジウム講演後刷り集, The 20th Workshop on JAXA Astrodynamics and Flight Mechanics}, month = {Mar}, note = {第20回アストロダイナミクスシンポジウム (2010年7月26-27日. 宇宙航空研究開発機構宇宙科学研究所), 相模原市, 神奈川県, The 20th Workshop on JAXA Astrodynamics and Flight Mechanics 2010 (July 26-27, 2010. Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency (JAXA)(ISAS)), Sagamihara, Kanagawa Japan, 土星周辺の探査機の軌道計算などでリングもしくはディスクの重力加速度場が必要となるが、球面調和関数による多重極展開は収束が悪いため定義に立ち戻って計算しなければならない。一様リングのポテンシャルは第一種完全楕円積分K(m)で表わされる(Kellogg 1929)が加速度場の具体的な表現は文献には見当たらないので導出するとK(m)と第二種完全楕円積分E(m)の線形結合で書けることがわかる。しかし、この表現は対称軸近傍で激しい桁落ちに直面するため、補助的な完全楕円積分B(m)およびS(m)を導入する必要が生じる。同様のことは一様ディスクの加速度場にも言える。ポテンシャルの表現はLass & Blitzer(1983)に、また加速度場の表現はKrogh et al.(1982)に既に与えられているが、この加速度表現も桁落ちが激しく使い物にならない。これも上記の補助的積分を使えば精度良く計算できる。実際の軌道計算にはB(m)およびS(m)の数値計算が必要であるが、前者の高速計算法はFukushima(2011)に、また後者の高速計算法はFukushima(2010b)に詳しい。なおK(m)およびE(m)の高速計算法についてはFukushima(2009b)を参照されたい。, We describe a precise numerical method to evaluate the acceleration vector of an inverse-square force law caused by a uniform riong or disk. The method consists of two parts. One is rewriting of the analytic expressions of the vector to reduce the loss of information in its computation. The other is a procedure to compute a special linear combination of the complete elliptic integrals of the first and the second kind appeared in the rewritten expressions. These details are already published in Fukushima (2010b). As for the fast computation procedures of the complete elliptic integrals and their auxiliaries are found in Fukushima (2009b, 2010b, 2011)., 資料番号: AA0065115018}, pages = {98--103}, publisher = {宇宙航空研究開発機構宇宙科学研究所 (JAXA)(ISAS), Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency (JAXA)(ISAS)}, title = {Precise Computation of Acceleration due to Uniform Ring or Disk}, volume = {20}, year = {2011} }