arxiv:math/ v1 [math.na] 18 Jan 2007

Transcription

1 MIXED FINITE ELEMENT METHODS FOR LINEAR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY arxiv:mat/ v1 [mat.na] 18 Jan 2007 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Abstract. In tis paper, we construct new inite element metods or te approximation o te equations o linear elasticity in tree space dimensions tat produce direct approximations to bot stresses and displacements. Te metods are based on a modiied orm o te Hellinger Reissner variational principle tat only weakly imposes te symmetry condition on te stresses. Altoug tis approac as been previously used by a number o autors, a key new ingredient ere is a constructive derivation o te elasticity complex starting rom te de Ram complex. By mimicking tis construction in te discrete case, we derive new mixed inite elements or elasticity in a systematic manner rom known discretizations o te de Ram complex. Tese elements appear to be simpler tan te ones previously derived. For example, we construct stable discretizations wic use only piecewise linear elements to approximate te stress ield and piecewise constant unctions to approximate te displacement ield. 1. Introduction Te equations o linear elasticity can be written as a system o equations o te orm (1.1) Aσ = ǫu, divσ = in Ω. Here te unknowns σ and u denote te stress and displacement ields engendered by a body orce acting onalinearly elastic bodywic occupies a regionω R 3. Ten σ takes values in te space S := Rsym 3 3 o symmetric matrices and u takes values in V := R3. Te dierential operator ǫ is te symmetric part o te gradient, te div operator is applied row-wise to a matrix, and te compliance tensor A = A(x) : S S is a bounded and symmetric, uniormly positive deinite operator relecting te properties o te body. I te body is clamped on te boundary Ω o Ω, ten te proper boundary condition or te system (1.1) is u = 0 on Ω. For simplicity, tis boundary condition will be assumed ere. Te issues tat arise wen oter boundary conditions are assumed (e.g., te case o pure traction boundary conditions σn = g) are discussed in [9]. Te pair (σ,u) can alternatively be caracterized as te unique critical point o te Hellinger Reissner unctional (1 (1.2) J(τ,v) = 2 Aτ : τ +divτ v v) dx. Ω 2000 Matematics Subject Classiication. 65N30, 74S05. Key words and prases. mixed metod, inite element, elasticity. October 31, Revised September 1, Te work o te irst autor was supported in part by NSF grant DMS Te work o te second autor was supported in part by NSF grant DMS Te work o te tird autor was supported by te Norwegian Researc Council. 1

2 2 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Te critical point is sougt among all τ H(div,Ω;S), te space o square-integrable symmetric matrix ields wit square-integrable divergence, and all v L 2 (Ω;V), te space o square-integrable vector ields. Equivalently, (σ,u) H(div,Ω;S) L 2 (Ω;V) is te unique solution to te ollowing weak ormulation o te system (1.1): (Aσ : τ +divτ u)dx = 0, Ω (1.3) τ H(div,Ω;S), divσ vdx = vdx, v Ω Ω L2 (Ω;V). A mixed inite element metod determines an approximate stress ield σ and an approximate displacement ield u as te critical point o J over Σ V were Σ H(div,Ω;S) and V L 2 (Ω;V) are suitable piecewise polynomial subspaces. Equivalently, te pair (σ,u ) Σ V is determined by te weak ormulation (1.3), wit te test space restricted to Σ V. As is well known, te subspaces Σ and V cannot be cosen arbitrarily. To ensure tat a unique critical point exists and tat it provides a good approximation o te true solution, tey must satisy te stability conditions rom Brezzi s teory o mixed metods [12, 13]. Despite our decades o eort, no stable simple mixed inite element spaces or elasticity ave been constructed. For te corresponding problem in two space dimensions, stable inite elements were presented in [10]. For te lowest order element, te space Σ is composed o piecewise cubic unctions, wit 24 degrees o reedom per triangle, wile te space V consists o piecewise linear unctions. Anoter approac wic as been discussed in te two dimensional case is te use o composite elements, in wic V consists o piecewise polynomials wit respect to one triangulation o te domain, wile Σ consists o piecewise polynomials wit respect to a dierent, more reined, triangulation [5, 21, 23, 31]. In tree dimensions, a partial analogue o te element in [10] as been proposed and sown to be stable in [1]. Tis element uses piecewise quartic stresses wit 162 degrees o reedom per tetraedron, and piecewise linear displacements. Because o te lack o suitable mixed elasticity elements, several autors ave resorted to te use o Lagrangian unctionals wic are modiications o te Hellinger Reissner unctional given above [2, 4, 6, 27, 28, 29, 30], in wic te symmetry o te stress tensor is enorced only weakly or abandoned altogeter. In order to discuss tese metods, we consider te compliance tensor A(x) as a symmetric and positive deinite operator mapping M into M, were M is te space o 3 3 matrices. In te isotropic case, or example, te mapping σ Aσ as te orm Aσ = 1 ( λ σ 2µ 2µ+3λ tr(σ)i), were λ(x), µ(x) are positive scalar coeicients, te Lamé coeicients. A modiication o te variational principle discussed above is obtained i we consider te extended Hellinger Reissner unctional (1.4) J e (τ,v,q) = J(τ,v)+ τ : qdx over te space H(div,Ω;M) L 2 (Ω;V) L 2 (Ω;K), were K denotes te space o skew symmetric matrices. We note tat te symmetry condition or te space o matrix ields is now enorced troug te introduction o te Lagrange multiplier, q. A critical point (σ, u, p) Ω

3 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 3 o te unctional J e is caracterized as te unique solution o te system (Aσ : τ +divτ u+τ : p)dx = 0, Ω τ H(div,Ω;M), (1.5) divσ vdx = vdx, v Ω Ω L2 (Ω;V), σ : qdx = 0, q Ω L2 (Ω;K). Clearly, i (σ,u,p)is asolutiono tissystem, tenσ issymmetric, i.e., σ H(div,Ω;S), and tereore te pair (σ,u) H(div,Ω;S) L 2 (Ω;V) solves te corresponding system (1.3). On teoterand, i(u,p)solves(1.3), tenu H 1 (Ω;V)and, iwesetptoteskew-symmetric part o gradu, ten (σ,u,p) solves (1.5). In tis respect, te two systems (1.3) and (1.5) are equivalent. However, te extended system (1.5) leads to new possibilities or discretization. Assume tat we coose inite element spaces Σ V Q H(div,Ω;M) L 2 (Ω;V) L 2 (Ω;K) and consider a discrete system corresponding to (1.5). I (σ,u,p ) Σ V Q is a discrete solution, ten σ will not necessary inerit te symmetry property o σ. Instead, σ will satisy te weak symmetry condition σ : qdx = 0, or all q Q. Ω Tereore, tese solutions will in general not correspond to solutions o te discrete system obtained rom (1.3). Discretizations based on te system (1.5) will be reerred to as mixed inite element metods wit weakly imposed symmetry. Suc discretizations were already introduced by Fraejis de Veubeke in [21] and urter developed in [4]. In particular, te so called PEERS element proposed in [4] or te corresponding problem in two space dimensions used a combination o piecewise linear unctions and cubic bubble unctions, wit respect to a triangulation o te domain, to approximate te stress σ, piecewise constants to approximate te displacements, and continuous piecewise linear unctions to approximate te Lagrange multiplier p. Prior to te PEERS paper, Amara and Tomas [2] developed metods wit weakly imposed symmetry using a dual ybrid approac. Te lowest order metod tey discussed approximates te stresses wit quadratic polynomials plus bubble unctions and te multiplier by discontinuous constant or linear polynomials. Te displacements are approximated on boundary edges by linear unctions. Generalizations o te idea o weakly imposed symmetry to oter triangular elements, rectangular elements, and tree space dimensions were developed in [28], [29], [30] and [24]. In [29], a amily o elements is developed in bot two and tree dimensions. Te lowest order element in te amily uses quadratics plus te curls o quartic bubble unctions in two dimensions or quintic bubble unctions in tree dimensions to approximate te stresses, discontinuous linears to approximate te displacements, and discontinuous quadratics to approximate te multiplier. In addition, a lower order metod is introduced tat approximates te stress by piecewise linear unctions augmented by te curls o cubic bubble unctions plus a cubic bubble times te gradient o local rigid motions. Te multiplier is approximated by discontinuous piecewise linear unctions and te displacement by local rigid motions. Morley [24] extends PEERS to a amily o triangular elements, to rectangular elements, and to tree dimensions. In addition, te multiplier is approximated by nonconorming rater tan continuous piecewise polynomials.

