Mirzaei Majid, Zeinali Ahad, Razmjoo Arash, Nazemi Majid
Department of Mechanical Engineering, Tarbiat Modares University, Tehran, Iran.
J Biomech. 2009 Aug 7;42(11):1584-91. doi: 10.1016/j.jbiomech.2009.04.042. Epub 2009 May 19.
This paper presents an effective patient-specific approach for prediction of failure initiation and growth in human vertebra using the general framework of the quantitative computed tomography (QCT)-based finite element method (FEM). The studies were carried out on 13 vertebrae (lumbar and thoracic), excised from 3 cadavers with the average age of 42 years old. Initially, 4 samples were QCT scanned and the images were directly converted into voxel-based 3D finite element models for linear and nonlinear analyses. The equivalent plastic strains obtained from the nonlinear analyses were used to predict the occurrence of local failures and development of the failure patterns. In the linear analyses, the strain energy density measure was used to identify the critical elements and predict the failure patterns. Subsequently, the samples were destructively tested in uniaxial compression and the experimental load-displacement diagrams were obtained. The plain radiographic images of the tested samples were also examined for observation of the failure patterns. In continuation, the presence of osteolytic defects in vertebrae was simulated by creation of artificial cavities within 9 remaining samples using a computer numerical control (CNC) milling machine. The same protocol was followed for scanning, modeling, and destructive testing of these samples. A strong correlation was found between the predicted and measured strengths. Finally, a typical vertebroplasty treatment was simulated by injection of low-viscosity bone cement within 3 compressed samples. The failure patterns and the associated load levels for these samples were also predicted using the QCT voxel-based FEM.
本文提出了一种有效的针对个体患者的方法,利用基于定量计算机断层扫描(QCT)的有限元方法(FEM)的通用框架来预测人体椎骨中失效的起始和扩展。研究在从3具平均年龄为42岁的尸体上切除的13块椎骨(腰椎和胸椎)上进行。最初,对4个样本进行了QCT扫描,图像被直接转换为基于体素的三维有限元模型,用于线性和非线性分析。从非线性分析中获得的等效塑性应变被用于预测局部失效的发生和失效模式的发展。在线性分析中,应变能密度度量被用于识别关键单元并预测失效模式。随后,对样本进行单轴压缩破坏试验,并获得实验载荷-位移图。还检查了测试样本的X线平片图像以观察失效模式。接着,通过使用计算机数控(CNC)铣床在9个剩余样本中创建人工空洞来模拟椎骨中溶骨性缺损的存在。对这些样本进行扫描、建模和破坏试验时遵循相同的方案。在预测强度和测量强度之间发现了很强的相关性。最后,通过在3个压缩样本中注入低粘度骨水泥来模拟典型的椎体成形术治疗。还使用基于QCT体素的有限元方法预测了这些样本的失效模式和相关的载荷水平。