The work performance of piston-cylinder liner system is affected by the lubrication condition and the secondary motion of the piston. Therefore, more and more attention has been paid to the secondary motion and lubrication of the piston. In this paper, the Jakobson-Floberg-Olsson (JFO) boundary condition is employed to describe the rupture and reformation of oil film. The average Reynolds equation of skirt lubrication is solved by the finite difference method (FDM). The secondary motion of piston-connecting rod system is modeled; the trajectory of the piston is calculated by the Runge-Kutta method. By considering the inertia of the connecting rod, the influence of the longitudinal and horizontal profiles of piston skirt, the offset of the piston pin, and the thermal deformation on the secondary motion and lubrication performance is investigated. The parabolic longitudinal profile, the smaller top radial reduction and ellipticities of the middle-convex piston, and the bigger bottom radial reduction and ellipticities can effectively reduce the secondary displacement and velocity, the skirt thrust, friction, and the friction power loss. The results show that the connecting rod inertia, piston skirt profile, and thermal deformation have important influence on secondary motion and lubrication performance of the piston.
National Natural Science Foundation of China51775428Open Project of State Key Laboratory of Digital Manufacturing Equipment and TechnologyDMETKF2017014Natural Science Foundation of Shaanxi Province2014JM2-5082Education Department of Shaanxi Province15JS0681. Introduction
The piston-cylinder liner system is a key subsystem of the internal combustion engines (ICEs); the performance of the piston-cylinder liner system influences the efficiency and performance of ICEs. The work performance of piston-cylinder liner system is affected by the lubrication condition and the secondary motion of the piston. Therefore, more and more attention has been paid to the piston skirt lubrication and the secondary motion of the piston by many scholars recently [1–6]. Meng et al. [7] introduced the oil film inertia force effects to the Reynolds lubrication equation; the effects of oil film inertia on film pressure, hydrodynamic friction force, and transverse displacement were investigated. A nonlinear model of the piston considering the reciprocating, lateral, and rotational degree of freedom was established by Tan and Ripin [8]. In their work, the piston secondary motion and the induced vibration behavior were studied; the proposed model was validated by the experimental data. Tan and Ripin [9] identified the intermittent transient contact of the piston skirt with the cylinder liner by measuring the tilting angle and secondary displacement. The impacting force was computed from the equation of motion and the frictional force ratio of the total friction force and torque was obtained. Murakami et al. [10] proposed a piston secondary motion analysis method coupling the structure analysis and the multibody dynamics analysis. The presented method was verified by comparing the numerical results with the measurement data of the piston skirt stress, secondary motion, and the vibrational acceleration produced on the cylinder. Besides the above works, more researches on the piston secondary motion have been carried out by taking some factors into consideration, such as surface roughness and lubrication condition. Meng et al. [11] investigated the influence of cylinder vibration on the lateral motion and tribological performances of the piston. The fluctuations of the cylinder vibration can be reduced by increasing mass, stiffness, and damping of the cylinder. Narayan [12] implemented a modeling of the piston secondary motion. The dynamic characteristics of the part of piston skirt and cylinder liner were analyzed, and the proposed model was validated by comparing the simulated results with the experimental data. Obert et al. [13] developed a reciprocating model test to investigate the influence of the temperature and oil supply on the tribological behavior of the piston ring-cylinder liner contacts. The real engine situation at fired top dead center was simulated by carefully adjusting the friction, scuffing, wear, oil supply, and temperature. Mazouzi et al. [14] presented a numerical secondary motion model to investigate the influence of the piston design parameters on the dynamic characteristics of a piston-cylinder contact. In their work, it was verified that the modification of clearance of the piston skirt-liner system and the axial location of piston pin can reduce the friction and piston slap. Fang et al. [15] modeled the piston skirt-liner lubrication and crank-rod-piston multibody system by taking the microgrooves on the piston skirt surface. Based on the presented model, the tribodynamic performance of the piston with skirt surface grooves was obtained and analyzed.
Regarding the influencing factors, the piston skirt profile has extreme significance for the secondary motion of the piston. Hence, the study of the effects of the piston skirt profile on the secondary motion and piston dynamics is essential. Littlefair et al. [16] analyzed the piston dynamics, thermoelastic distortion, and transient elastohydrodynamics of the lightweight piston with preferentially compliant short partial skirt. The greater entrainment wedge at bottom skirt and better conforming contact geometry at top skirt were achieved, and therefore, the load-carrying capacity was improved and the friction was reduced. McFadden and Turnbull [17] presented a model of the interface between a barrelled piston skirt and cylinder. In their work, the secondary motion and the pressure distribution of the oil film for changing the skirt profile were investigated. The research revealed that the wear of the cylinder arising from the rotation of the piston can be decreased by appropriate positioning of piston skirt. Gunelsu and Akalin [18] developed a dynamics model to investigate the secondary motion of piston for the different radius of curvature around the bulge of the barrel, the ellipse, and the aspect ratio. The secondary motion and frictional performance of the piston with the barrel and elliptical form of the skirt were analyzed. The results showed that the frictional performance can be improved via optimizing the piston skirt profiles. He et al. [19] established a numerical model of coupling lubrication and dynamic motion. The effects of the piston skirt parameters on dynamic performance were investigated, including the offset of piston pin, length of piston skirt, curvature parameter, and ellipticity. Furthermore, an asymmetrical skirt profile was proposed to reduce the thrust force on the liner. In order to investigate the influence of the design parameters of the piston skirt on the slap noise and the lubrication performance, the lubrication model of piston skirt-liner system and the dynamics model of multibody system were coupled by Zhao et al. [20].
Apart from the piston skirt profiles, the influence of the connecting rod inertia on lubrication performance of the piston skirt-liner system cannot be neglected. Therefore, more detailed dynamic analyses of the secondary motion by considering the inertia of the connecting rod were conducted [21–25]. Zhang et al. [21] developed a mathematical model of piston-connecting-rod-crankshaft system. Based on the coupling model of the secondary motion and the fluid dynamic lubrication, the influence of the variation in the system inertia on the piston secondary motion and the piston side force was investigated. The center-of-mass position was confirmed to be significant for the design of the connecting rod. Meng et al. [22, 23] conducted a numerical analysis of the piston dynamics, the oil film, and the friction loss of the piston skirt-liner system by taking the connecting rod inertia into consideration. The numerical results revealed that the influence of the connecting rod inertia on the side force, the oil film thickness, the frictional force, and the piston dynamics cannot be ignored, in particular at high speed. Furthermore, Meng et al. [24] presented a transient tribodynamic model by coupling the tribological performance and the dynamic characteristics of the piston skirt-liner system. By considering the angular acceleration of the crankshaft and the transient variable engine speeds, the secondary motion of the piston and engine performance were investigated. Zhu et al. [25] calculated the stress distribution of the connecting rod with consideration of the film lubrication or not. The minimum oil film thickness and maximum oil film pressure were calculated by the Reynolds equation. The results showed that the stress concentration can be reduced using smoother rounded fillet and longer fillet radius.
The thermal deformation is also an important influencing factor of the piston secondary motion and frictional characteristics. Therefore, the thermal deformation was increasingly introduced to relevant studies of the lubrication performance of the piston skirt. Pelosi and Ivantysynova [26] built a fully coupled simulation model for the piston/cylinder interface of axial piston machines and validated the model by comparing the prediction values with measurements of the fluid film boundary temperature distribution. The influence of the heat transfer and solids-induced thermal elastic deformation on the piston skirt-cylinder lubricating interface behavior was analyzed. Ning et al. [27] built a comprehensive lubrication model of the piston skirt-liner system. By considering the thermal and the elastic deformations of the piston skirt and the cylinder liner, the piston secondary motion was investigated.
The above-mentioned researches were carried out based on the Reynolds boundary. It is reasonable on the rupture boundary of oil film, but not on the reformation boundary. Jacoboson et al. [28, 29] presented the boundary conditions of conservation of mass on rupture and reformation boundaries of oil film, that is, JFO boundary condition. Since then, JFO boundary condition was considered as important factor in the investigation of the tribological properties. Zhao et al. [30] investigated the cavitation effects on the dynamic axial stiffness and damping coefficients of spiral-groove rotary seals for the Reynolds and JFO boundary conditions. By contrasting the theoretical value with the experimental data, the verification result was obtained that computation values of JFO boundary condition were in better agreement with the experimental values than the values using the Reynolds model. Kango et al. [31] conducted a mathematical model to investigate the influence of the surface microtexture on the performance parameters of the bearing by JFO boundary condition. The realistic results can be obtained based on JFO boundary condition in comparison to the Reynolds boundary condition. Liu et al. [32] investigated the influence of the cylinder liner dimples on the lubrication performance of the compression ring-cylinder liner system based on JFO boundary condition. For different dimple area density, radius, and depth, the minimum oil film thickness and friction force were compared and analyzed under the engine-like conditions. The results showed that the spherical dimples can improve the tribological performance.
In this paper, the average Reynolds equation of the piston skirt lubrication is solved by FDM, and the secondary motion trajectory of the piston is calculated by the Runge-Kutta method. The difference of the secondary motion of the piston under JFO and the Reynolds boundary conditions is significant. In order to approach the actual working condition, JFO boundary condition is employed, and the variable ϕ and cavitation index F are introduced to describe the full and rupture areas of oil film. By considering the connecting rod inertia, the effects of the longitudinal profiles of the piston skirt, the horizontal profiles, the offset of the piston pin, and the thermal deformation on the secondary motion and lubrication performance are investigated.
2. Secondary Motion Model
The piston makes reciprocating movement along cylinder liner; the crankshaft makes nonuniform rotation caused by piston. Figure 1 shows the geometrical model of piston-connecting rod mechanism.
Geometrical model of piston-connecting rod mechanism.
In Figure 1, Fg is the thrust of the gas on the top surface of the piston, F~IC and F~IP are the inertia forces of the piston and piston pin caused by reciprocating motion of piston, respectively, FIC, FIP, and MIP are the inertia forces of piston, piston pin, and the inertia moment of piston caused by the secondary motion of piston, respectively, Ft and Mt are the total lateral thrust force and lateral thrust moment of piston skirt, respectively, Ff and Mf are the friction and friction torque of piston skirt, respectively, FBX and FBY are the forces of connecting rod acting on piston in the x and y directions, respectively, g is the gravitational acceleration, κ is the angle between gravity direction and y-axis, dCM and yCM are the distances from center of mass of piston to center line of piston and to top skirt, respectively, dP and yP are the distances from center of piston pin hole to center line of piston and to top skirt, respectively, L is the length of the piston skirt, and α is the angle between the connecting rod and the center line of piston.
According to Figure 1, the force equation and moment equation of piston are established(1)Ft+FIC+FIP+FBX+mpgsinκ+mpingsinκ=0,where mp and mpin are the mass of piston and piston pin, respectively.(2)Fg+Ff+F~IC+F~IP+FBY+mpgcosκ+mpingcosκ=0Mt+Mf+MIC-Fgdp-mpgcosκdCM+dp+mpgsinκyP-yCM+FICyP-yCM-F~ICdCM+dp=0.
The angle α between the center line of piston and the connecting rod is(3)α=tan-1dc+ecsinθ×lcr2-dc+ecsinθ2-0.5,where ec is the length of the crankshaft, dc is the eccentricity of the centers of crankshaft and piston pin, and θ is the crank angle.
The inertia forces of piston and piston pin caused by reciprocating motion of piston are(4)F~IC=-mpy¨BF~IP=-mpiny¨B,where y¨B is the acceleration of the reciprocating motion of piston.
As shown in Figure 2, the piston swings around the center of pin hole between the major thrust side (the side acted upon by the expansion thrust, φ=0) and the minor thrust side (the side acted upon by the compression thrust, φ=π). In order to describe the motion of the piston, the displacement et of top skirt and displacement eb of bottom skirt are introduced to express the displacement, velocity, and tilting angle of piston secondary motion. In Figure 2, φ is the circumferential coordinate of piston, φl is the loading angle of piston skirt, γ is the tilting angle of piston around center of pin, and R is the radius of piston.
Schematic diagram of the piston-cylinder system.
The inertia forces of piston and piston pin and the inertia moment of piston caused by secondary motion are, respectively, (5)FIC=-mpe¨t+yCMLe¨b-e¨t,where e¨b, e¨t are the acceleration values of the top skirt and bottom skirt.(6)MIC=-Ipe¨t-e¨bL,where Ip is the inertia moment of piston around center of pin.(7)FIP=-mpine¨t+yPLe¨b-e¨t.
By considering the inertia of connecting rod, the equilibrium equation of connecting rod is obtained.(8)FAX-FBX+mcrgsinκ-mcrx¨cr=0FAY-FBY+mcrgcosκ-mcry¨cr=0,where FAX and FAY are the component forces caused by the crank acting on the connecting rod in the x and y directions, respectively, mcr is the mass of the connecting rod, and x¨cr and y¨cr are the components of acceleration of center of mass of the connecting rod in the x and y directions, respectively.
The moment equilibrium equation of the connecting rod is(9)-FBX1-jlcrcosα-FAXjlcrcosα+FBY1-jlcrsinα+FAYjlcrsinα-α¨Icr=0,where Icr is the moment of the connecting rod and j is the ratio of distances from center of mass and big side to length of connecting rod, as follows:(10)j=lAClAB.
According to the equilibrium equation of force and moment of the piston and connecting rod, the equations are written as follows:(11)A11A12A21A22e¨te¨b=Fs+Fss+FtMs+Mf+Mt,where Mt is the lateral thrust moment of piston skirt, A11, A12, A21, and A22 are relevant coefficients, Ft is the lateral thrust force of piston skirt, Fs and Ms are the physical quantities of the parameters of the piston structure, and Fss is the physical quantity of the parameters of the connecting rod; the expressions are as follows: (12)A11=mpin1-yPL+mp1-yCML+j2mcr1-yPLA12=mpinyPL+mpyCML+j2mcryPLA21=mp1-yCMLyP-yCM+IpLA22=mpyP-yCMyCML-IpLFs=-Fg-Ff-F~IC-F~IP-mpgcosκ-mpingcosκtanα+mpgsinκ+mpingsinκMs=-Fgdp-mpgcosκdCM+dp+mpgsinκyP-yCM-F~ICdCM+dPFss=jmcry¨cr-jmcrgcosκtanα+jmcrgsinκ+j1-jmcrecω2sinθ-α¨Icrlcrcosα.
3. Profile Design of Piston Skirt3.1. Horizontal Profile Design
The top surface of piston is acted upon by a larger gas pressure when the ICE operates; the piston expands outward along the pin direction because of the span between the right and left pin bosses. Therefore, the horizontal profile of piston skirt is designed as ellipse. In Figure 3, d1, d2 are the diameters of the major and minor axis of ellipse, φ is the circumferential angle, and Δ is the radial reduction of ellipse.
Horizontal profile of piston skirt.
Because of the special structure and working condition of piston, the external expansion deformation of bottom skirt is more than top skirt. Therefore, the top ellipse profile of piston is designed to be slightly smaller than the bottom one.
3.2. Longitudinal Profile Design
In order to improve the lubrication condition of piston skirt, the longitudinal profile of piston skirt is generally designed as dolioform to reduce the friction force of the motion and strengthen the guidance effect. Meanwhile, it is beneficial to form the hydrodynamic oil film. The longitudinal profile is shown in Figure 4; δ is the radial reduction of longitudinal profile of the piston skirt, δt is the radial reduction of top skirt, δb is the radial reduction of bottom skirt, and L0 is the distance from the middle-convex point to top skirt. The dolioform profiles of skirt are normally designed as parabola, ellipse, and circle curve.
Longitudinal profile of piston skirt.
4. Mixed-Lubrication Model of Piston Skirt4.1. Hydrodynamic Lubrication Model
By considering the effect of lubricated surface roughness on flow of oil film, the oil film pressure between piston skirt and cylinder liner is solved using two-dimensional average Reynolds equation based on the theory proposed by Patir and Cheng [33, 34]. (13)∂∂xϕxh3∂p∂x+∂∂yϕyh3∂p∂y=6μU∂hT∂x+σ∂ϕs∂x+12μ∂hT∂t,where hT is the average value of actual oil film thickness, p is the average hydrodynamic pressure, h is the theoretical oil film thickness, μ is the dynamic viscosity of lubricant, σ is the roughness of friction surface, ϕx and ϕy are the pressure flow factors, respectively, ϕs is the shear flow factor, and U is the relative velocity between the piston skirt and cylinder liner.
The pressure flow factors ϕx, ϕy [33] and shear flow factor ϕs [34] can be calculated as follows:(14)ϕx=ϕy=1-0.9e-0.56H′ϕs=1.899H′0.98·exp-0.92H′+0.05H′2H″≤51.126·exp-0.25H′H′>5,where H′=h/σ is the ratio of oil film thickness, σ is the combined surface roughness, and Ωp is the surface waviness of the piston.
Define touch factor ϕc=∂h-T/∂h [35]; while 0≤H′<3, ϕc can be calculated from(15)ϕc=exp-0.6912+0.782H′-0.304H′2+0.0401H′3,while H′≥3, ϕc=1.
The Reynolds equation can be written as (16)∂∂xϕxh3∂p∂x+∂∂yϕyh3∂p∂y=6μUϕc∂h∂x+6μUσ∂ϕs∂x+12μϕc∂h∂t.h can be calculated from(17)hφ,y,t=c+ettcosφ+ebt-ettyLcosφ+fφ,y.
By considering the elastic deformation of the piston, the oil film thickness can be rewritten as (18)hφ,y,t=c+ettcosφ+ebt-ettyLcosφ+fφ,y+dx,y,where d(x,y) is the elastic deformation of the piston caused by the oil film pressure; d(x,y) is obtained by the following formula [36]:(19)dx,y=1πE′∫0L∫0bpx,ydxdyx-x12-y-y12,where E′ is the comprehensive elastic modulus of contact surface, which is related to elastic moduli E1, E2 and Poisson ratios of piston and cylinder liner ν1, ν2.(20)E′=21-ν12E1+1-ν22E2.
By considering the symmetry of the structure of pistons skirt, the oil film force is calculated in a half of area (φ=0~π). On major thrust side (φ=0) and minor thrust side (φ=π), the boundary conditions can be expressed as(21)∂p∂φφ=0=∂p∂φφ=π=0.On the top skirt (y=0) and the bottom skirt (y=L), the boundary conditions can be written as(22)pφ,0=pφ,L=0.
Because the squeeze effect of the oil film decreases gradually with the increase of the loading angle, it has smaller influence on the lubrication performance and the secondary motion. Taking certain value of loading angle, the oil film pressure approximates to zero outside of the loading angle area; the corresponding boundary condition is(23)p=0φ1≤φ≤φ2.According to JFO boundary condition, the rupture boundary condition of oil film is(24)∂p-∂n=0;the reformation boundary condition of oil film is(25)h212μ∂p∂n=Un21-ρρ0,where n is the flowing direction of the oil film.
In order to couple the full oil film zone and the cavitation zone and determine the rupture and reformation boundary of the oil film, the dimensionless independent variable ϕ and cavitation index F are introduced:(26)Fϕ=p-ϕ≥0 on full oil film regionρρ0-1=1-Fϕon rupture region of oil film,where ρ is the fluid density in the full film region and ρ0 is the mix-density of the gas-fluid in the cavitation region.
The cavitation index F is(27)F=1ϕ≥00ϕ<0.
The Reynolds equation is converted into dimensionless equation and solved by the finite difference method.
As shown in Figure 5, the area of the oil film is divided into M×N elements. In the φ direction, the number of the node is M and the grid step is Δφ; in the y direction, the number of the node is N and the grid step is Δy.
Grid of oil film region of the piston skirt.
The grid independence is verified for three sets of different numbers of the node in the φ and y directions. Table 1 shows the comparison of the oil film pressures for different numbers of the grid (the pressures are at three points. Point 1 for φ=0°, y=0.0414; Point 2 for φ=28.8°, y=0.0207; Point 3 for φ=165°, y=0.0621; crankshaft angle is 90°). It can be seen that the changes of the oil film pressures are small with the increase of the grid number; the deviations among the oil film pressures are about 1% that can be neglected. The grid is independent of the calculations, but the calculation cost will increase when the grid with a bigger number is chosen. Therefore, the node numbers M=36 and N=20 are chosen as calculation grid number.
The numbers of the node in φ and y directions.
M(φ)
N(y)
Oil film pressure at point 1
Oil film pressure at point 2
Oil film pressure at point 3
36
20
554091 MPa
155160 MPa
103411 MPa
50
30
559202 MPa
156861 MPa
103376 MPa
80
60
561570 MPa
157100 MPa
103729 MPa
The tangential pressure of oil film on the skirt is calculated by (28).(28)τh=-μUhΦf+Φfs+Φfph2∂p∂y,where Φf, Φfs, and Φfp are shear stress factors related to lubrication surface roughness [37].
4.2. Solid-to-Solid Contact Model
Because the instantaneous gas pressure of piston is variable, the lubricant film thickness between the piston skirt and cylinder liner may get thinner during strenuous movement. Therefore, the solid-solid contact may be caused by the collision of the surface asperity by considering the surface roughness. In this situation, the total lateral thrust force consists of the lateral thrust force of oil film and lateral thrust force of solid-solid contact and the total friction force consists of the oil film shear stress and asperity shear stress.
The ratio of oil film thickness H′ is normally used to evaluate the lubrication condition; when H′≥4, the piston skirt-cylinder liner system is under hydrodynamic lubrication condition; when H′<4, the piston skirt-cylinder liner system is under mixed-lubrication condition, and asperity contact [38] is inescapable.
Based on the rough surface contact theory presented by Greenwood and Tripp [39], the solid-solid contact pressure pc of the piston skirt and cylinder liner is calculated by (29)pc=16215πηβσ2E′σβF5/2H′,where η is the asperity density of the rough surface, β is the radius of the asperity curvature, and ηβσ = 0.03~0.05,σ/β = 10-4~10−2.(30)αc=π2ηβσ2F2H′,where αc is the asperity contact area.
The shear force τc of solid-solid contact of piston skirt and cylinder liner is(31)τc=μf·pc,where μf is the frictional coefficient of solid-solid contact.
4.3. Calculation of Force and Moment
The lateral thrust force Fh caused by hydrodynamic lubrication is(32)Fh=R∫0L∫02πpcosφdφdy.The lateral thrust moment Mh caused by hydrodynamic lubrication is(33)Mh=R∫0L∫02πpcosφyP-ydφdy.The friction force Ffh caused by hydrodynamic lubrication is(34)Ffh=-signUR∫0L∫02πτhdφdy,where sign(U) is sign function.
The friction torque Mfh caused by hydrodynamic lubrication is(35)Mfh=signUR∫0L∫02πτhRcosφdφdy.The lateral thrust Fc caused by solid-solid contact is(36)Fc=R∫0L∫02πpccosφdφdy.The lateral thrust moment Mc caused by solid-solid contact is(37)Mc=R∫0L∫02πpccosφyP-ydφdy.The friction force Ffc caused by solid-solid contact is(38)Ffc=-signUR∫0L∫02πτcdφdy.The friction moment Mfc caused by solid-solid contact is(39)Mfc=signUR∫0L∫02πτcRcosφdφdy.The total lateral thrust force Ft of piston skirt is(40)Ft=Fh+Fc.The total friction force Ff of piston skirt is(41)Ff=Ffh+Ffc.The total lateral thrust moment Mt of piston skirt is(42)Mt=Mh+Mc.The total friction torque Mf of piston skirt is(43)Mf=Mfh+Mfc.
5. Thermal Analysis Model of Piston
The heat of the piston side is mainly transferred to the oil film and gas film in the cylinder clearance through the piston junk, piston ring, and piston skirt surface. Then the heat is transferred to coolant outside of cylinder liner through cylinder liner. The coolant is considered as heat conduction terminal. Figure 6 shows the structure of the piston ring groove.
Structure of piston ring groove area.
The heat transfer of the piston and cylinder liner is simplified to steady state heat conduction in multilayer flat wall; the heat transfer of cylinder liner and coolant is simplified to convective heat conduction state.(44)Φ=tw1-tw4δ1/λ1A+δ2/λ2A+δ3/λ3A+δ4/λ4A,where Φ is the heat flow of the heat conduction, tw1 and tw4 are the starting and terminal temperature, respectively, δ1 is the thickness of gas or lubricant, δ2 is the thickness of piston ring, δ3 is the thickness of cylinder wall, λ1 is the heat conduction coefficient of gas λg or heat conduction coefficient of lubricant λo, λg=0.023 W/(m2·K), λo=0.4 W/(m2·K), λ2 is the heat conduction coefficient of piston ring, λ2=0.2 W/(m2·K), λ3 is the heat conduction coefficient of cylinder liner, λ4 is the convective heat conduction coefficient of cylinder liner and coolant, and A is the heat conduction area; the heat conduction coefficient a of different positions of the piston side is deduced from deformation formula of (44).(45)a=1δ1/λ1+δ2/λ2+δ3/λ3+1/λ4.
The piston side consists of the piston junk, piston ring groove, piston ring land, piston pin area, and piston skirt. Because the structures are different in different areas, the calculation of heat conduction coefficients is different.
The heat conduction of the coolant and cylinder liner is simplified to convection heat transfer from tube banks in crossflow. According to strong turbulence formula and Nusselt number definition, the convective heat conduction coefficients of cylinder liner and coolant are obtained(46)λ4=0.023Ref0.8PrfmkL,where λ4 are the convective heat conduction coefficients of cylinder liner and coolant, Ref is within 104~1.2 × 105, Ref=11000, Prf is within 0.7~160, Prf = 100, when the fluid is heated by pipe, m=0.4, k is the heat conduction coefficient of static fluid, k=0.0036 W/(m·K), and L is the geometric characteristic length, L=0.005 m.
When the inner of cylinder liner contacts with the oil-gas mixture, the convective heat conduction coefficient is calculated by the following equation:(47)ac=1δc/λc,where ac is the convective heat conduction coefficient of the inner of the cylinder liner, δc is the thickness of the oil-gas mixture of the inner of the cylinder liner, and λc is the heat conduction coefficient of oil-gas mixture.
6. Model Solutions
Equation (11) is solved by the Runge-Kutta method. Taking the four-stroke ICE as an example, the piston motion is periodic motion. The simulation of the secondary motion is conducted in an engine cycle (i.e., the crankshaft angle is from 0° to 720°), the time step is Δt=2×10-6 s, the initial values et,eb,e˙t,e˙b,e¨t, and e¨b at the moment ti are provided, and the oil film pressure and solid-solid contact pressure of the hydrodynamic lubrication are calculated by coupling the motion parameters of the crankshaft connecting rod system. The shear stress of oil film and solid-solid contact shear stress are solved by (28) and (31), and the force and moment are calculated using integral formula. The values et, eb, e˙t, e˙b, e¨t, and e¨b at the moment ti+1 are calculated by the Runge-Kutta method. Repeat the above procedure until meeting convergence criteria.(48)et-ett+4πω<10-5eb-ebt+4πω<10-5e˙t-e˙tt+4πω<10-5e˙b-e˙bt+4πω<10-5,where ω is the rotational angular velocity of the crankshaft. The convergence can be obtained after calculating 2-3 cycles.
7. Numerical Examples
The simulation of the piston secondary motion is carried out in a certain-type four-stroke four-cylinder turbocharged diesel engine. The parameters used in the simulation are listed in Tables 2 and 3. Table 2 shows the structure parameters of engine. Table 3 shows the material characteristic parameters of piston and cylinder liner. The variation curve of the gas pressure for different crankshaft angles is shown in Figure 7, when the rotation speed is 1680 rpm.
Structure parameters of engine.
Symbol
Definition
Value
R
Piston radius
55 mm
L
Piston skirt length
82.8 mm
r
Crankshaft radius
72 mm
lAB
Connecting rod length
223 mm
lAC
Connecting rod center-of-mass location
72 mm
yCM
Piston center-of-mass location
10.8 mm
yP
Piston pin axial location
37.3 mm
mP
Piston mass
1.39 kg
Ip
Piston inertia moment
0.003 kg·m2
mcr
Connecting rod mass
1.7 kg
Ib
Connecting rod inertia moment
0.017 kg·m2
mc
Crankshaft mass
23.5 kg
Ic
Crankshaft moment of inertia
2 kg·m2
Ω
Crankshaft rotation speed
1680 rpm
c
Theoretical piston skirt clearance
50 μm
σs
Piston surface roughness
0.20 μm
Ωp
External waviness root-mean-square value of piston
20 μm
σb
Inner hole surface roughness
0.50 μm
σ
Combined roughness
0.54 μm
E
Modulus of elasticity (Al)
79 GPa
μ
Viscosity (100°, SAE 10W40)
0.0056 Pa·s
μf
Solid-solid friction coefficient
0.1
dc
Offset of crankshaft
0 mm
dp
Offset of pin
0 mm
dCM
Offset of center of mass
0 mm
κ
Cylinder axial tilting angle
0°
Characteristic parameters of piston material.
Physical quantity
Piston (eutectic Al Si alloy)
Cylinder liner (boron cast iron)
Thermal conductivity (λ/W·m−1·K−1)
145
42
Density (ρ/kg/m−3)
2700
7370
Specific heat capacity (cp/J·kg−1·K−1)
880
470
Elasticity modulus (E/GPa)
71
159
Poisson ratio (ν)
0.34
0.3
Coefficient of linear expansion (α/K−1)
2 × 10−6
12 × 10−6
Gas pressure inside of cylinder.
Figure 8 shows the comparison of the secondary motions of the piston with and without considering elastic deformation. When the elastic deformation of the piston is considered, the maximum secondary displacement and velocity increase, respectively, about 0.15% and 0.37%, the maximum of the minimum oil film thickness increases about 0.20%, and the friction force and friction power loss increase, respectively, about 0.015% and 0.016%. From Figure 8, it can be seen that the effect of elastic deformations is small within the selected data set, and therefore the piston was assumed to be rigid. So the elastic deformation of the piston is not considered as the influence factor in the current work.
Comparison of secondary motions of the piston with and without elastic deformation.
Secondary displacement
Secondary velocity
Minimum oil film thickness
Friction force
Friction power loss
7.1. Effect of the Boundary Condition on Secondary Motion
Figure 9 shows the comparison of secondary motions of the piston under the Reynolds and JFO boundary conditions. By considering the rupture and reformation boundaries of oil film, the secondary displacement, velocity, and the tilting angle under JFO boundary condition decrease obviously. Compared with the Reynolds boundary condition, because of the change of the oil film pressure, the maximum of the minimum oil film thickness on major thrust side increases at 135° and 495° and decreases at 270° and 630°; the corresponding friction force and friction power loss increase to the peak at 270° and 630°. The maximum friction force and friction power loss under JFO boundary condition increase, respectively, about 48.4% and 44.7%. The results show that the reformation of oil film has effect on piston secondary motion, so it is necessary to investigate the lubrication performance of the piston using JFO boundary condition.
Comparison of secondary motions of the piston under the Reynolds and JFO boundary conditions.
Displacement of the center of the piston pin
Friction force
Minimum oil film thickness
Friction power loss
7.2. Effect of Connecting Rod Inertia on Secondary Motion
Figure 10 shows the curve of secondary motion of the piston by considering the inertia of the connecting rod. The connecting rod inertia leads to the obvious change of the oil film pressure during the compression stroke and exhaust stroke, and the change of the pressure has an important influence on the secondary motion. In a cycle, the maximum secondary displacement and velocity and the maximums of the minimum oil film thicknesses on major and minor thrust sides increase slightly, the maximum thrust and friction forces increase, respectively, about 15.5% and 7.9%, the piston tilting angle has obvious increase at 360° and 720°, and the maximum friction power loss increases 6%. By considering the effect of the connecting rod inertia, the secondary motion is closer to the actual working condition, so the inertia of connecting rod is taken into consideration in the following numerical examples.
Effect of connecting rod inertia on secondary motion.
Secondary displacement
Secondary velocity
Displacement of the center of the piston pin
Piston tilting angle
Lateral thrust force
Friction force
Minimum oil film thickness
Friction power loss
7.3. Effect of Longitudinal Profile on Lubrication Performance
Figure 11 shows the profiles of three different longitudinal profile skirts. Profile 1 is the profile of parabola (the deviations of top and bottom skirts are, resp., 80 μm and 130 μm), profile 2 is the profile of ellipse (the deviations of top and bottom skirts are, resp., 94.95 μm and 135.55 μm), and profile 3 is the circular curve (the deviations of top and bottom skirts are, resp., 88.89 μm and 101.77 μm). Figure 12 shows the effects of different profiles on lubrication performance of the piston-cylinder liner system. The secondary displacements and velocities of top and bottom skirts and the maximum of the minimum oil film thickness of profile 1 decrease obviously contrasting with profiles 2 and 3, and the maximum tilting angle of profile 1 decreases. Due to the decrease of the minimum oil film thickness, the friction force and friction power loss increase slightly. By considering the secondary motion of three profiles comprehensively, profile 1 is chosen as longitudinal profile.
Profiles of different longitudinal profile skirts.
Effect of longitudinal profile on lubrication performance.
Secondary displacement of top skirt
Secondary displacement of bottom skirt
Secondary velocity of top skirt
Secondary velocity of bottom skirt
Piston tilting angle
Friction force
Minimum oil film thickness
Friction power loss
7.4. Effect of Radial Reduction on Lubrication Performance
Figure 13 shows the skirt profiles of different radial reductions on the top and bottom skirts, the deviations of top and bottom skirts of profile 1 are, respectively, 80 μm and 130 μm, the deviations of profile 4 are both 100 μm, and the deviations of profile 5 are, respectively, 130 μm and 80 μm. Figure 14 shows the effects of different radial reductions of top and bottom skirts on lubrication performance of the piston-cylinder liner system. The secondary displacement and velocity of top skirt of profile 1 decrease, while the secondary displacement and velocity of bottom skirt increase. The change of secondary displacements and velocities is beneficial to decrease the friction force and friction power loss. The maximum piston tilting angle of profile 1 decreases significantly compared to profiles 4 and 5, and the friction power loss decreases, respectively, about 3.7% and 5%. In Figure 14, the radial reduction of profile 1 is a more reasonable design parameter.
Profiles of different radial reduction skirts.
Effect of radial reduction on lubrication performance.
Secondary displacement of top skirt
Secondary displacement of bottom skirt
Secondary velocity of top skirt
Secondary velocity of bottom skirt
Piston tilting angle
Friction force
Minimum oil film thickness
Friction power loss
7.5. Effect of Middle-Convex Point on Lubrication Performance
Figure 15 shows the skirt profiles of different middle-convex points, the middle-convex point of profile 1 is in the middle of piston, that is, L0=40 mm, the middle-convex point of profile 6 shifts to top skirt, that is, L0=30 mm, and the middle-convex point of profile 7 shifts to bottom skirt, that is, L0=50 mm. Figure 16 shows the effects of different middle-convex points on lubrication performance of the piston-cylinder liner system. The secondary displacement and velocity of top skirt of the profile 1 increase contrasting with profile 6, the secondary displacement and velocity of bottom skirt of the profile 1 increase contrasting with profile 7. The maximum piston tilting angle of profile 1 increases contrasting with profile 6, and the maximum friction force and friction power loss increase slightly. Compared with profile 1, the major minimum oil film thickness of profile 6 decreases obviously for 100°~250° and 460°~610°, and the minor minimum oil film thickness of profile 7 decreases obviously for 225°~315° and 585°~675°. By considering the secondary displacements, velocities, and minimum oil film thicknesses, profile 1 is more reasonable to obtain better lubrication performance.
Profiles of different middle-convex point skirts.
Effect of middle-convex point on lubrication performance.
Secondary displacement of top skirt
Secondary displacement of bottom skirt
Secondary velocity of top skirt
Secondary velocity of bottom skirt
Piston tilting angle
Friction force
Minimum oil film thickness
Friction power loss
7.6. Effect of Horizontal Profile on Lubrication Performance
Figure 17 shows the horizontal profiles of different ellipticities. The ellipticities of top and bottom skirts of profile 1 are, respectively, 90 μm and 150 μm. The ellipticities of profile 8 are both 120 μm. The ellipticities of profile 9 are, respectively, 150 μm and 90 μm. Figure 18 shows the effects of horizontal profiles of different ellipticities on lubrication performance of the piston-cylinder liner system. Contrasting with profiles 8 and 9, the secondary displacement and velocity of the top skirt of profile 1 are smaller, and the secondary displacement and velocity of bottom skirt are bigger. The friction force and friction power loss of profile 1 decrease slightly, and the piston tilting angle of profile 1 increases compared with profiles 8 and 9. The results indicate that the lubrication performance is better when the horizontal profile is designed as the smaller ellipticity of top skirt and bigger ellipticity of bottom skirt.
Horizontal profiles of different ellipticities.
Effect of horizontal profile on lubrication performance.
Secondary displacement of top skirt
Secondary displacement of bottom skirt
Secondary velocity of top skirt
Secondary velocity of bottom skirt
Piston tilting angle
Friction force
Minimum oil film thickness
Friction power loss
7.7. Effect of the Offset of the Piston Pin Hole on Skirt Lubrication Performance
Figure 19 shows the effect of the offset of the piston pin hole on lubrication performance of the piston-cylinder liner system. By contrasting with the nonoffset of piston pin hole, the secondary displacements and velocities of top and bottom piston skirts increase when the piston pin hole shifts to the major thrust side. The trend of curve is contrary when the piston pin hole shifts to the minor thrust side. The reason of the difference is that the position of the piston pin hole determines the swing direction of the piston. Contrasting with the nonoffset of piston pin hole, the maximum displacement of the piston pin center decreases about 8% while the pin hole shifts to the major thrust side, the maximum friction force decreases about 19.6%, and the friction power loss decreases about 21.6%.
Effect of piston pin offset on skirt lubrication performance.
Secondary displacement of top skirt
Secondary displacement of bottom skirt
Secondary velocity of top skirt
Secondary velocity of bottom skirt
Displacement of piston pin center
Friction force
Minimum oil film thickness
Friction power loss
7.8. Effect of Thermal Deformation on Lubrication Performance
Figure 20 shows the temperature distribution nephogram of the piston, and the maximum temperature is 520.11 K and the minimum temperature is 364.43 K. Figure 21 shows the profiles of piston skirt and cylinder liner before and after thermal deformation, when the piston and cylinder liner are in the expansion stroke. Figure 22 shows the comparison of the secondary motions of piston, piston tilting angles, minimum oil film thicknesses, friction forces, and friction power losses before and after thermal deformation. As Figure 22 shows, the maximum secondary displacement of the piston pin center decreases, and the piston tilting angle increases after thermal deformation. The minimum oil film thickness decreases greatly due to the distortion of the piston skirt profile caused by the thermal deformation, and the friction force and friction power loss increase correspondingly. From Figure 22, the effect of the thermal deformation on lubrication performance of the piston-cylinder liner system is obvious.
Distribution nephogram of the piston temperature.
Profiles of piston skirt and cylinder liner before and after thermal deformation.
Profiles of piston skirt
Profiles of cylinder liner
Effect of thermal deformation on lubrication performance.
Secondary displacement of piston pin center
Piston tilting angle
Friction force
Minimum oil film thickness
Friction power loss
8. Conclusion
In this paper, JFO boundary condition is employed to describe the rupture and reformation of oil film of the piston-connecting rod system. By considering the inertia of the connecting rod, the influence of the longitudinal and horizontal profiles of piston skirt, the offset of the piston pin, and the thermal deformation on the secondary motion and lubrication performance is investigated. The results are as follows.
(1) Contrasting with the Reynolds boundary condition, the secondary motion of the piston under JFO boundary condition has a significant difference. The tribological performance of the piston-cylinder liner system under JFO boundary condition is closer to the realistic condition. By considering the connecting rod inertia, based on JFO boundary condition, the piston dynamic behaviors and mixed-lubrication performance of piston-cylinder liner system are analyzed. The results show that the effects of the connecting rod inertia, piston skirt profile, and thermal deformation on piston secondary motion and lubrication performance of piston skirt are significant, and the effect of elastic deformation is small within the selected data set.
(2) When the piston skirt profile is designed as parabola, the smaller the radial reduction of top skirt and the bigger the radial reduction of bottom skirt, the smaller the ellipticity of top skirt and the bigger the ellipticity of bottom skirt; the lubrication performance of piston skirt is superior.
(3) When the piston pin slightly shifts to the major thrust side, the secondary displacement and velocity of piston increase, the friction force and friction power loss decrease, and the friction characteristics of piston skirt are better.
(4) By considering thermal deformation of the piston skirt, the minimum oil film thickness obviously decreases; furthermore, the piston tilting angle, friction force, and friction power loss significantly increase. The influence of the thermal deformation of the piston skirt on the lubrication performance is a nonnegligible factor, so further research on the thermal deformation needs to be conducted in the following work.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The research was supported by National Natural Science Foundation of China (Grant no. 51775428), Open Project of State Key Laboratory of Digital Manufacturing Equipment and Technology (Grant no. DMETKF2017014), Natural Science Foundation of Shaanxi Province of China (Grant no. 2014JM2-5082), and Scientific Research Program of Shaanxi Provincial Education Department of China (Grant no. 15JS068).
MansouriS. H.WongV. W.Effects of piston design parameters on piston secondary motion and skirt-liner friction200521964354492-s2.0-3014443168510.1243/135065005X34026J03604LivanosG. A.KyrtatosN. P.Friction model of a marine diesel engine piston assembly20074010–121441145310.1016/j.triboint.2007.01.0202-s2.0-34548149879KirnerC.HalbhuberJ.UhligB.OlivaA.GrafS.WachtmeisterG.Experimental and simulative research advances in the piston assembly of an internal combustion engine2016991591682-s2.0-8496217761510.1016/j.triboint.2016.03.005WangY.LiY.YangJ.SunC.LiuK.FengR.Effects of anodic oxidation for combustion chamber on heat transfer of the piston in an aero-engine201730146652-s2.0-8498842898110.1080/08916152.2015.1135203ShahmohamadiH.MohammadpourM.RahmaniR.RahnejatH.GarnerC. P.Howell-SmithS.On the boundary conditions in multi-phase flow through the piston ring-cylinder liner conjunction2015901641742-s2.0-8492908727010.1016/j.triboint.2015.04.025ZhangJ.MengY.Direct observation of cavitation phenomenon and hydrodynamic lubrication analysis of textured surfaces20124621471582-s2.0-8486076535410.1007/s11249-012-9935-6MengF. M.ZhangY. Y.HuY. Z.WangH.Thermo-elasto-hydrodynamic lubrication analysis of piston skirt considering oil film inertia effect2007407108910992-s2.0-3394711578210.1016/j.triboint.2006.10.006TanY.-C.RipinZ. M.Analysis of piston secondary motion201333220516251762-s2.0-8487995688010.1016/j.jsv.2013.04.042TanY.-C.RipinZ. M.Technique to determine instantaneous piston skirt friction during piston slap2014741451532-s2.0-8489696516010.1016/j.triboint.2014.02.014MurakamiH.NakanishiN.OnoN.KawanoT.New three-dimensional piston secondary motion analysis method coupling structure analysis and multi body dynamics analysis20125142502-s2.0-8485932913710.4271/2011-32-0559MengF. M.WangX. F.LiT. T.ChenY. P.Influence of cylinder liner vibration on lateral motion and tribological behaviors for piston in internal combustion engine201522921511672-s2.0-8492128323310.1177/1350650114545140NarayanS.Effects of Various Parameters on Piston Secondary Motion20152015-March2-s2.0-8493837466510.4271/2015-01-0079ObertP.MüllerT.FüßerH.-J.BartelD.The influence of oil supply and cylinder liner temperature on friction, wear and scuffing behavior of piston ring cylinder liner contacts—a new model test20169430631410.1016/j.triboint.2015.08.0262-s2.0-84942275038MazouziR.KellaciA.KarasA.Effects of piston design parameters on skirt-liner dynamic behavior20166822502582-s2.0-8495952087010.1108/ILT-07-2015-0103FangC.MengX.XieY.A piston tribodynamic model with deterministic consideration of skirt surface grooves20171102322512-s2.0-8501366583410.1016/j.triboint.2017.02.026LittlefairB.De La CruzM.TheodossiadesS.MillsR.Howell-SmithS.RahnejatH.Dwyer-JoyceR. S.Transient tribo-dynamics of thermo-elastic compliant high-performance piston skirts201453151702-s2.0-8489291890210.1007/s11249-013-0243-6McFaddenP. D.TurnbullS. R.A study of the secondary piston motion arising from changes in the piston skirt profile using a simplified piston skirt model2013227138472-s2.0-8487584458210.1177/0954406212444509GunelsuO.AkalinO.The effects of piston skirt profiles on secondary motion and friction201413662-s2.0-8489326957010.1115/1.4026486062503HeZ.XieW.ZhangG.HongZ.ZhangJ.Piston dynamic characteristics analyses based on FEM method Part I: Effected by piston skirt parameters20147568852-s2.0-8490315873810.1016/j.advengsoft.2014.05.003ZhaoB.DaiX.-D.ZhangZ.-N.XieY.-B.A new numerical method for piston dynamics and lubrication analysis2016943954082-s2.0-8494515774210.1016/j.triboint.2015.09.037ZhangZ.XieY.ZhangX.MengX.Analysis of piston secondary motion considering the variation in the system inertia200922345495632-s2.0-6654908682710.1243/09544070JAUTO1078MengX.XieY.A new numerical analysis for piston skirtliner system lubrication considering the effects of connecting rod inertia2012472352432-s2.0-8485625504110.1016/j.triboint.2011.12.013MengX.NingL.XieY.WongV. W.Effects of the connecting-rod-related design parameters on the piston dynamics and the skirt-liner lubrication201322768858982-s2.0-8488052155910.1177/0954407012464025MengX.FangC.XieY.Transient tribodynamic model of piston skirt-liner systems with variable speed effects2016946406512-s2.0-8494694383010.1016/j.triboint.2015.10.034ZhuJ.-C.ZhuH.-H.FanS.-D.XueL.-J.LiY.-F.A study on the influence of oil film lubrication to the strength of engine connecting rod components201663941052-s2.0-8495920950210.1016/j.engfailanal.2016.02.006PelosiM.IvantysynovaM.Heat transfer and thermal elastic deformation analysis on the Piston/Cylinder interface of axial piston machines201213442-s2.0-8486562910310.1115/1.4006980041101NingL.MengX.XieY.Incorporation of deformation in a lubrication analysis for automotive piston skirt-liner system201322766546702-s2.0-8488453663510.1177/1350650112466768JacobosonB.FlobergL.The finite journal bearing considering vaporization195719072OlssonK. O.Cavitation in dynamically loaded bearing196530860ZhaoY.WeiC.YuanS.HuJ.Theoretical and experimental study of cavitation effects on the dynamic characteristic of spiral-groove rotary seals (SGRSs)2016643, article no. 5010.1007/s11249-016-0784-62-s2.0-84994890057KangoS.SharmaR. K.PandeyR. K.Thermal analysis of microtextured journal bearing using non-Newtonian rheology of lubricant and JFO boundary conditions20146919292-s2.0-8488433303210.1016/j.triboint.2013.08.009LiuC.LuY.-J.ZhangY.-F.LiS.MüllerN.Numerical study on the lubrication performance of compression ring-cylinder liner system with spherical dimples20171272-s2.0-8502545211310.1371/journal.pone.0181574e0181574PatirN.ChengH. S.An average flow model for determining effects of three-dimensional roughness on partial hydrodynamic lubrication19781001121710.1115/1.34531032-s2.0-0017825563PatirN.ChengH. S.Application of average flow model to lubrication between rough sliding surfaces1979101222023010.1115/1.34533292-s2.0-0018457141WuC.ZhengL.Average Reynolds equation for partial film lubrication with a contact factor198911111881912-s2.0-002430447810.1115/1.3261872YangQ.KeithT. G.Two-dimensional piston ring lubrication—part i: rigid ring and liner solution19963947577682-s2.0-003026804710.1080/10402009608983593ZhuD.ChengH.AraiT.HamaiK.Numerical analysis for piston skirts in mixed lubrication. Part I: Basic modeling199211435535622-s2.0-002689657710.1115/1.2920917YinB.LiX.FuY.YunW.Effect of laser textured dimples on the lubrication performance of cylinder liner in diesel engine20122472933122-s2.0-8486757310810.1002/ls.1185GreenwoodJ. A.TrippJ. H.The contact of two nominally flat rough surfaces2016185162563310.1243/PIME_PROC_1970_185_069_02