磁束はどのように計算しますか？

磁束Φは、磁束密度Bと面積要素ベクトルdAの内積を面積全体で積分して計算します（Φ = ∫ B・dA）。面と磁場の方向が角度θずれると、実質的な有効面積はA cosθとなるため、計算時には方向を考慮する必要があります。

磁束の単位ウェーバー(Wb)とは？

1ウェーバー(Wb)は、1平方メートルの面積に1テスラ(T)の一様な磁場が垂直に貫くときの磁束量です。例えば、表面磁束密度0.3Tのフェライト磁石が0.0001m²（1cm²）の面積を垂直に貫く場合、磁束は0.00003Wb（30μWb）となります。

磁束は保存されますか？

はい、磁束保存の法則（∇・B = 0）により、磁束の源（磁気単極子）は存在せず、磁力線は必ず閉じたループを形成します。つまり、ある領域から出る磁束と入る磁束の総和は常にゼロとなり、磁束は保存されます。

磁束 — CAE用語解説

Q: 磁束とは何ですか？

磁束は、ある面積を貫く磁場の総量を表す物理量です。磁束密度Bと面積Aの内積の面積分で定義され、単位はウェーバー(Wb)です。例えば、1平方メートルの面に1テスラの磁場が垂直に貫くと、磁束は1Wbになります。

カテゴリ: 用語集 | 2026-01-15

CAE visualization for magnetic flux - technical simulation diagram

磁束

🧑‍🎓

先生、磁束Φって何ですか？磁束密度Bとの違いがよく分からなくて…

理論と物理

磁束の定義と基本概念

🧑‍🎓

磁束って、教科書には「磁場が貫く量」と書いてありますが、具体的に何を数えているんですか？電荷の数みたいなイメージですか？

🎓

電荷の数とは違いますね。磁束は、ある面積を貫く「磁場の線」の総本数、という比喩が近いです。物理的には、磁束密度Bと面積Aの内積の面積分で定義されます。具体的には、1平方メートルの面に1テスラ(T)の一様な磁場が垂直に貫くとき、磁束は1ウェーバー(Wb)です。例えば、一般的なフェライト磁石の表面では磁束密度が0.3T程度なので、1cm²（0.0001 m²）の面積を貫く磁束は

$$ \Phi = B \cdot A = 0.3 \times 0.0001 = 3 \times 10^{-5} \, \text{Wb} $$

となります。

🧑‍🎓

面積分ということは、面の向きが重要なんですか？コイルを斜めに置くとどうなりますか？

🎓

その通りです。磁束は

$$ \Phi = \int_S \mathbf{B} \cdot d\mathbf{A} $$

で計算されます。ここでdAは面積要素ベクトルで、面の法線方向を向いています。コイル面と磁場の方向が角度θだけずれている場合、実質的に有効な面積はA cosθになります。例えば、先ほどの例でコイルを30度傾けると、磁束は

$$ \Phi = 0.3 \times 0.0001 \times \cos 30^\circ \approx 2.6 \times 10^{-5} \, \text{Wb} $$

と約13%減少します。これがファラデーの電磁誘導則で起電力が生じる磁束変化の元になる量です。

🧑‍🎓

磁束は保存されると聞きましたが、磁石から出た磁束はどこへ行くんですか？消えたりしないんですか？

🎓

良い質問です。マクスウェル方程式の一つ、磁束保存の法則

$$ \nabla \cdot \mathbf{B} = 0 $$

がそれを表しています。これは「磁場の発散はゼロ」、つまり磁束の源（モノポール）は存在せず、磁力線は必ずループを形成することを意味します。永久磁石の場合、N極から出た磁力線は空中や周辺の磁性体を通ってS極に戻り、閉じたループを形成します。CAEで磁界解析をすると、この磁束の連続性を満たすようにソルバーが計算します。磁束が「消える」のは、計算領域の境界条件の設定が不適切な時で、それは物理的ではなく数値的なエラーです。

数値解法と実装

CAEでの磁束計算手法

🧑‍🎓

CAEソフトで磁束を求める時、内部ではどのような計算が行われているんですか？単にBを面積分しているだけですか？

🎓

基本的には面積分ですが、FEM（有限要素法）では主に2つのアプローチがあります。一つは「磁束密度Bを直接積分」する方法。もう一つは、より一般的で精度が高い「磁気ベクトルポテンシャルA」を使う方法です。後者では、要素節点で求めたAの値から、ストークスの定理を用いて磁束を計算します：

$$ \Phi = \oint_C \mathbf{A} \cdot d\mathbf{l} = \int_S (\nabla \times \mathbf{A}) \cdot d\mathbf{A} $$

この方が、特にコイルの巻線全体を通る全磁束（リンク磁束）を求める時に便利で、Ansys MaxwellやJMAGでは標準的にこの方法が採用されています。

🧑‍🎓

「リンク磁束」とは何ですか？普通の磁束とどう違うのですか？

🎓

コイルがN回巻かれている場合、1本の導線が鎖交する磁束を鎖交磁束（リンク磁束）と呼び、実用的にはこれが重要です。1ターンのコイルを貫く磁束をΦとすると、Nターンのコイル全体の鎖交磁束λは

$$ \lambda = N \Phi $$

です。ただし、これは磁束が各ターンに均一に鎖交する理想的な場合です。実際の3次元コイルでは、端のターンと内側のターンで鎖交する磁束が異なるため、CAEでは各導体要素ごとに磁束を計算し、総和を取る必要があります。これが誘起される電圧

$$ V = -d\lambda / dt $$

を決める直接的な量です。

🧑‍🎓

メッシュの粗さは磁束の計算結果にどのくらい影響しますか？特にエアギャップ部分はどうすればいいですか？

🎓

影響は非常に大きいです。磁束密度Bは空間微分（∇×A）で求まるため、メッシュが粗いと勾配を正しく捉えられません。エアギャップは磁気抵抗が集中し、磁束密度が急変する領域なので、特に細かくメッシュ分割する必要があります。経験則では、エアギャップ長さの1/3から1/5のサイズのメッシュを切ります。例えば、ギャップが0.5mmの小型モーターなら、少なくとも0.1mm程度のメッシュサイズをエアギャップ周辺に設定します。Ansysの場合は「Inflation Layer」を設定して、境界層のように細かいメッシュを自動生成させるのが定石です。

実践ガイド

磁束評価のワークフロー

🧑‍🎓

モーター設計でトルクを計算する時、磁束は具体的にどう使われますか？どの断面の磁束を見ればいいですか？

🎓

ブラシレスDCモーターやIPMモーターのトルクTは、大まかに

$$ T \propto \Psi I $$

で表されます。ここでΨは1相あたりの鎖交磁束の実効値、Iは相電流です。CAEでは、ロータ位置を少しずつ回転させながら（例えば電気角で5度ステップ）、各相コイルを貫く鎖交磁束の瞬時値を計算します。その波形からフーリエ展開して基本波成分を取り出し、それがΨになります。見るべき「断面」は、各相のコイル導体全体を囲む任意の面ですが、ソフトウェアでは「Coil Terminal」や「Winding」としてコイルを定義すれば、自動で鎖交磁束を計算してくれます。

🧑‍🎓

磁束の計算結果を検証するには、どういうチェックをすれば信頼性が高まりますか？

🎓

最低限、以下の3点をチェックします。1. **磁束保存の確認**: 鉄心の任意の閉曲面で磁束の正味の出入りを計算し、ゼロに近いか確認。理論上はゼロですが、数値誤差で1e-6 Wb程度の残差は許容範囲です。2. **結果のメッシュ依存性**: メッシュを1.5倍細かくして、注目する磁束の値がどの程度変化するか確認。変化が1%未満なら収束と判断。3. **簡易計算との比較**: 例えば、エアギャップの磁束を

$$ \Phi \approx \frac{NI \cdot \mu_0 A_g}{l_g} $$

のような簡易式（NI:アンペアターン、Ag:ギャップ面積、lg:ギャップ長さ）で見積もり、CAE結果とオーダーが一致するか確認。これで大きな設定ミスを防げます。

ソフトウェア比較

各ソフトの磁束出力機能

🧑‍🎓

Ansys MaxwellとAbaqus（電磁気モジュール）で磁束を取得する操作はどう違いますか？

🎓

根本的にアプローチが異なります。**Ansys Maxwell**は専用の電磁界解析ソルバーで、「Excitations」設定でコイルを定義すると、結果として「Flux Linkage」が直接出力されます。プローブ機能で任意の面を定義し、貫通磁束も計算可能です。一方、**Abaqus/CAE**の電磁気解析は構造解析との連成が主目的で、磁束そのものを直接出力する標準的な結果変数はありません。代わりに、磁気ベクトルポテンシャル（AVOL）や磁束密度（B）を出力し、ユーザーが後処理でField Outputを利用して面積分を定義するか、Pythonスクリプトで計算する必要があります。磁束を頻繁に見る設計なら、MaxwellやJMAGの方が圧倒的に使いやすいです。

🧑‍🎓

無償や低価格のCAEソフト（例えばFEMMやSimscale）でも磁束は正確に計算できますか？

🎓

用途次第です。**FEMM**は2次元静磁界解析に特化しており、コイルの鎖交磁束や任意ラインを通過する磁束の計算機能があります。2D平面内では非常に正確で、多くのモーター初期設計に使われています。ただし3D効果は考慮できません。**Simscale**などのクラウドCAEプラットフォームは、バックエンドにOpenFOAMやCalculiXなどのオープンソースソルバーを使っています。磁界解析（magneto-statics）も可能ですが、専用ソフトに比べると前処理（特に複雑なコイル定義）や磁束に特化した後処理機能は限定的です。試作前の詳細設計では、やはり専用ソフトが有利です。

トラブルシューティング

磁束関連の解析エラー

🧑‍🎓

解析は収束したのに、計算したコイルの磁束が明らかに小さすぎる（または大きすぎる）値になります。どこを疑えばいいですか？

🎓

まず疑うべきは「コイルの定義」と「材料特性」です。1. **導体の登録漏れ**: コイルを構成する全ての導体（例えば、U相の全ターン）が「Winding」グループに正しく登録されているか確認。1本でも漏れると磁束は激減。2. **電流値の単位**: 設定した電流値の単位が[A]か[mA]か。10[A]のつもりが10[mA]だと1000分の1の磁束に。3. **材料のB-Hカーブ**: 特に鉄心材料。デフォルトの「Linear」のままだったり、飽和特性が正しく入力されていないと、透磁率が高すぎて異常に大きな磁束が出ることがあります。軟鉄の飽和磁束密度は約2.0Tであることを念頭に、結果のB分布も確認しましょう。

🧑‍🎓

過渡解析で磁束の時系列波形に不自然なノイズ（ギザギザ）が乗ります。原因と対策は？

🎓

主に2つの原因が考えられます。1. **時間ステップが粗すぎる**: 磁束変化が急峻な部分（コイルの通電ON/OFF時）で、時間ステップが大きいと数値微分の誤差が大きくなりノイズとして現れます。Ansys Maxwellの場合、「Time Step」設定を、少なくとも駆動電流の立ち上がり時間の1/10以下に設定します。2. **メッシュの移動に伴う補間誤差**: ロータが回転するモーター解析では、メッシュが再分割される「Remeshing」のタイミングで、変数値の補間誤差が生じます。これを軽減するには、「Band」領域のメッシュを細かくする、または「Sliding Mesh」手法が使えるソルバーではそれを選択します。また、結果を後処理で移動平均フィルタ（例えば1電気角周期で平均）をかけて実用的な値を抽出する方法もあります。

🧑‍🎓

「磁束が計算領域の境界を通過しています」という警告が出ました。これは問題ですか？

🎓

重大な問題を示す警告です。磁束保存則により、磁力線は閉じていなければならず、計算領域の境界で突然途切れることは物理的にあり得ません。この警告が出る原因は、**境界条件の設定が不適切**な場合がほとんどです。例えば、無限遠を模倣する「Balloon Boundary」や「Open Boundary」を設定すべき所に、誤って「磁気壁（Neumann条件）」や「対称境界」を設定していると発生します。対策は、モデルの対称性を正しく理解し、磁力線が自然に流出・流入できるように境界条件を見直すことです。特に外部磁場がない場合、モデル全体を囲む領域に「真空」や「空気」の材質を割り当て、その外側境界を「Balloon」に設定するのが安全です。