特許 7121111 磁場解析装置、解析方法、及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-08-08

(45)【発行日】2022-08-17

(54)【発明の名称】磁場解析装置、解析方法、及びプログラム

(51)【国際特許分類】

   G06F 30/20 20200101AFI20220809BHJP

   G16Z 99/00 20190101ALI20220809BHJP

【ＦＩ】

G06F30/20

G16Z99/00

【請求項の数】 12

(21)【出願番号】P 2020504776

(86)(22)【出願日】2018-11-12

(86)【国際出願番号】 JP2018041831

(87)【国際公開番号】W WO2019171660

(87)【国際公開日】2019-09-12

【審査請求日】2021-05-19

(31)【優先権主張番号】P 2018040400

(32)【優先日】2018-03-07

(33)【優先権主張国・地域又は機関】JP

(73)【特許権者】

【識別番号】000002107

【氏名又は名称】住友重機械工業株式会社

(74)【代理人】

【識別番号】100105887

【弁理士】

【氏名又は名称】来山 幹雄

(72)【発明者】

【氏名】宮崎 修司

【審査官】岡本 俊威

(56)【参考文献】

【文献】特開２００３－００６１８１（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｆ ３０／２０

Ｇ１６Ｚ ９９／００

(57)【特許請求の範囲】

【請求項1】

仮想空間を複数の体積要素に分割し、前記体積要素の各々に磁化を付与する磁化付与部と、
前記仮想空間の複数の観測点の各々について、前記観測点の周りの複数の前記体積要素に付与された磁化に基づいて、前記体積要素ごとに前記観測点に発生する磁場を演算し、複数の前記体積要素ごとの演算結果に基づいて、複数の前記観測点の各々に発生する磁場を求める磁場演算部と
を有し、
前記磁化付与部は、
前記仮想空間を、形状及び寸法が同一の複数のボクセルに分割し、
複数の前記ボクセルの各々に磁化を付与し、
第１方向に連続する複数の前記ボクセルをまとめて１つの前記体積要素とする粗視化を行うことにより、少なくとも一部の前記体積要素について、前記第１方向に関する寸法を、他の第２方向に関する寸法より大きくし、
粗視化された前記体積要素に含まれる複数の前記ボクセルに付与されている磁化に基づいて、粗視化された前記体積要素に磁化を付与する磁場解析装置。

【請求項2】

前記体積要素に付与される磁化と、周りの複数の前記観測点に発生する磁場との対応関係を、前記体積要素と前記観測点との相対位置ごとに保持する対応関係保持部を、さらに有し、
前記磁場演算部は、前記対応関係を用いて、前記体積要素に付与されている磁化が前記観測点に発生させる磁場を求める第１機能を持つ請求項１に記載の磁場解析装置。

【請求項3】

前記対応関係保持部は、前記体積要素の大きさごとに前記対応関係を保持しており、
前記磁場演算部は、前記体積要素に付与されている磁化が前記観測点に発生させる磁場を計算する際に、磁場を発生させるソースとなる前記体積要素の大きさに対応する前記対応関係を用いて磁場を計算する請求項２に記載の磁場解析装置。

【請求項4】

前記磁場演算部は、複数の前記ボクセルに付与されている磁化をそれぞれ磁気双極子で代表して、前記磁気双極子に基づいて前記観測点の磁場を求める第２機能を有する請求項２または３に記載の磁場解析装置。

【請求項5】

前記第２機能は、複数の前記ボクセルの前記磁気双極子によってそれぞれ前記観測点に発生する磁場を前記ボクセルごとに求め、前記ボクセルごとに求められた磁場を足し合わせる機能を含む請求項４に記載の磁場解析装置。

【請求項6】

前記第２機能は、一方向に一列に並ぶ複数の前記ボクセルの各々の前記磁気双極子を、１つの多重極展開点に多重極展開し、前記多重極展開点に対する前記観測点の相対位置に基づいて前記観測点の磁場を求める機能を含む請求項４に記載の磁場解析装置。

【請求項7】

１つの前記多重極展開点に多重極展開される複数の前記磁気双極子の大きさ及び向きは同一であり、１つの前記多重極展開点に多重極展開される複数の前記磁気双極子に対応する前記ボクセルの寸法は同一である請求項６に記載の磁場解析装置。

【請求項8】

前記磁場演算部は、磁場を発生させるソースとなる前記体積要素から前記観測点までの距離が第１距離以下の少なくとも一部の前記観測点に発生する磁場を、前記第１機能を用いて求め、前記第１距離より長い前記観測点に発生する磁場を、前記第２機能を用いて求める請求項４乃至７のいずれか１項に記載の磁場解析装置。

【請求項9】

前記磁場演算部は、磁場を発生させるソースとなる前記体積要素に付与された磁化を含む物理量を、前記体積要素の表面について面積分することにより、前記体積要素の周りの前記観測点に発生する磁場を求める第３機能を有する請求項８に記載の磁場解析装置。

【請求項10】

前記磁場演算部は、磁場を発生させるソースとなる前記体積要素から前記観測点までの距離が前記第１距離より短い第２距離以下の前記観測点に発生する磁場を、前記第３機能を用いて求める請求項９に記載の磁場解析装置。

【請求項11】

仮想空間を複数の体積要素に分割し、前記体積要素の各々に磁化を付与し、
前記体積要素に付与された磁化が前記仮想空間内の複数の観測点に発生させる磁場を前記体積要素ごとに計算し、複数の前記体積要素ごとの計算結果に基づいて、複数の前記観測点の各々に発生する磁場を求める磁場解析方法であって、
前記仮想空間を前記複数の体積要素に分割し、前記体積要素の各々に磁化を付与するにあたり、
前記仮想空間を、形状及び寸法が同一の複数のボクセルに分割し、
複数の前記ボクセルの各々に磁化を付与し、
第１方向に連続する複数の前記ボクセルをまとめて１つの前記体積要素とする粗視化を行うことにより、少なくとも一部の前記体積要素について、前記第１方向に関する寸法を、他の第２方向に関する寸法より大きくし、
粗視化された前記体積要素に含まれる複数の前記ボクセルに付与されている磁化に基づいて、粗視化された前記体積要素に磁化を付与する磁場解析方法。

【請求項12】

仮想空間を複数の体積要素に分割し、前記体積要素の各々に磁化を付与する機能と、
前記体積要素に付与された磁化が前記仮想空間内の複数の観測点に発生させる磁場を前記体積要素ごとに計算し、複数の前記体積要素ごとの計算結果に基づいて、複数の前記観測点の各々に発生する磁場を求める機能と、
前記仮想空間を前記体積要素に分割するにあたり、少なくとも一部の前記体積要素について、前記仮想空間内の第１方向に関する寸法を、他の第２方向に関する寸法より大きくする機能と
をコンピュータに実現させるためのプログラムであって、
前記仮想空間を前記複数の体積要素に分割し、前記体積要素の各々に磁化を付与する機能は、
前記仮想空間を、形状及び寸法が同一の複数のボクセルに分割し、
複数の前記ボクセルの各々に磁化を付与し、
第１方向に連続する複数の前記ボクセルをまとめて１つの前記体積要素とする粗視化を行うことにより、少なくとも一部の前記体積要素について、前記第１方向に関する寸法を、他の第２方向に関する寸法より大きくし、
粗視化された前記体積要素に含まれる複数の前記ボクセルに付与されている磁化に基づいて、粗視化された前記体積要素に磁化を付与する機能を含むプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、磁場解析装置、解析方法、及びプログラムに関する。

【背景技術】

【0002】

近年、コンピュータの処理能力向上に伴い、モータ等の電気機器の設計開発の現場で、磁場解析を取り入れたシミュレーションが行われる。シミュレーションを行うと、実際にプロトタイプを作製しなくてもある程度の評価が可能となるので、設計開発のスピードが向上し得る。

【0003】

古典力学や量子力学等を基にコンピュータを用いて物質科学全般の現象を探るための方法として、分子動力学法を、マクロスケールの系を扱えるように発展させた繰り込み群分子動力学法に基づくシミュレーションが知られている（例えば、特許文献１）。分子動力学法のような粒子法は、静的な現象だけではなく流れ等の動的な現象も取り扱えるため、主として静的な現象を解析対象とする有限要素法等に代わるシミュレーション手法として注目されている。

【0004】

例えば、各粒子に磁気モーメントを付与し、各粒子の球対称性に基づく厳密解を利用して磁気的な物理量を演算対象とすることで、比較的精度の高いシミュレーション結果を高速に得ることができる磁気ビーズ法が提案されている（例えば、特許文献２）。磁気ビーズ法では、対象物を複数の要素に分割し、各要素を球形状の粒子とみなして球対称性に基づく厳密解を利用して演算を行う。このとき、対象物を複数の粒子の集合体で表すと、粒子間に隙間が生じるため、その隙間や球形状の粒子に起因する誤差が生じ得る。計算精度を高めるためには、対象物を立方体等のボクセルで隙間なく分割し、ボクセルの形状に応じた厳密解に基づく面積分を実行する必要がある。この面積分を高精度で行おうとすると計算時間が膨大になってしまう。

【0005】

磁場演算の計算時間の短縮化を図ることが可能な手法が提案されている（例えば、特許文献３）。この手法では、磁気モーメントが付与される基準粒子が発生させる磁場の数値解を、基準粒子の周りに設定される複数の格子点ごとに予め算出し、基準粒子に付与される磁気モーメントのベクトル量と、磁場の数値解との対応関係を格子点ごとに保持しておく。磁気モーメントが付与される各粒子が発生させる磁場を、この対応関係を用いて演算する。各粒子が発生させる磁場を計算する度に面積分の計算を行う必要がないため、計算時間を短縮することができる。

【先行技術文献】

【特許文献】

【0006】

【文献】特開２００９－３７３３４号公報

【文献】特開２０１５－１１１４０１号公報

【文献】特開２０１７－１９４８８４号公報

【発明の概要】

【発明が解決しようとする課題】

【0007】

計算時間をより短縮する手法が望まれている。本発明の目的は、計算時間をより短縮することが可能な磁場解析装置、解析方法、及びプログラムを提供することである。

【課題を解決するための手段】

【0008】

本発明の一観点によると、
仮想空間を複数の体積要素に分割し、前記体積要素の各々に磁化を付与する磁化付与部と、
前記仮想空間の複数の観測点の各々について、前記観測点の周りの複数の前記体積要素に付与された磁化に基づいて、前記体積要素ごとに前記観測点に発生する磁場を演算し、複数の前記体積要素ごとの演算結果に基づいて、複数の前記観測点の各々に発生する磁場を求める磁場演算部と
を有し、
前記磁化付与部は、
前記仮想空間を、形状及び寸法が同一の複数のボクセルに分割し、
複数の前記ボクセルの各々に磁化を付与し、
第１方向に連続する複数の前記ボクセルをまとめて１つの前記体積要素とする粗視化を行うことにより、少なくとも一部の前記体積要素について、前記第１方向に関する寸法を、他の第２方向に関する寸法より大きくし、
粗視化された前記体積要素に含まれる複数の前記ボクセルに付与されている磁化に基づいて、粗視化された前記体積要素に磁化を付与する磁場解析装置が提供される。

【0009】

本発明の他の観点によると、
仮想空間を複数の体積要素に分割し、前記体積要素の各々に磁化を付与し、
前記体積要素に付与された磁化が前記仮想空間内の複数の観測点に発生させる磁場を前記体積要素ごとに計算し、複数の前記体積要素ごとの計算結果に基づいて、複数の前記観測点の各々に発生する磁場を求める磁場解析方法であって、
前記仮想空間を前記複数の体積要素に分割し、前記体積要素の各々に磁化を付与するにあたり、
前記仮想空間を、形状及び寸法が同一の複数のボクセルに分割し、
複数の前記ボクセルの各々に磁化を付与し、
第１方向に連続する複数の前記ボクセルをまとめて１つの前記体積要素とする粗視化を行うことにより、少なくとも一部の前記体積要素について、前記第１方向に関する寸法を、他の第２方向に関する寸法より大きくし、
粗視化された前記体積要素に含まれる複数の前記ボクセルに付与されている磁化に基づいて、粗視化された前記体積要素に磁化を付与する磁場解析方法が提供される。

【0010】

本発明のさらに他の観点によると、
仮想空間を複数の体積要素に分割し、前記体積要素の各々に磁化を付与する機能と、
前記体積要素に付与された磁化が前記仮想空間内の複数の観測点に発生させる磁場を前記体積要素ごとに計算し、複数の前記体積要素ごとの計算結果に基づいて、複数の前記観測点の各々に発生する磁場を求める機能と、
前記仮想空間を前記体積要素に分割するにあたり、少なくとも一部の前記体積要素について、前記仮想空間内の第１方向に関する寸法を、他の第２方向に関する寸法より大きくする機能と
をコンピュータに実現させるためのプログラムであって、
前記仮想空間を前記複数の体積要素に分割し、前記体積要素の各々に磁化を付与する機能は、
前記仮想空間を、形状及び寸法が同一の複数のボクセルに分割し、
複数の前記ボクセルの各々に磁化を付与し、
第１方向に連続する複数の前記ボクセルをまとめて１つの前記体積要素とする粗視化を行うことにより、少なくとも一部の前記体積要素について、前記第１方向に関する寸法を、他の第２方向に関する寸法より大きくし、
粗視化された前記体積要素に含まれる複数の前記ボクセルに付与されている磁化に基づいて、粗視化された前記体積要素に磁化を付与する機能を含むプログラムが提供される。

【発明の効果】

【0011】

体積要素の第１方向に関する寸法を第２方向に関する寸法より大きくすることで、体積要素の第２方向の寸法を一辺の長さとする正六面体のボクセルの個数より、体積要素の個数の方が少なくなる。このため、ボクセルごとに演算を行う場合に比べて、体積要素ごとに演算を行う方が、演算負荷が低くなる。また、磁場の空間的な変化が緩やかな方向を第１方向に選ぶと、計算精度の低下を抑制することができる。

【図面の簡単な説明】

【0012】

【図1】図１は、本実施例による磁場解析装置のブロック図である。

【図2】図２Ａは、対象物の一部分を複数のボクセルに分割した状態の斜視図であり、図２Ｂは、粗視化後の対象物の一部分の斜視図である。

【図3】図３Ａは、仮想空間内に定義される座標系を示す図であり、図３Ｂは、１つの体積要素の内部の磁化が一定であるとしたときに、体積要素の磁化によって観測点に発生する磁場Ｈ（ｒ_ｉ）を計算する方法を説明する図である。

【図4】図４は、ソース点となる体積要素と、その周りの複数の基準点との位置関係を示す図である。

【図5】図５は、磁化が付与されたボクセルと観測点との位置関係を示す図である。

【図6】図６は、ソース点となる体積要素と観測点Ｐｏとの位置関係を示す図である。

【図7】図７Ａ及び図７Ｂは、ソース点と観測点との位置関係を示す図である。

【図8】図８Ａ及び図８Ｂは、それぞれ解析対象の磁石電動機に用いられる永久磁石の平面図及び側面図であり、図８Ｃは、実施例による方法で磁場解析を行った結果、及び比較例による方法で磁場解析を行った結果を示すグラフである。

【図9】図９は、多重極展開の方法を説明するために用いる座標系を示す図である。

【図10】図１０Ａは、個々のボクセルについて双極子近似を行い、観測点Ｐｏに発生する磁場Ｈ（ｒ_ｉ）を計算により求める方法を説明する図であり、図１０Ｂは、多重極展開手法を用い、観測点Ｐｏに発生する磁場Ｈ（ｒ_ｉ）を計算により求める方法を説明する図である。

【図11】図１１は、磁場Ｈｏに対する磁場Ｈｃの誤差を、規格化距離の関数として示すグラフである。

【図12】図１２は、多重極展開を行うことによる計算の高速化効果を示すグラフである。

【発明を実施するための形態】

【0013】

古典力学や量子力学等を基に計算機を用いて物質科学全般の現象を探るための方法として、分子動力学法（Molecular Dynamics Method、以下ＭＤ法と称す。）や、量子分子動力学法（Quantum Molecular Dynamics Method）や、マクロスケールの系を扱えるようにＭＤ法を発展させた繰り込み群分子動力学法（Renormalized Molecular Dynamics、以下ＲＭＤ法と称す。）に基づくシミュレーションが知られている。正確にはＭＤ法やＲＭＤ法は運動論的手法（物理量の算出には統計力学を使う）であり、粒子法は、連続体を記述する微分方程式を離散化する手法であり、両者は別ものであるが、ここではＭＤ法やＲＭＤ法も粒子法と呼ぶ。

【0014】

粒子法は静的な現象だけでなく流れなどの動的な現象をも取り扱えるので、主に静的な現象を解析対象とする有限要素法などに代わるシミュレーション手法として注目されている。

【0015】

粒子法には、連続体を粒子で離散化することにより解析対象の粒子系を得るという微分的な見方がある。例えば、粒子法において流体を扱う場合、ナビエストークス（Navier-Stokes）方程式を粒子で離散化することが多い。

【0016】

一方、粒子法の別の見方として、多くの粒子を集めて連続体を形成するという積分的な見方もある。これは例えば、小さな鉄の粒を集めて固めて大きな鉄球を形成するという見方である。

【0017】

一般に、多くの磁気モーメントが存在する空間内のある点（観測点）の磁場を求める場合、重ね合わせの原理により、各磁気モーメントが観測点に発生させる磁場を磁気モーメントに亘って足し合わせる。本発明者は、このような磁気モーメントの集まりから磁場を求める手法と、積分的な見方をした場合の粒子法との親和性を独自に見い出し、粒子法における粒子に磁気モーメントを付与することに想到した。これにより、対流や大変形の解析に強いという粒子法の利点を維持しつつ、粒子法の適用範囲を磁場解析にまで広げることが可能となる。

【0018】

対象物を球形の粒子の集合体で表し、各粒子の球対称性に基づく厳密解を利用することで、高速かつ高精度の演算結果を得ることができる。しかしながら、連続体である対象物を球形の粒子の集合体で表すと粒子間に隙間が生じるため、その隙間に起因する計算誤差が生じてしまう。特に、粒子間の距離が近い場合の粒子近傍の磁場を計算する場合に誤差の影響が大きくなる。計算精度を高めるためには、対象物を立方体やボロノイ多面体といったボクセルに分割し、これらのボクセルにより隙間なく対象物を分割することが好ましい。この場合、磁場を求めるために、ボクセルごとに面積分を実行する必要がある。この面積分を実行するには、計算時間が膨大となる。

【0019】

以下に説明する実施例による磁場解析装置では、計算時間の短縮化を図ることが可能である。
図１は、本実施例による磁場解析装置３０のブロック図である。本実施例による磁場解析装置３０は、処理装置２０、記憶装置２５、入力装置２８、及び出力装置２９を含む。処理装置２０は、解析情報取得部２１、磁化付与部２２、磁場演算部２３、及び出力制御部２４を含む。記憶装置２５は、データ保持部２６及び対応関係保持部２７を含む。

【0020】

図１に示す各ブロックは、ハードウェア的には、コンピュータの中央処理ユニット（ＣＰＵ）をはじめとする素子や機械装置で実現することができ、ソフトウェア的にはコンピュータプログラム等によって実現することができる。図１では、ハードウェア及びソフトウェアの連携によって実現される機能ブロックが示されている。従って、これらの機能ブロックは、ハードウェア及びソフトウェアの組み合わせによって、種々の態様で実現することが可能である。

【0021】

処理装置２０は入力装置２８及び出力装置２９と接続される。入力装置２８は、処理装置２０で実行される処理に関係するユーザからのコマンド及びデータの入力を受ける。入力装置２８として、例えばユーザが操作を行うことにより入力を行うキーボードやマウス、インターネット等のネットワークを介して入力を行う通信装置、ＣＤ、ＤＶＤ等の記録媒体から入力を行う読取装置等を用いることができる。

【0022】

解析情報取得部２１は、入力装置２８を介して磁場解析情報を取得する。磁場解析情報には、磁場解析に必要な種々の情報が含まれる。例えば、仮想空間内に定義される解析対象物の形状、解析対象物をボクセルに分割する情報、ボクセルを粗視化する情報、仮想空間内に存在する磁場の空間的な変化が相対的に緩やかな方向を定義する情報等が含まれる。

【0023】

磁化付与部２２は、磁場解析情報に基づき仮想空間内の対象物を複数のボクセル（正六面体要素）に分割する。ボクセル５０の各々に、１つのボクセル５０を特定するためのボクセル識別子（ボクセルＩＤ）が付与される。

【0024】

図２Ａは、対象物の一部分を複数のボクセル５０に分割した状態の斜視図である。複数のボクセル５０が、ｘ軸方向、ｙ軸方向、及びｚ軸方向に並んでいる。複数のボクセル５０の形状及び寸法は同一であり、複数のボクセル５０は、単純立方単位胞構造を持つ。本実施例において、ｚ軸方向に関する磁場の変化が相対的に緩やかであり、ｘ軸方向及びｙ軸方向に関する磁場の変化が相対的に急峻である例について説明する。ここで、磁場には、外部磁場、及び対象物が磁化されることによって発生する磁場が含まれる。

