JP7293905B2 - friction compensator - Google Patents

friction compensator Download PDF

Info

Publication number
JP7293905B2
JP7293905B2 JP2019116155A JP2019116155A JP7293905B2 JP 7293905 B2 JP7293905 B2 JP 7293905B2 JP 2019116155 A JP2019116155 A JP 2019116155A JP 2019116155 A JP2019116155 A JP 2019116155A JP 7293905 B2 JP7293905 B2 JP 7293905B2
Authority
JP
Japan
Prior art keywords
friction
model
speed
driven body
force
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2019116155A
Other languages
Japanese (ja)
Other versions
JP2021002248A (en
Inventor
金塀 徐
貴之 金子
崇 林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuji Electric Co Ltd filed Critical Fuji Electric Co Ltd
Priority to JP2019116155A priority Critical patent/JP7293905B2/en
Publication of JP2021002248A publication Critical patent/JP2021002248A/en
Application granted granted Critical
Publication of JP7293905B2 publication Critical patent/JP7293905B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position Or Direction (AREA)

Description

本発明は、被駆動体の位置や姿勢等の制御量を目標値に追従させるサーボ制御システムにおいて、特に、被駆動体の速度(運転方向)反転時に発生する機械系摩擦による影響を抑制するための摩擦補償技術に関するものである。 The present invention is intended for a servo control system in which a controlled variable such as the position and attitude of a driven object follows a target value, in particular to suppress the influence of mechanical friction that occurs when the speed (operating direction) of the driven object is reversed. Friction compensation technology.

サーボ制御システムにおいて、被駆動体の位置指令と位置検出値との誤差を補正するフィードバック制御系では、位置誤差を検出した後にその誤差に応じた駆動力が制御対象に与えられるため、被駆動体は位置誤差の発生時刻から遅れて応答することになる。また、サーボ制御システムの駆動系に内在している摩擦力は、特に被駆動体の速度(運転方向)反転時の特性が非線形となる。
これらの現象は、サーボ制御システムの運動制御精度を大きく低下させる原因となっている。
In the servo control system, in the feedback control system that corrects the error between the position command of the driven body and the position detection value, after the position error is detected, a driving force corresponding to the error is applied to the controlled body. responds with a delay from the time of occurrence of the position error. In addition, the frictional force inherent in the drive system of the servo control system exhibits non-linear characteristics particularly when the speed (operating direction) of the driven body is reversed.
These phenomena cause a significant deterioration in the motion control accuracy of the servo control system.

高精度な位置制御が要求される被駆動体については、速度反転時の摩擦による外乱が制御精度に与える影響を少なくすることが重要であり、この問題を解決する技術としては、従来から、フィードバック摩擦補償技術やモデルベースフィードフォワード摩擦補償技術が存在する。
例えば、外乱オブザーバを用いて摩擦を含む外乱を推定し、フィードバック制御により外乱抑圧補償を行って摩擦の影響を抑制することが可能であり、この種の外乱オブザーバによる非線形摩擦補償は速度反転時にも有効であることが知られている。
また、複数の状態量に対応する外乱モデルを用いて外乱推定値をそれぞれ演算し、これらの外乱推定値の和から外乱補正量を生成する技術も一般的に使用されている。
For a driven body that requires high-precision position control, it is important to reduce the influence of frictional disturbances on control accuracy during speed reversal. Friction compensation techniques and model-based feedforward friction compensation techniques exist.
For example, it is possible to estimate disturbances including friction using a disturbance observer and perform disturbance suppression compensation by feedback control to suppress the effects of friction. known to be effective.
A technique is also generally used in which disturbance estimates are calculated using disturbance models corresponding to a plurality of state quantities, and a disturbance correction amount is generated from the sum of these disturbance estimate values.

近年、モデルベース開発はサーボ業界に拡がりを見せており、特にフィードフォワード摩擦補償技術が多く採用されている。
例えば、図14は、速度のフィードバック制御系と位置のフィードフォワード補償器及び摩擦補償装置を備えたサーボ制御システムの機能ブロック図であり、特許文献1に開示されている。
同図において、210は位置制御器、220は電流指令iを生成する速度制御器、230はフィードフォワード補償器、240は摩擦力(摩擦トルク)Fを推定する摩擦補償装置、241はクーロン摩擦力等の各種パラメータを同定するパラメータ同定部、242は各種パラメータに基づいて摩擦力Fを推定する摩擦推定部、250は摩擦力Fを電流補償指令idisに換算する電流指令換算部、300は電力変換器等からなるドライバ(駆動装置)、400は電動機や被駆動体を含む制御対象としての工作機械を示す。
In recent years, model-based development has spread to the servo industry, and feedforward friction compensation technology has been widely adopted.
For example, FIG. 14 is a functional block diagram of a servo control system that includes a velocity feedback control system, a position feedforward compensator, and a friction compensator, and is disclosed in Patent Document 1.
In the figure, 210 is a position controller, 220 is a speed controller that generates a current command i * , 230 is a feedforward compensator, 240 is a friction compensator that estimates a friction force (friction torque) F, 241 is a Coulomb friction A parameter identification unit that identifies various parameters such as force, a friction estimation unit 242 that estimates the frictional force F based on various parameters, a current command conversion unit 250 that converts the frictional force F into a current compensation command idis , and a 300 A driver (driving device) 400 composed of a power converter or the like indicates a machine tool as a controlled object including an electric motor and a driven body.

図14におけるパラメータ同定部241では、電流指令i と、工作機械400の被駆動体の位置とに基づいて実機摩擦特性を生成し、この実機摩擦特性から、クーロン摩擦力F、最大静止摩擦力F、各要素の最大摩擦力の比率α、ばね定数kをパラメータとして同定する。そして、同定したパラメータを摩擦モデルに採用し、この摩擦モデルと被駆動体の速度とを用いて被駆動体の位置と摩擦力との関係を示す検証用摩擦特性を演算する。
そして、検証用摩擦特性と実機摩擦特性との誤差が許容範囲になるまでクーロン摩擦力F、最大静止摩擦力F等のパラメータを調整し、上記誤差が許容範囲内になった時の全ての同定パラメータk,α,F,Fを摩擦推定部242に出力する。
The parameter identification unit 241 in FIG. 14 generates the actual machine friction characteristics based on the current command i c * and the position of the driven body of the machine tool 400. From the actual machine friction characteristics, the Coulomb friction force F c , the maximum static The frictional force F s , the maximum frictional force ratio α i of each element, and the spring constant k i are identified as parameters. Then, the identified parameters are employed in the friction model, and the friction model and the speed of the driven body are used to calculate the verification friction characteristics indicating the relationship between the position of the driven body and the frictional force.
Then, parameters such as the Coulomb friction force F c and the maximum static friction force F s are adjusted until the error between the verification friction characteristic and the actual machine friction characteristic falls within the allowable range. to the friction estimator 242 .

摩擦推定部242では、入力された最新の同定パラメータを用いて摩擦力Fを推定し、電流指令換算部250は摩擦力Fを電流補償指令idisに換算して出力する。
この電流補償指令idisとフィードフォワード補償器230の出力とを用いて速度制御器220の出力iを補正することにより、摩擦補償された電流指令i がドライバ300に与えられることとなる。
The friction estimator 242 estimates the frictional force F using the latest identification parameter input, and the current command conversion unit 250 converts the frictional force F into a current compensation command i_dis and outputs it.
By correcting the output i * of the speed controller 220 using this current compensation command i dis and the output of the feedforward compensator 230, the friction-compensated current command i c * is given to the driver 300. .

特許第6214948号公報([0022]~[0029]、[0042]~[0046]、図1~図3,図11等)Japanese Patent No. 6214948 ([0022] to [0029], [0042] to [0046], Figures 1 to 3, Figure 11, etc.)

前述した特許文献1に係る先行技術には、以下のような問題がある。
まず、位置や速度の動的な特性を考慮すると、特許文献1に数式(1)~(4)として記載された摩擦モデルは構造が複雑であるため、高精度のモデル化及びパラメータ同定が困難になる恐れがある。
一方、一定条件のもとで同定したパラメータに基づいて計算した検証用摩擦特性と実機摩擦特性との間には、システムの同定方法、位置・速度検出器の精度、更にはノイズ等の影響によって必ず誤差が存在する。特許文献1には、上記の誤差が許容範囲内になるように各パラメータを調整する旨、記載されているが、その具体的な手順は何ら開示されていない。
The prior art according to Patent Literature 1 described above has the following problems.
First, considering the dynamic characteristics of position and speed, the structure of the friction model described as formulas (1) to (4) in Patent Document 1 is complicated, so it is difficult to model with high accuracy and identify parameters. there is a risk of becoming
On the other hand, the difference between the verification friction characteristics calculated based on the parameters identified under certain conditions and the actual machine friction characteristics depends on the system identification method, the accuracy of the position/speed detector, noise, etc. There will always be errors. Although Patent Literature 1 states that each parameter is adjusted so that the above error is within the allowable range, it does not disclose any specific procedure.

また、特許文献1では、速度や加速度に対するクーロン摩擦力の依存特性を表現できないため、運転状態の変化に伴って制御性能の再現性等のロバスト性を保証できない場合がある。しかも、速度が反転して被駆動体が反対方向に運動するような特殊条件のもとでは、クーロン摩擦力や最大静止摩擦力を調整しても、検証用摩擦特性と実機摩擦特性との誤差が許容範囲内にならない恐れがある。
更に、工作機械の駆動系等に応答遅れがあることを考え併せると、パラメータを調整して上記の誤差が許容範囲内になったとしても、摩擦の影響による制御精度の低下を常に解決できるとは限らない。
Further, in Patent Document 1, since the dependence characteristics of the Coulomb friction force on speed and acceleration cannot be expressed, there are cases where robustness such as reproducibility of control performance cannot be guaranteed with changes in operating conditions. Moreover, under special conditions where the speed is reversed and the driven body moves in the opposite direction, even if the Coulomb friction force and maximum static friction force are adjusted, the difference between the verification friction characteristics and the actual machine friction characteristics may not fall within the allowable range.
Furthermore, considering that there is a response delay in the drive system of the machine tool, even if the above error is within the allowable range by adjusting the parameters, it is possible to always solve the decrease in control accuracy due to the influence of friction. is not limited.

そこで、本発明の解決課題は、極力簡易な摩擦モデルを用いて摩擦静特性及び動特性の両者を同時に考慮し、被駆動体の速度(運動方向)が反転する場合でも摩擦力を高精度に推定して適切な摩擦補償を実現可能とした摩擦補償装置を提供することにある。 Therefore, the problem to be solved by the present invention is to use a friction model that is as simple as possible to simultaneously consider both static friction characteristics and dynamic characteristics, so that the friction force can be calculated with high accuracy even when the speed (movement direction) of the driven body is reversed. To provide a friction compensator capable of realizing appropriate friction compensation by estimation.

上記課題を解決するため、請求項1に係る発明は、被駆動体の位置や姿勢の制御量を目標値に追従させるサーボ制御システムにおいて前記被駆動体に作用する摩擦力を補償するための摩擦補償装置であって、
前記被駆動体の速度と、この速度に対する摩擦力の立ち上がり係数と、クーロン摩擦力と、粘性摩擦係数と、を含む摩擦モデルを用いて、速度フィードバック制御系の出力に加算する摩擦力を摩擦補償量として推定する摩擦補償量推定手段を備え、
前記摩擦補償量推定手段は、
前記サーボ制御システムの入出力情報と前記摩擦モデルとを用いて、前記立ち上がり係数、前記クーロン摩擦力及び前記粘性摩擦係数を前記摩擦モデルのパラメータとして同定し、同定したパラメータを初期値として異なる加速度条件のもとで前記パラメータを更新して出力するモデルパラメータ演算部と、
前記モデルパラメータ演算部から出力された前記パラメータを用いて前記摩擦補償量を初期化し、前記被駆動体の速度反転時近傍における前記制御量の誤差が最小となるように、前記クーロン摩擦力を更新すると共に、複数の速度テスト信号にそれぞれ対応する複数の加速度と前記立ち上がり係数との関係を取得し、この関係と前記被駆動体の速度及び加速度と前記クーロン摩擦力とに基づいて前記摩擦補償量を推定するモデルベース摩擦補償部と、
を備えたことを特徴とする。
In order to solve the above-mentioned problems, the invention according to claim 1 provides a friction control system for compensating for the frictional force acting on the driven body in a servo control system that causes the control amount of the position and attitude of the driven body to follow a target value. A compensator,
Friction compensation for the friction force added to the output of the speed feedback control system using a friction model including the speed of the driven body, the rise coefficient of the friction force with respect to this speed, the Coulomb friction force, and the viscous friction coefficient Equipped with friction compensation amount estimation means for estimating as an amount,
The friction compensation amount estimating means is
Using the input/output information of the servo control system and the friction model, the rise coefficient, the Coulomb friction force, and the viscous friction coefficient are identified as parameters of the friction model, and the identified parameters are used as initial values for different acceleration conditions. a model parameter calculator that updates and outputs the parameters under
The friction compensation amount is initialized using the parameter output from the model parameter calculation unit, and the Coulomb friction force is updated so that an error in the control amount near the speed reversal of the driven body is minimized. Then, obtain the relationship between a plurality of accelerations corresponding to a plurality of speed test signals and the rise coefficient, and determine the friction compensation amount based on this relationship, the speed and acceleration of the driven body, and the Coulomb friction force. a model-based friction compensator that estimates
characterized by comprising

請求項2に係る発明は、請求項1に記載した摩擦補償装置において、前記モデルパラメータ演算部は、前記被駆動体の速度または加速度と摩擦力との関係を示す実機摩擦特性に基づいて前記摩擦モデルのパラメータを同定することを特徴とする。 The invention according to claim 2 is the friction compensating device according to claim 1, wherein the model parameter calculation unit calculates the frictional force based on the actual machine frictional characteristic indicating the relationship between the speed or acceleration of the driven body and the frictional force. It is characterized by identifying the parameters of the model.

請求項3に係る発明は、請求項1または2に記載した摩擦補償装置において、前記摩擦モデルが、前記被駆動体の速度が正から負、または、負から正に変化する不連続性を連続関数により近似したモデルであることを特徴とする。 The invention according to claim 3 is the friction compensating device according to claim 1 or 2, wherein the friction model continuously detects a discontinuity in which the speed of the driven body changes from positive to negative or from negative to positive. It is characterized by being a model approximated by a function.

請求項4に係る発明は、請求項1~3の何れか1項に記載した摩擦補償装置において、前記摩擦モデルのパラメータとして、前記被駆動体の加速度と前記立ち上がり係数との関係におけるヒステリシス幅を含むことを特徴とする。 The invention according to claim 4 is the friction compensation device according to any one of claims 1 to 3, wherein the hysteresis width in the relationship between the acceleration of the driven body and the rise coefficient is used as a parameter of the friction model. characterized by comprising

本発明によれば、運転状況に応じて摩擦補償量を動的に決定し、モデルベースフィードフォワード摩擦補償を行うことにより、速度反転時の摩擦の影響による制御誤差を低減させて適切な摩擦補償を実現することができる。 According to the present invention, the amount of friction compensation is dynamically determined according to the operating conditions, and model-based feedforward friction compensation is performed, thereby reducing control errors due to the effects of friction during speed reversal and providing appropriate friction compensation. can be realized.

本発明の実施形態に係るサーボ制御システムの一構成例を示す機能ブロック図である。1 is a functional block diagram showing one configuration example of a servo control system according to an embodiment of the present invention; FIG. 本発明の実施形態における速度-摩擦力特性の一例を示す図である。FIG. 4 is a diagram showing an example of velocity-frictional force characteristics in an embodiment of the present invention; 図1におけるモデルベース摩擦補償部の機能ブロック図である。FIG. 2 is a functional block diagram of a model-based friction compensator in FIG. 1; 本発明の実施形態において、ヒステリシス特性を有する速度-摩擦力特性を示す図である。FIG. 4 is a diagram showing velocity-frictional force characteristics with hysteresis characteristics in an embodiment of the present invention; 図1におけるモデルパラメータ演算部の機能ブロック図である。2 is a functional block diagram of a model parameter calculator in FIG. 1; FIG. 本発明の実施形態における加速度と立ち上がり係数との関係を示す図である。It is a figure which shows the relationship between the acceleration and a rise coefficient in embodiment of this invention. 本発明の実施形態における速度テスト信号の説明図である。FIG. 4 is an explanatory diagram of a speed test signal in the embodiment of the present invention; 本発明の実施形態におけるパラメータチューニング動作を示すフローチャートである。4 is a flow chart showing parameter tuning operation in the embodiment of the present invention; 本発明の実施形態における加速度と立ち上がり係数との関係を示す図である。It is a figure which shows the relationship between the acceleration and a rise coefficient in embodiment of this invention. 本発明の実施形態における摩擦補償量の演算手順を示すフローチャートである。4 is a flow chart showing a procedure for calculating a friction compensation amount in the embodiment of the present invention; 実機を対象とした、速度反転時近傍における加速度に応じた摩擦力の測定結果を示す図である。FIG. 10 is a diagram showing measurement results of frictional force according to acceleration near the time of speed reversal, for an actual machine; 本発明の実施形態を適用した場合の、速度反転時近傍における加速度に応じた摩擦力(摩擦補償量)の推定結果を示す図である。FIG. 7 is a diagram showing estimation results of frictional force (friction compensation amount) according to acceleration near the time of speed reversal when the embodiment of the present invention is applied; 本発明の実施形態による摩擦補償有りの場合及び摩擦補償無しの場合の、速度追従性のシミュレーション結果を示す図である。FIG. 5 is a diagram showing simulation results of velocity followability with and without friction compensation according to the embodiment of the present invention; 特許文献1に記載されたサーボ制御システムの機能ブロック図である。1 is a functional block diagram of a servo control system described in Patent Document 1; FIG.

以下、図に沿って本発明の実施形態を説明する。なお、本発明の技術的範囲は以下に説明する実施形態によって何ら限定されるものではない。
まず、図1は、この実施形態の摩擦補償装置を含むサーボ制御システムの一構成例を示す機能ブロック図であり、例えば、被駆動体を有する工作機械及びその駆動装置を含む制御対象30(伝達関数Gとする)に電流指令等を与えて被駆動体の位置を制御するシステムを想定している。
An embodiment of the present invention will be described below with reference to the drawings. In addition, the technical scope of the present invention is not limited at all by the embodiments described below.
First, FIG. 1 is a functional block diagram showing one configuration example of a servo control system including the friction compensator of this embodiment. A system is assumed in which a current command or the like is given to a function G) to control the position of the driven body.

