Abstract
The distribution of collagen fibers across articular cartilage layers is statistical in nature. Based on the concepts proposed in previous models, we developed a methodology to include the statistically distributed fibers across the cartilage thickness in the commercial FE software COMSOL which avoids extensive routine programming. The model includes many properties that are observed in real cartilage: finite hyperelastic deformation, depth-dependent collagen fiber concentration, depth- and deformation-dependent permeability, and statistically distributed collagen fiber orientation distribution across the cartilage thickness. Numerical tests were performed using confined and unconfined compressions. The model predictions on the depth-dependent strain distributions across the cartilage layer are consistent with the experimental data in the literature.