arxiv: 2604.17908 · v1 · submitted 2026-04-20 · 🧮 math.NA · cs.NA

Recognition: unknown

A Coupling Method of Mixed and Lagrange Finite Elements for Linear Elasticity Problem

Wei Chen , Jun Hu , Limin Ma , Mingyan Zhang

Authors on Pith no claims yet

Pith reviewed 2026-05-10 04:27 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords linear elasticitymixed finite elementsLagrange finite elementsdomain couplingstress concentrationa priori error estimateswell-posednessfinite element method

0 comments

The pith

A coupled mixed-Lagrange finite element method for linear elasticity stays well-posed and optimally convergent even when the mixed region shrinks to size O(h).

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

The paper develops a finite element scheme that applies conforming mixed elements only in small zones where stress concentrations occur and standard Lagrange elements in the rest of the domain. This hybrid formulation is shown to remain stable at the interface without extra constraints or penalties. The analysis proves that the discrete problem is well-posed and that the error bounds achieve the same optimal rates as a pure mixed method, even if the mixed subregion measures only O(h) in diameter. Numerical tests confirm the predicted convergence and illustrate reduced computational cost compared with using mixed elements everywhere.

Core claim

The central claim is that the coupled mixed-Lagrange formulation for the linear elasticity system admits a unique solution whose error in suitable norms is bounded by the best-approximation error of each element type, and that these optimal a priori estimates continue to hold when the diameter of the mixed-element subdomain is allowed to be as small as O(h).

What carries the argument

The conforming interface coupling between the mixed finite element subregion and the Lagrange subregion, which preserves the inf-sup stability of the mixed formulation without additional constraints.

If this is right

The coupled problem inherits the same stability constant as the pure mixed method restricted to the small subdomain.
Optimal a priori estimates hold simultaneously in the stress and displacement variables across both subregions.
The method yields stress accuracy comparable to a global mixed discretization while using fewer degrees of freedom outside the concentration zone.
Convergence rates remain independent of the precise size of the mixed subregion down to O(h).

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same interface treatment could be applied to other saddle-point problems where local high-accuracy regions are desirable.
An adaptive algorithm that selects the mixed region on the basis of a posteriori indicators would follow naturally from the analysis.
The O(h) tolerance suggests that the method remains robust even when the interface cuts through a few elements.

Load-bearing premise

The interface between the mixed and Lagrange subdomains can be aligned so that conformity and stability are preserved without extra constraints or loss of approximation quality.

What would settle it

A counter-example in which the discrete system loses uniqueness or the observed convergence rate drops below the optimal order when the mixed subregion is reduced to diameter O(h) would falsify the claim.

read the original abstract

This paper proposes a finite element method that couples mixed and Lagrange finite elements to efficiently capture stress concentrations in elasticity problems. The method employs conforming mixed finite elements in regions with stress concentration, while standard Lagrange elements are used elsewhere, achieving a balance between stress accuracy and computational efficiency. The well-posedness of the coupled formulation and optimal a priori error estimates are established, even when the size of the mixed finite element subregion is $O(h)$. Numerical experiments are presented to verify the theoretical convergence rates and to demonstrate the effectiveness and efficiency of the proposed method.

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

This paper shows how to confine mixed finite elements to an O(h) patch around stress concentrations, couple them stably to Lagrange elements outside, and keep optimal convergence.

read the letter

The main contribution is a coupling formulation that puts conforming mixed elements only in a small O(h)-sized region where stresses concentrate and switches to standard Lagrange elements elsewhere. They prove the coupled problem stays well-posed with a global inf-sup condition that does not blow up as the mixed region shrinks, and they derive a priori error estimates that remain optimal in the usual norms. Numerical tests are included to check the rates and show the expected computational savings. That localization result is the part that is not immediate from existing mixed or Lagrange theory, and the analysis appears to go through without extra stabilization or mesh restrictions beyond standard assumptions. The interface transmission conditions are handled so that conformity is preserved, which is the non-obvious step. The experiments confirm the theory on the test problems they chose. The interface proof is the place to look hardest, since any hidden regularity requirement or special alignment between the two subdomains could affect how small the mixed region can actually be in practice. The abstract is brief on the exact mixed space and the test cases, so the full write-up needs to make those choices explicit. This is aimed at people who already work with finite-element methods for elasticity and want a way to improve local stress accuracy without a full mixed discretization. It is solid enough on the central claim to go to referees rather than be rejected at the desk.

