2. 哈尔滨工程大学 自动化学院, 黑龙江 哈尔滨 150001
2. College of Automation, Harbin Engineering University, Harbin 150001, China
由于光纤陀螺具有精度高、可靠性好等优点,相比于传统机械转子陀螺仪,其具有无转动部件、结构简单、成本低、便于批量生产等特性[1-2]。在军事和民用范围具有广泛的应用价值以及应用前景。光纤陀螺的硬件部分由光路和电路2部分组成。光路部分的光纤环作为敏感光纤陀螺角速率输出的重要元件,当环境温度发生变化时,光纤环的性能会受到影响[3],为降低温度变化导致的非互易性,一般采用对称绕法[4]。除此之外,在绕制光纤环时由于工艺原因还会在光路中引入应力双折射点,产生非互易性误差,导致光纤陀螺精度不能达到预期要求。文献[5]中建立光纤陀螺光路数学模型时,默认光纤环中仅存在一个偏振串扰点,故建立模型与实际光纤环有区别,与实际得到的相位误差不符。文献[6]在建立光纤陀螺光路的数学模型时未考虑不同温度情况下光路中器件的变化,与工程应用中的实际环境不符,不具有一般性。
本文充分考虑光路中的熔接点、光学器件参数,通过相干矩阵表示输入光信号特性、琼斯矩阵表示光学器件的传输特性以及熔接点。利用麦克尔逊干涉仪检测光纤环中偏振串扰点位置和强度。完善光纤陀螺光路模型,使之适用于工程环境下。为验证模型的准确性,本文通过理论仿真分析和实验结果分析温度对光纤陀螺静态零偏的影响,证实本文提出的方法能够作为光纤环质量评估手段。
1 光路部件及光纤陀螺光路建模光纤陀螺误差源有2大部分:1)由于光路中各器件尾纤相接时,存在对轴误差角度,导致光信号传输至熔接点处时将发生偏振串扰现象;2)由于在绕环时工艺水平的限制以及操作不当引入的偏振串扰点。除此之外当温度环境发生变化时,这些偏振串扰点的位置和强度也随之变化。
对光纤陀螺光路中的光学器件进行建模[7],系统分析偏振串扰误差对光纤陀螺输出精度的影响。设定光波向光纤环传播的方向为顺时针,反之为逆时针,不考虑电路系统对光路的影响。光纤陀螺光路部分包括SLD宽带光源、光纤耦合器、光纤环、Y波导和PIN-FET探测器。根据光传播的途径对光学器件建模。文献[8]可以得到与光源偏振相关的光源相干矩阵模型;文献[9]可以得到光纤耦合器的琼斯矩阵;文献[10]可以得到与Y波导中起偏器的消光比ε相关的光纤偏振器琼斯矩阵模型以及加±π/2的相位调制时相位调制器的琼斯矩阵模型。
在目前相关的文献中,建立光纤环琼斯矩阵时,一般假设光纤环中仅存在一个串扰点[11],但在实际情况中光纤环内存在大量随机分布的偏振串扰点,且当环境温度变化时串扰点的位置和强度均发生变化。因此以往文献在通过光路模型分析光纤环静态零偏时与实际情况有所偏差。
本文通过麦克尔逊干涉仪测量光纤环中存在N个串扰点,并得到这些串扰点的位置和强度信息[12]。根据串扰点的位置信息,将2个串扰点之间的光纤长度等效为2个熔接点之间的尾纤长度;根据串扰点的强度信息等效为角度,将串扰点转化为熔接点。根据上述信息最终建立与温度相关的光纤环琼斯矩阵表达式:
$ \begin{array}{l} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\mathit{\boldsymbol{T}}_{{\rm{Coil\_CW}}}} = \prod\limits_{i = 1}^N {\left[ {\begin{array}{*{20}{c}} 1&0\\ 0&{{\rm{exp}}\left[ { - {\rm{i}}\frac{{2\pi \Delta {n^\prime }(T){l_i}}}{\lambda }} \right]} \end{array}} \right]} \cdot \\ \left[ {\begin{array}{*{20}{c}} {{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}}&{{\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}}\\ { - {\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}}&{{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}} \end{array}} \right] \cdot \left[ {\begin{array}{*{20}{c}} 1&0\\ 0&{ - {\rm{i}}\frac{{2\pi \Delta {n^\prime }(T)({l_{i + 1}} - {l_i})}}{\lambda }} \end{array}} \right] \cdot \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} {{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}}&{{\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}}\\ { - {\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}}&{{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}} \end{array}} \right] \cdot \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} 1&0\\ 0&{{\rm{exp}}\left[ { - {\rm{i}}\frac{{2\pi \Delta {n^\prime }(T)(L - {l_N})}}{\lambda }} \right]} \end{array}} \right] \end{array} $ | (1) |
$ \begin{array}{l} {\kern 1pt} {\mathit{\boldsymbol{T}}_{{\rm{Coil\_CCW}}}} = \left[ {\begin{array}{*{20}{c}} 1&0\\ 0&{{\rm{exp}}\left[ { - {\rm{i}}\frac{{2\pi \Delta {n^\prime }(T)\left( {L - {l_N}} \right)}}{\lambda }} \right]} \end{array}} \right] \cdot \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \prod\limits_{i = N}^1 {\left[ {\begin{array}{*{20}{c}} 1&0\\ 0&{ - {\rm{i}}\frac{{2\pi \Delta {n^\prime }(T)({l_{i + 1}} - {l_i})}}{\lambda }} \end{array}} \right] \cdot } \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} {{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}}&{{\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}}\\ { - {\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}}&{{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _{i + 1}}} \end{array}} \right] \cdot \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} 1&0\\ 0&{{\rm{exp}}\left[ { - {\rm{i}}\frac{{2\pi \Delta {n^\prime }(T){l_i}}}{\lambda }} \right]} \end{array}} \right] \cdot \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} {{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}}&{{\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}}\\ { - {\rm{sin}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}}&{{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\theta _i}} \end{array}} \right] \end{array} $ | (2) |
式中:TCoil_CW和TCoil_CCW表示光信号在光纤环中顺时针和逆时针传播时的琼斯矩阵;
Download:
|
|
由干涉仪的工作原理可知在最大干涉条纹包络2侧对称位置各出现一个由两波列产生的干涉条纹,由此计算出偏振光的耦合强度为:
$ h = 10 {\rm{lb}} \left( {\frac{{{I_{{\rm{cm}}}}}}{{{I_{\rm{m}}}}}} \right) $ | (3) |
式中:Icm为干涉次峰包络的幅值;Im为干涉主峰包络的幅值。
且当麦克尔逊干涉仪补偿引入的相位误差后,对应的干涉包络,耦合强度分别为:
$ \left\{ {\begin{array}{*{20}{l}} {{h_1} = \frac{1}{4}{{({\rm{cos}}{\kern 1pt} {\kern 1pt} 2\beta {\rm{sin}}{\kern 1pt} {\kern 1pt} (2\theta ))}^2}}\\ {{h_2} = \frac{1}{4}{{({\rm{sin}}{\kern 1pt} {\kern 1pt} 2\beta {\rm{cos}}{\kern 1pt} {\kern 1pt} (2\theta ))}^2}}\\ {{h_3} = \frac{1}{4}{{({\rm{sin}}{\kern 1pt} {\kern 1pt} 2\beta {\rm{sin}}{\kern 1pt} {\kern 1pt} (2\theta ))}^2}} \end{array}} \right. $ | (4) |
式中β表示起偏角度。当β=0时,干涉峰只剩下h1项,由式(4)可知,耦合光的强度是由对轴角度θ决定的,因此能够通过串扰强度得到光纤环内串扰点的对轴角度θ。
根据光纤陀螺光路原理图中光信号传输路径以及上述光路中串联的各光学器件的琼斯矩阵,建立光信号在光路中顺向和逆向传播时的琼斯矩阵分别为:
$ \begin{array}{l} {\mathit{\boldsymbol{G}}_{CW}} = \mathit{\boldsymbol{K}}_{C1}^{\rm{T}}{\mathit{\boldsymbol{P}}^{\rm{T}}}\mathit{\boldsymbol{K}}_{C2}^{\rm{T}}{\mathit{\boldsymbol{M}}_{CW}}{{\rm{e}}^{{\rm{j}}{\phi _s}/2}}{\mathit{\boldsymbol{T}}_{{\rm{Coil\_CW}}}}{\mathit{\boldsymbol{K}}_{T2}}\mathit{\boldsymbol{P}}{\mathit{\boldsymbol{K}}_{{T1}}} = \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{l}} {{S_{11}}}&{{S_{12}}}\\ {{S_{21}}}&{{S_{22}}} \end{array}} \right]{{\rm{e}}^{{\rm{j}}{\phi _s}/2}} \end{array} $ | (5) |
$ \begin{array}{l} {\mathit{\boldsymbol{G}}_{CCW}} = \mathit{\boldsymbol{K}}_{C1}^{\rm{T}}{\mathit{\boldsymbol{P}}^{\rm{T}}}\mathit{\boldsymbol{K}}_{T2}^{\rm{T}}{{\rm{e}}^{{\rm{ - j}}{\phi _s}/2}}{\mathit{\boldsymbol{T}}_{{\rm{Coil\_CCW}}}}{\mathit{\boldsymbol{M}}_{CCW}}{\mathit{\boldsymbol{K}}_{C2}}\mathit{\boldsymbol{P}}{\mathit{\boldsymbol{K}}_{T1}} = \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{l}} {{T_{11}}}&{{T_{12}}}\\ {{T_{21}}}&{{T_{22}}} \end{array}} \right]{{\rm{e}}^{{\rm{ - j}}{\phi _s}/2}} \end{array} $ | (6) |
式中:
$ \begin{array}{l} \mathit{\boldsymbol{G}} = {\mathit{\boldsymbol{G}}_{CW}} + {\mathit{\boldsymbol{G}}_{CCW}} = \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{l}} {{S_{11}}{{\rm{e}}^{{\rm{j}}{\phi _s}/2}} + {T_{11}}{{\rm{e}}^{ - {\rm{j}}{\phi _s}/2}}}&{{S_{12}}{{\rm{e}}^{{\rm{j}}{\phi _s}/2}} + {T_{12}}{{\rm{e}}^{ - {\rm{j}}{\phi _s}/2}}}\\ {{S_{21}}{{\rm{e}}^{{\rm{j}}{\phi _s}/2}} + {T_{21}}{{\rm{e}}^{ - {\rm{j}}{\phi _s}/2}}}&{{S_{22}}{{\rm{e}}^{{\rm{j}}{\phi _s}/2}} + {T_{22}}{{\rm{e}}^{ - {\rm{j}}{\phi _s}/2}}} \end{array}} \right] \end{array} $ | (7) |
则陀螺输出光信号的相干矩阵为:
$ \mathit{\boldsymbol{J}} = \mathit{\boldsymbol{G}}{\mathit{\boldsymbol{J}}_{{\rm{source}}}}{\mathit{\boldsymbol{G}}^{\rm{H}}} $ | (8) |
探测器检测到光强为:
$ {\mathit{\boldsymbol{I}}_{{\rm{OUT}}}} = {\rm{Tr}} [{\mathit{\boldsymbol{J}}_{{\rm{OUT}}}}] $ | (9) |
若通过方波调制,即设置相位调制器施加±π/2相位。当光信号传播至相位调制器时,调制相位满足:
$ {\phi _m}(t) - {\phi _m}(t - \tau ) = \pm \pi /2 $ | (10) |
系统光路输出光强表示为:
$ {I_{{\rm{ out }}}} = {I_0} + N{\rm{sin}}{\kern 1pt} {\kern 1pt} {\phi _s} + M{\rm{cos}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\phi _s} = {I_0} + {I_P}{\rm{sin}}({\phi _s} + {\varphi _e}) $ | (11) |
式中:I0表示输出光强的常数项; Ip表示偏振串扰导致的光强变化; φe表示由偏振串扰导致的相位误差,根据式(11)偏振相位误差可以表示为:
$ \begin{array}{*{20}{l}} {{\varphi _e} = {\rm{arctan}}\left( {\frac{M}{N}} \right) = {\rm{arg}}[({T_{11}}S_{11}^* + {T_{21}}S_{21}^*){C_{xx}} + }\\ {({T_{12}}S_{11}^* + {T_{22}}S_{21}^*){C_{yx}} + ({T_{11}}S_{12}^* + {T_{21}}S_{22}^*){C_{xy}} + }\\ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ({T_{12}}S_{12}^* + {T_{22}}S_{22}^*){C_{yy}}]} \end{array} $ | (12) |
本文以建立的光路模型为依据,结合式(11)和(12),分析各光学器件的性能对偏振相位误差的影响。
2.1 光源偏振度P与偏振相位误差若不考虑部分损耗和耦合,设光纤耦合器分光比为1:1,偏振器消光比ε=0.025,光源假定为线偏振光,椭圆方位角ψ=0.1 rad。光源偏振度P与偏振串扰相位误差φe的关系如图 1所示。
根据图 2仿真结果,当光源偏振度P=0时,对应相位误差φe=0,且相位误差φe与光源偏振度P呈近似正比的关系。为减小由于光源偏振度导致的相位误差可以考虑选用非偏振光源。
Download:
|
|
当不考虑光纤陀螺光路中耦合且假设各器件尾纤熔接点无失准角。由于实际情况下Y波导中集成偏振器的消光比ε不可能为零。假设
如图 2所示,在同一的转速下,消光比ε增大,输出光强Iout也增大且在ε逐渐增大的过程中Iout的增长速度变快。不同转速下,消光比ε与输出光强Iout变化规律大体相同。
2.3 光纤环内偏振串扰点等效通过麦克尔逊干涉仪测量出光纤环中偏振串扰信息,当建立不同温度下光纤陀螺光路模型时,可以将光纤环中的偏振串扰点等效为纤与纤之间的熔接点,通过已知的偏振串扰强度计算两段光纤偏振面旋转角度即对轴角度θ,再根据偏振面旋转角度对偏振相位误差的影响进行分析,最终得到偏振相位误差与偏振串扰强度之间的关系。
通常情况下β=0,故式(3)中串扰强度主要由h1项组成,则h1项串扰强度和对轴角度α之间的关系如图 3所示。
Download:
|
|
根据建立的光路模型和偏振相位误差表达式(12),推导出光纤环中偏振串扰强度等效的对轴角度θ与光纤陀螺偏振相位误差的关系如图 4所示。
Download:
|
|
图 4中的曲线表明,陀螺输出相位误差随偏振串扰强度等效的对轴角度θ的变化,呈近似正比关系,当对轴角度θ=0时相位误差也为0,综上光纤环中的偏振串扰强度越小,陀螺输出偏振相位误差越小。
3 温度实验结果及分析取实验室绕制的2个光纤环设置编号为#1和#2,进行高低温实验。参考GJB2 426 A-2004[13],得到在不同温度时光纤陀螺光路中由于偏振串扰相位误差而引起的陀螺零偏:
$ {B_{\rm{s}}} = \frac{1}{K}{\left[ {\frac{1}{{n - 1}}\sum\limits_{i = 1}^n {{{({\phi _i} - \bar \phi )}^2}} } \right]^{1/2}} $ | (13) |
式中:Bs表示光纤陀螺的零偏,(°)/h;n为采样点数;
将待测光纤环放置在温箱内,设置温度为-40 ℃~60 ℃匀速升温,每隔20 ℃,待温箱温度稳定后,通过麦克尔逊干涉仪测试并记录待测光纤环中偏振串扰点位置和强度。表 1和表 2为统计光纤环1#和2#在不同温度情况下偏振串扰强度的均值、最大值以及串扰强度大于-40 dB的个数。
对比表 1和表 2统计的测试结果,在不同温度点下光纤环1#和光纤环2#偏振串扰强度的均值和最大值均出现规律性变化,在温度范围端点光纤环中出现高串扰点的数目偏高,尤其在低温-40 ℃时,偏振串扰强度大于-40 dB的个数明显高于其他温度。出现该现象是由于,光纤环内部固胶在高温或低温情况下发生不可预测的体积变化,导致纤与纤之间出现挤压,光纤受到应力作用产生高串扰点。
除光纤环外,建立的光纤陀螺光路模型中各光学器件的琼斯矩阵中的元素,根据实验室用光纤陀螺的参数进行设置,光源偏振度P=1;输入光强度为1,Y波导中消光系数ε=0.025,消光比为35 dB;各器件间的尾纤为1 m;光源的波长为1 550 nm;设置光纤环偏振串扰测试温度为-40,-20,0,20,40,60,将测得的偏振串扰强度和位置的数据,代入光纤环琼斯矩阵中进行计算。
最终,根据模型得到光纤环1#和2#在不同温度点下的静态零偏预估值如表 3所示,可以认为光纤环2#在全温下的性能优于1#。
为了验证温度对陀螺性能的影响,将光纤环1#和2#接入光纤陀螺系统。采集常温(20 ℃)状态下时陀螺的静态输出值,通过式(13)计算光纤环1#和2#的静态零偏。
设置温箱温度变化范围-40 ℃~60 ℃~-40 ℃,变化速率为20 ℃/h。图 6表示温度匀速变化的情况下陀螺输出数据随温度变化的曲线,通过实验证明随着温度的变化陀螺输出呈趋势性变化。
对比表 4陀螺静态零偏,在变温情况下2个光纤环的静态零偏明显高于常温情况下。通过温度实验证明在建立光纤环琼斯矩阵模型时,引入温度变化量是正确且必要的。
Download:
|
|
设置起始温度为-40 ℃,每升温20 ℃待温箱温度稳定后,采集该温度下2 h的光纤陀螺静态输出数据,直至升温到60 ℃。根据获得的陀螺输出数据计算出光纤环1#和2#接入的陀螺系统的静态零偏,如表 5所示。
对比表 3和表 5的数据,光纤环1#和2#预估静态零偏值和实验测试静态零偏值,说明根据建立的光路模型能够得到与实验测试静态零偏较为一致的数据。因此认为根据建立的光路模型得到的预估结果能够大致反映光纤环接入系统后的实际测试结果,证明了模型具有合理性。根据光纤陀螺温度实验得到的数据,光纤环中偏振串扰的均值低,高串扰点少,由此光纤环组成的光纤陀螺静态性能好,反之静态性能差。考虑到在光纤陀螺系统中,由于偏振串扰导致的相位误差仅仅是导致零偏的一部分,预估值小于实测值是合理的。
4 结论1) 光纤陀螺静态性能在温度范围的端点会有所降低;通过光路模型解得的预估值与通过温度实验得到的测量值,结果具有线性关系。本文中提到的光纤环琼斯矩阵模型不仅仅适用于光纤陀螺用光纤环,其他领域用光纤环同样适用。
2) 由于光纤环测量技术的限制,局限于在固定温度下测量光纤环,因此本文仅针对固定温度点下光纤陀螺零偏进行分析。缺少光纤陀螺在温度变化时,由于偏振串扰点位置和强度的变化对陀螺精度的影响分析。
3) 建立的光路模型具有准确性和适用性,在实际工程中可以作为一种预先筛选光纤环的手段,用以提高光纤陀螺精度。
4) 实验结论该模型具有准确性和通用性能够成为评估和优化光纤陀螺静态性能的理论依据。而要改进光纤环性能,则需选择保偏光纤、改进绕环技术等方向努力。
[1] |
VALI V, SHORTHILL R W. Fiber ring interferometer[J]. Applied optics, 1976, 15(5): 1099-1100. DOI:10.1364/AO.15.001099 (0)
|
[2] |
王新, 陈淑芬, 方伟. 干涉式光纤陀螺光路偏振特性的理论分析[J]. 光学技术, 2005, 31(2): 187-189. WANG Xin, CHEN Shufen, FANG Wei. Analysis of optical path polarization characteristic of interferometric fiber-optic gyroscope[J]. Optical technique, 2005, 31(2): 187-189. (0) |
[3] |
HAN Mei, SETHI A, HUA Wei, et al. A detection-based multiple object tracking method[C]//Proceedings of 2014 International Conference on Image Processing. Singapore, 2004: 3065-3068.
(0)
|
[4] |
李绪友, 张春梅, 刘华兵, 等. 光纤环十六极对称绕法温度性能的仿真与分析[J]. 中国惯性技术学报, 2016, 24(6): 780-785. LI Xuyou, ZHANG Chunmei, LIU Huabing, et al. Simulation and analysis on temperature performance of fiber ring by 16-polar symmetrical winding method[J]. Journal of Chinese inertial technology, 2016, 24(6): 780-785. (0) |
[5] |
李绪友, 张琛, 张勇. 单轴光纤陀螺光学结构的研究与实现[J]. 光电工程, 2009, 36(4): 52-55. LI Xuyou, ZHANG Chen, ZHANG Yong. Research and realization on optical system of single axis fiber optic gyroscope[J]. Opto-electronic engineering, 2009, 36(4): 52-55. (0) |
[6] |
周超超, 霍光, 卫炀, 等. 非理想集成光学器件对光纤陀螺的影响[J]. 激光杂志, 2016, 37(1): 44-48. ZHOU Chaochao, HUO Guang, WEI Yang, et al. Influence of Nonideal integrated optic circuit on FOG performance[J]. Laser journal, 2016, 37(1): 44-48. (0) |
[7] |
陈世同, 李玉深, 吴鹏. 光纤陀螺光路传输模型仿真[J]. 辽宁工程技术大学学报(自然科学版), 2010, 29(2): 267-270. CHEN Shitong, LI Yushen, WU Peng. Transmission model simulation of fiber optic gyroscope[J]. Journal of Liaoning Technical University (natural science), 2010, 29(2): 267-270. (0) |
[8] |
赵冰, 陈军波. 光纤陀螺温度建模与补偿方法分析[J]. 海军航空工程学院学报, 2016, 31(1): 12-16. ZHAO Bing, CHEN Junbo. Analysis on temperature modeling and compensation mothed of FOG[J]. Journal of naval aeronautical and Astronautical university, 2016, 31(1): 12-16. (0) |
[9] |
DING Zhenyang, MENG Zhuo, YAO X S, et al. Accurate method for measuring the thermal coefficient of group birefringence of polarization-maintaining fibers[J]. Optics letters, 2011, 36(11): 2173. DOI:10.1364/OL.36.002173 (0)
|
[10] |
PAVLATH G A, SHAW H J. Birefringence and polarization effects in fiber gyroscopes[J]. Applied optics, 1982, 21(10): 1752-1757. DOI:10.1364/AO.21.001752 (0)
|
[11] |
BU Lingbing, SA Rina, KIM D. Calibration method of polarization Lidar based on Jones matrix[J]. Optik, 2016, 130: 834-839. (0)
|
[12] |
杨军, 苑勇贵, 喻张俊, 等. 光纤偏振器件与组件的分布式串音测量研究进展[J]. 光电工程, 2018, 45(9): 170625. YANG Jun, YUAN Yonggui, YU Zhangjun, et al. Recent progress of accurate measurement for distributed polarization crosstalk of fiber optic polarization component and device[J]. Opto-Electronic engineering, 2018, 45(9): 170625. (0) |
[13] |
GJB 2426A-2015, 光纤陀螺仪测试方法[S].北京: 中国标准出版社, 2015.
(0)
|