JP2008141835A - モータの制御方法及びそれを利用するモータ制御装置 - Google Patents

モータの制御方法及びそれを利用するモータ制御装置 Download PDF

Info

Publication number
JP2008141835A
JP2008141835A JP2006324077A JP2006324077A JP2008141835A JP 2008141835 A JP2008141835 A JP 2008141835A JP 2006324077 A JP2006324077 A JP 2006324077A JP 2006324077 A JP2006324077 A JP 2006324077A JP 2008141835 A JP2008141835 A JP 2008141835A
Authority
JP
Japan
Prior art keywords
motor
current
axis
phase
winding
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2006324077A
Other languages
English (en)
Other versions
JP5167631B2 (ja
Inventor
Masayuki Nashiki
政行 梨木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Denso Corp
Original Assignee
Denso Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Denso Corp filed Critical Denso Corp
Priority to JP2006324077A priority Critical patent/JP5167631B2/ja
Priority to US11/987,365 priority patent/US7772790B2/en
Priority to DE102007057499.3A priority patent/DE102007057499B4/de
Priority to CN2007101946877A priority patent/CN101192808B/zh
Publication of JP2008141835A publication Critical patent/JP2008141835A/ja
Application granted granted Critical
Publication of JP5167631B2 publication Critical patent/JP5167631B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K1/00Details of the magnetic circuit
    • H02K1/06Details of the magnetic circuit characterised by the shape, form or construction
    • H02K1/12Stationary parts of the magnetic circuit
    • H02K1/14Stator cores with salient poles
    • H02K1/145Stator cores with salient poles having an annular coil, e.g. of the claw-pole type
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K1/00Details of the magnetic circuit
    • H02K1/06Details of the magnetic circuit characterised by the shape, form or construction
    • H02K1/22Rotating parts of the magnetic circuit
    • H02K1/24Rotor cores with salient poles ; Variable reluctance rotors
    • H02K1/246Variable reluctance rotors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K3/00Details of windings
    • H02K3/04Windings characterised by the conductor shape, form or construction, e.g. with bar conductors
    • H02K3/28Layout of windings or of connections between windings
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/06Rotor flux based control involving the use of rotor position or rotor speed sensors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • H02P21/26Rotor flux based control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/08Reluctance motors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P27/00Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
    • H02P27/06Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/08Arrangements for controlling the speed or torque of a single motor
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/10Arrangements for controlling torque ripple, e.g. providing reduced torque ripple
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/09Motor speed determination based on the current and/or voltage without using a tachogenerator or a physical encoder
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Control Of Electric Motors In General (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

【課題】モータの磁気的な非線形性、磁気飽和に起因してモータのインダクタンスがその電流値により大きく変動する。その結果、モータの制御誤差が発生し、制御精度が劣化する。
【解決手段】モータをコンピュータを用いた非線形有限要素法などで解析し、各電流値における各巻線の磁束鎖交数を求め、この磁束鎖交数のデータテーブルを備え、このデータデーブルの情報を使用してモータの電圧制御を行い、モータの制御パラメータの適正化をモータ制御中にオンラインで制御することにより高精密なモータ制御を実現する。
【選択図】図4

Description

本発明は、自動車やトラック等に搭載されるモータの制御方法に関する。
従来から、ステータ磁極に各相のコイルが全節巻き、分布巻きに巻回されたモータが知られている。図34は、このような従来のモータの概略的な構成を示す縦断面図である。1はロータ軸、2はロータ、3は軸受け、4はステータコア、5はステータ巻線のコイルエンド部、6はモータのケースである。また、図37は図34のAA−AA線断面図である。これらの図には、4極、24スロットのシンクロナスリラクタンスモータが示されており、ステータの巻線構造は全節巻き、分布巻きであって、各ステータ磁極には各相のコイルが電気角で180°のピッチで分布状に巻回されている。35Jはステータのバックヨーク部で、35Hはステータの歯である。351、352はU相巻線で、357、358の−U相巻線と巻回されている。バックヨーク部に各巻線の相を付記している。355、356はV相巻線で、35B、35Cの−V相巻線へ巻回されている。359、35AはW相巻線で、35K、35Lの−W相巻線へ巻回されていて、残りの電気角360度については、同様の巻線である。なお、巻線の巻回方法は、種々の方法があり、同相の正符号と負符号とが巻回されていれば、前記説明の巻回方法から変更しても同一のモータ特性となる。
35Gはロータであり、電磁鋼板がロータ軸方向に積層されている。35Fはスリット状にプレス加工で打ち抜かれたフラックスバリアで、35Eはフラックスバリア35Fにより形成された細い磁路である。ロータの外周部がわずかにつながったブリッジにより、前記の細い磁路を支えている。このような構成とすることにより、ロータ座標系であるd、q軸で表現して、d軸方向の磁気インピーダンスを小さくし、q軸方向の磁気インピーダンスを大きくしている。図37に示すモータは、d、q軸の磁気インピーダンスの差によりリラクタンストルクを発生させ、駆動するモータである。
このようなモータを比較的精密に制御する制御装置のブロックダイアグラムの例を図39に示す。591はモータ、592はモータ591の回転位置θr、速度ωを検出するためのエンコーダである。593はエンコーダ592のインターフェイスであり、回転位置θr、速度ωを検出し、出力する。59Hは固定座標/回転座標変換手段であり、電流検出器で検出されたモータ591のU相電流iuとW相電流iwおよび前記回転位置θrを入力し、d軸電流idおよびq軸電流iqを計算し、出力する。
594は速度誤差検出手段であり、速度指令ωから検出速度ωを減算し、速度誤差を出力する。595は速度制御手段で速度誤差を入力し、比例、積分制御などの補償制御を行いトルク指令Tを電流指令手段596へ出力する。電流指令手段596は、トルク指令T、回転速度ωを入力としてd軸電流指令id、q軸電流指令iqを出力する。
597は電圧信号発生手段であり、d軸の電圧フィードフォワード指令FFd、q軸の電圧フィードフォワード指令FFqを出力する。59Bは電流制御ループゲインの設定手段で、d軸電流制御ループのループゲインGd、q軸電流制御ループのループゲインGqを、初期設定などで設定し、出力する。
598はd軸電流誤差検出手段であり、d軸電流指令idからd軸電流idを減算してd軸電流誤差をd軸電流制御手段599へ出力される。d軸電流制御手段599は、d軸電流誤差を比例、積分制御などの補償制御を行い、また、電流制御ループゲインGdに比例した値としてd軸電流制御電圧指令をd軸電圧制御手段59Aへ出力する。
同様に、59Cはq軸電流誤差検出手段であり、q軸電流指令iqからq軸電流iqを減算してq軸電流誤差をq軸電流制御手段59Dへ出力される。q軸電流制御手段59Dは、q軸電流誤差を比例、積分制御などの補償制御を行い、また、電流制御ループゲインGqに比例した値としてq軸電流制御電圧指令をq軸電圧制御手段59Eへ出力する。
d軸電圧制御手段59Aは加算機であり、前記d軸電流制御電圧指令とd軸の電圧フィードフォワード指令FFdとを加算し、d軸電圧指令Vdとして回転座標/固定座標座標変換手段59Fへ出力する。同様に、q軸電圧制御手段59Eは加算機であり、前記q軸電流制御電圧指令とq軸の電圧フィードフォワード指令FFqとを加算し、q軸電圧指令Vqとして回転座標/固定座標変換手段59Fへ出力する。
回転座標/固定座標変換手段59Fは、d、q軸座標でのd軸電圧指令Vdとq軸電圧指令Vqとを固定子座標のU、V、W相の3相電圧指令に変換し、3相のパルス幅変調PWMインバータ59Gへ出力し、モータ591の3相電圧、3相電流が制御される。3相インバータの主回路の具体例は、図38のような構成となっている。N95は直流電源で、N96、N97、N98N9A、N9B、N9Cは電力半導体素子であり、いわゆるIGBT、Power MOS FETなどである。N91は2相モータを表していて、N92、N93、N94は3相の巻線であり、インバータより3相電流N9D、N9E、N9Fの電流が供給される。このような制御装置の構成で、モータ591の電流、電圧、速度が制御されている。なお、前記の各ブロックで示した機能はマイクロプロセッサーを使用したソフトウェアで実現される。
図39に示すような制御装置は、例えば、特許文献1に開示されており、このように、マイクロプロセッサを使用し、回転座標系である、d、q軸座標系を基本座標として制御する制御装置が一般的に使用されている。
特開2004−289959号公報(第5頁、図1)
図37に示すシンクロナスリラクタンスモータを含め、永久磁石内蔵型モータなど種々のモータの制御は、図39に示すような制御装置で制御されることが多いが、種々の問題がある。大きな問題点の一つは、モータが理想的で線形なモータ特性ではないことに起因している。本発明の目的は、非線形な特性のモータをより正確に制御する制御装置を提供するものである。以下、現状の問題点について記述する。
モータは軟磁性体である電磁鋼板と巻線で構成され、必要に応じ、永久磁石が付加されている。モータの制御を困難にしている原因の一つは、軟磁性体の非線形性である。良く知られているように、電磁鋼板は非線形であり、磁気飽和特性となっている。そして、モータを小型化、低コスト化する必要があるため、ほとんどの場合、モータを磁気飽和する近傍を使用して制御されている。
そのようなモータの制御は、d軸インダクタンスLd、q軸インダクタンスLqなどの制御パラメータを、非線形であるにもかかわらず、ほぼ一定として制御されるため不正確であり、大きな制御的誤差が制御内部で発生するという問題がある。具体的には、図39のブロックダイアグラムにおいて、d軸電流制御手段599の適切な電流制御感度が電流の大きさにより変化することへの対応ができていないこと、q軸電流制御手段59Dの適切な電流制御感度が電流の大きさにより変化することへの対応ができていないこと、また、電圧信号発生手段597の出力であるd軸電圧フィードフォワード電圧FFd、q軸電圧フィードフォワード電圧FFqが不正確なことである。また、トルク特性が非線形であることから、電流指令手段596の電流指令id、iqが不正確であるという問題もある。その結果、モータ制御の特定領域については良好な制御が可能であるが、モータ制御の全速度領域、全トルク領域での高速、高応答な制御、正確な電圧制御ができないという問題が発生することが多い。
このような状況ではあるが、現実的には、電流検出、速度検出、位置検出などのフィードバック制御により補償して制御し、何とか目的を達成している。しかし、高速応答、高精度なモータ制御が求められる場合、あるいは、高速回転で制御する場合などでは、フィードバック制御で補償できる範囲には限界があり、問題が発生することが多い。また、ロータに永久磁石が付加された場合には、永久磁石の磁束も印加される起磁力により変化するなど、さらに、制御内部で複雑な挙動を示すことが多い。
前記の諸問題を解決するため、特許文献1に示すような種々の補償制御、モータ制御パラメータの領域ごとの設定などが行われているが、制御が複雑化する問題、対応策が個々の問題個別の対策になりがちで十分な改善が得られない問題、調整が複雑化する問題、制御装置も高価になりがちであるという問題等がある。また、モータの種類、例えば、全節巻き、分布巻きのモータといわゆる集中巻きのモータ、リラクタンスモータなどにより、制御の方法が変える必要があるという問題もある。高速回転では、漏れインダクタンスによる電圧降下に起因した制御誤差の問題、鉄損による最大効率運転の動作点が変動する問題等がある。
近年、モータ制御装置の低コスト要求および高信頼性化の要求から、ロータの位置を検出するエンコーダを排除した種々のセンサレス制御が広く使用されているが、センサレスによる位置検出誤差の問題、制御性能限界の問題がある。
本発明は上記問題の解決をその目的としている。
前記の制御装置の諸問題を解決する具体的な解決策として、モータの各電流動作点(id、iq)における各磁束鎖交数(Ψd、Ψq)を持ち、その動作点における適切な制御変数の計算を行い、より正確に制御することができる。各巻線の鎖交磁束数Ψが正確であれば、各巻線の電圧、電流が正確に制御できることになり、モータのトルク、インダクタンスなどの認識も可能となり、正確な制御を実現できる。
ロータに永久磁石を内蔵するような複雑な形状のモータの場合においても、永久磁石の有無、ロータ構造を意識することなく同一の手法で制御することができる。
前記の各磁束鎖交数(Ψd、Ψq)あるいは、このデータを元に変換したデータは、離散的に、テーブル化して格納し、使用する場合には、各電流動作点(id、iq)に応じて内挿計算して動作点のデータを、無視できる程度の誤差で、再現することができる。このテーブル化手法により、モータ特性を記録するデータのメモリ容量を許容量以下に低減することができる。
前記の各磁束鎖交数(Ψd、Ψq)のデータは、モータモデルを有限要素法などによりコンピュータで計算して求め、正確なデータを得ることができる。なお、モータを表現する座標系は広く使用されている回転子座標系である、d、q軸座標系が適切である。しかし、他の座標系で類似機能を実現することもできる。例えば、固定子座標系で3相の巻線情報を直接扱う方法は、複雑なデータの扱いとなるが、より実際のモータに近いモデルとなるのでより高精度にできる長所がある。
モータの種類としては、全節巻、分布巻きのモータ、いわゆる集中巻きのモータ、各相の巻線が概略ループ状に配置された構成のモータなど、種々モータへ類似の手法で適用することができる。
各磁束鎖交数(Ψd、Ψq)の替わりにその動作点の電流(id、iq)との比である、磁束鎖交数比例係数LLd(id1、iq1)で扱うこともできる。なお、永久磁石を含まないモータの場合は、前記磁束鎖交数比例係数LLd(id1、iq1)は、その電流動作点(id、iq)における平均インダクタンスである。
前記磁束鎖交数(Ψd、Ψq)のテーブルを利用して、モータの各トルク、回転数に応じた電流指令値(id、iq)を得るためのデータテーブルを、事前に作成することができる。また、モータを制御するときに、オンラインで、前記磁束鎖交数(Ψd、Ψq)のテーブルを利用して、モータの各トルク、回転数に応じた電流指令値(id、iq)を得ることも可能である。
前記磁束鎖交数(Ψd、Ψq)のテーブルを利用して、モータの電流動作点(id、iq)の電圧を計算し、制御装置のフィードフォワード電圧(FFd、FFq)とすることができる。また、電流動作点(id、iq)近傍の微少範囲のインダクタンスを計算し、電流制御ループのループゲイン定数(Gd、Gq)とすることができる。なお、これらのフィードフォワード電圧(FFd、FFq)、ループゲイン定数(Gd、Gq)は、前記磁束鎖交数(Ψd、Ψq)のテーブルデータから一部あるいは全てを計算し、計算したデータをテーブルに格納しておき、モータ制御時のプロセッサ負担を軽減することもできる。
データテーブルには、漏れインダクタンスあるいは漏れ磁束鎖交数データ、鉄損データなども格納しておき、より正確の制御を実現することができる。
モータの制御装置はマイクロプロセッサを使用して制御されることが多く、その制御速度の限界があり、そのサンプリング時間に比例した制御の遅れ時間Δtが発生する。この制御の遅れ時間に起因する誤差を補償するため、前記遅れ時間Δtの補正制御を行うことができる。
ロータの位置を検出するエンコーダを排除した種々のセンサレス制御において、前記磁束鎖交数(Ψd、Ψq)のテーブルを使用することにより、その動作点における固定子巻線電流、固定子巻線電圧、固定子巻線インダクタンスから、ロータ回転位置をより正確に認識することができる。
以下、本発明を適用した一実施形態のモータの制御装置について、図面を参照しながら詳細に説明する。図1は、3相、6極、36スロットのシンクロナスリラクタンスモータの横断面図である。ロータ111は、いわゆるマルチフラックスバリア型のロータで、112はロータ軸、113は空隙のフラックスバリア部、14はフラックスバリア部で形成された細い磁路である。ロータ外周部は、前記の細い磁路がそれぞれ分離されていて、いわゆるブリッジ部は無い。このような電磁鋼板がロータ軸方向に積層され、20枚に1枚の割合で非磁性ステンレス円板がロータ強化のために挿入されている。そして、これらの電磁鋼板、非磁性ステンレス円板は接着剤で相互に強固に固定されている。115はステータのバックヨーク部で、116は巻線を巻回するスロット、118は歯で、その先端部は部分的に拡大して図示しているように、軟磁性体粉末をプレス成形した部分的な歯である。
このモータは、特性を改善するため、フラックスバリア部113をできるだけ広くしてq軸インダクタンスを小さくし、ロータ外周部とステータの歯の先端部とが対向する部分の面積をできるだけ広くしてd軸方向磁気抵抗を小さくし、ステータとロータ間のエアギャップを0.13mmと小さな値としている。ステータ外形は172mm、ロータ外形は100mm、ステータコアの積厚は95mmである。図37のシンクロナスリラクタンスモータに対し、極力、d軸インダクタンスを大きくし、q軸インダクタンスを小さくする種々改良がなされている。前記の非磁性ステンレス鋼板は、電磁気的な影響は与えずにロータ構造を補強していて、ロータの軟磁性体部形状をモータ特性として有利な形状を実現するために使用されている。なお、このモータについては、電気学会論文誌D、Vol.126、No.2、pp.116-123、2006年、2月号「スリット回転子を用いたリラクタンスモータの力率・トルク向上に関する検討」に詳細を報告している。
この図1のモータは、図37などのモータに対し、力率、効率などを改善することができているが、そのd軸インダクタンスLd、q軸インダクタンスLqの非線形有限要素法FEMによる解析結果は、図2に示すように非線形な特性であり、当然、磁気飽和も発生している。横軸はモータ電流の連続定格電流に対する%表示、縦軸はインダクタンスである。また、非線形であるだけでなく、d、q軸の磁束鎖交数(Ψd、Ψq)はd軸電流idとq軸電流iqとの間の相互干渉も発生していて、大変複雑な関係となっている。
また、制御装置で電流を可変して制御する場合、その動作点近傍のインピーダンスの大きさ、すなわち、微少範囲のインダクタンスを知る必要がある。しかし、図3に示すように、電流iと磁束鎖交数Ψの関係が非線形であり、例えば、動作点が462の点であるとき、平均インダクタンスLaveは線464の勾配であり、微少範囲のインダクタンスLstは線463の勾配である。したがって、磁束鎖交数Ψ、巻線電圧v、トルクT等を求めるときには平均インダクタンスLaveが必要であり、電流制御を正確に行うためには微少範囲のインダクタンスLstが必要である。そして、動作点が変化すれば、Lave、Lstが変化し、電流id、iqの組み合わせによりそれぞれの値が相互干渉し、変化するという大変複雑な関係となっている。
また、モータのこのような磁気的非線形性は、巻線電流による起磁力の大きさと関係があるため、モータサイズとも関係することになる。概略的な話として、モータ直径が100mm以下の小さなモータの場合は巻線の最大アンペア・ターン数はそれほど大きくならず、極端な磁気飽和を考慮せずにモータを運転することができる。一方、モータ直径が150mm以上のモータの場合は巻線の最大アンペア・ターン数がモータの磁気回路を十分に磁気飽和する程度の量となり、磁気的な非線形性の問題が顕著となることがある。
(この実施例の制御方法の概略説明)
この実施例の制御方法は、このように複雑な挙動を示すモータの制御を、理論的にはほとんど誤差無く制御する方法、制御装置を提供するものである。その概略的な考え方は次の通りである。制御対象のモータについて、オフラインで、
(1)モータの各電流動作点(id、iq)における磁束鎖交数(Ψd、Ψq)をコンピュータを使用した非線形有限要素法などの解析技術によりもとめる。
(2)各電流動作点(id、iq)における磁束鎖交数(Ψd、Ψq)を表すデータテーブルData T1を作成する。例えば、図10に示すような各電流動作点(id、iq)に対する磁束鎖交数をデータテーブルとして表したものである。
(3)トルク指令Tおよび回転数ωに対応する電流指令値(id、iq)を表すデータテーブルData T2を作成する。
(4)モータを制御する制御装置において、オンラインで、トルク指令Tおよび回転数ωに対応した電流指令値(id、iq)を、データテーブルData T2を使用して、内挿計算して求める。
(5)電流指令値(id、iq)におけるフィードフォワード電圧(FFd、FFq)を求める。
(6)電流指令値(id、iq)における電流制御ループのループゲイン定数(Gd、Gq)を求める。
(7)その他については従来と同様の、モータの電圧、電流等の制御を行う。このような概略手順により、モータの各動作点における適切な制御を行うことにより、高速応答、高精度な制御を実現することができる。また、高速回転に対しても適切な制御が可能であり、ロータに永久磁石を含む場合においても永久磁石の有無を意識することなく同一の手法で制御することができる。それぞれの具体的な方法について、以下に説明する。
(図4のブロックダイアグラムの概略説明)
図4は本発明制御装置のブロックダイアグラムである。各制御信号である、変数、パラメータ等の演算は、マイクロプロセッサーのソフトウェアによって実行される。また、ゲートアレイなどのハードウェアによって実現することも可能である。図39に示す従来の制御装置と同一の構成要素については、同一の符号で示す。131は各電流動作点(id、iq)における磁束鎖交数(Ψd、Ψq)を表すデータテーブルData T1である。132はトルク指令値Tと回転数ωの条件から適切な電流指令値(id、iq)を求めるためのデータテーブルData T2である。134は電流指令値(id、iq)におけるフィードフォワード電圧(FFd、FFq)を回転数ωと磁束鎖交数(Ψd、Ψq)より求める電圧信号発生手段である。135は電流指令値(id、iq)近傍の磁束鎖交数(Ψd、Ψq)の情報から電流の微少変化の範囲におけるインダクタンスから電流制御ループのゲイン(Gd、Gq)を求める。
136はd軸電流制御手段で、d軸電流誤差信号を比例、積分制御などを行い、d軸電流制御ループのゲインGdを乗じてd軸電流制御電圧指令をd軸電圧制御手段59Aへ出力する。
137はq軸電流制御手段で、q軸電流誤差信号を比例、積分制御などを行い、q軸電流制御ループのゲインGqを乗じてq軸電流制御電圧指令をq軸電圧制御手段59Eへ出力する。
(シンクロナスリラクタンスモータのモータモデルにおける各電圧、トルク、磁束鎖交数との間の関係についての説明)
次に、図4に示した各構成要素の具体的な内容について、図1のモータを例として説明する。まず、モータの諸量は次のように表すことができる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
vdはd軸電圧、vqはq軸電圧、pは微分演算子、ωは回転角速度、Rは巻線抵抗、Tはトルク、Pnは極対数である。
モータの入力パワーPowerは、電圧と電流の積であり、数3で表され、また、数4の機械的出力ωTとモータ内部損失Plossで表される。ここで、巻線抵抗Rを無視し、モータ内部損失Plossが零であると仮定すると、数5となる。d軸磁束鎖交数Ψd、q軸磁束鎖交数Ψqは次の数6となる。
Figure 2008141835
Figure 2008141835
ここでは、これらのインダクタンスLd、Lqは、単純に各巻線の磁束鎖交数Ψと電流との比例係数であると定義する。例えば、インダクタンスを自己インダクタンスと相互インダクタンスに分類して考えることなども行わないこととする。
数1、数2に数6、数7を代入すると数8、数9が得られる。大型のモータの場合、巻線抵抗Rの値は小さいので、数8の代わりにその近似式としての数9を用いても良い。
Figure 2008141835
Figure 2008141835
Figure 2008141835
ここで,大型のモータの場合,巻線抵抗Rの値は小さく,数8の式は近似式として数9の式のように簡略化することもできる。改めて数1、数2の式について考えると,d軸インダクタンスLdとq軸インダクタンスLqが,図2に示すように非線形であって,定数ではないと言うことである。定数として扱うと,大半の電流領域で、数1、数2の式は大きな誤差を持つことになる。また,トルクTの数2、数10の式,巻線抵抗Rを無視し,モータ内部損失Plossが零であると仮定し,数1、数3、数5の式より導き出される。なお,本発明の説明では,モータのトルクを数2の式あるいは数9の式に示す近似式で説明するが,巻線抵抗R,鉄損,漏れインダクタンスを含めた厳密な計算によって求め,モータの制御情報を決定しても良く,また,それらの中間的な厳密さで計算し,制御しても良く,それらも本発明に含むものである。
したがって、d軸インダクタンスをLd(id、iq)として、電流(id、iq)の値により変化する関数であると考えて定義すれば、非線形ではあるが、数1、数2の式を正確な式として扱うことができる。同様に、q軸インダクタンスもLq(id、iq)とする。そして、従来の制御装置は、d軸インダクタンスをLd、q軸インダクタンスもLqを定数として扱うことを前提として使用されている点に根本的な問題がある。なお、モータをインダクタンスがほぼ一定である領域だけを使用する場合には、従来方法の制御装置でもある程度精度良くモータを制御することができるとも言える。
図18はモータを模式的に表現した図である。モータの入力251は電圧v、電流iであり、モータは電流iと磁束鎖交数Ψと電圧vと回転角速度ωで数11により表現できる。電圧vは磁束鎖交数Ψの時間変化率である。
Figure 2008141835
機械出力252は、回転角速度ωとトルクTである。この図18で、インダクタンスLがほぼ一定の値である領域で制御する場合には、磁束鎖交数Ψは数6、数7の式で表すことができるので、磁束鎖交数が解らなくても制御することが可能である。しかし、インダクタンスが非線形であるモータを扱う場合は、各電流動作点(id、iq)における磁束鎖交数Ψ(id、iq)、あるいは、インダクタンスL(id、iq)を知る必要がある。
ここで、磁束鎖交数Ψ(id、iq)、あるいは、インダクタンスL(id、iq)のどちらかが解ればよいのであるが、インダクタンスが定数でないのであれば、そもそも比例定数としてのインダクタンスの存在価値は低く、電圧方程式もインダクタンスを使用せず数8又は数9の式で表すことができる。トルクも数2の式ではなく、インダクタンスLを使用しない数10の式で表現することもできる。
したがって、従来とは異なるモータ開発手法であって、従来とは異なるモータ制御装置であり、モータの電流(id、iq)、磁束鎖交数(Ψd、Ψq)とが非線形有限要素法による磁場解析などで求まる場合は、数6、数7の式でインダクタンスLd(id、iq)およびLq(id、iq)を求める必要は必ずしもない。磁束鎖交数Ψd(id、iq)、Ψq(id、iq)あるいはインダクタンスLd(id、iq)、Lq(id、iq)のどちらかが解っていれば、モータの状態を表現でき、制御も可能である。磁束鎖交数Ψd(id、iq)、Ψq(id、iq)のデータテーブルの例を図10に、インダクタンスLd(id、iq)、Lq(id、iq)のデータテーブルの例を図9に示す。図9、図10に示すこれらのデータテーブルは、図4に示す制御装置のデータテーブル131に相当する。
図9のインダクタンスLd(id、iq)、Lq(id、iq)のデータテーブルの例は、横軸がd軸電流idで、縦軸がq軸電流iqとし、テーブルの各欄に各電流動作点(id、iq)におけるインダクタンスLd(id、iq)、Lq(id、iq)の値を格納している。テーブルの大きさ、数は、必要に応じて選択することができる。より精密な制御が必要な場合は、きめ細かなデータテーブルとすれば良い。また、テーブルは、現象を離散的にしか表現できないので、各電流動作点(id、iq)はテーブルの欄と欄との間の値になるので、各電流動作点(id、iq)におけるインダクタンスLd(id、iq)、Lq(id、iq)は周囲の欄のインダクタンスデータより内挿計算して、容易に得ることができる。内挿計算は単純な直線内挿計算あるいはより誤差を低減できるように、スプライン関数などのような曲線的な内挿計算を行っても良い。また、これらの推測計算は、一部、外挿計算であっても良い。また、モータの電磁気特性が比較的線形な領域については、前記の内挿計算の誤差が小さいので、前記のインダクタンステーブルは電流のピッチを粗くしても高精度なモータ制御が実現できる。逆に、電流が小さい領域では電磁鋼板のB−H特性が非線形であることが多く、電流が大きい領域においても電磁鋼板が磁気飽和しやすいのでB−H特性が非線形であることが多く、これらの領域においては前記データテーブルの電流のピッチを粗くして高精度なモータ制御が保てるようにすべきである。この工夫により、モータ制御を高精度に保ちながら、データテーブルのサイズを小さくすることができ、制御の簡素化、制御装置のコスト低減が可能である。また、前記の電磁鋼板等の軟磁性体の一般的な特性の観点から、電流の代表値は、原点領域、電流の小さい領域、磁気飽和にさしかかる領域、ある程度の磁気飽和特性となりモータ制御電流の最大値の近傍領域の4つの領域であり、前記データテーブルのサイズは、4行4列のサイズに簡素化できる。もちろん、データテーブルのサイズはそれ以上に大きい方が、高精度なモータ制御が可能であるが、内挿誤差を低減できる内挿手法との組み合わせによりデータテーブルのデータサイズを縮小することが可能である。
この結果、このテーブルを使用すれば、各電流動作点(id、iq)におけるモータの特性が正確に分かるので、モータの電圧と電流と数2の式で表されるトルクとを正確に制御できることになる。また、図10の磁束鎖交数Ψd(id、iq)、Ψq(id、iq)のデータテーブルも、図9のインダクタンスのデータテーブルと同様である。また、このデータテーブルには、各電流動作点(id、iq)における他のモータ制御に有用なデータも格納しておくこともできる。
以下、このモータ制御法の特徴の一部を説明する。
まず第1に、この実施形態では、モータ動作点における状態を磁束鎖交数Ψのd、q軸電流をパラメータとしたデータテーブルとして表現する(図10参照されたい)。たとえば、d、q軸電流を変数とした磁束鎖交数の関数である磁束鎖交数関数Ψ1、Ψ2の値とd、q軸電流の各値とのセットを各電流値ごとにデータテーブルとして保持する。
要するに、この実施形態では、磁束鎖交数Ψとしての磁束鎖交数関数Ψ1、Ψ2とd、q軸電流との関係によりモータモデルを表現し、このモータモデルを用いて必要なモータの状態量を求める。
第2に、この実施形態では、上記した磁束鎖交数Ψの値あるいは磁束鎖交数を変換した値は有限要素法解析法(FEM法)により計算した値より求める。
第3に、この実施形態では、上記モータモデルは、dq回転座標系により表示する。すなわち、磁束鎖交数Ψも電流もdq回転座標系により表示する(図4参照)。上記モータモデルは、d、q軸座標系で表現されたデータテーブルとしてもよく、静止座標系で表現したデータテーブルとしても表現してもよい。
(磁束鎖交数(Ψd、Ψq)を表すデータテーブルData T1を求める方法)
次に、図4の131に示す各電流動作点(id、iq)における磁束鎖交数(Ψd、Ψq)を表すデータテーブルData T1の具体的な求め方の例について説明する。
制御対象のモータを図1に示すモータとする。コンピュータ利用して、このモータ内の各部を通る磁束を計算、解析することができる。計算の方法は、有限要素法が主に使用されていて、本発明のように非線形な特性を対象とする場合は非線形有限要素法が有効である。有限要素法の他にも種々のコンピュータを使用した計算手法がある。例えば、モータ各部の磁気インピーダンスを計算し、各部の計算の合計として、全体および各部の磁束、トルク、電圧などを計算する手法である。この時、計算に使用するモータ使用に固有の諸定数、固有データなどをデータベースとして作成しておき、個々のモータの磁束鎖交数計算を高速に計算するできることもある。有限要素法の磁場解析の場合、計算時間が無視できないほど長くなることがあり、計算を高速にできる手法はその点で有利である。本発明では、コンピュータを使用した種々の計算手法による磁束φ、あるいは、磁束鎖交数Ψのデータ活用を対象としている。
図1の6極、36スロットのモータを図5の2極、12スロットのモータモデルとして換算して考えることができる。図1のモータは、図5のモータが3個内蔵されていると考えても良い。水平軸をd軸、首位直軸をq軸としたd、q軸座標で考えることにする。各スロットには、スロット番号である丸付き数字1〜12により示される巻線(1)〜(12)が配置されている。ステータコアの外周部に13番目のスロットから36番目のスロットの関係も付記している。そして、ステータコアのバックヨーク部に各巻線の相U、−U、V、−V、W、−Wの関係を付記している。
図5の各相の巻線151、152、153と各相の電流iu、iv、iwとロータ154の回転位置とd、q軸の位置関係を図6に模式的に示す。155はフラックスバリア部で空隙であり、156は軟磁性体部である。d軸方向へは磁気インピーダンスが小さく、q軸方向へは磁気インピーダンスが大きい。
この時、任意の電流振幅iaと電流の制御位相角θcを想定すると、その時のd軸電流id1、q軸電流iq1は、下記の数12、数13の式となる。
Figure 2008141835
Figure 2008141835
これらの電流id1、iq1は、図7のd軸巻線601、q軸巻線602に流れる電流と考えることができる。そして、これらのd軸巻線601、q軸巻線602と電流id1、iq1とロータの回転位置との関係は、図7のように表される。
(データテーブルの作成、利用についての概略説明)
このような条件において、例えば、d軸電流idを10種類、q軸電流iqを10種類だけ選択する時、全ての組み合わせの電流(id、iq)である100種類の電流(id、iq)について計算し、求めることになる。図1のモータについて、これらの電流値(id、iq)における磁束鎖交数Ψd(id、iq)、Ψq(id、iq)を非線形有限要素法等で求めることができ、その結果を図10のテーブルデータとして得ることができる。この図10のテーブルを図4に示す制御装置の121のデータテーブルData T1として使用することができる。
なお、前記の100種類の電流(id、iq)についての計算の量は、現在のコンピュータの計算速度では、比較的短時間で計算できる量であり、開発、設計上負担となる量ではない。また、モータの制御装置にそれらのデータを格納する場合において、そのデータ量は現在のメモリ容量からはほとんど問題とならない量である。また、前記の各10種類のd軸電流、q軸電流の値は、等間隔の電流値である必要はない。電磁気的により非線形により、電流値によるインダクタンスの変化が急激である領域ではより密な電流間隔とし、電流値によるインダクタンスの変化が直線的な領域では疎な電流間隔であるようなデータテーブルとすることにより、より少ないデータにより精度の高いモータの制御を実現できることになる。
なお、具体的な求め方は、後に詳述する。
なお、座標軸については、d、q軸座標で説明したが、軸数について、3軸以上の座標軸で変換して表現することができ、また、回転子座標軸で説明したが、固定子座標軸に変換して表現することもできる。例えば、各電流値をia、ib、ic・・・などの2個以上の変数とし、磁束鎖交数も2個以上の関数であるΨ1(ia、ib、ic・・・)、Ψ2(ia、ib、ic・・・)、Ψ3(ia、ib、ic・・・)・・・などとして表し、各電流値の時の値をデータテーブルとして持つものである。データテーブルの容量は増加するが、d、q軸で扱うのと同様の考え方で、モータの制御装置以外のデータテーブル作成システムなども含む制御システムを構成することができる。
特に、U、V、W相の3相モータの場合、電流はiu、iv、iwであって、iw=−iu−iwの関係である場合には、磁束鎖交数Ψu(iu、iv)、Ψv(iu、iv)、Ψw(iu、iv)として表現することができる。
図1のモータは、ロータ内部に永久磁石を含まないモータ構造であるが、図13、図52、図53、図54、図55に示すように永久磁石を内蔵した各種のモータが、モータの効率向上、小型化を実現するため、使用されている。ここで、201はステータコア、202はフラックスバリア部の空隙、203はフラックスバリア部に挟まれた細い磁路、204は永久磁石である。361はロータ軸、362はロータの軟磁性体部、363、364は永久磁石である。371はロータの軟磁性体部、372は永久磁石であ。381はロータの軟磁性体部、382は永久磁石である。391はロータの軟磁性体部、392は永久磁石である。
永久磁石を内蔵するモータの場合、数1の式は次の数14の式あるいは数15の式にように表現することができる。
Figure 2008141835
Figure 2008141835
永久磁石の磁束成分でd軸巻線に鎖交する磁束φmd、q軸巻線に鎖交する磁束φmqとして、巻線の巻回数Nssとするとき、d軸の永久磁石に起因する磁束鎖交数はΨmd=Nss×φmd、q軸の永久磁石に起因する磁束鎖交数はΨmq=Nss×φmqである。そして、磁石の向きがd軸方向へ向いている場合には、永久磁石に起因するq軸磁束鎖交数Ψmqは零となる。なお、厳密には、鎖交磁束φmd、φmqが複数回巻回されたそれぞれの巻線の位置によって少しずつ異なる値となるが、その巻線に鎖交する磁束の平均値として考えることとする。
トルク式は永久磁石成分の磁束鎖交数が付加されるので、数2の式は数16あるいは数17となる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Ψrdは永久磁石が発生する磁束鎖交数Ψmd以外のd軸の磁束鎖交数Ψであり、Ψrqは永久磁石が発生する磁束鎖交数Ψmq以外のq軸の磁束鎖交数Ψである。このように、永久磁石成分の磁束鎖交数(Ψmd、Ψmq)が付加されていると考えて、各式、121のデータテーブルData、T1を扱うことができる。
しかし、ここで、永久磁石により励起される磁束成分についても、特に分離することなく、各電流により励起される磁束と同一の磁束として扱うことも可能であり、その場合には、考え方、取り扱う式、取り扱うデータを単純化することができる。すなわち、永久磁石を内蔵するモータの場合においても、数8、数10の式で電圧(vd、vq)とトルクTを表現することができる。そして、この時、磁束鎖交数Ψd(id、iq)、Ψq(id、iq)には永久磁石により励起された磁束の磁束鎖交成分も含ませることができる。また、永久磁石に起因する磁束成分、磁束鎖交数成分をその他の磁束成分、磁束鎖交数成分と分離しない方が、有限要素法によるデータの処理計算、モータ制御における計算が簡略化できる。
これらの両手法を比較して考えると、モータに使用される永久磁石はコスト的な観点からその磁石厚みをできるだけ小さく設計するなど、磁石量を低減する傾向にある。また、磁石をロータ内部に埋め込んで構成する、いわゆる埋込磁石型同期モータの磁気回路設計においても、軟磁性体部の非線形性、磁気飽和を逆に利用して設計している現状から考え、特に磁石磁束を分離して扱う必要性がないと考えられる場合も多い。
また、従来は数14の式あるいは類似式で電圧方程式を考えてきたが、インダクタンスLd、Lqが一定値であることを想定した電圧方程式の場合、磁石に起因する磁束鎖交数Ψmを電流に起因する磁束鎖交数Ψ=L×iと分離することにより数14の式の線形方程式を構成できた。
これに対し、インダクタンスLは非線形であることを前提とし、磁束鎖交数のデータは有限要素法などのコンピュータを活用した解析のデータを活用し、磁束鎖交数Ψあるいは磁束鎖交数Ψと電流iとの比LLを非線形な表現が容易なデータテーブルで与えると言う手法により、永久磁石に起因する磁束鎖交数Ψmと電流に起因する磁束鎖交数Ψiとを分離することなく一緒にして取り扱うことが可能となったと言える。
磁束鎖交数Ψに所定関数を乗じるか加算してデータテーブルとして表現する場合について以下に説明する。つまり、データテーブルには、磁束鎖交数Ψではなく、磁束鎖交数Ψを変数とする所定の関数、更に言えば磁束鎖交数Ψに相関を有する関数の各値を記載してもよい。この関数としては、たとえばインダクタンスを示す量であるΨ/iを用いることができる。その他、磁石磁束量を定数値として上記磁束鎖交数Ψに加算してデータテーブルに記載してもよい。更に、磁束鎖交数Ψに相関を有する関数として、鎖交磁束数と巻線の巻回数とをデータテーブルに記載すれば、磁束鎖交数を実質的に表現することができる。更に、磁束鎖交数Ψに相関を有する関数として、平均磁束密度をデータテーブルに記載しても良い。その他、Li、Lmを鏡相インダクタンスで表現してもよい。
図10に示した磁束鎖交数Ψd(id、iq)、Ψq(id、iq)のテーブルは、種々の値を乗算したり、種々の値を加減算するなど、変形して使用することもできる。具体的な使用方法として、制御装置の演算において図10のデータを用いて何かを計算して制御する場合、図10のテーブルデータに何らかの計算を行った結果を新たなデータテーブルとして作成し、その新たなデータテーブルを使用した方が都度の計算は簡単になる。数6、数7の式の関係から磁束鎖交数Ψd(id、iq)、Ψq(id、iq)のテーブルを図9に示すように、インダクタンスLd(id、iq)、Lq(id、iq)のテーブルとして置き換えることもできる。
数6、数7の式で表される図9のデータテーブルのインダクタンスは、その電流値における平均インダクタンスを表している。一方、モータの電流制御においては、各電流id、iqの微小変化に対するインピーダンスであるインダクタンスの値が必要となる。その目的で、各電流値id、iqの近傍における微小変化範囲のインダクタンスを、図9のデータテーブルと同じ形態で微少電流変化に対するインダクタンスのデータテーブルとして持つことは有用である。もちろん、同一のデータテーブルとして、各電流id、iqに対して、平均インダクタンスと微少電流変化に対するインダクタンスとの両方の値をデータとして持つようにしても良い。
また磁束鎖交数は、数18、数19の式で示すように、永久磁石に起因する磁束鎖交数とその他の磁束鎖交数に分離して扱うこともできるので、磁束鎖交数Ψmd、Ψmqを定数として扱い、その他の磁束鎖交数Ψrd、Ψrqを各電流に関するデータテーブルとして作成し使用することもできる。
また、ある巻線部分の磁束鎖交数Ψは、その鎖交磁束φと巻線の巻回数NS1の積であり、
Figure 2008141835
として示される。例えばd軸巻線の磁束鎖交数Ψd(id、iq)はd軸巻線の巻き回数Nで除した鎖交磁束のデータテーブルとして作成し、使用しても良い。d軸巻線の磁束鎖交数Ψd(id、iq)は、前記鎖交磁束に巻き回数Nを乗じて簡単に求めることができる。
さらには、鎖交磁束φは平均磁束密度Bと面積Sの積B×Sとして数21の式のように表現することもできるので、磁束密度のデータテーブルを作成して、モータ制御に使用することも可能である。
Figure 2008141835
また、d、q軸のインダクタンスLd、Lqは、モータモデルの電圧方程式などの都合、制御演算の利便性などの理由で、異なる形式で表現されることもある。具体的な例は、次の式のインダクタンスLi、Lmである。
Figure 2008141835
Figure 2008141835
これらのインダクタンスLi、Lmを鏡相インダクタンスと呼ぶこともある。図9のインダクタンスのデータテーブルは、これらのインダクタンスLi、Lmで置き換えたデータテーブルとして表現することも可能である。以上説明したように、図9、図10のデータテーブルを変形することにより、モータの磁束分布の状態を種々の形式のデータテーブルで表すことができ、それらは、本発明に含まれる。
また、磁束鎖交数Ψd(id、iq)、Ψq(id、iq)のテーブルを使用して、磁束鎖交数Ψd、Ψqを電流id、iqの関数あるいは近似関数として関数式化することも可能である。例えば、次の数24、数25の式のような形態の関数である。
Figure 2008141835
Figure 2008141835
A〜I、L、L、L、Lは定数である。これらの式のように、種々の形態の関数として表現することができる。そして、磁束鎖交数Ψd、Ψqをデータテーブル化するのではなく、これらのA〜I、L、L、L、Lをモータパラメータ、制御パラメータとして得、各電流動作点における磁束鎖交数あるいはインダクタンスを計算して求めるものである。ただし、数24、数25の式等で表現される各値は、非線形有限要素法などのコンピュータを活用し手モータモデルから計算したデータを元に作成されるものであることを前提としている。そして、数24、数25の式の少なくともどちらかが磁気的に非線形である、あるいは、電流id、iqの相互の影響を受ける点が特徴的である。
次に、図10に示すような121のデータテーブルData T1のデータを使用する方法について考えてみる。このようなテーブルデータは、電流(id、iq)に関して、例えば前記のように10×10=100通りの組み合わせの離散的なデータ群である。したがって、任意の電流(id、iq)に対する磁束鎖交数Ψd(id、iq)、Ψq(id、iq)を得ることはできない。この解決策として、任意の電流(id、iq)の近傍の電流値における磁束鎖交数データから内挿計算してその電流(id、iq)の磁束鎖交数の近似値を得ることができる。内挿計算の具体的な例は、動作点のd軸電流idの値の前後の列と動作点のq軸電流iqの前後の行に該当する磁束鎖交数Ψd(id、iq)、Ψq(id、iq)の4組のデータを選択し、電流(id、iq)に比例して内挿計算する方法である。一般的な、マップ方式の離散的なデータから該当する動作点のデータを近似計算する手法が本発明にも適用できる。
次に、図1、図5、図6、図7に示すような、全節巻き、分布巻きのモータの磁束鎖交数を求める例について説明する。今、図5に示す丸付き数字で示す各巻線(1)〜(12)に鎖交する磁束量を求めることにより、各相の巻線の磁束鎖交数を求めることができる。そして、図7に示すd、q軸表現の磁束鎖交数Ψd(id、iq)、Ψq(id、iq)を求めることができる。
次に、巻線に鎖交する磁束の求め方の例について説明する。図8の611は巻線、612は巻線の片端、613は巻線の他端、614は巻線に流れる電流、615は電流614に起因する磁束とその他の起磁力に起因する磁束とを合わせた総磁束を表している。図8のモデルは、図1あるいは図5に示すような2次元形状のモータモデルであって、ロータ軸方向には同一の形状で、コアの積厚tが巻線611の片端612から他端613であるようなモータモデルをイメージしている。
この巻線611の各点には、これら全ての場によって生成されるベクトルポテンシャルAの値があり、そのベクトルポテンシャルAを巻線611に沿って片端612から他端613まで積分することにより巻線611に鎖交する磁束φが求められる。ある巻線部分の磁束鎖交数Ψはその鎖交磁束φと巻線の巻回数NS1の積で、Ψ=φ×NS1である。図8の巻線は、2次元のモータモデルを想定しているので、閉回路となっていないが、コアの積厚の単位長さあたりのモデルを有限要素法で解析しておき、コアの積厚tの値を比例倍してベクトルポテンシャルの積分値を求めることができる。
なお、閉回路構成の巻線の場合には巻線に沿ってベクトルポテンシャルAを積分すれば、その閉回路の巻線に鎖交する磁束φを求めることができ、巻回数Nsの積として磁束鎖交数Ψ=φ×Nsを求めることができる。
図8の説明と同様の方法で、図5の丸付き数字1で示す巻線(1)と鎖交する磁束φ1は、巻線丸1のベクトルポテンシャルAを丸付き数字1で示す巻線(1)に沿って線積分することにより求めることができる。また、例えば、U相巻線は、丸付き数字1で示す巻線(1)から丸付き数字7で示す巻線(7)へ巻回された巻線と、丸付き数字2で示す巻線(2)から丸付き数字8で示す巻線(8)へ巻回された巻線である。U相巻線である(1)、(7)と鎖交する磁束は、巻線(1)に沿ったベクトルポテンシャルの線積分であるφ1から巻線(7)に沿ったベクトルポテンシャルの線積分であるφ7の差であり、(φ1−φ7)である。同様に、残りのU相巻線である巻線(2)、(8)と鎖交する磁束は、巻線(2)に沿ったベクトルポテンシャルの線積分であるφ2から巻線(8)に沿ったベクトルポテンシャルの線積分であるφ8の差であり、(φ2−φ8)である。各スロットの巻回数をNsとすると、例えば、U相巻線の磁束鎖交数Ψuは、次の数26の式となる。
Figure 2008141835
(磁束鎖交数Ψd(id、iq)、Ψq(id、iq)の求め方の説明)
次に、磁束鎖交数を有限要素法などの解析データから求める具体的な方法を、次の(a)〜(i)項に示す。
(a)モータモデルと解析条件は、図1、図5、図6、図7に示す全節巻き、分布巻きのモータである。
(b)固定座標系のd、q軸電流成分の計算は、評価する動作点におけるd、q軸の電流id1 、iq1を電流振幅iaと制御位相角θcの値から〔数12〕、〔数13〕の式となる。次に、固定座標系での3相電流振幅値に換算した電流成分id2 、iq2を数27、数28により求める。
Figure 2008141835
Figure 2008141835
(c)電流成分id2 、iq2を固定子電流振幅値からU相からV相へ全電流が流れる位相の電流値に換算する。すなわち、相電流が最大の点から位相が30°進んだ電流値に換算したid3 、iq3 を数29、数30により求める。
Figure 2008141835
Figure 2008141835
(d)有限要素法で解析するために図1のモータに通電すべき各電流は、電流振幅iaと制御位相角θcより図5における3相巻線電流iu、iv、iwを数30〜数32により換算する。
Figure 2008141835
Figure 2008141835
Figure 2008141835
これらの3相電流を図1のモータモデルへ通電した状態で有限要素法解析を行い,トルクT1と各巻線の鎖交磁束φ1〜φ36を単位厚さt0当たりの解析値として求める。なお,解析条件として,ロータは図5のd,q軸の回転位置に固定している。すなわち,ロータのd軸がU相巻線が巻回された丸付き数字で示すスロット(1)、(2)の中間の向きに向いている。ロータはここで,前記鎖交磁束は各巻線Wnに沿ったベクトルポテンシャルAの積分値として求められる。なお,3次元形状巻線の場合でも,巻線に沿ってベクトルポテンシャルAを積分して鎖交磁束を求めることができる。
(e)磁束鎖交数Ψd3 ,Ψq3 の計算は,3相スター結線巻線の内,d軸を向いた2相分の直列巻線の磁束鎖交数Ψd3を求める。この時,2端子間のインダクタンスLxyの1/2がd軸インダクタンスLd=1/2×Lxyとなる関係を利用して計算し,計測する。
なお,この関係は,ロータ回転位置を変化させる時,スター結線した3相巻線の2端子間のインダクタンスが最大となる時の値Lmaxとし,最小となるときの値をLminとするとき,d軸インダクタンスLdは,Ld=Lmax/2,q軸インダクタンスLqは,Lq=Lmin/2の関係となる一般的な関係を利用している。ここでは,インダクタンスが最大,最小となる回転位置での関係を利用しているが,インダクタンスが最大,最小とならない中間的な回転位置θxにおいても,3相の各巻線のインダクタンスLu,Lv,Lwとd軸インダクタンスLd,q軸インダクタンスLqとの間には回転位置θxに応じた特定の関係があり,その中間的な回転位置θxにおける関係を利用してd軸インダクタンスLd,q軸インダクタンスLqを求めることも,やや複雑な計算となるが,可能である。その中間的な回転位置θxの関係を利用して計算する手法は本発明で開示する手法の変形であるので,本発明に含むものである。
丸付き数字で図示するスロット(9)、(10)、(11)、(12)から丸付き数字で示すスロット(3)、(4)、(5)、(6)へ数29の式のid3が通電されていると想定している。
Figure 2008141835
ここで、NSはスロットに巻回される巻き回数、tC はコア積厚、Pnは極対数6/2=3である。
なお、トルクリップル低減対策のために、ロータ形状が各磁極対ごとに異なる場合には、数34の式と同様の計算を36スロット分行い、数34の式から極対数Pnを削除する必要がある。
同様に、Ψq3(id1 、iq1)を数25により求める。丸付き数字で示すスロット(12)、(1)、(2)、(3)から丸付き数字で示すスロット(6)、(7)、(8)、(9)へ数30の式のiq3が通電されていると想定している。
Figure 2008141835
(f)固定座標系での図6の端子間インダクタンスLd3、Lq3を数36、数37により求める。
Figure 2008141835
Figure 2008141835
(g)d、q軸のインダクタンスLd、Lqは、スター結線時の端子間インダクタンスの最大値、最小値の1/2の値となる関係を利用して数38、数39から求める。
Figure 2008141835
Figure 2008141835
なお、前記手順で求められたインダクタンスLd(id1 、iq1)、Lq(id1 、iq1)の計算は、その表計算ソフトを準備しておけば、FEM解析結果のデータを入力するだけで数27〜数39の式を容易に計算することができる。そして、各電流id、iqにおけるインダクタンスLd(id1 、iq1)、Lq(id1 、iq1)を、例えば図9に示すデータテーブルとして得ることができる。
(dq軸における磁束鎖交数Ψをインダクタンスを使用せず、直接計算する方法の説明)
d、q軸の磁束鎖交数Ψd(id1 、iq1)、Ψq(id1 、iq1)は次の数40、数41の式として、数38)、数39の式の値を使用して求めることができる。そして、図10のテーブルを作成することができる。
Figure 2008141835
Figure 2008141835
また、解析、計算した数34、数35の式の磁束鎖交数データΨd3(id1 、iq1)、Ψq3(id1 、iq1)から、インダクタンスLd(id1 、iq1)、Lq(id1 、iq1)を求めることなく、直接、d、q軸の磁束鎖交数Ψd(id1 、iq1)、Ψq(id1 、iq1)を求めることができる。その方法は、数38、数39)、数27)、数28)、数29、数30、数34の式を用いて求めた次の数42、数43の式で求めることもできる。
Figure 2008141835
Figure 2008141835
以上、説明した方法で、数27から数43の式を使用することにより、有限要素法で磁場解析したデータであるベクトルポテンシャルA、磁束φを使用して、各電流(id1 、iq1)における各磁束鎖交数Ψd(id1 、iq1)、Ψq(id1 、iq1)、あるいは、インダクタンスLd(id1 、iq1)、Lq(id1 、iq1)を求めることができる。そして、例えば、d軸電流とq軸電流の組み合わせが100通りある場合、100種類の電流の組み合わせについて有限要素法解析を行い、磁束鎖交数を計算し、図10の磁束鎖交数の10行、10列のデータテーブルを作成することができる。なお、この程度の計算量は、現状の非線形有限要素法の解析技術、一般的なコンピュータ計算速度では大きな負担とならない程度の、実用上問題とならない程度の計算量である。また、このデータテーブルをメモリ上でどのように記憶するか、そして、モータの制御を行うマイクロプロセッサーがそれらのデータをどのような方法で扱うかについては、特に限定するものではなく、従来から一般的に使用されている種々のソフトウェア技術で実現することができる。
(3相集中巻きの磁束鎖交数の求め方の説明)
図36に示すモータは、3相、4極、6スロットのモータであり、短節巻き、集中巻き、ノンオーバラッピング巻きのモータである。ノンオーバラッピング巻きとは、3相のU、V、W巻線が相互に分離された構成であって、各相の巻線が3相全節巻きの巻線の様に物理的に重ならないことを意味している。TBU1、TBU2はU相の歯、TBV1、TBV2はV相の歯、TBW1、TBW2はW相の歯であり、344はステータコアのバックヨークである。これらの歯のロータに対向する面の円周方向形状を直線展開した図を図35に示す。横軸を電気角で示している。U、V、Wはモータの3相端子で、Nは3相巻線をスター結線した中性点である。各相の巻線を各相の歯に集中的に巻回している。345はロータの永久磁石で、表面磁石型ロータの例を図示している。
この図36に示すようなモータにおいても、前記の図5のモータモデルの磁束鎖交数を求めた方法と類似した方法で磁束鎖交数のデータテーブルを作成することができる。図23のモータモデルは、図36のモータを2極のモータで表現したモデルである。301はU相のステータ磁極、302はV相のステータ磁極、303はW相のステータ磁極である。丸付き数字で実際には図示されている(13)、(14)はU相の巻線、(15)、(16)tはV相の巻線、(17)、(18)はW相の巻線である。ただし、ロータ304は、図1のモータ、図5の2極化モデルと合わせるため、マルチフラックスバリア型のロータで表現している。
図23のモータモデルを図5のモータモデルと比較して説明する。図5のU相巻線である丸付き数字で示される巻線(1)と(7)および丸付き数字である巻線(2)及び(8)に相当する図23のU相巻線は丸付き数字で示す巻線(13)、(14)であり、短節巻きではあるが、電流位相、電圧位相、磁束の位相は同一の関係である。
同様に、図5のV相巻線である丸付き数字で示す巻線(5)と(11および(6)と(12)に相当する図23のV相巻線は巻線(15)、(16)である。
同様に、図5のW相巻線である丸付き数字で示す巻線(9)と(3)および(10)と(4)に相当する図23のW相巻線は丸付き数字で示す巻線(17)と(18)である。そして、上記した数27から数41の式において、数34の式を数44の式へ変形する必要がある。磁束鎖交数Ψd3の計算は、d軸を向いた2相分の巻線の磁束鎖交数Ψd3を求める。
Figure 2008141835
一方、磁束鎖交数Ψq3の計算は、図23のモータモデルでは、90°の位相差を持つステータの歯、巻線が存在しないので、求める方法を等価となるように変形する工夫が必要である。種々の等価な方法が考えられるが、その1方法を示す。図6のモータモデルは、図23のモータモデルのロータ回転位置を90°進め、このロータ位置にしたがってd、q軸の位置も修正し、U相電流、V相電流、W相電流の電流位相を電流(id1 、iq1)に対して電気角で90°進めた電流とする。この時、d、q軸における電流は(id1 、iq1)であるが、固定子巻線の各電流は位相が90°進んだことになり、各固定子電流の総合的な起磁力方向はd、q軸に対して図23と同じ方向となっている。この状態での電流を(id1x 、iq1x)として示すことにする。そして、図24で図示する電流ベクトルia、磁束鎖交数ベクトルΨaの方向となっている。この状態での有限要素法を使用した磁界解析で、磁束鎖交数Ψq3は下記の数45の式として求められる。
Figure 2008141835
ここで、NSはスロットに巻回される各相の巻き回数、tC はコア積厚、Pnは極対数4/2=2である。その他の磁束鎖交数の求め方、インダクタンスの求め方、データテーブルの作成方法などは数27から数43の式と同じ関係である。
なお、図24のこれらの操作は、磁束鎖交数Ψq3(id1x 、iq1x)をコンピュータ上で非線形有限要素法などにより計算、計測するための工夫の例である。そして、その原因は、モータの構成が図23、図24に示すように、ステータ磁極301、302、303および丸付き数字で示す各巻線(13)、(14)、(15)、(16)、(17)、(18)、(19)が電気角で120°に離散的に配置されていることである。120°の離散性があれば、位相差90°の表現は難しくなる。また、ロータ304も細い磁路とスリット状の空隙部が交互に配置された形状であり、磁気インピーダンスに離散性があるので、離散性が大きくなるとその対応が必要である。
同様の離散性の問題は、先に説明した図5のような全節巻、分布巻きのモータにおいても発生する。d、q軸で計算する場合には電気角で90°の位相差の巻線が配置されていれば都合がよく、図5のモータモデルの場合、電気角360の間に12個のスロットと巻線が配置されているので、不都合無く説明することができた。同様に、スロット数が4の整数倍であれば、90°位相差の巻線が存在することになり、都合良く磁束鎖交数を計算、計測することができる。る。しかし、例えば、スロット数が6個の場合、あるいは、18個の場合には、丁度90°の位相差の巻線はなく、d軸成分の鎖交磁束とq軸成分の鎖交磁束とを同時に計算、計測し、表現することは難しく、何らかの工夫が必要である。
この問題に対応する具体策の一つは、図23、図24の前記説明と同じ手法である。ロータの回転位置を90度回転させ、3相電流の位相も90°回転することによりd軸電流、q軸電流を同一とし、該当する位置の巻線の鎖交磁束および磁束鎖交数を計算することができる。
次にその他の手法について示す。例えば図25に示すように、電気角360°の間のスロット数が6個で6巻線(W1)、(W2)、(W3)、(W4)、(W5)、(W6)の場合、各巻線は60°の離散性がある。図25に示すように、電流iaを通電するとき、磁束鎖交数Ψaが発生していると仮定して、d、q軸それぞれの磁束鎖交成分を計算し、計測する。d軸成分の磁束鎖交数Ψd3(id1 、iq1)は、図25のモータモデルと下記の数46の式で計算する。
Figure 2008141835
次に、q軸成分の磁束鎖交数Ψd3(id1 、iq1)は、図25のモータモデルに該当する位置に巻線が無く、図5で説明した方法では、同等の条件で計算、計測することができない。そこで、図26に示すように、q軸成分の磁束鎖交数の計算ができるように、ロータ回転位置を30°回転させ、3相電流の位相も30°回転することによりd軸電流、q軸電流を同一とし、該当する位置の巻線の鎖交磁束を計算することができる。この状態での電流を(id1Y 、iq1Y)として示すことにする。磁束鎖交数Ψq3は下記の数47の式として求められる。その他の磁束鎖交数の求め方、インダクタンスの求め方、データテーブルの作成方法などは数27から数43の式と同じ関係である。
Figure 2008141835
また、前記のd、q軸における磁束鎖交数あるいはインダクタンスLd、Lqを求める方法は、リラクタンス型のモータのロータ回転位置を変化させる時、スター結線した3相巻線の2端子間のインダクタンスが最大となる時の値Lmaxとし、最小となるときの値をLminとするとき、数48、数49の関係となることを利用している。
Figure 2008141835
Figure 2008141835
しかし、これらの固定子座標系である3相巻線の磁束鎖交数あるいはインダクタンスは、回転子座標系であるd、q軸巻線の磁束鎖交数Ψd、ΨqあるいはインダクタンスLd、Lqと特定の関係があり、その関係を利用して、種々のロータ回転位置での各固定子巻線の鎖交磁束データを活用して、効果的に、d、q軸巻線の磁束鎖交数Ψd、ΨqあるいはインダクタンスLd、Lqを計算、計測することもできる。それらの変形した手法も本発明に含むものである。
また、トルクリップル低減対策のために、ロータ形状が各磁極対ごとに異なる場合には、数34の式と同様の計算を、各磁極対の該当する各巻線のそれぞれについて行い、各磁極対の該当する磁束鎖交数を足し合わせ、一方、数34、数35の式から極対数Pnを削除する必要がある。
また、図5、図23のモータ構成について説明したが、本発明は種々のモータ構成に適用することができる。図23に示すマルチフラックスバリア型のロータは、3相、集中巻きのステータとの組み合わせのモータ特性として、トルクリップルが大きくなる傾向にあり、その点では問題がある。一方、図36の表面磁石型のロータと図23の3相、集中巻きのステータとの組み合わせのモータ特性では、比較的容易にトルクリップルを低減できる。モータの各構成の組み合わせは、トルクリップル、界磁弱め等による定出力特性、より簡素で低コストなどの特徴があり、それぞれの用途により使い分けられている。ロータの種類については、図36の表面磁石型のモータを始め、図52〜図56のような各種のロータへ適用することができる。相数、極数が変わる場合には、各式を等価的に変形して適用することができる。
(ループ状巻線のEXM-Crossモータにおける磁束鎖交数の求め方の説明)
次に、図40に示すモータ150は、3相交流で動作する8極モータであり、ロータ11、永久磁石12、ステータ14を含んで構成されていて、特徴的には、円周方向に、概略、ループ形状の巻線15、16、17、18を持つモータである。
ロータ11は、表面に配置された複数の永久磁石12を備えている。これらの永久磁石12は、表面に沿って円周方向にN極とS極とが交互に配置されている。図41は、ロータ11の円周方向展開図である。横軸は機械角を示しており、機械角で360°の位置は電気角で1440°となる。
ステータ14は、それぞれ4個のU相ステータ磁極19、V相ステータ磁極20、W相ステータ磁極21を備えている。各ステータ磁極19、20、21は、ロータ11に対して突極状の形状を有している。図43は、ロータ11側から見たステータ14の内周側形状の展開図である。4個のU相ステータ磁極19は同一円周上に等間隔に配置されている。同様に、4個のV相ステータ磁極20は同一円周上に等間隔に配置されている。4個のW相ステータ磁極21は同一円周上に等間隔に配置されている。4個のU相ステータ磁極19をU相ステータ磁極群、4個のV相ステータ磁極20をV相ステータ磁極群、4個のW相ステータ磁極21をW相ステータ磁極群と称する。また、これらの各ステータ磁極群の中で、軸方向に沿って端部に配置されたU相ステータ磁極群とW相ステータ磁極群を端部ステータ磁極群、それ以外のV相ステータ磁極群を中間ステータ磁極群と称する。
また、U相ステータ磁極19、V相ステータ磁極20、W相ステータ磁極21のそれぞれは、互いに軸方向位置と周方向位置がずらして配置されている。具体的には、各ステータ磁極群は、相対的に機械角で30°、電気角で120°の位相差となるように互いに円周方向にずらして配置されている。図43に示す破線は、対向するロータ11の各永久磁石12を示している。同極のロータ磁極(N極に永久磁石12同士あるいはS極の永久磁石12同士)のピッチは電気角で360°であり、同相のステータ磁極のピッチも電気角で360°である。
ステータ14のU相ステータ磁極19、V相ステータ磁極20、W相ステータ磁極21のそれぞれの間には、U相巻線15、V相巻線16、17、W相巻線18が配置されている。図45は、各相の巻線の円周方向展開図を示す図である。U相巻線15は、U相ステータ磁極19とV相ステータ磁極20との間に設けられており、周方向に沿ったループ形状を成している。ロータ11側から見て時計回り方向の電流を正とすると(他の相の相巻線についても同様とする)、U相巻線15に流れる電流Iu は負(−Iu )となる。同様に、V相巻線16は、U相ステータ磁極19とV相ステータ磁極20との間に設けられており、周方向に沿ってループ形状を成している。V相巻線16に流れる電流Iv は正(+Iv )となる。V相巻線17は、V相ステータ磁極20とW相ステータ磁極21との間に設けられており、周方向に沿ったループ形状を成している。V相巻線17に流れる電流Iv は負(−Iv )となる。W相巻線18は、V相ステータ磁極20とW相ステータ磁極21との間に設けられており、周方向に沿ったループ形状を成している。W相巻線18に流れる電流Iw は正(+Iw )となる。これら3種類の電流Iu 、Iv 、Iw は、3相交流電流であり、互いに位相が120°ずつずれている。また39は軸方向起磁力を打ち消すための巻線である。
次に、ステータ14の各相ステータ磁極形状と各相巻線形状の詳細について説明する。図42は、図40のステータ14の断面箇所を示す図であり、図42(a)にはAA−AA線断面図が、図42(b)にはAB−AB線断面図が、図42(c)にはAC−AC線断面図がそれぞれ示されている。これらの図に示すように、U相ステータ磁極19、V相ステータ磁極20、W相ステータ磁極21のそれぞれは、ロータ11に対して突極形状を成しており、それぞれが相対的に機械角で30°、電気角で120°の位相差を有するような位置関係となるように配置されている。
図44は、U相巻線15の概略的な形状を模式的に示す図であり、正面図と側面図がそれぞれ示されている。ただし、巻線形状は、ステータ磁気回路の漏れ磁束、磁気飽和などを低減するために蛇行したループ状の巻線など種々変形も可能である。U相巻線15は、巻き始め端子Uと巻き終わり端子Nを有している。なお、同様に、V相巻線16、17は巻き始め端子Vと巻き終わり端子Nを有し、W相巻線18は巻き始め端子Wと巻き終わり端子Nを有している。各相巻線を3相Y結線する場合は、各相巻線15、16、17、18の巻き終わり端子Nが接続される。各相巻線15、16、17、18に流れる電流Iu 、Iv 、Iw は、各相ステータ磁極19、20、21とロータ11の永久磁石12との間でトルクを発生する電流位相に制御される。また、Iu +Iv +Iw =0となるように制御される。
次に、各相電流Iu 、Iv 、Iw とこれらの各相電流により各相ステータ磁極19、20、21に付与される起磁力との関係について説明する。図47は、エアギャップ面側(ロータ11側)から見た各相ステータ磁極19、20、21の展開図(図43)に等価的な各相電流巻線を書き加えた図である。
U相巻線は、4個のU相ステータ磁極19に同一方向で直列に巻回されている。したがって、各U相ステータ磁極19は同一方向に起磁力が付与されている。例えば、図47の左から2番目のU相ステータ磁極19に巻回されているU相巻線は、括弧付き数字で示す導線(3)、(4)、(5)、(6)によって形成されており、U相ステータ磁極19の回りにこの順番でこれらの導線が複数回巻回されている。なお、括弧付き数字で示す導線(2)、(7)は隣接するU相ステータ磁極19間の渡り線であり、電磁気的作用はない。
このようなU相巻線に流れる電流Iu の各部分について詳細に見ると、括弧付き数字で示す導線(1)と(3)の電流の大きさは同一で逆方向に流れており、起磁力アンペアターンは相殺されているため、これらの導線は等価的に電流が流れていないときと同じ状態にあるといえる。同様に、括弧付き数字で示す導線(5)と(8)の部分の電流についても起磁力アンペアターンは相殺されており、これらの導線は等価的に電流が流れていないときと同じ状態にあるといえる。このように、U相ステータ磁極19間に配置される導線に流れる電流は常に相殺されるため、電流を流す必要がなく、その部分の導線は排除することが可能である。その結果、括弧付き数字で示す導線(10)、(6)に対応するようにステータ14の円周上にループ状に流れるU相電流Iu と、括弧付き数字で示す導線(4)、(9)に対応するようにステータ14の円周上にループ状に流れるU相電流−Iu とが同時に流れている状態と同じと考えることができる。
しかも、上述した括弧付き数字で示す導線(10)、(6)に対応するようにステータ14の円周上にループ状に流れるU相電流Iu は、ステータコアの外部でループ状に流れる電流であり、ステータコアの外部は空気等であって磁気抵抗が大きいことから、ブラシレスモータ15への電磁気的作用はほとんどない。このため、省略しても影響はなく、ステータコアの外部に位置するループ状の巻線を排除することができる(なお、上述した例ではこのループ状の巻線を省略しているが、省略せずに残すようにしてもよい)。結局、図47に示すU相巻線の作用は、図40、図45に示すループ状のU相巻線15と等価であるということができる。
また、図47に示したV相巻線は、U相巻線と同様に、4個のV相ステータ磁極20を周回するように直列に巻回されている。この中で、括弧付き数字で示す導線(11)と(13)に流れる電流は大きさが同じで方向が逆であり、起磁力アンペアターンが相殺されるため、この部分は等価的に電流が流れていないときと同じ状態にあるといえる。同様に、括弧付き数字で示す導線(15)、(18)の電流についても起磁力アンペアターンは相殺されている。その結果、括弧付き数字で示す導線(20)、(16)に対応するようにステータ14の円周上に沿ってループ状に流れるV相電流Iv と、括弧付き数字で示す導線(14)、(19)に対応するようにステータ14の円周上にループ状に流れるV相電流−Iv とが同時に流れている状態と同じと考えることができる。結局、図47に示すV相巻線の作用は、図40、図45に示すループ状のV相巻線16、17と等価であるということができる。
また、図47に示したW相巻線は、U相巻線と同様に、4個のW相ステータ磁極21を周回するように直列に巻回されている。この中で、括弧付き数字で示す導線(21)と(23)に流れる電流は大きさが同じで方向が逆であり、起磁力アンペアターンは相殺されるため、この部分は等価的に電流が流れていないときと同じ状態にあるといえる。同様に、括弧付き数字で示す導線(25)、(28)の電流についても起磁力アンペアターンは相殺されている。その結果、括弧付き数字で示す導線(30)、(26)に対応するようにステータ14の円周上にループ状に流れるW相電流Iw と、括弧付き数字で示す導線(24)、(29)に対応するようにステータ14の円周上にループ状に流れるW相電流−Iw とが同時に流れている状態と同じと考えることができる。
しかも、上述した括弧付き数字で示す導線(24)、(29)に対応するようにステータ14の円周上にループ状に流れるW相電流−Iw は、ステータコアの外部でループ状に流れる電流であり、ステータコアの外部は空気等であり磁気抵抗が大きいことから、ブラシレスモータ15への電磁気的作用はほとんどない。このため、省略しても影響はなく、ステータコアの外部に位置するループ状の巻線を排除することができる。結局、図47に示すW相巻線の作用は、図40、図45に示すループ状のW相巻線18と等価であるということができる。
以上説明したように、ステータ14の各相ステータ磁極19、20、21に電磁気的作用を付与する巻線及び電流はループ状の簡素な巻線で代替えすることができ、かつ、ステータ14の軸方向両端のループ状の巻線を排除することができる。その結果、ブラシレスモータ15に使われる銅の量を大幅に低減することができるので、高効率化、高トルク化が可能となる。また、同相の円周方向のステータ磁極間に巻線(導線)を配置する必要がないため、従来構造以上の多極化が可能となり、特に巻線構造が簡素であることから、モータの生産性を向上させることができ、低コスト化が可能となる。
なお、磁気的には、U、V、W相のステータ磁極を通る磁束φu 、φv 、φw がバックヨーク部で合流し、3相交流磁束の総和が零となるφu +φv +φw =0の関係となっている。また、図35、図36に示した従来構造は、図47に示した各相突極19、20、21を2個ずつ合計6個を同一円周上に並べた構造であり、個々の突極の電磁気的作用、トルク発生はブラシレスモータ150と同じである。但し、図35、図36に示すような従来のブラシレスモータは、その構造上、図40から図46に示すモータ150のように巻線の一部を排除したり、巻線の簡素化を行うことはできない。
前記のモータ150はこのような構成を有しており、次にその動作を説明する。図48は、モータ150の電流、電圧、出力トルクのベクトル図である。X軸が実軸に、Y軸が虚軸にそれぞれ対応している。また、X軸に対する反時計回り方向の角度をベクトルの位相角とする。
ステータ14の各相ステータ磁極19、20、21に存在する磁束φu 、φv 、φw の回転角度変化率を単位電圧と称し、Eu =dφu /dθ、Ev =dφv /dθ、Ew =dφw /dθとする。各相ステータ磁極19、20、21のロータ11(永久磁石12)に対する相対位置は、図37に示したように、電気角で120°ずつシフトしているので、各相巻線15〜18の1ターンに誘起される単位電圧Eu 、Ev 、Ew は、図42に示すような3相交流電圧となる。
今、ロータが一定回転dθ/dt=S1で回転し、各相巻線15〜18の巻き回数をWu 、Wv 、Ww とし、これらの値がWc に等しいとすると、巻線15〜18の各誘起電圧Vu 、Vv 、Vw は次のように表される。なお、各ステータ磁極の漏れ磁束成分を無視すると、U相巻線の磁束鎖交数はWu ×φu、V相巻線の磁束鎖交数はWv ×φv、W相巻線の磁束鎖交数はWw ×φwである。
ここで、具体的な巻線と電圧の関係は次の数50〜数52のようになる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
U相の単位電圧Eu は、図40および図45に示されるU相巻線15の逆向きの1ターンに発生する電圧である。U相電圧Vu は、U相巻線15の逆向きに発生する電圧である。V相の単位電圧Ev は、V相巻線16の1ターンとV相巻線17の逆向きの1ターンとを直列に接続したときに両端に発生する電圧である。V相電圧Vv は、V相巻線16と逆向きのV相巻線17とを直列に接続したときの両端の電圧である。W相の単位電圧Ew は、図40および図45に示されるW相巻線18の1ターンに発生する電圧である。W相電圧Vw は、W相巻線18の逆向きに発生する電圧である。
ブラシレスモータ150のトルクを効率良く発生させようとすると、各相電流Iu 、Iv 、Iw は、各相巻線の単位電圧Eu 、Ev 、Ew と同一位相に通電する必要がある。図48では、Iu 、Iv 、Iw とEu 、Ev 、Ew とがそれぞれ同一位相であるものとし、ベクトル図の簡素化のため、同相の電圧ベクトル、電流ベクトルを同一のベクトル矢で表現している。
ブラシレスモータ150の出力パワーPa 、各相のパワーPu 、Pv 、Pw は、数53〜数56の式で表される。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
となる。また、ブラシレスモータ150の出力トルクTa 、各相のトルクTu 、Tv 、Tw は、下記の数57、数58、数59、数60で表される。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
となる。なお、本実施形態のブラシレスモータ150の電圧、電流、トルクに関するベクトル図は、図35、図36に示した従来のブラシレスモータのベクトル図と同じである。
次に、図40および図45に示した各相巻線と電流について、より高効率化する変形手法について説明する。U相巻線15とV相巻線16は、U相ステータ磁極19とV相ステータ磁極20の間に隣接して配置されたループ状の巻線であり、これらを単一の巻線にまとめることができる。同様に、V相巻線17とW相巻線18は、V相ステータ磁極20とW相ステータ磁極21の間に隣接して配置されたループ状の巻線であり、これらを単一の巻線にまとめることができる。
図46は、2つの巻線を単一の巻線にまとめた変形例を示す図である。図46と図45とを比較すると明らかなように、U相巻線15とV相巻線16が単一のM相巻線38に置き換えられ、V相巻線17とW相巻線18が単一のN相巻線39に置き換えられている。また、U相巻線15の電流(−Iu )とV相巻線16の電流(Iv )とを加算したM相電流Im (=−Iu +Iv )をM相巻線38に流すことにより、M相巻線38によって発生する磁束の状態とU相巻線15とV相巻線16のそれぞれによって発生する磁束を合成した状態とが同じになり、電磁気的に等価になる。同様に、V相巻線17の電流(−Iv )とW相巻線18の電流(Iw )とを加算したN相電流In (=−Iv +Iw )をN相巻線39に流すことにより、N相巻線39によって発生する磁束の状態とV相巻線17とW相巻線18のそれぞれによって発生する磁束を合成した状態とが同じになり、電磁気的に等価になる。
図48にはこれらの状態も示されている。図48に示されたM相巻線38の単位電圧Em 、N相巻線39の単位電圧En は以下のようになる。
Em =−Eu =−dφu /dθ
En =Ew =dφw /dθ
また、各巻線の電圧V、パワーP、トルクTのベクトル算式は以下数61〜数70で示される。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
ここで、数60の式で示されたトルク式は3相で表現され、数68の式で示されたトルク式は2相で表現されている。これらのトルク式の表現方法は異なるが、数68の式を展開すると数69の式となり、これら両式は数学的に等価であることがわかる。特に、電圧Vu 、Vv 、Vw および電流Iu 、Iv 、Iw が平衡3相交流の場合は数60の式で示されるトルクTa の値は一定となる。このとき、数68の式で示されるトルクTb は、図48に示すように、Tm とTn との位相差であるKmn=90°となる正弦波の2乗関数の和として得られ、一定値となる。
また、数68の式は2相交流モータの表現形態であり、数60の式と数70の式は3相交流モータの表現形態であるが、これらの値は同じである。しかし、数68の式において、(−Iu +Iv )の電流Im をM相巻線38へ通電する場合と−Iu とIv の電流をそれぞれU相巻線15とV相巻線16へ通電するのとでは、電磁気的には同じでも、銅損は異なる。図48のベクトル図に示すように、電流Im の実軸成分はIm にcos30°を乗じた値に減少するため、M相巻線38に電流Im を通電する方が銅損が75%になり、25%の銅損が低減されるという効果がある。
このように隣接して配置されたループ状の巻線を統合することにより、銅損が低減するだけではなく、巻線構造がさらに簡素になることから、モータの生産性をより向上させることができ、いっそうの低コスト化が可能となる。
次に、図40に示すモータのステータ14の形状に関し、そのギャップ面磁極形状の変形例について説明する。ステータ14の磁極形状は、トルク特性に大きく影響し、かつ、コギングトルクリップル、通電電流により誘起されるトルクリップルに密接に関係する。以下では、各ステータ磁極群に存在する磁束の回転角度変化率である単位電圧の形状および振幅がほぼ同一で相互に電気角で120°の位相差を維持するように、各ステータ磁極群のそれぞれに対応するステータ磁極の形状を変形する具体例について説明する。
図49は、ステータ磁極の変形例を示す円周方向展開図である。22はU相ステータ磁極、23はV相ステータ磁極、24はW相ステータ磁極である。図43に示した各相のステータ磁極22、23、24は、ロータ軸11と平行に配置された基本形状を有している。各ステータ磁極は、各相について同一形状であって、相対的に電気角で120°の位相差をなすように配置されている。このような形状を有する各ステータ磁極22、23、24を用いた場合にはトルクリップルが大きくなることが懸念される。しかし、各ステータ磁極22、23、24のラジアル方向にかまぼこ形状の凹凸を形成することにより、境界部での電磁気的作用を滑らかにすることができ、トルクリップルの低減が可能になる。また、他の方法として、ロータ11の永久磁石12の各極の表面にかまぼこ形状の凹凸を形成することにより、円周方向に正弦波的な磁束分布を実現することができ、これによりトルクリップルを低減するようにしてもよい。なお、図49の水平軸に付された角度は円周方向に沿った機械角であり、左端から右端までの1周が360°である。
また、図49に示した各相のステータ磁極22、23、24は、円周方向にスキューした形状とし、トルクリップルを低減することもできる。
ところで、図49に示したステータ磁極形状を採用した場合には、ステータ磁極のエアギャップ面形状を実現するためには、各相の巻線15、16、17、18とエアギャップ部との間にその磁極形状を実現するために各相のステータ磁極の先端がロータ軸方向に出た形状となり、軸方向に出るための磁路のスペースが必要であり、そのスペース確保のためモータ外形形状が大きくなりがちであるという問題がある。
図50は、ステータ磁極の他の変形例を示す円周方向展開図であり、この問題を軽減するステータ磁極形状が示されている。ステータ14のU相ステータ磁極28に存在する磁束φu の回転角度変化率であるU相の単位電圧をEu (=dφu /dθ)、V相ステータ磁極29に存在する磁束φv の回転角度変化率であるV相の単位電圧をEv (=dφv /dθ)、W相ステータ磁極30に存在する磁束φw の回転角度変化率であるW相の単位電圧をEw (=dφw /dθ)とするとき、各相の単位電圧Eu 、Ev 、Ew が形状、振幅がほぼ同一で、位相が相互に電気角で120°の位相差を保つように各相のステータ磁極28、29、30の形状を変形した例が図50に示されている。これらのステータ磁極形状の特徴は、各ステータ磁極28、29、30のエアギャップ面の大半がそれぞれのステータ磁極の歯の中間部分に対して距離が短く、ロータ11からの磁束が各ステータ磁極表面を通り、歯の中間部分を通り、そしてステータ14のバックヨークへの磁路を介して磁束が容易に通過できる点である。したがって、図50に示したステータ磁極形状は、図49に示したステータ磁極形状に比べて、各相巻線15、16、17、18とエアギャップ部との間のステータ磁極のスペースを小さくできることになる。その結果、ブラレスモータの外形形状を小さくすることが可能になる。
図51に示すステータ磁極の形状は、ロータに対向する面の形状が図43と図50の中間的な形状のステータ磁極である。34はU相ステータ磁極、35はV相ステータ磁極、36はW相ステータ磁極である。各相のステータ磁極のロータに面する形状は、円周方向の面積分布が正弦波分布に近い形状として、コギングトルク、トルクリップルが低減する形状としている。また、ステータ磁極の先端からステータのバックヨークまでの磁路形状をほぼ直線状とし、磁路の太さが途中で細くなって磁気飽和が起き難い形状とし、トルクの向上を図っている。その結果、ループ状巻線38、39は、図40、図46のような円環状とすることはできず、円周上を波状の形状のループ状巻線として配置している。図49、図50のモータについても、ループ状巻線の形状は、図示する各歯の形状に沿って、円周上を波状の形状のループ状巻線としても良い。
このように、ステータ磁極の歯の形状、ステータ磁極のロータに対向する部分の形状は、図43、49、50、51の例に示したように、トルクの増大、トルクリップルの低減、ステータ磁極間の漏れ磁束の低減、歯の磁気飽和の低減、製作の容易さなどの目的により種々の形状をとることができる。また、図40から図51に示したようなモータは、多相化など種々の変形も可能である。
(ループ状巻線モータの磁束鎖交数を求める方法の説明)
次に、図40に示すようなループ状巻線を持つモータの磁束鎖交数のデータテーブルを作成する方法について示す。前記したように、図40から図51に示したようなモータは、等価的に、図35、図36に示すような3相の短節巻き、集中巻き、ノンオーバラッピング巻きのモータと電磁気的なモータモデルは等価な作用をするので、先に示した図40等のいわゆる集中巻きのモータと同じ方法で磁束鎖交数を求めることができる。そして、図23、図24のモータモデルで、数44、数45の式を使用して解析し、計算すればよい。
具体的に、図40に示すようなループ状巻線を持つモータの各相巻線の磁束鎖交数Ψd3(id1 、iq1)、Ψq3(id1 、iq1)の求め方を示す。図23に示す3相集中巻きのモータにおいては、U相ステータ磁極301の方向にロータのq軸を向け、電流(id1 、iq1)において、スター結線され、逆向きに直列に巻回されたV相ステータ磁極の丸付き数字で示す巻線(15)、(16)とW相ステータ磁極の丸付き数字で示す巻線(17)、(18)との鎖交磁束をd軸方向鎖交磁束成分Ψd3(id1 、iq1)として計算し、数44の式で求めた。図40のループ状の巻線を持つモータにおいても、数44の式と等価な方法でd軸方向鎖交磁束成分Ψd3(id1 、iq1)を求めれば良いことになる。
図40のモータは、図23のモータモデルで表現すると、丸付き数字で示す巻線13の近傍の破線で示すバックヨーク部分から丸付き数字で示す巻線18の近傍の破線で示すバックヨークまでの間のバックヨークが削除されたモータモデルに相当する。その結果、丸付き数字で示す巻線(13)、(18)と鎖交する磁束φ13、φ18は零であり、モータのトルク発生上は効果がないので丸付き数字で示す巻線(13)と巻線(18)とを削除している。そして、図40のU相ステータ磁極19、V相ステータ磁極20、W相ステータ磁極21は、それぞれ、図23のU相ステータ磁極301、V相ステータ磁極302、W相ステータ磁極303に相当する関係となっている。図23の丸付き数字で示す巻線(14)、(15)、(16)、(17)は、それぞれ、図40の巻線15、16、17、18に相当する。
このような相対関係になっているので、図40のロータのq軸をU相ステータ磁極19の方向へ向け、電流(id1 、iq1)に相当する固定座標系の3相電流を固定子巻線へ通電したモータモデルで非線形有限要素法解析を行い、数71の式でd軸方向鎖交磁束成分Ψd3(id1 、iq1)を求める。
Figure 2008141835
ここで、図40の巻線16に鎖交する磁束はφ16L、巻線17に鎖交する磁束はφ17L、巻線18に鎖交する磁束はφ18Lである。
次に、図24に示すように、図40のモータのU相ステータ磁極19の方向にロータのd軸を向け、このロータ位置における電流(id1 、iq1)に相当する固定座標系の3相電流を固定子巻線へ通電したモータモデルで非線形有限要素法解析を行い、数72の式でq軸方向鎖交磁束成分Ψq3(id1 、iq1)を求める。この時、電流(id1 、iq1)は同じ値であるが、ロータ位置が電気角で90°異なるので固定子の3相電流の値は変化しており、その差を示すために電流を(id1z 、iq1z)として示すことにする。
Figure 2008141835
その他の磁束鎖交数の求め方、インダクタンスの求め方、データテーブルの作成方法などは数27〜数43の式と同じ関係である。このような非線形有限要素法の計算、磁束鎖交数の計算を繰り返すことにより、各電流値におけるd、q軸の磁束鎖交数Ψd(id 、iq)、Ψq(id 、iq)を求めることができる。
なお、図40に示すモータの各巻線は図45の展開図に示したように4個のコイルで構成されているが、図46に示したように、2個のコイルとした方が銅損が少なく、生産性も良いので、現実的には図46の展開図で示すようなコイルとすることが多い。巻線38の電流Imと巻線39の電流Inは先に説明したように、数73〜数74の式となっている。
Figure 2008141835
Figure 2008141835
したがって、回転子座標系のd、q軸電流(id1 、iq1)は、3相電流Iu、iv、iwに変換し、その後、数73、数74の式で変換してモータモデルの実電流Im、Inを計算し、有限要素法解析の電流値としても良い。
また、固定座標系の3相電流と回転座標系のd、q軸の2軸座標の電流との変換関係から、固定座標系の3相電流をもとめ、図40の3相巻線に与えるモデルとしてもほぼ等価な特性となる。これは、図45の巻線と図46の巻線は電磁気的にほぼ等価なので、図46に示すような巻線を持つモータであっても図45のモータモデルに変換して、数73、数74の式を使用せずに有限要素法解析を行うことでもある。
(FEMトルクと式算出トルクとの精度比較)
以上、各種のモータについて、図9に示すインダクタンステーブル、あるいは、図10に示す磁束鎖交数テーブルの求め方について詳細に示した。次に具体的な形状のモータについて非線形有限要素法で計算した例について示す。すなわち、(a)〜(f)に記述したインダクタンス算定方法で求めたLd(id1 、iq1)、Lq(id1 、iq1)の妥当性について、図1に示したマルチフラックスバリア型のシンクロナスリラクタンスモータトルクの例で検証する。
各動作点におけるインダクタンスの妥当性評価の方法は、非線形有限要素法で求めたトルクTfemとインダクタンスLd、Lqを使用して数2の式で計算されるトルクを比較して評価する。なお、ここでは、FEMで求めたトルクTfemがほぼ正しい値であると仮定し、トルクの検証値とする。信頼のおける非線形有限要素法で求められるトルク精度は、解析対象モデルの形状精度、材料特性と条件設定等が正確であれば十分に信頼に足るものであることは周知である。
図1のモータの各トルク特性を図11示す。電流の位相角θcは電気角で66°である。Tfemは非線形有限要素法で求めたトルク、Tは(a)〜(f)に記述したインダクタンス算定方法で求めた各電流動作点のインダクタンスLd(id1 、iq1)、Lq(id1 、iq1)より数2の式で計算したトルク、Tはd軸とq軸でそれぞれに単独に電流と磁束鎖交数を求めて得たインダクタンスLd(id1 、0)、Lq(0 、iq1)より数2の式で計算したトルク、Tは連続定格電流におけるLd 、 Lqを全電流領域で比例定数として扱って数2の式で計算したトルクである。
各トルクT 、T、Tを非線形有限要素法で求めたトルクTfemと比較して各インダクタンスの精度を検証する。トルクTはTfemとほぼ一致した値であり、電流ia=180Aの時の誤差でも3.3%と小さい値である。d軸とq軸との相互干渉を無視したTはやや誤差が大きく、電流ia=180Aの時の誤差が13.9%となった。インダクタンスLd、Lqを定数として扱うTは、特に電流の大きな領域では実際のトルクとはかけ離れた値となることが分かる。
(界磁弱め領域におけるトルク計算例)
次に、図1のモータの界磁弱め領域について、前記方法で求めたインダクタンスから計算されるトルクの精度について検証する。図12は横軸を電流位相角θc 、縦軸をトルクTとし、q軸電流iq=10A、20A、30A、40Aをパラメータとした特性である。インダクタンスから数2の式で計算したトルクT25L、T50L、T75L、T100Lは、d軸電流の微少な領域においても、非線形有限要素法で求めたトルクT25fem、T50fem、T75fem、T100femに近い値が得られており、各動作点において正確なインダクタンスが得られていることが解る。
次に、図13に示すNdFeB系の永久磁石を内蔵するモータの例について、(a)〜(f)に記述したインダクタンス算定方法で求めた各電流動作点のインダクタンスから計算したトルク計算値TAを検証する。通常、永久磁石を使用したモータの電圧方程式は、電流に起因する電圧成分と磁石に起因する電圧成分とを分けて表現することが多いが、数10の式のように、トルクは電流ベクトルiaと磁束鎖交数ベクトルΨaの外積なので、磁石に起因する磁束とその他の磁束とを分離することなく取り扱うこともできる。ただしその時は、d、q軸の全ての磁束鎖交数(Ψd、Ψq)と電流(id 、iq)の比をインダクタンス(Ld 、Lq)であると仮定して数2の式よりトルクTを求めることにする。この条件で、(a)〜(f)に記述したインダクタンス算定方法で求めた各電流動作点のインダクタンス(Ld 、Lq)から計算したトルク計算値Tは、図14に示すように、Tfemとほぼ一致した値であり、最大誤差でも2.9%と小さい値となった。この結果、永久磁石をロータに内蔵するモータについても(a)〜(f)に記述したインダクタンス算定方法、あるいは、磁束鎖交数算定方法は有効であることを確認できた。
次に、回転子座標軸の定義と回転子座標軸の回転移動について説明する。回転子座標系の一つとして、ここまで、d軸、q軸で表現して説明した。多く使用されている座標軸定義の一つは、ロータの磁極の方向をd軸とし、電気角で90°進んだ方向をq軸とする方法である。しかし、磁石内蔵型のブラシレスモータにおいては、ロータの磁極の方向が曖昧で、電流の大きさ、位相により変化するモータも多い。
(a)〜(f)に記述したインダクタンス算定方法の妥当性を検証する一つの方法として、d、q軸座標を図5の座標から図15に示すd1、q1軸座標に回転方向位置を電気角で30°変えてみる。図1のモータについて、このd1、q1座標で各軸の電流成分、磁束成分を(a)〜(f)に記述したインダクタンス算定方法で求め、数2の式で計算したトルクを図16に示す。新しい座標でのインダクタンスから計算したトルクTACWはTfemとほぼ一致した値であり、電流ia=180Aの時の誤差でも2.7%と小さい値となった。なお、図5の回転しない座標で求め、計算したトルクTを図16に書くとTACWと殆ど同じ値であり、重なった。この結果、正確な磁束鎖交数あるいはインダクタンスを求めることができれば、座標軸を回転させても大きなトルク誤差が発生しないことを検証できた。
この時、d1、q1軸の磁束鎖交数(Ψd1、Ψq1)と電流(id1、iq1)との比を(Ld1、Lq1)とした特性を図17に示す。この例では、q1軸電流iq1によりd1軸の磁束鎖交数Ψd1の大半がd1軸電流id1とは逆方向に発生するため、Ld1が負の値になっている。しかし、このLd1は電流ベクトルiaと磁束鎖交数ベクトルΨaの外積計算の一部の値としては正しく、数2の式の答であるトルクTACWは図16に示すように正確である。
また、回転子座標系でのd軸とq軸の定義は、(a)〜(f)に記述したインダクタンス算定方法に関しては、どの方向に定義しても良いことを示している。また、本発明の説明において、d、q軸のロータ方向を一定の方向に固定して説明したが、前記のように、d、q軸のロータ方向の定義を変えることもできる。
また、モータモデルの軸の定義は、2軸の直交軸とは限らず、90°以外の80°、60°などの角度で交差した軸であっても良い。すなわち、電圧v、電流i、磁束鎖交数Ψを、2つの軸にベクトルを分解して表現すればよい。また必要に応じて、各変数の各ベクトルを部分的に合成して諸量の関係式を構成することもできる。
また、モータモデルの軸の数は、2とは限らず、3、4、5などの数を選択することもできる。そして、これらの軸で構成される座標は、ロータ側あるいはステータ側のどちらでも定義できる。いわゆる、回転子座標系と固定子座標系のどちらを選択することもできる。さらには、センサレス位置検出技術で使用されているように、回転子座標でも固定子座標でもない第3の仮想座標系で諸量を定義することも可能である。
次に、モータの模式的な表現、ベクトル図、および、磁束鎖交数を使用したモータ表現について示す。図5から図17に示したように、高調波を含むFEMデータを使用して、正弦波理論でかつ集中定数的扱いをするインダクタンスLd 、 Lqあるいは磁束鎖交数Ψd 、Ψqを前記の方法で求め、数2の式あるいは数10の式でトルクTを計算し、正確な値が得られることを検証した。
現在、良く使用されているモータモデル、その電圧方程式等は数1、数2の式に示すように、モータの磁束鎖交数ΨがインダクタンスLと電流iで表現されている。図18の破線で囲まれたモータの模式図におけるΨ=L×iの関係である。しかし、インダクタンスLが非線形で、さらには永久磁石を含む特性を考える時、磁束鎖交数Ψと電流iで表現する方が汎用的であるとも考えられる。また、モータの各動作点で特性が変化するので、この磁束鎖交数Ψのデータテーブルを使用したモータのモデル化が、非線形な特性にも柔軟に対応でき適切であると考えられる。
なお、データテーブルの内容は非線形有限要素法データから正確に、かつ、容易に作成することができる。また、図18において破線で示されるモータの入力は電圧vと電流i、出力は回転速度ωとトルクTである。
モータを磁束鎖交数Ψで表現する方法は、数6、数7の式のように、Ψd=Ldid 、Ψq=Lqiqと置き換えれば良く、データテーブルは図9から図10へ置き換えられる。電圧方程式は、同様に、数1の式を数8の式のように変更して表現できる。なお、1kW以上の容量のモータで効率が90%を超える場合、抵抗値Rが相対的に小さいので数8の式の第2項を省略して数9の式として簡素化して扱っても、制御的に影響が少ない場合もある。また、トルクTは数10の式で表される。ベクトル図は図19のように変更して表現できる。
このように、数2、数25〜数28の式および図10のデータテーブル、図19のベクトル図に示すように、インダクタンスLd 、 Lqを使用せず、図18の電流i、電圧v、磁束鎖交数Ψ、回転数ωでモータを表現することもできる。
(定トルク型モータの電流指令置の設定方法の説明)
次に、モータの界磁弱め制御は行わず、モータのトルクを主に制御する場合、いわゆる定トルク型制御におけるブラシレスモータの具体的な制御の方法例について示す。まず最初に、制御対象モータについて、非線形有限要素法解析と計算により図9に示すようなインダクタンスデータテーブルあるいは図10に示すような磁束鎖交数データテーブルを作成する。そして、図4のモータ制御ブロックダイアグラムで示されるような制御を行う。そして、この制御の特徴的な部分は、トルク指令値Tから、132のデータテーブルData T2を使用して、電流指令値id、iqを決定する方法である。
あるトルクTを生成できる電流id、iqの組み合わせの値は、数10の式の関係から解るように、多く存在するので、モータの銅損が大きくならず効率を適度に維持できる領域の電流指令値id、iqを決めれば良い。具体的には、コンピュータによる非線形有限要素法により電流id、iqとトルクTとの関係を、必要なモータ動作領域について適切なピッチで評価することにより、例えば、図20の様な特性を決定することができる。これは、トルク指令値Tが決められたとき、電流指令値id、iqが、図20の矢印で示すように、一義的に決定できるような特性である。図297の特性は、種々の目的によって異なった特性とすることができる。例えば、銅損が最も小さくなる電流指令値id、iqの決定方法、目的の用途の回転角速度ωとトルクTにおいて最も損失が小さくなる電流指令値id、iqの決定方法、d軸電流指令値idを一定としてq軸電流指令値iqだけを可変する様な電流指令値id、iqの決定方法、モータの応答性を重視した電流指令値id、iqの決定方法などが考えられる。
図20はトルク指令値Tと電流指令値id、iqの特性関係が曲線の図で示されているが、現実にはマイクロプロセッサーを使用し、ソフトウェアで制御するので、図20のような特性をデータテーブルで近似的にメモリ上に持ち、ソフトウェア的に図20で説明したような機能を実現することができる。当然、有限の大きさのデータテーブルであり、データテーブルのデータ間の値については内挿計算により近似することができる。また、他の方法として、図20のような特性を近似関数化してソフトウェア的に処理し、制御することもできる。
また、図4の制御装置において、制御対象のモータの磁束鎖交数データを一式記憶しておけば、電流指令値決定などの種々制御に必要な、例えば、データテーブル132のデータテーブルを制御装置上でオンラインで作成する、あるいは、必要なデータだけをオンラインで計算することも可能であり、本発明の変形と位置付け、本発明に含むものである。
(定出力型モータの電流指令置の設定方法の説明)
次に、界磁弱め制御を行う場合の、ブラシレスモータの具体的な制御の方法例について示す。界磁弱め制御はモータの磁束量を制御するので、電源電圧の大きさが限定されている場合においても、高速回転まで駆動することができる。そして、界磁弱め制御の領域ではモータの端子電圧が電源電圧を超えないように制御し、電流は一定電流を超えないように制御すると、その時、出力は電圧と電流との積が一定になり、定出力制御を実現することにもなる。これは、モータの端子電圧を上げることなく高速回転まで定出力制御できるので、インバータのトランジスタの電流容量を過大にすることなく制御できることを意味し、インバータのサイズとコストを低減できることにもなる。
このような界磁弱め制御による高速回転領域の制御、あるいは、電圧的な側面では定電圧制御であり、あるいは、出力の側面では定出力制御であるような制御を実現する方法の例について説明する。まず最初に、制御対象モータについて、非線形有限要素法解析と計算により図9に示すようなインダクタンスデータテーブルあるいは図10に示すような磁束鎖交数データテーブルを作成し、そのモータの非線形な電磁気的特性を数値化する。そして、図4のモータ制御ブロックダイアグラムで示されるような制御を行う。そして、この制御の特徴的な部分は、トルク指令値Tから、132のデータテーブルData T2を使用して、電流指令値id、iqを決定する方法である。例えば、132のデータテーブルData T2に、図21に示すようなデータを格納する。横軸のω、ω、・・・ω、・・・ω1A、は回転角速度ωであり、縦軸のT、T・・・Tn・・・Tはトルク指令値で、その横軸と縦軸とが交叉する位置には回転角速度ωとトルクTとが所望の値であって、磁束鎖交数Ψの大きさおよびモータの電圧vがインバータの電源電圧によって決められる許容電圧値を超えないような電流(id、iq)の値が格納されている。この電流値(id、iq)は、前記の図9に示すようなインダクタンスデータテーブルあるいは図10に示すような磁束鎖交数データテーブルと数1〜数10の式などにより、前記条件であるトルク指令値T、モータの回転角速度ω、インバータの電源電圧を満たす電流値(id、iq)を計算し、図21のテーブルへ格納しておく。各電流値(id、iq)におけるインダクタンス(Ld、Lq)あるいは磁束鎖交数(Ψd、Ψq)は図9あるいは図10のテーブルの値から内挿して計算することができ、数1の式あるいは数8の式で電圧(vd、vq)を計算することができる。トルクTは、数2の式あるいは数10の式で計算することができる。なお、ここで、各電流値(id、iq)におけるトルクTは、非線形有限要素法で求められたトルクTfemの値を使用しても良い。このように、図4に示すモータ制御ブロックダイアグラムのデータテーブル132を作成することにより、モータを界磁弱め制御し、定出力制御を実現することができる。図4に示す制御において、データテーブル132以外の部分の定出力制御の方法は、例えば電圧のフィードフォワードの方法、電流制御ループのゲインの決定方法など、それぞれについて説明する本発明の制御方式が適用できる。
図1、2等に示したモータについて、非線形有限要素法でモータの磁束鎖交数テーブルを作成し、図21のような前記界磁弱め制御を適用した例を図22に示す。図22の横軸は毎分当たりの回転数、左縦軸はトルク、右縦軸は電圧である。図22の例では、2000rpm以上の領域で、モータの合成電圧が破線で示すVpa=280V以下であって、最大のトルクTが得られる電流(id、iq)の値とし、定出力特性を描いている。この時に、図21に相当するデータテーブルから計算されたトルクがTであり、その電流(id、iq)の時の非線形有限要素法により求めたトルクがTfemである。そして、モータ電圧はでVpaあり、d軸電圧はVpd、q軸電圧はVpqである。回転数が大きい領域では、理想的なモータでは小さな電圧であって欲しいq軸電圧Vpqが、d軸電圧Vpdよりも大きくなっている。このように、インバータの電源電圧より小さなモータ電圧となるように電流(id、iq)の値を決定して制御することが可能である。なお、例えば、モータのトルク指令値Tが図22のトルクTより小さい場合は、一般的には、電流(id、iq)のq軸電流iqを低減すれば、モータ電圧がVpaより小さな値に保ちながら、モータ出力トルクを所望の値に小さくすることができる。
なお、界磁弱め制御の領域の定出力制御の必要条件は、回転角速度ωとトルクTとが所望の値となり、電圧vと電流iがインバータの最大容量を超えなければよいので、電圧vと電流iと電流位相とはある程度の自由度があり、その自由度の中で選択することができる。具体的には、前記したように、銅損が最も小さくなる電流指令値(id、iq)決定方法、目的の用途の回転角速度ωとトルクTにおいて最も損失が小さくなる電流指令値(id、iq)の決定方法、d軸電流指令値idを一定としてq軸電流指令値iqだけを可変する様な電流指令値(id、iq)の決定方法、モータの応答性を重視した電流指令値(id、iq)の決定方法などが考えられる。
(磁束鎖交数Ψでトルク、電圧等を表示する方法の説明)
図1、2等に示したマルチフラックスバリア型のシンクロナスリラクタンスモータのベクトル図は、図19の様に表される。横軸はd軸で縦軸はq軸であり、回転座標系で示しており、回転子の回転位置を基準に構成したベクトル図となっている。電流はid、iqであり、合成電流はiaである。磁束鎖交数はΨd、Ψqであり、合成した磁束鎖交数はΨaである。電圧vd、vqは数1の式あるいは数8の式で表され、磁束鎖交数の電圧成分−ωΨd、ωΨqと抵抗の電圧降下成分id×R、iq×Rとで表される。磁束鎖交数の電圧成分はVaで、抵抗の電圧降下R2×iaで、合成電圧はV2である。電流位相角はθcで、モータの力率βは、数75の式で表される。
Figure 2008141835
このようにモータの変数である諸量を図19に示すベクトル図で表現することができる。ただし、モータは、一般的に、磁気的に非線形な特性なので、図19のベクトル図から外挿あるいは内挿して考えることができない動作領域もあり、注意を要する。
さらに、図19のベクトル図において、インダクタンスLd、Lqは数6、数7の式に示すように、Ψd=Ld×id、Ψq=Lq×iqとして表現できるので、このベクトル図をインダクタンスLd、Lqを使用せずに、磁束鎖交数Ψd、Ψqを使用して表現することもできる。
なお、永久磁石を使用しているモータの場合には、モータの一表現方法として、数14〜数19の式のような関係として表現することもでき、永久磁石に起因する鎖交磁束φmあるいは磁束鎖交数Ψmを分割して表現することもできる。
次に、モータの概念的な表現方法として、図18のように表現することができる。左側の矢印251はモータの入力で、電圧ベクトルvと電流ベクトルiである。右側の矢印252はモータの機械的出力であり、回転角速度ωとモータトルクTで表現している。破線の円で示した部分はモータであり、電圧ベクトルvと電流ベクトルiと磁束鎖交数ベクトルΨであり、物理的な諸量はそれらの積などの形で表現される。入力電力は、電圧ベクトルvと電流ベクトルiとの内積である。モータのトルクは電流ベクトルiと磁束鎖交数ベクトルΨとの外積である。モータの出力パワーはトルクTと回転角速度ωとのスカラー積である。なお、モータの内部損失は、図18では、省略している。
ここで、インダクタンスLは電流iと磁束鎖交数Ψとの比例定数として定義され、磁束鎖交数ΨはインダクタンスLと電流iとの積として扱われることが多かった。しかし、インダクタンスLが磁気的に非線形であったり、磁気飽和がある場合には、インダクタンスLは電流に関する複雑な関数となり、電流と磁束鎖交数との比例定数としての利便性は損なわれていることになる。すなわち、インダクタンスLが電流に関する複雑な関数となる場合には、モータ制御における各表現式、制御パラメータ、制御変数において、インダクタンスLを使用せず、磁束鎖交数Ψで直接表現し、取り扱った方がより簡素となり、モータ制御のソフトウェアの簡素化、高速化も可能となる。具体的なモータ表現式としては、数8、数15、数17、数18、数19の式などである。インダクタンスLを使用する場合には、これらの式に相当する毎回の計算において、乗算を1回多く行わなければならないことになる。
なお、後述するように、さらに精密なモータモデルとしては、漏れインダクタンス、あるいは、漏れ磁束鎖交数Ψkd、Ψkqを付加することもできる。さらに、鉄損成分も付加することができる。一般的に、高効率なモータにおいては、漏れ磁束、鉄損、機械的損失などを無視しても、モータの運転特性などの動的特性評価およびモータ運転には支障がないことが多い。しかし、モータの形態によっては、漏れ磁束鎖交数、鉄損の比率が大きく、無視できないモータも少なくない。
(電流指令値決定用のテーブルの切り替えについての説明)
次にモータの制御モードの切り替え制御について説明する。モータ及びその制御に対する要求は、用途およびモータ特性により種々である。また一つの機器におけるモータの運転においても、種々の運転特性を、その時々において求められることも少なくない。例えば、応答性が優先される高速応答型制御が求められる場合、最大効率の制御が求められる場合、それらの中間的な制御が求められる場合、高速回転での運転あるいは定出力制御が求められる場合などがある。自動車の運転に例えると、大きな加速が必要な場合は性能優先であり、一定速運転および回生運転などでは効率優先の運転が求められる。
具体的にこれらの要求に応えるためには、運転モードを認識する手段、各種運転モードを実現するためのデータテーブル等とモータ制御ソフトウェア、モータ運転モードの切り替え手段の準備が必要である。運転モードを認識する手段は、外部指示入力を認識する方法が簡単ではあるが、モータトルクTと回転各速度ωから運転の状況を判断し、運転モードを決定することもできる。例えば、トルクTが大きいときには性能優先とし、他の時には効率優先の制御とすることができる。
各種運転モードを実現するためのデータテーブルは、図20に示したような定トルク制御型のトルクTと電流iとの特性を実現するデータテーブル、図21に示したような定出力型のデータテーブルなどを複数持っておくことができる。そして、それらのデータテーブルを使用したモータ運転を実現するモータ制御ソフトウェアを準備しておく。
モータ運転モードの切り替え手段は、運転モードを認識する手段の出力により運転モードを認識し、モータを運転するソフトウエアを選択して制御すればよい。このような運転モードの切り替え運転により、その時々の要求に応えることのできるモータ運転を実現することができる。
(電圧フィードフォワード量の計算方法の説明)
次に、モータの制御精度を向上させる方法について説明する。モータを精度良く制御するためには、より正確な電圧信号を計算し、その電圧をモータへ印加させることは大変効果的である。図4に示す制御装置の例では、d軸電圧信号はFFdであり、q軸電圧信号はFFqである。これらの信号FFd、FFqは、電流指令値(id、iq)が電流指令手段133によりデータテーブル132を使用して計算された後、電圧信号発生手段134により、データテーブル131の各電流値(id、iq)における磁束鎖交数データ(Ψd、Ψq)を使用し、数8の式により計算され、数76、数77の式となる。
Figure 2008141835
Figure 2008141835
一方、モータの3相の内のU相電流iuとW相電流iwが電流検出器により検出され、固定座標系の3相電流から回転座標系のd、q軸の2相電流id、iqへ固定座標/回転座標変換手段59Hにより変換される。d軸電流指令idとd軸電流idとの差分をd軸電流誤差検出手段598により検出し、d軸電流制御手段136により比例、積分制御等の補償制御がなされd軸補正電圧信号がd軸電圧制御手段59Aへ出力され、前記d軸電圧信号はFFdと加算され、d軸電流制御電圧信号Vdとして回転座標/固定座標変換手段59Fへ出力する。
また、q軸電流指令iqとd軸電流iqとの差分をq軸電流誤差検出手段59Cにより検出し、q軸電流制御手段137により比例、積分制御等の補償制御がなされq軸補正電圧信号がq軸電圧制御手段59Aへ出力され、前記q軸電圧信号はFFqと加算され、q軸電流制御電圧信号Vqとして回転座標/固定座標変換手段59Fへ出力する。
このような方法により、d、q軸の電圧を正確に計算し、モータへその電圧を印加することにより小さな電流から大きな電流まで精密に制御することができる。この計算に使われるデータは、非線形有限要素法などにより計算された動作点の磁束鎖交数データに基づいており、磁気的の非線形性、磁気飽和、電流の過渡電圧特性なども計算された精密な値となっている。従来の手法では、ほぼ一定のインダクタンスLd、Lqの値として計算されるなどの方法がとられ、ある領域の特性に合わせてモータパラメータを調節すると他の領域では不正確になるなど、種々の問題があった。
また、前記の説明では、数76、数77の式の値を得るための電流値を電流指令値(id、iq)とした例について説明したが、この電流値を電流検出器により計測された実電流(id、iq)として、データテーブル131と数8の式により計算することもできる。この手法の場合、現実に近い電流値で計算することができるので、インバータの非線形性、電圧飽和などの現象が重なるなどより、電流指令値(id、iq)と実電流(id、iq)との差が大きくなる場合には、後者の方が優れた特性となる場合もある。また、電流指令と計測電流との中間的な電流値で制御することも可能であり、また、インバータが電圧飽和を発生している時には計測電流を使用すると言うように切り替えて制御することもできる。なお、これらの制御を行う場合は、図4のブロックダイアグラムの該当部を修正する必要がある。
これらの電圧制御は、d、q軸の相互の電圧干渉の問題を盛り込んだ制御であることから、非干渉制御と言うこともある。また、これらの制御は、マイクロプロセッサなどによりデジタル演算制御で実現されることが多く、サンプリング制御であり、サンプリング時間の遅れ時間、電流検出などの遅れ時間の問題、そして、それに伴う制御誤差が発生する。この解消策の一つとして、各制御部の時間遅れを補償する制御、遅れを推測し補償する制御が効果的である。具体的には、モータの制御のための計算において、その計算がモータへ印加される電圧、電流に反映される時間を想定し、モータの回転位置、回転速度、電流などの値の計算を行い使用するものである。これらの処理を行うことにより、予測可能な時間遅れ成分についてはその補償制御を行ったことになり、制御精度を改善することができる。
(電流ループゲイン量の計算方法の説明)
次に、モータの制御精度を、電流のフィードバックループの制御感度を各電流動作点において適正化させることにより、向上させる方法について説明する。電流のフィードバックループのループゲインが高すぎると発振現象を起こし、電流を精度良く制御できない。また、発振現象により過大電流が流れ、電力半導体を破損することもある。逆に、電流のフィードバックループのループゲインが低すぎると、発振現象は起きないが、電流の指令に対する応答性が悪くなり、電流制御誤差が大きくなり、モータの制御性能が劣化する。モータを精度良く制御するためには、より正確なモータ制御モデルであると同時に、モータの非線形性、電源電圧の変動などに対応できるように電流のフィードバックループのゲインを適正化することが重要である。前記したように、図2の示すようなインダクタンス変化を示すモータにおいては、低電流領域と最大電流の領域とではインダクタンスに10倍くらいの大きさの差が発生することもあり、このような場合には、電流制御ループゲインを10倍くらい可変する必要がある。例えば、ハイブリッド自動車において、モータによる車庫入れの自動位置決め行う場合、加速時などの大電流領域で過電流とならないように電流制御ループゲインGd 、Gqを設定すると、車庫入れ位置決めなどの小電流領域ではLd 、 Lqが数倍大きくなるためゲインが不足して応答性が劣化する問題がある。また、高速回転領域での界磁弱めによる定出力制御を行う場合、電源電圧の飽和を回避するための調整、電流の安定性を確保する調整が複雑化するという問題もある。
これらの問題を解決する例は、図4の制御装置において、電流制御ループゲインGd 、Gqを電流値(id、iq)に応じて、オンラインで常時適正な値に変更しながら制御する方法である。そうすることにより、電流フィードバック制御ループのゲインGd 、Gqを電流値(id、iq)の大きさには影響されずに制御することが可能となる。
電流制御ループのゲインGd 、Gqの具体的な計算方法は、図10に示す磁束鎖交数(Ψd、Ψq)のデータテーブルを利用する。モータを制御するとき、図4の電流指令値(id、iq)が電流指令手段133により決定された後、この電流指令値(id、iq)の近傍の磁束鎖交数(Ψd、Ψq)を図10の磁束鎖交数のデータテーブルを読み出し、この電流指令値(id、iq)の近傍の磁束鎖交数の微少範囲の電流の変化率ΔΨd/Δid 、ΔΨq/Δiqを計算する。この微少範囲の電流の変化率ΔΨd/Δid 、ΔΨ/Δiは、電流の微少範囲におけるインダクタンスLrd 、Lrqでもある。
Figure 2008141835
Figure 2008141835
そして、d軸電圧制御手段136、q軸電圧制御手段137において、電流制御ループのゲインに関するパラメータとして前記Gd 、Gqを使用する。その結果、電流の制御ループのゲインとしては、モータのインダクタンスLrd 、Lrqが電流制御ループのゲインGd 、Gqによって打ち消されることになり、電流の大きさに依存することなく電流制御ループのループゲインを一定に保つことが可能となる。その結果、小さな電流領域から大きな電流領域まで精度良く電流をフィードバック制御できることになる。
なお、インバータのスウィッチング損失を低減させる、モータの鉄損を低減させるなどの目的でインバータのパルス幅変調方式のPWM周波数を変更する場合、あるいは、インバータの微少電流領域における上下トランジスタの短絡防止用の不感帯および電流感度低下問題の補償を行う場合などにおいても、前記の電流ループゲインGd、Gqを可変することは効果があり、このような他の目的で電流ループゲインGd、Gqを可変することは、前記のモータインダクタンスの変動に対する電流制御ループゲインの可変、補償とは別の観点で行われる。
また、前記説明では、図10の磁束鎖交数テーブルを使用した例を説明したが、図9のインダクタンステーブルを使用して同様の制御を実現することができる。電流値についても、電流指令値(id、iq)を使用して電流の動作点を決定する例について説明したが、実際のモータ電流の検出値を使用して制御することも可能である。特に、電流が電源変動などの外乱等により変化が大きい場合には、電流の検出値を使用した方が現実的である。また、電流指令値と電流検出値の両方を使う、あるいは、条件により使い分けるという方法もある。
(電圧フィードフォワード量、電流ループゲイン量、あるいは、途中の計算量をd、q軸電流をパラメータとしたデータテーブルとして表現する方法の説明)
次に、モータの制御における逐次の計算量を低減し、マイクロプロセッサのモータ制御時の実時間における計算量を軽減する方法について説明する。前記のモータ制御における制御例では、図9もしくは図10のデータテーブルの値を使用して、数76〜数78の式の値を逐次計算する方法について説明した。これらの計算を前もって計算し、各電流(id、iq)における数76〜数78の式の値を図9もしくは図10に示すデータテーブルの行列に加えて記憶しておけば、モータ制御時の実時間における計算量を軽減することができる。制御実時間の計算は、電流(id、iq)の動作点におけるd、q軸の前記電圧信号FFd、FFq、および、前記電流ループゲインGd、Gqを、動作電流(id、iq)の周辺の電流値のデータテーブルの値から内挿計算して求めることができる。また、近似的には、動作電流(id、iq)の近傍の電流値のデータテーブルの値FFd、FFq、Gd、Gqで代用することもできる。
(漏れインダクタンスをモータモデルに追加する方法の説明)
次に、モータ巻線の漏れ磁束鎖交数Ψkd、Ψkqをより正確に表現し、モータを制御する方法について示す。漏れ磁束鎖交数とは、巻線に鎖交する磁束成分ではあるが、トルク発生にはほとんど作用しないコイルエンド部の磁束成分などの、無効な磁束鎖交数成分である。モータの種類により、漏れ磁束鎖交数Ψkd、Ψkqが大きいモータ、小さいモータがあり、大きいモータの場合、無視できない問題となる。
前記したモータ電圧式である数14の式は、漏れ磁束鎖交数Ψkd、Ψkqを付加し、次の数80、数81の式の様に書ける。
Figure 2008141835
Figure 2008141835
この時、漏れ磁束鎖交数Ψkd、Ψkqと漏れインダクタンスLkd、Lkqとは、数82、数83の関係となる。
Figure 2008141835
Figure 2008141835
これらの漏れインダクタンスLkd、Lkqを図7の等価回路に付加すると、図27のように示すことができる
また、数80、数81の式は下記の数84の式のようにも書ける。
Figure 2008141835
このように、漏れ磁束鎖交数Ψkd、Ψkq、あるいは、漏れインダクタンスLkd、Lkqを付加して表現し、モータ制御に使用することにより、より正確なモータ制御を実現することができる。なお、トルク式は、数2あるいは数10の式あるいは数16、数17の式で表され、式の上では、直接には、漏れ磁束に影響されないことになる。ただし、実際には、漏れ磁束部には磁気エネルギーが蓄積されることになり、力率低下を招いたり、あるいは、漏れ磁束成分によりモータ内の磁気回路が磁気飽和する問題が発生することがあり、トルクが低下する原因、モータの出力限界の原因の一つになっている。
本発明で具体的に漏れ磁束鎖交数を数80〜数84の式等を使用して制御する方法は、図10あるいは図9の各電流値(id、iq)のデータテーブルに漏れ磁束鎖交数Ψkd、Ψkqあるいは漏れインダクタンスLkd、Lkqを記憶し、使用すればよい。また、当然、図10あるいは図9とは別の、漏れ磁束鎖交数Ψkd、Ψkqあるいは漏れインダクタンスLkd、Lkqのデータテーブルを作っても良い。
(鉄損情報をモータモデルに反映する方法の説明)
次に、本発明においてより高精度化が必要な場合には、さらに、モータの固定子に発生する鉄損の計算を数14の式に付加することにより、数85の式とすることができる。
Figure 2008141835
トルク式である数16は数86の式となる。
Figure 2008141835
rmが等価鉄損抵抗で、材料と形状に固有の値であり、Krm=rm/ωである。一般に、出力Poにおける効率ηは数87の式ように表される。
Figure 2008141835
ここで、WLは電気的損失であり、数88の式により表される。
Figure 2008141835
上記した数85、数86の式も図4のモータ制御装置へ反映することにより、より高精度なモータの制御を実現することができる。具体的には、トルクTと電流(id+iq)の関係、数76、数77の電圧式、数78、数79のゲインGd 、Gqなどの各値へ影響する。
(磁束鎖交数Ψ相関情報を用いたセンサレスの位置検出/速度検出方法の説明)
次に、モータのセンサレス位置検出、センサレス速度検出への応用について説明する。モータ制御システムにおいて、エンコーダなどによりモータの回転位置、速度を検出して制御することが一般的に行われている。一方、回転位置検出器、回転速度検出器などはその信頼性、コスト、大きさ、そして、配線および制御装置側のインターフェイス回路のコスト、スペース、信頼性などの観点から、センサーレス位置検出に関する技術が盛んに研究され、実用化されている。本発明においては、モータの非線形有限要素法解析によるモータの正確な情報を図9、図10などのデータテーブルとして備えており、これらのモータ情報を効果的に使用して、より正確で効率良いセンサレス位置検出が可能となる。
センサレス位置検出およびそのモータの制御の例は、図29に示すように、センサレス検出手段59Sにより回転位置推定値θsrと回転角速度推定値ωsを推定して制御するものである。センサレス検出手段59Sは、通常、マイクロプロセッサとソフトウェアにより実現されるが、その具体的な構成は種々の方法がある。たとえば、固定子電流の零クロス点から回転子位置を推定する方法、固定子電圧の3倍調波を検出する方法、固定子電流の変化率を用いる方法、固定子磁束を間接的に制御する方法、モータ駆動電圧とは別の電圧を注入する方法などである。モータの誘起電圧成分と電流変化に伴うインダクタンスの電圧成分との両方を使用して、数14、数15の式を展開した関数として検出する、いわゆる拡張誘起電圧方式のセンサレス位置検出方式もある。詳しくは、たとえば、電気学会技術報告第1020号、31ページ及びその参考文献にも記載されている。
しかし、従来のセンサレス位置検出方法で使用されてきたインダクタンス等のモータパラメータは、一定値であったり、せいぜい、電流に関する折れ線近似の特性であった。その原因は、非線形なインダクタンスなどの電磁気情報を正確に得られなかったからである。そして、モータの正確な非線形インダクタンスなどの電磁気情報が得られなくとも高精度なロータ回転位置が検出できるように、センサレス位置検出方法が複雑化している面もある。
本発明では、磁束鎖交数をモータ全電流の関数Ψ(id、iq)として表し、かつ、これらのデータがコンピュータを利用した非線形有限要素法等で求められたデータであり、また、その記録手法、手段として、前記のデータをデータテーブルとしてコンピュータのメモリに記憶させておくことにより実現する。このデータテーブルのサイズは、モータ制御に求められる制御精度により異なるが、前記したように、現在の半導体、ハードディスクなどの高集積化技術では、ほとんど問題にならない程度のデータ量である。このような構成とすることにより、モータの運転状態を表すデータ精度きわめて高精度にすることができ、高精度なセンサレス位置検出、あるいは、高精度のセンサレス角速度検出が可能となる。また、モータの正確な電磁気的情報が得られるので、従来では難しかったような手法でも、高精度なセンサレス位置検出方法が可能となっていると言える。
また、前記のデータテーブルの手法は、非線形有限要素法等で求められたデータをモータ制御へ活用するための一手法であり、これらのデータを形式的にはデータテーブルではないが、データテーブルと等価な記憶手段で実現することも可能であり、それらも本発明に含むものである。データテーブルと等価な一つの方法は、磁束鎖交数あるいはインダクタンスのデータを数24、数25の式の例のように、データテーブルを元にした近似関数で表現する方法である。また他の方法として、前記データテーブルの数値をモータ制御プログラムの中に組み込むことも可能である。これらの方法もデータテーブルの変形として考えることができる。
(座標系d、q、α、β、UVW座標系の説明)
次に、本発明の説明で使用するモータモデルの3種類の座標について説明する。いづれも良く知られた座標である。大別して、固定子を中心として考える固定座標系と回転子を中心として考える回転座標系に分けられる。固定座標系は、図31に示すように、U、V、W軸で表す固定3相座標系とα、β軸で表す固定2相座標系を使用して表現している。図31は2極のモータモデルで、701、702、703は3相の固定子巻線であるU相巻線、V相巻線、W相巻線である。そして、これらの3相巻線はモータモデル的に等価な、破線で示す2相巻線のα巻線704とβ巻線705で置き換えることができる。一方、回転座標系は、図5、図6、図7で示したように、回転子の回転方向位置のd、q軸で表現する座標系である。なお、これらの座標系以外に変換して、等価な制御を実現することもできる。
前記のd、q軸の電流id、iq、電圧vd、vqと、U、V、W軸の電流iu、iv、iw、電圧vu、vv、vw、そして、α、β軸の電流iα、iβ、電圧vα、vβは次の数89〜数98の式のように、相関関係を持たせて表現することができる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
ここで、iaは2相の電流絶対値、ioは固定座標系の3相電流の振幅である。また、電圧関係は、電流位相と電圧位相とが同一で、力率が1の場合は次の数99〜数108のようになる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
Figure 2008141835
ここで、vaは2相の電圧絶対値、voは固定座標系の3相電圧の振幅である。数89〜数108の式に示したように、モータ制御の各変数は、座標系を変換することができ、計算、あるいは、計測の都合などにより異なる座標系で計算して変換することもできる。
次に、本発明の手法を、d、q軸以外の座標系ZXを使用してモータを制御する場合について説明する。モータの電流、電圧をd、q軸座標ではない座標系ZXで行う場合、図9あるいは図10のデータテーブルを使用する方法として、その座標系ZXにおける電流値をd、q軸の座標の電流値(id、iq)に変換して電圧、電流、トルクなどの計算をし、その結果をもとの座標系ZXに変換して戻すことができる。
また、他の方法としては、図9あるいは図10のデータテーブルを使用している座標系の電流値と磁束鎖交数とのデータテーブルDTXを、オフラインで計算して変換し、作成しておくことにより、このデータテーブルDTXを使用してモータの制御を実現することができる。例えば、α、β軸座標の場合には、電流値(iα、iβ)に関するデータテーブルDTXを作ることが可能である。ここで、固定座標系の場合は、ロータの回転位置θrの値により磁束鎖交数Ψが変化するので、3次元のiα、iβ、θrに関するデータテーブルDTXとなり、各欄にはα巻線とβ巻線の磁束鎖交数Ψα(iα、iβ、θr)、Ψβ(iα、iβ、θr)を計算し格納することになる。この結果、この方法の場合には、データテーブルDTXの大きさが大きくなる点は不利である。一方で、モータの動作として使用される範囲は限定されることが多く、そのことを利用して、データテーブルDTXの大きさを縮小することも可能である。例えば、図5において、通常使用する電流制御位相角θcの範囲が限定されている場合には、使用しない領域のデータは省略することも可能である。
U、V、W座標の場合は、例えば3相の内の一つの電流iw=−iu−ivなので自由度は2であり、電流値(iu、iv)に関するデータテーブルDTXを作っておけば良い。前記のα、β軸座標と同様に考えることができ、3次元のiu、iv、θrに関するデータテーブルDTXとなり、各欄にはu巻線とv巻線とwの磁束鎖交数Ψu(iu、iv、θr)、Ψv(iu、iv、θr)、Ψw(iu、iv、θr)を計算し格納することになる。やはり、データテーブルDTXの大きさが大きくなる点は不利であるが、前期同様、データテーブルDTXの縮小も可能である。
4相以上のデータテーブルについても、データテーブルが大きくなるが、同様に製作が可能である。磁束鎖交数のデータテーブルのデータは、各巻線の磁束鎖交数Ψである。インダクタンステーブルの場合の値は、数6、数7の式に示すように、磁束鎖交数Ψを巻線電流で割れば良い。
スウィッチトリラクタンスモータのように、正弦波電流ではないモータ電流を通電する場合においても、相当するデータテーブルを準備することにより、高精度な制御を実現することができる。すなわち、ロータの回転位置θrと各巻線の電流に対する各巻線の磁束鎖交数Ψの値がデータテーブルとして準備されていれば、モータの電圧、電流、トルクを正確に制御することができる。この時、リラクタンストルクの場合には、モータ入力の内、磁気エネルギーとモータ出力トルクとの分離の計算が、磁気非線形性などのため精度が落ちる場合もある。この対応として、データテーブルの欄に、磁束鎖交数Ψの値と共にトルクTの値を格納しておくことにより、モータの入力パワーから磁気エネルギーとモータの出力トルクを分離することなく正確なトルク知ることができるので、より高精度な制御が可能である。
また、従来の位置センサーレスモータ制御の種々方法において、図9、10に示すようなインダクタンスLd(id、iq)、Lq(id、iq)のデータテーブルあるいは磁束鎖交数Ψd(id、iq)、Ψq(id、iq)のデータテーブルを持ち、これらのデータを使用して制御することにより、より正確な電圧、電流、トルク等を制御することができるようになり、本発明に含むものである。
(他のセンサレス位置検出方法の説明)
次に、センサレス位置検出の一つの方法を、図28に示す。図29はそのモータ制御の全体を示す図で、センサレス位置検出を行うセンサレス検出手段は59Sである。このセンサレス検出手段59Sの内部構成は種々の構成で実現することができる。なお、図29のセンサレス検出手段59Sには、その入力信号等の接続関係を図示していないが、種々の制御情報を使用してモータの推定回転位置θsr、推定角速度ωsを推定計算する。図29のその他の構成は、図4に示すモータ回転位置を検出するエンコーダ592を備えるモータ制御装置と同じ構成の例を示している。前記の座標系の変形等、種々変形が可能である。
図28の661はマイクロプロセッサの割り込み処理等のスタートである。
662はモータの電圧と電流を計測する手段であり、その時点での推定回転位置θsr、推定角速度ωsにおけるモータ制御電圧Vpd、Vpqとして、図29のvd、vqを計測する。これらの値から、回転座標/固定座標変換手段59Fの出力である、固定座標系の3相電圧vu1、vv1、vw1を計算することができる。また、モータの3相電圧を測定して得ることも可能である。また、同時に3相電流iu1、iv1、iw1を計測、記憶する。663はモータ制御電圧Vpd、Vpqが計測された時点の推定回転位置θsr、推定角速度ωsを設定する。
664は、数1の式あるいは数8の式等を使用し、かつ、図9あるいは図10のデータテーブルに示すようなコンピュータにより非線形有限要素法などで正確に得られたモータ情報、前記の推定回転位置θsr、推定角速度ωs、モータの3相電流iu1、iv1、iw1を使用して、モータの推定電圧Vsd、Vsqを計算する。これらのVsd、Vsqは、3相電圧vus1、vvs1、vws1に変換しても良い。特にこの制御ステップでは、図9あるいは図10のデータテーブルの値が各電流動作点における電磁気的状態がその非線形性、磁気飽和特性も含めて正確に計算された情報なので、もし推定回転位置θsr、推定角速度ωsの誤差が小さければ、きわめて正確な計算であるところに本方式の特徴がある。そして、665は前記の推定回転位置θsr、推定角速度ωsの推定がどの程度の精度であるかを評価するために誤差評価関数δsを計算する。この誤差評価関数δsの例は、数109、数110の式に示すようにモータの電圧と電流の推定値と実際の値との差とする構成である。
Figure 2008141835
Figure 2008141835
そして、誤差判定手段666で、誤差評価関数δsの値が誤差許容値δoより小さいかどうかを判定し、誤差許容値δo以下であれば668で推定回転位置θsr、推定角速度ωsがほぼ正しい値であると判定し、669でこの計算処理を終了する。誤差評価関数δsが誤差許容値δo以上であれば、667で誤差評価関数δsの値が小さくなる方向へ、回転位置θsr、角速度ωsの値を補正計算し、再設定する。そして、経路66Aで664へもどり、電圧Vsd、Vsq等の計算を再度行い、この計算処理を誤差評価関数δsが誤差許容値δo以下となるまで繰り返す。
なお、667では、誤差評価関数δsの収束が早くなるような工夫、あるいは、計算回数が多くなった場合は終了するなどの例外処理を行っても良い。具体的な工夫の例は、誤差評価関数の値が大きければ、推定回転位置θsr、推定角速度ωsの予測値が大きく異なるので、次の演算サイクルで推定する推定回転位置θsr、推定角速度ωsの値を大きめに修正する方法がある。また、推定角速度ωsは急には大きく変化しないので、先回のモータ制御サンプリング時の値を使用するように固定することもできる。また、推定回転位置θsrを0°から360°までの間、荒いピッチで計算してみて、回転位置θrの領域を絞り込み、絞り込んだ回転位置θrの領域についてきめ細かく計算することもできる。さらには、前記の荒いピッチでの計算等は、事前に計算してデータテーブルとして記憶しておく方法も可能である。
また、現在の推定回転位置θsr、推定角速度ωsとから推定される推定値と観測されるモータ実測値との差である誤差評価関数δsの大きさに応じて、前記推定回転位置θsr、推定角速度ωsに加えるべき補正値を計算し、図28のフローチャートの経路66Aを66Bとし、計算を単純化する方法とすることもできる。すなわち、誤差評価関数δsの値が誤差許容値δo以下となり、収束するまで計算を繰り返す方法と一度の計算で推定値を決めてしまう方法とがある。また、それらの中間的な方法とすることもできる。
また、数1〜数25の式を変形してロータ回転位置θrと電圧、電流の関係式として、計測データを代入し、ロータ回転位置推定値θsrを計算する方法もある。特に、数8、数9、数14、数15の式は非線形な量であるインダクタンスLを含まないので、より正確な計算式であり、より正確な関係式を導き出すことができる。
なお、数109、数110の式では、d、q軸の電圧として換算したモータ電圧と推定回転位置θsrの仮定に基づいたモータ電圧の推定値との差から評価する例を示したが、3相電圧計測値vu1、vv1、vw1と3相電圧推定値vus1、vvs1、vws1の差で評価しても類似の結果となる。
(鎖交磁束の電流微少変化率を用いたセンサレスの位置検出/速度検出方法の説明)
次に、他の方法のモータのセンサレス位置検出、センサレス速度検出の方法について説明する。この方法は、モータ電流の微小変化に伴う電圧変化からロータ回転位置推定値θsrおよびロータ回転速度推定値ωsを求めるものである。これは、モータ電流に強制的にロータ回転位置検出用の高周波電流を重畳する、あるいは、測定するタイミングでパルス的な電流などを重畳して電流変化を与え、電圧変化を計測し、一方、図9あるいは図10のデータテーブルに示すようなコンピュータにより非線形有限要素法などで正確に得られたモータ情報を利用して推定電圧、推定電流を計算し、両値の差が小さくなるような回転位置推定値θsrと回転角速度推定値ωsを求めるものである。
モータ制御の構成例は、図29に示すように、センサレス検出手段59Sにより回転位置推定値θsrと回転角速度推定値ωsを推定して制御するものである。このセンサレス検出手段59Sの機能について説明する。
モータの磁気特性は図2に示したように、磁気的に非線形で、磁気飽和の特性もあり、特に小型、高出力、低コストなモータはほとんど磁気飽和領域を限界まで活用して設計されている。このような領域のモータ巻線の電流iとモータ巻線の磁束鎖交数Ψの関係は、図3に示すような関係となっている。モータの電流と磁束鎖交数Ψの関係は特性461であり、モータの電流動作点が462で電流がix1の時、磁束鎖交数Ψ=L×iの関係であり、平均インダクタンスLaveは直線464の勾配として示される。しかし、電流動作点462の近傍での微少範囲のインダクタンスLstは、電流動作点462の接線の勾配、すなわち、直線463の勾配であり、次の数111の式で表される。
Figure 2008141835
数111の式で示される直線463の勾配Lstと平均インダクタンスLaveである直線464の勾配とは大きく異なる値であることが解る。なお、数111の式の値をインダクタンスとして扱わず、基準電流変化Δiに対する磁束鎖交数変化ΔΨとして扱っても良い。
したがって、この方法のセンサレス位置検出において回転位置検出用の小さな高周波電流をモータ電流に重畳する場合、検出される高周波電圧は微少電流変化範囲に対する電圧であり、そのインダクタンスは平均インダクタンスLaveではなく、微少電流範囲のインダクタンスLstである。この検出方式の場合、モータの出力トルクが零でモータの実効的な電流が零である場合にもロータの回転位置θr、回転角速度ωを検出することができる。また、モータの回転数が零近傍であって、モータの誘起電圧成分、すなわち、磁束鎖交数の実効的な時間変化率dψ/dtが零近傍であっても、ロータの回転位置θr、回転角速度ωを検出することができる。
この具体的なセンサレス位置検出の一つの方法は、図30に示す方法である。
681はマイクロプロセッサの割り込み処理等のスタートである。
682は、位置検出用の高周波電流でその振幅がΔi=(ix3−ix1)/2の3相電流iu5、iv5、iw5をモータの3相電流iu4、iv4、iw4へ重畳して通電し、その時のモータの3相電圧vu4、vv4、vw4の高周波成分vu5、vv5、vw5を計測し、高周波電流の周波数を計測できるフィルターなどの処理計算を行う。そして、その時の3相電流iu4、iv4、iw4における、微少電流範囲のインダクタンスLstu4、Lstv4、Lstw4を得る。なお、高周波電流は3相高周波電流を3相巻線の電流へそれぞれ重畳する方法を示したが、重畳する高周波電流の値はU相端子からV、W相端子へ重畳して流すなど、他の方法であっても良い。また、この高周波信号は、一定周波数であれば、検出時においても信号をフィルター処理できるので、精度の高い検出が可能である。また、高周波電流の代わりに高周波電圧を重畳して、高周波電流成分を検出しても良い。また、パルス信号的な外乱注入信号でも、同様に、微少電流範囲のインダクタンスLstを検出することが可能である。683は、モータ動作点のモータ制御電圧vu4、vv4、vw4が計測された時点の推定回転位置θsr、推定角速度ωsを計算し、設定する。684は、モータ電流(id、iq)、回転位置θsr、角速度ωsの時の微少電流範囲のインダクタンスLstを計算する。この時、図9あるいは図10のデータテーブルを使用して、モータ電流(id、iq)の近傍における微少範囲のインダクタンス(Lstd、Lstq)を図3、数111の式で示した方法で、d軸成分インダクタンスLstd、q軸成分インダクタンスLstqとして計算する。なお、図9のインダクタンスのデータテーブルを使用する場合は、Ψ=L×iの関係から図3、数111の式を換算する必要がある。
ここで、d、q軸座標におけるd軸に対する各位相角θxのインダクタンスLxは、d軸インダクタンスをLd、q軸インダクタンスをLqとすると、位相角θxに関して、数112、数113の式として表すことができる。これは、各巻線のインダクタンスが、ロータ磁極の向きd軸に対する回転角の位置、すなわち、位相角θxの大きさにより、電気角で180°周期の正弦波的増減を繰り返す特性であることを表す式となっている。
Figure 2008141835
Figure 2008141835
数112の式より換算して、モータ電流(id、iq)、ロータ回転位置θrにおける固定座標の3相インダクタンスLu、Lv、Lwを、図9あるいは図10のデータテーブルを使用して数114〜数116の式のように計算することができる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
同様に、数104、数105、数106の式のLd、Lqに微少電流範囲のインダクタンスLstd、Lstqを代入することにより、微少電流範囲のインダクタンスLstu5、Lstv5、Lstw5として計算することができる。この時、この微少電流範囲のインダクタンスは、図3、数111の式で示した方法で、動作点462の接線であって、直線463の勾配としての各相の微少電流範囲のインダクタンスLstu5、Lstv5、Lstw5を求めたことになる。
次に、685は前記の推定回転位置θsr、推定角速度ωsの推定がどの程度の精度であるかを評価するために誤差評価関数δsを計算する。この誤差評価関数δsの例は、微少電流範囲のインダクタンスLstu5、Lstv5、Lstw5とLstu4、Lstv4、Lstw4とのそれぞれの差の絶対値の式として数117〜数119として表すことができる。
Figure 2008141835
Figure 2008141835
Figure 2008141835
そして、誤差判定手段686で、誤差評価関数δstu、δstv、δstwの各値が誤差許容値δoより小さいかどうかを判定し、誤差許容値δo以下であれば688で推定回転位置θsr、推定角速度ωsがほぼ正しい値であると判定し、689でこの計算処理を終了する。誤差評価関数δstu、δstv、δstwが誤差許容値δo以上であれば、687で誤差評価関数δstu、δstv、δstwの値が小さくなる方向へ、回転位置θsr、角速度ωsの値を補正計算し、再設定する。そして、経路68Aで684へもどり、電圧Vsd、Vsq等の計算を再度行い、この計算処理を誤差評価関数δstu、δstv、δstwが誤差許容値δo以下となるまで繰り返す。このようにして、誤差評価関数δstu、δstv、δstwが誤差許容値δo以下に収束するまで収束計算を行うことにより、推定回転位置θsr、推定角速度ωsを高精度な推定値とすることができる。
なお、667では、誤差評価関数δsの収束が早くなるような工夫、あるいは、計算回数が多くなった場合は終了するなどの例外処理を行っても良い。また、図30に示す計算アルゴリズムは、図28の例と同様に、同一目的で、変形することもできる。また、図28の方法はある程度の回転数で回転している時に精度良く位置を検出することができ、図30の方法はモータが静止あるいは低速回転においても位置検出が容易であるという傾向がある。従って、一つの方法として、低速回転では図30の方法で制御し、ある回転以上の高速回転では図28の方法で制御することもできる。また、制御の切り替えは、中間の回転数では両方式の平均値を取るなど、検出結果の重み付けを徐々に変えていく方法も可能である。また、図28の方法と図30の方法とを重ね合わせたセンサレス位置検出を実現することもできる。
(磁束鎖交数あるいはインダクタンスのデータテーブルを用いた拡張誘起電圧方式によるセンサレスの位置検出/速度検出方法の説明)
次に、センサレス位置検出方法の他の方法を示す。図32は、図29のセンサレス位置検出を含むモータ制御装置の一部を変形し、より詳しく記載したブロックダイアグラムである。モータの座標は、センサレス位置検出のアルゴリズムをより簡素に記載するため、回転座標系であるd、q軸座標と固定座標系であるα、β軸座標、および、U、V、W軸座標を使用している。図32の左上部の入力信号vd、vqは、それぞれ、d軸電圧指令、q軸電圧指令であり、図29のそれらと同じものである。711は、d、q軸信号をα、β軸信号vα、vβに変換する座標変換手段である。712は、α、β軸信号vα、vβをU、V、W軸信号vu、vv、vwに変換する座標変換手段である。713は、3相電流iu、iv=−iu−iw、iwをα、β軸座標の電流iα、iβへ変換する座標変換手段である。714は、ロータの推定回転位置θsrを使用して、d、q軸座標の推定電流isd、isqへ変換する座標変換手段である。
715は、図29のセンサレス検出手段59Sの一例であり、オブザーバを構成し、回転位置推定値θsrと回転角速度推定値ωsを推定計算し、出力する。715の各信号は電圧ベクトルv、電流ベクトルiなどで示しているが、以下に、固定座標系であるα、β軸座標で説明する。717は電圧vで、vα、vβである。破線で囲んで示す716は、モータ591を示している。718は、モータ591の電流を検出したiα、iβである。ここで、モータ591の電圧方程式は、次の数120の式で示される。
Figure 2008141835
ここで、Kはq軸方向の磁束鎖交数である。また、d軸方向の磁束鎖交数を付加することもできる。この数120の式は、固定座標系のα、β軸座標に変換し、さらに変形すると次の数121の式で示される。
Figure 2008141835
数121の式の第1項はモータ回転位置θrを含まないインピーダンス項EIであり、第2項はモータ回転位置θrを含む拡張誘起電圧項EEである。
なお、数121の式の変形については、「拡張誘起電圧モデルに基づくシンクロナスリラクタンスモータのセンサレス制御とそれに適したインダクタンス測定法」、電気学会論文誌D(産業応用部門誌)、2005年、Vol.125、No.1、16〜25ページおよびその参考文献等に記載されている。
今、図32の716はモータで、入力電圧vから前記拡張誘起電圧項EEを加減算機719で減算すると、その出力は(121)式の第1項のインピーダンス項EIの電圧である。71Aが数121の式の第1項の前部に記載されたインピーダンス{(R+pL)I−ω(L−L)J)}であると、その出力は電流iであり、iα、iβであると考えることができる。
ここで、71Bが前記インピーダンス{(R+pL)I−ω(L−L)J)}の逆数であるとすると、71Bの出力は数121の式の第1項のインピーダンス項EIの電圧である。そこで、717が電圧vでvα、vβであり、718が検出された電流iでiα、iβであるとすると、加減算機71Cの出力71Fはモータ716の電圧EEであり、数121の式の第2項のモータ回転位置θrを含む拡張誘起電圧項EEが計算されたことになる。71Dはノイズ等を除去するフイルターである。
71Eは推定回転位置θsr、推定角速度ωsを推定計算する位置検出手段である。フイルター71Dの出力はα、β軸の電圧v1α、v1βであり、数121の式に相当することから次の数122〜数123の式で表すことができる。
Figure 2008141835
Figure 2008141835
モータの推定回転位置θsrは、数122、数123の式より次の数124の式で計算される。
Figure 2008141835
モータの推定角速度ωsは次式で計算されることになり、推定回転位置θsrの変化率から数125の式で計算される。
Figure 2008141835
理論的には、数110〜数125の式で計算されることになるが、モータの磁気特性が非線形であり、磁気飽和もすることから、その誤差に起因する検出誤差が問題となる。本発明では、コンピュータにより非線形有限要素法で計算された図9、図10に示すようなデータテーブルあるいはそのデータから作られた近似式を使用して、数110〜数125の式の計算を行うことにより、高精度なセンサレス位置検出を実現することができる。そして、センサレス位置検出が正確であると、モータ制御装置の位置、および、速度の制御ゲインを高くすることが可能となり、モータの制御応答性が改善される。
また、前記センサレス検出において、パラメータの収束演算などを付加する場合の収束回数も短くすることができる。数110〜数123の式に含まれるモータパラメータは、モータの回転位置θr、電流値id、iqの影響を受けるパラメータが多く有り、モータ制御において同一サンプリングタイミングにおいても、より正確な推定回転位置θsrで数110〜数123の式の計算を行った方が数124、数125の式の計算結果が正確になり、この計算結果が回数と共に収束する関係となることもある。
(電流のフィードバックを固定座標系と回転座標系とで切り替える場合の説明)
次に、モータの位置制御、速度制御、電流制御において、エンコーダあるいはレゾルバなどの回転位置検出器を使用せずに、センサレス位置検出制御によりモータの制御を行い、より安全に、より高信頼にモータを制御する方法について説明する。このような場合の従来の問題点は、センサレス位置検出の位置精度が不正確な状態においても、モータ電流の電流フィードバックの電流制御ループを回転座標系を使用していることに起因する。極端な例としては、センサレス位置検出の回転位置情報θsrの位置誤差θeが電気角で180°の誤差である場合は、各相の電流制御は負帰還ではなく、正帰還になってしまい、電流制御のフィードバックループは発散し、インバータの電流制限の値まで過大な電流が流れてしまうことになる。また、前記回転位置情報θsrの位置検出誤差θeが比較的小さな場合においても、その位置検出誤差θeに相当した電流誤差Δieが発生し、そして、電流制御ループのオーバーシュートなどの不安定さの問題を誘発することもある。本発明の目的は、これらの電流誤差Δieを軽減し、その弊害を低減することである。
図33に本発明の構成例を示す。図33は、図29の制御装置の例の一部を変更した構成となっており、記載されていない部分は図29と同じである。電圧信号発生手段134、電流制御ループの設定手段135についても記載を省略している。721はd軸電流誤差検出手段598の出力に基づいて補償制御しd軸電流制御電圧信号Vdを出力するd軸電圧制御手段であり、722はq軸電流誤差検出手段59Cの出力に基づいて補償制御しq軸電流制御電圧信号Vdを出力するq軸電圧制御手段である。723は固定座標系であるα軸の電圧指令Vα1とVα2とを加算しα軸電圧指令Vαを出力する加算器であり、724は固定座標系であるβ軸の電圧指令Vβ1とVβ2とを加算しβ軸電圧指令Vβを出力する加算器である。726は回転座標のd、q軸信号を固定座標のα、β軸信号に変換する座標変換器であり、α軸の電流指令iα、β軸の電流指令iβを出力する。727はα軸の電流指令iαと計測したα軸の電流iαとを減算するα軸電流誤差検出手段であり、728はα軸の電流指令iβと計測したβ軸の電流iβとを減算するβ軸電流誤差検出手段である。729はα軸電流誤差検出手段727の出力に基づいて補償制御しα軸電流制御電圧信号Vα2を出力するα軸電圧制御手段であり、72Aはβ軸電流誤差検出手段728の出力に基づいて補償制御しβ軸電流制御電圧信号Vβ2を出力するβ軸電圧制御手段である。
図33の制御装置の構成は、図29の制御装置の構成に比較し、固定座標であるα、β軸座標での電流指令と電流負帰還と補償制御およびそれらの出力であるα、β軸電流制御電圧信号Vα2、Vβ2のモータ制御電圧信号への付加、加算とが追加されている。その結果、図33の制御では、d、q軸座標での電流の電流負帰還と、α、β軸座標での電流負帰還とが2重に存在することになる。そして、両方の電流負帰還をそれらの長所短所から特性を分担し並列使用する方法、あるいは、モータの運転状況により前記の両方の電流負帰還を使い分ける方法などがある。
ここで、モータの電流制御特性と制御する座標との関係について考察する。交流同期モータは、その回転と共に、ロータへステータの3相電流あるいは多相電流で起磁力を与え、トルクを発生している。従って、d、q軸座標のような回転座標からその起磁力等の電磁的作用を観ると、直流モータの制御と同様に扱うことができ、制御電圧、制御電流の扱いを単純化できるという特徴がある。その具体的な特長例は、電流の比例、積分補償制御において、積分補償制御を直流モータと同様に簡便に扱うことができる点である。この積分補償制御とは、低周波数領域の補償制御でもあり、電流誤差が極力小さくなるように制御することができという言い方もできる。また、前記したように回転座標系では、センサレス位置検出制御においては、推定回転位置θsrの位置誤差θeの大きさが不定の状態も発生するので、電流制御がモータ回転位置θrの不明確な状態で、正確なモータ回転位置θr情報を必要とする回転軸座標系で扱われることについては、技術的矛盾があり、その取り扱いに注意を要する。一方、固定座標系であるα、β軸座標、あるいは、U、V、W軸座標における電流制御においては、電流の負帰還に回転位置情報を必要としないが、交流電流を扱うことになり、電流制御の負帰還、補償制御において、低周波数領域の補償制御、すなわち、回転周波数より低い周波数成分の補償制御を行うことは困難である。
前記の2組の電流負帰還を並列使用する具体的な方法の例について説明する。回転座標であるd、q軸座標では電流制御、補償制御の内の定周波成分の制御をd、q軸電圧制御手段721、722で行い、座標変換手段711でα、β軸座標の電圧指令信号vα1、vβ1に変換する。固定座標であるα、β軸座標では電流制御、補償制御の内の高周波成分の制御をα、β軸電圧制御手段729、72Aで行い、α、β軸座標の電圧指令信号vα2、vβ2を出力する。そして、両座標の計算結果を加算機723、724で加算してモータ制御電圧指令Vα、Vβとする方法である。このような方法とすることにより、それぞれの座標で優れた面の特性を発揮することができるので、電流ループの即応性、安定性が改善され、かつ、電流誤差に対しては定周波成分の積分制御による電流誤差の低減を図ることができる。また、前記の補償制御には、比例制御、積分制御、微分制御、それらの組み合わせなど種々の方法があり、それらの種々方法で実現可能である。一つの方法として、d、q軸電圧制御手段721、722では積分補償制御だけを行い、α、β軸電圧制御手段729、72Aでは比例制御だけを行う構成に単純化することもできる。
次に、推定回転位置θsrの位置誤差θeが大きい期間が発生する恐れがある場合の、前記の2組の電流負帰還を使用する具体的な方法の例について説明する。まず、特に図示しないが、前記推定回転位置θsrが不定である期間あるいは位置誤差θeが過大である期間を認識する位置誤差量認識手段を設ける。例えば、センサレス位置検出が困難な状況においては、前記位置誤差量認識手段が「位置誤差θeが過大」と認識する。そして、「位置誤差θeが過大」と認識された場合には、α、β軸電圧制御手段729、72Aの出力を零あるいは小さな値とする。同時に、特に図示しないが、固定座標の電流指令iα、iβをα、β軸電流誤差検出手段727、728へ与え、各相の電流を制御する。なお、この時の電流指令値iα、iβの値は、その時のモータ制御システムに対するモータ運転モードなどの指令により決められる。そして、「位置誤差θeが過大では無く、制御可能」と認識される場合には、前記の2組の電流負帰還を並列使用する方法で制御する。また、「位置誤差θeが過大では無く、制御可能」と認識される場合には、d、q軸座標の電流制御だけで電流制御し、α、β軸電圧制御手段729、72Aの出力を零あるいは小さな値とすることもできる。
なお、固定座標での電流制御をα、β軸座標での制御の例について示したが、U、V、W軸座標での電流制御でも可能であり、図20の制御ブロックダイアグラムを変形することも可能である。また、他の固定軸座標で電流制御を行うことも可能である。
以上、本発明に関する種々形態の例について説明したが、本発明を種々変形も可能であり、本発明に含むものである。例えば、相数については3相交流のモータについて説明したが、単相、2相、4相、5相、7相、さらに相数の大きい多相が可能である。小容量の機器においては、コストの観点から部品点数が少ないことが望ましく相数の少ない2相、3相が有利であるが、トルクリップルの観点あるいは大容量機器の場合の1相のパワーデバイスの最大電流制約の点等では相数が多い方が有利なこともある。極数についても限定するものではない。また、モータの構造、分布巻き、短節巻きなどの巻線の形態などの変形が可能である。
また、ロータの種類については、表面磁石型のロータについて多く説明したが、図52〜56に示すようなロータ、さらに、ロータに巻線を持った巻線界磁型ロータ、軸方向端に固定された界磁巻線を持ちギャップを介してロータに磁束作り出すいわゆるクローポール構造ロータなどの種々ロータへの適用も可能である。永久磁石の種類、形状についても限定するものではない。
図57はシンクロナスリラクタンスモータを改良したモータのロータの例であり、スリットS02が複数個設けられ、各スリットに挟まれた細い磁路S03が構成されている。ロータの各磁極へは各巻線S04〜S05、S06〜S07、S08〜S09、S0A〜S0Bが巻回されていて、各巻線の結線方向、接続方向は、図57に示すロータの磁極方向の各磁束を強める方向に巻回されている。そして、その各巻線の結線は直列、並列などが可能であるが、図58に示すように、巻線S0C、S0D、S0E、S0FおよびダイオードS0Gが直列に接続されていて、電流の方向が決められている。この時、これらの巻線に流れる電流は界磁磁束を作る界磁電流となる。ダイオードSOGは、各巻線を複数組に分割してダイオードを介して閉回路とする構成としても良い。このような構成のモータで、ステータ巻線の電流によりロータ側の界磁電流を誘起して、この界磁電流も使いながら制御するモータである。
図59は、図57のロータを2極のモデルとして図示した断面図であり、724はスリットで空隙あるいは磁気抵抗の大きい樹脂など、725は軟磁性体の細い磁路、721、722はロータ内の巻線である。ステータの電流をd、q軸電流にモデル化し、id、iq、iaとして示している。
一般的に、良く使用されている誘導電動機、シンクロナスリラクタンスモータなどを3相インバータと3相巻線で駆動する場合、一つの相の巻線に流れる電流は、「回転と共に界磁エネルギーをモータへ与え、回転と共に界磁エネルギーを電源へ返す」という動作を繰り返している。そして、3相の電流の合成量としては、モータ内部に一定の界磁エネルギーを保って回転させている。この電源とモータを行き来するエネルギーは、良く知られているように、無効電力である。モータの抵抗が零であるような理想的な状態ではエネルギー消費が無く有害ではないが、現実的にはインバータの電流負担となり、インバータ損失、モータの銅損を発生させ、その意味では有害であり、問題である。図57、59のシンクロナスリラクタンスモータを改良したモータは、ステータ巻線に行き来する無効電力を低減し、モータ効率を改善し、インバータ負担を小さくすることが可能なモータである。
図59のステータ電流はベクトル図で書くと図60となり、界磁磁束の励磁電流であるd軸電流idとトルク電流であるq軸電流iqと合成電流iaが通電され、合成電流の電圧はvaとなる。従って、力率はcos(θa)となり、モータ効率、モータを駆動するインバータ効率を改善するためには、ステータ電流と電圧との位相差を小さくすることが有効である。この対策として、モータへ一度供給した界磁エネルギーを回生せず、モータ内部へ保持させることにより電源とモータ間の磁気エネルギーの授受、すなわち無効電流を低減する。具体的には、図61に示す界磁エネルギを示す磁気回路731、ステータ巻線733、ロータ巻線734とダイオード735の等価モデルの関係で示される。ステータ巻線733の電流は図59のd軸電流成分idをイメージしていて、界磁磁束732が生成される磁気回路の磁気エネルギーを供給する。この磁気エネルギーは、電源側へ回生されることはなく、ステータ巻線733の電流が零となるあるいは減少すると、界磁磁束732が減少しないように、そして、ステータ巻線733の電流の減少を補うように、巻線734にダイオード735の方向に電圧が発生し電流が流れる。このような作用により、ステータの励磁電流が途切れた期間も界磁磁束732を一定に保つことができる。そして、このような構成となっているので、ステータ電流の通電方法に種々の工夫を加えることが可能となる。すなわち、インバータとステータ巻線で、例えば、ごく短時間tdに主に励磁電流idを通電して磁気エネルギーを与え、他の時間帯tqは主にトルク電流iqを与える。このようにすると、モータの平均トルクはトルク電流iqに相当するトルクに比較してtq/(td+tq)に減少するが、tdが相対的に小さい値であれば実用的に問題ない。一方、インバータの電流負担、モータのステータ銅損は、電流が合成電流iaではなく、iq以下に低減できるので低減する。結果として、インバータ容量を低減し、モータ効率を向上することが可能となる。単純には、モータ単体のコストは上がる様にも考えられるが、モータ効率が向上することによるモータの小型化の効果があり、モータコストの低減も不可能ではない。
なお、図57、59のモータのさらなる改良策として、このロータ内に永久磁石を配置して力率を改善するなどの効果により、モータ効率、ピークトルク、インバータ駆動効率などを向上することも可能である。永久磁石の配置位置は、限定されるわけではないが、スリット部724の一部あるいは前部に配置でき、磁石の方向はd軸方向、あるいはq軸方向、あるいはその間が可能であり、必要なモータ特性によって選択できる。
図57、59に示すようなモータの制御においては、特に正確なモータの磁気特性が必要であり、図9、図10に示すようなインダクタンス、磁束鎖交数のデータテーブルを使用して効果的な制御を行うことができる。
以上、本発明について説明したが、本発明の特徴、効果として次のように言える。
(a)モータの特性を正確に表すパラメータは、動作点のインダクタンスあるいは動作点の磁束鎖交数である。そして、非線形で磁気飽和を含むモータ特性は、これらのデータテーブルを使用して、全動作領域を正確にモデル化できる。その結果、正確な磁束データを持つことになり、高精度で高応答な制御を実現できる。
(b)コンピュータによる非線形有限要素法などによりモータ特性を計算し、そのデータを使用して正確かつ容易に前記データテーブルを作成できる。データテーブルの値は離散的なデータであるが、 各動作点(id 、iq)の値は内挿計算して求めることができる。データテーブルのサイズは、電流id 、iqを各10ポイントについて、磁束鎖交数Ψd 、Ψqを2バイトデータとして合計400バイト程度のメモリ容量であり、全く問題ない量である。漏れ磁束Ψleak、鉄損等のその他の制御情報も付加できる。モータを制御するマイクロプロセッサの処理時間は、処理内容が単純なので、従来と同等あるいは短縮できると推定される。
(c)永久磁石を内蔵するモータについても同じ方法が適用できる。種々の形式のモータへ適用できる。
(d)回転子のd、q軸の方向が曖昧な電動機においても、座標軸に制約を受けず、正確に表現できる。
(e)また、 電動機の開発段階で、電動機のより精密な制御シミュレーションが可能となり、電動機開発から制御、負荷運転までの連成、連携開発が容易となり、開発期間短縮とコスト低減が期待できる。
(f)前記データテーブルを使用したセンサレスモータ制御では、より精密なモータモデルで制御することができ、モータ制御精度、応答性を改善できる。
なお、本特許の基本的な考え方、各数式は、本発明者が日本電気学会へ投稿した論文「シンクロナスリラクタンスモータの磁束鎖交数を用いたインダクタンス算定法とモデル化の提案」(2006年3月13日付け)に記載した内容であり、同論文は日本電気学会論文誌Dの2007年2月号に掲載予定である。
なお、上記説明では、制御方法として説明したが、次に示すように制御装置として本発明の各技術形態(以下、請求項とも称する)を把握することも当然可能である。
[請求項1]
交流モータの電機子巻線に流れる複数の電流を変数とする関数であって前記交流モータの所定部位の前記巻線と鎖交する磁束量に相関を有する量である磁束鎖交数Ψをそれぞれ表す複数の磁束鎖交数関数の値と前記複数の電流の値との間の関係をデータとして予め保持すると共に前記データに基づいて前記磁束鎖交数Ψに関する情報である磁束鎖交数相関量を求める磁束鎖交数決定部と、
求めた前記磁束鎖交数相関量に基づいて前記交流モータのトルク、速度及び回転位置の少なくとも一つを所定の目標値に制御するべく前記複数の電流を制御する制御部と、
を備えることを特徴とするモータの制御装置。
[請求項2]
請求項1において、
前記交流モータは、永久磁石をもつロータを有するモータの制御装置。
[請求項3]
請求項1又は2において、
前記データは、
前記複数の磁束鎖交数関数の値に所定の定数値を加算乃至乗算した値、及び、前記複数の磁束鎖交数関数の値に対して前記複数の電流を変数とする所定の関数の値を加算乃至乗算した値のどちらかと、前記複数の電流の値との関係を示すモータの制御装置。
[請求項4]
請求項1〜3のいずれかにおいて、
前記データは、
前記複数の電流の値に対応する前記複数の磁束鎖交数関数の値の近似値と、前記複数の電流の値との間の関係を示す近似関数の値と、前記複数の電流の値との間の関係を示すモータの制御装置。
[請求項5]
請求項1〜4において、
前記磁束鎖交数関数は、前記複数の電流に基づいて有限要素法により算出されるモータの制御装置。
[請求項6]
請求項1〜3のいずれかにおいて、
前記磁束鎖交数関数は、前記データから内挿法により算出されるモータの制御装置。
[請求項7]
請求項1〜6のいずれかにおいて、
前記データは、dq軸回転座標系により記載されるモータの制御装置。
[請求項8]
請求項1〜6のいずれかにおいて、
前記データは、直交静止座標系又は三相座標系により記載されるモータの制御装置。
[請求項9]
請求項7において、
前記交流モータは、3相交流モータからなり、
前記磁束鎖交数決定部は、
d、q軸に換算した電流(id、iq)が所定値(id1、iq1)である時の前記交流モータの磁束分布から所定のd軸磁束量φnd及びq軸磁束量φnqを求め、
dq軸電流(id、iq)が所定値(id1、iq1)である時のd軸磁束量φndと、前記巻線のうちほぼd軸位置にあると見なせる巻線部分のターン数Nsdとの積φnd×Nsdを3相の前記巻線ごとに求め、
前記3相の各巻線ごとに求めた前記積φnd×Nsdの和を求めてd軸磁束鎖交数の和Σ(φnd×Nsd)とし、
前記磁束鎖交数の和Σ(φnd×Nsd)に所定の定数値を乗算して前記dq軸電流(id、iq)の所定値(id1、iq1)におけるd軸磁束鎖交数Ψd(id1、iq1)を前記磁束鎖交数関数の値として求め、
dq軸電流(id、iq)が所定値(id1、iq1)である時のq軸磁束量φnqと、前記巻線のうちほぼq軸位置にあると見なせる巻線部分のターン数Nsqとの積φnq×Nsqを前記3相の巻線ごとに求め、
前記3相の各巻線ごとに求めた前記積φnq×Nsqの和を求めてq軸磁束鎖交数の和Σ(φnq×Nsq)とし、
前記磁束鎖交数の和Σ(φnq×Nsq)に所定の定数値を乗算して前記dq軸電流(id、iq)の所定値(id1、iq1)におけるq軸磁束鎖交数Ψq(id1、iq1)を前記磁束鎖交数関数の値として求めるモータの制御装置。
[請求項10]
請求項7において、
前記交流モータは、
短節巻き又は集中巻きされて3相の前記巻線が交叉しない3相交流モータからなり、
前記磁束鎖交数決定部は、
d、q軸に換算した電流(id、iq)が所定値(id1、iq1)である時の前記交流モータの磁束分布から所定のd軸磁束量φnd及びq軸磁束量φnqを求め、
dq軸電流(id、iq)が所定値(id1、iq1)である時のd軸磁束量φndと、前記巻線のうちほぼd軸位置にあると見なせる巻線部分のターン数Nsdとの積φnd×Nsdを前記3相の各巻線ごとに求め、
前記3相の各巻線ごとに求めた前記積φnd×Nsdの和を求めてd軸磁束鎖交数の和Σ(φnd×Nsd)とし、
前記磁束鎖交数の和Σ(φnd×Nsd)に所定の定数値を乗算して前記dq軸電流(id、iq)の所定値(id1、iq1)におけるd軸磁束鎖交数Ψd(id1、iq1)を前記磁束鎖交数関数の値として求め、
dq軸電流(id、iq)が所定値(id1、iq1)である時のq軸磁束量φnqと、前記巻線のうちほぼq軸位置にあると見なせる巻線部分のターン数Nsqとの積φnq×Nsqを前記3相の各巻線ごとに求め、
前記3相の各巻線ごとに求めた前記積φnq×Nsqの和を求めてq軸磁束鎖交数の和Σ(φnq×Nsq)とし、
前記磁束鎖交数の和Σ(φnq×Nsq)に所定の定数値を乗算して前記dq軸電流(id、iq)の所定値(id1、iq1)におけるq軸磁束鎖交数Ψq(id1、iq1)を前記磁束鎖交数関数の値として求めるモータの制御装置。
[請求項11]
請求項7において、
前記交流モータは、
同相のステータ磁極を6極以上有するとともに前記巻線がティースに巻装されたループ状の巻線からなる3相交流モータからなり、
前記磁束鎖交数決定部は、
d、q軸に換算した電流(id、iq)が所定値(id1、iq1)である時の前記交流モータの磁束分布から所定のd軸磁束量φnd及びq軸磁束量φnqを求め、
dq軸電流(id、iq)が所定値(id1、iq1)である時のd軸磁束量φndと、前記ティースのうちほぼd軸位置にあると見なせる前記ティースに巻装された前記巻線のターン数Nsdとの積φnd×Nsdを前記3相の各巻線ごとに求め、
前記3相の各巻線ごとに求めた前記積φnd×Nsdの和を求めてd軸磁束鎖交数の和Σ(φnd×Nsd)とし、
前記磁束鎖交数の和Σ(φnd×Nsd)に所定の定数値を乗算して前記dq軸電流(id、iq)の所定値(id1、iq1)におけるd軸磁束鎖交数Ψd(id1、iq1)を前記磁束鎖交数関数の値として求め、
dq軸電流(id、iq)が所定値(id1、iq1)である時のq軸磁束量φnqと、前記ティースのうちほぼd軸位置にあると見なせる前記ティースに巻装された前記巻線のターン数Nsqとの積φnq×Nsqを前記3相の各巻線ごとに求め、
前記3相の各巻線ごとに求めた前記積φnq×Nsqの和を求めてq軸磁束鎖交数の和Σ(φnq×Nsq)とし、
前記磁束鎖交数の和Σ(φnq×Nsq)に所定の定数値を乗算して前記dq軸電流(id、iq)の所定値(id1、iq1)におけるq軸磁束鎖交数Ψq(id1、iq1)を前記磁束鎖交数関数の値として求めるモータの制御装置。
[請求項12]
請求項9〜11のいずれかにおいて、
前記磁束鎖交数決定部は、
前記d軸磁束鎖交数Ψd(id1、iq1)とd軸電流id1との比であるd軸の磁束鎖交数比例係数LLd(id1、iq1)と、前記q軸磁束鎖交数Ψq(id1、iq1)とq軸電流iq1との比であるq軸の磁束鎖交数比例係数LLq(id1、iq1)とを前記データとして保持するとともに、前記d軸の磁束鎖交数比例係数LLd(id1、iq1)及びq軸の磁束鎖交数比例係数LLq(id1、iq1)に基づいて前記磁束鎖交数関数を導出するモータの制御装置。
[請求項13]
請求項1〜6のいずれかにおいて、
前記制御部は、
前記磁束鎖交数Ψに基づいて作成したトルクと前記電流との関係を表すデータテーブルあるいは関数を用いて前記電流の指令値を決定するモータの制御装置。
[請求項14]
請求項1〜6のいずれかにおいて、
前記制御部は、
前記磁束鎖交数Ψに基づいて作成したトルクと回転数と前記電流との関係を表すデータテーブルあるいは関数を用いて前記電流の指令値を決定するモータの制御装置。
[請求項15]
請求項1〜6のいずれかにおいて、
前記制御部は、
前記磁束鎖交数Ψに基づいて、前記トルク、前記巻線の電圧及び前記磁束鎖交数Ψの電流変化率の少なくとも一つを算出するモータの制御装置。
[請求項16]
請求項1〜6、14及び16のうちのいずれかにおいて、
前記制御部は、
複数の制御モードを有するとともに、前記磁束鎖交数Ψに基づいて作成したトルクと前記電流との関係を表すデータテーブルあるいは関数を前記各制御モードごとにそれぞれ有し、
選択された前記制御モードに応じて選択した前記データテーブルあるいは関数を用いて前記電流の指令値を決定乃至内挿法にて算出するモータの制御装置。
[請求項17]
請求項1〜6のいずれかにおいて、
前記制御部は、
前記磁束鎖交数Ψ又は前記磁束鎖交数Ψに基づいて得た電流値又は前記巻線のインダクタンス値又はそれらの近似関数値を用いて前記モータへ出力するべき電圧の指令値を求めるモータの制御装置。
[請求項18]
請求項1〜6のいずれかにおいて、
前記制御部は、
d軸電流制御ループ及びq軸電流制御ループを有するフィードバック制御系を有するとともに前記d軸電流制御ループのループゲインGd及び前記q軸電流制御ループのループゲインGqを制御するに際して、電流動作点(id、iq)の近傍における前記ループゲインGdをd軸電流idに対するd軸磁束鎖交数Ψdの微少変化率ΔΨd/Δidに比例する値に設定し、電流動作点(id、iq)の近傍における前記ループゲインGqをq軸電流iqに対するq軸磁束鎖交数Ψqの微少変化率ΔΨq/Δiqに比例する値に設定する制御を行うモータの制御装置。
[請求項19]
請求項1〜6のいずれかにおいて、
d軸電流制御ループ及びq軸電流制御ループを有するフィードバック制御系、もしくはd軸電圧とq軸電圧とをフィードフォワード制御するフィードフォワード制御系を有するとともに、
前記d軸電流制御ループのループゲインGd及び前記q軸電流制御ループのループゲインGqもしくは前記フィードフォワード制御系のd軸電圧フィードフォワード量FFd及びq軸電圧フィードフォワード量を、前記複数の電流としてのd、q軸電流(id、iq)を用いた前記データに基づいて演算するモータの制御装置。
[請求項20]
請求項1〜6のいずれかにおいて、
前記磁束鎖交数決定部又は前記制御部は、
d、q軸の漏れ磁束成分として定義されるd軸漏れ磁束鎖交数Ψkd及びq軸漏れ磁束鎖交数Ψkq、あるいは、d軸漏れインダクタンスLkd及びq軸漏れインダクタンスLkqを求めて前記データの一部として用いるモータの制御装置。
[請求項21]
請求項1〜6のいずれかにおいて、
前記磁束鎖交数決定部又は前記制御部は、
電流動作点(id、iq)及び回転角速度ωに基づいて求めた鉄損に相関を有する電気量を前記データの一部として用いるモータの制御装置。
[請求項22]
請求項1〜6のいずれかにおいて、
前記制御部は、
d軸電流制御ループ及びq軸電流制御ループを有するフィードバック制御系を有するとともに前記フィードバック制御の遅れ時間Δtを補償するための電流制御を行うモータの制御装置。
[請求項23]
請求項1乃至6のいずれか記載のモータの制御装置において、
前記複数の電流の各値と前記磁束鎖交数相関量の各値とのペアを記載するデータテーブルを前記データとして有し、
前記制御部は、
前記データテーブルを用いてロータ位置又はその角速度の演算を行うモータの制御装置。
[請求項24]
請求項1〜6のいずれかにおいて、
予め求めた前記交流モータの磁束鎖交数Ψと前記複数の電流との関係を示すデータテーブルあるいは前記巻線のインダクタンスのデータテーブルからなり、前記交流モータの電圧vと電流iとの関係を示すデータDataSを前記データとして有し、
前記制御部は、
前記データに基づいて求めた前記交流モータの回転位置の推定値θsrに基づいて決定した電圧vと電流iとにより定義される電力を前記交流モータに給電し、
所定の前記回転位置の推定値θsr1と前記データDataSと外部から入力された電流指令i*とにより前記モータの電圧vsを計算し、
前記電流ipと前記電流指令i*との差、及び、前記電圧vpと前記モータの電圧vsとの差がそれぞれの小さくなる前記推定値θsr1の値を前記交流モータの回転位置θsrの推定値として最終的に決定するモータの制御装置。
[請求項25]
請求項1〜6のいずれかにおいて、
予め求めた前記交流モータの磁束鎖交数Ψと前記複数の電流との関係を示すデータテーブルあるいは前記巻線のインダクタンスのデータテーブルからなり、前記交流モータの電圧vと電流iとの関係を示すデータDataSを前記データとして有し、
前記制御部は、
前記データに基づいて求めた前記交流モータの回転位置の推定値θsrに基づいて決定した電流ipk、電圧vpで定義される電力を前記巻線に給電し、
前記電流ipkは、所定の高周波電流成分ikを含み、
計測時点における前記電流ipk中の前記高周波電流成分ikの大きさと前記電圧vpの高周波成分vkの大きさとから前記高周波電流成分ikの所定の微少電流範囲内における前記巻線のインダクタンスLst4を求め、
補正した所定の前記回転位置の推定値θsr1と電流指令i*と前記データDataSとにより前記電流指令i*の近傍の微少電流範囲における前記巻線のインダクタンスLst5を計算し、
前記インダクタンスLst4とLst5との差が小さくなる前記回転位置θsr1の値を前記交流モータの回転位置の推定値として最終決定することを特徴とするモータの制御装置。
[請求項26]
請求項1〜6のいずれかにおいて、
予め求めた前記交流モータの磁束鎖交数Ψと前記複数の電流との関係を示すデータテーブルあるいは前記巻線のインダクタンスのデータテーブルからなり、前記交流モータの電圧vと電流iとの関係を示すデータDataSを前記データとして有し、
前記制御部は、
前記電圧vを、ロータの回転位置θrに関係する電圧である電圧項EEと、前記回転位置θrに対して依存性が少ない電圧である電圧項EIとの和とし、
前記複数の電流の計測値と前記データから得た前記巻線のインピーダンス値とにより前記電圧項EIに相当する電圧EIsを求め、
前記電圧vから前記電圧EIsを差し引くことにより前記電圧項EEに相当する電圧EEsを求め、
前記電圧EEsを使用して前記交流モータの回転位置推定値θsrをもとめるモータの制御装置。
[請求項27]
請求項1〜6のいずれかにおいて、
前記制御部は、
前記交流モータの電圧、電流に基づいて前記交流モータの回転位置をセンサレスに推定するセンサレス位置検出手段と、
回転座標系で表示された信号を用いて前記交流モータの電流制御を行う回転座標制御手段と、
静止直交座標系乃至三相静止座標系で表示された信号を用いて前記交流モータの電流を制御する固定座標制御手段と
を備えるモータの制御装置。
シンクロナスリラクタンスモータの横断面図である。 インダクタンスLd、Lqの特性である。 モータの磁束鎖交数Ψの例である。 モータ制御のブロックダイアグラムである。 6極のモータを2極のモータへモデル化した図である。 3相巻線とロータ磁極との配置関係を示す図である。 d、q軸に換算した巻線とロータ磁極との配置関係を示す図である。 モータ巻線と磁束φとベクトルポテンシャルAとの関係を示す図である。 d、q軸電流(id、iq)と各巻線のインダクタンスLd(id、iq)、Lq(id、iq)との関係を示すデータテーブルである。 d、q軸電流(id、iq)と各巻線の磁束鎖交数Ψd(id、iq)、Ψq(id、iq)との関係を示すデータテーブルである。 図1のモータの電流とトルクの関係を示す図である。 図1のモータの、特に定出力領域の、電流とトルクの関係を示す図である。 シンクロナスリラクタンスモータのロータ内部へ永久磁石を付加したモータの横断面図である。 図13のモータの電流とトルクの関係を示す図である。 図5のモータモデルのd、q軸座標をCWの方向へ30°傾けたモータモデルの図である。 図15のモータモデルの電流とトルクの関係を示す図である。 図15のモータモデルにおける電流とインダクタンスLd、Lqを示す図である。 モータおよびその入出力を模式的に示した図である。 シンクロナスリラクタンスのベクトル図である。 モータを定トルクモードで制御する時の、トルクと各電流id、iqの関係を示す図である。 モータを定出力モードで制御する時の、回転数とトルクと電流の関係を示すデータテーブルである。 図1のモータを定出力モードで制御する時の回転数とトルクと電圧の関係を示す図である。 3相交流、集中巻きの構成のモータモデルにおいて、d軸インダクタンスLdを求める方法を示す図である。 3相交流、集中巻きの構成のモータモデルにおいて、q軸インダクタンスLqを求める方法を示す図である。 3相交流のモータで、電気角360°に6スロットが配置されたモータにおいて、d軸インダクタンスLdを求める方法を示す図である。 3相交流のモータで、電気角360°に6スロットが配置されたモータにおいて、q軸インダクタンスLqを求める方法を示す図である。 各巻線の有効な磁束鎖交数成分を表すインダクタンスと漏れ磁束成分を表すインダクタンスLkd、Lkqを示す図である。 モータのセンサレス位置検出を行うフローチャートである。 モータのセンサレス位置検出、モータの速度制御を行うブロックダイアグラムである。 モータのセンサレス位置検出を行うフローチャートである。 U、V、W軸座標のU相巻線、V相巻線、W相巻線、および、α、β軸座標のα相巻線、β相巻線を示す図である。 拡張誘起電圧方式のセンサレス位置検出を示すブロックダイアグラムである。 センサレス位置検出のモータ制御において、電流のフィードバックを安定化させるモータ制御のブロックダイアグラムである。 従来の集中巻きブラシレスモータの縦断面図である。 図34のステータの歯のエアギャップ部に面した形状を直線状に展開し、各相巻線の配置関係を付記した展開図である。 図34のブラシレスモータの断面AA−AAを示す横断面図である。 3相交流、4極で、ロータにスリットを持つシンクロナスリラクタンスモータの横断面図である。 3相インバータと3相交流モータの各巻線との接続を示す図である。 従来のモータ制御装置のブロックダイアグラムである。 ループ状の巻線を持つ3相モータ縦断面図である。 ロータの円周方向展開図である。 図40のモータのAA−AA断面図、AB−AB断面図、AC−AC断面図である。 ロータ側から見たステータの内周側形状の展開図である。 巻線の概略的な形状を示す図である。 各相の巻線の円周方向展開図を示す図である。 2つの巻線を単一の巻線にまとめた変形例を示す図である。 エアギャップ面側から見た各相ステータ磁極の展開図に等価的な各相電流巻線を書き加えた図である。 本実施形態のブラシレスモータの電流、電圧、出力トルクのベクトル図である。 ステータ磁極の変形例を示す円周方向展開図である。 ステータ磁極の他の変形例を示す円周方向展開図である。 ステータ磁極の他の変形例を示す円周方向展開図である。 ロータの横断面図である。 ロータの横断面図である。 ロータの横断面図である。 ロータの横断面図である。 ロータの横断面図である。 スリットを持つロータへ巻線とダイオードを付加したロータである。 ロータの各巻線とダイオードとの接続関係を示す図である。 図57のロータを2極のロータにモデル化した図である。 電流と電圧とのベクトル図である。 ステータのd軸巻線、d軸電流、界磁の磁気回路、ロータの巻線、ロータのダイオードをモデル化して示した図である。

