# Effect of Surface Morphology on Dynamic Characteristics of Cam-Follower Oblique Impact System

Authors: Wangqun Deng; Haibo Feng; Jianjun Wang; Yongfeng Yang
Journal: Shock and Vibration
 (2019)

---

## Abstract

The effects of the fractal characteristics on the dynamic characteristics of the oblique impact system are studied via considering the normal contact stiffness with fractal dimension parameter of the microstructure of the interface in this paper. The normal distribution is used to describe the trend of contact line length. The effects of the position parameters and the scale parameters are analyzed. The simulation results show that the increasing of the interface roughness can make the system response more stable. The system response eventually evolves into chaos as the fractal dimension increases. The system response changes between quasiperiod and chaos affected by the contact line length with a normal random distribution.


---

## Body

## 1. Introduction

Mechanical systems are usually modeled as multibody systems with nonsmoothness. Typical examples are the vibration and noise produced in engines, impact print hammers, railway brakes, or chattering of machine tools. These effects are due to the nonsmooth characteristics such as impacts, clearances, intermittent contacts, dry friction, unbalances, crack, or a combination of these effects [1–5]. In nonsmooth systems, the time evolution of displacements and velocities is not requested to be smooth. Especially for the impact, the velocities are even allowed to undergo jump at certain time instances in order to fulfill the kinematical restrictions.Cam-follower devices are very important in impacting systems and widely applied in actual engineering systems. The rotation of the cam at special constant speed provides the force to operate the follower. The most common example is the valve train of the internal combustion engine. The cam rotation imparts force through the follower to the engine valve, while a spring provides the restoring force necessary to maintain contact between the two components [6, 7].The impact on the cam-follower is a typical oblique impact. By using nonsmooth dynamics, oblique impact on rough bodies with an unsymmetrical configuration can result in self-locking or “jam” at the sliding contact if the coefficient of friction is large enough [8]. It has been observed that, under the variation of the cam-follower system parameters, the responses may exhibit complex behaviors including bifurcations and chaos [9, 10]. For example, Zhang [11] designed a new experimental device to verify the vibration suppression of a cantilever beam bond with a piezoelectric actuator by an adaptive controller. Ding et al. [12, 13] investigated the convergence of the Galerkin method of an elastic beam resting on a nonlinear foundation with viscous damping subject to a moving concentrated load. Alzate [14, 15] shows that a sudden transition to chaos is observed in an experimental cam-follower by ignoring the tangential friction effect. Considering the effect of tangential friction, Sundar et al. [16] studied the linear response of a cam-follower. In addition, the normal contact stiffness is an important parameter affecting the dynamic characteristics of the contact between the cam and follower when considering the surface topography [17]. In traditional research, the stiffness of the impact system is mostly studied as a deterministic linear or nonlinear stiffness. However, the mechanical roughness of the surface topography may affect the wear, friction, and contact deformation, which make the stiffness no longer determined [18]. Majumdar and Bhushan [19] found that the mechanical processing surface has self-similar fractal characteristics and proposed the contact fractal theory and the contact fractal model (MB model), which uses the fractal dimension D and the fractal dimension parameter G to characterize the mechanical roughness of the surface topography.The contact of the cam-follower system consists of a round bar and a cam, which subject to the two intersecting cylindrical contact model. When the two parts contact, the contact surface is equivalent to an ellipse and the length of the contact line is equal to the length of the ellipse’s long axis. During the operation of the cam-follower system, the length of contact line is constantly changing due to the continuous change of the contact force between the follower and cam. Therefore, the length of the contact line should be dynamically changed in the actual analysis.The motivation of this paper is to establish the normal contact stiffness model based on the contact fractal theory via considering the fractal dimensionD and the fractal dimension parameter G. Normal distribution will be used to describe the variation of the contact line length. The effects of the normal contact stiffness on the dynamic characteristics of the oblique impact system are studied.


## 2. Mechanical Model and Dynamic Equation

### 2.1. System Description

