The goal of this thesis is to develop a hybrid _nite element formulation to carry out stress analysis of arteries. To the best of our knowledge, a hybrid _nite element impel mentation of the Holzapfel-Ogden artery model has not been carried out before. Since arteries are thin `shell-type' structures, they are subjected to membrane, shear and volumetric locking in case when standard _nite elements are used. Since hybrid _nite elements are known to overcome these problems, we develop hybrid hexahedral element formulations (both lower and higher-order) for artery analysis. We demonstrate The better coarse mesh accuracy of hybrid elements, which are based on a two-_eld variational formulation, over conventional displacement based elements. Typically, wend that three or four extra levels of renement are required with conventional elements to achieve the ame accuracy as hybrid elements. The recently proposed Holzapfel-Ogden constitutive model for the artery and its implementation both within the conventional and hybrid _nite element frameworks is discussed. The numerical implementation is particularly challenging due to the presence of _bers which can only take tensile loads. The mathematically exact tangent stiness matrix that we have derived in this work is crucial in ensuring convergence of the numerical strategy.