Computational analysis of the nonlinear mechanical properties of anisotropic hyperelastic materials aims at a better understanding of its physiology and pathophysiology under different loading conditions. This has an important role in biomechanics, surgical, clinical diagnostic and design of medical devices. This study investigates the modeling of arterial tissues made of a four-fiber family by using an anisotropic hyperelastic model. This model is based on the theory of polynomial invariant and was implemented in the university finite element code FER. The convex property of the strain energy function is investigated as well as the positive definite nature of the tangent stiffness matrix used within the framework of a finite element analysis. This allows us to guarantee the invertibility of the linearized problem and the uniqueness of the solution computed at each step of the Newton–Raphson scheme used to solve nonlinear problems.