An Abstract Stabilization Method with Applications to Nonlinear Incompressible Elasticity

by   Qingguo Hong, et al.
Penn State University
NetEase, Inc

In this paper, we propose and analyze an abstract stabilized mixed finite element framework that can be applied to nonlinear incompressible elasticity problems. In the abstract stabilized framework, we prove that any mixed finite element method that satisfies the discrete inf-sup condition can be modified so that it is stable and optimal convergent as long as the mixed continuous problem is stable. Furthermore, we apply the abstract stabilized framework to nonlinear incompressible elasticity problems and present numerical experiments to verify the theoretical results.



page 1

page 2

page 3

page 4


Fully discrete finite element approximation for a family of degenerate parabolic mixed equations

The aim of this work is to show an abstract framework to analyze the num...

An extended mixed finite element method for elliptic interface problems

In this paper, we propose an extended mixed finite element method for el...

Conforming and Nonconforming Finite Element Methods for Biharmonic Inverse Source Problem

This paper deals with the numerical approximation of the biharmonic inve...

Functional sets with typed symbols: Framework and mixed Polynotopes for hybrid nonlinear reachability and filtering

Verification and synthesis of Cyber-Physical Systems (CPS) are challengi...

Mixed finite element approximation of periodic Hamilton–Jacobi–Bellman problems with application to numerical homogenization

In the first part of the paper, we propose and rigorously analyze a mixe...

Analysis of an abstract mixed formulation for viscoelastic problems

This study provides an abstract framework to analyze mixed formulations ...

A mixed method for 3D nonlinear elasticity using finite element exterior calculus

This article discusses a mixed FE technique for 3D nonlinear elasticity ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Many finite element schemes perform very well in terms of both accuracy and stability for linear elastic problems (see [1, 2, 3, 4]), including in highly constrained situations such as incompressible cases. However, it is well-known that though such schemes can be extended to nonlinear elasticity cases (for instance, large deformation elasticity problems) the same level of stability is by no means guaranteed (see [5, 6, 7, 8]).

In [9, 10], a stabilized discontinuous Galerkin method for nonlinear compressible elasticity was introduced. Because the stabilization term adapts to the solution of the problem by locally changing the size of a penalty term on the appearance of discontinuities, it is called an adaptive stabilization strategy. This same adaptive stabilization strategy can also be found in [11]. Although all three papers discuss compressible and nearly incompressible nonlinear elasticity problems, they do not address the exactly incompressible problems.

In [12], an isogeometric “stream function” formulation was used to exactly enforce the linearized incompressibility constraint. In fact, the exact satisfaction of the linearized incompressibility constraint leads to a numerical approximation of the stability range that converges to the exact one (i.e., the one for the continuum problem). As the stream function is given by a fourth-order PDE, the high regularity of the NURBS shape functions must be used there and the stream function in 3dimensions is not unique.

Brezzi, Fortin, and Marini [13] presented a modified mixed formulation for second-order elliptic equations and linear elasticity problems that automatically satisfied the coercivity condition on the discrete level. Using the same technique and applying stable Stokes elements to a modified formulation for Darcy-Stokes-Brinkman models, Xie, Xu and Xue [14] obtained a uniformly stable method.

In this paper, motivated by the modified formulation used in [13, 14] for second-order equations and the Darcy-Stokes-Brinkman models respectively, we devise a stabilization strategy for the classical mixed finite method designed for Stokes equations and obtain a modified method for nonlinear incompressible elasticity. As indicated in [5], a usual mixed finite element method that is stable for Stokes equations such as the MINI element is not stable for nonlinear incompressible elasticity even though the continuous problem is stable. We note that, for the usual mixed finite method, the unphysical instability is caused by the fact that the discrete solution is not exactly divergence-free. Hence, we rewrite the continuous problem. Based on that we design a stabilization strategy involving the divergence of the discrete solution. We prove that the modified mixed finite element method can remove the unphysical instability and lead to optimal convergence. Hence, all stable Stokes elements can also be made stable for discrete nonlinear elasticity problems.

The rest of the paper is organized as follows. In section 2, we present the finite strain incompressible elasticity problem and obtain the linearized formulation of the nonlinear elasticity problem. In section 3, we propose the stabilization strategy and derive the modified mixed finite element method in an abstract framework, and in section 4 we apply this modified method to the nonlinear incompressible elasticity problem. Furthermore, we obtain an optimal approximation of the modified finite element method. In section 5, we give two simple examples and the corresponding numerical experiments to verify the stability and accuracy of the modified finite element method. We present our conclusions in section 6.

2 Motivation: The finite strain incompressible elasticity problem

To study the finite strain incompressible elasticity problem, we adopt what is known as the material description in this paper. Given a reference configuration for a d-dimensional bounded material body , the deformation of can be described by the map defined by

where  denotes the coordinates of a material point in the reference configuration and

represents the corresponding displacement vector. Following the standard notation (see

[12]), we introduce the deformation gradient

and the right Cauchy–Green deformation tensor

by setting


where I is the second-order identity tensor and is the gradient operator with respect to the coordinate .

For a homogeneous neo-Hookean material, for example, latex and rubber, we define (see [15]) the potential energy function as


where and are positive constants, “:”represents the usual inner product for second-order tensors, and is the deformation gradient Jacobian.

When we introduce he pressure-like variable (or simply pressure) , the potential energy (2.2) can be equivalently written as the following function of and (still denoted as with a little abuse of notation):

When the body is subject to a given load per unit volume in the reference configuration, the total elastic energy functional reads as


According to the standard variational principles, the equilibrium is derived by searching for critical points of (2.3) in suitable admissible displacement and pressure spaces and . The corresponding Euler–Lagrange equations arising from (2.3) lead to this solution: Find  such that


where and are the admissible variation spaces for the displacements and pressures, respectively. Also note that in (2.4), the linearization of the deformation gradient Jacobian is

We now focus on the case of an incompressible material, which corresponds to taking the limit in (2.4). Hence, our problem becomes: Find  such that


or, in residual form: Find  such that




We now derive the linearization of problem (2.5) around a generic point . Observing that

we easily get the problem for the infinitesimal increment : Find  such that



Remark 2.1

Since problem (2.9) is the linearization of problem (2.5), it can be interpreted as the generic step of a Newton-like iteration procedure for the solution of the nonlinear problem (2.5).

Remark 2.2

Taking in (2.9), we immediately recover the classical linear incompressible elasticity problem for small deformations; i.e., we find  such that


where denotes the symmetric gradient operator.

We now note that the Piola identity and give . Hence, we have

Let , we then consider the following problem: Find  such that




Let be the compatible finite element spaces. The discrete problem of (2.11) reads: Find  such that


where  are defined by (2.12).

In some cases, it is possible that the continuous problem (2.11) is stable (or well-posed) but that the discrete problem (2.13) is not (see [5]). In such cases, we say that there is an unphysical instability caused by the discretization.

3 Abstract stabilization strategy for mixed approximation

As shown at the end of the previous section, it is necessary to establish stable discretization when the continuous problem (2.11) is stable. So in this section, we propose an abstract stability strategy framework for the mixed approximation.

Let and each be a Hilbert space, and assume that

are continuous bilinear forms. Let and . We denote both the dual pairing of with and that of with by . We consider the following problem: Find  such that


We associate a mapping with the form :


As is well-known, the continuous problem (3.1) is well-posed under the following assumptions:

  • the inf-sup condition, i.e., it holds that


    where  are norms on the spaces  and , respectively;

  • the coercivity on the kernel space, i.e., it holds that


Now we rewrite the problem (3.1) in an equivalent form as follows: Find  such that


where denotes the inner product on , and is a parameter to be chosen properly.

Let and be finite dimensional subspaces of and such that the following assumption is satisfied:

  •  discrete inf-sup condition holds that



where is independent of .

We consider the discretization of the problem (3.4) as follows: Find  such that

Remark 3.1

Note that when , (3.6) becomes the classical mixed finite element method. We call (3.6) the modified mixed finite element.

Lemma 3.1

(see [16]) is equivalent to the following condition, the operator is an isomorphism, and

Theorem 3.1

Under the assumptions and , there is a constant such that the discrete problem (3.6) is stable for any in the sense that there exists a positive constant such that

  • For any , noting that , we have , where and denotes the orthogonal of  in . By Lemma 3.1, we have

    By assumption , the continuity of and the Cauchy inequality, we obtain that

    Noting that and choosing , we have

    which completes the proof.  

Furthermore, noting that is continuous, by the classic Brezzi theory for the saddle point problem [16, 17, 18], we have

Theorem 3.2

Let be the solution of the problem (3.4), then under the assumptions and providing that is sufficiently large, the discrete problem (3.6) has a unique solution , which satisfies


where is a constant depending on and .

4 Application to the nonlinear elasticity problem

We now apply the abstract stabilized framework proposed in the previous section to nonlinear incompressible elasticity problems. Suppose that  is connected, we consider the d-dimensional Dirichlet boundary value problem (other boundary value problems are similar) of (2.11), which means that . We choose as finite element spaces such that is satisfied. Furthermore, . We rewrite (2.11) as: Find  such that


where ,  are defined by (2.12),  denotes the inner product on , and is a parameter chosen properly.