図1において、上位の制御装置から送られた位置指令と制御対象30の位置出力との偏差が減算手段81により演算され、位置フィードバック制御器50は上記偏差をゼロにするように動作して速度指令を出力する。制御対象30の位置出力は微分手段91(sはラプラス演算子)により微分されて速度検出値に変換され、減算手段82は速度指令と速度検出値との偏差を演算すると共に、速度フィードバック制御器20は上記偏差をゼロにするように動作して電流指令(またはトルク指令)を出力する。 In FIG. 1, the deviation between the position command sent from the host control device and the position output of the controlled object 30 is calculated by the subtracting means 81, and the position feedback controller 50 operates to make the deviation zero, thereby increasing the speed. Output directive. The position output of the controlled object 30 is differentiated by the differentiating means 91 (s is the Laplace operator) and converted into a detected speed value. The subtracting means 82 calculates the deviation between the speed command and the detected speed value, 20 operates to make the deviation zero and outputs a current command (or torque command).

また、前記位置指令は、フィードフォワード補償用の逆システム50(制御対象30の伝達関数Gに対して、G-1とする)に入力され、その出力が加算手段83により速度フィードバック制御器40の出力に加算される。
加算手段83の出力は、加算手段84、減算手段85、加算手段86を経て制御対象30に与えられる。なお、加算手段84では、後述するモデルベース摩擦補償部70の出力が加算手段83の出力に加算され、減算手段85では、加算手段84の出力から外乱オブザーバ40による推定外乱が減算されると共に、加算手段86では減算手段85の出力と外乱信号dとが加算される。
In addition, the position command is input to an inverse system 50 for feedforward compensation (assumed to be G −1 for the transfer function G of the controlled object 30 ), and its output is sent to the velocity feedback controller 40 by the addition means 83 . Added to the output.
The output of the adding means 83 is given to the controlled object 30 via the adding means 84, the subtracting means 85 and the adding means 86. The adding means 84 adds the output of the model-based friction compensator 70, which will be described later, to the output of the adding means 83, and the subtracting means 85 subtracts the disturbance estimated by the disturbance observer 40 from the output of the adding means 84. The adding means 86 adds the output of the subtracting means 85 and the disturbance signal d.

更に、前記位置指令は微分手段93及び二階微分手段92に入力されて速度指令及び加速度指令が演算され、これらが摩擦に起因する制御誤差を補償するためのモデルベース摩擦補償部70の摩擦モデルに入力される。また、二階微分手段92の出力である加速度指令はモデルパラメータ演算部60に入力され、この演算部60により同定された摩擦モデルのパラメータがモデルベース摩擦補償部70に入力されている。モデルベース摩擦補償部70では、同定されたパラメータと速度指令及び加速度指令に基づいて補償するべき摩擦力(摩擦補償量)を演算し、前記加算手段84に出力する。
ここで、モデルパラメータ演算部60及びモデルベース摩擦補償部70は、特許請求の範囲における摩擦補償量推定手段を構成している。
Furthermore, the position command is input to the differentiating means 93 and the second-order differentiating means 92 to calculate a velocity command and an acceleration command, which are applied to the friction model of the model-based friction compensator 70 for compensating for control errors caused by friction. is entered. The acceleration command, which is the output of the second-order differentiating means 92 , is input to the model parameter calculator 60 , and the parameters of the friction model identified by this calculator 60 are input to the model-based friction compensator 70 . The model-based friction compensator 70 calculates the frictional force (friction compensation amount) to be compensated based on the identified parameters, the speed command and the acceleration command, and outputs it to the adding means 84 .
Here, the model parameter calculator 60 and the model-based friction compensator 70 constitute friction compensation amount estimating means in the claims.

モデルパラメータ演算部60は、例えば、以下の数式1に示すような簡易な摩擦モデルを用いて複数のパラメータを同定し、モデルベース摩擦補償部70はこれらのパラメータに基づいて摩擦力Tを推定する。
[数式1]
=Erf(αω)Tcf+Dω
ここで、Erfは誤差関数、αは摩擦力の立ち上がりの傾きを表す立ち上がり係数、ωは被駆動体の速度、Tcfはクーロン摩擦力、Dは粘性摩擦係数である。
The model parameter calculator 60 identifies a plurality of parameters using, for example, a simple friction model as shown in Equation 1 below, and the model-based friction compensator 70 estimates the frictional force T f based on these parameters. do.
[Formula 1]
Tf =Erf(αω) Tcf +Dω
Here, Erf is the error function, α is the rise coefficient representing the slope of rise of the friction force, ω is the speed of the driven body, Tcf is the Coulomb friction force, and D is the viscous friction coefficient.

被駆動体の速度ゼロを境として、クーロン摩擦力は正から負、または負から正に不連続に変化するのに対して、数式1のように誤差関数Erfを用いて摩擦力Tを近似することにより、図2の速度-摩擦力特性に示すように、速度がゼロに近付くにつれて所定の傾きにより滑らかに変化する静的摩擦特性を得ることができる。
なお、モデルパラメータ演算部60の機能については、図5と共に後述する。
While the Coulomb frictional force changes discontinuously from positive to negative or from negative to positive with the zero speed of the driven body as a boundary, the frictional force Tf can be approximated using the error function Erf as shown in Equation 1. By doing so, as shown in the speed-friction force characteristics of FIG. 2, static friction characteristics that smoothly change with a predetermined slope as the speed approaches zero can be obtained.
Note that the function of the model parameter calculator 60 will be described later together with FIG.

図1のモデルベース摩擦補償部70は、図3の機能ブロック図に示すように、モデルベース摩擦補償量計算部71、位相調整部72、ゲイン調整及び出力制限部73、及び出力フィルタ74を備えている。 As shown in the functional block diagram of FIG. 3, the model-based friction compensator 70 of FIG. ing.

モデルベース摩擦補償量計算部71は、モデルパラメータ演算部60により同定された摩擦モデルのパラメータと速度情報(速度指令及び加速度指令)とに基づいて、補償するべき摩擦力を計算する。
位相調整部72は、例えば、図2の速度-摩擦力特性を図4のごとく所定のヒステリシスを有する特性に調整する。すなわち、速度が負から正に反転する際、及び、正から負に反転する際に、摩擦力がクーロン摩擦力Tcfに等しくなるように摩擦力の位相を調整し、速度-摩擦力特性を所定の傾きで平滑化する。これにより、サーボ制御システムの応答性を高めて実機の摩擦力のヒステリシス特性を補償すると共に、速度反転前後における摩擦補償量を安定化することができる。なお、ヒステリシス幅は、位相調整部72が摩擦力の位相を調整することで任意に設定可能であり、このヒステリシス幅を摩擦モデルのパラメータに含めても良い。
勿論、上記に代えて、ヒステリシス特性補償及び位相補償を考慮した摩擦モデリングを行っても良い。
The model-based friction compensation amount calculator 71 calculates the frictional force to be compensated based on the parameters of the friction model identified by the model parameter calculator 60 and speed information (speed command and acceleration command).
The phase adjustment unit 72 adjusts the velocity-frictional force characteristic of FIG. 2 to a characteristic having a predetermined hysteresis as shown in FIG. 4, for example. That is, when the speed reverses from negative to positive and from positive to negative, the phase of the frictional force is adjusted so that the frictional force becomes equal to the Coulomb frictional force Tcf , and the speed-frictional force characteristic is Smooth with a given slope. As a result, the responsiveness of the servo control system can be enhanced to compensate for the hysteresis characteristic of the frictional force of the actual machine, and the friction compensation amount before and after the speed reversal can be stabilized. The hysteresis width can be arbitrarily set by adjusting the phase of the frictional force by the phase adjustment unit 72, and the hysteresis width may be included in the parameters of the friction model.
Of course, instead of the above, friction modeling considering hysteresis characteristic compensation and phase compensation may be performed.

図3に戻って、位相調整部72に後続するゲイン調整及び出力制限部73は、補償するべき摩擦力の数値的な安定性を保つと共に過剰補償を防止するためのものである。具体的には、モデルパラメータ演算部60により同定した各パラメータを前述した数式1の摩擦モデルに代入することにより計算した摩擦力の初期値の複数倍、例えば3倍程度の大きさに変換する。
ゲイン調整及び出力制限部73に後続する出力フィルタ74は、例えばローパスフィルタ特性を有し、この出力フィルタ74を通して得られた摩擦補償量が前記加算手段84に入力され、速度フィードバック制御器20の出力(詳しくは加算手段83の出力)に加算されてその結果が前記加算手段85に入力される。
Returning to FIG. 3, the gain adjustment and output limiter 73 following the phase adjuster 72 is for maintaining numerical stability of the frictional force to be compensated and preventing overcompensation. Specifically, each parameter identified by the model parameter calculation unit 60 is substituted into the friction model of the above-described Equation 1 to convert the initial value of the calculated frictional force to a magnitude that is multiple times, for example, about three times.
An output filter 74 subsequent to the gain adjustment and output limiter 73 has, for example, a low-pass filter characteristic. (more specifically, the output of the adding means 83) and the result is input to the adding means 85.

次に、モデルパラメータ演算部60の機能について、図5を参照しつつ説明する。
モデルパラメータ演算部60では、図5に示すように、数式1の摩擦モデルと事前に取得したサーボ制御システムの入出力情報とを用いて、システム同定により摩擦モデルの各種パラメータを同定する(ステップS1)。具体的には、例えば、被駆動体の速度または加速度と摩擦力との関係を示す実機摩擦特性に基づいて各種パラメータを同定する。これらの同定したパラメータを初期値として、異なる加速度条件を与えてチューニングにより動的なパラメータを演算する(同S2)。次いで、加速度とパラメータとの関係を記録しておき(同S3)、実際の運転時の加速度情報を上述の記録済みの関係と照合して動的にパラメータを更新し(同S4)、更新されたパラメータをモデルベース摩擦補償部70に出力する。
Next, functions of the model parameter calculator 60 will be described with reference to FIG.
As shown in FIG. 5, the model parameter calculation unit 60 identifies various parameters of the friction model by system identification using the friction model of Equation 1 and previously acquired input/output information of the servo control system (step S1 ). Specifically, for example, various parameters are identified based on actual machine friction characteristics that indicate the relationship between the speed or acceleration of the driven body and the frictional force. Using these identified parameters as initial values, different acceleration conditions are given and dynamic parameters are calculated by tuning (S2). Next, the relationship between the acceleration and the parameter is recorded (S3), and the acceleration information during actual driving is collated with the recorded relationship to dynamically update the parameter (S4). The parameters obtained are output to the model-based friction compensator 70 .

ここで、数式1は、速度に対する摩擦力の静特性のみを考慮した摩擦モデルを対象としている。例えば、図4のように速度-摩擦力特性を調整したとしても、摩擦補償量の立ち上がり係数は一定値であり、実際の運転時において図6のように加速度が種々変化する場合には、速度反転時近傍の微小速度領域における摩擦力を正確に再現することができない。
このため、本実施形態では、速度反転時近傍の微小速度領域における摩擦力の動的な挙動を反映させるように、上述したモデルパラメータ演算部60のチューニングによって動的にパラメータを取得することとした。
Here, Equation 1 is intended for a friction model that considers only static characteristics of frictional force with respect to velocity. For example, even if the speed-friction force characteristics are adjusted as shown in FIG. 4, the rise coefficient of the friction compensation amount is a constant value. The frictional force in the very small speed region near the time of reversal cannot be reproduced accurately.
For this reason, in the present embodiment, parameters are dynamically obtained by tuning the model parameter calculation unit 60 described above so as to reflect the dynamic behavior of the frictional force in the micro-velocity region near the time of velocity reversal. .

このモデルパラメータ演算部60によるチューニング動作について、図7,図8を参照しつつ更に説明する。
例えば、図7に示す台形状の速度テスト信号を用い(±ωは最大速度、Tは加速時間)、サーボ制御システムの運転可能な速度範囲(図7における-ω~+ωの範囲)にわたって、加速度が異なる第一速度テスト信号及び第二速度テスト信号を事前に定義しておく。
The tuning operation by this model parameter calculator 60 will be further described with reference to FIGS. 7 and 8. FIG.
For example, using the trapezoidal speed test signal shown in FIG. 7 (±ω r is the maximum speed and T is the acceleration time), the operable speed range of the servo control system (range from −ω r to +ω r in FIG. 7) A first velocity test signal and a second velocity test signal with different accelerations are predefined.

そして、まず、同定したパラメータを用いて数式1の推定摩擦力Tを初期化する(図8のステップS11)。初期化が完了したら、第一速度テスト信号によりサーボ制御システムを駆動する(同S12)。次に、摩擦モデルのパラメータのうちクーロン摩擦力Tcfの設定値を調整しながら(同S13)、速度反転時近傍における制御量の誤差(例えば、被駆動体の位置目標値と位置検出値との誤差)が最小になるまで上記設定値の調整を繰り返し(同S14)、制御量の誤差が最小になったら(同S14Yes)、その時の設定値によりクーロン摩擦力を更新する(同S15)。
次に、摩擦補償量の立ち上がり係数αを調整する(同S16)。そして、速度反転時近傍における制御量の誤差が最小になるまで立ち上がり係数αの調整を繰り返し(同S17)、制御量の誤差が最小になったら(同S17Yes)、その時の立ち上がり係数αを第一係数αとして保存する(同S18)。
Then, first, the estimated frictional force T f of Equation 1 is initialized using the identified parameters (step S11 in FIG. 8). After the initialization is completed, the servo control system is driven by the first speed test signal (S12). Next, while adjusting the set value of the Coulomb friction force Tcf among the parameters of the friction model (at step S13), the control amount error (for example, the position target value of the driven body and the position detection value error) is minimized (S14), and when the control amount error is minimized (Yes at S14), the Coulomb friction force is updated with the set value at that time (S15).
Next, the rise coefficient α of the friction compensation amount is adjusted (S16). Then, the rise coefficient α is repeatedly adjusted until the error in the control amount near the speed reversal is minimized (S17), and when the error in the control amount is minimized (Yes in S17), the rise coefficient It is stored as the coefficient α1 (at S18).

続いて、第二速度テスト信号によりサーボ制御システムを駆動して立ち上がり係数αを調整し(同S19,S20)、速度反転時近傍における制御量の誤差が最小になるまで立ち上がり係数αの調整を繰り返し(同S21)、制御量の誤差が最小になったら(同S21Yes)、その時の立ち上がり係数αを第二係数αとして保存する(同S22)。
以上の処理によって保存した第一係数α及び第二係数αと立ち上がり係数の上下限値とを用いて、図9に示す如く、一次線形近似関数により、加速度aと立ち上がり係数αとの関係を求める(同S23)。
Subsequently, the servo control system is driven by the second speed test signal to adjust the rise coefficient α (at S19 and S20), and the adjustment of the rise coefficient α is repeated until the error in the control amount near the speed reversal is minimized. (S21), when the control amount error becomes minimum (Yes in S21), the rising coefficient α at that time is stored as the second coefficient α2 (S22).
Using the first coefficient α1 and the second coefficient α2 saved by the above processing and the upper and lower limits of the rise coefficient, the relationship between the acceleration a and the rise coefficient α is obtained by a first-order linear approximation function as shown in FIG. is obtained (same S23).

次に、図10は、図1のモデルベース摩擦補償部70によって摩擦補償量を推定する動作を示すフローチャートである。
まず、運転中のサーボ制御システムの現時点における速度や加速度等の運転情報を取得し(ステップS31)、これらの運転情報を用いて摩擦モデルの各パラメータを更新すると共に(同S32)、これらのパラメータに基づき数式1により摩擦補償量を計算して出力する(同S33,S34)。
具体的には、図8の処理により求めた第一係数α及び第二係数αと立ち上がり係数αの上下限値とから、図9に示すような加速度a(言い換えれば速度ω)と立ち上がり係数αとの関係が求められ、また、図8のステップS15によりクーロン摩擦力が求められているので、これらのパラメータを前述の数式1に代入して摩擦力(摩擦補償量)Tを求めることができる。
ここで、サーボ制御システムの経年的な特性変化や周囲の環境変化等がある場合には、前述したモデルパラメータ演算部60のチューニング動作による動的なパラメータの演算(図5のステップS2)を再度実施すれば対応可能である。
Next, FIG. 10 is a flow chart showing the operation of estimating the friction compensation amount by the model-based friction compensator 70 of FIG.
First, the operating information such as the current speed and acceleration of the servo control system in operation is acquired (step S31), each parameter of the friction model is updated using this operating information (step S32), and these parameters Based on this, the friction compensation amount is calculated by the formula 1 and output (at steps S33 and S34).
Specifically, from the first coefficient α1 and the second coefficient α2 obtained by the processing of FIG. Since the relationship with the coefficient α has been obtained and the Coulomb friction force has been obtained in step S15 of FIG . be able to.
Here, if there is a change in the characteristics of the servo control system over time or a change in the surrounding environment, dynamic parameter calculation (step S2 in FIG. 5) by the tuning operation of the model parameter calculation unit 60 described above is performed again. It can be handled if implemented.

なお、図11は、実機を対象とした、速度反転時近傍における加速度に応じた摩擦力の測定結果を示す図である。図11によれば、加速度が小さいほど、速度反転時近傍における摩擦力の変動が大きくなっていることがわかる。
これに対し、図12は、本発明の実施形態を適用した場合の、速度反転時近傍における加速度に応じた摩擦力(摩擦補償量)の推定結果を示す図である。本発明の実施形態では、図9に示した加速度aと立ち上がり係数αとの関係を用いて数式1により推定摩擦力Tを求め、この推定摩擦力Tにより摩擦補償を行うため、加速度の大きさに限らず滑らかに変化する摩擦補償量を動的に求めることができ、摩擦の影響による制御誤差を低減することができる。
11A and 11B are diagrams showing the measurement results of the frictional force according to the acceleration in the vicinity of the speed reversal for the actual machine. According to FIG. 11, it can be seen that the smaller the acceleration, the larger the fluctuation of the frictional force near the speed reversal.
On the other hand, FIG. 12 is a diagram showing the estimation result of the frictional force (friction compensation amount) according to the acceleration near the speed reversal when the embodiment of the present invention is applied. In the embodiment of the present invention, the estimated frictional force Tf is obtained by Equation 1 using the relationship between the acceleration a and the rising coefficient α shown in FIG. It is possible to dynamically obtain a smoothly changing friction compensation amount regardless of the magnitude, and to reduce control errors due to the influence of friction.

更に、図13は、本発明の実施による摩擦補償有りの場合、及び、摩擦補償無しの場合の、速度追従性のシミュレーション結果を示す図である。
本発明の実施形態によって摩擦補償を行った場合には、速度反転時の前後において摩擦の影響を受けずに速度が連続的に変化しているのに対し、摩擦補償を行わない場合には、摩擦によって速度変化を妨げられる期間が生じていることがわかる。
Furthermore, FIG. 13 is a diagram showing simulation results of velocity followability with and without friction compensation according to the implementation of the present invention.
When friction compensation is performed according to the embodiment of the present invention, the speed changes continuously without being affected by friction before and after the speed reversal. It can be seen that there is a period during which speed change is hindered by friction.

以上説明したように、本実施形態によれば、摩擦の静特性及び動特性を考慮したうえで、異なる加速度に対して摩擦補償量を高精度に推定することができ、再現性の高い位置制御や速度制御を行うことができる。 As described above, according to the present embodiment, the friction compensation amount can be estimated with high accuracy for different accelerations in consideration of the static and dynamic characteristics of friction. and speed control.

10:位置フィードバック制御器
20:速度フィードバック制御器
30:制御対象
40:外乱オブザーバ
50:逆システム
60:モデルパラメータ演算部
70:モデルベース摩擦補償部
81,82,85:減算手段
83,84,86:加算手段
91,93:微分手段
92:二階微分手段
10: Position feedback controller 20: Velocity feedback controller 30: Control object 40: Disturbance observer 50: Inverse system 60: Model parameter calculator 70: Model-based friction compensator 81, 82, 85: Subtracting means 83, 84, 86 : addition means 91, 93: differentiation means 92: second-order differentiation means

Claims (4)

被駆動体の位置や姿勢の制御量を目標値に追従させるサーボ制御システムにおいて前記被駆動体に作用する摩擦力を補償するための摩擦補償装置であって、
前記被駆動体の速度と、この速度に対する摩擦力の立ち上がり係数と、クーロン摩擦力と、粘性摩擦係数と、を含む摩擦モデルを用いて、速度フィードバック制御系の出力に加算する摩擦力を摩擦補償量として推定する摩擦補償量推定手段を備え、
前記摩擦補償量推定手段は、
前記サーボ制御システムの入出力情報と前記摩擦モデルとを用いて、前記立ち上がり係数、前記クーロン摩擦力及び前記粘性摩擦係数を前記摩擦モデルのパラメータとして同定し、同定したパラメータを初期値として異なる加速度条件のもとで前記パラメータを更新して出力するモデルパラメータ演算部と、
前記モデルパラメータ演算部から出力された前記パラメータを用いて前記摩擦補償量を初期化し、前記被駆動体の速度反転時近傍における前記制御量の誤差が最小となるように、前記クーロン摩擦力を更新すると共に、複数の速度テスト信号にそれぞれ対応する複数の加速度と前記立ち上がり係数との関係を取得し、この関係と前記被駆動体の速度及び加速度と前記クーロン摩擦力とに基づいて前記摩擦補償量を推定するモデルベース摩擦補償部と、
を備えたことを特徴とする摩擦補償装置。
A friction compensator for compensating for frictional force acting on a driven body in a servo control system that causes a control amount of the position and attitude of the driven body to follow target values,
Friction compensation for the friction force added to the output of the speed feedback control system using a friction model including the speed of the driven body, the rise coefficient of the friction force with respect to this speed, the Coulomb friction force, and the viscous friction coefficient Equipped with friction compensation amount estimation means for estimating as an amount,
The friction compensation amount estimating means is
Using the input/output information of the servo control system and the friction model, the rise coefficient, the Coulomb friction force, and the viscous friction coefficient are identified as parameters of the friction model, and the identified parameters are used as initial values for different acceleration conditions. a model parameter calculator that updates and outputs the parameters under
The friction compensation amount is initialized using the parameter output from the model parameter calculation unit, and the Coulomb friction force is updated so that an error in the control amount near the speed reversal of the driven body is minimized. Further, obtaining relationships between a plurality of accelerations corresponding to a plurality of speed test signals and the rise coefficient, and calculating the friction compensation amount based on the relationships, the speed and acceleration of the driven body, and the Coulomb friction force. a model-based friction compensator that estimates
A friction compensator comprising:
請求項1に記載した摩擦補償装置において、
前記モデルパラメータ演算部は、前記被駆動体の速度または加速度と摩擦力との関係を示す実機摩擦特性に基づいて前記摩擦モデルのパラメータを同定することを特徴とする摩擦補償装置。
In the friction compensator according to claim 1,
The friction compensator according to claim 1, wherein the model parameter calculator identifies the parameters of the friction model based on actual machine friction characteristics indicating the relationship between the velocity or acceleration of the driven body and the friction force.
請求項1または2に記載した摩擦補償装置において、
前記摩擦モデルが、前記被駆動体の速度が正から負、または、負から正に変化する不連続性を連続関数により近似したモデルであることを特徴とする摩擦補償装置。
In the friction compensation device according to claim 1 or 2,
A friction compensation device according to claim 1, wherein the friction model is a model that approximates a discontinuity in which the speed of the driven body changes from positive to negative or from negative to positive with a continuous function.
請求項1~3の何れか1項に記載した摩擦補償装置において、
前記摩擦モデルのパラメータとして、前記被駆動体の加速度と前記立ち上がり係数との関係におけるヒステリシス幅を含むことを特徴とする摩擦補償装置。
In the friction compensation device according to any one of claims 1 to 3,
A friction compensator according to claim 1, wherein a hysteresis width in relation between the acceleration of the driven body and the rise coefficient is included as a parameter of the friction model.
JP2019116155A 2019-06-24 2019-06-24 friction compensator Active JP7293905B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019116155A JP7293905B2 (en) 2019-06-24 2019-06-24 friction compensator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019116155A JP7293905B2 (en) 2019-06-24 2019-06-24 friction compensator

Publications (2)

Publication Number Publication Date
JP2021002248A JP2021002248A (en) 2021-01-07
JP7293905B2 true JP7293905B2 (en) 2023-06-20

Family

ID=73995095

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019116155A Active JP7293905B2 (en) 2019-06-24 2019-06-24 friction compensator

Country Status (1)

Country Link
JP (1) JP7293905B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113110305B (en) * 2021-04-30 2022-07-05 华中科技大学 Friction modeling method of electromechanical system and application thereof
CN114291051B (en) * 2022-01-05 2023-03-28 清华大学 Method and device for modeling internal friction force of motor servo type hydraulic line control brake system
CN114545864A (en) * 2022-03-07 2022-05-27 哈尔滨理工大学 Friction compensation control method of electro-hydraulic servo system
CN115356928A (en) * 2022-08-18 2022-11-18 北京华卓精科科技股份有限公司 Friction compensation method and device for linear guide rail displacement system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001231280A (en) 2000-02-10 2001-08-24 Tamagawa Seiki Co Ltd Method and apparatus for controlling position
JP2011186669A (en) 2010-03-05 2011-09-22 Niigata Machine Techno Co Ltd Pressure control device and method
JP2016164697A (en) 2015-03-06 2016-09-08 富士電機株式会社 Position control system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001231280A (en) 2000-02-10 2001-08-24 Tamagawa Seiki Co Ltd Method and apparatus for controlling position
JP2011186669A (en) 2010-03-05 2011-09-22 Niigata Machine Techno Co Ltd Pressure control device and method
JP2016164697A (en) 2015-03-06 2016-09-08 富士電機株式会社 Position control system

Also Published As

Publication number Publication date
JP2021002248A (en) 2021-01-07

Similar Documents

Publication Publication Date Title
JP7293905B2 (en) friction compensator
JP2005135186A (en) Reference model follow-up type control system and its method
CN102969971B (en) Motor control apparatus
JP5279299B2 (en) POSITION CONTROL APPARATUS, EXPOSURE APPARATUS, AND DEVICE MANUFACTURING METHOD EQUIPPED WITH AN ITS
JP5847338B2 (en) Motor control device
JP3506157B2 (en) Motor position control device
US20160363939A1 (en) Control device for machine apparatus and gain determination method for friction compensation
JP5780473B2 (en) Motor control device and motor system
KR101813615B1 (en) Apparatus and method for control of actuator
CN110504880B (en) Feedforward compensation control method for interference observation of flux switching permanent magnet linear motor
JPWO2014091840A1 (en) Servo control device
US8442692B2 (en) Method for performing adaptive friction compensation in an actuator accounting for variation in friction characteristics of wave gear drive accompanying change in temperature
JP2021047556A (en) Object moving device, machine tool, information processing apparatus, information processing method, and information processing program
JP6930868B2 (en) Servo control device, servo control method and system
JP5017984B2 (en) Servo control device and speed tracking control method thereof
JP7120821B2 (en) Control device, control method and program
JP2019221032A (en) Electric motor control device
CN111103793A (en) Control method and device of gun recoil buffer system based on magnetorheological damper
JP2005071034A (en) Servo controller
JP3454942B2 (en) Position control device
KR102240723B1 (en) Controller comprising position predict unit and a control method thereof
JP7020200B2 (en) Speed / position control system
CN116339177A (en) Robot control method based on dynamic friction compensation
JPH04112690A (en) Control method for servo motor
CN109067271B (en) DC motor servo control method based on robust disturbance compensation scheme

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220516

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20220927

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20221006

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230421

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230522

R150 Certificate of patent or registration of utility model

Ref document number: 7293905

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150