Referee Report

0 major / 1 minor

Summary. The paper presents a coupling method for linear elasticity that uses mixed finite elements in subdomains with stress concentrations and Lagrange finite elements in the remaining domain. It proves the well-posedness of this coupled formulation and establishes optimal a priori error estimates that remain valid even when the mixed finite element subregion has size O(h). The theoretical results are supported by numerical experiments demonstrating the convergence rates and the method's efficiency.

Significance. If the results hold, this method offers a practical way to improve computational efficiency in elasticity simulations by limiting the use of mixed elements to small critical regions without compromising accuracy or convergence rates. The O(h) result is particularly noteworthy as it allows the mixed region to be minimal, potentially leading to significant savings in problems with localized features. The combination of rigorous analysis and numerical validation makes this a valuable contribution to finite element methods for elasticity.

minor comments (1)

Consider adding a figure illustrating the interface between the mixed and Lagrange subdomains to clarify the coupling.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive and encouraging review, which accurately captures the main contributions of our work on the coupled mixed-Lagrange finite element method for linear elasticity. We are pleased that the referee recognizes the practical value of allowing the mixed-element subdomain to be as small as O(h) and still achieve optimal error estimates. The recommendation is to accept, and there are no major comments provided in the report. Therefore, no specific points to address point-by-point, and no revisions are needed based on this feedback. We appreciate the validation of our theoretical and numerical results.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper's central claims—well-posedness of the coupled mixed-Lagrange formulation and optimal a priori error estimates even for mixed subregions of diameter O(h)—rest on direct mathematical analysis of the interface transmission conditions, conformity, and global inf-sup stability using standard finite-element techniques. No step reduces by construction to a fitted parameter renamed as a prediction, a self-citation chain, or an ansatz smuggled from prior work by the same authors. The abstract and manuscript structure indicate the results are derived from the discrete spaces and variational formulation without self-referential definitions or load-bearing external citations that merely rename the target result.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard assumptions from mixed finite element theory for elasticity and on the ability to define a stable interface between element types; no new free parameters or invented entities are mentioned.

axioms (2)

domain assumption Well-posedness and approximation properties of conforming mixed finite elements for linear elasticity hold in the subregion
Invoked to establish the coupled system's stability and error bounds.
domain assumption The interface between mixed and Lagrange subdomains preserves conformity and does not degrade the global error estimates
Required for the O(h) subregion result to hold.

pith-pipeline@v0.9.0 · 5387 in / 1213 out tokens · 35636 ms · 2026-05-10T04:27:28.931427+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

12 extracted references · 2 canonical work pages

[1]

Amara and J

[2]M. Amara and J. M. Thomas,Equilibrium finite elements for the linear elastic problem, Numer. Math., 33 (1979), pp. 367–383. [3]A. Ambroziak,Locking effects in the finite element method, Shell structures: Theory and application, 3 (2013), p

1979
[2]

[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 (1984), pp. 1–22. [5]D. N. Arnold, R. S. Falk, and R. Winther,Mixed finite element methods for linear elasticity with weakly imposed symmetry, Math. Comp., 76 (2007), pp. 1699–1723. 22WEI CHEN, JUN HU, LIMIN MA, AN...

1984
[3]

[10]S. C. Brenner,Korn’s inequalities for piecewiseH 1 vector fields, Math. Comp., 73 (2004), pp. 1067–1087. [11]S. C. Brenner and L. R. Scott,The mathematical theory of finite element methods, vol. 15, Springer, New York, third ed.,

2004
[4]

[12]S. C. Brenner and L.-Y. Sung,Linear finite element methods for planar linear elasticity, Math. Comp., 59 (1992), pp. 321–338. [13]F. Brezzi and M. Fortin,Mixed and hybrid finite element methods, vol. 15, Springer-Verlag, New York,

1992
[5]

Chen and X

[14]G. Chen and X. Xie,A robust weak Galerkin finite element method for linear elasticity with strong symmetric stresses, Comput. Methods Appl. Math., 16 (2016), pp. 389–408. [15]L. Chen, J. Hu, and X. Huang,Fast auxiliary space preconditioners for linear elasticity in mixed form, Math. Comp., 87 (2018), pp. 1601–1633. [16]W. Chen, X. Du, and J. Hu,The Hu...

work page arXiv 2016
[6]

Cockburn, J

[19]B. Cockburn, J. Gopalakrishnan, and J. Guzm ´an,A new elasticity element made for en- forcing weak stress symmetry, Math. Comp., 79 (2010), pp. 1331–1349. [20]R. Courant,Variational methods for the solution of problems of equilibrium and vibrations, Bull. Amer. Math. Soc., 49 (1943), pp. 1–23. [21]P. R. B. Devloo, A. M. Farias, S. M. Gomes, and J. a. ...

2010
[7]

[26]G. Fu, B. Cockburn, and H. Stolarski,Analysis of an HDG method for linear elasticity, Internat. J. Numer. Methods Engrg., 102 (2015), pp. 551–575. [27]V. Girault and P.-A. Raviart,Finite element methods for Navier-Stokes equations, vol. 5, Springer-Verlag, Berlin,

2015
[8]

Gopalakrishnan and J

[29]J. Gopalakrishnan and J. Guzm ´an,Symmetric nonconforming mixed finite elements for linear elasticity, SIAM J. Numer. Anal., 49 (2011), pp. 1504–1520. [30]C. O. Horgan,Korn’s inequalities and their applications in continuum mechanics, SIAM Rev., 37 (1995), pp. 491–511. [31]J. Hu,Finite element approximations of symmetric tensors on simplicial grids in...

2011
[9]

Hu and R

[33]J. Hu and R. Ma,Nonconforming mixed finite elements for linear elasticity on simplicial grids, Numer. Methods Partial Differ. Equations, 35 (2019), pp. 716–732. [34]J. Hu and R. Ma,Partial relaxation ofC 0 vertex continuity of stresses of conforming mixed finite elements for the elasticity problem, Comput. Methods Appl. Math., 21 (2021), pp. 89–

2019
[10]

Hu and S

[35]J. Hu and Z.-C. Shi,Lower order rectangular nonconforming mixed finite elements for plane elasticity, SIAM J. Numer. Anal., 46 (2007), pp. 88–102. A COUPLING METHOD OF MIXED AND LAGRANGE ELEMENTS23 [36]J. Hu and S. Zhang,A family of conforming mixed finite elements for linear elasticity on triangular grids, arXiv:1406.7457 [math.NA], (2014). [37]J. Hu...

work page arXiv 2007
[11]

[42]R. H. Macneal and R. L. Harder,A proposed standard set of problems to test finite element accuracy, Finite Elem. Anal. Des., 1 (1985), pp. 3–20. [43]A. H. Niemi, I. Babu ˇska, J. Pitk¨aranta, and L. Demkowicz,Finite element analysis of the girkmann problem using the modern hp-version and the classical h-version, Eng. Comput., 28 (2012), pp. 123–134. [...

1985
[12]

Stenberg,Postprocessing schemes for some mixed finite elements, RAIRO Mod´ el

[51]R. Stenberg,Postprocessing schemes for some mixed finite elements, RAIRO Mod´ el. Math. Anal. Num´ er., 25 (1991), pp. 151–167. [52]B. Szab ´o, I. Babuˇska, J. Pitk¨aranta, and S. Nervi,The problem of verification with refer- ence to the girkmann problem, Eng. Comput., 26 (2010), pp. 171–183. [53]R. Verf ¨urth,A review of a posteriori error estimation...

1991