中国激光, 2019, 46 (1): 0106003, 网络出版: 2019-01-27   

光纤陀螺环热扩散延迟响应模型及补偿技术研究 下载: 659次

Response Model and Compensation Technology of Thermal Diffusion Delay in Fiber Optic Gyro Coil
作者单位
北京航空航天大学仪器科学与光电工程学院, 北京 100191
摘要
基于热扩散机理, 研究了光纤陀螺中光纤环热扩散非互易相位误差导致的零漂的延迟响应特性, 发现了延迟时间的温度相关性, 获得了光纤环零漂延迟补偿模型。搭建高精度光纤陀螺实验系统, 在不同温度段, 采用近似“方波”形台阶温变激励, 验证了热致非互易性相位误差的延迟特性, 并获得了实验光纤环的模型参数。进行了光纤陀螺全温变零偏测试, 实验结果表明采用改进模型可获得更好的零偏补偿效果, 验证了新建模型的正确性和有效性。
Abstract
The delay response characteristics of bias drift caused by the non-reciprocity phase errors of a fiber coil in the fiber optic gyroscope (FOG) are investigated. The correlation between temperature and delay time is revealed and the modified temperature drift compensation model is established. In addition, the high-precision FOG experimental system is built. The delay characteristics of the non-reciprocity phase errors are verified using the approximate "square wave" temperature variation at different temperature ranges, and the model parameters of the experimental fiber coil are obtained as well. The test of FOG bias drift in full-temperature is carried out. The experimental results show that the modified model can used to achieve a better compensation effect based on the thermal diffusion mechanism, and the correctness and validity of the new model are verified.

1 引言

光纤陀螺(FOG)因具有无运动部件、精度覆盖面广、动态范围大等优点。经过40余年的研究与发展,FOG已在海、陆、空、天等领域得到广泛应用[1]。陀螺的性能除了与本身的结构有关之外,与周围环境也有十分密切的关系。其中,影响最大的是光纤环周围温度场的变化,时变的温度场所引起的非互易性相移会使光纤陀螺的零偏产生漂移,即Shupe效应[2]。Shupe效应是光纤陀螺实际应用的主要误差源,被认为是光纤陀螺向高精度发展的最后障碍。

国内外研究学者采用多种方法来抑制Shupe误差。Mohr[3]建立了实用性好和可信度高的温度模型并证实了四极对称绕法是抑制热引起的相位误差的最佳方案;文献[ 4]证明了四极对称绕法可使Shupe误差减少三个数量级;文献[ 5]证实改善光纤陀螺结构和组成部件,采用绝缘材料包围的光纤,设计传感器的金属外壳可进一步减小误差;光纤陀螺电路部分各器件发热不同,会影响陀螺工作的热平衡状态;文献[ 6]提出了一种基于分级控制、分段控制和闭环控制思想的温控方案,取得了较好的误差抑制效果;除上述方法之外,对光纤陀螺温度特性进行实验研究,建立温度漂移模型并实施温度补偿也是工程上解决光纤陀螺温度漂移问题的有效手段[7-8]

本文主要在光纤陀螺光纤环Shupe效应产生机理的基础上,基于热扩散原理,建立光纤环热致非互易性相移误差分析模型,分析Shupe效应产生的非相移误差与热扩散延迟响应时间、外界环境温度的相关性,建立基于光纤环热扩散延迟响应时间的温度补偿模型。搭建高精度光纤陀螺实验系统,在不同温度下,引入近似“方波”台阶温变激励,测试相移误差的延迟和温度相关性,验证模型的准确性并获得不同温度下的温度补偿模型参数。进行全温变实验,对比分析传统温度补偿模型与本文新建模型的误差补偿效果。

2 光纤环热致非互易相位误差模型

图1为光纤环展开示意图,光纤环采用对称(如四极对称)绕环技术绕制,长度为L,J为分光点,M为环中点,光纤环关于线段JM对称。在温度场中,光纤环对称部分经历的温度近似相同,如图1中颜色梯度所示。AA'段为光纤环中任意一对对称的光纤段,长度为dz,以J点为起点,沿顺时针方向,AA'段的位置分别为zL-z, T·(z)和 T·(L-z)为AA'点在相同时刻的温变速率,在光纤环中由这两段对称光纤产生的热致非互易性相位误差为[2]

dφe=2πλαSiO2+ndndTL-2zcv/n[T·(z)-T·(L-z)]dz,(1)

式中:αSiO2为石英材料的热膨胀系数,一般为0.5×10-6/℃;cv为光速;n为光纤有效折射率;dn/dT表示折射率随温度变化的系数;λ为工作波长。(1)式为经典型Shupe效应模型。目前,相关的误差补偿模型都集中在(1)式的温度变化因子[ T·(z)-T·(L-z)]上,由此开发了热致非互易性相位误差的模型补偿算法。但是,光纤环实际还存在热应力效应[9],光纤环材料热扩散等参数可能与作用温度值相关,温度及温度变化会引起dn/dT、热膨胀系数αSiO2的改变。为获得更准确的补偿模型,基于热扩散理论,对光纤环热致非互易性相位特性进行研究。

图 1. 光纤环展开示意图

Fig. 1. Schematic of fiber coil expansion

下载图片 查看所有图片

图 2. 光纤环温度分布图。(a)整体分布;(b)等效半无限热扩散模型

Fig. 2. Temperature distribution of fiber coil. (a) Overall distribution; (b) equivalent semi-infinite thermal diffusion model

下载图片 查看所有图片

一般情况下,可认为光纤环内外的温度场是均匀的,如图2(a)所示,阴影部分代表四极对称缠绕的光纤环,热扩散方向可能是由外往里或由里向外。为了简化模型,可认为环内等直径位置处温度相等。温度沿环径向外扩散,光纤环可等效为半无限的热扩散模型,如图2(b)所示。设图1AA'光纤段对应图2(b)中扩散深度为d的位置,基于热扩散理论,热传导方程可表示为[10]

T(d,t)t-K2T(d,t)d2=0T(0,t)=T0,(2)

K=rc0ρ,(3)

式中,T0为环境初始温度,K为材料的热扩散系数,r为材料的导热系数,ρ为材料的密度,c0为比热容。t时刻扩散深度为d处的温度分布为

T(d,t)=T0·erfcd2Kt,(4)

式中,误差函数为

erfc(x)=2πx+exp-y2)dy(5)

由(4)式可知,光纤环内的温度变化为非线性响应,且与热扩散系数K有关。设t=0时,环境温度由T0变化至T1,光纤环内d位置的归一化温度值与时间t的关系可表示为

T(d,t)-T0=(T1-T0)1-erfcd2Kt(6)

定义τc=d2/K为环中位置d处扩散时间参数,表示该点温度变化为设定温度的一半时所需要的时间。

一般情况下,材料的导热系数、比热容等参数都与材料本身温度相关,导致材料的热扩散系数随材料温度的变化而变化,即材料扩散系数K不但与材料导热特性相关,还与材料的工作温度相关。文献[ 11]测量了不同温度下石英材料的导热系数,结果表明,石英材料的温度越高,热扩散系数越大。设光纤环在不同温度条件下的热扩散系数分别为K1K2K3,且K1>K2>K3,以光纤材料的热扩散系数为输入,基于(6)式,可仿真光纤环中位置d处温度的变化曲线,如图3所示。由图3可知,外界环境温度由T0阶跃变化至T1时(蓝色方波实线),光纤环内部的温度变化为非线性响应(T_Ki表示热扩散系数为Ki时的响应曲线),存在明显延迟,而且延迟时间与环境温度、光纤环热扩散系数密切相关,热扩散系数K越大,时间常数τc越短。

图 3. 不同扩散系数下光纤环位置d处的温度变化示意图

Fig. 3. Schematic of temperature change at position d in fiber coil under different diffusion coefficients

下载图片 查看所有图片

(1)式中 T·(z)和 T·(L-z)为AA'点在相同时刻的温变速率,实际的温度补偿技术中,无法直接测到光纤环结构内的温度,需要通过测量环外温度并建立数学模型进行补偿。目前,大部分温度补偿的建模思路都是基于IEEE标准(IEEE Std 952-1997)给出的单轴光纤陀螺仪输入角速率与陀螺仪输出之间的模型方程[12]

D=cT·T+cT··T·,(7)

式中:D为由环境温度变化引起的零偏漂移;T为环境温度; T·表示温变速率,即为dT/dt;cTcT·分别为温度敏感系数与温变速率敏感系数。为获得更精确的实用温度补偿模型,一般将光纤陀螺零偏与温度近似为线性关系,并加入环境温度与温变速率的交叉项[13],形成实际应用的模型方程:

D'(t)=[aT(t)T·(t)+bT(t)+cT·(t)+d]·D0(t),(8)

式中,D0为陀螺零偏的实测值,D'为陀螺零偏的估计值,T(t)为t时刻陀螺温度, T·(t)为t时刻的温变速率,abcd为模型参数。

显然,这些模型补偿技术未考虑光纤环温度响应的非线性和时间延迟,很难实现好的补偿效果。(6)式和图3表明,光纤环的温度变化属于热扩散过程,环内温度变化存在明显的延迟,温变引起的陀螺漂移相对环外温度变化也存在延迟。设陀螺零漂与环外温度变化的延迟时间为τd,在光纤环中由图1两对称光纤产生的热致非互易性相位误差可重写为

dφe(t)=2πλαSiO2+ndndT×L-2zcv/n[T·(z,t-τd)-T·(L-z,t-τd)]dz,(9)

式中 T·(z)和 T·(L-z)为由(6)式根据环外温度计算得到的AA'点的温变速率。对应的温度补偿模型[(8)式]也可重写为

D'(t)=[aT(t-τd)T·(t-τd)+bT(t-τd)+cT·(t-τd)+d]·D0(t)(10)

在实际光纤陀螺测试实验中,光纤环内位置d处的温度及温度的变化率很难用类似(6)式的表达式精确计算,实际的延迟时间常数τc和延迟时间τd以及模型系数abcd需通过温度实验的测试数据拟合获得。

3 温度补偿模型验证及参数标定

为验证上述算法模型,并对不同温度下的温度误差模型参数进行标定,搭建实验系统进行近似“方波”形温度激励实验。实验系统结构框图和照片如图4图5所示,实验系统为光电分离的高精度光纤陀螺结构,由信号处理单元(SPU)和Sagnac光纤环路单元(SIA)组成。SPU包括宽带光源、耦合器、检测器、数字闭环信号检测模块、温度测量模块和用于数据采集和分析的计算机;SIA由多功能集成光学调制器(MIOC)、保偏光纤环(PMFC)和Pt温度传感器组成,光纤环采用四极对称放式进行缠绕,长度约为3 km,平均直径为150 mm,SIA被封装在一个圆形铝盒中。SIA通过光电混合缆与SPU连接,形成光电分离闭环光纤陀螺系统。实验时,只有SIA盒被安装在带温箱的高精度转台上。实验过程中,在-45~65 ℃范围内,分别选取低、中、高三个温度段,进行温度激励。具体变温过程为-25 ℃ →-45 ℃→-25 ℃、10 ℃ →-10 ℃→10 ℃、45 ℃→65 ℃→45 ℃。

图 4. 实验系统结构框图

Fig. 4. Structural diagram of experimental system

下载图片 查看所有图片

图 5. 实验照片。(a)系统照片;(b)带温箱转台安装照片

Fig. 5. Experimental photos. (a) System photo; (b) installation photo of turntable with temperature chamber

下载图片 查看所有图片

在变温实验,即进行-25 ℃→-45 ℃→-25 ℃近“方波”温度激励实验时,温箱温度变化程序为(-25 ℃,180 min)→(-25 ℃→-45 ℃,20 min)→(-45 ℃,180 min)→(-45 ℃→-25 ℃,20 min)→(-25 ℃,180 min)。在升降温过程中温箱温变速率约为1.0 ℃/min,在高、低温段保持恒温时间较长,均为3 h,以获得稳态时的实验数据。在整个实验过程中,温箱上电工作,但转台不上电,避免转台激磁电流引起的磁场对陀螺零偏的影响。实时同步记录光纤陀螺的输出和SIA盒外环境和盒内光纤环温度,采样率为1 Hz。根据高精度光纤陀螺的数据处理和测试规范,对输出数据进行100 s平滑。测试曲线如图6所示,在同一曲线图中,显示设定(温箱)温度TS(蓝色点划线)、SIA盒内光纤环温度TC(橘色虚线)、陀螺输出ΩB(黑实线)、盒内温度的微分 T·C(t)(粉色点线)和拟合后带延迟的温度微分 T·C(t-τd)(绿色短虚线)。

图 6. -25 ℃→-45 ℃→-25 ℃,陀螺输出结果图

Fig. 6. Output diagrams of FOG for -25 ℃→-45 ℃→-25 ℃

下载图片 查看所有图片

盒内温度的微分 T·C(t)为第tt时刻的温度与第t时刻温度的差值与时间Δt之比,计算方程为

T·C(t)=TC(t+Δt)-TCtΔt(11)

分别以(8)式和(10)式为目标函数进行拟合,获得无延迟模型系数abcd和有延迟模型系数abcdτd,如表1表2所列。采用同样的实验步骤和数据处理方法,获得其他两个温度段的模型参数,拟合得到的模型参数在表1表2中列出。由表1和2可知,在同种情况下,模型系数abcd在量级上是一致的,模型参数具有较好的一致性。

表 1. 无延迟模型参数值

Table 1. Parameter values of no-delay model

T /℃abcd
-25→-45→-25-4.89×10-6-3.08×10-5-0.0002930.99998
10→-10→10-6.12×10-6-1.29×10-5-0.0002781.00012
45→65→45-1.24×10-5-3.22×10-5-0.0005311.00141
Average value-7.80×10-6-2.53×10-5-0.0003671.00017

查看所有表

表 2. 有延迟模型参数值

Table 2. Parameter values of delay model

T /℃abcdτd /(102 s)τ'd /(102 s)
-25→-45→-25-5.11×10-6-3.45×10-5-0.0003281.000318.638.12
10→-10→10-7.66×10-6-1.94×10-5-0.0003621.000435.475.24
45→65→45-1.95×10-5-3.75×10-5-0.0006781.001023.093.21
Average value-1.08×10-5-2.87×10-5-0.0004561.00059

查看所有表

图 7. -25 ℃→-45 ℃→-25 ℃,陀螺原始零偏和补偿后零偏

Fig. 7. Original zero bias and compensated zero bias of gyro for -25 ℃→-45 ℃→-25 ℃

下载图片 查看所有图片

表 3. 补偿前后的陀螺零漂

Table 3. Zero biases of gyro before and after compensation

T /℃Original biasdrift ΩB /[(°)·h-1]Bias driftwithoutdelay ΩC /[(°)·h-1]Bias driftwith delayΩD /[(°)·h-1]
-25→-45→-250.00710.00370.0023
10→-10→100.00420.00300.0021
45→65→450.01130.00550.0031

查看所有表

图 8. 热扩散延迟响应时间与温度关系曲线图

Fig. 8. Relationship between thermal diffusion delay response time and temperature

下载图片 查看所有图片

4 全温变实验验证

为了进一步对上述理论模型和算法进行验证,进行了-45~ 65 ℃的全温变光纤陀螺零偏漂移测试和补偿实验,采用相同的实验装置,将SIA盒安装在带温箱的高精度转台上,温箱温度变化程序为(-25 ℃,60 min )→(-25 ℃→-45 ℃,40 min)→(-45 ℃,120 min)→(-45 ℃→65 ℃,220 min)→(65 ℃,120 min) →(65 ℃→25℃,80 min)→(25 ℃,30 min)。在升降温过程中温箱温速率约为0.5 ℃/min,在-45 ℃和65 ℃的恒温时间设置为2 h,以验证恒温时的补偿效果。在整个实验过程中,实时同步记录光纤陀螺的输出和SIA盒外环境温度和盒内光纤环温度,采样率为1 Hz,同样对输出数据进行100 s平滑,测试曲线如图9所示。

图 9. 全温段内陀螺输出结果图

Fig. 9. Gyroscope output in full-temperature

下载图片 查看所有图片

分别以传统温度补偿模型和本研究提出的有延迟补偿模型对陀螺原始输出进行补偿,其中模型参数abcd分别采用表1表2中计算的平均值 a-b-c-d-,有延迟补偿模型中的不同温度下的延迟时间由(12)式计算。补偿得到的零偏漂移曲线如图10所示,图10ΩB(黑色实线)为陀螺的原始输出、ΩC(蓝色点线)为传统无延迟模型补偿输出、ΩD(红色虚线)为有延迟模型补偿输出。计算得到的标准偏差分别为0.0198,0.0093,0.0046 (°)/h。

图 10. 全温段原始陀螺输出与补偿后陀螺输出

Fig. 10. Original gyro output and compensated gyro output in full temperature

下载图片 查看所有图片

实验证明,热扩散延迟响应时间与温度呈相关性,带有热扩散延迟响应时间的模型算法可将3 km的光纤环的零漂控制在0.005 (°)/h以内。带有热扩散延迟响应时间的补偿模型能够更好地补偿温变引入的零漂误差。

5 结论

通过将热扩散机理引入光纤环温度响应过程中,研究了光纤陀螺中光纤环热致非互易相位误差导致的零偏漂移的延迟响应特性,发现了延迟时间的温度相关性,改进了光纤环零漂延迟补偿模型,提出了基于多项式拟合的带有热扩散延迟响应时间的零偏温度漂移补偿技术。搭建高精度光纤陀螺实验系统,并进行三种不同温度下的“方波”台阶激励实验,通过非线性拟合得到不同温度下的热扩散延迟响应时间,完成了最优模型参数和不同温度下延迟时间的辨识。进行了(-45~65 ℃,0.5 ℃/min)的全温变实验验证,对比了传统模型与新建模型的补偿效果,实验结果验证了新建模型的正确性和有效性,同时表明采用新建模型可获得更好的零偏补偿效果,可将3 km的光纤环的零漂控制在0.005 (°)/h以内,为高精度光纤陀螺提供了实用有效的温度误差补偿技术。

参考文献

[1] Bergh R, Lefevre H, Shaw H. An overview of fiber-optic gyroscopes[J]. Journal of Lightwave Technology, 1984, 2(2): 91-107.

[2] Shupe D M. Thermally induced nonreciprocity in the fiber-optic interferometer[J]. Applied Optics, 1980, 19(5): 654-655.

[3] Mohr F. Thermooptically induced bias drift in fiber optical Sagnac interferometers[J]. Journal of Lightwave Technology, 1996, 14(1): 27-41.

[4] Lofts C M. Investigation of the effects of temporal thermal gradients in fiber optic gyroscope sensing coils[J]. Optical Engineering, 1995, 34(10): 2856.

[5] Tirat O F J, Euverte J M. Finite element model of thermal transient effect in fiber optic gyro[J]. Proceedings of SPIE, 1996, 2837: 230-239.

[6] KatoH, SaitoH, Sakai SI. Shupe effect compensation of temperature controlled fiber optical gyroscope[C]∥2010 11th IEEE International Workshop on Advanced Motion Control (AMC), March 21-24, 2010, Nagaoka, Niigata, Japan. New York: IEEE, 2010: 262- 267.

[7] 赵建辉, 徐帆, 李帆. 随钻测斜仪中光纤陀螺光纤环Shupe效应仿真[J]. 光学学报, 2018, 38(5): 0506001.

    Zhao J H, Xu F, Li F. Simulation of Shupe effect in fiber optic gyroscope fiber coil with inclinometer while drilling[J]. Acta Optica Sinica, 2018, 38(5): 0506001.

[8] 李绪友, 凌卫伟, 许振龙, 等. 双柱型绕法对干涉式光纤陀螺温度性能的影响[J]. 光学学报, 2016, 36(8): 0806003.

    Li X Y, Ling W W, Xu Z L, et al. Effect of double-cylinder winding on temperature performance of interferometric fiber optic gyroscope[J]. Acta Optica Sinica, 2016, 36(8): 0806003.

[9] Minakuchi S, Sanada T, Takeda N, et al. Thermal strain in lightweight composite fiber-optic gyroscope for space application[J]. Journal of Lightwave Technology, 2015, 33(12): 2658-2662.

[10] Mohr F. Bias error in fiber optic gyroscopes due to elasto-optic interactions in the sensor fiber[J]. Proceedings of SPIE, 2004, 5502: 410-413.

[11] 金太权. 测定石英玻璃的导热系数[J]. 吉林工业学院学报, 1991, 12(1): 39-42.

    Jin T Q. Determination of thermal conductivity of quartz glass[J]. Journal of Jilin Institute of Technology, 1991, 12(1): 39-42.

[12] IEEE Std 952-1997. IEEE standard specification format guide and test procedure for single-axis interferometric fiber opticgyros Std 952-1997. IEEE standard specification format guide and test procedure for single-axis interferometric fiber opticgyros[S]. IEEE Aerospace and Electronic Systems Society, 1997

[13] 柳贵福, 邢艳丽, 张树侠. 光纤陀螺零偏稳定性的数据建模方法研究[J]. 中国惯性技术学报, 2001, 9(3): 48-52.

    Liu G F, Xing Y L, Zhang S X. Investigation on the modeling of fiber optic gyro bias instability[J]. Journal of Chinese Inertial Technology, 2001, 9(3): 48-52.

闫晗, 杨远洪, 杨福铃. 光纤陀螺环热扩散延迟响应模型及补偿技术研究[J]. 中国激光, 2019, 46(1): 0106003. Yan Han, Yang Yuanhong, Yang Fuling. Response Model and Compensation Technology of Thermal Diffusion Delay in Fiber Optic Gyro Coil[J]. Chinese Journal of Lasers, 2019, 46(1): 0106003.

本文已被 5 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!