JPH0885083A - Method for deriving and calibrating tool parameter of robot - Google Patents

Method for deriving and calibrating tool parameter of robot

Info

Publication number
JPH0885083A
JPH0885083A JP22155894A JP22155894A JPH0885083A JP H0885083 A JPH0885083 A JP H0885083A JP 22155894 A JP22155894 A JP 22155894A JP 22155894 A JP22155894 A JP 22155894A JP H0885083 A JPH0885083 A JP H0885083A
Authority
JP
Japan
Prior art keywords
tool
robot
parameter
coordinate system
deriving
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
JP22155894A
Other languages
Japanese (ja)
Other versions
JP2774939B2 (en
Inventor
Masatoshi Hida
正俊 飛田
Rikiya Senda
力哉 先田
Toshio Kamigaki
敏雄 神垣
Toshihiro Fujii
稔洋 藤井
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.)
Kobe Steel Ltd
Original Assignee
Kobe Steel Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kobe Steel Ltd filed Critical Kobe Steel Ltd
Priority to JP6221558A priority Critical patent/JP2774939B2/en
Publication of JPH0885083A publication Critical patent/JPH0885083A/en
Application granted granted Critical
Publication of JP2774939B2 publication Critical patent/JP2774939B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Manipulator (AREA)

Abstract

PURPOSE: To provide a method for deriving and calibrating the.tool parameter of a robot, by which the tool parameter of the robot can be more easily derived and calibrated at high precision in a short time. CONSTITUTION: A coordinate system of a flange 4 as an attaching part for a tool 3 is calculated on the basis of respective positioning data at the time of positioning of the tool 3 on the pin tip 5a with at least three attitudes, and at least two straight lines including the tool, position candidacy seen from optional two points on the flange coordinate system to be indicated by using a tool parameter are found, and the tool parameters is derived on the basis of the intersection of mutual straight lines and the flange coordinate system. The calibration value of the tool parameter may be estimated on the basis of the deviation in relation to the mean value of the tool position candidacy seen from the flange coordinate system to be indicated by using the tool parameter.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は,ロボットのツールパラ
メータ導出方法及び較正方法に係り,例えばロボットの
ツールを交換したときや,ツールが作業ワーク等に衝突
したときなどツールパラメータの変更が生じたときに,
ツールの先端位置及び方向を決定するツールパラメータ
を導出又は較正する方法に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a tool parameter deriving method and a calibrating method for a robot. For example, when the tool of the robot is exchanged or when the tool collides with a work or the like, the tool parameter is changed. sometimes,
The present invention relates to a method for deriving or calibrating tool parameters that determine the tip position and orientation of a tool.

【0002】[0002]

【従来の技術】従来,ロボットのツールパラメータを導
出又は較正する方法として,例えば特開昭64−787
74号公報が開示されている。ここでは,ツール先端点
を任意の1点に4回以上位置決めし,この時,ツール取
り付け部であるフランジフレームが球状に配置されるこ
とからツール先端点を求める方法である。別の方法とし
て,特開昭61−25206号公報及び特開昭61−2
5207号公報には,較正用治具を用いた方法が開示さ
れている。
2. Description of the Related Art Conventionally, as a method for deriving or calibrating a tool parameter of a robot, for example, Japanese Patent Laid-Open No. 64-787 has been used.
Japanese Patent Publication No. 74 is disclosed. In this method, the tool tip point is positioned at any one point four times or more, and at this time, the tool tip point is obtained because the flange frame, which is the tool mounting portion, is arranged in a spherical shape. As another method, Japanese Patent Laid-Open Nos. 61-25206 and 61-2
Japanese Patent No. 5207 discloses a method using a calibration jig.

【0003】[0003]

【発明が解決しようとする課題】上記したような従来の
ロボットのツールパラメータ導出方法及び較正方法で
は,次のような問題点があった。特開昭64−7877
4号公報に開示の方法では,作業者がロボットの先端に
取り付けたツール最先端が較正位置と一致するように位
置決めを行うが,この作業は作業者の技量に依存する。
即ち,作業者の熟練度等がその位置決め精度や,さらに
は,ツール最先端の導出精度に影響を及ぼす。また,こ
の方法ではツール最先端の位置のみしか導出できず,ツ
ールの方向は設計値を用いるかあるいは実測する必要が
あった。また,特開昭61−25206号公報,特開昭
61−25207号公報に開示の方法では,測定点数は
少なくてすむものの,較正治具の製作精度や,治具の据
え付けの為の段取り替え時間が問題となっていた。本発
明は,上記事情に鑑みてなされたものであり,その第1
の目的とするところは,より簡便に短時間にしかも精度
よくロボットのツールパラメータを導出することであ
る。また,第2の目的は,同様にロボットのツールパラ
メータを較正することである。
The conventional robot tool parameter deriving method and calibrating method described above have the following problems. JP-A-64-7877
In the method disclosed in Japanese Patent No. 4, the operator positions the tool so that the tip of the tool attached to the tip of the robot coincides with the calibration position, but this operation depends on the skill of the operator.
That is, the skill level of the operator affects the positioning accuracy and the leading-edge derivation accuracy of the tool. In addition, this method can only derive the position of the tip of the tool, and it was necessary to use the design value or actually measure the tool direction. In the methods disclosed in Japanese Patent Laid-Open Nos. 61-25206 and 61-25207, the number of measurement points is small, but the accuracy of calibration jig preparation and the setup change for jig installation. Time was a problem. The present invention has been made in view of the above circumstances.
The purpose of is to derive the tool parameters of the robot more simply and in a short time and with high accuracy. The second purpose is to calibrate the robot tool parameters as well.

【0004】[0004]

