Abstract

A higher order mixed finite element method is presented for compressible transversely isotropic finite hyperelasticity. The independent variables of the three-field formulation are; displacement, fibre tension and fibre stretch. The formulation admits the description of a fully coupled stress response which evolves from an almost isotropic one into a hyper-anisotropic simply nearly inextensible response with increasing fibre tension, as for example observed in soft tissue biomechanics. Standard displacement approximation of order p in H1(Ω) is used while the auxiliary variables are approximated element wise by square integrable functions of order p−1 in L2(Ω). For finite extensibility the auxiliary variables are statically condensed out yielding a pure displacement based method. It is implemented in an hp-adaptive finite element code. A matching residual based error estimation capability is added and exercised. Numerical evidence indicating stability of approximation is supplied resolving a boundary layer caused by almost inextensible fibres. Coupled and uncoupled stress responses are compared. A generalised compressible transversely isotropic Holzapfel–Gasser–Ogden model is developed for the formulation that intentionally avoids the volumetric–isochoric split. Solutions obtained with the mixed method compare favourably to the corresponding ones obtained with pure displacement formulation. The latter fails in certain strongly anisotropic cases.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.