Next Article in Journal
An Enhanced Oxidation of Formate on PtNi/Ni Foam Catalyst in an Alkaline Medium
Previous Article in Journal
Oxidation and Microstructural Behaviors at 1200 °C of 32.5 wt.% Cr–Containing Co–Based Alloys Strengthened by HfC Carbides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elastodynamics Field of Non-Uniformly Moving Dislocation: From 3D to 2D

Applied Mechanics Laboratory, School of Aerospace Engineering, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Crystals 2022, 12(3), 363; https://doi.org/10.3390/cryst12030363
Submission received: 24 January 2022 / Revised: 3 March 2022 / Accepted: 3 March 2022 / Published: 8 March 2022
(This article belongs to the Section Crystalline Metals and Alloys)

Abstract

:
Molecular dynamics (MD) and experiments indicate that the high-speed dislocations dominate the plasticity properties of crystal materials under high strain rate. New physical features arise accompanied with the increase in dislocation speed, such as the “Lorentz contraction” effect of moving screw dislocation, anomalous nucleation, and annihilation in dislocation interaction. The static description of the dislocation is no longer applicable. The elastodynamics fields of non-uniformly moving dislocation are significantly temporal and spatially coupled. The corresponding mathematical formulas of the stress fields of three-dimensional (3D) and two-dimensional (2D) dislocations look quite different. To clarify these differences, we disclose the physical origin of their connections, which is inherently associated with different temporal and spatial decoupling strategies through the 2D and 3D elastodynamics Green tensor. In this work, the fundamental relationship between 2D and 3D dislocation elastodynamics is established, which has enlightening significance for establishing general high-speed dislocation theory, developing a numerical calculation method based on dislocation elastodynamics, and revealing more influences of dislocation on the macroscopic properties of materials.

1. Introduction

High strain rate deformation has long been a central concern in the fields such as machine manufacturing, automotive, aviation, aerospace, and military [1,2,3]. Under such loading conditions accompanied by extremely high temperature and pressure, there are many micromechanisms of plastic deformation for materials, such as dislocation and deformation twins. Complex phenomena can also arise; for example, a void can lead to high-pressure phase transformation [4], and the barrierless melt nucleation will arise for some energetic nitramine [5]. However, numerous molecular dynamics (MD) simulations have demonstrated that the dislocation density will go up with the increase of impact pressure, while the nucleation and motion of dislocations still dominate the plastic deformation of the material under the high strain rate condition [6,7,8,9]. For example, when the strain rate is below 2.77 × 10 8 /s 1 and temperature is over 50 K, the plasticity is dislocation-mediated in the body-centered-cubic (BCC) metal tantalum [10]. Their simulations predict that only on reaching certain limiting conditions of stress (the uniaxial stress reaches about 8 GPa) will deformation twinning become the dominant mode of the dynamic response. Therefore, the strength and plasticity properties of a metal are dominated by dislocations in the vast majority of high strain rate deformation conditions.
According to Orowan’s equation, plastic strain rate ε ˙ = ρ b v ( ρ is the mobile dislocation density, b is the magnitude of Burger’s vector, and v is the average dislocation velocity). If ρ = 10 13 / m 2 , b = 0.25 nm, and ε ˙ = 5 × 10 6 /s, v is expected to be above 2 km/s, which is closed to the shear wave speed of some metals. For high-speed dislocation, its elastodynamics stress fields are significantly different from the static stress field of dislocations, which has been recognized since 1949 [11,12,13].
On one hand, due to the insufficient spatial and temporal resolution of experimental apparatus, it is very difficult to directly observe high-speed dislocation under such high strain rate conditions. The observation of the high-speed dislocation has extremely high requirements for laboratory equipment and samples, so far just having one report about the high-speed dislocation in a two-dimensional plasma crystal [14]. Measuring the high-speed dislocation in more materials still is very challenging and has not been successful. In order to study high-speed dislocation, most work about the 3D and 2D elastodynamics field is based on theoretical analysis, numerical calculation [15,16,17,18,19,20], and molecular dynamics simulations [21,22,23,24,25].
On the other hand, the underlying physics of the high-speed dislocation is substantially different from that of the quasi-static dislocations, and the high-speed dislocation is manifested in the distinctive elastodynamics field. For instance, the displacement field solution of the uniformly moving screw dislocation differs from the static solution by a “Lorentz contraction” [26]. With the increase of dislocation velocity, the stress field is contracted to the axis, perpendicular to the slip plane, on which the dislocation is located. Moreover, molecular dynamics (MD) simulations show that high-speed dislocation can enter the transonic and supersonic regimes [23,27]. Correspondingly, the special stress field of Mach cones appears [28,29,30,31]. In addition, the elastodynamics field is sensitive to the motion history of the dislocation. The elastodynamics stress fields of a uniformly expanding and shrinking dislocation loop are remarkably different, even though the final radius is both of 100b, where b is the magnitude of Burgers vector, and its velocity magnitude is the same as 0.3 c T ( c T is the shear wave speed of the material) [15].
The elastodynamics effect leads to distinct dislocation interactions compared with the classical quasi-static dislocations. When the dislocations move at the low velocity, two coplanar opposite-sign colliding edge or screw dislocations annihilate. As the dislocation velocity increases close to the shear wave speed, those two dislocations, moving against each other, fully overshoot each other, even generating new dislocations on near planes [22,24,32]. In addition, the annihilation of a dislocation from the free surface is generally observed when it moves with a low velocity. Unexpectedly, the rebounding of a edge or a screw dislocation from the free surface has been found when it moves at a high-velocity [33]. These abnormal mechanisms can lead to a rapid increase in dislocation density that is strongly correlated both spatially and temporally, affecting high-strain-rate plastic deformation from the microscopic level. However, to this day, the elastodynamics effect on high-speed dislocation is far from being well understood, because the full history-dependent nature and spatial–temporal coupling feature significantly complicate the analysis of the problem.
Among the elastodynamics effect on high-speed dislocation, a fundamental problem is that the elastodynamics solution of high-speed dislocation looks very different for two-dimensional (2D) and three-dimensional (3D) cases. Most previous focus has been put on the 2D infinitely long straight pure screw or edge high-speed dislocations [20,31,34]. For 2D screw dislocation, it is derived based on the assumption of anti-plane strain framework, and for 2D edge dislocation, it is derived based on the plane strain framework. For the 3D case, a general theory of non-singular dislocations loop is also developed for homogeneous and isotropic infinite material [15,35]. After the concept of dislocation appeared, the theory of quasi-static dislocations has been well developed [26]. The stress fields of the high-speed dislocation are substantially different from those of the quasi-static dislocations, which have been reported since 1949 [11,12]. From then on, dislocation elastodynamics theories were further developed [13]. In the last ten years, the fundamental time-dependent solutions for the injection and non-uniform motion of straight edge and screw dislocations are widely studied. The two-dimensional elastodynamics Green tensor is used to obtain the elastodynamics field of the subsonic, transonic, or supersonic dislocation motion [20]. Refs. [13,31] have made several important contributions for elastodynamics solutions for the non-uniform motion of 3D dislocations. To complete the 3D elastodynamics framework, Cui et al. have built explicit compact formula of high-speed 3D dislocations with non-uniform motion, which can be directly used for the numerical implementation [15]. The corresponding mathematical formula of the stress fields of three-dimensional (3D) and two-dimensional (2D) dislocations looks quite different (see Equations (9), (10), and (16), which will be discussed in detail later).
However, few scholars have paid attention to the connection and physical differences between the 3D and 2D of non-uniformly moving dislocation’s elastodynamics field, even though it is important to connect the 2D and 3D underlying mechanisms. To our knowledge, only for quasi-static dislocation, Lazar built up the 2D elastic fields from 3D elastostatics solution for anisotropic elasticity in Fourier space [36]. However, for high-speed dislocation, how to directly obtain the 2D solution from the 3D elastodynamics solution is seldomly discussed. In the present work, 3D and 2D elastodynamics Green tensors act as mediums for correlations between both elastodynamics fields, the 2D dislocation elastodynamics is successfully derived from the higher dimensional elastodynamics fields, and their intrinsic physical connections and distinctions are investigated.

