一种软组织形变仿真方法.pdf

(19)中华人民共和国国家知识产权局*CN102044086A*(12)发明专利申请(10)申请公布号CN102044086A(43)申请公布日2011.05.04(21)申请号201010565036.6(22)申请日2010.11.30(71)申请人华北水利水电学院地址450011河南省郑州市北环路36号华北水利水电学院(72)发明人刘雪梅皇甫中民向明森闫雒恒赵振国赵晶闫新庆吴慧欣郝爱民刘明堂孙新娟杨礼波石秋华刘欢(51)Int.Cl.G06T17/00(2006.01)G06F19/00(2011.01)权利要求书4页说明书8页附图4页(54)发明名称一种软组织形变仿真方法(57)摘要本发明涉及一种基于光滑粒子流体动力学的软组织形变仿真方法，属于图形处理技术领域，该方法选取光滑粒子流体动力学法，以黏弹性力学模型来反映软组织的生物力学特性，包含以下步骤：依据黏弹性模型，构建软组织形变仿真计算相关的一系列方程；选择合适的支持域搜索策略和光滑核函数，采用粒子近似法对方程组的各相关项进行近似计算，通过显示积分法计算各粒子的密度、位置、速度等随时间的变化值；动态将粒子模型每个时间步长的状态输出到屏幕上，并进行纹理光照的渲染，显示软组织器官受力情况下的实时形变过程。该方法无需繁琐的网格计算，可提高软组织形变仿真的准确性和实时性。CN102486ACCNN110204408602044089A权利要求书1/4页1.一种软组织形变仿真方法和技术，其特征在于包含以下步骤：步骤1)：采集软组织的数据信息；步骤2)：选择黏弹性模型，构建用于软组织形变仿真的方程组：采用一个弹簧和一个黏壶并联的Kelvin黏弹性模型；首先，构建三维格式的Kelvin黏弹性应力-应变本构方程：依据Kelvin模型，各向同性材料的应力张量σ可分解成它的球形张量和偏斜张量部分，应变张量ε可分离为体积形变和等体积的形状畸变两部分：其中，α，β＝x，y，z.σαβ为应力分量，εαβ为应变分量；δαβ为Kronecker符号，σkk＝σxx+σyy+σzz和εkk＝εxx+εyy+εzz分别为体积应力和体积应变；Sαβ和eαβ分别为偏应力张量和偏应变张量的分量；根据Kelvin模型，偏应力张量和偏应变张量之间、体积应力和体积应变之间的三维黏弹性本构关系可表示为：上式中，Sαβ和eαβ分别为偏应力张量和偏应变张量的分量；σkk和εkk分别为体积应力和体积应变；E为材料的弹性模量，η为黏性系数；分别为偏应变分量和体积应变对时间的导数，即应变率，t为时间；其次，构建应变-位移几何方程：其中u，v，w为位移在三个坐标方向的分量；εxx，εyy，εzz为正应变分量；εxy，εyz，εzx为剪应变分量；然后，构建用于加速度计算的动量方程，其形式如下：其中，v为速度向量，t为时间，ρ为粒子的密度，x为坐标向量；步骤3)：依据步骤1)采集的软组织的数据信息，利用各数据点的位置向量，构建一个没有网格连接的粒子模型，并初始化所述粒子模型中各粒子的位置、质量、速度、加速度、受作用力等信息，构建粒子模型的初始化状态；2CCNN110204408602044089A权利要求书2/4页步骤4)：定义空间栅格，采用链表搜索法搜索粒子的支持域，构建支持域内的光滑核函数；核函数选取三次B-Spline光滑函数：式中，Wij为由邻近粒子pj估计粒子pi运动信息的光滑核函数，为粒子pi与pj之间的相对距离，r为粒子pi与pj之间的距离，h为光滑长度；对核函数在各方向求导，得到核函数的一阶导数为：式中，上标β＝x，y，z表示坐标方向，s，r，h的含义同式(6)，和分别表示粒子pi与pj的位置坐标向量在各方向的分量；步骤5)：应用光滑核函数W及其导数对参考粒子的支持域内所有粒子函数的加权平均近似的方法来构建步骤2)中各方程的SPH格式：构建密度方程SPH格式其中，ρi为粒子pi的密度，mj为pi的支持域内邻近粒子pj的质量，N为粒子pi的支持域内粒子总数，Wij的含义同式(6)；将动量方程(5)转换为SPH格式其中，分别为粒子pi的运动速度在各坐标方向的分量；为粒子pi应力向量的正应力分量，为粒子pi应力向量的剪应力分量；为粒子pj应力向量的正应力分量，为粒子pj应力向量的剪应力分量；ρi，ρj分别为粒子pi和pj的密度；mj，Wij，N的含义同公式(8)；3CCNN110204408602044089A权利要求书3/4页构建几何方程的SPH格式：其中，为粒子pi应变向量的正应变分量，为粒子pi应变向量的剪应变分量；ui，vi，wi为粒子pi的运动位移在各坐标方向的分量，uj，vj，wj为粒子pj的运动位移在各坐标方向的分量；ρj为粒子pj的密度，mj，Wij，N的含义同式(8)；步骤6)：用显示积分法求解步骤5)