(57)【特許請求の範囲】
【請求項１】
解析対象物の物性情報を受け付ける入力部と、
前記入力部から受け付けた前記解析対象物の物性情報に基づいて、当該解析対象物から発生する電流ベクトル及び磁化ベクトルを含む方程式から磁化ベクトルの項を削除した、電流ベクトルによる第一の磁界を算出するための第一の方程式と、前記削除された前記磁化ベクトルの項を含む方程式を、ベクトルポテンシャルを未知数とし、透磁率を含む方程式を空気の透磁率：１に仮定して前記透磁率を一定にすることでポワソン方程式に変形し、当該解析対象物から発生する前記磁化ベクトルによる第二の磁界を算出するための第二の方程式とを、有限要素法および境界積分法により解析可能な形式で作成する方程式作成部と、
前記方程式作成部によって作成された方程式を用いて、前記第一の磁界および前記第二の磁界それぞれを算出する磁界計算部と、
を有したことを特徴とする磁界解析装置。
【請求項２】
前記入力部は、前記解析対象物が有する物体それぞれの位置情報をさらに受け付け、
前記方程式作成部は、前記入力部から受け付けた前記解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、前記第一の方程式に陰解法または陽解法を適用することで、前記第一の磁界の時系列に沿った値を算出するための複数の方程式である第一の連立方程式をさらに作成し、
前記磁界計算部は、前記方程式作成部によって作成された前記第一の連立方程式を用いて、前記第一の磁界を算出することを特徴とする請求項１に記載の磁界解析装置。
【請求項３】
前記入力部は、前記解析対象物が有する物体それぞれの位置情報をさらに受け付け、
前記方程式作成部は、前記入力部から受け付けた前記解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、前記第二の方程式に陰解法または陽解法を適用することで、前記第二の磁界の時系列に沿った値を算出するための複数の方程式である第二の連立方程式をさらに作成し、
前記磁界計算部は、前記方程式作成部によって作成された前記第二の連立方程式を用いて、前記第二の磁界を算出することを特徴とする請求項１に記載の磁界解析装置。
【請求項４】
前記磁界計算部は、算出した前記第一の磁界および前記第二の磁界を合算し、合算した合算磁界を前記解析対象物における磁界の解析結果と決定することを特徴とする請求項２または３に記載の磁界解析装置。
【請求項５】
解析対象物の物性情報を受け付ける入力手順と、
前記入力手順から受け付けた前記解析対象物の物性情報に基づいて、当該解析対象物から発生する電流ベクトル及び磁化ベクトルを含む方程式から磁化ベクトルの項を削除した、電流ベクトルによる第一の磁界を算出するための第一の方程式と、前記削除された前記磁化ベクトルの項を含む方程式を、ベクトルポテンシャルを未知数とし、透磁率を含む方程式を空気の透磁率：１に仮定して前記透磁率を一定にすることでポワソン方程式に変形し、当該解析対象物から発生する前記磁化ベクトルによる第二の磁界を算出するための第二の方程式とを、有限要素法および境界積分法により解析可能な形式で作成する方程式作成手順と、
前記方程式作成部によって作成された方程式を用いて、前記第一の磁界および前記第二の磁界それぞれを算出する磁界計算手順と、
をコンピュータに実行させることを特徴とする磁界解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、磁界解析装置および磁界解析プログラムに関する。
【背景技術】
【０００２】
従来より、磁気記録装置のヘッドやモータなどの電磁機器の設計に際して、磁性体や導体などにどのような磁界が発生するのかを解析する磁界シミュレーションが行われている。磁界シミュレーションでは、解析領域を小さな領域に分割（メッシュ分割）し、分割した一つ一つの領域における磁化状態を求めるための方程式が作成される。そして、磁界シミュレーションでは、作成した領域ごとの方程式を連立方程式として解くことにより、電磁機器全体の磁化状態の挙動が予測される。方程式の作成や連立方程式の解法としては、例えば、有限要素法などの数値解析法が用いられる。
【０００３】
ここで、上記した有限要素法を用いた数値解析法では、磁性体や導体などの隙間に存在する空間領域に対してもメッシュ分割を行い、分割した一つ一つの空間領域における磁化状態を求めるための方程式を含めた連立方程式を解く。従って、上記した有限要素法を用いた数値解析法は、磁性体や導体が移動することにより空間が変化する場合には、移動するごとに方程式を作成し、連立方程式を解く必要があるため、解析に長時間要する場合がある。
【０００４】
そこで、磁性体や導体が移動することにより空間が変化する場合でも、空間領域のメッシュ分割を行なうことなく、磁界シミュレーションに適用できる数値解析法として、境界要素法を用いた数値解析法が知られている（例えば、特許文献１または２参照）。具体的には、境界要素法を用いた数値解析法では、有限要素法と境界要素法とのハイブリッド法により、複数の磁性体の一つ一つを任意の方向に移動しながらの磁化状態が解析される。しかしながら、上記した境界要素法を用いた数値解析法では、多数の連立方程式を解く必要があり、解析に長時間要する場合がある。
【０００５】
そこで、多数の連立方程式が不要となり、容易に磁界シミュレーションが行える数値解析法として、有限要素法および境界積分法を用いた数値解析法が知られている。
【０００６】
例えば、有限要素法および境界積分法を用いた数値解析法では、磁性体などから発生する磁化ベクトルによる磁界と導体などから発生する電流ベクトルによる磁界とをそれぞれ有限要素法および境界積分法により算出する。そして、有限要素法および境界積分法を用いた数値解析法では、算出したそれぞれの磁界から磁化状態を解析する解析方法が知られている（例えば、非特許文献１参照）。しかしながら、上記した数値解析法では、磁化ベクトルによる磁界を算出する際に、スカラーポテンシャルを用いており、磁性体同士が相対的な運動を伴う場合に発生する渦電流や、電流が時間変化する場合に発生する渦電流を考慮することができない。
【０００７】
そこで、近年、渦電流を考慮した有限要素法および境界積分法による数値解析法が知られている（例えば、非特許文献２参照）。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００６−１８３９３号公報
【特許文献２】特開２００７−２１３３８４号公報
【非特許文献】
【０００９】
【非特許文献１】Hybrid method for computing demagnetizing fields.  IEEE transactions on Magnetics, vol.26, No.2, March (1990)
【非特許文献２】Three-dimensional micromagnetic finite element simulations including eddy currents.  J. Appl. Phys. 97, 10E311 (2005)
【発明の概要】
【発明が解決しようとする課題】
【００１０】
ところで、上記した従来技術では、磁化ベクトルによる磁界への影響を正確に考慮できないという課題があった。すなわち、上記した従来技術では、導体などによって生じる電流ベクトルによる磁界に対しては有限要素法および境界積分法が適用されているが、磁性体などによって生じる磁化ベクトルを含む方程式には有限要素法および境界積分法が適用できず、磁化ベクトルによる磁界への影響を正確に考慮できない。
【００１１】
具体的には、有限要素法および境界積分法を適用するためには、方程式の両辺の体積積分を表面積分に変形できなければならないが、磁化ベクトルを含む方程式では、表面積分に変形することができないために、有限要素法および境界積分法が適用できない。従って、上記した従来技術では、磁化ベクトルによる磁界への影響を正確に考慮できない。
【００１２】
そこで、この発明は、上述した従来技術の課題を解決するためになされたものであり、磁化ベクトルによる磁界への影響を正確に考慮することが可能となる磁界解析装置および磁界解析プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１３】
上述した課題を解決し、目的を達成するため、この装置は、入力部が、解析対象物の物性情報を受け付ける。そして、方程式作成部は、入力部から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する電流ベクトルによる第一の磁界を算出するための第一の方程式を有限要素法および境界積分法により解析可能な形式で作成する。また、方程式作成部は、入力部から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する磁化ベクトルによる第二の磁界を算出するための第二の方程式を有限要素法および境界積分法により解析可能な形式で作成する。そして、磁界計算部は、方程式作成部によって作成された方程式を用いて、第一の磁界および第二の磁界それぞれを算出する。
【発明の効果】
【００１４】
開示の装置は、磁化ベクトルによる磁界への影響を正確に考慮することが可能となる。
【図面の簡単な説明】
【００１５】
【図１】図１は、実施例１における磁界解析装置の構成を説明するためのブロック図である。
【図２】図２は、実施例２における磁界解析装置の構成を説明するためのブロック図である。
【図３】図３は、数値解析法の概要を説明するための図である。
【図４】図４は、解析対象モデルを説明するための図である。
【図５】図５は、解析条件データを説明するための図である。
【図６】図６は、実施例２における磁界解析装置による磁界解析処理の手順を説明するための図である。
【図７】図７は、解析例の結果を説明するための図である。
【図８】図８は、実施例２の磁界解析プログラムを実行するコンピュータを説明するための図である。
【発明を実施するための形態】
【００１６】
以下に添付図面を参照して、本願の開示する磁界解析装置および磁界解析プログラムの実施例を詳細に説明する。なお、以下では、本願の開示する磁界解析プログラムを実施する磁界解析装置を実施例として説明する。
【実施例１】
【００１７】
まず最初に、図１を用いて、実施例１における磁界解析装置の構成を説明する。図１は、実施例１における磁界解析装置の構成を説明するためのブロック図である。
【００１８】
図１に示すように、実施例１における磁界解析装置１０は、方程式作成部２２と、磁界計算部２４と、入力部４０とを有する。
【００１９】
入力部４０は、解析対象物の物性情報を受け付ける。
【００２０】
方程式作成部２２は、入力部４０から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する電流ベクトルによる第一の磁界を算出するための第一の方程式を、有限要素法および境界積分法により解析可能な形式で作成する。
【００２１】
また、方程式作成部２２は、入力部４０から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する磁化ベクトルによる第二の磁界を算出するための第二の方程式を、有限要素法および境界積分法により解析可能な形式で作成する。
【００２２】
磁界計算部２４は、方程式作成部２２によって作成された方程式を用いて、第一の磁界および第二の磁界それぞれを算出する。
【００２３】
［実施例１の効果］
上記したように、実施例１によれば、入力部４０は、解析対象物の物性情報を受け付ける。そして、方程式作成部２２は、入力部４０から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する電流ベクトルによる第一の磁界を算出するための第一の方程式を有限要素法および境界積分法により解析可能な形式で作成する。また、方程式作成部２２は、入力部４０から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する磁化ベクトルによる第二の磁界を算出するための第二の方程式を有限要素法および境界積分法により解析可能な形式で作成する。そして、磁界計算部２４は、方程式作成部２２によって作成された方程式を用いて、第一の磁界および第二の磁界それぞれを算出する。従って、磁化ベクトルを含む方程式が体積積分の式から表面積分の式に変換可能となり、その結果、磁化ベクトルを含む方程式に有限要素法および境界積分法を適用することができ、磁化ベクトルの影響を正確に考慮することが可能となる。
【実施例２】
【００２４】
実施例２では、上述した実施例１における磁界解析装置１０が実行する処理内容の具体的な一例について説明する。
【００２５】
［実施例２における磁界解析装置の構成］
まず最初に、図２〜５を用いて、実施例２における磁界解析装置の構成を説明する。図２は、実施例２における磁界解析装置の構成を説明するための図であり、図３は数値解析法の概要を説明するための図であり、図４は、解析対象モデルを説明するための図である。また、図５は、解析条件データを説明するための図である。
【００２６】
図２に示すように、実施例２における磁界解析装置１０は、処理部２０と、記憶部３０と、入力部４０と、出力部５０と、入出力制御Ｉ／Ｆ部６０とを有する。
【００２７】
入力部４０は、キーボードやマウスなどを有し、利用者からの各種情報を受け付ける。具体的には、入力部４０は、利用者が所望する解析対象物および解析条件などの各種情報を受け付け、受け付けた解析対象物および解析条件を後述する処理部２０および記憶部３０に転送する。
【００２８】
出力部５０は、モニタなどを有し、後述する処理部２０の処理結果をモニタに出力する。
【００２９】
入出力制御Ｉ／Ｆ部６０は、処理部２０と、記憶部３０と、入力部４０と、出力部５０との間で送受信されるデータを中継するインターフェースである。
【００３０】
処理部２０は、記憶部３０が記憶するデータを用いて、解析対象物の磁化状態を解析する。以下、記憶部３０が記憶するデータを用いて処理部２０が実行する解析内容について、数式を用いて詳細に説明する。
【００３１】
本実施例における磁界解析に用いる数値解析法を説明する。本実施例における数値解析法は、解析対象物の導体や磁性体の領域のみをメッシュ分割し、空間領域をメッシュ分割することなく、解析対象物全体の磁化状態を解析する。例えば、本実施例における数値解析法では、図３に示すように、空間を隔てて配置された３つの領域のみをそれぞれメッシュ分割し、３つの領域全体の磁化状態を解析することができる。なお、図３に示す「ベクトルＨcurrent、φ、σ、ベクトルｍ、ベクトルＩ」それぞれは、「電流による磁界、スカラーポテンシャル、導電率、磁化ベクトル、電流ベクトル」を示す。また、図３に示す領域１は「導体」、領域２は、「磁化ベクトルが存在する導体」、領域３は、「電流ベクトルが存在する導体」である。また、図３に示す「Γ１、Γ２、Γ３」はそれぞれ「領域１の最外周面、領域２の最外周面、領域３の最外周面」を示し、「全外周面：Γ」を「Γ１+Γ２+Γ３」として表す。
【００３２】
以下、図３に示す記号などを用いた数式より、本実施例における処理部２０が実行する数値解析法である有限要素法および境界積分法（以下、「有限要素―境界積分法」と記す）について説明する。なお、以下では、本実施例によらない解析方法で用いられていた数式について説明したのち、本実施例の解析方法で用いられる数式について説明する。
【００３３】
まず、有限要素―境界積分法を用いた本実施例によらない解析方法は、上記した３つの領域に対して、磁化ベクトルによる磁界の算出にスカラーポテンシャルを使用し、電流ベクトルによる磁界の算出にはベクトルＨcurrentを使用する計算方法である。スカラーポテンシャルを使用する場合、上記した３つの領域に対して有限要素―境界積分法を適用すると、以下に示す式（１）〜（５）が得られる。
【００３４】
【数１】
0005412982.tif 000002
【００３５】
【数２】
0005412982.tif 000003
【００３６】
【数３】
0005412982.tif 000004
【００３７】
【数４】
0005412982.tif 000005
【００３８】
【数５】
0005412982.tif 000006
【００３９】
ここで、式（１）は、有限要素―境界積分法を適用するためにφをφ１およびφ２に分離した式であり、式（２）は、磁化ベクトルに依存する項をもつφ１のポワソン方程式であり、式（３）は、φ１の境界積分よりφ２の境界値を計算する式である。また、式（４）は、φ２を未知数とするラプラス方程式であり、式（５）は、スカラーポテンシャルの勾配により磁界を計算する方程式である。
【００４０】
そして、上記した方程式（２）〜（５）を順番に解くことにより、磁化ベクトルの定常磁界を計算することができる。
【００４１】
また、ベクトルＨcurrentを使用する場合、上記した３つの領域に対して有限要素―境界積分法を適用すると、以下に示す式（６）〜（９）が得られる。
【００４２】
【数６】
0005412982.tif 000007
【００４３】
【数７】
0005412982.tif 000008
【００４４】
【数８】
0005412982.tif 000009
【００４５】
【数９】
0005412982.tif 000010
【００４６】
ここで、式（６）は、電流ベクトルに依存する項を右辺に持つＨ１のポワソン方程式であり、式（７）は、Ｈ１の境界積分よりＨ２の境界値を計算する境界積分の式であり、式（８）は、Ｈ２のラプラス方程式である。また、式（９）は、磁界ベクトルの和により電流磁界を計算する方程式であり、上記した方程式（６）〜（９）を順番に解くことで、電流ベクトルによる磁界を計算することができる。なお、Ｈ１およびＨ２は、有限要素―境界積分法を適用するためにＨを分離したものである。また、式（３）および式（７）は、以下に示す式（１０）および（１１）に書き表すことができ、節点未知数に対して、それぞれ以下の式（１２）および式（１３）に変換することができる。
【００４７】
【数１０】
0005412982.tif 000011
【００４８】
【数１１】
0005412982.tif 000012
【００４９】
【数１２】
0005412982.tif 000013
【００５０】
【数１３】
0005412982.tif 000014
【００５１】
すなわち、Ｂijを含んだ式に変換することにより、節点未知数に対する計算をより簡略化させることが可能である。以下、ＢijをＢマトリックスと記す。
【００５２】
しかしながら、上記した計算方法では、磁性体同士が相対的な運動を伴う場合や電流が時間変化する場合に発生する渦電流などの過渡的効果を考慮した磁界解析を行うことができない。別の解析方法では、過渡的効果を考慮した以下に示す式（１４）が用いられる。なお、式（１４）のベクトルＨeddyは渦電流による磁界である。
【００５３】
【数１４】
0005412982.tif 000015
【００５４】
ところで、上記した式（１４）を用いた磁界解析では、磁化ベクトルを含む両辺の体積積分を表面積分に変形することができないために、有限要素―境界積分法を適用することができず、磁化ベクトルによる磁界への影響を正確に考慮することができない。
【００５５】
そこで、本実施例の数値解析法では、式（１４）に有限要素―境界積分法を適用するために、まず、式（１４）から磁化ベクトルの項を削除した以下の式（１５）を設定する。
【００５６】
【数１５】
0005412982.tif 000016
【００５７】
そして、式（１５）に対して有限要素―境界積分法および陰解法を適用することにより、以下に示す、式（１６）〜（１９）が得られる。
【００５８】
【数１６】
0005412982.tif 000017
【００５９】
【数１７】
0005412982.tif 000018
【００６０】
【数１８】
0005412982.tif 000019
【００６１】
【数１９】
0005412982.tif 000020
【００６２】
ここで、式（１６）〜（１９）を順番に解くことにより、「時間：ｎ+１」における電流ベクトルによる過渡磁界が計算できる。なお、「ｎ」は現在の時間ステップであり、「ｎ+１」は求めたい次の時間ステップである。また、上記では、時間解法として陰解法を適用させた場合について説明したが、本実施例はこれに限定されるものではなく、例えば、陽解法を適用する場合であってもよい。
【００６３】
続いて、本実施例における数値解析法では、式（１４）から式（１５）を導出する際に削除した磁化ベクトルの項について、ベクトルポテンシャルを未知数とする以下の式（２０）から有限要素―境界積分法が適用可能な式を導出する。
【００６４】
【数２０】
0005412982.tif 000021
【００６５】
なお、「Ａ、μ」はそれぞれ「ベクトルポテンシャル、透磁率」である。
【００６６】
そして、式（２０）に対して「空気の透磁率：１」を仮定すると、空気中は透磁率μ０が一定であることから、μ０をローテーションの外に出すことが可能となる。そして、ベクトル公式およびクーロンゲージ（∇・Ａ＝０）を使用することにより、式（２０）が以下の式（２１）となり、式（２２）が得られる。なお、「空気の透磁率：１」は、磁性体を磁化ベクトルで表現することに相当する。
【００６７】
【数２１】
0005412982.tif 000022
【００６８】
【数２２】
0005412982.tif 000023
【００６９】
ここで、式（２０）を式（２２）のようなポワソン方程式に変形することで、有限要素―境界積分法を適用することが可能となる。
【００７０】
そして、式（２２）に対して有限要素―境界積分法および陰解法を適用することにより、以下に示す、式（２３）〜（２７）が得られる。
【００７１】
【数２３】
0005412982.tif 000024
【００７２】
【数２４】
0005412982.tif 000025
【００７３】
【数２５】
0005412982.tif 000026
【００７４】
【数２６】
0005412982.tif 000027
【００７５】
【数２７】
0005412982.tif 000028
【００７６】
ここで、式（１５）および式（２４）は、式（１１）と、以下に示す式（２８）とに書き表すことができ、節点未知数に対して、それぞれＢマトリックスを含んだ式（１３）と、式（２９）とに変換することができる。
【００７７】
【数２８】
0005412982.tif 000029
【００７８】
【数２９】
0005412982.tif 000030
【００７９】
ここで、式（２３）〜（２７）を順番に解くことにより、時間「ｎ+１」における磁化ベクトルによる磁界が計算できる。なお、上記では、時間解法として陰解法を適用させた場合について説明したが、本実施例はこれに限定されるものではなく、例えば、陽解法を適用する場合であってもよい。
【００８０】
そして、式（１９）と式（２７）とを合算する以下の式（３０）を用いることにより、解析対象物全体の磁化状態を解析することができる。
【００８１】
【数３０】
0005412982.tif 000031
【００８２】
すなわち、本実施例における磁界解析装置１０の処理部２０は、式（１５）〜式（３０）を用いて解析対象物の磁界を解析する。
【００８３】
図２に戻って、記憶部３０は、解析対象物をモデル化した解析対象モデル、解析条件データ、各種方程式および後述する処理部２０の処理結果を記憶し、モデルデータ記憶部３１、解析条件データ記憶部３２、方程式記憶部３３、計算結果記憶部３４および結果データ記憶部３５を有する。
【００８４】
モデルデータ記憶部３１は、入力部４０を介して利用者から受け付けた解析対象モデルを記憶する。具体的には、モデルデータ記憶部３１は、磁界解析の解析対象となる電磁機器などの解析対象モデルを、解析対象モデルが有する各材料の物性情報や位置情報とともに記憶する。例えば、モデルデータ記憶部３１は、図４に示すように、非磁性の導体である材料Ａおよびコイルである材料Ｂそれぞれを小さな領域に分割（メッシュ分割）した解析対象モデルを記憶する。また、モデルデータ記憶部３１は、解析対象モデルが有する各材料それぞれの形態および配置状態の３次元空間における３次元情報を記憶する（図４に示すＸ−Ｙ−Ｚの３軸参照）。なお、解析対象モデルのメッシュ分割は、数十〜数百ｎｍオーダのサイズで、利用者によって設定される。
【００８５】
解析条件データ記憶部３２は、解析対象モデルを解析する際に、入力部４０を介して利用者から受け付けた詳細な解析条件データを記憶する。具体的には、解析条件データ記憶部３２は、利用者が解析対象モデルに設定する種々のパラメータを記憶する。例えば、解析条件データ記憶部３２は、図５の（Ａ）に示すように、「コイル電流の向き」、「コイルと導体の端が一致」、Ｚ軸に対する「コイル長：１０μｍ」、Ｚ軸に対する「導体長：３０μｍ＝２０μｍ＋１０μｍ」など、解析対象モデルを解析する際に用いられるパラメータを記憶する。
【００８６】
また、解析条件データ記憶部３２は、図５の（Ａ）に示すように、「解析位置：Ｚ＝２５μｍ、１５μｍおよび５μｍ」と設定されたパラメータを記憶する。さらに、解析条件データ記憶部３２は、図５の(Ｂ)に示すように、「材料：Ａ、透磁率：１、導電率：ｌｅ７、電流：０．０、磁化：０．０」や「材料：Ｂ、透磁率：１、導電率：ｌｅ７、電流：１．０、磁化：０．０」と設定されたパラメータを記憶する。
【００８７】
なお、電流および磁化は、解析開始からの経過時間に対して任意に設定される。例えば、図５の（Ｃ）に示すように、電流は、時間〔ｓｅｃ〕に対して一定値Ｉ〔Ａ〕で設定される。
【００８８】
方程式記憶部３３は、後述する処理部２０によって用いられる方程式と、後述する処理部２０によって作成された方程式とを記憶する。なお、方程式記憶部３３が記憶する各方程式については、後に詳述する。
【００８９】
計算結果記憶部３４および結果データ記憶部３５は、それぞれ後述する処理部２０によって算出された計算結果および解析結果を記憶する。なお、計算結果記憶部３４および結果データ記憶部３５それぞれが記憶する計算結果および解析結果については後に詳述する。
【００９０】
処理部２０は、入力部４０から解析対象物および解析条件の情報を受信し、受信した各種情報に基づいて、磁界解析の各種処理を実行し、入力情報処理部２１、方程式作成部２２と、Ｂマトリックス計算部２３と、磁界計算部２４とを有する。
【００９１】
入力情報処理部２１は、解析対象モデルおよび解析条件データを記憶部３０から読み出し、読み出したデータを後述する方程式作成部２２に転送し、解析用モデル読み出し部２１ａおよび解析条件読み出し部２１ｂを有する。
【００９２】
解析用モデル読み出し部２１ａは、利用者によって指定された解析対象モデルをモデルデータ記憶部３１から読み出し、読み出した解析対象モデルを後述する方程式作成部２２に転送する。
【００９３】
解析条件読み出し部２１ｂは、利用者によって指定された各種パラメータである解析条件データを解析条件データ記憶部３２から読み出し、読み出した解析条件データを後述する方程式作成部２２に転送する。
【００９４】
方程式作成部２２は、入力情報処理部２１から受信した解析対象モデルおよび解析条件データに基づいて、解析対象モデルの磁界解析を行うための方程式を作成し、第一連立方程式作成部２２ａと第二連立方程式作成部２２ｂとを有する。
【００９５】
第一連立方程式作成部２２ａは、解析用モデル読み出し部２１ａから受信した解析対象モデルの物性情報に基づいて、解析対象モデルから発生する電流ベクトルによる磁界である第一磁界を算出するための第一方程式を作成する。具体的には、第一連立方程式作成部２２ａは、解析対象モデルの導体から発生する電流ベクトルによる第一磁界を算出するために、解析用モデル読み出し部２１ａから受信した解析対象モデルに対して、予め方程式記憶部３３が記憶する式（１４）から式（１５）を作成する。
【００９６】
そして、第一連立方程式作成部２２ａは、解析対象モデルが有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、第一方程式に対して、有限要素―境界積分法と、陰解法または陽解法とを適用した第一連立方程式を作成する。そして、第一連立方程式作成部２２ａは、作成した第一連立方程式を方程式記憶部３３に格納する。
【００９７】
例えば、第一連立方程式作成部２２ａは、式（１５）に対して、有限要素―境界積分法と陰解法とを適用して、式（１６）〜（１９）を作成し、作成した式（１６）〜（１９）を方程式記憶部３３に格納する。
【００９８】
第二連立方程式作成部２２ｂは、解析用モデル読み出し部２１ａから受信した解析対象モデルの物性情報に基づいて、解析対象モデルから発生する磁化ベクトルによる磁界である第二磁界を算出するための第二方程式を作成する。具体的には、第二連立方程式作成部２２ｂは、解析対象モデルの磁性体から発生する磁化ベクトルによる第二磁界を算出するために、解析用モデル読み出し部２１ａから受信した解析対象モデルに対して、予め方程式記憶部３３が記憶する式（１４）から式（２２）を作成する。
【００９９】
そして、第二連立方程式作成部２２ｂは、解析対象モデルが有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、第二方程式に対して、有限要素―境界積分法と、陰解法または陽解法とを適用した第二連立方程式を作成する。そして、第二連立方程式作成部２２ｂは、作成した第二連立方程式を方程式記憶部３３に格納する。
【０１００】
例えば、第二連立方程式作成部２２ｂは、式（２２）に対して、有限要素―境界積分法と陰解法とを適用して、式（２３）〜（２７）を作成し、作成した式（２３）〜（２７）を方程式記憶部３３に格納する。
【０１０１】
Ｂマトリックス計算部２３は、解析用モデル読み出し部２１ａおよび解析条件読み出し部２１ｂからそれぞれ解析対象モデルおよび解析条件を受信する。そして、Ｂマトリックス計算部２３は、受信した解析対象モデルおよび解析条件に基づいて、メッシュ分割した解析対象モデルの連続する各領域間の節点における変位マトリックス（以下、Ｂマトリックスと記す）を算出する。
【０１０２】
磁界計算部２４は、方程式記憶部が記憶する第一連立方程式および第二連立方程式を解法し、第一磁界計算部２４ａと、第二磁界計算部２４ｂと、合算磁界計算部２４ｃとを有する。
【０１０３】
第一磁界計算部２４ａは、方程式記憶部３３が記憶する第一連立方程式を解法することで、時系列に沿った第一磁界を計算し、計算結果を計算結果記憶部３４に格納する。例えば、第一磁界計算部２４ａは、方程式記憶部３３が記憶する式（１６）〜（１９）を順番に解法することにより、時系列に沿った第一磁界を計算し計算結果を計算結果記憶部３４に格納する。
【０１０４】
第二磁界計算部２４ｂは、方程式記憶部３３が記憶する第二連立方程式を解法することで、時系列に沿った第二磁界を計算し、計算結果を計算結果記憶部３４に格納する。例えば、第二磁界計算部２４ｂは、方程式記憶部３３が記憶する式（２３）〜（２７）を順番に解法することにより、時系列に沿った第二磁界を計算し、計算結果を計算結果記憶部３４に格納する。
【０１０５】
合算磁界計算部２４ｃは、計算結果記憶部３４が記憶する第一磁界の計算結果および第二磁界の計算結果を合算した解析結果を算出し、算出した解析結果を結果データ記憶部３５に格納すると共に、入出力制御Ｉ／Ｆ部６０を介して出力部５０に送信する。例えば、合算磁界計算部２４ｃは、予め方程式記憶部３３が記憶する式（３０）を用いて、第一磁界と第二磁界とを合算した解析結果を算出し、算出した解析結果を結果データ記憶部３５に格納すると共に、入出力制御Ｉ／Ｆ部６０を介して出力部５０に送信する。
【０１０６】
なお、方程式記憶部３３に格納された式（１９）、式（２６）および式（２７）については、格納された以降の他の解析においても、解析ごとに作成せずに、既に格納されている式が用いられる場合であってもよい。
【０１０７】
［実施例２における磁界解析装置による処理の手順］
次に、図６を用いて、実施例２における磁界解析装置による処理を説明する。図６は、実施例２における磁界解析装置による磁界解析処理の手順を説明するための図である。
【０１０８】
図６に示すように、まず、実施例２における磁界解析装置１０は、利用者によって解析開始コマンドが実行されると（ステップＳ１０１肯定）、解析用モデル読み出し部２１ａは、モデルデータ記憶部３１から解析対象モデルを読み出して、方程式作成部２２に送信する（ステップＳ１０２）。具体的には、利用者によって解析対象モデルの指定および解析条件の設定が実行されると、解析用モデル読み出し部２１ａが、モデルデータ記憶部３１から解析対象モデルを読み出す。そして、解析用モデル読み出し部２１ａは、読み出した解析対象モデルを方程式作成部２２に送信する。
【０１０９】
そして、解析条件読み出し部２１ｂは、解析用モデル読み出し部２１ａが読み出した解析対象モデルに対する解析条件データを、解析条件データ記憶部３２から読み出して、方程式作成部２２に送信する（ステップＳ１０３）。
【０１１０】
そして、第一連立方程式作成部２２ａおよび第二連立方程式作成部２２ｂは、それぞれ電流ベクトルによる磁界を計算するための第一連立方程式および磁化ベクトルによる磁界を計算するための第二連立方程式を作成し、方程式作成部３３に格納する（ステップＳ１０４）。
【０１１１】
そののち、Ｂマトリックス計算部２３は、解析用モデル読み出し部２１ａおよび解析条件読み出し部２１ｂからそれぞれ受信した解析対象モデルおよび解析条件に基づいて、Ｂマトリックスを算出する（ステップＳ１０５）。
【０１１２】
そして、第一磁界計算部２４ａおよび第二磁界計算部２４ｂそれぞれは、方程式記憶部３３が記憶する第一連立方程式および第二連立方程式に対して「時間：ｔ」を初期化（ｔ＝０）する（ステップＳ１０６）。
【０１１３】
そののち、第一磁界計算部２４ａは、利用者によって設定された解析条件に基づいて、求める時間ステップの電流を設定し（ステップＳ１０７）、現在の時間ステップから求める時間ステップにおける第一磁界の値を算出し、算出した第一磁界の値を計算結果記憶部３４に格納する（ステップＳ１０８）。
【０１１４】
なお、ステップＳ１０７と並行して、第二磁界計算部２４ｂは、利用者によって設定された解析条件に基づいて、求める時間ステップの磁化を設定する（ステップＳ１０９）。そして、ステップＳ１０８と並行して、第二磁界計算部２４ｂは、現在の時間ステップから求める時間ステップにおける第二磁界の値を算出し、算出した第二磁界の値を計算結果記憶部３４に格納する（ステップＳ１１０）。
【０１１５】
続いて、合算磁界計算部２４ｃは、計算結果記憶部３４が記憶する第一磁界の計算結果および第二磁界の計算結果を合算した解析結果を算出し、算出した解析結果を計算結果記憶部３４に格納する（ステップＳ１１１）。そして、合算時間計算部２４ｃは、利用者が設定したすべての時間ステップにおける解析が終了したか否かの判定を行う（ステップＳ１１２）。
【０１１６】
ここで、解析が終了していない場合（ステップＳ１１２否定）、第一磁界計算部２４ａおよび第二磁界計算部２４ｂは、それぞれステップＳ１０７およびステップＳ１０９に戻って、次の時間ステップの電流および磁化の設定を行う。
【０１１７】
一方、解析が終了した場合（ステップＳ１１２肯定）、合算磁界計算部２４ｃは、解析結果を結果データ記憶部３５に格納して（ステップＳ１１３）、処理を終了する。
【０１１８】
［実施例２の効果］
上記したように、実施例２によれば、入力部４０は、解析対象物の物性情報を受け付ける。そして、第一連立方程式作成部２２ａは、入力部４０から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する電流ベクトルによる第一磁界を算出するための第一の方程式を有限要素法および境界積分法により解析可能な形式で作成する。また、第二連立方程式作成部２２ｂは、入力部４０から受け付けた解析対象物の物性情報に基づいて、解析対象物から発生する磁化ベクトルによる第二磁界を算出するための第二の方程式を有限要素法および境界積分法により解析可能な形式で作成する。そして、第一磁界計算部２４ａは、第一連立方程式作成部２２ａによって作成された方程式を用いて、第一磁界を算出する。また、第二磁界計算部２４ｂは、第二連立方程式作成部２２ｂによって作成された方程式を用いて、第二磁界を算出する。従って、磁化ベクトルを含む方程式に対して有限要素法および境界積分法を適用することができ、磁化ベクトルの影響を正確に考慮することが可能となる。
【０１１９】
また、実施例２によれば、入力部４０は、解析対象物が有する物体それぞれの位置情報をさらに受け付ける。そして、第一連立方程式作成部２２ａは、入力部４０から受け付けた解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、第一の方程式に陰解法または陽解法を適用することで、第一磁界の時系列に沿った値を算出するための第一連立方程式をさらに作成する。そして、第一磁界計算部２４ａは、第一連立方程式作成部２２ａによって作成された第一連立方程式を用いて、第一磁界を算出する。従って、時系列に沿った解析を容易に行うことができ、非定常磁界を含めた第一磁界を、電流ベクトルによる影響を正確に考慮して簡易に算出することができる。
【０１２０】
また、実施例２によれば、入力部４０は、解析対象物が有する物体それぞれの位置情報をさらに受け付ける。そして、第二連立方程式作成部２２ｂは、入力部４０から受け付けた解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、第二の方程式に陰解法または陽解法を適用することで、第二磁界の時系列に沿った値を算出するための第二連立方程式をさらに作成する。そして、第二磁界計算部２４ｂは、第二連立方程式作成部２２ｂによって作成された第二連立方程式を用いて、第二磁界を算出する。従って、時系列に沿った解析を容易に行うことができ、非定常磁界を含めた第二磁界を、磁化ベクトルによる影響を正確に考慮して簡易に算出することができる。
【０１２１】
また、実施例２によれば、合算磁界計算部２４ｃは、第一磁界計算部２４ａおよび第二磁界計算部２４ｂによってそれぞれ算出された第一磁界および第二磁界を合算し、合算した合算磁界を解析対象物における磁界の解析結果と決定する。従って、電流ベクトルによる磁界と磁化ベクトルによる磁界とを考慮した磁界解析を行うことができる。
【０１２２】
以下、本実施例における数値解析法を用いた解析例の結果を説明する。本解析例は、図３に示す解析対象モデルに対して、図４の（Ａ）、（Ｂ）および（Ｃ）に示す解析条件で解析を行った結果である。なお、本解析例では、材料Ａ，Ｂともに移動しない条件で実施していることから、運動条件は、材料Ａ、Ｂいずれも初期座標Ｘ０＝（０、０、０）、移動速度Ｖ＝（０、０、０）である。
【０１２３】
図７に示すように、本解析例では、本実施例における数値解析法を用いた解析例である本手法に加え、比較対象として、有限要素法を用いた既知の磁場解析ソフトであるソフトＡおよびソフトＢによる上記解析対象モデルの解析を実施している。なお、図７は、解析例の結果を説明するための図である。また、ソフトＡおよびソフトＢは、いずれも有限要素法を用いた数値解析であり、空間領域をメッシュ分割する必要があるため、本実施例における数値解析法と比較して、解析結果の算出に要する時間が長時間必要となる数値解析法である。
【０１２４】
本解析例では、図７に示すように、材料ＡのＺ＝５μｍ、Ｚ＝１５μｍおよびＺ＝２５μｍ位置における磁束密度ベクトルのＺ成分（Ｂｚ）を各数値解析法により経時的に算出した結果を示している。なお、グラフの縦軸は、磁束密度ベクトルのＺ成分（Ｂｚ）を示し、横軸は、時間を示している。
【０１２５】
図７に示すように、本手法による解析結果は、ソフトＡによる解析結果およびソフトＢによる解析結果と比較して、いずれの位置においても、ほぼ同様の結果が得られている。従って、上記した解析例の結果から、本実施例における数値解析法は、精度の高い解析を短時間で実行することが可能であることが示された。
【実施例３】
【０１２６】
さて、これまで実施例１および実施例２について説明したが、上述した実施例１および実施例２以外にも、種々の異なる形態にて実施されてよいものである。そこで、以下では、種々の異なる実施例を（１）〜（３）に区分けして説明する。
【０１２７】
（１）解析対象モデルおよび解析条件
上記した実施例１および実施例２では、解析対象モデルの指定と解析条件の設定とを解析開始時に同時に実行する場合について説明したが、本実施例はこれに限定されるものではなく、例えば、解析条件の設定を解析中に適宜変更する場合であってもよい。
【０１２８】
（２）システム構成等
また、上記の実施例において説明した処理手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に更新することができる。例えば、本実施例は、モデルデータ記憶部３１が記憶する解析対象モデルおよび解析条件データ記憶部３２が記憶する解析条件データの更新および新規追加を任意に行うことができる。
【０１２９】
また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示された構成要素と同一であることを要しない。すなわち、各処理部および各記憶部の分散・統合の具体的形態（例えば、図２の形態）は図示のものとは限られない。例えば、第一連立方程式作成部２２ａと第二連立方程式作成部２２ｂとが統合されている場合であってもよい。さらに、各装置にて行なわれる各処理機能は、その全部または任意の一部が、ＣＰＵおよび当該ＣＰＵにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。
【０１３０】
（３）磁界解析プログラム
ところで、上記の実施例２では、ハードウェアロジックによって各種の処理を実現する場合を説明したが、本実施例はこれに限定されるものではなく、予め用意されたプログラムをコンピュータで実行するようにしてもよい。そこで、以下では、図８を用いて上記の実施例２に示した磁界解析装置１０と同様の機能を有する磁界解析プログラムを実行するコンピュータの一例を説明する。図８は、実施例２における磁界解析プログラムを実行するコンピュータを説明するための図である。
【０１３１】
図８に示すように、情報処理装置としてのコンピュータ１０００は、キーボード１０２０、モニタ１０３０、ＲＡＭ１０４０、ＨＤＤ１０５０、ＣＰＵ１０６０およびＲＯＭ１０７０を有する。そして、キーボード１０２０、モニタ１０３０、ＲＡＭ１０４０、ＨＤＤ１０５０、ＣＰＵ１０６０およびＲＯＭ１０７０は、バス１０１０などで接続される。なお、図示していないが、コンピュータ１０００は、ネットワークを介して、デコーダと接続されている。
【０１３２】
ＲＯＭ１０７０には、上記の実施例２に示した磁界解析装置１０と同様の機能を発揮する磁界解析プログラム、つまり、図８に示すように、解析用モデル読み出しプログラム１０７１、解析条件読み出しプログラム１０７２、第一連立方程式作成プログラム１０７３、第二連立方程式作成プログラム１０７４、Ｂマトリックス計算プログラム１０７５、第一磁界計算プログラム１０７６、第二磁界計算プログラム１０７７、合算磁界計算プログラム１０７８が予め記憶されている。なお、これらのプログラム１０７１〜１０７８については、図２に示した磁界解析装置１０の各構成要素と同様、適宜統合または分散してもよい。
【０１３３】
そして、ＣＰＵ１０６０が、これらのプログラム１０７１〜１０７８をＲＯＭ１０７０から読み出して実行することで、図８に示すように、各プログラム１０７１〜１０７８は、解析用モデル読み出しプロセス１０６１、解析条件読み出しプロセス１０６２、第一連立方程式作成プロセス１０６３、第二連立方程式作成プロセス１０６４、Ｂマトリックス計算プロセス１０６５、第一磁界計算プロセス１０６６、第二磁界計算プロセス１０６７、合算磁界計算プロセス１０６８として機能するようになる。なお、各プロセス１０６１〜１０６８は、図２に示した、解析用モデル読み出し部２１ａ、解析条件読み出し部２１ｂ、第一連立方程式作成部２２ａ、第二連立方程式作成部２２ｂ、Ｂマトリックス計算部２３、第一磁界計算部２４ａ、第二磁界計算部２４ｂ、合算磁界計算部２４ｃにそれぞれ対応する。
【０１３４】
また、ＨＤＤ１０５０には、図８に示すように、モデルデータ記憶データ１０５１、解析条件データ記憶データ１０５２、方程式記憶データ１０５３、計算結果記憶データ１０５４、結果データ記憶データ１０５５が設けられる。このモデルデータ記憶データ１０５１、解析条件データ記憶データ１０５２、方程式記憶データ１０５３、計算結果記憶データ１０５４、結果データ記憶データ１０５５は、それぞれ図２に用いたモデルデータ記憶部３１、解析条件データ記憶部３２、方程式記憶部３３、計算結果記憶部３４、結果データ記憶部３５に対応する。そして、ＣＰＵ１０６０は、モデルデータ記憶データ１０４１、解析条件データ記憶データ１０４２、方程式記憶データ１０４３、計算結果記憶データ１０４４、結果データ記憶データ１０４５をそれぞれモデルデータ記憶データ１０５１、解析条件データ記憶データ１０５２、方程式記憶データ１０５３、計算結果記憶データ１０５４、結果データ記憶データ１０５５に対して登録し、このモデルデータ記憶データ１０４１、解析条件データ記憶データ１０４２、方程式記憶データ１０４３、計算結果記憶データ１０４４、結果データ記憶データ１０４５を読み出してＲＡＭ１０４０に格納する。そして、ＣＰＵ１０６０は、ＲＡＭ１０４０に格納されたモデルデータ記憶データ１０４１、解析条件データ記憶データ１０４２、方程式記憶データ１０４３、計算結果記憶データ１０４４、結果データ記憶データ１０４５に基づいて磁界解析処理を実行する。
【０１３５】
なお、上記した各プログラム１０７１〜１０７８については、必ずしも最初からＲＯＭ１０７０に記憶させておく必要は無く、例えばコンピュータ１０００に挿入されるフレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯディスク、ＤＶＤディスク、光磁気ディスク、ＩＣカードなどの「可搬用の物理媒体」、または、コンピュータ１０００の内外に備えられるＨＤＤなどの「固定用物理媒体」、さらには、公衆回線、インターネット、ＬＡＮ、ＷＡＮなどを介してコンピュータ１０００に接続される「他のコンピュータ（またはサーバ）」などに各プログラムを記憶させておき、コンピュータ１０００がこれらから各プログラムを読み出して実行するようにしてもよい。
【０１３６】
以上の本実施例を含む実施形態に関し、更に以下の付記を開示する。
【０１３７】
（付記１）解析対象物の物性情報を受け付ける入力部と、
前記入力部から受け付けた前記解析対象物の物性情報に基づいて、当該解析対象物から発生する電流ベクトルによる第一の磁界を算出するための第一の方程式と、当該解析対象物から発生する磁化ベクトルによる第二の磁界を算出するための第二の方程式とを、有限要素法および境界積分法により解析可能な形式で作成する方程式作成部と、
前記方程式作成部によって作成された方程式を用いて、前記第一の磁界および前記第二の磁界それぞれを算出する磁界計算部と、
を有したことを特徴とする磁界解析装置。
【０１３８】
（付記２）前記入力部は、前記解析対象物が有する物体それぞれの位置情報をさらに受け付け、
前記方程式作成部は、前記入力部から受け付けた前記解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、前記第一の方程式に陰解法または陽解法を適用することで、前記第一の磁界の時系列に沿った値を算出するための複数の方程式である第一の連立方程式をさらに作成し、
前記磁界計算部は、前記方程式作成部によって作成された前記第一の連立方程式を用いて、前記第一の磁界を算出することを特徴とする付記１に記載の磁界解析装置。
【０１３９】
（付記３）前記入力部は、前記解析対象物が有する物体それぞれの位置情報をさらに受け付け、
前記方程式作成部は、前記入力部から受け付けた前記解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、前記第二の方程式に陰解法または陽解法を適用することで、前記第二の磁界の時系列に沿った値を算出するための複数の方程式である第二の連立方程式をさらに作成し、
前記磁界計算部は、前記方程式作成部によって作成された前記第二の連立方程式を用いて、前記第二の磁界を算出することを特徴とする付記１に記載の磁界解析装置。
【０１４０】
（付記４）前記磁界計算部は、算出した前記第一の磁界および前記第二の磁界を合算し、合算した合算磁界を前記解析対象物における磁界の解析結果と決定することを特徴とする付記２または３に記載の磁界解析装置。
【０１４１】
（付記５）解析対象物の物性情報を受け付ける入力手順と、
前記入力手順から受け付けた前記解析対象物の物性情報に基づいて、当該解析対象物から発生する電流ベクトルによる第一の磁界を算出するための第一の方程式と、当該解析対象物から発生する磁化ベクトルによる第二の磁界を算出するための第二の方程式とを、有限要素法および境界積分法により解析可能な形式で作成する方程式作成手順と、
前記方程式作成部によって作成された方程式を用いて、前記第一の磁界および前記第二の磁界それぞれを算出する磁界計算手順と、
をコンピュータに実行させることを特徴とする磁界解析プログラム。
【０１４２】
（付記６）前記入力手順は、前記解析対象物が有する物体それぞれの位置情報をさらに受け付け、
前記方程式作成手順は、前記入力部から受け付けた前記解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、前記第一の方程式に陰解法または陽解法を適用することで、前記第一の磁界の時系列に沿った値を算出するための複数の方程式である第一の連立方程式をさらに作成し、
前記磁界計算手順は、前記方程式作成部によって作成された前記第一の連立方程式を用いて、前記第一の磁界を算出することを特徴とする付記５に記載の磁界解析プログラム。
【０１４３】
（付記７）前記入力部手順は、前記解析対象物が有する物体それぞれの位置情報をさらに受け付け、
前記方程式作成手順は、前記入力部から受け付けた前記解析対象物が有する物体それぞれの位置情報の時系列に沿った相対的位置情報の変化に基づいて、前記第二の方程式に陰解法または陽解法を適用することで、前記第二の磁界の時系列に沿った値を算出するための複数の方程式である第二の連立方程式をさらに作成し、
前記磁界計算手順は、前記方程式作成部によって作成された前記第二の連立方程式を用いて、前記第二の磁界を算出することを特徴とする付記５に記載の磁界解析プログラム。
【０１４４】
（付記８）前記磁界計算手順は、算出した前記第一の磁界および前記第二の磁界を合算し、合算した合算磁界を前記解析対象物における磁界の解析結果と決定することを特徴とする付記６または７に記載の磁界解析プログラム。
【符号の説明】
【０１４５】
１０  磁界解析装置
２０  処理部
２１  入力情報処理部
２１ａ  解析用モデル読み出し部
２１ｂ  解析条件読み出し部
２２  方程式作成部
２２ａ  第一連立方程式作成部
２２ｂ  第二連立方程式作成部
２３  Ｂマトリックス計算部
２４  磁界計算部
２４ａ  第一磁界計算部
２４ｂ  第二磁界計算部
２４ｃ  合算磁界計算部
３０  記憶部
３１  モデルデータ記憶部
３２  解析条件データ記憶部
３３  方程式記憶部
３４  計算結果記憶部
３５  結果データ記憶部
４０  入力部
５０  出力部
６０  入出力制御Ｉ／Ｆ部
【図１】
0005412982.tif 000032
【図２】
0005412982.tif 000033
【図３】
0005412982.tif 000034
【図４】
0005412982.tif 000035
【図５】
0005412982.tif 000036
【図６】
0005412982.tif 000037
【図７】
0005412982.tif 000038
【図８】
0005412982.tif 000039