Simultaneous measurement sensor of temperature and stress based on UV lithography phase-shifted grating
-
摘要: 传统光纤光栅传感器存在温度应力交叉敏感的问题,无法同时测量被测物体的温度应力变化情况。针对这种情况,提出一种新型的利用紫外光刻写相移光栅的方法,在刻写光栅前用电极放电的方法去除掉极小一段光纤的光敏性,使光纤原有的均匀周期分布状态被破坏,从而形成相移光栅,并对其进行理论分析。此种相移光栅的透射光谱存在2个明显的谐振峰,利用2个峰对温度和应变灵敏度不同的性质,可以通过建立解调矩阵来实现温度与应力的同时测量。实验结果表明:利用此种方法制成的相移布拉格光栅能够较为精确地实现温度、应力的同时测量,所制作的传感器温度灵敏度最高可达9.51 pm/℃,灵敏度方差低于2.125×10−7,应变灵敏度最高可达0.767 pm/με,灵敏度方差低于2.156×10−10。Abstract: The traditional fiber grating sensors are cross-sensitive to temperature stress and cannot simultaneously measure the change of the temperature stress of the measured object. In view of this situation, a new method of lithographing phase-shifted gratings using ultraviolet light was proposed. Before lithographing the gratings, the method of electrode discharge was used to remove the photosensitivity of a very small section of the optical fiber, so that the original uniform periodic distribution of the fiber was destroyed to form a phase-shifted grating, and its theoretical analysis was carried out. There were two distinct resonance peaks in the transmission spectrum of this phase-shifted grating. Using the properties of two peaks with different sensitivity to temperature and dependent variable, the simultaneous measurement of temperature and stress could be achieved by establishing the demodulation matrix. The experimental results show that the phase-shifted Bragg grating made by this method can realize the simultaneous measurement of temperature and stress accurately. The maximum temperature sensitivity of the sensor can reach to 9.51 pm/℃, and the sensitivity variance is lower than 2.125×10−7. The maximum strain sensitivity can reach to 0.767 pm/με, and the sensitivity variance is lower than 2.156×10−10.
-
Keywords:
- fiber grating /
- simultaneous measurement /
- UV lithography /
- temperature /
- stress
-
引言
安防镜头的测试标准为在−40 ℃~80 ℃的高低温环境下保证成像稳定。考虑到镜头成本的因素,企业大多采用玻塑混合的设计方案,但是在高低温环境下,塑胶镜片的热变形远大于玻璃镜片的热变形,其面型、厚度等更容易受机械结构的影响而产生变化,直接导致光学系统性能下降。为保证安防镜头在高低温环境下成像稳定,在光学设计阶段会利用塑胶镜片的特性对光学系统进行光学消热设计,以获得较好的高低温稳定性。但是,消热设计所选塑胶镜片材料的热膨胀系数往往大于镜框的热膨胀系数,高温工况下塑胶镜片与镜框会相互挤压, Zemax中消热设计并不能模拟镜面受挤压的应变变形等,因此,需要对镜头进行光机热一体化分析。
光机热一体化分析广泛应用于光机结构设计过程中,通过设计-仿真-设计反复迭代,优化光机结构对高低温环境的适应能力[1-5]。BROOME B G等[6]设计的机械被动式无热化镜头釆用了不同热膨胀系数材料搭配,并配合不同接触面形状的方式来实现光学系统的无热化。长春理工大学与长春光学精密机械研究所的王平、张国玉等[7]在航空变焦镜头机械被动式无热化设计中采用了差动元件,保证了在高低温环境下光学系统中透镜间的间隔。东北电子技术研究所的陈德富、李相军等[8]采用机械被动与光学组合的无热化方式,实现了对镜头离焦量的补偿。上述研究虽然对光学系统在热力载荷下的成像质量进行了诸多探讨,但针对热力载荷下机械结构对光学系统消热设计的影响,以及使机械消热设计与光学消热设计相互抵消的研究较少[9-10]。
本文以某款安防定焦镜头为分析对象,使用有限元分析计算光机结构在高低温环境下的热变形,选择合适材料的底座来抵消镜头的热离焦量,分析镜框材料分别为聚碳酸酯(polycarbonate,PC)+30%玻璃纤维(fiberglass,FG)和PC+20% FG时,镜片与镜框相互作用产生的应变。通过Zernike多项式拟合以及光机热一体化仿真,分析光学系统调制传递函数(modulation transfer function,MTF)变化,明确了PC+20% FG材料的底座和PC+20% FG材料的镜框对镜头热补偿的有效性。利用光机热一体化分析技术有效提高了光机系统的高低温稳定性,为玻塑混合定焦光学系统的优化设计提供理论与技术参考。
1 光学与机械热补偿公式
在光学系统处于不同温度时,透镜间距随着镜筒热胀冷缩而改变,进而影响整个系统的成像质量。因此,以三分离薄透镜为基础建立了数学模型,3组薄透镜光焦度分别为 Φ1、Φ2、Φ3,间距分别为 d1、d2,后截距(BFL)为 Lf,相应的镜筒材料热膨胀系数分别为 αH1、αH2、αH3。第一近轴光线在3组薄透镜上的入射高度分别为 h1、h2、h3。
三分离薄透镜光学系统总光焦度为
$$ \varPhi ={{{\displaystyle\sum\limits_3^{i=1}}}}{h}_{i}{\varPhi }_{i} $$ (1) 后截距Lf为
$$ {L_{{f}}} = \frac{{1 - {\varPhi _1}\left( {{d_1} + {d_2}} \right) - {\varPhi _2}{d_2}\left( {1 - {\varPhi _1}{d_1}} \right)}}{\varPhi } $$ (2) 温度变化Δt后,各薄透镜的光焦度为
$$ {\varPhi '_i} = {\varPhi _i}\left( {1 + {T_i}\Delta t} \right) $$ (3) 式中Ti为各透镜光学热特性参数。
由于各透镜之间由隔圈与镜框连接,因此当温度变化时,透镜间隔因材料的热胀冷缩而改变,变化后的透镜间隔为
$$ {d'_i} = {d_i}\left( {1 + {\alpha _H}_i{{\Delta }}t} \right) $$ (4) 温度变化后的系统总光焦度 Φ′为
$$ \begin{split} \varPhi ' = &{\varPhi '_1} + {\varPhi '_2} + {\varPhi '_3} - {\varPhi '_1}{d'_1}\left( {{{\varPhi '}_2} + {{\varPhi '}_3}} \right) -\\ & {\varPhi '_3}{d'_2}\left( {{{\varPhi '}_1} + {{\varPhi '}_2} - {{d'}_1}{{\varPhi '}_1}{{\varPhi '}_2}} \right) \end{split}$$ (5) 温度变化后的光学系统后截距Lʹ f为
$$ {L'_{{{f}}}} = \frac{{1 - {{\varPhi '}_1}\left( {{{d'}_1} + {{d'}_2}} \right) - {{\varPhi '}_2}{{d'}_2}\left( {1 - {{\varPhi '}_1}{{d'}_1}} \right)}}{{\varPhi '}} $$ (6) 当后截距部分相应底座的热变化量与光学系统焦距的变化量恰好相等时(差值Δz = 0),系统实现无热化,即:
$$ \begin{split} {{\Delta }}{\textit{z}} =& {L_{{f}}}\left( {1 + {\alpha _{Hi}}{{\Delta }}t} \right) - {{L'}_{{{f}}}} = \\ &\frac{{\left[ {1 - {\varPhi _1}{d_1} - {d_2}\left( {{\varPhi _1} + {\varPhi _2} - {\varPhi _1}{\varPhi _2}{d_1}} \right)} \right] \cdot \left( {1 + {\alpha _{Hi}}{{\Delta }}t} \right)}}{\varPhi } - \\ & \frac{{1 - {{\varPhi '}_1}\left( {{{d'}_1} + {{d'}_2}} \right) - {{\varPhi '}_2}{{d'}_2}\left[ {1 - {{{{\varPhi '}}}_1}{{d'}_1}} \right]}}{{\varPhi '}} = 0 \\[-15pt] \end{split} $$ (7) 综上所述,定焦镜头消热过程就是使各镜片光焦度变化和机械结构变形相互补偿的过程,而镜片光焦度的变化量又可以通过选择不同热特性参数的透镜组合加以改变,使之与机械热差匹配。因此,消热的关键在于选取合适的光学材料和结构材料组合[10-13]。
2 镜头光学系统设计
根据设计指标要求,镜头总长为22.205 mm,光学系统后焦距为5.498 mm,F数1.68,波长0.436 µm~0.85 µm,像质要求MTF大于0.3,测试温度范围为−40 ℃~80 ℃。光学系统如图1所示。镜头采用1枚玻璃镜片搭配4枚塑胶镜片的形式,第1、2、4、5枚镜片为塑胶镜片,第3枚为玻璃镜片。
光学系统在初始设计时进行消热设计,镜框、隔圈材料的热膨胀系数根据经验确定,光学系统MTF如图2所示。从图2可以看出,镜头在常温下各视场MTF大于0.3,可以满足设计指标要求。高温和低温组态下光学系统中心视场离焦量分别为13.5×10−3 mm、13.0×10−3 mm,需要搭配合适材料的底座以抵消镜头本身的离焦量。
3 镜头结构设计
定焦镜头主要由镜框、镜片、隔圈、底座等组成,镜片和隔圈的设计需要考虑其承靠稳定性与加工成本,此镜头采用平面承靠的方式固定镜片。另外,隔圈在此镜头中还兼顾着光阑的作用,因此隔圈还设计了充当光圈的凸台。为补偿Zemax中模拟光学系统本身的热离焦量,底座材料初步选择为PC+20% FG,结构设计如图3所示。
4 定焦镜头光机热一体化仿真分析
4.1 热结构有限元分析
安防镜头的镜框、镜片和隔圈全部为回转体,对镜头模型不产生影响的特征进行简化,利用Hypermesh对三维模型进行手动六面体网格划分,保证曲面上网格节点精度,接触面最少划分3层网格,共划分38 566个高质量六面体网格,网格雅可比大于0.7的占有率为92.7%,有限元模型网格剖视图如图4所示。
根据镜头测试要求,在Ansys workbench中对镜头进行热结构耦合分析,镜头在测试环境所受载荷为稳态热,温度载荷设置为极限温度−40 ℃和80 ℃,模型参考温度为常温25 ℃,边界条件根据实际工况设置为固定底座的导柱。镜头各部件材料参数如表1所示,镜头位移和应变云图如图5所示。
表 1 镜头部件材料参数Table 1. Material parameters of lens parts材料 弹性模量/
MPa密度/10−6
kg·mm−3泊松比 热膨胀系数
/10−5 ℃径向/
周向轴向 PC+30%FG 9510 1.43 0.36 4.21 3.17 PC+20%FG 5910 1.43 0.37 4.82 3.57 Plastic lens-1 2410 1.01 0.38 4.56 Plastic lens-2 2710 1.25 0.35 7.01 Glass lens-1 72000 2.53 0.32 1.02 Plastic lens-3 2430 1.12 0.33 6.79 Plastic lens-4 2510 1.03 0.32 6.54 Brass 11000 8.5 0.31 1.55 2320UV-7A2 1680 1.72 0.38 2.52 由计算结果可知,极限测试温度载荷下,镜头在底座热变形带动下,相对于相机参考面分别移动了13.3×10−3 mm、13.7×10−3 mm,可以有效补偿Zemax中模拟光学系统本身的热离焦量,PC+20% FG材料底座满足补偿要求。
高温80 ℃载荷下,由于塑胶镜片的热膨胀系数大于镜框径向的热膨胀系数,此时镜片膨胀受到镜框挤压,第2枚塑胶镜片上最大等效应变为2.36×10−3 mm。
低温−40 ℃载荷下,塑胶镜片和玻璃镜片受机械结构作用力产生的应变很小,这是因为塑胶镜片的热膨胀系数大于镜框的径向热膨胀系数,不会产生挤压作用力;玻璃镜片的弹性模量大于镜框的弹性模量,而且隔圈的热膨胀系数小于玻璃镜片的热膨胀系数,镜框低温收缩产生的力主要由隔圈承受,从而在一定程度上保护了玻璃镜片,所以玻璃镜片的挤压应变也非常小。
4.2 基于Sigfit的面型处理
连续的光学表面都可以用Zernike多项式线性组合表示,其具有在单位圆上相互正交、旋转对称性、与初级像差存在着一定对应关系等特点,且光学设计软件如Zemax、Code V等均支持 Zernike面形,所以将其常作为有限元分析与光学分析之间理想的接口工具[14-17]。本文采用Zernike多项式拟合变形后的面型。
Zernike多项式极坐标数学描述为
$$ \begin{split} {Z_i}\left( {r,\theta } \right) &= Z_n^m\left( {r,\theta } \right) \\ &= \left\{ {\begin{array}{*{20}{c}} {R_n^m\left( r \right){\text{cos}}\left( {m\theta } \right)}&{\left( {m \geqslant 0} \right)} \\ {R_n^m\left( r \right){\text{sin}}\left( {m\theta } \right)}&{(m < 0)} \end{array}} \right. \end{split} $$ (8) 式中:r为归一化半径,0 ≤ r ≤ 1;θ为极角,且 0 ≤ θ ≤ 2π;$R_{{n}}^{{m}}\left( r \right)$为径向多项式,可以表示为
$$ R_n^m\left( r \right) = \mathop {\mathop \sum \limits^{\tfrac{{n - m}}{2}} }\limits_{s = 0} {( - 1)^s}\dfrac{{\left( {n - s} \right)!}}{{s!\left( {\dfrac{{n + m}}{2} - s} \right)!\left( {\dfrac{{n - m}}{2} - s} \right)!}}{r^{\left( {n - 2s} \right)}} $$ (9) 式中:n为多项式的阶数,n=0,1,…;m 为序号,其值恒与 n 同奇偶性,且 m ≤ n。
在Sigfit中输入镜片基本参数,包括波长、曲率半径、非球面系数和包含每个透镜前后表面节点变形的.asig文件。高温80 ℃载荷下,分离出的镜片刚体位移和旋转量如表2所示。将分离刚体位移后的面型变化量使用Zernike多项式拟合,将结果输入到Zemax中,获得安防镜头施加热载荷后光学系统MTF曲线,如图6所示。
表 2 定焦镜头镜面刚体位移和旋转量Table 2. Mirror rigid body displacement and rotation of fixed-focus lensmm 镜面序号 X 刚体位移 Y 刚体位移 Z 刚体位移 X旋转刚体位移 Y旋转刚体位移 Z旋转刚体位移 1 9.95×10−4 −6.24×10−4 −5.05×10−2 −2.04×10−6 −3.34×10−6 8.83×10−5 2 9.86×10−4 −6.21×10−4 −5.08×10−2 −1.45×10−6 −4.80×10−7 8.83×10−5 3 9.86×10−4 −6.22×10−4 −3.66×10−2 9.28×10−7 −2.68×10−6 8.84×10−5 4 9.81×10−4 −6.24×10−4 −2.47×10−2 1.31×10−6 −2.71×10−6 8.84×10−5 5 9.89×10−4 −6.37×10−4 −2.85×10−2 −1.15×10−6 −2.66×10−6 8.77×10−5 6 9.81×10−4 −6.35×10−4 −2.54×10−2 −6.19×10−7 −2.72×10−6 8.77×10−5 7 9.83×10−4 −6.15×10−4 −2.80×10−2 −1.42×10−6 −2.13×10−6 8.78×10−5 8 9.79×10−4 −6.09×10−4 −1.63×10−2 −3.23×10−6 −2.26×10−6 8.79×10−5 9 9.75×10−4 −6.04×10−4 −1.64×10−2 −1.77×10−6 −1.91×10−6 8.77×10−5 10 9.73×10−4 −6.02×10−4 −1.38×10−2 −3.00×10−6 −3.14×10−6 8.78×10−5 光机热一体化分析结果表明,镜片受到镜框的挤压力过大,导致镜片光焦度变化量与光学消热设计要求的变化量产生差异,MTF曲线下降过大,已无法满足成像要求,需要改变光机结构的补偿措施。
4.3 PC+20% FG镜框有限元分析
经光机热一体化分析可知,镜框对塑胶镜片挤压应变过大,高低温环境下镜头已经不能满足使用要求。重新选择镜框材料为PC+20% FG,而且PC+20% FG的材料弹性模量也有所降低,在一定程度上减小了对内部镜片的挤压力。有限元仿真计算结果如图7所示。
−40 ℃载荷下,第1枚镜片产生微小应变,这是因为塑胶镜片的热膨胀系数大于镜框的轴向热膨胀系数,第1枚镜片受第2枚镜片的推力和镜框收缩的拉力影响形成力矩,所以产生微小应变。第2枚镜片的最大等效应变为0.53×10−3 mm,光机热一体化分析此时光学系统MTF曲线如图8所示。
经光机热一体化分析可知,此光机系统搭配PC+20% FG材料的镜框和PC+20% FG材料的底座满足光学补偿条件。高低温载荷下安防镜头光学系统MTF曲线有一定程度下降,但是总体大于0.3,符合设计要求。
5 镜头法兰焦距测量试验
为了验证设计和仿真分析的有效性,利用ImageMaster® HR紧凑型高精度MTF测量仪(型号:ImageMaster HR TempControl VIS)进行了法兰焦距测量实验。将镜头放置在对应台具上,将镜头所处空间箱体抽真空,防止低温空气凝结水珠,影响光线传递。通过辐射传热的方式,给镜头加热或降温,由测量设备内置相机检测镜头焦点位置,测出不同温度条件下镜头的焦点位置的变化量。为防止设备本身热胀冷缩造成的误差,设置一固定参照物,测量相机在不同温度下与固定参照物的距离,消除设备本身热变形的误差。测量平台与试验示意图如图9所示。
测量结果如表3所示。从表3可知,镜头法兰焦距在高低温环境下的变化量在5 µm以内,虽然镜头法兰焦距稍有变化,但仍然满足成像要求,可间接反映设计和仿真结果的有效性。
表 3 镜头高低温法兰焦距测量结果Table 3. Measurement results of high and low temperature flange focal length of lens温度/℃ 焦点位置
/mm参照物参考
位置/mm参照物
距离 /mm法兰焦距
/mm19.997 35.08617 25 9.645075 −0.44110 −39.994 35.09468 25 9.652904 −0.44178 20.001 35.09342 25 9.652887 −0.44053 79.998 35.07347 25 9.626284 −0.44518 20.003 35.08694 25 9.646567 −0.44038 6 结论
对某玻塑混合安防定焦镜头进行了光机热一体化分析,设计并验证了其高低温补偿措施。有限元分析结果表明,PC+20% FG材料的底座可以有效补偿Zemax中模拟光学系统本身的热离焦量。镜框材料为PC+30% FG时,镜框径向热膨胀系数与镜片热膨胀系数相差过大,在高温环境下镜片受到镜框内壁挤压,产生不均匀变形,产生的等效应变为2.36×10−3 mm,镜面光焦度无法按照光学消热设计的要求变化,造成镜头像质下降超出使用范围。当镜框材料为PC+20% FG时,镜框径向热膨胀系数增大,弹性模量减小,此时镜框对塑胶镜片产生的等效应变减小为0.53×10−3 mm。光机热一体化分析结果表明, 镜框材料为PC+20% FG时,光学系统MTF略有下降,但总体MTF值大于0.3,像质基本稳定,满足使用要求。最后通过法兰焦距测量试验,验证了安防镜头温度适应性良好,结构设计合理、有效,光机热一体化分析结果与实验结果基本吻合,为安防镜头研发提供了一种高效、可靠的验证手段。
-
表 1 相同应力条件下相移光栅温度测量灵敏度
Table 1 Temperature measurement sensitivity of phase-shifted grating under same stress
温度/℃ 30 45 55 65 75 85 95 灵敏度(dip A) 0.000 767 0.000 769 0.000 768 0.000 755 0.000 748 0.000 747 0.000 734 灵敏度(dip B) 0.000 767 0.000 774 0.000 768 0.000 754 0.000 748 0.000 745 0.000 734 表 2 相同温度条件下相移光栅应力测量灵敏度
Table 2 Stress measurement sensitivity of phase-shifted grating under same temperature
应力/με 0 400 800 1 200 1 600 2 000 灵敏度(dip A) 0.009 534 0.009 177 0.009 087 0.008 683 0.008 543 0.008 454 灵敏度(dip B) 0.009 62 0.009 086 0.009 07 0.008 773 0.008 446 0.008 392 -
[1] 杨思玉, 万生鹏, 王浩宇, 等. 单点相移光纤光栅光谱特性的研究与应用[J]. 应用光学,2019,40(2):349-355. YANG Siyu, WAN Shengpeng, WANG Haoyu, et al. Research and application of spectrum characteristics of single phase-shifted fiber Bragg grating[J]. Journal of Applied Optics,2019,40(2):349-355.
[2] 张宝荣, 黄震, 张玉唯. 基于叠栅的Bragg光纤光栅综合方法的研究[J]. 应用光学,2019,40(4):699-703. ZHANG Baorong, HUANG Zhen, ZHANG Yuwei. Study of Bragg grating synthesis methods based on superimposed fiber grating[J]. Journal of Applied Optics,2019,40(4):699-703.
[3] 童杏林, 何为, 张翠, 等. 光纤光栅与光纤法珀传感器在航空航天领域的研究与应用进展[J]. 激光杂志,2018,39(7):1-7. TONG Xinglin, HE Wei, ZHANG Cui, et al. Research and application progress of fiber Bragg grating and Fabry-Perot sensors in the field of aeronautics and astronautics[J]. Laser Journal,2018,39(7):1-7.
[4] 徐国权, 熊代余. 光纤光栅传感技术在工程中的应用[J]. 中国光学,2013,6(3):306-317. XU Guoquan, XIONG Daiyu. Applications of fiber Bragg grating sensing technology in engineering[J]. Chinese Optics,2013,6(3):306-317.
[5] 李婧怡, 朱振华. 光纤光栅传感技术在航空航天领域中的应用与发展[J]. 中国新通信,2018,20(4):67-69. doi: 10.3969/j.issn.1673-4866.2018.04.057 LI Jingyi, ZHU Zhenhua. Application and development of fiber grating sensing technology in aerospace[J]. China New Telecommunications,2018,20(4):67-69. doi: 10.3969/j.issn.1673-4866.2018.04.057
[6] 王维俊. 基于FBG宏应变空间分布的梁桥监测研究[D]. 苏州: 苏州科技大学, 2019: 1-7. WANG Weijun. Study on the monitoring of beam bridge based on FBG macro strain spatial distribution[D]. Suzhou: Suzhou University of Science and Technology, 2019: 1-7.
[7] 王正方. 桥隧工程安全监测的光纤光栅传感理论及关键技术研究[D]. 济南: 山东大学, 2014: 1-9. WANG Zhengfang. Study on fiber Bragg grating sensing theory and key technology for bridges and tunnels engineering safety monitoring[D]. Jinan: Shandong University, 2014: 1-9.
[8] FOOTE P D. Fiber Bragg grating strain sensors for aerospace smart structure[J]. SPIE,1994,2361:162-166.
[9] WEIS R S, BEADLE B M. MWD telemetry system for coiled-tubing drilling using optical fiber grating modulators downhole[C]//12th International Conference on Optical Fiber Sensors, Williamsburg, Virginia. Washington, D. C.: OSA, 1997: 416-419.
[10] 姜德生, 何伟. 光纤光栅传感器的应用概况[J]. 光电子·激光,2002,13(4):420-430. doi: 10.3321/j.issn:1005-0086.2002.04.026 JIANG Desheng, HE Wei. Review of applications for fiber Bragg grating sensors[J]. Journal of Optoelectronics Laser,2002,13(4):420-430. doi: 10.3321/j.issn:1005-0086.2002.04.026
[11] RUPALI S, SWEE C T. Thermal compensation of Bragg grating for structural health monitoring[J]. SPIE,1998,4934:150-157.
[12] 赵洪霞, 程培红, 丁志群, 等. 双锥形光纤光栅实现温度、折射率和液位同时测量[J]. 中国激光,2016,43(10):231-237. ZHAO Hongxia, CHENG Peihong, DING Zhiqun, et al. Simultaneous measurement of temperature, refractive index and liquid level based on biconical fiber gratings[J]. Chinese Journal of Lasers,2016,43(10):231-237.
[13] CAVALEIRO P M, ARAUJO F M, FERREIRA L A, et al. Simultaneous measurement of strain and temperature using Bragg gratings written in germanosilicate and boron-codoped germanosilicate fibers[J]. IEEE Photonics Technology Letters,1999,11(12):1635-1637.
[14] JIANG Y J, YUAN, XU J, et al. Phase-shifted fiber Bragg grating inscription by fusion splicing technique and femtosecond laser[J]. Optics and Lasers in Engineering,2016,86:236-241.
[15] JIANG Y J, LIU C, LI D, et al. Simultaneous measurement of temperature and strain using a phase-shifted fiber Bragg grating inscribed by femtosecond laser[J]. Measurement Science and Technology,2018,29(4):045101.
[16] YANG W, GENG T, YANG J, et al. A phase-shifted long period fiber grating based on filament heating method for simultaneous measurement of strain and temperature[J]. Journal of Optics,2015,17(7):075801.