Figure1 shows the schematic diagram of the cam-follower system. The follower makes a fixed axis rotation around O. The cam rotates around a fixed point O′, R is the radius, and ω is the rotating angular velocity of the cam. The eccentric e is the distance between the geometric center C and the rotation point O′. The angular displacement of the cam is θc=∠CO′x′, and the angle between the follower and the horizontal axis is θf. At the initial time, A is the cam and follower contact point. When the system is rotating, Bf and Bc are the impact points on the follower and cam, respectively. iAj is a local coordinate system, and note δi as the distance BfBc→, which indicates the deformation of contact springs between the follower and cam. Obviously, the size of the δi can be used to determine whether the impact occurs or not, that is, 0<δi indicates the separate state, and 0≥δi indicates the contact state.Figure 1
Schematic diagram of the cam-follower system.


### 2.2. The Equation of Motion

During the operation of the cam, the follower will show two kinds of motion states, namely, the contact state and the separate state, and the topology of the system will change with time.Using the Lagrangian approach, the follower motion equation can be obtained:(1)Jθ¨f+mgl2cosθf=0,separate  state,Jθ¨f+mgl2cosθf=Fnl1+μfsgnvrtFn,contact  state,where J is the follower’s moment of inertia, μf is the friction factor between follower and cam, and Fn is the normal contact force, and the expression is(2)Fnt=−Knδi−cδ˙i,where Kn and c represent the equivalent stiffness and damping coefficient of the normal contact, respectively.Using the pure damping model which was proposed by Padmanabhan [20], the given damping coefficient is(3)c=βKnδi,where β is the impact damping factor.


### 2.3. Impact Law

In this model, the impact time is very short, the displacement of the various parts is constant, and the speed changes before and after the impact. By ignoring the plastic characteristic, the oblique impact problem studied here only involves the elastic impact. Considering the impact recovery coefficient, the impact process is divided into two stages: compression and recovery. The velocity relationship before and after the impact is(4)θ˙fti+=−ξθ˙fti−+1+ξ−x˙Bcsinθf+y˙Bccosθfl1,where ξ is the coefficient of impact recovery and x˙Bc and y˙Bc are given by the following expressions:(5)x˙Bc=−eθ˙csinθc−Rθ˙fccosθfc,y˙Bc=eθ˙ccosθc−Rθ˙fcsinθfc.


### 2.4. Normal Contact Stiffness Fractal Model

The contact stiffness is related to the material properties of the two contacting objects and the geometry of the contact surface. For a specific mechanical system, the material properties of the components have been determined, and only the contact surface geometry may affect the contact stiffness coefficient.The actual mechanical interface consists of two rough surfaces. According to the results of Reference [21], the contact area can be reduced to a rigid plane in contact with a rough surface. Figure 2 shows a schematic diagram of the contact between a microconvex body and a rigid plane. When the contact occurs, the real area of the rough surface accounts for only a small part of the nominal contact area. Therefore, the contact problem of the mechanical joint surface can be handled with the contact of the convex body. When the surface is stressed, the microconvex body will occur elastic or plastic deformation. We focus on the elastic deformation of the local components.Figure 2
Schematic diagram of the contact between a microconvex body and a rigid plane. (a) Unloaded state. (b) Loaded state.
(a)
(b)The single microconvex body on the rough surface is equivalent to the sphere, andR is recorded as its equivalent radius of curvature.When the normal loadp is applied to the contact point, δ is the normal contact deformation. The radius of the contact area is R. According to Reference [20], the relationship between p and δ is(6)p=43ER1/2δ3/2,where E is the modulus of elasticity of the two contact materials:(7)1E=1−ϑf2Yf+1−ϑc2Yc,where ϑ and Y represent Poisson’s ratio and Young’s modulus of the material, respectively. Subscripts f and c represent the follower and cam.According to equation (6), the normal contact stiffness of a single microconvex body contact with a rigid plane is(8)kn=2ER1/2δ1/2.According to the geometric relationship between the microconvex deformation and the typical value of the fractal roughness parameterG in Reference [22], R≫δ, and the expression of the cross-sectional area a′ of the microconvex body is(9)a′=2πRδ.Substituting equation (9) into (8), we can obtain(10)kn=2Ea′2π.According to the contact fractal theory [23], the distribution function of the contact point of the microcontact area a′ is(11)na′=D2al′D/2a′D/2+1,(12)al′=2−DDAr,where al′ is the actual contact area of the maximum contact point, Ar is the real contact area of the rough surface, and D is the fractal dimension of rough surface which reflects the irregularity and complexity of surface profile on all scales. Considering the two-dimensional fractal problem, the range D∈1,2. According to the Hertz contact theory, when contact occurs, the contact area of two objects is elliptical, and the real contact area Ar is about 10−5m2.Assume that the rough surface is isotropic in the microtopography. According to the contact fractal theory, the microconvex body contact deformation belongs to the plastic deformation whena′<ac′, and the contact deformation belongs to the elastic deformation when a′>ac′. ac′ is determined by Reference [24] as(13)ac′=G2kφ/22/D−1,where G is the fractal roughness parameter of the rough surface. Larger G means rougher contact surface. k=H/σy is a factor related to the softer material hardness H and the yield strength σy, and k is generally taken as 2.8. φ=σy/E is a factor related to σy and the equivalent elastic modulus E.Thus the normal contact stiffness of the joint surface is(14)Kn=∫ac′al′knna′da′.Substituting Equations (10) and (11) into (14), we can obtain(15)Kn=2ED2π1−Dal′D/2al′1−D/2−ac′1−D/2.According to the relationshipa′=2a, we have(16)al′=2al,ac′=2ac.Substituting Equation (16) into (15), the normal contact stiffness can be obtained(17)Kn=2EDπ1−DalD/2al1−D/2−ac1−D/2.For a definite rough surface,ac is constant, and the normal stiffness depends on the real contact area al of the maximum contact point.


### 2.5. Contact Line Length with Random Normal Distribution

Using the Hertz contact theory in Reference [24], the equivalent stiffness coefficient is given by(18)Kn=π4Elλ,where lλ represents the length of the contact line.Assume that the length of the contact line obeys the law of normal distribution, that is, the contact line length is controlled by the position parameterμ and the scale parameter σ. μ and σ, respectively, determine the position and amplitude of the normal distribution. Figure 3 shows the numerical simulation of the normal contact line length with random distribution. It is found that the amplitude range of the contact line length increases as the scale parameter σ increases. Obviously, the position parameter μ and the scale parameter σ affect the dynamic response of the system. σ=0 indicates that the normal contact line length is constant, and σ≠0 indicates that the length is a random parameter.Figure 3
Normally distributed random contact line length. (a)μ=0.01m,σ=0.001m. (b) μ=0.01m,σ=0.002m. (c) μ=0.01m,σ=0.003m.
(a)
(b)
(c)


## 2.1. System Description

Figure1 shows the schematic diagram of the cam-follower system. The follower makes a fixed axis rotation around O. The cam rotates around a fixed point O′, R is the radius, and ω is the rotating angular velocity of the cam. The eccentric e is the distance between the geometric center C and the rotation point O′. The angular displacement of the cam is θc=∠CO′x′, and the angle between the follower and the horizontal axis is θf. At the initial time, A is the cam and follower contact point. When the system is rotating, Bf and Bc are the impact points on the follower and cam, respectively. iAj is a local coordinate system, and note δi as the distance BfBc→, which indicates the deformation of contact springs between the follower and cam. Obviously, the size of the δi can be used to determine whether the impact occurs or not, that is, 0<δi indicates the separate state, and 0≥δi indicates the contact state.Figure 1
Schematic diagram of the cam-follower system.


## 2.2. The Equation of Motion

