A quasi-incompressible and quasi-inextensible element formulation for transversely isotropic materials

Creative Commons License

Dal H.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, vol.117, pp.118-140, 2019 (SCI-Expanded) identifier identifier

  • Publication Type: Article / Article
  • Volume: 117
  • Publication Date: 2019
  • Doi Number: 10.1002/nme.5950
  • Journal Indexes: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Page Numbers: pp.118-140
  • Keywords: anisotropy, hyperelasticity, mixed finite element design, mixed variational principles, quasi-incompressiblity, quasi-inextensibility, MIXED FINITE-ELEMENT, STABILIZATION TECHNIQUE, GRADIENT PLASTICITY, STRAIN METHODS, DEFORMATION, ELASTICITY, IMPLEMENTATION, ELLIPTICITY, DERIVATION
  • Middle East Technical University Affiliated: Yes


The contribution presents a new finite element formulation for quasi-inextensible and quasi-incompressible finite hyperelastic behavior of transeversely isotropic materials and addresses its computational aspects. The material formulation is presented in purely Eulerian setting and based on the additive decomposition of the free energy function into isotropic and anisotropic parts, where the former is further decomposed into isochoric and volumetric parts. For the quasi-incompressible response, the Q1P0 element formulation is outlined briefly, where the pressure-type Lagrange multiplier and its conjugate enter the variational formulation as an extended set of variables. Using the similar argumentation, an extended Hu-Washizu-type mixed variational potential is introduced, where the volume averaged fiber stretch and fiber stress are additional field variables. Within this context, the resulting Euler-Lagrange equations and the element formulation resulting from the extended variational principle are derived. The numerical implementation exploits the underlying variational structure, leading to a canonical symmetric structure. The efficiency of the proposed approached is demonstrated through representative boundary value problems. The superiority of the proposed element formulation over the standard Q1 and Q1P0 element formulation is studied through convergence analyses. The proposed finite element formulation is modular and exhibits very robust performance for fiber reinforced elastomers in the inextensibility limit.