【0025】

さらに、磁化付与部２２（図１）は、ボクセル５０の各々に磁化を付与する。磁化は初期条件として付与してもよいし、後述する磁場演算部２３により求められた磁場に基づいて各ボクセル５０に発生する磁化を計算し、計算によって求められた磁化を付与してもよい。なお、磁化は、ボクセル５０内の全磁気モーメントのベクトル和の単位体積当たりの値であるから、ボクセル５０ごとに磁化として１つの値を付与することは、ボクセル５０に磁気モーメントの１つの値を付与することと実質的に同一である。磁化付与部２２は、付与された磁化を、ボクセルＩＤと関連付けてデータ保持部２６に登録する。

【0026】

さらに、磁化付与部２２は、磁場の空間的な変化が緩やかな方向（本実施例においてｚ軸方向）に連続する複数のボクセル５０をまとめて１つの体積要素５１とする。複数のボクセル５０をまとめて１つの体積要素５１とする処理を粗視化という。粗視化された体積要素５１の各々のｚ軸方向の寸法は、ｘ軸方向及びｙ軸方向の寸法より大きくなる。

【0027】

図２Ｂは、粗視化後の対象物の一部分の斜視図である。ｚ軸方向に並ぶ複数のボクセル５０（図２Ａ）が粗視化されて（まとめられて）体積要素５１が形成されている。１つの体積要素５１に粗視化するボクセル５０の個数は、磁場の変化の度合いに応じて設定するとよい。例えば、磁場の変化が緩やかであるほど、１つの体積要素５１に粗視化するボクセル５０の個数を多くするとよい。図２Ｂでは、上下方向の中央部の近傍で磁場のｚ軸方向の変化が緩やかであり、上下端に近づくに従って、磁場の変化が徐々に大きくなる例を示している。このため、上下方向の中央部の近傍において、１個の体積要素５１に粗視化するボクセル５０の個数を相対的に多くし、上下端に近づくほど、１個の体積要素５１に粗視化するボクセル５０の個数を相対的に少なくしている。

【0028】

磁化付与部２２（図１）は、元のボクセル５０に付与されている磁化に基づいて体積要素５１に、それぞれ磁化を付与する。例えば、１つの体積要素５１に粗視化された元のボクセル５０に付与されている磁化のベクトル平均を、体積要素５１に付与するとよい。磁化付与部２２は、体積要素５１に付与された磁化を、体積要素識別子（体積要素ＩＤ）と関連付けてデータ保持部２６に登録する。

【0029】

磁場演算部２３（図１）は、ボクセル５０に付与された磁化、または体積要素５１に付与された磁化に基づいて、仮想空間内の複数の観測点における磁場を計算し、計算結果をデータ保持部２６に登録する。観測点に磁場を生じさせる元となる磁化を持つボクセル５０または体積要素５１をソースということとする。観測点は、例えば体積要素５１の各々の中心位置に配置する。磁場演算部２３は、磁場を計算する３種類の機能を備えている。第１機能（テーブル参照方式）は、ソースから観測点までの距離が短いときの計算に適しており、対応関係保持部２７に登録されている対応関係を利用する。第２機能（双極子近似方式）は、ソースから観測点までの距離が長いときの計算に適している。第３機能（逐次積分方式）は、観測点がソースに近接しているときの計算に適している。第１機能～第３機能の詳細については、後に詳しく説明する。

【0030】

出力制御部２４（図１）は、解析結果、例えばデータ保持部２６に登録されている観測点ごとの磁場の計算結果を、出力装置２９に出力する。

【0031】

［第１機能（テーブル参照方式）］
次に、図３Ａ～図４を参照して、磁場演算部２３の第１機能（テーブル参照方式）について説明する。
図３Ａは、仮想空間内に定義される座標系を示す図である。座標系の原点をＯで表し、仮想空間内の観測点Ｐｏの位置をｒ_ｉで表す。微小体積ｄＶ’の位置をｒ’で表す。ここで、太文字はベクトルを意味する。観測点Ｐｏの磁場Ｈ（ｒ_ｉ）は、微小体積ｄＶ’の磁化Ｍ（ｒ’）によって発生する磁場を全空間に亘って積分することにより得られ、以下の式で表される。

【数1】

ここで、∇’は、位置ｒ’における微分を意味する。

【0032】

磁性材料が分布する領域をＶで表し、表面をＳで表すと、式（１）は、以下のように体積分と面積分の項に分けることができる。

【数2】

ここで、ｄＳ’は、磁性材料が分布する領域の表面の微小領域であり、ベクトルｎ_ｓは、表面Ｓの単位法線ベクトルである。

【0033】

図３Ｂは、１つの体積要素５１の内部の磁化が一定であるとしたときに、体積要素５１に付与された磁化によって観測点Ｐｏに発生する磁場Ｈ（ｒ_ｉ）を計算する方法を説明する図である。ｊ番目の体積要素５１（位置ｒ_ｊ）に付与された磁化をＭ_ｊで表す。体積要素５１の表面の微小領域をｄＳ_ｊで表し、微小領域ｄＳ_ｊの単位法線ベクトルをｎ_ｓで表す。

【0034】

式（２）を変形して、位置ｒ_ｉの観測点Ｐｏに発生する磁場Ｈ（ｒ_ｉ）は、以下の式で表される。

【数3】

ここで、ｒ_ｊは、微小領域ｄＳ’の位置を表す。パラメータｊに関するΣは、Ｎ個の体積要素５１についての和をとることを意味する。表面Ｓ_ｊに関する面積分は、体積要素５１ごとに表面全域で面積分を行うことを意味する。

【0035】

図４は、ソースとなる体積要素５１Ｓと、その周りの複数の基準点Ｇ_ｋとの位置関係を示す図である。図４では、ソースとなる体積要素５１Ｓが、連続する４個のボクセル５０を含む例を示している。本実施例においては、式（３）の面積分の演算処理を高速化させるために、ソースとなる１つの体積要素５１Ｓの磁化Ｍ_ｊと、複数の基準点Ｇ_ｋ（位置ｒ_ｉ）に発生する磁場Ｈ_ｋ（ｒ_ｉ）との対応関係を、体積要素５１Ｓの表面における面積分を実行することにより予め算出しておく。基準点Ｇ_ｋは、例えば複数のボクセル５０の頂点の位置に対応して配置される。

【0036】

ソースとなる体積要素５１Ｓ（位置ｒ_ｊ）の磁化Ｍ_ｊによって基準点Ｇ_ｋに発生する磁場Ｈ_ｋ（ｒ_ｉ）は、式（３）から、以下のように表すことができる。

【数4】

ここで、面積分は、体積要素５１Ｓの表面の全域についての面積分を意味する。

【0037】

式（４）で表される磁場Ｈ_ｋ（ｒ_ｉ）は、面積分を利用した厳密解である。ここで、式（４）を変形し、磁場Ｈ_ｋ（ｒ_ｉ）と磁化Ｍ_ｊとの対応関係を定めるテンソルＴ_ｋを用いて下記式（５）のように表すことができれば、任意の磁化Ｍ_ｊを持つ体積要素５１によって各基準点Ｇ_ｋに発生する磁場Ｈ_ｋ（ｒ_ｉ）を容易に計算できるようになる。

