A Stabilized Finite Element Framework for Incompressible Hyperelastic Materials at Finite Strains: Analysis, Implementation, and Applications

ABSTRACT This paper presents a stabilized finite element method (FEM) for incompressible hyperelastic materials at finite strains, addressing the computational and implementational challenges of traditional inf‐sup‐stable mixed FEMs. By augmenting the weak formulation with Galerkin/Least‐Squares (GLS) stabilization terms, the method enables equal‐order Lagrange elements (), eliminating the need for complex element pairings like Taylor–Hood (). The stabilization parameter is defined as where denotes the finite‐element mesh size and is the shear modulus; this choice, derived via eigenvalue analysis‐ensures robustness against volumetric locking while preserving consistency. The nonlinear system of algebraic equations is solved via a Newton–Raphson scheme with a block‐preconditioned GMRES solver, achieving optimal displacement convergence rates (comparable to Taylor–Hood elements). Validated on benchmarks including Cook's membrane and a pressurized neo‐Hookean cylinder, the method demonstrates optimal displacement convergence and less than 2% error in reaction forces compared to Abaqus. It reduces degrees of freedom by 30%–40% and solve times by 35% relative to Taylor–Hood elements, as shown in weak scaling tests. An industrial case study on seal compression under 4 mm displacement highlights its capability to handle contact mechanics and geometric nonlinearities, with pressure oscillations suppressed below 5%. The framework's open‐source implementation, combined with its accuracy, efficiency, and robustness under mesh distortion (aspect ratios up to 5:1) and near‐incompressibility (), makes it a versatile tool for engineering applications. Key advancements include a theoretically grounded stabilization strategy, computational efficiency gains, and seamless integration into industrial workflows. The open‐source implementation is available at https://urldefense.com/v3/__https://github.com/ujwalwarbhe/‐A‐Stabilized‐Finite‐Element‐Framework‐for‐Incompressible‐Hyperelastic‐Materials__;!!N11eV2iwtfs!rMduZ7bgMZ55fDmnpDY9Q77FwehCIBLV0_0OWq‐NRskvdXNke5hTj4gOMj625K5MdVsvjPSCktULIRtRsQ$ .

Authors

Institutions

Publication Details

Journal
International Journal for Numerical Methods in Engineering
Published
2025-09-03
DOI
https://doi.org/10.1002/nme.70121
Primary Topic
Advanced Numerical Methods in Computational Mathematics
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

A Stabilized Finite Element Framework for Incompressible Hyperelastic Materials at Finite Strains: Analysis, Implementation, and Applications

Ujwal Warbhe
International Journal for Numerical Methods in Engineering
Advanced Numerical Methods in Computational Mathematics
article

A Stabilized Finite Element Framework for Incompressible Hyperelastic Materials at Finite Strains: Analysis, Implementation, and Applications

Ujwal Warbhe
article en

Abstract

ABSTRACT This paper presents a stabilized finite element method (FEM) for incompressible hyperelastic materials at finite strains, addressing the computational and implementational challenges of traditional inf‐sup‐stable mixed FEMs. By augmenting the weak formulation with Galerkin/Least‐Squares (GLS) stabilization terms, the method enables equal‐order Lagrange elements (), eliminating the need for complex element pairings like Taylor–Hood (). The stabilization parameter is defined as where denotes the finite‐element mesh size and is the shear modulus; this choice, derived via eigenvalue analysis‐ensures robustness against volumetric locking while preserving consistency. The nonlinear system of algebraic equations is solved via a Newton–Raphson scheme with a block‐preconditioned GMRES solver, achieving optimal displacement convergence rates (comparable to Taylor–Hood elements). Validated on benchmarks including Cook's membrane and a pressurized neo‐Hookean cylinder, the method demonstrates optimal displacement convergence and less than 2% error in reaction forces compared to Abaqus. It reduces degrees of freedom by 30%–40% and solve times by 35% relative to Taylor–Hood elements, as shown in weak scaling tests. An industrial case study on seal compression under 4 mm displacement highlights its capability to handle contact mechanics and geometric nonlinearities, with pressure oscillations suppressed below 5%. The framework's open‐source implementation, combined with its accuracy, efficiency, and robustness under mesh distortion (aspect ratios up to 5:1) and near‐incompressibility (), makes it a versatile tool for engineering applications. Key advancements include a theoretically grounded stabilization strategy, computational efficiency gains, and seamless integration into industrial workflows. The open‐source implementation is available at https://urldefense.com/v3/__https://github.com/ujwalwarbhe/‐A‐Stabilized‐Finite‐Element‐Framework‐for‐Incompressible‐Hyperelastic‐Materials__;!!N11eV2iwtfs!rMduZ7bgMZ55fDmnpDY9Q77FwehCIBLV0_0OWq‐NRskvdXNke5hTj4gOMj625K5MdVsvjPSCktULIRtRsQ$ .

International Journal for Numerical Methods in EngineeringVol. 126(17)
Indian Institute of Information Technology, Nagpur (IN)
Openalex Percentile: Top 75%
Advanced Numerical Methods in Computational Mathematics
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.