4 4 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Tere is a close connection between mixed inite elements or linear elasticity and discretization o an associated dierential complex, te elasticity complex, wic will be introduced in 3 below. In act, te importance o tis complex was already recognized in [10], were mixed metods or elasticity in two space dimensions were discussed. Te new ingredient ere is tat we utilize a constructive derivation o te elasticity complex starting rom te de Ram complex. Tis construction is described in Eastwood [18] and is based on te te Bernstein Geland Geland resolution, c. [11] and also [14]. By mimicking te construction in te discrete case, we are able to derive new mixed inite elements or elasticity in a systematic manner rom known discretizations o te de Ram complex. As a result, we can construct new elements bot in two and tree space dimensions wic are signiicantly simpler tan tose derived previously. For example, we will construct stable discretizations o te system (1.5) wic only use piecewise linear and piecewise constant unctions, as illustrated in te igure below. For simplicity, te entire discussion o te present paper will be given in te tree dimensional case. A detailed discussion in two space dimensions can be ound in [8]. Besides te metods discussed ere, we note tat by sligtly generalizing te approac o tis paper, one can also analyze some o te previously known metods mentioned above tat are also based on te weak symmetry ormulation (see [19] or details). Figure 1. Elements or te stress, displacement, and multiplier in te lowest order case in two dimensions and tree dimensions. An alternative approac to construct inite element metods or linear elasticity is to consider a pure displacement ormulation. Since te coeicient A in (1.1) is invertible, te stress σ can be eliminated using te irst equation in (1.1), te stress strain relation. Tis leads to te second order equation (1.6) diva 1 ǫu = in Ω or te displacement u. A weak solution o tis equation can be caracterized as te global minimizer o te energy unctional (1 E(u) = 2 A 1 ǫu : ǫu+ u ) dx Ω over te Sobolev space H0(Ω;V). 1 Here H0(Ω;V) 1 denotes te space o all square integrable vector ields on Ω, wit square integrable derivatives, and wic vanis on te boundary Ω. A inite element approac based on tis ormulation, were we seek a minimum over a inite element subspace o H0 1 (Ω;V) is standard and discussed in textbooks, (e.g., [16]). However, or more general models, arising, or example, in viscoelasticity and plasticity (c. [15]), te stress strain relation is not local and an elimination o te stress σ is impossible. For suc models, a pure displacement model is excluded, and a mixed approac seems to be

5 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 5 an obvious alternative. Te construction o stable mixed elements or linear elasticity is an important step in te construction o mixed metods or tese more complicated models. Anoter advantage o te mixed approac is tat we automatically obtain scemes wic are uniormly stable in te incompressible limit, i.e., as te Lamé parameter λ tends to ininity. Since tis beavior o mixed metods is well known, we will not ocus urter on tis property ere. A more detailed discussion in tis direction can, or example, be ound in [5]. An outline o te paper is as ollows. In 2, we describe te notation to be used, state our main result, and provide some preliminary discussion on te relation between stability o mixed inite element metods and discrete exact complexes. In 3, we present two complexes related to te two mixed ormulations o elasticity given by (1.3) and (1.5). In 4, we introduce te ramework o dierential orms and sow ow te elasticity complex can be derived rom te de Ram complex. In 5, we derive discrete analogues o te elasticity complex beginning rom discrete analogues o te de Ram complex and identiy te required properties o te discrete spaces necessary or tis construction. Tis procedure is our basic design principle. In 6, we apply te construction o te preceding section to speciic discrete analogues o te de Ram complex to obtain a amily o discrete elasticity complexes. In 7 we use tis amily to construct stable inite element scemes or te approximation o te mixed ormulation o te equations o elasticity wit weakly imposed symmetry. Finally, in 8, we sow ow a sligtly more complicated procedure leads to a simpliied elasticity element. 2. Notation, statement o main results, and preliminaries We begin wit some basic notation and ypoteses. We continue to denote by V = R 3 te space o 3-vectors, by M te space o 3 3 real matrices, and by S and K te subspaces o symmetric and skew symmetric matrices, respectively. Te operators sym : M S and skw : M K denote te symmetric and skew symmetric parts, respectively. Note tat an element o te space K can be identiied wit its axial vector in V given by te map vec : K V: vec 0 v 3 v 2 v 3 0 v 1 = v 1 v 2, v 2 v 1 0 i.e., vec 1 (v)w = v w or any vectors v and w. WeassumetatΩisadomaininR 3 witboundary Ω. Wesallusetestandardunction spaces, like te Lebesgue space L 2 (Ω) and te Sobolev space H s (Ω). For vector-valued unctions, we include te range space in te notation ollowing a semicolon, so L 2 (Ω;X) denotes te space o square integrable unctions mapping Ω into a normed vector space X. Te space H(div,Ω;V) denotes te subspace o (vector-valued) unctions in L 2 (Ω;V) wose divergencebelongstol 2 (Ω). Similarly, H(div,Ω;M)denotestesubspaceo(matrix-valued) unctions in L 2 (Ω;M) wose divergence (by rows) belongs to L 2 (Ω;V). Assuming tat X is an inner product space, ten L 2 (Ω;X) as a natural norm and inner product, wic will be denoted by and (, ), respectively. For a Sobolev space H s (Ω;X), we denote te norm by s and or H(div,Ω;X), te norm is denoted by v 3

6 6 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER v div := ( v 2 + divv 2 ) 1/2. Te space P k (Ω) denotes te space o polynomial unctions on Ω o total degree k. Usually we abbreviate tis to just P k. In tis paper we sall consider mixed inite element approximations derived rom (1.5). Tese scemes take te orm: Find (σ,u,p ) Σ V Q suc tat (Aσ Ω : τ +divτ u +τ : p )dx = 0, τ Σ, (2.1) divσ Ω vdx = vdx, v V Ω, σ Ω : qdx = 0, q Q, were now Σ H(div,Ω;M), V L 2 (Ω;V) and Q L 2 (Ω;K). Following te general teory o mixed inite element metods, c. [12, 13], te stability o te saddle point system (2.1) is ensured by te ollowing conditions: (A1) τ 2 div c 1(Aτ,τ) wenever τ Σ satisies (divτ,v) = 0 v V, and (τ,q) = 0 q Q, (A2) or all nonzero (v,q) V Q, tere exists nonzero τ Σ wit (divτ,v) + (τ,q) c 2 τ div ( v + q ), were c 1 and c 2 are positive constants independent o. Te main result o tis paper, given in Teorem 7.1, is to construct a new amily o stable inite element spaces Σ, V, Q tat satisy te stability conditions (A1) and (A2). We sall sow tat or r 0, te coices o te Nédélec second amily o H(div) elements o degree r +1 or Σ, c. [26], and o discontinuous piecewise polynomials o degree r or V and Q provide a stable inite element approximation. In contrast to te previous work described in te introduction, no stabilizing bubble unctions are needed; nor is interelement continuity imposed on te multiplier. In 8 we also discuss a somewat simpler lowest order element (r = 0) in wic te local stress space is a strict subspace o te ull space o linear matrix ields. Our approac to te construction o stable mixed elements or elasticity is motivated by te success in developing stable mixed elements or steady eat conduction (i.e., te Poisson problem) based on discretizations o te de Ram complex. We recall (see, e.g., [7]) tat tere is a close connection between te construction o suc elements and discretizations o te de Ram complex (2.2) R C (Ω) grad C (Ω;V) curl C (Ω;V) div C (Ω) 0. More speciically, a key to te construction and analysis o stable mixed elements is a commuting diagram o te orm (2.3) grad curl R C (Ω) C (Ω;V) C (Ω;V) Π 1 R W Πc Π d div C (Ω) 0 Π 0 grad U curl V div Q 0. Here, te spaces V H(div) and Q L 2 are te inite element spaces used to discretize te lux and temperature ields, respectively. Te spaces U H(curl) and W H 1 are additional inite element spaces, wic can be ound or all well-known stable element coices. Te bottom row o te diagram is a discrete de Ram complex, wic is exact wen

7 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 7 te de Ram complex is (i.e., wen te domain is contractible). Te vertical operators are projections determined by te natural degrees o reedom o te inite element spaces. As pointed out in [7], tere are many suc discretizations o te de Ram complex. A diagram analogous to (2.3), but wit te de Ram complex replaced by te elasticity complex deined just below, will be crucial to our construction o stable mixed elements or elasticity. Discretization o te elasticity complex also gives insigt into te diiculties o constructing inite element approximations o te mixed ormulation o elasticity wit strongly imposed symmetry, c. [8]. 3. Te elasticity complex We now proceed to a description o two elasticity complexes, corresponding to strongly or weakly imposed symmetry o te stress tensor. In te case o strongly imposed symmetry, relevant to te mixed elasticity system (1.3), te caracterization o te divergence-ree symmetric matrix ields will be needed. In order to give suc a caracterization, deine curl : C (Ω;M) C (Ω;M) to be te dierential operator deined by taking curl o eac rowotematrix. Ten deineasecondorder dierential operatorj : C (Ω;S) C (Ω;S) by (3.1) Jτ = curl(curlτ) T, τ C (Ω;S). It is easy to ceck tat div J = 0 and tat J ǫ = 0. In oter words, (3.2) T C (V) ǫ C (S) J C (S) div C (V) 0, is a complex. Here te dependence o te domain Ω is suppressed, i.e., C (S) = C (Ω;S), and T = T(Ω) denotes te six dimensional space o ininitesimal rigid motions on Ω, i.e., unctions o te orm x a+bx wit a V and B K. In act, wen Ω is contractible, ten (3.2) is an exact sequence, a act wic will ollow rom te discussion below. Te complex (3.2) will be reerred to as te elasticity complex. A natural approac to te construction o stable mixed inite elements or elasticity would be to extend te complex (3.2) to a complete commuting diagram o te orm (2.3), were (3.2) iste toprowandte bottomrowis adiscrete analogue. However, due to tepointwise symmetry requirement on te discrete stresses, tis construction requires piecewise polynomials o ig order. For te corresponding problem in two space dimensions, suc a complex was proposed in [10] wit a piecewise cubic stress space, c. also [8]. An analogous complex was derived in te tree dimensional case in [3]. It uses a piecewise quartic space, wit 162 degrees o reedom on eac tetraedron or te stresses. We consider te ormulation based on weakly imposed symmetry o te stress tensor, i.e., te mixed system (1.5). Ten te relevant complex is, instead o (3.2), (3.3) T C (V K) Here, (grad, I) C (M) J C (M) T = {(v,gradv) v T}, (div,skw) T C (V K) 0. and J : C (Ω;M) C (Ω;M) denotes te extension o te operator deined on C (Ω;S) by (3.1) suc tat J 0 on C (Ω;K). We remark tat J may be written (3.4) Jτ = curlξ 1 curlτ,

8 8 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER were Ξ : M M is te algebraic operator (3.5) Ξµ = µ T tr(µ)δ, Ξ 1 µ = µ T 1 2 tr(µ)δ, wit δ te identity matrix. Indeed, i τ is symmetric, ten curlτ is trace ree, and tereore te deinition (3.4) reduces to (3.1) on C (Ω;S). On te oter and, i τ is skew wit axial vector u, ten curlτ = Ξgradu, and so curlξ 1 curlτ = 0. Observe tat tere is a close connection between (3.2) and (3.3). In act, (3.2) can be derived rom (3.3) by perorming a projection step. To see tis, consider te diagram (3.6) T C (V K) π 0 (grad, I) C (M) π 1 J C (M) π 2 (div,skw) T C (V K) 0 π 3 T C (V) ǫ C (S) were te projection operators π k are deined by J C (S) div C (V) 0, π 0 (u,q) = u, π 1 (σ) = π 2 (σ) = sym(σ), π 3 (u,q) = u divq. We may identiy C (V) wit a subspace o C (V K), namely, {(u,q) : u C (V),q = skw(gradu)}. Under tis identiication, T C (V) corresponds to T C (V K). We identiy te C (V) on te rigt wit a dierent subspace o C (V K), namely, {(u,q) : u C (V),q = 0}. Wit tese identiications, te bottom row is a subcomplex o te top row, and te operators π k are all projections. Furtermore, te diagram commutes. It ollows easily tat te exactness o te upper row implies exactness o te bottom row. In te next section, we sall discuss tese complexes urter. In particular, we sow te elasticity complex wit weakly imposed symmetry, i.e., (3.3), ollows rom te de Ram complex (2.2). Hence, as a consequence o te discussion above, bot (3.2) and (3.3) will ollow rom (2.2). 4. From te de Ram to te elasticity complex Te main purpose o tis section is to demonstrate te connection between te de Ram complex (2.2) and te elasticity complex (3.3). In particular, we sow tat wenever (2.2) is exact, (3.3) is exact. Tis section serves as an introduction to a corresponding construction o a discrete elasticity complex, to be given in te next section. In te ollowing section, te discrete complex will be used to construct stable inite elements or te system (1.5). Te de Ram complex (2.2) is most clearly stated in terms o dierential orms. Here we briely recall te deinitions and properties we will need. We use a completely coordinate-ree approac. For a sligtly more expanded discussion and te expressions in coordinates see, e.g., [7, 4]. We let Λ k denote te space o smoot dierential k-orms on Ω, i.e. Λ k = Λ k (Ω) = C (Ω;Alt k V), were Alt k V denotes te vector space o alternating k-linear maps on V. I ω Λ k we let ω x Alt k V denote ω evaluated at x, i.e., we use subscripts to indicate te spatial dependence.