【課題を解決するための手段】上記第1の目的を達成す
るために第1の発明は,ロボットのアーム先端に取り付
けられたツールの位置姿勢を決定するツールパラメータ
を導出する方法において,上記ツールを空間上の1点に
少なくとも異なる3つの姿勢で位置決めしたときの各位
置決めデータに基づいて上記ツールの取り付け部の座標
系を演算し,ツールパラメータを用いて表現される上記
取り付け部の座標系の任意の2つよりみたツール位置候
補を含む直線を少なくとも2本求め,上記直線同士の交
点と取り付け部の座標系とに基づいて上記ツールパラメ
ータを導出してなることを特徴とするロボットのツール
パラメータ導出方法として構成されている。さらには,
上記ツールの位置決めにおける各姿勢のうち,少なくと
も1つは既知の方向とすることを特徴とするロボットの
ツールパラメータ導出方法である。また第2の目的を達
成するために第2の発明は,ロボットのアーム先端に取
り付けられたツールの位置姿勢を決定するツールパラメ
ータを較正する方法において,上記ツールを空間上の1
点に少なくとも異なる3つの姿勢で位置決めしたときの
各位置決めデータに基づいて上記ツールの取り付け部の
座標系を演算し,ツールパラメータを用いて表現される
上記取り付け部の座標系よりみたツール位置候補の平均
値に対する偏差に基づいて上記ツールパラメータの較正
値を推定してなることを特徴とするロボットのツールパ
ラメータ較正方法として構成されている。さらには,上
記ツールの位置決めにおける各姿勢のうち,少なくとも
1つは既知の方向とすることを特徴とするロボットのツ
ールパラメータ較正方法である。
In order to achieve the above first object, a first invention is a method for deriving a tool parameter for determining the position and orientation of a tool attached to the tip of an arm of a robot. The coordinate system of the mounting part of the tool is calculated based on the respective positioning data when positioning is performed on one point in space in at least three different postures, and the coordinate system of the mounting part of the mounting part expressed by using tool parameters is calculated. A tool parameter of a robot, characterized in that at least two straight lines including tool position candidates viewed from arbitrary two are obtained, and the tool parameters are derived based on the intersection of the straight lines and the coordinate system of the mounting portion. It is configured as a derivation method. Furthermore,
At least one of the postures for positioning the tool is a known direction, which is a method for deriving the tool parameter of the robot. In order to achieve the second object, a second invention is a method of calibrating tool parameters for determining the position and orientation of a tool attached to the tip of an arm of a robot, wherein
The coordinate system of the attachment part of the tool is calculated based on each positioning data when the point is positioned in at least three different postures, and the tool position candidate viewed from the coordinate system of the attachment part expressed using the tool parameter is selected. The tool parameter calibration method for a robot is characterized in that the calibration value of the tool parameter is estimated based on the deviation from the average value. Furthermore, in the tool parameter calibration method for a robot, at least one of the postures for positioning the tool is set to a known direction.

【0005】[0005]

【作用】第1の発明によれば,ロボットのアーム先端に
取り付けられたツールの位置姿勢を決定するツールパラ
メータを導出するに際し,上記ツールを空間上の1点に
少なくとも異なる3つの姿勢で位置決めしたときの各位
置決めデータに基づいて上記ツールの取り付け部の座標
系が演算される。ツールパラメータを用いて表現される
上記取り付け部の座標系の任意の2つより見たツール位
置候補を含む直線が少なくとも2本求められる。上記直
線同士の交点と取り付け部の座標系とに基づいて,上記
ツールパラメータが導出される。このように,ロボット
の姿勢を変化させる教示点を極力少なくし,且つ,専用
の較正治具を用いなくても,より簡便に短時間に精度よ
くロボットのツールパラメータを導出することができ
る。さらに,上記ツールの位置決めにおける各姿勢のう
ち,少なくとも1つを既知の方向とすれば,ツールパラ
メータをその方向成分を含めて導出することができる。
また,第2の発明によれば,ロボットのアーム先端に取
り付けられたツールの位置姿勢を決定するツールパラメ
ータを較正するに際し,上記ツールを空間上の1点に少
なくとも異なる3つの姿勢で位置決めしたときの各位置
決めデータに基づいて,上記ツールの取り付け部の座標
系が演算される。ツールパラメータを用いて表現される
上記取り付け部の座標系より見たツール位置候補の平均
値に対する偏差に基づいて,上記ツールパラメータの較
正値が推定される。このように,ロボットの姿勢を変化
させる教示点を極力少なくし,かつ,専用の較正治具を
用いなくとも,より簡便にに短時間に精度よくロボット
のツールパラメータを較正することができる。さらに,
上記ツールの位置決めにおける各姿勢のうち,少なくと
も1つを既知の方向とすれば,ツールパラメータをその
方向成分を含めて較正することができる。
According to the first aspect of the invention, when deriving the tool parameters for determining the position and orientation of the tool attached to the tip of the robot arm, the tool is positioned at one point in space in at least three different orientations. The coordinate system of the attachment part of the tool is calculated based on each positioning data at this time. At least two straight lines including tool position candidates viewed from any two of the coordinate systems of the mounting portion expressed using the tool parameter are obtained. The tool parameter is derived based on the intersection of the straight lines and the coordinate system of the mounting portion. In this way, the number of teaching points that change the posture of the robot can be minimized, and the tool parameters of the robot can be derived easily and accurately in a short time without using a dedicated calibration jig. Furthermore, if at least one of the postures for positioning the tool is a known direction, the tool parameter can be derived including the direction component.
According to the second invention, when the tool parameters for determining the position and orientation of the tool attached to the arm tip of the robot are calibrated, the tool is positioned at one point in space in at least three different orientations. The coordinate system of the attachment part of the tool is calculated based on each positioning data of. The calibration value of the tool parameter is estimated on the basis of the deviation from the average value of the tool position candidates viewed from the coordinate system of the mounting portion expressed using the tool parameter. In this way, the number of teaching points that change the posture of the robot can be minimized, and the tool parameters of the robot can be calibrated more easily and accurately in a short time without using a dedicated calibration jig. further,
If at least one of the postures for positioning the tool is a known direction, the tool parameter can be calibrated including the direction component.

【0006】[0006]

【実施例】以下添付図面を参照して,本発明を具体化し
た実施例につき説明し,本発明の理解に供する。尚,以
下の実施例は,本発明を具体化した一例であって,本発
明の技術的範囲を限定する性格のものではない。ここ
に,図1は第1の発明の一実施例(第1の実施例)に係
るロボットのツールパラメータの導出方法の概略手順を
示すフロー図,図2は第2の発明の一実施例(第2の実
施例)に係るロボットのツールパラメータの較正方法の
概略手順を示すフロー図,図3は上記第1,第2の実施
例に係る方法を適用可能な装置0の概略構成図,図4は
ツール座標系の設定の様子を示す説明図,図5はツール
パラメータ導出の原理を示す説明図,図6はツールパラ
メータ較正の原理を示す説明図である。図1に示すごと
く,第1の発明の一実施例(第1の実施例)に係るロボ
ットのツールパラメータ導出方法は,ロボットのアーム
先端に取り付けられたツールの位置姿勢を決定するツー
ルパラメータを導出するに際し,上記ツールを空間上の
1点に少なくとも異なる3の姿勢で位置決めしたとき
の,各位置決めデータに基づいて,上記ツールの取り付
け部であるフランジの座標系を演算し,ツールパラメー
タを用いて表現される上記フランジ座標系の任意の2つ
より見たツール位置候補を含む直線を少なくとも2本求
め,上記直線同士の交点とフランジ座標系とに基づいて
上記ツールパラメータを導出するように構成されてい
る。さらに,上記ツールの位置決めにおける各姿勢の
内,少なくとも1つは既知の方向とするようにしてもよ
く,図1では第1回目の位置決めはツールの方向を既知
の方向ベクトルgに合わせて行うものとしている。図3
は上記第1の実施例に係るツールパラメータ導出方法を
適用可能な装置0の概略構成を示す。図中,1はロボッ
ト本体,2はアーム,3はツール,4はフランジ,5は
ピン(5aはその先端),6はコントローラ,7はティ
ーチングペンダント,8はパーソナルコンピュータを示
す。又,パーソナルコンピュータ8は例えばCPU,C
RT,メモリ,キーボード等を備えている市販のものを
用いる。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Embodiments embodying the present invention will be described below with reference to the accompanying drawings for the understanding of the present invention. The following embodiments are examples of embodying the present invention and are not intended to limit the technical scope of the present invention. 1 is a flow chart showing a schematic procedure of a method of deriving a tool parameter of a robot according to an embodiment (first embodiment) of the first invention, and FIG. 2 is an embodiment of the second invention ( Second Embodiment) A flow chart showing a schematic procedure of a method for calibrating a tool parameter of a robot according to the second embodiment, and FIG. 3 is a schematic configuration diagram of an apparatus 0 to which the method according to the first and second embodiments can be applied. 4 is an explanatory view showing a state of setting a tool coordinate system, FIG. 5 is an explanatory view showing a principle of deriving tool parameters, and FIG. 6 is an explanatory view showing a principle of tool parameter calibration. As shown in FIG. 1, a method of deriving a tool parameter of a robot according to an embodiment (first embodiment) of the first invention derives a tool parameter that determines the position and orientation of a tool attached to the arm tip of a robot. In doing so, the coordinate system of the flange that is the mounting part of the tool is calculated based on each positioning data when the tool is positioned at one point in space in at least three different postures, and the tool parameters are used. At least two straight lines including tool position candidates viewed from any two of the expressed flange coordinate systems are obtained, and the tool parameters are derived based on the intersection of the straight lines and the flange coordinate system. ing. Further, at least one of the postures for positioning the tool may be set to a known direction. In FIG. 1, the first positioning is performed by aligning the tool direction with a known direction vector g. I am trying. FIG.
Shows a schematic configuration of an apparatus 0 to which the tool parameter derivation method according to the first embodiment can be applied. In the figure, 1 is a robot body, 2 is an arm, 3 is a tool, 4 is a flange, 5 is a pin (5a is its tip), 6 is a controller, 7 is a teaching pendant, and 8 is a personal computer. The personal computer 8 is, for example, a CPU, C
A commercially available product equipped with an RT, a memory, a keyboard, etc. is used.

【0007】以下,この装置0を用いて,上記第1の実
施例方法をさらに具現化する。図4はツール座標系の設
定の様子を示している。図中,フランジ4の中心にはフ
ランジ座標系が設定されている。フランジ座標系はロボ
ット1の各軸のデータよりコントローラ6にて計算され
る。ここでは,空間上の1点は位置設定が容易に行える
ようにピン5の先端5aとしている。オペレータは,ま
ずティーチングペンダント7によりある姿勢で位置決め
する。このときのピン先端5aの位置をロボット座標上
の位置でp(px ,py ,pz ),フランジの座標系を
1 (l 1x,l1y,l1z,m1x,m1y,m1z,n1x,n
1y,n1z,o1x,o1y,o1z),求めたいフランジ上の
ツール位置ベクトルをt(tx ,ty ,tz )とする
と,これらの関係を次式で表すことができる。この時の
関係を図5(a)に示した。
Hereinafter, by using this device 0, the first actual
The embodiment method is further embodied. Figure 4 shows the tool coordinate system setup.
It shows the fixed situation. In the figure, the flange 4 is located at the center of the flange.
The lunge coordinate system is set. The flange coordinate system is Robo
Is calculated by the controller 6 from the data of each axis of
It Here, the position of one point in space can be easily set.
Thus, the tip 5a of the pin 5 is used. The operator
Without teaching pendant 7, positioning in a certain posture
To do. The position of the pin tip 5a at this time is on the robot coordinates.
At the position p (px, Py, Pz), The coordinate system of the flange
F1(L 1x, L1y, L1z, M1x, M1y, M1z, N1x, N
1y, N1z, O1x, O1y, O1z), On the desired flange
Let the tool position vector be t (tx, Ty, Tz) And
And the relationship between them can be expressed by the following equation. At this time
The relationship is shown in FIG.

【数1】 次にオペレータは,ロボットの姿勢を変化させる。この
ときのフランジ座標系をF2 (l2x,l2y,l2z
2x,m2y,m2z,n2x,n2y,n2z,o2x,o2y,o
2z)およびF3 (l3x,l3y,l3z,m3x,m3y
3z,n3x,n3y,n3z,o3x,o3y,o3z)とする
と,次の2つの式が成立する。
[Equation 1] Next, the operator changes the posture of the robot. The flange coordinate system at this time is F 2 (l 2x , l 2y , l 2z ,
m 2x , m 2y , m 2z , n 2x , n 2y , n 2z , o 2x , o 2y , o
2z ) and F 3 (l 3x , l 3y , l 3z , m 3x , m 3y ,
m 3z , n 3x , n 3y , n 3z , o 3x , o 3y , o 3z ), the following two equations hold.

【数2】 [Equation 2]

【0008】今,これらの3個の座標系のうち,i番目
とj番目(i,j=1〜3;i≠j)の座標系を考え
る。上記(1)式〜(3)式は,いずれもロボット座標
からのピン先端5aの位置ベクトルp(px ,py ,p
z )を示していることから次式が成立する。
Now, consider the i-th and j-th (i, j = 1 to 3; i ≠ j) coordinate systems among these three coordinate systems. The above equations (1) to (3) are all the position vectors p (p x , py , p of the pin tip 5a from the robot coordinates.
Since z ) is shown, the following equation holds.

【数3】 一般的に原点が同一で姿勢の異なる座標系の間には,図
5(b)に示すようにあるベクトル方向κにθ回転する
座標変換行列によって関係づけることができる。ここ
で,原点が同じで,姿勢が異なる2つの座標系をそれぞ
れ(lix,liy,liz,mix,miy,miz,nix
iy,niz,0,0,0)および(ljx,ljy,ljz
jx,mjy,mjz,njx,njy,njz,0,0,0)と
する。この間の関係を表現する座標変換行列をCとする
と,次式が成立する。
(Equation 3) Generally, coordinate systems having the same origin and different attitudes can be related by a coordinate transformation matrix that rotates θ in a certain vector direction κ as shown in FIG. 5B. Here, the origin is the same, orientation different two coordinate systems, respectively (l ix, l iy, l iz, m ix, m iy, m iz, n ix,
n iy , n iz , 0, 0, 0) and (l jx , l jy , l jz ,
m jx , m jy , m jz , n jx , n jy , n jz , 0, 0, 0). When the coordinate transformation matrix expressing the relationship between them is C, the following equation holds.

【数4】 上記(5)式より変換行列Cは次のようになる。[Equation 4] From the above equation (5), the transformation matrix C is as follows.

【数5】 (Equation 5)

【0009】一般に知られているように,変換行列Cは
等価回転軸κ及び等価回転角θにより変換を表す等価回
転行列Rot(κ,θ)に置き換えることができる。
As is generally known, the transformation matrix C can be replaced with an equivalent rotation matrix Rot (κ, θ) that represents transformation by an equivalent rotation axis κ and an equivalent rotation angle θ.

【数6】 上記(6)及び(7)式より,等価回転角θ及び等価回
転軸κは次のような値となる。
(Equation 6) From the above equations (6) and (7), the equivalent rotation angle θ and the equivalent rotation axis κ have the following values.

【数7】 この等価回転行列Rot(κ,θ)は方向ベクトルであ
る等価回転軸κに対して何の変換も行わない。従って,
αを任意の実数とすると,次式が成立する。 Rot(κ,θ)ακ=αCκ=ακ …(10)
(Equation 7) This equivalent rotation matrix Rot (κ, θ) does not perform any conversion on the equivalent rotation axis κ which is a direction vector. Therefore,
If α is an arbitrary real number, the following equation holds. Rot (κ, θ) ακ = αCκ = ακ (10)

【0010】上記(5)式の両辺にακをかけ,上記
(10)式を用いて変形すると次式が得られる。
By multiplying both sides of equation (5) by ακ and transforming using equation (10), the following equation is obtained.

【数8】 このときのv(vx ,vy ,vz )はロボット座標系上
からみた等価回転軸ベクトルを意味する。いま,ここで
上記(4)式に上記(10)式を加えると,次式が得ら
れる。
[Equation 8] At this time, v (v x , v y , v z ) means an equivalent rotation axis vector viewed from the robot coordinate system. Now, when the above equation (10) is added to the above equation (4), the following equation is obtained.

【数9】 [Equation 9]

【0011】これは,2つのフランジ座標系上の直線上
の点t+ακは一致し,ロボット座標系上では直線p+
αv上の1点に変換されることを意味する。この直線上
の1点pij(pijx ,pijy ,pijz )(ツール位置候
補に相当)を,piji とp ijj の2つのベクトルの
外積が等価回転軸ベクトルvの方向で,且つ2つのベク
トルが成す角度がθとなる点に設定する。すると直線L
ijは次式で表すことができる。
This is on a straight line on the two flange coordinate systems.
Point t + ακ coincides with each other, and a straight line p + on the robot coordinate system
It means that it is converted into one point on αv. On this straight line
1 point pij(Pijx, Pijy, Pijz) (Tool position
Equivalent to), pijoiAnd p ijojOf the two vectors
The outer product is in the direction of the equivalent rotation axis vector v, and two vectors
Set to the point where the angle formed by the tor is θ. Then straight line L
ijCan be expressed by the following equation.

【数10】 ここでとりうるフランジ座標系の組(i,j)は(1,
2),(1,3),(2,3)の3組存在し,上記(1
3)式より直線は3個求めることができる。この3本の
直線は物理的な意味から1点に交わるので,適当に直線
を2本選ぶことによりその交点が演算できる。この交点
は,ツール先端を位置決めした同一点pであるので,3
つのフランジ座標系の内の1つFi (i=1〜3)を用
いて以下のようにツール位置ベクトルtを求める。
[Equation 10] The set of flange coordinate systems (i, j) that can be taken here is (1,
There are 3 sets of 2), (1, 3) and (2, 3).
From equation (3), three straight lines can be obtained. Since these three straight lines intersect at one point from the physical sense, the intersection can be calculated by appropriately selecting two straight lines. Since this intersection is the same point p where the tool tip is positioned, 3
The tool position vector t is obtained as follows using one of the two flange coordinate systems F i (i = 1 to 3).

【数11】 [Equation 11]

【0012】以上でツール位置ベクトルtを導出でき
た。以上の演算は,パーソナルコンピュータ8のCPU
により実行され,CRTに実行結果が表示される。この
ように,ロボットの姿勢を変化させる教示点を極力少な
くし,且つ専用の較正治具を用いなくても,より簡便に
短時間に精度よくロボットのツールパラメータ(tx
y ,tz )を導出することができた。次にツールの方
向を導出する原理について述べる。方向成分を加えたツ
ールパラメータ(tx ,ty ,tz ,ε1 ,ε2
ε3 )によってフランジ座標系でのツールの座標系は次
のように表現される。ただし,ε1 ,ε2 ,ε3 はいず
れもオイラー角を示す。
As described above, the tool position vector t can be derived. The above calculation is performed by the CPU of the personal computer 8.
And the result is displayed on the CRT. As described above, the number of teaching points for changing the posture of the robot is minimized, and the tool parameter (t x ,
ty , tz ) could be derived. Next, the principle of deriving the direction of the tool will be described. Tool parameters (t x , t y , t z , ε 1 , ε 2 ,
The tool coordinate system in the flange coordinate system is expressed by ε 3 ) as follows. However, ε 1 , ε 2 , and ε 3 all represent Euler angles.

