(57)【要約】【課題】電動機の低回転数から高回転数までの広い動作領域において永久磁石の温度の推定精度を向上させることができる同期電動機の磁石温度推定システムを提供する。【解決手段】同期電動機２を駆動する基本波の周波数とは異なる周波数の電圧または電流を少なくとも同期電動機２のｄ軸に重畳し、重畳した電圧または電流と、重畳することにより得られた電流または電圧から同期電動機２のインピーダンスを演算し、演算したインピーダンスに基づいて永久磁石の温度Ｔｍを推定する。【選択図】図１
2015133890.tif 000002
【特許請求の範囲】
【請求項１】
永久磁石を有する同期電動機の磁石温度推定システムにおいて、
前記同期電動機を駆動する基本波の周波数とは異なる周波数の電圧または電流を少なくとも前記同期電動機のｄ軸に重畳する重畳部と、
重畳した電圧または電流と、重畳することにより得られた電流または電圧から前記同期電動機のインピーダンスを演算する演算部と、
演算した前記インピーダンスに基づいて前記永久磁石の温度を推定する磁石温度推定部と
を備えることを特徴とする同期電動機の磁石温度推定システム。
【請求項２】
前記重畳部は、前記基本波の周波数とは異なる周波数の電圧または電流を前記同期電動機のｄ軸にのみ重畳することを特徴とする請求項１に記載の同期電動機の磁石温度推定システム。
【請求項３】
前記重畳部は、前記基本波の周波数とは異なる周波数の電圧または電流を所定の間隔で重畳することを特徴とする請求項１または２に記載の同期電動機の磁石温度推定システム。
【請求項４】
磁石温度推定部は、予め取得した前記インピーダンスと前記永久磁石の温度との関係と、演算した前記インピーダンスとに基づいて、前記永久磁石の温度を推定することを特徴とする請求項１〜３のいずれか１項に記載の同期電動機の磁石温度推定システム。
【請求項５】
前記同期電動機の固定子コイルの温度を計測する温度計測部をさらに備え、
前記重畳部が、前記基本波の周波数とは異なる周波数の電圧または電流を前記同期電動機に重畳した際に、前記演算部は、前記温度計測部で計測された前記固定子コイルの温度に応じて、前記インピーダンスに含まれる前記固定子コイルの抵抗値を補正し、前記磁石温度推定部は、補正後の前記固定子コイルの抵抗値を含む前記インピーダンスに基づいて前記永久磁石の温度を推定することを特徴とする請求項１〜４のいずれか１項に記載の同期電動機の磁石温度推定システム。
【請求項６】
前記演算部は、前記温度計測部で計測された前記固定子コイルの温度と、予め取得した前記固定子コイルの温度と前記固定子コイルの抵抗値との関係に基づいて、前記インピーダンスに含まれる前記固定子コイルの抵抗値を補正し、前記磁石温度推定部は、補正後の前記固定子コイルの抵抗値を含む前記インピーダンスに基づいて前記永久磁石の温度を推定することを特徴とする請求項５に記載の同期電動機の磁石温度推定システム。
【請求項７】
前記同期電動機は、電流型インバータで駆動され、
前記重畳部は、下記（１）式で表される電流を前記同期電動機に重畳し、
前記演算部は、下記（１）式で表される電流を前記同期電動機に重畳することにより得られた下記（２）式で表される電圧を用いて、下記（１）式と下記（２）式から求めた下記（３）式で表される前記インピーダンスを演算することを特徴とする請求項１〜６のいずれか１項に記載の同期電動機の磁石温度推定システム。
【数１】
2015133890.tif 000012
ここで、
ｉｄｓｃ＊：ｄ軸高調波電流指令値、
ｉｑｓｃ＊：ｑ軸高調波電流指令値、
Ｉｃ：ｄ軸高調波電流指令値の振幅、
ωｃ：ｄ軸高調波電流指令値の角周波数、
ｔ：は時間
である。
【数２】
2015133890.tif 000013
ここで、
ｖｄｓｃ：ｄ軸高調波電圧値、
ｖｑｓｃ：ｑ軸高調波電圧値、
Ｒｄ：コイル抵抗値と磁石抵抗値を含む回路抵抗値、
Ｌｄ：ｄ軸インダクタンス、
ωｃ：ｄ軸高調波電流指令値の角周波数、
である。
【数３】
2015133890.tif 000014
ここで
Ｚｄｓ：インピーダンス、
Ｒｄ：コイル抵抗値と磁石抵抗値を含む回路抵抗値、
Ｌｄ：ｄ軸インダクタンス、
ωｃ：ｄ軸高調波電流指令値の角周波数、
ｖｄｓｃ：ｄ軸高調波電圧値、
ψｚｄ：インピーダンスの相差角、
Ｉｃ：ｄ軸高調波電流指令値の振幅、
ｔ：時間
である。
【請求項８】
前記同期電動機は、電圧型インバータで駆動され、
前記重畳部は、下記（４）式で表される電圧を前記同期電動機に重畳し、
前記演算部は、下記（４）式で表される電圧を前記同期電動機に重畳することにより得られた下記（５）式で表される電流を用いて、下記（４）式と下記（５）式から求めた下記（６）式で表される前記インピーダンスを演算することを特徴とする請求項１〜６のいずれか１項に記載の同期電動機の磁石温度推定システム。
【数４】
2015133890.tif 000015
ここで、
ｖｄｓｃ＊：ｄ軸高調波電圧指令値、
Ｖｃ：ｄ軸高調波電圧指令値の振幅、
ωｃ：ｄ軸高調波電圧指令値の角周波数、
ｔ：時間
である。
【数５】
2015133890.tif 000016
ここで、
ｉｄｓｃ：ｄ軸高調波電流値、
ｉｑｓｃ：ｑ軸高調波電流値、
ｖｄｓｃ＊：ｄ軸高調波電圧指令値、
Ｒｄ：コイル抵抗値と磁石抵抗値を含む回路抵抗値、
Ｌｄ：ｄ軸インダクタンス、
ωｃ：ｄ軸高調波電圧指令値の角周波数
である。
【数６】
2015133890.tif 000017
ここで
Ｚｄｓ：インピーダンス、
Ｒｄ：コイル抵抗値と磁石抵抗値を含む回路抵抗値、
Ｌｄ：ｄ軸インダクタンス、
ωｃ：ｄ軸高調波電圧指令値の角周波数、
Ｖｃ：ｄ軸高調波電圧指令値の振幅、
ψｚｄ：インピーダンスの相差角、
ｉｄｓｃ：ｄ軸高調波電流値、
ｔ：時間
である。
【請求項９】
前記基本波と、前記基本波の周波数とは異なる周波数の電圧または電流を分離するバンドパスフィルター、及びバンドストップフィルターをさらに備え、
前記演算部は、分離した電圧または電流に基づいて前記インピーダンスを演算することを特徴とする請求項１〜８のいずれか１項に記載の同期電動機の磁石温度推定システム。
【請求項１０】
分離した電圧または電流から前記同期電動機の回転方向とは異なる側の周波数成分を抽出する第２のバンドストップフィルターをさらに備え、
前記演算部は、抽出した電圧または電流に基づいて前記インピーダンスを演算することを特徴とする請求項９に記載の同期電動機の磁石温度推定システム。
【請求項１１】
演算した前記インピーダンスを補償する補償演算部をさらに備え、
前記磁石温度推定部は、補償した前記インピーダンスに基づいて前記永久磁石の温度を推定することを特徴とする請求項１〜１０のいずれか１項に記載の同期電動機の磁石温度推定システム。
【請求項１２】
前記補償演算部は、予め取得したｄ軸基本波電流と演算した前記インピーダンスを補償するための補償量との関係、または予め取得したｑ軸基本波電流と演算した前記インピーダンスを補償するための補償量との関係のうち、少なくとも１つの関係を用いて演算した前記インピーダンスを補償することを特徴とする請求項１１に記載の同期電動機の磁石温度推定システム。
【請求項１３】
予め取得した前記インピーダンスと前記永久磁石の着磁量との関係と、演算した前記インピーダンスとに基づいて、前記永久磁石の着磁量を推定する着磁量推定部をさらに備えることを特徴とする請求項１〜１２のいずれか１項に記載の同期電動機の磁石温度推定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、同期電動機の磁石温度推定システムに関する。
【背景技術】
【０００２】
回転子に永久磁石を配置し、固定子により形成された回転磁界と永久磁石の相互作用により回転子が回転する永久磁石型同期電動機が知られている。永久磁石は、温度が上昇すると減磁現象によって磁束が減少する。そこで、永久磁石の温度を許容温度以下に抑制するように制御する必要がある。永久磁石の温度を検出するには、永久磁石に温度センサを設けなければならず、電動機の小型化の制約になる。そこで、温度センサを用いずに永久磁石の温度を推定する方法が特許文献１に開示されている。特許文献１では、電動機の誘起電圧から永久磁石の温度を推定している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００７−６６１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１に開示された方法では、誘起電圧が小さい低回転時には、誘起電圧の測定誤差が大きくなり、永久磁石の温度の推定精度が低くなるおそれがある。
【０００５】
本発明は、上記問題に鑑みてなされたものであり、その目的とするところは、電動機の低回転数から高回転数までの広い動作領域において永久磁石の温度の推定精度を向上させることができる同期電動機の磁石温度推定システムを提供することである。
【課題を解決するための手段】
【０００６】
本発明の一態様に係る同期電動機の磁石温度推定システムは、同期電動機を駆動する基本波の周波数とは異なる周波数の電圧または電流を少なくとも同期電動機のｄ軸に重畳し、重畳した電圧または電流と、重畳することにより得られた電流または電圧から同期電動機のインピーダンスを演算し、演算したインピーダンスに基づいて永久磁石の温度を推定する。
【発明の効果】
【０００７】
本発明によれば、電動機の低回転数から高回転数までの広い動作領域において永久磁石の温度の推定精度を向上させることができる。
【図面の簡単な説明】
【０００８】
【図１】図１は、本発明の第１の実施形態に係る永久磁石型同期電動機の制御装置のシステム構成図である。
【図２】図２は、本発明の第１の実施形態に係る永久磁石温度推定の原理を説明する図である。
【図３】図３は、本発明の第１の実施形態に係る磁石温度推定部の構成図である。
【図４】図４は、高調波インピーダンス実部と永久磁石温度との関係を説明するグラフである。
【図５】図５は、本発明の第１の実施形態に係る磁石温度推定部の変形例１の構成図である。
【図６】図６は、固定子コイル温度と固定子コイル抵抗値との関係を説明するグラフである。
【図７】図７は、本発明の第１の実施形態に係る磁石温度推定部の変形例２の構成図である。
【図８】図８は、本発明の第１の実施形態に係る磁石温度推定部の変形例３の構成図である。
【図９】図９は、ｄ軸基本波電流値と高調波インピーダンス実部を補償するための補償量との関係を説明するグラフである。
【図１０】図１０は、ｑ軸基本波電流値と高調波インピーダンス実部を補償するための補償量との関係を説明するグラフである。
【図１１】図１１は、本発明の第２の実施形態に係る永久磁石型同期電動機の制御装置のシステム構成図である。
【図１２】図１２は、本発明の第２の実施形態に係る磁石温度推定部の構成図である。
【図１３】図１３は、本発明の第３の実施形態に係る永久磁石型同期電動機の制御装置のシステム構成図である。
【図１４】図１４は、本発明の第３の実施形態に係る磁石着磁量推定部の構成図である。
【図１５】図１５は、永久磁石着磁量と高調波インピーダンス実部との関係を説明するグラフである。
【発明を実施するための形態】
【０００９】
以下、本発明の実施形態を図面に基づいて説明する。
【００１０】
［第１の実施形態］
図１は、本発明の第１の実施形態に係る永久磁石型同期電動機の制御装置のシステム構成図である。図１に示すように、永久磁石型同期電動機の制御装置１は、電動機２を制御するために、電流制御部４と、座標変換部６，１１と、電力変換部７と、バンドストップフィルター９と、バンドパスフィルター１０と、共振制御部１３と、磁石温度推定部１４とを備えている。なお、図面において、２本斜線、３本斜線は、それぞれ２次元、３次元のベクトルを示す。
【００１１】
電動機２は、図２に示すように、三相の永久磁石型同期電動機（ＰＭＳＭ：permanent magnet synchronous motor）であり、永久磁石３２を備えた回転子３１が固定子３０の内部を回転するように構成されている。固定子３０には複数のスロット３３が設けられており、このスロット３３に固定子コイル３４が巻回されている。この固定子コイル３４に三相電力を供給することにより電流磁束が発生し、永久磁石３２は磁石磁束を発生する。なお、温度センサ３５は、固定子コイル３４の温度を計測するものである。
【００１２】
次に、図２を用いて、本実施形態における永久磁石温度推定の原理について説明する。本実施形態では、高調波電圧Ｖｈを固定子コイル３４に重畳し、発生する高調波電流値Ｉｈによって変化する高調波インピーダンスＺｈの実部Ｒｄに基づいて永久磁石３２の温度推定を行う。固定子コイルの抵抗値をＲｃ、固定子コイルのインダクタンスをＬｃ、永久磁石３２の抵抗値をＲｍ、永久磁石３２のインダクタンスをＬｍとすると、電動機２の等価回路は、図２のように表すことができる。高調波電圧Ｖｈと、等価回路に高調波電圧Ｖｈを重畳することにより発生する高調波電流値Ｉｈから高調波インピーダンスＺｈを演算すると、Ｚｈ＝Ｖｈ／Ｉｈで表される。また、高調波インピーダンスＺｈの実部Ｒｄは、（１）式で表される。
【数１】
2015133890.tif 000003
ここで、Ｍは相互インダクタンス、ωは高調波電圧Ｖｈの角周波数である。また、永久磁石３２の抵抗値Ｒｍ及びインダクタンスＬｍは永久磁石３２の温度Ｔｍの関数である。したがって、永久磁石３２の温度Ｔｍが変化すると、永久磁石３２の抵抗値Ｒｍ及びインダクタンスＬｍの値が変化する。すなわち、高調波インピーダンスＺｈの実部Ｒｄが変化する。後述するように、高調波インピーダンスＺｈの実部Ｒｄと永久磁石３２の温度Ｔｍには、相関関係がある。したがって、高調波電圧Ｖｈを重畳することにより得られる高調波インピーダンスＺｈの実部Ｒｄに基づいて永久磁石３２の温度Ｔｍを推定することができる。
【００１３】
以下、順次、図１に基づいて、電動機２の制御装置１の各構成部の構成、機能及び動作について説明する。
差分器３は、トルク指令値に基づくｄ−ｑ軸上の電流指令値ｉｄｓｆ＊，ｉｑｓｆ＊から、電動機２に実際に流れている検出電流ｉｄｓ，ｉｑｓから高調波成分をカットした検出基本波電流ｉｄｓｆ，ｉｑｓｆをそれぞれ減算する。
【００１４】
電流制御部４は、電流指令値ｉｄｓｆ＊，ｉｑｓｆ＊と、検出基本波電流ｉｄｓｆ，ｉｑｓｆとの偏差がなくなるように比例積分制御して第一の電圧指令値ｖｄ０＊，ｖｑ０＊を出力する。
【００１５】
差分器１２は、外部から与えられる高調波電流指令値ｉｄｓｃ＊，ｉｑｓｃ＊から、電動機２に実際に流れている検出電流ｉｄｓ，ｉｑｓから基本波成分をカットした検出高調波電流値ｉｄｓｃ，ｉｑｓｃをそれぞれ減算する。高調波電流指令値ｉｄｓｃ＊，ｉｑｓｃ＊は（２）式で表される。本実施形態では、電動機２を駆動するための基本波に重畳する高調波信号を電流の指令値として与える。
【数２】
2015133890.tif 000004
ここで、Ｉｃはｄ軸高調波電流指令値ｉｄｓｃ＊の振幅、ωｃはｄ軸高調波電流指令値ｉｄｓｃ＊の角周波数、ｔは時間である。ｄ軸高調波電流指令値ｉｄｓｃ＊の周波数は、基本波と異なる周波数である。（２）式からわかるように、本実施形態では、電動機２のｄ軸成分に高調波電流を重畳する。また、ｄ軸高調波電流指令値ｉｄｓｃ＊の振幅Ｉｃは、基本波の振幅より小さい。
【００１６】
共振制御部１３は、差分器１２の出力の大きさに応じて、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊を生成する。共振制御部１３は、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊の振幅及び間隔を任意に設定することができる。共振制御部１３で生成される高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊は、（３）式で表される。本実施形態において、基本波に高調波電圧を重畳する際に、パルセイティング・ベクトル・インジェクション（Pulsating vector injection）方式を用いる。パルセイティング・ベクトル・インジェクション方式とは、基本波のｄ軸成分に対して、ｄ軸と平行する方向（正負の両方向）に高調波電圧を重畳する方式である。
【数３】
2015133890.tif 000005
ここで、Ｒｄは固定子のコイル抵抗値と永久磁石抵抗値を含む回路抵抗値である。なお、この回路抵抗値には、電動機２と電力変換部７との間の配線抵抗値が含まれる。Ｌｄはｄ軸インダクタンス、ωｒは回転子の角周波数である。本実施形態では、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊のうち、ｑ軸成分をカットしてｄ軸成分のみを電動機２のｄ軸成分に重畳する。ｑ軸成分をカットする理由は、（３）式において、ｑ軸高調波電圧指令値ｖｑｓｃ＊にωｒが含まれているからである。このωｒは回転子の角周波数であるため、ｑ軸高調波電圧指令値ｖｑｓｃ＊は、回転子の角周波数によって値が変化する。このような回転子の角周波数による影響を防ぐため、ｑ軸成分をカットする。また、ｑ軸成分をカットすることにより、ｑ軸インダクタンスによる影響を防ぐこともできる。
【００１７】
加算器５は、第一の電圧指令値ｖｄ０＊，ｖｑ０＊に、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊を重畳して、第二の電圧指令値ｖｄｓ＊，ｖｑｓ＊を出力する。
【００１８】
座標変換部６は、加算器５の出力である第二の電圧指令値ｖｄｓ＊，ｖｑｓ＊を座標変換して、三相電圧指令値ｖｕ＊，ｖｖ＊，ｖｗ＊を出力する。
【００１９】
電力変換部７は、例えばコンバータとインバータで構成される電力変換回路であり、三相電圧指令値ｖｕ＊，ｖｖ＊，ｖｗ＊に基づいて、電動機２に三相電圧を印加する。なお、インバータは、電圧型インバータまたは電流型インバータを用いることができる。電動機２は、電力変換部７からの三相電流によって駆動し、または回生する。
【００２０】
電流検出器８は、ホール素子などを用いて、電動機２に流れる三相電流ｉｕ，ｉｖ，ｉｗを検出する。
【００２１】
座標変換部１１は、電流検出器８において検出された三相電流を座標変換して、ｄ−ｑ軸上の検出電流ｉｄｓ，ｉｑｓを出力する。検出電流ｉｄｓ，ｉｑｓは、基本波成分と高調波成分を含んでいる。そこで、バンドストップフィルター９とバンドパスフィルター１０を用いて、検出電流ｉｄｓ，ｉｑｓに含まれる基本波成分と高調波成分とを分離する。バンドストップフィルター９は、検出電流ｉｄｓ，ｉｑｓから高調波成分をカットして検出基本波電流ｉｄｓｆ、ｉｑｓｆを出力する。また、バンドパスフィルター１０は、検出電流ｉｄｓ，ｉｑｓから基本波成分をカットして検出高調波電流値ｉｄｓｃ，ｉｑｓｃを出力する。
【００２２】
次に、図３に基づいて、磁石温度推定部１４の構成、機能及び動作について説明する。
磁石温度推定部１４は、フィルタ１５と、バンドストップフィルター１６，１７と、演算器１８と、磁石温度推定器１９とを備えている。
【００２３】
フィルタ１５は、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊のうち、ｑ軸成分をカットしてｄ軸成分であるｄ軸高調波電圧値ｖｄｓｃを出力する。ｄ軸高調波電圧値ｖｄｓｃは、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊を電動機２のｄ軸成分に重畳することにより得られる電圧である。ｄ軸高調波電圧値ｖｄｓｃは、（４）式で表される。
【数４】
2015133890.tif 000006
【００２４】
バンドストップフィルター１６は、電動機２の回転の向きに応じて、ｄ軸高調波電圧値ｖｄｓｃのうち、プラス側またはマイナス側の周波数成分をカットすることで、プラス側、マイナス側いずれかの周波数成分を抽出する。
【００２５】
バンドストップフィルター１７は、バンドストップフィルター１６と同様に、プラス側またはマイナス側の周波数成分をカットすることで、プラス側、マイナス側いずれかの周波数成分を抽出する。
【００２６】
演算器１８は、バンドストップフィルター１６，１７から出力された電圧値及び電流値から高調波インピーダンスＺｄｓを演算する。演算された高調波インピーダンスＺｄｓは、（５）式で表される。
【数５】
2015133890.tif 000007
ここで、ψｚｄは、高調波インピーダンスＺｄｓの相差角である。
【００２７】
続いて、演算器１８は、（５）式から得られる高調波インピーダンスＺｄｓの実部Ｒｄを出力する。
【００２８】
磁石温度推定器１９は、演算器１８から出力された高調波インピーダンスＺｄｓの実部Ｒｄを利用して、永久磁石の温度Ｔｍを推定する。図４に示すように、高調波インピーダンスＺｄｓの実部Ｒｄと永久磁石の温度Ｔｍには、相関関係がある。この関係は、予め実験やシミュレーションを通じて取得することができる。磁石温度推定器１９は、高調波インピーダンスＺｄｓの実部Ｒｄと永久磁石の温度Ｔｍとの関係を示すマップないしは比例係数を記憶しておき、演算器１８から出力された高調波インピーダンスＺｄｓの実部Ｒｄとマップを参照して、永久磁石の温度Ｔｍを出力する。
【００２９】
以上説明してきたように、本実施形態の磁石温度推定システムによれば、基本波のｄ軸成分のみにｄ軸高調波電流指令値ｉｄｓｃ＊をパルセイティング・ベクトル・インジェクション方式で重畳し、電動機２の回転数及びｑ軸インダクタンスに影響を受けないｄ軸高調波電圧値ｖｄｓｃを得る。このｄ軸高調波電流指令値及びｄ軸高調波電圧を用いて、高調波インピーダンスＺｄｓを演算し、演算した高調波インピーダンスＺｄｓの実部Ｒｄに基づいて、永久磁石の温度Ｔｍを推定する。これにより、電動機２の低回転数から高回転数までの広い動作領域において永久磁石の温度Ｔｍを精度よく推定することができる。
【００３０】
また、本実施形態の磁石温度推定システムによれば、ｄ軸高調波電流指令値ｉｄｓｃ＊を所定の間隔で重畳する。ｄ軸高調波電流指令値ｉｄｓｃ＊の周波数や振幅によっては、永久磁石の温度Ｔｍが上昇するおそれがある。そのため、本実施形態では、所定の間隔でｄ軸高調波電流指令値ｉｄｓｃ＊を重畳する。これにより、永久磁石の温度Ｔｍの推定精度を高く保ちつつ、ｄ軸高調波電流指令値ｉｄｓｃ＊による電動機２の損失を抑えることができる。
【００３１】
また、本実施形態の磁石温度推定システムによれば、ｄ軸高調波電流指令値ｉｄｓｃ＊の振幅は、基本波の振幅より小さい。ｄ軸高調波電流指令値ｉｄｓｃ＊の振幅が基本波の振幅より大きい場合、電動機２の動作に影響を与えるだけでなく、永久磁石の温度Ｔｍが上昇するおそれがある。そのため、本実施形態では、ｄ軸高調波電流指令値ｉｄｓｃ＊の振幅を基本波の振幅より小さくする。これにより、永久磁石の温度Ｔｍの推定精度を高く保ちつつ、ｄ軸高調波電流指令値ｉｄｓｃ＊による電動機２の損失を抑えることができる。
【００３２】
また、本実施形態の磁石温度推定システムによれば、ｄ軸高調波電流指令値ｉｄｓｃ＊と、このｄ軸高調波電流指令値ｉｄｓｃ＊を重畳することにより得られたｄ軸高調波電圧値ｖｄｓｃから高調波インピーダンスＺｄｓを求めることができるため、容易に永久磁石の温度Ｔｍを推定することができる。
【００３３】
また、本実施形態の磁石温度推定システムによれば、バンドストップフィルター９、及びバンドパスフィルター１０を用いて、基本波と検出高調波電流値ｉｄｓｃ，ｉｑｓｃを分離する。これにより、高調波成分を高精度で抽出することができ、永久磁石の温度Ｔｍを精度よく推定することができる。
【００３４】
また、本実施形態の磁石温度推定システムによれば、バンドストップフィルター１６，１７を用いて、ｄ軸高調波電流指令値ｉｄｓｃ＊及びｄ軸高調波電圧値ｖｄｓｃから同期電動機の回転方向とは異なる側の周波数成分を抽出する。これにより、基本波の周波数とｄ軸高調波電流指令値ｉｄｓｃ＊の周波数が同じで場合でも、周波数の正負が異なるため、回転数域の広い電動機でも高調波の周波数を高くすることなく、永久磁石の温度Ｔｍを精度よく推定することができる。
【００３５】
次に、第１の実施形態の変形例１について説明する。
図５は、磁石温度推定部１４の変形例１の構成図である。図５に示すように、変形例１では、磁石温度推定部１４は、固定子のコイル温度Ｔｃを考慮して、永久磁石の温度Ｔｍを推定する。固定子のコイル温度Ｔｃは、例えば、図２に示すように固定子コイル３４に取り付けた温度センサ３５から取得することができる。固定子のコイル温度Ｔｃとコイル抵抗値Ｒｃには、図６に示すように、相関関係がある。この関係は、予め実験やシミュレーションを通じて取得することができる。磁石温度推定器１９は、固定子のコイル温度Ｔｃとコイル抵抗値Ｒｃとの関係を示すマップを記憶しておき、入力された固定子のコイル温度Ｔｃとマップを参照して、固定子のコイル抵抗値Ｒｃを推定することができる。固定子のコイル抵抗値Ｒｃは、高調波インピーダンスＺｄｓの実部Ｒｄに含まれる。そのため、固定子のコイル温度Ｔｃから固定子のコイル抵抗値Ｒｃを補正することにより、補正後のコイル抵抗値Ｒｃを含む高調波インピーダンスＺｄｓの実部Ｒｄを精度よく求めることができる。これにより、永久磁石の温度Ｔｍを精度よく推定することができる。
【００３６】
次に、第１の実施形態の変形例２について説明する。
図７は、磁石温度推定部１４の変形例２の構成図である。変形例１と変形例２が異なるのは、変形例２ではマップではなく、演算器２０の演算によって固定子のコイル抵抗値Ｒｃを推定することである。固定子のコイル温度がＴ０からＴ１に変化したとき、コイル温度がＴ０時の抵抗値をＲｃ、コイル温度がＴ１時の抵抗値をＲｃ’とすると、Ｒｃ’は、Ｒｃ’＝Ｒｃ（１＋α×（Ｔ１−Ｔ０））で表される。ここでαは抵抗温度係数である。このように変形例２では、変形例１のようにマップを用いることなく、固定子のコイル温度から固定子のコイル抵抗値Ｒｃを補正することにより、補正後のコイル抵抗値Ｒｃを含む高調波インピーダンスＺｄｓの実部Ｒｄを精度よく求めることができる。これにより、永久磁石の温度Ｔｍを精度よく推定することができる。
【００３７】
次に、第１の実施形態の変形例３について説明する。
図８は、磁石温度推定部１４の変形例３の構成図である。変形例３が、変形例１及び変形例２と異なるのは、図８に示すように、補償演算器２５を備えることである。補償演算器２５は、基本波電流値ｉｄ，ｉｑを用いて、高調波インピーダンスＺｄｓの実部Ｒｄを補償する。これにより、精度の高い高調波インピーダンスＺｄｓの実部Ｒｄ（ｃｏｍｐ）を求めることができる。変形例３では、基本波電流値ｉｄ，ｉｑとして、電流指令値ｉｄｓｆ＊，ｉｑｓｆ＊を用いるが、基本波電流値ｉｄ，ｉｑはこれに限定されるものではない。
次に、図９及び１０を参照して、基本波電流値ｉｄ，ｉｑを用いた高調波インピーダンスＺｄｓの実部Ｒｄを補償する方法を説明する。図９は、ｄ軸基本波電流値ｉｄと高調波インピーダンスＺｄｓの実部Ｒｄを補償するための補償量Ｃｄとの関係を示す。また、図１０は、ｑ軸基本波電流値ｉｑと高調波インピーダンスＺｄｓの実部Ｒｄを補償するための補償量Ｃｑとの関係を示す。図９及び１０に示すように、基本波電流値ｉｄ，ｉｑと高調波インピーダンスＺｄｓの実部Ｒｄを補償するための補償量Ｃｄ，Ｃｑには、相関関係がある。そこで、図９及び１０に示すグラフを線形補完すると、グラフの傾きは一定とみなすことができる。ここで、図９に示すグラフを線形補完した場合の傾きの大きさをαｄ、図１０に示すグラフを線形補完した場合の線分の傾きの大きさをαｑ、補償後の高調波インピーダンスＺｄｓの実部ＲｄをＲｄ（ｃｏｍｐ）とすると、Ｒｄ（ｃｏｍｐ）は、Ｒｄ（ｃｏｍｐ）＝Ｒｄ−Ｉｑ×αｑ＋Ｉｄ×αｄで表される。傾きαｄ，αｑは既知であるため、基本波電流値ｉｄ，ｉｑを入力すれば、補償後の高調波インピーダンスＺｄｓの実部Ｒｄ（ｃｏｍｐ）を求めることができる。
このように変形例３では、基本波電流値ｉｄ，ｉｑを用いることにより、高調波インピーダンスＺｄｓの実部Ｒｄを補償して、精度の高い実部Ｒｄ（ｃｏｍｐ）を求めることができる。そして、精度の高い高調波インピーダンスＺｄｓの実部Ｒｄ（ｃｏｍｐ）を用いることにより、永久磁石の温度Ｔｍを精度よく推定することができる。なお、変形例３において変形例１または２に示す固定子のコイル温度Ｔｃをさらに考慮してもよい。
【００３８】
［第２の実施形態］
次に、本発明の第２の実施形態について説明する。
第２の実施形態が第１実施形態と異なる点は、基本波のｄ軸成分に重畳するのが、高調波電流ではなく高調波電圧である点である。第１の実施形態と共通する点については説明を省略することとし、以下相違点を中心に説明を行う。
【００３９】
図１１は、第２の実施形態に係る永久磁石型同期電動機の制御装置１のシステム構成図である。
フィルタ２１は、バンドパスフィルター１０から出力される検出高調波電流値ｉｄｓｃ，ｉｑｓｃのうち、ｄ軸成分をカットしてｑ軸成分であるｑ軸検出高調波電流値ｉｑｓｃを出力する。
【００４０】
差分器１２は、外部から与えられるｑ軸高調波電流指令値ｉｑｓｃ＊からｑ軸検出高調波電流値ｉｑｓｃを減算する。外部から与えられるｑ軸高調波電流指令値ｉｑｓｃ＊は、（６）式で表される。
【数６】
2015133890.tif 000008
本実施形態では、ｑ軸高調波電流指令値ｉｑｓｃ＊を０とし、基本波のｑ軸成分に高調波電流が重畳しないようにしている。したがって、共振制御部１３で生成されるｑ軸高調波電圧指令値ｖｑｓｃ＊も０となる。加算器２２は、共振制御部１３から出力されるｑ軸高調波電圧指令値ｖｑｓｃ＊と、外部から与えられるｄ軸高調波電圧指令値ｖｄｓｃ＊とを加える。ｄ軸高調波電圧指令値ｖｄｓｃ＊は、（７）式で表される。
【数７】
2015133890.tif 000009
ここで、Ｖｃは振幅である。
【００４１】
加算器５は、第一の電圧指令値ｖｄ０＊，ｖｑ０＊に、高調波電圧指令値ｖｄｓｃ＊，ｖｑｓｃ＊を重畳して、第二の電圧指令値ｖｄｓ＊，ｖｑｓ＊を出力する。ｑ軸高調波電圧指令値ｖｑｓｃ＊は０であるため、加算器５は、基本波のｄ軸成分にのみｄ軸高調波電圧指令値ｖｄｓｃ＊を重畳する。これにより、ｑ軸インダクタンスによる影響を防ぐことができる。また、ｄ軸高調波電圧指令値ｖｄｓｃ＊の重畳方式は、第１の実施形態と同じように、パルセイティング・ベクトル・インジェクション方式である。これにより、電動機２の回転数による影響を抑えた高調波電流を得ることができる。基本波のｄ軸成分にのみｄ軸高調波電圧指令値ｖｄｓｃ＊を重畳することにより得られる高調波電流値ｉｄｓｃ，ｉｑｓｃは、バンドパスフィルター１０を介して出力され、（８）式で表される。
【数８】
2015133890.tif 000010
【００４２】
演算器１８は、図１２に示すように、上記（７）式と（８）式に基づいて、高調波インピーダンスＺｄｓを演算する。高調波インピーダンスＺｄｓは、（９）式で表される。
【数９】
2015133890.tif 000011
【００４３】
続いて、演算器１８は、（９）式から得られる高調波インピーダンスＺｄｓの実部Ｒｄを出力する。
【００４４】
磁石温度推定器１９は、演算器１８から出力された高調波インピーダンスＺｄｓの実部Ｒｄを用いて、永久磁石の温度Ｔｍを推定する。
【００４５】
以上説明してきたように、本実施形態によれば、基本波のｄ軸成分のみにｄ軸高調波電圧指令値ｖｄｓｃ＊をパルセイティング・ベクトル・インジェクション方式で重畳し、電動機２の回転数に影響を抑えたｄ軸高調波電流値ｉｄｓｃを得る。このｄ軸高調波電圧指令値ｖｄｓｃ＊及びｄ軸高調波電流値ｉｄｓｃを用いて、高調波インピーダンスＺｄｓを演算し、演算した高調波インピーダンスＺｄｓの実部Ｒｄに基づいて、永久磁石の温度Ｔｍを推定する。これにより、電動機２の低回転数から高回転数までの広い動作領域において永久磁石の温度Ｔｍを精度よく推定することができる。なお、第１の実施形態と同じように、固定子のコイル温度Ｔｃを用いてコイル抵抗値Ｒｃを補正してもよい。また、第１の実施形態と同じように、基本波電流値ｉｄ，ｉｑを用いて、高調波インピーダンスＺｄｓの実部Ｒｄを補償してもよい。
【００４６】
［第３の実施形態］
次に、本発明の第３の実施形態について説明する。
図１３は、第３の実施形態に係る永久磁石型同期電動機の制御装置１のシステム構成図である。第３の実施形態が第１の実施形態と異なる点は、永久磁石の温度Ｔｍではなく、永久磁石の着磁量Mag.stateを推定することである。第１の実施形態と共通する点については説明を省略することとし、以下相違点を中心に説明を行う。
【００４７】
図１３に示すように、本実施形態に係る制御装置１は、磁石着磁量推定部２３を備える。磁石着磁量推定部２３は、図１４に示すように、演算器１８が演算した高調波インピーダンスＺｄｓの実部Ｒｄから永久磁石の着磁量Mag.stateを推定する。図１５に示すように、永久磁石の着磁量Mag.stateと高調波インピーダンスＺｄｓの実部Ｒｄには、相関関係がある。この関係は、予め実験やシミュレーションを通じて取得することができる。磁石着磁量推定部２３は、永久磁石の着磁量Mag.stateと高調波インピーダンスＺｄｓの実部Ｒｄとの関係を示すマップを記憶しておき、演算器１８から出力された高調波インピーダンスＺｄｓの実部とマップを参照して、永久磁石の着磁量Mag.stateを出力する。本実施形態では、第１の実施形態と同様に、高調波インピーダンスＺｄｓの実部Ｒｄを、電動機２の回転数及びｑ軸インダクタンスの影響を受けずに演算することができる。本実施形態では、この高調波インピーダンスＺｄｓの実部Ｒｄに基づいて永久磁石の着磁量Mag.stateを推定することができる。これにより、電動機２の低回転数から高回転数までの広い動作領域において永久磁石の着磁量Mag.stateを精度よく推定することができる。なお、第１の実施形態と同じように、固定子のコイル温度を用いてコイル抵抗値を補正してもよい。また、第１の実施形態と同じように、基本波電流値ｉｄ，ｉｑを用いて、高調波インピーダンスＺｄｓの実部Ｒｄを補償してもよい。
【００４８】
欧州特許出願１４３８０００２．７（出願日２０１４年１月１３日）の全内容がここに援用される。
【００４９】
以上、実施例に沿って本発明の内容を説明したが、本発明はこれらの記載に限定されるものではなく、種々の変形及び改良が可能であることは、当業者には自明である。
【符号の説明】
【００５０】
１  制御装置
２  電動機
３、１２  差分器
４  電流制御部
５  加算器（重畳部）
６、１１  座標変換部
７  電力変換部
８  電流検出器
９、１６、１７  バンドストップフィルター
１０  バンドパスフィルター
１３  共振制御部
１４  磁石温度推定部
１５、２１  フィルタ
１８  演算器（演算部）
１９  磁石温度推定器
２０  演算器
２２  加算器
２３  磁石着磁量推定部
２４  磁石着磁量推定器
２５  補償演算器（補償演算部）
３０  固定子
３１  回転子
３２  永久磁石
３３  スロット
３４  固定子コイル
３５  温度センサ（温度計測部）
【図１】
2015133890.tif 000018
【図２】
2015133890.tif 000019
【図３】
2015133890.tif 000020
【図４】
2015133890.tif 000021
【図５】
2015133890.tif 000022
【図６】
2015133890.tif 000023
【図７】
2015133890.tif 000024
【図８】
2015133890.tif 000025
【図９】
2015133890.tif 000026
【図１０】
2015133890.tif 000027
【図１１】
2015133890.tif 000028
【図１２】
2015133890.tif 000029
【図１３】
2015133890.tif 000030
【図１４】
2015133890.tif 000031
【図１５】
2015133890.tif 000032