The discretization of the problem (4.1) is as follows: Find  such that

Lemma 4.1

([19], Corollary 2.3) We have

Let , where  denotes the duality paring between  and .

Lemma 4.2

([19], Corollary 2.4) Let  be connected. Then

  1. the operator grad is an isomorphism from  onto , and

  2. the operator div is an isomorphism from  onto .

Theorem 4.1

If defined by (2.12) satisfies

the discrete problem (4.2) is stable in the sense that there exists a positive constant such that

  • By the abstract framework, the proof is obvious.  

Also from the abstract framework in the previous section, we can get an approximate accuracy result similar to (3.9).

Theorem 4.2

Let be the solution of the problem (2.11), then under the assumptions and and providing that is sufficiently large, the discrete problem (4.2) has a unique solution , which satisfies


where the constant depends on .

5 Numerical experiment

In this section, we report our numerical experiments in regard to the performance of the stabilization strategy for two simple examples. We consider here simple problems using some mixed finite element formulations that are known to be optimal for Stokes equations. We first briefly present the finite element under consideration and then show the numerical results obtained using such element.

5.1 Two examples

In this subsection, we present two simple problems that will be used in subsection 5.3 to discuss the performance of the stabilization strategy proposed in Section 3. Using the usual Cartesian coordinates , we consider a square material body whose reference configuration is . We denote as the upper part of its boundary, while the remaining part of is denoted with . The total energy is assumed to be as in (2.3), where the external loads are given by the vertical uniform body forces: , where . The two problems differ in regard to the imposed boundary conditions. More precisely:

  • Problem 1. We set clamped boundary conditions on , but traction-free boundary conditions on .

  • Problem 2. We set vanishing normal displacements on , but traction-free boundary conditions on .

It is easy to see that both problems admit a trivial solution for every , where .

For the problems under investigation, the corresponding linearized problems (cf. (2.12)) can both be written as: Find  such that


where is the increment of the parameter of .

For these two different problems, the spaces and are defined as follows:

  • Problem 1. .

  • Problem 2. , where denotes the outward normal vector.

The stable discrete formulation reads as follows (see (4.2)): Find  such that


where ;
. Here is a parameter chosen depending on .

Remark 5.1

For Problem , it has been theoretically proved in [5] that the continuous problem (5.1) is stable when .

5.2 The MINI element

The considered scheme for (5.2) is the MINI element (see [18]). Let be a triangular mesh of with the mesh size . For the discretization of the displacement field, we take

where is the space of linear functions on , and is the linear space generated by , the standard cubic bubble function on . For the pressure discretization, we take

5.3 Numerical results

We now study the stability performance of the discretized model problems by means of the modified mixed finite element formulations briefly described above. It has been theoretically proved and numerically verified in [5] that for Problem the classical mixed finite element method is stable when , but unstable when . In this numerical experiment, we will demonstrate that the modified mixed finite element method is stable for both Problem and Problem when (which is the stability range for the continuous case of Problem in [5]) as predicated by our Theorem 3.1.

Noting Theorem 3.1

, we study the eigenvalues of the matrix induced by the bilinear form

for the problems under consideration. The first loads for which we find a negative eigenvalue are the critical ones. We start from for both positive and negative loading conditions, i.e., for and . We indicate as the critical loads and the first load values for which we find a negative eigenvalue. A subsequent bisection-type procedure is used to increase the accuracy of the critical load detections. The corresponding nondimensional quantities are denoted with and , respectively, where . Here, is some problem characteristic length, set equal to for simplicity, consistents also with the geometry of the model problems. If we do not detect any negative eigenvalue for extremely large values of the load multiplier , we set . Noting that the bound constant for the formulation is of order and the stable constant of is of order , we can choose by the proof of Theorem 3.1. Hence, in the codes, we set , then we adjust the parameter , and to investigate the stability performance. Furthermore, we take for Problem and for Problem .

Table 1: Stability Limits for Problem  
Table 2: Stability Limits for Problem  

To verify the convergence result (3.2), we set another data for (5.1) that and the true solution is .

Table 3: The Convergence of for Problem
Table 4: The convergence of for Problem

Table 1 and Table 2 show that the stabilization strategy together with the corresponding modified mixed finite element method proposed in this paper is effective. The stability performance can be improved obviously. In fact, these values of are competitive to the “exact” values claimed in [12]. Furthermore, we can see that the modified mixed finite element method is also locking-free by verifying the convergence results (3.2), which are shown in Table 3 and Table 4. In fact, the classical mixed finite element method is unstable for the Problem when , and hence the convergence fails. However, the modified method performs well.

6 Conclusions

