JPH06321442A - 昇降機の駆動装置における回転子の時定数適合方法 - Google Patents
昇降機の駆動装置における回転子の時定数適合方法Info
- Publication number
- JPH06321442A JPH06321442A JP6063672A JP6367294A JPH06321442A JP H06321442 A JPH06321442 A JP H06321442A JP 6063672 A JP6063672 A JP 6063672A JP 6367294 A JP6367294 A JP 6367294A JP H06321442 A JPH06321442 A JP H06321442A
- Authority
- JP
- Japan
- Prior art keywords
- rotor
- time constant
- torque
- reference value
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 26
- 230000001133 acceleration Effects 0.000 claims abstract description 53
- 230000006698 induction Effects 0.000 claims abstract description 48
- 230000004044 response Effects 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 description 24
- 230000004907 flux Effects 0.000 description 21
- 238000010586 diagram Methods 0.000 description 12
- 230000008859 change Effects 0.000 description 10
- 238000005259 measurement Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000012417 linear regression Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 239000003990 capacitor Substances 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- OVSKIKFHRZPJSS-UHFFFAOYSA-N 2,4-D Chemical compound OC(=O)COC1=CC=C(Cl)C=C1Cl OVSKIKFHRZPJSS-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011423 initialization method Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B66—HOISTING; LIFTING; HAULING
- B66B—ELEVATORS; ESCALATORS OR MOVING WALKWAYS
- B66B1/00—Control systems of elevators in general
- B66B1/24—Control systems with regulation, i.e. with retroactive action, for influencing travelling speed, acceleration, or deceleration
- B66B1/28—Control systems with regulation, i.e. with retroactive action, for influencing travelling speed, acceleration, or deceleration electrical
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/16—Estimation of constants, e.g. the rotor time constant
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
- H02P2207/01—Asynchronous machines
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- Elevator Control (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
定数を回転子の実時定数に適合させる。 【構成】 トルク電流基準値の2つの値と加速度基準値
の2つの値をサンプリングし、加速度基準値とトルク電
流基準値との間の無負荷時の線形関係の勾配を求める。
この勾配を回転子の実時定数に関する理想勾配から減算
して差動信号を得て、この信号に応答して回転子の推定
時定数を得る。
Description
動機において使用する回転子の時定数決定方法に関す
る。
置の多くはベクトル制御を利用している。ベクトル制御
を用いると電動機の動特性は向上し、昇降機速度はわず
かの偏差もなく基準値にぴったり追従するようになるた
め、快適な乗り心地で移動時間を短くすることができ
る。
には、制御対象となる電動機の電気的パラメータを知る
必要がある。ベクトル制御をするということは、裏を返
せば固定子での磁束制御とは全く独立してトルク電流を
制御するということである。このような独立制御は使用
する回転子の時定数値が正しい場合にしか行うことがで
きない。K.B.Nordin、D.W.Novoth
y、D.S.Zinger共著『The Influe
nce of Motor Parameter De
viations in Feedforward F
ield Orientation Drive Sy
stems』、IEEE IAS84:22B 第52
5頁(1984年)参照のこと。残念なことに、回転子
の時定数は一定ではなく回転子の抵抗に左右され、結果
として昇降機の負荷条件によって変化する回転子の温度
に左右されることになる。回転子の時定数を推定する一
般アルゴリズムは多数あるが、これらのアルゴリズムは
極めて複雑で高価なハードウェアを必要とするか、また
は昇降機の駆動装置では得られない入出力値を使用する
かのどちらかであるため昇降機の駆動装置には適さない
のが普通である。
信号を供給してこれを利用する類いのアルゴリズムもあ
る。このアルゴリズムでは、電動機の電圧、電流および
速度を測定し、これらの測定結果から回転子の時定数を
算出する。しかしながら電動機にノイズ信号を供給する
と急激な動きが起こって乗り心地が悪くなるため、実際
には昇降機用の電動機にノイズ信号を供給することはで
きずこのようなアルゴリズムは昇降機には適さない。
R.GabrielおよびW.Leonhard著『M
icroprocessor Control of
Induction Motor』、IEEE(198
2年)を参照のこと。
特別な動作モードを利用して回転子の時定数を推定する
方式のアルゴリズムもあるが、このような動作モードも
昇降機で利用することはできない。Proc.of t
he 18th Southeastern Symp
osium on System Theory収録、
M.SongおよびJ.Mescua著『On the
Identification of Induct
ion Motor Parameters』、IEE
E、Knoxville, USA(1986年4月)
参照のこと。
ternational Conference on
Power Electronics and Va
riable Speed Drive第287頁、
H.Schierling著『Self Commis
sioning − A Novel Feature
of Modern Inverter−Fed In
duction Motor Drives』、IEE
E、London(1988年7月)参照のこと。これ
らの文献では固定子の電圧を利用して回転子の時定数を
決めることを取り上げているが、ここで問題となるのは
昇降機ではコスト上の理由から固定子の電圧を測定する
ことは殆どなく、このようなアルゴリズムも無用の長物
になってしまうということである。
導電動機の回転子の推定時定数τ2を回転子の実時定数
に適合させることにある。
にベクトル制御式誘導電動機用のトルク基準値Trefに
対する電動機トルクTMのグラフは理想的には勾配1の
直線であり、第2に同ベクトル制御式誘導回転子に供給
される回転子の推定時定数τ'2によってこの直線はグラ
フの原点を中心に回転するということに注目している。
ref直線の理想的な1:1線形関係からの変動量を求め
ることができれば、この変動量を回転子の実時定数τ2
に応じて変わる回転子の推定時定数τ'2と結び付けて考
えることができるということにある。このため、TM対
Tref直線の1:1線形関係からの変動量を利用して、
TM対Tref直線の1:1線形関係を再現することのでき
る回転子の推定時定数τ'2を求める。
ためにはTMおよびTrefを求めなければならない。T
refは決まっているので、残っているのはTMだけであ
る。従って、当然TMセンサを用いたTMの算出が次のス
テップとなる。しかしながら、TMセンサを使用せずに
TMを算出してコストを抑えることもできる。本発明の
第4の着眼点はこの点にある。すなわち、より簡単かつ
低費用で得られる値、基準加速度arefおよび負荷トル
クTLOADに関する電動機トルクTMを数式TM=Qaref
+TLOADから求めることである。この場合、TLOADは近
似的にしか求めることができないが本発明に必要なのは
負荷トルクTLOADの符号のみであるので近似値で十分で
ある。よって、TMとTrefがどのように変化するかはT
refとQaref+TLOADがどのように変化するかに基づい
て決めることができる。
準値ITrefに正比例しているという点にも注目してい
る。
Irefとaref+TLOADの間の線形関係によって表すこと
もできるのである。電動機の速度ωmが例えば0など一
定であってITrefに加速度トルク成分が含まれないよう
な場合には、トルク電流基準値ITrefを測定することに
よってTLOADを求めてもよい。加速度トルクが全くなけ
ればTLOADはITrefに正比例するからである。従って、
この直線ITref対are fの1:1関係からの勾配とオフ
セットの変化は、回転子の理想時定数τ2からの回転子
の推定時定数τ'2の変化を示すものとなる。
refのオフセット(arefは0であり、加速度トルクT
ACCELERATIONも全くないので、このオフセットはTLOAD
に正比例する)に対する同一直線の勾配から、同勾配と
回転子の推定時定数τ'2との間の線形関係を証明できる
ということにある。TLOAD>0であって昇降機が上り方
向に加速しているかまたは下り方向に減速している時に
特徴的な線形関係が得られる。また、TLOAD<0であっ
て昇降機が下り方向に加速しているか上り方向に減速し
ている時にもう1つの特徴的な線形関係が得られる。こ
のためITref対arefの勾配を判定基準として利用して
回転子の推定時定数τ'2が回転子の実時定数τ2にうま
く適合しているか否かを判定できる。無負荷状態でこの
勾配を理想勾配から減算し、これによって得られる両者
の差を比例積分回路に供給して回転子の推定時定数τ'2
を生成する。この回転子の推定時定数は、TM対Tref直
線をベクトル制御式電動機において極めて重要な1:1
線形関係に戻すことのできるものである。
機の駆動装置におけるトルク電流基準値ITrefおよび基
準加速度arefをサンプリングして判定基準を展開し、
これによって回転子の推定時定数τ'2を回転子の実時定
数τ2に適合させる。
Trefおよび基準加速度arefをサンプリングした後、I
Tref対arefの勾配を負荷トルクTLOADに正比例する同
直線のオフセットに対してプロットする。このプロット
から無負荷勾配を求めて理想無負荷勾配から減算し、得
られた差を比例積分回路に供給して回転子の実時定数τ
2に適合した回転子の推定時定数τ'2を得る。
の値を示し、下付き文字2は回転子の値を示す。
し、高性能動的プレミアムサーボ機構を用いると仮定す
ると、電動機トルクTmはベクトル制御に必要な基準ト
ルクTrefに常に一致する。従って、理想的な例におい
て以下の数式1が得られる。
御ループを任意に設定することの他に回転子の時定数τ
2についての現在の値を正確に知ることができるという
ことを基本的な前提条件とする。
をプロットすると、理想的な適合例では原点で始まって
勾配1を有する数式1のグラフを得ることができる。こ
れを図1に示す。
転子の推定時定数τ'2がうまく適合していない場合に
は、関数Tref=f(Tm)の実際のプロファイルはこの
線形特性とは異なるものとなる。図1にはこの関数の非
適合時における基本特性も示してある。
(τ'2<τ2)における回転子の推定時定数τ'2の値を
小さくすると、速度制御装置(図示せず)によって算出
される基準トルクTrefは過比例的に増加して与えられ
た電動機トルクTmに達するが、昇降機系に加わるトル
クTmは簡単に測定パラメータとして得られるような類
いのものではないため、適合を目的としてTref=f
(Tm)の関係を直接評価することは不可能である。し
かしながら(インクリメンタルシャフトエンコーダによ
って)機械の速度ωmを測定すれば求められる間接測定
パラメータすなわち角加速度ωmを利用することはでき
る。角加速度ωmとトルクTmとの関係は(摩擦を無視す
れば)以下の数式2によって得られる。
フとなる。従って、角加速度ωmを機械のトルクTmにつ
いての測定値として用いることができる。ここで、基準
速度と実速度との間の誤差は小さいと仮定すると、ωm
の数字的に好ましくない微分値に近似値ωm≒ωrefを代
入することができる。ここで、ωrefはIMの基準角加
速度である。
めるとすると、(Tm=Tref)となる。ω≒ωrefとT
refとの関係は、上述の数式2から以下の数式3によっ
て得ることができる。
トTLOADの直線になる。理想適合状態τ'2=τ2という
仮定が成り立たない場合には、直線の勾配とオフセット
が変わるだけであとは図1と同じことが言える。
して使用するωrefと、速度制御装置における回転子の
推定時定数τ'2の理想値τ2からの偏差についての基準
トルクTrefとの基本的な関係を示す。
りにITrefおよびarefを制御系の内部パラメータとし
て用いる。換算基準加速度arefは、定数の係数(aref
=ωref×rT/u)以外は基準角加速度ωrefに対応す
る。同様に、トルク形成固定子電流成分ITrefはTref
に正比例する。トルク電流基準値ITrefは、負荷トルク
TLOADに正比例する静的な部分ITstatと、加速度トル
クTACCELERATIONまたは基準加速度arefに正比例する
動的な部分ITDYNとに分けることができる。ITrefとa
refとの関係から、速度制御装置において機械の実際の
回転子の時定数τ2の値はどの程度推定値τ'2と一致し
ているか分かる。判定基準を展開し、値ITrefおよびa
refを評価することによってτ'2のチューニング用の測
定値を得る。本発明による適合方法の適用例の構造を図
3に示す。
周囲(例えばamax/2付近など)において関数ITref
=f(aref)を近似することができる。
た動作点を中心とした(小さな)変化についての基準加
速度arefと電流ITrefの基準値との間の増幅係数とな
る。オフセットBは静負荷トルクTLOADによって変化す
るが、τ2=τ'2である場合に限ってBはTLOADに正比
例する。従って、言うまでもなく「A」は速度制御装置
の回転子の時定数τ'2のチューニングの関数となる。こ
のチューニング用の測定値として、比v=τ2/τ'2を
規定する。以下、数式5に示すような関係について定量
的に検討する。
究 arefとITrefとの一般的な関係を系の動きについての
式から得る。加速度トルクTACCELERATIONについては以
下の数式6が得られる。
る誘導電動機(図4に示す等価回路)のトルクの適応部
分は以下の数式7で示す値をとる。
別の仮定のもとで、電流および加速度トルクについての
実際の値に上述の式で得られた基準値を代入する。
TDYN(数式4参照)について以下の数式8を得る。
ける電流トルク関係からITstat=TLOAD/(k
Iflux)=Bとして直接求めることができる。
数式9によって電流ITrefを算出することができる。
の微分増幅係数と解釈することができる。ここで、比適
合状態(τ'2≠τ2)の場合を無視すると、arefのプリ
ファクターAは一定ではなくなり、数式9は∂Tm/∂
ITref=kIfluxを用いて以下の数式10のように表す
ことができる。
められる動作点である。動作点では「動的」電流成分I
TDYNの平均値を用いることができる。従って、数式4か
ら得られる係数Aを以下の数式11によって算出するこ
とができる。
Aへの影響を評価できるようにするために、トルクTm
に対する非適合状態による影響を考慮した一般化トルク
−電流関係を導き出さなければならない。
ける機械についての一般的なトルク式から説明する。
成分i1βおよびi1αとの関係は以下の数式13のよう
になる。
算出できる。
機械の滑り周波数ω2を特定した。
のようになる。
御装置v=τ2/τ'2における回転子の時定数τ2のチュ
ーニングの測定値からTmのみを決定できるようにする
ために、上記数式16における電流i1αにITrefおよ
びIfluxを代入しなければならない。
φ=ω2τ2からの固定子電流の基準値と実測値の大きさ
は等しいという仮定(「高速」電流制御および整流
器)、さらにtanφについての機械の滑り周波数に関
する数式15に基づいて、以下の数式17を導き出すこ
とができる。
について以下の数式18を得ることができる。
についての数式17から以下の数式20を得ることがで
きる。
に表すことができる。
係を示す最終的な形を得ることができる。これを数式2
2に示す。
るとすれば、τ2=τ'2であって数式22はTm=kI
fluxITrefになる。
推定時定数τ'2のチューニングによるTm=f
(ITref,τ2/τ'2)の特徴線の分枝の勾配への影響
をそのまま系の適合状態についての判定基準として使用
することはできない。しかしながら、数式11を用いれ
ば、制御系の内部パラメータとして得られる基準加速度
arefと基準電流値ITrefとの関係から上述の影響につ
いて間接的に評価することは可能である。係数Aは勾配
を示し、より正確にいうならば動作点APにおける関数
ITref=f(aref)の増幅係数を示す。
ついて以下の数式23を得る。
て、数式8から得られるすでに分かっているAについて
の式で数式23をまとめると、A=Jkv/(k
Iflux)=一定が得られる。
本発明のより一層簡単な形について説明する。上述した
形の方が高精度であるが、これは以下に述べる方法に比
べて処理オーバーヘッドに対して費用をかけているがた
めである。
ク図である。速度制御装置は、速度基準値、加速度基準
値、タコメータによって得られる誘導電動機の実速度、
誘導電動機回転子の時定数に応答する。速度制御装置
は、トルク電流基準値および磁束コマンドを独立に制御
して出力することによって誘導電動機の速度を制御す
る。トルク電流基準値ITrefおよび磁束コマンドIflux
は固定子電流の2つの成分である。これらの2つの成分
は固定子電流調節器に供給される。固定子電流制御部
は、固定子の三相電流のうちのIuとIvおよび回転子の
実測角度Phiに応答し、逆変換回路に固定子の三相電
圧Uu、Uv、Uwを供給する。逆変換回路は固定子の三
相電流に応答して誘導電動機にUu、Uv、Uwを供給す
る。
関係は線形かつ1:1である。図1を参照のこと。実際
には、このようなことは有り得ない。回転子の推定時定
数τ'2の値が変化し、これに対応する電動機トルクTm
およびトルク基準値Trefの値を測定すればTm対Tref
の直線は回転していることが分かる。従って、Tm対T
ref直線の理想的な1:1関係からの変動量を決定する
ことができれば、この変動量を回転子の理想時定数τ2
に対応して変化する回転子の推定時定数τ'2と関連付け
ることができる。よって、実測Tm対Tref直線と1:1
関係との間の変動量を利用して、Tm対Tref直線を1:
1線形関係に戻すことのできる回転子の推定時定数τ'2
を生成することができる。
めにはTmとTrefとを決めなければならない。Trefは
決まっているので、残っているのはTmだけである。従
って、当然次のステップはTmセンサを用いたTmの算出
になる。しかしながら、Tmセンサを使用せずにTmを算
出してコストを抑えることもできる。すなわち、より簡
単かつ低費用で得られる値、基準加速度arefおよび負
荷トルクTLOADに関する電動機トルクTmを以下の数式
24から求めることができる。
するかはTrefとAaref+TLOADでどのように変化する
かによって判断することができる。また、Trefはトル
ク電流基準値ITrefに正比例していることも分かる。簡
単に言えば、TmとTrefとの間の線形関係はIrefとa
ref+TLOADの間の線形関係によって表すこともできる
のである。電動機の速度が例えば0など一定であってト
ルク電流基準値ITrefに加速度トルク成分が含まれない
ような場合には、トルク電流基準値ITrefを測定するこ
とによってTLOADを求めてもよい。加速度トルクが全く
なければ、必要なトルクは負荷トルクTLOADに対するも
のだけになるため、負荷トルクTLOADはトルク電流基準
値ITrefに正比例する。定速では加速度トルクT
ACCELERATIONは全く必要ない。従って、この直線勾配と
オフセットの変化は、理想的な回転子の時定数τ2から
の回転子の推定時定数τ'2の変化を示すものとなる。図
7を参照のこと。
転(図1)と全く同じように選択する回転子の推定時定
数τ'2によって回転する。図7におけるオフセットB
は、加速度がなく加速度トルクも全くない状態ではトル
ク電流基準値ITrefの値と等しい。このオフセットB
は、数式5において示すようにτ2=τ'2である場合に
限って負荷トルクTLOADに正比例する。しかしながら、
BはTLOADと同一の符号を有する。aref=0であって
Tm=TLOADである場合には数式22および数式9を参
照のこと。τ2=τ'2であるとすると、TLOAD=KI
flux×ITrefになり、ITref=Bになる。
に対してプロットすると、図8および図9に示すグラフ
が得られる。この勾配Aは数式4の勾配Aと全く同じで
ある。図8は、回転子の推定時定数τ'2の様々な値につ
いてのかごの上り方向への加速時または下り方向への減
速時に得られる勾配A対オフセットBのグラフである。
図9は、回転子の推定時定数τ'2の様々な値についての
かごの下り方向への加速時または上り方向への減速時に
得られる勾配A対オフセットBのグラフである。負荷ト
ルクTLOADが0より大きい図8において、回転子の推定
時定数τ'2と勾配Aとの関係は特徴的なものとなる。図
9において、負荷トルクTLOADは0未満であり、回転子
の推定時定数τ'2と勾配Aとの関係は特徴的なものとな
る。よって、勾配AとオフセットBとの関係は、負荷ト
ルクTLOADが0より大きい場合と負荷トルクTLOADが0
未満である場合について、それぞれ以下の数式25およ
び26で表すことができる。
下り方向への減速時の図7に示す直線の勾配であり、A
2は下り方向への加速時または上り方向への減速時の図
7に示す直線の勾配;A0は回転子の推定時定数τ2=
τ'2である理想的な適合時の図7に示す直線の勾配;C
1およびC2はそれぞれ図8および図9に示す直線の線形
部分の勾配である。
に関連している。この慣性は上述の数式TM=Aaref+
TLOADにおけるAに関連している。すなわち、T=J×
dωr/dt+TLOAD(回転)である。ここでJは電動
機の駆動系の回転慣性である。この第2の式は上記の数
式の一般例である。
AとJとは等しくならない。誘導電動機の回転慣性に正
比例するのでC1およびC2はほぼ等しいと考えることが
できる。
じ直線のオフセットに対するグラフは、上り方向への加
速時または下り方向への減速時で負荷トルクが0より大
きい場合および下り方向への加速時または上り方向への
減速時で負荷トルクが0未満である時の勾配と回転子の
推定時定数τ'2との特徴的な関係を示す直線となる。従
って、この勾配を判定基準として利用して回転子の推定
時定数τ'2が回転子の実時定数τ2にうまく適合してい
るか否かを判定できる。勾配A0を理想勾配A0,idealか
ら減算し、これによって得られる両者の差を比例積分回
路に供給して回転子の推定時定数τ'2を生成する。この
回転子の推定時定数は、Tm対Tref直線を1:1線形関
係に戻すことのできるものである。
機能的に図6と同様の速度制御装置は、速度基準値、加
速度基準値、電動機速度、回転子の推定時定数に応答し
て、磁束とトルクに関連する固定子電流の2つの成分を
独立に制御する。速度制御装置は10ミリ秒のクロック
で動作して加速度基準値をサンプリングするとともに、
これに対応するトルク電流基準値ITrefを算出する。こ
の部分は本発明のなかで昇降機の動作時に実行しなけれ
ばならない唯一の部分である。
および加速度基準値arefの値は、ITrefとarefとの線
形関係の勾配AとオフセットBを算出する線形回帰ブロ
ックに供給されるまで格納される。A1は昇降機の上り
方向への加速時または下り方向への減速時における勾配
であり、A2は昇降機の下り方向への加速時または上り
方向への減速時における勾配である。Bは加速と減速の
両方について同一である。線形回帰ブロックは、ITref
対aref直線の無負荷状態での勾配を算出するブロック
にA1、A2、Bを供給する。このブロックは、まず加速
または減速の判定基準を使用すべきであるか否かを図8
および図9に基づいて決定する。表2はBの符号(=負
荷トルクTLOADの符号)、昇降機の移動方向および加速
度基準値とトルク電流基準値とのサンプル、加速相また
は減速相、使用する勾配の組み合わせを示す。
を算出する。C1およびC2は、線形回帰ブロックや、A
1、A2およびBからA0を算出するブロックでは得られ
ない。C1およびC2は後述するような適用時において得
られるものである。A0は次にノイズ抑制用のデジタル
フィルタに供給される。この種のデジタルフィルタはP
hilips著『Digital Control S
ystemAnalysis and Design』
第12章に記載されている。次に、後述するような形で
得られる回転子の時定数の実際の値に対応する理想値A
0,idealからITref対arefの直線の実際の無負荷勾配A
0を減算する。これによって得られた差を上述した文献
の第8章に記載されているような比例積分制御装置に供
給する。加算器からの差が負である場合には、回転子の
推定時定数τ′2は小さすぎるので比例積分制御装置か
ら供給される回転子の推定時定数τ′2の初期値を大き
くするかまたはその逆にする。回転子の推定時定数τ′
2を増減することによって、Tm対Trefを1:1に戻す
ことができ、図6に示すベクトル制御によって磁束とト
ルクとを独立して制御することができる。
0,idealの決定にもからんでいる。
ロック図である。まず、負荷トルクTLOADをゼロに設定
すなわちかご室の重量とかごの負荷を釣り合い重りと釣
り合わせる。例えば誘導電動機用の識別アルゴリズムな
どを利用して回転子の推定時定数τ′2をτ2の実際の値
に設定する。この実際の値は例えばH.Schierl
ing著、Proc.of Third Intern
ational Conference on Pow
er Electronics and Variab
le Speed Drive第287頁、『Self
−commissioning−−A Novel F
eature of Modern Inverter
−fed Induction Motor Driv
es』IEEE、London(1988年7月)など
に記載されているような自己規定手順において得ること
ができる。かごを移動させ、トルク電流基準値ITrefと
加速度基準値arefとをサンプリングする。Bは0であ
るので、数式25と26の間には差はなく、A0,ideal
は線形回帰ブロックから供給される。
想適合状態下にあるベクトル制御装置の理想値τ2に対
するブロック図である。ここで2つの初期化ランを取
る、すなわち、空のかごを0未満の負荷トルクで上昇さ
せてC2を決め、空のかごを0より大きい負荷トルクで
下降させてC1を決定する。トルク電流基準値ITrefお
よび加速度基準値arefをサンプリングし、A0,idealと
共にA1、A2、Bを算出ブロックに供給し、C1=(A1
−A0,ideal)/B(負荷トルクTLOADは0より大き
い)およびC2=(A2−A0,ideal)/B(負荷トルク
TLOADは0未満)を得る。
び固定子温度を昇降機の動作回数に対してプロットした
図である。回転子の時定数τ′2は回転子の温度に正比
例するため、この図は本発明が回転子の時定数を正しく
推定していることを証明する良い例となる。固定子の温
度は、定常状態では、回転子の温度と同じである。ま
た、回転子の温度を変えるためには固定子の温度を回転
子の温度の特性についての基準として用いることができ
る。
アルゴリズムに代わるアルゴリズムを以下に示す。
ontrol ofElectrical drive
s』、Springer−Verlag(Spring
er出版)、Berlin, Heidelberg,
New York, Tokyo(1985年)に記
載されているような周知の誘導電動機モデルとは若干異
なる誘導電動機モデルに基づいている。以下の数式は本
発明で使用するモデルを示すものである。これらの数式
から、軸(a)および(b)を有する座標系を使用した
固定子固定座標における静止状態での誘導電動機がどの
ようなものか分かる。ここで、(a)は直線R、S、T
を有する三相電動機についての相のうちの1つの電動機
相Rを示す。
動制御装置を調節するために利用することもできる。数
式27および28に示されるように、完全モデルを説明
するために使用するパラメータは4つある。これらのパ
ラメータは、固定子の抵抗R1、回転子の時定数τ2、全
漏れインダクタンスLσ、そして従来のモデルでは使わ
れていないがこれら3つのパラメータから誘導できるパ
ラメータLRである。LR=LH 2/R2であって、LHはメ
インインダクタンス、R2は回転子の抵抗である。等価
回路を図4に示す。
発明の一部ではないがH.Schierling著、博
士論文1987第61〜65頁『Self−Commi
ssioning−−A Novel Feature
of Modern Inverter−Fed I
nduction Motor Drives』、Si
emens AG(西ドイツ)およびH.Schier
ling著『Selbsteinstellendes
and selbstanpassendes An
triebsregelsystem fur die
Asynchronmaschine mit Pu
lsqechselrichter』(Technic
al University of Darmstad
t,西ドイツ)に記載されているアルゴリズムを利用し
て識別する。
と、以下の数式29が得られる。
電流の偏差から全漏れインダクタンスを算出することが
できる。さて、低電圧域におけるパルス制御式逆変換回
路の上述したような誤差をなくすために、>100Vの
範囲の高電圧を電動機に印加しなければならないので実
際の設定電圧は所望の値に対応する。電流は極めて急速
に高くなるので測定時間は短く、Ψ2a B=0はほぼ満た
される。もう1つの条件であるi1a=0を満たすため
に、電流のゼロ交点で測定を行う。
示すように印加電圧とこの電圧によって生じる電動機電
流の存在する時間内になる。
動機に印加される。時刻t1において電流は予め定めら
れた限界値に達する。この限界値は電動機の定格電流に
よって変わるが、安全上の理由から0.5Iratedまで
に制限されている。続いて電圧は0まで降下し、誘導機
の電流は周波数変換器の中間コンデンサを介して放電さ
れる。コンデンサにおける電圧の上昇は吸収エネルギす
なわち電動機の吸収エネルギE=0.5i2Lおよびコ
ンデンサの吸収エネルギE=0.5(U+dU)2Cに
基づいて算出することができる。
U=540Vであり、全漏れインダクタンスは、ここで
用いた電動機におけるC=1mFで最大30mHにな
る。また、電圧上昇Ud=2.78V=0.5%が得ら
れる。これは十分許容できる値である。
動機に印加される。電流はもう1つの予め定められた負
の限界値まで降下する。これによって電圧は再び0にな
り、インダクタは放電し、全漏れインダクタンスの識別
は終了する。時刻t2とt3の間に、電動機電流のゼロ交
点が検出される。このゼロ交点近辺のある時間dtにつ
いて電動機の電流を決定し、そこから全漏れインダクタ
ンスLσ=U2dt/di2を求める。
Rを求めるための識別アルゴリズムについても説明す
る。このアルゴリズムは、電動機の固定子電圧および固
定子電流が全て定常状態にある時に効を奏する。固定子
の電流および電圧が定常値にあるということは、単に電
流および電圧を印加してこれがなくなるまで遷移させる
には十分であると思われるだけの時間待っているだけで
簡単に確認することができる。
て電圧フェーザ(phaser)を予め設定することができ
る。従って、周波数変換器に供給される定格電圧フェー
ザを適当に変化させることで識別に必要な周波数ωの出
力周波数変換器ac電圧U1aを達成することができる。
しかしながら、周波数変換器内で発生する作用(最小起
動時間、むだ時間)の影響で、実際に設定される出力電
圧対所望の公称値の基本波に位相のシフトが起こる。低
電圧では、この位相シフトはここで使用した周波数変換
器については5〜10°に達する。したがって、くどい
ようではあるが公称電圧を識別用の基準値として使用す
ることは不可能である。
動機に図15に示すように正弦電圧U1aではなく方形波
を供給すればよい。
機が静止状態にある時に供給する。電動機がトルクを発
生せず止まったままでいるようにするために、もう1つ
の固定子電圧成分U1bはゼロに設定しておく。固定子の
電圧曲線は基本波ωと高調波とからなる。ωは時間T
(図1参照)からω=2π/Tによって導かれる。
られる電動機の固定子電流i1aを示す。固定子電流i1a
の曲線はu1aの正負によって増減される指数関数からな
る。
部を示す。角度φは電動機の固定子電圧の基本周波数と
固定子電流との間の位相角を意味する。tanφは本発
明のさらに他の部分で使用する。逆変換回路(3)は基
本波ω1(4)の所望の固定子電圧を誘導電動機(5)
に供給する。固定子電圧は図1に示すものと同一であ
る。A/D変換器(7)を利用して時間Tで電動機の固
定子電流(6)をサンプリングしてブロック(8)およ
び(9)に供給する。ブロック(8)は、以下の数式3
0に基づいて値ycを算出する。
ンプリングした値である。次にサンプリングされる値は
i1a[k+1]というふうになっていく。最初の値はi
1a[0]である。最後の値はi1a[T/tclock−1]
である。ブロック(8)の出力はycosであり、ycの最
終値に等しい。すなわち、以下の数式31が得られる。
間的に連続した式の離散形態になっている。
と同じようにして、数式33、34、35に基づいてy
sin の値を算出する(数式30、31、32参照)。
ックは数式36に基づいてtanφを算出する。
で、nは図1において説明したような固定子の電圧曲線
を電動機に供給する毎に得られるサンプルの数である。
基本周波数ωの値は毎回変化する。従って、全繰り返し
の結果は、基本周波数ωと適当なtanφの値のn対か
らなる集合であり、以下の数式37のように表すことが
できる。
よび数式37を使用して残りの電動機のパラメータ
R1、LRおよびτ2を推定するアルゴリズムを示すもの
である。定常状態の場合には、電動機モデル(数式2
7、28)から数式38のようにして固定子の複素イン
ピーダンスZ(ω)を導き出すことができる。
想成分は以下の数式39で示すようになる。
なる。
はゼロである(この量は実数である)ので、以下の数式
41が得られる。
よびX(ω)を数式27および28のパラメータに置き
換えると、非線形式として以下の数式42が得られる。
すようにして得たものである。
て書き替えることができる。書き替えによって得られる
結果は以下の数式58に示すようなn個の非線形式にな
る。
σ、tanφおよびωは分かっており、パラメータ
R1、LR、τ2が抜けている。これらの足りないパラメ
ータは数式42を満たさなければならない。実際には、
Lσおよびtanφには測定誤差があるため数式42に
示す条件f(R1,LR,τ2)=0を正確に満たすこと
は不可能である。
械モデルの設定時の抜けなども原因となる。従って、足
りないパラメータは数式42を正確に満足させるのでは
なく近似させるのである。有効な近似値についての判定
基準はf(R1,LR,τ2)の二乗の最小化である(最小
二乗法)。すなわち、以下の数式59のようになる。
集合(R1,LR,τ2)[k]から新しいパラメータ集
合(R1,LR,τ2)[k+1]を算出するための繰り
返しパターンを構築することができる。推定した初期値
(R1,LR,τ2)[0]がR1、LR、τ2を近似するの
に適したものであれば、このアルゴリズムは足りないパ
ラメータR1,LR,τ2に収束する。
昇降機の電動機駆動装置に供給する。ベクトル制御用に
は、回転子の時定数を速度制御装置に供給する。
囲を逸脱することなく様々な修正を加えることができ
る。
以下のような効果が得られる。
て移動時間を短縮し、かつ快適な乗り心地を達成でき
る。
時定数を得ることができる。
余計な移動時間の原因となる時間的な遅れはなくなる。
にほんのわずかな計算時間を加えるだけでよく、殆どの
計算は静止時に行うことができる。
る。
図である。
図である。
(℃)対昇降機の数を示すグラフである。
別に使用する固定子の電流および電圧波形を示すグラフ
である。
パラメータLRを識別するために使用する固定子の電圧
を示すグラフである。
パラメータLRを識別するために使用する固定子の電流
を示すグラフである。
ある。
図である。
Claims (15)
- 【請求項1】 トルク電流基準値の2つの値と加速度基
準値の2つの値をサンプリングし、 前記加速度基準値と前記トルク電流基準値との間の無負
荷時の線形関係の勾配を求め、 前記勾配を回転子の実時定数に関する理想勾配から減算
して差動信号を得、 前記差動信号に応答して回転子の推定時定数を得る昇降
機駆動用誘導電動機回転子の推定時定数供給方法。 - 【請求項2】 昇降機の駆動装置が回転子の時定数の実
際の値に応答する時には、前記理想勾配は前記勾配の値
になるため、前記理想勾配は前記回転子の実時定数に関
連している請求項1記載の誘導電動機用回転子の推定時
定数供給方法。 - 【請求項3】 電動機トルクを得、 トルク基準値を得、 前記電動機トルクと前記トルク基準値との間の線形関係
の勾配の回転子理想時定数に関する理想勾配からの変動
量を示す信号を得、 前記信号に応答して回転子の推定時定数を得る昇降機駆
動用誘導機回転子の推定時定数供給方法。 - 【請求項4】 トルク電流基準値によって前記トルク基
準値を近似的に得る請求項1記載の誘導電動機用回転子
の推定時定数供給方法。 - 【請求項5】 トルク電流基準値によって前記トルク基
準値を近似的に得る請求項3記載の誘導電動機用回転子
の推定時定数供給方法。 - 【請求項6】 加速度基準値および負荷トルクの符号に
応答して前記電動機トルクを近似的に得る請求項1記載
の誘導電動機用回転子の推定時定数供給方法。 - 【請求項7】 加速度基準値および負荷トルクの符号に
応答して前記電動機トルクを近似的に得る請求項3記載
の誘導電動機用回転子の推定時定数供給方法。 - 【請求項8】 前記勾配は誘導機が無負荷状態の時に得
られる請求項1記載の誘導電動機用回転子の推定時定数
供給方法。 - 【請求項9】 前記勾配は誘導機が無負荷状態の時に得
られる請求項3記載の誘導電動機用回転子の推定時定数
供給方法。 - 【請求項10】 昇降機の駆動装置が回転子の時定数の
実際の値に応答する時には、前記理想勾配は前記勾配の
値になるため、前記理想勾配は前記回転子の実時定数に
関連している請求項1記載の誘導電動機用回転子の推定
時定数供給方法。 - 【請求項11】 昇降機の駆動装置が回転子の時定数の
実際の値に応答する時には、前記理想勾配は前記勾配の
値になるため、前記理想勾配は前記回転子の実時定数に
関連している請求項3記載の誘導電動機用回転子の推定
時定数供給方法。 - 【請求項12】 負荷トルクが0より大きい時には前記
信号は第1の値となり、負荷トルクが0未満の時には前
記信号は第2の値となる請求項1記載の誘導電動機用回
転子の推定時定数供給方法。 - 【請求項13】 負荷トルクが0より大きい時には前記
信号は第1の値となり、負荷トルクが0未満の時には前
記信号は第2の値となる請求項3記載の誘導電動機用回
転子の推定時定数供給方法。 - 【請求項14】 トルク電流基準値をサンプリングする
ための手段と、 加速度基準値をサンプリングするための手段と、 前記誘導機が無負荷状態の時の前記トルク電流基準値と
前記加速度基準値との間の線形関係の勾配を得るための
手段と、 前記無負荷状態の勾配を回転子の実時定数に関連した理
想無負荷勾配から減算して差動信号を供給する加算器
と、 前記差動信号に応答して回転子の推定時定数を供給する
ための手段と、を有する誘導機用回転子の時定数供給装
置。 - 【請求項15】 前記手段は比例積分回路である請求項
14記載の装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US040,038 | 1993-03-31 | ||
| US08/040,038 US5476158A (en) | 1993-03-31 | 1993-03-31 | Rotor time constant adaptation for induction motor in vector controlled elevator drive |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06321442A true JPH06321442A (ja) | 1994-11-22 |
| JP3606899B2 JP3606899B2 (ja) | 2005-01-05 |
Family
ID=21908750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP06367294A Expired - Fee Related JP3606899B2 (ja) | 1993-03-31 | 1994-03-31 | 昇降機の駆動装置における回転子の時定数適合方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US5476158A (ja) |
| JP (1) | JP3606899B2 (ja) |
| DE (1) | DE4411128A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8917040B2 (en) | 2012-02-22 | 2014-12-23 | Denso Corporation | AC motor control apparatus |
Families Citing this family (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3230364B2 (ja) * | 1994-03-03 | 2001-11-19 | 株式会社明電舎 | 速度制御系における速度推定演算処理方法 |
| KR0186121B1 (ko) * | 1995-11-23 | 1999-04-15 | 이종수 | 엘리베이터의 진동보상을 위한 속도 제어장치 |
| US5929400A (en) * | 1997-12-22 | 1999-07-27 | Otis Elevator Company | Self commissioning controller for field-oriented elevator motor/drive system |
| US5880416A (en) * | 1997-12-22 | 1999-03-09 | Otis Elevator Company | Automatic calibration of motor speed loop gain for an elevator motor control |
| US5880415A (en) * | 1997-12-22 | 1999-03-09 | Otis Elevator Company | Automatic calibration of current regulator control compensation for an elevator motor drive with locked rotor |
| US5883344A (en) * | 1997-12-22 | 1999-03-16 | Otis Elevator Company | Automatic calibration of field-oriented elevator motor drive parameters using standstill motor measurements |
| US5896954A (en) * | 1997-12-22 | 1999-04-27 | Otis Elevator Company | Automatic fine tuning of rotor time constant in field-oriented elevator motor drive |
| US6316904B1 (en) | 2000-06-27 | 2001-11-13 | Ford Global Technologies, Inc. | Speed and rotor time constant estimation for torque control of an induction motor |
| JP2004173422A (ja) * | 2002-11-20 | 2004-06-17 | Fanuc Ltd | 誘導電動機の制御装置 |
| DE102004033074A1 (de) * | 2004-07-08 | 2006-02-02 | Dieter Miehlich | Trainingsgerät |
| ATE412915T1 (de) * | 2005-06-02 | 2008-11-15 | Abb Oy | Verfahren zur ermittlung der rotorzeitkonstanten einer asynchronmaschine |
| US20070105672A1 (en) * | 2005-10-18 | 2007-05-10 | Daren Luedtke | Variable speed transmission |
| JP2007164362A (ja) * | 2005-12-12 | 2007-06-28 | Fanuc Ltd | 時定数自動推定方法及び時定数自動推定装置並びにサーボモータ制御システム |
| US20080272596A1 (en) * | 2007-05-02 | 2008-11-06 | House Edward T | Wind turbine variable speed transmission |
| US8115441B2 (en) * | 2007-07-19 | 2012-02-14 | Hamilton Sundstrand Corporation | On-line measurement of an induction machine's rotor time constant by small signal d-axis current injection |
| US7769552B2 (en) * | 2008-05-16 | 2010-08-03 | Schneider Electric USA, Inc. | Method and apparatus for estimating induction motor rotor temperature |
| US8102140B2 (en) * | 2008-05-16 | 2012-01-24 | Schneider Electric USA, Inc. | Method and apparatus for estimating induction motor electrical parameters |
| FI124592B (fi) * | 2013-06-20 | 2014-10-31 | Kone Corp | Menetelmä ja laitteisto hissin sähkömoottorin ohjaamiseksi |
| US10698031B2 (en) | 2018-02-02 | 2020-06-30 | R. Gene Smiley | Systems and methods for obtaining and validating performance data on motors |
| CN120928212B (zh) * | 2025-10-14 | 2026-03-17 | 上海大恩芯源微电子有限公司 | 电池风险确定方法、装置、设备及介质 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6032440B2 (ja) * | 1979-12-13 | 1985-07-27 | 富士電機株式会社 | 誘導電動機の磁束ベクトル演算器 |
| US4677360A (en) * | 1986-03-13 | 1987-06-30 | General Electric Company | Field weakening induction drive |
| EP0274716A1 (de) * | 1987-01-09 | 1988-07-20 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zum Bestimmen des Flussvektors einer Drehfeldmaschine |
| US5076399A (en) * | 1990-09-28 | 1991-12-31 | Otis Elevator Company | Elevator start control technique for reduced start jerk and acceleration overshoot |
| US5334923A (en) * | 1990-10-01 | 1994-08-02 | Wisconsin Alumni Research Foundation | Motor torque control method and apparatus |
-
1993
- 1993-03-31 US US08/040,038 patent/US5476158A/en not_active Expired - Lifetime
-
1994
- 1994-03-30 DE DE4411128A patent/DE4411128A1/de not_active Withdrawn
- 1994-03-31 JP JP06367294A patent/JP3606899B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8917040B2 (en) | 2012-02-22 | 2014-12-23 | Denso Corporation | AC motor control apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| DE4411128A1 (de) | 1995-04-06 |
| JP3606899B2 (ja) | 2005-01-05 |
| US5476158A (en) | 1995-12-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3606899B2 (ja) | 昇降機の駆動装置における回転子の時定数適合方法 | |
| JP2780263B2 (ja) | 誘導電動機のベクトル制御方法と装置 | |
| KR100766665B1 (ko) | 유도 전동기의 속도 추정치 보정 방법 및 그 장치 | |
| JPH07110160B2 (ja) | 誘導電動機の制御装置 | |
| JPH06324127A (ja) | N個の誘導機電気パラメータを得る同定方法 | |
| JPS6042712B2 (ja) | 非同期機の運転装置 | |
| JP7055242B2 (ja) | 電力変換装置の制御装置、制御方法、及び電動機駆動システム | |
| EP0490024B1 (en) | Induction motor vector control | |
| JP7072728B2 (ja) | 電力変換装置の制御装置、及び電動機駆動システム | |
| JP2585376B2 (ja) | 誘導電動機の制御方法 | |
| US4510430A (en) | Vector control method and system for an induction motor | |
| EP1729405A1 (en) | Speed control apparatus of vector controlled alternating current motor | |
| JP2929344B2 (ja) | 電動機定数測定方法及びその装置 | |
| JP3099159B2 (ja) | 電動機定数測定方法及び装置 | |
| JPH10248300A (ja) | 電力変換装置 | |
| JP2816103B2 (ja) | 誘導電動機の制御装置 | |
| JP2579119B2 (ja) | 誘導電動機のベクトル制御装置 | |
| JP3229216B2 (ja) | ビートレス制御装置 | |
| JPH09238492A (ja) | 交流電動機の制御装置 | |
| JPS63185789A (ja) | エレベ−タ−の制御方法及び装置 | |
| JPS6159071B2 (ja) | ||
| JP3636901B2 (ja) | 誘導電動機の制御装置 | |
| JP3609098B2 (ja) | 誘導電動機のベクトル制御装置におけるモータ定数同定方法 | |
| JP2634959B2 (ja) | 速度センサレス速度制御方式 | |
| JP2002051594A (ja) | 主軸制御装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20040106 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20040109 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040218 |
|
| 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: 20040928 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041006 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081015 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081015 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091015 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091015 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101015 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |