arxiv: v1 [math.na] 11 Dec 2014

Transcription

1 DISCRETE KORN S INEQUALITY FOR SHELLS SHENG ZHANG arxiv: v1 [matna] 11 Dc 2014 Abstract W prov Korn s inqualitis for Nagdi and Koitr sll modls dfind on spacs of discontinuous picwis functions Ty ar usful in study of discontinuous finit lmnt mtods for slls Ky words Korn s inquality, Nagdi sll, Koitr sll, discontinuous finit lmnts Subjct classification 65N30, 46E35, 74S05 1 Introduction In t Nagdi sll modl, t strain nrgy is a sum of bnding, mmbran, and transvrs sar strain nrgis T tr strains ar xprssd in trms of t modl primary variabls comprising displacmnt of t sll mid-surfac and rotation of normal fibrs, in wic drivativs and function valus ar combind togtr in a complicatd mannr by t curvatur tnsors and Cristoffl symbols Korn s inquality for slls [5] stabliss an quivalnc btwn t strain nrgy norm of t primary variabls and t tir usual Sobolv norms, wic nsurs t wllposdnss of t t Nagdi sll modl in Sobolv spacs T situation for t Koitr sll modl is similar, wic xcluds t transvrs sar ffct in sll dformation and only uss mid-surfac displacmnt as t primary variabls, for wic a Korn s inquality is also availabl [5] Wn analyzing conforming finit lmnt mtods for slls [2, 4, 8], suc continuous vrsion of Korn s inqualitis play a fundamntal rol in tat ty attribut most of t numrical analysis and stimats to tat in Sobolv spacs To dal wit discontinuous Galrkin mtods for slls [16, 17], it is dsirabl to av appropriat gnralizations of t Korn s inqualitis tat can b applid to picwis dfind discontinuous functions It sms tat suc discrt Korn s inqualitis do not trivially follow from t Korn s inqualitis for slls Mtods of proving discrt Korn s inquality for plan lasticity [6] can not b asily adaptd to sll problms itr It is t purpos of tis papr to stablis suc discrt Korn s inqualitis for bot t Nagdi and Koitr sll modls, for t most gnral gomtry of sll mid-surfacs, witout making any additional assumption on t rgularity of sll modl solutions Dpartmnt of Matmatics, Wayn Stat Univrsity, Dtroit, MI (szang@wayndu) 1

2 2 SHENG ZHANG Lt Ω R 3 b t middl surfac of a sll of ticknss 2ǫ It is t imag of a twodimnsional domain Ω R 2 troug a mapping Φ tat as tird ordr continuous drivativs T Nagdi sll modl is a two-dimnsional modl dfind on Ω In t following, w us supr scripts to indicat contravariant componnts of tnsors, and us subscripts to indicat covariant componnts Grk supr and subscripts tak tir valus in {1, 2}, and Latin scripts tak tir valus in {1,2,3} Summation ruls wit rspct to rpatd and supr and subscripts will also b usd W us H α (τ) to dnot t L 2 basd Sobolv spac of ordr α of functions dfind on τ, in wic t norm and smi-norm ar dnotd by k,τ and k,τ, rspctivly W also us H 0 (τ) to dnot L 2 (τ) Wn τ = Ω, w simply writ t spac as H α T coordinats x α Ω furnis t curvilinar coordinats on Ω troug t mapping Φ W assum tat at any point on t surfac, along t coordinat lins, t two tangntial vctors a α = Φ/ x α ar linarly indpndnt T unit vctor a 3 = (a 1 a 2 )/ a 1 a 2 is normal to Ω T tripl a i furniss t covariant basis on Ω T contravariant basis a i is dfind by t rlations a α a β = δβ α and a3 = a 3, in wic δβ α is t Kronckr dlta It is obvious tat a α ar also tangnt to t surfac T mtric tnsor as t covariant componnts a αβ = a α a β, t dtrminant of wic is dnotd by a T contravariant componnts ar givn by a αβ = a α a β T curvatur tnsor as covariant componnts b αβ = a 3 β a α, wosmixdcomponnts arb α β = aαγ b γβ Tsymmtric tnsorc αβ = b γ α b γβ is calld t tird fundamntal form of t surfac T Cristoffl symbols ar dfind by Γ γ αβ = aγ β a α, wic ar symmtric wit rspct to t subscripts T Nagdi sll modl [15] uss t covariant componnts u α of t tangntial displacmnt vctor u = u α a α, cofficint w oft normal displacmnt wa 3 oft sll mid-surfac, and tcovariant componnts θ α of t normal fibr rotationvctor θ = θ α a α ast primary variabls T bnding strain tnsor, mmbran strain tnsor, and transvrs sar strain vctor associatd wit a dformation rprsntd by suc a st of primary variabls ar (11) ρ αβ (θ,u,w) = 1 2 (θ α β +θ β α ) 1 2 (bγ α u γ β +b γ β u γ α)+c αβ w, (12) γ αβ (u,w) = 1 2 (u α β +u β α ) b αβ w, (13) τ α (θ,u,w) = α w+b γ αu γ +θ α Hr, t covariant drivativ of a vctor u α is dfind by (14) u α β = β u α Γ γ αβ u γ T loading forcs on t sll body and uppr and lowr surfacs ntr t sll modl as rsultant loading forcs pr unit ara on t sll middl surfac, of wic t tangntial forc dnsity is p α a α and transvrs forc dnsity p 3 a 3 Lt t boundary Ω b dividd to

3 DISCRETE KORN S INEQUALITY FOR SHELLS 3 D Ω S Ω F Ω On D Ω t sll is clampd, on S Ω t sll is soft-simply supportd, and on F Ω t sll is fr of displacmnt constraint and subjct to forc or momnt only (Tr ar 32 diffrnt ways to spcify boundary conditions at any point on t sll boundary, of wic w considr t tr most typical) Lt H 1 = H 1 H 1 T sll modl is dfind in t Hilbrt spac (15) H = {(φ,v,z) H 1 H 1 H 1 ; v α and z ar 0 on D Ω S Ω, T modl dtrmins (θ,u,w) H suc tat (16) and θ α is 0 on D Ω} 1 a αβλγ ρ λγ (θ,u,w)ρ αβ (φ,v,z) adx 1 dx 2 3 Ω [ adx1 +ǫ 2 a αβλγ γ λγ (u,w)γ αβ (v,z)+κµ a αβ τ α (θ,u,w)τ β (φ,v,z)] dx 2 Ω Ω = (p α v α +p 3 z) ( ) adx 1 dx 2 + r α φ α d s+ q α v α +q 3 z +r α φ α d s Ω S Ω F Ω (φ,v,z) H Hr, q i and r α ar t forc rsultant and momnt rsultant on t sll dg [15] T factor κ is a sar corrction factor T last two intgrals on t sll dg is takn wit rspct to t arc lngt of t boundary of Ω T fourt ordr contravariant tnsor a αβγδ is t lastic tnsor of t sll, dfind by a αβγδ = µ(a αγ a βδ +a βγ a αδ )+ 2µλ 2µ+λ aαβ a γδ Hr, λ and µ ar t Lamé cofficints of t lastic matrial It satisfis t condition tat tr ar constants C 1 and C 2 tat only dpnd on t sll mid-surfac and t Lamé cofficints of t sll matrial suc tat for any tnsor ς αβ 2 α,β=1 ς αβ 2 C 1 a αβλγ ς αβ ς λγ, a αβλγ ς αβ ς λγ C 2 2 α,β=1 ς αβ 2 Also, tr ar constants C 1 and C 2 tat only dpnd on t sll mid-surfac suc tat for any vctor ς α, 2 2 ς α 2 C 1 a αβ ς α ς β, a αβ ς α ς β C 2 ς α 2 α=1 Ts inqualitis togtr wit a Korn s inquality [5] assurs tat t Nagdi sll modl (16) as a uniqu solution in t spac H T Korn s inquality stats tat tr is a C α=1

4 4 SHENG ZHANG tat could b dpndnt on t sll mid-surfac suc tat (17) θ H 1 + u H 1 + w H 1 [ C ρ αβ (θ,u,w) 2 L + γ 2 αβ (u,w) 2 L + τ 2 α (θ,u,w) 2 L 2 +f2 (θ,u,w) α,β=1 α,β=1 α=1 ] 1/2 (θ,u,w) H 1 H 1 H 1 Hr f is a continuous sminorm satisfying t rigid body motion condition tat if (θ,u,w) dfins a rigid body motion and f(θ,u,w) = 0 tn (θ,u,w) = 0 T displacmnt functions (θ,u,w) H 1 H 1 H 1 dfins a rigid body if and only if ρ αβ (θ,u,w) = 0, γ αβ (u,w) = 0, and τ α (θ,u,w) = 0 [5] W assum tat Ω is a boundd polygon Lt T b a sap rgular, but not ncssarily quasi-uniform, triangulation on Ω Sap rgularity of triangulations is a crucial notion in tispapr Itiswortwil torcall itsdfinitionr Considring atriangl, wltr andr b t radii of its inscribd circl and circumcircl, rspctivly Tn t ratio R/r is calld its sap rgularity constant, or simply sap rgularity For a triangulation, t maximum of sap rgularitis of all its triangls is calld t sap rgularity of t triangulation [9], dnotd by K W will nd to considr a (infinit) class of triangulations For a class, t sap rgularity K is t suprm of all t sap rgularitis of its triangulations For t triangulation T, w us T to dnot t st of all (opn) triangular lmnts of t partition, and us Ω to dnot t union of all t opn triangular lmnts W us E 0 to dnot t st of all intrior (opn) dgs and E all boundary dgs, and lt E = E 0 E W us τ to dnot t diamtr of an lmnt τ T and us to dnot t lngt of an dg E Lt H 1 b t spac of picwis H1 functions subordinatd to t triangulation T A function in H 1 is indpndntly dfind on vry lmnt τ T on wic it blongs to H 1 (τ) A function u in H 1 is crtainly in L2 (Ω ) On an dg E 0, a function u may av two diffrnt tracs from t two lmnts saring W us [u] to dnotd t absolut valu of diffrnc of t two tracs, wic is t jump of u ovr In t spac H 1, w dfin a norm (18) u H 1 := τ T u 2 1,τ + E 0 1 [u] 2 ds 1/2

5 DISCRETE KORN S INEQUALITY FOR SHELLS 5 Lt H 1 = H 1 H1 W prov t discrt Korn s inquality for Nagdi sll tat for all (θ,u,w) H 1 H1 H1 (19) θ 2 H + u 2 1 H + w 2 1 H [ 1 2 C ρ αβ (θ,u,w) 2 0,Ω + α,β=1 2 α,β=1 + E 0 γ αβ (u,w) 2 0,Ω + 1 α=1 2 τ α (θ,u,w) 2 0,Ω +f 2 (θ,u,w) α=1 2 ( [θα ] 2 + [u α ] 2) ds+ 1 E 0 [w] 2 ds Hr f is a sminorm suc tat f(θ,u,w) C( θ H 1 + u H 1 + w H 1 ) and it satisfis t rigid body motion condition tat if (θ,u,w) dfins a rigid body motion and f(θ,u,w) = 0 tn (θ,u,w) = 0 T constant C could b dpndnt on t sll mid-surfac and t sap rgularity K of t triangulation T, but otrwis is indpndnt of t triangulation W sall simply say tat suc a constant is indpndnt of T T Koitr sll modl [13] uss t covariant componnts u α of t tangntial displacmnt u α a α and t cofficint w of t normal displacmnt wa 3 of t sll mid-surfac as t primary variabl Suc a displacmnt dforms t surfac Ω and cangs its curvatur and mtric tnsors T linarizd cang in curvatur tnsor is t bnding strain tnsor It is xprssd in trms of t displacmnt componnts as (110) ρ K αβ (u,w) = 2 αβ w Γγ αβ γw +b γ α β u γ +b γ α u γ β +b γ β u γ α c αβ w T linarizd cang of mtric tnsor is t mmbran strain tnsor γ αβ (u,w) tat as t sam xprssion as in Nagdi modl (12) Tr is no transvrs sar in t Koitr modl Indd t bnding strain tnsor (110) can b obtaind from t bnding strain tnsor of t Nagdi modl (11) by liminating t rotation vctor θ using t zro-sar condition tat τ α (θ,u,w) = 0, and multiplying t rsult by 1 I, ρ K αβ (u,w) = ρ αβ(θ,u,w) wit θ α = α w b β α u β As in t Nagdi modl, w lt t boundary Ω b dividd to D Ω S Ω F Ω On D Ω t sll is clampd, on S Ω t sll is simply supportd, and on F Ω t sll is fr of displacmnt constraint and subjct to forc or momnt only (Tr ar 16 diffrnt ways to spcify boundary conditions at any point on t sll boundary, of wic w considr t tr most typical) T sll modl is dfind in t Hilbrt spac (111) H K = {(v,z) H 1 H 2 v α and z ar 0 on D Ω S Ω, and t normal drivativ of z is 0 on D Ω}

6 6 SHENG ZHANG T modl dtrmins (u,w) H K suc tat (112) 1 3 = Ω Ω a αβλγ ρ K λγ (u,w)ρk αβ (v,z) adx 1 dx 2 +ǫ 2 a αβλγ γ λγ (u,w)γ αβ (v,z) adx 1 dx 2 Ω (p α v α +p 3 z) ( adx 1 dx 2 + mdñzd s+ q α v α +q 3 z +mdñz ) d s S Ω F Ω (v,z) H K Hr, q i and marrsultant loading functions ont sll boundary, wic can bcalculatd fromforcrsultantsandmomntrsultantsontslldg[13] Tscalarz canbviwd as dfind on Ω W lt ñ = ñ α a α b t unit outward normal to Ω tat is tangnt to Ω T drivativ Dñz := ñ α α z is t dirctional drivativ in t dirction of ñ wit rspct to arc lngt T lastic tnsor a αβγδ is t sam as in Nagdi modl T wllposdnss of t Koitr modl (112) ings on a Korn s inquality for Koitr sll [5] tat tr is a constant C suc tat (113) u H 1 + w H 2 C [ α,β=1,2 ρ K αβ (u,w) 2 L 2 + α,β=1,2 γ αβ (u,w) 2 L 2 +f2 (u,w) ] 1/2 u H 1, w H 2 Hr f(u,w) is a smi-norm on H 1 H 2 tat satisfis t rigid body motion condition tat if (u,w) dfins a rigid body motion and f(u,w) = 0 tn u = 0 and w = 0 A displacmnt (u,w) H 1 H 2 dfins a rigid body motion of t sll mid-surfac if tr ar constant vctors canddsuc tat u α a α +wa 3 = c+d Φ(x 1,x 2 ) Tis isquivalnt toρ K αβ (u,w) = 0 and γ αβ (u,w) = 0 [5] For t Koitr modl, on t triangulation T, w also nd to considr picwis H 2 functions tat ar indpndntly dfind on ac lmnt, wit t norm dfind by (114) w H 2 := τ T w 2 2,τ + E 0 ( 1 ds) 1/2 [ α w] 2 ds+ 1 [w] 2 Lt f(u,w) b a smi-norm tat is continuous wit rspct to tis norm suc tat tr is a C only dpndnt on K of T and f(u,w) C( u H 1 + w H 2 ) (u,w) H 1 H2 W assum tat f satisfis t rigid body motion condition tat if (u,w) H 1 H 2 dfins a rigid body motion and f(u,w) = 0 tn u = 0 and w = 0 W av t discrt Korn s

7 DISCRETE KORN S INEQUALITY FOR SHELLS 7 inquality for Koitr sll tat for all u H 1 and w H 2 (115) u H 1 + w H 2 C + E 0 ( [ α,β=1,2 1 [u α ] 2 ds+ ρ K αβ (u,w) 2 0,Ω + 1 α,β=1,2 γ αβ (u,w) 2 0,Ω [ α w] 2 ds+ 1 [w] 2 ds ) +f 2 (v,z) Wit tis inquality stablisd, on may add t trm E 3 0 [w]2 ds to bot sids to obtain a nw inquality Lt [ s w] and [ n w] b t jumps in t tangntial and normal drivativs of w ovr and dg, rspctivly Tn w av t idntity 2 [ α w] 2 = [ s w] 2 + [ n w] 2 α=1 If w is a picwis polynomial, w av t invrs inquality [ s w] 2 C 2 [w]2 Tus for picwis polynomials u α and w, w av t following variant of discrt Korn s inquality for Koitr sll (116) u H 1 + w H2 C Hr + E 0 ( (117) w H2 := τ T 1 [ α,β=1,2 w 2 2,τ + E 0 ρ K αβ (u,w) 2 0,Ω + [u α ] 2 ds+ 1 ( 1 α,β=1,2 γ αβ (u,w) 2 0,Ω [ n w] 2 ds+ 3 [w] 2 ds ) ds) 1/2 [ n w] 2 ds+ 3 [w] 2 +f 2 (u,w) Tis inquality is usful for analysis of discontinuous Galrkin mtods for Koitr sll In proving ts discrt Korn s inqualitis, a important tool is a compact mbdding torm in t spac H 1, wic is provd in Sction 2 In Sction 3 w prov t discrt Korn s inqualitis (19) for t Nagdi sll, and Sction 4 is dvotd Koitr sll modl Trougout t papr, C will b a gnric constant tat may dpnd on t domain Ω, t mapping Φ tat dfins t sll mid-surfac, and sap rgularity K of a triangl, of a triangulation, or of a class of triangulations But otrwis, t constant is indpndnt of triangulations An intgral u(x τ 1,x 2 )dx 1 dx 2 or v(s)ds will b simply writtn as u or Ω τ v, in wic t intgration variabl and masur sould b clar from t contxt Ω 1/2 1/2

8 8 SHENG ZHANG 2 Compact mbdding in t spac of picwis H 1 functions W nd a discrt analogu of t Rllic-Kondracov compact mbdding torm, wic will play a fundamntal rol in proving t discrt Korn s inqualitis for slls Tr ar svral paprs rlvant to compact mbdding in picwis function spacs In [11], suc a rsult is statd undr t assumption tat t triangulation T is quasi-uniform In [10], a torm was provd for picwis polynomials undr t assumption tat t maximum ms siz tnds to zro In [7], tr is a sub-ms condition to b vrifid Altoug ts toris ar dvlopd in mor gnral sttings, tir rsults do not radily mt our nds W only assum tat T is a sap rgular triangulation of t polygon Ω Our proof also clarly sows tat t statmnt could brak down wn an intrior angl of Ω tnds to zro or 2π 21 A trac torm W first prov a tractorm forfunctions inh 1 Tis rsult itslf is a gnralization of a trac torm of Sobolv spac tory It will b usd in proving t discrt compact mbdding torm, and in vrifying t continuity of t sminorm dnotd by f in t rigt and sid of (19), (115), and (116) W will nd t following trac torm on an lmnt [3] Lmma 21 Lt τ b a triangl, and on of its dgs Tn tr is a constant C dpnding on t sap rgularity of τ suc tat [ 2 ] (21) u 2 C u 2 + α u 2 u H 1 (τ) 1 τ α=1 Torm 22 Lt T b a sap rgular, but not ncssarily quasi-uniform triangulation of Ω Tr xists a constant C suc tat (22) u L 2 ( Ω) C u H 1 u H 1 Proof Lt φ b a picwis smoot vctor fild on Ω wos normal componnt is continuous across any straigt lin sgmnt, and suc tat φ n = 1 on Ω On ac lmnt τ T, w av u 2 φ n = div(u 2 φ) = (2u u φ+u 2 divφ) τ τ Hr u is t gradint of u Summing up t abov quations ovr all lmnts of T, w gt Ωu 2 = [u 2 φ]+ (2u u φ+u 2 divφ) E 0 Ω If is t bordr btwn t lmnts τ 1 and τ 2 wit outward normals n 1 and n 2, tn [u 2 φ] = u 2 1φ 1 n 1 +u 2 2φ 2 n 2, wru 1 andu 2 arrstrictionsofuonτ 1 andτ 2, rspctivly It is notd tat altoug φ may b discontinuous across, it normal componnt is continuous, τ τ

9 DISCRETE KORN S INEQUALITY FOR SHELLS 9 i, φ 1 n 1 +φ 2 n 2 = 0 On t dg, w av [u 2 φ] [u 2 ] φ 0,,Ω Hr, [u 2 ] = u 2 1 u2 2 = 2 [u][{u}], wit [{u}] = (u 1 +u 2 )/2 bing t avrag W av (23) [u 2 φ] 2 φ 0,,Ω [ 1 [u] 2 ] 1/2 [ C φ 0,,Ω [ [{u}] 2 ] 1/2 1 [u] 2 ] 1/2 [ δ u δ u 2 ] 1/2 Hr, φ 0,,Ω is t Sobolv norm in t spac [W 0, (Ω)] 2 and C only dpnds on t sap rgularity of τ 1 and τ 2 W usd δ = τ 1 τ 2 to dnot t co-boundary of dg, and w usd t trac stimat (21) It tn follows from t Caucy Scwarz inquality tat u 2 L 2 ( Ω) C( φ 0,,Ω + divφ 0,,Ω ) u 2 L + u Ω E 0 1 [u] 2 Hr t constant C only dpnds on t sap rgularity of T T dpndnc on Ω of t C in (22) is iddn in t φ in t abov inquality 22 Compact mbdding in H 1 For δ > 0, w dfin a boundary strip Ω δ of widt O(δ) for t domain Ω Lt Ω 0 δ = Ω\Ω δ b t intrior domain T intrior domain Ω 0 δ as t proprty tat if a point is in Ω 0 δ, tn t disk cntrd at t point wit radius δ ntirly lis in Ω W first sow tat wn t strip is tin, t L 2 (Ω δ ) norm of t rstriction on Ω δ of a function in H 1 must b small Lmma 23 Tr is a constant C suc tat (24) u 2 Cδ u 2 H 1 Ω δ Hr Ω δ is a boundary strip of widt δ attacd to Ω u H 1 Proof W coos a picwis smoot a vctor fild φ wos normal componnt is continuous across any curv suc tat φ = 0 on t innr boundary of Ω δ, and divφ = 1 and φ Cδ on Ω δ (A construction of suc φ is givn blow) W tn xtnd φ by zro onto t ntir domain Ω T xtndd, still dnotd by φ, is a picwis smoot vctor fild wos normal componnts is continuous ovr any curv in Ω W tus av u 2 = u 2 divφ = u 2 divφ = 2u u φ+ u 2 φ n Ω δ τ T τ τ T τ τ T τ Ω T last trm can b writtn as u 2 φ n = τ T τ Ωu 2 φ n+ [u 2 φ] E 0

10 10 SHENG ZHANG Sinc t normal componnts of φ is continuous on dgs in E 0, w us t sam argumnt as in t proof of Torm 22, cf, (23), to gt [u 2 φ] C φ 0,,Ω u 2 L τ u [u] 2 τ τ T E 0 E 0 It tn follows from Torm 22 tat (25) u 2 C φ 0,,Ω u 2 L + u [u] 2 2 Ω δ Ω E 0 T proof is complt sinc φ 0,,Ω Cδ T constant C dpnds on Ω in trms of its intrior angls and xtrior angls at convx and concav vrtxs, rspctivly W dscrib a way to coos t boundary strip and construct t vctor fild φ tat was usd in t proof Tis fild can b constructd by picing togtr svral spcial vctor filds W nd som vctor filds on rctangls, wdgs, and circular disks On t xy-plan, w considr t vrtical rctangular strip R = (0,δ) (0,l) On tis strip, w considr ψ R = x,0 Tis vctor fild satisfis t condition tat divψ R = 1, ψ R = 0 on t lft sid, ψ R n = 0 on t top and bottom sids and t maximum of ψ R is δ tat is attaind on t rigt sid On a wdg W wit its vrtx at t origin, w considr t vctor fild ψ W = x,y /2 Tis ψ W satisfis t conditions tat divψ W = 1, ψ W n = 0 on t two sids of W, and ψ W = ρ/2 at a point in W wos distant from t origin is ρ On a circular disk C cntrd at t origin and of radius ρ, w considr t vctor fild ψ C = (1 ρ 2 /r 2 ) x,y /2 Hr r = (x 2 + y 2 ) 1/2 Tis vctor fild satisfis t condition tat divψ C = 1 on t disk xcpt at t cntr wr it is singular It points toward t cntr vry wr And it is zro on t boundary W us tis fild on a sctor of t circl C, on t two radial sids of wic w av ψ C n = 0 Wit ts spcial vctor δ θ θ Ω αδ Ω Figur 1 Boundary strip Ω δ nar a convx vrtx (lft) and a concav vrtx (rigt) filds, w can tn assmbl t φ on a boundary strip Ω δ Along t intrior sid of ac straigt sgmnt of Ω w coos a uniform strip of ticknss δ Ts strips ovrlap nar vrtxs If Ω is convx at a vrtx, w introduc a wdg wos vrtx is at t intrsction

11 DISCRETE KORN S INEQUALITY FOR SHELLS 11 of t intrior boundary of t uniform strips, and wos sids ar ortogonal to t mting straigt sgmnts, s t lft figur of Figur 1 If Ω is concav at a vrtx, w rsolv it by using a circular sctor, cntrd nar t vrtx and outsid of Ω T radius of t circl is sligtly biggr tan δ suc tat t arc is continuously connctd to t intrior dgs of t mting strips, and t two radial sids ar ortogonal to t mting boundary sgmnts, s t rigt figur in Figur 1 Wit suc tratmnt of t vrtics, t boundary strip Ω δ is composd of rctangular strips attacing to major portions of straigt sgmnts of Ω, portion of wdgs at convx vrtxs, and portion of circular sctors at concav vrtxs, s t sadd rgion in Figur 1 W tn transform ψ R, ψ W, and ψ C to various parts of Ω δ and assmbl a φ on Ω δ T vctor fild φ tus constructd is zro on t intrior boundary of Ω δ Its normal componnts ar continuous across any curv, and divφ = 1 on Ω δ T ticknss of Ω δ is t constant δ for t rctangular part It is maximizd to δ/sin θ at a 2 convx vrtx It is minimizd to αδ at t concav vrtx, wit 0 < α < 1, a valu can b cosn as, for xampl, 1/2 T norm φ as a maximum δ/sin θ at a convx vrtx wit θ 2 bing t intrior angl Tus wn θ is small, φ is big, and t stimat (25) dtriorats T norm φ also as a local maximum nar a concav vrtx It is boundd as φ δ 1 αsin θ 2 (1 α)sin θ 2 Wn t xtrior angl is small tis maximum would b big Also, on nds to coos α away from 1 and 0, to maintain a modrat ticknss of t strip and a rasonabl bound for φ wic affct t stimat (25) W rmark tat t constant C in t stimat (24) could tnd to infinity if an intrior angl of t polygonal domain Ω tnds to zro or 2π W tn prov tat functions in H 1 ar sift-continuous in L2, as statd in t following lmma W xtnd a function u H 1 to a function ũ on t wol R2 by zro Lmma 24 Tr is a constant C suc tat (26) [ũ(x+ρ) ũ(x)] 2 dx C ρ u 2 H 1 R 2 u H 1 Proof Bcaus for an lmnt τ T, smoot functions ar dns in H 1 (τ), w only nd to prov (26) for functions tat ar smoot on ac lmnt of T Lt u b suc a picwis smoot function Lt ρ b an arbitrary sort vctor W tak δ = ρ and coos a boundary strip Ω δ T intrior part Ω 0 δ of t domain as t proprty tat if x Ω0 δ tn t lin sgmnt [x,x+ρ] Ω W av [ũ(x+ρ) ũ(x)] 2 dx = R 2 Ω 0 δ [u(x+ρ) u(x)] 2 dx+ R 2 \Ω 0 δ [ũ(x+ρ) ũ(x)] 2 dx

12 12 SHENG ZHANG Using Lmma 23, w bound t scond trm as (27) [ũ(x+ρ) ũ(x)] 2 dx u 2 (x)dx C ρ u 2 H 1 Ω 2δ R 2 \Ω 0 δ W tn focus on t first trm Tis intgral can b takn on an qual masur subst of Ω 0 δ Tis subst is obtaind by rmoving a zro masur subst tat is composd of suc point x: x or x+ρ is on an opn dg E 0, or t closd straigt lin sgmnt [x,x+ρ] conncting x and x+ρ contains any vrtx of T, or [x,x+ρ] ovrlaps som dgs of E 0 By suc xclusion, for any x in t rmaining st, bot t nds of [x,x + ρ] ar in t intrior of som opn triangular lmnts, and [x, x + ρ] contains no vrtx T rstriction of u on [x,x+ρ] is a picwis smoot on dimnsional function, wic may av a finit numbr of jumping points in t opn straigt lin sgmnt (x,x+ρ) By t fundamntal torm of calculus, w av u(x+ρ) u(x) = 1 0 u(x+tρ) ρdt+ p [x,x+ρ] E 0 Not tat t intgrand in t intgral may mak no sns at t, if x+tρ E 0, wr u may jump Ts points ar xcludd from t intgration, wr t jumping ffct is rsolvd in t scond trm On t sgmnt [x,x+ρ], u may av a jump at p [x,x+ρ] E 0, wic is dnotd by [u] p tat is t valu of u from t sid of x minus tat from t sid of x+ρ W tus av 1 [u(x+ρ) u(x)] 2 ρ 2 u(x+tρ) 2 dt+ 0 p [x,x+ρ] E 0 [u] p [u] p Wn w tak intgral on Ω 0 δ (minus t aformntiond zro-masur subst), t first trm is boundd as follows 1 1 ρ 2 u(x+tρ) 2 dtdx = ρ 2 u(x+tρ) 2 dxdt ρ Ω 2 u 2 dx Ω 0 δ 0 0 To stimat t jumping rlatd scond trm, w writ [u] p = 1/2 1/2 [u] p if p [x,x+ ρ] E 0 is on t dg, and us t Caucy-Scwarz inquality to obtain t following stimat 2 [u] p 1 p [x,x+ρ] E 0 [x,x+ρ] Ω 0 δ [u] 2 [x,x+ρ] [x,x+ρ] 2

13 DISCRETE KORN S INEQUALITY FOR SHELLS 13 W sow blow tat tr is a C, dpnding on t domain Ω and t sap rgularity K of T, suc tat (28) C W tn av Ω 0 δ p [x,x+ρ] E 0 [u] p [x,x+ρ] 2 C Ω 0 δ [x,x+ρ] 1 [u] 2 [x,x+ρ]dx Evry trm in t rigt and sid is associatd wit a particular dg E 0 Eac dg E 0 is rlvant to at most t points in t paralllogram Ω in t Figur 2 Tus, by canging t ordr of sum and intgral, w av Ω 0 δ [x,x+ρ] 1 [u] 2 [x,x+ρ] dx E 0 1 E 0 Ω [u] 2 [x,x+ρ] dx 1 sin ρ, ρ [u] 2 ρ E 0 Hr, ρ, is t angl btwn t vctor ρ and t dg Trfor, w av 1 [u] 2 x Ω x + ρ Figur 2 A ρ sift and Ω for an dg Ω 0 δ[u(x+ρ) u(x)] 2 dx ρ 2 u 2 0,Ω + ρ E 0 1 [u] 2 Not tat t scond trm may carry t smallr cofficint ρ max{, ρ } suc tat t two trms ar closr in t ordr But w do not nd suc rfind stimats W tus provd (29) Ωδ[u(x+ρ) u(x)] 2 dx C ρ u 2 H u H 1 1 0

14 14 SHENG ZHANG T sift continuity (26) tn follows from (29) and (27) W av sown tat t st of zro xtndd functions is sift continuous in L 2 (R 2 ) Wgiv a prooffor tstimat (28) Lt l b astraigt lin cutting trougω Lt T b a sap rgular triangulation wit rgularity constant K Tn t sum of lngts of ms lin sgmnts intrscting l is boundd indpndnt of t triangulation Mor spcifically, w prov tat tr is a constant C, dpnding on t sap rgularity K, but otrwis indpndnt of t triangulation T suc tat (210) C E and l W sall us som facts tat follow from t sap rgularity assumption Tr is a minimum angl θ K among all angls of triangls of T T numbr of dgs saring a vrtx is boundd by a constant C tat only dpnds on K Lt 1 and 2 b two dgs saring a vrtx Tr ar constants C 1 and C 2 dpnding on K suc tat 1 C 1 2 and 2 C 2 1 Witout loss of gnrality, w assum l is orizontal W also assum tat l dos not pass any vrtx (Tis assumption was alrady mad in t contxt of (28) It can b rmovd by a sligt modification of t following argumnt) W first trim t st of intrscting A C b E c a B D F Figur 3 A lin cutting troug t triangulation dgs { E and l } to simplify t st witout significantly rducing its sum of lngts of all dgs Considr t lft most dg intrscting l, of wic only on nd vrtx is sard by som otr dgs intrscting l Lt A b tis nd vrtx, and assum it is abov l W xamin all t dgs intrscting l and saring A in t countrclockwis ordr W discard all suc dgs but t last on tat is AB in Figur 3 (T nxt dg saring A, as AC, dos not intrsct l) T dg BC must intrsct l Tr could b otr dgs intrscting l and saring t vrtx B W xamin all t dgs saring B and intrscting l in t clockwis ordr W discard all but t last on (It is BC in t figur) Now t vrtx C is in t sam situation as A, and w can dtrmin t dg CD using t sam rul as for AB Tn w dtrmin DE, EF, and so fort W do t trimming all t way to t rigt nd of l Tis procdur toucs vry dg intrscting l, by itr trimming an dg off or kping it Evry dg dltd as at last on vrtx-saring dg rtaind T rmaining dgs constitut a continuous picwis straigt pat as rprsntd by t

15 DISCRETE KORN S INEQUALITY FOR SHELLS 15 tick lin in t figur W dnot tis st by E l It follows from t aformntiond facts about t sap rgular triangulation tat tr is a constant C, dpnding on K only, suc tat E and l C E l W considr a typical triangl boundd by l and E l, as t sadd on in t figur, wos sids ar a, b, and c Lt t angl DEF b dnotd by θ Tn θ θ K Not tat c 2 = a 2 + b 2 2abcosθ If θ is obtus, tn a + b 2c Otrwis, w av c 2 = (a 2 + b 2 )(1 cosθ) +(a b) 2 cosθ (a 2 +b 2 )(1 cosθ) Tus a+b 2 cas, w av a+b 1 cosθ K c W tus provd 2 l Ω 1 cosθ K E l E0 2 1 cosθ c In any From tis, (210) follows Hr l Ω is t lngt of t lin sgmnt l Ω wic dos xcd t diamtr of Ω W can now prov t following compact mbdding torm Torm 25 Lt T i b a (infinit) class of sap rgular but not ncssarily quasi-uniform triangulations of t polygonal domain Ω, wit a sap rgularity constant K For ac i, lt H 1 i b t spac of picwis H 1 functions, subordinatd to t triangulation T i, quippd wit t norm (18) Lt {u i } b a squnc suc tat u i H 1 i for ac i and tr is a constant C, suc tat u i H 1 i C for all i Tn, t squnc {u i } as a convrgnt subsqunc in L 2 Proof It follows from (26) tat t squnc {u i } is a sift-continuous subst of L 2 Tis, togtr wit t stimat (22), vrifis t condition for a subst of L 2 to b compact, s Torm 212 in [1] 3 Discrt Korn s inquality for Nagdi sll In tis sction, w prov t discrt Korn s inquality (19) for t Nagdi sll modl W lt H = H 1 H 1 H 1, and dfin a norm in tis spac by [ ( ) 1/2 (31) (θ,u,w) H := θ α 2 H + u 1 α 2 H + w 2 1 H ] 1 Rcall tat t norm in H 1 is dfind by (18) Lt f(θ,u,w) b a smi-norm tat is continuous wit rspct to t H norm suc tat tr is a C and (32) f(θ,u,w) C (θ,u,w) H (θ,u,w) H

16 16 SHENG ZHANG W also assum tat f satisfis t condition tat if (θ,u,w) H 1 H 1 H 1 dfins a rigid body motion, as xplaind in t introduction, and f(θ,u,w) = 0 tn θ = 0, u = 0, and w = 0 W dfin a discrt nrgy norm on t spac H by (33) (θ,u,w) E := [ α,β=1,2 ( ραβ (θ,u,w) 2 0,Ω + γ αβ (u,w) 2 0,Ω ) + + E 0 1 ( T main rsult of tis sction is t following torm τ α (θ,u,w) 2 0,Ω ( [θα ] 2 + [u α ] 2) + [w] )+f 2 2 (θ,u,w) Torm 31 W assum tat Ω R 2 is a polygon, on wic T is a sap rgular but not ncssarily quasi-uniform triangulation wit a rgularity constant K Lt f(θ, u, w) b a sminorm on H satisfis t condition (32), and t rigid body motion condition Tr xists a constant C tat could b dpndnt on t sll mid-surfac and sap rgularity of t triangulation, but otrwis indpndnt of t triangulation, suc tat (34) (θ,u,w) H C (θ,u,w) E (θ,u,w) H To prov tis torm, w nd a discrt Korn s inquality for plan lasticity for picwis functions in H 1, s (121) of [6] It says tat tr is a constant C tat migt b dpndnt on t domain Ω and t sap rgularity K of t triangulation T, but otrwis indpndnt of T suc tat (35) u α 2 H C u 1 α 2 0,Ω + αβ (u) α,β=1,2 2 0,Ω + 1 [u α ] 2 E 0 Hr αβ (u) = ( β u α + α u β )/2 is t plan lasticity strain wic is t symmtric part of t gradint of u Proof of Torm 31 From (35), t dfinitions (11), (12), and (13) of ρ αβ, γ αβ, and τ α, and t dfinition of t discrt nrgy norm (33), w s tat tr is a constant C suc tat [ (36) (θ,u,w) 2 H C (θ,u,w) 2 E + ] ( ) θα 2 0,Ω + u α 2 0,Ω + w 2 0,Ω 1/2 (θ,u,w) H

17 DISCRETE KORN S INEQUALITY FOR SHELLS 17 Onafixd triangulationt, ittnfollowsfromtrllic Kondracov compact mbdding torm and Ptr s lmma (Torm 21, pag 18 in [12]) tat tr is a constant C T suc tat (θ,u,w) H C T (θ,u,w) E (θ,u,w) H W nd to sow tat for a class of sap rgular triangulations, suc C T as an uppr bound tat only dpnds on t sap rgularity K of t wol class If tis is not tru, tr would xist a squnc of triangulations {T n } and an associatd squnc of functions (θ n,u n,w n ) in H n suc tat (37) (θ n,u n,w n ) Hn = 1 and (θ n,u n,w n ) En 1/n It follows from Torm 25 tat tr is a subsqunc, still dnotd by (θ n,u n,w n ), suc tat (38) lim n (θ n,u n,w n ) = (θ 0,u 0,w 0 ) in L 2 L 2 L 2 W sow tat t limiting functions θα 0, u0 α, and w0 ar all in H 1, tis limit dfins a rigid body motion, and it is zro, wic will lad to a contradiction First, w sow tat w 0 is actually in H 1 and w av tat α w 0 +θα 0 + b β αu 0 β = 0 Sinc lim n (θ n,u n,w n ) = (θ 0,u 0,w 0 ) in L 2 L 2 L 2, in viw of t dfinition (13), w av lim τ α(θ n,u n,w n ) = τ α (θ 0,u 0,w 0 ) in H 1 n In t abov xprssions, t drivativs α w n and α w 0 ar undrstood in distributional sns For any compactly supportd smoot functions φ α D(Ω), for ac n, w av τ α (θ n,u n,w n ),φ α = w n α φ α + (θα n +bβ α un β )φα Ω n Ω n = ( α w n +θα n +b β αu n β)φ α [w n ] nα φ α Ω n Hr, n α ar t componnts of t unit normal n to t dg If is sard by τ 1 and τ 2 wit unit outward normals bing n 1 and n 2, on wic t rstrictions of w ar w 1 and w 2, tn [w] nα = w 1 n 1α +w 2 n 2α is t jump of w wit rspct to n Summation convntion is also usd in [w n ] nα φ α wit α bing viwd as a rpatd sub and supr scripts Using Höldr inquality and t trac inquality (21), w gt E 0 n τ α (θ n,u n,w n ),φ α τ α (θ n,u n,w n ) 0,Ωn φ α 0,Ω +C 1 [w n ] 2 E 0 n 1/2 2 φ α 2 0,Ω + 2 τ φ α 2 1,τ τ T n α=1 1/2

18 18 SHENG ZHANG Sinc (θ n,u n,w n ) En 0, w av lim τ α(θ n,u n,w n ) = 0 in H 1 n Trfor, in H 1, w av τ α (θ 0,u 0,w 0 ) = α w 0 +θα 0 +bβ α u0 β = 0 Sinc θ0 α and u0 α ar in L 2, tis quation sows tat t wak drivativs of w 0 ar in L 2 Tus w av w 0 H 1 and (39) τ α (θ 0,u 0,w 0 ) = α w 0 +θ 0 α +bβ α u0 β = 0 Nxt w sow u 0 α H1 and γ αβ (u 0,w 0 ) = 0 W av lim n γ αβ (u n,w n ) = γ αβ (u 0,w 0 ) in H 1, in wic t drivativs α u n β and αu 0 β ar undrstood in t distributional sns Lt φ αβ b an arbitrary symmtric tnsor valud function wit componnts in D(Ω) For ac n w av γ αβ (u n,w n ),φ αβ = u n α β φ αβ (Γ λ αβu n λ +b αβ w n )φ αβ Ω Ω = γ αβ (u n,w n )φ αβ [u n α] nβ φ αβ Ω n It follows from tis quation and t assumption lim n (θ n,u n,w n ) En = 0 tat w av lim n γ αβ (u n,w n ),φ αβ = 0 Tus γ αβ (u 0,w 0 ) = 0, in wic t drivativs ar undrstood in t distributional sns In viw of t dfinition (12), w av αβ (u 0 ) = Γ λ αβu 0 λ +b αβ w 0 Sinc u 0 α and w 0 ar in L 2, so αβ (u 0 ) L 2 In t sns of distribution, tr is t idntity tat [5, 14] αβ u 0 λ = β αλ (u 0 )+ α λβ (u 0 ) λ αβ (u 0 ) Fromtiswstat αβ u 0 λ H 1 ItfollowsfromaLmmaofJLLions(wosassumption on t domain is mt by our polygon, s pag 110 of [14] and pag 124 of [5]) and t fact β u 0 λ H 1 tat β u 0 λ L2 Trfor, w provd tat u 0 α H1, and w av (310) γ αβ (u 0,w 0 ) = 0 W tn sow tat θα 0 H1 and ρ αβ (θ 0,u 0,w 0 ) = 0 As in t abov, in t spac H 1, w av lim n ρ αβ (θ n,u n,w n ) = ρ αβ (θ 0,u 0,w 0 ), in wic t drivativs ar undrstood in t distributional sns For arbitrary φ αβ D(Ω), and for any n, w av t idntity ρ αβ (θ n,u n,w n ),φ αβ = ρ αβ (θ n,u n,w n )φ αβ [θα n ] nβ φ αβ + b γ α [un γ ] nβ φ αβ Ω n E 0 n E 0 n E 0 n Using t assumption lim n (θ n,u n,w n ) En 0 and tis quation, w gt lim ρ αβ(θ n,u n,w n ),φ αβ = 0 n

19 DISCRETE KORN S INEQUALITY FOR SHELLS 19 Trfor ρ αβ (θ 0,u 0,w 0 ) = 0, in wic t drivativs α θβ 0 ar undrstood in in t distributional sns In viw of t dfinition (11) of ρ αβ (θ 0,u 0,w 0 ), and t provd rgularity tat u 0 α H1, w s tat αβ (θ 0 ) ar in L 2 From tis and t argumnt usd abov, w s αβ 2 θ γ H 1 Using t Lmma of J L Lions again, w gt t rgularity θα 0 H 1, and t quation (311) ρ αβ (θ 0,u 0,w 0 ) = 0 According to Lmma 34 of [5], using t rgularitis tat θ 0 α, u0 α, and w0 ar all in H 1, and t quations (39), (310), and (311), w conclud tat t displacmnt functions (θ 0,u 0,w 0 ) dfins a rigid body motion Finally, w sow tat (θ 0,u 0,w 0 ) = 0 It follows from t bound (36), t assumption (37), and t convrgnc (38) tat lim n (θ n θ 0,u n u 0,w n w 0 ) Hn = 0 Sinc f is uniformly continuous wit rspct to t norm of H n and sinc f(θ n,u n,w n ) 0 (f is a part of t nrgy norm (33)), w s f(θ 0,u 0,w 0 ) = 0 Tus (θ 0,u 0,w 0 ) = 0 Trfor, lim n (θ n,u n,w n ) Hn = 0, wic is contradict to t assumption tat (θ n,u n,w n ) Hn = 1 As an xampl, w tak f(θ,u,w) = E D θα 2 + E S ED ( u 2 α + ) 1/2 w 2 It follows from Torm 22 tat tr is a C suc tat t continuity condition (32) is satisfid by tis f Undr t assumption tat t masur of D Ω is positiv, it is vrifid in [5] tat if (θ,u,w) dfins a rigid body motion and f(θ,u,w) = 0 tn θ = 0, u = 0, and w = 0 Wit tis f in t Korn s inquality (34), w could add boundary pnalty trm E D 1 θ 2 α + E S ED ( 1 u 2 α + 1 w 2 ) to t squars of bot sids of (34) to obtain an inquality tat is usful in analysis of discontinuous Galrkin mtods for Nagdi sll in wic t ssntial boundary conditions ar nforcd by Nitsc s consistnt boundary pnalty mtod

20 20 SHENG ZHANG 4 Discrt Korn s inquality for Koitr sll T wllposdnss of t Koitr modl (112) is basd on t Korn s inquality for Koitr sll tat tr is a constant C suc tat [ (41) u α 2 H 1 + w 2 H C ρ K 2 αβ (u,w) 2 L + ] γ 2 αβ (u,w) 2 L 2 +f2 (u,w) α,β=1,2 α,β=1,2 u H 1, w H 2 Hr f(u,w) is a smi-norm on H 1 H 2 tat satisfis t condition tat if u,w dfins a rigid body motion, as xplaind in t introduction, and f(u,w) = 0 tn u = 0 and w = 0 W gnraliz tis inquality to picwis functions on Ω For Koitr sll problms w also nd to us t discrt spac H 2 tat is composd of picwis H2 functions wit t norm dfind by (114) For (u,w) H 1 H2, w dfin a norm ( (42) (u,w) H K = 1/2 u α 2 H + w 2 1 H ) 2 Lt f(u,w) b a smi-norm tat is continuous wit rspct to tis norm suc tat tr is a C tat only dpnds on t sll mid-surfac and t rgularity K of t triangulation T (43) f(u,w) C (u,w) H K (u,w) H 1 H2 W assum tat f satisfis t rigid body motion condition tat if (u,w) H 1 H 2 dfins a rigid body motion and f(u,w) = 0 tn u = 0 and w = 0 W dfin t discrt nrgy norm on t spac H K (44) (u,w) E K = + E 0 [ α,β=1,2 ( ρ K αβ (u,w) 2 0,Ω + 1 [u α ] 2 + α,β=1,2 γ αβ (u,w) 2 0,Ω 1 [ α w] [w] )+f 2 2 (u,w) W tn av t following gnralization of t Korn s inquality for Koitr slls to picwis functions Torm 41 Tr xists a constant C tat could b dpndnt on t sll mid-surfac, and t sap rgularity K of T, but otrwis indpndnt of t triangulation, suc tat (45) (u,w) H K C (u,w) E K w H 2, u α H 1 1/2

21 DISCRETE KORN S INEQUALITY FOR SHELLS 21 Proof Koitr s sll modl is a rstriction of t Nagdi s modl on t subspac of zrosar dformations T discrt Korn s inquality (45) for Koitr sll can b drivd from t discrt Korn s inquality for Nagdi sll (34) For u α H 1 and w H2, w dfin picwis function θ α H 1 by (46) θ α = α w b γ α u γ W tn av ρ αβ (θ,u,w) = ρ K αβ (u,w) and τ α(θ,u,w) = 0 Tr ar constants C 1 and C 2 only dpnding on componnts of t mixd curvatur tnsor b γ α suc tat for all E 0 and for all τ T ( C 1 [u α ] 2 + ) [ α w] 2 [u α ] 2 + [θ α ] 2 C 2 ( [u α ] 2 + [ α w] 2 ), ( 2 ) C 1 u α 1,τ + w 2,τ α=1 2 2 θ α 1,τ + u α 1,τ + w 1,τ α=1 ( 2 ) C 2 u α 1,τ + w 2,τ α=1 From ts, w gt t quivalncs α=1 (47) C 1 (u,w) H K (θ,u,w) H C 2 (u,w) H K and (48) C 1 (u,w) E K (θ,u,w) E C 2 (u,w) E K T continuity (43) implis t condition (32) T inquality (45) tn follows from t Korn s inquality for Nagdi sll (34) As an xampl for t smi-norm satisfying t continuity condition (43), w tak f(u,w) = ( ) u 2 α + w 2 + 1/2 (D n w) 2 E S ED It follows from Torm 22 tat tr is a C only dpndnt on K suc tat t continuity condition (43) is satisfid by tis f Undr t assumption tat t masur of D Ω is positiv, it is vrifid in [5] tat if (u,w) H 1 H 2 dfins a rigid body motion and f(u,w) = 0 tn u = 0 and w = 0 E D

22 22 SHENG ZHANG Rfrncs [1] RA Adams, Sobolv spacs, Acadmic Prss, Nw York, 1975 [2] DN Arnold, F Brzzi, Locking fr finit lmnt mtods for slls, Mat Comp, 66 (1997), pp 1-14 [3] DN Arnold, F Brzzi, B Cockburn, LD Marini, Unifid analysis of discontinuous Galrkin mtods for lliptic problms, SIAM J Numr Anal, 39 (2002), pp [4] M Brnadou, Finit lmnt mtods for tin sll problms, Jon Wily & Sons - Masson, 1996 [5] M Brnadou, PG Ciarlt, B Miara, Existnc torms for two-dimnsional linar sll toris, J Elasticity, 34 (1994), pp [6] SC Brnnr, Korn s inqualitis for picwis H 1 vctor filds, Mat Comp, 73(2003) pp [7] A Buffa, C Ortnr, Compact mbddings of brokn Sobolv spacs and applications, IMA J Numr Anal, 29(2009) pp [8] D Capll, KJ Bat, T finit lmnt analysis of slls Fundamntals, Springr, 2011 [9] PG Ciarlt, T finit lmnt mtod for lliptic problms, Nort-Holland, 1978 [10] D A Di Pitro, A Ern, Discrt functional analysis tools for discontinuous Galrkin mtods wit application to t incomprssibl Navir Stoks quations, Mat Comp, 79(2010) pp [11] K Fng, On t tory of discontinuous finit lmnts, Mat Numr Sinica, 4 (1979), pp [12] V Girault, P-A Raviart, Finit lmnt mtods for Navir Stoks quations, tory and algoritms, Springr-Vrlag, 1986 [13] WT Koitr, On t foundations of t linar tory of tin lastic slls, Ndrl Akad Wtnsc Proc Sr B, 73(1970), pp [14] G Duvaut, JL Lions, Inqualitis in mcanics and pysics, Springr-Vrlag, 1976 [15] PM Nagdi, T tory of slls and plats, Handbuc dr Pysik Vol VIa/2, Springr-Vrlag, Brlin, 1972, pp [16] S Zang, Analysis of a discontinuous Galrkin mtod for Koitr sll, arxiv: [17] S Zang, A discontinuous Galrkin mtod for t Nagdi sll modl, arxiv: