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

    超聲空化對軟性磨粒流切削效率和質(zhì)量的影響

    2017-07-18 11:48:14計時鳴譚大鵬鄭晨亮
    農(nóng)業(yè)工程學報 2017年12期
    關(guān)鍵詞:硅片磨粒空化

    計時鳴,陳 凱,2※,譚大鵬,鄭晨亮

    ?

    超聲空化對軟性磨粒流切削效率和質(zhì)量的影響

    計時鳴1,陳 凱1,2※,譚大鵬1,鄭晨亮1

    (1. 浙江工業(yè)大學機械工程學院,杭州 310014;2. 麗水學院浙江省特色文創(chuàng)產(chǎn)品數(shù)字化設(shè)計與智能制造重點實驗室,麗水 323000)

    為提高軟性磨粒流切削的效能,研究了利用超聲空化作用加快材料去除的方法。在分析磨粒拋光試驗和計算流體動力學仿真結(jié)果的基礎(chǔ)上,提出并驗證磨粒碰撞壁面的動能是影響切削效率主要因素的假設(shè);根據(jù)對超聲空化氣泡潰滅,及其在外圍水域中沖擊波形成與傳播的研究結(jié)果,分析氣泡潰滅在推動磨粒改變其動能時所起的增強湍流效果的作用。在此基礎(chǔ)上搭建了超聲輔助磨粒流試驗裝置,以表面覆有氧化層的硅片為試驗對象,根據(jù)其表面氧化層的去除情況,分析超聲空化對硅片表面材料去除能力的影響,并在磨粒流加工無氧化層硅片表面的試驗中,觀察到超聲空化對磨粒流切削的顯著影響。驗證了超聲空化能使磨粒流中的磨粒的動能發(fā)生改變,從而起到提高切削效率和提高表面質(zhì)量的作用,以期為超聲波作用下磨粒流加工工藝的改進提供參考。

    超聲;空化;磨削;軟性磨粒流;三相流;磨粒運動

    0 引 言

    磨粒流加工是利用攜帶磨粒的流體對工件表面進行精整加工的一種加工方法,產(chǎn)生于20世紀60年代,它以流體中的磨粒作為切削工具,適合于對復(fù)雜結(jié)構(gòu)的零件進行加工。由于流體具有較好的仿形特性,因此在一些特定的加工場合有著不可比擬的優(yōu)勢[1-5]。

    用磨粒流對表面進行拋光,一般都要很長的時間,對鋼材表面進行加工時,表面粗糙度從Ra 0.9m降至0.3m要用35 h[5]。如何在不影響加工質(zhì)量甚至有所改善的前提下提高加工效率,是研究者非常感興趣且希望解決的問題。在先前的研究中,曾經(jīng)加入低強度的超聲波并觀察到了拋光效率和質(zhì)量改善的現(xiàn)象[6-8],當時認為是加載超聲波的流道出現(xiàn)了大量速度較小、壓力較大的空化區(qū)域,不斷發(fā)展撞擊壁面導(dǎo)致的結(jié)果。本文進一步分析了空化的物理過程,提出氣泡潰滅的沖擊波影響磨粒運動的觀點,對拋光效能改善的機理給出一種新的解釋。

    超聲波在液體中傳播時會發(fā)生空化現(xiàn)象,理論計算和實際觀察都表明,潰滅的氣泡會在其內(nèi)部產(chǎn)生極高的溫度和壓力,從而引發(fā)一系列復(fù)雜的物理化學現(xiàn)象[9-13]。研究表明,氣泡壓縮到最小時,氣泡壁壓強瞬間達到Pa以上,高壓以聲速向泡外傳播,沖擊波幅值大約以傳播距離倒數(shù)的速率衰減[14-16]。前人研究多關(guān)注于氣泡潰滅的微射流、沖擊波對壁面的直接作用,水中的超聲空化作用主要用于清洗,為防止氣蝕的發(fā)生,超聲波的強度不大,常見的喇叭型超聲換能器單個功率一般都在100 W左右。純水中的空化效果主要是由潰滅的氣泡實現(xiàn)的,如果在發(fā)生空化的水中有眾多高速運動的磨粒,磨粒對壁面的切削作用會大大超過氣泡對壁面的直接作用。合理加強超聲強度,利用空化沖擊波影響磨粒的運動,可以增強切削的效能。通過計算表明,沖擊波對磨粒所做的功的大小主要取決于潰滅點和磨粒的距離,以及受超聲波聲壓和水溫的影響很大的潰滅壓力。根據(jù)對磨粒流切削機理的分析,被加速的部分磨粒和隨機變化的運動方向,相當于起到了增加湍動能的效果,可以起到顯著增強切削效率甚至提高表面質(zhì)量的作用。

    1 磨粒流的切削機理

    從原子尺度上分析,軟性磨粒流拋光過程是被加工表面凹凸不平的峰谷結(jié)構(gòu)被大量的單個液體分子和少量的磨粒分子團撞擊的過程[17-21]。撞擊后工件表面的部分分子(或原子)獲得動能,使其產(chǎn)生晶格滑移或者脫離基體,處于“峰”部的分子團,在遭受碰撞后更容易滑移或脫離,首先是因為其受撞擊的概率大,其次是峰頂?shù)姆肿釉谧矒糁袝@得更多的能量而脫離原來的位置。經(jīng)過長時間大量的沖擊后,峰谷結(jié)構(gòu)趨向平緩導(dǎo)致表面粗糙度降低。一種被稱為納米顆粒膠體射流拋光的技術(shù),和磨粒流的切削機理非常類似[22-25]。張成光等[26]研究了磨粒水射流的材料去除機理,指出磨粒水射流加工時,高速高壓的水對磨粒起到的僅僅是載體的作用,實際的加工效果由磨粒實現(xiàn),兩種加工方法的主要區(qū)別是一個是對點加工,一個是對面加工。圖1是參考文獻[27]得到的試驗結(jié)果。

    圖1 納米顆粒膠體射流拋光試驗和仿真結(jié)果

    圖1a是用流道直徑為0.5 mm的普通錐形噴嘴對工件表面進行射流的垂直噴射2 h后,使用激光共聚焦顯微鏡觀察得到的圓環(huán)狀圖像,圓環(huán)最深處半徑在0.9 mm左右。值得注意的是,加工區(qū)域中心處的材料去除較少,去除最深處位于距拋光中心位置距離相等的圓環(huán)地帶。圖1b是使用輪廓儀對加工區(qū)域的檢測結(jié)果,非常明顯其截面形狀呈“W”形。圖1c是徑向截面上分析得到的動靜壓壓力曲線,截面直徑上的動壓壓力大小呈倒“W”型分布,射流中心動壓壓力趨于減小,動壓壓力高的區(qū)域集中在距中心一定距離的環(huán)形區(qū)域內(nèi)。試驗表明材料的去除效率和流體的動壓有關(guān),而跟靜壓無關(guān)。流體力學中動壓的定義為,磨粒流的密度在加工前后是不變的,這說明材料的去除效率和流體速度的平方成近似的線性關(guān)系。

    分析流體內(nèi)水、磨粒分子的情況運動,水分子雖然數(shù)量眾多,但由于尺寸相對于磨粒很小,在同樣速度的情況下,單個的水分子很難撞開工件表面的分子,切削作用主要靠磨粒來完成??梢哉J為水或其他流體主要起載體的作用,幫助磨粒獲得速度,以一定的角度撞擊加工表面。這一點可以從純水和攜帶磨粒的混合液的拋光效果對比就可以得到驗證。

    磨粒水射流加工試驗表明,在其它條件不變的情況下(加工時間,磨料的材料、體積分數(shù)、大小,被加工材料,環(huán)境溫度等),材料的去除量用某個位置點上的切削深度來表示的話,就有:

    該結(jié)論在軟性磨粒流拋光試驗也可以得到驗證,參見如圖2所示的磨粒流拋光試驗裝置,用泵向三根軟管注入磨粒和水的混合液,進入一個設(shè)計了特殊流道用來產(chǎn)生所需要流場的加工工具,如圖2a所示。在圓形的加工工具下方固定一個直徑一樣大小的鋁制圓盤,利用磨粒流對其上表面進行拋光。圖2b給出了磨粒流流動區(qū)域簡化模型,由于切削在工件壁面上進行,需要重點考察被加工壁面上方流場情況,流動區(qū)域最下面的大面積拋光區(qū)是分析的重點。

    圖2 磨粒流拋光試驗裝置及磨粒流流動區(qū)域簡化模型

    Fig.2 Abrasive flow polishing experiment device and simulation model of abrasive flow area

    拋光試驗中,液相為水,固相為碳化硅(SiC)顆粒,顆粒直徑為13m,密度為。選用鋁作為加工工件,工件初始表面粗糙度為。對鋁加工工件拋光試驗,并對該試驗條件下進行拋光模型仿真。

    利用Gambit 2.4.6軟件,對圖2b所示的試驗?zāi)チA髁鲃訁^(qū)域進行網(wǎng)格劃分后,再使用FLUENT軟件,建立基于Realizable 湍流模型和Mixture模型的液-固兩相磨粒流拋光加工模型,進行數(shù)值模擬。磨粒流3個入口直徑均為20 mm,磨粒流入口設(shè)置為速度入口,磨粒流入口速度為=60m/s,磨粒流出口設(shè)置為自由出口,體積分數(shù)為0.1(磨粒體積占混合液總體積的10%),其余邊界設(shè)置為固體邊界。

    圖3a為仿真得到的工件表面磨粒流速度矢量圖,在圖中明顯看到中間有一個流速較高的圓環(huán)區(qū)域。圖3b為仿真得到的工件表面磨粒流湍動能分布云圖,在圖中同樣可以看到中間有一個湍動能較高的圓環(huán)區(qū)域,說明該區(qū)域磨粒撞擊方向更混亂和隨機。圖3c為兩相流沖刷0.5 h后得到的工件表面效果圖,和仿真得到的湍動能云圖及速度矢量圖比較,可以明顯地看出:湍動能強,磨粒速度大的區(qū)域可以獲得較高的拋光效率。從磨粒速度矢量圖上,可以明顯看到被加工圓盤上的三個不同的速度區(qū)域,中間的圓環(huán)區(qū)域磨粒速度最高,用不同尺寸的磨粒進行仿真,得到的結(jié)果相近。

    FLUENT中使用湍流動能(TKE,turbulence kinetic energy)來描述湍流中單位質(zhì)量的平均動能。

    式中為湍動能,m2/s2;,,為三個直角坐標軸方向上的脈動速度,m/s。

    近壁面處的湍流動能可以反映出磨粒流在近壁面處攜帶的動能,也在一定程度上反應(yīng)出湍流的強度,因此也可以用湍流動能來估計磨粒撞向壁面方向的混亂程度,湍流強度足夠大時,可以認為各個方向的撞擊是均勻的。

    在圖3c中,也可以發(fā)現(xiàn)速度大、湍動能強的圓環(huán)形區(qū)域材料被去除的快,拋光效率高,比其他區(qū)域更早達到了較高的表面粗糙度。中間的圓形區(qū)域,由于磨粒速度最低,拋光效果最差。圓盤表面加工前粗糙度值為0.8m,用1 000目直徑13m的碳化硅(SiC)顆粒拋光鋁,2 h后材料去除最快區(qū)域的Ra從0.8m 降到0.6m。而用500目直徑25的碳化硅(SiC)顆粒拋光,1 h后就到0.6m左右??紤]到工件表面達到的粗糙度和材料去除量有關(guān),那么到達某粗糙度所用的時間也可以反映材料去除效率,可以得出:材料的去除效率和磨粒的大小有關(guān)。

    材料的去除效率和流體的中的磨粒的大?。ㄙ|(zhì)量)有關(guān),圖1的試驗又說明去除效率(相同時間內(nèi)的切削深度)和速度的平方成某種線性關(guān)系,因此可以合理地假設(shè),材料的去除效率和磨粒所攜帶的動能成近似的線性關(guān)系。即:磨粒所攜帶的動能和材料的去除率有關(guān)。

    不同大小的磨料最終所得到的表面粗糙度是不一樣的,但這不影響材料的去除效率和磨粒所攜帶的動能成近似線性關(guān)系的分析。

    2 空化沖擊波對磨粒動能的影響

    到目前為止對空化的機理研究尚未完全清楚。一般可以用改進的Rayleigh方程進行理論計算[28-29]。

    液體在超聲場中不再被當作不可壓縮流體,其質(zhì)點間的疏密隨聲波而變化,液體中氣泡的大小也隨超聲的疏密相而變化,如圖4所示,當其直徑小到為零可以認為其發(fā)生潰滅(閉合)。

    以氣泡壁面為研究對象,可以分析出氣泡的受力情況,泡內(nèi)有氣體壓力和蒸汽壓,泡外有液體靜壓力,表面張力以及超聲壓力。對氣泡大小變化的過程作適當?shù)暮喕诓豢紤]蒸氣壓力和輻射阻尼的情況下,經(jīng)后人改進的Rayleigh-plesset運動方程為[30]:

    其中為液體的表面張力系數(shù),N/m;為氣泡初始半徑,m;為超聲的聲壓,Pa;為液體靜壓,Pa;為液體的黏度,Pa· s;為氣體的多變指數(shù);為氣泡半徑,m。

    利用龍格-庫塔(Rung-Kutta)法[8]可解出該常微分方程,得到氣泡直徑在超聲場內(nèi)隨時間變化的規(guī)律,當氣泡半徑為零時,即發(fā)生氣泡潰滅。

    利用MATLAB計算單個氣泡在20 kHz聲場里半徑大小隨時間的變化可得如圖5的結(jié)果。

    圖5 水中不同大小的氣泡在超聲波作用下直徑變化的計算結(jié)果

    氣泡半徑在0.5m,500m時,穩(wěn)態(tài)空化(不潰滅);氣泡半徑在50m時,瞬態(tài)空化(潰滅)。計算結(jié)果說明單個氣泡在50m時,一接觸到超聲波在一個超聲周期內(nèi)(50s)就會潰滅。進一步可以分析出:成千上萬差不多大小的氣泡,只要在強度足夠的聲場里同樣會瞬間潰滅。

    可以用一個簡單而又令人信服的試驗證實這一分析:在盛水的容器內(nèi)加入大量由微納米氣泡泵所產(chǎn)生半徑在50m左右的氣泡,使之變成牛奶似的不透明的乳濁液,發(fā)射超聲波容器內(nèi)的水瞬間變透明。

    根據(jù)以上分析的結(jié)果,可以得出這樣一個結(jié)論:在一定頻率、一定強度的超聲波作用下,水中只有半徑小于最小閾值,或者半徑大于最大閾值的氣泡存在并作振蕩。小氣泡(氣核)在振蕩過程中因得到氣化補充不斷增大直至潰滅,大氣泡在震蕩過程上升直至脫離液體。

    通過計算還可以知道,超聲的頻率越低,越容易發(fā)生空化,對于一個變幅桿,換能器的機械尺寸已經(jīng)確定的超聲發(fā)生系統(tǒng),因為機電諧振頻率和機械共振一樣,是不可以隨意改變的,本文使用20 kHz這一頻率。

    假設(shè)瞬態(tài)空化泡的收縮(直到崩潰)過程是絕熱過程,則空化泡崩潰時泡內(nèi)的最高溫度()與最大壓力()可用下式表示[7]:

    (4)

    計算磨粒速度在氣泡潰滅時的變化,需要了解空化泡液體外圍的壓強??栈菀后w外圍壓強的分布,許多學者作了研究[14-16],文獻[16]給出了相關(guān)結(jié)果。

    假設(shè)氣泡保持一個球形變化,數(shù)值計算解氣泡外圍水中流體連續(xù)性方程和Navier-Stokes方程,在球?qū)ΨQ的條件下,形式為

    (6)

    采用與水蒸氣有關(guān)的改進的較精確的氣泡壁的運動RP方程為[16]

    計算發(fā)現(xiàn)[16],當氣泡接近它的最小半徑處,氣泡外圍形成很薄的高壓區(qū),這個高壓區(qū)域隨后向外傳播形成沖擊波。其研究結(jié)果顯示,氣泡外能觀測到向外的沖擊波,氣泡壓縮到最小時,氣泡壁壓強瞬間達到上千個以上大氣壓,高壓與氣泡內(nèi)有無沖擊波無關(guān)。尖銳的沖擊波脈沖在氣泡從最小半徑處反彈后形成,沖擊波形成后向外傳播,同時很快地衰減,沖擊波幅值大約以傳播距離倒數(shù)的速率衰減,沖擊波以聲速即1 480 m/s的速度向外傳播。

    利用這個模型,可以分析氣泡附近的磨粒因沖擊波而獲得的動能。為簡化數(shù)學上的處理,假設(shè)磨粒靜止,且受沖擊波作用的面積不變,磨粒距氣泡潰滅點的最近距離和最遠距離分別為和,如圖6所示。假設(shè)氣泡潰滅點的沖擊波壓強為,由于沖擊波幅值以傳播距離倒數(shù)的速率衰減即有,為離潰滅點距離等于處的壓強,為比例常數(shù),由于磨粒受沖擊波作用的面積也為常數(shù),那么沖擊波對磨粒做的功為

    磨粒的尺寸也會影響其獲得的動能,尺寸越大獲得的動能也越大,但其影響程度不如距離對獲得動能的影響,特別是在磨粒和潰滅點靠得很近的時候。

    如果磨粒原來靜止不動,并且相對于潰滅點的距離確定,就可以計算出磨粒所獲得的動能。但實際情況是潰滅點的位置無法預(yù)知,磨粒在氣泡潰滅前的運動方向無法預(yù)知(運動方向決定了磨??赡鼙患铀僖部赡鼙粶p速),磨粒與潰滅點的距離無法預(yù)知,很難對某一顆磨粒的運動進行計算。參照FLUENT對湍流的研究方法,可以從統(tǒng)計的角度來分析問題,先進行試驗,在試驗觀察的基礎(chǔ)上給出半經(jīng)驗的解釋和預(yù)測。

    3 超聲空化對磨粒流切削影響的試驗研究和分析

    根據(jù)空化機理,超聲強度顯然會影響空化強度??偟膩碚f強度越大效果越明顯,小功率超聲波只起到清洗的效果,大功率超聲波的空化能直接破壞被加工件表面[31-32]。拋光時要考慮被加工材料對超聲的耐受能力,不能發(fā)生氣蝕破壞表面。試驗發(fā)現(xiàn)不同的材料差別較大,鋁很脆弱,不銹鋼,玻璃,硅片,甚至塑料耐受能力很強。需要通過試驗找到合理的強度,留有安全裕度,讓空化時氣泡潰滅所造成的“沖擊波”和“微射流”能盡可能地推動磨粒加速運動又不會傷害被加工表面。

    3.1 合理超聲強度的確定

    通過大功率超聲波在液體中產(chǎn)生的空化作用,可以高速度模擬材料的氣蝕狀況。參照美國材料試驗協(xié)會編制的氣蝕試驗標準(ASTM)G32-06,使用一個頻率在20 kHz左右,功率可調(diào)的超聲波發(fā)生器,驗證氣蝕情況。

    超聲在水中的強度用聲強來描述,實際操作中超聲變幅桿發(fā)射端面和被加工表面只有1 mm左右的距離,不便于測量聲壓??梢允褂贸暟l(fā)射器功率,結(jié)合超聲發(fā)射面積來計算超聲強度。

    圖7為試驗裝置示意圖,在超聲換能器的兩電極上加一定頻率的交流電壓,壓電陶瓷因伸縮變形產(chǎn)生諧振,變幅桿將小振幅的振動放大,經(jīng)由下端面?zhèn)鞒觥:统R?guī)氣蝕試驗裝置不同的是,在超聲變幅桿的下方開孔,利用該孔由液壓泵壓入流動的水或者混合液,這樣可以比較不同的情況下的氣蝕程度。

    圖7 超聲空化磨粒流試驗裝置示意圖

    不論使用哪種材料試驗,被加工材料表面上都有圖8所示的一大一小2個圓環(huán)形的區(qū)域。小圓環(huán)稱之為兩相流切削顯著區(qū),通過FLUENT仿真可以得到該區(qū)域內(nèi)磨粒速度和湍動能都是最大的如圖9a所示,在水、磨粒兩相流動的狀態(tài)下,該區(qū)域的加工效率是最高的。由于壓電陶瓷是做成圓環(huán)形的,其振動產(chǎn)生的聲波經(jīng)變幅桿聚集后傳播到硅片上仍為圓環(huán)形區(qū)域,該大圓環(huán)區(qū)是超聲空化區(qū),如果發(fā)射超聲波,該區(qū)域會發(fā)生持續(xù)的氣泡潰滅現(xiàn)象,形成水、氣泡、磨粒的三相流動。

    試驗使用的硅片是經(jīng)過化學機械拋光的成品晶元(wafer),其表面經(jīng)氧化后有一層幾微米厚的SiO2,由于有少量雜質(zhì),在顯微鏡下呈黃色,氧化層之下是灰黑色的硅,如果表層脫落,可以很容易觀察到。

    圖8 被加工件表面的不同區(qū)域示意圖

    通過泵入純水進行試驗,找到合理的超聲波強度,在此強度下把純水換成不同尺寸磨粒的混合液,保證了既能觀測到空化效果,又不會破壞被加工表面。

    3.2 空化作用對磨粒切削能力影響的驗證試驗

    為便于觀察,選用表面附有一層SiO2的硅片(顯微鏡下SiO2黃色,硅黑色)進行試驗,可以通過觀察黃色SiO2的脫落情況,分析磨粒的切削效果。共進行了4次試驗:

    1)沒有觀察到超聲空化作用的顯著影響:使用100目,磨粒直徑150m左右的氧化鋁,經(jīng)過0.5 h的兩相加工后,可以在硅片表面的兩相流切削顯著區(qū)觀察到一個明顯的小圓環(huán),硅片表面其它區(qū)域沒有明顯的加工痕跡。加超聲波后,由于聲場對流場的影響,小圓環(huán)的形狀有所改變,其它區(qū)域表面沒有明顯的變化。圖9a為FLUENT仿真得到的磨粒速度場云圖,仿真使用了水、磨粒兩相混合模擬,標準的湍流模型,磨粒最大速度為21 m/s。圖9b為在兩相流切削顯著區(qū)得到的小圓環(huán)照片,實際加工得到的圖形和仿真得到的非常相似。圖9c為小圓環(huán)內(nèi)表面放大500倍的形貌照片,黃色的SiO2層被沖出最大5m左右的孔洞,底部可見黑色的硅,在這樣的速度下,150m左右的氧化鋁可以去除硅片表面的SiO2。

    2)觀察到超聲空化作用的顯著影響:使用1 000目,磨粒直徑13m左右的氧化鋁,仿真顯示磨粒速度沒有大的變化,由于磨粒質(zhì)量減小為原來的1/10,經(jīng)過1 h的兩相加工后,整個硅片表面沒有明顯的加工痕跡。加超聲波0.5 h后,由于氣泡潰滅對磨粒的推動作用,在超聲空化區(qū)出現(xiàn)一個大圓環(huán)。圖9d為在超聲空化區(qū)得到的大圓環(huán)照片。圖9e為兩相流切削顯著區(qū)表面放大500倍的形貌照片,全部顯黃色,基本沒有表面材料的去除。圖9f為超聲空化區(qū)表面放大500倍的形貌照片,有SiO2層被去除,底部露出黑色的硅。

    圖9 空化驗證試驗結(jié)果

    3)沒有觀察到超聲空化作用的顯著影響:使用2 500目,磨粒直徑5m左右的氧化鋁,由于磨粒尺寸再度減小,經(jīng)過兩相沖刷后,整個硅片同樣表面沒有明顯的加工痕跡。加超聲波半小時后,即使有氣泡潰滅對磨粒的推動作用,在超聲空化區(qū)出現(xiàn)的大圓環(huán)也不明顯。圖9g為原始照片表面放大500倍的形貌照片。圖9h為兩相流切削顯著區(qū)表面放大500倍的形貌照片,基本沒有表面材料的去除。圖9i為超聲空化區(qū)表面放大500倍的形貌照片,有SiO2層被去除的痕跡,但底部沒有露出黑色的硅。

    4)沒有觀察到超聲空化作用的顯著影響:使用3 000目,磨粒直徑2m左右的氧化鋁,由于磨粒尺寸太小,即使有氣泡潰滅的作用,磨粒速度還是不夠,無論有沒有加超聲波,硅片表面都沒有明顯的材料去除的痕跡。

    3.3 試驗結(jié)果分析

    由于拋光主要是靠磨粒的切削完成的,和切削有關(guān)的只是靠近被加工表面的那小部分,會和壁面碰撞的磨粒的運動狀況,遠處的流場情況對切削影響不大。特別需要強調(diào)的是,壁面處的靜壓和切削效果無關(guān),只要仿真得到被加工表面的流體速度(或湍動能),結(jié)合磨粒的大小,根據(jù)試驗數(shù)據(jù)就可以預(yù)期達到某個表面粗糙度需要的時間。可以采用不同壓力流量的泵,設(shè)計流道能讓被加工表面的速度,湍動能一樣,加工效率差別就不大。

    上述4組試驗,也驗證了磨粒流對材料的去除效率(或能力)和磨粒沖擊壁面的動能成正比,4組試驗除了使用磨粒的大小不同以外,其它條件都類似,F(xiàn)LUENT仿真結(jié)果顯示,4組試驗的磨粒速度是相近的,由于磨粒的大小不同,導(dǎo)致它們的動能差別很大,在兩相流切削顯著區(qū)(小圓環(huán)區(qū)域)大小磨粒的切削效果有非常明顯的區(qū)別。在21 m/s的速度下,小于13m的磨粒,短時間內(nèi)觀察不到明顯的切削痕跡。

    其次,試驗驗證了空化對磨粒切削能力的影響,試驗結(jié)果和理論分析能夠互相印證。

    空化氣泡對硅片壁面的作用主要表現(xiàn)為沖擊波和微射流,微射流是由緊靠壁面的氣泡所導(dǎo)致的,由于空化強度受到控制,微射流不足以對壁面產(chǎn)生切削。氣泡潰滅的沖擊波在整個流體范圍內(nèi)都存在,其壓力大小與潰滅點的距離成反比關(guān)系,磨粒受沖擊波的作用而導(dǎo)致的動能變化和其到潰滅點的距離有關(guān)。

    雖然某一顆磨粒與潰滅點的距離無法確定,但可以從統(tǒng)計學的角度來分析這個問題。假設(shè)磨粒在水中均勻分布,磨粒間的距離越小,就意味著磨粒離潰滅點的距離越近。在進行三相流試驗時,磨粒的體積濃度都是一樣的,即單位體積內(nèi)含有的磨粒體積(質(zhì)量)相同。由于小目數(shù)的顆粒直徑大,磨粒的數(shù)目就少,這意味著磨粒間的距離大,在空化沖擊波作用下的平均動能增加量就小。

    以試驗1中150m的磨粒為例,在速度21 m/s的情況下,0.5 h就會有明顯的切削效果,在空化區(qū)里磨粒的速度已經(jīng)不夠,沖擊波加速的幅度也不大,因此觀察不到切削效果。

    在試驗2中使用13m左右的磨粒,由于流速沒有大的變化,動能只有試驗1的十分之一左右,所以在兩相流切削顯著區(qū)觀察不到切削痕跡。在空化區(qū)磨粒的初始動能和應(yīng)沖擊波而增加的動能之和達到了可以去除表面材料的強度,所以可以觀察到SiO2的去除。

    在進一步減小磨粒尺寸的試驗3、4中,由于磨粒質(zhì)量減小,在兩相流切削顯著區(qū)自然觀察不到切削痕跡,在空化區(qū)雖然磨粒獲得的動能有增加,但磨粒初始動能減小,綜合效果是磨粒切削能力減弱,直至觀察不到。

    用空化的沖擊波機理,可以比較合理地解釋試驗中出現(xiàn)的現(xiàn)象。

    3.4 方案優(yōu)化后的拋光試驗

    上述試驗只是驗證了空化的效能卻沒有充分利用其作用,讓兩相流切削顯著區(qū)和超聲空化區(qū)重疊,才能得到顯著的效果。由于超聲空化區(qū)的位置難以變動,只有增大變幅桿內(nèi)孔的尺寸才能實現(xiàn)兩個區(qū)域的重合,但這樣會造成壁厚太薄而影響振動。由于超聲波可以穿過固體進入水中,采用未經(jīng)擴孔的變幅桿從硅片背面(非加工面)發(fā)射超聲波,同樣可以在硅片表面的流體中產(chǎn)生空化作用,再用內(nèi)孔增大的變幅桿向硅片正面(加工面)噴射磨粒流,觀察到了超聲空化對磨粒流切削的顯著影響。

    圖10是試驗結(jié)果照片,圖10a是未經(jīng)研磨拋光的硅片,測得原始表面粗糙度為Ra 0.33m。圖10b是用磨粒流加工20 min后的硅片表面,磨粒大小為80目,直徑180m左右。圖10c是同樣試驗條件下加入超聲波的情況,考慮到減小超聲波功率,空化的效果會降低,過大的功率又會對表面造成傷害,其功率為根據(jù)不同強度試驗后選定的一個較為合理的值。通過換能器上的電壓、電流、功率因數(shù)和發(fā)波面積,計算得超聲強度為92 W/cm2。

    圖10 磨粒直徑180 μm的黑剛玉沖刷硅片表面的試驗結(jié)果

    圖10d是原始硅片表面放大500倍的照片(Ra 0.33m)。圖10e是兩相流加工得到的硅片表面,顯微鏡下可見切削紋理,切削區(qū)表面粗糙度為Ra 0.20m。圖10f是加入超聲波后三相流加工得到的硅片表面,切削紋理消失,切削區(qū)表面粗糙度為Ra 0.17m。超聲作用下的磨粒流加工,效率顯著增加,表面質(zhì)量也有所提高。

    4 結(jié) 論

    超聲波作用下的液體會發(fā)生氣泡潰滅的空化作用,雖然氣泡在壁面的潰滅對壁面會產(chǎn)生一定的破壞作用,但是如果可以控制氣泡潰滅的強烈程度,則可以利用這一過程增強磨粒流中磨粒對壁面的微切削的效率。

    合理強度下的超聲空化,通過氣泡潰滅時在水中產(chǎn)生的沖擊波對磨粒運動的影響增強加工效能。假設(shè)磨粒和潰滅的氣泡在流體中均勻分布,潰滅的沖擊波雖不能增加磨粒的平均速度,但能改變磨粒運動方向,使其運動更加紊亂,有部分磨粒會得到加速,從而增強磨粒流的切削能力。

    理論分析和試驗表明,在一定的試驗條件下,(本文的試驗條件是:超聲頻率20 kHz,強度92 W/cm2,磨粒尺寸180m,初速度21 m/s),磨粒流的切削效率可以得到顯著的提高,并且有改善加工質(zhì)量的效果。要讓氣泡潰滅起作用需要一定的前提條件,磨粒的初始動能和沖擊加速的動能之和才是最終起作用的決定因素。

    超聲空化不是在任何情況下都起明顯的作用,針對不同的材料,要在綜合考慮磨粒大小,初速度,超聲波的頻率強度等因素,在合理設(shè)計的試驗條件下,才能顯著提高磨粒流去除表面材料的效能。

    [1] Li Chen, Ji Shiming, Tan Dapeng. Softness abrasive flow method oriented to tiny scale mold structural surface[J]. Int J Adv Manuf Technol, 2012, 61: 975-987.

    [2] 計時鳴,鐘佳奇,譚大鵬,等. 結(jié)構(gòu)化流道環(huán)境下不同顆粒濃度的磨粒群分布及其動力學特性[J]. 農(nóng)業(yè)工程學報,2012,28(4):45-53.Ji Shiming, Zhong Jiaqi, Tan Dapeng, et al. Distribution and dynamic characteristic of particle group with different concentration in structural flow passage[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(4): 45–53 (in Chinese with English abstract)

    [3] 計時鳴,付有志,譚大鵬.軟性磨粒流雙入口約束流場數(shù)值分析及加工試驗研究[J]. 機械工程學報,2012,48(19):177–185. Ji Shiming, Fu Youzhi, Tan Dapeng. Numerical analysis and processing experiment of double-inlet restraint flow field in the soft abrasive flow machining[J]. Journal of Mechanical Engineering, 2012, 48(19): 177-185 (in Chinese with English abstract)

    [4] 計時鳴,邱毅,蔡姚杰,等. 軟性磨粒流超聲強化機理及試驗研究[J]. 機械工程學報,2014,50(7):84-93.Ji Shiming, Qiu Yi, Cai Yaojie, et al. Research on mechanism of ultrasound enhancing and the experiment based on softness abrasive flow[J]. Journal of Mechanical Engineering, 2014, 50(7): 84-93 (in Chinese with English abstract)

    [5] Maurizio P, Cristian M, Alfredo S. Characterization of near-wall accumulation regions for inertial particles in turbulent boundary layers[J]. Physics of Fluids, 2005, 17(9): 1-4.

    [6] 吳曉霞,張華余,馬空軍. 超聲空化泡運動特性的研究進展[J]. 應(yīng)用聲學,2012,31(6):416-422.Wu Xiaoxia, Zhang Huayu, Ma Kongjun. Research progresses of movement characteristics of ultrasonic cavitation bubbles[J]. Applied Acoustics, 2012, 31(6): 416-422. (in Chinese with English abstract)

    [7] 崔方玲,紀威.超聲空化氣泡動力學仿真及其影響因素分析[J]. 農(nóng)業(yè)工程學報,2013,29(17):24-29. Cui Fangling, Ji Wei. Dynamic simulation of ultrasonic cavitation bubble and analysis of its influencing factors[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(17): 24-29 (in Chinese with English abstract)

    [8] 王含. 氣泡行為的數(shù)值研究[D]. 上海:復(fù)旦大學,2011. Wang Han. Numerical Investigation of Bubble Behavior [D]. Shanghai: Fudan University, 2011. (in Chinese with English abstract)

    [9] 張鵬利,林書玉,張濤. 空化泡液體外圍壓強的分布[J]. 應(yīng)用聲學,2012,31(2):67-70. Zhang Pengli, Lin Shuyu, Zhang Tao. Studies of the external pressure distribution of cavitation bubble[J]. Applied Acoustics, 2012, 31(2): 67-70. (in Chinese with English abstract)

    [10] 倪漢根,何子干.空泡潰滅沖擊波對流場攜帶顆粒的作用[J]. 大連理工大學學報, 1994,34(6):727-731. Ni Hangen, He Zigan. Impact of shock waves of bubble collapse on solid particles carried by flow fields[J]. Journal of Dalian University of Technology, 1994, 34(6): 727-731. (in Chinese with English abstract)

    [11] 劉海軍,安宇. 空化單氣泡外圍壓強分布[J]. 物理學報,2004,53(5):1406-1412. Liu Haijun, An Yu. Pressure distribution out side a single cavitation bubble[J]. Acta Phys. Sin, 2004, 53(5): 1406-1412. (in Chinese with English abstract)

    [12] Hashish M. A Model for abrasive-waterjet (AWJ) machining[J]. Journal of Engineering Materials and Technology ASME, 1989, 111(2): 154-162.

    [13] 楊曉京,陳子辰,樊瑜瑾,等. 磨粒磨損中微觀切削過程分子動力學模擬[J]. 農(nóng)業(yè)機械學報,2007,38(5):161-164. Yang Xiaojing, Chen Zichen, Fan Yujin, et al. Molecular dynanics simulation of micro cutting process of abrasive Wear[J]. Transactions of the Chinese Society for Agricultural Machinery, 2007, 38(5): 161-164. (in Chinese with English abstract)

    [14] Wang Xing, Zhang Yong, Zhang Feihu. Preliminary investigation on nanoparticle colloid hydrodynamic cavitation jet polishing technology[J]. Nanotechnology and Precision Engineering, 2011, 9(6): 483-487.

    [15] 馬占龍,王君林.微射流拋光機理仿真及實驗研究[J].光電工程,2012,39(5): 139-144. Ma Zhanlong, Wang Junlin. Simulation of mechanism and experimental study in micro fluid jet Polishing[J]. Opto- Electronic Engineering, 2012, 39(5): 139-144. (in Chinese with English abstract)

    [16] 王星. 納米膠體空化射流拋光及其關(guān)鍵技術(shù)研究[D]. 哈爾濱:哈爾濱工業(yè)大學,2014. Wang Xing. Study on the Key Techology of Nano Particle Colloid Hydrodynamic Cavitation Jet Polishing[D]. Harbin: Harbin Institute of Technology, 2014. (in Chinese with English abstract)

    [17] Rayleigh J W. On the Pressure developed in a liquid during the collapse of a spherical cavity[J]. Philos Mag, 1970, 34(200): 221-226.

    [18] Wang Guoyu, Senocak Inanc, Shyy Wei,et al. Dynamics of attached turbulent cavitating flows[J]. Progress in Aerospace Sciences, 2001, 37(6): 551-581.

    [19] 龔烈航,陸國勝,張曉鑫,等.三相磨料射流作用下材料的破壞機理研究[J]. 摩擦學學報,2002,22(4):283-286.Gong Liehang, Lu Guosheng, Zhang Xiaoxin, et al. Study on the damage mechanism under the three-phases abrasive jet[J]. Tribology, 2002, 22(4): 283-286 (in Chinese with English abstract)

    [20] Vyas B, Hansson L H. The cavitation erosion-corrosion of stainless steel[J]. Corrosion Science, 1990, 30(8): 761-770.

    [21] 宮雪,韓杰勝,李珍,等. 標準超聲振動氣蝕孕育期內(nèi)材料響應(yīng)與表面形貌[J]. 摩擦學學報,2013,33(6):550-555. Gong Xue, Han Jiesheng, Li Zhen, et al. Material response and surface morphology in incubation period of a standard ultrasonic vibration cavitation erosion test[J]. Tribology, 2013, 33(6): 550-555 (in Chinese with English abstract)

    [22] 王勇,林書玉,莫潤陽,等. 含氣泡液體中氣泡振動的研究[J]. 物理學報,2013,62(13):134304_1-134304_5. Wang Yong, Lin Shuyu, Mo Runyang, et al. Vibration of the bubble in bubbly liquids[J]. Acta Phys. Sin, 2013, 62(13): 134304_1-134304_5.(in Chinese with English abstract)

    [23] 陳逢軍,唐宇,苗想亮,等. 磨料射流表面拋光研究綜述[J]. 表面技術(shù),2015,44(11):119-127.Chen Fengjun, Tangyu, Miao Xiangliang, et al. Review on the abrasive jet surface polishing (AJP) technology[J]. Suface Technology, 2015, 44(11): 119-127. (in Chinese with English abstract)

    [24] 袁巨龍,吳喆,呂冰海,等. 非球面超精密拋光技術(shù)研究現(xiàn)狀[J]. 機械工程學報,2012,48(23):168-170. Yuan Julong, Wu Zhe, Lv Binghai, et al. Review on ultra-precision polishing technology of aspheric surface[J]. Journal of Mechanical Engineering, 2012, 48(23): 168-170. (in Chinese with English abstract)

    [25] Ke J H. Abrasive jet polishing of micro-channels using compound sic abrasives with compound additives[J]. Advanced Materials Research, 2013, 67(5/8): 1151-1159..

    [26] Reyssat E, Vandewalle N. Aging of an antibubble[J]. EPL (Europhysics Letters), 2005, 69(6): 966-970.

    [27] 張成光,張勇,張飛虎,等. 磨料水射流加工去除模型研究[J]. 機械工程學報,2015,51(7):189-194.Zhang Chengguang, Zhang Yong, Zhang Feihu, et al. Study on removal model of abrasive water jet machining[J]. Journal of Mechanical Engineering, 2015, 51(7): 189-194. (in Chinese with English abstract)

    [28] Beaucamp A, Namba Y. Super-smooth finishing of diamond turned hard X-ray molding dies by combined fluid jet and bonnet polishing[J]. CIRP Annals-Manufacturing Technology, 2013, 62(1): 315-318.

    [29] Dorbolo S, Caps H, Vandewalle N. Fluid instabilities in the birth and death of antibubbles[J]. New Journal of Physics, 2003, 5(161): 1-161.

    [30] Zou J, Ji C, Yuan G B, et al. Collapse of an antibubble[J]. Physical Review E, 2013, 87(6): 061002.

    [31] Anthony B, Yoshiharu N, Wilhelmus M. Surface Integrity of fluid jet polished tungsten carbiede[J]. Procedia CIRP, 2014, 13: 377-381.

    [32] 李宜燃. 軟性磨粒流拋光超聲波湍流強化發(fā)展過程及實驗研究[D]. 浙江:浙江工業(yè)大學,2012. Li Yiran. Ultrasonic Turbulence Strengthen and Expreimental Studies Based on The Sofe Abrasive Flow Machining[D]. Zhejiang: Zhejiang University of Technology, 2012.(in Chinese with English abstract)

    Effect of ultrasonic cavitation on maching efficiency and quality of soft abrasive flow

    Ji Shiming1, Chen Kai1,2※, Tan Dapeng1, Zheng Chenliang1

    (1.310014,; 2.323000,)

    The method of material removal by ultrasonic cavitation was studied to improve the efficiency of soft abrasive flow machining. Soft abrasive flow machining uses water as the carrier, and the actual processing effect is achieved by abrasive. The uneven peak and valley structure of the machined surface is impacted by a large number of abrasive particles. Those molecules at the high positions of the workpiece surface are more likely to slip or fall apart after collision. First of all, because of the probability of its impact, the molecule of the top in the collision will get more energy to escape from the original position. After a large amount of impact, the structure of the peak and valley tends to be smooth and the surface roughness decreases. The processing efficiency is determined by the kinetic energy of abrasive particles impinging on the wall of the workpiece. Based on the analysis of the abrasive particle polishing experiments and computational fluid dynamics simulation, it was proposed and verified that the kinetic energy of abrasive grain impacting the wall was the main factor affecting the cutting efficiency. When the ultrasonic wave passes through the water, the liquid is no longer a kind of incompressible fluid, and the density of particles varies with sound waves. Bubbles will appear in the water, and the size of bubble will be changed by the ultrasonic wave. The modified Rayleigh-plesset equation of motion can be used to calculate the variation of bubble diameter with time. Ultrasonic cavitation occurs when the bubble diameter is reduced to zero. The actual observation and theoretical calculation indicate that the collapse of the bubble in its interior will produce high temperature and pressure, which triggers a series of complex physical and chemical phenomena. According to the study of ultrasonic cavitation bubble collapse, and the research results of the peripheral water shock wave formation and its propagation, the role of bubble collapse in the change of abrasive particle kinetic energy was analyzed. When the bubble collapses, high-pressure area formed on the periphery of bubbles is thin, and the high-pressure area then spreads outward to form a shock wave. The amplitude of shock wave is attenuated with a speed of the reciprocal of propagation distance. The kinetic energy of the abrasive particles in the vicinity of the bubble due to shock wave can be analyzed. Strong ultrasonic wave can destroy the surface of parts. In order to make the cavitation bubble collapse caused by the "shock wave" and "micro jet" reasonably affect the machining process, the cavitation intensity must be controlled reasonably, so as to promote the acceleration of the abrasive particles as much as possible without hurting the machined surface. A large number of experiments have been carried out to obtain a reasonable ultrasonic intensity. The experimental device was set up to observe the effect of ultrasonic cavitation on the removal of SiO2on the surface of silicon wafer. Significant effect of ultrasonic cavitation was observed in the process of abrasive particle flow machining on the silicon wafer surface. The results show that the ultrasonic cavitation can change the kinetic energy of the abrasive particles in the particle flow, which can improve the cutting efficiency and improve the surface quality. When the frequency of ultrasonic wave is 20 kHz, the ultrasonic intensity is 92 W/cm2, the abrasive particle velocity reaches 21 m/s, and the size of the abrasive particles is above 13m, it can significantly improve the abrasive flow to remove the surface material under the action of ultrasonic wave.

    ultrasonics; cavitation; machining; soft abrasive flow; three phase flow; abrasive grain movement

    10.11975/j.issn.1002-6819.2017.12.011

    TG580.1

    A

    1002-6819(2017)-12-0082-09

    2016-12-01

    2017-05-14

    國家自然科學基金資助項目(51575494;51375446);浙江省自然科學基金重點項目(LZ14E050001);浙江省杰出青年科學基金項目(LR16E050001)

    計時鳴,男,浙江溫州人,浙江工業(yè)大學機械工程學院教授,博士,博士生導(dǎo)師,中國機械工程學會高級會員、中國農(nóng)業(yè)機械學會高級會員、IEEE會員。杭州浙江工業(yè)大學機械學院,310014。 Email:jishiming@zjut.edu.cn

    陳 凱,男,浙江麗水人,浙江麗水學院工學院,博士生,高級工程師。Email:ck207701@163.com

    計時鳴,陳 凱,譚大鵬,鄭晨亮. 超聲空化對軟性磨粒流切削效率和質(zhì)量的影響[J]. 農(nóng)業(yè)工程學報,2017,33(12):82-90. doi:10.11975/j.issn.1002-6819.2017.12.011 http://www.tcsae.org

    Ji Shiming, Chen Kai, Tan Dapeng, Zheng Chenliang. Effect of ultrasonic cavitation on maching efficiency and quality of soft abrasive flow[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 82-90. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.12.011 http://www.tcsae.org

    猜你喜歡
    硅片磨粒空化
    功率超聲作用下鋼液中空化泡尺寸的演變特性
    鋼鐵釩鈦(2023年5期)2023-11-17 08:48:34
    基于凸多面體碰撞檢測的虛擬砂輪建模研究
    超精密表面研拋磨粒的研究進展
    硅酸鹽通報(2021年3期)2021-04-18 11:01:42
    金剛石砂輪磨削貼膜硅片崩邊的研究*
    單個鐵氧體磨粒尺寸檢測電磁仿真
    光伏:硅片市場條件嚴峻 高效單晶需求回暖
    能源(2018年10期)2018-01-16 02:37:47
    微晶剛玉磨粒磨削20CrMnTi鋼的數(shù)值模擬研究
    三維扭曲水翼空化現(xiàn)象CFD模擬
    不同運動形式下水物相互作用空化數(shù)值模擬
    用于硅片檢測的精密直線運動平臺的設(shè)計與實現(xiàn)
    久久久久国内视频| 国产精品乱码一区二三区的特点| 麻豆久久精品国产亚洲av| 啦啦啦免费观看视频1| 国产欧美日韩精品亚洲av| 中国美女看黄片| 精品日产1卡2卡| 波多野结衣巨乳人妻| 在线观看www视频免费| 久久精品aⅴ一区二区三区四区| 午夜福利免费观看在线| 后天国语完整版免费观看| 欧美另类亚洲清纯唯美| 亚洲色图av天堂| 波多野结衣av一区二区av| 美女高潮到喷水免费观看| 午夜福利视频1000在线观看| 制服人妻中文乱码| 午夜两性在线视频| 欧美黄色片欧美黄色片| 18美女黄网站色大片免费观看| 老司机午夜福利在线观看视频| 色播在线永久视频| 日韩免费av在线播放| 岛国视频午夜一区免费看| 久久精品91无色码中文字幕| 国产成人啪精品午夜网站| 老司机深夜福利视频在线观看| 男人舔奶头视频| 一二三四社区在线视频社区8| 亚洲精品久久国产高清桃花| 国产又黄又爽又无遮挡在线| 一区二区三区激情视频| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影| 欧美色欧美亚洲另类二区| 看片在线看免费视频| 老汉色∧v一级毛片| 日韩欧美三级三区| 成人手机av| 日韩欧美在线二视频| 色婷婷久久久亚洲欧美| 国产高清激情床上av| 亚洲一码二码三码区别大吗| 久久伊人香网站| www日本在线高清视频| 无遮挡黄片免费观看| 亚洲av日韩精品久久久久久密| 久9热在线精品视频| 久久人妻av系列| 女同久久另类99精品国产91| 国产精品,欧美在线| 色精品久久人妻99蜜桃| 观看免费一级毛片| 日本黄色视频三级网站网址| 啦啦啦 在线观看视频| 亚洲色图 男人天堂 中文字幕| 观看免费一级毛片| 精品电影一区二区在线| 99国产综合亚洲精品| 免费在线观看亚洲国产| 亚洲国产精品合色在线| 国产免费av片在线观看野外av| 欧美在线一区亚洲| 中文在线观看免费www的网站 | 黑人巨大精品欧美一区二区mp4| 一区福利在线观看| 最近最新免费中文字幕在线| 国产aⅴ精品一区二区三区波| 麻豆成人av在线观看| 久久中文看片网| a在线观看视频网站| 久久午夜综合久久蜜桃| 国产精品亚洲一级av第二区| 天堂√8在线中文| 最近在线观看免费完整版| 亚洲在线自拍视频| 亚洲一区二区三区色噜噜| 亚洲国产欧洲综合997久久, | 国产男靠女视频免费网站| 亚洲av日韩精品久久久久久密| 免费在线观看亚洲国产| 成年版毛片免费区| 欧美午夜高清在线| 久久久国产成人免费| 国产精品电影一区二区三区| 2021天堂中文幕一二区在线观 | 一级黄色大片毛片| 亚洲熟妇熟女久久| xxxwww97欧美| 欧美不卡视频在线免费观看 | a级毛片在线看网站| 亚洲av电影不卡..在线观看| 亚洲午夜精品一区,二区,三区| 一区二区日韩欧美中文字幕| 欧美成狂野欧美在线观看| 人人澡人人妻人| 动漫黄色视频在线观看| 国产一区二区在线av高清观看| 中文字幕人妻熟女乱码| 亚洲国产欧美网| 国产野战对白在线观看| 日韩高清综合在线| 免费一级毛片在线播放高清视频| 啦啦啦免费观看视频1| 91麻豆精品激情在线观看国产| 一区二区三区精品91| 黄色a级毛片大全视频| 亚洲国产欧美日韩在线播放| 黑丝袜美女国产一区| 精品无人区乱码1区二区| 在线观看舔阴道视频| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利在线观看吧| 在线播放国产精品三级| 亚洲专区中文字幕在线| 狠狠狠狠99中文字幕| 成熟少妇高潮喷水视频| 久久久久久免费高清国产稀缺| 老汉色av国产亚洲站长工具| 十分钟在线观看高清视频www| 中文字幕精品免费在线观看视频| 亚洲精品色激情综合| 欧美色视频一区免费| 国产精品电影一区二区三区| 嫩草影视91久久| 成在线人永久免费视频| 日本三级黄在线观看| 中文亚洲av片在线观看爽| 深夜精品福利| 最新在线观看一区二区三区| 亚洲 欧美 日韩 在线 免费| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕精品亚洲无线码一区 | 狂野欧美激情性xxxx| 日本成人三级电影网站| 亚洲欧美精品综合一区二区三区| 在线观看免费视频日本深夜| 变态另类成人亚洲欧美熟女| 久久久久久免费高清国产稀缺| 久久国产精品男人的天堂亚洲| 免费av毛片视频| 久久久国产精品麻豆| 老汉色av国产亚洲站长工具| 亚洲黑人精品在线| 国产欧美日韩一区二区精品| 免费看日本二区| 啦啦啦韩国在线观看视频| 国产免费男女视频| 欧美日韩黄片免| 搡老岳熟女国产| 国产午夜精品久久久久久| 久久久久国内视频| 国产97色在线日韩免费| 中出人妻视频一区二区| 日本黄色视频三级网站网址| 亚洲成人国产一区在线观看| 国产高清有码在线观看视频 | 99国产综合亚洲精品| 麻豆国产av国片精品| 久久精品国产99精品国产亚洲性色| 久久热在线av| 窝窝影院91人妻| 999久久久精品免费观看国产| 一级a爱片免费观看的视频| 深夜精品福利| 黄片小视频在线播放| 51午夜福利影视在线观看| 夜夜夜夜夜久久久久| 激情在线观看视频在线高清| 亚洲熟妇中文字幕五十中出| 性欧美人与动物交配| 一级a爱视频在线免费观看| 亚洲精品美女久久久久99蜜臀| 美女 人体艺术 gogo| 老司机深夜福利视频在线观看| 国产私拍福利视频在线观看| 亚洲成国产人片在线观看| 日本熟妇午夜| 久久久国产精品麻豆| 国产精品久久电影中文字幕| 91字幕亚洲| 免费人成视频x8x8入口观看| 成人欧美大片| 男女做爰动态图高潮gif福利片| 午夜精品久久久久久毛片777| 日韩国内少妇激情av| 岛国视频午夜一区免费看| 老司机午夜十八禁免费视频| 日韩中文字幕欧美一区二区| 久久久国产精品麻豆| 亚洲男人的天堂狠狠| 精品一区二区三区视频在线观看免费| 亚洲性夜色夜夜综合| 美女扒开内裤让男人捅视频| 好男人电影高清在线观看| 免费人成视频x8x8入口观看| 91在线观看av| 深夜精品福利| 日韩欧美在线二视频| 日韩免费av在线播放| 亚洲黑人精品在线| 一本综合久久免费| 国产一区二区在线av高清观看| 特大巨黑吊av在线直播 | 可以在线观看的亚洲视频| 老司机深夜福利视频在线观看| 88av欧美| 国产激情久久老熟女| 丝袜在线中文字幕| 欧美成人一区二区免费高清观看 | 日本一本二区三区精品| 精华霜和精华液先用哪个| 亚洲精华国产精华精| 成人欧美大片| 亚洲,欧美精品.| 欧美色欧美亚洲另类二区| 在线观看免费视频日本深夜| 满18在线观看网站| 一级毛片女人18水好多| 久久精品aⅴ一区二区三区四区| 女警被强在线播放| 熟女电影av网| 嫩草影院精品99| 制服人妻中文乱码| 日本 av在线| x7x7x7水蜜桃| 国产精品综合久久久久久久免费| 久久精品成人免费网站| 欧美性长视频在线观看| 成在线人永久免费视频| 男女那种视频在线观看| 日韩精品免费视频一区二区三区| aaaaa片日本免费| 又紧又爽又黄一区二区| 51午夜福利影视在线观看| 99精品在免费线老司机午夜| 村上凉子中文字幕在线| 亚洲aⅴ乱码一区二区在线播放 | 国产精品99久久99久久久不卡| 日日摸夜夜添夜夜添小说| 此物有八面人人有两片| 一本久久中文字幕| 婷婷精品国产亚洲av| 19禁男女啪啪无遮挡网站| 日本在线视频免费播放| 亚洲国产毛片av蜜桃av| 欧美色欧美亚洲另类二区| 国产午夜福利久久久久久| 国产精品精品国产色婷婷| 国产又黄又爽又无遮挡在线| 久久婷婷人人爽人人干人人爱| 久久香蕉激情| 亚洲成人免费电影在线观看| 免费女性裸体啪啪无遮挡网站| 国产成人av激情在线播放| 男人舔女人下体高潮全视频| 免费看日本二区| 精品电影一区二区在线| 精品欧美国产一区二区三| 女性被躁到高潮视频| 欧美性猛交黑人性爽| 黄色丝袜av网址大全| 2021天堂中文幕一二区在线观 | 日韩成人在线观看一区二区三区| 国产色视频综合| 老熟妇乱子伦视频在线观看| 欧美日韩乱码在线| 搡老熟女国产l中国老女人| 黄片小视频在线播放| 亚洲av电影不卡..在线观看| 国产在线观看jvid| 夜夜夜夜夜久久久久| 少妇粗大呻吟视频| 男女视频在线观看网站免费 | 波多野结衣高清无吗| 18禁裸乳无遮挡免费网站照片 | 国产欧美日韩精品亚洲av| 日本在线视频免费播放| www日本黄色视频网| 免费在线观看黄色视频的| 又黄又粗又硬又大视频| 侵犯人妻中文字幕一二三四区| 国产色视频综合| 成人手机av| 特大巨黑吊av在线直播 | 午夜福利高清视频| 国产三级黄色录像| 免费在线观看亚洲国产| www日本在线高清视频| 午夜免费观看网址| 好男人在线观看高清免费视频 | 久久久久久久午夜电影| 免费在线观看成人毛片| 又黄又粗又硬又大视频| 国产免费av片在线观看野外av| 精品欧美国产一区二区三| 黑人操中国人逼视频| 在线观看一区二区三区| 久久精品成人免费网站| 窝窝影院91人妻| 在线观看一区二区三区| 日韩精品免费视频一区二区三区| 特大巨黑吊av在线直播 | 国内久久婷婷六月综合欲色啪| 亚洲 国产 在线| 听说在线观看完整版免费高清| 日韩欧美 国产精品| 神马国产精品三级电影在线观看 | 精品第一国产精品| 国产真实乱freesex| 午夜成年电影在线免费观看| 国产伦一二天堂av在线观看| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久国产精品人妻aⅴ院| 久久精品国产清高在天天线| 国产熟女午夜一区二区三区| 国产激情欧美一区二区| 麻豆久久精品国产亚洲av| 女生性感内裤真人,穿戴方法视频| 国产97色在线日韩免费| 动漫黄色视频在线观看| 极品教师在线免费播放| aaaaa片日本免费| 变态另类成人亚洲欧美熟女| 国产精品九九99| 人妻久久中文字幕网| 亚洲 欧美 日韩 在线 免费| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 国产精品久久久久久亚洲av鲁大| 欧美日韩亚洲综合一区二区三区_| 国产成人一区二区三区免费视频网站| 亚洲专区字幕在线| av在线播放免费不卡| 久久久久国内视频| 可以在线观看毛片的网站| 国产又爽黄色视频| 在线观看免费午夜福利视频| 精品福利观看| 亚洲无线在线观看| 最新美女视频免费是黄的| 美女免费视频网站| 日韩有码中文字幕| 国产精品国产高清国产av| 国产精品九九99| 一本精品99久久精品77| 国产成人欧美| 91麻豆精品激情在线观看国产| 久久亚洲真实| 精品久久久久久久久久免费视频| 97碰自拍视频| 国产精华一区二区三区| 国产精品乱码一区二三区的特点| 国产欧美日韩精品亚洲av| 法律面前人人平等表现在哪些方面| av天堂在线播放| 欧美成狂野欧美在线观看| 欧美乱色亚洲激情| 亚洲,欧美精品.| 国产激情久久老熟女| av中文乱码字幕在线| 国产精品 国内视频| 最新美女视频免费是黄的| 18禁黄网站禁片免费观看直播| 亚洲九九香蕉| 色av中文字幕| 久久久久久人人人人人| 在线观看午夜福利视频| 欧美乱色亚洲激情| 亚洲中文av在线| 色综合欧美亚洲国产小说| 变态另类成人亚洲欧美熟女| 女生性感内裤真人,穿戴方法视频| 久久草成人影院| 亚洲av成人一区二区三| 一a级毛片在线观看| 欧美日韩乱码在线| 伦理电影免费视频| 欧美三级亚洲精品| АⅤ资源中文在线天堂| 中文字幕精品亚洲无线码一区 | 制服丝袜大香蕉在线| 可以免费在线观看a视频的电影网站| 午夜两性在线视频| 国产精品日韩av在线免费观看| 国产精品美女特级片免费视频播放器 | 999精品在线视频| 亚洲精品美女久久av网站| 亚洲成a人片在线一区二区| 亚洲精品在线美女| 香蕉丝袜av| 国产精品一区二区免费欧美| 精华霜和精华液先用哪个| 免费在线观看黄色视频的| 黄网站色视频无遮挡免费观看| 久久中文字幕人妻熟女| 高潮久久久久久久久久久不卡| avwww免费| 亚洲成a人片在线一区二区| 999久久久精品免费观看国产| 免费看美女性在线毛片视频| 免费观看精品视频网站| 国产一区在线观看成人免费| 非洲黑人性xxxx精品又粗又长| 在线观看www视频免费| 亚洲七黄色美女视频| 可以在线观看毛片的网站| 免费av毛片视频| 午夜福利成人在线免费观看| 精品电影一区二区在线| 又大又爽又粗| 成人永久免费在线观看视频| 久久中文字幕人妻熟女| 激情在线观看视频在线高清| 黄片大片在线免费观看| 免费在线观看影片大全网站| 国产成人av教育| 午夜视频精品福利| 久久久久久久精品吃奶| 好看av亚洲va欧美ⅴa在| 国内揄拍国产精品人妻在线 | 成人三级做爰电影| 国产亚洲欧美在线一区二区| 欧美色视频一区免费| 侵犯人妻中文字幕一二三四区| 国产亚洲精品久久久久5区| 免费一级毛片在线播放高清视频| 国产精品久久久人人做人人爽| 国产97色在线日韩免费| 动漫黄色视频在线观看| 日本精品一区二区三区蜜桃| 在线观看www视频免费| 日韩欧美 国产精品| 国产精品久久视频播放| 757午夜福利合集在线观看| 亚洲av电影在线进入| 日本三级黄在线观看| www.自偷自拍.com| av超薄肉色丝袜交足视频| 热re99久久国产66热| 国产91精品成人一区二区三区| 国产精品美女特级片免费视频播放器 | 91在线观看av| 国产私拍福利视频在线观看| 狂野欧美激情性xxxx| 亚洲 国产 在线| 欧美色视频一区免费| 日本熟妇午夜| 色综合站精品国产| 精品第一国产精品| 美女午夜性视频免费| 99在线视频只有这里精品首页| 日韩精品中文字幕看吧| av有码第一页| 最新美女视频免费是黄的| 久久精品影院6| 久久精品国产综合久久久| 成人午夜高清在线视频 | 一区二区三区精品91| 999久久久国产精品视频| 日本撒尿小便嘘嘘汇集6| 美女大奶头视频| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美三级三区| 亚洲av片天天在线观看| 亚洲av熟女| 成人特级黄色片久久久久久久| 操出白浆在线播放| 久久国产精品男人的天堂亚洲| 一本精品99久久精品77| 亚洲av日韩精品久久久久久密| 欧美 亚洲 国产 日韩一| 18美女黄网站色大片免费观看| 悠悠久久av| 99久久精品国产亚洲精品| 国产精品香港三级国产av潘金莲| 91老司机精品| 欧美性猛交黑人性爽| 欧美精品啪啪一区二区三区| 给我免费播放毛片高清在线观看| 日本五十路高清| 特大巨黑吊av在线直播 | 欧美日本视频| 色精品久久人妻99蜜桃| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲国产一区二区在线观看| 男女床上黄色一级片免费看| 国产伦在线观看视频一区| 18美女黄网站色大片免费观看| 国产高清视频在线播放一区| 老司机午夜福利在线观看视频| 1024香蕉在线观看| 亚洲国产欧美网| 99久久无色码亚洲精品果冻| 不卡一级毛片| 最新在线观看一区二区三区| 宅男免费午夜| 在线观看舔阴道视频| 中亚洲国语对白在线视频| 国产爱豆传媒在线观看 | 色老头精品视频在线观看| 两性夫妻黄色片| 人人妻人人澡欧美一区二区| 天堂影院成人在线观看| av有码第一页| 老司机在亚洲福利影院| 久99久视频精品免费| 欧美av亚洲av综合av国产av| 亚洲天堂国产精品一区在线| 欧美性猛交黑人性爽| 白带黄色成豆腐渣| 最近在线观看免费完整版| 成人三级黄色视频| 我的亚洲天堂| 日本a在线网址| 黄网站色视频无遮挡免费观看| aaaaa片日本免费| 国产精品久久久人人做人人爽| 亚洲av中文字字幕乱码综合 | 亚洲精品在线美女| 久久久国产成人精品二区| 精品乱码久久久久久99久播| 久久精品国产综合久久久| 中国美女看黄片| 男人舔奶头视频| 欧美国产日韩亚洲一区| 91成年电影在线观看| 成年版毛片免费区| bbb黄色大片| 国产精品香港三级国产av潘金莲| 亚洲精品av麻豆狂野| 亚洲欧美一区二区三区黑人| 黄色 视频免费看| 久久久国产成人精品二区| 深夜精品福利| 欧美性长视频在线观看| 久9热在线精品视频| 欧美在线黄色| 国产激情欧美一区二区| 99久久无色码亚洲精品果冻| 自线自在国产av| 色尼玛亚洲综合影院| 国产精品爽爽va在线观看网站 | 中文字幕另类日韩欧美亚洲嫩草| 久久精品91蜜桃| 久久精品国产99精品国产亚洲性色| 深夜精品福利| 欧美精品啪啪一区二区三区| 亚洲人成电影免费在线| 婷婷六月久久综合丁香| 99久久99久久久精品蜜桃| 亚洲精品久久国产高清桃花| 午夜影院日韩av| 人人澡人人妻人| 免费看美女性在线毛片视频| 欧美绝顶高潮抽搐喷水| 成人手机av| av中文乱码字幕在线| 中文在线观看免费www的网站 | www国产在线视频色| 久久久久久久久久黄片| 中文字幕精品免费在线观看视频| 99在线视频只有这里精品首页| 国产不卡一卡二| 国产激情欧美一区二区| 亚洲精华国产精华精| 婷婷精品国产亚洲av在线| 亚洲欧洲精品一区二区精品久久久| 一个人免费在线观看的高清视频| 亚洲成人免费电影在线观看| 免费看日本二区| 免费在线观看完整版高清| 日韩国内少妇激情av| 中文字幕人妻熟女乱码| 国产免费av片在线观看野外av| 在线观看舔阴道视频| 欧美乱妇无乱码| 美女免费视频网站| 久久 成人 亚洲| 亚洲精品久久国产高清桃花| 精品久久久久久,| 啦啦啦韩国在线观看视频| 亚洲精华国产精华精| 2021天堂中文幕一二区在线观 | 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区视频在线观看免费| 人妻久久中文字幕网| 成人免费观看视频高清| 亚洲熟女毛片儿| 性欧美人与动物交配| 欧美激情极品国产一区二区三区| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色| 啦啦啦 在线观看视频| 亚洲专区国产一区二区| 国产精品亚洲av一区麻豆| 国产黄a三级三级三级人| 变态另类成人亚洲欧美熟女| 国产精品亚洲av一区麻豆| 19禁男女啪啪无遮挡网站| 女同久久另类99精品国产91| 成年女人毛片免费观看观看9| 国产午夜精品久久久久久| 在线观看免费视频日本深夜| 老司机午夜福利在线观看视频| 琪琪午夜伦伦电影理论片6080| 一级黄色大片毛片| 免费看日本二区| 亚洲av日韩精品久久久久久密| 亚洲国产高清在线一区二区三 | 色在线成人网| 日韩大尺度精品在线看网址| 国产精品久久久久久人妻精品电影|