CN103400023A - 软组织形变仿真方法 - Google Patents
软组织形变仿真方法 Download PDFInfo
- Publication number
- CN103400023A CN103400023A CN2013102675577A CN201310267557A CN103400023A CN 103400023 A CN103400023 A CN 103400023A CN 2013102675577 A CN2013102675577 A CN 2013102675577A CN 201310267557 A CN201310267557 A CN 201310267557A CN 103400023 A CN103400023 A CN 103400023A
- Authority
- CN
- China
- Prior art keywords
- particle
- soft tissue
- force
- deformation
- model
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种软组织形变仿真方法,包括以下步骤:建立软组织的生物力学模型、并对其中的各个质点进行初始化;力反馈设备对软组织施加作用力,进行碰撞检测;用改进的欧拉算法计算运动状态信息;将模型每个时间步长的状态输出到显示屏上,动态显示软组织形变过程;计算反馈力并输出触觉反馈。通过采用上述步骤,能够有效解决虚拟手术仿真中的实时性、精确性以及反馈力光滑性问题,提高了软组织形变仿真的精度和实时性,从而满足虚拟手术仿真的需要。
Description
技术领域
本发明涉及一种软组织形变仿真方法。
背景技术
虚拟手术仿真是虚拟现实技术的重要应用之一,其目标是利用各种医学影像数据以及虚拟现实技术在计算机中建立逼真的虚拟手术环境,使得医生或者受训者可以借助该虚拟环境进行手术训练。
在虚拟手术仿真中,软组织器官的形变模型决定了虚拟手术仿真的视觉效果和力反馈精度。常见的软组织形变计算模型分为两大类:基于几何的形变模型和基于物理的形变模型。基于几何的形变模型仅仅考虑了几何形态的变化,而忽略了软组织的实际力学本构方程以及形变过程中物体质量、力或其他物理现象的作用,因此不能真实的反映软组织的形变过程,该模型目前已较少使用;基于物理的形变模型则基于软组织的力学本构方程,通过相应的计算模型得出组织受力时的形变,能够更加真实的反映组织的形变,因此目前使用较多。目前,基于物理的形变模型主要有有限元模型和质点-弹簧模型。有限元模型的优点是参数连续化,在大形变情况下也能模拟组织形变,具有较高的精度和适应性,但其建模复杂,计算量大,需要大量的预处理工作,尤其在切割过程中难以实现快速的网格重构,因此,实时性较差。与有限元模型相比,质点-弹簧模型不需要将参数连续化,容易实现,运算速度较快。而质点-弹簧模型一般分为面模型和体模型两种,面模型的网格密度可以根据仿真的要求进行划分,计算效率高,可以满足仿真的实时性,但是,只能反映组织表面的形变特征,不能刻画三维实体的内部结构,是以牺牲精度为代价来保证实时性,体模型是对具有一定体积的实体组织进行建模,能够同时表现表面和内部结构,具有较高的精度,但实时性较差。
对于上述模型,目前求解其二阶微分方程的主要方法有欧拉法和龙格-库塔法。显式欧拉法计算量小,但是一般只一阶收敛,精度不高,为实现逼真的模拟效果,时间步长需要设置很小,从而导致整个形变过程延长,并且由于存在进退性冲击波,不适合用于质点-弹簧模型。隐式欧拉法精度高,稳定性好,但是计算量却较大。龙格-库塔法也具有较高的精度和稳定性,但是,在一步计算中需要计算四次函数的值,计算量大,实时性较差。在虚拟手术仿真中,光滑的触觉反馈对实时性的要求更高。因此,如何设计出一种支持光滑触觉反馈的快速且逼真的软组织形变仿真方法,尽可能的满足虚拟手术对实时性以及精确性的要求,已成为虚拟手术面临的首要问题。
发明内容
为克服以上现有技术的不足,本发明要解决的技术问题是提供一种软组织形变仿真方法,能够有效构建软组织生物力学模型,并对其中质点的速度、位移以及受力信息进行高效求解,从而有效解决虚拟手术仿真中的实时性、精确性以及反馈力光滑性问题,进而满足虚拟手术仿真的需要。
本发明的技术方案是:
一种软组织形变仿真方法,包含以下步骤:
步骤1):采集软组织的数据信息,并采用基于四面体的质点-弹簧体模型建立软组织的生物力学模型,该生物力学模型由n个质点组成,并且对于其中任意一个质点i,满足以下方程:
其中,mi为质点i的质量,xi表示质点i的位置矢量,表示与质点i相连的所有质点j(包括j1至jq)施加于质点i的内力,xj表示质点j的位置矢量,表示质点i所受外力,|xj-xi|表示两个质点间位置矢量之差,|xj-xi|0是发生形变前弹簧的长度,vi为质点i的速度,vj为质点j的速度,ai为质点i的加速度,μ为弹簧的弹性系数,η为阻尼器的阻尼系数;
步骤2):依据步骤1)所建的生物力学模型,对其中的各个质点进行初始化,构建模型的初始状态,并且计算该模型中的每个质点与相连质点间的弹簧初始长度;
步骤3):外接的力反馈设备对软组织施加作用力,进行碰撞检测,确定软组织上发生碰撞的质点和受力发生形变的区域,以及软组织被按压或者拉伸的长度;
步骤4):用改进的欧拉算法计算发生形变区域内的各个质点的运动状态信息;
步骤5):依据步骤3)和步骤4)计算反馈力,并将该反馈力输出至力反馈设备;
步骤6):循环执行步骤4)至步骤5),计算每个时间步长中各质点的运动状态,并在显示器上动态显示软组织的形变过程。
上述软组织形变仿真方法,其中步骤2)中的初始化包括初始化各个质点的位置、质量、速度、加速度和受力信息。
上述软组织形变仿真方法,其中步骤4)中的运动状态信息包括发生形变区域内的各个质点的位置、速度和受力随时间变化的信息。
上述软组织形变仿真方法,其中步骤4)的计算过程具体为:
步骤4.2):计算质点r的加速度:
步骤4.3):用改进的欧拉算法求解质点r的速度和位移,其中用显式欧拉法对速度vr进行迭代求解,用隐式欧拉法对位置矢量xr进行求解:
步骤4.5):循环执行步骤4.1)至4.4),计算受力区域内其他质点的位置、速度及受力信息。
本发明的有益效果是:本发明通过有效构建软组织生物力学模型,利用改进的欧拉算法求解软组织形变过程,对其中质点的速度、位移以及受力信息进行高效求解,并实现了光滑反馈力的输出,从而有效解决虚拟手术仿真中的实时性、精确性以及反馈力光滑性问题,提高了软组织形变仿真的精度和实时性,进而满足了虚拟手术仿真的需要。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1是软组织形变仿真方法的总体流程图示意图;
图2是软组织形变仿真方法中的质点-弹簧体模型的拓扑结构示意图。
具体实施方式
如图1和图2所示,一种软组织形变仿真方法,包含以下步骤:
步骤1):采集软组织的数据信息,并采用基于四面体的质点-弹簧体模型建立软组织的生物力学模型,该生物力学模型由n个质点组成,并且对于其中任意一个质点i,满足以下该动力学方程:
本发明采用Kelvin黏弹性模型来反映软组织的力学特性,即质点与质点间由弹簧和阻尼器并联连接。依据Kelvin模型,当力F作用到弹簧上时,弹簧服从胡克定律,伸长量u与F成正比,即F=μu;当力F作用到阻尼器上时,阻尼器服从牛顿黏性定律,形变速度u′与力F成正比,即F=μu′。因此,当软组织受到外力作用发生形变时,两个相连质点之间的内力由弹簧产生的弹性力和阻尼器产生的粘滞力组成,即fin=μu+ηu′。因此,施加于某质点的内力的计算方法为:
并且,质点的加速度计算方法为:
上述公式中,mi为质点i的质量,xi表示质点i的位置矢量,表示与质点i相连的所有质点j(包括j1至jq)施加于质点i的内力,xj表示质点j的位置矢量,表示质点i所受外力,|xj-xi|表示两个质点间位置矢量之差,|xj-xi|0是发生形变前弹簧的长度,vi为质点i的速度,vj为质点j的速度,ai为质点i的加速度,μ为弹簧的弹性系数,η为阻尼器的阻尼系数。
步骤2):依据步骤1)所建的生物力学模型,对其中的各个质点进行初始化,可包括初始化各个质点的位置、质量、速度、加速度和受力信息,并构建模型的初始状态,计算该模型中的每个质点与相连质点间的弹簧初始长度。上述初始化信息均存入邻接链表中。
步骤3):外接的力反馈设备对软组织施加作用力,进行碰撞检测,确定软组织上发生碰撞的质点和受力发生形变的区域,以及软组织被按压或者拉伸的长度。
步骤4):用改进的欧拉算法计算发生形变区域内的各个质点的运动状态信息,其中运动状态信息可包括发生形变区域内的各个质点的位置、速度和受力随时间变化的信息。
步骤5):依据步骤3),当外接的力反馈设备对软组织施加作用力时,计算设备按压或者拉伸软组织的长度,并结合依据步骤4)得出的各个质点的运动状态信息,计算反馈力,并将该反馈力输出至力反馈设备。此时,用户可以通过操纵杆感受到力的反馈;
步骤6):循环执行步骤4)至步骤5),计算每个时间步长中各质点的运动状态,利用OpenGL技术,将每个时间步长内计算出的质点按照位移向量输出至显示器,即可得到软组织模型在外力作用下的形变过程。
对于上述软组织形变仿真方法,其中步骤4)的计算过程可具体包括为:
步骤4.2):计算质点r的加速度:
步骤4.3):用改进的欧拉算法求解质点r的速度和位移,其中用显式欧拉法对速度vr进行迭代求解,用隐式欧拉法对位置矢量xr进行求解:
步骤4.5):循环执行步骤4.1)至4.4),计算受力区域内其他质点的位置、速度及受力信息。
在上述步骤和方法中,通过有效构建软组织生物力学模型,利用改进的欧拉算法求解软组织形变过程,对其中质点的速度、位移以及受力信息进行高效求解,并实现了光滑反馈力的输出,从而有效解决虚拟手术仿真中的实时性、精确性以及反馈力光滑性问题,提高了软组织形变仿真的精度和实时性,进而满足了虚拟手术仿真的需要。
上面结合附图对本发明优选的具体实施方式和实施例作了详细说明,但是本发明并不限于上述实施方式和实施例,在本领域技术人员所具备的知识范围内,还可以在不脱离本发明构思的前提下作出各种变化。
Claims (4)
1.一种软组织形变仿真方法,其特征在于包含以下步骤:
步骤1):采集软组织的数据信息,并采用基于四面体的质点-弹簧体模型建立软组织的生物力学模型,该生物力学模型由n个质点组成,并且对于其中任意一个质点i,满足以下方程:
其中,mi为质点i的质量,xi表示质点i的位置矢量,表示与质点i相连的所有质点j(包括j1至jq)施加于质点i的内力,xj表示质点j的位置矢量,表示质点i所受外力,|xj-xi|表示两个质点间位置矢量之差,|xj-xi|0是发生形变前弹簧的长度,vi为质点i的速度,vj为质点j的速度,ai为质点i的加速度,μ为弹簧的弹性系数,η为阻尼器的阻尼系数;
步骤2):依据步骤1)所建的生物力学模型,对其中的各个质点进行初始化,构建模型的初始状态,并且计算该模型中的每个质点与相连质点间的弹簧初始长度;
步骤3):外接的力反馈设备对软组织施加作用力,进行碰撞检测,确定软组织上发生碰撞的质点和受力发生形变的区域,以及软组织被按压或者拉伸的长度;
步骤4):用改进的欧拉算法计算发生形变区域内的各个质点的运动状态信息;
步骤5):依据步骤3)和步骤4)计算反馈力,并将该反馈力输出至力反馈设备;
步骤6):循环执行步骤4)至步骤5),计算每个时间步长中各质点的运动状态,并在显示器上动态显示软组织的形变过程。
2.如权利要求1所述的软组织形变仿真方法,其特征在于:所述步骤2)中的初始化包括初始化各个质点的位置、质量、速度、加速度和受力信息。
3.如权利要求2所述的软组织形变仿真方法,其特征在于:所述步骤4)中的运动状态信息包括发生形变区域内的各个质点的位置、速度和受力随时间变化的信息。
4.如权利要求3所述的软组织形变仿真方法,其特征在于:所述步骤4)中的计算过程具体为:
步骤4.2):计算质点r的加速度:
步骤4.3):用改进的欧拉算法求解质点r的速度和位移,其中用显式欧拉法对速度vr进行迭代求解,用隐式欧拉法对位置矢量xr进行求解:
步骤4.4):计算质点r开始运动后与其相连质点间弹簧的长度及弹簧的形变量,并计算质点r所受内力
步骤4.5):循环执行步骤4.1)至4.4),计算受力区域内其他质点的位置、速度及受力信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310267557.7A CN103400023B (zh) | 2013-06-28 | 2013-06-28 | 软组织形变仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310267557.7A CN103400023B (zh) | 2013-06-28 | 2013-06-28 | 软组织形变仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103400023A true CN103400023A (zh) | 2013-11-20 |
CN103400023B CN103400023B (zh) | 2016-11-02 |
Family
ID=49563651
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310267557.7A Expired - Fee Related CN103400023B (zh) | 2013-06-28 | 2013-06-28 | 软组织形变仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103400023B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105528491A (zh) * | 2015-12-24 | 2016-04-27 | 北京农业信息技术研究中心 | 一种植物动态虚拟仿真方法及系统 |
CN105912859A (zh) * | 2016-04-11 | 2016-08-31 | 浙江工业大学义乌科学技术研究院有限公司 | 一种基于质点弹簧和流体力学的组织形变方法 |
CN106485028A (zh) * | 2016-11-03 | 2017-03-08 | 南京信息工程大学 | 支持实时按压的柔性体变形仿真建模方法 |
CN106570341A (zh) * | 2016-11-14 | 2017-04-19 | 南昌大学 | 一种基于无网格径向基数据拟合的软组织受力形变模型建模方法 |
CN106650251A (zh) * | 2016-12-14 | 2017-05-10 | 南京信息工程大学 | 一种针灸力反馈形变模型的建模方法 |
CN106682425A (zh) * | 2016-12-29 | 2017-05-17 | 天津瀚海星云数字科技有限公司 | 带阻尼的软体受力变形的模拟仿真方法 |
CN106781941A (zh) * | 2016-11-24 | 2017-05-31 | 北京理工大学 | 一种用于模拟微创穿刺手术的方法及其系统 |
CN106777658A (zh) * | 2016-12-12 | 2017-05-31 | 南京信息工程大学 | 基于放松lra约束算法的缝合模拟方法 |
CN108511074A (zh) * | 2018-03-26 | 2018-09-07 | 福建师范大学福清分校 | 一种基于空间核映射和子空间聚集的软组织形变方法 |
CN109036567A (zh) * | 2018-01-10 | 2018-12-18 | 福建江夏学院 | 一种基于子空间凝聚算法的软组织形变仿真方法 |
CN109492250A (zh) * | 2018-09-29 | 2019-03-19 | 南京信息工程大学 | 一种非线性各向异性质量弹簧模型的软组织变形模拟方法 |
CN109658796A (zh) * | 2018-11-14 | 2019-04-19 | 华中科技大学 | 一种软组织损伤及异常连接模拟试验方法 |
CN109993830A (zh) * | 2019-04-10 | 2019-07-09 | 福建师范大学福清分校 | 一种软组织表面刺穿形变仿真方法及终端 |
CN110289104A (zh) * | 2019-06-27 | 2019-09-27 | 南京信息工程大学 | 软组织按压和形变恢复的模拟方法 |
CN110931121A (zh) * | 2019-11-29 | 2020-03-27 | 重庆邮电大学 | 一种基于Hololens的远程手术指导装置及操作方法 |
CN111026269A (zh) * | 2019-12-04 | 2020-04-17 | 上海褚信医学科技有限公司 | 基于力反馈的生物组织结构的触觉反馈方法、装置和设备 |
CN114117835A (zh) * | 2021-07-28 | 2022-03-01 | 上海华模科技有限公司 | 一种眼科手术模拟训练器中的眼球建模方法 |
CN117950504A (zh) * | 2024-03-26 | 2024-04-30 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
CN117950504B (zh) * | 2024-03-26 | 2024-05-31 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07100126A (ja) * | 1993-10-01 | 1995-04-18 | Shiseido Co Ltd | 表情しわの模擬方法および装置 |
US20020183992A1 (en) * | 2001-02-08 | 2002-12-05 | Inria Institut National De Recherche En Informatique Et En Automatique | Apparatus for simulating the deformation of materials, notably of soft body tissues |
CN1975784A (zh) * | 2006-12-28 | 2007-06-06 | 上海交通大学 | 基于骨骼线体网格的质点弹簧形变仿真方法 |
CN102044086A (zh) * | 2010-11-30 | 2011-05-04 | 华北水利水电学院 | 一种软组织形变仿真方法 |
CN102262699A (zh) * | 2011-07-27 | 2011-11-30 | 华北水利水电学院 | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 |
CA2713007A1 (en) * | 2010-08-11 | 2012-02-11 | National Research Council Of Canada | Presentation-enhanced solid mechanical simulation |
-
2013
- 2013-06-28 CN CN201310267557.7A patent/CN103400023B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07100126A (ja) * | 1993-10-01 | 1995-04-18 | Shiseido Co Ltd | 表情しわの模擬方法および装置 |
US20020183992A1 (en) * | 2001-02-08 | 2002-12-05 | Inria Institut National De Recherche En Informatique Et En Automatique | Apparatus for simulating the deformation of materials, notably of soft body tissues |
CN1975784A (zh) * | 2006-12-28 | 2007-06-06 | 上海交通大学 | 基于骨骼线体网格的质点弹簧形变仿真方法 |
CA2713007A1 (en) * | 2010-08-11 | 2012-02-11 | National Research Council Of Canada | Presentation-enhanced solid mechanical simulation |
CN102044086A (zh) * | 2010-11-30 | 2011-05-04 | 华北水利水电学院 | 一种软组织形变仿真方法 |
CN102262699A (zh) * | 2011-07-27 | 2011-11-30 | 华北水利水电学院 | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 |
Non-Patent Citations (2)
Title |
---|
叶秀芬,乔冰,郭书祥,郭庆昌: "虚拟手术仿真中人体软组织形变技术的研究", 《计算机应用》, vol. 29, no. 2, 28 February 2009 (2009-02-28) * |
闫雒恒: "基于改进弹簧振子模型的软组织形变仿真", 《计算机仿真》, vol. 29, no. 7, 31 July 2012 (2012-07-31), pages 332 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105528491A (zh) * | 2015-12-24 | 2016-04-27 | 北京农业信息技术研究中心 | 一种植物动态虚拟仿真方法及系统 |
CN105912859A (zh) * | 2016-04-11 | 2016-08-31 | 浙江工业大学义乌科学技术研究院有限公司 | 一种基于质点弹簧和流体力学的组织形变方法 |
CN105912859B (zh) * | 2016-04-11 | 2018-07-17 | 浙江工业大学义乌科学技术研究院有限公司 | 一种基于质点弹簧和流体力学的组织形变模拟方法 |
CN106485028A (zh) * | 2016-11-03 | 2017-03-08 | 南京信息工程大学 | 支持实时按压的柔性体变形仿真建模方法 |
CN106485028B (zh) * | 2016-11-03 | 2019-11-22 | 南京信息工程大学 | 支持实时按压的柔性体变形仿真建模方法 |
CN106570341B (zh) * | 2016-11-14 | 2019-06-18 | 南昌大学 | 一种基于无网格径向基数据拟合的软组织受力形变模型建模方法 |
CN106570341A (zh) * | 2016-11-14 | 2017-04-19 | 南昌大学 | 一种基于无网格径向基数据拟合的软组织受力形变模型建模方法 |
CN106781941A (zh) * | 2016-11-24 | 2017-05-31 | 北京理工大学 | 一种用于模拟微创穿刺手术的方法及其系统 |
CN106777658A (zh) * | 2016-12-12 | 2017-05-31 | 南京信息工程大学 | 基于放松lra约束算法的缝合模拟方法 |
CN106777658B (zh) * | 2016-12-12 | 2020-02-18 | 南京信息工程大学 | 基于放松lra约束算法的缝合模拟方法 |
CN106650251A (zh) * | 2016-12-14 | 2017-05-10 | 南京信息工程大学 | 一种针灸力反馈形变模型的建模方法 |
CN106682425A (zh) * | 2016-12-29 | 2017-05-17 | 天津瀚海星云数字科技有限公司 | 带阻尼的软体受力变形的模拟仿真方法 |
CN109036567B (zh) * | 2018-01-10 | 2021-11-09 | 福建江夏学院 | 一种基于子空间凝聚算法的软组织形变仿真方法 |
CN109036567A (zh) * | 2018-01-10 | 2018-12-18 | 福建江夏学院 | 一种基于子空间凝聚算法的软组织形变仿真方法 |
CN108511074A (zh) * | 2018-03-26 | 2018-09-07 | 福建师范大学福清分校 | 一种基于空间核映射和子空间聚集的软组织形变方法 |
CN108511074B (zh) * | 2018-03-26 | 2021-11-09 | 福建师范大学福清分校 | 一种基于空间核映射和子空间聚集的软组织形变方法 |
CN109492250A (zh) * | 2018-09-29 | 2019-03-19 | 南京信息工程大学 | 一种非线性各向异性质量弹簧模型的软组织变形模拟方法 |
CN109658796A (zh) * | 2018-11-14 | 2019-04-19 | 华中科技大学 | 一种软组织损伤及异常连接模拟试验方法 |
CN109658796B (zh) * | 2018-11-14 | 2021-04-20 | 华中科技大学 | 一种软组织损伤及异常连接模拟试验方法 |
CN109993830B (zh) * | 2019-04-10 | 2022-05-03 | 福建师范大学福清分校 | 一种软组织表面刺穿形变仿真方法及终端 |
CN109993830A (zh) * | 2019-04-10 | 2019-07-09 | 福建师范大学福清分校 | 一种软组织表面刺穿形变仿真方法及终端 |
CN110289104A (zh) * | 2019-06-27 | 2019-09-27 | 南京信息工程大学 | 软组织按压和形变恢复的模拟方法 |
CN110289104B (zh) * | 2019-06-27 | 2022-12-23 | 南京信息工程大学 | 软组织按压和形变恢复的模拟方法 |
CN110931121A (zh) * | 2019-11-29 | 2020-03-27 | 重庆邮电大学 | 一种基于Hololens的远程手术指导装置及操作方法 |
CN111026269A (zh) * | 2019-12-04 | 2020-04-17 | 上海褚信医学科技有限公司 | 基于力反馈的生物组织结构的触觉反馈方法、装置和设备 |
CN111026269B (zh) * | 2019-12-04 | 2023-10-10 | 上海褚信医学科技有限公司 | 基于力反馈的生物组织结构的触觉反馈方法、装置和设备 |
CN114117835A (zh) * | 2021-07-28 | 2022-03-01 | 上海华模科技有限公司 | 一种眼科手术模拟训练器中的眼球建模方法 |
CN117950504A (zh) * | 2024-03-26 | 2024-04-30 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
CN117950504B (zh) * | 2024-03-26 | 2024-05-31 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103400023B (zh) | 2016-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103400023A (zh) | 软组织形变仿真方法 | |
CN102262699B (zh) | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 | |
Teran et al. | Finite volume methods for the simulation of skeletal muscle | |
CN102968811B (zh) | 微创血管介入手术中导丝的实时运动仿真方法 | |
CN103559741B (zh) | 虚拟手术中基于粒子的多相耦合方法 | |
CN102663197A (zh) | 一种基于运动捕获的虚拟手抓取的仿真方法 | |
CN107330972B (zh) | 模拟生物力学特性的实时软组织形变方法和系统 | |
CN101496028A (zh) | 使用几何推动式模型模拟可变形物体的方法 | |
CN110289104B (zh) | 软组织按压和形变恢复的模拟方法 | |
CN103699776A (zh) | 一种面向心血管介入手术仿真的导丝模拟方法 | |
CN103729555A (zh) | 一种模拟血流与血管壁作用的方法和装置 | |
WO2018213607A1 (en) | Computerized rendering of objects having anisotropic elastoplasticity for codimensional frictional contact | |
CN103389649B (zh) | 一种基于球面拼接网格的飞行器机动运动模拟方法 | |
CN103914872A (zh) | 一种基于简化模态分析法的树动画模拟方法 | |
Lee et al. | Volumetric Object Modeling Using Internal Shape Preserving Constraint in Unity 3D. | |
Ishikawa et al. | Sonic-boom prediction using Euler CFD codes with structured/unstructured overset method | |
CN104112375A (zh) | 一种虚拟微创手术训练系统 | |
CN102521856B (zh) | 电子毛笔书写笔形变形方法 | |
CN107146288A (zh) | 虚拟手术中支持实时按压形变的软组织模型建模方法 | |
CN105912859B (zh) | 一种基于质点弹簧和流体力学的组织形变模拟方法 | |
CN115310339A (zh) | 基于物质点法的具有表面张力效应的固液耦合模拟方法 | |
Yu et al. | Research of simulation in character animation based on physics engine | |
CN102810212A (zh) | 基于cuda架构的针织布料仿真方法 | |
Wang et al. | An unfixed-elasticity mass spring model based simulation for soft tissue deformation | |
CN106934189A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161102 Termination date: 20190628 |
|
CF01 | Termination of patent right due to non-payment of annual fee |