A computationally efficient numerical integration scheme for non-linear planestress/strain FEM applications using one-point constitutive model evaluation

This work presents a proposal for employing reduced numerical integration in the formulation of the 4-node quadrilateral solid finite element. The use of these low-order integration rules leads to numerical instabilities such as those producing the hourglass effect. The proposed procedure allows evaluating a given constitutive model only in one integration point, achieving an attractive computational cost reduction and, also, successfully controls the hourglass effect. A validation of the proposal is included and discussed throughout the paper. To show the efficiency of the proposal, several application examples of masonry structures are studied and discussed. To represent the non-linear mechanical behaviour of masonry a plastic-damage model is implemented within the application of this sub-integration scheme. Also, in order to have a full and computationally efficient strategy to determine the behaviour of masonry structures, involving its evolution to collapse, a homogenization technique with a macro-modeling approach is used. The methodology discussed throughout this paper demonstrates a substantial computational cost reduction and an improved approximation of the non-linear problem evidenced by a reduction of up to 85% of the computational time for some cases.