ハルバッハ配列 — 理論と支配方程式

カテゴリ: 電磁場解析 | 2026-01-15
halbach-array-theory
理論と物理の世界へ

ハルバッハ配列とは

🧑‍🎓

先生、ハルバッハ配列って普通の磁石配列と何が違うんですか?


🎓

磁化方向を段階的に回転させることで、片面に磁束を集中させる配列だ。理想的なハルバッハ配列では片面の磁束密度が2倍、反対面はゼロになる。


$$ B_{one\_side} = B_r \left(1 - e^{-2\pi d/\lambda}\right) \cdot \frac{\sin(\pi/M)}{\pi/M} $$

$d$: 磁石厚、$\lambda$: 磁化パターンの波長、$M$: 1波長あたりの磁石セグメント数。


🧑‍🎓

M=4が一般的(90度ずつ回転)ですね。


🎓

そう。M=4で$\sin(\pi/4)/(\pi/4) \approx 0.90$。Mを増やすほど理想に近づくが、組立精度と部品点数が増える。リニアモータ、磁気浮上、Wiggler磁石に使われる。


まとめ

🎓
  • 片面集中 — 一方の磁束密度を最大2倍に増強
  • $M$セグメント — M=4(90度回転)が標準構成
  • 用途 — リニアモータ、SPM型モータ、粒子加速器

Coffee Break よもやま話

電気自動車モータ開発と電磁界解析

テスラのModel 3のモータは、リラクタンストルクと磁石トルクの両方を使うIPMSM(埋込磁石型同期モータ)。この複雑な磁場分布を最適化するには数千回の電磁界FEA解析が必要です。1回の解析に数分としても、最適化ループ全体では数週間のCPU時間。それでも実機を何十台も試作するよりは圧倒的に速くて安い。

各項の物理的意味
  • 電場項 $\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仮定(断面解析):電流方向が一様で、端部効果を無視できる場合に有効
  • 等方性仮定:異方性材料(珪素鋼板の圧延方向等)では方向別の特性定義が必要
  • 適用外ケース:プラズマ(電離気体)、超伝導体、非線形光学材料では追加の構成則が必要
次元解析と単位系
変数SI単位注意点・換算メモ
磁束密度 $B$T(テスラ)1T = 1 Wb/m²。永久磁石: 0.2〜1.4T
磁場強度 $H$A/mB-Hカーブの横軸。CGS系のOe(エルステッド)との換算: 1 Oe = 79.577 A/m
電流密度 $J$A/m²導体の断面積と総電流から算出。表皮効果で不均一分布に注意
透磁率 $\mu$H/m$\mu = \mu_0 \mu_r$。真空中 $\mu_0 = 4\pi \times 10^{-7}$ H/m
導電率 $\sigma$S/m銅: 約5.96×10⁷ S/m。温度上昇で低下

数値例:表皮深さの周波数依存性(銅, σ=5.96×10⁷ S/m, μr=1)

表皮深さ δ = √(2/(ωμσ))。周波数が上がるほど電流は表面に集中します:

50 Hz(家庭用電源)δ = 9.22 mm
1 kHz(オーディオ帯域)δ = 2.06 mm
100 kHz(IH調理器)δ = 0.206 mm
1 GHz(携帯電話)δ = 0.002 mm

GHz帯では電流はほぼ表面のみ! 高周波回路で導体の表面粗さが性能に直結するのはこのためです。メッシュもδの1/3以下の要素サイズが必要です。

簡易計算ツール:電磁気学基礎

表皮深さ δ = √(2/(ωμσ)) を計算します。高周波電磁場の浸透深さの目安。

電磁界解析の精度と計算コストの両立は永遠の課題です。 — Project NovaSolverは、既存ワークフローの改善を目指す取り組みとして、この問題に向き合っています。

Project NovaSolver — CAE実務の課題に向き合う研究開発

「ハルバッハ配列をもっと効率的に解析できないか?」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。

進捗通知を受け取る →