Next Article in Journal
Modelling and Simulating the Digital Measuring Twin Based on CMM
Previous Article in Journal
Modeling of Human-Exoskeleton Alignment and Its Effect on the Elbow Flexor and Extensor Muscles during Rehabilitation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Second-Order Dynamic Friction Model Compared to Commercial Stick–Slip Models

Mechanical Engineering, OTH Regensburg, Galgenbergstr. 30, 93053 Regensburg, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Modelling 2023, 4(3), 366-381; https://doi.org/10.3390/modelling4030021
Submission received: 16 June 2023 / Revised: 13 July 2023 / Accepted: 8 August 2023 / Published: 11 August 2023

Abstract

:
Friction has long been an important issue in multibody dynamics. Static friction models apply appropriate regularization techniques to convert the stick inequality and the non-smooth stick–slip transition of Coulomb’s approach into a continuous and smooth function of the sliding velocity. However, a regularized friction force is not able to maintain long-term stick. That is why dynamic friction models were developed in recent decades. The friction force depends herein not only on the sliding velocity but also on internal states. The probably best-known representative, the LuGre friction model, is based on a fictitious bristle but realizes a too-simple approximation. The recently published second-order dynamic friction model describes the dynamics of a fictitious bristle more accurately. It is based on a regularized friction force characteristic, which is continuous and smooth but can maintain long-term stick due to an appropriate shift in the regularization. Its performance is compared here to stick–slip friction models, developed and launched not long ago by commercial multibody software packages. The results obtained by a virtual friction test-bench and by a more practical festoon cable system are very promising. Thus, the second-order dynamic friction model may serve not only as an alternative to the LuGre model but also to commercial stick–slip models.

1. Introduction

Friction has long been an important issue in multibody dynamics. A detailed survey and comparison of several friction force models for dynamic analysis of multibody mechanical systems is provided in [1]. Just as in the present work, the use of internal states defines the difference between dynamic and static friction models. It starts with different variations of the static Coulomb model, including viscous parts and the Stribeck effect, discusses several regularization approaches with finite slopes at zero velocity, and ends up with dynamic friction force models, including the probably best-known representative, the LuGre friction model.
Static friction models, where an appropriate regularization technique converts the stick inequality and the non-smooth stick–slip transition of Coulomb’s approach into a continuous and smooth function of the sliding velocity, are not able to maintain long-term stick. That is why commercial multibody software packages, like Adams, RecurDyn, and Simpack, also offer dynamic friction models. To reproduce stick–slip effects and maintain long-term stick, they rely on their own developments and do not simply use models known from the literature, e.g., [1]. Just the LuGre model, as published in [2], is implemented in the commercial software package Adams [3]. However, as demonstrated in [4], the LuGre approach does not appear to be an engineer’s first choice because it has too many drawbacks and is not able to reproduce a pre-defined friction characteristic in dynamic applications.
A second-order dynamic friction model (FrDyn2) was introduced in [4] as a reference to the LuGre and to a standard static friction model. The FrDyn2 model produced accurate and reliable results in standard stick–slip examples as well as in a more practical model of a festoon cable system.
The present paper compares this model to commercial stick–slip models as as presented and analyzed in [5]. It turns out that the concepts of Adams and Recurdyn are rather similar but different from the Simpack stick–slip model. That is why just the Adams and the Simpack stick–slip models are used here for the comparison with the FrDyn2 model.
Section 2 illustrates a typical standard and the shifted regularization technique and provides a short description of the FrDyn2 model, as defined in [4]. Section 3 demonstrates the long-term stick potential of the FrDyn2 model, achieved by an appropriate shift in the standard regularization. It also illustrates the poor performance of a static friction model which approximates stick by slow sliding. A virtual friction test-bench was set up in Matlab for that purpose. Section 4 compares the break-away and stick–slip transition of the FrDyn2 model with the stick–slip models of Adams and Simpack. For that purpose, the virtual friction test-bench was also set up in Adams and Simpack. The focus of Section 5 is on the dynamic response to step-like force excitations with different durations and a force amplitude close to the adhesion limit. Section 6 analyzes the performance of the friction models under consideration by the more practical model of a festoon cable system, with a setup in Matlab, Adams, and Simpack. Section 7 summarizes and discusses the results and and provides an outlook for future research.
Finally, Appendix A contains a Matlab script and the corresponding Matlab functions describing the virtual friction test-bench operated with the FrDyn2 model.

2. Static and Dynamic Friction Models

