(57)【特許請求の範囲】
【請求項１】
ＰＭモータの特性評価時にｄｑ軸等価回路定数を電磁界解析によって算出するものにおいて、
前記ｄｑ軸等価回路定数のｄ軸磁束λdテーブルとｑ軸磁束λqテーブルを有する磁束テーブル、及び電流Ｉとインバータの制御角γをパラメータとしたトルクテーブを設け、与えられたトルク及び制御角γを用いて電流Ｉを求め、求められた電流Ｉと制御角γから前記磁束テーブルからｄ軸磁束λd及びｑ軸磁束λqを求めると共に、前記電流Ｉと制御角γから演算手段を介して前記等価回路のｄ軸電流ｉd’、ｑ軸電流ｉq’を算出し、これら軸磁束λd、λqと軸電流ｉd’、ｉq’及び電気角速度ωを用いて各軸等価回路の軸電圧Ｖd’とＶq’を求め、この軸電圧Ｖd’、Ｖq’と前記軸電流ｉd’、ｉq’から線間電圧実効値Ｖ1と相電流実効値Ｉ1を求めることを特徴としたＰＭモータの特性演算方法。
【請求項２】
前記電流Ｉ、制御角γ、電気角速度ωの条件で負荷をかけて電磁界解析により相電圧を求め、得られた相電圧を周波数波数分析して誘起電圧と同位相の成分Ｖq''と誘起電圧より位相が９０°進んだ成分Ｖd''とに分離し、次式により磁束λdとλqを得て平均トルクを求め、この手段を繰り返すことにより電流Ｉ、制御角γをパラメータとした前記トルクテーブルを作成することを特徴とした請求項１記載のＰＭモータの特性演算方法。
【数１】
0005145634.tif 000017
【請求項３】
前記算出された線間電圧実効値Ｖ1と予め決められた最大電圧値と比較し、線間電圧実効値Ｖ1＞最大電圧値のときに前記制御角γに任意の値である
△γを加算し、線間電圧実効値Ｖ1＜最大電圧値時の電圧値を弱め界磁制御時の線間電圧実効値とすることを特徴とした請求項１又は２記載のＰＭモータの特性演算方法。
【請求項４】
ＰＭモータの特性評価時にｄｑ軸等価回路定数を電磁界解析によって算出するものにおいて、
前記ｄｑ軸等価回路定数のｄ軸磁束λdテーブルとｑ軸磁束λqテーブルを有する磁束テーブル、及び電流Ｉとインバータの制御角γをパラメータとしたトルクテーブをコンピュータのメモリに設け、与えられたトルク及び制御角γを用いて電流Ｉを求める電流演算手段と、この電流演算手段によって求められた電流Ｉと制御角γから前記磁束からｄ軸磁束λd及びｑ軸磁束λqを求める軸磁束演算手段と、前記電流Ｉと制御角γから前記ｄｑ軸等価回路のｄ軸電流ｉd’、ｑ軸電流ｉq’を算出する軸電流演算手段と、前記軸磁束λd、λqと軸電流ｉd’、ｉq’及び電気角速度ωを用いて各軸等価回路の軸電圧Ｖd’とＶq’を求める軸電圧演算手段と、この軸電圧Ｖd’とＶq’と前記軸電流ｉd’、ｉq’から線間電圧実効値Ｖ1と相電流実効値Ｉ1を求める手段としてコンピュータを機能させるためのＰＭモータの特性演算用プログラム。
【請求項５】
前記電流Ｉ、制御角γ、電気角速度ωの条件で負荷をかけた電磁界解析により相電圧を求め、得られた相電圧を周波数波数分析して誘起電圧と同位相の成分Ｖq''と誘起電圧より位相が９０°進んだ成分Ｖd''とに分離した後、次式により磁束λdとλqを得て平均トルクを求め、この手段を繰り返すことにより電流Ｉ、制御角γをパラメータとした前記トルクテーブルを作成する手段としてコンピュータを機能させるための請求項４記載のＰＭモータの特性演算用プログラム。
【数２】
0005145634.tif 000018
【請求項６】
前記算出された線間電圧実効値Ｖ1と予め決められた最大電圧値と比較し、線間電圧実効値Ｖ1＞最大電圧値のときに前記制御角γに任意の値である
△γを加算し、線間電圧実効値Ｖ1＜最大電圧値時の電圧値を弱め界磁制御時の線間電圧実効値の演算手段としてコンピュータを機能させるための請求項４又は５記載のＰＭモータの特性演算用プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＰＭモータの特性演算方法とプログラムに関するものである。
【背景技術】
【０００２】
ＰＭモータ（永久磁石型同期電動機）は、一般産業用等としてインバータ駆動で使用される。ＰＭモータの特性及び制御法を評価するためにはｄｑ軸等価回路が用いられるが、特にモータ設計時には正確なｄｑ軸等価定数の把握ができれば、度重なるモータ試作、制御性評価回数の低減が可能となる。そのためには電磁界解析によるｄｑ軸等価回路定数を算出することが重要となっている。
近年、その電磁界解析によるモータ特性をより高速に、且つ高精度に把握することが求められている。この特性把握に用いられる電磁界解析は、非特許文献１や非特許文献２に開示された電磁界解析技術が用いられている。
【０００３】
図６は非特許文献２の図７．１２で示されているような一般的なＰＭモータ特性の計算に使用される等価回路を示したものである。電圧方程を用いて、電磁場解析によりモータ端子電圧を求め、等価回路定数を逆算する手法で等価回路を計算することにより、ＰＭモータのトルク及び回転数から電圧及び電流等を計算している。すなわち、図６の等価回路において、
【０００４】
【数３】
0005145634.tif 000002
【０００５】
【数４】
0005145634.tif 000003
【０００６】
【数５】
0005145634.tif 000004
【０００７】
【数６】
0005145634.tif 000005
【０００８】
ｉd'、ｉq'は各軸のトルク分電流、ｉmd、ｉmqは鉄損分電流であり、次式が成り立つ。
【０００９】
【数７】
0005145634.tif 000006
【００１０】
軸電圧Ｖd’とＶq’は次式で表される。
【００１１】
【数８】
0005145634.tif 000007
【００１２】
また、トルクは次式で表される。
【００１３】
【数９】
0005145634.tif 000008
【００１４】
線間電圧実効値及び相電流実効値は次式で表される。
【００１５】
【数１０】
0005145634.tif 000009
【００１６】
【数１１】
0005145634.tif 000010
【００１７】
その他の計算方法としては、磁界解析ソフトと回路方程式ソフトとを組み合わせた方法のものも上述した文献で報告されている。
【非特許文献１】電気学会技術報告書、第７７６号「回転機のバーチャルエンジニアリングのための電磁界解析技術」  ２０００年３月  電気学会発行
【非特許文献２】電気学会技術報告書、第８５５号「回転機の三次元ＣＡＥのための電磁界解析技術」  ２００１年９月  電気学会発行
【発明の開示】
【発明が解決しようとする課題】
【００１８】
図６で示す等価回路に基づく特性演算では、鉄心の磁気飽和によるｄ軸とｑ軸のインダクタンスＬｄ、Ｌｑ及びトルクの変動が考慮されていないためトルク−電流の線形性が保たれている領域では比較的精度のよい特性演算は可能であるが、図７で示すようにトルク−電流特性が非線形となる領域では、実測値に対して等価回路により計算した電圧値及びトルク値に誤差が生じる。この誤差を縮小するために電磁界解析と等価回路方式とを連携させたソフトウェアを用いて特性演算を実行すると、精度は向上するが、演算時間が長くなる問題を有している。なお、図７の（ａ）はトルク−電流特性図、（ｂ）は電圧−電流特性図である。
【００１９】
また、特性演算時にはＰＭモータの形状等を含む綿密な条件を必要とするが、企業のノウハウ等の存在によって企業間の定数等のデータ授受が円滑に実行できないことに基づき、高精度な特性演算が困難となっている。
【００２０】
そこで本発明が目的とするところは、特性演算の高速化及び高精度化を実現するＰＭモータの特性計算方法とプログラムを提供することにある。
【課題を解決するための手段】
【００２１】
本発明は、ＰＭモータの特性評価時にｄｑ軸等価回路定数を電磁界解析によって算出するものにおいて、
前記ｄｑ軸等価回路定数のｄ軸磁束λdテーブルとｑ軸磁束λqテーブルを有する磁束テーブル、及び電流Ｉとインバータの制御角γをパラメータとしたトルクテーブを設け、与えられたトルク及び制御角γを用いて電流Ｉを求め、求められた電流Ｉと制御角γから前記磁束テーブルからｄ軸磁束λd及びｑ軸磁束λqを求めると共に、前記電流Ｉと制御角γから演算手段を介して前記等価回路のｄ軸電流ｉd’、ｑ軸電流ｉq’を算出し、これら軸磁束λd、λqと軸電流ｉd’、ｉq’及び電気角速度ωを用いて各軸等価回路の軸電圧Ｖd’とＶq’を求め、この軸電圧Ｖd’、Ｖq’と前記軸電流ｉd’、ｉq’から線間電圧実効値Ｖ1と相電流実効値Ｉ1を求めることを特徴としたものである。
【００２２】
請求項２に記載された本発明は、前記電流Ｉ、制御角γ、電気角速度ωの条件で負荷をかけて電磁界解析により相電圧を求め、得られた相電圧を周波数波数分析して誘起電圧と同位相の成分Ｖq''と誘起電圧より位相が９０°進んだ成分Ｖd''とに分離した後、磁束λdとλqを得て平均トルクを求め、この手段を繰り返すことにより電流Ｉ、制御角γをパラメータとした前記トルクテーブルを作成することを特徴としたものである。
【００２３】
請求項３に記載された本発明は、前記算出された線間電圧実効値Ｖ1と予め決められた最大電圧値と比較し、線間電圧実効値Ｖ1＞最大電圧値のときに前記制御角γに任意の値である△γを加算し、線間電圧実効値Ｖ1＜最大電圧値時の電圧値を弱め界磁制御時の線間電圧実効値とすることを特徴としたものである。
【００２４】
請求項４に記載された本発明は、ＰＭモータの特性評価時にｄｑ軸等価回路定数を電磁界解析によって算出するものにおいて、
前記ｄｑ軸等価回路定数のｄ軸磁束λdテーブルとｑ軸磁束λqテーブルを有する磁束テーブル、及び電流Ｉとインバータの制御角γをパラメータとしたトルクテーブをコンピュータのメモリに設け、与えられたトルク及び制御角γを用いて電流Ｉを求める電流演算手段と、この電流演算手段によって求められた電流Ｉと制御角γから前記磁束からｄ軸磁束λd及びｑ軸磁束λqを求める軸磁束演算手段と、前記電流Ｉと制御角γから前記ｄｑ軸等価回路のｄ軸電流ｉd’、ｑ軸電流ｉq’を算出する軸電流演算手段と、前記軸磁束λd、λqと軸電流ｉd’、ｉq’及び電気角速度ωを用いて各軸等価回路の軸電圧Ｖd’とＶq’を求める軸電圧演算手段と、この軸電圧Ｖd’とＶq’と前記軸電流ｉd’、ｉq’から線間電圧実効値Ｖ1と相電流実効値Ｉ1を求める手段としてコンピュータを機能させることを特徴としたものである。
【００２５】
請求項５に記載された本発明は、前記電流Ｉ、制御角γ、電気角速度ωの条件で負荷をかけた電磁界解析により相電圧を求め、得られた相電圧を周波数波数分析して誘起電圧と同位相の成分Ｖq''と誘起電圧より位相が９０°進んだ成分Ｖd''とに分離した後、磁束λdとλqを得て平均トルクを求め、この手段を繰り返すことにより電流Ｉ、制御角γをパラメータとした前記トルクテーブルを作成する手段としてコンピュータを機能させることを特徴としたものである。
【００２６】
請求項６に記載された本発明は、前記算出された線間電圧実効値Ｖ1と予め決められた最大電圧値と比較し、線間電圧実効値Ｖ1＞最大電圧値のときに前記制御角γに任意の値である△γを加算し、線間電圧実効値Ｖ1＜最大電圧値時の電圧値を弱め界磁制御時の線間電圧実効値の演算手段としてコンピュータを機能させることを特徴としたものである。
【発明の効果】
【００２７】
以上のとおり、本発明によれば、トルクテーブル及び磁束テーブルを用意し、任意のトルク及びインバータ制御角を用いてトルクテーブルを補間しながら電流を求めてから等価回路の軸電流、軸電圧を求め、その後に線間電圧、相電流の各実行値を求めるようにしたものである。このため、従来の等価回路による演算よりも正確となり、且つ磁界解析法と等価回路方式の連携解析よりも高速演算が可能となるものである。また、弱め界磁制御時のモータ特性も容易に求めることが出来るものである。
【発明を実施するための最良の形態】
【００２８】
図１は、本発明の実施例を示す説明図で、１はＰＭモータにおけるｄ軸等価回路、２はｑ軸等価回路、３は変換部で各軸の電流ｉd’ｉq’と各軸磁束λd、λq
（Ｉ、γ）相互間の変換を行う。４は電流Ｉとインバータの制御角γとの関係が格納された磁束テーブルで、この磁束テーブル４はｄ軸磁束λd（Ｉ、γ）用４ａとｑ軸磁束λq（Ｉ、γ）用４ｂを有している。５はトルクテーブルである。
ここで、ｄ軸電流ｉd’及びｑ軸電流ｉq’とＩ、γとの関係は次式で表される。
【００２９】
【数１２】
0005145634.tif 000011
【００３０】
【数１３】
0005145634.tif 000012
【００３１】
鉄損抵抗ｉmd及びｉmqは次式で表される。
【００３２】
【数１４】
0005145634.tif 000013
【００３３】
また、軸電圧Ｖd’とＶq’は次式で表される。
【００３４】
【数１５】
0005145634.tif 000014
【００３５】
図２はトルクテーブル５の例を示したもので、モータ特性を計算する場合にはこのテーブルを補間しながら、目的とする電流Ｉ、インバータ制御角γの時の磁束及びトルク値を求める。補間は当該値の中間値を等価回路の両端より計算して求める。
【００３６】
図３は上述に基づいてＰＭモータの特性を算出するためのフローチャートを示したものである。ステップＳ１で、演算装置に対して所望する任意のトルクＴと制御角γを与えると、制御部はトルク値Ｔと制御角γを用いてトルクテーブル５を参照し、且つこのデーブルを補間演算しながら電流Ｉを求める。変換部３は、ｄ軸用磁束テーブル４ａ及びｑ軸用磁束テーブル４ｂを参照し、求められた電流Ｉと与えられた制御角γに対応したｄ軸磁束λdとｑ軸用磁束λqを求めると共に、電流Ｉと制御角γを用いて１０式に基づいてｄ軸等価回路の電流ｉd’とｑ軸等価回路の電流ｉq’を演算する（ステップＳ２）。
【００３７】
ステップＳ３では、磁束λd、λqと電流ｉd’、ｉq’と電気角速度ωを用いて１３式の演算を実行して軸電圧Ｖd’とＶq’を求め、求められた軸電圧を用いて８式の演算をおこなって線間電圧実行値を求める。また、ステップＳ２で求まった軸電流ｉd’、ｉq’を用いて９式の演算を実行して相電流実効値を求める。
【００３８】
この実施例によれば、トルクテーブルを設け、このテーブルを補間しながら電流を求めるようにしたものであるから、ＰＭモータの特性を等価回路による計算よりも正確であり、磁界解析と等価回路方式との連携解析よりも高速に計算することが可能となる。
【００３９】
図４はトルクテーブル５を電磁界解析の手法によって求める場合の例を示したものである。
【００４０】
図１で示すように電流I、制御角γ、及び電気角速度ωの条件で負荷をかけた場合の電磁界解析を行って相電圧を求める。得られた相電圧の周波数波数分析を行ってsin(誘起電圧と同位相の成分)とcos(誘起電圧より９０°位相が進んだ成分Ｖd'')に分離する。そして、図４で示すように誘起電圧（相電圧波高値）と同位相の成分をＶq''とし、誘起電圧より９０°位相が進んだ成分をＶd''として次式によりλdとλqを求めた後、電磁界解析によりトルクを求める。
【００４１】
【数１６】
0005145634.tif 000015
【００４２】
【数１７】
0005145634.tif 000016
【００４３】
次いで、電流I、制御角γをパラメータとし、上記を繰り返すことによって図２で示すトルクテーブルを作成する。したがって、この手法を用いることによって、図１に示す磁束λdとλqのテーブルを電磁界解析の結果により容易に求めることが可能となる。
【００４４】
図５は弱め界磁制御されるPMモータの特性演算方法を示したフローチャートである。この場合も図１で示した磁界テーブルとトルクテーブルを用いた等価回路に基づき算出されるが、ここでは、制御角γに初期値を代入し、線間電圧実効値Ｖ1が最大の電圧値になるまで制御角γに任意の値Δγを加えて繰り返し実行される。Δγとしては、例えばテーブルの値が１０度間隔で、位相の精度が１０度毎でよい場合にはΔγ＝１０となり、例えば、１度毎の精度が欲しい場合にはΔγ＝１となる。このΔγが大きいと電圧精度は悪くなるが、計算時間が速くなり、Δγが小さい場合には、電圧精度は良くなるが、計算時間が遅くなる。
また、Δγの初期値は、弱め界磁制御のみを行うときには基本的に０となり、最大トルク制御を行う場合には、トルクテーブルを使用し、要求するトルクを最も小さい電流で発生する制御角を求めて初期値とする。
【００４５】
図５において、ステップＳ１１でトルクＴと制御角γを与えると、制御部はトルク値Ｔと制御角γを用いてトルクテーブル５を参照し、且つこのデーブルを補間演算しながら電流Ｉを求める。変換部３は、ｄ軸用磁束テーブル４ａ及びｑ軸用磁束テーブル４ｂを参照し、求められた電流Ｉと制御角γに対応したｄ軸磁束λdとｑ軸用磁束λqを求めると共に、電流Ｉと制御角γを用いて１０式に基づいてｄ軸等価回路の電流ｉd’とｑ軸等価回路の電流ｉq’を演算する（ステップＳ１２）。
【００４６】
ステップＳ１３では、磁束λd、λqと電流ｉd’、ｉq’と電気角速度ωを用いて１３式の演算を実行して軸電圧Ｖd’とＶq’を求め、求められた軸電圧を用いて８式の演算をおこなって線間電圧実行値Ｖ1を求める。また、ステップＳ１２で求まった軸電流ｉd’、ｉq’を用いて９式の演算を実行して相電流実効値Ｉ1を求める。ステップ１４では、求まった電圧Ｖ1が予め定められた最大電圧値と比較され、最大電圧値より線間電圧実行値Ｖ1が大きい場合にはＳ１０において任意の制御角Δγ分をγに加算し、再度Ｓ１１〜Ｓ１４を実行してこれをＶ1＜最大電圧となるまで繰り返す。Ｖ1＜最大電圧となったときの線間電圧実行値Ｖ1とＳ１３で求まった線間電圧実行値Ｖ1が当該弱め界磁時のモータ特性となる。
【００４７】
したがって、この実施例によれば、弱め界磁時のＰＭモータの特性を、従来の等価回路による計算よりも正確となり、且つ界磁解析と等価回路方式の連携解析よりも高速に演算可能となるものである。
【図面の簡単な説明】
【００４８】
【図１】本発明の実施形態を示す磁束λとトルクＴテーブルを用いた等価回路。
【図２】本発明に使用されるトルクテーブルの例図。
【図３】本発明のモータ特性演算のためのフローチャート。
【図４】ｄ軸磁束λd、ｑ軸磁束λqの計算説明図。
【図５】本発明の弱め界磁時のモータ特性演算のためのフローチャート。
【図６】ＰＭモータの等価回路図。
【図７】トルク、電圧特性図で、（ａ）はトルク−電流図、（ｂ）は電圧−電流図。
【符号の説明】
【００４９】
１…  ｄ軸等価回路
２…  ｑ軸等価回路
３…  変換部
４…  磁束テーブル
５…  トルクテーブル
【図１】
0005145634.tif 000019
【図２】
0005145634.tif 000020
【図３】
0005145634.tif 000021
【図４】
0005145634.tif 000022
【図５】
0005145634.tif 000023
【図６】
0005145634.tif 000024
【図７】
0005145634.tif 000025