著者
横山 隆明 ヨコヤマ タカアキ Takaaki YOKOYAMA
巻号頁・発行日
(Released:2010-02-22)

本研究では、レゴリスといわれる細かい砂に覆われている月や惑星表面への着陸時の<br />着陸衝撃力および沈下量を事前に算定する方法に関する研究を行った。<br /> レゴリスのような細かい砂で覆われた地盤から受ける反力や沈下量を求めようとする<br />時、砂のような粉粒体は変形時に連続体的挙動から非連続体的挙動を示すため、解析的<br />に解くのは難しい。そのため地盤と機械との相互作用について研究しているテラメカニ<br />ックス分野では、解析対象を模擬した実験を行い、得られた実験結果から地盤をモデル<br />化し半経験式を導き、地盤の変形および抵抗を求める方法が良く用いられる。<br /> まず既存の研究の調査を行った結果、テラメカニックス分野で、砂地盤への着陸現象<br />を扱っている研究は見当たらなかったが、着陸と似たような現象を取り扱っている研究<br />として、重錘落下による地盤の締め固め効果に関する研究があった。これらの研究では、<br />衝撃加速度および沈下量は落下速度に比例すること、地盤密度が密な場合は衝撃加速度<br />波形が1つのピークを持つこと、緩い場合は2つのピークを持つこと、質量が大きい場<br />合は衝撃加速度ピーク値が小さくなることなどが実験から明らかにされていることがわ<br />かった。<br /> そこで本研究においても、まずはテラメカニックス分野で行われている方法と同じく、<br />実験から着陸衝撃力および沈下量を求める事を試みた。実験の対象は月面への着陸とし、<br />月面の高真空環境を模擬するため実験は真空槽の中で行い、月面の1/6G環境を模擬する<br />ため1/6G相似則に基づき1/6スケールの試験体を使用した。試験体は平成14年度に提<br />案されたSELENE-Bの着陸パットの諸元に基づき作成した。実験結果から、テラメカニッ<br />クス分野の重錘落下実験結果と同じく、衝撃加速度および沈下量は落下速度に比例する<br />ことがわかったが、月面環境を模擬するための真空引きの際に制御できない地盤密度の<br />変化が発生し、落下質量の影響が明確に現れない、既存の実験では密度が低い場合に現<br />れる2つのピークを持つ衝撃加速度形状が、密な地盤に対する実験でも現れるなどの真<br />空を模擬した実験特有の難点があることがわかった。<br /> そのため、アポロ計画で使われた地盤モデル中のパラメータを本研究で行ったフット<br />パット試験体の直径に合わせ変更して数値解析をおこない衝撃加速度および沈下量の予<br />測を行った。この結果、衝撃加速度および沈下量は速度に比例し、また、実験は表せな<br />かった落下質量の影響も表すことができることがわかった。しかしアポロ計画で使われ<br />たモデルでは、実験結果に見られた2つのピークを持つ衝撃加速度を再現できない、沈<br />下量がピーク値を示した後急激に減少してしまうなど、実験値と良く合わない点もあり<br />注意が必要であることが明らかになった。その結果、アポロの地盤モデルの適用性とし<br />ては、衝撃加速度、沈下量ともにピーク値までの予測であれば充分に使用可能であるこ<br />とがわかった。<br /> 次に、圧縮性流体の解析法として開発された連続体から非連続体まで解析できるSPH<br />法を砂地盤への衝突問題へ適用することを試みた。SPH法は圧縮性流体の解析法として<br />銀河系の形成など天体物理の問題解決のため開発された方法で、その後、非圧縮性流体<br />解析、構造解析など適用範囲が広がっているが、砂地盤への衝突題にSPH法を適用した<br />例は未だなく本研究が世界初である。本研究では地盤をモール・クーロンの破壊基準を<br />降伏条件として持つビンガム流体と仮定し解析を行った。その結果、実験と同様に衝撃<br />加速度および沈下量は衝突速度に比例し、実験では明らかにならなかった衝突質量の影<br />響、アポロモデルでは実験値と大きく異なったピーク後の沈下量の推移、およびアポロ<br />地盤モデルで表せなかった2つのピークをもつ衝撃加速度形状を表せることが明らかに<br />なった。そしてSPH法による解析結果を分析することで第2番目のピークの発生原因も<br />特定された。すなわち、実験の制約により制限せざるを得なかった地盤厚さをSPH法に<br />より任意に変化させ解析を行った解析結果、衝突中の地盤の密度変化の状態のモニター<br />結果から、地盤底面からの反射波および地盤密度の増加が原因であることが確かめられ<br />た。<br /> 以上の検討からSPH法は実験を充分に補える解析法であることが確かめられ、本研究<br />成果は月面上のみならず他の惑星上への着陸にも適用可能であることが示された。<br />最後に、これまでに得られた知見を用い、SELEM-B月探査機の月面着陸時の衝撃加速<br />度および沈下量の予測をSPH法を用いて行った。垂直方向のみ3mの衝突速度がある場合<br />を計算した結果、最大衝撃加速度は約21G程度(4脚接地の場合)、最大沈下量は約7.4<br />cm程度(1脚接地の場合)であることがわかった。<br /> 本研究では、月および惑星上への着陸時の衝撃加速度の予測法として、月面を対象と<br />した実験、地盤モデルによる数値計算、SPH法による解析の3つの方法を試みた。その<br />結果、月面環境を模擬しなければならない実験には多くの制限があり、アポロの際用い<br />られた地盤モデルによる数値計算ではピーク値までしか信頼性がないことがわかった。<br />反面、SPH法での解析は、実験で見られた特徴を良く表しているため、実験を補える手<br />段として信頼性が高いことが示され、月惑星上砂地盤への着陸衝撃力算定方法として信<br />頼性が高い方法であることがわかった。<br /> SPH法は現在様々な現象に対して適用範囲が広がりつつある方法であり、その特徴と<br />して連続体的挙動から非連続体的挙動まで表すことができるため、例えば着陸機が月お<br />よび惑星上の砂地盤に対して斜めに衝突する際の機体の安定性から砂の飛散範囲までの<br />解析などにも幅広く適用可能である。また、例えば今後活発になると思われる月の利用<br />に係わる問題に関しても、月の砂の掘削、運搬方法など月の砂に係わる問題は数多い。<br />SPH法を用いた砂の解析は、今回適用性が示された衝突問題以外にも幅広い適用性を持<br />つと考えられ、今後の発展が期待される手法と言える。<br />