Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP5169014B2 - Phase estimation method for AC motors - Google Patents
[go: Go Back, main page]

JP5169014B2 - Phase estimation method for AC motors - Google Patents

Phase estimation method for AC motors Download PDF

Info

Publication number
JP5169014B2
JP5169014B2 JP2007125010A JP2007125010A JP5169014B2 JP 5169014 B2 JP5169014 B2 JP 5169014B2 JP 2007125010 A JP2007125010 A JP 2007125010A JP 2007125010 A JP2007125010 A JP 2007125010A JP 5169014 B2 JP5169014 B2 JP 5169014B2
Authority
JP
Japan
Prior art keywords
rotor
magnetic flux
phase
rotor magnetic
value
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
JP2007125010A
Other languages
Japanese (ja)
Other versions
JP2008263764A (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.)
Nidec Corp
Original Assignee
Nidec Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nidec Corp filed Critical Nidec Corp
Priority to JP2007125010A priority Critical patent/JP5169014B2/en
Publication of JP2008263764A publication Critical patent/JP2008263764A/en
Application granted granted Critical
Publication of JP5169014B2 publication Critical patent/JP5169014B2/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 generates a rotor phase estimation value or a rotor magnetic flux phase estimation value necessary for sensorless drive control of an AC motor including a synchronous motor such as a permanent magnet synchronous motor and a synchronous reluctance motor and an induction motor. It relates to a phase estimation method.

電源周波数と同一の周波数で同期回転する交流電動機(永久磁石同期電動機、同期リラクタンス電動機などの同期電動機)をして高い制御性能を発揮せしめるには、回転子位相(N極位相、負突極位相など)が必要である。また、回転子磁束を有する永久磁石同期電動機、誘導電動機などの交流電動機をして高い制御性能を発揮せしめるには、回転子磁束の位相が必要である。このため、エンコーダ、レゾルバーと言った位置速度センサを回転子に装着し、これら位相を直接的あるいは間接的に得ていた。しかし、位置速度センサの装着は、信頼性、コストなどの点で問題があり、電動機駆動用の固定子電圧、電流の相当値を用いて回転子位相を推定することが行なわれている。当業者には周知のように、永久磁石同期電動機のように回転子磁束を有する同期電動機においては、回転子位相と回転子磁束位相とは、同一と捕らえることができる。The rotor phase (N pole phase, negative salient pole phase) is required to achieve high control performance by using AC motors (synchronous motors such as permanent magnet synchronous motors and synchronous reluctance motors) that rotate synchronously at the same frequency as the power supply frequency. Etc.) is necessary. In addition, the phase of the rotor magnetic flux is necessary for an AC motor such as a permanent magnet synchronous motor or an induction motor having a rotor magnetic flux to exhibit high control performance. For this reason, position and velocity sensors such as encoders and resolvers are mounted on the rotor, and these phases are obtained directly or indirectly. However, the mounting of the position / speed sensor is problematic in terms of reliability, cost, and the like, and the rotor phase is estimated using the equivalent values of the stator voltage and current for driving the motor. As is well known to those skilled in the art, in a synchronous motor having a rotor magnetic flux such as a permanent magnet synchronous motor, the rotor phase and the rotor magnetic flux phase can be regarded as the same.

交流電動機の合理的な駆動制御には、同期電動機にあっては回転子位相が、誘導電動機にあっては回転子磁束位相が必要である。位置速度センサレスによる駆動では、これら位相を推定する必要がある。以下に、永久磁石同期電動機、同期リラクタンス電動機、誘導電動機を例にとり、最も伝統的な位相推定法を紹介する。For rational drive control of an AC motor, a rotor phase is required for a synchronous motor, and a rotor magnetic flux phase is required for an induction motor. In driving without a position / velocity sensor, it is necessary to estimate these phases. In the following, the most traditional phase estimation method will be introduced by taking permanent magnet synchronous motor, synchronous reluctance motor, and induction motor as examples.

非突極の永久磁石同期電動機の場合、固定子の電圧、電流と回転子磁束を関係を記述した数学モデル(通常、電圧モデルと呼ばれる)は、固定αβ座標系上では、次式となる。

Figure 0005169014
Figure 0005169014
Figure 0005169014
In the case of a non-salient permanent magnet synchronous motor, a mathematical model describing the relationship between the stator voltage and current and the rotor magnetic flux (usually called a voltage model) is expressed by the following equation on the fixed αβ coordinate system.
Figure 0005169014
Figure 0005169014
Figure 0005169014

ここに、2x1ベクトルν、iは、それぞれ固定子の電圧、電流を意味している。また、2x1ベクトルφは固定子反作用磁束を、2x1ベクトルφは回転子永久磁石に起因する回転子磁束を意味している。Rは固定子巻線の抵抗であり、Lは固定子(d軸)インダクタンスであり、sは微分演算子d/dtである。永久磁石同期電動機の場合には、回転子位相θαは回転子磁束φの位相と同一であり、本位相は、(2)、(3)式が示しているように回転子磁束に含まれる。従って、回転子磁束が推定できれば、回転子磁束推定値より直ちに回転子位相推定値を得ることができる。Here, 2 × 1 vectors ν 1 and i 1 mean the voltage and current of the stator, respectively. The 2 × 1 vector φ i means the stator reaction magnetic flux, and the 2 × 1 vector φ m means the rotor magnetic flux caused by the rotor permanent magnet. R 1 is the resistance of the stator winding, L d is the stator (d-axis) inductance, and s is the differential operator d / dt. In the case of a permanent magnet synchronous motor, the rotor phase θ α is the same as the phase of the rotor magnetic flux φ m , and this phase is included in the rotor magnetic flux as shown in equations (2) and (3). It is. Therefore, if the rotor magnetic flux can be estimated, the rotor phase estimated value can be obtained immediately from the rotor magnetic flux estimated value.

最も伝統的な回転子磁束推定法は、誘起電圧相当値を固定子の電圧・電流を用いて生成し、生成値を、次のように近似積分するものである。

Figure 0005169014
Figure 0005169014
式右辺の四角括弧の中が、誘起電圧相当値を意味している。The most traditional method of estimating the rotor magnetic flux is to generate an induced voltage equivalent value by using the voltage / current of the stator, and approximate integration of the generated value as follows.
Figure 0005169014
Figure 0005169014
The square brackets on the right side of the equation mean the induced voltage equivalent value.

誘起電圧相当値に対する処理は、原理的には、−π/2(rad)の位相遅れをもつ純粋積分とする必要がある。しかし、実際には、電圧、電流、パラメータ等を含む種々の誤差が回転子磁束推定値(特にこの位相)に外乱として作用し、この緩和のため純粋積分は利用できず、影響緩和の観点から、電動機の回転速度を考慮の上、大き目の一定なフィルタ係数ωを持つ近似積分が行なわれてきた。近似積分において、理論上−π/2(rad)の位相遅れを持たせるには、被処理信号の周波数は無限大としなければならない。概ね−π/2(rad)の位相遅れを持たせるには、被処理信号の周波数は、一定フィルタ係数の10倍以上でなくてならない。従って、低〜中速域駆動時での推定において、大き目の一定なフィルタ係数をもつ近似積分による場合には、適切な−π/2(rad)の位相遅れが得られず、位相遅れが不足する。位相遅れの大きな不足は、回転子位相推定値における位相進み方向への位相の大きな狂いを意味する。位相推定法における位相推定値の位相狂いは、致命的な欠陥である。位相推定値における位相進みを低減するには、上述の一定フィルタ係数を小さくすればよいが、この場合、上に説明した外乱影響緩和上の問題が顕在化する。実際には、本推定法は、位相推定値の位相狂いと外乱影響緩和との妥協を図りつつ、駆動速度域を限定することにより、限定的に利用されてきた。In principle, the processing for the induced voltage equivalent value needs to be a pure integration having a phase delay of −π / 2 (rad). However, in practice, various errors including voltage, current, parameters, etc. act as disturbances on the rotor flux estimate (especially this phase), and pure integration cannot be used for this relaxation. Considering the rotational speed of the motor, approximate integration with a large constant filter coefficient ω c has been performed. In the approximate integration, in order to theoretically have a phase delay of −π / 2 (rad), the frequency of the signal to be processed must be infinite. In order to have a phase delay of approximately −π / 2 (rad), the frequency of the signal to be processed must be at least 10 times the constant filter coefficient. Therefore, in the estimation at the time of driving in the low to medium speed range, if approximate integration with a large constant filter coefficient is used, an appropriate phase delay of −π / 2 (rad) cannot be obtained, and the phase delay is insufficient. To do. A large shortage of phase delay means a large phase shift in the phase advance direction in the rotor phase estimation value. The phase error of the phase estimation value in the phase estimation method is a fatal defect. In order to reduce the phase advance in the phase estimation value, the above-described constant filter coefficient may be reduced. However, in this case, the above-described problem of mitigating the disturbance influence becomes obvious. Actually, this estimation method has been used in a limited manner by limiting the driving speed range while compromising the phase deviation of the phase estimation value and the disturbance influence mitigation.

突極特性をもつ同期リラクタンス電動機の場合、固定子の関係を記述した電圧モデルは、固定αβ座標系上では、次式となる。

