Simultaneous analysis of curved Kirchhoff beams and Kirchhoff--Love shells embedded in bulk domains
Pith reviewed 2026-05-22 11:55 UTC · model grok-4.3
The pith
A mixed-hybrid Bulk Trace FEM lets standard C0 Lagrange elements from the bulk domain model multiple curved Kirchhoff beams and shells at once via level sets.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors propose a mixed-hybrid and higher-order accurate Bulk Trace FEM that enables the use of standard C0-continuous Lagrange elements with dimensionality of the bulk domain for simultaneous modeling of beams and shells on level sets, with the geometries of the structures implied by level sets of a scalar function over the bulk domain.
What carries the argument
Mixed-hybrid Bulk Trace FEM that approximates displacements and moments on level-set traces using bulk-domain Lagrange elements of arbitrary polynomial degree.
If this is right
- Standard C0-continuous Lagrange elements from the bulk domain can be used for higher-order accurate solutions of Kirchhoff beam and shell models.
- Multiple beams and shells lying on different level sets can be analyzed within a single bulk computation without generating separate surface meshes.
- The approach achieves optimal convergence rates in displacement and stress quantities as verified through numerical tests.
- Level-set geometry descriptions allow parametric variations of the embedded structures by simply modifying the scalar level-set function.
Where Pith is reading between the lines
- The single-bulk-mesh treatment could simplify coupled problems in which beams or shells exchange forces or heat with surrounding material inside the same domain.
- The benchmark cases introduced in the paper provide concrete reference solutions that other embedded or trace finite element schemes can be tested against.
- The mixed-hybrid relaxation of continuity requirements may carry over to time-dependent or moderately nonlinear Kirchhoff-type models while preserving the bulk discretization advantage.
Load-bearing premise
The geometries of the beams and shells are accurately implied by level sets of a scalar function over the bulk domain, and the Kirchhoff-Love assumptions of small displacements without shear deformations hold for all structures simultaneously.
What would settle it
A benchmark test with a known exact solution for a curved Kirchhoff shell in which the measured L2 or energy-norm error fails to decrease at the optimal rate for the chosen polynomial degree would falsify the higher-order accuracy of the mixed-hybrid formulation.
Figures
read the original abstract
A set of curved beams and shells is geometrically implied by level sets of a scalar function over some bulk domain. The mechanical model for each structure is based on the Kirchhoff--Love theory, that is, small displacements without shear deformations are considered. These models for individual geometries are extended to bulk models, simultaneously modeling the whole set of beams/shells on all level sets. A major focus is on the numerical analysis of such models. A mixed-hybrid and higher-order accurate Bulk Trace FEM is proposed that enables the use of standard $C^0$-continuous Lagrange elements with dimensionality of the bulk domain. That is, the higher-order continuity requirements of displacement-based formulations in context of the Kirchhoff--Love theory are successfully alleviated. Several numerical tests confirm the accuracy and higher-order convergence of the proposed methodology, also qualifying as benchmark test cases in future studies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a mixed-hybrid Bulk Trace FEM for simultaneous modeling of curved Kirchhoff beams and Kirchhoff-Love shells geometrically implied by level sets of a scalar function in a bulk domain. Individual Kirchhoff-Love models (small displacements, no shear) are extended to bulk models, and the method uses standard C0-continuous Lagrange elements of bulk dimensionality to alleviate higher-order continuity requirements. Numerical tests are claimed to confirm accuracy and higher-order convergence, positioning the work as providing benchmark cases.
Significance. If the stability and convergence claims hold, the approach would enable efficient implicit embedding of mixed 1D/2D thin structures in bulk meshes without specialized high-continuity elements, offering a practical tool for structural analysis with complex level-set geometries.
major comments (1)
- [Numerical analysis / stability discussion] The central claim of higher-order convergence with standard C0 Lagrange bulk elements requires uniform stability of the mixed-hybrid spaces (displacements plus hybrid multipliers for moments/rotations). No explicit inf-sup analysis or numerical check of the inf-sup constant is reported under mesh refinement or varying level-set curvatures/gradients, leaving the stability (and thus the convergence) conditional on unverified assumptions about the trace operators and hybridization.
minor comments (1)
- [Abstract] The abstract states that 'several numerical tests confirm the accuracy and higher-order convergence' but provides no details on the specific test geometries, curvature ranges, or error measures used; this should be expanded for reproducibility.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback on our manuscript. We address the major comment point by point below, providing clarifications and indicating planned revisions where appropriate.
read point-by-point responses
-
Referee: The central claim of higher-order convergence with standard C0 Lagrange bulk elements requires uniform stability of the mixed-hybrid spaces (displacements plus hybrid multipliers for moments/rotations). No explicit inf-sup analysis or numerical check of the inf-sup constant is reported under mesh refinement or varying level-set curvatures/gradients, leaving the stability (and thus the convergence) conditional on unverified assumptions about the trace operators and hybridization.
Authors: We acknowledge that an explicit inf-sup analysis would provide a stronger theoretical basis. The current work emphasizes the formulation of the mixed-hybrid Bulk Trace FEM and its practical numerical performance rather than a complete stability theory for arbitrary level sets. Section 5 reports multiple benchmark tests demonstrating optimal higher-order convergence for beams and shells under mesh refinement, including examples with curved geometries and different level-set gradients; these results show no degradation indicative of instability. To directly address the concern, the revised manuscript will add a subsection with numerical computations of the inf-sup constant on successively refined meshes for representative level-set configurations. revision: partial
Circularity Check
No circularity: Bulk Trace FEM derivation relies on standard FEM theory and numerical verification
full rationale
The paper extends Kirchhoff-Love beam/shell models to bulk domains via level-set geometry and proposes a mixed-hybrid discretization using standard C0 Lagrange elements. Central claims of higher-order accuracy and simultaneous modeling are supported by explicit numerical convergence tests on benchmark cases rather than by any self-referential definitions, fitted parameters renamed as predictions, or load-bearing self-citations. The derivation chain begins from classical thin-structure kinematics and trace operators, then applies hybridization to relax continuity requirements; these steps are independent of the target results and do not reduce to the inputs by construction. External benchmarks (convergence rates) are used to qualify the method, confirming self-contained analysis.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Kirchhoff-Love theory applies: small displacements without shear deformations for all structures on level sets.
- domain assumption Geometries of beams and shells are exactly implied by level sets of a scalar function over the bulk domain.
Reference graph
Works this paper leans on
-
[1]
Arnold, D.N.; Brezzi, F.: Mixed and nonconforming finite element methods: imple- mentation, postprocessing and error estimates.ESAIM Math. Model. Numer. Anal., 19, 7–32, 1985
work page 1985
-
[2]
Vieweg+Teubner Verlag, Braunschweig, 1985
Başar, Y.; Krätzig, W.B.:Mechanik der Flächentragwerke. Vieweg+Teubner Verlag, Braunschweig, 1985. REFERENCES33
work page 1985
-
[3]
InEncyclopedia of Computational Mechanics
Bischoff, M.; Ramm, E.; Irslinger, J.: Models and Finite Elements for Thin-Walled Structures. InEncyclopedia of Computational Mechanics. (Stein, E.; Borst, R.; Hughes, T.J.; Hughes, T.J., Eds.), John Wiley & Sons, Chichester, 2 edition, 2017
work page 2017
-
[4]
44,Springer Series in Computational Mathematics
Boffi, D.; Brezzi, F.; Fortin, M.:Mixed Finite Element Methods and Applications, Vol. 44,Springer Series in Computational Mathematics. Springer, Berlin, 2013
work page 2013
-
[5]
Brezzi, F.: On the existence, uniqueness and approximation of saddle-point problems arising from Lagrange multipliers.RAIRO Anal. Numér.,8, 129–151, 1974
work page 1974
-
[6]
Burger, M.: Finite element approximation of elliptic partial differential equations on implicit surfaces.Comput. Vis. Sci.,12, 87–100, 2009
work page 2009
-
[7]
Burman, E.; Hansbo, P.; Larson, M.G.: A stabilized cut finite element method for partial differential equations on surfaces: The Laplace-Beltrami operator.Comput. Methods Appl. Mech. Eng.,285, 188–207, 2015
work page 2015
-
[8]
Burman, E.; Hansbo, P.; Larson, M.G.; Massing, A.: Cut finite element methods for partial differential equations on embedded manifolds of arbitrary codimensions. ESAIM Math. Model. Numer. Anal.,52, 2247–2282, 2018
work page 2018
-
[9]
Cenanovic, M.; Hansbo, P.; Larson, M.G.: Cut finite element modeling of linear membranes.Comput. Methods Appl. Mech. Eng.,310, 98–111, 2016
work page 2016
-
[10]
Chapelle, D.; Bathe, K.J.: Fundamental considerations for the finite element analysis of shell structures.Comput. Struct.,66, 19–36, 1998
work page 1998
-
[11]
Cockburn, B.; Gopalakrishnan, J.: A Characterization of Hybridized Mixed Methods for Second Order Elliptic Problems.SIAM J. Numer. Anal.,42, 283–301, 2004
work page 2004
-
[12]
Cockburn, B.; Gopalakrishnan, J.; Lazarov, R.: Unified Hybridization of Discontin- uous Galerkin, Mixed, and Continuous Galerkin Methods for Second Order Elliptic Problems.SIAM J. Numer. Anal.,47, 1319–1365, 2009
work page 2009
-
[13]
Comodi, M.I.: The Hellan-Herrmann-Johnson Method: some new error estimates and postprocessing.Math. Comput.,52, 17–29, 1989
work page 1989
-
[14]
Delfour, M.C.; Zolésio, J.P.: A Boundary Differential Equation for Thin Shells.J. Differ. Equ.,119, 426–449, 1995. 34 REFERENCES
work page 1995
-
[15]
Delfour, M.C.; Zolésio, J.P.: Tangential Differential Equations for Dynamical Thin/Shallow Shells.J. Differ. Equ.,128, 125–167, 1996
work page 1996
-
[16]
Advances in Design and Control
Delfour, M.C.; Zolésio, J.P.:Shapes and geometries: Metrics, Analysis, Differential Calculus, and Optimization. Advances in Design and Control. SIAM, Philadelphia, PA, 2 edition, 2011
work page 2011
-
[17]
Dziuk, G.; Elliott, C.M.: An Eulerian approach to transport and diffusion on evolving implicit surfaces.Comput. Vis. Sci.,13, 17–28, 2008
work page 2008
-
[18]
Dziuk, G.; Elliott, C.M.: Eulerian finite element method for parabolic PDEs on im- plicit surfaces.Interfaces Free Bound.,10, 119–138, 2008
work page 2008
-
[19]
Dziuk, G.; Elliott, C.M.: Finite element methods for surface PDEs.Acta Numer.,22, 289–396, 2013
work page 2013
-
[20]
Echter, R.; Oesterle, B.; Bischoff, M.: A hierarchic family of isogeometric shell finite elements.Comput. Methods Appl. Mech. Eng.,254, 170–180, 2013
work page 2013
-
[21]
Federer, H.:Geometric measure theory. Springer, New York, 1969
work page 1969
-
[22]
Fries, T.P.: Higher-order surface FEM for incompressible Navier-Stokes flows on man- ifolds.Int. J. Numer. Methods Fluids,88, 55–78, 2018
work page 2018
-
[23]
Fries, T.P.; Kaiser, M.W.: On the Simultaneous Solution of Structural Membranes on all Level Sets within a Bulk Domain.Comput. Methods Appl. Mech. Eng.,415, 116223, 2023
work page 2023
-
[24]
Fries, T.P.; Kaiser, M.W.: Simultaneous, Dynamical Analysis of Structural Ropes and Membranes on all Level-sets.Proceedings of the 9th European Congress on Com- putational Methods in Applied Sciences and Engineering (ECCOMAS 2024), Lisbon, Portugal, 2024
work page 2024
-
[25]
Fries, T.P.; Neumeyer, J.; Kaiser, M.W.: A new concept for embedding sub-structures via level-sets.Proceedings of the 16th World Congress on Computational Mechanics (WCCM 2024), Vancouver, Canada, 2024
work page 2024
-
[26]
Fries, T.P.; Schöllhammer, D.: A unified finite strain theory for membranes and ropes. Comput. Methods Appl. Mech. Eng.,365, 113031, 2020. REFERENCES35
work page 2020
-
[27]
Gfrerer, M.H.: AC1-continuousTrace-Finite-Cell-Methodforlinearthinshellanalysis on implicitly defined surfaces.Comput. Mech.,67, 679–697, 2021
work page 2021
-
[28]
Gimena, L.; Gimena, F.N.; Gonzaga, P.: Structural analysis of a curved beam element defined in global coordinates.Eng. Struct.,30, 3355–3364, 2008
work page 2008
-
[29]
Hansbo, P.; Larson, M.G.: Finiteelementmodelingofalinearmembraneshellproblem using tangential differential calculus.Comput. Methods Appl. Mech. Eng.,270, 1–14, 2014
work page 2014
-
[30]
Hansbo, P.; Larson, M.G.; Larsson, K.: Variational formulation of curved beams in global coordinates.Comput. Mech.,53, 611–623, 2014
work page 2014
-
[31]
Jankuhn, T.; Olshanskii, M.A.; Reusken, A.: Incompressible fluid problems on em- bedded surfaces: Modeling and variational formulations.Interfaces Free Bound.,20, 353–377, 2018
work page 2018
-
[32]
Johnson, C.: On the convergence of a mixed finite-element method for plate bending problems.Numer. Math.,21, 43–62, 1973
work page 1973
-
[33]
Kabaria, H.; Lew, A.J.; Cockburn, B.: A hybridizable discontinuous Galerkin formu- lation for non-linear elasticity.Comput. Methods Appl. Mech. Eng.,283, 303–329, 2015
work page 2015
-
[34]
Kaiser, M.W.; Fries, T.P.: Curved, linear Kirchhoff beams formulated using tan- gential differential calculus and Lagrange multipliers.Proc. Appl. Math. Mech.,22, e202200042, 2023
work page 2023
-
[35]
Kaiser, M.W.; Fries, T.P.: Simultaneous analysis of continuously embedded Reissner- Mindlin shells in 3D bulk domains.Int. J. Numer. Methods Eng.,125, e7495, 2024
work page 2024
-
[36]
Kaiser, M.W.; Fries, T.P.: Simultaneous solution of implicitly defined curved, linear Timoshenko beams in two-dimensional bulk domains.Proc. Appl. Math. Mech.,24, e202400019, 2024
work page 2024
-
[37]
Kaiser, M.W.; Fries, T.P.: SimultaneoussolutionofincompressibleNavier-Stokesflows on multiple surfaces.Arch. Appl. Mech.,95, 2025. 36 REFERENCES
work page 2025
-
[38]
Kiendl, J.; Bletzinger, K.U.; Linhard, J.; Wüchner, R.: Isogeometric shell analysis with Kirchhoff-Love elements.Comput. Methods Appl. Mech. Eng.,198, 3902–3914, 2009
work page 2009
-
[39]
Kirchhoff, G.: Über das Gleichgewicht und die Bewegung einer elastischen Scheibe. J. reine angew. Math.,40, 51–88, 1850
-
[40]
The small free vibrations and deformation of a thin elastic shell
Love, A.E.H.: XVI. The small free vibrations and deformation of a thin elastic shell. Philos. Trans. R. Soc. A,179, 491–546, 1888
-
[41]
Mindlin, R.D.: Influence of Rotatory Inertia and Shear on Flexural Motions of Isotropic, Elastic Plates.ASME J. Appl. Mech.,18, 31–38, 1951
work page 1951
-
[42]
Academic press, San Diego, 1988
Morgan, F.:Geometric measure theory: a beginner’s guide. Academic press, San Diego, 1988
work page 1988
-
[43]
Neumeyer, J.; Kaiser, M.W.; Fries, T.P.: Higher-Order, Mixed-HybridFiniteElements for Kirchhoff–Love Shells.Int. J. Numer. Methods Eng.,126, e70175, 2025
work page 2025
-
[44]
Neunteufel, M.; Schöberl, J.: The Hellan–Herrmann–Johnson method for nonlinear shells.Comput. Struct.,225, 106109, 2019
work page 2019
-
[45]
Neunteufel, M.; Schöberl, J.: The Hellan–Herrmann–Johnson and TDNNS methods for linear and nonlinear shells.Comput. Struct.,305, 107543, 2024
work page 2024
-
[46]
InGeometrically Unfitted Finite Element Methods and Applications
Olshanskii, M.A.; Reusken, A.: Trace finite element methods for PDEs on surfaces. InGeometrically Unfitted Finite Element Methods and Applications. (Bordas, S.P.A.; Burman, E.; Larson, M.G.; Olshanskii, M.A.; Olshanskii, M.A., Eds.), Vol. 121,Lec- ture notes in computational science and engineering, Springer Nature, Cham, 211–258, 2017
work page 2017
-
[47]
Olshanskii, M.A.; Reusken, A.; Grande, J.: A finite element method for elliptic equa- tions on surfaces.SIAM J. Numer. Anal.,47, 3339–3358, 2009
work page 2009
-
[48]
Olshanskii, M.A.; Reusken, A.; Xu, X.: A stabilized finite element method for advection-diffusion equations on surfaces.IMA J. Numer. Anal.,34, 732–758, 2014
work page 2014
-
[49]
Osher, S.; Fedkiw, R.P.: Level set methods: an overview and some recent results.J. Comput. Phys.,169, 463–502, 2001. REFERENCES37
work page 2001
-
[50]
Osher, S.; Fedkiw, R.P.:Level Set Methods and Dynamic Implicit Surfaces. Springer, Berlin, 2003
work page 2003
-
[51]
Rafetseder, K.; Zulehner, W.: A decomposition result for Kirchhoff plate bending problems and a new discretization approach.SIAM J. Numer. Anal.,56, 1961–1986, 2018
work page 1961
-
[52]
Rafetseder, K.; Zulehner, W.: A new mixed approach to Kirchhoff–Love shells.Com- put. Methods Appl. Mech. Eng.,346, 440–455, 2019
work page 2019
-
[53]
Reissner, E.: The Effect of Transverse Shear Deformation on the Bending of Elastic Plates.ASME J. Appl. Mech.,12, A69–A77, 1945
work page 1945
-
[54]
Schöllhammer, D.; Fries, T.P.: Kirchhoff-Love shell theory based on tangential differ- ential calculus.Comput. Mech.,64, 113–131, 2019
work page 2019
-
[55]
Schöllhammer, D.; Fries, T.P.: Reissner-Mindlin shell theory based on tangential differential calculus.Comput. Methods Appl. Mech. Eng.,352, 172–188, 2019
work page 2019
-
[56]
Schöllhammer, D.; Fries, T.P.: A higher-order Trace finite element method for shells. Int. J. Numer. Methods Eng.,122, 1217–1238, 2021
work page 2021
-
[57]
Cambridge University Press, Cambridge, 2 edition, 1999
Sethian, J.A.:Level Set Methods and Fast Marching Methods. Cambridge University Press, Cambridge, 2 edition, 1999
work page 1999
-
[58]
Timoshenko, S.P.: LXVI. On the correction for shear of the differential equation for transverse vibrations of prismatic bars.Lond. Edinb. Dublin Philos. Mag. J. Sci.,41, 744–746, 1921
work page 1921
-
[59]
Butterworth-Heinemann, Oxford, 7 edition, 2013
Zienkiewicz, O.C.; Taylor, R.L.; Zhu, J.Z.:The Finite Element Method: Its Basis and Fundamentals. Butterworth-Heinemann, Oxford, 7 edition, 2013
work page 2013
-
[60]
ICE Publishing, London, 2 edition, 2017
Zingoni, A.:Shell Structures in Civil and Mechanical Engineering. ICE Publishing, London, 2 edition, 2017
work page 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.