CN105656449B - A kind of electric spark voltage across poles signal processing method based on kalman filter - Google Patents

A kind of electric spark voltage across poles signal processing method based on kalman filter Download PDF

Info

Publication number
CN105656449B
CN105656449B CN201510991604.1A CN201510991604A CN105656449B CN 105656449 B CN105656449 B CN 105656449B CN 201510991604 A CN201510991604 A CN 201510991604A CN 105656449 B CN105656449 B CN 105656449B
Authority
CN
China
Prior art keywords
mrow
noise
inter
kalman filter
electrode voltage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510991604.1A
Other languages
Chinese (zh)
Other versions
CN105656449A (en
Inventor
奚学程
赵万生
陈默
刘宏达
陈昊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiao Tong University
Original Assignee
Shanghai Jiao Tong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiao Tong University filed Critical Shanghai Jiao Tong University
Priority to CN201510991604.1A priority Critical patent/CN105656449B/en
Publication of CN105656449A publication Critical patent/CN105656449A/en
Application granted granted Critical
Publication of CN105656449B publication Critical patent/CN105656449B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0202Two or more dimensional filters; Filters for complex signals
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H2017/0072Theoretical filter design
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0202Two or more dimensional filters; Filters for complex signals
    • H03H2017/0205Kalman filters

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Electrical Discharge Machining, Electrochemical Machining, And Combined Machining (AREA)

Abstract

本发明涉及一种基于卡曼滤波器的电火花极间电压信号处理方法,假设极间电压是由白噪音通过一个线性滤波器产生,在这个假设的基础上,利用Yule‑Walker新息增益矩阵法对采集的极间电压信号进行处理,可以得出从白噪音到极间电压的一个传递函数,再由该传递函数转换成设计卡曼滤波器所需的状态方程。本发明能够大幅减少采集的极间电压信号中的噪音干扰,从而避免了由于噪音原因而影响电火花加工过程中的伺服运动,减少了一些不必要的伺服轴来回运动,提升电火花加工的效率。

The present invention relates to a kind of method of processing the electric spark inter-electrode voltage signal based on the Kalman filter. It is assumed that the inter-electrode voltage is generated by white noise through a linear filter. On the basis of this assumption, the Yule-Walker innovation gain matrix is used By processing the collected inter-electrode voltage signal, a transfer function from white noise to inter-electrode voltage can be obtained, and then the transfer function is converted into the state equation required for designing a Kalman filter. The invention can greatly reduce the noise interference in the collected inter-electrode voltage signal, thereby avoiding the influence of the noise on the servo movement in the electric discharge machining process, reducing some unnecessary back and forth movements of the servo shaft, and improving the efficiency of the electric discharge machining .

Description

一种基于卡曼滤波器的电火花极间电压信号处理方法A Kalman Filter-Based Processing Method for EDM Interpolar Voltage Signal

技术领域technical field

本发明涉及电火花过程控制,属于特种加工技术领域,具体是一种用于电火花极间电压检测的卡曼滤波器设计方法。The invention relates to electric spark process control, belongs to the technical field of special processing, in particular to a design method of a Kaman filter used for electric spark inter-electrode voltage detection.

背景技术Background technique

电火花加工是利用工件和电极之间的一系列火花放电,对工件材料进行去除的一种过程。电火花加工普遍被用于模具、航空航天、医疗器械等领域。电火花加工与传统铣削加工相比,有很多不同的地方。比如电火花加工的速度,与铣削加工中根据的给定的速度进给不同,是根据测量的极间状态决定加工方向向前还是向后,并且根据当前测量的极间电压和设定伺服电压的差值决定向前或向后的速度,因此电火花加工速度是由极间状态决定的。因为极间电压决定了伺服轴的运动方向和速度,所以对极间电压的检测与估计就十分重要。目前最常用的极间电压的检测方法是平均电压法和滑动平均电压法,就是对从过去几个周期到现在采样周期的极间电压进行滑动平均电压的计算,计算的结果跟预先设定伺服电压进行比较,如图1所示。EDM is a process of removing material from a workpiece by a series of spark discharges between the workpiece and the electrode. EDM is widely used in molds, aerospace, medical equipment and other fields. Compared with traditional milling, EDM has many differences. For example, the speed of electric discharge machining is different from the given speed feed in milling. It determines whether the machining direction is forward or backward according to the measured inter-electrode state, and according to the current measured inter-electrode voltage and the set servo voltage. The difference determines the forward or backward speed, so the speed of EDM is determined by the state between poles. Because the inter-electrode voltage determines the direction and speed of the servo axis, it is very important to detect and estimate the inter-electrode voltage. At present, the most commonly used inter-electrode voltage detection methods are the average voltage method and the sliding average voltage method, which is to calculate the sliding average voltage for the inter-electrode voltage from the past few cycles to the current sampling cycle. The calculated results are consistent with the preset servo voltages are compared as shown in Figure 1.

公开号为101791728A的中国发明专利申请《静电感应微细电火花加工的非接触极间电压检测方法及其电路设计》是提出了一种用于静电感应微细电火花加工的非接触极间电压检测方法。其所用的方法用于检测极间电压,但是总体上是通过硬件的方法实现的。而且此方法仍然无法解决极间电压信号中噪音干扰的问题。本发明为检测后续步骤,假设极间电压是由白噪音通过一个线性滤波器产生,在这个假设的基础上,建立卡曼滤波模型,对极间电压信号进行卡曼滤波处理,减小信号中的噪音干扰,因而其与本发明的内容是完全不相同的。The Chinese invention patent application with publication number 101791728A "Non-contact inter-electrode voltage detection method and circuit design for electrostatic induction micro-EDM" proposes a non-contact inter-electrode voltage detection method for electrostatic induction micro-EDM. The method used is used to detect the inter-electrode voltage, but it is generally realized by hardware. Moreover, this method still cannot solve the problem of noise interference in the inter-electrode voltage signal. In order to detect the follow-up steps, the present invention assumes that the inter-electrode voltage is generated by white noise through a linear filter. On the basis of this assumption, a Kalman filter model is established to perform Kalman filter processing on the inter-electrode voltage signal to reduce the signal noise interference, thus it is completely different from the content of the present invention.

由于电火花加工时产生的强烈电场和磁场干扰,采集的极间电压中包含有很强的噪音干扰,在这种情况下,采用卡曼滤波器是一个很好的选择。因为卡曼滤波器是一个最优的线性滤波器,对信号的估计误差可以认为是白噪音。卡曼滤波器的设计是基于系统的状态方程和信号的自协方差矩阵,有了这些信息设计才能够进行。但是对于电火花加工来说,状态方程和自协方差矩阵都必须从采集的极间电压中得出。Due to the strong electric field and magnetic field interference generated during EDM, the collected inter-electrode voltage contains strong noise interference. In this case, using a Kalman filter is a good choice. Because the Kalman filter is an optimal linear filter, the estimation error of the signal can be considered as white noise. The design of the Kalman filter is based on the state equation of the system and the autocovariance matrix of the signal, and the design can only be carried out with these information. But for EDM, both the state equation and the autocovariance matrix must be derived from the collected inter-electrode voltage.

发明内容Contents of the invention

为克服上述现有技术的不足,本发明提供一种基于卡曼滤波器的电火花极间电压信号处理方法,能够大幅减少采集的极间电压信号中的噪音干扰,避免由于噪音原因而影响电火花加工过程中的伺服运动,减少不必要的伺服轴来回运动,提升电火花加工的效率。In order to overcome the above-mentioned deficiencies in the prior art, the present invention provides a method for processing electric spark inter-electrode voltage signals based on a Kalman filter, which can greatly reduce the noise interference in the collected inter-electrode voltage signals, and avoid affecting the electric spark due to noise. The servo movement during the spark machining process reduces unnecessary back and forth movement of the servo axis and improves the efficiency of EDM.

本发明的技术解决方案如下:Technical solution of the present invention is as follows:

一种基于卡曼滤波器的电火花极间电压信号处理方法,该方法包含以下步骤:A method for processing electrical spark inter-electrode voltage signals based on a Kalman filter, the method comprising the following steps:

步骤一:极间电压信号的实时采集,然后传输到数控系统,用于后面传递函数、状态方程和稳态增益的计算;Step 1: Real-time acquisition of the inter-electrode voltage signal, and then transmitted to the numerical control system for subsequent calculation of transfer function, state equation and steady-state gain;

步骤二:假设极间电压由白噪音通过一个线性滤波器产生,建立白噪音与极间电压之间的传递函数,由该传递函数转换成规范可控型的状态方程;Step 2: Assuming that the inter-electrode voltage is generated by white noise through a linear filter, the transfer function between the white noise and the inter-electrode voltage is established, and the transfer function is converted into a normative controllable state equation;

步骤三:由状态方程及白噪音方差求出测量噪音和过程噪音之间的自协方差矩阵,和两者之间的互协方差矩阵,进而得到极间电压与估计电压之间的稳态增益;Step 3: Calculate the autocovariance matrix between the measurement noise and the process noise and the cross-covariance matrix between the two from the state equation and the white noise variance, and then obtain the steady-state gain between the inter-electrode voltage and the estimated voltage ;

步骤四:根据状态方程和协方差矩阵设计卡曼滤波器;Step 4: Design a Kalman filter according to the state equation and covariance matrix;

步骤五:利用卡曼滤波器模型对极间电压信号进行处理,滤波后的极间电压作为伺服控制器的反馈信号,能够大幅减少极间电压信号中的噪音,从而减少不必要的伺服轴运动,提升加工效率,为加工服务。不必要的伺服轴运动是指,由于噪音的干扰使得原本能正常加工的电压信号小于当前伺服运动设定的最小阀值而导致数控系统认为当前短路,产生伺服轴回退,或者超出最大阀值而导致数控系统认为当前开路,造成伺服轴过冲。Step 5: Use the Kalman filter model to process the inter-electrode voltage signal. The filtered inter-electrode voltage is used as the feedback signal of the servo controller, which can greatly reduce the noise in the inter-electrode voltage signal, thereby reducing unnecessary servo axis movement , Improve processing efficiency and serve processing. Unnecessary servo axis movement means that due to noise interference, the voltage signal that can be processed normally is less than the minimum threshold value set by the current servo movement, which causes the CNC system to think that the current short circuit causes the servo axis to retreat, or exceeds the maximum threshold value As a result, the CNC system thinks that the current is open, causing the servo axis to overshoot.

本发明原理如下:Principle of the present invention is as follows:

本发明假设极间电压是由白噪音通过一个线性滤波器产生,在这个假设的基础上,利用Yule-Walker自协方差法对采集的极间电压信号进行处理,可以得出从白噪音到极间电压的一个传递函数。由这个传递函数转化成状态方程。白噪音和极间电压之间的传递函数可以表示为The present invention assumes that the inter-electrode voltage is generated by white noise through a linear filter. On the basis of this assumption, the Yule-Walker autocovariance method is used to process the collected inter-electrode voltage signal, and it can be obtained from white noise to the pole-to-electrode voltage signal. A transfer function of the voltage between. This transfer function is transformed into an equation of state. The transfer function between white noise and inter-electrode voltage can be expressed as

这里q-1是时间移动操作符,e(k)是白噪音,ai,i=1…n是分母上的系数,b0是一个使稳态增益为1的常数,y(k)是极间电压。为了得到方程(1),采用Yule-Walker自协方差的方法来求解系数,具体如下:Here q -1 is the time shift operator, e(k) is white noise, a i , i=1...n is the coefficient on the denominator, b 0 is a constant that makes the steady-state gain 1, y(k) is voltage between poles. In order to obtain Equation (1), the Yule-Walker autocovariance method is used to solve the coefficients, as follows:

其中r(k)是时间偏移为k时的自协方差,k的取值范围是从0到M。方程(2)可以表示为更为简洁的的矩阵形式where r(k) is the autocovariance when the time offset is k, and the value of k ranges from 0 to M. Equation (2) can be expressed in a more compact matrix form

Rc=r (3)Rc=r (3)

其中R是自协方差矩阵,故可得到系数向量cwhere R is the autocovariance matrix, so the coefficient vector c can be obtained

c=R-1r (4)系数向量c和系数向量a的关系可以表示为c=R -1 r (4) The relationship between coefficient vector c and coefficient vector a can be expressed as

c=[c1,c2,…,cM]=[-a1,-a2,…,-aM] (5)c=[c 1 ,c 2 ,…,c M ]=[-a 1 ,-a 2 ,…,-a M ] (5)

在有了传递方程之后,希望对输入信号是白噪音的假设做一个测试。根据从白噪音到极间电压信号之间的专递函数,可以利用从白噪音到极间电压的传递函数的逆函数来检测传递函数的输入e(k)是否为白噪音。e(k)可以由Now that we have the transfer equation, we want to test the assumption that the input signal is white noise. According to the transfer function between white noise and inter-electrode voltage signal, the inverse function of the transfer function from white noise to inter-electrode voltage can be used To detect whether the input e(k) of the transfer function is white noise. e(k) can be given by

求得。白噪音可由以下条件来判断:Get it. White noise can be judged by the following conditions:

根据传递函数方程(1),采用标准可控形式,可以求出状态方程According to the transfer function equation (1), using the standard controllable form, the state equation can be obtained

y(k)=[(b1-b0a1)(b2-b0a2)…(bn-b0an)]x(k)+Hw(k)+v(k)y(k)=[(b 1 -b 0 a 1 )(b 2 -b 0 a 2 )…(b n -b 0 a n )]x(k)+Hw(k)+v(k)

其中w(k)是过程噪音,v(k)是测量噪音。方程(8)可以表示为更为简洁的矩阵形式:where w(k) is the process noise and v(k) is the measurement noise. Equation (8) can be expressed in a more compact matrix form:

x(k+1)=Ax(k)+Gw(k)x(k+1)=Ax(k)+Gw(k)

y(k)=Cx(k)+Hw(k)+v(k) (9)y(k)=Cx(k)+Hw(k)+v(k) (9)

其中, in,

C=[(b1-b0a1)(b2-b0a2)…(bn-b0an)],H=b0C=[(b 1 -b 0 a 1 )(b 2 -b 0 a 2 )...(b n -b 0 a n )], H=b 0 .

过程噪音w(k),测量噪音v(k)和白噪音e(k)可以分别表示为Process noise w(k), measurement noise v(k) and white noise e(k) can be expressed as

w(k)=e(k) (10)w(k)=e(k) (10)

v(k)=e(k) (11)v(k)=e(k) (11)

在卡曼滤波器的设计中,需要用到二个自协方差矩阵和一个互协方差矩阵,为了方便方差矩阵的推导,先定义复合过程噪音和复合测量噪音 In the design of the Kalman filter, two autocovariance matrices and a cross-covariance matrix are needed. In order to facilitate the derivation of the variance matrix, the compound process noise is first defined and composite measurement noise

有了复合过程噪音和复合测量噪音状态方程可以写为With compound process noise and composite measurement noise The state equation can be written as

根据状态方程(14),可以推导出复合过程噪音的自协方差矩阵According to the state equation (14), the autocovariance matrix of the composite process noise can be derived

复合测量噪音的自协方差矩阵可由下式得到:The autocovariance matrix of the composite measurement noise can be obtained by the following formula:

复合过程噪音和复合测量噪音之间的互协方差矩阵可以表示为:compound process noise and composite measurement noise The cross-covariance matrix between can be expressed as:

其中用到的三个方差矩阵Q(k),R(k)和N(k)分别定义为The three variance matrices Q(k), R(k) and N(k) used are defined as

Q(k)=E(w(k)wT(k)) (18)Q(k)=E(w(k)w T (k)) (18)

R(k)=E(v(k)vT(k)) (19)R(k)=E(v(k) vT (k)) (19)

N(k)=E(w(k)vT(k)) (20)N(k)=E(w(k) vT (k)) (20)

在有了状态方程,和三个方差矩阵的信息之后,就可以设计卡曼滤波器了。事后状态矢量可以表示为事前状态的一种修正。After having the state equation and the information of the three variance matrices, the Kalman filter can be designed. Post-event state vector ex ante state a correction of .

其中M(k)是新息增益矩阵,是预测误差。L(k)是卡曼增益矩阵。where M(k) is the innovation gain matrix, is the prediction error. L(k) is the Kaman gain matrix.

在计算中需要状态估计误差矩阵,可以表示为:The state estimation error matrix is required in the calculation, which can be expressed as:

在瞬间状态过程过去之后,P(k)就进入了稳态。为了减少实时计算量,在电火花加工时可以应用稳态P。稳态P可以有下式计算After the transient state process has passed, P(k) enters a steady state. In order to reduce the amount of real-time calculation, the steady state P can be applied during EDM. Steady state P can be calculated by the following formula

P(1|0)=pI (27)P(1|0)=pI (27)

其中,P(1|0)是P的初始值,I是单位矩阵,p是一个常量,可取为100。Among them, P(1|0) is the initial value of P, I is the identity matrix, and p is a constant, which can be 100.

在电火花加工过程中,测量的极间电压将经过卡曼滤波器的滤波处理。为了使处理后的滤波电压值的范围和采集的原始电压保持在同一个读数范围之内,必须求得从原始信号y(k)到滤波之后的信号之间的稳态增益。这样在求出之后,可以乘以一个稳态增益的倒数,从而保持的卡曼滤波器的输入和输出电压保持在相同的取值范围。During the EDM process, the measured inter-electrode voltage will be filtered by a Kalman filter. In order to keep the range of the processed filtered voltage value and the collected original voltage within the same reading range, it is necessary to obtain the signal from the original signal y(k) to the filtered signal The steady-state gain between. This is asking for Afterwards, it can be multiplied by an inverse of the steady-state gain, thus keeping the input and output voltages of the Kalman filter in the same range of values.

从y(k)到的稳态增益from y(k) to Steady state gain of

状态方程表达如下:The state equation is expressed as follows:

对方程(29)至(31)进行z变换,我们可以得到Applying z-transformation to equations (29) to (31), we can get

因此之间的传递函数为:therefore The transfer function between them is:

其中:因此把z=1代入之间的传递函数,就可以得到从y(k)到的稳态增益in: So substitute z=1 into The transfer function between, you can get from y(k) to Steady state gain of

与现有技术相比,原始极间电压信号图中有很强的噪音干扰,在这种情况下,并不能反映当前加工的实际状态,有可能由于噪音的干扰使得原本正常加工的电压信号小于当前伺服运动设定的最小阀值而导致数控系统认为当前加工短路,产生伺服轴短路回退,或者超出最大阀值而导致数控系统认为当前加工开路,造成伺服轴过冲,影响正常加工。本发明能够大幅减少采集的极间电压信号中的噪音干扰,从而避免了由于噪音原因而影响电火花加工过程中的伺服运动,减少了一些不必要的伺服轴来回运动,提升电火花加工的效率。Compared with the existing technology, there is strong noise interference in the original inter-electrode voltage signal diagram. In this case, it cannot reflect the actual state of the current processing. It is possible that the original normal processing voltage signal is less than The minimum threshold value set by the current servo motion causes the CNC system to think that the current processing is short-circuited, resulting in a short-circuit retreat of the servo axis, or exceeds the maximum threshold and causes the CNC system to consider the current processing to be open, causing the servo axis to overshoot and affect normal processing. The invention can greatly reduce the noise interference in the collected inter-electrode voltage signal, thereby avoiding the influence of the noise on the servo movement in the electric discharge machining process, reducing some unnecessary back and forth movements of the servo shaft, and improving the efficiency of the electric discharge machining .

附图说明Description of drawings

图1是本发明的间隙电压检测的伺服控制图。Fig. 1 is a servo control diagram of gap voltage detection in the present invention.

图2是本发明的实施平台示意图。Fig. 2 is a schematic diagram of the implementation platform of the present invention.

图3是本发明的一个具体实施例的极间电压信号采集图,其中图(a)是未经处理的原始图,图(b)为经过卡曼滤波。Fig. 3 is a graph of inter-electrode voltage signal acquisition in a specific embodiment of the present invention, wherein graph (a) is an unprocessed original graph, graph (b) is a Kalman filter.

图4是一个具体实施例中同样加工10mm深的孔,分别对极间电压信号采用不同处理方法的加工速率对比图。Fig. 4 is a comparison chart of processing rates of different processing methods for inter-electrode voltage signals in the same embodiment of processing a hole with a depth of 10 mm.

具体实施方式Detailed ways

本发明的具体实施例在上海汉霸机电有限公司生产的HE 70电火花成形机上进行,平台示意图如图2所示,整个机床的运动控制是通过PMAC运动控制卡来实现的。PC机设定的实时采样周期为2ms,PMAC运动控制卡的伺服周期为0.88ms,PC机与PMCA运动控制卡是通过PMAC运动控制卡上的DPRAM模块实现的。极间电压信号是通过与PMAC控制卡相连的ACC-28A卡来采集,并且通过DPRAM模块实时传输到PC机。The specific embodiment of the present invention is carried out on the HE 70 EDM machine produced by Shanghai Hanba Electromechanical Co., Ltd. The schematic diagram of the platform is shown in Figure 2. The motion control of the entire machine tool is realized by the PMAC motion control card. The real-time sampling period set by the PC is 2ms, and the servo period of the PMAC motion control card is 0.88ms. The PC and the PMCA motion control card are realized through the DPRAM module on the PMAC motion control card. The inter-electrode voltage signal is collected through the ACC-28A card connected to the PMAC control card, and is transmitted to the PC in real time through the DPRAM module.

实施例:Example:

首先假设极间电压是由白噪音通过一个线性滤波器产生,建立白噪音与极间电压的传递函数,并且求出式(1)中的相关参数,具体见表1:First assume that the inter-electrode voltage is generated by white noise through a linear filter, establish the transfer function of white noise and inter-electrode voltage, and obtain the relevant parameters in formula (1), see Table 1 for details:

表1间隙电压的动力学模型采样周期2msTable 1 The dynamic model sampling period of the gap voltage is 2ms

然后用基于以上模型的卡曼滤波器进行加工10mm深的小孔,与滑动平均滤波器进行对比实验,实验使用的是Ф10的棒状石墨电极,工件为工具钢。加工条件如下表2所示:Then use the Kalman filter based on the above model to process a small hole with a depth of 10mm, and conduct a comparative experiment with the moving average filter. The experiment uses a Ф10 rod-shaped graphite electrode, and the workpiece is tool steel. The processing conditions are as shown in Table 2 below:

表2小孔加工实验参数Table 2 Small hole processing experimental parameters

伺服电压servo voltage 峰值电流peak current 脉冲宽度Pulse Width 脉冲间隔pulse interval 抬刀时间Lifting time 抬刀周期Lifting cycle 49V49V 28A28A 60μs60μs 40μs40μs 1s1s 5s5s

所述的滑动平均滤波器定义为对PMAC采样周期采集到的极间电压信号进行如下处理:Described moving average filter is defined as following processing to the inter-electrode voltage signal that PMAC sampling period gathers:

其中,u(k)为当前采样周期的采样信号。Among them, u(k) is the sampling signal of the current sampling period.

具体实施例中卡曼滤波效果如图3所示,原始极间电压信号图中有很强的噪音干扰,在这种情况下,并不能反映当前加工的实际状态,有可能由于噪音的干扰使得电压信号小于当前伺服运动设定的最小阀值而产生不必要的短路回退,或者超出最大阀值而导致数控系统认为当前开路,造成过冲,影响正常加工。In the specific embodiment, the effect of the Kalman filter is shown in Figure 3. There is strong noise interference in the original inter-electrode voltage signal diagram. In this case, it cannot reflect the actual state of the current processing, and it may be caused by noise interference. If the voltage signal is less than the minimum threshold value set by the current servo motion, unnecessary short-circuit retreat will occur, or if it exceeds the maximum threshold value, the CNC system will think that the current circuit is open, causing overshoot and affecting normal processing.

具体加工结果的表3所示:Shown in table 3 of concrete processing result:

表3两种滤波方法加工结果比较Table 3 Comparison of processing results of two filtering methods

其中回退时间比定义如下:where the backoff time ratio is defined as follows:

其中,Tretract为回退时间,Ttotal为包含抬刀的总加工时间。Among them, T retract is the retracting time, and T total is the total processing time including lifting the tool.

由表3可以看出,无论是加工时间还是回退时间比,使用卡曼滤波得到的结果均好于滑动平均滤波,具体从1mm加工到10mm的速度对比图见图4。而且就加工效果来说,采用滑动平均滤波3次加工均产生积碳,而采用卡曼滤波3次加工均没产生积碳,也就是采用卡曼滤波的加工稳定性比采用滑动平均滤波要好。It can be seen from Table 3 that regardless of the processing time or the retraction time ratio, the results obtained by using the Kaman filter are better than the moving average filter. The specific speed comparison chart from 1mm to 10mm is shown in Figure 4. And as far as the processing effect is concerned, the three processings using the sliding average filter all produce carbon deposits, while the three processings using the Kaman filter do not produce carbon deposits, that is, the processing stability of the Kaman filter is better than that of the sliding average filter.

以上详细描述了本发明的较佳具体实施例。本领域的技术人员应当明白,在不脱离本发明范围的情况下,可以根据本发明的构思做出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过改变材料或者是特点情况可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。The preferred specific embodiments of the present invention have been described in detail above. Those skilled in the art should understand that many modifications and changes can be made according to the concept of the present invention without departing from the scope of the present invention. Therefore, all technical solutions that can be obtained by those skilled in the art based on the concept of the present invention by changing materials or characteristics on the basis of the prior art should be within the scope of protection defined by the claims.

Claims (2)

1. A method for processing an interelectrode voltage signal based on a Kalman filter is characterized by comprising the following steps:
the method comprises the following steps: acquiring an interelectrode voltage signal in real time and transmitting the interelectrode voltage signal to a numerical control system;
step two: assuming that the interelectrode voltage is generated by white noise through a linear filter, establishing a transfer function between the white noise and the interelectrode voltage, and converting the transfer function into a standard controllable equation of state;
the transfer function between the white noise and the interelectrode voltage is expressed as:
<mrow> <mover> <mi>G</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <msup> <mi>q</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>e</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>b</mi> <mn>0</mn> </msub> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>a</mi> <mn>1</mn> </msub> <msup> <mi>q</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>+</mo> <mo>...</mo> <mo>+</mo> <msub> <mi>a</mi> <mi>n</mi> </msub> <msup> <mi>q</mi> <mrow> <mo>-</mo> <mi>n</mi> </mrow> </msup> </mrow> </mfrac> </mrow>
in the formula, q-1Is a time-shift operator, e (k) is white noise, aiI is 1 … n, b is a coefficient on the denominator0Is a constant that gives a steady state gain of 1, and y (k) is the filtered white noise;
the transfer function adopts Yule-Walker autocovariance to solve the coefficient, and the formula is as follows:
wherein r (k) is the autocovariance at time offset k, k ranging from 0 to M;
step three: solving an autocovariance matrix between the measurement noise and the process noise and a cross covariance matrix between the measurement noise and the process noise by using a state equation and a white noise variance so as to obtain a steady gain between the interelectrode voltage and the estimated voltage;
step four: designing a Kalman filter according to a state equation and a covariance matrix;
step five: and processing the interelectrode voltage signal by using a Kalman filter model, wherein the filtered interelectrode voltage is used as a feedback signal of the servo controller.
2. The method for processing the electric spark interpolar voltage signal based on the kalman filter of claim 1, wherein the step one collects the interpolar voltage by an ACC-28A card and transmits the collected interpolar voltage to a numerical control system in real time by a DPRAM module of a PMAC.
CN201510991604.1A 2015-12-24 2015-12-24 A kind of electric spark voltage across poles signal processing method based on kalman filter Active CN105656449B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510991604.1A CN105656449B (en) 2015-12-24 2015-12-24 A kind of electric spark voltage across poles signal processing method based on kalman filter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510991604.1A CN105656449B (en) 2015-12-24 2015-12-24 A kind of electric spark voltage across poles signal processing method based on kalman filter

Publications (2)

Publication Number Publication Date
CN105656449A CN105656449A (en) 2016-06-08
CN105656449B true CN105656449B (en) 2018-05-08

Family

ID=56477927

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510991604.1A Active CN105656449B (en) 2015-12-24 2015-12-24 A kind of electric spark voltage across poles signal processing method based on kalman filter

Country Status (1)

Country Link
CN (1) CN105656449B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104729492A (en) * 2013-12-18 2015-06-24 广西大学 Optical fiber gyroscope signal processing method based on Kalman filtering
CN105092711A (en) * 2015-08-04 2015-11-25 哈尔滨工业大学 Steel rail crack acoustic emission signal detecting and denoising method based on Kalman filtering

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8045645B2 (en) * 2007-06-08 2011-10-25 Telefonaktiebolaget Lm Ericsson (Publ) Signal processor for estimating signal parameters using an approximated inverse matrix

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104729492A (en) * 2013-12-18 2015-06-24 广西大学 Optical fiber gyroscope signal processing method based on Kalman filtering
CN105092711A (en) * 2015-08-04 2015-11-25 哈尔滨工业大学 Steel rail crack acoustic emission signal detecting and denoising method based on Kalman filtering

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Evaluation of linear Kalman filter processing geodetic kinematic measurements";Sonja Bogatin等;《Measurement》;20080630;第41卷(第5期);第561-578页 *
"电火花加工放电状态的自适应滤波";周明等;《哈尔滨工程大学学报》;20151130;第36卷(第11期);第1522-1525页 *

Also Published As

Publication number Publication date
CN105656449A (en) 2016-06-08

Similar Documents

Publication Publication Date Title
CN100449933C (en) Servo motor drive controller
CN107037340B (en) A reciprocating wire electric discharge machine tool discharge position acquisition system
CN103302369B (en) The control device of the WEDM of the correction of machining path is carried out according to programmed instruction
CN103691559B (en) Fuzzy control method of high frequency and high voltage power supply for electrostatic dust removal
CN106513879B (en) A kind of spark discharge state recognition and detection method based on chaology
JP6227599B2 (en) Wire electrical discharge machine with constant distance between poles
JP2015042432A (en) Wire electric discharge machine including average discharge delay time calculation means
DE69010625T2 (en) JUMP CONTROL SYSTEM OF A SPARKING EDM MACHINE.
DE69025917T2 (en) Touch sensing device for positioning slidable elements
CN111331211A (en) On-line penetration detection method for EDM small hole machining
Huang et al. Online workpiece height estimation for reciprocated traveling wire EDM based on support vector machine
CN105656449B (en) A kind of electric spark voltage across poles signal processing method based on kalman filter
CN105159390B (en) Dynamic variable step Maximum power point tracing in photovoltaic system with interference rejection ability
Jiang et al. Adaptive control for small-hole EDM process with wavelet transform detecting method
DE102017110735B4 (en) Method for producing a metal component, in particular a blade component of a turbomachine
CN105930589B (en) Data processing method of multi-axis simultaneous EDM feed rate post-processor based on space mapping
DE112009005053T5 (en) Drahterodierbearbeitungs device
DE3808646C1 (en)
CN110018669A (en) The profile errors control method of five-axle number control machine tool decoupling
CN209787027U (en) Feedforward sampling circuit capable of effectively suppressing ripple voltage and switching power supply thereof
CN106015692B (en) The adaptive proportional valve amplifier of card insert type vibrating signal
JPS597523A (en) Wire-cut electric discharge machine
Jaspreet et al. Optimization of machining parameters of Electric Discharge Machining for 202 stainless steel
DE102020007192A1 (en) CONTROL DEVICE AND PROGRAM
CN105665851A (en) Nanosecond-scale pulse width pulse power supply used for electrochemical processing

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant