Chapter 5
Analytical
techniques and solutions for linear elastic solids
In
the preceding chapter, we solved some simple linear elastic boundary value
problems. The problems were trivial,
however, in that the stress, strain and displacement fields were axially or
spherically symmetric. Most problems
of practical interest involve fully 3D stress and displamement fields.
It
is extremely difficult to solve general boundary value problems. However, some of the best mathematicians
over the past 200 years have turned their attention to this matter, and have
developed several very elegant techniques.
None of these are completely general, but solutions derived using
these techniques have provided invaluable insight into the behavior of
deformable solids.
Sadly,
these days any fool with a PC and a finite element package can solve
virtually any linear elastic boundary value problem, so you will not be able
to make a living calculating exact elasticity solutions. Nevertheless, some exact solutions are of
fundamental practical importance. Examples include contact problems,
solutions for cracks, stress concentrations, thermal stress problems, and
problems involving defects such as dislocations in solids. It is worth seeing how such solutions were
derived.
In
addition, there are some very powerful theorems in linear elasticity (such as
the principle of minimum potential energy, and the reciprocal theorem), which
can be used to directly calculate quantities of interest without necessarily
having to solve all the governing equations directly.
In
this chapter, we present a very brief survey of the field of linear
elasticity. Specifically,
1. We will outline some important general features of
solutions to boundary and initial value problems;
2. We will discuss some solution techniques and present
solutions to selected boundary value problems of interest;
3. We will discuss energy methods for solving problems
involving linear elastic solids, including the principle of minimum potential
energy, the reciprocal theorem, and the Rayleigh-Ritz method for estimating
natural frequencies of vibrating elastic solids.
5.1 General Principles
5.1.1
Summary of the governing equations of linear elasticity
Static problems. We already listed the governing
equations of linear elasticity in our discussion of simple axisymmetric
problems. They are repeated here for convenience.
Given:
1. The shape of the solid in its unloaded condition
2. The initial stress field in the solid (we will take
this to be zero)
3. The elastic constants for the solid and its mass density
4. The thermal expansion coefficients for the solid,
and temperature change from the initial configuration
5. A body force distribution (per unit mass) acting on the solid
6. Boundary conditions, specifying displacements on a portion or tractions on a portion of the boundary of R
Calculate
displacements, strains and stresses satisfying the governing equations of
linear elastostatics
Dynamic problems Dynamic
problems are essentially identical, except that the boundary conditions must
be specified as functions of time, and the initial displacement and velocity
field must be specified. In this case
the governing equations are
5.1.2 Alternative form of the
governing equations the Navier equation
The
governing equations can be simplified by eliminating stress and strain from
the governing equations, and solving directly for the displacements. In this case the linear momentum balance
equation (in terms of displacement) reduces to
For
the special case of an isotropic solid with shear modulus and Poisson ratio and uniform temperature this equation reduces to
These are known as the
Navier (or Cauchy-Navier) equations of elasticity.
The boundary conditions
remain as given in the preceding section.
5.1.3
Superposition and linearity of solutions
The
governing equations of elasticity are linear. This has two important consequences:
1. The stresses, strains and displacements in a solid
are directly proportional to the loads (or displacements) applied to the
solid.
2. If you can find two sets of displacements, strains
and stresses that satisfy the governing equations, you can add them to create
more solutions.
These
principles can be illustrated clearly using some of the simple solutions
derived in Section 4.1. For example,
examine the solution to the pressurized sphere (Sect 4.1.4). As an example, the radial stress induced
by pressure on the interior, and zero pressure on the
exterior surface is
The radial stress induced by pressure on the exterior surface, with zero pressure
on the interior surface is
Note
that in both cases the stress is directly proportional to the pressure. In addition, to find the radial stress by
combined pressures on the interior and on the exterior surface, you can just add
these two solutions.
Further examples of superposition and linearity will be
given in subsequent sections.
5.1.4
Uniqueness and existence of solutions to the linear elasticity equations
The following results are
useful:
1. If only displacements are prescribed on the boundary
of the solid, the governing equations of linear elasticity always have a
solution, and the solution is unique.
2. If mixed boundary conditions are specified, a static
solution exists and is unique if the displacements constrain rigid motions. A
dynamic solution always exists and is unique, provided the velocity field and
displacement field at time t=0
are known.
3. If only tractions are prescribed on the boundary, a
static solution exists only if the tractions are in equilibrium. In this case, the stresses and strains are
unique, but the displacements are not.
A dynamic solution always exists and is unique, again, providing
initial conditions are known.
5.1.5 Saint-Venant’s principle
Saint-Venant’s
principle is often invoked to justify approximate solutions to boundary value
problems in linear elasticity. For
example, when we solve problems involving bending or axial deformation of slender
beams and rods in elementary strength of materials courses, we only specify
the resultant forces acting on the ends of a rod, or the magnitudes of point
forces acting on a beam, we don’t specify the distribution of traction in
detail. We rely on Saint Venant’s
principle to justify this approach. In this context, the principle states the
following.
The
stresses, strains and displacements far from the ends of a rod or beam
subjected to end loading depend only on the resultant forces and moments acting
on its ends, and do not depend on how the tractions themselves are
distributed.
Although
SVP is widely used, it turns out to be remarkably difficult to prove
mathematically. The difficulty is
partly that it is not easy to state the principle itself precisely enough to
apply any mathematical machinery to it.
A rigorous statement is given by Sternberg (Q. J. Appl. Mech 11 p. 393 1954), among several other
versions. Here, we will just
illustrate the most common applications of the principle through specific
examples.
One
version of SVP can be stated as follows.
Suppose
that we calculate the stress, strain and displacement induced in a solid by
two different traction distributions and that act on some small region of a solid
with characteristic size a. If the
tractions exert the same resultant force and moment, then the stresses,
strains and displacements induced by the two traction distributions at a
distance r from the loaded region are identical for large .
In practice `large’ usually means .
This
principle can be illustrated using a simple example.
Consider a large solid with a flat surface, as shown
in the picture. It is possible to
calculate formulas for the stresses and displacements induced by various
pressure distributions acting on the flat surface the procedure to do this will be outlined
later. For now, we will compare the
stresses induced by
1. A uniform pressure
2. A parabolic pressure
You can verify for yourself that both pressure
distributions exert a resultant force P
acting in the vertical direction on the surface, and exert zero moment about
the origin. The variation of stress
down the axis of symmetry ( ), expressed in cylindrical-polar
coordinates, can be derived as
Case 1: Uniform pressure
Case 2: Parabolic pressure
Now, to demonstrate SVP, we
want to show that the stresses are equal for large z/a. We can do this
graphically the figures below compare the variation of
vertical and radial stress down the axis of symmetry with z/a.
The
stresses induced by the two different pressures are clearly indistinguishable
for . This example helps to quantify what we mean
by a `large’ distance.
The
second commonly used application of SVP is a rather vague statement that
A localized geometrical feature with
characteristic size R in a large solid only influences the stress in a region
with size approximately 3R surrounding the feature.
This
is more a rule of thumb than a precise mathematical statement. It can be illustrated by looking at
specific solutions. For example, the
figure below shows the Mises stress contours surrounding a circular hole in a
thin rectangular plate that is subjected to extensional loading (calculated
using the finite element method). Far from the hole, the stress is
uniform. The contours deviate from the
uniform solution in a region that is about three times the hole radius.
|