PERIDO, JOHN PAULO M.

AENG 260 -EXPERIMENTAL STRESS ANALYSIS

MSAE – FPM PROBLEM SET IV
2018-68796 DR. DELFIN C. SUMINISTRADO

3.31. A steel (E = 207 GPa and v = 0.30) ring is shrunk onto another steel ring, as shown in the

figure below. Determine the maximum interference possible without yielding one of the rings if

the yield strength of both rings is 900 MPa and the dimensions of the rings are as follows:

θ Case 1 Case 2 Case 3 Case 4

𝒂𝒊 , 𝒎𝒎 100 100 100 0
𝒃𝒊 , 𝒎𝒎 125 115 150 100
𝒃𝒐 , 𝒎𝒎 150 200 175 150

Interference between two different materials

𝒓 𝒓𝒐 𝟐 + 𝒓𝟐 𝒅 𝒓𝟐 + 𝒓𝒊 𝟐
𝜹 = 𝝈[ ( 𝟐 + 𝒗 𝒐 ) + ( − 𝒗𝒊 )]
𝑬 𝒐 𝒓𝒐 − 𝒓𝟐 𝑬 𝒊 𝒓𝟐 − 𝒓𝒊 𝟐

Where:

𝛿 - interference, mm

𝜎 - pressure generated with the designed connection or yield strength, MPa

𝐸𝑜 - Modulus of elasticity of the outer materials, GPa

𝐸𝑖 - Modulus of elasticity of the inner materials, GPa

 𝑟𝑜 - outer radius of the larger ring, mm

𝑟𝑖 - inner radius of the smaller ring, mm

𝑟 - inner radius of the larger ring or outer radius of the smaller ring, mm

𝑣𝑖 - Poisson’s ratio of the inner ring

𝑣𝑜 - Poisson’s ratio of the outer ring

Interference of same the material

𝟐𝝈𝒓𝟑 𝒓𝒐 𝟐 − 𝒓 𝒊 𝟐
𝜹= [ ]
𝑬 (𝒓𝒐 𝟐 − 𝒓𝟐 )(𝒓𝟐 − 𝒓𝒊 𝟐 )

Where:

𝛿 - interference, mm

𝜎 - pressure generated with the designed connection or yield strength, MPa

𝐸𝑖 - Modulus of elasticity of the inner materials, GPa

𝑟𝑜 - outer radius of the larger ring, mm

𝑟𝑖 - inner radius of the smaller ring, mm

𝑟 - inner radius of the larger ring or outer radius of the smaller ring, mm

CASE 1.

𝟐(𝟗𝟎𝟎 𝑴𝑷𝒂)(𝟏𝟐𝟓 𝒎𝒎)𝟑 (𝟏𝟓𝟎 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐

𝜹= [ ]
𝟐𝟎𝟕 𝑮𝑷𝒂 ((𝟏𝟓𝟎𝒎𝒎)𝟐 − (𝟏𝟐𝟓 𝒎𝒎)𝟐 )((𝟏𝟐𝟓 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐 )
 𝜹 = 𝟓. 𝟒𝟖𝟗𝟕 𝒎𝒎

CASE 2.

𝟐(𝟗𝟎𝟎 𝑴𝑷𝒂)(𝟏𝟏𝟓 𝒎𝒎)𝟑 (𝟐𝟎𝟎 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐

𝜹= [ ]
𝟐𝟎𝟕 𝑮𝑷𝒂 ((𝟐𝟎𝟎𝒎𝒎)𝟐 − (𝟏𝟏𝟓 𝒎𝒎)𝟐 )((𝟏𝟏𝟓 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐 )

𝜹 = 𝟒. 𝟓𝟗𝟒𝟕 𝒎𝒎

CASE 3.

𝟐(𝟗𝟎𝟎 𝑴𝑷𝒂)(𝟏𝟓𝟎 𝒎𝒎)𝟑 (𝟏𝟕𝟓 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐

𝜹= [ ]
𝟐𝟎𝟕 𝑮𝑷𝒂 ((𝟏𝟕𝟓 𝒎𝒎)𝟐 − (𝟏𝟓𝟎 𝒎𝒎)𝟐 )((𝟏𝟓𝟎 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐 )

𝜹 = 𝟓. 𝟗𝟓𝟗𝟗 𝒎𝒎

CASE 4.

𝟐(𝟗𝟎𝟎 𝑴𝑷𝒂)(𝟏𝟎𝟎 𝒎𝒎)𝟑 (𝟏𝟓𝟎 𝒎𝒎)𝟐 − (𝟎 𝒎𝒎)𝟐

𝜹= [ ]
𝟐𝟎𝟕 𝑮𝑷𝒂 ((𝟏𝟓𝟎 𝒎𝒎)𝟐 − (𝟏𝟎𝟎 𝒎𝒎)𝟐 )((𝟏𝟎𝟎 𝒎𝒎)𝟐 − (𝟎 𝒎𝒎)𝟐 )

𝜹 = 𝟏. 𝟓𝟔𝟓𝟐 𝒎𝒎

Reference

Press Fit Pressure Calculator – Optimize Your Interference/Transition Fit Design (2012, October
0000021). Retrieved from https://www.meadinfo.org/2009/07/press-fit-pressure- calculator-
00000optimize.html.
3.16. Determine the stresses in the uniformly loaded cantilever beam shown in the figure below

by using Airy’s stress function approach. Assume a unit thickness.

Boundary conditions

At x = 0;

𝜎𝑥 = 0

𝜏𝑥𝑦 = 0

At x = L;

ℎ
∫ 𝜏𝑥𝑦 𝑏𝑑𝑦 = 𝑞𝐿
−ℎ

ℎ
∫ 𝜎𝑥 𝑏𝑑𝑦 = 0
−ℎ

ℎ
1 2
∫ 𝜎𝑥 𝑦𝑏𝑑𝑦 = 𝑞𝐿
−ℎ 2

At y = h;

𝑞
𝜎𝑦 = −
𝑏

𝜏𝑥𝑦 = 0
At y = -h;

𝜎𝑦 = 0

𝜏𝑥𝑦 = 0

According to Budynas (1999), Condition 𝜎𝑥 = 0 will be ignored however, at 𝑥 = 0, the

ℎ ℎ
conditions that ∫−ℎ 𝜎𝑥 𝑏𝑑𝑦 = 0 and ∫−ℎ 𝜎𝑥 𝑦𝑏𝑑𝑦 = 0 must hold for equilibrium. Some observations

should be made before attempting to establish the stress function. Since the cross section is
𝑞
symmetric about the y axis and 𝜎𝑦 = − 𝑏 at y = h and 𝜎𝑦 = 0 at y = -h, 𝜎𝑦 should be an odd

function of y. Thus, the stress function should contain odd functions of y. This is substantiated by

the fact that 𝜎𝑦 should also contain odd functions of y since the net axial force is zero (which

ℎ
satisfy condition ∫−ℎ 𝜎𝑥 𝑏𝑑𝑦 = 0). Second, since 𝜎𝑦 is constant as a function of x on the top and

bottom faces, the stress function should not contain powers of 𝑥 greater than 𝑥 2 . Thus the

following stress function will be tried:

𝛷 = 𝐴𝑥𝑦 + 𝐵𝑥 2 + 𝐶𝑥 2 𝑦 + 𝐷𝑦 3 + 𝐸𝑥𝑦 3 + 𝐹𝑥 2 𝑦 3 + 𝐺𝑦 5

Since ∇4 𝛷 = 0,

∇4 𝛷 = 24𝐹𝑦 + 120𝐺𝑦 = 0

Thus

𝐹 = −5𝐺

Equations of Stresses

𝜕 2𝛷
𝜎𝑥 = = 6𝐷𝑦 + 6𝐸𝑥𝑦 − 30𝐺𝑥 2 𝑦 + 20𝐺𝑦 3
𝜕𝑦 2
 𝜕 2𝛷
𝜎𝑦 = = 2𝐵 + 2𝐶𝑦 − 10𝐺𝑦 3
𝜕𝑥 2

𝜕 2𝛷
𝜏𝑥𝑦 = = −(𝐴 + 2𝐶𝑥 + 3𝐸𝑦 2 − 30𝐺𝑥𝑦 2 )
𝜕𝑥𝜕𝑦

For 𝜏𝑥𝑦 at 𝑥 = 0,

𝐴 + 3𝐸𝑦 2 = 0

this can be true for all values of y when,

𝐴=𝐸=0

For 𝜏𝑥𝑦 at 𝑦 = ±ℎ and equating to zero,

0 = −(2𝐶𝑥 − 30𝐺ℎ2 𝑥)

𝐶 = 15𝐺ℎ2

Substituting to equation of 𝜎𝑦 ,

𝜎𝑦 = 2𝐵 + 30𝐺ℎ2 𝑦 − 10𝐺𝑦 3

𝑞
For 𝜎𝑦 = − 𝑏 and 𝜎𝑦 = 0 at 𝑦 = ℎ and 𝑦 = −ℎ,

𝑞
− = 2𝐵 + 30𝐺ℎ3 − 10𝐺ℎ3 = 2𝐵 + 20𝐺ℎ3
𝑏

And

0 = 2𝐵 − 30𝐺ℎ3 + 10𝐺ℎ3 = 2𝐵 − 20𝐺ℎ3

Solving simultaneously,
 𝑞
𝐵=−
4𝑏

𝑞
𝐺=−
40𝑏ℎ3

Substituting to the equation G,

𝑞 2
3𝑞
𝐶 = 15 (− ) ℎ = −
40𝑏ℎ3 8𝑏ℎ

Substituting E and G to

𝑞 𝑞
𝜎𝑥 = 6𝐷𝑦 + 6(0)𝑥𝑦 − 30(− 3
)𝑥 2 𝑦 + 20(− )𝑦 3
40𝑏ℎ 40𝑏ℎ3

3𝑞 2 𝑞
𝜎𝑥 = 6𝐷𝑦 + 3
𝑥 𝑦− 𝑦3
4𝑏ℎ 2𝑏ℎ3

ℎ 1
For ∫−ℎ 𝜎𝑥 𝑦𝑏𝑑𝑦 = 2 𝑞𝐿2 at 𝑥 = 𝐿,

ℎ
3𝑞 2 𝑞 1
∫ (6𝐷𝑦 + 3
𝐿 𝑦− 3
𝑦 3 ) 𝑦𝑏𝑑𝑦 = 𝑞𝐿2
−ℎ 4𝑏ℎ 2𝑏ℎ 2

ℎ
3𝑞 2 2 𝑞 1 2
∫ (6𝐷𝑦 2 + 𝐿 𝑦 − 𝑦 4
) 𝑏𝑑𝑦 = 𝑞𝐿
−ℎ 4𝑏ℎ3 2𝑏ℎ3 2

𝑞 2 3 𝑞 1 2
(2𝐷𝑦 3 + 𝐿 𝑦 − 𝑦 5
) 𝑏 = 𝑞𝐿
4𝑏ℎ3 10𝑏ℎ3 2

𝑞 2 3 𝑞 1
{2𝐷[ℎ3 − (−ℎ3 )] + 3
𝐿 [ℎ − (−ℎ3 )] − 3
[ℎ5 − (−ℎ5 )]}𝑏 = 𝑞𝐿2
4𝑏ℎ 10𝑏ℎ 2

𝑞 2 3 𝑞 1
{2𝐷[ℎ3 − (−ℎ3 )] + 3
𝐿 [ℎ − (−ℎ3 )] − 3
[ℎ5 − (−ℎ5 )]}𝑏 = 𝑞𝐿2
4𝑏ℎ 10𝑏ℎ 2

𝑞 𝑞ℎ2 1 2
4𝐷𝑏ℎ3 + 𝐿2 − = 𝑞𝐿
2 5 2
 𝑞ℎ2
3
4𝐷𝑏ℎ =
5

𝑞
𝐷=
20𝑏ℎ

Substituting the constants back to the equations for stresses,

𝑞 𝑞 𝑞
𝜎𝑥 = 6 ( ) 𝑦 + 6(0)𝑥𝑦 − 30 (− 3
) 𝑥 2 𝑦 + 20 (− ) 𝑦3
20𝑏ℎ 40𝑏ℎ 40𝑏ℎ3

3𝑞 3𝑞 2 𝑞
𝜎𝑥 = 𝑦+ 3
𝑥 𝑦− 𝑦3
10𝑏ℎ 4𝑏ℎ 2𝑏ℎ3

𝟑𝒒𝒚 𝒒
𝝈𝒙 = (𝟐𝒉𝟐 + 𝟓𝒙𝟐 ) 𝟑
− 𝟑
𝒚𝟑
𝟐𝟎𝒃𝒉 𝟐𝒃𝒉

𝑞 3𝑞 𝑞
𝜎𝑦 = 2 (− ) + 2 (− ) 𝑦 − 10 (− ) 𝑦3
4𝑏 8𝑏ℎ 40𝑏ℎ3

𝑞 3𝑞𝑦 𝑞𝑦 3
𝜎𝑦 = − − +
2𝑏 4𝑏ℎ 4𝑏ℎ3

𝒒
𝝈𝒚 = − (𝟐𝒉𝟑 + 𝟑𝒉𝟐 𝒚 − 𝒚𝟑 )
𝟒𝒃𝒉𝟑

3𝑞 𝑞
𝜏𝑥𝑦 = −(0 + 2 (− ) 𝑥 + 3(0)𝑦 2 − 30 (− 3
) 𝑥𝑦 2 )
8𝑏ℎ 40𝑏ℎ

3𝑞𝑥 3𝑞𝑥𝑦 2
𝜏𝑥𝑦 = −
4𝑏ℎ 4𝑏ℎ3

𝟑𝒒𝒙 𝟐
𝝉𝒙𝒚 = (𝒉 − 𝒚𝟐 )
𝟒𝒃𝒉𝟑

Equations of stresses

𝟑𝒒𝒚 𝒒
𝝈𝒙 = (𝟐𝒉𝟐 + 𝟓𝒙𝟐 ) 𝟑
− 𝒚𝟑
𝟐𝟎𝒃𝒉 𝟐𝒃𝒉𝟑
 𝒒
𝝈𝒚 = − (𝟐𝒉𝟑 + 𝟑𝒉𝟐 𝒚 − 𝒚𝟑 )
𝟒𝒃𝒉𝟑

𝟑𝒒𝒙 𝟐
𝝉𝒙𝒚 = (𝒉 − 𝒚𝟐 )
𝟒𝒃𝒉𝟑

2
Substituting the area moment of inertia is 𝐼𝑧 = 3 𝑏ℎ3 , the equation becomes

𝒒𝒚 𝒒 𝟑
𝝈𝒙 = (𝟐𝒉𝟐 + 𝟓𝒙𝟐 ) − 𝒚
𝟏𝟎𝑰𝒛 𝟑𝑰𝒛

𝒒
𝝈𝒚 = − (𝟐𝒉𝟑 + 𝟑𝒉𝟐 𝒚 − 𝒚𝟑 )
𝟔𝑰𝒛

𝒒𝒙 𝟐
𝝉𝒙𝒚 = (𝒉 − 𝒚𝟐 )
𝟐𝑰𝒛