A variational r-adaption and shape-optimization method for finite-deformation elasticity

This paper is concerned with the formulation of a variational r-adaption method for finite-deformation elastostatic problems. The distinguishing characteristic of the method is that the variational principle simultaneously supplies the solution, the optimal mesh and, in problems of shape optimization, the equilibrium shapes of the system. This is accomplished by minimizing the energy functional with respect to the nodal field values as well as with respect to the triangulation of the domain of analysis. Energy minimization with respect to the referential nodal positions has the effect of equilibrating the energetic or configurational forces acting on the nodes. We derive general expressions for the configurational forces for isoparametric elements and non-linear, possibly anisotropic, materials under general loading. We illustrate the versatility and convergence characteristics of the method by way of selected numerical tests and applications, including the problem of a semi-infinite crack in linear and non-linear elastic bodies; and the optimization of the shape of elastic inclusions