【数5】

【0038】

テンソルＴ_ｋの各要素は、以下のように表すことができる。

【数6】

ここで、ｎ_ｘ、ｎ_ｙ、ｎ_ｚは、それぞれ単位法線ベクトルｎ_ｓのｘ成分、ｙ成分、ｚ成分である。ｘ_ｉｊ、ｙ_ｉｊ、ｚ_ｉｊは、それぞれソースとなる体積要素５１Ｓの微小領域ｄＳ_ｊを始点とし基準点Ｇ_ｋを終点とするベクトルｒ_ｉｊのｘ成分、ｙ成分、ｚ成分である。

【0039】

本実施例では、ソースとなる体積要素５１Ｓに対して、複数の基準点Ｇ_ｋについて式（６）の面積分を実行し、テンソルＴ_ｋを複数の基準点Ｇ_ｋについて予め算出している。テンソルＴ_ｋは、ソースとなる体積要素５１Ｓに付与された磁化Ｍ_ｊと、磁化Ｍ_ｊによって複数の基準点Ｇ_ｋの各々に発生する磁場Ｈ_ｋとの対応関係を示している。この対応関係は、ソースとなる体積要素５１Ｓと観測点との相対位置ごとに準備されている。さらに、対応関係を示すテンソルＴ_ｋは、体積要素５１（図２Ｂ）の大きさごとに準備されている。例えば、粗視化するボクセル５０の最大個数を１０個とすると、粗視化するボクセルの個数が１個の場合から１０個の場合までのそれぞれについてテンソルＴ_ｋが準備される。

【0040】

対応関係を示すテンソルＴ_ｋは、対応関係保持部２７（図１）に記憶されている。テンソルＴ_ｋを予め計算する具体的な方法の一例が、特開２０１７－１９４８８４号公報に説明されている。なお、式（６）の面積分を予め実行しておく方法は、特開２０１７－１９４８８４号公報に示された例に限定されない。

【0041】

位置ｒ_ｉの観測点Ｐｏ（図３Ｂ）の磁場Ｈ（ｒ_ｉ）は、観測点Ｐｏを内部に含むボクセル５０の８個の頂点に相当する基準点Ｇ_ｋの磁場Ｈ_ｋを平均化することにより算出することができる。

【0042】

テーブル参照方式は、ソースとなる体積要素５１Ｓから観測点Ｐｏまでの距離が、ある基準距離以下のときに利用される。従って、ソースとなる体積要素５１Ｓからの距離がこの基準距離以下の基準点Ｇ_ｋについてテンソルＴ_ｋを準備しておけばよい。テンソルＴ_ｋを準備しておくべき基準点Ｇ_ｋの個数は、仮想空間内の体積要素５１の個数より少ない。このため、全ての体積要素５１が観測点Ｐｏ（図３Ｂ）に発生させる磁場を、その都度面積分を用いて計算する方式に比べて、対応関係を示すテンソルＴ_ｋを利用して磁場を計算するテーブル参照方式を適用する方が、演算負荷が低くなる。

【0043】

テーブル参照方式では、粗視化後の１つの体積要素５１が１つのソースとなり、かつ１つの観測点となる。このため、ｎ個のボクセル５０をまとめて１つの体積要素５１とすると、ソース及び観測点の個数が共に１／ｎになる。ソースと観測点との組み合わせの個数は１／ｎ^２になるため、磁場の演算時間がほぼｎ^２に比例して減少する。

【0044】

［第２機能（双極子近似方式）］
次に、図５を参照して、磁場演算部２３（図１）の第２機能（双極子近似方式）について説明する。

【0045】

図５は、磁化が付与されたボクセル５０と観測点Ｐｏとの位置関係を示す図である。観測点Ｐｏは、例えば体積要素５１の各々の中心に設定される。磁場を計算する元となる磁化として、粗視化前のボクセル５０に付与された磁化を用いる。

【0046】

式（１）を部分積分して近似を行うと、以下の式（７）が得られる。

【数7】

【0047】

双極子近似方式は、粗視化前のボクセル５０に付与された磁化Ｍ_ｊを、ボクセル５０の中心に配置された磁気双極子で代表させて磁場を求めることと等価である。ボクセル５０の体積をΔＶ_ｊで表すと、磁気モーメントはΔＶ_ｊＭ_ｊと表すことができる。双極子近似方式においては、ソース点となる体積要素５１（図２Ａ）の各々について面積分を実行する必要が無いため、演算負荷が軽減される。

【0048】

上述の例では、体積要素５１を構成する複数のボクセル５０の各々に付与された磁化Ｍ_ｊを用いて磁場を計算したが、粗視化後の体積要素５１の中心位置に、複数のボクセル５０に付与された磁気双極子を多重極展開して、多重極子に基づいて観測点Ｐｏの磁場を計算してもよい。

【0049】

双極子近似方式では、粗視化後の１つの体積要素５１に対応して１つの観測点Ｐｏが配置され、粗視化前の１つのボクセル５０が１つのソースとなる。このため、ｎ個のボクセル５０をまとめて１つの体積要素５１に粗視化すると、観測点Ｐｏの個数が１／ｎになり、ソースの個数は不変である。その結果、ソースと観測点との組み合わせの個数が１／ｎになることにより、磁場の演算時間がほぼｎに比例して減少する。

【0050】

［テーブル参照方式と双極子近似方式との比較］
上述のテーブル参照方式と双極子近似方式とを比較すると、計算精度の点ではテーブル参照方式の方が優れており、演算負荷の点では双極子近似方式の方が優れている。なお、ソースとなるボクセル５０から遠く離れた領域の観測点の磁場を計算する際には、ボクセル５０内に分布する磁化を磁気双極子で近似しても十分な計算精度が得られる。

【0051】

従って、ソースとなるボクセル５０と観測点Ｐｏとの距離に基づいて、テーブル参照方式と双極子近似方式とのどちらを適用するか決定するとよい。例えば、ソースを中心としたある半径の球の内部の観測点Ｐｏの磁場の計算にはテーブル参照方式を適用し、この球の外側の観測点Ｐｏの磁場の計算には双極子近似方式を適用するとよい。

【0052】

［第３機能（逐次積分方式）］
次に、図６を参照して、磁場演算部２３（図１）の第３機能（逐次積分方式）について説明する。

【0053】

図６は、ソースとなる体積要素５１Ｓと観測点Ｐｏとの位置関係を示す図である。観測点Ｐｏは、例えば体積要素５１の各々の中心に設定される。ソースとなる体積要素５１Ｓ内の位置ｒ’の磁化をＭ（ｒ’）で表す。観測点Ｐｏの位置をｒ_ｉで表すと、ソースとなる体積要素５１Ｓに付与された磁化Ｍ（ｒ’）が観測点Ｐｏの位置に発生させる磁場Ｈ（ｒ_ｉ）は、式（２）に基づいて計算することができる。ここで、式（２）の体積分は、ソースとなる体積要素５１Ｓ内の微小体積要素ｄＶ’について実行され、面積分は、ソースとなる体積要素５１Ｓの表面の微小領域ｄＳ’について実行される。

【0054】

