超声焊接全状态频率跟踪算法.pdf

第 43 卷 第 1 期Vol.43,No.12024 年 1 月Journal of Applied AcousticsJanuary,2024 研究报告 超声焊接全状态频率跟踪算法姚 震1罗子伦1刘凯捷1丁荣杰1莫毅强1何 建2乌斯别斯基亚历山大3(1 广东工业大学广州510006)(2 广州市科普超声电子技术有限公司广州510520)(3 白俄罗斯国家科学院国际技术转移中心明斯克220012)摘要：超声金属焊接负载重且变化剧烈，易导致换能器出现无阻性点和电源频率误跟踪。基于梅森等效模型，提出一种能同时计算出谐振频率与反谐振频率，避免误跟踪的全状态频率跟踪算法。存在阻性点时，利用3个不同发波频率及其对应的相位差，算出跟踪目标频率；当换能器处于无阻性点状态时，把跟踪目标变为相位差最小点，从而实现全状态频率跟踪。通过仿真软件对算法进行仿真，算法能有效避免误跟踪。焊接试验结果表明，能较易实现谐振频率和反谐振频率的频率跟踪，在无阻性点的工况下焊接过程平稳。关键词：频率跟踪；全状态控制；误跟踪；超声电源；超声金属焊接中图法分类号:TB552文献标识码:A文章编号:1000-310X(2024)01-0190-08DOI:10.11684/j.issn.1000-310X.2024.01.022Full-state frequency tracking algorithm of ultrasonic weldingYAO Zhen1LUO Zilun1LIU Kaijie1DING Rongjie1MO Yiqiang1HE Jian2USPENSKIY Alexander3(1 Guangdong University of Technology,Guangzhou 510006,China)(2 Guangzhou Kepu Ultrasonic Electric Technology Co.,Ltd.,Guangzhou 510520,China)(3 International Technology Transfer Center of the National Academy of Sciences of Belarus,Minsk 220012,Republic of Belarus)Abstract:Ultrasonic metal welding loads are heavy and vary drastically,which can easily lead to non-resistivefrequency state in the transducer and supply frequency mis-tracking.This paper proposes an full-state fre-quency tracking algorithm based on the Mason equivalent model that can simultaneously calculate the reso-nant and anti-resonant frequencies to avoid mis-tracking.When the transducer has resonant frequency andanti-resonant frequency,the target frequency is calculated by using three different frequencies and their corre-sponding phase differences in welding process;when transducer is in non-resistive frequency state,the trackingtarget is changed to the point with the minimum phase difference,thus achieving full-state frequency tracking.The algorithm was simulated by simulation software and it was concluded that the algorithm can avoid thefrequency mis-tracking phenomenon.The welding test results show that the frequency tracking of resonantfrequency and anti-resonant frequency can be easily achieved,and the welding process is smooth when thetransducer is in the non-resistive frequency state.Keywords:Frequency tracking;Full-state control;Mis-tracking;Ultrasonic generator;Ultrasonic metal weld-ing2022-09-06收稿;2023-02-28定稿广东省科技计划项目(2021A0505030068)作者简介:姚震(1976),男,广东韶关人,高级实验师,博士,硕士生导师,研究方向:超声加工技术。通信作者 E-mail:第43卷 第1期姚震等：超声焊接全状态频率跟踪算法1910 引言超声焊接被应用于动力电池的极耳及线束焊接生产中，其焊接质量直接影响产品的整体性能、良品率和使用寿命1。动力电池极耳超声波金属焊接频率通常为20 kHz，焊接功率达数千瓦，焊接过程负载变化剧烈且焊接时间短2，需要更快的频率跟踪速度。当换能器失谐时，焊头的输出功率波动较大，容易出现虚焊或过焊等现象3。超声换能器失谐工作时无功功率增加，换能器的电流电压升高，导致频率进一步漂移，容易损坏换能器。换能器焊接失谐一般指工作在无阻性点状态或频率误跟踪。超声焊接系统为非线性时变系统，其传递函数建模过程复杂。数字超声电源往往采用变步长跟踪4、PID5、模糊PID6算法等，不需建立具体模型的算法，优点是控制易实现，满足一般情况下的工况。但也存在一些缺点，PID参数整定困难，出现误跟踪无法自动复位，在无阻性点状态时容易频率大幅波动。本文基于梅森等效电路，通过数学公式推导，提出了一种全状态的频率跟踪算法，能在有阻性点时选择谐振频率或反谐振频率作为追频目标；在无阻性点时工作在换能器电压电流最小相位差频率下。新算法能同时计算出谐振频率与反谐振频率的值，判断换能器是否存在阻性点，并能实现误跟踪的自复位。1 换能器无阻性点状态与频率误跟踪超声换能器主要完成电能转换为机械能的过程，通常采用机电等效法对其工作过程进行分析。图1为压电换能器的梅森等效模型(忽略介电损耗)。?CCLR图1压电换能器等效电路图Fig.1Equivalent circuit of the piezoelectrictransducer在图1中，L1为换能器的动态电感，表征其等效质量；R1为换能器的动态电阻，表征其机械内阻；C1为换能器的动态电容，表征其柔度7。三者串联的支路为换能器的动态支路。C0为换能器的静态电容，是换能器的夹持电容，C0所在的支路为换能器的静态支路。换能器的复阻抗Z 和发波角频率之间的关系为Z()=R1+j(L1 1/C1)(1 2L1C0+C0/C1)+jC0R1,(1)式(1)中，角频率=2f，f 为换能器的驱动频率。一般情况下，Z()=0存在两个解，分别是谐振频率fr和反谐振频率fa。在焊接过程，这两个频率随时间不断变化，而频率跟踪则是在这两个频率点中选择一个作为跟踪目标，使换能器始终在谐振状态工作，保证焊接效果。1.1换能器的无阻性点状态换能器的动态电阻R1主要受负载的影响，换能器的负载与焊接压力、模头纹路尺寸和被焊接的材料有关。通过式(1)可知：当R1满足式(2)时8，换能器处于无阻性点状态。R1(C0+C1C0)L1C0C1.(2)换能器的相位角与负载大小密切相关，图2为不同负载下换能器相位曲线变化图，其中负载2大于负载1。-90-60-300306090?/()?/Hz?20 kg?30 kg?1?2Hz?frfrfrfafafa图2不同负载下换能器相频特性曲线比较Fig.2Relationship between the frequency andphase of the transducer in operation under differ-ent loads由图2可知，随着负载加重，换能器的fr和fa在逐渐靠近，换能器的感性区也在逐渐缩小，最终突破临界点后，换能器的阻性点消失。此时，无论以何种频率发波，换能器都无法在谐振状态工作。在无阻性点状态工作时，输出功率会存在无功功率。1922024 年 1 月1.2频率的误跟踪原因分析误跟踪是指追踪的频率点发生变化。图3为换能器的相频特性曲线，依据换能器振动频率和阻抗特性的不同，整个曲线可分为4个区域。设换能器的当前工作频率为f，当f fr为第一个区域，此时换能器阻抗呈现容性，当增加时，换能器容性会越来越弱，最终呈现阻性；当fr f(f2r+f2a)/2为第二个区域，此时换能器阻抗呈现感性，当增加时，换能器感性会越来越强，最终达到整个换能器感性状态的峰值；当(f2r+f2a)/2 f fa为第四个区域，此时换能器阻抗呈现容性，当增加时，换能器容性会越来越强。传统的频率跟踪算法通过对起始跟踪频率所处区域来确定跟踪目标。当起始频率处于第一、第二区域时，跟踪目标为fr；当起始频率处于第三、第四区域时，跟踪目标为fa。在小负载的情况下，第二、第三区域比较宽，传统跟踪算法根据跟踪目标来设定不同的跟踪频率上下限，避免误跟踪的发生。-90-60-300306090?/()?/Hz?+2?fafrfr2fa2+图3换能器的相频特性曲线Fig.3 Schematics for the phase and frequency ofthe transducer由式(1)可以推导出fr、fa两个频率的间距为|fa fr|=12(C0+C1C0)2L1C0C1R21L21.(3)在超声金属焊接中，负载重且波动大，导致跟踪频率上下限难以合理设置。假设目标跟踪频率为fa，根据扫频结果，起始频率定位在第三、第四区域。焊接加载时，换能器的R1会随着负载增加而增加，连续的焊接会使换能器发热，导致换能器的C0上升。由式(3)可知，假设L1、C1不变，当C0、R1增加时，fr、fa的间距减小8，即第二、第三区域缩小。负载抖动和频率跟踪的延迟有可能导致电源的发波频率从第三区域跑到了第二区域，引发误跟踪(追fa变为追fr)。现有的变步长频率跟踪、PID、模糊PID等各种频率跟踪控制算法，一旦发生误跟踪就难以找回目标跟踪频率。图4(a)和图4(b)均为同次焊接时换能器电压电流波形，深蓝色线为换能器电压波形，浅绿色线为换能器电流波形。电源匹配电路为LC匹配，跟踪算法为变步长频率跟踪算法4，在焊接的全过程换能器均处于谐振状态。图4(a)中为跟踪到反谐振频率点时的电压电流波形，换能器电压达1380 V，电流是663 mA。电压高电流低是换能器处于fa的特征，此时焊接功率为915 W。图4(b)中为跟踪到谐振频率点的电压电流波形，电压降为844 V，电流上升到1780 mA。电压低电流高是换能器处于fr的特征，此时焊接功率上升到1502 W，功率上升近65%。由此可知，换能器的跟踪频率由fa误跟踪到了fr，出现这种情况容易导致功率波动，影响焊接一致性。?(a)?(b)?图4金属焊接频率误跟踪换能器电压电流波形图Fig.4 Waveforms of the voltage and current ap-plied on the transducer during the process of weld-ing the metal when its frequency was tracked in-accurately第43卷 第1期姚震等：超声焊接全状态频率跟踪算法1932 全状态频率跟踪算法的推导与实现2.1全状态频率跟踪算法的推导假设当前电源输出频率为f，角频率为，其对应的换能器两端电压电流相位差为。由于相位差角度与负载阻抗角相同4，由式(1)可以推导出tan与之间的关系为tan=1R1(L11C1)C0(R21+(L11C1)2).(4)在式(4)的基础上进一步化简可以得到以下表达式：tan=C0L213R1+2C0L1+C1L1 R21C1C0R1C1C0+C1R1C1.(5)为使表达式变得简洁，现进行以下等效变换：K1=C0L1R1,K2=2C0L1+C1L1 2C0L1R1C1,K3=C0+C1R1C21.(6)把式(6)代入式(5)后等式两边同时乘，可得 tan=K14+K22+K3.(7)不难看出，式(7)有4个解，由于负频率没有物理意义，舍去两个负数解后可得fr=12K2K22 4K1K32K1,fa=12K2+K22 4K1K32K1.(8)当阻性点不存在时，目标频率计算公式如下：f=12K22K1.(9)在式(7)和式(8)中，是当前的发波频率，tan则是换能器相位差的正切值，可选fa或fr作为目标跟踪频率。因此，只要求出K1、K2、K3三个未知数，即可同时得到换能器的谐振频率和反谐振频率。值得注意的是，当K224K1K3?fnfnfnnnn?(7),?fcfn/fcn l fn/fnsSnfcfnfmaxfn/fnfmaxfc?fc=real(fc图5全状态频率跟踪算法流程图Fig.5 Algorithm flow for the transducer was generated at any frequency1942024 年 1 月每间隔一段时间采集相位差，实时计算出驱动频率。图5中的f1为单次焊接的首个发波频率，同理fn则为第n次改变后的发波频率，其对应的换能器电压电流相位差为n。流程图中还存在有两个控制参数，分别为抖动系数s和最大限定步长fmax。两个控制参数会共同影响算法的跟踪速度和稳态误差。一般通过扫频获取换能器空载时的fa作为起始发波频率，并依照fn=fn1+s n%3的方式得出算法起始所需的3个f 值和值。后续驱动频率根据本次和前两次发波的频率及其值进行计算得出，计算结果与上次输出的差值再与fmax比较，选择较低的输出。当n小于阈值l时，fn将在fa附近进行小范围波动。2.3抖动系数和最大限定步长对算法的影响设置抖动系数和最大限度步长虽然影响算法跟踪速度和稳态误差，但增加了算法稳定性。从仿真中发现，当接近0时，每一次的跟踪步长将会变得很小，就会导致计算的3个发波频率非常接近，计算的目标频率会出现严重的跳变。另外运算器的精度及计算误差也会造成目标频率的跳变。抖动系数就是人为加入微小偏差，消除计算频率的跳变。通常换能器的R1是时变的，每次采样值与理想值之间存在误差，可通过最大限度步长来限制这种误差的影响范围。3 仿真与实验验证3.1仿真设计为了验证算法的可行性，本文采用仿真软件对全状态频率跟踪算法进行理论验证。仿真设计的需求为：(1)换能器负载以连续的方式变化。(2)反馈需要以离散的方式获得。(3)结果则需要以曲线的方式较为直观地展现出来。选取20 kHz超声焊接换能器，使用阻抗分析仪测试其等效参数，如表1所示。表1换能器等效参数Table 1 Equivalent parameters of transduceC0/nFC1/nFL1/mHR1/10.27790.2208288.98626.1007忽略焊接时的抖动影响，换能器R1总体是以先快后慢的方式进行变化的9，为了更接近实际负载对换能器的影响，设置R1以式(10)方式变化，其中参数K 可用于调节R1的变化速度。R1=50+K ln(1+1000000 t).(10)仿真设计思路如下：(1)假设换能器符合梅森等效模型，其等效参数初始值如表1所示。焊接开始后，换能器的R1随时间以式(10)规律变化。(2)设定超声电源的起始频率，每个控制周期输出一个目标值，此目标值为下一个控制周期电源的发波频率，把这个频率输入到模型中，根据当前的等效参数计算并返回换能器相位差值。(3)记录整个仿真过程中发波频率、fa和fr的值，将换能器相差角随时间的变化绘制成曲线，观察其在理想情况下，算法的跟踪速度、跟踪稳定性和跟踪精度。3.2仿真结果分析仿真时算法的控制周期为100 s，焊接时间选择0.05 s，所需要用到的反馈值tan则基于梅森等效电路通过理论计算获得，从而得到该发波频率下对应的值。仿真结果如图68所示。图6是以fr为跟踪目标时，换能器两端电压电流相位差、fr、fa和发波频率随时间的变化曲线。为了验证算法能在误跟踪时自复位，将跟踪的起始频率设置为20170 Hz(远离fr值，跨过fa值)。可以看出，发波频率在焊接开始后在5 ms左右跟踪到fr，并在短暂的抖动后完全跟上。在之后的焊接时间里，跟踪稳定且几乎没有稳态误差，跟踪精度达0.1 Hz。00.010.020.030.040.05?/s-80-60-40-2002040?/()1.981.992.002.012.022.03?/(104 Hz)?fafr图6以fr为跟踪目标的仿真结果Fig.6 Simulation results with fras the trackingtarget第43卷 第1期姚震等：超声焊接全状态频率跟踪算法195图7是以fa为跟踪目标时，换能器两端电压电流相位差、fr、fa和发波频率随时间的变化曲线。这里把将跟踪的起始频率设置为19830 Hz(远离fa值，跨过fr值)。虽然发波频率也可以在焊接开始后5 ms左右跟踪到fa，完成误跟踪的自恢复，但相对于跟踪目标为fr情况，其跟踪精度与跟踪速度均有下降，跟踪精度只有2 Hz。?/()00.010.020.030.040.05?/s-80-60-40-20020401.981.992.002.012.022.03?/(104 Hz)?fafr图7以fa为跟踪目标的仿真结果Fig.7 Simulation results with faas the tracking target图8为换能器在无阻性点状态时，换能器两端电压电流相位差、fr、fa和发波频率随时间的变化曲线。在后半段，随着R1的增长换能器进入无阻性点状态，fr和fa合并成一条曲线，为换能器容性最低点频率的曲线。在换能器进入无阻性点状态后，算法自动更换控制目标，在短暂的抖动后完成频率跟踪。但是，相对于阻性点存在的情况，其跟踪精度有所下降，只有4 Hz。在R1很高的情况下，换能器的相位随频率变化的灵敏度会下降，在相位曲线上仍具有较高的稳定性。?/()?/s?/(104 Hz)00.020.040.060.080.10-80-60-40-20020401.981.992.002.012.022.03?fafr图8无阻性点状态时算法跟踪仿真结果Fig.8Algorithm tracking simulation results inthe non-resistance point state3.3焊接实验验证为了验证算法的实际效果，搭建实验平台进行验证。其中电源的控制芯片为TMS320F28377S，主电路采用LC匹配10，焊件为3层0.2 mm铜镀镍金属片，焊接功率约为1100 W。实验中的电源驱动频率和相位差数据由超声电源通过串口导出，与TEKTRONIX示波器MDO3104测得数据一致。假定换能器电压电流相位差小于5为谐振状态，则驱动频率不改变。以谐振频率fr为跟踪目标的结果如图9所示，图9(a)为焊接时换能器两端电压电流波形图，图9(b)为驱动频率和相位差随焊接时间的变化曲线。跟踪的启动频率设置为20010 Hz，气压为0.2 MPa。-80-65-50-35-20-5102540552000020020200402006020080201002012020140201602018000.090.180.270.360.45?/Hz?/s?/s?00.10.20.30.40.5-2000-1500-1000-5000500100015002000?/V-3000-2000-1000010002000?/mA(a)?(b)?/(O)图9以fr为跟踪目标的实验结果Fig.9The experimental result when the algo-rithm uses fras tracking target图10是在跟踪目标为反谐振频率fa时的实验结果，其中图10(a)为焊接时换能器两端电压电流波形图，图10(b)为驱动频率和相位差随焊接时间的变化曲线。焊接工况与跟踪目标为fr时相同。为了验证算法的自复位效果，跟踪的启动频率同样设置为20010 Hz。1962024 年 1 月?/Hz-80-70-60-50-40-30-20-1001020200002002020040200602008020100201202014020160201802020000.090.180.270.360.45?/s?/s?00.10.20.30.40.5-2500-2000-1500-1000-50005001000150020002500-3000-2000-1000010002000?/V?/mA(a)?(b)?/(O)图10以fa为跟踪目标的实验结果Fig.10 The experimental result when the algorithm uses faas tracking target-80-70-60-50-40-30-20-10020000200302006020090201202015020180202102024000.090.180.270.360.45?/Hz?/s?/s?/(O)0.10.20.30.40.5-1500-1000-500050010001500-3000-2000-10000100020003000?/V?/mA(a)?(b)?图11无阻性点状态时算法跟踪实验结果Fig.11 The experimental results when the transducer is in the non-resistance point state把焊接时的气压从0.2 MPa提高到0.4 MPa，使得换能器处于无阻性点状态，测试算法在换能器处于无阻性点时是否能较为稳定的跟踪。图11为换能器处于无阻性点状态时的频率跟踪效果，其中图11(a)为焊接时换能器两端电压电流波形图，图11(b)为电压发波频率和换能器的相位差随焊接时间的变化曲线。上述实验结果表明，算法能够较易实现谐振频率或反谐振频率的频率跟踪；在无阻性点状态时仍能维持较小相位差。存在阻性点时可完成负载动态变化(持续约80 ms)的跟踪，随后的相位差基本保持在8以内，在换能器处于无阻性点时不会出现严重的频率失锁问题，跟踪效果稳定。4 结论本文主要研究超声焊接的频率跟踪技术，考虑了换能器的负载特性，提出全状态频率跟踪的算法。超声电源可以根据焊接对象选择工作在谐振或反谐振频率点。仿真结果和焊接实验均表明，算法能够避免误跟踪；启动阶段的起始频率与目标频率相差较大时，追频时间大约5 ms；在重载情况下，换能器谐振频率消失时，仍能够找到最佳工作频率稳定工作。本文提出的控制算法可应用于负载变化快的场合，也可为动力电池上连接件、线束等较厚工件的超声焊接提供参考。参考文献1 Zwicker M F R,Moghadam M,Zhang W,et al.Auto-motive battery pack manufacturingA review of batteryto tab joiningJ.Journal of Advanced Joining Processes,2020,1:100017.2 李正中,周光平,梁召峰,等.全数字控制超声波金属焊接电源的设计J.声学技术,2012,31(4):403406.Li Zhengzhong,Zhou Guangping,Liang Zhaofeng,et al.Design of digital power source for ultrasonic metal weld-ingJ.Technical Acoustics,2012,31(4):403406.第43卷 第1期姚震等：超声焊接全状态频率跟踪算法1973 彭呈祥,段发阶,蒋佳佳,等.基于频率跟踪的超声驱动电源研制J.电力电子技术,2019,53(5):15.Peng Chengxiang,Duan Fajie,Jiang Jiajia,et al.Devol-opment of ultrasonic driving power based on frequencytrackingJ.Power Electronics Synopsis,2019,53(5):15.4 侯光华,杜贵平,罗杰.超声波电源的改进频率跟踪方法 J.电源学报,2019,17(1):152158.Hou Guanghua,Du Guiping,Luo Jie.Improved frequencytracking method for ultrasonic power supplyJ.Journal ofPower Supply,2019,17(1):152158.5 刘丽晨,杨明,李世阳,等.超声换能器并联谐振频率的复合式跟踪方法研究J.应用声学,2015,34(1):4550.Liu Lichen,Yang Ming,Li Shiyang,et al.Parallelresonant frequency composite tracking of piezoelectrictransducerJ.Journal of Applied Acoustics,2015,34(1):4550.6 刘宁庄,张远宝,许龙,等.基于模糊 PID 控制的超声电源频率跟踪设计J.电力电子技术,2021,55(10):8082,86.Liu Ningzhuang,Zhang Yuanbao,Xu Long,et al.Designof ultrasonic power frequency tracking based on fuzzy PIDcontrolJ.电力电子技术,2021,55(10):8082,86.7 刘礼平,张朋朋,朱学明.基于负载变化的超声振子电学性能试验研究J.航空制造技术,2020,63(Z1):7380,86.Liu Liping,Zhang Pengpeng,Zhu Xueming.Experimen-tal study on electrical performance of ultrasound vibra-tor based on load changeJ.Aeronautical ManufacturingTechnology,2020,63(Z1):7380,86.8 武剑,董惠娟,张广玉.压电换能器锁相环频率跟踪的失效分析与解决 J.华南理工大学学报(自然科学版),2010,38(3):123128.Wu Jian,Dong Huijuan,Zhang Guangyu.Analysisand solution for frequency tracking degradation of phase-locked loop of piezoelectric transducersJ.Journal ofSouth China University of Technology(Natural ScienceEdition),2010,38(3):123128.9 姚震,莫毅强,罗子伦,等.非线性工作状态压电换能器匹配技术J.应用声学,2023,42(2):380386.Yao Zhen,Mo Yiqiang,Luo Zilun,et al.Matchingtechnology of piezoelectric transducer in nonlinear work-ing stateJ.Journal of Applied Acoustics,2023,42(2):380386.10 韩丽轩,于保华,胡小平.功率超声压电换能器阻抗匹配电路参数化设计J.压电与声光,2015,37(4):713716,720.Han Lixuan,Yu Baohua,Hu Xiaoping.Parametric designof impedance matching circuit for power ultrasonic piezo-electric transducersJ.Piezoelectrics and Acoustooptics,2015,37(4):713716,720.