Next Article in Journal
A Generalized Viscosity Inertial Projection and Contraction Method for Pseudomonotone Variational Inequality and Fixed Point Problems
Previous Article in Journal
Dissipativity of Fractional Navier–Stokes Equations with Variable Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New High Accuracy Analysis of a Double Set Parameter Nonconforming Element for the Clamped Kirchhoff Plate Unilaterally Constrained by an Elastic Obstacle

School of Mathematics and Statistics, Zhengzhou University, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(11), 2038; https://doi.org/10.3390/math8112038
Submission received: 26 October 2020 / Revised: 6 November 2020 / Accepted: 13 November 2020 / Published: 16 November 2020
(This article belongs to the Section Computational and Applied Mathematics)

Abstract

:
In this paper, a non- C 0 double set parameter finite element method is presented for the clamped Kirchhoff plate with an elastic unilateral obstacle. A new high accuracy error estimate with order O ( h 2 ) in the broken energy norm is derived by use of a series of novel approaches, including some special features of the element and an incomplete biquadratic interpolation operator. At the same time, some experimental results are provided to verify the theoretical analysis.

1. Introduction

Let Ω R 2 be a bounded convex polygonal domain with a Lipschitz continuous boundary Γ . We consider a thin flat plate Ω × ( d / 2 , d / 2 ) clamped on its lateral boundary Γ × ( d / 2 , d / 2 ) , where d > 0 is the thickness of the plate and assumed to be small. Assume that the plate is subject to a load f L 2 ( Ω ) ,   f 0   a . e .   in   Ω and constrained unilaterally by an elastic obstacle ψ L 2 ( Ω ) , then displacement function u can be obtained by solving the weak formulation (cf. [1]):
Find     u V ,   such   that     a ( u , v ) + ( u , v ) = ( f , v ) ,     v V ,
where V = H 0 2 ( Ω ) , a ( w , v ) = Ω [ Δ w Δ v + ( 1 ν ) ( 2 w x y v x y w x x v y y w y y v x x ) ] d x d y , ( u , v ) = Ω κ ( u ψ ) v d x d y ,   t = min { t , 0 } , ( f , v ) = Ω f v d x d y , ν ( 0 , 1 / 2 ) is Poisson’s ratio and κ κ 0 > 0 describes the stiffness of the obstacle.
It was shown in [1] that the displacement function u can also be derived by the minimization problem:
Find     u V ,   such   that   v V     J ( u ) J ( v )
with the total energy J ( v ) = 1 2 a ( v , v ) + 1 2 Ω κ [ ( v ψ ) ] 2 d x d y ( f , v ) . Furthermore, the problem (2) is equivalent to the following fourth-order variational inequality of the second kind:
Find     u V ,   such   that a ( u , v u ) + Ω κ 2 [ ( v ψ ) ] 2 d x d y Ω κ 2 [ ( u ψ ) ] 2 d x d y ( f , v u ) , v V .
In particular, in the limit κ , the obstacle becomes rigid, and the problem reduces to the constrained minimization:
Find     u K * ,   such   that u = arg min v K * [ 1 2 a ( v , v ) ( f , v ) ]
with K * = { v V : v ψ   in   Ω } , which is the classical displacement obstacle problem of the clamped Kirchhoff plate (cf. [2]) and is equivalent to the following fourth-order variational inequality of the first kind:
Find     u K * ,   such   that a ( u , v u ) ( f , v u ) ,     v K * .
It is well known that nonconforming finite element methods (FEMs) are an attractive option for solving high order differential equations since the smoothness requirement on finite element functions is weakened. As we know, nonconforming FEMs for the problem (4), the plate obstacle problem with a rigid obstacle, have achieved fruitful results [3,4,5,6,7,8,9], but very few for the elastic obstacle problem (2), except [10]. Compared with the conforming FEMs studied in [1,11], C 0 continuous and discontinuous nonconforming (i.e., C 0 and non- C 0 nonconforming) FEMs were discussed in [10], and convergence analyses were established. On the other hand, the high accuracy analysis of nonconforming FEMs has been an active research area in practical computations, and many high accuracy results have been derived for variational inequalities and the boundary value problem [12,13,14,15,16,17,18,19,20]. However, for the plate obstacle problem (4) with a rigid obstacle, the exact solution u only belongs to H 3 ( Ω ) C 2 ( Ω ) instead of H l o c 4 ( Ω ) [21,22]; thus, the lack of H 4 regularity makes it impossible to develop high accuracy analysis. Fortunately, it was shown in [10] that the solution is more regular if the obstacle is elastic. In fact, the solution of the problem (2) belongs to H 4 ( Ω ) when the largest interior angles of the domain Ω are smaller than 126.283696 ° , and thus, the high accuracy analysis of FEMs for the problem (2) is possible and worth exploring.
In this paper, we attempt to present a high accuracy analysis of nonconforming double set parameter FEMs for the obstacle problem (2). The double set parameter method presented in [23,24] is a useful method to construct unconventional elements, and it involves two sets of parameters. The first set is chosen to meet the convergence requirements according to the generalized patch test [25] or F-E-M-test [26], while the second set is selected to be simple so that the total number of unknowns in the resulting discrete system is small. Up to now, several nonconforming double set parameter plate elements have been successfully applied to deal with the two-sided displacement obstacle problem of the clamped plate [8], plate bending problems [16,27], the linear elasticity problem [28], the fourth-order elliptic singular perturbation problem [29,30,31] and so on. In this paper, a non- C 0 nonconforming double set parameter plate element is employed for the elastic obstacle problem (2). We develop a series of novel approaches including some special features of the element and an incomplete biquadratic interpolation operator so as to get the high accuracy result with order O ( h 2 ) in the broken energy norm. Furthermore, we carry out a numerical experiment to show the performance of the proposed method.

2. A Twelve Parameter Double Set Parameter Element and Its Typical Properties

In the beginning, we introduce a twelve parameter double set parameter element briefly. The readers are referred to [27] for details.
Assume that T h is a regular rectangular decomposition of Ω . For a given K T h , let its center be ( x K , y K ) , four vertices be a i ( x i , y i ) , four edges be F i = a i a i + 1 ¯   ( i = 1 , 2 , 3 , 4 ( m o d 4 ) ) , and the edge length be 2 h x and 2 h y , respectively. Moreover, we denote h = max K T h { h x , h y } . Besides, assume that K ^ = [ 1 , 1 ] × [ 1 , 1 ] is the reference element on the ξ η plane with a center point ( 0 , 0 ) , four vertices a ^ 1 = ( 1 , 1 ) , a ^ 2 = ( 1 , 1 ) , a ^ 3 = ( 1 , 1 ) , and a ^ 4 = ( 1 , 1 ) , and four edges F ^ 1 = a ^ 1 a ^ 2 ¯ , F ^ 2 = a ^ 2 a ^ 3 ¯ , F ^ 3 = a ^ 3 a ^ 4 ¯ , and F ^ 4 = a ^ 4 a ^ 1 ¯ . Then, there exists an affine mapping J : K ^ K :
x = x K + h x ξ y = y K + h y η
satisfying J ( K ^ ) = K , v ( x , y ) = v ^ ( ξ , η ) .
The shape function space on the element K is taken as:
P ( K ) = P 3 ( K ) { x 4 , y 4 } = S p a n { p 1 , , p 12 } ,
where p 1 = 1 4 ( 1 ξ ) ( 1 η ) ,     p 2 = 1 4 ( 1 + ξ ) ( 1 η ) ,     p 3 = 1 4 ( 1 + ξ ) ( 1 + η ) ,     p 4 = 1 4 ( 1 ξ ) ( 1 + η ) ,     p 5 = ( 1 ξ 2 ) ,     p 6 = ( 1 η 2 ) ,     p 7 = ( 1 ξ 2 ) η ,     p 8 = ( 1 η 2 ) ξ ,     p 9 = ( 1 ξ 2 ) ξ ,     p 10 = ( 1 η 2 ) η , p 11 = ( 1 ξ 2 ) ξ 2 ,     p 12 = ( 1 η 2 ) η 2 , and P m ( K ) denotes the space consisting of piecewise polynomials of degree m on element K.
The degrees of freedom are selected as:
D ( v ) = ( d 1 ( v ) , d 2 ( v ) , , d 12 ( v ) ) ,
where:
d i ( v ) = v ( a i ) , i = 1 , 2 , 3 , 4 ,
d 5 ( v ) = 1 h x F 1 v d s = 1 1 v ^ ( ξ , 1 ) d ξ ,     d 6 ( v ) = 1 h y F 2 v d s = 1 1 v ^ ( 1 , η ) d η ,
d 7 ( v ) = 1 h x F 3 v d s = 1 1 v ^ ( ξ , 1 ) d ξ ,         d 8 ( v ) = 1 h y F 4 v d s = 1 1 v ^ ( 1 , η ) d η ,
d 9 ( v ) = h y h x F 1 v n d s = 1 1 v ^ η ( ξ , 1 ) d ξ ,     d 10 ( v ) = h x h y F 2 v n d s = 1 1 v ^ ξ ( 1 , η ) d η ,
d 11 ( v ) = h y h x F 3 v n d s = 1 1 v ^ η ( ξ , 1 ) d ξ ,     d 12 ( v ) = h x h y F 4 v n d s = 1 1 v ^ ξ ( 1 , η ) d η .
For v P ( K ) , let:
v = i = 1 12 β i p i = ( P , b ) ,
where b = ( β 1 , β 2 , · · · , β 12 ) R 12 , P = ( p 1 , p 2 , · · · , p 12 ) .
Substituting (9) into (8), we get:
D ( v ) = M b
with:
M = 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 1 0 0 4 3 0 4 3 0 0 0 4 15 0 0 1 1 0 0 4 3 0 4 3 0 0 0 4 15 0 0 1 1 4 3 0 4 3 0 0 0 4 15 0 1 0 0 1 0 4 3 0 4 3 0 0 0 4 15 1 2 1 2 1 2 1 2 0 4 4 3 0 0 4 0 4 1 2 1 2 1 2 1 2 4 0 0 4 3 4 0 4 0 1 2 1 2 1 2 1 2 0 4 4 3 0 0 4 0 4 1 2 1 2 1 2 1 2 4 0 0 4 3 4 0 4 0 .
It is easy to check that | M | = 2 24 3 4 * 5 2 0 ; thus b = M 1 D ( v ) , i.e., it holds that:
v P ( K ) ,   v = ( P , M 1 D ( v ) ) .
Next, we take nodal parameters as:
Q ( v ) = ( v 1 , v 1 x , v 1 y , v 2 , v 2 x , v 2 y , v 3 , v 3 x , v 3 y , v 4 , v 4 x , v 4 y ) ,
where v i , v i x , v i y are the function values of v and its first derivatives at vertices a i , respectively.
Then, approximating D ( v ) by a linear combination of the nodal parameters Q ( v ) as follows:
d i ( v ) = v i ,   ( i = 1 , 2 , 3 , 4 ) ,
for d i ( v )   ( i = 5 , 6 , 7 , 8 ) , the numerical integrating values of cubic Hermite interpolation polynomials on the corresponding sides result in:
d 5 ( v ) = v 1 + v 2 + h x 3 ( v 1 x v 2 x ) + O ( h 4 | v | 4 , K , ) , d 6 ( v ) = v 2 + v 3 + h y 3 ( v 2 y v 3 y ) + O ( h 4 | v | 4 , K , ) , d 7 ( v ) = v 3 + v 4 + h x 3 ( v 4 x v 3 x ) + O ( h 4 | v | 4 , K , ) , d 8 ( v ) = v 1 + v 4 + h y 3 ( v 1 y v 4 y ) + O ( h 4 | v | 4 , K , ) ,
for d i ( v )   ( i = 9 , 10 , 11 , 12 ) , the trapezoidal rule of numerical integration giving:
d 9 ( v ) = h y ( v 1 y + v 2 y ) + O ( h 3 | v | 3 , K , ) , d 10 ( v ) = h x ( v 2 x + v 3 x ) + O ( h 3 | v | 3 , K , ) , d 11 ( v ) = h y ( v 3 y + v 4 y ) + O ( h 3 | v | 3 , K , ) , d 12 ( v ) = h x ( v 1 x + v 4 x ) + O ( h 3 | v | 3 , K , ) .
The above discretizations can be rewritten in matrix form as:
D ( v ) = G Q ( v ) + E ( v ) ,
where:
G = 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 h x 3 0 1 h x 3 0 0 0 0 0 0 0 0 0 0 1 0 h y 3 1 0 h y 3 0 0 0 0 0 0 0 0 0 1 h x 3 0 1 h x 3 0 1 0 h y 3 0 0 0 0 0 0 1 0 h y 3 0 0 h y 0 0 h y 0 0 0 0 0 0 0 0 0 0 h x 0 0 h x 0 0 0 0 0 0 0 0 0 0 0 0 h y 0 0 h y 0 h x 0 0 0 0 0 0 0 0 h x 0 .
and:
E ( v ) = ( 0 , 0 , 0 , 0 , ε 1 ( v ) , ε 1 ( v ) , ε 1 ( v ) , ε 1 ( v ) , ε 2 ( v ) , ε 2 ( v ) , ε 2 ( v ) , ε 2 ( v ) )
with ε 1 ( v ) = O ( h 4 | v | 4 , K , ) , ε 2 ( v ) = O ( h 3 | v | 3 , K , ) . From (18), it is obvious that the term E ( v ) does not affect the convergence properties, and we neglect it and introduce a new set of parameters:
Q ¯ ( v ) = ( q 1 ( v ) , q 2 ( v ) , · · · , q 12 ( v ) )
such that:
D ( v ) = G Q ¯ ( v ) ,
where q 1 ( v ) = v ( a 1 ) ,   q 4 ( v ) = v ( a 2 ) ,   q 7 ( v ) = v ( a 3 ) ,   q 10 ( v ) = v ( a 4 ) ,   q 2 ( v ) = v ¯ x ( a 1 ) ,   q 5 ( v ) = v ¯ x ( a 2 ) , q 8 ( v ) = v ¯ x ( a 3 ) ,   q 11 ( v ) = v ¯ x ( a 4 ) ,   q 3 ( v ) = v ¯ y ( a 1 ) ,   q 6 ( v ) = v ¯ y ( a 2 ) ,   q 9 ( v ) = v ¯ y ( a 3 ) ,   q 12 ( v ) = v ¯ y ( a 4 ) , and v ¯ x ( a ) and v ¯ y ( a ) have perturbations on v x ( a ) and v y ( a ) , respectively. The readers are referred to [23,24] for details.
Then, using (10) and (19), we obtain:
b = M 1 G Q ¯ ( v ) ,
which in conjunction with (9) leads to the real shape function v. Thus, the unknowns of the discrete system are Q ¯ ( v ) , and the real shape function space P ¯ ( K ) is a subspace of P ( K ) :
P ¯ ( K ) = { v P ( K ) ; v = ( P , M 1 G q ) ,   q R 12 } .
It is easy to check that r a n k ( G ) = 10 , i.e., the dimension of P ¯ ( K ) is 10, and v P ¯ ( K ) can be expressed as v = i = 1 12 β i p i , where:
β 1 = q 1 , β 2 = q 4 , β 3 = q 7 , β 4 = q 10 , β 5 = h x 8 ( q 2 q 5 q 8 + q 11 ) , β 6 = h y 8 ( q 3 + q 6 q 9 q 12 ) , β 7 = h x 8 ( q 2 + q 5 q 8 + q 11 ) , β 8 = h y 8 ( q 3 + q 6 q 9 + q 12 ) , β 9 = 1 8 ( q 1 + q 4 + q 7 q 10 ) + h x 8 ( q 2 q 5 q 8 q 11 ) + h y 24 ( q 3 + q 6 q 9 + q 12 ) , β 10 = 1 8 ( q 1 q 4 + q 7 + q 10 ) + h x 24 ( q 2 + q 5 q 8 + q 11 ) + h y 8 ( q 3 q 6 q 9 q 12 ) , β 11 = 0 , β 12 = 0 .
Obviously, v P ¯ ( K ) is a cubic polynomial. In other words, P ¯ ( K ) is a completely cubic polynomial space. We denote the corresponding FE space by V h :
V h = { v ; v | K = ( P , M 1 G q ) , q R 12 , K T h , v ( a ) = v ¯ x ( a ) = v ¯ y ( a ) = 0 , for   any   nodes   a Γ } .
Lemma 1.
For any v h V h , the following properties hold:
( R 1 ) v h is continuous at the vertices of elements and is zero at the vertices on Γ.
( R 2 ) F v h d s is continuous across the element edge F and is zero on F Γ .
( R 3 )   p ( s ) P 1 ( F ) , i.e., p ( s ) is a linear polynomial function defined on the element edge F, and F p ( s ) v h n d s is continuous across the element edge F and is zero on F Γ .
( R 4 ) · h is a norm over V h , where · h = ( K T h | · | 2 , K 2 ) 1 2 .
Proof. 
From (8) and (19), we know that D ( v h ) is continuous between elements and:
d 1 = q 1 , d 2 = q 4 , d 3 = q 7 , d 4 = q 10 , d 5 = q 1 + q 4 + h x 3 ( q 2 q 5 ) , d 6 = q 4 + q 7 + h y 3 ( q 6 q 9 ) , d 7 = q 7 + q 10 + h x 3 ( q 11 q 8 ) , d 8 = q 1 + q 10 + h y 3 ( q 3 q 12 ) , d 9 = h y ( q 3 + q 6 ) , d 10 = h x ( q 5 + q 8 ) , d 11 = h y ( q 9 + q 12 ) , d 12 = h x ( q 2 + q 11 ) .
Furthermore, it follows from (22) that:
F 1 s v h n d s = h x 3 ( q 3 q 6 ) , F 2 s v h n d s = h y 3 ( q 8 q 5 ) , F 3 s v h n d s = h x 3 ( q 9 q 12 ) , F 4 s v h n d s = h y 3 ( q 2 q 11 ) ,
which together with (24) and the definition of d i   ( i = 1 , 2 , , 12 ) yields ( R 1 ) , ( R 2 ) , and ( R 3 ) .
At the same time, suppose v h h = 0 , then for any K T h , it holds that v h K P 1 ( K ) . From ( R 1 ) and ( R 2 ) , we can immediately get v h 0 ; thus, ( R 4 ) holds true.  □

3. High Accuracy Analysis

We consider the discrete approximation form of the variational inequality (3) as:
Find     u h V h ,   such   that   v h V h a ( u h , v h u h ) + K T h K κ 2 [ ( v h ψ ) ] 2 d x d y K T h K κ 2 [ ( u h ψ ) ] 2 d x d y ( f , v h u h ) ,
where a h ( w h , v h ) = K T h K [ Δ w h Δ v h + ( 1 ν ) ( 2 w h x y v h x y w h x x v h y y w h y y v h x x ) ] d x d y .
It was shown in [1] that the problem (26) is equivalent to the discrete approximation of the plate problem:
Find   u h V h ,   such   that   v h V h , a h ( u h , v h ) + K T h K κ ( u h ψ ) v h d x d y = ( f , v h ) .
which is also equivalent to the minimization problem:
Find     u h V h ,   such   that   v h V h     J h ( u h ) J h ( v h )
where the energy functional J h ( v h ) = 1 2 a h ( v h , v h ) + 1 2 K T h K κ [ ( v h ψ ) ] 2 d x d y ( f , v h ) . Since a h ( · , · ) is coercive on V h , we can get the lower bound for J h ( v h ) and J h ( v h ) (as v h h ). Then, using the elementary inequality ( t s ) ( t s ) 0 and an argument similar to that used in [10], we have:
Theorem 1.
The problem (26) has a unique solution u h . Moreover, u h h and ( u h ψ ) 0 are uniformly bounded independently of h.
In what follows, we will give error estimates for (26).
Theorem 2.
Assume that u and u h are the solutions of (3) and (26), respectively, and the largest internal angle of Ω is ω < 126.283696 ° . Then, we have:
u u h h C h 2 ( u 4 + f 0 + ( u ψ ) 0 ) .
Here and later, positive constant C is independent of h and may be different at each appearance.
Proof. 
Employing the triangle inequality yields:
u u h h u Π h u h + Π h u u h h ,
where Π h is the associated interpolation operator on V h . Let w h = Π h u u h . We have by (27) that:
w h h 2 C a h ( w h , w h ) = C a h ( Π h u u , w h ) + a h ( u , w h ) a h ( u h , w h ) C Π h u u h w h h + C a h ( u , w h ) + K T h K κ ( u h ψ ) w h d x d y ( f , w h ) .
Now, we focus on the estimate of the second term on the right-hand side of (31). Obviously,
a h ( u , w h ) ( f , w h ) + K T h K κ ( u h ψ ) w h d x d y = K T h K M ( u ) w h n d s + ϵ K T h K 2 u n s w h s d s       K T h K ( Δ u ) · w h d x d y ( f , w h ) + K T h K κ ( u h ψ ) w h d x d y ,
where M ( u ) = Δ u ( 1 ν ) 2 u s 2 and ϵ = ( 1 ν ) . In order to estimate (32), we notice that the solution u of the problem (3) satisfies (see [10]):
Ω Δ u · v d x d y + Ω κ ( u ψ ) v d x d y = ( f , v ) ,     v H 0 1 ( Ω ) .
However, from the construction of the element and ( R 1 ) ( R 3 ) in Lemma 1, we know that w h H 1 ( Ω ) . In this situation, we need to introduce a piecewise interpolation polynomial of w h as follows:
w h I K P 2 ( K ) { x 2 y , x y 2 } ,     w h I K ( a i ) = w h ( a i ) , F w h I d s = F w h d s ,     F K .
Then, from ( R 1 ) and ( R 2 ) in Lemma 1, we have w h I H 0 1 ( Ω ) , which in conjunction with (32) and (33) yields:
a h ( u , w h ) ( f , w h ) + K T h K κ ( u h ψ ) w h d x d y = K T h K M ( u ) w h n d s + ϵ K T h K 2 u n s ( w h w h I ) s d s K T h K ( Δ u ) · ( w h w h I ) d x d y   K T h K ( Δ u ) · w h I d x d y ( f , w h ) + K T h K κ ( u h ψ ) w h d x d y = K T h K M ( u ) w h n d s + ϵ K T h K 2 u n s ( w h w h I ) s d s K T h K ( Δ u ) · ( w h w h I ) d x d y   + ( f , w h I w h ) Ω κ ( u ψ ) w h I d x d y + K T h K κ ( u h ψ ) w h d x d y = K T h K M ( u ) w h n d s + ϵ K T h K 2 u n s ( w h w h I ) s d s K T h K ( Δ u ) · ( w h w h I ) d x d y   + ( f , w h I w h ) K T h K κ ( u ψ ) ( w h I w h ) d x d y + K T h K κ [ ( u h ψ ) ( u ψ ) ] w h d x d y j = 1 6 ( E r ) j .
Now, we start to estimate ( E r ) j one by one for j = 1 , 2 , , 6 .
Firstly, for ϕ C 0 ( Ω ) , let P 1 F ϕ be the piecewise linear Lagrange interpolation on edge F and P 0 F w h n = 1 F F w h n d s with F = F 1 d s , then the corresponding remainders are:
R 1 F ϕ = ϕ P 1 F ϕ ,     R 0 F w h n = w h n P 0 F w h n .
Applying ( R 3 ) in Lemma 1 leads to:
K T h F K F ( P 1 F M ( u ) ) w h n d s = 0 ,
K T h F K F ( P 1 F M ( u ) ) ( P 0 F w h n ) d s = 0 ,     K T h F K F M ( u ) ( P 0 F w h n ) d s = 0 ,
which follows from:
| ( E r ) 1 | = | K T h F K F ( R 1 F M ( u ) ) ( R 0 F w h n ) d s | | K T h F K ( F ( R 1 F M ( u ) ) 2 d s ) 1 2 ( F ( R 0 F w h n ) 2 d s ) 1 2 | C h 2 u 4 w h h .
Here, the Cauchy–Schwarz inequality and interpolation theorem [32,33] are used.
Secondly, for K T h , we have:
F K F ( P 1 F 2 u n s ) ( w h w h I ) s d s = F 1 ( P 1 F 1 2 u n s ) ( x , y K h y ) ( w h w h I ) x d x F 3 ( P 1 F 3 2 u n s ) ( x , y K + h y ) ( w h w h I ) x d x + F 2 ( P 1 F 2 2 u n s ) ( x K + h x , y ) ( w h w h I ) y d y F 4 ( P 1 F 4 2 u n s ) ( x K h x , y ) ( w h w h I ) y d y .
Employing integration by parts and (34), the first term of (37) can be estimated as:
F 1 ( P 1 F 1 2 u n s ) ( x , y K h y ) ( w h w h I ) x d x = ( P 1 F 1 2 u n s ) ( w h w h I ) | a 2 a 1 F 1 d ( P 1 F 1 2 u n s ) ( x , y K h y ) d x ( w h w h I ) d x = 0 .
Similarly, it holds that:
F K F ( P 1 F 2 u n s ) ( w h w h I ) s d s = 0 ,
which together with interpolation theorem yields:
| ( E r ) 2 | = | ϵ K T h F K F ( 2 u n s P 1 F 2 u n s ) ( w h w h I ) s d s | C h 2 u 4 w h h .
Thirdly, it follows from the definition of w h I in (34) that:
K ( w h w h I ) d x d y = F K F ( w h w h I ) d s = 0 .
Then, let P 0 K ( Δ u ) = 1 K K ( Δ u ) d x d y and K = K 1 d x d y . With a similar argument as ( E r ) 1 and ( E r ) 2 , we obtain:
| ( E r ) 3 | = | K T h K ( ( Δ u ) P 0 K ( Δ u ) ) · ( w h w h I ) d x d y | C h 2 u 4 w h h .
As for ( E r ) 4 and ( E r ) 5 , from the Cauchy–Schwarz inequality and interpolation theorem, it holds that:
| ( E r ) 4 | = | f ( w h I w h ) | f 0 w h w h I 0 C h 2 f 0 w h h
and:
| ( E r ) 5 | κ ( u ψ ) 0 w h w h I 0 C h 2 ( u ψ ) 0 w h h .
Finally, employing the elementary inequality ( t s ) ( t s ) 0 and Theorem 1, we have:
( E r ) 6 = K T h K κ [ ( u h ψ ) ( u ψ ) ] ( Π h u u + u u h ) d x d y K T h K κ [ ( u h ψ ) ( u ψ ) ] ( Π h u u ) d x d y C ( ( u h ψ ) 0 + ( u ψ ) 0 ) Π h u u 0 C ( u ψ ) 0 Π h u u 0 .
Combining (31), (35), and the above bounds of ( E r ) 1 ( E r ) 6 results in:
Π h u u h h 2 C h 4 ( u 4 2 + f 0 2 + ( u ψ ) 0 2 ) + C ( Π h u u h 2 + ( u ψ ) 0 Π h u u 0 ) .
Therefore, the desired result (29) follows from (30), (46) and the interpolation theorem immediately.  □
Remark 1.
w h I plays a key role in the estimations of ( E r ) 2 and ( E r ) 3 in Theorem 2. If using the error estimate method in [10], we can only obtain the convergence result with the order of O ( h ) . In addition, it should be pointed out that the analysis presented herein is also valid for the elastic obstacle problem with j ( v ) = Ω κ [ ( v ψ ) + ] 2 d x d y , t + = max { t , 0 } , and f L 2 ( Ω ) ,   f 0   a . e .   i n   Ω .
Remark 2.
Compared with the work in [10], we obtain the same high accuracy result O ( h 2 ) with a non- C 0 nonconforming element instead of the C 0 nonconforming ACM’s rectangular element. Therefore, our work is an extension of [10], and the requirement for the nonconforming plate element’s smoothness from C 0 continuity to mean C 1 continuity is further reduced ( ( R 1 ) ( R 3 ) implies that the element discussed in this paper is of the mean C 1 type). It can be checked that the high accuracy result in Theorem 2 is no longer true for many classical non- C 0 nonconforming plate elements, such as the Morley triangle element, the Veubeke triangle element, and so on.

