一种针灸力反馈形变模型的建模方法技术

本发明专利技术公开了一种针灸力反馈形变模型的建模方法，首先重建软组织模型；然后计算用于记录软组织模型的全局形变信息的全局形变矩阵；接着计算针刺点形变前的初始坐标和针刺点坐标之间的距离，并据此判定针尖是否刺破表皮；如果针尖未刺破表皮，实时产生虚拟弹簧计算皮肤表层的力反馈；如果针尖已刺破表皮，则计算软组织模型的内部力反馈，同时对针刺后的皮肤进行复原控制。本发明专利技术通过虚拟弹簧和体元模型可以分别计算针在皮肤内部和外部的力反馈，在确保计算简单的同时确保高精度的变形模拟，在交互过程中，操作者可以感受到触觉信息的实时性和真实感。

Modeling method of acupuncture force feedback deformation model

The present invention discloses a kind of acupuncture force feedback deformation model modeling method, the first reconstruction of the soft tissue model; then calculate the global deformation matrix global deformation information recording soft tissue model; then calculate the initial coordinates and coordinate between the acupuncture point acupuncture point deformation before the distance, and then determine whether needle punctured skin; if the needle not pierce the skin, feedback real-time virtual spring force calculation the surface of the skin; if the tip has punctured the skin, soft tissue model feedback calculation of internal force, and the acupuncture skin after restoration control. The virtual spring and the body element model can feedback the needle in the skin of internal and external force are calculated, in order to ensure the simple calculation and to ensure high precision deformation simulation, in the interactive process, real-time and realistic operator can feel tactile information.

【技术实现步骤摘要】
一种针灸力反馈形变模型的建模方法
本专利技术涉及针灸模拟，尤其涉及一种针灸力反馈形变模型的建模方法。
技术介绍
针灸是传统中国养生治病的神奇医术，在现代养生医学中也有很广泛的应用。然而人体穴位比较难找，对于刚接触针灸的新手来说是一个很大的挑战，几乎无法实战练习，而扎错穴位也会对人体产生或大或小的影响，因此训练至关重要。传统的针灸培训成本高、风险高，虚拟针灸模拟系统可以帮助中医新手练习针灸，有效地解决这一问题。柔性体表面形变是虚拟触觉交互过程的关键。常用的物理模型中，弹簧-质点模型具有建模简单、计算灵活、速度快的优点，但该模型精度有限，稳定性差。有限元模型精度高，但计算相对复杂。因此，具有良好的交互功能和形变效果的柔体模型对于虚拟手术系统来说至关重要。
技术实现思路
本专利技术所要解决的技术问题是针对
技术介绍
的不足，提供了一种针灸力反馈形变模型的建模方法。本专利技术为解决上述技术问题采用以下技术方案：一种针灸力反馈形变模型的建模方法，包含以下步骤：步骤1)，重建软组织模型：采用体元模型重建软组织内部，采用网格模型重建软组织表皮；所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元；进行针灸模拟时，将模型针头在软组织内部所经过路径上的内部体元替换成针体元；步骤2)，根据以下公式计算用于记录软组织模型的全局形变信息的的全局形变矩阵A：其中，n是软组织模型的网格质点数，mi是软组织模型中质点i的质量，i为大于等于1小于等于n的自然数，pi是软组织模型中质点i的坐标，pcm是形变后软组织模型的质心坐标，是质点i的初始坐标，是初始状态下软组织模型的质心坐标；步骤3)，计算针刺点形变前的初始坐标和针刺点坐标pc之间的距离当且针尖坐标x与针刺点坐标pc重合时，判定皮肤表层为未被刺破的状态，执行步骤4)，其中，dmax为预设的皮肤最大形变量；当且针尖坐标x与针刺点坐标pc不重合时，判定皮肤表层为被刺破的状态，执行步骤5)；步骤4)，根据以下公式实时产生虚拟弹簧计算皮肤表层的力反馈Fout：其中，σ(i)为软组织模型中受形变影响的网格质点集合,k为预设的虚拟弹簧系数，为软组织模型中质点i的初始坐标和当前坐标pi之间的距离，Nout为以针尖为顶点沿针方向的单位向量；步骤5)，根据针刺深度计算软组织模型的内部力反馈，同时对针刺后的皮肤进行复原控制：步骤5.1)，根据针刺点坐标pc和针尖坐标x的距离din计算软组织模型的内部力反馈Fin：Fin＝ρ·din·Nin其中，ρ是预设的的皮下组织密度，din是针刺点坐标pc和针尖坐标x的距离，Nin是从针尖坐标x到针刺点坐标pc的单位向量；步骤5.2)，根据以下公式计算软组织模型中各个质点目标位置的坐标：其中，为质点i目标位置的坐标，R是旋转矩阵，S是对称系数，矩阵AT是矩阵A的转置矩阵；步骤5.3)，根据以下公式优化软组织模型中各质点目标位置的坐标：其中，ηi为预设的恢复系数，ηi∈[0,1]；步骤5.4)，对于软组织模型中的各个质点i，控制其在预设的速率Rr下以ηi＝ηi-Rr的速度向其目标位置移动。本专利技术可以在保证较好地表现针刺形变和恢复效果的前提下实现快速、逼真、实时的力反馈计算，与现有技术相比，具有以下技术效果：1.使用体元和网格模型建模从而保证模型的稳定性，并且方便力反馈和形变的计算。2.计算全局形变矩阵A以保证模型的全局稳定性，防止模型走样，并未之后的恢复算法做准备。3.计算目标位置gi实现皮肤被针刺破后的恢复，使针灸模型更为逼真。4.插入恢复系数ηi以控制表面的恢复速率，进一步模拟真实情景下的针灸形变，保证形变恢复的稳定性和真实性。5.由针头x和针刺点pc之间的位置以及针刺点pc和针刺点初始位置的距离为根据进行刺破判定，可以快速找到针头刺破皮肤的时间点，并以此为分割线进行表面形变的模拟和力反馈的模拟，逼真地展现在两种不同状态下的针灸模型。6.未刺破前使用虚拟弹簧计算外部力反馈，使用虚拟弹簧计算力反馈方便快捷，计算时通过形变位置各顶点与其初始位置间的距离的改变，保证了力反馈的刷新。7.刺破后考虑皮下组织密度，进行内部力反馈的计算，让模型的力反馈更为逼真，力反馈随针刺入的深度实时刷新，满足了系统的实时性。附图说明图1是针灸力反馈形变模型的仿真流程图；图2是改良的形状匹配算法示例图；图3是刺破前手臂横截面示例图；图4是刺破后手臂横截面示例图。具体实施方式下面结合附图对本专利技术的技术方案做进一步的详细说明：如图1所示，本专利技术公开了一种针灸力反馈形变模型的建模方法，包含以下步骤：步骤1)，采用体元模型重建软组织内部，采用网格模型重建软组织表皮；所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元；进行针灸模拟时，将模型针头在软组织内部所经过路径上的内部体元替换成针体元。步骤2)，计算全局形变矩阵A，全局形变矩阵A记录软组织的全局形变信息，也是形变恢复算法的基础，全局形变矩阵A需要获取软组织模型初始状态下的各质点坐标及质心坐标，在整个针灸模拟的过程中实时计算更新，根据以下公式计算软组织内部的全局形变矩阵A：其中，n是软组织模型的网格质点数，mi是软组织模型中质点i的质量，i为大于等于1小于等于n的自然数,pi是软组织模型中质点i的坐标，pcm是形变后软组织模型的质心坐标，是质点i的初始坐标，是初始状态下软组织模型的质心坐标。步骤3)，进行刺破判定，计算针刺点形变前的初始坐标和针刺点坐标pc之间的距离其中，dmax为预设的皮肤最大形变量；当且针尖坐标x与针刺点坐标pc重合时，视为皮肤表层未被刺破的状态，执行步骤4)；当且针尖坐标x与针刺点坐标pc不重合时，视为皮肤表层被刺破的状态，执行步骤5)。步骤4)，皮肤未刺破，根据皮肤形变实时产生虚拟弹簧计算皮肤表层的力反馈，如图3所示，皮肤表层的力反馈Fout计算如下：其中,σ(i)为软组织模型受形变影响的网格质点集合,k为预设的虚拟弹簧系数，为软组织模型中质点i初始坐标和当前坐标pi之间的距离，Nout为以针尖为顶点沿针方向的单位向量。步骤5)，皮肤已刺破，根据针刺深度计算软组织模型的内部力反馈，同时对针刺后的皮肤形变进行优化，使皮肤在刺破后一定程度地恢复原状，如图4所示。步骤5.1)，计算软组织模型的内部接触力，通过针刺点坐标pc和针尖坐标x的距离din的改变生成持续且变化的力反馈，内部力反馈Fin计算如下：Fin＝ρ·din·Nin其中，ρ是预设的皮下组织密度，din是针刺点坐标pc和针尖坐标x的距离，Nin是从针尖坐标x到针刺点坐标pc的单位向量。计算刺破后的皮肤形变，当皮肤组织被刺破后，皮肤会一定程度地恢复形变，表现为各质点i向其对应的目标位置gi移动，从而恢复初始形状，如图2所示。步骤5.2)，目标位置gi是皮肤形变恢复的关键，通过全局形变矩阵A求解目标位置gi的坐标根据以下公式计算各质点的目标位置gi的坐标其中R是旋转矩阵，S是对称系数，R和S都是求解目标位置gi的坐标所需要的重要参数，矩阵AT是矩阵A的转置矩阵；在实际操作时各质点对应的目标位置gi位置固定，各质点i会快速回复形变，造成模拟的失真，为了防止软组织表面过快恢复，得到逼真的针灸形变效果，插入恢复系数ηi，ηi∈[0,1]，如图4所示本文档来自技高网...

【技术保护点】
一种针灸力反馈形变模型的建模方法，其特征在于，包含以下步骤：步骤1)，重建软组织模型：采用体元模型重建软组织内部，采用网格模型重建软组织表皮；所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元；进行针灸模拟时，将模型针头在软组织内部所经过路径上的内部体元替换成针体元；步骤2)，根据以下公式计算用于记录软组织模型的全局形变信息的的全局形变矩阵A：

【技术特征摘要】
1.一种针灸力反馈形变模型的建模方法，其特征在于，包含以下步骤：步骤1)，重建软组织模型：采用体元模型重建软组织内部，采用网格模型重建软组织表皮；所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元；进行针灸模拟时，将模型针头在软组织内部所经过路径上的内部体元替换成针体元；步骤2)，根据以下公式计算用于记录软组织模型的全局形变信息的的全局形变矩阵A：其中，n是软组织模型的网格质点数，mi是软组织模型中质点i的质量，i为大于等于1小于等于n的自然数，pi是软组织模型中质点i的坐标，pcm是形变后软组织模型的质心坐标，是质点i的初始坐标，是初始状态下软组织模型的质心坐标；步骤3)，计算针刺点形变前的初始坐标和针刺点坐标pc之间的距离当且针尖坐标x与针刺点坐标pc重合时，判定皮肤表层为未被刺破的状态，执行步骤4)，其中，dmax为预设的皮肤最大形变量；当且针尖坐标x与针刺点坐标pc不重合时，判定皮肤表层为被刺破的状态，执行步骤5)；步骤4)，根据以下公式实时产生虚拟弹簧计算皮肤表层的力反馈Fout：

【专利技术属性】
技术研发人员：张小瑞，段佳骊，孙伟，王澎湃，夏薇薇，宋爱国，牛建伟，朱利丰，刘佳，徐慧，舒雅，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏,32