Due to the complicated structure of the elastic fiber network in annulus fibrosus, existing in-silico studies either simplified or just overlooked its distribution pattern. Nonetheless, experimental and simulation results have proven that elastic fibers are of great importance to maintaining the structural integrity of annulus fibrosus and therefore to ensuring the load-bearing ability of intervertebral discs. Such needs call for a fine model. This work aims at developing a biphasic annulus fibrosus model by incorporating the accurate distribution pattern of collagen and elastic fibers. Both the structural parameters and intrinsic mechanical parameters were successfully identified using single lamella and inter-lamella microscopy anatomy and micromechanical testing data. The proposed model was then used to implement finite element simulations on various anterior and posterolateral multi-lamellae annulus fibrosus specimens. In general, simulation results agree well with available experimental and simulation data. On this basis, the effects of elastic fibers on the integrity of annulus fibrosus were further investigated. It was found that elastic fibers significantly influence the free swelling, radial stretching and circumferential shear performances of annulus fibrosus. Nonetheless, no significant effects were found for the circumferential stretching capability. The proposed biphasic model considers for the first time the distribution characteristics of elastic fibers at two scales, including both the principal orientations of all fiber families and the detailed distribution pattern within each family. Better understandings on the functions of collagen and elastic fibers can therefore be realized. To further enhance its prediction capability, the current model can be extended in the future by taking the fiber–matrix interaction as well as progressive damages into consideration.