Figure 0005169014
Figure 0005169014
In the case of a synchronous reluctance motor having salient pole characteristics, a voltage model describing the relationship of the stator is expressed by the following equation on the fixed αβ coordinate system.
Figure 0005169014
Figure 0005169014

ここに、L,Lは固定子dqインダクタンスと次の関係を有する同相インダクタンス、鏡相インダクタンスである。

Figure 0005169014
また、i、iは固定子のd軸、q軸電流であり、ω2nは回転子の電気速度である。また、Jは次式で定義された2x2交代行列である。
Figure 0005169014
Here, L i and L m are in-phase inductance and mirror-phase inductance having the following relationship with the stator dq inductance.
Figure 0005169014
Further, i d and i q are the d-axis and q-axis currents of the stator, and ω 2n is the electric speed of the rotor. J is a 2 × 2 alternating matrix defined by the following equation.
Figure 0005169014

同期リラクタンス電動機においては、(6)式の磁束φは、固定子反作用磁束の一部である。しかし、(5)、(6)式と(1)、(2)式との比較より明瞭のように、(6)式の磁束は、モデル上では、(2)式の回転子磁束と等価に扱うことができる。本明細書では、本等価性を考慮し、以降では、(6)式の磁束を擬似回転子磁束、または、簡単に回転子磁束と呼称する。In synchronous reluctance motor, (6) the magnetic flux phi m of Formula is part of the stator reaction magnetic flux. However, as is clear from the comparison between Equations (5) and (6) and Equations (1) and (2), the magnetic flux of Equation (6) is equivalent to the rotor magnetic flux of Equation (2) on the model. Can be handled. In the present specification, in consideration of this equivalence, hereinafter, the magnetic flux of the formula (6) is referred to as a pseudo rotor magnetic flux or simply a rotor magnetic flux.

(5)、(6)式におけるθαは、回転子位相(回転子の負突極または正突極の位相)を意味している。すなわち、設計者がd軸を負突極に選定すれば負突極の位相を、d軸を正突極に選定すれば正突極の位相を意味する。何れの場合も、設計者が選定したd軸の位相を意味する。(5)、(6)式と(1)、(2)式との比較より明白なように、同期リラクタンス電動機の電圧モデルは、永久磁石同期電動機の電圧モデルと高い類似性を有している。本類似性より理解されるように、同期リラクタンス電動機の回転子位相推定に関しては、永久磁石同期電動機のための(4)式の関係が利用可能である。(4)式の関係を利用する限りにおいては、位相推定値における位相の狂いの問題も、同様に発生する。なお、同期リラクタンス電動機の場合には、(4)式の右辺括弧に対応する値(すなわち、被積分値あるいは被近似積分値)は、拡張誘起電圧相当値となる。(5), which means the theta alpha in equation (6), the rotor phase (negative salient pole or positive salient pole of the rotor phase). That is, if the designer selects the d-axis as a negative salient pole, it means the phase of the negative salient pole, and if the designer selects the d-axis as a positive salient pole, it means the phase of the positive salient pole. In either case, it means the d-axis phase selected by the designer. As is clear from the comparison between the equations (5) and (6) and the equations (1) and (2), the voltage model of the synchronous reluctance motor has a high similarity to the voltage model of the permanent magnet synchronous motor. . As understood from this similarity, the relationship of the equation (4) for the permanent magnet synchronous motor can be used for the rotor phase estimation of the synchronous reluctance motor. As long as the relationship of equation (4) is used, the problem of phase shift in the phase estimation value also occurs. In the case of a synchronous reluctance motor, the value corresponding to the right parenthesis in equation (4) (that is, the integrated value or approximated integrated value) is an extended induced voltage equivalent value.

誘導電動機の固定子の関係を記述した電圧モデルは、固定αβ座標系上では、次式となる。

Figure 0005169014
ここに、φは誘導電動機における回転子磁束であり、φは固定子反作用磁束である。L1tは固定子の総合漏れインダクタンスである。誘導電動機のセンサレス駆動においては、回転子磁束の位相を把握する必要があるが、当業者には周知にように、本把握は回転子磁束の推定を通じ直ちに行うことができる。(9)式と(1)、(2)式との比較より明白なように、誘導電動機の電圧モデルは、永久磁石同期電動機の電圧モデル(特に、非突極の電圧モデル)と高い類似性を有している。本類似性より理解されるように、誘導電動機の回転子磁束位相推定に関しては、永久磁石同期電動機のための(4)式の関係が利用可能である。
(4)式の関係を利用する限りは、位相推定値における位相の狂いの問題も、同様に発生する。なお、誘導電動電動機の場合には、(4)式の右辺括弧に対応する値(すなわち、被積分値あるいは被近似積分値)は、誘起電圧相当値となる。The voltage model describing the relationship of the induction motor stator is expressed by the following equation on the fixed αβ coordinate system.
Figure 0005169014
Here, φ m is the rotor magnetic flux in the induction motor, and φ i is the stator reaction magnetic flux. L 1t is the total leakage inductance of the stator. In the sensorless drive of the induction motor, it is necessary to grasp the phase of the rotor magnetic flux, but as is well known to those skilled in the art, this grasp can be immediately performed by estimating the rotor magnetic flux. As is clear from the comparison between the expression (9) and the expressions (1) and (2), the voltage model of the induction motor is highly similar to the voltage model of the permanent magnet synchronous motor (particularly, the voltage model of a non-salient pole). have. As understood from this similarity, the relationship of the formula (4) for the permanent magnet synchronous motor can be used for the rotor magnetic flux phase estimation of the induction motor.
As long as the relationship of equation (4) is used, the problem of phase shift in the phase estimation value also occurs. In the case of an induction motor, the value corresponding to the right parenthesis in equation (4) (that is, the integrated value or approximated integrated value) is an induced voltage equivalent value.

位相推定値における位相の狂いは、致命的な欠陥である。致命的欠陥をなおすには、位相の狂いを補正する以外にはない。位相補正は、制御器設計でよく利用されているように、動的な位相補正器を利用するのが一般である。この種の動的位相補正器は、「ある周波数帯域で概ね位相進み」と言ったような概略的な位相補正には有用であるが、回転子位相推定あるいは回転子磁束位相推定のような精密な位相推定においては、不向きである。事実、この種の位相推定法においては、動的位相補正器を利用した位相補正法の活用は、殆んど行われていないようである。A phase shift in the phase estimate is a fatal defect. The only way to fix a fatal defect is to correct the phase error. Phase correction generally uses a dynamic phase corrector, as is often used in controller design. This type of dynamic phase compensator is useful for rough phase correction such as “generally phase advance in a certain frequency band”, but it is not as precise as rotor phase estimation or rotor flux phase estimation. It is not suitable for accurate phase estimation. In fact, in this type of phase estimation method, it seems that the phase correction method using a dynamic phase corrector is hardly utilized.

発明が解決しようとする課題Problems to be solved by the invention

本発明は上記背景の下になされたものであり、その目的は、同期電動機、誘導電動機を含む交流電動機のセンサレス駆動制御において必要とされる回転子位相あるいは回転子磁束位相を、位相の狂いが出ない状態でかつ外乱影響を抑圧した状態で推定できる新規な位相推定法を提供することである。The present invention has been made under the above-mentioned background, and its purpose is to change the rotor phase or rotor magnetic flux phase required for sensorless drive control of an AC motor including a synchronous motor and an induction motor. It is to provide a novel phase estimation method that can be estimated in a state that does not appear and in which a disturbance influence is suppressed.

課題を解決するための手段Means for solving the problem

上記目的を達成するために、請求項1に記載の発明は、回転子が固定子電源周波数と同期して回転する交流電動機の回転子の位相推定法であって、固定子の電圧と電流との相当値を少なくとも用いて生成した、回転子の回転に起因する誘起電圧または拡張誘起電圧の相当値を、回転子の回転に応じて通過帯域幅を変化させる速応帯域特性と回転子の非ゼロ回転速度に対応した周波数で-π/2(rad)の位相遅れ特性と高域で−40(dB/dec)以上の減衰特性とをもつ高次フィルタで処理し、回転子の非ゼロ回転速度に対応した周波数の基本波成分を主とする被フィルタ処理信号を生成し、被フィルタ処理信号から回転子の位相推定値を生成することを特徴とする。
In order to achieve the above object, the invention described in claim 1 is a method for estimating a phase of a rotor of an AC motor in which a rotor rotates in synchronization with a stator power supply frequency. the corresponding value was at least using generated, and the equivalent value of the induced voltage or the extended electromotive force caused by the rotation of the rotor, speed応帯frequency response and the rotor of the non-changing the pass bandwidth in accordance with the rotation of the rotor Processing with a high-order filter having a phase lag characteristic of -π / 2 (rad) at a frequency corresponding to zero rotation speed and an attenuation characteristic of -40 (dB / dec) or higher in the high band, and non-zero rotation of the rotor A filtered signal mainly including a fundamental wave component having a frequency corresponding to the speed is generated, and a phase estimation value of the rotor is generated from the filtered signal.

