Non-Metricity and the Nonlinear Mechanics of Distributed Point Defects

Transcription

1 Non-Metricity and the Nonlinear Mechanics of Distributed Point Defects Arash Yavari Alain Goriely May 24 Abstract We discuss the relevance of non-metricity in a metric-affine manifold (a manifold equipped with a connection and a metric) and the nonlinear mechanics of distributed point defects. We describe a geometric framework in which one can calculate analytically the residual stress field of nonlinear elastic solids with distributed point defects. In particular, we use Cartan s machinery of moving frames and construct the material manifold of a finite ball with a spherically-symmetric distribution of point defects. We then calculate the residual stress field when the ball is made of an arbitrary incompressible isotropic solid. We will show that an isotropic distribution of point defects cannot be represented by a distribution of purely dilatational eigenstrains. However, it can be represented by a distribution of radial eigenstrains. We also discuss an analogy between the residual stress field and the gravitational field of a spherical mass. Contents Introduction 2 Non-iemannian Geometries and Anelasticity 2 3 Point Defects in an Incompressible Isotropic Ball 4 3. Construction of the flat Weyl material manifold Calculation of the residual stress field Incompressible isotropic solids Compressible isotropic solids An analogy between the point defect metric and the Schwarzschild metric Singularity in the material metric in the case of a single point defect Exterior residual stress field of a ball of point defects Isotropic distribution of point defects and pure dilatational eigenstrains Conclusions Introduction The first mathematical study of line defects in solids goes back to the work of Volterra [97] more than a century ago. The close connection between the mechanics of solids with distributed defects and non-iemannian geometries was independently discovered in the 95s by Kondo [955a,b] and Bilby, et al. [955]; Bilby and Smith [956]. Defects influence many of the mechanical properties of solids and have been the focus of intense research in the last few decades. Motivated by applications of metals in industry and the need to take into To appear in. J. Knops, G.-Q. Chen, and M. Grinfeld (Editors), Differential Geometry and Continuum Mechanics, Springer Proceedings in Mathematics & Statistics (POMS). School of Civil and Environmental Engineering & The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA 3332, USA. arash.yavari@ce.gatech.edu. OCCAM, Mathematical Institute, University of Oxford, Oxford, OX 3LB, UK.

