CT-based finite element analysis (FEA) of human bones helps estimate fracture risk in clinical practice by linking bone ash density (ρash) to mechanical parameters. However, phase field models for fracture prediction require the heterogeneous fracture toughness GIc, which can be derived from the critical stress intensity factor KIc, determined through various experimental methods. Due to a lack of standards for determining cortical bone’s KIc, an experimental campaign is presented using 53 cortical specimens from two fresh frozen femurs to investigate whether a correlation exists between KIc and ρash. We investigated various experimental techniques for correlating KIc with ρash. We conducted FEAs employing the phase field method (PFM) to determine the most suitable correlation among the five possible ones stemming from the experimental methods. The ASTM standard using displacement at force application point was found to be the recommended experimental method for the estimation of KIc perpendicular to osteons’ direction KIc[MPam]=1.89ρash[gr/cc]1.88R2=0.5374.The corresponding statistical critical energy release rate bounds were determined: GIc[N/m]=321.94(ρash[gr/cc])1.69×exp(±2SD),with a standard deviation SD=0.30 representing a 95.4% confidence interval. The average GIc resulted in good correlations between the predicted fracture force by PFM-FEA of four representative specimens and experimental fracture forces. The proposed correlations will be used in CT-based PFM FEA to estimate the risk of hip and humeral fractures.
Read full abstract