留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

多芯光纤光栅形状传感性能与重构误差研究

周勇 胡文彬 程普 叶泓蕤 郭东来 杨明红

周勇, 胡文彬, 程普, 叶泓蕤, 郭东来, 杨明红. 多芯光纤光栅形状传感性能与重构误差研究[J]. 红外与激光工程, 2023, 52(3): 20220485. doi: 10.3788/IRLA20220485
引用本文: 周勇, 胡文彬, 程普, 叶泓蕤, 郭东来, 杨明红. 多芯光纤光栅形状传感性能与重构误差研究[J]. 红外与激光工程, 2023, 52(3): 20220485. doi: 10.3788/IRLA20220485
Zhou Yong, Hu Wenbin, Cheng Pu, Ye Hongrui, Guo Donglai, Yang Minghong. Research on shape sensing performance and reconstruction error of multi-core fiber grating[J]. Infrared and Laser Engineering, 2023, 52(3): 20220485. doi: 10.3788/IRLA20220485
Citation: Zhou Yong, Hu Wenbin, Cheng Pu, Ye Hongrui, Guo Donglai, Yang Minghong. Research on shape sensing performance and reconstruction error of multi-core fiber grating[J]. Infrared and Laser Engineering, 2023, 52(3): 20220485. doi: 10.3788/IRLA20220485

多芯光纤光栅形状传感性能与重构误差研究

doi: 10.3788/IRLA20220485
基金项目: 湖北省技术创新专项重大项目 (2018AAA016)
详细信息
    作者简介:

    周勇,男,硕士生,主要研究方向为光纤传感技术

  • 中图分类号: TN253

Research on shape sensing performance and reconstruction error of multi-core fiber grating

Funds: Hubei Province Technology Innovation Special Major Project (2018AAA016)
  • 摘要: 多芯光纤光栅形状传感技术利用空分复用以及应变监测的优势,结合不同的栅点布设方案,实现待测对象的连续曲率和形状传感。首先介绍了多芯光纤光栅曲率和挠率传感原理,提出采用齐次矩阵变换的三维重构算法实现光纤的三维形状重构。为了探究不同光栅密度对实验精度的影响,利用算法编程模拟了不同光栅间距下的三维形状重构精度,依据模拟仿真的结果,建立了不同光栅间距与三维重构误差之间的关系。三维形状传感实验使用光栅间距为10 cm和5 cm的七芯光纤光栅串。实验结果表明,最大误差出现在尾点处,分别为2.56 cm和1.15 cm,占全长的3.2%和1.4%,平均误差为1.32 cm和0.62 cm,占全长的1.7%和0.8%。实验结果与仿真值比较接近,说明可以依据仿真结果对不同光栅间距下的三维形状误差进行预测。结合具体的应用场景合理配置测点资源,在较低的成本范围内实现高性能的检测。
  • 图  1  三维形状重建技术路线图

    Figure  1.  3D shape reconstruction technology flow chart

    图  2  七芯光纤传感原理示意图。(a)曲率;(b)挠率

    Figure  2.  Schematic diagram of the seven-core optical fiber sensing principle. (a) Curvature; (b) Torsion

    图  3  三维形状重构示意图

    Figure  3.  Schematic diagram of 3D shape reconstruction

    图  4  不同光栅间距的模拟结果。(a)三维重构效果图;(b)空间位置误差;(c)光栅间距与误差曲线斜率之间的关系

    Figure  4.  Simulation results of different grating spacings. (a) Diagram of 3D reconstruction effect; (b) Spatial position error; (c) Relationship between grating spacing and error curve slope

    图  5  三维形状传感实验装置。(a)系统简图;(b)七芯光纤端面;(c)光纤并带结构;(d)三维螺旋曲率槽

    Figure  5.  Experimental setup for 3D shape sensing. (a) Schematic diagram of the system; (b) End face of seven-core fiber; (c) Optical fiber with ribbon structure; (d) 3D helical curvature groove

    图  6  (a)各个纤芯的光栅光谱;(b) 光栅串的光谱

    Figure  6.  (a) Grating spectrum of each fiber core; (b) Spectrum of grating string

    图  7  三维形状传感实验。(a)重构曲线与真实曲线;(b)误差

    Figure  7.  3D shape sensing experiment. (a) Reconstructed curve and true curve; (b) Error

    表  1  实验值与模拟值对比

    Table  1.   Comparison of experimental and simulated values

    Grating pitch,
    L/cm
    Maximum error
    /cm
    Error percentageSimulation error percentage
    51.151.4%0.6%
    102.563.2%2.7%
    下载: 导出CSV
  • [1] Lai T, Cheng P, Yan C, et al. 2D and 3D shape sensing based on 7-core fiber Bragg gratings [J]. Photonic Sensors, 2020, 10(4): 306-315.
    [2] Hou Maoxiang, Yang Kaiming, He Jun, et al. Two-dimensional vector bending sensor based on seven-core fiber Bragg gratings [J]. Optics Express, 2018, 26(18): 23770-23781. doi:  10.1364/OE.26.023770
    [3] Luo Shujun, Wang Xiaoliang, Chen Daru, et al. Multi-core fiber and its application for bending sensor [J]. Acta Photonica Sinica, 2016, 45(2): 0206005. (in Chinese)
    [4] Roesthuis R J, Kemp M, van den Dobbelsteen J J, et al. Three-dimensional needle shape reconstruction using an array of fiber bragg grating sensors [J]. IEEE/ASME Transactions on Mechatronics, 2014, 19(4): 1115-1126. doi:  10.1109/TMECH.2013.2269836
    [5] Hoehn K, Olsson A, Arkwright J W. High capacity torque and compression measurements using fibre optic sensors [C]//2017 2nd International Conference for Fibre-optic and Photonic Sensors for Industrial and Safety Applications (OFSIS), IEEE, 2017: 39-44.
    [6] 机器人末端力测量方法研究[J]. 红外与激光工程, 2021, 50(07): 219-227.

    Tang Chen, He Yanlin, Zhu Lianqing, et al. Research on end force measurement method of soft robot based on fiber Bragg grating[J]. Infrared and Laser Engineering, 2021, 50(7): 20200386. (in Chinese)
    [7] Zhang Zhili, Fu Zhulin, Zhao Bing, et al. FBG sensing technology application in detection of bastion safety [J]. Infrared and Laser Engineering, 2011, 40(3): 492-496. (in Chinese) doi:  10.3969/j.issn.1007-2276.2011.03.023
    [8] Burgner-Kahrs J, Rucker D C, Choset H. Continuum robots for medical applications: A survey [J]. IEEE Transactions on Robotics, 2015, 31(6): 1261-1280. doi:  10.1109/TRO.2015.2489500
    [9] Swaney P J, York P A, Gilbert H B, et al. Design, fabrication, and testing of a needle-sized wrist for surgical instruments [J]. J Med Device, 2017, 11(1): 0145011-0145019.
    [10] Shi C, Luo X, Qi P, et al. Shape sensing techniques for continuum robots in minimally invasive surgery: A survey [J]. IEEE Transactions on Biomedical Engineering, 2017, 64(8): 1665-1678. doi:  10.1109/TBME.2016.2622361
    [11] Heiden M, Henken K R, Chen L K, et al. Accurate and efficient fiber optical shape sensor for MRI compatible minimally invasive instruments [C]//Proceedings of SPIE, 2012, 8550: 85500L.
    [12] Hoffman J, Waters D H, Khadka S, et al. Shape sensing of polymer core composite electrical transmission lines using FBG sensors [J]. IEEE Transactions on Instrumentation and Measurement, 2019, 69(1): 249-257.
    [13] Duncan R G, Froggatt M E, Kreger S T, et al. High-accuracy fiber-optic shape sensing [C]//Proceedings of SPIE, 2007, 6530: 65301S.
    [14] Li Zhengying, Sun Wenfeng, Wang Honghai. Research on the ultra-weak reflective fiber Bragg grating sensing technology based on optical frequency domain reflection technology [J]. Acta Optica Sinica, 2015, 35(8): 0806003. (in Chinese)
    [15] Khan F, Denasi A, Barrera D, et al. Multi-core optical fibers with Bragg gratings as shape sensor for flexible medical instruments [J]. IEEE Sensors Journal, 2019, 19(14): 5878-5884. doi:  10.1109/JSEN.2019.2905010
    [16] Al-Ahmad O, Ourak M, Roosbroeck J V, et al. Improved FBG-based shape sensing methods for vascular catheterization treatment [J]. IEEE Robotics and Automation Letters, 2020, 5: 4683-4694.
    [17] Zhao Shiyuan, Cui Jiwen, Chen Mengmeng. Review on optical fiber shape sensing technology [J]. Optics and Precision Engineering, 2020, 28(1): 10-29. (in Chinese)
    [18] 郭启全. 计算机图形学教程[M]. 北京: 机械工业出版社, 2003.
    [19] Paloschi D, Bronnikov K A, Korganbayev S, et al. 3D shape sensing with multicore optical fibers: Transformation matrices vs Frenet-Serret equations for real-time application [J]. IEEE Sensors Journal, 2020, 21(4): 4599-4609.
  • [1] 解宇恒, 裴丽, 何倩, 常彦彪, 郭智君, 王建帅, 郑晶晶, 宁提纲, 李晶.  多芯光纤折射率与内应力分布重构技术(特邀) . 红外与激光工程, 2022, 51(1): 20210758-1-20210758-6. doi: 10.3788/IRLA20210758
    [2] 刘炳锋, 董明利, 孙广开, 何彦霖, 祝连庆.  浮空器柔性复合蒙皮形变光纤布拉格光栅传感器应变传递特性 . 红外与激光工程, 2021, 50(5): 20200315-1-20200315-9. doi: 10.3788/IRLA20200315
    [3] 吕佳豪, 董明利, 何彦霖, 孙广开, 周康鹏.  引入曲率与角度校正的柔性机构三维形状多芯光纤重构方法 . 红外与激光工程, 2021, 50(5): 20200453-1-20200453-7. doi: 10.3788/IRLA20200453
    [4] 张学强, 孙博, 贾静.  光纤布拉格光栅温度传感增敏特性的实验研究 . 红外与激光工程, 2019, 48(11): 1118003-1118003(7). doi: 10.3788/IRLA201948.1118003
    [5] 张雄雄, 宋言明, 孙广开, 孟凡勇, 董明利.  柔性复合基体光纤布拉格光栅曲率传感器 . 红外与激光工程, 2019, 48(2): 222001-0222001(7). doi: 10.3788/IRLA201948.0222001
    [6] 刘畅, 王双, 梁应剑, 江俊峰, 梅运桥, 刘琨, 齐晓光, 李鑫, 李元耀, 刘铁根.  光纤法珀压力传感系统设计与风洞初步实验 . 红外与激光工程, 2018, 47(7): 722002-0722002(7). doi: 10.3788/IRLA201847.0722002
    [7] 裴丽, 王建帅, 郑晶晶, 宁提纲, 解宇恒, 何倩, 李晶.  空分复用光纤的特性及其应用研究 . 红外与激光工程, 2018, 47(10): 1002001-1002001(12). doi: 10.3788/IRLA201847.1002001
    [8] 吕珊珊, 耿湘宜, 张法业, 肖航, 姜明顺, 曹玉强, 隋青美.  基于菱形光纤布拉格光栅传感阵列的声发射定位技术 . 红外与激光工程, 2017, 46(12): 1222005-1222005(5). doi: 10.3788/IRLA201746.1222005
    [9] 鞠恒, 林成新, 张佳琪, 刘志杰.  Fe-Mn-Si形状记忆合金涂层残余应力模拟与测量 . 红外与激光工程, 2017, 46(10): 1017009-1017009(10). doi: 10.3788/IRLA201777.1017009
    [10] 廖天河, 刘伟, 高穹.  不同形状激光波束在大气中传输的湍流效应 . 红外与激光工程, 2015, 44(S1): 41-45.
    [11] 常金勇, 强希文, 胡月宏, 宗飞, 李志朝, 封双连.  激光传输光束抖动效应的数值模拟 . 红外与激光工程, 2015, 44(S1): 46-49.
    [12] 戴毅斌, 樊玉杰, 李明尧, 刘晓宇, 蒋彭胜, 殷开婷.  多点激光微冲击成形的数值模拟研究 . 红外与激光工程, 2015, 44(S1): 50-56.
    [13] 张金玉, 孟祥兵, 杨正伟, 王冬冬, 陶胜杰.  红外锁相法涂层测厚数值模拟与分析 . 红外与激光工程, 2015, 44(1): 6-11.
    [14] 周倩, 宁提纲, 温晓东, 李超.  一种双包层半径光纤布拉格光栅传感器 . 红外与激光工程, 2015, 44(3): 1024-1027.
    [15] 许兆美, 刘永志, 周建忠, 蒋素琴, 王庆安, 汪通, 洪宗海.  基于生死单元的激光铣削温度场数值模拟与验证 . 红外与激光工程, 2014, 43(6): 1755-1760.
    [16] 陈林, 杨立, 范春利, 吕事桂, 石宏臣.  线性调频激励的红外无损检测及其数值模拟 . 红外与激光工程, 2014, 43(5): 1385-1389.
    [17] 孙杜娟, 胡以华, 李乐.  气溶胶“沉降—扩散”联合动态数值模拟 . 红外与激光工程, 2014, 43(2): 449-453.
    [18] 潘洪刚, 童峥嵘, 张爱玲, 宋殿友, 薛玉明.  基于级联多模布拉格光栅和高精细度滤波器的可调谐双波长窄线宽掺铒光纤激光器 . 红外与激光工程, 2014, 43(12): 3912-3917.
    [19] 唐力铁, 于志闯, 赵乐至, 尹飞, 郭士波, 谈斌.  DF/HF化学激光器喷管总压损失的数值模拟 . 红外与激光工程, 2013, 42(5): 1194-1197.
    [20] 宋民青, 侯尚林, 张保侠, 黎锁平, 刘延君.  光子晶体光纤布拉格光栅慢光的研究 . 红外与激光工程, 2013, 42(6): 1547-1552.
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  248
  • HTML全文浏览量:  38
  • PDF下载量:  41
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-07-13
  • 修回日期:  2022-09-12
  • 刊出日期:  2023-03-25

多芯光纤光栅形状传感性能与重构误差研究

doi: 10.3788/IRLA20220485
    作者简介:

    周勇,男,硕士生,主要研究方向为光纤传感技术

基金项目:  湖北省技术创新专项重大项目 (2018AAA016)
  • 中图分类号: TN253

摘要: 多芯光纤光栅形状传感技术利用空分复用以及应变监测的优势,结合不同的栅点布设方案,实现待测对象的连续曲率和形状传感。首先介绍了多芯光纤光栅曲率和挠率传感原理,提出采用齐次矩阵变换的三维重构算法实现光纤的三维形状重构。为了探究不同光栅密度对实验精度的影响,利用算法编程模拟了不同光栅间距下的三维形状重构精度,依据模拟仿真的结果,建立了不同光栅间距与三维重构误差之间的关系。三维形状传感实验使用光栅间距为10 cm和5 cm的七芯光纤光栅串。实验结果表明,最大误差出现在尾点处,分别为2.56 cm和1.15 cm,占全长的3.2%和1.4%,平均误差为1.32 cm和0.62 cm,占全长的1.7%和0.8%。实验结果与仿真值比较接近,说明可以依据仿真结果对不同光栅间距下的三维形状误差进行预测。结合具体的应用场景合理配置测点资源,在较低的成本范围内实现高性能的检测。

English Abstract

    • 光纤光栅(Fiber Bragg Gratings, FBG)传感技术广泛应用于各种2D和3D形状传感中[1-4]。相较于传统的电学形状传感器,光纤光栅形状传感器具有体积小、质量轻、灵敏度高、抗电磁干扰等显著优点,使其在建筑,工业,医疗等领域有着不可替代的优势[5-12]。基于光栅的多芯光纤形状传感技术经过多年的发展也愈加成熟。目前,基于多芯光纤光栅的形状传感技术主要分为两类,一类是以LUNA公司为代表的基于全同弱光栅的光频域反射(Optical Frequency Domain Reflectometry, OFDR)解调技术。2007年,Duncan等[13]人采用OFDR技术对长度为110 cm的多芯光纤光栅阵列进行了二维和三维形状重构实验,实验误差分别为1.9 mm和5.8 mm。基于OFDR解调技术空间应变传感分辨率较高,但由于其解调方式决定了其采集速度相对较慢,解调系统设备较复杂,成本较高[14]。另一类是基于波长解调方式的多芯光栅或者集成光栅形状传感技术,2019年,Khan等[15]人将四根长度为11.8 cm的多芯光纤集成到一个导管内,进行形状重构,实验结果表明其最大平均误差和最大绝对误差分别为0.44 mm和1.05 mm。2020年,Omar Al-Ahmad等[16]人提出了一种扭转补偿的方法,将长度为170 mm的多芯光纤嵌入导管内对于静态几何形状进行测量,结果表明其平均误差为0.7 mm,最大空间误差为2.52 mm。基于波长解调方式的多芯光栅阵列传感系统具有体积小巧、信噪比高、采集速度快,解调实时性高等[17]优点。鉴于实际应用领域中存在多种多样的形状传感需求,采用波长解调方式的多芯光栅阵列多芯传感系统可根据实际需求进行光栅数量和间距的灵活配置,可实现具有更灵活的传感距离、更多样的形状变化跨度的待测对象形状监测,具有更广泛的应用前景。然而由于多芯光栅测点在空间上无法实现连续,故无法获取测点之间空白无栅区域的形状,存在测点盲区。在实际应用中,为了提高形状传感的准确性,需要研究待测对象的形状变形特点,在兼顾关键点可测和待测区全范围可测的前提下,设计合理的光栅分布配置方案。

      文中基于七芯光纤光栅的应变传感特性,采用了一种基于齐次变换矩阵[18-19]的三维形状重构方法,对两种不同光栅间距的七芯光纤光栅进行形状重构,通过光栅点波长相对变化值求解有栅测点所在区域的曲率和挠率,通过三次样条插值法求得FBG测点之间空白无栅区域的曲率和挠率,最后将所有点的曲率和挠率整合到同一坐标系下,实现待测对象的三维形状重构。为了探究光栅间距对形状重构精度的影响,利用算法编程模拟了基于该算法原理下不同光栅间距下的三维形状重构误差,通过搭建一套三维形状传感系统完成实验验证,探讨了实验中的误差和模拟误差的合理性,为今后多芯光栅形状传感技术应用中的光栅配置设计提供一定的参考依据。

    • 多芯光栅阵列三维形状重构方法的前提在于有栅测点区域的曲率和挠率的精确监测,继而采用插值与坐标解算实现整体形状重构,解算过程共包括五个阶段,如图1所示。首先通过解调仪获取所有有栅测点未形变前与形变后的中心波长数据,根据每个FBG的波长漂移量进行曲率和挠率的计算,然后对这些FBG点之间的空白无栅区域的曲率和挠率利用三次样条插值法进行插值,再利用齐次矩阵法进行连续坐标解算,并整合至同一坐标系,最后绘制出光纤的整体三维形状。以下将重点对曲率和挠率计算以及齐次矩阵法进行阐述。

      图  1  三维形状重建技术路线图

      Figure 1.  3D shape reconstruction technology flow chart

    • 图2(a)所示,当七芯光纤朝某一个方向弯曲时,边芯(顺时针编号2~7)和中芯(编号1)两纤芯的中心之间的距离为d,一般称为芯间距。当光纤发生弯曲后,假设光纤弯曲的曲率半径为RN为垂直于弯曲方向的中性轴,则任意边芯的轴向应变ε可以表示为:

      $$ \varepsilon =\frac{d}{R}\cdot \mathrm{cos}\theta $$ (1)

      式中:θ为横截面上偏芯和弯曲方向之间的夹角。处于中性面的外侧的边芯为拉伸状态,处于中性面的内侧的边芯为压缩状态。由于中间芯位于中性面上,其中心波长可视为不受光纤弯曲的影响,可以进行轴向应变和温度的补偿。光纤弯曲曲率与波长漂移之间的线性关系可以表示为:

      $$ \Delta {\lambda }_{n}-\Delta {\lambda }_{1}=k\cdot d\cdot \rho \cdot {\lambda }_{n}\cdot \mathrm{cos}{\theta }_{n}\text{,}n=2,3,\cdot \cdot \cdot ,7 $$ (2)

      式中:Δλn为偏芯波长漂移量;Δλ1为中芯光栅中心波长漂移量;ρ=1/R为弯曲的曲率;k=(1−Pε)为标准石英光纤纵向应变敏感系数,通常为0.784;λn为多芯光纤光栅中偏芯在自由平直无弯曲状态下的初始波长。实验过程中,当检测到波长漂移时,利用中芯和任意两个边芯的数据分别代入公式(2)联立起来,就可以推导出每个芯与弯曲方向的夹角以及弯曲的曲率。如图2(b)所示,以边芯2芯为例,可以求出相邻光栅的2芯与弯曲方向的夹角$ {\theta }_{i} $,再根据公式(3)就可以求出各相邻FBG点之间的挠率大小。

      图  2  七芯光纤传感原理示意图。(a)曲率;(b)挠率

      Figure 2.  Schematic diagram of the seven-core optical fiber sensing principle. (a) Curvature; (b) Torsion

      $$\begin{split} \\ {\tau }_{12}={\theta }_{2}-{\theta }_{1}\text{;}{\tau }_{23}={\theta }_{3}-{\theta }_{2}\cdots \end{split} $$ (3)

      根据公式(1)、(2)分析七芯光纤光栅的波长漂移量,可以计算出FBG处的弯曲曲率与弯曲方向,再根据公式(3)可以计算出挠率的大小。显然由于七芯光纤的六个边芯是对称分布的,所有方向上都具有很好的弯曲分辨能力,同时相较于只有两个边芯且互呈90°分布的多芯光纤,七芯光纤多出来的边芯可以用作取平均值,这样可以减小偶然误差对实验结果的影响,提高测量精度。

    • 根据1.1节中计算出的曲率与挠率,利用齐次变换矩阵法(公式(4)~(6))将离散的点坐标解算出来。

      $$ \tilde S_{i + 1}^0 = M_1^0M_2^1 \cdot \cdot \cdot M_i^{i - 1}M_{(i + 1)}^i\tilde S_{i + 1}^{i + 1} = M_{i + 1}^0\tilde S_{i + 1}^{i + 1} $$ (4)
      $$ \begin{split} \tilde S_{i + 1}^i = \left[ {\begin{array}{*{20}{c}} {s_{i + 1,x}^i} \\ {s_{i + 1,y}^i} \\ {s_{i + 1,z}^i} \\ 1 \end{array}} \right] = M_{i + 1}^i\tilde S_{i + 1}^{i + 1} = \left[ {\begin{array}{*{20}{c}} {R_{i + 1}^i}&{{\text{o}}_{i + 1}^i} \\ O&1 \end{array}} \right]\tilde S_{i + 1}^{i + 1} \end{split}$$ (5)
      $$ M_{i + 1}^i = \left[ {\begin{array}{*{20}{c}} {\cos {\alpha _i}\cos {\varphi _i}}&{ - \sin {\alpha _i}}&{\cos {\alpha _i}\sin {\varphi _i}}&{{R_i}(1 - \cos {\varphi _i}){\cos} {\alpha _i}} \\ {\sin {\alpha _i}\cos {\varphi _i}}&{\cos {\alpha _i}}&{\sin {\alpha _i}\sin {\varphi _i}}&{{R_i}(1 - \cos {\varphi _i}){\sin} {\alpha _i}} \\ { - \sin {\varphi _i}}&0&{\cos {\varphi _i}}&{ {R _i}\sin {\varphi _i}} \\ 0&0&0&1 \end{array}} \right] $$ (6)

      式中:$ {\tilde{S}}_{i+1}^{i} $是以FBGi为坐标系来表示的FBGi+1的齐次坐标点;$ {s}_{i+1,k}^{i} $中的k={x, y, z}是$ {\tilde{S}}_{i+1}^{i} $的位置分量;${{M}}_{{i+1}}^{{i}}$是从坐标系i+1到坐标系i的齐次变换矩阵;$ {\tilde{S}}_{i+1}^{i+1} $是FBGi+1在FBGi+1坐标系下的齐次坐标,即$ {{[0 \; 0\; 0\; 1]}}^{\text{T}} $$ {\alpha }_{i} $为弯曲方向角的变化量;$ {\phi }_{i} $为光纤圆弧所对应的圆心角度数。齐次变换矩阵可以看成由四个矩阵组合而成,每个矩阵发挥不同的作用,其中${{R}}_{{i+1}}^{{i}}$是一个3×3的旋转矩阵,其作用为坐标系对齐,${{o}}_{{i+1}}^{{i}}$i+1坐标系的原点在i坐标系下的坐标,${{O}}_{\text{3}}^{{\text T}}$是一个1×3的零向量。因此,在形状重构期间将采用公式(6)这种坐标变换方法来计算以起始坐标系表示的所有FBG光栅点位置。齐次变换矩阵将每个独立坐标系中的点转换到同一坐标系下,将这些点按照顺序连接起来即可实现三维空间曲线的重构,如图3所示。

      图  3  三维形状重构示意图

      Figure 3.  Schematic diagram of 3D shape reconstruction

    • 以半径为10 cm、螺距为9 cm、总长为160 cm的圆柱螺旋线为对象,基于上述三维形状重构算法进行编程计算,探究七芯光纤光栅不同光栅间距对三维螺旋线形状重构精度的影响。选取光栅间距L为12.5 cm、10 cm、8 cm、5 cm和1 cm,绘制出这些光栅间距下的三维重构曲线与真实曲线的图形,如图4(a)所示。该模拟计算均假设有栅测点的曲率和挠率测量误差为0。从图4(b)中可以看,出当光栅间距L为12.5 cm、10 cm、8 cm、5 cm和1 cm,空间位置误差随着光纤长度呈现一种逐渐上升趋势,光栅间距越大,误差也会越大,其重建最大误差均落在模拟长度的尾点,分别为7.75 cm、4.35 cm、2.63 cm、0.94 cm和0.25 cm,占全长的4.8%、2.7%、1.6%、0.6%和0.16%。其原理是由于三维形状的每个坐标点都由前一个点经平移旋转变化而来,所以会由于算法的迭代引起发散误差。为了探究误差变化率与不同长度和间距的总体变化趋势,采用线性拟合方式对这几条误差曲线进行拟合,拟合的斜率分别为0.048、0.027、0.016、0.006和0.001,改斜率即对应于不同光栅间距的三维形状的误差百分率。通过建立间距和线性拟合的斜率之间的变化关系,可探讨不同光栅间距与三维重建误差之间的对应关系。如图4(c)所示,对这六个点进行指数拟合,拟合相关系数R2达0.999。因此,可以根据这一拟合曲线函数对任一光栅间距的三维形状重构的误差进行预测。

      图  4  不同光栅间距的模拟结果。(a)三维重构效果图;(b)空间位置误差;(c)光栅间距与误差曲线斜率之间的关系

      Figure 4.  Simulation results of different grating spacings. (a) Diagram of 3D reconstruction effect; (b) Spatial position error; (c) Relationship between grating spacing and error curve slope

    • 为了验证上述仿真结果的可靠性,设计了一套实验装置。图5(a)为实验系统图,该系统主要由七芯光纤光栅、七芯光纤扇入扇出装置、光纤光栅解调仪和上位机构成。实验所选用的七芯光纤是长飞光纤光缆股份有限公司生产的,由一个中芯和六个边芯组成,芯间距为42 μm,如图5(b)所示。通过掩模版法,在七芯光纤上写入不同中心波长的光栅。光纤光栅解调仪是由武汉高斯光电科技有限公司生产的,其分辨率为10 pm。光纤总长80 cm,选取了两种不同光栅间距的光栅串,分别为10 cm和5 cm,栅区的长度为1 cm,各个纤芯的光栅光谱以及光栅串的光谱如图6(a)、(b)所示。为了避免光纤的自然扭转效应引起的误差,将光纤做成并带的形式,即将多根光纤用并带胶粘在一起形成带状结构,如图5(c)所示,这样就可以避免光纤在测量过程中发生扭转。然后将并带结构放入二维曲率槽内,再用公式(3)进行预先的扭转角标定。标定完后,需要将七芯光纤光栅笔直地放置实验平台上,在无任何应力应变的干扰下测量得到光栅的初始中心波长,初始波长是否准确将会影响后续算法的计算结果,所以必须十分准确地测量才能获得高精度的重构形状。随后将光纤放置在3D打印的圆柱三维螺线凹槽内,如图5(d)所示。圆柱的半径为100 mm,螺旋槽的螺距为90 mm。

      图  5  三维形状传感实验装置。(a)系统简图;(b)七芯光纤端面;(c)光纤并带结构;(d)三维螺旋曲率槽

      Figure 5.  Experimental setup for 3D shape sensing. (a) Schematic diagram of the system; (b) End face of seven-core fiber; (c) Optical fiber with ribbon structure; (d) 3D helical curvature groove

      图  6  (a)各个纤芯的光栅光谱;(b) 光栅串的光谱

      Figure 6.  (a) Grating spectrum of each fiber core; (b) Spectrum of grating string

    • 根据波长漂移量计算出各光栅点的曲率和挠率,实验先后进行了五次测量,取其平均值,利用三次样条插值法对空白处的曲率和挠率进行插值,再结合齐次变换矩阵重建算法得到了三维重建效果图,如图7(a)所示。其中红色曲线为光栅间距为10 cm的实验重建三维螺旋线,蓝色曲线为光栅间距为5 cm的实验重建三维螺旋线,黑色曲线为真实三维螺旋线。实验结果表明,重构出的三维形状与真实形状符合得较好。图7(b)是对重建图形的误差分析,从图中可以看出,整体误差随着光纤长度的增加呈现出一种增长的趋势,最大误差出现在尾点处。表1为实验值与模拟值的数据对比,从表中可以看出,光栅间距5 cm时的最大误差为1.15 cm,占全长的1.4%,平均误差为0.62 cm,占全长的0.8%。光栅间距10 cm时的最大误差为2.56 cm,占全长的3.2%,平均误差为1.32 cm,占全长的1.7%。第3节得出模拟值的最大误差占全长的2.7%和0.6%,实验值与模拟值符合得较好,两组实验值与模拟值的相对误差分别为0.5%和0.8%。误差产生的主要原因有:(1)通过显微镜观察发现,实验所用的七芯光纤的六个边芯并非完全按正六边形排布的,而是有一定的误差,这对计算弯曲曲率会产生一定的影响,从而导致实验误差;(2)实验中采用的是分辨率为10 pm的解调仪解调FBG的峰值波长,当波长变化小于10 pm时无法检测到,存在系统误差。

      图  7  三维形状传感实验。(a)重构曲线与真实曲线;(b)误差

      Figure 7.  3D shape sensing experiment. (a) Reconstructed curve and true curve; (b) Error

      表 1  实验值与模拟值对比

      Table 1.  Comparison of experimental and simulated values

      Grating pitch,
      L/cm
      Maximum error
      /cm
      Error percentageSimulation error percentage
      51.151.4%0.6%
      102.563.2%2.7%
    • 文中探究了不同光栅间距对七芯光纤光栅三维形状传感精度的影响关系,通过模拟计算五组不同光栅间距下的三维螺旋线形状的重构误差,建立了误差与光栅间距之间的关系,然后在选取两组光栅间距进行实验,并将光纤做成并带的形式以防止光纤扭转带来的误差。结果表明,最大误差均出现在尾点处,分别为2.56 cm和1.15 cm,占全长的3.2%和1.4%,平均误差为1.32 cm和0.62 cm,占全长的1.7%和0.8%。无论是模拟仿真还是实验验证,误差值与待测对象相对长度的整体趋势符合线形增长的趋势。最大误差点均发生于尾点处。通过建立光栅间距和线性拟合的斜率之间的变化关系,可探讨不同光栅间距与三维重建误差之间的对应关系。根据该关系可以对相似的形状下的任一光栅间距、任意长度光纤的三维形状重构误差进行预测,这样就可以结合具体的应用场景选择合适的光栅间距以及解调方法,合理配置测点资源,在较低的成本范围内提高检测性能。此外,提高光纤光栅解调仪的分辨率将有利于获得更精确的波长数据,得出更精确的曲率值;该三维重构算法原理是三维形状的每个坐标点都由前一个点经平移旋转变化而来,所以会由于算法的迭代引起发散误差,可以通过增加等距校正点、可收敛发散误差有效抑制尾点处误差。该方法有望为不同距离待测对象的光栅布置方案与栅点间距参数提供一定的理论依据与设计参考。

参考文献 (19)

目录

    /

    返回文章
    返回