W2 A New Model For Control of Systems With Friction
W2 A New Model For Control of Systems With Friction
Abstract-In this paper we propose a new dynamic model for very much resembles that of a connection with a stiff spring
friction. The model captures most of the friction behavior that has with damper and is sometimes referred to as the Dahl effect.
been observed experimentally. This includes the Stribeck effect, Later studies (see, e.g., [l], [2]) have shown that a friction
hysteresis, spring-like characteristics for stiction, and varying
break-away force. Properties of the model that are relevant to model involving dynamics is necessary to describe the friction
control design are investigated by analysis and simulation. New phenomena accurately.
control strategies, including a friction observer, are explored, and A dynamic model describing the spring-like behavior during
stability results are presented. stiction was proposed by Dahl [3]. The Dahl model is essen-
tially Coulomb friction with a lag in the change of friction
force when the direction of motion is changed. The model has
I. INTRODUCTION many nice features, and it is also well understood theoretically.
~ ?
T 1 1 1 1
1 ' '
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.
, ,,,
420 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 40, NO. 3, MARCH 1995
- 1 7 1 1 -
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.
CANUDAS DE WIT et al.: CONTROL OF SYSTEMS WITH FRICTION 42 1
The derivative is negative when IzI > g(u). Since g ( u ) Fig. 2. Presliding displacement as described by the model. The simulation
is strictly positive and bounded by a, we see that the set was started with zero initial conditions.
I’ z(r)u(~)d.r 2 V ( t )- v ( 0 ) .
experimental results [l]. The stiffness 00 was chosen to give
a presliding displacement of the same magnitude as reported
in various experiments. The value of the damping coefficient
Proofi It follows from (1) that <
01 was chosen to give a damping of = 0.5 for the linearized
equation (9) with a unit mass. The Coulomb friction level FC
corresponds to a friction coefficient p M 0.1 for a unit mass,
and Fs gives a 50% higher friction for very low velocities.
dz The viscous friction u2 and the Stribeck velocity U , are also
2 z-.
dt of the same order of magnitude as given in [I].
Hence The different behaviors shown in the following subsections
cannot be attributed to single parameters but rather to the
L t z ( 7 ) u ( r )d r 2 L t z ( T ) F d r 2 V ( t )- V(0). behavior of the nonlinear differential equation (1) and the
shape of the function g. The presliding displacement and
the varying break-away force are due to the dynamics. This
Linearization in Stiction Regime
behavior is also present in the Dah1 model. The Stribeck shape
To get some insight into the behavior of the model in the of g together with the dynamics give rise to the type of
stiction regime we will consider a mass m in contact with a hysteresis observed in the subsection on frictional lag.
fixed horizontal surface. Let x be the coordinate of the mass,
i.e., U = d x / d t . The equation of motion becomes A. Presliding Displacement
m-
d 2 x = -F = -o0z - g l -d z - d x Courtney-Pratt and Eisner have shown that friction behaves
dt2 dt dt (7) like a spring if the applied force is less than the break-away
where z is given by (1). Linearizing (1) around z = 0 and force. If a force is applied to two surfaces in contact there will
U = 0 we get
be a displacement. A simulation was performed to investigate
if our model captures this phenomenon. An external force was
applied to a unit mass subjected to friction. The applied force
was slowly ramped up to 1.425 N which is 95% of Fs. The
Inserting (8) into (7) gives force was then kept constant for a while and later ramped down
to the value -1.425 N, where it was kept constant and then
d2x dx
m-
dt2
(CTI + +
~ 2 ) -
dt
+
uox = 0. (9) ramped up to 1.425 N again. The results of the simulation
are shown in Fig. 2 where the friction force is shown as a
This shows that the system behaves like a damped second- function of displacement. The behavior shown in Fig. 2 agrees
order system. Notice that the bristle stiffness, CO, is usually qualitatively with the experimental results in [ 131.
7-
~ ~
I’ 1 ‘I 7
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.
422 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 40,NO. 3, MARCH 1995
Y ,
1.2.
Fig. 5. Experimental setup for stick-slip motion.
1Position [m]
1.
1 Velocity [ds]
0 I. io-) 2. 3.G
Fig. 3 . Hysteresis in friction force with varying velocity. The velocity
variation with the highest frequency shows the widest hysteresis loop.
k4
0
13J-Friction force [NI
io m
Time [SI
c
Y Velocity [ds] F
03
Time [SI
0 10 m
Fig. 6 . Simulation of stick-slip motion.
1.2-
did depend on the rate of increase of the force but not on the
0
dwell-time; see also [161. Simulations were performed using
0
1.1 - 0
our model to determine the break-away force for different rates
0
of force application. Since the model is dynamic, a varying
1-
0
0 0 0 0
break-away force can be expected. A force applied to a unit
mass was ramped up at different rates, and the friction force
Force rate [Nk] when the mass started to slide was determined. Note that since
the model behavior in stiction is essentially that of a spring,
there will be microscopic motion, i.e., velocity different from
zero, as soon as a force is applied. The break-away force was
therefore determined at the time where a sharp increase in the
B. Frictional Lag velocity could be observed. Fig. 4 shows the force at break-
Hess and Soom [14] studied the dynamic behavior of fric- away as a function of the rate of increase of the applied force.
tion when velocity is varied during unidirectional motion. They The results agree qualitatively with the experimental results
showed that there is hysteresis in the relation between friction in [151 and [16].
and velocity. The friction force is lower for decreasing veloci-
ties than for increasing velocities. The hysteresis loop becomes [Link]-Slip Motion
wider at higher rates of the velocity changes. Hess and Soom Stick-slip motion is a typical behavior for systems with
explained their experimental results by a pure time delay in friction. It is caused by the fact that friction is larger at rest than
the relation between velocity and friction force. Fig. 3 shows
during motion. A typical experiment that may give stick-slip
a simulation of the Hess-Soom experiment using our friction motion is shown in Fig. 5. A unit mass is attached to a spring
model. The input to the friction model was the velocity which with stiffness k = 2 N/m. The end of the spring is pulled
was changed sinusoidally around an equilibrium. The resulting with constant velocity, i.e., d y l d t = 0.1 m/s. Fig. 6 shows
friction force is given as a function of velocity in Fig. 3. Our results of a simulation of the system based on the friction
model clearly exhibits hysteresis. The width of the hysteresis model in Section 11. The mass is originally at rest and the
loop also increases with frequency. Our model thus captures force from the spring increases linearly. The friction force
the hysteretic behavior of real friction described in [14]. counteracts the spring force, and there is a small displacement.
When the applied force reaches the break-away force, in this
C . Varying Break-Away Force case approximately aog(O), the mass starts to slide and the
The break-away force can be investigated through exper- friction decreases rapidly due to the Stribeck effect. The spring
iments with stick-slip motion. In [15] it is pointed out that contracts, and the spring force decreases. The mass slows
in such experiments the dwell-time when sticking and the down and the friction force increases because of the Stribeck
rate of increase of the applied force are always related and effect and the motion stops. The phenomenon then repeats
hence the effects of these factors cannot be separated. The itself. In Fig. 6 we show the positions of the mass and the
experiment was therefore redesigned so that the time in stiction spring, the friction force and the velocity. Notice the highly
and the rate of increase of the applied force could be varied irregular behavior of the friction force around the region where
independently. The results showed that the break-away force the mass stops.
- 1 I I 1 -1
I -
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.
~
Fig. 7. Block diagram for the servo problem with PID controller
I I
Fig. 9. Block diagram for the position control problem using a friction
observer.
V. FEEDBACK
CONTROL - 14
d i =U--i-ke, k>O
dt dv)
To further illustrate the properties of our friction model we
will investigate its application to some typical servo problems.
First we will use it to show that it predicts limit cycle
oscillations in servos with PID control. We will then use it and the following control law
to design observer based friction compensators.
u = -H(s)e + F + m-d2xd
dt2
A. Limit Cycles Caused by Friction where e = x - X d is the position error and X d is the desired
It has been observed experimentally that friction may give reference which is assumed to be twice [Link] term
rise to limit cycles in servo drives where the controller has ke in the observer is a correction term from the position error.
integral action; for references see [2]. This phenomenon is The closed-loop system is represented by the block diagram
often referred to as hunting. in Fig. 9. With the observer based friction compensation, we
Consider the linear motion of a mass m at position x. The achieve position tracking as shown in the following theorem.
equation of motion is Theorem 1 : Consider system (10) together with the friction
model (1) and (3), friction observer (12) and (13), and control
d2x
m-=u-~ law (14). If H ( s ) is chosen such that
dt2
where d x / d t = U is the velocity, F the friction force given by
(3) and U the control force which is given by the PID controller
is strictly positive real (SPR) then the observer error, F - F,
and the position error, e, will asymptotically go to zero.
Proof: The control law yields the following equations
A block diagram of the system is shown in Fig. 7. In Fig. 8 1 - OlS+(TO
we show the results of a simulation of the system. The friction e= (-F) = (-5) = -G(s)i
parameters are given by Table I, m = 1 and the controller
+
ms2 H ( s ) +
ms2 H ( s )
parameters are K , = 6, K p = 3, and Ki = 4. The reference
position is chosen as X d = 1. The model clearly predicts limit
cycles as have been observed experimentally in systems of where F = F - F and z" = z - 1. Now introduce
this type.
52
V=<TP<+-
B. Friction Compensation k
as a Lyapunov function and
Only linear feedback from the position was used in the PID
control law (1 1). Knowledge about friction was not used. It "
- = A< + B(-5)
is of course more appealing to make a model-based control dt
that uses the model to predict the friction to compensate for e = C<
-r I 1'I 1
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.
424 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 40, NO. 3, MARCH 1995
i
- -W) -
e
is strictly positive real, then the observer error, F - F , and
the velocity error will asymptotically go to zero.
Proofi The theorem is proven in the same way as Theo-
rem 1 after observing that the control law yields the following
error equations
1 -
(-F) =
OlS+(TO
(-5) = - G ( s ) ~
Fig. 10. The block diagram in Fig. 9 redrawn with e and Z as outputs of a
linear and a nonlinear block, respectively.
w - v d =
+
ms H(s) ms H(s) +
which is a state-space representation of G(s). Since G(s) is
SPR [17] it follows from the Kalman-Yakubovitch Lemma This is again an interconnection of a dissipative system with 2
[17] that there exist matrices P = PT > 0 and Q = QT > 0 as its output and a linear SPR system with w - V d as its output.
such that To assume that the friction model and its parameters are
known exactly is of course a strong assumption. Investigation
+
A ~ PP A = -Q of the sensitivity of the results to these assumptions is an
P B = CT. interesting problem that is outside the scope of this paper. The
accuracy required in the velocity measurement is a similar
Now
problem.
VI. CONCLUSIONS
A new dynamic model for friction has been presented. The
model is simple yet captures most friction phenomena that
are of interest for feedback control. The low velocity friction
characteristics are particularly important for high performance
I -IT&<. pointing and tracking. The model can describe arbitrary steady-
state friction characteristics. It supports hysteretic behavior due
The radial unboundedness of V together with the semi- to frictional lag, spring-like behavior in stiction and gives a
definiteness of dV/dt implies that the states are bounded. We varying break-away force depending on the rate of change of
<
can now apply LaSalle’s theorem to see that 4 0 and 2 + 0 the applied force. All these phenomena are unified into a first-
which means that both e and F tends to zero and the theorem order nonlinear differential equation. The model can readily
is proven. be used in simulations of systems with friction.
The theorem can also be understood from the following Some relevant properties of the model have been investi-
observations. By introducing the observer we get a dissipative gated. The model was used to simulate position control of a
map from e to 2 and by adding the friction estimate to servo with a PID controller. The simulations predict hunting
the control signal, the position error will be the output of as has been observed in applications of position control with
a linear system operating on 5. This means that we have integral action. The model has also been used to construct
an interconnection of a dissipative system and a linear SPR a friction observer and to perform friction compensation for
system as seen in Fig. 10. Such a system is known to be position and velocity tracking. When the parameters are known
asymptotically stable. the observer error and the control error will asymptotically go
Velocity Control: The same type of observer-based control to zero. Sensitivity studies, parameter estimation and adapta-
can be used for velocity control. For this control problem the tion are natural extensions of this work.
controller is changed to
.* dud REFERENCES
~ = - H ( s ) e + F + m -dt
[ 11 B. Armstrong-Htlouvry, Control of Machines with Friction. Boston,
MA: Kluwer, 1991.
[2] B. Armstrong-Htlouvry, P. Dupont, and C. Canudas de Wit, “A survey
of models, analysis tools and compensation methods for the control of
machines with friction,” Automatica, vol. 30, no. 7, pp. 1083-1138,
1994.
[3] P. Dahl, “A solid friction model,” Aerospace Corp., El Segundo, CA,
where w - V d is the velocity error and V d the desired velocity Tech. Rep. TOR-0158(3 107-1 8)- 1, 1968.
[4] P.-A. Bliman, “Mathematical study of the Dahl’s friction model,”
which is assumed to be differentiable. Velocity tracking is European J . Mechanics. AISolids, vol. 11, no. 6, pp. 835-848, 1992.
achieved as shown in the following theorem. [5] P.-A. Bliman and M. Sorine, “Friction modelling by hysteresisoperators.
Application to Dahl, Sticktion, and Stribeck effects,” in Proc. Con6
Theorem 2: Consider system (10) together with friction Models of Hysteresis, Trento, Italy, 1991.
model (1) and (3) and observer based control law (15). If [6] C. Walrath, “Adaptive bearing friction compensation based on re-
H ( s ) is chosen such that cent knowledge of dynamic friction,” Automatica, vol. 20, no. 6, pp.
717-727, 1984.
G(s) =
O1S +00 [7] N. Ehrich Leonard and P. Krishnaprasad, “Adaptive friction compen-
ms +H(s) sation for bi-directional low-velocity position tracking,” in Proc. 31st
Conf Decis Contr, 1992, pp. 267-273.
-r n -1- I-
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.
CANUDAS DE WIT ef al.: CONTROL OF SYSTEMS WITH FRICTION 425
[8] J. R. Rice and A.L. Ruina, “Stability of steady frictional slipping,” J . Henrik Olsson (S’91) received the [Link]. degree
Applied Mechanics, vol. SO, no. 2, 1983. in electrical engineering from Lund Institute of
191 P. Dupont, “Avoiding stick-slip through PD control,” lEEE Trans. Technology, Lund, Sweden in 1989.
Automat. Contr., vol. 39, pp. 1094-1097, 1994. He spent the academic year 1989-1990 at the
[IO] D. Haessig and B. Friedland, “On the modeling and simulation of Department of Electrical Engineering at University
friction,” in Proc. I990 Amer. Conrr. Conf., San Diego, CA, 1990, pp. of Califomia, Santa Barbara. Since 1990 he has
1256-1261, been with the Department of Automatic Control at
[I 11 P. Dahl, “Solid friction damping of spacecraft oscillations,” in PI-OC. Lund Institute of Technology where he is currently
AIAA Guidance and Conrr. C o n f , Boston, MA, Paper 75-1 104, 1975. completing the Ph.D. degree. His main research
[ 121 J. Willems, “Dissipative dynamical systems part I: General theory,” interest is in control of nonlinear servosystems.
Arch. Rational Mech. Anal., vol. 45, pp. 321-51, 1972.
[I31 J. Courtney-Pratt and E. Eisner, “The effect of a tangential force on the
contact of metallic bodies,” in Proc. Royal Society, vol. A238, 1957,
pp. 529-550.
[ 141 D. P. Hess and A. Soom, “Friction at a lubricated line contact operating
at ocillating sliding velocities,” J . Tribology, vol. 112, pp. 147-152, Karl Johan Astrom (M’71-SM’77-F’79) receievd
1990. the Ph D. degree in automatic control and mdthemat-
[ 151 V. I. Johannes, M. A. Green, and C. A. Brockley, “The role of the rate ics from the Royal Institue of Technology (KTH),
of application of the tangential force in determining the static friction Stockholm, in 1960
coefficient,” Wear, vol. 24, no. 381-385, 1973. He has been Professor of Automatic Control at
[I61 R. S. H. Richardson and H. Nolle, “Surface friction under time- Lund Institute of Technology/Lund University since
dependent loads,” Wear, vol. 37, no. 1, pp. 87-101, 1976. 1965 His research interests include broad aspect\ of
1171 H. K. Khalil, Nonlinear Systems. New York: Macmillan, 1992. automatic control, stochastic control, system iden-
tification, adaptive control, computer control, and
computer-aided control engineering.
Dr. Astrom has published five books and many
papers. He is a member of the Royal Swedish Academy ot Sciences, and
the Royal Swedish Academy of Engineering Sciences (IVA) He has received
many awards among them the Quazza medal from IFAC in 1987 and the
IEEE Medal of Honor in 1993
Authorized licensed use limited to: TU Delft Library. Downloaded on June 02,2025 at [Link] UTC from IEEE Xplore. Restrictions apply.