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

1. Introduction

Friction has long been an important issue in multibody dynamics. Today, static and dynamic friction force models are applied in multibody dynamics [1]. Static friction models usually 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 the last 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 [2]. A second order dynamic friction model describes the dynamics of a fictitious bristle more accurately. It was introduced in [2] as a reference and performed well in standard stick-slip examples as well as in a more practical model of a festoon cable system. In the present paper the second order dynamic friction model is compared to stick-slip friction models, developed and launched not long ago by commercial multbody software packages. They are presented and analyzed in [3]. It turns out that the concepts of ADAMS and Recurdyn are quite similar but different from the Simpack stick-slip model. That is why, just the ADAMS and the Simpack approaches are used here for the comparison with the second order dynamic friction model.

2. Static and Dynamic Friction Models

The idealized friction model of Coulomb simply distinguishes between sticking and sliding, Figure 1a. The friction force $F_R$ depends on the sliding velocity v and is realized by by combining an inequality with a simple relation

$$\left|F_R\right|\le {\mu }_s{F}_N\mathrm{if}v=0\mathrm{and}{F}_R=\frac{v}{\left|v\right|}{\mu }_d{F}_N\mathrm{if}\left|v\right|>0$$

(1)

The friction force is proportional to the normal force $F_N$ and characterized by the parameters ${\mu }_s$ and ${\mu }_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 because of the inequality representing stick. That is why, dry friction is usually approximated by an unambiguous function, Figure 1b. The regularized characteristics $F_R=F_R\left(v\right)$ introduces the fictitious velocity $v_s$, which defines the width of the regularization interval $-v_s\le v\le +v_s$ and the dynamic velocity $v_d$, that characterizes at $\left|v\right|\ge v_d$ the range of full sliding, where the simple relation $F_R=\pm {\mu }_d{F}_N$ applies. The transitions from the value pairs $(-v_d,-{\mu }_d{F}_N)\to (-v_s,-{\mu }_s{F}_N)$ and $(-v_s,-{\mu }_s{F}_N)\to (+v_s,+{\mu }_s{F}_N)$ as well as $(+v_s,+{\mu }_s{F}_N)\to (+v_s,+{\mu }_s{F}_N)$ 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 [2]. Even much simpler regularizations, assuming by $\mu ={\mu }_s={\mu }_d$ just one unique friction value, are in good conformity to dynamic measurements [4]. However, static friction models cannot maintain long-term stick. They describes the friction force just as a function of the sliding velocity $F_R=F_R\left(v\right)$. Figure 1b defines a commonly used regularized friction characteristics $F_R=F_R\left(v\right)$ by a set of “static” (${\mu }_s,v_s$) and “dynamic” (${\mu }_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 [6].

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 [5].

A software capable formulation of a friction force model takes also 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 [2], 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, which was introduced in [2] as a reference 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\ddot{z}=F_R-F_B=F_R(v_C-\dot{z})-\left({\sigma }_{0}z+{\sigma }_1\dot{z}\right)$$

(2)

where $m_b\ddot{z}$ approximates the inertia force of the fictitious bristle, $F_R=F_R(v_C-\dot{z})$ describes the friction force, $F_B={\sigma }_{0}z+{\sigma }_1\dot{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, ${\sigma }_0$ and ${\sigma }_1$ describe the stiffness and the damping of the fictitious bristle. The fictitious mass of the bristle is defined by

$$m_b=\frac{{\sigma }_1^2}{4{\sigma }_0}$$

(3)

which represents the aperiodic case of the homogenous second order differential equation (2) thus avoiding unwanted oscillations of the fictitious bristle. The second order dynamic bristle model implies with $s=\left[z,\dot{z}\right]$ two internal states. It is provided as a Matlab function in the Appendix by Listing 4.

The stick-slip models of Adams and RecurDyn describe the dynamic friction force just as a function of the contact point velocity $F_{Rd}=F_N{\mu }_d\left(v_C\right)$ and approximate the static friction force by a two-dimensional function $F_{Rs}=F_N{\mu }_s(v_C,x)$, where x is a fictitious displacement [3]. A smooth function of the contact point velocity $v_C$, which is not explained in detail in the user manual, describes the transition from ${\mu }_s$ to ${\mu }_d$. The fictitious displacement serves as an internal state and generates the static friction force $F_{Rs}$ or the static friction coefficient ${\mu }_s$ as a nonlinear function of x by adding a viscous damping term.

Simpack provides a stick-slip model which realizes, like the Coulomb’s approach in Figure 1a, a sudden drop from the static ${\mu }_s$ to the dynamic friction coefficient ${\mu }_d$ [3]. 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 ($\ddot{z}=0$ and $\dot{z}=v_C$) of the fictitious bristle dynamics (2) provides the required sticking force as

$$F_0=F_R\left(0\right)={\sigma }_{0}z+{\sigma }_1{v}_C$$

(4)

The second order dynamic friction model describes here the transition from the static to the dynamic friction force $({\mu }_s,F_s)\to ({\mu }_d,F_d)$ by a cubic polynomial and defines the friction force in the regularization range $-v_s\le v\le +v_s$ by a parabolic function. Then

$$\Delta v=\frac{F_0}{\left|F_0\right|}{v}_s{w}_F\mathrm{with}{w}_F=\left\{\begin{array}{cccc}1-\sqrt{1-\left|F_0\right|/F_s}\hfill & & \multicolumn{2}{c}{\mathrm{if}\left|F_0\right|\le F_s}\\ 1\hfill & & \multicolumn{2}{c}{\mathrm{elsewhere}}\end{array}\right.$$

(5)

delivers the corresponding horizontal shift. The static friction force is defined by $F_s=F_N{\mu }_s$ and $F_0/\left|F_0\right|$ adjusts the horizontal shift to the sign of the required sticking force value.

Figure 2 provides a simple test-bench which is used here to demonstrate the long-term sticking quality of the second order dynamic friction model and in the following for a comparison to the stick-slip models of ADAMS and Simpack.

A body of unit mass $m=1\mathrm{kg}$ is in contact to a horizontal rough plate. It is exposed to a horizontal force $F=F\left(t\right)$, which is continously increased in the interval $0\le t\le t_E$ from $F=0$ to the final value $F=F_E$. The friction parameters ${\mu }_s$, $v_s$ and ${\mu }_d$, $v_d$ model a regularized friction characteristics as defined in Figure 1b. The bristle parameters are adjusted to the body and the friction parameters by estimation a reference friction force of $F_{Rref}\approx 5\mathrm{N}$ and defining a reference bristle deflection of $z_{ref}=1·{10}^{-6}\mathrm{m}$. Then, ${\sigma }_0=F_{Rref}/z_{ref}=5·{10}^6\mathrm{N}/\mathrm{m}$ provides the bristle stiffness. The reference friction force corresponds with a reference mass of $m_{ref}=F_{Rref}/g=0.51\mathrm{kg}$ which provides the damping of the fictitious bristle as ${\sigma }_1=2\eta \sqrt{m_{ref}{\sigma }_0}=3192.75\mathrm{N}/(\mathrm{m}/\mathrm{s})$ if a viscous damping rate of $\eta =1$ is assumed hereby.

The Matlab script, provided in the Appendix by Listing 1, performs simulations with different step durations. It applies the standard implicit solver ode15s where the default tolerances are changed to RelTol=1.e-6 and AbsTol=1.e-9 because the reference bristle deflection was chosen very small in this example. The Matlab function dyn_fr_test_bench, provided in the Appendix by Listing 2, computes the dynamics of the simple friction test bench including the second order dynamic friction model as a set of first order differential equations. The Matlab functions Step3 and FrDyn2, defined in the Appendix by Listings 3 and 4 provide the step input and the dynamics of the second order friction model. The simulation results, plotted in Figure 3 demonstrate that the second order dynamic friction model perfectly maintains long-term stick.

The excitation force is slowly increased in the time interval $0\mathrm{s}\le t\le 30\mathrm{s}$ from $F=0$ to 99 % of the static friction force $F_s$, which is given here by $F_s={\mu }_s{F}_N=0.6·1\mathrm{kg}·9.81\mathrm{m}/{\mathrm{s}}^2=5.886\mathrm{N}$. In the subsequent time interval $30\mathrm{s}\le t\le 300\mathrm{s}$ the force $F_E=0.99·5.886\mathrm{N}=5.827\mathrm{N}$ is kept constant. The friction force $F_R\left(t\right)$, generated by the second order dynamic friction model counteracts the excitation force $F_R=F$, upper left plot in Figure 3. Note 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_{st}=u_{st}=F_E/{\sigma }_0=5.827\mathrm{N}/5.0·{10}^6\mathrm{N}/\mathrm{m}=1.1654·{10}^{-6}\mathrm{m}$. The friction force diagram $F_R\left(v\right)$ in the upper right plot of Figure 3, shows that the second-order dynamic friction 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\left(v\right)$. 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, Figure 4.

In the regularization range, the friction force characteristics $F_R=F_R\left(v\right)$ 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_{st}=0.99·F_s$ at the velocity of $v_{st}=0.9·v_s=0.9·{10}^{-3}\mathrm{m}/\mathrm{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 friction test-bench, defined in Figure 2, consists of a body in contact to a horizontal rough plate. The excitation force $F=F\left(t\right)$ is now slowly increased within $10\mathrm{s}$ from $F=0$ to the final value $F=F_E=1.05F_s=6.180\mathrm{N}$, which exceeds the adhesion limit of $F_s$ by 5 %. The excitation force $F=F\left(t\right)$ modeled by a third order polynomial reaches the adhesion limit $F=F_s=5.886\mathrm{N}$ at time $t=t_B=8.6808\mathrm{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 second order dynamic friction model (FrDyn2), the ADAMS stick-slip model, and the Simpack stick-slip model.

At first ($0\le t\le t_B$), the body remains in a quasi-static equilibrium, where the slowly increasing excitation force $F\left(t\right)$ is perfectly counteracted by a friction force $F_R\left(t\right)=F\left(t\right)$ generated by each of the friction models, as demonstrated in Figure 5a where the lines for $F\left(t\right)$ and $F_R\left(t\right)$ perfectly coincide in the time interval $0\le t\le 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 $\dot{z}$ as internal states $s=\left[z,\dot{z}\right]$. In the quasi-static equilibrium mode the velocity of the body and the time derivative of the bristle deflection are negligible small $v\approx 0$ and $\dot{z}\approx 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={\sigma }_{0}z+{\sigma }_1\dot{z}$. In the quasi-static equilibrium mode the body acceleration $a=\ddot{u}$ is negligible small too, as indicated in Figure 5b by the time history $a=a\left(t\right)$ in the time interval $0\le t\le 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\to {\sigma }_{0}z={\sigma }_{0}u$. The quasi-static force $F_R={\sigma }_{0}u$ equals the static friction force $F_s$ at the reference displacement of

$$u=x_R=\frac{F_s}{{\sigma }_0}=\frac{5.886\mathrm{N}}{5·{10}^6\mathrm{N}/\mathrm{m}}=1.1772·{10}^{-6}\mathrm{m}$$

(6)

The adhesion range is characterized by a vanishing sliding velocity ($v\approx 0$) and extends here to displacements in the range of $0\le u\le 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, 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$, Figure 5a and Figure 5d in particular. The transition from the static to the dynamic friction force $F_s\to F_d$ are modeled in the FrDyn2 and the ADAMS stick-slip model 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 characteristics (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 quite similar (solid black line). Most likely, ADAMS models the transition $F_s\to F_d$ by a 5th order polynomial. As a consequence, the FrDyn2 and the ADAMS stick-slip models produce slightly and slightly more delayed drops in the time histories of the computed friction forces, Figure 5d.

Figure 5b and Figure 5c illustrate the break-away effect at $t=t_B$ by the time histories of the body acceleration $\ddot{u}=a=a\left(t\right)$ and the body displacement $u=u\left(t\right)$. All friction models under consideration approximate sliding at $v\ge v_d$ by a constant friction force $F_R(v\ge v_d)=F_d=2.943\mathrm{N}$. Viscous components in the friction force are not considered here. The free body diagram in Figure 2b delivers the linear momentum

$$m\ddot{u}=F-F_R$$

(7)

for the body of mass $m=1\mathrm{kg}$. At $t>10\mathrm{s}$ which includes $t>t_B$ the applied force is defined by $F=F_E=1.05F_s=6.180\mathrm{N}$ and the friction force is represented by its dynamic value $F_R=F_d=2.943\mathrm{N}$. Then, the maximum acceleration of the body is defined by

$$a_{max}=\ddot{u}(t>10\mathrm{s})=(F_E-F_d)/m=(6.180\mathrm{N}-2.943\mathrm{N})/1\mathrm{kg}=3.237\mathrm{m}/{\mathrm{s}}^2$$

(8)

which is exactly reproduced by the friction models, Figure 5b.

5. Dynamic Response

A pulse load excitation, performed in [2], revealed the tendency of dynamic friction models to produce dynamic overshoots in the friction force time histories. That is why, the test-bench defined in Figure 2 is now exposed to excitation forces where the amplitude $F_E=0.95F_s=5.592\mathrm{N}$ is 5 % less than the static friction force $F_s=5.886\mathrm{N}$ and the step duration is varied from $t_E=0.1\mathrm{s}$ to $t_E=0.0001\mathrm{s}$. The corresponding simulation results are shown in Figure 6. The solid thin grey line represents the excitation force $F=F\left(t\right)$, the dashed blue line, the solid black line, and the dotted red line mark the friction forces $F_R=F_R\left(t\right)$ 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\left(t\right)=F\left(t\right)$ at a step duration of $t_E=0.1\mathrm{s}$, 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=\dot{u}$ but only on the displacement of the body u, Figure 7a and Figure 7b. In case of the FrDyn2 and the Simpack stick-slip models $F_R={\sigma }_{0}u$ holds, which provides the friction force $F_R=F_E=5.592\mathrm{N}$ at the steady state displacement $u=x_{st}^{D2}=x_{st}^S=5.592\mathrm{N}/5·{10}^6\mathrm{N}/\mathrm{m}=1.12·{10}^{-6}\mathrm{m}$, as indicated in Figure 7b and Figure 7c 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\mathrm{N}$ at the steady state displacement $u=x_{st}^A=0.71·{10}^{-6}\mathrm{m}$, Figure 7b and Figure 7c. As expected from the time histories $F_R=F_R\left(t\right)$ plotted in Figure 6a, the time histories of the body displacement reach their steady state values $u=x_{st}^{D2}=x_{st}^S$ and $u=x_{st}^A$ at $t>t_E$ without any overshoots, 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 simple friction test bench simplifies to

$$m\ddot{u}=F\left(t\right)-{\sigma }_{0}u\mathrm{or}m\ddot{u}+{\sigma }_{0}u=F\left(t\right)$$

(9)

where the quasi-static friction force is generated by the bristle compliance $F_R={\sigma }_{0}z$ and $z=u$ holds in addition. The simplified equation of motion (9) is characterized by the eigen-frequency ${\omega }_0=\sqrt{{\sigma }_0/m}$ and delivers the value and its of corresponding oscillation period as

$${\omega }_0=\sqrt{\left(5·{10}^6\mathrm{N}/\mathrm{m}\right)/\left(1\mathrm{kg}\right)}=2236{\mathrm{s}}^{-1}\mathrm{and}T=\frac{2\pi }{{\omega }_0}=0.003\mathrm{s}$$

(10)

As a consequence, even a rather short step duration of $t_E=0.01\mathrm{s}$ will still represent a subcritical excitation of the simple friction test bench. The time histories of the friction forces $F_R\left(t\right)$ exhibit just a tiny overshoot at $t>t_E=0.01\mathrm{s}$, Figure 6b.

The situation becomes complicated for step durations $t_E<T$, Figure 6c and Figure 6d. The Simpack stick-slip model (dotted red line) generates now significant overshoots, which amount to

$${\left.F_{Rmax}\right|}_{Simpack}^{t_E=0.001\mathrm{s}}=6.61\mathrm{N}\mathrm{and}{\left.F_{Rmax}\right|}_{Simpack}^{t_E=0.0001\mathrm{s}}=6.75\mathrm{N}$$

(11)

The values exceed the steady state value $F_{Rst}=F_E=5.592\mathrm{N}$ by $18.2$ % and $20.7$ % and even the static friction value $F_s=5.886\mathrm{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 $\left|F_R\right|\le F_s$, as expected from friction models, in general.

The friction models generate now friction forces which strongly depend on the body velocity $v=\dot{u}$ and the body displacement u, Figure 8a and Figure 8b.

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_{st}=1.12·{10}^{-6}\mathrm{m}$ with a small overshoot shortly after the step duration of $t_E=0.001\mathrm{s}$, Figure 8c. The FrDyn2 model overshoots and partly slides resulting in displacements at $t>t_E$ which exceed with $u_{st}^{D2}=2.27·{10}^{-6}\mathrm{m}$ the quasi steady state value of $x_{st}=1.12·{10}^{-6}\mathrm{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\left(t\right)$ shows a strange behavior for the ADAMS stick-slip model, solid black line in Figure 8c. At first ($t<3.5·{10}^{-3}\mathrm{s}$), it approaches the quasi steady state value $u=x_{st}^{D2}=x_{st}^S=1.12·{10}^{-6}\mathrm{m}$ of the FrDyn2 and Simpack solution and then ($4·{10}^{-3}\mathrm{s}<t\le 10·{10}^{-3}\mathrm{s}$) it starts to decrease very slowly but the simulated time interval $0\le t\le 10·{10}^{-3}\mathrm{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\left(t\right)$ as plotted in Figure 9,

where the section shown in Figure 8c has a light grey background. It seems that ADAMS applies in its stick-slip model different time constants for the increase and decrease of the body displacement $u=u\left(t\right)$. The force excitation with a step duration of $t_E=0.001\mathrm{s}$ is much faster than the dynamics of the friction test bench. The body displacement reaches its maximum value at $t\approx 0.0035\mathrm{s}$ which corresponds to oscillation period computed in (10). The decay from the maximum displacement to the steady state value takes about $0.05\mathrm{s}$, which is fourteen times as much. This strange behavior was also reported in [3] where 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.8F_s$, the body returned to its initial position. However, in the present example a small but permanent deviation of $\delta u=u_{st}^A-x_{st}^A=0.12·{10}^{-6}\mathrm{m}$ remains, Figure 9. This indicates that the ADAMS stick-slip model also tends to slip partly, 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 [2] to asses different friction models in a more practical example. The model consists of three cable $m_C$ and two trolley $m_T$ masses, Figure 10.

At the beginning $(t=0)$ the towing trolley is fixed at $u_{TT}\left(0\right)=1.5529\mathrm{m}$. The equilibrium position of the cable system places the trolleys at $x_{T1}\left(0\right)=u_{TT}\left(0\right)/3$, $x_{T2}\left(0\right)=2u_{TT}\left(0\right)/3$, and locates the cable masses at $x_{C1}\left(0\right)=u_{TT}\left(0\right)/6$, $x_{C2}\left(0\right)=u_{TT}\left(0\right)/2$, $x_{C3}\left(0\right)=5u_{TT}\left(0\right)/6$, as well as $z_{C1}\left(0\right)=z_{C2}\left(0\right)=z_{C3}\left(0\right)=0.96698\mathrm{m}$. The non-holonomic constraint ${\dot{u}}_{TT}=v_{TT}$ relates the movements of the towing trolley to a pre-defined velocity profile $v_{TT}=v_{TT}\left(t\right)$. The velocity profile, defined by the solid grey lines in Figure 11a, models an extension maneuver, which moves the towing trolley from the initial position $u(t=0)=u_{TT}\left(0\right)=1.5529\mathrm{m}$ to a final position of $u_{TT}(t>5\mathrm{s})=5.553\mathrm{m}$.

The Matlab simulation with the second order dynamic friction model (FrDyn2) generates the output at every simulation step. It applies the ode15s solver with error tolerances of RelTol=1.e-6 and AbsTol=1.e-9. The ADAMS and Simpack simulations were performed with an output step size of $\Delta t=1·{10}^{-4}\mathrm{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\mathrm{s}$. After that the trolleys perform to and fro motions which at $t\ge 5\mathrm{s}$ are indicated by sign changes in the time histories of trolley velocities $v_{T1}\left(t\right)$ and $v_{T2}\left(t\right)$. An arrow pointing from $v_{T1}(t\approx 6.47\mathrm{s})=0$ over $F_{N1}$ down to $F_{R1}$ highlights such an event, in particular. The dynamic motions of the cable masses induce variations in the normal forces $F_{N1}$ and $F_{N2}$ acting between the trolleys and the rail, Figure 11b. The time histories of the velocities $v_{T1}\left(t\right)$, $v_{T2}\left(t\right)$ and the normal forces $F_{N1}\left(t\right)$, $F_{N2}\left(t\right)$ generated with FrDyn2 and the stick-slip models of ADAMS and Simpack match nearly perfectly. However, the time histories of the friction forces $F_{R1}\left(t\right)$, $F_{R2}\left(t\right)$ exhibit some discrepancies, 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 of the trolley velocity $v_{T1}$ at $t\approx 6.47\mathrm{s}$ and a sticking period of trolley 2 in the time interval $7.2\mathrm{s}\le t\le 7.8\mathrm{s}$.

In the very short time interval $6.45\mathrm{s}\le t\le 6.48\mathrm{s}$ the normal force between trolley 1 and the rail amounts to $F_{N1}=21.2\mathrm{N}$ as indicated in Figure 11b. The friction values defined in Figure 10 provide in this case a static friction force of $F_s={\mu }_s{F}_{N1}=0.08·21.2=1.696\mathrm{N}$ and a dynamic friction force of $F_d={\mu }_d{F}_{N1}=0.05·21.2=1.06\mathrm{N}$. At times $t=6.45\mathrm{s}$ and $t=6.48\mathrm{s}$ the first trolley is in a full sliding mode, indicated in Figure 12a by the friction forces $F_{R1}(t=6.45\mathrm{s})=+F_d=1.06\mathrm{N}$ and $F_{R1}(t=6.48\mathrm{s})=-F_d=-1.06\mathrm{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_{T1}\to 0$ with the dynamic force value and does not reproduce a potential velocity dependent increase of 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\mathrm{s}\le t\le 6.47\mathrm{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 $\left|F_{R1}\right|\le F_s$, which results in $F_{R1}^{D2}=F_{R1}^A=-1.696\mathrm{N}$ at $t\approx 6.47\mathrm{s}$. However, The Simpack stick-slip model overshoots and produces the peak-value of $F_{R1}^S=-2.685\mathrm{N}$ which exceeds by nearly 60 % the static friction force $F_s$ or $-F_S$, respectively.

The sticking period $7.2\mathrm{s}\le t\le 7.8\mathrm{s}$ is represented quite similar by the friction models under consideration, 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\approx 7.2\mathrm{s}$. However, the small time delay of the peak values visible in Figure 12a is not noticeable because of the large time interval applied in this plot. The Simpack stick-slip model is based on the 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 model (FrDyn2) with the commercial stick-slip models of ADAMS and Simpack. The FrDyn2 model was introduced in [2] as a reference and tested versus the LuGre and a simple regularized friction model. The stick-slip models of ADAMS and Simpack as well as RecurDyn are described and investigated in [3].

The comparison is performed here with a simple 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 characteristics $F_R=F_R\left(v\right)$ 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 to a suitable alternative to commercial stick-slip models.