MECH2610 - S2 - Unit 2 - Lecture Notes - Full
MECH2610 - S2 - Unit 2 - Lecture Notes - Full
MECH2610 - S2 - Unit 2 - Lecture Notes - Full
2.1 Introduction:
In this unit, we will extend the stress analysis of symmetric beams under pure bending to that of the
asymmetric bending. If a bending moment about the z axis (M(z)) is applied on a symmetric beam, the
beam only bends in the xoy plane. Thus, the direct stress along the beam (in the longitudinal direction,
x) is caused by this bending deflection. However, if a pure bending moment is applied on an
asymmetric beam, the beam bends in both xoy and xoz planes, Moreover, even if a beam is
geometrically symmetric, when it is subjected to a bending moment about the z axis (M(z)) and another
about the y axis (M(y)), then the beam will also bend in both xoy and xoz planes. The bending in both
planes will result in the direct stress along the beam (x). Therefore, calculating this longitudinal
direct (normal) stress along the beams under asymmetric bending is more difficult, and the aim of
this unit is to solve this problem.
In this unit, we will first revise the stress analysis of symmetric beams under pure bending. Then the
same idea will be extended to asymmetric bending. The key idea is that the total direct stress along
the beam is the sum of the two stresses (both in the longitudinal direction) due to the deformation
in both planes. We will end up with an equation that describes the distribution of the direct stress on
the cross section of the beams subjected to asymmetric bending. Using this equation, the maximum
tensile and compressive direct stresses in the beam, as well as their positions, can be analysed.
2.2.1 Neutral axis: Consider a beam subjected to a pure bending moment M (no forces applied on
the beam!) about a symmetric axis (the z axis in this case). In Figure 2.1, the cross section is symmetric
(about the y axis) and the applied moment is about the z axis. The upper surface (AB) shortens and
thus is in compression. The lower surface stretches and therefore in tension. Hence there must be an
xoz plane in between in which the longitudinal deformation is zero. This plane is termed the neutral
plane. A transverse axis lying in the neutral plane is the neutral axis.
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.1
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
Figure 2.1
2.2.2 Longitudinal strain: Consider the deformation between two sections AC and BD (Figure 2.1).
The initial distance between AC and BD is dx. A longitudinal fibre GH, which is y distance below the
neutral axis, initially has the same length as the neutral axis, dx. During bending, GH becomes GH,
but the length of the neutral axis remains the same. If R is the radius of curvature of the neutral axis
and d is the angle between AC and BD, then Rd is the length of the neutral axis after deformation.
Because the length of the neutral axis remains the same during bending, Rd = dx.
GH = (R+y) d
Then the longitudinal strain of fibre GH is:
GH − GH (𝑅 + 𝑦)d − d𝑥 𝑅d + 𝑦d − 𝑅d
𝒙 = = =
𝐺𝐻 d𝑥 𝑅d (2.1)
𝒚
=
𝑹
We know from equation (2.1) that under pure bending, the strain of a symmetric beam
(1) is distributed linearly across the section;
(2) is zero at the neutral plane;
(3) has the maximum value at the outer surfaces.
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.2
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
(3) has the maximum value at the outer surfaces: being maximum tension and compression at the
two outer surfaces where y is maximum (Figure 2.2).
Figure 2.2
Figure 2.3
Consider an element of area, dA, at a distance y from the neutral plane (Figure 2.3). The internal force
on the element caused by the longitudinal direct stress in the x direction is
dFx = 𝑥 d𝐴
Therefore, the total longitudinal force on the cross section is
𝐹𝑥 = ∫ 𝑥 d𝐴
𝐴
𝐹𝑥 = ∫ 𝑥 d𝐴 = 0
𝐴
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.3
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
From equation (1.9), we can calculate the position of the centroid of the cross section
1
𝑦̅ = ∫ 𝑦d𝐴 = 0
𝐴
In other words, the neutral axis must pass through the centroid of the cross section.
∫ 𝑦𝑥 d𝐴
𝐴
To satisfy the equilibrium of the moment, this total internal resisting moment must balance the applied
external moment M.
∫ 𝑦𝑥 d𝐴 = 𝑀
𝐴
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.4
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
the beam is not symmetric. External moments My and Mz are applied about an arbitrary set of centroidal
axes. Therefore, the beam will bend in both xoy and xoz planes. The bending in both planes will
cause strain and stress in the longitudinal (x) direction.
Figure 2.4
Figure 2.5
From section 2.2, we have known that the longitudinal direct strain ( 𝑥 ) and stress (𝑥 ) due to the
bending in the xoy plane is 𝑦⁄𝑅𝑧 and 𝐸𝑦⁄𝑅𝑧 , respectively. Similarly, the longitudinal strain and stress
due to the bending in the xoz plane is 𝑧⁄𝑅𝑦 and 𝐸𝑧⁄𝑅𝑦 , respectively. Therefore, the total longitudinal
direct stress on the cross section should be the superposition of the stresses caused by the bending
in both planes:
𝑬 𝑬
𝒙 = 𝒚+ 𝒛 (2.6)
𝑹𝒛 𝑹𝒚
2.3.2 Equilibrium of moment
Consider an element of area dA on the cross section (Figure 2.4). The moment of the axial force on the
element about the z axis is
ydFx = y𝑥 d𝐴
Then the resultant total internal resisting moment about the z axis is
∫ 𝑦𝑥 d𝐴
𝐴
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.5
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
To satisfy the equilibrium of moment, this total internal resisting moment must balance the applied
external moment Mz.
∫ 𝑦𝑥 d𝐴 = 𝑀𝑧
𝐴
where ∫𝐴 𝑦 2 d𝐴 is the second moment of area about the z axis (Iz). ∫𝑨 𝒚𝒛 𝐝𝑨 is the product moment
of area of the cross section, Iyz. More details about Iyz can be found in sections 2.3.6 and 2.4.
So
𝑬𝑰𝒛 𝑬𝑰𝒚𝒛
+ = 𝑴𝒛 (2.7)
𝑹𝒛 𝑹𝒚
Consider the element of area dA again (Figure 2.4). The moment of the longitudinal force on the
element about the y axis is
𝑧d𝐹𝑥 = 𝑧𝑥 d𝐴
Then the resultant total internal resisting moment about the y axis is
∫ 𝑧𝑥 d𝐴
𝐴
To satisfy the equilibrium of moment, this total internal resisting moment must balance the applied
external moment My.
∫𝐴 𝑧𝑥 d𝐴 = 𝑀𝑦
Substitute 𝑥 by equation (2.6)
𝐸 𝐸
∫ 𝑧( 𝑦+ 𝑧)d𝐴 = 𝑀𝑦
𝐴 𝑅𝑧 𝑅𝑦
𝐸 𝐸
∫ 𝑧𝑦d𝐴 + ∫ 𝑧 2 d𝐴 = 𝑀𝑦
𝑅𝑧 𝐴 𝑅𝑦 𝐴
Where ∫𝐴 𝑧 2 d𝐴 is the second moment of area about the y axis (Iz). So
𝑬𝑰𝒚𝒛 𝑬𝑰𝒚
+ = 𝑴𝒚 (2.8)
𝑹𝒛 𝑹𝒚
2.3.3 𝑬⁄𝑹𝒛 and 𝑬⁄𝑹𝒚
Solving equations (2.7) and (2.8) gives
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.6
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
𝑬 𝑴𝒛 𝑰𝒚 − 𝑴𝒚 𝑰𝒚𝒛
= (2.9)
𝑹𝒛 𝑰𝒚 𝑰𝒛 − 𝑰𝒚𝒛 𝟐
𝑬 𝑴𝒚 𝑰𝒛 − 𝑴𝒛 𝑰𝒚𝒛 (2.10)
= 𝟐
𝑹𝒚 𝑰𝒚 𝑰𝒛 − 𝑰𝒚𝒛
𝑀𝑧 and 𝑀𝑦 are the moments applied on the beam, which are normally known. 𝐼𝑦 , 𝐼𝑧 , and 𝐼𝑦𝑧 are the
area properties of the cross section of the beam, which are determined by the dimensions of the cross
section. So 𝐸 ⁄𝑅𝑧 and 𝐸 ⁄𝑅𝑦 can be calculated if we know the moments applied and the geometry of
the beam.
Note: (1) if the beam is asymmetric, but only 𝑀𝑧 is applied (𝑀𝑦 = 0), then both 𝐸 ⁄𝑅𝑧 and 𝐸 ⁄𝑅𝑧 is
not zero. It is still asymmetric pure bending. (2) if the beam is symmetric, then 𝐼𝑦𝑧 is zero. However,
if both 𝑀𝑧 and 𝑀𝑦 are applied, then both 𝐸 ⁄𝑅𝑧 and 𝐸 ⁄𝑅𝑧 is not zero. This is also an asymmetric pure
bending case. (3) if the beam is symmetric but only 𝑀𝑧 is applied (𝑀𝑦 = 0), then 𝐸 ⁄𝑅𝑦 = 0. This is
equivalent to the symmetric pure bending in Section 2.2.
Figure 2.6
∫ 𝑦𝑧 d𝐴 = ∬ 𝑦𝑧 d𝑧d𝑦
𝐴
Because this is a rectangle (z and y are independent of each other), the above equation can be written
as
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.7
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
𝑧max 𝑦max
∬ 𝑦𝑧 d𝑥d𝑦 = ∫ 𝑧 d𝑧 ∫ 𝑦 d𝑦
𝑧min 𝑦min
Integrating it gives
1
𝐼𝑦𝑧 = (𝑦max 2 − 𝑦min 2 )(𝑧max 2 − 𝑧min 2 ) (2.11)
4
Figure 2.7
Step 1: Calculating the centroid.
This is important because we need Iz, Iy, and Iyz to calculate 𝐸 ⁄𝑅𝑧 and 𝐸 ⁄𝑅𝑦 . Iz, Iy, and Iyz are all
about the centroidal axes. Therefore, in order to calculate Iz, Iy, and Iyz, we need to know the position
of the centroid.
Using equating (1.10),
𝑦̅ ∑ 𝐴 = ∑ 𝐴𝑦̃ ; 𝑧̅ ∑ 𝐴 = ∑ 𝐴𝑧̃
90 10
Component A: 𝐴𝐴 = 9010 mm2, 𝑦
̃𝐴 = = 45 mm, 𝑧̃𝐴 = = 5 mm
2 2
60
Component B: 𝐴𝐵 = 6010 mm2 , 𝑦̃
𝐵 = 90 + 5 = 95 mm, 𝑧
̃𝐵 = = 30 mm
2
𝑦̅ (9010 + 6010) = (901045) + (601095) (𝑦̃ was calculated from the top of the cross
section. Of course it could also be calculated from the bottom.)
𝑧̅(9010 + 6010) = (𝑦 ̃𝐴 = 45 mm 90105) + (601030) (𝑧̃ was calculated from the left of
the cross section. It could also be calculated from the right.)
(901045) + (601095)
𝑦̅ = = 65 mm
9010 + 6010
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.8
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
(90105) + (601030)
𝑧̅ = = 15 mm
9010 + 6010
Step 2: Calculating Iz, Iy, and Iyz.
2a: After we obtain the position of the centroid, we locate a global coordinate system (yoz) at the
centroid to help us calculate the second moment of area. The following calculation is all in this new
global coordinate system (we don’t use the reference point at the left top corner anymore after
we got the position of the centroid).
2b: Then we find the coordinates (𝑦̂, 𝑧̂) of the local centroids of the two components in the global
coordinate system. The coordinates of the local centroids of the two components are (4565, 515)
and (9565, 3015), respectively.
2c: 𝑰𝒛 = 𝑰𝑨𝒛 + 𝑰𝑩
𝒛
2 10903
𝐼𝑧𝐴 = 𝐼𝑧𝐴′ + 𝐴𝐴 𝑦
̂𝐴 = + 9010(45 − 65)2 = 967500 mm4
12
60103
̂𝐵 2 =
𝐼𝑧𝐵 = 𝐼𝑧𝐵′ + 𝐴𝐵 𝑦 + 6010(95 − 65)2 = 545000 mm4
12
𝐼𝑧 = 𝐼𝑧𝐴 + 𝐼𝑧𝐵 = 1512500 mm4
2d: 𝑰𝒚 = 𝑰𝑨𝒚 + 𝑰𝑩
𝒚
90103
𝐼𝑦𝐴 = 𝐼𝑦𝐴′ + 𝐴𝐴 𝑧̂𝐴 2 = + 9010(5 − 15)2 = 97500 mm4
12
10603
2
𝐼𝑦𝐵 = 𝐼𝑦𝐵′ + 𝐴𝐵 𝑧̂
𝐵 = + 6010(30 − 15)2 = 315000 mm4
12
𝐼𝑦 = 𝐼𝑦𝐴 + 𝐼𝑦𝐵 = 412500 mm4
𝐴 𝐵
In order to use equation (2.11) to calculate 𝐼𝑦𝑧 and 𝐼𝑦𝑧 , we need to find the ymin, ymax, zmin and zmax for
both components. These are the points on the component that are furthest from the centroidal
axes. These values take into account +/- signs depending on their position compared to the
centroidal axes. Consider ymin and ymax of component A. We are considering their position relative
to the centroidal axes. The uppermost point on the component is 65 mm in the Y- direction and
the lowermost point is 25mm in the Y+ direction. Therefore the value for ymin = 65mm and the
value for ymax = +25mm. Similarly, zmin = 15 mm ; zmax = 5 mm.
For component B: ymin = 25 mm; ymax = 35 mm; zmin = 15 mm ; zmax = 45 mm
𝐴
1 1
𝐼𝑦𝑧 = (𝑦max 2 − 𝑦min 2 )(𝑧max 2 − 𝑧min 2 ) = (252 − 652 )(52 − 152 ) = 180000 mm4
4 4
𝐵
1 1
𝐼𝑦𝑧 = (𝑦max 2 − 𝑦min 2 )(𝑧max 2 − 𝑧min 2 ) = (352 − 252 )(452 − 152 ) = 270000 mm4
4 4
𝐴 𝐵
𝐼𝑦𝑧 = 𝐼𝑦𝑧 + 𝐼𝑦𝑧 = 450000 mm4
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.9
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
Calculate the stresses at these two points and compare them to get the maximum compressive stress.
At the first point,
𝑥 = 1.96(−65) − 2.14(−5) = −116.7 MPa
At the second point,
𝑥 = 1.9625 − 2.1445 = −47.3 MPa > −116.7 MPa
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.10
MECH2610 – Engineering Mathematics –2nd Semester Unit 2
So the maximum compressive stress occurs at Point 1 in the above figure (top right corner of
component A).
Note: a common mistake here is to choose (y = 65 mm, z = 45mm) to calculate the maximum
compressive stress. But is this point on the beam?
Dr Qingen Meng Copyright2023 University of Leeds UK. All rights reserved. 2.11