�第22卷�第1期爆炸与冲击Vol.22,No.1��2002年1月EXPLOSIONANDSHOCKWAVESJan.,2002�文章编号:1001�1455(2002)01�0061�05光斑尺寸对K9玻璃近红外激光损伤阈值的影响罗�福1,杜祥琬2,孙承纬1(1.中物院流体物理研究所,四川绵阳�621900;2.北京计算物理与应用数学研究所,北京�100088)��摘要:基于热弹性模型,用有限元法对不同光斑近红外激光辐照下K9玻璃样品中的温度和压力分布进行了计算。对厚2mm、半径10mm的圆片样品的计算结果表明,K9玻璃的损伤受表面环向拉伸应力的控制,光斑大小对损伤阈值有较大的影响,对光束总入射功率为28kW的激光辐照,当光斑半径从0.4mm增加到3mm时,损伤所需时间从14ms增加到1s。在一定光斑尺寸范围内,样品损伤(表面拉伸断裂)所需的激光功率密度在同一量级内变化,在更小光斑时,损伤时的最高温升趋于565K。��关键词:K9玻璃;激光光斑效应;温度;热应力;损伤��中图分类号:O346���文献标识码:A�1�引�言��K9玻璃广泛用作光学仪器的窗口、棱镜以及反射镜、滤光片的基体。K9玻璃的近红外激光损伤,是光学系统抗激光辐照中的典型问题。在高功率激光辐照下其损伤机制为缺陷、夹杂物吸收引起的多光子吸收、自由电子加热及自聚焦效应,导致雪崩电离击穿破坏[1]或夹杂热爆破坏[2]等。由于大气的限制,缺陷、夹杂物吸收的激光能量不足以引发上述过程,其吸收的能量将转变为热能,扩散到玻璃介质中而表现为光斑内的体加热。夏晋军等[3]测量了K9玻璃在氧碘化学激光辐照下的破坏阈值,认为损伤的机制是热应力炸裂,但他们只考虑了一种光斑尺寸。对非均匀加热的情形,光学材料的损伤模式与激光束辐照样品的几何参数有关,几何参数包括样品尺寸、样品表面处激光光斑尺寸、热扩散距离等尺度的相对大小以及样品所受的边界条件,它对温度和热应力的时空分布有重要影响。如果辐照时间不长,边界影响不起作用,此问题通常被忽略。这是因为在一般的应用中,玻璃的激光损伤问题往往出现在高功率激光器件设计和激光加工领域,在器件设计方面可以用耐激光辐照性能更好的材料如熔石英等,而在激光加工中则会选用吸收效率高的激光如CO2激光来加工。��当激光较弱、辐照时间较长时,样品的边界条件和辐照光斑大小对玻璃样品的损伤机制和阈值有着不可忽视的影响。本文中,我们不用解析模型而用数值计算来分析这个问题。2�理论模型与计算模型��在较低功率和较长激光作用时间时,受辐照样品中的温度分布和应力分布可以解耦处理,即可先由导热微分方程求解温度场,再将各个时刻的温度场分别代入热弹性平衡方程求解应力场。在应变点(K9玻璃对应粘度1013.5Pa�s)以下,玻璃表现为弹性体,这时使用热弹性本构关系。��考虑厚l、半径rs的圆板样品,取柱坐标,坐标原点在样品的激光入射表面中心,z轴与激光照射方向一致,Gauss型空间分布的激光束垂直照射到样品表面(z=0),t=0时刻激光开始照射样品,初始时�收稿日期:2001�01�04;�修回日期:2001�03�15���作者简介:罗�福(1965��),男,博士,副研究员。时样品温度均匀,记为T0,样品上、下表面为冷却边界,样品侧面(r=rs)处绝热,并且为自由面,但环线r=rs、z=l的轴向位移为0。样品内部的温度场T(r,z,t)满足热传导方程1��T�t=�2T�r2+1r�T�r+�2T�z2+F0��g(t)e-(r2/R2)-�z(1)初始条件和边界条件为T(r,z,0)=T0��-k�T�n|z=0=h(T-T0)-k�T�n|z=l=h(T-T0)��-k�T�n|r=rs=0(2)式中:�为热扩散率,k为热导率,�为吸收系数,F0为进入玻璃的激光光斑中心的通量,R为Gauss光束半径,T0为初始温度,n为表面法向,g(t)为归一化的随时间变化的激光波形函数,这里为单位阶跃函数。h为表面换热系数,在真空环境中,玻璃只发生热辐射损失,故h=��B(T4-T40)/(T-T0)(3)式中:�为黑度,对K9玻璃,��0.8;�B为Stefan�Boltzman常数。��轴对称热弹性问题的平衡微分方程为��r�r+��rz�z+�r-��r=0����z�z+��rz�r+�rzr=0(4)本构关系为�r=E1+�(�1-2�e+�r)����=E1+�(�1-2�e+��)�z=E1+�(�1-2�e+�z)���rz=E2(1+�)�rz(5)几何变形关系为�r=�ur�r+�(T-T0)���=urr+�(T-T0)�z=�uz�z+�(T-T0)��rz=�ur�z+�uz�r(6)...
