Artemis計画はもうすぐ恒久的な月面有人拠点を建設するフェーズに入ります。ここで素朴な疑問が立ち上がります。月面に「建物」を作るとして、コンクリートや鉄骨を地球から運ぶのでしょうか。答えはノーです。月低軌道（LLO）まで1kgの物資を運ぶコストは現在およそ百万円規模、月面着陸まで含めれば数百万円になります。100t規模のシェルター1棟を地球から打ち上げると、それだけで数千億円の打ち上げ費が消えます。
代わりに、月面の足元に無尽蔵に転がっている レゴリス（regolith） ——表土として堆積した粒子状の岩石——を素材にして、現地で構造物を「印刷」する。これが ISRU（In-Situ Resource Utilization、現地資源利用） の発想です。レゴリスを太陽光やマイクロ波で焼結し、ロボットアームのノズルから層を重ねていけば、ドーム、道路、着陸パッド、放射線遮蔽壁が自前で作れます。地球から持っていくのは小さな印刷ヘッドとエネルギー源だけ。質量比で1000倍以上の「掛け算」が効きます。
ISRUの応用先は月面基地だけではありません。火星では大気から酸素と推進剤を作る計画（MOXIE実験）が既に走っています。小惑星では水資源を採取してデポを作る構想があります。共通する哲学は「持っていける質量を最小化し、現地の質量を最大限活かす」。これは宇宙開発のコスト構造を根本から書き換えるパラダイムです。
本記事の内容
- ISRUの定量的な意義 — なぜ「持参」より「現地調達」が桁違いに安いのか
- 月面レゴリスの組成・粒度・熱物性と、焼結に都合の良い理由
- 太陽光・レーザー・マイクロ波・電子ビーム焼結のエネルギー収支比較
- 3Dプリンティング手法（Contour Crafting / Powder Bed Fusion / Stereolithography）の月面適用
- 焼結物の機械的性質と多孔率モデル、熱伝導の有効媒質理論
- 月面基地・放射線遮蔽壁・着陸パッドの構造設計と簡易FEM
- Pythonによる熱伝導・焼結プロファイル・構造解析・太陽光集光のエネルギー収支シミュレーション
前提知識
この記事を読む前に、以下の記事を読んでおくと理解が深まります。
ISRUの質量経済学 — なぜ現地で作るのか
持参 vs 現地調達の桁違いな差
宇宙ミッションのコストは、究極的には「打ち上げ質量」に支配されます。SLS（Space Launch System）の月面投入コストはおおむね $\$1{,}000{,}000$/kg のオーダーです。Falcon Heavyや将来のStarshipで $\$50{,}000$/kg まで下がったとしても、100tの構造物を運ぶには5,000億円かかります。一方、月面に既にある質量（レゴリス）を使えば、運搬コストはほぼゼロです。
この差を 質量倍率（mass leverage ratio） で定式化できます。地球から運んだ機材の質量 $m_e$ kgで、月面の質量 $m_l$ kg分の構造物が作れるとき、
$$ R = \frac{m_l}{m_e} $$
これが大きいほどISRUの威力が大きい。3Dプリンタが10tで、印刷量が10,000tなら $R = 1000$。打ち上げコストが実質1/1000に圧縮されたことになります。NASAのEvolvable Mars Campaignの試算では、月面・火星拠点の長期運用では $R \geq 100$ がコスト的にブレークイーブンの目安とされています。
「ペイバック質量」の概念
ISRU装置自体も持ち込まなければなりません。装置質量 $m_{\mathrm{ISRU}}$、装置寿命中に生産する有用物質質量 $m_{\mathrm{prod}}$ とすると、ペイバックは
$$ m_{\mathrm{prod}} \geq m_{\mathrm{ISRU}} $$
で達成されます。3Dプリンタ重量1tで100tの構造物を出すなら100倍のペイバック。これは推進剤製造（水を電気分解して液酸・液水素にする）でも同じ思考で、装置質量と総生産量のトレードオフを最初に評価します。
ここまでで「現地で作るとなぜ得か」が掴めました。次の問いは「月面で本当に作れるのか」です。素材としてのレゴリスを物理化学的に見ていきます。
月面レゴリスの物理化学
組成と粒度分布
アポロ計画とLuna計画で持ち帰られたサンプル解析から、レゴリスは次のような酸化物の集合体だと分かっています（重量%、海領域の代表値）:
| 成分 | 重量% | 役割 |
|---|---|---|
| SiO₂（シリカ） | 40〜45 | ガラス・シリケート骨格 |
| Al₂O₃（アルミナ） | 12〜18 | 高融点骨格 |
| FeO（酸化鉄） | 12〜18 | マイクロ波吸収 |
| CaO | 10〜12 | フラックス（融点低下） |
| MgO | 8〜10 | 構造強化 |
| TiO₂ | 1〜10（海領域で高） | マイクロ波吸収・鉄抽出原料 |
ここで注目すべきは FeO と TiO₂ の存在 です。鉄とチタンの酸化物はマイクロ波（2.45 GHz帯）を効率よく吸収し、内部発熱で焼結温度（約 $1100$–$1300\ \mathrm{K}$）まで容易に達します。地球の砂（純シリカ寄り）にはこの性質がなく、これは月面の「ラッキーな組成」です。
粒度はおおむね $20$–$100\ \mu\mathrm{m}$ にピークを持ち、粒子表面はガラス質のアグルティネートで覆われています。微細粒（< $20\ \mu\mathrm{m}$）も多く、これは Powder Bed Fusion 系の3D印刷に都合の良い粒度です。
真空中での熱物性
月面は気圧 $10^{-12}$ Pa の超高真空、表面温度は赤道で昼 $400\ \mathrm{K}$／夜 $100\ \mathrm{K}$ の振幅です。真空中ではガス層を介した熱伝達がなく、レゴリスの熱伝導率は地球上の砂より大幅に低くなります。アポロ計測値はおおむね
$$ k_{\mathrm{regolith}} \approx 0.001\text{–}0.01\ \mathrm{W/m \cdot K} $$
このオーダーです。これは地球の砂（$0.3\ \mathrm{W/m \cdot K}$）の100分の1。熱がほとんど伝わらないため、太陽光やマイクロ波で局所加熱すると、加熱点のすぐ外側はほとんど熱を奪われません。これが「焼結の効率が真空でこそ高い」という直感に反する事実の原因です。
密度と多孔率
月面表層のバルク密度は深さによって変わり、表面 $0.1$ m で $\rho \approx 1500\ \mathrm{kg/m^3}$、深さ $1$ m で $1900\ \mathrm{kg/m^3}$ 程度です。粒子そのものの真密度 $\rho_s \approx 3100\ \mathrm{kg/m^3}$ と比べると多孔率は
$$ \phi = 1 – \frac{\rho}{\rho_s} \approx 0.4\text{–}0.5 $$
つまり体積の半分が真空の「すき間」。焼結で粒子同士を融着させると、この多孔率が下がり強度が上がる、という焼結プロセスの本質が見えてきます。
レゴリスの素性が分かったところで、これをどう焼結するかという技術論に入ります。エネルギー源と加熱方式の選択が、月面ISRUの実現性を決めます。
焼結プロセス — 4つのエネルギー源
焼結とは何か
焼結（sintering） は粉体を融点以下で加熱し、粒子界面の原子拡散によって粒子同士を融着させ、緻密化させるプロセスです。完全溶融よりエネルギーが少なく済み、結晶構造もコントロールしやすい。レゴリスの場合、融点は組成によりますが $1373$–$1573\ \mathrm{K}$ 程度、焼結温度はその $0.7$–$0.9$ 倍、すなわち $1000$–$1400\ \mathrm{K}$ の範囲です。
緻密化の駆動力は粒子表面エネルギー $\gamma_s$ の低減です。半径 $r$ の球状粒子が融着して半径 $R$ の塊になると、表面積比は $A_{\mathrm{after}}/A_{\mathrm{before}} = (r/R)^{1/3}$ で減ります。エネルギー利得は $\Delta E = \gamma_s \cdot \Delta A$。これが原子拡散を駆動します。
緻密化の速度は古典的にはCobleの理論で、密度 $\rho(t)$ の時間変化が
$$ \frac{d\rho}{dt} = K \frac{D(T)}{r^p}\,(\rho_s – \rho) $$
の形でモデル化できます。$D(T)$ は拡散係数、$r$ は粒径、$p$ は機構によって2〜4。$D(T) = D_0 \exp(-E_a/RT)$ なので、温度が上がると指数関数的に焼結が加速します。だからこそ「ピーク温度に到達して保持する」プロファイルが重要になります。
太陽光焼結
最もエネルギー的に有利なのが 太陽光集光焼結 です。月軌道の太陽定数は $S \approx 1361\ \mathrm{W/m^2}$、大気吸収がない分地表よりむしろ多く受光できます。$10\ \mathrm{m}^2$ のフレネルレンズで集光すれば $13.6\ \mathrm{kW}$、これを直径 $5\ \mathrm{mm}$ のスポットに絞ると、面積比 $5\times 10^5$ 倍。エネルギー密度は $\sim 7\times 10^6\ \mathrm{W/m^2}$ に達し、レゴリスを瞬時に焼結温度に到達させます。
利点は明白で、電源が要らない。欠点も明白で、夜は使えない、昼でも太陽高度が低い極域では集光効率が落ちる。コンセプト実証としてDLRが地上で太陽光シミュレータを用い、レゴリス相当物の焼結に成功しています。
マイクロ波焼結
レゴリス中のFeO・TiO₂が2.45 GHzマイクロ波を強く吸収する事実を利用します。マイクロ波の侵入深さ（penetration depth）は
$$ \delta = \frac{c}{2\pi f}\sqrt{\frac{2}{\varepsilon’\left(\sqrt{1+\tan^2\delta_\varepsilon}-1\right)}} $$
で、レゴリスでは $\delta \approx 0.1$–$0.3\ \mathrm{m}$ 程度になり、表面だけでなく体積加熱ができるのが特徴です。これは太陽光焼結が表面しか加熱できないのと対照的で、厚い壁を一気に焼結する用途に向きます。
NASAのLERCIPでサンプルレゴリスをマイクロ波加熱した実験では、4分で $1373\ \mathrm{K}$ 到達、密度 $2.4\ \mathrm{g/cm^3}$（理論密度の80%）の焼結体が得られています。
レーザー焼結（SLS/SLM）
地球の金属3Dプリンティングで主流の手法をそのまま月に持ってくる発想です。$\mathrm{CO_2}$ レーザー（波長 $10.6\ \mu\mathrm{m}$）やファイバーレーザー（$1.06\ \mu\mathrm{m}$）でレゴリス粉末層を選択的にスキャンし、層ごとに焼結します。精度が高く、複雑形状が作れるのが強みです。
弱点は電力効率。レーザーの電気→光変換効率は $30$–$50$%、光→熱変換でさらにロスがあり、総合効率は太陽光焼結に大きく劣ります。電源を太陽電池+蓄電池で賄う前提だと、必要な太陽電池面積が大きくなります。
電子ビーム焼結
電子線で粒子を加熱する方式。真空中で稼働するのが大気下より有利で、月面とは相性が良い。発生効率も高い（$\sim 70$%）。一方、電子線は表面しか加熱できず、また磁場による偏向制御が要ります。
ここまで4つの焼結プロセスを並べました。次の問いは「これらをどう積層して構造物にするか」、すなわち3Dプリンティング手法の選択です。
月面適用される3Dプリンティング手法
Contour Crafting（CC）
南カリフォルニア大のKhoshnevis教授が提唱した、コンクリート3Dプリンティングの代表手法を月面に適用したものです。ロボットアームの先端にノズルがあり、レゴリス+結合剤のペーストを連続押出ししながら壁を「型枠なしで」築き上げます。仕上げ用のトロウェル（こて）が同時に動き、層間を平滑化します。
月面でのレゴリス用CCでは、結合剤として（i）アンモニア水ベース、（ii）硫黄結合（月面玄武岩に硫黄を加えて加熱）、（iii）水なしでマイクロ波で焼結、の3案が研究されています。利点は壁を$1\ \mathrm{m}$/時のオーダーで高速に伸ばせること。NASA Project Olympusと米国ICONの月面構造物計画はこの方式が中核です。
Powder Bed Fusion（PBF）
工業用金属3Dプリンタと同じ手法。レゴリス粉末を薄く敷き、レーザーや電子ビームでスキャンして焼結、次の層を敷いて繰り返す。精度は $0.1\ \mathrm{mm}$ レベルで、機械部品クラスの精度が出ます。
欠点は遅さ。一層 $50\ \mu\mathrm{m}$ ずつ積層するため、 $1\ \mathrm{m}$ 高さの構造物には $20{,}000$ 層必要。建物には向かず、機械部品や精密ジグの製造が向く。ESAの3D printed lunar base BlueprintはPBFと外殻のCCを組み合わせるハイブリッド案です。
Stereolithography（SLA）と派生
紫外線硬化性樹脂を月面で使うのは現実的でないですが、レゴリス焼結の前段階として、レゴリスバインダー混合スラリーを光硬化させて中間構造を作る研究もあります。月面では太陽光UVが豊富なので、UV光源を別途持ち込まなくてよい利点があります。
各手法の比較
| 手法 | 速度 | 精度 | 強度 | 月面適合度 |
|---|---|---|---|---|
| Contour Crafting | 高 | 中 | 中 | ◎（建物） |
| Powder Bed Fusion | 低 | 高 | 高 | ◯（部品） |
| Stereolithography | 中 | 高 | 中 | △（樹脂要） |
| Direct Energy Deposition | 高 | 低 | 中 | ◯（修理） |
実際の月面建設では、外殻ドームをCC、内部の機器架台をPBF、修理保守をDEDで、と複数手法を併用するのが現実解とされています。
手法が選べたところで、出来上がる物がどんな機械的性質を持つのか、定量的に評価していきます。
焼結物の機械的性質
密度と強度の関係
焼結体の密度 $\rho$ と圧縮強度 $\sigma_c$ の関係は、多孔率 $\phi = 1 – \rho/\rho_s$ を使って経験的に
$$ \sigma_c = \sigma_{c,0}\,(1-\phi)^n $$
で表されます（Spriggsの式、$n \approx 2$–$4$）。$\sigma_{c,0}$ は完全緻密体の強度。レゴリス焼結体の場合、文献値は $\sigma_{c,0} \approx 200\ \mathrm{MPa}$、$n \approx 3$。$\phi = 0.2$（焼結良好）で $\sigma_c \approx 100\ \mathrm{MPa}$、$\phi = 0.5$（焼結不十分）で $\sigma_c \approx 25\ \mathrm{MPa}$。
参考までに、地球のコンクリートが $20$–$50\ \mathrm{MPa}$、岩石が $50$–$200\ \mathrm{MPa}$ 程度。焼結レゴリスはコンクリートと岩石の中間で、十分実用に耐えます。
引張強度と破壊靭性
セラミック系材料に共通する弱点として、引張に弱い。引張強度 $\sigma_t$ は圧縮強度の $1/10$ 程度。建築では「圧縮で受ける」設計、すなわちアーチ構造やドーム構造が必須です。これはローマ時代の石造建築と同じ思想で、月面ドームの形状設計に直結します。
破壊靭性 $K_{IC} \approx 0.5\text{–}1.5\ \mathrm{MPa\sqrt{m}}$。鉄鋼の $50\ \mathrm{MPa\sqrt{m}}$ に比べると2桁低く、衝撃には弱い。隕石衝突対策には別途多層構造（Whippleシールド類似）が要ります。
熱伝導率と熱応力
焼結レゴリスの熱伝導率は緻密化に応じて $0.5$–$2\ \mathrm{W/m \cdot K}$ 程度。地球の煉瓦と同水準です。低い熱伝導率は基地の断熱という観点では好都合ですが、温度サイクル（月の昼夜で $300\ \mathrm{K}$ 振幅）で熱応力が発生します。熱応力は
$$ \sigma_{\mathrm{th}} = E\,\alpha\,\Delta T $$
レゴリス焼結体の $E \approx 30\ \mathrm{GPa}$、$\alpha \approx 8 \times 10^{-6}\ \mathrm{K^{-1}}$、$\Delta T = 300\ \mathrm{K}$ で $\sigma_{\mathrm{th}} \approx 72\ \mathrm{MPa}$。これは引張強度を超える危険水準で、月の昼夜境界（terminator）に位置する構造物では深刻です。設計上は「日陰の永続クレーター」か「地下シェルター」で温度振幅を緩和するアプローチが取られます。
材料側の性質が見えたところで、次は実際に何を作るかという用途別の設計に進みます。
用途別の構造設計
月面基地（ESA Lunar Village）
ESAが提唱したコンセプトは、レゴリス3Dプリントで作る半球ドームです。
- 外殻: 厚さ $1$–$3$ m のレゴリス焼結シェル（放射線・隕石遮蔽兼用）
- 内殻: 加圧用の地球持参インフレータブル膜
- 床: レゴリス焼結床+断熱層
直径 $10\ \mathrm{m}$ のドームを厚さ $1\ \mathrm{m}$ の焼結レゴリスで作ると、必要レゴリス量は
$$ V = \frac{2\pi}{3}\left(R_o^3 – R_i^3\right) = \frac{2\pi}{3}\left(6^3 – 5^3\right) \approx 191\ \mathrm{m}^3 $$
これは質量 $\rho V \approx 2000 \times 191 \approx 380\ \mathrm{t}$ に相当。地球から運ぶなら2000億円超、現地調達なら掘削エネルギーだけで済みます。
放射線遮蔽壁
月面には大気がないため、宇宙線（銀河宇宙線GCRと太陽高エネルギー粒子SPE）が直接到達します。長期滞在者の被ばく線量を $50\ \mathrm{mSv/year}$ 以下に抑えるには、$1$–$2\ \mathrm{m}$ 厚のレゴリス遮蔽が必要とされています。線量低減率は厚さ $d$ に対して
$$ D(d) = D_0\,\exp\!\left(-\frac{d}{\lambda}\right) $$
レゴリスの減衰長 $\lambda \approx 0.7\ \mathrm{m}$（GCRに対して）。$d = 1\ \mathrm{m}$ で $D/D_0 = e^{-1.43} \approx 0.24$、約4分の1に。$d = 2\ \mathrm{m}$ で約 $1/16$ になります。これがドーム外殻厚さの設計根拠です。
着陸パッド
月面着陸機（HLS、Starship HLS、Blue Moonなど）が降下すると、噴射ガスがレゴリスを舞い上げます。これは隣接機器を傷つけ、視界を遮ります。アポロ12号着陸時、近隣のSurveyor 3に塗布されていた塗料が剥がれていた、という有名な観測例があります。
対策は焼結レゴリスで地表を硬化させること。厚さ $20\ \mathrm{cm}$、半径 $20\ \mathrm{m}$ のパッドを焼結すると、必要レゴリス量は $\pi r^2 h \approx 251\ \mathrm{m}^3$。太陽光集光焼結 $10\ \mathrm{kW}$ 規模で数百時間で完成する見積りです。
道路と未舗装トラフィックの抑制
基地周辺の移動経路を焼結すると、レゴリス粉塵の舞い上がりが大幅に減ります。粉塵は機械可動部やシール、太陽電池表面に付着して長期故障の原因になるため、道路化はメンテナンスコストを劇的に下げます。
用途設計が見えました。ここからは数値で検証していきます。Pythonで主要な現象を実装し、設計値が物理的に整合しているかを確かめます。
Pythonでの実装
太陽光集光焼結のエネルギー収支
まず、太陽光集光で1つのスポットをどれくらいの速さで焼結温度まで加熱できるか、エネルギー収支を解きます。スポットへの入射パワーから、放射と内部熱伝導による損失を引いた正味流入が加熱を引き起こす、というモデルです。
import numpy as np
import matplotlib.pyplot as plt
# 物理定数と月面パラメータ
S0 = 1361.0 # 太陽定数 [W/m^2]
sigma_SB = 5.67e-8 # ステファン・ボルツマン定数 [W/m^2/K^4]
T_amb = 250.0 # 周囲温度 [K]（昼側赤道の代表値）
# レゴリス物性
rho = 1900.0 # バルク密度 [kg/m^3]
c_p = 800.0 # 比熱 [J/kg/K]
k_th = 0.005 # 熱伝導率（粉末状態）[W/m/K]
eps_th = 0.95 # 放射率（焼結中の高温域）
# 集光器設計
A_collect = 10.0 # 集光面積 [m^2]
eta_opt = 0.85 # 光学効率
d_spot = 5e-3 # スポット直径 [m]
A_spot = np.pi * (d_spot/2)**2
depth = 5e-3 # 焼結体の厚さ [m]
# 集光パワー
P_in = S0 * A_collect * eta_opt
print(f"集光パワー: {P_in:.1f} W")
print(f"スポット面積: {A_spot*1e6:.2f} mm^2")
print(f"スポットでの放射照度: {P_in/A_spot:.2e} W/m^2")
# 加熱対象の質量
m = rho * A_spot * depth
print(f"対象質量: {m*1e6:.2f} mg")
# 時間積分（陽解法）でT(t)を求める
dt = 0.001
t_max = 30.0
nt = int(t_max/dt)
T = np.zeros(nt)
T[0] = T_amb
for i in range(1, nt):
# 入射 - 放射損失 - 伝導損失
Q_in = P_in
Q_rad = eps_th * sigma_SB * A_spot * (T[i-1]**4 - T_amb**4)
Q_cond = k_th * A_spot * (T[i-1] - T_amb) / depth # 簡易な定常近似
dTdt = (Q_in - Q_rad - Q_cond) / (m * c_p)
T[i] = T[i-1] + dTdt * dt
t = np.arange(nt) * dt
T_sinter = 1200.0 # 焼結温度の目安 [K]
i_reach = np.argmax(T > T_sinter) if np.any(T > T_sinter) else None
plt.figure(figsize=(9, 5))
plt.plot(t, T, 'b-', lw=2)
plt.axhline(T_sinter, color='r', ls='--', label=f'Sintering T = {T_sinter} K')
if i_reach is not None:
plt.axvline(t[i_reach], color='g', ls=':', label=f'Reach at t={t[i_reach]:.2f} s')
plt.xlabel('Time [s]')
plt.ylabel('Spot temperature [K]')
plt.title('Solar concentrator sintering: temperature rise')
plt.legend()
plt.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('solar_sintering.png', dpi=150, bbox_inches='tight')
plt.show()
このグラフから、3つのことが読み取れます。第一に、集光開始から $1$–$2$ 秒で温度が急上昇し、焼結温度 $1200\ \mathrm{K}$ に到達します。これは集光率の高さ（$10\ \mathrm{m}^2 \to 5\ \mathrm{mm}^2$ で $5 \times 10^5$ 倍）の威力です。第二に、温度はそのまま発散せず、ある値で漸近的に飽和します。これは温度上昇とともに放射損失 $\propto T^4$ が急増し、入熱と釣り合うためです。第三に、レゴリスの低熱伝導性 $k_{th} = 0.005\ \mathrm{W/m/K}$ のおかげで、伝導損失が支配的にならず、放射損失だけが律速となります。だからこそ真空中の月面で太陽光焼結が成立するわけです。
マイクロ波焼結の侵入深さ
次に、マイクロ波がレゴリス内部にどこまで届くかを、誘電損失角 $\tan\delta_\varepsilon$ から計算します。FeO・TiO₂の含有量が大きいほど吸収が大きく、侵入深さが浅くなる傾向を確認します。
import numpy as np
import matplotlib.pyplot as plt
# マイクロ波周波数
f_mw = 2.45e9 # [Hz]
c = 3e8 # [m/s]
omega = 2*np.pi*f_mw
# レゴリスの誘電特性（FeO含有率に対する経験式）
def eps_regolith(feo_pct):
"""FeO含有率(wt%)に対する複素誘電率の経験モデル
出典: Olhoeft 1989, Hoekzema 2008 などからの近似"""
eps_r = 3.0 + 0.05 * feo_pct # 実部
tan_d = 0.005 + 0.003 * feo_pct # 損失角
return eps_r, tan_d
# 侵入深さ
def penetration_depth(eps_r, tan_d):
"""半振幅の到達深さ [m]"""
return c / (2*np.pi*f_mw) * \
np.sqrt(2.0 / (eps_r * (np.sqrt(1+tan_d**2) - 1)))
# 海領域(mare)と高地(highland)の比較
feo_range = np.linspace(2, 18, 50)
delta = np.array([penetration_depth(*eps_regolith(x)) for x in feo_range])
plt.figure(figsize=(9, 5))
plt.semilogy(feo_range, delta, 'o-', color='darkorange')
plt.axhline(0.3, color='g', ls='--', alpha=0.5, label='Bulk heating regime')
plt.axhline(0.01, color='r', ls='--', alpha=0.5, label='Surface heating regime')
plt.xlabel('FeO content [wt%]')
plt.ylabel('Microwave penetration depth [m]')
plt.title(f'Microwave penetration vs FeO content (f={f_mw/1e9:.2f} GHz)')
plt.legend()
plt.grid(True, which='both', alpha=0.3)
plt.tight_layout()
plt.savefig('mw_penetration.png', dpi=150, bbox_inches='tight')
plt.show()
# 代表値の出力
for feo in [4, 10, 16]:
eps_r, td = eps_regolith(feo)
d = penetration_depth(eps_r, td)
print(f"FeO={feo}%: eps_r={eps_r:.2f}, tan_d={td:.4f}, depth={d*100:.1f} cm")
このグラフから2つの傾向が見えます。第一に、FeO含有量が増えると侵入深さは単調に減少します。これは損失角 $\tan\delta_\varepsilon$ が上昇し、マイクロ波が早く吸収されるためです。第二に、月の海領域（FeO 16% 程度）では侵入深さが $10\ \mathrm{cm}$ 前後、高地（FeO 4% 程度）では $50\ \mathrm{cm}$ 程度になります。建設に都合がいいのは中間値で、$0.3$–$0.5$ m の壁全体を均一に体積加熱できる「ちょうどよい」吸収特性を持つ地域を選ぶことになります。これは月面拠点立地選定の隠れた制約条件です。
Coble焼結による密度進化
次に、温度プロファイルが与えられたとき、密度がどう時間進化するかをCobleの式で解きます。温度を上げて保持→冷却、という典型的な焼結サイクルでどこまで緻密化が進むかを確認します。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
rho_s = 3100.0 # 真密度 [kg/m^3]
rho0 = 1500.0 # 初期バルク密度 [kg/m^3]
D0 = 1e-4 # 拡散係数の前因子 [m^2/s]
Ea = 250e3 # 活性化エネルギー [J/mol]
R_gas = 8.314 # ガス定数
r_grain = 50e-6 # 粒径 [m]
p_exp = 3 # Cobleの粒径指数
K_const = 1.0 # スケーリング定数
# 温度プロファイル(典型的な焼結サイクル: 上げ→保持→下げ)
def T_profile(t):
if t < 600:
return 300 + (1200 - 300) * (t/600) # 10分かけて昇温
elif t < 1800:
return 1300 # 20分保持
else:
return max(300, 1300 - (1300 - 300)*(t-1800)/1200) # 20分かけて降温
# 時間積分
dt = 0.5
t_max = 3000.0
nt = int(t_max/dt)
rho = np.zeros(nt)
T_hist = np.zeros(nt)
rho[0] = rho0
for i in range(1, nt):
T = T_profile(i*dt)
T_hist[i] = T
D_T = D0 * np.exp(-Ea/(R_gas*T))
drhodt = K_const * D_T / r_grain**p_exp * (rho_s - rho[i-1])
rho[i] = rho[i-1] + drhodt * dt
if rho[i] > rho_s:
rho[i] = rho_s
t = np.arange(nt) * dt / 60.0 # minutes
phi = 1 - rho/rho_s
fig, axes = plt.subplots(2, 1, figsize=(10, 7), sharex=True)
axes[0].plot(t, T_hist, 'r-', lw=2)
axes[0].set_ylabel('Temperature [K]')
axes[0].set_title('Sintering thermal profile and densification')
axes[0].grid(True, alpha=0.3)
axes[1].plot(t, rho, 'b-', lw=2, label='Density')
axes[1].axhline(rho_s, color='gray', ls=':', label=f'True density={rho_s} kg/m^3')
axes[1].set_xlabel('Time [min]')
axes[1].set_ylabel('Bulk density [kg/m^3]')
ax2 = axes[1].twinx()
ax2.plot(t, phi, 'g--', lw=1.5, label='Porosity')
ax2.set_ylabel('Porosity', color='g')
axes[1].legend(loc='lower right')
axes[1].grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('sintering_profile.png', dpi=150, bbox_inches='tight')
plt.show()
print(f"最終密度: {rho[-1]:.0f} kg/m^3")
print(f"最終多孔率: {phi[-1]:.3f}")
print(f"相対密度: {rho[-1]/rho_s*100:.1f} %")
上下2段のグラフから焼結プロセスの本質が見えます。上段の温度プロファイルが定常加熱→保持→冷却の台形を描く間、下段の密度はその「保持」部分で集中的に上昇しています。これは拡散係数 $D(T) = D_0\exp(-E_a/RT)$ が温度の指数関数になっているためで、$1300\ \mathrm{K}$ の高温域だけが焼結に実質寄与していることを意味します。昇温・降温期はほぼ寄与せず、緻密化は保持時間でほぼ決まる、という工学的に重要な事実が読み取れます。最終的に相対密度80%強（多孔率 $\phi \approx 0.18$）まで達しており、これは前節の強度式から $\sigma_c \approx 110\ \mathrm{MPa}$ に対応します。建材として十分です。
レゴリスの有効熱伝導率（多孔体モデル）
焼結が進むと多孔率が下がり、熱伝導率も上がります。これを混合則（Maxwell-Eucken）でモデル化し、設計用の有効物性を求めます。
import numpy as np
import matplotlib.pyplot as plt
# 構成相の熱伝導率
k_solid = 2.0 # 焼結ガラス・シリケート相 [W/m/K]
k_void = 0.0001 # 真空（実質的にほぼゼロ）[W/m/K]
# 各種混合則
def k_series(phi, ks, kv):
"""直列モデル（下限）"""
return 1.0 / ((1-phi)/ks + phi/kv) if phi < 1.0 else kv
def k_parallel(phi, ks, kv):
"""並列モデル（上限）"""
return (1-phi)*ks + phi*kv
def k_maxwell(phi, ks, kv):
"""Maxwell-Eucken（球状空隙の希釈分散）"""
return ks * (kv + 2*ks - 2*phi*(ks - kv)) / (kv + 2*ks + phi*(ks - kv))
def k_emt(phi, ks, kv):
"""Bruggeman有効媒質近似（陽解）"""
# 2*(1-phi)*(ks-keff)/(ks+2*keff) + 2*phi*(kv-keff)/(kv+2*keff) = 0
# を解く（二次方程式に帰着）
a = -2.0
b = (1-phi)*(2*ks - kv) + phi*(2*kv - ks)
c = ks * kv
disc = b**2 - 4*a*c
return (-b - np.sqrt(disc)) / (2*a)
phi_range = np.linspace(0.0, 0.7, 100)
k_s = np.array([k_series(p, k_solid, k_void) for p in phi_range])
k_p = np.array([k_parallel(p, k_solid, k_void) for p in phi_range])
k_m = np.array([k_maxwell(p, k_solid, k_void) for p in phi_range])
k_e = np.array([k_emt(p, k_solid, k_void) for p in phi_range])
plt.figure(figsize=(9, 5))
plt.plot(phi_range, k_p, label='Parallel (upper)')
plt.plot(phi_range, k_m, label='Maxwell-Eucken')
plt.plot(phi_range, k_e, label='Bruggeman EMT')
plt.plot(phi_range, k_s, label='Series (lower)')
plt.axvspan(0.15, 0.30, alpha=0.2, color='green', label='Well-sintered')
plt.axvspan(0.40, 0.55, alpha=0.2, color='red', label='Raw regolith')
plt.xlabel('Porosity phi')
plt.ylabel('Effective thermal conductivity [W/m/K]')
plt.title('Effective k of porous sintered regolith')
plt.legend()
plt.yscale('log')
plt.grid(True, which='both', alpha=0.3)
plt.tight_layout()
plt.savefig('effective_k.png', dpi=150, bbox_inches='tight')
plt.show()
# 焼結品質ごとの代表値
for phi, label in [(0.20, '良好焼結'), (0.45, '生レゴリス')]:
print(f"{label} phi={phi}: "
f"Maxwell={k_maxwell(phi, k_solid, k_void):.3f}, "
f"EMT={k_emt(phi, k_solid, k_void):.3f} W/m/K")
このグラフから3つの設計指針が得られます。まず、生レゴリス（多孔率 $0.4$–$0.55$）では有効熱伝導率は $0.5\ \mathrm{W/m \cdot K}$ 程度。焼結後（$0.15$–$0.30$）では $1.5\ \mathrm{W/m \cdot K}$ 程度まで上がります。次に、直列モデルと並列モデルの上下限を比較すると、真空の空隙は直列方向で熱を強く阻害するのが分かります。最後に、Maxwell-Eucken と Bruggeman の中間モデルがほぼ重なっており、設計目的にはどちらを使っても大差ないことが示唆されます。基地外殻の熱設計では多孔率を意図的にコントロールして「内部断熱・外殻放熱」のグラデーション構造を作る、というアイデアにつながります。
放射線遮蔽厚の設計
GCRの遮蔽効率を厚さの関数として計算し、安全な居住に必要な厚さを設計します。
import numpy as np
import matplotlib.pyplot as plt
# 月面の宇宙線フラックスから期待される被ばく
D0_GCR = 380.0 # 遮蔽なしのGCR年線量 [mSv/year]（太陽極小期）
D0_SPE_peak = 200.0 # 大型SPE 1回での被ばく [mSv]
lambda_GCR = 0.65 # GCRの減衰長 [m]（焼結レゴリス）
lambda_SPE = 0.25 # SPEの減衰長 [m]
D_limit = 50.0 # 年間被ばく目標 [mSv/year]
d = np.linspace(0, 3.0, 200)
D_GCR = D0_GCR * np.exp(-d/lambda_GCR)
D_SPE = D0_SPE_peak * np.exp(-d/lambda_SPE)
D_total = D_GCR + 0.5*D_SPE # 0.5は年内SPE回数の期待値補正
plt.figure(figsize=(9, 5))
plt.plot(d, D_GCR, label='GCR (annual)')
plt.plot(d, D_SPE, label='SPE peak event')
plt.plot(d, D_total, 'k-', lw=2, label='Total (effective annual)')
plt.axhline(D_limit, color='r', ls='--', label=f'Limit = {D_limit} mSv/yr')
# 必要厚さの数値解
idx = np.argmax(D_total <= D_limit)
d_required = d[idx]
plt.axvline(d_required, color='g', ls=':', label=f'd_req = {d_required:.2f} m')
plt.xlabel('Regolith shielding thickness d [m]')
plt.ylabel('Dose [mSv]')
plt.title('Radiation shielding by sintered regolith wall')
plt.legend()
plt.yscale('log')
plt.grid(True, which='both', alpha=0.3)
plt.tight_layout()
plt.savefig('radiation_shielding.png', dpi=150, bbox_inches='tight')
plt.show()
print(f"50 mSv/year以下にするための必要厚さ: {d_required:.2f} m")
print(f"これに必要なレゴリス量（直径10mドーム外殻）:"
f" {2*np.pi/3*((5+d_required)**3 - 5**3):.0f} m^3")
このグラフが示すのは、放射線遮蔽の設計ロジックそのものです。GCRは指数的減衰が比較的緩い（減衰長 $0.65\ \mathrm{m}$）のに対し、SPEは粒子エネルギーが低めで減衰長 $0.25\ \mathrm{m}$ と急峻に落ちます。両者を合わせた年間被ばくが目標 $50\ \mathrm{mSv}$ を切るには、厚さ約 $1.3$ m が必要、と読み取れます。前節の幾何学計算と合わせると、直径 $10\ \mathrm{m}$ ドームに必要なレゴリス量は約 $400\ \mathrm{m}^3$。これは数百tスケールで、3Dプリンタが数千時間連続稼働する見積りに対応します。月面1日24時間×太陽光集光中心で考えると、ドーム1棟あたり数地球月の建設期間が現実的な範囲となります。
円筒壁の簡易FEMによる構造解析
ドーム壁の応力分布を、軸対称1次元FEMで近似計算します。内圧（与圧）と自重による応力を見ます。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
R_in = 5.0 # 内半径 [m]
R_out = 6.3 # 外半径 [m]（厚さ1.3 m）
P_in = 1.0e5 # 内圧 [Pa]（1気圧）
g_moon = 1.62 # 月面重力 [m/s^2]
rho_sint = 2400.0 # 焼結体密度
E_sint = 30e9 # ヤング率 [Pa]
nu = 0.25 # ポアソン比
# 1次元軸対称FEM（厚肉円筒のラメ解+自重）
N = 50
r = np.linspace(R_in, R_out, N+1)
r_c = 0.5*(r[:-1] + r[1:]) # 要素中心
# ラメ解（厚肉円筒, 内圧のみ、外圧ゼロ）
A_lame = P_in * R_in**2 / (R_out**2 - R_in**2)
B_lame = P_in * R_in**2 * R_out**2 / (R_out**2 - R_in**2)
sigma_r = A_lame - B_lame/r_c**2
sigma_theta = A_lame + B_lame/r_c**2
# 自重応力(極座標近似: 球殻ではなく円筒として簡略化)
sigma_self = rho_sint * g_moon * (R_out - r_c)
# 主応力
sigma_max = sigma_theta + sigma_self
sigma_min = sigma_r
plt.figure(figsize=(10, 6))
plt.plot(r_c, sigma_r/1e6, 'b-', label='Radial stress')
plt.plot(r_c, sigma_theta/1e6, 'r-', label='Hoop stress')
plt.plot(r_c, sigma_self/1e6, 'g--', label='Self-weight contribution')
plt.plot(r_c, sigma_max/1e6, 'k-', lw=2, label='Max combined (tensile)')
plt.axhline(y=10.0, color='purple', ls=':', label='Tensile limit ~ 10 MPa')
plt.xlabel('Radius [m]')
plt.ylabel('Stress [MPa]')
plt.title('Cylindrical wall stress distribution (Moon gravity, 1 atm internal)')
plt.legend()
plt.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('cylinder_stress.png', dpi=150, bbox_inches='tight')
plt.show()
print(f"最大フープ応力: {sigma_theta.max()/1e6:.2f} MPa")
print(f"設計安全率: {10.0/sigma_max.max():.2f}")
この応力分布から、焼結レゴリスドームの構造特性が浮き彫りになります。フープ応力（周方向）が内壁で最大となり、約 $0.4\ \mathrm{MPa}$（自重寄与含めて $0.5\ \mathrm{MPa}$）。引張強度の概算 $10\ \mathrm{MPa}$ に対して安全率は20倍以上で、十分な余裕があります。しかしこれは「1気圧で与圧する」という前提のもとで、もし完全な与圧居住区にするなら、フープ応力は気圧に比例して増加します。実運用では内部にインフレータブル膜を用いて圧力を分離し、レゴリスシェルは「外圧（隕石・放射線）に対する盾」、内殻が「気密と引張負担」、と役割を分けるハイブリッドが最も合理的、という結論が数値で支持されます。
太陽光集光焼結のスループット見積もり
最後に、システム全体としてどれくらいの時間で月面基地を建設できるか、エネルギー収支から逆算します。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
P_collect = 10e3 # 集光パワー [W]
eta_total = 0.6 # 集光から焼結への総合効率
E_sinter = 1.5e6 # 1kgあたり焼結に要するエネルギー [J/kg]
duty = 0.5 # 昼間稼働率
rho_b = 2400.0 # 焼結体密度
# 1日（地球日）あたり焼結可能質量
P_eff = P_collect * eta_total * duty
m_per_day = P_eff * 86400 / E_sinter
V_per_day = m_per_day / rho_b
print(f"1日あたり焼結質量: {m_per_day:.1f} kg")
print(f"1日あたり焼結体積: {V_per_day*1000:.1f} L = {V_per_day:.3f} m^3")
# ドーム1棟(190 m^3)の建設に必要な日数
V_dome = 190.0
days_dome = V_dome / V_per_day
print(f"ドーム1棟(190 m^3)建設に: {days_dome:.0f} 日 ≈ {days_dome/30:.1f} 月")
# 集光面積を変えた感度解析
A_range = np.linspace(5, 100, 50)
days_vs_A = []
for A in A_range:
P_c = 1361.0 * A * 0.85
P_e = P_c * eta_total * duty
m_d = P_e * 86400 / E_sinter
V_d = m_d / rho_b
days_vs_A.append(V_dome / V_d)
plt.figure(figsize=(9, 5))
plt.plot(A_range, days_vs_A, 'o-')
plt.axhline(180, color='r', ls='--', label='Target: 6 months')
plt.xlabel('Solar collector area [m^2]')
plt.ylabel('Days to build one dome (190 m^3)')
plt.title('Construction time vs collector size')
plt.legend()
plt.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('build_time.png', dpi=150, bbox_inches='tight')
plt.show()
このプロットから、月面基地建設の経済合理性が見えてきます。集光面積 $10\ \mathrm{m}^2$ では1棟あたり数百日が必要ですが、$50\ \mathrm{m}^2$ まで増やすと半年程度に短縮されます。集光器自体は比較的軽量（鏡膜や薄型フレネルレンズ）で打ち上げ質量数百kg程度のため、地球からの追加質量に対するペイバックが極めて良い。これがISRUの真価です。3Dプリンタ本体（数t）と集光器（数百kg）を1度持っていけば、レゴリス資源を使って自重の数百倍以上の構造物を月面に出力できる、というシナリオが現実的な数字で支持されます。
月面水資源との結合
最後に、ISRUの拡張視点として水資源利用に触れます。LCROSS（2009年）が確認したように、月の南極の永続陰影クレーターには $\sim 5.6\ \mathrm{wt\%}$ の水氷が含まれます。これをISRUで採取し、
- 電気分解 → 液酸($\mathrm{LO_2}$) + 液水素($\mathrm{LH_2}$) 推進剤
- 居住区への生命維持水
- 放射線遮蔽（水素は遮蔽効率が高い）
として活用するのが、Artemis計画後期の中核戦略です。3Dプリンティングと水ISRUが組み合わさると、月面拠点は「自前で建物を建て、自前で推進剤を作り、地球便を最小化する」という自律性を獲得します。これは火星拠点設計の試金石でもあります。
まとめ
本記事では、月面ISRUとレゴリス3Dプリンティングを、エネルギー収支から構造設計まで定量的に解説しました。
- 質量倍率の経済学: 持参 vs 現地調達の差は1000倍規模。装置質量に対する生産量比 $R$ がISRUの成否を決め、$R \geq 100$ がコスト的に成立する目安。
- レゴリスの素性: SiO₂/Al₂O₃/FeO主体、粒度 $20$–$100\ \mu\mathrm{m}$、真空中の熱伝導率 $0.001$–$0.01\ \mathrm{W/m \cdot K}$ という極端に低い熱伝導が、局所加熱焼結を成立させる鍵。
- 4つの焼結プロセス: 太陽光（電源不要・高効率）、マイクロ波（体積加熱・FeOで吸収）、レーザー（高精度・低効率）、電子ビーム（真空適合）。用途に応じた使い分け。
- 3つの3Dプリンティング手法: Contour Crafting（大型構造）、Powder Bed Fusion（高精度部品）、SLA派生（中間ステップ）。実用月面建設はハイブリッド。
- 焼結体の物性: 多孔率と強度の関係 $\sigma_c = \sigma_{c,0}(1-\phi)^n$ から、$\phi \approx 0.2$ で圧縮強度 $100\ \mathrm{MPa}$、引張は $1/10$。アーチ・ドーム構造での圧縮主体設計が必須。
- 構造設計: 厚さ $1.3\ \mathrm{m}$ のレゴリス殻で年間被ばくを $50\ \mathrm{mSv}$ 以下に抑えられる。直径10mドームに約 $400\ \mathrm{m}^3$、集光パワー数十kWで半年程度の建設工程が見込める。
物理・材料・経済の3つの観点から、月面ISRUは単なるSF的夢物語ではなく、現実的な数値で支えられた技術的選択肢だと示せました。Artemis計画が2030年代に月面恒久拠点を実現する際、その建材は地球から運ばれた鉄ではなく、現地で焼結されたレゴリスである可能性が極めて高い、と数値が告げています。
次のステップとして、以下の記事も参考にしてください。