CN106777584B - 一种模拟骨折愈合过程的仿真系统 - Google Patents
一种模拟骨折愈合过程的仿真系统 Download PDFInfo
- Publication number
- CN106777584B CN106777584B CN201611086920.5A CN201611086920A CN106777584B CN 106777584 B CN106777584 B CN 106777584B CN 201611086920 A CN201611086920 A CN 201611086920A CN 106777584 B CN106777584 B CN 106777584B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- cell
- bone
- tissue
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
Abstract
一种模拟骨折愈合过程的仿真系统,涉及生物医学工程领域。本发明用来预测骨折愈合的复杂过程,为医生制定良好的骨折愈合手术方案提供指导意见。所述系统包括三维几何建模模块、有限元建立模块、骨折部位生物力学建模模块、骨折部位细胞进化建模模块、骨痂生长建模模块和显示模块。通过骨折部位生物力学建模模块和骨折部位细胞进化模块对细胞浓度及组织材料属性进行更新;通过骨痂生长建模模块对骨痂几何形状进行更新;通过显示模块判断骨折愈合的效果。本发明可以更加准确的模拟骨折愈合过程,可以无限多次重复进行实验模拟,从而对制定最佳的骨折愈合手术方案提供有益帮助。
Description
技术领域
本发明涉及生物医学工程领域,尤其涉及一种模拟骨折愈合过程的仿真系统。
背景技术
肌骨系统在人体中起到重要的作用,其主要功能包括:运动功能、保护功能、支持功能、造血功能、贮存功能。骨折的高发性以及骨折后对患者造成的痛苦和对社会造成的经济负担使得骨折愈合的研究显得尤为迫切与重要。尽管关于骨折愈合的研究一直备受关注,但仍有5%~10%的骨折因各种原因发生延迟愈合甚至是不愈合。
骨折愈合是一个复杂的过程,包括骨髓间充质干细胞的迁移,增殖,分化和凋亡。骨髓间充质干细胞向成纤维细胞、软骨细胞、骨细胞分化,这些细胞最终形成纤维结缔组织,软骨组织和骨组织。这些细胞活动由骨折区域所受的力学刺激决定。骨折愈合的过程中,骨痂的形状以及尺寸也会随着骨折区域的细胞活动发生改变,从而影响骨折区域所受到的力学刺激。尽管关于骨折愈合的研究一直在进行,但仍存在许多缺点与不足:
1.没有建立针对患者的个体化模型;
2.关于骨折愈合的研究很多只停留在组织层面,没有对骨折区域中的细胞活动进行研究;
3.力学因素与细胞分化没有一个确定性关系;
4.没有考虑骨折愈合过程中骨痂形状和尺寸的改变;
5.骨折部位的模型和生物学材料设置过于简化。
发明内容
本发明的目的是提供一种模拟骨折愈合过程的仿真系统,用来预测骨折愈合的复杂过程,从而寻找最佳的骨折愈合的手术规划方案。
本发明的目的通过下述技术方案实现:一种模拟骨折愈合过程的仿真系统,其特征为骨折部位三维几何建模模块、骨折部位有限元建模模块、骨折部位生物力学建模模块、骨折部位细胞进化建模模块、骨痂生长建模模块和显示模块;
骨折部位三维几何建模模块,用于根据导入的二维断层扫描图像数据,经过图像预处理后进行骨折部位的三维表面几何模型的建立;
骨折部位有限元建模模块,用于对得到的骨折部位三维表面几何模型进行网格划分,实现连续的三维几何模型的离散化,得到节点坐标和单元坐标;还用于将细胞浓度和组织体积分数存储到单元节点中,与节点坐标共同组成节点信息;单元编号和节点信息共同构成本发明所需要的骨折部位有限元模型;
骨折部位生物力学建模模块,用于建立骨折部位生物力学模型,求解骨痂单元力学刺激;
骨折部位细胞进化建模模块,用于对骨折部位细胞的迁移,增殖,分化,凋亡进行建模,分析骨折部位的细胞进化。
骨痂生长建模模块,用于对骨痂形状和尺寸随时间的变化关系进行建模,分析骨痂随时间的变化情况。
显示模块,用于观察骨折的愈合情况,实现用户与系统之间的视觉交互。
本发明的有益效果是:
1.本发明所提出的模拟骨折愈合过程的仿真系统是基于windows开发语言平台来开发软件,通过自主编程实现骨折愈合过程的模拟,基于对话框的形式,易于操作,培训周期短;
2.将骨痂在骨折愈合过程中尺寸和形状的变化加入模拟骨折愈合过程的仿真系统中,能够更加真实的反映骨折愈合过程,是仿真结果更加精确;
3.将骨折区域设置为双相多孔弹性模型,相比单相模型,更加符合骨折区域的生物特性,使仿真结果更加精确;
4.通过仿真细胞的扩散、增殖、分化和凋亡来模拟骨折区域的细胞活动,更加符合骨折区域组织分化的实质,是仿真结果更加精确;
5.通过构建骨折愈合仿真方法,可以对医生制定最优的手术方案提供指导,进而提高手术成功率,提高骨折愈合质量,减少骨折骨折不愈合和延迟愈合的情况。
6.通过构建骨折愈合仿真方法,可以对建立的仿真模型进行多次重复实验,减少真实的生物实验,节省时间,提高效率,节省费用,避免人道主义争议。
综上,本发明的仿真方法克服了现有技术的缺点与不足。
附图说明
图1为模拟骨折愈合仿真系统中各模块之间的关系示意图;
图2为骨折部位三维几何模型建立流程图;
图3为有限元模型建立流程图;
图4为骨痂单元力学刺激求解流程图;
图5为骨折部位细胞进化建模模块、骨痂生长模块与骨痂单元力学刺激之间的关系示意图。
具体实施方式
具体实施方式一:如图1所示,本实施方式所述的一种模拟骨折愈合过程的仿真系统包括:骨折部位三维几何建模模块、骨折部位有限元建模模块、骨折部位生物力学建模模块、骨折部位细胞进化建模模块、骨痂生长建模模块和显示模块;
所述的骨折部位三维几何建模模块,用于根据导入的二维断层扫描图像数据,经过图像预处理后进行骨折部位的三维表面几何模型的建立;
所述的骨折部位有限元建模模块,用于对得到的骨折部位三维表面几何模型进行网格划分,实现连续的三维几何模型的离散化,得到节点坐标和单元编号;还用于将细胞浓度和组织体积分数存储到单元节点中,与节点坐标共同组成节点信息;单元编号和节点信息共同构成所需要的骨折部位有限元模型;
所述的骨折部位生物力学建模模块,用于建立骨折部位生物力学模型,求解骨痂单元力学刺激;
所述的骨折部位细胞进化建模模块,用于对骨折部位细胞的迁移,增殖,分化,凋亡进行建模,分析骨折部位的细胞进化。
所述的骨痂生长建模模块,用于对骨痂形状和尺寸随时间的变化关系进行建模,分析骨痂随时间的变化情况。
所述的显示模块,用于观察骨折的愈合情况,实现用户与系统之间的视觉交互。
具体实施方式二:如图2所示,本实施方式中,所述的骨折部位三维几何建模模块,实现其功能的具体过程为:
采用基于分割的三维医学影像表面重建算法对图像进行三维表面重构,通过阈值筛选、交互式分割和三维重建过程得到三维表面几何模型;
所述的影像由影像设备CT得到,数据存储格式为DICOM。
本实施方式的其他组成与连接关系与具体实施方式一相同。
具体实施方式三:如图3所示,本实施方式中,所述的骨折部位有限元建模模块,实现其功能的具体过程为:
将建立的三维表面几何模型进行网格划分,使连续的三维几何模型离散化,得到节点坐标和单元编号;将细胞浓度和组织体积分数存储到单元节点中,与节点坐标共同构成节点信息;由节点信息和单元编号组成所需要的骨折部位三维有限元模型;
所述的网格划分包括面网格划分和体网格划分两个步骤;面网格划分过程用于将三维表面模型进行优化,包括:表面模型优化,平滑处理,修补漏洞;表面模型的优化通过减小表面模型的三角面片来实现,该过程只需将相邻的两个顶点合并到一个新的顶点上,并延续原有的拓扑关系;平滑处理的过程中,对三维的面网格模型进行去噪;修补漏洞的过程中,通过将模型当中的空洞提取成空间多边形,然后对空洞多边形进行三角化的方法实现;体网格划分的过程是将面网格模型进行拉伸、旋转步骤来实现的。
本实施方式的其他组成与连接关系与具体实施方式一或二之一相同。
具体实施方式四:如图4所示,本实施方式中,所述的骨折部位生物力学建模模块,实现其功能的具体过程为:
将骨折部位看成是多孔弹性材料,则有如下关系式:
a.固体基质、液相和总的应力应变关系如下:
σs=-φspI+σE (1)
σf=-φfpI (2)
σt=σs+σf=-pI+σE (3)
式中,σs、σf、σt分别为固相、液相和总的应力张量;p为液体压力;φs、φf分别为固相和液相体积分数;σE为有效应力张量;I为单位张量;
线弹性材料的有效应力张量可以表示为:
σE=Cε (4)
式中,σE为有效应力张量;C为刚度张量;ε为总的弹性应变;
刚度张量由下式表示:
式中,E为弹性模量;υ为泊松比;
b.考虑到两相均不可压缩性和各向同性,多孔弹性模型的连续性方程为:
式中,φf为液相的体积分数;vs、vf分别为固相和液相的速度向量;
c.固相和液相的动量方程如下:
式中,πs、πf分别为固相和液相的体力向量;φf为液相体积分数;k为渗透率;vf、vs分别为固相和液相的速度向量;σs、σf、σE分别为固相、液相和有效的应力张量;p为液体压力;
d.偏应变张量的第二不变量可由下式表示:
式中,Ψ(x,t)为骨痂单元力学刺激;J2为偏应变张量第二不变量;ε1、ε2、ε3为主应变,εoct为八面体应变;
其中,八面体应变εoct可由下式求得:
通过有限元法求解上述方程,得到偏应变张量的第二不变量J2,并以此作为骨痂单元的力学刺激Ψ(x,t)。
本实施方式的其他组成与连接关系与具体实施方式一至三之一相同。
具体实施方式五:如图5所示,本实施方式中,所述的骨折部位细胞进化建模模块,实现其功能的具体过程为:
细胞的进化包括细胞的迁移、增殖、分化和凋亡。骨折愈合过程中所涉及到的细胞包括:骨髓间充质干细胞、成纤维细胞、软骨细胞和骨细胞。涉及到的组织包括肉芽组织、纤维结缔组织、软骨组织和骨组织。
细胞进化建模模块分为骨髓间充质干细胞进化建模子模块、骨细胞形成建模子模块和材料属性求解子模块。
1)骨髓间充质干细胞进化建模子模块
在骨折愈合初期,骨折部位组织被破坏,骨髓间充质干细胞侵入骨折部位并进行增殖和分化。根据骨痂单元力学刺激的不同,骨髓间充质干细胞可分化为成纤维细胞、软骨细胞和骨细胞。当骨痂单元力学刺激过高时,骨髓间充质干细胞发生凋亡。当骨髓间充质干细胞的浓度达到饱和状态,细胞浓度增加。
因此,骨髓间充质干细胞的进化模型可由下式表示:
式中,c干细胞为骨髓间充质干细胞浓度,单位为细胞数目/mm3;V破坏为被破坏组织的体积分数;Ψ为骨痂单元力学刺激;α增殖、Ψ增殖为骨髓间充质干细胞分化常数;D(V破坏)为扩散系数,单位为mm2/天;
其中扩散系数D(V破坏)下式表示:
式中,D0为常数;V破坏为被破坏组织的体积分数;
骨髓间充质干细胞的分化可由下式表示:
式中,o分化(Ψ,t)、c分化(Ψ,t)、l分化(Ψ,t)为骨髓间充质干细胞向骨细胞、软骨细胞、成纤维细胞分化;-c干细胞为骨髓间充质干细胞发生凋亡;为每种细胞类型的成熟时间;Ψ极限、Ψ骨、Ψ软骨、Ψ纤维、Ψ凋亡为每种细胞类型的力学刺激极限;
骨细胞成熟时间可由下式表示:
当Ψ极限<Ψ<Ψ骨 (16)
式中,t成熟为成熟时间常数;
软骨细胞成熟时间可由下式表示:
当Ψ骨<Ψ<Ψ软骨 (17)
成纤维细胞成熟时间可由下式表示:
当Ψ软骨<Ψ<Ψ纤维 (18)
式中,t极限为时间极限。
2)骨细胞形成建模子模块
骨细胞的形成过程包括膜内骨化和软骨骨化两种途径。膜内骨化是指间充质干细胞直接分化为骨细胞。软骨骨化需经过软骨生成之后才能进行软骨骨化过程;
所以骨细胞的形成可由下式表示:
式中,c骨为骨细胞浓度;o分化(Ψ,t)为由膜内骨化形成的骨细胞;e分化(Ψ,t)为由软骨骨化形成的骨细胞;
其中,由膜内骨化形成的骨细胞o分化(Ψ,t)可由下式表示:
式中,D骨为扩散常数,单位为mm2/天;c骨为骨浓度,单位为细胞数目/mm3;Ψ极限、Ψ骨为膜内骨化过程中骨痂单元力学刺激极限值;为膜内骨化过程中最小骨浓度;为骨髓间充质干细胞分化为骨细胞所用成熟时间;
由软骨骨化形成骨细胞e分化(Ψ,t)可由下式表示:
式中,c骨为骨细胞浓度,单位为细胞数目/mm3;c软骨为软骨细胞浓度,单位为细胞数目/mm3;p为被矿化的矿化的软骨的百分比;p最小、为常数;
软骨细胞浓度与骨髓间充质干细胞的分化和软骨细胞的骨化有关,可由下式表示:
式中,g分化(Ψ,t)为由骨髓间充质干细胞分化得到的软骨细胞;g骨化(Ψ,t)为软骨骨化过程;
其中,由骨髓间充质干细胞分化得到的软骨细胞g分化(Ψ,t)可由下式表示:
软骨骨化过程g骨化(Ψ,t)可由下式表示:
g骨化=-c分化 (24)
骨髓间充质干细胞除了可以分化为骨细胞和软骨细胞,还可以分化为成纤维细胞。如下式所示:
3)材料属性求解子模块
骨折愈合过程中,随着组织分化的进行,骨折区域组织的材料属性会发生变化。对软骨组织而言,新生成的软骨组织和钙化后的软骨组织材料属性发生改变。对骨组织而言,不成熟的骨组织和成熟的骨组织材料属性也有所不同。所以需对组织的材料属性进行求解。不同的组织由以下主要成分组成:水、矿物质、基质、Ⅰ型胶原、Ⅱ型胶原和Ⅲ型胶原。所以,被破坏的组织,不成熟骨组织,软骨组织,钙化的软骨组织,纤维结缔组织和肉芽组织的材料属性可由下式表示:
E=2000p矿化+430pⅠ型胶原+200pⅡ型胶原+100pⅢ型胶原+0.7p基质 (26)
υ=0.33p矿化+0.48(pⅠ型胶原+pⅡ型胶原+pⅢ型胶原)+0.49p基质 (27)
式中,p矿化为所含矿物质比例,pⅠ型胶原为所含Ⅰ型胶原比例;pⅡ型胶原为所含Ⅱ型胶原比例;pⅢ型胶原为所含Ⅲ型胶原比例;p基质为所含基质比例。
成熟骨材料属性由表观密度表示,如下式所示:
E骨=2014·ρ2.5,υ骨=0.2当(ρ≤1.2g/cc) (28)
E骨=1763·ρ3.2,υ骨=0.32当(ρ≥1.2g/cc) (29)
本实施方式的其他组成与连接关系与具体实施方式一至四之一相同。
具体实施方式六:如图5所示,本实施方式中,所述的骨痂生长建模模块,实现其功能的具体过程为:
在骨折愈合过程中,骨痂的尺寸和形状随时间发生改变。影响骨痂尺寸和形状的原因有两个:骨髓间充质干细胞增殖和软骨细胞肥大化。可由下式表示:
式中,为骨髓间充质干细胞增殖过程中骨痂生长速率;v为骨痂生长速率;为软骨细胞肥大化过程中骨痂生长速率;
其中,骨髓间充质干细胞增殖所引起的骨痂尺寸和形状的改变由下式表示:
式中,为最大骨髓间充质干细胞浓度;
软骨细胞肥大化引起的骨痂尺寸和形状的改变由下式表示:
式中,c软骨为软骨浓度,单位为软骨数目/mm3;
其中,f生长(Ψ,t)可由下式表示:
式中,k肥大、Ψ钙化、为与软骨细胞肥大化相关的常数。
本实施方式其他组成与连接关系与具体实施方式一至五之一相同。
具体实施方式七:本实施方式中,所述的显示模块,实现其功能的具体过程为:
将计算机的显示设备连接到计算机,用于进行视觉上的交互,进而判断骨折愈合的状态。
本实施方式中其他组成与连接关系与具体实施方式一至六之一相同。
Claims (1)
1.一种模拟骨折愈合过程的仿真系统,其特征在于,所述系统包括:
骨折部位三维几何建模模块(1)、骨折部位有限元建模模块(2)、骨折部位生物力学建模模块(3)、骨折部位细胞进化建模模块(4)、骨痂生长建模模块(5)和显示模块(6);
骨折部位三维几何建模模块(1),用于根据导入的二维断层扫描图像数据,经过图像预处理后进行骨折部位的三维表面几何模型的建立;
所述的骨折部位三维几何建模模块(1),实现其功能的具体过程为:
采用基于分割的三维医学影像表面重建算法对图像进行三维表面重构,通过阈值筛选、交互式分割和三维重建过程得到三维表面几何模型;
所述的影像由影像设备CT得到,数据存储格式为DICOM;
骨折部位有限元建模模块(2),用于对得到的骨折部位三维表面几何模型进行网格划分,实现连续的三维几何模型的离散化,得到节点坐标和单元编号;还用于将细胞浓度和组织体积分数存储到单元节点中,与节点坐标共同组成节点信息;单元编号和节点信息共同构成所需要的骨折部位有限元模型;
所述的骨折部位有限元建模模块(2),实现其功能的具体过程为:
将建立的三维表面几何模型进行网格划分,使连续的三维几何模型离散化,得到节点坐标和单元编号;将细胞浓度和组织体积分数存储到单元节点中,与节点坐标共同构成节点信息;由节点信息和单元编号组成所需要的骨折部位三维有限元模型;
所述的网格划分包括面网格划分和体网格划分两个步骤;面网格划分过程用于将三维表面模型进行优化,包括:表面模型优化,平滑处理,修补漏洞;表面模型的优化通过减小表面模型的三角面片来实现,该过程只需将相邻的两个顶点合并到一个新的顶点上,并延续原有的拓扑关系;平滑处理的过程中,对三维的面网格模型进行去噪;修补漏洞的过程中,通过将模型当中的空洞提取成空间多边形,然后对空洞多边形进行三角化的方法实现;体网格划分的过程是将面网格模型进行拉伸、旋转步骤来实现的;
骨折部位生物力学建模模块(3),用于建立骨折部位生物力学模型,求解骨痂单元力学刺激;
所述的骨折部位生物力学建模模块(3),实现其功能的具体过程为:
将骨折部位看成是多孔弹性材料,则有如下关系式:
a.固体基质、液相和总的应力应变关系如下:
σs=-φspI+σE (1)
σf=-φfpI (2)
σt=σs+σf=-pI+σE (3)
式中,σs、σf、σt分别为固相、液相和总的应力张量;p为液体压力;φs、φf分别为固相和液相体积分数;σE为有效应力张量;I为单位张量;
线弹性材料的有效应力张量可以表示为:
σE=Cε (4)
式中,σE为有效应力张量;C为刚度张量;ε为总的弹性应变;
刚度张量由下式表示:
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&upsi;</mi>
<mo>)</mo>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>&upsi;</mi>
<mo>)</mo>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mfrac>
<mi>&upsi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mi>&upsi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mi>&upsi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mfrac>
<mi>&upsi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mi>&upsi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mi>&upsi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>&upsi;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>&upsi;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>&upsi;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&upsi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,E为弹性模量;υ为泊松比;
b.考虑到两相均不可压缩性和各向同性,多孔弹性模型的连续性方程为:
▽·vs+▽·(φf(vf-vs))=0 (6)
式中,φf为液相的体积分数;vs、vf分别为固相和液相的速度向量;
c.固相和液相的动量方程如下:
<mrow>
<msup>
<mi>&pi;</mi>
<mi>s</mi>
</msup>
<mo>=</mo>
<mo>-</mo>
<msup>
<mi>&pi;</mi>
<mi>f</mi>
</msup>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&phi;</mi>
<mi>f</mi>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>k</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>v</mi>
<mi>f</mi>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mi>s</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
▽·σs+πs=0 (8)
▽·σf+πf=0 (9)
▽·σE+▽p=0 (10)
式中,πs、πf分别为固相和液相的体力向量;φf为液相体积分数;k为渗透率;vf、vs分别为固相和液相的速度向量;σs、σf、σE分别为固相、液相和有效的应力张量;p为液体压力;
d.偏应变张量的第二不变量可由下式表示:
<mrow>
<mi>&Psi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Ψ(x,t)为骨痂单元力学刺激;J2为偏应变张量第二不变量;ε1、ε2、ε3为主应变,εoct为八面体应变;
其中,八面体应变εoct可由下式求得:
<mrow>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
通过有限元法求解上述方程,得到偏应变张量的第二不变量J2,并以此作为骨痂单元的力学刺激Ψ(x,t);
骨折部位细胞进化建模模块(4),用于对骨折部位细胞的迁移,增殖,分化,凋亡进行建模,分析骨折部位的细胞进化;
所述的骨折部位细胞进化建模模块(4),实现其功能的具体过程为:
细胞的进化包括细胞的迁移、增殖、分化和凋亡;骨折愈合过程中所涉及到的细胞包括:骨髓间充质干细胞、成纤维细胞、软骨细胞和骨细胞;涉及到的组织包括肉芽组织、纤维结缔组织、软骨组织和骨组织;
细胞进化建模模块分为骨髓间充质干细胞进化建模子模块、骨细胞形成建模子模块和材料属性求解子模块;
1)骨髓间充质干细胞进化建模子模块
在骨折愈合初期,骨折部位组织被破坏,骨髓间充质干细胞侵入骨折部位并进行增殖和分化;根据骨痂单元力学刺激的不同,骨髓间充质干细胞可分化为成纤维细胞、软骨细胞和骨细胞;当骨痂单元力学刺激过高时,骨髓间充质干细胞发生凋亡;当骨髓间充质干细胞的浓度达到饱和状态,细胞浓度增加;
因此,骨髓间充质干细胞的进化模型可由下式表示:
式中,c干细胞为骨髓间充质干细胞浓度,单位为细胞数目/mm3;V破坏为被破坏组织的体积分数;Ψ为骨痂单元力学刺激;α增殖、Ψ增殖为骨髓间充质干细胞分化常数;D(V破坏)为扩散系数,单位为mm2/天;
其中扩散系数D(V破坏)下式表示:
式中,D0为常数;V破坏为被破坏组织的体积分数;
骨髓间充质干细胞的分化可由下式表示:
式中,o分化(Ψ,t)、c分化(Ψ,t)、l分化(Ψ,t)为骨髓间充质干细胞向骨细胞、软骨细胞、成纤维细胞分化;c干细胞为骨髓间充质干细胞发生凋亡;为每种细胞类型的成熟时间;Ψ极限、Ψ骨、Ψ软骨、Ψ纤维、Ψ凋亡为每种细胞类型的力学刺激极限;
骨细胞成熟时间可由下式表示:
当Ψ极限<Ψ<Ψ骨 (16)
式中,t成熟为成熟时间常数;
软骨细胞成熟时间可由下式表示:
当Ψ骨<Ψ<Ψ软骨 (17)
成纤维细胞成熟时间可由下式表示:
当Ψ软骨<Ψ<Ψ纤维 (18)
式中,t极限为时间极限;
2)骨细胞形成建模子模块
骨细胞的形成过程包括膜内骨化和软骨骨化两种途径:膜内骨化是指间充质干细胞直接分化为骨细胞;软骨骨化需经过软骨生成之后才能进行软骨骨化过程;
所以骨细胞的形成可由下式表示:
式中,c骨为骨细胞浓度;o分化(Ψ,t)为由膜内骨化形成的骨细胞;e分化(Ψ,t)为由软骨骨化形成的骨细胞;
其中,由膜内骨化形成的骨细胞o分化(Ψ,t)可由下式表示:
式中,D骨为扩散常数,单位为mm2/天;c骨为骨浓度,单位为细胞数目/mm3;Ψ极限、Ψ骨为膜内骨化过程中骨痂单元力学刺激极限值;为膜内骨化过程中最小骨浓度;为骨髓间充质干细胞分化为骨细胞所用成熟时间;
由软骨骨化形成骨细胞e分化(Ψ,t)可由下式表示:
式中,c骨为骨细胞浓度,单位为细胞数目/mm3;c软骨为软骨细胞浓度,单位为细胞数目/mm3;p为被矿化的矿化的软骨的百分比;p最小、为常数;
软骨细胞浓度与骨髓间充质干细胞的分化和软骨细胞的骨化有关,可由下式表示:
式中,g分化(Ψ,t)为由骨髓间充质干细胞分化得到的软骨细胞;g骨化(Ψ,t)为软骨骨化过程;
其中,由骨髓间充质干细胞分化得到的软骨细胞g分化(Ψ,t)可由下式表示:
软骨骨化过程g骨化(Ψ,t)可由下式表示:
g骨化=-c分化 (24)
骨髓间充质干细胞除了可以分化为骨细胞和软骨细胞,还可以分化为成纤维细胞;如下式所示:
3)材料属性求解子模块
骨折愈合过程中,随着组织分化的进行,骨折区域组织的材料属性会发生变化;对软骨组织而言,新生成的软骨组织和钙化后的软骨组织材料属性发生改变;对骨组织而言,不成熟的骨组织和成熟的骨组织材料属性也有所不同;所以需对组织的材料属性进行求解;不同的组织由以下主要成分组成:水、矿物质、基质、Ⅰ型胶原、Ⅱ型胶原和Ⅲ型胶原;所以,被破坏的组织,不成熟骨组织,软骨组织,钙化的软骨组织,纤维结缔组织和肉芽组织的材料属性可由下式表示:
E=2000p矿化+430pⅠ型胶原+200pⅡ型胶原+100pⅢ型胶原+0.7p基质 (26)
υ=0.33p矿化+0.48(pⅠ型胶原+pⅡ型胶原+pⅢ型胶原)+0.49p基质 (27)
式中,p矿化为所含矿物质比例,pⅠ型胶原为所含Ⅰ型胶原比例;pⅡ型胶原为所含Ⅱ型胶原比例;pⅢ型胶原为所含Ⅲ型胶原比例;p基质为所含基质比例;
成熟骨材料属性由表观密度表示,如下式所示:
E骨=2014·ρ2.5,υ骨=0.2当(ρ≤1.2g/cc) (28)
E骨=1763·ρ3.2,υ骨=0.32当(ρ≥1.2g/cc) (29)
骨痂生长建模模块(5),用于对骨痂形状和尺寸随时间的变化关系进行建模,分析骨痂随时间的变化情况;
所述的骨痂生长建模模块(5),实现其功能的具体过程为:
在骨折愈合过程中,骨痂的尺寸和形状随时间发生改变;影响骨痂尺寸和形状的原因有两个:骨髓间充质干细胞增殖和软骨细胞肥大化;可由下式表示:
式中,为骨髓间充质干细胞增殖过程中骨痂生长速率;v为骨痂生长速率;为软骨细胞肥大化过程中骨痂生长速率;
其中,骨髓间充质干细胞增殖所引起的骨痂尺寸和形状的改变由下式表示:
式中,为最大骨髓间充质干细胞浓度;
软骨细胞肥大化引起的骨痂尺寸和形状的改变由下式表示:
式中,c软骨为软骨浓度,单位为软骨数目/mm3;
其中,f生长(Ψ,t)可由下式表示:
式中,k肥大、Ψ钙化、为与软骨细胞肥大化相关的常数;
显示模块(6),用于观察骨折的愈合情况,实现用户与系统之间的视觉交互;
所述的显示模块(6),实现其功能的具体过程为:
将计算机的显示设备连接到计算机,用于进行视觉上的交互,进而判断骨折愈合的状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611086920.5A CN106777584B (zh) | 2016-12-01 | 2016-12-01 | 一种模拟骨折愈合过程的仿真系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611086920.5A CN106777584B (zh) | 2016-12-01 | 2016-12-01 | 一种模拟骨折愈合过程的仿真系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777584A CN106777584A (zh) | 2017-05-31 |
CN106777584B true CN106777584B (zh) | 2017-12-08 |
Family
ID=58914089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611086920.5A Expired - Fee Related CN106777584B (zh) | 2016-12-01 | 2016-12-01 | 一种模拟骨折愈合过程的仿真系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777584B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107610781B (zh) * | 2017-08-28 | 2019-04-09 | 哈尔滨理工大学 | 一种基于组织氧气环境以及力学环境的骨折愈合仿真方法 |
CN108511076B (zh) * | 2018-04-09 | 2019-10-25 | 哈尔滨理工大学 | 一种基于力学刺激和生物联合刺激的骨折愈合仿真系统 |
CN108565027B (zh) * | 2018-04-09 | 2019-12-10 | 哈尔滨理工大学 | 一种模拟骨折愈合过程的仿真系统 |
CN108536985A (zh) * | 2018-04-20 | 2018-09-14 | 河海大学常州校区 | 基于骨折愈合过程的内固定参数优化治疗个性化建模方法 |
CN113361182A (zh) * | 2021-07-02 | 2021-09-07 | 哈尔滨理工大学 | 一种基于免疫系统作用的骨折愈合仿真方法 |
CN116741387A (zh) * | 2023-06-26 | 2023-09-12 | 强联智创(北京)科技有限公司 | 由计算机实施的介入耗材的模拟方法及相关产品 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020186875A1 (en) * | 2001-04-09 | 2002-12-12 | Burmer Glenna C. | Computer methods for image pattern recognition in organic material |
CN104915519B (zh) * | 2015-06-30 | 2018-05-11 | 中国人民解放军第三军医大学第二附属医院 | 一种颅颌面骨模型建立方法及装置 |
CN105550461B (zh) * | 2015-12-30 | 2018-10-02 | 哈尔滨理工大学 | 一种基于断端微动和血供的骨折愈合仿真系统 |
-
2016
- 2016-12-01 CN CN201611086920.5A patent/CN106777584B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106777584A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777584B (zh) | 一种模拟骨折愈合过程的仿真系统 | |
Keeve et al. | Anatomy-based facial tissue modeling using the finite element method | |
CN105303605B (zh) | 一种基于力反馈的骨外科手术仿真系统 | |
CN102087676B (zh) | 一种基于孔隙网络模型的仿生骨支架设计方法 | |
CN105550461B (zh) | 一种基于断端微动和血供的骨折愈合仿真系统 | |
CN106777583B (zh) | 一种基于细胞活动的骨折愈合仿真方法 | |
CN106227993B (zh) | 一种基于生物学机理的骨折愈合动态过程仿真方法 | |
CN106777582B (zh) | 一种基于组织分化的长骨骨折愈合仿真系统 | |
CN105913718A (zh) | 一种埋线提升整形手术仿真系统 | |
CN108511076B (zh) | 一种基于力学刺激和生物联合刺激的骨折愈合仿真系统 | |
CN110443005A (zh) | 一种全膝关节有限元建模方法 | |
CN107610781A (zh) | 一种基于组织氧气环境以及力学环境的骨折愈合仿真方法 | |
Guo et al. | Tensor-mass Model based real-time simulation of vessel deformation and force feedback for the interventional surgery training system | |
Szczerba et al. | Computational model of flow–tissue interactions in intussusceptive angiogenesis | |
Szczerba et al. | Macroscopic modeling of vascular systems | |
CN108565027B (zh) | 一种模拟骨折愈合过程的仿真系统 | |
Szczerba et al. | Simulating vascular systems in arbitrary anatomies | |
Qiu et al. | Biomechanical analysis of reduction malarplasty with L-shaped osteotomy | |
CN113223724A (zh) | 一种腰椎后纵韧带六面体有限元模型的建模方法 | |
Wang et al. | A surrogate mechanostatistical microstructural model to inform whole hip cortical bone remodelling | |
Park | Topology Optimization to Design Bone Replacement Shapes in Craniofacial Reconstructive Surgery: Design, Simulation and Experimental Validation | |
Lei et al. | Study on circumpelvic muscle deformation and dynamic simulation of pelvic fracture reduction | |
Alim | Simulation of fracture strength improvements of a human proximal femur using finite element analysis. | |
Ho et al. | The Congruent Facelift | |
Xu et al. | Biomechanical Modeling and Simulation of Head |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171208 Termination date: 20191201 |