In this thesis, the numerical simulation of nonlinear behavior of hyperelastic materials in large deflections using the well known finite element method, is discussed. Firstly, the mechanics of rubber and rubber-like materials and the fundamentals of continuum mechanics corresponding to the hyperelastic materials are presented. Moreover, the most familiar models introduced in the litriture and used by hyperelasticity researchers are discussed and their theories, pros and cons are expressed. Additionally the implementation procedure for a new hyperelastic constitutive equation in ABAQUS software has been presented. Since, some hyperelastic models are available in this software, the present thesis provides the reader with the possibility of adding any kind of potential energy function to the material library of the software. One of the hyperelastic energy function wich has been recently proposed by Attatrd, has been selected and implemented in ABAQUS. The results for compressible and incompressible cases are in good agreement with those obtained from theory. Subsequently, the code is used to predict the response of an elastomer bushing. The results are compared with experimental results and some other models available in ABAQUS and a good agreement has been observed.