2. Eastodynamics Solution of Non-Uniformly Moving Dislocation

The elastodynamics stress field of a 3D non-uniformly moving dislocation loop is as follows:
σ i j ( x , t ) = C i j k l β k l E ( x , t ) = t 0 t d t L ( t ) ϵ r q l C i j k l C m n p q G k m , n ( x x , t t ) b p ξ r ( x , t ) d L t 0 t d t L ( t ) C i j k l ρ G ˙ k p ( x x , t t ) b p ϵ l s r V s ( x , t ) ξ r ( x , t ) d L
where C i j k l is the elastic modulus tensor, β k l E is the elastic distortion solution, ρ is the material density, G k p is the elastodynamics Green tensor at the 3D, the subscript using Greek alphabet ranges from 1 to 3 for 3D case. ( ) , n denotes differentiation with respect to x n . A superposed dot means a time derivative. t 0 and t represent the initial and current time, x is a position vector of the observer point. x is the integration point position at time t along the dislocation line L . b is the Burgers vector, V is the dislocation velocity vector, and ξ is the unit tangent direction vector of the dislocation line. The deriviation of Equation (1) is referred to [13,15,37,38,39].
Equation (1) clearly show that the elastodynamics solution is spatial–temporal coupled and depends on the full history of the dislocation motion. In order to obtain the 3D and 2D solution in the same framework, the first key problem is to obtain the corresponding elastodynamics Green tensor and its degradation for 2D problems.

3. Degradation of Elastodynamics Green Tensor

This section will demonstrate how to degrade the 3D elastodynamics Green tensor to the 2D counterpart, corresponding to the problems of infinitely long straight dislocation lines.
The 3D elastodynamics Green’s function in an isotropic crystal is derived in [40,41],
G i j ( x , t ) = 1 4 π ρ x 1 c T 2 δ i j x i x j x 2 δ t x c T + 1 c L 2 x i x j x 2 δ t x c L + 3 x i x j x 2 δ i j 1 x 2 x / c L x / c T τ δ t τ d τ
where x and x i represent the magnitude and the ith component of vector x , respectively. δ i j is the Kronecker delta. c T and c L are the transverse and longitudinal wave speed, respectively.
Without loss of generality, an infinite-long straight dislocation is assumed to be along the x 3 direction (see Figure 1a). Therefore, the degradation of elastodynamics Green tensor needs to integrate Equation (2) for x 3 from to . c is used to express a specific elastic wave velocity c, which may be c T , c L , or the elastic wave speed between c T and c L . Then, the basic term in Equation (2) can be expressed as f ( x , t ) δ t x c . Its integration can be calculated as follows,
f ( x , t ) δ t x c d x 3 = f ( x , t ) 1 c d x d x 3 | x = c t = 2 f ( x , t ) 1 c x 3 x | x = c t = 2 c t f ( x , t ) t 2 l 2 / c 2 | x = c t
where l 2 = x 1 2 + x 2 2 . According to Equation (3), many basic equations can be obtained to calculate the 2D elastodynamics Green tensor and its time derivative, as listed in Appendix A. Combining Appendix A and Equation (2), the 2D elastodynamics Green tensor becomes,
G α β ( x , t ) = 1 2 π ρ δ α β l 2 t 2 H ( t 2 l 2 / c T 2 ) t 2 l 2 / c T 2 t 2 l 2 / c L 2 H ( t 2 l 2 / c L 2 ) + x α x β l 4 ( 2 t 2 + l 2 / c T 2 ) H ( t 2 l 2 / c T 2 ) t 2 l 2 / c T 2 + ( 2 t 2 l 2 / c L 2 ) H ( t 2 l 2 / c L 2 ) t 2 l 2 / c L 2
G 33 ( x , t ) = 1 2 π μ t 2 l 2 / c T 2
G α 3 ( x , t ) = 0
where the Roman subscripts α and β range from 1 to 2.
For isotropic material, C i j k l = λ δ i j δ k l + μ ( δ i k δ j l + δ i l δ j k ) , where λ and μ are Lamé’s first parameter and shear modulus, respectively. According to Equation (1),
β k l E ( x , t ) = t 0 t d t L ( t ) [ λ ϵ r q l G k m , m ( s , t t ) b q ξ r ( s , t ) + μ ϵ r q l G k p , q ( s , t t ) b p ξ r ( s , t ) + μ ϵ r q l G k q , p ( s , t t ) b p ξ r ( s , t ) ] d L t 0 t d t L ( t ) ρ G ˙ k p ( s , t t ) b p ϵ l s r V s ( s , t ) ξ r ( s , t ) d L .
For screw dislocation, the dislocation line vector ξ = [ 001 ] , so the Burger vector b = [ 00 b ] . It is easy to see from Equation (7) that ϵ r q l b q ξ r = 0 , so only G k 3 , p , G k p , 3 , and G ˙ k 3 , t contribute. Since G α β , 3 = 0 , G α 3 = 0 , only G 33 plays a role.
With respect to edge dislocation, the dislocation line vector ξ = [ 001 ] , and the Burger vector b = [ b 00 ] . According to Equation (7), only G k m , m , G k 1 , q , G k q , 1 , and G ˙ k 1 contribute, where subscript q cannot be 3. Since G 31 = 0 , G α 3 = 0 , and G 33 , 3 = 0 , only G α β plays a role.
The obtained 2D Green tensor is consistent with previous work, which is derived for a plane-strain and anti-plane strain problem with respect to edge dislocation and screw dislocation, respectively [31].

4. Different Spatial–Temporal Decoupling Strategies at 2D and 3D

4.1. Spatial–Temporal Decoupling Strategy at 2D

The degradation of elastodynamics Green tensor indicates that for an infinitely long straight dislocation in infinitely extended material, the spatial–temporal coupling integration is firstly done on a spatial scale, because the geometry of the dislocation line remains unchanged, as indicated by the dislocation configuration at different times ( t 0 , t 1 , t 2 , t) in Figure 1a. When the dislocations are assumed to be always straight, along x 3 all the time, the intersection point of the dislocation line and plane O x 1 x 2 can well describe the position of dislocation at different times. The straight line integrals along the dislocation are clear and definite at any time. Depending on the atomic displacement caused by the screw and edge dislocation, it will be the anti-plane strain problem of screw dislocation and plane strain problem of edge dislocation. Therefore, the 2D treatment is able to represent the 3D information of the infinitely long straight dislocation, which is similar to previous work, such as [17,31]. Afterward, this problem can be simplified as a 2D problem. As described in Figure 1b, each infinite long straight dislocation line is represented through a point, which is marked by the red nodes at different times. Then, only time integration is required to obtain the elastodynamics solution.
As discussed in Section 3, for screw dislocation, only G 33 works. The time derivative of G 33 is,
G ˙ 33 ( x x , t t ) = G 33 ( x x , t t ) t + G 33 ( x x , t t ) x η x η t = G 33 , t ( x x , t t ) G 33 , η ( x x , t t ) V η = t ˜ H ( t ˜ 2 l 2 / c T 2 ) 2 π μ ( t ˜ 2 l 2 / c T 2 ) 3 ( x η x η ) H ( t ˜ 2 l 2 / c T 2 ) 2 π c T 2 μ ( t ˜ 2 l 2 / c T 2 ) 3 V η
where t ˜ = t t , l 2 = ( x 1 x 1 ) 2 + ( x 2 x 2 ) 2 , η = 1 , 2 , G 33 | t = t = 0 , V k | t = 0 = 0 . Combining Equations (7) and (8), we can obtain the strain field function of non-uniformly moving screw dislocation,
β 31 E = 0 t [ μ G 33 , 2 ( x x , t t ) b 3 ξ 3 + ρ G ˙ 33 ( x x , t t ) V 2 b 3 ξ 3 ] d t = 0 t μ G 33 , 2 ( x x , t t ) b 3 ξ 3 μ G 33 , η ( x x , t t ) V η c T 2 V 2 b 3 ξ 3 d t 0 t ρ G 33 , t ( x x , t t ) V 2 b 3 ξ 3 d t = 0 t μ G 33 , 2 ( x x , t t ) b 3 ξ 3 μ G 33 , η ( x x , t t ) V η c T 2 V 2 b 3 ξ 3 d t + 0 t ρ G 33 ( x x , t t ) V 2 ˙ b 3 ξ 3 d t ρ G 33 ( x x , t t ) V 2 b 3 ξ 3 | t = 0 t = t = 0 t μ G 33 , 2 ( x x , t t ) b 3 ξ 3 μ G 33 , η ( x x , t t ) V η c T 2 V 2 b 3 ξ 3 d t + 0 t ρ G 33 ( x x , t t ) V 2 ˙ b 3 ξ 3 d t = 0 t μ G 33 , 2 ( x x , t t ) b 3 ξ 3 μ G 33 , 1 ( x x , t t ) V 1 c T 2 V 2 b 3 ξ 3 μ G 33 , 2 ( x x , t t ) V 2 c T 2 V 2 b 3 ξ 3 d t + 0 t ρ G 33 ( x x , t t ) V 2 ˙ b 3 ξ 3 d t = b 3 ξ 3 2 π c T 2 0 t V 2 ˙ t ˜ 2 l 2 / c T 2 + 1 V 2 2 c T 2 ( x 2 x 2 ) ( t ˜ 2 l 2 / c T 2 ) 3 V 1 V 2 c T 2 ( x 1 x 1 ) ( t ˜ 2 l 2 / c T 2 ) 3 H ( t ˜ 2 l 2 / c T 2 ) d t
β 32 E = 0 t [ μ G 33 , 1 ( x x , t t ) b 3 ξ 3 + ρ G ˙ 33 ( x x , t t ) V 1 b 3 ξ 3 ] d t = b 3 ξ 3 2 π c T 2 0 t V 1 ˙ t ˜ 2 l 2 / c T 2 + 1 V 1 2 c T 2 ( x 1 x 1 ) ( t ˜ 2 l 2 / c T 2 ) 3 V 1 V 2 c T 2 ( x 2 x 2 ) ( t ˜ 2 l 2 / c T 2 ) 3 H ( t ˜ 2 l 2 / c T 2 ) d t
β 11 E = β 12 E = β 13 E = β 21 E = β 22 E = β 23 E = β 33 E = 0
For screw dislocation considered with constant speed V, at time t = 0.15 ns, the numerical results of Equations (9) and (10) can be seen in Figure 2. The 2D elastodynamics stress fields derived from the degradation of 3D are consistent with 3D-discrete dislocation elastodynamics solutions in [15].
For edge dislocation, we can obtain,
G ˙ α β ( x x , t t ) = G α β , t ( x x , t t ) G α β , γ ( x x , t t ) V γ
where α , β , γ = 1 , 2 , G α β | t = t = 0 , V γ | t = 0 = 0 .
According to Equation (4), the derivative of t 2 l 2 / c 2 is given by,
t 2 l 2 / c 2 , γ = l c 2 t 2 l 2 / c 2 l x γ H ( t 2 l 2 / c 2 ) = x γ c 2 t 2 l 2 / c 2 H ( t 2 l 2 / c 2 )
The derivative of 1 l m t 2 l 2 / c 2 is,
1 l m t 2 l 2 / c 2 , γ = m x γ H ( t 2 l 2 / c 2 ) l m + 2 t 2 l 2 / c 2 + x γ H ( t 2 l 2 / c 2 ) c 2 l m ( t 2 l 2 / c 2 ) 3
where m 0 , c = c T , c L . Using Equations (13) and (14) and x α / x γ = δ α γ , the derivative of the 2D Green function is expressed as
G α β , γ ( x , t ) = 1 2 π ρ δ α γ x β + δ β γ x α l 4 4 x α x β x γ l 6 t 2 t 2 l 2 / c L 2 + t 2 l 2 / c L 2 H ( t 2 l 2 / c L 2 ) t 2 t 2 l 2 / c T 2 + t 2 l 2 / c T 2 H ( t 2 l 2 / c T 2 ) + x α x β x γ l 4 t 2 / c L 2 H ( t 2 l 2 / c L 2 ) ( t 2 l 2 / c L 2 ) 3 H ( t 2 l 2 / c L 2 ) c L 2 t 2 l 2 / c L 2 t 2 / c T 2 H ( t 2 l 2 / c T 2 ) ( t 2 l 2 / c T 2 ) 3 + H ( t 2 l 2 / c T 2 ) c T 2 t 2 l 2 / c T 2 + 2 δ α β x γ l 4 t 2 l 2 / c L 2 H ( t 2 l 2 / c L 2 ) t 2 H ( t 2 l 2 / c T 2 ) t 2 l 2 / c T 2 + δ α β x γ l 2 H ( t 2 l 2 / c L 2 ) c L 2 t 2 l 2 / c L 2 + t 2 H ( t 2 l 2 / c T 2 ) c T 2 ( t 2 l 2 / c T 2 ) 3
Substituting Equations (12) and (15) into Equation (7), we can obtain the elastodynamics field of non-uniformly moving edge dislocation,
β α ι E = 0 t ϵ n ς ι G α β , γ ( x x , t t ) C β γ κ ς ( t ) b κ ξ n ϵ ι ς n ρ G ˙ α β ( x x , t t ) V ς b β ξ n d t = 0 t ϵ n ς ι G α β , γ ( x x , t t ) C β γ κ ς ( t ) b κ ξ n + ρ G ˙ α β ( x x , t t ) V ς b β ξ n d t = 0 t ϵ n ς ι G α β , γ ( x x , t t ) ( C β γ κ ς ρ V γ V ς δ β κ ) b κ ξ n ρ G α β , t ( x x , t t ) V ς b β ξ n d t = 0 t ϵ n ς ι G α β , γ ( x x , t t ) ( C β γ κ ς ρ V γ V ς δ β κ ) b κ ξ n + ρ G α β ( x x , t t ) V ˙ ς b β ξ n d t ϵ n ς ι ρ G α β ( x x , t t ) V ς b β ξ n | t = 0 t = t = 0 t ϵ n ς ι G α β , γ ( x x , t t ) ( C β γ κ ς ρ V γ V ς δ β κ ) b κ ξ n + ρ G α β ( x x , t t ) V ˙ ς b β ξ n d t = 1 2 π ρ 0 t δ α γ x β ˜ + δ β γ x α ˜ l 4 4 x α ˜ x β ˜ x γ ˜ l 6 t ˜ 2 t ˜ 2 l 2 / c T 2 + t ˜ 2 l 2 / c T 2 H ( t ˜ 2 l 2 / c T 2 ) t ˜ 2 t ˜ 2 l 2 / c L 2 + t ˜ 2 l 2 / c L 2 H ( t ˜ 2 l 2 / c L 2 ) + x α ˜ x β ˜ x γ ˜ l 4 t ˜ 2 / c L 2 H ( t ˜ 2 l 2 / c L 2 ) ( t ˜ 2 l 2 / c L 2 ) 3 H ( t ˜ 2 l 2 / c L 2 ) c L 2 t ˜ 2 l 2 / c L 2 t ˜ 2 / c T 2 H ( t ˜ 2 l 2 / c T 2 ) ( t ˜ 2 l 2 / c T 2 ) 3 + H ( t ˜ 2 l 2 / c T 2 ) c T 2 t ˜ 2 l 2 / c L 2 + 2 δ α β x γ ˜ l 4 t ˜ 2 l 2 / c L 2 H ( t ˜ 2 l 2 / c L 2 ) t ˜ 2 H ( t ˜ 2 l 2 / c T 2 ) t ˜ 2 l 2 / c T 2 + δ α β x γ ˜ l 2 H ( t ˜ 2 l 2 / c L 2 ) c L 2 t ˜ 2 l 2 / c L 2 + t ˜ 2 H ( t ˜ 2 l 2 / c T 2 ) c T 2 ( t 2 l 2 / c T 2 ) 3 ϵ n ς ι ( C β γ κ ς ρ V γ V ς δ β κ ) b κ ξ n d t + ρ δ α β l 2 t ˜ 2 H ( t ˜ 2 l 2 / c T 2 ) t ˜ 2 l 2 / c T 2 t ˜ 2 l 2 / c L 2 H ( t ˜ 2 l 2 / c L 2 ) + x α ˜ x β ˜ l 4 ( 2 t ˜ 2 + l 2 / c T 2 ) H ( t ˜ 2 l 2 / c T 2 ) t ˜ 2 l 2 / c T 2 + ( 2 t ˜ 2 l 2 / c L 2 ) H ( t ˜ 2 l 2 / c L 2 ) t ˜ 2 l 2 / c L 2 ϵ n ς ι V ˙ ς b β ξ n d t
where t ˜ = t t , x α ˜ = x α x α ( α = 1 , 2 ), l 2 = x 1 ˜ 2 + x 2 ˜ 2 . The 2D elastodynamics results of Equations (9), (10), and (16) are in agreement with the previous works [31].

4.2. Spatial–Temporal Decoupling Strategy at 3D

However, the spatial–temporal decoupling strategies at 2D dislocation discussed above are not applied to arbitrarily moving 3D dislocation, because their 3D geometry and curvature keep changing during plastic deformation. Therefore, the alternative strategy is to carry out the temporal integration firstly. As a result of the existence of the Dirac delta function δ ( t t x x c ) in Green’s function, the retarded solution used in Equation (1) is obtained as follows by performing temporal integration.
L ( t R / c ) = x such that x x ( t ret ) = c ( t t ret )
where x is a source position on the retarded dislocation and L ( t R / c ) represents the collection of retarded points x on an arbitrarily moving 3D dislocation. x is the observer point. t ret is retarded time corresponding to a specific elastic wave velocity c. For retarded dislocation, we make use of a simple example of subsonic 3D dislocation, sliding in plane O x 1 x 3 , with constant velocity V to elaborate the case. For retarded dislocation, we make use of a simple example of subsonic 3D dislocation, sliding in plane O x 1 x 3 with constant speed V. As shown in Figure 1c, for the arbitrary observer point P ( x 1 P , x 2 P , 0 ) in plane O x 1 x 2 , after the temporal integration, the points ( V t r , 0 , x 3 r ) on the retarded dislocation line satisfy ( x 1 P V t r ) 2 + ( x 2 P ) 2 + ( x 3 B ) 2 = c T 2 ( t t r ) 2 , where t r is the retarded time corresponding to the shear wave ( c T ). More examples of retarded dislocations can be seen in [15]. For 3D moving dislocation with arbitrary curve configuration, complex curve integrals along the retarded dislocation contain the 3D space information of elastodynamics fields.
This technique of using retarded potentials simplifies the calculation of coupled spatial and temporal integration to a spatial integral over retarded dislocation positions. This retarded potential method is standard in electrodynamics [42]. This method extends the technique of retarded potentials, which was originally used to describe the electrodynamics of charged particles moving near the speed of light. The retarded functions are the consequence of the finite speed of the wave and with respect to the effect of retardation. This idea is applicable to any dynamical theory with retardation, including elastodynamics [39]. The simplest example is the well-known Doppler effect in acoustics.
Accordingly, the elastodynamics stress field of arbitrarily moving 3D dislocation is as follows [15],
σ ( x , t ) = μ 4 π s ( x , t ) + s T ( x , t ) ,
where the auxiliary tensor s ( x , t ) is given by,
s ( x , t ) = L ( t R / c T ) 1 R 2 R R · V / c T 2 R ( b × ) + ν 1 2 ν b · ( R × ) I 3 b 4 R · b R 2 R R × d L t L ( t R / c T ) 1 c T R R R · V / c T ν 1 2 ν b · ( R × ) I + b 2 R · b R 2 R ( R × ) d L + L ( t R / c L ) 1 2 ν ( 1 ν ) R 2 R R · V / c L 1 3 ν 1 2 ν R ( b × ) + ν 2 ( 1 2 ν ) 2 b · ( R × ) I + b 6 R · b R 2 R R × d L t L ( t R / c L ) 1 2 ν ( 1 ν ) c L R R R · V / c L ν 1 2 ν R ( b × ) ν 2 1 2 ν 2 b · ( R × ) I + R · b R 2 R ( R × ) d L + 1 / c L 1 / c T L ( t R κ ) 6 c T 2 κ R 2 R R · V κ R ( b × ) + b 5 R · b R 2 R R × d L d κ + t L ( t R / c T ) 1 c T 2 R R · V / c T ν 1 2 ν b R · b R 2 R · ( V × ) I + b R · b R 2 R ( V × ) d L + t L ( t R / c L ) 1 c L 2 R R · V / c L ν 1 2 ν R · b R 2 R · ( V × ) I + R · b R 2 R ( V × ) d L 1 / c L 1 / c T L ( t R κ ) 1 R R R · V κ ν 1 2 ν b 3 R · b R 2 R · ( V × ) I + b 3 R · b R 2 R ( V × ) d L d κ + L ( t R / c T ) 1 c T R R R · V / c T ν 1 2 ν b 3 R · b R 2 R · ( V × ) I + b 3 R · b R 2 R ( V × ) d L L ( t R / c L ) 1 c L R R R · V / c L ν 1 2 ν b 3 R · b R 2 R · ( V × ) I + b 3 R · b R 2 R ( V × ) d L
and s T ( x , t ) is its transpose. Here, R is the magnitude of R , which is the vector connecting the source point to the observer point, namely, R = x x .

4.3. Consistency

Different spatial and temporal integration strategies lead to different formulas of the elastodynamics solution (see Equations (16) and (19)). The case of an infinitely long straight dislocation line described in Figure 1a–c schematically shows the difference of the integration strategies. Figure 1b needs to perform temporal integration over each point at different times, while Figure 1c needs to perform spatial integration over the retarded dislocation lines. For 3D dislocation, carrying out the space integration along retarded dislocation is capable of completing calculation of the 3D elastodynamics fields. For 2D dislocation, it is assumed to be along x 3 all the time; the intersection point of the dislocation line and plane O x 1 x 2 can well record the position of dislocation. Performing the straight line integration along the straight dislocation is specific. Thereby, performing the temporal integration can calculate the 2D elastodynamics fields after space integration. It is not apparent to see their physical connections from these figures.
To simplify the discussion and clarify their consistency, the spacial and temporal integration of a basic term is derived as follows, using the integration strategies described in Section 4.1 and Section 4.2. If the spatial integration is firstly done,
t 0 t f ( x , t ) δ t t x c d x 3 d t = t 0 t 2 x f ( x , t ) H ( t 2 l 2 / c 2 ) ( t t ) 2 l 2 / c 2 | x = c ( t t ) d t
where x 3 ϵ [ c 2 ( t t 0 ) 2 l 2 , c 2 ( t t 0 ) 2 l 2 ] , l 2 = ( x 1 x 1 ) 2 + ( x 2 x 2 ) 2 .
If the temporal integration is firstly done,
t 0 t f ( x , t ) δ t t x c d t d x 3 = f ( x , t ) t = t x c d x 3
where x ϵ [ 0 , c ( t t 0 ) ] , therefore t = t x c ϵ [ t , t 0 ] , mathematically.
Furthermore, d t = d ( t x c ) = d x c d x 3 d x 3 = x 3 c x d x 3 . Consequently, Equation (21) becomes
t 0 t f ( x , t ) δ t t x c d t d x 3 = t 0 t 2 x f ( x , t ) H ( t 2 x 2 / c 2 ) ( t t ) 2 l 2 / c 2 | x = c ( t t ) ) d t
Equation (22), which firstly calculates the temporal integration, is exactly the same as Equation (20), which firstly carries out the spatial integration. Therefore, for more complex formulas associated with elastodynamics solutions, they are exactly the same.

