(57)【要約】【課題】回転機の磁界分布を、特定部位における特定高調波でフィルタリングする回転機の磁界解析法を提供する。【構成】解析空間内の第１の特定部位に与えられた電磁気的な変数に関する第１の境界条件に基づき、解析空間の磁界解析を行う第１の磁界解析手順１と、第１の磁界解析手順１により得られた解析空間内の第２の特定部位における電磁気的な変数を高調波成分毎に展開する展開手順２と、展開した電磁気的な変数の高調波成分のうちの少なくとも一つの高調波成分に関する電磁気的な変数を第２の特定部位に第２の境界条件として与える第２の境界条件設定手順と、第１及び第２の境界条件に基づいて解析空間の磁界解析を行う手順３とをコンピュータに実行させる。【選択図】  図１
2004245814.tif 000002
【特許請求の範囲】
【請求項１】
解析空間内の第１の特定部位に与えられた電磁気的な変数に関する第１の境界条件に基づき、解析空間の磁界解析を行う第１の磁界解析手順と、
前記第１の磁界解析手順により得られた前記解析空間内の第２の特定部位における電磁気的な変数を高調波成分毎に展開する展開手順と、
展開した前記電磁気的な変数の高調波成分のうちの少なくとも一つの高調波成分に関する電磁気的な変数を前記第２の特定部位に第２の境界条件として与える第２の境界条件設定手順と、
前記第１及び第２の境界条件に基づいて解析空間の磁界解析を行う手順と
をコンピュータに実行させるための回転電機の磁界解析プログラム。
【請求項２】
請求項１において、前記展開手順では前記電磁気的な変数を空間的な分布に関する高調波成分毎に展開することを特徴とする回転電機の磁界解析プログラム。
【請求項３】
請求項１において、前記展開手順では前記電磁気的な変数を時間的な分布に関する高調波成分毎に展開することを特徴とする回転電機の磁界解析プログラム。
【請求項４】
請求項１において、前記電磁気的な変数はポテンシャルまたは起磁力であることを特徴とする回転電機の磁界解析プログラム。
【請求項５】
請求項１において、前記第１の磁界解析手順の前に前記解析空間の透磁率分布を記憶する記憶手順を含み、
前記第１の磁界解析手順及び前記第２の磁界解析手順においては、前記記憶手順により記憶された透磁率分布に基づき解析空間の磁界解析を行うことを特徴とする回転電機の磁界解析プログラム。
【請求項６】
請求項１において、前記解析空間は固定子に対応する解析空間と回転子に対応する解析空間を有し、
前記第２の特定部位は前記固定子に対応する解析空間と前記回転子に対応する解析空間の境界であることを特徴とする回転電機の磁界解析プログラム。
【請求項７】
解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えて、回転機全体の磁界解析を行う全体磁界解析手順と、
前記全体磁界解析手順により得られた磁性体内部の透磁率分布を記憶する手順と、
前記全体磁界解析により得られた回転子と固定子との間のスライド面上における磁界分布を記述するポテンシャルを、回転軸方向にはある分布を有し回転方向に関して定数あるいは周期的な分布をもつ複数の高調波成分に展開する手順と、
前記全体磁界解析手順における境界条件と同じ境界条件を与え、さらに、前記スライド面上にはディリクレ条件を与え、記憶した前記透磁率分布を解析空間全体に与えて固定し、前記全体磁界解析手順で与えた電流或いは磁石磁力を与えて、回転子を含む回転子空間と、固定子を含む固定子空間と、のそれぞれに線形の磁界解析を実施する手順と、
前記全体磁界解析手順における境界条件と同じ境界条件を与え、さらに、前記高調波成分を順次前記スライド面上に境界条件として与え、さらに、前記透磁率分布を磁性体全体に与え、電流或いは磁石磁力を零にして、回転子を含む回転子空間と固定子を含む固定子空間とのそれぞれに線形の磁界解析を実施する手順とをコンピュータに実行させるための回転電機の磁界解析プログラム。
【請求項８】
三次元の回転機構造体を分析する場合に、
解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えて回転機の磁界解析を実施する手順において、ゲージフリー又は部分的にゲージフリーの条件のもとで得られたポテンシャルを、ゲージ条件を満足するポテンシャルに変換する手順を、前記時間高調波に展開する手順の前に行うことを特徴とする請求項７に記載の回転電機の磁界解析プログラム。
【請求項９】
さらに、上記一連の回転子空間と固定子空間と、のそれぞれ個別の解析により得られた磁界分布を、与えた前記境界条件別に個別に表示させる手順を含むことを特徴とする請求項７又は８に記載の回転電機の磁界解析プログラム。
【請求項１０】
解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えて、一連の時間ステップにおける回転機全体の磁界解析を行う全体磁界解析手順と、
前記全体磁界解析手順により得られた前記一連の時間ステップにおける磁性体内部の透磁率分布を記憶する手順と、
前記全体磁界解析手順により得られた回転子と固定子との間のスライド面上における磁界分布を記述するポテンシャルを、回転軸方向にはある分布を有し回転方向に関して定数あるいは周期的な分布をもつ複数の高調波成分に展開し記憶する手順と、
一連の時間ステップにおける上記解析から得られた、回転方向に関して定数あるいは周期的な分布をもつ複数の高調波成分を時間軸方向に高調波展開して、時間・空間高調波成分を抽出する手順と、
前記全体磁界解析における境界条件と同じ境界条件を与え、さらに、前記スライド面上にはディリクレ条件を与え、前記全体磁界解析で与えた電流又は磁石磁力を与え、前記一連の時間ステップにおける前記磁性体内部の前記透磁率分布を時間ステップ毎に与え、回転子を含む回転子空間と、固定子を含む固定子空間と、のそれぞれにおいて線形の磁界解析を実施する手順と、
前記全体磁界解析手順における境界条件と同じ境界条件を与え、前記時間・空間高調波成分を順次前記スライド面上に境界条件として与え、さらに、一連の時間ステップにおける磁性体内の前記透磁率分布を時間ステップ毎に固定して与え、電流或いは磁石磁力を零にして、回転子を含む回転子空間と固定子を含む固定子空間とのそれぞれに線形の磁界解析を実施する手順と
をコンピュータに実行させるための回転電機の磁界解析プログラム。
【請求項１１】
三次元の回転機構造体を分析する場合に、
解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えて回転機の磁界解析を実施する手順において、ゲージフリー又は部分的にゲージフリーの条件のもとで得られたポテンシャルを、ゲージ条件を満足するポテンシャルに変換する手順を、前記時間高調波に展開する手順の前に行うことを特徴とする請求項１０に記載の回転電機の磁界解析プログラム。
【請求項１２】
さらに、上記一連の回転子空間と固定子空間と、のそれぞれ個別の解析により得られた磁界分布を、与えた前記境界条件別に個別に表示させる手順を含むことを特徴とする請求項１０又は１１に記載の回転電機の磁界解析プログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、回転機の磁界解析技術に関し、特に回転機の磁界分布を構成する複数の高調波成分と漏れ磁界成分に分解することができる回転機の磁界解析技術に関するものである。
【０００２】
【従来の技術】
これまでの回転機の磁界解析法において、回転子を小刻みに回転させながら逐次解析する時間ステップ法が採用されている（例えば、非特許文献１参照）。
【０００３】
【非特許文献１】
Ｔ．Ｗ．Ｐｒｅｓｔｏｎ，  Ａ．Ｂ．Ｊ．Ｒｅｅｃｅ  ａｎｄ  Ｐ．Ｓ．Ｓａｎｇｈａ：  “Ｉｎｄｕｃｔｉｏｎ  ｍｏｔｏｒ  ａｎａｌｙｓｉｓ  ｂｙ  ｔｉｍｅ−ｓｔｅｐｐｉｎｇ  ｔｅｃｈｎｉｑｕｅｓ”，  ＩＥＥＥ  Ｔｒａｎｓ．  ｏｎ  Ｍａｇｎｅｔｉｃｓ，  ｖｏｌ．２６，  Ｎｏ．２，  ｐｐ．５５１−５５４，  １９９０
【０００４】
【発明が解決しようとする課題】
振動、騒音、鉄損、渦電流損などは、特定の高調波成分によって引き起こされる。ところが、上記従来技術による回転機の磁界解析方法では、解析結果として得られる磁界分布に複数の高調波成分が含まれ、高調波成分単独での磁界解析ができなかった。すなわち、従来の磁界解析結果として得られる磁界分布は複数の高調波成分全てが重畳された状態で得られるため、問題となる高調波を特定したり、振動、騒音、鉄損、渦電流損などの低減方法を策定したりする場合に有用な特定高調波毎の磁界分布が得られないという問題があった。さらに、解析対象を２次元解析のみではなく３次元解析にまで拡張できれば便利である。
【０００５】
本発明の目的は、特定の空間高調波又は時間高調波の磁界分布を表示することにより、回転機の設計者が振動・騒音等を分析するのに有益な情報を得ることできる回転機の磁界解析技術を提供することにある。
【０００６】
【課題を解決するための手段】
本発明の一観点によれば、解析空間内の第１の特定部位に与えられた電磁気的な変数に関する第１の境界条件に基づき、解析空間の磁界解析を行う第１の磁界解析手順と、前記第１の磁界解析手順により得られた前記解析空間内の第２の特定部位における電磁気的な変数を高調波成分毎に展開する展開手順と、前記分析した電磁気的な変数の高調波成分のうちの少なくとも一つの高調波成分に関する電磁気的な変数を前記第２の特定部位に第２の境界条件として与える第２の境界条件設定手順と、前記第１及び第２の境界条件に基づいて解析空間の磁界解析を行う手順とをコンピュータに実行させるためのプログラムが提供される。
【０００７】
また、解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えて、回転機全体の磁界解析を行う全体磁界解析手順と、前記全体磁界解析手順により得られた磁性体内部の透磁率分布を記憶する手順と、前記全体磁界解析により得られた回転子と固定子との間のスライド面上における磁界分布を記述するポテンシャルを、回転軸方向にはある分布を有し回転方向に関して定数あるいは周期的な分布をもつ複数の高調波成分に展開する手順と、前記全体磁界解析手順における境界条件と同じ境界条件を与え、さらに、前記スライド面上にはディリクレ条件を与え、記憶した前記透磁率分布を解析空間全体に与えて固定し、前記全体磁界解析手順で与えた電流或いは磁石磁力を与えて、回転子を含む回転子空間と、固定子を含む固定子空間のそれぞれにおいて線形の磁界解析を実施する手順と、前記全体磁界解析手順における境界条件と同じ境界条件を与え、さらに、前記高調波成分を順次前記スライド面上に境界条件として与え、さらに、前記透磁率分布を磁性体全体に与え、電流或いは磁石磁力を零にして、回転子を含む回転子空間と固定子を含む固定子空間とのそれぞれに線形の磁界解析を実施する手順とをコンピュータに実行させるための回転電機の磁界解析プログラムが提供される。
【０００８】
全体磁界解析により得られた回転子と固定子との間のスライド面上における磁界分布を記述するポテンシャルを、回転軸方向にはある分布を有し回転方向に関して定数或いは周期的な分布をもつ高調波成分に展開する手順を含めることにより、解析対象を３次元に拡張することができる。
【０００９】
三次元の回転機構造体を分析する場合に、解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えて回転機の磁界解析を実施する手順において、ゲージフリー又は部分的にゲージフリーの条件のもとで得られたポテンシャルをゲージ条件を満足するポテンシャルに変換する手順を、ポテンシャルを、回転軸方向にはある分布を有し回転方向に関して定数あるいは周期的な分布をもつ複数の高調波成分に展開する手順の前に行うのが好ましい。
【００１０】
【発明の実施の形態】
本明細書においては、定数成分や基本波成分も含めて高調波と称する。
以下、本発明に係る回転機の磁界解析技術について図面に参照しつつ説明を行う。まず、本発明の第１の実施の形態による磁界解析方法について図１を参照して説明する。図１に示すように、本実施の形態による磁界解析方法では、第１ステップとして、従来技術に従った全体磁界解析１を実施する。このとき、回転子は計算条件として設定される回転数で回転させながら、逐次、時間ステップ法で解析を進める。磁界のソースとしては、電流、磁石などの起磁力を与える。また、非線形磁気特性を考慮した非線形解析とする場合は、時刻ステップ毎に、各要素の透磁率データ５を格納するか又は記憶していく。第２ステップとして、第１ステップの解析結果として得られたポテンシャルのうち、特定部位におけるポテンシャルの高調波展開２を行う。高調波としては、空間的な分布に含まれる空間高調波と、時間変化に含まれる時間高調波があるが、どちらか一方を展開（分離）する場合、或いは双方とも展開する場合が含まれる。高調波成分への展開の具体的な手法としては、公知のフーリエ変換（フーリエ展開）あるいはウェーブレット変換等の手法を用いることができる。
【００１１】
第３ステップとして、成分磁界解析３を実施する。成分磁界解析３は、ソース成分磁界解析３ａと、ポテンシャル高調波を境界条件とする成分磁界解析３ｂとを有する。成分磁界解析３では、全体解析１で得られた各要素の透磁率データ５を、時刻ステップ毎に各要素の透磁率として読み込んで与える。ソース成分磁界解析３ａでは、磁界のソースとして、電流、磁石磁力を与えた線形磁界解析を実施する。この際、電流、磁石磁力としては、全体磁界解析１で与えた電流、磁石磁力をそのものを与えて一回の磁界解析を行っても良いし、総和が全体磁界解析１で与えた電流、磁石磁力となるように成分分離させた各電流、磁石磁力成分を与えて複数回の磁界解析を行っても良い。
【００１２】
ポテンシャル高調波を境界条件とする成分磁界解析３ｂでは、特定部位におけるポテンシャルの高調波展開２で得られたポテンシャルの高調波成分を特定部位上に境界条件として与え、磁界を解析する。高調波成分は複数あるため、成分磁界解析３ｂも複数の高調波成分に対して複数回実施する。第４ステップは、任意成分磁界解析結果の重畳４を実施するステップであり、複数の特定高調波成分からなる磁界分布を知りたい場合に頭に、任意に実行することができる。
【００１３】
図２は本発明の第２の実施の形態による磁界解析法の流れを示し、図３は、本実施の形態による磁界解析法の補足説明図である。図２に示すように、本実施の形態による磁界解析法では、外周側解析空間１２と内周側解析空間１３とからなる全体解析空間１１があり、高調波としては、空間高調波と時間高調波の双方を展開（分離）する。全体解析空間１１は、軸方向断面の二次元空間を想定し、磁界解析の変数は、磁気ベクトルポテンシャルの軸方向成分を想定して説明する。外周側解析空間１２、内周側解析空間１３は、それぞれ固定子側、回転子側、或いは、それぞれ回転子側、固定子側のいずれかとなる。外周側解析空間１２、内周側解析空間１３の境界となるｂｂ’間の円弧は、固定子と回転子の間のエアギャップ中に位置させる。
【００１４】
全体磁界解析２１の方法を説明する。全体解析において、全体解析空間１１の境界条件は、外周側円弧上のｓｓ’間と内周円弧上のｒｒ’間ポテンシャルを、
【００１５】
【数１】
2004245814.tif 000003
とおき、周方向端部の境界ｓｒ間とｓ’ｒ’間は互いに周期境界条件となるようにする。全体解析空間１１が一周にわたる解析領域として設定されている場合、周期境界条件は用いない。磁界解析では、磁気特性データに基づいた非線形解析とし、磁界のソースとなる電流、磁石磁力等を入力し、円弧ｂｂ’上において回転子を移動させながら、逐次解析する時間ステップ法を実施する。全体解析空間１１内に導体がある場合は、導電率を入力して渦電流を解析すればよい。尚、透磁率を一定とした線形解析、渦電流を考慮しない解析であってもよい。ここで、各時刻での各要素の透磁率は、各要素の透磁率データ５に格納しておく。
【００１６】
ポテンシャルａの空間高調波、時間高調波分析２２では、全体磁界解析２１で得られた円弧ｂｂ’上でのポテンシャルを、
【００１７】
【数２】
2004245814.tif 000004
のように空間高調波と時間高調波の成分に展開する。尚、ポテンシャルａの高調波成分を展開する円弧ｂｂ’は、内周側解析空間１３を回転子と見立てて移動させる場合、図３に示すように、静止系の座標に円弧ｂｂ’として固定する場合と、回転系の座標に円弧ｃｃ’として固定する場合とが考えられるが、静止系の座標から回転系の座標へは座標変換で変換できるため何れの方法を用いても良い。言い換えれば、静止系の座標に円弧ｂｂ’上で展開して、回転系の座標に円弧ｃｃ’上に座標変換することもできるし、回転系の座標に円弧ｃｃ’上で分析して、静止系の座標に円弧ｂｂ’上に座標変換することもできる。
【００１８】
次いで、固定子、回転子漏れ磁束の磁界解析２３を実施する。固定子、回転子漏れ磁束の磁界解析２３では、全体磁界解析２１で得られた各要素の透磁率データベース５を、各時刻ステップ毎に各要素の透磁率として読み込んで与える。境界条件としては、全体磁界解析２１と同じく、全体解析空間１１の外周側円弧上のｓｓ’間と内周円弧上のｒｒ’間ポテンシャルは、上記（１）式とおき、周方向端部の境界ｓｒ間とｓ’ｒ’間は互いに周期境界条件となるようにする。また、外周側解析空間１２側の円弧ｂｂ’上、内周側解析空間１３側の円弧ｃｃ’上には、上記（１）式と同じ、
【００１９】
【数３】
2004245814.tif 000005
の境界条件を設定することで、外周側解析空間１２内および内周側解析空間１３内で漏れ磁束のみを発生させるようにする。磁界のソースとしては、全体磁界解析２１で与えた電流、磁石磁力を与え、線形磁界解析とし、全体磁界解析２１導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を実施すればよい。このとき電流、磁石磁力としては、全体解析１で与えた電流、磁石磁力をそのものを与えて１回の磁界解析を行っても良いし、総和が全体解析１で与えた電流、磁石磁力となるように成分分離させた各電流、磁石磁力成分を与えて複数回の磁界解析を行っても良い。また、この固定子、回転子漏れ磁束の磁界解析２３では、外周側解析空間１２と内周側解析空間１３の境界となる、円弧ｂｂ’上、円弧ｃｃ’上にそれぞれ上記（３）式の境界条件を与えるため、外周側解析空間１２と内周側解析空間１３とは、個別に解析することが可能になる。
【００２０】
固定子、回転子漏れ磁束の磁界解析２３と平行して、特定空間高調波、時間高調波磁界解析２４を実施する。特定空間高調波、時間高調波磁界解析２４では、固定子、回転子漏れ磁束の磁界解析２３と同様に、全体磁界解析２１で得られた各要素の透磁率データ５を、各時刻ステップ毎に各要素の透磁率として読み込んで与える。境界条件としては、固定子、回転子漏れ磁束の磁界解析２３と同様に、全体磁界解析２１と同じく、全体解析空間１１の外周側円弧上のｓｓ’間と内周円弧上のｒｒ’間ポテンシャルは（１）式とおき、周方向端部の境界ｓｒ間とｓ’ｒ’間は互いに周期境界条件となるようにする。外周側解析空間１２側の円弧ｂｂ’上、内周側解析空間１３側の円弧ｃｃ’上には、固定子、回転子漏れ磁束の磁界解析２３と異なり、（２）式で分析したポテンシャルのうち、特定のｍ次空間高調波、ｌ次時間高調波のみを抽出した、
【００２１】
【数４】
2004245814.tif 000006
を境界条件として設定する。磁界のソースは、全体磁界解析２１で与えた電流、磁石磁力を与えずに、線形磁界解析とし、全体磁界解析２１において導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を行えばよい。なお、円弧ｂｂ’上と円弧ｃｃ’上では、空間高調波次数が異なることになるが、これは一方の円弧上の展開結果に対し、静止系と回転系の間での座標変換を実施して他方の円弧上に与えることができる。
【００２２】
固定子、回転子漏れ磁束の磁界解析２３と同様に、特定空間高調波、時間高調波磁界解析２４でも、外周側解析空間１２と内周側解析空間１３の境界となる、円弧ｂｂ’上、円弧ｃｃ’上にそれぞれ（３）式の境界条件を与えるため、外周側解析空間１２と内周側解析空間１３とは、個別に解析することができる。
【００２３】
また、固定子、回転子漏れ磁束の磁界解析２３と、特定空間高調波、時間高調波磁界解析２４は、順不同で、任意に実施可能であり、特定空間高調波、時間高調波磁界解析２４における特定のｍ次空間高調波、ｌ次空間高調波は、解析者の欲する次数のみ着目して解析すればよい。
【００２４】
図４は本発明の第３の実施の形態による磁界解析法の流れを示す図である。図４に示す解析例は、図２に示す解析例に対し、取り扱う高調波を空間高調波のみとしたものである。全体磁界解析２１は図２と同じであり、各要素の透磁率データ５も同じように格納する。ポテンシャルａの空間高調波分析２５では、全体磁界解析２１で得られた円弧ｂｂ’上でのポテンシャルを、
【００２５】
【数５】
2004245814.tif 000007
のように空間高調波の成分に展開する。固定子、回転子漏れ磁束の磁界解析２３は、図２と同じである。
【００２６】
特定空間高調波磁界解析２６では、全体磁界解析２１で得られた各要素の透磁率データ５を、各時刻ステップ毎に各要素の透磁率として読み込んで与える。境界条件としては、全体解析空間１１の外周側円弧上のｓｓ’間と内周円弧上のｒｒ’間ポテンシャルを図２と同じく（１）式とおき、周方向端部の境界ｓｒ間とｓ’ｒ’間は互いに周期境界条件とする。外周側解析空間１２側の円弧ｂｂ’上、内周側解析空間１３側の円弧ｃｃ’上には、（５）式により分析したポテンシャルのうち、特定のｍ次空間高調波のみを抽出した、
【００２７】
【数６】
2004245814.tif 000008
を境界条件として設定する。磁界のソースは、全体磁界解析２１で与えた電流、磁石磁力を与えずに、線形磁界解析とし、全体磁界解析２１導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を実施すればよい。その後、解析した特定空間高調波の磁界分布のみを表示手段に表示する。本実施例によれば、特定空間高調波の磁界分布を表示して回転機の設計者に示すことにより、設計者は振動・騒音発生等の要因を分析するために有益な情報を得ることができる。
【００２８】
図５は本発明の第４の実施の形態による磁界解析法の流れを示す図である。図５に示す解析例は、図２に示す解析例に対し、取り扱う高調波を時間高調波のみとしたものである。全体磁界解析２１は図２と同じであり、各要素の透磁率データ５も同じように格納する。ポテンシャルａの時間高調波分析２５では、全体磁界解析２１で得られた円弧ｂｂ’上でのポテンシャルを、
【００２９】
【数７】
2004245814.tif 000009
のように時間高調波の成分に分析する。固定子、回転子漏れ磁束の磁界解析２３は、図２と同じである。
【００３０】
特定時間高調波磁界解析２６では、全体磁界解析２１で得られた各要素の透磁率データ５を、各時刻ステップ毎に各要素の透磁率として読み込んで与える。境界条件としては、全体解析空間１１の外周側円弧上のｓｓ’間と内周円弧上のｒｒ’間ポテンシャルを図２と同じく（１）式とおき、周方向端部の境界ｓｒ間とｓ’ｒ’間は互いに周期境界条件とする。外周側解析空間１２側の円弧ｂｂ’上、内周側解析空間１３側の円弧ｃｃ’上には、（７）式により分析したポテンシャルのうち、特定の１次時間高調波のみを抽出した、
【００３１】
【数８】
2004245814.tif 000010
を境界条件として設定する。磁界のソースは、全体磁界解析２１で与えた電流、磁石磁力を与えずに、線形磁界解析とし、全体磁界解析２１導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を実施すればよい。その後、解析した特定時間高調波の磁界分布のみを表示手段に表示する。
【００３２】
本実施の形態によれば、特定時間高調波の磁界分布を表示して回転機の設計者に示すことにより、設計者は振動・騒音発生等の要因を分析するために有益な情報を得ることができる。
【００３３】
次に、磁界解析を３次元に拡張する３次元解析法について説明する。図６は、本発明の第５の実施の形態による磁界解析法の流れを示す図である。第１ステップとして、従来技術に沿った全体磁界解析１０１を実施する。この際、電流、磁石などの起磁力を与え、回転子は計算条件として設定される回転数で回転させながら、逐次、時間ステップ法で解析を進める。また、非線形磁気特性を考慮した非線形解析とする場合は、時間ステップ毎に磁性体の各要素の透磁率データ１０５を記憶し格納する。第２ステップとして、第１ステップの解析結果として得られたポテンシャルのうち、特定部位におけるポテンシャルの高調波展開１０２を行う。この場合、空間高調波と時間高調波のどちらか一方を展開することもできるし、双方とも展開することもできる。
【００３４】
第３ステップとして、成分磁界解析１０３を実施する。成分磁界解析１０３は、ソース成分磁界解析１０３ａとポテンシャル高調波を境界条件とする成分磁界解析１０３ｂとを含む。成分磁界解析１０３では、全体解析１０１で得られた各時間ステップにおける磁性体内の各要素の透磁率データ１０５を、磁性体内の各要素の透磁率として固定して割り当て、線形磁界解析を実行する。
【００３５】
ソース成分磁界解析１０３ａでは、磁界のソースとして電流又は磁石磁力の少なくともいずれか一方を与えた線形磁界解析を実施する。この際、電流、磁石磁力としては、全体磁界解析１０１で与えたコイル電流、磁石磁力そのものを与えて一回の磁界解析を行っても良いし、総和が全体磁界解析１０１で与えた電流、磁石磁力となるように成分分離させた各電流、磁石磁力成分を与えて個別に磁界解析を行っても良い。３次元解析に拡張するために、ポテンシャル高調波を境界条件とする成分磁界解析１０３ｂでは、特定部位におけるポテンシャルの高調波分析１０２で得られた回転方向に関する定数成分（回転軸方向にはある分布をもち、回転方向には定数である成分）あるいは高調波成分（回転軸方向にはある分布をもち、回転方向に関しては周期的な分布をもつ複数の空間高調波モード成分）を同じ特定部位上に境界条件として与え、コイル電流、磁石磁力を全て零にして磁界解析する。高調波成分は複数あるため、成分磁界解析１０３ｂも複数の高調波成分に対応させて複数回行う。第４ステップは、任意成分磁界解析結果の重畳１０４を実施する処理を行うものであり、複数の特定高調波成分からなる磁界分布を知りたい場合に、任意に実行すればよい。
【００３６】
図７は、本発明の第６の実施の形態による回転機の磁界解析法の流れを示す図である。図７に示す解析例は、外周側解析空間１１２と内周側解析空間１１３とを有する全体解析空間１１１があり、高調波としては、空間高調波と時間高調波の双方を展開する例である。全体解析空間１１１は回転軸方向に垂直な断面における二次元空間を想定し、磁界解析の変数は、例として辺要素有限要素法において各要素の辺上に配置された磁気ベクトルポテンシャルの辺上への射影線積分値ａを想定して説明する。ポテンシャルを対象とするため、以下の実施例で説明するポテンシャルに、ある定数ベクトルを一様に加えて、辺上の線積分で得られる変数ａを用いても良いことは言うまでもない。
【００３７】
ここでは、簡単のため全体の解析空間の外周面における周期境界以外の境界面での境界条件として、磁界が境界面に平行に沿うディリクレ条件としてａ＝０を設定する場合について説明する。磁界強度が弱い外周面では、このように通常ディリクレ条件としてａ＝０を設定する場合が多い。外周側解析空間１１２、内周側解析空間１１３は、それぞれ固定子側、回転子側、或いは、それぞれ回転子側、固定子側のいずれかとなる。外周側解析空間１１２、内周側解析空間１１３の境界となるスライド面ｂｂ’は、固定子と回転子の間のエアギャップ中に位置させる。全体磁界解析１２１の方法を説明する。全体解析において、全体解析空間１１１の境界条件は、外周面ｓｓ’と内周面ｒｒ’のポテンシャルを０とおき、周方向端部の境界面ｓｒと面ｓ’ｒ’は互いに周期境界条件で接続される。全体解析空間１１１が一周にわたる解析領域として設定されている場合、周期境界条件は用いない。
【００３８】
磁界解析では、磁気特性データに基づいた非線形解析とし、磁界のソースとなる電流、磁石磁力等を入力し、スライド面ｂｂ’上に沿って回転子を回転しながら、時間ステップ法を用いて逐次解析する。全体解析空間１１１内に導体がある場合は、導電率を入力して渦電流を解析すればよい。尚、透磁率を時間的に一定とした線形解析や、渦電流を考慮しない解析であってもよい。ここで、各時刻での磁性体内各要素の透磁率は、透磁率データ１０５に格納しておく。ポテンシャルａの空間高調波、時間高調波成分への展開１２２では、全体磁界解析１２１で得られたスライド面ｂｂ’上でのポテンシャルを、
【００３９】
【数９】
2004245814.tif 000011
のように空間高調波並びに時間高調波に成分展開する。この場合、静止系で成分展開しても良いし、回転子にのった回転座標系で成分展開しても良い。
【００４０】
次いで、固定子、回転子漏れ磁束の磁界解析１２３を実施する。固定子、回転子漏れ磁束の磁界解析１２３では、全体磁界解析１２１で得られた磁性体内各要素の透磁率データ１０５を、時間ステップ毎に磁性体内各要素に割り当てる。境界条件としては、全体磁界解析１２１と同じく、全体解析空間１１１の外周面ｓｓ’と内周面ｒｒ’のポテンシャルを０とおき、周方向端部の境界面ｓｒともう一方の境界面ｓ’ｒ’は互いに周期境界条件で接続される。また、外周側解析空間１１２側のスライド面ｂｂ’上、内周側解析空間１１３側のスライド面ｃｃ’上のポテンシャルは、外周面ｓｓ’及び内周面ｒｒ’と同様に０に設定することで、外周側解析空間１１２内および内周側解析空間１１３内で漏れ磁束のみを発生させるようにする。磁界のソースとしては、全体磁界解析１２１で与えた電流、磁石磁力を与え、線形磁界解析とし、全体磁界解析１２１で導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を実施すればよい。このとき電流、磁石磁力としては、全体磁界解析１２１で与えた電流、磁石磁力そのものを与えて１回の磁界解析を行っても良いし、総和が全体磁界解析１２１で与えた電流、磁石磁力となるように成分展開させた各電流、磁石磁力成分を与えて個別に磁界解析を行っても良い。
【００４１】
また、この固定子、回転子漏れ磁束の磁界解析１２３では、外周側解析空間１１２と内周側解析空間１１３の境界となるスライド面ｂｂ’上、スライド面ｃｃ’上の境界条件を０に設定するため、外周側解析空間１１２と内周側解析空間１１３とは、個別に解析することが可能になる。固定子、回転子漏れ磁束の磁界解析１２３とは別に、特定空間高調波、時間高調波磁界解析１２４を実施する。特定空間高調波、時間高調波磁界解析１２４では、固定子、回転子漏れ磁束の磁界解析２３と同様に、全体磁界解析１２１で得られた磁性体内各要素の透磁率データ１０５を、時刻ステップ毎に磁性体内各要素に割り当てる。境界条件としては、固定子、回転子漏れ磁束の磁界解析１２３と同様に、全体磁界解析１２１と同じく、全体解析空間１１１の外周面ｓｓ’と内周面ｒｒ’において境界条件を０に設定し、周方向端部の境界面ｓｒと境界面ｓ’ｒ’は互いに周期境界条件で接続される。外周側解析空間１２側のスライド面ｂｂ’上、内周側解析空間１１３側のスライド面ｃｃ’上には、固定子、回転子漏れ磁束の磁界解析１２３とは異なり、（９）式で展開したポテンシャルのうち、特定のｍ次空間高調波、１次時間高調波のみを抽出した、
【００４２】
【数１０】
2004245814.tif 000012
を境界条件として設定する。外周側解析空間１１２が固定子ならば、スライド面ｂｂ’の境界条件としての（２）式は静止座標系のものを用い、外周側解析空間１１２が回転子ならば、回転座標系のものを用いる。また同様にして、内周側解析空間１１３が固定子ならば、スライド面ｃｃ’の境界条件としての（１０）式は静止座標系のものを用い、内周側解析空間１１３が回転子ならば、回転座標系のものを用いる。こうした上で、全体磁界解析１２１で磁界ソースとして与えた電流・磁石磁力は与えずに、線形磁界解析を実行する。全体磁界解析１２１で導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を実施すればよい。
【００４３】
特定空間高調波、時間高調波磁界解析１２４では、外周側解析空間１１２と内周側解析空間１１３の境界となる、スライド面ｂｂ’上、スライド面ｃｃ’上にそれぞれ（１０）式の境界条件を与えるため、外周側解析空間１１２と内周側解析空間１１３は、個別に解析することができる。また、固定子、回転子漏れ磁束の磁界解析１２３と、特定空間高調波、時間高調波磁界解析１２４は、順不同で、任意に実施可能であり、特定空間高調波、時間高調波磁界解析１２４における特定のｍ次空間高調波、１次時間高調波は、解析者の欲する次数のみ着目して解析すればよい。
【００４４】
尚、ここで磁気ベクトルポテンシャルＡを用いた解析について説明したが、二次元空間の解析であれば、自動的にクーロンゲージ
【００４５】
【数１１】
2004245814.tif 000013
が成立するが、三次元空間の解析では一般に高速解析を行うためにゲージ（ｄｉｖＡ）を固定しないでゲージフリーあるいは部分的にゲージフリーで計算するのが好ましい。このため、時間ステップを追って磁気ベクトルポテンシャルＡを解析する場合、各時間ステップでｄｉｖＡがゆらぐため、磁気ベクトルポテンシャルＡが関する物理変数を時間高調波成分に展開する場合に問題が生じる。そこで、各時間ステップにおいて、求まった磁気ベクトルポテンシャルＡに対して、
【００４６】
【数１２】
2004245814.tif 000014
を用いてクーロンゲージを課した磁気ベクトルポテンシャルＡ’に変換するのが望ましい。この場合、スカラーポテンシャルΨは、
【００４７】
【数１３】
2004245814.tif 000015
を用いて求めることができる。このようにして求めた各節点上のスカラーポテンシャルΨに対して、
【００４８】
【数１４】
2004245814.tif 000016
に示した新たな変数ａ’を変数ａの代わりに用いれば、ゲージのゆらぎに関する問題は解消される。ここで、Ψ１およびΨ２は、それぞれ、変数ａの定義における磁気ベクトルポテンシャルＡの辺上の線積分における辺上の始点および終点に相当する節点上のスカラーポテンシャルΨの値である。
【００４９】
図８は、本発明の第７の実施の形態による磁界解析法の流れを示す図である。図８の実施例は、図７の実施例に対し、取り扱う高調波を空間高調波のみとしたものである。全体磁界解析１３１は図７と同じであり、磁性体内各要素の透磁率データ１０５も同じように格納する。ポテンシャルａの空間高調波分析１３２では、全体磁界解析１３１で得られたスライド面ｂｂ’上でのポテンシャルを、
【００５０】
【数１５】
2004245814.tif 000017
のように空間高調波の成分に展開する。固定子、回転子漏れ磁束の磁界解析１３３は、図７と同じである。特定空間高調波磁界解析１３４では、全体磁界解析１３１で得られた磁性体内各要素の透磁率データ１０５を、各時刻ステップ毎に磁性体内各要素に割り当てる。境界条件としては、全体解析空間１１１の外周面ｓｓ’と内周面ｒｒ’上のポテンシャルを図２と同じく０とおき、周方向端部の境界面ｓｒと境界面ｓ’ｒ’は互いに周期境界条件で接続される。外周側解析空間１１２側のスライド面ｂｂ’上、内周側解析空間１１３側のスライド面ｃｃ’上には、（１５）式で分析したポテンシャルのうち、特定のｍ次空間高調波（ｍ＝０，１，２、…）のみを抽出した、
【００５１】
【数１６】
2004245814.tif 000018
を境界条件として設定する。磁界のソースは，全体磁界解析１３１で与えた電流、磁石磁力を与えずに線形磁界解析とし、全体磁界解析１３１導電率を与えて渦電流を計算していれば同様に導電率を与え、渦電流を考慮した解析を実施すればよい。
【００５２】
図９は、本発明の第８の実施の形態による磁界解析法の流れを示す。図９に示す解析法は、図７に示す解析法に対して、取り扱う高調波を時間高調波のみとしたものである。全体磁界解析１４１は図７の場合と同じであり、磁性体内各要素の透磁率データ１０５も同じように格納する。ポテンシャルａの時間高調波分析１４２では、全体磁界解析１４１で得られたスライド面ｂｂ’上でのポテンシャルを、
【００５３】
【数１７】
2004245814.tif 000019
のように時間高調波の成分に展開する。固定子、回転子漏れ磁束の磁界解析１４３は、図７と同じである。特定時間高調波磁界解析１４４では、全体磁界解析１４１で得られた磁性体内各要素の透磁率データ１０５を、時刻ステップ毎に磁性体内各要素に割り当てる。境界条件としては、全体解析空間１１１の外周面ｓｓ’と内周面ｒｒ’のポテンシャルを図２と同じく０とおき、周方向端部の境界面ｓｒと境界面ｓ’ｒ’は互いに周期境界条件で接続される。外周側解析空間１１２側のスライド面ｂｂ’上、内周側解析空間１１３側のスライド面ｃｃ’上には、（１７）式で分析したポテンシャルのうち、特定の１次時間高調波のみを抽出した、
【００５４】
【数１８】
2004245814.tif 000020
【００５５】
を境界条件として設定する。全体磁界解析１４１で与えた磁界のソースとしてのコイル電流・磁石磁力を与えずに、線形磁界解析を実行する。全体磁界解析１４１で導電率を与えて渦電流を計算していれば、同様に導電率を与え、渦電流を考慮した解析を実施すればよい。
【００５６】
以上に説明したいずれの実施例においても、図１０に示すように、計算装置１５１で一連の回転子空間・固定子空間個別の解析で得られた磁界分布を、視覚的にわかるように、与えた境界条件別に個別に表示装置１５２に表示させることができる。
【００５７】
尚、以上に説明した解析法は、回転半径が実質的に無限大であるリニアモータにも適用できるのは言うまでもない。この場合には、例えば以下の手順により解析を行うことができる。
【００５８】
１）解析空間の外側境界面に境界条件を与え、電流又は磁石磁力による起磁力を与えてリニアモータ全体の磁界解析を実施する全体磁界解析手順。２）前記全体磁界解析手順で得られた磁性体内部の透磁率分布を記憶する手順。３）前記全体磁界解析手順で得られた、可動子と固定子との間のスライド面上における磁界分布を記述するポテンシャルを、移動方向と垂直な方向にはある分布を有し、移動方向に関して定数あるいは周期的な分布をもつ複数の局在化した高調波成分に展開する手順。４）前記全体磁界解析手順における境界条件と同じ境界条件を与え、さらに、前記スライド面上にはディリクレ条件を与え、前記透磁率分布を解析空間全体に与え、前記全体磁界解析により与えた電流又は磁石磁力を与えて、可動子を含む可動子空間と固定子を含む固定子空間とのそれぞれに線形の磁界解析を実施する手順。５）前記全体磁界解析における境界条件と同じ境界条件を与え、さらに、分析した前記高調波成分を順次前記スライド面上に境界条件として与え、さらに前記透磁率分布を磁性体全体に与えて固定し、電流又は磁石磁力を零にして、可動子を含む可動子空間と、固定子を含む固定子空間と、のそれぞれに線形の磁界解析を実施する手順。この場合において、さらに、上記一連の可動子空間と固定子空間と、のそれぞれ個別の解析により得られた磁界分布を、与えた前記境界条件別に個別に表示させる手順を行っても良い。
【００５９】
以上、本実施の形態に沿って説明したが、本発明はこれらの例に限定されるものではなく、種々の変形が可能であるのは言うまでもない。例えば、上記解析プログラムにおける手順に基づいて行われる解析方法又は解析装置も本発明の範疇に入ることは言うまでもない。
【００６０】
【発明の効果】
以上説明した本発明に関る回転機の磁界解析の方法によれば、特定の高調波成分に起因する、振動、騒音、鉄損、渦電流損などの諸問題の解決が迅速に実施できる。すなわち、特定の高調波成分に対する磁界分布や、漏れ磁束のみの磁界分布を得ることができるため、振動、騒音、鉄損、渦電流損などの寄与成分が特定でき、磁束流線から特定の高調波磁界を低減する方策が容易に策定できるため、高調波成分に起因する、振動、騒音、鉄損、渦電流損などの諸問題の考察、低減のための検討が大幅に加速される効果がある。また、２次元解析のみならず３次元解析にも適用可能であるため、より詳細な分析が可能になる。
【図面の簡単な説明】
【図１】本発明の第１の実施の形態による回転機の磁界解析法の流れを示す図である。
【図２】本発明の第２の実施の形態による磁界解析法の流れを示す図である。
【図３】本発明の実施の形態による磁界解析法の補足説明図である。
【図４】本発明の第３の実施の形態による磁界解析法の流れを示す図である。
【図５】本発明の第４の実施の形態による磁界解析法の流れを示す図である。
【図６】本発明の第５の実施の形態による磁界解析法の流れを示す図である。
【図７】本発明の第６の実施の形態による磁界解析法の流れを示す図である。
【図８】本発明の第７の実施の形態による磁界解析法の流れを示す図である。
【図９】本発明の第８の実施の形態である磁界解析法の流れを示す図である。
【図１０】本発明の回転機の磁界解析を実施する装置を示す図である。
【符号の説明】１…全体磁界解析、２…特定部位におけるポテンシャルの高調波成分への展開、３…成分磁界解析、４…任意成分磁界解析結果の重畳、５…各要素の透磁率データ、１１…全体解析空間、１２…外周側解析空間、１３…内周側解析空間、２１…全体磁界解析、２２…ポテンシャルａの空間高調波、時間高調波成分展開、２３…固定子、回転子漏れ磁束の磁界解析、２４…特定空間高調波、時間高調波磁界解析、２５…ポテンシャルａの空間高調波分析、２６…特定空間高調波磁界解析、２３…固定子、回転子漏れ磁束の磁界解析、１０１…全体磁界解析、１０２…高調波成分展開、１０３…成分磁界解析、１０３ａ…ソース成分磁界解析、１０３ｂ…成分磁界解析、１０５…透磁率データ。
【図１】
2004245814.tif 000021
【図２】
2004245814.tif 000022
【図３】
2004245814.tif 000023
【図４】
2004245814.tif 000024
【図５】
2004245814.tif 000025
【図６】
2004245814.tif 000026
【図７】
2004245814.tif 000027
【図８】
2004245814.tif 000028
【図９】
2004245814.tif 000029
【図１０】
2004245814.tif 000030