JP2004219178A

JP2004219178A - 磁区構造解析装置、磁区構造解析方法、磁区構造解析プログラム及び記録媒体

Info

Publication number: JP2004219178A
Application number: JP2003005012A
Authority: JP
Inventors: Keiji Iwata; 圭司岩田; Masahiro Fujikura; 昌浩藤倉; Muneyuki Imafuku; 宗行今福
Original assignee: Nippon Steel Corp
Current assignee: Nippon Steel Corp
Priority date: 2003-01-10
Filing date: 2003-01-10
Publication date: 2004-08-05
Anticipated expiration: 2023-01-10
Also published as: JP4057922B2

Abstract

【課題】容易に精度が高い磁区構造のシミュレーションを行うことができる磁区構造解析装置、磁区構造解析方法、磁区構造解析プログラム及び記録媒体を提供する。
【解決手段】コンピュータ等の計算機を用いて、電磁鋼板の静磁エネルギＥｓｔ、磁区内の結晶磁気異方性エネルギＥａｎ及び磁壁エネルギＥｗを求め、これらの総和（総エネルギ）を求める。磁壁エネルギＥｗは磁壁内の結晶磁気異方性エネルギ及び交換エネルギから構成される。そして、総エネルギが最小となるような磁化の方向、磁化が生成するベクトルポテンシャル及び磁区幅を、有限要素法又は差分法等により求める。この有限要素法に当たっては、磁区の大きさを変化させる際に、磁区幅が変化する方向のスケールを磁区幅の変化に応じて変化させる。このとき、全領域の再区画を行うのではなく、磁区幅が変化する方向のスケールを変化させる。
【選択図】図３

Description

【０００１】
【発明の属する技術分野】
本発明は、電磁鋼板等の磁性材料の開発に好適な磁性材料の磁区構造解析装置、磁区構造解析方法、磁区構造解析プログラム及び記録媒体に関する。
【０００２】
【従来の技術】
磁性材料の特性、例えばエネルギ損失及び磁歪等は、一般的に、磁区構造に起因すると考えられている。従って、磁性材料の特性向上は、磁区構造を最適化することにより達成される。磁区構造を変化させる要因は多数存在し、例えば材料自体の形状、表面の状態、結晶方位及び歪等が挙げられる。
【０００３】
そして、磁区構造を変化させることを目的として、トランス等に使用される電磁鋼板等では、その表面に溝等の外乱が付与されている。このような外乱の付与により、磁区が細分化される。一方、磁性材料の損失には、大別して、ヒステリシス損と渦電流損とが存在する。渦電流損は、更に古典渦電流損と異常渦電流損とに分類される。上記の磁区の細分化には、異常渦電流損を低減する効果がある。
【０００４】
磁区構造を解析する方法としては、例えば、コロイド法、磁区ＳＥＭ、カー効果又はファラデー効果を利用した磁気光学法等がある。但し、これらの方法では、使用の変化に応じて、その都度試料を作る必要があり、多大な手間及び時間が必要とされる。また、試料の表面近傍の磁区構造の解析は可能であるが、試料の内部の磁区構造を解析することが困難である。これは、磁性材料の磁気特性には、形状異方性があるため、試料の切断等を行うと、本来の磁区構造が変化してしまうためである。
【０００５】
そこで、計算機を使用した磁区構造のシミュレーション方法が研究されている。例えば、文献（ＲｉｃｃａｒｄｏＨｅｒｔｅｌａｎｄＨｅｌｍｕｔＫｒｏｎｍｕｌｌｅｒ、「ＰｈｙｓｉｃａｌＲｅｖｉｅｗＢ」、１９９９年、６０巻、ｐ．７３６６−７３７８）には、パーマロイ薄膜における磁区構造の解析を、Ｌａｎｄａｕ−Ｌｉｆｓｈｉｔｚ−Ｇｉｌｂｅｒｔ（ＬＬＧ）方程式を用いて行う方法が記載されている。この解析で対象とされているパーマロイ薄膜は、磁気記録装置の磁気ヘッド等に使用されており、そのサイズはナノスケールである。このため、磁区の大きさと磁壁の大きさとの差は小さく、数値計算における最小スケールで磁壁内の磁化分布を表現することが可能である。従って、パーマロイ薄膜については、ＬＬＧ方程式を直接適用してその磁区構造をシミュレーションすることが可能である。
【０００６】
また、方向性珪素鋼板についての磁区細分化効果の解析方法が文献（石田昌義、中野恒、本田厚人、佐藤圭司、「方向性珪素鋼板に導入した線状溝による磁区細分化効果の解析」、日本応用磁気学会誌、１９９４年、１８巻、ｐ．８０９−８１３）に記載されている。一般に、電磁鋼板については、上述のように、その表面に溝等の外乱を導入して磁区細分化を行い、異常渦電流を下げる開発が行われている。そして、上記の文献には、溝による磁区細分化の効果を極めて単純な計算モデルを用いて調べる方法が記載されている。
【０００７】
【非特許文献１】
ＲｉｃｃａｒｄｏＨｅｒｔｅｌａｎｄＨｅｌｍｕｔＫｒｏｎｍｕｌｌｅｒ、「ＰｈｙｓｉｃａｌＲｅｖｉｅｗＢ」、１９９９年、６０巻、ｐ．７３６６−７３７８
【非特許文献２】
石田昌義、中野恒、本田厚人、佐藤圭司、「方向性珪素鋼板に導入した線状溝による磁区細分化効果の解析」、日本応用磁気学会誌、１９９４年、１８巻、ｐ．８０９−８１３
【０００８】
【発明が解決しようとする課題】
しかしながら、アモルファス材料を含む電磁鋼板のような軟質磁性材料においては、磁区の大きさが数百μｍであるのに対し、磁壁の厚さが数十ｎｍ程度であり、これらの大きさの相違が極めて大きい。従って、パーマロイ薄膜と同様にして、ＬＬＧ方程式を直接適用すると、分割数、即ち未知変数の数が増大し、膨大な計算量が必要とされ、非現実的である。
【０００９】
また、上記文献に記載された極めて単純な計算モデルを用いた方法は、解析的な計算に基づくものであるため、矩形の溝が形成されているモデル以外に適用することができない。更に、磁化の変化や結晶磁気異方性も考慮されていないため、還流磁区を扱うことができず、磁気エネルギを正確に見積もることができていない。この結果、このシミュレーションでも、実測値と計算値との差が大きく、この方法に基づいて磁区構造を最適化することは極めて困難である。
【００１０】
本発明は、かかる問題点に鑑みてなされたものであり、容易に精度が高い磁区構造のシミュレーションを行うことができる磁区構造解析装置、磁区構造解析方法、磁区構造解析プログラム及び記録媒体を提供することを目的とする。
【００１１】
【課題を解決するための手段】
本願発明者は、鋭意検討の結果、以下に示す発明の諸態様に想到した。
【００１２】
本願の第１の発明に係る磁区構造解析装置は、磁性材料内の磁区の構造を段階的に変化させ、段階毎に前記磁性材料の磁気エネルギを算出して、最も磁気エネルギが小さくなる段階を求める磁区構造解析装置であって、磁化を考慮しながら磁区内の静磁エネルギを算出する静磁エネルギ算出手段と、磁化を考慮しながら磁区内の結晶磁気異方性エネルギを算出する結晶磁気異方性エネルギ算出手段と、磁壁内の磁壁エネルギを算出する磁壁エネルギ算出手段と、前記静磁エネルギ、前記結晶磁気異方性エネルギ及び前記磁壁エネルギの総和を前記磁気エネルギとして算出する総エネルギ算出手段と、を有することを特徴とする。
【００１３】
本願の第２の発明に係る磁区構造解析方法は、計算機を用いて、磁性材料内の磁区の構造を段階的に変化させ、段階毎に前記磁性材料の磁気エネルギを算出して、最も磁気エネルギが小さくなる段階を求める磁区構造解析方法であって、磁化を考慮しながら磁区内の静磁エネルギを算出する静磁エネルギ算出工程と、磁化を考慮しながら磁区内の結晶磁気異方性エネルギを算出する結晶磁気異方性エネルギ算出工程と、磁壁内の磁壁エネルギを算出する磁壁エネルギ算出工程と、前記静磁エネルギ、前記結晶磁気異方性エネルギ及び前記磁壁エネルギの総和を前記磁気エネルギとして算出する総エネルギ算出工程と、を有することを特徴とする。
【００１４】
本願の第３の発明に係る磁区構造解析プログラムは、コンピュータに、磁性材料内の磁区の構造を段階的に変化させ、段階毎に前記磁性材料の磁気エネルギを算出して、最も磁気エネルギが小さくなる段階を求める処理を実行させる磁区構造解析プログラムであって、前記コンピュータに、磁化を考慮しながら磁区内の静磁エネルギを算出する静磁エネルギ算出処理と、磁化を考慮しながら磁区内の結晶磁気異方性エネルギを算出する結晶磁気異方性エネルギ算出処理と、磁壁内の磁壁エネルギを算出する磁壁エネルギ算出処理と、前記静磁エネルギ、前記結晶磁気異方性エネルギ及び前記磁壁エネルギの総和を前記磁気エネルギとして算出する総エネルギ算出処理と、を実行させることを特徴とする。
【００１５】
本願の第４の発明に係るコンピュータ読み取り可能な記録媒体は、上記の磁区構造解析プログラムを記録したことを特徴とする。
【００１６】
【発明の実施の形態】
以下、本発明の実施形態に係る磁性材料の磁区構造解析装置、磁区構造解析方法、磁区構造解析プログラム及びコンピュータ読み取り可能な記録媒体について、添付の図面を参照して詳細に説明する。
【００１７】
先ず、磁性材料の磁区構造の解析方法について説明する。本実施形態においては、コンピュータ等の計算機を用いて、磁性材料、例えば電磁鋼板の静磁エネルギＥｓｔ、磁区内の結晶磁気異方性エネルギＥａｎ及び磁壁エネルギＥｗを求め、これらの総和（総エネルギ）を求める。ここで、磁壁エネルギＥｗは、例えば磁壁内の結晶磁気異方性エネルギ及び交換エネルギから構成される。また、この磁壁エネルギＥｗとして、例えば１８０度磁区を区画する磁壁の磁壁エネルギを考慮する。そして、総エネルギが最小となるような磁化の方向、磁化が生成するベクトルポテンシャル及び磁区幅を、有限要素法又は差分法等により求める。
【００１８】
ここで、静磁エネルギＥｓｔ、磁区内の結晶磁気異方性エネルギＥａｎ及び磁壁エネルギＥｗを求める方法について説明する。
【００１９】
先ず、静磁エネルギＥｓｔを求める方法について説明する。静磁エネルギＥｓｔは、例えば有限要素法又は差分法等により求めることができ、本実施形態では、有限要素法を採用する。具体的には、先ず、図１に示すように、シミュレーションの対象となる電磁鋼板１が含まれる空間を、任意の形状の複数の微小な要素（領域）２に区画する。次に、各要素２について、ベクトルポテンシャル、磁化方位及び結晶方位を設定する。このとき、電磁鋼板１の周囲の空気に属する要素２には、ベクトルポテンシャルのみを設定する。
【００２０】
磁化は、｜ｍ（ベクトル）｜＝Ｍｓ（自発磁化の大きさ）の条件の下で、任意の方向を向くことができるように設定する。例えば、要素２の集合に含まれる要素ｉの磁化は、例えば下記数式１のように極座標表示で表す。但し、θ_ｉ及びφ_ｉは、未知変数である。
【００２１】
【数１】

