CN108920875A - 一种基于磁性氧化模型的虚拟切割算法 - Google Patents
一种基于磁性氧化模型的虚拟切割算法 Download PDFInfo
- Publication number
- CN108920875A CN108920875A CN201810853504.6A CN201810853504A CN108920875A CN 108920875 A CN108920875 A CN 108920875A CN 201810853504 A CN201810853504 A CN 201810853504A CN 108920875 A CN108920875 A CN 108920875A
- Authority
- CN
- China
- Prior art keywords
- model
- formula
- soft tissue
- equation
- force
- 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.)
- Granted
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 23
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 62
- 210000004872 soft tissue Anatomy 0.000 claims abstract description 50
- 229910052742 iron Inorganic materials 0.000 claims abstract description 31
- 238000005520 cutting process Methods 0.000 claims abstract description 20
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 239000000126 substance Substances 0.000 claims abstract description 18
- 210000001519 tissue Anatomy 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 238000009877 rendering Methods 0.000 claims abstract description 7
- 230000003647 oxidation Effects 0.000 claims description 29
- 238000007254 oxidation reaction Methods 0.000 claims description 29
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N iron oxide Inorganic materials [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 claims description 11
- NDLPOXTZKUMGOV-UHFFFAOYSA-N oxo(oxoferriooxy)iron hydrate Chemical compound O.O=[Fe]O[Fe]=O NDLPOXTZKUMGOV-UHFFFAOYSA-N 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000001125 extrusion Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 8
- 206010003694 Atrophy Diseases 0.000 claims description 7
- 230000037444 atrophy Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 208000015181 infectious disease Diseases 0.000 claims description 5
- JEIPFZHSYJVQDO-UHFFFAOYSA-N iron(III) oxide Inorganic materials O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 230000001052 transient effect Effects 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000006479 redox reaction Methods 0.000 claims description 3
- 239000002994 raw material Substances 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 10
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 102000016942 Elastin Human genes 0.000 description 1
- 108010014258 Elastin Proteins 0.000 description 1
- 206010028289 Muscle atrophy Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229920002549 elastin Polymers 0.000 description 1
- 230000002779 inactivation Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- YOBAEOGBNPPUQV-UHFFFAOYSA-N iron;trihydrate Chemical compound O.O.O.[Fe].[Fe] YOBAEOGBNPPUQV-UHFFFAOYSA-N 0.000 description 1
- 230000020763 muscle atrophy Effects 0.000 description 1
- 201000000585 muscular atrophy Diseases 0.000 description 1
- 230000001590 oxidative effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Epidemiology (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Prostheses (AREA)
Abstract
本发明提出了一种基于磁性氧化模型的虚拟切割算法，具体包括如下步骤：步骤1)建立基于无网格的弹性软组织控制方程；步骤2）通过准线性黏弹性模型并入非线性粘弹性，得到应力应变行为的非线性描述；步骤3）加入磁性氧化模型的网格模型：在网格与手术刀接触面的背侧纳入铁单质；在网格与手术刀接触面实时产生网格模型，并在网格模型下检测切割条件产生切口；步骤4）对所述网格模型进行实时状态检测，得到接触面上点的位移与实时点的位置图；步骤5）实时对组织衰败部分进行适当的渲染渲染。有益效果：在现有网格模型的基础上纳入了磁性氧化模型，本方法虚拟的手术刀的挤压力更有利于手术医师把握手感。
Description
技术领域
本发明属于医疗模拟技术领域，尤其涉及一种基于磁性氧化模型的虚拟切割算法。
背景技术
近年来，随着虚拟现实技术的不断发展，通过虚拟现实平台模拟手术训练得到可能。虚拟手术给予医护人员极大的便利，他们可以使用该技术进行反复的模拟练习，提升自己的技术。目前的虚拟手术通过建模、渲染、计算将模拟手术的过程呈现在平台上，使用了有限元模型、无网格模型等。然而，目前虚拟手术技术考虑的方面仍不全面，真实性仍旧不够。在做力的反馈时，血肉对刀具的挤压力以及该过程中的粘弹性对手术的操作有着较大的影响；而血肉与空气接触过久，容易引发感染，在虚拟现实平台中引入这一点，可以使医生把握手术速度。
发明内容
本发明目的在于克服现有技术的不足，本发明提出一种基于磁性氧化模型的网格模型的虚拟切割算法，该方法能够在不降低变形精度的情况下，通过改进后的网格模型更好地实现力的反馈并能使手术医生把控手术进程，所得到的模型更加具有真实性，具体由以下技术方案实现：
所述基于磁性氧化模型的虚拟切割算法，具体包括如下步骤：
步骤1)建立基于无网格的弹性软组织控制方程；
步骤2)通过准线性黏弹性模型并入非线性粘弹性，得到如式(1)的应力应变行为的非线性描述；
其中，σ为应力，ε是应变张量，t表示某一时刻，是考虑与应变相关的弹性非线性的函数，E(t-τ)为松弛函数，τ为一常数。
步骤3)加入磁性氧化模型的网格模型：在网格与手术刀接触面的背侧纳入铁单质；在网格与手术刀接触面实时产生网格模型，并在网格模型下检测切割条件产生切口；
步骤4)对所述网格模型进行实时状态检测，得到接触面上点的位移与实时点的位置图；
步骤5)实时对组织衰败部分进行适当的渲染渲染。
所述基于磁性氧化模型的虚拟切割算法的进一步设计在于，所述步骤1)中当施加力时，用X＝(x,y,z)T表示点的原始材料坐标，设定将所述点移动到新的位置X'＝(x',y',z')T，导致软组织变形。
所述基于磁性氧化模型的虚拟切割算法的进一步设计在于，所述软组织变形的方式替换为通过位移场U＝(u(x,y,z),v(x,y,z),w(x,y,z))实现软组织的变形，原始位置与变形得到的新位置之间的关系、相互间作用力的表示分别如式(2)、式(3)：
U＝X'-X (2)
其中，U为位移场，X'表示软组织变形得到的新位置，X表示软组织变形的原始位置，f为相互间的作用力，E为弹性系数，是位移场的梯度。
所述基于磁性氧化模型的虚拟切割算法的进一步设计在于，所述步骤2)中设定在t＝t0的时刻一外力作用于软组织，引起瞬时应变，该应变用阶跃函数表示如式(4)，
ε(t)＝H(t)ε0 (4)
其中，H(t)是单位阶跃函数，ε0表示在t0时刻的应变张量；
阶跃函数的导数是脉冲函数如式(5)：
其中，δ(t)是单位脉冲函数；
根据解决奇异函数的卷积积分的策略，将t＝t0分解为t＝t0 -和t＝t0 +，且当t≠t0时，脉冲函数为零，因此从0到t0 -的脉冲函数的卷积积分为零，结合式(1)得到弹性模量的平衡值与瞬态模量，如式(6)：
其中，为Prony系列的松弛函数，pi和τi是可以从实验数据中获得的两个系数；
设T的模拟时间被分成n个时间片，即结合式(1)，得到t+Δt时刻的等式如式(7)：
根据式(8)在非线性粘弹性软组织的任意点i上施加力，
根据式(9)使用Leap Frog方案计算位移：
对材料给予系数后，得到非线性粘弹性软组织的应力应变响应，并使用时间增量算法计算基于前一次已知解的应力，作用在任意点上的力将从应力中获得；
最后，遵循牛顿的运动定律，建立软组织的运动控制方程。
所述基于磁性氧化模型的虚拟切割算法的进一步设计在于，所述步骤3)所述网格模型中使用碰撞检测得到的点适应切割面，切割面的获取具体为：
设定有n个碰撞点Pi(xi,yi,zi)，设定Pi(xi,yi,zi)构成的切割面符合要求，如式(10)：
式(10)中定义z＝a0x+a1y+a2，A、B、C、D均为系数，且C≠0，满足Ax+By+Cz+D＝0；
当时，式(10)达到最小值，因此得到式(11)、式(12)：
求解式(12)，得到系数a0,a1,a2，进而得到切割面。
所述基于磁性氧化模型的虚拟切割算法的进一步设计在于，所述步骤4)中通过对刀的排斥力模拟生物软组织对手术刀的挤压力，通过铁生锈模拟生物组织受外界因素影响而引发的感染腐烂。
所述基于磁性氧化模型的虚拟切割算法的进一步设计在于，通过式(13)、式(14)分别限定排斥力、组织萎缩所产生的形变力与质量的关系，
F1＝k1m1 (13)
F2＝k2m2 (14)
其中，F1为当前铁单质对手术刀的排斥力，k1为由实验得出的常数，表示铁单质质量与产生的排斥力之间的线性关系，m1表示当前铁单质的质量；F2为组织萎缩所产生的形变力，k2为由实验得出的常数，表示三氧化二铁质量与萎缩力之间的线性关系，m2表示当前三氧化二铁的质量；
根据铁氧化的化学反应方程式得到铁单质的氧化还原反应随时间的转化率：
α＝f(t) (15)
通过所述转化率，计算出铁单质与三氧化二铁的质量，并结合式(13)、式(14)、式(8)以及式(9)得到接触面上点的位移，得到实时点的位置图。
本发明的优点如下：
(1)本发明的方法在现有网格模型的基础上纳入了磁性氧化模型，由于生物软组织对虚拟手术刀的挤压力是生物软组织的基本特性，因此存在挤压力进行模拟切割更加符合实际，本发明虚拟手术刀的挤压力更有利于手术医师把握手感。
(2)该方法以铁单质随时间缓慢氧化，即生锈，模拟因产生切口而裸露的生物软组织接触空气等因素引起的感染，更加符合实际，利于手术医师把握手术进度，降低伤口感染组织萎缩的几率。同时，可以有效地对伤亡人员进行急救。
(3)该方法提出的磁性氧化模型具有较强的鲁棒性，适用于绝大多数虚拟手术，有利于手术医师在虚拟现实系统中磨练手术技能。
附图说明
图1是切割示意图。
图2是磁性氧化网格模型图的主视图。
图3是磁性氧化网格模型图的左视图。
具体实施方式
结合具体实施例与附图对本发明的技术方案进一步说明。
如图1，本实施例提供的基于磁性氧化模型的虚拟切割算法，其主要步骤主要包括建立基于无网格的弹性软组织控制方程、并入非线性粘弹性、检测切割条件产生切口、实时状态检测、实时渲染，具体如下：
步骤1)建立基于无网格的弹性软组织控制方程。
在三维空间中，连续弹性理论被用来表示连续物体如软组织的行为。遵循无网格法的原理，软组织典型地通过一系列点离散化。当施加力时，用X＝(x,y,z)T表示点的原始材料坐标，上述点将移动到新的位置X'＝(x',y',z')T，导致软组织变形。
本实施例还可以用位移场U＝(u(x,y,z),v(x,y,z),w(x,y,z))来描述软组织的变形。原始位置与变形得到的新位置之间的关系是：
U＝X'-X (1)
力可以被表示为：
其中，U为位移场，X'表示软组织变形得到的新位置，X表示软组织变形的原始位置，f为相互间的作用力，E为弹性系数，是位移场的梯度。
步骤2：并入非线性粘弹性。
软组织显示出良好的弹性(伸展15％无损伤)和粘度，因为它主要由两种类型的蛋白质：胶原蛋白和弹性蛋白组成。在生物力学中，有几种模型描述了软组织的粘弹行为，最常用的模型是准线性黏弹性(QLV)模型。由于时间依赖性与QLV中的应变或应力无关，应力应变行为的非线性描述如下：
其中，σ为应力，ε是应变张量，t表示某一时刻，是考虑与应变相关的弹性非线性的函数，E(t-τ)为松弛函数，τ为一常数。
假设在t＝t0的时刻有一个外力作用于软组织，则会引起瞬时应变，该应变用阶跃函数表示：
ε(t)＝H(t)ε0 (4)
其中，H(t)是单位阶跃函数，ε0表示在t0时刻的应变张量。阶跃函数的导数是脉冲函数，有：
其中，δ(t)是单位脉冲函数。将方程(7)代入方程(3)，同时根据解决奇异函数的卷积积分的策略，将t＝t0分解为t＝t0 -和t＝t0 +，且当t≠t0时，脉冲函数为零，所以从0到t0 -的脉冲函数的卷积积分为零，得到：
其中，为Prony系列的松弛函数，pi和τi是可以从实验数据中获得的两个系数。
如方程(6)下侧式子所示，第一项是弹性模量的平衡值，第二项是瞬态模量。此外，应力由应变的稳定响应和应变的瞬态响应组成。
与显式方案相比，隐式方案是无条件稳定的。这需要使用增量形式来计算增加的时间间隔之后的压力。假设T的模拟时间被分成n个时间片，即结合(3)，得到t+Δt时刻的等式(Δt趋近于0)：
在非线性粘弹性软组织的任意点i上施加的力为：
使用Leap Frog方案计算位移：
对材料给予适当的系数后，可以得到非线性粘弹性软组织的应力应变响应。使用时间增量算法来计算基于前一次已知解的应力。另外，作用在任意点上的力将从应力中获得。最后，遵循牛顿的运动定律，建立软组织的运动控制方程。
步骤3：检测切割条件产生切口，参见图1(图1中1表示接触面，2表示虚拟手术刀)，并实时产生改进的网格模型，参见图2(图2中3表示铁单质，2表示网格)。
在真正的切割中，在手术刀穿透软组织的过程中发生以下两个重要的相互作用。首先，当切割力增加时，软组织的变形不断发生。直到施加到软组织上的切割力超过阈值水平才会发生分离。其次，当施加在软组织上的应力达到断裂应力时，软组织将被切开，并随着切割力的急剧下降而开始快速膨胀。一旦软组织被切开，软组织将以与切割平面中心相反的方向变形。因此，本实施例把切割过程分成三个阶段。采用了包含非线性粘弹性的变形模型，更切合实际地表示了变形。
手术刀通常只是沿着刀片的方向切割软组织，所以在计算软组织变形之前，必须确定切割面。由于切割面是虚构的，使用碰撞检测得到的点来适应切割面。在几何中，一个平面的方程可以表示为：
Ax+By+Cz+D＝0 (10)
当C≠0时，有：
为了方便，定义有：
z＝a0x+a1y+a2 (12)
假设有n个碰撞点Pi(xi,yi,zi)，设定Pi(xi,yi,zi)构成的切割面符合要求：
当时，上式达到最小值。因此：
求解上式，可以得到系数a0,a1,a2，然后得到切割面。
加入磁性氧化模型的网格模型：在网格-手术刀接触面的背侧纳入适当质量(实验得出)的铁单质；在网格-手术刀接触面实时产生网格模型。
步骤4：进行实时状态检测。
生物软组织在手术刀切割时，会对手术刀产生一个较小的挤压力；而切割产生的开口接触空气等容易发生感染，即组织萎缩。在虚拟手术系统中加入这两点，能更好的模拟真实，使手术医师能更好的磨练技术。
本文通过铁对刀(金属间)的排斥力，铁生锈分别模拟生物软组织对手术刀的挤压力与生物组织受外界因素影响而引发的感染腐烂。
因排斥力与组织萎缩所产生的形变力往往比较小，故将两种力与质量的关系视为线性关系，定义：
F1＝k1m1 (16)
其中，F1为当前铁单质对手术刀的排斥力，k1为由实验得出的常数，表示铁单质质量与产生的排斥力之间的线性关系，m1表示当前铁单质的质量。
F2＝k2m2 (17)
其中，F2为当前由于铁单质被氧化，生成三氧化二铁，模拟局部组织腐烂，背离切口施加该力模拟产生的肌肉萎缩状态，k2为由实验得出的常数，表示三氧化二铁质量与萎缩力之间的线性关系，m2表示当前三氧化二铁的质量。
不考虑铁单质与三氧化二铁互相掺杂，两者之间有明显的界限，产生三氧化二铁的部分无挤压力，模拟组织失活。
铁氧化涉及的化学反应方程式：
2Fe+O2+2H2O＝2Fe(OH)2
4Fe(OH)2+O2+2H2O＝4Fe(OH)3
2Fe(OH)3＝Fe2O3+3H2O
通过实验，得到铁单质的氧化还原反应随时间的转化率：
α＝f(t) (18)
通过实时转化率，可以计算出铁单质与三氧化二铁的质量。
将(16)(17)代入(8)(9)可以得到接触面上点的位移，得到实时点的位置图。
步骤5：实时渲染。组织衰败部分进行适当的渲染，达到腐烂效果。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (7)
1.一种基于磁性氧化模型的虚拟切割算法，其特征在于具体包括如下步骤：
步骤1)建立基于无网格的弹性软组织控制方程；
步骤2)通过准线性黏弹性模型并入非线性粘弹性，得到如式(1)的应力应变行为的非线性描述；
其中，σ为应力，ε是应变张量，t表示某一时刻，是考虑与应变相关的弹性非线性的函数，E(t-τ)为松弛函数，τ为一常数。
步骤3)加入磁性氧化模型的网格模型：在网格与手术刀接触面的背侧纳入铁单质；在与手术刀接触面实时产生网格模型，并在网格模型下检测切割条件产生切口；
步骤4)对所述网格模型进行实时状态检测，得到接触面上点的位移与实时点的位置图；
步骤5)实时对组织衰败部分进行适当的渲染渲染。
2.根据权利要求1所述的基于磁性氧化模型的虚拟切割算法，其特征在于所述步骤1)中当施加力时，用X＝(x,y,z)T表示点的原始材料坐标，设定将所述点移动到新的位置X'＝(x',y',z')T，导致软组织变形。
3.根据权利要求2所述的基于磁性氧化模型的虚拟切割算法，其特征在于所述软组织变形的方式替换为通过位移场U＝(u(x,y,z),v(x,y,z),w(x,y,z))实现软组织的变形，原始位置与变形得到的新位置之间的关系、相互间作用力的表示分别如式(2)、式(3)：
U＝X'-X (2)
其中，U为位移场，X'表示软组织变形得到的新位置，X表示软组织变形的原始位置，f为相互间的作用力，E为弹性系数，是位移场的梯度。
4.根据权利要求2所述的基于磁性氧化模型的虚拟切割算法，其特征在于所述步骤2)中设定在t＝t0的时刻一外力作用于软组织，引起瞬时应变，该应变用阶跃函数表示如式(4)，
ε(t)＝H(t)ε0 (4)
其中，H(t)是单位阶跃函数，ε0表示在t0时刻的应变张量；
阶跃函数的导数是脉冲函数如式(5)：
其中，δ(t)是单位脉冲函数；
根据解决奇异函数的卷积积分的策略，将t＝t0分解为t＝t0 -和t＝t0 +，且当t≠t0时，脉冲函数为零，因此从0到t0 -的脉冲函数的卷积积分为零，结合式(1)得到弹性模量的平衡值与瞬态模量，如式(6)：
其中，为Prony系列的松弛函数，pi和τi是可以从实验数据中获得的两个系数；
设T的模拟时间被分成n个时间片，即结合式(1)，得到t+Δt时刻的等式如式(7)：
根据式(8)在非线性粘弹性软组织的任意点i上施加力，
根据式(9)使用Leap Frog方案计算位移：
对材料给予系数后，得到非线性粘弹性软组织的应力应变响应，并使用时间增量算法计算基于前一次已知解的应力，作用在任意点上的力将从应力中获得；
最后，遵循牛顿的运动定律，建立软组织的运动控制方程。
5.根据权利要求1所述的基于磁性氧化模型的虚拟切割算法，其特征在于所述步骤3)所述网格模型中使用碰撞检测得到的点适应切割面，切割面的获取具体为：
设定有n个碰撞点Pi(xi,yi,zi)，设定Pi(xi,yi,zi)构成的切割面符合要求，如式(10)：
式(10)中定义z＝a0x+a1y+a2，A、B、C、D均为系数，且C≠0，满足Ax+By+Cz+D＝0；
当时，式(10)达到最小值，因此得到式(11)、式(12)：
求解式(12)，得到系数a0,a1,a2，进而得到切割面。
6.根据权利要求1所述的基于磁性氧化模型的虚拟切割算法，其特征在于所述步骤4)中通过对刀的排斥力模拟生物软组织对手术刀的挤压力，通过铁生锈模拟生物组织受外界因素影响而引发的感染腐烂。
7.根据权利要求6所述的基于磁性氧化模型的虚拟切割算法，其特征在于通过式(13)、式(14)分别限定排斥力、组织萎缩所产生的形变力与质量的关系，
F1＝k1m1 (13)
F2＝k2m2 (14)
其中，F1为当前铁单质对手术刀的排斥力，k1为由实验得出的常数，表示铁单质质量与产生的排斥力之间的线性关系，m1表示当前铁单质的质量；F2为组织萎缩所产生的形变力，k2为由实验得出的常数，表示三氧化二铁质量与萎缩力之间的线性关系，m2表示当前三氧化二铁的质量；
根据铁氧化的化学反应方程式得到铁单质的氧化还原反应随时间的转化率：
α＝f(t) (15)
通过所述转化率，计算出铁单质与三氧化二铁的质量，并结合式(13)、式(14)、式(8)以及式(9)得到接触面上点的位移，得到实时点的位置图。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810853504.6A CN108920875B (zh) | 2018-07-29 | 2018-07-29 | 一种基于磁性氧化模型的虚拟切割方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810853504.6A CN108920875B (zh) | 2018-07-29 | 2018-07-29 | 一种基于磁性氧化模型的虚拟切割方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108920875A true CN108920875A (zh) | 2018-11-30 |
| CN108920875B CN108920875B (zh) | 2023-01-24 |
Family
ID=64417681
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810853504.6A Active CN108920875B (zh) | 2018-07-29 | 2018-07-29 | 一种基于磁性氧化模型的虚拟切割方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108920875B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007140814A (ja) * | 2005-11-17 | 2007-06-07 | Ritsumeikan | アダプティブメッシュを用いた物体の変形シミュレーション方法及び切断シミュレーション方法並びにアダプティブメッシュを用いた物体の変形シミュレーションプログラム及び切断シミュレーションプログラム |
| CN101404039A (zh) * | 2008-03-28 | 2009-04-08 | 华南师范大学 | 一种虚拟手术的方法及其装置 |
| CN106485000A (zh) * | 2016-10-11 | 2017-03-08 | 南京信息工程大学 | 基于虚拟弹簧和球面调和函数力反馈模型建模方法 |
| CN107080586A (zh) * | 2017-03-24 | 2017-08-22 | 南京信息工程大学 | 虚拟手术中用于软组织实时形变的交叉弹簧模型 |
| CN107146288A (zh) * | 2017-05-16 | 2017-09-08 | 南京信息工程大学 | 虚拟手术中支持实时按压形变的软组织模型建模方法 |
| US20180204378A1 (en) * | 2017-01-13 | 2018-07-19 | Beihang University | Method for real-time cutting of digital organ based on metaball model and hybrid driving method |
-
2018
- 2018-07-29 CN CN201810853504.6A patent/CN108920875B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007140814A (ja) * | 2005-11-17 | 2007-06-07 | Ritsumeikan | アダプティブメッシュを用いた物体の変形シミュレーション方法及び切断シミュレーション方法並びにアダプティブメッシュを用いた物体の変形シミュレーションプログラム及び切断シミュレーションプログラム |
| CN101404039A (zh) * | 2008-03-28 | 2009-04-08 | 华南师范大学 | 一种虚拟手术的方法及其装置 |
| CN106485000A (zh) * | 2016-10-11 | 2017-03-08 | 南京信息工程大学 | 基于虚拟弹簧和球面调和函数力反馈模型建模方法 |
| US20180204378A1 (en) * | 2017-01-13 | 2018-07-19 | Beihang University | Method for real-time cutting of digital organ based on metaball model and hybrid driving method |
| CN107080586A (zh) * | 2017-03-24 | 2017-08-22 | 南京信息工程大学 | 虚拟手术中用于软组织实时形变的交叉弹簧模型 |
| CN107146288A (zh) * | 2017-05-16 | 2017-09-08 | 南京信息工程大学 | 虚拟手术中支持实时按压形变的软组织模型建模方法 |
Non-Patent Citations (7)
| Title |
|---|
| 王党校等: "计算机触觉:虚拟现实环境的力触觉建模和生成", 《计算机辅助设计与图形学学报》 * |
| 程强强: "虚拟手术仿真系统中软组织切割模型研究", 《中国优秀博硕士学位论文全文数据库(博士)医药卫生科技辑》 * |
| 臧爱云等: "基于真实切割的虚拟手术中力觉模型建模方法研究", 《中国体视学与图像分析》 * |
| 谢叻等: "虚拟手术中的力学变形和力觉感知", 《医用生物力学》 * |
| 赵沁平: "虚拟现实综述", 《中国科学（F辑:信息科学）》 * |
| 郭煜等: "虚拟手术中软组织形变仿真研究综述", 《集成技术》 * |
| 魏迎梅等: "手术仿真中的碰撞检测研究", 《计算机研究与发展》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108920875B (zh) | 2023-01-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6049788B2 (ja) | 仮想道具操作システム | |
| Okamura et al. | The haptic scissors: Cutting in virtual environments | |
| CN101320526B (zh) | 一种手术预测和训练的设备及其方法 | |
| Mahvash et al. | Modeling the forces of cutting with scissors | |
| Jeon et al. | Haptic simulation of breast cancer palpation: A case study of haptic augmented reality | |
| KR20100099507A (ko) | 햅틱 인터페이스를 이용한 임플랜트 시뮬레이션 시스템 | |
| CN114462265A (zh) | 一种用于可变形物体剪切断裂的仿真方法及材料仿真方法 | |
| CN108920875B (zh) | 一种基于磁性氧化模型的虚拟切割方法 | |
| CN109598799B (zh) | 一种基于CycleGAN的虚拟切割方法 | |
| Liu et al. | Real-time simulation of virtual palpation system | |
| WO2008072756A1 (ja) | 反力提示方法、および力覚提示システム | |
| Gick et al. | Categorical variation in lip posture is determined by quantal biomechanical-articulatory relations | |
| KR101514237B1 (ko) | 광대뼈 절단 가이드 템플릿 및 그 제작방법 | |
| Siregar et al. | Numerical Solution of Delayed SIR Model of Tuberculosis with Combination of Runge Kutta Method and Taylor Series Approach | |
| CN113705039B (zh) | 一种融入生物特性的高保真度肺部变形模型及其建模方法 | |
| Tillier et al. | Three-dimensional finite element modelling for soft tissues surgery | |
| CN114693892A (zh) | 一种虚拟手术中切割流血模型的建模方法 | |
| Voß et al. | Virtual cutting of anatomical structures | |
| Nikolaev | Non-linear mass-spring system for large soft tissue deformations modeling | |
| De Casson et al. | Simulating 2D tearing phenomena for interactive medical surgery simulators | |
| CN109033641A (zh) | 一种基于硅胶愈合模型的虚拟切割算法 | |
| Trantina et al. | Elastic-plastic fracture mechanics analysis of small cracks | |
| Wood | Structure-Property Relationships and Morphometric Effects of Different Shark Teeth on Shearing Performance | |
| Mathiyalagan et al. | Wave propagation in human tissue from transient high frequency vibrations | |
| Kim et al. | An empirical nonlinear viscoelastic model of reflective force by a layer of soft tissue |
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 |