Skip to main content

Cut Bogner-Fox-Schmit elements for plates

Abstract

We present and analyze a method for thin plates based on cut Bogner-Fox-Schmit elements, which are \(C^1\) elements obtained by taking tensor products of Hermite splines. The formulation is based on Nitsche’s method for weak enforcement of essential boundary conditions together with addition of certain stabilization terms that enable us to establish coercivity and stability of the resulting system of linear equations. We also take geometric approximation of the boundary into account and we focus our presentation on the simply supported boundary conditions which is the most sensitive case for geometric approximation of the boundary.

Introduction

The Bogner-Fox-Schmit (BFS) element [6] is a classical \(C^1\) thin plate element obtained by taking tensor products of cubic Hermite splines and removing the interior degrees of freedom that are zero on the boundary. In this paper, we consider a variant where we retain these degrees of freedom to obtain a \(C^1\) version of the \(Q_3\) approximation [23]. This element is only \(C^1\) on tensor product (rectangular) elements, which is a serious drawback since it severely limits the applicability of the resulting finite element method. However, on geometries allowing for tensor product discretization it is generally considered to be one of the most efficient elements for plate analysis, cf. [24, p. 153]. It is also a reasonably low order element for plates which is very simple to implement, in contrast with triangular elements which either use higher order polynomials, such as the Argyris element [1], or macro element techniques, such as the Clough–Tocher element [10]. The construction of curved versions of these elements for boundary fitting can also be cumbersome, see, e.g., [5]. It should be noted that the use of straight line segments for discretizing the boundary is not to be recommended, not only because of accuracy issues but also due to Babuška’s paradox for simply supported plates, see [3].

An alternative to \(C^1\) approximations for Kirchhoff plates is to either use discontinuous Galerkin methods [8, 15, 17, 22], or to use mixed finite elements for the Reissner–Mindlin model with small plate thickness [2, 4, 12]. These \(C^0\) methods alleviate the problem of boundary approximation. In this paper we present an alternative idea where \(C^1\) continuity is retained: we develop a cut finite element version, allowing for discretizing a smooth boundary which may cut through the tensor product mesh in an arbitrary manner. Adding stabilization terms on the faces associated with elements that intersect the boundary, we obtain a stable method with optimal order convergence. We prove a priori error estimates which also take approximation of the boundary into account. The focus of the analysis is on simply supported boundary conditions, the computationally most challenging case.

The paper is organized as follows. In “The Kirchhoff plate” section, we recall the thin plate Kirchhoff model; in “The finite element method” section we formulate the cut finite element method; in “Error estimates” section we present the analysis of the method starting with a sequence of technical results leading up to a Strang Lemma and an estimate of the consistency error and finally a priori error estimates in the energy and \(L^2\) norms. In “Numerics” section, we present some numerical illustrations, and in “Concluding remarks” section some concluding remarks are included.

The Kirchhoff plate

Consider a simply supported thin plate in a domain \(\Omega \subset {\mathbb {R}}^2\) with smooth boundary \(\partial \Omega \). The displacement \(u:\Omega \rightarrow {\mathbb {R}}\) satisfies

$$\begin{aligned} \nabla \cdot (\sigma (\nabla u) \cdot \nabla ) = f \end{aligned}$$
(1)

where \(\sigma (\nabla v)\) is the stress tensor

$$\begin{aligned} \sigma (\nabla v) = \kappa (\epsilon (\nabla v) +\nu (1-\nu )^{-1} (\nabla \cdot (\nabla v) ) I = \kappa ( \nabla \otimes \nabla v + \nu (1-\nu )^{-1} (\Delta v) I)) \end{aligned}$$
(2)

with \(\epsilon (\nabla v)\) the strain tensor

$$\begin{aligned} \epsilon (\nabla v) = ( (\nabla v) \otimes \nabla + \nabla \otimes (\nabla v) )/2 = \nabla \otimes \nabla v \end{aligned}$$
(3)

and \(\kappa \) the parameter

$$\begin{aligned} \kappa = \frac{E t^3}{12(1+\nu )} \end{aligned}$$
(4)

with E the Young’s modulus, \(\nu \) the Poisson’s ratio, and t the plate thickness. Since \(0 \le \nu \le 0.5\) both \(\kappa \) and \(\nu (1-\nu )^{-1}\) are uniformly bounded.

We shall focus our presentation on simply supported boundary conditions

$$\begin{aligned} u = 0 \quad \hbox { on}\ \partial \Omega , \qquad M_{nn}(u) = 0 \quad \hbox { on}\ \partial \Omega \end{aligned}$$
(5)

where the moment tensor M is defined by

$$\begin{aligned} M(u) = \sigma (\nabla u) \end{aligned}$$
(6)

and \(M_{ab} = a\cdot M \cdot b\) for \(a,b \in {\mathbb {R}}^2\). Other conditions, such as clamped boundaries, can be handled using the same techniques as in the following, cf. [15].

The weak form of (1) and (5) takes the form: find \(u \in V = \{ v \in H^2(\Omega ) : v=0 \text { on } \, \partial \Omega \}\) such that

$$\begin{aligned} a(u,v) = l(v) \qquad \forall v \in V \end{aligned}$$
(7)

where

$$\begin{aligned} a(v,w) = (\sigma (\nabla v),\epsilon (\nabla w))_\Omega = \kappa (\nabla \otimes \nabla v,\nabla \otimes \nabla w)_{\Omega }+ \nu (1-\nu )^{-1} (\Delta v,\Delta w)_\Omega ) \end{aligned}$$
(8)

and \(l(v) = (f,v)_\Omega \). The form a is symmetric, continuous, and coercive on V equipped with the \(H^2(\Omega )\) norm and it follows from the Lax-Milgram theorem that there exists a unique solution in V to (7). Furthermore, for smooth boundary and \(f \in L^2\) we have the elliptic regularity

$$\begin{aligned} \Vert u \Vert _{H^4(\Omega )} \lesssim \Vert f \Vert _\Omega \end{aligned}$$
(9)

The finite element method

The mesh and finite element space

We begin by introducing the following notation.

  • Let \({\widetilde{{\mathcal {T}}}}_h\), \(h\in (0,h_0]\), be a family of partitions of \({\mathbb {R}}^2\) into squares with side h. Let \({\widetilde{V}}_h\) be the Bogner-Fox-Schmit space consisting of tensor products of cubic Hermite splines on \({\widetilde{{\mathcal {T}}}}_h\). Note that \({\widetilde{V}}_h|_T = Q_3(T)\), with \(Q_3(T)\) the tensor product \(P_3(I_1) \otimes P_3(I_2)\) of cubic polynomials where \(T=I_1 \times I_2 \subset {\mathbb {R}}^2\).

  • Let \(\rho \) be the signed distance function, positive on the outside and negative on the inside, associated with \(\partial \Omega \) and let \(U_\delta (\partial \Omega ) = \{x\in {\mathbb {R}}^2 : |\rho (x)|< \delta \}\) be the tubular neighborhood of \(\partial \Omega \) of thickness \(2\delta \). Then there is \(\delta _0>0\) such that the closest point mapping \(p:U_{\delta _0}(\partial \Omega ) \rightarrow \partial \Omega \) is a well defined function of the form \(p(x) = x - \rho (x) n(p(x))\), cf. [14, Section 14.6].

  • Let \(\{\Omega _h, h \in (0,h_0]\}\) be a family of approximations of \(\Omega \) such that \(\partial \Omega _h \subset U_{\delta _0}(\partial \Omega )\) is piecewise smooth and

    $$\begin{aligned} \Vert \rho \Vert _{L^\infty (\partial \Omega _h)}\lesssim & {} h^4 \end{aligned}$$
    (10)
    $$\begin{aligned} \Vert n(p) - n_h\Vert _{L^\infty (\partial \Omega _h)}\lesssim & {} h^3 \end{aligned}$$
    (11)

    Furthermore, we assume that for each element T such that \(\partial \Omega _h\) intersects the interior of T, i.e. \(\text {int}(T) \cap \partial \Omega _h \ne \emptyset \), the curve segment \(\partial \Omega _h \cap T\) is smooth and intersect the boundary \(\partial T\) of T in precisely two different points. Let \({\mathcal {X}}_h\) be the set of all points where \(\partial \Omega _h\) is not smooth and note that the number of elements \(|{\mathcal {X}}_h|\) in \({\mathcal {X}}_h\) satisfies \(|{\mathcal {X}}_h| \lesssim h^{-1}\).

  • Let \({\mathcal {T}}_h = \{T \in {\widetilde{{\mathcal {T}}}}_h : T \cap \Omega _h \ne \emptyset \}\) be the active mesh and \(\mathcal {F}_h\) the set of interior faces in \(|\mathcal {T}_h\). Let \({\mathcal {T}}_{h,I}\) be the set of elements such that \(T\subset \Omega \cap \Omega _h\) and let \({\mathcal {F}}_{h,I}\) be the set of interior faces in \({\mathcal {T}}_{h,I}\). Let \({\mathcal {T}}_{h,B} = {\mathcal {T}}_h \setminus {\mathcal {T}}_{h,I}\) and \({\mathcal {F}}_{h,B} = {\mathcal {F}}_h \setminus {\mathcal {F}}_{h,I}\). For simplicity we assume that

    $$\begin{aligned} \Omega \subset O_h := \cup _{T\in {\mathcal {T}}_h} T \end{aligned}$$
    (12)

    We can always satisfy this assumption by enlarging the active mesh \({\mathcal {T}}_h\) if necessary.

  • Let \(V_h\) be the restriction of \({\widetilde{V}}_h\) to \({\mathcal {T}}_h\). Let \({\mathcal {K}}_{h}= {\mathcal {T}}_{h}\cap \Omega _h\) be the intersection of the active elements T with \(\Omega _h\).

We illustrate some of these quantities in Fig. 1. Here, \(\partial \Omega _h\) is indicated by the dotted line; \({\mathcal {T}}_h\) is the union of all elements shown; \({\mathcal {T}}_{h,I}\) consists of the white elements and \({\mathcal {T}}_{h,B}\) the grey elements. \({\mathcal {K}}_{h}\) consists of the elements and cut parts of elements inside \(\partial \Omega _h\) and \({\mathcal {F}}_{h,B}\) consists of all element sides on the grey elements excluding those without neighbouring elements.

Fig. 1
figure 1

Illustration of mesh definitions

The finite element method

The method reads: find \(u_h \in V_h\) such that

$$\begin{aligned} \boxed { A_h(u_h,v) = l_h(v) \qquad \forall v \in V_h } \end{aligned}$$
(13)

The forms are defined by

$$\begin{aligned} A_h(v,w)&= a_h(v,w)+\beta s_h(v,w) \end{aligned}$$
(14)
$$\begin{aligned} a_h(v,w)&= (\sigma (\nabla v), \epsilon (\nabla w ) )_{\Omega _h} +(T(v),w)_{\partial \Omega _h} +(v,T(w))_{\partial \Omega _h} \nonumber \\&\quad + \gamma h^{-3} (v,w)_{{\partial \Omega _h}} \end{aligned}$$
(15)
$$\begin{aligned} s_h(v,w)&= h ( [ \nabla ^2_{n_F} v], [ \nabla ^2_{n_F} w])_{{\mathcal {F}}_{h,B}} + h^3 ( [ \nabla ^3_{n_F} v], [ \nabla ^3_{n_F} w])_{{\mathcal {F}}_{h,B}} \end{aligned}$$
(16)
$$\begin{aligned} l_h(v)&= (f, v)_{\Omega _h} \end{aligned}$$
(17)

where

$$\begin{aligned} T = (M\cdot \nabla )_n + \nabla _t M_{nt} \end{aligned}$$
(18)

with sub-indices n and t indicating scalar product with the normal and tangent of \(\partial \Omega _h\), and \(\beta , \gamma \) are positive parameters which are proportional to \(\kappa \). Here \(s_h\) is a stabilization form involving jumps of the second and third normal derivatives in the direction \(n_F\), the normal to the element face, with

$$\begin{aligned}{}[v] := \lim _{\epsilon \downarrow 0} \left( v (x+\epsilon n_F)-v (x-\epsilon n_F))\right) \end{aligned}$$
(19)

which provides necessary stability at the cut elements, see (23). The bilinear form, apart from the stabilization terms, stems from Nitsche’s method [20], first analyzed for plates in a discontinuous Galerkin setting in [15].

Error estimates

Basic properties of \(\varvec{A}_{\varvec{h}}\)

The energy norm

Define the following energy norm on \(V+V_h\), with \(V = H^4(O_h)\),

$$\begin{aligned} |||v |||_h^2&=|||v |||^2_{\Omega _h} + \beta \Vert v \Vert ^2_{s_h} + h^3 \Vert T(v ) \Vert ^2_{ \partial \Omega _h} + h^{-3}\Vert v \Vert ^2_{{\partial \Omega _h}} \end{aligned}$$
(20)

where

$$\begin{aligned} |||v |||^2_{\Omega _h} = {\left( \sigma (\nabla v),\epsilon (\nabla v)\right) _{\Omega _h} } \end{aligned}$$
(21)

and we employ the standard notation \(\Vert v \Vert ^2_{s_h} = s_h(v,v)\). In view of (8) we have \(\kappa \Vert \nabla ^2 v \Vert ^2_{\Omega _h} \lesssim |||v |||^2_{\Omega _h}\), where \(\nabla ^j v\) is the tensor of all j:th order derivatives.

Stabilization

The stabilization term provides us with the following bound

$$\begin{aligned} \boxed { \Vert \nabla ^j v \Vert ^2_{{\mathcal {T}}_{h}} \lesssim \Vert \nabla ^j v \Vert _{{\mathcal {T}}_{h,I}} + h^{2(2-j)} \Vert v \Vert ^2_{s_h}, \qquad j=0,1,2,3 } \end{aligned}$$
(22)

which follows from the standard estimate

$$\begin{aligned} \Vert \nabla ^j v \Vert ^2_{T_2 }\lesssim {\Vert \nabla ^jv \Vert ^2_{T_1 } } + \sum _{k=j}^p h^{2(k-j)} \Vert [ \nabla ^k v ] \Vert ^2_{F} \end{aligned}$$
(23)

where \(T_1\) and \(T_2\) are elements that share the face F, and \(v|_{T_i} \in P_p(T_i)\), the space of polynomials of order p. See for instance [16, 19] for further details.

Continuity and coercivity

The form \(A_h\) is continuous

$$\begin{aligned} A_h(v,w) \lesssim |||v |||_h |||w |||_h \qquad v,w \in V + V_h \end{aligned}$$
(24)

which follows directly from the Cauchy-Schwarz inequality, and for \(\gamma \) large enough coercive

$$\begin{aligned} |||v |||_h^2 \lesssim A_h(v,v) \qquad v \in V_h \end{aligned}$$
(25)

Verification of (25)

We first recall the cut trace inequality

$$\begin{aligned} \Vert v\Vert ^2_{\partial \Omega _h \cap T} \lesssim h^{-1} \Vert v\Vert ^2_{T} + h \Vert \nabla v\Vert ^2_{T} \qquad v \in H^1(T) \end{aligned}$$
(26)

see [25] for a derivation. For \(v \in V_h |_T = V_h(T)\) we have the standard inverse inequality

$$\begin{aligned} \Vert \nabla ^k v \Vert _{T} \lesssim h^{l-k} \Vert \nabla ^l v \Vert _{T} \qquad v \in V_h(T), \qquad k \ge l \end{aligned}$$
(27)

which combined with (26) give the cut inverse trace inequality

$$\begin{aligned} \Vert v\Vert ^2_{\partial \Omega _h \cap T} \lesssim h^{-1} \Vert v\Vert ^2_{T} \qquad V_h(T) \end{aligned}$$
(28)

Now, using the inverse trace inequality (28), the inverse inequality (27), followed by the stabilization estimate (22) we obtain

$$\begin{aligned} \kappa ^{-1} h^3 \Vert T(v) \Vert ^2_{\partial \Omega _h}&\lesssim \kappa h^2 \Vert \nabla ^3 v \Vert ^2_{{\mathcal {T}}_{h}(\partial \Omega _h)} \lesssim \kappa \Vert \nabla ^2 v \Vert ^2_{{\mathcal {T}}_{h}(\partial \Omega _h)} \end{aligned}$$
(29)
$$\begin{aligned}&\lesssim \kappa (\Vert \nabla ^2 v \Vert ^2_{\Omega _h} + \Vert v \Vert ^2_{s_h}) \lesssim |||v |||^2_{\Omega _h} + \kappa \Vert v \Vert ^2_{s_h} \end{aligned}$$
(30)

and thus there is a constant \(C_*\) such that

$$\begin{aligned} \kappa ^{-1} h^3 \Vert T(v) \Vert ^2_{\partial \Omega _h} \le C_* ( |||v |||^2_{\Omega _h} + \kappa \Vert v \Vert ^2_{s_h} ) \end{aligned}$$
(31)

As a consequence, \(|||v |||_h\) and \(|||v |||_{h,*}\), where

$$\begin{aligned} |||v |||_{h,*}^2&=|||v |||^2_{\Omega _h} + \beta \Vert v \Vert ^2_{s_h} + h^{-3}\Vert v \Vert ^2_{{\partial \Omega _h}} \end{aligned}$$
(32)

are equivalent norms on \(V_h\). We then have

$$\begin{aligned} A_h(v,v)&= |||v |||^2_{\Omega _h} + \beta \Vert v \Vert ^2_{s_h} {+ 2(T(v),v)_{\partial \Omega _h}} + \gamma h^{-3} \Vert v \Vert ^2_{\partial \Omega _h} \end{aligned}$$
(33)
$$\begin{aligned}&\ge |||v |||^2_{\Omega _h} + \beta \Vert v \Vert ^2_{s_h} - \delta \kappa ^{-1}h^3 \Vert T(v)\Vert ^2_{\partial \Omega _h} + (\gamma - \delta ^{-1} \kappa ) h^{-3} \Vert v\Vert ^2_{\partial \Omega _h} \end{aligned}$$
(34)
$$\begin{aligned}&\ge (1 - C_* \delta ) |||v |||^2_{\Omega _h} + (\beta - \kappa C_* \delta ) \Vert v \Vert ^2_{s_h} + (\gamma - \delta ^{-1} \kappa )h^{-3} \Vert v\Vert ^2_{\partial \Omega _h} \end{aligned}$$
(35)

and we find that taking \(\delta \) small enough to guarantee that \(1 - C_* \delta \ge m>0\), \(\beta \) large enough to guarantee that \(\beta - \kappa C_* \delta \ge m\), and \(\gamma \) large enough to guarantee that \(\gamma -\delta ^{-1}\kappa \ge m\) leads to \(A_h(v,v) > rsim |||v |||_{{h,*}}^2 > rsim |||v |||_{h}^2\).

Poincaré inequality

We have the following Poincaré inequality

$$\begin{aligned} \boxed { \Vert v \Vert _{H^2({\mathcal {T}}_h)} \lesssim |||v |||_h \qquad v \in V_h } \end{aligned}$$
(36)

Verification of (36)

Using the stabilization estimate (22) and the fact that \({\mathcal {T}}_{h,I}\) is covered by \(\Omega _h\) we have

$$\begin{aligned} \Vert \nabla ^2 v\Vert ^2_{O_h} \lesssim \Vert \nabla ^2 v \Vert ^2_{\Omega _h} + \Vert v \Vert ^2_{s_h} \lesssim |||v |||_h^2 \end{aligned}$$
(37)

Next to estimate \(\Vert v \Vert ^2_{O_h}\) we again use the stabilization estimate (22) and the fact that \({\mathcal {T}}_{h,I}\) is covered by \(\Omega \),

$$\begin{aligned} \Vert v \Vert ^2_{H^1(O_h)} \lesssim \Vert v \Vert ^2_{H^1(\Omega )} + \Vert v \Vert ^2_{s_h} \end{aligned}$$
(38)

