专利名称：基于声音特征提取的超高压水射流速度优化控制方法
技术领域：
本发明涉及超高压水射流切割技术，具体是超高压水射流速度优化控制方法。
背景技术：
超高压水射流切割机在切割过程中具有无化学变化、无热变形、切缝窄、精度高、切面光洁、清洁无污染等优点，随着高科技的迅速发展，水射流切割机更趋向复杂化、更高压力、高性能、大负荷、高集成度、机电一体化、依赖性强等特点，因此，如何提高切割产品的光滑度和切割效率是目前超高压水射流切割优化控制的主要研究内容。
由水射流基本理论得到影响水射流切割光滑度
K工艺参数主要有喷嘴移动速度6、水压P、流量g、喷嘴结构参数(内径、长径比、材
料、形状)、靶距“冲击角度5、切割次数《、总功率况、磨料参数(种类、尺寸、流量)、磨
料混合方式(前置式、后置式)、产品参数(厚度、材料)等20多个，因此水射流系统光滑度&
与这些参数是强非线性映射关系，对这些参数进行全面定量的研究和优化控制目前存在很大的困难。然而，在水射流切割工业现场中，受生产条件限制，并不是上述所有的工艺参数
都可以调整，可以调整的主要的工艺参数有喷嘴移动速度G、水压P、冲击角度0等，其中
喷嘴移动速度i是主要优化控制因素。由生产经验得知不同厚度的产品切穿持续时间明显不同，切割时间越长，切割质量越好，对应喷嘴的横移速度越慢；通过水射流设备的数控系统上层指令可以方便调整，可行性好，喷嘴移动速度对提高切割的光滑度 < 的贡献率在18%。因此，如何根据不同的材料、不同厚度、不同光滑度要求自适应优化最佳的喷嘴移动速度z是关键。现有的超高压水射流速度控制仍然处于简单的人工经验控制，效率低，通过在数控系统中手工设置最大值和最小值速度，按照直线减速方式进行硬切换，难以满足柔性制造的高光滑度自适应控制要求。
发明内容
本发明正是为了解决上述问题，提供一种基于声音特征提取的自适应优化控制超高压水射流喷嘴横移速度的方法，提高切割产品表面光滑度和切割效率，满足实际工业生
产的需要。本发明解决其技术问题所采用的技术方案是包括如下步骤(I)在没有切割工件时采集多段背景噪声信号，对背景噪声信号进行分析获得时域特征和频域特征，背景噪声能量集中在0-8500HZ滤波频带基准内，背景噪声信号均方差为; (2)获取#组水射流
切割现场声音信号7 U)，选取*>令)的￡个离散序列/ )= 是正整数，对7ft)进行五个尺度的正交小波分解，获得各尺度的小波分解系数a = 1,2,...5为分解尺度，^=0,2,...2^为分解个数；根据小波分解系数的平方和计算各尺度上的分解能量，建立分解尺度和分解能量分布图，并利用所述背景噪声的均方差I和所述滤波频带基
准，基于最小距离原则确定尺度阈值A = 3，令大于该尺度阈值的各尺度的小波系数为零，其他尺度分解小波系数不变，得到新的小波系数，对所述新的小波系数进行信号小波逆变换，得到纯粹的水射流切割反射声音/ 的近似反射声音信号/4),求出/4)的振幅包络;当振幅包络,(X)的绝对值大于设定的幅值阈值门限A = 0.35时，获取对应的水射流切割开始时刻％ ;当振幅包络,￠)的绝对值小于设定的幅值阈值门限A =0.35时，获取对应的水射流切割结束时刻h ;得到#个水射流反射声音切穿时间HI，计算信号能量
^ = ElZi^f ; (3)将切穿时间L、能量5特征作为输入量、以人工经验获得的最佳喷嘴横 移速度为输出量建立最小二乘支持向量机模型；(4)切割工件时，在线采集切割现场声音信号，提取其切穿时间和信号能量，将切穿时间和信号能量输入到所述最小二乘支持向量机模型，获得最佳喷嘴横移速度。当超高压水射流冲击不同材质和几何形状参数的目标物体时，产生声音的特征值不同，因此本发明利用超高压水射流反射声音识别是否切穿物体，从而确定合适的优化速度。本发明具有很强的实用性，能够对不同材料、不同厚度、不同光滑度要求基于声音特征自适应优化最佳的喷嘴移动速度，可提高水射流切割的工作效率以及切割产品的品质，减轻操作者的劳动强度。
图I是声音信号计算机采集和控制系统应用 图2是图I中的计算机9的结构 图3是背景噪声信号 图4是现场声音信号 图5是现场声音信号的频谱 图6是纯粹的水射流切割反射声音信号包络 图7是本发明基于声音特征提取的超高压水射流速度优化控制流程；
图8是切割声音时间与最优切割速度的匹配关系图。
具体实施例方式参见图1，首先在现有的超高压水射流切割机上设置声音信号计算机采集和控制系统，现有的超高压水射流切割机包括升压泵系统I、加砂部分2和水射流喷嘴3，加工工件4固定在工作台5，声音信号计算机采集和控制系统由声音传感器6、声音信号传输线7、信号放大调理模块8以及计算机系统9依次相串接组成，其中，声音传感器6设置在加工工件4上，本发明可采用多个声音传感器6。参见图2，计算机系统9是基于TMS320DM642 DSP为核心的声音采集计算机系统，信号放大调理模块8采用AWA14604型前置放大器组成前端测试单元，通过信号传输线和高速A/D数据采集卡送入计算机I/O接口处理，信号采集频率为32000Hz。声音传感器6米集的现场声音信号由纯粹的水射流切割反射声音和背景噪声组成，本发明要对背景噪声进行有效抑制，再对纯粹的水射流切割反射声音进行特征提取。具体是
参见图7，首先，采用图I的声音信号计算机采集和控制系统提取背景噪声。背景噪声主要包括切割机升压泵系统I工作的噪音、工作台5的低频振动噪声等等。当水射流喷嘴3没有切割加工工件4时，通过声音传感器6采集多段背景噪声信号，采集的背景噪声信号参见图3。对图3的背景噪声信号进行分析，获得时域特征和频域特征。频域特征是经傅里叶变换获得噪声的频带，背景噪声能量集中在低频段的0-8500HZ内，为下面的声音信号处理提供信号滤波频带基准。时域特征是通过常规的自相关统计和直方图技术，可获取背景噪声信号均方差<7 时域特征。参见图4-6，现场进行切割不同材料和不同厚度工件的工作，且利用图I的声音信 号计算机采集和控制系统获取多组水射流现场声音信号。现场试验需要考虑两种情况一种情况是材料固定，厚度变化；另一种情况是材料变化，厚度变化。由于超高压水射流切割现场声音非常嘈杂，对提取纯粹的水射流切割反射声音产生了严重的影响，而只有纯粹的水射流切割反射声音才反映切割工件的特征。因此，提取纯粹的水射流切割反射声音/￠)特征时，先采用图I的声音信号计算机采集和控制系统获取超高压水射流切割现场
声音信号_7(0，再通过与背景噪声信号均方差^比较判断是否有纯粹的水射流切割反射
声音信号/￠)。切割现场声音信号频谱参见图5，高频段8500-15000HZ内为有用信号。对于优化速度控制，纯粹的水射流切割反射声音信号/￠)最主要的特征是切割时间长度*和信号能量￡。h-t, = 是水射流切割开始时刻，t,是水射流切割结束时刻。如何精确提取纯粹的水射流切割反射声音信号的特征切割时间长度^和信号能量s是关键，需要精确确定奇异跳变的两个时刻L ^ ,而切割现场声音信号_7(0受到噪声影响，难以精确提取特征量。采集的切割现场声音信号的声音数据；中含有较多噪声，具有非平稳性特征，传统的频域滤波难以准确恢复纯粹的水射流切割反射声音/￠)，本发明采用小波信号消噪，通过小波多尺度分析得到奇异跳变精确时刻，消除了噪声影响。小波信号消噪的主要步骤如下
I )选取片)的Z个离散序列:7 )，L = 2M ,为离散数据样点数，《是正整数。对.KU进行多尺度正交小波分解(设尺度为5)，获得各尺度的小波分解系数其中《 = 1,2,...5为分解尺度，先二 0,2，...2&1为分解样点数，Fajt为小波分解系数。II )计算小波分解后各个尺度上的分解能量值。对于正交小波，用小波系数表示
2"T-i^
信号能量，则各个尺度上信号分布的分解能量可表示为P, = ZPU，^ = 1-5 , p表示
A-O匕 a尺度的分解能量。需要注意的是水射流背景噪声的影响在较高尺度上，在较低尺度上信号的分解能量主要反映真实信号的变化。这样，建立尺度和分解能量分布图，并利用前面获
得背景噪声的均方差I和滤波频带基准，基于最小距离原则确定尺度阈值A = 3，令大于该尺度阈值的各个尺度的小波系数为零，即= O',s = jr,…,5 ;其他尺度分解小波系数不
变，二 =，这样得到新的小波系数矛u。III )用新的小波系数承#进行信号小波逆变换，得到纯粹的水射流切割反射声·音/ 的近似反射声音信号/4)。对近似反射声音信号^ 进行Hilbert变换，求出其振幅包络,它反映了近似反射声音信号J1(X)信号的整体变化趋势，便于识别奇异特征。获取近似反射声音信号/#)信号和振幅包络『(<)后，可以精确确定声音信号奇异跳变时刻《，G。利用特征值搜索程序搜索当振幅包络的绝对值大于设定的幅值阈值门限A时，即1,(01 >▲，获取对应的水射流切割开始时刻& ;当振幅包络,( 的绝对值小于设定的幅值阈值门限时，即| U| <A ,获取对应的水射流切割结束时刻I2。设定的幅值阈值门限为A = CG5，这是经验值，信号能量5=2^1/4)12 ;计算#个水射流反射声音切穿时间4。
计算水射流反射声音切穿时间t - 。将提取的将L,丑作为水射流反射声音特征。将切穿时间信号能量5*特征作为最小二乘支持向量机模型的输入量，同时考虑了现场工作环境参数的水射流水压P、光滑度八对优化速度的影响，将这四个物理量作为增广的输入向量，第t个向量为& = [mt],下标t为试验样本数目。以人工经验获得的最佳喷嘴横移速度为输出量，即％= [&],这里的&是经过大
量现场试验和试凑获得的。考虑训练样本集为#个，越多越好，为后续的最小二乘支持向量机的建模提供完备的输入输出信息。将输出量与增广的输入向量对应，建立#个非线性输入输出向量。获取非线性输入输出数据{% 4 } J = I..I后，求解满足下面的方程参数h
0Ir rI -] 「O _
~(I)
1O 4- T IZ Z
其中，Z= [Z1, Z2,…, f , I = [l,。《* =[珥，0；2，…, ]7" , CIfjfc 为权值变量，Ij 为偏
置。n= K(P^xj) , i,j = l2,-,N > [( ) = exp(-为非线性映射核函数，将输入数据映射到高维特征空间，参数A人工选择，C为正规化参数。由于
高压水射流系统的最佳喷嘴横移速度与高压水射流的切穿时间、信号能量、水压、光滑度等参数是强非线性映射关系，这里建立了简易的最小二乘支持向量机模型，即基于式(I)
获得^ 4建立的最小二乘支持向量机模型
权利要求
1.一种基于声音特征提取的超高压水射流速度优化控制方法，其特征是包括如下步骤(1)在没有切割工件时采集多段背景噪声信号，对背景噪声信号进行分析获得时域特征和频域特征，背景噪声能量集中在0-8500HZ滤波频带基准内，背景噪声信号均方差为 (2)获取#组水射流切割现场声音信号_fU)，选取的￡个离散序列/电)，L = 2M是正整数，对>电)进行五个尺度的正交小波分解，获得各尺度的小波分解系数,认,a = l,2,...5为分解尺度，k= 0,2,...2*-1为分解个数；根据小波分解系数的平方和计算各尺度上的分解能量，建立分解尺度和分解能量分布图，并利用所述背景噪声的均方差^和所述滤波频带基准，基于最小距离原则确定尺度阈值A = 3，令大于该尺度阈值的各尺度的小波系数为零，其他尺度分解小波系数不变，得到新的小波系数，对所述新的小波系数进行信号小波逆变换，得到纯粹的水射流切割反射声音的近似反射声音信号M),求出/4)的振幅包络妒(1);当振幅包络,(I)的绝对值大于设定的幅值阈值门限A =0.35时，获取对应的水射流切割开始时刻4 ;当振幅包络『的的绝对值小于设定的幅值阈值门限A = 时，获取对应的水射流切割结束时刻得到#个水射流反射声音切穿时间& = I1 - ^，计算信号能量￡ = Z U I2 ; (3)将切穿时间i、能量特征作为输入量、以人工经验获得的最佳喷嘴横移速度为输出量建立最小二乘支持向量机模型； (4)切割工件时，在线采集切割现场声音信号，提取其切穿时间和信号能量，将切穿时间和信号能量输入到所述最小二乘支持向量机模型，获得最佳喷嘴横移速度。
全文摘要
本发明公开一种基于声音特征提取的超高压水射流速度优化控制方法，在没有切割工件时采集多段背景噪声信号，分析获得时域特征和频域特征，对选取的N组水射流切割现场声音信号的个离散序列进行五个尺度的正交小波分解，得到近似反射声音信号，求出的振幅包络；获取对应的水射流切割开始时刻和水射流切割结束时刻，得到切穿时间，计算信号能量；将切穿时间、信号能量作为输入量、以人工经验获得的最佳喷嘴横移速度为输出量建立最小二乘支持向量机模型；切割工件时，在线采集切割现场声音信号，提取其切穿时间和信号能量输入到最小二乘支持向量机模型，获得最佳喷嘴横移速度；可提高水射流切割的工作效率及切割产品的品质。
文档编号G05B13/02GK102672626SQ201210168068
公开日2012年9月19日 申请日期2012年5月28日 优先权日2012年5月28日
发明者孙月平, 张军, 李雪峰, 赵德安, 陈波 申请人:江苏大学