发明名称 一种模拟柔性体可旋转变形的发条弹簧模型
摘要 本发明提出了一种模拟柔性体可旋转变形的发条弹簧模型,所述发条弹簧模型由多个发条弹簧依次串接组成,在交互过程中,输出回馈为采用发条弹簧模型计算出来的反应在扭矩作用下柔性体实时变形仿真的力触觉信息的信号,该发条弹簧模型中所有圈上产生扭转变形量之和的叠加对外等效为柔性体表面的变形。本发明发条弹簧模型每圈发条弹簧的扭转变形量计算方法相同,计算简单,加快了扭转变形计算的速度;通过调节发条弹簧可旋转芯轴半径,发条弹簧的厚度和悬挂截面宽度,就可模拟不同类型的柔性体,适用性广;可应用于虚拟外科手术仿真、遥控操作机器人控制、远程医疗等领域。
申请公布号 CN103714205A 申请公布日期 2014.04.09
申请号 CN201310713728.4 申请日期 2013.12.23
申请人 南京信息工程大学 发明人 张小瑞;朱建栋;孙伟;宋爱国;牛建伟;杨松;朱哲夫
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 1.一种模拟柔性体可旋转变形的发条弹簧模型,其特征在于,包括如下步骤:步骤1,建立空间直角坐标系,确定柔性体表面上任意点处安放的发条弹簧,其过程如下:步骤1.1,建立空间直角坐标系,在给定扭矩M作用下,在柔性体表面任意点O处安放一发条弹簧,在距离任意点O为r处设置一发条弹簧可旋转芯轴外侧点A<sub>0</sub>,即r为发条弹簧可旋转芯轴半径,以任意点O为原点,射线OA<sub>0</sub>所在方向为X轴正方向,建立XYZ空间坐标系;步骤1.2,按圈设置发条弹簧,其中,在第i-1圈发条弹簧外侧A<sub>i-1</sub>点,悬挂厚度为h<sub>i</sub>=h<sub>1</sub>+(i-1)d的第i圈发条弹簧,其中h<sub>1</sub>为第1圈发条弹簧的厚度,d为给定的常数且d&gt;0,即第i圈发条弹簧是以原点O(0,0,0)为圆心,以<img file="FDA0000444032470000011.GIF" wi="264" he="139" />为半径的圆,i=1,2,3,……S,S为自然数,每圈发条弹簧的宽均为b、截面惯性矩<img file="FDA0000444032470000012.GIF" wi="194" he="137" />弹性模量E取决于柔性体的材质,且柔性体材质均相同;步骤2,确定任一圈发条弹簧所消耗的扭矩;设定给定扭矩M的作用线与发条弹簧可旋转芯轴外侧点A<sub>0</sub>处所在的圆相切,且在给定扭矩M作用下,如果柔性体中共有前P圈发条弹簧产生扭转变形,则第P圈发条弹簧称为变形截止圈;且P≤S,也就是发条弹簧的圈数至少等于P;根据发条弹簧特性,设定:前P-1圈发条弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力相等且均为F<sub>0</sub>;第P圈发条弹簧上任意一点在给定扭矩M作用下,所消耗的最大扭转力<img file="FDA0000444032470000013.GIF" wi="64" he="75" />均相等,且不大于F<sub>0</sub>;步骤2.1,确定前P-1圈发条弹簧中,第i圈上所有点共消耗的扭矩M<sub>i</sub>为:<![CDATA[<math><mrow><msub><mi>M</mi><mi>i</mi></msub><mo>=</mo><msub><mi>F</mi><mn>0</mn></msub><mo>&CenterDot;</mo><msub><mi>R</mi><mi>i</mi></msub><mo>=</mo><msub><mi>F</mi><mn>0</mn></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>r</mi><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>h</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>步骤2.2,确定第P圈发条弹簧所有点共消耗的扭矩<img file="FDA0000444032470000026.GIF" wi="78" he="63" />为:<![CDATA[<math><mrow><msubsup><mi>M</mi><mi>P</mi><mo>&prime;</mo></msubsup><mo>=</mo><mi>M</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>P</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>M</mi><mi>i</mi></msub></mrow></math>]]></maths>第P圈发条弹簧上任意一点在给定扭矩M作用下,消耗的最大扭转力<img file="FDA0000444032470000027.GIF" wi="57" he="65" />为:<![CDATA[<math><mrow><msubsup><mi>F</mi><mi>P</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><msubsup><mi>M</mi><mi>P</mi><mo>&prime;</mo></msubsup><msub><mi>R</mi><mi>P</mi></msub></mfrac><mo>=</mo><mfrac><msubsup><mi>M</mi><mi>P</mi><mo>&prime;</mo></msubsup><mrow><mi>r</mi><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>P</mi></munderover><msub><mi>h</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths>步骤3,计算第i圈发条弹簧上所有点,在共消耗的扭矩M<sub>i</sub>作用下,产生的扭转圈数n<sub>i</sub>;<![CDATA[<math><mrow><msub><mi>n</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><msub><mi>M</mi><mi>i</mi></msub><msub><mi>l</mi><mi>i</mi></msub></mrow><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>EI</mi><mi>i</mi></msub></mrow></mfrac></mtd><mtd><mi>i</mi><mo>=</mo><mn>1,2,3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>P</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mfrac><mrow><msubsup><mi>M</mi><mi>P</mi><mo>&prime;</mo></msubsup><msub><mi>l</mi><mi>i</mi></msub></mrow><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>EI</mi><mi>i</mi></msub></mrow></mfrac></mtd><mtd><mi>i</mi><mo>=</mo><mi>P</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,l<sub>i</sub>为第i圈发条弹簧的有效长度,其表达式如下:<![CDATA[<math><mrow><msub><mi>l</mi><mi>i</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>R</mi><mn>1</mn></msub><mo>+</mo><msub><mi>R</mi><mn>2</mn></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>R</mi><mi>i</mi></msub><mo>)</mo></mrow></mtd><mtd><mi>i</mi><mo>=</mo><mn>1,2,3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>P</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>R</mi><mn>1</mn></msub><mo>+</mo><msub><mi>R</mi><mn>2</mn></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>R</mi><mrow><mi>P</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><mfrac><msubsup><mi>M</mi><mi>P</mi><mo>&prime;</mo></msubsup><msub><mi>F</mi><mn>0</mn></msub></mfrac></mtd><mtd><mi>i</mi><mo>=</mo><mi>P</mi></mtd></mtr></mtable><mo>.</mo></mfenced></mrow></math>]]></maths>
地址 210044 江苏省南京市浦口区宁六路219号