Pakdel Amirreza, Fialkov Jeffrey, Whyne Cari M
Sunnybrook Research Institute, Toronto, ON, Canada; Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, ON, Canada.
Sunnybrook Research Institute, Toronto, ON, Canada; Department of Surgery, University of Toronto, Toronto, ON, Canada.
J Biomech. 2016 Jun 14;49(9):1454-1460. doi: 10.1016/j.jbiomech.2016.03.015. Epub 2016 Mar 16.
Accurate finite element (FE) modeling of complex skeletal anatomy requires high resolution in both meshing and the heterogeneous mapping of material properties onto the generated mesh. This study introduces Node-based elastic Modulus Assignment with Partial-volume correction (NMAP) as a new approach for FE material property assignment to thin bone structures. The NMAP approach incorporates point spread function based deblurring of CT images, partial-volume correction of CT image voxel intensities and anisotropic interpolation and mapping of CT intensity assignment to FE mesh nodes. The NMAP procedure combined with a derived craniomaxillo-facial skeleton (CMFS) specific density-isotropic elastic modulus relationship was applied to produce specimen-specific FE models of 6 cadaveric heads. The NMAP procedure successfully generated models of the complex thin bone structures with surface elastic moduli reflective of cortical bone material properties. The specimen-specific CMFS FE models were able to accurately predict experimental strains measured under in vitro temporalis and masseter muscle loading (r=0.93, slope=1.01, n=5). The strength of this correlation represents a robust validation for CMFS FE modeling that can be used to better understand load transfer in this complex musculoskeletal system. The developed methodology offers a systematic process-flow able to address the complexity of the CMFS that can be further applied to create high-fidelity models of any musculoskeletal anatomy.
复杂骨骼解剖结构的精确有限元(FE)建模需要在网格划分以及将材料属性非均匀映射到生成的网格方面具备高分辨率。本研究引入基于节点的弹性模量分配与部分体积校正(NMAP),作为一种将FE材料属性分配到薄骨结构的新方法。NMAP方法包括基于点扩散函数的CT图像去模糊、CT图像体素强度的部分体积校正以及CT强度分配到FE网格节点的各向异性插值和映射。将NMAP程序与推导得出的颅颌面骨骼(CMFS)特定密度 - 各向同性弹性模量关系相结合,应用于生成6个尸体头部的特定标本FE模型。NMAP程序成功生成了具有反映皮质骨材料属性的表面弹性模量的复杂薄骨结构模型。特定标本的CMFS FE模型能够准确预测在体外颞肌和咬肌加载下测得的实验应变(r = 0.93,斜率 = 1.01,n = 5)。这种相关性的强度代表了对CMFS FE建模的有力验证,可用于更好地理解这个复杂肌肉骨骼系统中的载荷传递。所开发的方法提供了一个系统的流程,能够解决CMFS的复杂性,可进一步应用于创建任何肌肉骨骼解剖结构的高保真模型。