2 2 Non-iemannian Geometries and Anelasticity 2 account plastic deformations, the micro mechanism of plasticity, i.e., dislocations have been studied by many researchers but mostly in the framework of linearized elasticity. Other line defects, e.g. disclinations have also been the subject of many investigations. However, point defects have not received much attention even in the linearized setting after the original works of Love [927] and Eshelby [954]. In particular, Love [927] calculated the stress field of a single point defect in an infinite linear elastic solid and observed a /r 3 singularity. ecently, we [Yavari and Goriely, 22b] showed how one can calculate the residual stress field of a nonlinear elastic solid with distributed point defects. Non-metricity proved to be an essential geometric entity in describing the zero-stress configuration (material manifold) of solids with point defects. Here, material manifold (B, G, ) is a flat, torsion-free 3-manifold with non-metricity in which the body is stress free. We should mention that many researches have known the relevance of non-metricity to the mechanics of point defects [Falk, 98; de Wit, 98; Grachev, et al., 989; Kröner, 99; Miri and ivier, 22]. However, there has not been a concrete use of non-metricity in the literature for calculating residual stresses. The geometric framework discussed here has been recently used by the authors in the analysis of distributed dislocations and disclinations as well [Yavari and Goriely, 22a, 23b]. In this book chapter we review the results of [Yavari and Goriely, 22b], extend the residual stress calculation to arbitrary incompressible isotropic solids, and make some new observations. In particular, we discuss an analogy with relativity and the gravitational field of a spherical ball of mass m in an infinite empty space-time. Another problem that can benefit from geometric ideas is the stress analysis of inclusions in nonlinear elastic solids [Yavari and Goriely, 23a]. In [Yavari and Goriely, 23a] we showed that collapsing a small spherical inclusion with pure dilatational eigenstrain while keeping the strength of the inclusion fixed one recovers the stress field of a single point defect in a linear elastic solid. Earlier in [Yavari and Goriely, 22b] we had shown that using the nonlinear solution one can recover the classical linear solution for small strength point defect distributions supported in a small ball. Now one may be tempted to think that any isotropic distribution of point defects can be represented by pure dilatational eigenstrains. We will show that this is not the case and that material metric calculated in [Yavari and Goriely, 22b] is equivalent to a distribution of radial eigenstrains with no circumferential eigenstrains. 2 Non-iemannian Geometries and Anelasticity To make this book chapter self contained, in the following we tersely review the necessary geometric background. iemann-cartan manifolds. On a manifold B a linear (affine) connection is an operation X (B) X (B) X (B), where X (B) denotes the set of vector fields on B. In a local coordinate chart {X A }, A B = Γ C AB C, where Γ C AB are Christoffel symbols of the connection and A = are the natural bases for the tangent x A space. is compatible with a metric G of the manifold if G =. The torsion of is defined by T (X, Y) = X Y Y X [X, Y], where [.,.] is the commutator of vector fields. is symmetric if it is torsion-free, i.e. X Y Y X = [X, Y]. In (B,, G) the curvature is a map X (B) X (B) X (B) X (B) defined by (X, Y)Z = X Y Z Y X Z [X,Y] Z. A iemann-cartan manifold (B,, G) is a metric-affine manifold in which the metric and the connection are compatible. Cartan s moving frames. Let us consider a frame field {e α } N α= that forms a basis for the tangent space of B everywhere. We assume that this frame is orthonormal, i.e. e α, e β G = δ αβ. {e α } N α= is, in general, a noncoordinate basis for the tangent space, i.e. it is not necessarily induced from a coordinate chart. The frame field {e α } naturally defines the co-frame field {ϑ α } N α= such that ϑ α (e β ) = δβ α. The connection -forms are defined by e α = e γ ω γ α. The corresponding connection coefficients are defined as eβ e α = ω γ α, e β e γ = ω γ βαe γ, i.e. ω γ α = ω γ βαϑ β. Similarly, ϑ α = ω α γϑ γ, and eβ ϑ α = ω α βγϑ γ. In an orthonormal frame, the metric has the simple representation G = δ αβ ϑ α ϑ β. Non-metricity. Given a metric-affine manifold (B,, G), the non-metricity is a map Q X (B) X (B) X (B) defined as Q(U, V, W) = U V, W G + V, U W G U[ V, W G ]. In other words, Q = G. In the frame {e α }, Q γαβ = Q(e γ, e α, e β ). Non-metricity -forms are defined by Q αβ = Q γαβ ϑ γ. One can show that In a metric-affine manifold the torsion, curvature, and non-metricity, are, in general, non-vanishing.