【００２２】
そして、図２に示すモデルに基づいて、静磁エネルギＥｓｔを下記数式２のような汎関数で表現する。数式２以降の式中の太字は、その文字が表すものが方向をもっていることを示す。
【００２３】
【数２】

【００２４】
ここで、数式２中の太字の「Ａ」は、ベクトルポテンシャルであり、下記数式３のように表すことができる。
【００２５】
【数３】

【００２６】
ここで、Ａ_ｋは辺上の未知変数であり、このベクトルポテンシャルでは、要素の辺上に射影した値を未知変数とする補間関数Ｎ_ｋ（ベクトル）が利用されている。この関数は、特許第２８１９２２０号公報及び文献（岩田圭司、「辺要素法三次元磁場解析における電流連続入力方法」、日本応用磁気学会誌、１９９４年、１８巻、ｐ．６０７−６１２）に記載されているように、マクロな電磁場解析において高速計算が可能であると知られているが、磁性材料の磁化等を扱うマイクロマグネティクスの分野への適用は、これまで行われていない。
【００２７】
なお、ベクトルポテンシャルの表し方は、上記の方法に限定されるものではなく、例えば未知変数を要素の節点上において補間関数を定めてもよい。また、補間関数も数式３で表されるものに限定されるものではない。例えば、数式３で表される補間関数は１次要素用（ｋ＝１〜１２）であるが、高次要素用の補間関数を用いてもよい。
【００２８】
次に、磁区内の結晶磁気異方性エネルギＥａｎを求める方法について説明する。結晶磁気異方性は、各要素に結晶方位を設定することで定義することができる。要素内の磁化は、静磁エネルギが存在しない場合には、磁化容易軸と一致したときに最も安定するが、実際には、静磁エネルギの存在により、ある特定の方向で落ち着くこととなる。例えば、電磁鋼板（９７質量％Ｆｅ＋３質量％Ｓｉ）には、磁化容易軸が３方向存在し、結晶構造が立方晶であれば、この結晶磁気異方性エネルギＥａｎは、下記数式４で表すことができる。
【００２９】
【数４】

【００３０】
ここで、α_１、α_２及びα_３は、磁化ｍ（ベクトル）の磁化容易軸に対する方向余弦であり、Ｋは磁気異方性定数である。
【００３１】
次に、磁壁エネルギＥｗを求める方法について説明する。この磁壁エネルギＥｗには、前述のように磁壁内の結晶磁気異方性エネルギ及び交換エネルギが含まれている。そして、磁気異方性定数Ｋ及び交換定数Ａから、π（２ＡＫ）^１／２で表される値をε_ｗとし、磁区の幅をｄすると、磁壁エネルギＥｗは、下記数式５で表すことができる。
【００３２】
【数５】

【００３３】
このように、総エネルギは、汎関数Ｅ（Ａｋ、θ_ｉ、φ_ｉ、・・・）で表され、ベクトルポテンシャルの未知変数Ａｋ並びに磁化の向きθ_ｉ及びφ_ｉ等から構成される。
【００３４】
次に、総エネルギが最小となる磁化方向、磁化が生成するベクトルポテンシャル及び磁区幅を有限要素法により求める方法について説明する。
【００３５】
本実施形態では、磁区の大きさ及び要素２の大きさを変化させながら、上記の磁化方向等を求める。このとき、磁区の大きさを変化させる際に、磁区幅が変化する方向のスケールを磁区幅の変化に応じて変化させる。本実施形態では、静磁エネルギＥｓｔを用いており、初期状態では、図３（ａ）に示すように、互いに磁化５の向きが相反する２つの磁区３間に１８０度磁壁４が存在し、磁区３内に複数の要素２が存在している。そして、磁化方向等を求めるに当たって、磁区３の大きさを変化させるが、このとき、本実施形態では、図３（ｂ）に示すように、全領域の再区画を行うのではなく、図３（ｃ）に示すように、磁区幅が変化する方向のスケールを変化させる。
【００３６】
そして、本実施形態では、各未知変数（Ａｋ、θ_ｉ、φ_ｉ）に対する勾配、及びこの勾配の勾配を用いた最適化法、例えばＮｅｗｔｏｎ−Ｒａｐｈｓｏｎ法等により、総エネルギの最小値を求める。Ｎｅｗｔｏｎ−Ｒａｐｈｓｏｎ法に基づく計算を行うためには、総エネルギの各未知変数に対する勾配、及びこの勾配の勾配から構成される多次元連立方程式を解く必要があるが、本実施形態では、この計算を収束計算により行う。
【００３７】
このような方法で総エネルギの最小値を求める場合、要素分割を行う回数は１回のみであるため、極めて効率よく計算を行うことが可能となる。また、有限要素法を用いることにより、複雑な形状の磁性材料にも容易に対応することが可能である。つまり、電磁鋼板では、前述のように、磁区の細分化のために、鋼板に溝を形成することが多いが、本実施形態によれば、いかなる形状の溝にも容易に対応することが可能である。
【００３８】
更に、本実施形態では、ベクトルポテンシャルの補間関数に要素の辺上に未知変数が設定される関数を用いているので、多次元連立方程式を解くために使用する行列が対角優勢の疎行列となり、計算が容易である。
【００３９】
次に、実際に、本願発明者が行ったシミュレーション及び実測の結果について説明する。この測定では、電磁鋼板として珪素鋼板（Ｆｅ−３質量％Ｓｉ）を対象とし、導入溝の深さと最適な磁区幅との関係を解析した。試料の形状については、珪素鋼板の厚さは０．２３ｍｍとし、溝の幅は０．１ｍｍとし、溝のピッチは５ｍｍとした。また、シミュレーションによる解析の条件として、鋼板と鋼板を取り囲む空間を含んだ解析領域全体の分割要素数は３５４００とした。また、鋼板部分に関しては、図１に示すように、周期境界条件及び鏡面対称条件を設定した。
【００４０】
また、実測の試料として、鋼板の表面の面方位と容易磁化軸［１００］とが直交しておらず、図４に示すように、容易磁化軸が１．５度だけ珪素鋼板の表面から傾斜しているものを使用した。シミュレーションでもこの傾斜を考慮した。また、シミュレーションでは、珪素鋼板には容易磁化軸が３軸存在し、その結晶は立方晶であるため、結晶磁気異方性を求めるに当たり、数式４を用いた。ここで、磁気異方性定数Ｋは、４．８×１０^４Ｊ／ｍ^３とした。更に、自発磁化Ｍｓは２．１６Ｔとし、１８０度磁壁の単位面積当たりのエネルギε_Ｗは３．９×１０^−３Ｊ／ｍ^２とした。
【００４１】
この結果を、図５に示す。図５中の●は本発明の上記の実施形態に係るシミュレーションの結果を示し、四角印は従来の方法に基づくシミュレーションの結果を示し、×は実測値を示す。図５に示すように、本実施形態によれば、従来のものと比較して、実測値に極めて近い結果が得られた。従来の方法は、磁化ｍ（ベクトル）が動かない固定された方法であるため、静磁エネルギが過大に見積もられ、このエネルギを下げるために、結果として磁区幅が小さくなっている。一方、本実施形態では、磁化の方向を未知変数とすることができているので、磁化の変化によるエネルギの利得を考慮することができ、静磁エネルギ等の磁気エネルギを高い精度で見積もることができる。この結果として、還流磁区構造を再現することが可能となり、実測値に近い計算結果が得られている。
【００４２】
このような本発明は、特に、電磁鋼板に好適である。これは、電磁鋼板では、磁区の大きさが数百μｍであるのに対し、磁壁の厚さが数十ｎｍ程度であるため、これらのスケール差が１万倍程度にまで及び、従来のパーマロイ薄膜等に適用しているシミュレーション方法では、空間を区画して方程式を解くに当たり、空間分割数に比例して未知変数が増大し、計算時間が膨大なものとなるためである。
【００４３】
なお、本発明の実施形態は、上述のように、コンピュータがプログラムを実行することによって実現することができる。また、プログラムをコンピュータに供給するための手段、例えばかかるプログラムを記録したＣＤ−ＲＯＭ等のコンピュータ読み取り可能な記録媒体又はかかるプログラムを伝送するインターネット等の伝送媒体も本発明の実施形態として適用することができる。また、上記のプログラムも本発明の実施形態として適用することができる。上記のプログラム、記録媒体、伝送媒体及びプログラムプロダクトは、本発明の範疇に含まれる。
【００４４】
【発明の効果】
以上詳述したように、本発明によれば、珪素鋼板及びアモルファス等の電磁鋼板のように、磁区の大きさ（幅）と磁壁の幅との間に大きなスケール差が存在する磁性材料に対しても、高い精度で速やかにその磁区の構造の解析を行うことができる。このため、試行錯誤に基づく実験を行わずとも、シミュレーションにより好ましい磁区構造を見つけ出し、設計及び仕様の変更を容易に行うことができる。
【図面の簡単な説明】
【図１】空間を区画する方法を示す模式図である。
【図２】静磁エネルギを求めるためのモデルを示す模式図である。
【図３】本発明の実施形態における、磁区の大きさを変化させた場合の要素の変化を示す模式図である。
【図４】試料とその容易磁化軸との関係を示す模式図である。
【図５】シミュレーション及び実測の結果を示すグラフである。
【符号の説明】
１；電磁鋼板
２；要素
３；磁区
４；磁壁
５；磁化

