CN103678806A - 一种模拟柔性体扭转变形的盘香弹簧模型 - Google Patents
一种模拟柔性体扭转变形的盘香弹簧模型 Download PDFInfo
- Publication number
- CN103678806A CN103678806A CN201310680358.9A CN201310680358A CN103678806A CN 103678806 A CN103678806 A CN 103678806A CN 201310680358 A CN201310680358 A CN 201310680358A CN 103678806 A CN103678806 A CN 103678806A
- Authority
- CN
- China
- Prior art keywords
- coil spring
- incense coil
- circle
- torsion
- flexible body
- 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
- Developing Agents For Electrophotography (AREA)
- Springs (AREA)
Abstract
本发明提出了一种模拟柔性体扭转变形的盘香弹簧模型,所述盘香弹簧模型由多个盘香弹簧依次串接组成,在交互过程中,输出回馈为采用盘香弹簧模型计算出来的反应在扭矩作用下柔性体实时变形仿真的力触觉信息的信号,该盘香弹簧模型中所有圈上产生扭转变形量之和的叠加对外等效为柔性体表面的变形。本发明盘香弹簧模型每圈盘香弹簧的扭转变形量计算方法相同,计算简单,加快了扭转变形计算的速度;通过调节盘香弹簧可旋转芯轴半径,盘香弹簧的厚度和悬挂截面宽度等,就可模拟不同类型的柔性体,适用性广;可应用于虚拟外科手术仿真、遥控操作机器人控制、远程医疗等领域。
Description
技术领域
本发明属于计算机仿真技术领域,尤其涉及一种模拟柔性体扭转变形的盘香弹簧模型。
背景技术
虚拟手术作为一种新式的手术训练手段,即在计算机建立的虚拟环境中进行实时手术仿真已成为近几年研究的热点。虚拟手术仿真系统除在视觉上给使用者营造一种沉浸感、真实感之外,还要提供一种力触觉反馈机制,使人们在看见虚拟对象的同时,还可以触摸到它的几何形状和物理特征,视觉和触觉的有机结合使虚拟现实、虚拟仿真更加真实可信。
准确性和实时性是虚拟手术系统至关重要的性能指标,良好的建模方法将为此提供必要的保障。目前,主要有质点-弹簧模型,具有简单易行、计算量较小、易于实现等优点,但存在拓扑分析困难等弊端;有限元建模方法,虽然可对物体的形变进行精确和定量的模拟,但涉及大量繁杂的计算,实时交互性能差;边界元建模方法虽对建模方法的边界进行离散,简化了计算,但在稳定性方面却存在一定的不足。因此,深入开展用于柔性体扭转变形仿真技术建模的研究无疑对推动虚拟手术的发展具有重要意义。
发明内容
本发明所要解决的技术问题在于克服现有技术的不足,提供一种模拟柔性体扭转变形的盘香弹簧模型。
为了解决上述技术问题,本发明所采用的技术方案是:
一种模拟柔性体扭转变形的盘香弹簧模型,包括如下步骤:
步骤1,建立空间直角坐标系,确定柔性体表面上任意点处安放的盘香弹簧,其过程如下:
步骤1.1,建立空间直角坐标系,
在给定扭矩M作用下,在柔性体表面任意点O处设置一盘香弹簧;在距离任意点O为r处,选取盘香弹簧可旋转芯轴外侧点A0,r为盘香弹簧可旋转芯轴半径;以任意点O为原点,以射线OA0所在方向为X轴正方向,建立XYZ空间坐标系;
步骤1.2,按圈设置盘香弹簧,
其中,第i圈盘香弹簧是以原点O(0,0,0)为圆心,以Ri=r+ih为半径的圆,h为第i圈盘香弹簧的厚度,i=1,2,3,……S,S为自然数;
步骤2,确定任一圈盘香弹簧所消耗的扭矩;
设定给定扭矩M的作用线与盘香弹簧可旋转芯轴外侧点A0所在的圆相切,且在给定扭矩M作用下,如果柔性体中共有前P圈盘香弹簧产生扭转变形,则第P圈盘香弹簧称为扭转变形截止圈;P≤S,也就是盘香弹簧的圈数至少等于P;
根据盘香弹簧特性,设定:
前P-1圈盘香弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力相等且均为F0;第P圈盘香弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力FP′均相等,且不大于F0;
步骤2.1,确定前P-1圈盘香弹簧中,第i圈上所有点共消耗的扭矩Mi为:
Mi=li·F0=2πRi·F0
其中,li为第i圈盘香弹簧的有效长度,Ri为第i圈盘香弹簧的半径,
步骤2.2,确定第P圈盘香弹簧所有点共消耗的扭矩MP′为:
第P圈盘香弹簧在给定扭矩M作用下,任意一点消耗的最大扭转力FP′为:
步骤3,计算第i圈盘香弹簧上所有点,在共消耗的扭矩Mi作用下,产生的扭转圈数ni;
其中,li为第i圈盘香弹簧的有效长度,其表达式如下:
所述模拟柔性体扭转变形的盘香弹簧模型,所述前P圈盘香产生扭转变形总计需要的时延时间不大于1ms,也即满足刷新频率不低于1000Hz的要求。
所述模拟柔性体扭转变形的盘香弹簧模型,其所有盘香弹簧圈数中,每圈产生扭转变形需要的时延时间构成等比数列,即满足:
ti=qi-1t1
其中,ti表示第i圈盘香弹簧产生扭转变形需要的时延时间,t1为第1圈盘香弹簧产生扭转变形需要的时延时间,1≤i≤P。
所述模拟柔性体扭转变形的盘香弹簧模型,任一圈盘香弹簧上任意一点消耗的扭转力达到最大扭转力后,其下一圈开始产生扭转变形。
本发明的有益效果是:本发明提出了一种模拟柔性体扭转变形的盘香弹簧模型,所述盘香弹簧模型由多个盘香弹簧依次串接组成,在交互过程中,输出回馈为采用盘香弹簧模型计算出来的反应在扭矩作用下柔性体实时变形仿真的力触觉信息的信号,该盘香弹簧模型中所有圈上产生扭转变形量之和的叠加对外等效为柔性体表面的变形。本发明盘香弹簧模型每圈盘香弹簧的扭转变形量计算方法相同,计算简单,加快了扭转变形计算的速度;通过调节盘香弹簧可旋转芯轴半径,盘香弹簧的厚度和悬挂截面宽度等,就可模拟不同类型的柔性体,适用性广;可应用于虚拟外科手术仿真、遥控操作机器人控制、远程医疗等领域。
附图说明
图1是盘香弹簧模型示意图。
图2是柔性体变形仿真流程图。
图3是盘香弹簧模型构建方法流程图。
图4是盘香弹簧模型扭矩、扭转圈数与时延时间关系示意图。
具体实施方式
下面结合附图所示流程对本发明提出的一种模拟柔性体扭转变形的盘香弹簧模型进行详细说明:
如图1所示的盘香弹簧模型示意图。一种模拟柔性体扭转变形的盘香弹簧模型,其具体步骤如下:
步骤1,建立空间直角坐标系,确定柔性体表面上任意点处安放的盘香弹簧,其过程如下:
步骤1.1,建立空间直角坐标系,
在给定扭矩M作用下,在柔性体表面任意点O处设置一盘香弹簧;在距离任意点O为r处,选取盘香弹簧可旋转芯轴外侧点A0,r为盘香弹簧可旋转芯轴半径;以任意点O为原点,以射线OA0所在方向为X轴正方向,建立XYZ空间坐标系;
步骤1.2,按圈设置盘香弹簧,
其中,第i圈盘香弹簧是以原点O(0,0,0)为圆心,以Ri=r+ih为半径的圆,h为第i圈盘香弹簧的厚度,i=1,2,3,……S,S为自然数;
步骤2,确定任一圈盘香弹簧所消耗的扭矩;
设定给定扭矩M的作用线与盘香弹簧可旋转芯轴外侧点A0所在的圆相切,且在给定扭矩M作用下,如果柔性体中共有前P圈盘香弹簧产生扭转变形,则第P圈盘香弹簧称为扭转变形截止圈;P≤S,也就是盘香弹簧的圈数至少等于P;
根据盘香弹簧特性,设定:
前P-1圈盘香弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力相等且均为F0;第P圈盘香弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力FP′均相等,且不大于F0;
步骤2.1,确定前P-1圈盘香弹簧中,第i圈上所有点共消耗的扭矩Mi为:
Mi=li·F0=2πRi·F0
其中,li为第i圈盘香弹簧的有效长度,Ri为第i圈盘香弹簧的半径,
步骤2.2,确定第P圈盘香弹簧所有点共消耗的扭矩MP′为:
第P圈盘香弹簧在给定扭矩M作用下,任意一点消耗的最大扭转力FP′为:
步骤3,计算第i圈盘香弹簧上所有点,在共消耗的扭矩Mi作用下,产生的扭转圈数ni;
其中,li为第i圈盘香弹簧的有效长度,其表达式如下:
所述模拟柔性体扭转变形的盘香弹簧模型,所述前i圈盘香产生扭转变形总计需要的时延时间不大于1ms,也即满足刷新频率不低于1000Hz的要求。
所述模拟柔性体扭转变形的盘香弹簧模型,其所有盘香弹簧圈数中,每圈产生扭转变形需要的时延时间构成等比数列,即满足:
ti=qi-1t1
其中,ti表示第i圈盘香弹簧产生扭转变形需要的时延时间,t1为第1圈盘香弹簧产生扭转变形需要的时延时间,1≤i≤P。
所述模拟柔性体扭转变形的盘香弹簧模型,任一圈盘香弹簧上任意一点消耗的扭转力达到最大扭转力后,其下一圈开始产生扭转变形。
下面以虚拟手和虚拟血管模型为例,列举本发明技术方案的具体实施方式。
本实例中所有虚拟手和虚拟血管模型都直接采用从3DS MAX2013软件中导出的OBJ格式,以599个质点,1181个三角网格构成的虚拟手和2864个质点,6122个三角网格构成的虚拟血管模型为例来进行变形仿真,实验过程中模型获取和修改非常方便;操作系统为Windows2000,以3DS MAX2013、OpenGL图形库为基础,在VC++2012软件开发平台上进行仿真。
当检测到虚拟手碰撞到虚拟血管表面上任何一点时,在给定逆时针扭矩M=26×10-1N·m作用下,虚拟手与虚拟血管交互的局部区域内部填充盘香弹簧模型,在交互过程中,输出反馈为采用盘香弹簧模型计算出来的反应在给定扭矩作用下虚拟血管实时变形仿真的力触觉信息的信号,如图2所示;
在给定逆时针扭矩M=26×10-1N·m作用下,当虚拟手碰撞到虚拟血管表面上任意点O时,在任意点O处安放一盘香弹簧,在距离任意点O点为r=2×10-2m处设置一盘香弹簧可旋转芯轴外侧点A0,r=2×10-2m为盘香弹簧可旋转芯轴半径,以任意点O为原点,以射线OA0所在方向为X轴正方向,建立XYZ空间坐标系;
按圈设置盘香弹簧,如图1和图3所示,其中,第i圈盘香弹簧是以原点O(0,0,0)为圆心,以Ri=r+ih为半径的圆,h=1×10-2m为任意第i圈盘香弹簧的厚度,i=1,2,3,……S,S为自然数;
即以O(0,0,0)为圆心,
以R1=r+h=2×10-2+1×10-2=3×10-2m为半径画圆,形成第1圈盘香弹簧;
以R2=r+2h=2×10-2+2×1×10-2=4×10-2m为半径画圆,形成第2圈盘香弹簧;
以R3=r+3h=2×10-2+3×1×10-2=5×10-2m为半径画圆,形成第3圈盘香弹簧;
以R4=r+4h=2×10-2+4×1×10-2=6×10-2m为半径画圆,形成第4圈盘香弹簧;
以R5=r+5h=2×10-2+5×1×10-2=7×10-2m为半径画圆,形成第5圈盘香弹簧;
第i圈盘香弹簧的半径Ri形成以r=2×10-2m为首项,以h=1×10-2m为公差的等差数列;
计算的中间过程、最后的数据均按四舍五入法保留小数点后3位。
第i圈盘香弹簧的悬挂截面宽度均为b=2×10-3m、截面惯性矩均为 弹性模量均为E=3.09×107Pa,取决于柔性体的材质,且柔性体材质均相同;
假定给定逆时针扭矩M的作用线和盘香弹簧可旋转芯轴外侧点A0所在的圆相切,且在给定逆时针扭矩M作用下,如果柔性体中共有前P圈盘香弹簧产生扭转变形,则第P圈盘香弹簧称为扭转变形截止圈;P≤S,也就是盘香弹簧的圈数至少等于P;
根据盘香弹簧特性设定:如图4所示,前P-1圈盘香弹簧上任意一点在给定逆时针扭矩M作用下,所消耗的最大扭转力相等并均为F0=2.06N,第P圈盘香弹簧上任意一点在给定逆时针扭矩M作用下,所消耗的最大扭转力FP′均相等,且不大于F0=2.06N;
前i圈盘香弹簧上产生扭转变形总计需要的时延时间,满足刷新频率不低于1000Hz的要求;且所有盘香弹簧圈数中,每圈产生扭转变形需要的时延时间构成以第1圈盘香弹簧产生扭转变形需要的时延时间t1=10-5s为首项、以q=1.2为公比的等比数列;在盘香弹簧任意一圈上的任意一点消耗的扭转力达到最大扭转力后,其下一圈开始产生扭转变形。
假定力触觉再现刷新频率为1200Hz,则力触觉再现刷新频率的倒数
若第1圈盘香弹簧上任一点在给定逆时针扭矩M作用下,所消耗的最大扭转力相等并均到达F0,则第1圈盘香弹簧上所有点共消耗的扭矩M1为:
第1圈盘香弹簧的有效长度l1=2πR1,
M1=l1·F0=2πR1·Fo=2×3.142×3×10-2×2.06=3.884×10-1N·m
M1<M;
第1圈盘香弹上所有点,在共消耗的扭矩M1作用下,与其产生的扭转圈数n1之间满足:
第1圈盘香弹簧上产生扭转变形需要的时延时间T1=t1=10-5s<T;
因此,第1圈香弹弹簧上所有点共消耗的扭矩M1<M,且第1圈盘香弹簧产生扭转变形总计需要的时延时间T1=10-5s<T,满足刷新频率的要求;只有当第1圈香弹弹簧上任意一点消耗的扭转力到消耗的最大扭转力F0后,第2圈盘香卷弹簧才开始产生扭转变形。
若第2圈盘香弹簧上任一点在给定逆时针扭矩M作用下,所消耗的最大扭转力相等并均到达F0,则第2圈盘香弹簧上所有点共消耗的扭矩M2为:
第2圈盘香弹簧的有效长度l2=2πR2,
M2=l2·F0=2πR2·Fo=2×3.142×4×10-2×2.06=5.178×10-1N·m
前2圈盘香弹簧上所有点共消耗的扭矩之和为:
M1+M2=(3.884+5.178)×10-1=9.062×10-1N·m<M
第2圈盘香弹上所有点,在共消耗的扭矩M2作用下,与其产生的扭转圈数n2之间满足:
前2圈盘香弹簧产生扭转变形需要的时延时间T2=t1+t2=(1+q)t1=(1+1.2)×10-5=2.2×10-5s<T,这里T为力触觉再现刷新频率的倒数,
因此,前2圈盘香弹簧上所有点共消耗的扭矩M1+M2<M,且前2圈盘香弹簧产生扭转变形总计需要的时延时间T2=2.2×10-5s<T,满足刷新频率的要求;只有当第2圈盘香弹簧上任意一点消耗的扭转力达到消耗的最大扭转力F0后,第3圈盘香弹簧才开始产生扭转变形。
若第3圈盘香弹簧上任一点在给定逆时针扭矩M作用下,所消耗的最大扭转力相等并均到达F0,则第3圈盘香弹簧上所有点共消耗的扭矩M3为:
第3圈盘香弹簧的有效长度l3=2πR3,
M3=2πR3·Fo=2×3.142×5×10-2×2.06=6.473×10-1N·m
前3圈盘香弹簧上所有点共消耗的扭矩之和为:
M1+M2+M3=(3.884+5.178+6.473)×10-1=15.535×10-1N·m<M
第3圈盘香弹上所有点,在共消耗的扭矩M3作用下,与其产生的扭转圈数n3之间满足:
前3圈盘香弹簧产生扭转变形需要的时延时间T3=t1+t2+t3=(1+q+q2)t1=(1+1.2+1.22)×10-5=3.64×10-5s<T,这里T为力触觉再现刷新频率的倒数,
因此,前3圈盘香弹簧上所有点共消耗的扭矩M1+M2+M3<M,且前3圈盘香弹簧产生扭转变形总计需要的时延时间T3=3.64×10-5s<T,满足刷新频率的要求;只有当第3圈盘香弹簧上任意一点消耗的扭转力达到消耗的最大扭转力F0后,第4圈盘香弹簧才开始产生扭转变形。
若第4圈盘香弹簧上任一点在给定逆时针扭矩M作用下,所消耗的最大扭转力相等并均到达F0,则第4圈盘香弹簧上所有点共消耗的扭矩M4为:
第4圈盘香弹簧的有效长度l4=2πR4,
M4=2πR4·Fo=2×3.142×6×10-2×2.06=7.767×10-1N·m
前4圈盘香弹簧上所有点共消耗的扭矩之和为:
M1+M2+M3+M4=(3.884+5.178+6.473+7.767)×10-1=23.302×10-1N·m<M
第4圈盘香弹上所有点,在共消耗的扭矩M4作用下,与其产生的扭转圈数n4之间满足:
前4圈盘香弹簧产生扭转变形需要的时延时间T4=t1+t2+t3+t4=(1+q+q2+q3)t1=(1+1.2+1.22+1.23)×10-5=5.368×10-5s<T,这里T为力触觉再现刷新频率的倒数,
因此,前4圈盘香弹簧上所有点共消耗的扭矩M1+M2+M3+M4<M,且前4圈盘香弹簧产生扭转变形总计需要的时延时间T4=5.368×10-5s<T,满足刷新频率的要求;只有当第4圈盘香弹簧上任意一点消耗的扭转力达到消耗的最大扭转力F0后,第5圈盘香弹簧才开始产生扭转变形。
若第5圈盘香弹簧上任一点在给定逆时针扭矩M作用下,所消耗的最大扭转力相等并均到达F0,则第5圈盘香弹簧上所有点共消耗的扭矩M5为:
第5圈盘香弹簧的有效长度l5=2πR5,
M5=2πR5·Fo=2×3.142×7×10-2×2.06=9.062×10-1N·m
前5圈盘香弹簧上所有点共消耗的扭矩之和为:
M1+M2+M3+M4+M5=(3.884+5.178+6.473+7.767+9.062)×10-1=32.364×10-1N·m>M=26×10-1N·m
因此,前5圈盘香弹簧上所有点共消耗的扭矩M1+M2+M3+M4+M5>M,前5圈盘香弹簧上所有点共消耗的扭矩之和不小于给定逆时针扭矩M,则第5圈盘香弹簧为变形截止圈,不需要再判断是否满足刷新频率的要求;
第5圈盘香弹簧即变形截止圈上所有点共消耗的扭矩M5′为:
变形截止圈即第5圈盘香弹簧上任意一点在给定逆时针扭矩M作用下,消耗的最大扭矩F5′为:
根据盘香卷弹簧特性设定:前4圈盘香弹簧上任意一点在给定逆时针扭矩M作用下,消耗的最大扭转力相等并都假定为F0=2.06N,变形截止圈也即第5圈盘香弹簧在给定逆时针扭矩M作用下,任意一点消耗的最大扭转力F5′=0.613N均相等并都不大于F0=2.06N;
第5圈盘香簧即变形截止圈上所有点,在共消耗的扭矩M5′作用下,与其产生的扭转圈数n5之间满足:
前5圈盘香弹簧在给定逆时针扭矩M=26×10-1N·m作用下,共产生的扭转转圈数n为:
n=n1+n2+n3+n4+n5=2.828+5.026+7.854+11.309+1.365=28.382
注意:在采用柔性体变形仿真的盘香弹簧模型,来计算在给定扭矩作用下柔性体实时变形仿真的过程中,若r、h、b这些参数选取的过大,则柔性体变形仿真的盘香弹簧模型的变形圈数就少,计算量小,实时性好,但变形仿真效果不佳;如果r、h、b这些参数选取的过小,则柔性体变形仿真的盘香弹簧模型的变形圈数就越多,计算量大,实时性不佳,但变形仿真效果较好;另外在设置t1和ti之间的比例关系时,要考虑到程序运行时计算机本身的硬件配置,故在调试整个程序的过程中,要折中选择这些参数,不断反复调试,从而使变形效果更加逼真。
为验证本发明的实施效果,操作者通过PHANTOM OMNI手控器端部的手柄来触摸、感知和控制虚拟手对虚拟血管进行扭转的变形仿真,并将交互过程中产生的力触觉信息实时反馈给操作者。实验结果表明,该模型算法简单、计算速度快、力触觉反馈平稳,形变效果逼真,能够满足虚拟手术仿真交互系统实时性和稳定性的要求,从而使操作者对虚拟环境的感知和交互更加准确真实。
Claims (4)
1.一种模拟柔性体扭转变形的盘香弹簧模型,其特征在于,包括如下步骤:
步骤1,建立空间直角坐标系,确定柔性体表面上任意点处安放的盘香弹簧,其过程如下:
步骤1.1,建立空间直角坐标系,
在给定扭矩M作用下,在柔性体表面任意点O处设置一盘香弹簧;在距离任意点O为r处,选取盘香弹簧可旋转芯轴外侧点A0,r为盘香弹簧可旋转芯轴半径;以任意点O为原点,以射线OA0所在方向为X轴正方向,建立XYZ空间坐标系;
步骤1.2,按圈设置盘香弹簧,
其中,第i圈盘香弹簧是以原点O(0,0,0)为圆心,以Ri=r+ih为半径的圆,h为第i圈盘香弹簧的厚度,i=1,2,3,……S,S为自然数;
步骤2,确定任一圈盘香弹簧所消耗的扭矩;
设定给定扭矩M的作用线与盘香弹簧可旋转芯轴外侧点A0所在的圆相切,且在给定扭矩M作用下,如果柔性体中共有前P圈盘香弹簧产生扭转变形,则第P圈盘香弹簧称为扭转变形截止圈;P≤S,也就是盘香弹簧的圈数至少等于P;
根据盘香弹簧特性,设定:
前P-1圈盘香弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力相等且均为F0;第P圈盘香弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力FP′均相等,且不大于F0;
步骤2.1,确定前P-1圈盘香弹簧中,第i圈上所有点共消耗的扭矩Mi为:
Mi=li·F0=2πRi·F0
其中,li为第i圈盘香弹簧的有效长度,Ri为第i圈盘香弹簧的半径,
步骤2.2,确定第P圈盘香弹簧所有点共消耗的扭矩MP′为:
第P圈盘香弹簧在给定扭矩M作用下,任意一点消耗的最大扭转力FP′为:
步骤3,计算第i圈盘香弹簧上所有点,在共消耗的扭矩Mi作用下,产生的扭转圈数ni;
其中,li为第i圈盘香弹簧的有效长度,其表达式如下:
2.根据权利要求1所述的一种模拟柔性体扭转变形的盘香弹簧模型,其特征在于,所述前P圈盘香产生扭转变形总计需要的时延时间不大于1ms。
3.根据权利要求1所述的一种模拟柔性体扭转变形的盘香弹簧模型,其特征在于,所有盘香弹簧圈数中,每圈产生扭转变形需要的时延时间构成等比数列,即满足:
ti=qi-1t1
其中,ti表示第i圈盘香弹簧产生扭转变形需要的时延时间,t1为第1圈盘香弹簧产生扭转变形需要的时延时间,1≤i≤P。
4.根据权利要求1或2或3所述的一种模拟柔性体扭转变形的盘香弹簧模型,其特征在于,所述盘香弹簧模型,在任一圈盘香弹簧上任意一点消耗的扭转力达到最大扭转力后,其下一圈开始产生扭转变形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310680358.9A CN103678806B (zh) | 2013-12-12 | 2013-12-12 | 一种模拟柔性体扭转变形的盘香弹簧模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310680358.9A CN103678806B (zh) | 2013-12-12 | 2013-12-12 | 一种模拟柔性体扭转变形的盘香弹簧模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103678806A true CN103678806A (zh) | 2014-03-26 |
CN103678806B CN103678806B (zh) | 2016-07-06 |
Family
ID=50316342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310680358.9A Expired - Fee Related CN103678806B (zh) | 2013-12-12 | 2013-12-12 | 一种模拟柔性体扭转变形的盘香弹簧模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103678806B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104615879A (zh) * | 2015-01-28 | 2015-05-13 | 南京信息工程大学 | 一种增强力触觉再现的圆柱扭转弹簧虚拟模型 |
-
2013
- 2013-12-12 CN CN201310680358.9A patent/CN103678806B/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
孙伟等: "Real-Time haptic model for soft tissue deformation in surgery simulation", 《高技术通讯:英文版》 * |
崔桐等: "一种用于力觉再现的柔性体变形仿真弹簧-质点模型", 《东南大学学报(自然科学版)》 * |
张小瑞等: "用于实时人机交互的增强力触觉模型", 《高技术通讯》 * |
张小瑞等: "用于软组织变形仿真的层状菱形连接模型", 《计算机辅助设计与图形学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104615879A (zh) * | 2015-01-28 | 2015-05-13 | 南京信息工程大学 | 一种增强力触觉再现的圆柱扭转弹簧虚拟模型 |
Also Published As
Publication number | Publication date |
---|---|
CN103678806B (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101587372B (zh) | 一种用于虚拟现实人机交互的增强力触觉建模方法 | |
CN101976298B (zh) | 增强力触觉再现的对称式板弹簧虚拟模型的建模方法 | |
CN101964023B (zh) | 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 | |
CN101853072B (zh) | 一种用于柔性体变形仿真的力触觉建模方法 | |
CN103400023B (zh) | 软组织形变仿真方法 | |
CN101976299B (zh) | 柔性力触觉再现的对称式板弹簧虚拟模型的建模方法 | |
CN103869983B (zh) | 一种用于力触觉人机交互的柔性物体变形仿真方法 | |
CN105184030A (zh) | 基于弯扭复合弹簧质点模型柔性线缆位姿模拟方法及装置 | |
CN101496028A (zh) | 使用几何推动式模型模拟可变形物体的方法 | |
CN107992672A (zh) | 一种基于虚拟弹簧的软组织形变建模方法 | |
CN106709203A (zh) | 一种基于Texgen纬编织物模型的三维模拟方法 | |
CN104063057B (zh) | 一种支持增强力触觉反馈的涡卷弹簧仿真系统 | |
CN103745058B (zh) | 一种任意形状软组织表皮上受拉力/变形的模拟方法 | |
CN103793552A (zh) | 一种软组织形变的局部质点弹簧模型的实时动态生成方法 | |
CN103729879A (zh) | 基于力觉计算的虚拟手稳定抓取方法 | |
CN107146288A (zh) | 虚拟手术中支持实时按压形变的软组织模型建模方法 | |
CN103678806A (zh) | 一种模拟柔性体扭转变形的盘香弹簧模型 | |
CN103699753B (zh) | 一种支持增强力触觉反馈的多股螺旋弹簧虚拟模型 | |
CN108428256A (zh) | 一种基于柔软度的自适应网格细化的软组织形变仿真方法 | |
CN101901565B (zh) | 支持力触觉反馈的虚拟柔性体变形手术仿真系统 | |
CN107080586A (zh) | 虚拟手术中用于软组织实时形变的交叉弹簧模型 | |
CN103714205B (zh) | 一种模拟柔性体可旋转变形的发条弹簧模型 | |
CN106485028B (zh) | 支持实时按压的柔性体变形仿真建模方法 | |
CN103699741A (zh) | 一种增强模拟柔性体可旋转变形的发条弹簧模型 | |
CN104063540A (zh) | 一种支持增强力触觉反馈的涡卷弹簧模型 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 Termination date: 20181212 |