Finite element approximation of field dislocation mechanics

A tool for studying links between continuum plasticity and dislocation theory within a field framework is presented. A finite element implementation of the geometrically linear version of a recently proposed theory of field dislocation mechanics (J. Mech. Phys. Solids 49 (2001) 761; Proc. Roy. Soc. 459 (2003) 1343; J. Mech. Phys. Solids 52 (2004) 301) represents the main idea behind the tool. The constitutive ingredients of the theory under consideration are simply elasticity and a specification of dislocation velocity and nucleation. The set of equations to be approximated are non-standard in the context of solid mechanics applications. It comprises the standard second-order equilibrium equations, a first-order div-curl system for the elastic incompatibility, and a first-order, wave-propagative system for the evolution of dislocation density. The latter two sets of equations require special treatment as the standard Galerkin method is not adequate, and are solved utilizing a least-squares finite element strategy. The implementation is validated against analytical results of the classical elastic theory of dislocations and analytical results of the theory itself. Elastic stress fields of dislocation distributions in generally anisotropic media of finite extent, deviation from elastic response, yield-drop, and back-stress are shown to be natural consequences of the model. The development of inhomogeneity, from homogeneous initial conditions and boundary conditions corresponding to homogeneous deformation in conventional plasticity, is also demonstrated. To our knowledge, this work represents the first computational implementation of a theory of dislocation mechanics where no analytical results, singular solutions in particular, are required to formulate the implementation. In particular, a part of the work is the first finite element implementation of Kröner's linear elastic theory of continuously distributed dislocations in its full generality.