Claims

磁性材料内の磁区の構造を段階的に変化させ、段階毎に前記磁性材料の磁気エネルギを算出して、最も磁気エネルギが小さくなる段階を求める磁区構造解析装置であって、
磁化を考慮しながら磁区内の静磁エネルギを算出する静磁エネルギ算出手段と、
磁化を考慮しながら磁区内の結晶磁気異方性エネルギを算出する結晶磁気異方性エネルギ算出手段と、
磁壁内の磁壁エネルギを算出する磁壁エネルギ算出手段と、
前記静磁エネルギ、前記結晶磁気異方性エネルギ及び前記磁壁エネルギの総和を前記磁気エネルギとして算出する総エネルギ算出手段と、
を有することを特徴とする磁区構造解析装置。
前記磁壁エネルギは、磁壁内の結晶磁気異方性エネルギ及び交換エネルギからなることを特徴とする請求項１に記載の磁区構造解析装置。
前記各磁壁は、１８０°磁壁であることを特徴とする請求項１又は２に記載の磁区構造解析装置。
前記静磁エネルギ算出手段は、有限要素法又は差分法を用いて前記静磁エネルギを算出することを特徴とする請求項１乃至３のいずれか１項に記載の磁区構造解析装置。
前記静磁エネルギ算出手段は、前記磁性材料内の磁区の構造を段階的に変化させるにあたり、
前記磁性材料内の磁区の構造を最初に規定する際に前記磁性材料を複数の領域に区画し、
その後の前記磁区の構造の段階的な変化の度に、前記複数の領域の大きさを変化させながら前記静磁エネルギを算出することを特徴とする請求項１乃至４のいずれか１項に記載の磁区構造解析装置。
前記磁性材料は、電磁鋼板であることを特徴とする請求項１乃至５のいずれか１項に記載の磁区構造解析装置。
計算機を用いて、磁性材料内の磁区の構造を段階的に変化させ、段階毎に前記磁性材料の磁気エネルギを算出して、最も磁気エネルギが小さくなる段階を求める磁区構造解析方法であって、
磁化を考慮しながら磁区内の静磁エネルギを算出する静磁エネルギ算出工程と、
磁化を考慮しながら磁区内の結晶磁気異方性エネルギを算出する結晶磁気異方性エネルギ算出工程と、
磁壁内の磁壁エネルギを算出する磁壁エネルギ算出工程と、
前記静磁エネルギ、前記結晶磁気異方性エネルギ及び前記磁壁エネルギの総和を前記磁気エネルギとして算出する総エネルギ算出工程と、
を有することを特徴とする磁区構造解析方法。
前記磁壁エネルギは、磁壁内の結晶磁気異方性エネルギ及び交換エネルギからなることを特徴とする請求項７に記載の磁区構造解析方法。
前記各磁壁は、１８０°磁壁であることを特徴とする請求項７又は８に記載の磁区構造解析方法。
前記静磁エネルギ算出工程において、有限要素法又は差分法を用いて前記静磁エネルギを算出することを特徴とする請求項７乃至９のいずれか１項に記載の磁区構造解析方法。
前記静磁エネルギ算出工程において、前記磁性材料内の磁区の構造を段階的に変化させるにあたり、
前記磁性材料内の磁区の構造を最初に規定する際に前記磁性材料を複数の領域に区画し、
その後の前記磁区の構造の段階的な変化の度に、前記複数の領域の大きさを変化させながら前記静磁エネルギを算出することを特徴とする請求項７乃至１０のいずれか１項に記載の磁区構造解析方法。
前記磁性材料は、電磁鋼板であることを特徴とする請求項７乃至１１のいずれか１項に記載の磁区構造解析方法。
コンピュータに、磁性材料内の磁区の構造を段階的に変化させ、段階毎に前記磁性材料の磁気エネルギを算出して、最も磁気エネルギが小さくなる段階を求める処理を実行させる磁区構造解析プログラムであって、
前記コンピュータに、
磁化を考慮しながら磁区内の静磁エネルギを算出する静磁エネルギ算出処理と、
磁化を考慮しながら磁区内の結晶磁気異方性エネルギを算出する結晶磁気異方性エネルギ算出処理と、
磁壁内の磁壁エネルギを算出する磁壁エネルギ算出処理と、
前記静磁エネルギ、前記結晶磁気異方性エネルギ及び前記磁壁エネルギの総和を前記磁気エネルギとして算出する総エネルギ算出処理と、
を実行させることを特徴とする磁区構造解析プログラム。
前記磁壁エネルギは、磁壁内の結晶磁気異方性エネルギ及び交換エネルギからなることを特徴とする請求項１３に記載の磁区構造解析プログラム。
前記各磁壁は、１８０°磁壁であることを特徴とする請求項１３又は１４に記載の磁区構造解析プログラム。
前記静磁エネルギ算出処理において、有限要素法又は差分法を用いて前記静磁エネルギを算出することを特徴とする請求項１３乃至１５のいずれか１項に記載の磁区構造解析プログラム。
前記静磁エネルギ算出処理において、前記磁区の構造の段階的な変化の度に、前記磁性材料を段階的に複数段階の大きさの領域に、その総数を一定としながら区画し、前記領域毎に静磁エネルギを算出することを特徴とする請求項１３乃至１６のいずれか１項に記載の磁区構造解析プログラム。
前記磁性材料は、電磁鋼板であることを特徴とする請求項１３乃至１７のいずれか１項に記載の磁区構造解析プログラム。
請求項１３乃至１８のいずれか１項に記載の磁区構造解析プログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。