本发明涉及计算机辅助技术及脊柱外科治疗领域，更具体地说，涉及一种胸腰椎爆裂性骨折仿真模型建立方法及装置。
背景技术：
胸腰椎爆裂性骨折是脊柱外科常见的损伤性疾病，约占所有脊柱骨折的50％左右，常见于车祸及高处坠落等高能量损伤，患者常合并后凸畸形及骨折块突入椎管，压迫神经根和马尾，造成脊髓神经功能损伤。患者可出现不同程度的感觉障碍和运动功能丧失，更严重的会导致瘫痪，极大的影响了患者的身体和心理健康。解除骨折占位，恢复椎管容积，减轻脊髓压迫获得了外科医生的广泛共识。
随着医疗技术的发展和医疗器械的改善，经椎弓根伤椎置钉治疗胸腰椎爆裂性骨折得到了广泛关注和应用，并取得了良好的临床效果。与保守治疗和直接复位相比，经椎弓根伤椎置钉复位骨折，避免了骨折块复位不完全、出血量大损伤严重等问题，具有独特的优势。
目前对胸腰椎爆裂性骨折复位的生物力学研究主要通过动物脊柱、人的尸体脊柱等的体外实验实现。选取胸腰椎节段，模拟暴力打击脊柱，获得胸腰椎爆裂性骨折模型，从而进行力学分析和手术入路研究，其优点是：能直接、精确测量脊柱各节段的运动，同时可直接观察到脊柱的生理病理变化和组织形态学改变。缺点是：实验方法较为昂贵，标本难以获取、各标本之间存在一定差异性，可重复性差，同时，对于标本受力的测取，主要集中在皮质表面，对于内部的应力和应变，则明显受限。有限元法的优势则在于可以模拟各种复杂工况，通过参数设计，可无限次重复实验。因此，可利用有限元法与体外实验相结合，相互支持补充，以期可以更好的研究胸腰椎爆裂性骨折及相关手术治疗方法。
2016年第三军医大学硕士学位论文“横连接在胸腰椎骨折椎弓根螺钉固定中作用的有限元分析及临床研究”，利用64排ct对志愿者t9-l3节段进行扫描，得到的图像以dicom格式保存。然后用mimics软件读取已保存的dicom图像。根据骨与软组织具有不同灰度值的原理，进一步阈值化及擦除操作处理目标图像，然后剔除骨质旁软组织图像，利用mask的编辑功能填充椎间盘，得到目标节段的几何模型。最后利用rapidform软件对得到的几何模型进行光滑处理。再利用truss单元模拟前、后纵韧带、黄韧带、关节囊，棘间及棘上韧带等，将上述模型导入abaqus6.9.1软件中进行组合得到正常胸腰椎(t9-l3)的有限元模型。然后利用布尔运算法模拟t12椎体前方压缩1/2高度，椎间盘、椎体后方骨性结构及后方韧带复合体均完整，最终建立了t12椎体爆裂骨折的有限元模型。
期刊《中国矫形外科杂志》2015年第23卷第10期，刊出的论文“胸腰椎爆裂性骨折短节段伤椎固定三维有限元模型构建及生物力学比较研究”，建立了胸腰椎爆裂性骨折后路不同固定方式的尸体标本和三维有限元模型，具体取1具尸体为正常胸腰椎模型，用螺旋ct以层厚2mm(0.6mm居多)连续扫描获得横断面图像，dicom格式导入mimics软件，得到标本的三维信息图。手动完成t11-l3椎体模型后，在mimics的fea模块中对模型进行网格划分，得到较为满意的三维图像后导入ansys12.0软件中。然后对各单元进行赋值，包括弹性模量、泊松比等材料常数及特征值。整个模型共有结点675个，不同类型单元860个。
期刊《中国生物医学工程学报》2017年第36卷第6期，刊出的论文“后路固定术治疗胸腰椎爆裂骨折的生物力学研究”，建立了t10-l2节段正常脊柱模型、爆裂性骨折模型。具体地，将ct扫描获得的水平位图像以dicom格式导入医学图像重建软件mimics10.01中，以得到t10-l2节段的三维几何模型；再通过geomagicstudio软件处理，使其成为由nurbs曲面表达的实体结构，用于后续有限元的仿真模拟；然后经过solidworks2012的装配以及在hypermesh中进行六面体网格划分，最后导入abaqus6.13进行分析。每个椎体分为皮质骨和松质骨两部分。每个椎间盘模型包含3个部分：髓核、纤维环基质和纤维环纤维层，纤维环基质从内到外分为7层，分割纤维环的纤维层也分为7层，嵌入在纤维环基质中，胶原纤维这些层的角度设置为30°。此模型包含了7种韧带，这些韧带被设置为线性拉伸弹簧元件。胸椎t12节段采用楔形截骨法切除部分椎骨，保留后部结构以建立不稳定的胸腰椎爆裂骨折的有限元模型。
2018年南京中医药大学博士学位论文“3d打印外固定支具治疗单纯性胸腰椎骨折的有限元分析及临床应用研究”，建立了胸腰椎骨折的三维有限元模型(t12-l4)，具体地，选取一名l1骨折患者，应用螺旋ct对胸腰椎进行薄层连续扫描，扫描范围从胸椎t11至骶骨s1，共获得断层图像序列810张，以dicom格式存储于计算机中。然后导入mimics医学图像处理软件中，利用软件三维重建功能建立t12-l4椎体骨各个椎体模型。再将上述模型导入到3-matic软件，以椎体骨的上下缘为边界，提取椎骨的上下曲面，用造型工具建立各节段椎间盘和上下终板模型。接下来，利用geomagic软件建立实体模型。将在geomagic软件组装完成的模型，导入hypermesh软件中，模型的各部分分别建立为一个元件。对每个部分经过的必要的quickedit处理后，划分网格。网格选用10节点四面体som187单元网格。对划分完成的体网格，进行必要的质量检查，以供后续有限元分析使用。t12-l4节段有限元模型共生成870869个结点，506589个单元。再建立韧带模型。最后在hypermesh软件中对模型的各部分进行材料属性赋值，导入到ansys软件中，模拟人体站立状态，对模型施加边界条件和载荷。
然而，上述文献缺乏精细模拟椎管内占位性骨折块的步骤，仅通过修改椎体高度等方法简单模拟骨折，难以真实反映骨折块与周边骨组织以及软组织的关系，无法进一步研究骨折块复位等手术操作。因此，本专利能更加准确地模拟胸腰椎爆裂性骨折，为骨折机制研究及后期手术研究提供更加精细的仿真模型。
技术实现要素：
本发明的目的是针对现有技术中的不足，为更加有效地研究胸腰椎爆裂性骨折及相关复位方法，通过分析设计胸腰椎爆裂性骨折的几何形状和有限元网格划分，模拟骨折生物力学特性，提供一种更加精确的用于研究胸腰椎爆裂性骨折及相关手术治疗的方法和装置。
第一方面，本发明提供了一种胸腰椎爆裂性骨折仿真模型的建立方法，包括：
步骤s1：筛选健康成人的胸腰段进行ct扫描，获取dicom数据；
步骤s2：重建胸腰椎三维实体模型；
步骤s3：依据正常解剖结构，按照点、线、面的原则，画出胸腰椎外形特征的曲线，通过封装连接得到nurbs曲面模型；
步骤s4：应用有限元建模软件，在第一腰椎上分割骨折块，对骨折块及其他骨组织进行网格划分，进一步划分椎间盘为纤维环和髓核，添加前纵韧带、后纵韧带、棘上韧带、棘间韧带、囊韧带、黄韧带以及横突间韧带，建立骨折块有限元模型和胸腰椎有限元模型，再对骨折块和周边骨组织设定边界条件，组装获得所述胸腰椎爆裂性骨折仿真模型；其中，所述在第一腰椎上分割骨折块的步骤中，骨折块几何边界的确认如下：骨折块在矢状面的长度拟定为椎体中柱的范围，即椎体前后径的1/3，骨折块在冠状位的长度设定为椎弓根基底部内侧缘，骨折块的高度设定自上终板至椎基静脉孔水平；
步骤s5：设定材料属性、载荷、边界条件。
作为本发明的一个优选例，所述材料属性为：皮质骨和松质骨为各向同性材料属性，骨组织材料属性如下：
髓核为不可压缩的流体材料属性，椎间盘材料属性如下：
韧带为非线性材料属性。
作为本发明的另一优选例，所述胸腰椎有限元模型包括皮质、松质，在椎骨皮质骨划分过程中，选择为1mm厚的六面体网格，松质骨为四面体网格；定义小关节为面面接触；终板厚度为0.5mm，髓核区域和纤维环区域分别为椎间盘横截面积的40％和60％。
作为本发明的另一优选例，所述骨折块有限元模型包括外侧的皮质骨和内侧的松质骨。
作为本发明的另一优选例，在确定所述胸腰椎实外形特征的曲线时，只画出几条重要的关键线；对于椎体的上下面，使其由平滑曲面构成。
作为本发明的另一优选例，网格划分步骤先划分一半椎体，通过对称，从而得到整个椎体。
第二方面，本发明提供了如上任一所述建立方法所建立的胸腰椎爆裂性骨折仿真模型在研究胸腰椎爆裂性骨折及手术治疗方法中的用途。
第三方面，本发明提供了一种胸腰椎爆裂性骨折仿真模型的建立装置，包括：
图像获取单元，用于获取胸腰椎骨组织ct扫描图像；
三维重建单元，用于重建胸腰椎的三维数字图像；
胸腰椎实体模型生成单元，用于将三维数字图像转化为实体模型；
椎骨有限元模型生成单元，用于利用胸腰椎实体模型生成胸腰椎椎骨三维有限元模型；
软组织有限元模型生成单元，用于利用胸腰椎实体模型生成胸腰椎椎间盘及韧带三维有限元模型；
骨折块有限元模型生成单元，用于生成胸腰椎骨折块三维有限元模型；
模型组合单元，将椎骨三维有限元模型与椎间盘及韧带三维有限元模型整合，将骨折块三维有限元模型与周边椎骨三维有限元模型整合，设置面面接触精细模拟，获得所述胸腰椎爆裂性骨折仿真模型；
模型参数设定单元，用于对胸腰椎三维有限元模型及骨折块三维有限元模型进行参数设定，所述参数包括材料属性、小关节突的接触、载荷以及边界条件等；
计算分析单元，对所述胸腰椎爆裂性骨折仿真模型进行有限元仿真计算，获得所述胸腰椎爆裂性骨折仿真模型的应力应变云图及椎间活动度参数，提供生物力学数据。
作为本发明的一个优选例，所述椎骨有限元模型生成单元包括：
胸11椎体模型生成子单元，通过有限元网格划分，用于对椎骨三维有限元模型加载力矩和垂直应力；
胸12、腰1、腰2椎体模型生成子单元，用于生成椎骨三维有限元模型活动单元；
腰3椎体模型生成子单元，用于施加限制椎骨三维有限元模型的边界位移。
作为本发明的另一优选例，所述骨折块有限元模型生成单元包括：
骨折轮廓线确定子单元，用于确定腰1椎体骨折的实体模型的几何特征；
骨折块划分子单元，用于对骨折块进行网格划分。
本发明优点在于：
本发明提供的胸腰椎爆裂性骨折仿真模型建立的方法，通过ct获取胸腰椎的图像生成三维数字图像，进一步建立三维实体模型，在此基础上分割骨折块，建立骨折块三维有限元模型和周边胸腰椎三维有限元模型，并进行组合得到胸腰椎爆裂性骨折三维有限元模型，最后对所述模型进行参数设定，包括设定材料属性、小关节突的接触、载荷以及边界条件等。由于本模型包括了椎骨、椎间盘、韧带、小关节等结构，符合胸腰椎的解剖结构，骨折块模型几何特性设置合理，与周围椎骨紧密组合，因此建立的胸腰椎爆裂性骨折三维有限元模型可用于其他类似创伤有限元测试分析，仿真性高，临床应用性更强。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。下面描述中的附图仅仅是本申请的实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为实施例1公开的一种胸腰椎爆裂性骨折仿真模型建立方法的流程图。
图2为实施例1建立的骨折块示意图，深灰色非网格区域表示骨折块。
图3为实施例1建立的胸腰椎三维有限元模型。
图4为屈曲工况下计算所得椎间活动度与文献中标本实验测量所得数据比对。
图5为背伸工况下计算所得椎间活动度与文献中标本实验测量所得数据比对。
图6为侧弯工况下计算所得椎间活动度与文献中标本实验测量所得数据比对。
图7为轴向旋转工况下计算所得椎间活动度与文献中标本实验测量所得数据比对。
图8为实施例2公开的一种胸腰椎爆裂性骨折仿真模型建立方法的流程图。
图9为实施例3公开的一种胸腰椎爆裂性骨折仿真模型建立装置的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
实施例1
请参见图1，图1是实施例1公开的一种胸腰椎爆裂性骨折仿真模型建立方法的流程图。所述方法包括以下步骤：
步骤s1：筛选健康成人志愿者的胸腰段进行ct扫描，获取dicom数据。
步骤s2：将数据导入三维建模软件mimics重建胸腰椎三维实体模型，并保存为stl格式文件。
步骤s3：将stl格式文件导入逆向工程软件geomagic，依据正常解剖结构，按照点、线、面的原则，先设计数条表达主要外形特征的曲线，再封装连接得到nurbs曲面模型，保存为iges文件。
需要说明的是，由于椎骨结构多由不规则的曲面构成，因此在确定胸腰椎外形线时，只需要画出几条重要的关键线；对于椎体的上下面，尽可能由平滑曲面构成，便于后续步骤中终板和椎间盘的划分。通过这种方式，就可以简便快捷地得到胸腰椎的实体模型。
步骤s4：利用hypermesh软件，分割骨折块及周围椎骨组织，建立骨折块有限元模型和胸腰椎有限元模型，进一步组装成胸腰椎骨折三维有限元模型。
具体地，包括：
步骤s401：将上述iges文件导入hypermesh软件，在第一腰椎上分割骨折块(图2)。
需要说明的是，骨折块几何边界的确认根据本申请发明人在本领域长期研究经验获得。基于脊柱三柱理论，中柱受到破坏后导致椎体高度丢失，并导致破裂的中柱骨折块进入椎管，故本方法中骨折块在矢状面的长度拟定为椎体中柱的范围，即椎体前后径的1/3；根据脊椎骨的微结构研究，骨小梁起源于椎弓根基部的内侧角，椎弓根基部为几何结构上变化较大的部位，故骨折块在冠状位的长度设定为椎弓根基底部内侧缘；椎基静脉孔是胸腰椎爆裂性骨折块后缘线的起点，故骨折块的高度设定自上终板至椎基静脉孔水平。
步骤s402：对骨折块及其余椎骨组织进行网格划分，建立骨折块有限元模型和胸腰椎有限元模型(图3)。
胸腰椎有限元模型包括皮质、松质等结构。在椎骨皮质骨划分过程中，选择为1mm厚的六面体网格，松质骨为四面体网格；定义小关节为面面接触；终板厚度为0.5mm，髓核区域和纤维环区域分别为椎间盘横截面积的40％和60％；根据解剖结构添加相关韧带，包括前纵韧带、后纵韧带、黄韧带、棘间韧带、棘上韧带、囊韧带、横突间韧带，进行韧带有限单元模拟。
骨折块有限元模型分为外侧的皮质骨和内侧的松质骨。
需要说明的是，为了简化划分步骤，并且保障对称性，可选择先划分一半椎体，通过对称，从而得到整个椎体。
步骤s403：将骨折块与周边椎骨组合，设定边界条件，设置为面面接触，组装成胸腰椎爆裂性骨折三维有限元模型。
步骤s5：导入abaqus软件，设定材料属性、载荷、边界条件等。
需要说明的是，皮质骨和松质骨为各向同性材料属性，髓核为不可压缩的流体材料属性(表1-2)，韧带为非线性材料属性。
表1骨组织材料属性
表2椎间盘材料属性
步骤s6：计算获得椎间活动度，与文献中的标本实验结果比对，验证模型，进行生物力学分析。
本次脊柱模型的加载处理条件参考以往的相关文献。载荷方面固定l1椎体的下表面，使其在各个方向上的活动度为0，在t11椎体上表面加载5nm的纯力矩，模拟模型前屈、后伸、侧弯和旋转工况下，t11t2、t12l1、l1l2和l2l3的腰椎椎间活动度，计算结果与既往文献biomechanicalanalysisofthethoracolumbarspineunderphysiologicalloadings:experimentalmotiondatacorridorsforvalidationoffiniteelementmodels的生物力学标本试验数据比对，仿真数值也均在文献的实验标准差范围内，验证了模型的合理性和科学性。
实施例2
请参见图6，图6是实施例2公开的一种胸腰椎爆裂性骨折仿真模型建立方法的流程图。所述方法包括以下步骤：
步骤s1：筛选健康成人志愿者的胸腰段进行ct扫描，获取dicom数据。
步骤s2：将数据导入三维建模软件mimics重建胸腰椎三维实体模型，并保存为stl格式文件。
步骤s3：将stl格式文件导入逆向工程软件geomagic，依据正常解剖结构，按照点、线、面的原则，先设计数条表达主要外形特征的曲线，再封装连接得到nurbs曲面模型，保存为iges文件。
需要说明的是，由于椎骨结构多由不规则的曲面构成，因此在确定胸腰椎外形线时，只需要画出几条重要的关键线；对于椎体的上下面，尽可能由平滑曲面构成，便于后续步骤中终板和椎间盘的划分。通过这种方式，就可以简便快捷地得到胸腰椎的实体模型。
步骤s4：利用hypermesh软件，分割骨折块及周围椎骨组织，建立骨折块有限元模型和胸腰椎有限元模型，进一步组装成胸腰椎骨折三维有限元模型。
具体地，包括：
步骤s401：将上述iges文件导入hypermesh软件，在第一腰椎上分割骨折块。
需要说明的是，骨折块几何边界的确认根据本申请发明人在本领域长期研究经验获得。基于脊柱三柱理论，中柱受到破坏后导致椎体高度丢失，并导致破裂的中柱骨折块进入椎管，故本方法中骨折块在矢状面的长度拟定为椎体中柱的范围，即椎体前后径的1/3；根据脊椎骨的微结构研究，骨小梁起源于椎弓根基部的内侧角，椎弓根基部为几何结构上变化较大的部位，故骨折块在冠状位的长度设定为椎弓根基底部内侧缘；有文献显示椎基静脉孔是胸腰椎爆裂性骨折块后缘线的起点，故骨折块的高度设定自上终板至椎基静脉孔水平。
步骤s402：对骨折块及其余椎骨组织进行网格划分，建立骨折块有限元模型和胸腰椎有限元模型。
胸腰椎有限元模型包括皮质、松质等结构。在椎骨皮质骨划分过程中，选择为1mm厚的六面体网格，松质骨为四面体网格；终板厚度为0.5mm，髓核区域和纤维环区域分别为椎间盘横截面积的40％和60％；根据解剖结构添加相关韧带，包括前纵韧带、后纵韧带、黄韧带、棘间韧带、棘上韧带、囊韧带、横突间韧带，进行韧带有限单元模拟。
骨折块有限元模型分为外侧的皮质骨和内侧的松质骨。
需要说明的是，为了简化划分步骤，并且保障对称性，可选择先划分一半椎体，通过对称，从而得到整个椎体。
步骤s403：将骨折块与周边椎骨组合，设置为面面接触，组装成胸腰椎骨折三维有限元模型。
步骤s5：导入abaqus软件，设定材料属性、载荷、边界条件等。
需要说明的是，皮质骨和松质骨为各向同性材料属性，髓核为不可压缩的流体材料属性(表1-2)，韧带为非线性材料属性。
表1骨组织材料属性
表2椎间盘材料属性
实施例3
下面对本实施例提供的胸腰椎爆裂性骨折模型建立装置进行描述，其与上文描述的胸腰椎爆裂性骨折模型建立方法可相互对应参照。
所述胸腰椎爆裂性骨折模型建立装置，包括：
图像获取单元1，用于获取胸腰椎骨组织ct扫描图像；
三维重建单元2，用于重建胸腰椎的三维数字图像；
胸腰椎实体模型生成单元3，用于将三维数字图像转化为实体模型；
椎骨有限元模型生成单元4，用于利用胸腰椎实体模型生成胸腰椎椎骨三维有限元模型；
软组织有限元模型生成单元5，用于利用胸腰椎实体模型生成胸腰椎椎间盘及韧带三维有限元模型；
骨折块有限元模型生成单元6，用于生成胸腰椎骨折块三维有限元模型；
模型组合单元7，将椎骨三维有限元模型与椎间盘等软组织有限元模型整合，将骨折块三维有限元模型与周边椎骨三维有限元模型整合，设置面面接触精细模拟；
模型参数设定单元8，用于对胸腰椎三维有限元模型及骨折块三维有限元模型进行参数设定，所述参数包材料属性、小关节突的接触、载荷以及边界条件等；
计算分析单元9，对所述胸腰椎爆裂性骨折模型进行有限元仿真计算，获得所述胸腰椎爆裂性骨折模型的应力应变云图及椎间活动度等参数，为临床工作提供生物力学数据。
所述椎骨有限元模型生成单元4可进一步包括：
胸11椎体模型生成子单元，通过有限元网格划分，用于对椎骨三维有限元模型加载力矩和垂直应力；
胸12、腰1、腰2椎体模型生成子单元，用于生成椎骨三维有限元模型活动单元；
腰3椎体模型生成子单元，用于施加限制椎骨三维有限元模型的边界位移；
所述骨折块有限元模型生成单元6可进一步包括：
骨折轮廓线确定子单元，用于确定腰1椎体骨折的三维实体模型的几何特征；
骨折块划分子单元，用于对骨折块进行网格划分。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员，在不脱离本发明方法的前提下，还可以做出若干改进和补充，这些改进和补充也应视为本发明的保护范围。
1.一种胸腰椎爆裂性骨折仿真模型的建立方法，其特征在于，包括：
步骤s1：筛选健康成人的胸腰段进行ct扫描，获取dicom数据；
步骤s2：重建胸腰椎三维实体模型；
步骤s3：依据正常解剖结构，按照点、线、面的原则，画出胸腰椎实外形特征的曲线，通过封装连接得到nurbs曲面模型；
步骤s4：应用有限元建模软件，在第一腰椎上分割骨折块，对骨折块及其他骨组织进行网格划分，进一步划分椎间盘为纤维环和髓核，添加前纵韧带、后纵韧带、棘上韧带、棘间韧带、囊韧带、黄韧带以及横突间韧带，建立骨折块有限元模型和胸腰椎有限元模型，再对骨折块和周边骨组织设定边界条件，组装获得所述胸腰椎爆裂性骨折仿真模型；其中，所述在第一腰椎上分割骨折块的步骤中，骨折块几何边界的确认如下：骨折块在矢状面的长度拟定为椎体中柱的范围，即椎体前后径的1/3，骨折块在冠状位的长度设定为椎弓根基底部内侧缘，骨折块的高度设定自上终板至椎基静脉孔水平；
步骤s5：设定材料属性、载荷、边界条件。
2.根据权利要求1所述的建立方法，其特征在于，所述材料属性为：皮质骨和松质骨为各向同性材料属性，骨组织材料属性如下：
髓核为不可压缩的流体材料属性，椎间盘材料属性如下：
韧带为非线性材料属性。
3.根据权利要求1所述的建立方法，其特征在于，所述胸腰椎有限元模型包括皮质、松质，在椎骨皮质骨划分过程中，选择为1mm厚的六面体网格，松质骨为四面体网格；定义小关节为面面接触；终板厚度为0.5mm，髓核区域和纤维环区域分别为椎间盘横截面积的40％和60％。
4.根据权利要求1所述的建立方法，其特征在于，所述骨折块有限元模型包括外侧的皮质骨和内侧的松质骨。
5.根据权利要求1所述的建立方法，其特征在于，在确定所述胸腰椎外形特征的曲线时，只画出几条重要的关键线；对于椎体的上下面，使其由平滑曲面构成。
6.根据权利要求1所述的建立方法，其特征在于，网格划分步骤先划分一半椎体，通过对称，从而得到整个椎体。
7.如权利要求1-6任一所述建立方法所建立的胸腰椎爆裂性骨折仿真模型在研究胸腰椎爆裂性骨折及手术治疗方法中的用途。
8.一种胸腰椎爆裂性骨折仿真模型的建立装置，其特征在于，包括：
图像获取单元，用于获取胸腰椎骨组织ct扫描图像；
三维重建单元，用于重建胸腰椎的三维数字图像；
胸腰椎实体模型生成单元，用于将三维数字图像转化为实体模型；
椎骨有限元模型生成单元，用于利用胸腰椎实体模型生成胸腰椎椎骨三维有限元模型；
软组织有限元模型生成单元，用于利用胸腰椎实体模型生成胸腰椎椎间盘及韧带三维有限元模型；
骨折块有限元模型生成单元，用于生成胸腰椎骨折块三维有限元模型；
模型组合单元，将椎骨三维有限元模型与椎间盘及韧带三维有限元模型整合，将骨折块三维有限元模型与周边椎骨三维有限元模型整合，设置面面接触精细模拟，获得所述胸腰椎爆裂性骨折仿真模型；
模型参数设定单元，用于对胸腰椎三维有限元模型及骨折块三维有限元模型进行参数设定，所述参数包括材料属性、小关节突的接触、载荷以及边界条件等；
计算分析单元，对所述胸腰椎爆裂性骨折仿真模型进行有限元仿真计算，获得所述胸腰椎爆裂性骨折仿真模型的应力应变云图及椎间活动度参数，提供生物力学数据。
9.根据权利要求8所述的建立装置，其特征在于，所述椎骨有限元模型生成单元包括：
胸11椎体模型生成子单元，通过有限元网格划分，用于对椎骨三维有限元模型加载力矩和垂直应力；
胸12、腰1、腰2椎体模型生成子单元，用于生成椎骨三维有限元模型活动单元；
腰3椎体模型生成子单元，用于施加限制椎骨三维有限元模型的边界位移。
10.根据权利要求8所述的建立装置，其特征在于，所述骨折块有限元模型生成单元包括：
骨折轮廓线确定子单元，用于确定腰1椎体骨折的实体模型的几何特征；
骨折块划分子单元，用于对骨折块进行网格划分。
技术总结