Filonenko-Borodich - Theory of Elasticty
Filonenko-Borodich - Theory of Elasticty
Filonenko-Borodich - Theory of Elasticty
FILONENKO-BORODICH
TEXT FLY WITHIN
THE BOOK ONLY
<f)> 00 ^
OU 164095
PEACE
PUBLISHERS
M. M. 4>HJIOHEHKO BOPOHH4
ynpyrociH
POCy^APCTBEHHOE H3#ATEJIbCTBO
<I>H3MKO-MATEMATMMECKOft JIMTEPATypbl
MocKea
of
ELASTICITY
by
M. FILONENKO-BORODICH
by
MARINA KONYAEVA
Introduction 9
16. General 72
17. Stfaihs' Expressed in Terms of Stresses 75
18. Stresses Expressed in Terms of Strains . .78
*19. Work Done by Elastic Forces in a Solid . 82
*20. Potential of Elastic Forces, . , 83
CONTENTS
*
Chapter IX. More General Methods of Solving Elasticity Problems ... 265
*
Chapter XI. Variational Methods of the Theory of Elasticity '. . 350
77. Variational Principles of the Theory of Elasticity. Fundamen-
tal Integral Identity 350
78. Lagrange's Variational Equation 352
79. Ritz-Timoshenko Method 358
80. Castigliano's Variational Equation 364
81. Application of
Castigliano's Variational Equation to the
Problem of Torsion of a Prismatical Rod 368
82. First Problem of the Theory of Elasticity; Second Theorem
of Minimum Energy 373
83. Approximate Method Based on Variational Equation (11.61) . .375
84. Lame's Problem for an Elastic Rectangular Prism 379
References 387
not involve internal forces since they are pairwise mutually equi-
librated on the basis of Newton's third law of equality of action
and counteraction. The theory of elasticity sets forth the problem
of determining internal forces in a solid body. These forces re-
present interaction between molecules; they ensure the existence of
a solid body as such, its strength. They also act when no external
forces are applied to the body; these forces are not in themselves
the object of study in the theory of elasticity; under the action of
external forces the body deforms, the mutual position of molecules
changes and so do the distances between
them; the action of the external forces that
produce deformation gives rise to addi-
tional internal forces. Statics deals with
external forces applied to a body; it is
necessary therefore first to convert the
internal forces in question into external
Forces. This is achieved by the so-called
method of sections which consists in the
following (Fig. 1).
Let us imagine a body cut into two
l
Darts by a plane this plane intersects
\
1
A curved surface is often taken instead of a plane.
12 INTRODUCTION
lim -r- =P
is called stress of internal forces or simply the stress of the point
M of the body on the area lying in the plane of the section made.
The notion of stress in the theory of elasticity is a fundamental
measure of the intensity of internal forces brought about by the
method of sections. The stress is characterised by its magnitude
2
having the dimension force/length 1, direction and position, i.e.,
|
The passage, made above, from the average stress across the
area to the stress at the point is connected with the imaginary
7
process of reducing the dimensions of the area A/ to zero, which
is necessary for the application of the analysis of infinitesimals.
The validity and justification of such a formal process, as already
stated above, had long been questioned and discussed by scientists;
however, the application of the resulting basic equations of the
theory of elasticity to the solution of problems of physics demon-
strated fairly soon the efficiency of the methods worked out and
brought a series of remarkable results confirmed by experience;
this pertains first and foremost to the study of vibrations and
propagation of waves (e.g., sound waves) in elastic bodies; some
of the simpler problems of this kind are treated in Chapters IV
and IX. The middle of the 19th century was especially rich in
achievements as regards the development of the theory of elasticity
and the solution of problems of importance for physics and
technical mechanics; the work of the eminent French investigator
Saint-Venant and his pupils played a dominant part here. Under
these circumstances, doubts concerning the physical justification
of the method of elasticity dealing, as it were, with a continuous
solid medium gradually disappeared. From this point of view, it
is sometimes said that the theory of elasticity is based on the
Theory of Stress
Fig. 3
Fig. 4
the part A after the section has been made, the outward normal to
the same area dF will be directed opposite; we denote it, therefore,
by v. The projections of the total stress on it (Fig. 3) have to
be denoted as: X. Vl Y~ v Z- v it is obvious that
, \
y -v
s+ YV
'* Y -v Y V* 7
** v
7
z"z/
Vy y*
Fig. 5
Xx ,
Y x Zx
, \ X y,
Yr Zy \ X
z, Y z Zz
,
.
AY*> YV yi ^z
7
in whose notation both letters are the same, are normal stresses.
The remaining six stresses are
shearing stresses.
In order to gain a better orienta-
tion in these new notations the reader
is recommended to do the following
exercises.
Exercises.
<
?. A cylindrical body(Fig. 7) is twisted by forces applied to the end cross
sections. Designate the stresses at any point k of the cross section In and
write down the condition that the normal stress is absent. Designate the stress
at any point m
of the lateral surface and write
down the conditions trjat this surface is not
loaded.
3. A
cylindrical rod (Fig. 8) is bent by
applied to the end cross section pq.
forces
n Write down the conditions that the lateral sur-
face is free of load (see the preceding exercise),
designate those stresses at the section mn (nor-
mal and shearing) which are determined in
strength of materials.
ig. 7 Fig. 8
so that
= T^y > *
2
==T y2 ^Z == ^XZ"
Fig. 9
z)\
(1-1)
z).
2-1013
18 THEORY OF STRESS [Ch f
= 0, (1.2)
(1.4)
as indicated in parentheses.
After cancelling out like terms in (1.4) and dividing the re-
sulting equation by dx dy dz = dr we obtain it in the final form
dXx dXy dXz
dy
'
dz
2*
20 THEORY OF STRESS [Ch.
+ Kp dx dy dz -~ ;
Zp dx dy dz '- .
Examining Fig. 11, we find only two forces producing the mo-
ment of the third order; they are indicated by thicker lines m
'S^l**
dx dy
~dT
-
d H-pJ" =o = P ^ (1.5)-
dZ r
dy
Xy
shows that three of these nine functions are pairwise equal to
three others. Consequently, we have three differential equations
(1.5) containing six unknown functions
y, z),
(1.7).
z).
dx dx
A
^
Fig. 12 Fig. 13
Xv = ^cos (vx) + X y
cos (vy) -f- J^ cos (vz).
vz)
= n. (1.9)
Xv = Xx l -f-
Xytn -f- X n\
z
>pXdt = pX *
I
*
,
while all terms of the
equation of (1.8)
first
J
Xv dS + f XP d* = 0.
(S) (t)
Sec. 4] ANALYSIS OF STATE OF STRESS AT GIVEN POINT IN A BODY 25-
Here the first integral is extended over the surface of the body,
and the second integral over its volume. However, the first integral
can also be transformed into a volume integral by the use of the
Green-Ostrogradsky formula
l
(S) (S)
'
dX r dXv d>
(0
The integral on the left-hand side must vanish for any volume r
isolated in the body, but this is possible only if the integrand is
identically equal to zero at all points of the volume t:
Let us suppose that an element Oabc has been cut out within
an elastic body and that X v Y v Z v are the projections of the
, ,
total stress P v over the oblique area abc on the axes of a random
co-ordinate system Oxyz (Fig. 13). Knowing the projections of
the total stress (X v Y v Z v ), we readily find the normal and
, ,
1
The Green-Ostrogradsky formula has the form
dQ . d
[P cos (vx) Q cos (vy) ^ cos
f (4r -j- -f-
(S)
Fig. 14
(1.11)
Obviously (Fig. 14), the normal stress on the area abc will be
obtained as the sum of the projections of stresses X v Yv Z v on
, ,
the normal v.
Vv =
(1.12)
Sec. 4] ANALYSIS OF STATE OF STRESS AT GIVEN POINT IN A BODY 27
+ Yz (
X
(1.13)
(/yi 3
/l
2 + ^2^3) + ^ (^3^2 + ^2^) +
Formulas (1.12) and (1.13) are considerably simplified in the case
of the so-called state of plane stress when all forces are parallel
to a single plane and distributed uniformly in a direction perpen-
dicular to this plane.
Fig. 15
U = (1.14)
(a) (b)
;28 THEORY OF STRESS [Ch. I
Vv = Xx sin a + Y
2
y
cos 2 a + 2X y
sin a cos a, }
2
sin a). J
X, y,
Y
1
X">
Y V Y Z>
' *
(1.16)
V ifv* W
v v> u wv
over any fourth area (at the point 0) determined by its outward
normal v\ in view of .this the quantities (1.16) are sometimes
called elements of the state of stress in a body. One can proceed
further and construct two more areas with outward normals u
and w at the point O\ the stresses over these areas will be
expressed by formulas of the same form as (1.12) and (1.13),
for instance
1
/
1
+ 2X y
l
1
m l
.
After this has been done, we obtain nine formulas [three of the
type (1.12) and six of the type (1.13)] expressing the stress
components over the areas of the trihedral Ouvw
UU V
-
tt
W
U
Uv Vv v W
u wv (1.17)
Uw Vv w W
^ ww
in terms of stresses (1.16) over the areas of the former trihedral
Oxyz.
Nine components (1.16) or (1.17) define the stress tensor and
from this point of view they are called components of the stress
tensor. The formulas of the type (1.12) and (1.13) (as we have
said, there are nine of them) determine the transformation of
the tensor from one co-ordinate system to another. The stress
tensor (1.16) is symmetric, since the components symmetric with
Sec. 5] STRESS DISTRIBUTION AT GIVEN POINT 29
where /, m, n are the cosines of the angles that the normal makes
with the co-ordinate axes. The normal stress on the area in ques-
tion according to (1.12) is
Vv = X x l 2
+ y m + Z,n + 2y mn + 2Z nl + 2X lm.
y
2 2
2 x y (1.19)
Xx ? + Y tf + Zp + 2Y rt + 2Zx
y z tt H- 2X y ty = p Vv
2
.
(1 .20)
the +
sign should be taken on the right-hand side if Vv > 0, i;e M
if the normal stress on the area is tensile, and the is
sign if it
compressive. Then
or ^=' (1-21)
O(5, TJ, )= c2 .
(1.23)
= (X x + X m + A =
2/? l
y
= 2p (Yx + Y m + Y n) =
l
y z
= 2(Z, ,C)
= (Z
2/? ,n)
= 2pZ v .
stress over the area concerned on the co-ordinate axes. Hence the
conclusion: the total stress MP
on the area is perpendicular to the
tangent plane 5Sj to the surface; knowing its direction, we obtain
1
Since there are no terms of the first degree in |, r).
in equation (1.23),
Sec. 5] STRESS DISTRIBUTION AT GIVEN POINT 31
its magnitude P v = MP
by drawing NP _L MN. Now, of course,
we can easily find the shearing stress MT. Thus, Cauchy's quadric
enables the stress distribution at the given point of the body to M
be fully investigated.
It is known from analytic geometry that by rotating the co-ordi-
nate axes the equation of the surface of the second order (1.23)
can be transformed so as to eliminate the terms containing prod-
ucts of co-ordinates in a new co-ordinate system (u, v, w). Then,
we shall obviously have
4- Nff =
2
c .
(1 .24)
The final form of this surface depends on the signs of the prin-
cipal stresses N\, 2 N N
3 and on the associated sign of the right-
,
\N 3
This is the equation of an ellipsoid with the semi-axes
e c c* (1.26)
( *\ ( *\ (
W7 UVi/ \Wl
A^ 2
+ A^ 2
-|Af 3 |t
2 = 0. (1.28)
This equation gives the asymptotic cone separating the above two
hyperboloids and tending to approach them at infinity (Fig. 16).
If the end of the vector representing the normal stress on the
area appears to be on the hyperboloid of one sheet (1.26), this
stress is positive, i.e., tensile; if, on the other hand, it appears
to be on the hyperboloid of two sheets, the stress is compressive.
In the intermediate case it may be directed along the generator
of the asymptotic cone (1.28); in this case the length of the vec-
tor becomes infinite and according to (1.21) Vv = 0. Hence, on
areas normal to the generators of the asymptotic cone there act
only shearing stresses.
If a given area abc (Fig. 14) is a principal area, the total
stress P v acting on it is directed along the normal v and is a
principal stress: we denote it by /V; its projections on the co-
ordinate axes are
Sec 5] STRESS DISTRIBUTION AT GIVEN POINT 33
N) m + Y> =
]
Yx l 4- (Y y 0, (
1
.29)
Zx l + Z m + (Z N)n =
y z Q. t
Adding here the basic relation among the direction cosines of the
normal v
>2_|_^2 =li (1.30)
YyY z
H=
' 1
xY y y
1 7
Y z^z
(1.33)
Y x 1Y y 1Y z
1
z x z y zz
All three roots of equation (1.32) are real they give the val-
1
;
1
This is proved in analytic geometry in investigating the principal axes
of the central surface of the second order.
3-1013
34 THEORY OF STRESS [Ch. I
we find the quantities //, m'i, n't proportional to the direction co-
sines of the normal to the unknown principal area
z N,.
Then equations (1.8a) are simplified and take the form
(1.35)
(1.36)
1
The third equation will be a consequence of Jhe other two, according to
condition (1.31).
Sec. 6] MAXIMUM SHEARING STRESSES 35
Fig. 17
given point are principal and the stresses on them are equal;
this will be the case, for instance, under all-round compression
or tension.
3*
36 _ _
(1.35),
The normal
If we
has the projections
X v = N,l\
THEORY OF STRESS
cerned,
Pl=Vl+Tl
whence we find
+ Nln - (M/ + N^ + A^
2
)
2 2 2
)
.
(1 .38)
(1.39)
dl ~v dl
*" * w '
dm dm
It may be well to note that, in a'ddition to the required solu-
tions, we shall then obtain the extraneous solution 7^ which =
will obviously give a principal area; this solution will easily be
found and dropped later. Writing out the preceding conditions
and using (1.39), we obtain two equations for the determination
of / and m
-N,)m = 0.
Consider the most general case when Afi, N% and Af 3 are all
different;then the preceding equations can be divided through
by the differences of the principal stresses entering in them to
Sec.
We
_
6]
some manipulation, as
ingly, we shall obtain three solutions that are easy to find* The
first and the simplest, in which
/ = =
/tt 0, n= 1,
(N l
N 3 ) (1 2/ 2 ) = 0, (1 .41)
whence
/== m==0: n=
In the second case, from the second equation of (1.40) we obtain
similarly
/ = 0; m=~\ J/2
/i=4=-
J^2
T
1 2
+
the remaining two solutions give similarly
(1.42)
7*,= "'-";
_4
'a- .
Fig. 18
T
'2
+
~
If, for instance,
Sec. 8] SPHERICAL TENSOR AND STRESS DEVIATOR 39
(1.45)
1
Here and subsequently an asterisk indicates the chapters, sections and
parts thereof which the reader may easily omit in his first acquaintance with
ihe subject-matter without impairing continuity.
40 THEORY OF STRESS [Ch. 1
in this case the three principal stresses are the same, as are there-
fore the stresses over all the areas at the given point; Cauchy's
stress surface and Lame's stress ellipsoid appear to be spheres and
the tensor (1.45) is thus called spherical. For the sake of conven-
ience we introduce the concept of addition (subtraction) of ten-
sors. The sum (difference) of two tensors is defined as a tensor
whose components are equal to the sums (differences) of the
respective components of these tensors.
Let us take an arbitrary stress tensor
(1.46)
(1.47)
(1.48)
Sec 8] SPHERICAL TENSOR AND STRESS DEVIATOR 41
H = (Y y - N ) (Z,
- N + (Xx - N (Z - N +
) ) z )
+ (Xx - N ) (Y y
- N - Y Z - ZX X - Y X = Y Z + X XZ +
) z y Z x y y z Z
-f- XX Y y 3A/o Y Zx X =
X XZ Z -f- z y
(1.49)
X x = Nt r y = Af 2 ; Zz = ^ X y = Y = Zx = Q<z
fA = ^
[(Ni
~ ^o
represents the average cubic deviation of a given state of stress
(Ni, Af 2 MB) from the mean hydrostatic stress
,
N = 3- (^ + ^+ ^3). (1.53)
42 THEORY OF STRESS [Ch. I
A* = -|H. d.55)
cepted
Furthermore,
*, = />, cos (/>,*); K jr
=P jr cos(P jr y);
1
The interpretation of the quadratic invariant as the average quadratic
deviation was first given by V. V. Novozhilov (1952) and S. D. Ponomarev
(1953).
Sec. 9] GENERALISATION OF LAW OF RECIPROCITY OF STRESSES 43
thus,
XV =VX \
(1.56)
Fig. 19
areas are equal to each other (see Fig. 19). The law of reci-
procity (1.6) now follows as a special case if the areas are
mutually perpendicular.
Exercises to Sections 2-4.
1. In the wall of the dam represented in Fig. 6 there are induced the
following stresses:
x
y
y
__ AVz Z*7 r = Z7z = AU.
Substitute these expressions in equations (1.5) and check in what conditions
they will be satisfied. Find the load acting on the face Oy (i.e., for *=0);
by using conditions (1.8) find also the load on the face OA (for *=*/tanp).
2. In the circular bar of diameter 2r (Fig. 7), there occur the following
stresses:
0; :
GTAT.
Check whether they satisfy equations (1.5) and in what conditions. Prove
by the use of equations (1.8) that the lateral surface is free of load (i.e.,
/,-y._Z.-0).
Hint: for a point on the lateral surface
3. In the straight bar under pure bending (with the system of co-ordinate
axes as shown in Fig. 8) there occur the stresses
Xy = Yg - = 0.
Fig. 20
h3
*y- r(
J
X =2 =t 2T =0
Check whether these stresses satisfy equations (1.5); find the loads acting
on the faces of the rod.
II
Geometrical Theory
of Strain
u, v, w\ (2.1)
y, *); (2.2)
y, z). .
ig. 21 Fig. 22
c'
B'
-M u j?
Fig. 23
~
e * x __
6 (dx) _
~ du
dx dx'
in the same way, for the edge parallel to the axis Oz (i.e., the
edge M3 in Fig. 22) we have
dw
A B"
'
~&x _ ~dx
'
dx + ^-dx
dx
1
1 + 4^-
]
dx
Since we have confined ourselves to the case of very small de-
(2.4)
48 GEOMETRICAL THEORY OF STRAIN [Ch. II
We now can easily find the shearing strain, i.e., the distortion
of the right angle BAG:
\
dv , du
Fig. 24
dy
(2.6)
dw dv
dz
du dw
"tydZ"
'
'
dy
dv
<T
2e
e y* *=?-.
dy "dT 1
(2.6a)
rT~ dU dw
'
dx
right angle between the 4-x and +y axes decreases. The same
rule applies to the shearing strains in the other two planes. Sum-
marising what has been said above, we obtain the following sign
rule for deformations:
(a) elongations along the co-ordinate axes correspond to posi-
tive linear deformations e xx e yy e zz and contractions to negative;
, , ,
negative.
Formulas (2.6) show that six functions e xx e yv e zzy e yx e zx e xy , , , , ,
dv dv dv
dx
'
dy
'
dz
'
(2.7)
dw dw dw
' '
dx dy
'
dz
The components of the principal diagonal of this matrix are
elongations [see the first three relations of (2.6)]. The remaining
components represent the angles of rotation of the edges of the
elementary parallelepiped M123, shown in Fig. 22, about the x,
y, z axes; this can be seen, for instance, from formulas (2.3) and
(2.4); in the notation accepted therein matrix (2.7) is written as
e xx* a
y*'
a ^ ;
ajry ^yy
a
zy' (2.8)
a w e*z-
41013
50 GEOMETRICAL THEORY OF STRAIN [Ch. II
Fig. 25
and DBDi are then subtended by the common arc DD^ hence
Sec. 10] DISPLACEMENT COMPONENTS AND STRAIN COMPONENTS 51
and p;<0 1
.
-g->0,
In the case of Simple shear along the y axis (Fig. 25 6) we find
in a similar way
ft '/_ 1 dv
?* 2"^T'
**)
du dw
-ted^'' (2.9)
dv du
-5r 57-
(2.10)
4*
52 GEOMETRICAL THEORY OF STRAIN [Ch II
CM=P\ c^ =q
we write down (2.10) as
ry\
i
pz; 1
(2.11)
it dx 4- v dy + w dz
is the total differential of a certain function
<(*, y, z),
i. e.,
~~dxT
= dy^ dx^ 2 '
GEOMETRICAL THEORY OF STRAIN [Ch. II
'
dy
2
dx 2 dx dy \dy *
dx / dxdy
dy
2 '+ Ar2
= ~dx~dy*
x '
^_
d 2 e yy d 2 e zz d 2
ev
[yf_. (2.13)
dz 2 ^ v22 '
'dy dydz
2
d e fr
dx 2 *
dz 2 dz dx
~ oZ dxd wdy
_ *
.
de*y
'
dx '
dy dz
dx dy dz ''Wdy'
we get
d fozx d**z.
dz
(
\
deyz
dx
,
dy
s^ =2 dxdy
'
(2.14)
_i/^i-t
de*x de
*\_ o
<te (
~^~
dx dy dz 1~ c^Arc^y
'
d*zx
'
dyl
i
^
*
^
|
^JC
y^
dy dz dx
1
From this point of view equations (2.15) are sometimes called Saint-
Venant's identities.
GEOMETRICAL THEORY OF STRAIN [Ch. II
dv dv dv
(2.16)
dw dw dw
1
See the above-cited N. I. Muskhelishvili's book.
Sec. 11] COMPATIBILITY EQUATIONS
The quantities
dv du_
dx dy
are the angles of rotation of its edges with respect to their orig-
inal positions.
Equalities (2.16) may be regarded as differential equations in a,
v, w. If the components of strain and rotation are given at all
points of a body, the displacements themselves of the points of the
body u, v, w can be found by integrating equations (2.16). Car-
rying out this integration, we shall, in effect, convince ourselves-
of the necessity of satisfying the compatibility conditions (2.6).
Let us turn to the equations in the first line of (2.16) from which
we shall seek the function u(x, y, z}\ we note that, according to
the meaning of these equations, their right-hand members must be
partial derivatives of one and the same function; hence, the fol-
lowing three conditions should be satisfied:
dexx
dy dz
'
dx dz
de K y de v
===
i
A-
dz A*
dz
~
dz (2.19)
and we have then
de xx
dx
(2.20)
'
dz
.58 GEOMETRICAL THEORY OF STRAIN [Ch. II
equalities but they will contain only the second derivatives of the
strain components. Only six of these equalities are different and
they will coincide with equations (2.15). Thus, Saint-Venant's
compatibility equations (2.15) are, in effect, the necessary condi-
tions for displacements to be determined from (2.6) according to
the given strain components.
Saint-Venant's conditions are simultaneously sufficient for this
purpose, if we consider a simply connected body that has no
through cavities. In the case of a multiply connected body, Saint-
Venant's conditions also permit the determination of displacements
u, v, w by integrating Cauchy's equations (2.6); these displace-
ments, however, may now be represented by multiple-valued func-
tions of x, y, z; besides Saint-Venant's conditions, it is necessary
to introduce some additional conditions to make displacements
single-valued, as it is required by the physical meaning of the
problem.
The case of deformation of a multiply connected body is treated
in detail in Chapter I and Appendix II of the above-cited
N. I. Muskhelishvili's book.
nates of the initial point A of the segment and those of the ter-
minal point B. In the process of deformation these points will
Fig. 26
occupy new positions AI and fli; their co-ordinates are also shown
in Fig. 26. The projections deformed segment A Bi on the
of the {
(2.21 a)
7)!
= + ij 8<y; d = C+8w. (2.215)
gj g = 8g; Tj! 7]
= 8i|; ^ = 8t (2.22)
dx = 5; dy = TJ; dz = C.
<60 GEOMETRICAL THEORY OF STRAIN [Ch. II
dx ^ dz' '
'
dy
* dv * dv dv *
* .
dx ~dz^
^
dw
dv
dw
dv
dw
..
(2.23)
~dx~
(2.24)
e^ n
du
dx
= du du
*
= dv
' dv_ r
dv
(2.25)
dw dw dw
~w "~dz"
since, in general,
dv du dw dv du dw
dx "37" "aF T5T
These inequalities are transformed into equalities only in the ab-
sence of rotation at a given point in a body (co x = 0)^ = 0)2 = 0), i.e.,
in the case of pure deformation (see end of Sec. 10).
The left-hand members of equalities (1.8a) and (2.25) are also
identical in meaning; they are projections on the
co-ordijjate
axes: in the first case the projections of the total stress on the
area with outward normal v(l, m, n)\ in the second case the
quantities proportional to the elongations of the projections of
a segment having the direction cosines /, m, n.
Now consider the elongation of the segment AB=p itself;
since
p p ?
'
dv
-3
dy
m *
z
H ,
'
?
dz
/ dv
n* +- ,
'
du \
\
-r
dx
,
H^
'
dy )
m 2
+ en 2
#/. (2.27)
*xx ex
and ~e (2.28)
62 GEOMETRICAL THEORY OF STRAIN [Ch. II
tion*^ (2.27) and in the second of matrices (2.28) (cf. 2.6a); rela-
tion (2.27) has then assumed a form quite similar to (1.19). The
first of matrices (2.28) was called the tensor of relative displace-
P
2
*
P
= exx $ + e yym? + e, n? + 2^ z y fy + 2e ^ -+ 2,.
yz (2.27')
P 2e
= p, i.e., p
= |/^
_ -
(2.29)
p
where
F&\ C)=6 2
,
(2.30)
F& ^ 9= ^ 2
+ * 1 + *t + 2^
yy
2 2
(2.31)
-xy *
-yy
e *y* = 0,
-yz
(2.32a)
^~^ -yz
Ti
(2.32b)
'xy
(2.32 c)
can write down the equation of strain surface (2.30) in the prin-
cipal axes
> = 62 .
(2.33)
The form of this surface, just as that of the stress surface, depends
on the signs of the principal elongations e it e 2t e 3 If all three .
64 GEOMETRICAL THEORY OF STRAIN [Ch II
i
e xx e yy i
&xx e zz \
e yy e zz ~f~ ^xx^yy^zzr
The four terms in parentheses, which are small quantities
last
of the second and third order, are disregarded; then, taking into
account that dxdydz = di, we have
*x '
yy '
*
Sec. 14] STRAIN DEVIATOR AND ITS INVARIANTS 65
+ *, (2-34)
or
(2.35)
Thus, the dilatation coincides exactly with the first, linear, invari-
ant of the strain tensor (2.32a). It is constructed, just as the
other two invariants, similarly to the respective invariant of the
stress tensor; therefore, all we have said in Sec. 5 concerning the
invariants of the stress tensor may apply to the invariants of the
strain tensor; this will reduce to the formal replacement of the
notations
Xx i
Yy ,
y,
/ z,
-yz (2.36)
-yy
e e yz (2.37)
tyz
where
1
(2.38)
5-1013
66 GEOMETRICAL THEORY OF STRAIN [Ch. II
(2.39)
&xx ^ ^xy
&
f> pv
c (2.40)
(2.41)
The resolution of (2.37) is not merely a formal operation; it
reflects the physical properties of the phenomenon of deformation,
since real materials behave differently, in the sense of resistance
to the change of volume and to the change of shape.
Consider the invariants of the strain deviator; they are con-
structed in the same way as the invariants of the stress and
strain tensors, with an appropriate replacement of notations.
The first, linear, invariant is zero, according to (2.41). The
second, quadratic, invariant, by analogy with (1.49), has the form
+e zx )*}\ (2.42)
-e e xy
-xy
e yy -e (2.43)
Sec
is
or,
15]
_
simpler form
X'
according to (1.52),
= (e,
FINITE STRAIN
e) (e 2 e) (e^
_ it
e)
is referred to the
similarly to (1.51),
*)
3
1-
67
(2-44)
xx
du
xx === -3
dx
a
dv
-3
J
dx
= '
'
Q
P
r
= du
-3
dy
we have disregarded small but finite quantities. One must take
it into consideration in .evaluating all the subsequent conclusions;
they will be sufficiently accurate in the case of small deforma-
tions (when all the strains e xx e yy e zx e xy are very small as
, , . . . , ,
5*
68
_ _
more exact
GEOMETRICAL THEORY OF STRAIN
does not affect their accuracy. It will now be shown how to obtain
relations corresponding to finite deformations in place
of relations (2.6).
Starting with formulas (2.23), we replace their left-hand
bers according to (2.2 la, b) and rewrite them as
[Ch. II
mem-
,
dv
+-3
dw ..
,
dw ,
p_
J pi
(2.45)
(2 46)
.
Sec. 15] FINITE STRAIN 69
dw dw
note that
~"
here e= pl p
represents the unit elongation of the 'segment p.
JC /*/), (2.47)
where for brevity we have introduced the notation
*
2
dv 2
du
P
xx ~~ dx
^L ^ 1
\j \ ~
r \^)
i ( \
"""
i
t
2 [ \dx ) (dx
dv , du . du du . dv dv . dw dw (2.48)
dy dz dz dy
du ,
dw ,
du dw ,
"
dv dv ,
dw dw
fry n AY '
"i
d^ dx f\y dx
fry '
dz dx*
/^^ /) I
/}y
dz dx
/)v-
70 GEOMETRICAL THEORY OF STRAIN [Ch II
(2.49)
/l+2e jrjr
1 + e,,,
of the rod at the origin is fixed (does not rotate in the plane Oxz)\ the same
for the element dz in the plane of the cross section.
2. In the case of the pure bending of the rod considered in the preceding
(2.5i;
2P
It is
requirea:
(a) to find the general expressions of strains according to equations (2.6)
and (2.35) and to check whether they satisfy the compatibility equations (2.15);
(b) to write the equation of the deflected axis of the rod; for that substi-
tute the co-ordinates of a point of the axis before bending in expressions (2.51)
(point "a" of the preceding exercise);
(c) to check whether the fastening conditions at the left-hand end are sa-
tisfied (point "b" of the preceding exercise).
3. In an investigation of the torsion of a circular bar (Fig. 7) one obtains
the displacements
u tyz -f- #y + bz -f- c\ v = ixz ax -\- ez -f- /; w= bx ey-\- k.
It is required:
to select the coefficients a, b f c, e, f, k so that the end section 2=0
(a)
will be fixed in the following manner: the point O
must have no displacements;
the element of the axis dz rotates neither in the plane xOz nor in the plane
yOz\ the element of the axis dy does not rotate in the plane xOy\
(b) to find the magnitudes of the strains (2.6) and (2.35);
(c) to check whether equations (2.15) are satisfied.
Ill
16. GENERAL
In the preceding chapters we presented the theory of stresses
elucidating the static aspect of the problem, and also the theory
of displacements and strains elucidating the problem from the
geometrical point of view. These two theories alone cannot be
instrumental in solving the physical problems of the theory of
elasticity concerning deformations that are produced in an elastic
body under the action of external forces applied to it until stresses
and strains have been connected by a physical law. The physical
nature of this law makes it connect dissimilar features of the phe-
nomenon under investigation stresses and strains.
The theory of elasticity treats this law in the most general form.
Its analytic form, i.e., the general form of functions relating
stresses to strains, is ascertained in the first place:
===
-^x J\ \Pxxi ^yy ^zz> ^yz zx , ^ x y)t
' :=r::
y
, e yy , (3.1)
===
/6
of (3.1) as follows:
Xx = a u e xx -+ a l2 e yy + a^e + a zz }4 e yz +a l5
e zjc -f a l6e xy .
(3.2)
1.
a b
Fig. 27
stance, of mild steels with low content of carbon. The initial por-
tion OA of this diagram is straight, provided the test is conducted
thoroughly, an indication of the proportionality between stress and
74 GENERALISED HOOKE'S LAW [Ch. Ill
unit elongation
Xx = Eexx \ (3.3a>
e yy =e =
zz oe x
-
(3.3b>
or
1
e yy = -r\Y -^Z,+ X )Y
y JC
e zz = [Zz -a(Xx +Y y )}. (3.5)
>-
Fig. 28
we then obtain
(3.6)
where
-yz G
(3.8)
G
Formulas (3.5) or (3.6) in conjunction with formulas (3.8) give
the generalised Hooke's law for a homogeneous and isotropic
elastic body, i.e., a body whose elastic properties are identical in
all directions.
This leads us to a new basic group of equations of the theory
of elasticity
I r ** / J
_ (j
l Z
_ (3.9)
e zx * x
G E "
P J- y 2 + g)
y
Let us revert to equalities (3.6) and add them up by members:
yy
= \(\ + a) (Xx + Y +Z -
y 2)
hence we find
or
y* _
~ __ (1 + 2a)
a) (1
~r
1 -f. a xx*
XX = M -\-2v4w (3.11)
where
X== *= 3 12 )
(l + a)(l-2a) ;
Xy = Kry
Sec. 18] STRESSES EXPRESSED IN TERMS OF STRAINS 79
Xx +Y y + Z = 3X9 + 2p
z (exx +eyy + *).
or
e= + (3X 2tx)8, (3.14)
[A 1 2a
hence
__
G
and further
c-
= ar
Thus, we have the formulas expressing E and a in terms of Lame's
coefficients:
1
See Sec. 7.
80 GENERALISED HOOKE'S LAW [Ch. Ill
2a (3.17)
2a
(3.18)
Zx =
e yz
e zz e
(3.19)
we then obtain
N Q
= W, (3.20)
where
Exe-cises.
Xz dxdydz\
to the work done by this couple we must multiply its
obtain
moment by the angle of rotation; this angle of rotation will
obviously be given by the increment
displacement:
= (Xx *exx + Y y
*e
yy +Z z le zz +Y z *e yz +
The work per unit volume of the body at the point where the
parallelepiped has been isolated is
.
(3.22)
form of work done by the forces when the deformations disappear >.
The potential energy is due to deformations and is caused by them
If we
give elastic virtual displacements to a body, then, accord-
ing to the hypothesis introduced, the work (3.22) done by internal
forces on these displacements must be transformed completely into
elastic energy and give an increment >W:
%A = 8 W. (3.24)
y\ Jf
"
-N > * ^ \ /^ y ~~V.
\ x yy zz
dW . .
dW (3.27)
dexy
a l4 e yz
Yy = a^exx + a 22 e yy
Yz = a^e xx + a 4<2
e yy
(3.28)
-f a 6 ^y .
3. If we
limit ourselves, as before, to a case of very small
deformations, we shall show that whatever the analytic form of
relations (3.1) may be, they can be replaced approximately by
linear functions; in fact, expanding the right-hand members of
equations (3.1) in a Maclaurin series, we obtain, to small quanti-
ties of the second order,
, (3-31).
where the zero index in the function / t and its derivatives shows
that one must put
C xx = 6yy ==...= 6x = y
U.
1
If we solve equations (3.28) with respect to strains, we shall also obtain,
linear expressions for them in terms of stresses.
GENERALISED HOOKE'S LAW [Ch. Ill
Fig. 29
1
This reasoning holds true unless all the coefficients of the strains in
equation (3.31) vanish, in which case one must take into consideration small
terms of the second order in the Maclaurin series; no such case is possible,
however, since in the case of simple tension the existence of first order terms
is readily ascertained (Hooke's law).
Sec. 21) STRESS-STRAIN RELATIONS 87
= _A: + _y + _
"
x
dF dF dF . .
nF(x, y, z, ...)
1
In the case considered, for instance, the magnitude of stresses depends
upon the size of the sectorial cut-out ACBD.
$8
_ _ GENERALISED HOOKE'S LAW
(3.33)
III
or
e zz
(3.34)
since X>0 and |ii>0, as has been shown at the end of Sec. 18.
=a '
(3 - 35 )
Sec. 23]
Now
order
find the same second
dW
ISOTROPIC BODY
_
dezx de yy
== de yy
_
derivative by differentiating in reverse
52<
coefficients a mn and a nm
amn = anm- (3-37>
Qfi A
which are located along the diagonal and ?p
=15 among the
the left-hand member of the equation; the first three terms in the
right-hand member (containing elongations) will remain unaltered;
the fifth term will not change either; the fourth and sixth terms
preserving the magnitude will change the sign; hence, the
equality between the left-hand and right-hand members will be
violated; this will not happen only when 0i4 = 0, 0ie = 0, i.e., these
two coefficients must vanish. If we reverse the direction of the
axis 02, we shall find in the same way that 0u = 0i5 = 0; thus, we
draw the general conclusion that au = a 15 = a 16 = 0, i.e., the nor-
mal stress Xx is not connected with shear deformations. In
view of the complete equality of the co-ordinate axes in an
isotropic body this conclusion also applies to the other two normal
y
j
Fig. 30
a41 =a =a =
42 43 0,
a 51 = #53 =
conclude that
and further
Y = (3.38)
Y y = a u e yy +a l2 (e zz
z 11 zz\ 12 V xx i
~yyn /Q QQ\
Y =a e (o-oy)
Zx = Q>\t zx ,
Y j
an a, *44'
92 GENERALISED HOOKE'S LAW [Ch. Ill
/
y (3.40)
X = y
Xy
= 0,446x3
= 0.
Hence
(3.41)
uv a^e .
(
an a ia) exx sin 2a = 2an aue uv (3 - 42)
Sec. 23] ISOTROPIC BODY
and then the shear angle e uv as the distortion of the right angle
ABC. We have
OA = 1 sin a; OB 1 cos a,
and
AA l
1 sin a (
e yy ), BB l
= 1 cos a e xx .
A^ = AA 2 l
=
cos a 1 sin a cos a ( yy ),
BBc, = BB l
sin a = 1 sin a cos a e xx .
r vy )smacosa,
or
e xx
-sin2a.
94 GENERALISED HOOKE'S LAW [Ch. Iir
.
,
*uv
= P Pi
= for* e yy ) sin 2<x,
au a l2 = 2a. (3.43)
(3.44)
Sec. 23] ISOTROPIC BODY 95
or 6 =
Substituting this value of in the first of equalities (3.44), we
find the expression of e xx in terms of stresses
e ;
Solution
of the Elasticity Problem
in Terms of Displacements
A. Static equations
dZ x
~dx~
Xv = Xx cos (vx) + X y
cos (vy) + X cos (vz)\
t
Yv = Y x cos (vx) + Y y
cos (vy) + Y cos
z (II)
Z& = Zx cos +Z
(I;A:) y
cos (^y) + Z^ cos (vz).
'Sec. 24] COMPENDIUM OF BASIC EQUATIONS OF THEORY OF ELASTICITY 97
B. Geometrical equations
3. Relation between displacements and strains (Cauchy's equa-
tions) (2.6):
du dw , dv
xx X yz i
fi z
'
(jy
l/V ^
fM-
i
dw
t/UC/
yy dv
*
^x /"I i
/) v
The
*
(III)
c^w dv . du
zz dz
'
*y d* '
dy
4. Compatibility equations (Saint-Venant's equations) (2.15):
e xx d 2 yy d 2 exy d /de yz de zx de xy \
(
' '
dy dx 2 dx dy dz \ dx '
dy dz / dxdy
d 2 e yy
_
2
_
L ~~~ ^yz
'
dx dx '
(IV)
dy dz dz dy dz
'
dy \ dy
2
de yz
****
dz 2 dz
'
dz ~dx~~~
^ dezx \
d e yy
dz dx
d;t dy \
C. Physical equations
5. Generalised Hooke's law (3.9), (3.13):
yy
= ^iy y (Z,+ A,)]; ^ = ^-Z r jr ;
[ (V)
! e*y
= ir X r
(Va)
(V)
(V'a)
On the basis of these groups of equations it is possible to
proceed directly to the solution of the general problem of the
theory of elasticity for the stresses and strains produced in an
isotropic elastic body under the action of external forces.
Itshould be noted once again that all the foregoing equations
have been derived on the assumption of* very small deformations;
7-1013
98 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF DISPLACEMENTS [Ch IV
'
**v* v* ^v'
/ d 2u \ \
dx dy
+ '
+ '
dz
'
.
dY v d2
h dzL + r =
i
y i
dZ x dZ
~~dx~ ^ J?L
.
dy
+ + Z = (^P~^/'
^ z
P
\
Xv = cos y
cos cos (vz)\
= ^T cos cos (II)
Zv = cos cos cos
dv
(4.3)
du
dX x , dO ,
d 2u d*u
100 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF DISPLACEMENTS [Ch. IV
==
We note that
d 2u ,
d*v d*w _
~ d Ida .
~^
dv ,
"^ a^
do; \
_
~ dQ
'
()A: (dx dy J dx
(4.5)
(VI)
the body forces pX, pK, pZ are absent, equations (VI) will
If
be homogeneous. In the case of equilibrium, their right-hand
members are equal to zero; dividing all the equations by [i and
introducing the elastic constants M, and k (3.17), we reduce them
to
, <?0
(VI')
COS ^
>
x'
W~~dx dV "dy' dV ~^~
~dF ~dV~ ~~dx
da da
4-
,
jf
COS (Vy)
v
+ -gj
,
COS
.
(VZ).
.
Reasoning in the same way with regard to the other two equa-
tions (II), we obtain the definitive form of the conditions on the
surface:
cos
Yv =
(Via)
cos fr -
cos ( <y2:
17
IF
102 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF DISPLACEMENTS [Ch IV
(4.7)
and, moreover, suppose that the body forces are absent, i.e.,
y, z\ consequently, if we
consider points which in the
(Fig. 32)
absence of motion are located plane P normal to the axis Ox,
in the
all these points will be displaced equally and simultaneously; in
other words, the plane P will move in the direction of the axis
Ox without deforming.
Indeed, the equation of the plane P at rest will be X = X Q It .
6 = du ae
dx
ao _ d 2u a9 _ d 2u
dy dx dy dz a^a^
d2w
dt 2
= 0.
d2u d2 u
or
d2u
dt 2 (4.8)
where
(4.9)
1
Since the distance d between any two planes XXQ and x**x\ depends
on the time t
during the motion,
(4.10)
w w(x, t).
In this case all displacements take place parallel to the axis Oz.
By reasoning as above, we easily deduce that all points of any
plane P (Fig. 32) displace equally and simultaneously, remaining
at a constant distance from the plane zOy. If several such parallel
planes are considered, they will move in a vertical direction. In
the case of a periodic motion we shall have uniform transverse
vibration along the axis Oz. Let us check the possibility of such
vibrations; proceeding from equations (4.10), we have
u ~~ *!L
dx ^ dy _i_^
+*L ^ dz u
o-'
dx~ dy ~ dz U
^_^L__^L_o- '
V W ~
ma)- dx*
-
'
~ ~U
^L-^L-Q
dt* dt*
'
or
dt*
__ dx*
'
where
2
(4.12)
(4.13)
or
(4.14)
or
x t tn / A 1 s*\
where m is an integer.
Equation (4.16) shows that the abscissa x of this point moves
uniformly as time passes; the velocity of this motion is
i/
dx /
V
~~~dt~ T
*
-jg-
sec > T > 50)000
sec
where cp and
\|)
are arbitrary functions; substituting this expres-
sion in equation (4.8), we see that it is satisfied identically. It
can easily be shown that the right-hand member of equation
(4.18) represents two vibrations travelling with the velocity a in
the positive and negative directions of the axis Ox.
Take, for example, the particular solution
u = y(x at). (4.19)
w = B sin 2ic y-
= B sin (* bt\
j^
where
* = 7- (4-21)
108 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF DISPLACEMENTS [Ch IV
(4.22)
\V \<\V\, l (4.23)
f*
(4.24)
v ^
+ 2f*
Fig. 34
(4.24), we obtain
Sec. 28) LONGITUDINAL VIBRATIONS OF A BAR 109
'
r r
p ,p
kg/cm
/ 2
, p
= 7
7.85
^9 kg/cm
1Q33 X981
10 81
3
,
a = 0.3; substitut-
~
remam ver y
Fig. 35
pF dx -^ =
d 2u
-~
dX
dxF, (4.25)
or
r
~~ '
dt 2 dx
and, consequently,
dXx ~~ F
L _ d 2u
*
dx dx 2
1
This equality coincides, obviously, in the problem under consideration with
the first ofNavier's equations (I).
Sec. 28] LONGITUDINAL VIBRATIONS OF A BAR HI
dx 92
a2 T dt 2
~ X dx 2
'
Vu;
The left-hand member of this equation depends only on t, while
the right-hand member depends only on x\ since both of these
variables are independent, equation (4.30) can be satisfied identi-
cally only in the case if its left-hand and right-hand members
are equal separately to one and the same constant number; desig-
2
nate it by A, then equation (4.30) is decomposed into two or-
;
'
dt 2 dx 2
Their general solutions are
A, B
C, D, A, may have their particular values in each term of
y
C = 0;
then
u = Z(A cos Xctf -f B sin Xorf) sin X* '.
(4.33)
Let us suppose that the upper end of the bar is free and the elon-
gation is therefore equal to ** = =
at any instant
for x=l
-^~
of time t.
Keeping this in mind, we obtain cos>J=0 from equa-
tion (4.33). Hence
M = -g"
,
t
in
or
.
\ = jj-,in
where /
may be equal to any odd number: i=l, 3, 5, . . .
,
oo and,
therefore,
00
a=
2 (A- cos -g- a/ + .
sin
-g-a/)
sin -- jc.
(4.34)
(4.35)
But
A, A- sin
(4.37)
Fig. 36
sin
-i
(4.38)
-I
8-1013
114 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF DISPLACEMENTS [Ch IV
dx
dx '
dy
'
1
W
dZ x t
dZ y ,
dZ z
dx
will be inadequate for this purpose, since the number of un-
knowns exceeds the number of equations (I); it will, there-
(5.1)
fore, be necessary to resort to compatibility conditions (IV).
With that end in view we shall further transform equations (IV),
substituting stresses for strains [in accordance with equations (V)
of Sec. 24].
Let us now consider the particular case when stresses (5.1) are
expressed by functions of the first degree (linear functions) of the
co-ordinates of a point or when stresses are constant. By proceed-
ing from equations (V), we can easily show that the second de-
rivatives of strains are always expressed by linear functions of
the second derivatives of stresses (5.1), for instance,
d*Z
dy
2
~ d*Yy
+ *MI
.
dy* I J
'
etc.
11) SOLUTION OF ELASTICITY PROBLEM IN TERMS OF STRESSES [Ch. V
Fig. 37
bases are acted upon by external forces that produce only shearing
stresses and reduce *to two opposite couples; the moment of such
a couple is called the twisting moment (torque). According to the
theory of torsion elaborated by C. A. Coulomb at the end of the
18th century, the deformation of a bar consists in rotations of
the plane cross sections of the bar with respect to one another,
involving no distortion (warping) of cross sections; at all points
of the cross section there occur only shearing stresses T z directed
normally to radius-vectors of the points. These shearing stresses
are assumed to be distributed similarly over all sections (includ-
ing the ends of the cylinder); the system of forces acting at the
Sec. 30] TORSION OF A CIRCULAR BAR 117
M = fT rdF
z z (5.2>
F
(5.3)
X=
z
T2 sina. = T2 2r ; Y 2 = Tz cos a = T z ~
or, according to (5.3),
= (5>4)
Xv = 0\ K w = 0; Zv = 0.
i.e., the lateral surface is free of load.
At the end cross section we have
cos (vx) cos (vy) = 0,
cos (vz) 1 .
< 5 5>
'
ff
X dF=
g
ff
since the origin of co-ordinates is placed at the centroid of the
section. Consequently, the forces applied to the end section reduce
to a couple; we thus really have the problem of the torsion of
a bar.
The moment of the torsional couple about the axis Oz is
where / = -g
TlR 4
p the polar moment of inertia of a circular sec-
is
tion.
Fig. 38
forces. However
great these forces may be (e.g., in cutting
a with pliers), they will obviously cause
wire no appreciable
stresses outside the small region A outlined by the dotted
line.
This postulate being accepted, we can pass on to Saint-Venant's
principle formulated above by means of a simple static transfor-
mation. Our reasoning will again be illustrated by a concrete
example (Fig. 39 a, b) of a rod bent in one case by a weight P
suspended from the end of the rod from below, and in the other
by the same force but pressing from above.
It is evident that the stresses in the region A will be essentially
different in these two cases, but it can easily be seen that the
-difference will be insignificant outside the region A.
: / :
A ,
Fig. 39
the case of thin shells and the so-called thin-walled bars, such
as bars of I-section (Fig. 40ft). Here the dimensions characterising
the shape of a body may differ widely; for instance, the overall
dimensions AB and BC of the cross section may be small as com-
pared with the length of the bar AE, but considerably greater
than the web thickness 6. Take two bars (Fig. 40) built into a
wall in the same manner, one of them with a rectangular sec-
tion and the other with an I-section, the dimensions AB, BC and
122 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF STRESSES [Ch. V
du n dv n dw n
~T~~ \Jy
"
-x U V U (5.7)
dw
w
da
dv
dz
dw
(5.8)
dz
w= ty(x, y),
(5.9)
dy
dz (5.10)
(5.11)
dy dz
n^ dx dy
124 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF STRESSES (Ch. V
- - .
(5.12)
v '
dydz dxdy
Thus, equations (5.11) and (5.12) give the following table of the
.second derivatives of functions (5.9):
() y
dx*
'
dxdz
= T .' y
2
==n .
(5.13)
(5.14)
*
= g* + hy + k. )
\ve obtain
h -f- ix +- e = ix,
or
uf=
v = y = txz ax^-ez-\-f, I
(5.15)
w = = bx
The obtained expressions of displacements (5.15) involve six
arbitrary constants whose presence here can easily be explained.
Indeed, in Sec. 10 we agreed to fix the elastic body under inves-
tigation in order to eliminate its motion dealt with in mechanics
of absolutely rigid bodies and defined above, at the end of Sec. 10,
as a rigid body displacement; in the present problem, however,
we have not yet realised such fixing; the linear trinomials in
formulas (5.15) indicate that it is still possible for the rod to
rotate as a whole about the co-ordinate axes through arbitrary
small angles e, 6, a and to have translatory displacements
c, /,k along these axes [cf. formulas (2.11)].
Let the fixing conditions be realised at the upper end of the
bar: fix the centroid of the section 2 = 0, requiring that
for x=y=zQ
This condition, however, is inadequate as, in spite of it, the bar
can rotate about the origin (Fig. 41); in order to eliminate the
possibility of rotation let us require that two of three elements
dx, dy and dz near the origin of co-ordinates O remain immov-
able. In Sec. 10 [see formulas (2.3) and (2.4) and matrix (2.8)]
we already obtained the expressions of the angles of rotation of
elementary segments; by using them, we can write down the
required conditions.
126 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF STRESSES [Ch. V
Fig. 41
the assumption that the axis of the bar OB will bend along a curve
lying in the xz plane and that plane cross sections of the bar
will remain plane and normal to the deflected axis. From simple
geometrical considerations (as set forth in the courses on strength
of materials) it may be concluded that longitudinal "fibres" of
the bar parallel to its axis receive unit elongations proportional
to their distances from the yz plane containing the axis of the
bar, which is called the neutral plane
Xx = Y, = X = Y = Z x = ^
y z Z2 = .
(5.20)
moment
1
The reader is recommended to confirm for himself that the resultant oi
forces Z z dF is equal to zero.
128 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF STRESSES [Ch V
'
t
'
~~~ (^ 9.]}
dx ~~jT "dy" p dz p
dw i
dv ~. du . dw UV I uli r\
'
dx dx (5.22)
dv '
dz dz ' '
dy
V = -
+ cp (X, Z), (5.23)
dz
' (5.24)
dz '
dx p
dy df <iy
"dx '
~3y~ p
dx 2 dy
2
dxdy
= 0.
(5.25)
= xz ,
bx ey +, ,
k.
The linear trinomials included here represent the rigid body dis-
placement of the bar, just as in the problem of torsion. To deter-
mine them let us fix the bar, for example, at its left-hand end
JZ?
Fig. 42
.(Fig. 42). In the first instance we fix the centroid of the section,
requiring that
for x=y =z=
We fix further the element dz of the bar axis eliminating its rota-
tions in the planes Oxy and Oyz, i.e., requiring that
for ,c = y = *-0 4^ = 0.
91013
130 SOLUTION OF ELASTICITY PROBLEM IN TERMS OF STRESSES [Ch V
(5.26)
w=
Each point (x\, t/i, of the bar moves to a new position after
deformation:
z=z + w = z l l
-
r
(5.27)
Consequently,
(5.28)
*= -4-* + p. (5.29)
Sec. 33] PURE BENDING OF PRISMATICAL BAR 131
tan a = / du \
(-du_
dz )z=z
,
zn
(5.30)
Fig. 44
*=z ;
y=Y;
after deformation these equations become
~2+u =
h h
x= , . ,
TT +
,
(5.32)
This the equation of two parabolas (the upper and lower sides).
is
The radius of curvature of the deflected axis (5.27) and the ra-
dius of curvature of parabolas (5.32) are approximately equal to
pand-, respectively.
p and .
Yz = Zx = X = ^ y (5.33)
Sec. 34] PRISM STRETCHED BY ITS OWN WEIGHT 133
Fig. 45
(5.35)
, z), (5.36)
'
dy dx
dz (5.37)
:
+4r=o.
The method of determining them remains similar to that used
in Sections 32 and 33 and, therefore, without repeating it, we
present the final result:
GO or
v ~- zy ax + ez + /rc, (5.38)
w= -ij;\z + c(x +y 2 2 2
)]
bx ey + k.
for x= Q, y = 0, z =
Jiere must be
dv
dz Jx~
u
// -.- .
OOP-
'o
-p-
yy
xz, v
fj\
>o
(JpjOr
-p- yz,
w
(5.39)
du
1
One may require instead that = 0.
Sec. 34] PRISM STRETCHED BY ITS OWN WEIGHT 135
1
(1
_
a P^
ZQ
^
j
,
(5.40)
z = z = const.
Q
or
ao/r \2
2
where r2 = x +y
2 2
,
rg=x + yg; thus each circle of radius r at the
cross section with centre on the axis of the prism is transformed
into a circle of radius r with r<r The surface into which the .
z z -
-1-.UL f z 2 P -\ - - r2 \
% _L
-Z ^ _. Iz _ ^^or
2 I
2
o
/
*
2
-I- a ( 2 -4-V 2M
i^oJJ
r 2o*
r _(i
r=r n --_2r~ P \
j/-
r0>
2
z ZA _J_ IfL
^01^ 2
f
\
Zo2 _ 72 I
i
0/-2\
or
oj'
(5.44)
v , Y Zv (5.46)
and body forces
X, Y, Z (5.47)
two different systems of stresses are possible:
X'x ,
Y'y , Z;, Y' Z'x , X'y (5.480
and
x\ V2
* ^Y y (5.48")
Both of .these systems must satisfy basic equations (I) and (II):
dxx dx' dX'z
~~dx~
dYi dY'z
(5.490
~JT ~dz~
dZ' dZ' dZ'
dx ^ dy
"^ dz -i-^ ^,
X = X'x cos
v (vx) + Xy cos (vy) + X' cos z (vz),
(?A: dy
<)K
_ar; (5.49")
~5F
az: dZ"
dAT
dx dy dz
dx dy dz
= 0,
dy
= (X X"x] cos ^) cos (vy) + (5.51)
Xx = X =X
y y , etc.
7
Hence, both systems of stresses (5.48 ) and (5.48") are coinci-
dent, which was to be proved.
If the principle of superposition does not hold in a given
problem, we again encounter a multiple-valued solution. This
can be illustrated by Euler's problem of the buckling of a bar.
Indeed, let us take two loads, PI and P 2 each of which is slightly
,
dX,
dx ^ OXy
dy
^ dX,
dz n-K"
dY x dYy dYz
~dT ~^~ ~d^~
~T~
~d~z
hp/ ~~ (I)
dZ x dZy dZ z
dx ' '
dz
:=o
dy
are inadequate, and to solve the problem one has also to make
appeal to the compatibility conditions (IV). In addition, the
conditions (II) must, of course, be satisfied on the surface:
Xv Xx cos (vx) + X y
cos (vy) +X z cos (vz)\
Yv = Yx cos (vx) + Yy cos (vy) ~\- Y z cos (vz)\
Zv = Zx cos (vx) + Z cos (vy) y -f- Z z cos (<vz).
Thus the problem is reduced to the integration of nine equations
(I) and (IV) with six unknown functions. Arbitrary functions
entering in the general solution of these equations must be deter-
mined from the surface conditions (II).
Since the conditions of continuity (IV) connect the components
of strain
,
d*e . )
+ <
(VII)
+ dydz -0;
: +
0+' =-0.
_
dx 1
dx 2 \dx /
r
dy'
2
\ dx ,
dx
"dt (5.54)
1
This derivation is omitted here; it may be found in S. P. Timoshenko,
Theory of Elasticity, New York, McGraw-Hill Book Co., 1934.
Sec 36] BELTRAMI-MICHELL EQUATIONS 141
adding up the obtained results and using (5.53) and (5.54), we get
(^ + ^5 + ^^ = 0. (5.56)
(5.57)
or
(5-58)
(5.59)
In our case (when there are no body forces or they are constant)
six equations (5.56) -(5.59) connecting the second derivatives of
strains of an isotropic elastic body are equivalent to six Saint-
Venant's equations (VI) and can replace them. Now it remains
to convert equations (5.56) -(5.59) to stresses on the basis of
Hooke's law (V) and (V'a) we have
= 0. (5.60)
"T"
2
u ~
3X + z '
2(x dy dz
body forces are not constant. Then the relations become more
complicated and, for instance, the first of equations (VII) is
replaced by
v A -1-
1
~~ _1
d^__ /* dY , dZ\
J = f(X v Z v w) dS.
(S)
On
the basis of the Green-Ostrogradsky formula presented in
Sec. 3 (footnote on p. 25) this integral can be transformed into
an integral extended over the volume of the body t:
'-/(!? (t)
++)*
We calculate the derivatives here
dx
~
^_ dx
__
dx r- dx --
~T-X dx
-- Y
Kx
dx
--
~x dx
+X x e. + Y e + Z e + Y e + Zx e
y yy z zz z yz z
(5.62)
Sec. 37) THREE KINDS OF PROBLEMS OF THEORY OF ELASTICITY 145
XV _
- y'V y"V ** """""
* >
vV _
* - v'
*
V
"~"~ v"V
' 7^V_
- 7' 7" -"U ~"U i
XX = X X - X X \
Yy = Yy
- Yy \
Xy = X - X y y ,
(5.63)
it a' u!'\ v v' v"\
xy xy
pX = 9 X' =
//
P Jf 0,
etc.
Thus, we obtain
J J (Xv u + Y v v + Zv w) dS
= 2 J / J W fa, (5.64)
(S) (T)
where X V9 Y v Z Vy
u, v, w are taken from the first and third lines
,
Xv = 0, Y v = 0, 7V = 0; (5.65)
CO
W>0
at all points of the body; the last equality, therefore, is possible
only if all arguments of the function are zero, i.e.,W
*- =O f
etc.
* = 0; * y , = 0; e zx =Q\
* J rjr
= <Plt*. y)5 *yy
=? 2 (*. ^ j
(6-2)
10*
148 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
Fig. 46
axis Oz; let us suppose that the load does not vary along the
axis Oz in all these cases.
If we isolate a thin element with a load acting upon it in such
a long prismatical body far from its ends by two sections a short
distance apart and parallel to the plane Oxy and imagine that
it works as a separate elastic body, we note that there must occur
or
(6.3)
m ii
37
P q
m n
Fig. 47
directed along their span. The cross sections of the beams distort
in bending (Fig. 476), as has been shown previously (Sec. 33),
and there result elongations and contractions e zz If, however, .
becomes
e xx = ^[XX - oY - y
<? (Xx + Y = 11 - Xx - oK
y }\ [(1 a) y| f
__ .... (6 4) '
where
.
= T^ :
i=Tb- < 6 5)
-
X, = 0; Z, = 0. (6.6)
and therefore
e xy = 2(1 + gl)
x .
>1
y
y (6.40
dYx dY,
"""
dx dy
The surface conditions (II) are simplified as follows:
Xv = X x cos (vx) + X y
cos (vy);
\
(Hp>
'
dy
e = du
dy
'
As can easily be seen, there remains only the first of the six
conditions of continuity (IV)
""
.
2 "^
dx 2 (IVp)
dy dxdy
Hooke's law (V) becomes
e
^xx -( p \'
"yy (Vp)
v
y
Equations (Ip)-(Vp) make it possible to proceed to the solution
of the problem of plane strain.
long as the load does not exceed a certain limit (critical value) ,.
152 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
will be stable and the plate will not deflect in the direction of
the axis Oz. This case is called the generalised plane stress.
The faces of the plate (parallel to the plane Oxy) are assumed
to be free of stress, i.e., everywhere on them
Zz== 0; *z = 0; Yz = 0.
But since the distance between these faces (the width of the
plate) is very small, the stresses can be considered to be vanish-
for the same reason (small thickness along the axis Oz) and,
therefore, we shall assume them to be independent of the co-or-
dinate z. Thus in the problem under consideration the following
-conditions are satisfied approximately for stresses:
(6.8)
They differ from the respective conditions (6.6) and (6.7) for
plane strain only by the condition Z = 0.
2
lel to Oxy.
Thestrain e zz will entail some distortion of the plane faces
of the plate. In view of the small thickness of the plate, how-
ever, this distortion will be very slight. Let us take, for instance,
the case of the pure bending of a rectangular beam considered
in Sec. 33. The load-free lateral edges of the cross section incline
through an angle y (Fig- 446); one can see, however, from for-
mula (5.31) of Sec. 33 that this angle is proportional to the
width of the beam b\ therefore, if b is small, i.e., if the beam
has the shape of a thin plate, angles y WJ"U a s be very small '
1
an assumption is justifiable with regard to normal stresses Z z
Such \
surface is distributed uniformly along the generator (as was supposed in the
case of plane strain). Let us generalise the problem and assume that the load
is distributed symmetrically with respect to the middle plane of the plate
O 0, but nonuniformly (Fig. 48a); then shearing stresses Z x or Z y may
Sec. 39] GENERALISED PLANE STRESS 153
P ( Y (
yy E ^ y (V'P)
x y
.
- {' =
0;
1 ( Z y dz
h J h J
Fig. 48
the stress tensor are replaced by their average values across the thickness of
the plate. For simplicity we shall assume in future that surface loads are
distributed uniformly along generators and, therefore, we shall assume that
throughout the plate
Xz = 0.
154 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch VI
(Up), (IIIp) and (IVp). Hooke's law, however, has different ex-
pressions for these problems: equations (Vp) for plane strain, and
equations (V'p) for plane stress. It is important, however, to
note that the form of these equations is the same in both cases;
the difference is only in the values of the elastic constants which
are expressed in terms of E and a by formulas (6.5) in the case
of plane strain.
Henceforth we choose stresses X Xl Yy ,
X y =Yx as unknowns;
it is, necessary to transform the equation of compat-
therefore,
ibility (IVp) replacing strains by their expressions from equa-
tions (Vp) or (V'p), depending on whether we have plane strain
or plane stress.
Let us differentiate equations (Vp) as required by equation
(IVp) 2 2
d exx 2
1 / d X, d Yy
~~
dy
2
E l \ dy
2 l
dy
2
d 2e
yy
1 /d 2 Y y d2 X (6.9)
U y2
f) A, p
L, |
\
\
f) y'i
1/./V
\
/)
L/^V y%
dy
where p is the weight per unit volume of the body; now equa-
tion (6.11) becomes
Y
A* A
_ Ayx fity
~~ ^
/ /)2 Y
y y
* ""T 2
dx dy \ dx 2
dy
2
dy E>
dx 2 ~"T
~~
^jc dy E l
'
d'*yy ^jry ~~ 1 _ [
^2 ^ ^2
^y &** ^y 1 ~
_
'
dx 2 dx dy E l L dx 2 ""
();c
2 '
r)>;
2 '
dy
2
\
', L dx 2 '
dy
2
dy
2
Hence we obtain
-f J%) = 0, (IV'p)
(6.14)
dx
'
X =-dy
y
' Yy = -L
dx
'
/
6 IQX
10.18)
dx
or --4--
dx ^
dy dy
T
,
dv
^=-^1-;
dy
' y
A = --dx
dv
-I-;'
158 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch VI
px '
+ |J)=0. (IX')
= -- cos (**)
(X)
dxdy
v
_
~ ~*
d 2(
P dy
~
. p dx _
d / dcp \
'
'
dx dy ds ~~~ds (~dy)
d 2 y dx d (6.19)
dy
dxdy ds ~dx*~ds 'ds dx
-^-
dx
~ A
" Y v as
f r
J
ds
(6.20)
160 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
dx dy
Fig. 49
Keeping in mind that the surface loads X v and Y v act per unit
arc length of the contour, we note that the quantities and YW X&
in the right-hand members of equalities (6.20) represent the sums
of the projections of forces applied to the part 5 5 of the bar on
the axes Ox and Oy. If we take new axes ON
and OT (Fig. 49)
directed parallel to the normal and the tangent to the contour
at the point 5 in place of the axes Ox and Oy, formulas (6.20)
in the new co-ordinates n and / will be written as
(6.21)
(6.22)
Sec. 39] GENERALISED PLANE STRESS 161
where 7>> is the longitudinal force at the point S of the bar and
is the transverse force at the point S of the bar.
In formula (6.21) the quantity
represents the derivative of the stress function q> along the nor-
mal to the contour; similarly
dt
~~'
we may consider that
cp
= Af (5)
, (0.23)
dy
~dn
11-1013
162 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
Xx = c\ Y y = a\ Xy = bpx. (6.25)
(6.26)
(6.27)
From equations (6.25) one can see that terms of the first degree need not
1
be included in function (6.24) as they will not affect the magnitude of stresses.
Sec. 41] BENDING OF A CANTILEVER 163
Fig. 50
? = iry
k ~
3
-
Yy = / _
y
= Xy = X _ = y
0;
11*
164 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
(6.29)
y y
1 *
yd
^y jb 2J
0\
Fig. 51
dy*
Y
*
y
=0u = -i
dx*
'
(6.30)
(6.31)
Sec. 41) BENDING OF A CANTILEVER 1$5
>''(*)
= o, (6.32)
Hence
E=-C; D = f--
The stress function (6.31) becomes
(6.310
(6.29a)
a> ? B
dxdy 2
This shows that stresses of the form (6.29) or of the more gener-
al form (6.30) satisfy the equations of the theory of elasticity.
We shall attempt now to select the indefinite coefficients A, B
and E so as to satisfy the conditions on the surface (Fig. 51).
1. On the upper and lower faces
166 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch VI
for x = l
~*"2
requiring that the forces at the end cross section reduce to the
given load Q.
The reader is recommended to carry out the necessary compu-
tations as an exercise and to confirm for himself that the stresses
will be represented by expressions (6.29).
Thus, in this problem the expressions for stresses (6.29) obtained
on the basis of the hypothesis of plane sections are supported
from the viewpoint of the theory of elasticity as well, provided
the load Q is distributed over the end section in accordance
with the law indicated by the last of formulas (6.29). When the
load Q is applied in a different manner, the expressions for stresses
will also be different, but a considerable numerical difference will
be only near the loaded right-hand end of the plate, on the basis
of Saint-Venant's principle (Sec. 31, Fig. 39).
In this problem it is of interest to investigate strains and
displacements and to compare the result with that given by the
hypothesis of plane sections in strength of materials. Since we.
shall now make use of Hooke's law, let us first decide which of
the two cases (plane strain or plane stress) we shall dwell on;
to obtain full agreement with an analogous problem in strength
of materials let us treat a state of plane stress and use Hooke's
law in (V'p) (Sec. 39). Then from (IIIp), (V'p) and
the form
(6.29) we
obtain the following system of equations for determin-
ing displacements:
..
du_ Q_
dx EJ ^ y
dv (6.34)
dy
dv ,
,
du (l-f-a)Q / h2 9
2
\ ,- orv
6 35 >
-57+17^-^7^11-- y )- <
-
Sec 41] BENDING OF A CANTILEVER 167
^
W /
\~2
ly
2
__ xy
2
2~~)
\
,
~r"Ej
_Q_
(6.36)
dv a
~~
EJ
- 2 EJ i'
dy \ / )
(6.38)
The bracketed functions depend: the first one on x alone, the
second one on y alone; but since x and y are arbitrary and
independent of each other, equality (6.38) can exist only if the
bracketed expressions are equal to constant numbers m and n, so
that
= (l+a)l~. (6.39)
A (x) = ~2
tx*
g-
(6.40)
168
We
_ PLANE PROBLEM IN CARTESIAN CO-ORDINATES
lya-My + p.
_
I
[Ch. VI
^ ^
for x=y= u = 0; v = 0.
Now from equations (6.41) we immediately have a
= p = 0.
Itis necessary next to fix (Fig. 51) the supporting section
(against rotation about the point 0); for this purpose the condi-
tion that the initial tangent be horizontal is most frequently
specified:
( X=Q dv
A TT^
y=
for ;
\
ox
1
Cf. this with formula (2.3) in Sec. 10.
Sec 41] BENDING OF A CANTILEVER
or
(6.45)
y]
+ q)A y
(2 + q)
y
,
4 6
du 3Q
4EJ
Fig. 52 Fig. 53
One can easily find the general formula for the inclination of the section
1
with respect to the axis at all its points; this inclination is obviously equal
(Fig. 52) to the shear
dv .du
~EJ
J
f ~~
3EJ
|
"t"
4EJ
~~
_ Q/
3
|
3Q/ ,
fid7 v
"y 'x 27
(6.48)
and check whether they satisfy all the equations of the plane
problem. The stresses Y v are usually neglected; this circumstance
permits us to state a priori that the system of stresses (6.48) in
172 PLANE PROBLEM IN CARTESIAN CO-ORDINATES (Ch. VI
T
h -*-,r
Fig. 54
First we select the general form of the stress function for stresses
(6.48); according to equations (VIII) (Sec. 39) we have
-
dx dy
x), (6.50)
B
-JT-
x 2 y-, + f l (x) y -f/2 (x) ,
(6.51)
where and
/ 2 are arbitrary functions of x. Differentiating equa-
/t
tion with respect to x and substituting the
(6.50) resulting value
in the second equation of
(6.49), we have
hence we obtain
'
V
and equation (IX) becomes (assuming /2 (x) =0)
(6. 5 3>
(6.54)
2 3
The term -yj
JC y in the function 0(3 may be omitted since equality
D
2
(6.52) already contains the term of a similar form -g-.x: y
3
; then-
(6>55)
(6.56)
1
The term Ey, which does not affect the stresses, is omitted.
174 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
(6.57)
(6.58)
24
3
L
2
.
r
is2\-'i 4
Hence we obtain
(6.59)
X = Y x = - Bxf + Cx.
y
__ 24 '
2
_ ,
'
,
^' 24 2
_,,_Q
"
= -f; (6.60)
^-
/,3
Ch = q. (6.61)
^p-+-C
= 0; (6.62)
"
yy
r ^L\JL
h* \ 3 4
y
> ^
I
12 (6.64)
for*= Y
(6.65)
/
" A
2
The last condition requires that the shearing stresses on the end
section reduce to the supporting reaction ( ; one can easily
-^-j
see that it is already satisfied, since the last of equations (6.64)
coincides exactly with the prescribed expression of the shearing
stress [formulas (6.48)]. At the same time we note that the first
i.e., the normal stresses over the end section reduce to a couple
with the moment
4
I==
^TT S~ ^"90"' ( 6<68 )
176 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
HHH I Hill M
Fig. 55
20
(6.69)
If we set here
(6.70)
12-1013
178 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
of the third degree (6.26) retaining the first four terms in it, i.e.,
setting
a =b=c= 0.
Y y = dx + ey\
= X = - ex y fy - ?gx. (6.74>
Here pg is the weight per unit volume of the wall; for brevity
designate it by p, i.e.,
Fig. 57
(Sec. 38) and introduce in them the values of the cosines of the
angles between the outward normal v and the axes (Fig. 57):
Xx cos p X y
sin p = 0; Y x cos p Y y sin p
= 0.
(e -t P) tan
2
p
== 7; (d tan p -f- *) tan p = - -
(e +p) tan p.
Sec 43] TRIANGULAR AND RECTANGULAR RETAINING WALLS 179
*
~-
tan 3
(6.75)
_ *V _ .
tan 2 a
,
x -
/Y
(6.76,
The reader can easily prove that these stresses satisfy the follow-
ing conditions on the surface (Fig. 59):
1. Forjc =
2. For x = a
x_ x = w r_ = o.
3. Fory =
^= 0; K,=0.
12*
180 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VB
is not satisfied completely and on the upper face of the wall there
remain shearing stresses
Yy diagram
Xy diagram
Fig. 58
more so as they act in the upper, less stressed part of the wall
where the stress analysis is not usually required.
Solving this problem by the elementary method, we would
obtain the folloVing expression for the stresses Y y :
_N -/T-T-
,
Af /
\x
a\_
Tj
IZiZ
2 3 a
M
'
V
I
Fig. 60
-=0 (6.77)
182 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
? (x, y)
= XY, (6.78)
A'
(IV)
K + 2X"Y" + AT = (IV)
0. (6.79)
X (IV}
= a ^; 4
X" =- X
2
*, (6.80)
* [IV1
= X *; 4
X" = -\ X. 2
(6.81)
K fiv] _ 2X 2
Sec. 44] BENDING OF RECTANGULAR STRIP 183-
y = es \
(
5 2__X2)2^Q, (G.82a)
cp (x, y) (K cos \x
l + KZ sin \x) (A cosh \y -f- B sinh \y +
4-CycoshXy-f-DysinhXy). (6.8 )
K K A, B, C, D and X.
: = and x=l
no normal stresses occur; in other words, we require that
for x = (
There will still remain shearing stresses at the ends, and the
loads on the long sides of the strip will be balanced by the
forces Qo and Qi applied at the ends of the beam; these forces
*can be considered to be the reactions of the supports of a simple
beam of length /.
Conditions (6.87) will be satisfied if we set in the particular
solution (6.86)
m=l, 2, 3, ..., /i
1
:
, (x) dx = 0, f
t2 (x) dx = 0. (6.89)
1
The value m= (J is omitted for the time being and will be considered later.
Sec. 44] BENDING OF RECTANGULAR STRIP 185
-y- 4-
\
J
.
t. . <
D n A cosh
2^.]
= /,(*);
-
mcos + C^ J =
Considering the second and fourth of these equations, we see
the necessity of the restriction imposed on the loads ti and t 2
by conditions (6.89); indeed, the corresponding integrals of the
left-hand members of these equalities taken between the limits
(0, /) vanish. In any case, the above equalities show that the
functions appearing in the right-hand members must be repre-
sented approximately in the interval (0, /) by Fourier series in
imtx , tnnx ,. f _,, , ,,
sin and cos The greater the number of
j -y- respectively.
terms in the stress function (6.88), the more exact this repre-
sentation. To determine the coefficients of the series we follow
the general rule: we multiply both members of the equality by
b-jry kllJC
sin
y
or
j- respectively
and integrate from
cos to /; in
27 .
sin
Knrt knh
KT\* .
^^^
+ C^ (cosh -y-
, x>
H- A smh
/ i- i i L.
/
(6.90b)
.
mh-^
,
h cosh
/JTC/l
-y-
= |-
. , ^7T/Z\
1
2
J
C 2
t{
/ \
(x) cos
k^X
-y-
~J
rf^:;
Fig. 61
the stress function we can find stresses at any point of the strip
according to the formulas
*y x ~ _. dy*
'
vy _
~ jl- dx*
J
dxdy
'
.
(6.91)
K (1V)
=O
their general solutions are
(6.92)
will now be
(6.93)
for y =
(6 ' 94>
for x= Xx
for x=I X x = (a/ + p) (Ay + 5), )
(6 95)
'
*im (*)
= a \m cosh z + b }m sinh z + c }m z cosh z + d lm z sinh 2,
(*)
= a cosh z + 6 2m sinh z + c^z cosh z + d2m z sinh
<2m z,
= a 3m cosh z-f 63m sinh z + c3m z cosh z + d^z sinh z,
(6.97)
(6.98)
(z), (6.99)
where A m Bm
, , Cm, Dm are new arbitrary constants. The stress
function (6.88) becomes
n
? = 2 sln T L}
2 -
'ni' (6.100)
(6.101)
a km">
where the index k is given the values 1, 2, 3, 4 in succession.
These equations can readily be solved in general form; the fol-
lowing table presents the expressions of the coefficients thus
obtained:
'1m
192 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. V
(6.102
I
/mJt\2 -r- .
sin =
mnx
cos
(6.104)
Sec. 45] ONE MODIFICATION OF FILON'S METHOD 193
(6.105)
sin
(6.106)
Z ~
_mny
I
= - = .
sin
cos
13-1013
194 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch. VI
(6.108)
2 1 4,*;*
BT 14,
(6.109)
x
where Ffe m (2) (k=\, 2, 3, 4) denote the following antideriva-
tives of the functions t|) ftm (z) [see (6.97)]:
lp
*m (?) = -
I(**m c km ) cosh z + (akm -f- d t j sinn z +
Sec 45] ONE MODIFICATION OF FILON'S METHOD 195
v = 21 VI 1
(6.110)
= c- sin
(6.111)
cos
13*
196 _ PLANE PROBLEM IN CARTESIAN CO-ORDINATES
which
(6.112)
cp (x, y)
= sin lxY l (Xy) -f cos X*K 2 (Xy), (6.1 14)
where
Y l (Xy)
=A (X) +B
cosh Xy (X) sinh Xy +
+C (X) Xy cosh Xy +D (X) Xy sinh Xy,
(Xy)
== ^ cosh Xy +
(X) , (X) sinh Xy + (6.115)
+ 00
X, = jp =
'
2 2
+ 00 + 00
yy =1^.
= Joo X2cosXjcK,(Xy)dX- J
X 2 sinXjcK2 (Xy)dX, (6.117)
oo
+ 00 +00
AT =
2 2
X sinXjcr,'(Xy)rfX-- X cosXjcK2(Xy)rfX.
y jj|y= jf f I
Here single and double primes denote the derivatives of the func-
tions Y\(hy) and Y^y) with respect to the argument "ky.
If loads <]i(x), ti(x), </ 2 (*), t 2 (x) are given on the upper and
lower faces (Fig. 60), the boundary conditions will be of the
form similar to (6.90a) (the limits of integration are omitted
for simplicity in writing):
for y = 0:
for y = h:
[(6.118)
Y^= 9l (x)
= Jf X cos lxY (XA) 2
l
rfX
Jf
X 2 sin Xxr 2 (XA) rfX,
X =
y <, (jc)
= J X sin XjcKj (XA) dX
2
f
2
X cos XxK2 (XA) dX.
C(X), D(X):
(6.119)
I, C(X),
198 PLANE PROBLEM IN CARTESIAN CO-ORDINATES [Ch VF
rfX. (6.120)
00
/ X
( )
27 J COS Xjc J / (a) COS
[
Xa
flfaj
d\ +
sin X * / (*) sin Xa dL (6- 120a>
Setting now f(x) = ti(x) and comparing with the fourth formula
of (6.118), we obtain
4-co
\ r . ../
If
27lX 2
4-00
* ^
ZTtA^ J J
CO 00
In the same way, by the use of the first two formulas of (6.118)
we find four boundary values of the functions
Sec. 46] STRIP OF INFINITE LENGTH 199
Fig. 63 Fig. 64
Fig. 65 Fig. 66
r + (e, + ^ dr) +
-f- (r dr) d6
- 9,rd9 = 0.
Cancelling out the obtained equations and disregarding infini-
tesimal quantities of the third order, we 'get
" Q
dr r d0
2/? (Ipp)
T "W "7"
r2
(VHIp)
'^
dr\r c
= = yM*, y);
Fig. 67 Fig. 68
which must be satisfied by the stress function cp(r, 0); by the use
it we calculate the stresses
of according to formulas (VIIIp).
We now fix attention on deformations and displacements
(Fig. 67). Let us designate the displacements of points along
the r axis by u, and along the 9 axis by v, the elongation along
the r axis by e rr along the 6 axis by e^
\ and shear, i.e., the
distortion of the right angle bad, by e r^
Let us now examine the shear. One sees from the drawing
(Fig. 68) that the shear is expressed as
*r , = (p-a) + 7 ;
(7.4)
next we note that
dv da du
dv v 1 du (Hipp)
r
Hooke's law (V'p) for the case of generalised plane stress has
the previous form, only the notations of stresses and strains vary:
(Vpp)
"
2(1 + )
E
In the case of plane strain E and a should be replaced by
and GI according to formulas (6.5).
~'' ~
(viirp)
cp (r)
=A 1 n r -f- Br 2 1 n r + Cr + D. 2
(7.5)
(7.6)
for r =a _r = -}-p l
or
for r= b (7.7)
b2 a2
Sec 48] PROBLEMS IN WHICH STRESSES ARE INDEPEN. OF POLAR ANGLE 207
(7.11)
Hence
I u
- U 3)
i
n
*V -=E\ -
1
I du u
<s\ \ dr
+-
1
r dr
dr
-=< 1
+<>) j5l ^i -
<
7 - 13a >
dr dr
whence
_
dR r
__l ,
2(/? r _ _
-C.)_
for r =a R r 0;
for r == b R r
= 0;
b b
J
e (J
dr = 0;
fe>rdr
= M.
a a
. _
C= -JSJ-
(b
2
a2 +2 (b
2
In b a2 In a)],
where
N = (b* a 2 )2 4a 262
(in ^ .
4M ( a*b* t
b
(7.14)
14-1013
210 PLANE PROBLEM IN POLAR CO-ORDINATES [Ch. VII
i
___
----
a P '
v
h '
'
Effect of Effect of
elongation err elongation egg
Fig- 71
dr-\- e rr dr
fi
716):
but
W _ = e^ + dd' -
dr (r + dr) d6 e^
therefore,
8(rfO)
= 51 (rf
straight rod the stresses X,/, corresponding to the stresses R r of the above case,
were absent [formula (6.28) of Sec. 40]; consequently, there is no mutual
pressure of individual fibres in a straight rod under pure bending
Sec. 49] EFFECT OF CONCENTRATED FORCE 211
(7.15)
14*
212 PLANE PROBLEM IN POLAR CO-ORDINATES [Ch. VII
indefinite
form-g-;
but the point C and even the small semicircle
f /?_,pcosflde =A (7.16)
*
It may be shown that stresses (7.15) are obtainable according
to formulas (VIIIp) by the use of the stress function
<f
= j-r^s\n^. (7.19)
dr*
"-"' dr
cp(6)
= -
and then
cp
= rf (6)
==
-|
/-0 sin + Ar cos + Br sin 0.
n 2P cosO^
(7.20>
but
1214 PLANE PROBLEM IN POLAR CO-ORDINATES [Ch. VI [
therefore
cos
(7.21)
A 7
8 ^
m C /7
Fig, 72 Fig. 73
Yx =R r
sin0 cos 0, (7.22)
j
Yy =R r s\n*Q. )
we have
x
^ _. y2)2
2P
y
22
2
)
'
(7.24)
1
We may recall that the points inside the semicircle run are excluded.
2
Cf. formulas (1.15), of Sec. 4.
Sec. 49] EFFECT OF CONCENTRATED FORCE 215-
Fig. 74 shows
the diagrams of these stresses. In practical cal-
culations, itthe custom to distribute the action of the load
is P
at a given depth x under a certain angle; in the conditions of the
problem considered we have
Thus, for calculation one may assume conditionally that the load
P is distributed uniformly over an area about 1.6* long which
corresponds to an angle of 0~38 (Fig. 75).
diagram
diagram
Fig. 74 Fig. 75
Fig. 76
that
rdft
= cos 2
and similarly
The meaning of the angle 9 is clear from the drawing (Fig. 76).
If load extends from point a(9 = 9i) to point &(9 = 9 2 ), by
the
.summing up the stresses due to its individual elements, we obtain
the final expressions of stresses:
(7.26)
(7.27)
Sec. 50] WEDGE LOADED AT THE VERTEX 217
Fig. 77
term satisfies them too, since it is obtainable from the first term
by rotating the polar axis through 90, i.e., by replacing 9 by
are chosen, just as in Sec. 49, from the condition that the normal
forces
a
-a
-fa
/
-fa -f-a
k cos 6
J cos 2 dQ +k sin 6
J sin cos 6 rfO = P cos p,
k cos 6
Ja cos sin dQ +k sin
J sin 2 rf6 == P sin p.
of
-fa
*-// (7.30)
B
then =f= P, and this direction is not coincident in the general
case with the direction of the load P. If we put 0. = -^, we shall
Fig 78
-n>
Hence
2P cos(0
R r
R
Here (6(3) is (Fig. 77) the polar angle measured from the
direction the load P. Thus, solution (7.20) can immediately
of
be applied to the case when the load P (Fig. 72) acts not nor-
mally to the plane AB
but at any angle; one chould only re-
member that the polar angle 6 is always measured from the di-
rection of the load P. Consequently, this solution can be em-
ployed in the case of a load applied in the direction AB\ the pre-
vious solution (Sec. 49) also applies to the problem of a con-
tinuous tangential load over the plane AB. Combining the cases of
normal and tangential loadings, we can solve the problem in-
volving any inclined continuous load.
As illustrative examples of the analysis of a wedge consider
two cases: compression of a wedge (Fig. 78a) and bending of
a wedge by a force applied at the end (Fig. 786).
220 _
In the first
PLANE PROBLEM IN POLAR CO-ORDINATES
case
J3
=
_ [Ch. VII
where
X* X
= (2<x + sin2a)
2/>
jt
'
* li i-=-*
But (y is the polar angle measured from the direction of
0j
Ihe load; consequently, solution (7.15), where k ha* the above
Sec. 50] WEDGE LOADED AT THE VERTEX 221
value, holds true in this case as well. The stresses over the sec-
tion mn, as previously, will be found similarly to (7.22) by the
use of equations (7.15) but bearing in mind that the angle 61,
as we have just said, should be measured from the direction of
the load as shown in Fig. 786. The stresses over the cross sec-
tion mn in this case are expressed by
(2a sin2a)
(7.32)
sm<,
.
more, since
= 0, (7.34)
cp
= r/?9, (7.36)
V2<? = --0"+(r#"+3/?' + 4) e -
<7
- 37>
or
V2<p
= ^-f <]> 2,
where
<h=>>e", (7.38)
,fc
= r(/?" + 3-. + (7.39)
)e.
Hence
V 2 ( V 2 cp) =V +V 2
<h
2
<j; 2
.
(7.40)
V2 and
Sec. 51] GENERAL SOLUTION OF PLANE PROBLEM IN POLAR CO-ORDINATES 223
R r
R
we replace R by /?" +3 -f--^- leaving 6 unaltered.
In this way we obtain
/?0
IV
+ 2AQ" + BQ = 0, (7.41)
(7 ' 42)
e+.e=o; (7.43)
2
(i)
V +(!)' e = -
(7-44)
f
\
6. we
(A -n-}
get
2 -
224 PLANE PROBLEM IN POLAR CO-ORDINATES (Ch. VII
e"-f/tt
2
e = 0; (7.46)
(4)'-
2m *
(4)'
=0 -
( 7 46a >
-
/?' V + 6/?'V + (5 2m 2
) /? r
2 2
(
1 + 2/n 2
) /? V+ ( 1 m 2 2
) /? = 0. (7.48)
/? =r n
.
(7.49)
Now consider the cases when (1) m = 0, (2) m=l. In the first
case equation (7.46) is replaced by 6" = 0, whence
e =D +D 01
6 02
.
r\ rlnr;
i lit
;
ylnr.
Its general solution is
/? = C r+C01 02 r In r +C 03 y+C 04
~ In r.
According to (7.36) the function cp is
To
= (<V + C 2
02 r
2
In r +C +C 03 04 In r) (D01 8 +D 02 ). (7.52)
Ti
= (C u r + C y + C r + C
3
12 13 14
r In r) (D n cos fl
+D 12 sin 0). (7.53)
(4)
/=0 or A=c & (
7 -54)
15-1013
226 PLANE PROBLEM IN POLAR CO-ORDINATES [Ch. VII
V
(B =0,
I
i.e.,
B = c2 & (7.55)
cp
= 0.
A ?!/?
= 0; B CiR = 0,
"V 3 + 5/?"r 2
R'
Here Ci and c 2 must be chosen so that equations (7.57) will have
solutions in common, and later use can be made, of course, only
of these common solutions; let us seek them in the form
R = rm .
(7.58)
-le = 0, (7.59)
^V + tf'r-^/^O,
V /?V + = 2 2
1) ]/? 0.
\
( W >
[1 (/rc j
The solutions
r
m and r~
m
(7.61)
will satisfy both equations (7.60) (since the latter are not altered
when m is replaced by m). The solution of equation (7.59) is
whence
2m = (m l)*
= (m I)
2
/
2
,
Sec. 51] GENERAL SOLUTION OF PLANE PROBLEM IN POLAR CO-ORDINATES 227
If
k=(m 1)1. (7.62)
'
(/I
2
1)2
= (/H* 1)2; J
(7 ' 65 >
<p
= (D 01 cos 6 + Z)02 sin +D 03 cos +D 04 sin 0) (C01 r +C 02
r In r).
To
= (Q)ir + C 02 r In r) (D 3 cos +D 04 sin 0). (7.68)
We proceed to the second case m=l. Equations (7.65) have the
roots rti=l, n 2 = 1 which
give two solutions:
r,
7. (7.69)
15*
228 PLANE PROBLEM IN POLAR CO-ORDINATES [Ch. VII
e= D n cos 20 + D 12 sin 20 -f Dj + DH .
(7.70)
<P
= (^V + C 2
13 ) ( Ai cos 29 + D 12
sin 20 +D +A
13
6 4 )-
<P
= 1<PO + + + + +
2fo l<Pl 2<Pl 1?0 2
where
cp
=C r +C r01
2
02
2
In r +C +C 03 04 In r,
!CPO
= (Coir + Co r 2
2
2
In r + Co + Co 3 4 In r) 0,
Ti
=
To
= (
coi r + C^r In r) cos 0, (7.72)
ffo =(
rlere one can set C 3=C 1 3=C / 13 = 0, since the expressions of the
unction <p corresponding to these constants yield no stresses. The
expression of function (7.71) was found by J. H. Michell except
Sec. 51] GENERAL SOLUTION OF PLANE PROBLEM IN POLAR CO-ORDINATES 229
For cp
= Co4 In r0
= Co4 Inr 1
For cp
=C 02 r In r sin
* =
/? e = C 02
(y
sin +7 cos
0)
.
For cp
=C 02 r In r0 cos
->
?r
= C^/ln
(7 COS
02
6
r\ 2 Inr .
Sin
n\
0J
,
/?o
^ Q2 (7 cos - sin
0)
.
i<Po
=C ol r 2 4- C02 r 2 In r + C04 In r.
In problems of Sections 48-50 the functions icpo and 2 <po are used
partly [cf. formula (7.19)]. The remaining terms of function
(7.71) enable us to solve many important problems related to
the analysis of a wedge and closed ring The stresses correspond- l
.
1
See S. P. Timoshenko, Theory of Elasticity, New York, McGraw-Hill
Book Co., 1934 and P. F. Papkovitch, Theory of Elasticity, Moscow, Oborongiz,
1939.
VIII
Torsion of Prismatical
Bars and Bending
u = ry
1} = TJC2.
(8.1)
accordingly,;
we assume that
, y), (8.2)
1
This will be proved at the end of the present section.
232 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
The first two of equations (VI) are satisfied identically while the
last equation becomes
0. (8.4)
Yv = Yx cos (vx) + Y y
cos (vy) + Y, cos (vz), (8.5)
Zv = Zx cos (vx) + Z y
cos (vy) -f- Zz cos (vz). <
!L
i *-
o
^-M
1
i
(8.6)
*.-<(*
(8.7)
now (8.6), we observe that the first two of conditions (8.5) are
satisfied identically and the third assumes the form
Fig. 80
this point; this condition follows directly from the law of reciproc-
ity of shearing stresses since the lateral surface of the bar is
free of forces.
Let us now express the condition on the contour of the section
in terms of the torsion function <p by the use of (8.7). Multiplying
by the differential element of arc of the contour ds we obtain l
,
(3*
But from Fig. 80 we have
ds cos (vx) = dy\ ds cos (vy) = dx,
The arc is assumed to grow,
1
i.e., ds>Q, when moving to the left as
indicated by an arrow in Fig. 80.
234 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
or
The condition on the contour of the cross section (8.8) may be put in the
alternate form; we write it as
-
cos (vx) + -- cos (vy) = y cos (vx) x cos (vy).
The left-hand member represents the derivative of the function (p along the
normal to the contour [cf. derivation of condition (Via) in Sec. 25]. To trans-
form the right-hand member we observe that from Fig. 80.
cos (vx) = cos (ty)\ cos (vy) = cos (tx),
y = r cos (ry),
where r is the radius-vector of a point M of the contour.
Therefore,
*L = rcos(rO. (8.10)
This condition defines the value of the derivative along the normal
at each point of the contour; thus, the problem of finding the torsion function
cp(#, y)} is reduced to the classical Neumann problem of determining a harmon-
ic function according to the values of its normal derivative prescribed on
the contour.
'
dy dx
or
By integrating, we get
*
y
<|>(jc, y)= **~^ +C, (8.12)
<|>
= -
+ C, (8.13)
*.-,*($-,).
236 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
^_ i i __, v| .
do
(F) x, y"
since (*/i) and G)(J/O) the values of the stress function at two
points of the contour are identical according to (8.17). In
exactly the same way it can be shown that the other projection
Y T of the resultant force vector vanishes. This proves that the
above solution of the problem corresponds to the pure torsion
of the bar.
In order to complete the solution of the problem it is neces-
sary to express the twisting moment in terms of the tangen- M t
M=t
f(Y z x-Xz y)dF=- fx^-dF- fy^dF. (8.18)
(F) (F) (F)
x ^ dF =
.~ f
x
JdO J
Sec. 52] TORSION OF PRISMATICAL BARS 237
(8-19)
lf y) and *(* ,
y)
Af,
= + 2 / ftodF, (8.20)
(/=)
or
2 2
(8.21)
* = W> <
8 22 )
'
where
K= f (2<|>
jc
2
y
2
) dF (8.23)
(F)
hence
H= ;
f=*
y)
= const = C.
Substituting this in the equation of the section contour (8.12),
we get
x 2 -f-y 2 const. =
Consequently, the cross section is a circle. Coulomb's theory is
not applicable to any other cross sections and the cross sections
will warp in torsion.
dU^^dV^.
dx '
dU _ dV
dx
dy dy
W=tU(x> y)
_u (x, y) = * + +C
2
y2 1
(8.25)
D2 =D2
2
2 (x + /y)
2 =D 2 (x
2 - y + i2D xy.
2
) 2
Here
Assume that
or
or, ultimately,
1
It is taken into account here that if the function U+iV is analytic,
V iU is also an analytic function, since the Cauchy-Riemann conditions are
identical for them.
240 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
Setting here
-%- = a 2
,
^= b\ (8.28)
i-D,
we obtain
4 +"F
-&=l,
L
(8-29)
Hence
w = 2D 2 y = il=L T*y. (8.30)
(8.31)
Finally, we correlate the stresses over the section with the twist-
ing moment M t :
(8.32)
But
f f x dF= J^
. 2
f fy2 dF=Jjc
where J x and J y are the* moments of inertia of the ellipse about
its axes:
*
Jx _ nab 4
3
.
.
Jy
GJp*
where J p is the polar moment of inertia of the section. In order
to compare these expressions let us correlate the quantity
Jp
_ /'x
-- I
\
/
Jy
_ "
3 3
F* a2 +b 2
a 6
i.e.,
T74
WJp 39AJp
and, consequently,
^ = -^r39.4y p
.
(8.33)
^r + ir. (8.34)
Hence one can see that the stress T, as a function of x and */,
(8.35)
16-1013
242 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
^ b*
~ '
(8.36)
'
a, \a t a
and the ellipse (8.36) will be more elongated than the ellipse
(8.29) of the section contour. Equation (8.36) gives a family of
h /|2
stant), the stress (8.35), which is constant for all points of the
Fig. 82
At the points A and B the stress (being the same at all points
of the corresponding ellipse /) will be smaller; we have for them
a,
= a == a?C\ C = -L ;
2
T r*/" fl^ ^. T
Sec. 53] SAINT-VENANT'S METHOD. SPECIAL CASES 243
and assume
cp
= V= C 3 (3* y
2
y
3
). (8.38)
Then, according to (8.25), we obtain the equation of the contour
or, setting
C3 = |/3 3 ; C = y>i, (8.39)
we have
D 3 (x
3
3*y
2
) + x* + y = D v 2
(8.40)
This equation represents a complete hyperbola of the third
order with three asymptotes whose directions will be obtained
in accordance with the general rule
by equating to zero a set
of the terms of the third degree x 3 3xy 2 x (x -f- y 3) (x = Y
y 1^3) =
0; hence the equations of three straight lines parallel
to the asymptotes are
;c 0; = r
* 4- y 1/3 0; x y |/ j=0. =
Let us take three straight lines parallel to these:
b. (8.41)
+ D b (2a + b)x D ab = z 3
2
0. (8.44)
+ 2fr)=l; D 6 (2a + 6) = 0;
3
3D 3a=l; 0^ = 0,
16*
244 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
Hence we find
4_
3 (8.45)
(x a) [(x + 2af 3y
2
]
=
or
(x 0)C* + 2aH- yl/3)(jcH-2a y |/1j)
= 0,
w = -^- (3x 2
y y
3
), (8.46)
Fig. 83 Fig. 84
or
which are the three altitudes of the triangle of the section, re-
main in the original plane of the section, since w=z=Q for all
and hence
X" = ,
\
I
(8 48)
'
cp
= jcy. (8.49>
'
b
(8.50>
.
. X, = or - = y; (8.52)
Y2 = or
1*-
= -*. (8.53)
COS I
-y-j
w-D m smh-jr
m u mnc mnx
cm -b r^
-nine nk' \ -
s-> ( i _i_
-(Dm cosh
sin
) -^.
m=l, 3, 5, ...
where
(8 .55)
mT =-2
n bn T
p *s
J
m +""
l
or m = -o-o-( ir"2
u
,
3
w ,
3
cosh
m=l,3,5,...
(8.56>
Further, by using the Cauchy-Riemann conditions (8.11), we seek
its conjugate function ty(x, y)\ we have
,e
'"
sin cosh =_
m 2o cosh i
-J
(1) 2
=r cos -^- cosh
m 3
cosh
m 3 cosh
1
Here the substitution is made
m-i
^~- = (-l) 2
for m-1, 3.5, ...
"248 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
<2 (y)
=i
C is an arbitrary constant.
We obtain finally
-* ,
W vi, ,,
/w +l
L
cos cosh
_
m ,
3
COSh
z^
(8.57)
We calculate now the stress function according to formula
C = -j-to make
b2
(8.15) (we take here this function vanish on
the contour) :
.
1X 9 1 mnx
- mny /0 rov
(-D 2
,
cos cosh-- i
.
(8.58)
/
m=l,3, 5, ...
^^^1
-2T
m = l,3, 5, ... J
= |itWc T ___ ^
Aq V f 1 64 6 *
* u ^^^ 1
.
tanh-jy-
m=l L f 3, 5, ... J
Hence
*= < 8 59 ) -
lifer-
where
1_^1A
m=l, 3, 5, ...
=
For a very narrow rectangle we assume approximately that
= 0, and then
a = ~o 0.333 ... .
^-g-SHD" m ^'os^slnh^, 2
cosh -jr-
J
( 8 61 >
'
26u
m+ i
__^^
OU rwix
tr
Kz at
r
\
1
r\
2;c + Ur ^J
i
i>, (1)
2
X. 1
/
V
-t \
/
O
2
1
mTCC
sin
. til
-r cosh
i
"'. (8.62)
9 . mite
[
M2cosh
_J
or by using (8.59)
Tm^ = (8.63)
-^>
where
u
m o
2
cosh -j
The stress at the point B (Fig. 85) will be obtained from (8.61),
setting A; = O, y = 4-:
'
2b
or
(8.64)
250 TORSION OF PRISMATICAL BARS AND BENDING [Ch. VIII
where
m-f 1
(8.66)
Cyy
-C - zz U, Cx y -
dw dw
dx
dexz
_
dy
(8.70)
Pi Pa
dy*
'
p, p2
membrane:
Jx*
= -&. r
_ 77'
P
(8.71)
dy
(8.73)
then
z = <$>\
(8.74)
Fig. 86
t
(T}
=X z cos (xT) +- Y z cos (yT)\
254 TORSION OF PR1SMATICAL BARS AND BENDING [Ch. VIII
but
cos (xT) = 4- cos (yv),
cos (yT) = cos (xv).
t
( T)
= d* CQS ( ^ ) _j_ CQS (yiy)i
or
(8.75)
111
Fig. 87 Fig. 88
Taking into account that the gradient (the tangent of the max-
imum angle of inclination of the surface of the membrane at
a given point to its initial plane) is proportional to the close-
ness of horizontals, we see that the system of horizontals gives
a visual picture of stress distribution over the section of a twist-
ed rod. Hence, one can draw, for instance, the following rather
general conclusions.
(1) If the section is bounded by a simple contour which is
everywhere convex, the closeness of horizontals increases as the
contour is approached, and it is there that one should expect the
greatest stresses.
(2) If the contour indicated in (1) above has two axes of sym-
metry and is somewhat elongated one of these
in the direction of
axes, the closeness of horizontals will be maximal at the ends of
the minor axis as, for instance, in the case of an ellipse or a
rectangle. (At the ends of the major axis the possibility of the
membrane buckling is restrained by its smaller width, measured
perpendicular to the major axis.)
(3)the cross section has the form of a very narrow rectangle
If
where side b is considerably smaller than side c (Fig. 88), then
at points remote from the short sides the surface of the corre-
sponding membrane can be assumed to be approximately cylindri-
cal; in view of the uniformity of the load p on the membrane,
the section of this surface by a plane perpendicular to the y axis
will be a parabola
}
C I b*
j \-^
v
Y ~ ~~ M
d ^ t
x
T-ryrrrl
WM/2\
Fig. 89
Fig. 90
(8.78)
17-1013
258 TORSION OF PRISttATICAL BARS AND BENDING [Ch. VIH
or (8.79)
M = 2C F - 2
t Q Cm Fm - 2 f f 3> ds> (8.80)
7 = M X = Q(lZ) Xi y
I z yy
I yx
y\ y n
U.
/L/2 j- j y\y
(8.81)
XX = Yy = Xy ~ 0.
The remaining stresses X Z X and Y z = Z y will be found
=Z from
the surface conditions so as to satisfy equations (I) and (VII) as
'
Fig- 91
well. Let us see what form these equations assume under our
assumptions (8.81). The equilibrium equations (I) will be
written as
^= 0;
dYz
=0; (8.82)
dZ K
Qx_
J (8.83)
17*
260 _
difference
loaded end B.
We now
in
TORSION OF PRISMATICAL BARS AND BENDING
stress
proceed
distribution
to equations
will
(VII),
be essential only near
observing first
[Ch
that
VIII
the
in
-- -
our case
O 7z _ <?(*-*) X.
y.
j
It can easily be seen that the first, second, third and sixth of
or
(8.86)
or
JL vV 2
co
V 0-
V>
-
v 2 cp
V ?
q(? _f J
dy /(I
- 2J >
ds dx ds
-
dy
d<pdx .
dy dy _ dy
'
dx ds '
dy ds ds
we get
dy
requiring that the stress function <p remain the same on the con-
tour of the section. It remains to consider the condition on the
end cross section. By projecting the forces applied to this section
on the axis Oy we have
262 TORSION OF PRISMATICAL BARS AND BENDING ICh. VIH
__ 2
4.
r
~
2
+C. (8.92)
'
y = m(r 2
x* y
2
) y.
_ 1 l+2a Q
x __ Q 3 + 2a /
2
*2
l-2a
"** 27 4(l + a) V
'** *
32 j y 11 '
(8<94)
Y = --*-.* }
It
TORSION OF PRISMATICAL BARS AND BENDING
~~
z 27 4(1 + a)
for o = -g-
we find from equations (8.94)
v
max X,-^ T r2 =
o
w -^r = T 7
- - Q7R
-== 11.375
T ,
/QQ^\
(8.95)
max*,, =
v 4Q
355-
=
0,
(9.1)
y z
~~ **
_-dz dx '
dxdy* dydz'
These formulas are obtained from one another by cyclic changes
between
(x, y, z); (cp,, cp 2 , cp 3 ).
We shall see later that there are other possible ways of prescrib-
ing the shearing stresses, but in any case their formulas should
not be altered in form in interchanging the letters in the nota-
tion of stresses, as required by the law of reciprocity of shearing
stresses, for instance
Y ~~
dxdy
Substituting the expressions of Xy and X z from (9.3) in the first
of equations (9.1), we obtain
dXx __~ d*<? 3 d2?2
">" '
dx dx dy* dz 2 dx
< 9 4>
-
M~
Z =-2H"
Functions (9.2) in this solution are called Maxwell's stress func-
tions.
Sec 56] SOLUT. OF DIFFEREN. EQUAT. OF EQUILIB. IN TERMS OF STRESSES 267
(Xt Yx Zx \
(9.5)
<h
= ffXydxdy; <ti
= f fy z dydz;
cp 2
= J J Xg dx dz.
dX y d^ 2 dXz
dx dy dz
y. *); <M* y. *) ( 9 6)
-
'
'
dz dx dy dz
dYz d3 <\> 2
== '
(9.8)
~T~~dz~ dxdydz
dYz
+ dy dx dy dz
dx dy dx dy dz dx
t
"+"
d^ 2
~ d^\
dz
'
\ dy )
whence we have
~ i
2
a
dz
/ah
( dx "^ dy
.
*h ah\
dz )'
d 1 / dfy[
' A^ry == ~'2"~dz'\Jx'^ ^4^2 ^4^3
dF
dydz dy'
dzdx '
z 2 dx ( dy
~l~
dz **
dx (9.9)
dxdy
'
f
'
JC
= ~'2'Jy\~dz"^
d ^^3 1
r ~dx'
I C^l
' '
2
dy dz
2
dx dy dz 9
d'F
dz 2 dx 2
' ^-a
Yz ~~
"dydzdx 2
'
(9.10)
Zz = (a + b) dx 2 dy f',
Zx = b
dzdxdy 2
Substituting these values of the components in the equations of
equilibrium (9.1), we see that they are identically satisfied. Let
us set up a sum of solutions of the form (9.10):
:=2X
m-l
dy
2 Xy
m =1
m dxdydz 2 '
y
1 v
-y
7i '
2 ^ '
(9.11)
dz 2 dx 2
r ^H^ dydzdx 2
m =1
^=^(am dy
2 '
__^ m =1
m dzdxdy 2 '
m =1 =1
(9.12)
1
Ju. A. Krutkov, The Stress-function Tensor and General Solutions in Stat-
ics of the Theory of Elasticity. Moscow, Izd. AN SSSR, 1949; V. I. Blokh>
Prikl. matem. mekhan., i960, Vol 14, No 4, pp. 415-422
i
270 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
+ Cm = am am = J (Pm + Tm + J;
(9.13)
and set
m= 1 m=
S 1
(9.H)
m=l
Zz \ R Z = ZT \
= er (9.15)
Sec. 57J EQUATIONS OF EQUILIBRIUM IN CYLINDRICAL CO-ORDINATES 271
dr
(9.16)
The first terms of each of
three
these equations can be formally ob-
tained from equations (9.1) by re-
placing the notations:
x by r, dx by dr\
y by 6; dy by rdft.
tions:
(9.17)
r
*-
<?r
(9.18)
v*
r <te r c>r dz
Z fl
= 4r c)0^ r2
1 (
A? r dr dz
r ^ 2 J
1 r)
_L _^_ ( f _ f) /i \ .
'
(9.19)
dr r f)^ v2 c)^ /
7 =
^ 1
^(r/ 3)
72--
dz*
whence we find
'
dr
~5T+ (te 2
'
^T ""^"'
2
(9.20)
'
r dr r2 dr
*
- 72"
~~~
dz '
--.
~~
r* dr
'
18-1013
274 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch. IX
where r= |/"jc
2
-f-
2
y -\-z
2
is the radius-vector from the origin to a
point (x, y, z).
Let us first show that function (9.23) is harmonic. Note first
that from the relation
r2 = +y jc
2 2 - 2
(9.24)
we obtain
dx
or
dr
~dx
~x
in the same way
9 26 '
77=7 -7- <
'
dx\r) dr ( r } dx r* r r3
1 . r3 _ Q-2 f^ r
r
-r + .
(9.27,
2
dy
(9.28)
= o
3 +3
-- _, o - =
which was to be proved. Hence function (9.23) is Laplace's po-
tential. The field of this function is the entire unbounded space
except for a single point the origin, where function (9.23) be-
comes infinite.
The gradient of function (9.23) has the following projections
on the co-ordinate axes [cf. formulas (9.27)]:
dx \ r )
~ JL
r3
'
dy
(9.29)
-L '
r3
18*
276 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
v
,
y), (9.30)
-^rj^-TrCosCr,
- =X T-
_X (9.33)
dy dy
,)
= 2 -&,
(9.34)
'
=== *"^
"7T~ ~H \~7T~ '
y] *
(9.37)
,
z)].
"278 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
or
= -/ cos (-,r). (9.38)
But
d
WO
"35~ VW ')~~dr
'
0> = A:-^4-v
dx ^^ y dy
-4-z
^^ dz
^r-^-.
dr
(9.39)
'
1
This follows from the definition of the derivative of a function of a point
-along a given direction.
Sec. 59] BIHARMONIC EQUATION 279*
Let us now show that if \|) 4 , \(?2, ^3 are some harmonic functions,
the functions
(9.43).
where cpi, q> 2 , cps, t|?i, t|?2, % are any harmonic functions. Instead
of these functions, one can, of course, take any of their partial
280 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch. IX
(9.44)
where cpi, cp 2 , cp 3
and \|)
are harmonic functions.
Let us now show that a function of the form
r^ (9.45)
dx
" dx Y-T'
2r = 2*, (9.46)
consequently,
-. T
2
r 2 -T- .
dx '
dx
Further,
Similarly, we get
(9.47>
' '
dx dy dz
(9.49)
dz
d<\>
a* _ '
dy dz
(9.50)
282 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch. IX
(VI)
(9.51)
V2 w + (k+} )*L = 0t
Z, =
(9.52)
Xy =
Thus, equations (9.51) and (9.52) are expressed in terms of two
elastic constants jn and k.
Turning now to Lame's equations (9.51), we recall that [Sec. 36,
equation (5.55)] in our case, when there are no body forces, the
dilatational strain satisfies Laplace's equation, i.e., is a har-
Sec. 60] REDUCTION OF BASIC EQUATIONS TO BIHARMONIC EQUATIONS 283-
monic function
V2 = 0.
The derivatives
<36 d6 ()Q
' '
dx dy dz
(9.53>
= 0.
(9.54>
(9.55)
(9.56)
We shall not use solutions of the type (9.50) here and shall only
indicate that they are needed for the problem of equilibrium of
an elastic sphere. We should note, however, the following: in-
stead of Lame's equations (9.51), we have obtained, by differen-
tiation, equations (9.53) of higher order but simpler form and
found the solutions of these equations (9.55) and (9.56). One
cannot say, however, that (9.55) and (9.56) will necessarily be
solutions of Lame's equations (9.51): the solutions of Lame's
equations will always satisfy equations (9.53), since the latter
are a consequence of Lame's equations; but because of the higher
order of equations (9.53), they will have a wider class of solu-
tions than equations (9.51). We
shall thus have to ascertain in
\vhat conditions (9.55) and (9.56) will satisfy equations (9.51); to
this end it is obviously necessary to substitute (9.55) and (9.56)
in (9.51) >.
0; W, VV,
We have
du _
~~
a ?l
dx dx dx 2 '
dv dy> 2
~dy~
dw ~" 3
~" Z ~* '
"
'
Further we obtain
(9.58)
or
The same computations for the second and third equations yield
<9 60 >
-
'
dy
(9.61)
where q>i, <p2, q>3 are arbitrary harmonic functions while the har-
monic function is determined from equation
\|> (9.60). In this
way one
\|>
can obviously obtain a great number of particular
solutions, there being available the arbitrariness of the choice of
harmonic functions cpj, (p 2 , 93.
1
Indeed, without violating the generality of equation (9.59), we can set
C = D(/z-f3), where D is a new arbitrary constant; we denote \|?
Dz=*tyi will
obviously be a harmonic function, just as ty; then (9.59) becomes
Sec. <51J BOUSSINESQ'S METHOD 287
but
d
dx
/>
\
^
dx}~
\
d < r2 )
dx
d*
dx~r
I
r2 ^
dx*
'
we find
(9.64)
(9.67)
since
1
See formula (9.39).
288 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
(9.68)
^ '
d*ty
dx 2
= dx
d I
\
a|/\
dx } dx '
T - \t
y
y y
dxdy dx \ dy /
'
JLfrJW
(*)
.
dxdz dx
dxdy dxdz dx dx
(9.69)
or
Hence we conclude
(k + 1) 6 -f 2<|. H- 2 (* + 3) $=
Sec. 61] BOUSSINESQ'S METHOD 289
-
(9.70)
r
"W "^ k +3 ^ =~ 2(6 + 3) '
or
or
rUV + (rU + *IS) l/= r
p0 .
(9.73)
1
This constant can be incorporated in either of functions 0o or t|)
in the
left-hand member of the equation.
19-1013
290 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch. IX
On the basis of (9.74) and (9.75) equation (9.73) takes the form
dV
or
the origin.
Z, = IMx, y); Z je
=^ 2 (x, y); Zy = O3 (.*, y). (9.77)
Since we are
solving problem the in terms of displacements
(Chapter IV, Sec. 24), we express the left-hand members of these
equalities in terms of displacements u, v, w by using Hooke's
law (9.52) and Cauchy's equations (III) of Sec. 24:
> --
(2
(9.78)
such that on the boundary 2=0 they turn into prescribed func-
tions (9.77), apart from a constant factor n which we introduce
19*
292 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch. IX
V, 0)
=* 2 (x f
y), (9.79)
y, Q)
= 3 (x, y).
*
3
_ 2 d^ *
dy
'
dz dz 1 dz 1 dz
(9.80)
k + l ~dz~~
dz (9.81)
dz "^
dy
= 3
.
dy
These equations, as already stated, hold on the boundary of the
half-space. Both their left-hand and right-hand members, how-
ever, are harmonic functions; it is known that if two harmonic
functions coincide on the boundary of a region, they are identical
in the whole region It follows from this that equations
l
.
(9.81)
are valid in the whole half-space 2<0. This conclusion enables
us, by eliminating the three functions cpi, cp 2 and cp 3 from (9.81)
and (9.60), to express ty in terms of coi, <o 2 co 3 which are assumed ,
_ "
k ~~ 1 ^ ^~ I / *?! I
dx "^ dy "^ dz
<*?2 I <??3 \ _
~~ d*l "^
<toi
do) 3
dx ~^ dy
|
+
'
k 1 dz* dz \ ) dz
In order to eliminate the functions <pi, 92, cpa from this we make
use of equation (9.60) and obtain finally
After this has been done, the procedure of solving the problem
is as follows:
1. According to the prescribed functions cp lt d> 2 <J> 3 [see equa- ,
placements u, v, w.
5. From Hooke's law (9.52) we find the stresses at an arbi-
trary point of the half-space.
Subsequently we shall confine ourselves to such cases of loading
on the boundary where the resultant is a finite quantity. Then,
as z increases beyond all bounds (i.e., for an infinitely removed
horizontal section of the half-space), all the displacements and
stresses must tend to zero; hence, on the basis of equations (9.61),
we note that as z-> oo the function \j? and its first derivatives
must likewise approach zero. We shall make use of this remark
in determining the function t]> from equation (9.82). By integrat-
ing it twice with respect to z from oo to an arbitrary value of 2,
we obtain
CO OO
294 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
/1=/2 = 0,
and we finally obtain
* T / oo
/(* +
oo
+
The first stage of the foregoing procedure of solving our prob-
lem makes it necessary to solve the Dirichlet boundary value prob-
lem for harmonic functions coi, o) 2 003 [see formulas (9.79)]. As ,
J
Zv dS = Zvm f dS = Z v
where Z vm is a certain average stress on the surface of the
hemisphere. Hence
and, finally,
(9.84)
Fig. 93
1
In addition, as pointed out above with respect to equations (9.81), it
must remain finite everywhere in the half-space, i.e., for 2<0, including the
case of x, y, 2 tending to infinity.
296 .
MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch. IX
(J)
i
= -7r.
Cz
* r /
oo
/***
oo
T
But in Sec. 58 we saw that [cf. (9.27)]
JL-J-
r
~~
3
dz
Therefore,
2 9^1 c
2 dty _
~ Cz . 2 C z _
~"
k+2 Cz
dz 1
k + 1 dz r3
"t"
* + 1 2 r2 AJ
+ 1 r3
'
but
Therefore, we have
foi _
~ ^3 ^_
~ C(l2a) A:
dz dx dx 2 r3
r = yV-K 2
, where 2 =
Sec. 63] EFFECT OF CONCENTRATED FORCE NORMAL TO BOUNDARY 297
( *L ~ A t tdt
+
~~~
*
~~
J r* ^J (*
2 2 8
5 ) *
2
+ 6
a
r(r *)
'
and, therefore,
(9.90)
In addition we shall
the expression of the volume
derive
expansion required seeking the
in stresses. For this purpose it
is easiest to make use of equation (9.80) by introducing therein
the value of the function \|> from (9.86):
= -r = ^r = ^ /t
c (l-2^r-
N /rk rnv
(9-91)
omit arbitrary constants in view of the fact that functions cpi, <p2 cpa must ,
vanish as z~>o.
298 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
K,
= Cp{-2gL-(l-2
(9.92)
Zr =
+(i-: r (/-)*
It now remains to C involved in the
determine the constant
expressions of stresses so as to make the load on the boundary (at
the origin) reduce to the given force P. For that it is sufficient
to require that the resultant of the normal forces Z z dS (Fig. 936)
over any horizontal section of the half-space at a constant depth
z= h be equal to P (P is the magnitude of the force which
is assumed to be compressive). Hence we obtain the condition
or
v2
Sec 63] EFFECT OF CONCENTRATED FORCE NORMAL TO BOUNDARY 299
p<*p _
~~
1
f <*(P
2
+ /*
2
) _
~~
1 1
3
'
1 2 J 5. 2 3/z
and hence
C|x-2ic = />; C=
-^p (9.94)
1
See E. Trefitz, Handbuch der Physik, Band VI. "Mechanik der elastischen
Korper", Berlin, 1928.
300 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
where
is the distance from the load (9.95) to the point of the bound-
ary (x, y) at which the displacement is being determined. The
expression of dw must be summed up with respect to the coordi-
nates and rj for all points of the loaded area; we obtain the
expression of the total vertical displacement of the point (x, y) of
the boundary:
(9.96)
^L ~ i
*'
dz
dx dy
(9.98)
df _ df dz _ df dz
dx dz dx dz dx
dx <)z dz
and similarly (9.99)
dy dz
dx dy
(9.100)
dx
^ l ~~~
dz
'
dy
~dx 2 ~l)z 2
dzdz d~z 2
d2f d2f
dz dz dz 2
t"79 J!
v / = T7
O J
dx
O
+ TT
2
I J
2
:===4
A
:Hr
dz dz
r
" J
-
/C\ 1
(9.101)
^ 1 \
dy
dzdz
and, therefore, we immediately find its general solution in the
form
, (9.102)
where Fi and F2 are arbitrary analytic functions of the variables
z and z (respectively). The correctness of solution (9.102) is
verified by differentiation.
Any harmonic function can thus be represented as the sum of
two functions of complex variables z and z, according to (9.102).
This fact underlies several methods of application of a complex
variable to the plane problem of the theory of elasticity. The most
significant contributions in this respect have been made by
G. V. Kolosoff and especially by N. I. Muskhelishvili who devel-
l
Complex Variable,
to the Plane Problem of Elasticity. Jurjev, 1909.
Sec. 65] FILON'S METHOD 303
(9.103)
where
0=^-4-4^-;
dx*dy
V2 ==-^
dx2+ -^-.
'
dy
2 (9.104)
^.IVT/
(9-105)
U=u ii).
(9.106)
}
_dU . dU
'
(9.107)
dz d7
dzdz
(9.108)
d2 U
dzdz
j!L
+ / Jli = 2 4i = 2 - (
dU
I
d[
L\
dx dy dz dz \ dz dz /
(9.109)
dx
/ =2 dz
=2 / 4-
dy dz \ dz dz ,
dz
Sec 65] FILON'S METHOD 305
Tg- 1* -i- o)
1
j l (Z) (ft -t- u/ 2 i*jj,
(9.110)
(9.111)
=i
2i
In this way we obtain
(9.112)
20-1013
306 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
%-(*+*) (9.114)
dv ,
da
= 7(Ff2) I
2 (k
_ _ _ (9-115)
1 )
[z/2
*
(z) 4- zf\ (z)}
- cp,'
(z)
- ?2 "
(z)}
,
1
I
|
(9.117)
where /i(2) and 91(2) are the conjugate expressions to f\(z) and 1
(9.118')
then
/2 (2) = /i (2)
=x (*, y) /t (je, y).
(9.118")
<p 2 (2)
=f , (2)
= xi (Jt, y) *<h (x, y).
Sec. 65] FILON'S METHOD 307
(9.119)
(9.120)
k +2 dx '
' '
dx dy dx dy
^_ ' '
(9.121)
dy dx dy dx"
20*
308 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
2lO,
z
---du
= dx
dv
TT .
dy
2.=4tffc-/;<*>]=
Thus,
(* + 2) 8 + 2/o>
2 =/;(*). (9.122)
where o>i (^) and 0)2(2) are arbitrary functions; it .is easy to verify
that the bracketed expression is a real quantity; in order that the
Sec. 65] FILON'S METHOD 309*
(9.125>
and if we set
_ _.
k+l
All our reasoning pertains to the case of plane strain; for the-
case of generalised plane stress Hooke's law takes the form
Hence
2
310 MORE GENERAL METHODS OF SOLVING ELASTICITY PROBLEMS [Ch IX
where
V = e xx + e yr
Therefore we have
Xx =
yy
i
where
, 2* 2a
*
1 a
Thus in the case of generalised plane stress all formulas Hold good
if the elastic constant k is replaced in them by k* '.
66 WAVE EQUATIONS
In Sections 26-28 of Chapter IV we established that Lame's
dynamic equations
(9.127)
Ox \~di*~
-.
dt* \dx>
the other two equalities of the same form are obtained replacing
x by y and by z; we get
(9.129)
or
2
V 2 0, (9.1290
-^Lr=a
where, as in (4.9), we designate
\~~dy ~dz)
\
p W r) /
\~dy ~dz)
\ m
'
&= -.
(9.132>
(9-133)
-*/?
The second case x
o).
= o) = o) =
1/ z occurs when the displacements
u, v, w have a potential (see the end of Sec. 10 of Chapter II)
fl
_ du ,
~T~
dv dw _
d**
"^ dz ~~ dx*
.
"*" "^
__ - 2<1>
~ v .
'
dx dy dy*
2
therefore,
1
See N. E. Kochin, Vector Analysis, Moscow, Objed. natichno-tekhn. izd.,
.1938, p. 223.
Sec. 67] SOME PARTICULAR SOLUTIONS OF WAVE EQUATION
(9.134)
w
d 2
j-
= a V w.9T79
2 2
where
dx + 4!( + ^/;
V 2/==
J |V 2 2 2 (9
v .136>
dy
'
dz '
dt 2
L_
-~ c 2
d *f
I)* 2
'
,
(
(9.140)
o-
hence
is
"where r=
_
a solution of equation
I/"* -f- y
2 2 2
(9.135), as
Consider an extremely important particular solution of the form
+
/ = /o(r,
z represents the radius-vector of an arbi-
trary point with respect to the origin. It is obvious that solution
(9.141) corresponds to a symmetrical distribution of factor /
with respect to the origin, so that the quantity / has the same
we
').
set out to prove.
'
(9-141)
x
~~
**
,
f)/
~~ ~~ df, r*-3 *
dx 2 dr 2 r2 '
dr r2 dr 2 r2 dr r
and z, we get
2
d /o
~~ _ 2
f x 2
.
df r 2
-x 2
dx 2 dr 2 r2 "^ dr r3
'
^ J
Ar22
dr
J
r^22
~ dr
I
I
'
<>f r*-y*
r3
~ dfdr r 2
-z 2
dz 2 dr 2 r2 r3
we obtain
/0 ~~
V 2f _^ /o 2
"^
.
2 d/ _
~ 1 ^
^r 2 r dr r dr 2 V^/o;-
or
r
~w '
r ~~
d< 2 a* 2
and the last equation takes the form of the simplest wave equa-
tion [cf. (9.137)]:
of a spherical wave:
(9-142)
f=F(x,y,z)-T(t). (9.143)
r<0 F(*,y,z)
Avhere W is an arbitrary constant. This leads to two differential
equations:
0, (9.144)
y, z)
= 0. (9.145)
1
See, for instance, L. S. Leibenzon, Theory of Elasticity, 2nd ed., Moscow.
Costekhizdat, 1947.
X
Bending of a Plate
68. GENERAL
A plate is a body of prismatical or cylindrical form whose
height h is small compared with the dimensions of the faces
(Fig. 94).
mark off the middle plane of the plate which divides
Let us
its height in half. This plane plays the same role in the theory
of plates as the axis and the neutral plane in the theory of bend-
ing of rods. Let us specify the system of co-ordinates, placing the
Fig. 94
(10.1)
Fig. 95
u = 2 sin a
,
2 tan a = dw
2-^-,
where a is the angle of inclination of the tangent to the trace of
the middle plane; we take the minus sign in the equality because,
for instance, in the case represented in Fig. 95 z>(\ tan a>0,
u<0. Thus
= -*-" (10-2)
xx
_ du _
dx 2 '
dx
'
~ dy
~ Z
dy
2 (10.4)
~ ~ d 2
w
*" '
dx dy dx dy
E ,
l Otfyy),
Yy =
j
E
4. a) 2(1
r,& Ez
dx* )
(10.5)
dxdy'
The remaining three components X z Yz Z
, , z can now be obtained
from the equations of equilibrium (I) which have not yet been
utilised; the first of them gives
dXz dX dX
Ez d
[ a2 Ox
where
dz
~ 1 s'
d
dy
2
w '
Sec. 69] BASIC EQUATIONS OF BENDING AND TORSION OF PLATE 321
and -
(10.6)
The loads and the supporting reactions of the plate are assumed
to be normal to its upper and lower planes and, therefore, on these
planes, i.e., for z = y there should be
dYz
dz dx dy
4
T-
~.
(10.10)
Here [see formula (6.5)]
Eh*
a'
_E^~' 110.11)
or for b=\ B=
In the expanded form equation (10.10) is written down as
+4
M dy=
l
f Xx zdydz = dy j
Xx zdz\ M^dx = dx f Y y dzz
_A A a
2 "2 f
~2
21
324 BENDING OF A PLATE [Ch. X
d*w (10.12)
M,=-D( dy*
In addition to the bending moments, the faces of the isolated ele-
ment are acted upon by shearing forces due to the presence of
Fig. 96 Fig. 97
t. ti
^ f Zx dydz\ N 2 dx = f Z y dxdz.
N l
= D^^w\ N = ~-D~V*w. 2 (10.13)
V
following forms:
d I dw \ d
dy \ Ox
where tan a* and tan a y are the angular coefficients of the tan-
gents to the aforementioned sections; thus the mixed derivative
is the rate of change of the slope of the tangent to the section when
the latter moves in the direction normal to its plane; in this sense
the mixed derivative characterises the "torsion" of the middle
surface of the plate. The shearing stresses Yx according to (10.5),
,
= J Y x zdydz,
H^-Dd)^. (10.14)
1
The moment considered positive if, when viewing from the side of the
is
positive x axis, we
see the rotation in the clockwise direction; the mixed deriv-
ative, which represents the "torsion" of a surface, is positive when the tan-
gent rotates counterclockwise [see above the second form of formula (a)].
This is taken into account in the sign of the right-hand member of formula
(10.14).
22-1013
326
to
Summarising
BENDING OF A PLATE _
the arguments set forth in this section, we come
the conclusion that the deformation of the plate under the
action of the load q(x, y) applied to its upper plane is deter-
[Ch. X
M =M =N =N =
l 2 l 2 0; H= const. (10.15)
-^ H
00-16)
This equation indicates that the middle surface of the plate under
pure torsion is a hyperbolic paraboloid (or the so-called oblique
plane). It has been ascertained above that the mixed derivative of
the deflection w represents the torsion of a surface which is equal
to the angle of twist per unit length designated in Chapter VII by
T. If we also take into account that
H
Sec. 70J
We
//=
, ANALYSIS OF RESULTS OBTAINED
by c, then
_ 327
T _ M t
'
ych*
M, = N = //=
2 0; Af ! = M = M *= const.
2 (10. 19)
a ^ ~~ ~D
~Jx* ~+~ ~Jy 2
~
'
d w
^ ddx*w ~ MD
2 2
,
- -
(3 '
2
dy
It follows from this that the middle plane of the plate under pure
bending is converted into a paraboloid of revolution. If the approx-
imate values of the curvatures were replaced by the exact ones in
22*
328 BENDING OF A PLATE [Ch X
equations (10.20),
d2w d2w 1
we would find directly from these equations that the middle plane
is bent into a spherical surface of radius
~~
M
The pure bending of a plate occurs only in exceptional cases, for
instance, in a circular plate loaded along its contour by uniformly
distributed bending moments of intensity M
per unit length of the
contour.
bending moment M^ \
twisting moment H, I
(10.22)
shearing force N 19 j
dx *
~
r o^
|
_ 329
Fig. 98 Fig. 99
for y = 0, w= 2 (10.23)
dy
In these two cases no conditions have been imposed on shearing
forces and twisting moments; but, being applied to the edge
resting on the support, they will have some effect only on the mag-
nitude of the supporting reactions of the plate and on the stress
distribution near this edge.
3. The right-hand edgeis free. This case presents the greatest
difficulties. No
restrictions are imposed here on the deflection;
however, all three forces MI, //, NI must be equal to zero and there
\0\
Fig. 100
arm dy. Having done this for all areas of the face concerned we
see that the forces applied at points and n on the boundary m
of the two areas reduce to a single force -jr-dy per length dy.
Sec. 72] ELLIPTIC PLATE CLAMPED AT EDGE 331
4^-
dy
(10.24
v
and two concentrated shearing forces H\ and //fat the ends, i.e.
at the corners of the plate; the contour is smooth and withoui
if
dx*
or (10.25
d3w
;
(io.26;
E. Reissner, J. Math, and Phys., 1944, Vol. 23, pp. 184-191; L. Bolle
1
23*
332 BENDING OF A PLATE [Ch. X
(10.27)
a
we have
(10.28)
^= -
+ {1-
1=0. (10.29)
dx dx d oy
dx*
= 8-rX
O C
+ 4-5-
9
2 I A C
dy*
= 84- i-^r*; (10.30)
dy*
c
dx* dx 1 dy* dy* IF
Substituting the values of the derivatives in equation (10.10'), we
find the constant c, which is not yet determined:
=:=
q
c
16
M ,Hi=H2,Ni
2 and N 2.
Calculate the maximum bending moments Mi and
3. 2 at the M
ends of the axes of the contour (moments of clamping) and at the
centre of the plate,
4. Setting b = a (circular plate), find the quantities indicated in
the problems above.
r sm
.
** sin
m-\ n=l
If we denote for the sake of brevity
A mn sln^- S \ n
^- = Umn , (10.32)
m=l /i=l
(10.33)
1
Indeed, on the contour, i.e., at *~0, x=a, #=0, yb we have w = Q,
__..-0, -3-5-
= 0, and the conditions for the simply supported edge will be
satisfied.
334 BENDING OF A PLATE [Ch. X
sin sin
or more compactly:
= q(x, ,
y),
.
/1A
(I0.34a)
, v
n= 1 /i = 1
where
(10.35)
/n=l n=l
(10.36)
sin -
sin
/sir
(10.37)
Sec. 731 RECTANGULAR PLATE. NAVIER'S SOLUTION 325
oo b
sl " *" y
j S Cmk
m =l
sin
^T~ =/?(* y) d y- (io.38)
. o
/.sm
s
r* 4
( fq(x,
2 +-\\'
b 2
J
b
f f q(X '
y} Sin ^ Sin
^ dX dy -
(1 '
39)
(10.40)
tively; for m and n even they vanish. On this basis, integral (10.40)
is different from zero only for m and n odd and then it is equal to
4qab
where m=
l, 3, 5, . . , oo; n=l, 3, 5, .
,
oo. . .
16
/f- m, ^
+
Sec. 73] RECTANGULAR PLATE. NAVIER'S SOLUTION 337
formula as
1
338 BENDING OF A PLATE [Ch. X
max
or q=. (10.43)
4/>sin^^-sin
jT\mn
m2
Sec 74] RECTANGULAR PLATE. LEVY'S SOLUTION 339
mite ,
rind
4P
n2 \
2 (10.44)
2
b )
Fig. 101
V2V2 ^ _ ? (* y)
(10.45)
Levy's method is valid for all cases where (Fig. 101) two op-
posite edges of a rectangular plate, for instance and AC, are OB
simply supported, while the other two, OA and BC, may be sup-
ported in any manner or be free. Bearing in mind these condi-
tions, we seek the solution of equation (10.45) in the form
(10.46)
m=l
340
Since
we have
_
d2w
gpr
=
BENDING OF A PLATE
rc
2
F i "* /m sm
= (M d z*_ n
=a\
V^
W~ a
o
2 ^
2
/
(y)
\
_ mnx
,
[Ch. X
(10.47)
or more compactly
mnx =q
V D*m (y) sin
,
'
,
(10.48)
ra=l
where
2
D^==/UV)(y)_ 2 (
J
where F m (y) are the coefficients in the series which are depend-
ent on y since the prescribed load q (x, y) is generally dependent
on it; substituting this in the right-hand member of (10.48), we
equate the coefficients of the respective terms on the left-hand
and right-hand sides:
account that
J sin - -
mnx .
sin
nnx
dx
,
A
j-n\r\
IV) / \
(y)-2(
n / H \
)
2
fn(y) +
* / \ i /
(
nn \
)
4
*
fn(y)
=
/ \
D-a
2
>
(10.50)
Upon integration in the right-hand member of this equation,
it will obviously become a function of y alone.
The general solution (10.50) will be composed of
of equation
the general solution of the homogeneous equation which coincides
with equation (6.82) of Sec. 44, if we set in the latter
(10.51)
Therefore, the general solution of equation (10.50) is
(10.52)
Jt
dy
s-
= 0. (10.54)
S*/(0)sm-^- =
//-v\ mix ^
0;
/IA nnx r\
(10.55)
( 10 .56)
7(y). (10.57)
f tt (Q)
= Ql /*(0) = 0. (10.58)
Sec. 74] RECTANGULAR PLATE. LEVY'S SOLUTION 343
A n = D n = Q. (10.59)
{t
<
CO sh
fflib
/-
~ftt / , \
(10.60)
2I::lsinh 2
a a
n 2 (rw) (10.62)
v '
(10.64)
e
mi
smh a
\
I a (x, 7])sm d
dx^dri. (I0.64a)
x
JJ I
1
344 BENDING OF A PLATE [Ch. X
y2 W_
~~
d2w
dx*
"" d w
2
.
_
~ d*w
dr 2
.
'7"
1 dw 1
. d2w
dr ~T~"7*' ~dP"
'
dy*
and the basic equation (10.10) becomes [cf. (IXp) in Sec. 47]
2 , 1 d 1 d2 \ / d2w 1 dw , 1 d*w
(10.65)
Here we shall deal only 'with the case where the load is in-
dependent of 6, i.e., when it is distributed similarly in all direc-
tions from the centre of the plate. In this case the deflected middle
surface must obviously be a surface of revolution, the deflection
w does not depend on and equation (10.65) is simplified [cf.
Sec
al
75]
__
equation (IX'p) in Sec.
I d*w .
do;
nn
345
< ia66
fifi
.
>
Cj In r w= +
C 2 r 2 In r C 3 C 4 r 2 (10.68) + + ;
w9TL
~ 64D '
(10.70)
"5y2"~7~5T
1
dx dy
~ *
(10.71)
24-1013
346 _ _
moments over
solution (10.70);
BENDING OF A PLATE
infinite
X
at
the centre of the plate, i.e., the deflection and curvature of the
plate will be equal to infinity. Hence we obtain
(10.72)
M = -~- [a
}
i
16 l
2
(1
v + a)' ;
(10.74)
(10.76)
Sec. 76] MEMBRANE ANALOGY. MARCUS'S METHOD 347
24*
348 BENDING OF A PLATE [Ch. X
DV 2w = 3K. (10.81)
The expressions of shearing forces (10.13) can be written down
^r
dx \dx* ~dy*')
d d*w d*w (10.82)
i \
--I-
J
or more compactly
V 2 * =--:. (10.84a)
= -|-, (10.83)
= jjr. (10.848)
Sec. 76] MEMBRANE ANALOGY. MARCUS'S METHOD
1
See H. Marcus, Die Theorie elastischer Gewebe und ihre Anwendung auf
die Berechnung biegsamer Flatten, Berlin, Springer, 1924.
XI
Variational Methods
of the Theory of Elasticity
~5T+~^ ]r
~dT + X = ?~WV
^yy~
dv'
dy
'
"**
,
:
===
-^r
dz
du'
+ dw'
dx
.
r
-
'
(11.4)
dy
(S)
where
'
dx
dX y dYy dZ y ' '
dw
_
'
dQ du dv
f ,
y Y 7
dw'
(11.8)
352 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch XI
'
+X e'
x xx +Y y
e'
yy +Z z
e'
zz +Y
(11.9)
We introduce the value (11.9) in (11.7) and replace / in the
left-hand member of the latter by its expression (11.5); we then
obtain
f 9 [(x
4S) (t)
(1UO)
This results from formal computations; nevertheless,
equality
attaching one or another meaning, on the one hand, to the state
of stress
X Y r Zz Y , Zx X y
,
(11. 11 a)
,p_^.
r=:
p^^.
=ip ~- = Oj and that the actual displacements
a, 0, w are taken as a', v', w'\ then equality (11.10) is transformed
into (5.62) and yields Clapeyron's theorem.
we
do
_
displacements, and we denote
LAGRANGE'S VARIATIONAL EQUATION
u!
imposed on these
= 8# ; i)
9
= 8-0; w f
l
Iw.
_ 35$
(S) CO
(11.13)
\-
dv' . du'
Ba =s tv = %w = 0.
and hence the conditions follow which are imposed on the variations at this
point:
354 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch. XI
(S) (t)
= *fV*, (11-14)
(-0
(S) (T)
t.
(11.15)
are not varied. Bearing this in mind, we put the variation sign
before the integral signs in the left-hand member of equation
(11.15) and write it down as
(S) (T)
where
(S) (T)
Sec. 78] LAGRANGE'S VARIATIONAL EQUATION 355-
Wdi = 0. (11.1Gb).
["-/
The quantity
variations should be equal to zero; if, at the same time, the body
forces are absent, the variation of the work done by external
forces is zero:
components.
However, in investigating the bending of bars, plates and
shells of small thickness the hypotheses of plane sections and
linear elements introduced there permit the elastic energy to be
calculated as the work of the bending and twisting moments
and shearing forces. For instance, in the bending of a straight
bar we isolate its element by two adjacent sections; then, neg-
lecting the work done by the shearing force, we have the energy
of bending of this element
1 JI/T dX
where l
's the curvature of the axis of the bar.
~~~T"'
The energy of bending per unit length of the bar (dx=l) is
2
w=v(l-} .
P dx*
we are in a position to express the elastic energy as a function
of the displacement v:
Jt. (11.20)
here
Pi
t
,
2
* r
ar e the curvatures of the middle surface,
1
= W
t 2 is the torsion.
d d
+20 )
forjc = ^(jc) = 0,
-^j-
= 0. (11.23)
for*-/
(11.25)
for*-/
Sec. 79] RITZ-TIMOSHENKO METHOD 359
(x)dx = ^-^cm ^; m = l, 3, 5,
^-2; c> -H
:
4
m
360 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch. XI
m
Hence
If we
retain only the first terms in formula (11.27), then, taking
into account that the error in an alternating series is less than
the first omitted term, we see that this error is less than
-^
= - = 0.0041 = 0.4%
and
J ~ __
^EJ
~ 382.5/
instead of the exact value
304% j-
F=
In the case of the bending of a rectangular plate having the
sides a and b long and simply supported along the contour, its
leflection can be expressed as
.
nny
1
,
., om
Here we satisfy not only geometrical but also statical conditions on the
1
a
-
b
--+- 4 \ a2 '
b2
\C mn
)
Sln
mr.l .
Sin
n^
-f n 2 !2 ) 2 ~~a~ ~T '
mrc; .
/ZTC7]
~"fl~
Sin
T" sin
. mnx
-
.
sin
nny
r
If we assumehere that m=
1, 2, 3, .,
oo and n= 1, 2, 3, oo, . . . . . ,
Dn* (a
2
62) 2
W max f = D^ ^ 2
_|_ 2)2
(11.34)
fied.
slants a m ,
fc m ,
c m are varied:
8w = 5ty m U, y, z)8c m .
(S) (T)
:
=0 f (11.36)
(S) CO
da m
(S) (t)
dV
(11.37)
(S) (T)
dV
Oc m
(S)
Ky X6 (11.38
OX Z
dx Oy ~0z~ ^" p
= '
dYx dK,
dx dy ^-+r p =o, (I)
dZ x
dx
4 '
<?y
^
-f
rfZ,
^ =
are inadequate for determining the stresses
*,, ry ,
z,, A V, K., z,. (11.40)
(11-41)
cos cos
25-1013
366 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch. XI
since we
are dealing with the equilibrium of a body.
We
take for the displacements ', v w' and the corresponding
f
',
\
>
(11.46)
(V)
c, y tii
Z : X y Y z Z x the variations
, , ,
^ V^ 5* I/" 5* *7
Oy\ v , O/ , 02-r,y
f (?>Xv u +W v + *Z w) = f (8*,^ + W
v v rfS y
e yy +
(S) <*)
,J</T. (11.48)
Sec
We may
80)
_
(Xx e xx
CASTIGLIANO'S VARIATIONAL EQUATION
e + Zzzz
+ Y yyy
y yy e xy + Y zyz
e + Xyxyz zz e +Z e y z yz
_ x zx ) d*
367
(T)
We
now introduce the expressions of strains in terms of
stresses (V) in the integrand on the right-hand side of (11.10);
this integral will then be expressed as
+ Z W + XJ
z y
= 8 (Xx Yy ), etc.
etc.
(t)
t r o n\. -i 1
(11.50)
(t)
25*
368 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch XI
f (8A> + 8K + 8Z ,w)</S =
t ,t> t
8
/ Wdt. (11.51)
(S) (T)
Consequently, .
{formula (8.16)]:
Xx = Y = Z> = Xy == Oj
y z
X = du Y = (11.52)
f f Udxdy. (11.53)
Let us assume that the lower base 2 = is fixed in its plane; then
u=v= and the corresponding integral becomes zero. It remains
to evaluate the integral over the upper base z = l where / is the
=
length of the rod; therefore, putting z l in (11.54), we write
down equation (11.51) as
(860 * dS.
(S) (S)
(11.55)
Consider the first of the integrals on the right-hand side:
(11.56)
y
fy(W)dy = ylU\- y.
flUdy.
y.
=- flUdy.
fy-jy(W)dy
Vi
2(1 + ) _ l
Sec. 81] APPLICATION OF CASTIGLIANO'S VARIATIONAL EQUATION 371
Fig. 102
-a (H.58,
But in the substitution outside the sign of the inner integral the
values of Xi and x2 refer to the points of the contour at which,
372 _ VARIAT10NAL METHODS OF THEORY OF ELASTICITY
J J
Similarly, we obtain
"*"
Substituting this in (11.58) and taking f*U out of the brackets,
we have
,(x,y)] 9
(11.60]
where Ci, C2 , ,
C n are arbitrary constants, fi(x, y), f2 (x, y),
. . . . . ,
. . .
-f- /IQI^I ~T~ ^02^2 I
0. (11.61)
of a
for
function of several independent variables in a
example, the first differential of a function
_ [Ch. XI
*(*, y) is
= ^J
I- *)[(*,
It will be noted that the expression in the braces represents
the elastic energy per unit volume produced by the stresses
82 / ( Wdt\>0. (11.62)
\(t) /
~~ -
J EF GFoi ~T~ GC
here C
the torsional rigidity of the bar; the energy is related
is
to unitlength of the axis of the bar and not to unit volume;
therefore, the variational equation (11.61) requires in this case a
minimum of the integral
N 2
Q
2
M
taken throughout the length 5 of the axis of the bar; this condition,
as is known, permits the solution of statically indeterminate prob-
lems. The general problem of the theory of elasticity for a body on
the surface of which there are given loads is also statically inde-
terminate; to solve it one may apply an approximate method
based on the variational equation (11.61) and analogous to the
Ritz-Timoshenko method presented in Sec. 79.
y
(
(II .66)
extremely
difficult, but it is not obligatory, as we have seen, in Castigliano's
method; at the same time the state of stress (11.65) subject to
conditions (11.66), as such, is impossible, since it would give
initial stresses ("self-stressed state"), whose absence is assumed
in advance (see Sec. 21) according to the hypothesis of the nat-
ural state. We
shall not, however, consider this state of stress by
itself but, using (11.64) and (11.65), we shall construct a new
state of stress
Yy Yy 4" C m Ky '>
ZX ZX (11.67)
:===
i"r($S) i
v/'"* *7( fft) y y(
Z
'
1
If stresses (11.64) satisfied the compatibility conditions, they would give
a complete solution of the problem.
Sec. 83] APPROXIMATE METHOD BASED ON VARIATIONAL EQUATION 377
the exact solution and the more closely we come to the satis-
faction of the compatibility conditions.
If we substitute the values of stresses (11.67) in the expression
of the elastic energy
4-2(l (1 1
.68)
.
-\-A mn C m C n -\-...
(11.69)
dW (11.70)
~dC^
:
+ A nmC
We note that from the first equation of (11.70) we have
therefore,
"mn "n (11.71)
378 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch XI
M=M +C M +C M +
((
1 1 2 2
... +CN M N = M V + 2Cm M m (
, (11.72)
where M() is the moment of a given load; MI, Af 2 ..., M N are the moments
,
The coefficients Ai$ in (11.69) and (11.70), which are integrals of squares
and pairwise products of functions involved in expressions (11.67), represent
a generalisation of known Mohr's integrals
MM' .
The analogy between (11.67) and (11.72) shows that the state of stress
(11.64), which satisfies the equations of equilibrium and the boundary condi-
tions but does not satisfy the compatibility conditions, corresponds, in the
case of the bar system, to the state of stress in the basic statically determinate
system of the force method (when the redundant constraints are removed in
the given system). The state of stress (11.65) is analogous to the "unit" states
of stress of the basic system which are due to the action of unit loads C m l. =
Consequently, these quantities may be interpreted as generalised forces in
formulas (11.67).
All the arguments of this section rest upon the condition that
the stresses on the surface of a body are not varied since it is
assumed that they are given. However, in the case where the
semi-inverse method is applied the stress distribution sometimes
is not given on certain parts of the surface, only the resultant
force vector (or the resultant) and the resultant moment vector
being given on these parts of the surface. For instance, in Chap-
ter VIII, where we considered the problems of torsion and bend-
ing of a prismatical bar, we prescribed on its bases: in bending,
the load Q subject to the condition that the moment of the
tangential tractions producing it be equal to zero; in torsion, the
twisting moment z M
subject to the condition that the resultant
force vector of the tangential tractions producing it be equal to
Sec. 84] LAME'S PROBLEM FOR ELASTIC RECTANGULAR PRISM 379
zero. The
stress distribution in all cross sections of the bar ap-
pears be
to identical; thus, varying the stresses throughout the
bar, we must admit their variation on its bases as well. In such
cases it becomes necessary to resort to the variational equation
of the general form (11.51) instead of (11.61). The next section
is concerned with the application of Castigliano's method to the
Fig. 103
)
= cos -cos = 0, I. 2, ...),
=0, 1, 2, ...),
(11.74)
= 0, 1, 2, ...)
2
They possess the following properties .
pp. 193-208.
Sec. 84] LAME'S PROBLEM FOR ELASTIC RECTANGULAR PRISM 381
the equality
*^:
Fig. 104
<p 3 and write down the expressions of the components in the form
of (9.4) of Sec. 56. The stress functions themselves will be repre-
sented as triple sums composed of functions (1 1.74):
(x) P n (y) Pp
m n p
<P 2
= 2 S 2 Bmnp Pm (x) P n (y) Pp (z), (11.75)
26-1013
382 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch. XI
+ C mnfP^(x)P,,(y)P p (z)},
(11.76)
=-sss
m n p
,
-
= 2SS (*).
2 = 0; z c
C ,
Cj, C2 , ..., CN _ l (11.78)
26*
384 _ VARIAT1ONAL METHODS OF THEORY OF ELASTICITY
The
application of Castigliano's method to Lame's problem
will be illustrated by the simplest example when the prism is
compressed or stretched by equal loads distributed arbitrarily
over two opposite faces, say, over the faces 2 = and z=c\ then
the boundary conditions on these faces will be
X = Y = 0; Z = F(x,
Z 2 z y) lorz =c
X_ 2 = r^ 2 = 0' J Z_ 2 = F(x, y) forz = 0;
See the
1
author's paper, Prikl. matem. i
niekhan., 1951, Vol. 15, No. 2,
pp. 137-148.
Fig. 105 shows the surface of distribution of stresses Z z at the
section passing through the middle of the altitude of the prism.
In spite of a sharp nonuniformity of
the load acting on the bases, the dis-
tribution of stresses at the mid-
section is sufficiently close to uni-
form distribution, in full accord with
Saint-Venant's principle.
Castigliano's variational method
made it possible to obtain the solu-
tion of Lame's problem for a prism
in othermore complicated cases of
loading. V. P. Netrebko treated thel
1
Vestnik Mosk. Univ., 1954, No. 12, pp. 15-26; 1956, No. 6, pp. 11-25.
2
In Investigations on the Theory of Structures, Moscow, GosstroyizdaU
1954, No. 6, pp. 455-468; 1957, No. 7, pp. 437-466.
,386 VARIATIONAL METHODS OF THEORY OF ELASTICITY [Ch. XI
dies can be carried out in a very general form through the use,
for instance, of formulas (9.18) which contain three stress func-
tions. Several problems of this kind were solved in the papers
of V. N. lonov, V. M. Ljubimov and E. R. Miroshnichenko 2 .
1
The author's paper. Prikl. matem. i
mekhan., 1953, Vol. 17, No. 4,
pp. 465-469.
2
V. N. Investigations on the Theory of Structures, Moscow,
lonov, in
Gosstrojizdat, pp. 413-436; E. R. Miroshnichenko, The Problem
1957,
7, No.
of a Cylinder Compressed Between Rigid Plates Without Sliding, Moscow,
Jviosk. lesotekhn. inst., 1957; V. M., Ljubimov. Prikl. matem. i mekhan., 1957,
Vol. 21, No. 4.
REFERENCES
5. A. Foppl and L. Foppl, Drang und Zwang, 2 Aufl, Berlin, Bd. 1, 1924,
Bd. 2, 1928
8. P. Pfeiffer, Handbuch der Physik. Band VI. Mechanik der elastischen Kor~
per, Berlin, 1928.
B G
Beltrami, E., 139
Galerkin, B. G., 301, 338, 344
Bernoulli, J., 9, 126
Galilei, G., 9
Beizeno, C. B., 229
Germain, Sophie, 9, 322
Blokh, V. L, 269 229
Grammel, R.,
Bolle, L, 331
Green, G., 10
Bredt, R., 256
Boussinesq, J., 211, 278, 284, 290
H
Hertz, H., 300
Castigliano, A., 273, 364 Hooke, R., 9, 73
Cauchy, A. L., 10, 29, 48, 235
Clapeyron, B. P. E., 10, 143
Coker, E. G., 188
Coulomb, C. A., 9, 116, 238
lonov, V. N., 386
D
D'Alembert, J. le R., Ill, 313, 354
Dirichlet, P. G. L., 198, 235 10
Jasinsky, F. S.,
Jastrzbski, N. F., 10
I N
Identities, Saint-Venant's, 55 Navier's equations, 21, 96, 110, 202,
Initial stress, 86 270
Integrals, Mohr's 378 Neumann's problem, 234
Invariants of Neutral axis, 127
strain deviator, 65 Nonhomogeneity, 73
strain tensor, 65 Number, Poisson's 74
stress deviator, 40
stress tensor, 34
Isotropy, 73, 89 O
Oblique planes, 326
Octahedral areas, 39
Operator, Laplacian, 74, 203-204
Kirchhoff's hypothesis of linear ele-
ments, 318
Stress, tensor, 28 U
invariants of, 34
Uniqueness of solution of
String, equation of transverse vibra- elasticity
tion equations, 136, 143
Unit elongation, 47
of, 111
Strip, in bending, 189
rectangular, 181
on absolutely rigid foundation, 189,
195 Variation of stress, 365
Superposition principle of, 72, 137 Variational equation, Castigliano's, 364
Symmetric matrix, 363 Lagrange's, 352
Velocity, of wave propagation, 106-
108
Vibrations, elastic, 19
harmonic, 104
Tensile test diagram, 73 longitudinal, of bar, 109
Tensor, basic, 380 uniform, 104
transverse of string, 111
correcting, 380
finite-strain, 67, 70 uniform, 104
of relative displacements, 56 Virtual displacements, principle of, 354
spherical, 37
strain, 62 W
stress, 29
Theorem, Bredt's, 256 Wave equation, 311
Clapeyron's, 143, 352 Wave length, 105
Euler's on homogeneous functions, Waves, equivoluminal, 313
87 Love, 316
Stokes's, 256 of dilatation, 312
Theorems of minimum elastic energy, irrotational, 313
355, 373 of distortion, 311
Torque, see Twisting moment, 116, Rayleigh, 316
236, 325 reflected, 316
Torsion, Coulomb's theory of, 116, 238 spherical, 314
of plates, 318, 326 standing, 316
of prismatical bars, 231, 368 velocity of, 106-108
of particular forms of section: Wedge, bending of, 219
circular, 71, 116 compression of, 219
elliptic, 239 loaded at vertex, 217
rectangular, 245 Work, done by elastic forces, 82
triangular, 243 by external loads, 143
of rod in form of body of revolu- by surface tractions, 145
tion, 273 least, principle of, 375
Torsion function, 232, 250
Torsional rigidity, 238
Trajectory of equal stresses, 214
Twisting moment, 116, 237, 325 Young's modulus, 74
To the Reader