# Modeling and Nonlinear Response of the Cam-Follower Oblique-Impact System

Authors: Yong-Feng Yang; Yu Lu; Ting-Dong Jiang; Na Lu
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

In order to quickly and accurately analyze the complex behavior of cam-follower oblique-impact system, a mathematical model which can describe separation, impact, and contact was established in this paper. The transient impact hypothesis was extended, and the oblique collision model was established by considering the tangential slip. Moreau time-stepping method was employed to solve the linear complementarity problem which transformed by the oblique-impact equations. The simulation results show that the cam and follower kept permanent contact when the cam rotational speed was low. With the increase of the cam rotational speed, the cam and follower would be separated and then impact under the gravity action. The system performance shows very complex nonlinear characteristics.


---

## Body

## 1. Introduction

Cam-follower devices are important classical impacting systems which are used in a wide range of applications. The rotation of the cam at some constant speed provides the force to operate the follower. The most common example is valve trains of internal combustion engines, where the cam rotation imparts through the follower the proper motion to the engine valves while a spring provides the restoring force necessary to maintain contact between the components [1]. The impact of cam-follower is a typical oblique impact. In analyses using nonsmooth dynamics, oblique impact of rough bodies in an unsymmetrical configuration can result in self-locking or “jam” at the sliding contact if the coefficient of friction is sufficiently large: this has been termed Painleve’s paradox [2]. Because of this, the follower may exhibit very complex behavior including bifurcations and chaos.Alzate et al. [3] show that a sudden transition to chaos is observed in an experimental cam-follower. In their mathematical model, the tangential friction effect is ignored. Considering the effect of tangential friction, Sundar et al. [4] study the linear response of a cam-follower. But they do not discuss the nonlinear characteristic. Some researches studied the beam’s analysis and experimental methods that are similar and can be used in our study. For example, Ding et al. [5, 6] investigates the convergence of the Galerkin method for the dynamic response of an elastic beam resting on a nonlinear foundation with viscous damping subject to a moving concentrated load. Zhang et al. [7] present experimental verifications of vibration suppression for a cantilever beam bond with a piezoelectric actuator by an adaptive controller. In our research, the transient impact hypothesis is extended, and the oblique collision model is established by considering the tangential slip. Moreau time-stepping method is employed to solve the linear complementarity problem which transforms by the oblique-impact equations. The complex behavior including bifurcations and chaos is shown by the simulation.


## 2. Linear Complementarity Problem

A linear complementarity problem (LCP) is a problem of the formy
=
M
x
+
b, y
≥
0, x
≥
0, and y
T
x
=
0, where M is an n
×
n matrix andx,y, andb are an n-dimensional vector [8].LCP equations may have unique solution, no solution, or multiple solutions. Many problems in scientific computing and engineering applications may lead to solutions of LCP. For example, the LCP may arise from application problems such as the linear and quadratic programming, the economies with institutional restrictions upon prices, the optimal stopping in Markov chain, and the free boundary problems. Several computational methods have been developed for solving LCPs. The two basic functions of the LCP are Upr function and Sgn function:(1)
U
p
r
=
0
x
>
0
-
∞
,
0
x
=
0
⌀
x
<
0
,
S
g
n
=
∂
x
=
1
x
>
0
-
1,1
x
=
0
-
1
x
<
0
.In mechanics, Upr function is used to simulate geometry and dynamics of unilateral constraint, and Sgn function is used to simulate various kinds of dry friction.


## 3. Oblique-Impact Model

A cam-follower with oblique-impact on the contact point is studied in Figure1. The summary of the general notations used in Figure 1 is shown in Notations. By extending the transient impact hypothesis and considering the tangential slip, the oblique collision model is established.Figure 1
Sketch of cam-follower oblique-impact model.In the present work, the normal contact between rigid bodies is characterized by a set-valued force law called Signorini’s condition [9]. Figure 2 shows two convex rigid bodies apart from each other by a relative normal gap or distance denoted by g
N. The relative normal gap is nonnegative due to bodies’ impenetrability condition, being the two bodies in contact with each other when g
N
=
0. According to Signorini’s condition, the normal contact conditions are(2)
g
N
=
0
∧
λ
N
≥
0
c
l
o
s
e
d
c
o
n
t
a
c
t
,
g
N
>
0
∧
λ
N
=
0
o
p
e
n
c
o
n
t
a
c
t
.Equation (2) represents a complementarity behavior; always zero is the product of the relative normal gap and normal contact force; that is, g
N
λ
N
=
0.Figure 2
Sketch of rigid body contact: (a) normal distanceg
N and (b) tangential and normal force.
(a)
(b)The classical Coulomb friction law is another typical example that can be considered as a set-valued force law [9]. For the two conditions sliding and sticking, the tangential forces are(3)
-
λ
T
=
μ
λ
N
S
g
n
γ
T
,
γ
T
≠
0
s
l
i
d
i
n
g
,
-
μ
λ
N
≤
λ
T
≤
μ
λ
N
s
t
i
c
k
i
n
g
,where μ is the friction coefficient; λ
N is the normal contact force; and γ
T is the relative tangential velocity. These two conditions can be summarized by a set-valued force law as [9](4)
-
λ
T
∈
μ
λ
N
S
g
n
γ
T
.Tangential Coulomb’s contact states can be expressed by Sgn function. It also represents a complementarity behavior.


## 4. Dynamics Equations of the Follower

For the model shown in Figure1, the dynamics of follower under different speed conditions are considered. There are three states of the follower and cam: separate, contact, and impact. The dynamic equations of the follower can be expressed as [10](5)
M
u
˙
+
h
=
0
s
e
p
a
r
a
t
e
,
(6)
M
u
˙
+
h
-
w
N
λ
N
-
w
T
λ
T
=
0
c
o
n
t
a
c
t
,
(7)
M
u
+
-
u
-
-
w
N
Λ
N
-
w
T
Λ
T
=
0
i
m
p
a
c
t
,where u is generalized velocity and here is the follower’s angular velocity θ
˙
f; M
=
J
o
=
m
l
2
/
3
+
m
r
2
/
4 is the follower’s moment of inertia; h
=
m
g
l
cos
⁡
(
θ
f
)
/
2 is the follower’s external force; w is the generalized coefficient in contact point; Λ is the force of the follower, and u
+ and u
- are the velocity before and after the impact point, respectively. Notice that λ
N is a continuous boundary parameter when there is no impact. λ
T is a noncontinuous parameter when there is viscous-to-slip transition or change in the direction of relative speed. N denotes normal direction and T denotes tangential direction.The relationship betweenw
N, w
T and contact point’s relative speeds is(8)
γ
N
=
w
N
u
+
∂
w
N
∂
t
,
γ
T
=
w
T
u
+
∂
w
T
∂
t
.In order to get the unknown parametersw
N and w
T, we need to calculate the distance function in Figure 1:(9)
g
N
=
x
0
+
e
cos
⁡
θ
c
s
i
n
θ
f
-
e
sin
⁡
θ
c
cos
⁡
θ
f
-
r
-
R
,
g
T
=
x
0
+
e
cos
⁡
θ
c
c
o
s
θ
f
+
e
sin
⁡
θ
c
sin
⁡
θ
f
.According to the results of [11],(10)
g
˙
=
w
q
˙
+
∂
w
∂
t
.Therefore,(11)
w
N
=
x
0
+
e
cos
⁡
θ
c
c
o
s
θ
f
+
e
sin
⁡
θ
c
s
i
n
θ
f
,
w
T
=
-
x
0
+
e
cos
⁡
θ
c
s
i
n
θ
f
+
e
sin
⁡
θ
c
c
o
s
θ
f
.There are three equations for the cam-follower impact system. In order to quickly and accurately analyze the complex behavior of cam-follower oblique-impact system, it is better to establish one mathematical equation which can describe separation, impact, and contact motion. Equation measure is introduced here. Equation (6) is multiplied by the Lebesgue measure d
t and (7) is multiplied by the atoms measure d
η. Then add the two equations:(12)
M
u
˙
d
t
+
u
+
-
u
-
d
η
+
h
d
t
-
w
N
λ
N
d
t
+
Λ
N
d
η
-
w
T
λ
T
d
t
+
Λ
T
d
η
=
0
.This equation can be abbreviated as(13)
M
d
u
+
h
d
t
-
w
N
d
P
N
-
w
T
d
P
T
=
0
,where d
u
=
u
˙
d
t
+
u
+
-
u
-
d
η is velocity measure. It contains two parts: one is the continuous part, Lebesgue measure u
˙
d
t; the other is the noncontinuous part, Dirac measure. When the cam-follower is in the nonimpact state, the velocity measure d
u
=
u
˙
d
t. The impact measure which relates to the Lagrange multiplier also contains two parts: d
P
=
λ
d
t
+
Λ
d
η. The separation, impact, and contact equations are expressed as the unified form in (13) [10].Moreau time-stepping method [10] is employed to discretize (13). We can get(14)
M
M
u
E
-
u
A
+
h
M
Δ
t
-
w
N
M
p
N
-
w
T
M
p
T
=
0
,where A and E represent starting and ending time. M is the middle time of A and E. For the normal impact parameter p
N, apply the Upr function:(15)
-
p
N
∈
U
p
r
ξ
N
⟺
p
N
≥
0
,
ξ
N
≥
0
,
p
N
ξ
N
=
0
,where ξ
N
=
γ
N
+
+
ε
N
γ
N
-.In order to get the expectation LCP condition, we need to separate Sgn function to two Upr functions. For the tangential impact parameterp
T, we can get(16)
-
p
T
∈
μ
p
N
S
g
n
ξ
T
⟺
∃
ξ
R
,
ξ
L
s
.
t
.
μ
p
N
+
p
T
≥
0
,
ξ
R
≥
0
,
μ
p
N
+
p
T
ξ
R
=
0
μ
p
N
-
p
T
≥
0
,
ξ
L
≥
0
,
μ
p
N
-
p
T
ξ
L
=
0
ξ
T
=
ξ
R
-
ξ
L
.For the sake of brevity, we definep
R
≔
μ
p
N
+
p
T and p
L
≔
μ
p
N
-
p
T, where “≔” means “is defined as.” The second and third expressions of (16) are P
R
≥
0
,
ξ
R
≥
0
,
P
R
ξ
R
=
0 and P
L
≥
0
,
ξ
L
≥
0
,
P
L
ξ
L
=
0. The LCP condition is(17)
0
≤
ξ
N
ξ
R
p
L
⊥
p
N
p
R
ξ
L
≥
0
.The full LCP equation of follower can be expressed as(18)
ξ
N
ξ
R
p
L
=
w
N
M
M
M
-
1
w
N
M
-
w
T
M
μ
w
N
M
M
M
-
1
w
T
M
0
w
T
M
M
M
-
1
w
N
M
-
w
T
M
μ
w
N
M
M
M
-
1
w
T
M
1
2
μ
-
1
0
p
N
p
R
ξ
L
+
w
N
M
M
M
-
1
h
M
Δ
t
+
1
+
ε
N
γ
N
A
w
T
M
M
M
-
1
h
M
Δ
t
+
1
+
ε
N
γ
T
A
0
.From the solution of (17) and (18), we can get p
N and p
R. Then the generalized velocity can be calculated by(19)
u
E
=
M
M
-
1
w
N
M
-
w
T
M
μ
p
N
+
M
M
-
1
w
T
M
p
R
+
M
M
-
1
h
M
Δ
t
+
u
Aand the generalized displacement can be gotten by the definition of Moreau time time-stepping method.


## 5. Simulation of the Dynamics Response

MATLAB is employed for simulation of the equations. In order to get a proper accuracy, the time step of Moreau midpoint method must be very small. We choose 0.0002 s as the time step, and the values of other parameters can be found in Table1.Table 1
Cam-follower oblique-impact numerical simulation parameters.
Parameter
Value
Length of the follower,L
0.8 m
Radius of the follower,r
0.01 m
Mass of the follower,m
1.975 kg
Eccentricity of the cam,e
0.036 m
Radius of the cam,R
0.1 m
Coefficient of restitution in normal direction
0.42
Coefficient of restitution in tangential direction
0.4
Friction coefficient,μ
0.2
Distance between bear fix point of the follower and rotation center of the cam,x
0
0.4 m
Initial angular position of the cam,θ
0
0 radWith the increasing of the cam rotational speed, the state of the cam-follower is transformed from the initial permanent contact to the separation and impact. The system performs very complex nonlinear characteristics, such as period, quasiperiod, and chaos response. The bifurcation diagrams are shown in Figure3.Figure 3
Bifurcation diagram with the change of cam speed: (a) angular velocity and (b) angular displacement.
(a)
(b)


### 5.1. Permanent Contact

When the cam rotational speed is low, the constraining force between cam and follower is less than the follower’s restoring force. Therefore, the cam and follower keep in contact. The angular displacement of time history diagram and phase diagram is shown in Figure4 when ω
=
12 rad/s. The motion of cam and follower remained the same (there are a solid line and a dashed line in Figure 4(a)), and the phase plane shows a smooth circle.Figure 4
Cam and follower permanent contact whenω
=
12 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)


### 5.2. Period Motion

With the increasing of the cam speed, the cam and follower would be separated. When the cam and follower impact, there are minor collisions after the first impact, and then follower would contact the cam’s surface until the next separation. As shown in Figure5, the motion of the follower is periodic 1 or the orbit is P
(
∞
,
1
).Figure 5
P
(
∞
,
1
) motion when ω
=
16 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)When the cam speed goes higher, as shown in Figures6 and 7, the responses are P
1,1
,
P
2,2
,
P
4,4
,
…. From the time history of these P
(
n
,
n
) responses, the common point is that there is only 1 impact in each cam period.Figure 6
P
(
1
,
1
) motion when ω
=
17.4 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)Figure 7
P
(
4
,
4
) motion when ω
=
18.5 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)There also existP
(
∞
,
n
) motions. The rotation speeds of these motions are between two chaos areas in Figure 3. Figure 8 shows the P
(
∞
,
2
) motion. For these motions, the period of impact is greater than the cam. When the collision occurs, the normal velocity of cam and follower is close. Therefore, collision could not change the follower’s angular velocity direction.Figure 8
P
(
∞
,
2
) motion when ω
=
22.5 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)


### 5.3. Chaos Motion

Generally, for vibroimpact system, there are two ways for Poincaré section selection:(
1
) to select the fixed phase section of trajectory as Poincaré and (
2
) to select the collision surface as Poincaré section. In this paper, the collision surface is used:(20)
Σ
=
θ
,
θ
˙
,
t
∈
R
2
×
R
+
∣
ξ
N
=
0
,
γ
N
A
<
0
.
ξ
N
=
0 represents the collision point, and the mirror collision can be excluded by γ
N
A
<
0.Figure9 shows the response when ω
=
18.84 rad/s. The collision does not follow the periodic variation and became chaotic. Figure 10 shows the response when ω
=
21.7 rad/s. Compared with Figure 9, the chaos response shows more characteristics of slip phenomenon with 21.7 rad/s.Figure 9
Chaos motion whenω
=
18.84 rad/s: (a) time history; (b) phase diagram; and (c) Poincaré map.
(a)
(b)
(c)Figure 10
Chaos motion whenω
=
21.7 rad/s: (a) time history; (b) phase diagram; and (c) Poincaré map.
(a)
(b)
(c)


## 5.1. Permanent Contact

When the cam rotational speed is low, the constraining force between cam and follower is less than the follower’s restoring force. Therefore, the cam and follower keep in contact. The angular displacement of time history diagram and phase diagram is shown in Figure4 when ω
=
12 rad/s. The motion of cam and follower remained the same (there are a solid line and a dashed line in Figure 4(a)), and the phase plane shows a smooth circle.Figure 4
Cam and follower permanent contact whenω
=
12 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)


## 5.2. Period Motion

With the increasing of the cam speed, the cam and follower would be separated. When the cam and follower impact, there are minor collisions after the first impact, and then follower would contact the cam’s surface until the next separation. As shown in Figure5, the motion of the follower is periodic 1 or the orbit is P
(
∞
,
1
).Figure 5
P
(
∞
,
1
) motion when ω
=
16 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)When the cam speed goes higher, as shown in Figures6 and 7, the responses are P
1,1
,
P
2,2
,
P
4,4
,
…. From the time history of these P
(
n
,
n
) responses, the common point is that there is only 1 impact in each cam period.Figure 6
P
(
1
,
1
) motion when ω
=
17.4 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)Figure 7
P
(
4
,
4
) motion when ω
=
18.5 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)There also existP
(
∞
,
n
) motions. The rotation speeds of these motions are between two chaos areas in Figure 3. Figure 8 shows the P
(
∞
,
2
) motion. For these motions, the period of impact is greater than the cam. When the collision occurs, the normal velocity of cam and follower is close. Therefore, collision could not change the follower’s angular velocity direction.Figure 8
P
(
∞
,
2
) motion when ω
=
22.5 rad/s: (a) time history of θ
f (solid line) and θ
f
c (dashed line) and (b) phase diagram.
(a)
(b)


## 5.3. Chaos Motion

Generally, for vibroimpact system, there are two ways for Poincaré section selection:(
1
) to select the fixed phase section of trajectory as Poincaré and (
2
) to select the collision surface as Poincaré section. In this paper, the collision surface is used:(20)
Σ
=
θ
,
θ
˙
,
t
∈
R
2
×
R
+
∣
ξ
N
=
0
,
γ
N
A
<
0
.
ξ
N
=
0 represents the collision point, and the mirror collision can be excluded by γ
N
A
<
0.Figure9 shows the response when ω
=
18.84 rad/s. The collision does not follow the periodic variation and became chaotic. Figure 10 shows the response when ω
=
21.7 rad/s. Compared with Figure 9, the chaos response shows more characteristics of slip phenomenon with 21.7 rad/s.Figure 9
Chaos motion whenω
=
18.84 rad/s: (a) time history; (b) phase diagram; and (c) Poincaré map.
(a)
(b)
(c)Figure 10
Chaos motion whenω
=
21.7 rad/s: (a) time history; (b) phase diagram; and (c) Poincaré map.
(a)
(b)
(c)


## 6. Conclusions

The cam and follower oblique collision model is established in this paper. In order to simulate the response, the equations of the separation, impact, and contact motion are transformed to a linear complementarity problem. The main conclusions are as follows:(1)
The contact and tangential impact law represent a complementarity behavior. LCP equation method could be used to solve the oblique collision problem.
(2)
The system performs very complex nonlinear characteristics, such as period, quasiperiod, and chaos response. For every motion, there is a special collision phenomenon.



---

*Source: 6142501-2016-07-19.xml*