[go: up one dir, main page]

JP6726390B2 - Controller for permanent magnet type synchronous motor - Google Patents

Controller for permanent magnet type synchronous motor Download PDF

Info

Publication number
JP6726390B2
JP6726390B2 JP2016133789A JP2016133789A JP6726390B2 JP 6726390 B2 JP6726390 B2 JP 6726390B2 JP 2016133789 A JP2016133789 A JP 2016133789A JP 2016133789 A JP2016133789 A JP 2016133789A JP 6726390 B2 JP6726390 B2 JP 6726390B2
Authority
JP
Japan
Prior art keywords
current
command value
value
torque
electric motor
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
JP2016133789A
Other languages
Japanese (ja)
Other versions
JP2018007473A (en
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.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric 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 Fuji Electric Co Ltd filed Critical Fuji Electric Co Ltd
Priority to JP2016133789A priority Critical patent/JP6726390B2/en
Publication of JP2018007473A publication Critical patent/JP2018007473A/en
Application granted granted Critical
Publication of JP6726390B2 publication Critical patent/JP6726390B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Motors That Do Not Use Commutators (AREA)

Description

本発明は、磁極位置検出器を持たない永久磁石形同期電動機の制御装置に関し、詳しくは、電動機の速度範囲に応じて各種の制御方式を切替える際のショックを低減して電動機をスムーズに加減速運転するための制御装置に関するものである。 The present invention relates to a controller for a permanent magnet type synchronous motor that does not have a magnetic pole position detector, and more specifically, it reduces shock when switching between various control methods according to the speed range of the motor to smoothly accelerate and decelerate the motor. The present invention relates to a control device for driving.

永久磁石形同期電動機の制御装置をコストダウンするため、回転子の磁極位置を検出するための磁極位置検出器を使用しないで運転する、いわゆる、センサレスベクトル制御が実用化されている。
センサレスベクトル制御は、電動機の端子電圧や電流の情報から回転子の磁極位置と速度とを演算し、これらに基づいて電流制御を行うことでトルク制御や速度制御を実現するものである。なお、センサレスベクトル制御の具体的な方法は良く知られているため、ここでは詳述を省略する。
In order to reduce the cost of the control device for the permanent magnet type synchronous motor, so-called sensorless vector control has been put into practical use, which operates without using a magnetic pole position detector for detecting the magnetic pole position of the rotor.
The sensorless vector control realizes torque control and speed control by calculating the magnetic pole position and speed of the rotor from information on the terminal voltage and current of the electric motor and performing current control based on these. Since a specific method of sensorless vector control is well known, its detailed description is omitted here.

しかしながら、センサレスベクトル制御は、始動時を始めとする低速域における安定性に問題があることから、電動機の電流の振幅を一定とし、電流の角速度を速度指令値に制御することで、電動機の回転子を電流に引き込んで運転する技術が適用されることがある。このような運転方式を、以下では「電流引込制御」と呼ぶ。 However, since the sensorless vector control has a problem in stability in the low speed range such as at the time of starting, the rotation speed of the electric motor is controlled by keeping the amplitude of the electric current of the electric motor constant and controlling the angular velocity of the electric current to the speed command value. A technique for drawing a child into an electric current to drive the child may be applied. Hereinafter, such an operation method is referred to as “current pull-in control”.

ここで、例えば特許文献1には、電流引込制御とセンサレスベクトル制御との相互間で制御方式を切り替える際の電圧や電流の急変を防止してショックを低減するようにした永久磁石形同期電動機の制御装置が開示されている。特に、特許文献1の請求項6や明細書の段落[0083]〜[0093]、図6等には、運転中の電動機を減速して停止する場合に、センサレスベクトル制御から電流引込制御へショックレスにて切り替える技術が記載されている。 Here, for example, in Patent Document 1, there is disclosed a permanent magnet type synchronous motor which is configured to prevent a sudden change in voltage or current when switching control methods between current pull-in control and sensorless vector control to reduce shock. A controller is disclosed. Particularly, in claim 6 of Patent Document 1, paragraphs [0083] to [0093] of the specification, FIG. 6 and the like, when decelerating and stopping an electric motor during operation, shock is transferred from sensorless vector control to current drawing control. It describes the technology to switch between less.

特許第5277724号公報(請求項6、段落[0083]〜[0093]、図6等)Japanese Patent No. 5277724 (claim 6, paragraphs [0083] to [0093], FIG. 6, etc.)

上述した特許文献1では、センサレスベクトル制御から電流引込制御への切替前後で電動機の発生トルクが変化しないように、切替時におけるq軸電流指令値の初期値を数式23により演算すると共に、d軸電流指令値の初期値を、数式23により求めたq軸電流指令値の初期値と電流指令値振幅とを用いて数式24により演算している。
しかしながら、数式23はあくまで近似式であるため、その近似誤差に起因して電流指令値が適切な値にならず、結果として、電流引込制御への切替時に電動機の発生トルクにショックが発生する場合があった。
In Patent Document 1 described above, the initial value of the q-axis current command value at the time of switching is calculated by Equation 23 and the d-axis is set so that the generated torque of the electric motor does not change before and after switching from sensorless vector control to current drawing control. The initial value of the current command value is calculated by Expression 24 using the initial value of the q-axis current command value obtained by Expression 23 and the current command value amplitude.
However, since the formula 23 is only an approximate formula, the current command value does not become an appropriate value due to the approximation error, and as a result, a shock occurs in the torque generated by the electric motor when switching to the current drawing control. was there.

そこで、本発明の解決課題は、センサレスベクトル制御から電流引込制御への切替時における電流指令値の初期値を適切に演算することにより、ショックレスにて制御方式の切替を可能にした永久磁石形同期電動機の制御装置を提供することにある。 Therefore, a problem to be solved by the present invention is to permanently calculate the initial value of the current command value at the time of switching from the sensorless vector control to the current pull-in control, so that the control method can be switched without a shock. It is to provide a control device for a synchronous motor.

上記課題を解決するため、請求項1に係る発明は、磁極位置検出器を持たない永久磁石形同期電動機を電力変換器により運転するための制御装置において、
前記電動機の電流及び端子電圧をベクトルとしてとらえ、
前記電動機の電流指令値の角速度を速度指令値に制御する第1の運転モードと、
前記電動機の電流及び端子電圧から演算した回転子の磁極位置推定値及び速度推定値を用いて前記電動機の速度を速度指令値に制御する第2の運転モードと、を有し、
前記第2の運転モードから前記第1の運転モードへ切替えるときに、
前記電動機のトルク指令値と前記電流指令値の振幅とから、前記電流指令値と回転子の磁極位置との角度差を演算する手段と、
前記角度差から電気角補正値を求める手段と、
前記電気角補正値を用いて前記電流指令値の角度を初期化する手段と、
前記電気角補正値を用いて、前記電動機の端子電圧が急変しないように前記電動機の端子電圧指令値を初期化する手段と、
前記速度指令値と前記速度推定値との偏差から求めた角速度補償値を用いて、前記電流指令値の角速度が急変しないように前記電流指令値の角速度を補正する手段と、
を備え
前記電流指令値と前記磁極位置との角度差を演算する手段は、
前記電流指令値の振幅から、トルクが前記角度差の増加関数となる、前記角度差の区間を演算する手段と、前記角度差の区間から、前記トルク指令値に前記電動機のトルクが一致するように前記角度差を演算する手段と、を有するものである。
本発明によれば、第2の運転モードであるセンサレスベクトル制御から第1の運転モードである電流引込制御へ切替える際のショックを従来よりも低減することができ、永久磁石形同期電動機をスムーズに運転可能であると共に、上述した制御切替時の演算を簡略化することができる。
In order to solve the above-mentioned problems, the invention according to claim 1 is a control device for operating a permanent magnet type synchronous motor having no magnetic pole position detector by a power converter,
Capturing the current and terminal voltage of the motor as a vector,
A first operation mode for controlling an angular velocity of a current command value of the electric motor to a speed command value;
A second operation mode in which the speed of the electric motor is controlled to a speed command value using an estimated magnetic pole position value and an estimated speed value of the rotor calculated from the electric current and the terminal voltage of the electric motor,
When switching from the second operation mode to the first operation mode,
Means for calculating an angular difference between the current command value and the magnetic pole position of the rotor from the torque command value of the electric motor and the amplitude of the current command value;
Means for obtaining an electrical angle correction value from the angle difference,
Means for initializing the angle of the current command value using the electrical angle correction value;
Means for initializing a terminal voltage command value of the electric motor so that the terminal voltage of the electric motor does not suddenly change using the electrical angle correction value;
A means for correcting the angular velocity of the current command value so that the angular velocity of the current command value does not change suddenly using the angular velocity compensation value obtained from the deviation between the speed command value and the estimated speed value;
Equipped with
Means for calculating the angle difference between the current command value and the magnetic pole position,
A means for calculating a section of the angle difference in which torque is an increasing function of the angle difference from the amplitude of the current command value, and a torque of the electric motor matches the torque command value from the section of the angle difference. And means for calculating the angle difference .
According to the present invention, it is possible to reduce the shock at the time of switching from the sensorless vector control which is the second operation mode to the current drawing control which is the first operation mode, as compared with the prior art, and to make the permanent magnet type synchronous motor smooth. with luck a rolling Friendly ability can be simplified calculation at the time of control switching as described above.

本発明によれば、磁極位置検出器を持たない永久磁石形同期電動機の制御方式をセンサレスベクトル制御から電流引込制御へ切り替える際のショックを低減し、例えば減速時におけるスムーズな運転を実現することができる。 According to the present invention, it is possible to reduce the shock when switching the control method of the permanent magnet type synchronous motor having no magnetic pole position detector from the sensorless vector control to the current drawing control, and to realize smooth operation during deceleration, for example. it can.

本発明の実施形態の構成を示すブロック図である。It is a block diagram showing composition of an embodiment of the present invention. γ,δ軸及びd,q軸の定義を示す図である。It is a figure which shows the definition of (gamma), (delta) axis, and (d) and (q) axis. センサレスベクトル制御から電流引込制御への切替処理の原理を示すベクトル図である。It is a vector diagram which shows the principle of the switching process from sensorless vector control to current drawing control. 図1における角速度補償器の動作原理を示す図である。It is a figure which shows the operating principle of the angular velocity compensator in FIG. 電流引込制御におけるd軸と電流ベクトルとの角度差δと、トルクτとの関係を示す図である。It is a figure which shows the angle difference (delta) i of the d-axis and a current vector in current drawing control, and the relationship with torque (tau). 角度差δの上限値δimaxの演算原理を説明するためのベクトル図である。It is a vector diagram for explaining the calculation principle of the upper limit value δ imax of the angle difference δ i . 角度差δの下限値δiminの演算原理を説明するためのベクトル図である。It is a vector diagram for explaining the calculation principle of the lower limit value δ imin of the angle difference δ i .

以下、図に沿って本発明の実施形態を説明する。
図1は、この実施形態に係る制御装置を主回路と共に示したブロック図であり、永久磁石形同期電動機の制御方式を、拡張誘起電圧を利用したセンサレスベクトル制御から電流引込制御へショックレスにて切替えるためのものである。
Embodiments of the present invention will be described below with reference to the drawings.
FIG. 1 is a block diagram showing a control device according to this embodiment together with a main circuit. The control method of a permanent magnet type synchronous motor is changed from sensorless vector control using extended electromotive voltage to current drawing control without a shock. It is for switching.

制御装置における演算は、図2に示すように角速度ωで回転するγ,δ軸からなる直交回転座標系により行う。なお、図2において、d軸は回転子の磁極方向の軸、q軸はd軸と直交する方向の軸、θerrは位置推定誤差である。 The calculation in the control device is performed by an orthogonal rotation coordinate system composed of γ and δ axes rotating at an angular velocity ω 1 , as shown in FIG. In FIG. 2, the d-axis is the axis in the direction of the magnetic pole of the rotor, the q-axis is the axis in the direction orthogonal to the d-axis, and θ err is the position estimation error.

まず、この実施形態の電流引込制御時の動作を、制御装置の構成と共に説明する。
図1に戻って、切替スイッチ23の入力側を「S」に設定し、γ,δ軸の角速度ωを加算器124から出力される速度指令値ωに制御する。電気角演算器24は、角速度ωを積分してγ,δ軸の角度θを演算する。
First, the operation at the time of current drawing control of this embodiment will be described together with the configuration of the control device.
Returning to FIG. 1, the input side of the changeover switch 23 is set to “S 1 ”, and the angular velocity ω 1 of the γ and δ axes is controlled to the velocity command value ω * output from the adder 124. Electrical angle calculator 24, gamma integrates the angular velocity omega 1, calculates the angle theta 1 of the δ-axis.

一方、切替スイッチ21a,21bの入力側を「S」に設定する。第1の電流指令演算部41は、γ,δ軸電流指令値iγ ,iδ を数式1のように制御する。
[数1]
γ =Iapull
δ =0
ここで、
apull >0
On the other hand, the input sides of the changeover switches 21a and 21b are set to "S 1 ". The first current command calculation unit 41 controls the γ, δ axis current command values i γ * , i δ * as in Expression 1.
[Equation 1]
i γ * =I apull *
i δ * =0
here,
I apull * >0

ローパスフィルタ22a,22bは、γ,δ軸電流指令値iγ ,iδ の低周波成分iγf ,iδf をそれぞれ演算する。
電流座標変換器14は、u相電流検出器11u、w相電流検出器11wによりそれぞれ検出した相電流検出値i,iを、γ,δ軸の角度θに基づいてγ,δ軸電流検出値iγ,iδに座標変換する。この角度θは、回転子の位置推定値に相当する。
Low pass filter 22a, 22b is, gamma, [delta] -axis current value i γ *, i δ * of low-frequency component i .gamma.f *, calculates i delta] f *, respectively.
The current coordinate converter 14 calculates the phase current detection values i u and i w detected by the u-phase current detector 11 u and the w-phase current detector 11 w, respectively, based on the angle θ 1 of the γ and δ axes. The coordinates are converted into the detected current values i γ and i δ . This angle θ 1 corresponds to the estimated position value of the rotor.

減算器19a,19bにより、γ軸電流指令値iγf とγ軸電流検出値iγとの偏差、及び、δ軸電流指令値iδf とδ軸電流検出値iδとの偏差をそれぞれ演算し、電流調節器20は、上記偏差をそれぞれゼロにするようにγ軸電圧指令値vγ 、δ軸電圧指令値vδ を演算する。
これらの電圧指令値vγ ,vδ は、電圧座標変換器15により、γ,δ軸の角度θに基づいて相電圧指令値v ,v ,v に変換される。
The subtracters 19a and 19b respectively calculate the deviation between the γ-axis current command value i γf * and the γ-axis current detection value i γ and the deviation between the δ-axis current command value i δf * and the δ-axis current detection value i δ. Then, the current controller 20 calculates the γ-axis voltage command value v γ * and the δ-axis voltage command value v δ * so as to make the above deviations zero.
These voltage command values v γ * , v δ * are converted by the voltage coordinate converter 15 into phase voltage command values v u * , v v * , v w * based on the angle θ 1 of the γ and δ axes. It

三相交流電源50の電圧は整流回路60により直流電圧に変換され、インバータ等の電力変換器70に供給される。PWM回路13は、相電圧指令値v ,v ,v 、及び、入力電圧検出回路12により検出した電力変換器70の直流入力電圧Edcから、電力変換器70の出力電圧を相電圧指令値v ,v ,v に制御するための駆動信号(ゲート信号)を生成する。
電力変換器70はこのゲート信号に基づいて内部のIGBT等の半導体スイッチング素子を制御することにより、永久磁石形同期電動機(PMSM)80の端子電圧を相電圧指令値v ,v ,v に制御する。
The voltage of the three-phase AC power supply 50 is converted into a DC voltage by the rectifier circuit 60 and supplied to the power converter 70 such as an inverter. The PWM circuit 13 outputs the output voltage of the power converter 70 from the phase voltage command values v u * , v v * , v w * , and the DC input voltage E dc of the power converter 70 detected by the input voltage detection circuit 12. Drive signals (gate signals) for controlling the phase voltage command values v u * , v v * , v w * .
The power converter 70 controls a semiconductor switching element such as an internal IGBT based on the gate signal to change the terminal voltage of the permanent magnet synchronous motor (PMSM) 80 to the phase voltage command values v u * , v v * , Control to v w * .

