Stress Analysis and Weight Reduction of A One-Cyli
Stress Analysis and Weight Reduction of A One-Cyli
Stress Analysis and Weight Reduction of A One-Cyli
net/publication/282462725
CITATION READS
1 1,929
2 authors, including:
János Égert
Széchenyi István University, Gyor
15 PUBLICATIONS 32 CITATIONS
SEE PROFILE
All content following this page was uploaded by János Égert on 17 November 2015.
Abstract: The paper deals with the weight reduction and the achievement of a more
homogeneous stress distribution of a one-cylinder engine’s crankshaft,
developed for the Formula Student series. In order to achieve this goal, a
series of static FEM analysis were performed on the crankshaft. For the
FEM analyses, the critical loading and kinematic boundary conditions of
the crankshaft were calculated with two different methods. The loads were
determined analytically and were verified by Creo Mechanism software.
The calculations were done for different engine speeds and - as a result -
the critical engine speed, where the highest loads appear could be
determined. An essential weight reduction and a more homogeneous stress
distribution were achieved on the crankshaft as the result of FEM analysis.
1. Introduction
The paper shows the stress analysis of the crankshaft of a one-cylinder engine,
designed by the SZEngine team (the stundent’s engine development team of Széchenyi
István University) for the Formula Student (FS) race series and gives recommendations
for the weight reduction of the crankshaft. The FS race tracks feature many curves,
slaloms and short, straight sections. These features must also be taken into consideration
in the course of designing and developing the proper engine for the race car. The main
objective of the engine design is the creation of an engine with the possible lowest
weight and the possible highest performance. In order to be able to reduce the weight of
the crankshaft, a series of static FEM analysis were performed on the crankshaft. On the
basis of the results gained by the analysis, one could reduce the weight of the
crankshaft, make the stress distribution more homogenous, and fix the emerging failures
in the construction. For the analysis, it is necessary to exactly know the maximum loads
acting on the crankshaft therefore – as the first step – the paper deals with this problem.
The second part shows the results of the stress analyses and the weight reduction.
201
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
The analyzed crankshaft (Fig. 1.) is a construction made of three parts, consisting of a
crankpin and two other parts, which includes a web, a counterweight and a half-shaft.
The three parts are assembled by an interference fit. The material of the crankshaft is
42CrMo4 alloy steel, while for balancing the crankshaft, tungsten hard metal rods are
used, which are also connected by an interference fit into the bores of the
counterweight. The material properties of the crankshaft are summarized in Table 1.
202
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
203
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
In Fig. 2. the vector r2 represents the connecting rod and r1 models the crank. On the
basis of Fig. 2., the following vector-equation can be written:
r1 r2 r3 0 . (1)
Differentiating once the above vector-equation with respect to time, we get the angular
velocity of the connecting rod, while differentiating it twice, the angular acceleration of
the connecting rod can be achieved.
1r1 cos
2 , (2)
r2 cos
The acceleration of point A is known therefore the acceleration of the arbitrary point P
of the connecting rod can be determined by using the angular velocity and the angular
acceleration:
aP a A 2 rAP 22 rAP , (4)
From the evaluation of each term (4), the acceleration of an arbitrary point P is:
aP (r112 sin 2 rAP cos 22 rAP sin ) ex
(5)
(r112 cos 2 rAP sin 22 rAP cos ) ey
Knowing the mass properties mh , J h of the connecting rod and the mass md of the
piston and the acceleration of an arbitrary point P specified above, the forces acting on
the connecting rod (Fig. 3.) can be determined. In Fig. 3. the notation S means the center
of gravity of the connecting rod.
204
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Applying the principle of angular momentum for the connecting rod rotating about the
point B:
J S 2 2 J S 2 rBS mh aS mh 2 rBS vB M B , (8)
where vB vB ey is the velocity of point B, which value can be determined in the
following way:
vB 1r1 sin 2 r2 sin . (9)
It is important to note that many of the literature sources, such as literature [3] gives
only the incomplete form of the principle of angular momentum for a rigid body, which
does not contain the fourth term on the left side of the equation (8). According to [6],
however, if the principle of angular momentum for a rigid body is applied to a moving
point B as in the present case – the principle must be completed by the fourth term
containing the velocity of the point B. In case of leaving this term, the x component of
the force acting in point A shows equation (10), while – taking into consideration this
term – it shows equation (11):
J h mh aSy (r2 uS ) sin mh aSx (r2 uS ) cos FAy r2 cos
FAx , (10)
r2 sin
Force acting on the crankpin of the crankshaft C will be the counter-force of the force
acting on point A of the connecting rod:
FCx FAx , FCy FAy . (13)
At 8000 rpm engine speed, the differently calculated tangential and radial components
of the force acting on the crankpin, were performed and compared with the results
achieved by Creo Mechanism program (Fig. 4. and 5.). In the course of comparison, a
proper match could be observed, but – as it can be seen in the figures – a slight
deviation could be noticed between the two analytical solutions, which are presented in
the tangential force diagram. It is interesting to observe that if the component in x
direction of force FA was calculated by applying the incomplete form of the principle
of angular momentum (10), an exact match with the result provided by Creo was to be
obtained.
205
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Figure 4. The tangential force acting on the crankpin in function of crank angle at
crankshaft speed of 8000 rpm with different calculation methods
Figure 5. The radial force acting on the crankpin in function of crank angle at
crankshaft speed of 8000 rpm with different calculation methods
In reality, the force does not act at a certain point, but on a cylindrical contact surface.
The crankpin is assumed to have a sinusoidal pressure distribution. Therefore, the force
acting on the crankpin – that was determined above – was transformed into a sinusoidal
pressure distribution acting over an angle of 120 (Fig. 6.).
206
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
The pressure curve along the contour of crankpin can be prescribed by the following
formula:
180 180
p( ) pmax sin max 90 , (14)
where the value 0o belongs to the t axis.
In (14) the pressure constant pmax can be achieved by integrating the pressure surface:
dC
360 wC
360 180FC
FC pmax sin d s dz ds
0 0 C
pmax
wC dC
, (15)
Figure 7. The loads acting on the crankpin in the function of crank-angle at both ends
of the operational range
The FEM analysis should be performed at the engine speed, at which the load acting on
the crankpin has the maximum value. Fig. 8. shows that the maximum occurs at 6500
rpm engine speed.
Figure 8. The maximum values of the resultant forces acting on the crankpin in the
function of engine speed
207
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
2.3. The forces of inertia resulting from the speed irregularity of the crankshaft
This section discusses the instantaneous angular acceleration, which is necessary for
the calculation of the body forces resulting from the rotation of the crankshaft, as well
as for calculating the torque originating from the flywheel.
The tangential force acting on the crankpin is continuously changing in the function
of the crank angle, as it can be seen in Fig. 4. Since the torque of the engine is provided
by this force, the torque of the engine will also be irregular, so it will depend on the
crank angle.
The inertia of the crankshaft is increased by the flywheel mounted at the end of the
crankshaft, which also serves as an energy storing unit. When the instantaneous torque
is higher than the medium torque, the crankshaft accelerates, while the flywheel picks
up energy. When the instantaneous torque is lower than the medium torque, the
crankshaft slows down, while the flywheel delivers energy back to the crankshaft. In the
majority of engines, the flywheel – which also includes the clutch – is mounted on one
side of the crankshaft. In case of the analyzed single-cylinder engine, the clutch and the
flywheel are mounted as two separate parts at both ends of the crankshaft.
In order to be able to determine the speed irregularity, as a first step the medium
torque must be calculated. This is a torque value supposed to be constant, which
performs the same amount of work as the varying (actual) torque. For this, the torque
curve must be integrated for a whole engine cycle and the height of a rectangle of the
same width provides the medium value on the basis of the equality of the surfaces.
Fig. 9. illustrates the instantaneous torque and the medium torque at crankshaft speed
of 6500 rpm. The work below the medium torque – marked by red stripes – is negative
and this decreases the angular velocity, while the work above the medium torque –
marked by green stripes – is positive thus increasing the angular velocity.
208
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
n
E n E 0 M M K d , (16)
0
where E is the energy, resulting from the rotation of the crankshaft, clutch and flywheel
and M K is the medium torque. Knowing the instantaneous total energy of the system,
the instantaneous angular velocity of the crankshaft can easily be determined:
2 E ( )
1 ( ) . (17)
J f Jt Jl
Figure 10. Variation of crankshaft speed at initial crankshaft speed of 6500 rpm
The flywheel – with its high inertia – tends to decrease the irregularity during the
work of the engine. This balancing effect causes high torque on the crankshaft. The
torque, resulting from the inertia of the flywheel can be calculated by multiplying the
inertia of the flywheel with the instantaneous angular acceleration, which is the
derivative of the instantaneous angular velocity with respect to time:
d 1
M t Jl1 Jl . (18)
dt
Figure 11. Variation of torque caused by the flywheel at crankshaft speed of 6500 rpm
209
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
210
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Figure 15. Maximum von Mises stresses occurring as a result of interference fits and
the thermal load
211
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Figure 16. Variation of maximum von Mises stresses caused by the individual loads in
function of crank angle
Fig. 17. shows the stress distribution of each load at the crank angle, where the
maximum von Mises stress occur. On the basis of the results, the following conclusions
could be drawn:
The maximum equivalent von Mises stress acting on the crankpin occurs at 371o
crank angle, this means around the ignition its magnitude is about 300 MPa. The
maximum von Mises stresses occur between the fillet of the web and at the main
journal (Fig. 17.a).
The maximum von Mises stress caused by the forces of inertia varies according to
the angular velocity. The stresses occurring here are much lower compared to
stresses caused by former two loads (max. 70 MPa). The highest stresses are
occurring between the main journal and the counterweight (Fig. 17.b).
The maximum von Mises stress occurring as a result of the torque caused by the
flywheel varies according to the variation of the torque. The highest stresses
occur on the surface of the half shafts (Fig. 17.c).
212
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Figure 18. The maximum von Mises and radial stresses in case of different loads and
different mesh density
213
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Figure 20. Variation of the minimum radial stress (maximum contact pressure) in
function of crank angle in case of different interferences
In case of interference of 0,02 mm, the minimum value of the contact pressure around
ignition reaches zero value, so there will be disconnection between the crankpin and the
bore, it is to be feared that the crankpin will become loose, therefore at least 0,03 mm
interference is necessary in order to be able to properly ensure the interference fit
(Fig. 20.). As there is no torsional or axial load between the crankpin and the bore of the
web, sliding won’t occur between connecting components, so we didn’t take it into
consideration.
214
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
215
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
Figure 22: Von Mises stress distribution of the modified crankshaft caused by critical
total load
Summary
For the analysis, the critical loads acting on the crankshaft were analytically
determined. In the course of this, we have stated that the highest load occurs on the
crankshaft at a speed of 6500 rpm, therefore the series of static FEM analysis were
performed at this engine speed. The results of the FEM analysis have shown that the
crankshaft has a constructional failure, because the equivalent von Mises stress caused
by the interference between the crankpin and the web was significantly higher than the
stresses caused by the combustion. As a first step, this failure was eliminated by
decreasing the extent of the interference. From the FEM analysis it can also be
concluded, that there are several areas on the crankshaft, at which very small stresses
occur, so material was removed from these areas. As a result of the constructional
modifications, the stress distribution became more homogenous furthermore the weight
of the crankshaft could be reduced by 19 %.
References
[1] Dezsényi Gy: Design and testing of internal combustion engines (in Hungarian).
Nemzeti Tankönyvkiadó, Budapest, 1989.
216
P. Horváth and J. Égert – Acta Technica Jaurinensis, Vol. 8, No. 3, pp. 201-217, 2015
217