Kumar Ajay, Shitole Pankaj, Ghosh Rajesh, Kumar Rajeev, Gupta Arpan
School of Engineering, Indian Institute of Technology Mandi, Mandi, India.
Proc Inst Mech Eng H. 2019 Aug;233(8):823-838. doi: 10.1177/0954411919853918. Epub 2019 Jun 6.
Stress intensity factor and energy release rate are important parameters to understand the fracture behaviour of bone. The objective of this study is to predict stress intensity factor and energy release rate using finite element method, element-free Galerkin method, and extended finite element method and compare these results with the experimentally determined values. For experimental purpose, 20 longitudinally and transversely fractured single-edge notched bend specimens were prepared and tested according to ASTM standard. All specimens were tested using the universal testing machine. For numerical simulations (finite element method, element-free Galerkin method, and extended finite element method), two-dimensional model of cortical bone was developed by assuming plane strain condition. Material properties of the cortical bone were considered as anisotropic and homogeneous. The values obtained through finite element method, element-free Galerkin method, and extended finite element method are well corroborated to experimentally determined values and earlier published data. However, element-free Galerkin method and extended finite element method predict more accurate results as compared to finite element method. In the case of the transversely fractured specimen, the values of stress intensity factor and energy release rate were found to be higher as compared to the longitudinally fractured specimen, which shows consistency with earlier published data. This study also indicates element-free Galerkin method and extended finite element method predicted stress intensity factor and energy release rate results are more close to experimental results as compared to finite element method, and therefore, these methods can be used in the different field of biomechanics, particularly to predict bone fracture.
应力强度因子和能量释放率是了解骨骼断裂行为的重要参数。本研究的目的是使用有限元法、无单元伽辽金法和扩展有限元法预测应力强度因子和能量释放率,并将这些结果与实验测定值进行比较。出于实验目的,制备了20个纵向和横向断裂的单边缺口弯曲试样,并按照ASTM标准进行测试。所有试样均使用万能试验机进行测试。对于数值模拟(有限元法、无单元伽辽金法和扩展有限元法),通过假设平面应变条件建立了皮质骨的二维模型。皮质骨的材料特性被视为各向异性和均匀的。通过有限元法、无单元伽辽金法和扩展有限元法获得的值与实验测定值和早期发表的数据得到了很好的证实。然而,与有限元法相比,无单元伽辽金法和扩展有限元法预测的结果更准确。在横向断裂试样的情况下,发现应力强度因子和能量释放率的值比纵向断裂试样更高,这与早期发表的数据一致。本研究还表明,与有限元法相比,无单元伽辽金法和扩展有限元法预测的应力强度因子和能量释放率结果更接近实验结果,因此,这些方法可用于生物力学的不同领域,特别是预测骨折。