CN112484652A - 一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法 - Google Patents
一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法 Download PDFInfo
- Publication number
- CN112484652A CN112484652A CN201910866719.6A CN201910866719A CN112484652A CN 112484652 A CN112484652 A CN 112484652A CN 201910866719 A CN201910866719 A CN 201910866719A CN 112484652 A CN112484652 A CN 112484652A
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- strain
- birefringent
- birefringence
- initial
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/161—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供一种一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法,通过任意一组相邻的波峰波谷波长λ'n、λn及其双折射光纤初始长度、初始双折射率和初始双折射应变系数可计算双折射光纤所受应变大小。根据干涉光谱波谷波长、波峰波长的相对位置蕴含着应变信息的特点,与初始相角无关,可以剔除外界干扰,提高测量精度。该方法无需人为判断,有助于实现计算机在线测量。且该方法需要的信息量较小,可以实现短Bi‑FLM传感器的在线测量。本申请研究结果对Bi‑FLM各类传感器实现计算机在线测量,提高测量精度具有指导意义。
Description
技术领域
本发明涉及一种通过一组相邻的波峰波谷波长即可实现双折射光纤环镜应变在线测量方法。
背景技术
双折射光纤环镜(Birefringence fiber loop mirror,Bi-FLM)传感器具有低成本,易于制造,极化独立,不受电磁干扰等优点,已成功应用于应变、振动、扭矩等各类传感器。目前,该传感器多通过波长解调的方法实现离线测量,即根据Bi-FLM干涉光谱波长的相对变化量来推算传感量的大小。由于干涉光谱是周期性信号,需要人为判断外界传感量的变化是导致干涉光谱左移还是右移,需要人为判断是较小外界传感量产生较小相角变化导致的干涉光谱平移,还是较大外界传感量产生较大相角的周而复始的干涉光谱,不利于实施传感器的在线测量。且在测试过程中外界干扰容易改变干涉光谱的初始相角,导致干涉光谱平移,从而改变波长的相对变化量,基于波长解调的方法无法区分是干扰还是外界传感量的变化导致的波长相对变化,致使测量精度下降。也有学者基于强度解调原理实现Bi-FLM传感器在线测量,即将Bi-FLM传感器的光信号强度通过光电转换器转换为电信号,通过监测电信号的变化反推光信号的变化,从而反推外界传感量的变化。由于强度解调受光源稳定性影响较大,因此该方法精度较低。
申请人在之前的专利申请(CN 201810656015)中公开了通过任意4个相邻波谷波长及双折射光纤初始条件计算应变大小的方法,该方法无需人为判断,有助于促进传感器与计算机有效对接,实现在线测量。且与初始相角无关,因此可以剔除外界干扰,提高测量精度。但该方法需要4个相邻波谷的波长,包含3个周期的干涉波形,需要的信息量较多。双折射光纤长度越短,则Bi-FLM干涉光谱周期越大,在一定的光源波长范围内,光纤长度较短的Bi-FLM可能没有足够的波谷点(4个),导致无法计算。
发明内容
本发明的目的是提供一种双折射光纤环镜应变传感器在线测量方法,通过任意一组相邻的波峰波谷波长及其双折射光纤初始长度、初始双折射率和初始双折射应变系数即可计算得到双折射光纤所受应变大小。
本发明的技术方案为:一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法,其步骤包括:
1)将光纤耦合器、双折射光纤、光隔离器以及光谱仪构建成双折射光纤环境轴向应变传感器测量系统,其中,双折射光纤的两端分别与光纤耦合器的两个输出臂相连接,入射光经光隔离器与光纤耦合器的输入臂相连接,光谱仪的输入端与光纤耦合器的干涉光谱输出端口相连接;
2)将双折射光纤粘贴在被测物体上,记下双折射光纤初始长度L0、双折射光纤初始双折射率B0和双折射应变系数k;
3)待被测物体承受轴向应变后,此时双折射光纤所受应变即为被测物体所受应变,通过光谱仪测量双折射光纤环镜干涉光谱;
4)通过计算机编程找到任意一组相邻的波峰波谷波长λ'n、λn;
5)将光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k代入公式(24)
计算得到双折射光纤受轴向应变后的绝对长度L';
本发明理论分析如下:
Bi-FLM传感器初始干涉光谱表达式为:
式中,λ为干涉光谱波长,T(λ)为干涉光谱强度,相角θ=2πL0B0/λ,L0为光纤初始长度, B0为光纤初始双折射率。
当双折射光纤受轴向应变后,相角变化量Δθ为:
式中,εz=ΔL/L0=(L'-L0)/L0为双折射光纤轴向应变,单位为ε,其中L'为双折射光纤受应变后的长度。k是双折射应变系数,单位为1/ε,即光纤受1ε后双折射率变化大小。
由(2)式可得通过L'表示的Δθ表达式为:
由(1)、(2)式可得通过εz表示的受轴向应变后的Bi-FLM干涉光谱表达式为:
(4)式是通过双折射光纤所受应变εz描述干涉光谱T'(λ),用于与本文下面推导的计算应变的理论表达式对比,以校验本文推导的表达式正确性。
由(1)式、(3)式可得通过L'表示的受轴向应变后的干涉光谱表达式为:
欲使(5)式对应干涉光谱T'(λ)的值最小,则:
式中,n为整数,λn为整数n对应的波谷波长,其它依此类推。由(6)式解出:
欲使(5)式对应干涉光谱T'(λ)的值最大,则:
式中,n为整数,λ'n为整数n对应的波峰波长,其它依此类推。
由(8)式解出:
由(7)式可得:
由(9)式可得:
由(10)式=(11)式可得:
由(12)代入(10)式可得:
根据(7)式、(9)式可得:
当(14)、(15)式n值相同时,波谷与波峰相邻。由(14)式可知,当双折射光纤所受应变εz为正时,因为波长必须为正数,故n为正,此时λn>λ'n。即εz>0时使n为正,满足λn>λ'n,应用(13)式。
当双折射光纤所受应变εz为负应变且足够使n为负数时,此时λ'n>λn。由(14)式可知,因为波长必须为正数,要使n为负数,则:
B0+(B0+k)εz<0(16)
由(16)式可得:
假设双折射光纤长度L0=0.1m,双折射率B0=2.6×10-4,双折射应变系数k=7.3×0.001/ε,经计算需满足εz<-0.0343915343915344ε,即εz<34391.5343915344με时,才能使n为负数,该应变是个极大的负应变,一般情况下难于产生如此大的负应变,故一般情况下n为正,此时λn>λ'n。下面就n为正时展开讨论,即当满足λn>λ'n时可应用(13)式。
由(13)式可知,双折射光纤受轴向应变后的绝对长度L'可以由任意一组相邻的波峰波长λ'n、波谷波长λn,光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k求出,再将L'代入εz=ΔL/L0=(L'-L0)/L0计算双折射光纤所受应变大小。该方法可以由任意一组相邻的波峰波长、波谷波长及光纤初始条件求出双折射光纤所受应变大小。由(13)式可知,双折射光纤受轴向应变后的绝对长度L'与波峰波长λ'n、波谷波长λn有关,即该方法与波谷波长、波峰波长的绝对位置有关,绝对位置与初始相角有关,故该方法无法区分是干扰还是外界传感量的变化导致的波长变化,不能提高测量精度。但该方法与选定的监测点无关,无需进行传感器的校准;该方法无需人为判断干涉光谱左移还是右移,无需人为判断是否是周而复始的干涉光谱,有助于实现在线测量;该方法只需包含1个周期的干涉波形,需要的信息量更少,可以实现双折射光纤长度更短的Bi-FLM传感器的在线测量。
欲使(5)式对应干涉光谱T'(λ)的值最大,也可令:
由(18)式可得:
由(19)式可得:
由(10)式=(20)式可得:
由(21)代入(10)式可得:
由(22)式可知,双折射光纤受轴向应变后的绝对长度L'可以由任意一组相邻的波峰波长λ'n、波谷波长λn,光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k求出,再将L'代入εz=ΔL/L0=(L'-L0)/L0计算双折射光纤所受应变大小。
由(19)式可得:
当(14)式、(23)式n值相同时,波谷与波峰相邻。当εz>0时使n为正,满足λ'n>λn,应用(22)式。
比较(13)式、(22)式可知:一般情况下n为正,满足λ'n<λn时使用(13)式;满足λ'n>λn时使用(22)式。若已知任意一组相邻的波峰波长、波谷波长,无需判断波峰波长λ'n、波谷波长λn谁大谁小,(13)式、(22)式统一为(24)式:
由式(24)可知,双折射光纤受轴向应变后的绝对长度L'可以由任意一组相邻的波峰波长λ'n、波谷波长λn,光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k求出,再将L'代入εz=ΔL/L0=(L'-L0)/L0计算双折射光纤所受应变大小。该方法可以由任意一组相邻的波峰波长、波谷波长及光纤初始条件求出双折射光纤所受应变大小。由(13)式和(22) 式可知,双折射光纤受轴向应变后的绝对长度L'与波峰波长λ'n、波谷波长λn有关,即该方法与波谷波长、波峰波长的绝对位置有关,绝对位置与初始相角有关,故该方法无法区分是干扰还是外界传感量的变化导致的波长变化,不能提高测量精度。但该方法与选定的监测点无关,无需进行传感器的校准;该方法无需人为判断干涉光谱左移还是右移,无需人为判断是否是周而复始的干涉光谱,有助于实现在线测量;该方法只需包含1个周期的干涉波形,需要的信息量更少,可以实现双折射光纤长度更短的Bi-FLM传感器的在线测量。
附图说明
图1是双折射光纤环镜传感器原理图。
图2是1550nm波长附近的1组相邻波峰波长和波谷组合干涉光谱图。
具体实施方式
为了更好地理解本发明,下面用具体实例来详细说明本发明的技术方案,但是本发明并不局限于此。
实施例1
本双折射光纤环镜应变传感器在线测量方法,其步骤包括:
1)将3dB光纤耦合器5(3dB couple)、双折射光纤6(Bi-FLM)、光隔离器7(isolator) 以及光谱仪构建成双折射光纤环境轴向应变传感器测量系统,其中,双折射光纤的两端分别与光纤耦合器的两个输出臂相连接,入射光经光隔离器与光纤耦合器的输入臂相连接,光谱仪的输入端与光纤耦合器的干涉光谱输出端口相连接;如图1所示,入射光从端口1 经光隔离器7进入3dB光纤耦合器5,按1:1分成从端口3顺时针和端口4逆时针相向传输的两束光,最后汇聚在端口2,由于双折射光纤6具有双折射效应,汇聚在端口2的两束光发生干涉。当双折射光纤受到应变时,导致双折射光纤双折射率和双折射光纤长度发生变化,从而导致干涉光谱随之改变,以此实现应变测量;
2)将双折射光纤粘贴在被测物体上,记下双折射光纤初始长度L0、双折射光纤初始双折射率B0和双折射应变系数k;
3)待被测物体承受轴向应变后,此时双折射光纤所受应变即为被测物体所受应变,通过光谱仪测量双折射光纤环镜干涉光谱;
4)通过计算机编程找到任意一组相邻的波峰波谷波长λ'n、λn;
5)将光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k代入公式(24)
计算得到双折射光纤受轴向应变后的绝对长度L';
假设双折射光纤长度L0=0.1m,双折射率B0=2.6×10-4,双折射应变系数 k=7.3×0.001/ε,经计算需满足εz<-0.0343915343915344ε,即εz<34391.5343915344με时,才能使n为负数,该应变是个极大的负应变,一般情况下难于产生如此大的负应变,故通常情况下只考虑n为正数的情形。本申请方案也只适用于n为正数的情形。
假设双折射光纤长度L0=0.1m,双折射率B0=2.6×10-4,双折射应变系数k=7.3× 0.001/ε,选取波长范围为典型通讯波长1550nm附近,横坐标λ步长增量设置为0.0001nm,εz=200με时由(4)式绘制的干涉光谱如图2所示。图中包含连续的2个波谷点和2个波峰点,波长由小到大,分别以图中的λ'n+1/λn+1,λn+1/λ'n,λ'n/λn两两组合为一组,组成3组波峰波长、波谷波长,结合光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k,应用(24)式编程求解双折射光纤受轴向应变后的绝对长度L',再将L'代入εz=ΔL/L0=(L'-L0)/L0计算双折射光纤所受应变大小。计算应变及误差如表1所示。从表1可以看出,无需判断波峰波长λ'n、波谷波长λn谁大谁小,统一应用(24)式便可计算双折射光纤受轴向应变大小,计算应变与给定应变基本一致,存在一定误差是由于数字模拟计算的波峰波谷点不是严格意义上的波峰波谷点,而是近似的,从而导致存在一定误差。
定义表1中误差为:
表1 1550nm波长附近的1组相邻波峰波长和波谷波长计算应变的结果
Claims (1)
1.一种一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法,其步骤包括:
1)将光纤耦合器、双折射光纤、光隔离器以及光谱仪构建成双折射光纤环境轴向应变传感器测量系统,其中,双折射光纤的两端分别与光纤耦合器的两个输出臂相连接,入射光经光隔离器与光纤耦合器的输入臂相连接,光谱仪的输入端与光纤耦合器的干涉光谱输出端口相连接;
2)将双折射光纤粘贴在被测物体上,记下双折射光纤初始长度L0、双折射光纤初始双折射率B0和双折射应变系数k;
3)待被测物体承受轴向应变后,此时双折射光纤所受应变即为被测物体所受应变,通过光谱仪测量双折射光纤环镜干涉光谱;
4)通过计算机编程找到任意一组相邻的波峰波谷波长λ'n、λn;
5)将光纤初始长度L0、光纤初始双折射率B0和双折射应变系数k代入公式(24)
计算得到双折射光纤受轴向应变后的绝对长度L';
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910866719.6A CN112484652B (zh) | 2019-09-12 | 2019-09-12 | 一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910866719.6A CN112484652B (zh) | 2019-09-12 | 2019-09-12 | 一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112484652A true CN112484652A (zh) | 2021-03-12 |
CN112484652B CN112484652B (zh) | 2022-05-17 |
Family
ID=74920562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910866719.6A Active CN112484652B (zh) | 2019-09-12 | 2019-09-12 | 一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112484652B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113091781A (zh) * | 2021-03-29 | 2021-07-09 | 北京航空航天大学 | 用于偏光干涉型光纤传感器的解调方法、系统、设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5308973A (en) * | 1990-11-22 | 1994-05-03 | Hilti Aktiengesellschaft | Method and device for the measurement of force by a fiber optics system by evaluating phase shift of light waves |
CN101414093A (zh) * | 2008-11-28 | 2009-04-22 | 南开大学 | 一种基于功能材料填充微结构光纤的可调谐干涉仪 |
CN202066630U (zh) * | 2010-11-19 | 2011-12-07 | 中国计量学院 | 基于差分解调的HiBi-PCF-FLM应力传感器 |
CN103148794A (zh) * | 2013-02-07 | 2013-06-12 | 南京航空航天大学 | 一种基于双折射光纤环境的轴向应变测量方法 |
CN103322934A (zh) * | 2013-06-07 | 2013-09-25 | 南京航空航天大学 | 一种基于双折射光纤环镜的轴向应变测量系统 |
CN106568466A (zh) * | 2016-10-19 | 2017-04-19 | 暨南大学 | 细芯微结构光纤干涉仪传感器及其温度、应变检测方法 |
-
2019
- 2019-09-12 CN CN201910866719.6A patent/CN112484652B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5308973A (en) * | 1990-11-22 | 1994-05-03 | Hilti Aktiengesellschaft | Method and device for the measurement of force by a fiber optics system by evaluating phase shift of light waves |
CN101414093A (zh) * | 2008-11-28 | 2009-04-22 | 南开大学 | 一种基于功能材料填充微结构光纤的可调谐干涉仪 |
CN202066630U (zh) * | 2010-11-19 | 2011-12-07 | 中国计量学院 | 基于差分解调的HiBi-PCF-FLM应力传感器 |
CN103148794A (zh) * | 2013-02-07 | 2013-06-12 | 南京航空航天大学 | 一种基于双折射光纤环境的轴向应变测量方法 |
CN103322934A (zh) * | 2013-06-07 | 2013-09-25 | 南京航空航天大学 | 一种基于双折射光纤环镜的轴向应变测量系统 |
CN106568466A (zh) * | 2016-10-19 | 2017-04-19 | 暨南大学 | 细芯微结构光纤干涉仪传感器及其温度、应变检测方法 |
Non-Patent Citations (1)
Title |
---|
江莺等: "双折射光纤环镜应变灵敏度优化研究", 《光谱学与光谱分析》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113091781A (zh) * | 2021-03-29 | 2021-07-09 | 北京航空航天大学 | 用于偏光干涉型光纤传感器的解调方法、系统、设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112484652B (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dong et al. | High sensitivity optical fiber curvature sensor based on cascaded fiber interferometer | |
CN108895978B (zh) | 一种基于裸光纤的光纤应变传感器灵敏度标定方法 | |
CN102944253B (zh) | 基于偏振测量的光纤光栅横向压力和温度同时测量系统 | |
CN108132067B (zh) | 一种无源温度补偿的保偏光纤干涉仪及其补偿方法 | |
CN108180866B (zh) | 光纤光栅矢量弯曲识别器 | |
CN112525372B (zh) | 基于保偏光纤双臂异轴干涉仪的应变温度同时测量装置及方法 | |
CN110057307B (zh) | 一种提高光纤干涉仪应变灵敏度的方法及光纤干涉仪 | |
CN112484652B (zh) | 一组相邻的波峰波谷波长实现双折射光纤环镜应变在线测量方法 | |
CN209820413U (zh) | 一种光纤干涉仪 | |
CN101344439B (zh) | 一种光纤温度测量装置的构建方法 | |
Jin et al. | Analysis and correction method of axial strain error in multi-core fiber shape sensing | |
CN105486425A (zh) | 一种大量程的温度绝对值测量方法及测量装置 | |
CN110631500B (zh) | 双折射光纤环镜应变传感器在线测量方法 | |
CN110411354B (zh) | 光纤光栅宽量程位移监测装置及系统 | |
CN210036956U (zh) | 一种光纤光栅双边沿校正解调振动测量装置 | |
CN108827189B (zh) | 一种基于反射式微纳光纤耦合器的扭转传感器 | |
CN1851414A (zh) | 消除交叉敏感的光纤布拉格光栅传感系统 | |
CN114964029B (zh) | 一种基于光纤测量技术的轴扭力检测系统 | |
CN113494890B (zh) | 基于fpi干涉仪的光纤光栅应变传感器精度测量装置及方法 | |
Peng et al. | RETRACTED ARTICLE: Intelligent electrical equipment fiber Bragg grating temperature measurement system | |
Zhu et al. | Cross-Sensitivity-Free Highly Sensitive Torsion and Strain Sensor Based on Concatenated DTP-Customizable Helical Fiber Gratings | |
CN209495726U (zh) | 一种折射率和温度同时测量的全光纤传感器 | |
CN109163759B (zh) | 全光纤反射式硝酸根浓度和温度同时测量装置 | |
Zhao et al. | A fast test and compensation system for optical encoders based on extreme learning machine—Fourier Neural Network | |
CN200959040Y (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 |