知財求人 - 知財ポータルサイト「IP Force」
特許7260782電磁鋼板の磁気特性予測方法、磁気特性予測装置、およびコンピュータプログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-04-11
(45)【発行日】2023-04-19
(54)【発明の名称】電磁鋼板の磁気特性予測方法、磁気特性予測装置、およびコンピュータプログラム
(51)【国際特許分類】
G01R 33/12 20060101AFI20230412BHJP
G01N 27/72 20060101ALI20230412BHJP
C22C 38/00 20060101ALN20230412BHJP
C22C 38/06 20060101ALN20230412BHJP
H01F 1/147 20060101ALN20230412BHJP
【ＦＩ】
G01R33/12 Z
G01N27/72
C22C38/00 303U
C22C38/06
H01F1/147 175
【請求項の数】
8
(21)【出願番号】P 2019129413
(22)【出願日】2019-07-11
(65)【公開番号】P2021015032
(43)【公開日】2021-02-12
【審査請求日】2022-03-03
(73)【特許権者】
【識別番号】000006655
【氏名又は名称】日本製鉄株式会社
(74)【代理人】
【識別番号】100129838
【弁理士】
【氏名又は名称】山本 典輝
(74)【代理人】
【識別番号】100101203
【弁理士】
【氏名又は名称】山下 昭彦
(74)【代理人】
【識別番号】100104499
【弁理士】
【氏名又は名称】岸本 達人
(72)【発明者】
【氏名】村川 鉄州
(72)【発明者】
【氏名】富田 俊郎
【審査官】續山 浩二
(56)【参考文献】
【文献】特開２００７－３０４６８８（ＪＰ，Ａ）
【文献】特開２０１２－１７３１１５（ＪＰ，Ａ）
【文献】特開２０１２－１７３１１６（ＪＰ，Ａ）
【文献】特開２０１３－１３１０７２（ＪＰ，Ａ）
【文献】特開２０１３－０６４６１７（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｒ ３３／１２
Ｇ０１Ｎ ２７／７２
Ｃ２２Ｃ ３８／００
Ｃ２２Ｃ ３８／０６
Ｈ０１Ｆ １／１４７
(57)【特許請求の範囲】
【請求項1】
外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、前記平均磁化ベクトルＭａｖに基づいて前記電磁鋼板の磁束密度を求める磁束密度計算工程と
を有し、
前記平均磁化計算工程は、
前記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、
前記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、
外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、
平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、
下記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、前記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、
前記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて下記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、前記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、
前記平均磁化ベクトルＭａｖの計算値を前記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、前記平均磁化仮定工程、前記第１計算工程および前記第２計算工程を繰り返し行う反復計算工程と
を有することを特徴とする電磁鋼板の磁気特性予測方法。
【数1】
（上記式（１）において、第一項が静磁エネルギー、第二項が磁気異方性エネルギーを表す。ここで、前記各結晶粒は、前記結晶粒の結晶方位で決まる６種の磁化容易方向のうち、前記外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つ３つの磁区を有すると仮定する。また、Ｈｄは周囲の結晶粒が各結晶粒に及ぼす反磁界ベクトルを表し、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を用いて、下記式（３）により決まるものとする。
【数2】
【数3】
【請求項2】
前記平均磁化計算工程は、前記反磁界係数設定工程前に、前記外部磁場方向の反磁界係数Ｎ／／および前記外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を求める反磁界係数計算工程を有することを特徴とする請求項１に記載の電磁鋼板の磁気特性予測方法。
【請求項3】
前記外部磁場方向の反磁界係数Ｎ／／および前記外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を、Ｎ／／＜Ｎ⊥＜＜１とすることを特徴とする請求項１または請求項２に記載の電磁鋼板の磁気特性予測方法。
【請求項4】
前記平均磁化計算工程前に、前記集合組織を測定し、前記結晶粒情報を得る集合組織測定工程を有し、前記集合組織測定工程では、後方散乱電子回折法、Ｘ線回折法または中性子回折法により、前記集合組織を測定することを特徴とする請求項１から請求項３までのいずれかに記載の電磁鋼板の磁気特性予測方法。
【請求項5】
外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算手段と、前記平均磁化ベクトルＭａｖに基づいて前記電磁鋼板の磁束密度を求める磁束密度計算手段と
を有し、
前記平均磁化計算手段は、
前記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定手段と、
前記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定手段と、
外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定手段と、
平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定手段と、
下記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、前記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算手段と、
前記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて下記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、前記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算手段と、
前記平均磁化ベクトルＭａｖの計算値を前記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、反復計算する反復計算手段と
を有することを特徴とする電磁鋼板の磁気特性予測装置。
【数4】
【数5】
【数6】
【請求項6】
前記集合組織を測定し、前記結晶粒情報を得る集合組織測定手段を有し、前記集合組織測定手段は、後方散乱電子回折装置、Ｘ線回折装置または中性子回折装置であることを特徴とする請求項５に記載の電磁鋼板の磁気特性予測装置。
【請求項7】
電磁鋼板の磁気特性を予測することをコンピューターに実行させるためのコンピュータプログラムであって、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、
前記平均磁化ベクトルＭａｖに基づいて前記電磁鋼板の磁束密度を求める磁束密度計算工程と
をコンピューターに実行させ、
前記平均磁化計算工程は、
前記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、
前記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、
外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、
平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、
下記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、前記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、
前記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて下記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、前記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、
前記平均磁化ベクトルＭａｖの計算値を前記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、前記平均磁化仮定工程、前記第１計算工程および前記第２計算工程を繰り返し行う反復計算工程と
を有することを特徴とするコンピュータプログラム。
【数7】
【数8】
【数9】
【請求項8】
前記Ｋ１、前記Ｎ／／、及び前記Ｎ⊥について、前記電磁鋼板に含まれるＳｉの含有割合（質量％）をＳｉ、主方位の磁化容易方向である<１００>方向と板厚方向とがなす角の最少の角をφＮとしたとき、
Ｋ１＝４１５００－１６６６・Ｓｉ、
Ｎ／／＝０．０００６０・ｃｏｓφＮ―０．０００３５、及び、
Ｎ⊥＝０．００３５５・ｃｏｓφＮ―０．００１４５
で求められる値（ただし、Ｋ１について下二桁は四捨五入した値、Ｎ／／が０以下の時はＮ／／＝０とする）を用いることを特徴とする請求項７に記載のコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、電磁鋼板の磁気特性を集合組織から高精度に予測する磁気特性予測方法、磁気特性予測装置、およびコンピュータプログラムに関するものである。
【背景技術】
【0002】
従来、電動機、発電機、変圧器等の磁心に用いられる電磁鋼板の磁気特性を、材料の集合組織から求める方法が検討されている。しかしながら、その予測精度は十分ではなかった。また、集合組織から磁束密度を求められるにしても、特定の磁化力域における磁束密度が予測できるのみで、その精度も低かった。したがって、広い磁化力の範囲で高い予測精度の磁束密度予測法が求められていた。
【0003】
例えば非特許文献１には、ベクトル法による三次元解析データを用いて、電磁鋼板の集合組織から、その磁気特性を計算する方法が提案されている。非特許文献１では、（１）鉄（単結晶）が磁化されるのは＜１００＞方向のみであること、を前提とし、（２）それぞれの結晶粒は単磁区構造であること（磁化力５０００Ａ／ｍの外部磁界が加わると各結晶粒は単磁区化すること）、（３）それぞれの結晶粒において、外部磁界の向きとのなす角が最も近い磁化容易軸＜１００＞に平行に磁化されていること（各磁区の自発磁化は外部磁界に最も近い「磁化容易方向」を向くこと）、（４）隣接し合う結晶粒の間の相互作用は無視できる程度に十分小さいこと、の仮定の下で、電磁鋼板の磁気特性を計算している。
【0004】
この仮定に基づいて材料の集合組織から求められる磁束密度は、電磁鋼板のＢ５０（磁化力５０００Ａ／ｍにおける磁束密度）に近いことが知られている。しかしながら、この計算には外部磁界の強度は考慮されておらず、計算値が偶然Ｂ５０に近いといってもよい。また、上記以外の磁化力における磁束密度を正確に計算することはできない。
【先行技術文献】
【非特許文献】
【0005】
【文献】立野一郎、「無方向性電磁鋼板の集合組織に基づく磁化の異方性」、鉄と鋼、１９９０年、第７６巻、第１号、ｐ．８１－８８
【発明の概要】
【発明が解決しようとする課題】
【0006】
本発明は、上記問題点に鑑みてなされたものであり、電磁鋼板の集合組織から、広い磁化力の範囲において磁気特性を高精度に予測することができる電磁鋼板の磁気特性予測方法、磁気特性予測装置、およびコンピュータプログラムを提供することを主目的とする。
【課題を解決するための手段】
【0007】
上記目的を達成するために、本発明は、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算工程とを有し、上記平均磁化計算工程は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、下記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて下記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う反復計算工程とを有することを特徴とする電磁鋼板の磁気特性予測方法を提供する。
【0008】
【数1】
【0009】
上記式（１）において、第一項が静磁エネルギー、第二項が磁気異方性エネルギーを表す。ここで、上記各結晶粒は、上記結晶粒の結晶方位で決まる６種の磁化容易方向のうち、上記外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つ３つの磁区を有すると仮定する。また、Ｈｄは周囲の結晶粒が各結晶粒に及ぼす反磁界ベクトルを表し、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を用いて、下記式（３）により決まるものとする。
【0010】
【数2】
【0011】
ここで、上記式（３）において、（）／／と（）⊥は各々、括弧内のベクトルの外部磁場方向成分のベクトルおよび外部磁場方向に垂直な方向成分のベクトルを表す。上記式（１）において、μ０は真空中の透磁率、Ｈｅｘは外部磁場ベクトル、Ｖｉは上記各結晶粒中の磁区の磁区体積、Ｍｉは上記各結晶粒中の磁区の磁化ベクトル、Ｋ１は結晶磁気異方性定数、αｉ，１、αｉ，２、αｉ，３は上記自発磁化ベクトルＭｉの方向余弦を表す。ただし、Ｖ１＋Ｖ２＋Ｖ３＝１である。
【0012】
【数3】
【0013】
上記発明において、上記平均磁化計算工程は、上記反磁界係数設定工程前に、上記外部磁場方向の反磁界係数Ｎ／／および上記外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を求める反磁界係数計算工程を有していてもよい。
【0014】
また本発明においては、上記外部磁場方向の反磁界係数Ｎ／／および上記外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を、Ｎ／／＜Ｎ⊥＜＜１としてもよい。
【0015】
また本発明においては、上記平均磁化計算工程前に、上記集合組織を測定し、上記結晶粒情報を得る集合組織測定工程を有し、上記集合組織測定工程では、後方散乱電子回折法、Ｘ線回折法または中性子回折法により、上記集合組織を測定することが好ましい。
【0016】
また本発明は、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算手段と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算手段とを有し、上記平均磁化計算手段は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定手段と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定手段と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定手段と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定手段と、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算手段と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算手段と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、反復計算する反復計算手段とを有することを特徴とする電磁鋼板の磁気特性予測装置を提供する。
【0017】
また、本発明の電磁鋼板の磁気特性予測装置は、上記集合組織を測定し、上記結晶粒情報を得る集合組織測定手段を有していてもよく、上記集合組織測定手段は、後方散乱電子回折装置、Ｘ線回折装置または中性子回折装置であることが好ましい。
【0018】
また本発明は、電磁鋼板の磁気特性を予測することをコンピューターに実行させるためのコンピュータプログラムであって、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算工程とをコンピューターに実行させ、上記平均磁化計算工程は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う反復計算工程とを有することを特徴とするコンピュータプログラムを提供する。
【0019】
当該コンピュータプログラムにおいて、Ｋ１、Ｎ／／、及びＮ⊥について、電磁鋼板に含まれるＳｉの含有割合（質量％）をＳｉ、主方位の磁化容易方向である<１００>方向と板厚方向とがなす角の最少の角をφＮとしたとき、
Ｋ１＝４１５００－１６６６・Ｓｉ、
Ｎ／／＝０．０００６０・ｃｏｓφＮ―０．０００３５、及び、
Ｎ⊥＝０．００３５５・ｃｏｓφＮ―０．００１４５
で求められる値（ただし、Ｋ１について下二桁は四捨五入した値、Ｎ／／が０以下の時はＮ／／＝０とする）を用いてもよい。
Ｋ１＝４１５００－１６６６・Ｓｉ、
Ｎ／／＝０．０００６０・ｃｏｓφＮ―０．０００３５、及び、
Ｎ⊥＝０．００３５５・ｃｏｓφＮ―０．００１４５
で求められる値（ただし、Ｋ１について下二桁は四捨五入した値、Ｎ／／が０以下の時はＮ／／＝０とする）を用いてもよい。
【発明の効果】
【0020】
本発明は、電磁鋼板の集合組織から、電磁鋼板に要求される、１０００Ａ／ｍ～１００００Ａ／ｍ程度の広い磁化力の範囲において磁束密度を高精度に予測することができるという効果を奏する。したがって、所望の磁気特性を得る電磁鋼板の開発のリードタイムを短縮することができるだけでなく、オンラインまたはオフラインで集合組織を測定し、製品の磁気特性の高精度管理も可能となる。
【図面の簡単な説明】
【0021】
【図1】平均場近似を説明する模式図である。
【図2】結晶座標系の一例を示す模式図である。
【図3】結晶座標系の他の例を示す模式図である。
【図4】本発明の電磁鋼板の磁気特性予測方法の一例を示すフローチャートである。
【図5】本発明における結晶粒情報の一例を示す模式図である。
【図6】本発明の電磁鋼板の磁気特性予測装置の一例を示す模式図である。
【図7】実施例１におけるＥＢＳＤデータである。
【図8】実施例１における、圧延方向から０°～９０°方向の磁束密度を示すグラフである。
【図9】実施例２におけるＥＢＳＤデータである。
【図10】実施例２における、圧延方向から０°～９０°方向の磁束密度を示すグラフである。
【図11】実施例３におけるＥＢＳＤデータである。
【図12】実施例３における、圧延方向から０°～９０°方向の磁束密度を示すグラフである。
【図13】実施例４における無方向性電磁鋼板の断面組織の光学顕微鏡写真である。
【図14】実施例４におけるＯＤＦの断面図である。
【図15】実施例４における、圧延方向から０°～９０°方向の磁束密度を示すグラフである。
【発明を実施するための形態】
【0022】
以下、本発明の電磁鋼板の磁気特性予測方法、磁気特性予測装置、およびコンピュータプログラムについて図面に示す形態例を参照しつつ説明する。
【0023】
Ａ．電磁鋼板の磁気特性予測方法
本発明の電磁鋼板の磁気特性予測方法は、平均場近似による磁気特性予測方法である。
平均場近似とは、個々の局所的状態を決定するときに、周囲からの影響を平均場（個々の局所状態を平均した状態）として捉え、個々の局所状態を計算する方法である。ある平均場を仮定して局所状態が計算できれば、それらを平均することで平均場が求まるので、自己無撞着に計算が行える。したがって、本発明においては、図１に示すように、外部磁場ベクトルＨｅｘの方向に一様に磁化した平均場ベクトルＢａｖの中に、各結晶粒２０が埋め込まれているとして、各結晶粒の磁化容易方向および磁区状態を決定することを考える。つまり、ここでの平均場とは、それによって計算される各結晶粒の磁化容易方向および磁区状態を平均したものである。
本発明の電磁鋼板の磁気特性予測方法は、平均場近似による磁気特性予測方法である。
平均場近似とは、個々の局所的状態を決定するときに、周囲からの影響を平均場（個々の局所状態を平均した状態）として捉え、個々の局所状態を計算する方法である。ある平均場を仮定して局所状態が計算できれば、それらを平均することで平均場が求まるので、自己無撞着に計算が行える。したがって、本発明においては、図１に示すように、外部磁場ベクトルＨｅｘの方向に一様に磁化した平均場ベクトルＢａｖの中に、各結晶粒２０が埋め込まれているとして、各結晶粒の磁化容易方向および磁区状態を決定することを考える。つまり、ここでの平均場とは、それによって計算される各結晶粒の磁化容易方向および磁区状態を平均したものである。
【0024】
本発明においては、以下のような状態を考える。
（１）１０００Ａ／ｍ程度以上の磁化力域では、１８０°磁壁の移動は終了しており、外部磁場の方向に近い自発磁化を持つ３つの磁区が存在する。
（２）各磁区の自発磁化は、結晶磁気異方性および磁界で決められる角度だけ、磁化容易方向から回転している。
（３）周囲の結晶粒との相互作用を、境界に浸み出した磁極による反磁界の影響として捉え計算する。
（１）１０００Ａ／ｍ程度以上の磁化力域では、１８０°磁壁の移動は終了しており、外部磁場の方向に近い自発磁化を持つ３つの磁区が存在する。
（２）各磁区の自発磁化は、結晶磁気異方性および磁界で決められる角度だけ、磁化容易方向から回転している。
（３）周囲の結晶粒との相互作用を、境界に浸み出した磁極による反磁界の影響として捉え計算する。
【0025】
（１）の仮定は、１８０°磁壁が９０°磁壁に比べ動き易いので、１０００Ａ／ｍ程度を超える磁化力域では自然な仮定である。
（２）の仮定は、結晶磁気異方性定数から容易に評価できる。
（３）の仮定は、平均場近似に従えば、図１に示すように、埋め込まれた結晶粒２０が持つ磁化ベクトル（各磁区の持つ磁化ベクトルの和）をＭtotalとすると、Ｍtotalが作る反磁界ベクトルＨｄおよび外部磁界ベクトルＨｅｘの和で作られる磁界中でＭtotalが持つエネルギーが最小になるように各磁区の分量が決まるとすることである。
（２）の仮定は、結晶磁気異方性定数から容易に評価できる。
（３）の仮定は、平均場近似に従えば、図１に示すように、埋め込まれた結晶粒２０が持つ磁化ベクトル（各磁区の持つ磁化ベクトルの和）をＭtotalとすると、Ｍtotalが作る反磁界ベクトルＨｄおよび外部磁界ベクトルＨｅｘの和で作られる磁界中でＭtotalが持つエネルギーが最小になるように各磁区の分量が決まるとすることである。
【0026】
結晶磁気異方性エネルギーＥａは、下記式（４）で表される。
Ｅａ＝Ｋ１・（α１ ２・α２ ２＋α２ ２・α３ ２＋α３ ２・α１ ２） （４）
上記式（４）において、Ｋ１は結晶磁気異方性定数、α１、α２、α３は結晶座標系からみた自発磁化ベクトルの方向余弦である。
図２に結晶座標系を示す。α１、α２、α３はそれぞれ結晶座標系のＸ軸、Ｙ軸、Ｚ軸に対する自発磁化ベクトルＭの方向余弦であり、自発磁化ベクトルＭとＸ軸、Ｙ軸、Ｚ軸とのなす角をそれぞれθ１、θ２、θ３とすると、
α１＝ｃｏｓ（θ１）
α２＝ｃｏｓ（θ２）
α３＝ｃｏｓ（θ３）
で表される。
ここで、Ｋ１は電磁鋼板の化学組成に依存するパラメータである。この値は各種文献などによる値を採用することも可能であるが、本発明では、電磁鋼板に含まれるＳｉの比率をＳｉ（質量％）としたとき、
Ｋ１＝４１５００－１６６６・Ｓｉ
により決定される値を用いることが好ましい。この値を用いる際には下２桁を四捨五入した値とすることができる。
Ｅａ＝Ｋ１・（α１ ２・α２ ２＋α２ ２・α３ ２＋α３ ２・α１ ２） （４）
上記式（４）において、Ｋ１は結晶磁気異方性定数、α１、α２、α３は結晶座標系からみた自発磁化ベクトルの方向余弦である。
図２に結晶座標系を示す。α１、α２、α３はそれぞれ結晶座標系のＸ軸、Ｙ軸、Ｚ軸に対する自発磁化ベクトルＭの方向余弦であり、自発磁化ベクトルＭとＸ軸、Ｙ軸、Ｚ軸とのなす角をそれぞれθ１、θ２、θ３とすると、
α１＝ｃｏｓ（θ１）
α２＝ｃｏｓ（θ２）
α３＝ｃｏｓ（θ３）
で表される。
ここで、Ｋ１は電磁鋼板の化学組成に依存するパラメータである。この値は各種文献などによる値を採用することも可能であるが、本発明では、電磁鋼板に含まれるＳｉの比率をＳｉ（質量％）としたとき、
Ｋ１＝４１５００－１６６６・Ｓｉ
により決定される値を用いることが好ましい。この値を用いる際には下２桁を四捨五入した値とすることができる。
【0027】
自発磁化ベクトルＭが作る反磁界ベクトルＨｄおよび外部磁界ベクトルＨｅｘの和で作られる磁界中で自発磁化ベクトルＭが持つエネルギーを最小にするには、下記式（５）で表される結晶磁気異方性エネルギーと静磁エネルギーとの和の最小化を考えればよい。
【0028】
【数4】
【0029】
上記式（５）において、μ０は真空中の透磁率、Ｈｅｘは外部磁界ベクトル、Ｈｄは自発磁化ベクトルＭが作る反磁界ベクトル、Ｋ１は結晶磁気異方性定数、α１、α２、α３は結晶座標系からみた自発磁化ベクトルＭの方向余弦である。
【0030】
電磁鋼板は体心立方格子（ｂｃｃ）を有するため、磁化容易方向は［１００］［０１０］［００１］［－１００］［０－１０］［００－１］の６種である。このうち、外部磁場の方向に近い３種の磁化容易方向［１００］［０１０］［００１］を考える。自発磁化ベクトルは、結晶磁気異方性および磁界で決められる角度だけ、磁化容易方向から回転していると仮定することから、３種の磁化容易方向［１００］［０１０］［００１］について、磁化容易方向から所定の角度だけ回転している自発磁化ベクトルを持つ３つの磁区が存在すると仮定する。図３に結晶座標系を示す。結晶座標系のＸ軸、Ｙ軸、Ｚ軸はそれぞれ３種の磁化容易方向［１００］［０１０］［００１］のいずれかになる。自発磁化ベクトルＭ１、Ｍ２、Ｍ３はそれぞれ磁化容易方向、つまりＸ軸、Ｙ軸、Ｚ軸から所定の角度だけ回転しており、自発磁化ベクトルＭ１、Ｍ２、Ｍ３を持つ３つの磁区が存在する。すなわち、各結晶粒は、結晶粒の結晶方位で決まる６種の磁化容易方向のうち、外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つ３つの磁区を有すると仮定する。
【0031】
ここで、６種の磁化容易方向のうち、外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つとは、自発磁化ベクトルが磁化容易方向から所定の角度だけ回転していることを前提として、６種の磁化容易方向のうち、外部磁場の方向に近い３種の磁化容易方向から所定の角度だけ回転している自発磁化ベクトルＭｉを持つことを意味する。
【0032】
したがって、結晶粒が持つ磁化ベクトルＭtotalは、３つの磁区が持つ磁化ベクトルの和であり、下記式（６）で表される。
【0033】
【数5】
【0034】
ただし、Ｖ１＋Ｖ２＋Ｖ３＝１である。上記式（６）において、Ｖ１、Ｖ２、Ｖ３は結晶粒中の各磁区の磁区体積、Ｍ１、Ｍ２、Ｍ３は結晶粒中の各磁区の磁化ベクトルである。
【0035】
よって、下記式（１）で表される結晶磁気異方性エネルギーと静磁エネルギーとの和Ｅの最小化を考えればよい。
【0036】
【数6】
【0037】
上記式（１）において、μ０は真空中の透磁率、Ｈｅｘは外部磁界ベクトル、Ｈｄは自発磁化ベクトルＭｉ（ｉ＝１，２，３）が作る反磁界ベクトル、Ｖｉは結晶粒中の各磁区の磁区体積、Ｍｉは結晶粒中の各磁区の磁化ベクトル、Ｋ１は結晶磁気異方性定数、αｉ，１、αｉ，２、αｉ，３は結晶座標系からみた自発磁化ベクトルＭｉの方向余弦である。
【0038】
また、隣接する結晶粒の影響は、平均場近似を考え、下記式（７）に示すように平均場ベクトルＢａｖと結晶粒が持つ磁化ベクトルＭtotalとの差から発生する反磁界ベクトルＨｄで表される。
【0039】
【数7】
【0040】
上記式（７）において、Ｎは反磁界係数、Ｍtotalは結晶粒が持つ磁化ベクトル、Ｂａｖは平均場ベクトル、Ｖ１、Ｖ２、Ｖ３は結晶粒中の各磁区の磁区体積、Ｍ１、Ｍ２、Ｍ３は結晶粒中の各磁区の磁化ベクトルである。
【0041】
上述したように、ここでの平均場ベクトルＢａｖとは、それによって計算される各結晶粒の磁化容易方向および磁区状態を平均したものであり、つまり、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖに真空中の透磁率μ０を掛けたものである。
【0042】
【数8】
【0043】
したがって、周囲の結晶粒が各結晶粒に及ぼす反磁界ベクトルＨｄは、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を用いて、下記式（３）により決まるものとする。
【0044】
【数9】
【0045】
ここで、上記式（３）において、（）／／と（）⊥は各々、括弧内のベクトルの外部磁場方向成分のベクトルおよび外部磁場方向に垂直な方向成分のベクトルを表す。Ｎ／／は外部磁場方向の反磁界係数、Ｎ⊥は外部磁場方向に垂直な方向の反磁界係数、Ｍａｖは外部磁場下における電磁鋼板の材料全体の平均磁化ベクトル、Ｖｉは各結晶粒中の磁区の磁区体積、Ｍｉは各結晶粒中の磁区の磁化ベクトルである。
【0046】
よって、本発明においては、平均場近似に従って、上記式（１）で表される結晶磁気異方性エネルギーと静磁エネルギーとの和が最小になるように、平均場ベクトル、つまり外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める。この際、上記式（１）には求めるべき平均磁化ベクトルＭａｖが含まれているため、まず適当な平均磁化ベクトルＭａｖを仮定し、結晶磁気異方性エネルギーと静磁エネルギーとの和を最小化することで求められる平均磁化ベクトルＭａｖの値が一定値に収斂するまで反復計算する。
【0047】
すなわち、本発明の電磁鋼板の磁気特性予測方法は、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算工程とを有し、上記平均磁化計算工程は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、下記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて下記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う反復計算工程とを有することを特徴とする。
【0048】
【数10】
【0049】
上記式（１）において、第一項が静磁エネルギー、第二項が磁気異方性エネルギーを表す。ここで、上記各結晶粒は、上記結晶粒の結晶方位で決まる６種の磁化容易方向のうち、上記外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つ３つの磁区を有すると仮定する。また、Ｈｄは周囲の結晶粒が各結晶粒に及ぼす反磁界ベクトルを表し、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を用いて、下記式（３）により決まるものとする。
【0050】
【数11】
【0051】
ここで、上記式（３）において、（）／／と（）⊥は各々、括弧内のベクトルの外部磁場方向成分のベクトルおよび外部磁場方向に垂直な方向成分のベクトルを表す。上記式において、μ０は真空中の透磁率、Ｈｅｘは上記外部磁場ベクトル、Ｖｉは上記各結晶粒中の磁区の磁区体積、Ｍｉは上記各結晶粒中の磁区の磁化ベクトル、Ｋ１は結晶磁気異方性定数、αｉ，１、αｉ，２、αｉ，３は上記自発磁化ベクトルＭｉの方向余弦を表す。ただし、Ｖ１＋Ｖ２＋Ｖ３＝１である。
【0052】
【数12】
【0053】
なお、本願明細書においては、単位系としてＳＩ単位系を用いる。
【0054】
図４は本発明の電磁鋼板の磁気特性予測方法の一例を示すフローチャートである。まず、電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与えるステップＳ１（結晶粒情報設定工程）を行う。次に、電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与えるステップＳ２（外部磁場設定工程）を行う。次に、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与えるステップＳ３（反磁界係数設定工程）を行う。次に、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖの仮定値を与えるステップＳ４（平均磁化仮定工程）を行う。
次に、ステップＳ５において、結晶粒番号ｊとして「１」を設定する。次に、結晶粒番号ｊの結晶粒について、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、結晶粒番号ｊの結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるステップＳ６（第１計算工程）を行う。次に、ステップＳ７において、ステップＳ１で与えた全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めたか否かを判定する。全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めていない場合には、ステップＳ８に進む。ステップＳ８に進むと、結晶粒番号ｊとして「ｊ＋１」を設定する。そして、ステップＳ６に戻り、全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるまで、ステップＳ６～Ｓ８を繰り返し行う。
そして、ステップＳ７において、全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めたと判定されると、ステップＳ９（第２計算工程）に進み、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、各結晶粒の磁化ベクトルＭtotalを各結晶粒の体積に比例する重みを掛け平均して平均磁化ベクトルＭａｖの計算値を求める。
次に、ステップＳ１０において、ステップＳ４で与えた平均磁化ベクトルＭａｖの仮定値とステップＳ９で求めた平均磁化ベクトルＭａｖの計算値とを比較し、平均磁化ベクトルＭａｖの値が一定値に収斂したか否かを判定する。平均磁化ベクトルＭａｖの値が収斂していない場合には、ステップＳ１１に進む。ステップＳ１１に進むと、平均磁化ベクトルＭａｖの仮定値として、ステップＳ９で求めた平均磁化ベクトルＭａｖの計算値を与える。そして、ステップＳ５に戻り、平均磁化ベクトルＭａｖの値が収斂するまで、ステップＳ５～Ｓ１１を繰り返し行う。ここまでのステップＳ１～Ｓ１０が平均磁化計算工程である。
そして、ステップＳ１１において、平均磁化ベクトルＭａｖの値が収斂したと判定されると、ステップＳ１２（磁束密度計算工程）に進み、平均磁化ベクトルＭａｖの収斂値に基づいて電磁鋼板の磁束密度を求める。
次に、ステップＳ５において、結晶粒番号ｊとして「１」を設定する。次に、結晶粒番号ｊの結晶粒について、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、結晶粒番号ｊの結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるステップＳ６（第１計算工程）を行う。次に、ステップＳ７において、ステップＳ１で与えた全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めたか否かを判定する。全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めていない場合には、ステップＳ８に進む。ステップＳ８に進むと、結晶粒番号ｊとして「ｊ＋１」を設定する。そして、ステップＳ６に戻り、全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるまで、ステップＳ６～Ｓ８を繰り返し行う。
そして、ステップＳ７において、全ての結晶粒について、結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めたと判定されると、ステップＳ９（第２計算工程）に進み、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、各結晶粒の磁化ベクトルＭtotalを各結晶粒の体積に比例する重みを掛け平均して平均磁化ベクトルＭａｖの計算値を求める。
次に、ステップＳ１０において、ステップＳ４で与えた平均磁化ベクトルＭａｖの仮定値とステップＳ９で求めた平均磁化ベクトルＭａｖの計算値とを比較し、平均磁化ベクトルＭａｖの値が一定値に収斂したか否かを判定する。平均磁化ベクトルＭａｖの値が収斂していない場合には、ステップＳ１１に進む。ステップＳ１１に進むと、平均磁化ベクトルＭａｖの仮定値として、ステップＳ９で求めた平均磁化ベクトルＭａｖの計算値を与える。そして、ステップＳ５に戻り、平均磁化ベクトルＭａｖの値が収斂するまで、ステップＳ５～Ｓ１１を繰り返し行う。ここまでのステップＳ１～Ｓ１０が平均磁化計算工程である。
そして、ステップＳ１１において、平均磁化ベクトルＭａｖの値が収斂したと判定されると、ステップＳ１２（磁束密度計算工程）に進み、平均磁化ベクトルＭａｖの収斂値に基づいて電磁鋼板の磁束密度を求める。
【0055】
本発明においては、上述の（１）～（３）の仮定に基づいて、上記式（１）、（３）を使用して、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求めることにより、外部磁界の強度、磁区構造、および隣接する結晶粒の影響等を考慮して電磁鋼板の磁束密度を求めることができる。したがって、電磁鋼板の材料の集合組織から、電磁鋼板に要求される、１０００Ａ／ｍ～１００００Ａ／ｍ程度の広い磁化力の範囲において磁束密度を精度良く予測することが可能である。したがって、所望の磁気特性を得る電磁鋼板の開発のリードタイムを短縮することができるだけでなく、オンラインまたはオフラインで集合組織を測定し、製品の磁気特性の高精度管理も可能となる。
【0056】
以下、本発明の電磁鋼板の磁気特性予測方法における各工程について説明する。
【0057】
１．平均磁化計算工程
平均磁化計算工程では、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める。外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求めるに際しては、下記の結晶粒情報設定工程、外部磁場設定工程、平均磁化仮定工程、第１計算工程、第２計算工程、および反復計算工程を行う。
平均磁化計算工程では、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める。外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求めるに際しては、下記の結晶粒情報設定工程、外部磁場設定工程、平均磁化仮定工程、第１計算工程、第２計算工程、および反復計算工程を行う。
【0058】
（１）結晶粒情報設定工程
結晶粒情報設定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、また第２計算工程にて、各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求めるに際して、電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える。
結晶粒情報設定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、また第２計算工程にて、各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求めるに際して、電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える。
【0059】
ここで、結晶方位とは、試料座標系と結晶粒座標系との間の相対的な関係を表すものである。試料座標系は、圧延方向（ＲＤ）、板幅方向（ＴＤ）、板厚方向（ＮＤ）を軸とする三次元座標系であり、結晶座標系は、結晶主軸を軸とする三次元座標系である。
【0060】
結晶粒情報は、図５に例示するように、結晶粒番号ｊと、結晶粒番号ｊの結晶粒の結晶方位と、結晶粒番号ｊの結晶粒の体積とが相互に関連付けられたものである。なお、結晶粒番号ｊは整数であり、１から昇順に付けられるものとする。結晶粒番号ｊは、結晶粒を特定する番号である。
【0061】
各結晶粒の結晶方位および体積としては、測定した各結晶粒の結晶方位および体積のデータを用いてもよく、集合組織モデルを適用して、その集合組織モデルの各結晶粒の結晶方位および体積のデータを用いてもよい。
計算に使用する結晶粒の数は、電磁鋼板の種類や集合組織の集積度等により適宜選択される。例えば、方向性電磁鋼板のように集合組織の集積度が高い場合には、計算に使用する結晶粒の数は５０個～５００個程度とすることができる。また、無方向性電磁鋼板のように集合組織の集積度が低い場合には、計算に使用する結晶粒の数は３００個～３０００個程度とすることができる。結晶粒の数が上記範囲内であれば、電磁鋼板の集合組織から磁束密度を精度良く予測することができる。
計算に使用する結晶粒の数は、電磁鋼板の種類や集合組織の集積度等により適宜選択される。例えば、方向性電磁鋼板のように集合組織の集積度が高い場合には、計算に使用する結晶粒の数は５０個～５００個程度とすることができる。また、無方向性電磁鋼板のように集合組織の集積度が低い場合には、計算に使用する結晶粒の数は３００個～３０００個程度とすることができる。結晶粒の数が上記範囲内であれば、電磁鋼板の集合組織から磁束密度を精度良く予測することができる。
【0062】
測定した各結晶粒の結晶方位および体積のデータが用いられる場合、各結晶粒の結晶方位および体積は、後方散乱電子回折（ＥＢＳＤ；Electron Back Scattering Diffraction）法、Ｘ線回折法または中性子回折法により測定することができる。
ＥＢＳＤ法では、各結晶粒の結晶方位および体積を直接測定することができる。
また、Ｘ線回折法、中性子線回折法では、まず極点図を得て、その極点図をもとに結晶方位分布関数を得て、その結晶方位分布関数を用いて各結晶粒の結晶方位および体積を求めることができる。結晶方位分布関数を用いる場合、所望の結晶粒の数の代表方位を選定し、それらに結晶方位分布関数から求められる体積比率を付与して、各結晶粒の結晶方位および体積として使用する。
ＥＢＳＤ法では、各結晶粒の結晶方位および体積を直接測定することができる。
また、Ｘ線回折法、中性子線回折法では、まず極点図を得て、その極点図をもとに結晶方位分布関数を得て、その結晶方位分布関数を用いて各結晶粒の結晶方位および体積を求めることができる。結晶方位分布関数を用いる場合、所望の結晶粒の数の代表方位を選定し、それらに結晶方位分布関数から求められる体積比率を付与して、各結晶粒の結晶方位および体積として使用する。
【0063】
（２）外部磁場設定工程
外部磁場設定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える。
外部磁場の大きさについては、例えば、磁束密度としてＢ５０を求める場合には、５０００Ａ／ｍが設定される。また、外部磁場の方向と電磁鋼板の圧延方向とのなす角度は如何なる角度にも設定することができる。
外部磁場設定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える。
外部磁場の大きさについては、例えば、磁束密度としてＢ５０を求める場合には、５０００Ａ／ｍが設定される。また、外部磁場の方向と電磁鋼板の圧延方向とのなす角度は如何なる角度にも設定することができる。
【0064】
（３）反磁界係数設定工程
反磁界係数設定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える。
上記式（３）において、Ｎ／／およびＮ⊥はそれぞれ外部磁場方向の反磁界係数および外部磁場に垂直な方向の反磁界係数である。この値は各種文献などによる値を採用することも可能であるが、本発明者が検討したところ、Ｎ／／およびＮ⊥は、電磁鋼板の集合組織に比較的強い影響を受け、計算結果に影響を及ぼす。このため、本発明者らは以下のようにＮ／／およびＮ⊥を鋼板の集合組織の主方位から決定できる値として設定することを考え、以下のように決定した。
まず、集合組織が判明しているいくつかの鋼板において、実験値との比較によって、妥当と判断できるＮ／／およびＮ⊥を、次のように計算する。
最初に、磁束密度の実験値から求めた平均磁化ベクトルＭａｖの実験値と、Ｎ／／およびＮ⊥の仮定値とを与え、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求め、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める。平均磁化ベクトルＭａｖの計算値と平均磁化ベクトルＭａｖの実験値とを比較し、平均磁化ベクトルＭａｖの計算値が平均磁化ベクトルＭａｖの実験値に収斂するように、Ｎ／／およびＮ⊥の値を変更して反復計算し、最終的に相当数の事例の集合組織のそれぞれにおいて妥当と判断できるＮ／／およびＮ⊥の値を決定する。
そして、集合組織の主方位の変化とＮ／／およびＮ⊥の値の変化の関連から、Ｎ／／およびＮ⊥を鋼板の集合組織の主方位から決定できるパラメータとして設定する。得られた結果は以下の通りである。
反磁界係数設定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える。
上記式（３）において、Ｎ／／およびＮ⊥はそれぞれ外部磁場方向の反磁界係数および外部磁場に垂直な方向の反磁界係数である。この値は各種文献などによる値を採用することも可能であるが、本発明者が検討したところ、Ｎ／／およびＮ⊥は、電磁鋼板の集合組織に比較的強い影響を受け、計算結果に影響を及ぼす。このため、本発明者らは以下のようにＮ／／およびＮ⊥を鋼板の集合組織の主方位から決定できる値として設定することを考え、以下のように決定した。
まず、集合組織が判明しているいくつかの鋼板において、実験値との比較によって、妥当と判断できるＮ／／およびＮ⊥を、次のように計算する。
最初に、磁束密度の実験値から求めた平均磁化ベクトルＭａｖの実験値と、Ｎ／／およびＮ⊥の仮定値とを与え、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求め、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める。平均磁化ベクトルＭａｖの計算値と平均磁化ベクトルＭａｖの実験値とを比較し、平均磁化ベクトルＭａｖの計算値が平均磁化ベクトルＭａｖの実験値に収斂するように、Ｎ／／およびＮ⊥の値を変更して反復計算し、最終的に相当数の事例の集合組織のそれぞれにおいて妥当と判断できるＮ／／およびＮ⊥の値を決定する。
そして、集合組織の主方位の変化とＮ／／およびＮ⊥の値の変化の関連から、Ｎ／／およびＮ⊥を鋼板の集合組織の主方位から決定できるパラメータとして設定する。得られた結果は以下の通りである。
【0065】
本発明では、Ｎ／／およびＮ⊥は、Ｎ／／＜Ｎ⊥＜＜１であり、主方位の磁化容易方向である<１００>方向と板厚方向とがなす角の最少の角をφＮとする時、
Ｎ／／＝０．０００６０・ｃｏｓφＮ―０．０００３５
Ｎ⊥＝０．００３５５・ｃｏｓφＮ―０．００１４５
により決定できる。ただし、Ｎ／／が０以下の時はＮ／／は０とする。ここで、無方向性電磁鋼板の一般的な集合組織である{１１１}方位粒のｃｏｓφＮは３－０．５であり、理想の集合組織である｛１００｝のｃｏｓφＮは１である。
Ｎ／／＝０．０００６０・ｃｏｓφＮ―０．０００３５
Ｎ⊥＝０．００３５５・ｃｏｓφＮ―０．００１４５
により決定できる。ただし、Ｎ／／が０以下の時はＮ／／は０とする。ここで、無方向性電磁鋼板の一般的な集合組織である{１１１}方位粒のｃｏｓφＮは３－０．５であり、理想の集合組織である｛１００｝のｃｏｓφＮは１である。
【0066】
（４）平均磁化仮定工程
平均磁化仮定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、平均磁化ベクトルＭａｖの仮定値を与える。
平均磁化ベクトルＭａｖの仮定値としては、任意の値でよいが、収束を速めるために、解に近い値を与えることが好ましい。一般に、対象とする電磁鋼板と同様の試料の測定値を用いることで、速い収束が得られる。
平均磁化仮定工程では、第１計算工程にて、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求めるに際して、平均磁化ベクトルＭａｖの仮定値を与える。
平均磁化ベクトルＭａｖの仮定値としては、任意の値でよいが、収束を速めるために、解に近い値を与えることが好ましい。一般に、対象とする電磁鋼板と同様の試料の測定値を用いることで、速い収束が得られる。
【0067】
（５）第１計算工程
第１計算工程では、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める。
第１計算工程では、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める。
【0068】
各結晶粒中の磁区の磁化ベクトルＭｉの長さ（大きさ）は、飽和磁束密度Ｂｓを用いて、下記式（８）により求める。
【0069】
【数13】
【0070】
上記式（８）において、μ０は真空中の透磁率である。Ｂｓは飽和磁束密度であり、電磁鋼板の材料の組成によって決まる値である。
【0071】
各結晶粒は、上記結晶粒情報設定工程にて与えた結晶粒の結晶方位で決まる６種の磁化容易方向のうち、上記外部磁場設定工程にて与えた外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つ３つの磁区を有すると仮定する。例えば、計算に使用する結晶粒の数が５００個である場合、５００個の結晶粒それぞれについて、結晶粒の結晶方位で決まる６種の磁化容易方向のうち、外部磁場の方向に近い３種の磁化容易方向に近い方向に自発磁化ベクトルＭｉ（ｉ＝１，２，３）を持つ３つの磁区を有すると仮定する。
結晶粒は上記の３つの磁区を有すると仮定するため、結晶粒中の磁区体積Ｖｉ（ｉ＝１，２，３）は体積分率であり、Ｖ１＋Ｖ２＋Ｖ３＝１である。
結晶粒は上記の３つの磁区を有すると仮定するため、結晶粒中の磁区体積Ｖｉ（ｉ＝１，２，３）は体積分率であり、Ｖ１＋Ｖ２＋Ｖ３＝１である。
【0072】
上記式（１）において、αｉ，１、αｉ，２、αｉ，３は自発磁化ベクトルＭｉの方向余弦である。図３に示す結晶座標系において、αｉ，１、αｉ，２、αｉ，３はそれぞれ結晶座標系のＸ軸、Ｙ軸、Ｚ軸に対する自発磁化ベクトルＭｉの方向余弦であり、自発磁化ベクトルＭｉとＸ軸、Ｙ軸、Ｚ軸とのなす角をそれぞれθｉ，１、θｉ，２、θｉ，３とすると、
αｉ，１＝ｃｏｓ（θｉ，１）
αｉ，２＝ｃｏｓ（θｉ，２）
αｉ，３＝ｃｏｓ（θｉ，３）
で表される。具体的には、θｉ，１、θｉ，２、θｉ，３は下記式で表される。
αｉ，１＝ｃｏｓ（θｉ，１）
αｉ，２＝ｃｏｓ（θｉ，２）
αｉ，３＝ｃｏｓ（θｉ，３）
で表される。具体的には、θｉ，１、θｉ，２、θｉ，３は下記式で表される。
【0073】
【数14】
【0074】
上記式中、Ｍｉ，Ｘ、Ｍｉ，Ｙ、Ｍｉ，Ｚはそれぞれ自発磁化ベクトルＭｉのＸ軸、Ｙ軸、Ｚ軸方向成分である。
【0075】
自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３は、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和を最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める際に、同時に求めることができる。すなわち、自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３も、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和が最小になるように決められる。
【0076】
また、各結晶粒中の磁区の磁化ベクトルＭｉおよび自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３を求める際には、上記式（１）を簡略化することができ、各磁区について下記式（９）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅｉが最小になればよいことになる。この場合、各結晶粒中の磁区体積Ｖｉとは無関係に、各結晶粒中の磁区の磁化ベクトルＭｉおよび自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３を求めることができる。
【0077】
【数15】
【0078】
上記式（９）において、各記号は上記式（１）と同様である。
よって、まず、上記式（９）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅｉを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３を求め、次に、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和を最小化するように、各結晶粒中の磁区体積Ｖｉを求めることもできる。この場合には、上記式（９）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅｉを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３を求める際に、反磁界ベクトルＨｄ＝０とし、上記式（９）をさらに簡略化してもよい。
よって、まず、上記式（９）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅｉを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３を求め、次に、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和を最小化するように、各結晶粒中の磁区体積Ｖｉを求めることもできる。この場合には、上記式（９）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅｉを最小化するように、各結晶粒中の磁区の磁化ベクトルＭｉおよび自発磁化ベクトルＭｉの方向余弦αｉ，１、αｉ，２、αｉ，３を求める際に、反磁界ベクトルＨｄ＝０とし、上記式（９）をさらに簡略化してもよい。
【0079】
上記式（１）において、Ｋ１は結晶磁気異方性定数であり、電磁鋼板の材料の組成によって決まる値である。例えば、ケイ素鋼板であって、３％Ｓｉ－Ｆｅの組成の場合、結晶磁気異方性定数Ｋ１は３６５００Ｊ／ｍ３である。
【0080】
第１計算工程では、計算に使用する全ての結晶粒について、各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める。
【0081】
（６）第２計算工程
第２計算工程では、上記第１計算工程にて求めた各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める。
第２計算工程では、計算に使用する全ての結晶粒について、各結晶粒の磁化ベクトルＭtotalを求める。
平均磁化ベクトルＭａｖの計算値は、下記式（１０）に示すように、各結晶粒の磁化ベクトルＭtotalについて、各結晶粒の体積Ｗによる加重平均を行って求められる。
第２計算工程では、上記第１計算工程にて求めた各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める。
第２計算工程では、計算に使用する全ての結晶粒について、各結晶粒の磁化ベクトルＭtotalを求める。
平均磁化ベクトルＭａｖの計算値は、下記式（１０）に示すように、各結晶粒の磁化ベクトルＭtotalについて、各結晶粒の体積Ｗによる加重平均を行って求められる。
【0082】
【数16】
【0083】
上記式（１０）において、ｋは計算に使用する結晶粒の数、Ｍj,totalは結晶粒番号ｊの結晶粒の磁化ベクトルＭtotal、Ｗｊは結晶粒番号ｊの結晶粒の体積である。
各結晶粒の体積Ｗｊは、上記結晶粒情報設定工程にて与えた各結晶粒の体積である。
各結晶粒の体積Ｗｊは、上記結晶粒情報設定工程にて与えた各結晶粒の体積である。
【0084】
（７）反復計算工程
反復計算工程では、上記第２計算工程にて求めた平均磁化ベクトルＭａｖの計算値を上記平均磁化仮定工程での平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う。
平均磁化ベクトルＭａｖが収斂するまで反復計算を行う際には、仮定値に対する計算値の変化の割合が１％以下になるまで計算することが好ましく、より好ましくは０．５％以下になるまで計算する。なお、（仮定値に対する計算値の変化の割合）＝｛｜（計算値）－（仮定値）｜｝／（仮定値）×１００である。例えば、平均磁化ベクトルＭａｖの仮定値と計算値とが２桁目まで一致するように反復計算を行えばよい。具体的には、平均磁化ベクトルＭａｖの値は数回の反復計算で収斂する。
反復計算工程では、上記第２計算工程にて求めた平均磁化ベクトルＭａｖの計算値を上記平均磁化仮定工程での平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う。
平均磁化ベクトルＭａｖが収斂するまで反復計算を行う際には、仮定値に対する計算値の変化の割合が１％以下になるまで計算することが好ましく、より好ましくは０．５％以下になるまで計算する。なお、（仮定値に対する計算値の変化の割合）＝｛｜（計算値）－（仮定値）｜｝／（仮定値）×１００である。例えば、平均磁化ベクトルＭａｖの仮定値と計算値とが２桁目まで一致するように反復計算を行えばよい。具体的には、平均磁化ベクトルＭａｖの値は数回の反復計算で収斂する。
【0085】
（８）反磁界係数計算工程
平均磁化計算工程は、上記反磁界係数設定工程前に、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を求める反磁界係数計算工程を有していてもよい。なお、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥の計算については、上記反磁界係数設定工程に記載したので、ここでの説明は省略する。
平均磁化計算工程は、上記反磁界係数設定工程前に、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を求める反磁界係数計算工程を有していてもよい。なお、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥の計算については、上記反磁界係数設定工程に記載したので、ここでの説明は省略する。
【0086】
２．磁束密度計算工程
磁束密度計算工程では、上記反復計算工程にて求めた平均磁化ベクトルＭａｖの収斂値に基づいて電磁鋼板の磁束密度を求める。
磁束密度計算工程では、上記反復計算工程にて求めた平均磁化ベクトルＭａｖの収斂値に基づいて電磁鋼板の磁束密度を求める。
【0087】
３．集合組織測定工程
本発明においては、上記平均磁化計算工程前に、電磁鋼板の集合組織を測定し、各結晶粒の結晶方位および体積を含む結晶粒情報を得る集合組織測定工程を有してもよい。
集合組織の測定方法としては、後方散乱電子回折法、Ｘ線回折法、中性子回折法が挙げられる。なお、これらの方法については、上記結晶粒情報設定工程にて記載したので、ここでの説明は省略する。
本発明においては、上記平均磁化計算工程前に、電磁鋼板の集合組織を測定し、各結晶粒の結晶方位および体積を含む結晶粒情報を得る集合組織測定工程を有してもよい。
集合組織の測定方法としては、後方散乱電子回折法、Ｘ線回折法、中性子回折法が挙げられる。なお、これらの方法については、上記結晶粒情報設定工程にて記載したので、ここでの説明は省略する。
【0088】
４．電磁鋼板
本発明の電磁鋼板の磁気特性予測方法は、無方向性電磁鋼板および方向性電磁鋼板のいずれにも適用可能である。
本発明の電磁鋼板の磁気特性予測方法は、無方向性電磁鋼板および方向性電磁鋼板のいずれにも適用可能である。
【0089】
Ｂ．電磁鋼板の磁気特性予測装置
本発明の電磁鋼板の磁気特性予測装置は、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算手段と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算手段とを有し、上記平均磁化計算手段は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定手段と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定手段と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定手段と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定手段と、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算手段と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算手段と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、反復計算する反復計算手段とを有することを特徴とするものである。
本発明の電磁鋼板の磁気特性予測装置は、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算手段と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算手段とを有し、上記平均磁化計算手段は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定手段と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定手段と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定手段と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定手段と、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算手段と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算手段と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、反復計算する反復計算手段とを有することを特徴とするものである。
【0090】
図６は本発明の電磁鋼板の磁気特性予測装置の一例を示すブロック図である。図６において、電磁鋼板の磁気特性予測装置１は、平均磁化計算手段２と、磁束密度計算手段３と、集合組織測定手段４とを有しており、平均磁化計算手段２は、結晶粒情報設定手段１１と、外部磁界設定手段１２と、反磁界係数設定手段１３と、平均磁化仮定手段１４と、第１計算手段１５と、第２計算手段１６と、反復計算手段１７と有している。電磁鋼板の磁気特性予測装置１にはパーソナルコンピュータ等のコンピュータを用いることができる。
【0091】
電磁鋼板の磁気特性予測装置における各手段は、上記電磁鋼板の磁気特性予測方法における各工程を実行するものであるため、ここでの説明は省略する。
【0092】
Ｃ．コンピュータプログラム
本発明のコンピュータプログラムは、電磁鋼板の磁気特性を予測することをコンピューターに実行させるためのものであって、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算工程とをコンピューターに実行させ、上記平均磁化計算工程は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う反復計算工程とを有することを特徴とするものである。
本発明のコンピュータプログラムは、電磁鋼板の磁気特性を予測することをコンピューターに実行させるためのものであって、外部磁場下における電磁鋼板の材料全体の平均磁化ベクトルＭａｖを求める平均磁化計算工程と、上記平均磁化ベクトルＭａｖに基づいて上記電磁鋼板の磁束密度を求める磁束密度計算工程とをコンピューターに実行させ、上記平均磁化計算工程は、上記電磁鋼板の集合組織を構成する各結晶粒の結晶方位および体積を含む結晶粒情報を与える結晶粒情報設定工程と、上記電磁鋼板に外部から印加される外部磁場の大きさおよび方向を与える外部磁場設定工程と、外部磁場方向の反磁界係数Ｎ／／および外部磁場方向に垂直な方向の反磁界係数Ｎ⊥を与える反磁界係数設定工程と、平均磁化ベクトルＭａｖの仮定値を与える平均磁化仮定工程と、上記式（１）で表される静磁エネルギーおよび磁気異方性エネルギーの和Ｅを最小化するように、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉを求める第１計算工程と、上記各結晶粒中の磁区の磁化ベクトルＭｉおよび磁区体積Ｖｉに基づいて上記式（２）により各結晶粒の磁化ベクトルＭtotalを求め、上記各結晶粒の磁化ベクトルＭtotalを平均して平均磁化ベクトルＭａｖの計算値を求める第２計算工程と、上記平均磁化ベクトルＭａｖの計算値を上記平均磁化ベクトルＭａｖの仮定値として与え、平均磁化ベクトルＭａｖの値が収斂するまで、上記平均磁化仮定工程、上記第１計算工程および上記第２計算工程を繰り返し行う反復計算工程とを有することを特徴とするものである。
【0093】
例えば、ＣＰＵ、ＲＡＭ、ＲＯＭ、ＨＤＤ等を備えるコンピュータにおいて、ＣＰＵが、ＲＯＭやＨＤＤに記憶されているコンピュータプログラムをＲＡＭを使用して実行することによって、上記の各工程が実行される。
【0094】
【0095】
本発明は、上記形態に限定されるものではない。上記形態は例示であり、本発明の特許請求の範囲に記載された技術的思想と実質的に同一な構成を有し、同様な作用効果を奏するものは、いかなるものであっても本発明の技術的範囲に包含される。
【実施例】
【0096】
以下、実施例を挙げて本発明を具体的に説明する。
【0097】
［実施例１］
３％Ｓｉ－０．６％Ｍｎの組成の｛１００｝＜００１＞集合組織を持つ板厚０．３５ｍｍの二方向性電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を単板磁化測定装置を用いて測定した。
また、二方向性電磁鋼板の結晶組織と集合組織をＥＢＳＤによって求めた。組織の例を図７に示す。図７（ａ）は｛１００｝極点図、図７（ｂ）は結晶方位マップである。使用したＥＢＳＤデータは５３個の結晶粒を含み、結晶粒の平均粒径は０．３１ｍｍであった。観察から求めた各結晶粒の体積（面積）と結晶方位を用いて磁気特性の予測計算を行った。計算に用いた結晶粒の個数は約５０個とした。また、この計算では、Ｎ／／およびＮ⊥を各々０．０００２５および０．００５、結晶磁気異方性定数Ｋ１を３％Ｓｉ－Ｆｅの３６５００Ｊ／ｍ３、飽和磁束密度Ｂｓを２．０３Ｔとした。
計算値と実験値の比較を図８に示す。Ｂ８からＢ１００の広い磁化力の範囲で、計算値と実験値がよく一致した。
３％Ｓｉ－０．６％Ｍｎの組成の｛１００｝＜００１＞集合組織を持つ板厚０．３５ｍｍの二方向性電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を単板磁化測定装置を用いて測定した。
また、二方向性電磁鋼板の結晶組織と集合組織をＥＢＳＤによって求めた。組織の例を図７に示す。図７（ａ）は｛１００｝極点図、図７（ｂ）は結晶方位マップである。使用したＥＢＳＤデータは５３個の結晶粒を含み、結晶粒の平均粒径は０．３１ｍｍであった。観察から求めた各結晶粒の体積（面積）と結晶方位を用いて磁気特性の予測計算を行った。計算に用いた結晶粒の個数は約５０個とした。また、この計算では、Ｎ／／およびＮ⊥を各々０．０００２５および０．００５、結晶磁気異方性定数Ｋ１を３％Ｓｉ－Ｆｅの３６５００Ｊ／ｍ３、飽和磁束密度Ｂｓを２．０３Ｔとした。
計算値と実験値の比較を図８に示す。Ｂ８からＢ１００の広い磁化力の範囲で、計算値と実験値がよく一致した。
【0098】
［実施例２］
３％Ｓｉ－０．６％Ｍｎの組成の｛１００｝＜０１２＞集合組織を持つ板厚０．３５ｍｍの電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を単板磁化測定装置を用いて測定した。
また、電磁鋼板の結晶組織と集合組織をＥＢＳＤによって求めた。組織の例を図９に示す。図９（ａ）は｛１００｝極点図、図９（ｂ）は結晶方位マップである。ＥＢＳＤデータは９０個の結晶粒を含み、結晶粒の平均粒径は０．４ｍｍであった。観察から求めた各結晶粒の体積（面積）と結晶方位を用いて磁気特性の予測計算を行った。計算に用いた結晶粒の個数は約９０個とした。また、この計算では、実施例１と同様に、Ｎ／／およびＮ⊥を各々０．０００２５および０．００５、結晶磁気異方性定数Ｋ１を３％Ｓｉ－Ｆｅの３６５００Ｊ／ｍ３、飽和磁束密度Ｂｓを２．０３Ｔとした。
計算値と実験値の比較を図１０に示す。Ｂ１０からＢ１００の広い磁化力の範囲で、計算値と実験値がよく一致した。
３％Ｓｉ－０．６％Ｍｎの組成の｛１００｝＜０１２＞集合組織を持つ板厚０．３５ｍｍの電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を単板磁化測定装置を用いて測定した。
また、電磁鋼板の結晶組織と集合組織をＥＢＳＤによって求めた。組織の例を図９に示す。図９（ａ）は｛１００｝極点図、図９（ｂ）は結晶方位マップである。ＥＢＳＤデータは９０個の結晶粒を含み、結晶粒の平均粒径は０．４ｍｍであった。観察から求めた各結晶粒の体積（面積）と結晶方位を用いて磁気特性の予測計算を行った。計算に用いた結晶粒の個数は約９０個とした。また、この計算では、実施例１と同様に、Ｎ／／およびＮ⊥を各々０．０００２５および０．００５、結晶磁気異方性定数Ｋ１を３％Ｓｉ－Ｆｅの３６５００Ｊ／ｍ３、飽和磁束密度Ｂｓを２．０３Ｔとした。
計算値と実験値の比較を図１０に示す。Ｂ１０からＢ１００の広い磁化力の範囲で、計算値と実験値がよく一致した。
【0099】
［実施例３］
３％Ｓｉ－０．９％Ｍｎの組成の弱い｛１００｝集合組織を持つ板厚０．３５ｍｍの電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を単板磁化測定装置を用いて測定した。
また、電磁鋼板の結晶組織と集合組織をＥＢＳＤによって求めた。組織の例を図１１に示す。図１１は方位分布関数（Orientation Distribution Function；ＯＤＦ）のφ２＝４５°断面図である。使用したＥＢＳＤデータは３１０個の結晶粒を含んでいた。観察から求めた各結晶粒の体積（面積）と結晶方位を用いて磁気特性の予測計算を行った。計算に用いた結晶粒の個数は約３００個とした。また、この計算では、実施例１、２と同様に、Ｎ／／およびＮ⊥を各々０．０００２５および０．００５、結晶磁気異方性定数Ｋ１を３％Ｓｉ－Ｆｅの３６５００Ｊ／ｍ３とし、飽和磁束密度Ｂｓは２．００Ｔとした。
計算値と実験値の比較を図１２に示す。計算値と実験値がよく一致した。
また、立野モデルによる計算値も図１２に示すが、立野モデルの結果は、Ｂ５０の実験値に比較して０．０５～０．１Ｔ程度大きな値となった。
３％Ｓｉ－０．９％Ｍｎの組成の弱い｛１００｝集合組織を持つ板厚０．３５ｍｍの電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を単板磁化測定装置を用いて測定した。
また、電磁鋼板の結晶組織と集合組織をＥＢＳＤによって求めた。組織の例を図１１に示す。図１１は方位分布関数（Orientation Distribution Function；ＯＤＦ）のφ２＝４５°断面図である。使用したＥＢＳＤデータは３１０個の結晶粒を含んでいた。観察から求めた各結晶粒の体積（面積）と結晶方位を用いて磁気特性の予測計算を行った。計算に用いた結晶粒の個数は約３００個とした。また、この計算では、実施例１、２と同様に、Ｎ／／およびＮ⊥を各々０．０００２５および０．００５、結晶磁気異方性定数Ｋ１を３％Ｓｉ－Ｆｅの３６５００Ｊ／ｍ３とし、飽和磁束密度Ｂｓは２．００Ｔとした。
計算値と実験値の比較を図１２に示す。計算値と実験値がよく一致した。
また、立野モデルによる計算値も図１２に示すが、立野モデルの結果は、Ｂ５０の実験値に比較して０．０５～０．１Ｔ程度大きな値となった。
【0100】
［実施例４］
１．８％Ｓｉ－０．２％Ａｌを含有する、板厚０．５ｍｍの市販の無方向性電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を、単板磁化測定装置を用いて測定した。
また、電磁鋼板の結晶組織と集合組織を各々光学顕微鏡およびＸ線回折によって求めた。組織の例を図１３に示す。図１３は無方向性電磁鋼板の断面組織の光学顕微鏡写真である。上記の実施例１～３に比較して、結晶粒径は小さく、平均粒径は５０μｍ程度であった。
図１４は（ａ）板の表面、（ｂ）板の表面から板厚の１／４深さ（ｔ／４）、（ｃ）板厚の中心で測定した集合組織をＯＤＦの断面図として表したものである。左上から右下に向かって、φ２＝０°～９０°の範囲で、φ２を５°刻みの断面で表している。ＯＤＦは、Ｘ線回折によって求めた｛２００｝、｛１１０｝、｛２１１｝極点図から球調和関数展開法を用いて得た。このＯＤＦから１２００個の結晶方位をランダムに抽出し、各結晶方位のＯＤＦの値から、各結晶方位の結晶粒の体積を求め、それら結晶方位と結晶粒の体積を計算に用いた。ここでは、結晶粒の体積はＯＤＦの値にｓｉｎ（Φ）を掛けたものに比例するとし、１２００個の結晶粒の体積の和を１とした。ここで、φ１、Φ、φ２は結晶方位を表すオイラー角である。この際、板の表面から板厚の１／４深さ（ｔ／４）の影響は表面と板厚中心に比べ２倍あると仮定し、この深さの結晶粒の体積を２倍にして計算した。
また、Ｎ／／およびＮ⊥を各々０．０および０．００３５、結晶磁気異方性定数Ｋ１を３８５００Ｊ／ｍ３とし、飽和磁束密度Ｂｓは２．０９Ｔとした。
計算値と実験値の比較を図１５に示す。Ｂ２５の計算値が実験値をわずかに下回るものの、計算値と実験値はよく一致した。
１．８％Ｓｉ－０．２％Ａｌを含有する、板厚０．５ｍｍの市販の無方向性電磁鋼板を準備し、短冊状の試験片を圧延方向から種々の角度に切り出し、圧延方向に対して０°～９０°方向の磁気特性を、単板磁化測定装置を用いて測定した。
また、電磁鋼板の結晶組織と集合組織を各々光学顕微鏡およびＸ線回折によって求めた。組織の例を図１３に示す。図１３は無方向性電磁鋼板の断面組織の光学顕微鏡写真である。上記の実施例１～３に比較して、結晶粒径は小さく、平均粒径は５０μｍ程度であった。
図１４は（ａ）板の表面、（ｂ）板の表面から板厚の１／４深さ（ｔ／４）、（ｃ）板厚の中心で測定した集合組織をＯＤＦの断面図として表したものである。左上から右下に向かって、φ２＝０°～９０°の範囲で、φ２を５°刻みの断面で表している。ＯＤＦは、Ｘ線回折によって求めた｛２００｝、｛１１０｝、｛２１１｝極点図から球調和関数展開法を用いて得た。このＯＤＦから１２００個の結晶方位をランダムに抽出し、各結晶方位のＯＤＦの値から、各結晶方位の結晶粒の体積を求め、それら結晶方位と結晶粒の体積を計算に用いた。ここでは、結晶粒の体積はＯＤＦの値にｓｉｎ（Φ）を掛けたものに比例するとし、１２００個の結晶粒の体積の和を１とした。ここで、φ１、Φ、φ２は結晶方位を表すオイラー角である。この際、板の表面から板厚の１／４深さ（ｔ／４）の影響は表面と板厚中心に比べ２倍あると仮定し、この深さの結晶粒の体積を２倍にして計算した。
また、Ｎ／／およびＮ⊥を各々０．０および０．００３５、結晶磁気異方性定数Ｋ１を３８５００Ｊ／ｍ３とし、飽和磁束密度Ｂｓは２．０９Ｔとした。
計算値と実験値の比較を図１５に示す。Ｂ２５の計算値が実験値をわずかに下回るものの、計算値と実験値はよく一致した。
【0101】
以上のように、本発明によって、材料の集合組織から電磁鋼板の磁束密度を精度良く予測できることが確認された。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
