JP5155344B2 - 電動機の磁極位置推定装置 - Google Patents

電動機の磁極位置推定装置 Download PDF

Info

Publication number
JP5155344B2
JP5155344B2 JP2010007019A JP2010007019A JP5155344B2 JP 5155344 B2 JP5155344 B2 JP 5155344B2 JP 2010007019 A JP2010007019 A JP 2010007019A JP 2010007019 A JP2010007019 A JP 2010007019A JP 5155344 B2 JP5155344 B2 JP 5155344B2
Authority
JP
Japan
Prior art keywords
magnetic pole
pole position
current
motor
difference
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.)
Expired - Fee Related
Application number
JP2010007019A
Other languages
English (en)
Other versions
JP2011147287A (ja
Inventor
雅彦 秋山
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.)
Honda Motor Co Ltd
Original Assignee
Honda Motor Co Ltd
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 Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2010007019A priority Critical patent/JP5155344B2/ja
Priority to US12/968,296 priority patent/US8384323B2/en
Publication of JP2011147287A publication Critical patent/JP2011147287A/ja
Application granted granted Critical
Publication of JP5155344B2 publication Critical patent/JP5155344B2/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
    • 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
    • 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
    • H02P21/18Estimation of position or speed
    • 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
    • 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
    • H02P6/181Circuit arrangements for detecting position without separate position detecting elements using different methods depending on the speed
    • 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
    • H02P6/186Circuit arrangements for detecting position without separate position detecting elements using difference of inductance or reluctance between the phases
    • 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
    • H02P6/188Circuit arrangements for detecting position without separate position detecting elements using the voltage difference between the 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
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/01Motor rotor position determination based on the detected or calculated phase inductance, e.g. for a Switched Reluctance 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
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/05Determination of the rotor position by using two different methods and/or motor models

Landscapes

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

Description

