フェーズドアレイアンテナのCAEシミュレーション
理論と物理
概要 — フェーズドアレイとは
フェーズドアレイって戦闘機のレーダーに使われてるやつですか？
そう、F-35に搭載されているAN/APG-81のようなAESA（Active Electronically Scanned Array）レーダーが代表例だ。最近は5G基地局やStarlinkの衛星端末にも搭載されている。各素子の位相・振幅を電子的に制御してビームをミリ秒単位でスキャンできるのが最大の特徴だよ。
機械的にアンテナを回さなくていいってことですか？ それはすごいですね。
そういうこと。従来の回転式パラボラレーダーだと、ビーム走査は機械回転のスピードに制限される。フェーズドアレイなら電子制御だから、1秒間に数千回のビーム切り替えが可能だ。同時に複数の目標を追跡したり、レーダーと通信を切り替えたりもできる。イージス艦のSPY-1レーダーは約4,400素子で全方位をカバーしているし、5Gの基地局は64〜256素子のMassive MIMOで都市部のユーザーに個別ビームを向けている。
で、そのフェーズドアレイをCAEでシミュレーションするってことですよね。何がキーポイントになるんですか？
核心は3つだ。(1) アレイファクタ——素子配置と位相差が決めるビームの指向性パターン、(2) 素子パターン——個々の素子の放射特性、(3) 相互結合——隣接素子間の電磁的干渉。この3つを正確にモデリングできるかどうかで、シミュレーションの信頼性が決まる。
アレイファクタ（AF）
アレイファクタの数式を教えてもらえますか？
N個の等間隔素子からなるリニアアレイ（1次元配列）のアレイファクタはこう表される。
各変数の意味を整理しよう。
- $N$ — 素子数
- $a_n$ — 第 $n$ 素子の励振振幅（テーパリング重み）
- $k = 2\pi/\lambda$ — 自由空間の波数
- $d$ — 素子間隔
- $\theta$ — アレイ法線からの角度
- $\beta$ — 隣接素子間の位相差（ビームステアリング用）
等振幅 ($a_n = 1$) の場合は閉じた形に書けたりしますか？
いい質問だ。等方等振幅のユニフォームアレイなら幾何級数の公式で閉じて、
この式を見ると、$\psi = 0$ のとき $AF = N$ で最大値を取る。つまりメインビームは $\sin\theta_0 = -\beta / (kd)$ の方向に向く。ヌル（ゼロ）は $\psi = 2m\pi/N$ ($m \neq 0, N, 2N, \ldots$) で現れる。
素子パターンとパターン乗算の原理
アレイファクタだけで放射パターン全体が決まるんですか？
いや、それだけでは不十分だ。各素子自体にも固有の放射パターン（素子パターン）がある。全体の放射パターンはパターン乗算の原理で決まる。
例えばダイポール素子を使うアレイなら、ダイポール固有の $\cos\theta$ パターンがAFに掛け合わされる。パッチアンテナ素子なら $\cos^n\theta$ のようなパターンになる。素子パターンのエンベロープがAFのピークを包むので、例えばアレイファクタのグレーティングローブの位置が素子パターンのヌルに重なれば、事実上グレーティングローブが抑制されるケースもあるんだ。
つまり素子の選び方もアレイ全体の性能に大きく影響するってことですね。
その通り。だからCAEシミュレーションでは、まず単素子の放射パターンを正確に求めて、それからアレイ全体の解析に進むという2段階アプローチが基本になる。
ビームステアリングの原理
位相差 $\beta$ を変えるだけでビームの向きが変わるんですか？ 直感的に理解したいです。
波面のイメージで考えてみよう。全素子が同位相 ($\beta = 0$) で放射すると、波面はアレイ面に平行——つまり正面（ブロードサイド）方向にビームが出る。ここで隣接素子ごとに少しずつ位相を遅らせると、波面が傾く。これがビームステアリングだ。
ビームを $\theta_0$ 方向に向けたければ、上の式で $\beta$ を設定すればいい。例えば $d = \lambda/2$、$\theta_0 = 30°$ なら $\beta = -\pi \sin(30°) = -\pi/2 \approx -90°$。各素子を $90°$ ずつ遅らせるだけでビームが $30°$ に傾く。
めちゃくちゃシンプルですね。でもどの角度まで振れるんですか？
理論上は $\pm 90°$ まで振れるけど、実用上は $\pm 60°$ 程度が限界だ。理由は2つ。(1) ビーム幅が $1/\cos\theta_0$ に比例して広がり、利得が低下する（$\cos\theta_0$ ファクタ）。(2) 大角度ではグレーティングローブやスキャンブラインドネスが発生しやすくなる。
グレーティングローブ条件
グレーティングローブってよく聞くんですけど、なぜ $d < \lambda/2$ が必要なんですか？
グレーティングローブは、メインビーム以外にメインローブと同等の強度を持つ放射が現れる現象だ。光学の回折格子（diffraction grating）と同じ原理で発生する。
グレーティングローブが発生しない条件は、ビームを $\theta_0$ 方向にステアリングする場合、
ブロードサイド（$\theta_0 = 0°$）のみを考えるなら $d < \lambda$ だけど、全方向にステアリングしたい場合（$\theta_0 = 90°$）は $d < \lambda / 2$ が必要になる。実務では $d = 0.5\lambda$ を超えないことが大原則だ。
5Gのミリ波帯（28 GHz）だと $\lambda \approx 10.7$ mmだから、素子間隔は5 mm程度ですか。かなり密ですね。
そう。ミリ波帯では物理的なスペースの制約が厳しくなるから、TR（送受信）モジュールの小型化とともに、素子間隔の設計が極めて重要になる。逆にいえば、VHF帯（数百MHz）のレーダーだと $\lambda$ が1 m前後になるので、アレイ全体が数十メートル規模になる。
サイドローブ制御 — テーパリング
ユニフォームアレイだとサイドローブが高くなりますよね。どうやって抑えるんですか？
振幅テーパリング（amplitude tapering）を使う。中央の素子の励振を強く、端の素子を弱くすることでサイドローブを抑制する。代表的な窓関数とサイドローブレベル（SLL）の関係をまとめよう。
| 窓関数 | 第1サイドローブ (dB) | ビーム幅の増加 | 用途 |
|---|---|---|---|
| ユニフォーム | -13.3 | 基準 | 最大利得が必要な場合 |
| コサイン | -23 | ×1.36 | 中程度の抑圧 |
| ハミング | -42.8 | ×1.50 | 通信用アンテナ |
| テイラー ($\bar{n}$) | 任意設定可 | 任意 | レーダー（SLL/BW制御） |
| チェビシェフ | 等リップル | 任意 | 最小ビーム幅でSLLを保証 |
| ドルフ-チェビシェフ | -20〜-40 | ×1.1〜1.4 | 最適トレードオフ |
サイドローブを下げるとメインビームが太くなるんですね。トレードオフがある。
そうだ。これはフーリエ変換の不確定性原理と同じで、空間ドメインで窓を絞ると角度ドメインのサイドローブは下がるがメインローブは広がる。レーダーでは Taylor テーパリングがよく使われる。「最初の $\bar{n}$ 個のサイドローブを等しいレベルに保ち、残りを $1/u$ で減衰させる」という設計ができて、実用的なバランスが取れるんだ。
2Dプレーナアレイへの拡張
リニアアレイの理論を2次元に拡張すると、$M \times N$ の矩形平面アレイのアレイファクタは
$$ AF(\theta,\phi) = \sum_{m=0}^{M-1}\sum_{n=0}^{N-1} a_{mn}\, e^{\,j\,(m\psi_x + n\psi_y)} $$ここで $\psi_x = k d_x \sin\theta\cos\phi + \beta_x$, $\psi_y = k d_y \sin\theta\sin\phi + \beta_y$。分離可能な場合（$a_{mn} = a_m^{(x)} a_n^{(y)}$）は $AF = AF_x \cdot AF_y$ と書ける。実際の5G基地局アンテナ（例えば $8 \times 8$ パネル）はこの2D AFでビームを仰角・方位角の両方に制御する。
利得とディレクティビティ
等方性素子のユニフォームリニアアレイの最大ディレクティビティは
$$ D_{\max} = N \cdot \frac{d}{\lambda} \cdot 2 \quad (\text{ブロードサイド, } d = \lambda/2 \text{ のとき } D = N) $$$M \times N$ のプレーナアレイでは $D \approx \pi M N (d_x d_y / \lambda^2) \cdot 4$ と近似できる。ただし素子パターンの影響やエッジ効果で実効的な利得は理論値より2〜3 dB低くなることが多い。
次元解析と単位系
| 変数 | SI単位 | 備考 |
|---|---|---|
| 波長 $\lambda$ | m | $\lambda = c/f$。28 GHzで10.7 mm |
| 波数 $k$ | rad/m | $k = 2\pi/\lambda$ |
| 素子間隔 $d$ | m | $d/\lambda$ で正規化して議論 |
| 位相差 $\beta$ | rad | 度数法との混同に注意 |
| ディレクティビティ $D$ | 無次元 / dBi | $D_{\mathrm{dBi}} = 10\log_{10}(D)$ |
| 電界強度 $E$ | V/m | 遠方界: $1/r$ に比例して減衰 |
フェーズドアレイの起源 — 戦時中のレーダー競争
フェーズドアレイの概念は第二次世界大戦中に生まれた。ドイツのMAMMUT（マンモート）レーダーやイギリスのChain Homeが初期のアレイアンテナシステムだ。「機械的にアンテナを回すと遅すぎる」という切実な軍事ニーズが、位相制御ビームスキャニングの理論を一気に発展させた。冷戦期のBMEWS（弾道ミサイル早期警戒システム）でAN/FPS-115 PAVEPAWSが実用化され、現代のAESAレーダーにつながっている。アレイファクタの数式自体は平面波の重ね合わせという極めてシンプルな物理で、この同じ数式が70年後の5Gスマートフォンやテスラのオートパイロットレーダーにも使われている。
数値解法と実装
電磁界解析手法の比較
フェーズドアレイの電磁界シミュレーションにはどんな手法が使われるんですか？
主に4つの手法がある。それぞれ得意分野が違うから、問題に応じて使い分けることが重要だ。
| 手法 | 略称 | 支配方程式 | 得意な問題 | 代表ソフト |
|---|---|---|---|---|
| 有限要素法 | FEM | Maxwell方程式の弱形式 | 複雑形状・不均質媒質 | Ansys HFSS |
| 時間領域差分法 | FDTD | Maxwell curl方程式のYee格子離散化 | 広帯域・過渡解析 | CST, OpenEMS |
| モーメント法 | MoM | 積分方程式 (EFIE/MFIE) | 金属構造・開領域 | FEKO, ADS Momentum |
| 物理光学法 | PO/GO | 高周波近似 | 電気的に大きな構造 | Ansys SBR+, GRASP |
フルウェーブ解析のMaxwell方程式はどういう形になるんですか？
時間調和場（$e^{j\omega t}$）を仮定すると、Maxwell方程式からベクトル波動方程式が導かれる。
FEMではこれをGalerkin法で弱形式化し、辺要素（Nedelec要素）で離散化する。辺要素を使う理由は、ベクトル場の接線連続性を自動的に保証し、スプリアス（非物理的な偽の解）モードを排除できるからだ。節点要素で電界ベクトルを近似するとスプリアス解が大量に混入して使い物にならない。
無限アレイ近似（Floquet境界）
数百素子のアレイを全部モデリングしたら計算が爆発しそうですけど、どうやって対処するんですか？
非常にいいポイントだ。大規模アレイのシミュレーションで最もよく使われるのが無限アレイ近似だ。周期的な構造の1ユニットセル（1素子分）だけをモデリングし、周囲にFloquet（フロッケ）周期境界条件を適用する。
Floquet境界条件は、隣接ユニットセルとの電界の関係を次のように規定する。
ここで $\mathbf{d}$ は周期ベクトル、$\mathbf{k}_t$ は入射波の接線方向波数ベクトル（スキャン角度に対応）。1素子分の計算コストでスキャン角度ごとのアクティブインピーダンス、アクティブ素子パターン、反射係数 $S_{11}$ が得られる。HFSSの「Floquet Port」やCSTの「Unit Cell Boundary」がまさにこの機能だ。
でも無限アレイだとエッジ効果（端の素子の影響）は分からないですよね？
その通り。端から3〜5素子程度の範囲ではアクティブインピーダンスが中央部と異なる。この影響を把握するには有限アレイの解析が必要になる。実務では、まず無限アレイで素子設計を最適化してから、有限アレイで全体性能を検証する2段階アプローチが標準だ。
有限アレイのフルウェーブ解析
有限アレイの場合はどうやって計算コストを抑えるんですか？
いくつかのアプローチがある。
- ドメイン分割法（DDM）：アレイを素子ごとにサブドメインに分割し、各サブドメインを独立にFEMで解いた後、境界面で結合する。HFSSの「3D Component DDM」が典型。
- ハイブリッドFEM-BI法：素子近傍はFEMで精密に解き、遠方界と素子間の結合は境界積分（BI）で接続する。
- MLFMA加速MoM：マルチレベル高速多極子法で$O(N^2)$のMoMを$O(N\log N)$に加速。FEKOが強い。
- 埋め込み素子パターン法：有限アレイの中央素子のパターンを「アクティブ素子パターン」として抽出し、AFと合成する準解析的手法。
DDMだと各サブドメインを並列計算できるからHPCと相性が良さそうですね。
そうだ。HFSSのDDMはMPI並列でスケールするので、64素子アレイを64コアで計算すれば、理想的には1素子分の時間で全体が解ける。ただしサブドメイン間の結合反復が必要なので、実際にはその1.5〜3倍程度の時間がかかる。
素子間相互結合のモデリング
相互結合ってそんなに大きな問題なんですか？
大問題だ。素子間隔が $\lambda/2$ 程度だと、隣接素子間の結合係数 $S_{21}$ は $-10$ 〜 $-15$ dB程度になることが多い。これを無視すると、
- アクティブインピーダンスがスキャン角度で大きく変動する
- 放射パターンに予想外の歪みが生じる
- ヌル方向がずれる
- 最悪の場合、スキャンブラインドネス（特定角度で放射効率がゼロに近づく）が発生する
相互結合の影響はSパラメータ行列で定量化する。$N$ 素子アレイでは $N \times N$ のフルSパラメータ行列を求め、各素子のアクティブ反射係数を以下で計算する。
ここで $a_m$, $\phi_m$ は第 $m$ 素子の励振振幅と位相。スキャン角度を変えると $\phi_m$ が変わるので、$\Gamma_n^{\mathrm{active}}$ も変化する。ある角度で $|\Gamma_n^{\mathrm{active}}| \to 1$ に近づくと、ほぼ全反射になってスキャンブラインドネスが起きるんだ。
メッシュ戦略
アンテナ解析のメッシュって、構造解析とは勝手が違いそうですね。
全然違う。高周波電磁界解析のメッシュは波長に対するサイズが重要になる。基本ルールをまとめよう。
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| 最大要素サイズ | $\lambda / 6$ 〜 $\lambda / 10$ | $\lambda/10$ が安全側の基準 |
| 金属エッジ付近 | $\lambda / 20$ 以下 | 回折・電流集中領域 |
| 給電構造（via, slot） | $\lambda / 30$ 以下 | 局所的な電界集中 |
| 誘電体基板厚さ方向 | 3層以上 | 表面波の捕捉に必要 |
| Radiation Boundary距離 | $\lambda / 4$ 以上 | 放射面から境界までの距離 |
| PML層数 | 5〜8層 | 多いほど反射が減るが計算コスト増 |
$\lambda/10$ って28 GHzだと約1 mmですよね。かなり細かいメッシュが必要ですね。
そう。だからHFSSでは適応メッシュリファインメント（Adaptive Mesh Refinement）が標準で、Sパラメータの変化が収束するまで自動的にメッシュを細分化する。最初は粗いメッシュで解いて、電界勾配が大きい箇所を自動で特定・細分化していくから、人手でメッシュサイズを指定するよりずっと効率的だ。
周波数領域と時間領域の使い分け
FEM（HFSS）は周波数ごとに1回の行列解法が必要で、狭帯域の詳細解析に向く。一方FDTD（CST）は1回の時間進行シミュレーションで広帯域のSパラメータが得られるが、高いQ値の共振構造では収束に時間がかかる。フェーズドアレイでは「中心周波数でのパターン最適化」はFEM、「広帯域マッチング確認」はFDTDと使い分けるのが効率的だ。
実践ガイド
解析ワークフロー
フェーズドアレイの設計シミュレーションって、どういう手順で進めるんですか？
標準的なワークフローは5ステップだ。
- 仕様定義：動作周波数帯、スキャン範囲、サイドローブレベル、利得、偏波を決める
- 単素子設計：パッチ・スロット・ダイポール等の素子を選定し、単体での$S_{11}$・放射パターンを最適化する
- 無限アレイ解析：Floquet境界で素子間隔・基板パラメータを最適化。アクティブ$S_{11}$のスキャン角度依存性を確認
- 有限アレイ解析：DDMまたはフルモデルでエッジ効果・相互結合を検証。放射パターン・利得を確認
- 給電ネットワーク設計：位相シフタ・パワーデバイダ・BFN（Beam Forming Network）を含めた統合解析
最初から全部モデリングしちゃダメなんですか？
ダメではないけど非効率だ。例えば $16 \times 16 = 256$ 素子アレイを最初からフルモデルで解くと、メッシュ数は数千万〜数億要素になる。設計変更のたびに何十時間もかかるような解析では、パラメトリックスタディが回らない。まず単素子で高速にパラメータを最適化してから、無限アレイ → 有限アレイと段階的に複雑さを上げていくのが現実的だ。
Ansys HFSSでのセットアップ
HFSSで無限アレイ解析を組むときのポイントを教えてください。
HFSSのUnit Cell解析のセットアップ手順を具体的に説明しよう。
- ユニットセルの構築：1素子 + 基板 + グランドプレーンを作成。セル境界は素子の中心を基準に $d_x/2$ ずつ
- Floquet Port：アレイ面の上方に配置。Modal数は最低2（TE/TM基本モード）。高次モードの数は $d/\lambda$ に応じて増やす
- Master/Slave Boundary：セルの対向する面にペアで設定。位相差（Phase Delay）はスキャン角度パラメータに連動させる
- 解析設定：Maximum Delta S = 0.01（収束基準）、Maximum Passes = 20、Minimum Converged Passes = 2
- パラメトリックスイープ：$\theta_{\mathrm{scan}}$ を0°〜60°でスイープし、アクティブ $S_{11}$ の変動を確認
スキャン角度ごとに別々のシミュレーションが必要になるんですか？
そうだ。各スキャン角度でFloquet Portの位相条件が変わるから、$\theta_{\mathrm{scan}}$ ごとに1回のFEM解が必要になる。ただしHFSSのOptimetricsを使えば自動化できる。例えば $\theta$ を0°から60°まで5°刻みで13点、$\phi$ を0°と90°の2面で計26ケースを一括投入する。各ケースの適応メッシュは独立に収束させる。
CST Studio Suiteでのセットアップ
CSTの場合はどう違うんですか？
CSTはFDTD法がメインソルバーなので、アプローチが少し異なる。
- Unit Cell Boundary：HFSSのMaster/Slaveに相当。Floquetモードの位相差を指定
- Scan Angle：$\theta$, $\phi$ を直接指定可能。CSTはスキャン角度を内部でFloquet条件に変換する
- Time Domain Solver：広帯域解析に強い。1回の計算でSパラメータの周波数特性が得られる
- Array Factor Post-Processing：単素子の遠方界パターンを元にAFを合成する機能がある
CSTの注意点として、FDTD法は均一格子（Hexahedral Mesh）が基本なので、斜めのエッジやカーブした構造にはステアケースエラーが生じやすい。パッチアンテナのように直交するエッジが多い構造ならFDTDと相性が良いが、円形素子やスパイラルアンテナではFEM（HFSSやCSTのFrequency Domain Solver）の方が精度が出やすい。
V&V — 解析結果の検証
シミュレーション結果が正しいかどうか、どうやって確認するんですか？
V&V（Verification & Validation）は以下のチェックリストで進める。
| チェック項目 | 基準 | NG時の対処 |
|---|---|---|
| $S_{11}$ のメッシュ収束 | Delta S < 0.01（2パス連続） | 最大パス数を増やす / 初期メッシュを細分化 |
| エネルギー保存 | $|S_{11}|^2 + \eta_{\mathrm{rad}} \leq 1$ | PML / 吸収境界の設定見直し |
| ブロードサイド利得 | 理論値 $D = \pi N d^2/\lambda^2 \cdot 4$ と ±1 dB以内 | 相互結合 / 損失の影響を確認 |
| グレーティングローブ | $d \geq \lambda/2$ の場合に出現を確認 | 理論予測角度と実際の位置を照合 |
| スキャン時の利得低下 | $\cos\theta_0$ ファクタと概ね一致 | 素子パターンのcos依存性を確認 |
| 相互結合 | 隣接 $|S_{21}|$ < -10 dB（一般的な目安） | 素子間隔・基板パラメータを見直し |
実験結果との比較も必要ですよね？
もちろん。アンテナの場合は電波暗室（Anechoic Chamber）での放射パターン測定が最終的な検証になる。ニアフィールドスキャナーを使えば、遠方界パターンを近接場の振幅・位相測定からフーリエ変換で再構成できる。シミュレーションと測定の差が2 dB以内であれば、実務上は合格とみなすことが多い。
Starlinkが変えたフェーズドアレイの量産革命
SpaceXのStarlink衛星端末（Dishy McFlatface）は、一般消費者向けに量産されたフェーズドアレイアンテナだ。従来のフェーズドアレイは軍用レーダーや基地局向けで1台数千万〜数億円だったが、Starlinkは約5万円前後で提供している。これが可能になったのは、RFIC（高周波集積回路）の進歩で位相シフタとLNAをワンチップに集積できるようになったからだ。アンテナ設計のシミュレーションも、従来の「1素子に1ヶ月」から「AIアシスト最適化で1週間」に短縮されている。
吸収境界条件の選び方
開放空間をシミュレートするための吸収境界条件は、結果の精度に大きく影響する。HFSS/CSTで使えるRadiation Boundaryは $\lambda/4$ 以上離す必要がある。PML（Perfectly Matched Layer）なら $\lambda/8$ 程度でも十分だが、計算コストは若干増える。「放射面からPMLまでの距離を $\lambda/4$ → $\lambda/2$ に伸ばしたら $S_{11}$ が0.5 dBも変わった」というトラブルは現場で頻出する。PMLの距離は必ず感度解析で確認すべきだ。
ソフトウェア比較
主要ツール比較
フェーズドアレイのシミュレーションに使えるソフトを比較してもらえますか？
商用ツールの主要な比較ポイントをまとめた。
| 機能 | Ansys HFSS | CST Studio Suite | FEKO | COMSOL RF |
|---|---|---|---|---|
| 主要解法 | FEM (適応メッシュ) | FDTD / FEM | MoM / MLFMA | FEM |
| Floquet境界 | Floquet Port | Unit Cell BC | PBC + 無限アレイ | Floquet周期BC |
| DDM（大規模アレイ） | 3D Component DDM | Array Assembly | MLFMA加速 | なし |
| スクリプト自動化 | AEDT + Python API | VBA / Python | LUA / Python | COMSOL API (Java) |
| MPI並列スケーラビリティ | 良好（数百コア） | 良好 | 優秀（数千コア） | 普通 |
| レーダーRCS連携 | SBR+ / Savant | IE + PO | PO / UTD | なし |
| 開発元 | Ansys Inc. | Dassault Systemes SIMULIA | Altair | COMSOL AB |
| ライセンス形態 | ノードロック / フローティング | ノードロック / フローティング | フローティング | ノードロック / フローティング |
どれが一番いいですか？
「一番」は状況による。ざっくり言うと、
- HFSS：業界標準。素子設計〜小中規模アレイに強い。適応メッシュが優秀で初心者にも比較的扱いやすい
- CST：広帯域解析に強い。FDTD + FEM + MoMのマルチソルバーが一つのGUIで使える
- FEKO：MoM/MLFMA ベースで大規模アレイ（数百〜数千素子）の効率計算に強い。プラットフォームレベルのRCS/EMC解析に適する
- COMSOL：マルチフィジクス（電磁+熱+構造）の連成が得意。アンテナ+筐体の発熱解析などに
オープンソースの選択肢
学生なので商用ツールは手が出ないです。オープンソースで使えるものはありますか？
いくつかの選択肢がある。
| ツール | 手法 | 特徴 | 制限 |
|---|---|---|---|
| OpenEMS | FDTD | Octave/MATLAB連携、活発なコミュニティ | 適応メッシュなし、GUIが弱い |
| NEC2 / xnec2c | MoM | ワイヤアンテナに特化、非常に軽量 | 誘電体・マイクロストリップに非対応 |
| Meep (MIT) | FDTD | Python API充実、周期構造対応 | アンテナ特化機能なし |
| PyAEDT | HFSS API | HFSSのPythonラッパー（ライセンス要） | 商用ライセンスが必要 |
学生であればAnsysのStudent版（無償）が最も実用的だ。HFSS単素子設計とFloquet解析が使える。ただしメッシュ数やDDMに制限があるので、大規模アレイの解析は難しい。Pythonでアレイファクタの計算だけなら自分でコーディングするのも良い勉強になるよ。NumPy + Matplotlibで放射パターンの可視化まで30行程度でできる。
選定ガイドライン
結局どうやって選べばいいんですか？ 判断基準を教えてください。
以下の3つの問いに答えれば、ツールは自ずと決まる。
- 「アレイの規模は？」 — 4〜16素子ならほぼ全ソフトでフルモデル解析可能。64素子以上ならDDM対応のHFSSかMLFMA対応のFEKO。数千素子ならFEKO + HPC一択に近い
- 「帯域幅は？」 — 10%以下の狭帯域ならFEM (HFSS)。50%以上の超広帯域（UWBなど）ならFDTD (CST)
- 「連成は必要？」 — 電磁+熱連成（ハイパワーアレイの放熱設計）ならCOMSOL。電磁+構造連成（振動環境のパターン変動）ならHFSS + Ansys Mechanical
HFSSとCSTの「宗教戦争」
アンテナ業界では「HFSSかCSTか」は永遠の論争テーマだ。FEM派は「適応メッシュの精度が最高」と主張し、FDTD派は「1回で広帯域特性が取れるのが合理的」と反論する。実態としては、北米の防衛・宇宙企業はHFSS優勢、欧州の自動車・通信企業はCST優勢という地理的な色分けもある。最近はAnsysもCSTもPython APIを充実させているので、「使い慣れた方を使えばいい」というのが現場の正直な結論だ。
先端技術
5G ミリ波アレイの設計動向
5Gのフェーズドアレイって、今どんな方向に進化してるんですか？
5G NR（New Radio）のミリ波帯（FR2: 24.25〜52.6 GHz）では、端末側にも4〜8素子のフェーズドアレイが搭載されるようになった。3GPPのRelease 17以降では、基地局側は最大1024素子のMassive MIMOが規格化されている。設計の焦点は以下の3つだ。
- AiP（Antenna in Package）：アンテナとRFICをワンパッケージに統合。基板・モールド樹脂の誘電特性がアンテナ性能に直結するため、電磁界シミュレーションでパッケージ全体をモデリングする必要がある
- ビームマネジメント：コードブック設計（離散的なビーム方向のテーブル）のシミュレーション最適化。128ビーム×複数偏波の全パターンをスキャンして最適なコードブックを選定する
- EM-回路協調シミュレーション：アンテナの電磁界シミュレーションとRFIC/ビームフォーマーの回路シミュレーションを統合した協調解析
アンテナデジタルツイン
デジタルツインって最近よく聞きますけど、アンテナでもやるんですか？
やっている。特に衛星通信と防衛分野で活発だ。アンテナデジタルツインの典型的な構成は、
- 物理モデル：高精度なフルウェーブEM解析モデル（HFSS/CSTベース）
- データ駆動モデル：測定データと機械学習で構築した高速代理モデル
- リアルタイム推論：運用中のアンテナから取得するセンサーデータ（温度、振動、受信SNR）をフィードバックしてビーム再較正
例えば、人工衛星の太陽熱変形でアレイ面が0.1 mm反ったときのパターン劣化を即座に予測し、位相オフセットで補正する——こういったリアルタイム運用にデジタルツインが使われ始めている。
機械学習による最適化
AIでアンテナ設計を最適化するのは現実的ですか？
かなり実用化が進んでいる。代表的なアプローチをいくつか紹介しよう。
- サロゲートモデル最適化：数百回のEM解析で訓練データを作り、ガウス過程回帰やニューラルネットワークで代理モデルを構築。多目標最適化（帯域幅 vs. SLL vs. スキャン範囲）を代理モデル上で高速に実行する
- GAN/VAEによるトポロジー生成：生成モデルでアンテナの新しい形状を「発明」させる。人間が思いつかない非直感的な形状で良い特性が出るケースがある
- PINN（Physics-Informed Neural Network）：Maxwell方程式を損失関数に組み込んだニューラルネットで、物理整合性を保ちながら高速にフィールドを予測する
- 強化学習：ビームフォーミングのコードブック設計に適用。環境（電波伝搬チャネル）が動的に変化する場合に有効
AIがアンテナエンジニアの仕事を奪うかもしれないですね。
むしろ逆で、AIを「使える」エンジニアの生産性が飛躍的に上がるんだ。AIは設計空間の探索を加速するけど、物理的な洞察——なぜこの構造が良いのか、実製造で何が問題になるか——は人間の経験が不可欠だ。アレイファクタの理論をしっかり理解した上でAIツールを使いこなすのが、これからのアンテナエンジニアの姿だろう。
トラブルシューティング
収束失敗
HFSSで適応メッシュが全然収束しないんですけど、何が原因ですか？
アンテナ解析で収束しない原因は大きく3つある。
- Radiation Boundaryが近すぎる：放射面から $\lambda/4$ 未満だと近傍界がBCに干渉する。$\lambda/4$ 以上離すか、PMLに変更する
- Floquet Portのモード数不足：$d/\lambda > 0.5$ の場合、高次Floquetモードが伝搬する。Modal数を増やして再計算
- 金属構造の微小ギャップ：0.001 mm以下のギャップがあると、そこに電界が集中してメッシュが爆発する。CADのクリーンアップで微小ギャップを除去するか、意図的にマージする
Delta S = 0.01に収束しない場合、0.02に緩和してもいいですか？
設計探索フェーズなら0.02でも許容できる場合があるが、最終検証では0.01以下を死守すべきだ。収束しない場合は、まず初期メッシュの段階で「Mesh Operations > Length Based」を使って重要領域のメッシュを手動で細分化してから適応リファインを走らせる。これで初期メッシュの品質が上がり、収束が改善されることが多い。
グレーティングローブの見落とし
ブロードサイドでは問題なかったのに、ステアリングしたら急にパターンが崩れたんですけど。
典型的なグレーティングローブの見落としだ。$d = 0.6\lambda$ で設計した場合、ブロードサイドでは $d < \lambda$ なのでグレーティングローブは出ない。しかし $\theta_0 = 45°$ にステアリングすると、$d(1 + \sin 45°) = 0.6 \times 1.707 = 1.02\lambda > \lambda$ となり、グレーティングローブが出現する。対策は以下。
- 素子間隔を $d = 0.5\lambda$ 以下に変更する（根本解決）
- 最大スキャン角度を制限する（仕様レベルの妥協）
- 三角格子配列を採用する（矩形格子より密に詰められる）
相互結合の無視による誤差
アレイファクタの理論計算と、フルウェーブシミュレーションの結果が全然合わないんですが。
それは「理想的なパターン乗算」と「実際の相互結合」の差だ。パターン乗算の原理（$E_{\mathrm{total}} = E_{\mathrm{element}} \times AF$）は、全素子が孤立状態（isolated element pattern）で放射していることを仮定している。しかし実際のアレイでは隣接素子間の電磁結合で各素子のパターンが歪む（embedded element pattern）。
対策としては、
- 理論計算には「アクティブ素子パターン」（相互結合込み）を使う
- フルウェーブ解析と理論値の差が3 dB以上あれば、素子間結合が無視できないレベルと判断する
- アレイの中央素子と端素子でSパラメータを比較し、結合の影響範囲を特定する
スキャンブラインドネス
特定のスキャン角度でいきなり利得がガクッと落ちる現象が出たんですけど、これがスキャンブラインドネスですか？
その通り。スキャンブラインドネスは、アレイ基板の表面波がFloquet高次モードと結合する現象だ。特定の角度で表面波が強く励振されると、素子から放射されるべきエネルギーが基板内に閉じ込められ、$|S_{11}| \to 0$ dB（全反射）に近づく。
スキャンブラインドネスが発生しやすい条件は、
- 基板の誘電率が高い（$\varepsilon_r > 4$）
- 基板が厚い（$h > 0.05\lambda_0$）
- 素子間隔が $\lambda/2$ に近い
対策としては、(1) 低誘電率基板の使用、(2) 基板にビア（接地スタブ）を追加して表面波を抑制、(3) EBG（Electromagnetic Band Gap）構造の挿入がある。シミュレーションでは、アクティブ $S_{11}$ を $\theta$ = 0°〜75° のスキャン角度でスイープし、$|S_{11}|$ が急上昇する角度がないか事前に確認することが必須だ。
フェーズドアレイのシミュレーション、奥が深いですね。理論を理解した上でツールを使わないと、見当外れの結果になりかねないということがよく分かりました。
その感覚は正しい。CAEツールはブラックボックスではなく、物理を理解した上で使う「拡張された計算機」だ。アレイファクタの式を手計算で検証できる力があれば、シミュレーション結果が怪しいときにすぐ気づける。理論とシミュレーションの両輪を回せるエンジニアが、良い設計を生み出すんだ。
なった
詳しく
報告