Extensive research has been performed in the past decades to study the behavior of soft biological tissues in order to reduce the need for practical experiments. The applicability of these researches, particularly for skin, ligament, muscles and the heart, brings up its importance in various biological science and technology disciplines such as surgery and medicine. Softness and large deformation govern the behavior of soft materials and prohibit the use of small strains solutions in finite element method.In this work, the focus is set on a strain energy function which has the advantage of accurately representing the behavior of a variety of soft tissues with only a few parameters in a finite element approach. The numerical results are verified with a set of tensile experiments to demonstrate the performance of the proposed model. The parameters include the matrix and collagen bundles and their orientation. Different cases are analyzed and discussed for better prediction of different soft tissue responses.