Let \(P_{1,\Omega }:L^2(\Omega ) \rightarrow P_1(\Omega )\) be the \(L^2\) projection onto the space of linear functions on \(\Omega \). Then for \(v \in H^2(\Omega )\), and in particular for \(v \in V_h \subset H^2(\Omega )\), we have the Poincaré estimate

$$\begin{aligned} \Vert v - P_{1,\Omega } v \Vert _{H^1(\Omega )} \lesssim \Vert \nabla ^2 v \Vert _\Omega \end{aligned}$$
(39)

and using the trace inequality \(\Vert w \Vert _{\partial \Omega } \lesssim \Vert w \Vert _{H^1(\Omega )}\) with \(w = v - P_{1,\Omega } v\), we obtain

$$\begin{aligned} \Vert v - P_{1,\Omega } v \Vert _{\partial \Omega } \lesssim \Vert v - P_{1,\Omega } v \Vert _{H^1( \Omega )} \lesssim \Vert \nabla ^2 v \Vert _\Omega \end{aligned}$$
(40)

Note that the constants are independent of the mesh parameter since \(\Omega \) is fixed. We then have

$$\begin{aligned} \Vert v \Vert ^2_{H^1(\Omega )}&\lesssim \Vert (I- P_{1,\Omega } )v \Vert ^2_{H^1(\Omega )} + \Vert P_{1,\Omega } v \Vert ^2_{H^1(\Omega )} \end{aligned}$$
(41)
$$\begin{aligned}&\lesssim \Vert \nabla ^2 v \Vert ^2_\Omega + \Vert P_{1,\Omega } v \Vert ^2_{\partial \Omega _h} \end{aligned}$$
(42)
$$\begin{aligned}&\lesssim \underbrace{\Vert \nabla ^2 v \Vert ^2_{\Omega _h} + \Vert v \Vert ^2_{s_h} + \Vert v \Vert ^2_{\partial \Omega _h}}_{i \lesssim |||v |||_h^2} + \underbrace{\Vert (I- P_{1,\Omega }) v \Vert ^2_{\partial \Omega _h}}_{ii \lesssim |||v |||_{h_2}} \end{aligned}$$
(43)

which together with (37) proves (36). Here we estimated term i using the stabilization (22),

$$\begin{aligned} \Vert \nabla ^2 v \Vert ^2_\Omega \lesssim \Vert \nabla ^2 v \Vert ^2_{O_h} \lesssim \Vert \nabla ^2 v \Vert ^2_{\Omega _h} + \Vert v \Vert ^2_{s_h} \lesssim |||v |||_h^2 \end{aligned}$$
(44)

Finally, to estimate ii we recall the following technical estimate, which we prove in the “Appendix” of this paper see also the appendix of [9],

$$\begin{aligned} \Vert w \Vert ^2_{\partial \Omega _h} \lesssim \Vert w \Vert ^2_{\partial \Omega } + \delta ^{1/2} \Vert \nabla w \Vert ^2_{U_\delta (\partial \Omega )\cap O_h} \qquad w \in H^1(O_h) \end{aligned}$$
(45)

where \(\delta \sim h^4\), bounds the distance between \(\partial \Omega \) and \(\partial \Omega _h\), see (10). Note that the finite element functions are defined on \(O_h\), that contains both \(\Omega \) and \(\Omega _h\), and \(V_h \subset H^1(O_h)\). Setting \(w= (I- P_{1,h}) v\), using (40), the stabilization estimate (22), we obtain

$$\begin{aligned} \Vert (I- P_{1,h}) v \Vert ^2_{\partial \Omega _h}&\lesssim \Vert (I- P_{1,h}) v \Vert ^2_{\partial \Omega } + \delta ^{1/2} \Vert \nabla (I- P_{1,h}) v \Vert ^2_{U_\delta (\partial \Omega ) \cap O_h} \end{aligned}$$
(46)
$$\begin{aligned}&\lesssim \Vert \nabla ^2 v \Vert ^2_{\Omega } + h^2 \Vert \nabla (I- P_{1,h}) v \Vert ^2_{O_h} \end{aligned}$$
(47)
$$\begin{aligned}&\lesssim \Vert \nabla ^2 v \Vert ^2_{\Omega } + h^2 (\Vert \nabla (I- P_{1,h}) v \Vert ^2_{\Omega } + \Vert (I- P_{1,h}) v \Vert ^2_{s_h} ) \end{aligned}$$
(48)
$$\begin{aligned}&\lesssim \Vert \nabla ^2 v \Vert ^2_{\Omega } + h^2 (\Vert \nabla ^2 v \Vert ^2_{\Omega } + \Vert v \Vert ^2_{s_h} ) \end{aligned}$$
(49)
$$\begin{aligned}&\lesssim \Vert \nabla ^2 v \Vert ^2_{\Omega _h} + h^2 \Vert v \Vert ^2_{s_h} \end{aligned}$$
(50)
$$\begin{aligned}&\lesssim |||v |||_h^2 \end{aligned}$$
(51)

where we finally used (44) to estimate \(\Vert \nabla ^2 v \Vert _\Omega \) and the fact that \(\Vert (I - P_{1,\Omega }) v \Vert _{s_h} = \Vert v \Vert _{s_h}\). This completes the verification. \(\square \)

Interpolation

Let \(I_{h}:C^1({\mathbb {R}}^2) \rightarrow V_h\) be the standard element wise interpolant associated with the degrees of freedom in \(V_h\). Then we have the estimate

$$\begin{aligned} \Vert v - I_{h} v \Vert _{H^m(T)} \lesssim h^{4-m} \Vert v \Vert _{H^4(T)}\qquad m=0,1,2,3 \end{aligned}$$
(52)

To construct an interpolation operator for cut elements we recall that given \(v \in H^s(\Omega )\) there is an extension operator \(E:H^s(\Omega ) \rightarrow H^s({\mathbb {R}}^2)\) such that

$$\begin{aligned} \Vert E v \Vert _{H^s({\mathbb {R}}^2)} \lesssim \Vert v \Vert _{H^s(\Omega )} \end{aligned}$$
(53)

for all \(s>0\), cf. [21]. For simplicity we will often use the notation \(w = E w\) for the extension of \(w \in H^s(\Omega )\) to \({\mathbb {R}}^2\).

We define the interpolation operator

$$\begin{aligned} C^1(\Omega ) \ni v \mapsto I_{h} (E v ) = \pi _h v \in V_{h} \end{aligned}$$
(54)

Combining (52) with (53) we obtain the interpolation error estimate

$$\begin{aligned} \boxed {\Vert v - \pi _h v \Vert _{H^m({\mathcal {T}}_{h})} \lesssim h^{4-m} \Vert v \Vert _{H^4(\Omega )}\qquad m=0,1,2,3} \end{aligned}$$
(55)

For the energy norm we have the estimate

$$\begin{aligned} \boxed {|||v - \pi _h v |||_h \lesssim h^2 \Vert v \Vert _{H^4(\Omega )}} \end{aligned}$$
(56)

Verification of (56)

Let \(\eta = v - \pi _h v\) and recall that

$$\begin{aligned} |||\eta |||^2_h = |||\eta |||^2_{\Omega _h} + \Vert \eta \Vert ^2_{s_h} + h^3\Vert T(\eta ) \Vert ^2_{\partial \Omega _h} + h^{-3} \Vert \eta \Vert ^2_{\partial \Omega _h} \end{aligned}$$
(57)

The first term is directly estimated using the interpolation error estimate (55),

$$\begin{aligned} {\Vert \eta \Vert ^2_{\Omega _h}} \lesssim h^4 \Vert v \Vert ^2_{H^4(\Omega )} \end{aligned}$$
(58)

For the second term we employ the trace inequality

$$\begin{aligned} \Vert w \Vert ^2_{\partial T} \lesssim h^{-1} \Vert w \Vert ^2_T + h \Vert \nabla w \Vert ^2_T \end{aligned}$$
(59)

to conclude that

$$\begin{aligned} \Vert \eta \Vert ^2_{s_h}&= \sum _{j=2}^3 h^{2j - 3} \Vert [ \nabla _{n_F}^j \eta ] \Vert ^2_{{\mathcal {F}}_{h,B}} \nonumber \\&\lesssim \sum _{j=2}^3 h^{2j - 3} (h^{-1} \Vert \nabla _{n_F}^j \eta \Vert ^2_{{\mathcal {T}}_h({\mathcal {F}}_{h,B})} + h \Vert \nabla _{n_F}^{j+1} \eta \Vert ^2_{{\mathcal {T}}_h({\mathcal {F}}_{h,B})} \end{aligned}$$
(60)
$$\begin{aligned}&\lesssim \sum _{j=2}^3 h^{2j - 4} (\Vert \nabla ^j \eta \Vert ^2_{{\mathcal {T}}_{h}} + h^2 \Vert \nabla ^{j+1} \eta \Vert ^2_{{\mathcal {T}}_{h}} )\nonumber \\&\lesssim \sum _{j=2}^3 h^{2j - 4} h^{2(4-j)} \Vert v \Vert ^2_{H^4(\Omega )} \lesssim h^4 \Vert v \Vert ^2_{H^4(\Omega )} \end{aligned}$$
(61)

For the third term we use the cut trace inequality (26) and the interpolation estimate (55),

$$\begin{aligned} h^3 \Vert T(\eta )\Vert ^2_{\partial \Omega _h}&\lesssim h^3 ( h^{-1} \Vert \nabla ^3 \eta \Vert ^2_{{\mathcal {T}}_h(\partial \Omega _h)} + h \Vert \nabla ^4 \eta \Vert ^2_{{\mathcal {T}}_{h}(\partial \Omega _h)} ) \end{aligned}$$
(62)
$$\begin{aligned}&\lesssim h^2 \Vert \nabla ^3 \eta \Vert ^2_{{\mathcal {T}}_h} + h^4 \Vert \nabla ^4 \eta \Vert ^2_{{\mathcal {T}}_{h}} \lesssim h^4 \Vert v \Vert ^2_{H^4(\Omega )} \end{aligned}$$
(63)

where \({\mathcal {T}}_h({\mathcal {F}}_{h,B}) \subset {\mathcal {T}}_h\) is the set of elements with a face in \({\mathcal {F}}_{h,B}\). Finally, the fourth term is estimated in the same way as the third,

$$\begin{aligned} h^{-3}\Vert \eta \Vert ^2_{\partial \Omega _h}&\lesssim h^{-3} (h^{-1} \Vert \eta \Vert ^2_{{\mathcal {T}}_{h}(\partial \Omega _h)} + h \Vert \nabla \eta \Vert ^2_{{\mathcal {T}}_{h}(\partial \Omega _h)} ) \end{aligned}$$
(64)
$$\begin{aligned}&\lesssim h^{-4} \Vert \eta \Vert ^2_{{\mathcal {T}}_h} + h^{-2} \Vert \nabla \eta \Vert ^2_{{\mathcal {T}}_{h}} \lesssim h^4 \Vert v \Vert ^2_{H^4(\Omega )} \end{aligned}$$
(65)

which completes the verification of (56). \(\square \)

Consistency error estimate

Lemma 1

Let u be the exact solution to (1) with boundary conditions (5), and \(u_h\) the finite element approximation defined by (13), then

$$\begin{aligned} \boxed { |||u - u_h |||_{h} \lesssim |||u - \pi _h u |||_h + \sup _{v \in V_h\setminus \{0\}} {\frac{A_h(\pi _hu,v) - l_h(v)}{|||v |||_h}} } \end{aligned}$$
(66)

Proof

Adding and subtracting an interpolant we obtain

$$\begin{aligned} |||u - u_h |||_{h} \le |||u - \pi _h u |||_h + |||\pi _h u - u_h |||_h \end{aligned}$$
(67)

Using coercivity we can estimate the second term on the right hand side as follows

$$\begin{aligned} |||\pi _h u - u_h |||_{h}&\le \sup _{v \in V_h\setminus \{0\}} \frac{A_h(\pi _h u - u_h,v) }{|||v |||_h} \end{aligned}$$
(68)
$$\begin{aligned}&\le \sup _{v \in V_h\setminus \{0\}} \frac{A_h(\pi _h u - u,v) }{|||v |||_h} + \sup _{v \in V_h\setminus \{0\}} \frac{A_h(\pi _h u - u_h,v) }{|||v |||_h} \end{aligned}$$
(69)
$$\begin{aligned}&\le |||\pi _h u - u |||_h + \sup _{v \in V_h\setminus \{0\}} \frac{A_h(\pi _h u,v) - l_h(v)}{|||v |||_h} \end{aligned}$$
(70)

where we added and subtracted u in the numerator and for the first term used the estimate \(A_h(\pi _h u - u,v) \lesssim |||\pi _h u - u |||_h |||v |||_h\) and for the second used (13) to eliminate \(u_h\). Combining the estimates the desired result follows directly. \(\square \)

Lemma 2

Let \(\varphi \in H^4({\mathbb {R}}^2)\) and \(v \in V+V_h\), then

$$\begin{aligned} (\nabla \cdot (M(\varphi ) \cdot \nabla ),v )_{\Omega _h}&= (M(\varphi ) , \epsilon (\nabla v) )_{\Omega _h} - ( M_{nn}(\varphi ), \nabla _n v)_{\partial \Omega _h} \end{aligned}$$
(71)
$$\begin{aligned}&\quad + (T(\varphi ),v)_{\partial \Omega _h} + ([M_{nt}],v)_{{\mathcal {X}}_h} \end{aligned}$$
(72)

where, for \(x\in {\mathcal {X}}_h\) (the set of points where \(\partial \Omega _h\) is not smooth), \([M_{nt}]_x\) is defined by

$$\begin{aligned}{}[M_{nt}] |_x = M(x)_{n_h^+ t_h^+} - M(x)_{n_h^- t_h^-} \end{aligned}$$
(73)

with \(n_h^\pm \) and \(t_h^\pm \) the left and right limits to tangent and normal to the discrete boundary \(\partial \Omega _h\) at \(x \in {\mathcal {X}}_h\). In the case of \(C^1\) boundary \((v, [M_{nt}])_{{\mathcal {X}}_h}=0\).

Proof

Using the simplified notation \(M = M(\varphi )\) and \(T=T(\varphi )\) for brevity we obtain by integrating by parts

$$\begin{aligned} (\nabla \cdot (M\cdot \nabla ),v)_{\Omega _h}&=( (M \cdot \nabla )_n,v )_{\partial \Omega _h} - (M \cdot \nabla , \nabla v )_{\Omega _h} \end{aligned}$$
(74)
$$\begin{aligned}&= ((M \cdot \nabla )_n, v )_{\partial \Omega _h} - (M_n, \nabla v )_{\partial \Omega _h} + (M, \epsilon (\nabla v))_{\Omega _h} \end{aligned}$$
(75)

Splitting \(\nabla v\) in tangent and normal contributions on \(\partial \Omega _h\), we have the identity

$$\begin{aligned} (\nabla v, M_n )_{\partial \Omega _h \cap T}&= (\nabla _n v, M_{nn} )_{\partial \Omega _h \cap T} + (\nabla _t v, M_{nt} )_{\partial \Omega _h \cap T} \end{aligned}$$
(76)
$$\begin{aligned}&= (\nabla _n v, M_{nn} )_{\partial \Omega _h\cap T} - (v, \nabla _t M_{nt} )_{\partial \Omega _h \cap T} + (v, M_{nt} t \cdot \nu )_{ \partial (\partial \Omega _h \cap T )} \end{aligned}$$
(77)

where we integrated by parts along the curve segments \(\partial \Omega _h \cap T\), and \(\nu \) is the exterior unit tangent vector to \(\partial \Omega _h \cap T\). Summing over all elements that intersect \(\partial \Omega _h\), we obtain the identity

$$\begin{aligned} (\nabla v, M_n )_{\partial \Omega _h}&= (\nabla _n v, M_{nn} )_{\partial \Omega _h} - (v, \nabla _t M_{nt} )_{\partial \Omega _h} + (v, [M_{nt}])_{{\mathcal {X}}_h} \end{aligned}$$
(78)

Combining (75) and (78), we obtain

$$\begin{aligned} (v,\nabla \cdot (M\cdot \nabla ) )_{\Omega _h}&= (\epsilon (\nabla v) ,M)_{\Omega _h} - (\nabla _n v, M_{nn} )_{\partial \Omega _h} \end{aligned}$$
(79)
$$\begin{aligned}&\quad + (v, (M \cdot \nabla )_n +\nabla _t M_{nt} )_{\partial \Omega _h} - (v, [M_{nt}])_{{\mathcal {X}}_h} \end{aligned}$$
(80)

and setting \(T = (M \cdot \nabla )_n +\nabla _t M_{nt}\) we obtain the desired result. \(\square \)

Lemma 3

Let u be the exact solution to (1) with boundary conditions (5), then there is a constant such that for all \(v\in V_h\),

$$\begin{aligned} \boxed {A_h(u ,v) - l_h(v) \lesssim h^4 \Vert u \Vert _{H^4(\Omega )} |||v |||_{h,\bigstar } \lesssim h^{5/2} \Vert u \Vert _{H^4(\Omega )} |||v |||_{h}} \end{aligned}$$
(81)

where \(|||v |||_{h,\bigstar }\) is the norm

$$\begin{aligned} |||v |||^2_{h,\bigstar } = |||v |||_h^2 + \Vert T(v) \Vert ^2_{\partial \Omega _h} + h^{-6} \Vert v \Vert ^2_{\partial \Omega _h} \le (1+ h^{-3}) |||v |||^2_h \end{aligned}$$
(82)

Proof

Using the definition (14), the fact that \(s_h(u,v) = 0\) for \(u \in H^4(\Omega )\), and the partial integration identity (71) we obtain

$$\begin{aligned} A_h(u,v)- l_h(v)&= (M(u),\epsilon (\nabla v))_{\Omega _h} + (T(u),v)_{\partial \Omega _h} + (u,T(v))_{\partial \Omega _h} \end{aligned}$$
(83)
$$\begin{aligned}&\quad + \gamma h^{-3}(u,v)_{\partial \Omega _h} - (\nabla \cdot ( M(u) \cdot \nabla ) ,v)_{\Omega _h} \end{aligned}$$
(84)
$$\begin{aligned}&= (M_{nn}(u), \nabla _n v )_{\partial \Omega _h} + ([M_{nt}],v)_{{\mathcal {X}}_h} \end{aligned}$$
(85)
$$\begin{aligned}&\quad +(u, T(v))_{\partial \Omega _h} + \gamma h^{-3} (u,v)_{\partial \Omega _h} \end{aligned}$$
(86)
$$\begin{aligned}&=I + II + III + IV \end{aligned}$$
(87)

Before turning to the estimates of \(I-IV\) we first note that for \(w\in H^1_0(\Omega )\), with extension to \({\mathbb {R}}^2\) also denoted by w, we may apply (45) and the stability (53) of the extension that

$$\begin{aligned} \Vert w \Vert ^2_{\partial \Omega _h} \lesssim \delta \Vert w \Vert ^2_{H^1(U_\delta (\partial \Omega ))} \lesssim \delta \Vert w \Vert ^2_{H^1(U_{\delta _0}(\partial \Omega )\cup \Omega )} \lesssim \delta \Vert w \Vert ^2_{H^1(\Omega )} \end{aligned}$$
(88)

with \(\delta \sim h^4\). Note that here we do not need to restrict the norms to \(O_h\) as in (45) since the extended function is defined on \({\mathbb {R}}^2\). See [9, Appendix] for detailed derivations. Furthermore, for more regular functions such that \(w \in H^2(\Omega )\) with \(w =0\) on \(\partial \Omega \), we may strengthen the estimate as follows

$$\begin{aligned}&\Vert w \Vert ^2_{\partial \Omega _h} \lesssim \delta \Vert \nabla w \Vert ^2_{L^2(U_\delta (\partial \Omega ))}\lesssim \delta ^2 \sup _{t \in (-\delta ,\delta )} \Vert \nabla w \Vert ^2_{\partial \Omega _t} \lesssim \delta ^2 \Vert w \Vert _{H^2(\Omega )} \end{aligned}$$
(89)

where \(\partial \Omega _t = \{ x \in {\mathbb {R}}^2 | \rho (x) = t \}\) for \(t \in (-\delta _0,\delta _0)\), are level sets of \(\rho \), and again \(\delta \sim h^4\). In the last step we used a version of (45) to conclude that

$$\begin{aligned}&\Vert \nabla w \Vert ^2_{\partial \Omega _t} \lesssim \Vert \nabla w \Vert ^2_{\partial \Omega } + t \Vert \nabla ^2 w \Vert ^2_{U_t(\partial \Omega )}\nonumber \\&\quad \lesssim \Vert w \Vert ^2_{H^2(\Omega )} + \delta _0 \Vert \nabla ^2 w \Vert ^2_{U_{\delta _0}(\partial \Omega )} \lesssim \Vert w \Vert ^2_{H^2(\Omega )} \end{aligned}$$
(90)

where we used a trace inequality on \(\Omega \) and the stability (53) of the extension operator.

  1. I.

    Using Cauchy-Schwarz followed by (89) with \(w = M_{nn}(u)\), we get

    $$\begin{aligned}&(M_{nn}(u), \nabla _n v )_{\partial \Omega _h} \lesssim \Vert M_{nn}(u)\Vert _{\partial \Omega _h} \Vert \nabla _n v \Vert _{\partial \Omega _h} \lesssim \delta \Vert u \Vert _{H^4(\Omega )} |||v |||_h \end{aligned}$$
    (91)

    Here we used the estimate

    $$\begin{aligned} \Vert \nabla _n v \Vert _{\partial \Omega _h} \lesssim |||v |||_h \end{aligned}$$
    (92)

    which we derive by applying (45) with \(w=\nabla v\),

    $$\begin{aligned} \Vert \nabla _n v \Vert ^2_{\partial \Omega _h}&\lesssim \Vert \nabla v \Vert ^2_{\partial \Omega _h} \lesssim \Vert \nabla v \Vert ^2_{\partial \Omega } + \delta \Vert \nabla ^2 v \Vert ^2_{U_\delta (\partial \Omega ) \cap O_h} \end{aligned}$$
    (93)
    $$\begin{aligned}&\lesssim \Vert v \Vert ^2_{H^2(\Omega )} + \delta \Vert \nabla ^2 v \Vert ^2_{U_\delta (\partial \Omega ) \cap O_h} \lesssim \Vert v \Vert ^2_{O_h} \lesssim |||v |||_h^2 \end{aligned}$$
    (94)

    where we used the trace inequality \(\Vert \nabla v \Vert _{\partial \Omega } \lesssim \Vert \nabla v \Vert _{H^1(\Omega )} \lesssim \Vert v \Vert _{H^2(\Omega )}\), and at last the Poincaré inequality (36).

  2. II.

    Using the assumption on the accuracy of the discrete normal (11) we have for each \(x\in {\mathcal {X}}_h\),

    $$\begin{aligned} |[M_{nt}]| = M^+_{n_h t_h} - M^-_{n_h t_h}= M^+_{n_h t_h} -M_{n t} + M_{n t} - M^-_{n_h t_h} \end{aligned}$$
    (95)

    where the first term on the right hand side can be estimated as follows

    $$\begin{aligned} {|M^+_{n_h t_h} - M_{n t} |}\le |(n_h -n )\cdot M \cdot t_h |+ | n\cdot M \cdot (t_h - t)| \lesssim h^3 |M| \end{aligned}$$
    (96)

    We then have

    $$\begin{aligned} ([M_{nt}],v)_{{\mathcal {X}}_h}&\le \Vert [M_{nt}]\Vert _{{\mathcal {X}}_h} \Vert v\Vert _{{\mathcal {X}}_h} \lesssim h^3 \Vert M \Vert _{{\mathcal {X}}_h} \Vert v\Vert _{{\mathcal {X}}_h}\nonumber \\&\lesssim h^2 h^{1/2}\Vert M \Vert _{{\mathcal {X}}_h} h^{1/2} \Vert v\Vert _{{\mathcal {X}}_h} \end{aligned}$$
    (97)
    $$\begin{aligned}&\lesssim h^2 \Vert M \Vert _{L^\infty ({\mathcal {X}}_h)} h^2 |||v |||_h \lesssim h^4 \Vert u \Vert _{H^4(\Omega )} |||v |||_h \end{aligned}$$
    (98)

    where we used the fact that the number of elements, denoted by \(| {\mathcal {X}}_h |\), in \({\mathcal {X}}_h\) satisfies \(| {\mathcal {X}}_h | \sim h^{-1}\), and the Sobolev inequality [13] followed by the stability (53) of the extension operator to obtain

    $$\begin{aligned} h \Vert M \Vert ^2_{L^\infty ({\mathcal {X}}_h)} \lesssim \Vert u \Vert ^2_{W^2_\infty ({\mathbb {R}}^2)} \lesssim \Vert u \Vert ^2_{H^{3+\epsilon }({\mathbb {R}}^2)} \lesssim \Vert u \Vert ^2_{H^4({\mathbb {R}}^2)} \lesssim \Vert u \Vert ^2_{H^4(\Omega )} \end{aligned}$$
    (99)

    and the estimate

    $$\begin{aligned} h \Vert v\Vert ^2_{{\mathcal {X}}_h} \lesssim h^3 |||v |||^2_h\qquad v \in V_h \end{aligned}$$
    (100)

    To verify (100) we shall employ an inverse inequality locally using a linear approximation of the boundary. To that end consider \(x \in {\mathcal {X}}_h\) and let \(B_r(x)\) be a ball of radius \(r \sim h\) centred at x. Let \(T_x \in {\mathcal {T}}_h\) be one of the elements such that \(x \in \partial T_x \cap \partial \Omega _h\) and given \(v\in V_h\) let \(v_x \in Q_3({\mathbb {R}}^2)\) be the extension to \({\mathbb {R}}^2\) of \(v|_{T_x} \in Q_3(T_x)\) such that \(v_x|_{T_x} = v|_{T_x}\). Let \(\Gamma _x\) be a line through x such that \(\Vert \rho \Vert _{L^\infty (\Gamma _x \cap B_r(x))} \lesssim h^2\). Such a line exists since \(\partial \Omega \) is smooth and linear approximation is of second order locally. We then have

    $$\begin{aligned} h \Vert v \Vert ^2_{{\mathcal {X}}_h} \lesssim \sum _{x \in {\mathcal {X}}_h} h |v(x)|^2 \lesssim \sum _{x \in {\mathcal {X}}_h} \Vert v_x \Vert ^2_{\Gamma _x \cap B_r(x)} \end{aligned}$$
    (101)

    where we used the fact that \(|{\Gamma _x \cap B_r(x)}| \sim h\). In order to employ a local version of (45) we define the cylindrical tubular neighborhood

    $$\begin{aligned} U_{\delta ,x} = U_{\delta _1}(p(\Gamma _x \cap B_r(x))) \qquad \delta \in (0,\delta _0] \end{aligned}$$
    (102)

    over \(p(\Gamma _x \cap B_r(x))\subset \partial \Omega \), where \(p:U_{\delta _0}(\partial \Omega ) \rightarrow \partial \Omega \) is the closest point mapping. Then there is \(\delta _1 \sim h^2\) such that \(\Gamma _x \cap B_r(x) \subset U_{\delta _1,x}\) and \(\partial \Omega _h \cap U_{\delta _0,x} \subset U_{\delta _1,x}\), since \(\Gamma _x\) is \(O(h^2)\) accurate locally and \(\partial \Omega _h\) is \(O(h^4)\) accurate. We obtain

    $$\begin{aligned} \Vert v_x \Vert ^2_{\Gamma _x \cap B_r(x)}&\lesssim \Vert v_x \Vert ^2_{\partial \Omega _h \cap U_{\delta _1,x}} + \delta _1 \Vert \nabla v_x \Vert ^2_{ U_{\delta _1,x}} \end{aligned}$$
    (103)
    $$\begin{aligned}&\lesssim \Vert v_x \Vert ^2_{\partial \Omega _h \cap U_{\delta _1,x}}+ \delta _1 \Vert \nabla v_x \Vert ^2_{T_x} \end{aligned}$$
    (104)
    $$\begin{aligned}&\lesssim \Vert v \Vert ^2_{\partial \Omega _h \cap U_{\delta _1,x}} + \Vert v_x - v \Vert ^2_{\partial \Omega _h \cap B_r(x)} + \delta _1 \Vert \nabla v \Vert ^2_{T_x} \end{aligned}$$
    (105)

    then we added and subtracted v in the first term, used the triangle inequality, the inverse estimate \(\Vert \nabla v_x \Vert ^2_{U_{\delta _1,x}} \lesssim \Vert \nabla v_x \Vert _{T_x}\) which holds since \(v_x\) is a polynomial, and \(U_{\delta _1,x} \subset B_{r'}(x)\) for some ball of radius \(r'\sim h\), and finally we noted that \(v_x = v\) on \(T_x\). Inequality (103) is an application of (154). Inserting (105) in (101) we get

    $$\begin{aligned} h \Vert v \Vert ^2_{{\mathcal {X}}_h}&\lesssim \sum _{x \in {\mathcal {X}}_h} \Vert v \Vert ^2_{\partial \Omega _h \cap U_{\delta _1,x}} + \Vert v_x - v \Vert ^2_{\partial \Omega _h \cap U_{\delta _1,x}} + \delta _1 \Vert \nabla v \Vert ^2_{T_x} \end{aligned}$$
    (106)
    $$\begin{aligned}&\lesssim \Vert v \Vert ^2_{\partial \Omega _h} + \underbrace{\sum _{x \in {\mathcal {X}}_h} \Vert v_x - v \Vert ^2_{\partial \Omega _h \cap B_r(x)}}_{i \lesssim h^4 \Vert v \Vert ^2_{s_h}} + \underbrace{\sum _{x \in {\mathcal {X}}_h} \delta _1 \Vert \nabla v \Vert ^2_{T_x}}_{ii\lesssim h^3 \Vert v \Vert ^2_{H^2(\Omega )} + h^4 \Vert v \Vert ^2_{s_h}} \end{aligned}$$
    (107)
    $$\begin{aligned}&\lesssim \Vert v \Vert ^2_{\partial \Omega _h} + h^3 \Vert v \Vert ^2_{H^2(\Omega )} + h^4 \Vert v \Vert ^2_{s_h} \end{aligned}$$
    (108)
    $$\begin{aligned}&\lesssim h^3 |||v |||_h^2 \end{aligned}$$
    (109)

    which establishes (100). Here we used the fact that the number of cylinders \(U_{\delta _1,y}\), \(y\in {\mathcal {X}}_h\), that intersect \(U_{\delta _1,x}\) is uniformly bounded independent of \(x \in {\mathcal {X}}_h\) and \(h \in (0,h_0]\). We also used certain estimates of terms i and ii, which we verify next.

    1. i.

      Let \({\mathcal {T}}_{h,x} = {\mathcal {T}}_h(\partial \Omega _h \cap B_r(x))\), \({\mathcal {F}}_{h,x}\) be the interior faces in \({\mathcal {T}}_{h,x}\), and \(s_{h,x}\) be defined by (16) with \({\mathcal {F}}_{h,B}\) replaced by \({\mathcal {F}}_{h,x}\). Using the estimate

      $$\begin{aligned} \Vert v_x - v \Vert ^2_{{\mathcal {T}}_{h,x}} \lesssim h^4 \Vert v \Vert ^2_{s_{h,x}} \end{aligned}$$
      (110)

      which is a local version of (22) on the patch \({\mathcal {T}}_{h,x}\), we get

      $$\begin{aligned} i \lesssim \sum _{x\in {\mathcal {X}}_h} \Vert v_x - v \Vert ^2_{{\mathcal {T}}_{h,x}} \lesssim \sum _{x\in {\mathcal {X}}_h} h^4 \Vert v \Vert ^2_{s_{h,x}} \lesssim h^4 \Vert v \Vert ^2_{s_h} \end{aligned}$$
      (111)

      .

    2. ii.

      We first note that

      $$\begin{aligned} ii= \sum _{x \in {\mathcal {X}}_h} \delta _1 \Vert \nabla v \Vert ^2_{T_x}&\lesssim \delta _1 \Vert \nabla v \Vert ^2_{{\mathcal {T}}_h(\partial \Omega _h)}\nonumber \\&\lesssim \delta _1( \Vert \nabla v \Vert ^2_{{\mathcal {N}}_h({\mathcal {T}}_h(\partial \Omega _h)) \setminus {\mathcal {T}}_h(\partial \Omega _h)} + h^2\Vert v \Vert ^2_{s_h} ) \end{aligned}$$
      (112)

      where \({\mathcal {N}}_h({\mathcal {T}}_h(\partial \Omega _h))\) is the set of elements that share a node with an element in \({\mathcal {T}}_h(\partial \Omega _h)\), and we used the stabilization. We can then choose \(\delta _2 \sim h\) such that

      $$\begin{aligned} {\mathcal {N}}_h({\mathcal {T}}_h(\partial \Omega _h)) \setminus {\mathcal {T}}_h(\partial \Omega _h) \subset U_{\delta _2}(\partial \Omega ) \cap \Omega \end{aligned}$$
      (113)

      We now proceed in a similar way as in (89),

      $$\begin{aligned} \Vert \nabla v \Vert ^2_{{\mathcal {N}}_h({\mathcal {T}}_h(\partial \Omega _h)) \setminus {\mathcal {T}}_h(\partial \Omega _h)}&\lesssim \Vert \nabla v \Vert ^2_{U_{\delta _2}(\partial \Omega ) \cap \Omega }\nonumber \\&\lesssim \sup _{t \in (0,-\delta _2)} \delta _2 \Vert \nabla v \Vert ^2_{\partial \Omega _t} \lesssim \delta _2 \Vert v \Vert ^2_{H^2(\Omega )} \end{aligned}$$
      (114)

      where we used the estimate (90) in the last step. Inserting (114) in (112) and using that \(\delta _1 \sim h^2\) and \(\delta _2 \sim h\), we obtain

      $$\begin{aligned} ii \lesssim \delta _1 ( \delta _2 \Vert v \Vert ^2_{H^2(\Omega )} + h^2 \Vert v \Vert ^2_{s_h} ) \lesssim h^3 \Vert v \Vert ^2_{H^2(\Omega )} + h^4 \Vert v \Vert ^2_{s_h} \end{aligned}$$
      (115)
  3. III.

    Using (89) with \(w=u\) and recalling \(\delta \sim h^4\),

    $$\begin{aligned} (u, T(v))_{\partial \Omega _h}&\le \Vert u\Vert _{\partial \Omega _h} \Vert T(v)\Vert _{\partial \Omega _h} \lesssim \delta \Vert u \Vert _{H^4(\Omega )} \Vert T(v)\Vert _{\partial \Omega _h} \end{aligned}$$
    (116)
    $$\begin{aligned}&\lesssim \delta h^{-3/2} \Vert u \Vert _{H^4(\Omega )} h^{3/2} \Vert T(v)\Vert _{\partial \Omega _h} \lesssim h^{5/2} \Vert u \Vert _{H^4(\Omega )}|||v |||_h \end{aligned}$$
    (117)
  4. IV.

    Proceeding in the same way as for Term III,

    $$\begin{aligned} h^{-3}(u,v)_{\partial \Omega _h}&\lesssim h^{-3} \Vert u \Vert _{\partial \Omega _h} \Vert v \Vert _{\partial \Omega _h} \lesssim \delta h^{-3} \Vert u \Vert _{H^4(\Omega )} \Vert v \Vert _{\partial \Omega _h} \end{aligned}$$
    (118)
    $$\begin{aligned}&\lesssim \delta h^{-3/2} \Vert u \Vert _{H^4(\Omega )} |||v |||_h \lesssim h^{5/2} \Vert u \Vert _{H^4(\Omega )} |||v |||_h \end{aligned}$$
    (119)

    Combining the estimates we find that

    $$\begin{aligned} A_h(u,v) - l_h(v)&\lesssim h^4 \Vert u \Vert _{H^4(\Omega )} |||v |||_h + {h^{5/2} \Vert u \Vert _{H^4(\Omega )} |||v |||_h} \end{aligned}$$
    (120)
    $$\begin{aligned}&\lesssim h^{5/2} \Vert u \Vert _{H^4(\Omega )} |||v |||_h \end{aligned}$$
    (121)

    which completes the proof. \(\square \)

Error estimates

Theorem 1

The finite element solution defined by (13) satisfies

$$\begin{aligned} \boxed {|||u - u_h |||_h \lesssim h^2 \Vert u \Vert _{H^4(\Omega )}} \end{aligned}$$
(122)

Proof

Using the second bound of (81) in (66) followed by the interpolation error bound (56) we directly get the desired estimate. \(\square \)

Theorem 2

The finite element solution defined by (13) satisfies

$$\begin{aligned} \boxed {\Vert u - u_h \Vert _{\Omega _h} \lesssim h^4 \Vert u \Vert _{H^4(\Omega )}} \end{aligned}$$
(123)

Proof

Adding and subtracting an interpolant and using the interpolation error estimate (55) we have the estimate

$$\begin{aligned} \Vert u - u_h \Vert _{\Omega _h}&\le \Vert u - \pi _h u \Vert _{\Omega _h} + \Vert \pi _h u - u_h \Vert _{\Omega _h} \end{aligned}$$
(124)
$$\begin{aligned}&\lesssim h^4 \Vert u \Vert _{H^4(\Omega )} + \Vert \pi _h u - u_h \Vert _{\Omega _h} \end{aligned}$$
(125)

In order to estimate \(\Vert \pi _h u - u_h \Vert _{\Omega _h}\) we let \(\phi _h\in V_h\) be the solution to the discrete dual problem

$$\begin{aligned} (v,\psi )_{\Omega _h}&= A_h(v,\phi _h) \qquad \forall v \in V_h \end{aligned}$$
(126)

Setting \(v = \pi _h u - u_h\) we obtain the error representation

$$\begin{aligned} (\pi _h u - u_h,\psi )_{\Omega _h}&= A_h( \pi _h u - u_h, \phi _h) \end{aligned}$$
(127)
$$\begin{aligned}&= A_h( \pi _h u - u, \phi _h) + A_h( u - u_h, \phi _h) \end{aligned}$$
(128)
$$\begin{aligned}&= \underbrace{A_h( \pi _h u - u, \phi _h - \phi )}_{I} + \underbrace{A_h( \pi _h u - u,\phi )}_{II} \nonumber \\&\quad + \underbrace{A_h( u,\phi _h) - l_h(\phi _h)}_{III} \end{aligned}$$
(129)

Here \(\phi \) is the solution to the continuous dual problem

$$\begin{aligned} \nabla \cdot (\sigma (\nabla \phi ) \cdot \nabla ) = \psi \quad \hbox { in}\ \Omega , \qquad \phi = M_{nn}(\phi ) = 0 \quad {\hbox { on}\ \partial \Omega } \end{aligned}$$
(130)

extended to \({\mathbb {R}}^2\) using the stable extension operator, see (53).

  1. I.

    Since \(\phi _h\) is the finite element approximation of \(\phi \) we have the error estimate

    $$\begin{aligned} |||\phi - \phi _h |||_h \lesssim h^2 \Vert \phi \Vert _{H^4(\Omega )} \lesssim h^2 \Vert \psi \Vert _{\Omega } \end{aligned}$$
    (131)

    where we used elliptic regularity (9), which directly gives

    $$\begin{aligned} A_h( \pi _h u - u, \phi _h - \phi )&\le |||\pi _h u - u |||_h |||\phi _h - \phi |||_h \lesssim h^4 \Vert u \Vert _{H^4(\Omega )} \Vert \psi \Vert _{\Omega } \end{aligned}$$
    (132)
  2. II.

    Using the fact that \(s_h(\pi _h u - u,\phi ) = 0\) since \(\phi \in H^4(\Omega )\), the partial integration formula (71), the Cauchy-Schwarz inequality, and the interpolation error estimates we obtain

    $$\begin{aligned} A_h(\pi _h u - u,\phi )&= (\pi _h u - u, \psi )_{\Omega _h} -(\nabla _n ( \pi _h u - u),M_{nn}(\phi ))_{\partial \Omega _h} \end{aligned}$$
    (133)
    $$\begin{aligned}&\quad + (T(\pi _h u - u), \phi )_{\partial \Omega _h} + \gamma h^{-3} (\pi _h u - u, \phi )_{\partial \Omega _h} \end{aligned}$$
    (134)
    $$\begin{aligned}&\le \Vert \pi _h u - u\Vert _{\Omega _h} \Vert \psi \Vert _{\Omega _h} +\Vert \nabla _n ( \pi _h u - u)\Vert _{\partial \Omega _h} \Vert M_{nn}(\phi )\Vert _{\partial \Omega _h} \end{aligned}$$
    (135)
    $$\begin{aligned}&\quad + \Vert T(\pi _h u - u) \Vert _{\partial \Omega _h}\Vert \phi \Vert _{\partial \Omega _h} + \gamma h^{-3} \Vert \pi _h u - u\Vert _{\partial \Omega _h} \Vert \phi \Vert _{\partial \Omega _h}\nonumber \\ \end{aligned}$$
    (136)
    $$\begin{aligned}&\lesssim h^4 \Vert u\Vert _{H^4(\Omega _h)} \Vert \psi \Vert _{\Omega _h} + h^{5/2} \Vert u\Vert _{H^4(\Omega _h)} \Vert M_{nn}(\phi )\Vert _{\partial \Omega _h}\nonumber \\ \end{aligned}$$
    (137)
    $$\begin{aligned}&\quad + h^{1/2} \Vert u\Vert _{H^4(\Omega _h)} \Vert \phi \Vert _{\partial \Omega _h} + \gamma h^{-3} h^{7/2} \Vert u\Vert _{H^4(\Omega _h)} \Vert \phi \Vert _{\partial \Omega _h} \end{aligned}$$
    (138)
    $$\begin{aligned}&\lesssim \underbrace{(h^4 + h^{5/2} h^2 + h^{1/2} h^4)}_{\lesssim h^4} \Vert u\Vert _{H^4(\Omega _h)} \Vert \psi \Vert _{\Omega _h} \end{aligned}$$
    (139)

    Here we used (88) with \(\delta \sim h^4\) followed by the elliptic regularity (9) to conclude that

    $$\begin{aligned} \Vert M_{nn}(\phi ) \Vert _{\partial \Omega _h} \lesssim \delta ^{1/2} \Vert \phi \Vert _{H^3(U_\delta (\partial \Omega ))} \lesssim h^2\Vert \phi \Vert _{H^4(\Omega )}\lesssim h^2\Vert \psi \Vert _{\Omega } \end{aligned}$$
    (140)

    and using (89) we obtain

    $$\begin{aligned} \Vert \phi \Vert _{\partial \Omega _h} \lesssim \delta \Vert \phi \Vert _{H^4(\Omega )} \lesssim h^4 \Vert \psi \Vert _{\Omega _h} \end{aligned}$$
    (141)
  3. III.

    Using (81) we obtain

    $$\begin{aligned} |A_h( u,\phi _h) - (f, \phi _h)_{\Omega _h}|&\lesssim h^4 \Vert u \Vert _{H^4(\Omega )} |||\phi _h |||_{h,\bigstar } \lesssim h^4 \Vert u \Vert _{H^4(\Omega )} \Vert \psi \Vert _{\Omega _h} \end{aligned}$$
    (142)

    We used the estimate

    $$\begin{aligned} |||\phi _h |||^2_{h,\bigstar }&\lesssim |||\phi _h -\phi |||^2_{h,\bigstar } + |||\phi |||^2_{h,\bigstar } \end{aligned}$$
    (143)
    $$\begin{aligned}&\lesssim (1 + h^{-3} )|||\phi _h -\phi |||^2_{h} + |||\phi |||^2_{h} + \Vert T(\phi ) \Vert ^2_{\partial \Omega _h} + {h^{-6} \Vert \phi \Vert ^2_{\partial \Omega _h} } \end{aligned}$$
    (144)
    $$\begin{aligned}&\lesssim (1 + h^{-3})h^4 \Vert \phi \Vert ^2_{H^4(\Omega )} + |||\phi |||^2_{h} + ({1 + h^{-6} \delta ^2} )\Vert \phi \Vert ^2_{H^4(\Omega )} \end{aligned}$$
    (145)
    $$\begin{aligned}&\lesssim {\Vert \psi \Vert ^2_{\Omega _h}} \end{aligned}$$
    (146)

    where we used (89) to conclude that \(\Vert T(\phi ) \Vert ^2_{\partial \Omega _h} + h^{-6} \Vert \phi \Vert ^2_{\partial \Omega _h} \lesssim (1 + h^{-6} \delta ^2)\Vert \phi \Vert ^2_{H^4(\Omega )}\). Collecting the estimates of Terms IIII completes the proof. \(\square \)

Numerics

Implementation

We consider two higher order approximations of the boundary: a piecewise cubic \(C^0\) approximation or a piecewise cubic \(C^1\) approximation. The steps to create the approximate boundary are as follows.

  1. 1.

    The elements cut by the boundary are located, Fig. 2.

  2. 2.

    Straight segments connecting the intersection points between the boundary and the elements are established, and the geometrical object inside the domain is triangulated for ease if integration, Fig. 3.

  3. 3.

    The end points of the boundary segments and the inclinations at the endpoints (computed by use of tangent vectors) is used to obtain a \(C^1\) interpolant of the boundary, Fig. 4. (This step is skipped in the case of a \(C^0\) approximation of the boundary.)

  4. 4.

    The geometry is approximated by a cubic triangle, interpolating the exact boundary (\(C^0\) case) or the spline boundary (\(C^1\) case), Fig 5.

Note that the approximation of the boundary may partly land outside the element. In such cases, the basis functions of the element containing the straight segment is used also outside of the element.

Fig. 2
figure 2

Element intersected by the boundary (dashed)

Fig. 3
figure 3

Straight line approximation of the boundary (dotted) and triangulation for integration purposes

Fig. 4
figure 4

Cubic spline approximation of the boundary (solid line)

Example

We consider a circular simply supported plate under uniform load p. The plate is of radius \(R=1/2\) and has its center at \(x=1/2\), \(y=1/2\). Defining r as the distance from the midpoint we then have the exact solution

$$\begin{aligned} u=\frac{pR^4}{64\kappa }\left( 1-\left( \frac{r}{R}\right) ^2\right) \left( \frac{5+\nu }{1+\nu }-\left( \frac{r}{R}\right) ^2\right) \end{aligned}$$

see, e.g., [18]. The constitutive parameters were chosen as \(E=10^2\), \(\nu =0.3\), \(t=10^{-1}\), and the stabilization parameters as \(\beta = 10^{-2}\), \(\gamma = 10^2(2\kappa + 2\kappa \nu (1-\nu )^{-1})\).

We compare the convergence in normalized (\(|| u- u_h||/||u||\)) \(L_2\), \(H^1\) and piecewise \(H^2\) norms in Fig. 6. These norms are computed on the discrete geometry, for simplicity the straight segment geometry. The solid lines indicate second, third, and fourth order convergence, respectively from top to bottom, and we note that we observe a slightly higher than optimal rate of convercence of about \(O(h^{1/2})\) in all norms. We note that the continuity of the approximation of the boundary seems not to be crucial as the convergence curves are very close. In Fig. 7 we show an elevation of the solution on one of the meshes in the sequence.

Fig. 5
figure 5

Isoparametric cubic triangle approximation of the geometry

Fig. 6
figure 6

Convergence in normalized \(L_2\), \(H^1\), and piecewise \(H^2\) norms

We also show the influence of the parameter \(\beta \) on a fixed mesh (coarse, 386 active nodes). In Fig. 8 we show the condition number as \(\beta \) increases from its critical number, the number for which the system matrix is singular. For lower values of \(\beta \) there are negative eigenvalues in the system matrix. We note that as \(\beta \) increases, the condition number will eventually increase again after an initial drop. In Fig. 9 we show the \(H^1\) error which is more sensitive to the increase in \(\beta \). We remark that this effect, however, does not affect the convergence rate.

Fig. 7
figure 7

Elevation of the discrete solution on one of the meshes in the sequence

Fig. 8
figure 8

Variation of the condition number as a function of \(\beta \)

Fig. 9
figure 9

Variation of the \(H^1\) error as a function of \(\beta \)

Concluding remarks

We have proposed and analyzed a cut finite element method for a rectangular plate element, allowing for curved boundaries. The analysis shows that the method is optimally order convergent and stable. Two different approximations of the boundary have been tested, a standard cubic interpolation of the exact boundary and a cubic spline approximation leading to a continuously differentiable approximation of the boundary. Numerical results and theory indicate that the continuity of the boundary approximation is not crucial. With our method, the simple rectangular \(C^1\) element has greatly increased its practical applicability.

References

  1. Argyris JH, Fried I, Scharpf DW. The tuba family of plate elements for the matrix displacement method. Aeronaut J Roy Aeronaut Soc. 1968;72(692):701–9.

    Google Scholar 

  2. Arnold DN, Falk RS. A uniformly accurate finite element method for the Reissner-Mindlin plate. SIAM J Numer Anal. 1989;26(6):1276–90.

    MathSciNet  MATH  Google Scholar 

  3. Babuška I, Pitkäranta J. The plate paradox for hard and soft simple support. SIAM J Math Anal. 1990;21(3):551–76.

    MathSciNet  MATH  Google Scholar 

  4. Beirão da Veiga L, Niiranen J, Stenberg R. A family of \(C^0\) finite elements for Kirchhoff plates. II. Numerical results. Comput Methods Appl Mech Eng. 2008;197(21–24):1850–64.

    MATH  Google Scholar 

  5. Bernadou M, Boisserie J-M. Curved finite elements of class \(C^1\): implementation and numerical experiments. I. Construction and numerical tests of the interpolation properties. Comput Methods Appl Mech Eng. 1993;106(1–2):229–69.

    MATH  Google Scholar 

  6. Bogner FK, Fox RL, Schmit LA. The generation of interelement compatible stiffness and mass matrices by the use of interpolation formulae. In: Proc. Conf. Matrix Methods in Struct. Mech., AirForce Inst. of Tech., Wright Patterson AF Base, Ohio. 1965. p. 397–444.

  7. Brenner SC, Scott LR. The mathematical theory of finite element methods, volume 15 of Texts in Applied Mathematics. 3rd ed. New York: Springer; 2008.

    Google Scholar 

  8. Brenner SC, Sung L-Y. \(C^0\) interior penalty methods for fourth order elliptic boundary value problems on polygonal domains. J Sci Comput. 2005;22(23):83–118.

    MathSciNet  MATH  Google Scholar 

  9. Burman E, Hansbo P, Larson MG. A cut finite element method with boundary value correction. Math Comp. 2018;87(310):633–57.

    MathSciNet  MATH  Google Scholar 

  10. Clough RW, Tocher JL. Finite element stiffness matrices for analysis of plate bending. In: Matrix methods in structural mechanics, number AFFDL-TR-66-80. 1966. p. 515–45.

  11. Di Pietro DA, Ern A. Mathematical aspects of discontinuous Galerkin methods, volume 69 of Mathématiques & Applications (Berlin) [Mathematics & Applications]. Heidelberg: Springer; 2012.

  12. Durán R, Liberman E. On mixed finite element methods for the Reissner-Mindlin plate model. Math Comp. 1992;58(198):561–73.

    MathSciNet  MATH  Google Scholar 

  13. Ern A, Guermond J-L. Theory and practice of finite elements, volume 159 of Applied Mathematical Sciences. New York: Springer; 2004.

    MATH  Google Scholar 

  14. Gilbarg D, Trudinger NS. Elliptic partial differential equations of second order. Classics in Mathematics. Berlin: Springer; 2001. Reprint of the 1998 edition.

  15. Hansbo P, Larson MG. A discontinuous Galerkin method for the plate equation. Calcolo. 2002;39(1):41–59.

    MathSciNet  MATH  Google Scholar 

  16. Hansbo P, Larson MG, Larsson K. Cut finite element methods for linear elasticity problems. In: Geometrically unfitted finite element methods and applications, volume 121 of Lect. Notes Comput. Sci. Eng.. Cham: Springer; 2017. p. 25–63.

  17. Huang J, Huang X. A hybridizable discontinuous Galerkin method for Kirchhoff plates. J Sci Comput. 2019;78(1):290–320.

    MathSciNet  MATH  Google Scholar 

  18. Huston R, Josephs H. Practical stress analysis in engineering design. Boca Raton: CRC Press; 2008.

    Google Scholar 

  19. Massing A, Larson MG, Logg A, Rognes ME. A stabilized Nitsche fictitious domain method for the Stokes problem. J Sci Comput. 2014;61(3):604–28.

    MathSciNet  MATH  Google Scholar 

  20. Nitsche JA. Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind. Abh Math Univ Hamburg. 1971;36:9–15.

    MathSciNet  MATH  Google Scholar 

  21. Stein EM. Singular integrals and differentiability properties of functions. Princeton Mathematical Series, No. 30. Princeton: Princeton University Press; 1970.

  22. Wells GN, Dung NT. A \(C^0\) discontinuous Galerkin formulation for Kirchhoff plates. Comput Methods Appl Mech Eng. 2007;196(35–36):3370–80.

    MATH  Google Scholar 

  23. Zhang S. On the full \(C_1\)-\(Q_k\) finite element spaces on rectangles and cuboids. Adv Appl Math Mech. 2010;2(6):701–21.

    MathSciNet  MATH  Google Scholar 

  24. Zienkiewicz OC, Taylor RL. The finite element method. Vol. 2: solid mechaincs. 5th ed. Oxford: Butterworth-Heinemann; 2000.

    Google Scholar 

  25. Hansbo A, Hansbo P, Larson MG. A finite element method on composite grids based on Nitsche’s method. ESAIM: Math. Model. Numer. Anal. 2003;37(3):495–514.

    MathSciNet  MATH  Google Scholar 

Download references

Funding

Open access funding provided by Jönköping University. Funding was provided by Stiftelsen för Strategisk Forskning (AM13-0029), Vetenskapsrådet (2013-4708, 2017-03911, 2018-05262), Engineering and Physical Sciences Research Council (EP/P01576X/1) and Swedish Research Programme Essence.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Peter Hansbo.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix: Some inequalities

Appendix: Some inequalities

Let \(\omega \subset \partial \Omega \) and define the cylindrical tubular neighborhood

$$\begin{aligned} U_\delta (\omega ) = \{ x \in U_\delta (\partial \Omega ) : p(x) \in \omega \} \qquad \delta \in (0,\delta _0] \end{aligned}$$
(147)

Let \(\Gamma _1\) and \(\Gamma _2\) be two surface segments in \(U_\delta (\omega )\) with unit normals \(n_i\), such that the closest point mapping \(p:\Gamma _i \rightarrow \omega \) is a bijection with inverse \(p_i^{-1}: \omega \rightarrow \Gamma _i\) and there is constant such that

$$\begin{aligned} 1 \lesssim \min _{x \in \Gamma _i} |n_ i(x) \cdot n(x)| \end{aligned}$$
(148)

where \(n(x) = n\circ p(x)\). We can then define a bijection \(q:\Gamma _1 \rightarrow \Gamma _2\) by \(q=p_2^{-1} p\). For each \(x \in \Gamma _1\) let \(I_x\) be the line segment with endpoints \(x\in \Gamma _1\) and \(q(x) \in \Gamma _2\). We then have

$$\begin{aligned} v(x) = v(q(x)) + \int _{I_x} \nabla _n v ds \end{aligned}$$
(149)

where we integrate along I from q(x) to x. Using the Cauchy Schwarz inequality we get

$$\begin{aligned} v^2(x)&\lesssim v^2(q(x)) + \left( \int _{I_x} \nabla _n v \right) ^2 \lesssim v^2(q(x)) + \int _{I_x} 1 \int _{I_x} |\nabla _n v|^2 \nonumber \\&\lesssim v^2(q(x)) + \delta \int _{I_x} |\nabla _n v|^2 \end{aligned}$$
(150)

Integrating over \(\Gamma _1\),

$$\begin{aligned} \int _{\Gamma _1} v^2(x)&\lesssim \int _{\Gamma _1} v^2(q(x)) + \delta \int _{\Gamma _1} \int _{I_x} |\nabla _n v|^2 \end{aligned}$$
(151)
$$\begin{aligned}&\lesssim \int _{\Gamma _2} v^2(y_2) + \int _{S} |\nabla _n v |^2 \end{aligned}$$
(152)

Here \(S = \cup _{x\in \Gamma _1} I_x\) is the domain between the surfaces \(\Gamma _1\) and \(\Gamma _2\), and we changed coordinates to integration over \(\Gamma _2\) and S equipped with Euclidian measure. We conclude that

$$\begin{aligned} \boxed {\Vert v \Vert ^2_{\Gamma _1} \lesssim \Vert v \Vert ^2_{\Gamma _2} + \delta \Vert \nabla v \Vert ^2_{S}} \end{aligned}$$
(153)

In applications, it is of convenient to simply replace S by the larger domain \(U_\delta (\omega )\).

Typical applications include taking \(\Gamma _2 = \omega = \partial \Omega \) and \(\Gamma _1 = \partial \Omega _h\), which gives

$$\begin{aligned} \Vert v \Vert ^2_{\partial \Omega _h} \lesssim \Vert v \Vert ^2_{\partial \Omega } + \delta \Vert \nabla v \Vert ^2_{S} \end{aligned}$$
(154)

For \(v \in V_h\) we have \(S\subset O_h = \cup _{T \in {\mathcal {T}}_h} T\) since we assume the \(\Omega \cup \Omega _h \subset O_h\) and in particular \(S\subset U_\delta (\partial \Omega ) \cap O_h\) and thus we get

$$\begin{aligned} \Vert v \Vert ^2_{\partial \Omega _h} \lesssim \Vert v \Vert ^2_{\partial \Omega } + \delta \Vert \nabla v \Vert ^2_{U_\delta (\partial \Omega ) \cap O_h} \end{aligned}$$
(155)

For \(v\in H^1(\Omega )\), with extension to \({\mathbb {R}}^2\) also denoted by v, we have \(S \subset U_\delta (\partial \Omega )\) and we get

$$\begin{aligned} \Vert v \Vert ^2_{\partial \Omega _h} \lesssim \Vert v \Vert ^2_{\partial \Omega } + \delta \Vert \nabla v \Vert ^2_{U_\delta (\partial \Omega ) \cap O_h} \end{aligned}$$
(156)

In this paper we have recall that \(\partial \Omega _h \subset U_\delta (\partial \Omega )\), with \(\delta \sim h^4\), see (10).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Burman, E., Hansbo, P. & Larson, M.G. Cut Bogner-Fox-Schmit elements for plates. Adv. Model. and Simul. in Eng. Sci. 7, 27 (2020). https://doi.org/10.1186/s40323-020-00164-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40323-020-00164-3

Keywords