電磁場解析システム
【課題】軟磁性材料の透磁率や磁気抵抗率に対する、異方性と弾性曲げの影響を考慮に入れてマックスウェル方程式を解く電磁場解析システムにより、鉄損を評価する。
【解決手段】微小領域における軟磁性材料の予め定められた方向と磁束密度の方向との間の成す角度θおよび弾性曲げ曲率Kを異方性のパラメータとして、H−B曲線を格納するデータベースとW−B曲線を格納するデータベースと、微小領域においてマックスウェル方程式に基づき、前記角度θ、および、磁束密度の大きさBを決定する磁束密度ベクトル決定手段と、微小領域の鉄損を計算する鉄損計算手段と、前記微小領域の鉄損の総和を求める鉄損総和手段とを有することを特徴とする電磁場解析システム。
【解決手段】微小領域における軟磁性材料の予め定められた方向と磁束密度の方向との間の成す角度θおよび弾性曲げ曲率Kを異方性のパラメータとして、H−B曲線を格納するデータベースとW−B曲線を格納するデータベースと、微小領域においてマックスウェル方程式に基づき、前記角度θ、および、磁束密度の大きさBを決定する磁束密度ベクトル決定手段と、微小領域の鉄損を計算する鉄損計算手段と、前記微小領域の鉄損の総和を求める鉄損総和手段とを有することを特徴とする電磁場解析システム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、鋼板をはじめとする軟磁性材料の異方性および弾性曲げ変形による磁気特性劣化を考慮した変形・電磁場連成解析システムに関し、特に軟磁性材料の弾性曲げ変形に関する磁気特性の劣化を考慮した電気機器の鉄損を予測する方法に関する。
【背景技術】
【0002】
環境問題に対する意識の高まりにより、受配電に使用される変圧器ならびに空調機やハイブリット自動車に使用されるモータに対し、鉄損や銅損に代表される損失の低減の要求が高まっている。それに対し、軟磁性材料として実用上多く使用されている電磁鋼板は、鋼板圧延方向からの角度、励磁方向への応力および、弾性変形による曲げにより、励磁特性や鉄損特性が劣化する。このため、変圧器やモータ等の電気機器の性能を高精度に予測するためには、このような使用している材料の異方性や励磁方向への応力、そして、弾性変形による磁気特性の劣化を考慮することが欠かせない。
【0003】
従来は、励磁方向への応力や弾性変形による曲げを考慮した磁気特性を評価する場合、実機試作材料の実測を行いエプスタイン試験法による鉄損特性と単板磁気測定法(SST法)による鉄損特性との差を励磁方向への応力や弾性変形による応力に起因する磁気特性の劣化分として算出していた。この方法では、励磁方向への応力の影響による磁気特性の劣化分と弾性変形による応力に起因する磁気特性の劣化分を分離することができなかった。
【0004】
さらに、従来の弾性変形による曲げの影響評価方法は実機試作を前提としていたため、時間的、金銭的なコストがかかることに加え、形状を連続的に変更し試作することが困難であるため、電気機器設計時に損失最小となる最適形状を求めることが非常に困難であった。
【0005】
一方、従来よりマックスウエル方程式による電磁場解析法は、鋼板等の軟磁性材料の鉄損評価に用いられてきた。以下ここでは、軟磁性材料として実用上多く使用されている鋼板でもって代表的に説明することにする。非特許文献2に示されるように、電磁場解析には、コンピュータが使用され、鋼板の形状、計算のために分割された微小領域の大きさ、磁界Hに対する磁束密度B、周波数等は、コンピュータによる計算の解を求めるための物理量パラメータ条件として採用されている。即ち、これらの条件を考慮に入れて、マックスウエル方程式の計算機解が得られる。
【0006】
マックスウエル方程式の計算機解に、磁束密度Bに対応して測定された鉄損Wのデータ(W−B曲線)を考慮に加えて鉄損が求められている。例えば図1は、従来技術に基づいて実行される鉄損の数値計算ルーチンの流れを示すブロック図である。まず、鋼板の形状、微小領域分割、周波数等をパラメータとしてH−B曲線など、解を求めるための条件を与え、計算によってマックスウエル方程式の数値解を求める。磁束密度Bに対応して測定された、微小領域における鋼材の鉄損Wのデータ(W−B曲線)を上の数値解に与え、鋼板全体の損失Wを算出する。ここで、鋼板に作用する弾性曲げ変形の影響は無視され、鋼板に曲げ変形が加わっても磁気特性は変わらないと仮定している。
【0007】
特許文献1には、データ入力部と連成解析部と結果出力部より構成される磁気デバイスの磁気特性解析方法について記載されている。ここで、データ入力部より入力されるデータは、磁気デバイスを構成される材料の特性に関わるデータ、磁気デバイスの構造に関わるデータ、構造解析の境界条件に関わるデータ、磁場解析の境界条件に関わるデータである。連成解析部において、データ入力部から入力された前記境界条件に関わるデータに基づき区分された複数部分ごとの応力分布を求め、前記磁気デバイスの電磁場解析の境界条件に関わるデータ及び前記応力分布に基づき前記複数部分ごとの磁気特性を求め、前記複数部分ごとの磁気特性に基づき磁気デバイス全体の磁気特性を求めるものである。ここでは、磁気デバイスで用いられる磁性材料や非磁性材料の荷重や熱応力による歪の影響は考慮されているものの、磁性材料に作用する弾性曲げの影響は無視され、鋼板に曲げ変形が加わっても磁気特性は変わらないと仮定している。
【0008】
さらに、特許文献2には、磁性体の予め与えられた方向と磁束密度の方向との間の角度θおよび応力σを異方性のパラメータとして磁束密度と磁界とを関係づける解析式およびデータに基づくH−B曲線、および、磁束密度と鉄損とを関係づける解析式およびデータに基づくW−B曲線を用いて磁束密度ベクトル決定手段と鉄損総和手段により、磁気特性を求める方法について記載されている。ここでも同様に、磁性体である電磁鋼板の応力の影響は考慮されているものの、磁性材料に作用する弾性曲げの影響は無視され、鋼板に曲げ変形が加わっても磁気特性が変わらないと仮定している。
【0009】
ここで、マックスウエル方程式の解法について簡単に説明する。公知の多くの文献から明らかなように、マックスウエル方程式は次式で与えられる。
【0010】
【数1】
【0011】
ここで、B、H、D、E、Jはそれぞれ磁束密度、磁界、電束密度、電界、電流密度である。また、ρは電荷密度である。B、H、D、E、Jの間には、次の関係がある。
【0012】
【数2】
【0013】
ここで、μ、ε、σはそれぞれ、透磁率、誘電率、導電率である。
【0014】
一方、非特許文献2によれば、電磁界に関する解析が詳細に記載されている。同文献よれば、∂D/∂tは無視されている。磁束の発散は常に零であるので、連続であり、磁気ベクトルポテンシャルAが次式によって与えられている。
【0015】
【数3】
【0016】
これらの式から
【0017】
【数4】
【0018】
が得られている。従って、
【0019】
【数5】
【0020】
が得られる。ここで、−gradφ=E、J0は外部からの強制電流密度、Jeはうず電流密度、テンソル量で与えられる磁気抵抗率[ν]は、[ν]=1/[μ]である。(5)式は、ガラーキン法(Galerkin Method)により2次元的、及び3次元的に解かれる。実際には、透磁率は、一般的に弾性曲げ変形の影響を受ける。しかし、従来の鉄損評価法では、透磁率が弾性曲げ変形の影響を受けないと仮定して、数値解析が行われ、解が求められていた。
【0021】
また、H−B曲線は曲率に対して影響がなく、一定であると仮定している。このようにして、与えられた磁界強度Hに対する微小領域ごとに求められた、H−B曲線から、各微小領域における磁束密度の分布が求められる。一方、磁束密度Bと鉄損Wの関係は曲率Kに対しては一定としており、W−B曲線の形でデータベースに格納されている。従って、上記数値計算によって求められた磁束密度分布をデータベースに格納されたW−B曲線に適用すれば、与えられた磁界強度Hに対する鋼板全体の鉄損が数値計算によって求められる。上記従来技術では、鉄損を求めるための電磁場解析において、透磁率は等方性とし、また、弾性曲げ変形の影響を無視した透磁率でもって計算を行っているため、実際の鉄鋼材料において存在する透磁率の異方性や曲率の影響が無視され、そのために実測値の計算値からの乖離が無視できず、鉄損を十分に評価できないという欠点があった。
【特許文献1】特開平8−249621号公報
【特許文献2】特許3676761号公報
【非特許文献1】Nippon Steel Corporation:「Technical Data on Domain Refined ORIENTCORE・HI−B ZDKH Second Edition」
【非特許文献2】中田、高橋:「電気工学の有限要素法」(森北出版、1982)
【発明の開示】
【発明が解決しようとする課題】
【0022】
そこで、本発明は、このような従来技術の問題点を解決し、従来の解決方法を踏襲しながら、軟磁性材料の物性値パラメータに対する、異方性と弾性曲げの影響を考慮に入れた電磁場解析システムを提供して、正確な鉄損評価を実現することを課題とする。
【課題を解決するための手段】
【0023】
上記目的を達成するため、本発明による電磁場解析システムは、軟磁性材料の電磁場解析対象領域を複数の微小領域に分割する領域分割手段と、該微小領域における磁性体の予め定められた方向と磁束密度Bの方向との間の成す角度θおよび軟磁性材料の弾性曲げ曲率Kを異方性のパラメータとして、磁束密度Bと磁界Hとを関係付けるH−B曲線および磁束密度Bと鉄損Wとを関係付けるW−B曲線を格納するデータベースと、該データベースに格納されているH−B曲線を基にして、前記微小領域におけるマックスウエル方程式および弾性曲げ曲率Kに基づいて、前記角度θおよび前記磁束密度Bの大きさを決定する磁束密度ベクトル決定手段と、前記データベースに格納されているW−B曲線を基にして、前記微小領域の鉄損Wを計算する鉄損計算手段と、前記微小領域の鉄損Wの総和を求める鉄損総和手段とを有することを特徴とする。
【0024】
また、前記磁束密度Bおよび磁界Hは、該磁束密度Bと磁界Hとの位相差θBHを考慮して関係付けられることを特徴とする。
【0025】
また、前記磁束密度Bは、時間高調波成分を含む磁束密度であることを特徴とする。
【0026】
さらに、前記磁束密度ベクトル決定手段は、前記弾性曲げ曲率Kを算出する弾性曲げ曲率分布算出手段を備えることを特徴とする。
【0027】
さらに、前記軟磁性材料は、電気機器を構成する軟磁性材料であり、前記弾性曲げ曲率Kは、該電気機器の製作時あるいは自重により発生する弾性曲げ曲率であることを特徴とする。
【0028】
さらに、前記磁束密度ベクトル決定手段は、前記微小領域における磁束密度ベクトルを、その大きさBmaxと、前記の磁性体の予め定められた方向とのなす角度θとに分解し、前記H−B曲線から、前記角度θにおけるH−B曲線を導出し、前記微小領域における磁束密度およびその前記角度θ、磁界とが、その曲線上に存在するように、収束計算を行うことを特徴とする。
【0029】
さらに、前記H−B曲線およびW−B曲線は、前記弾性曲げ曲率Kをパラメータとして、磁束密度と磁界の関係より磁界増分係数α(=HK/HK=0)、磁束密度と鉄損値の関係より鉄損増分係数β(WK/WK=0)を求め、前記弾性曲げ曲率K=0のときのH−B曲線およびW−B曲線に前記磁界増分係数αと鉄損増分係数βとを乗じて得られるH−B曲線およびW−B曲線であることを特徴とする。
【発明の効果】
【0030】
本発明は、圧延による軟磁性材料の異方性と弾性曲げ曲率のH−B曲線およびW−B曲線への影響を考慮した鉄損の評価手段を見出し、電磁場解析システムとして具現化したものである。すなわち、本発明は鉄鋼材料の全領域を微小領域に分割し、分割された各微小領域に対して予め与えられた方向と磁束密度との成す角度と、例えば加工や自重の影響で生じた弾性変形による曲率とをパラメータとして考慮に入れた鉄損データをデータベースから求めて計算し、更に各微小領域で求められた鉄損データの総和をとることにより、異方性が弾性曲げの影響の大きい材料であっても、磁束密度の異方性・弾性曲げ影響や磁界に対する磁束密度の非直線性を意識しないで鉄損を評価計算することができる。また、磁束密度Bと磁界Hとの位相差θBH、時間高調波を考慮した電磁場解析を行うことが出来るなど産業上有用な、著しい効果を奏する。
【発明を実施するための最良の形態】
【0031】
以下、図面を参照して本発明を詳細に説明する。
【0032】
鉄鋼材料には鋼材の予め定められた方向、例えば圧延方向とそれに直角の方向とでは透磁率に異方性がある。このような鉄鋼材料の透磁率の異方性を考慮に入れれば、鉄損を精度よく評価することが可能となる。まず、電磁場解析の対象とする全領域を複数の微小領域に分割する。各微小領域の内部で透磁率が等方性を満足する場合は、各微小領域間で透磁率μ及び磁気抵抗率νや、予め定められた方向、例えば圧延方向と磁束密度Bの方向との間の成す角度θに差があっても容易に計算をすることができる。透磁率に異方性がある場合は、各微小領域にそれぞれ角度θとそれに対する透磁率を与えれば、コンピュータを利用してマックスウェル方程式の解を求めることができる。
【0033】
また、鉄鋼材料には、モータ製造工程の一つである焼嵌めや自重により、鋼板が弾性変形している場合が多く、弾性変形により、透磁率をはじめとした磁気特性が大きく変化する。このような鉄鋼材料の透磁率に対する弾性変形の影響を考慮に入れれば、鉄損を精度よく評価することが可能となる。この場合、電磁場解析の対象とする全領域を複数の微小領域に分割し、各微小領域の内部で弾性曲げ曲率が一定であるとすれば、その微小区間内では同一の磁気特性、透磁率を用いても差し支えない。このように、各微小領域にそれぞれと弾性曲げ曲率とそれに対する透磁率の関係を与えれば、コンピュータを利用してマックスウェル方程式の解を求めることができる。
【0034】
図2は、鉄鋼材料の対象とする全領域を格子状の複数の微小領域に分割した模様を示す説明図である。図2において、鉄鋼材料は、1,2,3……,i,i+1,i+2,……,j,j+1,j+2,……,nの微小領域に分割してある。例えば、i番目の微小領域の磁束密度がΔBiであるとし、鉄損がΔWiであるとする。
【0035】
図2において分割された各微小領域の内部においては、透磁率μ及び磁気抵抗率νは一定の値であるとする。このようにすれば、有限要素法において、各微小領域間の境界では不連続であっても、領域内では一様なパラメータをもっていると考えることができる。従って、異方性・弾性曲げの影響や磁界に対する磁束密度の非直線性を有する鉄鋼材料の鉄損の計算においても、予め計算された各微小領域の鉄損の総和を求めることによって全体鉄損を容易に計算することができる。図2では、全体の鉄損W(watts)は
【0036】
【数6】
【0037】
で与えられる。ここで、Δwiは微小領域内の鉄損である。また、ΔBiの方向θと大きさΔBiは、各微小領域によってそれぞれ異なった値をとる。鉄損Wは磁束密度Bが大きい程、大きな値をとるが、必ずしも直線関係になるわけではない。上記領域分割は有限要素法が適用されることを前提にして実施したものであるが、差分法或いはその他、類似の計算方式に適用可能であることは云うまでもない。
【0038】
図3は、磁束密度の異方性を示す説明図である。以下弾性曲げ曲率だけでなく磁束密度の異方性をも考慮した場合についても考えてみる。RDは鋼板の圧延方向を(rolling direction)、TDはそれに直角な方向(transversal direction)である。磁束密度Bの方向と鋼板の圧延方向RDとの間の角度をθとすれば、H−B曲線及びW−B曲線のθ依存性は、例えば図4(a)、(b)に示すような形状で与えられる。ここで、鋼板の圧延方向の代わりに、任意の予め与えられた方向とすることも可能である。図4(a)、(b)はいずれも角度θの依存性を表したものであるが、複数葉のデータは更に曲率Kに対する依存性も示している。
【0039】
図4(a)、(b)では、θ=0°,15°,30°……のデータがH−B関数及びW−B関数として与えられているが、その他の角度θではデータが関数として与えられていない。この場合、θ=20°において求める必要のある磁束密度と磁界との関係は、θ=15°及びθ=30°におけるデータから求められた近似関数について、以下で詳細に述べる補間内挿して求めたり、ニュートンラプソン法によるテイラー展開からマックスウェル方程式との連成により繰り返し計算して求められる。一方、磁束密度が与えられたときの鉄損は、θ=15°及びθ=30°におけるデータから補間内挿して求められる。ここで、曲率Kの逆数の変化が十分に無視できるほど小さい場合は、曲率は一定であると仮定してかまわない。もし、曲率Kの逆数が十分に大きく変化する場合には、それぞれの曲率に対応したH−B曲線及びW−B曲線を選択して繰り返し計算を実行する。ここで、H−B曲線及びW−B曲線は、データに基づいて求める方法について記したが、理論的に得られた解析式およびデータから得られた回帰式といった解析式を用いても構わない。
【0040】
次に図4(a)、(b)に含まれている曲率Kに対する依存性に関して少し詳細に説明する。例えば、図4(a)、(b)では、0[1/m]<K<20[1/m]の範囲のデータを取り扱い、サンプル化された曲率Kの値に対するデータとして示している。従って、データの表に直接与えられていない曲率に対する磁束密度や鉄損は、角度θの場合と同様にニュートンラプソン法によるテイラー展開から補間内挿して求められる。すなわち、(Bi,θi,ρ)によって微小領域iにおける鉄損wiが求められる。
【0041】
例えば、図4(b)に示すK=20[1/m]のデータにおいて、θ=15°及び30°の曲線、及びB=1.0T(Tesla)及び1.2Tによって切断され、点ABCDによって囲まれた領域を仮定し、点P(θ=20°、B=1.1T)における鉄損を求めるものとする。まず、θ=15°及び30°の曲線から内挿近似によって級数展開を行い、θ=20°のW−B曲線を求める。次にB=1.0T及び1.2Tの直線(曲線でもよい。)によって上記θ=20°のW−B曲線を切断、B軸上の内挿によってB=1.1Tにおける点Pを確定する。確定された点Pの縦軸から鉄損Wを読取ることができる。実際の処理では、図4(b)のグラフはデータベース内に格納されているので、コンピュータのソフトウエア処理によって内挿計算を行い、鉄損を求めることができる。
【0042】
曲率Kの影響によってH−B曲線が変化することは、図4によって示したとおりである。例えば、コンプレッサ等で使用される凹型鉄芯は鋼板をせん断して成形され、70℃から90℃くらいに加熱された軟鋼製の円筒リングに挿入されることで焼嵌めにより固定され、電気機器である回転電機に使用される。円筒リングが常温まで冷却収縮することにより、凹型鉄芯は、内部に圧縮応力を受けるとともに、弾性変形する。この弾性変形量を曲率Kにて表現する。図21は、凹型鉄芯における弾性変形による曲率の磁束密度への影響例を示した説明図である。図21において、xは鋼板のある方向としての圧延方向(RD)、Bは磁束の接線方向であり磁束密度を表す、θはx方向と磁束密度Bとのなす角度を示す。空間上の指定された位置を決定すると、磁束密度ベクトルBが与えられ、曲率Kの値が決定される。すなわち、指定された空間位置において、H−B曲線より、(θ、K)に対応する磁束密度Bが求められ、これによって、W−B曲線に(B、θ、K)を適用することにより鉄損Wが算出される。空間位置に対応する曲率は、構造力学の理論に従って有限要素法により算出することが可能である。
【0043】
さらに、電気機器の一つである変圧器においては、方向性電磁鋼板をせん断加工後、図22のように組み上げてコアを形成する。その際、コアの自重により鋼板が弾性変形する。弾性変形した際の取り扱いの方法は、前記の焼嵌め固定されたコアと同様である。
【0044】
図5は、磁束密度と磁界のベクトルの位相差を示す図である。図5において、磁界Hと磁束密度Bはベクトルで、その間に位相差θBHがある。このとき、磁気抵抗率νは簡略的に次式のようにモデル化できる。
【0045】
【数7】
【0046】
このとき、磁束密度ベクトルおよび磁界ベクトルは次式で表される。
【0047】
【数8】
【0048】
従って、磁気抵抗率νは、Bおよびθの関数として
【0049】
【数9】
【0050】
であらわされ、図6のようなH−B曲線から求めることができる。また、磁束密度Bと磁界Hとの位相差θBHも、Bおよびθの関数として
【0051】
【数10】
【0052】
で表されるので、図7のようなB−θBH曲線をデータベースとして用いることによりθBHを考慮した磁束密度ベクトルを決定することができる。
【0053】
図8は、磁束密度Bの時間的変化を示す図である。図8において、横軸は時間、縦軸は磁束密度を示している。例えば、磁界Hが正弦波であっても、磁束密度Bは正弦波とはならず時間高調波成分を含むため図8のような歪が生じているので、正確な電磁場解析を行うためにはこの時間高調波を考慮する必要がある。そこで、時間高調波成分を含んだ磁束密度Bをフーリエ級数展開し、その周波数ごとに図9のようなH−B曲線と図10のようなW−B曲線をデータベースとして用いることによって時間高調波成分を考慮した電磁場解析を行うことができる。このH−B曲線およびW−B曲線から磁気抵抗率νおよび鉄損Wは、B,θおよび周波数fの関数として、次式で与えられる。
【0054】
【数11】
【0055】
また、時間高調波成分を含んだ磁束密度Bをフーリエ級数展開すると、Bのx成分、y成分はそれぞれ次式で表される。
【0056】
【数12】
【0057】
これにより、高調波における第i番目の磁束密度ベクトルおよび、その位相差は次式によって定義できる。
【0058】
【数13】
【0059】
【数14】
【0060】
これによって、高調波成分における鉄損Wを次式によって求めることができる。
【0061】
【数15】
【0062】
このようにして、高調波成分を含む磁束密度ベクトルに基づいた鉄損Wを求めることができる。ここで求めたWは、高調波成分を考慮しない場合に比べて約30%大きな鉄損値となり、実機の鉄損との差も低減することから、より正確な電磁場解析が実現できる。
【0063】
図11に示すように、曲率K毎にH−B曲線およびW−B曲線を作成し、磁束密度Bおよび鉄損Wのデータベースを構築することにより、曲率Kを考慮した電磁場解析を行うことができる。このH−B曲線およびW−B曲線から、磁気抵抗率νおよび鉄損Wは、B,θ,ρの関数として、次式で与えられる。
【0064】
【数16】
【0065】
上記計算処理のルーチンを図12に示す。図12において、1は領域分割手段、2は磁束密度決定手段、3は磁束密度及び鉄損のデータベース、4は補間内挿計算手段、5は微小領域内鉄損計算手段、6は鉄損総和手段、7は曲率分布算出手段である。曲率分布算出手段7は、コンピュータを用いた構造解析による解析結果や三次元形状測定装置の測定結果により微小領域内aiの曲率分布Kiを算出する。データベース3はH−B曲線を表すデータベース、及びW−B曲線を表すデータベースより成り立つ。H−B曲線は周波数に依存し、磁気抵抗率νは周波数の増加に伴って増加する。従って、マックスウェル方程式では、νによって鋼板の周波数特性を含ませている。領域分割手段1は鉄鋼材料の対象領域全体を複数の微小領域(i=1〜n)に分割する。鉄鋼材料が一様な厚さの平面状板材であれば、微小領域の形状を三角形とすることができ、この場合には有限要素法による計算の適用が容易となる。磁束密度ベクトル決定手段2には、曲率分布算出手段7から算出された曲率分布データを入力し、一方では領域分割手段1からの領域分割の結果、及びデータベース3からの(H−B曲線)のデータを入力する。磁束密度ベクトル決定手段2では、分割された領域ごとにニュートンラプソン法を用い、マックスウェル方程式との連成により、テイラー展開から磁束密度ベクトルの収束計算を実行する。
【0066】
即ち、領域分割手段1からの入力データを基にして、磁束密度ベクトル決定手段は微小領域aiで、鋼板の圧延方向RDと磁束密度の方向との間の角度θ、及び磁束密度の大きさBとを決定して、補間内挿計算手段4に与える。ここで、圧延方向の代わりに、予め与えられた任意の方向とすることもできる。磁束密度に対する磁界(H−B曲線)及び鉄損(W−B曲線)のデータベース3は、曲率Kの値をパラメータとして表したデータ別に複数葉に分けられたデータの表を格納し、各データの表には角度θに対するH−B曲線及びW−B曲線を格納している。データベース3に格納されているH−B曲線及びW−B曲線は、それぞれサンプル化して与えられた角度θと曲率Kに対応して、磁束密度と磁界、及び磁束密度と鉄損の関係を与えるものである。これらのデータは、予め測定されたデータをサンプル化された角度θと曲率Kの関数の形で保持している。
【0067】
従って、任意の角度θと曲率Kに対する磁束密度Bと損失Wの関係は、補間内挿計算手段4により補間内挿法として、ニュートンラプソン法を使い、繰り返し計算を実行して求める。鉄鋼材料の鉄損データ(W−B曲線)は、データベースから補間内挿計算手段4に提供される。次にデータベース3からのW−B曲線のデータを利用し、補間内挿計算手段4によって求められた微小領域ai内の磁束密度Bと角度θを使い、微小領域鉄損計算手段5によって鉄損wiを計算する。そこで、鉄損総和手段6は、微小領域ai内の鉄損wiの総和Wを求める。このようにして、異方性をもった鉄鋼材料の鉄損が具体的に数値計算によって求められる。
【0068】
図13では、具体的実測に基づいたH−B曲線データについて、異方性を考慮した鉄鋼材料の鉄損計算ルーチンを示した。このルーチンによって実行されるマックスウェル方程式の解は、理論モデルにより次のように精度よく求められる。すなわち、磁気ベクトルポテンシャルAを用いて電流密度J0を表すと、次式が得られる。
【0069】
【数17】
【0070】
νをテンソル表示し、二次元場の式で表すと次式が得られる。
【0071】
【数18】
【0072】
これに対応した汎関数χは次式で与えられる。
【0073】
【数19】
【0074】
これを要素eにおける値で表し、要素eを構成する節点ieのポテンシャルで偏微分すると、
【0075】
【数20】
【0076】
が得られる。そこで、
【0077】
【数21】
【0078】
が有限要素法で解くべき方程式である。
【0079】
ここで、nuが未知節点の総数である。そこで、ニュートンラプソン法を適用するために書き直すと、次式が得られる。
【0080】
【数22】
【0081】
このマトリックスを解くことにより得られたk回目のδAjが求められれば、k+1回目の反復で得られる節点iでのポテンシャルの近似解
【0082】
【数23】
【0083】
が得られる。ただし、
【0084】
【数24】
【0085】
さらに解析を行うと、次式が得られる。
【0086】
【数25】
【0087】
上記の解析に従えば、テンソルで表された磁気抵抗率が解析モデルによる数値計算により求められる。具体的計算ルーチンの一例を図13に示す。図13に示したように、最初のステップS61において時刻を0にセットする。
【0088】
次に、S62で、
【0089】
【数26】
【0090】
の初期値を設定する。ここで、t=0のときは、
【0091】
【数27】
【0092】
t>0のときは、
【0093】
【数28】
【0094】
と、一つ前の時刻の値を用いる。次にS63で、式26に基づき、剛性マトリックス[K]、荷重ベクトル[F]の計算を行う。これにより、S64で、式29に基づき、δA(k,t)を計算する。これにより、S65で、k+1でのベクトルポテンシャル
【0095】
【数29】
【0096】
を計算する。
【0097】
S66で、ベクトルポテンシャルに基づき、式3より、磁束密度ベクトルBを求める。そこで、S67で、収束判定
【0098】
【数30】
【0099】
を行い、所定の収束性を満たしておれば、S71に進む。もし、収束性が不十分であれば、S68に進む。S68で、磁束密度ベクトルを、
【0100】
【数31】
【0101】
で、最大値とある方向の角度に分解する。そこで、S69で、予めデータベースにある異方性を考慮したH−B曲線より、その角度θ(k,t)における磁束密度と磁界の関係が分かり、それにより、磁気抵抗率νを求めることができる。ここでのH−B曲線は、その微小領域における曲率におけるH−B曲線であり、この曲率を考慮したH−B曲線を用いることで、異方性および曲率を考慮した磁束密度を得ることができる。
【0102】
これにより、
【0103】
【数32】
【0104】
の計算を行うことができる。これらのちは、S63で用いる剛性マトリックス、荷重ベクトルを求めるのに使用される。S71でkをひとつ進める。そして、S63に進み、計算を続ける。S63−S70までの計算は、すべての領域について計算を行う。この収束計算を行うことで、各微小領域において曲率を考慮したH−B曲線上にのった磁束密度ベクトルを計算することができ、曲率を考慮した計算ができる。S67で収束していたら、S71で時刻を一つ進めて、S62に向かう。この計算は、複数周期実行し、1周期分収束するまで計算を続ける。収束したら、S72に進み、鉄損の計算を行う。図13よりマックスウェル方程式の数値解から曲率を考慮した磁気抵抗率νがテンソル量として精度よく与えられているので、H−B曲線が精度よく解析される。
【0105】
次に、磁界増分係数αおよび鉄損増分係数βを用いてH−B曲線、W−B曲線を求める方法について述べる。そもそも、すべての鋼材、異方性を考慮した場合に対してH−B曲線、W−B曲線を求めることは通常大変なことであり、何らかの簡易計算が必要となる。ここでは、基本となる、H−B曲線、W−B曲線があった場合に、他の鋼材、異方性を考慮した特性を得る方法について述べる。
【0106】
図14は、実験から求められる磁束密度をパラメータとした曲率と磁界の関係図の一例である。磁束密度B1からB4に対して、曲率を0(平坦)から20[1/m]まで変更したときの磁界特性を示している。このとき、磁束密度Bをパラメータとして磁界増分係数α(=HK/HK=∞)(ここで、HK=0:曲率0のときの磁界、HK:曲率Kのときの磁界)が、図14上の各点に対して求めることができる。このデータに対し、曲率をパラメータにして、磁束密度−磁界増分係数αの関係を示したのが図15である。実測結果を尊重し、磁界増分係数αの磁束密度に対する関係は、測定生データを用いた。磁界増分係数αの場合、特性上B=0の近傍では、α=1となり、磁束密度Bが大きくなるにつれてαは大きくなるが、磁束密度が飽和する2T近傍では再びα=1となる。ここで、磁界増分係数αの磁束密度に対する関係式は生データを用いたが、測定生データのバラツキの影響を排除する意味合いから、測定生データの代わりに回帰式を用いても構わない。
【0107】
図16は、実験から求められる磁束密度をパラメータとした曲率と鉄損との関係図の一例である。磁束密度B1からB4に対して、曲率を0から20[1/m]までかけた時の鉄損特性を示している。このとき、磁束密度Bをパラメータにして鉄損増分係数β(=WK/WK=0)(ここで、WK=0:曲率0のときの鉄損、WK:曲率Kのときの鉄損)が、図16上の各点に対して求めることができる。これのデータに対して、曲率をパラメータとして磁束密度−鉄損増分係数βの関係にしたのが、図17である。
【0108】
ここで、鉄損増分係数βの磁束密度に対する関係は、測定生データを用いたが、測定生データのバラツキの影響を排除する意味合いから、測定生データの代わりに回帰式を用いても構わない。
【0109】
このとき、図18の点線のように、曲率0のときのH−B曲線が与えられたときに、点線上のある磁束密度に対して下記の式で、曲率Kとなったときの磁界Hρを求めることができる。
【0110】
【数33】
【0111】
これを、曲率0のときの、H−B曲線のすべての磁束密度に対して計算すれば、曲率ありのときのH−B曲線を求めることができる。これを図18上では、実線で示している。
【0112】
一方、図19の点線のように、曲率0のW−B曲線が与えられたときに、点線上のある磁束密度に対して下記の式で、曲率Kとなったときの鉄損Wρを求めることができる。
【0113】
【数34】
【0114】
これを、曲率0のときのW−B曲線のすべての磁束密度に対して計算すれば、曲率ありのときのW−B曲線を求めることができる。これを図19では、実線で示している。
【0115】
図20は、本発明のハードウエア構成を例示する図である。本発明における電磁場解析システムを実現するハードウエアは、スタンドアローン式のコンピュータと、十分な記憶容量を有するハードディスクでもよいが、多くのユーザが使用できる環境を整えるためには、図20のようなネットワークコンピュータと専用サーバから構成される方が好ましい。この構成により、ユーザ自身のコンピュータ端末からネットワークコンピュータに格納された本発明の電磁場解析システムに接続し、電磁場解析に必要な解析条件を入力することにより、ネットワークを通じて、本発明の電磁場解析システムが行った電磁場解析結果として、鉄損の評価情報を受け取ることができる。
【実施例】
【0116】
図23に示す外径φ115mm、内径φ60mm、積厚50mmのステータコア20に対して、外径φ59mm、積厚50mmの4極のロータを挿入したモータに対して、励磁電流500ATrms/相、50Hzの三相交流を励磁し、供試材は35A300とした条件での電磁場解析により、本発明の有効性を確認した。
【0117】
まず、比較例1として、ステータコア20の外部にSS400製の焼嵌めシェル(図示せず)による焼嵌めを行わないコアに対して電磁場解析を実行し、鉄損値を算出した。鉄損値は、電磁場解析の結果得られた磁束密度分布に対して、曲率0のW−B曲線と参照することにより求めた。
【0118】
次に、実施例1として、ステータコア20の外部に厚さ3mmのSS400製の焼嵌めシェル(図示せず)を常温でのステータコア外径と常温での焼嵌めシェル内径の差で示される焼きばめ代を直径で200μmとした条件で焼嵌めによりステータコアを固定した状況を有限要素法による応力解析を実行し、ステータコア20の弾性変形を算出した。その曲率半径分布を図24に示す。この曲率半径分布に従い、図4(a)に示すH−B曲線を用い、電磁場解析により磁束密度分布を算出し、その後、図4(b)に示すW−B曲線を用いて鉄損値を算出した。
【0119】
比較例1と実施例1の鉄損を比較例1の鉄損により正規化したものを図25に示す。焼嵌めによる弾性変形による曲率の変化を考慮することで、鉄損値が5%増加することが分かることから本発明の有効性を示すことができた。
【図面の簡単な説明】
【0120】
【図1】従来技術に基づいて実行される鉄損の数値計算ルーチンの流れを示すブロック図である。
【図2】鉄鋼材料の対象とする全領域を格子状の複数の微小領域に分割した模様を示す説明図である。
【図3】磁束密度の異方性を示す説明図である。
【図4】曲率に対する依存性を複数葉のグラフで与えることにより、H−B曲線及びW−B曲線の角度θ依存性及び曲率K依存性を示すグラフ図である。
【図5】磁束密度と磁界のベクトルの位相差を示す図である。
【図6】本発明に用いるH−B曲線を示す図である。
【図7】本発明に用いるB−θBH曲線を示す図である。
【図8】磁束密度Bの時間的変化を示す図である。
【図9】本発明に用いるH−B曲線を示す図である。
【図10】本発明に用いるW−B曲線を示す図である。
【図11】本発明に用いるH−B曲線ならびにW−B曲線を示す図である。
【図12】異方性を考慮した鉄鋼材料の鉄損計算ルーチンを示すブロック図である。
【図13】マックスウェル方程式により磁気抵抗率νをテンソルとして精度よく求めるルーチンを示す説明図である。
【図14】実験から求められる磁束密度をパラメータとした曲率と磁界の関係図の一例である。
【図15】図14より求めた曲率をパラメータとした磁束密度と磁界増分係数との関係図の一例である。
【図16】実験から求められる磁束密度をパラメータとした曲率と鉄損との関係図の一例である。
【図17】図16より求めた曲率をパラメータとした磁束密度と鉄損増分係数との関係図の一例である。
【図18】図15より求めた曲率ありの場合におけるH−B曲線の一例を示す図である。
【図19】図17より求めた曲率ありの場合におけるW−B曲線の一例を示す図である。
【図20】本発明のハードウエア構成を例示する図である。
【図21】鉄芯における曲率の磁束密度への影響例を示した説明図である。
【図22】変圧器鉄芯におけるコアの積み上げ状態を示した説明図である。
【図23】実施例で使用したステータコア概形図である。
【図24】本発明の実施例で使用した焼嵌め後の曲率半径分布を示す図である。
【図25】本発明による効果を示す鉄損比を示す図である。
【符号の説明】
【0121】
1 領域分割手段、
2 磁束密度決定手段、
3 磁束密度及び鉄損のデータベース、
4 補間内挿計算手段、
5 微小領域内鉄損計算手段、
6 鉄損総和手段
7 曲率分布算出手段
11 凹型鉄芯
12 円筒リング
20 ステータコア
【技術分野】
【0001】
本発明は、鋼板をはじめとする軟磁性材料の異方性および弾性曲げ変形による磁気特性劣化を考慮した変形・電磁場連成解析システムに関し、特に軟磁性材料の弾性曲げ変形に関する磁気特性の劣化を考慮した電気機器の鉄損を予測する方法に関する。
【背景技術】
【0002】
環境問題に対する意識の高まりにより、受配電に使用される変圧器ならびに空調機やハイブリット自動車に使用されるモータに対し、鉄損や銅損に代表される損失の低減の要求が高まっている。それに対し、軟磁性材料として実用上多く使用されている電磁鋼板は、鋼板圧延方向からの角度、励磁方向への応力および、弾性変形による曲げにより、励磁特性や鉄損特性が劣化する。このため、変圧器やモータ等の電気機器の性能を高精度に予測するためには、このような使用している材料の異方性や励磁方向への応力、そして、弾性変形による磁気特性の劣化を考慮することが欠かせない。
【0003】
従来は、励磁方向への応力や弾性変形による曲げを考慮した磁気特性を評価する場合、実機試作材料の実測を行いエプスタイン試験法による鉄損特性と単板磁気測定法(SST法)による鉄損特性との差を励磁方向への応力や弾性変形による応力に起因する磁気特性の劣化分として算出していた。この方法では、励磁方向への応力の影響による磁気特性の劣化分と弾性変形による応力に起因する磁気特性の劣化分を分離することができなかった。
【0004】
さらに、従来の弾性変形による曲げの影響評価方法は実機試作を前提としていたため、時間的、金銭的なコストがかかることに加え、形状を連続的に変更し試作することが困難であるため、電気機器設計時に損失最小となる最適形状を求めることが非常に困難であった。
【0005】
一方、従来よりマックスウエル方程式による電磁場解析法は、鋼板等の軟磁性材料の鉄損評価に用いられてきた。以下ここでは、軟磁性材料として実用上多く使用されている鋼板でもって代表的に説明することにする。非特許文献2に示されるように、電磁場解析には、コンピュータが使用され、鋼板の形状、計算のために分割された微小領域の大きさ、磁界Hに対する磁束密度B、周波数等は、コンピュータによる計算の解を求めるための物理量パラメータ条件として採用されている。即ち、これらの条件を考慮に入れて、マックスウエル方程式の計算機解が得られる。
【0006】
マックスウエル方程式の計算機解に、磁束密度Bに対応して測定された鉄損Wのデータ(W−B曲線)を考慮に加えて鉄損が求められている。例えば図1は、従来技術に基づいて実行される鉄損の数値計算ルーチンの流れを示すブロック図である。まず、鋼板の形状、微小領域分割、周波数等をパラメータとしてH−B曲線など、解を求めるための条件を与え、計算によってマックスウエル方程式の数値解を求める。磁束密度Bに対応して測定された、微小領域における鋼材の鉄損Wのデータ(W−B曲線)を上の数値解に与え、鋼板全体の損失Wを算出する。ここで、鋼板に作用する弾性曲げ変形の影響は無視され、鋼板に曲げ変形が加わっても磁気特性は変わらないと仮定している。
【0007】
特許文献1には、データ入力部と連成解析部と結果出力部より構成される磁気デバイスの磁気特性解析方法について記載されている。ここで、データ入力部より入力されるデータは、磁気デバイスを構成される材料の特性に関わるデータ、磁気デバイスの構造に関わるデータ、構造解析の境界条件に関わるデータ、磁場解析の境界条件に関わるデータである。連成解析部において、データ入力部から入力された前記境界条件に関わるデータに基づき区分された複数部分ごとの応力分布を求め、前記磁気デバイスの電磁場解析の境界条件に関わるデータ及び前記応力分布に基づき前記複数部分ごとの磁気特性を求め、前記複数部分ごとの磁気特性に基づき磁気デバイス全体の磁気特性を求めるものである。ここでは、磁気デバイスで用いられる磁性材料や非磁性材料の荷重や熱応力による歪の影響は考慮されているものの、磁性材料に作用する弾性曲げの影響は無視され、鋼板に曲げ変形が加わっても磁気特性は変わらないと仮定している。
【0008】
さらに、特許文献2には、磁性体の予め与えられた方向と磁束密度の方向との間の角度θおよび応力σを異方性のパラメータとして磁束密度と磁界とを関係づける解析式およびデータに基づくH−B曲線、および、磁束密度と鉄損とを関係づける解析式およびデータに基づくW−B曲線を用いて磁束密度ベクトル決定手段と鉄損総和手段により、磁気特性を求める方法について記載されている。ここでも同様に、磁性体である電磁鋼板の応力の影響は考慮されているものの、磁性材料に作用する弾性曲げの影響は無視され、鋼板に曲げ変形が加わっても磁気特性が変わらないと仮定している。
【0009】
ここで、マックスウエル方程式の解法について簡単に説明する。公知の多くの文献から明らかなように、マックスウエル方程式は次式で与えられる。
【0010】
【数1】
【0011】
ここで、B、H、D、E、Jはそれぞれ磁束密度、磁界、電束密度、電界、電流密度である。また、ρは電荷密度である。B、H、D、E、Jの間には、次の関係がある。
【0012】
【数2】
【0013】
ここで、μ、ε、σはそれぞれ、透磁率、誘電率、導電率である。
【0014】
一方、非特許文献2によれば、電磁界に関する解析が詳細に記載されている。同文献よれば、∂D/∂tは無視されている。磁束の発散は常に零であるので、連続であり、磁気ベクトルポテンシャルAが次式によって与えられている。
【0015】
【数3】
【0016】
これらの式から
【0017】
【数4】
【0018】
が得られている。従って、
【0019】
【数5】
【0020】
が得られる。ここで、−gradφ=E、J0は外部からの強制電流密度、Jeはうず電流密度、テンソル量で与えられる磁気抵抗率[ν]は、[ν]=1/[μ]である。(5)式は、ガラーキン法(Galerkin Method)により2次元的、及び3次元的に解かれる。実際には、透磁率は、一般的に弾性曲げ変形の影響を受ける。しかし、従来の鉄損評価法では、透磁率が弾性曲げ変形の影響を受けないと仮定して、数値解析が行われ、解が求められていた。
【0021】
また、H−B曲線は曲率に対して影響がなく、一定であると仮定している。このようにして、与えられた磁界強度Hに対する微小領域ごとに求められた、H−B曲線から、各微小領域における磁束密度の分布が求められる。一方、磁束密度Bと鉄損Wの関係は曲率Kに対しては一定としており、W−B曲線の形でデータベースに格納されている。従って、上記数値計算によって求められた磁束密度分布をデータベースに格納されたW−B曲線に適用すれば、与えられた磁界強度Hに対する鋼板全体の鉄損が数値計算によって求められる。上記従来技術では、鉄損を求めるための電磁場解析において、透磁率は等方性とし、また、弾性曲げ変形の影響を無視した透磁率でもって計算を行っているため、実際の鉄鋼材料において存在する透磁率の異方性や曲率の影響が無視され、そのために実測値の計算値からの乖離が無視できず、鉄損を十分に評価できないという欠点があった。
【特許文献1】特開平8−249621号公報
【特許文献2】特許3676761号公報
【非特許文献1】Nippon Steel Corporation:「Technical Data on Domain Refined ORIENTCORE・HI−B ZDKH Second Edition」
【非特許文献2】中田、高橋:「電気工学の有限要素法」(森北出版、1982)
【発明の開示】
【発明が解決しようとする課題】
【0022】
そこで、本発明は、このような従来技術の問題点を解決し、従来の解決方法を踏襲しながら、軟磁性材料の物性値パラメータに対する、異方性と弾性曲げの影響を考慮に入れた電磁場解析システムを提供して、正確な鉄損評価を実現することを課題とする。
【課題を解決するための手段】
【0023】
上記目的を達成するため、本発明による電磁場解析システムは、軟磁性材料の電磁場解析対象領域を複数の微小領域に分割する領域分割手段と、該微小領域における磁性体の予め定められた方向と磁束密度Bの方向との間の成す角度θおよび軟磁性材料の弾性曲げ曲率Kを異方性のパラメータとして、磁束密度Bと磁界Hとを関係付けるH−B曲線および磁束密度Bと鉄損Wとを関係付けるW−B曲線を格納するデータベースと、該データベースに格納されているH−B曲線を基にして、前記微小領域におけるマックスウエル方程式および弾性曲げ曲率Kに基づいて、前記角度θおよび前記磁束密度Bの大きさを決定する磁束密度ベクトル決定手段と、前記データベースに格納されているW−B曲線を基にして、前記微小領域の鉄損Wを計算する鉄損計算手段と、前記微小領域の鉄損Wの総和を求める鉄損総和手段とを有することを特徴とする。
【0024】
また、前記磁束密度Bおよび磁界Hは、該磁束密度Bと磁界Hとの位相差θBHを考慮して関係付けられることを特徴とする。
【0025】
また、前記磁束密度Bは、時間高調波成分を含む磁束密度であることを特徴とする。
【0026】
さらに、前記磁束密度ベクトル決定手段は、前記弾性曲げ曲率Kを算出する弾性曲げ曲率分布算出手段を備えることを特徴とする。
【0027】
さらに、前記軟磁性材料は、電気機器を構成する軟磁性材料であり、前記弾性曲げ曲率Kは、該電気機器の製作時あるいは自重により発生する弾性曲げ曲率であることを特徴とする。
【0028】
さらに、前記磁束密度ベクトル決定手段は、前記微小領域における磁束密度ベクトルを、その大きさBmaxと、前記の磁性体の予め定められた方向とのなす角度θとに分解し、前記H−B曲線から、前記角度θにおけるH−B曲線を導出し、前記微小領域における磁束密度およびその前記角度θ、磁界とが、その曲線上に存在するように、収束計算を行うことを特徴とする。
【0029】
さらに、前記H−B曲線およびW−B曲線は、前記弾性曲げ曲率Kをパラメータとして、磁束密度と磁界の関係より磁界増分係数α(=HK/HK=0)、磁束密度と鉄損値の関係より鉄損増分係数β(WK/WK=0)を求め、前記弾性曲げ曲率K=0のときのH−B曲線およびW−B曲線に前記磁界増分係数αと鉄損増分係数βとを乗じて得られるH−B曲線およびW−B曲線であることを特徴とする。
【発明の効果】
【0030】
本発明は、圧延による軟磁性材料の異方性と弾性曲げ曲率のH−B曲線およびW−B曲線への影響を考慮した鉄損の評価手段を見出し、電磁場解析システムとして具現化したものである。すなわち、本発明は鉄鋼材料の全領域を微小領域に分割し、分割された各微小領域に対して予め与えられた方向と磁束密度との成す角度と、例えば加工や自重の影響で生じた弾性変形による曲率とをパラメータとして考慮に入れた鉄損データをデータベースから求めて計算し、更に各微小領域で求められた鉄損データの総和をとることにより、異方性が弾性曲げの影響の大きい材料であっても、磁束密度の異方性・弾性曲げ影響や磁界に対する磁束密度の非直線性を意識しないで鉄損を評価計算することができる。また、磁束密度Bと磁界Hとの位相差θBH、時間高調波を考慮した電磁場解析を行うことが出来るなど産業上有用な、著しい効果を奏する。
【発明を実施するための最良の形態】
【0031】
以下、図面を参照して本発明を詳細に説明する。
【0032】
鉄鋼材料には鋼材の予め定められた方向、例えば圧延方向とそれに直角の方向とでは透磁率に異方性がある。このような鉄鋼材料の透磁率の異方性を考慮に入れれば、鉄損を精度よく評価することが可能となる。まず、電磁場解析の対象とする全領域を複数の微小領域に分割する。各微小領域の内部で透磁率が等方性を満足する場合は、各微小領域間で透磁率μ及び磁気抵抗率νや、予め定められた方向、例えば圧延方向と磁束密度Bの方向との間の成す角度θに差があっても容易に計算をすることができる。透磁率に異方性がある場合は、各微小領域にそれぞれ角度θとそれに対する透磁率を与えれば、コンピュータを利用してマックスウェル方程式の解を求めることができる。
【0033】
また、鉄鋼材料には、モータ製造工程の一つである焼嵌めや自重により、鋼板が弾性変形している場合が多く、弾性変形により、透磁率をはじめとした磁気特性が大きく変化する。このような鉄鋼材料の透磁率に対する弾性変形の影響を考慮に入れれば、鉄損を精度よく評価することが可能となる。この場合、電磁場解析の対象とする全領域を複数の微小領域に分割し、各微小領域の内部で弾性曲げ曲率が一定であるとすれば、その微小区間内では同一の磁気特性、透磁率を用いても差し支えない。このように、各微小領域にそれぞれと弾性曲げ曲率とそれに対する透磁率の関係を与えれば、コンピュータを利用してマックスウェル方程式の解を求めることができる。
【0034】
図2は、鉄鋼材料の対象とする全領域を格子状の複数の微小領域に分割した模様を示す説明図である。図2において、鉄鋼材料は、1,2,3……,i,i+1,i+2,……,j,j+1,j+2,……,nの微小領域に分割してある。例えば、i番目の微小領域の磁束密度がΔBiであるとし、鉄損がΔWiであるとする。
【0035】
図2において分割された各微小領域の内部においては、透磁率μ及び磁気抵抗率νは一定の値であるとする。このようにすれば、有限要素法において、各微小領域間の境界では不連続であっても、領域内では一様なパラメータをもっていると考えることができる。従って、異方性・弾性曲げの影響や磁界に対する磁束密度の非直線性を有する鉄鋼材料の鉄損の計算においても、予め計算された各微小領域の鉄損の総和を求めることによって全体鉄損を容易に計算することができる。図2では、全体の鉄損W(watts)は
【0036】
【数6】
【0037】
で与えられる。ここで、Δwiは微小領域内の鉄損である。また、ΔBiの方向θと大きさΔBiは、各微小領域によってそれぞれ異なった値をとる。鉄損Wは磁束密度Bが大きい程、大きな値をとるが、必ずしも直線関係になるわけではない。上記領域分割は有限要素法が適用されることを前提にして実施したものであるが、差分法或いはその他、類似の計算方式に適用可能であることは云うまでもない。
【0038】
図3は、磁束密度の異方性を示す説明図である。以下弾性曲げ曲率だけでなく磁束密度の異方性をも考慮した場合についても考えてみる。RDは鋼板の圧延方向を(rolling direction)、TDはそれに直角な方向(transversal direction)である。磁束密度Bの方向と鋼板の圧延方向RDとの間の角度をθとすれば、H−B曲線及びW−B曲線のθ依存性は、例えば図4(a)、(b)に示すような形状で与えられる。ここで、鋼板の圧延方向の代わりに、任意の予め与えられた方向とすることも可能である。図4(a)、(b)はいずれも角度θの依存性を表したものであるが、複数葉のデータは更に曲率Kに対する依存性も示している。
【0039】
図4(a)、(b)では、θ=0°,15°,30°……のデータがH−B関数及びW−B関数として与えられているが、その他の角度θではデータが関数として与えられていない。この場合、θ=20°において求める必要のある磁束密度と磁界との関係は、θ=15°及びθ=30°におけるデータから求められた近似関数について、以下で詳細に述べる補間内挿して求めたり、ニュートンラプソン法によるテイラー展開からマックスウェル方程式との連成により繰り返し計算して求められる。一方、磁束密度が与えられたときの鉄損は、θ=15°及びθ=30°におけるデータから補間内挿して求められる。ここで、曲率Kの逆数の変化が十分に無視できるほど小さい場合は、曲率は一定であると仮定してかまわない。もし、曲率Kの逆数が十分に大きく変化する場合には、それぞれの曲率に対応したH−B曲線及びW−B曲線を選択して繰り返し計算を実行する。ここで、H−B曲線及びW−B曲線は、データに基づいて求める方法について記したが、理論的に得られた解析式およびデータから得られた回帰式といった解析式を用いても構わない。
【0040】
次に図4(a)、(b)に含まれている曲率Kに対する依存性に関して少し詳細に説明する。例えば、図4(a)、(b)では、0[1/m]<K<20[1/m]の範囲のデータを取り扱い、サンプル化された曲率Kの値に対するデータとして示している。従って、データの表に直接与えられていない曲率に対する磁束密度や鉄損は、角度θの場合と同様にニュートンラプソン法によるテイラー展開から補間内挿して求められる。すなわち、(Bi,θi,ρ)によって微小領域iにおける鉄損wiが求められる。
【0041】
例えば、図4(b)に示すK=20[1/m]のデータにおいて、θ=15°及び30°の曲線、及びB=1.0T(Tesla)及び1.2Tによって切断され、点ABCDによって囲まれた領域を仮定し、点P(θ=20°、B=1.1T)における鉄損を求めるものとする。まず、θ=15°及び30°の曲線から内挿近似によって級数展開を行い、θ=20°のW−B曲線を求める。次にB=1.0T及び1.2Tの直線(曲線でもよい。)によって上記θ=20°のW−B曲線を切断、B軸上の内挿によってB=1.1Tにおける点Pを確定する。確定された点Pの縦軸から鉄損Wを読取ることができる。実際の処理では、図4(b)のグラフはデータベース内に格納されているので、コンピュータのソフトウエア処理によって内挿計算を行い、鉄損を求めることができる。
【0042】
曲率Kの影響によってH−B曲線が変化することは、図4によって示したとおりである。例えば、コンプレッサ等で使用される凹型鉄芯は鋼板をせん断して成形され、70℃から90℃くらいに加熱された軟鋼製の円筒リングに挿入されることで焼嵌めにより固定され、電気機器である回転電機に使用される。円筒リングが常温まで冷却収縮することにより、凹型鉄芯は、内部に圧縮応力を受けるとともに、弾性変形する。この弾性変形量を曲率Kにて表現する。図21は、凹型鉄芯における弾性変形による曲率の磁束密度への影響例を示した説明図である。図21において、xは鋼板のある方向としての圧延方向(RD)、Bは磁束の接線方向であり磁束密度を表す、θはx方向と磁束密度Bとのなす角度を示す。空間上の指定された位置を決定すると、磁束密度ベクトルBが与えられ、曲率Kの値が決定される。すなわち、指定された空間位置において、H−B曲線より、(θ、K)に対応する磁束密度Bが求められ、これによって、W−B曲線に(B、θ、K)を適用することにより鉄損Wが算出される。空間位置に対応する曲率は、構造力学の理論に従って有限要素法により算出することが可能である。
【0043】
さらに、電気機器の一つである変圧器においては、方向性電磁鋼板をせん断加工後、図22のように組み上げてコアを形成する。その際、コアの自重により鋼板が弾性変形する。弾性変形した際の取り扱いの方法は、前記の焼嵌め固定されたコアと同様である。
【0044】
図5は、磁束密度と磁界のベクトルの位相差を示す図である。図5において、磁界Hと磁束密度Bはベクトルで、その間に位相差θBHがある。このとき、磁気抵抗率νは簡略的に次式のようにモデル化できる。
【0045】
【数7】
【0046】
このとき、磁束密度ベクトルおよび磁界ベクトルは次式で表される。
【0047】
【数8】
【0048】
従って、磁気抵抗率νは、Bおよびθの関数として
【0049】
【数9】
【0050】
であらわされ、図6のようなH−B曲線から求めることができる。また、磁束密度Bと磁界Hとの位相差θBHも、Bおよびθの関数として
【0051】
【数10】
【0052】
で表されるので、図7のようなB−θBH曲線をデータベースとして用いることによりθBHを考慮した磁束密度ベクトルを決定することができる。
【0053】
図8は、磁束密度Bの時間的変化を示す図である。図8において、横軸は時間、縦軸は磁束密度を示している。例えば、磁界Hが正弦波であっても、磁束密度Bは正弦波とはならず時間高調波成分を含むため図8のような歪が生じているので、正確な電磁場解析を行うためにはこの時間高調波を考慮する必要がある。そこで、時間高調波成分を含んだ磁束密度Bをフーリエ級数展開し、その周波数ごとに図9のようなH−B曲線と図10のようなW−B曲線をデータベースとして用いることによって時間高調波成分を考慮した電磁場解析を行うことができる。このH−B曲線およびW−B曲線から磁気抵抗率νおよび鉄損Wは、B,θおよび周波数fの関数として、次式で与えられる。
【0054】
【数11】
【0055】
また、時間高調波成分を含んだ磁束密度Bをフーリエ級数展開すると、Bのx成分、y成分はそれぞれ次式で表される。
【0056】
【数12】
【0057】
これにより、高調波における第i番目の磁束密度ベクトルおよび、その位相差は次式によって定義できる。
【0058】
【数13】
【0059】
【数14】
【0060】
これによって、高調波成分における鉄損Wを次式によって求めることができる。
【0061】
【数15】
【0062】
このようにして、高調波成分を含む磁束密度ベクトルに基づいた鉄損Wを求めることができる。ここで求めたWは、高調波成分を考慮しない場合に比べて約30%大きな鉄損値となり、実機の鉄損との差も低減することから、より正確な電磁場解析が実現できる。
【0063】
図11に示すように、曲率K毎にH−B曲線およびW−B曲線を作成し、磁束密度Bおよび鉄損Wのデータベースを構築することにより、曲率Kを考慮した電磁場解析を行うことができる。このH−B曲線およびW−B曲線から、磁気抵抗率νおよび鉄損Wは、B,θ,ρの関数として、次式で与えられる。
【0064】
【数16】
【0065】
上記計算処理のルーチンを図12に示す。図12において、1は領域分割手段、2は磁束密度決定手段、3は磁束密度及び鉄損のデータベース、4は補間内挿計算手段、5は微小領域内鉄損計算手段、6は鉄損総和手段、7は曲率分布算出手段である。曲率分布算出手段7は、コンピュータを用いた構造解析による解析結果や三次元形状測定装置の測定結果により微小領域内aiの曲率分布Kiを算出する。データベース3はH−B曲線を表すデータベース、及びW−B曲線を表すデータベースより成り立つ。H−B曲線は周波数に依存し、磁気抵抗率νは周波数の増加に伴って増加する。従って、マックスウェル方程式では、νによって鋼板の周波数特性を含ませている。領域分割手段1は鉄鋼材料の対象領域全体を複数の微小領域(i=1〜n)に分割する。鉄鋼材料が一様な厚さの平面状板材であれば、微小領域の形状を三角形とすることができ、この場合には有限要素法による計算の適用が容易となる。磁束密度ベクトル決定手段2には、曲率分布算出手段7から算出された曲率分布データを入力し、一方では領域分割手段1からの領域分割の結果、及びデータベース3からの(H−B曲線)のデータを入力する。磁束密度ベクトル決定手段2では、分割された領域ごとにニュートンラプソン法を用い、マックスウェル方程式との連成により、テイラー展開から磁束密度ベクトルの収束計算を実行する。
【0066】
即ち、領域分割手段1からの入力データを基にして、磁束密度ベクトル決定手段は微小領域aiで、鋼板の圧延方向RDと磁束密度の方向との間の角度θ、及び磁束密度の大きさBとを決定して、補間内挿計算手段4に与える。ここで、圧延方向の代わりに、予め与えられた任意の方向とすることもできる。磁束密度に対する磁界(H−B曲線)及び鉄損(W−B曲線)のデータベース3は、曲率Kの値をパラメータとして表したデータ別に複数葉に分けられたデータの表を格納し、各データの表には角度θに対するH−B曲線及びW−B曲線を格納している。データベース3に格納されているH−B曲線及びW−B曲線は、それぞれサンプル化して与えられた角度θと曲率Kに対応して、磁束密度と磁界、及び磁束密度と鉄損の関係を与えるものである。これらのデータは、予め測定されたデータをサンプル化された角度θと曲率Kの関数の形で保持している。
【0067】
従って、任意の角度θと曲率Kに対する磁束密度Bと損失Wの関係は、補間内挿計算手段4により補間内挿法として、ニュートンラプソン法を使い、繰り返し計算を実行して求める。鉄鋼材料の鉄損データ(W−B曲線)は、データベースから補間内挿計算手段4に提供される。次にデータベース3からのW−B曲線のデータを利用し、補間内挿計算手段4によって求められた微小領域ai内の磁束密度Bと角度θを使い、微小領域鉄損計算手段5によって鉄損wiを計算する。そこで、鉄損総和手段6は、微小領域ai内の鉄損wiの総和Wを求める。このようにして、異方性をもった鉄鋼材料の鉄損が具体的に数値計算によって求められる。
【0068】
図13では、具体的実測に基づいたH−B曲線データについて、異方性を考慮した鉄鋼材料の鉄損計算ルーチンを示した。このルーチンによって実行されるマックスウェル方程式の解は、理論モデルにより次のように精度よく求められる。すなわち、磁気ベクトルポテンシャルAを用いて電流密度J0を表すと、次式が得られる。
【0069】
【数17】
【0070】
νをテンソル表示し、二次元場の式で表すと次式が得られる。
【0071】
【数18】
【0072】
これに対応した汎関数χは次式で与えられる。
【0073】
【数19】
【0074】
これを要素eにおける値で表し、要素eを構成する節点ieのポテンシャルで偏微分すると、
【0075】
【数20】
【0076】
が得られる。そこで、
【0077】
【数21】
【0078】
が有限要素法で解くべき方程式である。
【0079】
ここで、nuが未知節点の総数である。そこで、ニュートンラプソン法を適用するために書き直すと、次式が得られる。
【0080】
【数22】
【0081】
このマトリックスを解くことにより得られたk回目のδAjが求められれば、k+1回目の反復で得られる節点iでのポテンシャルの近似解
【0082】
【数23】
【0083】
が得られる。ただし、
【0084】
【数24】
【0085】
さらに解析を行うと、次式が得られる。
【0086】
【数25】
【0087】
上記の解析に従えば、テンソルで表された磁気抵抗率が解析モデルによる数値計算により求められる。具体的計算ルーチンの一例を図13に示す。図13に示したように、最初のステップS61において時刻を0にセットする。
【0088】
次に、S62で、
【0089】
【数26】
【0090】
の初期値を設定する。ここで、t=0のときは、
【0091】
【数27】
【0092】
t>0のときは、
【0093】
【数28】
【0094】
と、一つ前の時刻の値を用いる。次にS63で、式26に基づき、剛性マトリックス[K]、荷重ベクトル[F]の計算を行う。これにより、S64で、式29に基づき、δA(k,t)を計算する。これにより、S65で、k+1でのベクトルポテンシャル
【0095】
【数29】
【0096】
を計算する。
【0097】
S66で、ベクトルポテンシャルに基づき、式3より、磁束密度ベクトルBを求める。そこで、S67で、収束判定
【0098】
【数30】
【0099】
を行い、所定の収束性を満たしておれば、S71に進む。もし、収束性が不十分であれば、S68に進む。S68で、磁束密度ベクトルを、
【0100】
【数31】
【0101】
で、最大値とある方向の角度に分解する。そこで、S69で、予めデータベースにある異方性を考慮したH−B曲線より、その角度θ(k,t)における磁束密度と磁界の関係が分かり、それにより、磁気抵抗率νを求めることができる。ここでのH−B曲線は、その微小領域における曲率におけるH−B曲線であり、この曲率を考慮したH−B曲線を用いることで、異方性および曲率を考慮した磁束密度を得ることができる。
【0102】
これにより、
【0103】
【数32】
【0104】
の計算を行うことができる。これらのちは、S63で用いる剛性マトリックス、荷重ベクトルを求めるのに使用される。S71でkをひとつ進める。そして、S63に進み、計算を続ける。S63−S70までの計算は、すべての領域について計算を行う。この収束計算を行うことで、各微小領域において曲率を考慮したH−B曲線上にのった磁束密度ベクトルを計算することができ、曲率を考慮した計算ができる。S67で収束していたら、S71で時刻を一つ進めて、S62に向かう。この計算は、複数周期実行し、1周期分収束するまで計算を続ける。収束したら、S72に進み、鉄損の計算を行う。図13よりマックスウェル方程式の数値解から曲率を考慮した磁気抵抗率νがテンソル量として精度よく与えられているので、H−B曲線が精度よく解析される。
【0105】
次に、磁界増分係数αおよび鉄損増分係数βを用いてH−B曲線、W−B曲線を求める方法について述べる。そもそも、すべての鋼材、異方性を考慮した場合に対してH−B曲線、W−B曲線を求めることは通常大変なことであり、何らかの簡易計算が必要となる。ここでは、基本となる、H−B曲線、W−B曲線があった場合に、他の鋼材、異方性を考慮した特性を得る方法について述べる。
【0106】
図14は、実験から求められる磁束密度をパラメータとした曲率と磁界の関係図の一例である。磁束密度B1からB4に対して、曲率を0(平坦)から20[1/m]まで変更したときの磁界特性を示している。このとき、磁束密度Bをパラメータとして磁界増分係数α(=HK/HK=∞)(ここで、HK=0:曲率0のときの磁界、HK:曲率Kのときの磁界)が、図14上の各点に対して求めることができる。このデータに対し、曲率をパラメータにして、磁束密度−磁界増分係数αの関係を示したのが図15である。実測結果を尊重し、磁界増分係数αの磁束密度に対する関係は、測定生データを用いた。磁界増分係数αの場合、特性上B=0の近傍では、α=1となり、磁束密度Bが大きくなるにつれてαは大きくなるが、磁束密度が飽和する2T近傍では再びα=1となる。ここで、磁界増分係数αの磁束密度に対する関係式は生データを用いたが、測定生データのバラツキの影響を排除する意味合いから、測定生データの代わりに回帰式を用いても構わない。
【0107】
図16は、実験から求められる磁束密度をパラメータとした曲率と鉄損との関係図の一例である。磁束密度B1からB4に対して、曲率を0から20[1/m]までかけた時の鉄損特性を示している。このとき、磁束密度Bをパラメータにして鉄損増分係数β(=WK/WK=0)(ここで、WK=0:曲率0のときの鉄損、WK:曲率Kのときの鉄損)が、図16上の各点に対して求めることができる。これのデータに対して、曲率をパラメータとして磁束密度−鉄損増分係数βの関係にしたのが、図17である。
【0108】
ここで、鉄損増分係数βの磁束密度に対する関係は、測定生データを用いたが、測定生データのバラツキの影響を排除する意味合いから、測定生データの代わりに回帰式を用いても構わない。
【0109】
このとき、図18の点線のように、曲率0のときのH−B曲線が与えられたときに、点線上のある磁束密度に対して下記の式で、曲率Kとなったときの磁界Hρを求めることができる。
【0110】
【数33】
【0111】
これを、曲率0のときの、H−B曲線のすべての磁束密度に対して計算すれば、曲率ありのときのH−B曲線を求めることができる。これを図18上では、実線で示している。
【0112】
一方、図19の点線のように、曲率0のW−B曲線が与えられたときに、点線上のある磁束密度に対して下記の式で、曲率Kとなったときの鉄損Wρを求めることができる。
【0113】
【数34】
【0114】
これを、曲率0のときのW−B曲線のすべての磁束密度に対して計算すれば、曲率ありのときのW−B曲線を求めることができる。これを図19では、実線で示している。
【0115】
図20は、本発明のハードウエア構成を例示する図である。本発明における電磁場解析システムを実現するハードウエアは、スタンドアローン式のコンピュータと、十分な記憶容量を有するハードディスクでもよいが、多くのユーザが使用できる環境を整えるためには、図20のようなネットワークコンピュータと専用サーバから構成される方が好ましい。この構成により、ユーザ自身のコンピュータ端末からネットワークコンピュータに格納された本発明の電磁場解析システムに接続し、電磁場解析に必要な解析条件を入力することにより、ネットワークを通じて、本発明の電磁場解析システムが行った電磁場解析結果として、鉄損の評価情報を受け取ることができる。
【実施例】
【0116】
図23に示す外径φ115mm、内径φ60mm、積厚50mmのステータコア20に対して、外径φ59mm、積厚50mmの4極のロータを挿入したモータに対して、励磁電流500ATrms/相、50Hzの三相交流を励磁し、供試材は35A300とした条件での電磁場解析により、本発明の有効性を確認した。
【0117】
まず、比較例1として、ステータコア20の外部にSS400製の焼嵌めシェル(図示せず)による焼嵌めを行わないコアに対して電磁場解析を実行し、鉄損値を算出した。鉄損値は、電磁場解析の結果得られた磁束密度分布に対して、曲率0のW−B曲線と参照することにより求めた。
【0118】
次に、実施例1として、ステータコア20の外部に厚さ3mmのSS400製の焼嵌めシェル(図示せず)を常温でのステータコア外径と常温での焼嵌めシェル内径の差で示される焼きばめ代を直径で200μmとした条件で焼嵌めによりステータコアを固定した状況を有限要素法による応力解析を実行し、ステータコア20の弾性変形を算出した。その曲率半径分布を図24に示す。この曲率半径分布に従い、図4(a)に示すH−B曲線を用い、電磁場解析により磁束密度分布を算出し、その後、図4(b)に示すW−B曲線を用いて鉄損値を算出した。
【0119】
比較例1と実施例1の鉄損を比較例1の鉄損により正規化したものを図25に示す。焼嵌めによる弾性変形による曲率の変化を考慮することで、鉄損値が5%増加することが分かることから本発明の有効性を示すことができた。
【図面の簡単な説明】
【0120】
【図1】従来技術に基づいて実行される鉄損の数値計算ルーチンの流れを示すブロック図である。
【図2】鉄鋼材料の対象とする全領域を格子状の複数の微小領域に分割した模様を示す説明図である。
【図3】磁束密度の異方性を示す説明図である。
【図4】曲率に対する依存性を複数葉のグラフで与えることにより、H−B曲線及びW−B曲線の角度θ依存性及び曲率K依存性を示すグラフ図である。
【図5】磁束密度と磁界のベクトルの位相差を示す図である。
【図6】本発明に用いるH−B曲線を示す図である。
【図7】本発明に用いるB−θBH曲線を示す図である。
【図8】磁束密度Bの時間的変化を示す図である。
【図9】本発明に用いるH−B曲線を示す図である。
【図10】本発明に用いるW−B曲線を示す図である。
【図11】本発明に用いるH−B曲線ならびにW−B曲線を示す図である。
【図12】異方性を考慮した鉄鋼材料の鉄損計算ルーチンを示すブロック図である。
【図13】マックスウェル方程式により磁気抵抗率νをテンソルとして精度よく求めるルーチンを示す説明図である。
【図14】実験から求められる磁束密度をパラメータとした曲率と磁界の関係図の一例である。
【図15】図14より求めた曲率をパラメータとした磁束密度と磁界増分係数との関係図の一例である。
【図16】実験から求められる磁束密度をパラメータとした曲率と鉄損との関係図の一例である。
【図17】図16より求めた曲率をパラメータとした磁束密度と鉄損増分係数との関係図の一例である。
【図18】図15より求めた曲率ありの場合におけるH−B曲線の一例を示す図である。
【図19】図17より求めた曲率ありの場合におけるW−B曲線の一例を示す図である。
【図20】本発明のハードウエア構成を例示する図である。
【図21】鉄芯における曲率の磁束密度への影響例を示した説明図である。
【図22】変圧器鉄芯におけるコアの積み上げ状態を示した説明図である。
【図23】実施例で使用したステータコア概形図である。
【図24】本発明の実施例で使用した焼嵌め後の曲率半径分布を示す図である。
【図25】本発明による効果を示す鉄損比を示す図である。
【符号の説明】
【0121】
1 領域分割手段、
2 磁束密度決定手段、
3 磁束密度及び鉄損のデータベース、
4 補間内挿計算手段、
5 微小領域内鉄損計算手段、
6 鉄損総和手段
7 曲率分布算出手段
11 凹型鉄芯
12 円筒リング
20 ステータコア
【特許請求の範囲】
【請求項1】
軟磁性材料の電磁場解析対象領域を複数の微小領域に分割する領域分割手段と、該微小領域における磁性体の予め定められた方向と磁束密度Bの方向との間の成す角度θおよび軟磁性材料の弾性曲げ曲率Kを異方性のパラメータとして、磁束密度Bと磁界Hとを関係付けるH−B曲線および磁束密度Bと鉄損Wとを関係付けるW−B曲線を格納するデータベースと、該データベースに格納されているH−B曲線を基にして、前記微小領域におけるマックスウエル方程式および弾性曲げ曲率Kに基づいて、前記角度θおよび前記磁束密度Bの大きさを決定する磁束密度ベクトル決定手段と、前記データベースに格納されているW−B曲線を基にして、前記微小領域の鉄損Wを計算する鉄損計算手段と、前記微小領域の鉄損Wの総和を求める鉄損総和手段とを有することを特徴とする電磁場解析システム。
【請求項2】
前記磁束密度Bおよび磁界Hは、該磁束密度Bと磁界Hとの位相差θBHを考慮して関係付けられることを特徴とする請求項1に記載の電磁場解析システム。
【請求項3】
前記磁束密度Bは、時間高調波成分を含む磁束密度であることを特徴とする、請求項1または請求項2に記載の電磁場解析システム。
【請求項4】
前記磁束密度ベクトル決定手段は、前記弾性曲げ曲率Kを算出する弾性曲げ曲率分布算出手段を備えることを特徴とする請求項1〜3のいずれか1項に記載の電磁場解析システム。
【請求項5】
前記軟磁性材料は、電気機器を構成する軟磁性材料であり、前記弾性曲げ曲率Kは、該電気機器の製作時あるいは自重により発生する弾性曲げ曲率であることを特徴とする請求項1〜4のいずれか1項に記載の電磁場解析システム。
【請求項6】
前記磁束密度ベクトル決定手段は、前記微小領域における磁束密度ベクトルを、その大きさBmaxと、前記の磁性体の予め定められた方向とのなす角度θとに分解し、前記H−B曲線から、前記角度θにおけるH−B曲線を導出し、前記微小領域における磁束密度およびその前記角度θ、磁界とが、その曲線上に存在するように、収束計算を行うことを特徴とする請求項1〜5のいずれか1項に記載の電磁場解析システム。
【請求項7】
前記H−B曲線およびW−B曲線は、前記弾性曲げ曲率Kをパラメータとして、磁束密度と磁界の関係より磁界増分係数α(=HK/HK=0)、磁束密度と鉄損値の関係より鉄損増分係数β(WK/WK=0)を求め、前記弾性曲げ曲率K=0のときのH−B曲線およびW−B曲線に前記磁界増分係数αと鉄損増分係数βとを乗じて得られるH−B曲線およびW−B曲線であることを特徴とする請求項1〜6のいずれか1項に記載の電磁場解析システム。
【請求項1】
軟磁性材料の電磁場解析対象領域を複数の微小領域に分割する領域分割手段と、該微小領域における磁性体の予め定められた方向と磁束密度Bの方向との間の成す角度θおよび軟磁性材料の弾性曲げ曲率Kを異方性のパラメータとして、磁束密度Bと磁界Hとを関係付けるH−B曲線および磁束密度Bと鉄損Wとを関係付けるW−B曲線を格納するデータベースと、該データベースに格納されているH−B曲線を基にして、前記微小領域におけるマックスウエル方程式および弾性曲げ曲率Kに基づいて、前記角度θおよび前記磁束密度Bの大きさを決定する磁束密度ベクトル決定手段と、前記データベースに格納されているW−B曲線を基にして、前記微小領域の鉄損Wを計算する鉄損計算手段と、前記微小領域の鉄損Wの総和を求める鉄損総和手段とを有することを特徴とする電磁場解析システム。
【請求項2】
前記磁束密度Bおよび磁界Hは、該磁束密度Bと磁界Hとの位相差θBHを考慮して関係付けられることを特徴とする請求項1に記載の電磁場解析システム。
【請求項3】
前記磁束密度Bは、時間高調波成分を含む磁束密度であることを特徴とする、請求項1または請求項2に記載の電磁場解析システム。
【請求項4】
前記磁束密度ベクトル決定手段は、前記弾性曲げ曲率Kを算出する弾性曲げ曲率分布算出手段を備えることを特徴とする請求項1〜3のいずれか1項に記載の電磁場解析システム。
【請求項5】
前記軟磁性材料は、電気機器を構成する軟磁性材料であり、前記弾性曲げ曲率Kは、該電気機器の製作時あるいは自重により発生する弾性曲げ曲率であることを特徴とする請求項1〜4のいずれか1項に記載の電磁場解析システム。
【請求項6】
前記磁束密度ベクトル決定手段は、前記微小領域における磁束密度ベクトルを、その大きさBmaxと、前記の磁性体の予め定められた方向とのなす角度θとに分解し、前記H−B曲線から、前記角度θにおけるH−B曲線を導出し、前記微小領域における磁束密度およびその前記角度θ、磁界とが、その曲線上に存在するように、収束計算を行うことを特徴とする請求項1〜5のいずれか1項に記載の電磁場解析システム。
【請求項7】
前記H−B曲線およびW−B曲線は、前記弾性曲げ曲率Kをパラメータとして、磁束密度と磁界の関係より磁界増分係数α(=HK/HK=0)、磁束密度と鉄損値の関係より鉄損増分係数β(WK/WK=0)を求め、前記弾性曲げ曲率K=0のときのH−B曲線およびW−B曲線に前記磁界増分係数αと鉄損増分係数βとを乗じて得られるH−B曲線およびW−B曲線であることを特徴とする請求項1〜6のいずれか1項に記載の電磁場解析システム。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図25】
【図24】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図25】
【図24】
【公開番号】特開2008−170320(P2008−170320A)
【公開日】平成20年7月24日(2008.7.24)
【国際特許分類】
【出願番号】特願2007−4556(P2007−4556)
【出願日】平成19年1月12日(2007.1.12)
【出願人】(000006655)新日本製鐵株式会社 (6,474)
【Fターム(参考)】
【公開日】平成20年7月24日(2008.7.24)
【国際特許分類】
【出願日】平成19年1月12日(2007.1.12)
【出願人】(000006655)新日本製鐵株式会社 (6,474)
【Fターム(参考)】
[ Back to top ]