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
𝒒𝒚 𝒒 𝟑
𝝈𝒙 = (𝟐𝒉𝟐 + 𝟓𝒙𝟐 ) − 𝒚
𝟏𝟎𝑰𝒛 𝟑𝑰𝒛
𝒒
𝝈𝒚 = − (𝟐𝒉𝟑 + 𝟑𝒉𝟐 𝒚 − 𝒚𝟑 )
𝟔𝑰𝒛
𝒒𝒙 𝟐
𝝉𝒙𝒚 = (𝒉 − 𝒚𝟐 )
𝟐𝑰𝒛