【数12】 ここでは,図 1に示すようにツールの方向座標軸wを既
知の方向ベクトルg(gx ,gy ,gz )と一致させる
こととする。このとき,ツール座標系T及びフランジ座
標系Fi 及び方向ベクトルgの関係は以下のようにな
る。
[Equation 12] Here, as shown in FIG. 1, the direction coordinate axis w of the tool is made to coincide with the known direction vector g (g x , g y , g z ). At this time, the relationship between the tool coordinate system T, the flange coordinate system F i, and the direction vector g is as follows.

【数13】 [Equation 13]

【0013】フランジ座標系Fi 及び方向ベクトルgは
いずれも既知であるので,結局上記(17)式を変形す
ると次式が得られる。
Since the flange coordinate system F i and the direction vector g are both known, the following equation can be obtained by modifying the above equation (17).

【数14】 これより,オイラー角ε1 ,ε2 は以下のように求ま
る。 ε1 =atan2(rx ,ry ) …(19) ε2 =atan2(rz ,±√(rx 2 +ry 2 )) …(20) 以上より,上記(11)式,(16)式及び(17)式
よりツールの先端及びツールの方向を決定することがで
きる。この決定演算もパーソナルコンピュータ8のCP
Uにより実行され,CRTに実行結果が表示される。以
上によりツールパラメータをオイラー角表現したときの
(tx ,ty ,tz ,ε1 ,ε2 ,ε3 )の内の5つの
成分tx ,ty ,tz ,ε1 ,ε2 が求まる。残りの
1つの成分であるε3 はパーソナルコンピュータ8のキ
ーボード入力により設定値もしくはデフォルトを設定す
る。ロボットのツールパラメータを求めるには,上記し
たような第1の実施例のごとく,ロボットの姿勢を代え
てツール先端点を空間上の所定の位置5aに位置付けて
姿勢を変えて3回位置決めを行えば原理的には求まる。
しかし,実際にはロボットの姿勢を変え,ツール先端点
を所定点5aに近づける際に,正確に位置決めが行われ
ない場合があり得る。これらの場合を想定し,第2の発
明では3点以上の点を位置決めし,これら3点以上のフ
ランジ座標系の値より最小自乗法を用いてツールパラメ
ータの較正を行うこととした。
[Equation 14]From this, the Euler angle ε1, Ε2Is calculated as follows
It ε1= Atan2 (rx, Ry)… (19) ε2= Atan2 (rz, ± √ (rx 2+ Ry 2)) (20) From the above, the above equations (11), (16) and (17)
It is possible to determine the tip of the tool and the direction of the tool.
Wear. This decision calculation is also the CP of the personal computer 8.
It is executed by U and the execution result is displayed on the CRT. Since
When the tool parameters are expressed by Euler angles by the above
(Tx, Ty, Tz , Ε1, Ε2, Ε3) Of 5
Ingredient tx, Ty, Tz, Ε1, Ε2Is required. Remaining
One component, ε3Is the key for the personal computer 8.
Set value or default by board input
It To find the tool parameters of the robot, do the above.
As in the first embodiment described above, the posture of the robot is changed.
Position the tip of the tool at a predetermined position 5a in the space.
In principle, this can be obtained by changing the posture and performing positioning three times.
However, in reality, changing the posture of the robot,
When positioning is close to the predetermined point 5a, accurate positioning is performed.
There may not be. Assuming these cases, the second departure
In Ming, 3 or more points are positioned, and 3 or more points are positioned.
From the values in the lunge coordinate system, the tool parameter is
It was decided to calibrate the data.

【0014】引き続いて第2の発明についてのべる。図
2に示すごとく,第2の発明の一実施例(第2の実施
例)に係るロボットのツールパラメータ較正方法は,ロ
ボットのアーム先端に取り付けられたツールの位置姿勢
を決定するツールパラメータを較正するに際し,上記ツ
ールを空間上の1点に少なくとも異なる3つの姿勢で位
置決めしたときの各位置決めデータに基づいて上記ツー
ルの取り付け部であるフランジの座標系を演算し,ツー
ルパラメータを用いて表現される上記フランジ座標系よ
り見たツール位置候補の平均値に対する偏差に基づい
て,上記ツールパラメータの較正値を推定するように構
成されている。さらに,上記ツールの位置決めにおける
各姿勢のうち少なくとも1つは既知の方向とするように
してもよく,図2では第1回目の位置決めツールの方向
を既知の方向ベクトルgに合わせて行うものとしてい
る。本第2の実施例に係るツールパラメータ較正方法に
ついても,上記第1の実施例におけると同様の装置0を
適用可能である。以下,この装置0による上記第2の実
施例方法を更に具現化する。本第2の実施例では,q回
(q≧3)の位置決めを行い,その時のそれぞれのフラ
ンジ座標系F1 〜Fq よりツール位置ベクトルtを求め
る。この時のフランジ座標系F1 〜Fq とツール位置ベ
クトルtとの関係を図6に示した。図6では,正確に位
置決めが行われていないため,各フランジ座標系から等
価回転軸κにおろした法線の位置は一致せず,誤差εij
が生じる。ここでは,まず先に示した(4)式を変形し
以下の方程式を得る。
Next, the second invention will be described. As shown in FIG. 2, a method for calibrating tool parameters of a robot according to an embodiment (second embodiment) of the second invention calibrates tool parameters for determining the position and orientation of a tool attached to the tip of an arm of a robot. In doing so, the coordinate system of the flange that is the attachment part of the tool is calculated based on each positioning data when the tool is positioned at one point in space in at least three different postures, and it is expressed using tool parameters. The calibration value of the tool parameter is estimated based on the deviation from the average value of the tool position candidates viewed from the flange coordinate system. Further, at least one of the postures for positioning the tool may be set to a known direction. In FIG. 2, the direction of the first positioning tool is set to match the known direction vector g. . The same apparatus 0 as in the first embodiment can be applied to the tool parameter calibration method according to the second embodiment. Hereinafter, the method of the second embodiment by the device 0 will be further embodied. In the second embodiment, the positioning is performed q times (q ≧ 3), and the tool position vector t is obtained from the respective flange coordinate systems F 1 to F q at that time. FIG. 6 shows the relationship between the flange coordinate systems F 1 to F q and the tool position vector t at this time. In FIG. 6, since the positioning is not performed accurately, the positions of the normals drawn from the respective flange coordinate systems to the equivalent rotation axis κ do not match, and the error ε ij
Occurs. Here, first, the above equation (4) is modified to obtain the following equation.

【数15】 (Equation 15)

【0015】この場合,6変数tx ,ty ,tz
x ,py ,pz に対する方程式が6個導かれる。しか
し,これらの方程式だけでは解が一意に求まらない。逆
に,先に第1の実施例で述べた原理に従えば,少なくと
もフランジ座標系を3つ用いて方程式を作れば,解が一
意に求まることが明らかである。全てのフランジ座標系
を用いると以下の方程式が得られる。
In this case, the six variables t x , t y , t z ,
Six equations for p x , p y , and p z are derived. However, a solution cannot be uniquely obtained from these equations alone. On the contrary, according to the principle described in the first embodiment, it is apparent that a solution can be uniquely obtained by forming an equation using at least three flange coordinate systems. The following equation is obtained using all flange coordinate systems.

【数16】 今,ツール位置ベクトルtが誤差Δt(Δtx ,Δ
y ,Δtz )を含んでいるとき,上記(22)式より
以下の誤差方程式が導かれる。
[Equation 16] Now, the tool position vector t has an error Δt (Δt x , Δ
When t y , Δt z ) is included, the following error equation is derived from the above equation (22).

【数17】 [Equation 17]

【0016】同一点の誤差Δpi を以下のように導出さ
れる同一点の平均からの誤差と定義する。これはツール
位置候補の平均値に対する偏差に相当し,前述した誤差
εijに比例する。
The error Δp i at the same point is defined as the error from the average of the same points derived as follows. This corresponds to the deviation of the tool position candidate from the average value and is proportional to the error ε ij described above.

【数18】 上記(23)式の左辺のフランジ座標系よりなる行列J
(F1 ,F2 ,…,F q T の擬似逆行列J+ を(2
3)式の両辺にかけると,以下の式が得られる。
[Equation 18]The matrix J consisting of the flange coordinate system on the left side of the above equation (23)
(F1, F2,,, F q)TPseudo-inverse matrix J of+To (2
By multiplying both sides of equation 3), the following equation is obtained.

【数19】 これより,ツール位置ベクトルtの誤差Δtが求まる。
この誤差Δtを用いてツール位置ベクトルtの値を以下
のように更新し,対応する位置の誤差Δpj の総和がΣ
|Δpj |<εとなるまで繰り返す。ただし,εは予め
定めた値である。 tnew =told +Δt …(26) また,ツールの方向は前述した既知の方向ベクトルgを
例えば(0,0,−1)とおくことにより前記第1の実
施例と同様にして求めることができる。以上のように,
本第2の実施例によれば,ロボットの姿勢を変化させる
教示点を極力少なくし,且つ専用の較正治具を用いなく
ても,より簡便に短時間に精度よくロボットのツールパ
ラメータを較正することができた。
[Formula 19] From this, the error Δt of the tool position vector t is obtained.
Using this error Δt, the value of the tool position vector t is updated as follows, and the sum of the error Δp j at the corresponding position is Σ
Repeat until | Δp j | <ε. However, ε is a predetermined value. t new = t old + Δt (26) Further, the direction of the tool can be obtained in the same manner as in the first embodiment by setting the known direction vector g described above to (0, 0, -1), for example. it can. As mentioned above,
According to the second embodiment, the number of teaching points for changing the posture of the robot is reduced as much as possible, and the tool parameters of the robot can be calibrated more simply and accurately in a short time without using a dedicated calibration jig. I was able to.

【0017】[0017]

【発明の効果】本発明に係るロボットのツールパラメー
タ導出方法及び較正方法は,上記したように構成されて
いるため,ロボットの位置姿勢を変化させる教示点を極
力少なくし,且つ,専用の較正治具を用いなくても,よ
り簡便に短時間に精度よくロボットのツールパラメータ
を導出及び較正することができる。
Since the robot tool parameter deriving method and calibration method according to the present invention are configured as described above, the number of teaching points for changing the position and orientation of the robot is minimized, and a dedicated calibration jig is used. Without using tools, the robot tool parameters can be derived and calibrated easily and in a short time with high accuracy.

【図面の簡単な説明】[Brief description of drawings]

【図1】 第1の発明の一実施例(第1の実施例)に係
るロボットのツールパラメータの導出方法の概略手順を
示すフロー図。
FIG. 1 is a flowchart showing a schematic procedure of a method of deriving a tool parameter of a robot according to an embodiment (first embodiment) of the first invention.

【図2】 第2の発明の一実施例(第2の実施例)に係
るロボットのツールパラメータの較正方法の概略手順を
示すフロー図。
FIG. 2 is a flowchart showing a schematic procedure of a method for calibrating a tool parameter of a robot according to an embodiment (second embodiment) of the second invention.

【図3】 上記第1,第2の実施例に係る方法を適用可
能な装置0の概略構成図。
FIG. 3 is a schematic configuration diagram of an apparatus 0 to which the methods according to the first and second embodiments can be applied.

【図4】 ツール座標系の設定の様子を示す説明図。FIG. 4 is an explanatory diagram showing how the tool coordinate system is set.

【図5】 ツールパラメータ導出の原理を示す説明図。FIG. 5 is an explanatory diagram showing the principle of tool parameter derivation.

【図6】 ツールパラメータ較正の原理を示す説明図。FIG. 6 is an explanatory diagram showing the principle of tool parameter calibration.

【符号の説明】[Explanation of symbols]

0…ロボットのツールパラメータ導出及び較正装置 1…ロボット本体 2…アーム 3…ツール 4…フランジ(取り付け部に相当) 5a…ピン先端(空間上の1点に相当) 0 ... Robot tool parameter derivation and calibration device 1 ... Robot body 2 ... Arm 3 ... Tool 4 ... Flange (corresponding to mounting part) 5a ... Pin tip (corresponding to one point in space)

───────────────────────────────────────────────────── フロントページの続き (72)発明者 神垣 敏雄 愛知県豊橋市三弥町字中原1番地2 株式 会社神戸製鋼所豊橋FA・ロボットセンタ ー内 (72)発明者 藤井 稔洋 愛知県豊橋市三弥町字中原1番地2 株式 会社神戸製鋼所豊橋FA・ロボットセンタ ー内 ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Toshio Kangaki, Nakahara, Toyohashi City, Aichi Prefecture 2 1 Nakahara, Kobe Steel Co., Ltd. Toyohashi FA / Robot Center (72) Minoru Fujii Inami, Miyachi Town Nakahara No. 1 Address 2 Kobe Steel, Ltd. Toyohashi FA / Robot Center

Claims (4)

【特許請求の範囲】[Claims] 【請求項1】 ロボットのアーム先端に取り付けられた
ツールの位置姿勢を決定するツールパラメータを導出す
る方法において,上記ツールを空間上の1点に少なくと
も異なる3つの姿勢で位置決めしたときの各位置決めデ
ータに基づいて上記ツールの取り付け部の座標系を演算
し,ツールパラメータを用いて表現される上記取り付け
部の座標系の任意の2つよりみたツール位置候補を含む
直線を少なくとも2本求め,上記直線同士の交点と取り
付け部の座標系とに基づいて上記ツールパラメータを導
出してなることを特徴とするロボットのツールパラメー
タ導出方法。
1. A method of deriving a tool parameter for determining the position and orientation of a tool attached to the end of a robot arm, wherein each positioning data is obtained when the tool is positioned at one point in space in at least three different postures. Based on the above, the coordinate system of the mounting portion of the tool is calculated, and at least two straight lines including tool position candidates viewed from any two of the coordinate systems of the mounting portion expressed using the tool parameters are obtained, and the straight line is obtained. A method for deriving a tool parameter of a robot, characterized in that the above-mentioned tool parameter is derived based on an intersection of the two and a coordinate system of a mounting portion.
【請求項2】 上記ツールの位置決めにおける各姿勢の
うち,少なくとも1つは既知の方向とすることを特徴と
する請求項1記載のロボットのツールパラメータ導出方
法。
2. The method for deriving a tool parameter of a robot according to claim 1, wherein at least one of the postures for positioning the tool has a known direction.
【請求項3】 ロボットのアーム先端に取り付けられた
ツールの位置姿勢を決定するツールパラメータを較正す
る方法において,上記ツールを空間上の1点に少なくと
も異なる3つの姿勢で位置決めしたときの各位置決めデ
ータに基づいて上記ツールの取り付け部の座標系を演算
し,ツールパラメータを用いて表現される上記取り付け
部の座標系よりみたツール位置候補の平均値に対する偏
差に基づいて上記ツールパラメータの較正値を推定して
なることを特徴とするロボットのツールパラメータ較正
方法。
3. A method for calibrating tool parameters for determining the position and orientation of a tool attached to the tip of a robot arm, wherein each positioning data is obtained when the tool is positioned at one point in space in at least three different postures. The coordinate system of the mounting part of the tool is calculated based on the above, and the calibration value of the tool parameter is estimated based on the deviation from the average value of the tool position candidates viewed from the coordinate system of the mounting part expressed using the tool parameter. A method for calibrating a tool parameter of a robot characterized by the following.
【請求項4】 上記ツールの位置決めにおける各姿勢の
うち,少なくとも1つは既知の方向とすることを特徴と
する請求項3記載のロボットのツールパラメータ較正方
法。
4. The method for calibrating tool parameters of a robot according to claim 3, wherein at least one of the postures for positioning the tool has a known direction.
JP6221558A 1994-09-16 1994-09-16 Robot tool parameter derivation method and calibration method Expired - Lifetime JP2774939B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP6221558A JP2774939B2 (en) 1994-09-16 1994-09-16 Robot tool parameter derivation method and calibration method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP6221558A JP2774939B2 (en) 1994-09-16 1994-09-16 Robot tool parameter derivation method and calibration method

Publications (2)

Publication Number Publication Date
JPH0885083A true JPH0885083A (en) 1996-04-02
JP2774939B2 JP2774939B2 (en) 1998-07-09

Family

ID=16768614

Family Applications (1)

Application Number Title Priority Date Filing Date
JP6221558A Expired - Lifetime JP2774939B2 (en) 1994-09-16 1994-09-16 Robot tool parameter derivation method and calibration method

Country Status (1)

Country Link
JP (1) JP2774939B2 (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009006452A (en) * 2007-06-29 2009-01-15 Nissan Motor Co Ltd Method for calibrating between camera and robot, and device therefor
DE102010032840A1 (en) 2009-08-03 2011-02-10 FANUC LTD, Oshino-mura Apparatus and method for measuring the position of a tool center of a robot
JP2012020347A (en) * 2010-07-12 2012-02-02 Kobe Steel Ltd Jig used for deriving tool vector of robot
JP2013223901A (en) * 2012-04-20 2013-10-31 Mitsubishi Electric Corp Automatic setting method of tool vector
JP2014076517A (en) * 2012-10-11 2014-05-01 Meiwa E Tec:Kk Error detection method of robot hand
CN104325467A (en) * 2014-11-05 2015-02-04 南京熊猫电子股份有限公司 Device for testing motion accuracy of robots
CN104842353A (en) * 2015-05-19 2015-08-19 北京理工大学 Ultrasonic transducer calibration-based manipulator positioning method, system and calibrating tool
JP2018001332A (en) * 2016-06-30 2018-01-11 セイコーエプソン株式会社 Robot, control device, and robot system
US9884425B2 (en) 2015-03-30 2018-02-06 Seiko Epson Corporation Robot, robot control device, and robotic system
US10551821B2 (en) 2016-06-30 2020-02-04 Seiko Epson Corporation Robot, robot control apparatus and robot system
US10882189B2 (en) 2017-04-21 2021-01-05 Seiko Epson Corporation Control device and robot system
EP3885079A1 (en) 2020-03-27 2021-09-29 Seiko Epson Corporation Teaching method
DE112021003659T5 (en) 2020-07-08 2023-04-27 Fanuc Corporation robot control device
DE112021005322T5 (en) 2020-12-21 2023-09-07 Fanuc Corporation Training data generating device, machine learning device and robot joint angle estimating device
US11759955B2 (en) 2020-03-19 2023-09-19 Seiko Epson Corporation Calibration method
EP4286107A1 (en) 2022-05-30 2023-12-06 Seiko Epson Corporation Calibration method and robot system

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011224672A (en) * 2010-04-15 2011-11-10 Kobe Steel Ltd Deriving method and calibration method for tool vector of robot
JP5645760B2 (en) 2011-06-21 2014-12-24 株式会社神戸製鋼所 Robot tool parameter correction method

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05189017A (en) * 1992-01-14 1993-07-30 Yaskawa Electric Corp Tool constant setting method for industrial robot

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05189017A (en) * 1992-01-14 1993-07-30 Yaskawa Electric Corp Tool constant setting method for industrial robot

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009006452A (en) * 2007-06-29 2009-01-15 Nissan Motor Co Ltd Method for calibrating between camera and robot, and device therefor
DE102010032840A1 (en) 2009-08-03 2011-02-10 FANUC LTD, Oshino-mura Apparatus and method for measuring the position of a tool center of a robot
JP2011031346A (en) * 2009-08-03 2011-02-17 Fanuc Ltd Apparatus and method for measuring position of tool end point of robot
JP4763074B2 (en) * 2009-08-03 2011-08-31 ファナック株式会社 Measuring device and measuring method of position of tool tip of robot
US9050728B2 (en) 2009-08-03 2015-06-09 Fanuc Ltd Apparatus and method for measuring tool center point position of robot
DE102010032840B4 (en) * 2009-08-03 2015-06-18 Fanuc Corporation Apparatus and method for measuring the position of a tool center of a robot
JP2012020347A (en) * 2010-07-12 2012-02-02 Kobe Steel Ltd Jig used for deriving tool vector of robot
JP2013223901A (en) * 2012-04-20 2013-10-31 Mitsubishi Electric Corp Automatic setting method of tool vector
JP2014076517A (en) * 2012-10-11 2014-05-01 Meiwa E Tec:Kk Error detection method of robot hand
CN104325467A (en) * 2014-11-05 2015-02-04 南京熊猫电子股份有限公司 Device for testing motion accuracy of robots
US9884425B2 (en) 2015-03-30 2018-02-06 Seiko Epson Corporation Robot, robot control device, and robotic system
CN104842353A (en) * 2015-05-19 2015-08-19 北京理工大学 Ultrasonic transducer calibration-based manipulator positioning method, system and calibrating tool
JP2018001332A (en) * 2016-06-30 2018-01-11 セイコーエプソン株式会社 Robot, control device, and robot system
US10551821B2 (en) 2016-06-30 2020-02-04 Seiko Epson Corporation Robot, robot control apparatus and robot system
US10882189B2 (en) 2017-04-21 2021-01-05 Seiko Epson Corporation Control device and robot system
US11759955B2 (en) 2020-03-19 2023-09-19 Seiko Epson Corporation Calibration method
EP3885079A1 (en) 2020-03-27 2021-09-29 Seiko Epson Corporation Teaching method
US11712803B2 (en) 2020-03-27 2023-08-01 Seiko Epson Corporation Teaching method
DE112021003659T5 (en) 2020-07-08 2023-04-27 Fanuc Corporation robot control device
DE112021005322T5 (en) 2020-12-21 2023-09-07 Fanuc Corporation Training data generating device, machine learning device and robot joint angle estimating device
EP4286107A1 (en) 2022-05-30 2023-12-06 Seiko Epson Corporation Calibration method and robot system

Also Published As

Publication number Publication date
JP2774939B2 (en) 1998-07-09

Similar Documents

Publication Publication Date Title
JP2774939B2 (en) Robot tool parameter derivation method and calibration method
US9517563B2 (en) Robot system using visual feedback
CN102458779B (en) Robot calibration apparatus and method for same
JP3600878B2 (en) Robot Position Correction Method Using Laser Measuring Machine
EP2767370A2 (en) Robot system and method for controlling the same
EP0358765B1 (en) Automatic setting of tool tip point
JP2005300230A (en) Measuring instrument
CN101239469A (en) Calibration device and method for robot mechanism
KR20080088165A (en) Robot calibration method
JP2010149267A (en) Robot calibration method and device
JPH09128549A (en) Relative position attitude detecting method for robot system
JP2011224672A (en) Deriving method and calibration method for tool vector of robot
KR100301231B1 (en) The automatic compensation method for robot working path
JPS61133409A (en) Automatic correction system of robot constant
JPH06278064A (en) Method for converting coordinate of visual sensor
JP6628170B1 (en) Measurement system and measurement method
JPH0871968A (en) Method for drawing tool parameter of robot and method for calibrating the same
JPH012104A (en) Robot positioning error correction method
JP7264763B2 (en) calibration device
JP3507032B2 (en) Method for deriving machine differences of robots
JPH08141950A (en) Calibration method of articulated robot
JP2000033586A (en) Error compensation method of robot and device thereof
WO2021062694A1 (en) Method, device and system for controlling robot, storage medium and terminal
RU2761924C1 (en) Manipulator control method
JPH08185216A (en) Tool attitude parameter setting method and robot controller

Legal Events

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

Free format text: PAYMENT UNTIL: 20080424

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20090424

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20100424

Year of fee payment: 12

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

Free format text: PAYMENT UNTIL: 20100424

Year of fee payment: 12

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

Free format text: PAYMENT UNTIL: 20110424

Year of fee payment: 13

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

Free format text: PAYMENT UNTIL: 20120424

Year of fee payment: 14

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

Free format text: PAYMENT UNTIL: 20130424

Year of fee payment: 15

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

Free format text: PAYMENT UNTIL: 20130424

Year of fee payment: 15

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

Free format text: PAYMENT UNTIL: 20140424

Year of fee payment: 16

EXPY Cancellation because of completion of term