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

    非對(duì)稱表面織構(gòu)表征及動(dòng)壓潤滑效應(yīng)分析

    2022-03-31 05:48:24朱春霞王德全楊曉楠常永祥
    表面技術(shù) 2022年3期
    關(guān)鍵詞:動(dòng)壓織構(gòu)非對(duì)稱

    朱春霞,王德全,楊曉楠,常永祥

    非對(duì)稱表面織構(gòu)表征及動(dòng)壓潤滑效應(yīng)分析

    朱春霞,王德全,楊曉楠,常永祥

    (沈陽建筑大學(xué) 機(jī)械工程學(xué)院,沈陽 110168)

    建立具有非對(duì)稱表面織構(gòu)的粗糙平面的表征方程,并分析方程參數(shù)對(duì)其形貌的影響,依托表征方程重建的織構(gòu)平面研究非對(duì)稱表面織構(gòu)在動(dòng)壓潤滑條件下的潤滑特性。根據(jù)非對(duì)稱織構(gòu)平面的加工過程,基于W-M函數(shù)重構(gòu)粗糙表面,再疊加非對(duì)稱織構(gòu)形貌特征方程,實(shí)現(xiàn)不同尺寸非對(duì)稱表面織構(gòu)的數(shù)學(xué)表征,利用Matlab軟件對(duì)方程中表面織構(gòu)參數(shù)和粗糙表面參數(shù)對(duì)織構(gòu)平面形貌的影響進(jìn)行分析,并利用方程表征的織構(gòu)平面建立動(dòng)壓潤滑壓力油膜的實(shí)體模型,對(duì)壓力油膜進(jìn)行CFD計(jì)算,分別得出正向非對(duì)稱織構(gòu)、對(duì)稱織構(gòu)、反向非對(duì)稱織構(gòu)中的壓力分布規(guī)律和流線分布規(guī)律,實(shí)現(xiàn)非對(duì)稱表面織構(gòu)的動(dòng)壓潤滑效應(yīng)分析。表征方程的各個(gè)參數(shù)中,尺度系數(shù)和分形維數(shù)主要影響粗糙表面輪廓高度的變化范圍和變化頻率,微觀輪廓高度的變化范圍隨著尺度系數(shù)的減小以同樣的倍率減小,微觀輪廓高度變化隨著分形維數(shù)的逐漸增大越來越劇烈,變化頻率也越來越快,非對(duì)稱度、織構(gòu)深度等織構(gòu)相關(guān)參數(shù)主要影響織構(gòu)的形貌特征。潤滑壓力油膜CFD分析計(jì)算表明,正向非對(duì)稱織構(gòu)(=2)的進(jìn)出口壓力差為1747 Pa,明顯大于對(duì)稱織構(gòu)(=1,387 Pa)和反向非對(duì)稱織構(gòu)(=0.5,707 Pa),三者的流線分布云圖表明,正向非對(duì)稱織構(gòu)內(nèi)的回流效應(yīng)更加劇烈,流線分布更加復(fù)雜。所建立的非對(duì)稱織構(gòu)平面的表征方程可以清晰地描述其形貌特征,方程參數(shù)對(duì)形貌特征的影響機(jī)理明確,且構(gòu)型合理的非對(duì)稱織構(gòu)在提高承載能力和改善潤滑方面存在明顯的優(yōu)越性。

    非對(duì)稱;表面織構(gòu);表征方程;W-M函數(shù);CFD;動(dòng)壓潤滑效應(yīng)

    表面織構(gòu)技術(shù)是一種可以實(shí)現(xiàn)減摩、抗磨,能提高承載能力等功效的表面處理手段[1-3],隨著精加工技術(shù)的進(jìn)步和機(jī)械設(shè)備可靠性要求的提高,表面織構(gòu)技術(shù)在改善機(jī)械結(jié)合面特性方面起著越來越重要的作用[4-7]。

    在以往的表面織構(gòu)相關(guān)研究中,對(duì)于表面織構(gòu)輪廓形狀的描述往往通過二維或者三維的圖像實(shí)現(xiàn),缺少具體的表征方程,導(dǎo)致織構(gòu)動(dòng)壓潤滑研究中的幾何模型往往通過構(gòu)建理想光滑平面,這與實(shí)際加工條件下所得的結(jié)合面輪廓形貌存在較大誤差。在1991年,Majumdar和Bhushan[8]基于W-M函數(shù),創(chuàng)新性地提出了用于表征粗糙表面特征的數(shù)學(xué)模型,將粗糙結(jié)合面的表征模型問題帶離了統(tǒng)計(jì)學(xué)參數(shù)的束縛,為后續(xù)粗糙表面的多種數(shù)學(xué)表達(dá)模型提供了理論基礎(chǔ)。后來Yan和Komvopoulos[9]采用雙變量W-M函數(shù),將二維的模型延伸至三維形式,自此W-M函數(shù)逐漸成為粗糙表面表征最常用的理論基礎(chǔ)之一,這為研究非對(duì)稱織構(gòu)平面表征方程的建立指明了方向。文中表征非對(duì)稱織構(gòu)平面無織構(gòu)區(qū)域采用基于W-M函數(shù)的方程,可以得到接近于實(shí)際加工過程的基體粗糙平面,在此基礎(chǔ)上再疊加非對(duì)稱織構(gòu)的形貌特征方程,使得非對(duì)稱織構(gòu)平面的表征更加具體,更加符合實(shí)際加工場景。

    目前,表面織構(gòu)截面形狀的研究主要集中在對(duì)稱織構(gòu)方向上,例如常見的圓形凹坑、矩形凹槽、圓形凸起等[10-13],針對(duì)非對(duì)稱表面織構(gòu)的研究較少,且尚未形成完整的理論體系。現(xiàn)有部分非對(duì)稱織構(gòu)相關(guān)研究初步表明,合理打破截面對(duì)稱性的表面織構(gòu)可增強(qiáng)空化效應(yīng)和慣性效應(yīng), 進(jìn)而提高結(jié)合面的潤滑能力,減小摩擦,還可實(shí)現(xiàn)提升承載能力的效果[14]。Ewoldt[15]等通過對(duì)具有縱向不對(duì)稱深度剖面的非對(duì)稱表面織構(gòu)牛頓流體間隙進(jìn)行了控制實(shí)驗(yàn),證明在沒有空化的情況下,必須打破表面織構(gòu)深度的對(duì)稱才能產(chǎn)生額外的法向承載力,不對(duì)稱的表面織構(gòu)降低了剪應(yīng)力,產(chǎn)生了正載荷,降低了有效摩擦因數(shù)。王婷等[16]通過研究正弦狀非對(duì)稱表面織構(gòu)的動(dòng)壓潤滑性能,證明在特定工況下,構(gòu)型合理的非對(duì)稱織構(gòu)比對(duì)稱織構(gòu)的動(dòng)壓潤滑性能更加優(yōu)異,這是由于正弦狀織構(gòu)結(jié)合面間隙的特殊形貌使得潤滑介質(zhì)產(chǎn)生了有利于潤滑的正向壓力差。劉冬妮[17]通過研究具有非對(duì)稱表面織構(gòu)結(jié)合面之間的壓力分布、氣相分?jǐn)?shù)、摩擦因數(shù)等參數(shù),分析了影響空化效應(yīng)的因素,結(jié)果表明,非對(duì)稱表面織構(gòu)的形貌參數(shù)通過改變潤滑介質(zhì)壓力和渦流的分布規(guī)律影響了承載力和摩擦因數(shù)。鄧志強(qiáng)[18]通過切削試驗(yàn)證明,非對(duì)稱表面織構(gòu)應(yīng)用在硬質(zhì)合金刀具切削鈦合金時(shí),前刀面的磨損有了很大的改善,減少了崩邊現(xiàn)象的發(fā)生率,刀具的磨損量顯著減少,同時(shí)切屑性能有所提高??号卫鸞19]的研究中發(fā)現(xiàn),將非對(duì)稱織構(gòu)應(yīng)用在制動(dòng)盤的優(yōu)化設(shè)計(jì)中,證明形狀合理的非對(duì)稱表面織構(gòu)使制動(dòng)盤第2棱邊對(duì)閘片的刮擦效應(yīng)顯著,具有更加高效的制動(dòng)效果。從以上研究可以發(fā)現(xiàn),非對(duì)稱織構(gòu)相關(guān)研究熱度在不斷提高,但是缺乏系統(tǒng)性的研究過程,且這些研究建立的幾何模型全部在理想的純平光滑平面基礎(chǔ)上,潤滑油膜實(shí)體模型的建立過程不夠貼合實(shí)際加工工況,缺少關(guān)于建立非對(duì)稱織構(gòu)平面表征方程的研究,因此開展系統(tǒng)性的非對(duì)稱織構(gòu)研究具有重要的實(shí)際意義。

    通過上述研究現(xiàn)狀的分析發(fā)現(xiàn),目前有關(guān)非對(duì)稱織構(gòu)平面的數(shù)學(xué)表征處于匱乏狀態(tài),只能通過在理想狀態(tài)下勾勒絕對(duì)光滑的平面,這就造成非對(duì)稱織構(gòu)平面描述與實(shí)際加工獲得的粗糙平面存在較大出入,進(jìn)而影響動(dòng)壓潤滑分析過程中油膜實(shí)體模型的準(zhǔn)確性。由此,文中在較為成熟的W-M函數(shù)表征粗糙表面的理論基礎(chǔ)上,創(chuàng)建具有非對(duì)稱表面織構(gòu)的粗糙表面表征方程,并基于方程表征的織構(gòu)平面建立潤滑油膜的幾何模型,對(duì)傳統(tǒng)分析過程進(jìn)行優(yōu)化。

    1 具有非對(duì)稱表面織構(gòu)的粗糙平面表征

    潤滑油膜幾何模型的建立是流體動(dòng)壓潤滑效應(yīng)分析的重要環(huán)節(jié),更加符合實(shí)際工況的織構(gòu)平面,可以直接提高分析過程的可靠性。文中在進(jìn)行非對(duì)稱織構(gòu)的動(dòng)壓潤滑效應(yīng)研究之前,先分析具有非對(duì)稱表面織構(gòu)粗糙平面的輪廓特征,通過建立表征方程的方式,創(chuàng)造一種重構(gòu)織構(gòu)平面的方法,并且分析了方程中各個(gè)參數(shù)對(duì)具體輪廓細(xì)節(jié)的影響機(jī)理,避免直接使用理想光滑平面,為后續(xù)潤滑油膜幾何模型的建立提供理論基礎(chǔ)。

    1.1 建立表征方程

    表面織構(gòu)的加工過程一般在粗糙表面的基礎(chǔ)上通過減材制造技術(shù)獲得,因此想要準(zhǔn)確地描述具有非對(duì)稱表面織構(gòu)的平面,必須先將加工前的粗糙表面輪廓特征表征出來,且經(jīng)過加工后織構(gòu)平面的無織構(gòu)區(qū)域形貌不發(fā)生改變,再疊加非對(duì)稱表面織構(gòu)的形狀特征,即可得出更接近實(shí)際工況、具有非對(duì)稱表面織構(gòu)的粗糙表面。

    在任何粗糙度要求下加工獲得的平面,在微觀狀態(tài)下都是粗糙不平的,其輪廓呈現(xiàn)出的形貌特點(diǎn)較為復(fù)雜。分形理論的相關(guān)研究證明[20-22],粗糙表面的隨機(jī)輪廓具有統(tǒng)計(jì)自相似性和自仿射性的特征,隨著表面觀測技術(shù)的提高,研究證明該特性在原子尺度仍然存在[23]?;诜中卫碚撎岢龅腤-M函數(shù)廣泛應(yīng)用于具有自相似性的輪廓表征中[24-25],工程中常用的三維函數(shù)模型見式(1)。

    式中:()為粗糙表面沿軸方向的隨機(jī)高度;G為尺度系數(shù);為輪廓空間頻率(>1,常取=1.5);s為分形維數(shù)(2AB為相互獨(dú)立且都服從[0,2π]上均勻分布的隨機(jī)數(shù)。

    文中所采用的非對(duì)稱表面織構(gòu)形式為截面形狀是不等腰三角形的溝槽,其結(jié)構(gòu)示意如圖1所示。其中,ll分別表示織構(gòu)單元的長邊和短邊沿相對(duì)運(yùn)動(dòng)方向的寬度;表示織構(gòu)寬度;h表示無織構(gòu)區(qū)域油膜厚度;表示織構(gòu)底部到無織構(gòu)平面的垂直距離,稱為織構(gòu)深度;表示兩相鄰織構(gòu)單元間的距離,稱為織構(gòu)間距。

    圖1 非對(duì)稱表面織構(gòu)示意

    根據(jù)非對(duì)稱織構(gòu)平面的加工過程,在原有W-M函數(shù)表征的粗糙平面基礎(chǔ)上,疊加非對(duì)稱織構(gòu)的形貌特征函數(shù),得到的表征函數(shù)如式(2)所示。

    為增強(qiáng)表征方程的普適性和進(jìn)一步簡化公式,特引入一個(gè)可以描述該形狀下各種尺寸非對(duì)稱織構(gòu)的參數(shù)——非對(duì)稱度同時(shí)將12轉(zhuǎn)化為織構(gòu)寬度:

    將式(3)、式(4)帶入式(2)中,得到最終的非對(duì)稱織構(gòu)平面的表征方程,如式(5)所示。

    1.2 織構(gòu)平面形貌仿真分析

    為了驗(yàn)證表征方程的準(zhǔn)確性,并且將方程中的具體參數(shù)對(duì)非對(duì)稱織構(gòu)平面表面特征的影響效果加以分析,文中將上述得到的表征方程在Matlab中進(jìn)行仿真模擬,利用surf命令編寫程序,繪制出非對(duì)稱織構(gòu)平面的三維圖及其對(duì)應(yīng)的側(cè)視圖。根據(jù)表征方程的建立過程,其參數(shù)可以分為2個(gè)部分:織構(gòu)形貌參數(shù)和粗糙表面參數(shù)。

    1.2.1 織構(gòu)形貌參數(shù)對(duì)織構(gòu)平面形貌的影響

    有關(guān)織構(gòu)形狀的參數(shù)包括非對(duì)稱度、織構(gòu)深度、織構(gòu)寬度、織構(gòu)間距等,選取非對(duì)稱度、織構(gòu)深度等2項(xiàng)對(duì)非對(duì)稱織構(gòu)形貌影響較為明顯的參數(shù)進(jìn)行討論,設(shè)置的織構(gòu)形貌參數(shù)對(duì)比方案如表1所示。

    3組非對(duì)稱度值=0.5、=1、=2分別對(duì)應(yīng)的織構(gòu)形狀為反向非對(duì)稱織構(gòu)、對(duì)稱織構(gòu)和正向非對(duì)稱織構(gòu),見圖2。在其他織構(gòu)形貌參數(shù)不變的情況下,非對(duì)稱度的改變可以實(shí)現(xiàn)非對(duì)稱形式的直接改變,其中當(dāng)=1時(shí),方程表征的是截面形狀為等腰三角形的對(duì)稱溝槽織構(gòu)。同時(shí),2組非對(duì)稱度數(shù)值互為倒數(shù)的非對(duì)稱織構(gòu)可以作為對(duì)照,二者等同于同一非對(duì)稱織構(gòu)進(jìn)行潤滑介質(zhì)沿相反方向流動(dòng),可用于對(duì)比潤滑介質(zhì)流動(dòng)方向?qū)Ψ菍?duì)稱織構(gòu)潤滑效果的影響。

    表1 織構(gòu)形貌參數(shù)對(duì)比方案

    Tab.1 Comparison scheme of texture morphology parameters

    從圖3可以發(fā)現(xiàn),織構(gòu)深度參數(shù)直觀影響無織構(gòu)區(qū)域到織構(gòu)底部的垂直距離,但是對(duì)于非對(duì)稱織構(gòu)在粗糙平面上的分布等特征沒有任何影響。圖3中的細(xì)節(jié)顯示,隨著織構(gòu)深度的增大,粗糙表面的輪廓高度起伏變得相對(duì)不明顯,從而使得織構(gòu)平面的整體輪廓的高度變化看起來更加趨于平緩。如果在相同的尺度維數(shù)下,這種無織構(gòu)區(qū)域的輪廓差異便會(huì)消失,這也說明織構(gòu)參數(shù)與無織構(gòu)區(qū)域形貌無關(guān)。

    圖2 非對(duì)稱度組三維圖像及側(cè)視圖

    圖3 織構(gòu)深度組三維圖像及側(cè)視圖

    從以上織構(gòu)參數(shù)對(duì)于織構(gòu)平面形貌的影響可以發(fā)現(xiàn),織構(gòu)形貌參數(shù)直接影響織構(gòu)的截面形狀,且各個(gè)參數(shù)之間相互獨(dú)立,不發(fā)生耦合作用。由此,可以利用所建立的表征方程進(jìn)行不同尺寸非對(duì)稱織構(gòu)平面的重構(gòu),實(shí)現(xiàn)相同加工條件下不同形狀特征的織構(gòu)平面表征。

    1.2.2 粗糙表面參數(shù)對(duì)織構(gòu)平面形貌的影響

    粗糙表面的相關(guān)參數(shù)主要包括尺度系數(shù)G和分形維數(shù)s,為了清晰地研究二者具體對(duì)織構(gòu)平面形貌的影響,文中設(shè)置的粗糙表面參數(shù)對(duì)比方案如表2所示。

    從圖4中可以發(fā)現(xiàn),尺度系數(shù)G主要影響表面粗糙表面輪廓高度的變化范圍,隨著尺度系數(shù)的減小,微觀輪廓高度的變化幅度也以同樣的倍率減小。同時(shí),隨著微觀輪廓高度變化幅度的減小,整個(gè)織構(gòu)平面也更加趨于平滑,織構(gòu)的截面形狀也更加規(guī)則,這一點(diǎn)在模擬圖像的側(cè)視圖上表現(xiàn)得十分清楚。

    表2 粗糙表面參數(shù)對(duì)比方案

    Tab.2 Comparison scheme of rough surface parameters

    從圖5中可以發(fā)現(xiàn),分形維數(shù)s主要影響粗糙表面輪廓高度的變化頻率,隨著分形維數(shù)逐漸增大,粗糙表面輪廓高度變化越來越劇烈,變化頻率越來越快,表面的復(fù)雜程度也越來越高。同時(shí),也對(duì)織構(gòu)平面的細(xì)節(jié)描述產(chǎn)生了較大影響,使得織構(gòu)截面形狀更加接近實(shí)際的加工效果,這也顯示了分形維數(shù)在粗糙表面表征參數(shù)中的重要地位。

    圖4 尺度系數(shù)組三維圖像及側(cè)視圖

    圖5 分型維數(shù)組三維圖像及側(cè)視圖

    2 非對(duì)稱表面織構(gòu)動(dòng)壓潤滑效應(yīng)分析

    從流體力學(xué)角度分析,每個(gè)單獨(dú)的表面織構(gòu)單元都相當(dāng)于1個(gè)動(dòng)壓軸承,結(jié)合面間潤滑介質(zhì)的流動(dòng)空間沿著“窄—寬—窄”的趨勢變化,內(nèi)部壓力的波動(dòng)形成了壓力差,引發(fā)了空化效應(yīng),從而產(chǎn)生了額外的承載能力,降低了摩擦[26-27]。由于空化效應(yīng)的發(fā)生需要滿足特定的條件,包括潤滑介質(zhì)的運(yùn)動(dòng)速度、溫度、壓力等[28-29],因此當(dāng)部分條件無法滿足時(shí)需要通過打破織構(gòu)單元沿運(yùn)動(dòng)方向的對(duì)稱性,以保證表面織構(gòu)的高效性。

    2.1 潤滑油膜幾何模型

    建立潤滑油膜幾何模型實(shí)質(zhì)上就是描述結(jié)合面間潤滑介質(zhì)實(shí)際厚度的過程,無織構(gòu)結(jié)合面潤滑油膜厚度始終保持恒定,而非對(duì)稱表面織構(gòu)結(jié)合面的潤滑油膜具有由發(fā)散到收縮的變化趨勢,非對(duì)稱表面織構(gòu)的引入是結(jié)合面實(shí)現(xiàn)動(dòng)壓潤滑的動(dòng)力源泉。無織構(gòu)結(jié)合面和非對(duì)稱織構(gòu)結(jié)合面潤滑油膜結(jié)構(gòu)示意如圖6所示,()表示結(jié)合面間的實(shí)際油膜厚度,無織構(gòu)結(jié)合面的()恒等于0。

    在動(dòng)壓潤滑條件下,為了準(zhǔn)確地建立符合工程實(shí)際的非對(duì)稱表面織構(gòu)結(jié)合面潤滑油膜幾何模型,采用了逆向建模的思路,選取單獨(dú)非對(duì)稱織構(gòu)單元,將前面表征方程通過Matlab軟件建立的三維曲面轉(zhuǎn)換成stl三角網(wǎng)格文件,導(dǎo)入SolidWorks軟件構(gòu)建織構(gòu)平面實(shí)體模型,整個(gè)動(dòng)壓潤滑油膜模型的建立過程如圖7所示。

    圖6 潤滑油膜結(jié)構(gòu)示意

    圖7 動(dòng)壓油膜模型的建立

    2.2 不可壓縮Reynolds方程

    潤滑問題的實(shí)質(zhì)在于研究潤滑介質(zhì)的流動(dòng)規(guī)律。文中織構(gòu)潤滑模型的流動(dòng)形式可視為不可壓縮流體的層流運(yùn)動(dòng),僅考慮潤滑油膜壓力形成的動(dòng)壓效應(yīng),用于流體動(dòng)壓潤滑計(jì)算的Reynolds方程[30-32],如式(6)所示。

    式中:為潤滑介質(zhì)粘度;為油膜厚度;為織構(gòu)內(nèi)部壓力;為潤滑介質(zhì)速度。

    根據(jù)構(gòu)建的表征方程,建立潤滑油膜模型的膜厚方程見式(7)。

    對(duì)方程中的物理量進(jìn)行無量綱化處理,采用的無量綱參數(shù)如式(8—11)所示。

    得到上述Reynolds方程的無量綱形式如式(12)所示。

    2.3 潤滑效應(yīng)CFD分析

    利用Ansys的fluent模塊進(jìn)行流體CFD仿真計(jì)算,進(jìn)而分析非對(duì)稱表面織構(gòu)相比于對(duì)稱織構(gòu)的優(yōu)越性,以及運(yùn)動(dòng)方向與非對(duì)稱度之間的關(guān)系。選取三維油膜實(shí)體模型沿軸方向的中間截面(= 1.5 mm)作為二維簡化分析平面,潤滑油的粘度為1.06 Pa·s,潤滑油密度為889 kg/m3,對(duì)正向非對(duì)稱織構(gòu)(=2)、對(duì)稱織構(gòu)(=1)、反向非對(duì)稱織構(gòu)(=0.5)分別進(jìn)行壓力分布分析和流速流線分析,得到的油膜壓力分布云圖和油膜速度流線圖如圖8—9所示。

    圖8 油膜壓力分布云圖

    圖9 油膜速度流線圖

    2.3.1 非對(duì)稱織構(gòu)的優(yōu)勢分析

    在正向非對(duì)稱織構(gòu)(=2)與對(duì)稱織構(gòu)(=1)的壓力云圖(圖8a—b)對(duì)比中可以發(fā)現(xiàn),對(duì)稱織構(gòu)的潤滑油膜在進(jìn)入織構(gòu)和離開織構(gòu)的兩側(cè)壓力差遠(yuǎn)小于非對(duì)稱織構(gòu)。在油膜壓力曲線圖10中可以看出,兩者在=0.2 mm的織構(gòu)入口處和=0.8 mm的織構(gòu)出口處的壓力差分別為1747 Pa和387 Pa,正向非對(duì)稱織構(gòu)明顯高于對(duì)稱織構(gòu)。由此可見,正向非對(duì)稱織構(gòu)可以顯著改善結(jié)合面的潤滑性能,產(chǎn)生更好的承載能力,抵消外部載荷對(duì)結(jié)合面的壓力,從而減小摩擦。

    圖10 k=2和k=1織構(gòu)油膜壓力曲線

    從正向非對(duì)稱織構(gòu)和對(duì)稱織構(gòu)的流線分布(圖9a—b)中可以看出,正向非對(duì)稱織構(gòu)的流線分布更加復(fù)雜、更加紊亂,并且回流更加顯著,這表明正向非對(duì)稱織構(gòu)在潤滑介質(zhì)的慣性效應(yīng)方面更具有優(yōu)勢,而慣性效應(yīng)可以有效提高結(jié)合面的承載能力,因此,正向非對(duì)稱織構(gòu)相比對(duì)稱織構(gòu)在結(jié)合面承載力提升方面存在明顯優(yōu)勢。

    為了驗(yàn)證計(jì)算分析結(jié)果的準(zhǔn)確性,將所得非對(duì)稱表面織構(gòu)在改善潤滑效應(yīng)方面的優(yōu)勢結(jié)論與文獻(xiàn)[17]進(jìn)行了對(duì)比。由于目前沒有按照非對(duì)稱度區(qū)分該形狀非對(duì)稱表面的相關(guān)研究,因此選用了文獻(xiàn)[17]中角度為50.19°和71.57°的織構(gòu)形狀與文中=1和=2進(jìn)行對(duì)比,文獻(xiàn)通過摩擦磨損試驗(yàn)得到的摩擦因數(shù)變化趨勢與文中研究進(jìn)出口壓力差變化趨勢的對(duì)比結(jié)果如圖11所示,文中研究織構(gòu)進(jìn)出口壓力差的增大與文獻(xiàn)[17]摩擦因數(shù)的減小均可證明非對(duì)稱織構(gòu)相比對(duì)稱織構(gòu)存在優(yōu)勢結(jié)論的可靠性。

    2.3.2 運(yùn)動(dòng)方向?qū)櫥匦缘挠绊?/p>

    =2的正向非對(duì)稱織構(gòu)和=0.5的反向非對(duì)稱織構(gòu)在本質(zhì)上可以視為潤滑介質(zhì)沿相反的方向在非對(duì)稱織構(gòu)內(nèi)流動(dòng),因此二者的對(duì)比可以用來研究結(jié)合面的相對(duì)運(yùn)動(dòng)方向?qū)Ψ菍?duì)稱織構(gòu)動(dòng)壓潤滑效應(yīng)的影響。

    圖11 k=2和k=1織構(gòu)與文獻(xiàn)[17]實(shí)驗(yàn)結(jié)果對(duì)比

    在正向和反向非對(duì)稱織構(gòu)的壓力對(duì)比(圖8a、c)中可以發(fā)現(xiàn),結(jié)合面的相對(duì)運(yùn)動(dòng)方向?qū)?dòng)壓潤滑效應(yīng)的影響效果十分明顯。進(jìn)出口壓力差的對(duì)比如圖12所示,正向非對(duì)稱織構(gòu)(1747 Pa)大于反向(707 Pa),證明當(dāng)織構(gòu)內(nèi)部潤滑介質(zhì)的流動(dòng)方向與=2的情形保持一致時(shí),可以更加顯著地提高非對(duì)稱織構(gòu)的承載能力。同樣,對(duì)比二者流線圖(圖9a、c)可知,正向非對(duì)稱織構(gòu)內(nèi)部的回流程度更加劇烈,慣性效應(yīng)的效果也更加優(yōu)異。由此可見,對(duì)非對(duì)稱織構(gòu)進(jìn)行更加合理的結(jié)構(gòu)設(shè)計(jì)可以在最大程度上提高其對(duì)結(jié)合面潤滑性能的正面作用。

    圖12 k=2和k=0.5織構(gòu)油膜壓力曲線

    同理,選用了文獻(xiàn)[17]中角度為40.6°和71.57°的織構(gòu)形狀與文中研究=0.5和=2的非對(duì)稱織構(gòu)進(jìn)行對(duì)比,對(duì)正向和反向非對(duì)稱織構(gòu)的進(jìn)出口壓力差變化趨勢加以驗(yàn)證,對(duì)比結(jié)果如圖13所示。同樣證明,

    圖13 k=2和k=0.5織構(gòu)與文獻(xiàn)[17]實(shí)驗(yàn)結(jié)果的對(duì)比

    結(jié)合面相對(duì)運(yùn)動(dòng)方向同文獻(xiàn)[17]中角度為40.6°或文中研究=2的非對(duì)稱織構(gòu)形狀一致時(shí),對(duì)提高結(jié)合面潤滑性能方面更加顯著。

    3 結(jié)論

    以表面織構(gòu)加工過程為思路,在W-M函數(shù)表征粗糙表面的基礎(chǔ)上,疊加非對(duì)稱織構(gòu)的形狀特征方程,最終得到具有非對(duì)稱表面織構(gòu)粗糙平面的表征方程,并且通過Matlab模擬重構(gòu)三維織構(gòu)平面對(duì)方程的準(zhǔn)確性加以驗(yàn)證。

    通過分別仿真分析織構(gòu)形狀參數(shù)和粗糙表面參數(shù)對(duì)整個(gè)織構(gòu)平面形貌特征的影響發(fā)現(xiàn),非對(duì)稱度、織構(gòu)深度等織構(gòu)相關(guān)參數(shù)主要影響織構(gòu)的形狀特征,尺度系數(shù)主要影響粗糙表面輪廓高度的變化范圍,分形維數(shù)主要影響粗糙表面輪廓高度的變化頻率。

    通過對(duì)比動(dòng)壓潤滑條件下正向非對(duì)稱織構(gòu)和對(duì)稱織構(gòu)的潤滑油膜壓力分布、流線分布發(fā)現(xiàn),正向非對(duì)稱織構(gòu)(=2)的進(jìn)出口壓力差值更大,因織構(gòu)存在而產(chǎn)生的額外承載力更大,從而更加有利于增大承載和改善摩擦;正向非對(duì)稱織構(gòu)內(nèi)部的流線分布也更加復(fù)雜,織構(gòu)區(qū)域內(nèi)回流更加明顯,潤滑介質(zhì)的慣性效應(yīng)更加顯著,對(duì)于改善結(jié)合面的潤滑效應(yīng)更具有優(yōu)越性。

    通過對(duì)比動(dòng)壓潤滑條件下正向非對(duì)稱織構(gòu)和反向非對(duì)稱織構(gòu)的壓力油膜壓力分布和流線分布發(fā)現(xiàn),=2的正向非對(duì)稱織構(gòu)具有更大的進(jìn)出口壓力差,即產(chǎn)生的額外承載力更加符合預(yù)期潤滑效果。當(dāng)結(jié)合面相對(duì)運(yùn)動(dòng)方向與正向非對(duì)稱織構(gòu)工況一致時(shí),潤滑效果最佳,因此非對(duì)稱織構(gòu)的應(yīng)用要結(jié)合合理的結(jié)構(gòu)設(shè)計(jì)和場景選擇。

    [1] PODSIADLO P, WOLSKI M, STACHOWIAK G W. Directional Signatures of Surface Texture[J]. Tribology Letters, 2019, 67(4): 1-13.

    [2] 祁鵬浩, 仝哲, 劉奇, 等. 表面織構(gòu)化DLC涂層在脂潤滑下的摩擦學(xué)性能研究[J]. 表面技術(shù), 2021, 50(1): 296- 304.

    QI Peng-hao, TONG Zhe, LIU Qi, et al. Tribological Properties of DLC Coated Textured Surfaces under GreaseLubrication[J]. Surface Technology, 2021, 50(1): 296-304.

    [3] SHARMA S, JAMWAL G, AWASTHI R K. Enhance-ment of Steady State Performance of Hydrodynamic Journal Bearing Using Chevron-Shaped Surface Texture[J]. Pro-ceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2019, 233(12): 1833-1843.

    [4] GROPPER D, WANG Ling, HARVEY T J. Hydrody-namic Lubrication of Textured Surfaces: A Review of Modeling Techniques and Key Findings[J]. Tribology International, 2016, 94: 509-529.

    [5] 郝秀清, 宋曉路, 李亮. 表面織構(gòu)化刀具的研究現(xiàn)狀與進(jìn)展[J]. 表面技術(shù), 2016, 45(9): 170-181.

    HAO Xiu-qing, SONG Xiao-lu, LI Liang. Development and Perspective of Surface Texturing Tools[J]. Surface Technology, 2016, 45(9): 170-181.

    [6] TOMANIK E, PROFITO F J, ZACHARIADIS D C. Modelling the Hydrodynamic Support of Cylinder Bore and Piston Rings with Laser Textured Surfaces[J]. Tribo-logy International, 2013, 59: 90-96.

    [7] 王靜秋, 王曉雷. 表面織構(gòu)創(chuàng)新設(shè)計(jì)的研究回顧及展望[J]. 機(jī)械工程學(xué)報(bào), 2015, 51(23): 84-95.

    WANG Jing-qiu, WANG Xiao-lei. State of the Art in Innovative Design of Surface Texture[J]. Journal of Mechanical Engineering, 2015, 51(23): 84-95.

    [8] MAJUMDAR A, BHUSHAN B. Fractal Model of Elastic- Plastic Contact between Rough Surfaces[J]. Journal of Tribology, 1991, 113(1): 1-11.

    [9] YAN W, KOMVOPOULOS K. Contact Analysis of Elastic-Plastic Fractal Surfaces[J]. Journal of Applied Physics, 1998, 84(7): 3617-3624.

    [10] YUE Hong-zhi, DENG Jian-xin, GE Dong-liang, et al. Effect of Surface Texturing on Tribological Performance of Sliding Guideway under Boundary Lubrication[J]. Journal of Manufacturing Processes, 2019, 47: 172-182.

    [11] 王麗麗, 郭少輝, 魏聿梁, 等. 表面微織構(gòu)對(duì)45#鋼摩擦副表面摩擦學(xué)性能影響的實(shí)驗(yàn)研究[J]. 表面技術(shù), 2018, 47(12): 149-154.

    WANG Li-li, GUO Shao-hui, WEI Yu-liang, et al. Experimental Study on the Influence of Surface Micro-texture on the Tribological Properties of 45#Steel Friction Pairs[J]. Surface Technology, 2018, 47(12): 149-154.

    [12] MO Fan, SHEN Cong, ZHOU Jia, et al. Statistical Analysis of the Influence of Imperfect Texture Shape and Dimensional Uncertainty on Surface Texture Performance [J]. IEEE Access, 2017, 5: 27023-27035.

    [13] 屠春娟, 郭旭紅, 郭大林, 等. 不同形貌微織構(gòu)自潤滑陶瓷刀具切削性能的對(duì)比[J]. 機(jī)械工程材料, 2018, 42(11): 47-51.

    TU Chun-juan, GUO Xu-hong, GUO Da-lin, et al. Com-parison of Machinability of Self-Lubrication Ceramic Tools with Different Morphological Micro-Textures[J]. Materials for Mechanical Engineering, 2018, 42(11): 47-51.

    [14] 張瑜, 陳國定, 王琳, 等. 空化與慣性效應(yīng)耦合作用下的非對(duì)稱表面微織構(gòu)滑塊承載力分析[J]. 西北工業(yè)大學(xué)學(xué)報(bào), 2017, 35(6): 1026-1032.

    ZHANG Yu, CHEN Guo-ding, WANG Lin, et al. Analy-sis of Bearing Load-Carrying Capacity with Asymmetric Surface Textures under Coupling Effects of Cavitation and Inertia Effect[J]. Journal of Northwestern Polytech-nical University, 2017, 35(6): 1026-1032.

    [15] SCHUH J K, EWOLDT R H. Asymmetric Surface Textures Decrease Friction with Newtonian Fluids in Full Film Lubricated Sliding Contact[J]. Tribology Interna-tional, 2016, 97: 490-498.

    [16] 王婷, 傅連東, 湛從昌, 等. 正弦狀非對(duì)稱織構(gòu)的動(dòng)壓潤滑性能[J]. 機(jī)械設(shè)計(jì)與制造, 2020(3): 103-106.

    WANG Ting, FU Lian-dong, ZHAN Cong-chang, et al. Dynamic Lubrication Performance of the Sinusoidal Asy-mmetric Texture[J]. Machinery Design & Manufacture, 2020(3): 103-106.

    [17] 劉東妮. 非對(duì)稱微織構(gòu)對(duì)摩擦副潤滑特性的影響[D]. 大連: 大連海事大學(xué), 2019: 36-51.

    LIU Dong-ni. The Effect of Asymmetric Micro-Texture on Lubrication Characteristics of Friction Pairs[D]. Dalian: Dalian Maritime University, 2019: 36-51.

    [18] 鄧志強(qiáng). 非對(duì)稱形貌織構(gòu)刀具干切削鈦合金實(shí)驗(yàn)研究[D]. 湘潭: 湘潭大學(xué), 2018: 36-60.

    DENG Zhi-qiang. Experimental Study on Dry Cutting of Titanium Alloy with Asymmetric Texture Tool[D]. Xiangtan: Xiangtan University, 2018: 36-60.

    [19] 亢盼利. 制動(dòng)盤表面溝槽對(duì)制動(dòng)過程的動(dòng)力學(xué)影響研究及溝槽優(yōu)化[D]. 大連: 大連交通大學(xué), 2020: 7-33.

    KANG Pan-li. Research on the Dynamics Effect of Brake Disc Surface Groove on Braking and Groove Optimiza-tion[D]. Dalian: Dalian Jiaotong University, 2020: 7-33.

    [20] 陳國安, 葛世榮, 王軍祥. 分形理論在摩擦學(xué)研究中的應(yīng)用[J]. 摩擦學(xué)學(xué)報(bào), 1998, 18(2): 84-89.

    CHEN Guo-an, GE Shi-rong, WANG Jun-xiang. Applica-tion of Fractal Geometry in Tribology[J]. Tribology, 1998, 18(2): 84-89.

    [21] GANTI S, BHUSHAN B. Generalized Fractal Analysis and Its Applications to Engineering Surfaces[J]. Wear, 1995, 180(1/2): 17-34.

    [22] MAJUMDAR A, BHUSHAN B. Role of Fractal Geo-metry in Roughness Characterization and Contact Me-chanics of Surfaces[J]. Journal of Tribology, 1990, 112(2): 205-216.

    [23] SAYLES R S, THOMAS T R. Surface Topography as a Nonstationary Random Process[J]. Nature, 1978, 271(5644): 431-434.

    [24] 鄧可月, 劉政, 鄧居軍, 等. W-M函數(shù)模型下表面輪廓形貌的變化規(guī)律[J]. 機(jī)械設(shè)計(jì)與制造, 2017(1): 47-50.

    DENG Ke-yue, LIU Zheng, DENG Ju-jun, et al. Varia-tion of Surface Profile Topography Based on W-M Func-tion Model[J]. Machinery Design & Manufacture, 2017(1): 47-50.

    [25] 劉宇, 鄧宏盛, 張生芳, 等. 基于W-M分形函數(shù)的三維粗糙表面摩擦生熱研究[J]. 中國工程機(jī)械學(xué)報(bào), 2018, 16(3): 194-201.

    LIU Yu, DENG Hong-sheng, ZHANG Sheng-fang, et al. Research on Friction Heat Generation of Three Dimen-sional Rough Surface Based on W-M Fractal Function[J]. Chinese Journal of Construction Machinery, 2018, 16(3): 194-201.

    [26] 華??? 王蓉, 周萬, 等. 45鋼的黏結(jié)型激光微織構(gòu)表面摩擦學(xué)性能及固體潤滑機(jī)理分析[J]. 材料工程, 2015, 43(9): 39-45.

    HUA Xi-jun, WANG Rong, ZHOU Wan, et al. Analysis of Tribological Property and Mechanism of Solid Lubri-cation on 45 Steel Bond-Type Laser Micro-Textured Sur-face[J]. Journal of Materials Engineering, 2015, 43(9): 39-45.

    [27] WOS S, KOSZELA W, PAWLUS P, et al. Effects of Surface Texturing and Kind of Lubricant on the Coefficient of Friction at Ambient and Elevated Temperatures[J]. Tribology International, 2018, 117: 174-179.

    [28] 王麗麗, 郭少輝, 尹國校, 等. 計(jì)入空化效應(yīng)的表面微織構(gòu)滑動(dòng)軸承特性分析[J]. 潤滑與密封, 2018, 43(5): 65-69.

    WANG Li-li, GUO Shao-hui, YIN Guo-xiao, et al. Re-search on Micro-Textured Journal Bearing Characteristics Considering the Effect of Cavitation[J]. Lubrication Engineering, 2018, 43(5): 65-69.

    [29] HIRAYAMA T, SAKURAI T, YABE H. A Theoretical Analysis Considering Cavitation Occurrence in Oil- Lubricated Spiral-Grooved Journal Bearings with Experi-mental Verification[J]. Journal of Tribology, 2004, 126(3): 490-498.

    [30] ROM M, MüLLER S. A New Model for Textured Surface Lubrication Based on a Modified Reynolds Equation Including Inertia Effects[J]. Tribology International, 2019, 133: 55-66.

    [31] 韓靜, 方亮, 孫甲鵬, 等. Reynolds方程在紋理表面動(dòng)壓潤滑計(jì)算中的有效性評(píng)價(jià)[J]. 摩擦學(xué)學(xué)報(bào), 2014, 34(4): 348-356.

    HAN Jing, FANG Liang, SUN Jia-peng, et al. On the Validation of Reynolds Equation for Hydrodynamics Lubrication Simulation of Textured Surface[J]. Tribology, 2014, 34(4): 348-356.

    [32] ZHANG Hui, HUA Meng, DONG Guang-neng, et al. A Mixed Lubrication Model for Studying Tribological Behaviors of Surface Texturing[J]. Tribology International, 2016, 93: 583-592.

    Characterization of Asymmetric Surface Texture and Analysis of Hydrodynamic Lubrication Effect

    ,,,

    (School of Mechanical Engineering, Shenyang Jianzhu University, Shenyang 110168, China)

    The work aims to establish the characterization equation of the rough plane with asymmetric surface texture and analyze the effect of the equation parameters on the morphology and study the lubrication characteristics of the asymmetric surface texture under hydrodynamic lubrication conditions based on the texture plane reconstructed by the characterization equation. According to the machining process of the asymmetric texture plane, the rough surface was reconstructed based on the W-M function and then the morphology characteristic equation of the asymmetric texture was superimposed to achieve the mathematical characterization of the asymmetric surface texture of different sizes. The effects of surface texture parameters and rough surface parameters in the equation on the texture plane morphology were analyzed by Matlab software. The solid model of the hydrodynamic lubrication pressure oil film was established by the texture plane represented by the equation, and the CFD calculation of the pressure oil film was carried out. The pressure distribution law and streamline distribution law in the forward asymmetric, symmetric and reverse asymmetric textures were obtained respectively, and the hydrodynamic lubrication effect analysis of the asymmetric surface texture was realized. Among the parameters of the equation, the scale coefficient and fractal dimension mainly affected the variation range and frequency of the rough surface contour height. The microscopic contour height decreased at the same rate with the decrease of the scale coefficient, and the variation of the microscopic contour height became more and more violent with the gradual increase of the fractal dimension, and the frequency of variation became faster and faster. Some texture-related parameters, such as asymmetry rate and texture depth, specifically affected the morphology of the texture. The CFD analysis and calculation of lubrication pressure film showed that the pressure difference between the inlet and outlet of the forward asymmetric texture (=2) was 1747 Pa, which was significantly greater than the 387 Pa of the symmetric texture (=1) and the 707 Pa of the reverse asymmetric texture (=0.5). The streamline distribution nephogram showed that the reflux effect in the forward asymmetric texture was more intense and the streamline distribution was more complex. The established characterization equation of the asymmetric texture plane can clearly describe its morphology, and the effect mechanism of the equation parameters on the morphology is clear. Moreover, the asymmetrical texture with reasonable configuration has obvious advantages in increasing load capacity and improving lubrication.

    asymmetric; surface texture; characterization equation; W-M function; CFD; hydrodynamic lubrication effect

    2021-03-31;

    2021-08-10

    ZHU Chun-xia (1980—), Female, Doctor, Professor, Research focus: surface texture and related technology of mechanical joint surface.

    朱春霞, 王德全, 楊曉楠, 等. 非對(duì)稱表面織構(gòu)表征及動(dòng)壓潤滑效應(yīng)分析[J]. 表面技術(shù), 2022, 51(3): 66-75.

    TH117.2

    A

    1001-3660(2022)03-0066-10

    10.16490/j.cnki.issn.1001-3660.2022.03.006

    2021-03-31;

    2021-08-10

    國家自然科學(xué)基金面上項(xiàng)目(51575365);遼寧省“興遼英才計(jì)劃”項(xiàng)目(XLYC1807065)

    Fund:Suported by the National Natural Science Foundation of China (51575365); The Liaoning Provincial Natural "Xingliao Talent Plan" Project (XLYC1807065)

    朱春霞(1980—),女,博士,教授,主要研究方向?yàn)闄C(jī)械結(jié)合面表面織構(gòu)及相關(guān)技術(shù)。

    ZHU Chun-xia, WANG De-quan, YANG Xiao-nan, et al. Characterization of Asymmetric Surface Texture and Analysis of Hydrodynamic Lubrication Effect[J]. Surface Technology, 2022, 51(3): 66-75.

    猜你喜歡
    動(dòng)壓織構(gòu)非對(duì)稱
    冷軋壓下率對(duì)3104鋁合金織構(gòu)演變的影響
    國內(nèi)首個(gè)現(xiàn)代箔片氣體動(dòng)壓軸承技術(shù)培訓(xùn)班在長沙成功舉辦
    非對(duì)稱Orlicz差體
    研磨拋光表面微孔織構(gòu)的形成
    點(diǎn)數(shù)不超過20的旗傳遞非對(duì)稱2-設(shè)計(jì)
    織構(gòu)布置位置對(duì)滑動(dòng)軸承承載力的影響分析
    船海工程(2015年4期)2016-01-05 15:53:32
    南屯煤礦深部泵房硐室群動(dòng)壓失穩(wěn)機(jī)理及控制對(duì)策
    強(qiáng)烈動(dòng)壓巷道支護(hù)技術(shù)探討
    非對(duì)稱負(fù)載下矩陣變換器改進(jìn)型PI重復(fù)控制
    電測與儀表(2015年4期)2015-04-12 00:43:04
    Time constant of a hydraulic servo valve withdynamic pressure feedback
    亚洲欧洲精品一区二区精品久久久| 欧美日韩中文字幕国产精品一区二区三区| 亚洲中文字幕日韩| 9191精品国产免费久久| 一进一出抽搐动态| 久久精品夜夜夜夜夜久久蜜豆 | 欧美av亚洲av综合av国产av| 91九色精品人成在线观看| 精品熟女少妇八av免费久了| 久久久久久大精品| 2021天堂中文幕一二区在线观 | 91成人精品电影| 欧美黑人精品巨大| 国产精品乱码一区二三区的特点| 最近在线观看免费完整版| 国产欧美日韩精品亚洲av| 久久精品夜夜夜夜夜久久蜜豆 | 狂野欧美激情性xxxx| 最近最新中文字幕大全免费视频| 极品教师在线免费播放| 免费看a级黄色片| 少妇熟女aⅴ在线视频| √禁漫天堂资源中文www| 免费人成视频x8x8入口观看| 免费av毛片视频| 国产精品综合久久久久久久免费| 国产激情欧美一区二区| 很黄的视频免费| 99国产精品一区二区蜜桃av| 久久天堂一区二区三区四区| 每晚都被弄得嗷嗷叫到高潮| 看黄色毛片网站| 熟女少妇亚洲综合色aaa.| 精品午夜福利视频在线观看一区| 久久中文字幕一级| 草草在线视频免费看| 亚洲精品在线美女| 黄色a级毛片大全视频| 国产aⅴ精品一区二区三区波| 亚洲五月天丁香| 99在线人妻在线中文字幕| 波多野结衣av一区二区av| 大型黄色视频在线免费观看| 少妇熟女aⅴ在线视频| 久久久久久久久免费视频了| 精品熟女少妇八av免费久了| 日韩欧美免费精品| 他把我摸到了高潮在线观看| 国产爱豆传媒在线观看 | 可以免费在线观看a视频的电影网站| 日本a在线网址| 亚洲,欧美精品.| 免费在线观看影片大全网站| 99久久99久久久精品蜜桃| 亚洲avbb在线观看| 久久九九热精品免费| 丝袜美腿诱惑在线| 美女免费视频网站| 亚洲熟女毛片儿| 久久久久久免费高清国产稀缺| 亚洲在线自拍视频| 狂野欧美激情性xxxx| 国产亚洲av嫩草精品影院| 日韩欧美免费精品| 深夜精品福利| 日韩精品青青久久久久久| 国产乱人伦免费视频| 精品国内亚洲2022精品成人| 老司机午夜福利在线观看视频| 99久久无色码亚洲精品果冻| 变态另类丝袜制服| 日本 av在线| 午夜福利高清视频| 亚洲色图av天堂| 黄色毛片三级朝国网站| 亚洲人成电影免费在线| 国产亚洲精品久久久久5区| 夜夜爽天天搞| 哪里可以看免费的av片| 国产真人三级小视频在线观看| 婷婷精品国产亚洲av| 一级毛片高清免费大全| 午夜成年电影在线免费观看| 国产国语露脸激情在线看| 搡老妇女老女人老熟妇| 亚洲成av片中文字幕在线观看| 免费看a级黄色片| 亚洲精品国产区一区二| 天天一区二区日本电影三级| 看黄色毛片网站| 中文字幕久久专区| 国产激情久久老熟女| 久久香蕉激情| 亚洲国产精品成人综合色| 香蕉av资源在线| av天堂在线播放| 一二三四社区在线视频社区8| 国产免费av片在线观看野外av| 亚洲一区二区三区色噜噜| 免费在线观看视频国产中文字幕亚洲| 午夜久久久在线观看| 1024视频免费在线观看| 啦啦啦 在线观看视频| 黄色女人牲交| 亚洲精品久久成人aⅴ小说| 亚洲国产精品999在线| 少妇粗大呻吟视频| 欧美成人午夜精品| 国产成+人综合+亚洲专区| 亚洲午夜精品一区,二区,三区| 精品国产乱子伦一区二区三区| 国产单亲对白刺激| 国产成人一区二区三区免费视频网站| 欧美在线一区亚洲| 90打野战视频偷拍视频| 两个人免费观看高清视频| 好男人电影高清在线观看| 精品国产亚洲在线| 中文亚洲av片在线观看爽| www日本在线高清视频| 最新美女视频免费是黄的| 色综合欧美亚洲国产小说| 国产亚洲欧美98| bbb黄色大片| 亚洲一卡2卡3卡4卡5卡精品中文| 成人午夜高清在线视频 | 亚洲精品国产一区二区精华液| 女警被强在线播放| 久久久久九九精品影院| 国产成人av教育| 99在线视频只有这里精品首页| 夜夜爽天天搞| 亚洲一码二码三码区别大吗| 欧美性猛交黑人性爽| 亚洲熟妇熟女久久| 欧美黑人欧美精品刺激| 国产精品久久久久久精品电影 | 99精品久久久久人妻精品| 一区二区三区精品91| 午夜精品久久久久久毛片777| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕人妻丝袜一区二区| 草草在线视频免费看| 国产真人三级小视频在线观看| 久久精品人妻少妇| 身体一侧抽搐| 亚洲人成网站在线播放欧美日韩| 中文字幕精品免费在线观看视频| 两个人看的免费小视频| 校园春色视频在线观看| 可以在线观看毛片的网站| 91老司机精品| 欧美久久黑人一区二区| 中文在线观看免费www的网站 | 国产日本99.免费观看| 无人区码免费观看不卡| 亚洲色图 男人天堂 中文字幕| 亚洲五月天丁香| 亚洲精品久久国产高清桃花| 色精品久久人妻99蜜桃| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品粉嫩美女一区| 久久久久久久久久黄片| 99国产精品一区二区蜜桃av| 久久中文字幕一级| 久久中文字幕一级| 亚洲九九香蕉| 国产亚洲精品久久久久5区| 欧美激情高清一区二区三区| 听说在线观看完整版免费高清| 久久香蕉国产精品| АⅤ资源中文在线天堂| 亚洲精品粉嫩美女一区| 国产一区二区在线av高清观看| 美女高潮喷水抽搐中文字幕| 欧美性猛交╳xxx乱大交人| 午夜福利免费观看在线| 久久久久久国产a免费观看| 波多野结衣高清作品| 亚洲av成人av| 长腿黑丝高跟| 亚洲精品美女久久久久99蜜臀| 亚洲激情在线av| 久久天堂一区二区三区四区| 国产一卡二卡三卡精品| 亚洲激情在线av| 深夜精品福利| 一进一出好大好爽视频| 亚洲第一青青草原| 欧美中文日本在线观看视频| 一进一出抽搐动态| 亚洲欧美精品综合久久99| 淫妇啪啪啪对白视频| 亚洲在线自拍视频| 国产激情欧美一区二区| 亚洲中文日韩欧美视频| 久久久久免费精品人妻一区二区 | 亚洲人成网站在线播放欧美日韩| 777久久人妻少妇嫩草av网站| 欧美性长视频在线观看| 婷婷亚洲欧美| 国产亚洲av高清不卡| 欧美日韩黄片免| 欧美亚洲日本最大视频资源| 动漫黄色视频在线观看| 国产一区二区三区视频了| 国产爱豆传媒在线观看 | av免费在线观看网站| 国产亚洲精品久久久久久毛片| 在线观看免费视频日本深夜| 久久香蕉精品热| 成人av一区二区三区在线看| 男人舔奶头视频| 国产精品久久久av美女十八| 国产在线观看jvid| 黄色女人牲交| 国产精品一区二区三区四区久久 | 日本 欧美在线| 久久午夜亚洲精品久久| 中文字幕精品亚洲无线码一区 | 欧美一级a爱片免费观看看 | www国产在线视频色| 色av中文字幕| 久久人妻av系列| 国产精品电影一区二区三区| 亚洲午夜理论影院| 亚洲va日本ⅴa欧美va伊人久久| 成人永久免费在线观看视频| 曰老女人黄片| 麻豆成人av在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久欧美精品欧美久久欧美| 美女免费视频网站| 国产蜜桃级精品一区二区三区| 午夜a级毛片| 免费高清视频大片| a级毛片a级免费在线| 一夜夜www| 国产精品永久免费网站| 亚洲天堂国产精品一区在线| 99精品在免费线老司机午夜| 麻豆一二三区av精品| 欧美乱色亚洲激情| 国语自产精品视频在线第100页| 激情在线观看视频在线高清| 亚洲 欧美 日韩 在线 免费| 久久中文看片网| 女性生殖器流出的白浆| 最近最新中文字幕大全免费视频| a级毛片在线看网站| 18禁裸乳无遮挡免费网站照片 | 国产成人啪精品午夜网站| 国内精品久久久久精免费| netflix在线观看网站| 日韩精品中文字幕看吧| 亚洲精品中文字幕在线视频| 亚洲全国av大片| 日韩 欧美 亚洲 中文字幕| 日本精品一区二区三区蜜桃| 露出奶头的视频| 亚洲国产欧美一区二区综合| 别揉我奶头~嗯~啊~动态视频| 黄色片一级片一级黄色片| 久久草成人影院| 国产aⅴ精品一区二区三区波| 亚洲国产欧美网| www日本在线高清视频| 国产亚洲av高清不卡| 此物有八面人人有两片| avwww免费| 国产亚洲欧美在线一区二区| 国产高清有码在线观看视频 | 午夜福利在线在线| 国产一区二区在线av高清观看| 极品教师在线免费播放| a在线观看视频网站| av片东京热男人的天堂| 欧美日韩亚洲国产一区二区在线观看| 久久国产精品人妻蜜桃| 国产成+人综合+亚洲专区| 一区二区三区激情视频| 亚洲 国产 在线| 日本 av在线| or卡值多少钱| 国产成人影院久久av| 欧美黄色淫秽网站| 欧美久久黑人一区二区| 亚洲一区二区三区不卡视频| 亚洲第一电影网av| 免费一级毛片在线播放高清视频| 国产免费男女视频| 99久久久亚洲精品蜜臀av| 亚洲第一av免费看| 69av精品久久久久久| 精品无人区乱码1区二区| 国产亚洲精品综合一区在线观看 | 欧美日本亚洲视频在线播放| 日韩av在线大香蕉| 九色国产91popny在线| av中文乱码字幕在线| 男女视频在线观看网站免费 | 日韩一卡2卡3卡4卡2021年| 亚洲午夜精品一区,二区,三区| 婷婷六月久久综合丁香| 国产一区二区三区在线臀色熟女| 国产精品久久久久久人妻精品电影| 国产欧美日韩一区二区精品| 精品久久久久久成人av| 亚洲午夜理论影院| 久久久久久免费高清国产稀缺| 亚洲真实伦在线观看| 十八禁网站免费在线| 夜夜爽天天搞| 亚洲一区高清亚洲精品| 日日夜夜操网爽| 熟女电影av网| 美女午夜性视频免费| 一二三四社区在线视频社区8| 免费女性裸体啪啪无遮挡网站| 欧美乱妇无乱码| 国产精品 欧美亚洲| 精品卡一卡二卡四卡免费| 又紧又爽又黄一区二区| 亚洲第一av免费看| 国内毛片毛片毛片毛片毛片| 人人妻人人澡人人看| 午夜成年电影在线免费观看| 国产伦一二天堂av在线观看| 久久久久久人人人人人| 亚洲国产欧洲综合997久久, | 韩国av一区二区三区四区| 身体一侧抽搐| 国产私拍福利视频在线观看| 亚洲av熟女| 可以在线观看的亚洲视频| 黄片播放在线免费| 亚洲一区中文字幕在线| 嫩草影视91久久| 国产欧美日韩精品亚洲av| 这个男人来自地球电影免费观看| 国内精品久久久久久久电影| 久久青草综合色| 精品国产超薄肉色丝袜足j| 国产三级在线视频| 哪里可以看免费的av片| 午夜免费激情av| √禁漫天堂资源中文www| 午夜日韩欧美国产| 精品国产超薄肉色丝袜足j| 久久狼人影院| 女生性感内裤真人,穿戴方法视频| 别揉我奶头~嗯~啊~动态视频| 天堂影院成人在线观看| 成人国产综合亚洲| 国产精品久久久av美女十八| 国产亚洲欧美在线一区二区| 国产av又大| 男女那种视频在线观看| 777久久人妻少妇嫩草av网站| 午夜视频精品福利| 1024香蕉在线观看| 日本免费a在线| 精品无人区乱码1区二区| 久久 成人 亚洲| ponron亚洲| 亚洲久久久国产精品| 久久久久九九精品影院| 国产精品 欧美亚洲| 香蕉久久夜色| 国产野战对白在线观看| 亚洲国产欧洲综合997久久, | 日本成人三级电影网站| 亚洲熟妇中文字幕五十中出| 麻豆av在线久日| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 亚洲精品在线观看二区| 国产精品爽爽va在线观看网站 | 日韩欧美在线二视频| av免费在线观看网站| 1024香蕉在线观看| 老鸭窝网址在线观看| 亚洲欧美精品综合久久99| 欧美性长视频在线观看| 美女大奶头视频| 日日干狠狠操夜夜爽| 老司机福利观看| 激情在线观看视频在线高清| 亚洲激情在线av| 欧美色视频一区免费| 一级毛片女人18水好多| 12—13女人毛片做爰片一| 国产黄a三级三级三级人| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩欧美在线二视频| 一级a爱视频在线免费观看| 亚洲国产毛片av蜜桃av| 亚洲av第一区精品v没综合| 精品欧美一区二区三区在线| 亚洲最大成人中文| 男女下面进入的视频免费午夜 | 热99re8久久精品国产| 亚洲国产欧美网| 国产爱豆传媒在线观看 | 美女大奶头视频| 亚洲色图 男人天堂 中文字幕| 精品一区二区三区av网在线观看| 亚洲中文字幕一区二区三区有码在线看 | 久久精品夜夜夜夜夜久久蜜豆 | 国产精品久久久av美女十八| 久久中文字幕一级| 亚洲欧美精品综合久久99| 精品福利观看| 久久中文看片网| 99在线视频只有这里精品首页| 我的亚洲天堂| 老汉色∧v一级毛片| 亚洲精品美女久久久久99蜜臀| 91国产中文字幕| 黑人欧美特级aaaaaa片| 一区二区三区激情视频| 在线免费观看的www视频| 麻豆成人午夜福利视频| 一区二区三区国产精品乱码| 欧美日本亚洲视频在线播放| 狠狠狠狠99中文字幕| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 丝袜美腿诱惑在线| 热99re8久久精品国产| 中文资源天堂在线| 国产精品久久久av美女十八| 99久久无色码亚洲精品果冻| 午夜久久久久精精品| 精品一区二区三区视频在线观看免费| 亚洲 欧美一区二区三区| 国产一级毛片七仙女欲春2 | 91麻豆av在线| 麻豆成人av在线观看| 国产av不卡久久| 国产日本99.免费观看| 午夜两性在线视频| 亚洲熟妇熟女久久| 国产视频内射| 777久久人妻少妇嫩草av网站| 久久精品aⅴ一区二区三区四区| 1024视频免费在线观看| 欧美久久黑人一区二区| www.自偷自拍.com| 99国产综合亚洲精品| 熟女少妇亚洲综合色aaa.| 久久香蕉精品热| 亚洲精品国产精品久久久不卡| 热99re8久久精品国产| 久久久久久免费高清国产稀缺| 欧美黄色片欧美黄色片| 最近在线观看免费完整版| 亚洲午夜理论影院| 99re在线观看精品视频| 禁无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 村上凉子中文字幕在线| 一边摸一边抽搐一进一小说| 亚洲 欧美一区二区三区| 国产精品影院久久| 老鸭窝网址在线观看| 黑人欧美特级aaaaaa片| 亚洲一区二区三区色噜噜| 久久中文字幕一级| 久久午夜综合久久蜜桃| 每晚都被弄得嗷嗷叫到高潮| 性色av乱码一区二区三区2| 久久人妻av系列| 人人妻人人澡欧美一区二区| 亚洲精品久久国产高清桃花| 亚洲av成人一区二区三| 国产av一区二区精品久久| 后天国语完整版免费观看| 久久精品亚洲精品国产色婷小说| 男女下面进入的视频免费午夜 | 国产精品一区二区免费欧美| 18禁观看日本| 亚洲精品中文字幕一二三四区| 国产片内射在线| 女同久久另类99精品国产91| 露出奶头的视频| 日韩欧美国产在线观看| 亚洲真实伦在线观看| 午夜免费激情av| 亚洲国产欧美网| 午夜福利欧美成人| 欧美日韩福利视频一区二区| 婷婷亚洲欧美| 久久久久久久午夜电影| 久久久久国产一级毛片高清牌| 久久国产精品影院| 亚洲性夜色夜夜综合| 九色国产91popny在线| 精品久久久久久,| 黄片大片在线免费观看| 又黄又粗又硬又大视频| 两个人免费观看高清视频| 香蕉国产在线看| 亚洲一区二区三区不卡视频| 亚洲第一青青草原| 99久久国产精品久久久| 女同久久另类99精品国产91| 91av网站免费观看| 精品国内亚洲2022精品成人| 久久精品夜夜夜夜夜久久蜜豆 | 精品乱码久久久久久99久播| 久久精品国产99精品国产亚洲性色| 一边摸一边做爽爽视频免费| 一个人免费在线观看的高清视频| 亚洲欧美日韩高清在线视频| 19禁男女啪啪无遮挡网站| 色播在线永久视频| 欧美性猛交╳xxx乱大交人| 亚洲av中文字字幕乱码综合 | xxxwww97欧美| 午夜福利在线在线| 这个男人来自地球电影免费观看| 波多野结衣巨乳人妻| 一进一出好大好爽视频| www国产在线视频色| 午夜福利高清视频| 这个男人来自地球电影免费观看| 色综合婷婷激情| 天堂动漫精品| svipshipincom国产片| 亚洲精品一区av在线观看| 午夜福利视频1000在线观看| 亚洲最大成人中文| or卡值多少钱| 啦啦啦韩国在线观看视频| 美女午夜性视频免费| 99riav亚洲国产免费| 超碰成人久久| 国产黄a三级三级三级人| 正在播放国产对白刺激| 啦啦啦免费观看视频1| 国产麻豆成人av免费视频| 亚洲 欧美 日韩 在线 免费| 色播亚洲综合网| 日韩国内少妇激情av| 亚洲狠狠婷婷综合久久图片| 久久青草综合色| 国产私拍福利视频在线观看| 午夜福利成人在线免费观看| 女人高潮潮喷娇喘18禁视频| 大香蕉久久成人网| 757午夜福利合集在线观看| 天天添夜夜摸| 91在线观看av| 在线天堂中文资源库| 天堂影院成人在线观看| 午夜免费观看网址| 精品久久久久久成人av| 97人妻精品一区二区三区麻豆 | 两个人看的免费小视频| 国产又色又爽无遮挡免费看| 最新在线观看一区二区三区| 午夜福利高清视频| 1024手机看黄色片| 动漫黄色视频在线观看| 久久精品成人免费网站| 日韩精品中文字幕看吧| 午夜亚洲福利在线播放| 中文资源天堂在线| or卡值多少钱| 精品久久久久久久人妻蜜臀av| 一进一出好大好爽视频| 国产久久久一区二区三区| 在线观看一区二区三区| 亚洲成a人片在线一区二区| 亚洲全国av大片| 18禁黄网站禁片午夜丰满| 18禁观看日本| 这个男人来自地球电影免费观看| 天天一区二区日本电影三级| 中文字幕另类日韩欧美亚洲嫩草| 我的亚洲天堂| 1024香蕉在线观看| www.精华液| 精华霜和精华液先用哪个| 成年女人毛片免费观看观看9| 嫩草影视91久久| 美女 人体艺术 gogo| 天天躁夜夜躁狠狠躁躁| 一本一本综合久久| 久久精品91蜜桃| 欧美乱妇无乱码| 欧美日韩黄片免| 男人操女人黄网站| 亚洲一区二区三区色噜噜| 国产精华一区二区三区| 日韩中文字幕欧美一区二区| 在线观看免费日韩欧美大片| 午夜激情av网站| 一夜夜www| 婷婷亚洲欧美| 亚洲国产欧美一区二区综合| 午夜激情av网站| 色播亚洲综合网| 变态另类成人亚洲欧美熟女| 又大又爽又粗| 日本黄色视频三级网站网址| 久久九九热精品免费| 老汉色∧v一级毛片| 国产亚洲av嫩草精品影院| 色哟哟哟哟哟哟| 久9热在线精品视频| 男人舔女人的私密视频| 亚洲人成电影免费在线| www.自偷自拍.com|