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

    基于微量潤滑的多主元高熵合金銑削加工數(shù)值模擬

    2025-07-18 00:00:00吳巖偉劉寅孫興偉楊赫然董祉序張維鋒
    金剛石與磨料磨具工程 2025年3期

    1,2,2,2,2,1,2(1.,沈陽110870)(2.,沈陽110870)(3.,湖南湘潭411100)

    關(guān)鍵詞 高熵合金;熱力耦合;微量潤滑;材料去除機理;銑削力

    中圖分類號 TG58;TG71 文獻標(biāo)志碼A

    文章編號 1006-852X(2025)03-0332-10

    DOI碼 10.13394/j.cnki.jgszz.2023.0265

    收稿日期 2023-12-05修回日期2024-07-13

    多主元高熵合金(high-entropyalloy,HEA)一般是由大于或等于5種金屬或非金屬元素按照等原子比或近似等原子比,經(jīng)真空熔煉、粉末冶金或其他手段形成的一種具有單相固溶體和高熵值的合金[]。HEA具有高強度、高硬度、高耐磨性、強耐腐蝕性和耐回火軟化等性能特點,在特種材料制造領(lǐng)域中應(yīng)用前景廣闊,引起了眾多學(xué)者的廣泛關(guān)注。

    由于HEA具有較高的硬度和屈服強度,在切削過程中易出現(xiàn)加工硬化現(xiàn)象,因此其屬于難加工材料。DOAN等應(yīng)用分子動力學(xué)仿真方法模擬研究了Al-CrCuFeNi高熵合金在常規(guī)切削和超聲波橢圓振動輔助切削(ultrasonic elliptical vibration assisted cutting,UE-VAC)下的力學(xué)響應(yīng),發(fā)現(xiàn)UEVAC下工件的溫度明顯高于常規(guī)切削時工件的溫度,且其平均切削力隨振動頻率和振幅比的增大而減小。CONSTANTIN等[4分析了加工 Al0.6CoCrFeNi 高熵合金對切削刀具磨損的影響及加工區(qū)域溫度和硬度的演變行為,并與304不銹鋼的加工情況進行了對比,從切削力的角度分析,發(fā)現(xiàn)Al0.6CoCrFeNi 高熵合金的切削力比304不銹鋼的切削力高約 59% 。ZHANG等[5]采用高頻一維超聲振動輔助切削金剛石方法,對FeCrCoMnNi高熵合金進行了微納米級材料去除機理的研究,發(fā)現(xiàn)與傳統(tǒng)金剛石切削相比,其產(chǎn)生的切屑更薄,切削力更低,刀具磨損更輕微,工件加工表面質(zhì)量更好。以上研究均獲得了HEA較低的銑削力、更輕微的刀具磨損和更好的加工表面質(zhì)量,但成本較為高昂,能耗較高,未擺脫傳統(tǒng)澆注式潤滑的弊端,且UEVAC下加工區(qū)域溫度的大幅升高還會降低工件加工精度。

    為此,采用微量潤滑對HEA進行銑削實驗,以期獲得更綠色的加工環(huán)境,更小的銑削力以及更高的加工表面質(zhì)量。微量潤滑(minimalquantitylubrication,MQL)技術(shù)是一種新型的、高效的、低污染的切削潤滑冷卻方式,是將壓縮氣體與極微量的潤滑油混合并汽化后形成微米級的油霧液滴,然后將其高速噴向切削區(qū)域進行潤滑和冷卻[。MQL技術(shù)適用于難加工材料的切削加工,具有大幅降低切削液使用量,減少切削液對非切削區(qū)域污染,降低切削力、切削溫度,防止粘連,延長刀具壽命,提高加工表面質(zhì)量等優(yōu)點[-10]

    HEA擁有廣闊的應(yīng)用前景,但在改善其加工性能方面存在困難。為解決CoCrFeNiMn高熵合金的難加工問題,采用MQL技術(shù)對其進行銑削,從材料去除機理的角度對銑削力進行分析并與干式銑削時的對比,且對不同銑削參數(shù)對銑削力的影響規(guī)律進行研究。

    1CoCrFeNiMn銑削有限元模型

    1.1 四刃立銑刀建模

    實驗使用的四刃立銑刀有限元模型如圖1所示,銑刀模型的主要參數(shù)如表1所示。

    圖1銑刀有限元模型

    Fig.1Finiteelementmodelformillingcutter

    表1銑刀模型主要參數(shù)

    Tab.1Mainparametersof millingcuttermodel

    1.2 CoCrFeNiMn本構(gòu)模型

    本構(gòu)模型是反映物質(zhì)宏觀性質(zhì)的數(shù)學(xué)模型,在固體力學(xué)中用來描述材料應(yīng)力與應(yīng)變之間的關(guān)系。選用Johnson-Cook(J-C)本構(gòu)模型[1]可以較好地描述材料的變形行為。J-C本構(gòu)模型表達式為:

    式中: A 為材料的初始屈服應(yīng)力, B 為材料應(yīng)變硬化模量, ε 為準(zhǔn)靜態(tài)閾值, n 為材料應(yīng)變硬化指數(shù), c 為

    材料應(yīng)變率強化參數(shù),和 分別為應(yīng)變率和參考應(yīng)變率, m 為材料熱軟化指數(shù), T 、 Tm 和 Tr 分別為實時溫度、融化溫度和室溫。

    CoCrFeNiMn的J-C本構(gòu)模型參數(shù)[12]如表2所示。

    表2CoCrFeNiMn的J-C本構(gòu)模型參數(shù)

    Tab.2J-Cconstitutivemodel parametersofCoCrFeNiMn

    1.3材料失效分離準(zhǔn)則

    J-C模型的失效分離準(zhǔn)則屬于物理分離準(zhǔn)則,提供了材料到達失效點時等效塑性應(yīng)變的計算方法,是一種比較適用于材料高應(yīng)變率變形的分離準(zhǔn)則。J-C剪切失效模型根據(jù)單元積分點處的等效塑性應(yīng)變值是否達到失效應(yīng)變來判斷材料是否失效,即當(dāng)破壞參數(shù) ω> 1時,材料網(wǎng)格單元發(fā)生失效并被去除,其表達式為[13]:

    式中: 為初始應(yīng)變, 為等效塑性應(yīng)變增量, 為最大失效應(yīng)變。

    最大失效應(yīng)變表達式為:

    式中: p 為壓應(yīng)力, q 為米塞斯應(yīng)力, 為塑性應(yīng)變率, , d1~d5 為材料的失效參數(shù)。

    1.4接觸摩擦模型

    在銑削過程中,切削刃的前刀面會與切屑產(chǎn)生摩擦,后刀面則會與已加工表面產(chǎn)生摩擦。因此,確定的摩擦模型對切屑形貌、銑削力以及應(yīng)力場和溫度場的仿真分析有很大影響。由于刀尖周圍的熱量不易釋放,工件切屑容易黏附在其表面,所以刀尖周圍主要發(fā)生黏結(jié)摩擦,而遠離刀尖的部位會發(fā)生滑動摩擦。摩擦力 τ 的表達式為[14:

    式中: μ 為摩擦系數(shù), τcrit 為剪切力, σ1 為壓力,μτ<τcrit 表示處于滑動摩擦區(qū), μσ?τcrit 表示處于黏結(jié)摩擦區(qū)。

    在MQL銑削條件下,潤滑基本只作用于滑動摩擦區(qū),且相比于無潤滑的干式銑削,MQL能顯著降低滑動摩擦區(qū)的摩擦系數(shù) μ 。仿真過程中MQL銑削和干式銑削在滑動摩擦區(qū)的摩擦系數(shù) μ[15] 如表3所示。

    Tab.3 Friction coefficients underdifferentlubrication modes

    1.5熱傳導(dǎo)模型

    在銑削過程中,切削刃會對工件和切屑進行剪切、擠壓和摩擦,使得已加工表面和切屑發(fā)生彈性和塑性變形。大部分機械能會被轉(zhuǎn)化成熱能,產(chǎn)生大量的熱量,從而使銑削區(qū)域的溫度大幅提升。熱傳導(dǎo)模型是建立在導(dǎo)熱微分方程和熱力學(xué)第一定律基礎(chǔ)上的,由于后刀面與已加工表面摩擦產(chǎn)生的切削熱相較于其他熱量很少,因此可將其忽略。熱傳導(dǎo)模型微分方程為[16]:

    式中: T 為實時溫度; t 為時間; λ 為熱導(dǎo)率; ρ 為密度; ∣c∣ 為比熱容; , 為熱源在坐標(biāo)軸上的速度分量; Wh 為非彈性熱份額,一般取0.9; 為等效應(yīng)力;ε為等效應(yīng)變速率; J 為熱功當(dāng)量系數(shù)。

    初始條件為:

    式中, T0 為初始溫度。

    邊界條件為:

    式中: w(t) 為刀具與切屑接觸面的熱流密度, 為溫度在刀具與切屑接觸面垂直方向的偏微分。

    且有:

    式中: h 為工件或刀具與空氣接觸面的膜層散熱系數(shù), 為溫度在工件或刀具對流換熱面的偏微分, Ta 為工件或刀具表面溫度。

    相比于MQL銑削,干式銑削時工件和刀具的為0。在MQL銑削條件下,仿真過程中的熱傳導(dǎo)模型參數(shù)[15]如表4所示。

    表3不同潤滑方式下的摩擦系數(shù)

    表4熱傳導(dǎo)模型參數(shù)

    Tab.4 Parametersof heat conduction model

    1.6 銑削模型

    銑刀的旋轉(zhuǎn)方向和進給方向如圖2a所示,其側(cè)刃與端刃共同參與銑削。在有限元仿真中,MQL的作用主要體現(xiàn)在摩擦系數(shù)和膜層散熱系數(shù)上,且作用于與切削刃接觸的2個表面(圖2b)。對工件進行六面體網(wǎng)格劃分,其近似網(wǎng)格尺寸為 0.018mm ;對銑刀進行四面體網(wǎng)格劃分,其近似網(wǎng)格尺寸為 0.300mm ○

    圖2銑削模型示意圖

    Fig.2Schematicdiagramofmillingmodel

    2HEA銑削仿真結(jié)果及分析

    對HEA工件進行銑削仿真實驗,仿真時進給速度為 800.0mm/min ,主軸轉(zhuǎn)速為 45000r/min ,銑削深度為0.10mm 。銑削工件的規(guī)格為 2.5mm×0.3mm×0.3mm ○

    2.1干式銑削與MQL銑削時的米塞斯應(yīng)力場對比

    圖3為干式與MQL銑削時的米塞斯應(yīng)力云圖。如圖3所示:2種方式下等效應(yīng)力較大的區(qū)域主要集中在第一變形區(qū),即在切削刃前刀面即將形成切屑的部位。這是由于此處的工件材料即將發(fā)生剪切滑移,產(chǎn)生的塑性變形程度最大。

    圖3干式和MQL銑削時的米塞斯應(yīng)力云圖

    Fig.3Mises stress nephograms during dry and MQL milling

    高速運動的切削刃初始接觸工件材料時,對工件材料施加壓力,工件材料發(fā)生劇烈的彈性和塑性變形,其等效應(yīng)力集中于端刃和側(cè)刃周圍。隨后產(chǎn)生切屑時,等效應(yīng)力集中于第一變形區(qū),且切削刃前刀面上的切屑逐漸滿足J-C剪切失效準(zhǔn)則,逐漸脫離工件并在壓力的驅(qū)使下飛濺,但此時的切屑仍然存在殘余應(yīng)力。在切削刃離開切削部位時,已加工表面也仍然存在殘余應(yīng)力。

    干式銑削在切削刃初始接觸工件材料時,其等效應(yīng)力除了集中于切削刃周圍工件材料上外,還集中于后刀面靠近切削刃的工件材料上,其值可達 1000MPa 而MQL銑削時后刀面靠近切削刃的工件材料上的等效應(yīng)力僅為 400MPa 左右。這是由于干式銑削的后刀面與材料之間的摩擦系數(shù)較大。另外,MQL銑削在第一變形區(qū)的等效應(yīng)力值可達 1000MPa ,且其等效應(yīng)力集中范圍略大于干式銑削時的。這是由于干式銑削的第一變形區(qū)會使材料受到熱軟化效應(yīng),從而易于切削,

    產(chǎn)生較低的等效應(yīng)力。

    2.2干式銑削與MQL銑削時的溫度場對比

    由于HEA具有極高的強度和硬度,其切削過程中的大部分機械能會被轉(zhuǎn)化為熱能,使得切削區(qū)域的溫度大幅提升,加劇了刀具的磨損;且使工件材料發(fā)生過大的塑性變形,影響了其加工精度;而銑削區(qū)域產(chǎn)生的大量熱量也會反過來影響材料的加工性能,使得材料受到熱軟化效應(yīng),變得易于切削。2種銑削下的溫度云圖如圖4所示。

    由圖4可知:高速運動的切削刃初始接觸工件材料時,工件材料發(fā)生劇烈的彈性和塑性變形,大部分塑性功轉(zhuǎn)化為熱能,瞬時產(chǎn)生大量熱量并聚集于切削刃周圍;隨后產(chǎn)生切屑時,熱量主要集中于材料發(fā)生剪切滑移的位置(第一變形區(qū))及切屑上。這是由于第一變形區(qū)的材料受熱量影響發(fā)生塑性變形和網(wǎng)格單元畸變的程度最為劇烈,切屑與前刀面產(chǎn)生了劇烈的摩擦,且隨著切削的進行,第一變形區(qū)發(fā)生剪切滑移的材料會不斷轉(zhuǎn)變?yōu)榍行?;切屑滿足J-C剪切失效準(zhǔn)則后會斷裂并飛濺,帶走銑削產(chǎn)生的大部分熱量,但已加工表面仍然殘余部分熱量。

    圖4干式和MQL銑削時的溫度云圖

    MQL銑削時已加工表面和切屑的溫度顯著低于干式銑削時的,其降幅達 30%~50% ,這是由于已加工表面的大部分熱量被空氣、霧化的切削液以及切屑帶走。同時,MQL銑削時切屑的完整度更高。這是由于MQL銑削降低了前刀面與切屑間的摩擦系數(shù)[15],因而減小了切屑的塑性變形,并消除了加工過程中產(chǎn)生的部分加工硬化。

    2.3 銑削力仿真

    選取立銑刀端面中心點,通過仿真軟件輸出得到X,Y,Z3 個坐標(biāo)軸方向的作用力,由此繪制出仿真過程中 X,Y,Z3 個坐標(biāo)軸方向的銑削力圖。仿真時進給速度為 800.0mm/min ,主軸轉(zhuǎn)速為 45000r/min ,銑削深度為 0.10mm 。

    圖5為干式銑削下 X,Y,Z3 個坐標(biāo)軸方向的銑削力分布。如圖5所示: X 軸方向上的銑削力最大;同時,由于銑削總進給量較小,其本質(zhì)上為單刃切削,切削刃和工件材料的接觸存在真空期,所以銑削力呈近似周期性分布;四刃立銑刀在1個周期內(nèi),共產(chǎn)生了4次切削行為,切削力逐漸降低。這是由于首次切削產(chǎn)生的熱量已使材料熱軟化,且已去除部分網(wǎng)格單元;在后續(xù)切削中,較小的每齒進給量導(dǎo)致去除的網(wǎng)格單元少于首次切削時的,且切削刃接觸的為已受熱軟化的材料,所以產(chǎn)生的切削力逐漸降低。

    Fig.4Temperature nephogramsduring dryand MQL milling

    圖5干式銑削下 X,Y,Z3 個坐標(biāo)軸方向的銑削力分布 Fig.5Distribution of drymilling forces in ΛX,Y,Z threeaxisdirections

    2.4冷卻潤滑方式及銑削參數(shù)對銑削力的影響

    由于 X 軸方向上的銑削力最大,因而對此方向的銑削力進行單獨分析。網(wǎng)格單元畸變會導(dǎo)致銑削力發(fā)生瞬時突變,因此需將銑削總歷程時間內(nèi)均勻分布時間點上產(chǎn)生的 X 軸方向的銑削力數(shù)據(jù)進行處理,剔除極大值和極小值后取平均值,作為銑削力的研究對象。

    2.4.1干式銑削與MQL銑削時的銑削力對比

    MQL銑削可大幅降低切削刃與工件材料之間的摩擦系數(shù),使銑削更加順暢,從而起到降低銑削力的作用。同時,MQL銑削會大幅提升切削區(qū)域的膜層散熱系數(shù)[15],使工件材料和切削刃表面溫度大幅降低,避免高溫對工件材料產(chǎn)生的熱軟化效應(yīng),又起到增大銑削力的作用。

    圖6為不同銑削深度下干式銑削與MQL銑削時的 X 軸方向平均銑削力對比,銑削時進給速度為400.0mm/min ,主軸轉(zhuǎn)速為 9550r/min 。如圖6所示:在銑削深度為 0.15~0.20mm 時,MQL銑削與干式銑削的平均銑削力相當(dāng)。這是因為MQL銑削減小了刀具與工件材料之間的摩擦系數(shù),因而降低了銑削力,但MQL銑削提升了待去除材料的膜層散熱系數(shù),導(dǎo)致待去除材料受到比干式銑削時更弱的熱軟化效應(yīng),又增大了銑削力;二者相互抵消導(dǎo)致MQL銑削與干式銑削的平均銑削力基本相當(dāng)。同時,在銑削深度 >0.20mm 時,MQL銑削降低刀具與工件材料之間摩擦系數(shù)的因素占主導(dǎo)地位,其降低銑削力的能力隨著銑削深度的增加而增強。在銑削深度為 0.30mm 時,MQL銑削可降低約 30% 的銑削力。

    圖6不同銑削深度下的 X 軸方向平均銑削力對比 Fig. 6 Comparison ofaveragemilling forces in X -axis directionatdifferentmillingdepths

    2.4.2MQL銑削參數(shù)對銑削力的影響及實驗驗證

    圖7是進給速度對 X 軸方向平均銑削力的影響圖,銑削時主軸轉(zhuǎn)速為 45000r/min ,銑削深度為 0.10mm 。圖8是主軸轉(zhuǎn)速對 X 軸方向平均銑削力的影響圖,銑削時進給速度為 500.0mm/min ,銑削深度為 0.10mm 。由圖7和圖8可知: X 軸方向的平均銑削力隨進給速度的增大而增加,隨著主軸轉(zhuǎn)速的升高而降低。這是由于進給速度越大、主軸轉(zhuǎn)速越低,則每齒進給量越大,立銑刀側(cè)刃前刀面受到的來自工件網(wǎng)格單元的作用力也就越大。

    圖7進給速度對 X 軸方向平均銑削力的影響

    Fig.7Influencesof feed rates onaveragemillingforces in

    圖8主軸轉(zhuǎn)速對 X 軸方向平均銑削力的影響 Fig.8Influences of spindle speeds on average milling forces in X-axisdirection

    隨后對仿真實驗中的進給速度對 X 軸方向平均銑削力的影響規(guī)律進行實驗驗證。實驗使用的立式四軸聯(lián)動加工中心主軸轉(zhuǎn)速范圍為 50~8000r/min ,進給速度范圍為 0.01~5000.00mm/min? 。壓電式三分力測力儀量程為 50kN ,銑削力通過數(shù)據(jù)采集卡傳輸?shù)綌?shù)據(jù)采集系統(tǒng)。四刃硬質(zhì)合金立銑刀總長為 50mm ,刃長為6mm ,其他參數(shù)如表1所示,其整體材質(zhì)為硬質(zhì)合金,銑削刃表面附著有CVD涂層。三分力測量儀被牢固地吸附在真空吸盤上,實驗試件通過夾具固定在三分力測量儀的上方。實驗設(shè)備如圖9所示。

    實驗試件為等原子比CoCrFeNiMn高熵合金,其是一種典型的共晶面心立方結(jié)構(gòu)金屬,具有極高的硬度和屈服極限、良好的強塑性匹配能力、優(yōu)異的耐磨性能以及熱穩(wěn)定性,且在室溫下表現(xiàn)出良好的塑性變形能力,擁有廣闊的應(yīng)用前景。實驗試件如圖10所示。試件經(jīng)場發(fā)射掃描電鏡能譜分析的結(jié)果如圖11所示,測得試件中Co、 Cr Fe、 ΔNi 、Mn元素原子數(shù)量之比近似為 1:1:1:1:1:1 。

    圖9實驗設(shè)備

    試件銑削方式為順銑,銑削時主軸轉(zhuǎn)速為 5000r/min 進給速度為 0.4~1.2mm/min ,銑削深度為 1.00mm 。垂直于進給方向的坐標(biāo)軸方向銑削力最大,與仿真時的相同( X 軸方向,即切向方向)。將穩(wěn)定銑削段均勻分布時間點上產(chǎn)生的 X 軸方向銑削力數(shù)據(jù)進行濾波處理后,剔除其最大值和最小值后取平均值,作為銑削力的研究對象。每次銑削完后均更換一次刀具,每組銑削實驗重復(fù)5次,取平均值以提高結(jié)果的準(zhǔn)確性。圖12為進給速度對 X 軸方向平均銑削力的影響。如圖12所示: X 軸方向的平均銑削力隨進給速度的增大而增大,與仿真實驗的結(jié)論相符(圖7)。

    圖10實驗試件

    圖12進給速度對 X 軸方向平均銑削力的影響Fig.12Influencesof feed rates on averagemilling forcesinX-axisdirection

    Fig.9Test equipment

    Fig.10 Test specimen

    圖11實驗試件能譜分析

    Fig.11Energy spectrum analysis of test specimen

    圖13為不同銑削深度下每齒進給量對 X 軸方向平均銑削力的影響。如圖13所示:在每齒進給量一定的前提下, X 軸方向的平均銑削力隨著銑削深度的增大而增加;且隨銑削深度的增大,每齒進給量對 X 軸方向平均銑削力的影響程度逐漸增強。

    圖13不同銑削深度下每齒進給量對 X 軸方向平均銑削力的影響

    Fig.13Effectsof feed per tooth onaverage milling forces in X -axisdirectionatdifferentmillingdepths

    3結(jié)論

    建立了基于CoCrFeNiMn高熵合金和四刃立銑刀的熱力耦合三維單刃銑削模型,通過分析材料去除機理以研究MQL銑削與干式銑削在銑削力上的差異,隨后通過單因素實驗研究不同銑削參數(shù)對平均銑削力的影響規(guī)律。得出如下結(jié)論:

    (1)MQL銑削降低了刀具與工件材料間的摩擦系數(shù),減小了銑削力;但MQL銑削提升了待去除材料的膜層散熱系數(shù),導(dǎo)致待去除材料受到比干式銑削時更弱的熱軟化效應(yīng),增大了第一變形區(qū)的等效應(yīng)力值和聚集范圍,又增大了銑削力。二者相互抵消導(dǎo)致MQL銑削與干式銑削在銑削深度為 0.15~0.20mm 時的銑削力基本相當(dāng)。但在銑削深度 >0.20mm 時,MQL銑削降低銑削力的能力隨銑削深度的增大而增強;在銑削深度為 0.30mm 時,MQL銑削可降低約 30% 的銑削力??傮w來看,MQL銑削優(yōu)勢明顯,可在銑削深度 > 0.20mm 時顯著降低銑削力,同時顯著降低銑削部位的溫度,以此提高工件加工精度并規(guī)避了傳統(tǒng)澆注式潤滑的一系列缺陷。

    (2)平均銑削力隨進給速度的增大而增大,隨主軸轉(zhuǎn)速的升高而降低,即隨每齒進給量的增大而增加。平均銑削力隨銑削深度的增大而增大,且隨銑削深度增加,每齒進給量增加對銑削力的影響程度逐漸加劇。因而在銑削深度較大時,應(yīng)降低每齒進給量,以應(yīng)對銑削力對銑削深度的高敏感性。

    參考文獻:

    [1] 梁秀兵,魏敏,程江波,等.高熵合金新材料的研究進展[J].材料工程, 2009,37(12):75-79. LIANGXiubing,WEI Min, CHENG Jiangbo,etal.Research progress in advanced materials ofhigh-entropy alloys [J].Journal of Materials Engineering,2009,37(12):75-79.

    [2] 陳永星,朱勝,王曉明,等.高熵合金制備及研究進展[J].材料工程, 2017,45(11):129-138. CHENYongxing,ZHU Sheng,WANG Xiaoming,etal.Progress in preparationand research of high entropy alloys[J].Journal of Materials Engineering,2017,45(11):129-138.

    [3] DOAN D Q,F(xiàn)ANG T H,CHEN T H. Machining mechanism and deformation behavior of high-entropy alloy underelliptical vibration cutting[J].Intermetallics,2021,131(11): 107079.

    [4] CONSTANTING,BALANE,VOICULESCUI, etal.Cuttingbehavior ofAl0.6CoCrFeNihigh entropyalloy[J].Materials,2020,13(18):4181.

    [5] ZHANGL,HASHIMOTOT,YANJW.Machinability exploration for high-entropy alloyFeCrCoMnNi byultrasonic vibration-assisted diamond turning[J]. CIRP Annals,2021,70(5): 37-40.

    [6] 袁松梅,韓文亮,朱光遠,等.綠色切削微量潤滑增效技術(shù)研究進展[J]. 機械工程學(xué)報,2019,55(5):175-185. YUAN Songmei,HAN Wenliang,ZHU Guangyuan,et al.Recent progress on the efficiency increasing methods of minimum quantity lubrication technology in green cutting [J].Journal of Mechanical Engineering,2019,55(5): 175-185.

    L7」 哀松梅,朱光遠,土利.綠色切削微量潤滑技木潤滑劑特性研究進展 [J].機械工程學(xué)報,2017,53(17):131-140. YUAN Songmei, ZHU Guangyuan,WANG Li. Recent progress on lubricant characteristicsof minimum quantity lubrication (MQL) technology in green cutting [J]. Journal of Mechanical Engineering,2017, 53(17): 131-140.

    [8] WANGYG,LICH, ZHANGYB, et al.Experimental evaluation of the lubrication properties of the wheel/workpiece interface in minimum quantity lubrication (MQL) grindingusing different typesof vegetable oils[J].Journal ofCleanerProduction,2016,127: 487-499.

    [9] SUY,GONGL,LI B,etal.Performanceevaluationof nanofluidMQL with vegetable-based oil and ester oil as base fluids in turning [J].The International Journal of Advanced Manufacturing Technology,2016, 83(9): 2083-2089.

    [10]WANG Y G, LI C H, ZHANG Y B,et al. Experimental evaluation of the lubrication propertiesof the wheel/workpiece interface in MQL grinding with different nanofluids[J].Tribology International,20l6,99:198-10.

    [11]周琳,王子豪,文鶴鳴.簡論金屬材料JC本構(gòu)模型的精確性[J].高壓 物理學(xué)報,2019,33(4):3-16. ZHOU Lin,WANG Zihao,WEN Heming.On the accuracy of the Johnson-Cook constitutive model for metals[J]. Chinese Journal of High Pressure Physics,2019,33(4): 3-16.

    [12]許海濤.動態(tài)預(yù)壓縮對CoCrFeNiMn高熵合金微尺度壓入硬度的影 響[D].太原:太原理工大學(xué),2021. XUHaitao.Effect of dynamic pre-compression micro-scale indentation hardnessof CoCrFeNiMnhigh-entropyalloy [D].Taiyuan:Taiyuan University ofTechnology,2021.

    [13]尹曉珊,鐘建琳,彭寶營,等.Ti2AINb 基合金切削力有限元仿真[J].工 具技術(shù),2023,57(4):65-69. YINXiaoshan,ZHONG Jianlin,PENG Baoying,etal.Finite element simulationofcutting force of TiAINb-based alloy[J].Tool Engineering, 2023,57(4): 65-69.

    [14]ZHANG Y,HUANG X. Simulation and research of aerospace material miling based on ABAQUS [J].Journal of Physics: Conference Series, 2020,1653(1): 012070.

    [15]梁正一,楊赫然,劉寅,等.基于MQL 的鋯基塊體金屬玻璃鉆削仿真 研究[J].工具技術(shù),2023,57(7):90-96. LIANG Zhengyi, YANG Heran, LIU Yin, et al. Simulation research on driling of Zr-based bulk metal glass based on MQL [J].Tool Engineering,2023,57(7): 90-96.

    [16]奚欣欣,陳濤,丁文鋒.TiAl合金低壓渦輪葉片榫頭磨削溫度場研 究[J].金剛石與磨料磨具工程,2020,40(5):17-22. XI Xinxin, CHEN Tao,DING Wenfeng.Research on temperature field duringgrinding of low-pressure turbine blade tenon of TiAl alloys [J]. Diamond Abrasives Engineering,2020,40(5):17-22.

    作者簡介

    劉寅,男,1986年生,副教授、博士。主要研究方向:難加工材料精密加工技術(shù),微尺度加工技術(shù),綠色與清潔加工技術(shù)。E-mail:liuyin_neu@163.com

    (編輯:周萬里)

    Numerical simulation of multi-principal elements high-entropy alloy milling based on minimal quantity lubrication

    WU Yanwei13, LIU Yin12, SUN Xingwei12,YANG Heran12, DONG Zhixu12, ZHANG Weifeng12 (1.School ofMechanical Engineering,Shenyang University ofTechnology,Shenyang 110870,China) (2.Liaoning Province Key Laboratory of Complex Surface NC Manufacturing Technology, Shenyang 110870,China) (3. Jianglu Machinery Electronics Group Co.,Ltd., China North Industries Group Co.,Ltd.,Xiangtan 4111oo, Hunan, China)

    AbstractObjectives: High-entropy alloys (HEA) with multi-principal elements have many excellent characteristics suchashigh strength,highhardnessandhighwearresistance,whichhaveatracted widespreadattentionfromscholars. However, the mechanical procesing ofHEA isdificult.This studyutilizes minimalquantity lubrication (MQL)technologyfor millng HEA to improve their machining performanceand explore the influences of diffrent millng parameters on their miling force. Methods: The thermodynamic coupling miling model of HEA (CoCrFeNiMn) and a fouredge end millng cutter is establishedusing finite element simulation software.The difference in milling force between MQL millingand dry milling is studied byanalyzing the materialremoval mechanism,and the influencesofdiferent milingparameters onmiling forces are studied bysingle-factor experiments.Firstly,the three-dimensional modelof the four-edgeend milling cutter,the J-C constitutive model of HEA (CoCrFeNiMn),the heat conduction model and the contact friction modelareestablished,andthe material failure separationcriteriaaredeterminedandthemiling model is established.Then,the material removal mechanism under the thermodynamic coupling condition is analyzed,and the changes in milling forces between MQL milling and dry miling are analyzed from the perspectives of equivalent stress, contact friction and the thermal softening effct,and the advantagesofMQL technology formiling HEA (CoCrFeNiMn)are analyzed.Finally,the effects of feed speed,spindle speed and millng depthon milling force are obtained by single-factor tests.Results:Through comparative experiments between dry millng and MQL miling,it is found that: (1) The equivalent stress produced bythe two miling methods are concentrated in the first deformation Zone,and the equivalent stress ofdrymiling is also concentrated inthe position near the cuting edge.(2)The equivalent stress value of MQL milling in the firstdeformation zoneis slightly greater thanthatofdry millng.(3)The heat generated bythe two millng methods isconcentrated inthe first deformation zone andthe chips,andthe chips takeawaymostof the heat.Thechip temperature generatedby dry milling is significantly higher than that generated by MQL miling. (4) MQ miling significantly reduces the temperature at the cutting site and improves chip integrity. (5) When the milling depth is 0.15 to 0.20mm ,the milling force of MQL milling is basically the same as that of dry milling. When the milling depth is greater than 0.20mm , the ability of MQL milling to reduce milling force increases with the increase of millng depth.This is due to theuse ofMQL technology in the miling process, which compensates for the reduced millng forcecaused bylowfrictioncoeficientand the increased milling force caused by weak thermal softening effect. The MQL milling single-factor tests show that: (1)The milling force increases with the increaseoffeed speed,and decreases withthe increaseof spindle speed,that is,itincreases with the increase offeed rate per toth.(2)The milling force increases withthe increase of milling depths,and the effect offeed rate per tooth on the average milling force is graduallyintensified with the increase of milling depth. Conclusions: The MQL miling has obvious advantages,which can significantly reduce milling force at milling depth greater than 0.20mm ,and the ability to reduce milling force increases with the increase of milling depth.Inaddition,the useof MQL technology in the milling processcan significantly reduce the temperatureatthecuting position,improve HEA machiningaccuracy,andavoidaseries of defects of traditional pouring lubrication.The cuttng force in MQL milling conforms to the change law of cuttng force in most metal milling with the process parameters. In MQL milling,the miling force canbe further reduced by increasing the spindle speedand reducing the feedspeed.When the miling depth is large,the spindle speed should be further increased and the feed speed should bereduced to reduce the feed per toth,inorderto deal with thehigh sensitivity of the milling force to the milling depth.

    Key wordshigh-entropy alloy(HEA);thermo-mechanical coupling; minimal quantity lubrication (MQL);material removal mechanism; milling force

    天天操日日干夜夜撸| 老熟女久久久| 国产免费一区二区三区四区乱码| 亚洲国产中文字幕在线视频| 十八禁人妻一区二区| 午夜福利免费观看在线| 热99久久久久精品小说推荐| www.自偷自拍.com| 亚洲中文字幕日韩| 精品高清国产在线一区| 国产高清videossex| 热99re8久久精品国产| 大型av网站在线播放| a级毛片在线看网站| 中亚洲国语对白在线视频| 亚洲美女黄色视频免费看| 天天操日日干夜夜撸| 国产av精品麻豆| 日韩,欧美,国产一区二区三区| h视频一区二区三区| 午夜免费成人在线视频| 欧美日韩av久久| 高清欧美精品videossex| 国产1区2区3区精品| 黄色a级毛片大全视频| 亚洲欧洲日产国产| 首页视频小说图片口味搜索| 国产欧美日韩一区二区精品| 精品熟女少妇八av免费久了| 欧美日韩福利视频一区二区| 亚洲精品自拍成人| 国产精品一区二区在线观看99| 亚洲国产精品999| 两个人免费观看高清视频| 中文字幕制服av| 精品久久久精品久久久| 我的亚洲天堂| 一级毛片电影观看| av又黄又爽大尺度在线免费看| 精品免费久久久久久久清纯 | 狠狠狠狠99中文字幕| 亚洲av成人不卡在线观看播放网 | 亚洲精品美女久久av网站| 国产一卡二卡三卡精品| 日韩欧美一区二区三区在线观看 | 一个人免费在线观看的高清视频 | 欧美成狂野欧美在线观看| 欧美精品一区二区大全| 高清黄色对白视频在线免费看| 久久人妻熟女aⅴ| 国产成人精品久久二区二区91| 久久精品国产亚洲av高清一级| 免费高清在线观看日韩| 午夜福利免费观看在线| 亚洲国产日韩一区二区| 亚洲国产日韩一区二区| 欧美日本中文国产一区发布| 亚洲国产日韩一区二区| 亚洲中文av在线| 美女高潮到喷水免费观看| 19禁男女啪啪无遮挡网站| 99精品欧美一区二区三区四区| 亚洲男人天堂网一区| 91精品三级在线观看| 精品少妇黑人巨大在线播放| 宅男免费午夜| 人妻一区二区av| 久久久久精品人妻al黑| 捣出白浆h1v1| 成人亚洲精品一区在线观看| 90打野战视频偷拍视频| 俄罗斯特黄特色一大片| 日韩视频一区二区在线观看| 国产亚洲欧美在线一区二区| 新久久久久国产一级毛片| 777久久人妻少妇嫩草av网站| 午夜成年电影在线免费观看| 两人在一起打扑克的视频| 男人操女人黄网站| 91国产中文字幕| 日韩大码丰满熟妇| a级毛片在线看网站| 亚洲欧美一区二区三区黑人| 叶爱在线成人免费视频播放| 欧美日韩中文字幕国产精品一区二区三区 | 久久热在线av| 亚洲av电影在线观看一区二区三区| 欧美少妇被猛烈插入视频| av一本久久久久| 精品一区在线观看国产| 国产成人啪精品午夜网站| 男女高潮啪啪啪动态图| 国产成人啪精品午夜网站| 2018国产大陆天天弄谢| 18禁裸乳无遮挡动漫免费视频| 91麻豆av在线| 亚洲精品国产一区二区精华液| 亚洲精品国产一区二区精华液| 女人精品久久久久毛片| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕色久视频| 老司机午夜福利在线观看视频 | 亚洲国产日韩一区二区| 亚洲精品一二三| 中文字幕人妻丝袜一区二区| 久久久水蜜桃国产精品网| 亚洲av电影在线进入| 成年美女黄网站色视频大全免费| 亚洲一码二码三码区别大吗| 视频在线观看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 香蕉国产在线看| 欧美人与性动交α欧美精品济南到| 久久这里只有精品19| 国产亚洲欧美精品永久| 首页视频小说图片口味搜索| 99久久精品国产亚洲精品| 老司机福利观看| 一级毛片电影观看| 青春草亚洲视频在线观看| 国产成人av激情在线播放| 淫妇啪啪啪对白视频 | 高清欧美精品videossex| 岛国在线观看网站| 亚洲中文字幕日韩| 热re99久久精品国产66热6| 伊人久久大香线蕉亚洲五| 久久久久久久国产电影| 精品亚洲乱码少妇综合久久| 精品国产一区二区三区四区第35| 久久天堂一区二区三区四区| 中文字幕人妻丝袜一区二区| 亚洲av成人一区二区三| 欧美精品av麻豆av| 日本av手机在线免费观看| 在线观看免费高清a一片| 午夜福利影视在线免费观看| 欧美人与性动交α欧美软件| 五月天丁香电影| 一边摸一边抽搐一进一出视频| 日韩 亚洲 欧美在线| 亚洲精品在线美女| 在线亚洲精品国产二区图片欧美| 乱人伦中国视频| 男女无遮挡免费网站观看| 18禁国产床啪视频网站| 欧美日本中文国产一区发布| 国产免费现黄频在线看| 亚洲欧美精品自产自拍| 永久免费av网站大全| av网站免费在线观看视频| 国产精品二区激情视频| 色综合欧美亚洲国产小说| 中文字幕av电影在线播放| 99国产精品一区二区三区| 久久亚洲国产成人精品v| 国产欧美日韩一区二区三区在线| 老汉色∧v一级毛片| 天天躁夜夜躁狠狠躁躁| av在线app专区| 他把我摸到了高潮在线观看 | 999久久久国产精品视频| 乱人伦中国视频| 国产日韩欧美在线精品| 日本91视频免费播放| xxxhd国产人妻xxx| 9色porny在线观看| netflix在线观看网站| 在线十欧美十亚洲十日本专区| 老司机靠b影院| 亚洲国产毛片av蜜桃av| 国产伦理片在线播放av一区| 国产亚洲欧美在线一区二区| 最近中文字幕2019免费版| 亚洲七黄色美女视频| 中亚洲国语对白在线视频| 国产欧美日韩一区二区精品| 女人精品久久久久毛片| 老司机午夜福利在线观看视频 | 久久国产精品大桥未久av| 美女脱内裤让男人舔精品视频| 电影成人av| 老司机午夜十八禁免费视频| 欧美黄色淫秽网站| 国产高清videossex| 日日夜夜操网爽| 丝袜脚勾引网站| 男人舔女人的私密视频| 天天影视国产精品| 亚洲国产精品999| 亚洲欧美一区二区三区黑人| 男人操女人黄网站| 男女国产视频网站| 99国产精品99久久久久| 亚洲va日本ⅴa欧美va伊人久久 | 美女中出高潮动态图| 中文字幕人妻丝袜一区二区| 视频区欧美日本亚洲| 国产福利在线免费观看视频| 国产av又大| 各种免费的搞黄视频| 91大片在线观看| 国产亚洲精品一区二区www | 五月开心婷婷网| 90打野战视频偷拍视频| 亚洲av欧美aⅴ国产| 别揉我奶头~嗯~啊~动态视频 | 中文字幕制服av| 欧美激情极品国产一区二区三区| 日本a在线网址| 男人操女人黄网站| 国产男人的电影天堂91| 日日摸夜夜添夜夜添小说| 午夜日韩欧美国产| 日韩中文字幕欧美一区二区| 无遮挡黄片免费观看| 国产伦理片在线播放av一区| 叶爱在线成人免费视频播放| 一级,二级,三级黄色视频| 精品国产乱码久久久久久小说| 亚洲av电影在线观看一区二区三区| 中文字幕av电影在线播放| av免费在线观看网站| 91字幕亚洲| 久久久欧美国产精品| 国产男女超爽视频在线观看| 最近最新免费中文字幕在线| 黄色片一级片一级黄色片| 男女国产视频网站| 男女边摸边吃奶| 欧美成狂野欧美在线观看| 国产激情久久老熟女| 中文字幕精品免费在线观看视频| 亚洲中文字幕日韩| 精品亚洲成国产av| 欧美老熟妇乱子伦牲交| 免费高清在线观看日韩| 日本av免费视频播放| 中文精品一卡2卡3卡4更新| 丁香六月欧美| 亚洲精品一卡2卡三卡4卡5卡 | 日本五十路高清| 久久中文看片网| 在线看a的网站| 国产av又大| 欧美中文综合在线视频| 精品国内亚洲2022精品成人 | 日韩 欧美 亚洲 中文字幕| 女性生殖器流出的白浆| 又大又爽又粗| 免费久久久久久久精品成人欧美视频| 精品少妇内射三级| 激情视频va一区二区三区| 久久亚洲精品不卡| 老司机影院毛片| 波多野结衣一区麻豆| 正在播放国产对白刺激| 老司机深夜福利视频在线观看 | 国产精品久久久久久人妻精品电影 | 国产在线观看jvid| 国产精品九九99| 亚洲精品久久久久久婷婷小说| 首页视频小说图片口味搜索| 国产成人av激情在线播放| 亚洲自偷自拍图片 自拍| 美女视频免费永久观看网站| 国产免费av片在线观看野外av| 国内毛片毛片毛片毛片毛片| 欧美精品一区二区免费开放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲 欧美一区二区三区| 亚洲专区国产一区二区| 日本91视频免费播放| 黄色视频不卡| 亚洲天堂av无毛| 正在播放国产对白刺激| 欧美日韩精品网址| 国产在线观看jvid| 午夜福利免费观看在线| 欧美成狂野欧美在线观看| 亚洲国产欧美在线一区| 99久久综合免费| 动漫黄色视频在线观看| 1024香蕉在线观看| 天天操日日干夜夜撸| 两人在一起打扑克的视频| 啦啦啦视频在线资源免费观看| 天堂俺去俺来也www色官网| 国产淫语在线视频| 欧美在线一区亚洲| 首页视频小说图片口味搜索| 91精品伊人久久大香线蕉| 在线观看www视频免费| 日本欧美视频一区| 99国产综合亚洲精品| 欧美日韩视频精品一区| 国产精品一区二区在线不卡| 亚洲专区国产一区二区| 天天影视国产精品| 成人免费观看视频高清| 中亚洲国语对白在线视频| 成年女人毛片免费观看观看9 | 亚洲欧美精品综合一区二区三区| 99国产极品粉嫩在线观看| 高清视频免费观看一区二区| 国产成人精品久久二区二区91| 日韩一卡2卡3卡4卡2021年| 91老司机精品| 国产亚洲一区二区精品| 国产一卡二卡三卡精品| 男女下面插进去视频免费观看| 亚洲国产av影院在线观看| 啦啦啦视频在线资源免费观看| 久久久久久免费高清国产稀缺| 99精品欧美一区二区三区四区| 精品福利永久在线观看| 一本色道久久久久久精品综合| 精品久久久精品久久久| 99久久人妻综合| 国产精品.久久久| 中文字幕人妻熟女乱码| 99九九在线精品视频| 日韩欧美国产一区二区入口| 国产精品二区激情视频| 久久性视频一级片| 青春草亚洲视频在线观看| 男人爽女人下面视频在线观看| 精品视频人人做人人爽| 亚洲精品中文字幕一二三四区 | 超色免费av| 国产精品一区二区免费欧美 | 男女边摸边吃奶| 一本—道久久a久久精品蜜桃钙片| 亚洲色图综合在线观看| 热re99久久精品国产66热6| 久久精品aⅴ一区二区三区四区| 国产福利在线免费观看视频| 一本色道久久久久久精品综合| 女人爽到高潮嗷嗷叫在线视频| 国产日韩欧美亚洲二区| 美女脱内裤让男人舔精品视频| 免费高清在线观看视频在线观看| 久久亚洲精品不卡| 99精国产麻豆久久婷婷| 国产国语露脸激情在线看| 国产在视频线精品| 啦啦啦免费观看视频1| 免费女性裸体啪啪无遮挡网站| 777米奇影视久久| 在线 av 中文字幕| 国产淫语在线视频| 久久天躁狠狠躁夜夜2o2o| 欧美另类一区| av不卡在线播放| 人人妻,人人澡人人爽秒播| 亚洲性夜色夜夜综合| 狂野欧美激情性xxxx| 成年人午夜在线观看视频| 欧美黑人精品巨大| 亚洲精品日韩在线中文字幕| www.熟女人妻精品国产| 少妇 在线观看| 久久精品aⅴ一区二区三区四区| 一级,二级,三级黄色视频| 亚洲国产精品一区三区| 欧美 日韩 精品 国产| 亚洲精品美女久久久久99蜜臀| 美女主播在线视频| 亚洲中文日韩欧美视频| 97在线人人人人妻| 大型av网站在线播放| 男女国产视频网站| 久久久久久免费高清国产稀缺| 国产亚洲欧美在线一区二区| 操美女的视频在线观看| 色婷婷久久久亚洲欧美| 在线观看免费视频网站a站| 麻豆av在线久日| 9色porny在线观看| 999久久久国产精品视频| 国产精品自产拍在线观看55亚洲 | videos熟女内射| 男女下面插进去视频免费观看| 亚洲精品国产色婷婷电影| 老司机靠b影院| 黑丝袜美女国产一区| 夜夜夜夜夜久久久久| 欧美精品亚洲一区二区| 午夜福利一区二区在线看| 亚洲人成电影观看| 国产av国产精品国产| 亚洲 国产 在线| 亚洲精品中文字幕一二三四区 | 国产精品亚洲av一区麻豆| 国产亚洲精品一区二区www | 亚洲精品国产一区二区精华液| 侵犯人妻中文字幕一二三四区| 精品乱码久久久久久99久播| 久久精品久久久久久噜噜老黄| netflix在线观看网站| 免费黄频网站在线观看国产| 欧美精品一区二区大全| 丰满迷人的少妇在线观看| 久久精品国产a三级三级三级| 美女大奶头黄色视频| 久久久久久久精品精品| 日本五十路高清| 纵有疾风起免费观看全集完整版| 青青草视频在线视频观看| 91成人精品电影| 1024视频免费在线观看| 日日爽夜夜爽网站| 国产精品久久久av美女十八| 中国美女看黄片| 国精品久久久久久国模美| 日本91视频免费播放| 男女边摸边吃奶| 美女福利国产在线| 美女中出高潮动态图| 香蕉国产在线看| 精品久久久精品久久久| 蜜桃在线观看..| 亚洲av男天堂| 亚洲精品成人av观看孕妇| 国产一区二区三区在线臀色熟女 | 国产精品一区二区在线不卡| 国产精品一二三区在线看| 一个人免费在线观看的高清视频 | 色婷婷久久久亚洲欧美| 十八禁网站免费在线| 99久久精品国产亚洲精品| 日韩一区二区三区影片| 美女脱内裤让男人舔精品视频| 国产片内射在线| 国产亚洲av片在线观看秒播厂| 欧美成狂野欧美在线观看| 美国免费a级毛片| a级片在线免费高清观看视频| 亚洲精华国产精华精| 男女下面插进去视频免费观看| 色视频在线一区二区三区| 亚洲欧美色中文字幕在线| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区久久| 国产伦理片在线播放av一区| 亚洲成国产人片在线观看| 建设人人有责人人尽责人人享有的| 成在线人永久免费视频| 搡老乐熟女国产| 国产又爽黄色视频| 欧美日韩一级在线毛片| 老熟妇乱子伦视频在线观看 | av超薄肉色丝袜交足视频| 色播在线永久视频| 在线永久观看黄色视频| 女人精品久久久久毛片| 国产一级毛片在线| 男女下面插进去视频免费观看| 久久青草综合色| 国产真人三级小视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 大香蕉久久网| 国产男女内射视频| 精品人妻一区二区三区麻豆| 免费高清在线观看视频在线观看| 无限看片的www在线观看| 久久久久精品人妻al黑| 国产精品 欧美亚洲| 好男人电影高清在线观看| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av成人精品| 欧美黑人欧美精品刺激| 中文精品一卡2卡3卡4更新| 伊人久久大香线蕉亚洲五| 丝袜人妻中文字幕| 亚洲欧美一区二区三区久久| 1024香蕉在线观看| 各种免费的搞黄视频| 国产一级毛片在线| 老鸭窝网址在线观看| 国产av国产精品国产| 欧美黑人精品巨大| 爱豆传媒免费全集在线观看| 少妇粗大呻吟视频| 亚洲av国产av综合av卡| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人精品久久二区二区91| 窝窝影院91人妻| 久久精品国产亚洲av香蕉五月 | 巨乳人妻的诱惑在线观看| 一区在线观看完整版| 热99国产精品久久久久久7| 日韩三级视频一区二区三区| 国产精品欧美亚洲77777| av天堂在线播放| 亚洲天堂av无毛| 亚洲午夜精品一区,二区,三区| 国产av精品麻豆| av在线播放精品| 亚洲精品久久成人aⅴ小说| 成人国语在线视频| 12—13女人毛片做爰片一| a级毛片在线看网站| 免费在线观看视频国产中文字幕亚洲 | 在线看a的网站| 欧美乱码精品一区二区三区| 老汉色av国产亚洲站长工具| 亚洲综合色网址| 国产深夜福利视频在线观看| 中文欧美无线码| www日本在线高清视频| 亚洲国产日韩一区二区| av又黄又爽大尺度在线免费看| 久久久久国产精品人妻一区二区| 亚洲第一青青草原| 桃红色精品国产亚洲av| www.精华液| 桃红色精品国产亚洲av| 亚洲一码二码三码区别大吗| www.999成人在线观看| 天天添夜夜摸| 久久久久视频综合| 飞空精品影院首页| 欧美成狂野欧美在线观看| 亚洲精品成人av观看孕妇| 亚洲一码二码三码区别大吗| 黑人操中国人逼视频| 大香蕉久久网| 久久久久久久精品精品| 十八禁网站免费在线| 欧美久久黑人一区二区| 亚洲精品久久成人aⅴ小说| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av高清不卡| 久久人妻熟女aⅴ| 中文字幕高清在线视频| 国产欧美亚洲国产| 人人妻人人添人人爽欧美一区卜| 国产精品久久久人人做人人爽| 国产成人欧美在线观看 | 亚洲熟女精品中文字幕| 午夜两性在线视频| 最新的欧美精品一区二区| 亚洲欧美精品综合一区二区三区| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 19禁男女啪啪无遮挡网站| 一级毛片女人18水好多| 99re6热这里在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲,欧美精品.| 国产精品一区二区免费欧美 | 可以免费在线观看a视频的电影网站| 国产精品影院久久| 国产成人一区二区三区免费视频网站| 少妇粗大呻吟视频| 国产色视频综合| 999精品在线视频| 制服诱惑二区| 动漫黄色视频在线观看| 嫁个100分男人电影在线观看| 又紧又爽又黄一区二区| 国产黄色免费在线视频| 九色亚洲精品在线播放| 亚洲av日韩在线播放| 欧美精品啪啪一区二区三区 | 亚洲免费av在线视频| 涩涩av久久男人的天堂| 三级毛片av免费| 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 极品少妇高潮喷水抽搐| 亚洲免费av在线视频| 国产欧美日韩一区二区三区在线| av福利片在线| 伊人久久大香线蕉亚洲五| 日本wwww免费看| 亚洲精品一区蜜桃| 性少妇av在线| 国产日韩欧美亚洲二区| 成人国产av品久久久| 亚洲专区字幕在线| 少妇裸体淫交视频免费看高清 | www.999成人在线观看| 制服人妻中文乱码| 午夜福利,免费看| 色视频在线一区二区三区| 91精品三级在线观看| 亚洲 欧美一区二区三区| 成年动漫av网址| 一本—道久久a久久精品蜜桃钙片| 男女午夜视频在线观看| 首页视频小说图片口味搜索| 日韩 欧美 亚洲 中文字幕| 午夜成年电影在线免费观看| 成年动漫av网址| www.熟女人妻精品国产| 美女高潮喷水抽搐中文字幕| 精品少妇久久久久久888优播| 国产三级黄色录像| 9191精品国产免费久久| 免费观看人在逋| 日韩 亚洲 欧美在线| 麻豆av在线久日| 国产精品久久久久久精品古装| 免费在线观看黄色视频的| 日本五十路高清| 91麻豆av在线| 国产91精品成人一区二区三区 | 久久女婷五月综合色啪小说| 91九色精品人成在线观看| 久久人妻熟女aⅴ| 亚洲精品一卡2卡三卡4卡5卡 | 搡老乐熟女国产| 欧美日本中文国产一区发布|