Fast-sum method for the elastic field of three-dimensional dislocation ensembles

The elastic field of complex shape ensembles of dislocation loops is developed as an essential ingredient in the dislocation dynamics method for computer simulation of mesoscopic plastic deformation. Dislocation ensembles are sorted into individual loops, which are then divided into segments represented as parametrized space curves. Numerical solutions are presented as fast numerical sums for relevant elastic field variables (i.e., displacement, strain, stress, force, self-energy, and interaction energy). Gaussian numerical quadratures are utilized to solve for field equations of linear elasticity in an infinite isotropic elastic medium. The accuracy of the method is verified by comparison of numerical results to analytical solutions for typical prismatic and slip dislocation loops. The method is shown to be highly accurate, computationally efficient, and numerically convergent as the number of segments and quadrature points are increased on each loop. Several examples of method applications to calculations of the elastic field of simple and complex loop geometries are given in infinite crystals. The effect of crystal surfaces on the redistribution of the elastic field is demonstrated by superposition of a finite-element image force field on the computed results. copyright 1999 The American Physical Society