JP5534991B2 - Control device for synchronous motor - Google Patents

Control device for synchronous motor Download PDF

Info

Publication number
JP5534991B2
JP5534991B2 JP2010162169A JP2010162169A JP5534991B2 JP 5534991 B2 JP5534991 B2 JP 5534991B2 JP 2010162169 A JP2010162169 A JP 2010162169A JP 2010162169 A JP2010162169 A JP 2010162169A JP 5534991 B2 JP5534991 B2 JP 5534991B2
Authority
JP
Japan
Prior art keywords
current
unit
speed
phase
interference
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.)
Active
Application number
JP2010162169A
Other languages
Japanese (ja)
Other versions
JP2012023937A (en
Inventor
正樹 小原
Original Assignee
株式会社竹中製作所
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 株式会社竹中製作所 filed Critical 株式会社竹中製作所
Priority to JP2010162169A priority Critical patent/JP5534991B2/en
Publication of JP2012023937A publication Critical patent/JP2012023937A/en
Application granted granted Critical
Publication of JP5534991B2 publication Critical patent/JP5534991B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Ac Motors In General (AREA)

Description

本発明は、同期電動機の回転速度を制御するための制御装置に関する。 The present invention relates to a control equipment for controlling the rotational speed of the synchronous motor.

近年において、永久磁石を用いた同期電動機(永久磁石同期電動機、PMSM)は、その性能向上が著しく、小形かつ高効率であるため、家電製品を含めたあらゆる産業分野で用いられている。   In recent years, synchronous motors using permanent magnets (permanent magnet synchronous motors, PMSM) have been remarkably improved in performance, are small and highly efficient, and are used in all industrial fields including home appliances.

一般的なPMSMは、回転子に永久磁石を設け、固定子に電機子捲線を設けた回転界磁形の構造である。PMSMの制御には、回転子の回転角度位置を検出することが必要であるが、そのために位置センサを用いた場合には、小型化、耐ノイズ性、および価格などの点で不利になるため、位置センサを用いないセンサレス制御が行われている。   A general PMSM has a rotating field structure in which a permanent magnet is provided on a rotor and an armature winding is provided on a stator. For PMSM control, it is necessary to detect the rotational angle position of the rotor. However, if a position sensor is used for this purpose, it is disadvantageous in terms of downsizing, noise resistance, and cost. Sensorless control without using a position sensor is performed.

PMSMのセンサレス制御のために種々の方法が提案されている。   Various methods have been proposed for sensorless control of PMSM.

つまり、非特許文献1には、PMSMの固定子座標系で電動機の電流と回転磁束を推定する同一次元オブサーバを構成し、推定した電流値と実際の電動機に流れる電流値との差から回転速度を求める手法が開示されている。   That is, in Non-Patent Document 1, a PMSM stator coordinate system is used to construct a same-dimensional observer that estimates the current and rotating magnetic flux of the motor, and rotation is performed based on the difference between the estimated current value and the current value flowing through the actual motor. A technique for determining speed is disclosed.

また、非特許文献2には、PMSMの回転子が作る永久磁石の磁極軸を、d軸とd軸から電気角でπ/2[rad]進んだq軸とによるd−q軸座標系で電動機の電流と回転磁束を推定する同一次元オブサーバを構成し、推定した電流値と実際の電動機に流れる電流値との差から回転速度を求める手法が開示されている。   Non-Patent Document 2 discloses a dq-axis coordinate system in which the magnetic pole axis of a permanent magnet made by a PMSM rotor is d-axis and a q-axis advanced by π / 2 [rad] in electrical angle from the d-axis. A technique is disclosed in which a same-dimensional observer that estimates the electric current of the electric motor and the rotating magnetic flux is configured, and the rotational speed is obtained from the difference between the estimated current value and the current value flowing through the actual electric motor.

また、特許文献1には、PMSMのモデルに基づいて電機子電流を推定するオブザーバを用いた制御装置が開示されている。   Patent Document 1 discloses a control device using an observer that estimates an armature current based on a PMSM model.

特開2009−95135JP 2009-95135 A

楊耕・富岡理知子・中野求・金東海著、「適応オブザーバによるブラシレスDCモータの位置センサレス制御」、電気学会論文集D113巻5号1993年Saeko, Ritsuko Tomioka, Toru Nakano and Tokai Kin, “Position Sensorless Control of Brushless DC Motor with Adaptive Observer”, IEEJ Proceedings Vol. 113, No. 5, 1993 金原義彦著、「回転座標上の適応オブザーバを用いたPM電動機の位置センサレス制御」電気学会論文集D123巻5号2003年Yoshihiko Kanehara, “Position Sensorless Control of PM Motor Using Adaptive Observer on Rotating Coordinates” IEEJ Proceedings, Vol. 123, No. 5, 2003

しかし、上に述べた従来の方法においては、同期電動機の電流および回転磁束を推定するために推定変数が4次となる同一次元のオブサ−バを構成し、しかも推定した電流値と実際の電動機に流れる電流値との差から回転速度を求めているため、全体の演算量が増大するという問題がある。そのため、演算器の処理速度を高速にする必要があり、また記憶装置の記憶容量を大きくする必要があった。   However, in the conventional method described above, in order to estimate the current and rotating magnetic flux of the synchronous motor, an observer of the same dimension in which the estimation variable is fourth order is configured, and the estimated current value and the actual motor Since the rotation speed is obtained from the difference from the current value flowing through the current, there is a problem that the total amount of calculation increases. Therefore, it is necessary to increase the processing speed of the arithmetic unit and to increase the storage capacity of the storage device.

また、オブザーバを適切に動作させるためには、オブザーバにより推定した電動機の電流と実際の電流値との差分にある一定の値を掛けた量を調整し、しかも回転速度を求める演算においても適切に動作させるために演算値の調整が必要であり、調整に時間がかかるという問題があった。   In addition, in order to operate the observer appropriately, the amount obtained by multiplying the difference between the electric current of the motor estimated by the observer and the actual current value by a certain value is adjusted, and the calculation for calculating the rotation speed is also appropriate. In order to make it operate | move, adjustment of the calculated value was required and there existed a problem that adjustment took time.

本発明は、同期電動機の回転速度をセンサレスで制御するに際し、従来のようなオブザーバを用いることなく、演算量の増大を抑制することのできる装置を提供することを目的とする。 The present invention, when controlling the rotational speed of the synchronous motor without using a sensor, without using a conventional such observer, and an object thereof is to provide a equipment which can suppress the increase in the calculation amount.

本発明に係る装置は、同期電動機の回転速度を制御するための制御装置であって、速度指令に基づいて速度調節を行う速度調節部と、電流調節を行う電流調節部と、前記電流調節部から出力される電圧指令と回転速度の推定値とに基づいて干渉分を除去するための非干渉制御を行う非干渉制御部と、前記非干渉制御部の出力する電圧を座標変換する2相/3相座標変換部と、前記2相/3相座標変換部の出力に基づいて電機子捲線に電流を流して駆動するための電力変換部と、前記電機子捲線に流れる電流を検出する電流検出部と、前記電流検出部により検出された電流を座標変換する3相/2相座標変換部と、前記電圧指令に基づいて干渉分を除去した制御モデルにおける電流の目標値を算出する電流目標算出部と、前記電流検出部により検出された電流と前記電流目標算出部により算出された電流の目標値との差が漸近的に零になるように回転速度の推定値および磁極位置の推定値を求める速度位置推定部と、を有し、前記制御モデルは、その電圧電流方程式が、次の(1)式で示される。

Figure 0005534991
An apparatus according to the present invention is a control device for controlling the rotational speed of a synchronous motor, and includes a speed adjusting unit that adjusts a speed based on a speed command, a current adjusting unit that adjusts a current, and the current adjusting unit. A non-interference control unit that performs non-interference control for removing interference based on the voltage command output from the rotation speed estimation value and a two-phase / coordinate conversion of the voltage output from the non-interference control unit A three-phase coordinate conversion unit, a power conversion unit for driving the current through the armature winding based on the output of the two-phase / three-phase coordinate conversion unit, and current detection for detecting the current flowing through the armature winding Unit, a three-phase / two-phase coordinate conversion unit that performs coordinate conversion of the current detected by the current detection unit, and a current target calculation that calculates a target value of current in a control model that eliminates interference based on the voltage command And the current detection unit A speed position estimating unit for obtaining an estimated value of the rotational speed and an estimated value of the magnetic pole position so that a difference between the generated current and the target value of the current calculated by the current target calculating unit is asymptotically zero. Yes, and the control model, the voltage-current equation is expressed by the following equation (1).
Figure 0005534991

前記速度位置推定部は、例えば、The speed position estimation unit is, for example,
次の(2)式に基づいて、前記回転速度の推定値および磁極位置の推定値を求める。Based on the following equation (2), an estimated value of the rotational speed and an estimated value of the magnetic pole position are obtained.

Figure 0005534991
Figure 0005534991

本発明は、同期電動機の回転速度をセンサレスで制御するに際し、従来のようなオブザーバを用いることなく、演算量の増大を抑制することのできる装置を提供することができる。
The present invention, when controlling the rotational speed of the synchronous motor without using a sensor, without using the conventional kind of observer may provide equipment which can suppress the increase in the calculation amount.

本実施形態におけるPMSMの制御装置の構成の例を示すブロック図である。It is a block diagram which shows the example of a structure of the control apparatus of PMSM in this embodiment. PMSMの等価回路である。It is an equivalent circuit of PMSM. 本実施形態における制御装置の概略の制御内容を示すフローチャートである。It is a flowchart which shows the outline control content of the control apparatus in this embodiment. センサレスでないPMSMの制御装置の構成の例を示すブロック図である。It is a block diagram which shows the example of a structure of the control apparatus of PMSM which is not sensorless. 非線形フィードバック系を示すブロック図である。It is a block diagram which shows a nonlinear feedback system.

図1には、PMSMのセンサレスによる非干渉ベクトル制御のブロック図が示されている。また、図2には、PMSMの等価回路が示されている。   FIG. 1 is a block diagram of PMSM sensorless non-interference vector control. FIG. 2 shows an equivalent circuit of PMSM.

図1において、PMSM(Permanent Magnet Synchronous Motor)は、回転子(ロータ)に永久磁石が設けられ、固定子(ステータ)に3相の電機子捲線が設けられた回転界磁形の永久磁石同期電動機である。PMSMは、速度指令(回転速度指令)ω*に応じて周波数制御され制御装置1から出力される3相交流電力Poutによって、その周波数に同期した回転速度ωで回転するよう駆動される。   In FIG. 1, a PMSM (Permanent Magnet Synchronous Motor) is a rotating field type permanent magnet synchronous motor in which a rotor (rotor) is provided with a permanent magnet and a stator (stator) is provided with a three-phase armature winding. It is. The PMSM is driven to rotate at a rotational speed ω synchronized with the frequency by the three-phase AC power Pout that is frequency-controlled according to the speed command (rotational speed command) ω * and is output from the control device 1.

つまり、PMSMの基本モデルは、図2に示されるように、永久磁石のN極の方向をd軸とし、これより電気角でπ/2進んだ方向をq軸とする。つまり、d軸、q軸はモデル軸である。U相の電機子捲線を基準とし、これに対するd軸の進み角をθとする。つまり、d軸−q軸座標系は、U相の電機子捲線を基準としてこれより角度θだけ進んだ位置にある。   That is, in the PMSM basic model, as shown in FIG. 2, the direction of the N pole of the permanent magnet is the d-axis, and the direction advanced by π / 2 from this is the q-axis. That is, the d axis and the q axis are model axes. The lead angle of the d-axis relative to the U-phase armature winding is defined as θ. That is, the d-axis-q-axis coordinate system is at a position advanced by an angle θ from the U-phase armature winding as a reference.

なお、θは、U相の電機子捲線に対する磁極の実際の角度位置(磁極位置)を示す。角度位置θは、回転速度ωを積分することによっても求めることができる。   Θ represents the actual angular position (magnetic pole position) of the magnetic pole with respect to the U-phase armature winding. The angular position θ can also be obtained by integrating the rotational speed ω.

また、U相の電機子捲線を基準とし、推定された磁極位置の角度θ^に対応してγ軸が決定される。γ軸よりも電気角でπ/2進んだ位置がδ軸である。つまり、γ軸−δ軸座標系は、U相の電機子捲線を基準としてこれより角度θ^だけ進んだ位置にある。また、d軸−q軸座標系とγ軸−δ軸座標系とは極接近しており、θとθ^とはほぼ等しく、この平衡点近傍において適応制御が行われる。   The γ-axis is determined in correspondence with the estimated angle θ ^ of the magnetic pole position with reference to the U-phase armature winding. A position advanced by π / 2 in electrical angle from the γ axis is the δ axis. That is, the γ-axis δ-axis coordinate system is at a position advanced by an angle θ ^ with respect to the U-phase armature winding. Further, the d-axis-q-axis coordinate system and the γ-axis-δ-axis coordinate system are extremely close to each other, and θ and θ ^ are substantially equal, and adaptive control is performed in the vicinity of the equilibrium point.

なお、U相、V相、W相の各電機子捲線は、それぞれ、漏れインダクタンスl’および有効インダクタンスL’を含む自己インダクタンスL(Lu、Lv、Lw)、および捲線抵抗Raを持っている。   The U-phase, V-phase, and W-phase armature windings each have a self-inductance L (Lu, Lv, Lw) including a leakage inductance l 'and an effective inductance L', and a winding resistance Ra.

本明細書において、記号「*」は、各パラメータにおける指令または指令値などを示す場合に、記号「^」は、各パラメータの推定値または算出値などを示す場合に、それぞれ用いられる。   In this specification, the symbol “*” is used to indicate a command or a command value for each parameter, and the symbol “^” is used to indicate an estimated value or a calculated value for each parameter.

さて、制御装置1は、速度調節部11、電流調節部12(12a,12b)、非干渉制御部13、2相/3相座標変換部14、電力変換部15、電流検出部16、3相/2相座標変換部17、電流目標算出部18、および速度位置推定部19を備える。   The control device 1 includes a speed adjustment unit 11, a current adjustment unit 12 (12a, 12b), a non-interference control unit 13, a two-phase / three-phase coordinate conversion unit 14, a power conversion unit 15, a current detection unit 16, and a three-phase. / 2 phase coordinate conversion unit 17, current target calculation unit 18, and speed position estimation unit 19 are provided.

速度調節部11は、速度指令(速度指令値)ω*に基づいた速度調節を行う。つまり、外部から入力された速度指令値ω*と、速度位置推定部19から出力された速度推定値ω^との差(ω*−ω^)に基づいて、δ軸電流指令値iδ*を演算する。   The speed adjustment unit 11 performs speed adjustment based on a speed command (speed command value) ω *. That is, based on the difference (ω * −ω ^) between the speed command value ω * input from the outside and the speed estimated value ω ^ output from the speed position estimation unit 19, the δ-axis current command value iδ * is Calculate.

つまり、速度調節部11は、速度指令値ω*と速度推定値ω^との差(ω*−ω^)を偏差として、比例・積分制御(PI制御)などのフィードバック制御によって速度指令値通りの回転速度(回転速度実際値)ωとなるように(偏差が零になるように)、δ軸電流指令値(トルク電流指令)iδ*を作り出す。   That is, the speed adjustment unit 11 uses the difference (ω * −ω ^) between the speed command value ω * and the estimated speed value ω ^ as a deviation, and performs feedback control such as proportional / integral control (PI control) according to the speed command value. The δ-axis current command value (torque current command) iδ * is generated so that the rotation speed (actual rotation speed actual value) ω becomes (so that the deviation becomes zero).

電流調節部12a,12bは、それぞれ、γ軸電流指令値iγ*またはδ軸電流指令値iδ*に基づいて電流調節を行う。つまり、電流調節部12aは、γ軸電流指令値iγ*とγ軸電流検出値iγとの差(iγ*−iγ)に基づいて、γ軸電圧指令値Vγ*を演算する。電流調節部12bは、δ軸電流指令値iδ*とδ軸電流検出値iδとの差(iδ*−iδ)に基づいて、δ軸電圧指令値Vδ*を演算する。   The current adjusters 12a and 12b perform current adjustment based on the γ-axis current command value iγ * or the δ-axis current command value iδ *, respectively. That is, the current adjustment unit 12a calculates the γ-axis voltage command value Vγ * based on the difference (iγ * −iγ) between the γ-axis current command value iγ * and the γ-axis current detection value iγ. The current adjusting unit 12b calculates the δ-axis voltage command value Vδ * based on the difference (iδ * −iδ) between the δ-axis current command value iδ * and the δ-axis current detection value iδ.

つまり、電流調節部12a,12bは、それぞれ、差(iγ*−iγ)、(iδ*−iδ)を偏差として、フィードバック制御により指令値通りの電流値になるように(偏差が零になるように)、γ軸またはδ軸の電圧指令Vγ*、Vδ*を作り出す。   In other words, the current adjusting units 12a and 12b use the differences (iγ * −iγ) and (iδ * −iδ) as deviations, respectively, so that the current values are equal to the command values by feedback control (the deviations become zero). ) To generate voltage commands Vγ * and Vδ * for the γ-axis or δ-axis.

非干渉制御部13は、電流調節部12から出力される電圧指令(γ軸電圧指令値Vγ*、δ軸電圧指令値Vδ*)と、回転速度の推定値(速度推定値ω^)とに基づいて、干渉分を除去するための非干渉制御を行う。つまり、非干渉制御部13は、γ軸電圧指令値Vγ*、δ軸電圧指令値Vδ*、γ軸電流検出値iγ、δ軸電流検出値iδ、および速度推定値ω^に基づいて、γ軸電圧指令値Vγ**およびδ軸電圧指令値Vδ**を演算する。非干渉制御部13の出力するγ軸電圧指令値Vγ**およびδ軸電圧指令値Vδ**は、入力されるγ軸電圧指令値Vγ*およびδ軸電圧指令値Vδ*に対して、干渉分が除去されている。   The non-interference control unit 13 uses the voltage command (γ-axis voltage command value Vγ *, δ-axis voltage command value Vδ *) output from the current adjustment unit 12 and the estimated rotational speed (speed estimated value ω ^). Based on this, non-interference control for removing the interference is performed. That is, the non-interference control unit 13 determines that the γ-axis voltage command value Vγ *, the δ-axis voltage command value Vδ *, the γ-axis current detection value iγ, the δ-axis current detection value iδ, and the speed estimation value ω ^ An axis voltage command value Vγ ** and a δ axis voltage command value Vδ ** are calculated. The γ-axis voltage command value Vγ ** and the δ-axis voltage command value Vδ ** output from the non-interference control unit 13 interfere with the input γ-axis voltage command value Vγ * and δ-axis voltage command value Vδ *. Minutes have been removed.