During the operation of the cam, the follower will show two kinds of motion states, namely, the contact state and the separate state, and the topology of the system will change with time.Using the Lagrangian approach, the follower motion equation can be obtained:(1)Jθ¨f+mgl2cosθf=0,separate  state,Jθ¨f+mgl2cosθf=Fnl1+μfsgnvrtFn,contact  state,where J is the follower’s moment of inertia, μf is the friction factor between follower and cam, and Fn is the normal contact force, and the expression is(2)Fnt=−Knδi−cδ˙i,where Kn and c represent the equivalent stiffness and damping coefficient of the normal contact, respectively.Using the pure damping model which was proposed by Padmanabhan [20], the given damping coefficient is(3)c=βKnδi,where β is the impact damping factor.


## 2.3. Impact Law

In this model, the impact time is very short, the displacement of the various parts is constant, and the speed changes before and after the impact. By ignoring the plastic characteristic, the oblique impact problem studied here only involves the elastic impact. Considering the impact recovery coefficient, the impact process is divided into two stages: compression and recovery. The velocity relationship before and after the impact is(4)θ˙fti+=−ξθ˙fti−+1+ξ−x˙Bcsinθf+y˙Bccosθfl1,where ξ is the coefficient of impact recovery and x˙Bc and y˙Bc are given by the following expressions:(5)x˙Bc=−eθ˙csinθc−Rθ˙fccosθfc,y˙Bc=eθ˙ccosθc−Rθ˙fcsinθfc.


## 2.4. Normal Contact Stiffness Fractal Model

The contact stiffness is related to the material properties of the two contacting objects and the geometry of the contact surface. For a specific mechanical system, the material properties of the components have been determined, and only the contact surface geometry may affect the contact stiffness coefficient.The actual mechanical interface consists of two rough surfaces. According to the results of Reference [21], the contact area can be reduced to a rigid plane in contact with a rough surface. Figure 2 shows a schematic diagram of the contact between a microconvex body and a rigid plane. When the contact occurs, the real area of the rough surface accounts for only a small part of the nominal contact area. Therefore, the contact problem of the mechanical joint surface can be handled with the contact of the convex body. When the surface is stressed, the microconvex body will occur elastic or plastic deformation. We focus on the elastic deformation of the local components.Figure 2
Schematic diagram of the contact between a microconvex body and a rigid plane. (a) Unloaded state. (b) Loaded state.
(a)
(b)The single microconvex body on the rough surface is equivalent to the sphere, andR is recorded as its equivalent radius of curvature.When the normal loadp is applied to the contact point, δ is the normal contact deformation. The radius of the contact area is R. According to Reference [20], the relationship between p and δ is(6)p=43ER1/2δ3/2,where E is the modulus of elasticity of the two contact materials:(7)1E=1−ϑf2Yf+1−ϑc2Yc,where ϑ and Y represent Poisson’s ratio and Young’s modulus of the material, respectively. Subscripts f and c represent the follower and cam.According to equation (6), the normal contact stiffness of a single microconvex body contact with a rigid plane is(8)kn=2ER1/2δ1/2.According to the geometric relationship between the microconvex deformation and the typical value of the fractal roughness parameterG in Reference [22], R≫δ, and the expression of the cross-sectional area a′ of the microconvex body is(9)a′=2πRδ.Substituting equation (9) into (8), we can obtain(10)kn=2Ea′2π.According to the contact fractal theory [23], the distribution function of the contact point of the microcontact area a′ is(11)na′=D2al′D/2a′D/2+1,(12)al′=2−DDAr,where al′ is the actual contact area of the maximum contact point, Ar is the real contact area of the rough surface, and D is the fractal dimension of rough surface which reflects the irregularity and complexity of surface profile on all scales. Considering the two-dimensional fractal problem, the range D∈1,2. According to the Hertz contact theory, when contact occurs, the contact area of two objects is elliptical, and the real contact area Ar is about 10−5m2.Assume that the rough surface is isotropic in the microtopography. According to the contact fractal theory, the microconvex body contact deformation belongs to the plastic deformation whena′<ac′, and the contact deformation belongs to the elastic deformation when a′>ac′. ac′ is determined by Reference [24] as(13)ac′=G2kφ/22/D−1,where G is the fractal roughness parameter of the rough surface. Larger G means rougher contact surface. k=H/σy is a factor related to the softer material hardness H and the yield strength σy, and k is generally taken as 2.8. φ=σy/E is a factor related to σy and the equivalent elastic modulus E.Thus the normal contact stiffness of the joint surface is(14)Kn=∫ac′al′knna′da′.Substituting Equations (10) and (11) into (14), we can obtain(15)Kn=2ED2π1−Dal′D/2al′1−D/2−ac′1−D/2.According to the relationshipa′=2a, we have(16)al′=2al,ac′=2ac.Substituting Equation (16) into (15), the normal contact stiffness can be obtained(17)Kn=2EDπ1−DalD/2al1−D/2−ac1−D/2.For a definite rough surface,ac is constant, and the normal stiffness depends on the real contact area al of the maximum contact point.


