CN107146288A - 虚拟手术中支持实时按压形变的软组织模型建模方法 - Google Patents
虚拟手术中支持实时按压形变的软组织模型建模方法 Download PDFInfo
- Publication number
- CN107146288A CN107146288A CN201710344222.9A CN201710344222A CN107146288A CN 107146288 A CN107146288 A CN 107146288A CN 201710344222 A CN201710344222 A CN 201710344222A CN 107146288 A CN107146288 A CN 107146288A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- layer
- half elliptic
- soft tissue
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/28—Force feedback
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/44—Morphing
Landscapes
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明涉及虚拟手术中支持实时按压形变的软组织模型建模方法，包括如下步骤：对虚拟场景进行初始化；当检测到虚拟代理碰撞到虚拟软组织表面上任意一点时，在给定虚拟接触压力F作用下，虚拟代理与虚拟软组织交互的局部区域内部充满叶片弹簧力触觉计算模型，在交互过程中，软组织局部区域通过叶片弹簧力触觉计算模型进行变形计算和图形刷新，输出反馈为采用叶片弹簧力触觉计算模型计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息信号。本发明采用结构可靠的叶片弹簧，从而保证模型稳定性，方便力反馈和形变计算，不仅保证了建模的仿真性，而且叶片弹簧利用率高，变形计算速度快。
Description
技术领域
本发明涉及一种力触觉再现的建模方法，尤其涉及一种用于虚拟手术仿真过程中，人机交互的基于物理意义的虚拟手术中支持实时按压形变的软组织模型建模方法。
背景技术
如今二十一世纪科技迅猛发展，虚拟现实技术也随之崛起，虚拟手术应运而生，它是虚拟现实技术在医学上的一种重要应用。通过精确建模，来逼真地模拟人体组织在外力作用下的各种变形，并通过视觉显示和力反馈系统以及其它各种可能的感官形式，给用户提供真实的手术现场的感觉，让医生能够在计算机辅助的虚拟环境中进行手术训练和规划，根据医生需要提供不同的可以反复利用的实验对象，大大降低了外科医生手术训练的成本。虚拟手术能够模拟常见的手术操作过程，医生可在手术前对手术过程中涉及到的人体软组织进行虚拟手术仿真，减小手术风险，获得最好的手术效果。
虚拟手术基础中软组织仿真技术已成为当前的虚拟现实技术的重要组成部分。虚拟力触觉交互过程中，形变模型是手术成功的关键。准确快速的力触觉再现模型是虚拟手术成功的关键。采用基于物理意义的力触觉再现的建模方法，对软组织的变形进行实时仿真，已成为近几年研究热点。
现今常用的基于物理意义的柔性体变形模型中，有限元模型在建立软组织模型时只需少量模型参数，精确度高，但计算复杂，实时性有待提高。弹簧质点模型具有建模简单，计算速度较快，去除和增加点的操作较容易实现等优点，然而若质点周围约束不足，将会造成系统不稳定，若约束过多，则会减少变形范围，且精确度较低。边界元模型仅适用于简单、线性、各项同性软组织建模，对一些复杂、具有丰富细节的软组织建模失真较大。这些常用的软组织变形物理模型均存在计算较为繁杂和仿真精度不高等问题。
发明内容
本发明所要解决的技术问题是克服现有技术的不足，提供虚拟手术中支持实时按压形变的软组织模型建模方法，采用结构可靠的叶片弹簧，从而保证模型稳定性，并且方便力反馈和形变计算，不仅保证了建模的仿真性，而且叶片弹簧的利用率高，变形计算速度快。
为解决以上技术问题，本发明采取的一种技术方案是：
虚拟手术中支持实时按压形变的软组织模型建模方法，包括如下步骤：
步骤1：对虚拟场景进行初始化；
步骤2：当检测到虚拟代理碰撞到虚拟软组织表面上任意一点时，在给定虚拟接触压力F作用下，虚拟代理与虚拟软组织交互的局部区域内部充满叶片弹簧力触觉计算模型，在交互过程中，软组织局部区域通过叶片弹簧力触觉计算模型进行变形计算和图形刷新，输出反馈为采用叶片弹簧力触觉计算模型计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息信号，叶片弹簧力触觉计算模型的建模方法为：
(1)参数初始化，
(2)在给定虚拟接触压力F作用下，当虚拟代理碰撞到虚拟软组织表面上任意点时，在碰撞点下垂直悬挂一个伸直长度为L，叶片数为n，n大于等于3，厚度为a，宽度为b，弹性模量为E，弧高为H的半椭圆形对称式叶片弹簧，正中心下接另一相反方向相同规格的半椭圆形对称式叶片弹簧，叶片两端处相接触，以此构成软组织模型第一层，第一层下竖直悬挂第二层，第二层与第一层结构相似，半椭圆形对称式叶片弹簧伸直长度为2L，以此构成第二层，第三层也采用相似结构，半椭圆形对称式叶片弹簧伸直长度为3L，构成第三层，自第二层始，第i层的半椭圆形对称式叶片弹簧伸直长度都是第一层的i倍，而所有半椭圆形对称式叶片弹簧的叶片数为n，n大于等于3、厚度为a、宽度为b、弹性模量为E、弧高为H；
该半椭圆形对称式叶片弹簧的弹性模量E取决于人体虚拟软组织部位，且假设该部位软组织材质相同，即弹性模量E相同；
假设该半椭圆形对称式叶片弹簧受到给定虚拟接触压力F，且F的作用线和半椭圆形对称式叶片弹簧模型对称中心线一致，在虚拟接触压力F作用下，第一层先发生形变，第二层后发生形变，以此类推，假设在同一层的半椭圆形对称式叶片弹簧中，一层内上面的半椭圆形对称式叶片弹簧和下面的半椭圆形对称式叶片弹簧同时发生形变；当每层半椭圆形对称式叶片弹簧消耗的虚拟接触压力达到最大时，下一层半椭圆形对称式叶片弹簧才开始形变；假设在给定虚拟外力作用下，该层半椭圆形对称式叶片弹簧内的弹簧第i层弹簧消耗的虚拟接触压力未达到最大，则半椭圆形对称式叶片弹簧称为该层半椭圆形对称式叶片弹簧的变形截止层；
步骤3：确定半椭圆形对称式叶片弹簧的虚拟接触压力；
由于第一层中半椭圆形对称式叶片弹簧弧高为H，假设叶片弹簧可完全变形，则第一层半椭圆形对称式叶片弹簧的静挠度为H：
其中δ为挠度增大系数，为实际半椭圆形对称式叶片弹簧的挠度比理论截面梁挠度的增大系数；
则第一层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F1为：
第二层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F2为：
第三层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F3为：
则第i层半椭圆形对称式叶片弹簧消耗的虚拟接触压力Fi为：
确定该半椭圆形对称式叶片弹簧前i层共消耗的虚拟接触压力Fi′为：
步骤4：确定半椭圆形对称式叶片弹簧的形变量；
假设虚拟接触压力F满足Fi-1′＜F＜Fi′，Fi-1′、Fi′分别为前i层共消耗的虚拟接触压力和前i-1层共消耗的虚拟接触压力，第i层为变形截止层，则前i-1层的总变形量Hi-1′为：
Hi-1′＝(i-1)H
则第i层变形量为：
则软组织模型前i层的总变形量为：
优选地，半椭圆形对称式叶片弹簧前i层变形总计需要的时延时间不大于1ms。
优选地，在半椭圆形叶片弹簧一层中的上面半椭圆形叶片弹簧和下面半椭圆形叶片弹簧同时产生形变。
由于以上技术方案的采用，本发明与现有技术相比具有如下优点：
1、本发明和以往常用的基于物理意义的软组织变形仿真力触觉建模方法相比，该建模方法采用了结构可靠的叶片弹簧，从而保证模型的稳定性，并且方便力反馈和形变的计算，不仅保证了建模的仿真性，而且叶片弹簧的利用率高，变形计算速度快；
2、本发明每层的叶片弹簧的伸直长度、厚度、宽度、弹性模量、弧高、叶片数均相同，而叶片弹簧伸直长度不同，故每一层的受力不同中，与伸直长度的立方成正比，而每一层叶片弹簧的伸直长度为前一层的2倍，因此可以求得每一层的受力情况，计算较为简单；
3、本发明通过改变建模方法中第一层叶片弹簧的叶片数、厚度、宽度、弹性模量、弧高这些参数，就可对不同部位的虚拟软组织进行变形仿真，建模方法灵活、适用性广；
4、本发明建模方法用于在虚拟代理与虚拟柔性体交互的柔性体局部区域的按压变形计算、交互过程中，操作者感觉自然舒适、力触觉平稳、模拟效果逼真；
5、本发明可应用于虚拟外科手术仿真、深空探索、航空航天工业、军事应用、电子商务、远程医疗、导盲助残、虚拟游戏娱乐产业等领域，应用范围广。
附图说明
图1是本发明软组织变形仿真流程图；
图2是叶片弹簧力触觉虚拟模型的建模方法流程图；
图3是叶片弹簧力触觉虚拟模型示意图；
图4是叶片弹簧力触觉虚拟模型受力示意图。
具体实施方式
以下结合附图及具体实施例对本发明做进一步详细说明。
如附图1至附图4所示，本发明虚拟手术中支持实时按压形变的软组织模型建模方法，其特征在于：包括如下步骤：
步骤1：对虚拟场景进行初始化；
步骤2：当检测到虚拟代理碰撞到虚拟软组织表面上任意一点时，在给定虚拟接触压力F作用下，虚拟代理与虚拟软组织交互的局部区域内部充满叶片弹簧力触觉计算模型，在交互过程中，软组织局部区域通过叶片弹簧力触觉计算模型进行变形计算和图形刷新，输出反馈为采用叶片弹簧力触觉计算模型计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息信号，所述叶片弹簧力触觉计算模型建模方法为：
(1)参数初始化，
(2)在给定虚拟接触压力F作用下，当虚拟代理碰撞到虚拟软组织表面上任意点时，在碰撞点下垂直悬挂一个伸直长度为L，叶片数为n，n大于等于3，厚度为a，宽度为b，弹性模量为E，弧高为H的半椭圆形对称式叶片弹簧，正中心下接另一相反方向相同规格的半椭圆形对称式叶片弹簧，叶片两端处相接触，以此构成软组织模型第一层，第一层下竖直悬挂第二层，第二层与第一层结构相似，半椭圆形对称式叶片弹簧伸直长度为2L，以此构成第二层，第三层也采用相似结构，半椭圆形对称式叶片弹簧伸直长度为3L，构成第三层，自第二层始，第i层的半椭圆形对称式叶片弹簧伸直长度都是第一层的i倍，而所有半椭圆形对称式叶片弹簧的叶片数为n，n大于等于3、厚度为a、宽度为b、弹性模量为E、弧高为H；
该半椭圆形对称式叶片弹簧的弹性模量E取决于人体虚拟软组织部位，且假设该部位软组织材质相同，即弹性模量E相同；
假设该半椭圆形对称式叶片弹簧受到给定虚拟接触压力F，且F的作用线和半椭圆形对称式叶片弹簧模型对称中心线一致，在虚拟接触压力F作用下，第一层先发生形变，第二层后发生形变，以此类推，假设在同一层的半椭圆形对称式叶片弹簧中，一层内上面的半椭圆形对称式叶片弹簧和下面的半椭圆形对称式叶片弹簧同时发生形变；当每层半椭圆形对称式叶片弹簧消耗的虚拟接触压力达到最大时，下一层半椭圆形对称式叶片弹簧才开始形变；假设在给定虚拟外力作用下，该层半椭圆形对称式叶片弹簧内的弹簧第i层弹簧消耗的虚拟接触压力未达到最大，则半椭圆形对称式叶片弹簧称为该层半椭圆形对称式叶片弹簧的变形截止层；
步骤3：确定半椭圆形对称式叶片弹簧的虚拟接触压力；
由于第一层中半椭圆形对称式叶片弹簧弧高为H，假设叶片弹簧可完全变形，则第一层半椭圆形对称式叶片弹簧的静挠度为H：
其中δ为挠度增大系数，为实际半椭圆形对称式叶片弹簧的挠度比理论截面梁挠度的增大系数；
则第一层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F1为：
第二层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F2为：
第三层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F3为：
则第i层半椭圆形对称式叶片弹簧消耗的虚拟接触压力Fi为：
确定该半椭圆形对称式叶片弹簧前i层共消耗的虚拟接触压力Fi′为：
步骤4：确定半椭圆形对称式叶片弹簧的形变量；
假设虚拟接触压力F满足Fi-1′＜F＜Fi′，Fi-1′、Fi′分别为前i层共消耗的虚拟接触压力和前i-1层共消耗的虚拟接触压力，第i层为变形截止层，则前i-1层的总变形量Hi-1′为：
Hi-1′＝(i-1)H
则第i层变形量为：
则软组织模型前i层的总变形量为：
在本实施例中，半椭圆形对称式叶片弹簧前i层变形总计需要的时延时间不大于1ms。
在本实施例中，在半椭圆形叶片弹簧一层中的上面半椭圆形叶片弹簧和下面半椭圆形叶片弹簧同时产生形变。
以上对本发明做了详尽的描述，实施例的说明只是用于帮助理解本发明的方法及其核心思想，其目的在于让熟悉此领域技术的人士能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围之内。
Claims (3)
1.一种虚拟手术中支持实时按压形变的软组织模型建模方法，其特征在于：包括如下步骤：
步骤1：对虚拟场景进行初始化；
步骤2：当检测到虚拟代理碰撞到虚拟软组织表面上任意一点时，在给定虚拟接触压力F作用下，虚拟代理与虚拟软组织交互的局部区域内部充满叶片弹簧力触觉计算模型，在交互过程中，软组织局部区域通过叶片弹簧力触觉计算模型进行变形计算和图形刷新，输出反馈为采用叶片弹簧力触觉计算模型计算出来的反应在外力作用下软组织实时变形仿真的力触觉信息信号，所述叶片弹簧力触觉计算模型建模方法为：
(1)参数初始化，
(2)在给定虚拟接触压力F作用下，当虚拟代理碰撞到虚拟软组织表面上任意点时，在碰撞点下垂直悬挂一个伸直长度为L，叶片数为n，n大于等于3，厚度为a，宽度为b，弹性模量为E，弧高为H的半椭圆形对称式叶片弹簧，正中心下接另一相反方向相同规格的半椭圆形对称式叶片弹簧，叶片两端处相接触，以此构成软组织模型第一层，第一层下竖直悬挂第二层，第二层与第一层结构相似，半椭圆形对称式叶片弹簧伸直长度为2L，以此构成第二层，第三层也采用相似结构，半椭圆形对称式叶片弹簧伸直长度为3L，构成第三层，自第二层始，第i层的半椭圆形对称式叶片弹簧伸直长度都是第一层的i倍，而所有半椭圆形对称式叶片弹簧的叶片数为n，n大于等于3、厚度为a、宽度为b、弹性模量为E、弧高为H；
该半椭圆形对称式叶片弹簧的弹性模量E取决于人体虚拟软组织部位，且假设该部位软组织材质相同，即弹性模量E相同；
假设该半椭圆形对称式叶片弹簧受到给定虚拟接触压力F，且F的作用线和半椭圆形对称式叶片弹簧模型对称中心线一致，在虚拟接触压力F作用下，第一层先发生形变，第二层后发生形变，以此类推，假设在同一层的半椭圆形对称式叶片弹簧中，一层内上面的半椭圆形对称式叶片弹簧和下面的半椭圆形对称式叶片弹簧同时发生形变；当每层半椭圆形对称式叶片弹簧消耗的虚拟接触压力达到最大时，下一层半椭圆形对称式叶片弹簧才开始形变；假设在给定虚拟外力作用下，该层半椭圆形对称式叶片弹簧内的弹簧第i层弹簧消耗的虚拟接触压力未达到最大，则半椭圆形对称式叶片弹簧称为该层半椭圆形对称式叶片弹簧的变形截止层；
步骤3：确定半椭圆形对称式叶片弹簧的虚拟接触压力；
由于第一层中半椭圆形对称式叶片弹簧弧高为H，假设叶片弹簧可完全变形，则第一层半椭圆形对称式叶片弹簧的静挠度为H：
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>&delta;FL</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>Enba</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
其中δ为挠度增大系数，为实际半椭圆形对称式叶片弹簧的挠度比理论截面梁挠度的增大系数；
则第一层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F1为：
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>EHnba</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>&delta;L</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
第二层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F2为：
<mrow>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>EHnba</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mi>&delta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
第三层半椭圆形对称式叶片弹簧消耗的虚拟接触压力F3为：
<mrow>
<msub>
<mi>F</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>EHnba</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mi>&delta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
则第i层半椭圆形对称式叶片弹簧消耗的虚拟接触压力Fi为：
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>EHnba</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mi>&delta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
确定该半椭圆形对称式叶片弹簧前i层共消耗的虚拟接触压力Fi′为：
<mrow>
<msubsup>
<mi>F</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mi>&Sigma;</mi>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>EHnba</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mi>&delta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
步骤4：确定半椭圆形对称式叶片弹簧的形变量；
假设虚拟接触压力F满足Fi-1′＜F＜Fi′，Fi-1′、Fi′分别为前i层共消耗的虚拟接触压力和前i-1层共消耗的虚拟接触压力，第i层为变形截止层，则前i-1层的总变形量Hi-1′为：
Hi-1′＝(i-1)H
则第i层变形量为：
<mrow>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>F</mi>
<mo>-</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>L</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>Enba</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
则软组织模型前i层的总变形量为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msup>
<msub>
<mi>H</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>H</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>F</mi>
<mo>-</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>L</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>Enba</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>H</mi>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>&delta;L</mi>
<mn>3</mn>
</msup>
<mrow>
<mo>{</mo>
<mrow>
<mi>F</mi>
<mo>-</mo>
<mo>&Sigma;</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>EHnba</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mi>&delta;</mi>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>L</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>Enba</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
2.根据权利要求1所述的虚拟手术中支持实时按压形变的软组织模型建模方法，其特征在于：所述半椭圆形对称式叶片弹簧前i层变形总计需要的时延时间不大于1ms。
3.根据权利要求1所述的虚拟手术中支持实时按压形变的软组织模型建模方法，其特征在于：在所述半椭圆形叶片弹簧一层中的上面半椭圆形叶片弹簧和下面半椭圆形叶片弹簧同时产生形变。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710344222.9A CN107146288A (zh) | 2017-05-16 | 2017-05-16 | 虚拟手术中支持实时按压形变的软组织模型建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710344222.9A CN107146288A (zh) | 2017-05-16 | 2017-05-16 | 虚拟手术中支持实时按压形变的软组织模型建模方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107146288A true CN107146288A (zh) | 2017-09-08 |
Family
ID=59778157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710344222.9A Pending CN107146288A (zh) | 2017-05-16 | 2017-05-16 | 虚拟手术中支持实时按压形变的软组织模型建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107146288A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108550180A (zh) * | 2018-03-09 | 2018-09-18 | 南京信息工程大学 | 基于内点集域约束及高斯过程参数优化的血管建模方法 |
| CN108920875A (zh) * | 2018-07-29 | 2018-11-30 | 南京信息工程大学 | 一种基于磁性氧化模型的虚拟切割算法 |
| CN113343513A (zh) * | 2021-05-11 | 2021-09-03 | 南京信息工程大学 | 一种用于模拟软组织形变和路径切割的方法及装置 |
| CN117950504A (zh) * | 2024-03-26 | 2024-04-30 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976299A (zh) * | 2010-09-27 | 2011-02-16 | 南京信息工程大学 | 柔性力触觉再现的对称式板弹簧虚拟模型的建模方法 |
| WO2015102484A1 (en) * | 2013-12-31 | 2015-07-09 | Vrest B.V. | Method for generating a real-time haptic feedback signal for a haptic device of a virtual surgery simulator and virtual surgery simulator |
| CN106295085A (zh) * | 2016-10-17 | 2017-01-04 | 南京信息工程大学 | 柔性力触觉再现的叠合碟形弹簧虚拟模型的建模方法 |
| CN106528993A (zh) * | 2016-11-03 | 2017-03-22 | 南京信息工程大学 | 基于由碟形弹簧片构成的组合弹簧虚拟模型建模方法 |
| CN106560821A (zh) * | 2016-10-17 | 2017-04-12 | 南京信息工程大学 | 基于对合碟形弹簧的柔性体力触觉建模方法 |
-
2017
- 2017-05-16 CN CN201710344222.9A patent/CN107146288A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976299A (zh) * | 2010-09-27 | 2011-02-16 | 南京信息工程大学 | 柔性力触觉再现的对称式板弹簧虚拟模型的建模方法 |
| WO2015102484A1 (en) * | 2013-12-31 | 2015-07-09 | Vrest B.V. | Method for generating a real-time haptic feedback signal for a haptic device of a virtual surgery simulator and virtual surgery simulator |
| CN106295085A (zh) * | 2016-10-17 | 2017-01-04 | 南京信息工程大学 | 柔性力触觉再现的叠合碟形弹簧虚拟模型的建模方法 |
| CN106560821A (zh) * | 2016-10-17 | 2017-04-12 | 南京信息工程大学 | 基于对合碟形弹簧的柔性体力触觉建模方法 |
| CN106528993A (zh) * | 2016-11-03 | 2017-03-22 | 南京信息工程大学 | 基于由碟形弹簧片构成的组合弹簧虚拟模型建模方法 |
Non-Patent Citations (1)
| Title |
|---|
| 郑银环: "汽车钢板弹簧计算模型研究", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅱ辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108550180A (zh) * | 2018-03-09 | 2018-09-18 | 南京信息工程大学 | 基于内点集域约束及高斯过程参数优化的血管建模方法 |
| CN108550180B (zh) * | 2018-03-09 | 2021-11-19 | 南京信息工程大学 | 基于内点集域约束及高斯过程参数优化的血管建模方法 |
| CN108920875A (zh) * | 2018-07-29 | 2018-11-30 | 南京信息工程大学 | 一种基于磁性氧化模型的虚拟切割算法 |
| CN113343513A (zh) * | 2021-05-11 | 2021-09-03 | 南京信息工程大学 | 一种用于模拟软组织形变和路径切割的方法及装置 |
| CN113343513B (zh) * | 2021-05-11 | 2023-05-23 | 南京信息工程大学 | 一种用于模拟软组织形变和路径切割的方法及装置 |
| CN117950504A (zh) * | 2024-03-26 | 2024-04-30 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
| CN117950504B (zh) * | 2024-03-26 | 2024-05-31 | 北京航空航天大学 | 一种面向实时虚拟手术的软组织按压与夹取仿真方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107146288A (zh) | 虚拟手术中支持实时按压形变的软组织模型建模方法 | |
| CN103400023B (zh) | 软组织形变仿真方法 | |
| CN101964023B (zh) | 柔性力触觉再现的非对称形板弹簧虚拟模型的建模方法 | |
| CN102968811B (zh) | 微创血管介入手术中导丝的实时运动仿真方法 | |
| US20190108300A1 (en) | Methods for realistic and efficient simulation of moving objects | |
| CN203673431U (zh) | 一种运动轨迹虚拟装置 | |
| CN107992672A (zh) | 一种基于虚拟弹簧的软组织形变建模方法 | |
| US20220151701A1 (en) | Methods for realistic and efficient simulation of moving objects | |
| CN105302972A (zh) | 一种基于元球模型的软组织形变方法 | |
| CN109344500A (zh) | 一种非弹性质量弹簧模型及改进欧拉算法的软组织变形模拟方法 | |
| US20230061175A1 (en) | Real-Time Simulation of Elastic Body | |
| CN106485000A (zh) | 基于虚拟弹簧和球面调和函数力反馈模型建模方法 | |
| CN104063057B (zh) | 一种支持增强力触觉反馈的涡卷弹簧仿真系统 | |
| CN106528993A (zh) | 基于由碟形弹簧片构成的组合弹簧虚拟模型建模方法 | |
| CN107080586B (zh) | 虚拟手术中用于软组织实时形变的交叉弹簧模型 | |
| CN104268305A (zh) | 一种基于位置动力学的角度弯曲模型的仿真方法 | |
| CN101425188A (zh) | 带刚体核的广义弹簧振子形变仿真方法 | |
| CN103699753B (zh) | 一种支持增强力触觉反馈的多股螺旋弹簧虚拟模型 | |
| CN106295085A (zh) | 柔性力触觉再现的叠合碟形弹簧虚拟模型的建模方法 | |
| CN106485028B (zh) | 支持实时按压的柔性体变形仿真建模方法 | |
| Volino et al. | Stop-and-go cloth draping | |
| KR102273317B1 (ko) | 의상 패턴들의 오토 그레이딩 방법 및 장치 | |
| CN103699741B (zh) | 一种增强模拟柔性体可旋转变形的发条弹簧模型 | |
| CN110232745A (zh) | 一种基于可变形椭球体模型的人体建模方法 | |
| Tsiknis | Better cloth through unbiased strain limiting and physics-aware subdivision |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170908 |
|
| RJ01 | Rejection of invention patent application after publication |