高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

连续激光辐照三结GaAs太阳电池热应力场研究

谭宇 陆健

引用本文:
Citation:

连续激光辐照三结GaAs太阳电池热应力场研究

    作者简介: 谭宇(1993-), 男, 硕士研究生, 主要从事激光与物质相互作用方面的研究.
    通讯作者: 陆健, lujian@njust.edu.cn
  • 基金项目:

    上海航天科技创新基金资助项目 SAST20161113

  • 中图分类号: TN249

Study on three-junction GaAs solar cell thermal stress field by continuous wave laser irradiation

    Corresponding author: LU Jian, lujian@njust.edu.cn
  • CLC number: TN249

  • 摘要: 为了研究1070nm光纤连续激光对三结GaAs太阳电池的热力作用, 使用COMSOL软件建立物理模型进行了数值仿真, 得到了不同激光功率密度激光作用下的热应力场。为了验证了热应力计算方法的正确性, 利用栅线投影法测量了功率为16.7W、辐照半径为1mm、辐照时间为10s的激光作用下电池表面的形变, 实验结果与模拟结果基本吻合。结果表明, 辐照半径为1.5cm、功率密度为16.7W/cm2的激光辐照20s时, 底电池中心温度刚好超过电池使用温度; 底电池中心等效应力为96.6MPa, 刚好超过底电池材料的屈服极限; 根据这一结果可推测电池失效与热应力导致的结构损伤有关。该数值模拟结果与实验结果为激光辐照太阳电池的热力效应研究提供了一定的理论依据。
  • Figure 1.  Model of solar cell

    Figure 2.  Comparison curve between experimental and calculated values of battery center temperature under 5W/cm2 power density laser irradiation

    Figure 3.  Equivalent stress curve of center point under laser irradiation of different power densities

    Figure 4.  Temperature curve of center point under laser irradiation of different power densities

    Figure 5.  Displacement curve of center point of battery surface under laser irradiation with power of 16.7W, radius of 1mm and irradiation time of 10s

    Figure 6.  Beam path of solar cells under laser irradiation

    Figure 7.  Battery surface image before irradiation

    Figure 8.  Battery surface image before irradiation

    Figure 9.  Deformation diagram of battery central area under irradiation

    Table 1.  Forbidden band width and absorbable wavelength range of each layer of battery

    material forbidden bandwidth/eV absorption wavelength range/nm
    GaInP2 1.85 350~700
    GaAs 1.42 700~880
    Ge 0.67 880~1750
    下载: 导出CSV

    Table 2.  Thermal parameters of Ge

    parameter value
    density ρ/(g·cm-3) 5.49
    thermal conductivity κ/(W·cm·K-1) 0.59
    heat capacity C/(J·K-1) 0.31
    thermal radiation rate φ 0.1
    Stefan constant σ/(W·cm-2·K-4) 5.67×10-14
    coefficient of heat transfer h/(W·cm-2·K-1) 1×10-5
    下载: 导出CSV

    Table 3.  Mechanical parameters of Ge

    parameter value
    Yang’s module E/(N·cm-2) 1.31×107
    fracture module f/(N·cm-2) 9.3×103
    thermal expansion coefficient αt/(K-1) 5.92×106
    Poisson ration μ 0.21
    下载: 导出CSV
  • [1]

    ZHANG Zh W, LU J F, CHI W Y, et al. Techbique development and prospects analysis of GaAs solar cell [J].Aerospace Shanghai, 2003, 20(3):33-38(in Chinese). 
    [2]

    XIANG X B, DU W H, CHANG X L, et al. The study on high efficient GaAs/Ge solar cells [J]. Solar Energy Materials & Solar Cells, 2001, 68(1):97-103.
    [3]

    LANDIS G A. Space power by ground-based laser illumination[J]. IEEE Aerospace & Electronic Systems Magazine, 1991, 6(11):3-7. 
    [4]

    STEINSIEK F. Wireless power transmission experiment as an early contribution to planetary exploration missions[C]//54th International Astronautical Congress of the International Astronautical Federation. New York, USA: the International Academy of Astronautics and the International Institute of Space Law, 2003: R306.
    [5]

    MEYER J R, KRUER M R, BARTOLI F J. Optical heating in semiconductors: Laser damage in Ge, Si, In Sb, and GaAs[J]. Journal of Applied Physics, 1980, 51(10):5513-5522. doi: 10.1063/1.327469
    [6]

    QI H F, WANG Q, ZHANG X, et al. Theoretical and experimental study of laser induced damage on GaAs by nanosecond pulsed irradiation[J]. Optics & Lasers in Engineering, 2011, 49(2):285-291. 
    [7]

    XUE H Z, LI W, ZHANG H T, et al. Experimental study about laser induced damage to photovoltaic detectors in vacuum[J].Laser Technology, 2006, 30(5):494-497(in Chinese). 
    [8]

    QIU D D, WANG R, CHENG X A, et al. Mechanisms research on continuous wave laser induced damage to solar cells[J].Chinese Journal of Lasers, 2011, 38(3):302006(in Chinese). doi: 10.3788/CJL201138.0302006
    [9]

    TIAN X Q, XIAO S, TAO Sh H, et al. Damage threshold of femtosecond pulse laser for silicon solar cells[J]. Infrared & Laser Engineering, 2014, 43(3):676-680(in Chinese). 
    [10]

    HUANG F, NIU Y X, WANG X F, et al. Calculation of thermal and mechanical effect induced by laser in optical window materials[J]. Acta Optica Sinica, 2006, 26(4): 576-580(in Chinese). 
    [11]

    YANG H, LU J, ZHOU D Y, et al. Experimental study of 1070nm continuous wave fiber laser irradiation effect on three-junction GaAs solar cell[J]. Laser Technology, 2017, 41(3): 318-321(in Chinese).
    [12]

    SUN H, XU J M, ZHANG H C, et al. Simulation of three-junction Ga As solar cell temperature field by continuous wave laser irradiation[J]. Laser Technology, 2018, 42(2): 239-244(in Chinese). 
    [13]

    LI B B, LI X J. Numerical simulation of photovoltaic cell temperature field of laser power beaming[J]. Laser Technology, 2017, 41(4):537-537(in Chinese). 
    [14]

    SHANGHAI INSTITUTE OF SPACE POWER-SOURCES. Physical power technology[M].Beijing: Science Press, 2015:45-47 (in Chinese).
    [15]

    SUN W, DONG E L, YANG J F, et al. Measurement of 3-D profile using micro fringe projection correlation method[J]. Acta Metrologica Sinica, 2007, 28(1):10-13(in Chinese). 
  • [1] 李明海柳爱国宋耀祖 . 激光放大介质温度场和热应力场的数值模拟. 激光技术, 2002, 26(2): 86-89.
    [2] 张潇允张巍夏盛强马遥金光勇 . 高功率激光辐照CFRP的温度场和应力场的数值分析. 激光技术, 2021, 45(5): 636-641. doi: 10.7510/jgjs.issn.1001-3806.2021.05.017
    [3] 周益军张永康周建忠史建国 . 激光诱发板材热应力成形技术及其研究现状. 激光技术, 2006, 30(4): 425-428.
    [4] 管一弘樊刚陈铁力王洪纲 . 激光淬火热致残余应力的数值分析. 激光技术, 1998, 22(5): 313-317.
    [5] 杨建坤曹丁象郑万国贺少勃袁晓东於海武韩伟 . 热容模式下片状激光介质瞬态温度及热应力分析. 激光技术, 2007, 31(2): 196-199.
    [6] 刘子昂石伟汪诚 . 激光冲击强化残余应力的数值模拟研究. 激光技术, 2017, 41(1): 1-5. doi: 10.7510/jgjs.issn.1001-3806.2017.01.001
    [7] 许伯强刘洪凯徐桂东徐晨光李俊敏 . 基于应力-位移混合有限元法的激光超声数值模拟. 激光技术, 2014, 38(2): 230-235. doi: 10.7510/jgjs.issn.1001-3806.2014.02.018
    [8] 李贝贝李小将 . 激光输能光电池温度场数值模拟. 激光技术, 2017, 41(4): 537-544. doi: 10.7510/jgjs.issn.1001-3806.2017.04.016
    [9] 王军荣闵敬春宋耀祖 . 高热负荷激光介质热应力场的分析与优化研究. 激光技术, 2004, 28(6): 575-579.
    [10] 刘顺洪周龙早万鹏腾 . 激光诱发热应力成形的研究. 激光技术, 2001, 25(4): 258-268.
    [11] 孙浩徐建明张宏超杨欢陆健 . 连续激光辐照三结GaAs太阳电池温度场仿真. 激光技术, 2018, 42(2): 239-244. doi: 10.7510/jgjs.issn.1001-3806.2018.02.019
    [12] 苗心向袁晓东程晓锋贺少勃郑万国 . 表面污染物诱导熔石英损伤的热力学数值模拟. 激光技术, 2011, 35(2): 285-288. doi: 10.3969/j.issn.1001-3806.2011.02.039
    [13] 赵麒白忠臣周骅陆安江张正平刘桥 . 拉盖尔-高斯光束作用下熔石英温度及应力研究. 激光技术, 2018, 42(1): 121-126. doi: 10.7510/jgjs.issn.1001-3806.2018.01.024
    [14] 沈显峰王洋姚进汪法根杨家林 . 直接金属选区激光烧结热应力场有限元模拟. 激光技术, 2005, 29(4): 343-346.
    [15] 刘顺洪吉巧杰扬晶 . 钢管激光弯曲成形的数值模拟. 激光技术, 2006, 30(4): 355-359.
    [16] 杨洪亮金湘中修腾飞费鑫江叶颖 . 钢/铝异种金属光纤激光焊接数值模拟. 激光技术, 2016, 40(4): 606-609. doi: 10.7510/jgjs.issn.1001-3806.2016.04.031
    [17] 李晓锋周昕卢熙伍波杨泽后陈涌周鼎富侯天晋 . 激光在烟雾中传输特性的数值模拟分析. 激光技术, 2010, 34(3): 381-384. doi: 10.3969/j.issn.1001-3806.2010.03.027
    [18] 汪建敏周群立姜银方张梦蕾程科升万里赵龑 . 中空激光冲击金属板料变形的数值模拟. 激光技术, 2012, 36(6): 727-730. doi: 10.3969/j.issn.1001-3806.2012.06.004
    [19] 高立张永康 . 镍-钴合金杆料激光冲击强化的实验和数值模拟. 激光技术, 2006, 30(5): 507-510.
    [20] 陈根余陈晓锋周聪李时春廖生慧 . 激光深熔焊细长小孔数值模拟与试验研究. 激光技术, 2015, 39(2): 170-175. doi: 10.7510/jgjs.issn.1001-3806.2015.02.005
  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  5460
  • HTML全文浏览量:  3916
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-21
  • 录用日期:  2019-06-19
  • 刊出日期:  2020-03-25

连续激光辐照三结GaAs太阳电池热应力场研究

    通讯作者: 陆健, lujian@njust.edu.cn
    作者简介: 谭宇(1993-), 男, 硕士研究生, 主要从事激光与物质相互作用方面的研究
  • 南京理工大学 理学院, 南京 210094
基金项目:  上海航天科技创新基金资助项目 SAST20161113

摘要: 为了研究1070nm光纤连续激光对三结GaAs太阳电池的热力作用, 使用COMSOL软件建立物理模型进行了数值仿真, 得到了不同激光功率密度激光作用下的热应力场。为了验证了热应力计算方法的正确性, 利用栅线投影法测量了功率为16.7W、辐照半径为1mm、辐照时间为10s的激光作用下电池表面的形变, 实验结果与模拟结果基本吻合。结果表明, 辐照半径为1.5cm、功率密度为16.7W/cm2的激光辐照20s时, 底电池中心温度刚好超过电池使用温度; 底电池中心等效应力为96.6MPa, 刚好超过底电池材料的屈服极限; 根据这一结果可推测电池失效与热应力导致的结构损伤有关。该数值模拟结果与实验结果为激光辐照太阳电池的热力效应研究提供了一定的理论依据。

English Abstract

    • 随着各国航空航天事业的发展,性能越来越高的各类飞行器需要更加有效率的主电源,三结GaAs太阳电池因其具有耐高温、抗辐射能力强、光电转换效率高等优点被广泛应用[1-2]。相较于传统的太阳能供能,激光供能具有能量密度大、不受天气时间影响、转换效率高的特点,然而功率过大的激光会使电池温度升高产生较大的热应力,电池输出特性会由于热效应的产生而降低,因此,通过仿真计算电池在连续激光辐照下的温度场和热应力场具有重要研究意义,仿真结果为进一步研究激光辐照太阳电池致其损伤效应提供参考。在使用激光进行无线传输方面,1991年,LANDIS[3]等人探究了地面的空间激光输能方式;2006年, 日本KINKI大学[4]利用激光实现了风筝、机器人及无人机的供能实验。在激光对半导体的损伤方面,早在1980年,MEYER等人就研究了激光对Ge, GaAs, Si等材料的损伤[5]; 2011年, QI等人研究了纳秒脉冲激光对GaAs材料的损伤[6]。在太阳电池损伤研究方面,2006年,XUE等人实验研究了对光伏探测器的破坏[7]; 2011年, QIU等人[8]使用CO2连续激光辐照单晶硅太阳电池结合损伤形貌分析了太阳电池的损伤机理;2014年, TIAN等人[9]研究了飞秒脉冲激光作用下硅太阳电池的损伤阈值。在激光作用下的靶材温度场及应力场研究方面,2006年,HUANG等人做了激光辐照下材料的热力效应解析计算研究[10]; 南京理工大学的YANG[11]和SUN[12]等人在2017年分别从实验与数值模拟探究了温度对三结砷化镓太阳电池的性能影响; 同年,LI等人也做了激光输能光电池的温度场模拟[13]。由于应力测量难度较大, 在热应力对太阳电池的性能影响方面的研究相对较少。

      作者利用多物理场耦合仿真软件COMSOL构建了物理模型,利用热传导方程与热弹塑性方程,仿真计算了与实验结果基本吻合的温度场结果,并在此基础上仿真了对应的热应力场,由于热应力难以测量,作者通过实验测量电池辐照区域的形变量并与数值计算中的形变量对比,验证了数值计算方法的正确性。

    • 模型假设采用1070nm连续激光,辐照方向垂直于电池表面,横截面强度为高斯分布。由于光束为中心对称,所以模型建立为2维轴对称模型。图 1为电池物理模型和各层尺寸示意图。图 1中,z轴为对称轴,r方向为太阳电池径向,a为太阳电池厚度,b为电池半径,其值为1.5cm。

      Figure 1.  Model of solar cell

    • 柱坐标系下各向同性均匀材料的热传导偏微分方程为:

      $ \begin{array}{*{20}{c}} {{\rho _n}{C_n}\frac{{\partial {T_n}\left( {r,z,t} \right)}}{{\partial t}} = \frac{1}{r}\frac{\partial }{{\partial r}}\left[ {r{\kappa _n}\frac{{\partial {T_n}\left( {r,z,t} \right)}}{{\partial r}}} \right] + }\\ {\frac{\partial }{{\partial z}}\left[ {{\kappa _n}\frac{{\partial {T_n}\left( {r,z,t} \right)}}{{\partial z}}} \right] + {Q_n}\left( {T,r,z,t} \right)} \end{array} $

      (1)

      式中, n=1,2,3;ρnCnκnTn(r, z, t)分别表示在t时刻第n层的材料密度、热容、热导率和温度场分布; Qn为热源项, T为温度场。

      表 1所示,由于GaInP2层和GaAs层对1070nm激光的吸收系数相较于Ge层的吸收系数很小,故可以假设只有底电池吸收能量,并将其吸收的热能作为热源[14]。热源项写为:

      $ \begin{array}{*{20}{c}} {Q\left( {T,r,z,t} \right) = {I_0}\left( {1 - R - \eta } \right)\alpha \left( T \right) \times }\\ {I\left( {r,z} \right)g\left( t \right)\exp \left[ { - \alpha \left( T \right)z} \right]} \end{array} $

      (2)

      Table 1.  Forbidden band width and absorbable wavelength range of each layer of battery

      material forbidden bandwidth/eV absorption wavelength range/nm
      GaInP2 1.85 350~700
      GaAs 1.42 700~880
      Ge 0.67 880~1750

      式中, I0为激光辐照功率密度,R为电池表面反射率,η为光电转换效率, I(r, z)与g(t)分别为入射激光能量的空间分布函数与时间分布函数,α(T)为底电池的激光吸收系数。

      $ \begin{array}{*{20}{l}} {\alpha \left( T \right) = 1.4 \times {{10}^4}\exp \left[ {2.81 \times \left( {1.16 + } \right.} \right.}\\ {\left. {0.67 - 0.83 + 3.9 \times {{10}^{ - 4}}T - 1.17} \right)} \end{array} $

      (3)

      高斯光束的空间分布函数为:

      $ I\left( {r,z} \right) = \left\{ {\begin{array}{*{20}{l}} {\exp \left( { - 2{r^2}/a_0^2} \right),\left( {0 < r \le {a_0}} \right)}\\ {0,\left( {r > {a_0}} \right)} \end{array}} \right. $

      (4)

      式中,a0为激光辐照半径。连续激光的时间分布函数为:

      $ g\left( t \right) = 1,\left( {0 < t \le \infty } \right) $

      (5)

      初始温度场为:

      $ {\left. {T\left( {r,z,t} \right)} \right|_{t = 0}} = {T_0} = 303{\rm{K}} $

      (6)

      除电池上表面外,不考虑其它表面的热辐射与热对流。电池上表面边界条件为:

      $ \begin{array}{*{20}{c}} { - {{\left. {\kappa \frac{{\partial T\left( {r,z,t} \right)}}{{\partial r}}} \right|}_{z = 0}} = - h\left[ {T\left( {r,z,t} \right) - {T_0}} \right] - }\\ {\sigma \varepsilon \left[ {{T^4}\left( {r,z,t} \right) - T_0^4} \right]} \end{array} $

      (7)

      式中, κ为导热系数, φ为热辐射率, σ为Stefan常数, h为传热系数。从初始温度场开始,通过时间、空间的网格划分,可以逐步求得模型中各点的瞬态温度分布。表 2中给出了Ge的热学参量与公式中各常数的值。

      Table 2.  Thermal parameters of Ge

      parameter value
      density ρ/(g·cm-3) 5.49
      thermal conductivity κ/(W·cm·K-1) 0.59
      heat capacity C/(J·K-1) 0.31
      thermal radiation rate φ 0.1
      Stefan constant σ/(W·cm-2·K-4) 5.67×10-14
      coefficient of heat transfer h/(W·cm-2·K-1) 1×10-5
    • 热应力场的计算同样采用2维轴对称模型,在得到激光辐照太阳电池的温度场后,根据施加于靶材的温度载荷求解靶材的热应力场。根据von-Mises硬化准则,在材料弹性阶段材料的应力与应变为线性关系,满足Hooke定律,卸载后材料的变形与内部应力能逐渐恢复至初始值,当等效应力σ超过初始屈服强度时材料产生塑性变形。研究太阳电池在工作状态的热力效应时,激光功率并不高,所以电池温升相对较小,且变形量也较小,故忽略塑性变形,可按照热弹性模型计算。

      2维轴对称模型中平衡方程为:

      $ \left\{ {\begin{array}{*{20}{l}} {\frac{{\partial {\sigma _r}}}{{\partial r}} + \frac{{\partial {\tau _{zr}}}}{{\partial z}} + \frac{{{\sigma _r} - {\sigma _\theta }}}{r} = 0}\\ {\frac{{\partial {\sigma _z}}}{{\partial z}} + \frac{{\partial {\tau _{rz}}}}{{\partial r}} + \frac{{{\tau _{rz}}}}{r} = 0} \end{array}} \right. $

      (8)

      式中, σr, σθ, σz, τrzτzr分别为径向应力、环向应力、轴向应力、剪应力和切应力。假设材料为各向同性,其热弹性本构方程为:

      $ \left\{ \sigma \right\} = \left\{ {\begin{array}{*{20}{l}} {{\sigma _r}}\\ {{\sigma _\theta }}\\ {{\sigma _z}}\\ {{\tau _{zr}}} \end{array}} \right\} = {\mathit{\boldsymbol{D}}_{\rm{e}}}\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _r}}\\ {{\varepsilon _\theta }}\\ {{\varepsilon _z}}\\ {{\gamma _{rz}}} \end{array}} \right\} - \left\{ {\begin{array}{*{20}{c}} {\frac{1}{{1 - 2\mu }}}\\ {\frac{1}{{1 - 2\mu }}}\\ {\frac{1}{{1 - 2\mu }}}\\ 0 \end{array}} \right\}E{\alpha _\text{t}}\Delta T $

      (9)

      式中, De为弹性矩阵,E为材料弹性模量,μ为材料泊松比, εr, εθ, εz, γrz分别表示径向应变、环向应变、轴向应变及剪应变。{}表示列向量。表 3中列出了Ge材料的部分力学参量。表中,f为假设屈服强度,αt为热膨胀系数。

      Table 3.  Mechanical parameters of Ge

      parameter value
      Yang’s module E/(N·cm-2) 1.31×107
      fracture module f/(N·cm-2) 9.3×103
      thermal expansion coefficient αt/(K-1) 5.92×106
      Poisson ration μ 0.21

      当材料应力超过屈服强度时, 必须考虑材料的塑性变形,考虑材料的塑性变形后,一个单元中的总应变增量表达式为:

      $ \left\{ {\Delta \varepsilon } \right\} = \left\{ {\Delta {\varepsilon _{\rm{t}}}} \right\} + \left\{ {\Delta {\varepsilon _{\rm{e}}}} \right\} + \left\{ {\Delta {\varepsilon _{\rm{p}}}} \right\} $

      (10)

      式中, {Δε}为总应变增量;{Δεt}为热应变增量;{Δεe}为弹性应变增量;{Δεp}为塑性应变增量。

      弹性阶段中应力增量与应变增量关系可表示为:

      $ \left\{ {\Delta \sigma } \right\} = {\mathit{\boldsymbol{D}}_{\rm{e}}}\left( {\left\{ {\Delta \varepsilon } \right\} - \left\{ {\Delta {\varepsilon _{\rm{t}}}} \right\} - \left\{ {\Delta {\varepsilon _0}} \right\}} \right) $

      (11)

      其中,

      $ \left\{ {\Delta {\varepsilon _0}} \right\} = \frac{{\partial \mathit{\boldsymbol{D}}_2^{ - 1}}}{{\partial T}}\left\{ \sigma \right\}{\rm{d}}T $

      (12)

      当材料等效应力超过屈服强度时材料发生塑性屈服产生塑性应变增量,根据表 3中的参量假设Ge实际屈服强度为93MPa。再根据Prandtl-Reuss塑性流动增量理论,热弹塑性本构方程如下:

      $ \begin{array}{*{20}{c}} {\left\{ {\Delta \sigma } \right\} = {\boldsymbol{D}_{{\rm{ep}}}}\left( {\left\{ {\Delta \varepsilon } \right\} - \left\{ {\Delta {\varepsilon _{\rm{t}}}} \right\} - } \right.}\\ {\left. {\left\{ {\Delta {\varepsilon _0}} \right\}} \right) + \left\{ {\Delta {\varepsilon _0}} \right\}} \end{array} $

      (13)

      式中,Dep为弹塑性矩阵。

    • 根据上述模型,不同激光功率下电池上表面中心点如图 2所示。在特定功率密度激光辐照下,底电池上表面中心点的温度在辐照时间内逐渐升高,在辐照时间结束即20s时达到最大值,然后逐渐下降直至恢复至室温。随着激光功率密度越大,以下3个数值也越大:20s前的温度上升速率,20s时的温度最大值,辐照结束后温度下降速率。在功率密度为5W/cm2的激光辐照下,20s时中心点温度为343.9K, 温度场模拟结果与温度场测量实验结果[8]对比如图 2所示。此时电池的温升较小,现有研究表明,在该功率的激光作用下,太阳的电池的工作性能基本不受影响。

      Figure 2.  Comparison curve between experimental and calculated values of battery center temperature under 5W/cm2 power density laser irradiation

    • 研究表明,当激光功率密度超过5W/cm2时,电池工作性能会受到影响,因此本文中选取8.4W/cm2,11.7W/cm2和16.7W/cm2 3个激光功率密度研究电池热应力。中心点等效应力随时间变化曲线如图 3所示。在辐照过程中中心点应力不断增大,在20s达到最大值,辐照结束后随着温度降低中心应力也开始下降并逐渐趋向于零。而随着激光功率的增大,电池的热能积累更快,因此中心点热应力增长速率也更快。激光功率密度越高,任一时刻中心等效应力也越大。在8.4W/cm2和16.7W/cm2激光功率密度辐照下,20s时中心点应力分别为43.4MPa和96.6MPa,表现为压应力。由于锗晶体的断裂模量为93MPa,因此在16.7W/cm2激光功率密度辐照20s时将产生结构破坏,对比图 4中该功率下中心点的温度值,20s时中心点温度达到476.6K,基本刚好超过该电池的使用温度473K,此时电池的光电转换效率基本为0,因此推测当底电池中心应力超过其材料屈服强度时,电池彻底失效,这也进一步从侧面验证了模拟结果的正确性。

      Figure 3.  Equivalent stress curve of center point under laser irradiation of different power densities

      Figure 4.  Temperature curve of center point under laser irradiation of different power densities

    • 较低功率密度激光辐照下电池表面变形小于可测量值,为了使形变量便于测量以便验证数值计算结果,调整激光能量为16.7W,光斑半径为1mm,辐照时间为10s,在此条件下数值计算得到电池上表面中心点位移量如图 5所示。同热弹性模型相同,辐照结束后形变有一定的回弹,最终中心点位移稳定于约42μm。

      Figure 5.  Displacement curve of center point of battery surface under laser irradiation with power of 16.7W, radius of 1mm and irradiation time of 10s

    • 栅线投影法是3维形貌测量中的主要方法,具有较高的精度。本文中将提前设置的光栅通过投影仪投射到经过处理的电池表面,并从与投影方向呈一定角度的另一方向记录栅线条纹图像,激光辐照后条纹图像发生偏移,通过偏移量可得出物体上表面各点的高度变化量[15]

    • 本实验实验装置包括光纤激光器、扩束装置、测量设备、投影仪、拍摄设备。激光辐照系统如图 6所示。光纤连续激光器波长选取1070nm,输出功率连续可调。扩束装置为凸透镜L1与L2及可调光阑组成,L1直径为25mm,焦距为40mm;凸透镜L2的直径为50mm,焦距为150mm。通过调节L1与L2之间的距离可以连续调节入射到电池表面的光斑大小。投影仪与相机均为普通商用投影仪与相机。

      Figure 6.  Beam path of solar cells under laser irradiation

      为了得到清晰的投影图像,辐照前对电池表面喷涂了少量白漆。激光经扩束后经过光阑照射到样品上,光斑直径大小为2mm。将辐照前的平整电池表面光栅作为参考图像,如图 7所示。辐照后电池表面光栅如图 8所示。

      Figure 7.  Battery surface image before irradiation

      Figure 8.  Battery surface image before irradiation

    • 经激光照射后辐照区域内有明显凹陷,辐照区域边缘有轻微突起,经分析计算得到的电池辐照中心变形如图 9所示。辐照中心凹陷值为25μm,与数值模拟结果有一定差异,误差主要有以下来源:(1)半导体材料为各向异性材料,将材料考虑为各向同性的热弹塑性模型与实际情况有差异; (2)材料在温度上升到一定程度后,屈服极限、光电转换效率等参量会有较大变化,数值模拟过程中缺乏这些参量的实际数值变化规律; (3)激光的吸收虽然集中于底电池,但其它各层的形变对于上表面形变也会产生影响,数值模拟缺乏其它层的部分热力参量而与实际情况有差异; (4)喷漆在激光作用下对电池的表面位移产生了影响,并且喷漆对表面形变有一定掩盖; (5)实验过程中的环境,光栅图像的识别误差,激光入射功率的测定误差等其它因素引入的误差。

      Figure 9.  Deformation diagram of battery central area under irradiation

    • 利用COMSOL有限元软件构建了2维轴对称模型,在已有的温度场模拟结果与实验结果上,数值计算了不同功率密度的1070nm连续激光辐照下三结GaAs太阳电池应力场,并利用栅线投影法实验测量了电池中心的形变量,验证了数值计算方法的正确性。在激光功率密度为16.7W/cm2、辐照半径1.5cm的激光辐照20s时,电池中心温度超过使用温度,底电池中心等效应力也刚好超过材料屈服极限,可认为电池失效与热应力导致的结构破坏有关。通过实验测量更高功率激光辐照下电池表面的形变侧面验证了数值计算方法的正确性。该实验与模拟结果为进一步研究激光损伤太阳电池机理提供参考。

参考文献 (15)

目录

    /

    返回文章
    返回