ロケットのフェアリングの内径は、世界最大級のFalcon 9でも約5.2 m、H3でも5.4 m程度しかありません。ところが、現代の通信衛星や地球観測衛星が必要とする太陽電池パドル(Solar Array Paddle, SAP)は、展開後に20 mから40 m級にもなります。ISSに至っては片翼で約34 m、両翼で73 mの巨大な構造です。当然、これらは打ち上げ時にはフェアリング内に折り畳まれ、軌道投入後に展開されます。
しかし展開機構は衛星設計者にとって悩ましい存在です。1990年代のJAXA技術試験衛星VIII型(ETS-VIII)では、片翼面積54 m²の大型展開アンテナ実験が行われ、世界最高クラスの展開技術が実証された一方で、海外の通信衛星でも展開不具合による「ミッション失敗」が複数報告されてきました。Hubble宇宙望遠鏡では、太陽パドルの熱歪みが姿勢制御を乱し、初期運用で大きな問題となったのは有名です。Starlink初期型では大気抵抗を最小化するため衛星全体を「シャークフィン」状に立てる運用(太陽キール)が行われましたが、これも結局はパドル形状とポインティングのトレードオフから来た苦肉の策でした。
ではなぜ、これほどリスクのある展開機構を採用するのでしょうか。答えは単純で、展開せずに必要な電力(数kW〜数十kW)を得る方法が他に存在しないからです。フェアリング体積制約と電力要求のあいだに、展開機構が架け橋として割り込んでいるのです。
本記事では、太陽電池パドルの展開機構と、そこから派生するポインティング誤差を機械系の観点から解きほぐします。応用先としては、
- 通信衛星(GEO/MEO)の大型SAP設計、SADMのジッタ予算
- 地球観測衛星(LEO)の高ポインティング要求(ピクセルサブ秒級)と熱歪みの折り合い
本記事の内容
- なぜフェアリング体積が展開機構を要求するのか、設計駆動の出発点
- 展開機構の3方式 — パッシブ(バネ式ヒンジ + ラッチ)、アクティブ(モータ駆動)、SMA(形状記憶合金)
- 多体動力学(Lagrange方程式)による2リンクパドルの展開シミュレーション
- ラッチング瞬間の衝撃トルクと、その姿勢制御系への影響
- SADM(太陽電池駆動機構)が生む角度誤差・ジッタ・バックラッシ
- 大面積パドルの片持ち梁的固有モードと姿勢制御帯域との分離設計
- 入射角の余弦則 $P=P_\mathrm{max}\cos\theta$ による発電量とポインティング精度の関係
- Pythonでの展開シミュレーション、誤差予算集計、固有モード解析、発電量プロット
前提知識
この記事を読む前に、以下の記事を読んでおくと理解が深まります。
- ニュートンの運動の3法則を数式で理解する
- 人工衛星の姿勢と座標系 — 軌道座標系から機体座標系まで
- 三軸姿勢制御 — PD制御による衛星のポインティング
- 有限要素法(FEM)入門を基礎から徹底解説
- 熱応力の理論と計算
- 座屈(オイラーの座屈荷重)
直感 — なぜ展開機構が必要か
地上の屋根に取り付ける家庭用ソーラーパネルは「畳まれず広げたまま」設置されています。地上には体積制約がほとんどないからです。一方、宇宙へ運ぶ衛星はフェアリングという「打ち上げの傘」の中に押し込まれて軌道へ送られます。Falcon 9のフェアリング内径は5.2 m、有効長は約13 m。H3-22型でも内径5.4 m、有効長は約12 mです。
電力要求はミッションによって異なりますが、典型的な目安は次の通りです。
- CubeSat(数kg): 数Wから数十W
- 地球観測小型衛星(数百kg): 数百Wから1 kW
- 通信衛星(GEO, 数トン級): 5〜20 kW
- ISSのような大型有人施設: 約100 kW
太陽電池セル(GaAs三接合, 軌道上効率約30 %)を基準にとると、軌道上での太陽放射 $1366\ \mathrm{W/m^2}$ から得られる電力密度は約 $400\ \mathrm{W/m^2}$ です。10 kW を発電するには約25 m²、20 kWでは50 m²の有効面積が必要になります。これをフェアリング内径5 m強の円柱に押し込もうとすると、当然「平面のまま」では収まりません。たとえば50 m²なら、たとえ正方形にしても7 m × 7 mで、フェアリングを大きくはみ出します。
そこで「折り畳んで運び、軌道で展開する」という発想が必然的に現れます。畳の屏風や扇のように、複数のリジッドパネル(または薄膜ブランケット)をヒンジでつなぎ、打ち上げ時には何重にも折り重ねておく。軌道に到達したらヒンジを解除して、バネやモータの力で展開する。これが太陽電池パドルの基本概念です。
折り畳めば運べる、しかし畳めるということは「動く部分がある」ということでもあります。動く部分は失敗の温床です。展開ヒンジは1回だけ正確に動けば十分(after launch)ですが、その1回が失敗すれば衛星は寿命を失います。冗長設計、地上での組合せ試験、極低温/高温下での動作試験、すべてがこの「1回」のために行われます。次節では、その1回をどんな駆動源で実現するのか — 展開機構の方式論を見ていきます。
展開機構の方式
太陽電池パドルの展開機構は、駆動源で大まかにパッシブ式、アクティブ式、形状記憶合金(SMA)式の3つに分類されます。それぞれに利点と欠点があり、衛星の規模・信頼性要求・コスト制約に応じて使い分けられます。
パッシブ式 — バネ式ヒンジとラッチ機構
パッシブ式はもっとも古典的かつ広く使われている方式で、駆動エネルギーは事前に「バネに蓄えた弾性エネルギー」のみです。代表的な構成は次の3要素から成ります。
- トーションバネ(コイルバネ): ヒンジ軸に巻き付けられ、折り畳まれた状態でねじれエネルギーを蓄えています。
- 減速機構(ダンパ): 展開速度を制限し、ラッチング瞬間の衝撃を緩和します。粘性ダンパ(オイル系)や渦電流ダンパが用いられます。
- ラッチ(エンドストッパ): 展開完了角度で機械的に「ガチッ」と固定し、剛結合状態を作る金具です。
打ち上げ時はホールドダウン機構(火工品によるピンプラ、または非火工式リリース装置)でパドルが衛星本体に縛り付けられており、軌道投入後に解放信号でリリースが切られます。バネが解放されると、ヒンジは予め設計された角度まで一気に回転し、最後にラッチへ突入して固定されます。
パッシブ式の利点は、駆動回路や電源が一切不要という点に尽きます。電気的な故障源を排除できるため信頼性が高く、CubeSatから大型衛星まで広く採用されています。欠点は、展開後に「もう一度畳む」ことができない不可逆性です。地上試験でも、毎回バネを巻き直し、ラッチをリセットしなければなりません。
アクティブ式 — モータ駆動
アクティブ式はステッピングモータやブラシレスDCモータをヒンジに組み込み、コマンドで角度制御しながら展開する方式です。展開速度を自由にコントロールでき、必要なら途中で停止することもできます。
利点は次の3点です。
- 展開シーケンスを地上から監視・調整可能(ISSのSARJなど)
- 大型構造でも安定した低速展開ができ、ラッチ衝撃を抑えられる
- 双方向制御により、軌道上で「畳み直す」(リコンフィギュア)ことも理論上可能
欠点は、モータ・減速機・電源・制御回路がすべて単一故障点(SPF, Single Point of Failure)になりうることです。冗長化のため2系統モータを設けることが一般的ですが、その分質量と複雑度が増します。Hubbleの太陽パドル交換ミッション(STS-61, 1993)では、まさにモータ駆動展開機構の不調が背景にあり、宇宙飛行士による交換修理が行われました。
SMA(形状記憶合金)式 — ヒート駆動
近年、CubeSat分野で急速に普及しているのが形状記憶合金(SMA, Shape Memory Alloy)を使ったヒンジです。代表的な材料はNiTi合金(ニチノール)で、ある変態温度($A_f$点, 典型的に40〜80°C)を超えると「記憶した形」に戻ろうとする性質を持ちます。
機構は驚くほどシンプルです。
- 製造時にSMAリンクを「展開状態」で記憶させておく。
- 折り畳んでSMAを変形させ、低温(マルテンサイト相)に保つ。
- 軌道上でヒーター(または太陽熱)でSMAを加熱、オーステナイト相へ変態。
- SMAが記憶形状に戻る復元力でパドルを展開する。
火工品が不要(衝撃ゼロ)、駆動部が少ない(故障源が少ない)、繰返し動作も可能 — というメリットがあります。一方で、展開速度の制御が難しい、加熱に時間とエネルギーが要る、低温環境では加熱しきれない、といった課題もあります。SMA式は今後10〜20年でさらに普及すると見込まれており、形状記憶高分子(SMP)による薄膜パドルの研究も進んでいます。
どの方式を選ぶにせよ、共通する設計課題は「展開時のダイナミクスをどう予測・制御するか」です。多体構造が回転・連結・衝突するため、運動方程式は非線形になり、ラッチ衝撃は瞬間的に大きなトルクを生みます。次節では、解析力学のLagrange方程式を使ってこの展開ダイナミクスを定量化します。
展開ダイナミクス — 多体動力学
太陽電池パドルを「複数の剛体リンクが回転ヒンジでつながった多体系」とモデル化することで、展開過程の力学的挙動を解析できます。最小モデルとして2リンクのパドル(本体に近い側のリンク1、外側のリンク2)を考えましょう。
系の幾何
衛星本体に固定された座標系で、ヒンジ1(本体とリンク1の接続点)の位置は原点$O$、ヒンジ2(リンク1とリンク2の接続点)の位置は$P_1$、リンク2の先端は$P_2$とします。それぞれのリンク長を$\ell_1, \ell_2$、質量を$m_1, m_2$、重心位置をリンク基準で$\ell_{c1}, \ell_{c2}$(各リンクの根本から重心までの距離)、重心まわりの慣性モーメントを$I_1, I_2$とします。
ヒンジ1の回転角を$\theta_1$(本体基準からの相対角)、ヒンジ2の回転角を$\theta_2$(リンク1基準からの相対角)とします。折り畳まれた状態を $\theta_1 = \theta_2 = 0$、完全展開状態を $\theta_1 = \theta_2 = \pi$ と取ることが多いです(実際にはミッションに応じて90°展開や180°展開など様々です)。
リンク1の重心位置は
$$ \bm{r}_{c1} = (\ell_{c1}\cos\theta_1, \ell_{c1}\sin\theta_1) $$
リンク2の重心位置は、リンク1の先端 $P_1 = (\ell_1\cos\theta_1, \ell_1\sin\theta_1)$ から測って
$$ \bm{r}_{c2} = (\ell_1\cos\theta_1 + \ell_{c2}\cos(\theta_1+\theta_2), \ \ell_1\sin\theta_1 + \ell_{c2}\sin(\theta_1+\theta_2)) $$
となります。重要なのは$\theta_2$が「絶対角」ではなく$\theta_1$からの相対角である点で、これにより各リンクの「絶対角速度」は$\dot\theta_1$と$\dot\theta_1+\dot\theta_2$になります。
運動エネルギーとLagrange方程式
各重心の速度を時間微分から求め、運動エネルギーの総和$T$を立てます。リンク1の並進運動エネルギーと回転運動エネルギーを合わせると、
$$ T_1 = \frac{1}{2}m_1\ell_{c1}^2\dot\theta_1^2 + \frac{1}{2}I_1\dot\theta_1^2 $$
リンク2は少し複雑です。並進速度の二乗を展開すると、
$$ |\dot{\bm{r}}_{c2}|^2 = \ell_1^2\dot\theta_1^2 + \ell_{c2}^2(\dot\theta_1+\dot\theta_2)^2 + 2\ell_1\ell_{c2}\dot\theta_1(\dot\theta_1+\dot\theta_2)\cos\theta_2 $$
最後の交差項に$\cos\theta_2$が現れるのは、リンク1の運動とリンク2の運動が相対角$\theta_2$でカップリングするためです。これに回転運動エネルギーを加えると、
$$ T_2 = \frac{1}{2}m_2|\dot{\bm{r}}_{c2}|^2 + \frac{1}{2}I_2(\dot\theta_1+\dot\theta_2)^2 $$
軌道上(微小重力環境)では重力ポテンシャルは無視できるため、ポテンシャルエネルギーはバネ($\theta_i$の関数)とラッチ反力のみから生じます。バネを線形回転バネ$U_\mathrm{spring} = \frac{1}{2}k_1(\theta_1 – \theta_1^\mathrm{eq})^2 + \frac{1}{2}k_2(\theta_2 – \theta_2^\mathrm{eq})^2$としても良いですが、実際の展開バネは「定トルクバネ」に近く、$\theta$によらず一定トルク$\tau_\mathrm{spring}$を出します。ここでは後者を採用し、Lagrange方程式を
$$ \frac{d}{dt}\left(\frac{\partial T}{\partial \dot\theta_i}\right) – \frac{\partial T}{\partial \theta_i} = Q_i \quad (i=1,2) $$
の形で書きます。一般化力$Q_i$には、バネトルク$\tau_{s,i}$、ダンパトルク$-c_i\dot\theta_i$、ラッチトルク$\tau_{l,i}(\theta_i)$ がすべて含まれます。
具体的に計算すると、最終的に質量行列$\bm{M}(\bm{\theta})$と非線形項(コリオリ・遠心力)$\bm{C}(\bm{\theta},\dot{\bm{\theta}})\dot{\bm{\theta}}$を使って、
$$ \bm{M}(\bm{\theta})\ddot{\bm{\theta}} + \bm{C}(\bm{\theta},\dot{\bm{\theta}})\dot{\bm{\theta}} = \bm{\tau} $$
の形に整理されます。質量行列は
$$ \bm{M}(\bm{\theta}) = \begin{pmatrix} M_{11} & M_{12} \\ M_{12} & M_{22} \end{pmatrix} $$
の対称行列で、各成分は
$$ \begin{aligned} M_{11} &= I_1 + I_2 + m_1\ell_{c1}^2 + m_2(\ell_1^2 + \ell_{c2}^2 + 2\ell_1\ell_{c2}\cos\theta_2) \\ M_{12} &= I_2 + m_2(\ell_{c2}^2 + \ell_1\ell_{c2}\cos\theta_2) \\ M_{22} &= I_2 + m_2\ell_{c2}^2 \end{aligned} $$
となります。$M_{12}$に$\cos\theta_2$が入っているため、リンク2が90°近辺($\theta_2\to\pi/2$)では結合が弱まり、180°($\theta_2\to\pi$)でリンクが一直線になると最大の負結合となります。
ラッチング衝撃トルク
展開の最終局面、$\theta_i$がエンドストッパの角度$\theta_i^\mathrm{end}$に到達した瞬間、ラッチが「ガチッ」と噛み合います。理想化すると、ラッチは$\theta_i = \theta_i^\mathrm{end}$でステップ的に剛結合となり、ヒンジは事実上「速度ゼロ」に強制されます。
この衝撃を運動量変化として扱うと、ラッチ前の角速度$\dot\theta_i^-$をゼロにするのに必要な角力積は
$$ J_i = \int_{\mathrm{latch}} \tau_l\, dt = I_\mathrm{eq,i}|\dot\theta_i^-| $$
ここで$I_\mathrm{eq,i}$はそのヒンジから見た等価慣性モーメントです。ラッチ時間$\Delta t$は数ms以下と短いため、瞬間ピークトルクは$\tau_\mathrm{peak} \sim J_i/\Delta t$と非常に大きな値になりえます。例えば$I_\mathrm{eq} = 50\ \mathrm{kg\,m^2}$、ラッチ前角速度$\dot\theta = 0.5\ \mathrm{rad/s}$、ラッチ時間$\Delta t = 2\ \mathrm{ms}$なら、瞬間トルクは
$$ \tau_\mathrm{peak} \approx \frac{50 \times 0.5}{0.002} = 12{,}500\ \mathrm{N\,m} $$
にも達します。これは衛星本体に対する反作用として戻ってくるため、姿勢制御系にはインパルス外乱となります。ダンパで角速度を絞る、エラストマー製のソフトストッパでパルス幅を伸ばす、といった対策が採られます。
ここまでで「展開がどう進み、最後にどんな衝撃を生むか」がモデル化できました。しかし展開後も問題は続きます。展開状態のパドルは、姿勢制御系から見ると「動く太陽指向ターゲット」を保持しなければならず、ここで新たな誤差源 — SADMの角度誤差や熱変形 — が顔を出します。次節ではこれらを整理します。
SADMとポインティング誤差源
展開が完了したパドルは、太陽方向を常に向けるために衛星本体に対して相対回転を続けます。地球周回衛星なら軌道周期(90分〜24時間)に同期して回り、GEO衛星では1日1回転です。この回転を担うのがSADM (Solar Array Drive Mechanism, 太陽電池駆動機構)です。
SADMの構成
SADMは典型的に次の要素で構成されます。
- ステッピングモータ — 通常2相または5相、1ステップあたり0.1〜1°の分解能
- 減速機(ハーモニックドライブ) — 高減速比(50:1〜200:1)で大トルクを得る
- スリップリング — 回転継ぎ手で電力と信号を本体側に伝送
- 角度センサ(ポテンショメータ or レゾルバ) — 現在角度のフィードバック
SADMの角度誤差源は、大きく次の4つに分解できます。
- 指令分解能誤差 — ステッピングモータの最小ステップ角(典型的に$\pm$ステップ角の半分)
- 減速機バックラッシ — ハーモニックドライブでも数十秒角〜数分角のガタが残る
- 追尾アルゴリズム誤差 — 太陽位置推算と実際の差(典型的に数分角〜0.1°)
- 熱変形 — モータ筐体や減速機の温度変化による角度ドリフト
これらを統計的に二乗和(RSS, Root-Sum-Square)で合成したものがSADMポインティング誤差予算となります。例えば各成分が0.05°なら、合計はRSSで$\sqrt{4\times 0.05^2} = 0.1°$です。
SADMが生む姿勢ジッタ
意外と見落とされがちなのが、SADMの回転自体が衛星姿勢にジッタを与える問題です。ステッピングモータは離散的にステップするため、各ステップで微小な角加速度パルスが発生します。パドル側の慣性$I_\mathrm{SAP}$と本体側の慣性$I_\mathrm{bus}$の比に応じて、本体には反作用として
$$ \Delta\omega_\mathrm{bus} = -\frac{I_\mathrm{SAP}}{I_\mathrm{bus}}\Delta\omega_\mathrm{SAP} $$
の角速度変化が生じます。地球観測衛星(典型$I_\mathrm{bus} \approx 1000\ \mathrm{kg\,m^2}$, $I_\mathrm{SAP} \approx 100\ \mathrm{kg\,m^2}$)では比が0.1程度。1ステップ$\Delta\theta_\mathrm{step} = 0.1°$をパルス的に$\Delta t = 10\ \mathrm{ms}$で実行すると、本体姿勢には
$$ \Delta\theta_\mathrm{bus} \approx -0.1 \times \frac{I_\mathrm{SAP}}{I_\mathrm{bus}} \times 0.1° = -0.001° $$
の偏りが瞬間的に発生します。$0.001° \approx 17\ \mu\mathrm{rad}$。高分解能観測(GSD 0.5 m, 軌道高度500 km)では1ピクセル $\approx 1\ \mu\mathrm{rad}$なので、SADMジッタは17ピクセル分の像ブレに相当します。これは無視できません。
ISSのBeta Gimbal(SADMの一種)は2007年に振動異常を起こし、回転動作の停止を余儀なくされました。ハーモニックドライブの歯の摩耗が原因とされ、SADMが「衛星寿命を律速する」可能性すらあることを示した事例です。
構造熱変形
太陽電池パドルは「半分が太陽光に晒され($+150°\mathrm{C}$)、半分が地球陰に入る($-100°\mathrm{C}$)」という極端な温度勾配を頻繁に受けます。日陰突入(eclipse entry)の瞬間、表裏で200°C以上の温度差が生じることもあります。
熱応力の記事で扱ったように、温度差$\Delta T$、長さ$L$、線膨張係数$\alpha$の梁の片側だけが$\Delta T$変化すると、バイメタル原理で曲率半径$\rho$が生じ、先端のたわみ$\delta$は
$$ \delta \approx \frac{\alpha\Delta T L^2}{2t} $$
($t$はパネル厚)となります。例えば$\alpha = 16 \times 10^{-6}/\mathrm{K}$(アルミ合金), $\Delta T = 100\ \mathrm{K}$, $L = 5\ \mathrm{m}$, $t = 0.02\ \mathrm{m}$ なら、
$$ \delta \approx \frac{16\times 10^{-6} \times 100 \times 25}{0.04} = 1.0\ \mathrm{m} $$
という驚くべき値になります。実際にはCFRP製ヨーク(線膨張係数 $\sim 10^{-6}/\mathrm{K}$)を採用するなどして1〜2桁低減しますが、それでも姿勢制御系から見ると「太陽光の入射角度に応じて動く軟性外乱」となります。Hubble打ち上げ初期の太陽パドル振動(orbital flutter)はまさにこの問題で、後に高安定材質に交換されました。
熱変形は時間スケールが軌道周期(数十分)とゆっくりなので、フィードフォワード補正やパッシブ熱設計(等熱膨張材の選定、熱フィルム)で対処されます。一方、次に見る柔軟構造振動はミリ秒〜秒スケールの動的現象で、姿勢制御系の帯域と干渉する可能性があります。
柔軟構造との連成振動
完全展開された太陽電池パドルは、衛星本体に「片持ち梁」のように取り付けられた巨大な薄板です。質量に対して曲げ剛性が小さいため、軌道上ではわずかな外乱で振動します。
片持ち梁の固有振動
長さ$L$、断面二次モーメント$I$、断面積$A$、密度$\rho$、ヤング率$E$の片持ち梁(Euler-Bernoulli梁)の自由振動方程式は
$$ EI\frac{\partial^4 w}{\partial x^4} + \rho A\frac{\partial^2 w}{\partial t^2} = 0 $$
で、境界条件は根元(固定端)で$w = 0, \partial w/\partial x = 0$、先端(自由端)で曲げモーメント・せん断力ゼロです。解を$w(x,t) = \phi(x)e^{i\omega t}$と仮定して代入すると固有値問題になり、無次元固有値$\beta_n L$は特性方程式
$$ \cos(\beta_n L)\cosh(\beta_n L) + 1 = 0 $$
の解として求まります。数値的に解くと$\beta_1 L \approx 1.875$, $\beta_2 L \approx 4.694$, $\beta_3 L \approx 7.855, \ldots$ となり、各モードの固有振動数は
$$ \omega_n = (\beta_n L)^2\sqrt{\frac{EI}{\rho A L^4}} $$
で与えられます。1次モードは最低周波数で、典型的なSAP(長さ10 m, アルミハニカム板)では$f_1 = \omega_1/(2\pi) \approx 0.1\sim 1\ \mathrm{Hz}$程度に来ます。これは姿勢制御の帯域(典型的に0.01〜0.1 Hz)と十分近く、制御-構造連成 (Control-Structure Interaction, CSI)問題を引き起こします。
制御-構造連成と帯域分離
姿勢制御系がパドル先端の慣性をリアルタイムで「叩く」と、パドル1次モードが励起され、本体姿勢にフィードバックされます。閉ループ系で見ると、制御帯域内に構造共振があると不安定化のリスクが高くなります。
設計の鉄則は 「制御帯域 < 構造1次振動数 / 10」 という余裕を取ることです。例えばパドル1次が0.5 Hzなら、姿勢制御帯域は0.05 Hz以下に抑える、というものです。これでは応答が遅すぎる場合は、
- ノッチフィルタ — 構造共振点だけ位相利得を下げる
- モード補償制御 — モード座標を状態変数に含めるH∞/LQG設計
- 能動制振 (active damping) — パドル根元にピエゾアクチュエータを置き、振動を能動的に減衰させる
といった対策が必要です。ETS-VIIIの大型展開アンテナはまさにこの制御-構造連成設計の試金石となりました。
ここまでで、展開機構の方式・展開時のダイナミクス・展開後の誤差源・柔軟振動 — 太陽電池パドルを取り巻く主要な機械的課題を一通り見渡しました。次節では、これらの理論をPythonで実装し、実際の数値で挙動を確認していきます。
Python実装 — 展開シミュレーション
ここからは、これまで導出した運動方程式と誤差モデルをPythonで実装し、グラフで挙動を確認します。最初に2リンクパドルの展開ダイナミクスから始めましょう。
2リンク多体動力学のシミュレーション
scipy.integrate.solve_ivpで時間積分し、ラッチ機構を「角度がエンドストッパに到達したら強い反力バネ」で近似します。
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp
# パドル物理パラメータ
m1, m2 = 5.0, 5.0 # 各リンク質量 [kg]
l1, l2 = 2.0, 2.0 # 各リンク長 [m]
lc1, lc2 = 1.0, 1.0 # 各リンク重心位置(根本から) [m]
I1, I2 = m1*l1**2/12, m2*l2**2/12 # 重心まわり慣性 [kg m^2]
# 展開機構パラメータ
tau_spring = 2.0 # バネ定トルク [N m] (ヒンジ1,2 同値)
c_damp = 0.5 # 粘性ダンパ係数 [N m s/rad]
theta_end = np.pi # ラッチ角(完全展開)
k_latch = 5000.0 # ラッチ反力バネ定数 [N m/rad]
c_latch = 20.0 # ラッチダンパ [N m s/rad]
def mass_matrix(theta2):
"""質量行列 M(theta2)"""
M11 = I1 + I2 + m1*lc1**2 + m2*(l1**2 + lc2**2 + 2*l1*lc2*np.cos(theta2))
M12 = I2 + m2*(lc2**2 + l1*lc2*np.cos(theta2))
M22 = I2 + m2*lc2**2
return np.array([[M11, M12], [M12, M22]])
def coriolis(theta2, dtheta1, dtheta2):
"""コリオリ・遠心力項 C(theta,dtheta) dtheta"""
h = -m2*l1*lc2*np.sin(theta2)
C11 = h*dtheta2
C12 = h*(dtheta1 + dtheta2)
C21 = -h*dtheta1
C22 = 0.0
return np.array([C11*dtheta1 + C12*dtheta2,
C21*dtheta1 + C22*dtheta2])
def latch_torque(theta, dtheta):
"""エンドストッパでの反力トルク(片側のみ作用)"""
if theta < theta_end:
return 0.0
return -k_latch*(theta - theta_end) - c_latch*dtheta
def dynamics(t, y):
"""状態 y = [theta1, theta2, dtheta1, dtheta2]"""
theta1, theta2, dtheta1, dtheta2 = y
M = mass_matrix(theta2)
C = coriolis(theta2, dtheta1, dtheta2)
tau1 = tau_spring - c_damp*dtheta1 + latch_torque(theta1, dtheta1)
tau2 = tau_spring - c_damp*dtheta2 + latch_torque(theta2, dtheta2)
tau = np.array([tau1, tau2])
ddtheta = np.linalg.solve(M, tau - C)
return [dtheta1, dtheta2, ddtheta[0], ddtheta[1]]
# 初期条件: 完全折畳状態(theta=0), 静止
y0 = [0.01, 0.01, 0.0, 0.0] # 微小オフセットで特異点回避
t_span = (0.0, 8.0)
t_eval = np.linspace(*t_span, 4000)
sol = solve_ivp(dynamics, t_span, y0, t_eval=t_eval, rtol=1e-8, atol=1e-10)
# プロット
fig, axes = plt.subplots(2, 1, figsize=(10, 7), sharex=True)
axes[0].plot(sol.t, np.degrees(sol.y[0]), label=r'$\theta_1$ (deg)')
axes[0].plot(sol.t, np.degrees(sol.y[1]), label=r'$\theta_2$ (deg)')
axes[0].axhline(180, color='red', ls='--', alpha=0.5, label='Latch angle')
axes[0].set_ylabel('Hinge angle [deg]')
axes[0].legend(); axes[0].grid(alpha=0.3)
axes[0].set_title('Two-link SAP deployment')
axes[1].plot(sol.t, np.degrees(sol.y[2]), label=r'$\dot\theta_1$ (deg/s)')
axes[1].plot(sol.t, np.degrees(sol.y[3]), label=r'$\dot\theta_2$ (deg/s)')
axes[1].set_xlabel('Time [s]')
axes[1].set_ylabel('Angular velocity [deg/s]')
axes[1].legend(); axes[1].grid(alpha=0.3)
plt.tight_layout()
plt.savefig('sap_deployment.png', dpi=150, bbox_inches='tight')
plt.show()
このグラフから、2リンクパドルの展開挙動の特徴的な3段階が読み取れます。第1段階は0〜3秒の加速期で、両ヒンジの角速度がほぼ線形に増加し、バネトルクとダンパ抵抗が釣り合うまで加速します。第2段階は3〜6秒の準定常期で、各ヒンジが約30 deg/sの一定角速度で展開し、$\theta_1$と$\theta_2$が180°に近づきます。第3段階はラッチング時(7秒前後)で、ラッチ反力バネにより角速度が瞬間的に振動的に止まり、角度が180°でロックされます。$M_{12}$結合のため2つのヒンジは独立には動かず、$\theta_2$の挙動が$\theta_1$にもわずかな揺れとして影響していることがわかります。
ラッチング衝撃トルクの可視化
次に、ラッチ瞬間に発生するトルクパルスを可視化し、衛星本体への外乱としての大きさを評価します。
import numpy as np
import matplotlib.pyplot as plt
# 上のsolutionから各時刻のラッチトルクを再計算
latch_tau1 = np.array([latch_torque(sol.y[0,i], sol.y[2,i]) for i in range(len(sol.t))])
latch_tau2 = np.array([latch_torque(sol.y[1,i], sol.y[3,i]) for i in range(len(sol.t))])
fig, ax = plt.subplots(figsize=(10, 4))
ax.plot(sol.t, latch_tau1, label='Latch torque, hinge 1')
ax.plot(sol.t, latch_tau2, label='Latch torque, hinge 2')
ax.set_xlabel('Time [s]')
ax.set_ylabel('Latch torque [N m]')
ax.set_title('Latching impulse torque on satellite bus')
ax.grid(alpha=0.3); ax.legend()
plt.tight_layout()
plt.savefig('latch_torque.png', dpi=150, bbox_inches='tight')
plt.show()
# 衝撃トルクのピーク値
print(f"Peak latch torque (hinge 1): {np.min(latch_tau1):.1f} N m")
print(f"Peak latch torque (hinge 2): {np.min(latch_tau2):.1f} N m")
この時系列プロットでは、展開中はラッチトルクがほぼゼロで、エンドストッパ突入の瞬間に負方向の鋭いパルスが立ち上がり、その後数百ミリ秒の振動を経て減衰しているのが見えます。ピーク値は数百N mのオーダーで、これがそのまま衛星本体への反作用となります。本体慣性$I_\mathrm{bus} = 500\ \mathrm{kg\,m^2}$なら$500/500 = 1\ \mathrm{rad/s^2}$の瞬間角加速度を本体に与えうるため、ホイールアンローディングやスラスタ補正で吸収する必要があります。
片持ち梁の固有モード解析
展開後のパドルを片持ち梁として、1次〜3次の振動モードを数値的に求めます。
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import brentq
# 特性方程式 cos(beta L) cosh(beta L) + 1 = 0
def char_eq(x):
return np.cos(x)*np.cosh(x) + 1.0
# 1次〜5次の固有値 (beta_n L)
roots = []
intervals = [(1.5, 2.5), (4.0, 5.0), (7.0, 8.5), (10.0, 11.5), (13.5, 14.5)]
for a, b in intervals:
roots.append(brentq(char_eq, a, b))
roots = np.array(roots)
print("beta_n L =", roots)
# パドルパラメータ(典型例: アルミハニカム両面板)
L = 10.0 # 長さ [m]
E = 70e9 # ヤング率 [Pa] (Al合金)
rho = 50.0 # 等価密度 [kg/m^3] (ハニカムは軽い)
A = 1.0 # 幅 1 m, 厚さ0.05 m を仮定
I = 1.0 * 0.05**3 / 12 # 断面二次モーメント [m^4]
# 固有振動数
omegas = roots**2 / L**2 * np.sqrt(E*I / (rho*A))
freqs = omegas / (2*np.pi)
for n, f in enumerate(freqs, start=1):
print(f"Mode {n}: f = {f:.3f} Hz")
# モード形状
x = np.linspace(0, L, 200)
plt.figure(figsize=(10, 5))
for n, beta_L in enumerate(roots[:4], start=1):
beta = beta_L / L
sigma = (np.sin(beta_L) + np.sinh(beta_L)) / (np.cos(beta_L) + np.cosh(beta_L))
phi = np.cosh(beta*x) - np.cos(beta*x) - sigma*(np.sinh(beta*x) - np.sin(beta*x))
phi /= np.max(np.abs(phi))
plt.plot(x, phi, label=f'Mode {n}, f={freqs[n-1]:.2f} Hz')
plt.xlabel('Position along paddle x [m]')
plt.ylabel('Mode shape (normalized)')
plt.title('Cantilever mode shapes of solar paddle')
plt.legend(); plt.grid(alpha=0.3)
plt.tight_layout()
plt.savefig('paddle_modes.png', dpi=150, bbox_inches='tight')
plt.show()
このモード形状プロットから、典型的な片持ち梁モードの特徴が確認できます。1次モードは単純なC字曲げで、節点なし、先端で最大変位。2次モードは1個の節点を持つS字、3次モードは2個の節点、というように、高次になるほど節点が増え、変位の符号変化が増えます。固有振動数は約0.07 Hz, 0.42 Hz, 1.18 Hz, …と$\beta_n^2$に比例して急速に上昇します。1次モードが0.07 Hzと姿勢制御帯域に近いのが見て取れ、CSI問題の典型的な数値感覚が掴めます。
ポインティング誤差予算の集計
複数の誤差源をRSS合成し、システムレベルのポインティング誤差を試算します。
import numpy as np
import matplotlib.pyplot as plt
# 各誤差源の3シグマ値 [deg]
errors = {
'SADM step resolution': 0.05,
'Harmonic drive backlash': 0.03,
'Sun sensor noise': 0.02,
'Thermal deformation (CFRP yoke)': 0.04,
'Attitude control bias': 0.05,
'Structural flexure (1st mode)': 0.03,
'Alignment / manufacturing': 0.02,
}
# RSS合成
sigma_sum = np.sqrt(sum(v**2 for v in errors.values()))
# 寄与度プロット
fig, ax = plt.subplots(figsize=(10, 5))
names = list(errors.keys())
values = np.array(list(errors.values()))
contrib = values**2 / np.sum(values**2) * 100
bars = ax.barh(names, values, color='steelblue')
ax.set_xlabel('3-sigma error [deg]')
ax.set_title(f'Pointing error budget (RSS total = {sigma_sum:.3f} deg)')
ax.grid(axis='x', alpha=0.3)
for i, (b, c) in enumerate(zip(bars, contrib)):
ax.text(b.get_width()+0.001, i, f' {c:.0f}%', va='center')
plt.tight_layout()
plt.savefig('error_budget.png', dpi=150, bbox_inches='tight')
plt.show()
print(f"Total RSS pointing error: {sigma_sum:.3f} deg = {sigma_sum*60:.1f} arcmin")
この棒グラフから、ポインティング誤差予算の構造がはっきり見えます。SADM分解能と姿勢制御バイアスがそれぞれ25 %前後の寄与で支配的で、熱変形がそれに続きます。一方、アライメント誤差や柔軟構造の振動寄与は10 %以下に抑えられています。RSS合計は約0.11°(6.6 arcmin)で、典型的なGEO通信衛星(要求0.05°級)では足りませんが、地球観測ミッション(要求0.1°級)なら満たせる水準です。要求精度を満たすには、寄与の大きいSADM分解能と姿勢制御バイアスを優先的に低減すべきと示唆されます。これがエラーバジェットの実践的な使い方です。
発電性能とポインティング精度
最後に、ポインティング誤差が発電量にどう影響するかを定量化します。
余弦則と発電量
太陽電池に対する太陽光の入射角度を$\theta$(法線からの傾き)とすると、単位面積あたりの吸収パワーは余弦則に従って
$$ P = P_\mathrm{max}\cos\theta $$
となります。これは、面積$A$の平面に対し、法線方向から見た投影面積が$A\cos\theta$になるためです。
ここで重要なのは、$\theta$が小さい範囲では損失は二次のオーダーということです。Taylor展開すると、
$$ P/P_\mathrm{max} = \cos\theta \approx 1 – \frac{\theta^2}{2} $$
なので、$\theta = 5°$($\approx 0.087$ rad)でも損失は$0.087^2/2 \approx 0.38$ %にすぎません。これがSADM追尾精度を「分角」レベルで十分とできる理由です。一方、$\theta = 30°$なら$\cos 30° = 0.866$で13 %の損失、$60°$なら50 %損失と急速に悪化します。
Pythonで発電量と誤差を可視化
import numpy as np
import matplotlib.pyplot as plt
# 入射角の範囲
theta_deg = np.linspace(0, 90, 200)
theta_rad = np.radians(theta_deg)
power_ratio = np.cos(theta_rad)
# 季節変動による太陽-軌道角(beta angle)の効果
# 典型LEO: beta = -75° から +75° で変動
beta_angles_deg = np.linspace(-75, 75, 100)
beta_rad = np.radians(beta_angles_deg)
# SADM が1軸追尾なら、残りの軸誤差として |sin(beta)| が残る
power_with_1axis_sadm = np.cos(beta_rad) # 単純化モデル
fig, axes = plt.subplots(1, 2, figsize=(13, 5))
axes[0].plot(theta_deg, power_ratio*100, lw=2)
axes[0].axvline(5, color='green', ls='--', alpha=0.6, label='5 deg (loss<1%)')
axes[0].axvline(30, color='orange', ls='--', alpha=0.6, label='30 deg (~13% loss)')
axes[0].axvline(60, color='red', ls='--', alpha=0.6, label='60 deg (50% loss)')
axes[0].set_xlabel('Pointing error theta [deg]')
axes[0].set_ylabel('Power ratio [%]')
axes[0].set_title('Power vs pointing angle (cos rule)')
axes[0].legend(); axes[0].grid(alpha=0.3)
axes[1].plot(beta_angles_deg, power_with_1axis_sadm*100, lw=2)
axes[1].set_xlabel('Beta angle [deg]')
axes[1].set_ylabel('Power ratio [%]')
axes[1].set_title('Seasonal variation (1-axis SADM)')
axes[1].grid(alpha=0.3)
axes[1].axhline(np.cos(np.radians(75))*100, color='red', ls='--', alpha=0.5,
label=f'min at |beta|=75: {np.cos(np.radians(75))*100:.0f}%')
axes[1].legend()
plt.tight_layout()
plt.savefig('power_pointing.png', dpi=150, bbox_inches='tight')
plt.show()
左のプロットから、ポインティング誤差が10°以下なら発電損失は数%以下に抑えられ、SADM追尾の要求精度(典型的に1°級)なら損失はほぼ無視できるレベル(0.02 %以下)と読み取れます。これは、SADMの設計目標を「太陽追尾精度よりも、本体姿勢への外乱抑制(ジッタ低減)に振る」設計判断の正当化につながります。
右のプロットは、太陽-軌道面角度(beta angle, 季節変動)による発電量の変化を示しています。1軸SADM(衛星本体に対し1軸のみ回転)では、beta角が大きくなると太陽方向と回転軸面のなす角がそのまま残るため、cos則で発電量が減衰します。beta = $\pm 75°$では$\cos 75° \approx 0.26$、つまり通常時の26 %まで発電量が落ちる可能性があります。これがGEO衛星やSSO(太陽同期軌道)衛星で、季節ごとに電力収支を再計算する理由です。LEO衛星では2軸SADM(片方は本体内、もう片方はヨーク先で回転)を採用してこの変動を打ち消す例もあります。
まとめ
本記事では、太陽電池パドルの展開機構とポインティング誤差を、構造設計の視点から包括的に解説しました。
- フェアリング体積制約が展開機構の根本的な存在理由 — 内径5 m強の打ち上げ機で20 m級のパドルを運ぶには、折り畳んで展開する以外に手段がない
- 展開機構の3方式 — パッシブ式(バネ+ダンパ+ラッチ)、アクティブ式(モータ駆動)、SMA式(形状記憶合金) — それぞれ信頼性・制御性・コストのトレードオフ
- 展開ダイナミクスは多体系のLagrange方程式で記述され、質量行列が$\cos\theta_2$を介してリンク間結合を持つ。ラッチング時には$I\dot\theta/\Delta t$オーダーの瞬間トルクパルスが衛星本体に作用する
- SADMはステッピングモータと減速機の組合せで、太陽追尾精度(数分角)よりも回転動作が衛星姿勢に与えるジッタが問題になりやすい
- 熱変形は表裏100°C以上の温度差で先端が数十cm〜数mたわむ可能性があり、CFRP等の低熱膨張材で1〜2桁低減する
- 柔軟構造振動は1次モードで0.1〜1 Hzに来るため、制御帯域は構造1次の1/10以下に抑える「帯域分離」が原則。難しい場合はノッチフィルタや能動制振が必要
- ポインティング誤差予算はSADM分解能・バックラッシ・熱変形・姿勢制御バイアスのRSSで合成され、GEO通信衛星0.05°級、地球観測0.1°級が一般的な目標
- 発電量と入射角は余弦則 $P=P_\mathrm{max}\cos\theta$ に従い、誤差5°以下なら損失は1%未満。一方、季節変動(beta角変化)では1軸SADMでは$\cos\beta$分の発電低下が避けられない
太陽電池パドルは「広げないと発電できないが、広げると揺れて姿勢を乱す」という構造設計者にとって根本的なジレンマを体現した存在です。展開機構の選定、SADMの分解能とジッタのトレードオフ、熱設計、柔軟構造との連成 — これらをシステムレベルで折り合わせるのが衛星構造設計の核心です。
次のステップとして、以下の記事も参考にしてください。