上述した一連の制御により、振幅がIapull 、角速度が速度指令値ωに等しい電流ベクトルが発生し、電動機80の回転子が電流ベクトルに引き込まれて回転子速度ωを速度指令値ωに制御することができる。 By the series of controls described above, a current vector whose amplitude is I apull * and whose angular velocity is equal to the speed command value ω * is generated, and the rotor of the electric motor 80 is drawn into the current vector to change the rotor speed ω r to the speed command value ω Can be controlled to * .

次に、センサレスベクトル制御時の動作について説明する。以下では、電流引込制御と異なる部分を中心に説明し、重複する部分については省略する。
速度指令値ωと速度推定値ωrestとの偏差を減算器16により演算し、この偏差がゼロになるように速度調節器17がトルク指令値τを演算する。
Next, the operation during sensorless vector control will be described. In the following, a description will be given focusing on the parts different from the current drawing control, and overlapping parts will be omitted.
The subtractor 16 calculates the deviation between the speed command value ω * and the estimated speed value ω rest, and the speed controller 17 calculates the torque command value τ * so that the deviation becomes zero.

第2の電流指令演算部18は、トルク指令値τに基づき、トルク/電流が最大になる条件でd,q軸電流指令値i ,i を演算する。
切替スイッチ21a,21bの入力側を「S」に設定することにより、γ,δ軸電流指令値iγ ,iδ としてd,q軸電流指令値i ,i を設定する。
The second current command calculation unit 18 calculates the d and q axis current command values i d * , i q * based on the torque command value τ * under the condition that the torque/current is maximized.
By setting the input sides of the changeover switches 21a and 21b to "S 2 ", the d and q axis current command values i d * and i q * are set as the γ and δ axis current command values i γ * and i δ *. To do.

ローパスフィルタ22a,22b、減算器19a,19b、電流調節器20、電流座標変換器14、電圧座標変換器15の動作は、電流引込制御の場合と同じであり、これらの動作によって、γ,δ軸電流検出値iγ,iδが各指令値iγ ,iδ にそれぞれ一致するように相電圧指令値v ,v ,v が制御される。 The operations of the low-pass filters 22a and 22b, the subtractors 19a and 19b, the current regulator 20, the current coordinate converter 14, and the voltage coordinate converter 15 are the same as those in the current pull-in control. The phase voltage command values v u * , v v * , v w * are controlled so that the detected axial current values i γ , i δ match the command values i γ * , i δ * , respectively.

また、拡張誘起電圧演算器31は、γ,δ軸電圧指令値vγ ,vδ 、γ,δ軸電流検出値iγ,iδ、及び角速度ωを用いて、数式2により拡張誘起電圧を演算する。

Figure 0006726390
Further, the extended electromotive voltage calculator 31 uses the γ, δ-axis voltage command values v γ * , v δ * , the γ, δ-axis current detection values i γ , i δ , and the angular velocity ω 1 to perform expansion according to Equation 2. Calculate the induced voltage.
Figure 0006726390

なお、数式2の演算は、γ,δ軸電圧指令値vγ ,vδ の代わりに、図示されていない電圧検出回路を用いて電動機80の相電圧または線間電圧を測定し、これらの測定値及び角度(位置推定値)θから演算したγ軸電圧vγ、δ軸電圧vδを用いて行っても良い。 In addition, in the calculation of Expression 2, instead of the γ and δ axis voltage command values v γ * and v δ * , the phase voltage or line voltage of the electric motor 80 is measured using a voltage detection circuit (not shown), and these It is also possible to use the γ-axis voltage v γ and the δ-axis voltage v δ calculated from the measured value and the angle (position estimated value) θ 1 .

角度演算器32は、拡張誘起電圧の角度δeexestを数式3により演算する。

Figure 0006726390
The angle calculator 32 calculates the angle δ eexest of the extended electromotive voltage according to Expression 3.
Figure 0006726390

拡張誘起電圧の角度δeexestは位置推定誤差θerrに等しくなるので、位置推定誤差演算値θerrestを数式4とする。
[数4]
θerrest=δeexest
Since the angle δ eexest of the extended electromotive voltage becomes equal to the position estimation error θ err , the position estimation error calculation value θ errest is set to Expression 4.
[Equation 4]
θ errest = δ eexest

速度推定器33は、位置推定誤差演算値θerrestを入力とするPI調節器によって構成されており、速度推定値ωrestを数式5により演算する。

Figure 0006726390
The speed estimator 33 is configured by a PI adjuster that receives the position estimation error calculation value θ errest, and calculates the speed estimation value ω rest by Expression 5.
Figure 0006726390

切替スイッチ23の入力側を「S」に設定することにより、γ,δ軸の角速度ωを速度推定値ωrestに制御する。
電気角演算器24は、電流引込制御の場合と同様に、角速度ωを積分して回転子の位置推定値に相当する角度θを演算する。
By setting the input side of the changeover switch 23 to “S 2 ”, the angular velocity ω 1 of the γ and δ axes is controlled to the estimated velocity value ω rest .
Electrical angle calculator 24, like the case of the current pull-in control, calculates the angle theta 1 which integrates the angular velocity omega 1 corresponds to the position estimate of the rotor.

上述した拡張誘起電圧演算器31、角度演算器32、速度推定器33、電気角演算器24の動作により、γ,δ軸とd,q軸との角度差である位置推定誤差θerrを零に収束させることができ、回転子の速度及び磁極位置を正確に演算して電動機80のトルク及び速度を正確に制御することができる。 By the operations of the extended electromotive voltage calculator 31, the angle calculator 32, the speed estimator 33, and the electrical angle calculator 24 described above, the position estimation error θ err , which is the angular difference between the γ and δ axes and the d and q axes, is zero. The speed and magnetic pole position of the rotor can be accurately calculated to accurately control the torque and speed of the electric motor 80.

次に、センサレスベクトル制御から電流引込制御への切替処理の原理について、図3のベクトル図に基づいて説明する。
センサレスベクトル制御の場合、トルク/電流が最大になるように電流指令ベクトルiを制御するが、電流引込制御では、一般に、電流指令値を一定に制御し、電流指令ベクトルiの方向をγ軸と定義して制御演算を行う。
Next, the principle of the switching process from the sensorless vector control to the current drawing control will be described based on the vector diagram of FIG.
In the sensorless vector control, the current command vector i * is controlled so that the torque/current is maximized, but in the current pull-in control, generally, the current command value is controlled to be constant, and the direction of the current command vector i * is γ. Control calculation is performed by defining it as an axis.

そこで、センサレスベクトル制御から電流引込制御へ切替えるときに、切替前後でトルクが変化しないように電流指令値を演算し、d軸(磁極位置)と電流指令ベクトルiとの角度差δを演算し、この角度差δを使ってγ,δ軸の角度θを初期化する。これと同時に、電圧指令ベクトルvとローパスフィルタ22a,22bの出力である電流指令ベクトルi とが切替前後で連続になるように初期化する。 Therefore, when the sensorless vector control is switched to the current drawing control, the current command value is calculated so that the torque does not change before and after the switching, and the angle difference δ i between the d-axis (magnetic pole position) and the current command vector i * is calculated. Then, this angle difference δ i is used to initialize the angle θ 1 of the γ and δ axes. At the same time, the voltage command vector v * and the current command vector if * , which is the output of the low-pass filters 22a and 22b, are initialized to be continuous before and after the switching.

次に、センサレスベクトル制御から電流引込制御に切替える場合の演算処理について説明する。
図1の電流指令初期値演算器121は、センサレスベクトル制御から電流引込制御への切替直後に、電動機80のトルクをセンサレスベクトル制御時のトルク指令値τに制御するための角度差δを演算する。なお、電流指令初期値演算器121による演算内容は、後に詳述する。
Next, the calculation process when switching from sensorless vector control to current drawing control will be described.
Immediately after switching from the sensorless vector control to the current drawing control, the current command initial value calculator 121 of FIG. 1 determines the angle difference δ i for controlling the torque of the electric motor 80 to the torque command value τ * during the sensorless vector control. Calculate The details of the calculation by the current command initial value calculator 121 will be described later.

ここで、上記の角度差δから、電気角補正値δcompを数式6により演算する。
[数6]
δcomp=−δ
Here, the electrical angle correction value δ comp is calculated by the mathematical expression 6 from the above-mentioned angle difference δ i .
[Equation 6]
δ comp =−δ i

この電気角補正値δcompを用いて、電気角演算器24の出力である角度θを、数式7により初期化する。
[数7]
θ1(n)=θ1(n−1)−δcomp+ω
ここで、T:サンプル周期
なお、数式7において、添え字“n”はサンプル点を示し、“n”は今回値(切替後の値)、“n−1”は前回値(切替前の値)を示す。
Using this electrical angle correction value δ comp , the angle θ 1 which is the output of the electrical angle calculator 24 is initialized by the mathematical expression 7.
[Equation 7]
θ 1(n)1(n−1) −δ comp1 T s
Here, T s : sampling period In the formula 7, the subscript “n” indicates a sampling point, “n” is the current value (value after switching), and “n−1” is the previous value (before switching). Value).

数式7によってγ,δ軸の角度θをδcompだけ補正したことに伴い、ローパスフィルタ22a,22bの出力であるγ,δ軸電流指令値iγf ,iδf 、及び、電流調節器20の出力であるγ,δ軸電圧指令値vγ ,vδ を、それぞれ数式8,数式9により初期化する。

Figure 0006726390
Figure 0006726390
Since the angle θ 1 of the γ and δ axes is corrected by δ comp by Equation 7, the γ and δ axis current command values i γf * , i δf * , which are the outputs of the low-pass filters 22a and 22b, and the current regulator. The γ and δ axis voltage command values v γ * and v δ * which are the outputs of 20 are initialized by Equations 8 and 9, respectively.
Figure 0006726390
Figure 0006726390

数式8,9により、電動機80の電流と端子電圧とが急変しないように、γ,δ軸電流指令値iγf ,iδf 及びγ,δ軸電圧指令値vγ ,vδ を初期化することができる。 The γ and δ axis current command values i γf * , i δf * and the γ and δ axis voltage command values v γ * , v δ * are set so that the current and the terminal voltage of the electric motor 80 do not change suddenly by the mathematical expressions 8 and 9. Can be initialized.

角速度補償器123は、γ,δ軸の角速度ωが切替え時に急変しないように、電流引込制御開始から角速度補償値ω1compを演算し、加算器124は速度指令値ωと角速度補償値ω1compとを加算してγ,δ軸の角速度ωを演算する。
図4は、角速度補償器123の動作原理を示す図である。角速度補償値ω1compの初期値は、切替え直前に図1の減算器122により演算したγ,δ軸の角速度ω(速度推定値に等しい)と速度指令値ωとの偏差(ω−ω)とし、この角速度補償値ω1compを所定の変化率で零まで低減する。
The angular velocity compensator 123 calculates the angular velocity compensation value ω 1comp from the start of the current drawing control so that the angular velocity ω 1 of the γ and δ axes does not change suddenly at the time of switching, and the adder 124 calculates the velocity command value ω * and the angular velocity compensation value ω 1. 1comp is added to calculate the angular velocity ω 1 of the γ and δ axes.
FIG. 4 is a diagram showing the operating principle of the angular velocity compensator 123. The initial value of the angular velocity compensation value ω 1comp is the deviation (ω 1 −) between the angular velocity ω 1 (equal to the estimated velocity value) of the γ and δ axes calculated by the subtractor 122 in FIG. 1 immediately before switching and the velocity command value ω *. ω * ), and the angular velocity compensation value ω 1comp is reduced to zero at a predetermined change rate.

次に、電流指令初期値演算器121による演算内容について詳述する。
電流引込制御によって運転している場合、電動機80のトルクτと角度差δとの間には数式10の関係がある。
[数10]
τ=IaPULL Ψsinδ+IaPULL *2(L−L)sinδcosδ
ただし、Ψ:永久磁石磁束
角度差δは、上記のトルクτをトルク指令値τに置き換えた方程式の解から求める。数式10は非線形方程式であるため、角度差δは、ニュートンラプソン法や二分法などの数値計算によって求めれば良い。
Next, the contents of calculation by the current command initial value calculator 121 will be described in detail.
When operating by current drawing control, there is a relationship of Expression 10 between the torque τ of the electric motor 80 and the angle difference δ i .
[Equation 10]
τ=I aPULL * Ψ m sin δ i +I aPULL *2 (L d −L q )sin δ i cos δ i
However, Ψ m : Permanent magnet magnetic flux Angle difference δ i is obtained from the solution of the equation in which the torque τ is replaced with the torque command value τ * . Since the mathematical expression 10 is a non-linear equation, the angular difference δ i may be obtained by numerical calculation such as the Newton-Raphson method or the bisection method.

一例として、角度差δを二分法により求める方法について説明する。二分法は、ニュートンラプソン法などの他の計算方式と比べて演算を簡素化できる特徴がある。
図5は、数式10に示した角度差δとトルクτとの関係を示している。図5において、トルクτをトルク指令値τとする角度差δの解をδ(τ)とする。
As an example, a method of obtaining the angle difference δ i by the bisection method will be described. The bisection method has a feature that the operation can be simplified as compared with other calculation methods such as the Newton-Raphson method.
FIG. 5 shows the relationship between the angular difference δ i and the torque τ shown in Expression 10. In FIG. 5, the solution of the angle difference δ i where the torque τ is the torque command value τ * is δ i* ).

二分法を使って解を求める場合、解の区間を設定する必要がある。このため、角度差δの下限値及び上限値を、それぞれδimin,δimaxとする。電流引込制御を安定に行えるのは、トルクτが角度差δの増加関数となる場合である。そこで、角度差δの下限値δimin及び上限値δimaxは、角度差δiの区間[δimin,δimax]においてトルクτが増加関数になるように決める。 When the solution is obtained using the dichotomy, it is necessary to set the interval of the solution. Therefore, the lower limit value and the upper limit value of the angle difference δ i are set to δ imin and δ imax , respectively. The current pull-in control can be performed stably when the torque τ is an increasing function of the angle difference δ i . Therefore, the lower limit value [delta] imin and an upper limit [delta] imax of angular difference [delta] i is the angular difference δi interval [δ imin, δ imax] torque τ in decide to be an increasing function.

次に、角度差δの上限値δimaxの演算方法について説明する。図6は、角度差δの上限値δimaxの演算原理を説明するためのベクトル図である。
電流の振幅を電流指令値振幅IaPULL に制御し、角度差δが角度差上限値δimaxであるときの電流ベクトルを図6の電流ベクトル(2)とする。
ここで、電流の振幅Iとd,q軸電流i,iとの間には、数式11の関係がある。
[数11]
=i +i
数式11より、電流ベクトル(2)は、円軌跡I=IaPULL の上に制御される。
Next, a method of calculating the upper limit value δ imax of the angle difference δ i will be described. FIG. 6 is a vector diagram for explaining the calculation principle of the upper limit value δ imax of the angle difference δ i .
The current amplitude is controlled to the current command value amplitude I aPULL *, and the current vector when the angle difference δ i is the angle difference upper limit value δ imax is set as the current vector (2) in FIG. 6.
Here, there is a relationship of Expression 11 between the amplitude I a of the current and the d, q-axis currents i d , i q .
[Equation 11]
I a 2 = id 2 +i q 2
From Expression 11, the current vector (2) is controlled on the circular locus I a =I aPULL * .

一方、電動機のトルクτとd,q軸電流i,iとの間には、数式12の関係がある。
[数12]
τ=Ψ+(L−L)i
数式12により、トルクτがトルク指令値τに等しい電流ベクトルは、放物線の軌跡τ=τである。トルクτがトルク指令値τよりも大きくなるための条件から、電流ベクトル(2)の動作点は、τ=τの放物線よりも上側にある必要がある。
On the other hand, there is a relationship of Expression 12 between the torque τ of the electric motor and the d and q axis currents i d and i q .
[Equation 12]
τ=Ψ m i q +(L d −L q ) i d i q
From Expression 12, the current vector in which the torque τ is equal to the torque command value τ * is the parabolic locus τ=τ * . From the condition that the torque τ becomes larger than the torque command value τ * , the operating point of the current vector (2) needs to be above the parabola of τ=τ * .

更に、トルク/電流を最大とするd,q軸電流i,iは、数式13の関係にある。

Figure 0006726390
Furthermore, the d and q axis currents i d and i q that maximize the torque/current have the relationship of Expression 13.
Figure 0006726390

数式13より、トルク/電流が最大になる電流ベクトルの軌跡は、図6における原点を通る放物線の軌跡「トルク/電流最大条件」となる。トルクτが角度差δの増加関数となるための条件から、電流ベクトル(2)は、トルク/電流最大条件の放物線よりも左側にある必要がある。 From Expression 13, the locus of the current vector that maximizes the torque/current is the parabolic locus “torque/current maximum condition” that passes through the origin in FIG. 6. Due to the condition that the torque τ is an increasing function of the angle difference δ i , the current vector (2) needs to be on the left side of the parabola of the torque/current maximum condition.

角度差上限値δimaxは、以上の条件を満たすように、以下に述べる手順によって演算する。
まず、トルク/電流最大条件でトルクτをトルク指令値τに制御したときの電流ベクトルを、電流ベクトル(1)とする。この電流ベクトル(1)は、τ=τの放物線とトルク/電流最大条件の放物線の交点に達する。電流ベクトル(1)は、トルク/電流最大条件であるので、電流ベクトル(1)の振幅は、電流指令値振幅IaPULL よりも小さくなる。
The angle difference upper limit value δ imax is calculated by the procedure described below so as to satisfy the above conditions.
First, the current vector when the torque τ is controlled to the torque command value τ * under the torque/current maximum condition is set as the current vector (1). This current vector (1) reaches the intersection of the parabola of τ=τ * and the parabola of the torque/current maximum condition. Since the current vector (1) is the torque/current maximum condition, the amplitude of the current vector (1) becomes smaller than the current command value amplitude I aPULL * .

次に、電流ベクトル(2)を、電流振幅が電流指令値振幅IaPULL に等しく、d軸電流が電流ベクトル(1)のd軸電流idMTPAに等しい動作点に選ぶ。具体的には、数式14により演算する。

Figure 0006726390
数式14におけるIdMTPAは、トルク指令値τからテーブルを用いて演算する。 Next, the current vector (2) is selected as an operating point where the current amplitude is equal to the current command value amplitude I aPULL * and the d-axis current is equal to the d-axis current i dMTPA of the current vector (1). Specifically, the calculation is performed using Equation 14.
Figure 0006726390
I dMTPA in Expression 14 is calculated from the torque command value τ * using a table.

角度差上限値δimaxは、d軸と電流ベクトル(2)との角度差から、数式15により演算する。

Figure 0006726390
The angle difference upper limit value δ imax is calculated by Expression 15 from the angle difference between the d-axis and the current vector (2).
Figure 0006726390

電流ベクトル(2)の動作点は、τ=τの放物線よりも上側であるため、トルクτは、トルク指令値τよりも大きくなる。また、トルク/電流最大条件の放物線の左側であるので、トルクτは角度差δの増加関数となる。
これらのことから、d軸と電流ベクトル(2)との角度差から演算した角度差上限値δimaxは、必要とする条件を満たす。
Since the operating point of the current vector (2) is above the parabola of τ=τ * , the torque τ is larger than the torque command value τ * . Further, since it is on the left side of the parabola of the torque/current maximum condition, the torque τ is an increasing function of the angle difference δ i .
From these, the angle difference upper limit value δ imax calculated from the angle difference between the d-axis and the current vector (2) satisfies the necessary condition.

次に、角度差δの下限値δiminの演算方法を説明する。図7は、角度差δの下限値δiminの演算原理を説明するためのベクトル図である。
数式12より、埋込磁石永久磁石形同期電動機(IPMSM)のように、d軸インダクタンスLがq軸インダクタンスLよりも小さい電動機の場合、d軸電流iが大きくなるとトルクτがq軸電流iの増加関数にならなくなる。トルクτがq軸電流iの増加関数であるためのd軸電流iの上限値をIdlimpとすると、Idlimpは数式16によって演算することができる。

Figure 0006726390
特に、d軸電流iがIdlimpである場合、トルクτはq軸電流iによらず、零になる。 Next, a method of calculating the lower limit value δ imin of the angle difference δ i will be described. FIG. 7 is a vector diagram for explaining the calculation principle of the lower limit value δ imin of the angle difference δ i .
From Expression 12, in the case of an electric motor having a d-axis inductance L d smaller than the q-axis inductance L q , such as an embedded magnet permanent magnet synchronous motor (IPMSM), the torque τ becomes q-axis when the d-axis current i d increases. It does not become an increasing function of the current i q . Since the torque τ is an increasing function of the q-axis current i q , the upper limit value of the d-axis current i d is I dlimp , and I dlimp can be calculated by Expression 16.
Figure 0006726390
In particular, when the d-axis current i d is I dlimp , the torque τ becomes zero regardless of the q-axis current i q .

以上のことから、電流引込制御の場合には、電流指令値振幅IaPULL が大きく、角度差δが小さい軽負荷時にd軸電流iが上限値Idlimよりも大きくなりやすい。そこで、角度差の下限値δiminは、このことを考慮して設計する。 From the above, in the case of the current drawing control, the d-axis current i d is likely to be larger than the upper limit value I dlim at the time of a light load in which the current command value amplitude I aPULL * is large and the angle difference δ i is small. Therefore, the lower limit value δ imin of the angle difference is designed in consideration of this.

まず、電流指令値振幅IaPULL がd軸電流iの上限値Idlimpより小さいときは、角度差δが小さいときも、トルクτは角度差δの増加関数になる。そこで、数式17に示すように、角度差の下限値δiminを0とする。
[数17]
δimin=0
First, when the current command value amplitude I aPULL * is smaller than the upper limit value I dlimp of the d-axis current i d , the torque τ becomes an increasing function of the angle difference δ i even when the angle difference δ i is small. Therefore, as shown in Expression 17, the lower limit value δ imin of the angle difference is set to 0.
[Equation 17]
δ imin =0

一方、電流指令値振幅IaPULL がd軸電流iの上限値Idlimp以上の場合、トルクτを零とする電流ベクトルは、図7に示すように、振幅が電流指令値振幅IaPULL に等しく、d軸電流がその上限値Idlimpに等しい電流ベクトル(3)となる。この電流ベクトル(3)のd,q軸電流は、数式18となる。

Figure 0006726390
On the other hand, when the current command value amplitude I aPULL * is greater than or equal to the upper limit value I dlimp of the d-axis current i d , the current vector that sets the torque τ to zero has a current command value amplitude I aPULL * And the d-axis current is a current vector (3) equal to its upper limit value I dlimp . The d and q-axis currents of this current vector (3) are given by Equation 18.
Figure 0006726390

角度差の下限値δiminは、d軸と電流ベクトル(3)との角度差から、数式19により演算する。

Figure 0006726390
The lower limit value δ imin of the angle difference is calculated by Expression 19 from the angle difference between the d-axis and the current vector (3).
Figure 0006726390

以上説明した演算により、トルクτが増加関数になり、かつ、トルク指令値τが存在する角度差δの区間[δimin,δimax]を演算することができる。この結果、数式10に示した角度差δとトルクτとの関係式、及び、角度差δの区間[δimin,δimax]から、トルクτをトルク指令値τに制御する角度差δi(τ*)を二分法によって演算することが可能である。 By the calculation described above, the section [δ imin , δ imax ] of the angle difference δ i in which the torque τ becomes an increasing function and the torque command value τ * exists can be calculated. As a result, from the relational expression between the angle difference δ i and the torque τ shown in Formula 10, and the section [δ imin , δ imax ] of the angle difference δ i , the angle difference for controlling the torque τ to the torque command value τ *. It is possible to calculate δ i (τ*) by the bisection method.

11u:u相電流検出回路
11w:w相電流検出回路
12:入力電圧検出回路
13:PWM回路
14:電流座標変換器
15:電圧座標変換器
16,19a,19b:減算器
17:速度調節器
18:電流指令演算部
20:電流調節器
21a,21b,23:切替スイッチ
22a,22b:ローパスフィルタ
24:電気角演算器
31:拡張誘起電圧演算器
32:角度演算器
33:速度推定器
41:電流指令演算部
50:三相交流電源
60:整流回路
70:電流変換器
80:永久磁石形同期電動機
121:電流指令初期値演算器
122:減算器
123:角速度補償器
124:加算器
11u: u-phase current detection circuit 11w: w-phase current detection circuit 12: input voltage detection circuit 13: PWM circuit 14: current coordinate converter 15: voltage coordinate converter 16, 19a, 19b: subtractor 17: speed controller 18 : Current command calculator 20: Current controllers 21a, 21b, 23: Changeover switches 22a, 22b: Low-pass filter 24: Electric angle calculator 31: Extended induced voltage calculator 32: Angle calculator 33: Speed estimator 41: Current Command calculator 50: Three-phase AC power supply 60: Rectifier circuit 70: Current converter 80: Permanent magnet type synchronous motor 121: Current command initial value calculator 122: Subtractor 123: Angular velocity compensator 124: Adder

Claims (1)

磁極位置検出器を持たない永久磁石形同期電動機を電力変換器により運転するための制御装置において、
前記電動機の電流及び端子電圧をベクトルとしてとらえ、
前記電動機の電流指令値の角速度を速度指令値に制御する第1の運転モードと、
前記電動機の電流及び端子電圧から演算した回転子の磁極位置推定値及び速度推定値を用いて前記電動機の速度を速度指令値に制御する第2の運転モードと、を有し、
前記第2の運転モードから前記第1の運転モードへ切替えるときに、
前記電動機のトルク指令値と前記電流指令値の振幅とから、前記電流指令値と回転子の磁極位置との角度差を演算する手段と、
前記角度差から電気角補正値を求める手段と、
前記電気角補正値を用いて前記電流指令値の角度を初期化する手段と、
前記電気角補正値を用いて、前記電動機の端子電圧が急変しないように前記電動機の端子電圧指令値を初期化する手段と、
前記速度指令値と前記速度推定値との偏差から求めた角速度補償値を用いて、前記電流指令値の角速度が急変しないように前記電流指令値の角速度を補正する手段と、
を備え
前記電流指令値と前記磁極位置との角度差を演算する手段は、
前記電流指令値の振幅から、トルクが前記角度差の増加関数となる、前記角度差の区間を演算する手段と、
前記角度差の区間から、前記トルク指令値に前記電動機のトルクが一致するように前記角度差を演算する手段と、
を有することを特徴とする永久磁石形同期電動機の制御装置。
In a control device for operating a permanent magnet type synchronous motor having no magnetic pole position detector by a power converter,
Capturing the current and terminal voltage of the motor as a vector,
A first operation mode for controlling an angular velocity of a current command value of the electric motor to a speed command value;
A second operation mode in which the speed of the electric motor is controlled to a speed command value using an estimated magnetic pole position value and an estimated speed value of the rotor calculated from the electric current and the terminal voltage of the electric motor,
When switching from the second operation mode to the first operation mode,
Means for calculating an angular difference between the current command value and the magnetic pole position of the rotor from the torque command value of the electric motor and the amplitude of the current command value;
Means for obtaining an electrical angle correction value from the angle difference,
Means for initializing the angle of the current command value using the electrical angle correction value;
Means for initializing a terminal voltage command value of the electric motor so that the terminal voltage of the electric motor does not change suddenly using the electrical angle correction value;
A means for correcting the angular velocity of the current command value so that the angular velocity of the current command value does not change suddenly using the angular velocity compensation value obtained from the deviation between the speed command value and the estimated speed value;
Equipped with
Means for calculating the angle difference between the current command value and the magnetic pole position,
Means for calculating a section of the angular difference, where torque is an increasing function of the angular difference, from the amplitude of the current command value;
Means for calculating the angular difference so that the torque of the electric motor matches the torque command value from the angular difference section;
Controller for a permanent magnet synchronous motor and having a.
JP2016133789A 2016-07-06 2016-07-06 Controller for permanent magnet type synchronous motor Active JP6726390B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016133789A JP6726390B2 (en) 2016-07-06 2016-07-06 Controller for permanent magnet type synchronous motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016133789A JP6726390B2 (en) 2016-07-06 2016-07-06 Controller for permanent magnet type synchronous motor

Publications (2)

Publication Number Publication Date
JP2018007473A JP2018007473A (en) 2018-01-11
JP6726390B2 true JP6726390B2 (en) 2020-07-22

Family

ID=60948209

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016133789A Active JP6726390B2 (en) 2016-07-06 2016-07-06 Controller for permanent magnet type synchronous motor

Country Status (1)

Country Link
JP (1) JP6726390B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6896174B2 (en) * 2018-07-18 2021-06-30 三菱電機株式会社 Rotating machine control device, refrigerant compressor, refrigeration cycle device and air conditioner
TWI702787B (en) * 2019-04-15 2020-08-21 楊逸群 Control device for three-phase permanent magnet motor
CN110729942B (en) * 2019-11-04 2020-11-03 广州新科佳都科技有限公司 Servo control method and device for door leaf module, electronic equipment and storage medium
JP7390914B2 (en) * 2020-02-04 2023-12-04 ニデックインスツルメンツ株式会社 Motor control method, motor drive device, industrial robot control method, and industrial robot
JP7428527B2 (en) * 2020-02-04 2024-02-06 ニデックインスツルメンツ株式会社 Motor control method, motor drive device, industrial robot control method, and industrial robot
JP7556301B2 (en) 2021-02-04 2024-09-26 株式会社明電舎 POWER CONVERSION SYSTEM AND METHOD FOR CONTROLLING POWER CONVERSION SYSTEM

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5277724B2 (en) * 2008-05-23 2013-08-28 富士電機株式会社 Control device for permanent magnet type synchronous motor

Also Published As

Publication number Publication date
JP2018007473A (en) 2018-01-11

Similar Documents

Publication Publication Date Title
JP6726390B2 (en) Controller for permanent magnet type synchronous motor
JP5155344B2 (en) Electric motor magnetic pole position estimation device
KR102108911B1 (en) Drive system and inverter device
JP6367332B2 (en) Inverter control device and motor drive system
WO2011145334A1 (en) Control device without a rotation sensor
US20140049202A1 (en) Motor control apparatus and motor control method
JP2009142116A (en) Position sensorless controller of permanent magnetic motor
US10696141B2 (en) Synchronous motor control device and method of controlling synchronous motor
JP5277724B2 (en) Control device for permanent magnet type synchronous motor
JP2015223023A (en) Controller of synchronous motor
JP2008278595A (en) Control device for permanent magnet type synchronous motor
JPWO2016121237A1 (en) Inverter control device and motor drive system
JP5709932B2 (en) Synchronous machine controller
JP5428202B2 (en) Control device for permanent magnet type synchronous motor
JP2008220169A (en) Motor controller
JP2006158046A (en) Sensorless control method and apparatus of ac electric motor
KR102409792B1 (en) Control device of permanent magnet synchronization electric motor, microcomputer, electric motor system, and driving method of permanent magnet synchronization electric motor
JP6590196B2 (en) Power converter
JP6565484B2 (en) Power converter
JP5332305B2 (en) Control device for permanent magnet type synchronous motor
JP6680104B2 (en) Motor control device and control method
JP4038412B2 (en) Vector control inverter device
JP6497584B2 (en) Control device for permanent magnet type synchronous motor
US20140152206A1 (en) Control device for alternating current rotating machine and electric power steering device equipped with control device for alternating current rotating machine
JP2021022965A (en) Driving device of induction motor, driving method, and electric vehicle

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190613

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200311

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200311

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200403

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200610

R150 Certificate of patent or registration of utility model

Ref document number: 6726390

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: R3D02

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250