doi: 10.1093/imamat/hxh000 Mathematical and computational models of incompressible materials subject to shear

Numerical modeling of incompressible nonlinear elastic materials plays an increasing role in computa-tional science and engineering, particularly in the high-fidelity simulation of rubber-like materials and many biological tissues. Our present study focuses on the treatment of the incompressibility constraint in finite-element discretizations for a cube subject to “simple shear”. We demonstrate that this test prob-lem is not easily captured in three-dimensional mathematical and computational models, with challenges related to the incompressibility constraint that are unique to each approach. Specifically, we review the mathematical model, which presupposes the simple shear deformation and requires additional as-sumptions to determine the response. A computational model avoids these difficulties, but gives rise to questions regarding the mismatch between continuum and discrete representations of material incom-pressibility. We consider three distinct finite-element formulations to enforce (near) incompressibility, with particular attention paid to the competing goals of physical fidelity and computational efficiency. We demonstrate that some of these standard approaches fail to resolve incompressibility in a point-wise man-ner, despite it holding in an averaged sense. Numerical results indicate the maximum ratio of deformed to undeformed volumes grows sharply with mesh refinement; this is in contrast to the mathematical model, but occurs in a manner consistent with finite-element convergence theory