CN102968811B - 微创血管介入手术中导丝的实时运动仿真方法 - Google Patents
微创血管介入手术中导丝的实时运动仿真方法 Download PDFInfo
- Publication number
- CN102968811B CN102968811B CN201210413201.5A CN201210413201A CN102968811B CN 102968811 B CN102968811 B CN 102968811B CN 201210413201 A CN201210413201 A CN 201210413201A CN 102968811 B CN102968811 B CN 102968811B
- Authority
- CN
- China
- Prior art keywords
- seal wire
- guide wire
- represent
- carry out
- blood vessel
- 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.)
- Active
Links
Landscapes
- Surgical Instruments (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种微创血管介入手术中导丝的实时运动仿真方法,包括以下步骤:1)采用三角形表面网格对血管进行建模;2)采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力;3)进行导丝的渲染;4)采用K-DOPS树进行导丝与血管之间碰撞检测;5)采用非迭代的约束方法计算导丝与血管之间的接触力并进行碰撞响应;6)使用Verlet积分公式来迭代更新导丝的位置与速度;7)利用拉格朗日乘数和快速投影方法来实现导丝的不可伸缩特性;8)使用力反馈设备Phantom Omni来渲染导丝给用户的作用力。与现有技术相比,本发明具有稳定可靠、实现方便、真实高效、工程应用性强等优点。
Description
技术领域
本发明涉及一种虚拟手术的模拟方法,尤其是涉及一种微创血管介入手术中导丝的实时运动仿真方法。
背景技术
随着科学技术的高速发展,虚拟手术的应用为医疗的现代化提供了越来越多的帮助。虚拟的微创血管介入手术***能够帮助实习医生快速提高手术技能,让更多的心脑血管疾病患者能够进行微创介入治疗。导丝的实时运动仿真作为微创血管介入手术***中不可缺少的一部分,成为许多学者研究的课题。在导丝实时运动仿真的研究过程中,早期的一部分学者采用几何模型对导丝进行建模,这种方法完全忽略了导丝的物理性质,然而如果完全忽略了导丝的物理性质,与其他模型之间的交互将难以处理。也有一部分学者采用质点弹簧这样简单的物理模型来对导丝进行建模,这种方法虽然计算效率很高,但是由于模型过于简单,根本无法真实地模拟导丝在血管中的运动。本文实现的微创血管介入手术中导丝的实时运动仿真对导丝的模拟有很高的要求,不仅要求导丝要有高的真实性,还对导丝与血管之间交互的模拟要有可靠的物理模型,并且由于要应用于虚拟手术,这就要求模型***达到高的更新速度,并且能达到实时。但是微创血管介入手术中的导丝的实时运动仿真对血管的模拟要求较低。一个完整的微创血管介入手术中的导丝的实时运动仿真包括了血管模拟、导丝模拟、碰撞检测及响应、导丝的不可伸缩和力反馈渲染。
对于血管的模拟,由于血管在仿真的过程中几乎不产生形变,只需要采用三角形表面网格这样简单的表面模型对血管建模就足够了。不需要利用四面体网格这样复杂的立体模型来模拟血管。这样有效地降低了计算成本和存储空间。
导丝的模拟,是微创血管介入手术中导丝的实时运动仿真中的重点部分,算法设计的好坏关系到整个仿真***的真实性、准确性、高效性。目前,大部分的方法是多刚体运动模型,弹簧模型,有限元模型,B-Spline等理论,这其中弹簧模型和 B-Spline被常用于导丝的模拟,由于其效率高能达到实时。然而弹簧模型过于简单,忽略了导丝很多重要的物理性质而导致真实感很差。B-Spline的缺点是缺乏物理支持,只是从几何的角度来仿真导丝在血管中的形态,无法知道导丝的受力情况。因此无法为后面的力反馈设备提供真实的力反馈。有限元模型能够提供很强的真实性。但是有限元方法有一个致命的弱点:实时性差。不能满足在仿真过程中与用户进行实时交互的要求。为了同时达到真实性和实时性的目标,本文采用了基尔霍夫弹性杆对导丝进行物理建模。
主流的碰撞检测方法之一是,采用层次结构(hierarchy)为基础核心,例如Octree,k-d树,BSP树等,利用包围体的hierarchy结构也被广泛地应用于各种碰撞检测算法(同时这种方法也被应用于其他领域,如:光线追踪),包围体的选择,学者们也提出了许多方法,包括,球包围体,有轴向的长方体(AABBs)另外一些碰撞检测方法,包括采用space-time bounds以及采用四维几何体来包围物体,近些年,学者们也提出基于Voronoi图的碰撞检测方法。
导丝与血管之间的碰撞响应,目前的两种主流算法是惩罚法和约束法。惩罚法的真实性和稳定性差,然而约束法需要进行多次迭代才能确保导丝不刺穿血管壁,计算成本很高。结合惩罚法与约束法的优点,本文采用一种非迭代的约束方法来计算导丝与血管之间的接触力并进行碰撞响应。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种稳定可靠、实现方便、真实高效的微创血管介入手术中导丝的实时运动仿真方法。
本发明的目的可以通过以下技术方案来实现:
一种微创血管介入手术中导丝的实时运动仿真方法,其特征在于,该方法包括以下步骤:
1)由于血管在仿真的过程中不产生形变,采用三角形表面网格对血管进行建模;
2)采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力;
3)进行导丝的渲染,相邻的两导丝质点之间绘制一个圆柱体,同时在相邻的两圆柱体之间绘制一个圆球进行连接;
4)采用K-DOPS树进行导丝与血管之间碰撞检测;
5)采用非迭代的约束方法计算导丝与血管之间的接触力并进行碰撞响应;
6)使用Verlet积分公式来迭代更新导丝的位置与速度;
7)利用拉格朗日乘数和快速投影方法来实现导丝的不可伸缩特性;
8)使用力反馈设备Phantom Omni来渲染导丝给用户的作用力。
所述的采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力的具体实现过程为:
21)将导丝描述成一条被坐标化的曲线Г={t,m1,m2},其中t表示沿着曲线的单位切线,m1,m2表示曲线的两条单位法向量,该三个正交向量构成一个局部移动坐标系{t(s),m1(s),m2(s)},其中s表示曲线的弧长;
22)利用连续的基尔霍夫弹性杆模型计算被坐标化的导丝的连续弹性能量E(Г),其中E(Г)=E弯曲(Г)+E扭曲(Г);
23)利用离散的基尔霍夫弹性杆模型对连续的弹性能量进行离散化处理,推导出导丝的离散弹性能量的公式;
24)利用拉格朗日方程对导丝的离散弹性能量进行导丝位置的求导,计算出导丝的形变力。
所述的连续弹性能量E(Г)=E弯曲(Г)+E扭曲(Г)的具体公式为:
其中二维向量ω=(ω1,ω2)T表示导丝在局部移动坐标系中的曲率向量,m=m1′m2表示导丝在局部移动坐标系中的扭曲率;其中α表示导丝的弯曲模数,β表示导丝的扭曲系数;ω1,ω2分别表示导丝在局部移动坐标系m1,m2上的曲率,ds表示对导丝的弧长进行积分。
所述的离散弹性能量E(Г)=E弯曲(Г)+E扭曲(Г)的具体公式为:
其中i=1,2,3,…n表示将导丝进行离散分段后的段号,表示第i段导丝的长度, θi表示导丝的第i个顶点的旋转角度,ωi表示第i段导丝在局部移动坐标系中的曲率向量,mi表示第i段导丝在局部移动坐标系中的扭曲率。
所述的采用非迭代的约束方法的具体实现过程为:
首先计算导丝在血管中的穿刺向量d;然后根据穿刺向量计算接触力
最后利用步骤6)提到的Verlet积分公式来进行导丝的位置更新,确保导丝仅在血管内部运动,不穿透血管壁。
所述的快速投影方法的具体实现过程为:
71)构造导丝的长度约束函数C(p1,p2)=|p1,p2|-d,其中p1,p2为相邻两导丝质点的位置坐标,d为该相邻两质点的原始长度;
72)如果长度约束函数的值超过了最低阈值,则对质点p1,p2进行位置修正并且循环执行该过程直至长度约束函数的值小于最低阈值,其中修正公式为: 因此p1,p2的修正量Δp1,Δp2为:
所述的力反馈渲染的具体实现方法为:
导丝的仿真频率为50Hz,而力反馈设备的更新频率要达到1000Hz才能提供平滑的力触觉,为了实现两者的同步,采用线性插值的方法来计算力反馈设备的力输入参数。
与现有技术相比,本发明具有以下优点:
(1)采用基尔霍夫弹性杆模型,作为导丝的模拟基础,并且采用Verlet更新方法,有效地提高了数值的稳定性和精确性,实现了真实高效的更新;
(2)采用非迭代的约束方法来计算导丝与血管之间的接触力并进行碰撞响应,该方法避免了传统的约束方法需多次迭代的问题。
(3)采用了K-DOPS树,提高了碰撞检测的效率,提高了微创血管介入手术中导丝的实时运动仿真的真实性和有效性。
(4)利用拉格朗日乘数和快速投影方法来实现导丝的不可伸缩特性,该方法能够快速、有效地实现导丝的不可拉伸。
附图说明
图1为本发明的流程示意图;
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
本实施例在一台拥有Pentium(R)Dual-Core CPU、一个NVIDIA GeForce GT
240和2G内存的个人计算机中实现,整个虚拟手术训练***采用C++语言编写,采用QT编写界面。***在实现的过程中使用的第三方库有OpenGL、Boost和CGAL。
本实施例的实施流程如图1所示,一种微创血管介入手术中导丝的实时运动仿真方法,具体包括以下步骤:
第一步,采用三角形表面网格,对血管进行建模。所述的血管三角形表面网格,首先是从人体CT数据中分割出STL格式的血管数据文件,然后利用Netgen生成三角形表面网格,最后对三角形表面网格进行平滑处理生成最终***可用的三角形表面网格。
第二步,采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力。所述的采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力的具体实现过程为:
将导丝描述成被一条被坐标化的曲线Г={t,m1,m2},其中t表示沿着曲线的单位切线,m1,m2表示曲线的两条单位法向量,该三个正交向量构成一个局部移动坐标系{t(s),m1(s),m2(s)},其中s表示曲线的弧长。利用连续的基尔霍夫弹性杆模型计算被坐标化的导丝的连续弹性能量E(Г),其中E(Г)=E弯曲(Г)+E扭曲(Г)。然后再利用离散的基尔霍夫弹性杆模型对连续的弹性能量进行离散化处理,推导出导丝的离散弹性能量的公式。最后利用拉格朗日方程对导丝的离散弹性能量进行导丝位置的求导,推导出导丝的形变受力公式。
所述的连续弹性能量E(Г)=E弯曲(Г)+E扭曲(Г)的具体公式为:
其中二维向量ω=(ω1,ω2)T表示导丝在局部移动坐标系中的曲率向量,m=m′1m2表示导丝在局部移动坐标系中的扭曲率;其中α表示导丝的弯曲模数,β表示导丝的扭曲系数。ω1,ω2分别表示导丝在局部移动坐标系m1,m2上的曲率,ds表示对导丝的弧长进行积分。
所述的离散弹性能量E(Г)=E弯曲(Г)+E弯曲(Г)的具体公式为:
其中i=1,2,3,…n表示将导丝进行离散分段后的段号,表示第i段导丝的长度,
θi表不导丝的第i顶点的旋转角度。
第三步,进行导丝的渲染,相邻的两导丝质点之间绘制一个圆柱体,同时在相邻的两圆柱体之间绘制一个小圆球进行连接。
第四步,采用K-DOPS树进行导丝与血管之间的碰撞检测。碰撞检测中的基元是三角面片与直线段(三角面片用于模拟血管,直线段用于模拟导丝):对三角面片和直线段,在进行碰撞检测和响应时,视为三角棱柱和两头有相同半径的半球状的圆柱体,当线段进入圆柱体,即判定二者相互碰撞。
第五步,采用非迭代的约束方法计算导丝与血管之间的接触力并进行碰撞响应。所述的采用非迭代的约束方法的具体实现过程为:首先计算导丝在血管中的穿刺向量d,然后根据穿刺向量计算接触力最后利用Verlet积分公式来进行导丝的位置更新,确保导丝仅在血管内部运动,而不穿透血管壁。
第六步,使用Verlet积分公式来迭代更新导丝的位置与速度。本文采用的是速度Verlet积分公式,具体的速度Verlet积分公式如下:
第七步,利用拉格朗日乘数和一种快速投影方法来实现导丝的不可伸缩特性。所述的快速投影方法的具体实现过程为:
701)构造导丝的长度约束函数C(p1,p2)=|p1,p2|-d,其中p1,p2为导丝中相邻两质点的位置坐标,d为该相邻两质点的原始长度。
702)如果长度约束函数的值超过了最低阈值,则对质点p1,p2进行位置修正并且循环执行该过程直至长度约束函数的值小于最低阈值,其中修正公式为: 因此p1,p2的修正量ΔP1,ΔP2为:
第八步,力反馈渲染,使用力反馈设备Phantom Omni渲染导丝给用户的作用力。所述的力反馈渲染的具体实现方法为:导丝的仿真频率只有50Hz,而力反馈设备的更新频率只要要达到1000Hz才能提供平滑的力触觉,为了实现两者的同步,我们采用线性插值的方法来计算力反馈设备的力输入参数。即对***计算的相邻时间间隔(20ms)的力f1,f2进行线性插值(即20等分),则力反馈设备的力输入参数为fi=f1+(f2-f1)i/20。
Claims (5)
1.一种微创血管介入手术中导丝的实时运动仿真方法,其特征在于,该方法包括以下步骤:
1)由于血管在仿真的过程中不产生形变,采用三角形表面网格对血管进行建模;
2)采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力;
3)进行导丝的渲染,相邻的两导丝质点之间绘制一个圆柱体,同时在相邻的两圆柱体之间绘制一个圆球进行连接;
4)采用K-DOPS树进行导丝与血管之间的碰撞检测;
5)采用非迭代的约束方法计算导丝与血管之间的接触力并进行碰撞响应;
6)使用Verlet积分公式来迭代更新导丝的位置与速度;
7)利用拉格朗日乘数和快速投影方法来实现导丝的不可伸缩特性;
8)使用力反馈设备Phantom Omni来渲染导丝给用户的作用力;
采用非迭代的约束方法的具体实现过程为:
首先计算导丝在血管中的穿刺向量d;
然后根据穿刺向量计算接触力M为导丝的质量。
2.根据权利要求1所述的一种微创血管介入手术中导丝的实时运动仿真方法,其特征在于,所述的采用基尔霍夫弹性杆模型,对导丝进行建模,并利用拉格朗日方程计算导丝在仿真过程产生的形变力的具体实现过程为:
21)将导丝描述成一条被坐标化的曲线Γ={t,m1,m2},其中t表示沿着曲线的单位切线,m1,m2表示曲线的两条单位法向量,该三个正交向量构成一个局部移动坐标系{t(s),m1(s),m2(s)},其中s表示曲线的弧长;
22)利用连续的基尔霍夫弹性杆模型计算被坐标化的导丝的连续弹性能量E(Γ),其中E(Γ)=E弯曲(Γ)+E扭曲(Γ);
23)利用离散的基尔霍夫弹性杆模型对连续的弹性能量进行离散化处理,推导出导丝的离散弹性能量的公式;
24)利用拉格朗日方程对导丝的离散弹性能量进行导丝位置的求导,计算出导丝的形变力。
3.根据权利要求2所述的一种微创血管介入手术中导丝的实时运动仿真方法,其特征在于,所述的连续弹性能量E(Γ)=E弯曲(Γ)+E扭曲(Γ)的具体公式为:
其中二维向量ω=(ω1,ω2)T表示导丝在局部移动坐标系中的曲率向量,m=m′1m2表示导丝在局部移动坐标系中的扭曲率;其中α表示导丝的弯曲模数,β表示导丝的扭曲系数;ω1,ω2分别表示导丝在局部移动坐标系m1,m2上的曲率,ds表示对导丝的弧长进行积分。
4.根据权利要求3所述的一种微创血管介入手术中导丝的实时运动仿真方法,其特征在于,所述的离散弹性能量E(Γ)/=E弯曲(Γ)/+E扭曲(Γ)/的具体公式为:
其中i=1,2,3,…n表示将导丝进行离散分段后的段号,表示第i段导丝的长度,θi表示导丝的第i个顶点的旋转角度,ωi表示第i段导丝在局部移动坐标系中的曲率向量,mi表示第i段导丝在局部移动坐标系中的扭曲率。
5.根据权利要求1所述的一种微创血管介入手术中导丝的实时运动仿真方法,其特征在于,所述的力反馈渲染的具体实现方法为:
导丝的仿真频率为50Hz,而力反馈设备的更新频率要达到1000Hz才能提供平滑的力触觉,为了实现两者的同步,采用线性插值的方法来计算力反馈设备的力输入参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210413201.5A CN102968811B (zh) | 2012-10-25 | 2012-10-25 | 微创血管介入手术中导丝的实时运动仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210413201.5A CN102968811B (zh) | 2012-10-25 | 2012-10-25 | 微创血管介入手术中导丝的实时运动仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102968811A CN102968811A (zh) | 2013-03-13 |
CN102968811B true CN102968811B (zh) | 2015-03-04 |
Family
ID=47798936
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210413201.5A Active CN102968811B (zh) | 2012-10-25 | 2012-10-25 | 微创血管介入手术中导丝的实时运动仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102968811B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103280145B (zh) * | 2013-05-03 | 2016-01-13 | 上海交通大学 | 心血管介入虚拟手术力反馈*** |
CN103699776A (zh) * | 2013-12-01 | 2014-04-02 | 北京航空航天大学 | 一种面向心血管介入手术仿真的导丝模拟方法 |
CN104318840A (zh) * | 2014-10-22 | 2015-01-28 | 北京航空航天大学 | 一种基于弹簧质子模型的医学手术器械导丝的模拟方法 |
CN104463965A (zh) * | 2014-12-17 | 2015-03-25 | 中国科学院自动化研究所 | 一种面向微创心血管介入手术的训练场景模拟***及方法 |
CN105261069B (zh) * | 2015-11-05 | 2018-04-06 | 北京航空航天大学 | 基于gpu的软组织器官元球模型自动生成与碰撞检测的方法 |
CN105448171B (zh) * | 2015-11-06 | 2018-04-06 | 北京航空航天大学 | 一种基于一体化链甲模型的医学手术器械导丝的模拟方法 |
CN106067269B (zh) * | 2016-05-13 | 2019-03-15 | 中国科学院自动化研究所 | 虚拟心血管介入手术培训***中反馈力的确定方法及*** |
CN106709902B (zh) * | 2016-11-14 | 2020-06-19 | 上海交通大学 | 微创血管介入手术中导丝受血流作用的实时运动仿真方法 |
CN110729051B (zh) * | 2019-10-10 | 2022-11-22 | 中国科学院深圳先进技术研究院 | 一种介入手术中的导丝力学分析方法、***及电子设备 |
CN111243746B (zh) * | 2020-01-20 | 2023-05-26 | 上海奥朋医疗科技有限公司 | 一种血管介入手术机器人的手术仿真方法及*** |
CN111797506B (zh) * | 2020-06-10 | 2022-07-29 | 上海工程技术大学 | 一种主从导丝控制方法 |
CN112017516B (zh) * | 2020-08-26 | 2022-10-04 | 深圳市爱博医疗机器人有限公司 | 一种远程血管介入手术训练*** |
CN113823161B (zh) * | 2021-09-22 | 2023-02-28 | 山东静禾医疗科技有限公司 | 双通道式介入手术模拟装置及模拟控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3084098B2 (ja) * | 1991-08-30 | 2000-09-04 | 株式会社日立製作所 | コンピュータアニメーションの対話処理方式及び装置 |
CN101295409A (zh) * | 2008-06-05 | 2008-10-29 | 上海交通大学 | 虚拟手术***中形变物体的实时模拟*** |
CN101919697A (zh) * | 2010-09-07 | 2010-12-22 | 上海交通大学 | 介入导丝二自由度运动的非接触式检测方法 |
-
2012
- 2012-10-25 CN CN201210413201.5A patent/CN102968811B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3084098B2 (ja) * | 1991-08-30 | 2000-09-04 | 株式会社日立製作所 | コンピュータアニメーションの対話処理方式及び装置 |
CN101295409A (zh) * | 2008-06-05 | 2008-10-29 | 上海交通大学 | 虚拟手术***中形变物体的实时模拟*** |
CN101919697A (zh) * | 2010-09-07 | 2010-12-22 | 上海交通大学 | 介入导丝二自由度运动的非接触式检测方法 |
Non-Patent Citations (4)
Title |
---|
Analytical guide wire motion algorithm for simulation of endovascular interventions;Dr M.K.Konings等;《Medical and Biological Engineering and Computing》;20031101;第41卷(第6期);第689-700页 * |
Real-time guidewire simulation in complex vascular models;Vincent Luboz等;《The Visual Computer》;20090901;第25卷(第9期);第827-834页 * |
Virtual reality based three-dimensional guide wire propagation simulation for endovascular intervention;周正东等;《Transactions of Nanjing University of Aeronautics & Astronautics》;20100315;第27卷(第1期);第62-69页 * |
周正东等.三维血管介入手术模拟方法.《东南大学学报(自然科学版)》.2010,第40卷(第2期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102968811A (zh) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102968811B (zh) | 微创血管介入手术中导丝的实时运动仿真方法 | |
CN107330972B (zh) | 模拟生物力学特性的实时软组织形变方法和*** | |
CN103400023B (zh) | 软组织形变仿真方法 | |
CN106504329B (zh) | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 | |
CN109308739A (zh) | 一种基于虚拟现实的软组织微创手术训练方法 | |
CN104318056A (zh) | 基于位置动力学的软组织变形和切割模拟方法 | |
CN103699776A (zh) | 一种面向心血管介入手术仿真的导丝模拟方法 | |
CN1975784A (zh) | 基于骨骼线体网格的质点弹簧形变仿真方法 | |
CN105261069A (zh) | 基于gpu的软组织器官元球模型自动生成与碰撞检测的方法 | |
CN103729555A (zh) | 一种模拟血流与血管壁作用的方法和装置 | |
Zhang et al. | A fast and accurate vascular tissue simulation model based on point primitive method | |
Zhang et al. | A Real-time Cutting Model Based on Finite Element and Order Reduction. | |
CN103559741A (zh) | 虚拟手术中基于粒子的多相耦合方法 | |
Ishikawa et al. | Sonic-boom prediction using Euler CFD codes with structured/unstructured overset method | |
CN105955458A (zh) | 反馈力生成方法及装置 | |
CN102867336B (zh) | 一种基于热力学模型的固体燃烧过程模拟方法 | |
Zhang et al. | Real-time simulation of deformable soft tissue based on mass-spring and medial representation | |
CN106934189A (zh) | 外科手术软组织变形的仿真方法及装置 | |
CN115908664A (zh) | 人机交互的动画生成方法、装置、计算机设备、存储介质 | |
CN104318840A (zh) | 一种基于弹簧质子模型的医学手术器械导丝的模拟方法 | |
CN102855353A (zh) | 一种埋藏式缝合技术的实时模拟方法 | |
Duysak et al. | Fast simulation of deformable objects | |
KR101110342B1 (ko) | 유체 시뮬레이션 형상 제어 시스템 및 방법 | |
Qiao et al. | The research of soft tissue deformation based on mass-spring model | |
CN103714205B (zh) | 一种模拟柔性体可旋转变形的发条弹簧模型 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |