发明名称 一种增强力触觉再现的圆柱扭转弹簧虚拟模型
摘要 本发明提出了一种增强力触觉再现的圆柱扭转弹簧虚拟模型,所述圆柱扭转弹簧模型由多个圆柱扭转弹簧依次串接组成,在交互过程中,输出回馈为采用圆柱扭转虚拟模型计算出来的反应在给定虚拟扭矩作用下柔性体实时变形仿真的力触觉信息的信号,该圆柱扭转弹簧虚拟模型中所有层上扭转变形量之和的叠加对外等效为柔性体表面的变形。本发明圆柱扭转弹簧虚拟模型每层圆柱扭转弹簧扭转变形量计算方法相同,计算简单,加快了扭转变形计算的速度;通过调节圆柱扭转弹簧第一层半径,和圆柱扭转弹簧弹簧丝的直径,就可模拟不同类型的柔性体,适用性广;可应用于虚拟外科手术仿真、遥控操作机器人控制、远程医疗等领域。
申请公布号 CN104615879A 申请公布日期 2015.05.13
申请号 CN201510044754.1 申请日期 2015.01.28
申请人 南京信息工程大学 发明人 张小瑞;黄昆;孙伟;徐骁晗;朱建栋;宋爱国;牛建伟;庄颖
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 一种增强力触觉再现的圆柱扭转弹簧虚拟模型,其特征在于,包括如下步骤:步骤1,对虚拟场景进行初始化;步骤2,在检测到虚拟代理碰撞到虚拟柔性体之前,即在靠近虚拟柔性体的过程中,反馈输出散粒噪声,散粒噪声的噪声功率谱S为:S=2eI<sub>g</sub>式中,e为电子电量,I<sub>g</sub>为电流;其中,电流I<sub>g</sub>的表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>I</mi><mi>g</mi></msub><mo>=</mo><mfrac><msup><mrow><mn>2</mn><mi>e</mi></mrow><mn>2</mn></msup><mi>h</mi></mfrac><mi>V</mi><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>n</mi></msub></mrow>]]></math><img file="FDA0000663162970000011.GIF" wi="359" he="161" /></maths>式中,h为普朗克常量,V为施加电压,N为通道个数,T<sub>n</sub>为本征通道的投射系数;步骤3,当检测到虚拟代理碰撞到虚拟柔性体表面上任何一点时,在给定虚拟扭矩M作用下,虚拟代理与虚拟柔性体交互的局部区域内部填充圆柱扭转弹簧虚拟模型,在交互过程中,虚拟柔性体实时变形产生的力触觉信号与散粒噪声相叠加,作为反馈输出信号;所述圆柱扭转弹簧虚拟模型的建模方法为:步骤3‑1,建立空间直角坐标系,确定柔性体表面上任意点处安放的圆柱扭转弹簧,其过程如下:步骤3‑1‑1,建立空间直角坐标系;在给定虚拟扭矩M作用下,在柔性体表面任意点O处安放一圆柱扭转弹簧,以任意点O为原点,水平方向为X轴正方向,建立XYZ空间坐标系,即O的坐标为(0,0,0),在Z轴上距离原点O为d/2处记为O<sub>1</sub>点,即O<sub>1</sub>的坐标为(0,0,d/2),距离O<sub>1</sub>点为R处,设置圆柱扭转弹簧第一层弹簧丝的外侧点P<sub>0</sub>,即R为圆柱扭转弹簧的第一层半径;步骤3‑1‑2,按层设置圆柱扭转弹簧;在XYZ空间坐标系的Z轴方向,在圆柱扭转弹簧第一层弹簧丝下,依次设置圆柱扭转弹簧其它各层,且该圆柱扭转弹簧各层弹簧丝的直径均为d,圆柱扭转弹簧各层的半径均为R;该圆柱扭转弹簧的截面惯性矩为I,<img file="FDA0000663162970000012.GIF" wi="211" he="136" />弹性模量E取决于柔性体材质,且柔性体材质相同;步骤3‑2,确定圆柱扭转弹簧任一层消耗的扭矩;假定该圆柱扭转弹簧下端固定,在上端施加给定虚拟扭矩M,给定虚拟扭矩M的作用线与圆柱扭转弹簧第一层弹簧丝的外侧点P<sub>0</sub>处所在的圆相切,且在给定虚拟扭矩M作用下,如果柔性体中共有圆柱扭转弹簧的前C层产生扭转变形,则该圆柱扭转弹簧的第C层称为变形截止层;圆柱扭转弹簧的层数至少等于C;根据圆柱扭转弹簧特性,设定:圆柱扭转弹簧前C‑1层上任意一点在给定虚拟扭矩M作用下,所受到的作用力相等且均为F,圆柱扭转弹簧第C层上任意一点在给定虚拟扭矩M作用下,所消耗的最大扭转力F′<sub>C</sub>均相等,且不大于F;步骤3‑3,确定该圆柱扭转弹簧每层的螺旋角α,其表达式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&alpha;</mi><mo>=</mo><mi>arctan</mi><mfrac><mi>d</mi><mrow><mn>2</mn><mi>&pi;R</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000663162970000021.GIF" wi="324" he="134" /></maths>步骤3‑4,当i&lt;C时,确定圆柱扭转弹簧第i层消耗的扭矩M<sub>i</sub>,其表达式为:M<sub>i</sub>=iFR步骤3‑5,当i&lt;C时,确定圆柱扭转弹簧第i层所产生的扭转变形角<img file="FDA0000663162970000024.GIF" wi="77" he="66" />其表达式为:<img file="FDA0000663162970000022.GIF" wi="337" he="139" />其中,l<sub>i</sub>为该圆柱扭转弹簧第i层的工作长度,其表达式为:l<sub>i</sub>=2πiR步骤3‑6,确定变形截止层任意一点消耗的最大扭转力F′<sub>C</sub>,其表达式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>F</mi><mi>C</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><msubsup><mi>M</mi><mi>C</mi><mo>&prime;</mo></msubsup><mi>iR</mi></mfrac><mo>=</mo><mfrac><mrow><mi>M</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>C</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>M</mi><mi>i</mi></msub></mrow><mi>iR</mi></mfrac></mrow>]]></math><img file="FDA0000663162970000023.GIF" wi="482" he="206" /></maths>其中,M′<sub>C</sub>为变形截止层任意一点消耗的最大扭矩;步骤3‑7,确定变形截止层所产生的扭转变形角为<img file="FDA0000663162970000025.GIF" wi="90" he="60" />其表达式为:<img file="FDA0000663162970000031.GIF" wi="926" he="224" />其中,l<sub>C</sub>为该圆柱扭转弹簧变形截止层的工作长度,其表达式为:l<sub>C</sub>=2πCR步骤3‑8,将圆柱扭转弹簧每层的扭转角叠加,等效为虚拟柔性体表面的变形<img file="FDA0000663162970000033.GIF" wi="68" he="63" />其表达式如下:<img file="FDA0000663162970000032.GIF" wi="1059" he="291" />
地址 210044 江苏省南京市浦口区宁六路219号