(Hyperelastic analysis using MATLAB) Consider a unit cube as shown in Fig

(Hyperelastic analysis using MATLAB) Consider a unit cube as shown in Fig. 3.15. An eight-node solid element is used to model the cube. The positive X1 face (Face 4) is extended with a stretch ratio λ=6.0. The following boundary conditions are given: u₁=0 at Face 6, u₂=0 at Face 3, and u₃=0 at Face 1. Using NLFEA, calculate the relation between the stretch ratio and Cauchy stress. Use the incompressible Mooney-Rivlin hyperelastic material with A₁₀=80 MPa, A₀₁=20 MPa, and K=10⁷ MPa. Use 20 load increments. Compare the results withthe analytical solution.