5. Conclusions

The significant effects of elastodynamics on a dislocation stress field have been recognized for more than 70 years. However, the marked difference in elastodynamics stress expression for 3D and 2D dislocation is seldomly directly discussed. In the current work, the physical origin of their difference and connections is clearly discussed from the perspective of different spatial and temporal decoupling strategies.
For the 2D case, the dislocation configuration keeps straight along a specific direction, which makes it suitable to firstly carry out the spatial integration along its line direction. For the 3D case, the evolution of the dislocation configuration is very complicated, one cannot obtain generalized results if the spatial integration is firstly carried out. Therefore, the temporal integration is firstly performed. The consistency between these two cases is clearly demonstrated through the 2D and 3D elastodynamics Green tensor and a simple generalized formula. This work helps to understand the difference and connections between 2D and 3D dislocation elastodynamics, which leads to an enlightening opportunity for establishing general high-speed dislocation theory, developing a numerical calculation method based on dislocation elastodynamics, and understanding the collective behavior of high-speed dislocations under complicated shock-loading conditions.

Author Contributions

S.L.: Formal analysis, Verification, Writing—review and editing. Y.C.: Conceptualization, Supervision, Formal analysis, Writing—original draft, review, and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China under Grant No. 11972208, 11921002, National Key Laboratory of Shock Wave and Detonation Physics under Grant No. JCKYS2021212004 and Tsinghua University Initiative Scientific Research Program.

Data Availability Statement

Not applicable.

Acknowledgments

We greatly thank M. Lazar, Nasr Ghoniem, and Giacomo Po for the insightful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Calculating the space integral over 3D Green function in the x 3 direction and 2D elastodynamicss field, we used the following auxiliary formula, where x = x 1 2 + x 2 2 + x 3 2 , l = x 1 2 + x 2 2 , x i ( i = 1 , 2 , 3 ) , and x α ( α = 1 , 2 ) are the Cartesian coordinate components, and c T and c L are shear and longitudinal wave velocities of the material.
δ t x c x d x 3 = 1 x 1 c d x d x 3 | x = c t = 2 x 1 c x 3 x | x = c t = 2 t 2 l 2 / c 2 δ t x c x 3 d x 3 = 2 c 2 t 2 t 2 l 2 / c 2 x α x 3 δ t x c x 3 d x 3 = 2 x α c t 2 x 3 2 δ t x c x 3 d x 3 = 2 t 2 l 2 / c 2 t 2 x 3 2 δ t x c x 5 d x 3 = 2 t 2 l 2 / c 2 c 2 t 4 x 3 3 δ t x c x 5 d x 3 = 2 ( t 2 l 2 / c 2 ) c t 4 x 3 δ t x c x 2 d x 3 = 2 t
1 x 3 x / c L x / c T τ δ t τ d τ d x 3 = 1 x 1 / c L 1 / c T κ δ t x κ d x 3 d κ = 2 l 2 t 2 l 2 / c T 2 + 2 l 2 t 2 l 2 / c L 2 1 x 5 x / c L x / c T τ δ t τ d τ d x 3 = 4 3 l 4 t 2 t 2 l 2 c T 2 3 / 2 + 4 3 l 4 t 2 t 2 l 2 c L 2 3 / 2 2 l 2 c T 2 t 2 t 2 l 2 c T 2 + 2 l 2 c L 2 t 2 t 2 l 2 c L 2 1 x 7 x / c L x / c T τ δ t τ d τ d x 3 = 2 t 2 l 2 / c T 2 15 l 4 8 + 4 l 2 c T 2 t 2 + 3 l 4 c T 4 t 4 + 2 t 2 l 2 / c L 2 15 l 4 8 + 4 l 2 c L 2 t 2 + 3 l 4 c L 4 t 4
x α x 3 x 5 x / c L x / c T τ δ t τ d τ d x 3 = 2 x α 3 t 2 1 c T 3 1 c L 3 x 3 2 x 5 x / c L x / c T τ δ t τ d τ d x 3 = 2 3 l 2 t 2 ( t 2 l 2 / c T 2 ) 3 + 2 3 l 2 t 2 ( t 2 l 2 / c L 2 ) 3 x 3 2 x 7 x / c L x / c T τ δ t τ d τ d x 3 = 2 5 1 l 2 c T 2 t 4 + 2 3 l 4 t 2 t 2 l 2 c T 2 3 / 2 + 2 5 1 l 2 c L 2 t 4 + 2 3 l 4 t 2 t 2 l 2 c L 2 3 / 2 x 3 3 x 7 x / c L x / c T τ δ t τ d τ d x 3 = t 2 l 2 / c T 2 15 2 l 4 + 1 c T 2 l 2 t 2 3 c T 4 t 4 + c L 2 t 2 l 2 15 c L 5 l 4 t 4 2 c L 4 t 4 + c L 2 l 2 t 2 3 l 4
δ t x c x 2 d x 3 , t = 4 c 2 t + 2 l 2 t 2 ( c 2 t 2 l 2 ) 3 δ t x c x 4 d x 3 , t = 8 c 2 t 2 + 6 l 2 c 2 t 4 ( c 2 t 2 l 2 ) 3 1 x 3 x / c L x / c T τ δ t τ d τ d x 3 , t = 2 t l 2 t 2 l 2 / c T 2 + 2 t l 2 t 2 l 2 / c L 2 x 3 2 δ t x c x 4 d x 3 , t = 2 t 2 l 2 / c 2 c t 3 , t = 2 c t 2 t 2 l 2 / c 2 6 t 2 l 2 / c 2 c t 4 x 3 2 x 5 x / c L x / c T τ δ t τ d τ d x 3 , t = 4 3 l 2 t 3 ( t 2 l 2 / c T 2 ) 3 2 l 2 t t 2 l 2 / c T 2 4 3 l 2 t 3 ( t 2 l 2 / c L 2 ) 3 + 2 l 2 t t 2 l 2 / c L 2

References

  1. Gray, G.T., III. High-strain-rate deformation: Mechanical behavior and deformation substructures induced. Annu. Rev. Mater. Res. 2012, 42, 285–303. [Google Scholar] [CrossRef]
  2. Yaghoobi, M.; Voyiadjis, G.Z.; Sundararaghavan, V. Crystal Plasticity Simulation of Magnesium and Its Alloys: A Review of Recent Advances. Crystals 2021, 11, 435. [Google Scholar] [CrossRef]
  3. Ma, M.; Ding, H.; Huang, Y.; Tian, C.W.; Langdon, T.G. Microstructural and Hardness Evolution in a Duplex Stainless Steel Processed by High-Pressure Torsion. Crystals 2020, 10, 1138. [Google Scholar] [CrossRef]
  4. Javanbakht, M. High pressure phase evolution under hydrostatic pressure in a single imperfect crystal due to nanovoids. Materialia 2021, 20, 101199. [Google Scholar] [CrossRef]
  5. Roy, A.M. Barrierless melt nucleation at solid-solid interface in energetic nitramine octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine. Materialia 2021, 15, 101000. [Google Scholar] [CrossRef]
  6. Meyers, M.A. Dynamic Behavior of Materials; John Wiley & Sons: Hoboken, NJ, USA, 1994. [Google Scholar]
  7. Bisht, A.; Neogi, A.; Mitra, N.; Jagadeesh, G.; Suwas, S. Investigation of the elastically shock-compressed region and elastic–plastic shock transition in single-crystalline copper to understand the dislocation nucleation mechanism under shock compression. Shock Waves 2019, 29, 913–927. [Google Scholar] [CrossRef]
  8. Sichani, M.M.; Spearot, D.E. A molecular dynamics study of dislocation density generation and plastic relaxation during shock of single crystal Cu. J. Appl. Phys. 2016, 120, 045902. [Google Scholar] [CrossRef]
  9. Tanguy, D.; Mareschal, M.; Lomdahl, P.S.; Germann, T.C.; Holian, B.L.; Ravelo, R. Dislocation nucleation induced by a shock wave in a perfect crystal: Molecular dynamics simulations and elastic calculations. Phys. Rev. B 2003, 68, 165–168. [Google Scholar] [CrossRef]
  10. Zepeda-Ruiz, L.A.; Stukowski, A.; Oppelstrup, T.; Bulatov, V.V. Probing the limits of metal plasticity with molecular dynamics simulations. Nature 2017, 550, 492–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Eshelby, J. Uniformly Moving Dislocations. Proc. Phys. Soc. 1949, 62, 307. [Google Scholar] [CrossRef]
  12. Frank, F.C. On the Equations of Motion of Crystal Dislocations. Proc. Phys. Soc. Sect. A 1949, 62, 131–134. [Google Scholar] [CrossRef]
  13. Markenscoff, X.; Clifton, R. The nonuniformly moving edge dislocation. J. Mech. Phys. Solids 1981, 29, 253–262. [Google Scholar] [CrossRef]
  14. Nosenko, V.; Zhdanov, S.; Morfill, G. Supersonic dislocations observed in a plasma crystal. Phys. Rev. Lett. 2007, 99, 25002. [Google Scholar] [CrossRef] [Green Version]
  15. Cui, Y.; Po, G.; Pellegrini, Y.P.; Lazar, M.; Ghoniem, N. Computational 3-dimensional dislocation elastodynamics. J. Mech. Phys. Solids 2019, 126, 20–51. [Google Scholar] [CrossRef]
  16. Gurrutxaga-Lerma, B.; Balint, D.S.; Dini, D.; Eakins, D.E.; Sutton, A.P. A dynamic discrete dislocation plasticity method for the simulation of plastic relaxation under shock loading. R. Soc. Math. Phys. Eng. Sci. 2013, 469, 20130141. [Google Scholar] [CrossRef]
  17. Gurrutxaga-Lerma, B.; Balint, D.S.; Dini, D.; Sutton, A.P. Elastodynamic image forces on dislocations. R. Soc. Math. Phys. Eng. Sci. 2015, 471, 20150433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Blaschke, D.N. A general solution for accelerating screw dislocations in arbitrary slip systems with reflection symmetry. J. Mech. Phys. Solids 2021, 152, 104448. [Google Scholar] [CrossRef]
  19. Cui, Y.; Wang, T.; Luo, S.; Li, Z.; Li, Z. A discrete–continuous model of three-dimensional dislocation elastodynamics. Int. J. Plast. 2022, 152, 103221. [Google Scholar] [CrossRef]
  20. Pellegrini, Y.P.; Lazar, M. On the gradient of the Green tensor in two-dimensional elastodynamic problems, and related integrals: Distributional approach and regularization, with application to nonuniformly moving sources. Wave Motion 2015, 57, 44–63. [Google Scholar] [CrossRef] [Green Version]
  21. Bringa, E.M.; Rosolankova, K.; Rudd, R.E.; Remington, B.A.; Wark, J.S.; Duchaineau, M.; Kalantar, D.H.; Hawreliak, J.; Belak, J. Shock deformation of face-centred-cubic metals on subnanosecond timescales. Nat. Mater. 2006, 5, 805–809. [Google Scholar] [CrossRef]
  22. Chu, H.J.; Wang, J.; Beyerlein, I.J. Anomalous reactions of a supersonic coplanar dislocation dipole: Bypass or twinning? Scr. Mater. 2012, 67, 69–72. [Google Scholar] [CrossRef]
  23. Gumbsch, P.; Gao, H. Dislocations faster than the speed of sound. Science 1999, 283, 965–968. [Google Scholar] [CrossRef] [PubMed]
  24. Li, Q.J.; Li, J.; Shan, Z.W.; Ma, E. Strongly correlated breeding of high-speed dislocations. Acta Mater. 2016, 119, 229–241. [Google Scholar] [CrossRef] [Green Version]
  25. Kim, S.; Kim, H.; Kang, K.; Kim, S.Y. Relativistic effect inducing drag on fast-moving dislocation in discrete system. Int. J. Plast. 2019, 126, 102629. [Google Scholar] [CrossRef]
  26. Cottrell, A.H. Theory of dislocations. Prog. Met. Phys. 1968, 1, 77–96. [Google Scholar] [CrossRef]
  27. Jin, Z.; Gao, H.; Gumbsch, P. Energy radiation and limiting speeds of fast moving edge dislocations in tungsten. Phys. Rev. B 2008, 77, 094303. [Google Scholar] [CrossRef] [Green Version]
  28. Tsuzuki, H.; Branicio, P.S.; Rino, J.P. Molecular dynamics simulation of fast dislocations in copper. Acta Mater. 2009, 57, 1843–1855. [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. Lazar, M. The gauge theory of dislocations: A uniformly moving screw dislocation. R. Soc. Lond. A Math. Phys. Eng. Sci. 2009, 465, 2505–2520. [Google Scholar] [CrossRef]
  31. Lazar, M.; Pellegrini, Y.P. Distributional and regularized radiation fields of non-uniformly moving straight dislocations, and elastodynamic Tamm problem. J. Mech. Phys. Solids 2016, 96, 632–659. [Google Scholar] [CrossRef] [Green Version]
  32. Chen, H.; Levitas, V.I.; Xiong, L.; Zhang, X. Stationary Dislocation Motion at Stresses Significantly below the Peierls Stress: Example of Shuffle Screw and 60 Dislocations in Silicon. Acta Mater. 2021, 206, 116623. [Google Scholar] [CrossRef]
  33. Li, Q.J.; Li, J.; Shan, Z.W.; Ma, E. Surface Rebound of Relativistic Dislocations Directly and Efficiently Initiates Deformation Twinning. Phys. Rev. Lett. 2018, 117, 165501. [Google Scholar] [CrossRef] [Green Version]
  34. Markenscoff, X. Evolution equation of moving defects: Dislocations and inclusions. Int. J. Fract. 2010, 166, 35–40. [Google Scholar] [CrossRef] [Green Version]
  35. Lazar, M. The fundamentals of non-singular dislocations in the theory of gradient elasticity: Dislocation loops and straight dislocations. Int. J. Solids Struct. 2013, 50, 352–362. [Google Scholar] [CrossRef] [Green Version]
  36. Lazar, M.; Kirchner, H.O.K. Generalised plane strain embedded in three-dimensional anisotropic elasticity. Philos. Mag. 2021, 101, 2584–2598. [Google Scholar] [CrossRef]
  37. Mura, T. Continuous distribution of moving dislocations. Philos. Mag. 1963, 8, 843–857. [Google Scholar] [CrossRef]
  38. Callias, C.; Markenscoff, X.; Ni, L.Q. A singular asymptotic expansion for the field near a moving dislocation loop. Q. Appl. Math. 1990, 48, 113–132. [Google Scholar] [CrossRef] [Green Version]
  39. Lazar, M. On the non-uniform motion of dislocations: The retarded elastic fields, the retarded dislocation tensor potentials and the Liénard—Wiechert tensor potentials. Philos. Mag. 2013, 93, 749–776. [Google Scholar] [CrossRef] [Green Version]
  40. Eringen, A.C.; Suhubi, E.S.; Bland, D.R. Elastodynamics, vols. 1 and 2. Phys. Today 1977, 30, 65. [Google Scholar] [CrossRef]
  41. Achenbach, J.D. Reciprocity in Elastodynamics; Cambridge University Press: Cambridge, MA, USA, 2003. [Google Scholar]
  42. Jackson, J.D. Classical Electrodynamics; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
Figure 1. (a) Schematic showing a high-speed infinite-long straight dislocation line; (b) Schematic showing the integration strategy of a 2D dislocation line; (c) Schematic showing the integration strategy based on retarded dislocation.
Figure 1. (a) Schematic showing a high-speed infinite-long straight dislocation line; (b) Schematic showing the integration strategy of a 2D dislocation line; (c) Schematic showing the integration strategy based on retarded dislocation.
Crystals 12 00363 g001
Figure 2. Stress field of screw dislocation with V = 0 or 0.77 c T , at t = 0.15 ns in the O x 1 x 2 plane ( μ is the shear model of the material). (ad) Two-dimensional (2D) elastodynamics solution by Equations (9) and (10), (eh) 3D elastodynamics field in [15].
Figure 2. Stress field of screw dislocation with V = 0 or 0.77 c T , at t = 0.15 ns in the O x 1 x 2 plane ( μ is the shear model of the material). (ad) Two-dimensional (2D) elastodynamics solution by Equations (9) and (10), (eh) 3D elastodynamics field in [15].
Crystals 12 00363 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, S.; Cui, Y. Elastodynamics Field of Non-Uniformly Moving Dislocation: From 3D to 2D. Crystals 2022, 12, 363. https://doi.org/10.3390/cryst12030363

AMA Style

Luo S, Cui Y. Elastodynamics Field of Non-Uniformly Moving Dislocation: From 3D to 2D. Crystals. 2022; 12(3):363. https://doi.org/10.3390/cryst12030363

Chicago/Turabian Style

Luo, Shichao, and Yinan Cui. 2022. "Elastodynamics Field of Non-Uniformly Moving Dislocation: From 3D to 2D" Crystals 12, no. 3: 363. https://doi.org/10.3390/cryst12030363

APA Style

Luo, S., & Cui, Y. (2022). Elastodynamics Field of Non-Uniformly Moving Dislocation: From 3D to 2D. Crystals, 12(3), 363. https://doi.org/10.3390/cryst12030363

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