Nonlinear finite element calculations of residual stresses in dislocated crystals

A finite element algorithm for the determination of both residual stresses and lattice orientation in dislocated crystals is presented. The applied theoretical approach is based on the continuum theory of dislocations. Due to the compatibility condition, an initial plastic distortion field has been used to introduce an elastic incompatibility corresponding to the assumed dislocation density tensor field. The presented method gives the possibility to develop a numerical code suitable for determination of the stress field and lattice orientation. Two numerical examples are given: a low-angle dislocation wall and a subboundary structure. The method can be applied especially to dislocation sets composed of a few groups of monomial dislocations.