Comparative analysis of finite element formulations at plane loading of an elastic body

Abstract

The aim of the work - comparison of the results of determining the parameters of the stress-strain state of plane-loaded elastic bodies based on the finite element method in the formulation of the displacement method and in the mixed formulation. Methods. Algorithms of the finite element method in various formulations have been developed and applied. Results. In the Cartesian coordinate system, to determine the stress-strain state of an elastic body under plane loading, a finite element of a quadrangular shape is used in two formulations: in the formulation of the method of displacements with nodal unknowns in the form of displacements and their derivatives, and in a mixed formulation with nodal unknowns in the form of displacements and stresses. The approximation of displacements through the nodal unknowns when obtaining the stiffness matrix of the finite element was carried out using the form function, whose elements were adopted Hermite polynomials of the third degree. Upon receipt of the deformation matrix, the displacements and stresses of the internal points of the finite element were approximated through nodal unknowns using bilinear functions. The stiffness matrix of the quadrangular finite element in the formulation of the displacement method is obtained on the basis of a functional based on the difference between the actual workings of external and internal forces under loading of a solid. The matrix of deformation of the finite element was formed on the basis of a mixed functional obtained from the proposed functional by repla-cing the actual work of internal forces with the difference between the total and additional work of internal forces when loading the body. The calculation example shows a significant advantage of using a finite element in a mixed formulation.