ローレンツ力連成解析 — 電磁力と構造変形のカップリング手法

いやいや、構造屋こそ知っておくべきテーマだよ。例えばMRIの超伝導磁石——起磁時にコイルに流れる数百アンペアの電流が強磁場の中でローレンツ力を受けて、巻線が膨張しようとする。この力が数トンにもなるから、支持構造が耐えられるか構造解析で評価しないといけない。

もっと身近な例だと、変圧器の短絡事故。短絡電流が通常の10〜25倍に跳ね上がると、コイル間に数百kNの電磁力が瞬間的に加わる。これが巻線を座屈させたり、絶縁を破壊したりする。電力会社が「変圧器の短絡耐力試験」を重視するのはこのためだ。他にもリニアモーター、電磁成形、レールガン——電流と磁場がある場所には必ずローレンツ力がある。

その通り、出発点は荷電粒子に働く力だ。電場と磁場の両方を含む一般形はこう:

CAEで扱う導体の連成解析では、これを連続体に拡張した体積力密度の形で使う。電流密度 $\mathbf{J}$（A/m²）と磁束密度 $\mathbf{B}$（T）の外積だ:

そう、重力と同じように体積力としてFEMの右辺に入る。構造の運動方程式で書くと:

ここで $\boldsymbol{\sigma}$ は応力テンソル、$\mathbf{f}_{\text{em}} = \mathbf{J} \times \mathbf{B}$ が電磁体積力、$\rho$ は密度、$\ddot{\mathbf{u}}$ は加速度。静解析なら右辺は0になる。ざっくり言うと「電磁界解析で $\mathbf{J}$ と $\mathbf{B}$ を求めて外積を取り、その力を構造解析に食わせる」——これがローレンツ力連成の骨格だ。

いい質問だ。J × B は電流が流れている領域にしか力が定義できない。でも、永久磁石やフェライトコアのような磁性体には電流が流れていなくても磁場による力が働く。このとき使うのがMaxwell応力テンソルだ:

この応力テンソルを物体表面で面積分すると、物体に作用する電磁力の合力が得られる:

理論的には物体を囲む任意の閉曲面で同じ結果になるはずだけど、数値的にはメッシュの粗さで変わってしまう。実務では積分面を物体表面から少し離す「エアギャップ積分」が安定するね。COMSOLではデフォルトでこの手法が使われている。

鋭いね。磁歪（magnetostriction）はローレンツ力とは別の電磁-構造結合メカニズムだ。磁性体が磁化されると、磁区の回転に伴って結晶格子が微小に変形する。変圧器の鉄心は50Hz/60Hzの交流磁場で繰り返し伸縮するから、あの唸り音が出る。

磁歪は歪みテンソル $\varepsilon^{\text{ms}}$ として構造方程式に入る:

熱膨張の式 $\sigma = C(\varepsilon - \alpha \Delta T)$ と構造が同じだろう？ 磁歪歪みは典型的に $10^{-6}$〜$10^{-5}$ のオーダーで、電磁鋼板で約 $\lambda_s \approx 5 \times 10^{-6}$。微小に見えるけど、鉄心全体で足し合わさると変圧器の騒音や振動の主原因になる。

その通り。実際の問題ではこれらが重畳することも多い。例えば電動モーターでは、ステーター巻線にローレンツ力が働きつつ、鉄心では磁歪振動が起きる。騒音解析をやるなら両方とも考慮しないと実験と合わない。

多くの場合は弱連成で十分だよ。MRI磁石の静的評価とか、変圧器の短絡耐力解析とかは、構造変形が微小（mmオーダー以下）だから電磁場への影響は無視できる。

ただし強連成が必要なケースもある。代表例は電磁成形。アルミ板に渦電流を瞬間的に流して磁場との相互作用で高速変形させる加工法だけど、変形でコイルと被加工材の距離が変わると磁場分布も電流分布も大きく変わる。片方向だと力を過大評価してしまう。リニアモーターのエアギャップ変動問題も同様だね。

低周波の電磁界解析では、Maxwellの方程式から変位電流を無視した準静的近似を使う。磁気ベクトルポテンシャル $\mathbf{A}$ とスカラポテンシャル $\phi$ で定式化するA-φ法が標準的だ:

$\mu$ は透磁率、$\sigma$ は導電率、$\mathbf{J}_s$ は外部印加電流密度。ここから磁束密度 $\mathbf{B} = \nabla \times \mathbf{A}$ を求めて、導体中の誘導電流密度は $\mathbf{J} = -\sigma(\partial \mathbf{A}/\partial t + \nabla \phi)$ で得られる。これらの外積がローレンツ力だ。

まさにその通り。ただし重要な違いがある。電磁場のFEMではエッジ要素（Nedelec要素）を使うことが多い。通常の節点要素だと磁束密度の接線連続性が保証されないんだ。構造のラグランジュ要素とは別物だから、連成するときに要素タイプが異なることに注意が必要だよ。

主に3つの手法がある。それぞれ得意・不得意があるから、用途によって使い分ける:

概念は同じだよ。物体を仮想的に微小変位 $\delta s$ させたときの磁場エネルギー変化から力を求める:

JMAGやAnsys Maxwellのようなツールでは、この仮想仕事法をノーダルフォースとして自動計算してくれる。精度と安定性のバランスが良いから、トルク計算などではデファクトスタンダードだ。

ここが連成解析の肝だね。メッシュマッピングには3つの方法がある:

電磁側で計算した合力と、構造側に転送した後の合力が一致すること。これが崩れると、構造モデルに存在しない力が勝手に発生したり消えたりする。変圧器の短絡力解析で力の保存性が5%以上ずれると、応力評価が定量的に使い物にならなくなるから、転送後に必ず合力を確認する習慣をつけよう。

これが過渡連成の厄介なところだ。電磁場は表皮効果の時定数（ミリ秒〜マイクロ秒）で変化するのに、構造の固有振動周期はミリ秒〜秒のオーダー。時間スケールが2〜3桁違うことがざらにある。

対策としてはサブサイクリングが有効だ。電磁場は細かい時間刻み $\Delta t_{\text{em}}$ で解き、構造は粗い時間刻み $\Delta t_{\text{str}} = n \cdot \Delta t_{\text{em}}$ で解く。力の転送は構造の各ステップで、その区間の電磁力を時間平均して渡す。変圧器の短絡解析だと、電磁場を $\Delta t = 0.1 \text{ms}$、構造を $\Delta t = 1 \text{ms}$ くらいで回すのが典型的だね。

弱連成（片方向）の場合の標準ワークフローを説明するよ。多くの実務ケースはこれで対応できる:

よくあるパターンだよ。例えばJMAGで電磁場を解いて、力をNastranやAbaqusに渡すケース。JMAGにはNastran用の力データ出力機能がある。Ansys陣営なら、Ansys MaxwellとMechanicalをWorkbenchで繋ぐのが最も楽だ。異ソルバー間連携の場合は、力のファイルフォーマット（CSV、UNV、.bdf等）とメッシュマッピング方式に注意が必要だね。

電磁メッシュの最重要ポイントは表皮深さに対するメッシュ解像度だ。表皮深さ $\delta$ は:

導体表面から深さ $\delta$ 以内に電流の63%が集中するから、この領域に最低3層、できれば5層の要素を入れる必要がある。例えば銅の50Hz時の表皮深さは約9.4mmだけど、短絡事故の過渡初期は数kHzの高調波成分があるから $\delta < 1$ mm にもなる。

そう、構造屋がつまずきやすいポイントだね。電磁場は空気中にも存在するから、導体や磁性体だけでなく周囲空間もメッシュが必要。無限遠を有限領域で打ち切るから、境界を近づけすぎると磁束が歪んで力の計算精度が落ちる。対策として無限要素（Infinite Element）やPML（Perfectly Matched Layer）を使うこともある。

代表的なものを挙げよう:

ITERでは、プラズマが突然崩壊したとき（ディスラプション）にブランケットに数MNの電磁力が数ミリ秒で作用する。構造の過渡応答を評価しないと安全設計ができない。Abaqusの電磁-構造連成モジュールやCOMSOLが使われているよ。

主要なツールを比較してみよう。ポイントは「電磁場解析の質」と「構造への連成パスの成熟度」の両方だ:

その通り。COMSOLはマルチフィジクスを1つのGUI・1つのメッシュ・1つのソルバーで扱えるのが最大の強み。メッシュマッピングの問題が発生しないし、強連成もシームレスにできる。ただしソルバーの速度はAnsysやJMAGに比べると大規模問題で劣る傾向がある。

JMAGは電磁力を.bdfファイル（Nastran形式）や.csvで出力できるから、それをNastranやAbaqusに読み込ませる。JMAG-RTとの連携でMBDシミュレーション（Adams等）と繋げて、モーターの電磁加振力→ハウジング振動→騒音というフルチェーンの連成もできる。日本の自動車メーカーのモーター開発ではこの手法が主流だね。

注目すべきトレンドをいくつか紹介しよう:

教師あり学習でFEMのデータを使うアプローチもあるし、PINNsならデータなしで支配方程式の残差だけで学習できる。ただし現時点では精度がFEMに及ばないケースが多く、「設計初期段階のスクリーニング」や「実時間モニタリングのサロゲート」に使う研究が多いね。

まさに最前線のテーマだ。電力変圧器のデジタルツインでは、運転中の負荷電流を入力として電磁力→巻線応力をリアルタイムに推定し、累積疲労を監視する構想がある。フルFEMでは計算が間に合わないから、モード展開法やROM（Reduced Order Model）で電磁力→構造応答のマッピングを事前に学習しておく。クエンチ（超伝導消失）の検知と構造健全性評価を数秒以内に行うMRI用のデジタルツインも研究されている。

連成特有のトラブルは大きく3つに分類できる。順番に見ていこう。

あるある。電磁力の空間振動（ノイズ）は最も頻繁に遭遇するトラブルだ。原因はいくつかある:

対策：まず電磁メッシュを表皮深さの1/5以下に細分化してみる。それでもノイズが残る場合は、力の空間フィルタリング（ガウシアンスムージング）を適用する。COMSOLでは「Volume Average」オペレータ、Maxwellでは力のスムージングオプションが使える。ただし過度なスムージングは物理を歪めるから、合力の値が変わらないことを確認すること。

緩和は有効だけど、まず原因を特定しよう。双方向連成の発散パターンは主に2つ:

Aitken緩和の更新式はこうだ:

$\mathbf{r}_k$ は第 $k$ 反復の残差ベクトル。COMSOLのSegregated Solverにはこれが組み込まれている。Ansys Workbenchでは手動で緩和係数を設定する必要があるから、最初は保守的に $\omega = 0.3$ から始めて徐々に上げるのが安全だ。

よく遭遇するものをまとめるよ:

典型的なやつだね。連成解析では「手を動かす前に紙の上で力のオーダーエスティメーションをやる」ことが何より大事だ。ざっくり $F = J \times B \times V$（$V$ は導体体積）で合力の桁数を見積もってから解析結果と突き合わせる。桁が1つ以上違ったら何かがおかしい。