Search results

The purpose of this paper is to give a review on the newest developments of high-performance finite element methods (FEMs), and exhibit the recent contributions achieved by the authors’ group, especially showing some breakthroughs against inherent difficulties existing in the traditional FEM for a long time.

Three kinds of new FEMs are emphasized and introduced, including the hybrid stress-function element method, the hybrid displacement-function element method for Mindlin–Reissner plate and the improved unsymmetric FEM. The distinguished feature of these three methods is that they all apply the fundamental analytical solutions of elasticity expressed in different coordinates as their trial functions.

The new FEMs show advantages from both analytical and numerical approaches. All the models exhibit outstanding capacity for resisting various severe mesh distortions, and even perform well when other models cannot work. Some difficulties in the history of FEM are also broken through, such as the limitations defined by MacNeal’s theorem and the edge-effect problems of Mindlin–Reissner plate.

These contributions possess high value for solving the difficulties in engineering computations, and promote the progress of FEM.

The purpose of this paper is to propose an efficient low-order quadrilateral flat shell element that possesses all outstanding advantages of novel shape-free plate bending and plane membrane elements proposed recently.

By assembling a shape-free quadrilateral hybrid displacement-function (HDF) plate bending element HDF-P4-11β (Cen et al. 2014) and a shape-free quadrilateral hybrid stress-function (HSF) plane membrane element HSF-Q4θ-7β (Cen et al. 2011b) with drilling degrees of freedom, a new 4-node, 24-DOF (6 DOFs per node) quadrilateral flat shell element is successfully constructed. The trial functions for resultant/stress fields within the element are derived from the analytical solutions of displacement and stress functions for plate bending and plane problems, respectively, so that they can a priori satisfy the related governing equations. Furthermore, in order to take the influences of moderately warping geometry into consideration, the rigid link correction strategy (Taylor 1987) is also employed.

The element stiffness matrix of a new simple 4-node, 24-DOF quadrilateral flat shell element is obtained. The resulting models, denoted as HDF-SH4, not only possesses all advantages of original HDF plate and HSF plane elements when analyzing plate and plane structures, but also exhibits good performances for analyses of complicated spatial shell structures. Especially, it is quite insensitive to mesh distortions.

This work presents a new scheme, which can take the advantages of both analytical and discrete methods, to develop low-order mesh-distortion resistant flat shell elements.

A simple shape-free high-order hybrid displacement function element method is presented for precise bending analyses of Mindlin–Reissner plates. Three distortion-resistant and locking-free eight-node plate elements are proposed by utilizing this method.

This method is based on the principle of minimum complementary energy, in which the trial functions for resultant fields are derived from two displacement functions, F and f, and satisfy all governing equations. Meanwhile, the element boundary displacements are determined by the locking-free arbitrary order Timoshenko’s beam functions. Then, three locking-free eight-node, 24-DOF quadrilateral plate-bending elements are formulated: HDF-P8-23β for general cases, HDF-P8-SS1 for edge effects along soft simply supported (SS1) boundary and HDF-P8-FREE for edge effects along free boundary.

The proposed elements can pass all patch tests, exhibit excellent convergence and possess superior precision when compared to all other existing eight-node models, and can still provide good and stable results even when extremely coarse and distorted meshes are used. They can also effectively solve the edge effect by accurately capturing the peak value and the dramatical variations of resultants near the SS1 and free boundaries. The proposed eight-node models possess potential in engineering applications and can be easily integrated into commercial software.

This work presents a new scheme, which can take the advantages of both analytical and discrete methods, to develop high-order mesh distortion-resistant Mindlin–Reissner plate-bending elements.

This paper aims to propose a simple but robust three-node triangular membrane element with rational drilling DOFs for efficiently analyzing plane problems.

This new element is developed within the general framework of unsymmetric FEM. The element test functions are determined by using a conforming displacement field which is slightly different with the classical Allman’s interpolations, while a self-equilibrated stress field formulated based on the analytical airy stress solutions is adopted as the trial functions. To ensure the correctness between the drilling DOFs and the true rotations in elasticity, reasonable constraints are introduced through the penalty function method. Moreover, the special quadrature strategy is used for operating related integrations for future enrichment of element behavior.

Numerical benchmark tests reveal that this new triangular membrane element has exceptional prediction capabilities. In particular, this element can correctly reproduce a rigid body rotation motion and correctly undertake the external in-plane twisting moments; thus, it is a reasonable choice for being used to formulate flat shell elements or to be connected with other kind of elements with physical rotational DOFs.

This work provides a new approach for developing high-performance lower-order elements with simple formulations and good numerical accuracies.

The purpose of this paper is to present a two-dimensional (2D) model of the base force element method (BFEM) based on the complementary energy principle. The study proposes a model of the BFEM for arbitrary mesh problems.

The BFEM uses the base forces given by Gao (2003) as fundamental variables to describe the stress state of an elastic system. An explicit expression of element compliance matrix is derived using the concept of base forces. The detailed formulations of governing equations for the BFEM are given using the Lagrange multiplier method. The explicit displacement expression of nodes is given. To verify the 2D model, a program on the BFEM using MATLAB language is made and a number of examples on arbitrary polygonal meshes and aberrant meshes are provided to illustrate the BFEM.

A good agreement is obtained between the numerical and theoretical results. Based on the studies, it is found that the 2D formulation of BFEM with complementary energy principle provides reliable predictions for arbitrary mesh problems.

Due to the use of Lagrange multiplier method, there are more basic unknowns in the control equation. The proposed method will be improved in the future.

This paper presents a new idea and a new numerical method, and to explore new ways to solve the problem of arbitrary meshes.

The paper presents a 2D model of the BFEM using the complementary energy principle for arbitrary mesh problems.