本発明は、電動機の磁極位置推定装置に関する。
従来、例えば突極性を有する同期モータにおいて、実軸であるdq座標軸に対して位相差を有する制御軸(dc−qc座標軸)での電圧方程式に基づいて、誘起電圧が位相に応じて変化することを用いて位相差を算出し、この算出結果から制御軸の位相を算出する装置が知られている(例えば、特許文献1参照)。
特開2001−251889号公報
ところで、上記従来技術に係る装置においては、モータの停止時を含む回転速度が低い場合には位相差を算出することができないことから、他の演算方法を用いて位相差を算出する必要が生じる。
このため、例えば、複数の異なる演算方法により位相差を算出する複数の演算装置を備え、モータの回転速度に応じて複数の演算装置を切り替えて位相差を算出するように構成すると、装置構成が複雑化して、演算負荷が増大してしまい、複数の演算装置を切り替えるときにチャタリングが発生するという問題が生じる。
本発明は上記事情に鑑みてなされたもので、装置構成の複雑化および演算負荷の増大を抑制し、磁極位置の算出結果に変動が生じることを防止することが可能な電動機の磁極位置推定装置を提供することを目的とする。
上記課題を解決して係る目的を達成するために、本発明の第1態様に係る電動機の磁極位置推定装置は、電動機への通電に対する指令電圧を出力する指令電圧出力手段(例えば、実施の形態での電圧加算部45)と、前記電動機に通電される実電流を検出して検出結果を出力する電流検出手段(例えば、実施の形態での相電流センサ32)と、前記指令電圧と前記実電流とに基づいて、前記電動機の所定モデルの電圧方程式に応じたモデル電流を算出して、算出結果を出力するモデル電流算出手段(例えば、実施の形態でのモータモデル電流演算部61)と、前記モデル電流と前記実電流との差である電流差を算出する電流差算出手段(例えば、実施の形態でのγ軸電流偏差算出部62aおよびδ軸電流偏差算出部62b)と、前記電流差に基づいて、前記電動機の実際の磁極位置と推定又は指定による磁極位置との差である磁極位置差を、前記電動機の回転速度が所定値未満であるか否かに応じて、互いに異なる演算により算出する第1及び第2の磁極位置差算出手段(例えば、実施の形態での低速域磁極位置誤差演算部66及び高速域磁極位置誤差演算部65)と、前記第1の磁極位置差算出手段又は前記第2の磁極位置差算出手段により算出された前記磁極位置差に基づいて、前記電動機の磁極位置を演算する磁極位置演算手段(例えば、実施の形態での回転速度−磁極位置演算部49)とを備える。
さらに前記第1の磁極位置差算出手段は、前記電動機の回転速度が所定値未満である場合に、前記電動機のインダクタンスが前記磁極位置に応じて変化することに基づいて
dq座標系での前記電動機の電圧方程式においてd軸をγ軸に置換し、かつq軸をδ軸に置換し、かつ回転速度をゼロとしてから該電圧方程式を離散化して得られる離散時間状態方程式と、該離散時間状態方程式をdq座標系に対して位相差を有するγδ座標系での前記電動機の電圧方程式に適用してから回転速度をゼロとして得られる式とに基づき、今回の制御サイクルでの前記電流差および前記実電流と前回の制御サイクルでの前記実電流とを用いて前記磁極位置差を算出し、前記第2の磁極位置差算出手段は、前記電動機の回転速度が所定値以上である場合に、前記電動機の誘起電圧が前記磁極位置に応じて変化することに基づいて、前記γδ座標系での前記電動機の電圧方程式において磁極位置および回転速度に係る項を外乱として削除する前後において該電圧方程式を離散化して得られる2つの離散時間状態方程式に基づき、今回の制御サイクルでの前記電流差と前回の制御サイクルでの前記実電流および前記回転速度とを用いて前記磁極位置差を算出する。
さらに、本発明の第態様に係る電動機の磁極位置推定装置は、前記第1の磁極位置差算出手段は、前記磁極位置差が小さいこと基づいて得られる近似式により前記磁極位置差を算出する。
本発明の第1態様に係る電動機の磁極位置推定装置によれば、モデル電流算出手段は、電動機の回転速度に応じて互いに異なる演算により磁極位置差を算出する2つの第1の磁極位置差算出手段と第2の磁極位置差算出手段とに対して、共通化された状態方程式によりモデル電流を算出する。また、磁極位置演算手段は、2つの第1の磁極位置差算出手段と第2の磁極位置差算出手段とに対して、共通化された演算によって電動機の磁極位置を算出する。
これにより、例えば第1の磁極位置差算出手段と第2の磁極位置差算出手段とに対して、共通化されていない演算によってモデル電流および電動機の磁極位置を算出する場合に比べて、装置構成を簡略化することができ、演算負荷を低減することができる。さらに、電動機の回転速度が所定値未満であるか否かに応じて、第1の磁極位置差算出手段による磁極位置差の算出と、第2の磁極位置差算出手段による磁極位置差の算出との、何れか一方を選択して切り替えるときの安定性を向上させることができる。
さらに、電動機の回転速度が所定値未満であるか否かに応じて、第1の磁極位置差算出手段による磁極位置差の算出と、第2の磁極位置差算出手段による磁極位置差の算出との、何れか一方を選択して切り替えることにより、適切に電動機の磁極位置を算出することができる。第1の磁極位置差算出手段によれば、電動機の停止時であっても磁極位置を算出することができる。第2の磁極位置差算出手段によれば、所定の検出用電圧を印加する必要無しに磁極位置を算出することができる。
本発明の実施形態に係る電動機の磁極位置推定装置の構成図である。 本発明の実施形態に係る回転直交座標のγ−δ軸およびd−q軸の一例を示す図である。 本発明の実施形態に係る電動機の磁極位置推定装置の構成図である。 図3に示す磁極位置誤差推定部の構成図である。 本発明の実施形態に係る切替フラグの回転速度に応じた変化の一例を示す図である。 図3に示す回転速度−磁極位置演算部の構成図である。 本発明の実施形態の変形例に係る回転速度−磁極位置演算部の構成図である。
以下、本発明の電動機の磁極位置推定装置の実施形態について添付図面を参照しながら説明する。
この実施形態による電動機の磁極位置推定装置10(以下、単に、磁極位置推定装置10と呼ぶ)は、例えば3相交流のブラシレスDCモータ11(以下、単に、モータ11と呼ぶ)の磁極位置(つまり、所定の基準回転位置からの回転子の磁極の回転角度)を推定するものであって、このモータ11は、界磁に利用する永久磁石を有する回転子(図示略)と、この回転子を回転させる回転磁界を発生する固定子(図示略)とを備えて構成されている。
磁極位置推定装置10は、例えば図1に示すように、バッテリ12を直流電源とするインバータ13と、モータ制御装置14とを備えて構成されている。
この3相(例えば、U相、V相、W相の3相)交流のモータ11の駆動はモータ制御装置14から出力される制御指令を受けてインバータ13によりおこなわれる。
インバータ13は、スイッチング素子(例えば、MOSFET:Metal Oxide Semi-conductor Field Effect Transistor)を複数用いてブリッジ接続してなるブリッジ回路13aと平滑コンデンサCとを具備し、このブリッジ回路13aがパルス幅変調(PWM)された信号によって駆動される。
このブリッジ回路13aでは、例えば各相毎に対をなすハイ側およびロー側U相トランジスタUH,ULと、ハイ側およびロー側V相トランジスタVH,VLと、ハイ側およびロー側W相トランジスタWH,WLとがブリッジ接続されている。そして、各トランジスタUH,VH,WHはドレインがバッテリ12の正極側端子に接続されてハイサイドアームを構成し、各トランジスタUL,VL,WLはソースがバッテリ12の接地された負極側端子に接続されてローサイドアームを構成している。そして、各相毎に、ハイサイドアームの各トランジスタUH,VH,WHのソースはローサイドアームの各トランジスタUL,VL,WLのドレインに接続され、各トランジスタUH,UL,VH,VL,WH,WLのドレイン−ソース間には、ソースからドレインに向けて順方向となるようにして、各ダイオードDUH,DUL,DVH,DVL,DWH,DWLが接続されている。
インバータ13は、例えばモータ11の駆動時等においてモータ制御装置14から出力されて各トランジスタUH,VH,WH,UL,VL,WLのゲートに入力されるスイッチング指令であるゲート信号(つまり、PWM信号)に基づき、各相毎に対をなす各トランジスタのオン(導通)/オフ(遮断)状態を切り替える。これによって、バッテリ12から供給される直流電力を3相交流電力に変換し、モータ11の3相のステータ巻線への通電を順次転流させることで、各相のステータ巻線に交流のU相電流IuおよびV相電流IvおよびW相電流Iwを通電する。
モータ制御装置14は、後述するように、回転直交座標をなすγ−δ座標上で電流のフィードバック制御(ベクトル制御)を行うものであり、指令γ軸電流Iγc及び指令δ軸電流Iδcを演算し、指令γ軸電流Iγc及び指令δ軸電流Iδcに基づいて各相電圧指令Vu,Vv,Vwを算出し、各相電圧指令Vu,Vv,Vwに応じてインバータ13に対するゲート信号であるPWM信号を出力する。そして、実際にインバータ13からモータ11に供給される各相電流Iu,Iv,Iwをγ−δ座標上に変換して得たγ軸電流Iγ及びδ軸電流Iδと、指令γ軸電流Iγc及び指令δ軸電流Iδcとの各偏差がゼロとなるように制御をおこなう。
モータ制御装置14は、例えば相電流センサI/F(インターフェース)21と、制御装置22と、PWM信号生成部23と、検出用電圧印加部24とを備えて構成されている。
相電流センサI/F(インターフェース)21は、インバータ13のブリッジ回路13aとモータ11との間において、3相の各相電流のうち少なくとも何れか2相の各相電流(例えば、U相電流およびV相電流)を検出する各相電流センサ32に接続され、各相電流センサ32から出力される検出信号を受けて、検出信号を制御装置22に出力する。
制御装置22は、例えば図2に示すように、実際のモータ11が有する回転直交座標のd−q軸に対して、実際の回転角と推定または指定した回転角との差である位相差Δθeおよび回転速度ωeを有する回転直交座標のγ−δ軸を設定し、このγ−δ座標上で電流のフィードバック制御(ベクトル制御)をおこなう。
制御装置22は、指令γ軸電流Iγc及び指令δ軸電流Iδcを生成し、指令γ軸電流Iγc及び指令δ軸電流Iδcに基づいて各相電圧指令Vu,Vv,Vwを算出し、PWM信号生成部23に出力する。
また、制御装置22は、各相電流センサ32から出力される検出信号に応じた各相電流Iu,Iv,Iwをγ−δ座標上に変換して得たγ軸電流Iγ及びδ軸電流Iδと、指令γ軸電流Iγc及び指令δ軸電流Iδcとの各偏差がゼロとなるように電流のフィードバック制御(ベクトル制御)をおこなう。
なお、この制御装置22の動作の詳細は後述する。
PWM信号生成部23は、正弦波状の電流を3相のステータ巻線に通電するために、各相電圧指令Vu,Vv,Vwと、三角波などのキャリア信号とを比較して、インバータ13の各トランジスタUH,VH,WH,UL,VL,WLをオン/オフ駆動させるゲート信号(つまり、PWM信号)を生成する。そして、インバータ13において3相の各相毎に対をなす各トランジスタのオン(導通)/オフ(遮断)状態を切り替えることによって、バッテリ12から供給される直流電力を3相交流電力に変換し、3相のモータ11の各ステータ巻線への通電を順次転流させることで、各ステータ巻線に交流のU相電流IuおよびV相電流IvおよびW相電流Iwを通電する。
検出用電圧印加部24は、例えばモータ11の低速域での回転時における磁極位置誤差推定の処理において、後述する制御装置22の電流制御部43から出力されるγ軸電圧指令値Vγcおよびδ軸電圧指令値Vδcに重畳される電圧として、例えば所定周波数のパルス電圧または高周波正弦波電圧などからなるγ軸検出用電圧Vhγおよびδ軸検出用電圧Vhδを出力する。
制御装置22は、例えば図3に示すように、速度制御部41と、指令電流生成部42と、電流制御部43と、印加電圧切替部44と、電圧加算部45と、γδ−3相変換部46と、3相−γδ変換部47と、磁極位置誤差推定部48と、回転速度−磁極位置演算部49と、電気角−機械角変換部50とを備えて構成されている。
速度制御部41は、外部から入力される回転速度指令値ωrcに基づき、例えば電気角−機械角変換部50から出力される回転速度ωr(機械角)に応じたクローズループ制御により、トルク指令Tcを演算し、演算結果を出力する。
なお、制御装置22は、この速度制御部41の代わりに、モータ11から出力されるトルクに応じたトルク制御(クローズループ制御など)を行なうトルク制御部を備えて、トルク制御を実行してもよい。
指令電流生成部42は、速度制御部41から出力されるトルク指令Tcに基づき、指令δ軸電流Iδcおよび指令γ軸電流Iγcを演算して出力する。
電流制御部43は、指令電流生成部42から出力される指令γ軸電流Iγcと3相−γδ変換部47から出力されるγ軸電流Iγとの偏差ΔIγを算出し、指令電流生成部42から出力される指令δ軸電流Iδcと3相−γδ変換部47から出力されるδ軸電流Iδとの偏差ΔIδを算出する。そして、例えばPI(比例・積分)動作などにより、偏差ΔIγを制御増幅してγ軸電圧指令値Vγcを算出し、偏差ΔIδを制御増幅してδ軸電圧指令値Vδcを算出する。そして、γ軸電圧指令値Vγcおよびδ軸電圧指令値Vδcを出力する。
印加電圧切替部44は、回転速度−磁極位置演算部49から出力されるモータ11の回転速度推定値ωeに応じて、印加電圧切替部44から出力する印加電圧を切り替える。
印加電圧切替部44は、例えば回転速度推定値ωeが所定値以下であって、モータ11が低速域で回転している場合には、切替フラグのフラグ値を「0」として、検出用電圧印加部24から出力されるγ軸検出用電圧Vhγおよびδ軸検出用電圧Vhδを印加電圧として出力する。一方、例えば回転速度推定値ωeが所定値よりも大きく、モータ11が高速域で回転している場合には、切替フラグのフラグ値を「1」として、ゼロを印加電圧として出力する。
電圧加算部45は、電流制御部43から出力される各γ軸電圧指令値Vγcおよびδ軸電圧指令値Vδcと、印加電圧切替部44から出力される印加電圧(つまり、検出用電圧印加部24から出力されるγ軸検出用電圧Vhγおよびδ軸検出用電圧Vhδ、又は、ゼロ)とを加算して、各γ軸電圧指令値Vγおよびδ軸電圧指令値Vδとして、γδ−3相変換部46および磁極位置誤差推定部48に入力する。
γδ−3相変換部46は、回転速度−磁極位置演算部49から出力されるモータ11の磁極位置推定値θeにより、γ−δ座標上でのγ軸電圧指令値Vγおよびδ軸電圧指令値Vδを、静止座標である3相交流座標上での電圧指令値であるU相電圧指令VuおよびV相電圧指令VvおよびW相電圧指令Vwに変換する。
3相−γδ変換部47は、相電流センサI/F(インターフェース)21から出力される各相電流Iu,Ivの検出信号に基づき、同一タイミングでの各相電流の電流値の総和はゼロであることを用いて、2相の相電流(例えば、各相電流Iu,Iv)の電流値から、他の1相の相電流(例えば、W相電流Iw)の電流値を算出する。そして、回転速度−磁極位置演算部49から出力されるモータ11の磁極位置推定値θeにより、各相電流Iu,Iv,Iwを、γ−δ座標上でのγ軸電流Iγ及びδ軸電流Iδに変換する。
磁極位置誤差推定部48は、例えば、電圧加算部45から出力されるγ軸電圧指令値Vγ及びδ軸電圧指令値Vδと、3相−γδ変換部47から出力されるγ軸電流Iγ及びδ軸電流Iδとに基づき、モータ11の回転時にモータ11が発生する誘起電圧又はモータ11のインダクタンスが磁極位置によって変化することを利用して位相差Δθeを推定する。
磁極位置誤差推定部48は、例えば、高速域ではモータ11が発生する誘起電圧が磁極位置によって変化することを利用して、また、低速域では電流制御部43から出力されるγ軸電圧指令値Vγc及びδ軸電圧指令値Vδcに検出用電圧印加部24から出力されるγ軸検出用電圧Vhγ及びδ軸検出用電圧Vhδを印加した状態でインダクタンスが磁極位置により変化すること(いわゆるインダクタンスの突極性)を利用して位相差Δθeを推定する。
磁極位置誤差推定部48は、例えば図4に示すように、モータモデル電流演算部61と、γ軸電流偏差算出部62aおよびδ軸電流偏差算出部62bと、出力切替部63と、ヒステリシス設定部64と、高速域磁極位置誤差演算部65と、低速域磁極位置誤差演算部66とを備えて構成されている。
モータモデル電流演算部61は、電圧加算部45から出力されるγ軸電圧指令値Vγ及びδ軸電圧指令値Vδと、3相−γδ変換部47から出力されるγ軸電流Iγ及びδ軸電流Iδとを用いて、モータ11の所定モデルの電圧方程式によりγ軸モデル電流IγM及びδ軸モデル電流IδMを算出し、各モデル電流IγM,IδMを出力する。
モータモデル電流演算部61は、モータ11が低速域で回転しているときに低速域磁極位置誤差演算部66によって位相差Δθeが推定される場合と、モータ11が高速域で回転しているときに高速域磁極位置誤差演算部65によって位相差Δθeが推定される場合との両方に対して、モータ11の所定モデルの電圧方程式が共通化されるように構成されている。
先ず、以下に、モータ11の所定モデルの電圧方程式について説明する。
γ−δ座標系でのモータ11の電圧方程式(つまり、連続時間での電圧方程式)は、例えば、γ軸電圧指令値Vγ及びδ軸電圧指令値Vδと、巻線抵抗Rと、微分演算子pと、d軸インダクタンスLdおよびq軸インダクタンスLqと、位相差Δθeと、回転速度推定値ωeと、d軸電流Id及びq軸電流Iqと、永久磁石の磁束成分(誘起電圧定数)φとにより、例えば下記数式(1)に示すように記述される。
Figure 0005155344
上記数式(1)において、磁極位置および回転速度に係る項を外乱として削除すると、モータ11の簡易的な電圧方程式は、下記数式(2)に示すように記述される。
Figure 0005155344
連続時間での簡易的な電圧方程式である上記数式(2)を、例えばオイラー近似を用いて離散化すると、離散時間状態方程式は、キャリア信号の1周期Tsおよび任意の自然数nにより、例えば下記数式(3)に示すように記述される。
なお、下記数式(3)において、任意の自然数nは、例えば任意の制御周期でのキャリア信号の山側のキャリア頂点に対応するインデックスである。
また、下記数式(3)において、3相−γδ変換部47から出力されるγ軸電流Iγ及びδ軸電流Iδの(n)番目の値つまりγ軸電流Iγ[n]及びδ軸電流Iδ[n]を、(n)番目でのγ軸モデル電流IγM[n]及びδ軸モデル電流IδM[n]とした。
Figure 0005155344
モータモデル電流演算部61は、電圧加算部45から出力されるγ軸電圧指令値Vγ及びδ軸電圧指令値Vδの(n−1)番目の値つまりγ軸電圧指令値Vγ[n−1]及びδ軸電圧指令値Vδ[n−1]と、3相−γδ変換部47から出力されるγ軸電流Iγ及びδ軸電流Iδの(n−1)番目の値つまりγ軸電流Iγ[n−1]及びδ軸電流Iδ[n−1]とを用いて、上記数式(3)に示すモータ11の所定モデルの電圧方程式により、(n)番目でのγ軸モデル電流IγM[n]及びδ軸モデル電流IδM[n]を算出し、各モデル電流IγM[n],IδM[n]を出力する。
また、d−q座標系でのモータ11の電圧方程式(つまり、連続時間での電圧方程式)は、磁極位置に係る拡張誘起電圧Eexを用いて記述すると、例えば下記数式(4)に示すように記述される。
下記数式(4)において、モータ11が低速域で回転(あるいは、停止)している状態であれば、角速度をゼロ(つまり、回転速度推定値ωe≒0)として無視することができる。さらに、各軸電流および各軸電圧指令値においてd軸をγ軸に置換し、かつq軸をδ軸に置換し、磁極位置および回転速度に係る拡張誘起電圧Eexの項を外乱として削除すると、モータ11の簡易的な電圧方程式は、上述した実施の形態と同様に、上記数式(2)に示すように記述される。
これにより、モータ11が低速域で回転あるいは停止している状態でのモータ11の所定モデルの電圧方程式は、上記数式(3)と同様に記述される離散時間状態方程式となる。
Figure 0005155344
γ軸電流偏差算出部62aおよびδ軸電流偏差算出部62bは、例えば下記数式(5)に示すように、同一タイミング、例えばキャリア信号の山側のキャリア頂点の(n)番目で、3相−γδ変換部47から出力されるγ軸電流Iγ[n]及びδ軸電流Iδ[n]と、モータモデル電流演算部61から出力される各モデル電流IγM[n],IδM[n]との各電流偏差ΔIγM[n],ΔIδM[n]を算出し、各電流偏差ΔIγM[n],ΔIδM[n]を出力する。
Figure 0005155344
出力切替部63は、ヒステリシス設定部64から出力される切替フラグのフラグ値に応じて、高速域磁極位置誤差演算部65または低速域磁極位置誤差演算部66の何れか一方を選択して、γ軸電流偏差算出部62aおよびδ軸電流偏差算出部62bから出力された各電流偏差ΔIγM[n],ΔIδM[n]の出力先を切り替える。
出力切替部63は、切替フラグのフラグ値が「1」である場合には、高速域磁極位置誤差演算部65に各電流偏差ΔIγM[n],ΔIδM[n]を出力し、一方、切替フラグのフラグ値が「0」である場合には、低速域磁極位置誤差演算部66に各電流偏差ΔIγM[n],ΔIδM[n]を出力する。
ヒステリシス設定部64は、モータ11の回転速度(例えば、回転速度−磁極位置演算部49から出力されるモータ11の回転速度推定値ωe)に応じてヒステリシスを設けて、切替フラグのフラグ値を「0」または「1」に設定する。
ヒステリシス設定部64は、切替フラグのフラグ値をモータ11の回転速度が所定値以上であるか否かに応じて切り替えるときことに対してヒステリシスを設定している。例えば図5に示すように、切替フラグのフラグ値が「0」に設定されている状態でモータ11の回転速度が上昇する場合には、所定の第1回転速度ωaよりも大きい所定の第2回転速度ωb(>ωa)以上になったときに切替フラグのフラグ値を「0」から「1」へと切り替える。一方、切替フラグのフラグ値が「1」に設定されている状態でモータ11の回転速度が低下する場合には、所定の第1回転速度ωa(<ωb)以上になったときに切替フラグのフラグ値を「1」から「0」へと切り替える。
高速域磁極位置誤差演算部65は、例えば下記数式(6)により位相差Δθeを算出する。
この下記数式(6)は、上記数式(1)に示すモータ11の電圧方程式を離散化して得た離散時間状態方程式と、上記数式(3)に示すモータ11の所定モデルの電圧方程式(つまり、モータ11の簡易的な電圧方程式である上記数式(2)を離散化して得た離散時間状態方程式)とに基づき算出される。
下記数式(6)によれば、例えば、γ軸電流偏差算出部62aおよびδ軸電流偏差算出部62bから出力される(n)番目での各電流偏差ΔIγM[n],ΔIδM[n]と、(n−1)番目での回転速度推定値ωe[n−1]と、3相−γδ変換部47から出力される(n−1)番目でのγ軸電流Iγ[n−1]及びδ軸電流Iδ[n−1]と、予め既知とされるモータ11の電気的回路定数(つまり、巻線抵抗Rおよび各インダクタンスLd,Lqおよびキャリア信号の1周期Ts)とにより、(n)番目での位相差Δθeが算出される。
Figure 0005155344
低速域磁極位置誤差演算部66は、例えば下記数式(7)により位相差Δθeを算出する。
この下記数式(7)は、上記数式(3)に示すモータ11の所定モデルの電圧方程式(つまり、モータ11の簡易的な電圧方程式である上記数式(2)を離散化して得た離散時間状態方程式)と、上記数式(4)に基づく離散時間状態方程式と、γ−δ座標系でのモータ11の電圧方程式(下記数式(8))に基づく拡張誘起電圧Eexの数式とに基づき算出される。
なお、上記数式(4)に基づく離散時間状態方程式は、上記数式(4)において、モータ11が低速域で回転(あるいは、停止)している状態であるとして、角速度をゼロ(つまり、回転速度推定値ωe≒0)として無視し、さらに、各軸電流および各軸電圧指令値でのd軸をγ軸に置換し、かつq軸をδ軸に置換して得たモータ11の電圧方程式を、離散化して得た離散時間状態方程式である。
また、下記数式(8)の電圧方程式に基づく拡張誘起電圧Eexの数式は、下記数式(8)に対して、上記数式(4)において各軸電流および各軸電圧指令値でのd軸をγ軸に置換し、かつq軸をδ軸に置換して得た電圧方程式を適用し、さらに、モータ11が低速域で回転(あるいは、停止)している状態であるとして、角速度をゼロ(つまり、回転速度推定値ωe≒0)として無視して得た数式である。
Figure 0005155344
Figure 0005155344
また、低速域磁極位置誤差演算部66は、例えば位相差Δθeが小さい(Δθe≒0)場合には、上記数式(7)に係る近似式である下記数式(9)により位相差Δθeを算出してもよい。
Figure 0005155344
回転速度−磁極位置演算部49は、高速域磁極位置誤差演算部65または低速域磁極位置誤差演算部66から出力される位相差Δθeに基づき、PLL(Phase-locked loop)による位相同期処理をおこなう。
回転速度−磁極位置演算部49は、例えば図6に示すように、位相同期部71と、積分演算部72とを備えて構成されている。
位相同期部71は、位相同期処理として、例えばPI(比例・積分)動作を実行し、例えば下記数式(10)に示すように記述される伝達関数Ge(s)に基づき、比例ゲインKpおよび積分ゲインKiにより、例えば下記数式(11)に示すようにして、位相差Δθeから回転速度推定値ωeを演算して、演算結果を出力する。
Figure 0005155344
Figure 0005155344
積分演算部72は、例えば下記数式(12)に示すように、位相同期部71から出力される回転速度推定値ωeを積分して磁極位置推定値θeを演算し、演算結果を出力する。
Figure 0005155344
なお、回転速度−磁極位置演算部49は、PI動作に限定されず、例えば図7および下記数式(13)に示す第1変形例のように、位相差Δθeを入力値とする同一次元オブザーバによる追従演算処理を実行して、回転速度推定値ωe及び磁極位置推定値θeを演算してもよい。
Figure 0005155344
電気角−機械角変換部50は、モータ11の極対数qに応じて、回転速度−磁極位置演算部49から出力される回転速度推定値ωeを回転速度ωr(機械角)に変換し、回転速度ωr(機械角)を出力する。
上述したように、本実施形態による電動機の磁極位置推定装置10によれば、モータモデル電流演算部61は、モータ11の回転速度に応じて互いに異なる演算により位相差Δθeを算出する2つの高速域磁極位置誤差演算部65と低速域磁極位置誤差演算部66とに対して、共通化された状態方程式により各モデル電流IγM,IδMを算出する。また、回転速度−磁極位置演算部49は、2つの2つの高速域磁極位置誤差演算部65と低速域磁極位置誤差演算部66とに対して、共通化された演算によってモータ11の磁極位置を算出する。
これにより、例えば高速域磁極位置誤差演算部65と低速域磁極位置誤差演算部66とに対して、共通化されていない演算によって各モデル電流IγM,IδMおよび位相差Δθeを算出する場合に比べて、装置構成を簡略化することができ、演算負荷を低減することができる。さらに、モータ11の回転速度が所定値未満であるか否かに応じて、磁極位置誤差演算部65による位相差Δθeの算出と、高速域低速域磁極位置誤差演算部66による位相差Δθeの算出との、何れか一方を選択して切り替えるときの安定性を向上させることができる。
さらに、モータ11の回転速度が所定値未満であるか否かに応じて、低速域磁極位置誤差演算部66による位相差Δθeの算出と、高速域磁極位置誤差演算部65による位相差Δθeの算出との、何れか一方を選択して切り替えることにより、適切にモータ11の磁極位置(磁極位置推定値θe)を算出することができる。低速域磁極位置誤差演算部66によれば、モータ11の停止時であっても磁極位置を算出することができる。高速域磁極位置誤差演算部65によれば、所定のγ軸検出用電圧Vhγおよびδ軸検出用電圧Vhδを印加する必要無しに磁極位置を算出することができる。
さらに、モータ11の回転速度が所定値未満であるか否かに応じて、低速域磁極位置誤差演算部66による位相差Δθeの算出と、高速域磁極位置誤差演算部65による位相差Δθeの算出との、何れか一方を選択して切り替えることに対してヒステリシスを設けることで、切替時にチャタリングが発生することを防止することができる。
なお、上述した実施の形態においては、各相電流センサ32の代わりに、インバータ13のブリッジ回路13aとバッテリ12の負極側端子あるいは正極側端子との間においてインバータ13のブリッジ回路13aの直流側電流Idcを検出する直流側電流センサを設けてもよい。この第2変形例では、直流側電流センサから出力される検出信号と、PWM信号生成部23からインバータ13に入力されるゲート信号とに基づき各相電流を推定し、各相電流の推定値を3相−γδ変換部47に入力する。なお、この場合、上記数式(3)に示す電圧方程式は、各相電流の推定値が有するタイミングと同一のタイミングで同期化されていればよい。
10 電動機の磁極位置推定装置
11 モータ
13 インバータ
22 制御装置
32 相電流センサ(電流検出手段)
45 電圧加算部(指令電圧出力手段)
49 回転速度−磁極位置演算部(磁極位置演算手段)
61 モータモデル電流演算部(モデル電流算出手段)
62a γ軸電流偏差算出部(電流差算出手段)
62b δ軸電流偏差算出部(電流差算出手段)
65 高速域磁極位置誤差演算部(第2の磁極位置差算出手段)
66 低速域磁極位置誤差演算部(第1の磁極位置差算出手段)

Claims (2)

  1. 電動機への通電に対する指令電圧を出力する指令電圧出力手段と、
    前記電動機に通電される実電流を検出して検出結果を出力する電流検出手段と、
    前記指令電圧と前記実電流とに基づいて、前記電動機の所定モデルの電圧方程式に応じたモデル電流を算出して、算出結果を出力するモデル電流算出手段と、
    前記モデル電流と前記実電流との差である電流差を算出する電流差算出手段と、
    前記電流差に基づいて、前記電動機の実際の磁極位置と推定又は指定による磁極位置との差である磁極位置差を、前記電動機の回転速度が所定値未満であるか否かに応じて、互いに異なる演算により算出する第1及び第2の磁極位置差算出手段と、
    前記第1の磁極位置差算出手段又は前記第2の磁極位置差算出手段により算出された前記磁極位置差に基づいて、前記電動機の磁極位置を演算する磁極位置演算手段と
    を備え
    前記第1の磁極位置差算出手段は、前記電動機の回転速度が所定値未満である場合に、前記電動機のインダクタンスが前記磁極位置に応じて変化することに基づいて、
    dq座標系での前記電動機の電圧方程式においてd軸をγ軸に置換し、かつq軸をδ軸に置換し、かつ回転速度をゼロとしてから該電圧方程式を離散化して得られる離散時間状態方程式と、該離散時間状態方程式を前記dq座標系に対して位相差を有するγδ座標系での前記電動機の電圧方程式に適用してから回転速度をゼロとして得られる式とに基づき、今回の制御サイクルでの前記電流差および前記実電流と前回の制御サイクルでの前記実電流とを用いて前記磁極位置差を算出し、
    前記第2の磁極位置差算出手段は、前記電動機の回転速度が所定値以上である場合に、前記電動機の誘起電圧が前記磁極位置に応じて変化することに基づいて、
    前記γδ座標系での前記電動機の電圧方程式において磁極位置および回転速度に係る項を外乱として削除する前後において該電圧方程式を離散化して得られる2つの離散時間状態方程式に基づき、今回の制御サイクルでの前記電流差と前回の制御サイクルでの前記実電流および前記回転速度とを用いて前記磁極位置差を算出することを特徴とする電動機の磁極位置推定装置。
  2. 前記第1の磁極位置差算出手段は、前記磁極位置差が小さいこと基づいて得られる近似式により前記磁極位置差を算出することを特徴とする請求項1に記載の電動機の磁極位置推定装置。
JP2010007019A 2010-01-15 2010-01-15 電動機の磁極位置推定装置 Expired - Fee Related JP5155344B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2010007019A JP5155344B2 (ja) 2010-01-15 2010-01-15 電動機の磁極位置推定装置
US12/968,296 US8384323B2 (en) 2010-01-15 2010-12-15 Motor magnetic-pole-position estimating apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010007019A JP5155344B2 (ja) 2010-01-15 2010-01-15 電動機の磁極位置推定装置

Publications (2)

Publication Number Publication Date
JP2011147287A JP2011147287A (ja) 2011-07-28
JP5155344B2 true JP5155344B2 (ja) 2013-03-06

Family

ID=44277139

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010007019A Expired - Fee Related JP5155344B2 (ja) 2010-01-15 2010-01-15 電動機の磁極位置推定装置

Country Status (2)

Country Link
US (1) US8384323B2 (ja)
JP (1) JP5155344B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180120831A (ko) * 2017-04-27 2018-11-07 현대모비스 주식회사 전동식 파워 스티어링 시스템의 모터 제어 장치 및 방법

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8264860B2 (en) * 2009-08-10 2012-09-11 Emerson Climate Technologies, Inc. System and method for power factor correction frequency tracking and reference generation
US8358098B2 (en) * 2009-08-10 2013-01-22 Emerson Climate Technologies, Inc. System and method for power factor correction
US8344706B2 (en) * 2009-08-10 2013-01-01 Emerson Climate Technologies, Inc. System and method for rejecting DC current in power factor correction systems
US8264192B2 (en) * 2009-08-10 2012-09-11 Emerson Climate Technologies, Inc. Controller and method for transitioning between control angles
US8698433B2 (en) 2009-08-10 2014-04-15 Emerson Climate Technologies, Inc. Controller and method for minimizing phase advance current
US8406021B2 (en) 2009-08-10 2013-03-26 Emerson Climate Technologies, Inc. System and method for reducing line current distortion
US8508166B2 (en) 2009-08-10 2013-08-13 Emerson Climate Technologies, Inc. Power factor correction with variable bus voltage
US8476873B2 (en) * 2009-08-10 2013-07-02 Emerson Climate Technologies, Inc. System and method for current balancing
US8493014B2 (en) 2009-08-10 2013-07-23 Emerson Climate Technologies, Inc. Controller and method for estimating, managing, and diagnosing motor parameters
JP5413400B2 (ja) * 2011-04-20 2014-02-12 株式会社安川電機 交流電動機の制御装置
KR101923731B1 (ko) 2012-02-03 2018-11-29 한국전자통신연구원 Bldc 모터의 구동 장치
US9634593B2 (en) 2012-04-26 2017-04-25 Emerson Climate Technologies, Inc. System and method for permanent magnet motor control
EP2883302B1 (en) 2012-08-10 2020-09-30 Emerson Climate Technologies, Inc. Motor drive control using pulse-width modulation pulse skipping
GB2507322B (en) 2012-10-26 2020-02-12 Control Techniques Ltd Rotor control
CN103684137B (zh) * 2013-11-21 2016-08-17 南京航空航天大学 一种基于串联电感斜率阈值的电励磁双凸极电机低速运行无位置法
FR3016256B1 (fr) * 2014-01-07 2016-01-22 Leroy Somer Moteurs Procede pour determiner la polarite d'un pole de rotor de machine electrique tournante
KR101628145B1 (ko) * 2014-06-16 2016-06-09 현대자동차 주식회사 전동기의 센서리스 제어 방법 및 시스템
US9705437B2 (en) 2014-09-24 2017-07-11 Texas Instruments Incorporated Angular position estimation for PM motors
JP6490540B2 (ja) * 2015-08-25 2019-03-27 株式会社東芝 回転位置検出装置,空気調和機及び回転位置検出方法
JP6667407B2 (ja) * 2016-09-12 2020-03-18 ルネサスエレクトロニクス株式会社 制御装置
JP7052255B2 (ja) * 2017-08-25 2022-04-12 コニカミノルタ株式会社 画像形成装置
AU2018425573B2 (en) * 2018-05-30 2021-12-02 Mitsubishi Electric Corporation Permanent-magnet synchronous motor and ventilation blower
JP7099226B2 (ja) * 2018-09-26 2022-07-12 株式会社アドヴィックス モータ制御装置
JP7154987B2 (ja) * 2018-12-11 2022-10-18 株式会社東芝 永久磁石同期電動機の制御装置,マイクロコンピュータ,電動機システム及び永久磁石同期電動機の運転方法
CN110912489A (zh) * 2019-12-13 2020-03-24 西安锐驰电器有限公司 一种电机转速的检测方法
JP6977837B1 (ja) * 2020-09-03 2021-12-08 株式会社明電舎 電力変換システムおよび電力変換システムの制御方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5463299A (en) * 1989-06-07 1995-10-31 Hitachi, Ltd. Current controller for controlling a current flowing in a load using a PWM inverter and method used thereby
JP3381408B2 (ja) * 1993-10-26 2003-02-24 トヨタ自動車株式会社 電気角検出装置およびこれを用いた同期モータの駆動装置
JP2000023498A (ja) * 1998-07-03 2000-01-21 Yaskawa Electric Corp Ipmモータの制御方法及び制御装置
JP3411878B2 (ja) * 2000-03-06 2003-06-03 株式会社日立製作所 同期モータの回転子位置推定方法、位置センサレス制御方法及び制御装置
JP3888082B2 (ja) * 2001-06-08 2007-02-28 株式会社豊田自動織機 モータ装置およびその制御方法
JP4681453B2 (ja) * 2003-09-02 2011-05-11 日本精工株式会社 電動パワーステアリング装置の制御装置
JP4198162B2 (ja) * 2006-04-07 2008-12-17 三洋電機株式会社 モータ制御装置
JP4895703B2 (ja) 2006-06-28 2012-03-14 三洋電機株式会社 モータ制御装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180120831A (ko) * 2017-04-27 2018-11-07 현대모비스 주식회사 전동식 파워 스티어링 시스템의 모터 제어 장치 및 방법
KR102205254B1 (ko) 2017-04-27 2021-01-21 현대모비스 주식회사 전동식 파워 스티어링 시스템의 모터 제어 장치 및 방법

Also Published As

Publication number Publication date
JP2011147287A (ja) 2011-07-28
US20110175560A1 (en) 2011-07-21
US8384323B2 (en) 2013-02-26

Similar Documents

Publication Publication Date Title
JP5155344B2 (ja) 電動機の磁極位置推定装置
JP4674525B2 (ja) 磁極位置推定方法及びモータ制御装置
JP4746667B2 (ja) 電動機の相電流推定装置および電動機の磁極位置推定装置
JP5081131B2 (ja) 電動機の相電流推定装置および電動機の磁極位置推定装置
JP4168730B2 (ja) 3相交流電動機の制御装置
JP6183291B2 (ja) 同期モータの制御装置
US9112436B2 (en) System for controlling controlled variable of rotary machine
JP4926492B2 (ja) モータ制御装置
JP2007097263A (ja) 同期モータの磁極位置推定方法
JP5485232B2 (ja) スイッチング回路の制御装置
JP5709932B2 (ja) 同期機制御装置
JP2004048958A (ja) Dcブラシレスモータの制御装置
JP2015149875A (ja) 電動機制御装置
JP5165545B2 (ja) 電動機の磁極位置推定装置
JP4781933B2 (ja) 電動機の制御装置
JP2002051580A (ja) 同期モータの位置センサレス制御方法および位置センサレス制御装置。
JP2013146155A (ja) 巻線温度推定装置及び巻線温度推定方法
JP5186352B2 (ja) 電動機の磁極位置推定装置
KR102409792B1 (ko) 영구 자석 동기 전동기의 제어 장치, 마이크로 컴퓨터, 전동기 시스템 및 영구 자석 동기 전동기의 운전 방법
JP2010028981A (ja) 同期モータの回転子位置推定方法および同期モータの制御装置
JP2010136586A (ja) 電動機の磁極位置推定装置
JP2010136585A (ja) 電動機の制御装置
JP5456873B1 (ja) 同期機制御装置
JP2011147286A (ja) 電動機の磁極位置推定装置
JP2010136582A (ja) 電動機の磁極位置推定装置

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20111215

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20111220

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120215

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: 20121113

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121206

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

Free format text: PAYMENT UNTIL: 20151214

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees