Study of OCT weld depth curve fitting method for laser keyhole welding
-
摘要: 为了解决激光焊接熔深在线检测及光学相干光层析成像(OCT)熔深检测信号噪声多的问题，采用了一种百分位滤波加移动平均算法的熔深曲线拟合方法。通过光学相干光层析成像系统对304不锈钢匙孔焊过程进行了熔深在线检测，在使用百分位滤波从原始散点数据中提取出百分位熔深曲线过程中，发现该方法提取的百分位熔深曲线存在尖峰干扰噪声，通过移动平均算法对百分位熔深曲线进行低通滤波处理，最终从原始散点数据中提取拟合出OCT熔深测量曲线。结果表明，将OCT熔深测量曲线与焊缝纵抛面的实际熔深曲线进行对比，该方法获取OCT熔深测量曲线的测量精度最大提升了15%。该方法能从OCT原始散点数据中有效提取拟合出更为准确的OCT熔深测量曲线。Abstract: In order to solve the online detection of laser welding depth of fusion and to solve the problem of optical coherence tomography (OCT) depth of fusion detection signal noise, a percentile filtering plus moving average algorithm was used to fit the weld depth curve. The weld depth of 304 stainless steel keyhole welding process was detected online by optical coherence photomicrographic imaging system, in the process of extracting the percentile depth curves from the original scatter data using percentile filtering; it was found that there were spikes interfering noises in the percentile depth melting curves extracted by this method, and the percentile depth curves were processed by low-pass filtering through the moving average algorithm, and the OCT weld depth curves were finally extracted and fitted from the original scatter data. The results show that, comparing the OCT weld depth curve with the actual weld depth curve on the longitudinal surface of the weld, the measurement accuracy of the OCT depth of fusion curve obtained by this method has been improved by 15% at most. The method is effective in extracting and fitting more accurate OCT weld depth curves from the original scatter data.
-
0. 引言
空间相机的分辨率与口径成正比，随着科技发展，对于空间大口径反射镜的需求日趋明显。随着反射镜口径的增大，反射膜镀制的难度也急剧提升。热蒸发镀膜过程中存在晶格失配、热膨胀系数失配，这使得完成镀膜工作后，由镀膜温度冷却至室温过程中在“薄膜-基板”系统中存在热应力。当应力过大时会导致反射镜面形产生改变从而严重威胁到“薄膜-基板”系统的稳定性。另外，对于大口径反射镜来说，随着口径的增大，重力对于反射镜面形的影响变得更为明显，反射镜的面形发生变化将极大地增加反射膜的沉积难度。本文中针对热应力及重力对大口径反射镜镀膜的影响进行了仿真分析研究。
1. 生长热应力及重力对反射膜面形的影响
1.1 生长热应力产生的影响
为保证大口径光学元件的面形精度，并避免不必要的风险，镀膜过程中必须将基底温度控制在较低的水平[1-5]。在完成镀膜后，反射镜从镀膜温度冷却到室温后，仍会有一定的热应力存在于反射膜与反射镜基底之间[6]，它对反射膜的牢固度会产生很大的威胁。
一般将薄膜应力分为张应力和压应力，如图 1所示。当表现为张应力时，薄膜会产生收缩的趋势。当张应力超过薄膜本身的弹性限度时，会导致膜层破裂甚至脱落。压应力则与之相反，当压应力超过薄膜的弹性限度时，会导致膜层起皱[7-11]。对于Ag膜反射镜，采用热蒸发镀膜方式，从沉积温度冷却至室温的过程中，应力主要表现为张应力，具体表现是薄膜更容易在反射镜边缘部分出现脱落及破损的情况。
![图 1 薄膜受到张应力和压应力示意图]()
图 1 薄膜受到张应力和压应力示意图Figure 1. Schematic diagram of film subjected to tensile stress and compressive stress1.2 重力对反射膜面形的影响
到目前为止，大口径反射镜主流的镀膜方式仍是采用镜面向下、向上蒸发镀膜材料的方法[12-14]。通过优化反射镜背部固定点的位置，经研究发现，这种镀膜方式对口径在2 m左右的反射镜面形产生的影响仍在可控范围之内[15-18]，随着反射镜口径继续增大，由于重力的影响，反射镜面形出现明显的变化，若继续采用镜面向下的镀膜方式，会极大增加镀膜的难度[19]，同时严重影响反射膜的牢固度[20]。另外，完成镀膜工作后，反射镜需要进行吊装和镜面翻转[21]，在这个过程中，反射镜面以及沉积在镜面上的反射膜，受到的重力方向会产生改变，这种改变带来的反射镜面形变化对于大口径反射镜的反射膜伤害很大，容易造成膜层破裂以及脱落[22-23]。
2. 生长热应力对反射膜影响的仿真研究
2.1 反射膜热应力系统仿真建模
建模工况模拟真实镀膜情况，采用ZZS-2500型箱式真空热蒸发镀膜机，反射镜固定在镀膜机中的镀膜工装上，镀膜工装可绕中心轴旋转且高度可调。旋转速度提高有利于提高膜厚的均匀性，考虑到镀膜机旋转轴的承载能力，最大转速为15 r/min。镀膜时, 蒸发源距离反射镜垂直高度为1.4 m，镀膜工装以15 r/min速度转动，通过修正挡板对膜厚的均匀性进行调整。
针对不同口径反射镜反射膜生长热应力的分布及影响规律进行了以下仿真实验。Ag膜反射镜结构如图 2所示。反射镜基底口径分别设置为600 mm、1200 mm以及2400 mm，对应的基底厚度分别为60 mm、120 mm和240 mm(满足基底厚度∶反射镜口径=10 ∶1的要求)。Ag膜层厚度为100 nm，镀膜温度为200 ℃，环境温度为20 ℃，金属Ag膜的材料参数见表 1。
表 1 Ag膜材料参数Table 1. Material parameters of Ag filmdensity ρ/(kg·m-3) elastic modulus E/GPa Poisson’s ratio μ thermal expansion coefficient A/10-6K-1 coefficient of heat transfer K/(W·m-2·K-1) thermal conductivity κ/(W·m-1·K-1) 10500 73.2 0.38 19.5 429 411 为便于仿真计算，设定基底材料和Ag反射膜是均匀、各向同性的弹性材料，且不考虑周围环境影响。镀膜后，反射镜退火过程在真空镀膜机中完成，与外界无导热及对流。为研究热应力对反射膜本身所产生的面形变化趋势，针对反射膜本身进行仿真，Ag膜反射镜系统模型如图 3所示。基底下表面为z=0 mm平面，原点为基底下表面中心，图中金属灰色为单层银膜，深褐色为基底，划分网格后对模型进行有限元分析。
2.2 生长热应力仿真结果及分析
600 mm、1200 mm、2400 mm口径反射镜每个反射镜的吊装固定点数量为4个，分别设置为反射镜半径的中心位置处，反射镜基底温度设定为室温20 ℃。反射膜热应力仿真结果分别如图 4~图 6所示。从图中可以看出，靠近反射镜边缘部分的反射膜在冷却过程中温度梯度的变化更为剧烈，在冷却过程中释放出更多的应力，因此靠近反射镜边缘位置处，膜层更容易出现剥落和损坏的情况。
![图 4 600 mm口径反射镜反射膜热应力分布]()
图 4 600 mm口径反射镜反射膜热应力分布Figure 4. Thermal stress distribution of reflective film of 600 mm aperture reflector![图 5 1200 mm口径反射镜反射膜热应力分布]()
图 5 1200 mm口径反射镜反射膜热应力分布Figure 5. Thermal stress distribution of reflective film of 1200 mm aperture reflector![图 6 2400 mm口径反射镜反射膜热应力分布]()
图 6 2400 mm口径反射镜反射膜热应力分布Figure 6. Thermal stress distribution of reflective film of 2400 mm aperture reflector根据分析，原因在于应力的变化对反射膜的非球形形变有很大影响，非球形形变不仅发生在反射膜的边缘，而且布满整个镀膜区域，不论应力分布如何，发生在边缘位置的非球形形变，都是随着自由区的增大而增加的，反射镜口径的增大导致边缘部分自由区的增大，因此, 大口径反射镜更容易在边缘区出现膜层破裂、剥落等情况，仿真结果与参考文献[23]中的结论相符。
3. 重力对大口径反射镜基底形变影响的仿真分析
电子束镀膜过程中，反射镜镜面向下，一般采用多点吊挂或套圈吊挂方式。在镀膜实验研究中发现，基底材料不同的反射镜重力相差很大，反射镜自身重力所导致的面形变化会对镀膜工作带来很大影响。另外, 完成反射膜镀制工作后，反射镜的吊装、翻转过程使反射镜所受重力的方向发生变化，重力方向改变会引起反射镜面形产生变化。大口径反射镜由于其口径大、重量大，受重力影响反射镜面形的变化也更为明显，反射膜也有更大的可能出现损坏和脱落。为充分了解大口径反射镜镀膜及吊装、翻转过程中的重力分布，解决反射镜由于其自身重力导致的面形变化对镀膜工作带来的影响，对相同口径、不同基底材料的反射镜进行了重力仿真分析。
3.1 大口径反射镜-膜系统重力仿真建模
由于反射镜口径较大，且反射膜与基底之间的厚度相差过大，反射膜层本身有重力引起的面形变化依赖于基底面形由重力所导致的改变，因此重力仿真针对反射镜基底进行。
反射镜基底材料分别为石英、零膨胀(ultra-low expansion，ULE)玻璃、微晶这3种常用于制备空间反射镜基底的材料，材料参数如表 2所示。
表 2 基底材料及参数Table 2. Base material and its parametersserial number base material density ρ/(kg·m-3) Young’s modulus E/GPa Poisson’s ratio μ 1 quartz 2150 72 0.15 2 ULE 2210 67.6 0.17 3 crystallite 2500 95 0.25 以微晶玻璃基底反射镜为例，反射镜采用微晶玻璃作为基底材料，口径1320 mm、厚度132 mm。设基底下表面为z=0 mm平面，原点为基底下表面中心，反射镜背面外缘均匀分布6个吊挂点，反射镜基底模型如图 7所示。
![图 7 1320 mm口径微晶玻璃基底反射镜建模]()
图 7 1320 mm口径微晶玻璃基底反射镜建模Figure 7. Model of 1320 mm aperture reflector on glass ceramics substrate图 8为有限元模型网格划分后，对反射镜模型施加方向向下(与重力方向相同)的重力载荷并求解。
![图 8 网格划分后1320 mm口径反射镜有限元模型]()
图 8 网格划分后1320 mm口径反射镜有限元模型Figure 8. Finite element model of 1320 mm aperture reflector after grid division3.2 重力仿真结果及分析
3.2.1 反射镜镜面向下边缘吊挂方式镀膜
如图 9所示，6点吊挂方式，反射镜面形变化区域主要在镜子中心，由于重力作用，反射镜中心沿重力方向下移，反射镜基底产生明显形变。如果不消除重力影响，会对最终镀膜产品的面形、膜层牢固度等质量产生影响。
![图 9 1320 mm口径微晶玻璃基底反射镜6点吊挂重力影响仿真]()
图 9 1320 mm口径微晶玻璃基底反射镜6点吊挂重力影响仿真Figure 9. Simulation of 6 point hanging gravity influence of 1320 mm aperture reflector on glass ceramics substrate对1320 mm口径反射镜采用背部均布6个吊挂点的方式进行重力影响仿真，基底材料为空间反射镜常用的石英、微晶、ULE玻璃，重力仿真结果如图 10~图 15所示。反射镜背部形变主要集中在吊挂点位置处，该位置在镀膜过程中固定无位移，反射镜其余位置受重力影响发生形变。随着基底材料密度的增大，反射镜重量随之增大，反射镜背部边缘位置附近产生的形变也随之增大，镜面支撑点内部区域明显小于边缘处。
![图 10 1320 mm口径石英基底反射镜背部支撑点及反射镜面形变]()
图 10 1320 mm口径石英基底反射镜背部支撑点及反射镜面形变Figure 10. Back support point and specularity deformation of 1320 mm aperture reflector on quartz substrate![图 11 1320 mm口径石英基底反射镜背部支撑点及反射镜面应力变化]()
图 11 1320 mm口径石英基底反射镜背部支撑点及反射镜面应力变化Figure 11. Back support point and specularity stress change of 1320 mm aperture reflector on quartz substrate![图 12 1320 mm口径ULE基底反射镜背部支撑点及反射镜面形变]()
图 12 1320 mm口径ULE基底反射镜背部支撑点及反射镜面形变Figure 12. Back support point and specularity deformation of 1320 mm aperture reflector on ULE substrate![图 13 1320 mm口径ULE基底反射镜背部支撑点及反射镜面应力变化]()
图 13 1320 mm口径ULE基底反射镜背部支撑点及反射镜面应力变化Figure 13. Back support point and specularity stress change of 1320 mm aperture reflector on ULE substrate![图 14 1320 mm口径微晶基底反射镜背部支撑点及反射镜面形变]()
图 14 1320 mm口径微晶基底反射镜背部支撑点及反射镜面形变Figure 14. Back support point and specularity deformation of 1320 mm aperture reflector on glass ceramics substrate![图 15 1320 mm口径微晶基底反射镜背部支撑点及反射镜面应力变化]()
图 15 1320 mm口径微晶基底反射镜背部支撑点及反射镜面应力变化Figure 15. Back support point and specularity stress change of 1320 mm aperture reflector on glass ceramics substrate相同口径、不同基底材料反射镜重力引起的应力仿真结果与面形仿真结果类似。但图 14中微晶基底反射镜面应力分布与其余3种基底镜面应力分布略有不同，这与微晶材料的泊松比的值更大有关，导致微晶基底材料较其余3种基底材料更容易产生横向变形。
由重力仿真数值分析可知，受反射镜自身重力影响，采用6点边缘吊挂、镜面向下的镀膜方式时，镜面会在重力方向发生形变；这一形变与6点中间区域吊挂产生的形变可以互补。
3.2.2 反射镜镜面向上、底部支撑方式镀膜
大口径反射镜由于镜坯价格昂贵、加工周期长，镀膜过程中反射镜的安全是镀膜最基本要求。光学加工过程一般是镜面向上加工，如果采取镜面向下镀膜，需要进行翻转操作，反转过程会使镜体受力不均衡，如果设计或操作不当容易对镜体产生损坏。
对于大口径反射镜采用镜面向上的镀膜方式能够大大减小镜面翻转带来的风险。镀膜时可采取边缘套圈或多点支撑方式，图 16是3000 mm口径反射镜镜坯外环有12个区域支撑，重力环境下的镜坯变形与应力仿真图。由图可以看出，只进行边缘支撑，同样会引起镜坯产生形变，需要通过仿真设计，优化重力支撑点分布，减小重力影响，如图 17所示。
![图 16 3000 mm口径微晶基底反射镜背部支撑点及反射镜面应力变化]()
图 16 3000 mm口径微晶基底反射镜背部支撑点及反射镜面应力变化Figure 16. Back support point and specularity stress change of 3000 mm aperture reflector on glass ceramics substrate![图 17 3000 mm口径微晶基底反射镜重力消除背部支撑点分布示意图]()
图 17 3000 mm口径微晶基底反射镜重力消除背部支撑点分布示意图Figure 17. Distribution diagram of back support points for gravity elimination of 3000 mm reflector on microcrystalline substrate4. 结论
分别对600 mm、1200 mm、2400 mm不同口径的Ag膜反射镜进行了镀膜热应力有限元仿真，仿真结果表明，靠近反射镜边缘部分的反射膜在冷却过程中温度梯度的变化更为剧烈，在冷却过程中释放出更多的应力，因此在膜层靠近边缘的位置处更容易出现剥落和损坏的情况。另外，随着反射镜口径的增大，热应力也随之增大，因此, 相对于小口径反射镜来说，大口径反射镜的反射膜层更容易出现剥落和损坏的情况。
其次进行了重力对反射镜形变影响的仿真。仿真结果表明，对于大口径反射镜镀膜过程中，如果采取镜面向下的吊挂方式，随着反射镜口径增大，重力会引起反射镜形变，需要进行重力消除处理；通过对相同口径反射镜背部吊挂方式的重力影响仿真，设计了镀膜过程中的消应力背部吊挂装置，通过调整，可以达到减少重力对大口径反射镜镀膜过程基底形变的影响。应用重力仿真研究成果设计出1320 mm口径反射镜镀膜重力卸载工装。对于3000 mm口径的反射镜，在镀膜过程中，为提升膜层和镀膜后的面形精度，抵消重力影响，同样需要对支撑方式进行优化设计。
-
![]()
图 1 OCT熔深检测系统原理图
Figure 1. Schematic diagram of OCT laser weld depth monitoring system
![]()
图 4 样例A的焊缝纵抛面的金相图片
Figure 4. Metallographic picture of the longitudinal throwing surface of the weld of sample A
![]()
图 6 OCT熔深原始数据与实际熔深曲线的对比
Figure 6. Comparison of OCT raw data and actual depth of fusion curves
![]()
图 7 不同百分位数的熔深曲线与实际熔深的局部对比
Figure 7. Local comparison of weld depth curves at different percentiles with actual weld depth curves
![]()
图 8 百分位数p=5%的熔深曲线与实际熔深的全局对比
Figure 8. Global comparison of weld depth extracted by percentile filter p=5% with actual weld depth
![]()
图 9 基于百分位数p=5%和移动平均算法提取的熔深曲线
Figure 9. Weld depth curve extracted based on percentile filtering p=5% and moving average algorithm
![]()
图 10 拟合曲线与实际熔深曲线的差异评估
Figure 10. Evaluation of the difference between the fitted curve and the actual melt depth curve
![]()
图 11 样例B和样例C的金相图片、百分位滤波和基于百分位滤波及移动平均算法提取的熔深曲线
Figure 11. Metallographic images of sample B and sample C, weld depth curves extracted based on percentile filtering and weld depth curves extracted based on percentile filtering and moving average algorithm
表 1 3组实验样例的平均误差
Table 1 Average error of the three experimental samples
experimental sample percentile filter/% percentile filter+MA/% accuracy improvement/% sample A 4.93 4.21 15 sample B 4.97 4.34 13 sample C 4.53 4.09 10 
下载: 导出CSV
-
[1] XIAO X, LIU X, CHENG M, et al. Towards monitoring laser welding process via a coaxial pyrometer[J]. Journal of Materials Processing Technology, 2020, 277: 116409. DOI: 10.1016/j.jmatprotec.2019.116409
[2] SANDERS P G, LEONG K H, KESKE J S, et al. Real-time monitoring of laser beam welding using infrared weld emissions[J]. Journal of Laser Applications, Laser Institute of America, 1998, 10(5): 205-211. DOI: 10.2351/1.521853
[3] WANG L, MOHAMMADPOUR M, GAO X, et al. Adjustable ring mode (ARM) laser welding of stainless steels[J]. Optics and Lasers in Engineering, 2021, 137: 106360. DOI: 10.1016/j.optlaseng.2020.106360
[4] 李竹曼, 高向东, 张南峰. 大功率碟片激光焊接状态多特征融合分析法[J]. 激光技术, 2017, 41(5): 764-768. DOI: 10.7510/jgjs.issn.1001-3806.2017.05.029 LI Zh M, GAO X D, ZHANG N F. Analysis of high-power disk laser welding status based on multi-feature fusion[J]. Laser Technology, 2017, 41(5): 764-768(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2017.05.029
[5] FAN X, GAO X, LIU G, et al. Research and prospect of welding monitoring technology based on machine vision[J]. The International Journal of Advanced Manufacturing Technology, 2021, 115(11/12): 3365-3391.
[6] ZHANG X, TANG Z, WU Y, et al. Progress in in situ X-ray imaging of welding process[J]. Review of Scientific Instruments, Melville: AIP Publishing, 2022, 93(7): 071501. DOI: 10.1063/5.0074042
[7] HUANG D, SWANSON E A, LIN C P, et al. Optical coherence tomography[J]. Science, 1991, 254(5035): 1178-1181. DOI: 10.1126/science.1957169
[8] DUPRIEZ N D, DENKL A. Advances of OCT technology for laser beam processing: Precision and quality during laser welding[J]. Laser Technik Journal, 2017, 14(4): 34-38. DOI: 10.1002/latj.201700021
[9] 刘逸飞, 苏亚, 姚晓天, 等. OCT无创血糖检测图像处理最优化方法研究[J]. 激光技术, 2023, 47(2): 178-184. DOI: 10.7510/jgjs.issn.1001-3806.2023.02.004 LIU Y F, SU Y, YAO X T, et al. An optimization method of image processing for OCT non-invasive blood glucose detection[J]. Laser Technology, 2023, 47(2): 178-184(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2023.02.004
[10] MIYAGI M, KAWAHITO Y, KAWAKAMI H, et al. Dynamics of solid-liquid interface and porosity formation determined through X-ray phase-contrast in laser welding of pure Al[J]. Journal of Materials Processing Technology, 2017, 250: 9-15. DOI: 10.1016/j.jmatprotec.2017.06.033
[11] LOHAUS L, BAUTZE T, DIEPOLD K. Evaluation of optical sensors for laser welding in a technical cognitive environment[C]//International Congress on Applications of Lasers & Electro-Optics. Anaheim, California, USA: Laser Institute of America, 2010: 1541-1546.
[12] BLECHER J J, GALBRAITH C M, van VLACK C, et al. Real time monitoring of laser beam welding keyhole depth by laser interferometry[J]. Science and Technology of Welding and Joining, 2014, 19(7): 560-564. DOI: 10.1179/1362171814Y.0000000225
[13] DORSCH F, HARRER T, HAUG P, et al. Process control using capillary depth measurement[C]// International Congress on Applications of Lasers & Electro-Optics. San Diego, California, USA: Laser Institute of America Press, 2016: 1505.
[14] 谢冠明, 王三宏, 张跃强, 等. 基于光学相干层析的激光焊接熔深检测方法[J]. 光学学报, 2023, 43(11): 1114002. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202311014.htm XIE G M, WANG S H, ZHANG Y Q, et al. Laser welding depth monitoring method based on optical coherence tomography[J]. Acta Optica Sinica, 2023, 43(11): 1114002(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202311014.htm
[15] 黄威威, 游德勇, 高向东, 等. 基于相关分析和神经网络的激光焊接稳态识别[J]. 激光技术, 2022, 46(3): 312-319. DOI: 10.7510/jgjs.issn.1001-3806.2022.03.004 HUANG W W, YOU D Y, GAO X D, et al. Laser welding steady status recognition method based on correlation analysis and neural network[J]. Laser Technology, 2022, 46(3): 312-319(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2022.03.004
[16] SOKOLOV M, FRANCIOSA P, SUN T, et al. Applying optical coherence tomography for weld depth monitoring in remote laser welding of automotive battery tab connectors[J]. Journal of Laser Applications, Laser Institute of America, 2021, 33(1): 012028.
[17] DUIN R, HARINGA H, ZEELEN R. Fast percentile filtering[J]. Pattern Recognition Letters, 1986, 4(4): 269-272.
[18] BOLEY M, FETZER F, WEBER R, et al. Statistical evaluation method to determine the laser welding depth by optical coherence tomography[J]. Optics and Lasers in Engineering, 2019, 119: 56-64.
[19] MITTELSTÁDT C, MATTULAT T, SEEFELD T, et al. Novel approach for weld depth determination using optical coherence tomography measurement in laser deep penetration welding of aluminum and steel[J]. Journal of Laser Applications, 2019, 31(2): 022007.
[20] SCHMOELLER M, STADTER C, LIEBL S, et al. Inline weld depth measurement for high brilliance laser beam sources using optical coherence tomography[J]. Journal of Laser Applications, Laser Institute of America, 2019, 31(2): 022409.
计量
- 文章访问数: 422
- HTML全文浏览量: 5
- PDF下载量: 7