Within the framework of finite elasticity for incompressible materials, it is well known that the classical mixed finite element discretization can sometimes be unstable even though the continuous problem is stable. In this paper, we reformulated the continuous problem and proposed an abstract stabilization strategy based on the new continuous formulation and obtained a modified mixed finite element method. We proved theoretically in Section 3 that for a sufficiently large , the modified mixed finite element method is stable whenever the continuous problem is stable, and the method maintains the optimal convergence of the classical one. We verified by numerical experiments in Section 5 that the modified mixed finite element method is much more stable than the classical one and is also locking-free. The provided in Section 3

always overestimates the parameter used in practical problems. However, we can choose the parameter heuristically by analyzing the stability and continuity of continuous problems in the numerical experiments presented in Section



  • [1] B. Reddy, J. Simo, Stability and convergence of a class of enhanced strain methods, SIAM J. Numer. Anal. 32 (6) (1995) 1705–1728.
  • [2]

    D. Braess, C. Carstensen, B. Reddy, Uniform convergence and a posteriori error estimators for the enhanced strain finite element method, Numer. Math. 96 (3) (2004) 461–479.

  • [3] P. Houston, D. Schotzau, T. Wihler, An hp-adaptive mixed discontinuous galerkin fem for nearly incompressible linear elasticity, Comput. Methods Appl. Mech. Engrg. 195 (25-28) (2006) 3224–3246.
  • [4] P. Hansbo, M. Larson, Discontinuous Galerkin methods for incompressible and nearly incompressible elasticity by Nitsche’s method, Comput. Methods Appl. Mech. Engrg. 191 (17-18) (2002) 1895–1908.
  • [5] F. Auricchio, L. Beirao da Veiga, C. Lovadina, A. Reali, A stability study of some mixed finite elements for large deformation elasticity problems, Comput. Methods Appl. Mech. Engrg. 194 (9-11) (2005) 1075–1092.
  • [6] P. Wriggers, S. Reese, A note on enhanced strain methods for large deformations, Comput. Methods Appl. Mech. Engrg. 135 (3-4) (1996) 201–209.
  • [7] C. Lovadina, F. Auricchio, On the enhanced strain technique for elasticity problems, Comput.& Structures 81 (8-11) (2003) 777–787.
  • [8] D. Pantuso, K. Bathe, On the stability of mixed finite elements in large strain analysis of incompressible solids, Finite Elem. Anal. Des. 28 (2) (1997) 83–104.
  • [9] A. Ten Eyck, F. Celiker, A. Lew, Adaptive stabilization of discontinuous Galerkin methods for nonlinear elasticity: Analytical estimates, Comput. Methods Appl. Mech. Engrg. 197 (33-40) (2008) 2989–3000.
  • [10] A. Ten Eyck, F. Celiker, A. Lew, Adaptive stabilization of discontinuous Galerkin methods for nonlinear elasticity: Motivation, formulation, and numerical examples, Comput. Methods Appl. Mech. Engrg. 197 (45-48) (2008) 3605–3622.
  • [11] A. Ten Eyck, A.Eyck, A. Lew, An adaptive stabilization strategy for enhanced strain methods in non-linear elasticity, Internat. J. Numer. Methods Engrg. 81 (11) (2010) 1387–1416.
  • [12] F. Auricchio, L. Beirão da Veiga, C. Lovadina, A. Reali, The importance of the exact satisfaction of the incompressibility constraint in nonlinear elasticity: mixed FEMs versus NURBS-based approximations, Comput. Methods Appl. Mech. Engrg. 199 (5-8) (2010) 314–323.
  • [13] F. Brezzi, M. Fortin, L. Marini, Mixed finite element methods with continuous stresses, Math. Models Methods Appl. Sci. 3 (1993) 275–287.
  • [14] X. Xie, J. Xu, G. Xue, Uniformly stable finite element methods for Darcy-Stokes-Brinkman models, J. Comput. Math. 26 (3) (2008) 437–455.
  • [15] J. Bonet, R. Wood, Nonlinear continuum mechanics for finite element analysis, Cambridge University Press, New York, 1997.
  • [16] F. Brezzi, M. Fortin, Mixed and hybrid finite element methods, Springer-Verlag, New York, 1991.
  • [17] F. Brezzi, On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers, RAIRO Anal. Numer. 8 (2) (1974) 129–151.
  • [18] D. Arnold, F. Brezzi, M. Fortin, A stable finite element for the Stokes equations, Calcolo 21 (4) (1984) 337–344.
  • [19] V. Girault, P. Raviart, Finite element methods for Navier-Stokes equations:theory and algorithms, Springer-Verlag, Berlin-Heidelberg-New York-Tokyo, 1986.