Achieving Material Robustness via Symmetric Stress Finite Element Discretizations

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Brubeck, Pablo, Parker, Charles, Zerbinati, Umberto
Format: Preprint
Veröffentlicht: 2026
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866914583590993920
author Brubeck, Pablo
Parker, Charles
Zerbinati, Umberto
author_facet Brubeck, Pablo
Parker, Charles
Zerbinati, Umberto
contents When discretizing symmetric stress tensors in variational problems arising in continuum mechanics, one has to choose how to enforce the symmetry of the stress tensor: (i) strongly by requiring the discrete tensors to be pointwise symmetric or (ii) weakly by introducing a Lagrange multiplier. For $H(\mathrm{div})$-conforming finite element discretizations of Hellinger--Reissner elasticity and velocity--stress formulations of incompressible flow, where symmetry of the Cauchy stress tensor is tied to the conservation of angular momentum, we show that this choice may substantially impact the accuracy of the numerical scheme. Through a series of benchmark problems featuring anisotropic constitutive laws inspired by fiber reinforced material, liquid crystal polymer networks, and polar fluids, we show that schemes enforcing symmetry weakly can yield arbitrarily poor stress approximations -- even for zero-stress configurations. However, schemes enforcing symmetry strongly deliver accurate stress approximations independently of the constitutive law, a property we term material robustness. We present a unifying theory that rigorously explains this behavior.
format Preprint
id arxiv_https___arxiv_org_abs_2605_21425
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Achieving Material Robustness via Symmetric Stress Finite Element Discretizations
Brubeck, Pablo
Parker, Charles
Zerbinati, Umberto
Numerical Analysis
65N30 (Primary) 74S05, 76M10 (Secondary)
When discretizing symmetric stress tensors in variational problems arising in continuum mechanics, one has to choose how to enforce the symmetry of the stress tensor: (i) strongly by requiring the discrete tensors to be pointwise symmetric or (ii) weakly by introducing a Lagrange multiplier. For $H(\mathrm{div})$-conforming finite element discretizations of Hellinger--Reissner elasticity and velocity--stress formulations of incompressible flow, where symmetry of the Cauchy stress tensor is tied to the conservation of angular momentum, we show that this choice may substantially impact the accuracy of the numerical scheme. Through a series of benchmark problems featuring anisotropic constitutive laws inspired by fiber reinforced material, liquid crystal polymer networks, and polar fluids, we show that schemes enforcing symmetry weakly can yield arbitrarily poor stress approximations -- even for zero-stress configurations. However, schemes enforcing symmetry strongly deliver accurate stress approximations independently of the constitutive law, a property we term material robustness. We present a unifying theory that rigorously explains this behavior.
title Achieving Material Robustness via Symmetric Stress Finite Element Discretizations
topic Numerical Analysis
65N30 (Primary) 74S05, 76M10 (Secondary)
url https://arxiv.org/abs/2605.21425