CN114098837A - 基于光干涉及长光栅啁啾效应的力-形自感知穿刺针 - Google Patents
基于光干涉及长光栅啁啾效应的力-形自感知穿刺针 Download PDFInfo
- Publication number
- CN114098837A CN114098837A CN202111334357.XA CN202111334357A CN114098837A CN 114098837 A CN114098837 A CN 114098837A CN 202111334357 A CN202111334357 A CN 202111334357A CN 114098837 A CN114098837 A CN 114098837A
- Authority
- CN
- China
- Prior art keywords
- force
- grating
- fiber grating
- cavity
- needle
- 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
- 230000000694 effects Effects 0.000 title claims abstract description 18
- 239000000835 fiber Substances 0.000 claims abstract description 106
- 239000013307 optical fiber Substances 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000001228 spectrum Methods 0.000 claims description 27
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000013178 mathematical model Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 7
- KRHYYFGTRYWZRS-UHFFFAOYSA-N Fluorane Chemical compound F KRHYYFGTRYWZRS-UHFFFAOYSA-N 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 6
- 230000007935 neutral effect Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 4
- 238000003780 insertion Methods 0.000 claims description 4
- 230000037431 insertion Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 239000011347 resin Substances 0.000 claims description 4
- 229920005989 resin Polymers 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 239000003292 glue Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000002310 reflectometry Methods 0.000 claims description 3
- 238000003466 welding Methods 0.000 claims description 3
- 238000005530 etching Methods 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000001574 biopsy Methods 0.000 abstract description 2
- 230000004660 morphological change Effects 0.000 abstract description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
- A61B10/02—Instruments for taking cell samples or for biopsy
- A61B10/0233—Pointed or sharp biopsy instruments
- A61B10/0266—Pointed or sharp biopsy instruments means for severing sample
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B17/34—Trocars; Puncturing needles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B17/34—Trocars; Puncturing needles
- A61B17/3494—Trocars; Puncturing needles with safety means for protection against accidental cutting or pricking, e.g. limiting insertion depth, pressure sensors
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optical Transform (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于光干涉及长光栅啁啾效应的力‑形自感知穿刺针,包括针体,针体远离针尖的一端设有定位平面,针体侧壁上均匀设置有Ng个U形槽,U形槽内分别设置有第一至第Ng光纤,其中Ng≥3,每根光纤上均设置有一段啁啾光纤光栅与F‑P腔,啁啾光纤光栅设置于远离针尖一端,F‑P腔设置于靠近针尖一端,且啁啾光纤光栅与F‑P腔之间相隔一定距离。本发明实现了穿刺活检手术过程中穿刺针整体微量形态变化和尖端受力情况的精准监测。
Description
技术领域
本发明涉及医疗器械领域,具体涉及一种基于光干涉及长光栅啁啾效应的力-形自感知穿刺针。
背景技术
在临床诊断过程中,仅仅通过影像学检查难以对病灶定性,活体组织的病理检查才是诊断的金标准。经皮穿刺手术在临床操作时需要将穿刺针经过皮肤入针,经由人体组织到达病灶目标点进行取材。由于人体组织和器官生理结构错综复杂,给穿刺手术的精准入针带来极大困难,因此外科医生对实时监测穿刺针变形和针尖与人体组织的交互作用力情况具有迫切需求。目前穿刺针变形监测常用的技术方法有超声成像、计算机断层扫描以及磁共振成像。超声成像可以实时使用,但组织对比度低;计算机断层扫描需要高剂量的辐射,不能实时使用;核磁共振成像速度慢,也不能实时应用,且上述三种方法空间分辨率相对较低,难以准确监测针尖与穿刺路径间的微量偏差。目前用于交互作用力监测的传感器主要有压阻式、电容式、光电式等,以上传感器具有动态范围大,精度高等优势,但因为其单元体积较大,容易受到电磁干扰,处理电路复杂等问题,在穿刺针等小型医疗设备上应用受限。近年来,光纤光栅传感器因具有抗电磁干扰、良好的生物兼容性和微型化尺寸等优势而引起国内外学者的广泛关注,它能克服电类敏感元件的大部分局限,实现狭小空间下针体变形的精准测量。但是目前还只是通过光纤光栅分布式布置方式来进行变形监测,该种方式传感器的测量精度受光纤光栅位置及数量的影响,需要大量仿真或实验来优化传感器配置方案,设计过程需要消耗更多时间成本。
发明内容
本发明的目的在于提供一种基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,以实现穿刺活检手术过程中穿刺针整体微量形态变化和尖端受力情况的精准监测。
为解决上述技术问题,本发明提供了一种技术方案:基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,包括针体,针体远离针尖的一端设有定位平面,针体侧壁上均匀设置有Ng个U形槽,U形槽内分别设置有第一第Ng光纤,其中Ng≥3,每根光纤上均设置有一段啁啾光纤光栅与F-P腔,啁啾光纤光栅设置于远离针尖一端,F-P腔设置于靠近针尖一端,且啁啾光纤光栅与F-P腔之间相隔一定距离。
按上述方案,所述啁啾光纤光栅为线性纵向啁啾光纤光栅,其啁啾率为常数。
按上述方案,所述F-P腔为一单模光纤通过氢氟酸溶液腐蚀形成凹坑而后与另一单模光纤通过激光焊接而成。
按上述方案,所述第一至第Ng光纤在紧绷状态下通过树脂胶固定于U形槽内。
按上述方案,所述针体与国标16号穿刺针尺寸相同。
利用上文所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针实现的穿刺过程中三维力-形感知方法,其中:
受力感知方法为:
针体所受接触力分为轴向力Fz和横向力Fl,其中Fl在针体横截面上分解为沿x轴方向的力Fx和沿y轴方向的力Fy;接触力使三根光纤上的F-P腔的腔长发生不同程度的轴向变形和横向变形,F-P腔的腔长总变形分解为轴向变形和横向变形的线性叠加,通过分析针体的力学模型,得到在不同力的作用下各F-P腔的腔长变化量,进而得到此变化量与三维力的解耦矩阵,实现针尖处三维力大小和方向的测量;
形变感知方法为:
针体在受到轴向线性应变时,经解调得到的啁啾光纤光栅的反射光谱具有对称性,其峰值反射率对应的最大波长和线性应变函数的常系数存在映射关系,通过确定该波长值可计算得出线性应变函数表达式,结合悬臂梁模型的约束条件,对线性应变函数表达式进行两次积分得到针体的空间形变曲线。
按上述方案,所述受力感知方法具体为:
利用双解调算法将解调后形成的反射光谱分为F-P腔的干涉光谱与啁啾光纤光栅的宽带光谱;根据F-P腔干涉光谱,利用相位解调技术得出F-P腔的腔长变化量,得到腔长变化量与所受接触力之间的数学模型;随后得到全尺度啁啾光纤光栅的反射光谱与针体形变之间的数学模型,结合上述数学模型得到反射光谱信息与所受接触力的关系;
当力F作用在针尖处F-P腔上时,设第i个F-P腔腔长变化量δli(F)是由轴向受力形变δli,l和横向受力形变δli,z的线性叠加而成;
首先分析F-P腔轴向受力的形变模型,针体受轴向力Fz时,所有F-P腔受力情况相同,理论上对应的腔长变化量相同;
其中A为穿刺针的横截面积,E为穿刺针的杨氏模量;
当Ng=3时,在横截面上建立坐标系,三个F-P腔的坐标分别为(xi,yi)(i=1,2,3),有:
di=xicosα-yisinα
综合轴向受力和横向受力分析,F-P腔腔长变化量与针尖处所受三维力的关系表示为:
令Fx=Flcosα和Fy=Flsinα,则有:
令上式等号右侧的第一个矩阵为A,则有δl=AF,通过对穿刺针的标定实验获得矩阵A,进而可得穿刺针受力F=A-1δl。
按上述方案,所述形变感知方法具体为:
由于啁啾光纤光栅的结构特征,将一个长为L的啁啾光纤光栅看作是由N段长度为的均匀光纤光栅级联而成,若要获得啁啾光纤光栅上非均匀应变分布,首先需要得到每一段均匀光纤光栅上的均匀应变值;第n小段均匀光纤光栅的坐标进而得出啁啾光纤光栅第n小段的初始周期Λn,0:
其中Λ0为啁啾光纤光栅的初始周期,z为沿啁啾光纤光栅长度方向的坐标,C为啁啾率,L为光栅长度;
忽略所用光纤材料的吸收效应,根据耦合模理论和F矩阵理论,上述均匀光栅的传输特性表示为:
其中Ei(0)、Er(0)、Ei(L)、Er(L)分别为前后向传导模电场幅度,f11、f12、f21、f22分别为传输矩阵的四个矩阵元;为光纤光栅的输入矩阵,为光纤光栅的输出矩阵,F为光纤光栅的传输矩阵;其中为前后向传导模的耦合系数,λ为入射波长;为波数,为波数失谐量;φ为光栅相位;
由此,长为L的啁啾光纤光栅的传输矩阵表示为:
根据光纤光栅的边界条件:Ei(0)=0,Er(L)=0,得到啁啾光纤光栅反射光谱;
Ng=3时,当啁啾光纤光栅受到线性应变,结合针尖处应变为零的约束条件,设应变函数表达式为:
光纤光栅受到应力后的周期Λn,ε为:
其中pe为光纤纤芯的有效弹光系数;
通过上式通过针体上三根啁啾光纤光栅的宽带反射光谱得出针体受线性应变的函数表达,理论上三者的空间形态相同,所以其应变函数的系数值也相同,则针体中心轴的应变函数的系数认为是上述三个系数值的均值;
本发明的有益效果是:通过设置的啁啾光纤光栅和F-P腔实现了对针体形变和受力的精准感知;同时,设置的啁啾光纤光栅对针体全尺度的变形监测,相较于分布式监测,避免了复杂的试验和仿真分析优化。
附图说明
图1为本发明一实施例的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针整体结构图与光纤放大图;
图2为本发明一实施例的光纤布局示意图;
图3为本发明一实施例的针体横向与周向受力力学模型示意图;
图4为本发明一实施例的针体横截面参数示意图。
图中:1-定位平面,2-U形槽,3-树脂胶,4-啁啾光纤光栅,5-F-P腔,401-第一啁啾光纤光栅,402-第二啁啾光纤光栅,403-第三啁啾光纤光栅,501-第一F-P腔,502-第二F-P腔,503-第三F-P腔。
具体实施方式
为使本公开实施例的目的、技术方案和优点更加清楚,下面将结合本公开实施例的附图,对本公开实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本公开的一部分实施例,而不是全部的实施例。基于所描述的本公开的实施例,本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例,都属于本公开保护的范围。
参见图1,图2,基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,包括针体,针体远离针尖的一端设有定位平面1,针体侧壁上间隔120°均匀设置有三个U形槽2,U形槽2内分别设置有第一、第二和第三光纤,第一光纤上设置有第一啁啾光纤光栅401和第一F-P腔501,第而光纤上设置有第二啁啾光纤光栅402和第二F-P腔502,第三光纤上设置有第三啁啾光纤光栅403和第三F-P腔503,啁啾光纤光栅4设置于远离针尖一端,F-P腔5设置于靠近针尖一端,且啁啾光纤光栅4与F-P腔5之间相隔一定距离。
按上述方案,所述啁啾光纤光栅4为线性纵向啁啾光纤光栅,其啁啾率为常数。
按上述方案,所述F-P腔5为一单模光纤通过氢氟酸溶液腐蚀形成凹坑而后与另一单模光纤通过激光焊接而成。
按上述方案,所述第一、第二和第三光纤在紧绷状态下通过树脂胶3固定于U形槽2内。
按上述方案,所述针体与国标16号穿刺针尺寸相同。
利用上文所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针实现的穿刺过程中三维力-形感知方法,其中:
受力感知方法为:
针体所受接触力分为轴向力Fz和横向力Fl,其中Fl在针体横截面上分解为沿x轴方向的力Fx和沿y轴方向的力Fy;接触力使三根光纤上的F-P腔5的腔长发生不同程度的轴向变形和横向变形,F-P腔5的腔长总变形分解为轴向变形和横向变形的线性叠加,通过分析针体的力学模型,得到在不同力的作用下各F-P腔5的腔长变化量,进而得到此变化量与三维力的解耦矩阵,实现针尖处三维力大小和方向的测量;
形变感知方法为:
针体在受到轴向线性应变时,经解调得到的啁啾光纤光栅4的反射光谱具有对称性,其峰值反射率对应的最大波长和线性应变函数的常系数存在映射关系,通过确定该波长值可计算得出线性应变函数表达式,结合悬臂梁模型的约束条件,对线性应变函数表达式进行两次积分得到针体的空间形变曲线。
按上述方案,所述受力感知方法具体为:
利用双解调算法将解调后形成的反射光谱分为F-P腔5的干涉光谱与啁啾光纤光栅4的宽带光谱;根据F-P腔5干涉光谱,利用相位解调技术得出F-P腔5的腔长变化量,得到腔长变化量与所受接触力之间的数学模型;随后得到全尺度啁啾光纤光栅4的反射光谱与针体形变之间的数学模型,结合上述数学模型得到反射光谱信息与所受接触力的关系;
当力F作用在针尖处F-P腔5上时,设第i个F-P腔5腔长变化量δli(F)是由轴向受力形变δli,l和横向受力形变δli,z的线性叠加而成;
首先分析F-P腔5轴向受力的形变模型,针体受轴向力Fz时,所有F-P腔5受力情况相同,理论上对应的腔长变化量相同;
其中A为穿刺针的横截面积,E为穿刺针的杨氏模量;
在横截面上建立坐标系,三个F-P腔5的坐标分别为(xi,yi)(i=1,2,3),有:
di=xicosα-yisinα
综合轴向受力和横向受力分析,F-P腔5腔长变化量与针尖处所受三维力的关系表示为:
令Fx=Flcosα和Fy=Flsinα,则有:
令上式等号右侧的第一个矩阵为A,则有δl=AF,通过对穿刺针的标定实验获得矩阵A,进而可得穿刺针受力F=A-1δl。
按上述方案,所述形变感知方法具体为:
由于啁啾光纤光栅4的结构特征,将一个长为L的啁啾光纤光栅4看作是由N段长度为的均匀光纤光栅级联而成,若要获得啁啾光纤光栅4上非均匀应变分布,首先需要得到每一段均匀光纤光栅上的均匀应变值;第n小段均匀光纤光栅的坐标进而得出啁啾光纤光栅4第n小段的初始周期Λn,0:
其中Λ0为啁啾光纤光栅的初始周期,z为沿啁啾光纤光栅长度方向的坐标,C为啁啾率,L为光栅长度;
忽略所用光纤材料的吸收效应,根据耦合模理论和F矩阵理论,上述均匀光栅的传输特性表示为:
其中Ei(0)、Er(0)、Ei(L)、Er(L)分别为前后向传导模电场幅度,f11、f12、f21、f22分别为传输矩阵的四个矩阵元;为光纤光栅的输入矩阵,为光纤光栅的输出矩阵,F为光纤光栅的传输矩阵;其中为前后向传导模的耦合系数,λ为入射波长;为波数,为波数失谐量;φ为光栅相位;
由此,长为L的啁啾光纤光栅4的传输矩阵表示为:
根据光纤光栅的边界条件:Ei(0)=0,Er(L)=0,得到啁啾光纤光栅4反射光谱;
当啁啾光纤光栅4受到线性应变,结合针尖处应变为零的约束条件,设应变函数表达式为:
光纤光栅受到应力后的周期Λn,ε为:
其中pe为光纤纤芯的有效弹光系数;
通过上式通过针体上三根啁啾光纤光栅4的宽带反射光谱得出针体受线性应变的函数表达,理论上三者的空间形态相同,所以其应变函数的系数值也相同,则针体中心轴的应变函数的系数认为是上述三个系数值的均值;
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (8)
1.基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,其特征在于:包括针体,针体远离针尖的一端设有定位平面,针体侧壁上均匀设置有Ng个U形槽,U形槽内分别设置有第一至第Ng光纤,其中Ng≥3,每根光纤上均设置有一段啁啾光纤光栅与F-P腔,啁啾光纤光栅设置于远离针尖一端,F-P腔设置于靠近针尖一端,且啁啾光纤光栅与F-P腔之间相隔一定距离。
2.根据权利要求1所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,其特征在于:所述啁啾光纤光栅为线性纵向啁啾光纤光栅,其啁啾率为常数。
3.根据权利要求1所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,其特征在于:所述F-P腔为一单模光纤通过氢氟酸溶液腐蚀形成凹坑而后与另一单模光纤通过激光焊接而成。
4.根据权利要求1所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,其特征在于:所述第一至第Ng光纤在紧绷状态下通过树脂胶固定于U形槽内。
5.根据权利要求1所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针,其特征在于:所述针体与国标16号穿刺针尺寸相同。
6.利用权利要求1-5任一所述的基于光干涉及长光栅啁啾效应的力-形自感知穿刺针实现的穿刺过程中三维力-形感知方法,其特征在于:
受力感知方法为:
针体所受接触力分为轴向力Fz和横向力Fl,其中Fl在针体横截面上分解为沿x轴方向的力Fx和沿y轴方向的力Fy;接触力使三根光纤上的F-P腔的腔长发生不同程度的轴向变形和横向变形,F-P腔的腔长总变形分解为轴向变形和横向变形的线性叠加,通过分析针体的力学模型,得到在不同力的作用下各F-P腔的腔长变化量,进而得到此变化量与三维力的解耦矩阵,实现针尖处三维力大小和方向的测量;
形变感知方法为:
针体在受到轴向线性应变时,经解调得到的啁啾光纤光栅的反射光谱具有对称性,其峰值反射率对应的最大波长和线性应变函数的常系数存在映射关系,通过确定该波长值可计算得出线性应变函数表达式,结合悬臂梁模型的约束条件,对线性应变函数表达式进行两次积分得到针体的空间形变曲线。
7.根据权利要求6所述的穿刺过程中三维力-形感知方法,其特征在于:所述受力感知方法具体为:
利用双解调算法将解调后形成的反射光谱分为F-P腔的干涉光谱与啁啾光纤光栅的宽带光谱;根据F-P腔干涉光谱,利用相位解调技术得出F-P腔的腔长变化量,得到腔长变化量与所受接触力之间的数学模型;随后得到全尺度啁啾光纤光栅的反射光谱与针体形变之间的数学模型,结合上述数学模型得到反射光谱信息与所受接触力的关系;
当力F作用在针尖处F-P腔上时,设第i个F-P腔腔长变化量δli(F)是由轴向受力形变δli,l和横向受力形变δli,z的线性叠加而成;
首先分析F-P腔轴向受力的形变模型,针体受轴向力Fz时,所有F-P腔受力情况相同,理论上对应的腔长变化量相同;
其中A为穿刺针的横截面积,E为穿刺针的杨氏模量;
当Ng=3时,在横截面上建立坐标系,三个F-P腔的坐标分别为(xi,yi)(i=1,2,3),有:
di=xicosα-yisinα
综合轴向受力和横向受力分析,F-P腔腔长变化量与针尖处所受三维力的关系表示为:
令Fx=Flcosα和Fy=Flsinα,则有:
令上式等号右侧的第一个矩阵为A,则有δl=AF,通过对穿刺针的标定实验获得矩阵A,进而可得穿刺针受力F=A-1δl。
8.根据权利要求6所述的穿刺过程中三维力-形感知方法,其特征在于:所述形变感知方法具体为:
由于啁啾光纤光栅的结构特征,将一个长为L的啁啾光纤光栅看作是由N段长度为的均匀光纤光栅级联而成,若要获得啁啾光纤光栅上非均匀应变分布,首先需要得到每一段均匀光纤光栅上的均匀应变值;第n小段均匀光纤光栅的坐标进而得出啁啾光纤光栅第n小段的初始周期Λn,0:
其中Λ0为啁啾光纤光栅的初始周期,z为沿啁啾光纤光栅长度方向的坐标,C为啁啾率,L为光栅长度;
忽略所用光纤材料的吸收效应,根据耦合模理论和F矩阵理论,上述均匀光栅的传输特性表示为:
其中Ei(0)、Er(0)、Ei(L)、Er(L)分别为前后向传导模电场幅度,f11、f12、f21、f22分别为传输矩阵的四个矩阵元;为光纤光栅的输入矩阵,为光纤光栅的输出矩阵,F为光纤光栅的传输矩阵;其中为前后向传导模的耦合系数,λ为入射波长;为波数,为波数失谐量;φ为光栅相位;
由此,长为L的啁啾光纤光栅的传输矩阵表示为:
根据光纤光栅的边界条件:Ei(0)=0,Er(L)=0,得到啁啾光纤光栅反射光谱;
Ng=3时,当啁啾光纤光栅受到线性应变,结合针尖处应变为零的约束条件,设应变函数表达式为:
光纤光栅受到应力后的周期Λn,ε为:
其中pe为光纤纤芯的有效弹光系数;
通过上式通过针体上三根啁啾光纤光栅的宽带反射光谱得出针体受线性应变的函数表达,理论上三者的空间形态相同,所以其应变函数的系数值也相同,则针体中心轴的应变函数的系数认为是上述三个系数值的均值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111334357.XA CN114098837B (zh) | 2021-11-11 | 2021-11-11 | 基于光干涉及长光栅啁啾效应的力-形自感知穿刺针 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111334357.XA CN114098837B (zh) | 2021-11-11 | 2021-11-11 | 基于光干涉及长光栅啁啾效应的力-形自感知穿刺针 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114098837A true CN114098837A (zh) | 2022-03-01 |
CN114098837B CN114098837B (zh) | 2023-12-26 |
Family
ID=80378554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111334357.XA Active CN114098837B (zh) | 2021-11-11 | 2021-11-11 | 基于光干涉及长光栅啁啾效应的力-形自感知穿刺针 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114098837B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114795427A (zh) * | 2022-06-30 | 2022-07-29 | 北京微刀医疗科技有限公司 | 多方向电容力反馈穿刺针及穿刺设备 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070146699A1 (en) * | 2005-12-23 | 2007-06-28 | Bwt Property, Inc. | Raman spectroscopic apparatus utilizing internal grating stabilized semiconductor laser with high spectral brightness |
CN103153248A (zh) * | 2010-07-20 | 2013-06-12 | 约翰霍普金斯大学 | 用于外科手术仪器的干涉型力传感器 |
US8649847B1 (en) * | 2009-05-04 | 2014-02-11 | Intelligent Fiber Optic Systems, Inc. | Steerable shape sensing biopsy needle and catheter |
US20140053654A1 (en) * | 2012-08-22 | 2014-02-27 | U.S.A As Represented By The Administrator Of The National Aeronautics And Space Administration | Shape Sensing Using a Multi-Core Optical Fiber Having an Arbitrary Initial Shape in the Presence of Extrinsic Forces |
US20160008026A1 (en) * | 2013-03-04 | 2016-01-14 | The Board Of Trustees Of The Leland Stanford Junior University | Apparatus and methods involving elongated-medical instrument for sensing tissue interaction forces |
CN108652590A (zh) * | 2018-05-18 | 2018-10-16 | 武汉理工大学 | 一种集成光纤传感多参量测量的oct成像复合微探针及其制备方法 |
CN109765029A (zh) * | 2018-12-30 | 2019-05-17 | 北京信息科技大学 | 一种基于光纤光栅和f-p微腔的流场校测装置 |
US20190252845A1 (en) * | 2016-10-21 | 2019-08-15 | South China University Of Technology | Multi-wavelength narrow-linewidth single-frequency optical fiber laser source for laser radar system |
US20190247050A1 (en) * | 2006-11-21 | 2019-08-15 | David S. Goldsmith | Integrated system for the infixion and retrieval of implants |
CN112370124A (zh) * | 2020-11-12 | 2021-02-19 | 哈尔滨工业大学 | 一种用于机器人辅助穿刺手术的手动穿刺装置及系统 |
CN113208706A (zh) * | 2021-04-30 | 2021-08-06 | 哈尔滨工业大学 | 一种具有形状感知功能的可操纵穿刺针 |
-
2021
- 2021-11-11 CN CN202111334357.XA patent/CN114098837B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070146699A1 (en) * | 2005-12-23 | 2007-06-28 | Bwt Property, Inc. | Raman spectroscopic apparatus utilizing internal grating stabilized semiconductor laser with high spectral brightness |
US20190247050A1 (en) * | 2006-11-21 | 2019-08-15 | David S. Goldsmith | Integrated system for the infixion and retrieval of implants |
US8649847B1 (en) * | 2009-05-04 | 2014-02-11 | Intelligent Fiber Optic Systems, Inc. | Steerable shape sensing biopsy needle and catheter |
CN103153248A (zh) * | 2010-07-20 | 2013-06-12 | 约翰霍普金斯大学 | 用于外科手术仪器的干涉型力传感器 |
US20140053654A1 (en) * | 2012-08-22 | 2014-02-27 | U.S.A As Represented By The Administrator Of The National Aeronautics And Space Administration | Shape Sensing Using a Multi-Core Optical Fiber Having an Arbitrary Initial Shape in the Presence of Extrinsic Forces |
US20160008026A1 (en) * | 2013-03-04 | 2016-01-14 | The Board Of Trustees Of The Leland Stanford Junior University | Apparatus and methods involving elongated-medical instrument for sensing tissue interaction forces |
US20190252845A1 (en) * | 2016-10-21 | 2019-08-15 | South China University Of Technology | Multi-wavelength narrow-linewidth single-frequency optical fiber laser source for laser radar system |
CN108652590A (zh) * | 2018-05-18 | 2018-10-16 | 武汉理工大学 | 一种集成光纤传感多参量测量的oct成像复合微探针及其制备方法 |
CN109765029A (zh) * | 2018-12-30 | 2019-05-17 | 北京信息科技大学 | 一种基于光纤光栅和f-p微腔的流场校测装置 |
CN112370124A (zh) * | 2020-11-12 | 2021-02-19 | 哈尔滨工业大学 | 一种用于机器人辅助穿刺手术的手动穿刺装置及系统 |
CN113208706A (zh) * | 2021-04-30 | 2021-08-06 | 哈尔滨工业大学 | 一种具有形状感知功能的可操纵穿刺针 |
Non-Patent Citations (3)
Title |
---|
FOUZIA KHAN,等: "Pose Measurement of Flexible Medical Instruments Using Fiber Bragg Gratings in Multi-Core Fiber", 《IEEE SENSORS JOURNAL》, vol. 20, no. 18, XP011804892, DOI: 10.1109/JSEN.2020.2993452 * |
JIN SEOB KIM,等: "Shape Determination During Needle Insertion With Curvature Measurements", 《2017 IEEE/RSJ INTERNATIONAL CONFERENCE ON INTELLIGENT ROBOTS AND SYSTEMS》 * |
李天梁,等: "受弯件上粘贴型光纤布拉格光栅的应变传递规律", 《光学精密工程》, vol. 23, no. 5 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114795427A (zh) * | 2022-06-30 | 2022-07-29 | 北京微刀医疗科技有限公司 | 多方向电容力反馈穿刺针及穿刺设备 |
CN114795427B (zh) * | 2022-06-30 | 2022-09-02 | 北京微刀医疗科技有限公司 | 多方向电容力反馈穿刺针及穿刺设备 |
Also Published As
Publication number | Publication date |
---|---|
CN114098837B (zh) | 2023-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Su et al. | Fiber-optic force sensors for MRI-guided interventions and rehabilitation: A review | |
Lv et al. | A high-precision and miniature fiber Bragg grating-based force sensor for tissue palpation during minimally invasive surgery | |
CN208621784U (zh) | 一种柔性介入式医用导管空间弯曲检测的光栅光波导器件 | |
RU2628638C2 (ru) | Информация о вводе и выводе медицинского устройства с использованием распределенного измерения температуры оптическим волокном | |
CA2726367C (en) | Catheter with strain gauge sensor | |
EP3324850B1 (en) | Fiber-optic realshape sensor for enhanced doppler measurement display | |
US20110130645A9 (en) | System and method for reconstruction of the human ear canal from optical coherence tomography scans | |
Li et al. | Disposable FBG-based tridirectional force/torque sensor for aspiration instruments in neurosurgery | |
CN108593161A (zh) | 一种基于光纤光栅的微创外科手术机器人三维力传感器 | |
Moon et al. | Fiber-Bragg-grating-based ultrathin shape sensors displaying single-channel sweeping for minimally invasive surgery | |
Su et al. | A miniature MRI-compatible fiber-optic force sensor utilizing fabry-perot interferometer | |
Li et al. | Reaction force mapping by 3-axis tactile sensing with arbitrary angles for tissue hard-inclusion localization | |
US20140240713A1 (en) | Apparatuses and methods for imaging inside a vessel | |
Deng et al. | A miniature triaxial fiber optic force sensor for flexible ureteroscopy | |
Uzun et al. | Optical force sensor with enhanced resolution for MRI guided biopsy | |
Roodsari et al. | Fabrication and characterization of a flexible FBG-based shape sensor using single-mode fibers | |
CN114098837B (zh) | 基于光干涉及长光栅啁啾效应的力-形自感知穿刺针 | |
Leitão et al. | Optical fiber sensors for central arterial pressure monitoring | |
Xu et al. | Temperature-insensitive fiber-optic contact force sensor for steerable catheters | |
Chen et al. | Fabrication and shape detection of a catheter using fiber Bragg grating | |
Bützer et al. | Design and Evaluation of a Fiber‐Optic Grip Force Sensor with Compliant 3D‐Printable Structure for (f) MRI Applications | |
CN113520617A (zh) | 一种无源式三维力传感测头及光学式力传感器 | |
Parent et al. | 3D shape tracking of minimally invasive medical instruments using optical frequency domain reflectometry | |
Karpiouk et al. | Integrated optical coherence tomography and multielement ultrasound transducer probe for shear wave elasticity imaging of moving tissues | |
CN115290243A (zh) | 一种用于微创手术的三维力传感器及测量系统 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |