ISSN 2303-4521

PERIODICALS OF ENGINEERING AND NATURAL SCIENCES

Vol. 4 No. 1 (2016)

Available online at: http://pen.ius.edu.ba

Validation of Stresses with Numerical Method and Analytical

Method

Enes Akca

International University of Sarajevo, Department of Mechanical Engineering

Sarajevo, Bosnia and Herzegovina

[email protected], [email protected]

ABSTRACT
In this study, the goal is to confirm stresses under kind of loads, with analytical and
numerical methods. But each material is drawn in SolidWorks 3D design software, and
exported to ANSYS. Also the other purpose is to give information about type of stresses
and formulation of each stress with different loads.

Keywords:ANSYS, Stress, Shear stress, FEM, Analyticalsolution, Numerical solution

1. INTRODUCTION

Stress is defined as the strength of a material per unit area or unit strength. It is the force on a
member divided by area, which carries the force, formerly express in N/mm2or MPa. And there
are types of stresses and component of stresses; normal stress, shear or transversal stress,
maximum and minimum principle stresses. Besides, they will be defined in detail in chapter 2.
As computers have become more and more powerful, people have tended to use numerical
approaches to develop theoretical models to predict the effect of whatever is studied. This has
improved stress analysis and computer simulations. Numerical methods can potentially provide
more accurate solutions since they normally require much less restrictive assumptions. The finite
element method is very often used to analyze the stress state of an elastic body with complicated
geometries.
The finite element method is capable of providing this information, but the time needed to create
such a model is large. In order to reduce the modeling time, a 3D model created in solid

DOI: 10.21533/pen.v4i1.47 17
modeling software can be used. One such model is provided by SolidWorks. SolidWorks can
generate models of three-dimensional members easily. In SolidWorks, the geometry is saved as a
file and then it can be transferred from SolidWorks to ANSYS. ANSYS has pioneered the
development and application of simulation methods to solve the most challenging product
engineering problems. Simulation software enables organizations to confidently predict how
their products will operate in the real world [1].
First of all, it should be emphasized that the "numerical approach" is not automatically
equivalent to the "approach with use of computer", although we usually use numerical approach
to find the solution with use of computers. That is because of the high computer performance
incomparable to abilities of human brain. Numerical approach enables solution of a complex
problem with a great number of very simple operations. It can be distinguished two main
situations when numerical methods are used instead of analytical methods [2];
1. When analytical solution of the mathematically defined problem is possible but it is time-
consuming and the error of approximation we obtain with numerical solution is acceptable. In
this case the calculations are mostly made with use of computer because otherwise its highly
doubtful if any time is saved. It is also indivually to decide what we mean by "time-consuming
analytical solution". In my discipline even very simple mechanical problems are solved
numerically simply.
2. When analytical solution is impossible, this means that we have to apply numerical
methods in order to find the solution. This does not define that we must do calculations with
computer although it usually happens so because of the number of required operations.
Hence in the next chapters, there will be results of both numerical and analytical solutions, also
confirmation of 3 tests and one member is just consisting of numerical solution. On account of
the fact that it will be proved that ANYS results are reliable if 3 tests are confirmed.
2. STRESSES

The purpose is to find the normal and shear stresses acting on any inclined section. For uniaxial
load and pure shear, now the transformation relationships that give the stress components for any
orientation is need to be derived. This is referred as stress transformation[3].
When an element is rotated from one orientation to another, the stresses acting on the faces of
the element are different but they still represent the same state of stress, namely, the stress at the
point under consideration.

Axial load 𝜎 = 𝑃/𝐴

Torsional load in circular shaft 𝜏 = 𝑇𝑝 /𝐼𝑝

Bending moment and shear force in beam 𝜎 = 𝑀𝑦 /𝐼 𝜏 = 𝑉𝑄/𝐼𝑏

Plane stress occurs when the material at a point is subjected to two normal stress components
σxandσy and a shear stress τxy. Provided these components are known, then the stress components
acting on an element having a different orientation θ can be determined using the two force
equations of equilibrium or the equations of stress transformation [4].

𝜎𝑥 + 𝜎𝑦 𝜎𝑥 − 𝜎𝑦
𝜎𝑥 ′ = + 𝑐𝑜𝑠2𝜃 + 𝜏𝑥𝑦 𝑠𝑖𝑛2𝜃
2 2

18
 𝜎𝑥 − 𝜎𝑦
𝜏𝑥 ′ 𝑦 ′ = − 𝑠𝑖𝑛2𝜃 + 𝜏𝑥𝑦 𝑐𝑜𝑠2𝜃
2

Figure 1: Plane stress

It is important to determine the orientation of the element that produces the maximum principal
normal stresses and the maximum in-plane shear stress. Using the stress transformation
equations, it is found that no shear stress acts on the planes of principal stress.The principal
stresses are

𝜎𝑥 + 𝜎𝑦 𝜎𝑥 − 𝜎𝑦 2
𝜎1,2 = ± √( ) + 𝜏𝑥𝑦 2
2 2

The planes of maximum in-plane shear stress are oriented 45° from this orientation, and on these
shear planes there is an associated average normal stress.

𝜎𝑥 − 𝜎𝑦 2
𝜏 𝑚𝑎𝑥 = √( ) + 𝜏𝑥𝑦 2
𝑖𝑛−𝑝𝑙𝑎𝑛𝑒 2

𝜎𝑥 + 𝜎𝑦
𝜎𝑎𝑣𝑔 =
2

Figure 2: Planes of maximum in-plane shear stress

Mohr’s circle provides a semi-graphical method for finding the stress on any plane, the principal
normal stresses, and the maximum in-plane shear stress. To draw the circle, the σ and τ axes are
established, the center of the circle C [(σx+σy)/2, 0]and the reference point A (σx,τxy)are plotted.

19
The radius R of the circle extends between these two points and is determined from trigonometry
[4].

Figure 3: Mohr’s circle explanation

If σ1 and σ2 are of the same sign, then the absolute maximum shear stress will lie out of plane.
𝜎1
𝜏 𝑎𝑏𝑠 =
𝑚𝑎𝑥 2
In the case of plane stress, the absolute maximum shear stress will be equal to the maximum in-
plane shear stress provided the principal stresses σ1 and σ2 have the opposite sign [4].
𝜎1 − 𝜎2
𝜏 𝑎𝑏𝑠 =
𝑚𝑎𝑥 2

Figure 4: Absolute maximum shear stress

3. APPLICATION OF STRESSES
3.1.Tensile Test
Tensile testingis a fundamentalmaterials sciencetest in which a sample is subjected to a
controlledtensionuntil failure. The results from the test are commonly used to select a material
for an application, forquality control, and to predict how a material will react under other types
offorces. Properties that are directly measured via a tensile test areultimate tensile strength,
maximumelongationand reduction in area.From these measurements the following properties can
also be determined:Young's modulus,Poisson's ratio,yield strength, andstrain-
hardeningcharacteristics. And tension or tensile force will tend to lengthen the member [5].

20
 Figure 5: Drawing the bar
A bar that is applied the forces on the both edges of it.

Figure 6: 3D drawing of the bar

Figure 7: Stress distribution of the bar

Properties and dimensions of the member;
And numerical (ANSYS)result;

21
 Material Alloy steel
Tensile Yield Strength σy 250 MPa
Poisson's Ratio 0,3
Young's Modulus 210 GPa
Diameter 10 mm
Fn 10 kN
Normal Stress 127,46 MPa

Result from analytical calculation;

𝐹
𝜎=
𝐴
10𝑥103
𝜎= 𝜋∗(10∗10−3 )2
4

σ = 127 MPa
Hence, the results are close enough to confirm the analytical solution and numerical solution. It
is shown in the table below;

Analytical Result 127 Mpa

Numerical (ANSYS) Result 127,46 MPa

3.2.Pressure Test
The pressure is applied inside of a pipe, and it is seen that maximum stress will occur inside of
the pipe.

Figure 8: 3D model of the pipe

22
 .
Figure 9: Drawing of the pipe (dimensions are in mm)

Figure 10: Stress distribution of the pipe

Properties and dimensions of the member;
And numerical (ANSYS)result;

Material Alloy steel

Tensile Yield Strength σy 250 MPa
Poisson's Ratio 0,3
Young's Modulus 210 GPa
Out diameter 200 mm
Inner diameter 188 mm
Thickness 6 mm
Pressure 5 Mpa
Maximum Principal Stress 81.156 MPa

Result from analyticalcalculation;

OD = 200 mm
ID = 200-12= 188 mm
188
Ri = = 94 mm
2

23
 Ro = 100 mm
200−188
t= = 6 mm
2

𝑃 𝑥 𝑅𝑖
σt = 𝑡
≤σt,max
5 𝑥 94
σt = 6
≤σt,max

σt,max = 80 MPa
Hence, the results are close enough to confirm the analytical solution and numerical solution
(ANSYS). It is shown in the table below;

Analytical result 80 MPa

Numerical (ANSYS) result 81.156 MPa

3.3.Combined Load
Force A, B and C are applied to the surface of the L and on the other surface of the L is fixed
support with D that is shown in figure 11. Hence the purpose is to calculate the stresses on
different forces [6].

Figure 11: Forces on the L

24
 Figure 12: Dimension and the point of L (dimensions are in mm)

Properties of the member;

Material Alloy steel

Tensile Yield Strength σy 250 MPa
Poisson's Ratio 0,3
Young's Modulus 210 GPa
Diameter 20 mm
Force A 355 N
Force B 333 N
Force C 444 N

All the calculations are made by manually in below, and each of it represents types of moment
and type of forces that includes axial forces, transversal forces in 3D, torsion and bending
moment in 2D.
Fx =333 N (axial tension)
Fy =355 N (transverse shear in y direction)
Fz =- 444 N (transverse shear in –z direction)
Mx=- (7x10-2) x 355 = -24,85 Nm (torsion x-axis)
My=333 x 7 x 10-2- 444 x 20 x 10-2= -65,5 Nm (bending moment y-axis)
Mz= - 355 x 20 x 10-2 = -71 Nm (bending moment z-axis)
𝜋 𝑥 𝑑2
𝐴= = 3,14 𝑥 10−4 m2 (surface area)
4

𝜋 𝑥 𝑟4
𝐼= = 7,854 𝑥 10−9 m4 (moment of inertia)
4

𝜋 𝑥 𝑟4
𝐽= = 1,57 𝑥 10−8 m4 (polar moment of inertia)
2

𝐹 333
𝜎=𝐴= = 1,06 𝑀𝑃𝑎 (Axial)
3,14 𝑥 10−4

𝑀𝑥 𝑟 24,85 𝑥 0,01
𝜏𝑡𝑜𝑟 = = = 15,82 𝑀𝑃𝑎 (Torque)
𝐽 1,57 𝑥 10−8

4𝑉 4 𝑥 444
𝜏𝑥𝑦 = = = 1,88 𝑀𝑃𝑎 (Transversal shear)
3𝐴 3,14 𝑥 3 𝑥 10−4

𝑀𝑧 𝑦 71 𝑥 0,01
𝜎= = = 90,4 𝑀𝑃𝑎 (Bending)
𝐼 7,854 𝑥 10−9

Thus, with those calculation, it will be calculated normal stress, shear stress, maximum principle
stress, minimum principle stress in the next subtitles, and also the calculation of those stresses is
explained in 2nd chapter as theory.

25
3.3.1. Normal Stress

Figure 13: Normal stress

σx= 1,06 + 90,4 = 91,46 MPa

τ= 0
A stress is defined as the load divided by the area. For a normal stress, it is all loads
perpendicular to the surface.

3.3.2. Shear Stress

Figure 14: Shear stress

σ= 91,46 MPa
τ= τxy +τtor = 1,88 + 15,82 = 17,7 MPa
Shear stress is a stress state where the stress is parallel to the surface of the material, as opposed
to normal stress when the stress is vertical to the surface.

26
3.3.3. Maximum Principle Stress

Figure 15: Max principle stress

𝑅 = √(91,46⁄2)2 + 17,7 = 45,92

91,46
σ1 = 2
+ 45,92 = 91,65 𝑀𝑃𝑎

3.3.4. Minimum Principle Stress

Figure 16: Min principle stress

𝑅 = √(91,46⁄2)2 + 17,7 = 45,92

91,46
σ2 = 2
− 45,92 = − 0,19 𝑀𝑃𝑎

3.3.5. Results

Analytical (manual) Results

Normal Stress 91,46 Mpa
Shear Stress 17,7 MPa
Max Principle Stress 91,65 MPa

27
 Min Principle Stress -0,19 MPa
Numerical (Ansys) Results
Normal Stress 91,347 Mpa
Shear Stress 13 <τ< 20 MPa
Max Principle Stress 76 <σ1< 99 MPa
Min Principle Stress -14<σ2< 5 MPa

Hence, results which including analytical and numerical are compared, and it is confirmed in the
tables above.
The mesh size must depend on the element size and it is an important topic in a finite element
method because of its relationship to accuracy. As the element size increases, it is seen that the
results become more accurate and almost constant in table above [7].

4. CONCLUSION
Finite element modelling of the members in some forces was examined in detail in chapter 3, and
it is confirmed that the manual solutions and analytical solutions are verified. Hence it is proven
that the software (ANSYS) is reliable,and it is considered that ANSYS results are true since the
results approach each other and it is nearly constant. Also, in the chapter 4, it is tested only with
numerical solution, it is seen that results are nearly constant and more accurate when element
sizing is 40, 50 and 60 for both maximum principle stresses and minimum principle stresses.

REFERENCES
1. M.J.Fagan. Finite Element Analysis-Theory and Practices, Longman Group UK Ltd, 1992.
2. W. O. Jolley, J .J. Rencis, and H. T. Grandin, Jr., “A Module for Teaching Fundamentals of
Finite Element , Theory and Practice Using Elementary Mechanics of Materials”,
Proceedings of the 2003 American , Society for Engineering Education Annual Conference
& Exposition, Nashville, Tennessee, June 2003.
3. B.J. Hamrock, B. Jacobson, S. R. Schmid, “Fundamentals of Machine Elements”, 1999.
4. R. C. Hibbeler, “Mechanics of Materials”, 8th edition in SI units, Prentice-Hall, 2011.
5. Budynas, R. G., & Nisbett, J. K., “Shigley's Mechanical Engineering Design”, New York:
McGraw-Hill, 2008.
6. V.B.Bhandari, “Design of Machine Elements”, Tata McGraw Hill, New Delhi, 2010.
7. R. L. Norton, “Machine Design: An Integrated Approach”, 2nd edition, Prentice-Hall,New
Jersey, 2000.

28