Numerical analysis of a locking-free mixed finite element method for a bending moment formulation of Reissner-Mindlin plate model
|
|
- Oliver Bradley
- 5 years ago
- Views:
Transcription
1 Numerical analysis of a locking-free mixed finite element method for a bending moment formulation of Reissner-Mindlin plate model LOURENÇO BEIRÃO DA VEIGA Dipartimento di Matematica F. Enriques, Università degli Studi di Milano, Via Saldini 50, 2033 Milano, Italy. DAVID MORA Departamento de Matemática, Facultad de Ciencias, Universidad del Bío Bío, Casilla 5-C, Concepción, Chile. AND RODOLFO RODRÍGUEZ CI 2 MA, Departamento de Ingeniería Matemática, Universidad de Concepción, Casilla 60-C, Concepción, Chile. This paper deals with the approximation of the bending of a clamped plate modeled by Reissner-Mindlin equations. It is known that standard finite element methods applied to this model lead to wrong results when the thickness t is small. Here, we propose a new mixed formulation in terms of the bending moments, shear stress, rotations and transversal displacement. To prove that the resulting variational formulation is well posed, we use standard Babuška-Brezzi theory with appropriate t-dependent norms. The problem is discretized by standard finite elements and error estimates are proved with constants independent of the plate thickness. Moreover, these constants depend on norms of the solution that can be a priori bounded independently of the plate thickness, which leads to the conclusion that the method is locking-free. A local post-processing leading to H approximations of transversal displacement and rotations is introduced. Finally, we report numerical experiments confirming the theoretical results. Keywords: Reissner-Mindlin; bending moment formulation; locking-free finite elements; error analysis.. Introduction The Reissner-Mindlin theory is the most used model to approximate the deformation of an elastic thin or moderately thick plate. It is very well understood that standard finite element methods applied to this model lead to wrong results when the thickness t is small with respect to the other dimensions of the plate, due to the so called locking phenomenon. Nevertheless, adopting for instance a reduced integration or a mixed interpolation technique, this phenomenon can be avoided. Indeed, nowadays several families of methods have been rigorously shown to be free from locking and optimally convergent. We mention the recent monograph by Falk (2008) for a thorough description of the state of the art and further references. Among the existing methods, a large success has been shared by the mixed interpolation of tensorial components (MITC) methods introduced in Bathe & Dvorkin (985) or variants of them (Durán lourenco.beirao@unimi.it dmora@ubiobio.cl rodolfo@ing-mat.udec.cl
2 2 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ & Liberman (992)). Another solution is to write an equivalent formulation between the plate equations and an uncoupled system of two Poisson equations plus a rotated Stokes system, by means of a Helmholtz decomposition of the shear stress, as in Arnold & Falk (989). More recently, Amara et al. (2002) proposed and analyzed a conforming finite element method for the Reissner-Mindlin model satisfying various boundary conditions. In their analysis the bending moment is written in terms of three auxiliary variables belonging to classical Sobolev spaces. A mixed formulation in terms of these new variables is discretized by standard finite elements. Under regularity assumption on the exact solution, optimal error estimates were proved with constants independent of the plate thickness. Another approach is presented by Behrens & Guzmán (2009). In this case the plate bending problem is posed in terms of six variables lying in L 2 and H(div) spaces. A discretization in terms of discontinuous polynomials and enriched Raviart-Thomas elements is proposed. Error estimates with t-independent constants are proved. These estimates are quasi optimal in regularity, since they involve a norm of the shear stress which can not be a priori bounded independently of t. In the present paper we consider a bending moment formulation for the plate problem. We introduce these moments (which in practice usually represents the quantity of true interest in applications) as a new unknown, together with the shear stress, the rotations and the transversal displacement. We obtain a mixed variational formulation (Elasticity-like system) which, thanks to Babuška-Brezzi theory, is shown to be well posed and stable in appropriate t-dependent norms. For the approximation of bending moment and rotations we employ PEERS finite elements introduced by Arnold et al. (984), classical Raviart- Thomas elements are used for the shear stress and piecewise constants for the transversal displacement. We prove an uniform inf-sup condition with respect to the discretization parameter h and the thickness t. Moreover, the convergence rate is proved to be optimal O(h). The obtained estimates only depend on norms of the quantities which are known to be bounded independent of t. Therefore, the method turns out thoroughly locking-free. In addition, we propose a local post-processing procedure which gives piecewise linear rotations and transversal displacement which converge to the exact solution in a stronger H type discrete norm. The outline of this paper is as follows: In Section 2, we first recall the Reissner-Mindlin equations and some regularity results. Then, we prove the unique solvability and stability properties of the proposed formulation. In Section 3, we present the finite element scheme, prove a stability result and show the (linear) convergence of the method. In addition, we introduce and analyze a local post-processing procedure for transversal displacements and rotations. Finally, in Section 4 we report numerical tests which allow us to assess the performance of the proposed method. Throughout the paper we will use standard notations for Sobolev spaces, norms and seminorms. Moreover, we will denote with c and C, with or without subscripts, tildes, or hats a generic constant independent of the mesh parameter h and the plate thickness t, which may take different values in different occurrences. From now on, we use the following notation for any tensor field τ = (τ i j ) i, j =,2, any vector field η = (η i ) i =,2 and any scalar field v: divη := η + 2 η 2, rotη := η 2 2 η, v := ( ) ( ) v 2 v, curlv :=, 2 v v ( ) ( ) ( ) τ div τ := + 2 τ 2 2 η, Curlη := η η, η := 2 η, τ τ 22 2 η 2 η 2 η 2 2 η 2
3 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 3 of 20 Moreover, we denote 2. The plate model. τ t := (τ ji ), tr(τ) := I := 2 i= ( ) 0, J := 0 τ ii, τ a : τ b := ( ) i, j= τ a i jτ b i j. Consider an elastic plate of thickness t such that 0 < t, with reference configuration ( t 2, t 2), where is a convex polygonal domain of R 2 occupied by the midsection of the plate. The deformation of the plate is described by means of the Reissner-Mindlin model in terms of the rotations β = (β,β 2 ) of the fibers initially normal to the plate s midsurface, the scaled shear stress γ = (γ,γ 2 ), and the transversal displacement w. Assuming that the plate is clamped on its whole boundary, the following strong equations describe the plate s response to conveniently scaled transversal load g L 2 (): div(c(ε(β))) γ = 0 in, (2.) divγ = g in, (2.2) γ = κ t2( w β) in, (2.3) w = 0, β = 0 on, (2.4) where κ := Ek/2(+ν) is the shear modulus, with E being the Young modulus, ν the Poisson ratio, and k a correction factor usually taken as 5/6 for clamped plates, ε(β) := 2 ( β +( β)t ) is the standard strain tensor, and C is the tensor of bending moduli, given by (for isotropic materials) C τ := E 2( ν 2 ) (( ν)τ + νtr(τ)i) τ [L2 ()] 2 2. The tensor C is invertible with its inverse given by C τ := 2( ( ) ν2 ) E ( ν) τ ν ( ν 2 ) tr(τ)i τ [L 2 ()] 2 2. To write a variational formulation of the Reissner-Mindlin plate problem, we introduce the bending moment σ = (σ i j ), i, j =,2 as a new unknown defined by We can rewrite the equation above as: σ := C (ε(β)). C σ = β + ( ) 2 rotβ J, introducing the auxiliary unknown r := 2 rotβ. Multiplying by test function and then integrating by parts, we get C σ : τ + β div τ + r(τ 2 τ 2 ) = 0. (2.5)
4 4 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ Now, by testing (2.)-(2.3) with adequate functions, integrating by parts, using (2.5) and (2.4), and imposing weakly the symmetry of σ, we obtain the following mixed variational formulation: Find ((σ,γ),(β,r,w)) H Q such that C σ : τ + t2 κ γ ξ + for all ((τ,ξ),(η,s,v)) H Q, where with and where β (div τ + ξ)+ η (div σ + γ)+ r(τ 2 τ 2 )+ s(σ 2 σ 2 )+ H := H(div;) H(div;), Q := [L 2 ()] 2 L 2 () L 2 (), H(div;) := {τ [L 2 ()] 2 2 : div τ [L 2 ()] 2 }, H(div;) := {ξ [L 2 ()] 2 : divξ L 2 ()}. wdivξ =0, vdivγ = In the analysis, we will utilize the following t-dependent norm for the space H while for the space Q, we will use (τ,ξ) H := τ 0, + div τ + ξ 0, + ξ t,h(div;), ξ t,h(div;) := t ξ 0, + divξ 0,, (η,s,v) Q := η 0, + s 0, + v 0,. We note that for all (τ,ξ) H(div;) H(div;), where C is independent of t and (τ,ξ) H C( τ H(div;) + ξ H(div;) ), τ H(div;) + ξ H(div;) := τ 0, + div τ 0, + ξ 0, + divξ 0,. We rewrite problem (2.6) as follows: Find ((σ,γ),(β,r,w)) H Q such that gv, (2.6) a((σ,γ),(τ,ξ))+b((τ,ξ),(β,r,w)) = 0 (τ,ξ) H, b((σ,γ),(η,s,v)) = F(η,s,v) (η,s,v) Q, (2.7) where the bilinear forms a : H H R and b : H Q R, and the linear functional F : Q R, are defined by a((σ,γ),(τ,ξ)) := C σ : τ + t2 κ = 2( ( ν2 ) E ν γ ξ σ : τ ν ) ν 2 tr(σ)tr(τ) + t2 γ ξ, κ (2.8)
5 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 5 of 20 b((τ,ξ),(η,s,v)) := η (div τ + ξ)+ s(τ 2 τ 2 )+ vdivξ, (2.9) and F(η,s,v) := gv, for all (σ,γ),(τ,ξ) H and for all (η,s,v) Q. We recall the following regularity result for the solution of problem (2.7) (see Arnold & Falk (989)). PROPOSITION 2. Suppose that is a convex polygon and g L 2 (). Let ((σ,γ),(β,r,w)) be the solution of problem (2.7). Then, there exists a constant C, independent of t and g, such that w 2, + β 2, + γ H(div;) + t γ, + σ, + t div σ, + r, C g 0,. Now, we will prove that problem (2.7) satisfies the hypotheses of the Babuška-Brezzi theory, which yields the unique solvability and continuous dependence of this variational formulation. We first observe that the bilinear forms a, b, and the linear functional F are bounded with constants independent of plate thickness t. Let V := {(τ,ξ) H : b((τ,ξ),(η,s,v)) = 0 (η,s,v) Q} be the null space of the bilinear form b. From (2.9), we have that V := {(τ,ξ) H : ξ + div τ = 0, τ = τ t and divξ = 0 in }. The following lemma shows that the bilinear form a is V-elliptic, with a constant independent of the plate thickness t. LEMMA 2. There exists C > 0, independent of t, such that a((τ,ξ),(τ,ξ)) C (τ,ξ) 2 H (τ,ξ) V. Proof. Given (τ,ξ) V, using tr(τ) 2 2(τ : τ) τ [L 2 ()] 2 2, from (2.8) we obtain a((τ,ξ),(τ,ξ)) 2( ν) τ 2 0, E + t2 κ ξ 2 0,. Since div τ + ξ 0, = 0 and divξ 0, = 0, we get { 2( ν) a((τ,ξ),(τ,ξ)) min, } ( τ 2 0, E κ + div τ + ξ 2 0, + t2 ξ 2 0, + divξ 0,) 2. Thus where a((τ,ξ),(τ,ξ)) C (τ,ξ) 2 H, { 6( ν) C := min, }. E 2κ Therefore, a is V-elliptic, and we end the proof. In order to prove the corresponding inf-sup condition, we first prove the following lemma. LEMMA 2.2 There exists c > 0, independent of t, such that the following holds. For all s L 2 (), there exists τ s H(div;) such that (τ s 2 τs 2 ) = s, div τs = 0 in, and τ s H(div;) c s 0,.
6 6 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ Proof. Let s L 2 (), let us put P(s) := s, and consider λ := (s P(s)). We have that λ L 2 0 () := {v L2 () : v = 0} and moreover λ 0, c s 0,. Then, there exists v [H0 ()]2 such that divv = λ in and v, c λ 0,. We consider the following function ϕ := v+ P(s) ( ) x, 2 y which satisties divϕ = s and ϕ, v, + c s 0,. Now, we define ( τ s := Curlϕ = 2 v v 2 P(s) ) 2 v P(s) v 2 From this, we get that div τ s = 0, hence τ s H(div;). Moreover [L 2 ()] 2 2. and it is easy to check that (τ2 s τs 2 ) = v + P(s) 2 + 2v 2 + P(s) = divv+p(s) = λ + P(s) = s, 2 τ s H(div;) c s 0,. Thus, we end the proof. We are ready to prove the inf-sup condition for the bilinear form b. LEMMA 2.3 There exists C > 0, independent of t, such that b((τ,ξ),(η,s,v)) sup C (η,s,v) Q (η,s,v) Q. (τ,ξ) H (τ,ξ) H Proof. Let (η,s,v) Q. From Lemma 2.2, there exists τ s H(div;) such that div τ s = 0 in, (τ s 2 τs 2 ) = s in and τs H(div;) c s 0,. Then, b((τ,ξ),(η,s,v)) sup b((τs, 0),(η,s,v)) (τ,ξ) H (τ,ξ) H τ s 0, + div τ s 0, = τ s η div τ s + s(τ2 s τs 2 ) 0, = s 2 0, τ s 0, c s 0,. (2.0) Now, let τ := ε(z), where z [H 0 ()]2 is the unique weak solution (as a consequence of Korn s inequality and Lax-Milgram s lemma) of the following auxiliary problem { divε(z) = η in, z = 0 on.
7 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 7 of 20 We have that div τ = η [L 2 ()] 2, hence τ H(div;) and τ = τ t. Moreover, applying the continuous dependence result, there exists c > 0 depending only on such that τ 0, + div τ 0, c η 0,. (2.) Therefore, b((τ,ξ),(η,s,v)) sup b(( τ, 0),(η,s,v)) (τ,ξ) H (τ,ξ) H τ 0, + div τ 0, = τ 0, + div τ η div τ + s( τ 2 τ 2 ) 0, = η 2 0, τ 0, + div τ 0, c η 0,, (2.2) where the last inequality follows from (2.). Finally, let ξ := z, where z H0 () is the unique weak solution (as a consequence of Poincaré inequality and Lax-Milgram s lemma) of the following auxiliary problem { z = v in, z = 0 on. As before, there exists ĉ > 0 depending only on such that ξ 0, + div ξ 0, ĉ v 0,. It follows that b((τ,ξ),(η,s,v)) b((0, ξ),(η,s,v)) sup (τ,ξ) H (τ,ξ) H (+ t) ξ 0, + div ξ 0, 2( ξ 0, + div ξ η ξ + vdiv ξ 0, ) = 2( ξ 0, + div ξ η ξ + v 2 0, 0, ) ) ( v 20, 2( ξ 0, + div ξ 0, ) η ξ 2ĉ v 0, η 0,. Given we have proved that (cf. (2.2)) b((τ,ξ),(η,s,v)) A := sup, (τ,ξ) H (τ,ξ) H A c η 0,,
8 8 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ therefore, A 2ĉ(+ c) v 0,. (2.3) Thus, the proof follows from (2.0), (2.2) and (2.3). We are now in position to state the main result of this section which give the solvability of the continuous problem (2.7). THEOREM 2.2 There exists a unique ((σ,γ),(β,r,w)) H Q solution of the mixed variational formulation (2.7), and the following continuous dependence result holds where C is independent of t. ((σ,γ),(β,r,w)) H Q C g 0,, Proof. By virtue of Lemmas 2. and 2.3, the proof follows from a straightforward application of the abstract Theorem. in Chapter II of Brezzi & Fortin (99). 3. The Finite Element Scheme Let T h be a regular family of triangulations of the polygonal region by triangles T of diameter h T with mesh size h := max{h T : T T h }, and such that there holds = {T : T T h }. In addition, given an integer k 0 and a subset S of R 2, we denote by P k (S) the space of polynomials in two variables defined in S of total degree at most k, and for each T T h we define the local Raviart-Thomas space of order zero ( x where y ) is a generic vector of R 2. {( RT 0 (T) := span 0 ) ( 0, ) ( x, y On the other hand, for each triangle T T h we take the unique polynomial b T P 3 (T) that vanishes on T and is normalized by T b T =. This cubic bubble function is extended by zero onto the region T and therefore it becomes an element of H0 (). Hence, we define )}, B(T h ) := {τ h H(div,) : (τ ih,τ i2h ) Z(T),i =,2, T T h }, where Z(T) := span{curl(b T ),T T h }. Hence, we define the following finite element subspaces: H σ h := X h B(T h ), where X h is the global Raviart-Thomas space of lowest order, X h := {τ h H(div,) : τ h T [RT 0 (T) t ] 2, T T h }, H γ h := {ξ h H(div,) : ξ h T RT 0 (T), T T h }, Q w h := {v h L 2 () : v h T P 0 (T), T T h },
9 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 9 of 20 Q β h := {η h [L 2 ()] 2 : η h T [P 0 (T)] 2, T T h }, Q r h := { s h H () : s h T P (T), T T h }. At this point we recall that H σ h Q β h Qr h correspond to the PEERS-space given by Arnold, Brezzi and Douglas in Arnold et al. (984). Defining H h := H σ h H γ h and Q h := Q β h Qr h Qw h our mixed finite element scheme associated with the continuous formulation (2.7) reads as follows: Find ((σ h,γ h ),(β h,r h,w h )) H h Q h such that a((σ h,γ h ),(τ h,ξ h ))+b((τ h,ξ h ),(β h,r h,w h )) = 0 (τ h,ξ h ) H h, b((σ h,γ h ),(η h,s h,v h )) = F(η h,s h,v h ) (η h,s h,v h ) Q h. (3.) Our next goal is to show the corresponding discrete versions of Lemmas 2. and 2.3 to have the solvability and stability of problem (3.). With this aim, we note that the discrete null space of the bilinear form b reduces to: { V h := (τ h,ξ h ) H h : η h (div τ h + ξ h )+ s h (τ 2h τ 2h )+ v h divξ h = 0 (η h,s h,v h ) Q h }. Let (τ h,ξ h ) V h, taking (0,0,v h ) Q h and using that (divξ h ) T is a constant, since v h T is also a constant, we conclude that divξ h = 0 in. Now, taking (η h,0,0) Q h, since div τ h = 0 in τ h B(T h ), we have that (div τ h ) T is a constant vector. Moreover, since divξ h = 0, we have that ξ h T is also a constant vector. Therefore, since η h T is also a constant vector, we conclude that (div τ h + ξ h ) = 0 in. Thus, we obtain { V h = (τ h,ξ h ) H h : ξ h + div τ h = 0, s h (τ 2h τ 2h ) = 0 } s h Q r h and divξ h = 0 in. Note that the second condition in the above definition does not guarantee the symmetry of the tensors H σ h, as it was the case for the continuous kernel of b. Hence, we have that V h is not included in V, however the proof of Lemma 2. can be repeated (because in this proof we never used that τ = τ t ), to get the following result LEMMA 3. There exists C > 0 such that where the constant C is independent of h and t. a((τ h,ξ h ),(τ h,ξ h )) C (τ h,ξ h ) 2 H (τ h,ξ h ) V h, We introduce the Raviart-Thomas interpolation operator R : [H ()] 2 H γ h. Let us review some properties of this operator that we will use in the sequel: Let P be the orthogonal projection from L 2 () onto the finite element subspace Q w h. Then for all ξ [H ()] 2, we have divrξ = P(divξ). (3.2)
10 0 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ There exists c > 0, independent of h, such that ξ Rξ 0, ch ξ, ξ [H ()] 2. (3.3) Now, let Π h : [H ()] 2 2 X h be the usual interpolation operator defined as the cartesian product of Raviart-Thomas interpolation operator R, which satisfies (see Arnold et al. (984)): Let P be the orthogonal projection from [L 2 ()] 2 onto the finite element subspace Q β h. Then for all τ [H ()] 2 2, we have divπ h τ = P(div τ). There exists c > 0, independent of h, such that τ Π h τ 0, ch τ, τ [H ()] 2 2. Since X h H σ h, the operator Π h can be considered from [H ()] 2 2 into H σ h, with the same properties given above. Moreover, we let P : L 2 () Q r h the ortogonal projection. Then, we have s P (s) 0, Ch s,. (3.4) We continue with the following lemma establishing the discrete analogue of Lemma 2.3. LEMMA 3.2 There exists C > 0, independent of h and t, such that b((τ h,ξ h ),(η h,s h,v h )) sup C (η h,s h,v h ) Q (η h,s h,v h ) Q h. (τ h,ξ h ) H h (τ h,ξ h ) H Proof. Let (η h,s h,v h ) Q h. From Lemma 4.4 in Arnold et al. (984), we have that there exists τ h H σ h and c > 0 such that, η h div τ h + s h ( τ 2h τ 2h ) c { } η h 0, + s h 0,. (3.5) τ h 0, + div τ h 0, Thus, using (3.5), we have which yields = b((τ h,ξ h ),(η h,s h,v h ) sup b(( τ h, 0),(η h,s h,v h ) (τ h,ξ h ) H h (τ h,ξ h ) H τ h 0, + div τ h 0, η h div τ h + s h ( τ 2h τ 2h ) c { } η h 0, + s h 0,, τ h 0, + div τ h 0, sup (τ h,ξ h ) H h b((τ h,ξ h ),(η h,s h,v h )) (τ h,ξ h ) H c { η h 0, + s h 0, } (ηh,s h,v h ) Q h. (3.6)
11 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES of 20 Let now z be the unique weak solution (as a consequence of Poincaré inequality and Lax-Milgram s lemma) of the following problem: { z = vh in, z = 0 on. Since v h L 2 () and is a convex domain, a classical elliptic regularity result guarantees that z H 2 () H 0 () and there exists c > 0 such that z 2, c v h 0,. Now, we define ξ := z [H ()] 2. We note that div ξ = v h in, and ξ, = z, z 2, c v h 0,. (3.7) Let ξ h := R ξ. From (3.2) and the fact that div ξ = v h, we have that div ξ h = v h in. Hence, using the approximation property (3.3), we deduce that ξ h 0, + div ξ h 0, ξ h ξ 0, + ξ 0, + div ξ 0, ch ξ, + ξ,. From estimates (3.7), we obtain It follows that Given we have proved that (cf. (3.6)) and therefore we have that (+ t) ξ h 0, + div ξ h 0, 2( ξ h 0, + div ξ h 0, ) ĉ v h 0,. b((τ h,ξ h ),(η h,s h,v h )) sup b((0, ξ h ),(η h,s h,v h )) (τ h,ξ h ) H h (τ h,ξ h ) H (+ t) ξ h 0, + div ξ h 0, η h ξ h + v h div ξ h 2( ξ h 0, + div ξ h 0, ) 2( ξ h 0, + div ξ η h ξ h + v h 2 0, h 0, ) ) ( v 2( ξ h 0, + div ξ h 20, h 0, ) η h ξ h A h := ĉ v h 0, η h 0,. b((τ h,ξ h ),(η h,s h,v h )) sup, (τ h,ξ h ) H h (τ h,ξ h ) H A h A h c η h 0,, c (+ c)ĉ v h 0,. This allows us to conclude the proof. We are now in a position to establish the unique solvability, the stability, and the convergence properties of the discrete problem (3.).
12 2 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ THEOREM 3. There exists a unique ((σ h,γ h ),(β h,r h,w h )) H h Q h solution of the discrete problem (3.). Moreover, there exist C,C > 0, independent of h and t, such that and ((σ h,γ h ),(β h,r h,w h )) H Q C g 0,, ((σ,γ),(β,r,w)) ((σ h,γ h ),(β h,r h,w h )) H Q C inf ((σ,γ),(β,r,w)) ((τ h,ξ h ),(η h,s h,v h )) H Q, ((τ h,ξ h ),(η h,s h,v h )) H h Q h (3.8) where ((σ, γ),(β, r, w)) H Q is the unique solution of the mixed variational formulation (2.7). Proof. Is a direct application of the Theorem 2. in Chapter II of Brezzi & Fortin (99). The following theorem provides the rate of convergence of our mixed finite element scheme (3.). THEOREM 3.2 Let ((σ,γ),(β,r,w)) H Q and ((σ h,γ h ),(β h,r h,w h )) H h Q h be the unique solutions of the continuous and discrete problems (2.7) and (3.), respectively. If g H (), then, Proof. First, we note that (see (2.)) ((σ,γ),(β,r,w)) ((σ h,γ h ),(β h,r h,w h )) H Q Ch g,. (div σ + γ) (div σ h + γ h ) 0, = div σ h + γ h 0,. From the second equation of (3.), we have that η h (div σ h + γ h ) = 0 η h Q β h, hence Thus, also recalling (2.2), Therefore (div σ h ) T = γ h = T P(γ h ). T div σ h + γ h 2 0, = T T h div σ h + γ h 2 0,T = On the other hand, we have γ h P(γ h ) 2 0,T C h 2 T γ h 2,T T T h T T h C h 2 T divγ h 2 0,T = C h 2 T P(divγ) 2 0,T Ch2 g 2 0,. T T h T T h div σ h + γ h 0, Ch g 0,. (3.9) divγ divγ h 2 0, = T T h divγ divγ h 2 0,T = T T h divγ P(divγ) 2 0,T C T T h h 2 T divγ 2,T, which, using (2.2), yields divγ divγ h 0, Ch g,. (3.0)
13 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 3 of 20 Now, it is easy to check that (see Lemma 2.) ( τ 2 0, + t2 ξ 2 ) 0, Ca((τ,ξ),(τ,ξ)) (τ,ξ) H. In particular, taking (σ σ h,γ γ h ) H, we get ( σ σ h 2 0, + t2 γ γ h 2 ) 0, Ca((σ σ h,γ γ h ),(σ σ h,γ γ h )), and, using the first equation of (2.7), we obtain ( σ σ h 2 0, + t2 γ γ h 2 ) 0, Cb((σ h σ,γ h γ),(β,r,w)). (3.) Now, from the definition of the bilinear form b(, ) we get b((σ h σ,γ h γ),(β,r,w)) = β (div(σ h σ)+(γ h γ)) + r((σ 2h σ 2h ) (σ 2 σ 2 ))+ wdiv(γ h γ). (3.2) Subtracting the second equation of (3.) from the second equation (2.7), we have that η h (div(σ h σ)+(γ h γ)) = 0 η h Q β h, (3.3) s h ((σ 2h σ 2h ) (σ 2 σ 2 )) = 0 s h Q r h, (3.4) v h div(γ h γ) = 0 v h Q w h. (3.5) Considering P(β) Q β h, P (r) Q r h and P(w) Qw h, using (3.3), (3.4) and (3.5), we rewrite (3.2) as follow b((σ h σ,γ h γ),(β,r,w)) = (β P(β)) (div(σ h σ)+(γ h γ)) + (r P (r))((σ 2h σ 2h ) (σ 2 σ 2 )) + (w P(w))div(γ h γ). From (3.), the above equation and the Cauchy-Schwarz inequality, we obtain σ σ h 2 0, + t2 γ γ h 2 0, C β P(β) 0, div(σ h σ)+(γ h γ) 0, +C r P (r) 0, σ h σ 0, +C w P(w) 0, div(γ h γ) 0,. Applying the inequality pq 2 p2 + 2 q2, from the above bound it follows σ σ h 2 0, + t2 γ γ h 2 0, C β P(β) 0, div(σ h σ)+(γ h γ) 0, + 2 C2 r P (r) 2 0, + 2 σ h σ 2 0, +C w P(w) 0, div(γ h γ) 0,.
14 4 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ Using standard error estimates arguments, (3.4), (3.9) and (3.0), yields σ σ h 2 0, + t2 γ γ h 2 0, Ĉ( ch 2 β, g 0, + ch 2 r, + ch 2 w, g, ), and thus, from Proposition 2., we get Finally, using (3.8) we obtain ((σ,γ),(β,r,w)) ((σ h,γ h ),(β h,r h,w h )) H Q σ σ h 0, + t γ γ h 0, Ch g,. (3.6) C ((σ,γ),(β,r,w)) ((σ h,γ h ),( P(β),P (r),p(w))) H Q Ch g,, where in the last inequality we have used standard error estimates arguments, (3.4), (3.9), (3.0) and (3.6). We conclude the proof. 3. A post processing of transversal displacement and rotations In this section we present an element-wise post processing procedure which allows to build piecewise linear transversal displacement and rotations with improved approximation properties. In the following, we indicate with e a general edge of the triangulation and with E h the set of all such edges. Moreover, we indicate with h e the length of e E h and associate to each edge a unit normal vector n e, chosen once and for all. For each internal edge e of E h, we indicate with T + and T the two triangles of the mesh which have the edge e in common, where n e corresponds to the outward normal for T + and the opposite for T. Then, given any piecewise regular (scalar or vector) function v on, for each e E h we define the jump on internal edges [[v]] = v + e v e, where v ± is the restriction of v to T ±. On boundary edges, the jump is simply given by the value of v on the edge. We introduce the following H type discrete norm v 2,h = v 2 0,T + h e [[v]] 2 0,e, T T h e E h for all sufficiently regular (scalar or vector) functions v. Given, the discrete solution ((σ h,γ h ),(β h,r h,w h )), we define a post-processed transversal displacement w h L2 () as follows. For all T T h let w h P (T) such that Pw h = w h, w h = P(β h + t 2 κ γ h ). (3.7) It is immediate to check that w h is well defined and unique. Let w I = Pw, we start proving the following preliminary result. LEMMA 3.3 There holds w I w h,h Ch g,.
15 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 5 of 20 Proof. To prove the result, we will apply the following inf-sup condition: For all v h Q w h, there exists ξ h H γ h such that v h divξ h = v h 2,h, ξ h 0, C v h,h. (3.8) The simple proof of the above inf-sup condition will be shown briefly. Defining the degrees of freedom of ξ h by ξ h n e := h e [[v h ]] for all e E h, an element-wise integration by parts and the definition of the jump operator yield v h divξ h = h e [[v h ]] 2 0,e = v h 2,h. e E h The second bound in (3.8) follows easily by a scaling argument. Applying (3.8) to v h = (w h w I ), noting that divξ h is piecewise constant and finally using the discrete equations (3.), we obtain w I w h 2,h = (w I w h )divξ h = (w w h )divξ h = (β β h )ξ h + t 2 κ (γ γ h )ξ h. The proof then follows from the above equation using a Cauchy-Schwarz inequality, recalling Theorem 3.2 and using (3.8). We have the following improved convergence result for the post-processed transversal displacement. PROPOSITION 3.3 There holds Proof. We note that we can split w w h,h Ch g,. w h = w h + w h, w = w I + w, (3.9) where w h and w I, already defined above, are piecewise constant while w h and w have zero average on each element. Applying a scaled trace inequality on each triangle T and using that w h and w have zero average on each element, we have that h e e E h [[ w h w]] 2 0,e C (h 2 T T T h w h w 2 0,T + w h w 2,T) C ( w h w) 2 0,T. (3.20) T T h We now observe that, due to (3.9), there hold w h T = w h T and w T = w T for all T T h. Therefore, first due to definition (3.7) and (2.3), then using standard properties of the projector P, for all T T h we obtain ( w h w) 2 0,T = P(β h + t 2 κ γ h ) (β + t 2 κ γ) 2 0,T P(β h + t 2 κ γ h ) P(β + t 2 κ γ) 2 0,T + P(β + t 2 κ γ) (β + t 2 κ γ) 2 0,T (β h + t 2 κ γ h ) (β + t 2 κ γ) 2 0,T +Ch2 T (β + t2 κ γ) 2,T C( β h β 2 0,T + t4 γ h γ 2 0,T + h2 T β 2,T + h2 T t4 γ 2,T ). The above estimate, combined with Theorem 3.2 and Proposition 2., immediately yields T T h ( w h w) 2 0,T Ch2 g 2,. (3.2)
16 6 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ From (3.20), (3.2) and the definition of,h norm, we finally obtain w h w,h Ch g,. The above estimate, combined with Lemma 3.3 and a triangle inequality, finally gives the proof of the Proposition: w w h,h w I w h,h + w h w,h Ch g,. We define also a post-processed rotation field βh [L2 ()] 2 as follows: For all T T h let βh [P (T)] 2 such that Pβ h = β h, βh = P(C (3.22) σ h + r h J), where P is the L 2 projection onto the space of piecewise constant R 2 2 tensor fields. It is immediate to check that βh is well defined and unique. Moreover, the following results can be proved following the same lines shown above. PROPOSITION 3.4 There holds β β h,h Ch g,. Finally note that both post-processing procedures are fully local and therefore have a negligible computational cost. REMARK 3. Although the main purpose of this scheme is to compute a better approximation of the stresses, using this post-processing a piecewise linear approximation of transversal displacement and rotations, converging in a H type norm can be recovered. Note in particular that, from the definition of the norm,h and the fact that the jumps of w and β are null, it follows that at the limit for h 0 the post-processed discrete functions will also be continuous. 4. Numerical results We report in this section some numerical experiments which confirm the theoretical results proved above. The numerical method analyzed has been implemented in a MATLAB code. As a test problem we have taken an isotropic and homogeneous plate := (0,) (0,) clamped on the whole boundary for which the analytical solution is explicitly known (see Chinosi et al. (2006)). We analyze the convergence properties of the elements proposed here by considering different uniform decompositions as shown in Figure, and keeping the thickness t = Choosing the transversal load g as: E g(x,y) = [2y(y 2( ν 2 )(5x 2 5x+) ( 2y 2 (y ) 2 + x(x )(5y 2 5y+) ) ) the exact solution of problem (2.7) is given by w(x,y) = 3 x3 (x ) 3 y 3 (y ) 3 +2x(x )(5y 2 5y+) ( 2x 2 (x ) 2 + y(y )(5x 2 5x+) )], 2t2 [ ] y 3 (y ) 3 x(x )(5x 2 5x+)+x 3 (x ) 3 y(y )(5y 2 5y+), 5( ν)
17 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 7 of 20 FIG.. Square plate: uniform meshes. β (x,y) = y 3 (y ) 3 x 2 (x ) 2 (2x ), β 2 (x,y) = x 3 (x ) 3 y 2 (y ) 2 (2y ). The material constants have been chosen: E =, ν = 0.30 and the shear correction factor has been taken k = 5/6. In what follows, N denotes the number of degrees of freedom, namely, N := dim(h h Q h ). Moreover, we define the individual errors by: e(σ) := σ σ h 0,, e(σ,γ) := (div σ + γ) (div σ h + γ h ) 0,, e(γ) := γ γ h t,h(div;), e(r) := r r h 0,, e(β) := β β h 0,, e(w) := w w h 0,, where ((σ,γ),(β,r,w)) H Q and ((σ h,γ h ),(β h,r h,w h )) H h Q h are the unique solutions of problems (2.7) and (3.), respectively. Also, we define the experimental rates of convergence (rc) for the errors e(σ), e(σ, γ), e(γ), e(r), e(β) and e(w) by rc( ) := 2 log(e( )/e ( )) log(n/n, ) where N and N denote the degrees of freedom of two consecutive triangulations with errors e and e. Tables, 2 and 3 show the convergence history of the mixed finite element scheme (3.) applied to our test problem. Table. Errors and experimental rates of convergence for variables σ and (div σ + γ), computed on uniform meshes. N e(σ) rc(σ) e(σ,γ) rc(σ,γ) e e e e e e e e e e
18 8 of 20 L. BEIRÃO DA VEIGA, D. MORA AND R. RODRÍGUEZ Table 2. Errors and experimental rates of convergence for variables γ and r, computed on uniform meshes. N e(γ) rc(γ) e(r) rc(r) e e e e e e e e e e Table 3. Errors and experimental rates of convergence for variables β and w, computed on uniform meshes. N e(β) rc(β) e(w) rc(w) e e e e e e e e e e We observe there that the rate of convergence O(h) predicted by Theorem 3.2 is attained for all variables. Figure 2 shows the profiles of the discrete transversal displacement w h (left) and the first component of the discrete rotation vector β h (right) for t = 0.00, and the finest mesh. FIG. 2. Approximate transversal displacement (left) and the first component of the rotation vector (right). Figure 3 shows the profiles of the discrete shear stress γ h for t = 0.00, and the finest mesh.
19 FEM FOR BENDING MOMENT FORMULATION OF REISSNER-MINDLIN PLATES 9 of 20 F IG. 3. Approximate shear stress: first component (left) and second component (right). Figures 4 and 5 show the profiles of the bending moment tensor σ h = (σ i jh ), i, j =, 2, for t = 0.00, and the finest mesh. F IG. 4. Approximate bending moment: σ h (left) and σ 2h (right).
20 20 of 20 L. BEIRA O DA VEIGA, D. MORA AND R. RODRI GUEZ F IG. 5. Approximate bending moment: σ 2h (left) and σ 22h (right). Funding FONDAP and BASAL projects CMM, Universidad de Chile (Chile) to L.B. and R.R.; CONICYT (Chile) to D.M. R EFERENCES A MARA, M., C APATINA -PAPAGHIUC, D. & C HATTI, A. (2002) New locking-free mixed method for the ReissnerMindlin thin plate model, SIAM J. Numer. Anal., 40, A RNOLD, D. N., B REZZI, F. & D OUGLAS, J. (984) PEERS: A new mixed finite element for the plane elasticity, Japan J. Appl. Math.,, A RNOLD, D. N. & FALK, R. S. (989) A uniformly accurate finite element method for the Reissner-Mindlin plate, SIAM J. Numer. Anal., 26, BATHE, K. J. & DVORKIN, E. N. (985) A four-node plate bending element based on Mindlin-Reissner plate theory and a mixed interpolation, Int. J. Numer. Methods Engrg. 2, B EHRENS, E. M. & G UZM A N, J. (2009) A new family of mixed methods for the Reissner-Mindlin plate model based on a system of first-order equations, Scientific Computing Group Reports ( ), Brown University. B REZZI, F. & F ORTIN, M. (99) Mixed and Hybrid Finite Element Methods. New York: Springer. C HINOSI, C., L OVADINA, C. & M ARINI, L. D. (2006) Nonconforming locking-free finite elements for ReissnerMindlin plates, Comput. Methods Appl. Mech. Engrg., 95, D UR A N, R. & L IBERMAN, E. (992) On mixed finite elements methods for the Reissner-Mindlin plate model, Math. Comp., 58, FALK, R. (2008) Finite Elements for the Reissner-Mindlin Plate, in Mixed Finite Elements, Compatibility Conditions, and Applications, D. Boffi and L. Gastaldi, eds., Berlin: Springer,
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES
ANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES F. Auricchio Dipartimento di Meccanica Strutturale Università di Pavia, Italy C. Lovadina Dipartimento di Ingegneria Meccanica e Strutturale
More informationA UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE*
SIAM J. NUMER. ANAL. c 989 Society for Industrial and Applied Mathematics Vol. 26, No. 6, pp. 5, December 989 000 A UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE* DOUGLAS N. ARNOLD
More informationANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company ANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL DOUGLAS N. ARNOLD Department of
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationTomás P. Barrios 1, Gabriel N. Gatica 2,María González 3 and Norbert Heuer Introduction
ESAIM: M2AN Vol. 40, N o 5, 2006, pp. 843 869 DOI: 10.1051/m2an:2006036 ESAIM: Mathematical Modelling Numerical Analysis www.edpsciences.org/m2an A RESIDUAL BASED A POSERIORI ERROR ESIMAOR FOR AN AUGMENED
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationA NEW CLASS OF MIXED FINITE ELEMENT METHODS FOR REISSNER MINDLIN PLATES
A NEW CLASS OF IXED FINITE ELEENT ETHODS FOR REISSNER INDLIN PLATES C. LOVADINA Abstract. A new class of finite elements for Reissner indlin plate problem is presented. The family is based on a modified
More informationA-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements
A-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements Joint work with Lourenco Beirão da Veiga (Milano), Claudia Chinosi (Alessandria) and Carlo Lovadina (Pavia) Previous
More informationTwo Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate
Two Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate Pingbing Ming and Zhong-ci Shi Institute of Computational Mathematics & Scientific/Engineering Computing, AMSS, Chinese Academy of
More informationANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT METHOD FOR REISSNER-MINDLIN PLATES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 1, Pages 48 62 c 2017 Institute for Scientific Computing and Information ANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationNumerische Mathematik
Numer. Math. 2015 130:395 423 DOI 10.1007/s00211-014-0672-7 Numerische Mathematik Weakly over-penalized discontinuous Galerkin schemes for Reissner Mindlin plates without the shear variable Paulo Rafael
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationLeast-squares Finite Element Approximations for the Reissner Mindlin Plate
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl., 6, 479 496 999 Least-squares Finite Element Approximations for the Reissner Mindlin Plate Zhiqiang Cai, Xiu Ye 2 and Huilong Zhang
More informationLocking phenomena in Computational Mechanics: nearly incompressible materials and plate problems
Locking phenomena in Computational Mechanics: nearly incompressible materials and plate problems C. Lovadina Dipartimento di Matematica Univ. di Pavia IMATI-CNR, Pavia Bologna, September, the 18th 2006
More informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
More informationH(div) and H(curl)-conforming Virtual Element Methods
Noname manuscript No. (will be inserted by the editor) H(div) and H(curl)-conforming Virtual Element Methods L. Beirão da Veiga F. Brezzi L. D. Marini A. Russo the date of receipt and acceptance should
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationMIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS
MIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS DOUGLAS N. ARNOLD, RICHARD S. FALK, AND JAY GOPALAKRISHNAN Abstract. We consider the finite element solution
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationDISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES
DISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES MAR AINSWORTH, JOHNNY GUZMÁN, AND FRANCISCO JAVIER SAYAS Abstract. The existence of uniformly bounded discrete extension
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationA MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationDominique Chapelle 1, Anca Ferent 1 and Patrick Le Tallec 2
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o, 2003, pp. 43 58 DOI: 0.05/m2an:200305 THE TREATMENT OF PINCHING LOCKING IN 3D-SHELL
More informationIMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS. Maksim Frolov
IMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS Maksim Frolov Peter the Great St.Petersburg Polytechnic University AANMPDE11, FINLAND 06 10 August
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationA posteriori error analysis of an augmented mixed finite element method for Darcy flow
A posteriori error analysis of an augmented mixed finite element method for Darcy flow Tomás P. Barrios, J. Manuel Cascón and María González Abstract We develop an a posteriori error analysis of residual
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for Kirchhoff plate elements Jarkko Niiranen Laboratory of Structural
More informationA SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY
A SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY J. GOPALAKRISHNAN AND J. GUZMÁN Abstract. We presented a family of finite elements that use a polynomial space augmented
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
More informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationMORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
More informationarxiv: v1 [math.na] 19 Dec 2017
Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite
More informationA DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
More informationERROR ESTIMATES FOR LOW-ORDER ISOPARAMETRIC QUADRILATERAL FINITE ELEMENTS FOR PLATES
ERROR ESTIATES FOR LOW-ORDER ISOPARAETRIC QUADRILATERAL FINITE ELEENTS FOR PLATES RICARDO G DURÁN, ERWIN HERNÁNDEZ, LUIS HERVELLA-NIETO, ELSA LIBERAN, AND RODOLFO RODRíGUEZ Abstract This paper deals with
More informationA UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL
A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL C. CARSTENSEN Abstract. Residual-based a posteriori error estimates are derived within a unified setting for lowest-order conforming, nonconforming,
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationENHANCED STRAIN METHODS FOR ELASTICITY PROBLEMS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationDivergence free Virtual Elements for the Stokes problem on polygonal meshes
Divergence free Virtual Elements for the Stokes problem on polygonal meshes L. Beirão da Veiga, C. Lovadina and G. Vacca Abstract arxiv:1510.01655v1 [math.na] 6 Oct 2015 In the present paper we develop
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationANALYSIS OF SOME LOW ORDER QUADRILATERAL REISSNER-MINDLIN PLATE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718(XX)0000-0 ANALYSIS OF SOME LOW ORDER QUADRILATERAL REISSNER-MINDLIN PLATE ELEMENTS PINGBING MING AND ZHONG-CI SHI Abstract. Four
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationA Second Elasticity Element Using the Matrix Bubble
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2011 A Second Elasticity Element Using the
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationROBUST BDDC PRECONDITIONERS FOR REISSNER-MINDLIN PLATE BENDING PROBLEMS AND MITC ELEMENTS
ROBUST BDDC PRECONDITIONERS FOR REISSNER-MINDLIN PLATE BENDING PROBLEMS AND MITC ELEMENTS L. BEIRÃO DA VEIGA, C. CHINOSI, C. LOVADINA, AND L. F. PAVARINO Abstract. A Balancing Domain Decomposition Method
More informationMIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY
MIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY MARIE E. ROGNES AND RAGNAR WINTHER Abstract. Small deformations of a viscoelastic body are considered through the linear Maxwell
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationFinite element exterior calculus: A new approach to the stability of finite elements
Finite element exterior calculus: A new approach to the stability of finite elements Douglas N. Arnold joint with R. Falk and R. Winther Institute for Mathematics and its Applications University of Minnesota
More informationFinite element approximation on quadrilateral meshes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001; 17:805 812 (DOI: 10.1002/cnm.450) Finite element approximation on quadrilateral meshes Douglas N. Arnold 1;, Daniele
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationMauricio A. Barrientos 1, Gabriel N. Gatica 1 and Matthias Maischak 2
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 36, N o 2, 2002, pp. 241 272 DOI: 10.1051/m2an:2002011 A-POSTERIORI ERROR ESTIMATES FOR
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationAdvances in Engineering Software
Advances in Engineering Software 41 (21) 712 728 Contents lists available at ScienceDirect Advances in Engineering Software journal homepage: www.elsevier.com/locate/advengsoft The quadratic MITC plate
More informationA Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear Elasticity Equations
www.oeaw.ac.at A Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear lasticity quations B. Ayuso de Dios, I. Georgiev, J. Kraus, L. Zikatanov RICAM-Report 2011-29 www.ricam.oeaw.ac.at
More informationAdaptive approximation of eigenproblems: multiple eigenvalues and clusters
Adaptive approximation of eigenproblems: multiple eigenvalues and clusters Francesca Gardini Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/gardini Banff, July 1-6,
More information