专利名称：：一种用于虚拟现实人机交互的增强力触觉建模方法
技术领域：
：本发明涉及一种增强力触觉建模方法，尤其涉及一种用于虚拟现实人机交互的基于物理意义的软组织变形仿真力触觉建模方法。
背景技术：
：当前，对于虚拟现实技术而言，不断优化和完善的图形算法及其软件巳经能够惟妙惟肖地模拟三维虛拟场景，而人们对"可望而不可及"的不满又促使近年来利用力触觉交互设备在计算机所营造的虚拟世界中模拟软组织在外力作用下产生变形得到越来越多研究人员的关注与重视。其中，针对虚拟现实人机交互的基于物理意思的软组织变形仿真的力触觉建模方法已成为当前虚拟现实技术研究的前沿和热点。目前基于物理意义的软组织变形仿真力触觉建模方法主要有5种①弹簧-质点建模方法具有简单易行、计算量较小、易于实现等优点，但精度和稳定性有限，另外还存在拓扑分析困难等弊端。②有限元建模方法虽然较为准确，但涉及大量繁杂的计算，实时交互性能差。③Kim等人提出的形状保持链建模方法(shaperetainingchainlinkedmodel,S-ChainModel)具有实时性较好的优点，但链元素是以刚性体为假设的，精度不高。④边界元建模方法虽对建模方法的边界进行离散，简化了计算，但在稳定性方面却存在一定的不足。⑤长单元建模方法具有较高的触觉反馈刷新速度，求解方便，但由于建模方法抽象度较高，因此计算精度较低。以上说明，这些常用的软组织变形物理建模方法均存在计算较为繁杂和仿真精度不高等问题，继而影响了计算的实时性和有效性。因此，在保证变形建模方法准确性的同时，简化算法，提高计算的实时性是当前虚拟力触觉交互系统亟待解决的首要问题。人机工程学的最新研究表明人的力触觉感知存在着信号与噪声之间的协同行为，即适量噪声的加入会提高人的力触觉感知能力，这主要原因是因为人手皮5肤对一个持续作用力具有很强的适应性，而对噪声信号比较敏感，所以当加入合适的噪声信号，会产生比较真实的感觉。根据这一原理，为了提高虚拟现实人机力触觉交互过程的稳定性和感知真实性，提出了一种基于叠加噪声的增强力触觉建模方法。
发明内容本发明提出一种基于物理意义的叠加噪声的增强力触觉建模方法，并将其用于虚拟现实人机交互的软组织变形仿真。该建模方法能准确快速的计算变形量，实现对软组织的实时变形仿真，从而提高虚拟力触觉交互的逼真度。本发明采用如下技术方案一种叠加噪声的增强力触觉建模方法，该建模方法的建模方法其主要步骤如下步骤1对虚拟场景进行初始化；步骤2虚拟代理(如虚拟镊子、虚拟手指、虚拟探针等)在碰撞到虚拟软组织之前，即在靠近虚拟软组织的过程中，反馈输出与软组织表面距离有关的服从iV(O,(《_^_)2)分布的白噪声信号，该白噪声信号为(Tx=cr.^-2log(l-ra"c/()/32768)cos(2;r.層d()/32768)其中o"为白噪声的标准差，log()为自然对数，mwd()为0-32767间的一个随机整数，K为正比例系数，d为虚拟代理与虚拟软组织表面之间的距离，^为距离阈值，c/。>0;由式(l)可知，随着虚拟代理与虚拟软组织表面的距离越来越小，白噪声的标准差由小到大，逐渐递增。在此过程中，仅输出反馈与软组织表面距离有关的白噪声信号。步骤3当检测到虚拟代理碰撞到虚拟软组织表面上任何一点时，在给定虚拟接触力F作用下，虚拟代理与虚拟软组织交互的局部区域内部填充该等比例层6状并联力触觉建模方法，在交互过程中，反馈输出为:服从TV(O,(《.;)2)恒定的白噪声信号与采用等比例层状并联力触觉建模方法计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息线性叠加在一起的信号，所述等比例层状并联力触觉的建模方法为(1)参数初始化，(2)当虚拟代理碰撞到软组织表面上的任意点时，在碰撞点处下方悬挂一弹簧，形成第一层；在第一层4"个弹簧下悬挂十字杆，在十字杆的各个端部下各悬挂一弹簧，形成由^个枏同弹簧组成的第二层；在第二层的各个弹簧下，分别悬挂十字杆，每根十字杆的各个端部下各悬挂一弹簧，形成由42个相同弹簧组成的第三层；依此类推，在第/-1层的各个弹簧下，分别悬挂十字杆，每根十字杆的各个端部下各悬挂一弹簧，形成由4"'个相同弹簧组成的第f层，f=l，2，3，，见iV为自然数，处于同一层的弹簧系数及变形量均相同，相邻层的弹簧弹性系数之比为常数，即第/层任一弹簧的弹性系数^与第一层弹簧的弹性系数A呈如下式所示的指数倍关系!=一！.=1,2,3,……，iV(2)在这里令"为常数；假定第/层所有弹簧的原始长度、拉伸及压縮后的极限长度、弹簧的变形量分别为X,.。、JT,.max、Xmin、AZ,，它们均以第一层弹簧的对应量为首项，组成以《为公比的等比数列，k,-《。k—1弹簧被拉伸时/=2,3,......W(3)《腿-^。k'—1弹簧被压缩时所述的十字杆由长度相同的第一杆和第二杆组成，第一杆与第二杆垂直且其交点与第一杆和第二杆的质心相重合，(3)使虚拟接触力作用于软组织碰撞点，第/层对应的4"个弹簧同时开始被拉伸或压縮，如果前/层所有弹簧的变形量与对应的弹性系数的乘积之和小于给定的虚拟接触力，且第/层弹簧变形总计需要的时延时间满足刷新频1000Hz以上的要求，设共计经过时延7;,第/层的所有弹簧均被拉伸或压縮到极限长度，只有当第/层的所有弹簧均被拉伸或压縮到极限长度后，第/+1层对应的4'个弹簧才同时开始被拉伸或压縮，依此类推，直到前M层所有弹簧的变形量与对应的弹性系数的乘积之和不小于给定的虚拟接触力，或第M层弹簧变形总计需要的时延时间&不满足刷新频率的要求，则第M层为变形截止层，此时前M-1层的所有弹簧均被拉伸或压縮到极限长度，第M层的所有弹簧被拉伸或压縮的长度不大于极限长度，第M层任一弹簧对应的变形量为'仏l^-弹簧被拉伸时尸"=乂(4)~^-弹簧被压縮时其中A^^为第M层任一弹簧对应的变形量，F为虚拟接触力，A,.为第/层任一弹簧的弹性系数，《_为第一层弹簧被拉伸后的极限长度，《。为第一层弹簧的原始长度，1,_为第一层弹簧被压縮后的极限长度，《为公比系数，/^为第M层弹簧的个数，《为第Z'层弹簧的个数为；研究表明，在虚拟现实人机交互过程中，力触觉再现和视觉显示相比，对计算速度和刷新频率有更高的要求，实时的力计算是保证力触觉感觉平稳，图像流畅，无停顿感，模拟效果逼真的前提，要使其满足交互的实时性、稳定性和逼真度的综合要求，力触觉再现的满足刷新频率1000Hz以上的要求，即要求变形计算具有快速性的特点，用f,、7;分别表示第z'层的弹簧变形需要的时间、前/层的所有弹簧变形需要的时间，且令层间的时延时间满足,,=4'—、(5)8从虚拟代理碰撞接触到虚拟软组织表面算起，假定第/层的弹簧变形总计需要时延的时间7;,且其须满足i;〈r，其中;=+^+^+.+4"-l(6)r为力触觉再现刷新频率的倒数；等比例层状并联力触觉建模方法中每层各个弹簧变形量之和的叠加对外等效为物体表面的变形，与之相连的弹簧弹性力的合力等效为物体表面接触力。本发明的优点(1)和以往常用的基于物理意义的软组织变形仿真力触觉建模方法相比，该建模方法中第一层弹簧对应的质点作为与虚拟软组织表面直接交互的作用点的中心O，第二层对应的4个弹簧质点，在虚拟软组织表面构成以O为中心，r为半径的圆，第三层对应的16个弹簧质点，只有最外层的4个弹簧质点参与软组织的变形计算，其余12个弹簧质点仅起到分担力的作用，在虚拟软组织表面构成以o为中心，以r+丄r为半径的圆。依次类推，第/层对应的4"个弹簧质点，4只有最外层的4个弹簧质点参与软组织的变形计算，其余4"-4个弹簧质点仅起到分担力的作用，在虚拟软组织表面构成以o为中心，以r+丄(/-2)r为半径的圆，4整个建模方法构成一系列的同心圆，简化计算，加快了变形计算的速度。(2)等比例层状并联力触觉建模方法用于在交互过程中，虚拟代理与虚拟软组织交互的局部区域软组织的变形计算，在此过程中实时反馈给操作者的除了软组织实时变形仿真的力触觉信息的信号外，还有线性叠加上的人手皮肤对一个持续作用力比较敏感的白噪声信号，增强操作者的感知，在交互过程中，自然舒适、力触觉感觉平稳、模拟效果逼真。G)通过改变建模方法中第一层弹簧的原始长度、各层弹簧的公比系数g及各层对应同心圆的半径，就可对不同的软组织进行变形仿真，从而实现用一种建模方法对不同软组织进行仿真，适用性广。(4)可将其应用于虚拟外科手术仿真、远程医疗、遥操作机器人控制等领域。9图1是软组织变形仿真流程图2是基于叠加噪声的增强力触觉建模方法示意图；图3是虚拟代理与软组织交互过程中的增强力触觉建模方法流程图；图4是增强力触觉建模方法中，外力、变形层数与时延时间关系示意图；图5是等比例层状并联力触觉计算建模方法示意图，(a)是原始状态，(b)是拉伸后的状态(c)是压縮后的状态(d)为该建模方法各层与虚拟软组织表面的对应关系及各层对应边界点的受力等效图6镊子与肝脏交互过程渲染效果图，(a)是虚拟肝脏受压后受力变形渲染效果图，(b)是虚拟肝脏受拉后受力变形渲染效果图。具体实施方式-实施例1:一种用于虚拟现实人机交互的增强力触觉建模方法，其特征在于该建模方法的主要步骤如下步骤1对虚拟场景进行初始化；步骤2虚拟代理(如虚拟镊子、虚拟手指、虚拟探针等)在碰撞到虚拟软组织之前，即在靠近虚拟软组织的过程中，反馈输出与软组织表面距离有关的服从W(O,(^~^)2)分布的白噪声信号，该白噪声信号为cr':c=cr.V-2log(l—()/32768).cos(2;r.mwd()/32768)o"=AT--其中(T为白噪声的标准差，Iog()为自然对数，为0-32767间的一个随机整数，《为正比例系数，"为虚拟代理与虚拟软组织表面之间的距离，^为距离阈值，t/。>0;歩骤3当检测到虚拟代理碰撞到虚拟软组织表面上任何一点时，在给定虚拟接触力F作用下，虚拟代理与虚拟软组织交互的局部区域内部填充该等比例层10状并联力触觉建模方法，在交互过程中，反馈输出为:服从iV(O,(《.±)2)恒定的白噪声信号与采用等比例层状并联力触觉建模方法计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息线性叠加在一起的信号，所述等比例层状并联力触觉的建模方法为-(1)参数初始化，(2)当虚拟代理碰撞到软组织表面上的任意点时，在碰撞点处下方悬挂一弹簧，形成第一层；在第一层^个弹簧下悬挂十字杆，在十字杆的各个端部下各悬挂一弹簧，形成由W个相同弹簧组成的第二层；在第二层的各个弹簧下，分别悬挂十字杆，每根十字杆的各个端部下各悬挂一弹簧，形成由42个相同弹簧组成的第三层；依此类推，在第层的各个弹簧下，分别悬挂十字杆，每根十字杆的各个端部下各悬挂一弹簧，形成由4'—'个相同弹簧组成的第/层，/=1，2，3，，7V，W为自然数，处于同一层的弹簧系数及变形量均相同，相邻层的弹簧弹性系数之比为常数，即第f层任一弹簧的弹性系数^与第一层弹簧的弹性系数^呈如下式所示的指数倍关系5=^/=2,3,……，W(2)在这里令/为常数；假定第z'层所有弹簧的原始长度、拉伸及压縮后的极限长度、弹簧的变形量分别为《。、X,.max、Xmin、A《，它们均以第一层弹簧的对应量为首项，组成以《为公比的等比数列，'k'隨-^。k'—'弹簧被拉伸时线=《'誦—《。/=2,3，……邪)Ilmn-^。|《'—1弹簧被压縮时所述的十字杆由长度相同的第一杆和第二杆组成，第一杆与第二杆垂直且其交点与第一杆和第二杆的质心相重合，11(3)使虛拟接触力作用于软组织碰撞点，第Z层对应的4'-'个弹簧同时开始被拉伸或压縮，如果前/层所有弹簧的变形量与对应的弹性系数的乘积之和小于给定的虚拟接触力，且第/层弹簧变形总计需要的时延时间满足满足刷新频率1000Hz以上的要求，设共计经过时延7;,第/层的所有弹簧均被拉伸或压縮到极限长度，只有当第z'层的所有弹簧均被拉伸或压縮到极限长度后，第/+1层对应的4'个弹簧才同时开始被拉伸或压缩，依此类推，直到前M层所有弹簧的变形量与对应的弹性系数的乘积之和不小于给定的虚拟接触力，或第M层弹簧变形总计需要的时延时间&不满足刷新频率的要求，则第M层为变形截止层，此时前M-l层的所有弹簧均被拉伸或压縮到极限长度，第M层的所有弹簧被拉伸或压縮的长度不大于极限长度，第M层任一弹簧对应的变形量为尸-Zh(A鹏—不。)IV—、S弹簧被拉伸时A/AM(4)/=1弹簧被压縮时尸为虚拟接触力，A,为第/层其中AZW为第M层任一弹簧对应的变形邐任一弹簧的弹性系数，Z^为第一层弹簧被拉伸后的极限长度，《。为第一层弹簧的原始长度，J^,为第一层弹簧被压縮后的极限长度，《为公比系数，P^为第M层弹簧的个数，《为第i层弹簧的个数为；用A、i;分别表示第z'层的弹簧变形需要的时间、前z'层的所有弹簧变形需要的时间，且令层间的时延时间满足<formula>formulaseeoriginaldocumentpage12</formula>(5)从虚拟代理碰撞接触到虚拟软组织表面算起，假定第/层的弹簧变形总计需要时延的时间7;，且其须满足7]<2%其中<formula>formulaseeoriginaldocumentpage12</formula>(6)r为力触觉再现刷新频率的倒数；等比例层状并联力触觉建模方法中每层各个弹簧变形量之和的叠加对外等效为物体表面的变形，与之相连的弹簧弹性力的合力等效为物体表面接触力。具体实施例2:*1、构建虚拟肝脏模型和虛拟镊子模型，实现虚拟场景的初始化。本实例中所有虚拟肝脏和虚拟镊子模型都直接采用从3DSMAX9.0软件中导出OBJ的格式，以3116个质点，6228个三角网格构成的虚拟肝脏和174个质点，348个三角网格构成的虚拟医用镊子为例来进行变形仿真，模型获取和修改非常方便；运行平台为Windows2000，开发语言为VC++6.0;以3DSMAX9.0、OpenGL图形库为基础。2、产生与肝脏表面距离有关的服从iV(O，(《^^)2)分布的白噪声信号。VC++6.0环境下的ram/()产生0-32767间的一个随机整数，令正比例系数《=0.5l%m(这里N表示力的单位牛顿)，虚拟镊子与虚拟肝脏表面之间的初始距离"=200，距离阈值^^20an，在虚拟镊子未碰撞到虚拟肝脏之前，产生的白噪声的标准差C7的大小取决与虚拟镊子与虚拟肝脏之间不断变化的距离。3、当检测到虚拟镊子碰撞到虚拟肝脏表面上任何一点时，假定在给定虚拟接触力F-8N的作用下，虚拟镊子与虚拟肝脏交互的局部区域内部填充该等比例层状并联力触觉建模方法，在交互过程中，反馈输出为:服从iV(O,(《>;)2)分布的恒定的白噪声信号与采用等比例层状并联力触觉建模方法计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息线性叠加在一起的信号，对人体肝脏来说一般取/=2.0，A:严0.15N/cm，取e二2.71828，计算的最后数据按四舍五入法保留小数点后5位；假定第一层弹簧的原始长度、拉伸及压縮后的极限长度、弹簧的变形量分别为Z)。-lcm、Ilmax=l.lcm、Xlmm=0.9，、AZ1=0.1cw，第Z层所有弹簧的原始长度、拉伸及压缩后的极限长度、弹簧的变形量分别为X,。、I,max、X,mm、AX,，它们均以第一层弹簧的对应量为首项，组成以9=1.2为公比的等比数列；若在给定虚拟接触力的作用下，第一层的一个弹簧能被拉伸或压縮到极限长度的话，第一层的一个弹簧对应的变形量为^隨-A。I《1—1=0.k柳弹簧被拉伸时ilmm-A。k'—1=o.icw弹簧被压缩时(第一层一个弹簧的变形量)x(第一层弹簧对应的弹性系数)f0.1cmx0.15N/cm=0.015N弹簧被拉伸时廿，^^入^占桐秘^^「。、TH，《处、^「^"，其小于给定虚拟接触力^-8N;l0.kwx0.15N/cw=0.015N弹簧被压縮时假定第一层的弹簧需要变形的时间为A=10、；假定力触觉再现刷新频率为1200Hz，则力触觉再现刷新频率的倒数7"=-S;1200第一层的弹簧需要变形的时间为7;=6=itr、<r(这里r为力触觉在新刷新频率的倒数)，第一层弹簧的变形量与对应的弹性系数的乘积之和小于给定的虚拟接触力，且第一层弹簧变形总计需要的时延时间满足满足刷新频率1000Hz以上的要求，设共计经过时延《，第一层的所有弹簧均被拉伸或压缩到极限长度，只有当第一层的所有弹簧均被拉伸或压縮到极限长度后，第二层对应的^个弹簧才同时开始被拉伸或压縮。第二层弹簧对应的弹性系数为&=0.15N/cmxe^=0.24731N/on，若在给定虚拟接触力的作用下，第二层的4个弹簧均被拉伸或压縮到极限长度的话，第二层每个弹簧对应的变形量为AX，=^鹏-Zjf=0.1cmx1.2"=0.12麵弹簧被拉伸时《，—=0-lcwx1'2"=0.12cm弹簧被压縮时(第二层每个弹簧的变形量)x(第二层弹簧对应的弹性系数)x(第二层]站敏八wf0.12，x0.24731N/，x4=0.11871N弹簧被拉伸时弹簧个数)叫《处、，「,，』，刖两层所有0.12ox0.24731N/ax4=0.11871N弹簧被压缩时14弹簧的变形量与对应弹簧弹性系数的乘积之和f0.015N+0.11871N=0.13371N弹簧被拉伸时工、八站旌协+、他i『^^+仏,4w愁、+「力^，无论弹簧被拉伸或压缩时均小0.015N+0.1187IN=0.13371N弹簧被压縮时于给定虚拟接触力F-8N，前二层所有弹簧需要变形的时间为72=^+^=5^=5x10—、<r(这里r为力触觉在新刷新频率的倒数)，前二层弹簧的变形量与对应的弹性系数的乘积之和小于给定的虚拟接触力，且前两层弹簧变形总计需要的时延时间满足满足刷新频率1000Hz以上的要求，设共计经过时延r,，第二层的所有弹簧均被拉伸或压缩到极限长度，只有当第二层的所有弹簧均被拉伸或压縮到极限长度后，第三层对应的42个弹簧才同时开始被拉伸或压缩。第三层弹簧对应的弹性系数为&=0.15N/cmxe5=0.40774N/m,若在给定虚拟接触力的作用下，第三层的42个弹簧均被拉伸或压缩到极限长度的话，第三层每个弹簧对应的变形量为-线=A丽—^。Ig'—1=0.lcmx1.23—1=0.144c附弹簧被拉伸时|《min-^。k—^0.1cwxl.23-1=0.144cw弹簧被压縮时(第三层每个弹簧的变形量)x(第三层弹簧对应的弹性系数)x(第三层w敏人^f0.144cmx0.40774N/cnxl6=0.93943N弹簧被拉伸时共一p"弹費个数)=<|目!j二层所l0.144cwx0.40774N/cmxl6=0.93943N弹簧被压缩时一有弹簧的变形量与对应弹簧弹性系数的乘积之和]></math>id="icf0001"file="A2009100319350002C1.tif"wi="31"he="10"top="72"left="30"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>分布的白噪声信号，该白噪声信号为<mathsid="math0002"num="0002"><math><![CDATA[<mrow><mfencedopen='{'close=''><mtable><mtr><mtd><mi>σ</mi><mo>·</mo><mi>x</mi><mo>=</mo><mi>σ</mi><mo>·</mo><msqrt><mo>-</mo><mn>2</mn><mi>log</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>rand</mi><mo>(</mo><mo>)</mo><mo>/</mo><mn>32768</mn><mo>)</mo></mrow></msqrt><mo>·</mo><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>π</mi><mo>·</mo><mi>rand</mi><mo>(</mo><mo>)</mo><mo>/</mo><mn>32768</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>σ</mi><mo>=</mo><mi>K</mi><mo>·</mo><mfrac><mn>1</mn><mrow><mi>d</mi><mo>+</mo><msub><mi>d</mi><mn>0</mn></msub></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math></maths>其中σ为白噪声的标准差，log()为自然对数，rand()为0-32767间的一个随机整数，K为正比例系数，d为虚拟代理与虚拟软组织表面之间的距离，d0为距离阈值，d0＞0；步骤3当检测到虚拟代理碰撞到虚拟软组织表面上任何一点时，在给定虚拟接触力F作用下，虚拟代理与虚拟软组织交互的局部区域内部填充该等比例层状并联力触觉建模方法，在交互过程中，反馈输出为服从<mathsid="math0003"num="0003"><math><![CDATA[<mrow><mi>N</mi><mrow><mo>(</mo><mn>0</mn><mo>,</mo><msup><mrow><mo>(</mo><mi>K</mi><mo>·</mo><mfrac><mn>1</mn><msub><mi>d</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow></mrow>]]></math>id="icf0003"file="A2009100319350002C3.tif"wi="26"he="10"top="160"left="133"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>恒定的白噪声信号与采用等比例层状并联力触觉建模方法计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息叠加在一起的信号，所述等比例层状并联力触觉的建模方法为(1)参数初始化，(2)当虚拟代理碰撞到软组织表面上的任意点时，在碰撞点处下方悬挂一弹簧，形成第一层；在第一层40个弹簧下悬挂十字杆，在十字杆的各个端部下各悬挂一弹簧，形成由41个相同弹簧组成的第二层；在第二层的各个弹簧下，分别悬挂十字杆，每根十字杆的各个端部下各悬挂一弹簧，形成由42个相同弹簧组成的第三层；依此类推，在第i-1层的各个弹簧下，分别悬挂十字杆，每根十字杆的各个端部下各悬挂一弹簧，形成由4i-1个相同弹簧组成的第i层，i＝1，2，3，……，N，N为自然数，处于同一层的弹簧系数及变形量均相同，相邻层的弹簧弹性系数之比为常数，即第i层任一弹簧的弹性系数ki与第一层弹簧的弹性系数k1呈如下式所示的指数倍关系<mathsid="math0004"num="0004"><math><![CDATA[<mrow><mfrac><msub><mi>k</mi><mi>i</mi></msub><msub><mi>k</mi><mn>1</mn></msub></mfrac><mo>=</mo><msup><mi>e</mi><mfrac><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>β</mi></mfrac></msup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math></maths>在这里令β为常数；假定第i层所有弹簧的原始长度、拉伸及压缩后的极限长度、弹簧的变形量分别为Xio、Ximax、Ximin、ΔXi，它们均以第一层弹簧的对应量为首项，组成以q为公比的等比数列，所述的十字杆由长度相同的第一杆和第二杆组成，第一杆与第二杆垂直且其交点与第一杆和第二杆的质心相重合，(3)使虚拟接触力作用于软组织碰撞点，第i层对应的4i-1个弹簧同时开始被拉伸或压缩，如果前i层所有弹簧的变形量与对应的弹性系数的乘积之和小于给定的虚拟接触力，且第i层弹簧变形总计需要的时延时间满足满足刷新频率1000Hz以上的要求，设共计经过时延Ti，第i层的所有弹簧均被拉伸或压缩到极限长度，只有当第i层的所有弹簧均被拉伸或压缩到极限长度后，第i+1层对应的4i个弹簧才同时开始被拉伸或压缩，依此类推，直到前M层所有弹簧的变形量与对应的弹性系数的乘积之和不小于给定的虚拟接触力，或第M层弹簧变形总计需要的时延时间TM不满足刷新频率的要求，则第M层为变形截止层，此时前M-1层的所有弹簧均被拉伸或压缩到极限长度，第M层的所有弹簧被拉伸或压缩的长度不大于极限长度，第M层任一弹簧对应的变形量为其中ΔXM为第M层任一弹簧对应的变形量，F为虚拟接触力，ki为第i层任一弹簧的弹性系数，X1max为第一层弹簧被拉伸后的极限长度，X1o为第一层弹簧的原始长度，X1min为第一层弹簧被压缩后的极限长度，q为公比系数，PM为第M层弹簧的个数，Pi为第i层弹簧的个数为；用ti、Ti分别表示第i层的弹簧变形需要的时间、前i层的所有弹簧变形需要的时间，且令层间的时延时间满足ti＝4i-1t1(5)从虚拟代理碰撞接触到虚拟软组织表面算起，假定第i层的弹簧变形总计需要时延的时间Ti，且其须满足Ti＜T，其中<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><msub><mi>t</mi><mn>1</mn></msub><mo>+</mo><msub><mi>t</mi><mn>2</mn></msub><mo>+</mo><msub><mi>t</mi><mn>3</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow>]]></math></maths><mathsid="math0006"num="0006"><math><![CDATA[<mrow><mo>=</mo><mfrac><mrow><msup><mn>4</mn><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mn>1</mn></mrow><mn>3</mn></mfrac><msub><mi>t</mi><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math></maths>T为力触觉再现刷新频率的倒数；等比例层状并联力触觉建模方法中每层各个弹簧变形量之和的叠加对外等效为物体表面的变形，与之相连的弹簧弹性力的合力等效为物体表面接触力。全文摘要本发明公开了一种用于虚拟现实人机交互的基于物理意义的增强力触觉建模方法，其特征是虚拟代理未碰撞到虚拟软组织之前，仅反馈输出与软组织表面距离有关的白噪声信号。而交互过程中，反馈输出人手皮肤可敏感感知的白噪声信号与采用等比例层状并联力触觉建模方法计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息叠加在一起的信号；等比例层状并联力触觉建模方法中每层各个弹簧变形量之和的叠加对外等效为物体表面的变形，与之相连的弹簧弹性力的合力等效为物体表面接触力。该建模方法计算简单，能准确快速的计算变形，实现对软组织的实时变形仿真。在交互过程中，自然舒适、力触觉感觉平稳、模拟效果逼真。文档编号G06F3/01GK101587372SQ20091003193公开日2009年11月25日申请日期2009年7月3日优先权日2009年7月3日发明者伟孙,宋爱国,张小瑞,李佳璐,楠王,程盈盈申请人:东南大学