Strategies for the Computation of Configurational Forces in Dissipative Media

Configurational forces can be interpreted as driving forces on material inhomogeneities such as crack tips. In dissipative media the total configurational force on an inhomogeneity consists of an elastic contribution and a contribution due to the dissipative processes in the material. For the computation of discrete configurational forces acting at the nodes of a finite element mesh, the elastic and dissipative contributions must be evaluated at integration point level. While the evaluation of the elastic contribution is straightforward, the evaluation of the dissipative part is faced with certain difficulties. This is because gradients of internal variables are necessary in order to compute the dissipative part of the configurational force. For the sake of efficiency, these internal variables are usually treated as local history data at integration point level in finite element (FE) implementations. Thus, the history data needs to be projected to the nodes of the FE mesh in order to compute the gradients by means of shape function interpolations of nodal data as it is standard practice. However, this is a rather cumbersome method which does not easily integrate into standard finite element frameworks. An alternative approach which facilitates the computation of gradients of local history data is investigated in this work. This approach is based on the definition of subelements within the elements of the FE mesh and allows for a straightforward integration of the configurational force computation into standard finite element software. The suitability and the numerical accuracy of different projection approaches and the subelement technique are discussed and analyzed exemplarily within the context of a crystal plasticity model. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)