CN106446452B - 基于力触觉交互的拱梁弹簧形变模型的建模方法 - Google Patents

基于力触觉交互的拱梁弹簧形变模型的建模方法 Download PDF

Info

Publication number
CN106446452B
CN106446452B CN201610910611.9A CN201610910611A CN106446452B CN 106446452 B CN106446452 B CN 106446452B CN 201610910611 A CN201610910611 A CN 201610910611A CN 106446452 B CN106446452 B CN 106446452B
Authority
CN
China
Prior art keywords
arched girder
spring
girder
deformation
virtual
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
Application number
CN201610910611.9A
Other languages
English (en)
Other versions
CN106446452A (zh
Inventor
孙伟
吕查德
张小瑞
刘佳
朱利丰
宋爱国
牛建伟
王澎湃
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.)
NANJING LOONG SHIELD INTELLIGENT TECHNOLOGY CO.,LTD.
Original Assignee
Nanjing University of Information Science and Technology
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 Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201610910611.9A priority Critical patent/CN106446452B/zh
Publication of CN106446452A publication Critical patent/CN106446452A/zh
Application granted granted Critical
Publication of CN106446452B publication Critical patent/CN106446452B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Bridges Or Land Bridges (AREA)

Abstract

本发明公开了一种基于力触觉交互的拱梁弹簧形变模型的建模方法,首先对虚拟场景进行初始化;然后进行位置检测,当虚拟代理碰撞到虚拟柔性体表面时,在给定虚拟接触压力F的作用下,利用拱梁弹簧形变模型来填充虚拟代理与虚拟柔性体交互的局部区域内部;接着,利用拱梁弹簧形变模型计算虚拟柔性体局部区域变形量;最后根据计算出的虚拟柔性体表面的形变刷新图形,并反馈输出力触觉信息。本发明变形计算速度快,计算量少,且通过改变建模方法中相应的参数,就可对不同的柔性体进行变形仿真,建模方法适用性广。

Description

基于力触觉交互的拱梁弹簧形变模型的建模方法
技术领域
本发明涉及一种力触觉交互的形变模型的建模方法,尤其涉及一种用于虚拟手术仿真过程中,基于力触觉交互的拱梁弹簧形变模型的建模方法。
背景技术
虚拟手术作为虚拟现实技术在医学领域的重要应用之一,从上世纪八九十年代起越来越受到国内外学界及工业界的关注。由于真实的手术成本高、风险大,通过对手术过程进行计算机仿真,虚拟手术能够有效地测试手术计划的可行性,辅助手术医生及学徒演练实际的操作过程,与先进的图形技术相结合还能帮助病患与手术医生进行沟通。随着中国人口老龄化浪潮和信息时代的到来,虚拟手术的发展也被推到了一个关键的时刻。
尽管从虚拟手术概念的提出到现在已有二十多年,涵盖心肺、肝脏、五官、骨科、牙科等不同领域的手术仿真技术均被广泛研究,虚拟手术***在医学训练诊断中的应用仍未被推广,其主要难点在于精准的虚拟手术仿真需要消耗大量计算成本。尤其是针对软组织的虚拟手术***,欲达到手术仿真中形变的真实感需要消耗大量的建模资源和计算资源。因此,提出准确快速的力触觉交互的形变模型的建模方法是虚拟手术成功的关键。采用基于物理意义的建模方法,对柔性体的变形进行实时仿真,已成为近几年研究的热点。
目前常用的基于物理意义的柔性体形变模型的建模方法中,有限元虽然精确,但要求网格的节点数多、计算时耗多,实时性较差。弹簧-质点网络构建简单直观,计算量小,但由于弹簧系数设定的任意性,导致所仿真物体行为特性的优劣较难定量地鉴别。边界元通过对模型边界进行离散来降低问题维数,简化计算,但在稳定性方面却存在一定困难。近年来,为提高虚拟手术仿真的交互性与实时性,许多采用GPU加速的有限元和非线性有限元计算的工作也相继开展,前者是建立在线性描述,小变形的基础之上的,后者虽能达到仿真的实时性要求,但模拟效果还有待提高。以上说明,这些常用的柔性体形变模型的建模方法均存在计算较为繁杂和模拟效果有待提高等问题。
发明内容
本发明所要解决的技术问题是针对背景技术中所涉及到的缺陷,提供一种基于力触觉交互的拱梁弹簧形变模型的建模方法。
本发明为解决上述技术问题采用以下技术方案:
基于力触觉交互的拱梁弹簧形变模型的建模方法,包含以下步骤:
步骤1),对虚拟场景进行初始化;
步骤2),位置检测,当虚拟物品碰撞到虚拟柔性体表面时,在给定虚拟接触压力F的作用下,利用拱梁弹簧形变模型来填充虚拟物品与虚拟柔性体交互的局部区域内部;
所述拱梁弹簧形变模型包含N层拱梁结构和N-1层弹簧结构,N为大于等于2的自然数,其中:
第一层拱梁结构包含1个受力不发生形变且左右对称的刚性拱梁,且该刚性拱梁的顶点悬挂在虚拟物品和虚拟柔性体表面的碰撞点;
对于每一个大于等于2且小于等于N的自然数i,第i层拱梁结构包含i个受力不发生形变且左右对称的刚性拱梁,该i个刚性拱梁的底端依次相连形成波浪线,具有i个顶端和i+1个底端;第i-1层弹簧结构包含i个弹簧,一一对应垂直设置在第i-1层拱梁结构的i个底端和第i层拱梁结构的i个顶端之间;
第i-1层弹簧结构中第一个弹簧的弹性系数ki-1,1和第一层弹簧结构中弹簧的弹性系数k1,1满足下面公式:
第i-1层弹簧结构中第j个弹簧的弹性系数与该层第一个弹簧的弹性系数之间满足如下关系式:
其中,j为大于等于2小于等于i的自然数;
步骤3),利用拱梁弹簧形变模型计算虚拟柔性体局部区域变形量;
步骤3.1),计算每一层拱梁结构中受力最大的刚性拱梁所受的分力值:
步骤3.2),依次将每一层拱梁结构中受力最大的刚性拱梁的所受的分力值和预设的力量阈值f进行比较,直至第M层拱梁结构中受力最大的刚性拱梁的所受的分力值小于预设的力量阈值f,M为大于等于1小于等于N的自然数;
步骤3.3),依次计算第一层弹簧结构到第M-1层弹簧结构的形变;
步骤3.4),将第一层弹簧结构到第M-1层弹簧结构的形变进行累加,得到虚拟柔性体表面的形变;
步骤4),根据计算出的虚拟柔性体表面的形变刷新图形。
作为本发明一种基于力触觉交互的拱梁弹簧形变模型的建模方法进一步的优化方案,步骤3.1)中根据以下公式计算每一层拱梁结构中受力最大的刚性拱梁所受的分力值:
当n为偶数时,第n层拱梁结构中第个刚性拱梁所受分力最大,为:
当n为奇数时,第n层拱梁结构中第个刚性拱梁所受分力最大,为:
n为大于等于1小于等于N的自然数。
作为本发明一种基于力触觉交互的拱梁弹簧形变模型的建模方法进一步的优化方案,步骤3.3)中根据下列公式依次计算第一层弹簧结构到第M-1层弹簧结构的形变:
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)和以往常用的基于物理意义的柔性体变形仿真力触觉建模方法相比,该建模方法变形计算速度加快,减少了计算量。
(2)通过改变建模方法中相应的参数,就可对不同的柔性体进行变形仿真,建模方法适用性广。
附图说明
图1是柔性体变形仿真流程图;
图2是拱梁弹簧形变模型中各层刚性拱梁所受分力大小关系示意图;
图3是拱梁弹簧形变模型的原始状态示意图;
图4是拱梁弹簧形变模型的按压下的状态示意图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
如图1所示,本发明公开了一种基于力触觉交互的拱梁弹簧形变模型的建模方法,包含以下步骤:
步骤1),对虚拟场景进行初始化;
步骤2),位置检测,当虚拟物品碰撞到虚拟柔性体表面时,在给定虚拟接触压力F的作用下,利用拱梁弹簧形变模型来填充虚拟物品与虚拟柔性体交互的局部区域内部;
所述拱梁弹簧形变模型包含N层拱梁结构和N-1层弹簧结构,如图3所示,N为大于等于2的自然数,其中:
第一层拱梁结构包含1个受力不发生形变且左右对称的刚性拱梁,且该刚性拱梁的顶点悬挂在虚拟物品和虚拟柔性体表面的碰撞点;
对于每一个大于等于2且小于等于N的自然数i,第i层拱梁结构包含i个受力不发生形变且左右对称的刚性拱梁,该i个刚性拱梁的底端依次相连形成波浪线,具有i个顶端和i+1个底端;第i-1层弹簧结构包含i个弹簧,一一对应垂直设置在第i-1层拱梁结构的i个底端和第i层拱梁结构的i个顶端之间;
第i-1层弹簧结构中第一个弹簧的弹性系数ki-1,1和第一层弹簧结构中弹簧的弹性系数k1,1满足下面公式:
第i-1层弹簧结构中第j个弹簧的弹性系数与该层第一个弹簧的弹性系数之间满足如下关系式:
其中,j为大于等于2小于等于i的自然数;
步骤3),利用拱梁弹簧形变模型计算虚拟柔性体局部区域变形量;
步骤3.1),假设每一层拱梁结构中每个刚性拱梁受力如图4所示,所受的分力值如图2所示,且满足下列公式,根据公式计算每一层拱梁结构中每个刚性拱梁所受的分力值以及受力最大的刚性拱梁所受的分力值:
第n层拱梁结构中第m个刚性拱梁所受分力大小为:
当n为偶数时,第n层拱梁结构中第个刚性拱梁所受分力最大,为:
当n为奇数时,第n层拱梁结构中第个刚性拱梁所受分力最大,为:
n为大于等于1小于等于N的自然数,m为大于等于1小于等于n的自然数。
步骤3.2),依次将每一层拱梁结构中受力最大的刚性拱梁的所受的分力值和预设的力量阈值f进行比较,直至第M层拱梁结构中受力最大的刚性拱梁的所受的分力值小于预设的力量阈值f,M为大于等于1小于等于N的自然数;
步骤3.3),根据下列公式依次计算第一层弹簧结构到第M-1层弹簧结构的形变;
步骤3.4),根据下列公式将第一层弹簧结构到第M-1层弹簧结构的形变进行累加,得到虚拟柔性体表面的形变;
其中,Δx为前M-1层弹簧产生的形变之和,即为虚拟柔性体表面的形变,Δxp-1为第p-1层弹簧结构产生的压缩变形量,kp-1,1为第p-1层弹簧结构中第一个弹簧的弹性系数,p为大于等于1小于等于M的自然数。
步骤4),根据计算出的虚拟柔性体表面的形变刷新图形,并反馈输出力触觉信息。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.基于力触觉交互的拱梁弹簧形变模型的建模方法,其特征在于,包含以下步骤:
步骤1),对虚拟场景进行初始化;
步骤2),位置检测,当虚拟物品碰撞到虚拟柔性体表面时,在给定虚拟接触压力F的作用下,利用拱梁弹簧形变模型来填充虚拟物品与虚拟柔性体交互的局部区域内部;
所述拱梁弹簧形变模型包含N层拱梁结构和N-1层弹簧结构,N为大于等于2的自然数,其中:
第一层拱梁结构包含1个受力不发生形变且左右对称的刚性拱梁,且该刚性拱梁的顶点悬挂在虚拟物品和虚拟柔性体表面的碰撞点;
对于每一个大于等于2且小于等于N的自然数i,第i层拱梁结构包含i个受力不发生形变且左右对称的刚性拱梁,该i个刚性拱梁的底端依次相连形成波浪线,具有i个顶端和i+1个底端;第i-1层弹簧结构包含i个弹簧,一一对应垂直设置在第i-1层拱梁结构的i个底端和第i层拱梁结构的i个顶端之间;
第i-1层弹簧结构中第一个弹簧的弹性系数ki-1,1和第一层弹簧结构中弹簧的弹性系数k1,1满足下面公式:
第i-1层弹簧结构中第j个弹簧的弹性系数与该层第一个弹簧的弹性系数之间满足如下关系式:
其中,j为大于等于2小于等于i的自然数;
步骤3),利用拱梁弹簧形变模型计算虚拟柔性体局部区域变形量;
步骤3.1),计算每一层拱梁结构中受力最大的刚性拱梁所受的分力值:
步骤3.2),依次将每一层拱梁结构中受力最大的刚性拱梁的所受的分力值和预设的力量阈值f进行比较,直至第M层拱梁结构中受力最大的刚性拱梁的所受的分力值小于预设的力量阈值f,M为大于等于1小于等于N的自然数;
步骤3.3),依次计算第一层弹簧结构到第M-1层弹簧结构的形变;
步骤3.4),将第一层弹簧结构到第M-1层弹簧结构的形变进行累加,得到虚拟柔性体表面的形变;
步骤4),根据计算出的虚拟柔性体表面的形变刷新图形。
2.根据权利要求1所述的基于力触觉交互的拱梁弹簧形变模型的建模方法,其特征在于,步骤3.1)中根据以下公式计算每一层拱梁结构中受力最大的刚性拱梁所受的分力值:
当n为偶数时,第n层拱梁结构中第个刚性拱梁所受分力最大,为:
当n为奇数时,第n层拱梁结构中第个刚性拱梁所受分力最大,为:
n为大于等于1小于等于N的自然数。
3.根据权利要求1所述的基于力触觉交互的拱梁弹簧形变模型的建模方法,其特征在于,步骤3.3)中根据下列公式依次计算第一层弹簧结构到第M-1层弹簧结构的形变;
CN201610910611.9A 2016-10-19 2016-10-19 基于力触觉交互的拱梁弹簧形变模型的建模方法 Active CN106446452B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610910611.9A CN106446452B (zh) 2016-10-19 2016-10-19 基于力触觉交互的拱梁弹簧形变模型的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610910611.9A CN106446452B (zh) 2016-10-19 2016-10-19 基于力触觉交互的拱梁弹簧形变模型的建模方法

Publications (2)

Publication Number Publication Date
CN106446452A CN106446452A (zh) 2017-02-22
CN106446452B true CN106446452B (zh) 2019-03-29

Family

ID=58177139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610910611.9A Active CN106446452B (zh) 2016-10-19 2016-10-19 基于力触觉交互的拱梁弹簧形变模型的建模方法

Country Status (1)

Country Link
CN (1) CN106446452B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101976298A (zh) * 2010-09-27 2011-02-16 南京信息工程大学 增强力触觉再现的对称式板弹簧虚拟模型的建模方法
CN102207997A (zh) * 2011-06-07 2011-10-05 哈尔滨工业大学 基于力反馈的机器人微创手术仿真***
CN103869983A (zh) * 2014-03-26 2014-06-18 南京信息工程大学 一种用于力触觉人机交互的柔性物体变形仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060277466A1 (en) * 2005-05-13 2006-12-07 Anderson Thomas G Bimodal user interaction with a simulated object

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101976298A (zh) * 2010-09-27 2011-02-16 南京信息工程大学 增强力触觉再现的对称式板弹簧虚拟模型的建模方法
CN101976298B (zh) * 2010-09-27 2012-06-27 南京信息工程大学 增强力触觉再现的对称式板弹簧虚拟模型的建模方法
CN102207997A (zh) * 2011-06-07 2011-10-05 哈尔滨工业大学 基于力反馈的机器人微创手术仿真***
CN103869983A (zh) * 2014-03-26 2014-06-18 南京信息工程大学 一种用于力触觉人机交互的柔性物体变形仿真方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
手术仿真中基于导纳控制的力触觉形变模型;张小瑞、孙伟 等;《计算机辅助设计与图形学学报》;20150115;第114-119页

Also Published As

Publication number Publication date
CN106446452A (zh) 2017-02-22

Similar Documents

Publication Publication Date Title
CN110289104B (zh) 软组织按压和形变恢复的模拟方法
CN105302972A (zh) 一种基于元球模型的软组织形变方法
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.
Hou et al. A new model of soft tissue with constraints for interactive surgical simulation
CN111488670B (zh) 一种非线性的质点弹簧软组织形变仿真方法
Joldes et al. An efficient hourglass control implementation for the uniform strain hexahedron using the total Lagrangian formulation
Liu et al. Real-time simulation of virtual palpation system
CN108710735A (zh) 一种实时交互的无网格软组织形变模拟方法
CN103745058A (zh) 一种任意形状软组织表皮上受拉力/变形的模拟方法
Dehghan et al. A three‐dimensional large deformation model for soft tissue using meshless method
CN107146288A (zh) 虚拟手术中支持实时按压形变的软组织模型建模方法
CN107080586B (zh) 虚拟手术中用于软组织实时形变的交叉弹簧模型
Banihani et al. POD for real‐time simulation of Hyperelastic soft biological tissue using the point collocation method of finite spheres
CN106446452B (zh) 基于力触觉交互的拱梁弹簧形变模型的建模方法
Marinkovic et al. Towards real-time simulation of deformable structures by means of co-rotational finite element formulation
CN106202689A (zh) 一种软组织有限元模型的加速计算方法
CN109598799A (zh) 一种基于CycleGAN的虚拟切割方法
Wang et al. An unfixed-elasticity mass spring model based simulation for soft tissue deformation
Liu et al. Modelling and simulation of vascular tissue based on finite element method
Shi et al. Virtual surgical bleeding simulation with navier-stokes equation and modified smooth particle hydrodynamics method
JP2003141566A (ja) 3次元物体の切断シミュレーション方法
Wang et al. Real-time cutting method for soft tissue based on TLED algorithm
Ling et al. An improved meshless method in virtual surgery 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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211129

Address after: 210000 Building 1, 116 Shiyang Road, Qinhuai District, Nanjing City, Jiangsu Province

Patentee after: NANJING LOONG SHIELD INTELLIGENT TECHNOLOGY CO.,LTD.

Address before: 210044, No. 219, Ning six road, Pukou District, Jiangsu, Nanjing

Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY