Research Article  Open Access
Modeling and Nonlinear Response of the CamFollower ObliqueImpact System
Abstract
In order to quickly and accurately analyze the complex behavior of camfollower obliqueimpact 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 timestepping method was employed to solve the linear complementarity problem which transformed by the obliqueimpact 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.
1. Introduction
Camfollower 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 camfollower is a typical oblique impact. In analyses using nonsmooth dynamics, oblique impact of rough bodies in an unsymmetrical configuration can result in selflocking 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 camfollower. 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 camfollower. 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 timestepping method is employed to solve the linear complementarity problem which transforms by the obliqueimpact 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 form , , , and , where is an matrix and x, y, and b are an 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:
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. ObliqueImpact Model
A camfollower with obliqueimpact on the contact point is studied in Figure 1. 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.
In the present work, the normal contact between rigid bodies is characterized by a setvalued 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 . The relative normal gap is nonnegative due to bodies’ impenetrability condition, being the two bodies in contact with each other when . According to Signorini’s condition, the normal contact conditions areEquation (2) represents a complementarity behavior; always zero is the product of the relative normal gap and normal contact force; that is, .
(a)
(b)
The classical Coulomb friction law is another typical example that can be considered as a setvalued force law [9]. For the two conditions sliding and sticking, the tangential forces arewhere is the friction coefficient; is the normal contact force; and is the relative tangential velocity. These two conditions can be summarized by a setvalued force law as [9]
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 Figure 1, 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]where is generalized velocity and here is the follower’s angular velocity ; is the follower’s moment of inertia; is the follower’s external force; is the generalized coefficient in contact point; is the force of the follower, and and are the velocity before and after the impact point, respectively. Notice that is a continuous boundary parameter when there is no impact. is a noncontinuous parameter when there is viscoustoslip transition or change in the direction of relative speed. denotes normal direction and denotes tangential direction.
The relationship between , and contact point’s relative speeds is
In order to get the unknown parameters and , we need to calculate the distance function in Figure 1:
According to the results of [11],Therefore,
There are three equations for the camfollower impact system. In order to quickly and accurately analyze the complex behavior of camfollower obliqueimpact 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 and (7) is multiplied by the atoms measure . Then add the two equations:This equation can be abbreviated aswhere is velocity measure. It contains two parts: one is the continuous part, Lebesgue measure ; the other is the noncontinuous part, Dirac measure. When the camfollower is in the nonimpact state, the velocity measure . The impact measure which relates to the Lagrange multiplier also contains two parts: . The separation, impact, and contact equations are expressed as the unified form in (13) [10].
Moreau timestepping method [10] is employed to discretize (13). We can getwhere and represent starting and ending time. is the middle time of and . For the normal impact parameter , apply the Upr function:where .
In order to get the expectation LCP condition, we need to separate Sgn function to two Upr functions. For the tangential impact parameter , we can get
For the sake of brevity, we define and , where “” means “is defined as.” The second and third expressions of (16) are and . The LCP condition is
The full LCP equation of follower can be expressed as
From the solution of (17) and (18), we can get and . Then the generalized velocity can be calculated byand the generalized displacement can be gotten by the definition of Moreau time timestepping 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 Table 1.

With the increasing of the cam rotational speed, the state of the camfollower 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 Figure 3.
(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 Figure 4 when 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.
(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 Figure 5, the motion of the follower is periodic 1 or the orbit is .
(a)
(b)
When the cam speed goes higher, as shown in Figures 6 and 7, the responses are . From the time history of these responses, the common point is that there is only 1 impact in each cam period.
(a)
(b)
(a)
(b)
There also exist motions. The rotation speeds of these motions are between two chaos areas in Figure 3. Figure 8 shows the 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.
(a)
(b)
5.3. Chaos Motion
Generally, for vibroimpact system, there are two ways for Poincaré section selection: to select the fixed phase section of trajectory as Poincaré and to select the collision surface as Poincaré section. In this paper, the collision surface is used: represents the collision point, and the mirror collision can be excluded by .
Figure 9 shows the response when rad/s. The collision does not follow the periodic variation and became chaotic. Figure 10 shows the response when rad/s. Compared with Figure 9, the chaos response shows more characteristics of slip phenomenon with 21.7 rad/s.
(a)
(b)
(c)
(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.
Notations
:  Radius of the follower 
:  Bear fix point of the follower 
:  Length of the follower 
:  Angular position of follower with counterclockwise sense of rotation 
:  Angular position of cam with counterclockwise sense of rotation 
:  Eccentricity of the cam 
:  Rotation speed of the cam 
:  Rotation center of the cam 
:  Geometric center of the cam 
:  Initial angular position of the cam 
:  Displacement between cam and follower 
:  Displacement between follower’s fix point and the point which is nearest to the cam. 
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grants nos. 11272257 and 11672251) and NPU Aoxiang New Star.
References
 R. Alzate, M. Di Bernardo, U. Montanaro, and S. Santini, “Experimental and numerical verification of bifurcations and chaos in camfollower impacting systems,” Nonlinear Dynamics, vol. 50, no. 3, pp. 409–429, 2007. View at: Publisher Site  Google Scholar
 Y. Shen and W. J. Stronge, “Painleve paradox during oblique impact with friction,” European Journal of Mechanics. A. Solids, vol. 30, no. 4, pp. 457–467, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 R. Alzate, P. T. Piiroinen, and M. Di Bernardo, “From complete to incomplete chattering: a novel route to chaos in impacting camfollower systems,” International Journal of Bifurcation and Chaos, vol. 22, no. 5, Article ID 1250102, 2012. View at: Publisher Site  Google Scholar
 S. Sundar, J. T. Dreyer, and R. Singh, “Rotational sliding contact dynamics in a nonlinear camfollower system as excited by a periodic motion,” Journal of Sound and Vibration, vol. 332, no. 18, pp. 4280–4295, 2013. View at: Publisher Site  Google Scholar
 H. Ding, L.Q. Chen, and S.P. Yang, “Convergence of Galerkin truncation for dynamic response of finite beams on nonlinear foundations under a moving load,” Journal of Sound and Vibration, vol. 331, no. 10, pp. 2426–2442, 2012. View at: Publisher Site  Google Scholar
 H. Ding and L.Q. Chen, “Galerkin methods for natural frequencies of highspeed axially moving beams,” Journal of Sound and Vibration, vol. 329, no. 17, pp. 3484–3494, 2010. View at: Publisher Site  Google Scholar
 T. Zhang, H. G. Li, Z. Y. Zhong, and G. P. Cai, “Hysteresis model and adaptive vibration suppression for a smart beam with time delay,” Journal of Sound and Vibration, vol. 358, pp. 35–47, 2015. View at: Publisher Site  Google Scholar
 M. B. Meingast, M. Legrand, and C. Pierre, “A linear complementarity problem formulation for periodic solutions to unilateral contact problems,” International Journal of NonLinear Mechanics, vol. 66, pp. 18–27, 2014. View at: Publisher Site  Google Scholar
 R. I. Leine and C. Glocker, “A setvalued force law for spatial CoulombContensou friction,” European Journal of Mechanics—A/Solids, vol. 22, no. 2, pp. 193–216, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 P. Flores, R. Leine, and C. Glocker, “Application of the nonsmooth dynamics approach to model and analysis of the contactimpact events in camfollower systems,” Nonlinear Dynamics, vol. 69, no. 4, pp. 2117–2133, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 C. Glocker, “On frictionless impact models in rigidbody systems,” Philosophical Transactions of the Royal Society of London B: Biological Sciences, vol. 359, no. 1789, pp. 2385–2404, 2001. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2016 YongFeng Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.