9 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 9 Using te inner product on Alt k V inerited rom te inner product on V (see equation (4.1) o [7, 4]), we may also deine te Hilbert space L 2 Λ k (Ω) = L 2 (Ω;Alt k V) o square integrable dierential orms wit norm denoted by, and also te mt order Sobolev space H m Λ k (Ω) = H m (Ω;Alt k V), consisting o square integrable k-orms or wic te norm ω m := ( α ω 2) 1/2 α m is inite (were te sum is over multi-indices o degree at most m). Tus, 0-orms are scalar unctions and 1-orms are covector ields. We will not empasize te distinction between vectors and covectors, since, given te inner product in V, we may identiy a 1-orm ω wit te vector ield v or wic ω(p) = v p, p V. In te 3-dimensional case, we can identiy a 2-orm ω wit a vector ield v and a 3-orm µ wit a scalar ield c by ω(p,q) = v p q, µ(p,q,r) = c(p q r), p,q,r V. Te exterior derivative d = d k : Λ k Λ k+1 is deined by k+1 (4.1) dω x (v 1,...,v k+1 ) = ( 1) j+1 vj ω x (v 1,...,ˆv j,...,v k+1 ), ω Λ k,v 1,...,v k+1 V, j=1 were te at is used to indicate a suppressed argument and v denotes te directional derivative in te direction o te vector v. It is useul to deine HΛ k = {ω L 2 (Ω;Alt k V) dω L 2 (Ω;Alt k+1 V)}, wit norm given by ω 2 HΛ = ω 2 + dω 2. Using te identiications given above, te d k correspond to grad, curl, and div or k = 0,1,2, respectively, and te HΛ k correspond to H 1, H(curl), H(div), and, or k = 3, L 2. Te de Ram complex (2.2) can ten be written (4.2) R Λ 0 d Λ 1 d Λ 2 d Λ 3 0. It is a complex since d d = 0. A dierential k-orm ω on Ω may be restricted to a dierential k-orm on any submaniold M Ω: at eac point o M te restriction o ω is an alternating linear orm on tangent vectors. Moreover i dimm = k, te integral ω is deined. M I X is a vector space, ten Λ k (X) = Λ k (Ω;X) reers to te k-orms wit values in X, i.e., Λ k (X) = C (Ω;Alt k (V;X)), were Alt k (V;X) are alternating k-linear orms on V wit values in X. Given an inner product on X, we obtain an inner product on Λ k (X). Obviously te corresponding complex (4.3) X Λ 0 (X) d Λ 1 (X) d Λ 2 (X) d Λ 3 (X) 0, is exact wenever te de Ram complex is. We now construct te elasticity complex as a subcomplex o a complex isomorpic to te de Ram complex wit values in te six-dimensional vector space W := K V. First, or any x R 3 we deine K x : V K by K x v = 2skw(xv T ). We ten deine an operator K : Λ k (Ω;V) Λ k (Ω;K) by (4.4) (Kω) x (v 1,...,v k ) = K x [ω x (v 1,...,v k )].

10 10 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Next, we deine an isomorpism Φ : Λ k (W) Λ k (W) by wit inverse given by Φ(ω,µ) = (ω +Kµ,µ), Φ 1 (ω,µ) = (ω Kµ,µ). Next, deineteoperatora : Λ k (W) Λ k+1 (W)byA = ΦdΦ 1. Insertingteisomorpisms Φ in te W-valued de Ram sequence, we obtain a complex (4.5) Φ(W) Λ 0 (W) A Λ 1 (W) A Λ 2 (W) A Λ 3 (W) 0, wic is exact wenever te de Ram complex is. TeoperatorAasasimpleorm. UsingtedeinitionoΦ, weobtainor(ω,µ) Λ k (W), A(ω,µ) = Φ d(ω Kµ,µ) = Φ(dω dkµ,dµ) = (dω Sµ,dµ), were S = S k : Λ k (V) Λ k+1 (K), k = 0,1,2, is given by S = dk Kd. Using te deinition (4.1) o te exterior derivative, te deinition (4.4) o K, and te Leibniz rule (4.6) d(ω µ) = dω µ+( 1) k ω dµ, ω Λ k, µ Λ l, we obtain k+1 (4.7) (Sω)(v 1,...,v k+1 ) = ( 1) j+1 K vj [ω(v 1,...,ˆv j,...v k+1 )], ω Λ k (Ω;V). j=1 Note tat te operator S is purely algebraic, and independent o x. Since d 2 = 0 we ave ds = d 2 K dkd = (dk Kd)d or (4.8) ds = Sd. Noting tat (S 1 µ)(v 1,v 2 ) = K v1 [µ(v 2 )] K v2 [µ(v 1 )] = 2skw[v 1 µ(v 2 ) T v 2 µ(v 1 ) T ], we ind, using te identity (4.9) a b = 2vecskwab T, tat S 1 is invertible wit µ Λ 1 (Ω;V), v 1,v 2 V, (S 1 1 ω)(v 1) v 2 v 3 = 1 2 [vec( ω(v 2,v 3 ) ) v 1 vec ( ω(v 1,v 2 ) ) v 3 +vec ( ω(v 1,v 3 ) ) v 2 ], We now deine te desired subcomplex. Deine ω Λ 2 (Ω;K),v 1,v 2,v 3 V. Γ 1 = {(ω,µ) Λ 1 (Ω;W) dω = S 1 µ}, Γ 2 = {(ω,µ) Λ 2 (Ω;W) ω = 0}, wit projections π 1 : Λ 1 (Ω;W) Γ 1 and π 2 : Λ 2 (Ω;W) Γ 2 given by π 1 (ω,µ) = (ω,s 1 1 dω), π2 (ω,µ) = (0,µ+dS 1 1 ω).

11 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 11 Using (4.8), it is straigtorward to ceck tat A maps Λ 0 (W) into Γ 1 and Γ 1 into Γ 2, and tat te diagram (4.10) Φ(W) Λ 0 (W) Λ 1 (W) id Φ(W) Λ 0 (W) A A Λ 2 (W) π 1 π 2 A Λ 3 (W) 0 id A Γ 1 A Γ 2 A Λ 3 (W) 0 commutes, and tereore te subcomplex in te bottom row is exact wen te de Ram complex is. Tis subcomplex is, essentially, te elasticity complex. Indeed, by identiying elements (ω,µ) Γ 1 wit ω Λ 1 (K), and elements (0,µ) Γ 2 wit µ Λ 2 (V), te subcomplex becomes (4.11) Φ(W) Λ 0 (K V) (d 0, S 0 ) Λ 1 (K) d 1 S 1 1 d 1 Λ 2 (V) ( S 2,d 2 ) T Λ 3 (K V) 0. Tis complex may be identiied wit (3.3). As an initial step o tis identiication we observe tat te algebraic operator Ξ : C (M) C (M) appearing in (3.3) via (3.4) and te operator S 1 : Λ 1 (V) Λ 2 (K) are connected by te identity (4.12) Ξ = Υ 1 2 S 1Υ 1, were Υ 1 : C (M) Λ 1 (V) and Υ 2 : C (M) Λ 2 (K) are given by Υ 1 F(v) = Fv and Υ 2 F(v 1,v 2 ) = vec 1 F(v 1 v 2 ) or F C (M). In act, using (4.9), we ave or any v 1,v 2 V On te oter and, S 1 Υ 1 F(v 1,v 2 ) = 2skw[v 1 (Fv 2 ) T v 2 (Fv 1 ) T ] = vec 1 (v 2 Fv 1 v 1 Fv 2 ). Υ 2 ΞF(v 1,v 2 ) = vec 1 [ΞF(v 1 v 2 )], and ence (4.12) ollows rom te algebraic identity ΞF(v 1 v 2 ) = v 2 Fv 1 v 1 Fv 2, wic olds or any F M. We may urter identiy te our spaces o ields in (3.3) wit te corresponding spaces o orms in (4.11) in a natural way: (u,p) C (V K) (vec 1 u,vecp) Λ 0 (K V). F C (M) ω Λ 1 (K) given by ω(v) = vec 1 (Fv). F C (M) µ Λ 2 (V) given by µ(v 1,v 2 ) = F(v 1 v 2 ). (u,p) C (V K) (ω,µ) Λ 3 (K V) given by ω(v 1,v 2,v 3 ) = p(v 1 v 2 v 3 ), µ(v 1,v 2,v 3 ) = u(v 1 v 2 v 3 ). Under tese identiications te we ind tat d 0 : Λ 0 (K) Λ 1 (K) corresponds to te row-wise gradient C (V) C (M). S 0 : Λ 0 (V) Λ 1 (K) corresponds to te inclusion o C (K) C (M). d 1 S1 1 d 1 : Λ 1 (K) Λ 2 (V) corresponds to J = curlξ 1 curl : C (M) C (M). d 2 : Λ 2 (V) Λ 3 (V) corresponds to te row-wise divergence C (M) C (V). S 2 : Λ 2 (V) Λ 3 (K) corresponds to te operator 2skw : C (M) C (K).

12 12 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Tus, modulo tese identiications and te (unimportant) constant actor in te last identiication, (3.3) and (4.11) are identical. Hence we ave establised te ollowing result. Teorem 4.1. Wen te de Ram complex (2.2) is exact, (i.e., te domain is contractible), ten so is te elasticity complex (3.3). To end tis section, we return to te operator S : Λ k (V) Λ k+1 (K) deined by S = dk Kd. Let K : Λ k (K) Λ k (V) be te adjoint o K (wit respect to te Euclidean inner product on V and te Frobenius inner product on K), wic is given by (K ω) x (v 1,...,v k ) = 2ω x (v 1,...,v k )x. Deine S : Λ k (K) Λ k+1 (V) by S = dk K d. Recall tat te wedge product : Λ k Λ l Λ k+l is given by (ω µ)(v 1,...,v k+l ) = (signσ)ω(v σ1,...,v σk )µ(v σk+1,...,v σk+l ), ω Λ k,µ Λ l,v i V, were te sum is over te set o all permutations o {1,...,k+l}, or wic σ 1 < σ 2 < < σ k and σ k+1 < σ j+2 < < σ k+l. Tis extends as well to dierential orms wit values in an inner product space, using te inner product to multiply te terms inside te summation. Using te Leibniz rule (4.6), we ave (4.13) (Sω) µ = ( 1) k ω S µ, ω Λ k (V), µ Λ l (K). We tus ave dkω µ = ( 1) k+1 Kω dµ+d(kω µ) = ( 1) k+1 ω K dµ+d(ω K µ), and Kdω µ = dω K µ = ( 1) k+1 ω dk µ+d(ω K µ), Subtracting tese two expressions gives (4.13). For later reerence, we note tat, analogously to (4.7), we ave k+1 (4.14) (S ω)(v 1,...,v k+1 ) = 2 ( 1) j+1 ω(v 1,...,ˆv j,...v k+1 )v j, ω Λ k (Ω;K). j=1 5. Te discrete construction In tis section we derive a discrete version o te elasticity sequence by adapting te construction o te previous section. To carry out te construction, we will use two discretizations o te de Ram sequence. For k = 0,1,2,3, let Λ k denote a inite-dimensional spaceohλ k orwicdλ k Λk+1, andorwictereexistprojectionsπ = Π k : Λk Λ k wic make te ollowing diagram commute: (5.1) R Λ 0 d Λ 1 d Λ 2 d Λ 3 0 Π Π Π Π R Λ 0 d Λ 1 d Λ 2 d Λ 3 0 Tis is simply te diagram (2.3) written in te language o dierential orms. We do not make a speciic coice o te discretization yet, but, as recalled in 2, tere exist many suc discrete de Ram complexes based on piecewise polynomials. In act, as explained in [7], or eac polynomial degree r 0 we may coose Λ 3 to be te space o all piecewise polynomial 3-orms wit respect to some simplicial decomposition o Ω, and construct our

13 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 13 suc diagrams. We make te assumption tat P 1 (Ω) Λ 0, wic is true in all te cases mentioned. Let Λ k beasecond set oinitedimensional spaces witcorresponding projectionoperators Π enjoying te same properties, giving us a second discretization o te de Ram sequence. Supposing a compatibility condition between tese two discretizations, wic we describe below, we sall construct a discrete elasticity complex. We start wit te complex (5.2) K V Λ 0 (K) Λ 0 (V) d d Λ 3 (K) Λ 3 (V) 0 were Λ k (K) denotes te K-valued analogue o Λk and similarly or Λ k (V). For brevity, we enceort write Λ k (W) or Λk (K) Λ k (V). As a discrete analogue o te operator K, we deine K : Λk (V) Λ k (K) by K = Π K were Π is te interpolation operator onto Λ k (K). Next deine S = S k, : Λ k (V) Λk+1 (K) by S = dk K d, or k = 0,1,2. Observe tat te discrete version o (4.8), (5.3) ds = S d, ollows exactly as in te continuous case. From te commutative diagram (5.1), we see tat S = dπ K Π Kd = Π (dk Kd) = Π S. Continuing to mimic te continuous case, we deine te automorpism Φ on Λ k (W) by Φ (ω,µ) = (ω +K µ,µ), and te operator A : Λ k (W) Λk+1 (W) by A = Φ dφ 1, wic leads to A (ω,µ) = (dω S µ,dµ). Tanks to te assumption tat P 1 Λ 0, we ave Φ (W) = Φ(W). Hence, inserting te isomorpisms Φ into (5.2), we obtain (5.4) Φ(W) Λ 0 (W) A Λ 1 (W) A Λ 2 (W) A Λ 3 (W) 0. In analogy to te continuous case, we deine Γ 1 = {(ω,µ) Λ1 (W) dω = S 1,µ}, Γ 2 = {(ω,µ) Λ2 (W) ω = 0}. As in te continuous case, we can identiy Γ 2 wit Λ 2 (V), but, unlike in te continuous case, we cannot identiy Γ 1 wit Λ1 (K), since we do not require tat S 1, be invertible (it is not in te applications). Hence, in order to derive te analogue o te diagram (4.10) we require a surjectivity assumption: (5.5) Te operator S 1, maps Λ 1 (V) onto Λ 2 (K). Under tis assumption, te operator S = S 1, as a rigt inverse S mapping Λ2 (K) into Λ 1 (V). Tis allows us to deine discrete counterparts o te projection operators π1 and π 2 by π 1 (ω,µ) = (ω,µ S S µ+s dω), π2 (ω,µ) = (0,µ+dS ω),

