基于Lamb理论的激光陀螺光强调谐线型研究 下载: 1095次
1 引言
激光陀螺是捷联式惯性制导与导航系统的理想器件,具有启动时间短、动态范围大、标度因数稳定等优点,在导弹制导、航空器飞行控制以及航天遥感卫星姿态控制等领域占有重要的地位[1-2]。环形激光器是激光陀螺的核心[3-4],为了保持窄频带、低功耗、高稳定工作,激光陀螺环形激光器一般采用He-Ne气体作为增益介质[5]。研究表明,当环形激光器工作纵模满足一定的条件时,激光陀螺呈现自偏频状态[6-7]。
研究调制环形气体激光器光强调谐曲线线型的方法,为实现这种激光器强弱双纵模满足特定的振荡强度比例关系(如:Istrong∶
本文采用Lamb理论和等离子体色散函数,研究增益介质Ne同位素频率分裂、同位素配比、充气压力等因素对环形激光器光强调谐曲线线型的影响。建立环形激光器光强调谐曲线物理模型,对Ne20和Ne22充气比例分别为1∶1和7∶3、充气压力为400Pa条件下,环形激光器的光强调谐曲线进行仿真分析。实验上搭建了激光陀螺光强调谐曲线调制实验平台,验证了理论分析的正确性。最后,在上述理论基础上,研究了Ne双同位素的组分、配比等因素对环形激光器光强调谐曲线线型的影响,该研究为双纵模自偏频激光陀螺稳频工作建立了基础。
2 环形激光器
环形激光器是激光陀螺的核心器件,其具体结构如
环形激光器由四块反射镜构成环形光路。其中,两块反射镜为平面镜,另两块反射镜为凹面镜[14]。该环形激光器以He-Ne混合气体作为增益介质。按照生产工艺要求,为保证腔内He-Ne增益介质纯净、无杂气,环形激光器充入增益气体需要逐次经历第一次充气、老练、真空排气、第二次充气等工序。第一次充入腔体的气体与第二次充入气体的同位素配比完全相同。第一次充气后,环形激光器在高电压激励下连续老练,一周后排除腔内气体,并进行真空处理。然后,再次充入纯净的He-Ne混合气体。这样做的目的是充分排除微晶玻璃腔体分子结构中可能包裹的杂气粒子,避免陀螺长期工作时因增益气体污染而造成失效。此外,在平面镜对侧沉孔内放置吸气剂,工作时实时吸附激光器产生的杂气粒子。上述措施的目的是最大限度地保证激光陀螺增益介质的纯净性,避免增益介质中混入杂质引起激光器输出光强调谐曲线线型在激光陀螺额定寿命内出现畸变,进而影响陀螺的稳频精度。增益介质中逐渐混入杂气,将导致环形激光器阈值升高、限模结构变化、稳频精度降低,甚至造成激光器“光发白”等失效性质的故障。接下来,针对环形激光器光强调谐曲线线型与增益介质气体配比、充气压力等因素的关系进行研究。
3 Ne同位素与环形激光器光强调谐曲线线型
3.1 Lamb理论与环形谐振腔
激光Lamb理论利用麦克斯韦方程组描述光频电磁场,量子力学描述物质粒子,进而分析光和物质的相互作用[15-16]。该理论可以较好地解释激光频率间的耦合效应对激光增益的影响,适合研究Ne同位素频率分裂以及频率耦合造成的激光光强调谐曲线线型变化等问题。
根据激光Lamb理论,含有增益介质激励源项与非均匀损耗项的环形谐振腔行波场方程为
式中:σ(r)为阻尼项,唯象地表征激光场分布在环形谐振腔内受到各种损耗所造成的衰减;E(r,t)和P(r,t)分别表示激活增益介质的电场强度和电极化强度,它们是空间坐标r、时间坐标t的函数;c为光速;
式中:PI(t)和PII(t)分别为电极化强度P(r,t)的空间傅里叶分量;φ为电场强度的相位;k为波矢,k1和k2分别表示正反行波的波矢;ω1和ω2为与转动方向同、反向的角频率。
将(2)式代入(1)式,保留一阶项,如
式中:α为单程净增益系数;β为自饱和系数;θ为互饱和系数;σ为模牵引系数;ρ为模自推斥系数;τ为模互饱和系数;Ω0为理想空腔Sagnac效应频率;ψ为行波相位差,
当环形激光器工作在稳定条件下
令
将(7)、(8)两式代入(9)式,得
令
将(11)式代入(10)式中,并略去二阶小量,得到准稳态条件下,单纵模工作环形激光器稳定工作时的平均光强为
式中:R为辐射捕获修正量,它的取值与增益介质的充气压力相关。(12)式即为激光器输出光强度随环形谐振腔振荡频率变化的曲线I(ν)-ν,即激光器光强调谐曲线。该曲线是激光器稳频的依据,曲线的形貌以及斜率直接影响稳频精度。接下来,借助等离子体色散函数,继续求解He-Ne激光器光强调谐曲线。
3.2 等离子体色散函数理论
等离子体色散函数Z(ξ,η)通过热平衡状态下的麦克斯韦-玻尔兹曼分布,描述增益介质粒子(经典粒子)的等离子体状态[17]。环形激光器工作在稳定状态下,利用Lamb理论求解环形激光器中与转动方向同、反向的行波间,受频率耦合效应影响的各自光强随时间变化的特性时,需要求解(12)式中对应的各阶Lamb系数。而气体环形激光器的各阶Lamb系数与频率耦合效应、原子热运动所产生的多普勒加宽、气体均匀展宽等因素密切相关,所以Lamb系数可以用等离子体色散函数表示。等离子体色散函数的表达式可以写为
式中:ξ为频率参量;η为均匀展宽和非均匀展宽之比;L为洛伦兹函数。
等离子体色散函数是复函数,(13)式的实部、虚部分别表示为
由(14)、(15)式可知,等离子体色散函数为奇异积分函数,因此,不存在精确的显函数表达式。但可以通过数值积分,求出不同激光振荡频率对应的等离子体色散函数的值。利用Matlab软件,对等离子体色散函数Z(ξ,η)的实部与虚部以及它们的一阶导数进行数值计算,其计算结果如
图 2. 等离子体色散函数数值计算结果。(a)色散函数实部;(b)色散函数虚部;(c)色散函数实部的一阶导数;(d)色散函数虚部的一阶导数
Fig. 2. Numerical results of plasma dispersion function. (a) Real part of dispersion function; (b) imaginary part of dispersion function; (c) first derivative of dispersion function real part; (d) first derivative of dispersion function imaginary part
3.3 Ne双同位素与光强调谐曲线线型
Ne20和Ne22同位素各自的增益中心频率位置存在261MHz的频率分裂[18]。自然界中存在的Ne元素同位素成分复杂、各同位素配比存在差异。由于Ne同位素原子质量差异等因素,将导致He-Ne气体增益介质的激活发光特性存在一定的区别。假设Ne20和Ne22的多普勒展宽分别为ku、
按照等离子体色散函数,(3)式和(4)式中的Lamb系数α、β、θ可以表示为
式中:假设增益介质为Ne20、Ne22、He4混合气体;F为同位素中Ne20所占比例,
(19)式为单纵模工作环形激光器光强调谐曲线的表达式。其中,除等离子体色散函数外,其他参数均可实验测得。等离子体色散函数是ξ、η的函数,因此,光强调谐曲线仅是ξ、η的函数,代入等离子体色散函数可以数值求解。
将色散函数及其一阶导数的表达式代入(19)式,求解获得光强调谐曲线和ξ、η的关系,其中ξ与双同位素组分、配比相关;η与充气压力相关。单纵模工作激光陀螺光强调谐曲线仿真结果如
图 3. 环形激光器光强调谐曲线仿真结果。(a) Ne20∶Ne22=1∶1,充气压力400 Pa,光强调谐曲线;(b) Ne20∶Ne22=7∶3,充气压力400 Pa,光强调谐曲线
Fig. 3. Simulation results of intensity tuning curves in ring laser. (a) Intensity tuning curves under Ne20∶Ne22=1∶1, 400 Pa pressure; (b) intensity tuning curves of Ne20∶Ne22=7∶3, 400 Pa pressure
由
图 4. 环形激光器光强调谐曲线仿真结果。(a) Ne20∶Ne27=1∶1,充气压力400 Pa,光强调谐曲线;(b) Ne20∶Ne27=7∶3,充气压力400 Pa,光强调谐曲线
Fig. 4. Simulation results of intensity tuning curves in ring laser. (a) Intensity tuning curves of Ne20∶Ne27=1∶1, 400 Pa pressure; (b) intensity tuning curve of Ne20∶Ne27=7∶3, 400 Pa pressure
由
4 实验结果
搭建激光陀螺光强调谐曲线线型检测系统如
图 5. 激光陀螺光强调谐曲线线型检测系统
Fig. 5. Linear type detection system of laser gyro intensity tuning curve
法布里-珀罗(F-P)腔扫描获得环形激光器纵模状态如
图 7. 激光陀螺光强调谐曲线检测方法示意图。(a)纵模扫描及光强调谐曲线获得;(b)扫模电压
Fig. 7. Schematic of detection method of intensity tuning curve in laser gyro. (a) Longitudinal mode is swept and intensity tuning curve is obtained; (b) voltage of mode sweeping
选取同一个待测环形谐振腔,保持增益介质充气压力400 Pa不变,
图 8. 单纵模工作激光陀螺光强调谐曲线测试结果。(a) Ne20∶Ne22=7∶3测试结果;(b)自然Ne同位素混合气体测试结果
Fig. 8. Test results of intensity tuning curve of single longitudinal mode operated laser gyro. (a) Test result of Ne20∶Ne22=7∶3; (b) test result of natural Ne isotope mixture gas
5 结论
激光陀螺输出的光强调谐曲线线型与气体增益介质的同位素组分、配比、充气压力密切相关。光强调谐曲线的线型直接影响着陀螺的稳频精度。本文采用Lamb理论、等离子体色散理论研究了He-Ne增益介质中Ne元素同位素中心频率分裂以及频率耦合对光强调谐曲线的影响。实验建立了激光陀螺输出光强调谐曲线检测实验系统,采用天然Ne同位素混合气体、Ne20和Ne22双同位素按照7∶3的比例混合,以400 Pa气压充入Ne-Ne混合气体,实验获得了陀螺输出光强调谐曲线,实验结果验证了理论分析的正确性。该研究为实现一种环形激光器双纵模光强调谐曲线频率非对称位置稳频提供理论基础。
[1] Hokmabadi M P, Schumer A, Christodoulides D N, et al. Non-Hermitian ring laser gyroscopes with enhanced Sagnac sensitivity[J]. Nature, 2019, 576(7785): 70-74.
[2] Hinuber E L V, Reimer C, Schneider T, et al. INS/GNSS integration for aerobatic flight applications and aircraft motion surveying[J]. Sensors, 2017, 17(5): 941.
[3] Chow W W, Gea-Banacloche J, Pedrotti L M, et al. The ring laser gyro[J]. Reviews of Modern Physics, 1985, 57(1): 61-104.
[4] Nejad SM, PourmahyabadiM. The performance modeling of ring laser gyro in inertial navigation[C]//2006 13th IEEE International Conference on Electronics, Circuits and Systems. December 10-13, 2006, Nice, France.New York: IEEE Press, 2006: 318- 321.
[5] 姜亚南. 环形激光陀螺[M]. 北京: 清华大学出版社, 1985: 214- 216.
Jiang YN. Ring laser gyro[M]. Beijing: Tsinghua University Press, 1985: 214- 216.
[6] Liu J N, Jiao M X, Jiang J B, et al. Self-biasing phenomenon in prism laser gyro operating in double-longitudinal-mode state[J]. Optics Express, 2018, 26(24): 32353-32364.
[7] 刘健宁, 焦明星, 马家君, 等. 棱镜式激光陀螺双纵模自偏频现象研究[J]. 兵工学报, 2017, 38(6): 1113-1119.
Liu J N, Jiao M X, Ma J J, et al. Study of the double longitudinal mode operation and the self-biasing phenomenon in prism laser gyro[J]. Acta Armamentarii, 2017, 38(6): 1113-1119.
[8] Hurst R B, Mayerbacher M, Gebauer A, et al. High-accuracy absolute rotation rate measurements with a large ring laser gyro: establishing the scale factor[J]. Applied Optics, 2017, 56(4): 1124-1130.
[9] 任莉娜, 刘健宁, 焦明星. 温度扰动对棱镜式激光陀螺特性的影响[J]. 光子学报, 2018, 47(6): 0614002.
[10] 马家君, 蒋军彪, 刘健宁. 全反射棱镜式激光陀螺双纵模稳频技术[J]. 中国激光, 2014, 41(9): 0902011.
[11] 贾惠霞, 张英敏, 张琼, 等. 一种棱镜式激光陀螺光强控制系统研究[J]. 光学学报, 2014, 34(8): 0814001.
[12] 金世龙, 龙兴武, 王飞, 等. 激光陀螺克服高、低温冲击的技术研究[J]. 光学学报, 2006, 26(3): 409-414.
[13] 蒋军彪, 马家君, 刘建宁. 折射率非均匀性对激光陀螺背向散射的影响分析[J]. 激光与光电子学进展, 2016, 53(6): 061402.
[14] 韩宗虎, 胡晓东. 激光陀螺反射镜散射检测方法[J]. 中国惯性技术学报, 2015, 23(4): 540-543.
Han Z H, Hu X D. Method of mirror scattering test for laser gyro[J]. Journal of Chinese Inertial Technology, 2015, 23(4): 540-543.
[15] Izmailov E A, Kukhtievitch S E, Tikhomirov V V, et al. Analysis of laser gyro drift components[J]. Gyroscopy and Navigation, 2016, 7(1): 24-28.
[16] 周炳琨, 高以智, 陈倜嵘, 等. 激光原理[M]. 6版.北京: 国防工业出版社, 2009: 278- 281.
Zhou BK, Gao ZY, Chen TR, et al.Principles of laser[M].6th ed. Beijing: National Defense Industry Press, 2009: 278- 281.
[17] 徐茵. 等离子体色散函数的计算方法[J]. 计算物理, 1989, 6(1): 9-18.
[18] Aronowitz F. Single-isotope laser gyro[J]. Applied Optics, 1972, 11(2): 405-412.
[19] Audi G, Bersillon O, Blachot J, et al. The Nubase evaluation of nuclear and decay properties[J]. Nuclear Physics A, 2003, 729(1): 3-128.
Article Outline
刘健宁, 刘玉洁, 焦明星, 蒋军彪, 翁浚, 郭兰朋. 基于Lamb理论的激光陀螺光强调谐线型研究[J]. 光学学报, 2020, 40(21): 2114001. Jianning Liu, Yujie Liu, Mingxing Jiao, Junbiao Jiang, Jun Weng, Lanpeng Guo. Study on the Intensity Tuning Curve in a Laser Gyro Based on Lamb Theory[J]. Acta Optica Sinica, 2020, 40(21): 2114001.