JPH0484201A - Automatic tuning controller - Google Patents

Automatic tuning controller

Info

Publication number
JPH0484201A
JPH0484201A JP19838990A JP19838990A JPH0484201A JP H0484201 A JPH0484201 A JP H0484201A JP 19838990 A JP19838990 A JP 19838990A JP 19838990 A JP19838990 A JP 19838990A JP H0484201 A JPH0484201 A JP H0484201A
Authority
JP
Japan
Prior art keywords
tuning
plant
control
auto
attenuation rate
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
Application number
JP19838990A
Other languages
Japanese (ja)
Other versions
JP2865109B2 (en
Inventor
Shinichi Kitagawa
伸一 北川
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.)
Shimadzu Corp
Original Assignee
Shimadzu 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 Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP2198389A priority Critical patent/JP2865109B2/en
Publication of JPH0484201A publication Critical patent/JPH0484201A/en
Application granted granted Critical
Publication of JP2865109B2 publication Critical patent/JP2865109B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
(57) [Summary] This bulletin contains application data before electronic filing, so abstract data is not recorded.

Description

【発明の詳細な説明】 (イ)産業上の利用分野 この発明は、プロセス制御で用いられる調節計、特に、
その制御パラメータ(P I D定数)を自動調節する
オートチューニング調節計に関する。
DETAILED DESCRIPTION OF THE INVENTION (a) Industrial Application Field This invention relates to a controller used in process control, in particular,
This invention relates to an auto-tuning controller that automatically adjusts its control parameters (PID constants).

(ロ)従来の技術 従来、制御パラメータを自動調節する方法として、 ■プラントのパルス応答やインデシャル応答の定数から
ジーグラ・ニコルスの公式で決定する方法。
(B) Conventional technology Conventionally, as a method for automatically adjusting control parameters, there are two methods: - Determination using the Ziegler-Nichols formula from the constants of the plant's pulse response and individual response.

■比例制御で振動が減衰しない限界振動を発生させ、そ
の振動周期等をジーグラ・ニコルスの公式に当てはめて
決定するいわゆる限界感度法。
■The so-called limit sensitivity method uses proportional control to generate a limit vibration that does not dampen, and determines its vibration period by applying the Ziegler-Nichols formula.

■プラントに同定信号を印加し、得られる応答を統計的
に処理してパラメータを推定し、それに適するパラメー
タを決定する方法。
■A method of applying identification signals to the plant, statistically processing the responses obtained, estimating parameters, and determining appropriate parameters.

■制御系に非線型要素を挿入して振動を発生させて決定
する、いわゆるリミットサイクル法。
■The so-called limit cycle method, which determines by inserting nonlinear elements into the control system and generating vibrations.

等がある。etc.

(ハ)発明が解決しようとする課題 上記した従来の■の方法はパラメータの調整時に、制御
系を一旦開ループにする必要があるので非制御状態とな
る上、ステップやパルス入力を加えるので外乱となって
系を乱す。■の方法はパラメータ調整時に限界振動を発
生するので制御系を乱す上、限界振動はプラント応答の
発散寸前であり、安全上も好ましくない。■の方法はパ
ラメータ決定のためのプラントの同定を正確にするのに
、大振幅の同定信号を長時間印加する方が有利であるが
、そのために大きな外乱を与えるとう矛盾がある。■の
方法は■の方法と同様の問題点がある。
(c) Problems to be Solved by the Invention In the conventional method (2) described above, when adjusting parameters, it is necessary to open the control system once, resulting in an uncontrolled state, and addition of step and pulse inputs causes disturbances. and disturbs the system. Method (2) generates limit vibrations during parameter adjustment, which disturbs the control system, and limit vibrations are on the verge of divergence of the plant response, which is unfavorable from a safety standpoint. In the method (2), it is advantageous to apply a large-amplitude identification signal for a long time in order to accurately identify the plant for parameter determination, but there is a contradiction in that this causes a large disturbance. Method (2) has the same problems as method (2).

この発明は、上記各問題点に着目してなされたものであ
って、プラントを非制御状態とすることなく、またプラ
ントに余分な外乱を与えることなく、通常の制御状態で
安定して調整できるオートチューニング調節計を提供す
ることを目的としている。
This invention was made by focusing on each of the above-mentioned problems, and it is possible to stably adjust the plant in a normal control state without putting the plant in an uncontrolled state or giving any unnecessary disturbance to the plant. The purpose is to provide an auto-tuning controller.

(ニ)課題を解決するための手段及び作用この発明のオ
ートチューニング調節計は、被制御出力を計測する手段
と、被制御出力の振幅減衰率を算出する手段と、行き過
ぎ量を算出する手段と、算出された振幅減衰率と行き過
ぎ量とから制御パラメータを更新するための更新係数を
算出する手段と、算出された更新係数でパラメータを更
新する手段を持ち、PIまたはPID制御を実行しなが
ら、更新係数が所定値となるまで振幅減衰率算出、行き
過ぎ量算出、更新係数算出、パラメータ更新を繰り返す
調整手段から構成されている。
(d) Means and operation for solving the problems The auto-tuning controller of the present invention includes means for measuring the controlled output, means for calculating the amplitude attenuation rate of the controlled output, and means for calculating the overshoot amount. , a means for calculating an update coefficient for updating a control parameter from the calculated amplitude attenuation rate and an overshoot amount, and a means for updating the parameter with the calculated update coefficient, and while executing PI or PID control, It is comprised of an adjusting means that repeats amplitude attenuation rate calculation, excess amount calculation, update coefficient calculation, and parameter update until the update coefficient reaches a predetermined value.

この発明のオートチューニング調節計は、プラントを非
制御状態とすることなく、PIまたはPID制御を実行
しながらプラントの応答を計測し、制御中に起こる種々
の外乱で生じる振動の振幅減衰率の算出、行き過ぎ量の
算出、振幅減衰率と行き過ぎ量からのパラメータ更新係
数の算出、更新係数による制御パラメータ(比例帯およ
び積分時間)の更新を行う。更新係数が1となるまで以
上の手順を繰り返すことにより、希望するプラント応答
を得る制御パラメータが決定される。応答の計測は、振
動のピークが一定時間以上こなければ打ち切る。また応
答の計測中に振幅減衰率が1以上となることを予測する
と、ただちにオートチューニングを中止し、予め定めた
安全値へ制御パラメータを更新する。
The autotuning controller of this invention measures the response of the plant while executing PI or PID control without bringing the plant into a non-controlled state, and calculates the amplitude attenuation rate of vibrations caused by various disturbances that occur during control. , calculates the amount of overshoot, calculates a parameter update coefficient from the amplitude attenuation rate and the amount of overshoot, and updates the control parameters (proportional band and integral time) using the update coefficient. By repeating the above procedure until the update coefficient becomes 1, control parameters that obtain the desired plant response are determined. Response measurement is stopped if the peak of vibration does not occur for a certain period of time or longer. Furthermore, if it is predicted that the amplitude attenuation rate will be 1 or more during response measurement, auto-tuning is immediately stopped and the control parameters are updated to predetermined safe values.

被制御量の計測の対象は、通常の制御状態で種々の外乱
により往じる振動であり、プラントの同定や制御パラメ
ータの算出のために信号を印加してプラントに余計な外
乱を与えることがない。
The target of controlled variable measurement is the vibration that occurs due to various disturbances in normal control conditions, and it is possible to apply signals to identify the plant and calculate control parameters to cause unnecessary disturbances to the plant. do not have.

チューニングは常に制御状態で行われるためプラントが
安定に保たれる。応答の計測は一定時間以内にピークが
こなければ打ち切られるので必要以上の長時間を要しな
い。チューニング中にプラントの応答が万一発散の方向
に向かうと、振幅減衰率が1以上であることを予測した
時点で直ちにオートチューニングを中止し、制御パラメ
ータを安全値に変更するので、プラントの安全性が保た
れる。
Tuning is always performed under control, so the plant remains stable. Response measurement is aborted if a peak does not occur within a certain period of time, so it does not take an unnecessarily long time. If the response of the plant should diverge during tuning, auto-tuning will be immediately stopped as soon as the amplitude attenuation rate is predicted to be 1 or more, and the control parameters will be changed to safe values, ensuring plant safety. Maintains sexuality.

(ホ)実施例 以下、実施例により、この発明をさらに詳細に説明する
(E) Examples The present invention will be explained in more detail with reference to Examples below.

第1図は、この発明の一実施例を示すオートチューニン
グ調節計の構成を示すプロ・7り図である。このオート
チューニング調節計1は、比例演算(P)部2、積分演
算(1)部3、微分演算(D)部4及びオートチューニ
ングアルゴリズム(制御)部5から成り、比例演算部2
、積分演算部3、微分演算部4でPID制御を実行しな
がらプラント6の応答を計測し、プロセス量(P V)
、設定量(SV)から求まる偏差(y)を基に制御パラ
メータ(比例帯P、積分時間TI)を自動的に調節する
FIG. 1 is a schematic diagram showing the configuration of an auto-tuning controller according to an embodiment of the present invention. This auto-tuning controller 1 consists of a proportional calculation (P) section 2, an integral calculation (1) section 3, a differential calculation (D) section 4, and an auto-tuning algorithm (control) section 5.
, the integral calculation unit 3 and the differential calculation unit 4 measure the response of the plant 6 while executing PID control, and calculate the process amount (PV).
, control parameters (proportional band P, integral time TI) are automatically adjusted based on the deviation (y) found from the set value (SV).

第2図で、プラントの閉ループ応答(被制御量の偏差の
波形)と振幅減衰率aと行き過ぎ量すを定義する。すな
わち振幅減衰率aは、 yゴ −yt ’lr    Vz 行き過ぎ量すは、b””  Y z / ’I l 、
最大待ち時間をTとして、T+’ 、Tt 、TsはT
以内である。
In FIG. 2, the closed loop response of the plant (the waveform of the deviation of the controlled variable), the amplitude attenuation rate a, and the overshoot amount are defined. That is, the amplitude attenuation rate a is: ygo -yt 'lr Vz The amount of overshoot is b"" Yz / 'I l ,
When the maximum waiting time is T, T+', Tt, and Ts are T
within

yl、y2、y3はそれぞれ第1〜第3ピークのyの値
を示し、正負の符号を持つ。yoはオートチューニング
の処理が開始するための条件を与えるyの正の値で、y
の絶対値がyoを越えるとオートチューニングの処理を
開始する。ただしオートチューニング開始から第1ピー
クまでの時間TI、あるいは各ピーク間の時間Tz 、
T3、・・・が予め設定した最大待ち時間Tよりも大き
い場合(すなわちピークが検出されない場合)、後述す
るフローチャートに示すように別途aとbの値を規定す
る。
yl, y2, and y3 each indicate the y value of the first to third peaks, and have positive and negative signs. yo is a positive value of y that provides the conditions for starting the auto-tuning process, and y
When the absolute value of exceeds yo, auto-tuning processing is started. However, the time TI from the start of auto-tuning to the first peak, or the time Tz between each peak,
When T3, . . . is larger than the preset maximum waiting time T (that is, when no peak is detected), the values of a and b are separately defined as shown in the flowchart described later.

次に、第3図(A)、第3図(B)に示すフロー図によ
り、実施例オートチューニング調節計の動作を説明する
Next, the operation of the embodiment auto-tuning controller will be explained with reference to flowcharts shown in FIGS. 3(A) and 3(B).

先ず、予め設定したyoに対し、lyl>y。First, for a preset yo, lyl>y.

か否か判定する〔ステップST (STという〕■〕。It is determined whether or not [step ST (referred to as ST]■)].

yl>yoとなると制御状態が良くないとみなし、処理
を開始する。PID制御は継続しながら応答を計測し、
ピークの値を求める。つまり、T時間内にピークありか
否か判定する(ST2)、一定時間内にピークが検出で
きない場合はa、 bを規定値とする(ST3)。ST
2で、T時間内にピーク値があると、このピーク値y1
を第1ピーク値として記憶する(ST4)。そして、さ
らに以後T時間内にピークありか否か判定しく5T5)
、ピーク値がなければ、a、bを規定値としく5T6)
、ピーク値があれば、このピーク値y2を第2ピーク値
として記憶する。
When yl>yo, it is assumed that the control condition is not good and processing is started. Measure the response while continuing PID control,
Find the peak value. That is, it is determined whether or not there is a peak within time T (ST2). If a peak cannot be detected within a certain time, a and b are set as specified values (ST3). ST
2, if there is a peak value within time T, this peak value y1
is stored as the first peak value (ST4). Then, it is further determined whether there is a peak within T time 5T5)
, if there is no peak value, set a and b to the specified values 5T6)
, if there is a peak value, this peak value y2 is stored as the second peak value.

)’+    ’lz 制御開始時は、この判定がNOであり、次にT時間内に
ピークありか否か判定しく5T9)、N。
)'+'lz At the start of control, this determination is NO, and then it is determined whether or not there is a peak within time T (5T9), NO.

であれば、yの現在値をy、とし、ピークがあれば、ピ
ーク値y sを第3のピーク値として記憶する(STI
I)。また、ピーク値がなければ、現在値yを第3のピ
ーク値y、とする(STIO)。
If so, set the current value of y to y, and if there is a peak, store the peak value ys as the third peak value (STI
I). Furthermore, if there is no peak value, the current value y is set as the third peak value y (STIO).

そして、これらピークイ直(’/I% Yz 、)’s
 )からa、bを算出しく5T12)、a、bの値の範
囲に応じて更新係数に、 、kzを求める(STI3〜
5T19)そして、kl=kz=1か否か判定しく5T
20)、Noであれば、制御パラメータP(比例ゲイン
)、TI(積分定数)にそれぞれに、 、k、を乗じた
値を新たなパラメータ値としてP、TIを更新する(S
T21)。以下、応答計測からパラメータ更新までの処
理を繰り返し実行する。そして更新係数に、 、k2を
算出の際、k、=に2 =1となった時点で繰り返しを
やめ、パラメータの更新を停止して、オートチューニン
グを完了する。完了後でもlyl>yOとなれば再びチ
ューニングを開始する。なお応答計測中に現在の値yに
対して(y−y2)/ (yl−y2)の値が1となっ
た時点で、ST8の判定がYESとなり、aが1以上に
なると判断し、ただちにオートチューニングを中止して
、制御パラメータを予め設定した安全値に変更する(S
T22)。
And these peaky direct('/I% Yz,)'s
), calculate a and b from
5T19) Then, it is necessary to judge whether kl=kz=1 or not.5T
20), if No, update P and TI by setting the values obtained by multiplying the control parameters P (proportional gain) and TI (integral constant) by , k as new parameter values (S
T21). Thereafter, the process from response measurement to parameter update is repeatedly executed. Then, when calculating the update coefficient, k2, when k,=2=1, the repetition is stopped, parameter updating is stopped, and auto-tuning is completed. Even after completion, if lyl>yO, tuning is started again. During response measurement, when the value of (y-y2)/(yl-y2) becomes 1 for the current value y, the judgment in ST8 becomes YES, and it is determined that a is greater than or equal to 1. Cancel auto-tuning and change the control parameters to preset safe values (S
T22).

これにより万一チューニングが良好に行われず、プラン
トの応答が発散傾向になった場合でも、危険が未然に回
避される。
As a result, even if tuning is not performed properly and the response of the plant tends to diverge, danger can be avoided.

なお、上記実施例ではPID制御としているが、これを
PI制御で置き替えてもよい。
Note that although PID control is used in the above embodiment, this may be replaced with PI control.

また、上記実施例にある最大待ち時間Tは単に応答計測
の時間的打ち切り限界を定めるものであり、ピーク間の
時間すなわち振動周期自体をチューニングの必要条件と
するわけではない。
Furthermore, the maximum waiting time T in the above embodiment merely defines the time limit of response measurement, and does not make the time between peaks, that is, the vibration period itself, a necessary condition for tuning.

また、a、bの場合分けの範囲を示す不等式およびに+
、kzを算出する式は適用されるシステムにより上記実
施例とは異なるものであうでもよい。
In addition, there are also inequalities that indicate the range of cases for a and b, and +
, kz may be different from those in the above embodiment depending on the system to which it is applied.

また、振幅減衰率1以上を予測した時点でオートチュー
ニングを中止してパラメータを安全値にするのは、万一
のための安全機能であり、チューニングの必要条件では
ない。
Furthermore, stopping auto-tuning and setting parameters to safe values when an amplitude attenuation rate of 1 or more is predicted is a safety function in case of an emergency, and is not a necessary condition for tuning.

(へ)発明の効果 この発明によれば、従来の方式に比べ、プラントに余分
な外乱を与えることなく通常の制御状態で安定してチュ
ーニングができる。
(F) Effects of the Invention According to the present invention, tuning can be performed stably in a normal control state without causing unnecessary disturbance to the plant, compared to conventional methods.

実施例では一定時間以上ピークを検出しないと応答計測
を打ち切り、更新係数の算出に移るので、必要以上の長
時間を要しない。また応答計測中にプラントが発散傾向
になると、直ちに制御パラメータを安全値に変更するの
で、安全生が確保される。
In the embodiment, if a peak is not detected for a certain period of time or more, the response measurement is aborted and the process moves on to calculating the update coefficient, so that an unnecessarily long time is not required. Furthermore, if the plant tends to diverge during response measurement, the control parameters are immediately changed to safe values, so safety is ensured.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は、この発明の一実施例を示すオートチューニン
グ調節計の構成を示すブロック図、第2図は、プラント
の閉ループ応答を示す図、第3図(A)、第3図(B)
は、実施例オートチューニング調節計の動作を説明する
ためのフロー図である。 2:比例演算部、  3:積分演算部、4:微分演算部
、 5:オートチューニングアルゴリズム部。
FIG. 1 is a block diagram showing the configuration of an autotuning controller according to an embodiment of the present invention, FIG. 2 is a diagram showing the closed loop response of the plant, FIG. 3 (A), FIG. 3 (B)
FIG. 2 is a flowchart for explaining the operation of the auto-tuning controller according to the embodiment. 2: Proportional calculation section, 3: Integral calculation section, 4: Differential calculation section, 5: Auto-tuning algorithm section.

Claims (1)

【特許請求の範囲】[Claims] (1)PIまたはPID制御を行う調節計において、 被制御出力を計測する手段と、被制御出力の振幅減衰率
を算出する手段と、行き過ぎ量を算出する手段と、算出
された振幅減衰率と行き過ぎ量とから制御パラメータを
更新するための更新係数を算出する手段と、算出された
更新係数でパラメータを更新する手段を持ち、PIまた
はPID制御を実行しながら、更新係数が所定値となる
まで振幅減衰率算出、行き過ぎ量算出、更新係数算出、
パラメータ更新を繰り返す調整手段を持つことを特徴と
するオートチューニング調節計。
(1) In a controller that performs PI or PID control, a means for measuring the controlled output, a means for calculating the amplitude attenuation rate of the controlled output, a means for calculating the overshoot amount, and the calculated amplitude attenuation rate. It has a means for calculating an update coefficient for updating the control parameter from the overshoot amount and a means for updating the parameter with the calculated update coefficient, and while executing PI or PID control, until the update coefficient reaches a predetermined value. Amplitude attenuation rate calculation, excess amount calculation, update coefficient calculation,
An auto-tuning controller characterized by having an adjustment means that repeatedly updates parameters.
JP2198389A 1990-07-26 1990-07-26 Auto tuning controller Expired - Lifetime JP2865109B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2198389A JP2865109B2 (en) 1990-07-26 1990-07-26 Auto tuning controller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2198389A JP2865109B2 (en) 1990-07-26 1990-07-26 Auto tuning controller

Publications (2)

Publication Number Publication Date
JPH0484201A true JPH0484201A (en) 1992-03-17
JP2865109B2 JP2865109B2 (en) 1999-03-08

Family

ID=16390316

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2198389A Expired - Lifetime JP2865109B2 (en) 1990-07-26 1990-07-26 Auto tuning controller

Country Status (1)

Country Link
JP (1) JP2865109B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5691896A (en) * 1995-08-15 1997-11-25 Rosemount, Inc. Field based process control system with auto-tuning
US5812428A (en) * 1995-09-22 1998-09-22 Rosemount Inc. Process controller having non-integrating control function and adaptive bias
US5818714A (en) * 1996-08-01 1998-10-06 Rosemount, Inc. Process control system with asymptotic auto-tuning
JP2007213483A (en) * 2006-02-13 2007-08-23 Idemitsu Kosan Co Ltd Optimal adjustment system and optimum adjustment method for PID controller

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4003733B2 (en) * 2003-10-22 2007-11-07 株式会社日立製作所 PLANT SYSTEM, CONTROL DEVICE AND CONTROL METHOD

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01205202A (en) * 1988-02-12 1989-08-17 Hitachi Ltd auto tuning controller

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01205202A (en) * 1988-02-12 1989-08-17 Hitachi Ltd auto tuning controller

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5691896A (en) * 1995-08-15 1997-11-25 Rosemount, Inc. Field based process control system with auto-tuning
US5812428A (en) * 1995-09-22 1998-09-22 Rosemount Inc. Process controller having non-integrating control function and adaptive bias
US5818714A (en) * 1996-08-01 1998-10-06 Rosemount, Inc. Process control system with asymptotic auto-tuning
JP2007213483A (en) * 2006-02-13 2007-08-23 Idemitsu Kosan Co Ltd Optimal adjustment system and optimum adjustment method for PID controller

Also Published As

Publication number Publication date
JP2865109B2 (en) 1999-03-08

Similar Documents

Publication Publication Date Title
JPS62502433A (en) How to adjust digital regulator settings
JP2954660B2 (en) Model predictive controller
KR102451225B1 (en) Adjustment support device of the control unit and the adjustment device of the control unit
JP2865109B2 (en) Auto tuning controller
US5796608A (en) Self controllable regulator device
JP6511483B2 (en) Control device with error compensation with learning ability
JP2019530103A (en) Control device with adjustable control behavior
JP3470197B2 (en) Self-adjusting controller
JP2641855B2 (en) Adaptive control device
JP2651152B2 (en) PID control device
JPH04138501A (en) auto tuning controller
JPH0934503A (en) Adjustment method for pid controller
JPH11161301A (en) PID controller with auto tuning function
JP3774376B2 (en) Method and apparatus for identifying limit gain and transfer function of control system
JPH06208403A (en) Control device
JPH01276202A (en) auto tuning controller
JP2003189653A (en) Motor control device and gain setting method thereof
RU2393520C1 (en) Method to stabilise oscillating object output signal
JPH043202A (en) PID control method
JP2839679B2 (en) Automatic tuning device for control parameters
JP3779512B2 (en) Controller design equipment
JPH02201605A (en) Adjustor
JPS61226803A (en) Process control device
JPS641801B2 (en)
JPH0395602A (en) Control method for controller

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071218

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20081218

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20091218

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20091218

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20101218

Year of fee payment: 12

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101218

Year of fee payment: 12