Finite Element Analysis of a Rocker Arm

Finite Element Analysis of a Rocker Arm

 

Component Description

 

1.1.            Component Function

The rocker arm is an oscillating,  two-arm lever that provides a means of actuating the valves in the combustion chamber of an internal combustion engine. It translates the radial motion of the profile of the cam lobe through a fulcrum into linear motion for opening and closing the intake and exhaust valves.  It also provides a means of multiplying the lift ratio.

During operation, the rocker arm experience stresses and undergo deflection. Severe rocker arm deflection causes inefficient engine performance, and often results in metal fatigue leading to increased wear and friction in the valve train and eventually engine failure.

Figure 1: Diesel engine valve train showing rocker arm (1)

1.2.            Component Geometry

The rocker arm is a two-arm lever that pivots about a fulcrum. One end is connected to the push rod which rests over cams on the camshaft, while the other acts on the spring-loaded valve stem and pivoted on the rocker shaft.

Figure 2:(From left) 3-D model of the rocker arm and the 2-D geometry

 

1.3.            Service Loading

The rocker arm is subjected to compressive load at the fulcrum on the rocker shaft during the opening and closing of the valves.  The forces acting at the valve end (FE) include; the gas back pressure on the valve, the spring force and the force due to valve acceleration. There is also the load at the cam end, (FC) which is transmitted to the rocker arm through the push rod.

Figure 3: Free Body Diagram of Rocker arm depicting the location of the acting forces: FE and FC

Methodology (154)

         Two-Dimensional Idealization

In a rocker arm, the constraint at the pivot hole and the applied loads at both ends, all act in a plane, parallel to the cross-section plane of the rocker arm.

To estimate the stress concentration of the geometry, it is idealized as a shell and shaped in the profile of the planar cross section of the rocker arm under the assumption that the rocker arm is completely solid with only the holes drilled through as seen in Figure 4. The two-dimensional idealization makes it easier to deal with and the main stresses are obtained with reasonable accuracy.

Figure 4: 2D geometrical Idealization using Autodesk Inventor

2.2.            Loading, Boundary Conditions and Constraints

The rocker arm has two major loads applied on either end respectively:

Force due to exhaust valve loading; FE

Force due to cam action through the push rod; FC

The main forces will be calculated based on the Diesel Engine specifications below.

Type of Engine

Turbocharged 6-cylinder, Diesel Engine, 2523CCm2DiCR Engine.

Number of Cylinders

6

Engine capacity

2523 CC

Maximum Engine Power

46.3 kW @ 3200 rpm

Maximum Torque

195 Nm @ 1440 – 2200 rpm.

Table 1: Diesel Engine Specification (2)

Mass of the valve mv

0.09kg

Diameter of the valve head dv

40mm

Lift of the valve   h

9.4mm

Cylinder pressure Pc

0.4N/mm2

Maximum suction pressure Ps

0.02N/mm2

Diameter of fulcrum pin d1

22mm

Diameter of boss D1

34mm

Rocker arm ratio

1.64

Engine Speed

3200 RPM

Angle action of cam (Ɵ)

110o

Spring rate k

23 N/mm

Spring Preload P1

249.5N

Weight of associates parts with valve w

0.882 N

Acceleration of valve a

1550 m/s2

Table 2: Input data for calculation of applied load on rocker arm (2)

Total load on valve:
Pt=Pg+w=π4×d12×Pc+w=502.4N

Initial spring force:
Fs
=  
π4×d12×Pc–w
= 24.249 N

Force due to valve acceleration with valve weight:
Fa=m×a–w=140.38N
  

Maximum load on the rocker arm for exhaust valve:
Fe=P+Fs+Fa+(P1+k×h)=1133.87N

Load on push rod side of the rocker arm:
 Fc=Fe×Rocker arm ratio=1859.54N

A pinned constraint is applied on the pivoting hole in all degrees of freedom.

Finite Element Representation

         Modelling and Justification of Shell type

The material for the rocker arm; carbon steel is assumed to be homogeneous and isotropic with details in Table 3.

Young’s Modulus (GPa)

200

Poisson’s ratio

0.29

Density (kg/m3)

7850

Yield Strength (MPa)

415

Ultimate Strength (MPa)

585

Table 3: Carbon Steel Material Properties

Plane stress idealization is assumed because the thickness in the z-direction is smaller than the length in the x-y plane. Given the 2D geometry, it is idealized as a plane stress model and shell elements must be used. The shape of the element can either be triangular or quadrilateral and the element order is either linear or parabolic.

Element Size

Maximum von-Mises Stress

(MPa)

(Triangular element)

Maximum von-Mises Stress (MPa)

(Quadrilateral element)

Linear

Parabolic

Linear

Parabolic

1.0mm

97.51

140.7

141.9

167.8

0.8mm

109.2

149.0

137.6

167.9

0.6mm

118.9

162.3

143.6

176.5

0.4mm

144.9

182.3

170.5

179.4

0.2mm

162.6

185.6

178.6

185.4

Table 4: Comparison of Element type and order

The element type chosen for the meshing is the second order quadrilateral element. This was due to the accuracy and consistency of the von Mises stress values obtained when the Force FE is applied on the rocker arm in comparison to other element types. Therefore, it has proven to be the best for this geometry.

To accurately represent the pinned pivot constraint, a pinned, rigid body connector was applied at the centre point and connected to the pivot hole so all the nodes at the hole can better represent the motion.  For each loading scenario, one hole must have a fixed constraint the other has the load applied for the solution to be solved. This can be seen in Figure 5 and Figure 6.

Figure 5: Mesh diagram for FC loading scenario

Figure 6: Mesh diagram for FE loading scenario

Contour Plots

After the finite element model was solved, the maximum Von-Mises stress was recorded. The finite element contour plot for each loading scenario can be seen in Figure 7 and Figure 8.

Figure 7: FE loading scenario

Figure 8: FC loading scenario

When Force FC was applied, it predicted the largest maximum von-Mises stress compared to when Force FE was applied. Hence, FC is the worst-case loading scenario and the stress obtained from this scenario will be evaluated.

4.1.            Justification of results parameter

The von Mises stress is chosen as the main result parameter to be evaluated. This is required to evaluate the failure rate of the rocker arm, since it is made of a ductile material and is largely likely to fail as a result of shearing action. The von Mises stress yield criterion predicts failure for a ductile material, when the von Mises stress becomes equal to the yield strength.

The contour plot predicts the maximum von Mises stress value as 203.3MPa, occurring at the neck of the rocker arm.  This value is within the Yield strength for carbon steel of 415MPa.

Mesh Convergence

In order to explore the effects of mesh density on the prediction from the finite element model, several finite element models were generated, varying the number of elements from 39 (coarsest mesh) to 12069 (finest mesh) from element sizes ranging from 0.2mm to 20.5 mm. The von Mises stress from the worst-case scenario were calculated and compared, to study the influence of mesh density on the stresses predicted, as shown below.

Number of elements

39

84

500

1019

2364

3045

3970

5368

9189

12069

Von Mises stress (MPa)

107.4

122.5

189.6

195.6

205.6

203.1

205.7

204.5

204.4

203.3

Table 5: Effect of mesh density

From Figure 9, the stress converges to the maximum possible value of 203.3 MPa.

Conclusion

In this case, them model produced acceptable results. The maximum stress was 203.3 MPa and was localized mostly and the neck of the rocker arm which is where the stress was expected to be most concentrated. The convergence of the results shows that the result is not significantly dependent on the choice of discretization and is an indication of the validity of the accuracy of the prediction. Also, the result obtained under both loading conditions, indicating the same location (the neck of the rocker arm) as where the maximum von Mises stress is predicted and where failure is likely to occur provides further validation of the result. However, the component was assumed to be completely solid except for the holes. In reality, the 3D geometry (Figure 2) is a better representation of the component and the aforementioned assumption is a relatively poor one. Apart from the applied loads, the component would be subjected to thermal stresses. Furthermore, though plane stress was assumed, the thickness of 19 mm was relatively large compared to the length of 50.8 mm. To validate the model, plane strain assumptions may be made during the next analysis in which a displacement-controlled test can be used to identify the strain that causes failure.

References

1.) Diesel Engine Valve Train. Nuclear Power Training. [Online] [Cited: November 14th, 2018 at 21:34.] http://nuclearpowertraining.tpub.com/h1018v1/css/Figure-10-Diesel-Engine-Valve-Train-31.htm.

2.) DESIGN AND STATIC STRUCTURAL ANALYSIS OF ROCKER ARM IN I.C. Bacha, Sachin et al. s.l. : INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH, 2018, INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH, p. 9.