Xiang-Rong Fu

Analytical trial function method for development of new 8‐node plane element based on the variational principle containing Airy stress function

Purpose − The purpose of this paper is to propose a novel and simple strategy for construction of hybrid‐“stress function” plane element. Design/methodology/approach − First, a complementary energy functional, in which the Airy stress function is taken as the functional variable, is established within an element for analysis of plane problems. Second, 15 basic analytical solutions (in global Cartesian coordinates) of the stress function are taken as the trial functions for an 8‐node element, and meanwhile, 15 unknown constants are then introduced. Third, according to the principle of minimum complementary energy, the unknown constants can be expressed in terms of the displacements along element edges, which are interpolated by element nodal displacements. Finally, the whole system can be rewritten in terms of element nodal displacement vector. Findings − A new hybrid element stiffness matrix is obtained. The resulting 8‐node plane element, denoted as analytical trial function (ATF‐Q8), possesses excellent...

A novel hybrid stress-function finite element method immune to severe mesh distortion

This paper introduces a hybrid stress-function finite element method proposed recently for developing 2D finite element models immune to element shapes. Deferent from the first version of the hybrid-stress element constructed by Pian, the stress function of 2D elastic or fracture problem is regarded as the functional variable of the complementary energy functional. Then, the basic analytical solutions of are taken as the trial functions for finite element models, and meanwhile, the corresponding unknown stress-function constants are introduced. By using the principle of minimum complementary energy, these unknown stress-function constants can be expressed in terms of the displacements along element edges. Finally, the complementary energy functional can be rewritten in terms of element nodal displacement vector, and thus, the element stiffness matrix of such hybrid-function element can be obtained. As examples, two (8- and 12-node) quadrilateral plane elements and an arbitrary polygonal crack element are constructed by employing different basic analytical solutions of different stress functions. Numerical results show that, the 8- and 12-node plane models can produce the exact solutions for pure bending and linear bending problems, respectively, even the element shape degenerates into triangle and concave quadrangle; and the crack element can also predict accurate results with very low computational cost in analysis of stress-singularity problems.