Claims (27)

  1. それぞれモータの複数の電流(ia、ib・・・)を変数とする関数であって磁束鎖交数(Ψ)にそれぞれ相当する複数の磁束鎖交数(Ψ1(ia、ib・・・)、Ψ2(ia、ib・・・))の各値と前記各電流値との関係をデータテーブル等のデータとして備え、前記データを使用して前記モータの電流動作点(ia、ib・・・)近傍における各巻線の電圧あるいは巻線インピーダンス等を求め、モータのトルク、速度、位置等を制御することを特徴とするモータの制御方法。
  2. 請求項1において、
    制御対象とするモータが永久磁石を備えるモータである
    ことを特徴とするモータの制御方法。
  3. 請求項1、2において、
    前記データテーブルの各値は、各磁束鎖交数Ψ1(ia、ib・・・)、Ψ2(ia、ib、・・・)、の各電流値の時の値へ定数を加算した値、定数を乗算した値、動作点の電流値に関する関数を加算した値、あるいは、動作点の電流値に関する関数を乗算した値等のように各磁束鎖交数を変換した値であることを特徴とするモータの制御方法。
  4. 請求項1〜3において、
    モータの各電流値ia、ib・・・における磁束鎖交数を、変数がia、ib・・・である近似関数である関数FΨ1(ia、ib・・・)、FΨ2(ia、ib・・・)・・・として備えることを特徴とするモータの制御方法。
  5. 請求項1〜4において、
    前記磁束鎖交数Ψの値あるいは磁束鎖交数を変換した値は有限要素法解析により計算した値より得られたものであることを特徴とするモータの制御方法。
  6. 請求項1〜3、5において、
    各電流の動作点(ia、ib・・・)の値に対する磁束鎖交数Ψ1(ia、ib・・・)、Ψ2(ia、ib・・・)・・・を前記データテーブルから内挿計算して求めることを特徴とするモータの制御方法。
  7. 請求項1〜6において、
    d、q軸座標系で表現されたデータテーブルであることを特徴とするモータの制御方法。
  8. 請求項1〜7において、
    固定軸の3相の座標系で表現されたデータテーブルであることを特徴とするモータの制御方法。
  9. 請求項7において、
    3相交流モータであって、
    d、q軸に換算した電流が(id1、iq1)の時のモータ内部の磁束の分布を求め、
    d軸方向へ向いた磁束成分であって、各スロットの巻線へ鎖交する磁束φndと各スロットに巻回された巻線数Nsの積φnd×Nsを各巻線ごとに求め、
    それらの磁束鎖交数の和Σ(φnd×Ns)を求め、
    前記磁束鎖交数の和Σ(φnd×Ns)の1/1.4142倍をその電流(id1、iq1)の時のd軸磁束鎖交数Ψd(id1、iq1)=Σ(φnd×Ns)/1.4142とし、
    同様に、q軸方向についても求め、
    q軸方向へ向いた磁束成分であって、各スロットの巻線へ鎖交する磁束φnqと各スロットに巻回された巻線数Nsの積φnq×Nsを各巻線ごとに求め、
    それらの磁束鎖交数の和Σ(φnq×Ns)を求め、
    前記磁束鎖交数の和Σ(φnq×Ns)の1/1.4142倍をその電流(id1、iq1)の時のq軸磁束鎖交数Ψq(id1、iq1)=Σ(φnq×Ns)/1.4142とする
    ことを特徴とするモータの制御方法。
  10. 請求項7において、
    3相交流の短節巻き,集中巻きで,3相の各巻線が交叉しない構造のモータの制御方法であって,
    ロータのq軸をU相ステータ磁極の方向の回転位置へ固定し、3相電流のd,q軸に換算した電流が(id1,iq1)の時のモータ内部の磁束の分布を求め,
    d軸方向へ向いた磁束成分であって,各歯に巻回された巻線へ鎖交する磁束φndと各巻線数Nsの積φnd×Nsを各巻線ごとに求め,
    それらの磁束鎖交数の和Σ(φnd×Ns)を求め,
    前記磁束鎖交数の和Σ(φnd×Ns)の1/1.4142倍をその電流(id1,iq1)の時のd軸磁束鎖交数Ψd(id1,iq1)=Σ(φnd×Ns)/1.4142とし,
    同様に,q軸方向についても,ロータのd軸をU相ステータ磁極の方向の回転位置へ固定し,3相電流のd,q軸に換算した電流が(id1,iq1)の時のモータ内部の磁束の分布を求め,
    q軸方向へ向いた磁束成分であって,各歯に巻回された巻線へ鎖交する磁束φnqと各巻線数Nsの積φnq×Nsを各巻線ごとに求め,
    それらの磁束鎖交数の和Σ(φnq×Ns)を求め,
    前記磁束鎖交数の和Σ(φnq×Ns)の1/1.4142倍をその電流(id1,iq1)の時のq軸磁束鎖交数Ψq(id1,iq1)=Σ(φnq×Ns)/1.4142とすることを特徴とするモータの制御方法。
  11. 請求項7において、
    ループ状巻線を持ち,同一相のステータ磁極を同一円周上に持つ6極以上の3相モータであって,
    ロータのq軸をU相ステータ磁極の方向の回転位置へ固定し,3相電流のd,q軸に換算した電流が(id1,iq1)の時のモータ内部の磁束の分布を求め,
    d軸方向へ向いた磁束成分であって,各歯に巻回された巻線へ鎖交する磁束φndと各巻線数Nsの積φnd×Nsを各巻線ごとに求め,
    それらの磁束鎖交数の和Σ(φnd×Ns)を求め,
    前記磁束鎖交数の和Σ(φnd×Ns)の1/1.4142倍をその電流(id1,iq1)の時のd軸磁束鎖交数Ψd(id1,iq1)=Σ(φnd×Ns)/1.4142とし,
    同様に,q軸方向についても,ロータのd軸をU相ステータ磁極の方向の回転位置へ固定し,3相電流のd,q軸に換算した電流が(id1,iq1)の時のモータ内部の磁束の分布を求め,
    q軸方向へ向いた磁束成分であって,各歯に巻回された巻線へ鎖交する磁束φnqと各巻線数Nsの積φnq×Nsを各巻線ごとに求め,
    それらの磁束鎖交数の和Σ(φnq×Ns)を求め,
    前記磁束鎖交数の和Σ(φnq×Ns)の1/1.4142倍をその電流(id1,iq1)の時のq軸磁束鎖交数Ψq(id1,iq1)=Σ(φnq×Ns)/1.4142とする
    ことを特徴とするモータの制御方法。
  12. 請求項9、10、11において、
    3相交流モータであって、
    d、q軸に換算した電流が(id1、iq1)の時のモータ内部の磁束の分布を求め、
    d軸方向へ向いた各巻線へ鎖交する磁束φndと各巻線数Nsの積φnd×Nsを各巻線ごとに求め、
    それらの磁束鎖交数の和Σ(φnd×Ns)を求め、
    前記磁束鎖交数の和Σ(φnd×Ns)の1/1.4142倍をその電流(id1、iq1)の時のd軸磁束鎖交数Ψd(id1、iq1)=Σ(φnd×Ns)/1.4142とし、
    前記d軸磁束鎖交数Ψd(id1、iq1)とd軸電流id1の比をd軸の磁束鎖交数比例係数LLd(id1、iq1)とし、
    同様に、q軸方向についても求め、
    q軸方向へ向いた各巻線へ鎖交する磁束φnqと各巻線数Nsの積φnq×Nsを各巻線ごとに求め、
    それらの磁束鎖交数の和Σ(φnq×Ns)を求め、
    前記磁束鎖交数の和Σ(φnq×Ns)の1/1.4142倍をその電流(id1、iq1)の時のq軸磁束鎖交数Ψq(id1、iq1)=Σ(φnq×Ns)/1.4142とし、
    前記q軸磁束鎖交数Ψq(id1、iq1)とq軸電流iq1の比をq軸の磁束鎖交数比例係数LLq(id1、iq1)とする
    ことを特徴とするモータの制御方法。
  13. 請求項1〜6において、
    前記磁束鎖交数Ψのデータテーブルより作成されたトルクとモータ電流の関係を表すデータテーブルあるいは関数を備えることを特徴とするモータの制御方法。
  14. 請求項1〜6において、
    前記磁束鎖交数Ψのデータテーブルより作成されたトルクと回転数とモータ電流の関係を表すデータテーブルあるいは関数
    を備えることを特徴とするモータの制御方法。
  15. 請求項1〜6において、
    トルク、モータ電圧、巻線の磁束鎖交数の電流変化率等を前記磁束鎖交数Ψを基礎データとして計算あるいは制御する
    ことを特徴とするモータの制御方法。
  16. 請求項1〜6、14、16において、
    モータの電流指令値を求めるための複数組の電流情報から、モータの制御モード指令に対応する一組を選択して使用するかあるいは前記の複数の電流情報を用いて内挿計算することによりモータの電流指令値を決定することを特徴とするモータの制御方法。
  17. 請求項1〜6において、
    モータの電圧信号を、電流の値に応じた非線形なインダクタンスデータ、あるいは、電流の値に応じた磁束鎖交数データ、あるいは、電流に関するそれらの近似関数を使用して求めることを特徴とするモータの制御方法。
  18. 請求項1〜6において、
    d軸電流制御ループのループゲインをGd、q軸電流制御ループのループゲインGq、電流動作点(id、iq)の近傍におけるd軸磁束鎖交数Ψdのd軸電流idに関する微少変化率をΔΨd/Δid、電流動作点(id、iq)の近傍におけるq軸磁束鎖交数Ψqのq軸電流iqに関する微少変化率をΔΨq/Δiqとした場合に、前記ゲインGd、Gqをそれぞれ磁束鎖交数の電流変化率ΔΨd/Δid 、ΔΨq/Δiqに比例した値とすることを特徴とするモータの制御方法。
  19. 請求項1〜6において、
    d、q軸の電圧のフィードフォワード量FFd、FFq、あるいは、d、q軸電流制御ループのループゲイン定数Gd、Gqをd、q軸電流をid、iqをパラメータとして作成したデータテーブル
    を備えることを特徴とするモータの制御方法。
  20. 請求項1〜6において、
    d、q軸の漏れ磁束成分であるd軸漏れ磁束鎖交数Ψkd、q軸漏れ磁束鎖交数Ψkq、あるいは、d軸漏れインダクタンスLkd、q軸漏れインダクタンスLkqを付加して制御することを特徴とするモータの制御方法。
  21. 請求項1〜6において、
    電流の動作点(id、iq)とその回転角速度ωとに対応した鉄損成分を含めて制御することを特徴とするモータの制御方法。
  22. 請求項1〜6において、
    制御の遅れ時間Δtを補償した電流制御を行うことを特徴とするモータの制御方法。
  23. 請求項1〜6において、
    モータの電流の値ごとの磁束鎖交数情報を記録した前記データテーブル、あるいは、同様のインダクタンスのデータテーブル、あるいは、これらを近似関数化した情報を使用するセンサレス位置検出、センサレス角速度検出を行うことを特徴とするモータの制御方法。
  24. 請求項1〜6において、
    モータの電圧vと電流iの関係を、あらかじめコンピュータ計算などにより得られたモータの磁束鎖交数のデータテーブルあるいはインダクタンスのデータテーブル等のデータDataSとして備え、
    モータの回転位置の推定値θsrに基づいてモータへ電流ip、電圧vpを印加し、
    モータの回転位置の推定値θsr1と前記データDataSと電流指令iよりモータ電圧vsを計算し、
    前記モータ電流ip、電圧vpと前記モータ電流指令i、モータ電圧vsとのそれぞれの差が小さくなるような前記推定値θsr1をモータの回転位置推定値とすることを特徴とするモータの制御方法。
  25. 請求項1〜6において、
    モータの電圧vと電流iの関係を、あらかじめコンピュータ計算などにより得られたモータの磁束鎖交数のデータテーブルあるいはインダクタンスのデータテーブル等のデータDataSとして備え、
    モータの回転位置の推定値θsrに基づいて高周波電流成分ikを重畳した電流ipk、電圧vpを印加し、
    計測タイミングでの電流ipkの高周波成分ikと電圧vpの高周波成分vkとからその高周波電流成分の微少電流範囲におけるインダクタンスLst4を計測、計算し、
    モータの回転位置の補正した推定値θsr1と電流指令iより、前記データDataSを使用して、モータ電流指令iの近傍の微少電流範囲におけるインダクタンスLst5を計算し、
    前記インダクタンスLst4とLst5との差が小さくなるような回転位置θsr1をモータの回転位置の推定値とすることを特徴とするモータの制御方法。
  26. 請求項1〜6において、
    モータの電圧vと電流iの関係を、あらかじめコンピュータ計算などにより得られたモータの磁束鎖交数のデータテーブルあるいはインダクタンスのデータテーブル等のデータDataSとして備え、
    モータの電圧方程式等をモータの回転位置θrに関係する電圧項EEとモータの回転位置θrに対する依存性に低い電圧項EIに分離し、
    モータの計測した電流値とモータのインピーダンスより前記電圧項EIに相当する電圧EIsを求め、
    モータの電圧vから前記電圧EIsを差し引くことにより前記電圧項EEに相当する電圧EEsを求め、
    この電圧EEsを使用してモータの回転位置推定値θsrをもとめることを特徴とするモータの制御方法。
  27. 請求項1〜6において、
    モータの回転位置の検出を行うエンコーダ等の位置検出手段の代わりに、モータの電圧、電流等の情報からモータの回転位置を検出するセンサレス位置検出手段と、
    d相、q相等の回転座標で制御するモータ電流の回転座標制御手段と、
    α相、q相、あるいは、U相、V相、W相等の固定座標で制御するモータ電流の固定座標制御手段と
    を備えることを特徴とするモータの制御装置。
JP2006324077A 2006-11-30 2006-11-30 モータの制御方法及びそれを利用するモータ制御装置 Expired - Fee Related JP5167631B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2006324077A JP5167631B2 (ja) 2006-11-30 2006-11-30 モータの制御方法及びそれを利用するモータ制御装置
US11/987,365 US7772790B2 (en) 2006-11-30 2007-11-29 System and method for controlling motor using parameter associated with magnetic flux
DE102007057499.3A DE102007057499B4 (de) 2006-11-30 2007-11-29 System und Verfahren zum Steuern eines Motors unter Verwendung eines mit einem magnetischen Fluss zusammenhängenden Parameters
CN2007101946877A CN101192808B (zh) 2006-11-30 2007-11-30 使用磁通量相关参数来控制电机的系统和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006324077A JP5167631B2 (ja) 2006-11-30 2006-11-30 モータの制御方法及びそれを利用するモータ制御装置

Publications (2)

Publication Number Publication Date
JP2008141835A true JP2008141835A (ja) 2008-06-19
JP5167631B2 JP5167631B2 (ja) 2013-03-21

Family

ID=39339162

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006324077A Expired - Fee Related JP5167631B2 (ja) 2006-11-30 2006-11-30 モータの制御方法及びそれを利用するモータ制御装置

Country Status (4)

Country Link
US (1) US7772790B2 (ja)
JP (1) JP5167631B2 (ja)
CN (1) CN101192808B (ja)
DE (1) DE102007057499B4 (ja)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010041839A (ja) * 2008-08-06 2010-02-18 Fuji Electric Systems Co Ltd 永久磁石形同期電動機の制御装置
WO2010029879A1 (ja) * 2008-09-11 2010-03-18 本田技研工業株式会社 電動パワーステアリング装置
EP2290806A2 (en) 2009-08-27 2011-03-02 Hitachi, Ltd. Control apparatus for synchronous motor
CN102013870A (zh) * 2010-11-11 2011-04-13 江苏大学 五自由度无轴承同步磁阻电机逆系统解耦控制器
DE112010001465T5 (de) 2009-03-30 2012-05-16 Hitachi, Ltd. Wechselstrommotor-Steuervorrichtung und Wechselstrommotor-Treibersystem
KR101328192B1 (ko) * 2009-06-23 2013-11-08 가트너-일렉트로닉-디자인 게엠베하 자석과 함께 고정자 및 회전자를 구비한 브러시리스 직류 모터(bldc) 내의 부하율을 보상하기 위한 장치
WO2015132991A1 (ja) * 2014-03-05 2015-09-11 三菱電機株式会社 シンクロナスリラクタンスモータ
JP2015192495A (ja) * 2014-03-27 2015-11-02 株式会社豊田中央研究所 情報処理装置、情報記憶装置、及び回転電機の制御装置
JP2016005427A (ja) * 2014-06-16 2016-01-12 現代自動車株式会社Hyundaimotor Company 電動機のセンサレス制御方法およびシステム
JP2016140139A (ja) * 2015-01-26 2016-08-04 株式会社デンソー 回転電機の制御装置
KR20170008381A (ko) * 2015-07-13 2017-01-24 현대모비스 주식회사 모터 속도 측정 장치 및 방법
JP2017085879A (ja) * 2015-10-13 2017-05-18 シュネーデル、トウシバ、インベーター、ヨーロッパ、ソシエテ、パル、アクション、セプリフエSchneider Toshiba Inverter Europe Sas 可変リラクタンス同期電気モータのインダクタンス値を識別するための制御プロセス
JP2018085851A (ja) * 2016-11-24 2018-05-31 日産自動車株式会社 可変磁束モータの電流制御方法、及び電流制御装置
JP2018085852A (ja) * 2016-11-24 2018-05-31 日産自動車株式会社 可変磁束型回転電機の電流制御方法、及び電流制御装置
KR20190016377A (ko) * 2017-08-08 2019-02-18 엘지전자 주식회사 모터 구동 장치 및 이를 구비하는 공기조화기
JP2019068601A (ja) * 2017-09-29 2019-04-25 日産自動車株式会社 モータの制御方法、及び、モータの制御装置。
JP2019068600A (ja) * 2017-09-29 2019-04-25 日産自動車株式会社 モータの制御方法、及び、モータの制御装置
CN110445439A (zh) * 2018-05-04 2019-11-12 美的集团股份有限公司 永磁同步电机的控制方法和装置
CN110460284A (zh) * 2018-05-04 2019-11-15 美的集团股份有限公司 食物料理机、电机驱动系统、感应电机的控制方法和装置
CN110829919A (zh) * 2019-10-29 2020-02-21 西安石油大学 一种连续波泥浆脉冲器的控制系统及其开发方法
WO2020075620A1 (ja) * 2018-10-10 2020-04-16 日本電産株式会社 モータ制御装置およびパワーステアリング装置
US11296633B2 (en) 2018-01-12 2022-04-05 Mitsubishi Electric Corporation Rotary machine control device
KR102500019B1 (ko) * 2022-12-29 2023-02-16 주식회사 에스엠전자 2승저감부하에 연결된 모터의 운영전력효율이 개선된 제어관제시스템

