CN101529487B - 柔性的对象模拟器 - Google Patents

柔性的对象模拟器 Download PDF

Info

Publication number
CN101529487B
CN101529487B CN200780038656.0A CN200780038656A CN101529487B CN 101529487 B CN101529487 B CN 101529487B CN 200780038656 A CN200780038656 A CN 200780038656A CN 101529487 B CN101529487 B CN 101529487B
Authority
CN
China
Prior art keywords
node
fragment
flexible elongated
nodes
elongated objects
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
Application number
CN200780038656.0A
Other languages
English (en)
Other versions
CN101529487A (zh
Inventor
M·伊基茨
D·尼尔森
C·J·乌尔里希
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Immersion Corp
Original Assignee
Immersion Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Immersion Corp filed Critical Immersion Corp
Publication of CN101529487A publication Critical patent/CN101529487A/zh
Application granted granted Critical
Publication of CN101529487B publication Critical patent/CN101529487B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B23/00—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
    • G09B23/28—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
    • G09B23/30—Anatomical models
    • G09B23/32—Anatomical models with moving parts

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Medicinal Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Medical Informatics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Educational Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Instructional Devices (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Media Introduction/Drainage Providing Device (AREA)

Abstract

一种模拟柔性的对象例如导管的方法包括将柔性的对象建模为多个片段，其中每个片段具有由一个或多个边连接的多个节点。该方法进一步包括：将每个片段建模为由节点和边形成的多个四面体元，以产生四面体有限元模型。该方法进一步包括：索引每个节点，以便第一节点的第一索引值和第二节点的第二索引值之间的最大差值是有界的并且约等于s＊k，其中s是在一个节点处连接的片段的最大个数，并且k是属于一个片段的节点的最大个数。

Description

柔性的对象模拟器
政府许可权
美国政府具有本发明的已付许可，并且根据合理的条款，如由国家科学基金会(National Science Foundation)授予的第DMI-0239344号基金所提供的条款，在特定情况下，美国政府具有要求专利所有人许可他人使用的权利。
技术领域
本发明的一个实施方式涉及柔性的对象的模拟器。更具体地，本发明的一个实施方式涉及用于模拟导管和其他线状结构的医疗装置模拟器。
背景技术
微创治疗的发展引起了对用于训练的基于计算机的模拟(computer-based simulation)的潜力的日益增加的兴趣。已设计出许多基于计算机的医疗模拟器来教授这些新技术。这些成果大多集中在腹腔镜手术，其中刚性设备被***到患者的腹腔，并且由产生高分辨率彩色图像的内窥相机来提供视觉反馈。这些模拟器是主要的技术训练工具(skills trainer)，其强调有助于腹腔镜手术的二维的手眼协调实践。
介入心脏病学与其他的微创治疗具有同样的使其适用于基于模拟器的学习的特征：其要求从二维显示和良好的手眼协调中综合理解三维的解剖构造。与手术一样，由不正确地操作心导管***引起的并发症可具有恶劣的后果，因此存在对于在尝试真正的过程之前进行基于模拟器的训练的强烈需要。
然而，介入心脏病学模拟具有其特有的挑战。首先，视觉反馈不是由可见光提供，而是由荧光检查提供，荧光检查必须被实时模拟，同时当荧光镜在病人周围移动时允许视点的变化。第二，导管、引导线和支架是柔性装置，并且因此必须将其建模为可变形的对象(deformable object)，而刚性的腹腔镜设备却不是这种情况。医师可以仅推、拉或扭转装置的近端来控制导管或引导线在血管网络内的移动。因为这种装置被限制在患者的脉管***内部，所以是输入力和接触力的组合允许装置朝目标移动。
模拟模型试图获得的线状结构(wire-like structure)或柔性的对象的主要特征包括几何非线性、高拉伸强度和对弯曲的低阻抗。然而，许多已知的柔性的对象模拟模型不适合实时应用
在医疗模拟环境中使用的用于柔性的对象的已知模型包括关节体方法(articulated body method)(“ABM”)，其将对象表示为由旋转和扭转弹簧连接的一组刚性片段(见，如Dawson等人的“Designing a Computer-Based Simulator for Interventional CardiologyTraining”，Catheterization and Cardiovascular Intervention51：522-527(2000))。然而，因为这些方法使用显式积分，所以其不能提供如许多应用所要求的必要的硬度和速度。此外，稳定性受到最小片段的长度的影响。
其他用于柔性的对象的已知模型使用梁有限元(见，如Cotin等人的“New Approaches to Catheter Navigation for InterventionalRadiology Simulation”，MICCAI(2005))。然而，由于使用显式积分和迭代解技术，这些方法呈现了与关节体方法类似的问题。
基于上述内容，存在对用于基于计算机的模拟的对柔性的对象进行建模和模拟柔性的对象的***和方法的需要。
发明内容
本发明的一个实施方式是模拟柔性的对象例如导管的方法。该方法包括：将柔性的对象建模为多个片段，其中每个片段具有由一个或多个边连接的多个节点。该方法进一步包括：将每个片段建模为由节点和边形成的多个四面体元，以产生四面体有限元模型。该方法进一步包括：索引每个节点，以便第一节点的第一索引值和第二节点的第二索引值之间的最大差值是有界的并且约等于s＊k，其中s是在一个节点处连接的片段的最大个数，并且k是属于一个片段的节点的最大个数。
附图说明
图1是根据本发明的一个实施方式用于模拟柔性的对象的***的透视图。
图2示出待被模拟的柔性的对象的一个片段。
图3示出如何通过***其他的边将每个棱柱进一步分解为三个四面体元。
图4示出根据本发明的一个实施方式的两个连续的片段和相应的节点，并且示出节点索引方案(node indexing scheme)。
图5是根据本发明的一个实施方式由图1的计算机执行来模拟柔性的对象例如导管的功能的流程图。
图6用图形示出根据本发明的一个实施方式的***矩阵A的带状结构(banded structure)。
图7示出根据本发明的一个实施方式的所得到的模拟的四面体模型。
图8示出根据本发明的一个实施方式的使用者所见的实际的模拟的导管。
具体实施方式
本发明的一个实施方式是用于柔性的对象(例如那些在医疗过程中使用或遇到的对象)的基于物理的建模和交互模拟的方法和***。使用带有线性弹性材料模型和共转弯曲(corotational warping)的有限元来对对象进行建模。
图1是根据本发明的一个实施方式用于模拟柔性的对象的***100的透视图。在医疗过程期间，***100被用来模拟导管，但其可用来模拟任何柔性的对象。***100包括人/计算机界面102、电子界面104和计算机106。
导管108由使用者操控，并且，响应于这种操控，将虚拟现实图像显示在计算机106的显示器110上。计算机106可以是任何类型的包括处理器和存储器的通用的或专门的计算机，其中存储器用来储存由处理器执行的指令。
除了导管108，人/计算机界面102包括屏障112和“中心线”114，导管108通过中心线114***到“身体”中。用屏障112来代表覆盖患者身体的皮肤部分。在一个实施方式中，屏障112由人体模型或其他身体或身体部分(如躯干、上肢或下肢)的仿真代表物(life-likerepresentation)而形成。中心线114***到屏障112中来为导管108提供从屏障112进入或移除点，并且允许操控在患者体内的导管108的远处部分，同时将组织损害减到最低。虽然因为导管108的末端对于医疗模拟不是必须的而在一个实施方式中除去了导管108的末端，以防止对人或道具(property)的任何潜在伤害，但是导管108可以是任何商业可得的导管。
导管108包括把手或“手柄”116和杆118。手柄116可以是任何用来操控导管108的常规装置，或手柄116本身可包括杆118。杆118是细长的柔性的对象，并且，具体地，是细长的圆柱形对象。为了模拟导管108，***100在三维空间中追踪杆118的运动，其运动被约束以便杆118只具有二、三或四个运动自由度。这是导管的典型用法的合理模拟，因为一旦导管被***到患者体内，在某种情况下导管将沿其长度被限制到约两个自由度上。
触觉界面120接收杆118并且在杆118上施加触觉反馈，触觉反馈可由使用者感受并且为使用者提供导管108进入实际身体的感觉。在一个实施方式中，触觉界面120包括一个或多个致动器和其他产生触觉反馈的装置。触觉界面120可以是用于在杆118上产生触觉反馈的任何已知装置，包括在美国专利第5,821,920号中公开的触觉界面。触觉界面120也确定导管108在模拟身体内的位置，包括使用者是否在手柄116处推、拉或扭转导管108。
电子界面104通过电缆122从触觉界面120接收位置信息，并通过电缆124将信息传输至计算机106。如在以下更详细地公开的，作为响应，计算机106对导管108的位置进行建模，并且在显示器110上产生模拟的图形图像。进一步地，计算机106基于导管108的位置产生所需要的触觉效果，并且向触觉界面120提供信号来产生由使用者感受的触觉效果。
在模拟柔性的对象例如导管108时，本发明的一个实施方式最先使用了用于生成有限元模型的四面体分解和节点索引方法。然后，本发明的实施方式使用直接数值解法(direct numerical solver)来执行快速模型模拟。
在一个实施方式中，把被模拟的柔性的对象建模为形成分段的空间曲线的连接的片段的串联链(serial chain)。除了限定片段相对于链中的前面片段的取向的两个角度外，每个片段还具有长度、半径和与其相关的材料参数。虽然在其他实施方式中，可对不同形状进行建模，但是在一个实施方式中，假设对象具有圆形横截面。
图2示出被模拟的柔性的对象的一个片段200。沿着对象的中心线或空间曲线220通过片段节点201、202将每个片段连接到其他片段。将其他周围的节点210-215放置在片段节点201、202周围来形成片段。然后，根据图2所示的样式，通过连接边和节点，将每个片段分解为三棱柱。图3示出如何通过***其他的边301-303将每个棱柱进一步分解为三个四面体元。
图4示出根据本发明的一个实施方式的两个连续的片段400和401以及相应的节点，并且示出节点索引方案。索引节点，以便由边连接的两个节点的索引值之间的差值保持由一常数界定，该常数相对于对象中的节点总数来说是小的。在一个实施方式中，以增加的顺序来为片段逐个分配索引值。如图4所示，如果将k-1个节点放置在每个中心线节点周围，那么就将索引值ki、ki+1、…、ki+k-1分配到属于片段i的节点。索引方法产生由边连接的两个节点的索引值之间的最大差值2k-1。可为非管状对象开发类似的方案。
在一个实施方式中，两个节点的索引值之间的最大差值由s＊k界定，其中s是在一节点处连接的片段的最大个数，并且k是属于一片段的节点的最大个数。因此，对于导管和其他长的柔性的对象，范围是2＊k。
在一个实施方式中，通过连通的无向图“G”来对片段的连通性进行建模，其中每个片段由图边表示，并且连接两个片段的所有对象节点由图节点表示。索引每个对象节点，以便第一对象节点的第一索引值和第二对象节点的第二索引值之间的最大差值由(1+b^l)＊k界定，其中第一对象节点和第二对象节点由一对象边连接。“l”是连通性图的生成树“T”的分枝层次(branching level)的个数，“b”是生成树T中一节点的子节点(children)的最大个数，以及“k”是属于一片段的对象节点的最大个数。在该实施方式中，以下参数可应用于被模拟的对象的实施方式。
(1)线性连通性：b＝1，l＝0-＞2＊k；
(2)单个环：b＝2，l＝1-＞3＊k；
(3)在一节点处连接的两个环：b＝4，l＝1-＞5＊k；
(4)在一边处连接的两个环：b＝3，l＝2-＞10＊k(由于其也由5＊k界定，因此10＊k不是紧密界限(tight bound))；
(5)完全3层次二元树：b＝2，l＝3-＞9＊k。
本发明的一个实施方式通过持续计算每个片段周围的节点的位置，如图4中所示的节点，来模拟柔性的对象的移动。对象的移动遵循牛顿第二定律：
Ma＝f(x，v)
其中x、v和a是包含节点的位置、速度和加速度的矢量，M是表示对象中质量分布的矩阵，以及f包括作用在节点处的内力和外力。在一个实施方式中，根据以高度交互的速度(如，＞30Hz)作用在对象上的力来更新节点中每一个的位置。在一个实施方式中，这是使用半隐式时间步方法(semi-implicit time-stepping method)通过对运动方程进行数值积分来完成的，例如在Baraff和Witkin的“Large Stepsin Cloth Simulation”SIGGRAPH(1998)(“Baraff”)中公开的方法。
图5是根据本发明的一个实施方式由图1的计算机106执行来模拟柔性的对象例如导管108的功能的流程图。在一个实施方式中，图5的功能是由储存在存储器中和由处理器执行的软件来实现的。在其他的实施方式中，功能可由硬件或任何硬件和软件的组合来执行。
在500处，为被模拟的柔性的对象建立四面体有限元模型和索引。在一个实施方式中，如结合以上图2-4所公开的来建立四面体有限元模型。
在520处，计算元素旋转。在一个实施方式中，使用共转弯曲方法计算旋转，例如在Mueller和Gross的“Interactive VirtualMaterials”，Graphics Interface(2004)中公开的方法。
在530处，计算节点力(f)和力微分(和
Figure G2007800386560D00072
)。在一个实施方式中，使用半隐式时间步方法来计算节点力和力微分，例如Baraff公开的方法。
在540处，为给定步长h按 b = h ( f + h ∂ f ∂ x v ) 计算“b”。在一个实施方式中，使用半隐式时间步方法来计算b，例如Baraff公开的方法。
在550处，为给定步长h按 A = [ M - h ∂ f ∂ v - h 2 ∂ f ∂ x ] 计算“A”。在一个实施方式中，使用半隐式时间步方法来计算A，例如Baraff公开的方法。
在560处，根据计算的A值和b值来求解下面的线性***：AΔv＝b。在一个实施方式中，使用直接解法来求解线性***。在一个实施方式中，直接解法是带状Cholesky矩阵分解算法，例如在Gene H.Golub和Charles F.Van Loan的“Matrix Computations”，约翰霍布金斯大学出版社(Johns Hopkins University Press)，第156页(1996)所公开的算法。在570处，使用在560处的解的结果，Δv，来更新模拟循环中的节点位置和速度。然后，循环回到520。
如图5所公开的柔性的对象的建模和模拟提供了许多相对现有技术来说的优点。就一个优点来说，隐式积分提供了高频动态响应所需要的必要的稳定性和性能。而用显式积分方法不可能获得类似水平的稳定性和性能。此外，相对于使用线性FE模型来说，共转弯曲提供了更准确的体积守恒。体积守恒能导致所模拟对象的更真实的物理性状。
此外，结合以上图2-4所公开的四面体有限元模型和索引的结果是***矩阵A的带宽为2k-1。因此，根据一个实施方式，该模型生成带有低带宽的非常稀疏的矩阵，其允许使用直接数值解法。图6用图形示出根据本发明的一个实施方式的***矩阵A的带状结构。如图所示，所有的非零矩阵元都在距对角线600有限的距离内。***矩阵A的带状结构是节点编号方法的直接结果。带状解法使用矩阵的带状结构来更快地求解***AΔv＝b。
在现有技术的方法中，通常通过迭代方法来求解线性***AΔv＝b，例如共轭梯度方法。相反地，本发明的实施方式使用带状Cholesky矩阵分解算法来高效地求解从图2-4所示的模型中获得的线性***。由于索引节点而产生的网络的连通性，直接解法优胜于现有技术的迭代解法。在一个实施方式中，算法的运行时间复杂度是O(m2n)，其中n是片段的个数以及m＝6k-3是***矩阵A的带宽。在一个实施方式中，***矩阵的非零元素储存在计算机存储器中的线性阵列中，还储存在限定矩阵中元素的位置的二维索引值的阵列中，以及储存在指标(pointer)的二维阵列中，其非零指标是指一维阵列中的元素。在一个实施方式中，修改了带状Cholesky算法，因此将在执行该算法期间生成的零元素从二维指标阵列中移除。
图7示出根据本发明的一个实施方式的所得到的模拟的四面体模型。模型710与模拟的身体器官部分720如心脏、静脉、动脉等并排示出。在图8中示出了根据本发明的一个实施方式由使用者看见的实际的模拟导管800。
如所公开的，本发明的一个实施方式通过建立四面体有限元模型和索引节点并且然后使用直接数值解法来执行快速模型模拟，来模拟柔性的对象，例如导管。其结果是有效并真实的对象模拟。
在此具体地说明和/或描述了本发明的几个实施方式。然而，应理解，本发明的修改和变化包含在以上所教授的内容中，并且在所附权利要求的范围之内而不偏离本发明的精神和预期的范围。

Claims (26)

1.一种模拟柔性的细长的对象的方法，其包括以下步骤：
将所述柔性细长对象建模为多个片段的串联链，每个片段具有由一个或多个边连接的多个节点；
将每个片段建模为由所述节点和所述边形成的多个四面体元，以产生四面体有限元模型；以及
索引每个节点，以便第一节点的第一索引值和第二节点的第二索引值之间的最大差值是有界的并且约等于s*k，其中s是在一个节点处连接的片段的最大个数，并且k是属于一个片段的节点的最大个数；
通过持续计算每个片段周围的节点的位置来模拟所述柔性细长对象的移动，包括：
为所述四面体有限元模型计算第一时间期的元素旋转；
为第一时间期和每个节点计算力；以及
在显示器上显示所模拟的柔性细长对象的图形图像。
2.如权利要求1所述的方法，其中所述柔性细长对象包括分段的空间曲线，并且沿着所述空间曲线通过片段节点连接每个片段。
3.如权利要求2所述的方法，其中每个片段包括所述片段节点周围的多个周围的节点。
4.如权利要求3所述的方法，其中将至多k-1个周围的节点放置在属于片段i的每个片段节点周围，并且将索引值ki、ki+1、...、ki+k-1分配到所述周围的节点。
5.如权利要求1所述的方法，其中所述计算步骤包括共转弯曲。
6.如权利要求1所述的方法，其中所述计算步骤包括求解AΔv＝b，并且其中使用半隐式时间步方法来计算A和b。
7.如权利要求6所述的方法，其中对于给定步长h按
Figure FSB0000118456510000021
计算b，并且对于所述给定步长h按
Figure FSB0000118456510000022
计算A。
8.如权利要求6所述的方法，其中所述求解步骤包括直接解法。
9.如权利要求8所述的方法，其中所述直接解法是带状Cholesky矩阵分解算法。
10.如权利要求9所述的方法，其中所述直接解法被修改以移除在执行期间生成的零元素。
11.如权利要求6所述的方法，其中A的带宽为m＝6k-3。
12.如权利要求1所述的方法，其中所述柔性细长对象是导管。
13.如权利要求1所述的方法，其中所述柔性细长对象是距近似分段的空间曲线有限距离的可变形体。
14.如权利要求1所述的方法，其中所述柔性细长对象可在几何上分解为一组层，其中每一层是距近似分段的空间曲线有限距离的可变形体。
15.如权利要求1所述的方法，其中所述柔性细长对象是距一组分段的空间曲线有限距离的可变形体，其中曲线片段的子集被连接以形成图。
16.一种用来模拟柔性的细长的对象的***，其包括：
用于将所述柔性细长对象建模为多个片段的串联链的装置，每个片段具有由一个或多个边连接的多个节点；
用于将每个片段建模为由所述节点和所述边形成的多个四面体元以产生四面体有限元模型的装置；
用于索引每个节点的装置，以便第一节点的第一索引值和第二节点的第二索引值之间的最大差值是有界的并且约等于s*k，其中s是在一个节点处连接的片段的最大个数，并且k是属于一个片段的节点的最大个数；
用于通过持续计算每个片段周围的节点的位置来模拟所述柔性细长对象的移动的装置，包括：
用于为所述四面体有限元模型计算第一时间期的元素旋转的装置；
用于为第一时间期和每个节点计算力的装置；以及
用于在显示器上显示所模拟的柔性细长对象的图形图像的装置。
17.如权利要求16所述的***，其中所述柔性细长对象包括分段的空间曲线，并且沿着所述空间曲线通过片段节点连接每个片段。
18.如权利要求17所述的***，其中每个片段包括所述片段节点周围的多个周围的节点。
19.如权利要求18所述的***，其中将至多k-1个周围的节点放置在属于片段i的每个片段节点周围，并且将索引值ki、ki+1、...、ki+k-1分配到所述周围的节点。
20.一种用于模拟柔性的细长的对象的***，其包括：
处理器；
存储器，其耦合到所述处理器；以及
显示器，其耦合到所述处理器；
其中所述存储器储存指令，当所述处理器执行所述指令时，所述指令使所述处理器：
将所述柔性细长对象建模为多个片段的串联链，每个片段具有由一个或多个边连接的多个节点；
将每个片段建模为由所述节点和所述边形成的多个四面体元，以产生四面体有限元模型；
索引每个节点，以便第一节点的第一索引值和第二节点的第二索引值之间的最大差值是有界的并且约等于s*k，其中s是在一个节点处连接的片段的最大个数，并且k是属于一个片段的节点的最大个数；
用于通过持续计算每个片段周围的节点的位置来模拟所述柔性细长对象的移动，包括：
为所述四面体有限元模型计算第一时间期的元素旋转；
为第一时间期和每个节点计算力；以及
在显示器上显示所模拟的柔性细长对象的图形图像。
21.如权利要求20所述的***，其中所述柔性细长对象包括分段的空间曲线，并且沿着所述空间曲线通过片段节点连接每个片段。
22.如权利要求21所述的***，其中每个片段包括所述片段节点周围的多个周围的节点。
23.如权利要求22所述的***，其中将至多k-1个周围的节点放置在属于片段i的每个片段节点周围，并且将索引值ki、ki+1、...、ki+k-1分配到所述周围的节点。
24.如权利要求20所述的***，其中所述柔性细长对象是导管。
25.一种模拟柔性的细长的对象的方法，其包括以下步骤：
将所述柔性细长对象建模为多个片段的串联链，每个片段具有由一个或多个对象边连接的多个对象节点；
将每个片段建模为由所述对象节点和所述对象边形成的多个四面体元，以产生四面体有限元模型；
通过连通的无向图G对所述片段的连通性进行建模，其中每个片段由图边表示，并且连接两个片段的所有对象节点由图节点表示；
索引每个对象节点，以便第一对象节点的第一索引值和第二对象节点的第二索引值之间的最大差值由(1+b^I)*k界定，其中所述第一对象节点和所述第二对象节点由对象边连接，I是所述图G的生成树T的分枝层次的个数，b是所述生成树T中一个节点的子节点的最大个数，并且k是属于一个片段的对象节点的最大个数；以及
通过持续计算每个片段周围的节点的位置来模拟所述柔性细长对象的移动，包括：
为所述四面体有限元模型计算第一时间期的元素旋转；
为第一时间期和每个节点计算力；以及
在显示器上显示所模拟的所述柔性细长对象的图形图像。
26.一种模拟柔性的细长的对象的方法，其包括以下步骤：
将所述柔性细长对象建模为多个片段的串联链，每个片段具有由一个或多个边连接的多个节点；
将每个片段建模为由所述节点和所述边形成的多个四面体元，以产生四面体有限元模型；以及
索引每个节点，以便第一节点的第一索引值和第二节点的第二索引值之间的最大差值是有界的并且大约小于所述多个节点的数量的一半，其中所述第一节点和所述第二节点由边连接；
通过持续计算每个片段周围的节点的位置来模拟所述柔性细长对象的移动，包括：
为所述四面体有限元模型计算第一时间期的元素旋转；
为第一时间期和每个节点计算力；以及
在显示器上显示所模拟的柔性细长对象的图形图像。
CN200780038656.0A 2006-10-16 2007-10-10 柔性的对象模拟器 Expired - Fee Related CN101529487B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11/549,834 2006-10-16
US11/549,834 US20080088578A1 (en) 2006-10-16 2006-10-16 Flexible object simulator
PCT/US2007/080913 WO2008048831A2 (en) 2006-10-16 2007-10-10 Flexible object simulator

Publications (2)

Publication Number Publication Date
CN101529487A CN101529487A (zh) 2009-09-09
CN101529487B true CN101529487B (zh) 2014-04-09

Family

ID=39314978

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200780038656.0A Expired - Fee Related CN101529487B (zh) 2006-10-16 2007-10-10 柔性的对象模拟器

Country Status (5)

Country Link
US (1) US20080088578A1 (zh)
EP (1) EP2080184A2 (zh)
JP (1) JP5144669B2 (zh)
CN (1) CN101529487B (zh)
WO (1) WO2008048831A2 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007030603A2 (en) 2005-09-08 2007-03-15 Wms Gaming Inc. Gaming machine having display with sensory feedback
WO2007117418A2 (en) 2006-03-31 2007-10-18 Wms Gaming Inc. Portable wagering game with vibrational cues and feedback mechanism
US8005659B2 (en) * 2007-07-06 2011-08-23 Immersion Medical, Inc. Simulation of coupled objects
US8442806B2 (en) 2010-03-03 2013-05-14 Immersion Medical, Inc. Systems and methods for simulations utilizing a virtual coupling
US20120100515A1 (en) * 2010-10-20 2012-04-26 Northwestern University Fluoroscopy Simulator
US9058714B2 (en) 2011-05-23 2015-06-16 Wms Gaming Inc. Wagering game systems, wagering gaming machines, and wagering gaming chairs having haptic and thermal feedback
US9449456B2 (en) 2011-06-13 2016-09-20 Bally Gaming, Inc. Automated gaming chairs and wagering game systems and machines with an automated gaming chair
CA2869079A1 (en) * 2013-10-29 2015-04-29 The Royal Institution For The Advancement Of Learning/Mcgill University Finite element methods and systems
CN111437033B (zh) * 2020-04-03 2021-03-02 天津理工大学 一种用于血管介入手术机器人***的虚拟传感器

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6929481B1 (en) * 1996-09-04 2005-08-16 Immersion Medical, Inc. Interface device and method for interfacing instruments to medical procedure simulation systems

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5821920A (en) * 1994-07-14 1998-10-13 Immersion Human Interface Corporation Control input device for interfacing an elongated flexible object with a computer system
US5594651A (en) * 1995-02-14 1997-01-14 St. Ville; James A. Method and apparatus for manufacturing objects having optimized response characteristics
JP3712451B2 (ja) * 1995-10-30 2005-11-02 株式会社ナムコ 画像合成方法及び画像合成装置
US7084884B1 (en) * 1998-11-03 2006-08-01 Immersion Corporation Graphical object interactions
AU2001270265A1 (en) * 2000-06-29 2002-01-14 Object Reservoir, Inc. Method and system for coordinate transformation to model radial flow near a singularity
JP3845682B2 (ja) * 2002-02-18 2006-11-15 独立行政法人理化学研究所 シミュレーション方法
US6965370B2 (en) * 2002-11-19 2005-11-15 Immersion Corporation Haptic feedback devices for simulating an orifice
WO2004063988A1 (en) * 2003-01-15 2004-07-29 Koninklijke Philips Electronics N.V. Image processing method for automatic adaptation of 3-d deformable model onto a subtantially tubular surface of a 3-d object
JP4169263B2 (ja) * 2003-04-07 2008-10-22 株式会社バンダイナムコゲームス 画像生成装置及び情報記憶媒体
JP2005027587A (ja) * 2003-07-08 2005-02-03 Shimano Inc 竿体及び釣竿
JP4266783B2 (ja) * 2003-11-14 2009-05-20 関東自動車工業株式会社 柔軟物の変形解析装置
JP2006164219A (ja) * 2004-11-09 2006-06-22 Phifit Kk 有限要素解析用インターフェース、有限要素解析用インターフェースのプログラムを記録したコンピュータ読み取り可能な記録媒体、高付加価値通信網を利用した有限要素解析方法、有限要素解析用並列処理計算機および有限要素解析用計算機
JP4601405B2 (ja) * 2004-11-29 2010-12-22 株式会社日本総合研究所 ハイドロフォーミングシミュレーションシステム及びハイドロフォーミングシミュレーション用プログラム

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6929481B1 (en) * 1996-09-04 2005-08-16 Immersion Medical, Inc. Interface device and method for interfacing instruments to medical procedure simulation systems

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
A hybrid elastic model for real-time cutting,deformations, and force feedback for surgery training and simulation;Cotin et al.;《The Visual Computer》;20001231;第16卷(第8期);第441页第4.2部分、附图2 *
Baraff et al..Large Steps in Cloth Simulation.《SIGGRAPH, Computer Graphics Proceedings, Annual Conference Series》.1998,第49页第5.2部分,第44页第1.2部分.
Cotin et al..A hybrid elastic model for real-time cutting,deformations, and force feedback for surgery training and simulation.《The Visual Computer》.2000,第16卷(第8期),第441页第4.2部分、附图2.
Cotin et al..New Approaches to Catheter Navigation for Interventional Radiology Simulation.《MICCAI(2005),LNCS 3750》.2005,第536页第2.1-2.2部分、附图1,5.
Large Steps in Cloth Simulation;Baraff et al.;《SIGGRAPH, Computer Graphics Proceedings, Annual Conference Series》;19981231;第49页第5.2部分,第44页第1.2部分 *
New Approaches to Catheter Navigation for Interventional Radiology Simulation;Cotin et al.;《MICCAI(2005),LNCS 3750》;20051231;第536页第2.1-2.2部分、附图1,5 *

Also Published As

Publication number Publication date
JP2010507133A (ja) 2010-03-04
JP5144669B2 (ja) 2013-02-13
WO2008048831A2 (en) 2008-04-24
US20080088578A1 (en) 2008-04-17
WO2008048831A3 (en) 2008-07-10
CN101529487A (zh) 2009-09-09
EP2080184A2 (en) 2009-07-22

Similar Documents

Publication Publication Date Title
CN101529487B (zh) 柔性的对象模拟器
Lenoir et al. Interactive physically-based simulation of catheter and guidewire
Chentanez et al. Interactive simulation of surgical needle insertion and steering
Cotin et al. New approaches to catheter navigation for interventional radiology simulation
Duriez et al. New approaches to catheter navigation for interventional radiology simulation
Peterlik et al. Constraint-based haptic rendering of multirate compliant mechanisms
US20140272865A1 (en) Physics Engine for Virtual Reality Surgical Training Simulator
US8005659B2 (en) Simulation of coupled objects
Wang et al. A robust and fast approach to simulating the behavior of guidewire in vascular interventional radiology
Luo et al. A robust and real-time vascular intervention simulation based on Kirchhoff elastic rod
Kuroda et al. Interaction model between elastic objects for haptic feedback considering collisions of soft tissue
US20090130643A1 (en) Method for simulating a manual interventional operation by a user in a medical procedure
CN104574503A (zh) 造影剂扩散过程模拟方法和装置
Wang et al. A computer-based real-time simulation of interventional radiology
Cai et al. VR simulated training for less invasive vascular intervention
Cotin et al. Volumetric deformable models for simulation of laparoscopic surgery
Li et al. A Novel FEM‐Based Numerical Solver for Interactive Catheter Simulation in Virtual Catheterization
Guo et al. Tensor-mass Model based real-time simulation of vessel deformation and force feedback for the interventional surgery training system
Zhang et al. Real-time simulation of deformable soft tissue based on mass-spring and medial representation
Nakao et al. Haptic reproduction and interactive visualization of a beating heart for cardiovascular surgery simulation
Cheng et al. A fast and stable guidewire model for minimally invasive vascular surgery based on Lagrange multipliers
Li et al. Haptics-equiped interactive PCI simulation for patient-specific surgery training and rehearsing.
Zuo et al. Position based catheterization and angiography simulation
Aloisio et al. Computer-based simulator for catheter insertion training
Westwood Physics-based models for catheter, guidewire and stent simulation

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
C56 Change in the name or address of the patentee

Owner name: YIMEISEN COMPANY

Free format text: FORMER NAME: IMMERSION CORP.

CP01 Change in the name or title of a patent holder

Address after: American California

Patentee after: IMMERSION CORPORATION

Address before: American California

Patentee before: Immersion Corp.

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: 20140409

Termination date: 20191010