Published: 15 November 2017

Dynamic analysis of reciprocating compressor with clearance and subsidence

Shungen Xiao1
Shulin Liu2
Shouguo Cheng3
Xiaogang Xue4
Mengmeng Song5
Xin Sun6
1, 2, 3, 4, 6School of Mechatronical Engineering and Automation, Shanghai University, Shanghai, People’s Republic of China
3Department of Mechanical and Electrical Engineering, Jiangyin Polytechnic College, Jiangsu, People’s Republic of China
5, 1College of Information, Mechanical and Electrical Engineering, Ningde Normal University, Ningde, People’s Republic of China
Corresponding Author:
Shulin Liu
Views 958
Reads 406
Downloads 1717

Abstract

In this paper, the reciprocating compressor with clearance and subsidence of the crosshead is studied, considering the varying cylinder pressure. To investigate the nonlinear dynamical behaviors of the reciprocating compressor precisely, a novel dynamic model of the reciprocating compressor of single cylinder with clearance and subsidence for the varying cylinder pressure is established based on the Lagrangian approach. It is demonstrated that different clearance sizes, subsidence sizes and crankshaft angular velocities can affect the dynamic response of the crosshead in the horizontal direction, respectively. The analysis results of dynamic response illustrate that with the increase of clearance sizes and subsidence sizes, the influence of crosshead displacement, velocity and acceleration increase, but the effect of clearance is much stronger than subsidence. Fast Fourier transformation is also applied to analyze the frequency spectrum of the system response. Moreover, the stability of the motion is studied using phase trajectories, maximum Lyapunov exponents and Poincaré sections. Numerical solution results show that clearance sizes, subsidence sizes and crankshaft angular velocities change the dynamical behavior of the system. Through calculation of the maximum Lyapunov exponents and Poincaré portraits, the results reveal that the reciprocating compressor system is characterized by chaotic behavior with the faults of clearance and subsidence. Furthermore, with the increase of the clearance sizes and the crankshaft angular velocities, the shapes of chaotic phase trajectories have changed greatly, and varied inconspicuously with the increase of subsidence sizes in reciprocating compressor system.

1. Introduction

Reciprocating compressors are one of the most popular machines used in petroleum and chemical production processes, such as gas compression, petroleum transportation and natural gas transportation [1, 2]. In practice, with reciprocating compressor working on a period of time, as the result of manufacture tolerance, material deformation, and wear, some clearances in its joints and subsidence of crosshead commonly exist, which are inevitable. In the case of oversized joint clearances and subsidence, contact forces generate impulsive effect, which causes increased vibration and noise, and reduces system reliability, stability, life and precision. So, clearance and subsidence play a significant role in the prediction of kinematic and dynamic behavior of reciprocating compressor [3].

In the past few decades, most studies have focused on the fault diagnosis method of reciprocating compressor with signal processing. For example, aiming at the bearing clearance fault of reciprocating compressor, a series of fault detection approaches have been used extensively in the literatures [4-6]. In addition, there are also many other studies on the fault diagnosis of reciprocating compressor valves [7-10]. However, few researchers studied the failure mechanism of reciprocating compressor with joint clearances fault and subsidence fault of crosshead, and so on. Zhao et al. [11] studied the dynamic behavior of a reciprocating compressor transmission mechanism with joint clearance. In their study, taking cylinder pressure, amplitude of clearance, crank speed and flexibility of connecting rod as influence factors, the results showed that these factors play a significant role in dynamic behavior of mechanisms. Zhao et al. [12] also presented a parameter optimization approach for planar joint clearance model and its application for dynamics simulation of reciprocating compressor. The dynamics response experimental test verified the effectiveness of this application.

Reciprocating compressor of single cylinder is a slider-crank mechanism, which is shown in Fig. 1. The crankshaft rotates one revolution, driving the connecting rod, the crosshead and the piston to reciprocate. The cylinder can realize the four processes of expansion, suction, compression and exhaust. Although working environment of reciprocating compressor has its particularity, kinematics and dynamics analysis of slider-crank mechanism with joint clearance can be used for reference. There are several publications on the dynamical analysis of slider-crank mechanisms with clearance: with regard to continuous contact force model with clearance, Earles [13], Zhang [14] and Erkaya [15] introduced a model based on permanent contact condition to study clearances in the joints, respectively. In the model, clearance is replaced by a massless virtual link that connects the journal center to the bearing center. Moreover, Feng [16] et al. performed an optimization method to control the inertia forces by redistribution of masses of the moving links in planar mechanisms, in the presence of clearances at joints. In their study, clearance is also modeled with a virtual massless link that connects the journal center to the bearing center.

For contact-separation model with clearance, Dubowsky and Freudenstein [17, 18] formulated an impact pair model to predict the dynamic response of an elastic mechanical joint with clearance. In fact, in the recent decades, there are a lot of researches [19-25] using the impact model with clearance. The contact or separation between the bushing and pin is determined by the relative eccentricity between joint elements. If the eccentricity distance is greater than the radial clearance of a joint, the bushing and pin will contact each other in the direction of the eccentricity. Conversely, if the eccentricity is less than the radial clearance of the joint, the bushing and pin will remain separated from each other [26].

For chaotic behavior with joint clearance case, Farahan [27] analyzed the nonlinear dynamic behavior of a four-bar mechanism with joint clearance at the connection between the coupler and rocker. In their study, Numerical simulations display both periodic and chaotic motions in the system behavior. Wang [28] et al. carried out a study on chaotic control of a flexible multibody system with uncertain joint clearance. In their research, with the purpose to keep the continuous contact of the clearance joint, a modified extended delayed feedback control is used to stabilize the chaotic motion of the flexible multibody system.

Fig. 1Schematic diagrams of reciprocating compressor of single cylinder

Schematic diagrams of reciprocating compressor of single cylinder

Based on the research results of the previous description, the main purpose of this study is to conduct a study based on dynamic response of single cylinder reciprocating compressor with clearance joints and subsidence of the crosshead. The paper is organized as follows: The dynamic model of single cylinder reciprocating compressor with clearance and subsidence is established in Section 2. The effects of joint clearance size, subsidence size and angular velocity are discussed in Section 3, respectively. Phase trajectories of the reciprocating compressor crosshead are analyzed, and the chaotic behavior of the motion is studied using Lyapunov exponents and Poincaré sections with the change of clearance size, subsidence size and crankshaft angular velocity under varying cylinder pressure in Section 4. Dynamic response and chaotic behavior of the reciprocating compressor are summarized in Section 5. In addition, the conclusions of this paper are also given in Section 5.

2. System model of reciprocating compressor with clearance and subsidence

After the reciprocating compressor works for a period of time, the revolute joints clearance gradually increases and the crosshead may generate subsidence due to wear. In order to study the combined effects of clearance and subsidence on the motion of the crosshead, we build the dynamical mathematical model of the reciprocating compressor system.

2.1. Model of cylinder pressure

With the crankshaft of reciprocating compressor rotating a circle, the cylinder is accordingly carried out a working cycle, including expansion, suction, compression and exhaust processes, as shown in Fig. 2. In a circle, when the crankshaft is located at zero angle, the piston is located at right end with the inlet and exhaust vales closing at the same time. As the crankshaft turns clockwise, the crosshead and piston move from right to left, working volume that is located at the right side of the piston gradually increases, cavity gas gradually expands and cylinder pressure gradually decreases, which is expansion process. When the cylinder internal pressure decreases to slightly less than the external pressure of the cavity, the inlet valve is opened, until the crosshead and piston move to the far left, which is suction process. In the suction process, cylinder pressure is almost constant. The suction process is accomplished at the far left. The compression and exhaust processes are mostly the opposite of expansion and suction. Combining with the indicator diagram of reciprocating compressor, cylinder pressure can be expressed using the following equation:

1
P=200000μ,

where P is the cylinder pressure, μ is a cylinder pressure coefficient given by:

2
μ=sinπ2-θ1, 0°+2nπθ17π16+2nπ n is integer Expansion,sinπ16, 7π16+2nπθ1π+2nπ n is integer Suction,-sinθ1+π16, π+2nπθ123π16+2nπ n is integer Compression,1, 23π16+2nπθ12π+2nπ n is integer Exhaust.

2.2. Equations of motion

In this subsection, motion equations are obtained using Lagrangian approach for the reciprocating compressor mechanism with subsidence of the crosshead and revolute clearance joint between the crankshaft and the connecting rod. In normal working condition, the schematic of kinematic diagram of reciprocating compressor is showed in Fig. 3(a). It is assumed that there are both a clearance in B position between the crankshaft and the connecting rod, and a subsidence of the crosshead produces because of wear.

Fig. 2Expansion, suction, compression and exhaust processes

Expansion, suction, compression and exhaust processes

Fig. 3Schematic of kinematic diagram of reciprocating compressor

Schematic of kinematic diagram of reciprocating compressor

a) Schematic of reciprocating compressor in normal working condition

Schematic of kinematic diagram of reciprocating compressor

b) Hypothetical rod in revolute pair with clearance

Schematic of kinematic diagram of reciprocating compressor

c) Schematic of reciprocating compressor with clearance in B position and subsidence of the crosshead

According to the positional relationship shown in Fig. 3, the displacement equations of crankshaft, connecting rod and crosshead in centroid position can be written as:

3
x1 y1x2 y2x3 y3=rs 0 0r ls rcr l rccosθ1 sinθ1cosθ2 sinθ2cosθc sinθc.

In which xi and yi i=1,2,3 are the centroid position of crankshaft, connecting rod and crosshead, respectively. θ1, θc and θ2 are the angle of the crankshaft, the clearance and connecting rod with the x-axis, respectively. r is the length of crankshaft, l is the length of connecting rod and rc is clearance radius in B position. rs is the length of A joint to the centroid position of crankshaft and ls is the length of B joint to the centroid position of connecting rod.

According to the geometric relationship of the crank-slider mechanism, θ2 can be calculated using:

4
θ2=-sin-1d+rsinθ1+rcsinθcl.

In this expression, d is the depth of crosshead subsidence, namely, d=-y3. It is clear that θ2 is a function of θ1 and θc. Therefore, θ˙2 is given by:

5
θ˙2=θ2θ1θ˙1+θ2θcθ˙c.

Obviously, the Eq. (4) can become:

6
rsinθ1+rcsinθc+lsinθ2+d=0.

In Eq. (5), θ2/θ1 and θ2/θc are calculated by taking θ1 and θc partial derivatives of Eq. (6), respectively, as follows:

7
rcosθ1+lcosθ2θ2θ1=0,
8
rccosθc+lcosθ2θ2θc=0.

θ2/θ1 and θ2/θc are obtained in matrix form as:

9
θ2θ1θ2θc=-1lcosθ2rcosθ1rccosθc.

By taking time derivatives of xi and yi i=1,2,3, the velocity equations of crankshaft, connecting rod and crosshead in centroid position can be expressed as:

10
x˙iy˙i=θ˙1xiθ1yiθ1+θ˙cxiθcyiθc.

In Eq. (10), xi/θ1, yi/θ1, xi/θc and yi/θc are calculated by taking θ1 and θc partial derivatives of Eq. (3) as follows, respectively:

11
x1θ1y1θ1=rs-sinθ1cosθ1,
12
x2θ1y2θ1=r lsθ2θ1rlsθ2θ1-sinθ1cosθ1-sinθ2cosθ2,
13
x3θ1y3θ1=r lθ2θ1rlθ2θ1-sinθ1cosθ1-sinθ2cosθ2,
14
x1θcy1θc=00,
15
x2θcy2θc=rc lsθ2θcrclsθ2θc-sinθccosθc-sinθ2cosθ2,
16
x3θcy3θc=rc lθ2θcrclθ2θc-sinθccosθc-sinθ2cosθ2.

These Eqs. (11) to (16) can be expressed into a matrix as follow:

17
x1θ1y1θ1x2θ1 y2θ1x3θ1y3θ1x2θcy2θcx3θcy3θc=rs00rlsθ2θ10rlθ2θ100lsθ2θcrc0lθ2θcrc-sinθ1cosθ1-sinθ2cosθ2-sinθccosθc.

By taking time derivative of Eq. (5), the angular velocity equations of connecting rod can be given by:

18
θ¨2=2θ2θ12θ˙12+22θ2θ1θcθ˙1θ˙c+2θ2θc2θ˙12+θ¨cθ2θc+θ¨1θ2θ1.

In Eq. (18), 2θ2/θ12, 2θ2/θC2 and 2θ2/θcθ1 are calculated by taking θ1 and θc partial derivatives of Eq. (9), respectively, as follows:

19
2θ2θ12=rsinθ1+lsinθ2θ2θ12lcosθ2,
20
2θ2θ12=r c lsinθc(cosθ2)2+(rccosθc)2sinθ2l2(cosθc)3,
21
2θ2θcθ1=r.rc.sinθ2.cosθ1.cosθcl2(cosθ2)2.

The acceleration equations of crankshaft, connecting rod and crosshead in centroid position can be solved by taking time derivative of Eq. (10) as follow:

22
x¨iy¨i=θ˙122xiθ122yiθ12+2θ˙1θ˙c2xiθ1θc2yiθ1θc+θ¨cxiθcyiθc+θ˙c22xiθc22yiθc2, i=1,2,3,

where 2xi/θ12, 2yi/θ12, 2xi/θ1θc, 2yi/θ1θc, 2xi/θc2 and 2yi/θc2 can be determined using the following equations from Eqs. (23) to (31).

2x1/θ12 and 2y1/θ12 are calculated by taking θ1 derivative of Eq. (11) as follow:

23
2x1θ122y1θ12=rs-cosθ1-sinθ1.

2x2/θ12 and 2y2/θ12 are computed by taking θ1 derivative of Eq. (12) as follow:

24
2x2θ122y2θ12=-rls2θ22θ1lsθ2θ12r- ls2θ22θ1lsθ2θ12cosθ1sinθ1sinθ2cosθ2cosθ2sinθ2.

2x3/θ12 and 2y3/θ12 can be deduced by taking θ1 derivative of Eq. (13) as follow:

25
2x3θ122y3θ12=-rl2θ22θ1lθ2θ12r- l2θ22θ1 lθ2θ12cosθ1sinθ1sinθ2cosθ2cosθ2sinθ2.

2x1/θ1θc and 2y1/θ1θc are expressed by taking θc derivative of Eq. (11) as follow:

26
2x1θ1θc2y1θ1θc=00.

2x2/θ1θc and 2y2/θ1θc are defined by taking θc derivative of Eq. (12) as follow:

27
2x2θ1θc2y2θ1θc=-ls2θ2θ1θclsθ2θ1θ2θc- ls2θ2θ1θclsθ2θ1θ2θcsinθ2cosθ2cosθ2sinθ2.

2x3/θ1θc and 2y3/θ1θc are computed by taking θc derivative of Eq. (13) as follow

28
2x3θ1θc2y3θ1θc=-l2θ2θ1θclθ2θ1θ2θc- l2θ2θ1θclθ2θ1θ2θcsinθ2cosθ2cosθ2sinθ2.

2x1/θc2 and 2y1/θc2 are calculated by taking θc derivative of Eq. (14) as follow:

29
2x1θc22y1θc2=00.

2x2/θc2 and 2y2/θc2 can be deduced by taking θc derivative of Eq. (15) as follow:

30
2x2θc22y2θc2=-ls2θ22θclsθ2θc2rc- ls2θ22θclsθ2θc2 rcsinθ2cosθ2cosθ2sinθ2cosθcsinθc.

2x3/θc2 and 2y3/θc2 are calculated by taking θc derivative of Eq. (16) as follow:

31
2x3θc22y3θc2=-l2θ22θclθ2θc2rc- l2θ22θclθ2θc2 rcsinθ2cosθ2cosθ2sinθ2cosθcsinθc.

These Eqs. (23) to (31) can be written into a matrix as follow:

32
2x1θ12 2y1θ122x2θ122y2θ122x2θ1θc2y2θ1θc2x2θc22y2θc22x3θ122y3θ122x3θ1θc2y3θ1θc2x3θc22y3θc2=rs000r±ls2θ2θ12lsθ2θ1200±ls2θ2θ1θclsθ2θ1θ2θc00±ls2θ2θc2lsθ2θ12rcr±l2θ2θ12lθ2θ1200±l2θ2θ1θclθ2θ1θ2θc00±l2θ2θc2lθ2θ12rc×cosθ1sinθ1sinθ2cosθ2cosθ2sinθ2cosθcsinθc.

In Eq. (32), The meaning of the “±” is that “+” and “-” represent the symbol of the first column and the second column of the solved matrix, respectively.

In addition, it is easily obtained as follows:

33
2x1θ1θc=2y1θ1θc=2x1θc2 = 2y1θc2=x 1θc=y1θc
=y3θ1=y3θc=2y3θ12=2y3θ1θc=2y3θc2=0.

2.3. Equations of dynamics

Lagrange motion equation for an unconstrained dynamical system is as follow:

34
ddtEq˙j-Eqj+Uqj=Qj, j=1,2,

where E and U are the system kinetic and potential energies, respectively. Qj is the nonconservative generalized force corresponding to the generalized coordinate qj, and can be calculated as:

35
Qj=i=13Fi*.Viq˙j+Mi*.ωiq˙j.

In Eq. (35), Fi* is the resultant of external force acting at the center of mass. Mi* is the external torque acting on body i. Vi and ωi are the translational and rotational velocity for the mass center of body i, respectively.

Taking the crankshaft, connecting rod and crosshead of the reciprocating compressor as a system, the kinetic energy and potential energy of the system are:

36
E=i=13Ei=i=1312mix˙i2+y˙i2+i=1212Jiθ˙i2,U=i=12migyi.

Substituting Eq. (36) into Eq. (34) leads to:

37
i=12Jiθ¨iθiqj+i=13mix¨ixiqj+y¨iyiqj+gi=13miyiqj=Qj.

Substitute Eqs. (18) and (22) into Eq. (37). Because the crankshaft works for uniform rotation, the angular acceleration of crankshaft equals zero, Namely, θ¨1= 0. Therefore, differential equation about θ¨c can be written as:

38
Aθ¨c=Bθ˙c+Cθ˙c2+D,

where A, B, C and D are expressed as:

39
A=J2θ2θc2+i=13mixiθc2+yiθc2,
40
B=-2ωJ2θ2θc2θ2θcθ1+i=13mixiθc2xiθ1θc+yiθc2yiθ1θc,
41
C=-J2θ2θc2θ2θc2+i=13mixiθc2xiθc2+yiθc2yiθc2,
42
D=Q2-J22θ2θ1θ2θc+i=13mixiθc2xiθ12+yiθc2yiθ12ω2-gi=13miyiθc.

Obviously, Eq. (38) is a nonlinear equation.

It is supposed that the centroid of connecting rod is middle position, that is to say, l=2ls. Moreover, since the clearance size rc is very small, it can be considered the expression rc.x+yy, where x and y are two expressions that are roughly the same order of magnitude. Substituting Eqs. (9), (17), (19)-(21) and (32) into Eqs. (39)-(42), respectively. The expressions of A, B, C and D can be expressed as follows:

43
A=a1+a24cosθccosθ22+a2+a2cosθc(sinθctanθ2-cosθc)
+a3(tanθ2cosθc-sinθc)2,

where ai i=1, 2, 3 are constants, a1=J2rc/l2, a2=m2rc2, a3=m3rc2:

44
B=cosθ1cosθccos4θ2b1sinθ2cosθ1+b2sin(θ2+θc)+(b3+b4)sin(θc-θ2).

In which, bi i=1, 2, 3, 4 are constants, b1=-2ωJ2rc2r/l3,b2=-m2rrc2ω/4l,b3=3b2, b4=-2m3rrc2ω:

45
C=1cos2θ2c1sinθc+c2516cos2θ2sin2θc-14cos2θcsin2θ2+c3sin2θ2-c28sin2θc.

In which, ci i=1,2,3 are constants, c1=-J2rc2/2l2+3c2/16-c3,c2=m2rc2,c3=-m3rc2/2:

46
D=d1(-sinθc+tanθ2cosθc)-d2cosθccos4θ2(lcos2θ2sinθ1+rsinθ2cos2θ1) +d3sin(θc-θ2)cosθ1cosθ2+d4tanθ2sinθ1+rcos2θ1lcos3θ2 +d5(sinθccosθ1+sin(θc-θ1)),

where:

47
d0=ω2rrcd1=Prc, d2=d0J2l2, d3=d0m22+m3,
d4=d0m24+m3+m24sinθc, d5=d0m24.

3. Numerical simulation and dynamic responses analysis

In this previous section, the numerical solution of the clearance angle θc is solved by the Runge-Kutta method according to the Eqs. (38) and (4). Subsequently, the displacement, velocity and acceleration of the crosshead are also obtained accordingly in the horizontal direction. In this section, we discuss the influence of clearance, subsidence and rotation speed of the crankshaft. Taking the 2D12 model of reciprocating compressor as the research object, the simulation characteristics are presented in Table 1. In addition, the centroids of crankshaft and connecting rod are middle position.

Table 1Dimension for the 2D12 model of reciprocating compressor

Parts
Length (m)
Mass (kg)
Moment of inertia (kg.m2)
Crankshaft
0.12
1
0.0012
Connecting rod
0.6
5
0.15
Crosshead
1

3.1. Dynamic response of reciprocating compressor system with different clearance sizes

In this subsection, from the perspective of dynamic response, the effects of the reciprocating compressor with different clearance sizes are investigated under the conditions of without subsidence, crankshaft angular velocity 20 rad/s and variable cylinder load shown in Fig. 2. In this study, the revolute joint between the crankshaft and the linkage is used as the clearance for the research. The displacement, velocity, acceleration and acceleration spectrum of the crosshead without/with 0.1 mm, 0.5 mm and 5 mm clearance sizes are shown in Figs. 4 to 7, respectively. It is noteworthy that the clearance size is exaggerated to 5 mm so as to better investigate the dynamics response of the variational clearance sizes for the reciprocating compressor.

Fig. 4Dynamic responses of crosshead without clearance

Dynamic responses of crosshead without clearance

As can be seen from Figs. 4 to 7, with the increase of clearance sizes, the effects of displacement, velocity, and acceleration increase. It can be observed that the crosshead velocity and acceleration fluctuate more violently than the crosshead displacement. Furthermore, Fig. 4(d), 5(d), 6(d) and 7(d) show the acceleration spectrum using FFT transformation. In Fig. 4(d), it is clear that 3.109 Hz which is the rotational frequency of crankshaft and 6.217 Hz which is two times rotational frequency of crankshaft are displayed without clearance. In Fig. 5(d), In addition to 3.109 Hz and 6.217 Hz, there are very weak other frequency components with 0.1 mm clearance size, and it is distinct that the amplitudes of the other frequency components are lower than 6.217 Hz. In Fig. 6(d), except 3.109 Hz and 6.217 Hz, the other frequency components are obvious with 0.5 mm clearance size, and it is evident that the amplitudes of some other frequency components are higher than 6.217 Hz, but are lower than 3.109 Hz. It can be found that there are 113.5 Hz frequency component and approximate multiplier such as 223.1 Hz and 330.3 Hz in other frequency components. In Fig. 7(d), apart from 3.109 Hz and 6.217 Hz, the other frequency components are extremely obvious with 5mm clearance size, and appear 148.4 Hz frequency component and approximate multiplier such as 287.6 Hz and 448.4 Hz. It is noteworthy that the amplitudes of some other frequency components are higher than 3.109 Hz. Therefore, the results show that besides rotational frequency of crankshaft, the larger the clearance sizes are, the more obvious the other frequency components are, even beyond the amplitude of the rotational frequency components of the crankshaft.

Fig. 5Dynamic responses of crosshead with 0.1 mm clearance

Dynamic responses of crosshead with 0.1 mm clearance

Fig. 6Dynamic responses of crosshead with 0.5 mm clearance

Dynamic responses of crosshead with 0.5 mm clearance

The comparison of response curves of displacement, velocity and acceleration is demonstrated in Figs. 8 to 10, respectively. It can be seen from Figs. 8 to 10 that the variational clearance does not influence the crosshead displacement in a remarkable way, and the crosshead displacement and velocity of the reciprocating compressor without clearance are generally consistent with those of the reciprocating compressor with a different clearance size. When the clearance sizes increase from 0.1 mm to 1 mm, the corresponding maximal deviation values of the displacement increase from the 0.1 mm to 1.0 mm at the low dead point, and the maximal deviation values of velocity increase from 0.002 m/s to 0.042 m/s. In vivid contrast, the crosshead acceleration is dramatically influenced by the clearance sizes and the maximal peak values of crosshead acceleration increase from 77.68 m/s2 to 402.6 m/s2.

When the clearance size is exaggerated to 5 mm, the corresponding maximal deviation value of the displacement increases to 4.9 mm at the low dead point, and the maximal deviation value of velocity increases to 0.959 m/s, and the maximal peak value of crosshead acceleration increases to 1164 m/s2. Distinctly, 5 mm clearance size creates a large vibration response to the crosshead, which is likely to cause the reciprocating compressor system to collapse.

Fig. 7Dynamic responses of crosshead with 5 mm clearance

Dynamic responses of crosshead with 5 mm clearance

Fig. 8Comparison of crosshead displacement with different clearance sizes

Comparison of crosshead displacement with different clearance sizes

Fig. 9Comparison of crosshead velocity with different clearance sizes

Comparison of crosshead velocity with different clearance sizes

Fig. 10Comparison of crosshead acceleration with different clearance sizes

Comparison of crosshead acceleration with different clearance sizes

3.2. Influence of the crosshead subsidence

As reciprocating compressor works for a period of time, there may be wear between the crosshead and the side of the cylinder, which may cause the crosshead to sink. That is to say, the crosshead produces decentration. Obviously, the crosshead subsidence will exacerbate the vibration of the crosshead and transfer it to the reciprocating compressor body. In this subsection, the influences of subsidence sizes on the dynamic behavior of reciprocating compressor with subsidence of the crosshead are discussed under the conditions of clearance size 1mm, crankshaft angular velocity 20 rad/s and variable cylinder load shown in Fig. 2. The displacement, velocity, and acceleration of the crosshead without/with 0.1 mm, 1 mm and 5 mm subsidence sizes are shown in Figs. 11 to 13, respectively. It is noteworthy that in order to better research the dynamic response of the variational subsidence sizes to the reciprocating compressor, the subsidence size is also exaggerated to 5 mm.

Fig. 11Comparison of crosshead displacement with different subsidence sizes

Comparison of crosshead displacement with different subsidence sizes

As can be seen from Figs. 11 to 13, the existence of the subsidence brings a very weak impact for the crosshead displacement, and does not produce strongly effect on the crosshead velocity and acceleration with different subsidence sizes. With the increase in subsidence sizes, as the subsidence sizes increase from 0 mm to 1 mm, the corresponding maximal deviation values of the displacement increase from 0 mm to 0.4 mm at the low dead point, and the maximal deviation values of velocity increase from 0 m/s to 0.059 m/s, and the maximal peak values of crosshead acceleration increase from 402.6 m/s2 to 548.1 m/s2.

When the subsidence size is exaggerated to 5 mm, the corresponding maximal deviation value of the displacement increases to 0.7 mm at the low dead point, and the maximal deviation value of velocity increases to 0.312 m/s, and the maximal peak value of crosshead acceleration increases to 772.4 m/s2. Compared with the dynamic response of without subsidence size, the 5 mm subsidence size increases the slight fluctuation. Thus, the varying subsidence sizes do not have a significant impact on the reciprocating compressor system. The reason may be that in this paper, the flexibility of the connecting rod between the crosshead and the piston is not considered. So, it reduces the impact of the crosshead on the reciprocating compressor.

Distinctly, compared with the effect of the clearance sizes, the oscillating range of crosshead displacement, velocity and acceleration are smaller with subsidence sizes. That is to say, the influence of the clearance is much stronger than the subsidence of the crosshead in the horizontal direction.

Fig. 12Comparison of crosshead velocity with different subsidence sizes

Comparison of crosshead velocity with different subsidence sizes

Fig. 13Comparison of crosshead acceleration with different subsidence sizes

Comparison of crosshead acceleration with different subsidence sizes

3.3. Influence of the input crankshaft angular velocity

In this subsection, the influence of the input crankshaft speed on the dynamic behavior of the reciprocating compressor is reported. In addition, in order to fully understand the dynamic behavior of the crankshaft at the low speed, medium speed and high speed, dynamic responses of the crosshead from 10 rad/s to 150 rad/s are analyzed under the conditions of subsidence 1 mm, clearance size 1 mm and variable cylinder load shown in Fig. 2, as showed in Fig. 14 to 16, respectively.

It is can be seen from Figs. 14 to 16, the changes of the crankshaft angular velocity do not influence the crosshead displacement in a conspicuous manner, but produce evident effects for the crosshead velocity and acceleration. When the crankshaft angular velocities increase from 10 rad/s to 150 rad/s, the corresponding maximal values of the displacement increase from 480.6 mm to 481.1 mm at the low dead point, and the maximal values of the velocity increase from 1.267 m/s to 20.14 m/s, and the maximal peak values of crosshead acceleration increase from 172.7 m/s2 to 16010 m/s2. In addition, as can be seen from Fig. 14(d), Fig. 15(d) and Fig. 16(d), with the increase in the crankshaft angular velocity, besides rotational frequency of crankshaft, the lower the crankshaft angular velocities are, the more obvious the other frequency components are, even beyond the amplitude of the rotational frequency components.

Fig. 14Dynamic responses of crosshead with the crankshaft speed 10 rad/s

Dynamic responses of crosshead with the crankshaft speed 10 rad/s

Fig. 15Dynamic responses of crosshead with the crankshaft angular velocity 70 rad/s

Dynamic responses of crosshead with the crankshaft angular velocity 70 rad/s

When the crankshaft angular velocity is 10 rad/s, except 1.554 Hz which is the rotational frequency of crankshaft and 3.109 Hz which is the double rotational frequency, there are still some other remarkable frequency components that almost submerge the rotational frequency of crankshaft. One can observe 167.9 Hz frequency component and approximate multiplier such as 337.3 Hz in some other frequency components. When the crankshaft angular velocity is 70 rad/s, in addition to the rotational frequency 10.88 Hz and double rotational frequency 22.44 Hz, there are very weak other frequency components, and it is distinct that the amplitudes of the other frequency components are lower than 22.44 Hz. When the crankshaft angular velocity is 150 rad/s, the rotational frequency 24.04 Hz and double rotational frequency 48.04 Hz are more significant than the other frequency components which can be ignored.

Fig. 16Dynamic responses of crosshead with the crankshaft angular velocity 150 rad/s

Dynamic responses of crosshead with the crankshaft angular velocity 150 rad/s

4. Chaos analysis

Data obtained from a deterministic system can be classified as either periodic or non-periodic data [29]. Non-periodic data may correspond to a quasi-periodic, transient or chaotic motion. Lyapunov exponents produce a method that can distinguish the chaotic and non-chaotic (periodic or quasi-periodic) behavior. Periodic attractors show only negative and zero exponents which indicate convergence to a predictable motion, whereas there is at least one positive exponent for a chaotic system [30, 31]. Hence, one needs to determine the sign of Lyapunov exponents to characterize the behavior of a system. In this section, nonlinear dynamic analysis of the reciprocating compressor crosshead is investigated. We discuss the chaotic and non-chaotic behavior of the crosshead using phase trajectories, Lyapunov exponents and Poincaré section under the different clearance sizes, subsidence sizes and crankshaft angular velocities.

Fig. 17 shows the phase trajectories of displacement-velocity and velocity-acceleration for the crosshead when the reciprocating compressor has neither clearance nor subsidence. Obviously, under such ideal conditions reciprocating compressor has regular periodic motion.

Fig. 17The phase trajectories of displacement-velocity and velocity-acceleration without clearance and subsidence

The phase trajectories of displacement-velocity  and velocity-acceleration without clearance and subsidence

a)

The phase trajectories of displacement-velocity  and velocity-acceleration without clearance and subsidence

b)

Fig. 18 and Fig. 19 show that the phase trajectories of displacement-velocity and |velocity-acceleration for the crosshead under the conditions of subsidence size 1 mm, crankshaft angular velocity 20 rad/s and variable cylinder load with 0.1 mm, 0.5 mm,1 mm and 5 mm clearance sizes, respectively. It is noteworthy that in order to better investigate the dynamic behavior of the varying clearance size to the reciprocating compressor, the clearance size is exaggerated to 5 mm.

It is can be seen from Fig. 18 and Fig. 19, the change of the clearance sizes produces evident effect for the shape of phase trajectories. With the increase of the clearance sizes, the phase trajectories of velocity–acceleration are more and more like the -character. Fig. 24(a) shows all calculated Lyapunov exponents by Wolf [32] method for the crosshead acceleration with 0.1 mm, 0.5 mm, 1 mm and 5 mm clearance sizes. The largest Lyapunov exponents λi have been calculated for clearance sizes rc1= 0.1 mm, rc2= 0.5 mm, rc3= 1 mm and rc4= 5 mm. The values of the largest Lyapunov exponents are λ1= 1.57 to the corresponding rc1= 0.1 mm, λ2= 0.68 to the corresponding rc2= 0.5 mm, λ3= 0.62 to the corresponding rc3= 1 mm, and λ4= 0.59 to the corresponding rc4= 5 mm. One can see a decrease in values of the largest Lyapunov exponents as the clearance sizes increase, but the downward trend is gradually slowed down.

Fig. 18The phase trajectories of displacement-velocity with different clearance sizes: a) 0.1 mm, b) 0.5 mm, c) 1 mm, d) 5 mm

The phase trajectories of displacement-velocity with different clearance sizes:  a) 0.1 mm, b) 0.5 mm, c) 1 mm, d) 5 mm

Fig. 19The phase trajectories of velocity-acceleration with different clearance sizes: a) 0.1 mm, b) 0.5 mm, c) 1 mm, d) 5 mm

The phase trajectories of velocity-acceleration with different clearance sizes:  a) 0.1 mm, b) 0.5 mm, c) 1 mm, d) 5 mm

Fig. 20 and Fig. 21 show that the phase trajectories of displacement-velocity and velocity-acceleration for the crosshead under the conditions of clearance size 1 mm, crankshaft angular velocity 20 rad/s and variable cylinder load with 0 mm, 0.1 mm, 1 mm and 5 mm subsidence sizes, respectively. It is noteworthy that in order to better study the dynamic behavior of the varying subsidence size to the reciprocating compressor, the subsidence size is also exaggerated to 5 mm.

As can be seen from Fig. 20 and Fig. 21, the changes of the subsidence sizes produce weak influence for the shape of phase trajectory. It is visible that the phase trajectory of the crosshead vibration does not close for each revolution. Fig. 24(b) shows all calculated Lyapunov exponents for the crosshead acceleration with 0 mm, 0.1mm, 1 mm and 5 mm subsidence sizes. The largest Lyapunov exponents λi have been computed for subsidence sizes d1= 0 mm, d2= 0.1 mm, d3= 1 mm and d4= 5 mm. The values of the largest Lyapunov exponents are λ1= 0.51 to the corresponding d1= 0 mm, λ2= 0.53 to the corresponding d2= 0.1 mm, λ3= 0.62 to the corresponding d3= 1 mm, and λ4= 0.52 to the corresponding d4= 5 mm.

Fig. 20The phase trajectories of displacement-velocity with different subsidence sizes: a) 0 mm, b) 0.1 mm, c) 1 mm, d) 5 mm

The phase trajectories of displacement-velocity with different subsidence sizes:  a) 0 mm, b) 0.1 mm, c) 1 mm, d) 5 mm

Fig. 21The phase trajectories of velocity-acceleration with different subsidence sizes: a) 0 mm, b) 0.1 mm, c) 1 mm, d) 5 mm

The phase trajectories of velocity-acceleration with different subsidence sizes:  a) 0 mm, b) 0.1 mm, c) 1 mm, d) 5 mm

Fig. 22The phase trajectories of displacement–velocity with different angular velocities: a) 10 rad/s, b) 20 rad/s, c) 70 rad/s, d) 150 rad/s

The phase trajectories of displacement–velocity with different angular velocities: a) 10 rad/s, b) 20 rad/s, c) 70 rad/s, d) 150 rad/s

Fig. 23The phase trajectories of velocity-acceleration with different angular velocities: a) 10 rad/s, b) 20 rad/s, c) 70 rad/s, d) 150 rad/s

The phase trajectories of velocity-acceleration with different angular velocities: a) 10 rad/s, b) 20 rad/s, c) 70 rad/s, d) 150 rad/s

Fig. 22 and Fig. 23 show that the phase trajectories of displacement-velocity and velocity-acceleration for the crosshead under the conditions of clearance size 5 mm, subsidence size 1 mm and variable cylinder load with 10 rad/s, 20 rad/s, 70 rad/s and 150 rad/s crankshaft angular velocities, respectively. It is can be seen from Fig. 22 and Fig. 23, the changes of the crankshaft angular velocity produce remarkable effects on the shape of phase trajectories,the higher the crankshaft angular velocity, the more regular the trajectory. In addition, with the decrease of the crankshaft angular velocities, the phase trajectories of velocity-acceleration are more and more like the -character. Fig. 24(c) shows all calculated Lyapunov exponents for the crosshead acceleration with 10 rad/s, 20 rad/s, 70 rad/s and 150 rad/s crankshaft angular velocities. The largest Lyapunov exponents λi have been computed for crankshaft angular velocities ω11= 10 rad/s, ω12= 20 rad/s, ω13= 70 rad/s and ω14= 150 rad/s. The values of the largest Lyapunov exponents are λ1= 0.97 to the corresponding ω11= 10 rad/s, λ2= 0.59 to the corresponding ω12= 20 rad/s, λ3= 0.23 to the corresponding ω13= 70 rad/s, and λ4= 1.5 to the corresponding ω14= 150 rad/s.

Fig. 24Lyapunov exponents with different clearance sizes, subsidence sizes and angular velocities

Lyapunov exponents with different clearance sizes, subsidence sizes and angular velocities

a) Lyapunov exponents with different clearance sizes

Lyapunov exponents with different clearance sizes, subsidence sizes and angular velocities

b) Lyapunov exponents with different subsidence sizes

Lyapunov exponents with different clearance sizes, subsidence sizes and angular velocities

c) Lyapunov exponents with different angular velocities

Obviously, the signs of the largest Lyapunov exponents are positive for different clearance sizes, subsidence sizes and angular velocities, indicating the exponential separation of nearby trajectories as time evolves, that is to say, the reciprocating compressor system is characterized by chaotic behavior with the faults of clearance and subsidence. Therefore, one can deduce at this point the chaotic behavior of the reciprocating compressor system. The trajectories shown in figure are associated with the motion around the attractor and show the classical characteristics of chaotic motion.

In addition, in order to better analyze the chaos of the reciprocating compressor system, the Poincaré section method is applied in this study. Figs. 25 to 28 are the Poincaré portraits corresponding to Figs. 17, 19, 21 and 23, respectively. Whether there is chaos in a system can be judged by observing the intercept point on the Poincare section. When the Poincare section has and has only one fixed point or a few discrete points, the movement is periodic. When the Poincare section is a closed curve, the movement is quasi-periodic. When the Poincare section is a piece of dense points with a fractal structure, the movement is chaos.

As illustrated in Fig. 25, When the reciprocating compressor system does not have clearance and subsidence, the Poincare sections have only one fixed point, which indicates that the reciprocating compressor is a periodic system in the absence of fault. However, as can be seen from Figs. 26 to 28, obviously, the reciprocating compressor system of the single cylinder with clearance and subsidence has chaotic behavior, which is consistent with the judgment result of the maximum Lyapunov index.

Fig. 25The Poincaré portraits of displacement-velocity and velocity-acceleration without clearance and subsidence

The Poincaré portraits of displacement-velocity  and velocity-acceleration without clearance and subsidence

a) Crosshead displacement (m)

The Poincaré portraits of displacement-velocity  and velocity-acceleration without clearance and subsidence

b) Crosshead velocity (m/s)

Fig. 26The Poincaré portraits for (velocity, acceleration) with different clearance sizes: a) 0.1 mm, b) 0.5 mm, c) 1 mm, d) 5 mm

The Poincaré portraits for (velocity, acceleration) with different clearance sizes:  a) 0.1 mm, b) 0.5 mm, c) 1 mm, d) 5 mm

Fig. 27The Poincaré portraits for (velocity, acceleration) with different subsidence sizes: a) 0 mm, b) 0.1 mm, c) 1 mm, d) 5 mm

The Poincaré portraits for (velocity, acceleration) with different subsidence sizes: a) 0 mm, b) 0.1 mm, c) 1 mm, d) 5 mm

Fig. 28The Poincaré portraits for (velocity, acceleration) with different angular velocities: a) 10 rad/s, b) 20 rad/s, c) 70 rad/s, d) 150 rad/s

The Poincaré portraits for (velocity, acceleration) with different angular velocities: a) 10 rad/s, b) 20 rad/s, c) 70 rad/s, d) 150 rad/s

5. Conclusions

In this work, the nonlinear dynamics of reciprocating compressor of single cylinder with subsidence of the crosshead and clearance between crankshaft and connecting rod are investigated. Dynamics differential equations are established based on the Lagrangian approach. Runge-Kutta method is used to solve the dynamics differential equations by MATLAB software.

Under the condition of the varying cylinder pressure in a cycle, the dynamics responses of reciprocating compressor cylinder are analyzed in the horizontal direction with the changes of clearance sizes, subsidence sizes and crankshaft angular velocities. The numerical simulation results show that the displacement, velocity and acceleration of the crosshead are slightly effected with the increase of the subsidence sizes, but the displacement, velocity and acceleration of the crosshead are significantly affected as the increase of the clearance sizes and crankshaft angular velocities. In addition, the response spectrums of the crosshead acceleration are drawn by using the fast Fourier transform. The results display that besides the rotation frequency of crankshaft, some other frequency components are more and more obvious and even submerge the rotation frequency components in the spectrum with the increase of the clearance sizes. The change in crankshaft angular velocities is just the opposite of the change in the clearance. That is to say, as the crankshaft angular velocities increase, some other frequency components are becoming weaker than the rotation frequency components.

Using Lyapunov exponents and Poincaré section, the stability of the reciprocating compressor system with clearance and subsidence is investigated. The signs of the largest Lyapunov exponents calculated by Wolf method are positive, and the Poincare section is a piece of dense points with a fractal structure for this nonlinear system with clearance and subsidence, indicating the existence of strange attractors and chaos phenomena. In addition, the mechanism exhibits different chaotic motion trajectories with the changes of clearance, subsidence and crankshaft angular velocity, in which the shapes of the phase trajectories significantly change with the variation of clearance and the crankshaft angular velocity, and inconspicuously vary as the change of subsidence.

References

  • Elhaj M., Gub F., Ballb A. D. Numerical simulation and experimental study of a two-stage reciprocating compressor for condition monitoring. Mechanical Systems and Signal Processing, Vol. 22, 2008, p. 374-389.
  • Almasi A. A new study and model for the mechanism of process reciprocating compressors and pumps. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering, Vol. 224, 2010, p. 143-148.
  • Dupac M., Beale D. G. Dynamic analysis of a flexible linkage mechanism with cracks and clearance. Mechanism and Machine Theory, Vol. 45, 2010, p. 1909-1923.
  • Chen G. J., Zou L. Q., Zhao H. Y. An improved local mean decomposition method and its application for fault diagnosis of reciprocating compressor. Journal of Vibroengineering, Vol. 18, Issue 3, 2016, p. 1474-1485.
  • Zhao H. Y., Xu M. Q., Wang J. D. Local mean decomposition based on rational hermite interpolation and its application for fault diagnosis of reciprocating compressor. Journal of Mechanical Engineering, Vol. 51, Issue 1, 2015, p. 83-9.
  • Zhao H. Y., Wang J. D., Hui H. A feature extraction method based on HLMD and MFE for bearing clearance fault of reciprocating compressor. Measurement, Vol. 89, 2016, p. 34-43.
  • Farzaneh-Gord M., Khoshnazar H. Valve fault detection for single-stage reciprocating compressors. Journal of Natural Gas Science and Engineering, Vol. 35, 2016, p. 1239-1248.
  • Pichler K., Lughofer E., Pichler M. Fault detection in reciprocating compressor valves under varying load conditions. Mechanical Systems and Signal Processing, Vol. 70, Issue 71, 2016, p. 104-119.
  • Tran V. T., AlThobiani F., Ball A. An approach to fault diagnosis of reciprocating compressor valves using Teager–Kaiser energy operator and deep belief networks. Expert Systems with Applications, Vol. 41, 2014, p. 4113-4122.
  • Qin Q. Q., Jiang Z. N., Feng K. A novel scheme for fault detection of reciprocating compressor valves based on basis pursuit, wave matching and support vector machine. Measurement, Vol. 45, Issue 5, 2012, p. 897-908.
  • Zhao H. Y., Ji Y. D., Wang J. D. Dynamic analysis of reciprocating compressor transmission mechanism with joint clearance. Fluid Machinery, Vol. 41, Issue 3, 2013, p. 15-33.
  • Zhao H. Y., Xu M. Q., Wang J. D. A parameters optimization method for planar joint clearance model and its application for dynamic s simulation of reciprocating compressor. Journal of Sound and Vibration, Vol. 344, 2015, p. 416-433.
  • Earles S. W. E., Wu C. L. S. Motion analysis of a rigid link mechanism with clearance at a bearing using Lagrangian mechanics and digital computation. Mechanisms, 1973, p. 83-89.
  • Zhang X. C., Zhang X. M. A comparative study of planar 3-RRR and 4-RRR mechanisms with joint clearances. Robotics and Computer-Integrated Manufacturing, Vol. 40, 2016, p. 24-33.
  • Selçuk Erkaya, İbrahim Uzmay Modeling and simulation of joint clearance effects on mechanisms having rigid and flexible links. Journal of Mechanical Science and Technology, Vol. 28, Issue 8, 2014, p. 2979-2986.
  • Feng B., Morita N., Torii T. A new optimization method for dynamic design of planar mechanism with clearances at joints. ASME, Vol. 124, Issue 1, 2002, p. 68-73.
  • Dubowsky S., Freudenstein F. Dynamic analysis of mechanical systems with clearances. Part 1: Formulation of dynamic model. Journal of Engineering for Industry, Vol. 93, Issue 1, 1971, p. 305-309.
  • Dubowsky S., Freudenstein F. Dynamic analysis of mechanical systems with clearances. Part 2: Dynamic response. Journal of Engineering for Industry, Vol. 93, Issue 1, 1971, p. 310-316.
  • Flores P., Ambrosio J., Claro J.C.P., Lankarani H.M. Dynamic behavior of planar rigid multi-body systems including revolute joints with clearance. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, Vol. 221, Issue 2, 2007, p. 161-174.
  • Flores P., Lankarani H. M. Dynamic response of multibody systems with multiple clearance joints. Journal Compute Nonlinear Dynamics, Vol. 7, Issue 3, 2011, p. 636-647.
  • Machado M., Costa J., Seabra E., Flores P. The effect of the lubricated revolute joint parameters and hydrodynamic force models on the dynamic response of planar multibody systems. Nonlinear Dynamics, Vol. 69, Issues 1-2, 2012, p. 635-654.
  • Tian Q., Zhang Y. Q., Chen L. P., Yang J. Z. Simulation of planar flexible multibody systems with clearance and lubricated revolute joints. Nonlinear Dynamics, Vol. 60, 2010, p. 489-511.
  • Flores P. A parametric study on the dynamic response of planar multibody systems with multiple clearance joints. Nonlinear Dynamics, Vol. 61, Issue 61, 2010, p. 633-653.
  • Flores P., Ambrósio J., Claro J. C. P., Lankarani H. M., Koshy C. S. Lubricated revolute joints in rigid multibody systems. Nonlinear Dynamics, Vol. 56, Issue 3, 2009, p. 277-295.
  • Flores P., Lankarani H. M. Spatial rigid-multibody systems with lubricated spherical clearance joints: modeling and simulation. Nonlinear Dynamics, Vol. 60, Issue 1, 2010, p. 99-114.
  • Xu L. X., Han Y. C. A method for contact analysis of revolute joints with noncircular clearance in a planar multibody system. Proceedings of the Institution of Mechanical Engineers Part K: Journal Multi-body Dynamics, Vol. 230, Issue 4, 2016, p. 589-605.
  • Farahan S. B., Ghazavi M. R., Rahmanian S. Bifurcation in a planar four-bar mechanism with revolute clearance joint. Nonlinear Dynamics, Vol. 87, 2017, p. 955-973.
  • Wang Z., Tian Q., Hu H. Nonlinear dynamics and chaotic control of a flexible multibody system w ith uncertain joint clearance. Nonlinear Dynamics, Vol. 86, Issue 3, 2016, p. 1571-1597.
  • Nayfeh A. H., Balachandran B. Applied Nonlinear Dynamics. John Wiley and Sons, Inc., New York, 1995.
  • Dupac M., Beale D. G. Dynamic analysis of a flexible linkage mechanism with cracks and clearance. Mechanism and Machine Theory, Vol. 45, Issue 12, 2010, p. 1909-1923.
  • Dupac M., Marghitu D. B. Nonlinear dynamics of a flexible mechanism with impact. Journal of Sound and Vibration, Vol. 289, Issues 4-5, 2006, p. 952-966.
  • Wolf A., Swift J. B., Swinney H. L., Vastano J. A. Determining Lyapunov exponents from a time series Physica D, Vol. 16, 1985, p. 285-317.

Cited by

Nonlinear analysis of complex mechanisms with multi-clearances considering dry friction and lubricated joints
Zhimin Wang | Guoguang Jin | Dong Liang | Zhan Wei | Boyan Chang | Yang Zhou
(2023)
Nonlinear dynamics of coupling rub-impact of double translational joints with subsidence considering the flexibility of piston rod
Shungen Xiao | Shulin Liu | Hongzhou Wang | Yang Lin | Mengmeng Song | Hongli Zhang
(2020)
Positive Displacement Machines
Ibrahim A. Sultan | Truong H. Phung
(2019)
Mechanism Design of Equivalent Simulation of Cylinder Pressure Load in Reciprocating Compressor Fault System
Shun Gen Xiao | Meng Meng Song | Xiao Hu Sang | Li Xia Huang | Hong Long Chen | Tao Qiu
(2019)
Nonlinear dynamic response of reciprocating compressor system with rub-impact fault caused by subsidence
Shungen Xiao | Shulin Liu | Feng Jiang | Mengmeng Song | Shouguo Cheng
(2019)
Dynamic analysis of the moving mechanism of the reciprocating compressor with clearance joints
Qingfeng Xiao | Shulin Liu | Ang Nie | Shungen Xiao | Shouguo Cheng
(2018)

About this article

Received
16 June 2017
Accepted
20 August 2017
Published
15 November 2017
SUBJECTS
Fault diagnosis based on vibration signal analysis
Keywords
reciprocating compressor
dynamics
clearance
subsidence
chaos
Acknowledgements

This work is supported by the Natural Science Foundation of China (Grant No. 51575331 and Grant No. 51175316), Fujian Natural Science Foundation (Grant No. 2015J01643), and Talents Cultivation Program for Outstanding Young Scientists in Fujian Universities (Grant No. MIN Education (2015) 54). These supports are gratefully acknowledged.

Author Contributions

In this paper, Shungen Xiao is responsible for the writing of the entire paper, the realization of the core ideas and so on. Shulin Liu is the advocate of the core idea. Shouguo Cheng helps the first author to debug the programs of the dynamics equation. Xiaogan Xue assists the first author to establish a mathematical model of cylinder pressure. Mengmeng Song helps the first author to analyze the influence of the clearance size. Xin Sun helps the first author to debug the programs of the Lyapunov exponents.