Stress Application Example
Stress Application Example
Stress Analysis
As you have learned from CVE 220 and/or MCE 301, when an elastic body is
subjected to applied loadings, stresses are created inside the body. In general
these stresses often vary in complicated ways from point to point and from plane
to plane within the structure. To help characterize this situation, stresses are
normally defined with respect to a given coordinate system. The illustration
below shows that at a typical point within a loaded body, the state of stress can
be characterized on a small cube of material defined with respect to a Cartesian
coordinate system. These nine components are called the stress components,
with x, y, z referred to as normal stresses and xy, yx, yz, zy, zx, xz called the
shearing stresses.
P3 y
P2
yz yx
xy
zy
y x
zx xz
z
p
x
P1
(Externally Loaded Body) z (State of Stress at a Point)
x xy xz
[ ] yx y yz (1)
zx zy z
and since the shearing stresses have the equalities xy = yx , yz =zy , zx = xz , the
stress matrix is symmetric.
Again from your previous courses, you studied how the stresses vary by
orientation of the plane of action. For example, as we change the orientation of a
particular plane the normal stress component (say x) will vary. There exists a
special orientation where the normal stress will be a maximum, and these are
called principal planes and the normal stresses acting on them are called the
principal stresses. It turns out that for the general three-dimensional case, the
theory to determine principal stresses and the planes on which they act is
formulated by the eigenvalue problem
where [] is the stress matrix, {n} is the principal direction vector and (the
eigenvalue) is the principal stress. Thus solving the eigenvalue problem will
determine up to three distinct principal stresses and the corresponding three
principal directions. It turns out for this application (3 x3, symmetric real matrix)
the principal directions are mutually orthogonal. Again you should remember
that the shear stress components will vanish on these three principal planes and
so for a coordinate system that is aligned with the principal directions the stress
matrix takes on the simplified diagonal form
1 0 0
[ ] 0 2 0 (3)
0 0 3
y
y 2
2
1
yz yx
xy 1
zy
x
zx xz
3
z x
z 3
(General Coordinate System) (Principal Coordinate System)
y
P
x
2c
2L
From mechanics of materials theory the in-plane stress components for the right
half of the beam (x > 0) are given by
3P 3P 2
x 3
( L x ) y , y 0 , xy 3
(c y 2 ) (4)
4c 8c
Thus the stress matrix for this problem reduces to the 2x2 form
3 3P
x ( L x ) y ( c 2
y 2
)
xy 2c 3
8c 3
[] 3P
xy y 3 (c y )
2 2
0
(5)
8c
We can now employ MATLAB to calculate the principal stresses and make a
contour plots of the distribution in the region 0 < x < L, -c < y < c. Contour plots
are shown below along with the MATLAB code.
First Principal Stress Contours, 1
0.4
0.2
-0.2
-0.4
0.4
0.3
0.2
0.1
0
-0.1
-0.2
-0.3
-0.4