HEATS 发表于 2010-9-12 11:02:26

刀具磨损的小波检测

<H1><FONT size=2>1 引言</FONT></H1>
<DD>刀具磨损的检测方法通常可分为直接测量法和间接测量法两类。直接测量法就是直接测量刀具后刀面磨损带中间部分的平均磨损量;而间接测量法则是通过测量与刀具磨损有关的物理量(如切削力、声发射等),并根据它们与刀具磨损之间的对应关系来判别刀具磨损情况。
<DD>常用的一种间接测量法是选取切削分力的比值和该比值的变化作为判别刀具磨损的特征量,一旦这两个特征量超过某一阈值,即判定刀具已经磨损。但在应用该方法对刀具磨损进行实际检测时,由于振动和随机噪声的干扰,很容易出现误判。此外,由于刀具的正常磨损与异常磨损之间的界限具有不确定性,因此对判别阈值的预先设定比较困难,必须依靠经验和多次切削试验才能确定。另外该方法在不同的应用场合性能不够稳定。为此,本文提出一种新的间接测量方法,即利用小波变换对切削力信号进行分析,通过计算切削力信号的奇异性指数来实现对刀具磨损的检测。
<H1><FONT size=2>2 奇异性指数与刀具磨损的对应关系</FONT></H1>
<DD>若信号f(t)在某点出现突变或某阶导数不连续,则称信号在该点具有奇异性,一般用Lipschitz 指数(简称Lip指数)来描述信号f(t)的奇异性程度。设有非负整数n,n≤<FONT face=symbol>a</FONT>≤n+1,如果存在常数A>0以及n 次多项式p<SUB>n</SUB>(t),对于t∈(t<SUB>0</SUB>-<FONT face=symbol>d</FONT>,t<SUB>0</SUB>+<FONT face=symbol>d</FONT>)时下式成立:<BR>
<TABLE width=550>
<TBODY>
<TR>
<TD align=middle><FONT size=2>|f(t)-p<SUB>n</SUB>(t-t<SUB>0</SUB>)|≤A|t-t<SUB>0</SUB>|<FONT face=symbol>a</FONT></FONT></TD>
<TD width=10><FONT size=2>(1)</FONT></TD></TR></TBODY></TABLE>则称f(t)在t<SUB>0</SUB>处是Lip<FONT face=symbol>a</FONT>的。如果f(t)在t<SUB>0</SUB>的Lip指数<FONT face=symbol>a</FONT><1,则称t<SUB>0</SUB>为f(t)的奇异点。
<DD>Lip指数给出了信号f(t)在t<SUB>0</SUB>点可导性的精确信息。若n≥1,则f(t)在t<SUB>0</SUB>处n阶可导,第n阶导数在t<SUB>0</SUB>点具有奇异性,且p<SUB>n</SUB>(t-t<SUB>0</SUB>)为f(t)在t<SUB>0</SUB>处的Taylor 级数的前n+1 项;当n=0 时,p<SUB>n</SUB>(t-t<SUB>0</SUB>)=f(t<SUB>0</SUB>),则由(1)式可知阶跃信号在阶跃点处<FONT face=symbol>a</FONT>=0。若信号f(t)无限次可导,则该信号不具有奇异性。
<DD>在切削加工中,当刀具磨损时,切削力会随之增大,且切削力在不同磨损阶段的差异可反映刀具的实际磨损状态。如果切削力信号突然上升或下降,则预示着刀具可能折断。由此可见,当刀具磨损时,相应的切削力信号存在奇异性,且奇异性指数的大小与刀具磨损状态具有一定对应关系,即刀具的磨损状态决定着切削力信号变化的剧烈程度,而奇异性指数正是用于描述这种剧烈程度的一个指标。刀具磨损越严重,相应的切削力信号变化越剧烈,Lip指数也越小;当刀具折断时,切削力信号会发生突变,此时对应的Lip指数<FONT face=symbol>a</FONT>=0;Lip指数越接近于1,则说明切削力越平稳,刀具的磨损程度也越轻。由此可知,通过Lip指数对切削力信号变化程度的描述,可以真实反映刀具磨损状态。基于上述分析,本文将刀具磨损分为5种状态,每一种状态对应于某一范围的切削力信号Lip指数(见表1),这样就在刀具磨损与Lip指数<FONT face=symbol>a</FONT>之间建立了直接联系,实际检测时通过计算切削力信号的Lip指数,根据表1即可确定刀具的磨损状况。
<TABLE cellSpacing=0 borderColorDark=#ffffff cellPadding=0 align=center bgColor=#e5ebba borderColorLight=#006600 border=1>
<CAPTION><FONT size=2><STRONG>表1 Lip指数与刀具磨损状态的对应关系</STRONG></FONT></CAPTION>
<TBODY>
<TR align=middle>
<TH><FONT size=2>磨损状态</FONT></TH>
<TH><FONT size=2>折断或破损</FONT></TH>
<TH><FONT size=2>剧烈磨损</FONT></TH>
<TH><FONT size=2>磨损较严重</FONT></TH>
<TH><FONT size=2>轻微磨损</FONT></TH>
<TH><FONT size=2>正常磨损</FONT></TH></TR>
<TR align=middle>
<TD><FONT size=2>Lip指数</FONT></TD>
<TD><FONT size=2>0.0~0.2</FONT></TD>
<TD><FONT size=2>0.2~0.4</FONT></TD>
<TD><FONT size=2>0.4~0.6</FONT></TD>
<TD><FONT size=2>0.6~0.8</FONT></TD>
<TD><FONT size=2>0.8~1.0</FONT></TD></TR></TBODY></TABLE>
<H1><FONT size=2>3 信号奇异性的小波检测原理</FONT></H1>
<DD>小波变换是检测信号奇异性的有效工具。小波&Oslash;(t)在L<SUP>2</SUP>(R)中满足容许性条件∫<SUP>R</SUP>&Oslash;(t)dt=0,小波变换可定义为信号f与小波伸缩系的卷积,即<BR>
<TABLE width=550>
<TBODY>
<TR>
<TD align=middle>
<TABLE>
<TBODY>
<TR align=middle>
<TD rowSpan=3><FONT size=2>Wf(s,x)=f&Oslash;s(x)=</FONT></TD>
<TD><FONT size=2>1</FONT></TD>
<TD style="FONT-SIZE: 18pt" rowSpan=3><FONT size=2>∫</FONT></TD>
<TD rowSpan=2><FONT size=2></FONT></TD>
<TD rowSpan=3><FONT size=2>f(t)&Oslash;</FONT></TD>
<TD><FONT size=2>(x-t)</FONT></TD>
<TD rowSpan=3><FONT size=2>dt</FONT></TD></TR>
<TR bgColor=#006600 height=1>
<TD><FONT size=2></FONT></TD>
<TD><FONT size=2></FONT></TD></TR>
<TR align=middle>
<TD><FONT size=2>s</FONT></TD>
<TD><SUB><FONT size=2>R</FONT></SUB></TD>
<TD><FONT size=2>s</FONT></TD></TR></TBODY></TABLE></TD>
<TD width=10><FONT size=2>(2)</FONT></TD></TR></TBODY></TABLE>式中,s为小波的尺度参数。
<DD>对于小波的奇异性检测效果,有些文献给出了如下定理:若f(t)在t<SUB>0</SUB>点的Lip指数为<FONT face=symbol>a</FONT>,且0<<FONT face=symbol>a</FONT><1,在t<SUB>0</SUB>的邻域(t<SUB>0</SUB>-<FONT face=symbol>d</FONT>,t<SUB>0</SUB>+<FONT face=symbol>d</FONT>)除t<SUB>0</SUB>点外均可导,则t<SUB>0</SUB>点是f(t)的一个孤立奇点,且f(t)的小波变换满足<BR>
<TABLE width=550>
<TBODY>
<TR>
<TD align=middle><FONT size=2>|Wf(s,x)|≤As<SUP><FONT face=symbol>a</FONT></SUP></FONT></TD>
<TD width=10><FONT size=2>(3)</FONT></TD></TR></TBODY></TABLE>
<DD>当尺度s→0 时,Wf(s,x)在t<SUB>0</SUB>点趋于零的速度为s<SUP><FONT face=symbol>a</FONT></SUP>,小于Wf(s,x)在其它点趋于零的速度s,因此,当尺度s较小时,f(t)的小波变换在t<SUB>0</SUB>点取得极大值。
<DD>对于两个尺度s<SUB>i</SUB>,s<SUB>i+1</SUB>,由式(3)两端取对数并相减,可得<BR>
<TABLE width=550>
<TBODY>
<TR>
<TD align=middle>
<TABLE>
<TBODY>
<TR align=middle>
<TD rowSpan=3><FONT size=2><FONT face=symbol>a</FONT>=tan</FONT></TD>
<TD style="FONT-SIZE: 18pt" rowSpan=3><FONT size=2>|</FONT></TD>
<TD><FONT size=2>Wf(s<SUB>i+1</SUB>,x)</FONT></TD>
<TD style="FONT-SIZE: 18pt" rowSpan=3><FONT size=2>|</FONT></TD>
<TD rowSpan=3><FONT size=2>/lg(</FONT></TD>
<TD><FONT size=2>s<SUB>i+1</SUB></FONT></TD>
<TD rowSpan=3><FONT size=2>)</FONT></TD></TR>
<TR bgColor=#006600 height=1>
<TD><FONT size=2></FONT></TD>
<TD><FONT size=2></FONT></TD></TR>
<TR align=middle>
<TD><FONT size=2>Wf(s<SUB>i</SUB>,x)</FONT></TD>
<TD><FONT size=2>s<SUB>i</SUB></FONT></TD></TR></TBODY></TABLE></TD>
<TD width=10><FONT size=2>(4)</FONT></TD></TR></TBODY></TABLE>
<DD>式(4)表明,以lgs为横坐标、lg|Wf|为纵坐标的模极大值曲线上相邻两点的斜率就是信号在奇异点处的Lip 指数。<BR>
<TABLE align=right>
<TBODY>
<TR>
<TD align=middle><FONT size=2><IMG src="http://www.chmcw.com/upload/news/RCL/13220_d0vnuf2007798509.gif"><BR><B>图1 刀具磨损时的切削力信号及其小波分解结果</B></FONT></TD></TR></TBODY></TABLE>
<H1><FONT size=2>4 试验与分析</FONT></H1>
<DD>在CM6140 车床上进行刀具磨损检测试验。试验条件为:工件材料:45 钢(调质);刀具:YT15 硬质合金刀片(<FONT face=symbol>g</FONT><SUB>0</SUB>=25°,<FONT face=symbol>a</FONT><SUB>0</SUB>=5°,<FONT face=symbol>K</FONT><SUB>r</SUB>=45°,<FONT face=symbol>l</FONT><SUB>s</SUB>=0°);切削用量:v<SUB>c</SUB>=1.65m/ s,a<SUB>p</SUB>=3mm,f=0.2mm/r;测力仪:瑞士KISTLER公司9257A1型压电晶体测力仪,采用5006 型电荷放大器,采样频率为100Hz。
<DD>图1a所示为试验测得的刀具磨损时的一个切削力信号,由于背景噪声较大,因此切削力的变化曲线较平滑。采用小波变换则可将切削力的变化从背景噪声中检测出来。设N(x)是一个方差为<FONT face=symbol>s</FONT><SUP>2</SUP> 的白噪声,它实际上是一种Lip指数为-0.5-<FONT face=symbol>e</FONT>(<FONT face=symbol>e</FONT>>0)的随机分布,其小波变换的平均功率为E|Wf(s,x)|<SUP>2</SUP>=<FONT face=symbol>s</FONT><SUP>2</SUP>‖<FONT face=symbol>y</FONT>‖/ s,它随着尺度的增大而迅速减小,这也是小波变换的另一优点。
<DD>以Gauss小波为小波基,对图1a所示切削力信号在4个尺度上进行变换,得到各尺度上的信号分别如图1b、c、d、e所示。可以看出,当尺度较小时,模极大值的数目急剧增加,此时的模极大值主要由测量过程中的噪声所支配,因此难以据此分析切削力的真实变化情况。但随着尺度的增大,模极大值的数目迅速减少。当尺度s=16时(如图1e所示),噪声所产生的模极大值已基本消失,仅在采样点n=61,116,147,196,251,320,389,448处还存在模极大值,且这些采样点处的模极大值在各个尺度上均存在,因此可以确定这是由信号本身的奇异性所产生的。由图1e可见,切削力在n=61处已开始波动,然后在n=116,147处进一步加强,最终在n=196处发生剧烈变化;随后又在n=251处骤然下降,在n=320,389,448等处反复波动,这种变化反映了刀具的磨损过程。根据这些点在各尺度上的小波变换值,利用式(4)可计算出它们的Lip指数(见表2)。
<TABLE cellSpacing=0 borderColorDark=#ffffff cellPadding=0 align=center bgColor=#e5ebba borderColorLight=#006600 border=1>
<CAPTION><FONT size=2><STRONG>表2 部分模极大值点的Lip指数</STRONG></FONT></CAPTION>
<TBODY>
<TR align=middle>
<TH><FONT size=2>模极大值点</FONT></TH>
<TH><FONT size=2>n=61</FONT></TH>
<TH><FONT size=2>n=116</FONT></TH>
<TH><FONT size=2>n=147</FONT></TH>
<TH><FONT size=2>n=196</FONT></TH>
<TH><FONT size=2>n=251</FONT></TH></TR>
<TR align=middle>
<TD><FONT size=2>Lip指数</FONT></TD>
<TD><FONT size=2>0.7635</FONT></TD>
<TD><FONT size=2>0.6216</FONT></TD>
<TD><FONT size=2>0.6138</FONT></TD>
<TD><FONT size=2>0.5820</FONT></TD>
<TD><FONT size=2>0.5926</FONT></TD></TR></TBODY></TABLE>
<DD>根据表1所示Lip指数与刀具磨损状态的对应关系,由表2所示Lip指数可判断出刀具的磨损有一个积累过程,并可判定刀具发生了较严重的磨损。对该刀具后刀面磨损带的直接测量结果(磨损量为0.382mm)证明本试验的分析结果与实际情况相符,正是这种磨损导致了切削力的波动。
<H1><FONT size=2>5 结论</FONT></H1>
<DD>刀具磨损时的切削力信号存在奇异性,奇异性的大小与刀具磨损的剧烈程度有关。利用小波变换可检测出这种奇异性,并能计算出奇异性指数的大小,从而可实现对刀具磨损状态的正确判定。 </DD>
               
页: [1]
查看完整版本: 刀具磨损的小波检测

中国磨床技术论坛
论 坛 声 明 郑重声明:本论坛属技术交流,非盈利性论坛。本论坛言论纯属发表者个人意见,与“中国磨削技术论坛”立场无关。 涉及政治言论一律删除,请所有会员注意.论坛资源由会员从网上收集整理所得,版权属于原作者. 论坛所有资源是进行学习和科研测试之用,请在下载后24小时删除, 本站出于学习和科研的目的进行交流和讨论,如有侵犯原作者的版权, 请来信告知,我们将立即做出整改,并给予相应的答复,谢谢合作!

中国磨削网