• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    單顆金剛石切削無(wú)氧銅的聲發(fā)射關(guān)聯(lián)維特征?

    2020-09-24 03:09:04吳海勇
    應(yīng)用聲學(xué) 2020年4期
    關(guān)鍵詞:相空間切削速度維數(shù)

    吳海勇

    (1 漳州職業(yè)技術(shù)學(xué)院機(jī)械工程學(xué)院 漳州 363000)

    (2 華僑大學(xué)制造工程研究院 廈門(mén) 361021)

    (3 福建龍溪軸承(集團(tuán))股份有限公司博士后科研工作站 漳州 363000)

    0 引言

    無(wú)氧銅是一種高純度、高電導(dǎo)性和高延展性的重要工業(yè)金屬材料,在真空電子器件、電纜電線、電機(jī)制造等電子電力領(lǐng)域有著廣泛應(yīng)用。利用金剛石工具可實(shí)現(xiàn)無(wú)氧銅的精密和超精密加工[1?2],這主要是因?yàn)榍邢鳠o(wú)氧銅過(guò)程中,金剛石工具的機(jī)械磨損量甚微[3?5],并且金剛石與無(wú)氧銅的親和性差,不產(chǎn)生化學(xué)磨損[6?7],金剛石工具的磨損基本上不會(huì)影響到無(wú)氧銅的表面加工質(zhì)量。

    聲發(fā)射技術(shù)是無(wú)氧銅切削加工特征監(jiān)測(cè)和評(píng)價(jià)的主要方法之一[8?9],可用于評(píng)價(jià)切削加工過(guò)程中工件塑性變形和裂紋變化等特性[10?12],通過(guò)對(duì)切削加工過(guò)程中聲發(fā)射信號(hào)的特征參數(shù)、頻譜和功率譜等特征參量的分析和對(duì)比,可實(shí)現(xiàn)對(duì)無(wú)氧銅切削加工過(guò)程特征的監(jiān)測(cè)和評(píng)價(jià)[13]。

    根據(jù)文獻(xiàn)顯示,無(wú)氧銅切削過(guò)程的聲發(fā)射時(shí)域信號(hào)屬于非線性一維時(shí)間序列函數(shù)[5,7?8,14],具有標(biāo)度不變性的分形結(jié)構(gòu)特征,可利用混沌分形理論對(duì)無(wú)氧銅切削過(guò)程的聲發(fā)射信號(hào)進(jìn)行分析,其中分形維數(shù)法是混沌分形理論中的一種重要分析方法,可以應(yīng)用于解析混沌吸引子自相似性特征和定量描述非線性系統(tǒng)動(dòng)力學(xué)行為。在分形維數(shù)法中,基于Grassberger-Procaccia(G-P)算法的關(guān)聯(lián)維數(shù)是用于計(jì)算非線性時(shí)間序列數(shù)據(jù)的一種常見(jiàn)重要方法,能較好反映出動(dòng)力學(xué)系統(tǒng)分形結(jié)構(gòu)的不規(guī)則度和復(fù)雜程度,在監(jiān)測(cè)金屬加工形貌特征描述[15?16]及聲發(fā)射信號(hào)分析[17?18]等領(lǐng)域有著廣泛應(yīng)用。然而,無(wú)氧銅在不同切削參數(shù)工況下的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)特征尚未有相關(guān)報(bào)道。

    基于上述分析,本文跟蹤監(jiān)測(cè)單顆金剛石旋轉(zhuǎn)切削無(wú)氧銅過(guò)程的聲發(fā)射信號(hào),分析不同切削工況下的聲發(fā)射時(shí)域信號(hào),采用G-P算法分析計(jì)算聲發(fā)射信號(hào)的關(guān)聯(lián)維數(shù)特征,旨在為無(wú)氧銅切削機(jī)理提供理論參考。

    1 試驗(yàn)條件及方法

    試驗(yàn)在Planomat HP-480 高速精密平面磨床上開(kāi)展,見(jiàn)圖1。單顆金剛石工具通過(guò)螺紋連接安裝固定在鋁盤(pán)基體(?360 mm×50 mm)上,主軸帶動(dòng)鋁盤(pán)基體和單顆金剛石工具以一定切削速度旋轉(zhuǎn)切削無(wú)氧銅,無(wú)氧銅通過(guò)夾具安裝固定在磁力吸盤(pán)工作臺(tái)上,并隨工作臺(tái)作橫向進(jìn)給運(yùn)動(dòng)。聲發(fā)射傳感器安裝在夾具側(cè)邊,通過(guò)信號(hào)放大器和采集卡收集、提取和處理分析聲發(fā)射信號(hào)。

    采用Element Six 公司SDB1125 磨料級(jí)金剛石,磨粒平均粒徑為0.725 mm;利用NiCrBSi 釬料制備釬焊單顆金剛石磨粒工具。試驗(yàn)所采用的無(wú)氧銅物理機(jī)械性能如表1所示。

    圖1 單顆金剛石切削無(wú)氧銅試驗(yàn)裝置圖Fig.1 The experimental setup diagram for single diamond cutting oxygen free copper

    表1 無(wú)氧銅的物理機(jī)械性能Table1 Mechanical and physical performance of oxygen free copper

    表2 聲發(fā)射信號(hào)采集參數(shù)設(shè)置Table2 The parameter settings of AE signals acquisition

    采用的聲發(fā)射檢測(cè)系統(tǒng)是Physical Acoustics Corporation 公司的MICRO-II2型,利用PCI-2-PAC 四通道信號(hào)采集系統(tǒng)采集無(wú)氧銅切削過(guò)程的聲發(fā)射信號(hào),試驗(yàn)前,先進(jìn)行斷鉛筆芯試驗(yàn),以保證聲發(fā)射傳感器安裝位置和接觸狀態(tài)的合理性,并通過(guò)空載試驗(yàn)測(cè)定聲發(fā)射信號(hào)門(mén)檻電壓等主要采集參數(shù)值,主要參數(shù)設(shè)置如表2所示。

    試驗(yàn)采用單因素變量法監(jiān)測(cè)不同切削參數(shù)下的聲發(fā)射信號(hào),變量因素有3個(gè),分別是切削速度vs、切削深度ap和進(jìn)給速度vw,每個(gè)變量因素有4個(gè)參數(shù)值,變量因素的具體參數(shù)值如表3所示。表中帶下劃線的參數(shù)值表示在試驗(yàn)該變量因素時(shí)其余兩個(gè)變量因素的取值,即當(dāng)切削深度為試驗(yàn)變量時(shí),切削速度和進(jìn)給速度取固定值分別為20 m/s和10 m/min,其余試驗(yàn)參數(shù)值以此類(lèi)推。

    表3 切削試驗(yàn)參數(shù)表Table3 Experimental parameters of cutting tests

    2 試驗(yàn)結(jié)果

    2.1 劃痕形貌

    單顆金剛石磨粒切削無(wú)氧銅的劃痕形貌(以vs=20 m/s、vw=10 m/min、ap=10 μm為例)見(jiàn)圖2,金剛石旋轉(zhuǎn)劃擦切削無(wú)氧銅的劃痕呈現(xiàn)出“滑擦-耕犁-切削-耕犁-滑擦”形貌。在金剛石切入無(wú)氧銅時(shí),金剛石切刃剛接觸無(wú)氧銅,切削深度由零逐漸增大,金剛石并未實(shí)際參與切削無(wú)氧銅,在無(wú)氧銅表面滑動(dòng)摩擦并使其產(chǎn)生彈性變形,形成滑擦劃痕。當(dāng)金剛石繼續(xù)切入工件,無(wú)氧銅表面產(chǎn)生塑性變形,無(wú)氧銅受到金剛石切刃的擠壓向兩邊塑性流動(dòng),在金剛石兩側(cè)形成微微隆起的耕犁溝槽。隨著切削深度的不斷增大,金剛石繼續(xù)擠壓無(wú)氧銅,使無(wú)氧銅在金剛石前部產(chǎn)生剪切滑移形成切屑,并在金剛石兩側(cè)形成工件材料的塑性隆起。在達(dá)到最大切深之后,金剛石切出無(wú)氧銅,切出過(guò)程劃痕的形成與切入過(guò)程相反。

    圖2 無(wú)氧銅表面的劃痕形貌Fig.2 The scratching morphology of oxygen free copper surface

    2.2 聲發(fā)射特征參數(shù)

    不同切削參數(shù)對(duì)聲發(fā)射特征參數(shù)值的影響如表4所示。從表中可知,隨著切削速度的增加,上升時(shí)間、均方根值(Root mean square,RMS)和能量基本呈上升趨勢(shì),而計(jì)數(shù)、持續(xù)時(shí)間和幅值的變化規(guī)律并不明顯。隨著切削深度的增加,上升時(shí)間、能量、持續(xù)時(shí)間、幅值和RMS值均呈現(xiàn)出顯著增加的趨勢(shì),而計(jì)數(shù)值的變化較為不明顯,可見(jiàn),切深對(duì)聲發(fā)射特征參數(shù)值的影響是較為顯著的。隨著進(jìn)給速度的增加,各個(gè)特征參數(shù)值均呈現(xiàn)出在一定范圍內(nèi)波動(dòng)變化,但并未呈現(xiàn)出明顯的變化規(guī)律。由此可見(jiàn),切削深度和切削速度對(duì)聲發(fā)射特征參數(shù)的影響相對(duì)較為顯著,而進(jìn)給速度對(duì)聲發(fā)射特征參數(shù)的影響較不顯著。

    2.3 切削速度對(duì)聲發(fā)射信號(hào)的影響

    切削速度對(duì)聲發(fā)射時(shí)域信號(hào)的影響見(jiàn)圖3。金剛石切削無(wú)氧銅的聲發(fā)射時(shí)域信號(hào)振動(dòng)幅值變化范圍在0~0.2 V之間,切削速度由小增大時(shí),其最大振幅分別達(dá)到了0.2 V、0.12 V、0.11 V和0.2 V。另外,可以發(fā)現(xiàn),不同切削速度時(shí),聲發(fā)射時(shí)域信號(hào)包絡(luò)線的疏密程度是有所不一樣的,切削速度為10 m/s(圖3(a))和25 m/s(圖3(d))時(shí),信號(hào)包絡(luò)線較為密集,且其信號(hào)振幅值亦較大;切削速度為15 m/s(圖3(b))和20 m/s(圖3(c))的信號(hào)包絡(luò)曲線較為稀疏,振幅值也相對(duì)較小。

    表4 切削參數(shù)對(duì)聲發(fā)射特征參數(shù)值的影響Table4 The influence of AE characteristic parameters on cutting parameters

    圖3 不同切削速度時(shí)的聲發(fā)射時(shí)域信號(hào)Fig.3 The AE time-domain signals with different cutting speeds

    2.4 切削深度對(duì)聲發(fā)射信號(hào)的影響

    金剛石以不同切深切削無(wú)氧銅的聲發(fā)射時(shí)域信號(hào)如圖4所示。從圖中可見(jiàn),切削深度的變化對(duì)聲發(fā)射時(shí)域信號(hào)有著顯著的影響,切深與聲發(fā)射信號(hào)振幅呈現(xiàn)出顯著的正效應(yīng)關(guān)系,即隨著切深的增加,信號(hào)振蕩幅值亦隨之逐漸增大,最大振幅值分別達(dá)到了0.02 V、0.04 V、0.08 V和0.21 V;聲發(fā)射時(shí)域信號(hào)包絡(luò)曲線密集度和包絡(luò)面積亦隨著切深的增加而顯著增加。

    圖4 不同切深時(shí)的聲發(fā)射時(shí)域信號(hào)Fig.4 The AE time-domain signals with different cutting depths

    2.5 進(jìn)給速度對(duì)聲發(fā)射信號(hào)的影響

    進(jìn)給速度對(duì)聲發(fā)射時(shí)域信號(hào)的影響如圖5所示??梢?jiàn),進(jìn)給速度由小增大時(shí),聲發(fā)射時(shí)域信號(hào)最大振幅值分別達(dá)到了0.1 V、0.09 V、0.1 V和0.12 V;時(shí)域信號(hào)包絡(luò)曲線的疏密度和包絡(luò)面積并未隨著進(jìn)給速度的增加而呈現(xiàn)出明顯變化規(guī)律。由此可見(jiàn),進(jìn)給速度對(duì)聲發(fā)射時(shí)域信號(hào)的影響并不顯著。

    圖5 不同進(jìn)給速度時(shí)的聲發(fā)射時(shí)域信號(hào)Fig.5 The AE time-domain signals with different feeding speeds

    3 關(guān)聯(lián)維分析

    3.1 G-P算法

    混沌系統(tǒng)是一個(gè)復(fù)雜的動(dòng)力學(xué)系統(tǒng),可用豪斯道夫維數(shù)、信息維數(shù)和相似維數(shù)等方法對(duì)其進(jìn)行解析,但是這些方法所描述的對(duì)象系統(tǒng)往往屬于簡(jiǎn)單規(guī)則的分析系統(tǒng)。從上述試驗(yàn)結(jié)果可見(jiàn),無(wú)氧銅切削聲發(fā)射信號(hào)是一個(gè)復(fù)雜的混沌動(dòng)力學(xué)系統(tǒng),因此,本文采用基于經(jīng)典G-P算法的關(guān)聯(lián)維數(shù)法分析無(wú)氧銅切削聲發(fā)射信號(hào)。G-P算法是利用嵌入理論和相空間重構(gòu)思想,從一維非線性時(shí)間序列動(dòng)力系統(tǒng)中直接計(jì)算出關(guān)聯(lián)維數(shù)的計(jì)算方法[19]。金剛石切削無(wú)氧銅的聲發(fā)射信號(hào)屬于一維非線性時(shí)間序列{xn,n=1,2,3,···,N},需對(duì)信號(hào)進(jìn)行相空間的重構(gòu),其中時(shí)間延遲法是一種行之有效的計(jì)算方法,利用該方法可將金剛石切削無(wú)氧銅聲發(fā)射信號(hào)重構(gòu)在長(zhǎng)度為N、維數(shù)為m的歐氏空間Rm中,并使該空間中的數(shù)據(jù)點(diǎn)維持原聲發(fā)射信號(hào)中未知吸引子的拓?fù)涮匦裕?jīng)時(shí)間延遲法處理后,可得到相空間向量Yi:

    其中,?t為切削聲發(fā)射信號(hào)的采樣時(shí)間間隔,i=1,2,3,···,N?m。則可得到N?m個(gè)相空間向量序列

    若相空間軌跡矩陣Y可描述為

    金剛石切削無(wú)氧銅的聲發(fā)射信號(hào)是非線性變化的,切削參數(shù)的變化影響了聲發(fā)射信號(hào)相空間狀態(tài)向量的重構(gòu)。試驗(yàn)采樣頻率為10 MHz,單次旋轉(zhuǎn)切削采樣點(diǎn)數(shù)達(dá)到15360個(gè),滿足聲發(fā)射信號(hào)相空間重構(gòu)對(duì)一維時(shí)間序列點(diǎn)數(shù)量的要求。在確定合適的時(shí)延參數(shù)和嵌入維數(shù)之后,即可得到無(wú)氧銅切削聲發(fā)射信號(hào)的關(guān)聯(lián)函數(shù)C(r):

    式(2)中,H(x)為階躍函數(shù),r為相空間的測(cè)度距離,rij為相空間中Yi到Y(jié)j之間的測(cè)度距離。對(duì)于點(diǎn)Yi而言,共有M?1個(gè)測(cè)度距離,關(guān)聯(lián)函數(shù)是統(tǒng)計(jì)無(wú)氧銅切削聲發(fā)射信號(hào)相空間中距離小于或是等于測(cè)度值r的數(shù)據(jù)點(diǎn)所占的百分比,若r值過(guò)大,C(r)趨于1;若r值過(guò)小,C(r)趨于0。在取合適測(cè)度值時(shí),C(r)=rυ,其中υ可認(rèn)為是對(duì)關(guān)聯(lián)維維數(shù)D的一種逼近計(jì)算,無(wú)氧銅切削聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)D的計(jì)算可如式(3)所示:

    因此,重構(gòu)無(wú)氧銅切削聲發(fā)射時(shí)域信號(hào)的相空間,計(jì)算出其關(guān)聯(lián)函數(shù)與時(shí)延參數(shù)的雙對(duì)數(shù)曲線關(guān)系(即lnC(τ)-lnτ),分析lnC(τ)-lnτ曲線中無(wú)標(biāo)度區(qū)間的斜率,即可得到無(wú)氧銅切削聲發(fā)射信號(hào)的關(guān)聯(lián)維數(shù)D。

    3.2 時(shí)間延遲的確定

    時(shí)延參數(shù)τ的選擇直接影響到關(guān)聯(lián)維數(shù)計(jì)算的準(zhǔn)確性,在確定合適時(shí)間延遲τ時(shí)應(yīng)充分考慮到以下因素:(1)若時(shí)間延遲過(guò)小,則任意兩個(gè)相鄰狀態(tài)矢量在數(shù)值上太過(guò)于接近,聲發(fā)射信號(hào)相空間軌跡趨同,致使相空間信息過(guò)于冗余;(2)若時(shí)間延遲過(guò)大,又將使相鄰兩個(gè)狀態(tài)矢量距離過(guò)大,將導(dǎo)致向量之間喪失關(guān)聯(lián)性;(3)時(shí)間延遲應(yīng)是試驗(yàn)采樣時(shí)間間隔的倍數(shù)。

    通過(guò)計(jì)算無(wú)氧銅切削聲發(fā)射信號(hào)時(shí)間序列的自相關(guān)函數(shù)得到相空間時(shí)間延遲參數(shù)。自相關(guān)函數(shù)顯示了無(wú)氧銅切削聲發(fā)射信號(hào)及其時(shí)間延遲之間的相似程度,并合理計(jì)算出兩者之間從不相關(guān)到冗余的折中度量。確定無(wú)氧銅切削聲發(fā)射信號(hào)相空間時(shí)間延遲的自相關(guān)函數(shù)可如式(4)所示:

    其中,x為聲發(fā)射信號(hào)數(shù)據(jù)序列的均值,σ為標(biāo)準(zhǔn)差。自相關(guān)函數(shù)曲線首次達(dá)到值1/e時(shí)的τ確定為最優(yōu)時(shí)延值[20?21]。利用Matlab 軟件,可以計(jì)算得到聲發(fā)射信號(hào)相空間的R(τ)-τ之間的關(guān)系,如圖6所示。從圖中可見(jiàn),不同切削參數(shù)時(shí)的自相關(guān)函數(shù)R(τ)與時(shí)延參數(shù)τ之間的變化特征是各異的,且各個(gè)時(shí)間延遲值亦各不相同。切削參數(shù)由小增大時(shí),不同切削速度下時(shí)間延遲的最佳取值分別為23、72、47和36;不同切削深度下時(shí)間延遲的最佳取值分別為21、27、62和48;不同進(jìn)給速度下時(shí)間延遲的最佳取值分別為47、23、38和26。

    圖6 不同工況下的R(τ)-τ的變化關(guān)系Fig.6 The relationship of R(τ)-τ for AE signals with different cutting conditions

    3.3 切削參數(shù)對(duì)關(guān)聯(lián)維數(shù)的影響分析

    結(jié)合式(2)可知,關(guān)聯(lián)積分C(r)的取值與測(cè)度r息息相關(guān),在選擇測(cè)度r時(shí),若取值過(guò)小,則C(r)=0,系統(tǒng)噪聲信號(hào)將干擾聲發(fā)射有用信號(hào);若取值過(guò)大,則C(r)=1,導(dǎo)致聲發(fā)射有用信號(hào)被淹沒(méi)。因此,測(cè)度r的取值應(yīng)涵蓋整個(gè)相空間的所有無(wú)標(biāo)度區(qū)間。本文采用等指數(shù)比例遞增法計(jì)算聲發(fā)射信號(hào)相空間關(guān)聯(lián)函數(shù)。利用3.2節(jié)所確定的時(shí)延參數(shù)值,可進(jìn)一步計(jì)算得到不同切削參數(shù)下的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)D。根據(jù)相關(guān)文獻(xiàn)可知[19,21],通過(guò)關(guān)聯(lián)維數(shù)值可以進(jìn)一步確定信號(hào)系統(tǒng)的動(dòng)力學(xué)特性,具體如表5所示。

    表5 關(guān)聯(lián)維數(shù)與系統(tǒng)動(dòng)力學(xué)特性Table5 The relationship of correlation dimension and system dynamic characteristic

    因此,根據(jù)計(jì)算得到無(wú)氧銅在不同切削參數(shù)下的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)D,利用表5可進(jìn)一步確定該聲發(fā)射信號(hào)系統(tǒng)屬于何種運(yùn)動(dòng)狀態(tài)特性。

    3.3.1 切削速度對(duì)關(guān)聯(lián)維數(shù)的影響分析

    切削速度對(duì)聲發(fā)射信號(hào)雙對(duì)數(shù)曲線數(shù)的影響如圖7所示。隨著嵌入維數(shù)的增大,雙對(duì)數(shù)曲線的lnC(r)取值逐漸增大,不同切削速度的雙對(duì)數(shù)曲線變化特征較為相似,不同切削速度的雙對(duì)數(shù)曲線均呈階段性增長(zhǎng)趨勢(shì),并逐漸收斂于飽和狀態(tài)。當(dāng)m較小時(shí),曲線斜率較小;隨著m的增大,曲線斜率亦隨之逐漸增大,分形區(qū)域無(wú)標(biāo)度區(qū)間不斷縮小并趨于飽和;在D值達(dá)到飽和之后,隨著m的增加,lnC(r)值基本不再發(fā)生變化。另外,在切削速度較小時(shí)(vs=10 m/s,圖7(a)),lnC(r)的取值范圍在?1.1~0之間,隨著切削速度的增大,其他3個(gè)切削速度(圖7(b)、7(c)和7(d))的lnC(r)的取值范圍則擴(kuò)大到?7.9~0之間。通過(guò)對(duì)無(wú)標(biāo)度區(qū)間的斜率計(jì)算,可得切削速度與關(guān)聯(lián)維數(shù)的關(guān)系。

    切削速度對(duì)關(guān)聯(lián)維數(shù)的影響關(guān)系見(jiàn)圖8(a),不同切削速度對(duì)聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)的影響變化趨勢(shì)大致相同,隨著m的增大,不同切削速度的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)曲線先呈現(xiàn)快速下降的趨勢(shì),并在m=4 之后趨于平緩。不同切削速度的關(guān)聯(lián)維數(shù)變化范圍是有所不同,切削速度由小增大,取值范圍依次為0.0282~0.0656、0.0488~0.0687、0.0389~0.0631、0.0303~0.0583,趨于平穩(wěn)時(shí)的取值分別為0.028、0.0493、0.0393和0.0305,由此可見(jiàn),不同切削速度時(shí)的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)均屬于大于零的正分?jǐn)?shù)。另外,分析相同嵌入維數(shù)時(shí)切削速度對(duì)聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)影響可知,在關(guān)聯(lián)維數(shù)較小(m=1)時(shí),關(guān)聯(lián)維數(shù)隨著進(jìn)給速度的增加基本呈現(xiàn)出線性負(fù)效應(yīng)關(guān)系,提取m=1時(shí)切削速度和關(guān)聯(lián)維數(shù)之間的關(guān)系可得圖8(b),關(guān)聯(lián)維數(shù)與切削速度基本呈現(xiàn)線性負(fù)效應(yīng)關(guān)系,對(duì)其進(jìn)行線性擬合,可得擬合直線D=?0.00051vs+0.007235,擬合相關(guān)系數(shù)為0.66595,顯示出較好的線性擬合關(guān)系。因此,不同切削速度的聲發(fā)射信號(hào)系統(tǒng)均具有混沌動(dòng)力學(xué)特性,利用切削速度與關(guān)聯(lián)維數(shù)的線性負(fù)效應(yīng)關(guān)系,可實(shí)現(xiàn)無(wú)氧銅切削過(guò)程中切削速度特征的評(píng)價(jià)和表征。

    圖7 不同切削速度時(shí)聲發(fā)射信號(hào)的雙對(duì)數(shù)lnC(r)-lnr 圖Fig.7 The double logarithmic lnC(r)-lnr for AE signals with different cutting speeds

    圖8 切削速度對(duì)聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)的影響Fig.8 The influence of cutting speeds on the correlation dimension of AE signals

    3.3.2 切深對(duì)關(guān)聯(lián)維數(shù)的影響分析

    不同切深下聲發(fā)射信號(hào)的雙對(duì)數(shù)曲線如圖9所示。隨著m的增加,雙對(duì)數(shù)曲線的lnC(r)值逐漸增大,不同切深下的雙對(duì)數(shù)曲線變化趨勢(shì)較為相似。當(dāng)m較小時(shí),雙對(duì)數(shù)曲線均呈現(xiàn)出了平穩(wěn)變化特征,隨著m的增大,lnC(r)值呈線性平緩增加趨勢(shì),曲線斜率亦隨之逐漸增大,分形區(qū)域無(wú)標(biāo)度區(qū)間不斷縮小,這種階段性增長(zhǎng)趨勢(shì)最后逐漸收斂于飽和狀態(tài);之后,隨著m的增加,D值基本不再發(fā)生變化。另外,可以發(fā)現(xiàn)的是,隨著切深的增加,lnC(r)取值范圍呈現(xiàn)逐漸增加的趨勢(shì),即從切深為10 μm時(shí)的?6.1~0(圖9(a))逐漸增大到切深為40 μm時(shí)的?7.9~0(圖9(d))。計(jì)算不同切深下的雙對(duì)數(shù)曲線無(wú)標(biāo)度區(qū)間斜率,可得到不同切深時(shí)的關(guān)聯(lián)維數(shù)。

    圖9 不同切深時(shí)聲發(fā)射信號(hào)的雙對(duì)數(shù)lnC(r)-lnr 圖Fig.9 The double logarithmic lnC(r)-lnr for AE signals with different cutting depths

    切深對(duì)聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)的影響特性如圖10(a)所示。隨著m的增大,不同切削深度下的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)均呈現(xiàn)了快速下降的趨勢(shì),并且切深越大,D取值越小;在m為4時(shí),曲線呈現(xiàn)拐點(diǎn)變化,之后,隨著m的增加曲線呈現(xiàn)出趨于平緩且略有所上升的變化趨勢(shì)。另外,不同切深的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)變化范圍是有所不同的,切深由小增大,關(guān)聯(lián)維數(shù)的取值范圍依次為0.0387~0.0826、0.0517~0.0772、0.0439~0.0682、0.0389~0.0631,可見(jiàn),不同切深時(shí)的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)均屬于大于零的正分?jǐn)?shù)。另外,對(duì)比不同嵌入維數(shù)時(shí)聲發(fā)射信號(hào)的關(guān)聯(lián)維數(shù),可以發(fā)現(xiàn),在嵌入維數(shù)取較小值時(shí)(m=1),關(guān)聯(lián)維數(shù)值隨著切深的增大而逐漸減小,關(guān)聯(lián)維數(shù)與切深呈現(xiàn)出明顯的負(fù)效應(yīng)關(guān)系。提取m=1時(shí)的關(guān)聯(lián)維數(shù)與切深的變化情況,即可得圖10(b),可見(jiàn),在m=1時(shí),關(guān)聯(lián)維數(shù)與切深呈現(xiàn)顯著的線性負(fù)效應(yīng)關(guān)系,對(duì)其進(jìn)行線性擬合,可得到擬合直線D=?0.000675ap+0.08965,擬合相關(guān)系數(shù)值為0.98156,顯示出了很好的線性擬合關(guān)系。因此,不同切深下的聲發(fā)射信號(hào)系統(tǒng)均具有混沌動(dòng)力學(xué)特性,可將切深與關(guān)聯(lián)維數(shù)的線性負(fù)效應(yīng)關(guān)系,應(yīng)用于無(wú)氧銅切削過(guò)程中切深特征的評(píng)價(jià)和表征。

    圖10 切深對(duì)聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)的影響Fig.10 The influence of cutting depths on the correlation dimension of AE signals

    3.3.3 進(jìn)給速度對(duì)關(guān)聯(lián)維數(shù)的影響分析

    不同進(jìn)給速度時(shí)的聲發(fā)射信號(hào)雙對(duì)數(shù)曲線如圖11所示。當(dāng)m較小時(shí),雙對(duì)數(shù)曲線斜率較小;隨著m的增加,lnC(r)取值逐漸增大,雙對(duì)數(shù)曲線的斜率亦隨之逐漸增大,分形區(qū)域無(wú)標(biāo)度區(qū)間不斷縮小,雙對(duì)數(shù)曲線呈現(xiàn)階段性增加的趨勢(shì),最后都逐漸收斂于飽和狀態(tài),其中進(jìn)給速度為20 m/min(圖11(c))和25 m/min(圖11(d))時(shí),雙對(duì)數(shù)曲線所呈現(xiàn)的階段性增長(zhǎng)趨勢(shì)更為顯著;在關(guān)維數(shù)達(dá)到飽和之后,隨著m的增加,聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)不再發(fā)生變化。另外,進(jìn)給速度對(duì)lnC(r)取值變化趨勢(shì)的影響并不顯著,計(jì)算不同進(jìn)給速度下的雙對(duì)數(shù)曲線無(wú)標(biāo)度區(qū)間斜率,可得進(jìn)給速度與關(guān)聯(lián)維數(shù)之間的關(guān)系。

    圖11 不同進(jìn)給速度時(shí)聲發(fā)射信號(hào)的雙對(duì)數(shù)lnC(r)-lnr 圖Fig.11 The double logarithmic lnC(r)-lnr for AE signals with different feeding speeds

    進(jìn)給速度對(duì)關(guān)聯(lián)維數(shù)的影響特性如圖12(a)所示,不同進(jìn)給速度對(duì)關(guān)聯(lián)維數(shù)的影響變化趨勢(shì)大致相似,m的增加使D呈現(xiàn)快速下降趨勢(shì),并且進(jìn)給速度越大,D取值越大;m大于4 之后,D的變化逐漸趨于平穩(wěn)。另外,不同進(jìn)給速度時(shí)的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)變化范圍是有所不同的,進(jìn)給速度由小增大,D的取值范圍依次為0.0389~0.0631、

    0.024~0.0747、0.0698~0.0939、0.0861~0.1361,可見(jiàn),不同進(jìn)給速度時(shí)的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)均為大于零的正分?jǐn)?shù)。另外,在嵌入維數(shù)較小(m=1)時(shí),關(guān)聯(lián)維數(shù)隨著進(jìn)給速度的增加而增加,故而提取m=1時(shí)兩者之間的關(guān)系,如圖12(b)所示,可見(jiàn),關(guān)聯(lián)維數(shù)與進(jìn)給速度基本呈線性正效應(yīng)關(guān)系,對(duì)其進(jìn)行線性擬合,擬合直線為D=0.00476vw+0.00858,擬合相關(guān)系數(shù)值為0.88033,具有良好的線性擬合關(guān)系。由此可見(jiàn),不同進(jìn)給速度下的聲發(fā)射信號(hào)同樣是具有混沌動(dòng)力學(xué)特性,根據(jù)進(jìn)給速度與關(guān)聯(lián)維數(shù)的線性正效應(yīng)關(guān)系特征,可應(yīng)用于無(wú)氧銅切削過(guò)程中進(jìn)給速度特征的評(píng)價(jià)和表征。

    圖12 進(jìn)給速度對(duì)聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)的影響Fig.12 The influence of feeding speeds on the correlation dimension of AE signals

    4 討論

    根據(jù)彈塑性切削理論可知[5,13],無(wú)氧銅在切削過(guò)程中的塑性形變功率W可表述為

    其中,σ為材料應(yīng)力,ε為材料應(yīng)變率,V為材料變形體積。根據(jù)聲發(fā)射原理可知,若無(wú)氧銅的σ與ε為常數(shù),無(wú)氧銅塑性形變的聲發(fā)射能W可簡(jiǎn)化為

    經(jīng)推導(dǎo)計(jì)算可知,單顆金剛石以一定的幾何切削參數(shù)切削無(wú)氧銅所產(chǎn)生的聲發(fā)射信號(hào)總功率Wt可描述為

    式(7)中,τ為剪切應(yīng)力,b為切削寬度,α為切削前角,β為剪切角,L為切屑變形長(zhǎng)度。從式(7)可知,無(wú)氧銅切削過(guò)程中的聲發(fā)射信號(hào)能量主要取決于切削過(guò)程中所產(chǎn)生的應(yīng)力、無(wú)氧銅材料的應(yīng)變率以及材料的形變體積。因此,綜合上述金屬切削理論、聲發(fā)射原理、試驗(yàn)結(jié)果與關(guān)聯(lián)維分析,金剛石切削無(wú)氧銅切削的聲發(fā)射信號(hào)機(jī)制可進(jìn)一步分析如下:

    (1)劃痕形貌與聲發(fā)射信號(hào)。從劃痕形貌圖(圖2)可知,金剛石切削無(wú)氧銅經(jīng)歷了“滑擦-耕犁-切削-耕犁-滑擦”的變化過(guò)程,在此過(guò)程中,應(yīng)變率ε和材料應(yīng)力σ將隨著金剛石與無(wú)氧銅接觸形式的變化而變化,其數(shù)值呈現(xiàn)出由小增大再減小的變化歷程,應(yīng)力與應(yīng)變的變化歷程直接映射反饋到切削過(guò)程中的聲發(fā)射時(shí)域信號(hào)振幅(圖3、圖4和圖5)以及能量、RMS 等特征參數(shù)值(表4)等的變化特征上,與聲發(fā)射信號(hào)所呈現(xiàn)出“橄欖型”變化特征是基本相吻合。

    (2)切削參數(shù)對(duì)聲發(fā)射的影響。不同切削參數(shù)工況與聲發(fā)射信號(hào)之間具有一定映射關(guān)系:(a)切削深度。切削深度的增加,金剛石的切削寬度b和無(wú)氧銅切屑變形長(zhǎng)度L亦隨之增加,無(wú)氧銅材料去除體積V得到了明顯的增加,塑性變形程度增大,使其聲發(fā)射時(shí)域信號(hào)振幅增大(圖4),聲發(fā)射信號(hào)能量、振幅和RMS 等特征參數(shù)值(表4)亦隨之增加,切深與聲發(fā)射信號(hào)特征參數(shù)與時(shí)域信號(hào)振幅呈現(xiàn)出顯著的正效應(yīng)關(guān)系。(b)切削速度。結(jié)合式(7)可知,切削速度的增加,無(wú)氧銅材料切屑形變速率增大,瞬時(shí)應(yīng)變率ε增大,聲發(fā)射信號(hào)的總功率隨之增大,能量、RMS等特征參數(shù)值亦隨之有所增加(如表4所示),切削速度與聲發(fā)射信號(hào)總功率基本呈現(xiàn)正效應(yīng)關(guān)系。(c)進(jìn)給速度。進(jìn)給速度的增加,使無(wú)氧銅材料切削瞬時(shí)應(yīng)變率在一定程度上略有所增加,但不及切削深度和切削速度的影響那么顯著,因此,進(jìn)給速度的增加,使聲發(fā)射特征參數(shù)值和時(shí)域信號(hào)振幅在一定范圍內(nèi)波動(dòng)變化,其對(duì)聲發(fā)射信號(hào)的影響并不顯著。

    (3)切削參數(shù)對(duì)關(guān)聯(lián)維數(shù)的影響。無(wú)氧銅切削過(guò)程的聲發(fā)射信號(hào)屬于非線性一維混沌動(dòng)力學(xué)函數(shù),在較小嵌入維數(shù)時(shí)(m=1),不同切削參數(shù)下的聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)分形無(wú)標(biāo)度區(qū)的區(qū)別較大;嵌入維數(shù)的增加,聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)區(qū)域平穩(wěn),分形無(wú)標(biāo)度區(qū)不斷縮小并逐漸趨于飽和狀態(tài),噪聲亦隨之充滿重構(gòu)相空間的每一維,進(jìn)而增加了不必要的數(shù)據(jù)污染。因此,無(wú)氧銅切削聲發(fā)射信號(hào)的嵌入維數(shù)應(yīng)選擇較小值,不同切削參數(shù)對(duì)聲發(fā)射信號(hào)的影響程度不同,關(guān)聯(lián)維數(shù)與不同切削參數(shù)之間的線性相關(guān)程度亦有所差別,其中切削深度和切削速度與關(guān)聯(lián)維數(shù)呈現(xiàn)線性負(fù)效應(yīng)關(guān)系,進(jìn)給速度與關(guān)聯(lián)維數(shù)呈現(xiàn)線性正效應(yīng)關(guān)系。

    5 結(jié)論

    (1)切削速度和進(jìn)給速度對(duì)金剛石切削無(wú)氧銅的聲發(fā)射信號(hào)影響較不顯著,而切削深度的影響則較為顯著,切削深度與聲發(fā)射信號(hào)振幅呈正效應(yīng)關(guān)系。

    (2)金剛石切削無(wú)氧銅的聲發(fā)射信號(hào)雙對(duì)數(shù)曲線均呈現(xiàn)階段性增加趨勢(shì),并逐漸收斂于飽和狀態(tài);嵌入維數(shù)的增加,使聲發(fā)射信號(hào)關(guān)聯(lián)維數(shù)呈現(xiàn)出先快速下降后趨于平穩(wěn)的變化特征。

    (3)金剛石切削無(wú)氧銅的聲發(fā)射信號(hào)具有混沌運(yùn)動(dòng)變化特性,在較小嵌入維數(shù)時(shí),關(guān)聯(lián)維數(shù)與切削深度和切削速度呈現(xiàn)線性負(fù)效應(yīng)關(guān)系,與進(jìn)給速度呈現(xiàn)線性正效應(yīng)關(guān)系。

    猜你喜歡
    相空間切削速度維數(shù)
    β-變換中一致丟番圖逼近問(wèn)題的維數(shù)理論
    束團(tuán)相空間分布重建技術(shù)在西安200 MeV質(zhì)子應(yīng)用裝置的應(yīng)用
    切削速度對(duì)高硬合金高速切削工件溫度和系統(tǒng)勢(shì)能的影響
    切削速度對(duì)鈦合金切屑形貌和剪切帶的影響研究
    一類(lèi)齊次Moran集的上盒維數(shù)
    關(guān)于齊次Moran集的packing維數(shù)結(jié)果
    非對(duì)易空間中的三維諧振子Wigner函數(shù)
    涉及相變問(wèn)題Julia集的Hausdorff維數(shù)
    基于相空間重構(gòu)的電磁繼電器電性能參數(shù)預(yù)測(cè)研究
    基于刀-屑摩擦與切削速度關(guān)聯(lián)模型的切削力數(shù)值分析*
    2018国产大陆天天弄谢| 天美传媒精品一区二区| 大片电影免费在线观看免费| 午夜av观看不卡| 亚洲精品乱码久久久久久按摩| 日韩人妻高清精品专区| 高清黄色对白视频在线免费看| 亚洲第一区二区三区不卡| 成年美女黄网站色视频大全免费 | 国产高清有码在线观看视频| 春色校园在线视频观看| 日本欧美视频一区| 高清黄色对白视频在线免费看| 亚洲av.av天堂| 一区二区三区四区激情视频| 欧美成人午夜免费资源| 你懂的网址亚洲精品在线观看| 视频区图区小说| 国产黄片视频在线免费观看| 成人二区视频| 18+在线观看网站| 交换朋友夫妻互换小说| 国产精品免费大片| 99热全是精品| 两个人免费观看高清视频| 三级国产精品欧美在线观看| 男女边吃奶边做爰视频| 日韩一区二区三区影片| 久久亚洲国产成人精品v| 亚洲av国产av综合av卡| 日本爱情动作片www.在线观看| 在线观看www视频免费| 国产精品人妻久久久久久| 成人毛片a级毛片在线播放| 在线观看三级黄色| 国产亚洲精品久久久com| av在线app专区| 国产精品一区二区在线观看99| 少妇被粗大猛烈的视频| videossex国产| 高清毛片免费看| 大香蕉久久网| 国产精品无大码| 亚洲av电影在线观看一区二区三区| 少妇人妻精品综合一区二区| 国产成人一区二区在线| 国产精品成人在线| 国产免费视频播放在线视频| 精品亚洲成a人片在线观看| 国产精品一区二区在线观看99| 亚洲av成人精品一二三区| 一区二区三区四区激情视频| av在线app专区| 精品国产一区二区三区久久久樱花| 69精品国产乱码久久久| 亚洲,一卡二卡三卡| 亚洲国产成人一精品久久久| 精品国产一区二区三区久久久樱花| 99热网站在线观看| 一本久久精品| 在线观看人妻少妇| 日韩人妻高清精品专区| 黄色怎么调成土黄色| 蜜桃国产av成人99| 国产免费现黄频在线看| 亚洲av.av天堂| 免费看av在线观看网站| 青春草视频在线免费观看| 免费观看av网站的网址| 男女高潮啪啪啪动态图| 在线免费观看不下载黄p国产| 久久毛片免费看一区二区三区| 亚洲激情五月婷婷啪啪| 中国国产av一级| 丰满饥渴人妻一区二区三| 欧美激情 高清一区二区三区| 国产精品国产三级专区第一集| a级毛片免费高清观看在线播放| 狂野欧美白嫩少妇大欣赏| 久久毛片免费看一区二区三区| 男女边摸边吃奶| 男女国产视频网站| 国产精品蜜桃在线观看| 高清午夜精品一区二区三区| 亚洲内射少妇av| 色婷婷久久久亚洲欧美| 亚洲精华国产精华液的使用体验| 91国产中文字幕| 国产综合精华液| 秋霞在线观看毛片| 国产精品一区二区三区四区免费观看| 国产一区亚洲一区在线观看| 51国产日韩欧美| 少妇被粗大的猛进出69影院 | 亚洲色图 男人天堂 中文字幕 | 成人无遮挡网站| 国产日韩欧美亚洲二区| 三级国产精品欧美在线观看| 多毛熟女@视频| 最近2019中文字幕mv第一页| 亚洲精品美女久久av网站| 精品国产乱码久久久久久小说| 中文字幕最新亚洲高清| 香蕉精品网在线| 亚洲中文av在线| 少妇人妻久久综合中文| 国产精品久久久久久精品电影小说| 男人添女人高潮全过程视频| 国产精品成人在线| 18禁在线无遮挡免费观看视频| 新久久久久国产一级毛片| 欧美xxⅹ黑人| 男女边吃奶边做爰视频| 久久久久久久大尺度免费视频| 伦理电影大哥的女人| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色综合www| 欧美97在线视频| 午夜日本视频在线| 美女福利国产在线| 在线观看免费视频网站a站| 最新的欧美精品一区二区| 黑人欧美特级aaaaaa片| 十八禁高潮呻吟视频| 永久网站在线| 免费av中文字幕在线| 国产亚洲一区二区精品| 久久久久久久久久人人人人人人| 欧美国产精品一级二级三级| 女的被弄到高潮叫床怎么办| 满18在线观看网站| 久久99热这里只频精品6学生| 亚洲精品国产av蜜桃| 精品国产一区二区三区久久久樱花| 精品午夜福利在线看| 欧美日韩av久久| 国产伦理片在线播放av一区| 色婷婷av一区二区三区视频| 大香蕉97超碰在线| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 性色avwww在线观看| 成人综合一区亚洲| 久久精品国产自在天天线| 午夜91福利影院| 男人添女人高潮全过程视频| 少妇被粗大猛烈的视频| 黄色视频在线播放观看不卡| 精品酒店卫生间| 97超视频在线观看视频| 亚洲精品第二区| 九草在线视频观看| 国国产精品蜜臀av免费| 乱人伦中国视频| 成人毛片60女人毛片免费| 欧美性感艳星| 99久久人妻综合| 午夜福利在线观看免费完整高清在| 国产精品三级大全| 精品99又大又爽又粗少妇毛片| 国产av码专区亚洲av| √禁漫天堂资源中文www| 在线免费观看不下载黄p国产| www.av在线官网国产| 亚洲国产av新网站| 99热这里只有是精品在线观看| 久久精品人人爽人人爽视色| 男女高潮啪啪啪动态图| av在线app专区| 日韩av不卡免费在线播放| 国产高清有码在线观看视频| 两个人免费观看高清视频| av免费观看日本| 精品久久久精品久久久| 十八禁网站网址无遮挡| 色哟哟·www| 蜜桃国产av成人99| 欧美精品一区二区大全| 亚洲国产最新在线播放| 99国产精品免费福利视频| 亚洲国产精品专区欧美| 啦啦啦啦在线视频资源| 欧美精品一区二区大全| 99九九在线精品视频| 久久久精品免费免费高清| 亚洲精品日韩av片在线观看| 男人爽女人下面视频在线观看| 国产毛片在线视频| 午夜激情久久久久久久| 中文字幕亚洲精品专区| 午夜日本视频在线| 亚洲av免费高清在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美激情 高清一区二区三区| 精品国产一区二区久久| 久久精品夜色国产| 国产av码专区亚洲av| 国产精品一区二区在线不卡| av在线老鸭窝| 视频在线观看一区二区三区| 免费av中文字幕在线| 午夜福利影视在线免费观看| 亚洲四区av| 嘟嘟电影网在线观看| 五月开心婷婷网| 最新中文字幕久久久久| 热re99久久国产66热| 国产亚洲精品第一综合不卡 | 亚洲怡红院男人天堂| 午夜免费鲁丝| 另类精品久久| 午夜免费男女啪啪视频观看| 成年美女黄网站色视频大全免费 | 飞空精品影院首页| 中国美白少妇内射xxxbb| 亚洲av日韩在线播放| 亚洲一级一片aⅴ在线观看| 18禁在线无遮挡免费观看视频| 国产精品久久久久久精品古装| 伦精品一区二区三区| 国产精品成人在线| 高清av免费在线| 亚洲四区av| 在线观看一区二区三区激情| 肉色欧美久久久久久久蜜桃| 国产极品粉嫩免费观看在线 | 成人二区视频| 久久精品人人爽人人爽视色| 亚洲av.av天堂| 青青草视频在线视频观看| 少妇人妻精品综合一区二区| 国产毛片在线视频| 欧美日韩在线观看h| 最后的刺客免费高清国语| 精品亚洲乱码少妇综合久久| 亚洲精品一二三| 免费高清在线观看视频在线观看| 国产高清国产精品国产三级| 久久久精品区二区三区| 亚洲久久久国产精品| 哪个播放器可以免费观看大片| 久久午夜综合久久蜜桃| 久久久久久久大尺度免费视频| 欧美人与善性xxx| 国产成人freesex在线| 波野结衣二区三区在线| 边亲边吃奶的免费视频| 亚洲av男天堂| 国产日韩一区二区三区精品不卡 | 18+在线观看网站| 免费观看av网站的网址| 久久久久人妻精品一区果冻| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 久久99精品国语久久久| 搡老乐熟女国产| 国产69精品久久久久777片| 亚洲三级黄色毛片| 美女大奶头黄色视频| 高清av免费在线| 欧美精品亚洲一区二区| 高清视频免费观看一区二区| 亚洲精品,欧美精品| 91久久精品国产一区二区成人| 黑人巨大精品欧美一区二区蜜桃 | 大陆偷拍与自拍| 国产成人91sexporn| 国产片特级美女逼逼视频| 国产精品麻豆人妻色哟哟久久| 这个男人来自地球电影免费观看 | a 毛片基地| 亚洲丝袜综合中文字幕| 午夜福利网站1000一区二区三区| 国国产精品蜜臀av免费| 亚洲欧美日韩另类电影网站| 18+在线观看网站| 国产一区二区三区av在线| 丝袜美足系列| 大又大粗又爽又黄少妇毛片口| 久久99精品国语久久久| 国产综合精华液| 在线观看www视频免费| 人妻制服诱惑在线中文字幕| www.色视频.com| 一级黄片播放器| 精品人妻在线不人妻| 插逼视频在线观看| 久久久国产一区二区| 99久久精品国产国产毛片| 免费黄网站久久成人精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 少妇的逼好多水| 在线 av 中文字幕| 亚洲国产成人一精品久久久| 精品一品国产午夜福利视频| 亚洲av中文av极速乱| 国产免费视频播放在线视频| 精品久久久噜噜| 91国产中文字幕| 一区二区三区免费毛片| 久久狼人影院| 亚洲美女黄色视频免费看| 精品卡一卡二卡四卡免费| 夜夜爽夜夜爽视频| 中文欧美无线码| 亚洲欧美清纯卡通| 久久久欧美国产精品| 国产一区二区在线观看日韩| 国产精品嫩草影院av在线观看| av专区在线播放| 亚洲精品视频女| 欧美日韩视频高清一区二区三区二| av女优亚洲男人天堂| 国产不卡av网站在线观看| 一区二区三区精品91| 丰满少妇做爰视频| 精品人妻一区二区三区麻豆| 22中文网久久字幕| a级毛片免费高清观看在线播放| 婷婷成人精品国产| 欧美最新免费一区二区三区| 一级毛片 在线播放| 美女国产高潮福利片在线看| 在线播放无遮挡| 国产伦精品一区二区三区视频9| 曰老女人黄片| 80岁老熟妇乱子伦牲交| 一区二区av电影网| 春色校园在线视频观看| 久久久精品免费免费高清| freevideosex欧美| 18禁在线无遮挡免费观看视频| 一区二区日韩欧美中文字幕 | 亚洲国产日韩一区二区| 只有这里有精品99| 国产亚洲精品久久久com| 午夜免费男女啪啪视频观看| 看非洲黑人一级黄片| 国产精品熟女久久久久浪| 午夜免费男女啪啪视频观看| freevideosex欧美| 日韩成人av中文字幕在线观看| 国产爽快片一区二区三区| 老司机影院毛片| 高清av免费在线| 在线观看www视频免费| 日韩熟女老妇一区二区性免费视频| 亚洲国产欧美在线一区| 免费黄色在线免费观看| 精品酒店卫生间| 国产亚洲午夜精品一区二区久久| 中文字幕免费在线视频6| 成人毛片60女人毛片免费| 久久精品久久精品一区二区三区| 国产色爽女视频免费观看| 欧美亚洲 丝袜 人妻 在线| 日韩欧美精品免费久久| 人妻一区二区av| 亚洲一级一片aⅴ在线观看| 久久精品夜色国产| 日韩,欧美,国产一区二区三区| 精品久久久久久电影网| 亚洲美女黄色视频免费看| 欧美变态另类bdsm刘玥| 亚洲伊人久久精品综合| 性色av一级| 男的添女的下面高潮视频| 少妇猛男粗大的猛烈进出视频| 热99久久久久精品小说推荐| 七月丁香在线播放| 999精品在线视频| 国产成人精品无人区| 久久精品久久精品一区二区三区| 亚洲美女搞黄在线观看| 亚洲国产精品一区二区三区在线| 啦啦啦视频在线资源免费观看| 亚洲精品一二三| 久久久久久伊人网av| 五月天丁香电影| 国产毛片在线视频| 国产一区二区在线观看av| 亚洲人与动物交配视频| 国产在线免费精品| 97在线视频观看| 免费少妇av软件| 日本vs欧美在线观看视频| 久久久精品区二区三区| 两个人的视频大全免费| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 97超视频在线观看视频| h视频一区二区三区| 有码 亚洲区| 午夜视频国产福利| 狂野欧美激情性xxxx在线观看| 18+在线观看网站| 一级二级三级毛片免费看| 欧美+日韩+精品| 亚洲精品久久午夜乱码| .国产精品久久| 国语对白做爰xxxⅹ性视频网站| 全区人妻精品视频| 成年人午夜在线观看视频| 午夜福利网站1000一区二区三区| 亚洲欧美日韩另类电影网站| 亚洲熟女精品中文字幕| 乱人伦中国视频| 我的女老师完整版在线观看| 日韩精品免费视频一区二区三区 | 成人国产av品久久久| 中国三级夫妇交换| 男女无遮挡免费网站观看| 欧美国产精品一级二级三级| 伦精品一区二区三区| 国精品久久久久久国模美| 国产精品一区二区在线观看99| 午夜老司机福利剧场| 免费观看av网站的网址| 狂野欧美白嫩少妇大欣赏| 国产精品偷伦视频观看了| 久热久热在线精品观看| 日本vs欧美在线观看视频| 久久久久久久亚洲中文字幕| 男女高潮啪啪啪动态图| 精品人妻熟女毛片av久久网站| 一二三四中文在线观看免费高清| 高清午夜精品一区二区三区| 九色亚洲精品在线播放| 如日韩欧美国产精品一区二区三区 | 精品亚洲乱码少妇综合久久| 日韩三级伦理在线观看| 久久毛片免费看一区二区三区| 欧美性感艳星| 看十八女毛片水多多多| 国产精品熟女久久久久浪| 男女啪啪激烈高潮av片| 亚洲国产欧美日韩在线播放| 亚洲国产av影院在线观看| 国产一区有黄有色的免费视频| 高清av免费在线| 国产精品久久久久久精品古装| 国产黄频视频在线观看| 九九爱精品视频在线观看| 日韩大片免费观看网站| 亚洲精品久久成人aⅴ小说 | 精品一品国产午夜福利视频| 欧美日韩视频高清一区二区三区二| av有码第一页| 色婷婷久久久亚洲欧美| 日本与韩国留学比较| 久久久久久久大尺度免费视频| 美女国产视频在线观看| 青春草视频在线免费观看| 中文字幕最新亚洲高清| 日韩熟女老妇一区二区性免费视频| 日韩免费高清中文字幕av| 亚洲精品aⅴ在线观看| 国产乱来视频区| 国产一区二区在线观看日韩| 国产精品麻豆人妻色哟哟久久| 免费av中文字幕在线| 精品少妇黑人巨大在线播放| 一级毛片我不卡| 午夜福利网站1000一区二区三区| 另类精品久久| 日本vs欧美在线观看视频| 男女无遮挡免费网站观看| 在线观看国产h片| 女性生殖器流出的白浆| 中文乱码字字幕精品一区二区三区| 免费高清在线观看视频在线观看| 欧美成人午夜免费资源| 日韩av免费高清视频| 亚洲国产精品999| 大片免费播放器 马上看| 午夜福利在线观看免费完整高清在| 亚洲欧美日韩另类电影网站| 国产极品粉嫩免费观看在线 | 精品一区二区免费观看| 精品一品国产午夜福利视频| 校园人妻丝袜中文字幕| 中文字幕人妻丝袜制服| a级毛色黄片| 亚洲怡红院男人天堂| 成年女人在线观看亚洲视频| 人人妻人人添人人爽欧美一区卜| 亚洲经典国产精华液单| 熟女人妻精品中文字幕| 国产熟女午夜一区二区三区 | 一级二级三级毛片免费看| 成人亚洲精品一区在线观看| 久久久久精品性色| 欧美人与性动交α欧美精品济南到 | 国产精品欧美亚洲77777| 国产探花极品一区二区| 亚洲四区av| 日本与韩国留学比较| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 久久女婷五月综合色啪小说| a 毛片基地| 3wmmmm亚洲av在线观看| 成人综合一区亚洲| 免费观看性生交大片5| 国产精品无大码| 久久久久久久久久久丰满| 久久久国产一区二区| 亚洲在久久综合| 中文字幕精品免费在线观看视频 | 观看美女的网站| 你懂的网址亚洲精品在线观看| 国产日韩欧美视频二区| 少妇人妻久久综合中文| 18+在线观看网站| 大片免费播放器 马上看| 天天躁夜夜躁狠狠久久av| 亚洲精品中文字幕在线视频| 各种免费的搞黄视频| 色5月婷婷丁香| 王馨瑶露胸无遮挡在线观看| 伊人亚洲综合成人网| 纯流量卡能插随身wifi吗| 一级片'在线观看视频| 国产一区二区三区综合在线观看 | 中文字幕久久专区| 国产av国产精品国产| 日本黄大片高清| 国产69精品久久久久777片| 日韩大片免费观看网站| 亚洲性久久影院| 夜夜骑夜夜射夜夜干| 91久久精品电影网| 久久久久久人妻| 一边亲一边摸免费视频| 九草在线视频观看| 国产精品一区www在线观看| 99热国产这里只有精品6| 男女无遮挡免费网站观看| 亚洲精品av麻豆狂野| 国产免费一级a男人的天堂| 国产精品蜜桃在线观看| 色哟哟·www| 色94色欧美一区二区| 少妇人妻精品综合一区二区| 久久久久久伊人网av| 欧美激情 高清一区二区三区| 日韩免费高清中文字幕av| 黄片无遮挡物在线观看| 国产免费又黄又爽又色| 啦啦啦视频在线资源免费观看| 欧美精品人与动牲交sv欧美| 国产深夜福利视频在线观看| 91久久精品国产一区二区三区| 99精国产麻豆久久婷婷| 国产69精品久久久久777片| 亚洲精品色激情综合| 国产成人freesex在线| 亚洲国产精品成人久久小说| 人人妻人人澡人人看| 久久毛片免费看一区二区三区| 午夜91福利影院| 成人手机av| 又黄又爽又刺激的免费视频.| 精品少妇久久久久久888优播| 我的老师免费观看完整版| 中文字幕制服av| 男女边吃奶边做爰视频| 黑丝袜美女国产一区| 欧美丝袜亚洲另类| 国产成人精品婷婷| 亚洲av二区三区四区| a级毛色黄片| 精品久久蜜臀av无| 十八禁网站网址无遮挡| 午夜激情av网站| 国产成人91sexporn| 伦精品一区二区三区| 精品一区二区三区视频在线| 精品国产一区二区久久| 尾随美女入室| 国产片特级美女逼逼视频| 草草在线视频免费看| 亚洲无线观看免费| 美女国产高潮福利片在线看| 国产伦精品一区二区三区视频9| 日韩成人伦理影院| 日韩亚洲欧美综合| 又大又黄又爽视频免费| 另类亚洲欧美激情| 我要看黄色一级片免费的| 少妇被粗大猛烈的视频| 18+在线观看网站| 国内精品宾馆在线| 亚洲欧洲国产日韩| 色5月婷婷丁香| 国产精品久久久久久精品电影小说| 亚洲怡红院男人天堂| 伦精品一区二区三区| a 毛片基地| 我的女老师完整版在线观看| 久久鲁丝午夜福利片| 亚洲精品第二区| 国产精品人妻久久久影院| 日韩熟女老妇一区二区性免费视频| 日韩一区二区三区影片| 啦啦啦在线观看免费高清www| 国产精品一区二区在线观看99| 久久久精品区二区三区| 九色成人免费人妻av| 99国产综合亚洲精品| 国产精品蜜桃在线观看| 国产精品一区二区在线不卡| 久久99一区二区三区| 乱人伦中国视频| 校园人妻丝袜中文字幕|