14 14 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER and obtain te discrete analogue o (4.10): (5.6) Φ(W) Λ 0 (W) A Λ 1 (W) A Λ 2 (W) A Λ 3 (W) 0 id π 1 π 2 id Φ(W) Λ 0 (W) A Γ 1 A Γ 2 A Λ 3 (W) 0 It is straigtorward to ceck tat tis diagram commutes. For example, i (ω,µ) Λ 0 (W), ten π 1 A (ω,µ) = π 1 (dω S µ,dµ) = (dω S µ,dµ S S dµ+s d[dω S µ]) = (dω S µ,dµ S [S dµ+ds µ]) = A (ω,µ), were te last equality ollows rom (5.3). Tus te bottom row o (5.6) is a subcomplex o te top row, and te vertical maps are commuting projections. In particular, wen te top row is exact, so is te bottom. Tus we ave establised te ollowing result. Teorem 5.1. For k = 0,...,3, let Λ k be a inite dimensional subspace o HΛk or wic dλ k Λk+1 and or wic tere exist projections Π = Π k : Λk Λ k tat make te diagram (5.1) commute. Let Λ k be a second set o inite dimensional spaces wit corresponding projection operators Π k wit te same properties. I S 1, := dπ 1 K Π2 Kd maps Λ 1 (V) onto Λ 2 (K), and te bottom row o (5.1) is exact or bot sequences Λk and Λ k, ten te discrete elasticity sequence given by te bottom row o (5.6) is also exact. Te exactness o te bottom row o (5.6) suggests tat te ollowing coice o inite element spaces will lead to a stable discretization o (2.1): Σ Λ 2 (V), V Λ 3 (V), Q Λ 3 (K). In te next section we will make speciic coices or te discrete de Ram complexes, and ten veriy te stability in te ollowing section. For use in te next section, we state te ollowing result, giving a suicient condition or te key requirement tat S 1, be surjective. Proposition 5.2. I te diagram (5.7) Λ 1 (V) Π 1 S 1 Λ 2 (K) Π 2 commutes, ten S 1, is surjective. Λ 1 (V) S 1, Λ 2 (K) 6. A amily o discrete elasticity complexes In tis section, we present a amily o examples o te general discrete construction presented in te previous section, by coosing speciic discrete de Ram complexes. Tese urnis a amily o discrete elasticity complexes, indexed by an integer degree r 0. In te next section we use tese complexes to derive inite elements or elasticity. In te lowest

15 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 15 order case, te metod will require only piecewise linear unctions to approximate stresses and piecewise constants to approximate te displacements and multipliers. We begin by recalling te two principal amilies o piecewise polynomial spaces o dierential orms, ollowing te presentation in [7]. We enceort assume tat te domain Ω is a contractible polyedron. Let T be a triangulation o Ω. Let T be a triangulation o Ω by tetraedra, and set P r Λ k (T ) = {ω HΛ k (Ω) ω T P r Λ k (T) T T }, P + r Λk (T ) = {ω HΛ k (Ω) ω T P + r Λk (T) T T }. Here P + r Λk (T) := P r Λ k (T)+κP r Λ k+1 (T) were κ : Λ k+1 (T) Λ k (T) is te Koszul dierential deined by (κω) x (v 1,,v k ) = ω x (x,v 1,,v k ). Te spaces P r + Λ 0 (T ) = P r+1 Λ 0 (T ) correspond to te usual degree r+1 Lagrange piecewise polynomial subspaces o H 1, and te spaces P r + Λ 3 (T ) = P r Λ 3 (T ) correspond to te usual degree r subspace o discontinuous piecewise polynomials in L 2 (Ω). For k = 1 and 2, te spaces P r + Λk (T ) correspond to te discretizations o H(curl) and H(div), respectively, presented by Nédélec in [25], and te spaces P r Λ k (T ) are te ones presented by Nédélec in [26]. An element ω P r Λ k (T ) is uniquely determined by te ollowing quantities: (6.1) ω ζ, ζ P + r d 1+k Λd k (), d (T ), k d 3. Here d (T ) is te set o vertices, edges, aces, or tetraedra in te mes T, or d = 0,1,2,3, respectively, and or r < 0, we interpret P r + Λ k (T) = P r Λ k (T) = 0. Note tat or ω Λ k, ω naturally restricts on te ace to an element o Λ k (). Tereore, or ζ Λ d k (), te wedge product ω ζ belongs to Λ d () and ence te integral o ω ζ on te d-dimensional ace o T is a well-deined and natural quantity. Using te quantities in (6.1), we obtain a projection operator rom Λ k to P r Λ k (T ). Similarly, an element ω P r + Λ k (T ) is uniquely determined by (6.2) ω ζ, ζ P r d+k Λ d k (), d (T ), k d 3, and so tese quantities determine a projection. I X is a vector space, we use te notation P r Λ k (T ;X) and P r + Λk (T ;X) to denote te corresponding spaces o piecewise polynomial dierential orms wit values in X. Furtermore, i X is an inner product space, te corresponding degrees o reedom are given by (6.1) and (6.2), but were te test spaces are replaced by te corresponding X valued spaces. To carry out te construction described in te previous section we need to coose te two sets o spaces Λ k and Λ k or k = 0,1,2,3. We ix r 0 and set Λk = P+ r Λk (T ), k = 0,1,2,3, and Λ 0 = P r+2λ 0 (T ), Λ 1 = P+ r+1λ 1 (T ), Λ 2 = P r+1λ 2 (T ), and Λ 3 = P rλ 3. As explained in [7], bot tese coices give a discrete de Ram sequence wit commuting projections, i.e., a diagram like (5.1) makes sense and is commutative. We establis te key surjectivity assumption or our coice o spaces, by veriying te commutativity o (5.7).

16 16 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Lemma 6.1. Let Λ 1 (V) = P+ r+1λ 1 (T ;V) and Λ 2 (K) = P+ r Λ 2 (T ;K) wit projections Π 1, Π 2 deined via te corresponding vector valued moments o te orm (6.1) and (6.2). I S 1, = Π 2 S 1 ten (6.3) S 1, Π1 = Π 2 S 1, and so S 1, is surjective. Proo. We must sow tat (Π 2 S 1 S 1, Π1 )σ = 0 or σ Λ 1 (V). Deining ω = (I Π 1 )σ, te required condition becomes Π 2 S 1ω = 0. Since Π 1 ω = 0, we ave (6.4) ω ζ = 0, ζ P r+2 d Λ d 1 (;V), d (T ), 2 d 3, (in act (6.4) olds or d = 1 as well, but tis is not used ere). We must sow tat (6.4) implies (6.5) S 1 ω µ = 0, µ P r+2 d Λ d 2 (;K), d (T ), 2 d 3. From (4.13), we ave S 1 ω µ = ω ζ, were ζ = S d 2 µ P r+2 dλ d 1 (;V) or µ P r+2 d Λ d 2 (;K), as is evident rom (4.14). Hence (6.5) ollows rom (6.4). 7. Stable mixed inite elements or elasticity Based on te discrete elasticity complexes just constructed, we obtain mixed inite element spaces or te ormulation (2.1) o te elasticity equations by coosing Σ, V, and Q as te spaces o matrix and vector ields corresponding to appropriate spaces o orms in te K- and V-valued de Ram sequences used in te construction. Speciically, tese are (7.1) Σ P r+1 Λ 2 (T ;V), V P r Λ 3 (T ;V), Q P r Λ 3 (T ;K). In oter terminology, Σ may be tougt o as te product o tree copies o te Nédélec H(div) space o te second kind o degree r +1, and V and Q are spaces o all piecewise polynomials o degree at most r wit values in K and V, respectively, wit no imposed interelement continuity. In tis section, we establis stability and convergence or tis inite element metod. Te stability o te metod requires te two conditions (A1) and (A2) stated in 2. Te irst condition is obvious since, by construction, divσ V, i.e., dp r+1 Λ 2 (T ;V) P r Λ 3 (T ;V). Te condition (A2) is more subtle. We will prove a stronger version, namely (A2 ) or all nonzero (v,q) V Q, tere exists nonzero τ Σ wit divτ = v, 2Π Q skwτ = q and τ div c( v + q ), were Π Q is te L 2 projection into Q and c is a constant. Recalling tat Γ 2 = 0 P r+1λ 2 (T ;V) and A (0,σ) = ( S 2, σ,dσ), and tat te operator S 2 corresponds on te matrix level to 2skw, we restate condition (A2 ) in te language o dierential orms.

17 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 17 Teorem 7.1. Given (ω,µ) P r Λ 3 (T ;K) P r Λ 3 (T ;V), tere exists σ P r+1 Λ 2 (T ;V) suc tat A (0,σ) = (ω,µ) and (7.2) σ HΛ c( ω + µ ), were te constant c is independent o ω,µ and. Beore proceeding to te proo, we need to establis some bounds on projection operators. We do tis or te corresponding scalar-valued spaces. Te extensions to vector-valued spaces are straigtorward. First we claim tat (7.3) Π 2 η c η 1 η H 1 Λ 2, Π 3 ω c ω 0 ω H 1 Λ 3. Here te constant may depend on te sape regularity o te mes, but not on te messize. Te second bound is obvious (wit c = 1), since Π 3 is just te L2 projection. Te irst bound ollows by a standard scaling argument. Namely, let ˆT denote te reerence simplex. For any ˆβ P r+1 Λ 2 (ˆT), we ave (7.4) ˆβ 0,ˆT c( ˆ ˆµ ˆ ˆβ ˆµ + ˆζ ˆT ˆβ ˆζ ), were ˆ ranges over te aces o ˆT, ˆµ over a basis or P + r Λ 0 (ˆ), and ˆζ over a basis or P r 1Λ + 1 (ˆT). Tis is true because te integrals on te rigt and side o (7.4) orm a set o degrees o reedom or ˆβ P r+1 Λ 2 (ˆT) (see (6.1)), and so we may use te equivalence o all norms on tis inite dimensional space. We apply tis result wit ˆβ = ˆΠ 2 ˆη, were ˆΠ 2 is te projection deined to preserve te integrals on te rigt and side o (7.4). It ollows tat ˆΠ 2 ˆη 0,ˆT c( ˆ ˆη ˆµ + ˆζ ˆη ˆζ ) c ˆη 1,ˆT, ˆµ ˆ ˆT were we ave used a standard trace inequality in te last step. Next, i T is an arbitrary simplex and η H 1 Λ 2 (T), we map te reerence simplex ˆT onto T by an aine map ˆx Bˆx+b, and deine ˆη H 1 Λ 2 (ˆT) by ˆηˆx (ˆv 1,ˆv 2 ) = η x (Bˆv 1,Bˆv 2 ), or any x = Bˆx+b T and any vectors ˆv 1,ˆv 2. It is easy to ceck tat Π2 η = ˆΠ 2 ˆη, and tat Π 2 η 0,T c ˆΠ 2 ˆη 0,ˆT c ˆη 1,ˆT c( η 0,T + η 1,T ) c η 1,T. Squaring and adding tis over all te simplices in te mes T gives te irst bound in (7.3). WealsoneedaboundonteprojectionoaorminH 1 Λ 1 into Λ 1 = P+ r+1λ 1 (T ). However, te projection operator Π 1 is not bounded on H1, because its deinition involves integrals over edges. A similar problem as arisen beore (see, e.g., [10]), and we use te same remedy.

18 18 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER Namely we start by deining an operator Π 1 0 : H1 Λ 1 P r+1λ + 1 (T ) by te conditions (7.5) Π 1 0ω ζ = ω ζ, ζ P r 1 Λ 2 (T), T T, T T (7.6) Π 1 0ω ζ = ω ζ, ζ P r Λ 1 (), 2 (T ), (7.7) Π 1 0 ω ζ = 0, ζ P r+1λ 0 (e), e 1 (T ). e Note tat, in contrast to Π 1, in te deinition o Π 1 0, we ave set te troublesome edge degrees o reedom to zero. Let ˆΠ 1 0 : H1 Λ 1 (ˆT) P r+1 + Λ1 (ˆT) be deined analogously on te reerence element. Now or ˆρ H 1 Λ 1 (ˆT), ˆdˆΠ 1 0ˆρ P r+1λ 2 (ˆT), so ˆdˆΠ 1 0ˆρ 0,ˆT c( ˆ ˆdˆΠ 1 0ˆρ ˆµ + ˆζ ˆdˆΠ 1 0ˆρ ˆζ ), ˆµ ˆ ˆT were again ˆ ranges over te aces o ˆT, µ over a basis o P r + Λ0 (ˆ), and ζ over a basis o P r 1Λ + 1 (ˆT). But ˆdˆΠ 1 0ˆρ ˆµ = ˆΠ 1 0ˆρ ˆdˆµ = ˆρ ˆdˆµ, ˆ ˆ were we ave used Stokes teorem and te act tat te vanising o te edge quantities in te deinition o ˆΠ 1 0 to obtain te irst equality, and te ace degrees o reedom entering te deinition o ˆΠ 1 0 to obtain te second. Similarly, ˆdˆΠ 1 0ˆρ ˆζ = ˆΠ 1 0ˆρ ˆdˆζ + ˆΠ 1 0ˆρ ˆζ = ˆρ ˆdˆζ + ˆρ ˆζ = ˆdˆρ ˆζ. It ollows tat ˆT ˆT ˆT ˆdΠ 1 0ˆρ 0,ˆT c ˆρ 1,ˆT, ρ H 1 Λ 1 (ˆT). Wen we scale tis result to an arbitrary simplex and add over te mes, we obtain dπ 1 0 ρ c( 1 ρ + ρ 1 ), ρ H 1 Λ 1 (Ω). To remove te problematic 1 in te last estimate, we introduce te Clement interpolant R mapping H 1 Λ 1 into continuous piecewise linear 1-orms (still ollowing [10]). Ten Deining Π 1 : H1 Λ 1 P + r+1λ 1 by ρ R ρ + ρ R ρ 1 c ρ 1, ρ H 1 Λ 1. (7.8) Π1 = Π 1 0(I R )+R, we obtain d Π 1 ρ d Π 1 0(I R )ρ + dr ρ c( 1 (I R )ρ + (I R )ρ 1 + dr ρ ) c ρ 1. Tus we ave sown tat (7.9) d Π 1 ρ c ρ 1, ρ H 1 Λ 1. ˆT ˆ ˆT ˆT

19 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 19 Having modiied Π 1 to obtain te bounded operator Π 1, we now veriy tat te key property (6.3) in Lemma 6.1 carries over to (7.10) S 1, Π1 = Π 2 S 1, were we now use te vector-valued orms o te projection operators. It ollows easily rom (7.8), (7.5), and (7.6) tat (6.4) olds wit ω = (I Π 1 )σ, so tat te proo o (7.10) is te same as or (6.3). We can now give te proo o Teorem 7.1. Proo o Teorem 7.1. Given µ P r Λ 3 (T ;V) tere exists η H 1 Λ 2 (V) suc tat dη = µ wit te bound η 1 c µ (since d maps H 1 Λ 2 onto L 2 Λ 3 ). Similarly, given ω P r Λ 3 (T ;K) tere exists τ H 1 Λ 2 (K) wit dτ = ω + S 2, Π2 η wit te bound τ 1 c ω +S 2, Π2 η. Let ρ = S 1 1 τ (recall tat S 1 is an isomorpism) and set σ = d Π 1 ρ+ Π 2 η. We will now sow tat A (0,σ) = (ω,µ). From te deinition o σ, we ave S 2, σ = S 2, d Π 1 ρ S 2, Π2 η. Ten, using (5.3), (7.10), and te commutativity dπ = Π d, we see S 2, d Π 1 ρ = ds 1, Π 1 ρ = dπ2 S 1ρ = dπ 2 τ = Π 3 dτ = Π 3 (ω +S 2, Π2 η) = ω S 2, Π2 η. Combining, we get S 2, σ = ω as desired. Furter, rom te commutativity d Π = Π d and te deinition o η, we get dσ = d Π 2 η = Π 3 dη = Π 3 µ = µ, and so we ave establised tat A (0,σ) = (µ,ω). It remains to prove te bound (7.2). Using (7.3), we ave S 2, Π2 η = Π 3 S 2 Π 2 η c S 2 Π 2 η c Π 2 η c η 1 c µ. Tus ρ 1 c τ 1 c( ω + µ ). Using (7.9), wetenget d Π 1 ρ c ρ 1 c( ω + µ ), and, using (7.3), tat Π 2 η c η 1 c µ. Tereore σ c( ω + µ ), wile dσ = µ, and tus we ave te desired bound (7.2). We ave tus veriied te stability conditions (A1) and (A2), and so may apply te standard teory o mixed metods (c. [12], [13], [17], [20]) and standard results about approximation by inite element spaces to obtain convergence and error estimates. Teorem 7.2. Suppose (σ,u,p) is te solution o te elasticity system (1.5) and (σ,u,p ) is te solution o discrete system (2.1), were te inite element spaces Σ, V, and Q are given by (7.1) or some integer r 0. Ten tere is a constant C, independent o, suc

20 20 DOUGLAS N. ARNOLD, RICHARD S. FALK, AND RAGNAR WINTHER tat σ σ div + u u + p p C in ( σ τ div + u v + p q ), τ Σ,v V,q Q σ σ + p p + u Π n n 1 u C( σ Π σ + p Π n p ), n 1 u u C( σ Π σ + p Π n p + u Π n u ), I u and σ are suiciently smoot, ten d(σ σ ) = dσ Π n dσ. σ σ + u u + p p C r+1 u r+2, div(σ σ ) C r+1 divσ r+1. Remark. Note tat te errors σ σ and u Π n u depend on te approximation o bot σ and p. For te coices made ere, te approximation o p is one order less tan te approximation o σ, and tus we do not obtain improved estimates, as one does in te case o te approximation o Poisson s equation, were te extra variable p does not enter. 8. A simpliied element Recall tat te lowest order element in te stable amily described above, or a discretization based on (1.5), is o te orm Σ P 1 Λ 2 (T ;V), V P 0 Λ 3 (T ;V), Q P 0 Λ 3 (T ;K). Te purpose o tis section is to present a stable element wic is sligtly simpler tan tis one. More precisely, te spaces V and Q are uncanged, but Σ will be simpliied rom ull linears to matrix ields wose tangential normal components on eac two dimensional ace o a tetraedron are only a reduced space o linears. We will still adopt te notation o dierential orms. By examining te proo o Teorem 7.1, we realize tat we do not use te complete sequence (5.2) or te given spaces. We only use te sequences (8.1) P 0 + Λ 2 (T ;K) P 0 Λ 3 (T ;K) 0, P 1 + Λ1 d (T ;V) P 1 Λ 2 d (T ;V) P 0 Λ 3 (T ;V) 0. Te purpose ere is to sow tat it is possible to coose subspaces o some o te spaces in (8.1) suc tat te desired properties still old. More precisely, compared to(8.1), te spaces P + 1 Λ 1 (T ;V) and P 1 Λ 2 (T ;V) are simpliied, wile te tree oters remain uncanged. I we denote by P + 1, Λ1 (T ;V) and P 1, Λ 2 (T ;V) te simpliications o te spaces P + 1 Λ1 (T ;V) and P 1 Λ 2 (T ;V), respectively, ten te properties we need are tat: (8.2) P + 1, Λ 1 (T ;V) d P 1, Λ 2 (T ;V) d d P 0 Λ 3 (T ;V) 0 is a complex and tat te surjectivity assumption (5.5) olds, i.e., S = S 1, maps te space P + 1, Λ1 (T ;V) onto P + 0 Λ2 (T ;K). Note tat i P + 0 Λ2 (T ;V) P 1, Λ 2 (T ;V), ten d maps P 1, Λ 2 (T ;V) onto P 0 Λ 3 (T ;V). We irst sow ow P + 1, Λ 1 (T ;V) can be constructed as a subspace o P + 1 Λ 1 (T ;V). Since te construction is done locally on eac tetraedron, we will sow ow to construct a space

21 MIXED METHODS FOR ELASTICITY WITH WEAKLY IMPOSED SYMMETRY 21 P 1, + Λ1 (T;V) as a subspace o P 1 + Λ1 (T;V). We begin by recalling tat te ace degrees o reedom o P 1 + Λ1 (T;V) ave te orm ω µ, µ P 0 Λ 1 (,V). We ten observe tat tis six dimensional space can be decomposed into P 0 Λ 1 (;V) = P 0 Λ 1 (;T )+P 0 Λ 1 (;N ), i.e. into orms wit values in te tangent space to, T or te normal space N. Tis is a dimensional decomposition. Furtermore, P 0 Λ 1 (;T ) = P 0 Λ 1 sym(;t )+P 0 Λ 1 skw(;t ), were µ P 0 Λ 1 (;T ) is in P 0 Λ 1 sym (;T ) i and only i µ(s) t = µ(t) s or ortonormal tangent vectors s and t. Finally, we obtain a 3+3 dimensional decomposition were P 0 Λ 1 (;V) = P 0 Λ 1 sym (;T )+P 0 Λ 1 skw (;V), P 0 Λ 1 skw(;v) = P 0 Λ 1 skw(;t )+P 0 Λ 1 (;N ). In more explicit terms, i µ P 0 Λ 1 (F;V) as te orm µ(q) = (a 1 t+a 2 s+a 3 n)q t+(a 4 t+a 5 s+a 6 n)q s, were t and s are ortonormal tangent vectors on, n is te unit normal to, and q is a tangent vector on, ten we can write µ = µ 1 + µ 2, wit µ 1 P 0 Λ 1 sym(;v) and µ 2 P 0 Λ 1 skw (;V), were ( µ 1 (q) = a 1 t+ a ) ( ) 2 +a 4 a2 +a 4 s q t+ t+a 5 s q s, 2 2 ( ) ( ) a2 a 4 a4 a 2 µ 2 (q) = s+a 3 n q t+ t+a 6 n q s. 2 2 Te reason or tis particular decomposition o te degrees o reedom is tat i we examine te proo o Lemma 6.1, were equation (6.3) is establised, we see tat te only degrees o reedom tat are used or an element ω P 1 + Λ1 (T;V) are te subset o te ace degrees o reedom given by: ω (S 0 ν), ν P 0Λ 0 (;K). However, or ν P 0 Λ 0 (;K), µ = S 0ν is given by µ(q) = νq. Since te general element ν P 0 Λ 0 (K) can be written in te orm b 1 (ts T st T ) + b 2 (nt T tn T ) + b 3 (ns T sn T ), νq = ( b 1 s+b 2 n)q t+(b 1 t+b 3 n)q s or q a tangent vector, and tus µ P 0 Λ 1 skw (;V). Hence, we ave split te degrees o reedom into tree on eac ace tat we need to retain or te proo o Lemma 6.1 and tree on eac ace tat are not needed. Te reduced space P + 1, Λ 1 (T;V) tat we now construct as two properties. Te irst is tat it still contains te space P 1 Λ 1 (T;V) and te second is tat te unused ace degrees o reedom are eliminated (by setting tem equal to zero). We can acieve tese conditions by irst writing an element ω P + 1 Λ 1 (T;V)asω = Π ω+(i Π )ω,wereπ denotesteusual projectionoperatorinto