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

    細(xì)徑金剛石線鋸鋸縫內(nèi)切削液流場(chǎng)仿真研究

    2024-12-06 00:00:00陳佳虎葛培琪
    金剛石與磨料磨具工程 2024年6期
    關(guān)鍵詞:計(jì)算流體力學(xué)

    摘要 隨著工件大尺寸化及鋸絲細(xì)線化,鋸切加工過(guò)程中的鋸縫越來(lái)越深且窄,切削液在鋸切過(guò)程中無(wú)法充分發(fā)揮作用,對(duì)切片質(zhì)量影響較大?;谟?jì)算流體力學(xué)(computational fluid dynamics,CFD)數(shù)值模擬,通過(guò)建立 CFD 鋸縫模型,對(duì)金剛石線鋸鋸切加工材料時(shí)鋸縫內(nèi)切削液流場(chǎng)進(jìn)行分析研究。仿真分析發(fā)現(xiàn):在小尺寸鋸縫內(nèi),隨著走絲速度增大至 25 m/s,切削液更能充分進(jìn)入鋸縫,在鋸絲與工件接觸區(qū)域及非接觸區(qū)域充滿液體后,接觸區(qū)域流體壓力在 0.179 0 MPa 左右,非接觸區(qū)域流體壓力在 0.159 0 MPa 左右;切削液黏度和表面張力在一定范圍內(nèi)的降低,有利于保證鋸縫內(nèi)切削液的相對(duì)飽和與穩(wěn)定,同時(shí)可以使鋸縫內(nèi)切削液壓力分布更為穩(wěn)定。

    關(guān)鍵詞 金剛石線鋸;計(jì)算流體力學(xué);小尺寸鋸縫;切削液流場(chǎng)

    中圖分類號(hào) TG58; TG74; O357 文獻(xiàn)標(biāo)志碼 A

    文章編號(hào) 1006-852X(2024)06-0781-08

    DOI 碼 10.13394/j.cnki.jgszz.2023.0235

    收稿日期 2023-11-07 修回日期 2023-12-21

    切削液在金剛石線鋸鋸切加工過(guò)程中具有非常重要的作用,其潤(rùn)滑冷卻作用可以減少鋸絲與工件的摩擦,吸收加工過(guò)程中產(chǎn)生的熱量,降低局部加工區(qū)域溫度,防止切片熱變形,提升鋸絲使用壽命[1] 。隨著鋸絲逐漸細(xì)線化,鋸切大尺寸工件時(shí),鋸縫越來(lái)越深且窄,切削液不能大量進(jìn)入鋸切區(qū)域,這使得金剛石線鋸鋸切加工過(guò)程中,切削液的潤(rùn)滑冷卻效果變差,進(jìn)而導(dǎo)致切片表面質(zhì)量下降[2] 。因此,為充分發(fā)揮切削液在鋸切過(guò)程中的作用從而保證切片質(zhì)量,需要對(duì)金剛石線鋸鋸切過(guò)程中切削液在鋸縫內(nèi)的流場(chǎng)進(jìn)行進(jìn)一步研究。

    M?LLER 等[3]提出一種二維模型分析漿料的彈性流體動(dòng)力學(xué)行為,發(fā)現(xiàn)漿液的動(dòng)力學(xué)行為及其與鋸絲的相互作用是影響鋸切效果的重要因素。ZHU 等[4]將鋸絲振動(dòng)特性引入模型中,探究了鋸絲振動(dòng)與流體動(dòng)壓效應(yīng)的相互影響。GE 等[5]引入了自柔順系數(shù),通過(guò)有限差分法對(duì)一維模型進(jìn)行數(shù)值求解,分析了不同工藝下的彈性流體動(dòng)壓力分布和膜厚分布。ISHIKAWA等[6]采用高速攝像的方法分析游離磨粒線切割中的流體,觀察了漿料在線網(wǎng)上的成膜狀態(tài)以及在鋸縫內(nèi)的流動(dòng)狀態(tài)。NASSAUER 等[7]通過(guò)實(shí)驗(yàn)觀測(cè),發(fā)現(xiàn)隨著鋸絲的運(yùn)動(dòng),鋸切通道內(nèi)會(huì)出現(xiàn)氣泡的產(chǎn)生和潰滅現(xiàn)象。林麟[8]利用計(jì)算流體力學(xué)(computational fluid dy-namics,CFD)數(shù)值模擬,系統(tǒng)探究了切削液施加于鋸絲及鋸縫 2 種不同供液方式下,線鋸鋸切加工過(guò)程中的流體行為。除此之外,鄭忠利[9]以增加進(jìn)入鋸切區(qū)域有效切削液為基本思想,采用靜電噴霧供液方式代替?zhèn)鹘y(tǒng)的澆注式供液,發(fā)現(xiàn)靜電噴霧供液方式獲得的硅片質(zhì)量比澆注式供液獲得的硅片質(zhì)量好。QIU 等[10]采用將硅錠浸入水箱中(水浴鋸)的方法鋸切硅片,并與傳統(tǒng)澆注式供液方法在相同鋸切參數(shù)下比較,發(fā)現(xiàn)改進(jìn)的金剛石線鋸冷卻潤(rùn)滑方法能獲得表面質(zhì)量更好的硅片。施郁虎[11]將微量潤(rùn)滑技術(shù)應(yīng)用于金剛石線鋸鋸切加工中,發(fā)現(xiàn)硅片表面粗糙度下降,加工效率提高。此外,寧培桓等[12]研究了線鋸鋸切加工過(guò)程中切削液的掛線能力與其黏度和表面張力的關(guān)系,發(fā)現(xiàn)具有合適表面張力及黏度的切削液在線鋸鋸切加工過(guò)程中有助于提高切片效率,切出表面質(zhì)量更好的切片??梢?,在游離與固結(jié)磨料線鋸鋸切加工過(guò)程中,切削液在鋸縫內(nèi)的流動(dòng)狀態(tài)與鋸切效果密切相關(guān)。隨著鋸絲細(xì)線化以及工件大尺寸化,實(shí)際鋸切加工中將出現(xiàn)更小的流體尺度,但當(dāng)下對(duì)于細(xì)徑金剛石線鋸鋸切加工過(guò)程中,切削液充分進(jìn)入鋸縫以發(fā)揮潤(rùn)滑冷卻作用及其在鋸縫內(nèi)的分布狀態(tài)的研究較少。因此,在金剛石線鋸鋸切加工過(guò)程中,對(duì)切削液在小尺寸鋸縫內(nèi)流場(chǎng)的研究具有重要意義。本研究中通過(guò)建立 CFD 鋸縫模型,仿真分析了金剛石線鋸鋸切過(guò)程中切削液在小尺寸鋸縫中的流場(chǎng),為研究金剛石線鋸鋸切加工過(guò)程中切削液的冷卻潤(rùn)滑效果提供了理論基礎(chǔ)。

    1

    鋸縫流場(chǎng)仿真幾何模型

    基于切削液噴射至鋸絲且由鋸絲運(yùn)動(dòng)攜帶切削液進(jìn)入鋸縫內(nèi)的供液方式建立鋸縫流場(chǎng)仿真幾何模型,如圖 1 所示。

    圖 1 所示為帶有單根鋸絲的鋸縫流場(chǎng)仿真幾何模型。金剛石線鋸平推式鋸切加工工件時(shí),鋸絲在加工中形成的弓角很小,為簡(jiǎn)化分析往往忽略線弓的影響,因此將狹長(zhǎng)的鋸絲縫隙簡(jiǎn)化為平直縫隙[13] 。鋸縫寬度、鋸絲線徑的設(shè)定,參考使用芯線直徑為 36 μm 的電鍍金剛石線鋸鋸切 210 mm × 210 mm 單晶硅片時(shí)的參數(shù),鋸絲與工件間的接觸高度參考鋸絲表面磨粒的平均出刃高度,鋸縫流場(chǎng)仿真幾何模型參數(shù)如表 1 所示。

    將鋸縫流場(chǎng)仿真幾何模型導(dǎo)入 ICEM-CFD 軟件,劃分為結(jié)構(gòu)化六面體網(wǎng)格,如圖 2 所示。通過(guò)在不同網(wǎng)格數(shù)量下仿真計(jì)算金剛石線鋸鋸切過(guò)程中切削液在鋸絲與工件接觸一側(cè)中間長(zhǎng)度位置的壓力,對(duì)網(wǎng)格進(jìn)行獨(dú)立性驗(yàn)證,結(jié)果如表 2 所示。

    由不同網(wǎng)格數(shù)量下的計(jì)算結(jié)果可知:隨著網(wǎng)格數(shù)量的增多,切削液在鋸絲與工件接觸一側(cè)中間長(zhǎng)度位置的壓力趨于穩(wěn)定,以網(wǎng)格數(shù)量 3 518 680 條件下的計(jì)算結(jié)果為基準(zhǔn),計(jì)算得到網(wǎng)格數(shù)量 5 581 660 條件下的計(jì)算結(jié)果誤差僅為 1.67%,因此為縮短計(jì)算時(shí)間,后續(xù)計(jì)算時(shí)采用的網(wǎng)格數(shù)量為 3 518 680。

    5

    仿真結(jié)果

    對(duì)于各方程的求解方法,動(dòng)量方程由基于壓力的求解器求解,通過(guò)隱式時(shí)間離散化求解連續(xù)性方程,對(duì)于壓力速度耦合,使用 PISO 算法,得到的結(jié)果如下。

    5.1

    走絲速度的影響

    隨著工件大尺寸化及鋸絲細(xì)線化的發(fā)展,鋸縫尺寸越來(lái)越小,因此需要研究小尺寸鋸縫條件下如何使切削液充分進(jìn)入鋸切區(qū)域。進(jìn)入鋸縫的切削液流體主要是剪切流,在剪切流體中,影響流體運(yùn)動(dòng)狀態(tài)的因素主要為走絲速度,因此在小尺寸鋸縫(D w =60 μm)以及切削液入口供液量 Q=2.25 × 10?5 L 的條件下,對(duì)走絲速度的影響進(jìn)行了分析計(jì)算。

    因鋸縫狹長(zhǎng),為能更好地對(duì)切削液在鋸縫內(nèi)的分布情況進(jìn)行觀察,在后續(xù)計(jì)算獲取切削液在鋸縫內(nèi)體積分?jǐn)?shù)分布云圖時(shí),重點(diǎn)觀察了鋸絲軸向方向平面上的鋸絲進(jìn)線端、鋸縫中部及鋸絲出線端的切削液體積分?jǐn)?shù)分布云圖。圖 3 為所觀察鋸縫位置的示意圖,圖 4為不同走絲速度條件下切削液在鋸絲與工件接觸一側(cè)及非接觸一側(cè)的體積分?jǐn)?shù)分布云圖,并且提取了最終進(jìn)入鋸縫內(nèi)部的切削液量,圖中用 V 表示。

    鋸縫內(nèi)鋸絲與工件的接觸區(qū)域一側(cè)以及非接觸區(qū)域一側(cè)均未完全充滿液體,且在接觸區(qū)域液相體積分?jǐn)?shù)<100%,區(qū)域內(nèi)存在空氣層,這將對(duì)切片質(zhì)量產(chǎn)生不利影響,該結(jié)果也與文獻(xiàn) [7] 中實(shí)驗(yàn)觀測(cè)到的現(xiàn)象相吻合。隨著走絲速度的增大,進(jìn)入鋸縫的切削液也增多,當(dāng)v w gt;25 m/s 時(shí),接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)均充滿液體,切削液的充分進(jìn)入將會(huì)提升切片質(zhì)量。

    單晶硅等高硬脆材料的斷裂韌性越大,單位截面積所受的鋸切力也越大,因此降低鋸切力可以有效改善切片的表面質(zhì)量[17] 。金剛石線鋸鋸切加工過(guò)程中,切削液附著在鋸絲上并形成液體膜,由鋸絲運(yùn)動(dòng)帶入鋸縫,鋸縫內(nèi)部的切削液壓力會(huì)與鋸切力部分抵消,因此,為探究鋸縫中切削液的壓力分布,在鋸絲與工件接觸側(cè)和非接觸側(cè)進(jìn)行了壓力分布捕捉。如圖 5 所示,切削液的壓力將沿著線 1 及線 2 測(cè)量。

    圖 6 及圖 7 為不同走絲速度條件下,切削液在鋸絲與工件接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)的壓力分布曲線圖。

    由圖 6 及圖 7 可知:鋸縫內(nèi)切削液壓力在鋸絲與工件的接觸區(qū)域一側(cè)和非接觸區(qū)域一側(cè)隨著走絲速度的增大總體上呈增大趨勢(shì),兩側(cè)的壓力差也總體上呈增大趨勢(shì),兩側(cè)壓力分布更加平穩(wěn),這將更有利于提高金剛石線鋸鋸切加工后切片的表面質(zhì)量。鋸縫兩側(cè)均充滿液體后,接觸區(qū)域一側(cè)壓力在 0.179 0 MPa 左右,非接觸區(qū)域一側(cè)壓力在 0.159 0 MPa 左右。接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)的切削液壓力分布與液體在鋸縫中的分布狀態(tài)密切相關(guān),當(dāng)切削液在鋸切區(qū)域分布不均時(shí),鋸絲附近出現(xiàn)空氣層,這使得氣體存在區(qū)域與液體存在區(qū)域的壓力存在差值,從而引起壓力波動(dòng),對(duì)鋸切加工后的切片表面質(zhì)量產(chǎn)生不利影響。

    5.2

    切削液物理屬性的影響

    切削液的物理屬性與其在鋸縫內(nèi)的流動(dòng)狀態(tài)密切相關(guān)。為了研究切削液物理屬性對(duì)金剛石線鋸鋸切加工過(guò)程中鋸縫內(nèi)切削液流場(chǎng)的影響,進(jìn)行了 5 種切削液(見表 3)在 v w =10 m/s 條件下的物理性質(zhì)仿真研究。

    圖 8 為不同切削液物理屬性條件下,切削液在鋸絲與工件接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)的體積分?jǐn)?shù)分布云圖。

    由圖 8 可知:隨著切削液黏度及表面張力的逐漸減小,切削液逐漸充滿鋸縫,但當(dāng)切削液的黏度及表面張力過(guò)小時(shí),進(jìn)入鋸縫的切削液會(huì)有減少的趨勢(shì)。文獻(xiàn) [12] 發(fā)現(xiàn),切削液黏度及表面張力過(guò)大或過(guò)小都不利于切削液的掛線性能,合適的黏度及表面張力更有利于切削液黏附在鋸絲上。在金剛石線鋸鋸切加工過(guò)程中,切削液通過(guò)黏附在鋸絲上,由鋸絲運(yùn)動(dòng)攜帶進(jìn)入鋸縫,所以進(jìn)入鋸縫的切削液量越多,說(shuō)明切削液黏附在鋸絲上的量越多,因此仿真結(jié)果與文獻(xiàn) [12] 實(shí)驗(yàn)得到的不同黏度及表面張力液體掛線性能的規(guī)律相吻合,由此說(shuō)明本次仿真研究結(jié)果的正確性。

    圖 9 及圖 10 為不同切削液物理屬性條件下,切削液在鋸絲與工件接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)的壓力分布曲線圖。由圖 9 和圖 10 可知:鋸縫內(nèi)切削液壓力在鋸絲與工件的接觸區(qū)域一側(cè)和非接觸區(qū)域一側(cè)隨著切削液黏度及表面張力的減小總體呈增大趨勢(shì),且壓力分布更加平穩(wěn)。但黏度及表面張力過(guò)小時(shí),壓力也會(huì)產(chǎn)生一定的波動(dòng),這是因?yàn)榇藭r(shí)切削液不能很好地黏附在鋸絲上,進(jìn)入鋸縫內(nèi)的切削液減少,進(jìn)而會(huì)影響切片質(zhì)量。接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)的壓力分布與鋸縫內(nèi)切削液分布密切相關(guān),對(duì)比分析 5 種切削液,在切削液 C 的物理屬性下,接觸區(qū)域一側(cè)及非接觸區(qū)域一側(cè)的切削液壓力分布更為平穩(wěn),接觸區(qū)域一側(cè)的最大壓力在 0.016 9 MPa 左右,非接觸區(qū)域一側(cè)的最大壓力在 0.013 2 MPa 左右。

    6

    結(jié)"論

    (1)在小尺寸鋸縫下,走絲速度較小時(shí),切削液難以充分進(jìn)入鋸切區(qū)域發(fā)揮作用,隨著走絲速度的增大(v w gt;25 m/s),鋸絲與工件接觸區(qū)域及非接觸區(qū)域內(nèi)逐漸充滿液體,接觸區(qū)域與非接觸區(qū)域的切削液壓力及壓力差整體呈增大趨勢(shì),當(dāng)接觸區(qū)域及非接觸區(qū)域內(nèi)均充滿切削液時(shí),鋸縫內(nèi)切削液壓力分布較為穩(wěn)定,接觸區(qū)域的壓力在 0.179 0 MPa 左右,非接觸區(qū)域的壓力在 0.159 0 MPa 左右。

    (2)液體黏度和表面張力在一定范圍內(nèi)的降低,有利于保證鋸縫內(nèi)切削液的相對(duì)飽和與穩(wěn)定,同時(shí)可以使鋸縫內(nèi)切削液壓力分布更為穩(wěn)定。綜合對(duì)比 5 組切削液的物理屬性,金剛石線鋸鋸切加工過(guò)程中切削液 C 的物理屬性(密度為872.5 kg?m?3 ,黏度為1.15 mPa?s,表面張力為 34.02 mN?m?1 ,與單晶硅表面接觸角為 20 °,與鋸絲表面接觸角為 33 °)更有利于其進(jìn)入鋸縫。

    參考文獻(xiàn):

    [1]葛培琪, 陳自彬, 王沛志. 單晶硅切片加工技術(shù)研究進(jìn)展 [J]. 金剛石與磨料磨具工程,2020,40(4):12-18.

    GE Peiqi, CHEN Zibin, WANG Peizhi. Review of monocrystalline

    silicon slicing technology [J]. Diamond amp; Abrasives Engineering,2020,40(4):12-18.

    [2]BHAGAVAT S, KAO I. A finite element analysis of temperature

    variation in silicon wafers during wiresaw slicing [J]. International

    Journal of Machine Tools and Manufacture,2008,48(1):95-106.

    [3]M?LLER H J. Basic mechanisms and models of multi-wire sawing [J].

    Advanced Engineering Materials,2004,6(7):501-513.

    [4]ZHU L Q, KAO I. Galerkin-based modal analysis on the vibration of

    wire–slurry system in wafer slicing using a wiresaw [J]. Journal of Sound

    and Vibration,2005,283(3/4/5):589-620.

    [5]GE P Q, SANG B, GAO Y F. The numerical analysis on action

    mechanism of slurry in free abrasive wiresaw slicing [J]. Key

    Engineering Materials,2007(359/360):455-459.

    [6]ISHIKAWA K I, SUWABE H, ITOH S I, et al. A basic study of the

    behavior of slurry action at multi-wire saw [J]. Key Engineering

    Materials,2003(238/239):89-92.

    [7]NASSAUER B, HESS A, KUNA M. Numerical and experimental

    investigations of micromechanical processes during wire sawing [J].

    International Journal of Solids and Structures,2014,51(14):2656-2665.

    [8]林麟. 磨粒線切割中的流體行為研究 [D]. 廈門: 華僑大學(xué), 2020.

    LIN Lin. Study on fluid behavior in abrasive wire cutting [D]. Xiamen:

    Huaqiao University, 2020.

    [9]鄭忠利. 靜電噴霧金剛石線鋸切割液吸附特性及切割實(shí)驗(yàn)研究 [D].杭州: 浙江工業(yè)大學(xué), 2020.

    ZHENG Zhongli. Adsorption characteristics and cutting experiments of

    electrostatic spray diamond wire saw cutting fluid [D]. Hangzhou:

    Zhejiang University of Technology, 2020.

    [10]QIU J, LI X F, ZHANG S B. Research on an improved bath cooling and

    lubrication method for diamond wire sawing [J]. The International

    Journal of Advanced Manufacturing Technology,2021,112:1-10.

    [11]施郁虎. 基于金剛石線鋸的微量潤(rùn)滑與工藝研究 [D]. 鎮(zhèn)江: 江蘇科技大學(xué), 2020.

    SHI Yuhu. Micro lubrication and process research based on diamond

    wire saw [D]. Zhenjiang: Jiangsu University of Science and Technology,2020.

    [12]寧培桓, 周建偉, 劉玉嶺, 等. Si 單晶片切削液掛線性能的研究 [J]. 半導(dǎo)體技術(shù),2008,33(11):981-984.

    NING Peihuan, ZHOU Jianwei, LIU Yuling, et al. Study on the

    capability of adhesion of monocrystalline Si wafer cutting fluid [J].

    Semiconductor Technology,2008,33(11):981-984.

    [13]馮勇, 王曉宇, 徐振欽. 單晶硅料擺輔助多金剛線切片的鋸切力模型研究 [J]. 機(jī)械工程學(xué)報(bào),2021,57(19):260-272.

    FENG Yong, WANG Xiaoyu, XU Zhenqin. Research on cutting force

    model of diamond wire saw with single crystal silicon swing [J]. Journal

    of Mechanical Engineering,2021,57(19):260-272.

    [14]章梓雄, 董曾南. 黏性流體力學(xué) [M]. 北京: 清華大學(xué)出版社, 2011.

    ZHANG Zixiong, DONG Zengnan. Viscous fluid mechanics [M].

    Beijing: Tsinghua University Press, 2011.

    [15]任玉新, 陳海昕. 計(jì)算流體力學(xué)基礎(chǔ) [M]. 北京: 清華大學(xué)出版社, 2006.

    REN Yuxin, CHEN Haixin. Fundamentals of computational fluid

    dynamics [M]. Beijing: Tsinghua University Press, 2006.

    [16]L·普朗特. 流體力學(xué)概論 [M]. 郭永懷, 陸士嘉, 譯. 北京: 科學(xué)出版社,1984.

    L PRANDTL. Introduction to fluid mechanics [M]. GUO Yonghuai, LU

    Shijia, translated. Beijing: Science Press, 1984.

    [17]雷娟棉, 譚朝明. 基于 Transition SST 模型的高雷諾數(shù)圓柱繞流數(shù)值研究 [J]. 北京航空航天大學(xué)學(xué)報(bào),2017,43(2):207-217.

    LEI Juanmian, TAN Zhaoming. Numerical simulation for flow around

    circular cylinder at high Reynolds number based on Transition SST

    model [J]. Journal of Beijing University of Aeronautics and Astronautics,2017,43(2):207-217.

    作者簡(jiǎn)介

    通信作者: 葛培琪,男,1963 年生,教授。主要研究方向:金剛石線鋸技術(shù)。

    E-mail: pqge@sdu.edu.cn

    (編輯:李利娟)

    Simulation"study"of"cutting"fluid"flow"field"in"kerf"of"fine"diameter

    diamond"wire"saw

    CHEN Jiahu

    1 , GE Peiqi 1,2

    (1. School of Mechanical Engineering, Shandong University, Jinan 250061, China)

    (2. Key Laboratory of High Efficiency and Clean Mechanical Manufacture of Ministry of Education,

    Shandong University, Jinan 250061, China)

    Abstract

    Objectives:"Electroplated diamond wire saws are widely used in the field of slicing hard and brittle materi-als such as monocrystalline silicon and sapphire. Cutting fluid should give full play to its role in the sawing process,which is conducive to the improvement of wafer quality. As the size of the wafer increases and the diameter of the wiresaw decreases, the kerf in the sawing process becomes deeper and narrower, and the cutting fluid cannot enter the kerf inlarge quantities, resulting in worse lubrication and cooling effects during the sawing process, which leads to the declineof the surface quality of the wafer. Based on computational fluid dynamics (CFD) numerical simulation methods, thecutting fluid flow field in the cutting seam of the diamond wire saw was analyzed and studied. Methods:"In this paper,based on CFD numerical simulation methods, the cutting fluid flow field in the sawing seam of the diamond wire saw isanalyzed and studied. Firstly, according to the actual situation of the diamond wire saw cutting process, a 3D simulationgeometric model is established based on the liquid supply mode, where cutting fluid flows along the sawing wire and isbrought into the sawing area by the motion of the sawing wire. Heat transfer is not considered in this study, and the cut-ting fluid is assumed to be a viscous incompressible fluid. The governing equations of fluid flow include the continuityequation and the momentum equation. It is found from the equations that the main factors affecting the flow field distri-bution of cutting fluid in the kerf are wire speed and cutting fluid density. By calculating Reynolds number and Webernumber, the fluid model studied in this paper is selected as the Transition SST model. The VOF method is determined tocharacterize the fluid state of cutting fluid in the saw joint, and the CSF model is introduced into the VOF method tocharacterize the influence of surface tension. Considering the influence of physical properties of cutting fluid, the dens-ity, viscosity, surface tension, and wall contact angle of cutting fluid are measured experimentally. The momentumequation is solved by a pressure-based solver, the continuity equation is solved by implicit time discretization, and thePISO method is used for pressure-velocity coupling. Result: With the increase of chip size and the decrease of wire sawdiameter, the size of the saw seam is getting smaller and smaller. The main fluid entering the saw seam is shear flow,and the main factor affecting the fluid motion state is the wire speed. Under the condition of small-size sawing, when thewire speed is low (v w ≤25 m/s), both the contact area and the non-contact area of the saw wire in the sawing joint are notcompletely filled with liquid, and the liquid volume fraction in the contact area is lt; 100%, with an air layer in the area.With the increase of wire speed, more and more cutting fluid enters the sawing joint. With the increase of wire speed,more and more cutting fluid enters the saw seam. When v w gt;25 m/s, both the contact area and the non-contact area arefilled with liquid. The cutting fluid pressure in the saw seam increases with the increase of wire speed on both the con-tact area side and the non-contact area side, and the pressure difference on both sides also increases generally. The pres-sure distribution on both sides becomes more stable. After both sides of the saw seam are filled with liquid, the pressureon one side of the contact area is about 0.179 0 MPa, and the pressure on one side of the non-contact area is about 0.159 0MPa. With the gradual reduction of the viscosity and surface tension of the cutting fluid, the cutting fluid gradually fillsthe saw joint. Howver, when the viscosity and surface tension of the cutting fluid are too small, the cutting fluid enter- ing the saw joint will tend to decrease. The cutting fluid pressure in the saw seam increases with the decrease of the cut-ting fluid viscosity and surface tension in both the contact area and the non-contact area, and the pressure distributionbecomes more stable. However, when the viscosity and the surface tension are too small, the pressure will also fluctuate.Under the physical properties of the cutting fluid C, with a density of 872.5 kg?m?3 , viscosity of 1.15 mPa?s and surfacetension of 34.02 mN?m?1 , the pressure distribution of the cutting fluid on the contact area side and the non-contact areaside is more stable. The maximum pressure on the contact area side is about 0.016 9 MPa, and the maximum pressure onthe non-contact area side is about 0.013 2 MPa. Conclusions: (1) Under small sawing sizes, when the wire speed is low,the cutting fluid is difficult to fully enter the sawing area to play its role. With the increase of the wire speed (v w gt;25m/s), the contact area and non-contact area between the saw wire and the workpiece are gradually filled with liquid, andthe cutting fluid pressure and pressure difference between the contact area and the non-contact area show an overall in-creasing trend. When the contact area and non-contact area are full of cutting fluid, the cutting fluid pressure distribu-tion in the saw joint is relatively stable, with the pressure in the contact area being about 0.179 0 MPa and the pressurein the non-contact area being about 0.159 0 MPa. (2) The reduction of liquid viscosity and surface tension within a cer-tain range is conducive to ensuring the relative saturation and stability of the cutting fluid in the saw joint, and at thesame time, it can make the pressure distribution of the cutting fluid in the saw joint more stable. A comprehensive com-parison of the physical properties of the 5 groups of cutting fluids shows that the physical properties of the cutting fluidC during the diamond line saw cutting process are more conducive to its entry into the saw joint.

    Key"words

    diamond wire saw;computational fluid dynamics;small size kerf;cutting fluid flow field

    猜你喜歡
    計(jì)算流體力學(xué)
    熱等離子體煤裂解制乙炔的轉(zhuǎn)化規(guī)律與定向控制
    《工程流體力學(xué)》教學(xué)方法探討
    基于預(yù)條件技術(shù)的風(fēng)力機(jī)葉片計(jì)算方法研究
    風(fēng)力發(fā)電機(jī)組機(jī)艙結(jié)構(gòu)散熱性能分析方法
    汽車前端冷凝器進(jìn)風(fēng)不均勻性分析
    ORVR系統(tǒng)加油管液封性能數(shù)值模擬
    考慮懸移質(zhì)效應(yīng)的橋墩動(dòng)床沖刷精細(xì)化分析方法
    CFD數(shù)值模擬在計(jì)算流體力學(xué)項(xiàng)目驅(qū)動(dòng)實(shí)踐教學(xué)中的應(yīng)用
    CFD軟件及源程序輔助工程流體力學(xué)課程教學(xué)的探索與實(shí)踐
    民用飛機(jī)靜壓源位置誤差修正設(shè)計(jì)研究①
    科技資訊(2015年17期)2015-10-09 21:02:59
    18禁裸乳无遮挡动漫免费视频| 久久久国产精品麻豆| 在线观看免费视频网站a站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品国产精品| 婷婷色av中文字幕| 最近的中文字幕免费完整| 国产老妇伦熟女老妇高清| 天天躁日日躁夜夜躁夜夜| 99久久综合免费| 9191精品国产免费久久| 91在线精品国自产拍蜜月| 一区二区三区四区激情视频| 久久影院123| 国产亚洲欧美精品永久| 亚洲av男天堂| 亚洲欧洲国产日韩| 日日撸夜夜添| 最近最新中文字幕免费大全7| 国产精品无大码| 免费人妻精品一区二区三区视频| 欧美日韩成人在线一区二区| 男的添女的下面高潮视频| 交换朋友夫妻互换小说| 日韩av在线免费看完整版不卡| 天天躁夜夜躁狠狠久久av| 精品国产乱码久久久久久男人| 午夜福利在线观看免费完整高清在| 永久免费av网站大全| 美女大奶头黄色视频| 久久99一区二区三区| 叶爱在线成人免费视频播放| av有码第一页| 久久久国产精品麻豆| 免费女性裸体啪啪无遮挡网站| 寂寞人妻少妇视频99o| 国产成人精品无人区| 电影成人av| 日产精品乱码卡一卡2卡三| 只有这里有精品99| 五月伊人婷婷丁香| 久久久久久伊人网av| 成人国语在线视频| 一级毛片电影观看| 亚洲国产欧美日韩在线播放| 午夜福利影视在线免费观看| 久久久精品区二区三区| 黄色一级大片看看| av片东京热男人的天堂| 亚洲四区av| 亚洲少妇的诱惑av| 亚洲成人av在线免费| 亚洲综合色惰| 欧美人与性动交α欧美软件| 日韩一区二区视频免费看| 欧美精品人与动牲交sv欧美| 精品午夜福利在线看| 午夜福利视频精品| 一个人免费看片子| 日韩一本色道免费dvd| 国产精品免费大片| 亚洲一区中文字幕在线| 久久精品熟女亚洲av麻豆精品| 永久网站在线| 国产亚洲最大av| 97精品久久久久久久久久精品| 80岁老熟妇乱子伦牲交| 亚洲人成电影观看| 国精品久久久久久国模美| 如何舔出高潮| 熟女电影av网| 色哟哟·www| 下体分泌物呈黄色| 1024香蕉在线观看| 中文字幕人妻熟女乱码| 好男人视频免费观看在线| 亚洲图色成人| 在线看a的网站| 亚洲国产欧美在线一区| 青春草国产在线视频| 日韩在线高清观看一区二区三区| 国产精品一国产av| av国产精品久久久久影院| 香蕉精品网在线| 日韩熟女老妇一区二区性免费视频| 久久久久精品久久久久真实原创| 亚洲精品国产av成人精品| 欧美av亚洲av综合av国产av | 免费高清在线观看视频在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美日韩在线播放| 精品国产露脸久久av麻豆| 色吧在线观看| 亚洲欧美中文字幕日韩二区| 中文字幕人妻丝袜一区二区 | 亚洲 欧美一区二区三区| 黄色 视频免费看| www.av在线官网国产| 黄片无遮挡物在线观看| 一二三四中文在线观看免费高清| 国产精品久久久av美女十八| 国产精品 国内视频| 国产片内射在线| 国产乱来视频区| 欧美+日韩+精品| 国产欧美日韩综合在线一区二区| 卡戴珊不雅视频在线播放| 一二三四在线观看免费中文在| 亚洲伊人色综图| 天天躁夜夜躁狠狠躁躁| 少妇被粗大猛烈的视频| 中文字幕人妻熟女乱码| 日韩av免费高清视频| 一级毛片电影观看| 一区二区三区激情视频| 久久这里有精品视频免费| 老司机影院毛片| 国产有黄有色有爽视频| 赤兔流量卡办理| 麻豆乱淫一区二区| 丝瓜视频免费看黄片| 久久精品国产亚洲av天美| 飞空精品影院首页| 2021少妇久久久久久久久久久| 一级毛片电影观看| 日本色播在线视频| 亚洲精品美女久久久久99蜜臀 | 日韩中文字幕欧美一区二区 | 边亲边吃奶的免费视频| 丝袜喷水一区| 国产亚洲一区二区精品| 日本爱情动作片www.在线观看| 丰满迷人的少妇在线观看| 日日撸夜夜添| 亚洲一区二区三区欧美精品| 五月开心婷婷网| 少妇精品久久久久久久| 丝瓜视频免费看黄片| 久久久国产欧美日韩av| 午夜精品国产一区二区电影| 婷婷色av中文字幕| 丰满饥渴人妻一区二区三| 国产1区2区3区精品| 亚洲精品乱久久久久久| 精品一区二区三区四区五区乱码 | 另类精品久久| 春色校园在线视频观看| 免费不卡的大黄色大毛片视频在线观看| 国产成人精品婷婷| 啦啦啦在线免费观看视频4| 卡戴珊不雅视频在线播放| 啦啦啦啦在线视频资源| 一区二区日韩欧美中文字幕| 久久99一区二区三区| 亚洲色图 男人天堂 中文字幕| 国产亚洲欧美精品永久| 久久 成人 亚洲| 秋霞在线观看毛片| 欧美97在线视频| 黑丝袜美女国产一区| 久久99精品国语久久久| 国产精品免费大片| 人妻人人澡人人爽人人| 亚洲国产成人一精品久久久| 国产成人午夜福利电影在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲av.av天堂| 一区二区三区四区激情视频| 性色av一级| 一区福利在线观看| 天堂中文最新版在线下载| 美女主播在线视频| 啦啦啦中文免费视频观看日本| 九九爱精品视频在线观看| 天堂中文最新版在线下载| 亚洲国产精品一区三区| 97在线人人人人妻| 免费高清在线观看视频在线观看| 免费日韩欧美在线观看| 黄色配什么色好看| 欧美中文综合在线视频| 国产精品久久久久久久久免| 欧美xxⅹ黑人| 91精品伊人久久大香线蕉| 一区福利在线观看| 一级,二级,三级黄色视频| 一区二区av电影网| 在线精品无人区一区二区三| 亚洲av中文av极速乱| 午夜激情久久久久久久| 欧美日韩国产mv在线观看视频| 久久久久久人人人人人| 1024香蕉在线观看| 母亲3免费完整高清在线观看 | 99久久综合免费| 国产爽快片一区二区三区| av有码第一页| 人人妻人人澡人人看| 两个人看的免费小视频| 欧美日韩精品网址| 少妇精品久久久久久久| 亚洲成色77777| 99热国产这里只有精品6| 久久影院123| 超碰97精品在线观看| 香蕉丝袜av| 国产人伦9x9x在线观看 | 欧美另类一区| 国产一区有黄有色的免费视频| 国产精品无大码| 不卡av一区二区三区| 久久久久精品人妻al黑| √禁漫天堂资源中文www| 免费黄网站久久成人精品| 国产有黄有色有爽视频| 性少妇av在线| 春色校园在线视频观看| 午夜免费观看性视频| 如日韩欧美国产精品一区二区三区| 91精品伊人久久大香线蕉| 国产乱人偷精品视频| 丝袜美腿诱惑在线| 少妇被粗大猛烈的视频| 在线观看免费高清a一片| 国产色婷婷99| 只有这里有精品99| 99久久综合免费| 亚洲国产看品久久| 国产福利在线免费观看视频| 成人漫画全彩无遮挡| 丰满饥渴人妻一区二区三| 亚洲中文av在线| 一本—道久久a久久精品蜜桃钙片| 久久国内精品自在自线图片| 久久毛片免费看一区二区三区| 午夜日本视频在线| av国产精品久久久久影院| 老熟女久久久| 男女下面插进去视频免费观看| 成人漫画全彩无遮挡| 伊人久久大香线蕉亚洲五| 国产成人精品久久久久久| 97在线视频观看| a 毛片基地| 赤兔流量卡办理| 2022亚洲国产成人精品| 久久人人爽av亚洲精品天堂| 久久精品久久久久久久性| 亚洲精品美女久久av网站| 亚洲av成人精品一二三区| 日韩免费高清中文字幕av| 国产精品久久久av美女十八| 寂寞人妻少妇视频99o| 日韩av免费高清视频| 日韩 亚洲 欧美在线| 久久国产精品男人的天堂亚洲| 你懂的网址亚洲精品在线观看| 久久久国产一区二区| 亚洲欧美精品综合一区二区三区 | 精品国产乱码久久久久久小说| 午夜日韩欧美国产| 免费看不卡的av| 国产成人免费观看mmmm| 丝袜在线中文字幕| 黄片无遮挡物在线观看| 国产精品免费视频内射| 成人毛片60女人毛片免费| 美女福利国产在线| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 久久这里有精品视频免费| 久久久久精品久久久久真实原创| 国产精品熟女久久久久浪| 热99久久久久精品小说推荐| 午夜av观看不卡| 亚洲四区av| 中文字幕人妻丝袜一区二区 | 久久影院123| 韩国精品一区二区三区| 婷婷色麻豆天堂久久| 精品国产国语对白av| 少妇的逼水好多| 一区二区av电影网| 9热在线视频观看99| 亚洲精品美女久久久久99蜜臀 | 看十八女毛片水多多多| 亚洲av.av天堂| 亚洲国产精品成人久久小说| 亚洲,欧美精品.| 日韩精品有码人妻一区| 啦啦啦在线观看免费高清www| 在线免费观看不下载黄p国产| 免费在线观看完整版高清| 国产国语露脸激情在线看| 免费日韩欧美在线观看| 尾随美女入室| 日本vs欧美在线观看视频| 久久久久精品性色| 99re6热这里在线精品视频| 亚洲国产av新网站| 久久久久精品性色| 黄色配什么色好看| 久久久国产一区二区| 婷婷色av中文字幕| 精品国产露脸久久av麻豆| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| 午夜久久久在线观看| 丁香六月天网| av免费在线看不卡| 国产日韩欧美亚洲二区| 免费在线观看黄色视频的| 亚洲精品在线美女| 亚洲美女黄色视频免费看| 在线免费观看不下载黄p国产| 国产成人精品久久二区二区91 | 久久精品夜色国产| 日韩人妻精品一区2区三区| 亚洲欧美清纯卡通| 国产av码专区亚洲av| 久久国产亚洲av麻豆专区| 午夜福利,免费看| 午夜久久久在线观看| 国产又爽黄色视频| 高清av免费在线| 久久久久久久亚洲中文字幕| xxx大片免费视频| 老女人水多毛片| 日韩不卡一区二区三区视频在线| 伊人久久大香线蕉亚洲五| 五月伊人婷婷丁香| 18禁观看日本| 日韩制服骚丝袜av| 国产野战对白在线观看| 一本色道久久久久久精品综合| 精品人妻一区二区三区麻豆| 精品久久久久久电影网| av在线app专区| 国产麻豆69| 久久久久久久大尺度免费视频| 国产一区二区三区av在线| 麻豆精品久久久久久蜜桃| 亚洲av中文av极速乱| 国产精品99久久99久久久不卡 | 熟女少妇亚洲综合色aaa.| 免费黄频网站在线观看国产| 久久久久国产一级毛片高清牌| 日日啪夜夜爽| 一级毛片 在线播放| 亚洲国产精品国产精品| 色吧在线观看| 日本欧美视频一区| 国产成人精品婷婷| 80岁老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 国产不卡av网站在线观看| 91久久精品国产一区二区三区| 日韩电影二区| 欧美中文综合在线视频| 国产又色又爽无遮挡免| 国产精品女同一区二区软件| 亚洲精品美女久久av网站| 看免费av毛片| 在线天堂中文资源库| 亚洲 欧美一区二区三区| 捣出白浆h1v1| av卡一久久| 国产精品一区二区在线观看99| 热re99久久国产66热| 亚洲成色77777| 中文字幕av电影在线播放| 不卡视频在线观看欧美| 日本爱情动作片www.在线观看| 精品人妻熟女毛片av久久网站| 丰满乱子伦码专区| 18禁国产床啪视频网站| 日本爱情动作片www.在线观看| 国产在线视频一区二区| 91精品伊人久久大香线蕉| 欧美激情极品国产一区二区三区| 在线观看免费高清a一片| 国产免费现黄频在线看| 一区二区三区乱码不卡18| 在线天堂最新版资源| 丰满少妇做爰视频| 制服丝袜香蕉在线| 狂野欧美激情性bbbbbb| 欧美日韩成人在线一区二区| 如何舔出高潮| 国产一区二区在线观看av| 亚洲伊人色综图| av免费在线看不卡| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 日韩在线高清观看一区二区三区| 黄色怎么调成土黄色| 日本黄色日本黄色录像| 99久久综合免费| 国产精品av久久久久免费| 欧美 亚洲 国产 日韩一| 午夜福利网站1000一区二区三区| 欧美精品亚洲一区二区| 日产精品乱码卡一卡2卡三| 久久久久久免费高清国产稀缺| av免费观看日本| 在线观看www视频免费| 免费观看a级毛片全部| 18禁观看日本| 亚洲一级一片aⅴ在线观看| 欧美在线黄色| 亚洲精品国产色婷婷电影| av视频免费观看在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲成av片中文字幕在线观看 | 黄片小视频在线播放| 我的亚洲天堂| 久久久久久久久久久免费av| 国产精品久久久久成人av| 男女免费视频国产| 麻豆乱淫一区二区| 日韩伦理黄色片| 伊人久久国产一区二区| 一本—道久久a久久精品蜜桃钙片| 男人添女人高潮全过程视频| 男女国产视频网站| 亚洲精品美女久久久久99蜜臀 | 青青草视频在线视频观看| 91精品三级在线观看| 国产1区2区3区精品| 99精国产麻豆久久婷婷| 美女国产高潮福利片在线看| 亚洲精品,欧美精品| 蜜桃国产av成人99| 久久久久久伊人网av| 欧美在线黄色| 久久综合国产亚洲精品| 国产精品.久久久| 亚洲精品国产av成人精品| 成人毛片a级毛片在线播放| 麻豆精品久久久久久蜜桃| 欧美日韩精品成人综合77777| 精品亚洲乱码少妇综合久久| 亚洲内射少妇av| 美女中出高潮动态图| 久久久久人妻精品一区果冻| 九草在线视频观看| 午夜老司机福利剧场| 男女下面插进去视频免费观看| 狂野欧美激情性bbbbbb| 美女主播在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费av中文字幕在线| 亚洲,欧美,日韩| 亚洲国产欧美网| 免费观看在线日韩| 久久久久国产精品人妻一区二区| 日本午夜av视频| 精品亚洲成a人片在线观看| 午夜激情久久久久久久| 欧美人与善性xxx| 国产成人免费无遮挡视频| 90打野战视频偷拍视频| www日本在线高清视频| 最新的欧美精品一区二区| 999精品在线视频| 搡女人真爽免费视频火全软件| 国产精品麻豆人妻色哟哟久久| 激情五月婷婷亚洲| 五月开心婷婷网| 亚洲av免费高清在线观看| 午夜日本视频在线| www.自偷自拍.com| 有码 亚洲区| 久久精品国产亚洲av涩爱| 男男h啪啪无遮挡| a级片在线免费高清观看视频| 午夜免费鲁丝| 一级毛片黄色毛片免费观看视频| 亚洲精品第二区| 精品国产乱码久久久久久小说| 少妇精品久久久久久久| 亚洲av电影在线观看一区二区三区| 久久鲁丝午夜福利片| 午夜老司机福利剧场| 国产精品 国内视频| 飞空精品影院首页| xxx大片免费视频| 午夜福利,免费看| 王馨瑶露胸无遮挡在线观看| 天天操日日干夜夜撸| 激情视频va一区二区三区| 日韩制服骚丝袜av| 赤兔流量卡办理| 亚洲,一卡二卡三卡| 成年美女黄网站色视频大全免费| 丰满饥渴人妻一区二区三| 精品久久久久久电影网| 欧美变态另类bdsm刘玥| 热re99久久国产66热| 欧美精品一区二区大全| 国产成人欧美| 777米奇影视久久| 国产毛片在线视频| 一本久久精品| 18在线观看网站| 国产爽快片一区二区三区| av网站在线播放免费| 欧美人与善性xxx| 亚洲精品国产av成人精品| 亚洲精品久久久久久婷婷小说| 久久精品国产亚洲av天美| 夫妻午夜视频| 亚洲经典国产精华液单| 一级a爱视频在线免费观看| 欧美日韩一级在线毛片| 亚洲国产看品久久| 日日爽夜夜爽网站| 欧美精品人与动牲交sv欧美| 国产爽快片一区二区三区| 熟女少妇亚洲综合色aaa.| 日韩,欧美,国产一区二区三区| 亚洲av免费高清在线观看| 香蕉国产在线看| 国产一区有黄有色的免费视频| 国产精品女同一区二区软件| 香蕉丝袜av| 国产亚洲一区二区精品| 一本色道久久久久久精品综合| 亚洲国产精品999| 最近手机中文字幕大全| 欧美精品一区二区免费开放| 国产精品麻豆人妻色哟哟久久| 天堂8中文在线网| 丝袜美腿诱惑在线| 亚洲精品在线美女| 亚洲国产欧美在线一区| 看免费成人av毛片| 人妻少妇偷人精品九色| 亚洲欧美成人综合另类久久久| 一区二区日韩欧美中文字幕| 99久久综合免费| 亚洲av中文av极速乱| 亚洲国产精品国产精品| 亚洲av电影在线观看一区二区三区| 99国产综合亚洲精品| 三级国产精品片| 亚洲成国产人片在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线观看人妻少妇| 大码成人一级视频| 欧美变态另类bdsm刘玥| 校园人妻丝袜中文字幕| 99国产综合亚洲精品| 国产精品嫩草影院av在线观看| 日韩一本色道免费dvd| 久久精品久久久久久久性| 欧美人与性动交α欧美精品济南到 | 精品少妇久久久久久888优播| 国产欧美亚洲国产| 久久国产亚洲av麻豆专区| 精品国产一区二区三区久久久樱花| 99久国产av精品国产电影| 美女大奶头黄色视频| 青春草国产在线视频| av福利片在线| 97在线人人人人妻| 国产成人精品婷婷| 1024香蕉在线观看| 亚洲,一卡二卡三卡| 精品国产一区二区三区久久久樱花| 不卡av一区二区三区| 1024视频免费在线观看| 婷婷成人精品国产| 视频区图区小说| 亚洲国产毛片av蜜桃av| 亚洲内射少妇av| 午夜福利,免费看| 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| 夜夜骑夜夜射夜夜干| 欧美精品一区二区大全| 韩国高清视频一区二区三区| 精品国产乱码久久久久久小说| 国产野战对白在线观看| 熟女少妇亚洲综合色aaa.| 亚洲经典国产精华液单| 建设人人有责人人尽责人人享有的| 最近最新中文字幕免费大全7| 久久久久精品人妻al黑| 九色亚洲精品在线播放| 99久久中文字幕三级久久日本| 国产探花极品一区二区| a级毛片黄视频| 精品一区二区免费观看| 中文字幕人妻丝袜一区二区 | 欧美变态另类bdsm刘玥| 一本久久精品| 日韩伦理黄色片| 18禁国产床啪视频网站| 老司机影院成人| 精品一区二区免费观看| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 一级毛片 在线播放| 看免费成人av毛片| 精品少妇一区二区三区视频日本电影 | 亚洲色图 男人天堂 中文字幕| 在线观看免费视频网站a站| 精品久久久精品久久久| 五月伊人婷婷丁香| 国产不卡av网站在线观看| 老司机亚洲免费影院| 人人妻人人澡人人看| 深夜精品福利| 最近2019中文字幕mv第一页|