The aim of this paper is to give a simple, introductory presentation of the extension of the Virtual Element Method to the discretization of H(div)-conforming vector fields (or, more generally, of (n − 1) − Cochains). As we shall see, the methods presented here can be seen as extensions of the so-called BDM family to deal with more general element geometries (such as polygons with an almost arbitrary geometry). For the sake of simplicity, we limit ourselves to the 2-dimensional case, with the aim of making the basic philosophy clear. However, we consider an arbitrary degree of accuracy k (the Virtual Element analogue of dealing with polynomials of arbitrary order in the Finite Element Framework).

The present paper proposes and analyzes a general locking free mixed strategy for computing the deformation of incompressible three dimensional structures placed inside flexible membranes. The model involves as in Chapelle and Ferent [Math. Models Methods Appl. Sci.13 (2003) 573–595] a bending dominated shell envelope and a quasi incompressible elastic body. The present work extends an earlier work ofArnold and Brezzi [Math Comp.66 (1997) 1–14] treating the shell part and proposes a global stable finite element approximation by coupling optimal mixed finite element formulations of the different subproblems by mortar techniques. Examples of adequate finite elements are proposed. Convergence results are derived in two steps. First a global inf-sup condition is proved, deducedfrom the local conditions to be satisfied by the finite elements used for the external shell problem, the internal incompressible 3D problem, and the mortar coupling, respectively. Second, the analysis ofArnold and Brezzi [Math. Comp.66 (1997) 1–14] is extended to the present problem and least to convergence results for the full coupled problem, withconstants independent of the problem's small parameters.