Enhanced mixed interpolation XFEM formulations for discontinuous Timoshenko beam and Mindlin‐Reissner plate

Shear locking is a major issue emerging in the computational formulation of beam and plate finite elements of minimal number of degrees of freedom as it leads to artificial overstiffening. In this paper, discontinuous Timoshenko beam and Mindlin‐Reissner plate elements are developed by adopting the Hellinger‐Reissner functional with the displacements and through‐thickness shear strains as degrees of freedom. Heterogeneous beams and plates with weak discontinuity are considered, and the mixed formulation has been combined with the extended finite element method (FEM); thus, mixed enrichment functions are used. Both the displacement and the shear strain fields are enriched as opposed to the traditional extended FEM where only the displacement functions are enriched. The enrichment type is restricted to extrinsic mesh‐based topological local enrichment. The results from the proposed formulation correlate well with analytical solution in the case of the beam and in the case of the Mindlin‐Reissner plate with those of a finite element package (ABAQUS) and classical FEM and show higher rates of convergence. In all cases, the proposed method captures strain discontinuity accurately. Thus, the proposed method provides an accurate and a computationally more efficient way for the formulation of beam and plate finite elements of minimal number of degrees of freedom.     

Limit analysis and convex programming: A decomposition approach of the kinematic mixed method

This paper proposes an original decomposition approach to the upper bound method of limit analysis. It is based on a mixed finite element approach and on a convex interior point solver using linear or quadratic discontinuous velocity fields. Presented in plane strain, this method appears to be rapidly convergent, as verified in the Tresca compressed bar problem in the linear velocity case. Then, using discontinuous quadratic velocity fields, the method is applied to the celebrated problem of the stability factor of a Tresca vertical slope: the upper bound is lowered to 3.7776—value to be compared with the best published lower bound 3.772—by succeeding in solving non‐linear optimization problems with millions of variables and constraints. Copyright © 2008 John Wiley & Sons, Ltd.     

A new hybrid-mixed variational approach for Reissner–Mindlin plates. The MiSP model

A valuable variational approach for plate problems based on the Reissner–Mindlin theory is presented. The new MiSP (Mixed Shear Projected) approach is based on the Hellinger–Reissner variational principle, with a particular representation of transversal shear forces and transversal shear strains. The approximations of the shear forces are derived from those of the bending moments using the corresponding equilibrium relations. The shear strains are defined in terms of the edge tangential strains that are projected on the element degrees of freedom. Two finite elements are developed on the MiSP approach basis: 3-node triangular element MiSP3 and 4-node quadrilateral element MiSP4. Both elements can be considered as the most simple among the existent mixed elements. A modified MiSP model with a derived 4-node element is also presented. Numerical experiments are presented which show that the MiSP elements do not exhibit shear locking and give excellent results for thick and thin plates. They also pass the patch test for a general triangle and quadrilateral. © 1998 John Wiley & Sons, Ltd.     

On the performance of non‐conforming finite elements for the upper bound limit analysis of plates

In this paper, the upper bound limit analysis of thin plates in bending is addressed using various types of triangular finite elements for the generation of velocity fields and second‐order cone programming for the minimization problem. Three different C¹‐discontinuous finite elements are considered: the quadratic six‐node Lagrange triangle (T6), an enhanced T6 element with a cubic bubble function at centroid (T6b) and the cubic Hermite triangle (H3). Through numerical examples involving Johansen and von Mises yield criteria, it is shown that cubic elements (H3) give far better results in terms of convergence rate and precision than fully conforming elements found in the literature, especially for problems involving clamped boundaries. Copyright © 2013 John Wiley & Sons, Ltd.     

A new eight-node quadrilateral shear-bending plate finite element

A new eight-node quadrilateral shear-bending Reissner–Mindlin plate finite element for the very thin and thick plates without locking and spurious zero-energy modes is presented. The element has very good convergence characteristics both for thin and thick plates, is hardly insensitive to mesh distortions, and passes the patch tests. The formulation of the element is derived from a displacement variational principle and some general criteria to compute inconsistent transverse shear strains. These criteria have been applied with success to four- and eight-node quadrilateral plate finite elements and could be applied to construct triangular elements. The eight-node quadrilateral shear-bending plate finite element proposed has been found to be very efficient.     

A stabilized discrete shear gap finite element for adaptive limit analysis of Mindlin–Reissner plates

This paper presents a numerical formulation for computation of collapse load of Mindlin–Reissner plates that uses stabilized discrete shear gap finite elements and second‐order cone programming. Displacement fields are approximated using the discrete shear gap in combination with a stabilized strain smoothing technique, ensuring that shear‐locking problem can be avoided and that accurate solutions can be obtained. The underlying optimization problem is formulated in the form of a standard second‐order cone programming, so that it can be solved using highly efficient primal‐dual interior‐point algorithm. An error indicator based on plastic dissipation will be used in the adaptive refinement scheme. Various plates with arbitrary geometries and boundary conditions are examined to illustrate the performance of the proposed procedure. Copyright © 2013 John Wiley & Sons, Ltd.     

A stabilized discrete shear gap extended finite element for the analysis of cracked Reissner–Mindlin plate vibration problems involving distorted meshes

In this work, a new approach in the framework of the local partition of unity finite element method (XFEM) to significantly improve the accuracy of natural frequency in free vibration analysis of cracked Reissner–Mindlin plates is presented. Different from previous approaches, the present formulation is expected to be more accurate and effective in modeling cracked plates by integrating the stabilized discrete shear gap (DSG) into the XFEM setting. Intensive numerical results at low frequency demonstrated that the novel DSG-based XFEM approach possesses the following desirable properties: (1) the awkwardness of transverse shear-locking phenomenon can be overcome easily; (2) the DSG-based XFEM can be applicable to both moderately thick and thin plates straightforwardly; (3) the representation of cracks is independent of finite element mesh; (4) mesh distortion is insensitive and controllable; (5) the accuracy of natural frequency obtained by the present method is high and (6) the present method uses three-node triangular elements that can be much easily generated automatically for problems even with complicated geometry. These properties of the DSG-based XFEM are confirmed through several numerical examples of cracked plates with different boundary conditions.     

Lower bound static approach for the yield design of thick plates

The present work addresses the lower bound limit analysis (or yield design) of thick plates under shear‐bending interaction. Equilibrium finite elements are used to discretize the bending moment and the shear force fields. Different strength criteria, formulated in the five‐dimensional space of bending moment and shear force, are considered, one of them taking into account the interaction between bending and shear resistances. The criteria are chosen to be sufficiently simple so that the resulting optimization problem can be formulated as a second‐order cone programming problem (SOCP), which is solved by the dedicated solver MOSEK . The efficiency of the proposed finite element is illustrated by means of numerical examples on different plate geometries, for which the thin plate solutions as well as the pure shear solutions are accurately obtained as two different limit cases of the plate slenderness ratio. In particular, the proposed element exhibits a good behavior in the thin plate limit. Copyright © 2014 John Wiley & Sons, Ltd.     

Inelastic transient dynamic analysis of Reissner–Mindlin plates by the D/BEM

A direct domain/boundary element method (D/BEM) for dynamic analysis of elastoplastic Reissner–Mindlin plates in bending is developed. Thus, effects of shear deformation and rotatory inertia are included in the formulation. The method employs the elastostatic fundamental solution of the problem resulting in both boundary and domain integrals due to inertia and inelasticity. Thus, a boundary as well as a domain space discretization by means of quadratic boundary and interior elements is utilized. By using an explicit time‐integration scheme employed on the incremental form of the matrix equation of motion, the history of the plate dynamic response can be obtained. Numerical results for the forced vibration of elastoplastic Reissner–Mindlin plates with smooth boundaries subjected to impulsive loading are presented for illustrating the proposed method and demonstrating its merits. Copyright © 2000 John Wiley & Sons, Ltd.     

Eight‐node Reissner–Mindlin plate element based on boundary interpolation using Timoshenko beam function

A new eight‐node Reissner–Mindlin plate element is developed with a special interpolation within the element. This special interpolation is an extension of the element boundary interpolation that employs Timoshenko beam function for the boundary segment interpolation. The element function can effectively capture the structural behaviours of thick plates and achieve high precision in the analysis of thick plates. Patch tests and numerical investigations are conducted. It can be seen that the proposed element successfully passes all the patch tests. The results of the numerical investigation show that the proposed element is free of the shear locking phenomenon and possesses a higher accuracy in the analyses, as compared to the earlier research in the literature. Copyright © 2006 John Wiley & Sons, Ltd.     

Dispersion analysis of stabilized finite element methods for acoustic fluid interaction with Reissner–Mindlin plates

The application of stabilized finite element methods to model the vibration of elastic plates coupled with an acoustic fluid medium is considered. A complex‐wavenumber dispersion analysis of acoustic fluid interaction with Reissner–Mindlin plates is performed to quantify the accuracy of stabilized finite element methods for fluid‐loaded plates. Results demonstrate the improved accuracy of a recently developed hybrid least‐squares (HLS) plate element based on a modified Hellinger–Reissner functional, consistently combined with residual‐based methods for the acoustic fluid, compared to standard Galerkin and Galerkin gradient least‐squares plate elements. The technique of complex wavenumber dispersion analysis is used to examine the accuracy of the discretized system in the representation of free waves for fluid‐loaded plates. The influence of fluid and coupling matrices resulting from consistent implementation of pressure loading in the residual for the plate equation is examined and clarified for the different finite element approximations. Copyright © 2001 John Wiley & Sons, Ltd.     

C0 REISSNER–MINDLIN MULTILAYERED PLATE ELEMENTS INCLUDING ZIG-ZAG AND INTERLAMINAR STRESS CONTINUITY

Concerning composites plate theories and FEM (Finite Element Method) applications this paper presents some multilayered plate elements which meet computational requirements and include both the zig-zag distribution along the thickness co-ordinate of the in-plane displacements and the interlaminar continuity (equilibrium) for the transverse shear stresses. This is viewed as the extension to multilayered structures of well-known C⁰ Reissner–Mindlin finite plate elements. Two different fields along the plate thickness co-ordinate are assumed for the transverse shear stresses and for the displacements, respectively. In order to eliminate stress unknowns, reference is made to a Reissner mixed variational theorem. Sample tests have shown that the proposed elements, named RMZC, numerically work as the standard Reissner–Mindlin ones. Furthermore, comparisons with other results related to available higher-order shear deformation theories and to three-dimensional solutions have demonstrated the good performance of the RMZC elements.     

A new method for derivation of locking-free plate bending finite elements via mixed/hybrid formulation

The shear-locking phenomenon in discrete bending analysis of Mindlin/Reissner plates is investigated. Mixed/hybrid variational principles are introduced which, unlike the rigorous displacement model, allow systematic derivation of locking-free finite elements. This is achieved by satisfaction of an auxiliary condition, having the clear physical interpretation of shear-force elimination on account of equilibrium. An example, using competitive techniques, demonstrates the applicability of the idea.     

Refined 9‐Dof triangular Mindlin plate elements

Based on the Mindlin/Reissner plate theory, two refined triangular thin/thick plate elements, the conforming displacement element DKTM with one point quadrature for the part of shear strain and the element RDKTM with the re‐constitution of the shear strain, are proposed. In the formulations the exact displacement function of the Timoshenko's beam is used to derive the element displacements of the refined elements. Numerical examples are presented to show that the present models indeed possess properties of high accuracy for thin and thick plates, is capable of passing the patch test required for Kirchhoff thin plate elements, and does not exhibits extra zero energy modes. The element RDKTM is free of locking for very thin plate analysis and its convergence can be ensured theoretically. However, the element DKTM is not free of shear locking when the thickness/span ratios less than 10^?2. Copyright © 2001 John Wiley & Sons, Ltd.     

A stabilized one‐point integrated quadrilateral Reissner–Mindlin plate element

A new quadrilateral Reissner–Mindlin plate element with 12 element degrees of freedom is presented. For linear isotropic elasticity a Hellinger–Reissner functional with independent displacements, rotations and stress resultants is used. Within the mixed formulation the stress resultants are interpolated using five parameters for the bending moments and four parameters for the shear forces. The hybrid element stiffness matrix resulting from the stationary condition can be integrated analytically. This leads to a part obtained by one‐point integration and a stabilization matrix. The element possesses a correct rank, does not show shear locking and is applicable for the evaluation of displacements and stress resultants within the whole range of thin and thick plates. The bending patch test is fulfilled and the computed numerical examples show that the convergence behaviour is better than comparable quadrilateral assumed strain elements. Copyright © 2004 John Wiley & Sons, Ltd.     

Displacement‐based finite elements with nodal integration for Reissner–Mindlin plates

An assumed‐strain finite element technique is presented for shear‐deformable (Reissner–Mindlin) plates. The weighted residual method (reminiscent of the strain–displacement functional) is used to enforce weakly the balance equation with the natural boundary condition and, separately, the kinematic equation (the strain–displacement relationship). The a priori satisfaction of the kinematic weighted residual serves as a condition from which strain–displacement operators are derived via nodal integration, for linear triangles, and quadrilaterals, and also for quadratic triangles. The degrees of freedom are only the primitive variables: transverse displacements and rotations at the nodes. A straightforward constraint count can partially explain the insensitivity of the resulting finite element models to locking in the thin‐plate limit. We also construct an energy‐based argument for the ability of the present formulation to converge to the correct deflections in the limit of the thickness approaching zero. Examples are used to illustrate the performance with particular attention to the sensitivity to element shape and shear locking. Copyright © 2009 John Wiley & Sons, Ltd.     

Generalized mixed variational methods for Reissner plate and its applications

 Based on the mechanism of shear locking phenomenon and potential functional of Reissner plate bending problem, the generalized mixed variational principle for Reissner plate analysis is presented by parameterized Lagrange multiplier method. The proposed variational functional contains splitting factors which are able to adjust the shear potential energy and shear complementary energy components in it. The generalized mixed finite element formulation of bilinear quardrilateral element for Reissner plate bending analysis is established in terms of the new variational principle. The stiffness of the finite element model can be changed by the alteration of the splitting factors. Thus both the free of shear locking and higher accuracy are obtained by the choice of appropriate splitting factors. The most important is that this paper gives one self-adaptative way to choose the splitting factors for thin and moderately thick plates. This results in the comparative order of magnitude between the bending stiffness and shear stiffness for the arbitrary thickness. In the application of two-by-two exact Gaussian integration scheme to the proposed mixed element model, numerical examples show that free of locking is obtained even in the thin plate limit and high accuracy is given for moderately thick plate. Received: 15 January 2002 / Accepted: 10 September 2002 This work is partially supported by the National Nature Science Fund in China under Award No. 53978376     

A cell‐based smoothed discrete shear gap method using triangular elements for static and free vibration analyses of Reissner–Mindlin plates

The cell‐based strain smoothing technique is combined with discrete shear gap method using three‐node triangular elements to give a so‐called cell‐based smoothed discrete shear gap method (CS‐DSG3) for static and free vibration analyses of Reissner–Mindlin plates. In the process of formulating the system stiffness matrix of the CS‐DSG3, each triangular element will be divided into three subtriangles, and in each subtriangle, the stabilized discrete shear gap method is used to compute the strains and to avoid the transverse shear locking. Then the strain smoothing technique on whole the triangular element is used to smooth the strains on these three subtriangles. The numerical examples demonstrated that the CS‐DSG3 is free of shear locking, passes the patch test, and shows four superior properties such as: (1) being a strong competitor to many existing three‐node triangular plate elements in the static analysis; (2) can give high accurate solutions for problems with skew geometries in the static analysis; (3) can give high accurate solutions in free vibration analysis; and (4) can provide accurately the values of high frequencies of plates by using only coarse meshes. Copyright © 2012 John Wiley & Sons, Ltd.     

Hybrid displacement function element method: a simple hybrid‐Trefftz stress element method for analysis of Mindlin–Reissner plate

In order to develop robust finite element models for analysis of thin and moderately thick plates, a simple hybrid displacement function element method is presented. First, the variational functional of complementary energy for Mindlin–Reissner plates is modified to be expressed by a displacement function F, which can be used to derive displacement components satisfying all governing equations. Second, the assumed element resultant force fields, which can satisfy all related governing equations, are derived from the fundamental analytical solutions of F. Third, the displacements and shear strains along each element boundary are determined by the locking‐free formulae based on the Timoshenko's beam theory. Finally, by applying the principle of minimum complementary energy, the element stiffness matrix related to the conventional nodal displacement DOFs is obtained. Because the trial functions of the domain stress approximations a priori satisfy governing equations, this method is consistent with the hybrid‐Trefftz stress element method. As an example, a 4‐node, 12‐DOF quadrilateral plate bending element, HDF‐P4‐11 β, is formulated. Numerical benchmark examples have proved that the new model possesses excellent precision. It is also a shape‐free element that performs very well even when a severely distorted mesh containing concave quadrilateral and degenerated triangular elements is employed. Copyright © 2014 John Wiley & Sons, Ltd.     

An enhanced cell‐based smoothed finite element method for the analysis of Reissner–Mindlin plate bending problems involving distorted mesh

In this work, an enhanced cell‐based smoothed finite element method (FEM) is presented for the Reissner–Mindlin plate bending analysis. The smoothed curvature computed by a boundary integral along the boundaries of smoothing cells in original smoothed FEM is reformulated, and the relationship between the original approach and the present method in curvature smoothing is established. To improve the accuracy of shear strain in a distorted mesh, we span the shear strain space over the adjacent element. This is performed by employing an edge‐based smoothing technique through a simple area‐weighted smoothing procedure on MITC4 assumed shear strain field. A three‐field variational principle is utilized to develop the mixed formulation. The resultant element formulation is further reduced to a displacement‐based formulation via an assumed strain method defined by the edge‐smoothing technique. As the result, a new formulation consisting of smoothed curvature and smoothed shear strain interpolated by the standard transverse displacement/rotation fields and smoothing operators can be shown to improve the solution accuracy in cell‐based smoothed FEM for Reissner–Mindlin plate bending analysis. Several numerical examples are presented to demonstrate the accuracy of the proposed formulation.Copyright © 2013 John Wiley & Sons, Ltd.