Sパラメータ解析
理論と物理
Sパラメータとは
先生、Sパラメータって何を表す量ですか？
高周波回路の入出力特性を反射波と透過波の比で表す。2ポートの場合：
$a_i$: 入射波、$b_i$: 反射波。$S_{11}$: 反射係数、$S_{21}$: 透過係数。
$|S_{11}|$が小さいほどインピーダンス整合が良いんですね。
そう。$S_{11} = -20$ dBなら反射電力は1%。$S_{21} = -3$ dBなら透過電力は半分（3 dBロス）。Sパラメータは周波数の関数で、VNA（ベクトルネットワークアナライザ）で測定する。
まとめ
Sパラメータの誕生——散乱行列が「ポート間関係」を変えた
Sパラメータ（散乱パラメータ）の概念はK. KurokiとD. M. Pozarらによって整備され、ZパラメータやYパラメータでは扱いが難しかった「マイクロ波回路の入射波・反射波・透過波の関係」を統一的に記述した。特に伝送線路上で測定可能な「反射係数Γ（=S₁₁）」「透過係数S₂₁」は実験と理論を直結させ、ネットワークアナライザの普及とともに高周波設計の共通言語となった。CAEでのSパラメータ計算は固有モード展開→ポートモード規格化のプロセスで実現される。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ（発電機）は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波（GHz帯のアンテナ等）では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$：電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き（電荷）から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ（キャパシタ）の設計では、電極間の電場分布をこの法則で計算する。ESD（静電気放電）対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$：磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定：透磁率・誘電率が磁場・電場強度に依存しない（飽和領域では非線形B-Hカーブが必要）
- 準静的近似（低周波）：変位電流項を無視可能（$\omega \varepsilon \ll \sigma$）。渦電流解析で一般的
- 2D仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
FEMでのSパラメータ抽出
FEMからSパラメータをどう抽出しますか？
1. ポートにモードパターン（TE10等）を設定
2. 1つのポートから入射波を励振
3. 各ポートの反射波・透過波を計算
4. $S_{ij} = b_i/a_j$で算出
HFSSの適応メッシュは$\Delta S$（Sパラメータの変化量）を収束判定に使用。
マルチポートの場合は？
$N$ポートなら$N \times N$のSマトリクス。各ポートから順に励振して$N$回解く（Direct Solver）。または全ポート同時に解いて行列演算で$S$を抽出（Fast Frequency Sweep）。
まとめ
De-embedding——「測定値からポートの影響を除去する」技術
実際のVNAでSパラメータを測定するとき、コネクタ・ケーブル・治具の影響が混入する。これを除去する「De-embedding（脱埋め込み）」は、既知の参照構造（Open/Short/Thru）を測定してポートモデルを逆行列演算で除くプロセスだ。CAEでも同様のDe-embedding概念が使われ、解析領域内のポート波形から実際のデバイス特性を切り出す。TRL（Thru-Reflect-Line）校正がVNAでの標準手法であり、CST・HFSSのポート設定はこの手続きを自動化している。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務での活用
フィルタ設計、コネクタ評価、パッケージのSI解析、アンテナ給電回路の設計が代表的。
実務チェックリスト
「S₁₁が改善したらS₂₁が悪化した」——整合設計のトレードオフ
アンテナや増幅器の入力整合設計では「S₁₁（反射損）を下げると帯域が狭くなる」フロービッツの限界（Bode-Fano限界）がある。これは「反射係数の周波数積分は一定値を超えられない」という理論的制約で、反射を特定帯域で徹底的に下げれば他の帯域で必ず悪化する。CAEでは複数周波数点でSパラメータを最適化する際にこの制約を暗黙に踏んでしまうことが多い。最適化前にBode-Fano限界から「達成可能な帯域幅とS₁₁の組み合わせ」を確認する習慣が設計効率を上げる。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| Ansys HFSS | Sパラメータの業界標準。Touchstone出力 |
| CST Studio Suite | 時間/周波数領域の両方でSパラ計算 |
| Keysight ADS | 回路+電磁界の統合。Sパラのデエンベッディング |
| Cadence Clarity 3D | PCB/パッケージのSパラメータ抽出 |
Sパラメータ解析ツール——Keysight ADS vs ANSYS HFSS
Sパラメータ解析の専用ツールとしてKeysight ADS（Advanced Design System）とANSYS HFSSが代表格だ。ADSは回路レベルのSパラメータ解析・マッチング回路設計・EMCoSimulation（EM/回路協調）が得意で、RF回路設計のデファクト。HFSSは3D構造のポート解析・固有モード抽出・アダプティブメッシュでSパラメータを高精度に算出し、アンテナ・導波管・RFICパッケージに強い。Touchstoneファイル（.s2p/.snp）を介した連携が業界標準で、両ツールは補完関係にある。
選定で最も重要な3つの問い
- 「何を解くか」：Sパラメータ解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術
Sパラメータの「因果性」——パッシブ・因果・安定の三条件
Sパラメータが物理的に実現可能か判定する条件として「パッシブ性（エネルギー非増幅）・因果性（未来の入力に応答しない）・安定性」の三つがある。CAEで計算したSパラメータがこれらを満たさないとき、回路シミュレータに取り込んだ際に非因果的な過渡応答や「負の損失」が生じる。強制エンフォースメント法（パッシビティ保正アルゴリズム）で修正する手法が標準的で、ANSYS SIwaveやKeysight ADS内にこの機能が実装されている。測定精度と因果性違反の鑑別が実務での重要スキルだ。
トラブルシューティング
トラブル
「Sパラメータがポートモードで変わる」——ポート定義の落とし穴
FEM/FDTDシミュレーションでSパラメータを計算する際、ポートの「モード規格化インピーダンス」設定によって値が変わるという混乱が多い。50 Ω規格化と100 Ω規格化ではS₁₁の値が異なり、どちらが正しいかはシステム仕様に依存する。HFSS・CSTでは「Renormalization（規格化変換）」機能でポート後から変換が可能だが、多ポート問題では行列変換が必要だ。トラブル回避の鉄則は「CAEのポートインピーダンスとVNAの規格化インピーダンスを統一する」こと、これだけで測定・解析の乖離の8割は防げる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Sパラメータ解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告