Families Citing this family (125)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0724837D0 (en) * 2007-12-20 2008-01-30 Edwards Ltd vacuum pump
US8115437B2 (en) * 2008-02-22 2012-02-14 Rockwell Automation Technologies, Inc. Torque limit of PM motors for field-weakening region operation
US7847501B2 (en) * 2008-03-07 2010-12-07 Tesla Motors Varying flux versus torque for maximum efficiency
EP2151918A1 (en) * 2008-08-07 2010-02-10 Bombardier Transportation GmbH Operating a synchronous motor having a permanent magnet rotor
CN102171923B (zh) * 2008-08-28 2014-03-12 Thk株式会社 线性同步电机控制装置
US8228013B2 (en) * 2008-09-10 2012-07-24 Rockwell Automation Technologies, Inc. Method and apparatus for automatically identifying electrical parameters in a sensor-less PMSM
EP2343799A4 (en) * 2008-10-29 2017-12-27 Mitsubishi Electric Corporation Control device for permanent magnet synchronization electric motor
CA2756719C (en) * 2009-03-26 2014-09-09 Mitsubishi Electric Corporation Controller for ac rotary machine
US8296657B2 (en) * 2009-05-19 2012-10-23 Sony Corporation Random image selection without viewing duplication
JP5414420B2 (ja) * 2009-08-21 2014-02-12 ジェコー株式会社 電流センサ及びその製造方法
DE102009039672B4 (de) * 2009-09-02 2024-03-07 Sew-Eurodrive Gmbh & Co Kg Verfahren zur Bestimmung der Rotorlage einer feldorientiert betriebenen Synchronmaschine
KR101628385B1 (ko) * 2010-03-31 2016-06-08 현대자동차주식회사 영구자석 동기모터의 제어방법
ES2464565T3 (es) 2010-06-29 2014-06-03 Vestas Wind Systems A/S Identificación de perfil de inductancia de generador de imán permanente
EP2453248B1 (en) 2010-11-15 2013-06-26 ABB Oy Method and arrangement for determining inductances of synchronous reluctance machine
US9187100B2 (en) * 2010-12-20 2015-11-17 Cummins Inc. Hybrid power train flexible control integration
KR101209965B1 (ko) * 2010-12-30 2012-12-07 엘에스산전 주식회사 전기자동차의 유도 전동기의 토크 제어 시스템 및 그 방법
JP2013021840A (ja) * 2011-07-12 2013-01-31 Honda Motor Co Ltd リラクタンスモータのロータ
EP2555420B1 (en) * 2011-08-01 2019-10-23 ABB Schweiz AG Self-commissioning procedure for inductance estimation in an electrical machine
KR20130026873A (ko) * 2011-09-06 2013-03-14 현대자동차주식회사 하이브리드 차량의 제어방법
JP5438081B2 (ja) * 2011-09-21 2014-03-12 日立オートモティブシステムズ株式会社 ブラシレスモータの駆動装置
KR101549283B1 (ko) * 2011-10-12 2015-09-01 엘에스산전 주식회사 영구자석 동기 전동기 구동 시스템의 파라미터 추정장치
CA2800792C (en) * 2012-01-06 2016-10-25 Sears Brands, Llc Programmable portable power tool with brushless dc motor
JP5598767B2 (ja) * 2012-02-22 2014-10-01 株式会社デンソー 交流電動機の制御装置
JP5621998B2 (ja) 2012-02-22 2014-11-12 株式会社デンソー 交流電動機の制御装置
JP5483218B2 (ja) 2012-02-22 2014-05-07 株式会社デンソー 交流電動機の制御装置
JP5534252B2 (ja) 2012-02-22 2014-06-25 株式会社デンソー 交流電動機の制御装置
JP5488845B2 (ja) 2012-02-22 2014-05-14 株式会社デンソー 交流電動機の制御装置
DE102012103941B4 (de) * 2012-05-04 2017-01-26 Antriebstechnik Katt Hessen Gmbh Wechselstromwicklung für eine rotierende elektrische Maschine und Verfahren zu ihrer Herstellung aus Runddrähten
GB2503039B (en) * 2012-06-15 2020-05-27 Danfoss Drives As Method for controlling a synchronous reluctance electric motor
GB2503040B (en) * 2012-06-15 2020-05-06 Danfoss Drives As Variable torque angle for electric motor
US10234165B2 (en) * 2012-07-21 2019-03-19 Zhongshan Broad-Ocean Motor Co., Ltd. HVAC control system for household central air conditioning
JP5626306B2 (ja) * 2012-10-09 2014-11-19 トヨタ自動車株式会社 回転電機制御システム
US9871418B2 (en) 2012-11-01 2018-01-16 General Electric Company Sensorless electric machine
US9641033B2 (en) 2013-09-06 2017-05-02 General Electric Company Electric machine having offset rotor sections
US9941775B2 (en) 2012-11-01 2018-04-10 General Electric Company D-ring implementation in skewed rotor assembly
US9906108B2 (en) 2012-11-01 2018-02-27 General Electric Company Sensorless electric machine
US9093878B2 (en) * 2012-11-01 2015-07-28 General Electric Company Sensorless electric machine
GB201301259D0 (en) * 2013-01-24 2013-03-06 Rolls Royce Plc Method of controlling an ac machine and controller for controlling an ac machine
EP3007352B1 (en) * 2013-05-28 2019-04-03 Fuji Electric Co., Ltd. Device for detecting magnetic pole position of permanent magnet-type synchronous motor
US10521519B2 (en) * 2013-07-23 2019-12-31 Atieva, Inc. Induction motor flux and torque control with rotor flux estimation
US9906082B2 (en) 2013-09-06 2018-02-27 General Electric Company Electric machine having reduced torque oscillations and axial thrust
US9641119B2 (en) * 2013-12-10 2017-05-02 Mcmaster University Extended-speed low-ripple torque control of switched reluctance motor drives
CN204068791U (zh) * 2014-03-25 2014-12-31 睿能机电有限公司 一种无刷直流电机的电磁转矩脉动抑制装置
KR101601444B1 (ko) * 2014-07-04 2016-03-21 현대자동차주식회사 모터 구동 시스템의 인버터 6-스텝 제어 장치 및 방법
DE102014214392A1 (de) * 2014-07-23 2016-01-28 Ksb Aktiengesellschaft Herstellungsverfahren für einen Rotor einer Reluktanzmaschine sowie Rotor für eine Reluktanzmaschine
ES2784177T3 (es) * 2014-08-05 2020-09-22 Ego Elektro Geraetebau Gmbh Aparato electrodoméstico
JP6390446B2 (ja) * 2015-01-22 2018-09-19 株式会社デンソー 回転電機の制御装置
US9692327B2 (en) * 2015-04-22 2017-06-27 Infineon Technologies Ag Multi-phase machine current control
US10439535B2 (en) * 2015-04-27 2019-10-08 Mitsubishi Electric Corporation Control device of AC rotating machine and electric power steering device
DE102015005478A1 (de) * 2015-04-28 2016-11-03 Stöber Antriebstechnik GmbH & Co. KG Servoantrieb, Verfahren zur Bestimmung der Position eines Rotors eines Permanentmagnetservomotors sowie Rotor zur Verwendung beim Servoantrieb und beim Verfahren
KR101698775B1 (ko) * 2015-08-11 2017-01-23 엘지전자 주식회사 홈 어플라이언스
US10581274B2 (en) 2015-06-03 2020-03-03 Lg Electronics Inc. Home appliance
KR101684538B1 (ko) * 2015-06-18 2016-12-08 현대자동차 주식회사 하이브리드 차량의 인버터 제어 방법
US10608560B2 (en) * 2015-08-08 2020-03-31 Nidec Corporation Motor controller and a motor system
KR101663520B1 (ko) * 2015-08-11 2016-10-07 엘지전자 주식회사 모터 구동장치 및 이를 구비하는 홈 어플라이언스
JP6088604B1 (ja) * 2015-08-27 2017-03-01 ファナック株式会社 磁束制御部を有する電動機制御装置、ならびに機械学習装置およびその方法
FR3048139B1 (fr) * 2016-02-18 2018-01-26 Schneider Toshiba Inverter Europe Sas Procede et systeme de commande pour une installation de commande de moteur electrique
US10615728B2 (en) * 2016-02-25 2020-04-07 Panasonic Intellectual Property Management Co., Ltd. Motor control device and motor control method
US9831808B1 (en) * 2016-07-20 2017-11-28 Semiconductor Components Industries, Llc Rotor position sensing system for three phase motors and related methods
US9831809B1 (en) * 2016-07-20 2017-11-28 Semiconductor Components Industries, Llc Rotor position sensing system for three phase motors and related methods
CN109642918B (zh) * 2016-08-24 2021-03-09 东芝三菱电机产业系统株式会社 脉冲宽度调制转换器用输入滤波器的上电评价试验装置
CN106438325B (zh) * 2016-11-29 2018-07-17 芜湖美智空调设备有限公司 电器设备及其压缩机状态检测方法、服务器和智能终端
DE102017010523A1 (de) 2016-12-14 2018-06-14 Sew-Eurodrive Gmbh & Co Kg Verfahren zur Bestimmung der Winkellage des Rotors eines von einem Wechselrichter gespeisten Synchronmotors und einer Vorrichtung zur Durchführung des Verfahrens
US9876452B1 (en) 2016-12-16 2018-01-23 Semiconductor Components Industries, Llc Rotor position sensing system for permanent magnet synchronous motors and related methods
WO2018152177A1 (en) * 2017-02-14 2018-08-23 KSR IP Holdings, LLC Systems and methods for harmonic compensation
WO2018173587A1 (ja) * 2017-03-23 2018-09-27 日本電産株式会社 モータ制御方法、モータ制御システムおよび電動パワーステアリングシステム
FR3066276B1 (fr) * 2017-05-10 2019-04-26 Schneider Toshiba Inverter Europe Sas Procede d'identification de la resistance electrique du rotor d'un moteur electrique
RU2658741C1 (ru) * 2017-05-25 2018-06-22 Акционерное общество "Научно-исследовательское предприятие общего машиностроения" (АО "НИПОМ") Способ управлением током возбуждения синхронного электродвигателя в послеаварийных процессах энергосистемы
JP6748896B2 (ja) * 2017-06-16 2020-09-02 パナソニックIpマネジメント株式会社 電動工具の制御装置
KR101939476B1 (ko) * 2017-07-11 2019-01-16 엘지전자 주식회사 모터 구동 장치
TWI661296B (zh) * 2017-11-17 2019-06-01 財團法人資訊工業策進會 監測系統及監測方法
TWI661316B (zh) * 2017-11-17 2019-06-01 財團法人資訊工業策進會 評分系統及評分方法
KR102410944B1 (ko) * 2017-12-01 2022-06-20 현대자동차주식회사 3상 전동기의 임피던스 파라미터 추출 장치 및 방법
CN108155833B (zh) * 2017-12-16 2020-09-18 南京理工大学 考虑电气特性的电机伺服系统渐近稳定控制方法
US11499537B2 (en) * 2017-12-17 2022-11-15 Microchip Technology Incorporated Closed loop torque compensation for compressor applications
CN108199541B (zh) * 2018-01-31 2019-12-06 山东大学 抑制分数槽永磁同步电机固有轴电压的极槽配合选择方法
US10903771B2 (en) * 2018-03-13 2021-01-26 Board Of Regente, The University Of Texas System System, apparatus and method for high torque density permanent magnet switched reluctance machine
US10944352B2 (en) 2018-03-19 2021-03-09 Tula eTechnology, Inc. Boosted converter for pulsed electric machine control
US11623529B2 (en) 2018-03-19 2023-04-11 Tula eTechnology, Inc. Pulse modulated control with field weakening for improved motor efficiency
US10742155B2 (en) 2018-03-19 2020-08-11 Tula eTechnology, Inc. Pulsed electric machine control
DE102018211882A1 (de) * 2018-07-17 2020-01-23 Ziehl-Abegg Se Verfahren zur Erfassung von Betriebsparameterdaten eines Motorsystems mit Elektromotor und entsprechendes Motorsystem
CN109150055B (zh) * 2018-09-14 2020-04-10 清华大学 永磁同步电机的i/f控制中电磁转矩反馈控制系统及方法
DE102018008689A1 (de) * 2018-10-24 2020-04-30 Kostal Drives Technology Gmbh Verfahren und Vorrichtung zur lastfreien Bestimmung lastabhängiger Lagezuordnungsparameter einer Synchronmaschine ohne Lagegeber
CN109387781B (zh) * 2018-12-04 2023-11-24 北京交通大学 开关器件损伤评估方法、装置和设备
EP3896826A4 (en) * 2018-12-14 2022-08-31 NSK Ltd. VIBRATION CONTROL DEVICE
CN109921720B (zh) * 2019-03-07 2021-09-07 Oppo广东移动通信有限公司 电机的控制装置、电子设备及其控制方法
CN109818541B (zh) * 2019-03-13 2020-10-02 东南大学 一种用于磁链观测的记忆电机绕组复用控制方法及系统
JP7308631B2 (ja) * 2019-03-15 2023-07-14 株式会社日立製作所 情報連携システムおよび情報管理方法
CN109782173B (zh) * 2019-03-25 2021-07-16 中车青岛四方车辆研究所有限公司 异步电机励磁互感曲线测量系统及其测量方法
GB2583947B (en) * 2019-05-15 2023-11-29 Trw Ltd An electric motor apparatus
CA3141827A1 (en) * 2019-05-30 2020-12-03 Magna International Inc. Motor drive optimization system and method
CN110429895B (zh) * 2019-07-26 2021-01-15 江苏大学 混合动力车用开关磁阻bsg优化线性控制器的构造方法
CN112468047B (zh) * 2019-09-06 2022-07-05 中车株洲电力机车研究所有限公司 一种电机参数测试方法以及一种直线电机的控制方法
CN110707968B (zh) * 2019-09-10 2021-11-30 苏州安驰控制系统有限公司 一种单锭控制系统的控制方法、系统及计算机存储介质
EP3809578A1 (en) * 2019-10-15 2021-04-21 ABB Schweiz AG Method for starting a synchronous motor
TWI723605B (zh) * 2019-10-31 2021-04-01 國立中央大學 馬達控制系統及其最大功率因數控制器的控制方法
CN111106780B (zh) * 2019-12-04 2021-04-30 中国科学院深圳先进技术研究院 一种电机控制方法、装置、终端设备及存储介质
CN111092582B (zh) * 2019-12-31 2021-10-08 潍柴动力股份有限公司 一种位置识别方法和相关装置
DE102020201413A1 (de) 2020-02-05 2021-08-05 Zf Friedrichshafen Ag Verfahren und Steuergerät zum Betreiben einer elektrischen Maschine
CN111510040B (zh) * 2020-04-10 2023-10-10 江苏大学 用于五相永磁同步电机的无权重系数的双虚拟电压矢量预测转矩控制方法
WO2021216228A1 (en) * 2020-04-21 2021-10-28 Tula eTechnology, Inc. Pulse modulated control with field weakening for improved motor efficiency
CN111585493B (zh) * 2020-06-08 2021-12-21 福州大学 一种基于无永磁磁链观测器的定子永磁型记忆电机电流控制方法
CN112039374B (zh) * 2020-08-05 2022-04-08 青岛大学 考虑输入饱和的永磁同步电动机命令滤波离散控制方法
US11239780B1 (en) * 2020-09-24 2022-02-01 Rockwell Automation Technologies, Inc. Iterative flux identification
CN114594704B (zh) * 2020-12-07 2023-08-01 山东新松工业软件研究院股份有限公司 电机内环控制方法、装置及控制器
CN112636659A (zh) * 2020-12-17 2021-04-09 宁波拓普集团股份有限公司 一种基于电机控制指令曲线的电机控制方法
US11515824B2 (en) * 2020-12-21 2022-11-29 Caterpillar Inc. Switched reluctance self sensing active pulse torque compensation
CN112737160B (zh) * 2020-12-29 2022-04-19 东南大学 一种集中绕组外转子游标电机提升功率因数的方法
US11628730B2 (en) 2021-01-26 2023-04-18 Tula eTechnology, Inc. Pulsed electric machine control
DE102021104785A1 (de) 2021-03-01 2022-09-01 Audi Aktiengesellschaft Läufer für eine Synchronmaschine
WO2022197742A1 (en) 2021-03-15 2022-09-22 Tula eTechnology, Inc. Methods of optimizing waveforms for electric motors
CN112821835B (zh) * 2021-03-30 2022-09-09 宁波奥克斯电气股份有限公司 参数确定方法、装置及空调器
CN113098358B (zh) * 2021-04-26 2022-07-26 湖南科技大学 基于导通相实时电感定位的开关磁阻电机无位置传感器控制方法
CN113452287B (zh) * 2021-06-07 2023-03-28 江苏科技大学 一种水下航行器的多永磁同步电机的控制方法及控制系统
WO2022265986A1 (en) 2021-06-14 2022-12-22 Tula eTechnology, Inc. Electric machines with efficient torque transitions
WO2023278139A1 (en) 2021-06-28 2023-01-05 Tula eTechnology, Inc. Selective phase control of an electric machine
US20230006575A1 (en) * 2021-06-30 2023-01-05 Texas Instruments Incorporated Method for determining a position of a rotor at standstill
US11557996B1 (en) 2021-07-08 2023-01-17 Tula eTechnology, Inc. Methods of reducing vibrations for electric motors
US11345241B1 (en) 2021-08-12 2022-05-31 Tula eTechnology, Inc. Method of optimizing system efficiency for battery powered electric motors
WO2023038760A1 (en) 2021-09-08 2023-03-16 Tula eTechnology, Inc. Electric machine torque adjustment based on waveform integer multiples
WO2023069131A1 (en) 2021-10-18 2023-04-27 Tula eTechnology, Inc. Mechanical and electromechanical arrangements for field-weakening of an electric machine that utilizes permanent magnets
AT525623A1 (de) * 2021-11-12 2023-05-15 Avl List Gmbh Kontrollverfahren für eine Kontrolle von Betriebsparametern eines Elektromotors, insbesondere zum Antrieb eines Fahrzeugs
WO2023147088A1 (en) * 2022-01-28 2023-08-03 Tau Motors, Inc. Motor control using piecewise affine model
US11888424B1 (en) 2022-07-18 2024-01-30 Tula eTechnology, Inc. Methods for improving rate of rise of torque in electric machines with stator current biasing
CN115356624B (zh) * 2022-08-15 2023-07-18 小米汽车科技有限公司 电机铁耗确定方法、装置、车辆、存储介质及芯片

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5752392A (en) * 1980-09-10 1982-03-27 Fuji Electric Co Ltd Current control for ac machine
JPH07322700A (ja) * 1994-05-27 1995-12-08 Meidensha Corp 誘導電動機の最大効率制御方法
JPH11150996A (ja) * 1997-11-13 1999-06-02 Toyota Motor Corp モータ制御装置
JPH11205914A (ja) * 1998-01-12 1999-07-30 Yamaha Motor Co Ltd 電動車両の出力制御装置
JP2001161099A (ja) * 1999-11-30 2001-06-12 Meidensha Corp 同期電動機の制御方式
JP2001204200A (ja) * 2000-01-17 2001-07-27 Toyo Electric Mfg Co Ltd 永久磁石形回転電機の制御方法
JP2002186271A (ja) * 2000-12-13 2002-06-28 Toyota Central Res & Dev Lab Inc 電気系解析方法、記録媒体および電気系解析システム
JP2003070280A (ja) * 2001-08-28 2003-03-07 Hitachi Ltd モータ制御装置
JP2005328635A (ja) * 2004-05-14 2005-11-24 Toyo Electric Mfg Co Ltd スイッチトリラクタンスモータの制御装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0982844B1 (en) * 1996-08-19 2009-06-03 Daikin Industries, Limited Synchronous motor driving device
JP3336870B2 (ja) * 1996-09-04 2002-10-21 三菱電機株式会社 回転磁石形多相同期電動機の制御方法及びその装置
JP3813637B2 (ja) * 1998-09-03 2006-08-23 三菱電機株式会社 交流電動機の制御装置
GB2357201B (en) * 1999-08-20 2003-07-09 Mitsubishi Electric Corp Method and apparatus for controlling a synchronous motor
EP1211798B1 (en) * 2000-11-22 2018-01-10 Nissan Motor Co., Ltd. Motor control apparatus and motor control method
JP3771544B2 (ja) 2003-03-24 2006-04-26 株式会社日立製作所 永久磁石形同期電動機の制御方法及び装置
JP3928575B2 (ja) * 2003-04-07 2007-06-13 日産自動車株式会社 モーター制御装置
KR100655702B1 (ko) * 2004-12-20 2006-12-11 현대자동차주식회사 영구자석 동기 모터 제어방법

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5752392A (en) * 1980-09-10 1982-03-27 Fuji Electric Co Ltd Current control for ac machine
JPH07322700A (ja) * 1994-05-27 1995-12-08 Meidensha Corp 誘導電動機の最大効率制御方法
JPH11150996A (ja) * 1997-11-13 1999-06-02 Toyota Motor Corp モータ制御装置
JPH11205914A (ja) * 1998-01-12 1999-07-30 Yamaha Motor Co Ltd 電動車両の出力制御装置
JP2001161099A (ja) * 1999-11-30 2001-06-12 Meidensha Corp 同期電動機の制御方式
JP2001204200A (ja) * 2000-01-17 2001-07-27 Toyo Electric Mfg Co Ltd 永久磁石形回転電機の制御方法
JP2002186271A (ja) * 2000-12-13 2002-06-28 Toyota Central Res & Dev Lab Inc 電気系解析方法、記録媒体および電気系解析システム
JP2003070280A (ja) * 2001-08-28 2003-03-07 Hitachi Ltd モータ制御装置
JP2005328635A (ja) * 2004-05-14 2005-11-24 Toyo Electric Mfg Co Ltd スイッチトリラクタンスモータの制御装置

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010041839A (ja) * 2008-08-06 2010-02-18 Fuji Electric Systems Co Ltd 永久磁石形同期電動機の制御装置
WO2010029879A1 (ja) * 2008-09-11 2010-03-18 本田技研工業株式会社 電動パワーステアリング装置
US8606462B2 (en) 2008-09-11 2013-12-10 Honda Motor Co., Ltd. Electric power steering device
DE112010001465T5 (de) 2009-03-30 2012-05-16 Hitachi, Ltd. Wechselstrommotor-Steuervorrichtung und Wechselstrommotor-Treibersystem
US8742704B2 (en) 2009-03-30 2014-06-03 Hitachi, Ltd. AC motor control device and AC motor driving system
KR101328192B1 (ko) * 2009-06-23 2013-11-08 가트너-일렉트로닉-디자인 게엠베하 자석과 함께 고정자 및 회전자를 구비한 브러시리스 직류 모터(bldc) 내의 부하율을 보상하기 위한 장치
EP2290806A2 (en) 2009-08-27 2011-03-02 Hitachi, Ltd. Control apparatus for synchronous motor
CN102013870A (zh) * 2010-11-11 2011-04-13 江苏大学 五自由度无轴承同步磁阻电机逆系统解耦控制器
WO2015132991A1 (ja) * 2014-03-05 2015-09-11 三菱電機株式会社 シンクロナスリラクタンスモータ
JP5866074B1 (ja) * 2014-03-05 2016-02-17 三菱電機株式会社 シンクロナスリラクタンスモータ
US9825515B2 (en) 2014-03-05 2017-11-21 Mitsubishi Electric Corporation Synchronous reluctance motor having radial-direction widths slit configuration on a q-axis for improved power factor
KR101819848B1 (ko) 2014-03-05 2018-01-17 미쓰비시덴키 가부시키가이샤 싱크로너스 리럭턴스 모터
JP2015192495A (ja) * 2014-03-27 2015-11-02 株式会社豊田中央研究所 情報処理装置、情報記憶装置、及び回転電機の制御装置
JP2016005427A (ja) * 2014-06-16 2016-01-12 現代自動車株式会社Hyundaimotor Company 電動機のセンサレス制御方法およびシステム
JP2016140139A (ja) * 2015-01-26 2016-08-04 株式会社デンソー 回転電機の制御装置
KR20170008381A (ko) * 2015-07-13 2017-01-24 현대모비스 주식회사 모터 속도 측정 장치 및 방법
KR102419489B1 (ko) 2015-07-13 2022-07-12 현대모비스 주식회사 모터 속도 측정 장치 및 방법
JP2017085879A (ja) * 2015-10-13 2017-05-18 シュネーデル、トウシバ、インベーター、ヨーロッパ、ソシエテ、パル、アクション、セプリフエSchneider Toshiba Inverter Europe Sas 可変リラクタンス同期電気モータのインダクタンス値を識別するための制御プロセス
JP2018085852A (ja) * 2016-11-24 2018-05-31 日産自動車株式会社 可変磁束型回転電機の電流制御方法、及び電流制御装置
JP2018085851A (ja) * 2016-11-24 2018-05-31 日産自動車株式会社 可変磁束モータの電流制御方法、及び電流制御装置
KR20190016377A (ko) * 2017-08-08 2019-02-18 엘지전자 주식회사 모터 구동 장치 및 이를 구비하는 공기조화기
KR101990444B1 (ko) * 2017-08-08 2019-06-18 엘지전자 주식회사 모터 구동 장치 및 이를 구비하는 공기조화기
JP2019068601A (ja) * 2017-09-29 2019-04-25 日産自動車株式会社 モータの制御方法、及び、モータの制御装置。
JP2019068600A (ja) * 2017-09-29 2019-04-25 日産自動車株式会社 モータの制御方法、及び、モータの制御装置
US11296633B2 (en) 2018-01-12 2022-04-05 Mitsubishi Electric Corporation Rotary machine control device
CN110445439A (zh) * 2018-05-04 2019-11-12 美的集团股份有限公司 永磁同步电机的控制方法和装置
CN110460284A (zh) * 2018-05-04 2019-11-15 美的集团股份有限公司 食物料理机、电机驱动系统、感应电机的控制方法和装置
WO2020075620A1 (ja) * 2018-10-10 2020-04-16 日本電産株式会社 モータ制御装置およびパワーステアリング装置
CN110829919A (zh) * 2019-10-29 2020-02-21 西安石油大学 一种连续波泥浆脉冲器的控制系统及其开发方法
KR102500019B1 (ko) * 2022-12-29 2023-02-16 주식회사 에스엠전자 2승저감부하에 연결된 모터의 운영전력효율이 개선된 제어관제시스템