4. Numerical Experiment

In this section, we consider the elastic obstacle problem (2) with Ω = [ 0 , 1 ] 2 , f = 10 ,   ν = 0.25 , and an elastic obstacle defined by the function ψ = 0         if     ( x , y ) [ 0.3 ; 0.7 ] 2 , 1         otherwise . Since the exact solution is unknown, we refer the numerical solution on rectangles with h = 1 / 128 as the exact solution to show the errors. The domain Ω is divided into N × N rectangles, and the cases κ = 10 j ,   j = 0 , 1 , 3 , 4 with N = 8 , 16 , 32 , 64 are computed, respectively. The numerical results are given in the following Table 1. Moreover, Figure 1, Figure 2 and Figure 3 illustrate the numerical solution u h with N = 64 , and Figure 4 presents the errors in the logarithm scales. It is obvious that the errors in the energy norm are convergent at order O ( h 2 ) , which coincides with the theoretical analysis in Theorem 2.

Author Contributions

Conceptualization, L.P.; software, L.P.; methodology, L.P. and D.S.; writing, original draft preparation, L.P.; writing, review and editing, D.S. All authors read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 11701523, No. 11671369).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tosone, C.; Maceri, A. The clamped plate with elastic unilateral obstacles: A finite element approach. Math. Model. Methods Appl. Sci. 2003, 13, 1231–1243. [Google Scholar] [CrossRef]
  2. Glowinski, R. Numerical Methods for Nonlinear Variational Inequality Problems; Springer: New York, NY, USA, 1984. [Google Scholar]
  3. Wang, L.H. Some strongly discontinuous nonconforming finite element approximations of a fourth order variational inequality with displacement obstacle. Acta. Numer. Math. 1992, 14, 98–101. [Google Scholar]
  4. Shi, D.Y.; Chen, S.C. Quasi conforming element approximation for a fourth order variational inequality with displacement obstacle. Acta Math. Sci. 2003, 23B, 61–66. [Google Scholar] [CrossRef]
  5. Shi, D.Y.; Chen, S.C. General estimates on nonconforming elements for a fourth order variational problem. Numer. Math. Sin. 2003, 25, 99–106. [Google Scholar]
  6. Brenner, S.C.; Sung, L.Y.; Zhang, Y. Finite element methods for the displacement obstacle problem of clamped plates. Math. Comp. 2012, 81, 1247–1262. [Google Scholar] [CrossRef]
  7. Brenner, S.C.; Sung, L.Y.; Zhang, H.C.; Zhang, Y. A Morley finite element method for the displacement obstacle problem of clamped Kirchhoff plates. J. Comput. Appl. Math. 2013, 254, 31–42. [Google Scholar] [CrossRef]
  8. Shi, D.Y.; Pei, L.F. Double set parameter finite element method for two-sided displacement obstacle problem of clamped plate. J. Math. Anal. Appl. 2016, 436, 203–216. [Google Scholar] [CrossRef]
  9. Shi, D.Y.; Pei, L.F. A new error analysis of Bergan’s energy-orthogonal element for two-sided displacement obstacle problem of clamped plate. J. Math. Anal. Appl. 2016, 442, 339–352. [Google Scholar] [CrossRef]
  10. Han, W.M.; Hua, D.Y.; Wang, L.H. Nonconforming finite element methods for a clamped plate with elastic unilateral obstacle. J. Integral Eq. Appl. 2006, 18, 267–284. [Google Scholar] [CrossRef]
  11. Gustafsson, T.; Stenberg, R.; Videman, J. A stabilised finite element method for the plate obstacle problem. Br. Numer. Math. 2018, 2018, 1–28. [Google Scholar] [CrossRef] [Green Version]
  12. Li, M.X.; Lin, Q.; Zhang, S.H. Superconvergence of finite element method for the Signorini problem. J. Comput. Appl. Math. 2008, 222, 284–292. [Google Scholar] [CrossRef] [Green Version]
  13. Shi, D.Y.; Ren, J.C.; Gong, W. Convergence and superconvergence analysis of a nonconforming finite element method for solving the Signorini problem. Nonlinear Anal. Theor. 2012, 75, 3493–3502. [Google Scholar] [CrossRef]
  14. Shi, D.Y.; Guan, H.B.; Guan, X.F. Superconvergence analysis of finite element method for a second-type variational inequality. J. Appl. Math. 2012, 2012, 1–12. [Google Scholar] [CrossRef]
  15. Xu, C.; Shi, D.Y. Superconvergence analysis of low order nonconforming finite element methods for variational inequality problem with displacement obstacle. Appl. Math. Comp. 2019, 348, 1–11. [Google Scholar] [CrossRef]
  16. Chen, S.C.; Yin, L.; Mao, S.P. An anisotropic, superconvergent nonconforming plate finite element. J. Comput. Appl. Math. 2008, 220, 96–110. [Google Scholar] [CrossRef] [Green Version]
  17. Shi, D.Y.; Wu, Y.M. Uniform superconvergent analysis of a new mixed finite element method for nonlinear Bi-wave singular perturbation problem. Appl. Math. Lett. 2019, 93, 131–138. [Google Scholar] [CrossRef]
  18. Hu, J.; Shi, Z.C.; Yang, X.Q. Superconvergence of both two and three dimensional rectangular Morley elements for biharmonic equations. Mathematics 2015, 132, 491–509. [Google Scholar]
  19. Ruggieri, M.; Speciale, M.P. Approximate analysis of a nonlinear dissipative model. Acta Appl. Math. 2014, 132, 549–559. [Google Scholar] [CrossRef]
  20. Shi, D.Y.; Mao, S.P.; Chen, S.C. On the anisotropic accuracy analysis of ACM’s nonconforming finite element. J. Comput. Math. 2005, 23, 185–198. [Google Scholar]
  21. Frehse, J. On the regularity of the solution of the biharmonic variational inequality. Manuscripta Math. 1973, 9, 91–103. [Google Scholar] [CrossRef]
  22. Schild, B. A regularity result for polyharmonic variational inequalities with thin obstacles. Ann. Scuola Norm. Sup. Pisa Cl. Sci. 1984, 11, 87–122. [Google Scholar]
  23. Chen, S.C.; Shi, Z.C. Double set parameter method of constructing stiffness matrices. Numer. Math. Sin. 1991, 3, 286–296. [Google Scholar]
  24. Chen, S.C.; Shi, D.Y. The error estimates for double set parameter elemet. J. Eng. Math. 2005, 22, 599–605. [Google Scholar]
  25. Stummel, F. The generalized patch test. SIAM J. Numer. Anal. 1979, 16, 449–471. [Google Scholar] [CrossRef]
  26. Shi, Z.C. The F-E-M-test for convergence of nonconforming finite element. Math. Comp. 1987, 49, 391–405. [Google Scholar]
  27. Shi, D.Y. Studies on Nonconforming Elements. Ph.D. Thesis, Xi’an Jiaotong University, Xi’an, China, 1997. [Google Scholar]
  28. Mao, S.P.; Chen, S.C. A quadrilateral nonconforming finite element for linear elasticity problem. Adv. Comput. Math. 2008, 28, 81–100. [Google Scholar] [CrossRef]
  29. Chen, S.C.; Zhao, Y.C.; Shi, D.Y. Non C0 nonconforming elements for elliptic fourth order singular perpurbation problem. J. Comput. Math. 2005, 23, 185–198. [Google Scholar]
  30. Xie, P.L.; Shi, D.Y.; Li, H. A new robust C0-type nonconforming triangular element for singular perturbation problems. Appl. Math. Comput. 2010, 217, 3832–3843. [Google Scholar] [CrossRef]
  31. Chen, S.C.; Liu, M.F.; Qiao, Z.H. An anisotropic nonconforming element for fourth order elliptic singular perturbation problem. Int. J. Numer. Anal. Model. 2010, 7, 766–784. [Google Scholar]
  32. Ciarlet, P.G. The Finite Element Method for Elliptic Problems; North-Holland: Amsterdam, The Netherlands, 1978. [Google Scholar]
  33. Brenner, S.C.; Scott, L.R. Mathematics Theory Finite Element Methods; Springer: New York, NY, USA, 1994. [Google Scholar]
