[发明专利]一种模拟骨折愈合过程的仿真系统

权利要求书

1.一种模拟骨折愈合过程的仿真系统，其特征在于，所述系统包括：

骨折区域分析模型建立模块(1)、血运重建子系统(2)、骨折愈合过程建模子系统(3)和程序终止判断模块(4)；

其中，所述的骨折区域分析模型建立模块(1)用于建立骨折区域几何模型和有限元分析模型；

所述的骨折区域分析模型建立模块(1)，实现其功能的具体过程为：

1)骨折区域三维表面几何模型的建立

采用基于分割的三维医学影像表面重建算法对图像进行表面重构，通过阈值筛选、交互式分割和三维重建过程得到三维表面几何模型；

所述的影像由影像设备CT得到，数据存储格式为DICOM；

2)骨折区域有限元分析模型的建立

将骨折区域三维表面几何模型进行网格划分，使连续的几何模型离散化，得到骨折区域有限元分析模型；

所述的网格划分包括面网格划分和体网格划分两个步骤：面网格划分过程用于将三维表面几何模型进行优化，包括：表面模型优化、平滑处理、修补漏洞；表面模型的优化通过减少表面模型的三角面片来实现，该过程只需将相邻的两个顶点合并到一个新的顶点上，并延续原有的拓扑关系；平滑处理的过程中，对三维的面网格模型进行去噪；修补漏洞的过程中，通过将模型当中的空洞提取成空间多边形，然后对空间多边形进行三角化的方法实现；体网格划分的过程是将面网格模型进行拉伸、旋转步骤实现的；

通过网格划分得到的骨折区域有限元分析模型包括单元编号和节点坐标两部分；

单元编号包含四列数据，四列数据分别为每个单元的四个节点的节点序号；

节点坐标包含三列数据，三列数据分别为每个节点的空间坐标值；

其中，所述的血运重建子系统(2)包括：

细胞内部分子生理活动建模模块，用于从细胞内分子层次中相关分子生理活动对骨折愈合过程中的血运重建过程进行建模；

所述的相关分子包括血管生成细胞生长因子受体，Notch1蛋白，Dll4蛋白，激活的血管生成细胞生长因子受体，激活的Notch1蛋白，有效激活的血管生成细胞生长因子受体，有效激活的Notch1蛋白和肌动蛋白；

细胞生理活动建模模块，用于从细胞层次中血管生成细胞生理活动对骨折愈合过程中的血运重建过程进行建模；

其中，所述的骨折愈合过程建模子系统(3)包括：

相关细胞浓度建模模块，用于建立骨折愈合过程中相关细胞浓度随时间的变化关系；

所述的相关细胞包括：间充质干细胞、成纤维细胞、软骨细胞和骨细胞；

相关生长因子浓度建模模块，用于建立骨折愈合过程中相关生长因子浓度随时间的变化关系；

所述的相关生长因子包括：软骨细胞生长因子、骨细胞生长因子和血管生成生长因子；

相关细胞胞外基质浓度建模模块，用于建立骨折愈合过程中相关细胞胞外基质浓度随时间的变化关系；

所述的相关细胞胞外基质包括：成纤维细胞胞外基质、软骨细胞胞外基质和骨细胞胞外基质；

采用耦合偏微分方程组的方法对上述相关细胞浓度、相关生长因子浓度和相关细胞胞外基质浓度随时间的变化关系进行建模；

其中，所述的程序终止判断模块(4)用于判断程序是否终止，若不满足终止条件，则程序进行下一循环步；若满足终止条件，则程序结束，并记录骨折愈合时间，绘制相关细胞浓度、相关生长因子浓度和相关细胞胞外基质浓度随时间的变化云图；

所述的程序终止判断模块(4)，实现其功能的具体过程为：

判断当前分析步中，骨折区域单元骨浓度是否达到最大值并保持不变，若不满足上述条件，则程序进入下一循环步；若满足上述条件，则程序终止，并记录骨折愈合时间，绘制间充质干细胞浓度、成纤维细胞浓度、软骨细胞浓度、骨细胞浓度、软骨细胞生长因子浓度、骨细胞生长因子浓度、血管生成细胞生长因子浓度、成纤维细胞胞外基质浓度、软骨细胞胞外基质浓度和骨细胞胞外基质浓度随时间的变化云图。