寰枢关节带椎动脉有限元建模及牵引治疗模拟-影像医学与核医学专业论文.docxVIP

万方数据 万方数据 寰枢关节带椎动脉有限元建模及牵引治疗模拟 中文摘要 目的：1、建立寰枢关节带椎动脉有限元模型，并验证。 2、研究寰枢关节各向运动时寰枢关节各部位的应力分布。 3、研究寰枢关节旋转运动对椎动脉第 2、4 弯曲的形态影响。 4、通过此模型进行临床牵引治疗模拟，观察寰枢关节各部位的应力分布 及预期危害性。 方法：从我院 PACS 数据库头颈部 CTA 数据中，寰枢关节及双侧椎动脉未见明显 异常的检查者作为研究对象，且要求体位标准、图像清晰、椎动脉显示清楚、并 排除颈部疾病。图像以 DICOM 格式保存，将其导入 Mimics 软件中建立表面几何 模型，将几何模型在 Geomagic Studio 软件中进行网格优化处理，优化后的模型 导入 Hypermesh 中进行体网格及单元划分，生成三维实体，最后实体模型通过 Abaqus 软件给予材料赋值，生成寰枢关节有限元模型，对其设定边界条件、施 加载荷、提交运算，最后观察运算结果。 结果：1、成功建立寰枢关节带椎动脉的有限元模型，模型与三维 CT 图像一致， 模型在 5N、2N、0.5N 载荷下旋转角度依次为 19.5°、8.9°、2.5°，寰椎外侧 关节面第 874 号节点位移依次为 22.1mm、8.9mm、2.2mm，寰椎后弓应力值依次 为 3100-4200 N/mm2、500-1000 N/mm2、200-400 N/mm2。 2、模型在 5N 力载荷下前屈、后伸、侧屈、旋转运动时寰椎后弓应力集中 现 象 明 显 ， 各 向 运 动 时 寰 椎 后 弓 与 两 侧 块 连 接 部 应 力 值 分 别 为 前 屈 4500-5500N/mm2 ，后伸 4300-5400N/mm2 ， 侧 屈 7000-10000N/mm2 ，旋转 3100-4200N/mm2。各椎体横突因无明显应力作用，应力值接近为零。 3、左、右侧椎动脉第 2 弯曲宽度在左、右旋转位及中立位时分别为 7.6mm、 mm、7.5mm，7.6mm、7.3mm、7.1mm；左、右侧椎动脉第 2 弯曲高度在左、右 旋转位及中立位时分别为 3.7 mm、3.2mm、3.8mm，4.8mm、5.3mm、5.5mm；左、 右侧椎动脉第 4 弯曲宽度在左、右旋转位及中立位时分别为 10.2mm、10.1mm、 10.3mm， 9.0mm、9.2mm、9.2mm；左、右侧椎动脉第 4 弯曲高度在左、右旋转位 及中立位时分别为 4.2mm、4.1mm、4.2mm， 3.9mm、4.0mm、3.9mm。左、右椎动 III 脉倾角在左、右旋转位及中立位时分别为-17.7°、+16.9°、-2.9°，+3.8°、 -29.3°、-13. 1°。 4、模型在施以 5N 牵引力时，寰椎前后弓、枢椎齿突、枢椎椎弓根在前屈 30°、15°、中立位、后伸 15°、30°时的应力值分别为 2500-4000 N/mm2、 1000-2000 N/mm2、600-1200 N/mm2、800-1100 N/mm2、900-1200 N/mm2，6000-9000 N/mm2、1500-2500 N/mm2、2000-3000 N/mm2、4700-6500 N/mm2、5000-8000 N/mm2， 4000-7000 N/mm2、1800-3000 N/mm2、1200-1800 N/mm2、200-700 N/mm2、1200-1700 N/mm2，5500-7000 N/mm2、1200-2000 N/mm2、1100-1600 N/mm2、3500-4800 N/mm2、 4500-7000 N/mm2。模型在施以 2N 牵引力时，寰椎前后弓、枢椎齿突、枢椎椎弓 根在前屈 30°、15°、中立位、后伸 15°、30°时的应力值分别为 600-1000N/mm2、 500-800 N/mm2、400-700 N/mm2、200-600 N/mm2、300-700 N/mm2，1500-3000 N/mm2、 600-1300 N/mm2、500-1000 N/mm2、1600-2500 N/mm2、2000-3300 N/mm2，1300-1800 N/mm2、700-1200 N/mm2、500-800 N/mm2、150-500 N/mm2、200-700 N/mm2，1200-2000 N/mm2、500-800 N/mm2、400-700 N/mm2、1300-2000 N/mm2、1800-2800 N/mm2。 结论：1、基于 CTA 数据，寰枢关节带椎动脉的有限元模型外观逼真，几何相似 性好，可应用于生物力学模拟研究。 2、寰枢关节在各向运动时应力集中在寰椎后弓，尤其为后弓与两侧块连 接部，相当于椎动脉沟