Figure 1. The numerical solution u h with κ = 10 .
Figure 1. The numerical solution u h with κ = 10 .
Mathematics 08 02038 g001
Figure 2. The numerical solution u h with κ = 10 3 .
Figure 2. The numerical solution u h with κ = 10 3 .
Mathematics 08 02038 g002
Figure 3. The numerical solution u h with κ = 10 4 .
Figure 3. The numerical solution u h with κ = 10 4 .
Mathematics 08 02038 g003
Figure 4. Errors in the energy norm with different k.
Figure 4. Errors in the energy norm with different k.
Mathematics 08 02038 g004
Table 1. Errors of u u h h and the CPU time.
Table 1. Errors of u u h h and the CPU time.
κ N × N 4 × 4 8 × 8 16 × 16 32 × 32 64 × 64
κ = 10 0 13.500990360726113.603349640630920.887029283354530.208879901833670.04536489180482
t i m e ( s ) 0.0484820.0737260.2070290.66716951.440274
κ = 10 1 13.483451256013873.595354822628660.884516022038010.208289478802330.04523833745855
t i m e ( s ) 0.1084800.3947702.42188110.93627250.397260
κ = 10 3 11.445822593614402.882588358691550.698854362057020.166104000893760.03616945545106
t i m e ( s ) 5.30754728.444150187.683438793.3383763411.901540
κ = 10 4 11.2412620376426602.873232163905850.624999323349380.150262002401870.03545375844482
t i m e ( s ) 17.87694699.593113919.9725095191.45773719,710.033405
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shi, D.; Pei, L. New High Accuracy Analysis of a Double Set Parameter Nonconforming Element for the Clamped Kirchhoff Plate Unilaterally Constrained by an Elastic Obstacle. Mathematics 2020, 8, 2038. https://doi.org/10.3390/math8112038

AMA Style

Shi D, Pei L. New High Accuracy Analysis of a Double Set Parameter Nonconforming Element for the Clamped Kirchhoff Plate Unilaterally Constrained by an Elastic Obstacle. Mathematics. 2020; 8(11):2038. https://doi.org/10.3390/math8112038

Chicago/Turabian Style

Shi, Dongyang, and Lifang Pei. 2020. "New High Accuracy Analysis of a Double Set Parameter Nonconforming Element for the Clamped Kirchhoff Plate Unilaterally Constrained by an Elastic Obstacle" Mathematics 8, no. 11: 2038. https://doi.org/10.3390/math8112038

APA Style

Shi, D., & Pei, L. (2020). New High Accuracy Analysis of a Double Set Parameter Nonconforming Element for the Clamped Kirchhoff Plate Unilaterally Constrained by an Elastic Obstacle. Mathematics, 8(11), 2038. https://doi.org/10.3390/math8112038

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop