A mixed finite element for laminated composite plates based on the use of bubble functions

A new laminated composite plate finite element is proposed that is numerically stable and accurate in displacements and stresses, including transverse shear stress. The formulation is based on the Hellinger—Reissner principle with Mindlin kinematics. All stress components are given independent approximations and do not satisfy equilibrium conditions a priori. A novel feature of the formulation is the additive decomposition of the displacement field into two parts corresponding to nodal interpolations and independent local basis functions. The additional basis functions and their associated parameters play an important role in characterizing the accuracy of the element. These functions eliminate shear locking in the lower order elements and provide additional variational constraints on the stresses, leading to very accurate results. A 4‐node and a 9‐node version are developed and it is shown that both elements pass the patch test suggested by Zienkiewicz et al. and are stable in the sense of the Babuška—Brezzi condition. The special structure of the element flexibility matrix provides computational efficiency approaching that of displacement based formulations.