The fatigue damage behavior of bone has attracted significant attention in both the mechanical and orthopedic fields. However, due to the complex and hierarchical structure of bone, describing the damage process quantitively or qualitatively is still a significant challenge for researchers in this area. In this study, a nonlinear bi-modulus gradient model was proposed to quantify the neutral axis skewing under fatigue load in a four-point bending test. The digital image correlation technique was used to analyze the tensile and compressive strains during the fatigue process. The results showed that the compressive strain demonstrated an obvious two-stage ascending behavior, whereas the tensile strain revealed a slow upward progression during the fatigue process. Subsequently, a theoretical model was proposed to describe the degradation process of the elastic modulus and the movement of the neutral axis. The changes in the bone properties were determined using the FEM method based on the newly developed model. The results obtained from two different methods exhibited a good degree of consistency. The results obtained in this study are of help in terms of effectively exploring the damage evolution of the bone materials.
Keywords: FEM method; bi-modulus; fatigue; functionally graded beams; neutral axis.