体積要素５１Ｓ内の磁化Ｍ（ｒ’）が一定であると仮定できる場合、または体積要素５１Ｓ内における磁化Ｍ（ｒ’）の空間的な変化がほぼ無視できる場合、式（２）の右辺の体積分の項が０になるため、式（２）の右辺の面積分を実行することにより、磁場Ｈ（ｒ_ｉ）を求めることができる。

【0055】

逐次積分方式においては、テーブル参照方式で利用した基準点Ｇ_ｋの位置の磁場を平均化することなく、観測点Ｐｏの磁場を直接計算する。このため、テーブル参照方式に比べて、磁場の計算精度を高めることができる。特に、ソース点が観測点に近接しているときに、逐次積分方式を適用する効果が顕著である。

【0056】

また、テーブル参照方式では、ソースとなる体積要素５１Ｓ（図４）内の磁化Ｍ_ｊが一定であるという条件の下で磁場が計算される。逐次積分方式は、体積要素５１内の磁化が空間的に変化している場合にも適用することができる。

【0057】

逐次積分方式では、粗視化後の１つの体積要素５１が１つのソースとなり、かつ１つの観測点となる。このため、ｎ個のボクセル５０をまとめて１つの体積要素５１に粗視化すると、ソース及び観測点の個数が共に１／ｎになる。その結果、ソースと観測点との組み合わせの個数が１／ｎ^２になるため、磁場の演算時間がほぼｎ^２に比例して減少する。

【0058】

次に、図７Ａを参照して、観測点Ｐｏに発生する磁場を計算する一例について説明する。
図７Ａは、ソースＰｓと観測点Ｐｏ１、Ｐｏ２との位置関係を示す図である。ソースＰｓを中心とした半径ａ１の球の内側の第１領域Ｒ１に存在する観測点Ｐｏ１に発生する磁場の計算には、テーブル参照方式を適用する。これにより、計算精度を確保しつつ、演算負荷の軽減を図ることができる。第１領域Ｒ１の外側の第２領域Ｒ２に存在する観測点Ｐｏ２に発生する磁場の計算には、双極子近似方式を適用する。これにより、計算精度を確保しつつ、演算負荷をさらに軽減させることができる。

【0059】

次に、図７Ｂを参照して、観測点Ｐｏに発生する磁場を計算する他の例について説明する。
図７Ｂは、ソースＰｓと、観測点Ｐｏ２、Ｐｏ３、Ｐｏ４との位置関係を示す図である。図７Ｂに示した例では、第１領域Ｒ１が、ソースＰｓを中心とする半径ａ２（ａ２＜ａ１）の球の内側の第３領域Ｒ３と外側の第４領域Ｒ４とに区分されている。第４領域Ｒ４内の観測点Ｐｏ４に発生する磁場の計算には、図７Ａの例と同様にテーブル参照方式を適用する。第３領域Ｒ３内の観測点Ｐｏ３に発生する磁場の計算には、逐次積分方式を適用する。これにより、計算精度をより高めることができる。

【0060】

実施例による磁場解析装置を用いることにより得られる効果を確認するために、実際に磁場解析を行った。次に、図８Ａ～図８Ｃを参照して、この解析結果について説明する。

【0061】

図８Ａ及び図８Ｂは、それぞれ解析対象の磁石電動機に用いられる永久磁石の平面図及び側面図である。ｚ軸を中心軸とする円柱の側面に沿って、円周方向に１０個の永久磁石６０が等間隔に配置されている。半径方向をｘ軸及びｙ軸と定義する。永久磁石６０の各々の平断面は円弧状であり、半径方向に磁化されている。円周方向に隣り合う２つの永久磁石６０の磁化方向は反対向きである。磁場解析において、永久磁石６０を、複数のボクセル５０（図２Ａ）に分割する。図８Ｂには、粗視化後の体積要素５１のｚ軸方向の境界を示しており、ｘ軸方向及びｙ軸方向の境界は示していない。

【0062】

永久磁石６０によって形成される磁場のｚ軸方向の変化は、ｘ軸方向及びｙ軸方向の変化に比べて緩やかである。このため、ｚ軸方向に連続する複数のボクセル５０（図２Ａ）をまとめて体積要素５１（図２Ｂ）とする粗視化を行う。また、ｚ軸方向の中央近傍において磁場の変化が相対的に緩やかであり、ｚ軸方向の両端近傍では磁場の変化が急峻になる。このため、１つの体積要素５１にまとめるボクセル５０の個数を、永久磁石６０のｚ軸方向の中心近傍において相対的に多くし、ｚ軸方向の両端近傍において相対的に少なくした。この磁場解析では、ボクセル５０の個数を２３４，９５２個とし、粗視化後の体積要素５１の個数を２６，０６４個とした。比較のために、粗視化しない比較例、すなわち、すべての体積要素５１の各々が１つのボクセル５０で構成される例についても、磁場解析を行った。

【0063】

また、磁場解析において、図７Ｂの第３領域Ｒ３と第４領域Ｒ４との境界となる球面の半径ａ２を、ボクセル５０の一辺の長さの４倍とし、第４領域Ｒ４と第２領域Ｒ２との境界となる球面の半径ａ１を、ボクセル５０の一辺の長さの１６倍とした。

【0064】

図８Ｃは、実施例による方法で磁場解析を行った結果、及び比較例による方法で磁場解析を行った結果を示すグラフである。横軸は、ｚ軸を中心とした回転方向の角度を単位「ラジアン」で表し、縦軸は磁場の強さを任意単位で表す。図８Ｃのグラフには、ｚ軸方向の中間における永久磁石６０の表面における磁場の強さを示している。グラフ中の実線は実施例による磁場解析結果を示し、丸記号は、比較例による磁場解析結果を示す。

【0065】

図８Ｃに示すように、実施例と比較例との磁場解析結果がよく一致している。これは、粗視化しても十分な計算精度が得られていることを意味する。比較例による方法で磁場解析を行ったときの計算時間は、実施例による方法で磁場解析を行ったときの計算時間の３３．９倍であった。粗視化することにより、演算負荷が軽減され、磁場解析時間の短縮化を図ることが可能であることが確認された。

【0066】

［双極子近似方式への多重極展開の適用］
次に、図９～図１２を参照して、双極子近似方式を用いた磁場の計算に多重極展開を適用して計算時間の短縮を図る方法について説明する。

【0067】

図９は、多重極展開の方法を説明するために用いる座標系を示す図である。磁場の変化が相対的に緩やかな方向に一列に並ぶ複数の立方体のボクセル５０を粗視化して１つの体積要素５１が得られている。多重極展開する複数のボクセル５０の個数をＮで表す。粗視化された体積要素５１の中心の位置をｒ_ｃで表す。位置ｒ_ｃを多重極展開点Ｃということとする。多重極展開点Ｃを原点とし、磁場の変化が相対的に緩やかな方向をｚ軸方向とするｘｙｚ直交座標系を定義する。多重極展開点Ｃに多重極展開される対象となる磁気双極子を持つ複数のボクセル５０はｚ軸上に位置する。ｘｙｚ直交座標系において、ｊ番目のボクセル５０の位置はｒ_ｊ－ｒ_ｃで表され、観測点Ｐｏの位置はｒ_ｉ－ｒ_ｃで表される。ｊ番目のボクセル５０に付与された磁化をＭ_ｊで表し、ｊ番目のボクセル５０の体積をΔＶ_ｊで表す。ｚ軸と、多重極展開点Ｃから観測点Ｐｏに向かうベクトルｒ_ｉ－ｒ_ｃとのなす角度をθ_ｉで表す。

【0068】

ｊ番目のボクセル５０と観測点Ｐｏとの間の距離ｒ_ｉｊは、以下の式（８）で表される。

【数8】

ここで、ｒ_ｉｃ、ｒ_ｊｃは、それぞれ多重極展開点Ｃから観測点Ｐｏまで距離、及び多重極展開点Ｃからｊ番目のボクセル５０までの距離である。Ｐ_Ｌ（ｃｏｓθ_ｉ）は、ｃｏｓθ_ｉを引数とするＬ次のルジャンドル多項式である。なお、式（８）は、J. D. JACKSON, CLASSICAL ELECTRODYNAMICS THIRD EDITION, p.103に説明されている。

【0069】

ベクトルｒ_ｉ－ｒ_ｃのｘ成分であるｘ_ｉｃで式（８）を偏微分すると、以下の式（９）が得られる。

【数9】

ここで、ｒ_ｉｃ＝ｒ_ｉ－ｒ_ｃ、ｒ_ｊｃ＝ｒ_ｊ－ｒ_ｃである。ｘ_ｊｃは、位置ベクトルｒ_ｊｃのｘ成分である。

【0070】

ｊ番目のボクセル５０はｚ軸上に位置するため、ｘ_ｊｃ＝ｙ_ｊｃ＝０である。この条件を式（９）に適用すると、以下の式（１０）が導出される。

【数10】

ここで、ｇ_１、ｇ_２は以下の式で表される。

【数11】

また、ｚ_ｊｃは、ベクトルｒ_ｊｃのｚ成分である。

【0071】

式（１０）をｘ_ｉｃで偏微分した後、－ｘ_ｉｃで割ると、以下の式が得られる。

【数12】

ここで、ｆ_１、ｆ_２、ｆ_３は、以下の式で表される。

【数13】

【0072】

磁気モーメントｍ_ｊを持つｊ番目のボクセル５０が遠方の観測点Ｐｏに作る磁場Ｈ（ｒ_ｉ）は、以下の式で近似することができる。

【数14】

ここで、Ｍ_ｊ及びΔＶ_ｊは、それぞれｊ番目のボクセル５０に付与された磁化及びｊ番目のボクセル５０の体積である。Ｎ個のボクセル５０が粗視化された体積要素５１が遠方の観測点Ｐｏに作る磁場は、式（１４）をＮ個のボクセル５０について足し合わせることにより求まる。式（１４）をＮ個のボクセル５０について足し合わせると、式（７）が得られる。

【0073】

式（１４）に式（１０）及び式（１２）を代入して整理すると、以下の式が得られる。

【数15】

ここで、中カッコ内の１段目、２段目、及び３段目の式は、それぞれｘ成分、ｙ成分、及びｚ成分を表している。ｍ_ｊｚは、磁気モーメントｍ_ｊのｚ成分である。

【0074】

粗視化前のＮ個のボクセル５０に対応する磁気モーメントｍ_ｊが等しいと仮定する。このとき、ｍ_１＝ｍ_２＝・・・＝ｍ_Ｎ＝ｍと表すことができる。また、粗視化前のＮ個のボクセル５０の一辺の長さｄは等しいと仮定する。ｊ番目のボクセル５０の一辺の長さをｄ_ｊで表すと、ｄ_１＝ｄ_２＝・・・＝ｄ_Ｎ＝ｄと表すことができる。観測点ＰｏへのＮ個のボクセルからの寄与を考えると、式（１５）から以下の式が導出される。

【数16】

ここで、ｍ_ｚは、磁気モーメントｍのｚ成分である。係数ｓ_１ ^Ｌ、ｓ_２ ^Ｌ、ｓ_３ ^Ｌは、以下の式で定義される。

【数17】

【0075】

粗視化の対象となる複数のボクセル５０がｚ軸方向にのみ配列しており、一辺の長さｄ_ｊがすべてのボクセル５０において同一である場合、式（１７）に示す係数ｓ_１ ^Ｌ、ｓ_２ ^Ｌ、ｓ_３ ^Ｌは、ボクセル５０を粗視化して体積要素５１を定義する段階で計算することが可能である。式（１６）には、位置ベクトルとしてｒ_ｉｊが現れておらず、ｒ_ｉｃのみが現れている。このため、粗視化後の体積要素５１が観測点Ｐｏに作る磁場は、多重極展開点Ｃの位置ｒ_ｃと、観測点Ｐｏの位置ｒ_ｉとを用いて計算することが可能になる。これは、観測点Ｐｏの磁場を計算するにあたり、粗視化前のＮ個のボクセル５０を個別に扱う必要がないことを意味する。従って、Ｎ個のボクセル５０の各々の位置ｒ_ｊごとに磁場を計算する手法と比べて、計算時間を短縮することが可能である。

【0076】

従って、磁場の変化が相対的に緩やかな一方向（ｚ軸方向）に一列に配列した複数のボクセル５０の各々の磁気双極子を多重極展開点Ｃに多重極展開する手法を用いることにより、計算時間を短縮することが可能である。なお、計算時間の短縮を図るために、１つの体積要素５１に粗視化されるＮ個の複数のボクセル５０の一辺の長さｄ_ｊが等しく、Ｎ個のボクセル５０にそれぞれ付与される磁化Ｍ_ｊが等しいという条件を満たすようにシミュレーション条件を設定することが好ましい。

【0077】

係数ｓ_１ ^Ｌ、ｓ_２ ^Ｌ、ｓ_３ ^Ｌの値は、実際にシミュレーション計算を行うときに用いるルジャンドル多項式の展開次数まで予め計算しておき、記憶装置２５（図１）に記憶させておくとよい。

【0078】

次に、多重極展開手法を用いることの効果を確認するために行った評価実験の結果について説明する。多重極展開の手法を用いた磁場計算と、個々のボクセル５０について双極子近似を適用した磁場計算とを実際に行って両者の精度及び計算時間を比較した。

【0079】

図１０Ａは、個々のボクセル５０について双極子近似を行い、観測点Ｐｏに発生する磁場Ｈ（ｒ_ｉ）を計算により求める方法を説明する図である。１００個のボクセル５０がｚ軸方向に並んで配置されている。ボクセル５０の各々に磁化Ｍが付与されている。個々のボクセル５０が観測点Ｐｏに発生させる磁場を求め、１００個のボクセル５０について足し合わせることにより、磁場Ｈ（ｒ_ｉ）を求めた。この計算によって求められた磁場をＨｏと表記することとする。

【0080】

図１０Ｂは、多重極展開手法を用い、観測点Ｐｏに発生する磁場Ｈ（ｒ_ｉ）を計算により求める方法を説明する図である。１００個のボクセル５０（図１０Ａ）を粗視化して体積要素５１を定義し、体積要素５１の中心である多重極展開点Ｃに、１００個のボクセル５０の磁気モーメントを多重極展開した。立方体のボクセル５０の各々の一辺の長さをｄで表す。このとき、粗視化スケールは１００ｄになる。式（１７）を用いて係数ｓ_１ ^Ｌ、ｓ_２ ^Ｌ、ｓ_３ ^Ｌを予め計算しておき、式（１６）を用いて観測点Ｐｏの磁場Ｈ（ｒ_ｉ）を求めた。この計算によって求められた磁場をＨｃと表記することとする。磁場Ｈｃの計算は、ルジャンドル多項式の展開次数を０次、２次、４次とした３つの場合について行った。

【0081】

図１１は、磁場Ｈｏに対する磁場Ｈｃの誤差を、規格化距離の関数として示すグラフである。規格化距離は、多重極展開点Ｃから観測点Ｐｏまでの距離を粗視化スケール１００ｄで規格化したものである。誤差は、（Ｈｃ－Ｈｏ）／Ｈｏで定義される。図１１のグラフに示した細い実線、破線、及び太い実線は、それぞれルジャンドル多項式の展開次数を０次、２次、及び４次とした場合の誤差を示している。

【0082】

ルジャンドル多項式の展開次数を高くするほど、誤差が小さくなっていることがわかる。また、規格化距離が長くなるほど、誤差が小さくなっていることがわかる。図７Ａ及び図７Ｂに示したように、ソースＰｓから観測点Ｐｏまでの距離が短い範囲で逐次積分方式やテーブル参照方式を用いる場合は、ソースＰｓからある程度遠い範囲において双極子近似方式が適用される。このため、ルジャンドル多項式の展開次数を低くしても、十分高い精度が得られると考えられる。

【0083】

図１２は、多重極展開を行うことによる計算の高速化効果を示すグラフである。個々のボクセル５０について双極子近似を行う手法（図１０Ａ）を適用した場合の計算時間をＴ_０で表し、多重極展開手法を適用した場合の計算時間をＴ_Ｃで表したとき、高速化効果をＴ_０／Ｔ_Ｃと定義する。例えば、計算時間Ｔ_Ｃが計算時間Ｔ_０の１／２になった場合、高速化効果は２倍である。

【0084】

ルジャンドル多項式の展開次数を低くするごとに、高速化効果が高くなっていることがわかる。ルジャンドル多項式の展開次数を２次までとした場合、約２５倍の高速化効果が得られている。また、ルジャンドル多項式の展開次数を４次までとした場合でも、約５倍の高速化効果が得られている。

【0085】

図１１及び図１２に示した評価実験により、多重極展開手法を適用することによって、所望の精度を維持しつつ、計算時間を短縮することが可能であることが確認された。ルジャンドル多項式の展開次数を低くすると計算時間は短くなるが、精度は低下する。求められる精度と、許容される計算時間に応じて、ルジャンドル多項式の展開次数を決めるとよい。

【0086】

図１１に示したように、規格化距離が長くなるとルジャンドル多項式の展開次数を低くしても高い精度が得られる。このため、規格化距離に応じてルジャンドル多項式の展開次数を変化させてもよい。規格化距離が長くなるに従って、ルジャンドル多項式の展開次数を低くするとよい。この手法により、求められる精度を維持しつつ、計算時間をさらに短くすることが可能になる。

【0087】

次に、上記実施例の変形例について説明する。
上記実施例では、図２Ｂに示したｚ軸方向に関する磁場の変化が、ｘ軸方向及びｙ軸方向に関する磁場の変化より緩やかである場合を示した。この例は、磁化されている物体がｚ軸方向に長い棒状または筒状の形状を持つような場合に相当する。その他に、磁化されている物体がｘｙ面内に広がる板状の形状を持つような場合には、ｘ軸方向及びｙ軸方向に関する磁場の変化が、ｚ軸方向に関する磁場の変化より緩やかになる。このような場合には、体積要素５１のｘ軸方向及びｙ軸方向の寸法を、ｚ軸方向の寸法より大きくするとよい。

【0088】

上記実施例による解析手法は、ボクセル５０（図２Ａ）に内接する球形のビーズを配置し、ビーズに磁化（磁気モーメント）を付与し、ビーズの集合体からなる磁性材料の運動をシミュレーションする方法（磁気ビーズ法）にも適用することができる。磁気ビーズ法については、特開２０１５－１１１４０１号公報、特開２０１７－１９４８８４号公報に説明されている。

【0089】

式（７）に示した体積ΔＶ_ｊを磁気ビーズの体積と考えると、ボクセル５０の体積は、磁気ビーズの充填率γ_ｊを用いて、ΔＶ_ｊ／γ_ｊで表される。従って、本実施例による手法を磁気ビーズ法に適用する場合、式（７）のΔＶ_ｊに代えて、ΔＶ_ｊ／γ_ｊを用いるとよい。

【0090】

図８Ａ～図８Ｃに示したシミュレーションでは、磁石電動機の永久磁石が磁化されている例を取り扱ったが、上記実施例による解析手法は、その他の電磁アクチュエータの磁場解析に利用することができる。また、対象物が永久磁石ではなく、鉄等の強磁性材料である場合にも、上記実施例による解析手法を適用することができる。この場合には、各ボクセル５０（図２Ａ）には、外部磁場によって生じる磁化を付与し、この磁化による磁場を求めればよい。

【0091】

上記実施例では、解析対象をボクセル５０（図２Ａ）に分割し、連続する複数のボクセル５０を粗視化して体積要素５１を形成したが、ボクセル５０を経ることなく、解析対象を直接体積要素５１に分割してもよい。

【0092】

上記実施例では、１つの体積要素５１（図２Ｂ）内に１つの観測点を設定したが、１つの体積要素５１内に複数の観測点を設定してもよいし、ボクセル５０ごとに観測点を設定してもよい。１つの体積要素５１内に複数の観測点を設定し、観測点に発生する磁場に基づいて磁化を生じさせると、１つの体積要素５１内に磁化の分布が生じる。テーブル参照方式を適用する場合には、１つの体積要素５１内の磁化の分布を平均化した磁化Ｍ_ｊを用いて磁場を計算するとよい。双極子近似方式を適用する場合には、観測点の磁場から、ソースとなるボクセル５０（図５）ごとの磁化Ｍ_ｊを求めるとよい。逐次積分方式を適用する場合には、ソースとなる体積要素５１内の均一ではない磁化Ｍ（ｒ’）に基づいて、観測点の磁場を計算すればよい。

【0093】

上述の実施例は例示であり、異なる例で示した構成の部分的な置換または組み合わせが可能であることは言うまでもない。複数の例の同様の構成による同様の作用効果については実施例ごとには逐次言及しない。さらに、本発明は上述の実施例に制限されるものではない。例えば、種々の変更、改良、組み合わせ等が可能なことは当業者に自明であろう。

【符号の説明】

【0094】

２０ 処理装置
２１ 解析情報取得部
２２ 磁化付与部
２３ 磁場演算部
２４ 出力制御部
２５ 記憶装置
２６ データ保持部
２７ 対応関係保持部
２８ 入力装置
２９ 出力装置
３０ 磁場解析装置
５０ ボクセル
５１ 体積要素
５１Ｓ ソース点となる体積要素
６０ 永久磁石

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】