## 2.5. Contact Line Length with Random Normal Distribution

Using the Hertz contact theory in Reference [24], the equivalent stiffness coefficient is given by(18)Kn=π4Elλ,where lλ represents the length of the contact line.Assume that the length of the contact line obeys the law of normal distribution, that is, the contact line length is controlled by the position parameterμ and the scale parameter σ. μ and σ, respectively, determine the position and amplitude of the normal distribution. Figure 3 shows the numerical simulation of the normal contact line length with random distribution. It is found that the amplitude range of the contact line length increases as the scale parameter σ increases. Obviously, the position parameter μ and the scale parameter σ affect the dynamic response of the system. σ=0 indicates that the normal contact line length is constant, and σ≠0 indicates that the length is a random parameter.Figure 3
Normally distributed random contact line length. (a)μ=0.01m,σ=0.001m. (b) μ=0.01m,σ=0.002m. (c) μ=0.01m,σ=0.003m.
(a)
(b)
(c)


## 3. Numerical Simulation

Considering the system dynamic equations contain strong nonlinear time-varying parameters, the analytic solutions are difficult to obtain. The numerical iterative method, fourth-order Runge–Kutta, is employed to solve the dynamic equations and calculates the system’s response over 20 seconds. From the numerical simulation results, it is found that the first several periods’ responses are generally affected by the initial integral parameters and we give up these responses to eliminate the transient effect.At the initial time, the follower and cam keep in touch, and the parameter values of the cam-follower oblique impact system are shown in Table1. In the simulation process, the system response is calculated by equation (1), and the follower angular velocity value is calculated by the impact law of equation (4).Table 1
Values of the system parameters.
Parameters
Values
Follower length,l
0.6
m
Follower radius,r
0.005
m
Follower mass,m
0.3676
kg
Cam radius,R
0.09
m
Cam eccentricity,e
0.04
m
Poisson’s ratio,ϑ
0.3
Acceleration of gravity,g
9.81
m
/
s
2
Friction coefficient,μf
0.3
Impact damping factor,β
3.25
s
/
m
Impact recovery coefficient,ξ
0.42
Modulus of elasticity,Y
2.07
∗
10
11
Pa
Cam speed,ω
28.7
rad
/
s


### 3.1. Effect of Contact Line Length

The length of the contact line affects the nonlinear stiffness coefficient of the contact surface. Taking the contact line lengthlλ as the control parameter, lλ varies from 0.001 m to 0.018 m in steps of 0.00025 m, and the global bifurcation diagram of the follower angle position θf is shown in Figure 4. It can be found from the bifurcation diagram that the cam-follower system exhibits complex dynamic behaviors, including periodic motion, quasiperiodic motion, and chaotic motion. The system response changes from the initial chaotic motion first to periodic 2, then to quasiperiod 2, and then finally to chaos.Figure 4
Global bifurcation diagram ofθf.


### 3.2. Effect of Joint Surface Fractal

Takek = 2.8, φ = 0.01, Ar = 10−5m2, D∈1.36,1.56, and G takes 10−9m, 10−10m, and 10−11m, respectively. Figure 5 shows the variation of contact stiffness Kn with fractal parameters G and D. It can be found that Kn increases with the increase of the fractal dimension D and decreases with the increase of the scale parameter G, which indicates that increasing the roughness of the joint surface is benefit to the reduction of the normal contact stiffness.Figure 5
Variation of contact stiffnessKn with G and D.Taking the fractal dimensionD as the control parameter, the global bifurcation diagram of the follower angle position θf for the different scale parameters G is shown in Figure 6. It can be found from Figure 6 that the system response is more sensitive to the fractal parameters D and G. Smaller parameter variety will affect the response state. When G=10−10m, the system first enters the period 2, soon moves into the quasiperiod 2, and finally moves into large amplitude chaos. In addition, G=10−9m corresponds to the larger bandwidth of the relatively stable response interval of the system. Smaller G means smaller joint normal contact stiffness, and it is benefit to the stability of the system response.Figure 6
Global bifurcation diagram ofθf with D for different G. (a) G=10−9m. (b) G=10−10m.
(a)
(b)


### 3.3. Effect of Random Contact Line Length

From the bifurcation diagram of Figure4, it can be found that the system performs quasiperiod 2 motion when lλ=0.008m. Figure 7 shows the system response with different contact line length distributions for σ∈0.001m,0.003m and μ=0.008m. θfc is the angle between the rod and the horizontal direction in the persistent contact. It can be found that as the random deviation σ increases, the system response evolves from quasiperiodic to chaos, which indicates that the contact line length lλ with a normal random distribution can lead the system unsteady.Figure 7
System response with different contact line length distributions (μ=0.008m). (a) Quasiperiod 2 (σ=0.001m; left frame: red represents θf, blue represents θfc). (b) Quasiperiod 2 (σ=0.002m; left frame: red represents θf, blue represents θfc). (c) Chaos (σ=0.004m; left frame: red represents θf, blue represents θfc).
(a)
(b)
(c)In Figure4, when lλ=0.015m, the system response is chaos. Figure 8 shows the system response for different σ when μ=0.015m and lλ=0.015m. It can be seen that all systems exhibit quasiperiod 2 motions, which indicates that the random deviations can effectively suppress the chaos in the system.Figure 8
System response for the differentσ (μ=0.015m). (a) Quasiperiod 2 (σ=0.001m; left frame: red represents θf, blue represents θfc). (b) Quasiperiod 2 (σ=0.002m; left frame: red represents θf, blue represents θfc). (c) Quasiperiod 2 (σ=0.003m; left frame: red represents θf, blue represents θfc).
(a)
(b)
(c)The determinate response is shown in Figure9 when lλ=0.023m, and it is chaos. Figure 10 shows the system response after applying a random deviation when μ=0.023m. It can be found that the system also exhibits quasiperiod 2, which means applying random deviation can effectively suppress the chaos in the system.Figure 9
Chaotic response whenlλ=0.023m (μ=0.023m,σ=0).Figure 10
System response after applying a random deviation (μ=0.023m). (a) Quasiperiod 2 (σ=0.001m; left frame: red represents θf, blue represents θfc). (b) Quasiperiod 2 (σ=0.002m; left frame: red represents θf, blue represents θfc). (c) Quasiperiod 2 (σ=0.003m; left frame: red represents θf, blue represents θfc).
(a)
(b)
(c)


## 3.1. Effect of Contact Line Length

The length of the contact line affects the nonlinear stiffness coefficient of the contact surface. Taking the contact line lengthlλ as the control parameter, lλ varies from 0.001 m to 0.018 m in steps of 0.00025 m, and the global bifurcation diagram of the follower angle position θf is shown in Figure 4. It can be found from the bifurcation diagram that the cam-follower system exhibits complex dynamic behaviors, including periodic motion, quasiperiodic motion, and chaotic motion. The system response changes from the initial chaotic motion first to periodic 2, then to quasiperiod 2, and then finally to chaos.Figure 4
Global bifurcation diagram ofθf.


## 3.2. Effect of Joint Surface Fractal

Takek = 2.8, φ = 0.01, Ar = 10−5m2, D∈1.36,1.56, and G takes 10−9m, 10−10m, and 10−11m, respectively. Figure 5 shows the variation of contact stiffness Kn with fractal parameters G and D. It can be found that Kn increases with the increase of the fractal dimension D and decreases with the increase of the scale parameter G, which indicates that increasing the roughness of the joint surface is benefit to the reduction of the normal contact stiffness.Figure 5
Variation of contact stiffnessKn with G and D.Taking the fractal dimensionD as the control parameter, the global bifurcation diagram of the follower angle position θf for the different scale parameters G is shown in Figure 6. It can be found from Figure 6 that the system response is more sensitive to the fractal parameters D and G. Smaller parameter variety will affect the response state. When G=10−10m, the system first enters the period 2, soon moves into the quasiperiod 2, and finally moves into large amplitude chaos. In addition, G=10−9m corresponds to the larger bandwidth of the relatively stable response interval of the system. Smaller G means smaller joint normal contact stiffness, and it is benefit to the stability of the system response.Figure 6
Global bifurcation diagram ofθf with D for different G. (a) G=10−9m. (b) G=10−10m.
(a)
(b)


## 3.3. Effect of Random Contact Line Length

From the bifurcation diagram of Figure4, it can be found that the system performs quasiperiod 2 motion when lλ=0.008m. Figure 7 shows the system response with different contact line length distributions for σ∈0.001m,0.003m and μ=0.008m. θfc is the angle between the rod and the horizontal direction in the persistent contact. It can be found that as the random deviation σ increases, the system response evolves from quasiperiodic to chaos, which indicates that the contact line length lλ with a normal random distribution can lead the system unsteady.Figure 7
System response with different contact line length distributions (μ=0.008m). (a) Quasiperiod 2 (σ=0.001m; left frame: red represents θf, blue represents θfc). (b) Quasiperiod 2 (σ=0.002m; left frame: red represents θf, blue represents θfc). (c) Chaos (σ=0.004m; left frame: red represents θf, blue represents θfc).
(a)
(b)
(c)In Figure4, when lλ=0.015m, the system response is chaos. Figure 8 shows the system response for different σ when μ=0.015m and lλ=0.015m. It can be seen that all systems exhibit quasiperiod 2 motions, which indicates that the random deviations can effectively suppress the chaos in the system.Figure 8
System response for the differentσ (μ=0.015m). (a) Quasiperiod 2 (σ=0.001m; left frame: red represents θf, blue represents θfc). (b) Quasiperiod 2 (σ=0.002m; left frame: red represents θf, blue represents θfc). (c) Quasiperiod 2 (σ=0.003m; left frame: red represents θf, blue represents θfc).
(a)
(b)
(c)The determinate response is shown in Figure9 when lλ=0.023m, and it is chaos. Figure 10 shows the system response after applying a random deviation when μ=0.023m. It can be found that the system also exhibits quasiperiod 2, which means applying random deviation can effectively suppress the chaos in the system.Figure 9
Chaotic response whenlλ=0.023m (μ=0.023m,σ=0).Figure 10
System response after applying a random deviation (μ=0.023m). (a) Quasiperiod 2 (σ=0.001m; left frame: red represents θf, blue represents θfc). (b) Quasiperiod 2 (σ=0.002m; left frame: red represents θf, blue represents θfc). (c) Quasiperiod 2 (σ=0.003m; left frame: red represents θf, blue represents θfc).
(a)
(b)
(c)


## 4. Conclusions

In this paper, the fractal model of the normal contact stiffness has been established by considering the microtopography of the joint surface. The effects of the fractal dimensionD and the scale parameter G on the dynamic characteristics of the cam-follower oblique impact system have been studied. The effects of the location and scale parameters of the system dynamic characteristics have been analyzed. The main conclusions can be listed as follows:(1)
Increasing the roughness of the joint surface is a benefit to the reduction of the normal contact stiffness.(2)
The system response is more sensitive to the fractal parameters,D and G. The response state changes as a minor parameter varies, and the system response eventually turns into chaos with the increase of D.(3)
The system response changes between quasiperiod and chaos affected by the contact line length with a normal random distribution.



---

*Source: 3956169-2019-01-22.xml*