Elastic Media With Microstructure
Elastic Media With Microstructure
With 20 Figures
Springer-Verlag
Berlin Heidelberg New York Tokyo 1983
Professor Isaak A. Kunin
Department of Mechanical Engineering, University of Houston
Houston, TX 77004, USA
Guest Editor:
Professor Dr. Ekkehart Kroner
Institut flir Physik, Universitat Stuttgart
D-7000 Stuttgart 80, Fed. Rep. of Germany
Series Editors:
Professor Dr. Manuel Cardona
Professor Dr. Peter Fulde
Professor Dr. Hans-Joachim Queisser
Max-Planck-Institut flir Festk6rperforschung, Heisenbergstrasse 1
D-7000 Stuttgart 80, Fed. Rep. of Germany
ISBN-13:978-3-642-81962-9 e-ISBN-13:978-3-642-81960-5
DOl: 10.1007/978-3-642-81960-5
Library of Congress Cataloging in Publication Data. (Revised for volume 2). Kunin, I. A. (Isaak
Abramovich), 1924- Elastic media with microstructure. (Springer series in solid-state sciences; ) Rev. and
updated translation of: Teoriia uprugikh sred s mikrostrukturoI. Bibliography: v. 1, p. Includes index.
Contents: 1. One-dimensional models. - v. 2. Three-dimensional models. 1. Elasticity. I. Title. II. Series.
QA931.K8913 1982 531'.3823 81-18268
This work is subject to copyright. All rights are reserved, whether the whole or part of the material is
concerned, specifically those of translation, reprinting, reuse of illustrations, broadcasting, reproduction
by photocopying machine or similar means, and storage in data banks. Under § 54 ofthe German Copyright
Law, where copies are made for other than private use, a fee is payable to "Verwertungsgesellschaft Wort",
Munich.
© by Springer-Verlag Berlin Heidelberg 1983
Softcover reprint of the hardcover 1st edition 1983
The use of registered names, trademarks, etc. in this publication does not imply, even in the absence
of a specific statement, that such names are exempt from the relevant protective laws and regulations
and therefore free for general use.
2153/3130-543210
Preface
Crystals and polycrystals, composites and polymers, grids and multibar systems
can be considered as examples of media with microstructure. A characteristic
feature of all such models is the existence of scale parameters which are con-
nected with microgeometry or long-range interacting forces. As a result the cor-
responding theory must essentially be a nonlocal one.
This treatment provides a systematic investigation of the effects of micro-
structure, inner degrees of freedom and non locality in elastic media. The prop-
agation of linear and nonlinear waves in dispersive media, static, deterministic
and stochastic problems, and the theory of local defects and dislocations are
considered in detail. Especial attention is paid to approximate models and lim-
iting transitions to classical elasticity.
The book forms the second part of a revised and updated edition of the
author's monograph published under the same title in Russian in 1975. The
first part (Vol. 26 of Springer Series in Solid-State Sciences) presents a self-
contained theory of one-dimensional models. The theory of three-dimensional
models is considered in this volume.
I would like to thank E. Kroner and A. Seeger for supporting the idea of
an English edition of my original Russian book. I am also grateful to E. Borie,
H. Lotsch and H. Zorski who read the manuscript and offered many sugges-
tions.
1. Introduction ............................................ .
2. Medium of Simple Structure 6
2.1 Quasicontinuum .................................... . 6
2.2 Equations of Motion ................................ . 12
2.3 Elastic Energy Operator ............................... . 13
2.4 Symmetric Stress Tensor and Energy Density ............. . 17
2.5 Homogeneous Media ................................ . 22
2.6 Approximate Models ................................ . 27
2.7 Cubic Lattice ...................................... . 30
2.8 Isotropic Homogeneous Medium, ........................ . 33
2.9 Debye Quasicontinuum .............................. . 37
2.10 Boundary-Value Problems and Surface Waves ............. . 41
2.11 Notes ............................................ . 44
6. Dislocations 144
6.1 Elements of the Continuum Theory of Dislocations ........ . 144
6.2 Some Three-Dimensional Problems .................... . 150
6.3 Two-Dimensional Problems 152
6.4 Screw Dislocations ................................. . 154
6.5 Influence of Change of the Force Constants in Cores of Screw
Dislocations ...................................... . 158
6.6 Edge Dislocations .................................. . 160
6.7 Notes ........................................... . 163
References 249
Bibliography 255
In recent years, new physical and mathematical models of material media, which
can be considered far-reaching generalizations of classical theories of elasticity,
plasticity, and ideal and viscous liquids, have been developed intensively. Such
models have appealed for a number of reasons. Primary among them are the
use of new construction materials in extreme conditions and the intensification
of technological processes. The increasing tendency toward rapprochement of
mechanics with physics is closely connected with these factors. The internal logic
of the development of continuum mechanics as a science is also important.
This treatment is devoted to the study of models of elastic media with mi-
crostructure and to the development of the nonlocal theory of elasticity. Starting
from such models as a crystal lattice and simple discrete mechanical systems,
we develop the theory and its applications in a systematic way.
The Cosserat continuum was historically one of the first models of elastic
media which could not be described within the scope of classical elasticity [1.1]
However, the memoirs of E. and F. Cosserat (1909) remained unnoticeJ for a
long time, and only around 1960 did the generalized models of the Cosserat
continuum start to be developed intensively. They are known as oriented media,
asymmetric, multipolar, micromorphic, couple-stress, etc., theories. For short
we shall call them couple-stress theories. Essential contributions to the develop-
ment of couple-stress theories were made, for example, by Aero, Eringen, Green,
Grioli, Gunther, Herrmann, Koiter, Kuvshinsky, Mindlin, Naghdi, Nowacki,
Palmov, Rivlin, Sternberg, Toupin, and Wozniak; their fundamental works
are listed in the Bibliography. The Bibliography is far from being complete. A
survey of works before 1960 can be found in the fundamental treatment of
Truesdell and Toupin [1.2]; later ones are quoted in papers by Wozniak [1.3],
Savin and Nemish [1.4], Iliushin and Lomakin [1.5], as well as in monographs by
Misicu [1.6] and Nowacki [1.7].
From the very beginning of the development of the generalized Cosserat
models, attention was turned to their connections with the continuum theory of
dislocations. In 1967 a symposium was organized by the International Union of
Theoretical and Applied Mechanics, which had great significance in summing up
the ten-year period of development [1.8]. In the symposium a new trend closely
connected with the theory of the crystal lattice was also presented which con-
tained the above-indicated models as a long-wave approximation, namely, a
nonlocal theory of elasticity. The nonlocal theory of elasticity was also devel-
oped in works of Edelen, Eringen, Green, Kroner, Kunin, Laws and others,
2 1. Introduction
Let us return to nonlocal models. They can be divided into two classes:
discrete and continuous. Discrete structure of a medium could be taken into
account in the usual way, for example, as is done in the theory of the crystal
lattice. However, the apparatus of discrete mathematics is most cumbersome;
therefore, we shall also use the mathematical model of quasicontinuum for an
adequate description of the discrete medium. Its essence is an interpolation of
functions of discrete argument by a special class of analytical functions in such
a way that the condition of one-to-one correspondence between quasicontinuum
and the discrete medium is fulfilled. The advantages of such an approach con-
sist in an ability to describe discrete and continuous media within the scope of a
unified formalism and, in particular, to generalize correctly such concepts of
continuum mechanics as strain and stress. It is to be emphasized that the model
of quasicontinuum is applicable not only to crystal lattices but also to macrosys-
terns.
We shall also distinguish media of simple and complex structures. In the first
case, the displacement vector is the only kinematic variable and it determines a
state of the medium completely. Body forces are the corresponding force vari-
able. To describe a medium of complex structure, a set of microrotations and
microdeformations of different orders characterizing the internal degrees of
freedom and the corresponding force micromoments is additionally introduced.
The difference between media of simple and complex structures, generally
speaking, is conserved in the approximation of weak nonlocality, but this is
displayed only for high enough frequencies of the order of the natural fre-
quencies of the internal degrees of freedom. At low frequencies, the internal
degrees of freedom can be excluded from the equations of motion so that they
will contribute to the effective characteristics of the medium only. The difference
between the quasicontinuum and the continuous medium completely disappears
in the approximation of weak nonlocality.
In the zeroth long-wave approximation, at not very high frequencies, a
complete identification of different models of a medium with microstructure
takes place: all of them are equivalent to the classical model of elastic con-
tinua which was obtained on the basis of general phenomenological postulates.
Only effective elastic moduli "know" about the structure of the initial mic-
romodel, but this information cannot, of course, be derived from them. It fol-
lows that an explicit consideration of miclOstructure effects and, in particular,
of the internal degrees of freedom is possible only with the simultaneous con-
sideration of nonlocality, i.e., a consistent theory of elastic medium with
microstructure must necessarily be nonlocal.
Schematically, the connections between different theories are shown in
Fig. 1.1.
Our main purpose is the investigation of the effects of microstructure and
nonlocality. In addition, we wish to elucidate the domain of applicability of
different theories of media with microstructure. Such theories are considerably
more complex than the usual theory of elasticity, although they reduce to this in
4 1. Introduction
Nonlocal theory of
Fig. 1.1. Illustrating
elasliclly the interconnections
of the various theories
Weakly nonlocal
certain limits; their application is, as a rule, reasonable only when they describe
qualitatively new effects which are not derivable from the local theory.
We find it advisable to consider a number of simple models of media with
microstructure in order to acquire some nonlocal intuition before proceeding to
generalizations.
For these reasons, the treatment is divided into two parts. One-dimensional
models, for which cumbersome tensor algebra is not needed, are studied in the
first volume [1.12]. At the same time, one can trace a number of distinctions of
the nonlocal theory from the local one already in one-dimensional models.
These distinctions have both physical and methodological character. Particu-
larly, one has to analyze critically the possibility of using in nonlocal theories
such habitual notions as stress, strain, energy, density and flux.
In this volume, three-dimensional models of media with microstructure
are considered, the main attention being payed to specific three-dimensional
effects. The second and third chapters are devoted to the general theory of
media with simple and complex structures. Local defects in media with mi-
crostructure as well as in elasticity are considered in Chapter 4. Internal stress
and dislocations in the nonlocal elasticity are considered in Chapters 5 and 6.
Chapter 7 (written by S.K. Kanaun) is devoted to random fields of inhomo-
geneities.
The author has tried to avoid complex mathematical methods in the first
part. It is assumed that the reader is acquainted with the Fourier transform and
has some skill in working with a-functions (though it is easy to acquire it in
the process of reading the book).
The necessity of simultaneous use of space-time and frequency-wave re-
presentations is one of the peculiarities of the formalism used in the nonlocal
theory. This way of thinking is quite habitual for physicists, but could be, to
some extent, unusual for engineers. Because of this, the rate of presenting the
material is slow in the beginning and then speeds up gradually. For the same
1. Introduction 5
purpose, a number of results are presented in the form of problems, which are
considered a part of the text and they are referred to in the later text.
In conclusion, let us note that a number of important questions are omitted.
In particular, this is related to the thermodynamics of the nonlocal models,
which was contained in the original plan of the book. Unfortunately, the author
has not succeeded in representing this problem in a sufficiently simple and phy-
sically motivated form. An axiomatic approach to the nonIocal thermody-
namics was developed by Edelen [1.13].
2. Medium of Simple Structure
2.1 Quasicontinuum
IHere and in the following material by identical superscripts and subscripts we contemplate
summation.
2.1 QuasicontinuuID 7
(2.1.2)
Then k = kf3ef3 and kx = kaxa. Note that the last expression is also under-
stood as an ordinary scalar product (k, x), if between the spaces E3 (with fixed
origin) and E~ an identification is established with the help of the metric tensor
gafJ' i.e., the operation of raising and lowering of indices. For our purpose such
an identification is not convenient and we shall consider separately the "physi-
cal" x-space E3 and the dual k-space E~.
Let us construct a parallelepiped B { - 11: :s; kfJ :s; 11:}, in the k-space, whose
edges are parallel to vectors of the reciprocal basis ef3 and set
o(x) ~ _1_ Je
B (211:)3 B
ixk dk = _1~ IT
11:3Vo fJ=l
sin(xef3 )
xe f3 '
(2.1.3)
Defining u(k) to be zero outside B, and taking into account (2.1.3) we find
for the inverse Fourier transform u(x)
u(n) ~ u(x) ~ u(k) is well guaranteed by the condition of truncating the Fou-
rier-transform u(k) of the function u(x) and the uniqueness of the expansion
(2.1.4).
Let us proceed to the general case of the space N'(B), whose elements are
functions u(n), which increase no faster than some power of Inl as Inl --+ 00.
The space N'(B) can also be regarded as a space of generalized functions on the
space N(8), the latter being the space of rapidly decreasing functions.
Recall that in the one-dimensional case [1.12], when interpolating functions
u(n) E N'(B), it turned out to be necessary to identify the ends of the segment B,
transforming it into a circle. Analogously, let us identify the opposite faces of
the present parallelepiped B or, in other words, let us transform it into a three-
dimensional torus.
The series (2.1.4) now converges only in the sense of the generalized func-
tiom [2.1], and u(k) belongs to the space of the generalized functions K'(B),
defined on the space K(B) of infinitely differentiable test functions with sup-
ports in B.
Let us denote the inverse Fourier transforms of the spaces K(B) and K'(B)
by X(B) and X'(8), respectively. The inverse Fourier transform u(x) E X'(B) of a
function u(k) E K'(8) is defined, as usual, with the help of the Parseval equality
(2.1.6)
where q>(x) E X(B) and q>(k) E K(B) are connected with each other by the usual
Fourier transform.
It can be proved [2.1] that u(x) E X'(B) are regular generalized functions (i.e.
conesponding to ordinary functions), which increase no faster than a power
of Ixl as Ixl --+ 00 and which can be analytically continued to a complex region
as entire functions of the exponential type :5; Te. The proof of the fact that the
function u(x) has values u(n) at the knots, can be performed analogously to the
one-dimensional case [l.12].
The uniqueness of interpolation is ensured by the identification of the
opposite faces of the parallelepiped B. This excludes from X'(B) functions of
the type P(x) sin[Te(xe)~], [P(x) being a polynomial], which vanish at all the
knots.
Thus we have:
Theorem. The formulae (2.1.4, 6) fix linear isomorphisms between the spaces
We may now consider u(n), u(x) and u(k) as representations of one and the
same function in different functional bases. Of course, to this end it is also
2.1 Quasicontinuum 9
q(k) = 1: q(n)e-ink,
n
kEB,
q(x) = 1: q(n)oB(x -
n
n) . (2.1.9)
Note that, as distinct from (2.1.4, 5), in the correspondence q(n) +-+ q(x) +-+
q(x) the factor Vo is absent. This is connected with the fact that in the further
consideration the quantities q(n) will have the meaning of forces acting at the
knots of a lattice. Thus q(x) can be interpreted as the corresponding density of
body forces, q(k) as its Fourier transform and the functional (qlu) as the work,
done by the forces q on the displacement u.
As an important example, consider the functional
Problem 2.1.2. Show that the following identities are generated by the func-
tional with the kernel I (x) :
Here
= _1 _ SS "'(x
G1tP ~,
x')e-iCkx-k'x'ldx dx',
, (2.1.14)
One can establish the correspondences between the discrete and continuous
operations of multiplication, convolution, differentiation and so on, distinc-
tions compared with the one-dimensional case being small.
Thus, the functions u(n) E N'(B) constitute a ring with respect to the opera-
tion of multiplication. At the same time, the ordinary product of two functions
from X'(B), generally speaking, does not belong to X'(B). Nevertheless, it is
possible to introduce in X'(B) such an operation of multiplication, which cor-
responds to the multiplication of the inverse transforms in N'(B) and whose
result is still contained in X'(B).
In fact, on the torus B we can define an operation of displacement and,
hence, an integral convolution
The model for a medium with simple structure will be based on Born's model
of a simple lattice with the basis ea in the harmonic approximation. For the
elastic energy rp we have, analogously to the case of a one-dimensional chain
[1.12],
where rpaf3 (n, n') is a given tensor of the force constants, and uaCn) is the displace-
ment of the particles from the equilibrium position, i.e. it is assumed that initial
external forces are absent.
This model may be interpreted, if necessary, as a system of pointwise parti-
cles situated at the knots of the lattice and interacting by means of linear elastic
bonds of the most general nature.
The Lagrangian of this system has the form
L = i ~ gafJ m(n)uaCn)uf3(n)
where men) is the mass of a particle at the point n, and qa(n) is the external
force acting on the particle. The dependence of the field parameters on time is
not indicated explicitly.
In the usual way we find the equations of motion
where rpafJ(x, x') and rpaf3 (k, k') are expressed in terms of the force constants
(2.2.7)
(2.2.10)
The properties of the operator <pa~(x, x') are in many aspects analogous to the
properties of the operator <P(x, x'), which were considered in details in [Ref.
1.12, Sect. 2.4]. We shall therefore concentrate on features specific to the
three-dimensional case.
Hermiticity. From (2.2.4), it follows that
From this follows the possibility of representing cpa~(k, k') in the form
(2.3.5)
where c).aJ1.~(k, k') is an analytic function of the same type as cpa~(k, k'); without
loss of generality it can be considered to be symmetric in the indices inside the
first and second pairs and to be Hermitian with respect to permutation of the
pairs, i.e.
(2.3.6)
Local approximation: If the field parameters vary slowly over distances of the
order of I, then
2Here and subsequently indices which are contained in brackets (parentheses) are assumed
to be subject to the operation of alternation (symmetrization).
2.3 Elastic Energy Operator 15
Media with binary interaction: One of the possible specific properties is the
following one. The forces appearing due to a change of the distance between
points x and X', are proportional to this change and are directed along the
line, which connects these two points. It is easy to show that such a medium can
be described by an operator of elastic bonds 1Jfup with the kernel
where 1Jf(x, x') = 1Jf(x', x) is a scalar characteristic of the elastic bond which
connects the points x and x', such that 1Jf(x, x') = 0 when Ix - x'i > L.
For the energy operator rp we may write
The conditions of invariance with respect to rotation are here satisfied au-
tomatically.
Homogeneous media: In this case, which is important for applications
(2.3.15)
(2.3.18)
(2.3.19)
The static Green's tensor Ga{l is defined as an operator which is inverse with
respect to the energy operator (/jafJ and hence
Galx, x') satisfying the corresponding boundary conditions as x, x' --+ 00.
For a homogeneous medium,
(2.3.21 )
i.e. in the k representation, the kernel of the Green's tensor can be obtained by
the purely algebraic operation of inverting the matrix (/jafJ(k).
It is easy to show that from the definition of the Green's tensor, the elastic
energy can be represented in the static case, as
(2.3.22)
where the qa's are body forces. Conversely, this relation can be taken as a de-
finition of the static Green's tensor.
It is obvious that the kernels Ga{l(x, x'), Ga{l(k, k') possess the symmetry
properties (2.3.1, 2). For the quasicontinuum, due to the condition k E B, the
kernel Galx, x') is an entire function of the corresponding arguments.
The dynamic Green's tensor GafJ is defined as an operator inverse to the
operator
We leave it to the reader to write the analogs of the equations (2.3.20) for the
kernels GafJ(x, x', t), GafJ(x, x', ev) and Galk, k', ev).
2.4 Symmetric Stress Tensor and Energy Density 17
where the analytic function AJ.(k) is uniquely determed by the condition that
coefficients of its Taylor series
00
(2.4.3)
(2.4.4)
holds.
Lemma 3. If an analytic function Aa(k) satisfies the conditions
(2.4.5)
(2.4.6)
where the analytic function BJ.a(k) is symmetric with respect to the indices A, a.
In order to prove this let us write the function Aa(k) in the form
(2.4.7)
18 2. Medium of Simple Structure
The constant tensor A&a appearing here is symmetric because of the second
of the conditions (2.4.5); the tensor function A}.a"(k) is symmetric with respect
to the indices Afl- and, according to Lemma 1, is defined uniquely. This function
can be replaced by a tensor function A{Aa")(k) from Lemma 2, which is sym-
metric with respect to the indices A, a, since the factor k}.k" provides the opera-
tion of symmetrization needed in the identity (2.4.4). As a result, we find
(2.4.8)
(2.4.10)
The requirement concerning invariance of the energy (/J with respect to the
transformation Ua -> Ua + u£, (2.2.4) and the Hermiticity of the operator
(/Ja{3 lead to the following equation
which is to be satisfied for any u£ and u{3' From this one can obtain conditions
which must be satisfied by the kernel of the operator (/Ja{3'
First, let us make one remark. In the case of an infinite medium, the ex-
pression for the energy (2.2.4) and hence the second term in (2.4.11) are not a
priori defined on displacements ug of the form (2.4.9) or (2.4.10). Consequently
let us consider next a bounded medium and then perform the transition to the
limit.
Taking into account the arbitrariness of u{3' we conclude that (2.4.11) is
equivalent to the condition
Ju£(x)(/Ja{3(x, x')dx = 0 (2.4.12)
2.4 Symmetric Stress Tensor and Energy Density 19
or in the k representation
Using (2.4.9), aa being arbitrary, gives two equivalent conditions for f/Ja~
Due to the assumption about the boundedness of the medium the function
x') has finite support and hence the function f/Jap(k, k') is analytic with
f/Ja~(x,
respect to the pair of points k, k'. Taking into account also the conditions of
Hermiticity (2.3.1) we find that the relations (2.4.14) are equivalent to the
possibility of representing f/Ja~(k, k') in the form
we have
(2.4.18)
k ).f/JAa[p.(31(k, 0) = O. (2.4.19)
We now apply Lemma 3 to the function f/Ja~(k, k') for each argument k and
k'. This leads to the required result:
(2.4.20)
Here, the tensor cAap.~(k, k') satisfies the symmetry conditions (2.3.6) and
is connected with f/JAap.~(k, k') by
(2.4.21)
20 2. Medium of Simple Structure
Problem 2.4.1. Check that c Aap{3(k, k') satisfies the conditions (2.3.6).
It was assumed above that the medium is finite. However, the derived
result about the representability of l/Ja{3(k, k') in the form of (2.3.5) with the
tensor c Aap{3(k, k') possessing the symmetry (2.3.6) depends neither upon the
shape nor upon the size of the elastic body. Hence this remains valid also
for the limiting case of an infinite medium. According to the same reasoning,
if the second term in (2.4.11), vanishes for an abritrary finite medium, it can
be considered to be equal to zero also in the limiting case of an infinite medium
(natural regularization).
Let us now introduce the strain tensor
(2.4.22)
(2.4.23)
or in the k-representation
(2.4.25)
(2.4.27)
or
From the above expressions it follows that O'Aa(x) and rp(x) may be inter-
preted as the stress tensor and the elastic energy density, respectively, both
quantities being invariant with respect to translation and rotation, as they
should be according to their physical meaning.
2.4 Symmetric Stress Tensor and Energy Density 21
(1 = Co:, (2.4.29)
Problem 2.5.1. Show that a necessary condition for this invariance is the re-
presentability of (/JafJ(k) in the form
(2.5.2)
[Hint: Expand (/JafJ(k) in a series, write (2.4.12) and (2.4.16) in the k-represen-
tation and take into account the condition of Hermiticity.]
(2.5.3)
The energy density 'Po(x) differs from q>Q(x) by divergence terms, which must
2.5 Homogeneous Media 23
Problem 2.5.2. Verify that the most general expression for the divergence
part, which is quadratic with respect to the displacement and invariant with
respect to translation, has the form
(2.5.5)
(2.5.6)
where ~}J~ and b).al'~ possess the above-described symmetries. The conditions
of invariance with respect to rotation yield
(2.5.7)
Regarding (2.5.7) as an equation for b).a}J~, let us find conditions for its
solvability. For this purpose it is necessary to employ a simple but somewhat
cumbersome procedure: to symmetrize (2.5.7) with respect to Af-t and a{3
separately, then to take the difference of these expressions and, finally, to
symmetrize the difference with respect to Aa. The final result will contain
only two components of the tensor (/loa}J~ with permuted pairs of indices,
~}J~ = (/Yij).~}J. (2.5.8)
This is, in fact, the desired additional condition which ensures the invariance
of <po(x) with respect to rotation.
Problem 2.5.4. Show that if the condition (2.5.8) is satisfied, (2.5.7) has only
one solution:
b).a}J~ = (/J{f}J~ - rpArf"~, (2.5.9)
which possesses the necessary symmetry.
24 2. Medium of Simple Structure
(2.5.10)
where
(2.5.11)
(2.5.12)
(2.5.13)
(2.5.14)
(2.5.15)
These expressions enable us to express explicitly C$a p(3 in terms of the force
constants and to interpret (2.5.8) as restrictions on the force constants, due
to the requirements of invariance of the energy with respect to rotation.
Thus, we have completed the first part of the problem. Now it is necessary
to show that (2.5.2) admits a representation in the form
(2.5.16)
where the operator of the elastic moduli c),ap(3(k) possesses the symmetry (2.3.19),
and to find a connection between (/JJ.a p(3(k) and cJ.ap(3(k).
Let us now write down (2.2.4) for the elastic energy (/J in the form
(2.5.17)
2.5 Homogeneous Media 25
(2.5.18)
which is symmetric with respect to the indices ila because of the symmetry
properties of (/)).apfJ(k) (Lemma 2). This enables us to write (2.5.17) in the form
(2.5.19)
(2.5.20)
(2.5.21)
or
(2.5.23)
(2.5.24)
where
(2.5.25)
r.'apfJ _
U2 - ~
1 (a).a(pfJ)
1
+ a(p{3),a)
1, (2.5.26)
(2.5.27)
26 2. Medium of Simple Structure
<c).ala~a,ut3/c,ut3> ,
where
(2.5.28)
(2.5.29)
(2.5.30)
satisfies the symmetry conditions (2.3.19) and is connected with <pAa,utS(k) by the
above-described expressions, which implicitly determine the operator S.
Hooke's law in the operator form (2.4.23, 24) takes now the form
(2.5.32)
Problem 2.5.5. Show that the real (imaginary) part of c Aa ,uf3(k) is an even (odd)
function of k and hence, for an isotropic medium and for a medium with central
symmetry, c Aa,uf3(k) is a real function.
The equation of motion (2.2.8) for a homogeneous medium has the form
(2.5.33)
(2.5.34)
Taking into acount (2.5.16), the equations of motion may be also written
in the following form
J
pua(x) - 0). cAa,utS(x - x')0,uuix') dx' = qa(x) , (2.5.35)
Wi = w,{k), i = 1, 2, 3 (2.5.38)
give the dependence of the frequencies on the wave vector for the corresponding
mode of free vibrations. At the same time, according to (2.5.2) all the modes
pass through the points w = 0, k = O.
To real values of wand k, correspond nondecaying waves of the type
Generally speaking, the group velocity depends upon the modulus of the wave
vector k, i.e. spatial dispersion takes place. The absence of nonzero real roots
of the equation wi(k) = 0 is a necessary condition of stability of the medium.
In fact, such roots would correspond to static displacements different from
translation and rotation, which do not contribute to the elastic energy; this is
not possible for a stable elastic medium.
However, as in the one-dimensional case, it is possible to have complex
roots of the function w,(k), which are caused by the action-at-a-distance. To
these correspond exponentially increasing (or decreasing) displacements ua(x),
which cannot be considered to be admissible functions for an infinite medium.
At the same time such solutions of the static equation should be taken into ac-
count in boundary-value problems, where they can make a considerable con-
tribution to the elastic energy by means of the forces which act in the boundary
region.
As a rule, for high frequencies only complex values of the wave vector can
satisfy (2.5.37). Solutions having such frequencies must decay, while propagat-
ing, but this is not connected with dissipation of energy, since for such waves
the energy flux is equal to zero.
where the cn's are real constant tensors, which can be explicitly expressed in
terms of force constants of the micro-model.
For example, in the zeroth approximation, taking into account that
(2.6.3)
which coincides with the expression for the elastic moduli, known in the theory
of a crystal lattice [2.2].
Thus, in the zeroth approximation we obtain the equations of the classical
theory of elasticity
(2.6.4)
with the elastic-constants tensor c8al'fJ. Note that this approximation is universal
in the sense that it is physically correct for an infinite medium, as well as for
boundary-value problems.
As we saw in the examples of one-dimensional models, when constructing
successive approximations, it is usually necessary to distinguish the cases of
infinite and finite media. In this section, we consider infinite media only. As
distinct form the classical theory of elasticity, for nonlocal and weakly nonlocal
theories this case has a fundamental significance.
Recall that in models of media with weak dispersion, which correspond to
the long-wavelength approximation, only the fields u(x) whose Fourier trans-
forms u(k) are localized in the region of small (as compared with [-I) wave num-
bers are strictly speaking, to be considered. This enables one to restrict oneself
to a finite number of terms in (2.6.1), or, equivalently in the expansion of the
functions of bounded support CAal'fJ(x) in a multipole series.
To a first approximation
(2.6.5)
2.6 Approximate Models 29
(2.6.6)
and the equations of motion (2.5.35) are differential equations of third order
with respect to space derivatives.
It is known [2.2] that such models in crystal optics describe the specific
phenomenon of a rotation of the polarization plane of a propogating electro-
magnetic wave. Media for which this effect takes place are called gyrotropic.
An anologous phenomenon is possible also for an elastic medium, when cfPi3< #-
O. Unfortunately, as far as the author knows, no experimental investigations of
rotating of the polarization plane of elastic waves in a gyrotropic medium have
been carried out up to now. At the same time, such an investigation must be of
interest, because gyrotropy is one of not very many qualitatively new phenomena
which arise already in taking into account weak spatial dispersion. This is con-
nected with the fact that, despite the small value of the second term in (2.6.5),
the effect accumulates over a sufficiently large path and hence the plane of
polarization rotates through a finite angle.
For media with central symmetry, as well as for the isotropic medium, we
have
and, as can be easily shown (Problem 2.5.5), Cfpf3< = 0, i.e. such media are not
gyrotropic. For such media the next representation after the zeroth one is the
second approximation
(2.6.8)
C2 = /2C2' (2.6.10)
where C2 has the dimension of ordinary elastic moduli and / is the characteristic
scale parameter. The latter can also be considered as a small parameter with
respect to which the expansion is carried out. The equations of motion (2.5.35)
are now written in the form
(2.6.11)
of the order of I. The connection between the elastic moduli and rpap(x, x')
or cJ.app(x, x') can be found, as in [Ref. 1.12, Sect. 2.4].
Thus the equations for the medium with weak dispersion belong to the
well studied class of equations with a small coefficient of higher derivatives.
It is necessary to take this coefficient into account when constructing solutions
as well as when determining the field of applicability. In particular, such equa-
tions are also connected with phenomena of the boundary-layer type.
In an obvious way, the equations for higher-order approximations can be
written down, but they have no significance, since they do not describe qualita-
tively new phenomena and may provide at most a small correction.
Equation (2.6.11) of a medium of simple structure with weak spatial
dispersion does not differ in essence from the equations of the so-called couple-
stress theory of elasticity with constrained rotation, which was developed
by a number of authors [2.3-5]. A superficial distinction is only in the form
in which the equations of motion are written down, since in the couple-stress
theory of elasticity a nonsymmetric stress tensor usually appears. Such a
distinction is connected with the different definitions of the stress tensor and
has mainly a methodological character: the equations of motion expressed
in terms of displacements coincide.
However, the above-accepted physical interpretation of the approximate
model (2.6.11) differs fundamentally from the usual interpretation of the
couple-stress theory of elasticity. The latter is introduced purely phenomenolo-
gically in analogy with the classical theory of elasticity, and the parameter I
is not regarded as small. 3 As a result, typical problems, which are often solved
in connection with couple-stress theories, concern calculations of stress con-
centration on inhomogeneties having dimensions of the order of or even less
than I. It is natural that values of the stress concentration can differ strongly
from values obtained on the basis of the classical theory of elasticity, and
this creates the illusion of a theory being significantly improved. The preceeding
analysis shows that these results should rather be considered as a brave ex-
trapolation of the approximate model (2.6.11) in a field in which, generally
speaking, only the strongly nonlocal model is applicable.
to the binary interaction (with the bond of the type of a spring) since such
a lattice would be unstable; it would be possible to fold up such a lattice without
changing the distances between neighboring knots. Therefore, for a lattice
with nearest-neighbor interaction, it is necessary to introduce bonds of the
rod type. In the simplest case, each of these rods would be characterized by
two parameters: the longitudinal and the transverse stiffnesses A; and B;,
respectively (i = 1, 2, 3).
Taking into account the above reasoning let us give, in the lattice coordinate
system, the matrix of force constants in the form (see also [2.6]):
°° 0,0, 1, ° 0,0,
0, 0, 1
-1, °
n 1, 0,
I -1, 0, -1
Al ° ° °° °° B2 B3
(2.7.1)
- vol!j)afJ(n)
°° ° ° °° ° ° °° ° °
Bl
Bl
A2
B2
B3
A3
The first condition from (2.3.17) gives (here and subsequently only nonzero
components are presented and there is no summation over identical indices):
(2.7.2)
From (2.3.15) it follows that the matrix !j)afJ(k) is diagonal. For its com-
ponents in the crystallographic system of coordinates we find (a =1= (3 =1= r)
From these expressions, according to (2.5.2) we easily find the tensor !j)}.ap.{3(k),
which is symmetric with respect to the indices )" fJ-
(2.7.4)
For the components (2.5.20) which are symmetric and scewsymmetric with
respect to the second pair of indices we obtain
(2.7.7)
(2.7.8)
(2.7.10)
capap = B[S2(ka) + s2(k p) - 1],
Note that all the computations here were carried out in the lattice coordi-
nate system and therefore the components of the tensor c"app have the dimen-
sion [Nm-6]. In order to obtain the usual dimension, i.e. [Nm-2], it is neces-
sary to consider the mixed components d~p. Since gap '" 12, where 1 is the
length of a rod, equal to the parameter of the cubic lattice, we have
(2.7.11)
where the last equality is valid only in the lattice coordinate system.
For the elastic-constants tensor in the zeroth approximation, we find
where
Note that the different signs of the Lame constants Ao and f1.0 additionally
indicate the artificial nature of the model and, in particular, the fact that the
model does not correspond to a system of rods. In fact, in the latter, a side
thrust should not result from tension along a lattice axis, and this is possible
only if Ao = O. In the next chapter we shall see that in a correct consideration
of rod systems in terms of media of complex structure, such contradictions
do not appear.
Problem 2.7.1. Write down the elastic constants in the second approximation.
Problem 2.7.2. Taking into account (2.7.3), derive an expression for the dis-
persive modes w,{k) (j = 1, 2, 3).
As a start, let us define our notation in the case of isotropic media. The position
vectors of the physical and Fourier spaces will be denoted by rand k, respec-
tively. Their moduli are denoted by rand k, i.e. r = Irl, k = Ikl. Note also
34 2. Medium of Simple Structure
that in the three-dimensional case, a function <p(r) and its Fourier transform
<p(k) are connected by the following relations
(2.8.1)
<p(r) = - 4~2r r
-00
k<p(k)e irk dk , (2.8.2)
where )'(k), fJ,(k) are real scalar functions. When the action-at-a-distance is
finite, ), and fJ, are analytic functions of the vector k, and their expansions with
respect to k can contain only terms with k 2 (as a consequence of isotropy), i.e.
)'(k) and fJ,(k) are even analytic functions of k.
It can be easily seen that (2.8.3) for the operator <J)afJ(k) can be associated
with the corresponding operator of the elastic moduli
(2.8.4)
which has the same structure as the tensor of elastic moduli in the local theory
of elasticity. However, as distinct from the local theory, this representation for
clal'fJ(k) is not unique. In fact, without destroying the symmetry, one can add
to (2.8.4) terms of the form
which would only change the connection between the scalar coefficients of the
tensors <J)afl(k) and clal'fl(k). Thus, we have here one more example of nonuni-
queness of the operator of elastic moduli (and hence the stress tensor) in a non-
local theory.
The equations of motion (2.5.33) in direct notation have the form
which is identical with the equation for an isotropic medium in the local theory,
but with substitution of the Lame constrants ),0, fJ,0 by scalar operators A, M
2.8 Isotropic Homogeneous Medium 3S
with the kernels il(r), p(r). The latter are inverse Fourier transforms of il(k),
p(k). Note that the nonuniqueness of the operator of elastic moduli cAap{J does
not influence the equations of motion whatsoever: these are fully determined by
the elastic-energy operator f/Ja{J.
In order to obtain solutions of the equations of motion it is sufficient to
construct the Green's tensor G{r, t). Formally this can be acheived by replacing
the Lame constants ih the Green's tensor of the local theory of elasticity by
the corresponding operators. For example, for the static Green's tensor we
have
However, this representation is not very suitable for the investigation of specific
nonlocal phenomena. Consequently, we shall consider another approach.
An arbitrary vector function uP(k) can be decomposed into its longitudinal
(I) and transverse (t) components, as a result of the projection operators
(2.8.7)
(2.8.9)
where
(2.8.10)
Applying the operators 1C and {) to the equations of motion we find that these
equations are decomposed into independent equations with scalar operators
(2.8.12)
(2.8.13)
(2.8.14)
where
1 . 1 (2.8.15)
-P-'[(V"7(7ik~)---(V"'2] , J = ,t.
where, as (V -+ 0,
p(V2 2 p(V2
kr,o «(V) -+ ). +2 ,kt,o«(V) -+ - - . (2.8.17)
o flo flo
The inverse functions Gi(k, (V) = (/Jjl(k, (V) can be represented in the form
([Ref. 1.12., Eq. (2.9.17)])
Thus for the isotropic medium, the structure of the solutions is determined
2.9 Debye Quasicontinuum 37
by the distribution of roots of the functions l/Jj(k, m) as was also true in the one-
dimensional case.
When proceeding to models with weak nonlocality, we used the long-wave-
length approximation. However, for isotropic media, in a number of cases, it
may be more suitable to use an approximation with respect to the first roots
(Sect. 2.10). In this approximation, in the static case (m = 0)
(2.8.20)
We shall also need the Green's functions of the Laplacian g(r) and bi-
Laplacian h(r), i.e. the inverse Fourier transforms of - k 2 and k4, which satisfy
the condition
/:::..2h(r) = /:::..g(r) = air) . (2.9.2)
Simple computations yield
38 2. Medium of Simple Structure
g(r) = - - 12
2 Si(Kr), (2.9.3)
'ltr
h( r ) = - - 1
42
( .
KrSl(Kr) sin,.,
+- -+ cosKr ) , (2.9.4)
'It K Kr
(2.9.5)
where Vb VI are the logitudinal and transverse velocities of sound in the cor-
responding local model of the elastic continuum. It is evident that this model
does not take into account the spatial dispersion.
In a more realistic Born-Karman model, a dispersion law is accepted, which
is analogous to the case of a chain with interaction between nearest neighbors
2KV j 'ltk
(J)' (k)
.
J
= -~
'It
SIn -
2K
. (2.9.6)
(2.9.7)
(2.9.8)
(2.9.9)
is shown in Fig. 2.1. The straight line corresponds to the Debye model where
w~ = vi" is the Debye frequency. The dotted curve corresponds to the Bom-
Karman model and, as is clear from the picture, practically it coincides with the
curve for rj = 0.5.
1.0 ~
wJwf· ,//
//
I} 0
/1/
/'"
,--
If y --r
/
f/j/
j
1/..-- y(1O
I
If Fig. 2.1. Dispersion curves
o k/x !.O
The parameter n may be connected with the ratio of the boundary frequency
Wj(") to the Debye frequency w~, i.e.,
(2.9.10)
(2.9.11)
(2.9.12)
(2.9.13)
40 2. Medium of Simple Structure
where
(2.9.15)
(2.9.16)
where Qa = const, and a(r) can also be replaced by air) since due to the trunca-
tion of the Fourier spectrum of the Green's tensor, the final result is the same.
This reasoning enables us, in what follows, to make no distinction between the
"true" a-function and the a-function with the truncated Fourier spectrum,
when we consider the quasicontinuum model.
For the energy of the elastic field of a concentrated force, taking into account
(2.3.16), we have
(2.9.17)
where the (Xi are dimensionless constants, which depend on the form of the
functions oJj{k),
(2.9.18)
For the Debye model and for the model (2.9.9) we find
hj{r) ~ - 8r
10
(1 - 22Ti2 ).
ICr
(2.9.19)
Substitution in (2.9.12) and (2.9.13) yields the dominant terms of the asymptotic
form of Gap(r) to be of the order of ,-1 and ,-3, which, as it should, coincides
with the asymptotic form of the Green's tensor for the elastic continuum with
action-at-a-distance examined in the preceeding section. As IC ~ 00, the Green's
tensor is transformed into the Green's tensor of the local theory of elasticity
and (2.9.17) diverges linearly with IC.
The results obtained in [Ref. 1.12, Sects. 2.6, 11] can be transferred to the three-
dimensional case but here some specific three-dimensional effects arise; in par-
ticular, special types of surface waves appear.
Let us restrict ourselves to a consideration of the simplest type of model of
a medium with binary interaction. In this case the operator of the elastic bonds
1jfap has the form of(2.3.8). Repeating almost literally the arguments carried out
earlier for the one-dimensional case, one can introduce the three-dimensional
operators (/Jr!, raP(w), etc., write down the Green's formula and formulate the
basic boundary value problems. In particular, in the obvious notation (see the
notation [Ref. 1.12, Sect. 2.6]), the Green's formula has the form
(2.10.1)
(2.10.2)
(2.10.3)
(2.10.6)
(2.10.7)
where the functions Iffb Iff t and Iff are connected by the relations
I
/fft(r) = Sr-1/ff(r) dr . (2.10.8)
r
Problem 2.10.2. Show that the interaction of binary type implies the equality
of the Lame constants
271: I
AO = flo = -[5 Sr 4/ff(r) dr . (2.10.9)
o
Let us now assume that the medium occupies the upper half space and
choose the coordinate system so that the z-axis is directed normalIy to the
surface. The boundary conditions, corresponding to the absence of forces in
the surface layer S(z), have the form
Let us consider a plane surface wave propagating along the x-axis, i.e. let
us seek a solution of the form
where v(z) --+ 0 as z --+ 00, and p is a real parameter. It is easy to see that in
this case two types of surface waves can propagate independently in the medi-
um. For a wave of the first type the displacement vector lies in the plane, which
is drawn in the direction of propagation perpendicularly to the surface, and
which satisfies the boundary conditions
For a wave of the second type, the displacement vector is parallel to the z-axis
and satisfies the condition
(2.10.13)
Jl/fa2(r) dy = x a J ~ l/f(r) dy
r
= 0 (0: -:F 2) . (2.10.14)
au,
az = 0
. (2.10.15)
u, = O. (2.10.16)
Thus the wave of the second type does not have any classical analogy.
Let us proceed to the approximation with respect to the first root. The
equations of motion in this approximation take the form
(2.10.17)
(2.1O.IS)
At small values of lV, for the roots of the above equation we obtain
The root "0 corresponds to the usual Rayleigh wave and the roots "1,2
correspond to the new surface waves, which are absent in the classical theory.
With p ,.., lV/C, we have
2.11 Notes
The medium of complex structure is considered in detail with the main at-
tention being focused on specific three-dimensional effects.
Starting from the harmonic model of the crystal lattice with a basis, we
introduce the collective variables: the mass center displacement, micro rotations
and microdeformations. In terms of these variables, the general nonlocal equa-
tions of the medium of complex structure are obtained. The physical require-
ments of stability and invariance with respect to rigid body motions lead to the
existence of acoustical and optical modes of vibrations. It is shown that in the
corresponding acoustical region the microrotations and microdeformations
can be excluded from the equations of motion. As a result, the medium of com-
plex structure can be described in the acoustical region in terms of the medium
of simple structure with space and time dispersion (dynamical elastic moduli).
The detailed comparison of the theory with approximate couple-stress
theories is provided. It is shown that the latter have a very restricted range of
applicability. As an example of the medium of complex structure the multibar
system is considered in detail.
ea and the metric tensor ga{3 = ea ·e{3' the vector n has the integer component
In the absence of initial forces, the Lagrangian of the lattice has the form
2L = ga{3 ~m/.vaCn,j)win,j)
n;
where wa(n, j) is the displacement of the ph particle in the nth cell, and fa(n, j)
is the external force which acts on it.
Let us generalize the model in such a way that it would be possible to take
into consideration both discrete and continuous periodic structures. To this
end, let us introduce the density p(~) of the mass distribution within the cell
and a force matrix (/Ja{3(n, n', ~, e) which determines the interaction. Here,
~ = ~aea is the position vector of a point of a cell with the number n in the
local system of coordinates with the origin situated at the center of mass of
the cell. Then the Lagrangian is given by
- 1:
nn'
SJ wa(n, ~)(/Ja{3(n, n', ~, ~')w{3(n', e) d~ d~'
Using the formalism of the quasicontinuum let us write down the Lagrangian
(3.1.2) in the x representation
(3.1.6)
Here, s denotes the multi index ih ... As and the parentheses denote the cor-
responding scalar product.
Let us denote by P2; the tensor inverse to P2s and introduce the two quasi-
diagonal matrices
(3.1.7)
(3.1.8)
(3.1.9)
Problem 3.1.2. Taking into account the fact that the origin of the coordinate
system coincides with the mass center of a cell, show that for the first two
elements of the bi-orthogonal basis, the following relations are valid:
eO(';) = 1,
(3.1.10)
Let us expand the functions of the variables .;, .;' entering (3.1.4) in terms of
the bi-orthogonal basis
For the coefficients of the expansion, taking into account (3.1.8) we find
In the new set of variables, the Lagrangian (3.1.4) can be written in a form
which is invariant with respect to the n, x and k representations:
(3.1.13)
Jpas'fJ(x, x') ws,p(x') dx' + JJf/)saS'f3(x, x') WS'fJ(x') dx' = fsa(x) , (3.1.14)
where, in the general case, pas'f3(x, x') is the kernel of a nonlocal moment-of-
inertia operator.
Let us now consider the physical meaning of the collective cell variables
was(x) andfsa(x). In view of (3.1.10,12) we find that WOa is the displacement
of the mass center of the cell and fOa is the average body force density. It is
clear that these quantities will have a considerable significance when the medium
is described macroscopically. Hence, let us introduce the special notation
(3.1.15)
"pa
'Yl
-def Wpa'rIIPa - , p = 1" 2
deffPa (3.1.16)
Let
(3.1.18)
where p.().a) is the average density of force dipoles and m},a = - ma). is the aver-
age density of micromoments. It is clear that the p.().a) and m).a are generalized
forces, which correspond to the generalized displacements e5.a and Q;". The
quantities 1}pa and p'pa, for p > 1, describe the microdeformations and micro-
moments of higher orders.
Note that the number of independent elements of the basis is determined
by the number of degrees of freedom of a cell. If the number of degrees of
3.2 Energy Operator 49
freedom does not exceed 12 (in the case of a discrete structure, this means that
the number of particles in a cell does not exceed 4), then the basis is automati-
cally restricted to the first two elements. In this case, the role of kinematic
variables will be played by displacements of the center of mass u, micro-
deformation e' and micro rotation Q'; and the role of force variables will be
assumed by the body force density q, the dipole force density J1- and the micro-
moments density m.
If the cell consists of one particle with a finite moment of inertia, then the
internal degrees of freedom are restricted to a micro rotation Q'. This simplest
model of the medium of complex structure, known as the Cosserat medium,
will be considered in Sect 3.5.
In terms of the new variables, the equations of motion (3.1.14) take the form
where the operator coefficients are connected with the operators in (3.1.14)
in an obvious way.
From (3.1.13, 14) it follows that the properties of the medium are com-
pletely determined by the kinetic energy operators pas'fJ and the elastic energy
operators (/Jsas'{J. In this section we consider the structure of these operators.
According to a reasoning analogous to that of [Ref. 1.12, Sect, 3.3] let us
assume that the operator [sas'fJ is determined by a set of constant moment-of-
inertia tensors, the following ones having the principal significance:
(3.2.2)
(3.2.5)
(3.2.6)
(3.2.7)
(3.2.9)
Problem 3.2.1. Obtain (3.2.8, 9) from the laws of transformation of waCx, ~),
taking into account (3.1.10, 12).
</J0asf3(O, k') = 0,
(3.2.12)
ia[},</J0al sf3(O, k') + </J[},alsfl (0, k') = 0.
3.2 Energy Operator 51
This enables us to represent the matrix rpsas'fJ(k, k') in (3.1.19) in the form
[Ref. 1.12, Eq. (3.3.15)]
Problem 3.2.2. Prove (3.2.13). [Hint: Use the Hermitian adjointness and the
condition (3.2.12).]
The equations of motion (3.1.12), after taking into account (3.2.13), become
(3.2.17)
Here, rpo and rpl are the elastic energies, corresponding to the displacements
of the centers of masses and to the internal degrees of freedom and rpint is
the interaction energy. It is clear that the terms rpo and rpl are to be positive
definite forms.
Homogeneous Medium. As we know, in this case all the operators have dif-
ference kernels in the x representation and are reduced to a product of functions
of k in the k representation [Ref. 1.12, Eq. (3.9.7)]. The matrix (3.2.12) becomes
52 3. Medium of Complex Structure
(3.2.19)
When writing down the above matrix, the conditions of invariance of the
energy with respect to translation and rotation were used. However, additional
conditions, analogous to (3.2.15), are also to be fulfilled. One of them is obtained
immediately from the second relation in (3.2.15)
(3.2.20)
The simple form of the limiting transition is here due to the fact that the
tensors XPal'i3 and r Pa l'[3 do not possess any symmetry, which would be destroyed
by alternating with respect to the indices /1(3.
The analog of the first relation of (3.2.15) produces more difficulties. When
passing to the limiting case of a homogeneous medium, a new symmetry with
respect to the indices A/1 arises in the tensor r).al'i3, which does not commute
with the alternation with respect to /1(3. Therefore, as in the case of a medium
with simple structure, we have to obtain this condition directly for a homo-
geneous medium. The corresponding discussion was carried out in details in
Sect. 2.5. This enables us to restrict ourselves to a few principal points. Hence,
it may be worthwhile, with a view to brevity, to stress only the most salient
features.
We start with a requirement of invariance of the energy density of the
zeroth approximation ({Jo(x) with respect to rotation. Taking into account
(3.2.19), the density of the Lagrangian in the zeroth approximation can be
represented in the form
At last, analogously to the treatment given in [Ref. 1.12, Sect. 1.5], it can
be shown that the condition of solvability of (3.2.23) with respect to bJ.ap.{J
has the form
where
[Aa, ,u,81 = r~ap.[3 + (X~ap.[3 + X~{J).a + r~ap.(3)().p.) . (3.2.25)
Under the condition (3.2.24) the solution of (3.2.23) is unique
(3.2.26)
Problem 3.2.3. Write the three-dimensional analog of [Ref. 1.12, Eq. (3.2.35)]
for the case, in which the elementary cell has a centre of inversion.
Problem 3.2.4. Formulate the three-dimensional analog of the conditions
which were considered in 1.3.7 concerning the existence of an equivalent
model of medium of simple structure.
54 3. Medium of Complex Structure
In the general case the matrix (/)sas'~(k) does not possess any symmetry
with respect to a, (3. However, if this symmetry is present, then it can be given
a simple interpretation. Let us call the interaction in a discrete model binary
(trinary), if the elastic energy can be represented in the form of a sum, each
term of which depends on the difference of displacements of two (three) par-
ticles l . The generalization to the case of a continuous medium is obvious.
In [3.1] it was shown that, in the harmonic approximation, the trinary inter-
action is the most general form of interaction in a crystal lattice. It can be
proved [3.2] that the symmetry of (/)a~(n - n', ~, ~') and hence also of (/)sas'~
with respect to a, (3, is a necessary and sufficient condition of the interaction
being binary (in above-described sense). From this, in particular, it follows
that in anyone-dimensional model the interaction is always binary, and we
saw this also directly. In the three-dimensional models of simple structure,
invariance with respect to inversion
(/)a~(x) = (/)a~( - x) = (/)~a(x) , (3.2.28)
insures that the conditions of the interaction being binary, are fulfilled.
Let us point out an interesting case, when the conditions of invariance
of the equations of motion with respect to a homogeneous deformation
(3.2.29)
are satisfied. To this end, as follows from (3.2.27), it is necessary and sufficient
to satisfy the following conditions:
(3.2.30)
Together with (3.2.20) they give [Ref. 1.12, Eq.(3.8.7)]
(3.2.31)
In Sect. 3.4, we shall see that this condition considerably simplifies the com-
putation of the macroscopic tensor of elastic constants.
IThis definition is more general than the one given in Sect. 2.3, since the invariance of
the terms of the sum with respect to rotation is not required. In connection with this, the
definition has a formal nature.
3.3 Approximate Models and Comparison with Couple-Stress Theories 55
Let us expand the functions r(k), X(k) and F(k) entering (3.2.22) in powers
of k and consider only a finite number of terms. With k = 0, we have the
zeroth approximation, which is considered in the next section. When proceeding
to successive approximations, it is expedient to distinguish two cases (Sect.
2.6). For a gyrotropic medium, the next approximation is the first one and the
operators (3.2.17) have the form
r Aal'{3 = r6al'{3 + r~al'{3Ton Xpal'{3 = xPt,,{3 + Xfa,,{3TOr ,
(3.3.1)
FPap'{3 = F1fP'1' + F1ap'{3Tor .
As was pointed out above, this corresponds to the rotation of the polariza-
tion plane of a wave, which is typical of a gyrotropic medium. If the medium is
non gyro tropic then the approximation following the zeroth, is the second
order approximation where
r Aal'{3 = r6al'{3 + r~al'{3IlTOIlOT' XPa l'{3 = X~aJ1{3 + XP{I'{3IlTOIlOT ,
(3.3.2)
FPap'{3 = F~aP'{3 + FPz"p'{3llrOIlOT .
By the indices 1, 2 in (3.3.1, 2) we denote the corresponding coefficients in the
expansions of the functions r(k), X(k) and F(k) in powers of ik, the coefficients
obviously being real tensors.
Replacing the integral operators (3.2.17) by the differential operators (3.3.1,
2)[i.e. transforming (3.2.17) into (3.3.1)] or by analogous operators of higher
order, corresponds to the transition to a model of a medium in the approxima-
tion of weak dispersion. In other words, this means taking into account non-
local effects. From considerations of a dimensional analysis it follows that the
coefficients of the higher-order derivatives contain indirectly a scale parameter
I, which, due to the approximate character of the model, should be considered
small in comparison with a characteristic wave length or the size of the body.
Hence we conclude that the "resolving power" of the model is limited by the
parameter I, i.e. solutions of the approximate equations, in the best case, can
pretend to give a qualitatively correct description of non local effects in regions
of the order of I, but, generally speaking, should not be extrapolated to smaller
regions.
The equations of a medium with weak dispersion contain, as a particular
case, the equations of the Cosserat continuum, of an oriented medium, and of
couple-stress, multi poles, etc. theories of elasticity. The above considerations
point out the possible regions of applicability of these theories.
Usually, such theories are derived on phenomenological grounds, roughly
according to the following scheme. It is postulated, that the Lagrangian de-
pends on the strain and its derivatives and on a set of additional kinematic
variables such as microrotations and microstrain and their derivatives. By
analogy with the classical theory of elasticity the new kinematic variables and
their derivatives are put into correspondence with moment stresses of various
orders, for which equations of equilibrium (conservation laws) and equations of
56 3. Medium of Complex Structure
state of the Hooke's law type, are postulated. This enables us to obtain a closed
system of equations and a set of boundary conditions. Further formulation of
boundary-value problems usually does not differ from the formulation of an-
alogous problems in the classical theory of elasticity. For example, many papers
are devoted to investigations, in the scope of the couple-stress theory of elas-
ticity, concerning stress concentration around holes. In this situation, one would
expect that significant differences from the classical theory of elasticity arise for
holes with diameters of the order of or less than I.
Let us consider some of the difficulties which arise in such an approach.
1) All the couple-stress theories approximately account for the nonlocal
effects, and the parameter 1 which is contained in them has to be considered
to be small. Therefore, when extrapolating the results to the regions of the order
of or much smaller than I, one should be very cautious.
2) The presence of a small parameter in terms with higher-order derivatives,
as is well-known, leads to effects of the boundary-layer type. In connection with
this, it is hardly reasonable to consider boundary-value problems for com-
plicated regions in the scope of couple-stress theories. Moreover, it is also
unknown how to interpret the boundary conditions physically and from where
to take the values of the moments on the boundary.
3) The theories contain new constants (sometimes in large numbers) but
usually it is not clear, from what sort of experiments these could be determined.
It is necessary to know at least the qualitative relationship between these con-
stants and the parameters of the micromodel.
4) In an axiomatic construction of the theories without taking into account
their approximate nature, the equations of motion are often written down with
different degrees of accuracy with respect to different variables. We shall con-
sider this aspect with the example of the Cosserat model in Sect 3.5.
5) The connection between the couple-stress theories and the classical
theory of elasticity is usually discussed in a simplified manner: in order to pro-
ceed to the theory of elasticity, it is necessary to set the constants contained in
the equation equal to zero. As we have observed in the example of the one-
dimensional model of a medium of complex structure, the actual connection
is more involved: for sufficiently long wavelengths and for sufficiently low fre-
quencies the internal degrees of freedom are excluded from the equations of
motion, but they do make some contribution to the effective macroparameters.
In the next section this question will be considered for the three-dimensional
medium.
6) We observed that, in the nonlocal theory of a medium of simple structure,
the stress tensor is not a physically uniquely determined quantity. This holds
even more strongly for the couple stress in the medium of complex structure.
For its formally unique determination, it is necessary to postulate conserva-
tion laws,2 which do not have any group basis. As is well-known in field theory
[3.3], the conservation laws are a result of the invariance of the Lagrangian
with respect to the corresponding groups and in our case they are the groups
of translations and rotations. However, only displacements and microrotations
are essentially with this fundamental group; this leads to a special divergence
form of the equations of motion, which ensures the equilibrium of forces and
moments. The remaining internal degrees of freedom (microdeformations) are
invariant with respect to the fundamental group, and additional laws for couple
stress have to be postulated rather arbitrarily.
Note also that couple stresses do not satisfy the natural principle of cor-
respondence, i.e. they do not become the usual stresses in the limiting case of
long waves, when a transition to the conventional theory of elasticity should
take place.
Taking into account these arguments, one should not be surprised that the
rare attempts of experimental determination of couple stress in an nonhomo-
geneous medium have not led to any results.
Let us sum up. The couple-stress theories of elasticity and, generally, models
of a medium with weak dispersion approximately take into account compara-
tively delicate effects of nonlocality and of internal degrees of freedom. Their
application is expedient, first of all, when they yield not merely small quantita-
tive corrections, but describe qualitatively new phenomena. Therefore the
physically justified field of applicability of these theories is narrower than that
of the classical theory of elasticity, though the latter is the limiting case.
As in the one-dimensional case, we have the dispersion equation for free vibra-
tions
det { - (Ji /,as'f3 + (/)sas'f3(k)} = 0, (3.4.1)
which connects the frequency wand the wave vector k. However this equation
now yields three acoustic and 3N-3 optical modes w(k), where 3N is the number
of degrees of freedom of a cell.
As pointed out in [Ref. 1.12, Chap. 3], in the acoustic region offrequencies
the main kinematic variable is the displacement of the center of the masses of
a cell Ua = WOa. The distinctive role of this variable is connected with the fact
that in the acoustic region the internal degrees of freedom 1jpa can be excluded
from the equations of motion by transforming the latter into an equation,
which contains only one variable U a , and into equations which express explicitly
the remaining variables in terms of displacements. Recall that the possibility of
such an operation is connected with the special law of transformation of the
kinematic variables under translation and with the stability of the medium.
58 3. Medium of Complex Structure
Let us consider, for the beginning, the equations of the zeroth approximation
[Ref. 1.12, Eq. (3.6.1)]
Problem 3.4.1. Show that the optical frequencies at k = 0 are determined by the
equation
det( - w2JPaP'(> + FPryP'f') = 0, (3.4.3)
and, moreover, from the conditions of stability, the optical frequencies cannot
be zero.
We see that the matrix (- w2J + r 0) has an inverse in the region of acoustic
frequencies. This means that in the acoustic region the second equation of the
system (3.4.2) can be solved with respect to 7Jpl(>, and 7Jn can be eliminated from
the first equation. From here it also follows that there exists a finite neighbor-
hood of the origin of the k, w-space, in which the internal degrees of freedom
can be eliminated in such a manner.
In the direct matrix notations (3.2.22) coincide with [Ref. l.12, Eq. 3.3.24)]
for the one-dimensional medium. This enables us to describe the final results
without going into details.
Let us introduce the matrix
The neighborhood of the origin of the k, w-space, in which A(k, w) exists, will
be called an admissible acoustic region.
With the help of the matrix A(k, w) the system (3.2.22) can be transformed
into the form
(3.4.5)
(3.4.6)
on lV, or, in other words, there exists a time dispersion besides with the spatial
dispersion, the former being unconnected with the dissipation of energy.
Problem 3.4.2. Show that, when k and lV are small, the time dispersion appears
only beginning with the second approximation.
(3.4.8)
where cAal'P(k, lV) possesses the symmetry of the tensor of elastic moduli (2.3.19)
and is connected with rJ>AaI'P(k, lV) by the operator S. This operator is defined
by the set of operations described in Sect. 2.5
(3.4.10)
and write the equations of motion (3.4.8) in the (x, t)-representation in the usual
form
pu - div (J = q* . (3.4.11)
As distinct from the case of a medium with simple structure, the operator
Hooke's law has the form
(JAa(x, t) = H cAal'P(x - x', t - tl)CI'P(X', t') dx' dt', (3.4.12)
i.e. the strain tensor cl'P and the stress tensor (JAa have nonlocal connections not
only with respect to the spatial variables but also with respect to time.
It is important to emphasize that in the zeroth long-wavelength approxi-
mation, we obtain the equations of the classical theory of elasticity, and the
stress tensor (JAa transforms into the usual stress tensor, which is connected with
strain by the tensor of the elastic constants c«,I'P. The latter, according to (2.5.4),
is given by the expression3 [Ref. 1.12, Eq. (3.6.12)]
3For the crystal lattice, this expression coincides with that obtained by Born and Huang
[3.5, 6].
60 3. Medium of Complex Structure
Problem 3.4.3. Show that for the tensor cgal'(3 to possess the required symmetry,
it is necessary and sufficient to have symmetry of (/>}.al'(3 with respect to the
indices A.p. and aj3, as well as with respect to the permutation of these pairs of
indices, and this symmetry is ensured by (3.2.20).
m).al'(3 -
'Vo -
r).al'(3 _
0
1 (rl!.(3}.a
"2 0
+ r).(3tJa)
0, (3.4.14)
i.e. for the determination of caal'(3, it is not necessary to invert the matrix roo
Problem 3.4.5. On the basis of (3.2.17), show that in the long-wave length ap-
proximation the elastic energy (/) can be written in the usual form (with p. = 0)
(3.4.15)
(3.4.16)
Eq. (3.4.13) and the analogous formulae for the elastic constants of higher
approximations enable us to express these constants in terms of the parameters
of the micromodel. In the following section we shall illustrate this point with
the example of the Cosserat model.
When proceeding to the higher approximations in (3.4.8) it is expedient to
distinguish the cases of gyrotropic and nongyrotropic media. However, it is
necessary to note that the property of gyrotropy in the above sense is not
invariant with respect to the transformation of the equations carried out. A
medium can be gyrotropic with respect to the initial system of equations (3.2.22)
but then the first approximation in (3.4.8) may coincide with the zeroth one.
(3.4.17)
Here / and T are the characteristic length and time, Ao and po are the correspond-
ing coefficients of the expansion of the functions A(k, (1) and p(k, (1), and J can
be considered to be an inertial operator characteristic of the medium.
In the static case the equations for an isotropic medium coincide with the
equations for a medium of simple structure and, hence, for solving them, one
can use the static Green's tensor constructed in Sect. 2.8.
Thus, in the admissible acoustic legion, the equations of a medium with
complex structure are reduced to the equations for an equivalent medium with
simple structure with spatial and time dispersions. In the long-wavelength
approximation this corresponds to the transition to a macroscopic description
of the medium with microstructure, when the internal degrees of freedom are
excluded from explicit consideration, but contribute to the effective elastic
constants.
Note that, conversely, in the region of high frequencies it is possible to elimi-
nate the displacements of the mass centers from the equations, analogously to
what was done in the one-dimensional case [Ref. 1.12, Sect. 3.6].
Let us consider a medium for which the internal degrees of freedom are reduced
to microrotation. The latter can be characterized by an antisymmetric tensor
Q~~ or by the pseudo vector
(3.5.1)
62 3. Medium of Complex Structure
where in the cartesian coordinate system crf.'~ is the unit anti symmetric pseu-
dotensor. 4 Analogously, the force micromoments are given by the antisym-
metric tensor rn).a or by the pseudovector rn r • Eq. (3.2. 17), after obvious simpli-
fications, takes the form
More precisely, in order to obtain the usual equations of the Cosserat conti-
nuum [3.8], it is necessary to equate the constants {32' 72' Xz to zero in (3.5.3),
but to preserve the nonzer0 5 constant F z. The equations thus obtained
appear to be slightly inconsistent with respect to the order of approximation.
The first equation is of the zeroth order with respect to the scale parameter
I, while the second equation is of second order. This inconsistency is a con-
sequence of the "symmetrical" assumption, that the Lagrangian depends on
the displacements u, microrotations 0' and their first order derivatives.
In reality, the variables u and 0' enter the equations nonsymmetrically,
in view of the differences in the laws of their transformations with respect to
the translation group. Therefore, in order to have consistent accuracy, it is
necessary to assume that, in the Lagrangian, the displacement derivatives have
an order higher by one than the order of the derivatives of microrotation.
We shall not continue the general investigation of the Cosserat medium,
but shall instead proceed to a concrete model which is interesting for its own
sake. Let us consider a system of rods which forms a cubic lattice and assume
that, at the knots of the lattice, masses connected by the weightless rods are
located, and that the masses possess finite moments of inertia. Nearest-
neighbor interaction is assumed for simplicity.
Let us introduce the following characteristics of the rods being parallel
to the axis x a of the lattice coordinate system: length I, longitudinal stiffness
4In an arbitrary affine or lattice coordinate system this is the corresponding tensor density
[B8.251.
50ther possible variants are connected with different definitions of an isotropic medium.
3.5 Cosserat Model 63
+ 1: (/)ff~(n - n')uf3(n', t)
n'
(I)
P
= 12EJ ( 1 -
/3
1 )
-1I - _/2
1
i
()' ,
(u) (3.5.5)
2 6
where 1 and p. are the corresponding force and moment at the clamped end.
Taking into account analogous relations for tension and twisting, we
obtain explicit expressions for the force constants in terms of the parame-
ters of the rods. The nonzero components of the force matrices have the
form
64 3. Medium of Complex Structure
o 1
I~
1 0 0 0 001
-1 0 0 o -1 0 o 0-1
(3.5.6)
-B
+ I +B I
±B1 1 1
':fB z - ZI ±B3
(3.5.7)
Problem 3.5.1. Verify that the force constants of self-action, found from the
conditions of invariance of (3.5.4) with respect to translation and rotation,
see (3.2.10) and (3.2.11), have the form (do not use the summation convention!)
Problem 3.5.2. Write down the conditions of invariance with respect to rotation
(3.2.20, 24) taking into account (3.5.1), and verify that these conditions are
fulfilled for the system under consideration.
Problem 3.5.3. Write down the equations of motion (3.5.2) in the zeroth and
second approximations expressing the tensor coefficients in terms of the pa-
rameters of the system.
°
the form (3.4.5, 6). The tensor if) Aal'{3(k, w) entering this equation possesses
for k = and w = 0, the required symmetry with respect to the indices and
in particular, satisfies the conditions of invariance with respect to the rotation
(2.5.8). Recall that this condition was not fulfilled in the consideration of the
rod system (Sect. 2.7), in the framework of a medium of simple structure
(without taking into account microrotations).
Let us now present the final expressions for the operator of the elastic
moduli c Aal'{3(k, w) confining ourselves (for simplicity) to the case of identical
rods and setting Aa = A, Ba = B, C a = C. For the essential components of
cAal'{3(k, w) we have (a #- f3 #- r)
(3.5.12)
BB
aa - A a' CfIlP -
B a fJ co·aa·PfJ -- 0 (a..J.
.,.. t'f.l) (3.5.13)
+ B'
CO··aa - U··afJ - .
a fJ
For the rod system, let us consider the dispersion equation (3.4.1), which
connects the frequency wand the wavevector k. For some characteristic direc-
tions of the cubic lattice, in particular along the diagonals (1, 1, 1) and (1, 1,0)
and along the axis (1, 0, 0), it is possible to obtain the solution w(k) of the
dispersion equation explicitly. In all these cases degeneracy takes place: there
exists one acoustical and one optical mode, these being the analogs of longi-
tudinal waves, as well as two coinciding acoustical and two coinciding optical
modes, these being the analogs of transverse waves.
Let us write down the dispersive modes for the direction (1, 0, 0,), i.e.
for the wave vector (k h 0, 0). The acoustic modes are given by (p = VOl m)
w2 = wi(k 1) = :~ sin 2 ~ ,
(3.5.14)
w2 = w!,lk1) = ~ [I(k 1) - ,.; I2(k 1) - g(k 1)] ,
(3.5.15)
Here
(3.5.16)
Problem 3.5.5. Write down the solutions of the dispersion equation for the
directions along the diagonals.
3.6 Notes
Let 1/)0 be a linear operator for which the Green's function (inverse operator)
Go is known and satisfies certain boundary conditions, i.e., Gol/)o = I, where 1
is the identity operator. If 1/)1 is a perturbation of the operator 1/)0 such that a
Green's function G of the operator I/) = 1/)0 + 1/)1 exists, then one can show that
the Green's function G has the representation
G = Go - GoPG o , (4.1.1)
(4.1.2)
In fact, substituting Pinto (4.1.1) and multiplying by I/) from the left, we have
I/)G = (1/)0 + 1/)1) [Go - G 01/)1(1/)1 + 1/)1G OI/)1)-1 1/)1G O]
= 1 + 1/)1GO - 1/)1(1/)1 + 1/)1GOI/)1)-11/)1GO - 1/)1G01/)1(1/)1 + 1/)1GOI/)1)-1 1/)1GO
= 1 + [I - (1/)1 + 1/)1Go1/)1) (1/)1 + 1/)1GOI/)1)-1] 1/)1GO = I.
4.1 General Scheme 69
(4.1.3)
(4.1.4)
u = Uo - GoPuo, (4.1.5)
where the second term is the perturbation caused by the defects. It is obvious
that this equation is valid also in the case when the external field Uo is given
directly.
For the elastic energy (/), according to (2.3.22), we have
The first term corresponds to the energy (/)0 of the external field in the me-
dium without defects and the second term is the energy (/)int of interaction of
the defects with the external field and among themselves. This enables us to in-
terpret P as an operator for interaction energy.
Problem 4.1.1. Show that from the invariance of the energy with respect to tran-
slation, it follows that the operator P can be represented in the form
P = VPV. (4.1.7)
Substituting Pinto (4.1.6) we obtain the expression for the energy of inter-
action in terms of the deformation eo of the external field and the operator 15:
70 4. Local Defects
(4.1.8)
where the summation is carried out over the region of perturbation and the
matrix plp(mm') is defined by
plp(m, m') = {[rpfP(m, m')]-l + (flJ.a(m - m')}-l . (4.2.3)
Since the points m, m' here belong to the finite region of perturbed bonds, the
problem is reduced to the solution of a finite system of linear algebraic equations.
When inverting the matrices, it is necessary to take into account that, due
to the conditions of invariance with respect to translation and rotation, the in-
verse matrices exist only in some special vector subspaces, which can be distin-
guished by constructing the corresponding projection operators. Also, in the
general case, expressions ofthe type (4.1.2) should be used instead of (4.2.3).
Thus, an effective analytical solution of the problem can be constructed only
for the simplest models, which we shall not dwell upon. In more complicated
situations, numerical methods are needed.
In analogous manner (by numerical methods) the problem of the dynamic
Green's tensor Gain, n', ev) can be solved for the case of a single defect as well
as for the case of finite set of defects.
In order to obtain tractable analytical expressions, it is expedient to accept
another model for the defect, which will be considered in the next section.
quasicontinuum. This enables us to use the a-function model for the description
of local inhomogeneities, introducing a-functions and their derivatives directly
into the coefficients of the equations. Such a procedure would be incorrect for
an elastic continuum.
Suppose that there is a single defect at a point x = Xo. Let us first consider
the static case. The simplest model of a defect of the type of local inhomogeneity
of the operator for the elastic moduli in the neighborhood of the point Xo can
be described by an expression of the form
(4.3.1)
Here, the constant tensor cia i'{3 has the dimension and symmetry of the elastic-
constants tensor and characterizes the variation of the elastic properties in the
neighborhood of the point Xo. A more complicated model of a defect can be
obtained if we add additional terms with derivatives of a-functions.
Straightforward computations according to the scheme of Sect. 4.1 give for
the Green's tensor Ga{3(x, x') of the medium with the defect, the expression
(4.3.2)
The constant tensor PVi').~ appearing here has the symmetry of the tensor of
elastic constants and is determined by
P Vi').~ -- _ (....0
l> Vi').~
+ v-lb
0 Vi').~
)-1 , (4.3.3)
(4.3.4)
The existence of g2i').~ follows from the analyticity of G~ix) for the quasi-
continuum. From here it also follows that Ga{3(x, x') is an analytic function of
x and x'. However Gaix, x') depends on the parameter a nonanalytically.
Problem 4.3.1. Show that gO '" VOl, P '" Vo and, as a ~ 0, the second term in
(4.3.2) is of the order of Vo for the points (x "# xo, x' "# xo), of the order of a
for points (x = xo, x' "# xo) or (x "# xo, x' = xo), and tends to infinity as {l1
for the point (x = x' = xo) .
The construction of the dynamic Green's tensor for a single defect can be
carried out in an analogous way. Here it is reasonable to consider two cases.
If the density p(x) = Po = const and if the defect is caused only by a change
of the elastic properties of the medium, then the expression for the dynamic
Green's tensor Gaix, x', w) coincides with (4.3.2), after replacing G~ix) by
72 4. Local Defects
the dynamic Green's tensor for the homogeneous medium G~J3(x, w). Similarly,
p = P(w) is determined by (4.3.3) with gO = gO(w).
In the other case, when c(x, x') = co(x - x'), and the density has the form:
In the general case of a defect of mass and elastic moduli, the expression
for the Green's tensor is constructed in an analogous way, but has a more
complicated structure.
Let us now consider a perturbation of the external field u~(x), caused by the
defect, and let us confine ourselves to the static case. Taking into account (4.3.2),
we find
where u~(x) is the perturbation which coincides, as can be easily proved, with
the displacement in the homogeneous medium, caused by a force dipole at
the point Xo with density
(4.3.9)
(4.3.10)
(4.3.11)
introduce naturally the concept of a force f, with which the external field acts
on the defect. By definition
(4.3.12)
Cl = 2P.lE1 + AIE2,
Ela!'fJ def ~ (aafJaJ.!, + aaJ.afJ!'), E~afJ!' def aJ.aa!'fJ . (4.3.13)
Using the general expression for the Green's tensor of the Debye quasicon-
tinuum with spatial dispersion (2.9.11) and taking into account (4.3.3, 4 and
13) we obtain
(4.3.15)
where
r(3/ I 1-1
PI = V't 3(Ao + 2p.o) + 3Al + 2P.l J '
(4.3.16)
P _ v I.l.. ( 2{3/ + 3{3t) + I 1-1
2 - 't15 Ao + 2p.o Po
2p.1 J·
k4 .
(3j
3v?
= .1t'~ Jwr(k)
K
dk (J = I, t) . (4.3.17)
For the Debye model and for the model (2.9.8) we have, respectively,
(4.3.18)
Let us assume in addition that with the defect there is associated a force
dipole with the defect there is associated a force dipole with the density 1
1 Recall that for an isotropic medium we use the notations rand r instead of x and 1x I.
74 4. Local Defects
For the moment of the effective dipole, according to (4.3.10) and taking
into account (4.3.15, 16 and 20), we have
(4.3.21)
The expression for the total deformation c = cO + c1 now has the form
(4.3.22)
(4.3.23)
(4.3.24)
(4.3.25)
Va be equal to Ao + Al and f.J.a + f.J.l and let, on its boundary, radial forces with
the surface density Q/vo be applied. Then it can be shown that for the Debye
model (i.e., when spatial dispersion is absent) the values of e(O) and 0'(0) coincide
with the strain and stress within the sphere and the asymptotic values of e(r)
and O'(r) coincide with the strain and stress outside the sphere. The total elastic
energies for these cases also coincide exactly.2
For the interaction energy <Pint of the defect with the external filed egir),
we find, according to (4.3.11, 15)
<Pint = J
~l [ Tr{eO(O) } + ~{Dev{eO(O) } J' (4.3.26)
where Tr{eO} and Dev{eO} are the spherical and deviatoric parts of eO.
From this it follows that the force which acts upon the defect is equal to
Omitting simple but lengthy computations, we present the final result for
the static Green's tensor
The matrix P ij with tensor components entering the above expression, is the
inverse of the matrix
(4.4.3)
where
(4.4.4)
2Such a coincidence is, of course, caused by the high symmetry of the model. In the general
case, the correspondence between point defects and local inhomogeneities is more complicated.
We return to this question later.
76 4. Local Defects
The dynamic Green's tensor for the system of point mass defects has an
analogous structure.
In many cases, the dominant term of the asymptotics of the perturbed
elastic and or interaction energy is of principal interest. It may be shown that,
in such a situation, the problem is reduced to the binary interaction of defects.
In this connection, let us consider in more detail the case of two defects, which
is by itself of interest.
When i andj = 1,2, the computation of the components P ij is conveniently
carried out according to (no summation)
Pi .. = (R jj - RijRj/Rj,)-l, (i :f,j),
(4.4.5)
Pij = (Rji - RjjRj} Rjj)-l
p.ip) = P7 + O(p-6) ,
(4.4.6)
p .. /p) '" p-3 + O(p-9) ,
where the tensor PP corresponds to an isolated defect at the point Xj and is de-
fined by (4.3.3).
For the isotropic Debye quasicontinuum, the matrix P../p) can be computed
explicitly. Let us also present the generalization of (4.3.8) to the case of a
system of defects in the external field u~(x):
Here, u~(x) is the displacement due to a defect at the point Xj' However, the
moments of the equivalent force dipoles
Q"I'
• = .L.,.l'J ",0 (x.)
" P¥f!.).< ,,).< J (4.4.8)
J
now depend on the values of the field at all points x j as well as on the elastic
characteristics of all defects. Therefore, in the general case, the perturbation is
not the superposition of perturbations caused by individual defects, i.e. the
matrix Pij depends on the perturbations of the elastic constants of individual
defects in a complicated manner. In the limiting case of large distances between
the defects it follows from formulae of the type of (4.4.6) that in the zeroth
approximation, the perturbations of individual defects can be linearly super-
posed; in the next approximation we have the superposition of binary inter-
actions.
4.5 Local Inhomogeneity in an Elastic Medium 77
(4.4.9)
The force, which acts on the point defect at a point Xk is given by the
expression
rk def 0 (/J
J,l ~ - oxt int
For defects situated far enough from each other, the second term becomes
vanishingly small and the first term coincides with the force which acts on the
isolated defect at the point Xk' Conversely in a homogeneous external field
cO = const, only the second term contributes to the force.
(4.5.1)
where cff{l,lf.l is the tensor of the elastic constants of the homogeneous medium,
c7{l,lu is the perturbation caused by the defect localized in a (small) region V;.
The case Cj -+ 00, corresponds to a rigid inclusion and the case C; -+ - Co
corresponds to a cavity. Let us consider, for simplicity, an infinite elastic
medium.
Using the general scheme described in Sect. 4.1, one can show without diffi-
culty that the Green's tensor can be represented in the form
G = GO - GOVPVGO, (4.5.3)
The matrix pjj and its inverse matrix Rij have the structure analogous to (4.4.3):
It can easily be seen that the operator P is self-adjoint; its kernel satisfies the
symmetry conditions
(4.5.6)
(4.5.8)
Taking into account that P(y, y') implicitly depends on the coordinates of
the centers of mass of the defects, we have for the force fj, which acts on the
kth defect, due to the external field and other defects,
(4.5.9)
In an analogous manner one can find the moment which acts on a defect.
In the general case, the kernel P(y, y') may be found by numerical methods,
and use of the latter is facilitated by the fact that, unlike G(x, x'), P(y, y') is
localized in a bounded region.
The problem is considerably simplified and, in a number of cases, can be
solved analytically if the defects are located at distances, which are large in
comparison with their size and if only the asymptotic form of the perturbed
field is of interest. Eventually, this is equivalent to replacing local inhomo-
geneities by equivalent point defects.
First, let us consider approximations with respect to the distance r jj =
IXj - xjl between defects. In the zeroth approximation the defects do not in-
4.5 Local Inhomogeneity in an Elastic Medium 79
teract and Pij = Piou j , where pi is the operator of the i-th defect. One can show
that in the first approximation the problem is reduced to binary interaction.
Therefore let us restrict ourselves to the case of two defects. For two defects,
according to (4.5.5), we have
(4.5.10)
(4.5.11)
'" - -2
YJ'int -
1 EOa,8 pa{3A"(r
()() ij)EOA,.., P 00 -_ 4~ pij
., 00· (4.5.13)
From here it follows that only the asymptotics of the matrix ptIJ(r.) contri-
butes to the dominant term of the asymptotics of ~int in a homogeneous ex-
ternal field. In the case of an arbitrary external field, it is necessary to take into
account also the contribution to the asymptotics of ~int from the asymptotics
of the diagonal components of the matrices pM(rij) and pf6(r.). This contribu-
tion is equal to zero if the defects possess central symmetry or if the external
field changes slowly over a distance of the order of the size of a defect.
Let us sum up. For the description of the interaction of defects in an elastic
medium, it is necessary to know, in some approximation, the Green's function
for the medium with defects, or, equivalently, one has to be able to solve the
80 4. Local Defects
The equations of classical elasticity were written above in the form of (2.6.4).
In this section we confine ourselves to statics. Then the equations become
(4.6.1)
u = Gq, (4.6.2)
where
(4.6.5)
For the isotropic medium, the tensor of elastic moduli becomes, see (2.8.4),
where E1 and E2 are constant isotropic tensors defined in (A. 1), and )., p. are
4.6 Homogeneous Elastic Medium 81
the Lame constants. The expression for G(k) coincides with (2.8.6), where A
and p. are taken to be constants.
For an arbitrary anisotropic medium
where na = kallkl and gap(n) is a function defined on the unit sphere in k-space,
which is a determined explicitly by (4.6.4, 5).
The explicit expression for Gap(x) is known only for an isotropic medium
and media with transversly isotropic and hexagonal symmetries [4.1]. However,
we shall see that in a number of important cases it is sufficient to know general
properties of Gap(x) and its expression in the k-representation.
In particular, note that it follows from (4.6.4, 5) that
(4.6.8)
Let us consider briefly the behavior of the solution (4.6.3) for a few important
types of the density of external forces qa(x). Let qa(x) = ija(x) V(x), where V(x)
is the characteristic function of a bounded region V with a sufficiently smooth
boundary (), and ija(x) is a continuous function. Then uix) is a continuous
function, which is twice differentiable inside V, analytic outside V and tends to
zero as ,-1, when r = Ixl --+ 00. Its first derivatives and consequently, the strain
and stress tensors, have a finite discontinuity on ().
If qa(x) = qffo(x), then the solution is proportional to Garf.-x) and has, as
follows from (4.6.6), a singularity of the order of ,-1 as x --+ O. From here, it
is easy to find also the behavior of uix) for external forces concentrated on
surfaces.
Proofs of all these statements can be obtained in a manner analogous to
that used in the potential theory, we shall not therefore dwell on this.
We shall see later, that in a number of applications, in particular for a
medium with defects, it is convenient to have an equivalent form of the equations
of elasticity and to have their solutions not in terms of the displacement ua ,
but in terms of the strain eap or stress (laP. For this purpose, let us rewrite (4.6.1)
in the form of the usual system
(4.6.9)
(4.6.11)
where pal' = pl'a is the density of force moments, the latter being defined non-
uniquely. After some additional conditions are imposed, the correspondence
qa _ pal' becomes one-to-one, and it is not difficult to find pal' explicitly. It is to
be emphasized that, in problems concerning inhomogeneities in an external
field, pal', and not qa is the naturally given quantity (see below).
Substituting into (4.6.3) the expression (4.6.11) for qa and applying the opera-
tor defto both sides of the equality, we obtain
(4.6.12)
where
(4.6.13)
In operator notations,
From the theory of generalized functions, it follows [4.2] that the inverse
Fourier transform K(x) of the generalized homogeneous function of the zeroth
degree K(k), is defined with the help of the corresponding regularization. The
latter is reduced to a decomposition of K(x) into a regular part and a singular
one, which is proportional to o(x). Properties of the generalized function K(x)
are considered in detail in Appendix A.3.
From (4.6.12), it follows that the desired strain field is obtained as a con-
volution of the function K(x) with the density of force moments, p(x), i.e. K
is the Green's operator for strain. Let us consider a few features ofthis operator.
From (4.6.14), it is obvious that the operator K is annihilated by all opera-
tors which annul the operator def. In particular, this is true for the operator Rot,
4.6 Homogeneous Elastic Medium 83
RotK= O. (4.6.17)
Further, taking into account (4.6.4, 14) and the symmetry of the tensor of
elastic constants, we have
or
KCK= K. (4.6.19)
U= KC, li = CK (4.6.20)
(4.6.21)
while
li+ = U. (4.6.22)
m =Bp, (4.6.23)
i.e. the projections U and li give one of the representations of the Green's
operators for strain and stress, respectively.
This representation is closely connected with the elastic energy (/). In fact,
taking into account (4.6.19), we have
The connection of this operator with internal stress will be found in the
next chapter. Here, we consider only the relations which are analogous to those
for the operator K presented above.
SBS = S. (4.6.30)
e = SB, e = BS (4.6.31)
e+ = e. (4.6.32)
(4.6.35)
Thus, the strain and the stress can be obtained by projecting the field
sources onto the corresponding subspaces of the Hilbert space. To the pro-
4.6 Homogeneous Elastic Medium 85
jections li and ff, there correspond the usual stress a (caused by external forces
q, i.e. diva = -q) and the usual strain c (which satisfies the compatibility
conditions Rot c = 0 and is represented in terms of the displacement, as c =
def u ). Conversely, the internal stress a and the internal strain c, for which
diva == 0, but Rot c#-O correspond to the projections and e e. In this case,
the deformation c is not represented in terms of a displacement. A detailed
consideration of the internal stress and the corresponding field sources will
be given in the next two chapters.
Unfortunately, the projection operators introduced above are not self-
adjoint, i.e. the corresponding projections are not orthogonal. This leads to
some mathematical complications, since arbitrary projections do not possess
a number of nice properties inherent to orthogonal ones. This difficulty can
be easily overcome by a renormalization of the field variables, after which
the projections will become orthogonal. The loss of the usual physical meaning
of the field variables, is an inevitable price.
Positive definite tensors such as C and B = C-l, have uniquely determined
self-adjoint square roots C' and B'. Let
7: = Hp., (4.6.37)
where
HZ = H, H+ = H, (4.6.39)
(4.6.40)
ez = e, e+ = e, (4.6.42)
while
86 4. Local Defects
1=0+8. (4.6.43)
Thus, the solution of the problems of elasticity for both external and internal
stress is reduced to the construction of the orthogonal projections 0 and 8
in the Hilbert space with the scalar product, defined by the elastic energy
(4.6.40). The projections were explicitly expressed above by means of the
Green's operator G for displacement. It can be shown that the converse is
also valid, namely the operator G can be expressed in terms of 0 and 8.
It is essential that if one admits the definition of the operator K in the
form of (4.6.14), then all the relations, starting with (4.6.17), remain valid
for an arbitrary inhomogeneous elastic medium. In connection with this,
a closed system of equations for the projections is of interest. This system will
be deduced in Sect. 4.8.
Before we proceed to boundary value problems for elastic media with in-
homogeneities, let us consider a fundamental local problem about a discon-
tinuity of stress and strain at a point on the interface between two media.
Let two elastic media with tensors of elastic moduli cSIlA,u and callA,u be joined.
Assume that, on the (sufficiently smooth) interface {} between them, the
usual conditions of continuity of displacement and of the normal stress vector
(4.7.1)
are satisfied. Let us find an expression for the discontinuity of the stress and
strain at a point on the interface {}.
From the continuity of u on {}, follows the continuity of the tangential
component of the tensor Vu on {}. Introducing the projection operators onto
the normal and onto the tangent plane
(4.7.2)
(4.7.3)
(4.7.4)
where replacing CAl' by ihu,u is possible due to the symmetry of the tensor of
elastic moduli. Decomposing VtJ! into the sum of normal and tangent com-
ponents, we obtain
4.7 The Interface of Two Media 87
(4.7.6)
where
(4.7.7)
(4.7.8)
Adding this to (4.7.3) and symmetrizing with respect to the indices A, fl., we
have
(4.7.9)
where
or, equivalently,
K - K o + K[ C]Ko = 0. (4.7.13)
Let us present explicit formulae for the discontinuity [c] of the strain on the
interface between two isotropic media for the case, when the x 3-axis coincides
with the normal:
In this simplest case the formulae for [e] can be obtained directly, by obvious
reasoning. The continuity of the tangent component of the tensor Vu follows
from the first condition in (4.7.1). In the coordinate system chosen, this means
the continuity of six components of oau{3 with indices la and 2f3 and, hence,
of the three components of strain ell, e12 and e22' The continuity of stress
components 0'13,0'23 and 0'33 follows from the second condition in (4.7.1). Using
the derived conditions on 0' and e and also Hooke's law for an isotropic medium,
we arrive at (4.7.14).
Calculations of such a kind, in a fixed basis, have an essential defect, namely
they do not enable one to find an invariant meaning of the coefficients in the
formulae for the discontinuities. We shall see how useful it is to have a physical
interpretation of these coefficients.
Notice. The formulae obtained for the discontinuity of the strain are a
consequence of the interface conditions (4.7.1), but are not completely equiv-
alent to the latter, since they are valid under weaker assumptions. In particular,
for the case of the isotropic media considered, the component u3 can have a
discontinuity on the interface that will lead to a a-function term in e33' However,
this a-function will not affect the discontinuity in the limiting values of e33'
It is not difficult to obtain formulae for the discontinuity of the stress
from the formulae for the discontinuity in the strain. We present these formulae
in a form, which is completely symmetric to (4.7.11, 13),
S - So + S[B]So = O. (4.7.17)
Here
(4.7.20)
where
(4.7.21)
Let us agree that, when considering a defect in an inhomogeneous medium,
we shall denote by C and Co the moduli of the medium with and without the
4.8 Integral Equations for an Inhomogeneous Medium 89
(4.7.22)
(4.7.23)
holds.
Note, that C 1 = - Co or Bl -+ 00 corresponds to a cavity whereas C 1 -+ 00
or Bl = - Bo corresponds to a rigid inclusion.
In the notations adopted above,
(4.7.24)
Problem 4.7.1. Verify that the analogs of (4.6.19, 30) are particular cases of
(4.7.13, 17). (Hint: Set formally C = -Co.)
Problem 4.7.2. Show that (4.7.13, 17) are equivalent to the relations
where the projection operators are defined by (4.6.20, 31, 38 and 41).
(4.8.1)
(4.8.2)
(x E V) (4.8.4)
(x E V). (4.8.5)
where
(4.8.7)
(4.8.8)
where Uo = Goq is an external field which vanishes at infinity, i.e. the field
which would be present in the homogeneous medium under the action of
the force q. By virtue of our assumptions, u(x) and uo(x) are continuously
differentiable functions.
Let us apply the operator def to both sides of (4.8.8). Taking into account
the symmetry of Cb we obtain the equation for the strain t = def u
(4.8.9)
is the Green's operator for strain in the homogeneous medium, which was
introduced in Sect. 4.6.
(4.8.11)
So = Co - CoKoCo (4.8.12)
is the Green's operator for internal stress introduced in Sect. 4.6, see (4.6.29).
Let us consider now, how far one can extend the field of applicability of
(4.8.9, 11). We set temporarily
e=to+tl (4.8.13)
(4.8.14)
92 4. Local Defects
The right hand side contains only the restriction 4 of the external field
co to the support V of the function C1(x) with bounded support. Thus C1 will
not change, if one extends cb outside V in any (smooth) way. Therefore,
the external field co appearing in (4.8.9) can satisfy any conditions at infinity,
c - t co as Ixi - t 00 being automatically satistied. The same conclusion is also
valid for equation (4.8.11).
Let us now assume that C1(x) has a discontinuity on a certain surface D,
and that co is continuous on D. Then, by virtue of the properties of the operator
Ko (Appendix A.3) the solution also has a discontinuity, which satisfies the
conditions (4.7.11), i.e. the necessary conditions at the interface are fulfilled.
This enables us to extend (4.8.9) to piecewise continuous C1(x)'s, in particular,
to functions of the type of a product of a constant tensor and the characteristic
function of a region. Consequently, (4.8.11) can also be extended to piecewise
continuous functions B1(X),
Let us proceed to the extension of (4.8.9, 11) to the case of functions C1(x)
and B1(x) of the oscillatory type. We confine ourselves to the most interesting
case, when C1(x) and B1(x) are almost periodic functions, which are repre-
sentable as series of exponentials with arbitrary wave vectors. Let the external
fields be functions of the same type. Then, the solutions of the equations and,
consequently, the products C1c and B1c belong to the same class of almost-
periodic functions. Thus, we arrive at a problem of extending the operators
Ko and So to functions of the above-mentioned type.
Since the series are understood in the sense of generalized functions, the
action of the operators on exponentials can be taken to be termwise. The
action of an operator on an exponential with a nonzero wave vector is written
in the k representation in the form of a product of the corresponding homo-
geneous function and a a-function with support at a nonzero point. This is
a well-defined operation. Hence it remains only to complete the definition of
the action of the operators on constants, which corresponds to multiplying
the a-function concentrated at the zero point by homogeneous functions,
the latter being a quite nonunique operation.
In order to understand how to overcome this difficulty, let us consider a
model problem. Let an "inhomogeneity" be given by the tensors C1 or B1
which are constant throughout the whole space. Assume first that what is
really given is 0"0 = const. Then it is obvious that
(4.8.15)
It is natural to require that (4.8.9, 11) give the same result. This is possible only
if Ko and So act on constants as follows:
JKo(x - J
x') dx' = B o, So(x - x') dx' = 0, (4.8.16)
(4.8.19)
(4.8.20)
where K and S are naturally interpreted as the corresponding Green's operators
for a medium with characteristics Cor B while - C 1 and - Bl can be considered
as "perturbations."
Let us emphasize the correspondence between the pairs of integral equations
(4.8.9, 19) and (4.8.11, 20) on the one hand and the pairs of local relations at
the discontinuities (4.7.12, 16) on the other hand. This correspondence is
not accidental. The conditions for the jumps can be obtained from the integral
equations, making use of properties of the Green's operators. Conversely, the
conditions for the discontinuities imply a natural form of the integral equations. 4
3There is an analogy between these asymptotic conditions and the first and second bound-
ary value problems.
40 ne can say that there is a remarkable finite-dimensional representation of the global
Green's and projection operators by local ones.
94 4. Local Defects
K - Ko + KC1Ko = 0, (4.8.21)
(4.8.22)
(4.8.23)
are the global analog of (4.7.25). One can prove that these equations together
with the relations
cn = nc, ce = BC (4.8.24)
n+ = n, e+ = B (4.8.25)
and the analogous ones for no, eo, determine the projections nand e uniquely,
if the projections no and eo are given. For the latter, one can take the standard
projections Do, B o, which correspond to the case Co = l. Their explicit expres-
sion is presented in (B.7). Then (4.8.23-25) uniquely determine the Green's
operators and, hence the solutions of the corresponding problems for an in-
homogeneous medium.
Problem 4.8.3. Verify that the system (4.8.23-25) is equivalent to the system
see (4.7.26).
In conclusion, let us consider briefly the case, when C1(x), and consequently,
Bl(X) have a bounded support (a local inhomogeneity). Let V(x) be the charac-
teristic function of a bounded region V occupied by the inhomogeneity, and V
be the complement of V. Let us denote by + (-) the restriction of a field to
the region V (V) and introduce the operators
(4.8.28)
(4.8.29)
where the first equation uniquely determines the solution inside V and the
second uniquely determines its continuation to V. Analogously, for (4.8.11),
we have
(4.8.30)
(4.8.31)
Problem 4.8.4. Show that, for a cavity, (4.8.9, 11) take the form, using (4.8.4),
(4.8.32)
(4.8.33)
and for a rigid inclusion, (4.8.9, 11) take the form, using (4.8.5),
(4.8.34)
(4.8.35)
Note that, in the case of a cavity, c+ is finite and (1+ = O. Conversely, for a
rigid inclusion, (1+ is finite and c+ = O.
Let us now consider a field outside the inhomogeneity, e.g. for definiteness
the stress field. According to (4.8.31), for the perturbed field (1l(X), when x E V,
we have
(4.8.38)
(4.9.1)
It follows from the definition (4.8.27) of the operator KJ+) and from the
property (C.I) of the generalized function Ko(x) that KJ+) is decomposed into
two parts: the first one is proportional to the identity operator and the second
one is determined by the regularized part of Ko(x), an arbitrary region contained
in the regularization being naturally identified with V. As a result, it is sufficient
to consider the regularized part only. The latter acts on a function q>(x) which is
smooth within V according to the rule
where Ko(x) is the formal second derivative of Go(x) defined by (4.6.13). Taking
into account that Go(x) """ Ixl-1 and Go( -x) = Go(x) we can rewrite (4.9.2)
in the form
where
where
*
Let us set dx' = r 2drdwand integrate first over an elementary cone dw(n)
with vertex at the point x E V. Then
J(k)(X) =
{J)
Jgo(n)n(k)pk(x, n)dw, x E V, (4.9.8)
98 4. Local Defects
where
is the distance from the point x in the direction n to the boundary of the unit
sphere V. For pk we have
k
pk(X, n) = ~ (-I)ICL(x.n)l [.vI-x.x + (x.n)2] k-I. (4.9.10)
Since go(n) = go( -n) the contribution to J<kl(x) is due only to terms whose
product with n(kl is even. It is readily seen that they contain the radical in an
even degree and have the form
where
is a constant tensor.
Substituting (4.9.11) into (4.9.6) we find that F(ml(x) is a polynomial of
degree m. Finally, it is easy to show that the term x(ml has a nonzero coefficient.
Q.E.D.
Note 1. To extend the theorem to the cases of a cavity and a rigid inclusion
we have to make a further assumption. It follows from the proof that for the
cavity (C 1 = -Co) a weaker statement is valid: the operator 1- KJ+)Co
transforms a polynomial c:+(x) into a polynomial C:b(x) of the same degree but
equation (4.8.28) is not solvable for all c:b(x). This means that the theorem is
valid if and only if conditions of solvability are fulfilled. It is essential that for
the most important case, C:b = const, (4.8.28) always admits a solution as will
be shown below. The extension of the theorem to the case of a rigid inclusion
is performed in the same way but the corresponding field variables will be
stresses satisfying (4.8.30).
Note 2. It follows from the proof that even (odd) c:o(x) corresponds to even
(odd) c:+(x). In particular, if c:o(x) is linear then c:+(x) is also linear.
Note 3. Using the Green's functions for homogeneous elliptic operators of
the order 2p [4.2] it is easy to generalize the theorem to the case 2p < n where
n is the spatial dimension, n being arbitrary for p = 1.
4.9 Ellipsoidal Inhomogeneity 99
(4.9.13)
(4.9.14)
We define the mean value off(n) over the unit sphere n2 = I by the expression
<f(n) = ~ J
det a f(n) p3(n) dn . (4.9.16)
Problem 4.9.2. Show that <f(n) coincides with the mean value of f(x) over
the ellipsoid D defined by the expression
_ I f(x)
<f(x) - 4ndeta! p(x) dx, (4.9.17)
which permits one to interpret <f(n) as the mean value over the ellipsoid D.
Now, let us consider consequences of the theorem. Let the external field
co(or 0"0) be homogeneous, i.e. constant. It follows from the theorem that in
this case c+ = const and 0"+ = const. To calculate them explicitly we use (4.8.28)
and (4.8.30), rewriting them in the form
(4.9.18)
(4.9.19)
100 4. Local Defects
It is seen from here that we need to know the action of the operators KJ+) and
SJ+) on constants. Let us recall that these operators are represented in the form
where A and D are constant tensors and the tilde denotes the regularized parts;
these and the constants A and D depend on the choice of the regularization
region. In the case under consideration it is natural to identify the regulariza-
tion region with the ellipsoid. Then we have for the constants A and D (Ap-
pendix C)
(4.9.21)
where the integration is carried out over the unit sphere ()} = 1. Introducing
the new variable
(4.9.22)
and taking into account that Ko(k) and SoCk) are homogeneous functions of
zero degree, we find
It is readily seen from (4.9.2) that in the case e+(x) = const the contribution
of the operator KJ+) is equal to zero. Hence, the action of KJ+) on a constant
results in multiplication of this constant by A. In the same way the action of
Sd+) results in multiplication by D. Finally, we have
(4.9.24)
(4.9.25)
Thus, to find e+ and ()'+ we must invert the corresponding constant matrices.
It can be shown that the inverse matrices always exist including the limiting
cases of the cavity and the rigid inclusion.
Note that (4.9.23) permits one to rewrite the equations for e+ and ()'+ in the
form
(4.9.26)
(4.9.27)
A comparison with the discontinuity relations (4.7.12, 16) shows that the
coefficients in the equations admit an interesting interpretation: they are equal
4.9 Ellipsoidal Inhomogeneity 101
to the mean values over the ellipsoid of the coefficients in the corresponding
relations for jumps.
Now let us consider the case for which the external field is a linear function
of x
(4.9.28)
(4.9.29)
the tensors band d being linearly dependent. It is convenient to write this de-
pendence in the form
(4.9.30)
The solution of the problem for linear stresses can be deduced in the same way.
If the external field is a polynomial of the second or higher degree, then the
p-property permits one to find from (4.9.18, 19) recurrent algebraic equations
which connect the coefficients of the polynomials eo(x) and c+(x) or (7o(x) and
(7+(x). However, in contrast to the case ofthe homogeneous field these equations
are solvable only under additional conditions. The explanation is as follows:
formally there exist external fields (they do not satisfy the equilibrium condition
div (70 = 0) which are not perturbed by the ellipsoidal inhomogeneity. For ex-
ample, in the isotropic case this property has an external field of the type
ca{3 = (c,x)oa{3 where C is an arbitrary constant vector. As a consequence the
matrix A can be inverted only on the corresponding subspace.
Let us consider one more problem which is of interest for applications.
Let (7(x) be the limiting value, from outside, of the stress at the point x E D.
As above, instead of (7(x) we consider the functionf(n) of the normal n. The de-
pendence (7(n) upon the external field (7o(x) can be represented in the form
(4.9.32)
(4.9.33)
102 4. Local Defects
(4.9.34)
(4.9.35)
and consequently
(4.9.36)
In particular, for the cavity C 1 =- Co and S*(n) = So(n). Taking into account
(4.9.23) we have
It is interesting to note that this problem as well as the previous ones, were
solved for an arbitrary anisotropy without use of the explicit form for the
Green's tensor Go(n). It is sufficient to know its Fourier transform Go(k). In
other words, we find stress concentrations even though we have not been able
to solve the corresponding boundary value problem. This is a direct conse-
quence of the property of polynomial conservation.
The representation for F(n) obtained above is especially convenient for the
investigation of the dependence of the concentration coefficient on the shape
of the ellipsoid. Indeed, the local dependence of F(n) on the normal to the el-
lipsoid (or equivalently, on the point) is given by the first factor S*(n), which
is independent of the ellipsoid parameters and always remains finite. The second
factor <S*(n»-l does not depend on the point and becomes infinite in the case
of the passage to the limit of the elliptical crack.
It follows from (4.9.16) that the entire dependence of <S*(n» on the ellip-
soid parameters is concentrated in the scalar weight factor pen). This is essential
since it permits a complete investigation of the stress concentration in the limit-
ing cases of an ellipsoidal crack and a needle for an arbitrary anisotropic me-
dium by a series expansion of pen) for small parameter.
Note the interesting relationships
1
Dllll = To [1 - 81t (3111 + It)] ,
where
1 2po (4.9.41)
KO =
161tpo(1 - ]..10) , To = 1 - ]..10 '
the quantities
_ 3 2
_ 3
Ip - "2 V!
00 d~
(a~ + ~)Ll(~)' Ipq - "2 va p ! (a~ + ~)(a~ + ~)Ll(~)
00 d~
(4.9.43)
(4.9.44)
Let us consider two cases. First, let the external field O'ffP be the tension
along the coordinate axes, i.e. only the diagonal components 0'$" are nonzero.
Then the stress concentration is defined by the components
104 4. Local Defects
(4.9.45)
Secondly, let the external field off" be pure shear. Thd stress concentration
is determined by the components
(4.9.46)
(4.10.1)
where al ;::: a2 ;::: a3 are the ellipsoid's semiaxes. The case a <:: 1, ~ '" 1 cor-
responds to the needle, while ~ « 1, a '" 1 to the crack, and e « a « 1 to a
narrow crack. We note that in all cases we have e = ~a « 1.
It follows from (4.9.36), that the solution of the stress concentration prob-
lems for a needle and a crack reduces to the computation of the principal
terms in the expansion of the tensor D-l with respect to the corresponding
small parameters. From (4.9.23) and (4.9.16) we have
First we carry out the computation for a needle. Eq. (4.10.2) is transformed
to the spherical coordinates q;, () with the polar axis along the axis of the needle,
i.e. along Xl. We perform the change of variable cos () = t and we set
(4.10.3)
~ dtp
D - -2 J JD(q;, t)f1(t, al) dt ,
_ 21t 1
2 (4.10.4)
11: 0 cos q; + .;2·sm2 q; -1
where
(4.10.5)
For a needle a --+ 0, .; ,...., 1 andf1(t, al) --+ oCt). Consequently, the principal term
in the expansion of D with respect to a has the form
(4.10.8)
(4.10.9)
(4.10.10)
In all these cases the tensor DOl can be found in explicit form. For its
SHere and in the following we assume that the symmetry axes of the medium coincide
with the axes of the ellipsoid.
106 4. Local Defects
(4.10.11)
1
~(~A1VU1U2 + B1 ~ + C1vlu1D' (4.10.12)
(4.10.13)
Inserting DOl into (4.9.37) we can obtain, according to (4.9.32) the stress
q(n) at the surface of the needle for an arbitrary external field qo.
The expressions (4.10.12,14) give an obvious mechanism for the appearance
of the singularity in the conversion from a needle to a narrow crack, i.e. for
~ -+ 0: a singularity of order ~-1 appears only for external stresses qo con-
taining components with the index 3.
Now let us consider an arbitrary crack ~ « 1, a '" 1. In the same way
as in the case of a needle, we transform (4.10.2) to spherical coordinates 'P, 0,
but with the polar axis oriented along x 3• Setting
and making use of the definition (4.10.3) for D('P, t), we have
_ a 210 d'P
D - -2 J f
1
D('P, t)/2(t, ~1) dt (4.10.16)
'Jt 0 cos 'P + a sm 'P-l
2 2 . 2
with
(4.10.17)
In contrast to the case of the needle, for the computation of the principal
term in the expansion of D-l with respect to ~ we have to retain the first two
terms in D since the tensor Do in (4.10.17) does not have an inverse, i.e.
det Do = O. It is sufficient to show that all the components of the tensor Do
containing the index 3 are equal to zero. In fact, as ",ill be shown below, in the
integral for Do the tensor DC'P, 1), which corresponds to So(n) for nl = n2 = 0,
n3 = 1, occurs. But from the general expression (4.9.43) for So(n) it follows
that in the case of a rhombic structure, the corresponding components of
So(n) which contain the index 3 are equal to zero for nl = n2 = o. We should
expect that this is true also for an arbitrary anisotropy.
For the computation of Do and Db we consider12(t, ~l) as a generalized func-
tion on the segment It I :::;; 1 with identified points ± 1, i.e. on the circumference.
This is possible by virtue of the fact that D('P, t) is even and continuous. We
can verify that for ~l -. 0
D('P, ± 1)
D = ~21O
o !
2'Jt cos 2 'P + a 2 sin 2 'P d'P,
(4.10.18)
where the latter integral is written in regularized form. We note that the ex-
pressions obtained become considerably simpler for the circular (a = 1)
and for the narrow (a « 1) crack. In the isotropic case the integrals (4.10.18)
can be easily computed and the inversion of the tensor (4.10.17) allows us to
obtain explicit expressions for the components of the principal term of the
expansion of D-l
108 4. Local Defects
Here K(a) and E(a) are complete elliptic integrals of the first and second
kind, respectively. The remaining components «an be considered equal to
zero with the accuracy of 0(1). From this it follows that a contribution to
the singular stress O"(n) at the surface of the crack is given only by the com-
ponents of the external field 0"0 which have index equal to 3.
Passing to the narrow crack a --+ 0 yields
which coincides with (4.10.14) for ~ --+ O. In the case when the medium has
a rhombic symmetry, the tensor D-l has a similar structure and its components
can be expressed in the general case in terms of elliptic integrals. They become
considerably simpler for the narrow crack and have the form
(4.10.21)
which, obviously, in the particular case of the isotropic medium coincide with
(4.10.20).
We emphasize that, representing an independent interest, the case of the
narrow crack not only makes obvious the mechanism of the variation of the
stress concentration when passing from a needle to a crack, but also allows
us to pass to the plane problem if we consider the median section of the crack.
We now proceed to the direct investigation of the stress at the surface of
the crack. We start with a pure extension. Since the components 0"6 1 and 0"32,
according to (4.10.19), do not contribute to the singularity, we consider only
the extension 0"83 along the axis x 3 normal to the crack.
From (4.9.32, 37) and (4.10.19-21) we find
(4.10.22)
coordinates X"', connected at each point of the surface of the ellipsoidal crack
with the normal n. As a local basis we take
(4.10.23)
where ea is the unit vector of the coordinate system connected with the axes
of the ellipsoid, while the vector nO is normal to the boundary of the crack
at the point (nh n2, 0). It is obvious that for the points of the surface near
the edge, which are of primary interest, the axes XlI and x 2' are, respectively,
parallel and normal to the edge. With respect to the local axes, the tensor
O'alfJ'(n) is planar since all components with index 3 are equal to zero. This
follows from the equilibrium conditions.
First we consider the isotropic case. Transforming in (4.10.22) to local
coordinates and taking into account (4.9.43) and (4.10.19, 20), we obtain
O'1 /1'(n) = l.ioO'max(1 - n§), O'2 /2/(n) = O'max(1 - n§), O'1'2 /(n) = O. (4.10.24)
Here for the finite and the narrow crack we have, respectively,
tangential stress is 't"max = 10'21/2. It is clear from (4.10.24) that in the case of an
extension the largest value of the stress in any cross section normal to the edge
is attained for n3 = 0, i.e. on the edge itself.
In the anisotropic case the computations are similar but more cumbersome.
We give the final expressions for the stresses O'(n) at the surface of a narrow
crack in an orthotropic medium in the cross sections nl = 0 and n2 = O. In the
selected local basis, which in the present case coincides with the principle
axes of the tensor O'(n), we have for the cross section n2 = 0
Here (D-l)3333 is given by (4.10.21). The stresses in the cross section nl = 0 are
obtained by replacing in the right-hand sides the indices 1 and 5 by 2 and 4,
respectively.
110 4. Local Defects
We consider now a pure shear. The component (J6 2 does not contribute to the
singularity and the situation for (J6 3 and (J5 3 is similar. Therefore we assume that
the external field coincides with (J6 3.
In the isotropic case, we find from (4.9.32,37) and (4.10.19)
(4.10.28)
(4.10.29)
then the new axes are principal for (JaIW(n). For the principal stresses we find
(JI(n) .v
= (JI/I/(n) = (Jmax [(1 - ).io) (1 - nD (1 - n~) - (1 + vo) nln3] ,
.v
(J2(n) = (J2 /2/(n) = - (Jmax [(1 - vo) (1 - ni)(1 - n~) + (1 + vo) nln3] .
(4.10.30)
Here for the finite and the narrow crack we have, respectively,
(4.10.31)
The largest absolute value of the principal stresses 0"1 and (J2, equal to (Jmax,
occurs at the cross section n2 = 0 in the points
which are the maxima (minima) of the functions (JI(n) and (J2(n). In this cross
section (J1 and (J2 are equal to zero at the edge. It is essential that the maximum
point is situated at a very small distance from the edge of the crack, which
according to (4.9.14,30) is
(4.10.33)
We can also show that the characteristic width of the peak is of order c. By
varying the cross section from n2 = 0 to n1 = 0 the magnitude of the peak de-
creases while the point of the maximum approaches the edge coinciding with
it for n1 = O. For the cross sections in which (J1 and (J2 are of the same sign, the
character of the variation of '"max is evident. In particular, in the cross section
4.11 Crack in a Homogeneous Medium 111
n2 = 0 the quantity 't"max has a sharp peak of height O'rnax/2 at the point n* and
vanishes at the edge. In the case when 0'1 and 0'2 have different signs
This quantity reaches its maximum value equal to (1 - ).10) O'max at the edge in
the cross section n1 = O.
Thus, in the case of shear, a remarkable increase of stress takes place, which,
if not taken into account, wi11lead, in numerical computations, to qualitatively
irregular results on the stress concentration in the crack. This shows the neces-
sity of investigating the stresses not only at the characteristic points of the
contour of the crack, but also over its entire surface.
A similar phenomenon takes place in an anisotropic medium. The principal
stresses at the surface of a narrow crack in an orthotropic medium in the
cross section n2 = 0 are
Here a and (D-1)1313 are given by (4.10.27 and 21), respectively. In the cross
section n1 = 0
Obviously, as in the isotropic case, the increase effect takes place in the cross
section n2 = O.
From (4.11.1) one can find e+(x); then (4.11.2) determines the field O'(x) outside
the cavity.
Let us consider the limiting transition from the cavity to a crack, i.e. to a
cut along a smooth oriented surface D which is bounded by a closed contour
r which lies on the boundary of V. Let us choose a local coordinate system at
the point x' E D, such that its z-axis is directed along the normal n(x'). Let
h(x') be the transverse dimension of the cavity and Zl(X', h), Z2(X', h) be the co-
ordinates of the points of intersection of the z-axis with the boundary of V,
while Zh Z2 ~ 0, as h ~ 0. For a fixed point x ~ V, the kernel So(x - x') in
(4.11.2) is a smooth, bounded function. Hence,
Note that the quantity e+(x, h) can be interpreted as the coefficient of the
principal term of the expansion of the function e+(x) V(x) in multipoles, which
are concentrated on D,
(4.11. 7)
where
T(x, x') = n(x)So(x - x')n(x') . (4.11.10)
where r is the distance from the point x along the normal to a point Xo E r,
(3(xo) is a smooth function on r.
Let ut consider the asymptotics of the stress field O"(x) outside the crack in
the neighborhood of the edge r of the crack. Let Yh Y2, Y3 be the local cartesian
coordinates at the point Xo E r,
the Y2-axis being directed along the limiting nor-
mal to 0 at the point xo, and the Y3-axis being tangent to r.
Then the Yraxis
lies in the plane which is tangent to 0 at the point Xo. Taking into account
(4.11.12), we have
(4.11.13)
Using (4.11.8) and (4.11.3) we can write down the expression for 0" at the
point z = (- p cos (), - p sin (), 0) (p > 0)
Here ~ = P-1Yj and it is taken into account that So(x) is an even homogeneous
function of degree - 3. It is easy to show that, as p ~ 0, the integral tends to
a finite limit, and hence, the stress has a singularity of the order of p-ll2.
Let us introduce a tensor coefficient of stress intensity, which is of interest
for applications and defined by the expression
(4.11.16)
where
s(O)
00
o
00
-00
°
= - J -.I ~l d~l J So(cos + ~l> sin 0, ~3) d~3 (4.11.17)
that J~ooSo(~l> ~2' ~s> d~3 coincides, in essence, with the Green's tensor for stress
in the problem of plane deformation and complex shear (in dimensionless
variables ~) of a homogeneous medium with moduli Co, the normal to the
plane coinciding with the ~3-axis. Then J coincides with the stress at the point
~l = - cosO, ~2 = - sinO, when a discontinuity of the displacement vector
(J(xo) -.1ft is given along the positive semi-axis of ~l. In the isotropic case,
the function J(O) is easily found in an explicit form.
The tensor J can be represented as a sum of three tensors, which correspond
to the three components of the vector 1',
(4.11.19)
where
(4.11.20)
(do not sum over i!). The tensors J 1 and J 2 are found from the solution of
the corresponding planar problem, and J 3 is derived from the solution of an
anti-planar problem (complex shear). All these problems can be solved by
known methods for the case of an arbitrary anisotropy.
Note that in the case of linear fracture mechanics, the asymptotics of the
stress field in the neighborhood of the edge of a crack is always characterized
by the stress intensity factors K 1, K II , KIll. It can be shown that these
4.12 Elliptic Crack 115
(4.11.21)
(4.11.22)
where
(4.11.23)
This representation is convenient since s(O)n(xo) depends on a neighborbood
of the point Xo only, and the operator B(xo) depends on the whole surface Q.
Thus, the calculation of the tensor intensity coefficient J is finally reduced
to solving a certain planar problem and constructing the operator B. An
explicit expression for the latter, for the case of an elliptic crack, will be obtained
in the next section.
An elliptic crack is the simplest case which is of interest for applications and
admits an exact solution.
Let Q be a plane cut with an elliptic boundary in an elastic medium, and
b(x) be the discontinuity of the displacement vector at the cut. We have the
following analog of the theorem about polynomial conservation: for a poly-
nomial external field O'o(x) of degree m, the vector b(x) has the form
(4.12.1)
7This result was first obtained (in a very complicated way) in [4.4].
116 4. Local Defects
where X(Xlo X2) E D, n = const is the normal to D, eis a coordinate along the
normal and
(4.12.3)
(4.12.4)
where Q2k is a polynomial of a degree not higher than 2k. Hence it follows
that
(4.12.5)
where the integration is carried out over the whole plane, which contains D,
the vector b(x') is defined to be zero outside D, and T(x) is a generalized func-
tion, which is generated by the generalized function So(x) in three-space.
Substituting (4.12.1) for b(x) into (4.12.6) and taking into account (4.12.3),
we have
= J T(x)z(x)dx = 1
(2'Jt)Z J T(k)z(k) dk , (4.12.8)
where
(4.12.9)
sinlakl )
z(k h k z) = 2'Jta l aZ lakl- z ( lakl - cos lakl (4.12.10)
(4.12.11)
(4.12.12)
Problem 4.12.2. Verify that (4.12.1) for the vector b(x) has, in a neighborhood
of the point Xo E r,
the following asymptotics in the local coordinate system
Yh Yz (Sect. 4.11)
(4.12.14)
118 4. Local Defects
where
(4.12.15)
0)
a= (a1o az'
and B(xo) is the value of B(Xb xz) at the point Xo E r .
Thus, for a polynomial external field, the polynomial B(x) being considered
as a known, we have an explicit expression for {3(xo), and consequently, for the
tensor intensity coefficient
(4.12.16)
where
d 1 d 1 d-l
1 = - C3 + ).I(
Cz 2
- C3)' z= C3 + ).I(
C1 2
- C3)' 3 - C3'
C
1=
- Cz + 3C3' Cz = C3 - E(a) - K(a)
aZ - 1 E(a)
C _
,3 _ a Z'
and E(a) and K(a) are the complete elliptic integrals of the first and second
kinds.
(4.13.1)
Here v is the volume of the ellipsoid, C l is the tensor of perturbed elastic con-
stants, and A is a tensor which depends on parameters of the ellipsoid and on
the tensor of elastic constants Co of the external medium. In the general case,
the latter is given by the (4.9.21, 23), and for an isotropic medium this is given
by (4.9.39). The second equality in (4.13.1) is valid if CII exists.
For an ellipsoidal cavity, Cl = - Co and
(4.13.3)
(4.13.4)
.. (
PtiJ =
pao paoVVGO(Xl - X2)P&)
(4.13.5)
P50VVGO(X2 - Xl) pao P& '
where Xl> X2 are the position vectors of the centers of the ellipsoids. As is shown
in Sect. 4.5, the central symmetry of the inclusions enables us to ignore the
matrices PM u P{6.
Substituting PM in (4.5.13), we find an expansion of the interaction energy
<Pint in powers Of'12 = IXI - x21
120 4. Local Defects
(4.13.6)
where
(4.13.7)
(4.13.8)
It is obvious, that here </Jb is the self-energy of interaction of i-th defect with the
external field co, and </J 1 is a quadratic function of the external field, which
depends on Co and C .. (i = 1, 2), these depending on the defects and the unit
vector of the straight line which connects the defects.
We present below the explicit expressions for </J6 and </J 1 in the case, when
the external field is pure dilatation c~{l = cOoa{l and the inhomogeneities are two
spheroids with a common rotation axis x 2 and with isotropic elastic constants
2
"=
</Jo - ~
2Ko6...
(3P"1 + Pi)
2,
(4.13.9)
where
- 2[ - 12 + h2 - 8Te(1 - vo) ~. Vi ]2 .
fJ-i 1 + Vi
The quantities I p , Ipq contained in these expressions are given by (9.42), Ko
is expressed by (4.9.41), the index i is the number of a defect, and the lower
index of the quantities P:
is connected with the principal axes of the ellipsoids,
these taking the values 1, 2 only (due to the rotation symmetry). A limiting
case of a rigid inclusion is obtained as fJ-i ~ 00; the case of a cavity is obtained
as, fJ-i ~ - fJ-o, Vi ~ Vo·
If one of the ellipsoids, for example the second one, is a sphere, then P~ =
o and the expression for </J1 is greatly simplified
4.14 Notes 121
(4.13.11)
4.14 Notes
The content of Sects, 4.1-4.5 is based on [B6.l4, 20, 23]. As regards point
defects, see [B6.27, 39]. The Green's tensor for displacement in an anisotropic
medium was investigated in [B1.2] The Green's tensors for strain and stress
were introduced in [B6.l7, 24]. Elasticity in terms of projection operators was
formulated in [B6.25, 26]. Discontinuities of stress and strain on the boundary
of two media were obtained in [B6.24].
The polynomial conservation theorem for an ellipsoidal inclusion was
proved for an isotropic medium and homogeneous field in [B6.5-7] and for
an anisotropic medium and arbitrary polynomial field in [B6.22]. The com-
plete investigation of stress on an ellipsoidal inhomogeneity and its limiting
forms was carried out in [B6.2l, 24], see also [B6.l2, 31, 37, 42].
A survey on cracks can be found, for example, in [B6.8, 29]. As regards
integral equations for cracks, see [B6.43, 7.13, 8.6].
5. Internal Stress and Point Defects
In all the preceeding discussions, it was the displacement vector or the field of
the symmetric strain tensor of rank two, which is uniquely determined by
the former, that described a state of a medium of simple structure. Conversely,
if the strain field is given, satisfying the usual compatability conditions, then
the displacement field is determined to within a rigid translation and a rotation
of the medium. The deformation is connected with the stress tensor by Hooke's
law (local or nonlocal), the stress tensor being connected by means of the
equations of equilibrium, with the external body or boundary forces. In the
absence of external forces, both the strain and the stress are equal to zero.
According to such a scheme the classical as well as the nonlocal theory of
elasticity are constructed.
However, another situation is possible, occurring in a number of important
cases. In order to give an illustration, let us carry out a thought experiment.
Let us consider a medium which consists of particles connected by an
arbitrary system of elastic bonds of the type of springs. Let us also assume
that, in the initial state, the elastic energy of all the springs is equal to zero
and consequently forces of interaction between particles of the medium also
vanish. Let us take some of the springs, change their length (e.g. they are short-
ened) and return the springs to their original places. This, generally speaking,
will result in a deformation of the medium, and the latter will proceed to a
new state which corresponds to the minimum of the elastic energy. In this
state of equilibrium, as distinct from the initial state, the elastic energy of
bonds is no longer equal to zero. A convenient quantitative measure of inter-
action between particles of the medium is the stress tensor, which is not con-
nected at all with the external forces. To emphasize this fact, let us call this
tensor an internal stress tensor, while the usual stress, caused by external
forces, will be called the external stress.
5.1 Internal Stress in the Nonlocal Theory 123
(5.1.1)
Hence, it follows that the operator (/)la/lp is self-adjoint, i.e. its kernel satisfies
the following conditions
or in the k-representation
(5.1.4)
where ,,;,. and 'la are the external and internal distortions, respectively.
The external distortion ,,;,. is defined as a gradient of the displacement Ua ,
caused by external forces, i.e.
(5.1.5)
Thus, we have
124 5. Internal Stress and Point Defects
(5.1.6)
(5.1.7)
(5.1.8)
Hence, taking into account the arbitrariness of aJ.a and Xl'fl' we obtain a condi-
tion for C/JJ.al'fl [cf. 2.4.18)]
(5.1.9)
The energy operator C/JJ. al'f3 determines a scalar product in the functional
space of distortions, transforming the latter into a Hilbert space. This Hilbert
space can be decomposed into two orthogonal (with respect to energy) sub-
spaces of external and internal distortions, and the total energy can be re-
presented in the form of a sum of the external and internal energies.
We know from the results, obtained in the Chap. 2. that the most general
expression for the external elastic energy can be represented in the form
(5.1.10)
Here,
(5.1.11)
is the symmetric strain tensor and the Hermitian operator of the elastic moduli
cJ.al'fl is symmetric with respect to the indices Aa and pl3-
The external stress tensor (j' J.a is connected with the external strain c;a by
the nonlocal Hooke's law:
(5.1.12)
(5.1.13)
5.1 Internal Stress in the Nonlocal Theory 125
(5.1.15)
(5.1.16)
The quantity
will be called the tensor of internal stress. From (5.1.15), taking into account
(5.1.5) we find that (i.lrx satisfies the equation of equilibrium
(5.1.18)
(5.1.19)
lIn the case of a moving dislocation, the external forces perform work.
126 5. Internal Stress and Point Defects
(5.1.23)
where atfP{Jv is a constant tensor, which is anti-symmetric with respect to the
indices A)) and symmetric with respect to the indices p.fJ. In the presence of
a central symmetry and in an isotropic medium, the expansion of a.l.apfJ(k, k')
begins with terms of the second order in k, k'.
Here the following question can arise: why the tensor of internal stress
cannot be symmetrized with the help of transformations, analogous to those
carried out in Chap. 2. Recall, that there, besides conditions of invariance
with respect to translation and rotation, an essential use was made of the
representability of the external distortion in the form (5.l.5), which is ensured
by the condition (5.l.6). However this condition is obviously invalid in the
region, where the sources of internal stress are present.
It is assumed hereafter for the sake of simplicity that f/)AapfJ = cAapfJ. The
case aAapfJ =f. 0 can be considered in an analogous way but leads to more com-
plicated formulae.
In view of our assumption, (5.l.16) for the internal energy can now be
written in the form
(5.l.24)
which connects the symmetric tensor of internal stress with the symmetric
tensor of internal strain.
From the positive definiteness of the form (5. l. 24) , it follows that the
operator cAapfJ is non degenerate, i.e. there exists an inverse operator bpfJAa ,
which will be called the operator of elastic compliance. This enables us to
invert Hooke's law and write it down in the form
(5.l.26)
a = C e, e = Ba. (5.l.27)
Due to the assumption, concerning the equality of the operators (/J),ap.~ and
c),ap.f3, the role of the main kinematic variable, which determines the internal
energy (5.1.24), is now played not by the internal distortion (),a but by the
internal strain C),a' This enables us to simplify the geometric picture and to
confine ourselves to purely metric considerations.
Two metrics can be connected with the elastic medium. The external one
gives a distance ds' between points of the medium x a and x a + dx a in a de-
formed state, caused by the external forces. The corresponding expression
has the form
where g~{3(x) = gpa(x) is the external metric tensor, which determines the
external metric in the given state. Since we are working in the harmonic ap-
proximation, i.e. we are considering small deformations and small displace-
ments, the difference between Lagrangian and Eulerian systems of coordinates
is not essential. Therefore, the tensor g~f3(x) can be considered as a field quantity,
which is related to a fixed coordinate system x a in space, in which the deformed
medium is embedded. Without loss of generality, it can be assumed that the
coordinate system x a is cartesian, with the constant metirc tensor Oa~'
Let in the initial undeformed state the external metric tensor coincide
with oa{3' Then, as a measure of deformation we can take the difference
(5.2.2)
called the tensor of external strain. But the deformed state can be also described
in terms of the field of displacement Ua(x) , with which the tensor c~{3(x) is con-
nected by the relation (5.1.11). From the linear theory of elasticity it is well
known that, for the representability of c~ix) in the form of (5.1.11), it is neces-
sary and sufficient that the deformation c~{3(x) satisfies the Saint-Venant condi-
tions of compatibility. The latter can be written in the form
(5.2.3)
or briefly
Rot C = 0, (5.2.4)
2The operator Rot is sometimes denoted also by the symbol Ink. Attention should be
paid to the correct order of indices: when (5.2.5) is substituted in (5.2.3), the operators rot
are contracted with different indices.
5.2 Geometry of a Medium with Sources of Internal Stress 129
Problem 5.2.1. Verify that the operator Rot satisfies the identities
where the operators grad and def act on vector fields (recall that the operator
def is the symmetrized gradient).
Problem 5.2.2. Show that the operator Rot commutes with the operators of
symmeterization and alternation, and hence, conserves the symmetry of a
tensor of the second order.
Problem 5.2.3. Verify that analogously to the well-known identity for the
Laplacian
(5.2.11)
(5.2.12)
(5.2.13)
differs from zero only in the region, where there are sources of internal stress,
and it characterizes their density. This quantity is called an incompatibility.
From the properties of the operator Rot, it follows that 7]"p(x) is a symmetric
tensor, the divergence of which is equal to zero
(5.2.14)
Taking into account (5.1.18) and (5.1.27), let us write in abreviated notations,
the complete system of equations for the internal stress
Problem 5.2.4. Show that the incompatibility and the dislocation (or quasi-
dislocation) density are connected by the relation
(5.2.16)
the stress (J is given then r; is defined uniquely but a is defined with some
arbitrariness.
(5.2.18)
Thus, from the viewpoint of internal stress, the density of moments of dis-
locations is defined with even mOl e indeterminancy than the density of dis-
locations itself. In particular, the anti-symmetric component mCpl'l does not
affect r;AI' and, hence, (Ja{3 (but not al'». We consider this situation in more detail
in Section 6.1, with the example of a dislocation loop.
or briefly,
The tensor operators Z and G obviously, play here the role of the Green's oper-
ators for internal stress. From (5.2.18) it follows that they are connected by the
relation
G = ZRot. (5.3.3)
132 5. Internal Stress and Point Defects
where the Hermitian operator H is determined not uniquely. Later we shall use
this indeterminacy in order to obtain the simplest expression for H.
The internal elastic energy, which is defined by (5.1.24 or 28), can be now
represented in the forms
(5.3.8)
aIn the local limit, G coincides with the operator S of Chap. 4. Earlier, the symbol G
denoted the Green's operator for displacements, but the latter is not used in Chap. 5 and
6, which excludes a possibility of confusion.
4Here, as in Sect. 2.8, k and r denote the moduli of the vectors k and r.
5.3 Green's Tensors for Internal Stress 133
1 A(k) (5.3.9)
b(k) = 2fJ-(k) E1 - 2fJ-(k) [3A(k) + 2fJ-(k)] E2 .
Let us write the second equation of the system (5.3.4) in the k-representation
(5.3.10)
where
(5.3.11)
(5.3.12)
and have the required singularity, as k --* 0, which ensures the fulfillment of the
conditions for G(r) at infinity. Omitting calculations, we present the final result
(5.3.14)
(5.3.15)
(5.3.16)
134 5. Internal Stress and Point Defects
(5.3.18)
(5.3.19)
Thus, the problem is reduced to the determination of the inverse Fourier trans-
form of the functionf(k) = k- 4p(k)under the condition k::;; K. Calculations,
analogous to those carried out in Sect. 2.9 yield
p(b.) = Po P(b.) ,
(5.3.21)
(5.3.22)
(5.3.23)
is valid. Using the expansion for fer), it is easy to obtain the value of G(r) for
r = 0, which plays a significant role in the theory of point defects:
(5.3.26)
(5.3.27)
For the Debye model and for the approximation (5.18) we have, respec-
tively,
{3 = 1, (3 = 30 ~511 r . (5.3.28)
(5.3.29)
where flo f2 are dimensionless tensor functions of the unit vector na = xa/r
136 5. Internal Stress and Point Defects
and the Poisson coefficient ))0. For the Debye model, the second term in (5.3.30)
is absent and the asymptotics coincides with the zeroth approximation (with
a = 0). In the general case, the second term in (5.3.30) accounts for the devia-
tion of w(k) from the linear law and contains approximate information about
the presence of microstructure in the medium.
In conclusion, let us note that if the Green's tensor for an isotropic medium
is known, an explicit expression for the Green's tensor of a medium with weak
anisotropy can be constructed.
Let the operator of elastic compliance B have the form
B = Bo - BI . (5.3.31)
(5.3.32)
(5.4.1)
where bo(x - x') is the kernel of the operator Bo of the homogeneous medium,
bi is a constant tensor, which has the symmetry and dimension of the tensor of
elastic compliance.
5.4 Isolated Point Defect 137
Let us denote by Go(x - x') the Green's tensor of the homogeneous medium,
which corresponds to bo(x - x') and by G(x, x') the Green's tensor of the
medium with the defect, relating to b(x, x'). Omitting calculations, which are
analogous to those, carried out in Sect. 4.2, let us write down the expression
for G(x, x') in terms of Go(x - x')
The constant tensor P 11m has the symmetry of the tensor of elastic constants
and is the inverse of the tensor
(5.4.4)
where M)" is the effective moment of the defect, which depends on the dis-
persion law and the ratio of elastic constant of the defect and that of the medium.
Taking into account (5.4.4), after simple transformations one can obtain the
expression for M ;,,:
(5.4.6)
(5.4.8)
where
(5.4.9)
Let us consider the interaction between the defect and another field of stress.
We denote by -r(x) the unperturbed stress, i.e. the stress which would exist in the
medium in the absence of the defect. If -r(x) is caused by external forces, then
the defect behaves only as a local inhomogeneity. This case was considered in
138 5. Internal Stress and Point Defects
Sect. 4.2. Therefore, let us assume that -rex) is an internal stress, caused by some
distant sources with the moment density m (x). The resultant stress field will now
consist of the defect field O'(x), the field -rex) and a perturbed field rex). For the
latter, taking into account (4.3), it is easy to obtain
(5.4.10)
where
(5.4.11)
Thus, the perturbed field coincides with the field of a point source of internal
stress with the moment M.
For the energy of interaction Wint of the field -rex) with the field of the defect,
we have
(5.4.12)
where M~{3 is an effective moment of the defect. Hence, the force, which acts
on the defect, is equal to
(5.4.13)
Let us proceed to the case of an isotropic medium. The kernel of the operator
Bo is conveniently represented in the form 5
(5.4.14)
(5.4.15)
The operators Xo(Ll) and Po(Ll) are the analogs of the modulus of volume
compression Xo and the shear modulus PO' Assuming, as in Sect. 2.3, that the
Poisson coefficient ))0 is a constant, we have
PO(Ll) = Pop(Ll) ,
(5.4.16)
2Po(1 + )) 0) - (Ll) - (Ll)
Xo(Ll) = 3(1 _ 2))0) P = XoP ,
where P(Ll) is determined by the dispersion law. Recall (Sect. 2.3) that, for the
SAs above, in the case of an isotropic medium we introduce the notations rand k for the
position-vectors of points and r, k for their moduli.
5.4 Isolated Point Defect 139
Debye model, P(6.) = 1 and, for the model with a polynomial approximation
of the dispersion curve, the operator p(6.) is given by (5.3.21).
Set now B = Bo - B I , where BI is the operator of perturbed elastic com-
pliance. For better physical clarity it is expedient to express the kernel of the
operator BI in terms of the kernel of the operator of perturbed elastic moduli
C I . Confining ourselves to the case of a spherically symmetric defect, we obtain
P = PIll + pzlz,
-1 _ 2po(1+ ))0) (2~(1 - 2))0) _ Xo + Xl )
PI - v(1 - 2))0) 3(1 - ))0) Xl'
(5.4.18)
P21 = 2po ( ~(7 - 5))0) _ Po + PI) .
v 15(1 - ))0) PI
The formulae (5.4.3), (5.3.24) and (5.4.18) give an explicit expression for
the Green's tensor G(r, r') of an isotropic medium with a point defect.
Let the defect be a centre of dilatation, i.e. MAl' = M AI" Then, for the a
stress caused by the defect, we find from (4.15) (for convenience of notation
we assume r 0 = 0)
(5.4.21)
(5.4.22)
140 5. Internal Stress and Point Defects
(5.4.23)
However, at the same time, a vacancy is a source of internal stress, its asymp-
totics being of the form (5.4.21), and, as r -+ 0, the condition 0'(0) = 0 should
also be satisfied. The corresponding model in an elastic medium is a hollow
sphere, which is inserted with stress into a cavity.
It is easy to show that all the above conditions can be satisfied, if one
assumes that to the vacancy there corresponds the moment density
(5.4.25)
Then the stress caused by the vacancy takes the form (for the Debye model)
(5.4.26)
where Tr{ z-} = z-~, and the moment M' is determined by (5.4.20).
Let us consider a system of point sources of internal stress with the moment
density
(5.5.1)
In the general case, with the sources of internal stress there is associated a
local change of elastic constants, and hence, for the kernel of the operator
B we should set
The matrix pjj with tensor components is the inverse of the matrix
(5.5.4)
(5.5.6)
(5.5.7)
The first term 1/)0, as it should be, is equal to the sum of the self-energies of
each defect, which are given by (5.4.27). For the coefficients 1/)1 and 1/)2, we
have
(5.5.8)
where
mj=
(5.5.9)
, fl.j , Xj
fl.j= Xj= --'--'-'--
fl.o + fl.j , Xo + Xj
In the paper by Lifshitz and Tanatarov [5.11], the first term of the
expansion of the elastic energy of interaction for two identical centers of
dilatation in an ordinary isotropic elastic medium was calculated. Comparison
shows that it coincides with the term f/JIp-6 shown above for the case of identical
defects in the Debye model.
5.6 Notes
We suppose that the reader is familiar with the basic physical notion for dis-
locations as elementary carriers of plastic deformation in crystals. At the pre-
sent, elements of the theory of dislocations are standard tools for mechanical
engineers and stress analysts. Descriptions of the theory of dislocations, which
are intended for different groups of readers, can be found, for example, in
[6.1,2]. In this section we present some facts from the continuum theory of dislo-
cations, which are necessary for what follows.
Recall that a dislocation in a crystal is a linear defect of the crystal lattice
of a special type. In Fig. 6.1 an edge dislocation is schematically shown, its line
B/
shown, this being formed by the displacements of the surfaces of the cut ABCD
in the direction AB. Here, the segment DC belongs to the line of the dislocation.
The general case of a curvilinear dislocation can be considered as some sort of
combination of edge and screw components.
In Fig. 6.3, the edge dislocation L is surrounded by a closed contour r which
is called the Burgers circuit. Its image F' in a crystal without dislocation has a
discontinuity. The vector b required to close the contour, is called the Burgers
vector. For an edge dislocation the vector b is perpendicular to the dislocation
line L, while for a screw dislocation it is parallel to the latter.
12 11 10 8 8 7 1211 109 8 7
1.1 6 1.1 B
-c-" -
j..-- ~-
1"15 [j lq 5
18 .1
1718 19 1 2
I-- r
r'!-- f- -
15
15 b 0"
17 18 19 r- 1 Z
Fig. 6.3.
Burgers circuits
In the general case, the Burgers vector has an arbitrary angle with respect
to the dislocation line, this angle being able to change along the line for a
linear dislocation, but the Burgers vector remains constant. It is the existence
of a Burgers circuit with a non-zero constant vector b that is the definition of a
dislocation in a crystal. From the definition, it follows that the dislocation line
is always either closed or ends at the edge of the crystal (or goes to infinity
in the case of an unbounded crystal).
The plane, which passes through the Burgers vector and is tangent to the
line of dislocation, is called the slip plane. Dislocations move mostly along these
planes.
In the continuum theory of dislocations a real crystal is replaced by a con-
tinuous elastic medium and the dislocation is modelled by a singular line
L, which possesses the property that the integral
(6.1.1)
has a constant value when taken over any closed contour r enclosing L. Here
up. are displacements in the medium, and rand b are called the Burgers circuit
and the Burgers vector, respectively.
Besides an individual dislocation, there are also dislocations, which are
continuously distributed over a surface or a region. In such a case, the dis-
146 6. Dislocations
Analogously, for a surface S and region Vwe have O(S) and o(V) according to
So(S)g>(x) dx = Sg>(xs) dS ,
s
(6.1.3)
So(V)g>(x) dx = Sg>(xv) dV .
v
It is clear that o( V) coincides with the characteristic function of the region. The
above-introduced notation is due to the convenience in writing down the
following formulae.
Concerning L, S, and V it is assumed only that the corresponding integrals
in (6.1.2) and (6.1.3) exist.
The following representations are valid
Sdx g>(x) So(x - XL) dL = SdL So(x - xL)g>(x)dx = Sg>(XL) dL. (6.1.5)
L L L
where dL and d8 are the elements of the curve and the surface.
From Gauss' and Stockes' formulae, important relations for the vector 0-
functions follow.
Let S be the boundary of a region V. Then
In order to prove, for example, (6.1.8), let us write down the relations
(6.1.12)
where we have used the definitions of o-functions, the usual rules of operation
with derivatives of generalized functions and the analog of the Gauss theorem
(recall that we are working in an affine coordinate system). Then (6.1.12)
immediately implies (6.1.8).
(6.1.13)
148 6. Dislocations
or
(6.1.14)
Analogously, other a-functions with weight can also be defined. The deri-
vatives of a-functions with weight are defined, as usually, by a transposition of
the operation to the test function. In the cases when no confusion can arise,
we shall omit the brackets in the expressions for the a-functions with weight
and write, for example, rotl'.l.M.l.(xL)a(L), where the operator rot is applied to
the entire a-function with weight.
It can be easily seen that the following relations are valid:
(6.1.15)
where [.l. and n.l. are the unit tangent vector and the normal, respectively.
Note that formulae of the type (6.1.8) and (6.1.10) are applicable only to
a-functions with weight of the special form (6.1.15). More detailed information
on a-functions concentrated on lines and surfaces, may be found in [6.4, 5] .
Let us now proceed to the consideration of characteristics of a distribution
of dislocations. Let L be a closed (or going to infinity) line, which is the bound-
ary of the surface S, and let bv be a constant vector. The expression
(6.1.16)
(6.1.17)
A comparison with (5.2.17) shows that this dislocation density can be as-
sociated with the density of moments
(6.1.20)
where
(6.1.21)
mp"(x) = Mp"(xv)o(V) ,
al'"(x) = rot~p Mp"(xv)o(V) , (6.1.22)
7}AI'(X) = Rot~~12 Mp"(xv)o(V) ,
where Mp"(xv) is the volume density of moments of dislociations. If the dis-
locations are distributed over a surface S or a contour L, then o(V) is to be
replaced by o(S) or o(L), and Mp"(xv) by a surface moment density Mp"(xs)
or by a linear moment density Mp"(XL), respectively.
In the limiting case of an elementary dislocation loop, which can be obtained,
for example, by contracting the surface Sin (6.1.19)to a point Xo or, equivalent-
150 6. Dislocations
ly, by moving away far enough the point of observation x, the expression for the
density of moments takes the form
(6.1.23)
The computation of static fields and the interaction of dislocations within the
framework of the simplest model of an isotropic elastic continuum has been
dealt with in many papers, their results are given in the monographs and reviews
mentioned above. Computations, for specific crystals, carried out on com-
puters, showed that, in many cases, formulae of the continuum theory of dis-
locations remain valid even in the field where they seem to be inapplicable,
namely, for distances of the order of a few interatomic distances. However, at
the same time, it was also found that a number of effects essentially depend
on the discrete structure of the crystal, and for their theoretical description one
has to refute the continuum model.
An intermediate place between the crystal lattice and the elastic continuum
is occupied by the above investigated Debye's model of the quasi-continuum
with the space dispersion taken into account or not. As we know, this model
gives a qualitatively correct description of the discrete structure and it is not
too complicated for theoretical analysis. In Sections 6.4-6 within the framework
of this model, a few problems will be considered, in which typical effects of
microstructure appear. In particular, taking into account a change of force
constants in the core of the dislocation is one of such problems.
Let the operator of elastic compliance B for an isolated dislocation with
the line L, has the form
(6.2.1)
where B'(XL) is a tensor, which has dimension and symmetry of the tensor of
elastic compliance. It is obvious that B'(XL) must depend on the direction of
the Burgers vector and on the tangent vector to the dislocation line.
The construction of the Green's tensor G, which corresponds to (6.2.1),
the Green's tensor Go of the homogeneous medium being given, leads, in the
general case, to the solution of an integral equation on the dislocation line L.
If the perturbation of the force constants in the core' of the dislocation is small,
then the integral equation can be. effectively solved by the method of successive
approximations. In the first appoximation, we have
(6.2.3)
(6.2.4)
H(x, x')
= Ho(x - x') + J
s Ho(x - xL)RotB'(XL)RotHo(XL - x') dL .
L
(6.2.5)
(6.2.6)
where the tensor qiap(L) depends only on the line of the dislocation
(6.2.8)
152 6. Dislocations
Let us now consider the interaction of the dislocation with a stress field
o-(x) , which is caused by sufficiently remote sources of stress. The total force
q, which acts on the dislocation, can be represented as the sum of two com-
ponents
q = qo + q', (6.2.9)
where qo depends linearly on the Burgers vector b of the dislocation and on the
external field (J, and q' is quadratic with respect to the external field and does not
depend on b. Omitting not difficult but cumbersome calculations, we present
the final result.
The force qo can be written in a form, which coincides with the known ex-
pression of Peach-Koeller [6.1]
(6.2.10)
but with effective stress "C1'"' In the case of a small change of the force constants
in the core of the dislocation,
(6.2.12)
Beginning with this section, we consider the straight line dislocations, which
are parallel to the z-axis. In this case the stress and the strain depend only on
the two-dimensional position-vector r(x, y).
Let us decompose the 3-dimensional stress tensor (Jjk (i, k = 1, 2, 3) into a
two-dimensional tensor (Ja,B, vector aa and scalar (J
(Ja,B 1 (Ja )
(Jik = ( ..........:.......... , a, f3 = 1, 2. (6.3.1)
(Ja i (J
6.3 Two-Dimensional Problems 153
(6.3.2)
where lfJ, ¢ are scalars, and the operators Rot and rot are defined by the ex-
pressions
(6.3.3)
(6.3.4)
From the last equation, taking into account the boundary conditions at infinity,
it follows that c = 0.
The operator form of the Hooke's law can be now written in the form
Here
If the operator ($af>). is equal to zero, the system (6.3.6) is uncoupled. This
will take place in all cases, when, due to the symmetry conditions in the plane
x, y, there do not exist two-dimensional tensorial material characteristics of a
medium of an odd order. In particular, these conditions are satisfied in the
presence of central symmetry, cubic symmetry or in the isotropic case.
Let us represent the solution of the system (6.3.6) in the form
(6.3.9)
When the system (6.3.6) is uncoupled, the Green's operator is diagonal and
hence, according to (6.3.9), in the case of a homogeneous medium, edge dis-
locations do not interact with the screw ones. If the change of force constants
in the dislocation core is taken into account, then the above degeneration is
removed, since the diagonal components of the Green's operator depend on
the changes of force constants of all dislocations'!
Let us now proceed to a separate consideration of screw and edge dis-
locations.
Suppose there are only screw dislocations in the medium, i.e. ba = O. Let us
further assume that the conditions of symmetry are satisfied, ensuring the
uncoupling of the system (6.3.6). Then the elastic field of screw dislocations is
described by a scalar function cj;(r), which satisfies the equation
(6.4.1)
In the case under consideration, only the components (Ta(r) of the stress differ
from zero, and for them, according to (6.3.2), we have
Thus, the problem of determining elastic fields and of investigating the inter-
action between screw dislocations, is reduced to constructing the Green's
function G. In this section we consider the case of a homogeneous isotropic
medium without taking into account the change in the force constants in the
core of the dislocation. The kernel of the operator ~aJ3 and the Green's function
G are then of the difference type.
For an isotropic medium, taking into account (5.3.9), we obtain
(6.4.5)
where p is an operator with the kernel p(k). (6.4.2) takes the form
From the above discussion it follows that G(r) is the inverse Fourier trans-
form of p(k)k- 2, where k = Ikl. It is not difficult to prove that, in the two-
dimensional case, the inverse Fourier transform of the function k- 2 (k :::;; K) is
g(r) = -.l.
2rc
(7 Jo(,") -
0 '"
1 d," + c), (6.4.7)
where J o(,") is the Bessel function of the zeroth order, r = Irl. The constant c
should be chosen from the condition that, at large distances (Kr ~ 1), the func-
tion g(r) is to transform into the classical Green's function of the Laplace
operator
1 r
g(r) ~ - 2rc ln R' (6.4.8)
body. Taking into account the asymptotic behaviour of (6.4.7), we find that
c = In "R. The solution of (6.4.6) may now be represented in the form
(6.4.10)
where
(6.4.11)
and J 1("r) is the Bessel function of the first order. When" ..... 00 (a ..... 0), the
function air), becomes the ordinary two-dimensional a-function.
For the Debye's model
G(r) ~ - i~ In ~ , (6.4.13)
i.e., in the case of screw dislocations the space dispersion makes no additional
contribution to the asymptotics of G(r).
For a single screw dislocation at the point r = 0 with the Burgers vector b,
taking into account (6.4.3) we have
(6.4.15)
(6.4.16)
6.4 Screw Dislocations 157
wl;tere the constant G(O) for the model (5.3.18) and for the Debye's model is
equal to
respectively.
Consider now a system of screw dislocations with the Burgers vectors bj'
the dislocations being located at points r j' Then, for the stress we have
(6.4.18)
(6.4.19)
Thus, in the case of two dislocations, we find for the energy of interaction
(6.4.20)
(6.4.21)
1'\
\ 1'=-7
~
~ f5
~ !/ /'
-7~'" .L
=--
5
~
o1 pia 7
/ J
-OJ!
\) Fig. 6.4. Interaction force q as
function of distance p
158 6. Dislocations
(6.4.22)
respectively.
The asymptotics of (6.4.21, 22) coincide with the expression for the force of
interaction, which is known in the continuum theory of dislocations.
Since /JO(Kp)/ :s; 1, it follows from (6.4.22) that, in the Debye's model, the
sign of the force cannot change with decreasing distance between the dislocations.
At the same time, an investigation of (6.4.21) shows that, when the space dis-
persion is taken into account, the sign of the force can change, provided the
boundary value W(K) of the dispersion curve for transverse oscillations is greater
than the corresponding Debye frequency. The graphs of q(p) for different values
of r are presented in the Fig. 6.4. The dotted curve corresponds to the Debye
model.
where Bo(r, r') is the kernel of the operator of the homogeneous medium, the
B/s are constants which have dimensions of POl. For convenience of inter-
pretation of the constants B i , let us consider single screw dislocation within the
framework of the Debye model. The kernel of the operator p in this case can
be represented in the form
(6.5.2)
(6.5.3)
where
(6.5.4)
6.5 Influence of Change of the Force Constants in Cores of Screw Dislocations 159
Let us denote by G (r - r') the Green's function for the homogeneous med-
ium, which was constructed in the preceeding section. Then, for the Green's
function G(r, r') for the operator (6.5.1), we can derive the expression
(6.5.9)
Taking into account that oaGO(O) = 0, we find from (6.5.6, 8 and 9) that the
change of the force constants in the core of an isolated screw dislocation does
not contribute to the stress and the elastic self-energy.
Let us consider two screw dislocations with the Burgers vectors b l = b,
b2 = ± b. Assume also that BI = B2 = B. The elastic energy in this case has the
form
(6.5.10)
a[Go(p)]2
G(rb rl) = G(r2' rz) = Go(O) + a2 _ [G;;(p)]2 ,
(6.5.11)
_ [Go(p)]2G;;(p)
G(rb r2) - Go(p) + a2-[G;;(p)]2 '
where
160 6. Dislocations
The value Go(O) is given by (6.4.17), and the second term in (6.5.12), for the
Debye's model and the model (4.3.18), respectively, has the form
f - fJ.o,,2
[ 1 G'()]
poP p=o
_
-1- 81t'
(1 _2.8 T) .
fJ.o,,2
-9"1t
(6.5.13)
(6.5.14)
where (/)0 is the self-energy of the dislocations, (/)int is the interaction energy.
For the latter, taking into account (6.5.10,11), we obtain
d (6.5.16)
q(p) = - dp (/)int(P) .
q(p) ~ fJ.o b2
27rp
(± 1 + ~).
7rap2
(6.5.17)
From this expression, it follows that the change of the sign of the force of inter-
action at large distances, is possible only under simultaneous fulfillment of
two conditions: the values of the force constants in the core of the dislocation
increase, and the limiting frequency we,,) is higher than the Debye's frequency.
An analysis of (6.5.16) in the general case, shows that change of the sign of
the force at small distances (of the order offew interatomic distances) is possible
only if the spatial dispersion is taken into account.
Suppose that there are only edge dislocations in the medium, their lines being
parallel to the z-axis and their distribution being given by a density of the
vector baCr). As in the case of screw dislocations, we assume the fulfillment
6.6 Edge Dislocations 161
(6.6.1)
In this case, for non-zero components of stress, taking into account (6.3.2,
5), we have
(6.6.4)
where
(6.6.5)
It is not difficult to prove that the function her), which is the inverse Fourier
transform of k- 4, has the form
162 6. Dislocations
(6.6.9)
2
G(r) = -1--- ph(r) . (6.6.10)
- ))0
For the Debye's model and for the model (5.3.18) we have
2po
G(r) = T=-))o h(r) ,
(6.6.11)
G(r) = 1 2!0))0 [her) - p;;-2g(r) + 2r;- 1 ,,-40/ r )J,
respectively.
In the case of a single edge dislocation, the substitution of (6.6.11) into
(6.6.3) enables us to obtain closed expressions for the stress. For the Debye
model they coincide with the corresponding expressions obtained by Brailsford
[6.6].
Consider now two parallel dislocations at the points rl and r2 with the
Burgers vectors bI = (b, 0), b2 = (± b, 0) (the slip planes are parallel). Then,
for the energy of interaction I/>int' according to (6.6.6), we have
(6.6.12)
The projection of the force of interaction onto the slip plane for the model
(5.3.18) is equal to
(6.6.13)
For the Debye's model the last two terms in the above expression should be
omitted.
In the polar coordinate system the asymptotics of qx has the form
The first term coincides with the known expression for qx in the continuum
theory of dislocations. As distinct from the case of a screw dislocation, the
dispersion makes here the additional contribution of the order of p-3 to the
asymptotics. For the Debye model the asymptotics is given by the first term
only.
6.7 Notes 163
(6.6.15)
(6.6.17)
where
(6.6.18)
6.7 Notes
In this chapter the method of the effective field is applied to solve problems for
composites and cracked solids. Under the assumption of a random change of
the effective field from one particle to another the formulae for the first and
second moments of random stress-strain fields are presented.
7.1 Background
field for every particle, which is assumed equal to the external field, applies to
the whole body. Let the solution of the last problem be known and therefore
the function describing the state of every isolated particle can be expressed in
terms of its properties, the effective parameters of the medium and an external
field. After the summing of perturbations from all the particles an overall stress
and strain may be found and an effective elastic moduli tensor of the composite
is expressed in this way through the particle states functions.
As a result this way leads to a system of algebraic equations for the effective
parameters of elasticity (Sect. 7.3). Such a modification of SCM will be called
the method of an effective medium.
For the case of a composite medium another avenue for solving this problem
can be suggested. Every particle may be considered isolated in a homogeneous
medium-the matrix of the composite. The presence of surrounding particles is
accounted for through an effective external field applying to every particle. The
scheme which is based on such an assumption will be called the method of the
effective field.
It is commonly accepted that the effective field is constant (and the same
for all the particles). In this case the method of the effective field coincides
with a modification of SCM, which is widely used in quantum mechanics,
statistical physics and also in the mechanics of heterogeneous medium. The
assumption about the constant effective field can be confirmed by calculating
the effective parameters of an inhomogeneous medium. But in order to describe
the fluctuations of elastic fields in a composite we have to account for the
fluctuations of the effective field.
In this chapter the method of the effective field is applied to solve problems
for composites and cracked solids. Under the assumption of a random
change of the effective field from one particle to another the formulae for the
first and second moments of random stress-strain fields are presented. The
general scheme of the method is considered in Sect. 7.2-5. Sections 7.6-9 are
devoted to the first approximation of the method, coinciding with a modifica-
tion of SCM. Formulae for the effective constants of composites are obtained.
The construction of the second moment is described in Sect. 7.11 for the case
of point defects. The results are compared with the exact solutions for regular
composites and experiments. Estimates for the accuracy of the method is
given. Finally, the use of the results obtained for composite fracture is
shown in Sect. 7.12.
Let us consider an infinite elastic medium the elastic moduli tensor of which
is
(7.2.2)
where V,(x) is the characteristic function of the region occupied by the i-th
inclusion, the tensor Cli is constant within the i-th inclusion but, generally
speaking, varies from one inclusion to another.
The elastic compliance tensor of the inhomogeneous medium B = C-l
is represented in the same way as (7.2.1, 2),
where Bo = COl, Bli are the constant tensors within each inclusion.
If the fields of stress o-t(x) and strain ct(x) inside the inclusions are known,
the tensors I1(X) and c(x) at any point of the medium are represented in the
form (Sect. 4.8)
where co and 110 are the external strain and stress fields, respectively. Here
we denote
(7.2.5)
B l l1+(x) = 1:; Blil1t(X) V,{x) .
i
The kernels of the integral operators Ko and So in (7.2.4) are defined via
the Green's tensor for displacements of the homogeneous medium (Sect. 4.6).
It will be assumed further that the regions Vi occupied by inclusions consti-
tute a homogeneous random field and all the tensors Cli are independent ran-
dom variables having the same known density functions. Such a model involves
the most interesting composite materials. If an external field 110 (or co) is constant,
the functions C1c+(x) and B l l1+(x) for each realization of the random distribu-
tion of the regions Vi in space may be represented at the sum of an almost per-
iodic function and a function with bounded support. The operators Ko and So
are defined upon such a class of functions in Sect. 4.8 and Appendix A3. It is
essential that a unique definition of these operators is possible, if it is known
what kind of external field is given: the stress field 110 or the strain field co.
7.2 Formulation of the Problem 169
Further, for definiteness we assume that the external stress field 0'0 is always
fixed. Then the operators Ko and So are defined on constants via (Sect. 4.8)
The problem for a fixed external strain field can be treated similarly.
Let us consider separately the important case of a homogeneous medium
containing isolated cavities or cracks. For cavities C li = - Co, and Bli =
(Co + Cli)-l - COl tends to infinity while C li ~ - Co. Since
(7.2.7)
Note that the right-hand sides of (7.2.4 and 8) may be interpreted as the stress
and strain fields in a homogeneous medium, containing, in the regions V k ,
sources of dislocations moments with the density
(7.2.9)
(7.2.10)
in the case of cavities (cf. Sect. 5.3). It is essential that these dislocation moments
are induced in the medium by applying the external field 0'0 because the fields
O'+(x) and e+(X) are linear functions of 0'0'
Let us consider now the case in which defects are cracks or cuts along the
smooth, oriented surfaces Ok' If we pass to the limit from a cavity Vk to the
crack Ok, as it was shown in Sect. 4.11, we have et(X) ~ Mk(x)o(D k)
where
(7.2.11)
nk(x) is the normal to the surface Ok, b,.(x) is a vector field at Ok which can be
interpreted as a jump of displacements across the crack, 0(0.) is a delta-func-
tion concentrated at Ok (Sect. 6.1).
170 7. Elastic Medium with Random Fields of Inhomogeneities
As it follows from (7.2.6), the stress and strain fields in a cracked medium
are represented in the form
where we denote
(7.2.13)
The equations for the vectors bix) follow from the boundary conditions
on the surfaces of the cut and have the form
For the case of a single crack a similar equation is considered in Sect. 4.11.
X E Vj , (7.3.1)
where m,,(x) is defined by (7.2.9 or 10) for the case of inclusions or cavities and
by (7.2.11) for the case of cracks. The field O';(x) defined by this relation in the
region V j may be interpreted as the external field of the i-th defect. In the field
O',{x) this defect behaves as an isolated one.
Sometimes it is more convenient to consider an external strain field e,{x)
of the i-th defect. The connection between cSj(x) and O',{x) is represented by the
obvious relation
Let the solution of the elastic problem for isolated defects in an arbitrary
external field be known. This means that we have explicit expressions for the
7.3 The Effective Field 171
functions mk(x) in (7.3.1) henceforth denoted by mk(x, if k). Thus the system
of equations governing the fields O'k(X) for each of the interacting defects,
is obtained from (7.3.1) in the form
(7.3.3)
If the solution of this system has been found, the stress and strain fields in an
inhomogeneous medium are determined by
(7.3.4)
which are the consequences of (7.2.4, 8). Thus the fields ifk(X) can be considered
as the main unknowns of the problem.
If inhomogeneities constitute a random field, ifk(x) are random functions.
For the purpose of constructing the statistical moments of ifk(x) let us introduce
the following simplifying assumptions about the structure of O'k(X) (hypothesis
of the effective field method) :
HI) The field O'k(X) is practically constant inside the region Vk (or Ok) but,
generally speaking, these fields are different for the different defects.
H2) The random variables 0' k do not depend statistically on the geometric
characteristics and elastic constants of the defect Vk •
Note that the constants if k may be regarded as the average of the external
field O'k(X) over the volume Vk.
The field O'k(X) defined in the region V = UV k and equal to the constant
O'k in Vk will be called the effective field. Note that in the classical scheme of the
self-consisted method if is assumed constant and equal for each particle [7.3, 4].
In the present modification of the method this field is considered as a random
one and as it will be shown below, it enables us to describe more precisely the
interaction between defects in the various stochastic fields of inhomogeneities
The picture for the interactions corresponding to the hypotheses HI and H2
may be described qualitatively in the following way. For the typical defect the
external field O',.(x) which is represented as a sum of the external field 0"0 and the
fields caused by all the surrounding inhomogeneities, is approximately constant
in the region V .. occupied by this defect. Furthermore the contribution of every
individual defect to this field (0' .. ) is negligible. Without considering the realm
of application of these hypotheses, we proceed now to an analysis of the formal
consequences.
Further we shall assume that defects are inclusions of ellipsoidal shape, or
plane elliptical cuts (cracks). As follows from the solution for an isolated el-
lipsoidal inhomogeneity in a constant external field, the density of the disloca-
tion moment mk(x, O'k) induced in the region V k is represented in the form
172 7. Elastic Medium with Random Fields of Inhomogeneities
(7.3.5)
where
(7.3.6)
(7.3.7)
The linear transformation a,/ transforms V k into the unit sphere, Ko(k) is the
k-representation of the Green's function for strains, ill is the surface of the unit
sphere in k-space.
In the case of elliptical cracks the density mix, iJ k ) in (7.3.3, 4) has the form
(7.3.8)
where the expression for P k(X) is found from the solution of the problem for an
elliptical crack in the constant external field iJ k (Sect. 4.12) and has the form
(7.3.9)
( l 2)
h kX,X _ !Ii
- b
k
J (XI)2
1- -
ak
- (- X2 )2 ,
bk
(7.3.10)
where the axes Xl> X2 are connected with the main axes of the crack. Here ak and
bk are the half axes of the ellipse Ok, ak 2 b k. The constant tensor P k in (7.3.9)
is defined by
(7.3.11)
where nk is the unit normal to Ok and the tensor To has the form of(4.12.8).
In the case of an isotropic medium TOI is
(7.3.13)
Let us fix now a point Xo and introduce the region Vxo by the relation
if Xo E V
(7.3.14)
where Vj is the region occupied by j-th inclusion. Vxo denotes the complement
of Vxo to V. Evidently Vxo is the region Vj in which the point Xo is situated.
The characteristic functions (of argument x) of the regions Vxo and Vxo are
denoted by V(xo; x) and V(xo; x), respectively: hence
Let P(x) be an arbitrary continuous tensor field which coincides with Pi'
see (7.3.6), in the regions Vi. Let the field O'(x) be defined in the region V by
If the hypothesis HI is valid, the field O'(x) coincides with O'i in the regions
Vi and consequently with the effective field in V.
For the case of cracks let us introduce the region 0 of all surfaces of the
cracks and the region Oxo which is
(7.3.18)
Thus the equation governing the effective field in the region 0 is represented
in the form
where P(x) is an arbitrary continuous tensor field coinciding with P k(X) , see
(7.3.9), on the surfaces Ok.
Equations (7.3.17) and (7.3.19) represent the grounds for constructing
the statistical moments of the effective field 0'. It will be necessary for the
solution of this problem to calculate the various averages of the random
fields which have to be considered here. In the following section we shall
discuss in details the problem of construction of these averages.
174 7. Elastic Medium with Random Fields of Inhomogeneities
and the higher moments of f(x) are defined in the same manner.
In what follows, an important role will be played by the means of f(x)
calculated under the condition that the region V contains a fixed point Xl' For
this conditional mean, we introduce the commonly accepted notation 2 (f(x)lxI)'
By definition [7.5],
Analogously, the mean of the field f(x) under the condition that the region
V contains the points Xh ... , Xn has the form
n n
(f(x)lxh X2, ... , x n ) = (IT V(X;)-l S/(XI) IT V(x,-}df.l(V)' (7.4.4)
;=1 ;=1
Let now X E V, Xl E V x ' The mean of the fieldf(x) under these conditions
will be denoted by (f(x)Ix ; Xl), and the expression for this mean has the
form
2The symbol <I> is nowhere used in this chapter to denote the scalar product, therefore
no confusion should arise.
7.4 Several Mean Values of Homogeneous Random Fields 175
(7.4.6)
is valid.
In the general case, the sign";" separates variables which belong to the
region V x , from variables which belong to the region Vx • Thus the mean of
!(XI)!(X2) under the condition Xl E V, X2 E V:r' E VX1 is written in the form
(7.4.9)
and so on.
Let us note the equality which is a consequence of (7.3.16)
(7.4.10)
(7.4.11)
(7.4.12)
176 7. Elastic Medium with Random Fields of Inhomogeneities
where on the left-hand sides we have the means with respect to the ensemble
of realizations of the random function V(x), G is a region of the volume v, which
occupies the whole space at the limit.
Analogously, for conditional means, the ergodic property enables us to
write
where
Here
(7.4.16)
Then the mean of !(XI) under the condition that Xl> X2, ... , Xn ED, becomes
(7.4.18)
The mean of the field f(x) under the condition that xED, Xl> xz, ... , xn E
Q,,; x n+l> ... , xm ED" is represented in the form
n m m
= <D(x) IT Q(x; x,) IT D(x; Xk)-l S!(xI)D(x) IT D(x; xk)dp(x).
1=1 k=n+1 k=n+1
(7.4.20)
Note that it will always be clear from the context, which field (of cracks
or of inclusions) is considered; therefore the left-hand sides of (7.4.4 and 18)
have the same notations.
For ergodic fields of cracks, calculations of their means with respect to
the ensemble of realization can be substituted by averaging throughout the
space for a fixed typical realization D(x). In this case, the mean <f(XI)lxI;
Xl + xz) for example, can be found from the relation, see (7.4.15),
<f(XI)lxI; Xl + xz)
= <D(XI)D(XI; Xl + xz)-llim.l
V G
v-->oo
S!(XI)Q(XI)D(XI; Xl + Xz) dXI, (7.4.21)
where
(7.4.22)
In conclusion, let us note the formula for the mean <D(x)lxI), which is
analogous to (7.4.10),
(7.4.23)
For homogeneous fields of cracks, both terms on the right-hand side are
functions of the difference X - Xl, the first term being equal to zero when
X = Xl> and the support of the second one is concentrated in a finite neigh-
borhood of the origin, if the dimensions of the cracks are bound.
The construction of means of the type (7.4.22, 23) for specific stochastic
fields of cracks is considered in Appendix A4.
Let us begin with constructing the main statistical characteristics of the effective
field O'(x). We denote by O'n(Xl> ... , xn) an n-tuple moment of the effective field,
i.e. the mean of the tensor product O'(XI) ..... O'(xn) under the condition that the
178 7. Elastic Medium with Random Fields of Inhomogeneities
points XI. ... , Xn belong to V(or ()). In particular, the expectation and the two-
point moment of the effective field are conditional means of the form
(7.5.1)
In order to obtain an expression for ii', let us average both sides of(7.3.17)
with respect to the ensemble of realizations of the random field of inhomogen-
eities under the condition x E V
Since x is a fixed point here, So(x - x') is a determined kernel. Using the
hypothesis H2 (Sect. 7.3) about the statistical independence of the field ii(X)
at the point x of the elastic properties, the form and the dimension of the in-
clusion located at this point, let us represent the mean in the integrand in(7.5.2)
in the form of a product of means, analogously to (7.4.9). Then the expression
for iiI becomes
where the mean (ii(x')lx'; x) is calculated under the condition that x' E V and
XE Vx "
In order to construct the mean (ii(x')lx'; x) in the integrand of (7.5.3),
let us again consider (7.3.17). Averaging both sides of the equation under the
condition x E V, Xl E Vx and using the hypothesis H2 (Sect. 7.3) once more,
we have
The equations (7.5.3, 4) are not closed, since their right-hand sides contain
the mean values of the effective field, which are calculated under conditions,
different from those in the left-hand sides. Thus, a chain of equations arises
here, and to close it, one has to introduce additional assumptions about the
means in (7.5.3,4) or in analogous relations which connect more complicated
conditional means.
As it was pointed out in Sect. 7.3, the classical scheme of the method of self-
consistent field consists of assuming ii(x) to be constant and equal for all par-
ticles. Then ii' and the conditional mean (ii(X') lx' ; x) coincide with each other
and with ii, and the constant ii is given by (7.5.3).
In the scope of the modification (of the method) developed here, the equiv-
alent assumption which enables one to obtain equations for expectation of the
7.5 General Scheme for Constructing First Statistical Moments of the Solution 179
Substituting this relation into (7.5.3) we arrive at the equation with respect
to 0'1
The function €p(x' - Xl) is the mean value of the effective field at the point x'
under the condition that there is a defect at the point Xl> and this characterizes
the pair intt.raction in the system of interacting defects. It is obvious that,
as lx' - xli --+ 00, this function tends to a mean value of the effective field 0'1.
The equation for €p(x) follows from (7.5.4, 7) and has the form
€p(x) = 110 + f So(x - x') (P(x')V(x; x')lx; x', 0) qJ(x') dx' . (7.5.9)
To obtain the second moment 0'2(X1 - X2) of the effective field, let us take
a product of values of the field O'(x) (the relation (7.3.17) at different points
Xl and xz) and average the result under the condition Xl> X2 E V
The equation for the function 0'2(X1 - X2) will be obtained after splitting
the mean in the integrand with the help of the hypothesis H2 (Sect. 7.3)
SHere, the influence of a defect at the point x on the mean value of the effective field at
the point x' is not taken into account.
180 7. Elastic Medium with Random Fields of Inhomogeneities
(7.5.12)
(7.5.16)
B* = Bo - (P(x)V(x) A, (7.5.17)
7.5 General Scheme for Constructing First Statistical Moments of the Solution 181
where the fourth rank tensor A connects the external field with the first moment
jj1 of the effective field
(0'(XI)0'(X2» = 0'5
+ JSO(XI J
- x') dx' SO(X2 - x") (P(x')P(x") Vex') V(x") )jj2(X' - x") dx" ,
(E(XI)e(X2» = E5 + EO(P(X) V(x» jjl + (P(x) V(X»jjIEo
+ JKO(XI - x')Codx' KO(X2J
-x")Co(P(x')P(x")V(x')V(X"»jj2(X' - x") dx" . (7.5.20)
Analogously, one can write down the expression for mean values of the
energy density (rJJ) = (O'E)/2, of a composite medium, of the density (rJJint )
of defects interaction, etc.
When considering a random field of cracks, one has to substitute the func-
tions P(x') Vex') and P(x') Vex; x') into the preceding relations for P(X')Q(x')
and P(x')Q(x; x'), respectively, where now an arbitrary continuous function
P(x) coincides with the functions Pix), (7.3.9), on the surfaces Qk' Thus, to
construct the first two moments of the random fields O'(x) and E(X), one has
to find the statistical moments jj1 and jj2 of the effective field and then to con-
struct the means and to compute the integrals appearing in (7.5.17, 20).
In conclusion, let us compare the method of the self-consistent field with
another self-consistent method of constructing effective elastic constants,
namely, with the method of the effective medium. It was pointed out in Sect.
7.1 that the main assumption, on which the second is based, is that each inclu-
sion in a composite behaves as an isolated one in a homogeneous medium,
properties of the latter coinciding with the effective properties of the whole com-
posite. The external field for each inclusion is assumed to be 0'0'
Using these assumptions, we can transform (7.3.4) for stress and strain fields,
in the case of ellipsoidal inclusions, to the form
182 7. Elastic Medium with Random Fields of Inhomogeneities
where the tensors P k * are expressed through yet unknown effective elastic pa-
rameters B* of the composite
(7.5.22)
Here Clk * is a perturbation of the elastic modulus for the k-th inclusion with
respect to the effective tensor C* = B;l
(7.5.23)
(7.5.16)
Having averaged(7.5.1 l)for the stress and strain tensors and acting then on
the constant, which is the mean of the sum in the integrand, by the operators
Ko and So, we obtain from (7.5.16)
(7.5.24)
where
(7.5.25)
The relation (7.5.24) which connects the constants of the tensor C* with
each other, is to be considered an equation with respect to effective elastic
constants of the composite. The latter are contained on the right-hand side
through the tensor Clk *, (7.5.23), and A k *.
In the case of elliptic cracks, assumptions of the method of the effective
medium lead to expressions for stress and strain tensors, which coincide with
(7.3.13), if one sets if k = 0'0 and substitutes Pk(x) by Pk*(x), while Pk*(x) (an
analog to Pix) in (7.3.8)) is to be found from the solution of the problem
about an elliptic crack in a homogeneous medium with the elastic moduli
tensor C*. Using (7.5.16), one can obtain equations for the effective tensor
B* for a medium with cracks in the form
7.6 Random Field of Ellipsoidal Inhomogeneities 183
B* = Bo - (P *h(x)Q(x» , (7.5.26)
where
(7.5.27)
Let us now consider in more detail the method of the effective field for a medium,
which contains inclusions of the ellipsoidal form. In this section, the ex-
pectation 0'1 of the effective field and the tensor B* of the effective elastic
compliance in the case of composite and polycrystal media will be derived.
From (7.5.17, 18), it follows that the problem of determining the tensors
0'1 and B* is reduced to constructing explicit expressions for the means
(7.6.1)
and calculating the integral which represents the action of the generalized
functions So(x) on the function ?J!(x) defined by the last equality. The tensor
P" in (7.6.1) has the form of (7.3.6).
Let us start with constructing the means (7.6.1, 2). Due to the ergodic
property assumed, the mean (7.6.1) for homogeneous fields of nonintersecting
inclusions, has the form
(7.6.3)
where v" is the volume of the k-th inclusion, N is the number of inclusions
184 7. Elastic Medium with Random Fields of Inhomogeneities
in the region G. It is taken into account here that P k is a constant tensor inside
each of the inclusions.
If one considers the ensemble of realizations of the homogeneous random
field P(x) Vex), then the tensors PkVk will be stochastic variables with one and
the same distribution function for any k. Therefore, if both sides of (7.6.3)
are averaged with respect to the ensemble of realizations once more, then
where Va is the mean volume per an inclusion, v(a) is the volume of an ellipsoid
with semi axes aI, a2, a3, pea) is a tensor determined by the form of the ellipsoid
and by the elastic constants of the latter, see (7.3.6). On the right-hand side
of (7.6.4), we have already a mean value of the random variable v(a)P(a) ,
the distribution function of which is determir.ed by the statistical properties
of the ensemble of inclusions and is assumed to be known.
The problem of constructing the mean (7.6.2) for a homogeneous random
field of elliptic regions is considered in Appendix A4. Let us point out the
main properties of the function 1f!(x). This is a continuous, bound function,
its values being tensors of rank four. If the dimensions of all inclusions are
bound, then 1f!(x) can be represented as a sum of a constant component <1f!),
an oscilating function with the zero mean value and a function with bounded
<
support. The constant 1f!) is easily calculated, if one takes into account that
the constant components of the functions 1f!(x) and <P(x') V(x')lx) coincide,
since, for inclusions of bound dimensions, these functions differ from each
other in a finite neighborhood of the origin only. Using the ergodic property
and averaging later with respect to the ensemble, we obtain
where the definition of the conditional mean is taken into account (Sect. 7.4):
of the matrix being different from the symmetry of the effective constants
tensor of the inhomogeneous medium. For many stochastic composites, the
symmetry of the textUle connected with geometrical properties of the field
of inhomogeneities is not too low and this can be described with the help of
a second order tensor a. In such a simple case, a determines a linear transform,
which transforms 1f!(x) into a spherically symmetric function [1f!(ax) = 1f!'(lx'I)].
In the general case, such a linear transformation can not be found.
The function 1f!(x) is contained in (7.5.17, 18) for 0'1 and B* due to the
integral, which represents the action of the generalized function So(x) on 1f!(x).
Using properties of this generalized function (Appendix A3 and Sect. 4.6) and
the definition of the operator So acting on constants (7.2.6), we have
where
Substituting this relations into (7.5.17, 18) and taking into account (7.6.4),
we obtain expressions for the tensors 0'1 and B* in the form
(7.6.10)
where
(7.6.12)
where A is defined by (7.3.7) with dif3 = o«f3. It is easy to verify that the formal,
limiting transition, as p -+ 1, leads to a physically reasonable result: B* -+
Bo + Bl . Thus, the results obtained are noncontradictory in the entire region
of the concentration of inclusions, though, while p is close to 1, the hypotheses
of the method of an effective field lose sense, since one can not speak of isolated
inhomogeneities in this case.
Note one more important case. In the scope of the method presented, one
can consider not only effective elastic constants, but also more delicate char-
acteristics of a stress state in a composite. Indeed, by virtue of the hypothesis
HI (Sect. 7.3), each inclusion in the composite behaves as an isolated one
in a homogeneous matrix under a constant field iik, k being the number of
the inclusion. From the solution for an isolated ellipsoidal inhomogeneity
(Sect. 4.8), an expression follows for the stress tensor oJ inside each inclusion
(7.6.15)
In the first approximation, one can replace the tensor by its mean value iil,
namely (7.6.10). From here, it is easy to evaluate the stress concentration in
the matrix on the boundaries of inclusions, making use of (4.7.16) [7.6].
Let us apply the method of an effective field to consider themost important
class of inhomogeneous media, namely to one-phase polycrystals. In compar-
7.6 Random Field of Ellipsoidal Inhomogeneities 187
The tensor BI(x) has a constant value inside each crystallite, but changes
from one crystallite to another in accordance with a change of crystallographic
axes.
Following the scope of the method of an effective field, let us assume that
each crystallite k behaves as an isolated ellipsoidal inhomogeneity with the
elastic properties Bo + Blk in a homogeneous medium Bo and an external
field if k' the latter being the sum of an external field if k and a perturbation
caused by all other cristallites. In the first approximation, expressions for
ifl and the tensor of effective compliance B* still have the form of (7.5.18, 17).
Since now the inclusions (crystallites) fill the whole space,
(7.6.17)
where Vk(x) is the characteristic function of the region occupied by the k-th
crystallite. From here, there follow the equalities
where, on the right-hand side we have the mean of the stochastic variable
Pea), the latter still being of the form (7.6.12). It is determined by the distribu-
tion function of the orientation of the crystallites anisotropy axes and by the
form of the crystallites.
The properties of the function 7Jf(x - x') = <P(x') Vex, x')lx) coincide,
in the case under consideration, with all the properties of the same function
for composites, and its symmetry depends on the geometry of the random
field of crystallites, which form the polycrystal.
If there exists a linear transformation a, which converts the function 7Jf(x)
into a symmetric one 7Jf(ax) = 7Jf'(lxl), then, as before (7.5.17, 18), emply
the equalities
posing the tensor B(x) into the constant Bo and the perturbation B I (x), see
(7.6.16); however, the value ifi depends on this decomposition.
In conclusion, let us dwell on (7.5.24) for the effective elastic compliance,
obtained with the help of the method of an effective medium. We consider,
as above, the case of composite and polycrystalline media. Carrying out the
averaging in (7.5.24) according to the scheme presented in this section, we
have
(7.6.21)
As it was already pointed out, the last relation can be considered an equa-
tion with respect to effective elastic parameters of the medium. The right-
hand side contains components of the tensor B* through the tensors C I *
and A *(a) (Sect. 7.5). In the case of polycrystals, one has to assume additionally,
that v(a)Jvo = 1.
The equation (7.6.21) was investigated in [7.7-10]. Values of the effective
elastic compliance, obtained from its solution, satisfie rather delicate variation-
al extimates of Hashin-Shtrikman [7.11], and, in some cases, coincide well
with experimental data. However, for media, which contain cavities (pores),
the method of an effective medium leads to negative values of the effective
elastic moduli, when the concentration of pores is larger than in the case of
close packing. The method of an effective field has an explicit advantage in
this case: as it was pointed out, (7.6.11) for B* is physically non-contradictorary
up to a concentration of inclusions, which is equal to unity.
Let us apply the method of an effective field for calculating the elastic macro-
constants of composites. This problem is of particular interest, since in the
two-dimensional case, for a number of regular structures there are explicit
solutions, which enable us to estimate the error of the method.
Let identical ellipsoidal inclusions form a regular lattice in a homogeneous
matrix. The tensor of effective compliance of the composite is defined as before,
by
(7.7.1)
while averaging of the fields e(x) and o'(x) is performed here with respect to
an elementary cell of the lattice, that is obviously equivalent to averaging
with respect to the whole space. The means of e(x) and o{x) in turn coincide
with the means at a fixed point x under all possible translations of the regular
7.7 Regular Structures 189
where m = maea is the vector of the lattice formed by the centers of inclusions
(Sect. 2.1), the prime attached to the summation sign means excluding the
item with m = 0, and the function lex) is defined by the relation
From (7.7.4, 5), it follows that the function 1lf(x) is the sum of a periodic
function with the mean value pP(a) and a function -l(x)P(a) with bounded
support. The action of the generalized function on such functions is defined
by (7.6.7), from where it follows that
(7.7.8)
(7.7.11)
(7.7.12)
(7.7.13)
Here Po and ))0 are the shear modulus and the Poisson ratio of the matrix,
r
the components of the tensor are written in the basis ea of the lattice vectors,
and the tensors Ei are defined in Appendix AI. The scalar function a(p) in
(7.7.11) is illustrated by its graph if Fig. 7.1, which was obtained by numerical
summation of a converging series of integrals.
0.08
0.06
0.04
0.02
(7.6.11) in the case of stochastic fields of inclusions and of (7.7.10) for regular
structures. All tensors in these relations are to be replaced by their two-dimen-
sional analogs. The two-dimensional problem was considered in more details
in [7.12].
Let isotropic circular inhomogeneities form a regular lattice in an isotropic
(two-dimensational) matrix. The following function is the two-dimensional
analog of the function J(x), (7.7.6),
(7.7.15)
which coincides with the expression for B* in the case of an isotropic field of
inclusions, see (7.6.13). Then the tensor A for a two-dimensional stress state
if defined by
where E1, E2 are the basis tensors, defined in Appendix A4, but for the two-
dimensional case.
If inclusions form a square lattice, then the tensor B* is defined by a rela-
tion of the form of (7.7.11), where A is the tensor (7.7.16), and
(7.7.17)
+t
6 6
~ ~+E1,'
I
5 !..*
Eo ~ -.
"-0 5 l!...*
fto ~'.
~+E1=O
td
"
kt
4 4
b 10.1 b
3 3 I
2 2
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0.2 0.4 0.6 0.8 b/a 0.2 0.4 0.6 0.8 b/a
Fig. 7.2. Comparison of exact and effective Fig. 7.3. Comparison of exact and effective
Young's moduli shear moduli
tities for inclusions, E*, and f1.* are the values of the effective elastic modulus.
The solid curves represent the exact solutions, the dotted line with small circles
corresponds to (7.7.15)
For this case, let us consider the error !1 which arises when evaluating ef-
fective constants with the help of the method of the effective field
(7.7.18)
Here EI and E* are the exact values of the effective Young modulus and the one
calculated from (7.7.15), respectively. Calculations show that !1 is a function
of the dimensionless parameter
(7.7.19)
The dotted line is an extrapolated value of the error. The fact that, as p -+ I,
(7.7.15) gives a correct result for all Ch is taken into account.
The relative error for the effective shear modulus practically coincides with
A(X) for all values of the parameter X.
A comparison of the effective field approximation with exact solutions for a
square lattice of circular inclusions was carried out in [7.12]. The form of the
error of calculating the effective elastic constants is here the same as in the case
of a triangular lattice.
(7.8.1)
B* = Bo - (P(x)Q(x)A, (7.8.2)
where
Here the functions Q(x) and Q(x; x') are defined by (7.4.17, 19), the continuous
tensor field P(x) coincides with P,,(x), (7.3.9), on the surfaces of the cracks.
To construct the means in (7.8.2,8) let us use the same procedure as in the
determination of the means (7.6.1, 2) in the case of ellipsoidal inclusions. Using
194 7. Elastic Medium with Random Fields ofInhomogeneities
the ergodic property with the subsequent averaging with respect to the ensemble,
we obtain
The constant components (1ff) of the functions 1ff(x) and 1ff(x - x') =
(P(x')Q(x')lx) coincide, since their difference is a function, with bounded sup-
port, if the dimensions of the cracks are bound (Appendix A4). By analogy
with (7.6.5), we have
are taken into account, the latter being easily obtained, as (7.8.4) was. For cracks
of identical dimensions
21t a 3 (7.8.9)
(1ff) = -3 ~ (P(a, b),
Vo
i.e. ?P'(x) = const. By virtue of the definition (7.2.6) of the action of the operator
So on constants
(7.8.13)
(7.8.14)
where the tensors Ei are defined by (AI) and the averaging of these tensors
is carried out with respect to all possible orientations of cracks. If the distribu-
tion with respect to orientations is homogeneous, then
<a3 )
B* = Bo - - 8 (1 - ].10) --[].loB 2
- 2(5 - ].I0)E1], (7.8.15)
45 llo(2 - ].10) Vo
and, in particular, for the effective shear modulus Il* and the Poisson coeffi-
cient ].1*, we have
196 7. Elastic Medium with Random Fields oflnhomogeneities
_
fl* - flo
(1 + 32 (a 3 ) (1 - ).10) (5 - ).10)
45 -v;;- (2 - ).10)
)-1 , (7.8.16)
(7.8.19)
Note that the expression for B* will have the same form, if there exists a
linear transformation a, which makes lff(x) a spherically symmetric function
lff'(lxD = lff(ax). Then the tensor D depends on a and is defined by the relation
(7.6.8).
In the case of a spherically symmetric lff and of a homogeneous distribution
with respect to orientations, B* is an isotropic tensor
(7.8.20)
where
21t a3 (7.8.22)
(P(x)Q(x» = (lff) = -3 - P(a, b) ,
Vo
7.8 Fields of Elliptic Cracks 197
where P(a, b) is a constant tensor of the form of (7.3.11) and the function h(x)
has the form of (7.3.10) at each crack. In these relations, averaging is carried
out with respect to all translations, i.e. analogously to the case of regular
lattices of inclusions (Section 7.7).
As it is shown in AppendixA4, the mean (7.8.23) for regular lattices of elliptic
cracks is represented by
where m is a vector of the lattice formed by the centers of the cracks, the primed
summation sign means excluding the term with m = 0, Q;"(x) is the a-function
concentrated on a flat surface with its center at the point m and parallel to
all cracks. The form of the surface is an ellips with semiaxes, which are equal to
doubled semiaxes of the crack under consideration. The function J(x) is defined
by the relation (A 4.47)
where
(7.8.26)
and the xL and x 2-axes are directed along the principal axes of the elliptic cracks.
From (7.8.10) we can now find the result of the action of the generalized
function So(x) on 1/f(x) in the form of (7.8.24)
(7.8.28)
where D is defined by (7.6.8) with a«~ = aa~. Substituting (7.8.27) into (7.8.1,
2) for the tensors 0"1 and B*, we obtain
(7.8.29)
198 7. Elastic Medium with Random Fields of Inhomogeneities
where (W) has the form of (7.8.22). Note that, for simple lattices, iiI is the
real external field o{x) averaged over the surface of any crack.
If cracks form a complex lattice, then the expression for W(x) consists of
sums of the type of (7.8.24) over each sublattices. In this case, one can make
the method of solution more precise by assuming that the value iiI of the
effective field is different for different lattices. Equations for the set of quan-
tities iii (the values of iiI on the k-th sublattice) follow from an analog of
(7.5.6) for the case of cracks and become
n
iii + 2: Q}kl (Wi) iii = 0"0, k = 1, 2, ... , n , (7.8.31)
i=1
where
and, with i = k, the expression for Q)!'l has the form of (7.8.28), where the
quantities a, band m concern the k-th sublattice only. In these relations, m..
is a vector of the i-th sublattice, ai, bi are semi axes of elliptic cracks in the
i-th sublattice, ~ik is a vector, which connects knots of i-th and k-th sublattices
in an elementary cell, the tensor (Wi) is determined by the orientation and
the dimensions of cracks in the i-th sublattice and has the form, analogous
to (7.8.22).
Solving the system (7.8.31) with respect to the tensors ii} and substituting
the result into the expressions for the mean values of stress and strain, analogous
to (7.5.14), using then the definitions (7.5.16) for the tensor B* we obtain
n
B* = Bo - 2:
i=1
(Wi) Ai , (7.8.33)
where the tensors Ai connect ii} with the external field 0"0
(7.8.34)
where
and the radius of all cracks is a. The dimensionless scalar coefficients ao and
(30 depend only on the geometry of the lattice and are presented in the form
of converging series of the type of (7.8.28). Fig. 7.5, displays the graphs of
0.20
0.15
0.10
0.05
ao and {30 as functions of the parameter q = (a3)lvo for a regular cubic lattice
of cracks, all planes of the cracks being parallel to a fixed face of the cube.
The stress intensity factors on the contour of a crack are found from the
solution of the problem for a circular isolated crack in the constant field jji
of the form of (7.8.35). In the case of pure tension in the direction of the normal
200 7. Elastic Medium with Random Fields of Inhomogeneities
to the plane of cracks, lTff,B = lTonanf3, and the stress intensity factor K J has
the form [7.8]
(7.8.38)
1.2
(7.8.39)
calculated with the help of the method of effective field (dotted curve) is com-
pared with the maximun value of kJ on the edge of the cracks from [7.14]
(continuous curve).
In conclusion, let us consider the expression for the tensor of the effective
elastic compliance of a medium with cracks, to which the method of effective
medium leads. Eq. (7.5.26) for the effective compliance is represented by
(7.8.40)
where P * is determined from the solution for an elliptical crack in the medium
with the tensor B* of elastic compliance and has the form, which is analogous
to that of Ph in (7.3.8) (for more details, see the end of Sect. 7.5).
7.9 Two-Dimensional Systems of Rectilinear Cuts 201
B* = Bo - 2'Jt<
3 Vo a3 ) p * , where (7.8.41)
(7.8.43)
(7.8.44)
where the xLaxis of the cartesian coordinate system with the origin at the
center of the cut is directed along L; I is the half length of the cut. The constant
fourth order tensor P has the form
(7.9.4)
Here SoCk) is the k-representation of the Green's function for stress in the
two-dimensional problem, n is the normal to L, JI(k l ) is the Bessel function
of the first order, the Greek subscripts in this section take values 1, 2. For
an isotropic medium and for a plane state of stress
(7.9.5)
with x E L
(7.9.6)
with x E Li ,
where the continuous function P(X)if(X) = M(x) has the form of (7.9.2) on
each of the cuts L k. The expressions for the tensors l1(x) and c(x) are represented
by
7.9 Two-Dimensional Systems of Rectilinear Cuts 203
(7.9.10)
B* = Bo - (P(x)L(x»A, (7.9.11)
where
where c:vo is an average area per crack, the tensor P has the form of (7.9.3) and
depends on a stochastic orientation n of the crack. The mean in the integrand
in (7.9.12)
possesses all the properties of the mean (7.8.5) of the three-dimensional case.
The constant component of the function UJ(x) is equal to
(UJ) = ~ (/3P) ~
2c:vo ( I ) 2
1..(1tc:vo12) (P), (7.9.15)
while the last equality occurs if the dimensions of all the cracks are approxi-
mately identical.
A regularization of the integral in(7.9.12)is defined by a relation analogous
to (7.8.10). In the two-dimensional case, the tensor D becomes (Appendix A3)
where the integration is carried out over the unit circle in the two-dimensional
k-space. For an isotropic medium, the tensor D takes the form
204 7. Elastic Medium with Random Fields of Inhomogeneities
Here El and E2 are the two-dimensional analogs of the basis tensors, defined
by (AI).
Let us proceed to consider specific random fields of cracks on the plane.
1) A Poisson field of rectilinear cuts is defined analogously to a Poisson field
of elliptic cracks in the three-dimensional case (Sect. 7.8), and here
B* = Bo - -1
2
< 12
'It - P ) .
-
(vo
(7.9.20)
(7.9.21)
(7.9.22)
_ <'It /2)
q---. (7.9.23)
(VO
This implies
(7.9.24)
where Eo and E* are the Young moduli of the real and effective media, respec-
tively.
2) A model with a restriction on the intersections of cracks (Sect. 7.8). In
this case, the mean of W(x) defined by (7.9.14) is a continuous function, which
is equal to zero at x = O. If W(x) is also spherically symmetric, thus
<
where D has the form of (7.9.16), and 1Jf) has the form of (7.9.15). In this case,
the tensors 0'1 and B* are represented by
(7.9.28)
B* = Bo + 2 (1
Po
1
+ Vo) (1 - Tq
1 )q( 1 - 41)q [E1 + ~ q(E2 - E1) l
(7.9.29)
Hence, for the relative elastic modulus E*/ Eo and the Poisson coefficient v*/vo,
we have
(7.9.30)
(7.9.31)
(curve 3) are compared with experimental data [7.16]. The experiments were
carried out using thin rubber sheets, containing random fields of rectilinear
cuts (vo = 0.5). Results of the experiments are exhibited in Fig. 7.7, by the dot
line with small circles (E*/ Eo) and by the dot-and-dash line (v*/vo).
The statistical analysis shows that, to the field of cracks investigated in [7.16],
there corresponds the model with the restriction on cracks intersection. The
experimental data are described by (7.9.30, 31), obtained for this model, in the
best way. Therefore one can state that the approximation of an effective field
enables one to take into account correctly interaction between cracks in the
situation under consideration.
Let us present here expressions for E* and v*, which are given by the meth-
od of an effective medium in the case of an isotropic tensor of effective elastic
compliance. The analog of (7.8.41) in the two-dimensional case has the form
(7.9.32)
Thus, we have
E* = ~ = 1 _ q. (7.9.33)
Eo Vo
where m is a vector of the lattice formed by the centers of cracks, the primed
summation sign means the exclusion of the term with m = 0, I is half the length
of each crack, L;"(x) is the a-function concentrated on the segment of the length
41 with its center in the m-th knot and parallel to all the cracks. The function
lex) has the form
lex) =
2~ {(I - M) /1&(2 - 1&) + sin- (1-1&)
I Y ' I I
+~
2
1 if Ixl :::;; 21
11: 0 if Ixl > 21
(7.9.35)
The action of the generalized function So(x) on l/f(x) in the fOlm of (7.9.34)
is represented as in (7.8.27),
Here D has the form of(7.9.16) and the mean value of1ff(x) is
7C1 2
(1ff) = -2-P(n). (7.9.38)
000
Expressions for the tensors q1 and B* have in this case the form which is an-
alogous to (7.8.29, 30).
Let us consider a lattice of cracks of one orientation, which has two or-
thogonal symmetry axes. The expression for the vector q 1n (which, here, deter-
mines the state of all cracks) has the form
-Up
0' np = (1 _1 (3) Up (l:.t + 1_ a .l)
a _ {3 n np e..tnv •
0'0 (7.9.39)
The expression for the tensor of effective elastic compliance B* in terms of the
basis tensors Ei defined in Appendix AI, is
(7.9.40)
(7.9.41)
where 0'0 and 't"o are the values of external tensile stress and shear stress, re-
spectively,
1
kJ = 1 - a - (3 , k n = 1 - a (7.9.42)
1.0 2.0,--------------,0--,
1.5
OL-----~~------~
Fig. 7.8. Elastic moduli for triangle Fig. 7.9. Stress intensity factors for
lattice of cracks triangle lattice of cracks
b) A system of collinear cracks on one straight line. The values of the co-
efficients a and (3 in (7.9.42), which correspond to this case, can be obtained from
the solution for a rectangular lattice when one of the characteristic dimensions
of the latter tends to infinity. Then a and (3 are expressed as simple series [7.17].
The value of the coefficient kJ is in agreement with its exact value, which has
the form [7.19]
2 sin ~ A (7.9.43)
k9=-~~=~
.y' ATe sin TeA
(c) A row of parallel cracks. In Fig. 7.10, the values of the intensity factors
kI and kIl' obtained by the method of an effective field (dotted line with
small circles), are compared with results of numerical solution of the corres-
ponding integral equations [7.20] (solid curves).
kll k J
2.0 1.0
1.8
1.6
1.4
1.2 Fig. 7.10. Stress intensity factors for a row of
parallel cracks
0
As it can be seen from the above examples, in the two-dimensional case, the
first approximation of the method of an effective field gives good agreement
7.10 Random Field of Point Defects 209
with the results of more precise calculations of elastic constants and of in-
tensity factors, when the lengths of the cracks are smaller than the distances
between their centers. This statement appears to remain valid for other types of
regular systems of cracks as well.
When the centers of cracks in a lattice become closer (see case c), the ex-
ternal field q(x), in which each crack is contained, begins to differ more and
more from a constant one. Therefore, the error which is given here by the
method of the effective field, is connected with violating the hypothesis HI
(Sect. 7.3). However, a way of making the computational part of the method
more precise can be suggested.
Let the effective field q(x) be not constant and have the form of a polynomial.
The solution for an isolated inhomogeneity can be found, using results of
Sect. 4.9, 12. This will enable one to express the corresponding induced density
m(x) of dislocation moments on defects in terms of the coefficients of the
polynomial q(x).
Equations for these coefficients can be obtained, substituting q(x) and ex-
plicit expressions for m(x) into (7.3.3) for the effective field and averaging a
necessary number of times the expression obtained, over the region occupied
by defects. As weights for the averaging, a system of linearly independent func-
tions (polynomials, for example) is to be chosen. An analogous modification of
the method can be used also in the case of stochastic fields of ellipsoidal in-
clusions and cracks.
In the previous sections the effective field method was applied only to the
determination of the first moments of stress-strain fields in composite materials.
The calculation of the second moments is technically a more complicated
problem. Finding the averages in (7.5.9, 13), determined by the second mo-
ments of q(x), and the solution of the system of these equations is only possible
by using numerical methods, even for simple cases. However, many technical
difficulties disappear if we pass to the model of point inhomogeneities, where
inclusions of finite sizes are replaced by isolated point defects.
In this and the following section the general scheme of the effective field
method is applied to the description of the elastic medium containing a homo-
geneous random field of point defects. Note that in this case the application
of the method has a better foundation than for inhomogeneities of finite
sizes. As a matter of fact the main hypothesis HI (Sect. 7.3) about the constancy
of the effective field q(x) in the region occupied by an arbitrary inclusion,
turns into a trivial fact in case of point defects.
Let us begin with the problem of replacing an isolated finite inclusion by
a point defect. The basis for such a substitution is the following convention.
210 7. Elastic Medium with Random Fields of Inhomogeneities
The perturbed fields of the isolated inclusion and of the corresponding point
defect should have the same asymptotics at infinity. Thus in order to pass to
the point-defects model from the model of the composite with inclusions of
finite sizes, we should substitute the functions PkVk(X) and Pk(X)Ok(X) in the
previous relations by their first terms in the multipole series expansions
(7.10.1)
where ~k is the center of the region Vk • Note that PkVlx) characterizes the
density of dislocation moments induced in the region V k by applying an ex-
ternal field. The tensor P2 in (7.10.1) has the form
(7.10.2)
In case of cracks we deal with the function PiX)Ok(X) which is also re-
presented in the form of (7.10.1), where
(7.10.3)
Here Pk(x) and P k are defined by (7.3.9, 11), ak is the larger half-axis of the
ellips Ok'
Let X be a set of points ~k with the point defects specified by the set Xxo
and defined by
X
xo
= l x=
U ~i
i*k
U~i
'
by Xo E X
(7.10.4)
by the analogy with V Xo ' see (7.3.14). We consider the generalized functions
O'(x) = 0"0 + SSo(x - x') PO(x')O'(x') X(x; x') dx', x EX, (7.10.7)
7.10 Random Field of Point Defects 211
where Po(x) is a continuous function coinciding with P2, see (7.10.2, 3), at the
point x = ~k'
For the determination of the effective field's first moment 6'1 both parts
of (7.10.7) should be averaged with the convention x E X. Then the hypothesis
H2 (Sect. 7.3) about the statistical independence of O'(x) from PO(x')X(x; x')
at the point x together with an assumption similar to (7.5.5), should be used.
As a result the equation for 0'1 is presented by
where the operation of the average with the conventions x E X and Xl> •.. ,
Xn E Xx, is denoted by (Ix; Xl, ... , x n ), see also (Sect. 7.4).
In case of the point-defects model, (7.5.17) for the effective compliance
tensor B* is reduced to
B* = Bo - (PO(x)X(x» A , (7.10.9)
where P = (PO) is the average over the ensemble of random variables (7.10.2
or 3).
Let us consider two stochastic point-detects models in space (a Poisson
field and a random point lattice) and determine the tensors 6'1 and B* in these
cases.
1) The Poisson point field. An example is the Poisson field of the set of
the centers of the cracks described in Sect. 7.8.
In case of the Poisson field the averages on the right-hand sides of (7.10.11,
12) have the following forms, see (A4.6, 17),
1 (7.10.14)
B* = Bo - -P.
Vo
where the tensor D is defined by (7.6.5) and aa(3 = oa(3. This relation together
with (7.10.9) yields
B* = Bo - -1 P ( I
Vo
-1 DP
+ Vo )-1 . (7.10.16)
It coincides with (7.6.11) in the case of inclusions and with (7.8.19) in the
case of cracks.
Thus, the point-defects model with the necessary restrictions yields the
same expressions for the effective constants tensors as the model for finite
sizes of inhomogeneities.
7.10 Random Field of Point Defects 213
~m = m + Pm + r, (7.10.17)
where the primed summation sign means the omission of the term m = O.
The functions g(x) and f(x) are connected with each other by
Note that (lff) = l/v and this equality may be proved as (7.6.5). Thus
(7.10.21) is a consequence of (7.10.19) and (7.2.6).
In order to calculate the integral (7.10.21) let us pass to the k-representation
of the functions in the integrand. Using Parseval's formula we obtain
1
Q = S So(x)[volff(x) - 1] dx = (21t)3 S So(k)[volff(k) - (21tpo(k)] dk.
(7.10.22)
vo1/f(k) - (2rc)3 0(k) = exp( - J a2k 2)[(2rc)3 ~ o(k - 2rcm') -1], (7.10.24)
If the expectation of the random point field is the cubical lattice these
integrals become
(7.10.25)
where the components of the tensor r are defined in the system of the basic
vectors of the cubical lattice by
3
r(X~).1' = .t...J
'" u,uJ:qJ:l$o~ol!
t t t· (7.10.26)
;=1
The coefficients ao, al> a2 have the form of the convergent series (1 is the
lattice parameter)
a1 = t
;,j,k=-oo
(i2 + J('~~ (P
l
.i k 2) })~2+ i2) exp[-2( rcla )(i2 + p + k 2)]
+} + k
(7.10.27)
B* = Bo - -1 P ( I - -QP
Vo
1
Vo
)-1 (7.10.29)
6
5 1-'-*
4 fLo
3
2
1
0.8
0.6
0.4
0.2 Fig. 7.11. Shear modulus for square lattices of
holes (1) and rigid inclusions (2)
0
216 7. Elastic Medium with Random Fields of Inhomogeneities
(7.11.2)
The function
(7.11.4)
is the mean under the condition that Xl> X2 E X, where X is the set of points, at
which isolated defects are located. Obviously, if Xl i= X2, then in Xl and Xz
only two different defects can occur. Therefore, jj2(XI - X2) pertains to the
pair interaction in a random field of point defects.
Along with the two-point moment of the effective field (7.11.4), let us intro-
duce the mean of jj2 at a point Xl under the condition that there is a defect at X2
(7.11.5)
7.11 Correlation Functions in the Approximation by Point Defects 217
The limit of d(x), as Ixl -+ 00, will be denoted by doo- Since at large distances
the dependence on X2 in (7.11.5) disappears, we have
The expression for d(x) can be obtained analogously to (7.11.1) and has the
form
Our next aim is to construct closed equations for the functions 0'2(X), d(x)
and q;(x). Here we shall make use of (7.11.1, 2 and 7), uncoupling the means in
the integrands with the help of statistical hypotheses of the type of H2 (Sect.
7.3) and some additional assumptions concerning the structure of conditional
means.
Let X X1X2 be X with points ~j, which appeared at Xl or X2. If
(7.11.8)
then, by the virtue of the definition (7.10.6) of the function X(xt. x'), we have
the equality
the expression for each term on the right-hand side of (7.11.10) can be repre-
sented in the form
(7.11.13)
218 7. Elastic Medium with Random Fields of Inhomogeneities
where P is the mean with respect to the ensemble of the stochastic variables
(7.10.2,3).
Substituting these expressions into (7.11.10) and then the result into (7.11.1),
we find
where
(7.11.15)
The equation for the function 9'(x) follows from (7.11.3), if one represents
the mean in the integrand in the form
(PO(X') X(XI; x')q(X')IXI; X2) = P(X(Xl> X2; X')IXI; x2)(q(x')lx', Xl> X2)
+ o(x' - X2)P(q(X2)lxI; X2). (7.11.16)
The equations (7.11.14,17 and 18) form a closed system with respect to the
required three functions. A specific structure of the random field of defects is
contained in these equations through the function F(x', Xl> X2), which is de-
fied by (7.11.15). Let us consider this function in more detail.
From (7.11.15) and the definition of the conditional mean (Sect. 7.4), we
obtain
(7.11.19)
This expression for the function F is the point of departure for the construction
7.11 Correlation Functions in the Approximation by Point Defects 219
of the explicit expressions for random point fields of different types. Specific
examples of such constructions are carried out in Appendix A4.
In the case of a Poisson point field, the function F turns out to be equal to
the constant
If the set X is a random spatial lattice (Sect. 7.10), then F is determined from
(A4.30) and has the form
Here m, n are vectors of a fixed regular lattice (the expectation of the random
one under consideration), the primed summation signs mean exclusion of the
terms with m = 0 and n = 0, respectively, the function g(x) is defined by
(7.10.20).
The system of equations (7.11.14, 17 and 18) will be solved for the ex-
ample of the two-dimensional linear chain of defects. Let a system of parallel
rectilinear cuts of a stochastic length 21 located on one straight line be modeled
by point defects. The coordinates of the centers of cuts form a onedimensional
homogeneous random point field. Let us assume that the external stress 0"0 is
a tension along the normal to the line of cuts, i.e.
(7.11.22)
(7.11.24)
na
s,a(J).p( )
0 X np = Po 1... !lpJ.
21t(1 _ vo) x 2 u , (7.11.25)
220 7. Elastic Medium with Random Fields of Inhomogeneities
d(x ) -- 0'0\£>()
X + b2 ----.xz
0'2(x) + b2 SOO D('
(x _.1 X')2 L"' X ,x, 0)0'-2(X') d x,
'
r
-00
If the point-defects density tends to zero, then the integral terms in these
equations vanish, and the equations describe interaction between two isolate
defects on the straight line. At the same time, the solution of the system (7.11.26)
has the form
(7.11.28)
Here the expression for \D(x) is the normal component of the stress field, in which
there are located two point defects, the distance between them being x. It is
obvious, that \D(x) is the asymptotics, as x > b, ofthe solution for the interaction
of two rectilinear cracks, which lie on the same straight line. If x :::;; b, the solu-
tion (7.11.28) has no physical sense.
As pointed out in Sect. 7.10, in the case of finite defects, there is no analog
to the point field, in which defects can occur at arbitrary small distances. For
example, in the case under consideration, the centers of cracks must be not
closer to each other, than the sum of their half lengths. This circumstance can
be partially taken into account, if one makes use of the model of a one-dimen-
sional point field, which is analogous to a random lattice of defects in the space.
Let Xk be a coordinate of the k-th defect and the differences XHI - Xk be
independent stochastic variables (for different k's), which have the same normal
distribution
1 [ ( x - 10)2 ] (7.11.29)
f(x) = .v21t0' exp - 20'2 .
7.11 Correlation Functions in the Approximation by Point Defects 221
where the primed summation sign means exclusion of the term with k = 0,
and the double prime means the same for n = 0 and n = k,
3 [ ( x - xlo)2] (7.11.31)
fix) = .v21tlklo exp - 202 .
Let us substitute (7.11.30) into (7.11.26) and search for the solution of the
latter in the form
X2
9'(x) = (bI(X) X2 _ b2 . (7.11.32)
This choice of the structure of the solution is due to the following reasons.
The functions 9', ff2 and d characterize an interaction between two defects in a
random field of point defects. In the first approximation one can assume that
their form coincides with that of the corresponding functions for two isolated
defects (7.11.28), and the presence of other defects can be taken into account
by correcting the external field. Thus we immediately arrive at (7.11.32).
Numerical computations show that the functions ffi(x), dl(x) and 9'I(X) are
well approximated by constants, the values of the latter depending on the para-
meters p = 2b/lo and x = I/o. Let us denote these constants by ffMp, ,.), doo(p,,.)
and 9'oo(p, ,.), respectively; one can prove that
(7.11.33)
The dependence of 9'oo(p, ,.) and doo(p, ,.) - 9'!(p,,.) on the parameters is shown
in Figs 7.12 and 13.
Since 9'00 = ffl the mean value of the effective field, as it is seen from Fig.
7.12, attains its maximum for large relative dispersion of distance between
222 7. Elastic Medium with Random Fields of Inhomogeneities
cr 1
2.5
d=-\p~
~2-
0.8 0
0.6
0.4
P=l 0.2
0
0.8
-0.2
0.6
0.4 -0.4
10 'X
-0.6
Fig. 7.12. Dependence u,/q. on" andp Fig. 7.13. Dependence (d~ - Y~)M on "
andp
defects (I> = 2.5 ..;- 3.5). With I> = 5 ..;- 6, the quantity 0'1 no more coincides
with the value, which corresponds to a regular chain of defects (I> = (0) .
The difference doo - rp! is equal to the dispersion of the effective field 0'.
From Fig. 7.13, it is seen that the dispersion has its maximum for I> = 4 ..;- 5,
and for I> = 9 ..;- 10, is practically zero, which corresponds to the transition
to a regular structure. For I> < 3, the quantity doo - rp! becomes negative.
This loss of physical meaning is connected with an increase of the probability
of the defects approaching each other at distances less than b for small 1>.
As it is seen from (7.11.32) for 0'2(X), the correlation radius of the effective
field has the order of the average value of the dimension b of a defect and
weakly depends on the relative dispersion of the distance 1>-2 between defects.
Let us now proceed to calculating the second moment of the stress field
o-(x) in a medium with point defects. The relation (7.5.20) becomes
(0"(XVo"(X2» = o-ij
+ JSO(X1 - x')dx' JSO(X2 - x")(PO(x')PO(x")X(x')X(x")O'(x')O'(x"» dx'.
(7.11.34)
and the hypothesis H2 about the statistical independence (Sect. 7.3), we repre-
sent the mean in the integrand in (7.11.34) by
where
(7.11.40)
where
(7.11.41)
d
0'5 + /1t(x) J 1t{x -
00
where doo == d( 00) and a 2(x) have the form of (7.11.32), the function ¢(x) is
represented by
00
The primed summation sign means exclusion of the term with k = 0, and
fk{X) has the form of (7.11.31). The function 1t(x) is an analog of H(x) from
(7.11.37) and is defined by the relation
4It is taken into account here that at one and the same point there cannot be two different
defects; therefore, we have <a(x')a(x")lx, x") = <a(x')a(x')lx') = d« if x' = x".
224 7. Elastic Medium with Random Fields of Inhomogeneities
The graph of the continuous part of the function p(x) - O'~ is presented in
the Fig. 7.14. The presence of a singular part and of a singularity at x = bin
O~~k-~~~~~~~~~~
4 5 x/21
-1
-2
-3
Fig. 7.14. Dependence of correlation radius p on ,,/21
the correlation function O'nn(x) of the random field on the line of defects is con-
nected with replacing real defects by point inhomogeneities. For a random
field of defects with finite dimensions, the correlation function is smooth,
bound and has the minimum correlation radius of the order of the average
dimension of a defect. When the random field of defects tends to a regular
structure, the correlation radius of the stress field increases; this is easily seen
from Fig. 7.14.
7.12 Conclusions
Let us now estimate the accuracy of the method of an effective field and discuss
a possibility of using the first moments of the solution of the stochastic elastic
problem, for strength calculations of structures made of composites.
It was pointed out above that in the general case it is impossible to obtain
an explicit solution of the elastic problem for stochastically inhomogeneous
medium. Each statistical moment of the solution is expressed in terms of a set
of moments of higher orders, and to truncate the chain of equations, one has
to introduce additional assumptions, usually without rigourous proofs of their
validity.
There exists a number of approaches to solving such problems; surveys
can be found, for example, in [7.2, 21-24]. The method of an effective field
considered in this chapter, stands somewhere between self-consistent schemes
and methods of the type of smoothing approximation or stochastic phases
approximation [7.21, 27].
7.12 Conclusions 225
5The latter is, generally speaking, not obligatory. One can assume only an identical be-
havior of the field if in a neighborhood of each particle; see the end of Sect. 7.9.
226 7. Elastic Medium with Random Fidds of Inhomogeneities
the particles. This can occur when particles are situated close enough to each
other, i.e. when the concentration is high. On the other hand, the more the
defects affect each other, the more reasonable the hypothesis H2 is.
However, invalidity of the hypotheses HI and H2 can weakly influence
the first and the second moments of the solution, i.e. the quantities of most
interest for applications. These moments are extremely rough statistical char-
acteristics of a solution, and a part of the information about detailed behavior
of the particles is here lost. Therefore, a detailed consideration of the state
of each particle is not reasonable. Of course, this argument can be confirmed
only by comparing the derived results with exact solutions or experiments.
Since exact solutions of stochastic problems are unknown, comparisons
with experiments play an important role. Fr9m the results of Sect. 7.7 and 9,
it follows that the method of an effective field enables one to describe suffi-
ciently the results of experimental measurements of effective elastic constants
in the two-dimensional case, and, for regular composites, yields a fair agree-
ment with exact solutions.
In Sect. 7.7. a dimensionless scalar parameter X was introduced, with the
help of which the range of applicability of the method can be described. The
size of the region of values X, where the method gives an error more than 10 %,
is relatively small. This region contains, for example, X for the cases of circular
holes or of absolutely rigid inclusions at concentrations of p < 0.4. The esti-
mate seems to be valid not only in the case of a regular triangular lattice of
inclusions, but also for other types of structures, in which distributions of
inclusions differ weakly from the isotropic one. Analogous estimates for
stochastic fields of cracks follow from the results of Sect. 7.9.
In the spatial case, the method of an effective field has qualitatively the
same region of applicability. However, in the three-dimensional medium,
the potential of a particle decays more rapidly than in the planar case ( ,...., r- 1
and ,...., r- 2 , respectively). Another peculiarity of the three-dimensional problem
is, that here more particles can interact directly. All this can exhibit quantita-
tive boundaries for the applicability of the method.
Unfortunately, the absence of sufficient data in the literature, concerning
the determination of correlation functions of elastic fields in an inhomogeneous
medium, does not allow us to estimate the error of the method when calculating
the second moments of the solution.
In conclusion, let us dwell briefly upon the scheme of using the obtained
results for strength computations of structures made of composites. Such a
computation can be divided into two steps: computing the macro-stress and
estimating the level of micro-stress.
When computing the macro-stress, the composite is considered as a homo-
geneous medium, the elastic properties of which are determined by effective
constants. Here one can use the classical elasticity and consider the macro-
stress obtained as an external one when estimating the level of micro-stress
connected with an inhomogeneity of the material.
7.13 Notes 227
7.13 Notes
1
E!{3).!' = T(Oa),0fl!' + oa!'ofl),), E';fl).!' = oafl0),!',
where
(A 1. 3)
It is easy to verify that the tensors Eik can be expressed as linear combinations
of the basis elements Ei.
The linear tensor space, spanned by the basis Ei is closed with respect to the
product introduced, and hence forms an algebra which we denote A(n). We set
Ei.p = a~~jEj,
E4 1 E413E413E61 E41 E6 1 E6
E5 1 E51 E41 E6 1 E41 J(E5 + E6) 1 E6
E6 1 E61 E41 E61 E41 E6 1 E6
Ii = A~Ek, (A 1.5)
0 1 1 1 0 3
2 -2 -2 2
0 1 1 1 0 1
-2 2 2 2
0 0 -1 1 0 3 (A 1.6)
A~ = -2 2
0 0 -1 1 0 1
2 2
1 1 1 1 -2 1
-2 2 2 2
The inverse matrix a~ has the form
1 0 0 0 1 1
3 1 3 1 0 0
2 -2 -2 2
1 1 1 1 0 0
2 2 -2 -2
a~ = (A 1.7)
1 1 - 1 1 0 0
2 2
1 1 0 0 1 0
2 2 2
1 1 0 0 0 0
2 2
Al Fourth-Order Tensors of Special Structure 231
The multiplication table for the basis tensors Ji has a simpler form
I
J1 I J3 I J4 I J5 I J6
I
J2
JI JI J2 I J3 I J4 I 0 0
I I I
J2 J2 JI I-J4 I- J 3 I 0 0
I I I
J3 J3 J4 I Jl I J2 I 0 0 (A1.8)
I I I
J4
I
J4
I
J3 I-J2 1-J1 I 0 I
0
J5 0 0 0 0 J5 0
I I I I I I
J6 0 0 I 0 0 0 J6
I I I I I I
It is seen from the table that the algebra A(n) is split into the direct sum of two
subalgebras with the basis J5, J6 and JI, ... , J4, the second subalgebra being
isomorphic to a real form of the complex quaternion algebra (the usual real
quaternion algebra will be obtained if one substitutes J2 -+ J ~ = iJ2, J3 -+ J~ =
iJ3).
A tensors Qa~).1' is inverse to a tensor P a~).1' if these tensors satisfy the relation
(A 1.9)
Let a tensor P = Pk JR, and a tensor Q = q~k be inverse to P. Then, for the
components of the tensor Q, we have
PI _l!:L
( qb ... , q6) = ( T' 11 '
_l!L _l!.L
11 ' 11' is' ;6), (A 1. 10)
where
(A 1. 11)
It is easily seen that qh ... , q4 are connected with Ph ... , P4, as in the qua-
ternion case, and only difference being that "modulus squared" is not of con-
stant sign. The conditions of existence of the inverse tensor are: 11 =F 0, P5 =F 0,
P6 = o.
A special case, when the tensor P a~).1' is symmetric with respect to permuta-
tion of the pairs a{3 and Af1. is of interest, since the tensor of elastic constants
has this symmetry. One can show that in this case, the inverse tensor Qa~).1'
possesses the same symmetry.
232 Appendices
As usual, one can define a scalar product of two elements of the algebra
A(n) as the complete contraction of the two tensors. Then, the scalar product
of the unit tensor and an aribtrary one coincides with the trace of the latter. In
particular, for the basis elements Jk, we have
J1 I J2 I J3 I J41 J5
(AU2)
Tr Jk I 2 I 0 I 0 I 0 I 2
It is easily seen that Ell(n) belongs to the class of tensors considered in Appendix
AI. In fact, one can directly verify the identity
and taking into account the formulae (AU) for the basis Ei, we find
With the help of the matrix a~ given by the table (A1.7), we can rewrite
(A2.4) in the basis Ji of the algebra A(n)
men) = -.l
2
(Jl - J2) - J6 . (A2.5)
A2 Green's Operators of Elasticity 233
Note the following important identities for &len), being easily verified with
the help of the Table (A 1. 8) :
(A2.6)
Let ,(x) be a second-order tensor field, which vanishes at infinity, and let
,(k) be its Fourier transform. Then, as is easily verified, the following equiv-
alence relations
=0
div , '¢!> no, = 0 , (A2.8)
Rot, = 0 '¢!> eo, = 0
hold.
Let us proceed to the Green's operators of elasticity. The Green's operator
for displacements G(k) is defined by (4.6.4, 5). Its explicit expression in the
isotropic case is given by (2.8.6), where;' and f.i are to be taken as constants.
From the results of Sect. 4.6, it follows that the Green's operator for the
strain K(k) is connected with G(k) by the relation
(A2.9)
KCK= K. (A2.H)
satisfies
SBS = S. (A2.14)
8 0K = 0, DoS = O. (A2.I7)
In order to show the existence and uniqueness of the solution ofthe algebraic
(in the k-representation) system of equations (A2.I6, 17), let us introduce some
useful notations. Let L be an arbitrary linear operator, L1 and L2 its restrictions
to the subspaces of the projections Do and eo, respectively, and
(A2.I9)
where Cll and B22 are the inverses in the corresponding subspaces. Finally, the
existence of the inverse operators follows from the positive-definiteness of the
operators C and B. The Green's tensor G(k) for displacement is, in tum, ex-
pressed in terms of K or S (or D and e) with the help of simple algebraic opera-
tions. At the same time, obtaining explicit formulae in the x-representation is
possible only in the simplest cases, in particular for the isotropic or hexagonal
symmetries.
The Green's tensors and projection operators for the isotropic medium were
obtained in [B6.25]. In particular, it is shown that the projections Do and 8 0
admit an additional invariant decomposition
Do = m+ m, 8 0 = 80+8 ~, (A2.20)
where V(x) is the characteristic function of the region V. Further, taking into
account Parceval's equality (2.1.6), we find
1
A = (2n)3 SkG(k)kV(k) dk, (A3.5)
1
A = (21t)3 JK(k) V(k) dk . (A3.6)
Let us assume now that {) is an ellipsoid given by (4.9.13) and perform the
substitution x = ay. Then V'(y) = V(ay) is the characteristic function of the
unit sphere and V(k) = det{a}V'(ak). Substituting into (A3.6) and replacing
k by a-1k, we obtain
(A3.8)
and, hence,
(A3.9)
where the integration is carried out over the unit sphere Q)2 = 1.
Note that the constant A does not depend on the absolute dimensions of
the ellipsoid {), since K(k) is a homogeneous function of the zeroth degree.
Therefore the functional defined by (A3.2) does not depend on absolute dimen-
sions of the ellipsoid. This enables one to pass to the limit by contracting the
ellipsoid to a point; for definiteness, let this point be x = 0). Then the integral
over V vanishes, since K(x) - IxP, and the second integral is expressed as an
integral in the sense of principal value, which exists since (K, rp) is well defined.
Thus, finally,
where
D = C - CAC. (A3.15)
The regularization of K(x) and Sex) can be carried over to the plane case
in an obvious way, while
A = I
21t JK(tr 1l1) dll), D = I
21t JSea-ll1) dll), (A3.16)
(A3.19)
It is easily verified that the k-representation of the function T has the form
(A3.20)
where the integral over the plane is understood in the sense of principal value
and exists.
It is not difficult to verify that in the Fourier representation,
where b(x)o(D) is a o-function with a smooth weight b(x) (Sect. 6.1). The func-
tion lex) is defined for all x, other than xED, since, in the latter case, the in-
tegral (A3.23) formally diverges. Let us find a regular representation of the
integral lex) for x E D. For this purpose, we rewrite lex) in the form
From Gauss' formula and from the property div S = 0 of the function Sex),
follows the equality
This representation is valid for all x, while, for XED the existence of the integral
in the sense of a principal value follows from (A3.21).
Now let 0 in (A3.23) be a smooth, not closed, simply-connected surface
bound by a smooth contour r.This case is reduced to the previous one, if we
consider 0 to be a part of the closed surface Or. and hex) to be zero for all x E
{) = 0 1\0. In this case, the integral (A3.26) is represented by
With the help of Stokes' formula, the second integral is transformed into an
integral along the boundary r.
Here one has to take into account that, due to
(5.3.3)1
Let us consider the following averages of the functions X(x) and X(XI; x)
with respect to the ensemble of realizations of the point field X:
<X(X» = <eiEX
1: o(x - ~,.) , (A4.3)
(A4.4)
(A4.5)
Here G is the region of averaging, its limit coinciding with the whole space,
v is the volume of G, N is a number of points of a fixed realization of the point
fieId, which are contained in G, and Vo is an average volume per point.
Analogously,
to within a null-set. Note that all the equalities will be understood in this sense.
Let us proceed to the more complicated problem of constructing the two-
point moment of the function X(x). For ergodic fields
(A4.1O)
From here we find the expression for the numerator in the conditional mean
(A4.4). By virtue of the definition (A4.2) of the function X(XI; x), the equality
(A4.14)
then
(A4.15)
(X(Xl; x) I x) = -L . (A4.16)
Vo
In order to construct the conditional mean (A4.5), note that, due to the
absence of correlation between positions of the points (Poisson's field I), the
equalities
hold.
Calculating the right means analogously to the preceding and substituting
the result into (A4.5), we obtain
(A4.19)
f;m = m + pm + r , (A4.20)
(A4.21)
J
g(X) = (21t)-3 f(k)f( - k)eikx dk , (A4.23)
A4 Calculation of Certain Conditional Means 243
Let us consider the double sum in (A4.14) and proceed from two-index
quantities gmn to one-index ones, and rewrite this double sum as a single one.
The vectors m - n, with m, nEG, take on N2 values, some of which may be
identical. Each of these vectors coincides with one of the vectors of the fixed
lattice. The number of vectors, which are equal to a fixed one m, depends on
the relations between the length Iml of m and dimensions of G. If the dimen-
sions of G are much larger than Iml, then this number is '" N. Taking this
into account, (A4.14) can be represented by
where the primed summation sign means exclusion of the term with m = O.
From here and from (A4.6) follows the desired expression for the mean
(A4.4)
Let us proceed to constructing the mean (A4.5) for the random lattice.
The mean in the numerator of (A4.5) has the form
(A.4.27)
(A4.28)
244 Appendices
We substitute this result into (A4.27) and transform the double sums to
single ones. Using (A4.25), we find the expression for the mean (A4.5) in
the form
(A4.30)
where the primed summation sign means exclusion of the terms with m =
o or n = o.
We complete the consideration of random point fields and proceed to the
construction of analogous means for homogeneous random fields of noninter-
secting ellipsoids. Let us confine ourselves to the problem of constructing
the means of the characteristic functions Vex) and V(XI; x) of the form
(A4.33)
Note that JV,{x) V.(x + Xl) dx is the volume of intersection of two identical
ellipsoids, the centers of which differ by the vector Xl. From here, we have
(A4.36)
A4 Calculation of Certain Conditional Means 245
where
Starting with this expression, one can show that W(x) is a continuous
function with the mean value p. Since the different V; do not intersect W(O) =
O. For each realization it follows from the definition of the function V(x; Xl)
that V(x; x) = O.
The form of the function W(x) depends on specific structure of the random
field of inclusions. In particular, for regular lattices of identical inclusions,
the mean (V(x) V(Xl» is a periodic function of the difference X - Xl. In a
neighborhood of the origin, the form of this function is determined by the
first term in (A4.35). From here, it follows that W(x) can be represented in
the form
W(X) = p 1:'
m
J(x - m) , (A4.39)
(O(x» = lim .l
v-+oo
J1: 0; (x)dx =
V G ;
lim.l
v-+oo
t 1ta;b;,
V ;=1
(A4.40)
246 Appendices
where the aj, h j are semiaxes of the ellipse OJ. Since, for a homogeneous field,
all stochastic variables aj, h j have identical distributions with the same mean
value of <a h), then
The axes xl, x 2 coincide with the principal axes of the ellipse OJ.
Analogously to (A4.35), let us write (A4,42) in the form
<h(x)O(x)O(x + Xl»
= lim .l..
v-+ooV.
2: JG h;(x)O,{x)Oj(x + xI)dx + <~) U
?IJ'(XI) , (A4.44)
where
(A4,45)
where
A4 Calculation of Certain Conditional Means 247
(A4.48)
where m is a vector of the lattice, formed by the centers of the ellipses Dm,
the primed summation sign denoting omission of the term with m = O.
In Chap. 7, the simplified notation
was introduced.
References
Chapter 1
1.1 E. Cosserat, F. Cosserat: Theorie des Corps Deformables (Hermann, Paris
1909)
1.2 C. Truesdell, R.A. Toupin: "The Classical Field Theories", in: Principles
of Classical Mechanics and Field Theory, Encyclopedia of Physics, Vol.
III/I (Springer, Berlin, Gottingen, Heidelberg 1960) pp. 226-793
1.3 Cz. WoZniak: Dynamic models of certain bodies with discrete-con-
tinuous structure. Arch. Mech. Stosow. 21, 707 - 724 (1969)
1.4 G.N. Savin, Iu. N. Nemish: Investigations on stress concentration in
couple-stress elasticity (survey). Prikl. Mekh. 6, 1-17 (1968) [in Rus-
sian]
1.5 A.A. Il'jushin, V.A. Lomakin: "Couple-Stress Theories in the Mechanics
of Solid Deformable Bodies", in Prochnost i plastichnost (Strength and
Plasticity), ed. by A. Il'jushin (Nauka, Moscow 1971) pp. 54-61 [in
Russian]
1.6 M. Misicu: Mecanica Mediilor Deformabile. Fundamentele Eleaticitatii
Structurale (Continuum Mechanics. Foundations of Elasticity) (Acad.
R.S.R., Bucharest 1967) [in Rumanian]
1.7 W. Nowacki: Theory of Micropolar Elasticity (Wydawn. uczeln. poli-
techn., Poznan 1970) [in Polish]
W. Nowacki: Theory of Nonsymmetric Elasticity (Zaklad narod. im
Ossolinskich PAN, Wroclaw 1970) [in Polish]
1.8 E. Kroner (ed.): Mechanics of Generalized Continua. Proc IUTAM
Symposium on the Generalized Co sse rat Continuum and the Continuum
Theory of Dislocations with Applications, Freudenstadt and Stuttgart
1967 (Springer, Berlin, Heidelberg, New York 1968)
1.9 LA. Kunin, V.G. Kosilova (eds): Teoria Sred s Mikrostrukturoi (Theory
of Media with Microstructure, Bibliography) (Sib. Otdelenie AN SSSR,
Institut teplofiziki, Novosibirsk 1976)
1 .1 0 E. Kroner: On the physical reality of torque stresses in continuum
mechanics. Int. 1. Eng. Sci. 1 ,261-278 (1963)
1.11 LA. Kunin: Teoria Uprugikh Sred c Mikrostrukturoi (Theory of Elastic
Media with Microstructure) (Nauka, Moscow 1975)
l.12 LA. Kunin: Elastic Media with Microstructure lOne-Dimensional
Models, Springer Series in Solid-State Sciences, Vol. 26 (Springer,
Berlin, Heidelberg, New York 1982)
1.13 D.G.B. Edelen: Irreversible thermodynamics of nonlocal systems. Int. 1.
Eng. Sci. 12,607 -631 (1974)
250 References
Chapter 2
2.1 l.M. Gel'fand, G.E. Shilov: Generalized Functions (Academic, New York
1964)
2.2 M.Born, K. Huang: Dynamical Theory of Crystal Lattices (Clarendon
Press, Oxford 1954)
2.3 E.L. Aero, E.V. Kuvshinskii: Fundamental equations of the theory of
elastic media with rotationally interacting particles. SOy. Phys. Solid
State 2, 1272-1281 (1960)
2.4 G. Grio1i: Elasticita asimmetrica. Annali di Matematica Pura e Applicata
4,389-418 (1960)
2.5 R.D. Mindlin, H.E. Tiersten: Effects of couple-stresses in linear elasticity.
Arch. Rat. Mech. Anal. 11 ,415-448 (1962)
2.6 G. Leibfried: "Gittertheorie der mechanischen und thermischen Eigen-
schaften der Kristalle ", in Encyclopedia of Physics, Vol.7, Part 1
(Springer, Berlin 1955)
2.7 1.M. Ziman: Electrons and Phonons (University Press, Oxford 1960)
Chapter 3
3.1 V.S. Oskotskii, A.L. Efros: On the theory of crystal lattices with non-
central interatomic interaction. SOy. Phys.-Solid State 3, 448-457
(1961)
3.2 V.E.Vdovin, LA. Kunin: A theory of elasticity with spatial dispersion.
Three-dimensional complex structure. 1. Appl. Math. Mech. 30, 1272-
1281 (1966)
3.3 l.M. Gel'fand, S.V. Fomin: Calculus of Variations (Prentice-Hall, Engle-
wood Cliffs, NJ 1963)
3.4 V.E. Vdovin: On energy relations in the theory of media with micro-
structure. Dinamika Sploshnoi Sredy 7, 94-104 (1971)
3.5 M.Born, K. Huang: Dynamical Theory of Crystal Lattices (Clarendon
Press, Oxford 1954)
3.6 K. Huang: On the atomic theory of elasticity. Proc. R. Soc. A203, No.
1072, 178-194 (1950)
3.7 A.O. Gel'fond: Calcul des Differences Finies (Dunod, Paris 1963)
3.8 V.A. Pal'mov: Fundamental equations of the theory of asymmetric
elasticity. Prikl. Math. and Mech. 28, 401-408 (1964)
Chapter 4
4.1 l.M. Lifshits, L.N. Rosentsveig: Green's tensor for anisotropic unbounded
elastic medium. Zh. Exsp. Teor. Fiz.17, 783-791 (1947)
4.2 LM. Gel'fand, G.E. Shilov: Generalized Functions (Academic Press, New
York 1964)
4.3 G.l. Eskin: Boundary- Value Problems for Elliptic Pseudo-Differential
Equations (Math. Soc. of USA, N. Y. 1981)
References 251
4.4 J.R. Willis: The stress field around an elliptical crack in an anisotropic
elastic medium. Int. J. Eng. Sci. 6, 253-263 (1968)
4.5 J.D. Eshelby: "Elastic inclusions and inhomogeneities", in Prog. Solid
Mech., Vol. 2, ed. by LN. Sneddon, R. Hill (1961) pp. 88-140
4.6 J.D. Eshelby: The elastic field outside an ellipsoidal inclusion. Proc. R.
Soc. A252, 561-569 (1959)
4.7 J.D. Eshelby: The determination of the elastic field of an ellipsoidal
inclusion, and related problems. Proc. R. Soc. A241, 376-396 (1957)
4.8 I.M. Lifshifts, L.V. Tanatarov: On the elastic interaction of impurity
atoms in crystal. J. Phys. Metals Metallography 12, 331-338 (1961)
Chapter 5
5.1 E. Kroner: Kontinuumstheorie der Versetzungen und Eigenspannungen,
(Springer, Berlin, Gottingen, Heidelberg 1958)
5.2 B.A. Bilby: "Geometry and continuum mechanics", in [Ref. B. 2.1, pp.
180-199]
5.3 E. Kroner: Allgemeine Kontinuumstheorie der Versetzungen und Eigen-
spannungen, Arch. Rat. Mech. Anal. 4 , 273-334 (1960)
5.4 I.A. Kunin: "Methods of tensor analysis in the theory of dislocations".
A supplement to the Russian translation of J .A. Schouten: Tensor Ana-
lysis for Physicists (Nauka, Moscow 1965), pp. 374-443. (The supple-
ment is available in English from the U.S. Department of Commerce,
Clearing house for Federal Sci. and Techn. Information, Springfield,
VA 2215})
5.5 RAA G Memoirs of the Unifying Study of Basic Problems in Engineering
and Physical Sciences by Means of Geometry, ed. by K. Kondo (Gaku-
jutsu Bunken Fukyu-Kai, Tokyo, V.I. 1955, V. II 1958; V. III 1962,
V. IV 1968)
5.6 C.c. Wang: "On the geometric structure of simple bodies, a mathematical
fo'undation for the theories of continuous distributions of dislocations",
in [Ref. B2.1, pp. 247-250J
5.7 R. de Wit: "Differential geometry of a nonlinear continuum theory of dis-
locations", in [Ref. B2.1, pp. 251-261]
5.8 I.A. Kunin: Internal stresses in media with microstructure. J. Appl.
Math. Mech. 31, 898-906 (1967)
5.9 LA. Kunin: An algebra of tensor operators and its applications to elasti-
city. Int. J. Eng. Sci., 19, 1551-1561 (1981)
5.10 A.I. Lur'e: Three-dimensional Problems of the Theory of Elasticity (In-
terscience, New York 1964)
5.11 I.M. Lifshits, L.V. Tanatarov: On the elastic interaction of impurity
atoms in crystals. J. Phys. of Metals and Metallography 12, 331-338
(1961)
252 References
Chapter 6
6.1 J. Friedel: LesDislocations (Gauthier-Villars, Paris 1956)
6.2 E. Kroner: Kontinuumstheorie der Versetzungen und Eigenspannungen,
(Springer, Berlin, Gottingen, Heidelberg 1958)
6.3 LA. Kunin: Fields due to arbitrary distributions of dislocations in an ani-
sotropic elastic medium. J. Appl. Mech. Tech. Phys. 6, 76 (1955)
6.4 LA. Kunin: "Methods of tensor analysis in the theory of dislocations". A
supplement to the Russian translation of J.A. Schouten. Tensor Analysis
for Physicists (Nauka, Moscow 1965), pp. 374-443. (The supplement is
available in English from the U.S. Department of Commerce, Clearing
house for Federal Sci. and Techn. Information, Springfield, VA 22151)
6.5 LM. Gel'fand, G.E. Shilov: Generalized Functions (Academic Press, New
York 1964)
6.6 A.D. Brailsford: Stress field of a dislocation. Phys. Rev. 142,383-392
(1966)
Chapter 7
7.l 1.1. Gikhman, A.V. Skorokhod: Theory of Random Processes, Vol. 1,
(Nauka, Moscow 1971) [in Russian]
7.2 T.D. Shermergor: Theory of Elastic Inhomogeneous Media (Nauka,
Moskow 1977) [in Russian]
7.3 LA. Chaban: Self-consistent field approach to calculation of effective
parameters of microinhomogeneous media. Sov. Phys.-Acoust. 10, 298-
304 (1964)
7.4 E. Kroner: Bounds for effective elastic moduli of disordered materials. J.
Mech. and Phys. Solids, 25, 137 -15 5 (1977)
7.5 1.1. Gikhman, A.V. Skorokhod: Theory of Random Processes, Vol.2
(Nauka, Moscow 1971) [in Russian]
7.6 V.M. Levin: On the stress concentration in inclusions in composite
materials. J. Appl. Math. and Mech. 41, 753-761 (1977)
7.7 A.V. Hershey: The elasticity of an isotropic aggregate of anisotropic
cubic crystals. J. Appl. Mech. 21,236 (1954)
7.8 G. Kneer: Dber die Berechnung der Elastizitatmoduls vielkristaller Ag-
gregate mit Textur. Phys. Stat. Sol. 9, 825--838 (1965)
7.9 E. Kroner: Berechnung der elastischen Konstanten des Vielkristalls
aus den Konstanten des Einkristalls, Z. Phys.151, 504 (1958)
7.l 0 P.R. Morris: Elastic constants of polycrystals. Int. J. Eng. Sci. 8, 49
(1970)
7.ll Z. Hashin, S. Shtrikman: On some variational principles in anisotropic
and nonhomogeneous elasticity; J. Mech. Phys. Solids 10, 343 (1962)
7.12 S.K. Kanaun: Self-consistent field approximation for an elastic compo-
site medium. J. Appl. Mech. Tech. Phys. 18, 274-282 (1977)
7.13 E.L Grigoliuk, L.A. Fil'shtinskii: Perforated Plates and Shells (Nauka,
Moscow 1970) [in Russian]
References 253
7.14 A.E. Andreikiv, V.V. Panasiuk, M.M. Stadnik: Fracture of prismatic bars
with cracks. Problemy Prochnosti, No. 10. 10-16 (1972)
7.15 B. Budjanski, R.J. O'Connel: Elastic moduli of a cracked solid. Int. J.
Solids Struct. 12, (1976)
7.16 A.S. Vavakin, R.L. Salganik: On effective characteristics of inhomogenous
media with local inhomogeneities. Izv. Akad. Nauk SSSR, MTT No.3
(1975)
7.17 S.K. Kanaun, G.I. Iablokova: "Self-consistent field approximation in
plane problem for systems of interacting cracks", in Mekhanika Sterzhne-
vykh Sistem i Sploshnykh Sred, No.9, (Leningr. Inshenemostroitel'nyi
Institut, Leningrad 1976) pp. 194-203
7.18 L.A. Fil 'shtinskii: Interaction of a doubly-periodic system of rectilinear
cracks in an isotropic medium. J. Appl. Math. and Mech. 38 ,853-861
(1974)
7.19 P.e. Daris, G.C.M. Sih: "Stress analysis of cracks" in Fracture Toughness
Testing and its Applications (Am. Soc. for Testing and Materials,
1965) pp. 30-83
7.20 V.V.Panasiuk, M. Savruk, A.P. Datsishin: Stress Distribution Around
Cracks in Plates and Shells (Naukova Dumka, Kiev 1976)
7.21 G. Adomian: Linear random operator equations in mathematical physics.
J. Math. Phys. 12, 1944 (1971)
7.22 R.J. Elliott, J.A. Krumhansl, P.L. Leath: The theory and properties of
randomly disordered crystals. Rev. Mod. Phys. 46,465 (1974)
7.23 M. Hori: Statistical theory of effective electrical, thermal and magnetic
properties of random heterogeneous materials. J. Math. Phys. 14, 514-
523; 1942-1948 (1973); 15,2177-2185 (1974); 16,352-364; 1772-
1775 (1975); 18,487-501 (1977)
7.24 A.H. Nayfeh: Pertubation Methods (Wiley, New York 1970)
7.25 V.M. Levin: On calculation of effective moduli of composites. Dokl.
Akad. Nauk SSSR 220, No.5 (1975)
7.26 LT. Walpole: On bounds for the overall elastic moduli of inhomoge-
neous systems - I, II. J. Mech. Phys. Solids 14, 151 and 289 (1966)
7.27 M. Kac: "Mathematical mechanisms of phase transitions" in Statistical
Physics. Phase Transitions and Superf!uidity, Vol. 1 ed. by M. Chretien,
(Gordon & Breach, New York 1968) pp. 241-305
Bibliography
B2 Couple-Stress Theories
B2.l E. Kroner (ed.): Mechanics of Generalized Continua. Proc. IUTAM
Symposium on the Generalized Cosserat Continuum and the Continuum
Theory of Dislocations with Applications, Freudenstadt and Stuttgart
1967 (Springer, Berlin, Heidelberg, New York 1968)
B2.2 I.A. Kunin, V.G. Kosilova (eds.): Theory of Media with Microstructure.
Bibliography (Sib. Otdelenie AN SSSR, Institut Teplofiziki, Novosibirsk
1976)
B2.3 E.L. Aero, E.V. Kuvshinskii: Continuum theory of asymmetric elastici-
ty. The problem of internal rotation. SOy. Phys.-Solid State 5, 1892-
1897 (1963)
B2.4 E.L. Aero, E.V. Kuvshinskii: Continuum theory of asymmetric elastici-
ty. Equilibrium of an isotropic body. SOy. Phys.-Solid State 6,2141-
2148 (1964)
B2.S E.L. Aero, E.V. Kuvshinskii: Fundamental equations of the theory of
elastic media with rotationally interacting particles. SOy. Phys.-Solid
State 2, 1272-1281 (1960)
B2.6 A. Askar: Molecular crystals and the polar theories of the continua. Ex-
perimental values of material coefficients for KN0 3 . Int. J. Eng. Sci.
10,233-300 (1972)
256 Bibliography
B2.29 D. Rogula: Moment stresses and the symmetry of the stress tensor in
bodies with non-local structure. Bull Acad. Pol. Sci. Ser. Sci. Techn.18,
159-164 (1970)
B2.30 G.N. Savin, lu. N. Nemish: Investigations on stress concentration in
couple-stress elasticity (survey). Prikl. Mekhanika 6,1-17 (1968)
B2.31 H. Schaefer: The Cosserat continuum. Z. Angew. Math. und Mech. 47,
485-498 (1967)
B2.32 la. M. Shiriaev: Investigation of the influence of a scale factor on stress
concentrations around holes. Mekhanika Polimerov, No.3, 565 (Akad.
Nauk Latv. SSR, 1970)
B2.33 E. Sternberg: "Couple stresses and singular stress concentrations in
elastic solids", in [Ref. B2.1, pp. 95-108]
B2.34 R. Stojanovic: "Dislocations in the generalized elastic Cosserat con-
tinuum", in [Ref. B2.1, pp. 152-155]
B2.35 R.A. Toupin: Theories of elasticity with couple-stresses. Arch. Rat.
Mech. Anal. 17, 85-112 (1964)
B2.36 R.A. Toupin: "Dislocated and oriented media", in [Ref. B2.1, pp.
126-140]
B2.37 K. Wilmanski, Cz. Wozniak: On geometry of continuous medium with
microstructure. Arch. Mech. Stosowanej 19, 715 - 723 (1967)
B2.38 Cz. Wozniak: Dynamic models of certain bodies with discrete-contin-
uous structure. Arch. Mech. Stosowanej 21, 707-724 (1969)
B3 Nonlocal Theories
B3.1 D.G.B. Edelen, A.E. Green, N. Laws: Nonlocal continuum mechanics.
Arch. Rat. Mech. Anal. 43, 36-44 (1971)
B3.2 A.C. Eringen, D.G.B. Edelen: On non-local elasticity. Int. J. Engng. Sci.
10,233-248 (1972)
B3.3 A.C. Eringen: Nonlocal polar elastic continua. Int. J. Eng. ScilO, 1-6
(1972)
B3.4 A.C. Eringen: Linear theory of nonlocal microelasticity and dispersion
of plane waves. Lett. Appl. Sci. 1, 129-146 (1973)
B3.5 A.C. Eringen, C.G. Speziale, B.S. Kim: Crack-tip problem in nonlocal
elasticity. J. Mech. Phys. Solid. 25, 339-355 (1977)
B3.6 A.C. Eringen: "Nonlocal continuum mechanics and some applications",
in Nonlinear Equations in Physics and Mathematics, ed. by A.O. Barut
(Reidel, Dordrecht, Holland 1978) pp. 271-318
B3.7 F. Hehl, E. Kroner: Zum Materialgesetz eines elastichen Mediums mit
Momentenspannungen. Z. Naturforsch. 20, 336-350 (1965)
B3.8 S. Kaliski, C. Rymarz: Surface waves in continua with nonlocal interac-
tion. Biul. WAT J. Dabrowskiego 20,17-39 and 25-37 (1971)
B3.9 D. Kessel, E. Kroner: Nichtlokale Elastizitatstheorie der lonenkristalle
vom Typ A+Ir. Z. Naturforsch. 25. 1046-1053 (1970)
B3.1O V.G. Kosilova, I.A. Kunin: Dynamics of generalized Cosserat models.
Dinamika Sploshnoi Sredy No.4, 73-82 (1970)
258 Bibliography
B3.11 V.G. Kosilova, LA. Kunin: On a scattering problem for media with
dispersion. Chislennye Metody Mekhaniki Sploshnoi Sredy, 5, No.2,
51~72 (1974)
B3.12 V.G. Kosilova, LA. Kunin: Wave packets propagation in one-dimen-
sional medium with space dispersion. Dinamika Sploshnoi Sredy, No.6,
12'7~142 (1970)
B3.13 E. Kroner: On the physical reality of torque stresses in continuum
mechanics. Int. J. Eng. Sci. 1, 261~278 (1963)
B3.14 E. Kroner: Elasticity theory of materials with long range cohesive
forces. Int. J. Structures, 3, 731 ~ 742 (1967)
B3.15 E. Kroner: "Interrelations between various branches of continuum
mechanics", in [Ref. B2.1, pp. 330~340]
B3.16 E. Kroner: " Ober die Symmetrie-Eigenschaften elastischer Materialien
mit Kohasionskraften endlicher Reichweite," in "Problemy Gidrodina-
miki i Mekhaniki Sploshnoi Sredy" (Nauka, Moskva 1969) pp. 293~300
B3.17 E. Kroner, B.K. Datta: Nichtlokale Elastostatik: Ableitung aus der
Gittertheorie. Z. Physik 196. 203~211 (1966)
B3.18 E. Kroner, B.K. Datta: "Non-local theory of elasticity for a finite
inhomogeneous medium - a derivation from lattice theory", in Funda-
mental Aspects of Dislocation Theory, Vol. 2, ed. by J. Simmons and
R. de Wit (National Bureau of Standards, Washington 1970) pp. 737 ~
746
B3.19 J.A. Krumhansl: "Some considerations of the relation between solid
state physics and generalized continuum mechanics", in [Ref. B2.l,
pp. 298~311]
B3.20 LA. Kunin: Inhomogeneous elastic medium with nonlocal interaction.
J. App!. Mech. Tech. Phys. 8. 41 ~44 (1967)
B3.21 I.A. Kunin: Internal stresses in media with microstructure. J. Appl.
Math. Mech. 31, 898~906 (1967)
B3.22 LA. Kunin: Model of an elastic medium of simple structure with three-
dimensional dispersion. J. App!. Math. Mech. 30. 642~652 (1966)
B3.23 LA. Kunin: Theory of elasticity with spatial dispersion. One-dimen-
sional complex structure. J. App!. Math. Mech. 30, I 025~1 034 (1966)
B3.24 LA. Kunin: The Green's tensor for an elastic isotropic medium with
spatial dispersion. J. App!. Mech. Tech. Phys. 8, 103 (1967)
B3.25 LA. Kunin: Theories of elastic media with microstructure. Proc. Vibr.
Probl. 9, 323~336 (1968)
B3.26 LA. Kunin: "The theories of elastic media with microstructure and the
theory of dislocations", in [Ref. B2.1 , pp. 321 ~3 29]
B3.27 LA. Kunin: "The theory of elastic media with microstructure", in
Prochnost'i Plastichnost' (Nauka, Moskva 1971) pp. 65~ 70
B3.28 LA. Kunin, A.M. Vaisman: "On problems of the non-local theory of
elasticity", in: Fundamental Aspects of Dislocation Theory, Vo!. 2,
ed. by J. Simmons and R. de Wit (National Bureau of Standards,
Washington 1970) pp. 747~759
Bibliography 259
B4.51 H.C. Morris: A prolongation structure for the AKNS system and its
generalization. J. Math. Phys. 18, 533-536 (1977)
B4.52 H.C. Morris: Prolongation structures and a generalized inverse scat-
tering problem. J. Math. Phys.17, 1867-1869 (1976)
B4.53 S.P. Novikov: The periodic problem for the Korteweg-de Vries equa-
tion. Funct. Anal. Appl. 8,236-246 (1974)
B4.54 N. Ooyama, N. Saito: On the stability of lattice solitons. Supplement
of the Progress of Theoretical Physics No. 45,201-208 (1970)
B4.55 A.C. Scott, F.Y. Chu, D.W. Laughlin: The soliton: a new concept in
applied science. Proc. IEEE 61, 1443-1483 (1973)
B4.56 A.B. Shabat: On non-scattering potentials in the Sturm-Liouville
equation. Dinamika Sploshnoi Sredy, No.5, 130 (1970)
B4.57 A.B. Shabat: On the Korteweg-de Vries equation. Sov. Math. Dokl.
14, 1266-1270 (1973)
B4.58 Solitons, ed. by R. Bullough, P. Caudrey, Topics Current Phys. Vol. 17
(Springer, Berlin, Heidelberg, New York 1980)
B4.59 L.A. Takhtadzhyan: Exact theory of propagation of ultrashort optical
pulses in two-level media. Sov. Phys.-JETP 39, 228-233 (1974)
B4.60 LA. Takhtadzhyan, L.D. Faddeev: Essentially nonlinear one-dimen-
sional model of classical field theory. Theor. Math. Phys. 21, 1046-
1057 (1974)
B4.61 S. Tanaka: On the N-tuple wave solutions of the Korteweg-de Vries
equation. Publ. Res. Inst. Math. Sci. Kyoto University 8, 419-427
(1972/1973)
B4.62 Theory of Solitons, ed. by S.P. Novikov (Nauka, Moscow 1980) [in
Russian]
B4.63 M. Toda: Waves in nonlinear lattice. Supplement of the Progress of
Theoretical Physics No. 45,174-200 (1970)
B4.64 M. Toda: Theory of Nonlinear Lattices, Springer Ser. Solid-State Sci.,
Vol. 20 (Springer, Berlin, Heidelberg, New York 1981)
B4.65 M. Wadafi, H. Sanuki, K. Kondo: Relationships among inverse methods,
Blicklung transformation and an infinite number of conservation laws.
Progr. Theor. Phys. 53,419-436 (1975)
B4.66 M. Wadati, M. Toda: The exact N-solution of the Korteweg-de Vries
equation. J. Phys. Soc. Jpn. 32 (1972)
B4.67 H.D. Wahlquist, F.B. Estabrook: Prolongation structures of nonlinear
evolution equations. J. Math. Phys. 16, 1-7 (1975); 17, 1293-1297
(1976)
B4.68 Y. Watanabe: Note on gauge invariance and conservation laws for a
class of nonlinear partial differential equations. J. Math. Phys. 15, 453-
457 (1974)
B4.69 N.J. Zabusky: Solitons and bound states of the time-independent
SchrOdinger equation. Phys. Rev. 168, 124-128 (1968)
B4.70 V.E. Zakharov: Kinetic equation for solitons. Sov. Phys.-JETP 33,
538-541 (1971)
Bibliography 263
BS Crystal Lattice
B5.l V.M. Agranovich, V.L. Ginzburg: Crystal-Optics with Spatial Disper-
sion and the Theory of Excitons. Springer Ser. Solid-State Sci., Vol. 42
(Springer, Berlin, Heidelberg, New York 1983)
B5.2 M. Born, K. Huang: Dynamical Theory of Crystal Lattices (Clarendon
Press, Oxford 1954)
B5.3 L. Brillouin, M. Parodi: Propagation des On des dans les Milieux Perio-
diques (Masson, Paris 1956)
B5.4 K. Huang: On the atomic theory of elasticity. Proc. R. Soc. A203,
No. 1072, 178-194 (1950)
B5.5 P.N. Keating: On the sufficiency of the Born-Huang relations. Phys.
Lett. A25, 496-497 (1967)
B5.6 A.M. Kosevich: Fundamentals of Crystal Lattice Mechanics (Nauka,
Moscow 1972)
B5.7 G. Leibfried: "Gittertheorie der mechanischen und thermischen Eigen-
schaften der Kristalle", in Encyclopedia of Physics, Vol. 7, Part 1
(Springer, Berlin 1955)
B5.8 A.N. Luzin: Some problems in the dynamics of crystal lattices and of
the theory of elasticity. SOY. Phys.-JETP 24,615-620 (1966)
264 Bibliography
Tokyo
Systems. - Coherent Protential Approxima-
tion (CPA). - References. - Subject Index.