Extended penalty coefficients for elimination the locking effects in moderately thick beam and plate finite elements

The present paper is dedicated to presentation and energy verification of the methods of stabilization the strain energy by penalty coefficients. Verification of the methods is based on the consistency and ellipticity conditions to be satisfied by the finite elements. Three methods of stabilization are discussed. The first does not satisfy the above requirements. The second is consistent but cannot eliminate parasitic energy terms. The third method, proposed by the author, is based on the decomposition of the element stiffness matrix. The method can help to eliminate locking of the finite elements. For two-noded beam element with linear shape functions and exact integration a stabilized free of locking (and elliptical) element is received (equivalent to reduced integration element). Two plate finite elements are analyzed: four-noded rectangular element and DSG triangle. A new method of stabilization with the use of four independent parameters is proposed. The finite elements with this kind of stabilization satisfy the consistency condition. In the rectangular element it was not possible to eliminate one parasitic term of energy which appears during the procedure. For DSG triangle all parasitic terms of energy are eliminated. The penalty coefficients depends on the geometry of the triangle