[发明专利]一种基于组织氧气环境以及力学环境的骨折愈合仿真方法

摘要

本发明公开了一种基于组织氧气环境以及力学环境的骨折愈合仿真方法，具体包括以下步骤：建立骨折部位的三维几何模型；将得到的几何模型进行网格划分，建立骨和骨痂的有限元模型；对有限元模型进行分析解算，求解出骨折区域的偏应变；根据力刺激，建立骨折区域的血管生长模型；建立骨折区域各细胞增殖分化模型；建立骨折区域的氧气运输模型；计算组织内新的材料属性；建立力生物调节算法的骨折愈合迭代仿真过程；根据计算得到的组织材料属性值判断组织内的材料属性是否达到骨的材料属性值，若达到骨的材料属性值则仿真完成，否则更新组织材料属性值并进入下一迭代仿真。

主权项

1.一种基于组织氧气环境以及力学环境的骨折愈合仿真方法，包括如下步骤：步骤一：建立骨折部位的三维几何模型；步骤二：将得到的几何模型进行网格划分，建立骨和骨痂的有限元模型；步骤三：对骨痂内初始时刻的材料属性进行赋值；步骤四：对有限元模型进行分析解算，求解出骨折区域的偏应变；步骤五：根据计算得到的力刺激，建立骨折区域的血管生长模型；步骤六：建立骨折区域细胞增殖分化模型，包括间充质干细胞增殖分化模型、软骨细胞增殖分化模型、成纤维细胞增殖分化模型以及骨细胞增殖分化模型；步骤七：建立骨折区域的氧气运输模型；步骤八：计算组织内新的材料属性；步骤九：根据以上步骤建立力生物调节算法的骨折愈合迭代仿真过程，根据步骤八计算得到的组织材料属性值判断组织内的材料属性是否达到骨的材料属性值，若达到骨的材料属性值则仿真完成，否则更新组织材料属性值并进入下一迭代仿真；其特征在于：所述步骤五中血管生长模型如下：式中，Cvessel为血管浓度，Dvessel为骨痂内血管扩散系数,代表新血管在骨痂内生长的速度，t为时间，x，y，z分别为三维笛卡尔坐标系的x轴，y轴，z轴，血管生长的扩散源为髓腔和皮质骨膜处。