Mixed finite element formulations of strain-gradient elasticity problems Original Research Article

Theories on intrinsic or material length scales find applications in the modeling of size-dependent phenomena. In elasticity, length scales enter the constitutive equations through the elastic strain energy function, which, in this case, depends not only on the strain tensor but also on gradients of the rotation and strain tensors. In the present paper, the strain-gradient elasticity theories developed by Mindlin and co-workers in the 1960s are treated in detail. In such theories, when the problem is formulated in terms of displacements, the governing partial differential equation is of fourth order. If traditional finite elements are used for the numerical solution of such problems, then C1 displacement continuity is required. An alternative “mixed” finite element formulation is developed, in which both the displacement and the displacement gradients are used as independent unknowns and their relationship is enforced in an“integral-sense”. A variational formulation is developed which can be used for both linear and non-linear strain-gradient elasticity theories. The resulting finite elements require only C0 continuity and are simple to formulate. The proposed technique is applied to a number of problems and comparisons with available exact solutions are made.