3 2 Non-iemannian Geometries and Anelasticity 3 Q γαβ = ω βγα + ω αγβ dg αβ, e γ, where d is the exterior derivative. Thus, Q αβ = ω αβ + ω βα dg αβ = DG αβ, where D is the covariant exterior derivative. This is Cartan s zeroth structural equation. For an orthonormal frame G αβ = δ αβ and hence Q αβ = ω αβ + ω βα. In a metric-affine manifold with non-metricity, the Weyl -form is defined as Q = n Q αβg αβ. Therefore, Q αβ = Q αβ + QG αβ, where Q is the traceless part of the non-metricity. If Q = and if is torsion-free, (B,, G) is called a Weyl manifold. The torsion and curvature 2-forms are defined by These are, respectively, Cartan s first and second structural equations. T α = dϑ α + ω α β ϑ β, (2.) α β = dω α β + ω α γ ω γ β. (2.2) The compatible volume element on a Weyl manifold. A volume element on B is any non-vanishing n-form. In the orthonormal coframe field {ϑ α } the volume form is written as µ = hϑ... ϑ n, for some positive function h. In a coordinate chart {X A } we have µ = h det G dx... dx n. Divergence of an arbitrary vector field W on B is defined as (Div W)µ = L W µ, where L is the Lie derivative. Having a connection divergence can also be defined as the trace of the covariant derivative, i.e. Div W = W A A = W A,A + Γ A ABW B. The volume element µ is compatible with if L W µ = (W A A)µ, which is equivalent to D (h det G) = [Saa, 995]. Thus, dh h = d ln h = n Q. Note that this implies that dq =. Therefore, to be able to define a compatible volume 2 element the Weyl one-form must be closed. Geometric anelasticity. Let us briefly review geometric nonlinear elasticity. We identify a body B with a iemannian manifold B and a configuration of B is a mapping ϕ B S, where S is another iemannian manifold (B, G). It is assumed that the body is stress free in the material manifold. The deformation gradient is the tangent map of ϕ and is denoted by F = T ϕ. At every point X B, F is a linear map F(X) T X B T ϕ(x) S. Choosing local coordinate charts {x a } and {X A } on S and B, respectively, the components of F read F a A(X) = ϕa (X). (2.3) XA The transpose of F is defined by F T T x S T X B, FV, v g = V, F T v G, for all V T X B, v T x S. In components (F T (X)) A a = g ab (x)f b B(X)G AB (X), where g and G are metric tensors on S and B, respectively. The right Cauchy-Green deformation tensor C(X) T X B T X B is defined by C(X) = F(X) T F(X) = (g ab ϕ)f a AF b B. The relation between the iemannian volume element dv at X B and its corresponding deformed volume element at x = ϕ(x) S is dv = JdV, where J = det g/ det G det F is the Jacobian. The left Cauchy-Green deformation tensor is defined as B = ϕ (g ) with components B AB = (F ) A a(f ) B b g ab. The spatial analogues of C and B are c = ϕ (G), c ab = (F ) A a (F ) B b G AB, (2.4) b = ϕ (G ), b ab = F a AF b BG AB, (2.5) where ϕ and ϕ are the pull-back and push-forward by ϕ, respectively. b is called the Finger deformation tensor. Note that C and b have the same principal invariants denoted by I, I 2, and I 3 [Ogden, 984]. For an isotropic material the strain energy function W depends only on the principal invariants of b. One can show that for a compressible and isotropic material the Cauchy stress has the following representation [Doyle and Ericksen, 956; Simo and Marsden, 983] σ = 2 ( I 2 I 3 W I 2 + W I 3 ) g + 2 W I b 2 W I 2 b. (2.6) Similarly, for an incompressible and isotropic material the Cauchy stress has the following representation [Doyle

4 3 Point Defects in an Incompressible Isotropic Ball 4 and Ericksen, 956; Simo and Marsden, 983] σ = ( p + 2I 2 W I 2 ) g + 2 W I b 2 W I 2 b. (2.7) Material manifold. We assume that the defect-free body is stress-free in Euclidean space in the absence of external loads. This body may be made of a material with multiple stress-free configurations (corresponding to multiple wells of a strain-energy density) but we assume that all the material points are in the same energy well (that is we are not considering phase transformations). Next, we assume that some distribution of defects appears in this body and induces residual stresses (we are not considering nucleation or the work associated with the creation of defects). The stress-free configuration with defects (referred to as the material manifold) explicitly depends on the distribution of defects and their types but not on the constitutive equations. Of course, residual stresses explicitly depend on the choice of the constitutive equations. Point defects. The classical continuum picture of a single vacancy is the following. emove a small spherical ball from the body and identify all the points on its boundary sphere. A single interstitial or extra matter can be visualized by inserting a larger elastic ball inside the spherical cavity and letting the system relax. Consider a distribution of point defects or extra matter in a solid. If one imagines partitioning this body into a large number of small cubes of the same size and let them relax the resulting stress-free cubes will have different sizes (and hence volumes) depending on the distribution of point defects. The relaxed volumes are the local embedding of the underlying iemannian material manifold into the ambient (Euclidean) space. In other words, in the stress-free configuration of the defective-solid volume elements vary depending on the distribution of point defects. It is known that in a metric-affine manifold (B,, G) with non-metricity the iemannian volume element is not covariantly constant. More specifically, D det G = d det G ω α α det G = n 2 Q det G. This shows that the Weyl one-form Q is somehow related to the volume density of point defects. In the following example, we will use a semi-inverse method and start with a coframe field with some unknown function(s). This then implies that the material metric is known up to the unknown function(s). To relate this unknown function(s) to the volume density of point defects we will find a compatible volume element on the material manifold B, i.e., a volume element that is covariantly constant. 3 Point Defects in an Incompressible Isotropic Ball In [Yavari and Goriely, 22b] we considered a ball of radius o with a spherically-symmetric isotropic distribution of point defects. We constructed the material manifold and calculated the residual stress field for an incompressible neo-hookean solid. Here we first construct the material manifold and then calculate the residual stress field when the defective ball is made of an arbitrary incompressible isotropic solid. We also consider a special class of compressible solids for a particular example of distributed point defects. 3. Construction of the flat Weyl material manifold In a body with only point defects the material manifold is a flat Weyl manifold, i.e. the torsion and the curvature of the material connection both vanish. In order to find a solution, we start by an ansatz for the material coframe field. In the spherical coordinates (, Θ, Φ),, Θ π, Φ < 2π, we assume the following coframe field ϑ = f()d, ϑ 2 = dθ, ϑ 3 = sin Θ dφ, (3.) for some unknown function f() > to be determined. Assuming that the non-metricity is traceless and isotropic, i.e. Q αβ = 2δ αβ q()ϑ, the matrix of connection -forms has the following form ω = [ω α β] = ω ω 2 ω 3 ω 2 ω 2 2 ω 2 3, (3.2) ω 3 ω 2 3 ω 3 3