請求項2に記載の発明は、回転子磁束を有する交流電動機の回転子磁束の位相推定法であって、固定子の電圧と電流との相当値を少なくとも用いて生成した、回転子磁束の回転に起因する誘起電圧または拡張誘起電圧の相当値を、回転子磁束の回転に応じて通過帯域幅を変化させる速応帯域特性と回転子磁束の非ゼロ回転速度に対応した周波数で-π/2(rad)の位相遅れ特性と高域で−40(dB/dec)以上の減衰特性とをもつ高次フィルタで処理し、回転子磁束の非ゼロ回転速度に対応した周波数の基本波成分を主とする被フィルタ処理信号を生成し、被フィルタ処理信号から回転子磁束の位相推定値を生成することを特徴とする。
The invention according to claim 2 is a method for estimating the phase of a rotor magnetic flux of an AC motor having a rotor magnetic flux, wherein the rotation of the rotor magnetic flux is generated using at least the equivalent values of the voltage and current of the stator. -Π / 2 at the frequency corresponding to the non-zero rotation speed of the rotor flux and the speed response band characteristic that changes the passband width according to the rotation of the rotor magnetic flux. The fundamental wave component of the frequency corresponding to the non-zero rotational speed of the rotor magnetic flux is processed using a high-order filter having a phase delay characteristic of (rad) and an attenuation characteristic of −40 (dB / dec) or higher in the high frequency range. And a phase estimation value of the rotor magnetic flux is generated from the filtered signal.

請求項3の発明は、請求項1乃至請求項2のいずれか1項に記載の交流電動機のための位相推定法であって、速応帯域特性をもつ該高次フィルタを、高域で−40(dB/dec)の減衰特性をもつローパスフィルタとすることを特徴とする。The invention of claim 3 is the phase estimation method for an AC motor according to any one of claims 1 to 2, wherein the high-order filter having a quick response band characteristic is The low-pass filter has an attenuation characteristic of 40 (dB / dec).

次に本発明の作用を説明する。請求項1の発明は、回転子が固定子電源周波数と同期して回転する同期電動機を対象に、固定子の電圧と電流との相当値を少なくとも用いて生成した誘起電圧または拡張誘起電圧の相当値を、回転子の回転に応じて通過帯域幅を変化させる速応帯域特性と高域で−40(dB/dec)以上の減衰特性とをもつ高次フィルタで処理し、被フィルタ処理信号から回転子の位相推定値を生成するものである。Next, the operation of the present invention will be described. According to the first aspect of the present invention, for a synchronous motor whose rotor rotates in synchronization with the stator power supply frequency, an equivalent of an induced voltage or an extended induced voltage generated using at least the equivalent values of the stator voltage and current. The value is processed by a high-order filter having a rapid response band characteristic that changes the pass band width according to the rotation of the rotor and an attenuation characteristic of −40 (dB / dec) or higher in the high band, and from the signal to be filtered The phase estimation value of the rotor is generated.

説明の簡略化を図るため、本発明による速応帯域特性をもつ高次フィルタ(以下、速応帯域高次フィルタと略記)として、次の2次ローパスフィルタ考える。

Figure 0005169014
ここに、ζ,ωはフィルタ係数である。In order to simplify the description, the following second-order low-pass filter is considered as a high-order filter (hereinafter, abbreviated as a fast-response band high-order filter) having a rapid response band characteristic according to the present invention.
Figure 0005169014
Here, ζ and ω f are filter coefficients.

図1は、(10)式の周波数特性(振幅特性と位相特性)を、ζ=2、ω=100π(rad/s)を条件に示したものである。同図より容易に理解されるように、(10)式の2次ローパスフィルタにおいては、通過帯域幅は概ねωであり、ω以遠の高周波域では−40(dB/dec)の振幅減衰が得られ、また周波数ω=ωにおける位相遅れは正確に−π/2(rad)である。周波数特性における本特徴は、係数ζ,ωの値如何にかかわず不変である。ω=ωにおける位相遅れは純粋積分の位相遅れ−π/2(rad)と同一である点には、特に注意されたい。FIG. 1 shows the frequency characteristics (amplitude characteristics and phase characteristics) of the equation (10) under the conditions of ζ = 2 and ω f = 100π (rad / s). As can be easily understood from the figure, in the second-order low-pass filter of equation (10), the pass bandwidth is approximately ω f , and the amplitude attenuation is −40 (dB / dec) in the high frequency range beyond ω f. And the phase lag at the frequency ω = ω f is exactly −π / 2 (rad). This feature in the frequency characteristics, coefficient zeta, a value how involved not immutable omega f. phase lag in omega = omega f is the point is the same as the phase lag -π / 2 (rad) of the pure integration is to be particularly noted.

回転子が固定子電源周波数と同期して回転する同期電動機において、回転子磁束(等価回転子磁束を含む)以外は既知であるとして、請求項1の発明に従い、固定子の電圧と電流を用いて生成した誘起電圧または拡張誘起電圧の相当値を、速応帯域高次フィルタで処理すると、次の関係を得る。

Figure 0005169014
In the synchronous motor in which the rotor rotates in synchronization with the stator power supply frequency, it is assumed that other than the rotor magnetic flux (including the equivalent rotor magnetic flux) is known, and according to the invention of claim 1, the voltage and current of the stator are used. When the induced voltage or the equivalent value of the extended induced voltage is processed by the fast response band high-order filter, the following relationship is obtained.
Figure 0005169014

(11)式第3式のF(s.ω)sは、ω=ωを中心周波数とする通過帯域をもつバンドパス特性を意味し、中心周波数では次の性質をもつ。

Figure 0005169014
すなわち、中心周波数では、位相偏移(位相進み・位相遅れ)は発生しない。純粋積分処理と同一のローパスフィルタ位相遅れ−π/2(rad)が、回転子磁束に代わって誘起電圧あるいは拡張誘起電圧を処理することにより生じた等価な微分処理に伴う位相進み+π/2(rad)により、相殺された形となっている。(11) Equation F (s.ω f ) s in the third equation means a band pass characteristic having a pass band with ω = ω f as a center frequency, and has the following property at the center frequency.
Figure 0005169014
That is, no phase shift (phase advance / phase lag) occurs at the center frequency. The same low-pass filter phase lag -π / 2 (rad) as in the pure integration process is the phase lead + π / 2 (with the equivalent differential process generated by processing the induced voltage or the extended induced voltage instead of the rotor magnetic flux. rad) is offset.

図2は、図1に対応したF(s.ω)sの周波数特性を示したものである。回転子磁束に対しては、(12)式の特性に加え、ω以遠の高周波域では、−20(dB/dec)の振幅減衰が維持されていることが確認される。一方で、固定子の電圧・電流を用いて生成された誘起電圧相当値あるいは拡張誘起電圧相当値に含まれ、かつ高周波域に属する各種誤差成分に対しては、2次のローパスフィルタの効果があり、−40(dB/dec)の減衰を受けることになる(図1参照)。FIG. 2 shows the frequency characteristic of F (s.ω f ) s corresponding to FIG. For rotor flux, (12) was added to the characteristics of the type, in the high frequency range of omega f beyond, that amplitude attenuation of -20 (dB / dec) is maintained is confirmed. On the other hand, the effect of the second-order low-pass filter is effective against various error components that are included in the induced voltage equivalent value or the extended induced voltage equivalent value generated using the stator voltage / current and belong to the high frequency range. Yes, it will be attenuated by -40 (dB / dec) (see FIG. 1).

回転子磁束の基本波成分は、回転子の電気速度と同一であることを考慮するならば、フィルタ係数ωを電気速度絶対値と等しく選定する場合には(すなわち、ω=|ω2n|と選定する場合には)、図2より、回転子磁束基本波成分を、位相誤差なく検出できることがわかる。また、非正弦着磁に起因する高調波成分を含む回転子磁束を推定する場合には、図2の振幅特性は、回転子磁束推定値に含まれる高調波成分は減衰され、回転子磁束推定値はより正弦波的になることを意味する。具体的には、回転子磁束のn次高調波成分は、−20logn(dB)の減衰を受けることになる。例えば、この主成分である第3、第5高調波成分は、各々、−9.5(dB)、−14.0(dB)の減衰を受けることになる。(11)式を用いた以上の説明で既に明らかなように、高次フィルタによる被フィルタ処理信号は回転子磁束推定値(等価回転子磁束推定値を含む)のそのものとなる。If it is considered that the fundamental wave component of the rotor magnetic flux is the same as the electric speed of the rotor, the filter coefficient ω f is selected equal to the absolute value of the electric speed (that is, ω f = | ω 2n FIG. 2 shows that the rotor magnetic flux fundamental wave component can be detected without phase error. Further, when estimating the rotor magnetic flux including the harmonic component caused by non-sinusoidal magnetization, the amplitude characteristic of FIG. 2 is attenuated by the harmonic component included in the estimated value of the rotor magnetic flux. The value means more sinusoidal. Specifically, the nth harmonic component of the rotor magnetic flux is attenuated by −20 logn (dB). For example, the third and fifth harmonic components that are the main components are attenuated by −9.5 (dB) and −14.0 (dB), respectively. As is apparent from the above description using the expression (11), the signal to be filtered by the high-order filter is the rotor magnetic flux estimated value (including the equivalent rotor magnetic flux estimated value) itself.

回転子位相推定値は、高次フィルタによる被フィルタ処理信号である回転子磁束推定値の逆正接処理を通じ、直ちに決定される。このことを考慮するならば、このための回転子磁束推定は、基本波成分のみの位相偏移のない推定とすることが肝要である。ω=|ω2n|と選定する場合には、周波数ω=ωで位相偏移のない特性とω以遠の高周波域で減衰特性とを示す図2の周波数特性は、回転子磁束推定の上からは、特に好ましい特性である。本特性は、請求項1の本発明に従い、誘起電圧相当値あるいは拡張誘起電圧相当値を高次フィルタ(高域で−40(dB/dec)以上の振幅減衰をもつフィルタ)で処理することにより得られたものである点には、特に注意されたい。The rotor phase estimation value is immediately determined through the arc tangent process of the rotor magnetic flux estimation value, which is the filtered signal by the high-order filter. If this is taken into consideration, it is important to estimate the rotor magnetic flux for this purpose without estimating the phase shift of only the fundamental wave component. When ω f = | ω 2n | is selected, the frequency characteristic of FIG. 2 showing the characteristic without phase shift at the frequency ω = ω f and the attenuation characteristic in the high frequency range beyond ω f is the rotor magnetic flux estimation. From above, it is a particularly preferable characteristic. According to the present invention of claim 1, this characteristic is obtained by processing an induced voltage equivalent value or an extended induced voltage equivalent value with a high-order filter (a filter having an amplitude attenuation of −40 (dB / dec) or higher in a high band). Special attention should be paid to the fact that it was obtained.

上の例におけるω=|ω2n|の選定は、すなわち回転子電気速度ω2nに応じたローパスフィルタ係数ωの変更は、フィルタ係数ωとフィルタ通過帯域幅とが概ね同一であることを考慮すると、速度に応じてローパスフィルタの通過帯域幅が変化することを意味し、速応帯域性が確立されることを意味する。このように、請求項1の本発明によれば、高次フィルタは、回転子の回転に応じて通過帯域幅を変化させる速応帯域特性もつことになる。ひいては、位相推定値に外乱として作用しる電圧、電流、パラメータ等に含まれる種々の誤差を、回転速度に応じた速度で低減・抑圧する作用が得られる。The selection of ω f = | ω 2n | in the above example, that is, the change of the low-pass filter coefficient ω f according to the rotor electrical speed ω 2n is that the filter coefficient ω f and the filter pass bandwidth are substantially the same. Is taken into account, it means that the pass bandwidth of the low-pass filter changes according to the speed, which means that the quick response bandwidth is established. Thus, according to the first aspect of the present invention, the high-order filter has a quick response band characteristic that changes the pass band width according to the rotation of the rotor. As a result, the effect | action which reduces and suppresses various errors contained in the voltage, electric current, parameter, etc. which act as a disturbance to a phase estimated value at the speed according to a rotational speed is acquired.

以上の説明より明白なように、請求項1の発明によれば、交流電動機の中で、特に回転子が固定子電源周波数と同期して回転する同期電動機の回転子の位相推定において、回転子磁束(等価回転子磁束を含む)の基本波成分を、位相誤差のない状態で、さらには、電圧、電流、パラメータ等に含まれる種々の誤差を回転速度に応じ低減・抑圧した状態で、推定できると言う作用が得られる。As apparent from the above description, according to the invention of claim 1, in the phase estimation of the rotor of the synchronous motor in which the rotor rotates in synchronization with the stator power supply frequency, in the AC motor, the rotor Estimates the fundamental wave component of magnetic flux (including equivalent rotor magnetic flux) with no phase error, and with various errors in voltage, current, parameters, etc. reduced / suppressed according to the rotational speed The effect that it can be obtained.

続いて、請求項2の作用を説明する。請求項2の発明は、永久磁石同期電動機のように回転子磁束を有する同期電動機、誘導電動機を含む交流電動機を対象に、固定子の電圧と電流との相当値を用いて生成した誘起電圧または拡張誘起電圧の相当値を、回転子磁束の回転に応じて通過帯域幅を変化させる速応帯域特性と高域で−40(dB/dec)以上の減衰特性とをもつ高次フィルタで処理し、被フィルタ処理信号から回転子磁束の位相推定値を生成するものである。Next, the operation of claim 2 will be described. The invention of claim 2 is directed to a synchronous motor having a rotor magnetic flux, such as a permanent magnet synchronous motor, or an AC motor including an induction motor, and an induced voltage generated using an equivalent value of a voltage and a current of a stator. The equivalent value of the expansion induced voltage is processed by a high-order filter having a quick response band characteristic that changes the pass band width according to the rotation of the rotor magnetic flux and an attenuation characteristic of −40 (dB / dec) or higher in the high band. The phase estimation value of the rotor magnetic flux is generated from the filtered signal.

請求項2の発明の請求項1の発明対する唯一の違いは、高次フィルタが有する速応帯域特性を、回転子速度に代わって、回転子磁束速度に応じて持たせる点にある。簡明のため、(11)式の2次フィルタを用いて説明するならば、フィルタの帯域幅を定める可変フィルタ係数ωを、「回転子速度」に代わって、「回転子磁束の速度」に応じて変化させる点にある。他に関しては、請求項2の発明と請求項1の発明との違いはない。従って、請求項2の発明の作用の説明は、請求項1の発明に関する(10)式以降の説明において、「回転子速度」を「回転子磁束速度」と形式的に読み代えるだけで、そのまま成立する。当然のことながら、本読み替えにおいては、回転子速度を示したω2nは回転子磁束速度と読み替える必要がある。以上の読み替えより明白なように、請求項2の発明による作用は、適用可能な電動機の相違を除けば、請求項1の発明による作用と、基本的に同一である。請求項2の発明による場合にも、高次フィルタによる被フィルタ処理信号は回転子磁束推定値のそのものとなる。The only difference between the invention of claim 2 and the invention of claim 1 resides in that the high-speed filter has a quick response band characteristic according to the rotor magnetic flux speed instead of the rotor speed. For the sake of simplicity, if a description is given using the secondary filter of the equation (11), the variable filter coefficient ω f that defines the filter bandwidth is changed to “rotor magnetic flux speed” instead of “rotor speed”. The point is to change it accordingly. With respect to other aspects, there is no difference between the invention of claim 2 and the invention of claim 1. Therefore, the explanation of the operation of the invention of claim 2 can be made by simply replacing “rotor speed” with “rotor magnetic flux speed” formally in the explanation after the expression (10) relating to the invention of claim 1. To establish. As a matter of course, in this replacement, ω 2n indicating the rotor speed needs to be replaced with the rotor magnetic flux speed. As is clear from the above replacement, the operation of the invention of claim 2 is basically the same as the operation of the invention of claim 1 except for the difference in applicable electric motors. Also in the case of the invention of claim 2, the signal to be filtered by the high-order filter is the rotor magnetic flux estimation value itself.

上記の読み替えによる作用の妥当性は、交流電動機として、永久磁石同期電動機のように回転子磁束を有する同期電動機を考えれば、一層明瞭である。当業者には周知のように、回転子磁束を有する同期電動機においては、回転子位相と回転子磁束位相は同一であり、また、回転子速度と回転子磁束速度は同一である。従って、回転子磁束を有する同期電動機を対象とする場合には、請求項2の発明と請求項1の発明の違いは、無い。当然のことながら、両発明は同一の作用を示す。The validity of the action by the above replacement is clearer when a synchronous motor having a rotor magnetic flux such as a permanent magnet synchronous motor is considered as an AC motor. As is well known to those skilled in the art, in a synchronous motor having a rotor magnetic flux, the rotor phase and the rotor magnetic flux phase are the same, and the rotor speed and the rotor magnetic flux speed are the same. Therefore, when the synchronous motor having the rotor magnetic flux is targeted, there is no difference between the invention of claim 2 and the invention of claim 1. Naturally, both inventions exhibit the same effect.

以上より明らかなように、請求項2の発明によれば、回転子磁束を有する交流電動機の回転子磁束の位相推定において、回転子磁束の基本波成分を、位相誤差のない状態で、さらには、電圧、電流、パラメータ等に含まれる種々の誤差を回転速度に応じ低減・抑圧した状態で、推定できると言う作用が得られる。As apparent from the above, according to the invention of claim 2, in the phase estimation of the rotor magnetic flux of the AC motor having the rotor magnetic flux, the fundamental wave component of the rotor magnetic flux is further reduced in a state without phase error. Thus, an effect is obtained that estimation can be performed in a state where various errors included in the voltage, current, parameter, and the like are reduced / suppressed according to the rotation speed.

続いて、請求項3の発明の作用を説明する。請求項1の作用の説明、特に(10)式を用いた説明で明白なように、回転子磁束(等価回転子磁束を含む)に含まれる高調波成分を除去し、基本波成分を抽出・推定するには、最低限でも、速応帯域特性をもつ高次フィルタに高域で−40(dB/dec)の減衰特性をもたせる必要がある。換言するならば、高次フィルタに高域で−40(dB/dec)の減衰特性をもたせれば、回転子磁束(等価回転子磁束を含む)に含まれる高調波成分を除去し、基本波成分を抽出・推定することが可能となる。請求項3の発明によれば、速応帯域特性をもつ高次フィルタを、高域で−40(dB/dec)の減衰特性をもつローパスフィルタとするので、請求項1または請求項2の作用を、最小の処理負荷で得ることができると言う作用が得られる。Next, the operation of the invention of claim 3 will be described. As is apparent from the description of the operation of claim 1, particularly the description using the equation (10), the harmonic component contained in the rotor magnetic flux (including the equivalent rotor magnetic flux) is removed, and the fundamental wave component is extracted. In order to estimate, it is necessary to provide a high-order filter having a quick response band characteristic with an attenuation characteristic of −40 (dB / dec) in the high band at least. In other words, if the high-order filter has an attenuation characteristic of −40 (dB / dec) in the high band, the harmonic component contained in the rotor magnetic flux (including the equivalent rotor magnetic flux) is removed, and the fundamental wave It becomes possible to extract and estimate the components. According to the invention of claim 3, since the high-order filter having the quick response band characteristic is a low-pass filter having an attenuation characteristic of −40 (dB / dec) in the high band, the effect of claim 1 or claim 2 is achieved. Can be obtained with a minimum processing load.

以下、図面を用いて、本発明の実施形態を詳細に説明する。図3に、永久磁石同期電動機を対象として、本発明の位相推定法を利用したセンサレスベクトル制御システム(トルク制御モード)の1実施形態例を示した。図3においては、位相決定器2以外の機器は、通常のセンサ利用ベクトル制御システムの構成と同一である。これら機器の働きは、当業者には容易に理解されるので、説明を省略する。参考までに、同図面における機器の意味を、本明細書末尾の「符号の説明」につけたので、参照されたい。本発明の核心は、位相決定器2にある。位相決定器の目的は、固定子の電圧・電流相当値を利用して、固定αβ座標

Figure 0005169014
測定値を、固定子電圧相当値としては電圧指令値を利用する例となっている。Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. FIG. 3 shows an embodiment of a sensorless vector control system (torque control mode) using the phase estimation method of the present invention for a permanent magnet synchronous motor. In FIG. 3, the devices other than the phase determiner 2 have the same configuration as that of a normal sensor utilization vector control system. Since the functions of these devices are easily understood by those skilled in the art, description thereof is omitted. For reference, the meaning of the device in the drawing is given in “Explanation of Symbols” at the end of this specification, so please refer to it. The core of the present invention is the phase determiner 2. The purpose of the phase determiner is to use the fixed voltage and current equivalent values of the stator to determine the fixed αβ coordinates.
Figure 0005169014
In this example, the voltage command value is used as the measured value corresponding to the stator voltage equivalent value.

図4は、図3における位相決定器の内部構成の1例を示したものである。位相決定器2は、本発明による速応帯域高次フィルタ2a、速度推定器2b、係数生成器2cから構成されている。FIG. 4 shows an example of the internal configuration of the phase determiner in FIG. The phase determiner 2 includes a fast response band high-order filter 2a, a speed estimator 2b, and a coefficient generator 2c according to the present invention.

速応帯域高次フィルタ2aは、被フィルタ処理信号として固定αβ座標系上の回転子磁束推定値を出力し、この逆正接演算を通じて、回転子位相推定値(初期回転子位相推定値

Figure 0005169014
子電圧・電流の相当値としての固定子電流実測値と固定子電圧指令値に加えて、最終回転
Figure 0005169014
請求項1乃至請求項3のすべての発明に合致した、更には(10)〜(12)式の説明に利用した2次フィルタを採用した例を示している。本例では、固定子の電圧・電流相当値に対する処理を、次の(13)式のように行なっている。
Figure 0005169014
すなわち、固定子反作用磁束微分値の入力に代わって、固定子反作用磁束そのものを速応帯域高次フィルタに入力するようにしている。当然のことながら、この場合の固定子反作用磁束の入力すべき端子は、同図に明示したように、電圧相当値を入力する端子と異なることになる。The fast response band high-order filter 2a outputs a rotor magnetic flux estimated value on a fixed αβ coordinate system as a filtered signal, and through this arctangent calculation, a rotor phase estimated value (initial rotor phase estimated value).
Figure 0005169014
In addition to the actual stator current value and the stator voltage command value as the equivalent values of the child voltage and current, the final rotation
Figure 0005169014
An example is shown in which a second-order filter that matches all the inventions of claims 1 to 3 and that is used in the explanation of the expressions (10) to (12) is employed. In this example, processing for the voltage / current equivalent value of the stator is performed as shown in the following equation (13).
Figure 0005169014
That is, instead of inputting the stator reaction magnetic flux differential value, the stator reaction magnetic flux itself is input to the fast response band high-order filter. As a matter of course, the terminal to which the stator reaction magnetic flux is to be input in this case is different from the terminal to which the voltage equivalent value is input as clearly shown in FIG.

対象とする永久磁石同期電動機が突極特性をもつ場合には、固定子反作用磁束φは、最

Figure 0005169014
Figure 0005169014
ここに、Iは2x2単位行列であり、Q(・)は次式で定義された鏡行列である。
Figure 0005169014
図5では、最終位相推定値を利用した固定子反作用磁束生成の様子(ブロック2a−1)を破線の貫徹矢印で表現している。非突極の場合には、固定子反作用磁束φは、最終位相推定値の要なく、次の(16)式のように生成できる。
Figure 0005169014
When the target permanent magnet synchronous motor has salient pole characteristics, the stator reaction magnetic flux φ i is
Figure 0005169014
Figure 0005169014
Here, I is a 2 × 2 unit matrix, and Q (•) is a mirror matrix defined by the following equation.
Figure 0005169014
In FIG. 5, the state of the stator reaction magnetic flux generation using the final phase estimated value (block 2a-1) is expressed by a penetrating broken arrow. In the case of a non-salient pole, the stator reaction magnetic flux φ i can be generated as in the following equation (16) without the need for the final phase estimation value.
Figure 0005169014

速度推定器2bは、積分フィードバック形速度推定法に従い構成している。積分フィードバック形速度推定法の詳細は、下記文献に詳しく説明されているので、この説明は省略する。
(1)新中新二:「永久磁石同期モータの最小次元D因子状態オブザーバとこれを用いたセンサレスベクトル制御法の提案」、電気学会論文誌D分冊、Vol.123、No.12、pp.1446−1460(2003−12)
速度推定器2bは、初期回転子位相推定値を入力として受け、最終回転子位相推定値と同電気速度推定値を出力している(図4参照)。回転子電気速度推定値は、係数生成器2cへフィードバック利用されている。速度制御モードの場合には、位相決定器2の外部にも電気速度推定値を出力することになる。図4では、この様子を破線で示している。
The speed estimator 2b is configured according to an integral feedback type speed estimation method. The details of the integral feedback type speed estimation method are described in detail in the following document, and the description thereof will be omitted.
(1) Shinnaka Shinji: “Proposal of minimum dimension D-factor state observer of permanent magnet synchronous motor and sensorless vector control method using it”, IEEJ Transactions D Volume, Vol. 123, no. 12, pp. 1446-1460 (2003-12)
The speed estimator 2b receives the initial rotor phase estimated value as an input, and outputs the final rotor phase estimated value and the same electric speed estimated value (see FIG. 4). The estimated value of the rotor electrical speed is used as feedback to the coefficient generator 2c. In the case of the speed control mode, the estimated electric speed value is also output to the outside of the phase determiner 2. In FIG. 4, this state is indicated by a broken line.

係数生成器2cは、速度推定器から得た電気速度推定値を処理し、速応帯域高次フィルタの可変係数ωを、速応的に生成している。本発明による速応帯域高次フィルタで所期の作用を得るには、フィルタ係数ωに関しては、実質的にω=|ω2n|なる関係を達成する必要がある(なお、速応帯域の2次フィルタの場合には、係数ζは一定で問題ない)。係数

Figure 0005169014
リミッタ処理の3処理を施して得た最終処理信号をフィルタ係数ωとし、実質的に本関係を達成している。The coefficient generator 2c processes the electrical speed estimation value obtained from the speed estimator, and quickly generates the variable coefficient ω f of the fast response band high-order filter. In order to obtain the desired effect with the rapid response band high-order filter according to the present invention, it is necessary to substantially satisfy the relationship of ω f = | ω 2n | with respect to the filter coefficient ω f (note that the rapid response band) In the case of the second-order filter, the coefficient ζ is constant and has no problem). coefficient
Figure 0005169014
The final processed signal obtained by applying the 3 treatment of the limiter processing and the filter coefficient omega f, have achieved substantially the relationship.

Figure 0005169014
り、簡単には、次の一定係数をもつローパスフィルタによる処理でよい。
Figure 0005169014
このときのフィルタは次の1次のものでよい。
Figure 0005169014
脈動除去処理は、ローパスフィルタ処理に代わって、次の機能をもつマルチレイトサンプルホールダ処理でも達成できる。
Figure 0005169014
絶対値処理は、フィルタ係数ωが正であることを保証するためのものである。下限リミッタ処理は、フィルタ係数ωの下限値をある正の微小値に設定し、計算誤差が存在する場合にも、速応帯域高次フィルタの不安定化を防止するためのものである。
Figure 0005169014
In short, processing by a low-pass filter having the following constant coefficient may be used.
Figure 0005169014
The filter at this time may be the following primary filter.
Figure 0005169014
The pulsation removing process can be achieved by a multi-rate sample holder process having the following functions instead of the low-pass filter process.
Figure 0005169014
The absolute value processing is for assuring that the filter coefficient ω f is positive. Lower limiter processing is set to a positive small value with the lower limit value of the filter coefficients omega f, when the calculation error is present is also intended to prevent the destabilization of the speed応帯range order filter.

Figure 0005169014
からのフィルタ係数ωの生成例として、(17)式のローパスフィルタを用いた例を示している。また、同図では、絶対値処理と下限リミッタ処理とを、1つのブロックで表現している。
Figure 0005169014
As an example of generating the filter coefficient ω f from, an example using the low-pass filter of the equation (17) is shown. In the figure, the absolute value process and the lower limiter process are expressed by one block.

上記の実施形態例では、説明の簡明性と具体性とを確保するために、速応帯域高次フィルタは2次のローパスフィルタとした。一般には、速応帯域高次フィルタとして、n次ローパスフィルタを用いることが可能であることを指摘しておく。例えば、速応帯域高次フィルタを3次とする場合には、つぎのものが有用である。

Figure 0005169014
In the above embodiment, the fast response band high-order filter is a second-order low-pass filter in order to ensure the simplicity and concreteness of the description. In general, it should be pointed out that an n-order low-pass filter can be used as the fast response band high-order filter. For example, when the fast response band high-order filter is third order, the following is useful.
Figure 0005169014

上の実施形態例では、速応帯域高次フィルタを固定αβ座標系上で構成するものとした(図3参照)。本座標系上での速応帯域高次フィルタは、1入力1出力フィルタを2個並列に配したものである((11)式参照)。これに代わって、速応帯域高次フィルタは、回転子位相に位相差のない同期を目指した準同期γδ座標系上でも構成することができる。図7は、この1実施形態例である。位相決定器2が、ベクトル回転器6a、6bより左側に構成されている点が、図3の実施形態例と異なっている。In the above embodiment, the fast response band high-order filter is configured on the fixed αβ coordinate system (see FIG. 3). The fast response band high-order filter on this coordinate system is one in which two one-input one-output filters are arranged in parallel (see equation (11)). Alternatively, the fast response band high-order filter can also be configured on a quasi-synchronous γδ coordinate system aimed at synchronization without a phase difference in the rotor phase. FIG. 7 shows an example of this embodiment. The phase determiner 2 is different from the embodiment in FIG. 3 in that the phase determiner 2 is arranged on the left side of the vector rotators 6a and 6b.

図8は、図7における位相決定器の内部構成の1例を示したものである。位相決定器2は、本発明による速応帯域高次フィルタ2a、位相同期器2d、係数生成器2cから構成されている。FIG. 8 shows an example of the internal configuration of the phase determiner in FIG. The phase determiner 2 includes a fast response band high-order filter 2a according to the present invention, a phase synchronizer 2d, and a coefficient generator 2c.

準同期γδ座標系用の速応帯域高次フィルタは、固定αβ座標系用の速応帯域高次フィルタのラプラス演算子を形式的にD因子に置換することにより得ることができる。この結果、準同期γδ座標系上の速応帯域高次フィルタは、単一の2入力2出力D因子フィルタとなる。図5の固定αβ座標系用の速応帯域高次フィルタに対応した、準同期γδ座標系用の速応帯域高次フィルタを、図9に示した。同図のD因子は、以下のように定義されている2x2行列である。

Figure 0005169014
The fast response band high-order filter for the quasi-synchronous γδ coordinate system can be obtained by formally replacing the Laplace operator of the fast response band high-order filter for the fixed αβ coordinate system with the D factor. As a result, the fast response band high-order filter on the quasi-synchronous γδ coordinate system becomes a single 2-input 2-output D-factor filter. FIG. 9 shows a fast response band high-order filter for the quasi-synchronous γδ coordinate system corresponding to the fast response band high-order filter for the fixed αβ coordinate system of FIG. The D factor in the figure is a 2 × 2 matrix defined as follows.
Figure 0005169014

図9では、D因子は逆行列として使用されているが、この実現法は、上記文献(1)に既に詳しく紹介されているので、この説明は省略する。なお、図9においては、速応帯域高次フィルタへ入力される固定子反作用磁束は、座標系の準同期性を考慮し、以下のように生成されている。

Figure 0005169014
当然のことながら、非突極の場合には、L=L=Lであり、固定子反作用磁束は、固定子電流相当値にスカラーのインダクタンスを乗ずることにより、直ちに生成される((16)式参照)。準同期γδ座標系用の速応帯域高次フィルタは、準同期γδ座標系上での回転子磁束を推定しており、最終的には、被フィルタ処理信号である回転子磁束推定値を逆正接処理し、同位相を生成・出力している。本位相は、準同期γδ座標系上からみた回転子位相推定値となっている。In FIG. 9, the factor D is used as an inverse matrix, but since this realization method has already been introduced in detail in the above document (1), this description is omitted. In FIG. 9, the stator reaction magnetic flux input to the fast response band high-order filter is generated as follows in consideration of the quasi-synchronization of the coordinate system.
Figure 0005169014
Of course, in the case of a non-salient pole, L d = L q = L i and the stator reaction magnetic flux is immediately generated by multiplying the value corresponding to the stator current by the scalar inductance (( (See 16)). The fast response band high-order filter for the quasi-synchronous γδ coordinate system estimates the rotor magnetic flux on the quasi-synchronous γδ coordinate system, and finally reverses the rotor flux estimation value that is the filtered signal. Tangent processing is performed to generate and output the same phase. This phase is a rotor phase estimation value viewed from the quasi-synchronous γδ coordinate system.

位相同期器2dは、一般化積分形PLL法に従い構成している。位相同期器2dは、準同期γδ座標系上からみた回転子位相推定値を入力として、固定αβ座標系上の回転子の位相推定値と同電気速度推定値を出力している。一般化積分形PLL法の詳細は、上述の文献(1)に詳しく説明されているので、この説明は省略する。The phase synchronizer 2d is configured according to a generalized integral PLL method. The phase synchronizer 2d receives the rotor phase estimation value viewed from the quasi-synchronous γδ coordinate system, and outputs the rotor phase estimation value and the electric speed estimation value on the fixed αβ coordinate system. The details of the generalized integral PLL method are described in detail in the above-mentioned document (1), and thus the description thereof is omitted.

図3〜9を用いて説明した実施形態例は、永久磁石同期電動機を対象した場合の位相推定法の実施形態例であるが、同期リラクタンス電動機を対象とする場合にも、固定子反作用磁束の生成に関し若干の変更を行なうだけで、これを実質そのまま利用できる。具体的には、固定αβ座標系で速応帯域高次フィルタを構成する場合には、これへの入力である拡張誘起電圧の相当値を、固定子の電圧・電流相当値を用いて以下の関係に基づき生成すればよい。

Figure 0005169014
また、準同期γδ座標系で速応帯域高次フィルタを構成する場合には、これへの入力である拡張誘起電圧の相当値を、固定子の電圧・電流相当値を用いて以下の関係に基づき生成すればよい。
Figure 0005169014
The embodiment described with reference to FIGS. 3 to 9 is an embodiment of the phase estimation method for a permanent magnet synchronous motor. However, even in the case of a synchronous reluctance motor, the stator reaction magnetic flux This can be used as is, with only minor changes to the generation. Specifically, when a fast response band high-order filter is configured in a fixed αβ coordinate system, an equivalent value of an extended induced voltage that is an input to this is expressed as follows using a voltage / current equivalent value of a stator. What is necessary is just to produce | generate based on a relationship.
Figure 0005169014
Also, when configuring a fast response band high-order filter in a quasi-synchronous γδ coordinate system, the equivalent value of the extended induced voltage that is input to this is expressed as follows using the voltage / current equivalent value of the stator. It may be generated based on this.
Figure 0005169014

以上説明したように、図3〜9を用いて説明した実施形態例は、速応帯域高次フィルタへの入力を(23)式または(24)式へと変更するだけ、同期リラクタンス電動機の場合の実施形態例となる。なお、同期リラクタンス同期の場合には、(13)式左辺に従い、速応帯域高次フィルタへの入力信号を2分することなく印加した方が、都合がよい。As described above, the embodiment described with reference to FIGS. 3 to 9 is the case of the synchronous reluctance motor only by changing the input to the fast response band high-order filter to the equation (23) or (24). This is an example embodiment. In the case of synchronous reluctance synchronization, it is more convenient to apply the input signal to the fast response band high-order filter without dividing it into two according to the left side of equation (13).

続いて、誘導電動機を対象とする場合の実施形態例を説明する。誘導電動機を対象とした場合の回転子磁束推定法も、図3〜9を用いて説明した永久磁石同期電動機を対象した場合の位相推定法の実施形態例に対し、固定子反作用磁束の生成に関し若干の変更を行なうだけで、これを実質そのまま利用できる。具体的には、速応帯域高次フィルタの構成を示した図5、図9において、固定子反作用磁束生成ブロック2a−1に次の関係を適用し、固定子反作用磁束を生成するようにすればよい((9)式参照)。

Figure 0005169014
Subsequently, an embodiment example in the case of targeting an induction motor will be described. The rotor magnetic flux estimation method for the induction motor is also related to the generation of the stator reaction magnetic flux in comparison with the embodiment example of the phase estimation method for the permanent magnet synchronous motor described with reference to FIGS. You can use it as it is with just a few changes. Specifically, in FIGS. 5 and 9 showing the configuration of the fast response band high-order filter, the following relation is applied to the stator reaction magnetic flux generation block 2a-1 to generate the stator reaction magnetic flux. (See equation (9)).
Figure 0005169014

(25)式と(16)式との比較より明らかなように、誘導電動機を対象とする場合の実施形態例は、非突極の永久磁石同期電動機を対象とした場合の実施形態例と形式的には完全同一となる。ただし、図3〜9を用いて説明した永久磁石同期電動機を対象した場合の位相推

Figure 0005169014
回転子磁束速度推定値と読み替える必要がある。繰返し説明しているように、永久磁石同期電動においては、回転子速度も回転子磁束速度も同一である。誘導電動機においては、この同一性は一般には成立しない。回転子磁束の推定に必要な速度情報は、回転子磁束の速度である。As is clear from the comparison between the equations (25) and (16), the embodiment in the case of targeting an induction motor is the embodiment and the form in the case of a non-salient permanent magnet synchronous motor. It becomes completely the same. However, the phase estimation for the permanent magnet synchronous motor described with reference to FIGS.
Figure 0005169014
It should be read as the rotor flux velocity estimate. As described repeatedly, in the permanent magnet synchronous motor, the rotor speed and the rotor magnetic flux speed are the same. In an induction motor, this identity generally does not hold. The speed information necessary for estimating the rotor magnetic flux is the speed of the rotor magnetic flux.

以上、永久磁石同期電動機、同期リラクタンス電動機、誘導電動機を対象として、本発明による回転子位相推定法及び回転子磁束位相推定法を、固定αβ座標系上で構成する場合、準同期γδ座標系上で構成する場合について、詳しく説明した。本発明は、固定子の電気的関係を利用してこれら位相を推定するものであり、同一の固定子をもつ他の交流電動機にも適用可能であることを指摘しておく。同一の固定子をもつ他の交流電動機としては、例えば、巻線形同期電動機、回転子に巻線と永久磁石をもつハイブリッド界磁形同期電動機、回転子に同期用永久磁石と誘導用カゴをもつハイブリッド電動機などがある。As described above, when the rotor phase estimation method and the rotor magnetic flux phase estimation method according to the present invention are configured on a fixed αβ coordinate system for permanent magnet synchronous motors, synchronous reluctance motors, and induction motors, The case where it is configured with is described in detail. It should be pointed out that the present invention estimates these phases using the electrical relationship of the stator and is applicable to other AC motors having the same stator. Other AC motors having the same stator include, for example, a wound synchronous motor, a hybrid field synchronous motor having a winding and a permanent magnet in the rotor, and a synchronous permanent magnet and an induction cage in the rotor There are hybrid electric motors.

図3〜図9を用いて説明した実施形態例では、電圧相当値としては電圧指令値を利用したが、電圧実測値や他の近似値を利用してもよい。また電流相当値としては電流実測値を利用したが、電流指令値や他の近似値を利用してもよい。In the embodiment described with reference to FIGS. 3 to 9, the voltage command value is used as the voltage equivalent value, but an actual voltage measurement value or other approximate value may be used. Moreover, although the current measured value is used as the current equivalent value, a current command value or another approximate value may be used.

図3及び図7を用いた実施形態例では、回転子位相あるいは回転子磁束位相の推定が駆動制御の鍵となるトルク制御モードを対象に、本発明を説明した。当業者には容易に理解されるように、本発明は、トルク制御モードに限定されるものではなく、速度制御モードなど他のモードによる駆動にも適用可能である。In the embodiment using FIG. 3 and FIG. 7, the present invention has been described for the torque control mode in which the estimation of the rotor phase or the rotor magnetic flux phase is the key to the drive control. As will be readily understood by those skilled in the art, the present invention is not limited to the torque control mode, but can be applied to driving in other modes such as a speed control mode.

図3及び図7を用いた実施形態例では、電流制御はフィードバック的に行なっている。電流制御をフィードフォワード的に行なう場合にも、あるいは電流制御を実施しない場合にも、本発明は適用可能であることを指摘しておく。In the embodiment using FIG. 3 and FIG. 7, the current control is performed in a feedback manner. It should be pointed out that the present invention is applicable even when the current control is performed in a feedforward manner or when the current control is not performed.

本発明による位相決定器は、アナログ的に実現可能であるが、最近のディジタル技術の著しい進歩を考えるとディジタル的に構成することが好ましい。ディジタル構成はハードウェア的構成とソフトウェア的構成があるが、当業者にとっては既に自明のように本発明はいずれでも構成できる。以上、本発明に関し、各種の図を利用しつつ多様な実施形態例を詳しく説明した。Although the phase determiner according to the present invention can be implemented in an analog manner, it is preferably constructed digitally in view of the remarkable progress of recent digital technology. Although the digital configuration includes a hardware configuration and a software configuration, as will be apparent to those skilled in the art, any of the present invention can be configured. The present invention has been described in detail with reference to various embodiments using various drawings.

発明の効果Effect of the invention

以上の説明より明白なように、本発明は以下の効果を奏する。請求項1の発明によれば、交流電動機の中で、特に回転子が固定子電源周波数と同期して回転する同期電動機の回転子の位相推定において、回転子磁束(等価回転子磁束を含む)の基本波成分を、位相誤差のない状態で、さらには、電圧、電流、パラメータ等に含まれる種々の誤差を回転速度に応じ低減・抑圧した状態で、推定できると言う作用が得られた。ひいては、本作用の結果、請求項1の発明によれば、同期電動機の回転子磁束推定値(等価回転子磁束推定値を含む)から、回転子位相を、位相の狂いが出ない状態でかつ外乱影響を抑圧した状態で推定できると言う効果が得られる。As is clear from the above description, the present invention has the following effects. According to the first aspect of the present invention, in the estimation of the phase of the rotor of the synchronous motor in which the rotor rotates in synchronization with the stator power supply frequency among the AC motors, the rotor magnetic flux (including the equivalent rotor magnetic flux) is used. The fundamental wave component can be estimated in a state where there is no phase error, and further, in a state where various errors included in the voltage, current, parameter, and the like are reduced / suppressed according to the rotational speed. As a result, as a result of this action, according to the first aspect of the present invention, the rotor phase is calculated from the rotor magnetic flux estimated value (including the equivalent rotor magnetic flux estimated value) of the synchronous motor in a state where no phase deviation occurs. The effect that it can estimate in the state which suppressed disturbance influence is acquired.

次に、請求項2の本発明による効果を説明する。請求項2の発明によれば、回転子磁束を有する交流電動機の回転子磁束の位相推定において、回転子磁束の基本波成分を、位相誤差のない状態で、さらには、電圧、電流、パラメータ等に含まれる種々の誤差を回転速度に応じ低減・抑圧した状態で、推定できると言う作用が得られた。ひいては、本作用の結果、請求項2の発明によれば、回転子磁束を有する交流電動機の回転子磁束推定値から、回転子磁束の位相を、位相の狂いが出ない状態でかつ外乱影響を抑圧した状態で推定できると言う効果が得られる。Next, the effect of the present invention of claim 2 will be described. According to the invention of claim 2, in the phase estimation of the rotor magnetic flux of the AC motor having the rotor magnetic flux, the fundamental wave component of the rotor magnetic flux is in a state without phase error, and further, the voltage, current, parameter, etc. The effect that it is possible to estimate in a state where various errors included in is reduced / suppressed according to the rotation speed was obtained. As a result, as a result of this action, according to the invention of claim 2, from the estimated value of the rotor magnetic flux of the AC motor having the rotor magnetic flux, the phase of the rotor magnetic flux is affected in a state in which there is no phase deviation and the influence of disturbance. The effect that it can estimate in the suppressed state is acquired.

続いて、請求項3の本発明の効果を説明する。請求項3の発明によれば、請求項1または請求項2の作用を、最小の処理負荷で得ることができると言う作用が得られた。ひいては、本作用の結果、請求項3の本発明によれば、請求項1または請求項2の効果を、最小の処理負荷で得ることができると言う効果が得られる。Next, the effect of the present invention of claim 3 will be described. According to the invention of claim 3, the effect that the action of claim 1 or claim 2 can be obtained with a minimum processing load is obtained. As a result, according to the present invention of claim 3, as a result of this action, the effect of claim 1 or claim 2 can be obtained with a minimum processing load.

速応帯域高次フィルタの1周波数特性例Example of one-frequency characteristics of high-speed fast response band filter 速応帯域高次フィルタに微分処理を加えた場合の1周波数特性例Example of 1-frequency characteristics when differential processing is applied to a high-speed response band high-order filter 1実施形態例における駆動制御装置を含む駆動制御システムのブロック図1 is a block diagram of a drive control system including a drive control device according to an embodiment. 1実施形態例における位相決定器の基本構成を示すブロック図The block diagram which shows the basic composition of the phase determiner in the example of 1 embodiment 1実施形態例における速応帯域高次フィルタの構成例Configuration example of fast response band high-order filter in one embodiment 1実施形態例における係数生成器の構成例Configuration example of coefficient generator in one embodiment 1実施形態例における駆動制御装置を含む駆動制御システムのブロック図1 is a block diagram of a drive control system including a drive control device according to an embodiment. 1実施形態例における位相決定器の基本構成を示すブロック図The block diagram which shows the basic composition of the phase determiner in the example of 1 embodiment 1実施形態例における速応帯域高次フィルタの構成例Configuration example of fast response band high-order filter in one embodiment

符号の説明Explanation of symbols

1 交流電動機
2 位相決定器
2a 速応帯域高次フィルタ
2a−1 固定子反作用磁束生成ブロック
2b 速度推定器
2c 係数生成器
2d 位相同期器
3 電力変換器
4 電流検出器
5a 3相2相変換器
5b 2相3相変換器
6a ベクトル回転器
6b ベクトル回転器
7 余弦正弦信号発生器
8 電流制御器
DESCRIPTION OF SYMBOLS 1 AC motor 2 Phase determiner 2a Fast response band high-order filter 2a-1 Stator reaction magnetic flux generation block 2b Speed estimator 2c Coefficient generator 2d Phase synchronizer 3 Power converter 4 Current detector 5a Three-phase two-phase converter 5b 2-phase 3-phase converter 6a Vector rotator 6b Vector rotator 7 Cosine sine signal generator 8 Current controller

Claims (3)

回転子が固定子電源周波数と同期して回転する交流電動機の回転子の位相推定法であって、
固定子の電圧と電流との相当値を少なくとも用いて生成した、回転子の回転に起因する誘起電圧または拡張誘起電圧の相当値を、回転子の回転に応じて通過帯域幅を変化させる速応帯域特性と回転子の非ゼロ回転速度に対応した周波数で-π/2(rad)の位相遅れ特性と高域で-40(dB/dec)以上の減衰特性とをもつ高次フィルタで処理し、回転子の非ゼロ回転速度に対応した周波数の基本波成分を主とする被フィルタ処理信号を生成し、被フィルタ処理信号から回転子の位相推定値を生成することを特徴とする交流電動機のための回転子の位相推定法。
A phase estimation method for a rotor of an AC motor in which a rotor rotates in synchronization with a stator power supply frequency,
An equivalent value of the induced voltage or extended induced voltage caused by the rotation of the rotor, which is generated using at least the equivalent value of the stator voltage and current, is a speed response that changes the passband width according to the rotation of the rotor. It is processed with a high-order filter that has a phase lag characteristic of -π / 2 (rad) at a frequency corresponding to the non-zero rotation speed of the band characteristics and a rotor, and an attenuation characteristic of -40 (dB / dec) or higher at high frequencies. Generating a filtered signal mainly including a fundamental wave component having a frequency corresponding to a non-zero rotation speed of the rotor, and generating a rotor phase estimation value from the filtered signal. For rotor phase estimation.
回転子磁束を有する交流電動機の回転子磁束の位相推定法であって、
固定子の電圧と電流との相当値を少なくとも用いて生成した、回転子磁束の回転に起因する誘起電圧または拡張誘起電圧の相当値を、回転子磁束の回転に応じて通過帯域幅を変化させる速応帯域特性と回転子磁束の非ゼロ回転速度に対応した周波数で-π/2(rad)の位相遅れ特性と高域で-40(dB/dec)以上の減衰特性とをもつ高次フィルタで処理し、回転子磁束の非ゼロ回転速度に対応した周波数の基本波成分を主とする被フィルタ処理信号を生成し、被フィルタ処理信号から回転子磁束の位相推定値を生成することを特徴とする交流電動機のための回転子磁束の位相推定法。
A phase estimation method for a rotor magnetic flux of an AC motor having a rotor magnetic flux,
The pass band width is changed according to the rotation of the rotor magnetic flux, with the equivalent value of the induced voltage or the expansion induced voltage caused by the rotation of the rotor magnetic flux generated using at least the equivalent value of the stator voltage and current. High-order filter with fast response band characteristics and phase lag characteristics of -π / 2 (rad) at frequencies corresponding to non-zero rotational speed of rotor flux and attenuation characteristics of -40 (dB / dec) or higher at high frequencies To generate a filtered signal mainly including a fundamental wave component having a frequency corresponding to the non-zero rotation speed of the rotor magnetic flux, and to generate a phase estimation value of the rotor magnetic flux from the filtered signal. The rotor magnetic flux phase estimation method for AC motors.
速応帯域特性をもつ該高次フィルタを、高域で-40(dB/dec)の減衰特性をもつローパスフィルタとすることを特徴とする請求項1乃至請求項2のいずれか1項に記載の交流電動機のための位相推定法。 3. The high-order filter having a quick response band characteristic is a low-pass filter having an attenuation characteristic of −40 (dB / dec) in a high band. Phase estimation method for AC motors.
JP2007125010A 2007-04-09 2007-04-09 Phase estimation method for AC motors Active JP5169014B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007125010A JP5169014B2 (en) 2007-04-09 2007-04-09 Phase estimation method for AC motors

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007125010A JP5169014B2 (en) 2007-04-09 2007-04-09 Phase estimation method for AC motors

Publications (2)

Publication Number Publication Date
JP2008263764A JP2008263764A (en) 2008-10-30
JP5169014B2 true JP5169014B2 (en) 2013-03-27

Family

ID=39985838

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007125010A Active JP5169014B2 (en) 2007-04-09 2007-04-09 Phase estimation method for AC motors

Country Status (1)

Country Link
JP (1) JP5169014B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5828441B2 (en) * 2011-09-19 2015-12-09 新中 新二 Rotor phase estimator for synchronous motor
JP6691330B2 (en) * 2016-02-14 2020-04-28 有限会社シー・アンド・エス国際研究所 Phase velocity estimation device using multiphase low resolution detector signals

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3612636B2 (en) * 1996-09-18 2005-01-19 有限会社シー・アンド・エス国際研究所 Vector control method for synchronous motor
JP2000312499A (en) * 1999-04-27 2000-11-07 Meidensha Corp Vector controller for induction motor
JP2004080975A (en) * 2002-08-22 2004-03-11 Nissan Motor Co Ltd Motor control device

Also Published As

Publication number Publication date
JP2008263764A (en) 2008-10-30

Similar Documents

Publication Publication Date Title
CN1322667C (en) Sensorless controller of AC motor and control method
Lin et al. Improved rotor position estimation in sensorless-controlled permanent-magnet synchronous machines having asymmetric-EMF with harmonic compensation
CN102844979B (en) The rotor phase speed estimating device of alternating current motor
CN101682283B (en) Sensorless control device for permanent magnet synchronous motor
CN113866480B (en) Detection of offset errors in phase current measurements for motor control systems
JP4989075B2 (en) Electric motor drive control device and electric motor drive system
JP7218700B2 (en) motor controller
Malekipour et al. A closed-loop PMSM sensorless control based on the machine acoustic noise
CN114204865A (en) Method and device for estimating electrical angular velocity of permanent magnet synchronous generator
WO2020105204A1 (en) Power conversion device
JP4899509B2 (en) AC motor rotor phase estimation device
JP2015180130A (en) Control device of permanent magnet synchronous motor
Farahat et al. Back EMF harmonics and DC components rejection for improvement in sensorless FTC for multiphase PMSM based on SMO with enhanced cross-coupling second-order sequence filter
Giangrande et al. Modelling of linear motor end-effects for saliency based sensorless control
JP4120775B2 (en) Vector control method and apparatus for AC motor
JP5169014B2 (en) Phase estimation method for AC motors
JP5692572B2 (en) Synchronous motor drive control device
WO2012026568A1 (en) Single-phase signal input device and system interconnection device
JP5744151B2 (en) Electric motor driving apparatus and electric motor driving method
JP3735836B2 (en) Vector control method for permanent magnet synchronous motor
JP2008199868A (en) Permanent magnet synchronous motor drive control device
JP5851662B1 (en) AC rotating machine control device
JP5828441B2 (en) Rotor phase estimator for synchronous motor
JP6150212B2 (en) Digital rotor phase speed estimation device for AC motor
Pravica et al. Interior permanent magnet wind generator torque estimation considering low-pass filter phase shift compensation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100312

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20101111

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120221

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120222

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120423

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121217

R151 Written notification of patent or utility model registration

Ref document number: 5169014

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

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

Free format text: PAYMENT UNTIL: 20160111

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20160111

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

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