“…To overcome the numerical issues in solving saddle-point problems, the widely followed approach in computational solid mechanics is to impose the incompressibility constraint weakly using quasi-incompressible approximations. Some popular approaches for simulating quasiincompressible hyperelastic materials are F -bar formulation [6,7], enhanced-strain method [8], average nodal strain formulation [9,10], reduced integration method with hourglass control [11,12], enhanced assumed strain (EAS) methods [13,14], the two-field displacement-pressure formulation [3,4,15], the three-field formulation [16], mixed stabilised formulations [17][18][19][20][21][22][23][24][25], energy-sampling stabilisation [26,27], least-squares formulations [3,[28][29][30][31], F -bar projection methods [32]. Some noteworthy recent contributions to incompressible and quasi-incompressible computational solid mechanics are in [27,[33][34][35][36][37][38][39][40][41][42][43]…”