CN101964023B - 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 - Google Patents
柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 Download PDFInfo
- Publication number
- CN101964023B CN101964023B CN2010102991871A CN201010299187A CN101964023B CN 101964023 B CN101964023 B CN 101964023B CN 2010102991871 A CN2010102991871 A CN 2010102991871A CN 201010299187 A CN201010299187 A CN 201010299187A CN 101964023 B CN101964023 B CN 101964023B
- Authority
- CN
- China
- Prior art keywords
- asymmetry
- spring
- layer
- veneer
- veneer spring
- 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
Links
Images
Landscapes
- Veneer Processing And Manufacture Of Plywood (AREA)
Abstract
本发明公开了一种柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法,其特征是当检测到虚拟代理碰撞到虚拟柔性体表面上任何一点时,在给定虚拟接触拉力作用下,虚拟代理与虚拟柔性体交互的局部区域内部填充非对称形板弹簧虚拟模型,在交互过程中,输出反馈为采用非对称形板弹簧虚拟模型计算出来的反应在拉力作用下柔性体实时变形仿真的力触觉信息的信号;非对称形板弹簧虚拟模型中每层的非对称形单板弹簧变形量之和的叠加对外等效为虚拟柔性体表面的变形;该建模方法不仅计算简单,而且能够保证变形计算具有较高精度,仿真实验表明其能够提供实时的力触觉反馈,具有很强的真实感,在变形仿真过程中具有较好的稳定性和实时性。
Description
技术领域
本发明涉及一种力触觉再现的建模方法,尤其涉及一种用于虚拟手术仿真过程中,人机交互的基于物理意义的柔性体实时力触觉再现的非对称形板弹簧虚拟模型的建模方法。
背景技术
虚拟手术在实际应用中具有远程干预、合理制订手术方案、降低手术费用、等优点,并且提供了一种脱离实体的安全高效的学习和训练手段,随着计算机技术的不断提高和虚拟现实领域研究的不断深入,虚拟手术在仿真的逼真性和实时性上将不断提高,它必将在未来的手术仿真领域中得到更广泛的应用,并将发挥巨大的作用。
力触觉是人类感觉器官的重要组成部分,在虚拟手术中,引入力触觉反馈,可以使医生在手术仿真中不仅能够看到而且还能感觉到手术器官,同时医生在进入虚拟手术环境时,不仅能够通过手和手臂的运动与虚拟组织器官进行交互,从而形成对虚拟手术环境的一个完整认识,而且能像操作真实物体一样,真切感受到与虚拟对象交互时产生的力触觉,这无疑可使手术训练更真实、准确、可靠。虚拟手术作为一种新式的手术训练手段,准确快速的力触觉再现建模方法是虚拟手术成功的关键。采用基于物理意义的力触觉再现的建模方法,对柔性体的变形进行实时仿真,已成为该领域应用与研究的最新热点,并将对虚拟手术的发展产生深远的影响。
目前,常用的基于物理意义的柔性体实时力触觉再现的建模方法主要有如下5种:①弹簧-质点模型具有网格构建直观、简单、联合方程组易于求解的优点,但递归求解时容易引起不稳定的振荡现象,造成仿真结果不够精确。②有限元模型具有较高的仿真精度和真实度,然而计算量较大成为其应用于实时仿真的主要障碍。③锁链模型离散化精度依然不高,且形变由锁链的几何结构来决定,因而所仿真的物体不够真实,此外,该模型的实现不可避免地需要较高的计算量。④边界元模型在求解过程中,虽***仿真规模较小,但在稳定性方面却存在一定的困难。⑤形状保持链模型具有实时性较好的优点,但是基于链元素为刚性体这个假设,故精度不高。以上说明,这些常用的柔性体实时力触觉再现的建模方法,均存在计算较为繁杂和仿真精度不高等问题,继而影响了计算的实时性和有效性。再加上由于人体柔性体的复杂性、差异性和多样性,使现有常用的基于物理意义的柔性体实时力触觉再现的建模方法在变形效果上不够理想。
鉴于上述问题,为了使虚拟手术仿真过程中力触觉人机交互过程更加符合人们自身的习惯,提高交互的沉浸感和真实感,提出了柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法。
发明内容
本发明提出一种基于物理意义的柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法,并适用于虚拟人机交互的柔性体实时变形仿真过程中。该建模方法计算简单,能准确快速的计算变形量,实现对柔性体的实时变形仿真;且在交互过程中,操作者能实时自然地感知到仿真过程中的力触觉信息,提高了虚拟手术仿真过程中力触觉交互的逼真度。
本发明采用如下技术方案:
1、一种柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法,其特征在于该建模方法的主要步骤如下:
步骤1对虚拟场景进行初始化;
步骤2当检测到虚拟代理碰撞到虚拟柔性体表面上任何一点时,在给定虚拟接触拉力F作用下,虚拟代理与虚拟柔性体交互的局部区域内部填充非对称形板弹簧虚拟模型,在交互过程中,输出反馈为采用非对称形板弹簧虚拟模型计算出来的反应在拉力作用下虚拟柔性体实时变形仿真的力触觉信息的信号,所述非对称形板弹簧虚拟模型的建模方法为:
(1)参数初始化,
(2)在给定虚拟接触拉力F作用下,当虚拟代理碰撞到虚拟柔性体表面上的任意点时,在碰撞点处下悬挂一叶片宽度为b、厚度为h、长度为l1的非对称形单板弹簧,形成第一层,在第一层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为l2=l1q的非对称形单板弹簧,形成第二层,在第二层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为l3=l1q2的非对称形单板弹簧,形成第三层,依此类推,在第i-1层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为li=l1qi-1的非对称形单板弹簧,形成第i层,i=1,2,3,……N,N为自然数,各层非对称形单板弹簧长度构成以l1为首项,以q为公比的等比数列,q=1~2,用一宽度为B的簧箍,将组成非对称形板弹簧的所有非对称形单板弹簧夹紧在一起,且簧箍可在非对称形板弹簧上移动,b≤B≤b+0.1,
以虚拟接触拉力F的作用线为分界线,两侧非对称形板弹簧的叶片数相等,叶片厚度相同,分界线将第i层的非对称形单板弹簧长度li分成li′、li″左右两部分,且左右两部分长度的比值相同,即:
l1′、l1″,l2′、l2″,l3′、l3″,li′、li″分别表示分界线将第一、二、三、i层的非对称形单板弹簧长度l1、l2、l3、li分成左右两部分,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
采用簧箍夹紧时,组成柔性力触觉再现的非对称形板弹簧虚拟模型中所有非对称形单板弹簧均有一段无效部分,在计算弹簧刚度时,应考虑其影响并进行必要的修正,修正后第i层的非对称形单板弹簧的弹簧刚度为:
其中,li、Pi、B分别为第i层的非对称形单板弹簧的长度、第i层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
第i层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度Pi满足:
Pi=Pi′λ(1+λ) (3)
Pi′为第i层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,其满足:
其中,E、I、li′分别为弹性模量、截面惯性矩、第i层的非对称形单板弹簧长度li以分界线为界左侧的长度,
截面惯性矩I满足:
在虚拟接触拉力F作用下,若柔性体中共有M层的非对称形单板弹簧产生变形,则第M层称为变形截止层,
若给定的虚拟接触拉力F能使第一层的非对称形单板弹簧被拉伸时产生的变形量X1达到第一层的非对称形单板弹簧给定的挠度值这种情况下假定前M-1层的非对称形单板弹簧被拉伸时产生的变形量均与第一层的非对称形单板弹簧给定的挠度值相同,第M层的非对称形单板弹簧被拉伸时产生的变形量不大于第一层的非对称形单板弹簧给定的挠度值,
第一层的非对称形单板弹簧消耗的拉力F1为:
XC1表示与第一层的非对称形单板弹簧的宽度、厚度、长度均相同的对应的对称式单板弹簧给定的挠度值,λ表示每层的非对称形单板弹簧以分界线为界,分成左右两部分长度的比值,且每层分成左右两部分长度的比值均相同,
除第一层和变形截止层第M层外,其余各层的非对称形单板弹簧消耗的拉力Fj为:
j的取值范围为2~M-1,
lj、Pj、B分别为第j层的非对称形单板弹簧的长度、第j层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
lj=l1qj-1 (11)
第j层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度Pj满足:
Pj=Pj′λ(1+λ) (12)
Pj′表示第j层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,E、I、lj′分别为弹性模量、截面惯性矩、第j层的非对称形单板弹簧长度lj以分界线为界左侧的长度,
lj=l1qj-1 (15)
截面惯性矩I满足:
变形截止层第M层的变形量为:
lM、PM、B分别为变形截止层第M层的非对称形单板弹簧的长度、变形截止层第M层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
lM=l1qM-1 (19)
变形截止层第M层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度PM满足:
PM=PM′λ(1+λ) (20)
PM′表示第M层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,E、I、lM′分别为弹性模量、截面惯性矩、第M层的非对称形单板弹簧长度lM以分界线为界左侧的长度,
lM=l1qM-1 (23)
截面惯性矩I满足:
(3)使虚拟接触拉力F作用于虚拟柔性体碰撞点,第i层对应的非对称形单板弹簧被拉伸,如果前i层所有的非对称形单板弹簧消耗的拉力Fi之和小于给定的虚拟接触拉力F,且前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间满足刷新频率1000Hz以上的要求,设共计经过时延Ti后,第i层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值,只有当第i层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值后,第i+1层对应的非对称形单板弹簧才开始被拉伸,依此类推,直到前M层所有的非对称形单板弹簧消耗的拉力之和不小于给定的虚拟接触拉力,或第M层的非对称形单板弹簧产生拉伸变形总计需要的时延时间不满足刷新频率的要求;
用ti、Ti分别表示第i层的非对称形单板弹簧产生拉伸变形需要的时延时间、前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间,且令层间的时延时间满足以第一层的非对称形单板弹簧产生拉伸变形需要的时延时间t1为首项、以w为公比的等比数列:
ti=wi-1t1 (25)
从虚拟代理碰撞接触到虚拟柔性体表面算起,假定前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间Ti满足Ti<T,其中
T为力触觉再现刷新频率的倒数;
所述非对称形板弹簧虚拟模型中每层的非对称形单板弹簧变形量之和的叠加对外等效为虚拟柔性体表面的变形,
其中,X为前M层所有的非对称形单板弹簧产生的拉伸变形量之和,Xi为前M-1层中任意第i层的非对称形单板弹簧产生的拉伸变形量,XM为变形截止层第M层产生的拉伸变形量,F为给定的虚拟接触拉力,为修正后第一层的非对称形单板弹簧刚度。
本发明的优点:
(1)和以往常用的基于物理意义的柔性体变形仿真力触觉建模方法相比,该建模方法中,采用一宽度为B的簧箍,将组成非对称形板弹簧的所有非对称形单板弹簧夹紧在一起,且簧箍可在非对称形板弹簧上移动,从而使该模型的使用范围更加广泛。
(2)组成非对称形板弹簧的所有非对称形单板弹簧以虚拟接触拉力F的作用线为分界线,两侧非对称形板弹簧的叶片数相等,叶片厚度相同,每层的非对称形单板弹簧以分界线为界分成左右两部分,且每层分成左右两部分的长度的比值λ均相同,简化了计算过程。
(3)若给定的虚拟接触拉力F能使第一层的非对称形单板弹簧被拉伸时产生的变形量X1达到第一层的非对称形单板弹簧给定的挠度值这种情况下假定前M-1层的非对称形单板弹簧被拉伸时产生的变形量均与第一层的非对称形单板弹簧给定的挠度值相同,从而进一步减少了计算量,使变形计算速度加快。
(4)通过改变建模方法中组成非对称形板弹簧的叶片宽度b、厚度h、第一层的非对称形单板弹簧的长度l1、各层非对称形单板弹簧长度构成的等比数列的公比q、簧箍的宽度B,就可对不同的柔性体进行变形仿真。
附图说明
图1是采用非对称形板弹簧虚拟模型的建模方法的柔性体变形仿真流程图;
图2是非对称形板弹簧虚拟模型的建模方法流程图;
图3是柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法中,拉力、变形层数与时延时间关系示意图;
图4是非对称形板弹簧虚拟模型的建模方法示意图,(a)是原始状态,(b)是拉伸下的状态,(c)为该建模方法各层与虚拟柔性体表面的对应关系及各层对应边界点的受力等效图。
具体实施方式:
具体实施例1:
一种柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法,其特征在于该建模方法的主要步骤如下:
步骤1对虚拟场景进行初始化;
步骤2当检测到虚拟代理碰撞到虚拟柔性体表面上任何一点时,在给定虚拟接触拉力F作用下,虚拟代理与虚拟柔性体交互的局部区域内部填充非对称形板弹簧虚拟模型,在交互过程中,输出反馈为采用非对称形板弹簧虚拟模型计算出来的反应在拉力作用下虚拟柔性体实时变形仿真的力触觉信息的信号,所述非对称形板弹簧虚拟模型的建模方法为:
(1)参数初始化,
(2)在给定虚拟接触拉力F作用下,当虚拟代理碰撞到虚拟柔性体表面上的任意点时,在碰撞点处下悬挂一叶片宽度为b、厚度为h、长度为l1的非对称形单板弹簧,形成第一层,在第一层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为l2=l1q的非对称形单板弹簧,形成第二层,在第二层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为l3=l1q2的非对称形单板弹簧,形成第三层,依此类推,在第i-1层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为li=l1qi-1的非对称形单板弹簧,形成第i层,i=1,2,3,……N,N为自然数,各层非对称形单板弹簧长度构成以l1为首项,以q为公比的等比数列,q=1~2,用一宽度为B的簧箍,将组成非对称形板弹簧的所有非对称形单板弹簧夹紧在一起,且簧箍可在非对称形板弹簧上移动,b≤B≤b+0.1,
以虚拟接触拉力F的作用线为分界线,两侧非对称形板弹簧的叶片数相等,叶片厚度相同,分界线将第i层的非对称形单板弹簧长度li分成li′、li″左右两部分,且左右两部分长度的比值相同,即:
l1′、l1″,l2′、l2″,l3′、l3″,li′、li″分别表示分界线将第一、二、三、i层的非对称形单板弹簧长度l1、l2、l3、li分成左右两部分,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
采用簧箍夹紧时,组成柔性力触觉再现的非对称形板弹簧虚拟模型中所有非对称形单板弹簧均有一段无效部分,在计算弹簧刚度时,应考虑其影响并进行必要的修正,修正后第i层的非对称形单板弹簧的弹簧刚度为:
其中,li、Pi、B分别为第i层的非对称形单板弹簧的长度、第i层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
第i层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度Pi满足:
Pi=Pi′λ(1+λ) (3)
Pi′为第i层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,其满足:
其中,E、I、li′分别为弹性模量、截面惯性矩、第i层的非对称形单板弹簧长度li以分界线为界左侧的长度,
截面惯性矩I满足:
在虚拟接触拉力F作用下,若柔性体中共有M层的非对称形单板弹簧产生变形,则第M层称为变形截止层,
若给定的虚拟接触拉力F能使第一层的非对称形单板弹簧被拉伸时产生的变形量X1达到第一层的非对称形单板弹簧给定的挠度值这种情况下假定前M-1层的非对称形单板弹簧被拉伸时产生的变形量均与第一层的非对称形单板弹簧给定的挠度值相同,第M层的非对称形单板弹簧被拉伸时产生的变形量不大于第一层的非对称形单板弹簧给定的挠度值,
第一层的非对称形单板弹簧消耗的拉力F1为:
XC1表示与第一层的非对称形单板弹簧的宽度、厚度、长度均相同的对应的对称式单板弹簧给定的挠度值,λ表示每层的非对称形单板弹簧以分界线为界,分成左右两部分长度的比值,且每层分成左右两部分长度的比值均相同,
除第一层和变形截止层第M层外,其余各层的非对称形单板弹簧消耗的拉力Fj为:
j的取值范围为2~M-1,
其中,为修正后第j层的非对称形单板弹簧的弹簧刚度,
lj、Pj、B分别为第j层的非对称形单板弹簧的长度、第j层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
lj=l1qj-1 (11)
第j层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度Pj满足:
Pj=Pj′λ(1+λ) (12)
Pj′表示第j层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,E、I、lj′分别为弹性模量、截面惯性矩、第j层的非对称形单板弹簧长度lj以分界线为界左侧的长度,
lj=l1qj-1 (15)
截面惯性矩I满足:
变形截止层第M层的变形量为:
lM、PM、B分别为变形截止层第M层的非对称形单板弹簧的长度、变形截止层第M层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
lM=l1qM-1 (19)
变形截止层第M层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度PM满足:
PM=PM′λ(1+λ) (20)
PM′表示第M层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,E、I、lM′分别为弹性模量、截面惯性矩、第M层的非对称形单板弹簧长度lM以分界线为界左侧的长度,
lM=l1qM-1 (23)
截面惯性矩I满足:
(3)使虚拟接触拉力F作用于虚拟柔性体碰撞点,第i层对应的非对称形单板弹簧被拉伸,如果前i层所有的非对称形单板弹簧消耗的拉力Fi之和小于给定的虚拟接触拉力F,且前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间满足刷新频率1000Hz以上的要求,设共计经过时延Ti后,第i层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值,只有当第i层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值后,第i+1层对应的非对称形单板弹簧才开始被拉伸,依此类推,直到前M层所有的非对称形单板弹簧消耗的拉力之和不小于给定的虚拟接触拉力,或第M层的非对称形单板弹簧产生拉伸变形总计需要的时延时间不满足刷新频率的要求;
用ti、Ti分别表示第i层的非对称形单板弹簧产生拉伸变形需要的时延时间、前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间,且令层间的时延时间满足以第一层的非对称形单板弹簧产生拉伸变形需要的时延时间t1为首项、以w为公比的等比数列:
ti=wi-1t1 (25)
从虚拟代理碰撞接触到虚拟柔性体表面算起,假定前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间Ti满足Ti<T,其中
T为力触觉再现刷新频率的倒数;
所述非对称形板弹簧虚拟模型中每层的非对称形单板弹簧变形量之和的叠加对外等效为虚拟柔性体表面的变形,
其中,X为前M层所有的非对称形单板弹簧产生的拉伸变形量之和,Xi为前M-1层中任意第i层的非对称形单板弹簧产生的拉伸变形量,XM为变形截止层第M层产生的拉伸变形量,F为给定的虚拟接触拉力,为修正后第一层的非对称形单板弹簧刚度。
具体实施例2:
1、构建虚拟医用镊子模型和虚拟肾模型,实现虚拟场景的初始化。
本实例中所有虚拟医用镊子和虚拟肾模型都直接采用从3DS MAX 9.0软件中导出的OBJ格式,以440个质点,869个三角网格构成的虚拟医用镊子和2205个质点,4712个三角网格构成的虚拟肾模型为例来进行变形仿真,实验过程中模型获取和修改非常方便;操作***为Windows 2000,以3DS MAX 9.0、OpenGL图形库为基础,在VC++6.0软件开发平台上进行仿真。
2、当检测到虚拟医用镊子碰撞到虚拟肾表面上任何一点时,在给定虚拟接触拉力F=1.5N作用下,虚拟医用镊子与虚拟肾交互的局部区域内部填充非对称形板弹簧虚拟模型,在交互过程中,输出反馈为采用非对称形板弹簧虚拟模型计算出来的反应在拉力作用下虚拟肾实时变形仿真的力触觉信息的信号,
假定非对称形板弹簧虚拟模型中,叶片宽度b=2mm、厚度h=1.4mm、各层非对称形单板弹簧长度构成以l1=3mm为首项,以q=1.05为公比的等比数列,用一宽度为B=2.05mm的簧箍,将组成非对称形板弹簧的所有非对称形单板弹簧夹紧在一起,且簧箍可在非对称形板弹簧上移动,计算的中间过程、最后的数据均按四舍五入法保留小数点后6位。
任一层的非对称形单板弹簧的材质均相同,即取E=4.8204×105Pa,
因任一层的非对称形单板弹簧的叶片宽度、厚度均相同,故任一层的非对称形单板弹簧的截面惯性矩I为:
以虚拟接触拉力F的作用线为分界线,两侧非对称形板弹簧的叶片数相等,叶片厚度相同,每层的非对称形单板弹簧以分界线为界,被分成左右两部分,且每层的非对称形单板弹簧被分成左右两部分长度的比值λ=0.6均相同,
在虚拟接触拉力F作用下,若柔性体中共有M层的非对称形单板弹簧产生变形,则第M层称为变形截止层,
与第一层的非对称形单板弹簧的宽度、厚度、长度均相同的对应的对称式单板弹簧给定的挠度值XC1=0.76mm,
若给定的虚拟接触拉力F能使第一层的非对称形单板弹簧被拉伸时产生的变形量X1达到第一层的非对称形单板弹簧给定的挠度值这种情况下假定前M-1层的非对称形单板弹簧被拉伸时产生的变形量均与第一层的非对称形单板弹簧给定的挠度值相同,第M层的非对称形单板弹簧被拉伸时产生的变形量不大于第一层的非对称形单板弹簧给定的挠度值,
第一层的非对称形单板弹簧长度l1以分界线为界左侧的长度l1′为:
第一层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度P1′满足:
第一层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度P1满足:
P1=P1′λ(1+λ)=0.464493×0.6×1.6=0.445913N/mm (5)
采用簧箍夹紧时,组成柔性力触觉再现的非对称形板弹簧虚拟模型中所有非对称形单板弹簧均有一段无效部分,在计算弹簧刚度时,应考虑其影响并进行必要的修正,修正后第一层的非对称形单板弹簧的弹簧刚度为:
若给定的虚拟接触拉力F能使第一层的非对称形单板弹簧被拉伸时产生的变形量与第一层的非对称形单板弹簧给定的挠度值相同,则第一层的非对称形单板弹簧消耗的拉力F1为:
F1=0.674470N<F=1.5N,
假定第一层的非对称形单板弹簧产生拉伸变形需要的时延时间t1=10-5s;
第一层的非对称形单板弹簧产生拉伸变形需要的时延时间T1=t1=10-5s<T;
因此,第一层的非对称形单板弹簧消耗的拉力F1=0.674470N<F=2.6N,且第一层的非对称形单板弹簧产生拉伸变形需要的时延时间满足力触觉再现刷新频率1000Hz以上的要求;只有当第一层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值后,第二层对应的非对称形单板弹簧才开始被拉伸。
第二层的非对称形单板弹簧长度l2为:
l2=l1q=3×1.05=3.15mm (8)
第二层的非对称形单板弹簧长度l2以分界线为界左侧的长度l2′为:
第二层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度P2′满足:
第二层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度P2满足:
P2=P2′λ(1+λ)=0.401247×0.6×1.6=0.385197N/mm (11)
若给定的虚拟接触拉力F能使第二层的非对称形单板弹簧被拉伸时产生的变形量与第一层的非对称形单板弹簧给定的挠度值相同,则第二层的非对称形单板弹簧消耗的拉力F2为:
前二层的非对称形单板弹簧消耗的拉力之和为:
F1+F2=0.674470+0.561687=1.236157N<F=1.5N
层间的时延时间满足以第一层的非对称形单板弹簧产生拉伸变形需要的时延时间t1为首项,以w=1.1为公比的等比数列,
前两层的非对称形单板弹簧产生拉伸变形总计需要的时延时间:T2=t1+t2=(1+w)t1=2.1×10-5s<T,这里T为力触觉再现刷新频率的倒数,
因此,前两层的非对称形单板弹簧消耗的拉力之和小于给定的虚拟接触拉力,且前两层的非对称形单板弹簧产生拉伸变形总计需要的时延时间满足刷新频率1000Hz以上的要求;只有当第二层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值后,第三层对应的非对称形单板弹簧才开始被拉伸。
第三层的非对称形单板弹簧长度l3为:
l3=l1q2=3×1.052=3.3075mm (14)
第三层的非对称形单板弹簧长度l3以分界线为界左侧的长度l3′为:
第三层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度P3′满足:
第三层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度P3满足:
P3=P3′λ(1+λ)=0.346612×0.6×1.6=0.332748N/mm (17)
若给定的虚拟接触拉力F能使第三层的非对称形单板弹簧被拉伸时产生的变形量与第一层的非对称形单板弹簧给定的挠度值相同,则第三层的非对称形单板弹簧消耗的拉力F3为:
前三层的非对称形单板弹簧消耗的拉力之和为:
F1+F2+F3=0.674470+0.561687+0.468772=1.704929N>F=1.5N
由于前三层的非对称形单板弹簧消耗的拉力之和不小于给定的虚拟接触拉力,则第三层为变形截止层,不需要再判断是否满足刷新频率的要求,此时前二层的非对称形单板弹簧被拉伸时产生的变形量与第一层的非对称形单板弹簧给定的挠度值相同,变形截止层第三层的非对称形单板弹簧被拉伸时产生的变形量为:
故在给定的虚拟接触拉力F=2.6N作用下,柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法中,前三层的非对称形单板弹簧被拉伸时产生的变形量之和的叠加对外等效为虚拟柔性体表面的变形,变形量之和为:
X=X1+X2+X3
=XC1+XC1+X3 (21)
=0.76+0.76+0.427758
=1.947758mm
注意:在采用柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法来计算在给定虚拟接触拉力作用下柔性体实时变形仿真的过程中,若b、h、l1、q、B这些参数选取的过大,则柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法中变形的层数就少,计算量小,实时性好,但变形仿真效果不佳;如果b、h、l1、q、B这些参数选取的过小,则柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法中变形的层数就越多,计算量大,实时性不佳,但变形仿真效果较好;另外在设置t1和ti之间的比例关系时,要考虑到程序运行时计算机本身的硬件配置,故在调试整个程序的过程中,要折中选择这些参数,不断反复调试,从而使变形效果更加逼真。
为验证本发明的实施效果,操作者通过PHANTOM OMNI手控器端部的手柄来触摸、感知和控制虚拟医用镊子对虚拟肾进行拉拽的变形仿真,并将交互过程中产生的力触觉信息实时反馈给操作者。实验结果表明:该建模方法是有效的,在交互过程中,操作者可实时对虚拟肾进行拉拽操作,真实地感知到变形仿真过程中虚拟医用镊子与虚拟肾之间的自然交互操作。
Claims (1)
1.一种柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法,其特征在于该建模方法的主要步骤如下:
步骤1对虚拟场景进行初始化;
步骤2当检测到虚拟代理碰撞到虚拟柔性体表面上任何一点时,在给定虚拟接触拉力F作用下,虚拟代理与虚拟柔性体交互的局部区域内部填充非对称形板弹簧虚拟模型,在交互过程中,输出反馈为采用非对称形板弹簧虚拟模型计算出来的反应在拉力作用下虚拟柔性体实时变形仿真的力触觉信息的信号,所述非对称形板弹簧虚拟模型的建模方法为:
(1)参数初始化,
(2)在给定虚拟接触拉力F作用下,当虚拟代理碰撞到虚拟柔性体表面上的任意点时,在碰撞点处下悬挂一叶片宽度为b、厚度为h、长度为l1的非对称形单板弹簧,形成第一层,在第一层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为l2=l1q的非对称形单板弹簧,形成第二层,在第二层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为l3=l1q2的非对称形单板弹簧,形成第三层,依此类推,在第i-1层的非对称形单板弹簧下,设置一叶片宽度为b、厚度为h、长度为li=l1qi-1的非对称形单板弹簧,形成第i层,i=1,2,3,……N,N为自然数,各层非对称形单板弹簧长度构成以l1为首项,以q为公比的等比数列,q=1~2,用一宽度为B的簧箍,将组成非对称形板弹簧的所有非对称形单板弹簧夹紧在一起,且簧箍可在非对称形板弹簧上移动,b≤B≤b+0.1,
以虚拟接触拉力F的作用线为分界线,两侧非对称形板弹簧的叶片数相等,叶片厚度相同,分界线将第i层的非对称形单板弹簧长度li分成li′、li″左右两部分,且左右两部分长度的比值相同,即:
l1′、l1″,l2′、l2″,l3′、l3″,li′、li″分别表示分界线将第一、二、三、i层的非对称形单板弹簧长度l1、l2、l3、li分成左右两部分,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,li、Pi、B分别为第i层的非对称形单板弹簧的长度、第i层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
第i层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度Pi满足满足:
Pi=Pi′λ(1+λ) (3)
Pi′为第i层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,其满足:
其中,E、I、li′分别为弹性模量、截面惯性矩、第i层的非对称形单板弹簧长度li以分界线为界左侧的长度,
截面惯性矩I满足:
在虚拟接触拉力F作用下,若虚拟柔性体中共有M层的非对称形单板弹簧产生变形,则第M层称为变形截止层,
若给定的虚拟接触拉力F能使第一层的非对称形单板弹簧被拉伸时产生的变形量X1达到第一层的非对称形单板弹簧给定的挠度值这种情况下假定前M-1层的非对称形单板弹簧被拉伸时产生的变形量均与第一层的非对称形单板弹簧给定的挠度值相同,第M层的非对称形单板弹簧被拉伸时产生的变形量不大于第一层的非对称形单板弹簧给定的挠度值,
第一层的非对称形单板弹簧消耗的拉力F1为:
XC1表示与第一层的非对称形单板弹簧的宽度、厚度、长度均相同的对应的对称式单板弹簧给定的挠度值,λ表示每层的非对称形单板弹簧以分界线为界,分成左右两部分长度的比值,且每层分成左右两部分长度的比值均相同,
除第一层和变形截止层第M层外,其余各层的非对称形单板弹簧消耗的拉力Fj为:
j的取值范围为2~M-1,
lj、Pj、B分别为第j层的非对称形单板弹簧的长度、第j层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
lj=l1qj-1 (11)
第j层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度Pj满足:
Pj=Pj′λ(1+λ) (12)
Pj′表示第j层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,E、I、lj′分别为弹性模量、截面惯性矩、第j层的非对称形单板弹簧长度lj以分界线为界左侧的长度,
lj=l1qj-1 (15)
截面惯性矩I满足:
变形截止层第M层的变形量为:
lM、PM、B分别为变形截止层第M层的非对称形单板弹簧的长度、变形截止层第M层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度、簧箍宽度,
lM=l1qM-1 (19)
变形截止层第M层的非对称形单板弹簧在未考虑无效部分时的弹簧刚度PM满足:
PM=PM′λ(1+λ) (20)
PM′表示第M层的非对称形单板弹簧在未考虑无效部分时左侧的弹簧刚度,λ表示每层的非对称形单板弹簧以分界线为界分成左右两部分的长度的比值,且每层分成左右两部分长度的比值均相同,
其中,E、I、lM′分别为弹性模量、截面惯性矩、第M层的非对称形单板弹簧长度lM以分界线为界左侧的长度,
lM=l1qM-1 (23)
截面惯性矩I满足:
(3)使虚拟接触拉力F作用于虚拟柔性体碰撞点,第i层对应的非对称形单板弹簧被拉伸,如果前i层所有的非对称形单板弹簧消耗的拉力Fi之和小于给定的虚拟接触拉力F,且前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间满足刷新频率1000Hz以上的要求,设共计经过时延Ti后,第i层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值,只有当第i层的非对称形单板弹簧被拉伸到第一层的非对称形单板弹簧给定的挠度值后,第i+1层对应的非对称形单板弹簧才开始被拉伸,依此类推,直到前M层所有的非对称形单板弹簧消耗的拉力之和不小于给定的虚拟接触拉力,或第M层的非对称形单板弹簧产生拉伸变形总计需要的时延时间不满足刷新频率的要求;
用ti、Ti分别表示第i层的非对称形单板弹簧产生拉伸变形需要的时延时间、前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间,且令层间的时延时间满足以第一层的非对称形单板弹簧产生拉伸变形需要的时延时间t1为首项、以w为公比的等比数列:
ti=wi-1t1 (25)
从虚拟代理碰撞接触到虚拟柔性体表面算起,假定前i层所有的非对称形单板弹簧产生拉伸变形总计需要的时延时间Ti满足Ti<T,其中
T为力触觉再现刷新频率的倒数;
所述非对称形板弹簧虚拟模型中每层的非对称形单板弹簧变形量之和的叠加对外等效为虚拟柔性体表面的变形,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102991871A CN101964023B (zh) | 2010-09-27 | 2010-09-27 | 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102991871A CN101964023B (zh) | 2010-09-27 | 2010-09-27 | 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101964023A CN101964023A (zh) | 2011-02-02 |
CN101964023B true CN101964023B (zh) | 2012-02-22 |
Family
ID=43516894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102991871A Expired - Fee Related CN101964023B (zh) | 2010-09-27 | 2010-09-27 | 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101964023B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101996280B (zh) * | 2010-09-27 | 2012-07-11 | 南京信息工程大学 | 用于力触觉再现的板弹簧虚拟模型的建模方法 |
CN103699724B (zh) * | 2013-12-10 | 2016-07-06 | 南京信息工程大学 | 一种模拟柔性体扭转变形的蜗卷弹簧模型的方法 |
CN103699741B (zh) * | 2013-12-23 | 2016-08-17 | 南京信息工程大学 | 一种增强模拟柔性体可旋转变形的发条弹簧模型 |
CN103714205B (zh) * | 2013-12-23 | 2016-10-05 | 南京信息工程大学 | 一种模拟柔性体可旋转变形的发条弹簧模型 |
CN103699753B (zh) * | 2013-12-31 | 2016-08-17 | 南京信息工程大学 | 一种支持增强力触觉反馈的多股螺旋弹簧虚拟模型 |
CN104615879A (zh) * | 2015-01-28 | 2015-05-13 | 南京信息工程大学 | 一种增强力触觉再现的圆柱扭转弹簧虚拟模型 |
CN106295085A (zh) * | 2016-10-17 | 2017-01-04 | 南京信息工程大学 | 柔性力触觉再现的叠合碟形弹簧虚拟模型的建模方法 |
CN106560821A (zh) * | 2016-10-17 | 2017-04-12 | 南京信息工程大学 | 基于对合碟形弹簧的柔性体力触觉建模方法 |
CN106528993A (zh) * | 2016-11-03 | 2017-03-22 | 南京信息工程大学 | 基于由碟形弹簧片构成的组合弹簧虚拟模型建模方法 |
CN107978367B (zh) * | 2017-11-28 | 2021-07-13 | 电子科技大学 | 一种基于虚拟现实的远程诊疗机器人力触觉再现方法 |
CN111981997B (zh) * | 2020-07-21 | 2022-01-28 | 广东工业大学 | 一种基于大变形的ps-oct应变估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2905497B1 (fr) * | 2006-09-01 | 2008-11-14 | Michelin Soc Tech | Procede de simulation du comportement thermo-mecanique d'un pneu, et application |
CN101587372B (zh) * | 2009-07-03 | 2010-09-15 | 东南大学 | 一种用于虚拟现实人机交互的增强力触觉建模方法 |
-
2010
- 2010-09-27 CN CN2010102991871A patent/CN101964023B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101964023A (zh) | 2011-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101964023B (zh) | 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 | |
CN101976299B (zh) | 柔性力触觉再现的对称式板弹簧虚拟模型的建模方法 | |
CN101587372B (zh) | 一种用于虚拟现实人机交互的增强力触觉建模方法 | |
CN101976298B (zh) | 增强力触觉再现的对称式板弹簧虚拟模型的建模方法 | |
Gao et al. | Haptic B-spline surface sculpting with a shaped tool of implicit surface | |
CN101853072B (zh) | 一种用于柔性体变形仿真的力触觉建模方法 | |
CN101996280B (zh) | 用于力触觉再现的板弹簧虚拟模型的建模方法 | |
CN103869983A (zh) | 一种用于力触觉人机交互的柔性物体变形仿真方法 | |
CN107992672A (zh) | 一种基于虚拟弹簧的软组织形变建模方法 | |
CN110289104A (zh) | 软组织按压和形变恢复的模拟方法 | |
CN103745058A (zh) | 一种任意形状软组织表皮上受拉力/变形的模拟方法 | |
CN106528993A (zh) | 基于由碟形弹簧片构成的组合弹簧虚拟模型建模方法 | |
CN101901565B (zh) | 支持力触觉反馈的虚拟柔性体变形手术仿真*** | |
CN116631633A (zh) | 一种软组织拉伸撕裂模拟***、方法、设备、存储介质 | |
CN106295085A (zh) | 柔性力触觉再现的叠合碟形弹簧虚拟模型的建模方法 | |
CN106485028B (zh) | 支持实时按压的柔性体变形仿真建模方法 | |
Yang et al. | Haptic function evaluation of multi-material part design | |
Zhang et al. | Layered rhombus-chain-connected model for real-time haptic rendering | |
Salsedo et al. | Architectural design of the haptex system | |
CN103714205B (zh) | 一种模拟柔性体可旋转变形的发条弹簧模型 | |
Liu et al. | Modelling and simulation of vascular tissue based on finite element method | |
CN111103973A (zh) | 模型处理方法、装置、计算机设备和存储介质 | |
Wang et al. | A haptics–based virtual simulation system for product design | |
CN103678806B (zh) | 一种模拟柔性体扭转变形的盘香弹簧模型 | |
Leidinger et al. | Explicit isogeometric b-rep analysis on trimmed nurbs-based multi-patch cad models in ls-dyna |
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: 20120222 Termination date: 20140927 |
|
EXPY | Termination of patent right or utility model |