模拟，通过 改变断端微动距离和血供参数，从而得到不同状态下的骨折愈合情况，确定出最佳的骨折 愈合方案以及对骨折愈合产生有益影响的断端微动距离和血供状态。
【主权项】
1. 一种基于断端微动和血供的骨折愈合仿真系统，其特征在于，所述系统包括： 骨折愈合仿真子系统（1)、输入模块(2)和显示模块(3); 其中所述骨折愈合仿真子系统（1)包括： 几何建模模块：用于根据导入的二维断层扫描图像数据，经过图像预处理后进行骨折 处骨挪的Ξ维表面模型的建立； 网格划分模块：用于Ξ维表面模型的面网格划分，并用于将面网格模型生成体网格模 型，实现连续的Ξ维几何模型的离散化，保证能够进行后续的有限元计算，从而得到骨挪Ξ 维有限元模型； 环境初始化模块：用于对建立好的骨挪Ξ维有限元模型进行初始化设置，包括生物力 学环境初始化设置和生物学环境初始化设置； 生物力学计算模块:用于计算在骨折愈合过程中断端微动的距离，并通过有限元分析， 计算骨挪内单元的应变刺激;另外还包括不同愈合阶段，骨挪单元材料属性的计算； 生物学环境模块:用于在骨折愈合过程中，构建骨挪处的生化环境，完成骨挪内的组织 分化过程； 其中所述的输入模块(2)包括： 图像数据输入单元，用于读取患者的二维断层扫描图像数据； 用户操作单元，用于用户与骨折愈合仿真系统软件之间的信息传递与处理，信息包括 患者的个人基本信息和骨挪模型视图； 其中，所述的显示模块(3)，用于观察骨折愈合情况，实现用户与系统之间的视觉交互。2. 根据权利要求书1所述的一种基于断端微动和血供的骨折愈合仿真系统，其特征在 于，所述的骨折愈合仿真子系统（1)实现其功能的具体过程为： 1) 几何建模模块 通过阔值筛选、交互式分割W及Ξ维重建来实现骨挪的Ξ维几何建模;通过阔值筛选 的过程将患者的骨、肌肉、软骨、脂肪、皮肤组织分别进行初步的分割;通过交互式分割过程 手动擦除和修补单张图片来提高图像分割的质量和速度;通过Ξ维重建将在二维断层图像 中分割出来的区域生成Ξ维表面模型； 2) 网格划分模块 网格划分包括面网格划分和体网格划分两个步骤;面网格划分过程用于将Ξ维表面模 型进行优化，包括:表面模型优化、平滑处理、修补漏桐;表面模型的优化通过减小表面模型 的Ξ角面片来实现，该过程只需将相连的两个顶点合并到一个新顶点上，并延续原有的拓 扑关系；平滑处理的过程中，对Ξ维的面网格模型进行去噪;修补漏桐的过程中，通过将模 型当中的空桐提取成空间多边形，然后对空桐多边形进行Ξ角化的方法实现;体网格划分 的过程是将面网格模型进行拉伸、旋转步骤来实现； 3) 环境初始化模块 所述骨挪由Ξ种组织类型组成：结缔组织、软骨组织和骨组织；在骨折愈合初期，骨挪 主要由结缔组织组成，所W骨挪材料属性为结缔组织的材料属性； 在骨折愈合初期，骨挪处的血液供应有Ξ部分组成：皮质骨内的血液供应为100% ;由 于骨折断端处组织受到严重损伤，所W骨折断端处的血液供应为0;除了骨折处自身的血液 供应外，骨折周围的软组织也会向骨折处进行血液供应，将其设定为30%; 环境初始化模块还用于外加载荷的设置及皮质骨材料属性的赋值； 4)生物力学计算模块： 主要包括骨折断端微动距离的计算、骨挪单元应变的计算和骨挪单元材料属性的计 算； a. 骨折断端微动的计算： 骨折断端微动的计算包括当前骨挪硬度的计算和断端微动的计算； 假设骨挪具有线性，各向同性的性质;所W有如下关系： Fn 二 kcalWcal , 式中，Fn为骨挪受到的载荷，kca功当前骨挪硬度，μ。。功骨挪移动距离； 其中Fn可由下式求得： Fn = 〇A 式中，Fn为骨挪受到的载荷，σ为骨挪应力，A为受力面积； 骨挪应力计算公式为： σ = Εε 式中，σ为骨挪应力，Ε为骨挪弹性模量，ε为骨挪应变； 其中€ = 站为骨挪相对移动距离，L为骨折间距； 通过上述公式的计算，即可求得骨挪硬度kcal; 外加载荷、外部固定器和骨挪之间具有如下的关系： F(y)=Ffix(y)+kcaiy 式中，F为外加载荷，Ffix为固定器所承受的载荷，kcal为当前骨挪硬度，μ为断端微动距 离； 断端微动距离μ可由ρ(μ)的反函数μ(ρ)求得； b. 骨挪单元应力的计算 用体积应变和崎变应变来表示骨挪的单元应力，如下式所示：式中，ε〇为体积应变，γ 0为崎变应变； 利用单元主应变ε?，62，63分别求得体积应变ε〇和崎变应变丫 0，关系式如下所示：其中，ε〇为体积应变代表体积的改变；γ 0为崎变应变代表形状的改变； C.骨挪单元材料属性的计算 随着骨折愈合的进行，骨挪单元内的结缔组织、软骨组织和骨组织的浓度发生改变，从 而使骨挪单元的材料属性发生变化，计算公式如下：式中，Eele为单元弹性模量;E。。。为结缔组织弹性模量，C。。。为结缔组织浓度;Esart为软骨 组织弹性模量，Csart为软骨组织浓度;Ebone为骨组织弹性模量，CbDne为骨组织浓度； 白 ele -白 conCcon+白 cartCcart+白 boneCbone 式中，0ele为单元泊松比；Θ。。。为结缔组织泊松比，C。。。为结缔组织浓度；0Eart为软骨组织 泊松比，Csart为软骨组织浓度;0bDne为骨组织泊松比，CbDne为骨组织浓度； 由于骨挪是由结缔组织、软骨组织和骨组织的混合物组成，有下式： Ccon+Ccart+Cbone- 1 式中，C。。。为结缔组织浓度，Csart为软骨组织浓度，CbDne为骨组织浓度； 联立上述公式即可求得骨挪单元材料属性； 5)生物学环境模块 通过模糊控制来实现骨挪内的组织分化过程:其中模糊控制是一个关于屯个输入变量 和Ξ个输出变量的一种控制过程;屯个输入变量分别为:单元内当前血供、软骨浓度和骨浓 度，临近单元的血供和骨浓度W及当前单元的体积应变和崎变应变;Ξ个输出变量分别为： 当前单元血供变化、软骨浓度变化和骨浓度变化;所述的控制过程具体为:依据临床组织学 得到不同阶段骨挪内的组织分布情况获得输入变量的隶属度函数和输出变量的隶属度函 数，根据输入变量的隶属度函数可将输入变量模糊化;通过模糊规则W及输出变量的隶属 度函数，再通过反模糊化即可得到更新后的血供、软骨组织和骨组织状态； 由此可实现组织分化过程。3. 根据权利要求书1所述的一种基于断端微动和血供的骨折愈合仿真系统，其特征在 于:所述输入模块(2)实现其功能的具体过程为： 1) 图像数据输入单元 将二维断层扫描图像通过计算机输入设备导入到骨折愈合仿真系统软件中，为几何建 模做准备； 2) 用户操作单元 用于操作单元通过计算机输入设备输入手术相关信息，并通过操作输入设备与骨折愈 合仿真系统软件的操作界面进行交互。4. 根据权利要求书1所述的一种基于断端微动和血供的骨折愈合仿真系统，其特征在 于:所述显示模块(3)实现其功能的具体过程为： 将计算机的显示设备连接到计算机，与用于进行视觉上的交互。5. 根据权利要求书1所述的一种基于断端微动和血供的骨折愈合仿真系统，其特正在 于，所述骨折愈合仿真软件子系统（1)中的生物力学计算模块和生物学环境模块功能都是 由迭代循环实现的。
【专利摘要】一种基于断端微动和血供的骨折愈合仿真系统，涉及生物医学工程领域。本发明用来预测骨折愈合的复杂过程，探寻骨折愈合中尚无法确定的机理。所述系统包括骨折愈合仿真子系统、输入模块和显示模块；其中所述骨折愈合仿真子系统(1)包括：几何建模模块、网格划分模块、环境初始化模块、生物力学计算模块和生物学环境模块。生物学环境模块用于在骨折愈合过程中，构建骨痂处的生化环境，完成骨痂内的组织分化过程。本发明将血液供给作为动态变量引入到骨折愈合仿真系统中，能够更加准确的模拟骨折愈合过程。本发明可以无限多次重复进行实验模拟，从而对找到最佳的断端微动值和血液对骨折愈合的影响提供有益帮助。
【IPC分类】G06F17/50, G06T17/30
【公开号】CN105550461
【申请号】CN201511027837
【发明人】王沫楠, 杨宁
【申请人】哈尔滨理工大学
【公开日】2016年5月4日
【申请日】2015年12月30日