(57)【要約】モータ制御方法は、αβ固定座標系またはｄｑ回転座標系を基準とした、フェーザ表示による、合成磁束およびステータ電流を獲得するステップと、式（１）に基づいてトルク角δを演算するステップであって、【数１】2018159102.tif 000012 ここで、Ｌは、電機子インダクタンスであり、Ψｍは、ステータ磁石の磁束を示し、Ψｓは合成磁束の大きさを示し、ＩＳはステータ電流の大きさを示す、ステップと、トルク角δに基づいてモータを制御するステップと、を包含する。
2018159102.tif 000002
【特許請求の範囲】
【請求項１】
表面磁石型モータを制御するモータ制御方法であって、
αβ固定座標系またはｄｑ回転座標系を基準とした、フェーザ表示による、合成磁束およびステータ電流を獲得するステップと、
式（１）に基づいてトルク角δを演算するステップであって、
【数１】
2018159102.tif 000010
ここで、Ｌは、電機子インダクタンスであり、Ψｍは、ステータ磁石の磁束を示し、Ψｓは前記合成磁束の大きさを示し、ＩＳは前記ステータ電流の大きさを示す、ステップと、
前記トルク角δに基づいてモータを制御するステップと、
を包含するモータ制御方法。
【請求項２】
前記トルク角δに基づいてトルクＴを演算するステップをさらに包含し、
前記モータを制御するステップにおいて、前記トルクＴに基づいて前記表面磁石型モータを制御する、請求項１に記載のモータ制御方法。
【請求項３】
前記αβ固定座標系における前記合成磁束のα軸およびβ軸上の成分に基づいて位相角ρを演算し、かつ、前記トルク角δおよび前記位相角ρに基づいてモータのロータ角θを演算するステップをさらに包含し、
前記モータを制御するステップにおいて、前記ロータ角θおよび前記トルクＴに基づいて前記表面磁石型モータを制御する、請求項２に記載のモータ制御方法。
【請求項４】
表面磁石型モータと、
前記表面磁石型モータを制御する制御回路と、
を有し、
前記制御回路は、
αβ固定座標系またはｄｑ回転座標系を基準とした、フェーザ表示による、合成磁束およびステータ電流を獲得し、
式（２）に基づいてトルク角δを演算し、
【数２】
2018159102.tif 000011
ここで、Ｌは、電機子インダクタンスであり、Ψｍは、ステータ磁石の磁束を示し、Ψｓは前記合成磁束の大きさを示し、ＩＳは前記ステータ電流の大きさを示し、
前記トルク角δに基づいてモータを制御する、モータ制御システム。
【請求項５】
請求項４に記載のモータ制御システムを有する電動パワーステアリングシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、モータ制御方法、モータ制御システムおよび電動パワーステアリングシステムに関する。
【背景技術】
【０００２】
近年、電気駆動システムが様々な応用分野に広く用いられる。電気駆動システムとして、例えばモータ制御システムが挙げられる。モータ制御システムは、例えばベクトル制御を用いて電動モータ（以下、「モータ」と表記する。）を制御する。ベクトル制御には、例えば、電流センサおよび位置センサを用いる方式（以下、「センサ制御」と称する。）と、電流センサのみを用いる方式（以下、「センサレス制御」と称する。）と、がある。センサ制御では、位置センサの測定値に基づいてロータの位置（以下、「ロータ角」と称する。）が算出される。一方、センサレス制御では、ロータ角は、電流センサによって測定される電流などに基づいて推定される。  【０００３】
ベクトル制御には、一般にトルク情報が必要とされる。トルクは、例えばモータのトルク角に基づいて演算することが可能である。特に、センサレス制御では、トルク角に基づいてロータ角を推定することが求められる。このように、ベクトル制御の精度向上には、トルク角を正確に取得することが不可欠とされる。例えば、センサ制御において、トルク角は、ｄｑ回転座標系における変数を用いて演算できることが知られている。トルク角は負荷角とも称される。  【０００４】
特許文献１は、いわゆるオブザーバを用いてトルク角を推定するセンサレス制御を開示する。具体的に説明すると、オブザーバは、電流センサで測定された電流値に基づいてロータ角を推定し、さらに、推定されたロータ角に基づいてフィードバック・トルク角を推定する。特許文献２は、トルクの推定値に基づいてトルク角を求める演算式を開示する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２００７−５２５１３７号公報
【特許文献２】中国特許出願公開第１０３６８４１６９号明細書
【非特許文献】
【０００６】
【非特許文献１】Ghaderi, Ahmad, and Tsuyoshi Hanamoto. "Wide-speed-range sensorless vector control of synchronous reluctance motors based on extended programmable cascaded low-pass filters." IEEE Transactions on Industrial Electronics, Vol. 58, No. 6, (June 2011), p.2322-2333.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
センサ制御に利用される、ｄｑ回転座標系における変数に基づくトルク角の演算は、センサレス制御には適用できない場合がある。その理由は以下のとおりである。ｄｑ回転座標系は、ロータと共に回転する回転座標系であり、ロータ角および回転速度に基づいて設定される座標系である。一方で、センサレス制御では、ロータ角の推定にトルク角が必要とされることがある。その場合、センサレス制御において、ｄｑ回転座標系における変数に依存しない、トルク角を求める手法が求められる。  【０００８】
センサレス制御において、特許文献１に開示されているようなオブザーバを用いるトルク角の推定は、通常、モータに関する様々のパラメータ（例えば、電機子インダクタンスおよびリアクタンス）を必要とし、かつ、それらに強く影響を受ける。例えば非特許文献１で言及されているように、オブザーバを用いる推定は、特に初期値およびノイズ共分散行列に強く依存するとされている。その結果、それらの値および行列を誤って選択すると、モータ制御を不安定にさせる可能性がある。さらに、オブザーバによる推定には、より複雑な演算が必要とされる。そのため、コンピュータに対する演算負荷が増大するといった課題が生じる。以上の理由により、センサレス制御において、複雑な演算を特に必要としない、トルク角を推定するための手法が望まれる。  【０００９】
本開示の実施形態は、センサレス制御において、ｄｑ回転座標系における変数に依存せずにトルク角を推定することが可能な、新規なモータ制御方法、モータ制御システム、および、当該モータ制御システムを有する電動パワーステアリングシステムを提供する。
【課題を解決するための手段】
【００１０】
本開示の例示的なモータ制御方法は、表面磁石型モータを制御する方法であって、αβ固定座標系またはｄｑ回転座標系を基準とした、フェーザ表示による、合成磁束およびステータ電流を獲得するステップと、式（１）に基づいてトルク角δを演算するステップであって、    
【数１】
2018159102.tif 000003
ここで、Ｌは、電機子インダクタンスであり、Ψｍは、ステータ磁石の磁束を示し、Ψｓは前記合成磁束の大きさを示し、ＩＳは前記ステータ電流の大きさを示す、ステップと、前記トルク角δに基づいてモータを制御するステップと、を包含する。  【００１１】
本開示の例示的なモータ制御システムは、表面磁石型モータと、前記表面磁石型モータを制御する制御回路と、を有し、前記制御回路は、αβ固定座標系またはｄｑ回転座標系を基準とした、フェーザ表示による、合成磁束およびステータ電流を獲得し、式（２）に基づいてトルク角δを演算し、    
【数２】
2018159102.tif 000004
ここで、Ｌは、電機子インダクタンスであり、Ψｍは、ステータ磁石の磁束を示し、Ψｓは前記合成磁束の大きさを示し、ＩＳは前記ステータ電流の大きさを示し、前記トルク角δに基づいてモータを制御する。
【発明の効果】
【００１２】
本開示の例示的な実施形態によると、センサレス制御において、ｄｑ回転座標系における変数に依存せずにトルク角を求めることが可能な、新規なモータ制御方法、モータ制御システム、および当該モータ制御システムを有する電動パワーステアリングシステムが提供される。
【図面の簡単な説明】
【００１３】
【図１】図１は、実施形態１によるモータ制御システム１０００のハードウェアブロックを示すブロック図である。
【図２】図２は、実施形態１によるモータ制御システム１０００中のインバータ３００のハードウェア構成を示すブロック図である。
【図３】図３は、実施形態１の変形例によるモータ制御システム１０００のハードウェアブロックを示すブロック図である。
【図４】図４は、コントローラ１００の機能ブロック示す機能ブロック図である。
【図５】図５は、変数Ｉｓ、Ψｓ、ΦおよびＶｓを表示するフェーザ図である。
【図６】図６は、αβ固定座標系またはｄｑ回転座標系における合成磁束Ψｓを表示するフェーザ図である。
【図７】図７は、ロータ磁束Ψｍ、電機子磁束Ψａおよび合成磁束Ψｓを表すフェーザ図である。
【図８】図８は、所定期間内の、トルクの波形（上）、三相電流の波形（中間）、および、三相電圧の波形（下）を示すグラフである。
【図９】図９は、本開示の演算式を用いて推定された所定期間内のトルク角（度）、および、トルク角の実測値の波形を示すグラフである。
【図１０】図１０は、実施形態２によるＥＰＳシステム２０００の典型的な構成を示す模式図である。
【発明を実施するための形態】
【００１４】
以下、添付の図面を参照しながら、本開示のモータ制御方法、モータ制御システム、および当該モータ制御システムを有する電動パワーステアリングシステムの実施形態を詳細に説明する。但し、以下の説明が不必要に冗長になるのを避け、当業者の理解を容易にするため、必要以上に詳細な説明は省略する場合がある。例えば、既によく知られた事項の詳細説明や実質的に同一の構成に対する重複説明を省略する場合がある。  
【００１５】
（実施形態１）
〔モータ制御システム１０００の構成〕
図１は、本実施形態によるモータ制御システム１０００のハードウェアブロックを模式的に示す。
【００１６】
モータ制御システム１０００は、典型的に、モータＭと、コントローラ（制御回路）１００と、駆動回路２００と、インバータ（「インバータ回路」とも称される。）３００と、複数の電流センサ４００と、アナログデジタル変換回路（以下、「ＡＤコンバータ」と表記する。）５００と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）６００とを有する。モータ制御システム１０００は、モジュール化され、例えば、モータ、センサ、ドライバおよびコントローラを有するモータモジュールとして製造および販売され得る。本明細書では、構成要素としてモータＭを有するシステムを例に、モータ制御システム１０００を説明する。ただし、モータ制御システム１０００は、構成要素としてモータＭを有しない、モータＭを駆動するためのシステムであってもよい。  
【００１７】
モータＭは、表面磁石型（ＳＰＭ）モータであり、例えば表面磁石型同期モータ（ＳＰＭＳＭ）である。モータＭは、例えば三相（Ｕ相、Ｖ相およびＷ相）の巻線（不図示）を有する。三相の巻線は、インバータ３００に電気的に接続される。三相モータに限らず、五相、七相などの多相モータも本開示の範疇である。本明細書では、三相モータを制御するモータ制御システムを例に、本開示の実施形態を説明する。  
【００１８】
コントローラ１００は、例えばマイクロコントロールユニット（ＭＣＵ）である。または、コントローラ１００は、例えば、ＣＰＵコアが組み込まれたフィールドプログラマブルゲートアレイ（ＦＰＧＡ）によっても実現し得る。  
【００１９】
コントローラ１００は、モータ制御システム１０００の全体を制御し、例えばベクトル制御によってモータＭのトルクおよび回転速度を制御する。モータＭは、ベクトル制御に限らず、他のクローズドループ制御によっても制御され得る。回転速度は、単位時間（例えば１分間）にロータが回転する回転数（ｒｐｍ）または単位時間（例えば１秒間）にロータが回転する回転数（ｒｐｓ）で表される。ベクトル制御は、モータに流れる電流を、トルクの発生に寄与する電流成分と、磁束の発生に寄与する電流成分とに分解し、互いに直交する各電流成分を独立に制御する方法である。コントローラ１００は、例えば、複数の電流センサ４００によって測定された実電流値、および実電流値に基づいて推定されたロータ角などに従って目標電流値を設定する。コントローラ１００は、その目標電流値に基づいてＰＷＭ（Ｐｕｌｓｅ  Ｗｉｄｔｈ  Ｍｏｄｕｌａｔｉｏｎ）信号を生成し、駆動回路２００に出力する。  
【００２０】
駆動回路２００は、例えばゲートドライバである。駆動回路２００は、インバータ３００におけるスイッチング素子のスイッチング動作を制御する制御信号を、コントローラ１００から出力されるＰＷＭ信号に従って生成する。後述するように、駆動回路２００は、コントローラ１００に実装されていてもよい。  
【００２１】
インバータ３００は、例えば直流電源（不図示）から供給される直流電力を交流電力に変換し、変換された交流電力でモータＭを駆動する。例えば、インバータ３００は、駆動回路２００から出力される制御信号に基づいて、直流電力を、Ｕ相、Ｖ相およびＷ相の擬似正弦波である三相交流電力に変換する。この変換された三相交流電力でモータＭは駆動される。  
【００２２】
複数の電流センサ４００は、モータＭのＵ相、Ｖ相およびＷ相の巻線に流れる少なくとも２つの電流を検出する少なくとも２つの電流センサを有する。本実施形態では、複数の電流センサ４００は、Ｕ相およびＶ相に流れる電流を検出する２つの電流センサ４００Ａ、４００Ｂ（図２を参照）を有する。当然に、複数の電流センサ４００は、Ｕ相、Ｖ相およびＷ相の巻線に流れる３つの電流を検出する３つの電流センサを有していてもよいし、例えばＶ相およびＷ相に流れる電流またはＷ相およびＵ相に流れる電流を検出する２つの電流センサを有していてもよい。電流センサは、例えば、シャント抵抗、およびシャント抵抗に流れる電流を検出する電流検出回路（不図示）を有する。シャント抵抗の抵抗値は、例えば０．１Ω程度である。  
【００２３】
ＡＤコンバータ５００は、複数の電流センサ４００から出力されるアナログ信号をサンプリングしてデジタル信号に変換し、この変換したデジタル信号をコントローラ１００に出力する。コントローラ１００がＡＤ変換を行ってもよい。その場合、複数の電流センサ４００は、アナログ信号をコントローラ１００に直接出力する。  
【００２４】
ＲＯＭ６００は、例えば書き込み可能なメモリ（例えばＰＲＯＭ）、書き換え可能なメモリ（例えばフラッシュメモリ）または読み出し専用のメモリである。ＲＯＭ６００は、コントローラ１００にモータＭを制御させるための命令群を有する制御プログラムを格納する。例えば、制御プログラムはブート時にＲＡＭ（不図示）に一旦展開される。ＲＯＭ６００は、コントローラ１００に外付けされる必要はなく、コントローラ１００に搭載されていてもよい。ＲＯＭ６００を搭載したコントローラ１００は、例えば上述したＭＣＵであり得る。  
【００２５】
図２を参照して、インバータ３００のハードウェア構成を詳細に説明する。
【００２６】
図２は、本実施形態によるモータ制御システム１０００中のインバータ３００のハードウェア構成を模式的に示す。  
【００２７】
インバータ３００は、３個のローサイドスイッチング素子および３個のハイサイドスイッチング素子を有する。図示されるスイッチング素子ＳＷ＿Ｌ１、ＳＷ＿Ｌ２およびＳＷ＿Ｌ３がローサイドスイッチング素子であり、スイッチング素子ＳＷ＿Ｈ１、ＳＷ＿Ｈ２およびＳＷ＿Ｈ３が、ハイサイドスイッチング素子である。スイッチング素子として、例えば、電界効果トランジスタ（ＦＥＴ、典型的にはＭＯＳＦＥＴ）または絶縁ゲートバイポーラトランジスタ（ＩＧＢＴ）などの半導体スイッチ素子を用いることができる。スイッチング素子は、モータＭに向けて流れる回生電流を流す還流ダイオードを有する。  
【００２８】
図２に、Ｕ相およびＶ相に流れる電流を検出する２つの電流センサ４００Ａ、４００Ｂのシャント抵抗Ｒｓを示す。図示されるように、例えばシャント抵抗Ｒｓは、ローサイドスイッチング素子とグランドとの間に電気的に接続され得る。または、例えばシャント抵抗Ｒｓは、ハイサイドスイッチング素子と電源との間に電気的に接続され得る。  
【００２９】
コントローラ１００は、例えばベクトル制御に基づく三相通電による制御（以下、「三相通電制御」と表記する。）を行うことによってモータＭを駆動することができる。例えば、コントローラ１００は、三相通電制御を行うためのＰＷＭ信号を生成し、そのＰＷＭ信号を駆動回路２００に出力する。駆動回路２００は、インバータ３００中の各ＦＥＴのスイッチング動作を制御するゲート制御信号をＰＷＭ信号に基づいて生成し、各ＦＥＴのゲートに与える。  
【００３０】
図３は、本実施形態の変形例によるモータ制御システム１０００のハードウェアブロックを模式的に示す。  
【００３１】
図示されるように、モータ制御システム１０００は、駆動回路２００を有していなくてもよい。その場合、コントローラ１００は、インバータ３００の各ＦＥＴのスイッチング動作を直接制御することが可能なポートを有する。具体的に説明すると、コントローラ１００は、ゲート制御信号をＰＷＭ信号に基づいて生成することが可能である。コントローラ１００は、そのポートを介してゲート制御信号を出力し、そのゲート制御信号を各ＦＥＴのゲートに与えることができる。  
【００３２】
図３に示されるように、モータ制御システム１０００は、位置センサ７００をさらに有していてもよい。位置センサ７００は、モータＭに配置され、ロータ角を検出してコントローラ１００に出力する。位置センサ７００は、例えば磁気抵抗（ＭＲ）素子を有するＭＲセンサとセンサマグネットとの組み合わせによって実現される。位置センサ７００は、例えば、ホール素子を含むホールＩＣまたはレゾルバを用いても実現される。  
【００３３】
モータ制御システム１０００は、位置センサ７００の代わりに、例えば、速度センサまたは加速度センサを有し得る。コントローラ１００は、位置センサとして速度センサを用いる場合、回転速度信号または角速度信号に積分処理等を行うことによりロータ角、つまり、回転角を算出することができる。角速度は、１秒間にロータが回転する角度（ｒａｄ／ｓ）で表される。また、コントローラ１００は、位置センサとして加速度センサを用いる場合、角加速度信号に積分処理等を行うことにより回転角を算出することができる。
【００３４】
本開示のモータ制御システムは、例えば図１および２に示されるような、位置センサを有しない、センサレス制御を行うためのモータ制御システムに利用され得る。また、本開示のモータ制御システムは、例えば図３に示されるような、位置センサを有する、センサ制御を行うためのモータ制御システムにも利用され得る。  
【００３５】
以下、図４から図７を参照しながら、センサレス制御用のモータ制御システムを例に、そのシステムに用いられるモータ制御方法の具体例を説明し、トルク角の推定に用いる演算を主に説明する。本開示のモータ制御方法は、トルク角の推定が要求される、ＳＰＭモータを制御するための様々なモータ制御システムに利用され得る。  
【００３６】
〔モータ制御システム１０００の制御方法〕
モータ制御システム１０００の制御方法の概要は以下のとおりである。
【００３７】
まず、電流センサ４００で測定された三相電流Ｉａ、ＩｂおよびＩｃをαβ固定座標系におけるα軸およびβ軸上の電流Ｉα、Ｉβに変換する。次に、電流Ｉα、Ｉβに基づいて、位相角ρを演算し、かつ、ステータ電流Ｉｓおよび合成磁束Ψｓを演算する。次に、ステータ電流Ｉｓおよび合成磁束Ψｓに基づいてトルク角δを推定し、かつ、モータ制御に必要なトルクＴおよびロータ角θをトルク角δに基づいて決定する。最終的に、トルクＴおよびロータ角θに基づいてモータＭを制御する。  
【００３８】
本実施形態によるモータ制御方法を実現するためのアルゴリズムは、例えば特定用途向け集積回路（ＡＳＩＣ）またはＦＰＧＡなどのハードウェアのみで実現することもできるし、ハードおよびソフトウェアの組み合わせによっても実現することができる。  
【００３９】
図４は、トルク角δを推定するための、コントローラ１００の機能ブロックを模式的に示す。本明細書において、機能ブロック図における各ブロックは、ハードウェア単位ではなく機能ブロック単位で示される。モータ制御用ソフトウェアは、例えば、各機能ブロックに対応した特定の処理を実行させるためのコンピュータプログラムを構成するモジュールであり得る。そのようなコンピュータプログラムは、例えばＲＯＭ６００に格納される。  
【００４０】
図４に示されるように、コントローラ１００は、例えば、プレ演算ユニット１１０、トルク角演算ユニット１２０、位相角演算ユニット１３０、ロータ角演算ユニット１４０、トルク演算ユニット１５０およびモータ制御ユニット１６０を有する。コントローラ１００は、ステータ電流Ｉｓおよび合成磁束Ψｓに基づいてトルク角δを演算することができる。本明細書において、説明の便宜上、各機能ブロックをユニットと表記することとする。当然に、この表記は、各機能ブロックを、ハードウェアまたはソフトウェアに限定解釈する意図で用いられない。  
【００４１】
各機能ブロックがソフトウェアとしてコントローラ１００に実装される場合、そのソフトウェアの実行主体は、例えばコントローラ１００のコアであり得る。上述したように、コントローラ１００は、ＦＰＧＡによって実現され得る。その場合、全てまたは一部の機能ブロックは、ハードウェアで実現され得る。  
【００４２】
複数のＦＰＧＡを用いて処理を分散させることにより、特定のコンピュータの演算負荷を分散させることができる。その場合、図４に示される機能ブロックの全てまたは一部は、その複数のＦＰＧＡに分散して実装され得る。複数のＦＰＧＡは、例えば車載のコントロールエリアネットワーク（ＣＡＮ）によって互いに通信可能に接続され、データの送受信がなされる。  
【００４３】
例えば三相通電制御において、各相を流れる電流の総和は理想的にゼロになる。本明細書において、モータＭのＵ相の巻線に流れる電流をＩａ、モータＭのＶ相の巻線に流れる電流をＩｂ、および、モータＭのＷ相の巻線に流れる電流をＩｃとする。電流Ｉａ、ＩｂおよびＩｃの総和はゼロになる。  
【００４４】
コントローラ１００（例えばＣＰＵコア）は、電流Ｉａ、ＩｂおよびＩｃのうちの２つの電流を受け取って残りの１つの電流を演算により求める。本実施形態では、コントローラ１００は、電流センサ４００Ａで測定された電流Ｉａおよび電流センサ４００Ｂで測定された電流Ｉｂを取得する。コントローラ１００は、電流Ｉａ、ＩｂおよびＩｃの総和はゼロになる上記関係を用いて、電流Ｉａ、Ｉｂに基づいて電流Ｉｃを演算する。３つの電流センサを用いて電流Ｉａ、ＩｂおよびＩｃを測定し、それらをＡＤコンバータ５００を介してコントローラ１００に入力する構成を採用しても構わない。  
【００４５】
コントローラ１００は、ベクトル制御などに用いられるいわゆるクラーク変換を用いて、電流Ｉａ、ＩｂおよびＩｃを、αβ固定座標系における、α軸上の電流Ｉαおよびβ軸上の電流Ｉβに変換することができる。ここで、αβ固定座標系は静止座標系である。三相のうちの一相の方向（例えばＵ相方向）がα軸であり、α軸と直交する方向がβ軸である。  
【００４６】
コントローラ１００はさらに、クラーク変換を用いて、リファレンス電圧Ｖａ＊、Ｖｂ＊およびＶｃ＊を、αβ固定座標系における、α軸上のリファレンス電圧Ｖα＊およびβ軸上のリファレンス電圧Ｖβ＊に変換する。リファレンス電圧Ｖａ＊、Ｖｂ＊およびＶｃ＊は、インバータ３００の各スイッチング素子を制御するための、上述したＰＷＭ信号を表す。  
【００４７】
例えば、電流Ｉα、Ｉβ、リファレンス電圧Ｖα＊およびＶβ＊を求める演算は、コントローラ１００のモータ制御ユニット１６０によっても実行され得る。電流Ｉα、Ｉβ、リファレンス電圧Ｖα＊およびＶβ＊は、プレ演算ユニット１１０および位相角演算ユニット１３０に入力される。  
【００４８】
本実施形態によるモータ制御において、ステータ電流Ｉｓ、合成磁束Ψｓおよび位相角Φは、変数として与えられ、電機子抵抗Ｒ（ｍΩ）、電機子インダクタンスＬ（μＨ）およびロータ磁束Ψｍ（Ｗｂ）は、パラメータとして与えられる。ここで、ロータ磁束Ψｍは、ロータの永久磁石の磁束の大きさを示す。  
【００４９】
プレ演算ユニット１１０は、電流Ｉα、Ｉβ、リファレンス電圧Ｖα＊およびＶβ＊に基づいて、αβ固定座標系またはｄｑ回転座標系を基準とした、変数ＩｓおよびΨｓを獲得する。ｄｑ回転座標系は、ロータと共に回転する回転座標系である。プレ演算ユニット１１０は、後段のトルク角演算ユニット１２０に上記の変数を渡すためにプレ演算を行うためのユニットである。  
【００５０】
図５は、変数Ｉｓ、Ψｓ、ΦおよびＶｓを表示するフェーザ図である。図６は、αβ固定座標系またはｄｑ回転座標系における合成磁束Ψｓを表示するフェーザ図である。図示される変数はすべてフェーザ表示によって表される。以下、各変数をフェーザとして扱う。  
【００５１】
＜変数：ステータ電流Ｉｓ＞
プレ演算ユニット１１０は、式（１）に基づいてフェーザ図におけるステータ電流Ｉｓを演算する。
Ｉｓ＝（Ｉα２＋Ｉβ２）１／２      式（１）
【００５２】
＜変数：合成磁束Ψｓ＞
プレ演算ユニット１１０は、電流Ｉα、Ｉβ、リファレンス電圧Ｖα＊およびＶβ＊に基づいてフェーザ図における合成磁束Ψｓを演算する。具体的に説明すると、プレ演算ユニット１１０は、式（２）から（４）に基づいて合成磁束Ψｓを演算する。合成磁束Ψｓは、図５に示されるように、ロータ磁束Ψｍに電機子磁束Ψａ（＝Ｌ・Ｉｓ）を加算することにより得られる。  
【００５３】
プレ演算ユニット１１０は、例えば、式（２）に基づいて合成磁束Ψｓのα軸上の成分Ψαを演算する。プレ演算ユニット１１０は、式（３）に基づいて合成磁束Ψｓのβ軸上の成分Ψβを演算する。ここで、式（２）および（３）の中のＬＰＦは、ローパスフィルタによる処理を意味する。高調波を除去する目的で、例えばコントローラ１００が有する汎用ローパスフィルタを用いることができる。合成磁束Ψｓは、式（４）で表される。
Ψα＝ＬＰＦ（Ｖα＊−Ｒ・Ｉα）              式（２）
Ψβ＝ＬＰＦ（Ｖβ＊−Ｒ・Ｉβ）              式（３）
Ψｓ＝（Ψα２＋Ψβ２）１／２                  式（４）
【００５４】
プレ演算ユニット１１０は、変数Ｉｓ、Ψｓをトルク角演算ユニット１２０に出力する。コントローラ１００とは異なる他のハードウェア（例えば、ＦＰＧＡ）が変数Ｉｓ、Ψｓを演算してもよい。トルク角演算ユニット１２０は、他のハードウェアから変数Ｉｓ、Ψｓを受け取ることにより、それらを獲得してもよい。このような構成によると、コントローラ１００の演算負荷を低減することが可能となる。  
【００５５】
トルク角演算ユニット１２０は、パラメータＬ、Ψｍ、変数ＩｓおよびΨｓに基づいてトルク角δを演算する。トルク角δは、図６において、例えばｄｑ回転座標系における合成磁束Ψｓとｄ軸との間の角度で表され、反時計方向を正の方向とする角度である。  
【００５６】
図７は、ロータ磁束Ψｍ、電機子磁束Ψａおよび合成磁束Ψｓを表すフェーザ図である。  
【００５７】
図示されるように、Ψｍ、ΨａおよびΨｓを三辺とする三角形に、いわゆる余弦定理を適用すると、式（５）が得られる。式（５）を変形すると、式（６）が得られ、ｃｏｓ（δ）の逆三角関数を求めると、δに関する式（７）が得られる。Ψａ＝Ｌ・Ｉｓを式（７）に代入することにより、δに関する式（８）が最終的に得られる。
Ψａ２＝Ψｓ２＋Ψｍ２−２ΨｓΨｍｃｏｓ（δ）      式（５）
【数３】
2018159102.tif 000005
【数４】
2018159102.tif 000006
【数５】
2018159102.tif 000007
【００５８】
トルク角演算ユニット１２０は、トルク角δをトルク演算ユニット１５０およびロータ角演算ユニット１４０に出力する。式（８）に示されるように、トルク角δの推定に、ｄｑ回転座標系における変数および位相角Φは必要とされない。位相角Φは、図５に示されるように、例えばｄｑ回転座標系において、ステータ電流Ｉｓとステータ電圧Ｖｓとの間の角度で表され、反時計方向を正の方向とする角度である。ステータ電圧Ｖｓは、逆起電力電圧に対応した電圧である。本明細書では、逆起電力電圧をステータ電圧と呼ぶ。本実施形態によれば、パラメータＬ、Ψｍ、変数ＩｓおよびΨｓに基づいてトルク角δを演算することが可能となる。  
【００５９】
位相角演算ユニット１３０は、電流Ｉα、Ｉβ、リファレンス電圧Ｖα＊およびＶβ＊に基づいて位相角ρを推定する。位相角演算ユニット１３０は、プレ演算ユニットと同様に、例えば上記の式（２）および（３）に基づいて磁束成分Ψα、Ψβを演算する。位相角演算ユニット１３０は、さらに、例えば式（９）に基づいて位相角ρを演算する。位相角ρは、例えば図６に示されるように、αβ固定座標系において、合成磁束Ψｓとα軸との間の角度で表され、反時計方向を正の方向とする角度である。位相角演算ユニット１３０は、位相角ρをロータ角演算ユニット１４０に出力する。
ρ＝ｔａｎ−１（Ψβ／Ψα）                  式（９）
【００６０】
ロータ角演算ユニット１４０は、トルク角δおよび位相角ρに基づいてロータ角θを演算する。トルク角δ、位相角ρおよびロータ角θの関係は、図６に示されるとおりである。ロータ角演算ユニット１４０は、式（１０）に基づいてロータ角θを演算し、推定することができる。
θ＝ρ−δ      式（１０）
【００６１】
トルク演算ユニット１５０は、トルク角δに基づいてトルクＴを演算する。ＳＰＭモータを用いる場合、突極比（Ｌｄ／Ｌｑ）は１（つまり、Ｌ＝Ｌｄ＝Ｌｑ）となる。その場合、電機子に働くトルクの反作用として、トルクＴは式（１１）によって表されることが知られている。トルク演算ユニット１５０は、例えば式（１１）に基づいてトルクＴを演算することができる。    
【数６】
2018159102.tif 000008
ここでＰはモータ極対数を示すパラメータである。  
【００６２】
モータ制御ユニット１６０は、トルクＴおよびロータ角θに基づいてモータＭを制御することができる。モータ制御ユニット１６０は、例えば一般的なベクトル制御に必要な演算を行う。ベクトル制御は周知の技術であるので、その制御についての詳細な説明は省略する。  
【００６３】
本実施形態によると、センサレス制御において、ｄｑ回転座標系における変数に依存せずにトルク角を求めることが可能となる。また、トルク角の推定に複雑な演算は特に必要されないので、コンピュータに対する負荷を低減することが可能となり、かつ、メモリコストを低減することが可能となる。  
【００６４】
以下に、本開示によるトルク角δの演算の妥当性を、ｄＳＰＡＣＥ社の”ラピッドコントロールプロトタイピング（ＲＣＰ）システム”およびＭａｔｈＷｏｒｋｓ社のＭａｔｌａｂ／Ｓｉｍｕｌｉｎｋを用いて検証した結果を示す。この検証には、ベクトル制御により制御を受けるＳＰＭモータのモデルが用いられた。表１には、検証時の各種システムパラメータの値が示される。  
【００６５】
【表１】
2018159102.tif 000009
【００６６】
図８は、所定期間内（０．３５秒から０．３８秒までの０．０３秒）の、トルクの波形（上）、三相電流の波形（中間）、および、三相電圧の波形（下）を示す。図９は、本開示の演算式を用いて推定された所定期間内のトルク角（度）、および、トルク角の実測値の波形を示す。図８および図９の横軸は時間（ｍｓ）を表す。図８の縦軸は、上側から順番に、トルクの大きさ（Ｎ・ｍ）、電流値（ｍＡ）および電圧値（Ｖ）を表す。図９の縦軸は、トルク角の大きさ（度）を表す。  
【００６７】
図８のシミュレーション結果から、ベクトル制御が適切になされていることが分かる。また、図９のシミュレーション結果から、本開示の演算式を用いて推定されたトルク角δ、および、実測値は類似することが分かる。より詳細には、推定されたトルク角δと実測値との誤差は約１度である。センサレス制御において、一般に、その誤差の許容値は１０度程度とされている。本シミュレーション結果から得られた誤差は、その許容値の範囲に十分に収まる値である。  
【００６８】
以上のシミュレーション結果から、本明細書に提案する、トルク角を演算するための手法を用いることにより、センサレス制御においてトルク角を精度よく推定できることが分かる。  
【００６９】
本開示によるトルク角δの推定手法は、上述したとおり、センサレス制御に限らず、図３に示されるセンサ制御用のモータ制御システムにも好適に利用され得る。  
【００７０】
図３に示されるモータ制御システム１０００中のコントローラ１００は、ｄｑ回転座標系における変数に基づいてトルク角δを演算することができる。コントローラ１００は、例えば式（１２）に基づいてトルク角δを演算することが可能である（図５を参照）。
δ＝ｔａｎ−１〔（Ｖｄ−Ｒ・Ｉｄ）／（Ｖｑ−Ｒ・Ｉｑ）〕    式（１２）
ここで、Ｖｄは電機子電圧のｄ軸上の電圧成分であり、Ｖｑは電機子電圧のｑ軸上の電圧成分である。Ｉｄは電機子電流のｄ軸上の電流成分であり、Ｉｑは電機子電流のｑ軸上の電流成分である。
【００７１】
センサ制御において、位置センサが何らかの原因で破損した場合、ロータ角を測定することはできなくなる。そのため、センサ制御を継続することは困難となる。一方で、位置センサが故障した場合、モータ制御を、センサ制御からセンサレス制御に切替えることが可能である。そのセンサレス制御に、本開示によるトルク角の推定手法を適用することにより、位置センサが故障した場合でも、モータ制御を継続することが可能となる。  
【００７２】
（実施形態２）
図１０は、本実施形態によるＥＰＳシステム２０００の典型的な構成を模式的に示す。  
【００７３】
自動車等の車両は一般に、ＥＰＳシステムを有する。本実施形態によるＥＰＳシステム２０００は、ステアリングシステム５２０、および補助トルクを生成する補助トルク機構５４０を有する。ＥＰＳシステム２０００は、運転者がステアリングハンドルを操作することによって発生するステアリングシステムの操舵トルクを補助する補助トルクを生成する。補助トルクにより、運転者の操作の負担は軽減される。  
【００７４】
ステアリングシステム５２０は、例えば、ステアリングハンドル５２１、ステアリングシャフト５２２、自在軸継手５２３Ａ、５２３Ｂ、回転軸５２４、ラックアンドピニオン機構５２５、ラック軸５２６、左右のボールジョイント５５２Ａ、５５２Ｂ、タイロッド５２７Ａ、５２７Ｂ、ナックル５２８Ａ、５２８Ｂ、および左右の操舵車輪５２９Ａ、５２９Ｂ備える。  
【００７５】
補助トルク機構５４０は、例えば、操舵トルクセンサ５４１、自動車用電子制御ユニット（ＥＣＵ）５４２、モータ５４３および減速機構５４４を備える。操舵トルクセンサ５４１は、ステアリングシステム５２０における操舵トルクを検出する。ＥＣＵ５４２は、操舵トルクセンサ５４１の検出信号に基づいて駆動信号を生成する。モータ５４３は、駆動信号に基づいて操舵トルクに応じた補助トルクを生成する。モータ５４３は、減速機構５４４を介してステアリングシステム５２０に、生成した補助トルクを伝達する。  
【００７６】
ＥＣＵ５４２は、例えば、実施形態１によるコントローラ１００および駆動回路２００などを有する。自動車ではＥＣＵを核とした電子制御システムが構築される。ＥＰＳシステム２０００では、例えば、ＥＣＵ５４２、モータ５４３およびインバータ５４５によって、モータ制御システムが構築される。そのモータ制御システムとして、実施形態１によるモータ制御システム１０００を好適に用いることができる。  
【００７７】
本開示の実施形態は、トルク角の推定能力が求められる、シフトバイワイヤ、ステアリングバイワイヤ、ブレーキバイワイヤなどのエックスバイワイヤおよびトラクションモータなどのモータ制御システムにも好適に用いられる。例えば、本開示の実施形態によるモータ制御システムは、日本政府および米国運輸省道路交通安全局（ＮＨＴＳＡ）によって定められたレベル０から４（自動化の基準）に対応した自動運転車に搭載され得る。
【産業上の利用可能性】
【００７８】
本開示の実施形態は、掃除機、ドライヤ、シーリングファン、洗濯機、冷蔵庫および電動パワーステアリングシステムなどの、各種モータを有する多様な機器に幅広く利用され得る。
【符号の説明】
【００７９】
１００：コントローラ、１１０：プレ演算ユニット、１２０：トルク角演算ユニット、１３０：位相角演算ユニット、１４０：ロータ角演算ユニット、１５０：トルク演算ユニット、１６０：モータ制御ユニット、２００：駆動回路、３００：インバータ、４００、４００Ａ、４００Ｂ：電流センサ、５００：ＡＤコンバータ、６００：ＲＯＭ、７００：位置センサ、１０００：モータ制御システム、２０００：ＥＰＳシステム
【図１】
2018159102.tif 000013
【図２】
2018159102.tif 000014
【図３】
2018159102.tif 000015
【図４】
2018159102.tif 000016
【図５】
2018159102.tif 000017
【図６】
2018159102.tif 000018
【図７】
2018159102.tif 000019
【図８】
2018159102.tif 000020
【図９】
2018159102.tif 000021
【図１０】
2018159102.tif 000022
【国際調査報告】
2018159102.tif 000023
2018159102.tif 000024
2018159102.tif 000025