留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

激光二极管端泵Yb:YAG晶体的温度场及应力场

李昕阳 李隆 任嘉欣 贺政隆 石镨 宁江浩 张春玲

肖凯博, 蒋新颖, 袁晓东, 等. 间隔掺杂低温Yb: YAG叠片放大器的热效应优化 [J]. 红外与激光工程, 2016, 45(12): 77-84. doi:  10.3788/IRLA20230683
引用本文: 肖凯博, 蒋新颖, 袁晓东, 等. 间隔掺杂低温Yb: YAG叠片放大器的热效应优化 [J]. 红外与激光工程, 2016, 45(12): 77-84. doi:  10.3788/IRLA20230683
肖凯博, 蒋新颖, 袁晓东, 等. 间隔掺杂低温Yb: YAG叠片放大器的热效应优化 [J]. 红外与激光工程, 2016, 45(12): 77-84. doi:  10.3788/IRLA20230683
Citation: 肖凯博, 蒋新颖, 袁晓东, 等. 间隔掺杂低温Yb: YAG叠片放大器的热效应优化 [J]. 红外与激光工程, 2016, 45(12): 77-84. doi:  10.3788/IRLA20230683

激光二极管端泵Yb:YAG晶体的温度场及应力场

doi: 10.3788/IRLA20230683
基金项目: 国家自然科学基金项目(11705136)
详细信息
    作者简介:

    李昕阳,男,硕士生,主要从事激光技术方面的研究

  • 中图分类号: TN248.1

Temperature field and stress field of LD end-pumped Yb:YAG crystal

Funds: National Natural Science Foundation of China (11705136)
  • 摘要: 为了解决激光二极管端面泵浦Yb:YAG晶体引起的热效应问题,通过对晶体工作特点的分析,建立了周边冷却恒温、端面与空气存在热交换的有限元热模型。利用泊松方程,对Yb:YAG晶体的温度场、热应力场、热形变场和热透镜焦距进行了数值计算,并定量分析了激光二极管泵浦光的高斯阶次、光斑半径和泵浦功率对激光晶体温度场的影响。研究结果表明:若激光二极管泵浦功率为50 W,耦合到泵浦面的光斑半径为400 μm时,晶体尺寸为3 mm×3 mm×4 mm、掺杂浓度为5.0 at.%的Yb:YAG晶体端面的最高温升为59.2 K,最大热形变量为0.64567 μm,晶体内稳定时最大应力为2.380×108 N/m2,热透镜焦距为19.99 mm,该条件下激光器可正常运行。研究结果为全固态Yb:YAG激光器的设计提供了理论依据。
  • 图  1  LD端面泵浦Yb:YAG晶体激光器实验装置简图

    Figure  1.  Schematic diagram of LD end-pumped Yb:YAG crystal laser experimental device

    图  2  LD端面泵浦Yb:YAG晶体几何结构图

    Figure  2.  LD end-pumped Yb:YAG crystal geometry diagram

    图  3  三维晶体(a)、泵浦面(b)及z轴截面(c)温度场分布图

    Figure  3.  Three-dimensional crystal (a), (b) pump surface and (c) z axis section temperature field distribution diagrams

    图  4  泵浦面(a)及z轴截面(b)等温线分布图

    Figure  4.  Pump surface (a) and z axis section (b) isotherm distribution diagrams

    图  5  不同高斯阶次下Yb:YAG晶体z轴上的温度场分布对比图

    Figure  5.  Comparison of temperature field distribution on z axis of Yb:YAG crystal at different Gauss order

    图  6  不同光斑半径的Yb:YAG晶体z轴上的温度场分布对比图

    Figure  6.  Comparison of temperature field distribution on z axis of Yb:YAG crystal with different spot radius

    图  7  不同泵浦功率下的Yb:YAG晶体z轴上的温度场分布对比图

    Figure  7.  Comparison of temperature field distribution on z axis of Yb:YAG crystal at different pump power

    图  8  形变开始时(a)、热应力稳定后(b)三维晶体应力分布图;形变开始时(c)、热应力稳定后(d)泵浦面应力分布图

    Figure  8.  Three-dimensional crystal stress distribution at the beginning of deformation (a) and after thermal stress stabilization; pump surface stress distribution at the beginning of deformation (c) and after thermal stress stabilization (d)

    图  9  Yb:YAG晶体泵浦面热形变分布图

    Figure  9.  Yb:YAG crystal pump surface thermal deformation distribution

    图  10  不同泵浦功率下热透镜焦距的变化

    Figure  10.  Thermal lens focal length varying with pump power

  • [1] 肖凯博, 蒋新颖, 袁晓东, 等. 间隔掺杂低温Yb: YAG叠片放大器的热效应优化 [J]. 红外与激光工程, 2016, 45(12): 77-84.

    Xiao Kaibo, Jiang Xinying, Yuan Xiaodong, et al. Optimization of thermal effects in a cryogenically cooled Yb: YAG multislab amplifier with interlayers [J]. Infrared and Laser Engineering, 2016, 45(12): 1206004. (in Chinese)
    [2] 唐靓, 叶慧琪, 肖东. YAG: Ce, Yb近红外下转换荧光粉中的离子簇集 [J]. 红外与激光工程, 2019, 48(5): 295-299.

    Tang Liang, Ye Huiqi, Xiao Dong. Ion clustering in YAG: Ce, Yb near-infrared down-conversion phosphors [J]. Infrared and Laser Engineering, 2019, 48(5): 0521004. (in Chinese)
    [3] 李国会, 杜应磊, 徐宏来, 等. 双变形镜对Yb: YAG板条激光器光束质量校正技术 [J]. 红外与激光工程, 2022, 51(8): 165-175.

    Li Guohui, Du Yinglei, Xu Honglai, et al. Correction of beam quality correction of Yb: YAG laser with double deformable mirrors [J]. Infrared and Laser Engineering, 2022, 51(8): 20210800. (in Chinese)
    [4] 耿鹰鸽, 李隆, 潘晓瑞, 等. 脉冲LD端面抽运变热导率方片Yb: YAG晶体温度场 [J]. 激光与光电子学进展, 2016, 53(8): 167-174.

    Geng Yingge, Li Long, Pan Xiaorui, et al. Temperature field of pulsed LD end-pumped square Yb: YAG crystal with variable thermal conductivity [J]. Laser & Optoelectronics Progress, 2016, 53(8): 081401. (in Chinese)
    [5] 张德平, 吴超, 张蓉竹, 等. LD端面泵浦分离式放大器结构的热效应研究 [J]. 红外与激光工程, 2015, 44(8): 2250-2255.

    Zhang Deping, Wu Chao, Zhang Rongzhu, et al. Study on thermal effect of LD end-pumped separated amplifier structure [J]. Infrared and Laser Engineering, 2015, 44(8): 2250-2255. (in Chinese)
    [6] 蒙裴贝, 颜凡江, 李旭, 等. 热边界和泵浦结构对激光晶体热效应的影响 [J]. 红外与激光工程, 2015, 44(11): 3216-3222.

    Meng Peibei, Yan Fanjiang, Li Xu, et al. Influence of boundary condition and pump scheme on thermal effects of laser crystal [J]. Infrared and Laser Engineering, 2015, 44(11): 3216-3222. (in Chinese)
    [7] 邱海齐, 张彬, 祝颂军. 反射式 LD 端面泵浦 Nd: YAG陶瓷激光器的温度场分析 [J]. 激光杂志, 2008, 29(2): 8-10.

    Qiu Haiqi, Zhang Bin, Zhu Songjun. Analysis of temperature field distribution in reflective LD end-pumped Nd: YAG ceramic laser [J]. Laser Journal, 2008, 29(2): 8-10. (in Chinese)
    [8] 田晓, 齐兵. 空气导热作用下Nd: YAG晶体温场特性 [J]. 中国光学, 2019, 12(3): 686-692. doi:  10.3788/co.20191203.0686

    Tian Xiao, Qi Bing. Temperature field of Nd: YAG crystal under air heat transfer [J]. Chinese Optics, 2019, 12(3): 686-692. (in Chinese) doi:  10.3788/co.20191203.0686
    [9] 李隆, 聂建萍, 史彭, 等. 激光二极管端面抽运YAG-Nd: YAG复合晶体棒温度场特性 [J]. 中国激光, 2010, 37(4): 917-922. doi:  10.3788/CJL20103704.0917

    Li Long, Nie Jianping, Shi Peng, et al. Temperature field characteristic of YAG-Nd: YAG composite crystal rod end-pumped by laser diode [J]. Chinese Journal of Lasers, 2010, 37(4): 917-922. (in Chinese) doi:  10.3788/CJL20103704.0917
    [10] 岱钦, 李新忠, 王希军. YAG物理特性参量对LD泵浦固体激光器热效应的影响 [J]. 光子学报, 2007, 36(12): 2201-2204.

    Dai Qin, Li Xinzhong, Wang Xijun. Influence of YAG materials property on thermal effects in LD side pumping solid state lasers [J]. Acta Photonica Sinica, 2007, 36(12): 2201-2204. (in Chinese)
    [11] 甘安生, 李隆, 史彭. 激光二极管端面泵浦Yb: YAG薄片激光器的热效应 [J]. 光子学报, 2008, 37(4): 631-635.

    Gan Ansheng, Li Long, Shi Peng. Thermal effect of Yb: YAG slice laser by diode laser end-pumped [J]. Acta Photonica Sinica, 2008, 37(4): 631-635. (in Chinese)
    [12] 邵铭, 张立文. Yb: YAG激光器特性 [J]. 激光杂志, 2003, 24(5): 15-17.

    Shao Ming, Zhang Liwen. Performance of Yb: YAG laser [J]. Laser Journal, 2003, 24(5): 15-17. (in Chinese)
    [13] 杨鸿儒, 左铁钏. 全固体小型Yb: YAG激光器热效应及输出特性研究 [J]. 光子学报, 2003, 32(8): 907-910.

    Yang Hongru, Zuo Tiechuan. Investigation on Thermal effects and output of all solid state miniature Yb: YAG laser [J]. Acta Photonica Sinica, 2003, 32(8): 907-910. (in Chinese)
    [14] 毛小洁, 钟国舜, 王旭, 等. LD端面泵浦掺Nd晶体的热透镜焦距 [J]. 激光与红外, 2010, 40(12): 1316-1318.

    Mao Xiaojie, Zhong Guoshun, Wang Xu, et al. Thermal lens focus of the end-pumped Nd-doped crystallaser [J]. Lasers & Infrared, 2010, 40(12): 1316-1318. (in Chinese)
    [15] 肖娜娜. 二极管泵浦Tm: YAP激光器的热透镜效应研究[D]. 哈尔滨工程大学, 2016.

    Xiao Nana. Research on thermal lens efffects in the LD pumped Tm: YAP laser[D]. Harbin: Harbin Engineering University, 2016. (in Chinese)
  • [1] 任嘉欣, 李隆, 李昕阳, 杨恒鑫, 纪玉潇, 张春玲.  激光二极管端泵方形Tm:YAG复合晶体的热效应 . 红外与激光工程, 2024, 53(3): 20230717-1-20230717-10. doi: 10.3788/IRLA20230717
    [2] 孙佳宁, 王雨雷, 张雨, 齐瑶瑶, 丁洁, 颜秉政, 白振旭, 吕志伟.  LD端面泵浦Er:Yb:glass/Co:MALO晶体热效应分析 . 红外与激光工程, 2023, 52(8): 20230349-1-20230349-12. doi: 10.3788/IRLA20230349
    [3] 蔡家轩, 万乐, 石世宏, 吴捷, 杨强, 石拓, 程梦颖.  激光熔化沉积AlSi10Mg温度场对显微组织性能的影响 . 红外与激光工程, 2022, 51(5): 20210366-1-20210366-10. doi: 10.3788/IRLA20210366
    [4] 张玥婷, 谭毅, 王继红, 彭起, 杨智焜.  激光窗口杂散光固气耦合热效应对光束质量的影响 . 红外与激光工程, 2022, 51(9): 20210966-1-20210966-9. doi: 10.3788/IRLA20210966
    [5] 张鹏辉, 赵扬, 李鹏, 周志权, 白雪, 马健.  基于有限元法的激光声磁检测系统优化研究 . 红外与激光工程, 2022, 51(7): 20210533-1-20210533-9. doi: 10.3788/IRLA20210533
    [6] 万颖超, 杨保来, 奚小明, 张汉伟, 叶云, 王小林.  不同泵浦波长光纤激光器模式不稳定效应对比 . 红外与激光工程, 2022, 51(4): 20210256-1-20210256-8. doi: 10.3788/IRLA20210256
    [7] 李志杰, 孔庆庆, 张明栋, 金子蘅, 卞殷旭, 沈华, 朱日宏.  万瓦级激光光闸的热效应分析及处理 . 红外与激光工程, 2022, 51(2): 20210909-1-20210909-8. doi: 10.3788/IRLA20210909
    [8] 刘全喜, 任钢, 李轶国, 岳通, 王莉, 肖星, 邓翠, 李佳玲.  激光二极管端面抽运梯度浓度掺杂介质激光器热效应的有限元法分析 . 红外与激光工程, 2019, 48(11): 1105004-1105004(11). doi: 10.3788/IRLA201948.1105004
    [9] 徐志钮, 胡志伟, 赵丽娟, 李永倩, 张翼, 胡世勋.  光电复合海缆温度场建模分析及在在线监测中的应用 . 红外与激光工程, 2018, 47(7): 717003-0717003(9). doi: 10.3788/IRLA201847.0717003
    [10] 孟令武, 邵帅, 乔健.  导光镜镜面温度场的模拟分析 . 红外与激光工程, 2018, 47(11): 1106008-1106008(7). doi: 10.3788/IRLA201847.1106008
    [11] 蒋志刚, 王浟, 蔡和, 高明, 韩聚洪, 安国斐, 张伟, 薛亮平, 王顺艳, 王宏元, 周杰.  具有三级放大结构的铷蒸气激光系统的优化设计 . 红外与激光工程, 2016, 45(11): 1106003-1106003(7). doi: 10.3788/IRLA201645.1106003
    [12] 张耀平, 樊峻棋, 龙国云.  变形镜在激光辐照下热畸变有限元模拟 . 红外与激光工程, 2016, 45(11): 1136002-1136002(5). doi: 10.3788/IRLA201645.1136002
    [13] 肖凯博, 蒋新颖, 袁晓东, 郑建刚, 郑万国.  间隔掺杂低温Yb:YAG叠片放大器的热效应优化 . 红外与激光工程, 2016, 45(12): 1206004-1206004(8). doi: 10.3788/IRLA201645.1206004
    [14] Sekou Singare, 陈盛贵, 钟欢欢.  激光透射焊接聚碳酸酯的有限元分析 . 红外与激光工程, 2016, 45(2): 206005-0206005(6). doi: 10.3788/IRLA201645.0206005
    [15] 张德平, 吴超, 张蓉竹, 孙年春.  LD 端面泵浦分离式放大器结构的热效应研究 . 红外与激光工程, 2015, 44(8): 2250-2255.
    [16] 肖冬明, 何宽芳, 王迪.  基于多层有限元模型的激光选区熔化多层瞬态温度场演化规律研究 . 红外与激光工程, 2015, 44(9): 2672-2678.
    [17] 许兆美, 汪通悦, 裴旭, 蒋素琴, 李伯奎, 王庆安, 洪宗海.  Al2O3陶瓷激光多道铣削温度场有限元模拟 . 红外与激光工程, 2015, 44(2): 477-481.
    [18] 许艳军, 赵宇宸, 沙巍, 齐光.  大尺寸SIC空间反射镜离子束加工热效应分析与抑制 . 红外与激光工程, 2014, 43(8): 2556-2561.
    [19] 张东来, 李小将, 杨业伟.  激光辐照飞行靶三维温度场数值模拟 . 红外与激光工程, 2014, 43(9): 2883-2888.
    [20] 马依拉木·木斯得克, 姚建铨, 王鹏.  LD端面泵浦946 nm/473 nm连续Nd:YAG/LBO激光器 . 红外与激光工程, 2013, 42(11): 2931-2934.
  • 加载中
