Reduced cortical bone thickness increases stress and strain in the female femoral diaphysis analyzed by a CT-based finite element method: Implications for the anatomical background of fatigue fracture of the femur

The incidence of hip fractures is increasing in Japan and is high among women older than 70 years. While osteoporosis has been identified as one of the causative factors of fracture, atypical femoral fracture has emerged as a potential complication of bisphosphonate therapy. Atypical femoral fracture is prevalent among Asian women and has been attributed to morphological parameters. Age-related decreases in the morphological parameters of the femoral diaphysis, such as cortical bone thickness, cortical cross-sectional area, and the cortical index, were reported in Japanese women prior to bisphosphonate drugs being approved for treatment. Thus, in the present study, the relationships between biomechanical and morphological parameters were analyzed using a CT-based finite element method. Finite element models were constructed from 44 femurs of Japanese women aged 31–87 years using CT data. Loading conditions were set as the single-leg configuration and biomechanical parameters, maximum and minimum principal stresses, Drucker-Prager equivalent stress, maximum and minimum strains, and strain energy density were calculated in 7 zones from the subtrochanteric region to distal diaphysis. Pearson's correlation coefficient test was performed to investigate relationships with morphological parameters. While absolute stresses gradually decreased from the subtrochanteric region to distal diaphysis, absolute strains markedly declined in the proximal diaphysis and were maintained at the same levels as those in the distal regions. All types of stresses and minimum principal strain in the femoral diaphysis scored higher absolute values in the high-risk group (≥70 years, n = 28) than in the low-risk group (<70 years, n = 16) (p < 0.05). The distribution patterns of equivalent stress and strain energy density were similar to that of Young's modulus, except for the region of the linea aspera. All biomechanical parameters correlated with morphological parameters and correlation efficiencies, with the reciprocal of cortical bone thickness showing the strongest correlation. The present results demonstrated that biomechanical parameters may be predicted by calculating the cortical bone thickness of femurs not treated with bisphosphonates. Furthermore, strain appeared to be repressed at a low level despite differences in stress intensities among the regions by bone remodeling. This remodeling is considered to be regulated by Wolff's law driven by equivalent stress and strain energy densities from the proximal to distal femur. The present results will promote further investigations on the contribution of morphological parameters in the femoral diaphysis to the onset of atypical femoral fracture.