Massachusetts Institute of Technology Department of Mechanical Engineering Cambridge, MA 02139
Massachusetts Institute of Technology Department of Mechanical Engineering Cambridge, MA 02139
1 Objectives
The primary objectives of this lab are to introduce the concept of stress and strain
concentration factors in notched structural configurations. The notion of stress con
centration is experimentally explored qualitatively, using photoelasticity, and quan
titatively, using experimental, analytical, and numerical methods.
2 Lab Tasks
• review the concepts of linear elastic stress concentration factors associated with
geometric discontinuities such as notches, holes, and cutouts
• generalize the concept of stress concentration factor to finite bodies, and explore
the consequences of St. Venant’s principle for using analytical “infinite body”
stress concentration results to estimate stress concentration associated with
isolated discontinuities
• examine basic output numerical results from a linear elastic plane stress finite
element solution of the loaded notched specimen.
2. For the case of a circular hole in a finitewidth plate, and in particular, for the
geometry of the specimen in Fig 8a, obtain the value of the stress concentration
factor Ktg from the graph in Fig. 3. Compare the level of σθθ /σ at the edge
of the hole (r = a) in an infinite plate, computed in (1) above, to that of the
finitewidth plate. Explain the reasons for any differences between the two stress
concentration factor values.
3. During the lab session, local circumferential strain values, �θθ , were recorded at
various strain gauges placed along the ligament (θ = π/2). Plot these measured
strain values as a function of the radial coordinate r, measured from the center
of the hole.
4. An additional remote strain gauge measured the farfield axial strain. Using
this measurement, estimate the strain concentration factor at the hole, K� (see
eq. (13)). How does K� compare with Ktg obtained in (2) above? Comment on
the similarity/difference between the two values.
5. Using the analytical expressions for σrr , σθθ , and σrθ derived in (1) above, obtain
an analytical expression for the radial distribution of strain, �θθ , on the ligament
at θ = π/2. Use this expression to obtain an approximate prediction of the �θθ
profile along the ligament of the finite aluminum plate in Fig. 8a. The plate
is loaded with an axial force P = 40kN . In your calculations, assume a value
for the Young’s modulus of the aluminum plate as E = 72 GP a, and Poisson
ratio of ν = 0.32. Compare your predictions with the experimental values
obtained in (3) above, and comment on possible sources for differences between
the measurements and the prediction.
3
4 Background on Linear Elastic
Stress Concentration
The stress distributions around a central hole can be estimated for the simple case of
an infinitely wide plate subjected to tensile loading. The overall stress distributions
in the plate are given by (Figure 1)
�� � � � �
σ a2 a4 a2
σrr (r, θ) = 1 − 2 + 1 + 3 4 − 4 2 cos(2θ) (1a)
2 r r r
�� 2
� � 4
� �
σ a a
σθθ (r, θ) = 1 + 2 − 1 + 3 4
cos(2θ) (1b)
2 r r
�� 4 2
� �
σ a a
σrr (r, θ) = − 1 − 3 4 + 2 2 sin(2θ) (1c)
2 r r
∞
where “σ” is the magnitude of the remotelyapplied tensile stress, σxx .
s q s
x
. r
0
0
A B
sqq s
s rr
Tr q
R s rr
A B Tr q C
. C
q
. . .
a
r
.
.
The surface of the circular hole is r = a; the unit normal to this surface is n = −er ,
where er is the unit vector in the local radial direction. We understand that the state
of stress in the plate is approximately plane stress (providing plate thickness h � a),
so that σzz = σzr = σzθ = 0. In this case the traction vector lies in the plane, with
components
t = −(σrr er + σrθ eθ ) = 0, (2)
which is consistent with an evaluation of eq. (1) at r = a:
σrr (r = a, θ) = 0 (3a)
σrθ (r = a, θ) = 0 (3c)
For θ = π/2, the hoop stress in eq. (3b) attains its maximum value of σθθ = σmax =
3 σ. This corresponds to the peak of the stress distribution circumferential stress
distribution shown in Figure 2a. Hence we may say that the stress concentration
factor (the ratio of the maximum local stress [component] to the farfield stress
[component]) for this geometry is equal to 3. The concept of a stress concentration
factor will be further discussed in the following section. However, it is important
to note that stress near the hole greatly exceeds the farfield stress. Consequently,
failure processes 1 may initiate locally at the edge of the hole under values of farfield
stress which are themselves sufficiently small to preclude such failures from occurring
remotely.
Figure 2b, which shows the radial variation of σθθ along the ray θ = π/2, emphasizes
that the magnitude of the stress concentration associated with the hole decays rapidly
with increasing distance from the notch. This is a clear example of St. Venant’s
principle, which states that the perturbations in a linear elastic stress field due to
the presence of an isolated geometrical discontinuity of size “d” are localized within
a region of characteristic linear dimension ∼ 3d from the discontinuity. The stress
levels outside this region are therefore close to the nominal applied stress levels (un
perturbed).
1
e.g., yielding (failure to deform only elastically), fracture (failure to remain an intact body),
fatigue damage (failure to sustain loadcarrying capacity over repeated load cycles) , etc.
+
3s
-s
s . . . . q. s
(A)
s . . s
3s
(B)
Figure 2: Distribution of hoop stress component σθθ : (a) around the circumference
of circular hole in a large body, and (b) radial distribution along the ligament where
θ = π/2.
See also, Crandall, Dahl, & Lardner, Sections 5.9 and 5.7)
Stress concentration factors have been obtained for several geometries of engineering
significance. These are often tabulated in engineering handbooks (e.g., “[R. E.] Peter
son’s Stress Concentration Factors,” 2nd. Ed., John Wiley, New York, 1997, edited
by W. D. Pilkey). Two types of stress concentration factor definitions are found in
the literature, and it is important to apply the analysis consistent with the given
definition. The first type of definition of stress concentration factor is based on the
“gross” applied stress, σg (see equation (6) below for definition of σg ). This stress
concentration factor is given the symbol “Ktg ” according to
σmax
Ktg ≡ , (4)
σg
where σmax is the maximum local stress at the edge of the hole, and σg is the applied
farfield stress remote from the hole.
Similarly, we may also define the stress concentration factor based on the nominal (or
netsection average) applied stress, σnom , according to
σmax
Ktn ≡ . (5)
σnom
Here σnom is the nominal, or netsection average stress acting across the minimum net
area, Anet , of the plane containing the notch. The total tensile load, P , transmitted
along a finitesized member is finite and equal to the product of gross applied stress,
σg , and gross area, Agross ; we also require the product of nominal stress, σnom , and
netsection area, Anet , to equal the load:
so that
Agross
σnom = σg . (7)
Anet
For example, in a uniformthickness plate (thickness = h = constant) of width 2w
containing a central circular hole of radius a,Agross = 2wh and Anet = Agross −2ah; thus
Agross /Anet = 2wh/h(2w−2a) = 1/(1−a/w), so σnom = σg /(1−a/w). Evidently, since
σmax is unique, the stress concentration factors Ktg and Ktn differ by a corresponding
(but inverse!) area ratio:
� �
Anet
Ktn = Ktg = Ktg · (1 − a/w). (8)
Agross
Values for both stress concentration measures Ktg and Ktn for centrallylocated cir
cular holes in tensileloaded plates of finite width are given in the graph of Fig. 3; it
is readily verified from the figure that eq. (8) holds. It is also clear that in the limit
of a finite hole in a very large body, the distinction between “net” and “gross” stress
measures and stress concentration factors vanishes.
For the case of elliptical holes in an infinitely wide plate (e.g., plate width 2w → ∞,
b/w → 0) subject to remote tensile loading, the stress concentration factor Kt∞ can
be obtained analytically, and is given by:
b
Kt∞ = 1 + 2 , (9)
a
where 2a is the length of the ellipse diameter in a direction parallel to the remote
tension, and 2b is the length of the ellipse diameter in the direction perpendicular to
the remote tension. The location of the maximum stress is at the ends of the 2baxis
of the ellipse. Equation (9) is plotted as a solid line in Figure 4.
where ρ is the local radius of curvature of the edge of the ellipse near the 2baxis; for
an ellipse, at this location, ρ = a2 /b.
Equation (10) may be used to estimate the stress concentration factors for other
“nearellipse” isolated hole geometries in wide plates by introducing the idea of an
“equivalent” ellipse. Figure 5 illustrates the concept for three different nonelliptical
notch shapes of actual tip root radius ρ = r and actual length 2b projected perpen
dicular to the remote tensile axis; the “equivalent ellipse” for each notch is shown in
Use of the equivalent ellipse concept then allows an estimation of the maximum stress
at the root of the isolated nearelliptical notch, for remote applied stress “σ ∞ ,” as
� �
. ∞ ∞ b
σmax = σ × Kt∞ (b, aequivalent ) = σ × 1 + 2 . (11)
aequivalent
Image removed due to copyright considerations. See Chart 4.50 in Pilkey, Walter. Peterson's
Stress Concentrations, 2nd ed. New York: Wiley, 1992.
Image removed due to copyright considerations. See Chart 4.1 in Pilkey, Walter. Peterson's
Stress Concentrations, 2nd ed. New York: Wiley, 1992.
10
s
2b
. . .
r r 2a
. . .
s
Slit
Equivalent
Ellipse r
r 2a
. .
t t
2b
12
6 Introduction to Photoelasticity
Stress distributions in some transparent materials may also be measured using pho
toelastic techniques. These rely on illumination with plane polarized light obtained
by passing light rays through vertical slots that produce polarized light beams hav
ing waves that oscillate only along one plane perpendicular to the travel of the light
ray. When the polarized light passes through a material which is being stressed in
a direction parallel to the polarizing axis, then a fringe pattern is formed against a
light (bright field) background. Conversely, when the stress axis is perpendicular to
the polarizing axis, a “dark field”, or black image is formed.
In some materials, the application of stress may cause an incident planepolarized ray
to split into two coincident rays with directions that coincide with the directions of the
inplane principal stress directions. Since this phenomenon is only observed during
the temporary application of stress, it is known as “temporary birefringence”. Fur
thermore, the speeds of the rays are proportional to the magnitudes of the respective
principal stresses. Hence, the emerging rays are out of phase. They therefore produce
interference fringe patterns when they are recombined. If they are recombined at an
analyzer (shown in Fig. 6), then the amount of interference in the emerging rays is di
rectly proportional to the difference between the local inplane principal stress levels,
σp and σq . therefore, the amount of interference is related to the maximum inplane
shear stress, which is given by
1
τmax = (σp − σq ). (12)
2
The fringe patterns therefore provide a visual indication of the spatial variations in
the maximum inplane shear stress.
13
Model
Load
Polarised
light axis
(a) Unloaded-no birefringence
(polarised light passes through unaffected)
sp s p Plane
sq A sq Plane
Plane polarised
light
Figure e: Temporary birefringence. (a) Plane-polarised light directed onto an unstressed model passes
through unaltered. (b) When the model is stressed the incident plane-polarised light is split
into two component rays. The directions of the rays coincide with the directions of the principal
stresses, and the speeds of the rays are proportional to the magnitudes of the respective stresses in their
directions. The emerging rays are out of phase, and produce an interference pattern of fringes.
Two experiments will be performed in this lab. In the first experiment, an aluminum
plate of thickness h and width 2b, with a central circular hole of diameter 2a (see
Figure 8a), will be subjected to elasticlevel tensile load of magnitude P . Strain
distributions will be measured using strain gauges attached to different positions on
the plate. Local tangent strains, �θθ , will be measured at five positions (gauges #1,
2, 3, 5, and 8, as detailed in Figure 9) along the ligament. An additional, remotely
located strain gauge (#10 in Fig. 9), oriented parallel to the loading direction, will be
used to obtain measurements of axial farfield strain. Exact positions of the gauges
with respect to the hole are given in Fig. 9.
6" 6"
strain gages
R=0.5"
1.5"
24" 24"
1.5"
0.25"
(a) (b)
Figure 8: Geometry of the plates tested in this module. (a) Plate with a circular hole:
both an aluminum and a polycarbonate plate of this geometry will be tested. In the
figure, the positions of strain gauges applied to the aluminum specimen are schemat
ically indicated. Precise gauge locations are shown in Fig. 9. (b) Polycarbonate plate
with a slot.
15
10
11
37.8 mm
1 2 3 4 5 6 7 8 9
2.3 mm
6.5 mm
11.0 mm
14.7 mm
18.7 mm
22.6 mm
26.0 mm
30.8 mm
16
Today, most quantitative stress analysis is conducted using powerful computer tech
niques; the most widelyused of these is the finite element method. In applications
to isotropic linear elasticity, the method constructs an approximate (but [with suffi
cient care!] numerically accurate) solution to the equations governing elastic boundary
value problems. In practice, the approximate displacement field is constructed as a
piecewisecontinuous, loworder polynomial within small (but finite) “elements;” in
twodimensional models (e.g., plane stress), the elements are typically triangular or
quadrilateralshaped. Figure 10 shows a mesh of quadrilateral elements used here to
model the holeinplate tension specimen.
As part of the output of the approximate solution, we can obtain color contour plots of
the spatial variation of particular stress components. Figure 11 shows a detail of the
variation of stress component σ22 when the farfield load is P = 10 kN , leading to a
[remote] gross stress of σgross = P/2bh = 10.33 M P a (plate thickness is h = 6.35 mm;
total plate width is 2b = 152.4 mm; hole diameter is 2a = 38.1 mm); Figure 12 shows
a similar variation of shear stress component σ12 .
The contour plots can also be “sectioned” along a given path within the material, to
plot the variation of a particular variable along that path. Figure 13 plots the “radial”
variation (variation in direction x1 ) of normalized stress component σ22 across the
minimum ligament at the horizontal centerplane of the specimen. The local stress
17
3 1
Figure 10: Finite element mesh used to generate approximate stress analysis solution.
is normalized by the remote gross stress, σg = P/2bh. The peak value at the root of
the hole is 3.26, dropping to 0.904 at the edge of the specimen. Figure 14 plots the
normalized strain component �22 over the same region; here strain is normalized with
remote gross strain, �g = σg /E.
18
S22 VALUE
-9.24E+07
+2.22E-16
+3.18E+06
+6.36E+06
+9.55E+06
+1.27E+07
+1.59E+07
+1.91E+07
+2.23E+07
+2.55E+07
+2.86E+07
+3.18E+07
+3.50E+07
+2.04E+08
Figure 11: Contour plots of axial stress component σ22 . Remote (gross) stress level
∞
is σ22 = P/(2bh) = 10.33 M P a.
S12 VALUE
-1.06E+08
-1.00E+07
-8.18E+06
-6.36E+06
-4.55E+06
-2.73E+06
-9.09E+05
+9.09E+05
+2.73E+06
+4.55E+06
+6.36E+06
+8.18E+06
+1.00E+07
+1.06E+08
Figure 12: Contour plots of shear stress component σ12 . Remote (gross) stress level
∞
is σ22 = P/(2bh) ≡ σg = 10.33 M P a.
19
3.2
LIGSTRES
2.8
2.4
1.6
1.2
0.8
3.2
ELIGN
2.8
2.4
normalized ligament strain
2.0
1.6
1.2
0.8
20