图(10)
计量
  • 文章访问数:  51
  • HTML全文浏览量:  10
  • PDF下载量:  18
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-12-08
  • 修回日期:  2024-01-15
  • 刊出日期:  2024-03-21

激光二极管端泵Yb:YAG晶体的温度场及应力场

doi: 10.3788/IRLA20230683
    作者简介:

    李昕阳,男,硕士生,主要从事激光技术方面的研究

基金项目:  国家自然科学基金项目(11705136)
  • 中图分类号: TN248.1

摘要: 为了解决激光二极管端面泵浦Yb:YAG晶体引起的热效应问题,通过对晶体工作特点的分析,建立了周边冷却恒温、端面与空气存在热交换的有限元热模型。利用泊松方程,对Yb:YAG晶体的温度场、热应力场、热形变场和热透镜焦距进行了数值计算,并定量分析了激光二极管泵浦光的高斯阶次、光斑半径和泵浦功率对激光晶体温度场的影响。研究结果表明:若激光二极管泵浦功率为50 W,耦合到泵浦面的光斑半径为400 μm时,晶体尺寸为3 mm×3 mm×4 mm、掺杂浓度为5.0 at.%的Yb:YAG晶体端面的最高温升为59.2 K,最大热形变量为0.64567 μm,晶体内稳定时最大应力为2.380×108 N/m2,热透镜焦距为19.99 mm,该条件下激光器可正常运行。研究结果为全固态Yb:YAG激光器的设计提供了理论依据。

English Abstract

    • 在工业、医疗、国防等领域有着广泛应用的激光二极管(LD)端面泵浦Yb:YAG激光器,因其有着峰值功率高、寿命长等特点,成为了激光技术研究者们的关注对象[13]。与Nd:YAG晶体相比,Yb:YAG晶体在高掺杂情况下不会发生浓度猝灭现象,而且量子亏损很小,更容易产生高效高峰值功率的脉冲激光[4]。然而,激光二极管泵浦Yb:YAG晶体在产生激光辐射的同时也会产生大量的热。热的沉积会导致晶体内温度分布不均匀,使激光晶体产生热效应,严重影响了激光的输出功率及光束质量,减小热效应已经成为激光器设计中的关键步骤[57]

      以往的热效应研究常采用解析法计算端面绝热晶体热模型的温度场分布[810]。然而,这种模型忽略了空气与端面的热交换,分析结果存在一定误差。文中建立了晶体周边恒温冷却、通光面与空气热交换的热模型,采用有限元分析法,对LD端面泵浦的Yb:YAG方形晶体晶体内的温度场、应力场、形变场及热焦距进行了计算,并对不同高斯阶次、光斑半径和泵浦功率下的晶体温度场分布变化进行了定量分析。

    • LD发出的泵浦光经耦合镜聚焦后射入激光晶体端面中心。图1为LD端面泵浦Yb:YAG晶体激光器实验装置简图,图2为LD端面泵浦Yb:YAG晶体几何结构图。其中,晶体尺寸为a×b×c,晶体的通关面记作z=0、z=c,令泵浦光线从z=0的端面射入。在晶体端面镀有940 nm波长的增透膜(透射率大于99.8%)和1030 nm波长的增透膜(透射率大于99%)。作为谐振腔输入镜M1,镀有940 nm波长的增透膜(透射率大于99.8%),同时镀有1030 nm波长的高反膜(反射率大于98%);作为谐振腔输出镜M2,镀有1030 nm波长的高反膜(反射率为91%,透射率为9%)。此外,Yb:YAG晶体应涂有导热硅脂并被铟膜包裹,也要尽可能完全接触热沉,达到快速传热冷却的目的。

      图  1  LD端面泵浦Yb:YAG晶体激光器实验装置简图

      Figure 1.  Schematic diagram of LD end-pumped Yb:YAG crystal laser experimental device

      图  2  LD端面泵浦Yb:YAG晶体几何结构图

      Figure 2.  LD end-pumped Yb:YAG crystal geometry diagram

      泵浦光入射到单个光纤中时,其光纤出射光的功率分布为高斯分布,而激光器中由一捆光纤发射到光学耦合元件上的光可视为高斯光的叠加,其光强分布呈“大礼帽”状,称为超高斯分布[11],光强表达式为:

      $$ I\left( {x,y,z} \right) = {I_0}{\text{exp}}\left( {{{ - 2}}\frac{{{x^{2k}}{\text{ + }}{y^{2k}}}}{{\omega _0^{2k}}}} \right){\text{exp}}\left( { - \beta {{z}}} \right) $$ (1)

      式中:I0为泵浦功率;ω0为泵浦光高斯半径;k为高斯阶次;β为激光介质对泵浦光的吸收系数。

      晶体内部的荧光量子效应和内损耗对泵浦光的消耗远大于其他因素造成的光损耗,因此仅考虑量子效应和内损耗所带来的热效应。则z轴截面的热功率密度为:

      $$ {q_v}\left( {x,y,z} \right) = \beta \eta {I_0}\exp \left( { - 2\frac{{{x^{2k}} + {y^{2k}}}}{{\omega _0^{2k}}})} \right)\exp \left( {{{ - }}\beta z} \right) $$ (2)

      式中:η为热转换系数,其数值由量子效应及内损耗决定。有:

      $$ \eta = 1 - {\lambda _{\rm{p}}}/{\lambda _{\rm{L}}} $$ (3)

      式中:λp为泵浦光波长,940 nm;λL为出射激光波长,1030 nm。

      激光器工作一段时间后,其工作状态和热力学状态达到稳定,此时激光器输出功率稳定,Yb:YAG晶体内部温度场呈梯度分布。方形晶体Yb:YAG的热传导方程可以表示为:

      $$ \frac{{{\partial ^2}T\left( {x,y,z} \right)}}{{\partial {x^2}}}{\text{ + }}\frac{{{\partial ^2}T\left( {x,y,z} \right)}}{{\partial {y^2}}}{\text{ + }}\frac{{{\partial ^2}T\left( {x,y,z} \right)}}{{\partial {z^2}}} = - \frac{{{q_v}}}{\lambda } $$ (4)

      式中:T为温度函数;λ为晶体导热系数。

      涂有导热硅脂并被铟膜包裹的晶体与冷却铜块紧密接触,晶体侧面受到铜块限制,仅有与空气接触的端面可自由形变,该情况下z轴方向的应力起到引起形变的主导作用,故将晶体视为各向同性的线弹性材料。

    • 通常激光器稳定工作过程中会使用循环水冷,减小热能对激光工作的影响,因此初始温度设为冷却水温度:

      $$ T = {T_0} $$ (5)

      考虑到晶体周边恒温冷却、端面与空气热交换,提出以下边界条件。公式(6)描述的是周边恒温冷却条件,x=−a/2、x=a/2、y=−b/2、y=b/2,这四个面始终为冷却水温度。公式(7)、公式(8)为第三类边界条件,z=0、z=c,两个与空气相接触的端面始终与空气热交换。

      $$ \left\{\begin{array}{l} T(-a / 2, y, z)=T_0 ; T(a / 2, y, z)=T_0 \\ T(x,-b / 2, z)=T_0 ; T(x, b / 2, z)=T_0 \end{array}\right. $$ (6)
      $$ \lambda \frac{{\partial T}}{{\partial z}}{|_{z = 0}} = h(T{|_{z = 0}} - {T_0}) $$ (7)
      $$ - \lambda \frac{{\partial T}}{{\partial z}}{|_{z = c}} = h(T{|_{z = c}} - {T_0}) $$ (8)

      式中:h为空气传热系数。

      此外,激光器水冷装置中的铜块与涂有导热硅脂并被铟膜包裹的晶体紧密接触,故将晶体侧面视为受到铜块限制的固定约束;而端面不受约束,仅与空气接触,故视为可自由形变的端面。

    • 按照上述理论建立热模型,该模型尺寸为3 mm×3 mm×4 mm,掺杂浓度为5.0 at.%,对于940 nm泵浦光的吸收系数为5.6 cm−1[1113],导热系数为0.13 W·cm−1·K−1,高斯阶次为1,冷却温度T0为291 K,热膨胀系数为7.8×10−6 K−1

      当激光泵浦功率为50 W、泵浦面光斑半径为400 μm时,系统泵浦端面最高温升为59.2 K。图3为激光二极管端面泵浦的Yb:YAG晶体温度场分布图,图4为激光二极管端面泵浦的Yb:YAG晶体内部等温线分布图。

      图  3  三维晶体(a)、泵浦面(b)及z轴截面(c)温度场分布图

      Figure 3.  Three-dimensional crystal (a), (b) pump surface and (c) z axis section temperature field distribution diagrams

      图  4  泵浦面(a)及z轴截面(b)等温线分布图

      Figure 4.  Pump surface (a) and z axis section (b) isotherm distribution diagrams

    • 考虑到可以通过改变光纤疏密来调制光斑分布函数,期望通过对比不同高斯阶次的功率函数对激光晶体温度场的影响,获取更为合适的功率函数,对激光器设计提供依据。在其他条件与图3一致的情况下,高斯阶次分别选择1、2、3,其对应的泵浦端面最高温升分别为59.18、75.24、83.99 K。图5为不同高斯阶次下Yb:YAG晶体z轴上的温度场分布对比图,随着高斯阶次提升,温升也在提高,这是因为高斯阶次的提升会使光能的分布更加集中。

      图  5  不同高斯阶次下Yb:YAG晶体z轴上的温度场分布对比图

      Figure 5.  Comparison of temperature field distribution on z axis of Yb:YAG crystal at different Gauss order

    • 泵浦光与激光晶体内振荡模的匹配是提高Yb:YAG激光器输出光束质量的关键[14]。在其他条件与图3一致的情况下,只改变泵浦光斑半径,分别为300、350、400、450、500 μm时,其对应的泵浦端面最高温升为105.21、77.30、59.18、46.76、37.88 K。图6为不同光斑半径的Yb:YAG晶体z轴上的温度场分布对比图,可以看出,光斑半径越小,温升越大,这是因为光斑半径越小,其功率密度越大。

      图  6  不同光斑半径的Yb:YAG晶体z轴上的温度场分布对比图

      Figure 6.  Comparison of temperature field distribution on z axis of Yb:YAG crystal with different spot radius

    • 在符合实验要求的前提下选择适当的泵浦功率,能有效降低温度场的温度。在其他条件与图3一致的情况下,只改变泵浦功率,分别为50、60、70 W,其对应的最高温升为59.18、71.02、82.85 K。图7为不同泵浦功率下的Yb:YAG晶体z轴上的温度场分布对比图,可以看出,功率越大,温升越大。

      图  7  不同泵浦功率下的Yb:YAG晶体z轴上的温度场分布对比图

      Figure 7.  Comparison of temperature field distribution on z axis of Yb:YAG crystal at different pump power

    • Yb:YAG晶体在工作过程中不断产生热量,积累的热量引起温度场分布的变化,并使晶体内应力场与热形变场发生变化,在与图3条件一致的情况下,系统所受最大应力为2.380×108 N/m2,产生的最大形变为6.4567×10−4 mm。图8为Yb:YAG晶体应力场分布图,图9为Yb:YAG晶体热形变分布图。

      图  8  形变开始时(a)、热应力稳定后(b)三维晶体应力分布图;形变开始时(c)、热应力稳定后(d)泵浦面应力分布图

      Figure 8.  Three-dimensional crystal stress distribution at the beginning of deformation (a) and after thermal stress stabilization; pump surface stress distribution at the beginning of deformation (c) and after thermal stress stabilization (d)

      图  9  Yb:YAG晶体泵浦面热形变分布图

      Figure 9.  Yb:YAG crystal pump surface thermal deformation distribution

    • 产生激光晶体热透镜效应的原因主要是大量的泵浦能量以热的形式沉积在激光晶体中,而晶体在工作状态下产生热量的主要原因有:能级间非辐射跃迁所释放的能量被基质材料吸收;泵浦波长范围内的部分能量因与晶体泵浦能带不匹配而被基质材料吸收;晶体的量子亏损与内损耗。这些因素使晶体内热量沉积,并形成温度梯度、应力梯度及形变。因此,在计算LD端面泵浦Yb:YAG方形晶体的热透镜焦距时,考虑到的主要因素为温度梯度、端面形变和热致双折射[15]

      温度梯度、热致双折射、端面形变造成了热透镜效应,它们产生的热透镜焦距分别为:

      $$ {f_{TO}} = \frac{{\pi {K_c}\omega _p^2}}{{{P_{ph}}{\mathrm{d}}n/{\mathrm{d}}T}}\left[ {\frac{1}{{1 - \exp ( - \alpha l)}}} \right] $$ (9)
      $$ {f_{EO}} = \frac{{\pi {K_c}\omega _p^2}}{{2{\alpha _T}{n^3}{C_{\varphi ,r}}{P_{ph}}[1 - \exp ( - \alpha l)]}} $$ (10)
      $$ {f}_{TE}=\frac{\pi {K}_{c}{\omega }_{p}^{2}}{{\alpha }_{T}{r}_{o}{P}_{ph}[1-\mathrm{exp}(-\alpha {r}_{0})]2(n-1)} $$ (11)

      式中:Kc为晶体热导率;dn/dT为热光系数;Cφ,r为弹性系数;αT为热膨胀系数;α为吸收系数;Pph为泵浦光功率的致热功率,其与泵浦功率关系如下。

      $$ {P_{ph}} = {P_{in}}\xi $$ (12)

      式中:ξ为热转化效率。由公式(10)~(12)可得综合热透镜焦距为:

      $$ {f_{th}} = \frac{1}{{{1}/{{{f_{TO}}}} + {1}/{{{f_{EO}}}} + {1}/{{{f_{TE}}}}}} $$ (13)

      相对于温度梯度带来的热透镜效应,热应力双折射与热形变产生的热透镜效应十分微弱[15],因此,只关注温度分布时热透镜焦距随功率的变化,如图10所示,公式(14)为简化后只考虑温度梯度的热透镜效应。

      $$ {f_{th}} = \frac{{{\text{π }}{K_c}\omega _p^2}}{{{P_{in}}\xi {\mathrm{d}}n/{\mathrm{d}}T}}\left[ {\frac{1}{{1 - \exp ( - \alpha l)}}} \right] $$ (14)

      图  10  不同泵浦功率下热透镜焦距的变化

      Figure 10.  Thermal lens focal length varying with pump power

    • 考虑到晶体周边恒温冷却、通光面与空气热交换,建立了更符合激光二极管端面泵浦Yb:YAG晶体实际工作状态的热模型,利用泊松方程和有限元分析法,对激光二极管端面泵浦Yb:YAG方形晶体的温度场分布进行了分析,并进一步计算了Yb:YAG晶体内由非均匀温升引起的的应力场以及热形变场的变化情况,最后得出一定条件下热透镜焦距随功率变化的曲线图。当激光二极管泵浦功率为50 W、泵浦光高斯半径为400 μm时,晶体尺寸为3 mm×3 mm×4 mm、掺杂浓度为5.0 at.%的Yb:YAG晶体端面的最高温升为59.2 K,最大热形变量为0.64567 μm,晶体内稳定时最大应力为2.426×107 N/m2,热透镜焦距为19.99 mm。此外,最大温升会随着高斯阶次和功率的增大、光斑半径的减小而提高,这些量的变化实际上是功率密度的变化,故可得出最大温升与功率密度正相关的结论。与文献[4]中的实验测量结果进行比较,理论计算结果基本符合,所得结果为全固态Yb:YAG激光器的设计提供了理论依据。

参考文献 (15)

目录

    /

    返回文章
    返回