5 3.2 Calculation of the residual stress field 5 where ω = ω 2 2 = ω 3 3 = q()ϑ, for a function q() to be determined. We now need to enforce T α =. Note that dϑ =, dϑ 2 = f() ϑ ϑ 2, dϑ 3 = f() ϑ3 ϑ + cot Θ ϑ2 ϑ 3. (3.3) Cartan s first structural equations read This implies that T = ω 2 ϑ 2 ω 3 ϑ 3 =, (3.4) T 2 = [ f() + q()] ϑ ϑ 2 ω 2 ϑ + ω 2 3 ϑ 3 =, (3.5) T 3 = cot Θ ϑ2 ϑ 3 [ f() + q()] ϑ3 ϑ + ω 3 ϑ ω 2 3 ϑ 2 =. (3.6) ω 2 = [ f() + q()] ϑ2, ω 2 3 = cot Θ ϑ3, ω 3 = [ f() + q()] ϑ3. (3.7) It can be shown that = 2 2 = 3 3 = are trivially satisfied. The remaining Cartan s second structural equations read The first equation gives us the following ODE The solution is 2 = 2 = dω 2 + ω 3 ω 2 3 =, (3.8) 2 3 = 3 2 = dω ω 2 ω 3 =, (3.9) 3 = 3 = dω 3 + ω 2 3 ω 2 =. (3.) d f() d [ f() + q()] + f() [ + q()] =. (3.) f() Note that when q() =, we have f() = and hence C =. Therefore In this example the Weyl -form is written as f() + q() = C. (3.2) q() = [ f() ]. (3.3) Q = 2q()ϑ = 2 [ f() ] ϑ = 2(f() ) d. (3.4) The function f() is determined using the volume density of point defects n(x) and using the equation µ µ = nµ [Yavari and Goriely, 22b]. In the particular example of a defective ball µ = 2 sin Φd dθ dφ and µ = f()h()µ, and hence n() f() = 3 3y2 n(y)dy. (3.5) 3.2 Calculation of the residual stress field In this section we extend our previous calculation in [Yavari and Goriely, 22b] to arbitrary incompressible isotropic solids and a certain class of compressible isotropic solids. We consider a ball of radius o and assume that a point defect density n() is given.

6 3.2 Calculation of the residual stress field Incompressible isotropic solids The material metric has the following form G = f 2 () 2 2 sin 2 Θ, G = f 2 () 2 2 sin 2 Θ. (3.6) Having the underlying iemannian material manifold, we obtain the residual stress field by embedding it into the Euclidean ambient space, which is the Euclidean 3-space. We look for solutions of the form (r, φ, z) = (r(), Φ, Z), and hence det F = r (). Assuming an incompressible solid we have J = det g det G det F = r 2 f() 2 r () =. (3.7) Assuming that r() =, we obtain r() = ( 3x 2 f(x)dx) 3. (3.8) The Finger tensor b (b ab = F a AF b BG AB ) is found to be The principal invariants of b read b = 4 r 4 () 2 2 sin 2 Θ. (3.9) Now (b ) ab = c ab = g am g bm c mn has the following representation I = 2 r2 () r 4 (), I 2 = 2 2 r 2 () + r4 () 4. (3.2) b = r 4 () 4 2 r 4 () 2 r 4 () sin 2 Θ. (3.2) Therefore, the Cauchy stress can be written as σ = p + α 4 r 4 where α = 2 W I and β = 2 W I 2. The only non-trivial equilibrium equation reads + 2β 2 r 2 p + α + β ( 2 + r2 ) r 2 2 r 4 4 sin 2 Θ [ p r 2 + α 2 + β ( 2 r 4 + r2 4 )], (3.22) σ rr,r + 2 r σrr rσ θθ r sin 2 θ σ φφ =. (3.23) Or σ rr, + 2f2 r 2 ( r σrr rσ θθ ) =. (3.24)

