arxiv: 2604.26701 · v1 · submitted 2026-04-29 · 🧮 math.NA · cs.NA

Recognition: unknown

Explicit Planar Finite Element Elasticity Complexes and C¹ Elements on Barycentric Refinements

Chunyu Chen , Long Chen , Xuehai Huang

Authors on Pith no claims yet

Pith reviewed 2026-05-07 12:41 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords finite element elasticity complexbarycentric refinementsC1 finite elementsAiry potentialsmixed finite elementsplane elasticityHsieh-Clough-Tocher elementsexact sequences

0 comments

The pith

Enriching polynomial stresses with three explicit locally supported functions produces an exact finite element elasticity complex and new C1 elements on barycentric refinements.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper seeks to render the Arnold-Douglas-Gupta family of mixed finite elements for plane elasticity fully explicit by deriving closed-form expressions for the enrichment functions added to polynomial stresses on barycentric refined triangles. These enrichments are generated by specific Airy potentials, leading to an explicit C1 potential space of Hsieh-Clough-Tocher type that maps precisely to the desired symmetric stress space. Enforcing single-valued degrees of freedom across elements then assembles global spaces that form a complete finite element complex preserving the exact sequence on simply connected domains. A reader would care because this turns an abstract construction into implementable elements for higher-order accurate simulations of elastic deformations, including new quadratic through higher-order C1 elements.

Core claim

On each macro triangle the symmetric stress space is enriched by three locally supported functions whose closed-form formulas and Airy potentials are derived explicitly. The resulting potential space is a concrete Hsieh-Clough-Tocher type C1 space whose Airy image coincides with the Arnold-Douglas-Gupta stress space. Single-valued degrees of freedom produce global spaces and the full explicit finite element elasticity complex on simply connected domains, from which new families of C1 finite elements on barycentric refinements are obtained for quadratic and higher degrees.

What carries the argument

The three locally supported enrichment functions for the symmetric stress space on macro triangles, generated by explicit Airy potentials that ensure the exact sequence property under single-valued assembly.

Load-bearing premise

The three locally supported enrichment functions can be selected such that their single-valued global assembly exactly matches the Arnold-Douglas-Gupta stress space while preserving the exact sequence on simply connected domains.

What would settle it

Computing the dimension of the assembled stress space on a small triangulation and verifying whether it equals the known dimension of the Arnold-Douglas-Gupta space or whether the Airy map from the potential space is onto the stresses.

read the original abstract

The exact-sequence structure behind the Arnold--Douglas--Gupta family of higher-order mixed finite elements for plane elasticity on barycentric refinements is made explicit. On each macro triangle, the symmetric stress space is obtained by enriching polynomial stresses with three locally supported functions. We derive closed-form formulas for these enrichments and identify explicit Airy potentials that generate them. This leads to a concrete Hsieh--Clough--Tocher type $C^1$ potential space whose Airy image is exactly the Arnold--Douglas--Gupta stress space. By enforcing single-valued degrees of freedom, we obtain global spaces and a fully explicit finite element elasticity complex on simply connected domains. As a consequence, we construct a new family of $C^1$ finite elements on barycentric refinements, including quadratic, cubic, quartic, and higher-order elements.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives closed-form formulas for the three enrichment functions in the ADG stress space on barycentric refinements and ties them to an explicit HCT-style C1 potential space.

read the letter

The main point is that the authors derive explicit expressions for the three locally supported enrichment functions that complete the local polynomial symmetric stresses to the Arnold-Douglas-Gupta space on each macro-triangle, and they identify the corresponding Hsieh-Clough-Tocher type C1 space whose Airy potentials produce exactly those stresses. This turns an abstractly defined family into something with concrete bases and potentials that can be written down and coded directly. The byproduct is a new family of C1 elements on barycentric refinements at quadratic and higher orders, which removes one practical barrier for people who need to implement or analyze these elements. That explicitness is the real contribution here and it is useful within mixed finite-element methods for elasticity. The global construction works by adding the enrichments locally and then imposing single-valued degrees of freedom. For the elasticity complex to stay exact on simply connected domains, the enrichments must sit inside the target space, fill the dimension gap without remainder, and the continuity constraints must preserve the kernel and image properties of the Airy map. The abstract states that this holds after assembly, but the proofs must verify the global dimension count and rule out extra rigid modes or loss of surjectivity. If those checks are tight, the rest follows cleanly; any mismatch would break the exact sequence. This paper is for specialists in finite-element exterior calculus and mixed methods who already know the ADG family and want implementable bases rather than dimension arguments. A reader who codes these spaces or extends them to other refinements will get direct value from the formulas. It deserves a serious referee because the closed-form expressions address a concrete gap even if the global exactness step needs close checking in review.

Referee Report

3 major / 3 minor

Summary. The manuscript derives explicit closed-form expressions for the enrichment functions in the Arnold-Douglas-Gupta family of finite elements for plane elasticity on barycentric refinements. It identifies corresponding Airy potentials in a Hsieh-Clough-Tocher type C^1 space and shows that single-valued degree-of-freedom assembly yields global spaces forming an exact elasticity complex on simply connected domains. This construction is used to define a new family of C^1 finite elements of various polynomial degrees.

Significance. The explicit formulas and potential identification are significant contributions, as they make the abstract ADG complex concrete and implementable. This enables the construction of higher-order C^1 elements on refined meshes, which could have applications in plate bending and other problems requiring C^1 continuity. The work is grounded in finite element exterior calculus and provides a clear path from local to global exact sequences.

major comments (3)

[Local enrichment construction] The claim that the three locally supported enrichment functions lie within the target ADG stress space and complete it to the correct dimension requires explicit verification, such as showing they are divergence-free or satisfy the symmetry and polynomial reproduction properties of ADG.
[Global assembly and exactness (around the discussion of single-valued DOFs)] The argument that enforcing single-valued DOFs on the enrichments preserves the exactness of the elasticity complex on simply connected domains needs to address whether the rigid body modes and the surjectivity of the Airy operator are maintained globally. A dimension count for the global stress space would strengthen this.
[C^1 element family] For the new C^1 elements derived from the potentials, the manuscript should confirm that the degrees of freedom ensure C^1 continuity across macro-edges and provide the dimension formula for each order (quadratic, cubic, etc.).

minor comments (3)

The notation for the enrichment functions could be made more consistent throughout the paper.
Include a reference to the original Arnold-Douglas-Gupta paper in the introduction for context.
Some figures illustrating the support of the enrichment functions would aid clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the positive assessment of our work and for the constructive comments, which help clarify the presentation of the explicit constructions and global properties. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

Referee: The claim that the three locally supported enrichment functions lie within the target ADG stress space and complete it to the correct dimension requires explicit verification, such as showing they are divergence-free or satisfy the symmetry and polynomial reproduction properties of ADG.

Authors: We appreciate this suggestion for strengthening the local construction. The manuscript provides closed-form expressions for the three enrichment functions and asserts that they lie in the ADG space and complete its dimension. To address the request for explicit verification, the revised manuscript will include direct computations confirming that each enrichment is symmetric and divergence-free, together with a verification that the enriched space reproduces the required polynomials and matches the dimension of the target ADG stress space on each macro-triangle. revision: yes
Referee: The argument that enforcing single-valued DOFs on the enrichments preserves the exactness of the elasticity complex on simply connected domains needs to address whether the rigid body modes and the surjectivity of the Airy operator are maintained globally. A dimension count for the global stress space would strengthen this.

Authors: We agree that a global dimension count would make the exactness argument more transparent. The single-valued DOF assembly ensures that the global Airy potentials remain in the C^1 space, thereby preserving the kernel (rigid-body modes) and the surjectivity of the Airy map on simply connected domains. In the revision we will add an explicit global dimension count for the stress space that confirms the dimensions are consistent with exactness of the complex, including correct accounting for the rigid-body modes in the displacement space. revision: yes
Referee: For the new C^1 elements derived from the potentials, the manuscript should confirm that the degrees of freedom ensure C^1 continuity across macro-edges and provide the dimension formula for each order (quadratic, cubic, etc.).

Authors: We thank the referee for highlighting this point. The degrees of freedom are those of the underlying Hsieh-Clough-Tocher-type C^1 space; single-valued assembly of these DOFs guarantees C^1 continuity across macro-edges by construction. The revised manuscript will include an explicit statement confirming this continuity property and will provide the dimension formulas for the C^1 elements of quadratic, cubic, quartic, and higher orders on barycentric refinements. revision: yes

Circularity Check

0 steps flagged

Explicit closed-form enrichments and Airy potentials derived independently of the target ADG space.

full rationale

The paper starts from the known Arnold-Douglas-Gupta stress space on barycentric refinements and derives explicit local enrichment functions plus their Airy potentials from an HCT-type C^1 space. Enforcing single-valued DOFs to obtain global spaces is a standard finite-element assembly step that does not reduce the claimed exact-sequence property to a self-definition or fitted input. No load-bearing step collapses to a prior result by the same authors; the construction is self-contained once the three enrichment functions are exhibited in closed form. External benchmarks (dimension counts, exactness on simply-connected domains) remain independently verifiable.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper introduces no new free parameters or invented entities. It relies on standard polynomial spaces, Airy stress functions, and the exact-sequence framework of finite-element exterior calculus, plus the domain assumption of simple connectedness.

axioms (2)

standard math Standard properties of polynomial spaces and exact sequences in finite-element exterior calculus hold for the Arnold-Douglas-Gupta family.
Invoked throughout the construction of enriched stress spaces and their Airy potentials.
domain assumption The computational domain is simply connected.
Required for the global single-valued spaces to form an exact complex.

pith-pipeline@v0.9.0 · 5448 in / 1412 out tokens · 76029 ms · 2026-05-07T12:41:32.004488+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

37 extracted references · 3 canonical work pages

[1]

Ahmad, A

B. Ahmad, A. Alsaedi, F. Brezzi, L. D. Marini, and A. Russo. Equivalent projectors for virtual element methods.Comput. Math. Appl., 66(3):376–391, 2013. 13

2013
[2]

Argyris, I

J. Argyris, I. Fried, and D. Scharpf. The TUBA family of plate elements for the matrix displacement method. Aero. J. Roy. Aero. Soc., 72:701–709, 1968. 7

1968
[3]

D. N. Arnold.Finite element exterior calculus, volume 93 ofCBMS-NSF Regional Conference Series in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2018. 1

2018
[4]

D. N. Arnold, J. Douglas, Jr., and C. P. Gupta. A family of higher order mixed finite element methods for plane elasticity.Numer. Math., 45(1):1–22, 1984. 1

1984
[5]

D. N. Arnold, R. S. Falk, and R. Winther. Finite element exterior calculus, homological techniques, and applications.Acta Numer., 15:1–155, 2006. 1

2006
[6]

D. N. Arnold and K. Hu. Complexes from complexes.Found. Comput. Math., 21(6):1739–1774, 2021. 1

2021
[7]

D. N. Arnold and R. Winther. Mixed finite elements for elasticity.Numer. Math., 92(3):401–419, 2002. 2

2002
[8]

D. N. Arnold and R. Winther. Mixed finite elements for elasticity in the stress-displacement formulation. Contemporary Mathematics, 329:33–42, 2003. 2

2003
[9]

Beir ˜ao da Veiga and G

L. Beir ˜ao da Veiga and G. Manzini. A virtual element method with arbitrary regularity.IMA J. Numer. Anal., 34(2):759–781, 2014. 3, 13

2014
[10]

S. C. Brenner and L.-Y . Sung.C0 interior penalty methods for fourth order elliptic boundary value problems on polygonal domains.J. Sci. Comput., 22/23:83–118, 2005. 15

2005
[11]

Brezzi, J

F. Brezzi, J. Douglas, Jr., R. Dur ´an, and M. Fortin. Mixed finite elements for second order elliptic problems in three variables.Numer. Math., 51(2):237–250, 1987. 5

1987
[12]

Brezzi, J

F. Brezzi, J. Douglas, Jr., and L. D. Marini. Two families of mixed finite elements for second order elliptic problems.Numer. Math., 47(2):217–235, 1985. 5

1985
[13]

Brezzi and L

F. Brezzi and L. D. Marini. Virtual element methods for plate bending problems.Comput. Methods Appl. Mech. Engrg., 253:455–462, 2013. 3, 13

2013
[14]

C. Chen, L. Chen, X. Huang, and H. Wei. Geometric decomposition and efficient implementation of high order face and edge elements.Commun. Comput. Phys., 35(4):1045–1072, 2024. 12

2024
[15]

C. Chen, X. Huang, and H. Wei.H m-conforming virtual elements in arbitrary dimension.SIAM J. Numer. Anal., 60(6):3099–3123, 2022. 3, 13

2022
[16]

L. Chen, J. Hu, X. Huang, and H. Man. Residual-based a posteriori error estimates for symmetric conform- ing mixed finite elements for linear elasticity problems.Sci. China Math., 61(6):973–992, 2018. 2

2018
[17]

Chen and X

L. Chen and X. Huang. Geometric decompositions of the simplicial lattice and smooth finite elements in arbitrary dimension.arXiv preprint arXiv:2111.10712, 2021. 7

work page arXiv 2021
[18]

Chen and X

L. Chen and X. Huang. A finite element elasticity complex in three dimensions.Math. Comp., 91(337):2095–2127, 2022. 1

2095
[19]

Chen and X

L. Chen and X. Huang. Finite elements for div- and divdiv-conforming symmetric tensors in arbitrary di- mension.SIAM J. Numer. Anal., 60(4):1932–1961, 2022. 5 17

1932
[20]

Chen and X

L. Chen and X. Huang.H(div)-conforming finite element tensors with constraints.Results Appl. Math., 23:Paper No. 100494, 33, 2024. 4

2024
[21]

Chen and X

L. Chen and X. Huang. Finite element complexes in two dimensions.Sci. Sin. Math., 55(8):1593–1626,
[22]

Chen and X

L. Chen and X. Huang. Hybridizable symmetric stress elements on the barycentric refinement in arbitrary dimensions.Math. Comp., https://arxiv.org/abs/2501.02691, 2025. 4

work page arXiv 2025
[23]

Chen and X

L. Chen and X. Huang. Complexes from complexes: Finite element complexes in three dimensions.Math. Comp., 95(359):1083–1142, 2026. 1

2026
[24]

S. H. Christiansen, J. Gopalakrishnan, J. Guzm ´an, and K. Hu. A discrete elasticity complex on three- dimensional Alfeld splits.Numer. Math., 156(1):159–204, 2024. 1

2024
[25]

S. H. Christiansen, J. Hu, and K. Hu. Nodal finite element de Rham complexes.Numer. Math., 139(2):411– 446, 2018. 2

2018
[26]

S. H. Christiansen and K. Hu. Finite element systems for vector bundles: elasticity and curvature.Found. Comput. Math., 23(2):545–596, 2023. 1

2023
[27]

R. W. Clough and J. L. Tocher. Finite element stiffness matrices for analysis of plates in bending. InPro- ceedings of the Conference on Matrix Methods in Structural Mechanics, pages 515–545, Wright-Patterson Air Force Base, Ohio, 1965. 2, 3

1965
[28]

J. J. Douglas, T. Dupont, P. Percell, and R. Scott. A family of finite elements with optimal approxima- tion properties for various galerkin methods for 2nd and 4th order problems.RAIRO. Analyse num ´erique, 13(3):227–255, 1979. 3

1979
[29]

Hsieh.Analysis of Plate Bending by the Finite Element Method

J. Hsieh.Analysis of Plate Bending by the Finite Element Method. PhD thesis, University of California, Berkeley, USA, 1962. 2, 3

1962
[30]

J. Hu. Finite element approximations of symmetric tensors on simplicial grids inR n: the higher order case. J. Comput. Math., 33(3):283–296, 2015. 2, 4

2015
[31]

Hu and S

J. Hu and S. Zhang. A family of conforming mixed finite elements for linear elasticity on triangular grids. arXiv preprint arXiv:1406.7457, 2014. 2, 4

work page arXiv 2014
[32]

Huang and X

J. Huang and X. Huang. Local and parallel algorithms for fourth order problems discretized by the Morley- Wang-Xu element method.Numer. Math., 119(4):667–697, 2011. 15

2011
[33]

X. Huang. Nonconforming finite element Stokes complexes in three dimensions.Sci. China Math., 66(8):1879–1902, 2023. 15

1902
[34]

Johnson and B

C. Johnson and B. Mercier. Some equilibrium finite element methods for two-dimensional elasticity prob- lems.Numer. Math., 30(1):103–116, 1978. 2, 14

1978
[35]

T. K. Nilssen, X.-C. Tai, and R. Winther. A robust nonconformingH2-element.Math. Comp., 70(234):489– 505, 2001. 3

2001
[36]

P. Percell. On cubic and quartic Clough-Tocher finite elements.SIAM J. Numer. Anal., 13(1):100–103, 1976. 3

1976
[37]

Raviart and J

P.-A. Raviart and J. M. Thomas. A mixed finite element method for 2nd order elliptic problems. InMathe- matical aspects of finite element methods (Proc. Conf., Consiglio Naz. delle Ricerche (C.N.R.), Rome, 1975), pages 292–315. Lecture Notes in Math., V ol. 606. Springer, Berlin, 1977. 5 SCHOOL OFMATHEMATICALSCIENCES, PEKINGUNIVERSITY, BEIJING100871, CHIN...

1975