An anisotropic damage model for soft fibered tissue is presented in this paper, using a multi-scale scheme and focusing on the directionally dependent behavior of these materials. For this purpose, a micro-structural or, more precisely, a microsphere-based approach is used to model the contribution of the fibers. The link between micro-structural contribution and macroscopic response is achieved by means of computational homogenization, involving numerical integration over the surface of the unit sphere. In order to deal with the distribution of the fibrils within the fiber, a von Mises probability function is incorporated, and the mechanical (phenomenological) behavior of the fibrils is defined by an exponential-type model. We will restrict ourselves to affine deformations of the network, neglecting any cross-link between fibrils and sliding between fibers and the surrounding ground matrix. Damage in the fiber bundles is introduced through a thermodynamic formulation, which is directly included in the hyperelastic model. When the fibers are stretched far from their natural state, they become damaged. The damage increases gradually due to the progressive failure of the fibrils that make up such a structure. This model has been implemented in a finite element code, and different boundary value problems are solved and discussed herein in order to test the model features. Finally, a clinical application with the material behavior obtained from actual experimental data is also presented.