Mathematical Theory of Elasticity - 1892 PDF
Mathematical Theory of Elasticity - 1892 PDF
Mathematical Theory of Elasticity - 1892 PDF
ON THE
MATHEMATICAL THEORY
OF
ELASTICITY
BY
A. E. H. LOVE, M.A.
FELLOW AND LECTURER OF BT JOHN'S COLLEGE, CAMBRIDGE
VOLUME I.
CAMBRIDGE:
AT THE UNIVERSITY PRESS.
1892
(.dll Right. resm~td.]
PBIIIITED BY C. 1. CLAY, )I,J., A'ND 8011111,
totally different methods from those given in the text. They are
intended as an acknowledgement of indebtedness, and a suggestion
to the reader for further work. I have tried to avoid as far as
possible reference either to erroneous mathematics or to incon-
clusive experiments.
I have not thought it advisable to introdu.ce collections of
examples for practice. On the other hand a number of results
are stated without proof. These are generally either of historical
interest or else of importance in the development of the subject,
but the analysis necessary to prove them would involve no point
but such as will be found in the text, or may fairly be assumed to
be known. The student without previous acquaintance with the
subject is advised in all cases to provide the required proofs. It
is hoped that he will not then fail to understand the subject for
lack of examples, nor waste his time in mere problem grinding.
In conclusion I have to express my thanks to Prof. A. G.
Greenhill, Prof. Karl Pearson, and Mr J. Larmor for their kindness
in reading the proof-sheets and for many valuable criticisms, to
Mr A. Harker for his kind assistance in the revision of the articles
on Crystallography, and to Mr C. Chree for his very careful
revision of the proofs, and for the many suggestions he has made
for the improvement of the work during its passage through
the press.
A. E. H. LOVE.
HISTORICAL INTRODUCTION.
PAOB
I
Scope of History. Galilei's enquiry. Enunciation of Hooke's Law.
Marriotte. Galilei's problem in the 18th century. Coulomb on
flexure, on torsion, on rupture by shear. Young on shear. Young's
modulns. Value of researches before 1820. History of analysis of
strain and stress. Methods of obtaining general equations. Navier.
Cauchy. Poisson. Saint-Venant. Non-molecular methods. Lame.
Green. Stokes. Clebsch. Rari-constancy. Arguments on both
sides of controversy. Voigt's experiments. Conclnsion against
rari-constant theory. Thomson's thermodynamical investigation.
Present position of Hooke's Law. Elastic constants of isotropic
solids. Neumann and Voigt on elastic crystals. Saint-Venant's
researches on reolotropy. Alteration of constants by change of
temperature. Static and kinetic moduluses. The thermo-elastic
equations. Curvilinear coordinates. Wave-motion. General theory
of vibrations. General problems. Solutions for the plane. Bons-
sinesq's theory of local perturbations. Solutions for the sphere.
Tidal effective rigidity of the earth. Wangerin's theory for surfaces
of revolution. Betti's method of integration. Vibrations of a sphere.
Saint-Venant's semi-inverse method. Torsion of Prisms. Flexure
of Prisms. Extensions of Saint-Venant's theory.
Arts. 21-49. 68
Effects of loading. Set. Elastic Limits. State of ease. Hooke's Law.
Proofs of the Law. Isotropy. .lEolotropy. Elastic constants.
Moduluses. Isotropic solids. Resistance to compression. Rigidity.
Young's modulus. PoiBBon's ratio. Equations of equilibrium and
small motion. Table of constants. 1Eolotropic solids. Green's
reduction of the number of constants. Cauchy's reduction. Sketch
of Crystallography. Neumann's aBBUmptions. Energy-function for
various crystal-systems. Moduluses and Poisson's ratios for reolo-
tropic solids. Voigt's experimental results. Amorphous bodies.
Curvilinear distributions of elasticity. Note on double-suffix
notations.
PAGE
Equations of wave-propagation in isotropic solid. Solution of the equa-
tions. Stokes's interpretation. Plane waves. Propa.ga.tion of dis-
turbance in ooolotropic medium. Equation for the velocity. Equa-
tions of a. ra.y. Wave-surface.
Theory of the free vibrations of solids. Principal modes. Frequency.
Normal functions. Frequency-equation has real positive roots.
Conjugate property of normal functions. Effects of sudden ap-
plication or reversal of load.
~- 82--116. 146
Saint-Venant's senn-mverse method. Imposed stress-conditions. Dif-
ferential equations. Expressions for transverse displacement. Ex-
pression for longitudinal displacement. Formation of the boundary-
condition. Simplified form of solution. Character of solution. Ex-
tension. Uniform flexure. Distortion of shape of cross-section by
flexure. The problem of Torsion. Torsion-couple. Sa.int-Vena.nt's
torsion-factor. Torsional Rigidity. Hydrodyna.mical analogy. In-
troduction of conjugate functions. Boundary-condition. Second
form of hydrodyna.mical analogy. Effect of cylindrical fl.a.ws on
strength of shaft. Strength under torsion. Maximum shear a.t
surface. Effects of corners. Torsion of circular cylinder. Elliptic
cylinder, distortion of cross-sections. Equilateral triangle. Rect-
angle. Sectors. Saint-Vena.nt's approximate formula.. 1Eolotropic
rectangular beam. Voigt's approximate formula. Non-uniform
flexure. Character of the stress. Principle of equivalence of statically
equipollent systems of load. Bending by transverse force. Flexural
Rigidity. Displacements produced by flexure. Asymmetric loading.
Strength under flexure. Strength under combined strain. Distortion
of cross-sections into curved surfaces. Circular Cylinder. Hollow
circular beam. Elliptic cylinder. lEolotropic rectangular beam.
~- 117--132. 199
Orthogonal surfaces. The line-element. Vector-differentiation. The
rotations of the system of normals. Dupin's Theorem. Expressions
for strain-components. Stress-equations. Energy-method. Strain-
equations for isotropic solid. Systems of surfaces. Polar formuloo.
Cylindrical coordinates. Radial strain. Sphere compressed by its
own gravitation. Internal heterogeneity of the earth. Spherical
shell under pressure. Radial vibrations of spherical shell. Problems
on cylinders. Rotating circular disc. lEolotropic shells under
pressure.
XIV CONTENTS.
NoTEs 345
A. On shear and shearing stress. B. On ooolotropic bodies. C. On
Betti's method of integration.
INDEX 349
CORRIGENDA.
1-2
4 HISTORICAL INTRODUCTION,
noticed by writers on elastic theory. The fact that a medium is possible in which
it does not hold good appears to have been first noticed in connexion with Electro-
dynamics.
4 Lam~ and Clapeyron, • M~oire sur l'~quilibre int6rieur des corps solides
homogenes'. Mem ....par divers savans, IV. 1833. The date of the memoir is at
least as early as 1828,
HISTORICAL INTRODUCTION. 'I
ponents of finite strain are practically those of Green 1, and
Saint-Venant 1, but the latter was the first to consider them
minutely. To Saint-Venant more than anyone else belongs the
credit of the adequate discussion of shear • ; he was the first
mathematician to call attention to its importance as a specific
kind of strain ; previously to his time the quantities we should
now call shears made their appearance simply as mathematical
expressions. Sir W. Thomson further simplified the discussion
of strain by the introduction of his strain-ellipsoid', and the
kinematical theory reaches its highest development in Thomson
and Tait's Natural Philosophy, Part I.
To a modern reader it might appear that the analysis of stress
and strain is a necessary preliminary to a general theory of
elasticity, but historically this was not the order in which dis-
coveries were made. The investigation of the general equations
by Navier does not depend on any such analysis; Poisson's
investigation involves an analysis of stress, but not of strain,
Green's an analysis of strain, but not of stress. There are in fact
three fundamental methods of arriving at these equations. The
first consists in assuming a law as to the character of intermolecular
force, and deducing the differential equations of displacement from
the equations of equilibrium of a single displaced "molecule".
This is Navier's method. The second method consists in forming
differential equations of equilibrium of any element in terms of
the stresses exerted upon it by the surrounding matter, and then,
by means of relations between stress and relative displacement,
eliminating the stress-components from these equations. The
required relations may be assumed, as in Cauchy's first investi-
gation, or deduced from experiment, as by Sir G. Stokes, or
calculated from an assumed law of intermolecular force, as by
Poisson and Cauchy. The third method consists in writing down
an expression for the energy of the strained solid, and deducing
for the summations with respect to angular space about the given
molecule, but not for the summations with respect to distance from
this molecule. The equations thus obtained are identical with
Navier's. The principle, on which summations may be replaced
by integrals, has been explained as follows by Cauchy 1 :-If the
molecular distribution be such that the number of molecules in
any volume, which contains a very large number of molecules, and
whose dimensions are at the same time small compared with the
radius of the sphere of sensible molecular activity, be proportional
to the volume, then, making abstraction of the molecules in the
immediate neighbourhood of the one considered, the actions of all
the others, contained in one of the small volumes referred to, will
be equivalent to a force through the centroid of this volume,
which will be proportional to the volume and to a function of the
distance of the particular molecule from the centroid of the
volume. The action of the remoter molecules is said to be
"regular", and the action of the nearer ones, "irregular"; and
thus Poisson assumed that the irregular action of the nearer
molecules may be neglected, in comparison with the action of the
remoter ones, which is regular. This is Sir G. Stokes's 1 descrip-
tion of Poisson's assumption, and it is the text upon which he
founds his criticism of Poisson. Without making this assumption
Cauchy arrived at Poisson's results.
Among later investigations of the stress-strain relations, as
given by the molecular hypothesis, we must note especially those
of Saint-Venant 3• In the first place he gave an ingenious proof
that, if the elastic reactions arise from changes in the molecular
configuration, and the intermolecular forces are functions of the
intermolecular distances, then, for very small strains, the stresses
must be linear functions of the strains. For in this case the term
of any stress, that arises from the force between two molecules, is
the difference of the amounts of this force in the strained and
unstrained states; and, since the force is supposed a continuous
function of the molecular configuration, this difference must be
ultimately a linear function of the variations in the intermolecular
1 In his memoir first quoted. Eurcices de MatMmatiques, 1828, pp. 241-243
of the new edition. Paris, 1890.
1 Math. aud Phys. Papers, vol. I. pp. 116 sq. and Oamb. Phil. Soc. Tram. VIII.
1845.
8 See his ediiion of Moigno's Btatiqu.e, and of Navier's Lei}OflB, also the memoir
1 • Snr lee equations qui expriment lee conditions d'equilibre ou lee lois du
mouvement iut&ienr d'un corps solide '. Exercices, 1828. This memoir precedes
those iu which the same author made use of the molecular method.
s 'On the Equilibrium of Elastic Solids'. Edinburgh R. 8. Tram, xx. 1853.
12 HISTORICAL INTRODUCTION,
Scientiftc Paper1.
22 HISTORICAL INTRODUCTION.
1 Vorle.ungen 11btr die Theorie der Elalticitlit dtr festtn Karper und de1 Licht-
a Proc. Lcmd. Math. Soc. IV. 1873, and Theory of Sound, vol. I. 1877.
HISTORICAL INTRODUCTION. 27
1 '
Ricerche intorno all' equilibria de' corpi elastici isotropi '. Beak Accademia
dei Lincei, Rome, 1832.
HISTORICAL INTRODUCTION. 29
gravitation, and to any distribution of surface-traction. Lame'
commenced by transforming the equations to polar coordinates.
The equations of the problem in rectangular coordinates were first
solved by Sir W. Thomson •, who applied the results to the
consideration of astronomical problems relating to the elastic
equilibrium of the earth, deformed by tide-generating forces, or
centrifugal force. In the case of the tides it was shewn that the
degree of rigidity to be attributed to the solid, in order that ocean
tides upon it may be similar to those on the earth, is very con-
siderable, and the result discredits somewhat the geological
hypothesis of the internal fluidity of the earth. The application
of the problem to test this hypothesis is however beset with
difficulties which have not yet been surmounted. The spherical
harmonic solutions of the equations of elasticity have an extended
application to other problems besides that of the equilibrium of
the sphere. They are solutions in terms of integral powers of the
coordinates, and they have been considered in this light by Mr
Chree 1, who has shewn, by means of them, how to obtain a solution
of some problems relating to the equilibrium of ellipsoids, and has
also utilised them to verify Saint-Venant's solutions for the torsion
and flexure of beams. Another application of them which has been
recently made 4 is to investigate the effect of flaws in diminishing
the strength of structures, verifying for the simplest case the
factor of safety 2, allowed by engineers to guard against this form of
weakness. A different solution of Lame's problem has been given
by Borchardt 6• Instead of spherical harmonic series the displace-
ments are expressed in terms of definite integrals involving the
given surface-tractions, and a like solution has been given, by the
same writer, of the problem of the strain in a sphere deformed by
unequal heating 8, setting out from the thermo-elastic equations of
1 ' Ueber das Problem des Gleichgewichts elastischer Rotationskiirper ', Grunert's
ANALYSIS OF STRAIN 1•
Fig. 1.
1
By " extension " of a. line we sha.ll always mean the ratio of the increment of
lengih to the original length. Contraction will be treated as negative extension.
38 ANALYSIS OF STRAIN. [4
In the figure let B be the new position of A, and suppose AB
so chosen that the middle point of AB is on the axis of y, then
the amount of the shear is 2 tan !AOB. The angle !AOB is
called the angle of the shear, and its tangent is half the amount of
the shear.
4. Componenta of Strain.
The deformation of the figure will be completely known when
we know the new length of every line in it. Since parallel lines
are equally extended, we only need to know the new lengths of
lines drawn through the origin.
Let l, m, n be the direction-cosines of a line drawn through
the origin, and r its length, fll, y, z the coordinates of its other
extremity, so that
fll=lr, y=mr, z=nr.
After strain let the point (fll, y, z) come to (a;.,'!/~> z1). Then by
writing lr, mr, nr for fll, y, z in (1 ), squaring and adding we find
the new length r 1 of the line given by the equation
r;J = r2 [1 + 2 (lse1 + m 2es + n 2e3 + mns1 + nlss + lms8)] ••• (2),
where e1 =au+ i (au2 + an2 + aa12) } . . . . . . . . . . . . ( ),
3
&1 = ~ + aw + (~2~ + assa~ + a:JP.ss)
and Es, e3, Bs, 83 are to be found from these by cyclical interchanges
of the suffixes 1, 2, 3.
The deformation is thus completely determined by means of
the six quantities e~> e20 e3, s~> &2 , s3. We shall call these the
components of strain.
The meaning of the quantities e1 , Es, e3 is at once apparent, for
the extensions of lines parallel to the axes are
-/(1 + 2e1)- 1, -/(1 + 2es)- 1, -/(1 + 2ea) - 1.
To see the meaning of ~. s2> s3 , it is convenient to form an
expression for the cosine of the angle between the strained
positions of two lines through the origin.
Let (l, m, n), (l', m', n') be the primitive directions of the lines,
then the cosine of the angle between their strained positions will
be found by taking two points (fll, y, z), (fll', y', z') one on each line,
and supposing their strained positions to be (a;., Y~> z1), (a;.', y1', z1'),
the cosine of the angle between them after strain is
a;.a;.' + Yl'!// + z1z1'
5] COMPONENTS OF STRAIN. 39
Now ~. Y~> z1 are given by (1) in terms of x, y, z, and x1', y/, z/
are the same functions of a!, y', z' ; also x : y : z = l : m : n, and
al : y' : z' = l' : m' : n'. Hence we find, for the cosine in question,
the value
ll'(l +2•1)+mm'(l +2~2)+nn'(l +2• 3) +(mn' +m'n)11 +(nl' +n'l)s2 +(lm' +l'm)s3
J[l +2 (•u ••• •s• tsu tss• l-'Jlmn)Z] J[l +2 (•1• ••• •s• tsu tsu ¥Jl'm'n)2J
..................... (4).
In particular if the lines be the axes of x and y this reduces to
1
v(l + 2e1} .V(l + 2e2) ••• ••• • . . . . •••• • . . . . . (
5).
Thus s~> s2 , s3 depend upon the angles between the strained
positions of the lines initially coinciding with the axes.
Another way of looking at this matter is to suppose that the
strain consists of a simple shear, say of the planes y parallel to the
axiS X.
z' z. 'Ins
I '11'8
and let the components of the strain referred to the new axes be
e1', e2', e1', s1', s2', s8'. Then, since the new and old lengths of a
given line are independent of the system of axes, it follows that
the strain-quadric (7) will be transformed to
(e1', e;, e1', !-s1', !s1', -fss'ia;yz)1= k,
7] TRANSFORMATION OF STRAIN. 41
and thus e/ ... will be the coefficients in the transformed equation
of the quadric (7). We thus obtain the equations
•t' =•tlt1 +•s7nt1 +•an,.2 + 8 tfntn,. +•sn,.Zt +•s~fnt, }
•t'=2•tl.Js+2.sm,nts+2•sm2~+st(mzns+mzn2)+•2(n.Js+~ls)+•s(4ms+W
............ (11),
and the other components can be written down by cyclical inter-
changes of the suffixes of the l's, m's, and n's.
We remark that by a well-known theorem of Solid Geometry
the quantities
E1+ e, +e., &12+ Bs2+ s,•- 4 ( e.ea + eaet + e1e,),
4f!tEsf!a + S1Ss83 - EtSl- Es&12 - E,S11
are unaltered by the transformation of coordinates. These are
called invariants of the strain, they are the coefficients in the
discriminating cubic of the strain-quadric (7).
7. Examples. Eztension and Shear.
We may utilise the properties of the strain-quadric to discUBB
the components of strain in particular cases. Of this we shall give
two examples.
(1) Suppose the strain (e1 , e2 , e., s1 , &2 , s1) equivalent to a
simple extension. Let e be the amount of the extension, and
v'(1 + 2E1) -1 = e, then referred to its principal axes the strain-
quadric is
whence
Fig. 2.
In the figure let P be a point on one principal axis before
strain, and Q its strained position, ON the perpendicular on PQ.
Suppose ON= 1, then PQ = s = 2 tan a, and tan 8 = (1 + sin a) sec a.
Hence
NQ=cot 8 +s =cosa/(1 +sin a)+ 2 tan a= tan a+ sec a,
1 +sin a cos a
and angle POQ = tan-1 - - - - - tan-1 . =a.
cosa 1 +sma
Thus the principal axes of the shear are rotated through an
angle equal to the angle of the shear. It is clear that after the
strain the figure can be turned back through this angle without
any alteration of the length of any line in it, and the simple shear
combined with this rotation is called a pure shear.
In pure shear lines parallel to one of the principal axes of the
shear are extended, and lines parallel to the other principal axis
are contracted, and since the principal extensions are connected
by the relation
(1 + '1/1) (1 + 'IJt) = 1
given in equation (19), we get the following representation of pure
shear:
44 ANALYSIS OF STRAIN. [7
Let .ABCD be a rhombus, whose diagonals are in the ratio
1 + '1/1 : 1 + '1/2> and are in the direction of the principal axes of the
D' Fig. a.
shear, and let .A'B'C'D' be an equal rhombus, with its correspond-
ing diagonals at right angles to those of .A BCD; then by the pure
shear, consisting of contraction along .AC and extension along BD,
the first rhombus will be transformed into the second
The reader should find no difficulty in verifying the following
methods of producing any homogeneous strain :
(1) Any such strain can be produced in a figure by a shear
parallel to one axis of planes perpendicular to another, a
uniform extension perpendicular to the plane of the two axes, a
uniform extension of all lines of the figure, and a rotation.
(2) Any such strain can be produced by three shears each
of which is a shear parallel to one axis of planes perpendicular to
another, a uniform extension of all lines of the figure, and a
rotation.
8.Pure Strain.
In general a strain is said to be pure when the principal
axes of the strain-ellipsoid are lines which retain their primitive
8] PURE STRAIN. 45
X y z
E ll ~ ~
'1] l2 ms '712
--
~ l, m, na
and let Qh, ylJ z1 be the coordinates, after strain, of the point,
whose coordinates, before strain, are x, y, z; then
"h = ~E1 + l,'T]1 + la~~
= ~a.E + l,/3'1] + la'Y~
=~~+~+~+~~+~+~+~~+~+~
The coefficient of yin the expression for x1, i.e. the coefficient a 1s. is
a~~+ f3l2mz+ rylsm,,
and we should find the same value for ~· We should find in like
manner ass= as:a. an = au.
Thus, if the strain be pure, we have the relations
~ = as:a, am= au, au= an .................. (20).
9. The Elongation-Quadric.
The quadric (23) is called the Elongation-Quadric. Let P be
any point (ro, y, z) of the figure, which is transformed to P 1 (ro1, y1, z1)
by the strain, then, if we define the elongation of OP in direction
OP to be the projection of P 1P on OP, this is
e=-
au f_av OW
OX' -ay' g= oz'
Ow Ov au aw Ov ou
a= oy +az· b= oz +OX' c=-+-
ax ay' (36),
OW ov au -ow Ot1 au
2wl = oy - az' 2w2 = - -
oz ox' 2w1 = - - -
ax ay
we find that the component displacements of Q, when P is
regarded as held fixed, are
rotation &c. hold for the matter about any point. For convenience
we state here that
Ou OV OW
.:1 =ox+ oy + oz ........................(as)
is the cubical dilatation of the matter about the point (x, y, z).
If the strain be pure, 11T1 =0, 11T2 =0, 1D'1 =0, or we have
Ow OV ou Ow
oy = az · oz = ox •
so that there exists a function cf>, such that
ocf> ocf> ocf>
u = ox • v = oy • w = az ·
The function cf> is called the " displacement-potential." Its exist-
ence is confined to the case where the strain is pure.
If the displacements be finite, the deformation of the body
in the neighbourhood of any point P can still be expressed
by six components of strain. Let r be the unstrained length
of any short line PQ of the body, r 1 its length after strain, and
l, m, n the direction-cosines of PQ before the strain, then, as in
art. 4,
(r12 - r2)/r2 = 2 (l2e1 + m2e1 + n2e3 + mns1 + nls2 + lms3) ... (39),
where
ANALYSIS OF STRESS 1•
JJJ(pX +?~z +
0
~11 + O!z) rkdydz=O .......... (3).
Since this equation holds whatever the surface S may be, provided
only that it is closed and altogether within the solid, it holds when
S is reduced to a point, and thus we have the differential equation
of equilibrium
pX + oXz + oX11 + oXz=O ............... (4 ).
ox oy oz
In like manner, resolving parallel to y and z, we obtain the
equations
pY + ?~+ o:f11 + oYz= 0}
ox oy oz
oZz oZ oZz ............... (5).
11
pZ+a;-+ o9 + oz =O
In order that the part of the solid within S may be in equili-
brium the sum of the moments of all the forces applied to this part
~bout any axis must vanish. By taking moments about the axis
of x we get the equation
fffp (yZ -zY) rkdydz
+ JJ[y(lZz+mZ11 +nZz)-z(lYz+mY11 +nYz)] dS =0 .•. (6).
The surface-integral becomes by Green's transformation
1
For the symmetry see Note on Double Suffix Notations at the end of chapter m.
60 ANALYSIS OF STRESS. [14
Then the equations of equilibrium become
aP au aT
ax + ay + oz + pX = 0
au+ oQ +as +pY=O
ox Oy oz .................. (10).
aT+ as+ oR + pZ = 0
ax oy oz
In case the solid is in a state of small internal relative motion we
have to put instead of pX, p (X-~), and similar expressions in
that the moments of the bodily forces that act upon any element,
about some set of axes whose origin is in the element, are
ultimately the products of the mass of the element and small
quantities which vanish when the linear dimensions of the element
vanish, in other words that the bodily forces acting upon any
element are ultimately reducible to a single force in a line that
meets the element. This is true for such bodily forces as gravita-
tion, and we shall limit our enquiry to the cases for which it holds
good. In the theory of Magnetism and some other parts of
mathematical Physics it has to be supposed that the bodily forces
applied to an element of the medium reduce to a force through
the centre of the element and a couple 1 • When this is the case
we may take L, M, N for the components of the couple, and then
to equation (6) we have to add the term JJJLdxdydz, so that instead
of (7) and (8) we get
Zy- Yz +
Xz-Z.:+M=O .
L=Ol
Y.:-Xy+N=O
In the theory of Elasticity these considerations are unimportant,
and the equations (7) and (8) are true for the systems of forces
usually considered.
j
y z
·--- - · ---~
x'
I I l I '11~1 •••••••••••• (13).
y
z'
I '2
~,--11
I la
__
I mt
'111-:J
l1 U+~Q+~S.
hence
P'=l1 (~P+m1 U +~T)+m1 (l1 U + ~Q+~S)+~(l1 T+~S+~R),
and
U' = l2 (l1P +m1 U +n1 T)+1~(l1 U +~Q+~S)+1!..!(l]T + m1S + ~R);
or
P' =l12P+~'Q+nt'R+2~n 1 S+2n1~T+2l1n~ U, }
U'=~lsP+ntr~Q+nr_n~+(~ns+ms~)S+(~l~~)T+(~m,+l2~)U.
•..... (14),
and from these Q', R', and S', T' may be written down by cyclical
interchanges of the suffixes 1, 2, 3.
It is to be observed that, if the quadric
(P, Q, R, S, T, U"§..xy•)l=k ............... (15)
be transformed by the substitutions given by the scheme (13), it
becomes
(P', Q', R', S', T, U'"§..:rly'z')l=k.
This quadric is termed the Stress-Quadric, and it appears from the
above equations that the normal stress across any plane is inversely
as the square of the radius vector of the stress-quadric which is
normal to the plane.
Again the transformation may be such as to refer the quadric
to its principal planes. When this is so S', T', U' will vanish, and
the stress across any one of the coordinate planes is normal to the
plane. The stresses P', Q', R' are then called principal stresses,
and we learn that there are at any point of the solid three co-ortho-
gonal planes across each of which the stress is purely normal
These are called the principal planes of stress.
Fig. 4.
P=aw
oe' Q=aw
of' R=ow}
og
ow oW ow ············ (20),
8= oa, T=NJ, U=ac
and W is the potential energy of the strained elastic solid in the
condition denoted by (e,f, g, tl, b, c). We shall shew in chapter V.
that W exists when the solid is strained at constant temperature,
or is executing small vibrations.
5-2
CHAPTER III.
ELASTIC PROPERTIES OF MATTER-STRESB-STRAIN RELATIONS 1•
x'
y'
(X+ p.) oy
oA + p.V~ + pY = •p asv
ot •••••••••••• (12),
2
F= ZAa+2p.(:.+m1ll'a-n1ll's).
isotropic, but they are treated as such by the authorities referred to.
2 'Units and Physical Constants'.
78 STRESS-STRAIN RELATIONS. [31
that in some cases we give E, in others k, in others p., and in
others rr. The nnmbers in square brackets occasionally given in
some columns are obtained from those given in other columns by
application of known relations between elastic constants.
chapter.
32] CRYSTALLOGRAPHY. 79
Again we shall prove that, for a system of discrete particles
homogeneously arranged, whose action upon each other is such
that the force between any two is in the line joining them and
is a function of their distance, there is a further reduction in the
number of constants from 21 to 15, effected by the equations
c23 = c", c31 = C111, ~~ = c88 }
............... (20).
c1, = c1111, ~=eM, Cas= c415
This is Cauchy's reduction of the number of constants, and there
are many practical applications for which it is convenient to
suppose it to hold good. We shall not however expressly introduce
it into the general theory, as we have no sufficient ground for
supposing that it expresses a necessary property of solid bodies,
and it has not been verified by experiment.
32. Natural Orystall.
Further reductions in the number of constants take place, when
the solid exhibits any kind of structural symmetry. The theory of
the possible symmetries, and of the forms of bodies possessing these
symmetries, is the science of Crystallography 1• The internal
structure of crystalline bodies can be inferred from the forms in
which they crystallize. These forms are always bounded by
surfaaes which are very nearly plane, and deviations from
plane forms are treated as exceptions, crystals bounded by
planes being regarded as the rule. If any three edges of a
crystal be chosen as axes, the planes of the crystal may be re-
ferred to them, and any plane can be determined by its intercepts
on the axes. The law of crystal form experimentally ascertained
is that, for any crystal, these intercepts are rational numerical
multiples of three fixed lengths dependent on the form. The
ratios of these fixed lengths are in general irrational, and are
called axial ratios. Crystals are classified accoPding to their
symmetry, and their axial ratios, when referred to the most
convenient system of axes. Thus if a, b, c be three fixed inter-
cepts on three fixed axes, any crystal must be bounded by planes
parallel to the planes whose intercepts are afh, bjk, cfl, where h, k, l
are positive or negative integral numbers. The ratios a : b : c are
the axial ratios, and depend on the material. The symbol (h, k, l)
1 The system of notation adopted is that of Miller, see e.g., G. H. Williams,
Elements of Ory•tallography. Macmillan, London, 1890.
80 STRESS-STRAIN RELATIONS. [32
represents a family of parallel planes. If the crystal possess
symmetry with respect to the plane of two of the axes, (say the
b-axis and the c-axis), then the existence of a plane face, forming
one mem her of the family (h, k, l) requires the existence of a face
forming one member of the family (- h, k, l). The collection of all
the planes required by this law forms a complete or holohedral
crystal form. Of equal importance are the partial crystal forms
arrived at by the selection of certain planes from those of
any complete crystal form. If half the planes be selected the
resulting form is said to be hemihedral, if one quarter tetartohedral;
the half or quarter selected must however be chosen according to
certain rules, depending on the symmetry of the crystal. If, when
the axes are suitably chosen, any one of the axial ratios become
rational, it is clear that this ratio may be taken to be unity, and
the two axes concerned are said to be equivalent ; if further these
axes be normal to planes of symmetry, they are said to be equi-
valent axes of symmetry. The law of selection of planes to make a
hemihedral form is that only such planes can occur as intersect
equivalent axes of symmetry at the same distance from the origin,
at the same inclination, and in equal numbers 1• The selection of
half the planes of a complete crystal form may either include or
exclude pairs of parallel planes ; in the former case the resulting
form is said to be hemihedral with parallel faces, in the latter
hemihedral. with inclined faces.
In the theory of elastic crystals, it is convenient to introduce
two sets of rectangular axes. The axes of (a:, y, z) are perfectly
general, and the axes of (a:0 , y0 , z0) are parallel to lines to which it
is convenient to refer the faces of the crystal, (sometimes, but not
always, crystallographic axes). We shall denote the displacements,
stresses, and strains, referred to the latter system, by (u0 , Vo. w0),
(Po. Qo, Ro. So, To, Uo). and (eo, fo, Uo. ao. bo. Co); and the most
general system of elastic constants corresponding to (16), when
referred to the axes of (a:0 , y 0 , z0), will be denoted by a's with double
suffixes instead of c's.
1 An example will make this clearer. If no two axes of symmetry be equivalent,
'
'\
\\
\\
\
\
\\\
\
\
---------;:;1
/
, ... ""
/
:
I
I
/ I
,,' I
Fig. 5.
\\
\\
\
\\
\,\8'
I<~"'---------+---71L'·_/_/_/_~-~-~;~T···............ A
d
Fig. 6
6-2
84 STRESS-STRAIN RELATIONS. [36
planes of the form, and the complete crystal form is the octa-
hedron with rectangular diagonals of different lengths.
The types of the possible hemihedral forms are :
1°, The tetrahedron formed by the planes ABO, A'BO', A'B'O,
.A.B'O', where AA', BB' and 00' are the axes. This has no plane
of symmetry but the figure will come into a similar position with
respect to the axes after a rotation through two right angles about
either axis, hence for these inclined-faced hemihedral forms formula
(21) holds.
2°. The half-form whose planes are ABO, ABO', A'B'O,
A'B'O'. This is identical with the complete monoclinic form,
and may therefore be rejected.
3°. The half-form whose planes are ABO, A'BO, A'B'O, AB'O.
This has two planes of symmetry, and formula (21) holds.
We shall write (21) in the form
2W=(A, B, 0, F, G, H'§.eo,fo,go'f+Lao2 +Mbo'+Nco2 ···(21).
In this notation Cauchy's relations are
L=F, M=G, N=H.
36. Quadratic or Tetragonal System-(6 ConstantA).
This system has three rectangular planes of symmetry, and two
of the axes are equivalent; let these be the axes of :c0 and y0 , then
Po must be the same function of e0 that Q0 is ofj0 , and Po and Qo
must have the same term in g 0 ; also 8 0 must be the same function
of a 0 that T0 is of b0 ; we thus get the equalities au=~, au=~.
a"= a156 , and W is given by the equation
2W=(au, au, aas, a..~. ~'J'§.eo,fo,goJ +a" (ao9 + bo2) + aeaco2 . . . (22).
Let (h, Tc, l) denote any plane of the complete form. Then the
complete form will also contain the planes (± h, ± Tc, ± l) and
(±Tc, ±h, ±l). If Tc=h the figure is an octahedron with rectangular
diagonals two of which are equal in length. The hemihedral forms
derivable from the most general complete form are obtained as
follows:-
1 °, Selection by alternate planes. Of the 8 bounding planes,
4 are obtained from the form (h, Tc, l) by taking the signs all
positive, or two negative and one positive, and the other 4 are
obtained from these by interchanging h and Tc and at the same time
36] SYSTEMS OF CRYSTALS. 85
changing the sign of l. The resulting inclined-faced hemihedra.l
form has no plane of symmetry, but by rotating it round the z0
axis through 90° it will come into a similar position with respect
to the axes. We shall see presently that equation (22) does not
hold for this case (fig. 7).
Z'
Fig. 7.
In fig. 7 .AB and .A.'D are two sides of a rhombus, two of whose
corners are .A., .A.', and the other two are on oa equidistant from 0.
Ba and an are two sides of an equal rhombus, one of whose corners
is a, another on Oa opposite to a, and the other two on .A..A.'. The
complete form consists of two pyramids, vertices Z and Z' standing
on the 8 sided figure .A.BaD.A.'... The hemihedralforms obtained
by the first method of selection are such as that bounded by the
planes Z'.A.B, ZBO, Z'aD, ZD.A.' ...
2°. Selection by alternate pairs intersecting in the principal
plane of symmetry (ro0 y0). Of the 8 bounding planes, 4 are
obtained from (h, lc, l) by taking the first two letters either (h, lc)
or (- h, -lc) or (lc, - h) or ( -lc, h), and the other 4 from these by
changing the sign of l. The resulting parallel-faced hemihedral
form has three planes of symmetry, viz.: the plane (ro0 , Yo) and
86 STRESB-STRAIN RELATIONS. [36
planes parallel to the sides of the square base in this plane, and
formula (22) holds for these forms (fig. 8).
Z'
Fig. 8.
through 60° about the axis z0 • When the new axes are given by
the equations
~=!, ~=!.V3, n,.=O}
l 2 = - !.Y3, m 2 = t• ~ = 0 ............... (25),
l8 =0, ma= 0, na= 1
the equations of transformation of strain-components, given in (33)
of art. 10, become
e' = }e + !f + !.V3c, a'= !a- !.V3b, }
f = !e + l/ -l.V 3c, b' = !.V3a + !b, ... (26).
g' =g, c' =-!.V3e+!.V3/- ic
The equations of transformation of stress-components given in (14)
of art. 16, give P', Q', ... in terms of P, Q, ... If we write down
the corresponding formulre for P, Q, ... in terms of P', Q', ... we
shall get
P = }P' + !Q'- i.V3 U', 8 = !8' + !.V3T, }
Q =!P' +lQ' +!.V3U', T=-!.V38' + !T', ... (27).
R =R', U = !.V3P'- !.V3Q'- !U'
Now writing equation (21) in the form
2W =(.A, B, C, F, G, H'§.eo,fo, gof+Lao2 +Mbo2 +Nco2,
substituting for 8' and T' in the 8, T equations of (27), and
equating coefficients of a or b, we obtain L = M.
Substituting for R' in the R equation of (~7), and equating
coefficients of e or f, we obtain F = G.
Substituting for P', Q', U' in the P, Q equations of (27), and
equating coefficients of e in the P equation, and coefficients
off in the Q equation, we get fN = H-.A- -ArB- iH, and
!N = H-B- -&A - !H, from which .A = B, and N = ! (.A -H).
Thus the energy-function is given by
2W=(~l>~l> aaa, ~. ~. ~l_eo,Jo.gof+a"(ao2 +ho')+H~~-al2)co2
............ (28).
Note that this formula is unaltered by turning the axes of x0 and
Yo through any angle.
Beryl is an example of a crystal for which this formula holds.
39. Rhombohedral System-(6 Constants).
The most important hemihedral forms of the hexagonal system
are the rhombohedrons obtained from a hexagonal pyramid by
39] SYSTEMS OF CRYSTALS. 89
the selection of alternate planes. In figure 9 ABOA'B'O' is a
regular hexagon, and ZZ' a perpendicular axis, and the faces of
the rhombohedron are ZAB, Z'BO, ZOA', Z'A'B', ZB'O', and Z'O'A.
These forms are unaltered by rotation through 120° about the axis
z0 , and also by rotation through 180° about the line AA', which
we take for the axis rc0 • We have already seen (art. 34) that the
last property produces just the same simplification in the energy-
function as if the plane rc0 = 0 were a plane of symmetry, and we
may therefore set out from the form
2 W =(au, ~~ £Zaa, a" ... 2:eo, fo, 9o 1 ao'f +(a~~~~, atl81 ~9:bo, Co)2•
If we work out the conditions that this may be unaltered by a
Z'
Fig. 9.
90 STRESS-STRAIN RELATIONS. [40
rotation through 120° about the axis z0 , we shall find the energy-
function for this crystal system given by
2W =(au, au, asa, a_, aw, aaleo,J~. go)t +a"' (ao' + bo')+t(au -all) Co2
+ 2~,ao (eo-/o) + 2~,boeo............ (29).
After what we have just done for the holohedrons of the hexagonal
·system, the work may be left to the reader. Formula (29) holds
for Iceland Spar.
Quartz is an example of a tetartohedral form of the hexagonal
system, which is a hemihedral form of the rhombohedral system,
and has the property of being unaltered by the same rotations as
the rhombohedron. Formula (29) therefore holds for this mineral.
(; -
1
~) m...man . .na + (;. . -1) n2nsl2la + (~ -1) lJ,m...m J 8
2 2
+ (m...na + man-...) + (rt/-a of: nals) + (l2ma + lam2'f (35)
L M N ········· '
1 See the • Annotated Clebsch,' pp. 95 sq.
44] MODULUSES. 95
and the result (31) of the article preceding is that the resistance
to compression is the reciprocal of
1112 2 2111
E1 + E3 + Ea + F1 + F2 + Fs- L- M- N ......... (a 1)·
tTw = E2 ( 2~- j)
o-31 = E, (2~- jJ
0"311= Ea (_!_-
2L F
_!_)
1
1 Bee ta.ble, p. 77. The moduluses for different specimens of sieel vary
considerably.
2 These numbers shew that there really exist materials which poBBess what Sir
U=Nc J
where n is a number. These satisfy all the above conditions, and
are put forward by Saint-Venant as likely to prove more useful
in practical applications than the formulre for an isotropic solid.
(See the "Annotated Clebsch," p. 107.) In case the axes of x
and y are interchangeable, we can simplify these formulre by
putting L=M.
It is easy to shew that
E N•
;y1 = 8g (Ea lJ - 1) ;· (E,
~ - 1) .............. ( 47).
1
in 1694, but they were really introduced into practical treatises by Poncelet in 1839.
1 See Prof. Unwin's treatise on the Testing of Materials of Construction.
102 STRENGTH OF MATERIALS. [52
Then comes a stage in which the curve is generally concave down-
wards, so that the strain increases faster than it would do if it were
proportional to the traction; in this stage the strain is largely a
permanent set. As the traction increases there comes a region of
well-marked discontinuity, in which a small increase of traction
produces a large increase of set. The traction at the beginning of
this region is called the Yield-Point 1• If the traction be increased
above the yield-point, local How of the material generally takes
place at some point of the bar, determined apparently by accidental
circumstances, and the bar ultimately breaks.
52. Elaatic Limits.
The strain is in general wholly elastic, provided it be not
greater than a certain limit, i.e. if a certain limiting stress be not
exceeded, the strained body recovers its previous form and dimen-
sions on the removal of the load. The elastic limits are by no
means constant, even for the same kind of material, but appear to
depend very much upon the previous history of the body tested.
Cast iron is an example of a material which, in the natural state,
cannot be subjected to any finite stress without receiving a
permanent set. It is well known that the elastic limits may be
gradually raised by the application of loads that produce a perma-
nent set. Thus a bar of cast iron, after having received a set by
the application of a small finite load, may be strained by the
application of an equal or slightly greater load, without taking any
increase of set. It is possible that in cases of this kind the initial
set consists in the removal of a state of initial strain, or, in other
words, that the application of a load to the body is required to
reduce it to a state of ease. It is important to remember that
there are in general two elastic limits for the same kind of strain ;
thus, for most materials, the limiting elastic extension differs from
the limiting elastic compression. It is also important to remember
that the elastic limits may be considerably widened by the gradual
application and removal of loads that produce small permanent
sets. The experimental determination of elastic limits is beset by
difficulties, among which not the least important is that the elastic
1 The phenomenon appears to have been firs1 noticed by Dufour in Geneva in
or
and in like manner we may write down the values of r~, rv after
strain.
The new value of X (r) is
X (r) + erx' (r) ........................ (10).
The new value of p is
p [1-(e+ f + g)]=p' say .................. (11).
61] POINT-ATOM HYPOTHESIS. 113
Thus P, U, T become
8-2
116
remember that
ou
e= oa;' ... ,
= JJl
0
0~ SudS- JfJ~ ( 0~)auda:dydz ... ............. (29),
0
u= o, " = o, w= o, ou ~
az = o, a:z = o.
~
ax= o.
the expression for the displacements with given strains will be
unique. These particular equations indicate that a point of the
solid (the origin) retains its primitive position, that one line-
element of the solid (along the axis of z) retains its direction, and
that one plane-element of the solid (the zx plane) retains its
direction unaltered by the strain. It is manifestly possible, having
strained the solid in an arbitrary manner, to bring it back, by
translation and rotation, so that this point, this line, and this plane
shall recover their primitive positions.
(o) The strains cannot be given arbitrarily'.
1 This theorem was given by Saint-Venant in his edition of Navier's Ltt;OTI8,
Appendice III.
122 GENERAL THEOREM& [66
Suppose we have the equations
~ ~ ~ ~ ~ ~ ~ ~ ~
e=OX' I= oy' g = oz ' a= oy + oz' b= oz + OX ' c=ax+ oy <37),
and seek the conditions of compatibility.
The displacements u, v, w must be one-valued functions of
x, y, z, so that, if we integrate along any line drawn through the
origin, we must have
OU ou ou
I
u = Uo + OX dx + oy dy + oz dz ..................(38),
~= f'Ou) o +J~
oy \}Jy oy dx+ (oc ox dy+t(~
oy _of) oz + ob
oy _ oa)
ox dz} ...(41 ).
~=(ou)
oz oz o +Joedx+
oz "2" oz
1
(oc+~-~)d
oy oa; Y+(~-?f!)az
oz ox
The line integrals will be independent of the path of integra-
tion, if
66] UNIQUENESS OF SOLUTIONS. 123
to which may be added
a~f a2g _ o'a
oz' + olf- oyoz ........................ (43),
one of the equations obtained by operating in like manner upon v.
These six equations are the necessary and sufficient conditions
of compatibility of equations (37).
(e) If either the surface-displacements or the surface-tractions
be given, the solution of the general equations of equilibrium is
umque.
1°, Supposing the bodily forces and surface-tractions given,
then, taking W a quadratic function of the six strains, we have
ow
ae=P, ......
also the general equations of equilibrium are three such as
aP au aT
ax + ay + a-z + pX = 0 .................. <44),
and the boundary-conditions are three such as
lP+mU +nT=F ..................... (45).
If possible suppose there are two different solutions of these
sets of equations, and that the corresponding displacements are
'U:!, v1 , w1 and u2 , v,, w~ in the two solutions. Then, writing
u' =U:i -Us, v' = V1 w' = w1 -w1 ,
-v2>
we see that u', v', w' are a set of displacements which satisfy three
such differential equations as
aP' aU' aT'
OX + oy + (}z = 0 ..................... (46 ),
and three such boundary-conditions as
lP' +mU' +nT' = 0 ..................... (47),
where P', ... are the stresses corresponding to these displacements.
Now, by Green's transformation,
fff(e oW'
oe + ... + c'
oW')
oc' dxdydz = 0 .......... (48).
But, from the form of W', as a positive quadratic function, we
know that the expression under the integral sign is 2 W', so that
the integral is a sum of positive terms, which can vanish only
when e' = j' = ... = c' = 0. Thus the displacements (u', v', w') are
such as are possible for a rigid body, and the solution is only inde-
terminate to the extent of such displacements.
2°. Supposing the bodily forces and surface-displacements
given, we take as before two solutions u1 , v1 , w~> u 2 , V2 , W 2 , and
form their differences u', v', w', then u', v', w' satisfy stress-
equations like
oP' oU' oT'
+ oy + oz = 0•
ox
and boundary-conditions u' = 0, v' = 0, w' = 0, at the surface.
Thus we find that
- ffr , )
. (u- u
(oP1
oa: + oPe oP6) dxd d
oy + az y z
-(similar terms in v'-v and w' -w)+ ffff(e' -e)dxdydz ... (54),
and therefore, since the function f(e' -e) is necessarily positive, it
follows from (53) that V' is always greater than V, if the equation
JJJ [P (X1- a;) Us+ p ( Y1- ~;~) vs + p ( zl- a;;~) w2J rUJ;dydz
+ JJ(F1Us + G1vs + Hlws) dB
= JJJ [P (X~-~) Ut +p ( Y2- ~~~) vl + p (z~- a;;~) w~J rUJ;dydz
+ JJ< F2Ut + Gsvl + Hsw1) dB ......( 56),
the volume integrations being taken over all points within the
surface B of the solid.
To prove the theorem, let P1o Q11 ~. Bt, T1o U1 be the six
stresses, and e1 ,ft, g1 , ~. ht. Ct the six strains, corresponding to the
first set of displacements, and like expressions with suffix 2 those
corresponding to the second set, and let W1 be the energy-function
for the first set, and W2 that for the second set. Then, by the
equations of equilibrium,·the left hand side of (56) is
-JJJ[(~P~ + ~u~ + oT1) Us+ (aul + aQl + ?§.) v2
oa: ay az oa: ay az
+ e;: 1
+a:+ a:;) w 2] da:dydz + Jf(F1Us+ G1v1+ Hlws)dS
= fff<Ples+Q~~+~gs+B~as+Tlbs+ U~~)dP;dydz
- JJ[Us(lPl + mU1 + nT1 -F1)+ Vs(lU1+mQ1+nB1- G1)
+ W 2 (lT1 + mB1 + n~- H1)] dB.
The surface integral vanishes identically in virtue of the
boundary-conditions, and the volume integral is
ff.f [61 oW
oes + '"oWs
2
1 Cifs + g1oWs
J og, + ~ oW.
oas + boW,
1 ob, + c1oWs]
oc. .3-dydz, ua;
In order that P, may be unity and the rest zero, we must have
"A.+p. 1 }\. -tT
6:. = p. (3"A. + 2p.) =E. f, = g,=- 2p. (3"A. + 2p.) =E. as= bs= Cs =Oo
In order that S, may be unity and the rest zero, we must have
1
~ = -, e,=f, =g~=O, b,=c2 = Oo
p.
Hence, suppressing the suffix 1, we find the mean values of
the six strains in the following forms :
mean value of a=
2~V Jj(Hy + Gz) dS ..................... (63),
where Vis the volume of the body, and E is the Young's Modulus,
p. the Rigidity, and u the Poisson's ratio of the material, and these
three constants are connected by the relation
E = 2p. (1 + u ).
We notice that the mean value of the cubical compression is
L. 9
130
("A.+"')
oA
az + p.Vtw = P o-otw
1
1
See Poisson's memoir of April 14, 1828, Mem. Pari• .A.cad. 1829, and Lord
Rayleigh's Theory of Sound, vol. n. ch. uv.
70] PROPAGATION OF DISTURBANCE. 131
h~~li=~A
ot1 }
a~.ar ............ (71).
/CJV2fiT = ()tt , fiT= (fiT~> fiT2> fiT a)
a!=~~}
at-~
.............. ('74),
so that, if the initial values of~ and 00~ be denoted by ~0 and c/J 0,
sinh (atV) ·
cfJ = cosh (atV) c/Jo + t atV ~o ............ ('7 5),
where ~0 , ~ are functions of :c, y, z.
Observing that cosh (atV), and sin~~tV) are even functions
of V, we see that these are real operators, and the operations
indicated can be performed.
But there is another form into which the solution can be
thrown, in virtue of the theorem that the mean value of a function
+ over the surface of a sphere of radius R, whose centre is the
.
ongm, . t h e va1ue ofsinh
. IS -RV-(RV) ~r..
,.. at t he ongm.
..
9-2
132 GENERAL THEOREMS. (70
For, consider the function ff fJ"+IIIJ+ezdS, the integration ex-
tending over the surface of a sphere, whose centre is the origin,
and whose radius is R. Changing the axes, so that the new axis
of Z may be the normal to the plane aa; + by + cz = 0, we see that
z..; a1 + bS + as = aa; + by+ cz,
also dS = 2w-RdZ, and the integration for Z is taken between the
limits R and - R.
47TR11 •
=-RV smhRV,
if V 11 =a1 +h'+&.
Now suppose ,Y is any function of x, y, z, uniform within a
sphere whose centre is the origin and radius R,
where
is the mean value of 4>o over the surface of this sphere, and
sinh (atV) _,_
0
atV "'
is the mean value of rp0 over the surface of this sphere, thus
0 - _,_
c/> =at (tcf>o) + tl/lo •. •• .. ·••.• •·.•.....• ('77),
where i/)0 and ~ are the mean values of the initial cf> and 4> at all
points of the surface of a sphere whose radius is at and centre the
point at which the disturbance is to be estimated.
and therefore
wiT = l (~~ + 'Xsf'J + Ast'f
= H~. ~. X., As., ~. ~lE, ,, r1 ...............(100),
138 GENERAL THEOREMS. [76
where
~11 =Xl = (c,.l + c8m + c0n)1 = (Cn, c., C00 , 0118 , 0,.01 c111 Il, m, n)1 ••• (101),
~ ~=~=W+~+~~+~+~
=(C,.e.c.,cu, Heso+ce]. ![c,.,+cll8], Hc,.,+c.]Il, m, nf
......... (102),
and the other coefficients can be written down in like manner.
The function II is thus a complete quadratic function of ~. "''
~. a.nd, since W' is always positive, ID1II also is always positive.
Now from the equations (97) we find
piD~~ = IDs 00~ = ~n~ + "'A,ufJ + ~~.
and two similar equations.
Hence ID1 must satisfy the determina.ntal equation
2
~- ID p, ~~~ Xu
Xu, X.- ID1p, 'A.a = 0 ...... (103).
Xu, ~~ Aaa- ID2 p
Since the function ID1Il is always positive, the roots of this equation
are all real a.nd positive.
Thus there are in general three real values of ID, the velocity
with which the tangent pla.ne to the surface of discontinuity
advances, and these are functions of (l, m, n) the direction of the
tangent plane.
76. Equation• of a ray.
Let ao, /30 , 'Yo be any point on the surface 80 when t = 0, then
the parallel tangent plane at time t is
l (a: -ao) + m (y- /30) + n (z- ry0) = IDt,
and this contains the point
OID OID OID)
( IXo + t ol ' /3o + t om ' 'Yo + t On '
"f
l l OID OID +nOn
Ol +mOm OID =ID •••••••••••••••
(104)•
78. Wave-Motion,
We can now give a sketch of the propagation of the disturbance
through the medium. For this purpose we shall suppose that
initially the part of the medium outside a certain surface S 0 is
unstrained, and the medium within the surface is strained in a
given manner. If then we draw the central-surface corresponding
to time t for every point within S 0 , these surfaces will have an
envelope S, which will consist in general of six sheets, two for
each value of w. Fixing our attention on one value of w and the
corresponding sheets of S the motion of this type will be called a
wave. Three such waves are propagated. The parts of the
medium, not included between the two sheets of S, corresponding
to a wave are at rest and unstrained. Every element of the
medium when the wave reaches it takes suddenly the small dis-
placement-velocity propagated with the wave. After a time
depending on its position with respect to the original region of
disturbance (the space within S 0 ), the wave will have passed over
this element, and as the inner sheet of S passes over it the
element suddenly loses the small velocity that it had, and returns
to a position of rest and a configuration of no strain. The same
thing happens for each of the three waves. 'fhe element, if it be
far enough from S0 , is jerked into motion from rest, and returns
impulsively to rest from motion by the action of three separate
impulses, and its motion in each case lasts for a finite time
depending on the size of 80 • In every case the whole motion
depends simply on the form of the wave-surface and on the initial
state.
The particular case of an isotropic solid is an example of a
case in which the determinantal equation for w has two equal roots,
the lc of our previous work. The reader will find it an instructive
exercise to work out this case, and also the case of a medium
whose energy-function is of the form
l {A (e+ f +g)~+L (a9 -4jg) + M(b2 -4ge) +N(lf- 4ej)},
which leads to Fresnel's wave-surface 1 and a sphere as the general
wave-surface.
1 See Jlath. Paper• of the late George Ckeen, pp. 308--S05.
141
Jjj8Wda:dydz
+ JJJ{p ~; Su + p~ Sv + p: &w}da:dydz = 0 ....... (10'1).
If after performing the variations we put u = u,.tf>r, 'IJ = Vrtf>r,
w = Wrtf>r, and observe that
Q2u
ott= 'Urt/>rPr',
and so on, we see that tf>r will be a factor which can be removed
from the resulting equations, and the part that arises from 8W
will be the same as if we substitute (u,., Vr, wr) for (u, 'IJ, w) in the
expressions
aP au aT au aQ as aT as aR
~+~+~, ~+~+~, ~+~+~·
If Pr, Qr... Ur denote the values of P, Q, ... U when u,., 'IJr, Wr are
substituted for u, v, w the equations of vibration are
L. 10
CHAPTER VI.
1 Saint-Veuant began with a solid which has one plane of symmetry only,
perpendicular to the axis of the beam, but introduced the other two planes after-
wards to simplify the work. The student reading the subject for the first time is
advised to work over all the general theory for the case of an isotropic beam.
As a further example of the analysis in the next article it may be shewn that, if
Saint-Venant's stress conditions be imposed and the beam be supposed vibrating,
equation (12) will be satisfied, and equations (11) and (18) become
Ollw p (1 + 2cr) 0'10
o.zJ = 2~-' 11 + lr} oz at• .................................... (H),
10-2
148 SAINT-VENANT'S PROBLEJrL [83
two of which intersect in the axis of the beam, but we shall not at
:first take it to be isotropic.
83. Equatiou of the Problem.
Take then the axis of z parallel to the length of the beam, and
suppose that it is the line of centres of inertia of the normal
sections in the unstra.ined state, and suppose the energy-function
of the material when strained to be W, where
2W = (A,B, C, F, G,HJ...e,f,g)2 +La'+MI:r+Nc' ... (1),
so that the stresses are given by the equations
P=Ae+Hf+ Gg, S=La}
Q=He+ Bf+Fg, T=Mb ··············· (2).
R=Ge+Ff+Cg, U=Nc
The stress-conditions imposed by the semi-inverse method are
P=Q= U=0 ........................ (3).
01w trp OStD aaw
and ~=- 2#£(l+tr) azae•= oza,• ........................ (18),
whm·e tr is the Poisson's ratio of the material, supposed isotropic, and #' is the
rigidity, Equations (11), (12), and (18) cannot be satisfied UDless ~: is independent
of :e and y. The equation corresponding to (10) is
a~w aawa~to
&2 + ay• + 2 1 - ~ oz _ ae• ...........................(lO),
p OStD
TheBB are similar to the vibrations of water in a cylindrical tank whose curved
surface coincides with that of the beam.
84] GENERAL EQUATIONS. 149
oT =O oS=O}
;; +0~ +~~ =0 ..................... (4),
ox oy oz
and the only condition at the cylindrical boundary, which is not
satisfied identically, is
lT+mS=0 ........................... (5),
where l, m are the cosines of the angles which the normal to the
boundary drawn outwards makes with the axes of x andy.
We may also suppose the geometrical conditions satisfied at the
origin to be
u=O, v=O, w=O, au en, Ou
oz=O, oz=O, oy=()1 ... (6).
Then the problem consists in the discovery of the most general
solution of the equations (3), (4), (5), (6), and the determination of
the consequent amount and distribution of force over the plane
ends of the cylinder 1•
where 1T1 and 1T1 are the ratios of lateral contraction, parallel to x
andy respectively, to longitudinal extension parallel to z.
Also, since U = 0, we must have
Ou en,
oy + ox= 0 · .............. · · .. · .. · ....(8>·
1 These equations denote that the origin is supposed held fixed, that the
element of the axis of the beam at the origin retains iw primitive direction, and
the element of the plane through it and the am 11 retains its primitive direction.
H any other conditions be imposed at the origin the displacements consist of thoae
that we shall obtain combined with a suitable rigid body displaoement.
• The problem in this form was first considered by Clebsch. (See lnVoduoUon.)
150 SAINT-VENANT'~:; PROBLEM. [85
aT as
And, since oz = 0 and oz = 0, we must have
~~ =0}
+ ()Jw
::: : ........................ (9).
-+-=0
oz2 oyoz
The third equation of (4) is
O'u O'w) (
asv Qlw) O'u ()tv O'w
M ( ozo.x + ox' + L 'Oyoz + ays + G ozo.x + F o60y + O 2 = O, az
or by ('7),
asw asw
Jf ox' +Lays+ [0- (M +G) u1- (L + F)u,] az•
'O'w
=0. . (10).
+ !fJtY1(cr2- f)]
R= (0- Gcr1 - Fcr1) [tx + «1X + a,y + z ({3 + fJ.x+ {3.y)]
The condition (5) at the cylindrical boundary becomes
+ fJ, E-Lu
M (Ml:cy+tLma:')- Lmu1u:'] ...(24)
[ tLmu~~Y' +
1
'!/:
:c--tazt}
0 l ••.••••••••••••••••••••• (32),
0 y
Fig. 10.
If this figure be supposed raised through the distance necessary
to bring its centre on to the circle formed by the central-line, and
turned round the axis y until the axis a: produced passes through
the centre of this circle, we shall have a complete representation of
the changes of form and position involved, and the representation
applies equally well to the distortion of a beam of rectangular
section and to the distortion of a. rectangular portion of a. beam of
any section.
1 The curvature in the figure is much exaggerated.
92] TORSION-COUPLE. 15'7
u=-T'!JZ} ...........................(35),
::;a:
where cfJ satisfies the differential equation
!![M {a:L(~!-
0: TY )} +L:y {a:(~+
Tfl: )}] rkdy,
since (36) holds at all points to which the integration extends ;
and this can be transformed to
P jj•·(a: a<I>
a:iJ -y a<I>)
d.x d.xdy .................. (43),
where pis the density, and this is -Cclp(q-1)1.
If we suppose the vessel constructed of such material that its
moment of inertia about the axis of the cylinder is -pi, the whole
impulse required to start the motion will be - pCc1ql, so that,
identifying p and p., the impulse in the hydrodynamical problem
will be identical with the couple in the elastic problem. The
hydrodynamical problem is however no longer a real physical
problem as it involves a negative distribution of matter on the
surface of the cylinder.
The hydrodynamical analogy suggests the method to be
followed in the solution of the torsion-problem. We know that in
irrotational motion of a liquid in two dimensions there exists a
stream-function 'I', which is the conjugate-function of <I> with
respect to a; and y, and that the value of 'I' is given at the boundary,
and It. 18
. 1D
. genera1 Simp. Ier to SO1Ve t h e equatiOn· -QIV_$1 + o'V
oyi = 0
0
when the value of V is given at any boundary than when the value
of~~ (rate of variation in the direction of the normal) is given at
the boundary.
We shall accordingly suppose that cfJ and '+' are conjugate
functions of a; and y, so that cfJ + t'+' is a function of the complex
variable a;+ ''!I· then we know that cfJ and '+' satisfy the same
partial differential equation
o' Ql \
(ofli'l cfJ
+ oy2 ),. = 0•
and we have
ocfJ _ a,.
ocfJ _ a,.
-o:C- oy' oy-- ()a; ..................... (44).
We have to obtain the boundary-condition for '+'·
160 THE TORSION-PROBLEM. [94
Let ds be an element of arc of the bounding curve of a normal
section of the cylinder, measured in that direction in which the
Pig. 11.
dy ~t + dUJ oy =
ds oy ds oa:
'T (a: dtx + Y ~?!)
ds ds
.
We may obtain the integral of this equation in the form
y= tT(afJ+~)+ a..................... (46),
where a is an arbitrary constant.
Thus the problem is reduced to finding a solution of the
equation
tAfr 0~
oai' + oy' = o........................... (47)
at a.ll points of the section, subject to the condition
y. =iT (a:t + y2) +a..................... (48)
at a.ll points of the boundary.
96] STRENGTH UNDER TORSION. 161
95. Hydrodynamical Analogy, continued.
There is another form of the hydrodynamical analogy in which
the problem of torsion is compared with that of a liquid circulating
with uniform spin in a fixed cylindrical vessel coinciding with the
surface of the twisted prism. In this form the moment of
momentum of the liquid coincides with the torsional rigidity of
the prism, and the velocity of the liquid at any point represents
the shear at that point of the prism. In the first form of the
analogy the shear is represented by the velocity of the liquid
relative to the vessel.
A result of considerable importance 1 can be at once deduced
from the bydrodynamical analogy. Suppose a shaft transmitting
a couple to contain a cylindrical flaw whose axis is parallel to that
of the shaft. If the diameter of the cavity be small compared with
that of the shaft, and the cavity be at a distance from the surface
great compared with its diameter, the problem is very nearly the
same as that of liquid streaming past a cylinder. Now we know
that the velocity of liquid streaming past a cylinder bas a maximum
value equal to twice the velocity of the stream, and we may
therefore conclude that in the case of the shaft the shear near the
cavity is twice as great as that at a distance. The importance of
this result in connexion with the strength of materials bas been
previously pointed out (art. 58).
96. Strength or a beam under torsion.
With our notation the six strains of the material are
olf>
e=f=g=c=O, a= oy +TX, b= oa; -ry,
o4>
Fig. 12.
or
the stream-lines are therefore similar ellipses, and the distance
between consecutives is proportional to the perpendicular from the
centre on the tangent. The maximum velocity along a given
stream-line has place therefore on the minor axis. The velocity
at any point on the minor axis is proportional to the distance from
the centre, and therefore, in the torsion-problem, the shear is
99] TRIANGULAR PRISM, 165
greatest at the extremities of the minor axis of the ellipse, and is
there equal to 2Tat.bj(a2 + b2).
99. The equilateral triangular prism.
!I
Fig. 13.
The curves of equal distortion are shown in fig. 13. The triangle
is divided into 6 exactly similar and equal parts by the lines
joining the corners to the middle points of the sides, and, in these
parts, the section is alternately depressed and elevated.
In this, as in any case where there is a sharp corner which is
not re-entrant, the velocity of the fluid circulating in the prism
is zero at the corner, and it is obvious from .symmetry that its
166 THE TORSION-PROBLEM. [100
maximum in the case of the triangle has place at the middle
points of the sides. It may be verified by direct calculation that
the maximum shear in the twisted triangular prism has place at
the middle points of the sides.
provided
"" (2n+ 1)wy
T(~-b2 )=~An COS --~ ............ (63),
when b>y>-b.
Now between these limits we find, by expanding y 2 - b2 in
Fourier's series,
(-1)"
(2)8
An=- 4Th~ ;r (2n + l)3 . . . . . . . . . . . . . . . (64).
Thus
so that
q,=-TX'!f
sinh (~~±_l) wx ]
2( 2) 8 "" 2b (- . ~ 2n + 1) ?rY
1)"
+4Tb ;. ~ [ (2n+1)3 h(2n+1)7rasm- 2b ... ( 66 )·
cos 2b
The twisting couple is, by (39) of art. 93,
4p.Tab a
2
; b' -JJ-T JJ (:#- y') dxdy
+ •"'Tb' (~Y ff (ia!- y ~!) clxdy,
which is equal to
fp.TabB + 4~TIJ2 f
(~r J(X~~- Y ~:) dxdy .. ,. .. (67),
where
(2n + 1) 7rX • (2n + 1) 7rY
1)" inh 2b 2b
<I>-; (2n + f)a
_ 00 (- S SID
h (2n + 1) wa ... (GS).
cos 2b
A term of the double integra] is
(2n+l)7ra (-1)" ?r
sec h ---2b-- (2n + lf 2b
Jj [xsm. h (2n+2bl)7rx cos (2n+2bl)7ry
- y COS h
(2n + 1) wx . (2n +I) 7r!f] d d
2b SID
2b X '!J•
168 THE TORSION-PROBLEM. [100
Now
. h (2n +b1)
J
~ 7rX d
-a XSID 2 x
2b [ h(2n+1)7ra 2b . h(2n+1)7ra]
=--,(2,..-n-+-1) 7r 2a cos 2b - (2n + 1) 1r 2 SID · - 2b- '
and
so
J cos (2n+1)7ryd
fJ
-b 2b
=. 2b
y (2n + 1) 7r
2 (- 1 )11
'
J
~ h (2n+ 1)7rx ..]_ _ 2b . h (2n+ 1) 'Ira
-a cos 2b r.w;- (2n + 1) 7r
2 SID 2b
Hence the twisting couple is
i,u.Tabs +
2)s "" (2b)' 8tanh (2n ~!)'Ira
- 4,u.Tb' (;. ; 7r (2n + 1)6 '
which is equal to
siDce
and it is to be noticed that a tenn i~J-Tab3 is contributed by the
transcendental part of ¢.
The expression for ¢ must really be unaltered when x is
changed into y, and a into b, and the expression for the twisting
couple must be symmetrical in a and b. For an account of the
identities thus obtained the reader is referred to a paper by
Mr F. Purser in the Messenger of Mathematics, XI. 1882.
Saint-Venant has investigated 1 the fonns of the curves of equal
distortion given by ¢ = const. If we begin with the case of a
square prism cf> vanishes along the diagonals, and along the middle
lines of the square parallel to its sides. If we take one side great
compared with the adjacent side, then cf> vanishes only on the
middle lines and not on the diagonals. The limiting form between
1 See the great memoir on ' Torsion ' of 1855.
101] SECTORS. 169
rectangles which divide into 8 parts, in which cf> is alternately
positive and negative, and those which divide into 4 such parts, is
given by making the ratio of the sides equal to 1·4513.
101. Other sections.
When the section is not of one of the forms just considered
the problem can frequently be solved by means of conjugate
functions. Whenever we know a transformation by means of
conjugate functions, say
a+ t/3 = f(x + t.y),
such that the boundary consists partly of lines a= const. and
partly of ·lines /3 = const., the differential equation for 'ir can be
expressed in the form
(}2+ (}2+
oa2 + '0{32 = 0,
and Y. will be a given function of a along some part of the
boundary, and of /3 along the remainder. The simplest case is
that of a curvilinear rectangle bounded by two arcs of concentric
circles and two radii, including the case of a sector 1 of a circle of
any angle. The work in these cases may be left to the reader, we
give the results.
1°. For a sector of angle 2/3 in terms of polar coordinates
r, 8, so that the boundaries are
r=O, r=a, 8=±/3,
(r)(2n+1)2~ 1r8] ... (70),
'\fr = 2 cos
tT'r COS
28
213
+ Ta 2 «> [
; A.'l7Hl a COS (2n + 1)
213
where
A
mH- -
_ ( )n+1 [ 1 2 1
(2n + 1) 1r- 4/3- (2n + 1) 1r + (2n + 1) 1r + 4/3
J
......... (71).
If we write re' 9 =ax, then
af
'1/r- t.cf> = i-ra2 cos 2/3
-ra~ {af Ja;o ~dx---;:
- 2f)
a;~-3 4/3 ( ~) 1 Ja; a;l!~+l }
tan-1 a;2fJ + a;2 o ------;__& ... (72),
1+~ 1+#
where the modulus of a:: is ~ 1, and ta.n-1 (afo) is that branch
of the function which vanishes with a;,
1 Greenhill, Messenger of Mathematics, VIII. 1879, and x. 1881.
170 THE TORSION-PROBLEM. flO I
where w is the area and I the moment of inertia about the axis.
Now Saint-Venant 1 found that, in the case of all sections such
as those just described that had been worked out,
w•
G ="I p.T,
where the value of " varies only from ·0228 to ·026 while its mean
value is about ·025, or -iJ, and its value for an ellipse is ·02533 ... ,
and thus we may take with remarkable accuracy for most forms
of section likely to occur in practice
(1)4
G = irr y P.T •••••••••••••••••••••••• (77),
and the twisting couple for different prisms of the same material
is directly as the fourth power of the area of the cross-section, and
inversely as the moment of inertia of this section about the axis.
The theory of Coulomb and Navier made the couple directly
proportional to the moment of inertia I.
103. Tonlon of JEolotroptc Rectangular Beam.
On account of its importance we shall give the solution of the
problem of the torsion of the rectangular prism in the case where
the two principal rigidities of the material for shear of planes
through the axis are not the same.
We have to find a function 4> to satisfy the differential
equation
CJ24> o~q,
M aa,a +L oy2 = o .....................(78),
and the conditions
when y' = ± J( ~ L M) b
...... (83),
J(L~M)b=b'
y has to satisfy the conditions
~t = T Y 1 1
•• , , , , , . . . . . . . ••••••••••••••(84),
and ~~ = 1
7' /1J
1
............................. (85),
Then the angle turned through by the end at which the couple
G is applied is
3Gl
Moot' ( 1 - 3 flr) ......................(90).
16
Now when r is not < 3, I is remarkably nearly constant and equal
174 THE PROBLEM OF FLEXURE. [105
z{31E JJ xydxdy,
which vanishes identically when the stress parallel toy vanishes.
The couple about an axis in the section parallel to the axis y is
ff( Lx ocp
oy
-My ocb) dxdy
ax
-/31M ff[(E- Mul)(x2y/M- }y8/L)- !u1:x2y + uq] dxdy . .. (97),
and, as in the case of torsion, the first term of this can only be
calculated when cp is known.
We shall suppose that the principal axes of inertia of the
section of the beam coincide with its axes of symmetry of elastic
structure, and then we see that the resultant stress at any section
reduces to a transverse force parallel to· the axis x, a bending
couple about the axis y, and a torsional couple about the axis z.
By a combination of this solution with the previous solutions,
in which there was either simply a bending couple about the
axis of y, or simply a torsional couple about the axis of z, we
may obtain the solution for any resultant stress-system.
106. Practical Utility or the Solution.
Now although the resultant stress-system may consist of any
extending force, any transverse force, any twisting couple, and any
bending couple whatever, yet the mode found for the application
of these forces is unique and determinate. It is very unlikely
107] EQUIPOLLENT LOADS, 177
that in the pieces used in architecture or engineering the stresses
applied at an end of a bar are distributed over the end in the
manner here considered, and the application of the theory to
structures depends on the validity of a principle introduced by
Saint-Venant 1 which may be called the principle of the " elastic
equivalence of statically equipollent loads". He stated that any
system of surface-tractions applied to a small part of a surface of
a solid will produce at any point not very near to that surface the
same strain as any other system of forces having the same resultant
would produce. When the part of the surface to which the forces
are applied is vanishingly small the principle is obviously true,
but it becomes less and less exact as the surface subjected to
traction becomes larger. The principle is in fact equivalent to
this-the application of an equilibrating system of forces to a
part of the surface of a solid produces no strain in the interior-
a theorem which is in general obviously untrue, but which becomes
more and more exact as the surface subjected to traction is
diminished. It is on account of the assumed smallness of the
linear dimensions of the cross-section of the beam compared with
its length that the principle may be applied to obtain a theory
of the equilibrium of beams.
107. Bending by tranavene force.
We wish to obtain in the general case the solution for applied
transverse force X parallel to the axis a:.
Consider first the terms in #1 , these give rise to (1) a stress
E#11 1 parallel to the axis a:, (2) a couple - E#111z about the axis
y, and (3) a couple which we may call T about the axis z. (See
art. 97.) Tis constant along the beam.
The terms in a1 give rise to a couple - Ea1/ 1 about the axis y.
(See art. 89.)
The terms in #o or T give rise to a couple about the axis z (see
art. 92), and since T is at our disposal we may make this couple
equal to-T. This couple is constant all along the beam.
If then we take a1 + #1l = 0, the stress-system at z = l will
reduce to a force X where
X= E/3111 . . . . . . . . . . . . . . . . . . . . . . . . . . . (98).
1 See the memoir on ' Torsion ' of 1855.
L. 12
178 THE PROBLEM OF FLEXURE. [107
A z
Fig. 14.
the part BP on the part .AP, and the action of .AP on BP is equal
and opposite to this.
Now it may be observed that the equation
d2rc
EI1 dzt =X (l- z) ..................... (102),
which is derived at once from the figure by taking moments for
107] TRANSVERSE FORCE. 179
the equilibrium of BP about an axis through P parallel to the
axis y, is sufficient, with the conditions
d:r;
a:=O, a-z=O
at z = 0, to determine the form of the central-line.
We shall devote a subsequent chapter to the development of
this remark.
At present we notice that the defl.exion of the axis of the
beam is
t ~~ (l- }z) Z 2
........................ (103),
and the maximum defl.exion is at the end of the beam, and is
equal to
Xl3
! EI~····························(104).
The displacements are
X
u = t ]?/~ [(l- z) (ula:2- IT·J?l) + (l-}z) z 2] - T'!JZ
X
v = EI~ a-2 (l- z) rey + TaJZ •••• (105),
_ ~ X
w- z Ell
[<z - 27 _) a: -
2
~
E- Ma-1 -La-1
L rey
~]
+ .,..
,..
p
Fig. 15.
Psine
+ -Els (l- z) u 1tey- -ryz,
Paine
v= i Ela [(l- z)(u~2 - u1xt) + (l- !z) z2]
...... (107),
Pcose
+ EI (l-z) uatey + -r:cz,
1
_ 1 Pcose[<..!l l) E-Mu1-Lu1 ~J
W-1"-EJl ~-- 2 z :c L ·tey
+ 1l" Psine[(
--- Z 2
- 2l)
Z y-
E-Mu1-Lu2 yw-
_.]
Els M
+ c/>1 + cf>s
where Tis the twist corresponding to the twisting couple-T.
The equations of the central-line are
and the deflexion in the plane through the axis and the direction
of the force P is
l Ep (COS
2
1-;- + sine
2
e)
I;" (l-iz) z2 ............ (109).
The maximum defiexion is at the end of the beam and is
equal to
! Pl
3
E (cos
1
I;- +I;sin2 e e)
.................. (110).
Unless 11 = 12 the axis will not be bent into a plane curve in the
plane through P and the old position of the axis, but the plane in
which the curve lies is given by the equation
:csine_y.cose=O·
Is 11 '
the sine of the angle between this plane and the plane through P
and the axis is
. ( 1 I\ (sin1
smecose 11 - tJ I,' + 112
cos' e)-l e
•
182 THE PROBLEM OF FLEXURE. [109
The plane through the axis perpendicular to the plane in which
bending takes place cuts the cross-section in a line which is conju~
gate, with respect to the ellipse llf/11 + y'fl~ = 1, to the intersection
of the cross-section and the plane through the axis and the direc-
tion of the force P.
This is the theorem of Saint-Venant and Bresse that the
neutral line is conjugate with respect to the ellipse of inertia to
the trace of the load-plane on the cross-section.
109. Strength of a beam under flexure.
To simplify 1 the consideration of this question we shall suppose
that u 1 = u~ = u, and L = M = p.; also we shall suppose the section
such that the twisting couple T = 0, and thus reduce the displace-
ments to
u = !fJ1 [u (l-z) (x'- '!f)+ (l- !z)z'] }
v = {J1u (l- z) 1C!J
...... (111),
w=ifJ~ [<~- 2~)x- (~- 2u) ~ + cf> J
where 4> is a certain function of (x, y) to be determined by solving
the equation
!Jlcf> ()lcf>
aa;s + a'!f = 0•
subject to a. condition given at the boundary.
The six strains are
e = {J1u (l- z) x
f= fJ1u(l- z)x
g=- {11 (l-z)x
c=O
Now in the particular cases that we shall investigate we
shall find that ~ is at least of the order 3 in the linear
dimensions of the cross-section, and thus, if x and y be small
1 The suppositions of this article are not equiftlent to that of transverse
isotropy. Bee Note B at ihe end of ihe volume.
110] STRENGTH OF BEAMS. 183
compared with (l-z), we may for a first approximation neglect
a and b in comparison with e, f, g. It will then be seen that
the greatest principal extension has place at the highest point
on the contour and at the fixed end, a; being measured positive
downwards. According to the theory of Poncelet and Saint-
Venant (art. 57) this extension must not exceed a certain limit.
If T0 be the breaking stress of the material for pull in the
direction of the axis z, and t the thickness of the beam the limit
of safe loading X is given by the equation
ci>X = 2TJ1 j(lt) ..................... (113),
where ci> is a factor of safety.
In case the beam is not very long 1 in comparison with its
breadth and thickness, we cannot neglect the shears a, b, due to
flexure, and given in equations (112). The problem then becomes
very complicated. We have really to transform the elongation
quadric
(e,f, g, la, !b, !o1a;, y, z'f=k
.to its principal axes, so as to obtain the form
fhar +1I'!J1+ g1zl = k,
then the greatest positive value of fh, 11 , or g1 is the greatest
principal extension. The quadric can be found in particular cases
only after the function cf> has been determined.
110. Strength under combined ltrain.
Suppose the strains due to transverse force X (= fJ1El1) to be
as in the last article
e = fJ1rr (l-z)a;, I= fJ1u (l- z) a;, g =- fJ1 (l- z) a;,
and the strains due to torsional couple G to be as in art. 96
a -_ocf>+
oy ...,..,,
,. b-olf>_.,.,J
- oa; ..,,
and suppose o = 0.
The beam will be at once bent to curvature fJ1(l- z) at any
point, and twisted to torsion 'T.
The principal axes of the strain at any point are a line parallel
to the axis z and two rectangular lines in the plane of the section.
It is easy to see that the latter make half-right angles with the
tangent and normal to the stream-line of the circulating fluid
1 If this be not the case the forces whose reeultant is X must be regarded as
distributed over the end in a particular manner. See art. 106.
184 THE PROBLEM OF FLEXURE. [110
fa! ds=O,
the integral being taken right round the boundary.
Now
f+(a:', y')
(of) + (of)'
2
~~ =rll fJl [ial cos• fJ+ (~~-tal- a,) cos() sin51 eJwhen r =a
= a1fJ1 [ (iff -las) cos()+ ( !a1 +}a,- i!) cos J...(117). 3()
Now any function cf> which satisfies (115) within a circle whose
centre is the origin can be expressed in the form
4> = fJ1 [ (1!- !a.) a'rcos fJ+! (la. +}a,-~ ID ,..COB aeJ
= fJ1 [ (i ~ - fa,) a'a: + l ( !a1 + fa,- i ~) (a:'- 3a:y•) J. . (l18).
Fig. 16.
.Af5 we travel along the beam towards the end z = l the outer
circle expands. The double points move off farther from the
inner circle and the parts of the ovals within the inner circle
flatten out. If l > l-v'7 a 1, the outer circle becomes infinite for
z = l - -v'(l' -!a'), and, for greater values of z, a and /3 are
negative, and the outer circle disappears. The curves of equal
distortion may now be written
a:(.#+ys+p)=q,
1 See footnoie p. 183.
113] ISOTROPIC CIRCULAR BEAM. 191
where pis a constant depending on z, and q has the opposite sign
to the displacement tv. The curves of equal distortion are like
that shown in fig. 17.
The surfaces into which the cross-sections are distorted are
given by the equations
2X
z- Zo = E~ {zo (2l- z 0) - {a2 } rc (d + fJy' -1),
'IT'(})
.'f:
Fig. 17.
[(~!r + (~;rr.
Then we know that h-1 dE is the distance between the consecutive
surfaces Eand E+dE, and h-1 dn, is the distance between consecu-
tives "1 and "1 +an,; also we have
h-'J = tc
2
(cosh 2E - cos 2'1] ),
from which, by putting E=Eo, we find for the perpendicular p
from the centre of the elliptic section on any tangent the equation
hab = p ...........................(129).
The solutions cf> of the differential equation, which remain
finite continuous and one-valued within the elliptic boundary, are
of the form
cp = /31'I. [.An cosh nE cos n7J +Bnsinh nEsin n7J]. .. (130).
The boundary-condition (127) at E=Eo is
h ~~ = /31p [ {(2 + u) ~ + (1 -lu) ~}cosh Eo sinh2 Eo cos 7J sin2 7J
+ lu ~ cosh8 Eo cos3 '1/ J,
or of=fJ1ab [rl(2+u)a+(1-lu)a
of/> b~} t{cOS7J-COS37J)
or A 1 =!ac [(2+iu)a+(1-!u)~J,
and As=--hac2:z: 4b3 [(2+!u)a+(l-!u)~J.
115] ELLIPTIC CYLINDER. 195
Hence
- n.81 S& a: 462 [(2 + iu) a+ (1-iu) ~]cosh 3f cos 3'1] ... (132).
Now
cosh 3E cos 811 = IC [43;2 - 12Jf-
,;. 3&] ,
0
as is easily verified; hence
rp = t.81:c [(2 + iu) as+ (1 -lu) b2]
43;2 -12JI- 3&
- n.81:c 3& + 462 [(2 + !u) a1 + (1-iu) 7J2],
or
rp = 3~~ b [a 2
2
{2a2 + 62 + 2ua1} -l (:c<J- 3tl) {2a2 + b2 + !u (a'- 62)}]
............ (133).
We notice by way of verification that the value of rp for the
circle is deduced by putting a = b.
It is easy to see that the couple (97) about the axis z vanishes.
AP, in the case of the circle we find the displacements produced
by a load X at the end z = l, by substituting from (133) in (105),
in the forms :
where
a-
1
!{4a +262+u(a -b )}
1 1
and the curves of equal distortion are now ovals lying within an
ellipse when z is small, and curves like that in fig. 17 when z is
13-2
196 THE PROBLEM OF FLEXURE. [116
considerable. These are among the curves given in the frontis-
piece of Prof. Pearson's Elastical Researcha of Barre de Saint-
Venant.
116. l"laure of J.Eolotropic Rectangular Beam.
We shall suppose the sides of the rectangle parallel to the axes
of a; and y, and of lengths 2a and 2b.
When the two principal rigidities of the beam in the planes
(z, a;) and (z, y) are unequal we have to satisfy the equation
CJ2cf> CJ2cf>
M oa.a + L oy' = 0 ...................... (136),
within the rectangle, and the boundary-conditions
of/> L 2Lu, y ~J .......... (137).
oa: = if11 [ u1a.a + E- Mu1-
when a;= ± a and b > y > - b,
ocf> E- Mrr1
and oy = f11 L a;y ..................... (138),
oy'
Ox=o
at the y' boundary.
Assume
when b' > y' > - b'. Now between these limits we may expand y"
by Fourier's theorem in the form
y'2 =b'-
3
2
+4b'
2
«> l-(-)n
r- I1 -nc 2
os-
b' .
n1J
Hence
E b2
A =-a'-tr-
2M I 3,
And therefore
(
1/> = .811 (2~ a•- itrsb') a:- l(E- Mtrl) (~- ~')
\.
4b• "'
- __.. u I
(-t sinh mrx
-.
b 'WTr'!J] }
cos --b ....... (143).
w- 1 [ n h n7ra
cos -b-
Returning to the general case, it is easy to see by symmetry
that the integrals JJ :1: ~~ da;dy, and JJ y ~: da;dy vanish, and that
therefore the terms in fJ1contribute no couple about the axis z.
Now suppose the beam bent by a load X parallel tom applied
at the end • == l, we have, as in art. 107, to add the solutions for
fJ1 == iXf(Ea•b),
«1 == -!Xlf(Ea•b).
The displacements are
X
u ==i Ea•b [(l-z) (u1ufl-u1!l)- (l-}z) z2],
X
V == i Ea•b u, (l- z) my, (144),
X
W ==I Ea•b [z (2l- z) tc- (E- Mu 1- Lus) mysfL] + 4>
CURVILINEAR COORDINATES.
in which dE stands for ~! da: + ~% d/3 + ~ d"f, and 80~> 882 , 881
are the infinitesimal rotations of the three normals at (a, /3, 'Y)
about themselves necessary to bring them into coincidence with
the normals at (a +da, fJ+d/3, 'Y+ d"f).
2
Fig. 18.
Let P be the point (a,/3, 'Y); P~> (a+ da, /3, 'Y); Ps. (a, /3 + d/3, 'Y);
and Q, (a+ da, /3 + d/3, 'Y)·
88a will be a linear function of da,
It is clear that the rotation
d/3. To find the term in d/3, observe that the length PP1 is~,
~~·(iJa/3,
and this is the rotation of the system of axes from (1) towards
120] GEOMETRY. 203
(2) in passing from P to Pt, i.e. it is the term of 088 that contains
dfJ. In like manner the term in da. is
-~ 0~ (~)da.
Adding these we get the complete expression for BOa. We can
now put down the formulre
a c)
tlh--
y,=O, =r-
Ita Cry ~ da.
This meets the normal tory at (a, /3, 'Y) in the point
1
204 CURVILINEAR COORDINATES. [121
and thus the line of intersection of fJ = const. and ry = const. is a
line of curvature on ry. This proves Dupin's Theorem, and gives
for the principal curvature 1/p1 of the surface ry in the normal
section through the line drx
l
Since u, v, w are components of a vector, the changes in them
are by (3)
87£ = du- v888 +w801
Bv = dv - w881 + uBOa · · ••· ••.••. · · · · •· .••• (7).
Bw = dw- uM1 + v801
121] STRAIN-COMPONENTS. 205
Inserting the values of the 88's from (4), the first of these
equations becomes
ou ou Ou ou Ou Ou
+ + + +
oxl ~ oyl dyl ozl dzl = oa da. ofJ dfJ ().y d!y
- v [~:a (k) dfJ-~ 0~ (~) da] +w [ ha~ (k) da -~~ (k) d!yJ
•........... (8).
Remembering that dx1, dy1, dz1 are da./~. dfJJ~, dryfha we
obtain by equating the coefficients of da, df), d!y the results
Ou Ou 0 (1) + ~h.w ().y0 (1).
0~ = ~ oa + ~hav of) ~ h;.
ou - ~hsv _o_ (!\
au= hs of) ......... (9).
oyl oa 1i:J
ou =ha ou -~haw!(.!)
ozl ory oa hs
Similar results follow from the other equations of (7).
If as in ch. I. we write e, f, g, a, b, c for the six components of
strain
..... (10).
xl, yl> zl
being the components of the bodily force per unit of
mass parallel to d'nt, dn,, dna respectively.
In the above equation ~V = drx.dPd1/hthaha. and A1 = dPd1/Ma.
and so on.
Now PA1 is the x1-component of a vector quantity whose other
components are UA1 and TA1 , viz., this quantity is the resultant
stress across the surface-element A1 •
Hence
a(PAl)d o(PAt)d o(PA1)d
0 a;+-':l- Yt+-':l- z1
X1 U'/Jl UZt
In like manner
0 ('!_~~) dz = da.dfJdty [-_!__ oT + !'! (!) + !~ (.!) _!!~ (.!\]
oz 1
1
h k, o'Y ~ o'Y It,
1 h, o'Y ~ hs oa ha I
............ (16).
p (o2v
Ot2 - Y1) = h1 a
oau + h2 aQ as +haS [ 2h~ O"f
o/:3 + ha 07 a (~1) +~ a (~1 )]
07
+ ~ U [ 2~ ~ (t) + ha 0: (~) J
+ h~ [ ha (Q- R) 0~ (ij +~ (Q- P) 0~ (k) J...... (18),
and
2
ow ) aT as oR
P ( ot~ -Zl =hloa +h~ofJ+ha o'Y +~T 2h•oa
[ a (1)
ha +~ocia (1)]
~
~=e+f+g ,
We shall however be able to obtain the equations for u, v, w
as well as equations (17), (18), (19) more directly by using the
energy-method explained in the last chapter.
Iff BWdxdydz
= JJJ[ (Px1-Po;;)Bu+ (PY1-p ~)&+ (Pz1-p 00~)BwJa.a:dydz
+ JJJ<FBu+ GBv +HBw)dS....................................... (21).
123] APPLICATION OF ENERGY-METHOD. 209
.
fffs w dadSdtr
~hA
,
where sw-_aw
ae & + aw aw •
of sf+ ... + ac & ..........(21 >·
Also we have, by (10),
(1) (1)
o8u + h1h2 o{:Ja ~ 8v + ~ha orya ~ 8w .... ...(22).
Be = h1 a;
Thus we shall have to evaluate terms of the form
ha -ag - ~o (1)awJ
o (1)oW r,_ 1Ib +pXl = Po
2
u
- ~ha [ali ot~ ..... (27).
This is identical with equation (17), and the equations corre-
sponding to (18) and (19) can be written down by symmetry.
In like lll8nner by equating the terms in ~u, f>v, f>w in the
surface-integrals we obtain the boundary-conditions
lP+mU +nT=F,}
l U + mQ + nS = G, .................. (28).
lT+mS+nR=H
The strain-equations can be found as before by substituting
for the stresses their expression in terms of the strains, or by
beginning with the expression for W in terms of strain-com-
ponents.
124. Strain-equations for isotropic solid.
In the case of an isotropic solid, the strain-equations can be
put into a particularly simple form.
124] STRAIN -EQUATIONS, 211
(~~ ;~ -~: ;;) +... + ... =(a;l a;;l- a:~ aa~l)+ ... + ...
......... (31),
where the fixed system (x, y, z) is quite independent of the
directions of the (x~> y1 , z1) axes, which are the normals to the
surfaces a:, f:J, 'Y at any point.
Thus in varying the energy we have to find the variation of
the functions ~. 'GT11 fl1 2 , 'Gra. and of such quantities as
J
- ;z (00~1 8 wl) + 8 W1~~ dxdydz
+terms containing svl.
14-2
212 CURVILINEAR COORDINATES. [124
The parts :y ( 0~ ~ W
0 1
1) and - :z ( 0~ ~ W
0 1
1) contribute only
surface-integrals, and the other parts vanish identically. The
terms in ~V1 may similarly be shewn to contribute surface-
integrals only.
JJJ ~ {! (:\ + 2p.) fl.2 + 2p. (11112 + 'GT~I +'Gras)} d~~drt... (32),
and we can obtain the equations of equilibrium or small motion
in terms of u, v, w by variation of this integral.
The term in au proceeding from the variation of the term in
!l. is
1
and this is
and this is
2vl = __
rsme orfl ior (vrain
;!._[ow- e)] ,
2v2=~ [:r(ur)-:],
2va = ~- 8k-0 [oae<~w sin e)- 0~ (ur)J.
4t. The stress-equations are
1 oP 1 oU oT 3T+(P -Q)cot e _ fiJ~t x)
r-oe+ rsin tJ orfl +Or+ r - p \at 1 - 1 ,
(A+ 21') sin 8 ~!- 21' ~~ + 21' sin 8 :r (N111) = pr sin 8 (a;- X 1) ,
a.:1 a
(A+ 21') cosec (1 o<f>- 21'fu. (r•t) + 21' aB = pr
ova (Otv
ott- Y1 ) ,
2vl = !r [ow
oe - !oz (rv)J ,
2 "~~"t
au ow
= oz- Or,
[0
2v8 =-1 - (rv)- ;;n
r or ou
Ou] •
1
For applications of polar coordinates the reader is referred to Mr Chree's
paper on ' The equations of an isotropic solid in polar and cylindrical coordinates',
Camb. Phil. Soc. Trans. XIV., 1889.
126] RADIAL STRAIN, 217
4°. The stress-equations are
aP
0r
1 au
+;:
aT P- Q
ofJ + oz +-r-=p ot'-X
((Jtu 1
) '
au+ ! ~Q + as + 2 u = P (OSv _ y )
Or r ae az r ott ' 1
T-x =- a ( U',x)
x X~ + 2J£- -
r r or r'
or x x [ (A+ 2.u) oU'
T-=- ~.
r r . -+2A-
or r
oU' U'
Thus T=(A + 2p.) ~ + 2A- ................... (45).
ur r
This is the radial traction per unit area across any element of a
concentric sphere of radius r.
We shall now consider some examples of these formulae, and,
as we do not require U to denote a component of stress, we shall
suppress the accent on U'.
127] SPHERE-PROBLEMS. 219
B
18 U=Ar+]:2·
The radial stress at any point is
hence
and
( 3:\ + 2p.) A (_!_ -
rl
_!)
r0
3
=- ~+p
r r 1
3
1
8
0 '
•
from which
P.J.r1 - Poro r + _..!:_ ro r1 ( [JJ. -Po) .! ..... (53).
3 3 3 8
U= 1
3:\ + 2p. r 0 - r1
3 3
4p. r 08 - r18 rl
In particular if p 0 = 0, p 1 = p we have a spherical envelope
strained by internal pressure. The displacement is
3 1
pr1 r
-~--s
[ 3X 1 2 + T1 r:S 0]
.................. (54).
ro - rl + p. -slJI- .,-
The radial extension at any point is
8 2
pr1 [
- a - s 3:\
1 - -1 r_.0] .................. (55),
ro -r~ + 2p. 2p. .,-
which is greatest at the inner surface.
222 CURVILINEAR COORDINATES. [129
or !!___
dw
("w- (v- rc') tan w) = 0
v-rc'+vwtanw '
224 CURVILINEAR COORDINATES. [130
where xis written for ~ea.
This equation reduces to
{x2- v (3- v)} = 0,
afJ r;ec2 x
88 = 4 c:Z
21') [(2A +a,.,) a'- (2'A. +I') r].
and a tension in the direction of the axis of the cylinder
--
1111 = tJJ'p'A. [2"A. + 31' 2 - 2-"]
4 (>.. + 21') 'A. + I' a .,- ·
These are principal stresses, and the maximum of each is at the
axis, where ;; and 8i are > -;; and are each equal to
tJJ'p
4 ('A.+ 21') (2A + 31') at.
Thus Lame's condition of safety (art. 5'1) would be that
2"A. + 31'
T.o > tJJ2pa2 4("A.+ 21')'
or if Poisson's ratio be !,
To > frpspat.
Thus the maximum angular velocity for safety given by Lame's
method is less than that given by Poncelet's in the ratio vl
The maximum difference of greatest and least principal
stresses is the value of 8i--;; at the axis, and this is
tJJ2 pa1 f.' (2'A. + 31')
4 (>.. + p.)(>.. + 21')'
On the "stress-difference" theory (art. 5'1) this must be less
than T0 • The maximum angular velocity for safety according
to this theory is .Vi of that given by Poncelet's theory, Poisson's
ratio being t·
L. 15
226 CURVILINEAR COORDINATES. [130
Let the axis be the axis of z, and let r be the radius vector
to any point drawn perpendicular to this axis, and (} the angle
between the direction of r and a fixed plane through the axis ;
also let u and w be the displacements in the direction of the radius
and the axis of z. Then the strains are
ou ' t he extensiOn
e =or . a1ong r,
g = ow
oz , t h e extensiOn
. a1ong z,
ou ow
b = oz +or' the shear of the plane (r, z)
If the material be isotropic the stresses are 2
-
rr =P=:>..~+2p.or'
ou 'I
Bi=Q=~ +2p.~.
r
.................. (85).
1 Only the leading steps of the analysis are given, and the verification is left to
the reader.
2 See art. 49.
15-2
228 CURVILINEAR COORDINATES. [131
when r=a.
J Pdz
l
-l
= 0 ........................... (88),
where n• = l {1 + 8 A + ~- F},
and we can find, as in the cylinder-problem,
U= 1 {I'I"In+f-poron+fr"-i+(~r)-PI"lt--porot-
rosn-rr (n-i)a+2F 01 (n+t)a-2F r--1}
.....• (95),
which agrees with equation (53) in the case of isotropy.
The cubical dilatation of the spherical cavity is the value of
3 Ufr when r = r 1 , and this is
3rln-f {PI"ln+f- Pof'on+f "' Plr1t-- Porof-}
rosn-rl"' (n -l) a+ 2F + "o (n +l) a- 2F ... (9G).
This result is of importance in the theory of piezometer ex-
periments, for which a discrepancy appears to ha.ve been observed
between the results obtained and the dilatation that would have
place if the material were isotropic. The solution in (96) contains
3 independent constants and Saint-Venant 1 held that these could
be adjusted so as to explain the experiments in question.
1 See Pearson's Elaltical Btltarclu• of Barrl de Saiflt-VtMnt, p. 82.
CHAPTER VIII.
GENERAL SOLUTIONS
(A.+ J£)
aoyA + J£V'v + pY = 0, ............... (1),
(A.+ J£)
aozA + J£VIw+ pZ = 0
where A is the cubical dilatation given by
ou Ov ow
A ::;:: ax+ ay + az .. ···.................. ·<2)·
In order to solve these equations we seek first any set of
particular integrals in terms of X, Y, Z, and secondly the most
genera.l complementary solutions of the same equations with X,
Y, Z all equal to zero. The first set of particular integrals
obtained will not in general lead to values of the stresses or
displacements which satisfy the boundary-conditions. In that
case we have to determine the arbitrary functions or arbitrary
constants, that occur in the complementary solutions, so that the
complete solutions, consisting of particular integrals and com-
plementary functions, may satisfy these conditions.
134. Formulm for the Bodlly J'orce1 1•
Let X, Y, Z be the components of the bodily force, per unit of
mass, supposed finite continuous and one-valued functions of x, y, s
throughout the body; we seek to throw X, Y, Z into the forms
X=of _aw +av
OX Oy oz J
y = of _ au+ aw .....................(3),
ay oz ax •
Z=of _av +au
oz ox ay
where U, V, W, and/ are functions of x, y, z.
By differentiating the above equations with respect to x, y, s
and adding, we obtain
, _ax aY az _
vI- OX+ oy + oz- ci> say ................. (4)·
1 The subject-matter of this and the two following articles is due to Prof. Betti-
Teoria della Ekuticita. Il Nuo"o Cimento, 1872.
234 GENERAL SOLUTIONS. [134
Let 4>' be the value of ell at (a/, y', z), and r the distance of
(.21, y, z) from (m, y', z), then a particular integral of the equation
for I is the potential of a distribution whose density at (.21', y', z')
1
is - 71" ell', so that we may write
4
I= ,Y- 4~ J{j~' k'dy'dz' = ,Y + Fsay ........... (5),
where the integration extends throughout the solid, and + is a.
function which is finite continuous and one-valued within the body
and satisfies the equation
vsy. = o..............................(6);
we may complete the definition of ,Y by subjecting it to the
condition
ov oF
ov + ov -(lX +mY +nZ)=O ............... (7)
at the boundary, (l, m, n) denoting the direction-cosines of the
normal drawn outwards, and dv the element of this normal. Thus
the function f is completely determined.
Now let
of
X =a$+ G,
where, as in the case of ell', G' is the value of G at (.21', y', z'), and in
like manner let
135] DISCUSSION OF THE BODILY FORCES. 235
then we have
135. Interpretation.
Consider any surface tT drawn within the body. The surface-
integral of the normal component of the system of forces depending
on f is p JJJtl>dT, where dT is the element of the volume within
the surface tr, and, when the surface is contracted to a point,
we see that this system of forces tends to vary the volume of
an element.
The surface-integral of the normal component of the G, H, K
system is
.................. (15).
Then
and
oa:- oW
oy + oV')
oz = 0·
Hence we have a. solution in the form
tP = 47r
_p_ jjj[_r da:'dy'dz'' } ............(16),
(~ + 2p.)
L= 4~p.jjj~' rlAfdy'dz'
and similar forms for M and N, where a.s before f, U' are the
137] P A,RTICULAR UiTEGRALS. 237
values off, U at (x', y', z'). Hence we can write down u in the
form
u = 47T c:
+ 2JJ.)J Jf~ cos ndu;'dy'dz'
u = __!!__
47rp.
JJJ[X'r _2(A.A.++ 2p.)
p. ~ (x' 6r + Y' ar + Z' or)] da!dy'dz'
ox ox ()y oz '
p
v = 47Tp.
JJJ[ r--
Y' x + JJ. o (x' or
2(x + 2p.)oy
Y' ar Z' or)J -'d 'dz'
ox+ ay + oz Y •
.3
tWi
w = _e_ JJJ [~ - A.+ ~' ~ (x' or + Y' ar + Z' or)] dx' dy'dz'
47TJJ. r 2(X + 2JJ.) oz ox oy oz
...... (18).
Solutions equivalent to these are given in Thomson and Tait's
Natural Philusophy, Part II. art. 731.
0
()a; {(X+ 2p.) V2t/J + pf} = 0,
(X+ JJ-)
aA
ox + p.V'u + pp'u + p ofa: = 0 ............ (30),
0
where f satisfies V'f = 0, and has the time-factor e•Pt, we may put
_ 1 of _ 1 of 1 of
U-- lf ox' V--pt oy' W=--ps OZ ••••••.•. (31),
then these make
V•u=O, V'v=O, Vtw=O, A=O,
and we have a set of particular integrals.
&v G X [
011 = 2Jl- + ltiT 1 - m 2Jl- A -n•~o J ............ (34),
ew = HJl- -ltiTs+mtiT -n Xp. A 1
011 2 2
where (l, m, n) are the direction-cosines of the normal (dv) to the
boundary drawn outwards from the space occupied by the solid.
Thus we have to find u, v, w to satisfy equations of the form
14!1] THE DILATATION. 241
Vllu =a given function of (x, y, z), and ~:=a. given function at the
boundary.
Now Pro£ Betti has shewn that we can find the value of d,
at any point (uf, y', z'), so that the surfa.ce-displa.cement.a may
be given functions, if we can find systems of displacements
E, 'TJ, ~ which become equal at any point (a:, y, z) of the surface
()7·-1 ar-1 ar-1
to - ~ ' - oy , - oz ' where r is the distance between the points
(x, y, z) and (x', y', z'); and we can find d so that the surface-
tractions may be given functions if we can find displacements
(E, 'TJ, ~) auch that the surface-tractions that would produce them
are those that would occur if near the surface the displacements
ar-1
were - ox , - ar-1 or-1
ay , -az ; and he has given similar methods for
the determination of •1• •2• •a·
To calculate JJ{F' t 1
+ G' a;; +
1
H' ~~~} dS'
we may take the origin at (of, y', z'), and the surface S' a sphere
of small radius, whose centre is the origin. Then, remembering
that the normal to S' must be drawn towards its centre, we have,
by (15) of art. 29,
F' or--1 + G' Or--1 + H' Or-1
oa: oy oz
=~
r'
[u ~ + 2 auOr +
,.
P.
r oa:
P. ~ (av _ou) oz _aw)J
oy _P. z_r (ou oa:
+ two similar expressions,
=U+2p.(~Ou+l[O'V +:ew)
r' r' rer rOr rer.
141] THE DILATATION. 248
Again
F0 u + Go v + How = - 2p.
1 1 I . [ ()
u Clr
(()r-a:c1) + v ar() {rJr-1)
\.Ty + w Clr() ((),.-1)]
()z
u:c + vy+wz
=- 4p. r' .
Thus equation (37) becomes
J
where F, G, Hare the given surface-tractions.
142. Determination of the Rotatlon1.
To determine the rotation v 1 = i- (~:- ~), we take
or-
1
u, = Oy + Es· Va
or-1
As before, take (:rl, y1, s1) as origin, and 8' a small sphere
described round this point as centre, then, by (15) of art. 29,
Fl or--1 - Gl or--1 = 2p. ou- ~ ~) + e (~ -au)
(1Lrsar
oy ox rsar r ox ay
_ P. ~ [~ (aw _ov\ + ~
r r oy oil r
(auoz _aw)
OX
+ ~ (av _ au)]
r ox oy ,
Hence y- Gl ()r--1
F' ·ar-1 (F I G I HI )
i)X- aU+ at1+ aW
0
= 2p. [{l!.r ?~Or _uori (Ji.)}
r _{!!!ror ~ (!!!)}]
~ _" Or r + r~ (~om _ ou)
oy
+r~ [~r {0 (vz)-
oz
a(wz)) +~{a (wz)- 0 (uz)} +~{a (uz)- 0 (vz)}].
oy J r OX oz r oy OX
1
The integral of the last line over 8 vanishes identically.
The first line is
{o
2p. (ury) _? (vrx)~ + ~ (ov _ ou).
r' or or J r ox oy '
and, working as before, we find for the surface-integral the
value of
- 4'1Tp. (~- ~)
1
ox ay
at (x', y z'). ,
1 1 1
Hence at (x y , , Z ) we have
1
- &rp.Jf[(Fa+La)u+(G, +M,)v+ (H,+N,)w] d8 ... (44).
If the surface-displacements be given, we have to find E., 771 , ~.
a system of displacements to satisfy the general equations of
equilibrium and to make
ar-1 ()r--1
Ea=--ay, 'Ia== ox , ~1 =0 ..............(45)
246 BETTI'S METHOD. [142
at the surface ; then we calculate two sets of surface-tractions, viz.:
L 81 M1 , N 8 corresponding to f,, "'a.
t"a and F,, G,, Ha corresponding
• Or-1 Or"-1
to dJBplacements aij ,- ~ , 0.
When this is done the rotation 'GJ'a can be expressed in the
form
1
'GJ'a =- -Jf[(F,+L,)u+ (Ga+M,)v+ (H,+ N,) w] d8... (46),
87rp.
where u, v, w are the given surface-displacements.
If the surface-tractions be given, we have to find f 1 , .,,, t"a a
system of displacements to satisfy the general equations of equili-
brium, and to make the surface-tractious
F=-F, G=-G, H=-H1 •••••••••••• (47),
where F11 G,, H, are calculated as if the displacements were Ty,
ar-1
ar-1
- oa: ' 0; then the rotation 'GJ'a can be expressed in the form
{)lr-1
Ho=-2p. oz'
The system L 0 , M0 , N 0 is obtained by putting u' =Eo. v' = 'I'Jo,
w'=to· We get
I 4p. OSR-1
.:1 = X + 3p. OZ2 ' 'rill
I
=
(
1 + XX++3p.
,.,. ) 'jJtR:-1
oyoz I
I
'rllt =-
(1 +X+3p.
X+ ,.,. ) ozew··
Jt-1.
'jJt
This gives the value of .11 at (a/, y', z') in terms of the given
surface-displacements u, v, w.
145. DetermiaatioD of the cU.splacemeab.
We may now determine u, v, w at (x', y', z').
Let L, M, N denote the functions
where r is the distance from the origin to any point of the solid
Since the only (x, y, z) that occurs is the origin, we may suppress
the accents on (:c', y', z') and write
dm X+p. zx
u = 271' A + 3p. r' '
dm X+p. zy ............... (16).
'IJ=----
271' A + 3p. r' '
w=dm+dm;>..+p.~
2,.,. 2,. x + 3p. r
If dm be regarded as a small finite quantity the depression
near the origin is very great, and we must regard the origin as
excluded from the part of the solid whose deformation we in-
vestigate. The problem is that of a considerable depression near
a single point, and the above formulre shew how to find the
displacements at a distance from the point.
147] FIRST TYPE OF SIMPLE SOLUTIONS. 253
147. Elementary Dlacu.don of Particular Ezample.
Simple Solutloas of Fint Type.
On account of its importance we shall consider this solution
a priori.
It can be readily verified that the displacements
u =-a:;a:JJ
()t
dy', v =-a!~ JJp1rdrC dy',}
p1rda:'
2 ... (19)
w= - oz• JJp1rda:' dy' + ~:: V• JJp1rda:' dy'
is a solution, r being the distance of (a:, y, z) from the point (a:'y')
on the plane z =0, and p1 any function of a:', y'. Now we may
regard p1 as the surface-density of a distribution of matter on the
plane z = 0, and then p1rda:' ff dy'
is the "direct potential " of this
y,
distribution at (a:, z), and, since V'r = 2/r, tV•ffp
1rda:' is the dy'
"inverse potential" (i.e. the ordinary gravitation potential) of this
distribution.
J10 PlqdqdfJ'
2r foe
0 + (1 qt)!
= f7rPl·
A+ 2,u.
Hence H = A+ p. 4,u.7rPI ..................... (22).
The displacement at the surface is easily seen to be purely
normal and equal to
~ :3;JJ p~a.r;'dy' ......................... ( 23).
Now suppose PI to vanish at all points except near the origin,
and suppose that near the origin p1 becomes infinite in such a way
that ff p1dx' dy' is finite and equal to ,A + ,u. W .
"'+ 2,u. 47rp.
Then the part near the origin suffers a very great normal
displacement, and the resultant normal traction is W. If to fix
ideas we consider the plane z = 0 horizontal, and the axis z drawn
I49] NORMAL DISPLACEMENTS. 255
vertically downwards into the solid, the problem is that of finding
the deformations produced in the interior by very great normal
pressure distributed over a very small area so as to have a finite
resultant, and such tangential traction as will hold fixed the
parts of the bounding plane at a distance from the origin.
To obtain the displacements in this problem we have to
multiply the expressions (I7) by ...'A. + p. W .
"'+ 2p. 4'IT'J'
149. Weight supported at single point. Rest or 11111'·
race fb:ed.
The displacement can be analysed into:
. I displa.cement equaI to ..----
(I) a vert1ca W -,
'A. + 3p. -:;:----- I
"'+ 2J' -r'IT'J' r
(2) a radial displacement from the origin equal to
"A.+p. w
cose
'A.+ 2p. 47rl' 7 ,
e
where is the angle between the radius-vector and the vertical.
The stress exerted across any horizontal plane by the maUer
above it can be reduced to:
(I) a vertical pressure equal to
tt W cos
27r ('A. + 2p.) ,..
e(I + 3'A.+I' p. cost e) ,
(2) a radial tangential traction outwards from the axis z equa
p.W sine( "A.+p. 11 )
to
27r ('A.+ 2,U) r t I + 3 7 COB e.
At the surface these reduce to a radial tangential traction
p. w
'A.+ 2p. 27rr'
at all points at a finite distance from the origin. This is the
traction required to hold the surface fixed.
To find the strains we refer to polar coordinates (e, q,, r).
The displacements u', v', w' along the meridian, the parallel, and
the radius-vector are
, __ "A.+ 3p. w
sine ,_ ,_ cos w"·) e (
u - 'A.+ 2p. 47rp. r ' 11 - 0' w - 27rp. r .. · 2or •
Then, using the formulm of ch. VII. art. I25, we find that the
256 SOLID BOUNDED BY PLANE. [150
extension along the meridian is equal to that along the parallel,
and either of them is
l(X+p) Wr-1cos8/{7rp(X+2p)} ............ (25).
The contracticm along the radius vector is
i Wr-1cos 8/{JL7r) ..................... (26).
The cubical oompression is
t Wr-1 cos 8f{7r (X+ 2p)} ..................(27).
There is a. shear in the meridian plane of amount
t Wrsin 8/f7r(X + 2p)} ..................(28).
The axes of the elongation-quadric are in and perpendicular to
the meridian plane, and the two in the meridian plane can be
obtained by turning the tangent to the meridian (1) and the
radius vector (3) through an angle i tan-1 ( 2p tan 8/(3>.. + 5p )} in
the direction from (1) tow&rds (3).
150. Geaera.llN.tioa. Particular Integral for the
Bodlly Forcu 1•
The results of the preceding example are very important. We
see that if the mass ffp1da;' dy' be very small and be distributed
with a finite surface-density over a very small area, there will be
a finite normal surface-traction per unit area near the origin 1,
equal to
47Tp(X+ 2~)( __ _r____ d "t)
, sun~ eDSl y ,
"'+JL
and vanishing surface-traction elsewhere. The displacements
corresponding to this state of things are proportional to
'i1r _ o'r _ otr + >.. + 2p V"r
oxoz ' oyez ' o.e' x + JL ·
We also found that if Pl be the density of a surface-distribution
on the plane z = 0, the functions u, v, w given by
finite resultant for a very small area, and vanishing surface-traction elsewhere
see below, ari. 16ll.
150] GENERALISATION. 257
where r is the distance of any point (a:, y, .e) from the point (a:', y')
on the surface, are functions which satisfy the equations of
equilibrium at every point on either side of the surface, p1 being
finite.
It follows from this that, if p' be the volume-density at
(a:', y', .e') of a distribution of fictitious matter, and ~ be the
"direct potential" of this distribution given by
~ = fff p'rda:' dy' d.e',
where ~ = oa:
~ + ()y
011 + aw = _____!!'__!VI
oz 'X + 1'- ()z
fff p
'rda;' d I d.e'.
'!I
Observing that V' fff p'rda;' dy' dz' = - &rrp' when r = 0, we find
2
X'= 0, Y' = 0, pZ' =g.,.I'-~+ P. p', where Z' is the value of Z at
1\,+P.
(
I
a;, '!I'I z ' ).
Thus the displacements
~
u =- ()a;()z ,
U=k,U.
P Jf.J[X'r-2(~+2,u.)oa:
~+,u. a (x'ar Y'ar zar)].3-'dy,_,
Oil:+ O!f+ a8
'.3J u.u;
w= ~+ ,u. ~ ~ + ~ + 3,u. ~ !
~ + 2JA- 81r,u. r ~ + 2JA- s.,.JA- r
If in art. 149 we write !P for W, and take the vertical to
mean the direction of the force P, all the statements of that
article apply to this case.
152] LOCAL PERTURBATIONS. 259
1152. Local Perturbation•.
It is of great interest to enquire what will be the resulting
displacements when a system of forces, which acting on a rigid
body would produce equilibrium, is applied to a small part of a
solid. In the Theory of Beams we have seen that Saint-Venant
introduced a. principle, which we have called the" Principle of the
~qui valence of statically equipollent loads". This principle states
that the application of an equilibrating system of forces to a. small
part of the surface of a solid produces no sensible strain, except at
very small distances from the part subjected to the action of the
forces. M. Boussinesq brings this principle under a more general
one which he states thus :
" External forces which produce equilibrium being applied to
" an elastic solid at points within a. given sphere provoke no
"sensible displacement at distances from the sphere which are of
'' a certain order of magnitude in comparison with the radius."
M. BouBBinesq has given several examples of this principle, and
they lead to the conclusion that the application of forces to a
small part of a solid produces, at sensible distances from the part,
sensibly the same displacements as would be produced by the action
of any other system of forces equivalent to the same resultant
force and the same couple when applied to a. rigid body. Near
the region of application of the forces their mode of distribution
sensibly affects the result, and the displacements differ finitely
from those that would be produced by an equivalent set of forces
differently distributed ; but these deviations from the kind of
displacement that depends on resultant forces and moments are
practically confined to a small space near the region of application
of the forces, and they are called by M. BouBBinesq " Local Per-
turbations".
The student will find no difficulty in proving, by differentiating
the formulre of the last article with respect to z, that equal and
opposite forces, applied at points near together, in the same
straight line, produce a.t sensible distances displacements which
vary directly as the forces, and as the distance between their
points of application, and inversely as the square of the distance
from the point of application of one of them ; and that the re-
sulting strains are directly as the forces and the distance between
their points of application, and inversely as the cube of the dis-
17-2
260 SOLID BOUNDED BY PLANE. [153
tance from the point of application of one of them. Similar
results can also be proved in the case where the forces are applied
at points near together, but not in the same straight line. Such
systems produce then displacements which can be regarded simply
as local perturbations, insensible at sensible distances from the
region within which they are applied. In the case of a long thin
wire or rod strained by the application of forces at its ends, or a
very thin plate or shell strained by forces applied at its edge, the
falling off of the local perturbations at a little distance from the
region of application of the force is likely to be much more rapid.
The particular case of a very thin plate subjected to torsional
couple has been considered in Thomson and Tait's Natural Philo-
sophy, Part II., art. 728, where it is shewn that the local pertur-
bations diminish according to an exponential function of the
distance from the edge.
1 JJ(F 1fti
Or-
1 ar-1 ar-1) da:dy ...(33),
~= - 2.,. ('A. + 1-') + G oy + H ---az
where F, G, H are the given surface-tractions at z = 0.
-(z+f+RJ {1 + 2 z:it} ·
Hence
1 _ 3 (z+z' + 2R) __2(z+z')+R_ =O
V X- R(z+z' +R)- R(z+ z' + R'f R(z+z' +R'f .
262 SOLID BOUNDED BY' PLANE. [155
Also
OSx aR-1
(J:rfi)z=a;-'
OSx aR-1
aya-z= ay ·
OSx aR-1
az• = ----a;- .
Again, consider the function y defined by the equation
+ z
=(z + z') log (z + + R)- R ............ (35).
We find av- av-
az1 = (Jz =X.
and y, a;: and ~! satisfy Laplace's equation.
liSIJ. Determination of •ublldlary displacement. re-
quired in finding the rotation•.
To find v 1 we have, by art. 142, to find displacements f 11 '711 ~1
satisfying the equations of equilibrium, and such that the surface-
tractions that would produce them are the same as if the displace-
()r-1 ar-1
ments near the surface were 0,- az , Ty ; thus we have to find
displacements f 1 , '711 ~1 which satisfy the boundary-conditions
~(~~1+~)=~:;.
~ (o~oy1 + -d'7az 1)
= ~ ayt -
((Jtr- O'r-
1
azl '
1
)
......... (36).
a.,/
-x ( ax+ ay +
oE~~ a~1 )
1
a~~~ CJtR-1
and az + 2~ az =- 4~ ayaz ······<38>·
155] THE THREE ROTATIONS. 263
Now supposing that
E1, = oE"
oy • '1/1, = &r,''
oy . '1, = or'
ay ...............(39).
the third of these boundary-conditions becomes
(OE" o'IJ" oR-
oy)+<"A-+2~} o~"
1
'"=2•
o.O'x+.e~
a4' oz
where « and ,8 are constants : for with these values we find
oE"
oa: + 0'1}"
oy + a; ozl - (a- ,8) o-x
or' = 2 ~ OZ 1 •
or'
az =(2 + ,8) ~oz' + 2z o-x
ozl;
and, when z = 0, we have
oE" oo~" + or')
"A. ( oa: + oz + 2~ or'
oz = [( 2 - a + ,8) "A.+ 2 ( 2 + ,8) ~]
o-x
oz1 ;
and since ~ =0 ! 1
, the third of the boundary-conditions (40) is
satisfied if
-a"A.+ (,8 + 2) ("A.+ 2~)=- 4p..
The other two boundary-conditions become, when z = 0,
Hence a+,8+2=0.
=2zO'x_2~o'x_()tx
'Uyozl x + p. 'Uyoz oyoz
J
satisfy the boundary-conditions. It is easy to verify that they
also satisfy the differential equations of equilibrium.
Then since
ox=-~ «x=_ox ox=~
oal om • oy' oy • oz' oz
and V•x=O,
it follows that L, M, N, 4> all satisfy the equation
osv oV osv 1
0
om'• + oy'• + oz'1 =
at all points within the solid, and a.re finite, continuous, and one-
valued functions of m', y', z'.
Now the value of ~ given in art. 153 and the values found in
the last article for w~o w1 , w, can be re-written in the following
forms:
1 o4>
~ = 27r (A+ p.) oz"
2 A+ 2p. o4> 1 o (oM oL)
1
'111 = 27rp. (A+ p) oy' + 27rp. 'am' om' - 'Oy' ' ...... (47).
2
··=- A+2p. o4> 1 o (oM oL)
27rp.(A. +p.)om' + 27Tp.oy' om'- Oy' ,
when 8=0.
We have to add complementary solutions which make
aw H aw 1 a~
o.z' = - 21-' ' and o& = 4.,. (X + 1-') o&
respectively, when .z' = 0, and these are
1 aN ~
4.,.1-' o.z' ' and 4'11' (X + 1-')'
Hence the complete value of w is
1 aN ~ .z' a~
w =- 4'11'~-' az' + 4.,.(x + /-')- 4'11'~-' a&"" ........ (52).
, 1
oq,' 1 a roM' aL')
+ +rrp.oy' \a:t- 'Oy' '
aq,' 1 o (aM' oL')
l. . . (56).
U=-
w a; w ~a;
+ 47rp.r''
---
47r(~+p.)r(z+r)
w '!! w zy
=- 47r (~ + p.) r (z+ r) + 47rp. r•' · ········ (64)·
11
W= W(~+2p.) ~ + w ~
47rp. (~ + p.) r 47rp. r'
These correspond to displacements produced by a single weight.
W supported at the origin ; for the surface-tractions on any plane
z = const. are
F = 3W z'w
27r ,.. ' 27r ,. ' " .,- 2
G = 3 W zty H = 3 ~; ......... (6S).
2
: JJ~ P1 ck'dy' .....................(66),
where ff p1 ck'dy' = W. As in art. 148 we find H = P~> and the
resultant of H over the very small surface to which it is applied
is w.
The stress at any point across any plane parallel to the surface
is in the direction joining the point to the origin (where the
weight is supported), and is, as it were, a repulsion from that point
of amount
3W ~e
27rrlcos ,
where 8 is the angle the radius-vector makes with the vertical.
If we describe a sphere to pass through the origin and the
point (a;, y, z), and to have its centre on the axis of z, and if D be
the diameter of this sphere, we shall have r2 = Dz, and the stress
across horizontal planes will be the same at all points of such
a sphere, and its amount is
3W
27rP'
We notice that the expressions for the stresses (65) do not
contain any elastic constant, so that the transmission of force
across the horizontal planes is of the same character for all isotropic
solids.
162] WEIGHT AT SINGLE POINT. 271
The horizontal displacement is along the radius perpendicular
to the axis z, and is equal to
WsinO[008
4'1T'p.r
p. e 1
-'A.+p.1+cos0 '
J
where e is the angle the radius-vector from the origin makes with
the axis z.
Within the cone whose generators are given by the equation
COS1 e+cos e= p.f('A. + p.),
this displacement is from the axis, and without this cone it is
towards the axis. When X= p. the angle of this separating cone is
cos-1 !(.J3 - 1) or 68° 32' nearly.
The vertical displacement is
w ("A.+
2p.
-4- -..,.--+coste )•
'IT'p.r "'+P.
and is always downwards.
At the surface the vertical displacement downwards is
t Wr- (A.+ 2p.)f{p.'TT'(A. + p.)},
1
and-
4'1f' (~ + p.) of the displacements of the second type given by
(62), viz.:
oX oX oX
U= OX' V=-
oy. W= OZ'
272 SOLID BOUNDED BY PLANE. [163
we obtain displacements
u--
1 ax
-----
1 oa<1>
- 4-n-(X + J') o.x 4-n-1' ozo.x'
1 ax 1 eM>
t l = - 47r(X+J') oy- 4-n-l'ozoy' ... (67 ),
w= _ 1 ax_ ~ atc;I> + x + 2~-' v•<I>
47T (X+ I') oz 4"'1' ozl 4-n-J'(X +I')
where
X= ff /'llog (z + r) £k'dy',}
<l>=ffplrdofdy' .................. ( 6S);
a.nd these are the displacements produced by purely normal surface-
traction /'1 per unit area applied at z = 0.
M. Boussinesq has given several examples of the application of
these formulre to determine the displacements produced in a solid
bounded by a horizontal plane which supports a load distributed
in a given manner.
CHAPTER X.
LAM!f:'S PROBLEM 1•
where
ou ay
~=aa: +
Ov ow
+ az ........................ (2),
which are finite, continuous, and one-valued within a sphere of
radius a, and make u, v, w given functions of position on the
surface.
We may suppose the given surface-values of u, v, w expanded in
spherical surface-harmonics, and thus we may take at the surface
fi=QD ti=GO
u= I
tl=l
.A," w= I
•=1
o.......... (3),
where A., B., 0. are spherical surface-harmonics of order n.
We seek a solution of equations (1) expressed in terms of
spherical harmonics.
Differentiate equations (1) with respect to a:, y, z, add, and use
(2), and we find
(>..+ 2p.)V 2~ =0 ........................... (4).
Thus ~ satisfies Laplace's equation, and therefore, within a
sphere whose centre is the origin,~ may be expanded in a series
of spherical solid harmonics, so that we may write
ft=GO
A= I
ti=O
o• ........................... (5),
where On is a spherical solid harmonic of order n.
1'65] SURFACE-DISPLACEMENTS GIVEN. 275
in the form
>..+p. 1I=OO 1 aon
u=--2-rt l: - 2 1 ~ ........•...... (8),
p. •=O n + ufl:
and we have similar particular integrals of the equations for
v and w.
We have to add to the particular solutions complementary
solutions, so arranged that the complete expression of
oujoll: + ovfoy + owfoz
may be identical with A. Suppose these complementary solutions
are l:Un, ~ Vn, l: Wn, where Un, V", W" are spherical solid
harmonics of order n, then we have identically
~ [oU'*1 + oVn+l+ oWn+l] _ >..+ p.~_n_ e =I.O ... (9)
oll: oy oz p. 2n + 1 " " '
where we have picked out the terms containing spherical solid
harmonics of the same order n.
Thus, if we write
oUn+~ oVn+~ oWn+l .r~
--aiC + ---ay + ----az- = T n• • • • • • • • • • •••••
( )
10 ,
"tn will be a spherical solid harmonic of order n, and
>..n + p. (3n + 1)
Vn = lA' ( 2n + 1) On ...•••••••...•.•..(11 ),
and the complete expressions for u, v, w are of the forms
u= ~ ( Un- M"rs
0
t; 1
),
v= i (V.n -
1
M.nrt o"tn-1)
oy ' ......••••.•.•. (12)
18-2
276 EQUILIBRIUM: OF SPHERE. [165
I [( U,. -Mn+~a'
0
t:+l)- A,.;] =0 .....•.•. (16),
when r=a.
The left-hand side satisfies Laplace's equation within the
sphere of radius a, and vanishes at the surface, it is therefore
identically zero. There are three such equations as (16) which
are all true identically, and it is clear that the temls of any order
n separately vanish for all values of r.
Differentiating equations such as (16) with the sign I omitted
with respect to a:, y, z, and adding, we find
Yn-l =:a: (A,.~) + :y (B,. ;;) + :z (0,. ~)· .....(17),
which determines the function Yn-~> and in like manner all the
functions 'l[r are determined Then U,. is determined from (16),
and V,., W,. are given by similar equations. Thus we have finally
_ ·~"'
u- A rA
.. [ ..a.n---;+
•=1 a
( a1 -r') 11 OYn+I]
.Wn+t_-:>_
ua:
••••••••• (1 8'
)
and
M - A+p.
n+s-li{n+1)+~(3n+4) ............... (20).
166] GENERALISATION. 277
168. Dliplacement In any SoUd
Equations (12) express in terms of rational integral functions
of the coordinates quite general solutions of the equations of
equilibrium for a simply-connected region containing the origin.
For such a region ~ can always be expressed in a series of
spherical solid harmonics such as 0,., and the displacements consist
of particular integrals of the differential equations of the form given
by (8) and complementary solutions of the same equations ex-
pressible in a series of spherical solid harmonics, and the four sets
of harmonics thus introduced are connected by the set of relations
involved in the equation 'Oufo:c + 'OtJfoy + CJwfoz = ~. Mr Chree 1
has applied this method to the determination of general solutions
expressed in positive integral powers of the coordinates, and has
obtained by this means the displacements in a rotating ellipsoid.
He has also shewn that Saint-Venant's solution of the problem of
the flexure of an elliptic beam is the only possible solution which
contains no higher power than the third of the coordinates of a
point on the cross-section.
v-- == aa:a (.A,. a"r") + oya (B,. a"r") + oza (a,. anr") '
1
M. -l "A+p.
,.- "A.(n-1)+p.(3n-2)"
Thus if .A,., B,., a,.
could be expressed in terms of the surface-
tractions the problem would be solved.
We have to calculate the surface-tractions coiTesponding to
displacements such as (25).
For this purpose we first write down the value of t:,. ; it is
A - "<.' 2p.(2n -1) M. 41~
'-~>
- ._ ;>._ + p. n Tn-1·
Thus we obtain
1 o (rf"+a
- 2n + 1 ox a"'+l cf>--ra-t
) - En atn-11
rtt&+ o (atn-1
ox ,-.-1 Yn-1) · · .(32),
1 >.. (n + 2)- p. (n- 3)
where En= 2n + 1 >.. (n - 1) + p. (3n - 2) · · ··· · · · · ··· · ·<33).
Rearranging, and picking out terms that contain surface-
harmonics of order n, we have the surface-tractions required to
produce surface-displacements An, .•. expressed in terms of the
given surface-tractions Fn, ... by equations of the form
rn rtn+~ 0 (alrl-1 )
(n- 1) An a" - En atn-1 ox rt-1 Yn-1
1 0 (r"'+l )
Fnr"
- 2n + 1 tJX aln+l cf>-n-s - p.a.n-1 = 0 ............ (34),
whenr=a.
!80 EQUILIBRIUH OF SPHERE. {l67
The left-hand side of (34) satisfies Laplace's equation within
the sphere of radius a, and vanishes at the surface, it is therefore
identically zero. We have three such equations as (34) which are
all true identically, and these hold for each value of n and for all
values of r.
Introduce two new functions 'l'n-11 ~-n-~. which are spherical
solid harmonics of orders indicated hy their suffixes, and defined
by the equations
~R.n.
1
Fr = a:~Ru
oo a,
( '")" Gr == y7Rn
oo (r)"
a; , .. (r)"
Hr == &~Rn a ...(38).
Now the surface-tractions that must be applied to the surface
r ==a to produce the system of displacements (25) are given by
(34), and may be written in forms of which the type is given by
11-
= a• [ 2n + 1 ama (r-+~ ) 1 a ( r"+l)]
am+s c/1-A-~ + 2n + 3 am ~1 a-+ 1 '
and
These are the only .A's that occur, and the B's and O's can be
found by writing y and s respectively for a;.
169] SPHERICAL FLAW. 283
169. Spherical Cavity In lnflnlte 10Ud mua 1•
By analogy with the solutions for space inside a. sphere we
may write down those for space outside in forms of which the
type is U =I {!:i -knrl :a;(~~)} ............... (45),
where v and w are found by cyclical interchanges of the letters
( U, V, W}, (a:, y, .e), and 'tnH is given by the equation
(51).
P' 1
"=2n+l
From the equations of type (49), it is easy, by the method of
art. 167, to expreBS 'tn+l and X-1 in terms of given surface-
1 Only the leading steps of the analysis are given.
284 EQUILIBRIUM OF SPHERE. [169
tractions, and thence to infer the values of U.,., V.,., W.,.. We
shall consider a particular example.
Suppose that at a very great distance there is a finite shear s,
so that the displacement, when r = oo , can be expressed by
'U = sy, 'II= 0, w = 0 .................. (52).
In this case, it can be shewn, by the above analysis, that, if
the surface r ==a be free, the only harmonics y and X that occur
are proportional to xy,
and that the only harmonics U that occur
are ul and u.;also that ul and vl are respectively proportional
to '!I and a;, while wl = 0, and that u, is proportional to ,., a:(~)'
and V,r-1, W,r-1 are the same multiples of the differential co-
efficients of xyr-3 with respect toy and z.
We thus find for the forms of u, v, w
u = A 1L + B
,.. ~ (a;y)
ax rD
+ Or' ~ (fl!!i)
ax rD
+ sy
'
x
v=A-+B-
,.a
a - +Cr'-a -
ayr ()yr'
(u;y) (u;y) ...... (53),
.
The most 1mpo rtant pomt
. concerns t he shear - + ()v
~ . !tiS• au
oy uQ;
easy to shew that the value of this, when x=O, y=O, r=a, is
15X+30p.
9>.. + 1.4p.-8........................ (55),
1
See Phil. Mag. Jan. 1892, p. 77.
170] BODILY FORCES. 285
:!
(iii) Forces arising from the attraction of the harmonic in-
equalities and derivable from a potential which is I e,.Q-+1·
2 3
1 Cf. arts. 95, 96.
286 SPHERE UNDER BODILY FORCES. [171
3g
Let us write Yn+t = W,l+l + 2n + 3 enQn+~ ..................(57),
= [ Hr1 (5X+2p.)
0
~;1 +4p.H(n+1)a:QnH+(3X+2p.).A a~;HJ
......... (6'1).
Thus the part contributed to Fr is
1
a:[Ha2 (5X+ 6p.) +.A (3X+ 2p.)] [ 1 + (n: ) EnQnH]
0
+ 2HaenQnH (5X + 6p.) a:- aen [ Ha2 (5X + 2p.) ~~+1
+ 41p.H (n + 1) a:QnH + (3X + 2p.) .A
0
~:+1] ,
of which the first term vanishes identically by (64), and the
second is obtained by substituting for r its value a + enQn+t in
(66).
Thus, collecting the terms in a:Qn+1 and
0
~;+~ ,and transform-
ing by means of the identity (26), we have for the part contributed
to Fr by the strain produced by the radial forces
R [ sx- (2n- 4) p. r~ aQ,l+l 41 • aQ'*l
aen 2 2n + 3 ()a; + p.a ()a;
- 2
5X- (2n- 4) p. H o
(QA+l)]
2n + 3 -- ,an oa: ,an+a ..........( 68>·
173] GRAVITATIONAL TRACTIONS. 289
This is the typical term arising from the spherical harmonic term
enQnH in the equation of the surface.
Since we neglect e119, e..em, and enm2, we may take H to be
gp d .
-liT (~-.f.-2p,) a' an wnte a ax
2 oQnH £
or ,a -ax
oQn+I . h urfac
ill t e s
I
e-va ue
of (68), so that this surface-value becomes
__fff!_ I
1._ [5~ + 2 (n + 5) IJ. 2 oQnH
lr~+2p, en 2n+3 a ox
_ 5~- 2 (n- 2) 1J. ...-+a ~ (QnH)] ( )
2n + 3 r--. ole ,an+s • ••• • • 69 •
ox +~J.v2u+p---ax=
(X+iJ.) a~ aYn+l 0
............. (71),
where YnH is a spherical solid harmonic of order (n + 1), with a
small coefficient of the same order as e,..
Now such a particular integral can be found by assuming that
175] PARTICULAR INTEGRALS. 291
the strain throughout the sphere is irrotational, i.e. that there is a
displacement-potential cf> such that
ocf> of/> ocb
u =ox • v = oy • w = oz •
for then A = V2cf>, and the equations can be satisfied if
(X+ 2~) V1cf> + p Yn+1 = 0 .................. (72).
Just as in (8) of art. 165 we have a particular integral of this
equation in the form
p ,a
cf> =-X+ 2p, 2 (2n + 5) Yn+l•
Thus the particular integrals u, v, w of equations (71) are given
by three such equations as
- 1 p 0
u-- 2 (2n + 5) ~+ 2}£ ()$ (r-2YnH) ............ (73).
The terms contributed to~ (the product of the radius and the
radial displacement) are easily found from (73) to be
p n+3
-X+ 2J.£ 2 (2n + 5) r2Yn+l·
Thus the terms contributed to Fr by the particular integrals
(73) are found by using the formula (23) to be
U= ~ [An~ +(a 2
-r2)Mn+t ~~+l];
+ ~ [An~ +(a 2
-r2)Mn+2 °t;HJ. ........... (75),
in which
and the typical terms in the surface-tractions are this and the
terms given in (69) and (74). Since the surface is supposed free
we must add these terms together and equate the result to zero.
We thus obtain an equation which may be written
~[an oYnH+b
0/lJ n
r2"H~ (Yn+l)+an'e oQnH+b 'e ~t~ (Qn+l)
0/lJ r2"+3 n OllJ n n OllJ r2"+&
~ow Vn+I• (~)m+s cf>-n-2, and QnH are spherical solid har-
monics of order n + 1, and we have obtained in (82) and (83)
three equations which determine these in terms of W•+~· It is
clear that, if ~ w.+l be reduced to a single term, ~e.Qn+l will at
the same time be reduced to a single term containing the same
solid harmonic, and 'o/n+I and f"'l+scJ>-n-2 will be the only Y. and cJ>
functions that occur.
178. Determination or the unknown harmonics.
We may now suppose that the disturbing potential consists of
a single spherical solid harmonic Wn+l" Then they and cJ> func-
tions are determined, and likewise the harmonic inequality enQn+I•
and we seek to determine the unknown harmonics ~.A.., I.B.,
~c •.
From the equation (78) pick out the terms containing spherical
solid harmonics of order n, and of order n + 2. We find two equa-
tions
- f£ (n -
1) A r" 0Yn+I
n a" = an ax + an En
oQn+l
I
ax 0 II (rm+s
+ an OX am+l c/>-n-2 '
)
+bn "r2n+G~(Vn+I)
()a; rl"+a + ,..
u(n -1)A n a" = 0• r"]
when r = a ; and just as in art. 177 we find
X (n + 4)- p. (n- 1))
Vn+I { n + 1 + (n + 2) X (n + 1) + IJ. (3n + 4)f
1
= _n + ~ p (2n + 5) >.. + 2 (n -1- 2) IJ. W:
2n + 3 X + 2p. p. n+IJ
(!)m+a
a
cf>
--~
= _ n + 1 2n + 1_p_ >.. + p. (n + 2) W
2n 2n + 3 X + 2p. p. n+I
......... (85).
The functions An, AnH may then be written down by means
of equations (84).
In connexion with this problem we may notice in particular
179] ORAVITATION ANNULLED. 297
En
Yn+1
Qn+I = a [ 2n +5 -
(r)m+a
a 2n cf>--n-~
+ 1-
p n +3
X + 2p. 2 (2n +5)
W J
fl+l •
where
'Q 1 _ 1 (X + fJ-) (n + 1) ( 2 2) .lA
E'n n+I-2r0 X(n+1)+p.(3n+4) a -ro Tn+I
[ Yn+t (r)m+a
+ro 2n+S- a ra 0
2
2
cf>-n-t p n+3 W
2n+i-X+2p.2(2n+5) n+I
J
and we are to give to each spherical harmonic function that occurs
its value when r = r00
298 SPHERE UNDER BODILY FORCES. [180
Substituting for YnH and ~-n-2 their values from (85), we find
that the height of the harmonic inequality is given by
_ ap WnH n+ 1 (2n+ 3)X+(2n+2)p.
enQn+l- p. 2n (2n2 + 8n + 9) A+ (2n2 + 6n + 6) p. • • ·<86),
and the radial displacement is given by
e 'Q I- ro{J w!'+l n + 1 [ (2n + 3) A+ (2n + 2) Jl.
" nH- p. 2n (2n"+8n+9)X+(2n1 +6n+6)p.
+
a2 -r0•
-r::- (n + 1) (2n + (n+3)A+(n+2)JJ.
+ 9}i + (2n +
2 8n + 6)
J 2 6n p. ... (
87 )·
a.,
= _ _qpe..a}__ __ {5A. + 2 (n + 5) #£ _ 3 A+ (n + 2) #£}
(2n+3)(A.+2p.) 5 2n+3 '
1
fJ = _ _ gpe~- { (2n+ 5)A+(2n+4)f£ _ 5A- (2n- 4) #£}
3
:::=n~?K~~~>;~:;~;~+5) 5 .I
......... (91).
The three equations such as (90) give
n+l A(n+4)-#£(n-1)}
{n+ 2 + A(n+ 1) + #£(3n + 4) f'-+n+ 1
= 2n + 5 gpe..Q,.+t {a (2n. + 5) A.+ (2n + 4) p. _ 5A.- (2n- 4) #£}
2n+3 A+2p. (2n+3)(2n+5) 5 '
2
2n : 1 p. (~Y'~+a
cf>--s I
= (n + 1) gpe,.Qn+t {5A. + 2 (n + 5) p. _
3 A+ (n_-!:: 2) P.-} I
2n + 3 A+ 2J1- 5 2n + 3 )
.••••.•.. (92).
300 SPHERE UNDER BODILY FORCES. [180
- p. (n - 1) A n :t!
an = ,.
-n
- oQn+~ ~ (rt"+s""'t'-fi-'l)
ox + txn"OX }
rn+s
- p. ( n + 1) An+2 an+s = r2n
+a {QfJn OX
o (Q"+I) Q , o(Vn+I)}
rsn+s + fJn OX r2fl+3
(93).
H =n (X+g;IJ,) a'
A=- 5)..+6.u. -h gpa
3A. + 211- X + 2J.1- '
5)..
Ea = f 2).. + 7II- '
X+3p. ..••••......... (96),.
~=-!--
).. + 211- pa2,
5X + l2JJ-
~
I
b 1 _ - _L 5).. + 2JA-
l- n )..+ 211- gp
where the terms in a~ have been rejected from H and A as an
p. 289.
We shall consider the particular cases where 11-/A = 0 and
p.j).. = 1.
302 SPHERE UNDER BODILY FORCES. [182
182. Incompreutble material.
When the material of the sphere is incompressible we have X
infinite and the constants become
.A =0, H =0, Ea=-h_,
tlt = - }pa2, tlt' = }gpa2, b1 = }p, b1' = - }gp.
Hence the first two of equations (95) become
¥ 1-'"h = v p ( wll- fgelQt),
"' (~r c~>-~~ =- te <wll- tgelQll>;
and the third of equations (95) then gives
apWll
e1Ql = 2gpa 191-£ •
-5-+5
We shall write this result
EtQll =
15~
6~ + 19 gwll ..................... (9'1),
where ~=igpaf~-£,
Hence
y,- (~) 4>-a=» ~ W,- H- ~ e1Qs
7 a 3 ag a '
so that the third of equations (95) becomes
~ W:
-9t- ~ (W•
~ Q,-f- Q) ~) e1Qs
s = ( 1-25
· ag ,-ty.-el
a a -+tel
g - a.
From this equation we obtain the height of the harmonic
inequality
225~
elQ, = 275 + 93~ g
w, .....................(98).
The result (97) of the last article is that when the material is
incompressible the height of the harmonic inequality is
225~ w2
e1Q2 = 285 + 90~ g
The difference between the result which holds when the
material is incompressible and that which holds when the material
fulfils Poisson's condition is a very small fraction of either for any
the same value of~. so that in case it is uncertain which hypothesis
is the best to make no very great error can arise in our estimate
of the harmonic inequality if we assume the material to be in-
compressible.
In the applications that we shall make to problems relating to
the earth, considered as an elastic solid globe, we shall have to
assume the material incompressible to avoid the difficulties ex-
plained in art. 127.
r_ g
-a [ 1 - fP2 ora
1
+ ¥t :pa] ...............(99),
so that the ellipticity of the surface is
(J) a 1
2
i- ....................(100).
g 1+¥~
gpa
For a liquid sphere the ellipticity would be iro'afg.
If the globe be of the same mass and diameter as the Earth,
and of the rigidity of steel or iron, we have
a= 640 x 108 centimetres,
p=5"5,
p. = 780 x lOS grammes' weight per sq. centimetre 1•
The ellipticity of the surface due to the rotation is diminished
by the rigidity in the ratio
780
1 : 1 +.Y .1f640
or nearly 1 : 3.
If the rigidity were that of glass,
p. = 244 x 106 grammes' weight per square centimetre;
244
and the ratio is 1:1 +.Y lf640
or nearly 3 : 5.
The same numbers apply generally to a globe of the same
mass and diameter of the Earth whatever may be the forces whose
potential is W,, and we have
~ = J nearly for a rigidity equal to that of steel,
~ = 5 nearly for a rigidity equal to that of glass.
1 + ¥ __E:_
gpa
of the true equilibrium height, the material being incompressible.
This ratio is about ! when p. is the same as that for steel, and
about i when p. is the same as that for glass.
and
ou Ov Ow
A= ox+ oy + oz ........................ (2),
1 The following among other authorities may be consulted :
Jaerisch 'Ueber die elastisohen Bcihwingungen einer isotropen Kugel'. Crelle's
Journal, LUXVIII. 1880.
Lamb 'On the Vibrations of an elastic sphere'. Proc. Lond. Math. Soc. xm.
1882, and ' On the Vibrations of a spherical shell'. Proc. Lond. Math. Soc. XIV.
1888.
Love 'The free and forced Vibrations of an elastic spherical shell ... '. Proc.
Lond. Math. Soc. XIX. 1888,
Chree ' On the equations of an isotropic elastic solid in cylindrical and polar
coordinates'. Camb. Phil. Soc. Tram. XIV. 1887.
Rayleigh 'On Waves propagated along the plane surface of an elastic solid'.
Proc. Lond. Math. Soc, XVII. 1886.
310 vmRATIONS OF SPHERE. [190
which are simple harmonic functions of the time, are finite,
continuous, and one-valued within the boundary, and satisfy the
condition that the bounding surface is free from stress.
Suppose the solid performing free vibrations whose period
is 27r/p; then for a;; ... we may substitute - yu ..., and thus the
equations (1) become of the type
oA
(11. + p.) OX + p.VIu + pp"u = 0 ................. (3).
r:;s + " 2
) w = (1 - K'fh') ~~
Equations ( 6) can be satisfied by putting
1 oA 1 oA 1 oA
u = - fii OX ' v = - Jl2 Oy ' w = - fii oz ......... (7),
where A satisfies (5), and these satisfy (2).
Hence the complete solutions of the equations of vibration
consist of the sums of these solutions and the general solutions
of the equations
(V1 + K')u =0, (V 1 + K')v=O, (V' + "') w= 0,}
ou Ov Ow ...... (8).
-+-+-=0
ox oy oz
191. Description of method
Before proceeding we point out the kind of results to be
obtained. According to the theory explained in art. 79, there will
be an indefinite number of normal modes of oscillation, and the
oscillations of any normal mode can be executed independently.
191] NORMAL FUNCTIONS. 311
If the system be oscillating in a normal mode then at any instant
the displacements can be expressed in the form
u= u'.A cos(pt+ e), v =v'.A cos (pt+ e), w = w'.A cos (pt +e),
where 2w/p is the period, .A is a small arbitrary constant, and
u', v', w' are functions of X, y, z. These functions are called
normal function8, and the determination of the .vibrations of' any
elastic system is effected when the normal functions are known and
the frequency-equations have been formed and solved. In what
follows we shall first determine the forms of the normal functions;
and no confusion ought to arise if we denote them by u, v, w,
instead of u', v', w', and write A for 'iJufiJx + 'OtJfoy + ()wj(Jz, where
u, v, w are simply normal functions. In strictness each term of the
cubical dilatation also contains a factor of the form .A cos (pt+ e).
Among the vibrations of a sphere we shall find that for some
modes there are spherical surfaces at which the displacement
vanishes, just as in the vibrations of a string there may be one or
more nodal points. Such surfaces will be called nodal surfaces,
and their number and position are determined by the type of
vibration and the frequency, and, conversely, if the number and
position of these surfaces be given the type and the frequency are
determinate. We shall find also other modes for which there are
no surfaces at which the displacement vanishes, but there will
then be surfaces at which the radial displacement vanishes,
and we shall term such surfaces quasi-nodal. The number of the
quasi-nodal surfaces for a particular class of vibrations does not in
general determine the frequency or the type.
We proceed now to the consideration of the vibrations of an
isotropic elastic sphere.
192. Determination of the Dilatatibn.
We have to find a solution of the equation (V 1 +h'}A=0 in
a form adapted to satisfy boundary-conditions at the surface of
a sphere. We therefore suppose A, at the surface of the sphere,
expressible in spherical surface-harmonics, and we treat the typical
term A=RnS", where Sn is a spherical surface-harmonic, and Rn
is a function of r, defined by the equation
a• n(n+l)
o(hr)'(rRn)+rR,.- (hr'f (rRn)=0 ......... (9).
This is the case of Riccati's equation which is integrable in
312 VmRATIONS OF SPHERE. [192
terms of circular functions, and the solution which remains finite
in space containing the origin is
Rn =,... (~ ;r)" (s!nr~) ..................(10).
This function can be expanded in a convergent series of powers
of r, beginning with r'l, and, if we take such a multiplier as will
make the first coefficient unity, and write ,...B.= m,., we shall
have as the general form of A
a=ao
A= I ro,.,Y,.(hr) ..................... (ll),
•=0
where ro,. is a spherical solid harmonic of order n, and
,y,.(x) =(-).. 1.3. 5 ... (2n + 1) (~ :fxr (m:x) ..• (12).
We add a few properties of the functions V• (x) which admit
of ready verification :
The equations connecting consecutive y's are
d a;l
x do; '1/rn-1 (x) =- n
2
+ 1 ,Y,.(x) =(2n-l){,Y,._,(x)-y,._1(x)} ...(13).
The differential equation is
d"ts~) + 2(nx+ 1)dydx(x) +V",.(x)=O .....•... (H).
The series for y,. (x) is
2~
,Y,.(x) =I- 2. + 3 + 2 .4.(2n :3)(2n + 5)- ··· ( 15);
and thus 'tn (x) = l v'(2'11') 1. 3 ... (2n + l):r<n+ll Jn+l(x) ... (16),
where Jn+l (x) is the Bessel's function of order n + !.
The function y,. (z) of the complex variable z is a uniform
function in all parts of the plane of z which exclude the point
at infinity. This point is an essential critical point .of the
function.
193. Determination of the Displacement..
The forms of u, v, w, which satisfy (8 ), can be written down in
the same way ; thus
u=IU,.,Y,.(IC1"), v=IV,.,Y,.(IC1"), w=IW,.y,.(IC1") ...... (17),
where U,., V.u W,. are spherical solicl harmonics of order n, and
these have to be arranged to satisfy the condition
Ou ov Ow
ox +oy + oz =O ·····················<
18>·
193] DISPLACEMENTS. 313
Now this condition is
I [-t (..cr) (oU~
n ()a; + ~V
oy.. +oW,.)+
oz CJ't,.(..cr)
Or
(xU. . +yV.. +zW"')]
r
= 0 ....•.... (19).
If Xn be a spherical solid harmonic of order n, then the
forms
u.n = yax.. ~ v: _ ax.. ax.. w.___ ,.ax.. _ Y ax..
oz - z oy ' n - z 0/JJ - llJ a; ' .. ... oy 0/JJ
satisfy the equations
xU,.+yV,. +zW,.= 0, and oU.. fox+ oV,./oy + oW,.JOz= 0,
and U,., V ,., W,. are spherical solid harmonics of order n.
·
Agam ·r u.n
1
be ~ a (4>-~) and v n and
rm=l '
ocbn+l - a,.,.an+1 oa; w.n be
similar expressions with y and z respectively for u;, these will be
spherical solid harmonics of order n provided 4>n+I be one of order
n+ 1, and we shall have
au.. av.. aw.. < )
0/JJ + oy + Tz""=n 2n+ 1 a..4>-l.
and xU. . + y V,. + z W,. = (n + 1) 4>..+I + a,.'f'ln4>n-1 •
Thus the terms contributed to (19) by such functions U,., V,11
W,., will contain 4>n+I multiplied by
+ 1 Gyn
n-r- a..;. .+l + ( n + 2) ( 2n + 5) tln+s'tn+z•
Tr + tln+s (n + 2) r---ar-
where Vn is written for 'o/'n(K.r).
By using (13) the multiplier becomes
[- (n n:
+ 1) 2
3 + (n + 2) (2n + 5) 12n+s] Vn+I·
'l'his vanishes identically if we take
n+1 JC2
12n+t = n + 2 (2n + 3) (2n + 5) ·
Thus we have found solutions of (8) in the form
u =I [+. (..cr) (a4>n+1
ax + Y axaz.. _ z ~)
ay
n +1 ~+o o(4>n+t)]
- n + 2 (2n + 3)(2n + 5) Vn+s (..cr) ou; r-t• · · ·· .. ( 20),
where v and w are to be derived from this by cyclical interchanges
of the letters u;, y, z, and Xn and 4>n+1are spherical solid harmonics
of orders indicated by the suffixes.· This solution contains two
3lt VIBllATIO:SS OF SP'IIEBE. [193
unknown spherical hannonic functions, and therefore constitutes
the general 80lntion of 1ihe form required.
The complete expressions for the nonoal fonctioos will be found
by adding together the left-hand sides of the equations such as (20),
and the particular solutions given by (7) and (11). They depend
upon three sets of unknown spherical solid hannonics. e, x, ••
and we sball shew bow to determine the ratios of these from the
boundary-conditions.
For convenience of reference we state here the results 80 far
obtained. The cubical dilatation is
~A••+• (Jar) 008 (pt +E).
The displacement u is
~.A oos(pt+E) [- ~s!: {•a'ta(Ar)+ 't.(~) (y~• -s~•)
++-~ (~) U: -+•+1 (~) n:j (!n+~:+3)! (,!:~)]'
and the displacements " and to are to be derived by cyclical
interchanges of the letters a:, y, z. The summations extend to all
integral values of n, and to all values of p given by the frequency-
equations with the corresponding values o{ I& and "·
It is also convenient to state that the product o£ the radius-
vector and the radial displacement is given by the equation
r = u.a:+ "Y+ toz=- ~!.A cos{pt+ E)[{n-t.(l&r) +hr+.' (l&r)}•..]
+~A oos(pt+E)[n-t.(~) ..].
l
on the dilatation, Cor which, omitting the constant and the time-
Cactor,
~ = !.o>,.y.(hr),
u=-_!_o~ "= _ _!_a~ to=-_!_o~ ........(22).
h2 oa: • h• ()y • h2 az .
Using the identity
xj(a:, y, z)= 2n: 1 {~ _,-H ;a:(,L.)} ......(23),
19~] BOUNDARY-CONDITIONS. 315
we find
" -~ ["'"
a;~=.... 2n (hr)
+ 1 {OOJn 0 ( rm+~
()a; - , -H ()a; OJn )}] .........(24).
The terms of ua: + vy + we depending on a are - ~~ ~~
or - !~ IQ).. {nt" (hr) + r dvd;hr)} ............(25).
The terms contributed hereby to a: ( +ua; vy + wz) are
I [ Pn ( y ~-zaxn)
oz ay
+a aMn + b rsn+a! (~)
n oa: " aa: ,-+1
1 K' [ 2(n + 2) ,
bn=- h2 2n+ 1 V'n(ha)+ K'a 2 havn (ha) ,
J ... (33).
following article are reached the reader is referred to Prof. Lamb's paper in Proc.
Lond. Math. Soc. XIII., 1882.
199] P.ll\TICULAR MODES. 319
The first six roots are found to be
"a/'11' = ·7961,
2·2715, 3·3469, 4·3837, 5·4059, 6·4209.
The character of the vibration depends on the form of Xt· In
the case where X2 is the zonal harmonic 2z2 - afl- '!f, we have for
the normal functions
u=t~("r)yz, v=-V's("r)xz, w=O.
All the particles on the same parallel move along the parallel
through a small distance proportional to the sine of the latitude,
and the equatorial plane is nodal.
{vi (ha) +~~~hay/ (ha)} ;~I + l {vi (~ea) + ": y/ (~ea)} cf>I = 0,
and the frequency-equation, obtained by eliminating (I)I/h' and ~.
and supposing "= .j8h, is
4 'it, ("a)
1+ .\73 I J3 + ~t~j~~ = 0.
• r~
~ea'Yl
("a)
.;a "a'i'I (~ea)
where Q)" and n,. are spherical solid harmonics, and y.. (hr)
and 'I'" (hr) are defined by the equations
~ [<n+1)Vta+I(~a)cf>tt+l
- {(n + 1) "o/tt+1 (ha) +a d"o/ta+I
da(ha)} ~~,~-
7
J
0ta+I n + 1 W:ta+I e<.Pt ••• (53),
so that this is of the form Eta!;+I e.Pt, where Eta is a. number; and
the equation of the surface of the sphere at time t is
_ E. w'*l ~
r-a+ ap tf'. ••••••••••••• , •••••• (54).
......... (55).
Ct =- {~al v~ (~a)- 2'frl (~a>} 1
We give the arithmetical result for the special case when the
material is incompressible and so rigid that the fourth power of
tea may be neglected.
We have by (15)
ka = 0, V1 (ha) = 1, Vs (ha) = 1, 'fr,' (ha) = 0,
te"'a1 rat ~'
V'~(~a)=1- 10' Vs(~a)= 1-14 + 504'
o,=2(1-"';)~
2r { ra' 8 ( rat ~a')}
~ =- 15 1 -14 + tt'a' - T + 126 •
The height of the harmonic inequality becomes, when we reject
~' in the numerator and denominator of (56),
2W.
ap'
[-1 + -H + ~tasMJ
-H+tt'as-flg
tt'Pt,
-fi pa
,.,. Wstf'Pt ••••••••••••············(57),
which agrees with the result of the corresponding equilibrium
problem in art. 182.
328 VmRATIONS OF SOLID. [203
For a sphere of the mass and diameter of the earth, and of
the rigidity of steel or iron, executing vibrations of the species
considered with a semi-diurnal period, we have, in C.G.s. units
2w/p = 12 X 60 X 60, p = 5·6, p. = 800 X 10', a= 64o() X 1()8,
so that "a= t nearly.
It follows from this that the neglect of ("a)' would be fairly
justifiable in the case of such a body. We conclude that in the
case of an elastic solid earth the bodily tides would follow the
equilibrium law.
where ,-2 = / 1
+9 2
- h1 . . . . . . . . . . . . . . . . . . . . . . . . (59).
Supposing r real and positive, we have
~ = Pe-n ........................... (60),
where P ex: e•Cft+n+pt); and the particular solutions (7) become
if tg r
u =- h• Pe-n, v =- h' Pe-n, w = hs Pe-n ...... (61).
The complementary solutions satisfy equations such as
{)'.lu au av aw
a.zs- h = 0, ax+ ay + az = 0 ............... (62),
where s'=/'+g2 -r ........................... (63);
J Proc, Lond. Math. Soc. xvu. 1886.
203] SURFACE-WAVES. 329
so that we have
u=A6-u, v=Br", w = 06-81 .........(64),
where ifA+ ,gB -sO= 0 ...................... (65),
and A, B, 0 contain a factor 6d.ta:+n+ptJ.
Hence the expressions for the displacements
has one real root "' = ·9127 5, while the complex roots make the
real parts of r and s have opposite signs, so that they may be
rejected. We now have
"
2
= •91275 (f" + gt), r" = p + g", at= '08725 (f" + g•),
and hlu =if (- 6-n + ·54336-a) e•(pt+fz+n> }
h"v= '9(-6-n + ·5433r") 6'(pl+fz+n> ...... (71).
h"w = v(/ 2 + g1)(~ -1·840r") 6'(pt+j!r+gy)
For progressive waves whose fronts are parallel to the axis of
y, we have
u = U (rfz- ·54336-") sin (pt +fa:),} ......... (7!),
w = U(6-fz-1·840ru) cos(pt +fa:)
where U is a constant ; and the velocity of propagation is
plf = •9554 .J(p.fp),
which is slightly less than that of waves of distortion m an
330 WAVES AT SURFACE OF SOLID, [203
unlimited medium. The horizontal motion vanishes at a certain
depth. The motion at the surface is given by
u = ·4567 U sin (pt +fa:)}
o..t.O U cos ( p t +fi)
w = - 0"2'
0
a; .................. (73) j
so that the particles move in elliptic orbits whose axes are nearly
in the ratio 2 : 1.
Lord Rayleigh also considers the cases where X= p., or the
material fulfils Poisson's condition, where A= 0, or longitudinal
extension is unaccompanied by lateral contraction, and where
X= - JJ£, or the bulk-modulus vanishes. For A= J£ he finds
It'= •8453 (jt + g1), rt = ·7182 (jt +!f), ,. = •1547 (/1 +!f).
For a progressive wave
u = U (e-n- ·5773e-") sin (pt +fa:), }
W = U (•8475e-n- 1•4679e-") COB (pt +fa:) ...... (H),
and the ratio of the axes of the elliptic orbit, described by a
surface-particle, is reduced to about f.
Lord Rayleigh suggests that these surface-waves may play an
important part in earthquakes and in collision, as they diverge
from the source of disturbance in two dimensions only, and con-
sequently gain increasing relative importance at a considerable
distance.
CHAPTER XII.
The three rotations TlT1 , •!h TiTa can be most simply found by
putting
GO
8 2 =Mo+l:M~ ..................... (6),
1
00
8,=No+ IN,e~
1
P=!- (9 ~~
oa + Rav)
1D'a{3 ' 1D'
1J = "' V,e•14>
V0 +! ..................... (10).
1
Then we have
Wo = - f v ~~o da -v aiao d/3 •............... (11 ),
334 WANGEBIN'S PBOBLEJI. [205
which is the integral of a complete differential in virtue of the
differential equation for Y,.
Also it can be shewn that
u.·=a+•
aw, , M.}
Ul
o;,
UJ V, = ofJ - vL,
..................... (13),
and
A+21'
UJW,=Z,+zY,+ l--(r+zt)X,
I'
(:\ + 2,u.) ~!
0~
- 2,u. ~p = 0 }
0. ................ (16).
(:\ + 2,u.) of:J + 2p. oa = 0
These a.re found from (3'7) of ch. VII. by taking h.= 1, h1 =h.,
and remembering that fiT ( = 11T1) is the only one of the components
of rotation that occurs.
It is clear from the above equations that(:\+ 2,u.) ~and 2,u.fiT
are conjugate functions of a and {3 and therefore also of a; and y.
We have next to find the displacements u and v from the
equations
so that
336 PLANE STRAIN, [207
and a particular value of If> is the potential of a plane distribution
of density - f:l./27r, and likewise of+ for a distribution - vf'Tr.
This completes the solution in the general case, it will be seen
to be arbitrary in two ways viz. (X + 2p.) 1:1 + 2p.v£, and the com-
plementary (v + £U )/h are any functions of a+ tf:J.
The above includes as a. particular case the theory of solutions
in rectangular coordinates a: and y. In particular problems it is
generally better to use conjugate functions « and f:J, if it can be
arranged that curves a= const. and f:J = const. shall represent the
whole of the boundary.
1:1 = x1 2J£:£[e-(.A.ncosnf:J+Bnainnf:J)+e--(.A.,/cosnf:J+B,.'sinnf:J)]}
4'~=!
'I' p.
:£ [
1
4 (n + 1)
e(n+tl•(-B cosnt:J+.A. sinnt:J)
" ,.... " ,....
- 4 (n_!-1)
____ e-(n-tl•(B" 'cos nt:J
,.....
-.A. " 'sin ..,t:J)J ·
,.,..... '
208] POLAR COORDINATES. 337
1
+ 4 (n _ 1 )(.;:;:1' -~)e-(n-Sl•(.An'cosn#+Bn'sinn#)J. .. (22),
and
ve&== .,.. [ 1 (n+2 n )
~ 4(n+ 1) ~-A+ 2p. e<n+•l•(- B,cosn#+.A,sm n#)
.
l
Thus the particular integrals for u and v, so far as they depend
upon these terms, are
v= ( i A: 1
2f' + 2: a.) (B cos fJ- A sin /3)
...... (25).
209] POLAR COORDINATES. 339
209. General Formula.
Now taking ordinary polar coordinates r and 0 so that e- = r,
{:J = 0, we have for the general forms of u and "
U-
- (1 1}..+
+ log
"2"
r log r) ( 0 + B 0)
2p. + 21'
A
.L11
I
COS 1
I •
Sill
+~
"' ( n+ 21' -p.n) (n+
rn+1 .
2 }..+4 1 )(Ancosn0+Bnsmn8)
"' ( n-2
+:; A.+ 1'- P,
n) (nr-n+_ ) (An' cos nO+ Bn sin nO)
1
1
2 4 1
ao
+I [rn- 1 (-DncosnO+ OnsinnO) +r(n+li(Dn' cos n0-0n Sinn0)] 1
1
............ (26),
and
1
V= ( " 2 " - -
log
}..+ 2p.
r + 1-+2p.
-log-r) (B 1
I
COS
0 -.a.l
A I •
Sln
0)
+7(n--,;:-
ao+2
- :>.. +n21') 4 (nrn+~+ 1) (- Bn cos nO+ An sin nO)
"' (n - 2 n ) (n-l)(Bncosn0-An
+; -----;--}..+ 1'
r-n+ 1
smnO)
1 1 •
2 4
"' [rn-1 (On cos nO+ Dn sin nO)+ r(n+1) (0,: cos nO+ Dn1 sin nO)]
+I
1
............ (27).
In the same notation we have
~ = }..+
_!._ i[rn(Ancosn0+Bnsinn0)+rn(An COBn0+Bn BinnO)]}
2I' 1
1 1
2'111' = !
I'
i
1
[r" (- Bn COB nO+ An sin nO) + r-1l ( Bn' cos nO -An' sin nO)]
............ (28);
and the tractions at a cylindrical surface r = const. are given by
F=~+2p.~
ov J....................(29).
l
G =- 2p.'ID' + 2p. ar
This gives means for the complete analytical solution of any
problem of plane strain in a. solid bounded by coaxal circular
cylinders.
210. Particular Example.
As an example we may consider the case where there is a
cylindrical cavity of radius a in an infinite solid, and at an infinite
22-2·
340 PLANE STRAIN. [210
distance there is a distribution of shear. To represent this we may
take the displacements referred to the system of fixed axes of a;
and y to be U imd V, and suppose that at an infinite distance
U=sy, V=O.
In the notation of the last article the conditions at infinity
become
u cos 8 -vsin 8= sr sin 8,
usin8+vcos8=0;
or u = tsr sin 28,
v = !sr cos 28 - !sr.
We have already seen that u = 0, v =-tar satisfy the equa-
tions, and make the tractions F and G vanish a.t every cylindrical
surface, so that we shall have to add this solution to the solution for
u = tar sin 28, v = tar cos 28,
when r is very great.
From the general solutions (26) and (2'7) we have to keep the
terms in B~', 02 , 0~'· To satisfy the conditions when r is very
great we have to take 0~ = !s. The condition that there is no
traction across the surface r =a gives two relations among the
three constants by which B.; and 0,' are determined. The work
may be left to the reader, and the result is that .
u= ( X+
-- 2}-1- -+
}..+J.' r r'
a')
a~ ! r - t - ssin 28
'
}
...... (30).
v = (____!!:__ as + !r + t a') s cos 28 -.far
X+l' r r'
It may be as well perhaps to remark that this problem does
not, like the corresponding one in art. 169, yield a result in
connexion with the theory of torsion. In the case of torsion a
very important part of the shear consists as we know of a. shifl;ing
of the fibres of the twisted prism parallel to the axis of the prism,
and our work above, being confined to displacements in one plane,
does not take this into account.
211. Elliptic Coordinates.
We shall next consider the case of elliptic coordinates given by
a;+ £'!1 = c cosh (a+ t/3) ................... (31),
and suppose in the first place that the elastic medium extends to
infinity, and is bounded internally by an elliptic cylindrical surface
211] ELLIPI'IC COORDINATES. 341
of the family «, say ex= exo, which is deformed in a given manner.
Then, according to art. 206, we have to take for~ and 2v series of
the form
a= !Ane--e'nll,
2w£ =X+ 2J.Io IAne--e'nll,
p.
in which An is a complex constant ; and we may at the end keep
only the real part of the solution.
Now the displacements have to be found from the equations
a a (u) a (v)
hs = aex h. + a{3 h '
2w a (v) a
-h, = aa li - o/3 7i •
(u)
where h-» = -lei (cosh 2ex - cos 2/3).
The functions M-t, 2vh-t can each be expressed as sums of
terms of the forms
e-(n+l)•} c?s n/3 and {e--} c?s (n + 2) /3
{ smn/3' sm(n+2)/3'
and the equations for u and v can clearly be satisfied by assuming
for u and v sums of terms of these forms with suitable coefficients
These are the particular integrals of the equations for the dis-
placements, and the complementary functions will be found by
taking (v + tu)/h any function of (ex+ t/3) and therefore by taking
for vfh and ufh functions of the same forms as a and 2p.v/(X + 2p.).
Now suppose definitely that
~ = & (X! 2p.) Ie-- (.An cos n/3 + Bn sin n/3) l ...... (32);
;p. Ie-- (Bn cos n/3- An sin n/3)
2v =
J
then we can easily verify that
(X+ 2p.)~ =! [e-(n+ll•{(.An -An+s)cosnfJ+(Bn -Bn+s)sin n/3}
+ e-!n-t)a {(.An- An-t) cosn/3 + (Bn- Bn-t) sin n/3}]
l
2p. ~=I [e-!n+tl•{(Bn- Bn+t) cosn/3- (An- An+s) sin n/3} J
+ e- (n-t) • {(Bn- Bn-t) cos n/3- (.An -An-t) sin n/3} J
......... (33);
342 PLANE STRAIN. [211
and again we can easily verify that ufh and vfh are given by
u e-ltHI)a (n n+ 2) { }
Ji= I 4(n + 1 ) p,- }..+ 2p. (.An-.A.nH)cosn,8+(Bn-Bn+t)sinn,8
-- - - (n+2
e--Mo - - - -n-) (.An - .A.n+t)
4(n+ 1) p. }.. + 2p.
l
211] ELLIPTIC COORDINATES. 343
From the first two we get a difference-equation for the A's,
and from the second two we get a difference-equation for the
B's. When these are solved Dn and On are given by one of (37)
and one of (38).
As an example 1 suppose the cylinder «o, whose principal semi-
axes are a and b, turned through a small angle cf>. The boundary-
conditions are
h
h = e-.• (a+2 b)' cf> sm
u . 'II
2fJ, li = abcf>, w en a= «o.
All the A's vanish, Bo and all the odd B's vanish, and B, = B, = ... ,
all the D's vanish and all the CJs except 00 and 0,, and we find
p.
0, = - t (a + b)' A + 3,u. cf>,
B, = 20r"'" (A+ 2,u.),
and
v
h = abcf>
+Ha+ b)' cf> (e-'Jtl.o -e-.) (A+ 2p.)~3:,u.cos 2fJ.
THE term " shear " was first used by engineers to denote tangential stress,
and is so used in Rankine's Applied Mechanic&. The usage of it for sliding
strain in this work might be justified by reference to Sir W. Thomson, now
Lord Kelvin, and many other eminent authorities, theoretical and practical.
The kind of strain called shear has been considered in ch. I, and the kind of
stress called shearing stress has been considered in ch. IL The object of
this note is to insist more fully than is done in those chapters on the twofold
character of both shear and shearing stress as they occur in the mathematical
expressions. For simplicity we sha.ll limit our consideration to the case of
infinitesimal displa.cements.
The shears are represented by such expressions as 'Ow(Oy + ov(iJz. Now
this expression is the sum of two simple shears, viz. : a. simple shear ?no(Oy of
the planes y=const. pa.ra.llel to the axis z, and a. simple shear ov(iJz of the
planes z=const. parallel to the axis y. In like manner if we define the
(infinitesimal) shear of two initia.lly rectangular lines (1) and (2) to be the
cosine of the angle between them a.fter strain-a. definition which has been
shewn to coincide with the definition in terms of sliding motion-then this
shear will be made up of a. simple shear pa.ra.llel to (2) of the planes perpen-
dicular to (1), and a. simple shear pa.ra.llel to (1) of the planes perpendicular
to (2). The shears that occur in mathematical expreesions are in fact
generally the sums of two such simple shears which are not a.t first separated.
Thus in the energy-function the terms in a for example are just the same
whatever be the proportion in a of the simple shear pa.ra.llel to !I to that
pa.ra.llel to z.
Shearing stress also is of a. twofold character, but the like ambiguity does
not occur. Shearing stress consists of tangential stresses across two perpen-
dicular planes, but these are always equal.
We know that a. simple shear c is equivalent to equal extension and
contraction each ic, and conversely that equal extension and contraction
each e are equivalent to a. simple shear of amount 2e, and in the same way
the extension and contraction might be taken to be equivalent to two simple
shears each of amount e, which combine in the manner explained above ;
346 NOTES.
or again the same extension and contraction will be the equivalents of two
simple shears whose sum is 2e and whose ratio is anything whatever.
Equal pressure and tension each P are in like manner equivalent to a
shearing stress, but the amount of the shearing stress is P. This shearing
stress is really a stress-system consisting of equal tangential stresses P on two
perpendicular planes.
The above remarks appear to contain the secret of the "discrepant
reckonings of shear and shearing stress" to which Lord Kelvin has frequently
called attention. (See e.g. Thomson and Tait's Nat. Phil. Part II. art. 681,
and Lecture8 on Molecular Dynamic~ p. 176.) The discrepancy appears
to arise from the combination in a shear of two simple shears whose ratio
it is unnecessary to know, while the tangential stresses combined in a
shearing stress are always equal. Writing the discrepant statements in
parallel columns we have
Equal extension and contraction Equal pressure and tension each
each e are equivalent to two simple Pare equivalent to tangential stresses
shears of perpendicular planes ; the on two perpendicular planes ; each of
sum of the shears is 2e and their these is of amount P.
ratio may be anything whatever.
Finally we may note that the values of the two simple shears will be equal
if the strain be pure. It follows that, if we regard any small strain as
analysed into a small rotation and a small pure strain, then the extensions
and contractions to which the pure shears are equivalent are always obtained
from the simple shears by precisely the same rule as that by which the
pressures and tensions are obtained from the tangential stresses.
about a given point is extended, and the surface-tractions that must be applied
to hold the surface fixed when the same state of dilatation is produoed in the
spherica.l element.
Equation (40} on p. 244 shews that the dilatation produced at any point by
a given system of surface-displacements is proportional to the work done by
the tractions that must be applied to hold the Bllrface fixed, when there is
dilatation of the spherica.l element about the point, acting through the given
surface-displacements; and equation (41} on the same page shews that the
dilatation produced at any point by a given system of rmrface..tractions is
proportional to the work done by these tractions aoting through the displace-
ments that take place when the surface is free and there is dilatation of the
spherica.l element about the point.
There is a like interpretation of such equations as (48) and (46} on p. 246
for rotation about any given line in terms of the tractions that must be applied
to hold the surface fixed when a spherica.l element about a given point is made
to rotate about the line, and of the displacements that take place when the
surface is free and a similar rotation is effected at the point. In fact in the
above statements we have merely to read 'rotation about a given line' for
'dilatation'.
INDEX.
END OF VOL. I.