7 3.2 Calculation of the residual stress field 7 This gives us the differential equation p () = k(), where k() = α () 4 r 4 () + 2β () (β() + α() 2 r 2 () r 2 () ) [ 2 r 2 () f() 4 r 5 () f()r() ]. (3.25) 2 Suppose that at the boundary σ rr ( o ) = p. Thus 4 o 2 o p() = p + α( o ) r 4 ( o ) + 2β( o) r 2 ( o ) Once p() is known all the stress components are easily calculated from (3.22). o k(ξ)dξ. (3.26) Example 3.. In [Yavari and Goriely, 22b] we considered the following point defect distribution n i, n() = i <, (3.27) where i < o. Thus, one can see that This yields i f() =, (3.28) > i f() =. n ( i )3 (3.29) i r() =, (3.3) > i r() = [ 3 + n 3 i ln ( (/ i) 3 n n )] 3. (3.3) Note that for < i, λ = λ 2 = λ 3 = and hence I = I 2 = 3. Therefore, α = α and β = β are constants and consequently for < i, k() =. Thus i p() = p + α( o) 4 o r 4 ( o ) + 2β( o) 2 o r 2 ( o ) o k(ξ)dξ = p i, (3.32) i i o p() = p + α( o) 4 o r 4 ( o ) + 2β( o) 2 o r 2 ( o ) o k(ξ)dξ. (3.33) It is seen that pressure is uniform for < i. Now the radial stress has the following distribution i σ rr () = α + 2β p i = σ i, (3.34) i o σ rr () = α()4 r 4 () + 2β()2 r 2 () p α( o) 4 o r 4 ( o ) 2β( o) 2 o r 2 ( o ) + o k(ξ)dξ. (3.35) It is seen that the radial stress is uniform and equal to σ i for < i. components of stress are also equal to σ i for < i. Note that the other two physical emark 3.2. Note that the defective ball does not have to be homogenous. One can have different energy functions for < i and > i. In this case W = W (, I, I 2 ) and hence α and β will have jumps at = i. This will not affect the pressure for < i. However, for > i, one should add the term [[α + 2β]] i to the pressure field Compressible isotropic solids Next we consider a spherically-symmetric point defect distribution in a ball made of a compressible isotropic solid. For an isotropic solid instead of considering the strain energy density as a function of the principal

8 3.2 Calculation of the residual stress field 8 invariants of C one can assume that W explicitly depends on the principal invariants of U, i.e. W = Ŵ (i, i 2, i 3 ), where i = λ + λ 2 + λ 3, i 2 = λ λ 2 + λ 2 λ 3 + λ 3 λ, i 3 = λ λ 2 λ 3. (3.36) Carroll [988] rewrote the representation of the Cauchy stress for isotropic elastic solids in terms of the left stretch tensor. In our geometric framework Carroll [988] s Eq. (2.5) is rewritten as [Yavari and Goriely, 23a] σ = ( i 2 Ŵ + Ŵ ) g + Ŵ V Ŵ V. i 3 i 2 i 3 i 3 i i 2 (3.37) In components this reads σ ab = ( i 2 i 3 Ŵ i 2 + Ŵ i 3 ) g ab + i 3 Ŵ i V ab Ŵ i 2 (V ) ab, (3.38) where b ab = V am V bn g mn and c ab = (V ) am (V ) bn g mn. Carroll [988] considered a special class of compressible materials for which Ŵ (i, i 2, i 3 ) = u(i ) + v(i 2 ) + w(i 3 ), where u, v, and w are arbitrary C 2 functions. For this class of materials we have σ = ( i 2 i 3 v (i 2 ) + w (i 3 )) g + u (i ) i 3 V v (i 2 )V. (3.39) In the case of a ball with a spherically-symmetric point defect distribution, we have λ = r () f(), λ 2 = λ 3 = r(). (3.4) Thus A simple calculation gives us i = r () f() + 2r(), i 2 = 2r()r () + r2 (), i f() 2 3 = r ()r 2 (). (3.4) 2 f() V = r () f() r() r() sin 2 Θ, V = f() r () r 3 () r 3 () sin 2 Θ. (3.42) Hence, the non-zero stress components read σ rr () = u (i ) 2 r 2 () + v (i 2 ) 2 r() + w (i 3 ), (3.43) σ θθ () = u f() (i ) r 3 ()r () + v (i 2 ) ( r 3 () + f() r 2 ()r () ) + w (i 3 ) r 2 (), (3.44) σ φφ () = sin 2 Θ σθθ (). (3.45) The equilibrium equation (3.23) is simplified to read 2 r 2 du dr + 2 r dv dr + dw u + 2( f) ( dr r 2 r + v ) =. (3.46) rr We first consider a harmonic material [John, 96] for which v(i 2 ) = c 2 (i 2 3) and w(i 3 ) = c 3 (i 3 ), where c 2 and c 3 are constants (Class I materials according to Carroll [988]). In this case the above ODE is reduced to du dr u + 2( f) ( r + c 2r ) =. (3.47) 2 r

9 3.3 An analogy between the point defect metric and the Schwarzschild metric 9 Let us now consider the point defect distribution (3.27). For < i, f() = and hence du dr that i must be a constant and therefore =. This implies r() = C + C 2 2. (3.48) For r() to be bounded at the origin we must have C 2 =. Now for < i, i = 3C and hence the physical components of Cauchy stress read σ rr () = σ rr () = u (3C ) C 2 + 2c 2 C + c 3 = σ, (3.49) σ θθ () = r 2 ()σ θθ () = u (3C ) C 2 + 2c 2 + c 3 = σ, C (3.5) σ φφ () = σ θθ (), (3.5) i.e. the Cauchy stress inside the point defect sphere is uniform and hydrostatic. For Class II and III materials according to Carroll [988], u(i ) = c (i 3), w(i 3 ) = c 3 (i 3 ) and u(i ) = c (i 3), v(i 2 ) = c 2 (i 2 3), respectively. For < i, for Class II materials, we have r 2 () = C 2 + C2. Similarly, for Class III materials we have r 3 () = C 3 + C 2. Assuming that r() = in both cases C 2 = and hence for < i, we have r() = α, where α is a constant. This is identical to what we observed for harmonic materials. Therefore, the above result holds for materials of Types II and III as well. The unknown constant C is determined after one solves a nonlinear second-order ODE for r in the interval i < o and imposes the continuity conditions r(i ) = r(+ i ), σrr (i ) = σrr (i ), and the boundary condition σrr ( o ) = p. 3.3 An analogy between the point defect metric and the Schwarzschild metric Einstein s vacuum field equations can be solved exactly for a spherically-symmetric distribution of matter with gravitational mass m. The solution is called the Schwarzschild (exterior) solution. In the coordinates (t,, Θ, Φ) for space-time, the Schwarzschild metric reads [Hawking and Ellis, 973] ds 2 = ( 2m ) dt2 + ( 2m ) d dθ sin 2 ΘdΦ 2. (3.52) This metric represents the gravitational field outside of a ball of mass m. Note that this solution is valid only for > 2m. The interior solution can be determined using the energy-momentum tensor of the matter inside the ball. Note that when restricted to (, Θ, Φ) this metric looks very similar to our point defect metric if one replaces ( 2m ) by f(). Next we consider a single point defect and observe another interesting similarity between our metric and that of Schwarzschild. 3.4 Singularity in the material metric in the case of a single point defect We consider a single point defect with strength δv at the center of the ball. In this case n() = where δ() is the one-dimensional Dirac delta distribution. Therefore δv δ(), (3.53) 4π2 Hence f() = δv δ() 4π 2 3δv 4π 3 h() = 3δv 4π 3. (3.54) = 3 δv 4π δ() 3 3δv 4π = 3δv 4π 3. (3.55)

10 Note that f() > and hence this expression is meaningful only when 3.5 Exterior residual stress field of a ball of point defects 3.5 Exterior residual stress field of a ball of point defects > ( 3δv 4π ) 3. (3.56) In Yavari and Goriely [22b] we considered a finite ball of radius o with a uniform defect distribution n in a small ball of radius i and showed that the stress inside the defective ball is uniform for < i. Let us now assume that n() = for > i but is elsewhere arbitrary. The total volume of the point defects is Note that for > i, we have δv = i f() = 4π 2 n()d. (3.57) 3δv 4π 3. (3.58) It is seen that for > i, f() depends only on δv and not on the specific distribution of n() for < i. Therefore, r() and consequently all the stress components for > i depend only on δv. Thus, we have proved the following proposition. Proposition 3.3. Consider a ball of radius o made of an isotropic elastic material or a material with anisotropy respecting the spherical symmetry. Assume that the ball is defect free outside a ball of radius i < o. Then, the residual stress field for > i depends only on the total volume of the point defects in < i and is independent of the specific form of n() for < i. emark 3.4. Note that this result is similar to the effect of a spherical ball of matter with mass m on the gravitational field. The gravitational field of the space-time outside the ball depends only on m and not on the specific distribution of density inside the ball (as long as it is spherically symmetric). 3.6 Isotropic distribution of point defects and pure dilatational eigenstrains We know that in the limit of a vanishing inclusion with pure dilatational eigenstrain the linearized solution for a single point defect in an isotropic linear elastic solid is recovered as shown in [Yavari and Goriely, 22b] by fixing δv = 4π 3 i n /3 and in the limit of small i. Note that our point defect metric is equivalent to that of a distributed radial eigenstrain. A natural question is: can we represent an isotropic distribution of point defects by a pure dilatational eigenstrain distribution? We will see in the following that the answer is negative. Consider a coframe field of the following form ϑ = K()d, ϑ 2 = K()dΘ, ϑ 3 = K() sin Θ dφ, (3.59) for some unknown function K() to be determined. Assuming that the non-metricity is traceless and isotropic Q αβ = 2δ αβ q()ϑ, the matrix of connection -forms has the following form ω = [ω α β] = where ω = ω 2 2 = ω 3 3 = q()ϑ, for a function q() to be calculated. Note that dϑ =, dϑ 2 = ω ω 2 ω 3 ω 2 ω 2 2 ω 2 3, (3.6) ω 3 ω 2 3 ω 3 3 K() [ + K () K() ] ϑ ϑ 2, dϑ 3 = K() [ + K () K() ] ϑ3 ϑ + cot Θ K() ϑ2 ϑ 3. (3.6)

11 4 Conclusions Cartan s first structural equations read This gives us T = ω 2 ϑ 2 ω 3 ϑ 3 =, (3.62) T 2 = { K() [ + K () K() ] + q()} ϑ ϑ 2 ω 2 ϑ + ω 2 3 ϑ 3 =, (3.63) T 3 = cot Θ K() ϑ2 ϑ 3 { K() [ + K () K() ] + q()} ϑ3 ϑ + ω 3 ϑ ω 2 3 ϑ 2 =. (3.64) ω 2 = { K() [ + K () K() ] + q()} ϑ2, ω 2 3 = cot Θ K() ϑ3, ω 3 = { K() [ + K () K() ] + q()} ϑ3. (3.65) It can be shown that = 2 2 = 3 3 = are trivially satisfied. The remaining Cartan s second structural equations read The first equation gives us the following ODE with solution 2 = 2 = dω 2 + ω 3 ω 2 3 =, (3.66) 2 3 = 3 2 = dω ω 2 ω 3 =, (3.67) 3 = 3 = dω 3 + ω 2 3 ω 2 =. (3.68) d d { K() [ + K () K() ] + q()} + ( + K () K() ) { K() [ + K () ] + q()} =, (3.69) K() K() [ + K () K() ] + q() = C K(). (3.7) Note that when q() =, we have K() = and hence C =. Therefore q() = K() K() [ + K () K() ] = K () K 2 () = ( K() ). (3.7) Interestingly, the other two curvature 2-forms trivially vanish. The Weyl -form is written as The governing differential equation for the function h() reads Q = 2q()ϑ = 2K () d = 2d ln K(). (3.72) K() d ln h() = 3 Q = 3d ln K(). (3.73) 2 Therefore, h() = CK 3 (). Because for K() =, h() =, C =, and hence h() = K 3 (). Note that µ = hϑ ϑ 2 ϑ 3 = h()k 3 ()µ = µ, and hence n() =. This means that the metric (3.59) cannot represent a spherically-symmetric distribution of point defects. 4 Conclusions We discussed the relevance of non-metricity in the nonlinear mechanics of distributed point defects. An anelasticity problem is transformed to a classical nonlinear elasticity problem if one can construct the material manifold, i.e. a 3-manifold in which the defective body is stress-free by construction. The material manifold of a solid with distributed point defects is a flat Weyl manifold, i.e., a manifold with a connection and metric such that the non-metricity is traceless and both the torsion and the curvature tensors vanish. We revisited the problem

12 EFEENCES 2 of a finite ball with a spherically-symmetric and isotropic distribution of point defects. We constructed the material manifold and calculated the residual stress field when the ball is made of an arbitrary incompressible isotropic solid. We observed an interesting analogy between the residual stress field and the gravitational field of space-time with a ball made of matter. We also showed that an isotropic distribution of point defects cannot be represented by a distribution of pure dilatational eigenstrains. Acknowledgements. AY was partially supported by AFOS Grant No. FA and NSF Grant No. CMMI and CMMI AG is Wolfson/oyal Society Merit Award Holder and also acknowledges support from a eintegration Grant under EC Framework VII. eferences Bilby, B. A.,., Bullough, and E., Smith [955], Continuous distributions of dislocations: a new application of the methods of non-iemannian geometry. Proceedings of the oyal Society of London A23(85): Bilby, B. A., and E., Smith [956], Continuous distributions of dislocations. III. Proceedings of the oyal Society of London A236(27): Carroll, M.M. [988], Finite strain solutions in compressible isotropic elasticity. Journal of Elasticity 2(): de Wit,. [98], A view of the relation between the continuum theory of lattice defects and non-euclidean geometry in the linear approximation. International Journal of Engineering Science 9(2): Doyle, T. C. and Ericksen, J. L. [956]. Nonlinear Elasticity. in Advances in Applied Mechanics IV, Academic Press, New York. Eshelby, J.D. [954], Distortion of a crystal by point imperfections. Journal of Applied Physics 25(2): John, F. [96], Plane strain problems for a perfectly elastic material of harmonic type. Communications on Pure and Applied Mathematics 3(2): Falk, F. [98], Theory of elasticity of coherent inclusions by means of non-metric geometry. Journal of Elasticity (4): Grachev, A.V., A. I. Nesterov, et al. [989], The gauge theory of point defects. Physica Status Solidi (b) 56(2):43-4. Hawking, S.W. and Ellis, G.F.. [973], The Large Scale Structure of Space-Time. Cambridge University Press, Cambridge. Kondo, K. [955], Geometry of elastic deformation and incompatibility, Memoirs of the Unifying Study of the Basic Problems in Engineering Science by Means of Geometry, (K. Kondo, ed.), vol., Division C, Gakujutsu Bunken Fukyo-Kai, 955, pp Kondo, K. [955], Non-iemannian geometry of imperfect crystals from a macroscopic viewpoint, Memoirs of the Unifying Study of the Basic Problems in Engineering Science by Means of Geometry, (K. Kondo, ed.), vol., Division D-I, Gakujutsu Bunken Fukyo-Kai, 955, pp. 6-7 Kröner, E. [99], The differential geometry of elementary point and line defects in Bravais crystals. International Journal of Theoretical Physics 29(): Love, A.H. [927], Mathematical Theory of Elasticity. Cambridge University Press, Cambridge. Miri, M and ivier, N. [22], Continuum elasticity with topological defects, including dislocations and extramatter. Journal of Physics A - Mathematical and General 35: Ogden,.W. [984], Non-Linear Elastic Deformations, Dover, New York.

13 EFEENCES 3 Saa, A. [995], Volume-forms and minimal action principles in affine manifolds. Journal of Geometry and Physics 5:2-8. Simo, J.C. and J.E. Marsden[983], Stress tensors, iemannian metrics and the alternative representations of elasticity, Springer LNP 95: Teodosiu, C. [982], Elastic Models of Crystal Defects. Springer-Verlag, Berlin. Volterra, V. [97], Sur l équilibre des corps élastiques multiplement connexes. Annales Scientifiques de l Ecole Normale Supérieure, Paris 24(3):4-58. Yavari, A. and A. Goriely [22a] iemann-cartan geometry of nonlinear dislocation mechanics. Archive for ational Mechanics and Analysis 25():59-8. Yavari, A. and A. Goriely [22b] Weyl geometry and the nonlinear mechanics of distributed point defects, Proceedings of the oyal Society A 468: Yavari, A. and A. Goriely [23a] Nonlinear elastic inclusions in isotropic solids, Proceedings of the oyal Society A 469: Yavari, A. and A. Goriely [23b] iemann-cartan geometry of nonlinear disclination mechanics. Mathematics and Mechanics of Solids 8(): 9-2.