關(guān)佳亮 陳 玲 劉益嘉 朱生根

(①北京工業(yè)大學(xué)機(jī)電學(xué)院，北京100124;②北京工研精機(jī)股份有限公司，北京101312)

KDP(磷酸二氫鉀KH2PO4)晶體是一種性能優(yōu)良的非線性光學(xué)元件。目前，大口徑KDP 晶體是唯一可用作慣性約束核聚變系統(tǒng)中高功率激光驅(qū)動(dòng)裝置器的Pockels 盒的光學(xué)晶體材料［1］。表面粗糙度、表面波紋度和透射波是大口徑KDP 晶體3 項(xiàng)最難達(dá)到的技術(shù)指標(biāo)［2］，且表面波紋度是影響因素最復(fù)雜、最難控制的指標(biāo)。利用SPDT 技術(shù)對(duì)KDP 晶體進(jìn)行加工后，在已加工表面上沿著圓弧刀具的切削軌跡會(huì)形成一定的刀痕，該刀痕結(jié)構(gòu)會(huì)給KDP 晶體的透射波前添加小尺度周期性擾動(dòng)，這種小尺度的相位調(diào)制是強(qiáng)激光非線性增長(zhǎng)的噪聲源，在高功率情況下容易形成自聚焦破壞，是必須抑制的制造誤差［3］。而且，KDP 晶體表面波紋度具有很大的擾動(dòng)性，既與加工工藝參數(shù)的選擇有關(guān)，也受到機(jī)床系統(tǒng)的加工狀態(tài)及環(huán)境因素的影響［4-5］。由于KDP 晶體單點(diǎn)金剛石飛刀切削加工的進(jìn)給量很小，對(duì)于大尺寸KDP 晶體來(lái)講，需耗費(fèi)很長(zhǎng)的加工時(shí)間。因此，對(duì)KDP 晶體加工過(guò)程中表面波紋度的狀態(tài)進(jìn)行在線監(jiān)測(cè)對(duì)指導(dǎo)實(shí)際加工具有重要意義，若波紋度可以滿足要求則繼續(xù)加工，反之，則快刀走過(guò)，調(diào)整機(jī)床狀態(tài)或工藝參數(shù)后再進(jìn)行加工，利于節(jié)省加工、卸工件和檢測(cè)的時(shí)間，減少刀具的磨損，提高加工效率和表面質(zhì)量。

1 聲發(fā)射技術(shù)簡(jiǎn)介與信號(hào)采集

聲發(fā)射(Acoustic Emission，簡(jiǎn)稱AE)技術(shù)［6］檢測(cè)的原理如圖1 所示，從聲發(fā)射源發(fā)出的彈性波經(jīng)介質(zhì)傳播將最終到達(dá)被撿測(cè)工件的表面，引起工件表面的機(jī)械振動(dòng)。安裝在工件表面的傳感器將探測(cè)到的瞬態(tài)位移轉(zhuǎn)換為電信號(hào)，通過(guò)聲發(fā)射采集系統(tǒng)將該電信號(hào)放大、處理、記錄和顯示，形成特定的參數(shù)和波形數(shù)據(jù)。再對(duì)數(shù)據(jù)進(jìn)行處理和分析，得到材料產(chǎn)生聲發(fā)射的特性。

采用AE 技術(shù)對(duì)KDP 晶體不同表面波紋度狀態(tài)的聲發(fā)射信號(hào)進(jìn)行采集與分析，得到AE 信號(hào)特征與表面波紋度的映射關(guān)系，從而通過(guò)AE 信號(hào)來(lái)判斷加工時(shí)表面波紋度的狀態(tài)［7］。

2 實(shí)驗(yàn)數(shù)據(jù)的采集

2.1 聲發(fā)射采集系統(tǒng)的設(shè)置

將數(shù)字濾波器設(shè)置為直通，波形和參數(shù)門限均設(shè)置為30db。在正式切削時(shí)，將環(huán)境噪聲和主軸空轉(zhuǎn)的信號(hào)初步濾除。在試驗(yàn)中涉及5 種機(jī)床轉(zhuǎn)速，分別為100 r/min、275 r/min、450 r/min、625 r/min 和800 r/min，根據(jù)SPDT 原理、刀具旋轉(zhuǎn)半徑和工件的大小進(jìn)行計(jì)算，得到在各轉(zhuǎn)速下刀具與工件的接觸時(shí)間分別為131.25 ms、48 ms、19.167 ms、21 ms 和16.4 ms。在時(shí)域和頻域上初步濾除環(huán)境噪聲和主軸空轉(zhuǎn)的聲發(fā)射信號(hào)后，各工藝參數(shù)組合下的聲發(fā)射信號(hào)表現(xiàn)為刀具轉(zhuǎn)動(dòng)一周的周期函數(shù)，但只有在刀具與工件接觸時(shí)聲發(fā)射系統(tǒng)才有信號(hào)產(chǎn)生，因此采樣的最小時(shí)間應(yīng)大于刀具與工件的接觸時(shí)間。根據(jù)奈奎斯特采樣定理，采樣頻率應(yīng)至少為信號(hào)最高頻率的兩倍。本試驗(yàn)的采樣頻率設(shè)置3 333 kHz，采樣長(zhǎng)度根據(jù)上述各轉(zhuǎn)速下刀具與工件的接觸時(shí)間進(jìn)行選擇。根據(jù)實(shí)際情況在SAEU2S 聲發(fā)射系統(tǒng)軟件中選擇參數(shù)間隔、鎖閉時(shí)間、峰值間隔和存儲(chǔ)路徑等，對(duì)KDP 晶體加工過(guò)程中不同工藝參數(shù)下的聲發(fā)射信號(hào)進(jìn)行采集與存儲(chǔ)。

2.2 實(shí)驗(yàn)數(shù)據(jù)的采集

采用SPDT 和二次通用回歸旋轉(zhuǎn)組合方法［8］進(jìn)行31 組試驗(yàn)，最后測(cè)量到每組試驗(yàn)的表面波紋度如表1所示。

表1 二次通用回歸旋轉(zhuǎn)組合實(shí)驗(yàn)設(shè)計(jì)與試驗(yàn)結(jié)果數(shù)據(jù)

3 聲發(fā)射信號(hào)的分析

聲發(fā)射系統(tǒng)采集到的信號(hào)是工件表面的瞬態(tài)位移，包含了切削、機(jī)床振動(dòng)和環(huán)境噪聲等因素，數(shù)據(jù)量大，監(jiān)測(cè)信號(hào)無(wú)法直接用來(lái)判斷KDP 晶體表面波紋度的狀態(tài)。故需對(duì)采集到的聲發(fā)射信號(hào)進(jìn)行處理，得到與表面波紋度對(duì)應(yīng)關(guān)系明顯的信號(hào)特征，實(shí)現(xiàn)對(duì)加工過(guò)程中KDP 晶體表面波紋度的在線監(jiān)測(cè)。

根據(jù)實(shí)際情況對(duì)表1 中的表面波紋度試驗(yàn)結(jié)果按從小到大的順序分為好、較好、中、較差和差5 個(gè)等級(jí)，如表2 所示。

表2 KDP 晶體表面波紋度試驗(yàn)結(jié)果分組

3.1 聲發(fā)射信號(hào)的時(shí)域分析

分析加工過(guò)程中聲發(fā)射信號(hào)時(shí)域波形，得到信號(hào)隨時(shí)間的變化情況，判斷出部分信號(hào)的特性，及整個(gè)時(shí)域信號(hào)隨波紋度增大時(shí)的整體變化趨勢(shì)。圖2 為表面波紋度等級(jí)1 至4 所對(duì)應(yīng)的單幀聲發(fā)射時(shí)域信號(hào)。

各單幀聲發(fā)射信號(hào)顯示了加工過(guò)程中主軸旋轉(zhuǎn)一周，刀具與工件接觸的全過(guò)程信號(hào)。由圖2 可得，每個(gè)波形圖中的信號(hào)在開(kāi)始和結(jié)束時(shí)幅值均較大，這是由刀具切入和切出KDP 晶體時(shí)的能量突變所致，可為后續(xù)的頻域分析提供判斷依據(jù)。忽略兩端的跳躍信號(hào)時(shí)，對(duì)比上述4 個(gè)圖形可知:隨著波紋度值的增大，聲發(fā)射信號(hào)中間段的幅值存在整體變大的趨勢(shì)。

為了提高聲發(fā)射信號(hào)分析的可靠度，對(duì)每組工藝參數(shù)下的信號(hào)選取10 個(gè)采樣長(zhǎng)度進(jìn)行分析，并取其統(tǒng)計(jì)特征的平均值作為聲發(fā)射信號(hào)的特征量。根據(jù)時(shí)域分析方法，對(duì)表2 中25 個(gè)波紋度值所對(duì)應(yīng)的聲發(fā)射信號(hào)進(jìn)行時(shí)域特征量的分析計(jì)算，得到統(tǒng)計(jì)結(jié)果如表3所示。將時(shí)域統(tǒng)計(jì)結(jié)果繪制成折線圖，如圖3 所示，其中均值取其絕對(duì)值。

分析表3 和圖3 可知，在KDP 晶體SPDT 加工過(guò)程中，AE 信號(hào)的時(shí)域特征均值μx在KDP 晶體波紋度值增大時(shí)，沒(méi)有呈現(xiàn)規(guī)律性變化;方差σx2 隨著波紋度值的增大有整體增大趨勢(shì)，但存在間斷性跳躍;而均方值ψx2在排除個(gè)別特殊情況時(shí)，與波紋度值大小有著明顯的對(duì)應(yīng)關(guān)系，故可作為加工過(guò)程中監(jiān)測(cè)KDP 晶體表面波紋度狀態(tài)的重要特征量。

3.2 聲發(fā)射信號(hào)的頻域分析

3.2.1 聲發(fā)射信號(hào)功率譜分析

在對(duì)時(shí)域信號(hào)進(jìn)行頻域分析時(shí)，聲發(fā)射信號(hào)中的直流分量會(huì)造成頻譜泄漏，使直流信號(hào)泄漏到其他頻譜分量尤其是低頻分量上，從而影響頻譜分析的準(zhǔn)確性。因此在頻域分析前，需先在時(shí)域內(nèi)去除直流分量。然后利用MATLAB 中的FFT 函數(shù)對(duì)時(shí)域信號(hào)進(jìn)行傅立葉變化得到頻譜圖，再對(duì)聲發(fā)射信號(hào)進(jìn)行功率譜分析，圖4 依次列出了代表波紋度等級(jí)1 至4 的聲發(fā)射信號(hào)功率譜圖。

表3 聲發(fā)射信號(hào)時(shí)域特征統(tǒng)計(jì)

由功率譜圖4 可知，聲發(fā)射信號(hào)主要分布在300 kHz 以下，在0 ～300 kHz 內(nèi)存在兩個(gè)信號(hào)集中段:第一個(gè)信號(hào)集中段的幅值較大，并基本穩(wěn)定在50 kHz 附近;第二個(gè)信號(hào)集中段的頻率隨著KDP 晶體表面波紋度值的變化而波動(dòng)，且其幅度隨著波紋度值的增大明顯增加。

3.2.2 聲發(fā)射信號(hào)頻域特征提取

對(duì)聲發(fā)射信號(hào)進(jìn)行濾波分析，提取不同頻段的信號(hào)特征量。選擇0 ～50 kHz、50 ～100 kHz、100 ～150 kHz、150 ～200 kHz、200 ～250 kHz 和250 ～300 kHz 這6 個(gè)頻段進(jìn)行濾波分析。如圖5 為波紋度等級(jí)為2 功率譜濾波結(jié)果。

由圖5 的AE 信號(hào)時(shí)域波形圖和功率譜可知，刀具切入、切出的信號(hào)主要分布在0 ～50 kHz 以內(nèi)。通過(guò)對(duì)功率譜值在各濾波頻段上進(jìn)行積分，得到的平均能量作為該頻段的信號(hào)特征。取各頻段平均能量的10 個(gè)采樣長(zhǎng)度的平均值進(jìn)行統(tǒng)計(jì)，得到的結(jié)果如表4所示。通過(guò)折線圖來(lái)表示，如圖6 所示。

從表4 和圖6 中聲發(fā)射信號(hào)各頻段平均能量與波紋度的關(guān)系可知:隨著波紋度值的增大，50 ～100 kHz、100 ～150 kHz 和150 ～200 kHz 三個(gè)頻段的平均能量呈現(xiàn)出整體增大趨勢(shì)，其中100 ～150 kHz 頻段的平均能量與表面波紋度值的對(duì)應(yīng)關(guān)系最明顯。其他3 個(gè)頻段的AE 信號(hào)平均能量折線圖雜亂無(wú)章，與表面波紋度值沒(méi)有對(duì)應(yīng)關(guān)系。綜上得到100 ～150 kHz 可作為判斷KDP 晶體單點(diǎn)金剛石飛刀切削加工過(guò)程中表面波紋度狀態(tài)的特征頻段。

表4 AE 信號(hào)各濾波頻段的平均能量

4 結(jié)語(yǔ)

(1)對(duì)不同表面波紋度狀態(tài)下聲發(fā)射信號(hào)進(jìn)行時(shí)域特征統(tǒng)計(jì)，分析均值μx、方差σx2、均方值ψx2與波紋度的關(guān)系，得到均方值ψx2可作為KDP晶體加工過(guò)程中表面波紋度狀態(tài)在線監(jiān)測(cè)的重要時(shí)域特征量。

(2)基于MATLAB 軟件估計(jì)聲發(fā)射信號(hào)的功率譜圖，并對(duì)聲發(fā)射信號(hào)進(jìn)行6 個(gè)頻段的濾波，通過(guò)統(tǒng)計(jì)得到不同頻段的平均能量。分析各同頻段平均能量隨波紋度值增大的變化情況，得到150 ～200 kHz 可作為在線監(jiān)測(cè)KDP 晶體表面波紋度狀態(tài)的特征頻段。

［1］蘇根博，曾金波，賀友平，等.大截面KDP 晶體在激光核聚變研究中的應(yīng)用［J］.硅酸鹽學(xué)報(bào)，1997，25(6):717 -719.

［2］楊福興.KDP 晶體超精密加工技術(shù)的研究［J］. 制造技術(shù)與機(jī)床，2003(9):63 -65.

［3］楊力.先進(jìn)光學(xué)制造技術(shù)［M］.北京:科學(xué)出版社，2001:235 -286.

［4］姜偉，陳明君，李明全. KDP 晶體飛刀銑削加工參數(shù)對(duì)表面質(zhì)量影響的實(shí)驗(yàn)研究［J］. 航空精密制造技術(shù)，2009，45(5):4 -7.

［5］馬恩才.KDP 晶體超精密加工表面力學(xué)性能的研究［D］.哈爾濱:哈爾濱工業(yè)大學(xué)，2006:1 -5.

［6］國(guó)防科技工業(yè)無(wú)損檢測(cè)人員資格鑒定與認(rèn)證培訓(xùn)檢測(cè)編審委員會(huì).聲發(fā)射檢測(cè)［M］.北京:機(jī)械工業(yè)出版社，2005:1 -8.

［7］李孟源.聲發(fā)射檢測(cè)及信號(hào)處理［M］.北京:科學(xué)出版社，2010.

［8］王景賀，王洪祥，畢聽(tīng)，等.KDP 晶體單點(diǎn)金剛石切削表面粗糙度預(yù)測(cè)及實(shí)驗(yàn)研究［J］.光學(xué)技術(shù)，2006，32(2):267 -269.