(57)【要約】【課題】定常時と過渡時のモータ動作の把握を高精度に行うことのできるモータとその駆動回路のシミュレーション方法を提供する。【解決手段】モータとその駆動回路を回路方程式で記述し、電磁界解析により算出した速度起電力係数とトルク係数、インダクタンス列を利用して、定常時と過渡時のモータの動作解析を行うように構成した。【選択図】        図１
2004236392.tif 000002
【特許請求の範囲】
【請求項１】
モータとその駆動回路を回路方程式で記述し、電磁界解析にて算出した瞬時の速度起電力係数とトルク係数を利用して、定常時と過渡時のモータの動作解析を行うことを特徴とするモータとその駆動回路のシミュレーション方法。
【請求項２】
モータとその駆動回路を回路方程式で記述し、電磁界解析にて算出したインダクタンス列を利用して、定常時と過渡時のモータの動作解析を行うことを特徴とするモータとその駆動回路のシミュレーション方法。
【請求項３】
モータとその駆動回路を回路方程式で記述し、電磁界解析にて算出した瞬時の速度起電力係数とトルク係数、および、インダクタンス列を利用して、定常時と過渡時のモータの動作解析を行うことを特徴とするモータとその駆動回路のシミュレーション方法。
【請求項４】
請求項１ないし３において、前記モータとその通電状況にそれぞれ対応する複数の回路モデル（回路方程式）を備え、前記モータに応じて使用する回路モデルを択一的に切り替えて実行することを特徴とするモータとその駆動回路のシミュレーション方法。
【請求項５】
請求項１ないし４において、算出した瞬時電流の最大値から前記モータで発生する鉄損のシミュレーションを行うことを特徴とするモータとその駆動回路のシミュレーション方法。
【請求項６】
請求項１ないし５において、コンデンサ容量と抵抗値、交流電源周波数、交流電源の投入位相に係るデータに基づいて、交流整流とその後段に挿入されるコンデンサを考慮したシミュレーションを可能としたことを特徴とするモータとその駆動回路のシミュレーション方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、電磁界解析を用いて回路方程式に必要なパラメータのみを抽出することにより、より精度が高く、製品設計過程における試作品製作回数の低減が可能なモータと、その駆動回路のシミュレーション方法に関する。
【０００２】
【従来の技術】
従来から、モータを解析する方法としては種々の方式が提案されている。例えば、モータとモータ駆動回路をともに伝達関数で指定する方法や、モータとモータ駆動回路をともに回路モデル（回路方程式）で指定する方法などがある。
【０００３】
しかし、前述したように、モータとモータ駆動回路の両方を伝達関数で指定した場合、系の細かい振る舞いの把握は困難であり、また、前記モータとモータ駆動回路をともに回路モデルで指定した場合は、系の細かい振る舞いの把握は可能となるものの、シミュレーション（解析）に要する時間が増大してしまう問題があった。
【０００４】
そこで、このような問題を解決するために、モータを伝達関数で記述し、モータ駆動回路を回路モデル（回路方程式）で記述する方法や、モータの半分を伝達関数で記述し、前記モータの残りの半分とモータ駆動回路を回路モデル（回路方程式）で記述する方法が提案されている（特許文献１参照）。
【０００５】
この方法によれば、前述した問題点、すなわち、系の振る舞いの把握が困難であるといった問題や、シミュレーション時間が増大するといった問題を確実に解消できるとともに、モータ駆動回路の細かい設計が容易となる等の効果が得られる。
【０００６】
【特許文献１】
特開平６−２３３５７４号（第３頁、第１図）
【０００７】
【発明が解決しようとする課題】
然るに、モータを伝達関数で記述し、モータ駆動回路を回路モデルで記述した場合、モータの定常時の動作を把握することは可能であるが、過渡時の動作を把握することはできなかった。
【０００８】
また、モータの半分を伝達関数で記述し、前記モータの残りの半分とモータ駆動回路を回路モデルで記述した場合、定常時のみならず過渡時のモータの動作を把握することができるが、この場合、両者を結びつけるためのリンク手段を予め用意する必要があり、シミュレーションプログラムが複雑となるなど、あまり好ましくなかった。
【０００９】
さらに、前述した２つの方法においては、いづれも、解析に必要なモータの速度起電力係数とトルク係数を、種々の正弦波の合成関数を複数用意することにより与えているため、前記合成関数から速度起電力係数列の波形を求めた場合、速度起電力係数やトルク係数に含まれる歪み（高調波成分）をシミュレーション波形に反映させることが非常に困難であるといった問題を有していた。
【００１０】
しかも、前記合成関数から速度起電力係数とトルク係数を与える場合、一般的には、基本波だけでなく高次の波形も与える必要があるが、その次数の与え方が明示されていなかった。
【００１１】
そこで、本発明は、定常時と過渡時の両方のモータの動作解析を、短時間かつ高精度で行うことのできるモータとその駆動回路のシミュレーション方法を提供する。
【００１２】
【課題を解決するための手段】
請求項１記載のモータとその駆動回路のシミュレーション方法は、モータとモータ駆動回路を回路方程式で記述し、電磁界解析によって算出した瞬時の速度起電力係数とトルク係数を用いて、定常状態および過渡状態のモータのシミュレーションを行うことを特徴とする。
【００１３】
請求項２記載のモータとその駆動回路のシミュレーション方法は、モータとその駆動回路を回路方程式で記述し、電磁界解析によって算出したインダクタンス列を利用して、定常状態および過渡状態のモータのシミュレーションを行うことを特徴とする。
【００１４】
請求項３記載のモータとその駆動回路のシミュレーション方法は、モータとその駆動回路を回路方程式で記述し、電磁界解析によって算出した瞬時の速度起電力係数とトルク係数、および、インダクタンス列を利用して、定常状態および過渡状態のモータのシミュレーションを行うことを特徴とする。
【００１５】
請求項４記載のモータとその駆動回路のシミュレーション方法は、請求項１ないし３記載の前記モータとその通電状況にそれぞれ対応した複数の回路モデル（回路方程式）を備え、前記モータに応じて使用する回路モデルを択一的に切り替えて実行することを特徴とする。
【００１６】
請求項５記載のモータとその駆動回路のシミュレーション方法は、請求項１ないし４において、算出した瞬時電流の最大値から前記モータで発生する鉄損のシミュレーションを行うことを特徴とする。
【００１７】
請求項６記載のモータとその駆動回路のシミュレーション方法は、請求項１ないし５において、コンデンサ容量，抵抗値，交流電源周波数，交流電源の投入位相に関するデータに基づいて、交流整流とその後段のコンデンサを考慮したシミュレーションを可能としたことを特徴とする。
【００１８】
本発明のモータとその駆動回路のシミュレーション方法によれば、速度起電力係数とトルク係数を電磁界解析によって直接算出することができるので、前記速度起電力係数とトルク係数の高調波成分がシミュレーション波形に与える影響を把握するのに非常に有利である。
【００１９】
また、本発明のモータとその駆動回路のシミュレーション方法によれば、シミュレータに複雑なリンク手段を備えることなく、定常時および過渡時のモータ動作を正確に把握することができる。
【００２０】
さらに、本発明のモータとその駆動回路のシミュレーション方法によれば、解析の対象となるモータ毎に特化した回路モデル（回路方程式）を予め複数備えておき、解析時に複数の回路モデルの中から１つを択一的に選択してモータのシミュレーションを行うので、モータとモータ駆動回路をともに回路モデルで記述しても、シミュレーションに長時間を要することを確実に防止することができる。
【００２１】
【発明の実施の形態】
以下、本発明の実施の形態を図１ないし図１０により説明する。図１は本発明のモータとその駆動回路のシミュレーション方法において、シミュレーションの対象となるモータ毎に特化して生成した専用のシミュレータの構成要素を記載した説明図である。
【００２２】
まず、図１記載のモータとしては、本実施例においては、図２に示す三相永久磁石同期モータを例示する。前記三相永久磁石同期モータは、そのモータ本体１に、インバータ回路（３相インバータ回路）２と、該インバータ回路２を制御する制御装置３と、前記インバータ回路２の入力として電源（直流または交流）４を接続して概略構成されている。
【００２３】
図３は、前記三相永久磁石同期モータの回路モデルを示しており（図１記載のモデル化）、抵抗，自己インダクタンスおよび相互インダクタンス，速度起電力係数とモータ回転子の回転数の積に比例する速度起電力を等価に置き換えた電源を一相分の構成として、それが３つ集まって構成されている。
【００２４】
なお、前記自己インダクタンスおよび相互インダクタンスは、固定または回転子位置変動のどちらかを、また、前記トルク係数および速度起電力係数は、平均値または変動のどちらかを選択できるものとする。
【００２５】
次に、下記［数１］に図１記載の回路方程式を記述する。
【００２６】
【数１】
2004236392.tif 000003
【００２７】
この回路方程式は、前記三相永久磁石同期モータの三相のうちニ相に通電する場合について記述されている。また、前記回路方程式は、シミュレーションの対象となるモータ毎にその記述内容が異なり、複数種類のモータのシミュレーションを考慮した場合、各モータ毎に予めこれをモデル化して各々の回路方程式を作成する必要がある。
【００２８】
つまり、前記シミュレーションを行うモータが異なれば、複数用意した回路方程式の中から適用する回路方程式を１つ選択することにより、それぞれのモータに特化した専用のシミュレータを構築することができ、この結果、シミュレーションに要する時間を短縮することが可能となる。
【００２９】
また、図１に記す速度起電力係数とトルク係数は電磁界解析によって算出されるものであり、一方、運動方程式としては、下記［数２］に記すように、例えば、ラグランジェの運動方程式が利用される。
【００３０】
【数２】
2004236392.tif 000004
【００３１】
図４は、本発明のシミュレーションの手順を簡略化して示す模式図である。該シミュレーションを実行する場合、まず、［数１］に記述した回路方程式で算出される電流と電磁界解析より求めたトルク係数を積算して、シミュレーションの対象であるモータ（本実施例では、三相永久磁石同期モータ）のトルクを算出する。
【００３２】
次に、［数２］に記述したラグランジェの運動方程式に代表される運動方程式に、算出したトルクを代入することによって前記モータの回転数を算出し、つづいて、前記回転数と電磁界解析によって求めた速度起電力係数との積から速度起電力を算出する。
【００３３】
そして、この速度起電力を１つのパラメータとして［数１］記載の回路方程式に入力（フィードバック）することによって前述した電流を再度計算する。以上の動作を繰り返し行うことによって、本発明のモータのシミュレーションは実行される。なお、前記シミュレーションによる結果は、電流波形や電圧波形などの各種波形データとして出力することが可能である。
【００３４】
下記［表１］に前記シミュレーションを実行するために入力が必要なパラメータ（定数）を示す。各パラメータの説明は以下のとおりである。
【００３５】
【表１】
2004236392.tif 000005
【００３６】
ＴＤ：計算の時間間隔（単位は（ｓ））である。
【００３７】
ｔＥ１ｏｆｆ：計算開始時（時刻０）には、電源電圧はＥ１で与えられる値になっている。時刻０からこの時間が経過すると、Ｅ１が０になる。
【００３８】
ＴＳＴＯＰ：計算開始の時刻０から計算終了までの時刻（単位は（ｓ））。
【００３９】
ｓｔｒｔｃｕｔ：計算結果の出力ファイルへの書き出しを始める時刻（単位は（ｓ））。
【００４０】
ｓｔｏｐｃｕｔ：計算結果の出力ファイルへの書き出しを終える時刻（単位は（ｓ））。
【００４１】
ｐｈａｓｅ：通電位相角（正は進み位相，負は遅れ位相（単位は（ｒａｄ）））。
【００４２】
ｓｔａｒｔ−ｒｐ：時刻０における回転子の位置（単位は（ｒａｄ））。
【００４３】
ｓｔａｒｔ−ｗ：時刻０における回転子の速度（単位は（ｍｉｎ−１））。
【００４４】
極数：モータの極数。
【００４５】
Ｅ１：インバータへ入力する直流電圧値（単位は（Ｖ））。
【００４６】
Ｊ：モータの慣性モーメント（単位は（ｋｇ・ｍ２））。
【００４７】
ＧＭ：回転数に比例する損失の比例定数（回転制動係数）。
【００４８】
ｋｔ０：負荷トルクを示す定数（負荷が一定の場合の単位は（Ｎ・ｍ）、負荷が回転数の２乗に比例する場合の単位は（Ｎ・ｍ・ｓ２／（ｒａｄ）２）。
【００４９】
ｋｔ０ｎ：負荷が回転数の何乗にあたるかを示す定数。
【００５０】
線間Ｒ：モータの線間の抵抗値（初期温度で与えられる値であり、単位は（Ω））。
【００５１】
線間Ｌ：モータの線間のインダクタンス値（単位は（Ｈ））。
【００５２】
ＶＤ：ダイオードのオン電圧（単位は（Ｖ））。
【００５３】
ｋｅ：うず電流計算用定数（単位は（Ｗ／Ａ２））。
【００５４】
ｋｈ：ヒステリシス損計算用定数（単位は（Ｗ／Ａ））。
【００５５】
ｗｅｉｇｈｔ：固定子１脚あたりの重量（単位（ｋｇ））。
【００５６】
ｎｗ：１巻線の巻数。
【００５７】
ｒａｄｉｕｓ：鉄損失計算用の固定子１脚の換算半径（単位は（ｍ））。
【００５８】
Ｂ０：固定子材質の残留磁束密度（単位は（Ｔ））。
【００５９】
ＨＭ：固定子材質が飽和し始める磁界の強さ（単位は（Ｈ））。
【００６０】
ＢＭ：固定子の飽和磁束密度（単位は（Ｔ））。
【００６１】
線間ｋｅｗ：速度起電力係数を正弦波近似する場合の実行値（ｋｔｉｕｋｅｗ＝２の場合のみ有効であり、単位は（Ｖ・Ｓ／ｒａｄ））。
【００６２】
線間ｋｔｉ：トルク係数を正弦波近似する場合の実行値（ｋｔｉｕｋｅｗ＝２の場合のみ有効であり、単位は（Ｎ・ｍ／Ａ））。
【００６３】
ｋｔｉｕｋｅｗ：速度起電力係数とトルク係数の与え方を示す。１は外部ファイルからであり、２は内部で定数（線間ｋｅｗと線間ｋｔｉ）により計算する。
【００６４】
初期温度：巻線抵抗の初期温度（単位は（℃））。
【００６５】
使用温度：巻線抵抗の使用時の温度（単位は（℃））。
【００６６】
ＰＷＭｏｎ／ｏｆｆ：ＰＷＭを行うか行わないかを示しており、１は行う。２は行わない。
【００６７】
ＰＷＭｏｎ時間：  ＰＷＭの周期のうち通電している時間（単位は（ｓ））。
【００６８】
ＰＷＭ周期：ＰＷＭの周期（単位は（ｓ））。
【００６９】
ｍｄ１３５ｐｗｍ：ＰＷＭをかける場合、通電モードを１，３，５とした場合に上下どの素子を入切するかを示す。１は上下とも入り。２は上切，下入。３は上入，下切。
【００７０】
ｍｄ２４６ｐｗｍ：ＰＷＭをかける場合、通電モード２，４，６とした場合に上下どの素子を入切するかを示す。１は上下とも入り。２は上切，下入。３は上入，下切。
【００７１】
Ｌ変動：インダクタンス変動方法を指定するもの（０→ｏｆｆ，１→ｏｎ（ｓｉｎ），２→ｏｎ（ａｒｙ））。
０…線間Ｌの値のみを読み込み、他の値は無視する。Ｌａｍｐ，線間平均Ｍ，線間変動Ｍの値は無視される。
１…（線間Ｌ＋線間平均Ｍ）を一定値とし、（Ｌａｍｐ＋線間変動Ｍ）が正弦波状に変化するとして与える。
２…インダクタンス列を外部ファイル（デフォルトはｌａｒｙ．ｃｓｖ）から読み込む。Ｌ変動，Ｌａｍｐ，線間平均Ｍ，線間変動Ｍの値は無視される。角度のみを考慮したインダクタンス列となる。この場合のｌａｒｙ．ｃｓｖの例を下記の［表２］に示す。
【００７２】
【表２】
2004236392.tif 000006
【００７３】
３…インダクタンス列を外部ファイル（デフォルトはｌａｒｙ．ｃｓｖ）から読み込む。Ｌ変動，Ｌａｍｐ，線間平均Ｍ，線間変動Ｍの値は無視される。角度と電流を考慮したインダクタンス列となる。この場合のｌａｒｙ．ｃｓｖの例を下記の［表３］に示す。
【００７４】
【表３】
2004236392.tif 000007
【００７５】
Ｌａｍｐ（単位は（Ｈ））：Ｌ変動が１のときに有効な値であり、自己インダクタンスは下記［数３］にて与えられる。
【００７６】
【数３】
2004236392.tif 000008
【００７７】
［数３］において、θは回転子角度，ｓｉｎは正弦波である。また、自己インダクタンスの一定値分Ｌ０を線間Ｌとして与え、自己インダクタンスの変動分ＭＡをＬａｍｐとして与える。
【００７８】
線間平均Ｍ（単位は（Ｈ）），線間変動Ｍ（単位は（Ｈ））：Ｌ変動が１のときに有効な値である。相互インダクタンスは下記［数４］で与えられものとする。
【００７９】
【数４】
2004236392.tif 000009
【００８０】
［数４］において、θは回転子角度，ｓｉｎは正弦波である。また、相互インダクタンスの一定値分Ｍ０を線間平均Ｍとして与え、相互インダクタンスの変動分ＭＡを線間変動として与える。
【００８１】
電源選択：電源を指定するものであり、１を選択した場合はＥ１の値を直流電源とし、２を選択した場合はＥ１の値を交流正弦波の実行値とする。
【００８２】
コンデンサ容量：単位は（Ｆ）
抵抗：単位は（Ω）
交流電源周波数：単位は（Ｈｚ）
交流電源投入位相：単位は（電気ｄｅｇｒｅｅ）
上記パラメータは、前記電源選択で２を選択したときに有効であり、電源はコンデンサインプット型とする。図５に電源選択で２を選択したときの概念図を示す。
【００８３】
Ｌ＋ｋ・Ｍ選択（Ｌ＋Ｍ（１），Ｌ＋２Ｍ（２））：電圧平衡方程式中でＬ＋ＭとＬ＋２Ｍの係数がある。簡略化して扱うためにどちらかに統一する。
１…Ｌ＋Ｍに統一する。
２…Ｌ＋２Ｍに統一する。
【００８４】
以上記載したパラメータを入力することによって、図４に示すシミュレーションは実行される。また、図４に示す速度起電力係数を電磁界解析によって算出した結果を図６に示す。
【００８５】
図６（ａ）は埋め込み磁石型同期電動機（ＩＰＭ）の速度起電力係数をグラフ化したものであり、同図（ｂ）は表面磁石型同期電動機（ＳＰＭ）の速度起電力係数をグラフ化したものである。
【００８６】
［従来の技術］欄で説明した正弦波の合成関数から速度起電力係数を求めた場合、図６に示すような波形を得ることは困難であり、本発明の優位性を示すグラフであるといえる。なお、トルク係数については、前記速度起電力係数と同一値にできることは公知である。
【００８７】
図７は、本発明に係る電磁界解析によって算出したインダクタンス列を示す一例であり、前述の［表３］と対応している。即ち、図７は角度と電流をパラメータとしたインダクタンス列を電磁界解析によって算出し、それをグラフ化したものであり、電流値によって前記インダクタンス値が異なることがわかる。
【００８８】
そして、前記シミュレーションを行う場合は、電磁界解析によって算出したインダクタンス（図７）をテーブル化して利用したり、前記インダクタンスを一定としてシミュレーションを行う。
【００８９】
図８は、前記シミュレーションの結果、出力される波形の一例である。同図（ａ）に電流−時間波形を、同図（ｂ）に電圧−時間波形を示している。図８に示す波形が、オシロスコープ等を用いて実測する場合と非常に近い波形イメージであることがわかる。
【００９０】
図９は電流−時間波形について、これを実測した場合を（ａ）に、また、電磁界解析によって算出した速度起電力係数とトルク係数を用いて本発明のシミュレーションを実行した場合を同図（ｂ）に、さらに、［従来の技術］欄に記載した正弦波の合成関数から速度起電力係数とトルク係数を求めた場合を同図（ｃ）に示している。
【００９１】
図９（ａ）〜（ｃ）を互いに比較すると、本発明の電磁界解析を利用したシミュレーションによる電流−時間波形が、正弦波の合成関数から速度起電力係数とトルク係数を求めた場合に比べて、より実測波形に近似していることは明らかである。
【００９２】
つまり、図９は、本発明のシミュレーション方法によれば、従来技術と比較して、より精度の高い（実測値に近い）シミュレーションが可能であることを示している。
【００９３】
つづいて、本発明のシミュレーションによってモータの鉄損を算出する場合について説明する。前記鉄損を算出する場合は、まず最初に、前記シミュレーションによって算出した電流のピーク値にシミュレーションの対象であるモータのコイル巻数を乗じることによって、磁界の値を求める。
【００９４】
次に、モータコア材（ケイ素鋼板や焼結フェライト等）のヒステリシスカーブ（本実施例では、図１０に直線で近似した簡略的なヒステリシスカーブを例示するが、高次方程式を用いて、より精度の高いヒステリシスカーブに置き換えられることは当然である）上でその磁界に対応する最大磁束密度を読み取る。この最大磁束密度の値が判明すれば、下記の［数５］によって鉄損を容易に算出することができる。
【００９５】
【数５】
2004236392.tif 000010
【００９６】
以上説明したように、本発明のモータとその駆動回路のシミュレーション方法によれば、定常時のみならず過渡時におけるモータ動作をも確実に把握することができるとともに、その精度は、従来の正弦波の合成関数を利用したシミュレーション方法に比較して、より精度の高い解析を可能とした。
【００９７】
【発明の効果】
モータと駆動回路の両方を回路方程式で記述することにより、モータの瞬時および過渡時の波形をともに把握することができる。
【００９８】
モータと駆動回路の両方を回路方程式で記述することにより、複雑なリンク手段を予め用意する必要は一切なく、シミュレーションのプログラムを簡易に構成することができる。
【００９９】
電磁界解析によって算出した速度起電力係数とトルク係数を利用してモータのシミュレーションを実行することにより、より精度の高いシミュレーションを実現することができる。
【０１００】
シミュレーションの対象となるモータ毎に回路モデル（回路方程式）を予め作成し、対象が変更される都度、複数の回路モデル（回路方程式）を切り替えて、対象となるモータに特化した回路モデルを選択するので、モータと駆動回路の両方を回路モデルで記述しているにもかかわらず、そのシミュレーション時間を大幅に短縮することができる。
【０１０１】
モータの鉄損を容易に算出することができるとともに、コンデンサ容量，抵抗値，交流電源周波数，交流電源の投入位相のデータを基に、交流整流とコンデンサを具備したモータのシミュレーションを実現可能とした。
【図面の簡単な説明】
【図１】本発明のシミュレーション方法を実現する要素を記載した説明図である。
【図２】前記シミュレーション方法を適用するモータの一例を示す構成図である。
【図３】前記モータの回路モデルを示す回路図である。
【図４】前記シミュレーションの手順を示す模式図である。
【図５】前記シミュレーションの入力パラメータで電源選択を２としたときの概念図である。
【図６】（ａ）は埋め込み磁石型同期電動機（ＩＰＭ）の速度起電力係数を示すグラフであり、（ｂ）は表面磁石型同期電動機（ＳＰＭ）の速度起電力係数を示すグラフである。
【図７】電磁界解析により算出したインダクタンス列を示すグラフである。
【図８】（ａ）は、前記シミュレーションの結果として出力される電流−時間の関係を示すグラフであり、（ｂ）は電圧−時間の関係を示すグラフである。
【図９】（ａ）は、電流−時間の関係についての実測値を示すグラフであり、（ｂ）は電磁界解析によって算出した速度起電力係数とトルク係数を用いて前記シミュレーションを行った場合の電流−時間の関係を示すグラフであり、（ｃ）は正弦波の合成関数から速度起電力係数とトルク係数を求めた場合の電流−時間の関係を示すグラフである。
【図１０】モータコア材のヒステリシスカーブである。
【符号の説明】
１  モータ本体
２  インバータ回路
３  制御回路
４  電源
【図１】
2004236392.tif 000011
【図２】
2004236392.tif 000012
【図３】
2004236392.tif 000013
【図４】
2004236392.tif 000014
【図５】
2004236392.tif 000015
【図６】
2004236392.tif 000016
【図７】
2004236392.tif 000017
【図８】
2004236392.tif 000018
【図９】
2004236392.tif 000019
【図１０】
2004236392.tif 000020