The idealized friction model of Coulomb simply distinguishes between sticking and sliding, as seen in Figure 1a. The friction force F R depends on the sliding velocity v and is realized by combining an inequality with a simple relation
F R μ s F N if v = 0 and F R = v v μ d F N if v > 0
The friction force is proportional to the normal force F N and characterized by the parameters μ s and μ d , which specify the static and the dynamic coefficients of friction.
Coulomb’s dry friction approach, as defined in Equation (1), is practically unusable in general multibody dynamics due to the inequality representing stick. That is why dry friction is usually approximated by an unambiguous function, as seen in Figure 1b. The regularized characteristic F R = F R ( v ) introduces the fictitious velocity v s , which defines the width of the regularization interval v s v + v s and the dynamic velocity v d that characterizes at v v d the range of full sliding, where the simple relation F R = ± μ d F N applies. The transitions from the value pairs ( v d , μ d F N ) ( v s , μ s F N ) and ( v s , μ s F N ) ( + v s , + μ s F N ) as well as ( + v s , + μ s F N ) ( + v s , + μ s F N ) , which are usually modeled by sufficiently smooth functions, like polynomials or trigonometric functions.
Static friction models are widely used in multibody dynamics and control theory. They are able to reproduce stick–slip effects in standard and in a more practical application [4]. Even much simpler regularizations, assuming by μ = μ s = μ d just one unique friction value, are in good conformity to dynamic measurements [6]. However, static friction models cannot maintain long-term stick. They describe the friction force just as a function of the sliding velocity F R = F R ( v ) . Figure 1b defines a commonly used regularized friction characteristics F R = F R ( v ) by a set of “static” ( μ s , v s ) and “dynamic” ( μ d , v d ) friction parameters. However, the use of dynamic friction parameters is no essential criteria for a dynamic friction model, as erroneously assumed in [7]. Here, a static friction force defined by F f = F f ( μ s , v s , μ d , F n ) making use of the dynamic friction value μ d is already supposed to define a dynamic friction force.
Dynamic friction models are characterized by the use of internal states. The friction force is then defined by a more complex function F R = F R ( v , s ) , where s collects the internal states. Dynamic friction models are able to reproduce stick–slip effects and maintain long-term stick [8].
A software-capable formulation of a friction force model also takes the friction model parameters into account. Then, F R = F R ( v , p ) characterizes a static and F R = F R ( v , s , p ) a dynamic friction model, where p collects the friction model parameters.
The well-known LuGre friction model uses the displacement z of a fictitious bristle as an internal state. However, as demonstrated in [4], the LuGre approach represents just a first step approximation to the dynamics of a massless bristle, which results in several drawbacks of this dynamic friction model. If the mass of the fictitious bristle is also taken into account, this results in a second-order dynamic friction model (FrDyn2), which was introduced in [4] as a reference. The FrDyn2 model exhibits none of the LuGre drawbacks and performed well in standard stick–slip examples and in a more practical model of a festoon cable system.
The second-order bristle dynamics is defined by
m b z ¨ = F R F B = F R ( v C z ˙ ) σ 0 z + σ 1 z ˙
where m b z ¨ approximates the inertia force of the fictitious bristle, F R = F R ( v C z ˙ ) describes the friction force, F B = σ 0 z + σ 1 z ˙ models the bristle as a visco-elastic element, and v C represents the component of the contact point velocity that is perpendicular to the contact normal. Just as with the LuGre model, σ 0 and σ 1 describe the stiffness and the damping of the fictitious bristle. The fictitious mass of the bristle is defined by
m b = σ 1 2 4 σ 0
which represents the aperiodic case of the homogenous second-order differential Equation (2), thus avoiding unwanted oscillations in the fictitious bristle. The FrDyn2 model implies with s = z , z ˙ two internal states. It is provided as a Matlab function in Appendix A by Listing A4.
The stick–slip models of Adams and RecurDyn describe the dynamic friction force just as a function of the contact point velocity F R d = F N μ d ( v C ) and approximate the static friction force by a two-dimensional function F R s = F N μ s ( v C , x ) , where x is a fictitious displacement [5]. A smooth function of the contact point velocity v C , which is not explained in detail in the user manual, describes the transition from μ s to μ d . The fictitious displacement serves as an internal state and generates the static friction force F R s or the static friction coefficient μ s as a nonlinear function of x by adding a viscous damping term.
Simpack provides a stick–slip model which realizes, like Coulomb’s approach in Figure 1a, a sudden drop from the static μ s to the dynamic friction coefficient μ d [5]. In the adhesion region, the friction force is approximated by a visco-elastic element whose deflection again represents an internal state of this stick–slip model.

3. Long-Term Stick of the Second-Order Dynamic Friction Model

By applying an appropriate horizontal shift to the regularized friction characteristics, as indicated in Figure 1c, the second-order dynamic friction model can maintain long-term stick. The steady state solution ( z ¨ = 0 and z ˙ = v C ) of the fictitious bristle dynamics (2) provides the required sticking force as
F 0 = F R ( 0 ) = σ 0 z + σ 1 v C
The FrDyn2 model describes here the transition from the static to the dynamic friction force ( μ s , F s ) ( μ d , F d ) by a cubic polynomial and defines the friction force in the regularization range v s v + v s by a parabolic function. Then,
Δ v = F 0 | F 0 | v s w F with w F = 1 1 | F 0 | / F s if | F 0 | F s 1 elsewhere
delivers the corresponding horizontal shift. The static friction force is defined by F s = F N μ s and F 0 / | F 0 | adjusts the horizontal shift to the sign of the required sticking force value.
Figure 2 provides a virtual test-bench which is used here to demonstrate the long-term sticking quality of the FrDyn2 model and in the following for a comparison to the stick–slip models of Adams and Simpack.
A body of unit mass m = 1 kg is in contact to a horizontal rough plate, as seen in Figure 2a. It is exposed to a horizontal force F = F ( t ) , which is continously increased in the interval 0 t t E from F = 0 to the final value F = F E , as defined by Figure 2d. The friction parameters μ s , v s and μ d , v d , provided by Figure 2c, model a regularized friction characteristic as defined in Figure 1b. The FrDyn2 model uses the states z and z ˙ of a fictitious bristle, where m b denotes its fictitious mass and σ 0 and σ 1 characterize the visco-elastic properties of the bristle, as seen in Figure 2b. The bristle parameters, as defined by Figure 2c, are adjusted to the body and the friction parameters by estimating a reference friction force of F R r e f 5 N and defining a reference bristle deflection of z r e f = 1 · 10 6 m . Then, σ 0 = F R r e f / z r e f = 5 · 10 6 N / m provides the bristle stiffness. The reference friction force corresponds to a reference mass of m r e f = F R r e f / g = 0.51 kg . It provides in a first guess the damping of the fictitious bristle as σ 1 = 2 m r e f σ 0 = 3192.75 N / ( m / s ) , when the aperiodic case of a fictitious oscillator consisting of m r e f , σ 1 , and σ 0 is assumed hereby.
The Matlab script, provided in Appendix A by Listing A1, performs simulations with different step durations. It applies the standard implicit solver ode15s where the default tolerances are changed to RelTol = 1 · 10 6 and AbsTol = 1 · 10 9 . The reference bristle deflection was chosen very small in this example. The Matlab function dyn_fr_test_bench, provided in Appendix A by Listing A2, computes the dynamics of the virtual friction test-bench including the FrDyn2 model as a set of first-order differential equations. The Matlab functions Step3 and FrDyn2, defined in Appendix A by Listings A3 and A4, provide the step input and the dynamics of the FrDyn2 model. The simulation results, plotted in Figure 3, demonstrate that the FrDyn2 model perfectly maintains long-term stick.
The excitation force is slowly increased in the time interval 0 s t 30 s from F = 0 to 99% of the static friction force F s , which is given here by F s = μ s F N = 0.6 · 1 kg · 9.81 m / s 2 = 5.886 N . In the subsequent time interval 30 s t 300 s , the force F E = 0.99 · 5.886 N = 5.827 N is kept constant. The friction force F R ( t ) , generated by the FrDyn2 model, counteracts the excitation force F R = F , as seen in the upper-left plot in Figure 3. It is worth noting that the friction force, as defined by the free-body diagram in Figure 3b, points in the opposite direction of the contact point velocity. The tip of the bristle sticks to the ground because the excitation force does not exceed the static friction force. Then, the bristle deflection z coincides with the body displacement u and provides the friction force as a function of the bristle deflection z = u , as demonstrated by the lower-right plot in Figure 3. As a consequence, the body shifts slightly and comes to a stand-still at the steady state value of z s t = u s t = F E / σ 0 = 5.827 N / 5.0 · 10 6 N / m = 1.1654 · 10 6 m . The friction force diagram F R ( v ) , in the upper-right plot of Figure 3, shows that the FrDyn2 model can reproduce the ambiguous part of Coulomb’s friction law at vanishing contact point velocities v = 0 .
A static friction model describes the friction force just as a function of the sliding velocity, F R = F R ( v ) . A typical regularization without any horizontal shift is illustrated in Figure 1. A static friction model can reproduce the required friction force F R = F but definitely cannot maintain stick for longer time intervals, as seen in Figure 4.
In the regularization range, the friction force characteristic F R = F R ( v ) is described by a parabolic function. This pre-defined function, plotted in the right graph of Figure 4 by a thin dashed line, perfectly coincides with the computed friction force F R plotted by a solid thick line. The parabola delivers the required steady state friction force of F R = F s t = 0.99 · F s at the velocity of v s t = 0.9 · v s = 0.9 · 10 3 m / s . As a consequence, the body does not come to rest but continues to move inexorably at this velocity.

4. Break-Away and Stick–Slip Transition

The virtual friction test-bench, defined in Figure 2, consists of a body in contact to a horizontal rough plate. The excitation force F = F ( t ) is now slowly increased within 10 s from F = 0 to the final value F = F E = 1.05 F s = 6.180 N , which exceeds the adhesion limit of F s by 5%. The excitation force F = F ( t ) , modeled by a third-order polynomial, reaches the adhesion limit F = F s = 5.886 N at time t = t B = 8.6808 s .
The simulation results are plotted in Figure 5. The dashed blue line, the solid black line, and the dotted red line represent the results computed with the FrDyn2 model, the Adams stick–slip model, and the Simpack stick–slip model.
At first ( 0 t t B ), the body remains in a quasi-static equilibrium, where the slowly increasing excitation force F ( t ) is perfectly counteracted by a friction force F R ( t ) = F ( t ) generated by each of the friction models, as demonstrated in Figure 5a, where the lines for F ( t ) and F R ( t ) perfectly coincide in the time interval 0 t t B . In general, the friction models under consideration generate friction forces depending not only on the contact velocity but also on internal states. The FrDyn2 model uses the displacement z of a fictitious bristle and its time derivative z ˙ as internal states s = z , z ˙ . In the quasi-static equilibrium mode, the velocity of the body and the time derivative of the bristle deflection are negligibly small v 0 and z ˙ 0 . In this mode, the tip of the fictitious bristle sticks to the ground, which results in a bristle deflection that equals the body displacement z = u . The compliance of the fictitious bristle is modeled by a viscous force element, F B = σ 0 z + σ 1 z ˙ . In the quasi-static equilibrium mode, the body acceleration a = u ¨ is negligibly small too, as indicated in Figure 5b by the time history a = a ( t ) in the time interval 0 t t B . According to Equation (2), the friction force generated by the FrDyn2 model corresponds then to the elastic part of the bristle force, F R σ 0 z = σ 0 u . The quasi-static force F R = σ 0 u equals the static friction force F s at the reference displacement of
u = x R = F s σ 0 = 5.886 N 5 · 10 6 N / m = 1.1772 · 10 6 m
The adhesion range is characterized by a vanishing sliding velocity ( v 0 ) and extends here to displacements in the range of 0 u x R . In this range, the friction force is generated as a function of the displacement, where the FrDyn2 model and the stick–slip model of Simpack apply a linear and the stick–slip model of Adams a nonlinear digressive function, as seen in Figure 5f. The Adams manual does not specify the type of nonlinearity but, as indicated by the solid black line in Figure 5f, it approaches the limit value F R = F s at the reference displacement x R defined in (6) with a vanishing inclination.
The stick–slip model of Simpack is based on Coulomb’s approach, where the friction force drops in an instant from the static to the dynamic value as soon as the excitation force exceeds the static friction force at t = t B , as seen in Figure 5a,d in particular. The transition from the static to the dynamic friction force F s F d are modeled in the FrDyn2 and the Adams stick–slip models as functions of the velocity v controlled by the parameters v s and v d . The FrDyn2 model applies a cubic polynomial which is shifted in the horizontal direction to maintain stick at v = 0 , as indicated in Figure 1. As can be seen by inspecting Figure 5e, the FrDyn2 model generates a friction force characteristic (dashed blue line) which reproduces the pre-defined and horizontally shifted one (magenta colored circles) nearly perfectly. The friction characteristics produced by the Adams stick–slip model is rather similar (solid black line). Most likely, Adams models the transition F s F d by a fifth-order polynomial. As a consequence, the FrDyn2 and the Adams stick–slip models produce slightly more delayed drops in the time histories of the computed friction forces, as seen in Figure 5d.
Figure 5b,c illustrate the break-away effect at t = t B by the time histories of the body acceleration u ¨ = a = a ( t ) and the body displacement u = u ( t ) . All friction models under consideration approximate sliding at v v d by a constant friction force F R ( v v d ) = F d = 2.943 N . Viscous components in the friction force are not considered here. The free body diagram in Figure 2b delivers the linear momentum
m u ¨ = F F R
for the body of mass m = 1 kg . At t > 10 s which includes t > t B , the applied force is defined by F = F E = 1.05 F s = 6.180 N and the friction force is represented by its dynamic value F R = F d = 2.943 N . Then, the maximum acceleration of the body is defined by
a m a x = u ¨ ( t > 10 s ) = ( F E F d ) / m = ( 6.180 N 2.943 N ) / 1 kg = 3.237 m / s 2
which is exactly reproduced by the friction models, as seen in Figure 5b.

5. Dynamic Response

A pulse load excitation, performed in [4], revealed the tendency of dynamic friction models to produce dynamic overshoots in the friction force time histories. That is why the virtual test-bench defined in Figure 2 is now exposed to excitation forces where the amplitude F E = 0.95 F s = 5.592 N is 5% less than the static friction force F s = 5.886 N and the step duration is varied from t E = 0.1 s to t E = 0.0001 s . The corresponding simulation results are shown in Figure 6. The solid thin gray line represents the excitation force F = F ( t ) , the dashed blue line, the solid black line, and the dotted red line mark the friction forces F R = F R ( t ) computed with the FrDyn2, the Adams stick–slip, and the Simpack stick–slip models.
The time histories of the friction forces perfectly coincide with the excitation force F R ( t ) = F ( t ) at a step duration of t E = 0.1 s , as seen in Figure 6a. All friction models operate here in a quasi-static sticking mode where the friction forces are practically generated as a function of the body displacement, as already illustrated in Figure 5e, and in this specific case by Figure 7.
The forces F R generated by the friction models depend here practically not on the velocity v = u ˙ but only on the displacement of the body u, as seen in Figure 7a,b. In case of the FrDyn2 and the Simpack stick–slip models F R = σ 0 u holds, which provides the friction force F R = F E = 5.592 N at the steady state displacement u = x s t D 2 = x s t S = 5.592 N / 5 · 10 6 N / m = 1.12 · 10 6 m , as indicated in Figure 7b,c by thin dashed black lines. Adams models the friction force a quasi-static sticking mode by a strongly nonlinear and degressive function of the displacement. The Adams manual does not specify this function but the simulation results provide the friction force F R = F E = 5.592 N at the steady state displacement u = x s t A = 0.71 · 10 6 m , as seen in Figure 7b,c. As expected from the time histories F R = F R ( t ) plotted in Figure 6a, the time histories of the body displacement reach their steady state values u = x s t D 2 = x s t S and u = x s t A at t > t E without any overshoots, as seen in Figure 7c.
In a quasi-static mode, the tip of the fictitious bristle, which forms the basis of the FrDyn2 model, sticks to the ground. Then, the linear momentum (7) of the body in the virtual friction test-bench simplifies to
m u ¨ = F ( t ) σ 0 u or m u ¨ + σ 0 u = F ( t )
where the quasi-static friction force is generated by the bristle compliance F R = σ 0 z and z = u holds in addition. The simplified equation of motion (9) is characterized by the eigen-frequency ω 0 = σ 0 / m and delivers the value and its corresponding oscillation period as
ω 0 = 5 · 10 6 N / m / 1 kg = 2236 s 1 and T = 2 π ω 0 = 0.003 s
As a consequence, even a rather short step duration of t E = 0.01 s will still represent a subcritical excitation of the virtual friction test-bench. The time histories of the friction forces F R ( t ) exhibit just a small overshoot at t > t E = 0.01 s , as seen in Figure 6b.
The situation becomes complicated for step durations t E < T , as seen in Figure 6c,d. The Simpack stick–slip model (dotted red line) generates now significant overshoots, which amount to
F R m a x S i m p a c k t E = 0.001 s = 6.61 N and F R m a x S i m p a c k t E = 0.0001 s = 6.75 N
The values exceed the steady state value F R s t = F E = 5.592 N by 18.2 % and 20.7 % and even the static friction value F s = 5.886 N by 12.3 % and 14.5 %, which calls into question the physical basis of this stick–slip model. The time histories of the friction forces generated by the FrDyn2 and the Adams stick–slip models (dashed blue and solid black lines) differ somehow. But both models limit the friction force to the static value F R F s , as expected from friction models in general.
The friction models now generate friction forces which strongly depend on the body velocity v = u ˙ and the body displacement u, as seen in Figure 8a,b.
The Simpack stick–slip model (dotted red line) generates a time history of the body displacement which approaches the quasi steady state value u = x s t = 1.12 · 10 6 m with a small overshoot shortly after the step duration of t E = 0.001 s , as seen in Figure 8c. The FrDyn2 model overshoots and partly slides, resulting in displacements at t > t E which exceed with u s t D 2 = 2.27 · 10 6 m the quasi steady state value of x s t = 1.12 · 10 6 m significantly. Hence, the FrDyn2 model generates a dynamic break-away effect at high frequent excitation loads, which are close (here, 95%) to the static friction force. The time history of the body displacement u = u ( t ) shows a strange behavior for the Adams stick–slip model, as exemplified by the solid black line in Figure 8c. At first ( t < 3.5 · 10 3 s ), it approaches the quasi steady state value u = x s t D 2 = x s t S = 1.12 · 10 6 m of the FrDyn2 and Simpack solution and then ( 4 · 10 3 s < t 10 · 10 3 s ) it starts to decrease very slowly, but the simulated time interval 0 t 10 · 10 3 s is too short to indicate a limit value.
A simulation with the Adams stick–slip model over a longer time period results in the time history of the body displacement u = u ( t ) , as plotted in Figure 9.
The section shown in Figure 8c has a light gray background. It seems that Adams applies in its stick–slip model different time constants for the increase and decrease in the body displacement u = u ( t ) . The force excitation with a step duration of t E = 0.001 s is much faster than the dynamics of the virtual friction test-bench. The body displacement reaches its maximum value at t 0.0035 s , which corresponds to the oscillation period computed in (10). The decay from the maximum displacement to the steady state value takes about 0.05 s , which is fourteen times as much. This strange behavior was also reported in [5], wherein pulse loads are applied to a single mass resting on a horizontal plate. At the end of a series of impulse loads each of magnitude 0.8 F s , the body is returned to its initial position. However, in the present example, a small but permanent deviation of δ u = u s t A x s t A = 0.12 · 10 6 m remains, as seen in Figure 9. This indicates that the Adams stick–slip model also tends to partly slip, when high frequent excitation loads close to the static friction force are applied.

6. The Festoon Cable System Model

A planar model of a festoon cable system is used in [4] to asses different friction models in a more practical example. The model consists of three cable m C and two trolley m T masses, as seen in Figure 10.
At the beginning ( t = 0 ) , the towing trolley is fixed at u T T ( 0 ) = 1.5529 m . The equilibrium position of the cable system places the trolleys at x T 1 ( 0 ) = u T T ( 0 ) / 3 , x T 2 ( 0 ) = 2 u T T ( 0 ) / 3 , and locates the cable masses at x C 1 ( 0 ) = u T T ( 0 ) / 6 , x C 2 ( 0 ) = u T T ( 0 ) / 2 , x C 3 ( 0 ) = 5 u T T ( 0 ) / 6 , as well as z C 1 ( 0 ) = z C 2 ( 0 ) = z C 3 ( 0 ) = 0.96698 m . The non-holonomic constraint u ˙ T T = v T T relates the movements of the towing trolley to a pre-defined velocity profile v T T = v T T ( t ) . The velocity profile, defined by the solid gray lines in Figure 11a, models an extension maneuver, which moves the towing trolley from the initial position u ( t = 0 ) = u T T ( 0 ) = 1.5529 m to a final position of u T T ( t > 5 s ) = 5.553 m .
The Matlab simulation with the FrDyn2 model generates the output at every simulation step. It applies the Matlab standard solver for stiff differential equations ode15s with error tolerances of RelTol = 1 · 10 6 and AbsTol = 1 · 10 9 .
The Adams and Simpack simulations were performed with an output step size of Δ t = 1 · 10 4 s . The dashed blue, the solid black, and the dotted red lines mark the results obtained by the FrDyn2, the Adams, and the Simpack stick–slip models.
The movement of the towing trolley ends at t = 5 s . After that, the trolleys perform to and fro motions which at t 5 s are indicated by sign changes in the time histories of trolley velocities v T 1 ( t ) and v T 2 ( t ) . An arrow pointing from v T 1 ( t 6.47 s ) = 0 over F N 1 down to F R 1 highlights such an event, in particular. The dynamic motions of the cable masses induce variations in the normal forces F N 1 and F N 2 acting between the trolleys and the rail, as seen in Figure 11b. The time histories of the velocities v T 1 ( t ) , v T 2 ( t ) and the normal forces F N 1 ( t ) , F N 2 ( t ) generated with FrDyn2 and the stick–slip models of Adams and Simpack match nearly perfectly. However, the time histories of the friction forces F R 1 ( t ) and F R 2 ( t ) exhibit some discrepancies, as seen in Figure 11c. In particular, when the trolley velocities change their signs or during a sticking period of trolley 2.
The plots in Figure 12 focus on a sign change in the trolley velocity v T 1 at t 6.47 s and a sticking period of trolley 2 in the time interval 7.2 s t 7.8 s .
In the very short time interval 6.45 s t 6.48 s , the normal force between trolley 1 and the rail amounts to F N 1 = 21.2 N , as indicated in Figure 11b. The friction values defined in Figure 10 provide in this case a static friction force of F s = μ s F N 1 = 0.08 · 21.2 = 1.696 N and a dynamic friction force of F d = μ d F N 1 = 0.05 · 21.2 = 1.06 N . At times t = 6.45 s and t = 6.48 s , the first trolley is in a full sliding mode, as indicated in Figure 12a by the friction forces F R 1 ( t = 6.45 s ) = + F d = 1.06 N and F R 1 ( t = 6.48 s ) = F d = 1.06 N . These sliding modes are perfectly reproduced by the friction models under consideration. Shortly before the sign change in the trolley velocity, the friction forces computed by the FrDyn2 and the Adams stick–slip model make use of the Stribeck effect, which models a velocity-dependent transition from the static to the dynamic friction force and vice versa. The Simpack stick–slip model approaches v T 1 0 with the dynamic force value and does not reproduce a potential velocity-dependent increase in the friction force. The FrDyn2 and the Simpack stick–slip models describe the friction force at v = 0 by a linear spring, whereas the Adams stick–slip model uses a nonlinear approach. That is why, the FrDyn2 model (dashed blue line) corresponds in the time interval 6.46 s t 6.47 s more to the Simpack (dotted red line) than to the Adams stick–slip model (solid black line). The friction forces of the FrDyn2 and the Adams stick–slip models are limited to the static value F R 1 F s , which results in F R 1 D 2 = F R 1 A = 1.696 N at t 6.47 s . However, The Simpack stick–slip model overshoots and produces the peak value of F R 1 S = 2.685 N , which exceeds by nearly 60% the static friction force F s or F S , respectively.
The sticking period 7.2 s t 7.8 s is represented rather similarly by the friction models under consideration, as seen in Figure 12b. Again, the FrDyn2 and the Adams stick–slip models increase the friction forces from the dynamic to the static value when approaching stand-still at t 7.2 s . However, the small time delay of the peak values visible in Figure 12a is not noticeable due to the large time interval applied in this plot. The Simpack stick–slip model is based on Coulomb’s approach, which results in the discontinuities at the slip–stick and stick–slip transitions in the dotted red line.

7. Discussion

The present manuscript compares a second-order dynamic friction (FrDyn2) model with the commercial stick–slip models of Adams and Simpack. The comparison is performed here with a virtual friction test-bench and a more practical model of a festoon cable system.
All models can maintain long-term stick. The FrDyn2 model corresponds partly to the Adams and partly to the Simpack stick–slip models. The FrDyn2 and the Adams stick–slip models show dynamic break-away effects at high frequent excitation loads, which are close (here 95%) to the static friction force. The Simpack stick–slip model avoids dynamic break-away effects by overshoots in the friction force that far exceed the static friction force. Adams models the decay of a friction force overshoot much slower than the increase, whereas the FrDyn2 model considers in both cases the dynamics of the fictitious bristle.
The FrDyn2 model is based on a fictitious bristle characterized by its mass, stiffness, and damping. The fictitious mass of the bristle is automatically adjusted to the stiffness and damping parameters. The pre-defined friction force characteristic F R = F R ( v ) is described here by piecewise-defined polynomials but not limited to this. The bristle parameters can easily be derived from estimated reference friction forces and estimated bristle deflections.
The results obtained by the FrDyn2 model are reliable and based on the physical nature of the friction model approach, which makes the second order dynamic friction model a suitable alternative to commercial stick–slip models.
The dynamics of the FrDyn2 model are governed by the friction characteristics. A rapid transition from the static to the dynamic force, modeling the Stribeck effect, results in a stiff performance of the FrDyn2 model. But the shifted regularized friction characteristic is completely continuous and smooth, which makes it possible to apply any standard stiff ode solver for a multibody system which incorporates the FrDyn model.
Future works will implement the FrDyn2 model as an external force element in Adams and Simpack. Then, the run-time performance of this model can also be compared to the corresponding commercial stick–slip models. The influence of the velocities v s and v d , which model the regularization and the Stribeck effect on the results and on the run-time performance, will be studied in addition.

Author Contributions

Both authors validated the results, wrote, and reviewed the manuscript. G.R. developed the second-order dynamic friction model, performed the Matlab Simulations, and produced the figures. M.S. performed the simulations with the Adams and the Simpack stick–slip models. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Appendix A provides a Matlab script and functions which realize the simple friction test-bench including the second-order dynamic friction model.

Acknowledgments

The authors would like to thank the “Walter and Elfriede Niehoff Stiftung” for funding a research assistant position for Matthias Schuderer at the OTH Regensburg.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Virtual Friction Test-Bench Realized in Matlab

Listing A1. Virtual friction test-bench exposed to step-like force inputs.
Modelling 04 00021 i001Modelling 04 00021 i002
Listing A2. Virtual friction test-bench exposed to step-like force inputs.
Modelling 04 00021 i003Modelling 04 00021 i004
Listing A3. Virtual friction test-bench exposed to step-like force inputs.
Modelling 04 00021 i005
Listing A4. Virtual friction test-bench exposed to step-like force inputs.
Modelling 04 00021 i006Modelling 04 00021 i007

References

  1. Marques, F.; Paulo Flores, P.; Pimenta Claro, J.C.; Lankarani, H.M. A survey and comparison of several friction force models for dynamic analysis of multibody mechanical systems. Nonlinear Dyn. 2016, 86, 1407–1443. [Google Scholar] [CrossRef]
  2. Åström, K.J.; Canudas de Wit, C. Revisiting the LuGre friction model. IEEE Control Syst. Mag. 2008, 28, 101–114. Available online: https://hal.science/hal-00394988/document (accessed on 10 July 2023).
  3. Adams; Joint Friction. Available online: https://help.hexagonmi.com/de-DE/bundle/adams_2021.4/page/adams_help/Adams_Basic_Package/view/building_models_joints/TOC.Idealized.Joints.xhtml (accessed on 10 July 2023).
  4. Rill, G.; Schaeffer, T.; Schuderer, M. LuGre or not LuGre. Multibody Syst. Dyn. 2023, 1–28. [Google Scholar] [CrossRef]
  5. Schuderer, M.; Rill, G.; Schaeffer, T.; Schulz, C. Friction Modeling from a Practical Point of View. In Proceedings of the ECCOMAS Thematic Conference on Multibody Dynamics, Lisbon, Portugal, 24–28 July 2023; Available online: https://multibody2023.tecnico.ulisboa.pt/prog_MULTIBODY_WEB/MULTIBODY2023_PAPERS/ID_104_516_full_paper_Friction_Modeling.pdf (accessed on 8 August 2023).
  6. Pires, I.; Ayala, H.; Weber, H. Ensemble Models for Identification of Nonlinear Systems with Stick-Slip. In Proceedings of the ENOC 2020+2, Lyon, France, 17–22 July 2022; Available online: https://enoc2020.sciencesconf.org/386539/document (accessed on 10 June 2023).
  7. Chaturvedi, E.; Mukherjee, J.; Sandu, C. A novel dynamic dry friction model for applications in mechanical dynamical systems. Inst. Mech. Eng. Part K J. Multibody Dyn. 2023, 14644193231169325. [Google Scholar]
  8. Jing, Q.; Mi, N. Investigation of Selection Mechanism of Friction Models in Multibody Systems. In Proceedings of the 5th International Conference on Vehicle, Mechanical and Electrical Engineering (ICVMEE), Dalian, China, 28–30 September 2019; pp. 251–260. Available online: https://www.scitepress.org/Papers/2019/88736/88736.pdf (accessed on 10 June 2023).
Figure 1. Dry friction: (a) Coulomb’s approach, (b) regularized approximation, (c) shifted regularization.
Figure 1. Dry friction: (a) Coulomb’s approach, (b) regularized approximation, (c) shifted regularization.
Modelling 04 00021 g001
Figure 2. Virtual friction test-bench.
Figure 2. Virtual friction test-bench.
Modelling 04 00021 g002
Figure 3. Long-term stick potential of second-order dynamic friction model.
Figure 3. Long-term stick potential of second-order dynamic friction model.
Modelling 04 00021 g003
Figure 4. Slow sliding approximates sticking at a standard static friction model.
Figure 4. Slow sliding approximates sticking at a standard static friction model.
Modelling 04 00021 g004
Figure 5. Virtual friction test-bench simulation results with a force that exceeds the adhesion limit.
Figure 5. Virtual friction test-bench simulation results with a force that exceeds the adhesion limit.
Modelling 04 00021 g005
Figure 6. Dynamic friction forces resulting from step force excitation with different durations.
Figure 6. Dynamic friction forces resulting from step force excitation with different durations.
Modelling 04 00021 g006
Figure 7. Friction force characteristics and body displacements at a step duration of t E = 0.1 s .
Figure 7. Friction force characteristics and body displacements at a step duration of t E = 0.1 s .
Modelling 04 00021 g007
Figure 8. Friction force characteristics and body displacements at a step duration of t E = 0.001 s .
Figure 8. Friction force characteristics and body displacements at a step duration of t E = 0.001 s .
Modelling 04 00021 g008
Figure 9. Body displacements generated with Adams at a step duration of t E = 0.001 s .
Figure 9. Body displacements generated with Adams at a step duration of t E = 0.001 s .
Modelling 04 00021 g009
Figure 10. Multibody model of a crane festoon system as defined in [4].
Figure 10. Multibody model of a crane festoon system as defined in [4].
Modelling 04 00021 g010
Figure 11. Results of a simulated festoon extension maneuver.
Figure 11. Results of a simulated festoon extension maneuver.
Modelling 04 00021 g011
Figure 12. Friction forces in specific time intervals.
Figure 12. Friction forces in specific time intervals.
Modelling 04 00021 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rill, G.; Schuderer, M. A Second-Order Dynamic Friction Model Compared to Commercial Stick–Slip Models. Modelling 2023, 4, 366-381. https://doi.org/10.3390/modelling4030021

AMA Style

Rill G, Schuderer M. A Second-Order Dynamic Friction Model Compared to Commercial Stick–Slip Models. Modelling. 2023; 4(3):366-381. https://doi.org/10.3390/modelling4030021

Chicago/Turabian Style

Rill, Georg, and Matthias Schuderer. 2023. "A Second-Order Dynamic Friction Model Compared to Commercial Stick–Slip Models" Modelling 4, no. 3: 366-381. https://doi.org/10.3390/modelling4030021

APA Style

Rill, G., & Schuderer, M. (2023). A Second-Order Dynamic Friction Model Compared to Commercial Stick–Slip Models. Modelling, 4(3), 366-381. https://doi.org/10.3390/modelling4030021

Article Metrics

Back to TopTop