Plate spectral elements based upon Reissner–Mindlin theory

Plate bending spectral elements based upon the Reissner–Mindlin theory is developed for a rectangular plate. Numerical comparison studies for static and eigenvalue problems, demonstrate the excellent performance and robustness of the present method. The improvement in accuracy and computational efficiency in comparison with low order plate bending elements is remarkable and suggests that the proposed approach is an alternative to existing schemes.     

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.     

On stabilized finite element methods for the Reissner–Mindlin plate model

Stabilized finite element formulation for the Reissner–Mindlin plate model is considered. Physical interpretation for the stabilization procedure for low‐order elements is established. Explicit interpolation functions for linear and bilinear stabilized MITC elements are derived. Some numerical examples including buckling and frequency analyses are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

On the accuracy of p-version elements for the Reissner–Mindlin plate problem

This paper addresses the question of accuracy of p-version finite element formulations for Reissner–Mindlin plate problems. Three model problems, a circular arc, a rhombic plate and a geometrically complex structure are investigated. Whereas displacements and bending moments turn out to be very accurate without any post-processing even for very coarse meshes, the quality of shear forces computed from constitutive equations is poor. It is shown that significantly improved results can be obtained, if shear forces are computed from equilibrium equations instead. A consistent computation of second derivatives of the shape functions is derived. © 1998 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.     

On evaluating the inf–sup condition for plate bending elements

This paper addresses the evaluation of the inf–sup condition for Reissner–Mindlin plate bending elements. This fundamental condition for stability and optimality of a mixed finite element scheme is, in general, very difficult to evaluate analytically, considering for example distorted meshes. Therefore, we develop a numerical test methodology. To demonstrate the test methodology and to obtain specific results, we apply it to standard displacement-based elements and elements of the MITC family. Whereas the displacement-based elements fail to satisfy the inf–sup condition, we find that the MITC elements pass our numerical test for uniform meshes and a sequence of distorted meshes. © 1997 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.     

On the evaluation of the method of finite spheres for the solution of Reissner–Mindlin plate problems using the numerical inf–sup test

In this paper we use the numerical inf–sup test to evaluate both displacement‐based and mixed discretization schemes for the solution of Reissner–Mindlin plate problems using the meshfree method of finite spheres. While an analytical proof of whether a discretization scheme passes the inf–sup condition is most desirable, such a proof is usually out of reach due to the complexity of the meshfree approximation spaces involved. The numerical inf–sup test (Int. J. Numer. Meth. Engng 1997; 40 :3639–3663), developed to test finite element discretization spaces, has therefore been adopted in this paper. Tests have been performed for both regular and irregular nodal configurations. While, like linear finite elements, pure displacement‐based approximation spaces with linear consistency do not pass the inf–sup test and exhibit shear locking, quadratic discretizations, unlike quadratic finite elements, pass the test. Pure displacement‐based and mixed approximation spaces that pass the numerical inf–sup test exhibit optimal or near optimal convergence behaviour. Copyright © 2007 John Wiley & Sons, Ltd.     

An adaptive non‐conforming finite‐element method for Reissner–Mindlin plates

Adaptive algorithms are important tools for efficient finite‐element mesh design. In this paper, an error controlled adaptive mesh‐refining algorithm is proposed for a non‐conforming low‐order finite‐element method for the Reissner–Mindlin plate model. The algorithm is controlled by a reliable and efficient residual‐based a posteriori error estimate, which is robust with respect to the plate's thickness. Numerical evidence for this and the efficiency of the new algorithm is provided in the sense that non‐optimal convergence rates are optimally improved in our numerical experiments. Copyright © 2003 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.     

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.     

A Chebyshev collocation multidomain method to solve the Reissner–Mindlin equations for the transient response of an anisotropic plate subjected to impact

The transient response of an anisotropic rectangular plate subjected to impact is described through a Chebyshev collocation multidomain discretization of the Reissner–Mindlin plate equations. The trapezoidal rule is used for time-integration. The spatial collocation derivative operators are represented by matrices, and the subdomains are patched by natural and essential conditions. At each time level the resulting governing matrix equation is reduced by two consecutive block Gaussian eliminations, so that an equation for the variables at the subdomain corners has to be solved. Back-substitution gives the variables at all other collocation points. The time history as represented by computed contour plots has been compared with analytical results and with photos produced by holographic interferometry. The agreements are satisfactory. © 1997 John Wiley & Sons, Ltd.     

Automatic adaptive refinement for shell analysis using nine-node assumed strain element

An automatic adaptive refinement procedure for the analysis of shell structures using the nine-node degenerated solid shell element is suggested. The basic adaptive refinement principle and the effects of singularities and boundary layers on the convergence rate of the nine-node element used are discussed. A new stress recovery procedure based on the patch convective co-ordinate system concept is developed for the construction of a continuous smoothed stress field over the shell domains. The stress recovery procedure is easy to implement, requires a modest computational effort and needs only local patch information. It can be applied to shells with non-uniform thickness as well as to multi-layered shell structures. The smoothed recovered stress obtained is then used with the Zienkiewicz and Zhu error estimator for a posteriori error estimation during the adaptive refinement analysis. Numerical results which are in good agreement with theoretical predictions are obtained and they indicate that the current adaptive refinement procedure can eliminate the effect of singularities inside the problem domains so that a near-optimal convergence rate is achieved in all the numerical examples. This also indicates that the stress recovery procedure can produce an accurate stress field and as a result the error estimator can reflect the error distribution of the finite element solution. Even though in the current study only one type of element is used in the analysis, the whole adaptive refinement scheme can be readily applied to any other types of degenerated solid element. © 1997 John Wiley & Sons, Ltd.     

The partition of unity finite element method for elastic wave propagation in Reissner–Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four‐noded element for the Reissner–Mindlin plate is derived in this paper, which is free of shear locking. Such a locking‐free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency‐dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright © 2006 John Wiley & Sons, Ltd.     

Multilayered shell finite element with interlaminar continuous shear stresses: a refinement of the Reissner–Mindlin formulation

A finite element formulation for refined linear analysis of multilayered shell structures of moderate thickness is presented. An underlying shell model is a direct extension of the first‐order shear‐deformation theory of Reissner–Mindlin type. A refined theory with seven unknown kinematic fields is developed: (i) by introducing an assumption of a zig‐zag (i.e. layer‐wise linear) variation of displacement field through the thickness, and (ii) by assuming an independent transverse shear stress fields in each layer in the framework of Reissner's mixed variational principle. The introduced transverse shear stress unknowns are eliminated on the cross‐section level. At this process, the interlaminar equilibrium conditions (i.e. the interlaminar shear stress continuity conditions) are imposed. As a result, the weak form of constitutive equations (the so‐called weak form of Hooke's law) is obtained for the transverse strains–transverse stress resultants relation. A finite element approximation is based on the four‐noded isoparametric element. To eliminate the shear locking effect, the assumed strain variational concept is used. Performance of the derived finite element is illustrated with some numerical examples. The results are compared with the exact three‐dimensional solutions, as well as with the analytical and numerical solutions obtained by the classical, the first‐order and some representative refined models. Copyright © 2000 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.     

Implicit dynamic analysis using an isogeometric Reissner–Mindlin shell formulation

In isogeometric analysis, identical basis functions are used for geometrical representation and analysis. In this work, non‐uniform rational basis splines basis functions are applied in an isoparametric approach. An isogeometric Reissner–Mindlin shell formulation for implicit dynamic calculations using the Galerkin method is presented. A consistent as well as a lumped matrix formulation is implemented. The suitability of the developed shell formulation for natural frequency analysis is demonstrated by a numerical example. In a second set of examples, transient problems of plane and curved geometries undergoing large deformations in combination with nonlinear material behavior are investigated. Via a zero‐thickness stress algorithm for arbitrary material models, a J₂‐plasticity constitutive law is implemented. In the numerical examples, the effectiveness, robustness, and superior accuracy of a continuous interpolation method of the shell director vector is compared with experimental results and alternative numerical approaches. Copyright © 2016 John Wiley & Sons, Ltd.     

Static,free vibration,and buckling analysis of laminated composite Reissner–Mindlin plates using NURBS‐based isogeometric approach

This paper presents a novel numerical procedure based on the framework of isogeometric analysis for static, free vibration, and buckling analysis of laminated composite plates using the first‐order shear deformation theory. The isogeometric approach utilizes non‐uniform rational B‐splines to implement for the quadratic, cubic, and quartic elements. Shear locking problem still exists in the stiffness formulation, and hence, it can be significantly alleviated by a stabilization technique. Several numerical examples are presented to show the performance of the method, and the results obtained are compared with other available ones. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Reissner–Mindlin shell implementation and energy conserving isogeometric multi‐patch coupling

A shear‐flexible isogeometric Reissner–Mindlin shell formulation using non‐uniform rational B‐splines basis functions is introduced, which is used for the demonstration of a coupling approach for multiple non‐conforming patches. The six degrees of freedom formulation uses the exact surface normal vectors and curvature. The shell formulation is implemented in an isogeometric analysis framework for computation of structures composed of multiple geometric entities. To enable local model refinement as well as non‐matching domains coupling, a conservative multi‐patch approach using Lagrange multipliers for structured non‐uniform rational B‐splines patches is presented. Here, an additional border frame mesh is used to couple geometries during structural analyses. This frame interface approach avoids the problem of excessive constraints when multiple patches are coupled at one point. First, the shell formulation is verified with several reference cases. Then the influence of the frame interface discretization and frame penalty stiffness on the smoothness of the results is investigated. The effects of the perturbed Lagrangian method in combination with the frame interface approach is shown. In addition, results of models with T‐joint interface connections and perpendicular stiffener patches are presented. Copyright © 2016 John Wiley & Sons, Ltd.