(57)【要約】【課題】  コイル電流の変化に伴う非定常磁界解析を高速処理する解析法を提供すること。【解決手段】  非定常非線形磁界をNewton-Raphson法で解析する場合、前時刻の解を基点にして、それに対する補正項を求めるというプロセスで解を求める。もしくは、［磁束密度の変動量］／［磁界の変動量］で定義される差分透磁率を用いた変動場に関する方程式を解いて、磁気ベクトルポテンシャルの変動量を求め、求めた磁気ベクトルポテンシャルを基に、磁化曲線を用いて、新たに差分透磁率を求めなおし、求めなおす前後の２個の差分透磁率の線形結合で得られる別の差分透磁率を新たな差分透磁率とする。【選択図】図２
2006209273.tif 000002
【特許請求の範囲】
【請求項１】
非定常非線形磁界をNewton-Raphson法で解析する場合、前時刻の解を基点にして、それに対するソース項の変動による補正項を求めるというプロセスで解を求めることを特徴とする磁界解析法。
【請求項２】
［磁束密度の変動量］／［磁界の変動量］で定義される差分透磁率を用いた変動場に関する方程式を解いて、磁気ベクトルポテンシャルの変動量を求め、求めた磁気ベクトルポテンシャルを基に、磁化曲線を用いて、新たに差分透磁率を求めなおし、求めなおす前後の２個の差分透磁率の線形結合で得られる別の差分透磁率を新たな差分透磁率とすることを特徴とする磁界解析法。
【請求項３】
請求項１または２において、行列方程式を反復解法で解く場合に、変動場を表現する未知数ベクトルの初期値として、前時刻における解を用いることを特徴とする磁界解析法。
【請求項４】
請求項１または２において、行列方程式を反復解法で解く場合に、前時刻でのコイル電流等の磁界発生ソースの変動量に関する前記行列方程式の右辺ベクトルをｂとし、現時刻でのコイル電流等の磁界発生ソースの変動量に関する前記行列方程式の右辺ベクトルを
ｂ′として、変動場を表現する未知数ベクトルの初期値として、前時刻における解に因子ｇ＝ｂ・ｂ′／｜ｂ｜｜ｂ′｜をかけたものを用いることを特徴とする磁界解析法。
【請求項５】
コンピュータに請求項１から請求項４までのいずれかの処理を実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁界解析方法及び磁界解析プログラムに関するものである。
【背景技術】
【０００２】
従来の代表的な非線形磁界解析法として、有限要素法による解析法があり、ＩＣＣＧ法による反復解法ならびに透磁率を逐次修正するNewton-Raphson法を併用している。例えば、文献：中田高義・高橋則雄著：「電気工学の有限要素法」１９８６年（森北出版）等に示されている。
【０００３】
【非特許文献１】「電気工学の有限要素法」中田高義・高橋則雄著、森北出版発行、１９８６年
【発明の開示】
【発明が解決しようとする課題】
【０００４】
有限要素法で磁界解析する場合の一般的な方法として、磁気ベクトルポテンシャルＡを用いたＡ法で記述すると、変位電流の影響が無視できる準静的な場の基礎方程式は、
【０００５】
【数１】
2006209273.tif 000003
となる。ここに、μは透磁率、σは導電率、Ｊ0 はコイルの電流密度ベクトル、Ｍ0 は磁石の残留磁束密度ベクトルである。なお、磁石に関しては、便宜上、Ｂ＝μＨ＋Ｍ0 とした。透磁率μは場に依存するため、数（１）は非線形になり、通常Newton-Raphson法を用いて解を求めている。
【０００６】
数１を離散化すると、
【０００７】
【数２】
2006209273.tif 000004
の形に書ける。ここに、ｘj は磁気ベクトルポテンシャルＡの辺ｊに割り当てた未知数を意味している。
【０００８】
ここで、考え方を容易にするために、数２を次のような１次元問題に置き換える。
【０００９】
【数３】
2006209273.tif 000005
【００１０】
Ａijは、１／μを被積分項に含む積分形式で表現できるため、ｂとｘの関係は磁界Ｈと磁束密度Ｂの関係法に類似しており、ＢＨカーブと類似したｂｘカーブを形成する。以上の準備のもとに、Newton-Raphson法における解の求め方を図１を用いて説明する。
【００１１】
図に示す解曲線ｘ＝ｆ(ｂ)はＢ＝Ｂ(Ｈ)と類似した曲線形状を示し、図に示すようにｂの増加とともにｘの伸びが緩慢になる飽和特性をもつ。前時刻ｔ−Δｔにおけるソース項ｂ0 に対する解ｘ0 に対して、現時刻ｔにおけるソース項ｂ1 に対する解ｘ1 を求めることを考える。
【００１２】
既存のNewton-Raphson法では、前時刻における透磁率を用いてまず線形解を求める。図１の（ｂ，ｘ）空間では、勾配ｘ0 ／ｂ0 を用いて、ソース項ｂ1 に対する線形解を求めることに相当する。その後、Newton-Raphson 法により、非線形性による補正項を順次加算していく。しかし、この方法では、磁気飽和が顕著な場合、始めの線形解におけるｘの増加が顕著になり、収束のプロセスにおいて動作点（ｂ，ｘ）がｂの正の方向に大きく振られるため、解の収束性は大きく悪化することになり、多大な計算時間を要することになる。
【課題を解決するための手段】
【００１３】
上記目的を達成するために本発明に関る磁界解析法では、図２あるいは図３に示すように、前時刻の動作点から現時刻の動作点に直接近づくような解の求め方をすれば良い。その方法として、２つの方法が考えられる。
【００１４】
ひとつは、図２に示すように、前時刻の動作点から接線を引いて、近似解を求め、そこからNewton-Raphson法で解の精度を高めていく方法である。図に示すように、近似解は、現時刻の動作点の近傍にいるため、この場合のNewton-Raphson法による解の収束は極めて速い。
【００１５】
もうひとつの方法は、図３に示すように、現時刻の動作点をあらかじめ予測して、前時刻の動作点から現時刻の動作点あるいはその近傍に向かう線分を引き、これに相当する方程式を解いて、近似解を求め、この近似解により、前述の線分の傾きを修正して、再度方程式を解いて、解の精度を高める。このようなプロセスの反復により、高速に解を求めることができる。
【発明の効果】
【００１６】
本発明に関する磁界解析法によれば、磁性材料を含む非線形磁界を従来のNewton-
Raphson法よりも高速に解析できるという効果がある。
【発明を実施するための最良の形態】
【００１７】
本実施例は、磁性材料の非線形磁界解析の方法に係り、特に非線形磁界を高速に求めるための磁界解析法に関するものである。
【００１８】
用いる磁界解析法として、一般的な有限要素法をとりあげる。本発明の第一実施例を述べる。
【００１９】
第一の実施例として、Newton-Raphson法の高速化に関する方法について述べる。磁気ベクトルポテンシャルを有限要素法で解く場合、
【００２０】
【数４】
2006209273.tif 000006
のようにベクトル型内挿関数Ｎj で展開する。ここに、ａj は有限要素体系を構成する要素群の辺ｊに割り当てた未知変数である。解くべき方程式は、
【００２１】
【数５】
2006209273.tif 000007
【００２２】
【数６】
2006209273.tif 000008
【００２３】
【数７】
2006209273.tif 000009
【００２４】
【数８】
2006209273.tif 000010
【００２５】
【数９】
2006209273.tif 000011
と表現できる。ここに、νは透磁率の逆数である磁気抵抗率、σは導電率、Δｔは時間ステップ幅、ｂi は右辺ソース項である。この方程式を解くことにより、非線形性による補正項δａi を順次求めていく。ここで、補正前の近似解として、前時刻の解をそのまま使う。この操作は、図１に示したような接線を引いて近似解を求めるような操作に相当する。新しい近似解を求めたら、νとそのＢ2 に関する微分を求め直して、残差Ｈi の絶対値がある収束判定値よりも小さくなれば、収束完了ということになる。
【００２６】
本実施例によれば、従来のNewton-Raphson法による解析手順、すなわち、前時刻における磁気抵抗率νを用いて、Ｈi ＝０を解いてまず線形解を求め、その後、非線形性による補正項を求めて加算していくというプロセスに比べて、高速に収束解を得られるという効果がある。
【００２７】
つぎに第二の実施例について述べる。まず、差分透磁率μ′を次のように定義する。
【００２８】
【数１０】
2006209273.tif 000012
【００２９】
ここに、ΔＢおよびΔＨは、それぞれ磁束密度の変動量と磁界の変動量である。これを用いて、まず、磁気ベクトルポテンシャルの変動量に関する基礎方程式を導く。Maxwell方程式の一部をなすFaradayの式である
【００３０】
【数１１】
2006209273.tif 000013
【００３１】
【数１２】
2006209273.tif 000014
を出発点とする。ここに、添え字１，０は、それぞれ現時刻ｔ、並びに前時刻ｔ−Δｔにおける物理量であることを意味する。ここで、ΔＡ＝Ａ1−Ａ0として、Ａの時間微分項に関するNewmarkのβ法に基づく次式を用いる。
【００３２】
【数１３】
2006209273.tif 000015
【００３３】
ここに、βは０≦β≦１なる固定定数であり、安定した時間応答解析のために１／２≦β≦１なるβを用いる必要がある。ここでβ×数１１＋(１−β)×数１２に数１３を用いると、
【００３４】
【数１４】
2006209273.tif 000016
ここで、ΔＨ＝Ｈ1−Ｈ0，ΔＪ＝Ｊ1−Ｊ0を用いて変形すると、
【００３５】
【数１５】
2006209273.tif 000017
全体をβで割って適宜移項すると、
【００３６】
【数１６】
2006209273.tif 000018
を得る。ここで差分透磁率μ′を用いて、
【００３７】
【数１７】
2006209273.tif 000019
と書けるため、数１６は、
【００３８】
【数１８】
2006209273.tif 000020
となる。特にβ＝１では、
【００３９】
【数１９】
2006209273.tif 000021
数１８のrotＨ0の取り扱い方によって、次に示す３種類の解法が考えられる。
【００４０】
第一の方法は、rotＨ0にAmpereの式を用いる方法である。rotＨ0は数１２より
【００４１】
【数２０】
2006209273.tif 000022
と書けるので、これを数１８に代入すると、
【００４２】
【数２１】
2006209273.tif 000023
となる。右辺の∂Ａ0／∂ｔは数１３より各時間ステップで∂Ａ／∂ｔを更新することにより、求めることができる。特にβ＝１では、
【００４３】
【数２２】
2006209273.tif 000024
第二の方法は、透磁率μを用いる方法である。
【００４４】
【数２３】
2006209273.tif 000025
を数１１の右辺に代入すると、
【００４５】
【数２４】
2006209273.tif 000026
を得る。特にβ＝１では、
【００４６】
【数２５】
2006209273.tif 000027
第三の方法は、ΔＨの時間積算を用いる方法である。時刻ｔ＝０におけるＨ，Ｊ，Ａ、をそれぞれＨ0，Ｊ0，Ａ0とおくと、
【００４７】
【数２６】
2006209273.tif 000028
【００４８】
【数２７】
2006209273.tif 000029
の関係式を用いると、数１７から
【００４９】
【数２８】
2006209273.tif 000030
を得る。これを数１１に代入すると
【００５０】
【数２９】
2006209273.tif 000031
となる。特にβ＝１では、
【００５１】
【数３０】
2006209273.tif 000032
となる。
【００５２】
本発明における非線形動磁場解析のフローチャートを図４に示す。この図は第３時間ステップ目以降の解析の流れを示している。まず、前時刻におけるソース項の変動量Δｂ0 と求めた磁界の変動量ΔＨ0 から、現時刻におけるソース項の変動量Δｂ1 に対する磁界の変動量ΔＨ1 は、
【００５３】
【数３１】
2006209273.tif 000033
で予測できる。これより、前時刻における磁界の変動量ΔＨ0 から現時刻における磁界の変動量の第ゼロ近似値ΔＨ1(0)を推定する。次に初期磁化曲線により磁束密度の変動量
ΔＢ1(0)を求め、各要素における差分透磁率μ1′(0)を求める。近似的に求めた差分透磁率μ1(0)を用いて、次章に示す場の方程式を解くと、より正確な解（Ｈ1(1)，Ｂ1(1)）が求まる。この点と前時刻に得られた動作点（Ｈ0，Ｂ0）とを結ぶ線分の勾配として、より正確な差分透磁率μ1′(0)が求まる。そこで、２つの差分透磁率μ1′(0)とμ1′(1)を用いて、さらに正確な差分透磁率μ1′(1)を求める。ｎ回目の反復における差分透磁率
μ1′nは、次式のようになる。
【００５４】
【数３２】
2006209273.tif 000034
【００５５】
ここに、μdは前時刻に得られた動作点（Ｈ0，Ｂ0）における微分透磁率である。数
（８）におけるμ1′(n)とμ1′n-1の線形結合係数は、経験的に求めた値であり、
μ1′(n)よりも正確な差分透磁率を求める式になっている。
【００５６】
この係数は、数３２のものと多少ずれたものでも良い。この差分透磁率μ1′nならびに数２１，数２４，数２９のいずれかの式を用いて、有限要素法により、ΔＡ1(1)，
ΔＢ1(1)を求め、Ｂ1(1)＝Ｂ0＋ΔＢ1(1)より、磁束密度を求める。以上の操作を
｜ΔＢ1(n)｜＜ε（ε：収束判定値）を満足するまで反復する。
【００５７】
本実施例によれば、反復計算において、既に求めてある変動場に関する解を初期値にして、行列方程式を反復解法で計算するために、Newton-Raphson法よりも高速に解を求められるという効果がある。
【００５８】
図５に解析システムの一例の図を示す。本解析システムは、計算機１，表示装置２，記憶媒体３から構成される。ここでは、記憶媒体３を明示するために、計算機１の外に出しているが、計算機１内部に記憶媒体３を設置しても良い。計算機１には、上記に示した２つの実施例のいずれかのアルゴリズムに基づくプログラムが格納されており、計算機１で計算を実行し、計算結果を表示装置２に表示する。計算結果の一部は記憶媒体３に記憶される。時間ステップを刻みながらの解析になるため、記憶媒体３に記憶したデータを再利用しながら、計算を進めていくことになる。
【図面の簡単な説明】
【００５９】
【図１】従来のNewton-Raphson法による解の収束過程を示す図。
【図２】本発明の第一の実施例である新しいNewton-Raphson法による解の収束過程を示す図。
【図３】本発明の第二の実施例である差分透磁率を用いたときの解の収束過程を示す図。
【図４】本発明の第二の実施例である差分透磁率を用いたときの非線形動磁場解析のフローチャート計算を示す図。
【図５】本発明による解析システムの一例を示す図。
【符号の説明】
【００６０】
１…第１および第２の実施形態における計算機、２…第１および第２の実施形態における表示装置、３…第１および第２の実施形態における記憶媒体、１１…本発明の第二の実施例における計算ステップ１、１２…本発明の第二の実施例における計算ステップ２、
１３…本発明の第二の実施例における計算ステップ３、１４…本発明の第二の実施例における計算ステップ４、１５…本発明の第二の実施例における計算ステップ５、１６…本発明の第二の実施例における計算ステップ６、１７…本発明の第二の実施例における計算ステップ７。
【図１】
2006209273.tif 000035
【図２】
2006209273.tif 000036
【図３】
2006209273.tif 000037
【図４】
2006209273.tif 000038
【図５】
2006209273.tif 000039