Also Published As

Publication number Publication date
US7772790B2 (en) 2010-08-10
CN101192808B (zh) 2010-09-15
DE102007057499A1 (de) 2008-06-05
DE102007057499B4 (de) 2017-06-08
JP5167631B2 (ja) 2013-03-21
CN101192808A (zh) 2008-06-04
US20080129243A1 (en) 2008-06-05

Similar Documents

Publication Publication Date Title
JP5167631B2 (ja) モータの制御方法及びそれを利用するモータ制御装置
Kwon et al. Extending operational limit of IPMSM in signal-injection sensorless control by manipulation of convergence point
JP3467961B2 (ja) 回転電機の制御装置
CN102170256B (zh) 电机控制器
Wang et al. Large-scale optimization of synchronous reluctance machines using CE-FEA and differential evolution
JP3809783B2 (ja) モータ制御装置
JP2008295200A (ja) 同期モータの制御装置、及び同期モータ制御の最適化方法
Barcaro et al. Remarks on torque estimation accuracy in fractional-slot permanent-magnet motors
Nie et al. Deadbeat-direct torque and flux control for wound field synchronous machines
Hamiti et al. Modeling of a synchronous reluctance machine accounting for space harmonics in view of torque ripple minimization
JP2009268268A (ja) モータ制御装置及び発電機制御装置
JP2004032944A (ja) 同期電動機の制御装置及び同期電動機
JP6806983B2 (ja) 二重三相巻線永久磁石同期形電動機の駆動システム
JP6687228B1 (ja) 交流回転電機の制御装置
CN104094516A (zh) 交流电机的控制装置
US20220399838A1 (en) Induction motor and control device therefor
JP4996847B2 (ja) サーボモータにおける電流制御方法,電流制御プログラム,記録媒体,サーボモータ
Arafat et al. Fault tolerant control of five-phase permanent magnet assisted synchronous reluctance motor based on dynamic current phase advance
WO2017030055A1 (ja) 回転機の制御装置および制御方法
Yang et al. Torque ripple reduction strategy of model based predictive torque control for doubly salient permanent magnet synchronous machines
CN114865967A (zh) 交流旋转电机的控制装置
Oliveira et al. Extended DQ transformation for vectorial control applications of non-sinusoidal permanent magnet AC machines
CN101373945B (zh) 电动机控制装置以及电流相位的控制方法
JP5456873B1 (ja) 同期機制御装置
JP2018046723A (ja) 回転機制御装置及び回転機制御方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20081212

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110907

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110920

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20111111

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20120726

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121023

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20121030

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20121127

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121210

R151 Written notification of patent or utility model registration

Ref document number: 5167631

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160111

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees