Next Article in Journal
Movable Die and Loading Path Design in Tube Hydroforming of Irregular Bellows
Next Article in Special Issue
Effects of Manganese and Zirconium Dispersoids on Strain Localization in Aluminum Alloys
Previous Article in Journal
A Review of Microstructural Evolution and Modelling of Aluminium Alloys under Hot Forming Conditions
Previous Article in Special Issue
Quasi-Particle Approach to the Autowave Physics of Metal Plasticity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Revisiting the Application of Field Dislocation and Disclination Mechanics to Grain Boundaries

Université de Lorraine, CNRS, LEM3, 7 rue Félix Savart, 57070 Metz, France
*
Author to whom correspondence should be addressed.
Metals 2020, 10(11), 1517; https://doi.org/10.3390/met10111517
Submission received: 28 September 2020 / Revised: 29 October 2020 / Accepted: 9 November 2020 / Published: 16 November 2020
(This article belongs to the Special Issue Trends in Plasticity of Metals and Alloys)

Abstract

:
We review the mechanical theory of dislocation and disclination density fields and its application to grain boundary modeling. The theory accounts for the incompatibility of the elastic strain and curvature tensors due to the presence of dislocations and disclinations. The free energy density is assumed to be quadratic in elastic strain and curvature and has nonlocal character. The balance of loads in the body is described by higher-order equations using the work-conjugates of the strain and curvature tensors, i.e., the stress and couple-stress tensors. Conservation statements for the translational and rotational discontinuities provide a dynamic framework for dislocation and disclination motion in terms of transport relationships. Plasticity of the body is therefore viewed as being mediated by both dislocation and disclination motion. The driving forces for these motions are identified from the mechanical dissipation, which provides guidelines for the admissible constitutive relations. On this basis, the theory is expressed as a set of partial differential equations where the unknowns are the material displacement and the dislocation and disclination density fields. The theory is applied in cases where rotational defects matter in the structure and deformation of the body, such as grain boundaries in polycrystals and grain boundary-mediated plasticity. Characteristic examples are provided for the grain boundary structure in terms of periodic arrays of disclination dipoles and for grain boundary migration under applied shear.

1. Introduction

Disclinations and dislocations were simultaneously introduced by Volterra, more than a century ago [1]. Dislocations are crystal defects arising from translational lattice incompatibility, as measured by the Burgers vector. Disclinations are defects originating in the rotational incompatibility of the lattice characterized by the Frank vector [2]. Disclinations have long been neglected in the field theory of crystal defects, due to their diverging elastic energy, which precludes their occurrence as single defects [3]. However self-screened configurations, such as disclination dipoles, involve bounded elastic energy and may appear frequently [4,5]. Indeed, they have been proposed for describing grain boundaries in various polycrystals, by deriving the disclination density field from crystal orientation maps, as obtained by electron backscattered diffraction [6,7]. Disclination dipoles were also exhibited in planar configurations formed by C 60 molecular layers on W O 2 / W ( 110 ) substrates [8]. Generally, disclinations are liable for complementing dislocations in the description of the lattice structure when single-valued elastic rotation fields do not exist, as in polycrystals. Disclination dipoles were shown to provide a good description of grain boundaries [9,10,11,12].
In addition to obstacles to plastic deformation by dislocation glide and twinning, grain boundaries can also contribute to plastic deformation. In ultra-fine-grained metals or in materials lacking dislocation slip systems, they can become plasticity mediators, not only by supplying dislocation sources, but also via migration mechanisms. For example, when a shear stress is applied to a bicrystal, one crystal can grow at the expense of the other by the normal motion of the boundary, inducing plastic shear in the crystal traversed by the boundary. Such shear/migration coupling is characterized by the ratio β of the displacements parallel and normal to the boundary. Surface-dislocation-based models [13,14] captured well the variations of this coupling factor as a function of the misorientation from the Burgers vector’s content of the interface. Molecular dynamics methods [14,15,16] provided elementary mechanisms for boundary migration under stress. As considered here, field models involving dislocation and disclination densities are also candidates for the modeling of grain boundaries and their plasticity. The size of the investigated microstructures can be of the order of (but can be larger than) the dimensions typically considered in atomistic simulations. The field models allow analyzing the fine structure of the grain boundaries and the distribution of elastic energy at an atomic resolution length scale. It is also possible to model boundary migration by the motion of arrays of disclination dipoles, and to model dissipative and diffusive phenomena using disclination and dislocation transport equations [17].
In the present paper, we provide a review of the mechanical theory of dislocation and disclination fields [18] and show its interest in dealing with crystalline bodies subject to grain boundary-mediated plasticity. We limit the presentation to infinitesimal transformations. A more complete presentation of the model may be found in the framework of generalized disclinations [19]. By a field theory, we mean a theory of continuously distributed crystal defect densities using the tools of the mathematical theory of partial differential equations and boundary problem solving. As such, we capitalize on the earlier elasto-static theories of dislocation fields [20,21] and of dislocation and disclination fields [2]. Both theories are strongly connected: the latter reduces to the former when the disclination density field vanishes. We also benefit from the mechanical theory of dislocation fields [22], which additionally deals with dislocation-mediated plasticity. In the latter, the material displacement and dislocation density fields are derived from a set of partial differential equations complemented with initial and boundary conditions, provided constitutive information for the elastic behavior and dislocation motion is supplied [23]. This theory was shown to capture well complex features of elasto-plastic deformation such as size effects [24,25], directional hardening [24,26], plastic strain localization [27], dislocation cores [28,29] and patterning of dislocation ensembles [30,31].
When the dislocation velocity vanishes, the mechanical theory of dislocation fields [22] reduces to the elasto-static theory of dislocations [20,21]. The present framework of disclination and dislocation fields reduces to the mechanics of dislocation fields when the disclination density is ignored. It also reduces to the elasto-static crystal theory of dislocations and disclinations [2] when the velocities of dislocations and disclinations are ignored. A previous theory of plasticity accounting for both the translational and rotational aspects of lattice incompatibility was proposed by Kossecka and deWit [32,33]. Driving forces for disclination motion were discussed in [34]. The theory offered kinematic guidelines, but it was not used to solve realistic boundary value problems. In the present work, we try to use the field disclination and dislocation mechanics model in situations where the dynamic interplay between dislocations and disclinations is essential to the understanding of the plastic deformation of the material.
The outline of the paper is the following. Mathematical notations are provided in Section 2. A review of the incompatible elastic dislocation and disclination model [2] is given in Section 3. Static coupling between dislocations and disclinations occurs at this level, due to the continuity condition for dislocation densities. The equations of balance of momentum and moment of momentum are given in Section 4, which further includes the presentation of nonlinear nonlocal elastic constitutive laws that can be of interest for core problems. Section 5 is devoted to the transport properties of dislocations and disclinations. Dynamic coupling occurs at this stage, because disclination mobility yields a source/sink term for dislocations. In Section 6, guidance of the Clausius–Duhem inequality allows defining the driving forces for dislocation and disclination motion, and helps providing appropriate constitutive relations between the driving forces and defects velocities. Section 7 shows possible algorithms used for the solution of the model equations. A plane “edge-wedge” model is detailed in Section 8, and used in Section 9 and Section 10 to investigate the structure and mobility of grain boundaries seen as periodic arrays of disclination dipoles.

2. Mathematical Notations

A bold symbol is used for a tensor, as in: A . To avoid possible ambiguities, an upper arrow can be used for a vector: V . The transpose of a tensor A is A t . Tensor subscript indices are expressed with respect to a basis ( e i , i = 1 , 2 , 3 ) of a rectangular Cartesian coordinate system. Arrays of one/two dots denote contraction of the respective number of “adjacent” indices on two immediately neighboring tensors. For instance, the tensor A . B with components A i k B k j comes from the dot product of tensors A and B , and A : B = A i j B i j is the inner product. The cross product of a second-order tensor A and a vector V , the div and curl operators acting on second-order tensors are defined by:
( A × V ) i j = e j k l A i k V l
( div   A ) i = A i j , j
( curl   A ) i j = e j k l A i l , k .
where e j k l = e j · ( e k × e l ) is a component of the Levi–Civita tensor X . It is 1 if the j k l permutation is even, 1 if it is odd, 0 otherwise. The comma followed by a component index denotes a spatial derivative with respect to the Cartesian coordinate. A vector A is associated with tensor A through the inner product of A with the tensor X :
( A ) k = 1 2 ( X : A ) k = 1 2 e k i j A i j
( A ) i j = ( X · A ) i j = e i j k ( A ) k .
The trace of the second-order tensor A is t r ( A ) = A i i . The symmetric/skew-symmetric parts are A s y m and A s k e w . The hydrostatic part is A h y d = 1 3 t r ( A ) I , where I is the identity tensor, and its deviatoric part: A d e v = A A h y d . An upper dot represents a material time derivative.

3. Review of the Incompatible Elasto-Static Defect Theory

3.1. Incompatibility in the Dislocation Model

In the proposed framework, continuity of the displacement vector field u , along with continuity of its derivatives, is assumed at any point, including between atoms, in a simply connected body that is elasto-plastically deformed, except perhaps along interfaces where continuity of the derivatives may not hold. The rotation vector
ω = 1 2 curl   u
is also defined at any point. These statements assume in a dual way that the material is capable of transmitting stresses and couple stresses, possibly below inter-atomic length scales. Such continuity was very early proposed [35,36] when defining pressure and stresses in quantum systems. It was later used by [37] when defining the stress tensor in quantum mechanics. By definition, couple stresses (moments of stress) at a point come from mechanical stresses existing at distant material points. As such, couple stresses introduce nonlocality in the mechanical description of the material. The distortion tensor is the gradient of the displacement U = grad   u . It is curl-free:
curl   U = 0 .
If dislocations are present, the elastic U e and plastic U p components of U contain incompatible parts, which are not curl-free, but are div-free. Indeed, if dislocations thread a patch S in the body, a constant discontinuity b in the elastic displacement exists across S, which manifests itself as a closure defect along the circuit C surrounding S:
b = u e = C U e · dl .
b is referred to as the resulting Burgers vector of this dislocation ensemble. Applying Stoke’s theorem yields:
b = C U e · dl = S curl   U e · n d S ,
where n is the unit normal to S. The discontinuity b is also characterized in a pointwise continuous manner by introducing Nye’s dislocation density tensor α :
curl   U e = α
such that
b = S α · n d S .
Thus, dislocation lines threading S and having a non-zero resulting Burgers vector induce the existence of an incompatible part U e of the elastic distortion tensor. Applying the Stokes-Helmholtz decomposition of square-integrable tensor fields with square-integrable first order derivatives [38], we can uniquely define the square-integrable tensor and vector χ and z such that the elastic distortion U e reads:
U e = curl   χ + grad   z .
The curl of U e extracts curl   χ and removes grad   z while, conversely, taking the divergence extracts grad   z and removes curl   χ . Thus, Equation (10) involves only curl   χ , and
curl   U e = curl   curl   χ = α .
The term grad   z is the compatible (curl-free) part U e of the elastic distortion U e , and z is the compatible elastic displacement u e , up to a constant. However, Equation (13) is insufficient to identify U e from a given dislocation field α . To ensure completeness of the identification, U e must vanish if α = 0 . Following [38,39], we thus augment Equation (13) with the supplementary condition div   U e = 0 together with U e · n = 0 on the external boundary with unit normal n . This allows ensuring that the solution does not contain any gradient part. Applying the curl operator to Equation (13), we have:
curl   curl   U e = grad   div   U e div   grad   U e = curl   α .
Then, using the side condition div   U e = 0 , we get:
div   grad   U e = curl   α .
Equation (15) is a Poisson-type equation for the unique determination of the unknown tensor U e . We could have demonstrated in a similar way the existence of an incompatible part of the plastic distortion, U p , which is opposite to the incompatible elastic distortion U e to ensure the compatibility of the distortion U and the continuity of the displacement u . An additional compatible component U p may also exist, together with a compatible plastic displacement u p . We can then write the following equations:
U = U e + U p
U e = U e + U e
U p = U p + U p
0 = U e + U p .
From the above, we see that:
α = curl   U p
similarly to Equation (13). The continuity condition:
div   α = 0
derives naturally from Equations (10) and (13). The divergence-free character of α reflects the fact that dislocation lines cannot end within the material. They should either exit the material, or be closed lines. As will be shown later in Section 3.2, dislocation lines may possibly end on disclination lines. Now, in view of defining disclinations in the next section, we introduce curvatures and examine how they appear in the present dislocation model. The strain tensor ϵ is the symmetric part of U , the rotation tensor ω its skew-symmetric part, corresponding to the rotation vector ω through the relation:
ω = X · ω ,   ω = 1 2 ω : X .
With this, Equation (7) yields:
curl   ϵ + d i v ( ω ) I grad t ω = 0
where grad ω is the curvature tensor κ :
κ = grad   ω .
The trace of Equation (23) yields:
d i v ( ω ) = 0 ,
and implies that t r ( κ ) = 0 , meaning that the curvature tensor is deviatoric. Please note that transposing Equation (23) and then taking the curl yields:
curl   curl t   ϵ = 0 .
This is the Saint-Venant compatibility condition. Applying the above procedure to Equation (10), we find:
curl   ϵ e + d i v ( ω e ) I grad t ω e = α
and:
η = curl   curl t ϵ e = curl   ( α t 1 2 t r ( α ) I ) .
The incompatibility tensor η is a possible measure of the incompatibility due to the presence of dislocations. To terminate this section, we introduce the elastic curvature tensor κ e as the gradient of the elastic rotation ω e :
κ e = grad   ω e ,
which reflects both the continuity of the rotation vector and the curl-free character of the tensor κ e :
curl   κ e = 0 .
Using κ e , Equation (27) can be expressed as:
curl   ϵ e + t r ( κ e ) I κ e t = α .
It is important to note for what follows that the above equation holds even if κ e is not a gradient. Finally, relations similar to Equations (29) and (31) can be obtained by using Equation (19) for the plastic strain ϵ p , rotation ω p and curvature κ p :
κ p = grad   ω p
curl   ϵ p t r ( κ p ) I + κ p t = α .

3.2. Incompatibility in the Dislocation and Disclination Model

We now assume that elastic rotation discontinuities ω e exist in addition to the displacement discontinuities u e = b reflected by the dislocation densities. Thus, the elastic curvature tensor κ e is not a gradient anymore. This is the situation met along grain boundaries in polycrystals. Integrating κ e along closed circuits in the manner of Equation (9) leads to the angular closure defect ϕ known as the Frank vector:
ϕ = ω e = C κ e · dr = S curl   κ e · n d S .
We introduce a tensor field θ such that:
curl   κ e = θ
and
ϕ = S θ · n d S .
θ is the disclination density tensor [2]. In the presence of disclinations, the elastic displacement discontinuity reads:
u e = C [ ϵ e ( κ e t × r ) t ] · dr + ϕ × r 0 = b + ϕ × r 0 .
Thus, the discontinuity u e derives from the motion of a rigid body with translation b and rotation ϕ , in agreement with Weingarten’s theorem [40]. The Burgers vector can also be expressed as a function of dislocation and disclination density tensors:
b = S ( α ( θ t × r ) t ) · n d S .
When both dislocations and disclinations exist, the incompatibility tensor η is:
η = curl   curl t   ϵ e = curl   ( α t 1 2 t r ( α ) I ) + θ .
As similarly detailed above for U e , we must ensure that the incompatible part κ e of the curvature tensor is null if θ = 0 . In this aim, Equation (35) must be replaced with:
curl   κ e = θ ,
together with the side condition div   κ e = 0 and boundary condition κ e · n = 0 on external boundaries.
These conditions provide a unique solution κ e of the Poisson equation:
div   grad   κ e = curl   θ .
A compatible part κ e can also exist, but it is not involved in Equations (40) and (41). Therefore, one has the following equations:
κ e = κ e + κ e
κ p = κ p + κ p
κ = κ e + κ p
0 = κ e + κ p .
The incompatible plastic curvature is opposed to the incompatible elastic curvature to ensure that the total curvature is compatible and the material continuous. Then:
curl   κ p = θ
div   grad   κ p = curl   θ .
The continuity equation for disclinations:
div   θ = 0
derives from Equations (35) and (40). It mathematically expresses that disclination lines should not end within the body and form closed loops. The continuity equation for dislocations Equation (21) becomes:
div   α = θ : X .
It means that dislocations may now terminate on twist-disclination lines. Please note that if the disclination density is null, Equations (31), (35), (38) and (49) reduce to Equations (11), (21), (27) and (30), and the elasto-static model of dislocations and disclinations reduces to the above dislocation model.

4. Elasticity in Incompatible Media Containing Dislocations and Disclinations

4.1. Field Equations

We consider describing the mechanical equilibrium of a simply connected body of volume V and external surface S containing dislocations and disclinations, and submitted to external loading. Since the disclination densities are associated with the incompatibility of the elastic curvatures, we are interested in a framework involving, in addition to stresses, the work-conjugates of the curvatures, i.e., the couple stresses. We therefore follow the footsteps of [41] to derive the mechanical balance equations. Let n be the outward unit normal to S, t be the traction vector (force per unit area) and m the moment vector (couple per unit area) acting on S. As they are non-essential for the present purposes, body forces and couples are neglected for the sake of simplicity, as well as inertia forces. The conservation of momentum and moment of momentum of the body therefore imposes the following conditions:
S t d S = 0
S ( m + r × t ) = 0 .
Assuming that the tractions and moments at boundaries can be written as t = T · n and m = M · n , where T and M are respectively the stress and couple-stress tensors, we find by substituting in Equations (50) and (51) and applying the divergence theorem:
V div   T d v = 0
V ( div   M X : T ) d v = 0 .
Assuming sufficient continuity, the pointwise local equations are:
div   T = 0
div   M X : T = 0 .
Equations (54) and (55) are balance of momentum and moment of momentum equations where the stress tensor T is not necessarily symmetric. These two equations can be combined to form a higher-order balance equation coupling the symmetrical part of the stress tensor T s y m and the deviatoric part of the couple-stress tensor M d e v [41]:
div   T s y m + 1 2 curl   div   M d e v = 0 ,
but involving neither the skew-symmetric part T s k e w of the stress tensor nor the hydrostatic part M h y d of the couple-stress tensor. T s k e w may be obtained from Equation (55), but M h y d is altogether absent from the mechanical balance equations and is therefore indeterminate in such a couple-stress theory [41,42,43]. A more general presentation of the theory involving not only curvatures and couple stresses but also strain gradients and their work-conjugates (“hyper-stresses” or “double-stresses”) allows avoiding such indeterminacy. It features crystal defects more general than disclinations, and referred to as generalized disclinations. The differences between a couple-stress disclination theory involving only curvatures and a hyper-stress generalized disclination theory involving both curvatures and strain gradients are discussed in [44].
The mechanical dissipation in the material is the difference between the power of the applied forces and the rate of change of the stored energy:
D = S ( v · T + ω ˙ · M ) · n d S V ψ ˙ d v = S ( v i T i j + ω i ˙ M i j ) n j d S V ψ ˙ d v .
Here, v is the material velocity and ψ the specific free energy density. Using the balance Equations (54) and (55) together with the divergence theorem, we can write:
D = D [ ( v i T i j + ω ˙ i M i j ) , j v i T i j , j ω i ˙ ( M i j , j e i k l T k l ) ψ ˙ ] d v .
Then we use ω i = e i m n ω m n / 2 to substitute the rotation vector for the rotation tensor through Equation (22), and progressively derive:
D = D ( T i j v i , j + M i j ω ˙ i , j 1 2 e i k l e i m n ω ˙ m n T k l ψ ˙ ) d v
D = D ( T i j v i , j + M i j ω ˙ i , j 1 2 ( δ k m δ l n δ k n δ l m ) ω ˙ m n T k l ψ ˙ ) d v
D = D ( T i j v i , j + M i j ω ˙ i , j T i j ω ˙ i j ψ ˙ ) d v
D = D ( T i j ϵ ˙ i j + M i j ω ˙ i , j ψ ˙ ) d v
D = D ( T s y m : ϵ ˙ + M d e v : κ ˙ ψ ˙ ) d v .
The Clausius–Duhem inequality imposes that the dissipation be positive, or zero for reversible processes: D 0 . In Equation (63) the inner product between symmetric and skew-symmetric tensors is null. Hence the strain rate tensor ϵ ˙ extracts the symmetric part of the stress tensor T , and the deviatoric curvature rate tensor κ ˙ extracts the deviatoric part of the couple-stress tensor. A free energy density function ψ is taken as:
ψ = ψ ( ϵ e , κ e ) ,
with contributions from the elastic strain ϵ e and curvature κ e . A more detailed description of this free energy will be given below in Section 4.3. For the time being, it is sufficient to note that, differentiating Equation (64) on the one hand, and considering elastic reversible processes such that D = 0 in Equation (63) on the other hand, leads to identifying the symmetric part of the stress tensor T s y m and deviatoric part of the couple-stress tensor M d e v with the partial derivatives of the free energy density with respect to the elastic strain and curvature:
ψ ˙ = ψ ϵ e : ϵ ˙ e + ψ κ e : κ ˙ e
ψ ˙ = T s y m : ϵ ˙ e + M d e v : κ ˙ e
T s y m = ψ ϵ e
M d e v = ψ κ e .
Equations (67) and (68) will be the constitutive relationships describing the elastic behavior of the material.

4.2. Solving the Incompatible Elasticity Problem

Consider again a simply connected body of volume V containing dislocations and disclinations, submitted to tractions and moments on its external surface S with outward unit normal n . Tractions are imposed on a part S t of S and rotations/displacements on a part S u :
T · n = t d   on   S t
u e = u d   on   S u
ω e = ω d   on   S u .
We are interested in describing the state of stress and couple-stress in the body in the presence of dislocation and disclination density fields ( α , θ ) . The incompatible elastic strain and curvature ( ϵ e , κ e ) are assumed to be known from the solution of the Poisson Equations (15) and (41). They usually do not allow satisfying the balance Equation (56). Recalling the decomposition of the elastic strain and curvature in the sum of incompatible, ( ϵ e , κ e ) , and compatible, ( ϵ e , κ e ) , components, Equation (56) can be written as:
0 = div   T s y m ( ϵ e , κ e ) + 1 2 curl   div M d e v ( ϵ e , κ e ) + f ,
f = div   T s y m ( ϵ e , κ e ) + 1 2 curl   div M d e v ( ϵ e , κ e ) .
f is to be seen as a volumetric density of body forces associated with dislocations and disclinations. In Equations (72) and (73), the symmetric part of the stress tensor and the deviatoric part of the couple-stress tensor are obtained from the constitutive relations Equations (67) and (68). The solution of Equation (72) is the compatible elastic displacement field.

4.3. Elasticity in Dislocation/Disclination Cores

Elastic constitutive laws usually relate the elastic strain tensor to the stress tensor, both taken at the same material point. As such, they are referred to as “local” constitutive laws. They are valid in a statistical sense if the characteristic size of the volume element actually envisioned at the material point is large enough to include a large number of atoms in a perfectly ordered crystal. Local laws may fail at describing correctly elasticity if this characteristic size reduces to the order of a few atomic lattice spacings, particularly if the volume element involves crystal defects such as dislocations and disclinations. Since the mechanical behavior of a particular atom depends on its interactions with all neighboring atoms, whose location is strongly disturbed by the presence of these defects, defining the stress tensor in such a local manner may then become impossible. Nonlocal elastic laws are then needed. Nonlocal elastic laws are relationships using convolution integrals, where the stresses at a given material point depend on the deformation state in some neighborhood, whose extent is controlled by a nonlinear kernel function [45,46,47,48,49]. Such nonlocal elastic laws have been successful in retrieving phonon dispersion curves and regularizing elastic fields in the core region of crystal defects. Convolution integrals can be coupled with higher-order elastic laws to derive nonlocal elastic constitutive laws that are of interest in the core regions of dislocations and disclinations, while local laws may be sufficient in the non-defected regions. The occurrence of different elastic behavior in non-defected and defected regions was already assumed in grain and phase boundaries [50,51,52,53]. Following [54] and relying on Eringen-type convolution integrals applied to strains and curvatures, we propose below nonlocal laws describing the elastic behavior in the core region of dislocations and disclinations, and reducing to local laws in the perfectly ordered parts of the crystal.

Nonlocal Convolution Integrals

We introduce a free energy density functional at a point r in a body V in the form:
ψ ( r ) = 1 2 C i j k l ( r ) ϵ i j e ( r ) ϵ k l e ( r ) + 1 2 A i j k l ( r ) κ i j e ( r ) κ k l e ( r ) + . . . . . . + 1 2 ϵ i j e ( r ) V K ( R ) D i j k l ( r , r ) ϵ k l e ( r ) d V + . . . . . . + 1 2 κ i j e ( r ) V K ( R ) B i j k l ( r , r ) κ k l e ( r ) d V .
The first two pointwise terms on the r.h.s. include symmetric positive-definite quadratic forms with A and C given at point r . This leads to the linear elastic laws used and discussed in [2,18,54] in the framework of models for dislocations and disclinations. Nonlocal contributions are added in the form of Eringen-type convolutions, which are the last two terms of the r.h.s. [55]. In such terms, the location vectors r are points distinct from point r . The difference vectors R = ( r r ) and their norms R : = ( R . R ) 1 / 2 are also defined. The spatial range of nonlocality is controlled by the nonlocal kernel K ( R ) , while the elasticity tensors B and D depend on ( r , r ) and the elastic curvatures/strains are taken at points r . The proposed nonlocality is motivated by the strong coupling existing in the core region of defects between displacement/strain and rotation/curvature at distant locations r and r (see Figure 1). One can express the strain induced at point r by the curvature at point r :
ϵ k l e ( r ) = e k m n κ m l e ( r ) R n ,   ϵ e ( r ) = ( κ e t ( r ) × R ) t .
Note the presence of the above nonlocal term in the definition of the Burgers vector in Equation (37) above. Similarly, we can express the curvature at point r due to the strain at point r :
κ k l e ( r ) = e k m n ϵ m l e ( r ) R n R 2 ,   κ e ( r ) = 1 R 2 ( ϵ e ( r ) × R ) t .
The above nonlocal relations between strains and curvatures can be used to substitute curvatures and strains at points r for strains and curvatures at points r , respectively. This allows greatly simplifying the convolutions in the free energy functional Equation (74), in such a way that nonlocality applies only to elasticity tensors while all strains and curvatures are defined at point r (the reader is referred to ref [55] for more details):
ψ = 1 2 C i j k l ϵ i j e ϵ k l e + 1 2 A i j k l κ i j e κ k l e + 1 2 D ¯ i j k l ϵ i j e κ k l e + 1 2 B ¯ i j k l κ i j e ϵ k l e ,
D ¯ ( r ) = V K ( R ) ( D t × R ) t d V ,
B ¯ ( r ) = V K ( R ) ( B × R R 2 ) t d V .
By differentiating Equation (77) with respect to the elastic strain and curvature tensors, we find the stress and couple-stress tensors:
T i j s y m = ψ ϵ i j e = C i j k l ϵ k l e + 1 2 ( D ¯ i j k l + B ¯ k l i j ) κ k l e
M i j d e v = ψ κ i j e = A i j k l κ k l e + 1 2 ( B ¯ i j k l + D ¯ k l i j ) ϵ k l e .
Using finally E ¯ i j k l = 1 2 ( D ¯ i j k l + B ¯ k l i j ) , we can write:
T i j s y m = ψ ϵ i j e = C i j k l ϵ k l e + E ¯ i j k l κ k l e
M i j d e v = ψ κ i j e = A i j k l κ k l e + E ¯ k l i j ϵ k l e .
We end up with relationships similar to those proposed in [54], except that the present elastic tensors ( B ¯ , D ¯ , E ¯ ) involve nonlocal integrals. As proposed by Eringen [47,48], we assume D = C . We also impose B ¯ i j k l = D ¯ k l i j , which allows ensuring that the free energy functional is symmetric positive-definite and that the elastic solution is stable. Interestingly, stability of the elastic strain and curvature fields is expected at the convoluted level from the nonlocal character of the formulation, whereas instability would occur in the limit of a purely local formulation [55].

5. Transport of Dislocations and Disclinations

In this Section, we are interested in the kinematics of dislocations and disclinations in the context of the continuous setting set forth in Section 3 for their representation by smooth density fields. In this aim, we take advantage of the transport framework for the spatio-temporal evolution of the defects density fields that derives from the conservation of the Burgers and Frank vectors during their motion. Being unquestionable from a kinematic point of view, this framework, referred to as transport of dislocation and disclination densities, provides a natural basis for the dynamic description of plasticity through dislocation and disclination motion. Since disclinations are mainly distributed along grain boundaries in polycrystals, it provides in particular a basis for describing the contribution to plasticity of grain boundary motion.
Consider a material surface S delimited by a closed circuit C. Let f be the disclination flux used to measure the rate of inflow into S of disclination lines, carrying along with them their corresponding Frank vectors ϕ , through a line element dx of curve C. Let V θ be the velocity of the disclinations with respect to the material, and l the unit vector along the disclination lines. Then the disclination tensor is: θ = ϕ l . In the absence of a disclination source term, the conservation of the Frank’s vector content demands that the rate of change of the Frank’s vector of all disclination lines threading S be equal to the total disclination flux across C:
d d t S θ · n d S = C f · dx .
For small transformations, the pointwise statement corresponding to Equation (84) is:
θ ˙ = curl   f
where θ ˙ represents the time derivative of the disclination density tensor. The rate of inflow of Frank vectors across the surface l × dx is:
f . dx = ϕ ( V θ . ( l × dx ) ) .
Indeed, suppose first that the disclination velocity V θ is in the plane ( l , dx ) . Then, the disclination does not enter (or exit) surface S, and it is clear from Equation (86) that the flux f is zero as it should be. Conversely, f is always non-zero when V θ has a non-zero component normal to the plane ( l , dx ) , i.e., when the disclination enters or exits surface S, and it changes sign according to whether the disclination enters or exits S. We can then obtain from Equation (86):
f = ϕ l × V θ = θ × V θ .
Consequently, the local statement of balance Equation (85) becomes:
θ ˙ + curl   ( θ × V θ ) = 0 .
Equation (88) is a transport law for the disclination density tensor θ . It can be understood as an evolution equation for θ , provided the disclination velocity V θ is known as a function of the stress/couple-stress state from constitutive statements. Its meaning is that through the curl term, the incompatible part of the disclination flux incrementally feeds the disclination density. Comparing Equations (46) and (88) it follows, after time differentiation of Equation (46) that the cross product θ × V θ can be identified with the plastic curvature rate tensor κ ˙ p , up to a gradient:
κ ˙ p = θ × V θ + grad   ω ˙ p * .
At mesoscale, the term κ ˙ p * = grad   ω ˙ p * can be given the significance of a statistical plastic curvature rate, meaning that a non-vanishing plastic curvature rate may be obtained even when the disclination density θ vanishes at this scale [56]. Indeed, using space-time running averages of the disclination density tensor, disclination velocity V θ and plastic curvature rate tensor κ ˙ p over a domain of mesoscopic size, allows writing the mesoscopic plastic curvature rate κ ˙ p ¯ as:
κ ˙ p ¯ = θ × V θ ¯   =   θ ¯ × V θ ¯ + κ ˙ p * ,
where overbars indicate averaged variables. It is seen that κ ˙ p ¯ may be non-zero when the net disclination density vanishes at mesoscale ( θ ¯ = 0 ) , in which case it becomes:
κ ˙ p ¯ = κ ˙ p *   .
Dropping the overbars for convenience, the plastic curvature rate may be generally written as:
κ ˙ p = θ × V θ + κ ˙ p * .
As we shall discuss below, the compatible part κ ˙ p of κ ˙ p increments the history-dependent compatible plastic curvature produced by the motion of dislocations and disclinations. Equation (92) accounts for plasticity processes in relation to the mobility of grain boundaries that cannot be described by a pure dislocation-based model.
In the theory of dislocations, the dislocation density tensor α satisfies a transport equation similar to Equation (88), in the absence of a source term [22,57]:
α ˙ + curl   ( α × V α ) = 0 .
Here, V α denotes the dislocation velocity, and the cross product α × V α can be identified with the plastic distortion rate U ˙ p , up to a gradient. However, U ˙ p is undefined in a theory of disclinations, because its skew-symmetric part ω ˙ p has discontinuities. Hence, the dislocation transport Equation (93) does not account for curvature incompatibility and has to be modified in a disclination theory. For small transformations, a version correctly accounting for curvature incompatibility is readily obtained from the time derivative of Equation (33):
α ˙ = curl   ϵ p ˙ + κ ˙ p t t r ( κ ˙ p ) I = curl   ϵ p ˙ + s θ .
Here the plastic strain rate ϵ ˙ p is
ϵ ˙ p = ( α × V α ) s y m + ϵ ˙ p * ,
where ϵ ˙ p * is a gradient tensor acquiring significance at mesoscale, as discussed above for the mesoscopic plastic curvature rate tensor κ ˙ p * . The compatible part of the plastic strain rate, ϵ ˙ p , is not involved in Equation (94) but it nevertheless increments the compatible plastic strain produced by the motion of dislocations through the lattice. Note importantly in Equation (94) the presence of a dislocation source/sink term, s θ = κ ˙ p t t r ( κ ˙ p ) I , due to the mobility of disclinations. Such a term offers supplementary relaxation mechanisms [4,58], and will be shown to be useful later on in the prediction of shear-coupled boundary migration. When the dislocation/disclination velocities and mesoscopic contributions to plasticity ( κ ˙ p * , ϵ ˙ p * ) are known from the stress/couple-stress state through constitutive relations, Equations (88), (92), (94) and (95) can be used as evolution equations for the disclination and dislocation densities. Please note that Equations (88) and (94) are first order hyperbolic equations, which implies that their solutions have propagative character, as experimentally verified in [31]. Furthermore, hyperbolicity of these equations as a strong impact on the algorithms devoted to their numerical solution from finite elements or spectral approximations [59,60,61], as detailed below in Section 7.3.

6. Constitutive Relations for Dislocation and Disclination Mobility

In this Section, we account for thermodynamic requirements on positiveness of the mechanical dissipation to formulate constitutive relations for the mobility of dislocations and disclinations by following a procedure proposed in [62]. Substituting Equations (65) and (66) in Equation (63), and using the additivity of elastic and plastic strain rates and curvature rates tensors, we find the dissipation to be:
D = D ( T s y m : ϵ ˙ p + M d e v : κ ˙ p ) d v 0 .
To identify driving forces associated with dislocation and disclination velocities, we write down the plastic strain and curvature rate tensors ( ϵ ˙ p , κ ˙ p ) in terms of dislocation and disclination fluxes as in Equations (92) and (95) and substitute in Equation (96):
D = D ( T s y m : ( α × V α ) s y m + T s y m : ϵ ˙ p * + M d e v : ( θ × V θ ) + M d e v : κ ˙ p * ) d v ,
which reads in component form:
D = D ( 1 2 T i j s y m ( e j k l α i k V l α + e i k l α j k V l α ) + T i j s y m ϵ ˙ i j p , * + M i j d e v e j k l θ i k V l θ + M i j d e v κ ˙ i j p , * ) d v ,
or also:
D = D ( e j k l T i j s y m α i k V l α + T i j s y m ϵ ˙ i j p , * + e j k l M i j d e v θ i k V l θ + M i j d e v κ ˙ i j p , * ) d v .
Equivalently, Equation (99) is in intrinsic form:
D = D ( F α · V α + T s y m : ϵ ˙ p * + F θ · V θ + M d e v : κ ˙ p * ) d v ,
where we wrote:
F α = T s y m · α : X ;   F l α = e j k l T i j s y m α i k
F θ = M d e v , t · θ : X ;   F l θ = e j k l M i j d e v θ i k .
F α and F θ are the driving forces for dislocation and disclination motion. Using the dyadic notation α = b t for the dislocation density ( b is the Burgers vector, t the line vector), F α can be expressed as: F α = T s y m · b × t , a form similar to the Peach-Köhler force for discrete dislocations. Similarly for the disclination density tensor, θ = ϕ l ( ϕ is Frank’s vector, l the line vector), and F θ = M d e v , t · ϕ × l . This is a Peach-Köhler-type relationship for disclinations. Thus, F θ drives grain boundary motion in the present disclination theory. It is produced by the couple stresses, the stresses being not directly involved, and is normal to the disclination line. Positiveness of the mechanical dissipation rate Equation (99) can be satisfied by choosing the simple constitutive relations:
F α = B α V α ;   B α 0
F θ = B θ V θ ;   B θ 0
ϵ ˙ p * = λ ϵ T s y m ;   λ ϵ 0
κ ˙ p * = λ κ M d e v ;   λ κ 0 ,
where ( B α , B θ , λ ϵ , λ κ ) are positive material parameters. Equations (103)–(106) have the form of Newtonian viscous drag relationships, but ( B α , B θ , λ ϵ , λ κ ) are generally not simple constants. They may involve non-linearity, or/and hardening laws reflecting for example thermally assisted obstacle overcoming through adequate constitutive statements. Please note that both the glide and climb components of the Peach-Köhler force are included in Equation (103). λ ϵ describes plasticity through both dislocation glide and climb at mesoscale. Similarly, λ κ describes plasticity induced by grain boundary motion at mesoscale and driven by M d e v in Equation (106). Possible expressions for λ κ were proposed in [56] by averaging Equation (104).

7. Solution Algorithms

In the present Section, we propose two algorithms for the numerical solution of the model. The first version solves for the incompatible plastic strains and curvatures at each time step. The second version is a reduced, rate form of the latter, numerically faster, but where incompatible fields are not determined at each step.

7.1. Complete Algorithm

The complete model involves the following equations, gathered from the above sections:
curl   κ p = θ
curl   ϵ p = α + κ p t t r ( κ p ) I
T = C : ( ϵ ϵ p ) + D : ( κ κ p )
M d e v = A : ( κ κ p ) + B : ( ϵ ϵ p )
0 = div   T s y m + 1 2 curl   div   M d e v
κ ˙ p = θ × V θ + κ ˙ p *
ϵ ˙ p = ( α × V α ) s y m + ϵ ˙ p *
θ ˙ = curl   κ ˙ p
α ˙ = curl   ϵ p ˙ + κ ˙ p t t r ( κ ˙ p ) I
V α = 1 B α T s y m · α : X
V θ = 1 B θ M d e v , t · θ : X
ϵ ˙ p * = λ ϵ T s y m
κ ˙ p * = λ κ M d e v ,
with the boundary conditions:
T · n = t d   on   S t
u = u d   on   S u
ω = ω d   on   S u .
The unknowns are the displacement field u and the dislocation/disclination density fields ( α , θ ) . They are known at the initial time. The incompatible component κ p of κ p is solution of Equation (107), with the side conditions mentioned in Equation (40). ϵ p can be determined from Equation (108), with the side conditions div   ϵ p = 0 , ϵ p · n = 0 on S. Being history-dependent, κ p and ϵ p are initially set to zero without loss of generality. Furthermore, using the constitutive laws Equations (109) and (110), the balance of momentum and moment of momentum problem Equation (111) can be solved for the compatible displacement, as shown in Section 4.2. Then, the plastic strain rate and curvature rate ( ϵ ˙ p , κ ˙ p ) are found from Equations (112) and (113), where Equations (116)–(119) have been used to evaluate the dislocation/disclination velocities and the mesoscopic terms. Only the compatible parts ( ϵ ˙ p , κ ˙ p ) are employed for updating the plastic strain and curvature tensors. Finally, the dislocation/disclination densities are updated using Equations (114) and (115) and the procedure is iterated at the next time step.

7.2. Reduced Algorithm

Using the rate form of the model equations leads to a simpler and faster incremental algorithm. The main reason and difference is that the decomposition of plastic strains and curvatures into incompatible and compatible parts is possibly accounted for at the initial time, but not later. By taking the time derivative of Equations (109)–(111) and ignoring Equations (107) and (108), the algorithm uses the following equations:
0 = div   T ˙ s y m + 1 2 curl   div   M ˙ d e v
T ˙ = C : ( ϵ ˙ ϵ ˙ p ) + D : ( κ ˙ κ ˙ p )
M ˙ d e v = A : ( κ ˙ κ ˙ p ) + B : ( ϵ ˙ ϵ ˙ p )
κ ˙ p = θ × V θ + κ ˙ p *
ϵ ˙ p = ( α × V α ) s y m + ϵ ˙ p *
θ ˙ = curl   κ ˙ p
α ˙ = curl   ϵ ˙ p + κ ˙ p t t r ( κ ˙ p ) I
V α = 1 B α T s y m · α : X
V θ = 1 B θ M d e v , t · θ : X
ϵ ˙ p * = λ ϵ T s y m
κ ˙ p * = λ κ M d e v ,
with the boundary conditions:
T ˙ · n = t ˙ d     on   S t
u ˙ = u ˙ d   on   S u
ω ˙ = ω ˙ d   on   S u .
All fields ( u , α , θ ) are given at the initial time, as well as the plastic strain and curvature ( ϵ p , κ p ) arbitrarily set to zero without inconvenience. All fields can then be updated using the above rate equations. Clearly, the solutions obtained from these two algorithms may differ substantially because their continuity and differentiability properties are different.

7.3. Numerical Implementation

The finite element implementation couples a conventional Galerkin scheme for the solution of the equilibrium problem with a mixed Galerkin-Least Squares (GLS) scheme for the numerical approximation of the dislocation/disclination transport equations. Applying a standard Galerkin method to the convective transport equations leads to high frequency, propagating, numerical instabilities arising at sharp gradients. The Least-Squares (LS) finite element method, based on minimizing the norm of the residual, has the potential for eliminating such drawbacks, because it damps high frequency oscillations. Explicit LS methods have been considered so far [60,61]. However, damping is also a drawback of LS methods, because it leads to undue solution relaxation in the long-term. Therefore, an explicit Galerkin-Least-Squares (GLS) method involving convenient linear combinations of the Galerkin and LS algorithms has been proposed [60,61].
A constant drawback of finite element methods is that three-dimensional simulations are exceedingly demanding on computational resources. Spectral methods based on Fast Fourier Transform (FFT) algorithms were shown to be efficient at solving three-dimensional problems in periodic media and therefore appear as attractive alternatives to finite elements. However, the Gibbs oscillations they inherently produce in regions of strong spatial gradients reinforce the numerical instability of the solutions to the transport equations with respect to high frequency perturbations. Spectral low-pass filters in the Fourier space were shown to be efficient at damping these perturbations, and consequently at providing accurate and stable solutions to the transport problem [59]. FFT algorithms for the solution of the equilibrium problems are also available [55,63].

8. A Plane Edge-Wedge Model

As an illustration of the model and a basis for the next two Sections, we consider a bi-dimensional model with a distribution of wedge disclinations first discussed in [18]. In the orthonormal frame ( e 1 , e 2 , e 3 ) , we take the disclination tensor to be: θ = θ 33 e 3 e 3 . The continuity Equation (48) then implies that the disclination density θ 33 only depends on the coordinates ( x 1 , x 2 ) . Thus, the incompatibility Equation (46) reduces to:
θ 33 = κ 31 , 2 p , κ 32 , 1 p , ,
and the Poisson Equation (47) is:
κ 31 , 11 p , + κ 31 , 22 p , = + θ 33 , 2
κ 32 , 11 p , + κ 32 , 22 p , = θ 33 , 1 .
The involved plastic curvatures are only ( κ 31 p , κ 32 p ) and the disclination transport Equation (88) reads:
θ ˙ 33 = κ ˙ 31 , 2 p κ ˙ 32 , 1 p .
We ignore the mesoscopic plastic curvature rate κ ˙ p * in Equation (92) as we want to deal with nanoscale problems. Therefore, the plastic curvature rates are:
κ ˙ 31 p = θ 33 V 2 θ
κ ˙ 32 p = + θ 33 V 1 θ .
The constitutive relation Equation (104) for the disclination velocities then yields:
V 1 θ = + 1 B θ M 32 d e v θ 33
V 2 θ = 1 B θ M 31 d e v θ 33 ,
and the plastic curvature rates read:
κ ˙ 31 p = 1 B θ M 31 d e v θ 33 2
κ ˙ 32 p = 1 B θ M 32 d e v θ 33 2 ,
while the transport equation becomes:
θ ˙ 33 = 1 B θ ( M 31 d e v θ 33 2 ) , 2 1 B θ ( M 32 d e v θ 33 2 ) , 1 .
In the dislocation transport Equation (94), the source/sink term s θ generates only edge dislocations ( α 13 , α 23 ) . Ignoring also the mesoscopic plastic strain rate ϵ ˙ p * in Equation (95), the motion of these edge dislocations produces the plastic strain rates:
ϵ ˙ 11 p = α 13 V 2 α
ϵ ˙ 12 p = ϵ ˙ 21 p = 1 2 ( α 13 V 1 α α 23 V 2 α )
ϵ ˙ 22 p = + α 23 V 1 α ,
and the plastic rotation rate:
ω ˙ 3 p = 1 2 ( α 13 V 1 α + α 23 V 2 α ) .
The above relations Equations (148) and (150) show that climb of edge dislocations ( α 13 , α 23 ) produces the extension rates ( ϵ ˙ 11 p , ϵ ˙ 22 p ) , while their glide produces the shear ϵ ˙ 12 p in Equation (149). The dislocation transport Equation (94) therefore yields:
α ˙ 13 = ϵ ˙ 11 , 2 p ϵ ˙ 12 , 1 p + κ ˙ 31 p
α ˙ 23 = ϵ ˙ 21 , 2 p ϵ ˙ 22 , 1 p + κ ˙ 32 p .
Using Equations (141), (142), (148), (149) and (150), these edge dislocation generation rates read:
α ˙ 13 = ( α 13 V 2 α ) , 2 1 2 ( α 13 V 1 α α 23 V 2 α ) , 1 θ 33 V 2 θ
α ˙ 23 = ( α 23 V 1 α ) , 1 + 1 2 ( α 13 V 1 α α 23 V 2 α ) , 2 + θ 33 V 1 θ .
If all other dislocation densities are initially absent, the model therefore involves only α 13 and α 23 densities. The continuity condition Equation (49) then implies a planar distribution of these dislocations. The relation Equation (103) provides the dislocation velocities for climb:
V 1 α = + 1 B α T 22 s y m α 23
V 2 α = 1 B α T 11 s y m α 13 ,
and glide:
V 1 α = + 1 B α T 12 s y m α 13
V 2 α = 1 B α T 12 s y m α 23 .
The corresponding plastic strain rates are:
ϵ ˙ 11 p = 1 B α T 11 s y m α 13 2
ϵ ˙ 12 p = ϵ ˙ 21 p = 1 B α T 12 s y m ( α 13 2 + α 23 2 )
ϵ ˙ 22 p = 1 B α T 22 s y m α 23 2 ,
while the transport equation can be reformulated as:
α ˙ 13 = + 1 B α ( T 11 s y m α 13 2 ) , 2 1 2 B α ( T 12 s y m ( α 13 2 + α 23 2 ) ) , 1 + 1 B θ M 31 d e v θ 33 2
α ˙ 23 = 1 B α ( T 22 s y m α 23 2 ) , 1 + 1 2 B α ( T 12 s y m ( α 13 2 + α 23 2 ) ) , 2 + 1 B θ M 32 d e v θ 33 2 .
The relevant stress and deviatoric couple-stress components in the model are ( T 11 , T 12 , T 21 , T 22 ) and ( M 31 d e v , M 32 d e v ) . The balance of momentum Equation (56) read here:
T 11 , 1 s y m + T 12 , 2 s y m + 1 2 ( M 31 , 12 d e v + M 32 , 22 d e v ) = 0
T 21 , 1 s y m + T 22 , 2 s y m 1 2 ( M 31 , 11 d e v + M 32 , 21 d e v ) = 0 .
Thus, the loading applied at the boundaries of the domain may include shear/tension/ compression in the directions ( e 1 , e 2 ) . Couple stresses and rotations along e 3 can also be prescribed. The elastic constitutive laws Equations (82) and (83) yield:
T i j s y m = ψ ϵ i j e = C i j k l ϵ k l e + D i j k l κ k l e
M i j d e v = ψ κ i j e = A i j k l κ k l e + B i j k l ϵ k l e .
Substituting ( ϵ i j ϵ i j p , ω i , j κ i j p ) for ( ϵ i j e , κ i j e ) in the constitutive Equations (167) and (168), and introducing the resultant stresses and couple stresses in the balance Equations (165) and (166), one obtains two fourth-order partial differential equations for the material displacements ( u 1 , u 2 ) . As already mentioned, the boundary conditions comprise the prescription of tractions and moments, or/and displacements and rotations on the surface of the body. Please note that no boundary condition is specified on the disclination and dislocation densities.

9. Structure and Elastic Energy of Symmetric Tilt Boundaries

Copper bi-crystals with 37 ( 610 ) and 5 ( 310 ) symmetrical tilt boundaries of misorientation 18.9° and 36.9° were simulated in [64]. The bi-crystals are free of loads and periodic boundary conditions are imposed on the vertical boundaries in Figure 2. The characteristic size of the finite element mesh is 0.5Å. The wedge disclination density fields are shown in Figure 2, where they are superimposed on the relaxed atomic structure as predicted by atomistic calculations [14]. They consist of periodic arrays of wedge-disclination dipoles located at the edges of the structural units. Dilatations/contractions regions are associated with negative/positive disclinations. The local occurrence of large elastic strains suggests that a finite strain setting could be in order. Nevertheless, the predicted elastic fields are in good agreement with atomistic simulations and experiment. The elastic energy is essentially located in the dilatation/contraction areas in the core region of defects. In agreement with atomistic simulations, the energy is contained in the first atomic rows forming the grain boundary, essentially at structural units. Comparing the elastic fields of the 37 ( 610 ) boundary and of the 5 ( 310 ) boundary in Figure 2 shows that energy is more localized near the interface when the structural units are close to each other, due to elastic screening effects. Furthermore, the continuous description predicts features such as the shear deformation and curvature across/along the grain boundaries. As shown in Figure 3, there is a good agreement between quantitative estimates of the excess energy of grain boundaries using our model, and measurements from atomistic simulations and experiment [65,66,67]. Satisfactory agreement is obtained for all misorientations. Notably, the energy cusps predicted for the 5 ( 310 ) and 5 ( 210 ) grain boundaries originate in the self-screening of the elastic fields when the separation distance between structural units becomes small.

10. Shear-Coupled Migration of Symmetric Tilt Boundaries

As already suggested above, shear-induced normal grain boundary migration is a fascinating example of grain boundary-mediated plasticity. It was experimentally observed at moderate temperature in Al [68,69]. Following [17], we illustrate this mechanism by submitting the copper bicrystal featuring the 5 ( 310 ) symmetric <001> tilt boundary studied above to shear parallel to the boundary. In Figure 4, the wedge-disclination dipoles array is observed to move downward. This motion generates edge dislocation dipoles α 13 (line direction normal to the figure, Burgers vector along the boundary). Due to the shear stress T 12 these edge dipoles glide along the boundary, generate the plastic shear strain field ϵ 12 p seen in Figure 4c and annihilate at mid-distance, as Figure 4a,b show. Meanwhile, the downward migration of the disclination dipoles produces a plastic tilt rotation ω 3 p shown in Figure 4d. Striking similitude is observed between the heterogeneity of the plastic shear/rotation fields presented in Figure 4 and those obtained by post-processing atomistic simulations [16,70,71]. Because the boundary moves downward while positive shear is produced, the shear-coupling factor β , which corresponds to the ratio between the parallel translation of the two crystals and the normal migration of the boundary, is negative. The value β 0.94 measured in Figure 4b is in reasonable agreement with atomistic simulations and experiments [14,68]. The predicted variations of the shear-coupling factor in the full range of misorientations are presented in Figure 3e. The model retrieves the two <001> and <110> branches found in experiments [68] and obtained by molecular dynamics simulations and surface-dislocation-based models [13,14].
To further analyze the role of disclinations in plasticity, the model was used for simulating shear-coupled migration in forsterite. More specifically, we have used a periodic wedge-disclination array in the particular case of a (011)/<100> tilt boundary with misorientation 60°. This boundary was previously modeled at the atomic scale [72]. Figure 5 presents the disclination density field used to model the grain boundary structure [73]. Wedge-disclination densities are placed on the edges of structural units identified in [72]. The disclination arrangements result in a self-screening quadrupolar configuration, such that most of the energy is contained within the quadrupoles (see Figure 5b). The predicted excess energy of 1.3 J/m2 is in good quantitative agreement with the value 1.28 J/m2 measured from atomistic calculations [72]. When a shear stress is applied, the disclination array moves normal to the boundary, as shown in Figure 5. As detailed above for copper bi-crystals, the migration generates dislocation dipoles whose glide parallel to the boundary generates plastic shear within the region swept by the boundary. A more comprehensive study of the disclination structure, energy and mobility of grain boundaries as a function of their misorientation is clearly needed to eventually reach an accurate description of the rheology of olivine aggregates. EBSD-based disclination field imaging, as introduced in [6,7], applied to naturally deformed olivine aggregates should also be helpful in firmly establishing such models.

11. Concluding Remarks

The elasto-plastic model of dislocations and disclinations essentially consists of a set of partial differential equations, where the unknowns are the material displacement and dislocation/disclination density fields resolved at an atomic resolution length scale, and evolving at the time scale of the elementary dissipative mechanisms [18]. It is nonlocal in space and time in the conventional variables of compatible continuum mechanics and it provides a nonlocal incompatible generalization of the latter leading to well-posed problems in dislocation and disclination dynamics. It extends the elasto-plastic model of dislocations [22] by dealing, in addition to translational discontinuities and dislocations, with rotational discontinuities and disclinations. which allows grain boundary modeling in polycrystals. It also extends the elasto-static model of dislocations and disclinations [2], by dealing with plasticity mediated by transport of dislocations/disclinations within a consistent thermomechanical framework.
It has been sometimes believed that the assumption of continuity of the field variables is no longer valid at nanoscale and below, when the discrete atomistic nature of matter becomes apparent, and that continuum mechanics consequently fails to capture the physical phenomena at the scale of a few lattice spacings. In our opinion, the present results demonstrate instead that an incompatible, nonlocal continuum description involving fields of displacement and crystal defect densities that are smooth at this scale can be adequate for the purpose of describing such phenomena, in complement to atomistic representations. Indeed, the continuous disclination and dislocation theory proposed here was satisfactorily used to model the core structure of grain boundaries as well as the resulting elastic fields and excess energies. It was also shown to predict the complex interactions between disclinations/dislocations at work in the generation of plastic rotations/strains during shear-coupled boundary migration.
A continuous description of lattice defects may be more attractive than discrete approaches for various reasons. First, smoothness is desirable from the point of view of mathematical analysis and numerical computation, and because it allows dealing with defects’ core properties. When viewed at a small scale, lattice defects and the corresponding distributions of elastic strain and energy are certainly better described by smooth density fields than by singularities. Building on this point of view, recent work [16,70,71,74,75,76] has been devoted to the construction of field representations of the lattice deformation from the atomic displacements resulting from molecular dynamics simulations.
Secondly, such descriptions do not have to resolve atomic vibrations in the manner of atomistic simulations. The kinetic energy of atomic vibrations is time-averaged over periods of microseconds and rendered as dissipation, which results in higher computational efficiency. As a practical result, finite element or spectral FFT numerical approximations may allow considering the dynamics of crystal defect ensembles over time scales in the μ s or more, under realistic loading rates and stresses.
Finally, continuous approaches of crystal defects at nanoscale offer a convenient basis for seamless coarse-graining by using appropriate averaging techniques, which allows building mesoscale continuous theories of grain boundary-mediated plasticity in grain aggregates [56]. This ability could allow investigating alternative disclination-based mechanisms for plasticity, e.g., triple junction motion or emission/absorption of dislocations at grain boundaries, in media where dislocation-based plasticity is limited, either by scarcity of the slip systems, as in olivine, or because dislocation glide is restricted, as in polynanocrystals.

Author Contributions

Writing–review and editing, C.F. and V.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Volterra, V. Sur l’équilibre des corps élastiques multiplement connexes. Ann. Sci. Ecol. Norm. Supér. III 1907, 24, 401–517. [Google Scholar] [CrossRef] [Green Version]
  2. De Wit, R. Linear theory of static disclinations. In Fundamental Aspects of Dislocation Theory; Simmons, J.A., de Wit, R., Bullough, R., Eds.; National Bureau of Standards Special Publication: Washington, DC, USA, 1970; Volume 317, pp. 651–680. [Google Scholar]
  3. Friedel, J. Dislocations; Pergamon Press: Oxford, UK, 1964. [Google Scholar]
  4. Romanov, A.E.; Vladimirov, V.I. Disclinations in crystalline solids. In Dislocations in Solids; Nabarro, F.R.N., Ed.; Elsevier: Amsterdam, The Netherlans, 1992; Volume 9, p. 191. [Google Scholar]
  5. Romanov, A.E.; Kolesnikova, A.L. Application of disclination concept to solid structures. Prog. Mater. Sci. 2009, 54, 740–769. [Google Scholar] [CrossRef]
  6. Beausir, B.; Fressengeas, C. Disclination densities from EBSD orientation mapping. Int. J. Solids Struct. 2013, 50, 137–146. [Google Scholar] [CrossRef]
  7. Fressengeas, C.; Beausir, B. Tangential continuity of the curvature tensor at grain boundaries underpins disclination density determination from spatially mapped orientation data. Int. J. Solids Struct. 2019, 156–157, 210–215. [Google Scholar] [CrossRef]
  8. Bozhko, S.I.; Taupin, V.; Lebyodkin, M.; Fressengeas, C.; Levchenko, E.A.; Radikan, K.; Lübben, O.; Semenov, V.N.; Shvets, I.V. Disclinations in C60 molecular layers on WO2/W(110) surfaces. Phys. Rev. B 2014, 90, 214106. [Google Scholar] [CrossRef] [Green Version]
  9. Gertsman, V.Y.; Nazarov, A.A.; Romanov, A.E.; Valiev, R.Z.; Vladimirov, V.I. Disclination-structural unit model of grain boundaries. Philos. Mag. A 1989, 59, 1113–1118. [Google Scholar] [CrossRef]
  10. Hurtado, J.A.; Eliott, B.R.; Shodja, H.M.; Gorelikov, D.V.; Campbell, C.E.; Lippard, H.E.; Isabell, T.C.; Weertmann, J. Disclination grain boundary model with plastic deformation by dislocations. Mater. Sci. Eng. A 1995, 190, 1. [Google Scholar] [CrossRef]
  11. Li, J.C.M. Disclination model of high angle grain boundaries. Surf. Sci. 1972, 31, 12. [Google Scholar] [CrossRef]
  12. Shih, K.K.; Li, J.C.M. Energy of grain boundaries between cusp misorientations. Surf. Sci. 1975, 50, 109–124. [Google Scholar] [CrossRef]
  13. Berbenni, S.; Paliwal, B.; Cherkaoui, M. A micromechanical-based model for shear-coupled grain boundary migration in bicrystals. Int. J. Plast. 2013, 44, 68–94. [Google Scholar] [CrossRef] [Green Version]
  14. Cahn, J.W.; Mishin, Y.; Suzuki, A. Coupling grain boundary motion to shear deformation. Acta Mater. 2006, 54, 4953–4975. [Google Scholar] [CrossRef]
  15. Farkas, D.; Froseth, A.; Van Swygenhoven, H. Grain boundary migration during room temperature deformation of nanocrystalline Ni. Scr. Mater. 2006, 55, 695–698. [Google Scholar] [CrossRef]
  16. Tucker, G.J.; Zimmerman, J.A.; McDowell, D.L. Shear deformation kinematics of bicrystalline grain boundaries in atomistic simulations. Model. Simul. Mater. Sci. Eng. 2010, 18, 015002. [Google Scholar] [CrossRef]
  17. Taupin, V.; Capolungo, L.; Fressengeas, C. Disclination mediated plasticity in shear-coupled boundary migration. Int. J. Plast. 2014, 53, 179–192. [Google Scholar] [CrossRef]
  18. Fressengeas, C.; Taupin, V.; Capolungo, L. An elasto-plastic theory of dislocation and disclination fields. Int. J. Solids Struct. 2011, 48, 3499–3509. [Google Scholar] [CrossRef] [Green Version]
  19. Acharya, A.; Fressengeas, C. Continuum mechanics of the interactions between phase boundaries and dislocations in solids. In Differential Geometry and Continuum Mechanics; Chen, G.Q., Grinfeld, M., Knops, R.J., Eds.; Springer: Berlin/Heidelberg, Germany, 2015; Volume 137, pp. 125–168. [Google Scholar]
  20. Kröner, E. Kontinuumstheorie der Versetzungen und Eigenspannungen (Ergebnisse der Angewandten Mathematik, No. 5); Collatz, L., Lösch, F., Eds.; Springer: Berlin/Heidelberg, Germany, 1958. [Google Scholar]
  21. Kröner, E. Continuum theory of defects. In Physics of Defects; North Holland: Amsterdam, The Netherlands, 1980; pp. 218–314. [Google Scholar]
  22. Acharya, A. A model of crystal plasticity based on the theory of continuously distributed dislocations. J. Mech. Phys. Solids 2001, 49, 761. [Google Scholar] [CrossRef]
  23. Acharya, A. Driving forces and boundary conditions in continuum dislocation mechanics. Proc. R. Soc. Lond. 2003, A459, 1343–1363. [Google Scholar] [CrossRef] [Green Version]
  24. Puri, S.; Das, A.; Acharya, A. Mechanical response of multicrystalline thin films in mesoscale field dislocation mechanics. J. Mech. Phys. Solids 2011, 59, 2400–2417. [Google Scholar] [CrossRef]
  25. Taupin, V.; Varadhan, S.; Chevy, J.; Fressengeas, C.; Beaudoin, A.J.; Montagnat, M.; Duval, P. Effects of size on the dynamics of dislocations in ice single crystals. Phys. Rev. Lett. 2007, 99, 155507. [Google Scholar] [CrossRef]
  26. Taupin, V.; Varadhan, S.; Fressengeas, C.; Beaudoin, A.J. Directionality of yield point in strain-aged steels: The role of polar dislocations. Acta Mater. 2008, 56, 3002–3010. [Google Scholar] [CrossRef]
  27. Varadhan, S.; Beaudoin, A.J.; Fressengeas, C. Lattice incompatibility and strain-aging in single crystals. J. Mech. Phys. Solids 2009, 57, 1733–1748. [Google Scholar] [CrossRef]
  28. Gbemou, K.; Taupin, V.; Raulot, J.M.; Fressengeas, C. Building compact dislocation cores in an elasto-plastic model of dislocation fields. Int. J. Plast. 2016, 82, 241–259. [Google Scholar] [CrossRef]
  29. Zhang, X.; Acharya, A.; Walkington, N.J.; Bielak, J. A single theory for some quasi-static, supersonic, atomic, and tectonic scale applications of dislocations. J. Mech. Phys. Solids 2015, 84, 145–195. [Google Scholar] [CrossRef] [Green Version]
  30. Chevy, J.; Fressengeas, C.; Lebyodkin, M.; Taupin, V.; Bastie, P.; Duval, P. Characterizing short range vs. long range spatial correlations in dislocation distributions. Acta Mater. 2010, 58, 1837–1849. [Google Scholar] [CrossRef]
  31. Fressengeas, C.; Beaudoin, A.J.; Entemeyer, D.; Lebedkina, T.; Lebyodkin, M.; Taupin, V. Dislocation transport and intermittency in the plasticity of crystalline solids. Phys. Rev. B 2009, 79, 014108. [Google Scholar] [CrossRef] [Green Version]
  32. Kossecka, E.; de Wit, R. Disclination kinematics. Arch. Mech. 1977, 29, 633–651. [Google Scholar]
  33. Kossecka, E.; de Wit, R. Disclination dynamics. Arch. Mech. 1977, 29, 749–767. [Google Scholar]
  34. Das, E.S.P.; Marcinkowski, M.J.; Armstrong, R.W.; de Wit, R. The movement of Volterra disclinations and the associated mechanical forces. Philos. Mag. 1973, 27, 369–391. [Google Scholar] [CrossRef]
  35. Pauli, W. Handbuch der Physik, Band XXIV, Teil 1, Vol. V, Part 1; Springer: Berlin, Germany, 1933. [Google Scholar]
  36. Schrödinger, E. The energy-impulse hypothesis of material waves. Ann. Phys. 1927, 82, 265. [Google Scholar] [CrossRef]
  37. Nielsen, O.H.; Martin, R.M. Quantum-mechanical theory of stress and force. Phys. Rev. B 1985, 32, 3780–3791. [Google Scholar] [CrossRef]
  38. Jiang, B. The Least-Squares Finite Element Method; Springer: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  39. Acharya, A.; Roy, A. Size effects and idealized dislocation microstructure at small scales: Predictions of a phenomenological model of Mesoscopic Field Dislocation Mechanics: Part I. J. Mech. Phys. Solids 2006, 54, 1687–1710. [Google Scholar] [CrossRef] [Green Version]
  40. Weingarten, J. Sulle superficie di discontinuità nella teoria della elasticità dei corpi solidi. Rend. R. Accad. Lincei 1901, 10, 57–60. [Google Scholar]
  41. Mindlin, R.D.; Tiersten, H.F. Effects of couple-stresses in linear elasticity. Arch. Ratlat. Mech. Anal. 1962, 11, 415–448. [Google Scholar] [CrossRef]
  42. Eringen, A.C. Theory of micropolar elasticity. In Fracture; Liebowitz, H., Ed.; Academic Press: New York, NY, USA, 1968; Volume 2, pp. 621–729. [Google Scholar]
  43. Koiter, W.T. Couple stresses in the theory of elasticity. I II Proc. Ned. Akad. Wet. 1964, B, 17–44. [Google Scholar]
  44. Fressengeas, C.; Sun, X. On the theory of dislocations and generalized disclinations and its application to straight and stepped symmetrical tilt boundaries. J. Mech. Phys. Solids 2020, 143, 104092. [Google Scholar] [CrossRef]
  45. Eringen, A.C. A unified theory of thermomechanical materials. Int. J. Eng. Sci. 1966, 4, 179–202. [Google Scholar] [CrossRef]
  46. Eringen, A.C.; Edelen, D.G.B. Nonlocal elasticity. Int. J. Eng. Sci. 1972, 10, 233–248. [Google Scholar] [CrossRef]
  47. Eringen, A.C.; Kim, B.S. Relation between nonlocal elasticity and lattice dynamics. Cryst. Lattice Defects 1977, 7, 51–57. [Google Scholar]
  48. Eringen, A.C. Continuum Mechanics at the Atomic Scale; No. 77-SM-1; Princeton University: Princeton, NJ, USA, 1977. [Google Scholar]
  49. Eringen, A.C. Non Local Continuum Field Theories; Springer: New York, NY, USA, 2002. [Google Scholar]
  50. Dingreville, R.; Hallil, A.; Berbenni, S. From coherent to incoherent mismatched interfaces: A generalized continuum formulation of surface stresses. J. Mech. Phys. Solids 2014, 72, 40–60. [Google Scholar] [CrossRef]
  51. Spearot, D.E.; Jacob, K.I.; McDowell, D.L. Dislocation nucleation from bicrystal interfaces with dissociated structure. Int. J. Plast. 2007, 23, 143–160. [Google Scholar] [CrossRef]
  52. Spearot, D.E.; Capolungo, L.; Qu, J.; Cherkaoui, M. On the elastic tensile deformation of <100> bicrystal interfaces in copper. Comput. Mater. Sci. 2008, 42, 57–67. [Google Scholar]
  53. Tschopp, M.A.; Tucker, G.J.; McDowell, D.L. Atomistic simulations of tension-compression asymmetry in dislocation nucleation for copper grain boundaries. Comput. Mater. Sci. 2008, 44, 351–362. [Google Scholar] [CrossRef]
  54. Upadhyay, M.; Capolungo, L.; Taupin, V.; Fressengeas, C. Elastic constitutive laws for incompatible crystalline media: The contributions of dislocations, disclinations and G-disclinations. Philos. Mag. 2013, 93, 794–832. [Google Scholar] [CrossRef]
  55. Taupin, V.; Gbemou, K.; Fressengeas, C.; Capolungo, L. Nonlocal elasticity tensors in dislocation and disclination cores. J. Mech. Phys. Solids 2017, 100, 62–84. [Google Scholar] [CrossRef]
  56. Taupin, V.; Capolungo, L.; Fressengeas, C.; Upadhyay, M.; Beausir, B. A mesoscopic theory of dislocation and disclination fields for grain boundary-mediated crystal plasticity. Int. J. Solids Struct. 2015, 71, 277–290. [Google Scholar] [CrossRef]
  57. Mura, T. Continuous distributions of moving dislocations. Philos. Mag. 1963, 89, 843–857. [Google Scholar] [CrossRef]
  58. Kleman, M.; Friedel, J. Disclinations, dislocations, and continuous defects: A reappraisal. Rev. Mod. Phys. 2008, 80, 61–115. [Google Scholar] [CrossRef] [Green Version]
  59. Djaka, K.S.; Taupin, V.; Berbenni, S.; Fressengeas, C. A numerical spectral approach to solve the dislocation density transport equation. Model. Simul. Mater. Sci. Eng. 2015, 23, 065008. [Google Scholar] [CrossRef]
  60. Roy, A.; Acharya, A. Finite element approximation of field dislocation mechanics. J. Mech. Phys. Solids 2005, 53, 143–170. [Google Scholar] [CrossRef]
  61. Varadhan, S.; Beaudoin, A.J.; Acharya, A.; Fressengeas, C. Dislocation transport using an explicit Galerkin/least-squares formulation. Model. Simul. Mater. Sci. Eng. 2006, 14, 1245–1270. [Google Scholar] [CrossRef]
  62. Coleman, B.D.; Gurtin, M.E. Thermodynamics with internal state variables. J. Chem. Phys. 1967, 47, 597–613. [Google Scholar] [CrossRef]
  63. Djaka, K.S.; Villani, A.; Taupin, V.; Capolungo, L.; Berbenni, S. Field Dislocation Mechanics for heterogeneous elastic materials: A numerical spectral approach. Comput. Methods Appl. Mech. Eng. 2017, 315, 921–942. [Google Scholar] [CrossRef] [Green Version]
  64. Fressengeas, C.; Taupin, V.; Capolungo, L. Continuous modeling of the structure of symmetric tilt boundaries. Int. J. Solids Struct. 2014, 51, 1434–1441. [Google Scholar] [CrossRef]
  65. Hasson, G.; Boos, J.Y.; Herbeuval, I.; Biscondi, M.; Goux, C. Theoretical and experimental determinations of grain-boundary structures and energies—Correlations with various experimental results. Surf. Sci. 1972, 31, 115–137. [Google Scholar] [CrossRef]
  66. Tschopp, M.A.; Spearot, D.E.; McDowell, D.L. Influence of grain boundary structure on dislocation nucleation in FCC metals. In Dislocations in Solids; Hirth, J.P., Ed.; North Holland Publishing: Amsterdam, The Netherlands, 2008; Volume 14, pp. 43–139. [Google Scholar]
  67. Bachurin, D.V.; Murzaev, R.T.; Nazarov, A.A. Atomistic computer and disclination simulation of [001] tilt boundaries in nickel and copper. Phys. Met. Metallogr. 2003, 96, 555–561. [Google Scholar]
  68. Gorkaya, T.; Molodov, D.A.; Gottstein, G. Stress-driven migration of symmetric <001> tilt grain boundaries in Al bicrystals. Acta Mater. 2009, 57, 5396–5405. [Google Scholar]
  69. Mompiou, F.; Caillard, D.; Legros, M. Grain boundary shear-migration coupling—I. In situ TEM straining experiments in Al polycrystals. Acta Mater. 2009, 57, 2198–2209. [Google Scholar] [CrossRef]
  70. Tucker, G.J.; Tiwari, S.; Zimmerman, J.A.; McDowell, D.L. Investigating the deformation of nanocrystalline copper with microscale kinematic metrics and molecular dynamics. J. Mech. Phys. Solids 2012, 60, 471–486. [Google Scholar] [CrossRef]
  71. Zimmerman, J.A.; Bammann, D.J.; Gao, H. Deformation gradients for continuum mechanical analysis of atomistic simulations. Int. J. Solids Struct. 2009, 46, 238–253. [Google Scholar] [CrossRef] [Green Version]
  72. Adjaoud, O.; Marquardt, K.; Jahn, S. Atomic structures and energies of grain boundaries in Mg2SiO4 forsterite from atomistic modeling. Phys. Chem. Miner. 2012, 39, 749–760. [Google Scholar] [CrossRef]
  73. Cordier, P.; Demouchy, S.; Beausir, B.; Taupin, V.; Barou, F.; Fressengeas, C. Disclinations provide the missing mechanism for deforming olivine-rich rocks in the mantle. Nature 2014, 507, 51–56. [Google Scholar] [CrossRef] [PubMed]
  74. Sun, X.Y.; Taupin, V.; Fressengeas, C.; Cordier, P. Continuous description of the atomic structure of a grain boundary using dislocations and generalized disclination density fields. Int. J. Plast. 2016, 77, 75–89. [Google Scholar] [CrossRef]
  75. Sun, X.Y.; Cordier, P.; Taupin, V.; Fressengeas, C.; Jahn, S. Continuous description of a grain boundary in forsterite from atomic scale simulations: The role of disclinations. Philos. Mag. A 2016, 96, 1757–1772. [Google Scholar] [CrossRef] [Green Version]
  76. Sun, X.Y.; Fressengeas, C.; Taupin, V.; Cordier, P.; Combe, N. Disconnections, dislocations and generalized disclinations in grain boundary ledges. Int. J. Plast. 2018, 104, 134–146. [Google Scholar] [CrossRef]
Figure 1. Nonlocal strains and curvatures. (A) Nonlocal shear strain at location r induced by the rotation d ω e arising from a curvature tensor at location r ; (B) Nonlocal rotation and curvature at location r from the displacement du induced by a shear strain tensor at location r .
Figure 1. Nonlocal strains and curvatures. (A) Nonlocal shear strain at location r induced by the rotation d ω e arising from a curvature tensor at location r ; (B) Nonlocal rotation and curvature at location r from the displacement du induced by a shear strain tensor at location r .
Metals 10 01517 g001
Figure 2. Disclination model for a 37 ( 610 ) symmetrical <001> tilt boundary of misorientation 18.9° (left column) and of a 5 ( 310 ) symmetrical <001> tilt boundary of misorientation 36.9° (right column). The wedge-disclination density θ 33 (rad/m2) is shown on top of the relaxed atomic structure (a,f); Elastic dilatations/contractions (b,g); Shear stress field T 12 (Pa) (c,h); Elastic curvature κ 32 (rad/m) (d,i); Elastic energy density (J/m3) (e,j). Atomic structures are taken from reference [14].
Figure 2. Disclination model for a 37 ( 610 ) symmetrical <001> tilt boundary of misorientation 18.9° (left column) and of a 5 ( 310 ) symmetrical <001> tilt boundary of misorientation 36.9° (right column). The wedge-disclination density θ 33 (rad/m2) is shown on top of the relaxed atomic structure (a,f); Elastic dilatations/contractions (b,g); Shear stress field T 12 (Pa) (c,h); Elastic curvature κ 32 (rad/m) (d,i); Elastic energy density (J/m3) (e,j). Atomic structures are taken from reference [14].
Metals 10 01517 g002
Figure 3. Excess energy of <001> tilt grain boundaries in copper bi-crystals. Black diamonds: present work; Gray triangles/squares: experimental measurements/molecular statics predictions [65]; Gray circles: molecular statics predictions [66].
Figure 3. Excess energy of <001> tilt grain boundaries in copper bi-crystals. Black diamonds: present work; Gray triangles/squares: experimental measurements/molecular statics predictions [65]; Gray circles: molecular statics predictions [66].
Metals 10 01517 g003
Figure 4. Disclination model for shear-coupled boundary migration. (a,b): Wedge-disclination density θ 33 (rad/m2) (color coded) and edge dislocation density α 13 (m−1) (black/white contours), during migration of a 5 ( 310 ) symmetric <001> tilt grain boundary of misorientation 36.87°. Two snapshots are shown at 1 ms (a) and 10 ms (b); (c,d): Plastic shear strain ϵ 12 p (c) and plastic rotation ω 3 p ( ° ) (d) generated during the migration. In (ad), the black dashed-line shows the initial position of the grain boundary. (e): Shear-coupling factor as a function of grain boundary misorientation. The circles show values obtained with the disclination model while the solid lines show the <001> and <110> shear deformation modes proposed in [14].
Figure 4. Disclination model for shear-coupled boundary migration. (a,b): Wedge-disclination density θ 33 (rad/m2) (color coded) and edge dislocation density α 13 (m−1) (black/white contours), during migration of a 5 ( 310 ) symmetric <001> tilt grain boundary of misorientation 36.87°. Two snapshots are shown at 1 ms (a) and 10 ms (b); (c,d): Plastic shear strain ϵ 12 p (c) and plastic rotation ω 3 p ( ° ) (d) generated during the migration. In (ad), the black dashed-line shows the initial position of the grain boundary. (e): Shear-coupling factor as a function of grain boundary misorientation. The circles show values obtained with the disclination model while the solid lines show the <001> and <110> shear deformation modes proposed in [14].
Metals 10 01517 g004
Figure 5. Disclination model for a (011)/<100>, 60° tilt boundary in forsterite. (a) Initial wedge-disclination density (rad·m−2) superimposed on top of the atomic structure. Triangles show the structural units; (b) Elastic energy (J/m3); (c) Migration of the (011)/<100> tilt boundary. The wedge-disclination density θ 33 (rad·m−2) is shown in color. The plastic shear generated by migration is shown with gray contours.
Figure 5. Disclination model for a (011)/<100>, 60° tilt boundary in forsterite. (a) Initial wedge-disclination density (rad·m−2) superimposed on top of the atomic structure. Triangles show the structural units; (b) Elastic energy (J/m3); (c) Migration of the (011)/<100> tilt boundary. The wedge-disclination density θ 33 (rad·m−2) is shown in color. The plastic shear generated by migration is shown with gray contours.
Metals 10 01517 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fressengeas, C.; Taupin, V. Revisiting the Application of Field Dislocation and Disclination Mechanics to Grain Boundaries. Metals 2020, 10, 1517. https://doi.org/10.3390/met10111517

AMA Style

Fressengeas C, Taupin V. Revisiting the Application of Field Dislocation and Disclination Mechanics to Grain Boundaries. Metals. 2020; 10(11):1517. https://doi.org/10.3390/met10111517

Chicago/Turabian Style

Fressengeas, Claude, and Vincent Taupin. 2020. "Revisiting the Application of Field Dislocation and Disclination Mechanics to Grain Boundaries" Metals 10, no. 11: 1517. https://doi.org/10.3390/met10111517

APA Style

Fressengeas, C., & Taupin, V. (2020). Revisiting the Application of Field Dislocation and Disclination Mechanics to Grain Boundaries. Metals, 10(11), 1517. https://doi.org/10.3390/met10111517

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop