Computer Methods in Applied Mechanics and Engineering
A finite element method with discontinuous rotations for the Mindlin–Reissner plate model
Introduction
Plates are very common simplified models for thin structures in engineering practice. The most basic plate models are the Kirchhoff model, which is a fourth order partial differential equation, and the Mindlin–Reissner (MR) model, which is a system of second order partial differential equations. The Kirchhoff model can be seen as the limit of the MR model as the thickness of the plate tends to zero. Finite element approximations of plate models would seem to be easier to handle with the MR model, since then only C0 continuity is required, as opposed to the C1-continuous elements needed for the Kirchhoff model. However, in order for a finite element method to work asymptotically as t → 0 in the MR model, typically rather complicated approximations must be used.
In this paper, we will consider a family of simple continuous-discontinuous Galerkin finite element methods for the MR model, first proposed in [10], based on discontinuous piecewise polynomials for the discretization of the rotations and continuous piecewise polynomials of one degree higher for the transverse displacements.
When the thickness of the plate tends to zero we obtain the Kirchhoff plate and our scheme simplifies to the method proposed in [9]. In this context we also mention the discontinuous Galerkin methods for the Kirchhoff plate developed by Hansbo and Larson [12] and for the Mindlin–Reissner model by Bösing et al. [4].
Section snippets
The continuous problem
The Mindlin–Reissner plate model is described by the following partial differential equations:where u is the transverse displacement, θ is the rotation of the median surface, t is the thickness, assumed constant, t3g is the transverse surface load, andis the moment tensor. Here, 1 is the identity tensor and ε is the strain operator with componentsThe material constants are given by the relations
The finite element method
For simplicity, we shall consider the case of clamped boundary conditions. The transverse displacement and rotation vector are solutions to the following variational problem: find and such thatfor all . However, as is common in discontinuous Galerkin methods, we need to assume a higher regularity of the solution in order for all terms in our scheme to make sense with the exact solution inserted. In particular, we will need n · σ(θ
Stability estimates
For our analysis, we introduce the following edge norm:and mesh dependent energy-like normThe mesh dependent norm ⦀ · ⦀ can be used to bound the broken H1(Ω) norm on Θh, which is the statement of the following lemma. Lemma 3 There is a constant c, independent of h, μ, and λ such that Proof This is a discrete Korn-type inequality that results from the control of the rigid body rotations given by the
A priori error estimates
In this section, we will derive a priori error estimates for CDG methods in the case k = 2, and show that they hold uniformly in t. For higher order methods, edge effects will typically preclude global estimates because of the lack of regularity of the exact solution.
A duality-based posteriori estimate
For duality-based a posteriori error analysis, we consider the following variant of (19): find z and ψ such thatwith zero Dirichlet boundary conditions for z and ψ. With eθ = θ − θh and eu = u − uh we find, using Lemma 2, thatwhere πh and now denote arbitrary interpolants (or projections) onto the respective subspaces. Using the equilibrium equations we find that
Numerical examples
We apply Algorithm 1 to a set of simple model problems in order to: (1) exemplify the behavior of the adaptive procedure; (2) study convergence rates of the finite element method (6) with respect to meshsize and plate thickness; and (3) study how the choice of stability parameter affects the approximation on different types of meshes.
We remark that the solutions in Sections 7.1 An L-shaped membrane, 7.2 The unit square are not smooth enough to be covered by the a priori error analysis presented
Concluding remarks
We have presented a novel finite element method for the Mindlin–Reissner plate model, based on the discontinuous Galerkin approach. We show that our method does not lock as long as we make a proper choice of a free, but computable, parameter. Our approach avoids the current paradigm of projections of the rotations in the shear energy functional, which, at least from a conceptual point of view, requires a mixed implementation. We pay the price of having to use a higher number of degrees of
References (16)
- et al.
Continuous/discontinuous approximations of fourth-order elliptic problems in structural and continuum mechanics with applications to thin bending elements and strain gradient elasticity
Comput. Methods Appl. Mech. Engrg.
(2002) - et al.
Discontinuous Galerkin methods for incompressible and nearly incompressible elasticity by Nitsche’s method
Comput. Methods Appl. Mech. Engrg.
(2002) - et al.
A uniformly accurate finite-element method for the Reissner–Mindlin plate
SIAM J. Numer. Anal.
(1989) - et al.
A four-node plate bending element based on Mindlin–Reissner plate theory and mixed interpolation
Int. J. Numer. Methods Engrg.
(1985) - et al.
A feed-back approach to error control in finite element methods: basic analysis and examples
East–West J. Numer. Math.
(1996) - et al.
A new interior penalty discontinuous Galerkin method for the Reissner–Mindlin model
Math. Models Methods Appl. Sci.
(2010) Korn’s inequalities for piecewise H1 vector fields
Math. Comput.
(2004)- et al.
Error analysis of mixed-interpolated elements for Reissner–Mindlin plates
Math. Models Methods Appl. Sci.
(1991)
Cited by (27)
Embedded boundary conditions for shear-deformable plate bending
2023, Computer Methods in Applied Mechanics and EngineeringA locking-free weak Galerkin finite element method for Reissner–Mindlin plate on polygonal meshes
2020, Computers and Mathematics with ApplicationsC<sup>0</sup> interior penalty methods for a dynamic nonlinear beam model
2018, Applied Mathematics and ComputationA simple approach for finite element simulation of reinforced plates
2018, Finite Elements in Analysis and DesignCitation Excerpt :If the beams are located on one side of the plate the main approach is the same but we must include models for the membrane deformation of the plate and beams. The idea of using continuous/discontinuous approximations was first proposed by Engel et al. [6] and later analysed for Kirchhoff–Love and Mindlin–Reissner plates in Refs. [7,9,10], cf. also Wells and Dung [21]. Other boundary conditions for plates, for instance simply supported and free, can easily be included in the c/dG finite element method, see Ref. [8] for details.
C<sup>0</sup>-discontinuous Galerkin methods for a wind-driven ocean circulation model: Two-grid algorithm
2018, Computer Methods in Applied Mechanics and EngineeringA discontinuous Galerkin method for nonlinear shear-flexible shells
2016, Computer Methods in Applied Mechanics and EngineeringCitation Excerpt :The development of discontinuous Galerkin (DG) finite element methods has provided new means to attack locking problems, as evidenced by the formulations in [19–21] for nearly and fully incompressible elasticity and in [22] for Stokes flow. In particular, the last 10 years has seen the development of several locking-free DG methods for Reissner plates [23–26]. In [24–26], the flexibility in the choice of the approximation spaces afforded by discontinuous Galerkin methods is exploited to ensure satisfaction of the natural compatibility conditions of the thin plate limit.