Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

• Published in: Modelling and simulation in materials science and engineering Vol. 28; no. 7; pp. 75001 - 75034
• Main Authors: Queyreau, Sylvain, Hoang, Khiem, Shi, Xiangjun, Aubry, Sylvie, Arsenlis, Athanasios
• Format: Journal Article
• Language: English
• Published: United States IOP Publishing 01.10.2020
• Subjects: analytical model; dislocation dynamics simulations; dislocation induced tractions; dislocation interface interactions; ENGINEERING; hybrid dislocation dynamics finite elements methods simulations; mathematics and computing
• ISSN: 0965-0393; 1361-651X; 1361-651X
• DOI: 10.1088/1361-651X/aba736

An analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561-587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.