Recovery of strong equilibrium from displacement‐based finite element models of Reissner‐Mindlin plates

In this paper, we extend to Reissner‐Mindlin plate bending problems a technique, originally proposed in the context of two‐dimensional and three‐dimensional continua, for recovering fully equilibrated stresses from the solution of a compatible finite element model. The technique involves partition of unity functions and the analyses of overlapping star patches modelled with hybrid equilibrium plate elements. The patches are subjected to balanced systems of loads composed of partitioned and fictitious loads, where the latter are derived from the stresses of the compatible solution. The special case of assumed linear displacement fields of both deflection and rotation for the compatible model is included. This case requires additional fields of stress resultants to correct possible rotational imbalances of star patches, and these are derived elementwise. Other cases of nonconforming elements are briefly considered. Numerical examples are presented to illustrate the effectiveness of these techniques in terms of the deviation of the recovery, which compares the complementary strain energy of a recovered solution with that obtained by a global equilibrated analysis based on the same stress approximations.