なお、非干渉制御部13において、干渉分とは、PMSMの回転子の回転にともなって、電機子捲線のインダクタンスLなどに誘起される電圧成分(逆起電力)ωL、−ωL、ωφmなどである。インダクタンスLおよび回転子磁束φmは予め求めておくことができ、これらと回転速度の推定値(速度推定値)ω^から干渉分を求めることができる。   In the non-interference control unit 13, the interference component is a voltage component (back electromotive force) ωL, −ωL, ωφm, etc. induced in the inductance L of the armature winding along with the rotation of the PMSM rotor. is there. The inductance L and the rotor magnetic flux φm can be obtained in advance, and the interference can be obtained from these and the estimated value (speed estimated value) ω ^ of the rotational speed.

なお、非干渉制御それ自体は従来から用いられている技術であり、例えば、杉本英彦編著「ACサーボシステムの理論と設計の実際」(総合電子出版社、1999/05/08発行、77−80頁)を参照することができる。   Non-interference control itself is a technique that has been used in the past. For example, Hidehiko Sugimoto, “Theory and Design of AC Servo Systems” (general electronic publisher, 1999/05/08, 77-80) Page).

2相/3相座標変換部14は、非干渉制御部13の出力する電圧を座標変換する。つまり、非干渉制御部13からのγ軸電圧指令値Vγ**およびδ軸電圧指令値Vδ**を、電力変換部15を制御するための、u、v、wの3相の電圧指令値に変換する。   The two-phase / three-phase coordinate conversion unit 14 performs coordinate conversion of the voltage output from the non-interference control unit 13. That is, the γ-axis voltage command value Vγ ** and the δ-axis voltage command value Vδ ** from the non-interference control unit 13 are used to control the power conversion unit 15 in three-phase voltage command values u, v, and w. Convert to

電力変換部15は、PWM制御部およびインバータなどを備える。電力変換部15は、2相/3相座標変換部14からの出力に基づいてPWM(Pulse Width Modulation)信号を生成し、PWM信号をインバータのスイッチング素子の各ゲートに印加する。インバータからは、周波数および振幅が制御された3相交流電力Poutが出力され、これによってPMSMのu、v、wの電機子捲線に3相電流が流れ、回転子がそれに同期して回転駆動される。   The power conversion unit 15 includes a PWM control unit and an inverter. The power conversion unit 15 generates a PWM (Pulse Width Modulation) signal based on the output from the 2-phase / 3-phase coordinate conversion unit 14, and applies the PWM signal to each gate of the switching element of the inverter. The inverter outputs a three-phase AC power Pout with a controlled frequency and amplitude. This causes a three-phase current to flow through the PMSM u, v, and w armature windings, and the rotor is driven to rotate in synchronization therewith. The

電流検出部16は、PMSMの電機子捲線に流れる電流を検出する。つまり、PMSMに流れている実際の電流(電流検出値)を検出する。図1においては、u、v、wの各相に流れる電流iu、iv、iwを検出しているが、いずれか2相のみを検出して残りの1相については演算で求めるようにしてもよい。   The current detector 16 detects a current flowing through the armature winding of the PMSM. That is, the actual current (current detection value) flowing through the PMSM is detected. In FIG. 1, the currents iu, iv, and iw flowing in the u, v, and w phases are detected. However, only one of the two phases may be detected, and the remaining one phase may be obtained by calculation. Good.

3相/2相座標変換部17は、電流検出部16により検出された電流を座標変換する。つまり、電流検出部16からの3相の電流iu、iv、iwを、γ−δ軸座標系におけるγ軸方向の電流iγおよびδ軸方向の電流iδに変換する。電流iγ、iδは、変換されたものではあるが、PMSMに実際に流れる電流に対応しているので、これも電流検出値である。したがって、「電流検出値iγ、iδ」と記載することがある。   The three-phase / two-phase coordinate conversion unit 17 performs coordinate conversion of the current detected by the current detection unit 16. That is, the three-phase currents iu, iv, iw from the current detector 16 are converted into a current iγ in the γ-axis direction and a current iδ in the δ-axis direction in the γ-δ axis coordinate system. Although the currents i γ and i δ are converted, they correspond to the current that actually flows through the PMSM, and are also current detection values. Therefore, it may be described as “current detection values iγ, iδ”.

電流目標算出部18は、電流調節部12から出力されるγ軸電圧指令値Vγ*およびδ軸電圧指令値Vδ*に基づいて、干渉分を除去した制御モデルにおける電流の目標値を算出する。つまり、電流目標算出部18は、γ軸電圧指令値Vγ*およびδ軸電圧指令値Vδ*に基づいて、PMSMの回転速度ωに比例した逆起電力の影響を取り除いた理想的なPMSMの電流回路における、d軸およびq軸の電流目標値であるid^、iq^を演算する。   Based on the γ-axis voltage command value Vγ * and the δ-axis voltage command value Vδ * output from the current adjustment unit 12, the current target calculation unit 18 calculates a target value of current in the control model from which the interference is removed. In other words, the current target calculation unit 18 eliminates the influence of the counter electromotive force proportional to the rotational speed ω of the PMSM based on the γ-axis voltage command value Vγ * and the δ-axis voltage command value Vδ *. In the circuit, id ^ and iq ^, which are d-axis and q-axis current target values, are calculated.

なお、電流目標値id^、iq^は、干渉分を取り除いた理想的な制御モデルにおける目標値であり、これを電流目標算出部18において演算によって推定することから、「電流推定値」ということができる。したがって、「電流推定値id^、iq^」と記載することがある。   The current target values id ^ and iq ^ are target values in an ideal control model from which interference has been removed, and are estimated by calculation in the current target calculation unit 18 and are therefore referred to as "current estimated values". Can do. Therefore, it may be described as “current estimation value id ^, iq ^”.

また、電流目標算出部18は、干渉分を取り除いた理想的な制御モデル(適応モデル)における適応制御のための目標値を求めるものであり、電流目標算出部18において適応モデルを定義するものであるので、電流目標算出部18を「適応モデル部18」ということもできる。   The current target calculation unit 18 obtains a target value for adaptive control in an ideal control model (adaptive model) from which interference is removed, and the current target calculation unit 18 defines an adaptive model. Therefore, the current target calculation unit 18 can also be referred to as the “adaptive model unit 18”.

速度位置推定部19は、電流検出部16により検出された電流iγ、iδと、電流目標算出部18により算出された電流の目標値id^、iq^との差(id^−iγ)、(iq^−iδ)が、いずれも漸近的に零になるように、回転速度の推定値(速度推定値)ω^、および磁極位置の推定値(磁極推定値)θ^を求める。   The speed position estimator 19 compares the difference between the currents iγ and iδ detected by the current detector 16 and the target values id ^ and iq ^ calculated by the current target calculator 18 (id ^ -iγ), ( The estimated value of rotational speed (speed estimated value) ω ^ and the estimated value of magnetic pole position (magnetic pole estimated value) θ ^ are determined so that both iq ^ -iδ) become asymptotically zero.

なお、γ−δ軸座標系では、速度位置推定部19によって推定されるPMSMの磁極軸がγ軸であり、γ軸から電気角でπ/2だけ進んだ軸がδ軸である。   In the γ-δ axis coordinate system, the magnetic pole axis of the PMSM estimated by the velocity position estimation unit 19 is the γ axis, and the axis advanced by π / 2 in electrical angle from the γ axis is the δ axis.

速度位置推定部19によって求められた速度推定値ω^は、速度調節部11におけるδ軸電流指令値iδ*の演算に用いられる。また、磁極位置の推定値θ^は、2相/3相座標変換部14における座標変換に用いられる。   The estimated speed value ω ^ obtained by the speed position estimating unit 19 is used in the calculation of the δ-axis current command value iδ * in the speed adjusting unit 11. The estimated value θ ^ of the magnetic pole position is used for coordinate conversion in the two-phase / three-phase coordinate conversion unit 14.

電流目標算出部18および速度位置推定部19は、電流検出値iγ、iδおよび電圧指令値Vγ*、Vδ*に基づいて、制御モデルにおける電流の誤差が漸近的に零になるように適応制御を行うものであり、これらは適応制御部TSを構成する。   The current target calculation unit 18 and the speed position estimation unit 19 perform adaptive control based on the detected current values iγ and iδ and the voltage command values Vγ * and Vδ * so that the current error in the control model becomes asymptotically zero. These are performed, and these constitute the adaptive control unit TS.

なお、適応制御部TSは、誤差が漸近的に零になるように適応制御を行うものであるが、誤差が漸近的に零になるようにとは、制御が乱れないよう、急激な変化が起こることをできるだけ避けるように制御が行われることを意味する。   The adaptive control unit TS performs adaptive control so that the error is asymptotically zero. However, when the error is asymptotically zero, there is a sudden change so that the control is not disturbed. It means that control is done to avoid as much as possible.

本実施形態の制御装置1では、適応モデルと呼ぶPMSMの簡単な電流回路モデルを用い、その適応モデルにおいて推定した電流値とPMSMに実際に流れる電流値との差から、回転速度と磁極位置とを推定する。従来に用いられていた位置センサ、速度センサによる検出値に代えて、推定された速度推定値ω^および磁極推定値θ^を用いて制御を行う。   In the control device 1 of the present embodiment, a simple current circuit model of PMSM called an adaptive model is used. From the difference between the current value estimated in the adaptive model and the current value actually flowing through the PMSM, the rotational speed and the magnetic pole position are calculated. Is estimated. Control is performed using the estimated speed estimated value ω ^ and magnetic pole estimated value θ ^ instead of the detection values obtained by the position sensor and speed sensor conventionally used.

そこで、電流目標算出部18で用いられる適応モデルは、その電圧電流方程式を次の(7)式で示すことができる。   Therefore, the adaptive model used in the current target calculation unit 18 can show the voltage-current equation by the following equation (7).

Figure 0005534991
Figure 0005534991

また、速度位置推定部19は、次の(8)式に基づいて、速度推定値ω^、磁極推定値θ^を求めることができる。   Moreover, the speed position estimation part 19 can obtain | require speed estimated value (omega) ^ and magnetic pole estimated value (theta) ^ based on following (8) Formula.

Figure 0005534991
Figure 0005534991

また、速度位置推定部19は、次の(9)式に基づいて、速度推定値ω^、磁極推定値θ^を求めることもできる。   Moreover, the speed position estimation part 19 can also obtain | require speed estimated value (omega) ^ and magnetic pole estimated value (theta) ^ based on following (9) Formula.

Figure 0005534991
Figure 0005534991

速度位置推定部19が上の(8)式または(9)式を用いる場合に、非干渉制御部13は、後に述べる(16)式に基づいて非干渉制御を行う。   When the speed position estimation unit 19 uses the above equation (8) or (9), the non-interference control unit 13 performs non-interference control based on equation (16) described later.

また、速度位置推定部19は、次の(10)式に基づいて、速度推定値ω^、磁極推定値θ^を求めることもできる。   Moreover, the speed position estimation part 19 can also obtain | require speed estimated value (omega) ^ and magnetic pole estimated value (theta) ^ based on following (10) Formula.

Figure 0005534991
Figure 0005534991

速度位置推定部19が上の(10)式を用いる場合に、非干渉制御部13は、次の(11)式に基づいて非干渉制御を行う。   When the speed position estimation unit 19 uses the above equation (10), the non-interference control unit 13 performs non-interference control based on the following equation (11).

Figure 0005534991
Figure 0005534991

本実施形態の制御装置1においては、図3のフローチャートに示すように、PMSMの電機子捲線に流れる電流iγ、iδを検出し(#11)、電圧指令Vγ*、Vδ*に基づいて、干渉分を除去した制御モデルにおける電流の目標値id^、iq^を算出し(#12)、算出した電流の目標値と電流の検出値iγ、iδとの差(id^−iγ)、(iq^−iδ)が漸近的に零になるように、回転速度の推定値ω^および磁極位置の推定値θ^を求め(#13)、求めた回転速度の推定値ω^および磁極位置の推定値θ^に基づいて、速度制御および座標変換を行う(#14)。   In the control device 1 of the present embodiment, as shown in the flowchart of FIG. 3, the currents iγ and iδ flowing through the armature winding of the PMSM are detected (# 11), and interference is performed based on the voltage commands Vγ * and Vδ *. The current target values id ^ and iq ^ in the control model with the minutes removed are calculated (# 12), and the difference between the calculated current target value and the detected current values iγ and iδ (id ^ -iγ), (iq The estimated rotational speed value ω ^ and the estimated magnetic pole position θ ^ are obtained so that (^ -iδ) becomes asymptotically zero (# 13), and the obtained estimated rotational speed value ω ^ and the estimated magnetic pole position are obtained. Based on the value θ ^, speed control and coordinate conversion are performed (# 14).

本実施形態の制御装置1では、適応モデルと呼ぶPMSMの簡単な電流回路モデルを用い、その適応モデルにおいて推定した電流値とPMSMに実際に流れる電流値との差から回転速度と磁極位置とを推定するだけであるから、演算量がそれほど増えることなく抑制することができ、しかも調整が簡単である。   In the control device 1 of this embodiment, a simple PMSM current circuit model called an adaptive model is used, and the rotational speed and the magnetic pole position are calculated from the difference between the current value estimated in the adaptive model and the current value actually flowing through the PMSM. Since only the estimation is performed, the calculation amount can be suppressed without increasing so much, and the adjustment is easy.

次に、制御装置1についてさらに詳しく説明する。   Next, the control device 1 will be described in more detail.

図4には、本実施形態の制御装置1との比較のために、電流目標算出部18および速度位置推定部19を用いることなく、位置センサを用いた、つまりセンサレスでない制御装置80のブロック図が示されている。   FIG. 4 is a block diagram of a control device 80 that uses a position sensor, that is, is not sensorless, without using the current target calculation unit 18 and the speed position estimation unit 19 for comparison with the control device 1 of the present embodiment. It is shown.

図4において、制御装置80は、速度調節部81、電流調節部82a,82b、非干渉制御部83、2相/3相座標変換部84、電力変換部85、電流検出部86、3相/2相座標変換部87、位置検出器88、および速度検出部89を備える。   4, the control device 80 includes a speed adjustment unit 81, current adjustment units 82a and 82b, a non-interference control unit 83, a two-phase / three-phase coordinate conversion unit 84, a power conversion unit 85, a current detection unit 86, and a three-phase / A two-phase coordinate converter 87, a position detector 88, and a speed detector 89 are provided.

PMSMの回転子の磁極軸をd軸として、d−q軸座標系での電圧電流方程式は、次の(12)式で示すことができる。   With the magnetic pole axis of the PMSM rotor as the d axis, the voltage / current equation in the dq axis coordinate system can be expressed by the following equation (12).

Figure 0005534991
Figure 0005534991

上の(12)式を、電圧値を入力変数とし、電流値を状態変数とした状態方程式は、次の(13)式で示すことができる。   The equation (12) above can be expressed by the following equation (13) with the voltage value as the input variable and the current value as the state variable.

Figure 0005534991
Figure 0005534991

すなわち、d軸、q軸の実際の電流(電流検出値)id、iqは、入力電圧であるd軸電圧値Vdおよびq軸電圧値Vqが変化しなくても、その電流が流れているときに回転子が回転すると、回転子の回転速度ωに比例した逆起電力によって影響を受ける。   That is, the actual currents (current detection values) id and iq of the d-axis and q-axis are when the current flows even if the d-axis voltage value Vd and the q-axis voltage value Vq that are input voltages do not change. When the rotor rotates, it is affected by the counter electromotive force proportional to the rotational speed ω of the rotor.

つまり、実際のd軸電流idおよびq軸電流iqは、回転子による逆起電力によって干渉を受けることになり、これがPMSMの電流応答の劣化の原因となる。このような干渉をなくすために、干渉分を除去する非干渉制御が従来より行われているのである。   That is, the actual d-axis current id and q-axis current iq are subject to interference by the counter electromotive force generated by the rotor, which causes deterioration of the current response of the PMSM. In order to eliminate such interference, non-interference control for removing the interference is conventionally performed.

そこで、非干渉制御部83では、次の(14)式に示す演算を行う。   Therefore, the non-interference control unit 83 performs the calculation shown in the following equation (14).

Figure 0005534991
Figure 0005534991

すなわち、非干渉制御部83において、回転速度ωに比例する正負の逆起電力の成分を持った電圧指令値を発生し、これによって逆起電力の影響を打ち消す。これにより非干渉制御が行われる。   That is, the non-interference control unit 83 generates a voltage command value having positive and negative back electromotive force components proportional to the rotational speed ω, thereby canceling the influence of the back electromotive force. Thereby, non-interference control is performed.

その結果、電流調節部82の出力する電圧指令Vq*、Vd*を入力変数とし、電流値id、iqを状態変数とした状態方程式は、次の(15)式で示すことができる。   As a result, a state equation in which the voltage commands Vq * and Vd * output from the current adjustment unit 82 are input variables and the current values id and iq are state variables can be expressed by the following equation (15).

Figure 0005534991
Figure 0005534991

上の(15)式で示されるように、回転速度ωに比例した逆起電力ωL、−ωL、ωφmが消去されており、d−q軸の非干渉制御が実現されている。   As shown by the above equation (15), the counter electromotive forces ωL, −ωL, and ωφm proportional to the rotational speed ω are eliminated, and the dq axis non-interference control is realized.

図1に戻って、本実施形態の制御装置1はセンサレスによる非干渉制御を行っているため、非干渉制御部13は、実際の回転速度ωの代わりに、速度位置推定部19から出力される回転速度の推定値ω^を入力として、次の(16)式の演算を実行する。   Returning to FIG. 1, since the control device 1 of the present embodiment performs non-interference control without a sensor, the non-interference control unit 13 is output from the speed position estimation unit 19 instead of the actual rotational speed ω. The following equation (16) is calculated with the estimated rotational speed value ω ^ as an input.

Figure 0005534991
Figure 0005534991

非干渉制御部13の出力するδ軸電圧指令値Vδ**、γ軸電圧指令値Vγ**が、2相/3相座標変換部14に入力され、2相/3相座標変換部14において、速度位置推定部19から出力される磁極位置の推定値θ^を用いて、3相の交流電圧指令値に変換される。2相/3相座標変換部14の出力する電圧指令値によって、電力変換部15は3相交流電力Poutを出力し、これによってPMSMのに実際の電圧が加えられる。   The δ-axis voltage command value Vδ ** and the γ-axis voltage command value Vγ ** output from the non-interference control unit 13 are input to the 2-phase / 3-phase coordinate conversion unit 14, and the 2-phase / 3-phase coordinate conversion unit 14 The magnetic pole position estimated value θ ^ output from the speed position estimating unit 19 is used to convert it into a three-phase AC voltage command value. The power converter 15 outputs the three-phase AC power Pout according to the voltage command value output from the two-phase / three-phase coordinate converter 14, whereby an actual voltage is applied to the PMSM.

PMSMの電機子捲線に実際に流れる電流が、電流検出部16によって検出され、3相/2相座標変換部17に入力される。3相/2相座標変換部17において、速度位置推定部19から出力される磁極位置の推定値θ^を用いて、γ軸、δ軸の2相の直流電流の検出値iγ、iδに変換される。   The current that actually flows through the armature winding of the PMSM is detected by the current detection unit 16 and input to the three-phase / two-phase coordinate conversion unit 17. In the three-phase / two-phase coordinate conversion unit 17, the estimated value θ ^ of the magnetic pole position output from the velocity position estimation unit 19 is used to convert the detected values iγ and iδ of the two-phase DC currents of the γ-axis and δ-axis. Is done.

電流目標算出部18は、電流調節部12から出力される、γ軸電圧指令Vγ*、δ軸電圧指令Vδ*を入力値として、回転速度ωに比例した逆起電力の影響を取り除いた理想的なPMSMの電流回路の状態方程式である上の(15)式におけるγ軸、δ軸上の次の(17)式を演算する。   The current target calculation unit 18 uses the γ-axis voltage command Vγ * and the δ-axis voltage command Vδ * output from the current adjustment unit 12 as input values and removes the influence of the counter electromotive force proportional to the rotational speed ω. The following equation (17) on the γ-axis and δ-axis in the above equation (15), which is the state equation of the current circuit of the PMSM, is calculated.

Figure 0005534991
Figure 0005534991

電流目標算出部18は、上の(17)式の演算によって、d軸、q軸の電流推定値id^、iq^を出力する。   The current target calculation unit 18 outputs the d-axis and q-axis current estimation values id ^ and iq ^ by the calculation of the above equation (17).

速度位置推定部19は、電流目標算出部18から出力されるd軸、q軸の電流推定値id^、iq^、3相/2相座標変換部17から出力される実際の電流値iδ、iγに基づいて、回転速度ω^および磁極位置θ^を推定する。   The speed position estimator 19 includes the d-axis and q-axis current estimated values id ^ and iq ^ output from the current target calculator 18, the actual current values iδ output from the three-phase / two-phase coordinate converter 17, Based on iγ, the rotational speed ω ^ and the magnetic pole position θ ^ are estimated.

速度位置推定部19を実現するための1つの例は、上に述べた(8)式で示される。つまり、電流目標算出部18からのd軸、q軸の電流推定値id^、iq^と3相/2相座標変換部17からの電流値iδ、iγとの偏差から、回転速度ω^を推定する。磁極位置θ^は、推定された回転速度ω^を積分することによって求め、これを磁極位置の推定値とする。   One example for realizing the speed position estimation unit 19 is expressed by the above-described equation (8). That is, the rotational speed ω ^ is determined from the deviation between the d-axis and q-axis current estimation values id ^ and iq ^ from the current target calculation unit 18 and the current values iδ and iγ from the three-phase / two-phase coordinate conversion unit 17. presume. The magnetic pole position θ ^ is obtained by integrating the estimated rotational speed ω ^, and this is used as the estimated value of the magnetic pole position.

また、速度位置推定部19を実現するための他の1つの例は、上に述べた(9)式で示される。つまり、電流目標算出部18からのd軸、q軸の電流推定値id^、iq^と3相/2相座標変換部17からの電流値iδ、iγとの偏差から、回転速度ω^を推定する。   Another example for realizing the speed position estimating unit 19 is expressed by the above-described equation (9). That is, the rotational speed ω ^ is determined from the deviation between the d-axis and q-axis current estimation values id ^ and iq ^ from the current target calculation unit 18 and the current values iδ and iγ from the three-phase / two-phase coordinate conversion unit 17. presume.

この場合に、推定された回転速度ω^の回転方向を示す「sgnω^」は、「ω^」に置き代わっている。   In this case, “sgnω ^” indicating the rotational direction of the estimated rotational speed ω ^ is replaced by “ω ^”.

また、速度位置推定部19を実現するためのさらに他の1つの例は、上に述べた(10)式で示される。つまり、電流目標算出部18からのd軸、q軸の電流推定値id^、iq^と3相/2相座標変換部17からの電流値iδ、iγとの偏差から、回転速度ω^を推定する。   Still another example for realizing the speed position estimating unit 19 is expressed by the above-described equation (10). That is, the rotational speed ω ^ is determined from the deviation between the d-axis and q-axis current estimation values id ^ and iq ^ from the current target calculation unit 18 and the current values iδ and iγ from the three-phase / two-phase coordinate conversion unit 17. presume.

この場合に、推定された回転速度ω^の回転方向を示す「sgnω^」は、d軸とγ軸との角度の偏差(θ−θ^)の正弦である「sin(θ−θ^)」に置き代わっている。また、この場合において、非干渉制御部13は、上に述べた(16)式に代えて(11)式を演算する。   In this case, “sgnω ^” indicating the rotation direction of the estimated rotational speed ω ^ is “sin (θ−θ ^)”, which is a sine of an angular deviation (θ−θ ^) between the d axis and the γ axis. Has been replaced. In this case, the non-interference control unit 13 calculates the expression (11) instead of the expression (16) described above.

本実施形態の制御装置1は、1つまたは複数のCPUまたはDSP、ROM、RAM、周辺回路素子、インタフェ−ス回路素子、論理回路素子、スイッチング回路素子、IC回路素子、その他の半導体素子、受動回路素子、種々の電子機器、電子部品などを用いて構成することができる。上に説明した演算は、メモリに記憶されたプログラムをCPUまたはDSPが実行するソフトウエアにより、または論理回路などによるハードウエアにより、またはこれらの組み合わせにより、実現することができる。   The control device 1 of the present embodiment includes one or more CPUs or DSPs, ROM, RAM, peripheral circuit elements, interface circuit elements, logic circuit elements, switching circuit elements, IC circuit elements, other semiconductor elements, passive It can be configured using circuit elements, various electronic devices, electronic components, and the like. The operations described above can be realized by software executed by a CPU or DSP, hardware such as a logic circuit, or a combination of these programs stored in a memory.

上に述べたように、本実施形態の制御装置1によると、PMSMの回転速度をセンサレスで制御するに際し、従来のようなオブザーバを用いることなく、演算量の増大を抑制することができる。   As described above, according to the control device 1 of the present embodiment, when controlling the rotational speed of the PMSM without using a sensor, an increase in the amount of calculation can be suppressed without using a conventional observer.

制御装置1による適応制御について、さらに説明する。
〔モデル規範適応システムに基づくPMSMのセンサレス制御法〕
PMSMのセンサレス制御では、PMSMモデルのパラメータの正確さが、磁極位置推定の精度や性能に強く影響する。パラメータの変動に追従して同定しながら積極的にシステムのロバスト化を図る制御技術として適応システムがある。負荷の慣性モーメント同定にモデル規範適応システム(MRAS)が有効であることは既に確認されている(小原正樹著、「速度制御系のオートチューニング」平成元年電気学会全国大会論文集1689を参照できる)。ここに、抵抗、リアクトル、界磁定数などのパラメータは既知として、回転速度を未知パラメータとしたMRASを構築し、PMSMの磁極位置センサレス制御を実現する手法を提案する。以下の説明では、その適応同定アルゴリズムの導出を示す。
The adaptive control by the control device 1 will be further described.
[PMSM sensorless control method based on model reference adaptive system]
In PMSM sensorless control, the accuracy of PMSM model parameters strongly affects the accuracy and performance of magnetic pole position estimation. There is an adaptive system as a control technique for actively making the system robust while following and identifying parameter variations. It has already been confirmed that the model reference adaptive system (MRAS) is effective for identifying the moment of inertia of the load. ). Here, we propose a method for constructing MRAS with unknown parameters such as resistance, reactor, and field constant, and realizing PMSM magnetic pole position sensorless control with rotational speed as an unknown parameter. In the following description, the derivation of the adaptive identification algorithm is shown.

磁極位置の適応同定アルゴリズムについて説明する。   An adaptive identification algorithm for the magnetic pole position will be described.

まず、PMSMのd−q回転座標上における電圧電流方程式を(18)式に示す。   First, the voltage-current equation on the dsq rotation coordinate of PMSM is shown in equation (18).

Figure 0005534991
Figure 0005534991

図1にここで検討するモデル規範適応同定システムの構成が示されている。ここでは、実速度を一定と見なして速度を推定しているが、推定系の適応同定則に内部モデル原理を適用すれば、速度が変動する過度状態においても速度推定の誤差を少なくできると思われる。   FIG. 1 shows the configuration of the model reference adaptive identification system studied here. Here, the speed is estimated by assuming that the actual speed is constant, but if the internal model principle is applied to the adaptive identification rule of the estimation system, the error of speed estimation can be reduced even in an excessive state where the speed fluctuates. It is.

PMSMのセンサレス制御において、中・高速域用として、適応オブザ−バに基づく方式、拡張誘起電圧を用いた方式などがある。しかし、これらほとんどの方式は電動機の数学モデルを利用しており、センサレス制御の性能は使用するパラメータの精度に影響される。このため、停止や運転中での各種パラメータの計測技術の研究が盛んである。   In sensorless control of PMSM, there are a method based on an adaptive observer, a method using an extended induced voltage, and the like for medium and high speed regions. However, most of these methods use a mathematical model of the motor, and the performance of sensorless control is affected by the accuracy of the parameters used. For this reason, research on measurement techniques of various parameters during stoppage and operation is active.

この点に着目して、パラメータ同定にも拡張できるモデル規範適応システム(MRAS)に基づく磁極位置センサレス制御法を提案する。はじめに電流制御ループ内でのMRASに基づく非干渉センサレス制御の構成を述べる。次に磁極位置同定のための規範モデルの導出、安定性の証明、同定則の導出を述べる。   Focusing on this point, a magnetic pole position sensorless control method based on a model reference adaptive system (MRAS) that can be extended to parameter identification is proposed. First, the configuration of non-interference sensorless control based on MRAS in the current control loop will be described. Next, the derivation of the reference model for magnetic pole position identification, the proof of stability, and the derivation of the identification rule are described.

次に、MRASの構成について説明する。   Next, the configuration of MRAS will be described.

PMSMの位置・速度制御ではベクトル制御が一般的であり、しかもd−q回転座標間での速度起電力の影響を除くため非干渉制御が用いられている。ここで提案する磁極位置センサレス制御システムの構成が図1に示されている。図から明らかなように、これは従来のPMSM非干渉ベクトル制御システムにMRASの要素である適応モデルと速度・位置推定ブロックを付加した構造となっており、PMSMを駆動する推定回転座標(γ−δ軸)上で、しかも非干渉制御を付加した電流制御ループ内で、MRASを構成するのが特徴となっている。適応モデルは非干渉化されたPMSMの理想電流回路モデルであり、速度・位置推定は作成されたモデル電流id^、iq^と検出された実電流iγ、iδから磁極位置・回転速度を推定する。   In PMSM position / velocity control, vector control is generally used, and non-interference control is used in order to eliminate the influence of velocity electromotive force between dq rotation coordinates. The configuration of the magnetic pole position sensorless control system proposed here is shown in FIG. As is apparent from the figure, this is a structure in which an adaptive model and a velocity / position estimation block, which are elements of MRAS, are added to the conventional PMSM non-interference vector control system, and the estimated rotational coordinates (γ− It is characterized in that the MRAS is configured on the δ axis) and in a current control loop to which non-interference control is added. The adaptive model is a non-interfering PMSM ideal current circuit model, and the speed / position estimation is to estimate the magnetic pole position / rotation speed from the generated model currents id ^, iq ^ and the detected actual currents iγ, iδ. .

次に、磁極位置の適応同定について説明する。
(γ−δ回転軸での状態方程式)
PMSMのd−q回転座標上における電圧電流方程式が上の(18)式に示されている。
Next, adaptive identification of the magnetic pole position will be described.
(Equation of state at γ-δ rotation axis)
The voltage-current equation on the dsq rotation coordinate of PMSM is shown in the above equation (18).

上の(18)式を変形して電流を状態変数および出力変数、電圧を入力変数とするPMSMのd−q座標における状態方程式を求めると、次の(19)式となる。   When the above equation (18) is modified to obtain a state equation in the dq coordinate of PMSM using current as a state variable and output variable and voltage as an input variable, the following equation (19) is obtained.

Figure 0005534991
Figure 0005534991

一方、γ−δ推定回転座標上における状態方程式を求めると次の(20)式となる。   On the other hand, when the state equation on the γ-δ estimated rotation coordinate is obtained, the following equation (20) is obtained.

Figure 0005534991
Figure 0005534991

(非干渉制御と適応モデル)
上の(19)式ではd軸q軸それぞれお互いに干渉する速度起電力項が存在するので、
(Non-interference control and adaptive model)
In the above equation (19), there are velocity electromotive force terms that interfere with each other on the d-axis and the q-axis.

Figure 0005534991
Figure 0005534991

上の(21)式によって非干渉制御を施して干渉項を取り除き、求めた式を規範モデルとすると次の(22)式となる。   When the interference term is removed by performing non-interference control by the above equation (21) and the obtained equation is a reference model, the following equation (22) is obtained.

Figure 0005534991
Figure 0005534991

推定回転座標(γ−δ軸)の(20a)式に非干渉制御を施した後、確定的等価原理(CE原理)よりPMSMに加える制御則(電圧入力)を求める。   After performing non-interference control on the equation (20a) of the estimated rotational coordinate (γ-δ axis), a control law (voltage input) to be applied to the PMSM is obtained from a deterministic equivalent principle (CE principle).

Figure 0005534991
Figure 0005534991

上の(23)式を(20a)式に代入すると、推定回転座標(γ−δ軸)上での誘起電圧と推定誘起電圧を含んだ次の(24)式を求めることができる。   When the above equation (23) is substituted into the equation (20a), the following equation (24) including the induced voltage on the estimated rotational coordinate (γ-δ axis) and the estimated induced voltage can be obtained.

Figure 0005534991
Figure 0005534991

(誤差方程式)
規範モデルと実際値との誤差εγ,εδを以下の(25)式で定義する。
(Error equation)
Errors εγ and εδ between the reference model and actual values are defined by the following equation (25).

Figure 0005534991
Figure 0005534991

上の(22)式および(24)式から、抵抗、リアクトル、誘起電圧定数などのパラメータは既知とし、回転速度のみを未知パラメータとする誤差方程式を求める。(22)式から(24)式をひくと、   From the above equations (22) and (24), an error equation is obtained in which parameters such as resistance, reactor, and induced voltage constant are known and only the rotational speed is an unknown parameter. From equation (22) to equation (24),

Figure 0005534991
Figure 0005534991

となる。上の(26)式で推定軸上の誘起電圧   It becomes. Induced voltage on the estimated axis in equation (26) above

Figure 0005534991
Figure 0005534991

をそれぞれ次の(28)式のように近似できるので、 Can be approximated as in the following equation (28).

Figure 0005534991
Figure 0005534991

上の(28a)式および(28b)式と次の(29)式、 The above equations (28a) and (28b) and the following equation (29):

Figure 0005534991
Figure 0005534991

とを上の(26)式に代入すると、 Is substituted into the above equation (26),

Figure 0005534991
Figure 0005534991

上の(30)式において、磁極位置と回転速度の偏差が零となるように、つまりθ−θ^→0、ω−ω^→0となるように、MRASを構成する。磁極位置がθ−θ^→0ならば回転速度はω−ω^→0となるため、(30)式の回転速度偏差ω−ω^を|ω|(θ−θ^)に置き換えて安定性を論じても問題を生じない。したがって、誤差方程式は(30)式を変形した次の(31)式となる。   In the above equation (30), the MRAS is configured so that the deviation between the magnetic pole position and the rotation speed becomes zero, that is, θ−θ ^ → 0 and ω−ω ^ → 0. If the magnetic pole position is θ−θ ^ → 0, the rotation speed will be ω−ω ^ → 0, so the rotation speed deviation ω−ω ^ in equation (30) is replaced with | ω | (θ−θ ^) Discussing sex does not cause problems. Therefore, the error equation is the following equation (31) obtained by modifying equation (30).

Figure 0005534991
Figure 0005534991

(安定性と適応則)
MRASは非線形であり、ポポフの超安定論に基づいて安定性を証明する。なお、ここでは、速度位置推定部19が上の(8)式を用いた場合についての証明を行うが、速度位置推定部19が上の(9)式または(10)式を用いた場合についても、式を変形することによって同様に証明することが可能である。
(Stability and adaptation law)
MRAS is non-linear and proves stability based on Popov's theory of superstability. Here, although the proof that the speed position estimation unit 19 uses the above equation (8) is proved, the speed position estimation unit 19 uses the above equation (9) or (10). Can be proved in the same way by transforming the equation.

図5に非線形フィードバック系が示されている。上の(31)式において、   FIG. 5 shows a non-linear feedback system. In the above equation (31),

Figure 0005534991
Figure 0005534991

とおくと、 After all,

Figure 0005534991
Figure 0005534991

となる。(33)式は、入力u、出力eの関数に変形でき、図5の線形定常ブロック伝達関数G(s)に相当する。(33)式において行列Aが漸近安定行列であるため、G(s)は強正実となり、ポポフの第一条件を満たす。さらに図5のフィードバック系が安定であるためには、第二の条件であるポポフの積分不等式(45)を満足する必要がある。 It becomes. Expression (33) can be transformed into a function of input u and output e, and corresponds to the linear stationary block transfer function G (s) in FIG. In Equation (33), since the matrix A is an asymptotically stable matrix, G (s) is strongly positive and satisfies Popov's first condition. Furthermore, in order for the feedback system of FIG. 5 to be stable, it is necessary to satisfy Popov's integral inequality (45), which is the second condition.

Figure 0005534991
Figure 0005534991

上の(32)式においてw=−uであるので、   Since w = −u in the above equation (32),

Figure 0005534991
Figure 0005534991

次の(35)式を(34)式に代入すると、   Substituting the following equation (35) into equation (34),

Figure 0005534991
Figure 0005534991

となる。回転速度ωを一定と仮定すると磁極位置θはt の一次関数(1/s2 )となるので、内部モデル原理よりθの推定値に定常偏差を生じさせないためには適応則に1/s2 の成分を含む必要がある。ω−ω^の適応則を、 It becomes. Assuming that the rotational speed ω is constant, the magnetic pole position θ becomes a linear function (1 / s 2 ) of t. Therefore, in order to prevent a steady deviation from occurring in the estimated value of θ from the internal model principle, the adaptive law is 1 / s 2 It is necessary to include the ingredients. The adaptation law of ω-ω ^

Figure 0005534991
Figure 0005534991

Figure 0005534991
Figure 0005534991

とおき、上の(37)式の比例+積分を積分した、 And the proportional + integral of the above equation (37) is integrated,

Figure 0005534991
Figure 0005534991

を磁極位置誤差の適応則とする。上の(39)式より、 Is an adaptive law of magnetic pole position error. From equation (39) above,

Figure 0005534991
Figure 0005534991

となり、上の(36)式、(38)式、および(40)式より、 From the above equations (36), (38), and (40),

Figure 0005534991
Figure 0005534991

が求まる。このため、簡単な計算で次の(42)式が証明でき、 Is obtained. Therefore, the following equation (42) can be proved by a simple calculation,

Figure 0005534991
Figure 0005534991

ポポフの積分不等式(34)を満足するため、MRASは漸近安定となり、次の(43)式のように誤差は0に収束する。 Since Ms. Popov's integral inequality (34) is satisfied, MRAS becomes asymptotically stable, and the error converges to 0 as in the following equation (43).

Figure 0005534991
Figure 0005534991

以上より、推定磁極位置θ^と推定回転速度ω^の実際の適応則はωを一定と仮定し、ωとω^の回転方向は同一であるのでω^の極性(回転方向)を使用することにより、   From the above, the actual adaptive law of the estimated magnetic pole position θ ^ and the estimated rotation speed ω ^ assumes that ω is constant, and the rotation direction of ω and ω ^ is the same, so the polarity of ω ^ (rotation direction) is used. By

Figure 0005534991
Figure 0005534991

Figure 0005534991
Figure 0005534991

となる。
〔むすび〕
上に述べたように、本実施形態の制御装置1では、理想モデル(規範モデル)と実測値との差を求め、その偏差εγ(=id^−iγ)、εδ(=iq^−iδ)、εδ−εγsgnω^が安定的に零となるような推定回転速度ω^を求める。このような制御によると、少ない演算量で高速な制御を行うことができる。
It becomes.
[Musubi]
As described above, in the control device 1 of the present embodiment, the difference between the ideal model (normative model) and the actually measured value is obtained, and the deviations εγ (= id ^ -iγ), εδ (= iq ^ -iδ) , Εδ−εγsgnω ^ is obtained so that the estimated rotational speed ω ^ is stably zero. According to such control, high-speed control can be performed with a small amount of calculation.

上に述べた制御装置1において、速度調節部11は、偏差(ω*−ω^)が零になるように、比例・積分制御(PI制御)を行ったが、比例・積分・微分制御(PID制御)、その他の制御を行ってもよい。速度調節部11において、トルク指令値を求め、そのトルクを得るための電流指令値iγ*、iδ*を求めてもよい。速度調節部11から出力される電流指令値iγ*、iδ*を用いて、3相/2相座標変換部17から出力される電流検出値iγ、iδとの偏差(iγ*−iγ)、(iδ*−iδ)を演算し、これを電流調節部12a,12bに入力すればよい。   In the control device 1 described above, the speed adjustment unit 11 performs the proportional / integral control (PI control) so that the deviation (ω * −ω ^) becomes zero. PID control) and other controls may be performed. The speed adjusting unit 11 may obtain a torque command value and obtain current command values iγ * and iδ * for obtaining the torque. Deviations (iγ * −iγ) from the detected current values iγ and iδ output from the three-phase / two-phase coordinate conversion unit 17 using the current command values iγ * and iδ * output from the speed adjusting unit 11 ( i δ * −i δ) may be calculated and input to the current adjusting units 12a and 12b.

その他、PMSM、電流調節部12、非干渉制御部13、電流目標算出部18、速度位置推定部19、または制御装置1の各部または全体の構成、構造、回路、形状、方式、個数、演算内容、処理内容、処理順序などは、本発明の主旨に沿って適宜変更することができる。   In addition, PMSM, current adjustment unit 12, non-interference control unit 13, current target calculation unit 18, speed position estimation unit 19, or each part or entire configuration of control device 1, structure, circuit, shape, method, number, calculation contents The processing content, processing order, and the like can be changed as appropriate in accordance with the gist of the present invention.

1 制御装置
11 速度調節部
12 電流調節部
13 非干渉制御部
14 2相/3相座標変換部
15 電力変換部
16 電流検出部
17 3相/2相座標変換部
18 電流目標算出部
19 速度位置推定部
ω* 回転速度指令値
ω^ 回転速度推定値
ω 回転速度実際値
θ^ 推定磁極位置
θ 磁極位置実際値
iγ* γ軸電流指令値
iδ* δ軸電流指令値
iγ γ軸電流実際値
iδ δ軸電流実際値
id^ d軸推定電流値
iq^ q軸推定電流値
Vγ* γ軸電流調節部の指令値
Vδ* δ軸電流調節部の指令値
Vγ** γ軸電圧指令値
Vδ** δ軸電圧指令値
DESCRIPTION OF SYMBOLS 1 Control apparatus 11 Speed adjustment part 12 Current adjustment part 13 Non-interference control part 14 2-phase / 3-phase coordinate conversion part 15 Power conversion part 16 Current detection part 17 3-phase / 2-phase coordinate conversion part 18 Current target calculation part 19 Speed position Estimator ω * rotational speed command value ω ^ rotational speed estimated value ω rotational speed actual value θ ^ estimated magnetic pole position θ magnetic pole position actual value iγ * γ-axis current command value iδ * δ-axis current command value iγ γ-axis current actual value iδ δ-axis current actual value id ^ d-axis estimated current value iq ^ q-axis estimated current value Vγ * γ-axis current adjuster command value Vδ * δ-axis current adjuster command value Vγ ** γ-axis voltage command value Vδ ** δ-axis voltage command value

Claims (5)

同期電動機の回転速度を制御するための制御装置であって、
速度指令に基づいて速度調節を行う速度調節部と、
電流調節を行う電流調節部と、
前記電流調節部から出力される電圧指令と回転速度の推定値とに基づいて干渉分を除去するための非干渉制御を行う非干渉制御部と、
前記非干渉制御部の出力する電圧を座標変換する2相/3相座標変換部と、
前記2相/3相座標変換部の出力に基づいて電機子捲線に電流を流して駆動するための電力変換部と、
前記電機子捲線に流れる電流を検出する電流検出部と、
前記電流検出部により検出された電流を座標変換する3相/2相座標変換部と、
前記電圧指令に基づいて干渉分を除去した制御モデルにおける電流の目標値を算出する電流目標算出部と、
前記電流検出部により検出された電流と前記電流目標算出部により算出された電流の目標値との差が漸近的に零になるように回転速度の推定値および磁極位置の推定値を求める速度位置推定部と、を有し、
前記制御モデルは、
その電圧電流方程式が、次の(1)式で示され、
Figure 0005534991

前記速度位置推定部は、
次の(2)式に基づいて、前記回転速度の推定値および磁極位置の推定値を求める、
Figure 0005534991

ことを特徴とする同期電動機の制御装置。
A control device for controlling the rotational speed of the synchronous motor,
A speed adjusting unit for adjusting the speed based on the speed command;
A current adjustment unit for current adjustment;
A non-interference control unit that performs non-interference control to remove interference based on the voltage command output from the current adjustment unit and the estimated value of the rotation speed;
A two-phase / three-phase coordinate conversion unit for converting the voltage output from the non-interference control unit;
A power converter for driving the armature winding based on the output of the two-phase / three-phase coordinate converter;
A current detector for detecting a current flowing in the armature winding;
A three-phase / two-phase coordinate conversion unit that converts the current detected by the current detection unit;
A current target calculation unit for calculating a target value of a current in a control model from which interference is removed based on the voltage command;
Speed position for obtaining the estimated value of the rotational speed and the estimated value of the magnetic pole position so that the difference between the current detected by the current detector and the target value of the current calculated by the current target calculator is asymptotically zero. an estimation unit, the possess,
The control model is
The voltage-current equation is expressed by the following equation (1):
Figure 0005534991

The speed position estimation unit
Based on the following equation (2), an estimated value of the rotational speed and an estimated value of the magnetic pole position are obtained.
Figure 0005534991

The control apparatus of the synchronous motor characterized by the above-mentioned.
同期電動機の回転速度を制御するための制御装置であって、
速度指令に基づいて速度調節を行う速度調節部と、
電流調節を行う電流調節部と、
前記電流調節部から出力される電圧指令と回転速度の推定値とに基づいて干渉分を除去するための非干渉制御を行う非干渉制御部と、
前記非干渉制御部の出力する電圧を座標変換する2相/3相座標変換部と、
前記2相/3相座標変換部の出力に基づいて電機子捲線に電流を流して駆動するための電力変換部と、
前記電機子捲線に流れる電流を検出する電流検出部と、
前記電流検出部により検出された電流を座標変換する3相/2相座標変換部と、
前記電圧指令に基づいて干渉分を除去した制御モデルにおける電流の目標値を算出する電流目標算出部と、
前記電流検出部により検出された電流と前記電流目標算出部により算出された電流の目標値との差が漸近的に零になるように回転速度の推定値および磁極位置の推定値を求める速度位置推定部と、を有し、
前記制御モデルは、
その電圧電流方程式が、次の(1)式で示され、
Figure 0005534991

前記速度位置推定部は、
次の(3)式に基づいて、前記回転速度の推定値および磁極位置の推定値を求める、
Figure 0005534991

ことを特徴とする同期電動機の制御装置。
A control device for controlling the rotational speed of the synchronous motor,
A speed adjusting unit for adjusting the speed based on the speed command;
A current adjustment unit for current adjustment;
A non-interference control unit that performs non-interference control to remove interference based on the voltage command output from the current adjustment unit and the estimated value of the rotation speed;
A two-phase / three-phase coordinate conversion unit for converting the voltage output from the non-interference control unit;
A power converter for driving the armature winding based on the output of the two-phase / three-phase coordinate converter;
A current detector for detecting a current flowing in the armature winding;
A three-phase / two-phase coordinate conversion unit that converts the current detected by the current detection unit;
A current target calculation unit for calculating a target value of a current in a control model from which interference is removed based on the voltage command;
Speed position for obtaining the estimated value of the rotational speed and the estimated value of the magnetic pole position so that the difference between the current detected by the current detector and the target value of the current calculated by the current target calculator is asymptotically zero. an estimation unit, the possess,
The control model is
The voltage-current equation is expressed by the following equation (1):
Figure 0005534991

The speed position estimation unit
Based on the following equation (3), an estimated value of the rotational speed and an estimated value of the magnetic pole position are obtained.
Figure 0005534991

The control apparatus of the synchronous motor characterized by the above-mentioned.
同期電動機の回転速度を制御するための制御装置であって、
速度指令に基づいて速度調節を行う速度調節部と、
電流調節を行う電流調節部と、
前記電流調節部から出力される電圧指令と回転速度の推定値とに基づいて干渉分を除去するための非干渉制御を行う非干渉制御部と、
前記非干渉制御部の出力する電圧を座標変換する2相/3相座標変換部と、
前記2相/3相座標変換部の出力に基づいて電機子捲線に電流を流して駆動するための電力変換部と、
前記電機子捲線に流れる電流を検出する電流検出部と、
前記電流検出部により検出された電流を座標変換する3相/2相座標変換部と、
前記電圧指令に基づいて干渉分を除去した制御モデルにおける電流の目標値を算出する電流目標算出部と、
前記電流検出部により検出された電流と前記電流目標算出部により算出された電流の目標値との差が漸近的に零になるように回転速度の推定値および磁極位置の推定値を求める速度位置推定部と、を有し、
前記制御モデルは、
その電圧電流方程式が、次の(1)式で示され、
Figure 0005534991

前記速度位置推定部は、
次の(4)式に基づいて、前記回転速度の推定値および磁極位置の推定値を求める、
Figure 0005534991

ことを特徴とする同期電動機の制御装置。
A control device for controlling the rotational speed of the synchronous motor,
A speed adjusting unit for adjusting the speed based on the speed command;
A current adjustment unit for current adjustment;
A non-interference control unit that performs non-interference control to remove interference based on the voltage command output from the current adjustment unit and the estimated value of the rotation speed;
A two-phase / three-phase coordinate conversion unit for converting the voltage output from the non-interference control unit;
A power converter for driving the armature winding based on the output of the two-phase / three-phase coordinate converter;
A current detector for detecting a current flowing in the armature winding;
A three-phase / two-phase coordinate conversion unit that converts the current detected by the current detection unit;
A current target calculation unit for calculating a target value of a current in a control model from which interference is removed based on the voltage command;
Speed position for obtaining the estimated value of the rotational speed and the estimated value of the magnetic pole position so that the difference between the current detected by the current detector and the target value of the current calculated by the current target calculator is asymptotically zero. an estimation unit, the possess,
The control model is
The voltage-current equation is expressed by the following equation (1):
Figure 0005534991

The speed position estimation unit
Based on the following equation (4), an estimated value of the rotational speed and an estimated value of the magnetic pole position are obtained.
Figure 0005534991

The control apparatus of the synchronous motor characterized by the above-mentioned.
前記非干渉制御部は、
次の(5)式に基づいて非干渉制御を行う、
Figure 0005534991

請求項または記載の同期電動機の制御装置。
The non-interference control unit is
Perform non-interference control based on the following equation (5).
Figure 0005534991

The control apparatus for a synchronous motor according to claim 1 or 2 .
前記非干渉制御部は、
次の(6)式に基づいて非干渉制御を行う、
Figure 0005534991

請求項記載の同期電動機の制御装置。
The non-interference control unit is
Non-interference control is performed based on the following equation (6).
Figure 0005534991

The control apparatus for a synchronous motor according to claim 3 .
JP2010162169A 2010-07-16 2010-07-16 Control device for synchronous motor Active JP5534991B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010162169A JP5534991B2 (en) 2010-07-16 2010-07-16 Control device for synchronous motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010162169A JP5534991B2 (en) 2010-07-16 2010-07-16 Control device for synchronous motor

Publications (2)

Publication Number Publication Date
JP2012023937A JP2012023937A (en) 2012-02-02
JP5534991B2 true JP5534991B2 (en) 2014-07-02

Family

ID=45777688

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010162169A Active JP5534991B2 (en) 2010-07-16 2010-07-16 Control device for synchronous motor

Country Status (1)

Country Link
JP (1) JP5534991B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5798513B2 (en) * 2012-03-09 2015-10-21 株式会社竹中製作所 Method and apparatus for detecting initial magnetic pole position of permanent magnet synchronous motor, and control apparatus for permanent magnet synchronous motor
JP6397734B2 (en) * 2014-11-14 2018-09-26 株式会社竹中製作所 DC motor control device
CN105071719B (en) * 2015-07-27 2017-12-22 哈尔滨工业大学 Position-sensor-free evaluation method for low capacity bus capacitor driver

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3736287B2 (en) * 2000-05-22 2006-01-18 日産自動車株式会社 Vector control device for electric motor
JP5165545B2 (en) * 2008-12-08 2013-03-21 本田技研工業株式会社 Electric motor magnetic pole position estimation device

Also Published As

Publication number Publication date
JP2012023937A (en) 2012-02-02

Similar Documents

Publication Publication Date Title
JP5156352B2 (en) AC motor control device
JP6220941B2 (en) Induction motor parameter estimation method
JP5223109B2 (en) Control device for permanent magnet type synchronous motor
JP6776066B2 (en) Inverter controller and motor drive system
KR101485989B1 (en) Motor control device
JP2010035363A (en) Controller for permanent magnet type synchronous motor
JP6166601B2 (en) Motor control device and generator control device
JP6135713B2 (en) Motor control device, magnetic flux command generation device, and magnetic flux command generation method
JP5534991B2 (en) Control device for synchronous motor
JP2013150498A (en) Controller and control method of synchronous motor
JP5641774B2 (en) Method and apparatus for estimating rotor position and speed of stepping motor
JP2013146155A (en) Winding temperature estimating device and winding temperature estimating method
JP2011067066A (en) Controller for permanent magnet type synchronous motor
JP5332305B2 (en) Control device for permanent magnet type synchronous motor
JP6590196B2 (en) Power converter
JP5332301B2 (en) Control device for permanent magnet type synchronous motor
JP2023048833A (en) State estimation method for motor unit and state estimation device
JP7251424B2 (en) INVERTER DEVICE AND INVERTER DEVICE CONTROL METHOD
JP6108114B2 (en) Control device for permanent magnet type synchronous motor
JP5798513B2 (en) Method and apparatus for detecting initial magnetic pole position of permanent magnet synchronous motor, and control apparatus for permanent magnet synchronous motor
JP6848680B2 (en) Synchronous motor control device
JP5186352B2 (en) Electric motor magnetic pole position estimation device
JP2010028981A (en) Rotor position estimating method for synchronous motor, and controller for the synchronous motor
JP6497584B2 (en) Control device for permanent magnet type synchronous motor
JP2015192463A (en) Motor control device and power generator control device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130502

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140225

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140228

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140317

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140422

R150 Certificate of patent or registration of utility model

Ref document number: 5534991

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250