亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于多物理場耦合的航空發(fā)動機(jī)燃油泵液膜密封仿真分析與試驗(yàn)*

        2022-12-28 05:09:18任寶杰郝木明胡廣陽王銘章蘇志善李勇凡
        潤滑與密封 2022年12期
        關(guān)鍵詞:機(jī)械變形

        任寶杰 郝木明 胡廣陽 王銘章 蘇志善 李勇凡

        (1.東營海森密封技術(shù)有限責(zé)任公司 山東東營 257000;2.中國石油大學(xué)(華東) 新能源學(xué)院 山東青島 266580;3.中國航空發(fā)動機(jī)集團(tuán)航空發(fā)動機(jī)動力傳輸重點(diǎn)實(shí)驗(yàn)室 遼寧沈陽 110015;4.中國航發(fā)西安動力控制科技有限公司 陜西西安 710077)

        機(jī)械密封作為旋轉(zhuǎn)機(jī)械的重要部件,在工業(yè)上得到了廣泛應(yīng)用[1-2]。隨著現(xiàn)代工業(yè)的發(fā)展,機(jī)械密封的使用環(huán)境越來越苛刻,對機(jī)械密封的性能提出了更高的要求。對于接觸式機(jī)械密封,高參數(shù)工況下機(jī)械密封存在潤滑失效問題,導(dǎo)致整機(jī)故障[3]。

        非接觸式液膜密封因其兼具密封性和良好的潤滑特性,能夠?qū)崿F(xiàn)在高參數(shù)工況下的長壽命安全運(yùn)行。在實(shí)際運(yùn)轉(zhuǎn)中,非接觸式液膜密封端面的潤滑流體膜處于微米級,與密封環(huán)受到熱力影響產(chǎn)生的變形處于同一量級,因此在對高參數(shù)液膜密封進(jìn)行分析時,應(yīng)考慮熱-流-固多物理場的耦合效應(yīng)。QIU和KHONSARI[4]基于CFD-ACE+得到了螺旋槽機(jī)械密封副的溫度和壓力分布。WU等[5]采用FLUENT建立密封環(huán)、周圍沖洗液和其他密封組件的三維流固耦合模型,分析了機(jī)械密封的流場和溫度分布規(guī)律。楊丹丹等[6]建立了高溫油泵機(jī)械密封的三維力熱耦合分析模型,指出力熱耦合變形中熱變形占主導(dǎo)地位。朱維兵等[7]采用ANSYS Workbench對斜線槽上游泵送機(jī)械密封進(jìn)行單向流固耦合分析,分析了工況參數(shù)對密封環(huán)變形的影響。

        為提高非接觸式機(jī)械密封的性能,相關(guān)學(xué)者針對端面槽型優(yōu)化開展了大量研究工作。在單列螺旋槽的基礎(chǔ)上,LAI[8]提出了人字槽和Y形槽等雙列螺旋槽結(jié)構(gòu);丁雪興等[9]采用多目標(biāo)優(yōu)化理論對人字槽密封和螺旋槽密封進(jìn)行結(jié)構(gòu)優(yōu)化,結(jié)果表明人字槽密封的綜合性能高于螺旋槽密封;MENG等[10]采用三維熱流模型,分析了人字槽結(jié)構(gòu)參數(shù)對密封性能的影響??傮w而言,與單列螺旋槽液膜密封相比,人字槽液膜密封同時具有上游泵送和下游泵送的能力,具有更大的流體膜厚度及流體膜剛度,穩(wěn)定性優(yōu)于單列螺旋槽[11-12]。黃偉峰等[13-14]發(fā)現(xiàn)了機(jī)械密封在高轉(zhuǎn)速下的溫升和變形對密封性能的削弱效應(yīng),總結(jié)了機(jī)械密封特性參數(shù)變化規(guī)律并應(yīng)用于上游泵送機(jī)械密封的設(shè)計(jì)優(yōu)化。YANG等[15]基于熱彈流分析方法,對三角形織構(gòu)進(jìn)行數(shù)值計(jì)算后發(fā)現(xiàn)該織構(gòu)形槽型對于密封性能具有明顯提升作用。孫鑫暉等[16]在考慮熱流體動力學(xué)特性基礎(chǔ)上,對槽數(shù)、槽深等參數(shù)對于密封性能的影響進(jìn)行計(jì)算分析。

        本文作者以某型航空發(fā)動機(jī)燃油增壓泵燃油側(cè)機(jī)械密封為研究對象,針對燃油側(cè)機(jī)械密封溫升過高導(dǎo)致密封失效的問題進(jìn)行密封端面的結(jié)構(gòu)設(shè)計(jì)?;诹髯儗W(xué)、摩擦學(xué)與潤滑理論、材料學(xué)、傳熱學(xué)等學(xué)科的基本原理,建立了適用于航空發(fā)動機(jī)燃油增壓泵用機(jī)械密封多物理場仿真分析模型,并通過試驗(yàn)驗(yàn)證了人字槽機(jī)械密封的性能滿足燃油泵運(yùn)行要求。

        1 理論模型

        1.1 物理模型

        燃油側(cè)機(jī)械密封動靜環(huán)物理模型如圖1所示。

        圖1 機(jī)械密封結(jié)構(gòu)

        圖2所示為人字槽液膜密封端面結(jié)構(gòu)示意圖,人字槽型端面結(jié)構(gòu)特征參數(shù)為槽數(shù)N與槽深hg,槽數(shù)N預(yù)設(shè)計(jì)范圍為12~24,槽深hg預(yù)設(shè)計(jì)范圍為16.0~40.0 μm。

        圖2 密封端面結(jié)構(gòu)

        1.2 考慮粗糙度的雷諾方程

        考慮表面粗糙度的、描述機(jī)械密封端面液膜壓力分布的控制方程如式(1)所示[13]。

        (1)

        1.3 粗糙表面微凸體彈塑性接觸力模型

        隨微凸體法向變形量逐漸增大,將依次發(fā)生彈性變形、彈塑性變形和塑性變形,同時也產(chǎn)生相應(yīng)的彈性接觸力、彈塑性接觸力和塑性接觸力。描述包括所述3種接觸力階段的理論模型如式(2)所示。

        (2)

        式中:pt(h)、pet(h)、pept(h)、ppt(h)為總接觸壓力、彈性接觸壓力、彈塑性接觸壓力、塑性接觸壓力,Pa;w為微凸體法向變形量;w1、w2為初始屈服點(diǎn)臨界法向變形量、完全塑性變形的臨界法向變形量;ys為表面凸峰平均高度線和表面平均高度線之間的距離;z為微凸體的高度(相對于表面凸峰平均高度線);φs(z)為微凸體高度分布的概率密度函數(shù)(以表面凸峰平均高度線為基準(zhǔn));h為以表面平均高度線為基準(zhǔn)的兩表面平均距離;R為微凸體的曲率半徑;H為軟材料的硬度,Pa;η為微凸體面積密度。

        1.4 能量方程

        描述機(jī)械密封端面液膜溫度分布的控制方程(平均能量方程)如式(3)所示。

        (3)

        1.5 熱傳導(dǎo)方程

        非接觸式機(jī)械密封動環(huán)、靜環(huán)和液膜的溫度計(jì)算方法與接觸式機(jī)械密封動靜環(huán)溫度計(jì)算方法相一致,影響密封環(huán)溫度場的主要因素有2個:一是密封環(huán)旋轉(zhuǎn)生成熱量,二是密封靜動環(huán)與液膜及外界流體進(jìn)行熱量交換。熱量主要由兩部分組成:一是攪拌熱,由于主軸旋轉(zhuǎn)使得流體流動攪拌而產(chǎn)生;二是摩擦熱,因動環(huán)旋轉(zhuǎn)并且液體具有黏性而在端面間隙產(chǎn)生。由于對流傳熱是密封環(huán)同流體熱量交換的主要方式,故而影響環(huán)溫度分布的主要參數(shù)之一便是對流傳熱系數(shù)。

        在整個密封系統(tǒng)中主要涉及到的熱傳遞方式有以下幾種:

        (1)密封動靜環(huán)端面與間隙流體之間是流-固耦合傳熱;

        (2)密封動靜環(huán)內(nèi)部同種介質(zhì)之間通過熱傳導(dǎo)方式傳遞熱量;

        (3)密封環(huán)內(nèi)外側(cè)流體與環(huán)的內(nèi)外徑之間是對流傳熱。

        溫度場和變形的研究對象是動環(huán)和靜環(huán),極坐標(biāo)下彈塑性體的傳熱方程如下:

        (4)

        式中:T為溫度,℃;Q/υ為內(nèi)部生成熱;cp為比熱容,J/(kg·℃);k為傳熱系數(shù),W/(m·℃)。

        傳熱方程變?yōu)?/p>

        (5)

        密封環(huán)的周圍是密封介質(zhì),它們之間的傳熱方式是對流傳熱。對流邊界條件要求垂直于邊界的傳導(dǎo)熱流等于通過對流傳熱帶走的熱量。

        因此對流傳熱邊界滿足如下方程:

        (6)

        對于耦合端面,要求傳遞的總熱量與生成熱相等。則耦合邊界條件數(shù)學(xué)表達(dá)式為

        (7)

        式中:α為對流傳熱系數(shù),W/(m2·℃);T∞為環(huán)境溫度,℃;k1、k2為結(jié)構(gòu)1和2的傳熱系數(shù),W/(m·℃);T1、T2為結(jié)構(gòu)1和2的溫度,℃。

        1.6 對流傳熱系數(shù)

        機(jī)械密封的傳熱計(jì)算主要是指密封動、靜環(huán)與周圍介質(zhì)間傳熱效應(yīng)的計(jì)算,其對流傳熱系數(shù)α目前基本均采用經(jīng)驗(yàn)公式或半經(jīng)驗(yàn)公式進(jìn)行計(jì)算。

        (1)靜環(huán)外圓面與介質(zhì)之間:

        Nu=0.023εlRe0.8Pr0.4

        Nu=2αl/λl

        Re=2vl/υl

        (8)

        式中:εl為修正系數(shù),一般可取εl=2;l為靜環(huán)外圓面與密封箱內(nèi)周壁間的徑向間隙;λl為液體介質(zhì)的導(dǎo)熱系數(shù);υl為液體介質(zhì)的運(yùn)動黏度;v為靜環(huán)外周周圍介質(zhì)的軸向流速。

        (2)動環(huán)外圓面與介質(zhì)之間:

        (9)

        2 數(shù)值求解方法

        采用的仿真計(jì)算流程如圖3所示。該計(jì)算流程以開啟力與閉合力相平衡為最終的收斂判別依據(jù),同時包括液膜壓力場收斂、液膜溫度場收斂、密封環(huán)端面溫度分布收斂等3個子判別依據(jù);滿足上述各判別依據(jù)之后,端面間隙形狀將自行滿足收斂。

        圖3 仿真計(jì)算流程

        表1給出了文中所采用的密封結(jié)構(gòu)參數(shù)和工況參數(shù)。

        表1 密封面幾何結(jié)構(gòu)與工況參數(shù)

        采用數(shù)字式黏度計(jì)測量所購油品在不同溫度下的動力黏度,利用所測數(shù)據(jù)擬合介質(zhì)的黏溫關(guān)系方程(黏溫特性模型),以用于仿真計(jì)算。油品的實(shí)測數(shù)據(jù)及擬合公式如圖4所示。

        圖4 黏溫特性曲線

        表2給出了密封環(huán)溫度場及力熱變形計(jì)算相關(guān)物性參數(shù)。

        表2 密封動靜環(huán)物性參數(shù)

        機(jī)械密封密封環(huán)的力、熱邊界條件序號如圖5所示,密封環(huán)計(jì)算力、熱邊界條件序號所對應(yīng)的邊界類型如表3所示。

        圖5 密封環(huán)邊界序號

        表3 密封環(huán)力、熱邊界條件類型

        3 計(jì)算結(jié)果與分析

        設(shè)定轉(zhuǎn)速為9 210 r/min,密封外徑壓力設(shè)定0.60 MPa,環(huán)境溫度設(shè)定25 ℃,密封內(nèi)徑壓力設(shè)定0.10 MPa,計(jì)算所得燃油側(cè)平端面接觸式機(jī)械密封與人字槽液膜密封性能特征參數(shù)見圖6,人字槽液膜密封計(jì)算結(jié)果如圖7所示。

        由圖6可知,在計(jì)算工況下該人字槽型密封結(jié)構(gòu)泄漏量約0.138 mL/min,空化率約0.97%,端面平均溫度約117.6 ℃,計(jì)算泄漏量滿足相關(guān)設(shè)計(jì)要求,這與該密封端面結(jié)構(gòu)所具有的上游泵送效果密切相關(guān)。

        圖6 平端面密封與人字槽型密封性能特征參數(shù)

        與平端面接觸式機(jī)械密封進(jìn)行對比,雖然人字槽型液膜密封泄漏量大于平端面密封,但是人字槽型密封摩擦扭矩與摩擦功率均小于平端面接觸式機(jī)械密封,人字槽型密封端面平均溫度低于平端面密封端面平均溫度約50 ℃,端面溫度的降低可相當(dāng)程度減小密封介質(zhì)高溫變質(zhì)、積碳發(fā)生和密封環(huán)高溫?zé)崃?、皰疤的可能性。另外,與平端面接觸式機(jī)械密封相比,由于人字槽型液膜密封摩擦副之間存在一層薄膜流體,在運(yùn)行工況發(fā)生波動時,薄膜流體具有的剛度特性可實(shí)現(xiàn)端面膜厚的自動調(diào)整,以適應(yīng)外界工況的變化;在密封閉合力突然增大時,密封端面薄膜流體的存在可起到減沖緩阻作用,避免密封摩擦副可能出現(xiàn)的接觸撞擊,有利于密封運(yùn)行穩(wěn)定性的提高。

        分析圖7(a)可以得出,雙環(huán)帶人字槽型密封端面間流體膜的壓力分布有2個低壓區(qū)域,分別出現(xiàn)在人字槽的外徑側(cè)槽的槽根處及內(nèi)徑側(cè)槽的槽根處,最低壓力為流體在計(jì)算溫度下的飽和蒸汽壓;在人字槽型的岔口尖端處出現(xiàn)高壓區(qū),其最大壓力約為1.80 MPa。人字槽型液膜密封相當(dāng)于在同一密封端面開設(shè)不同的兩組螺旋槽,在動環(huán)旋轉(zhuǎn)以及內(nèi)外徑壓差的作用下,端面流體具有較高的動能和壓能,當(dāng)流體在隨著動環(huán)端面旋轉(zhuǎn)的情況下遇到人字槽岔口尖端壁面,降速升壓,故此處出現(xiàn)高壓區(qū)。由于端面高壓的存在使密封動、靜環(huán)端面間充滿一層厚度為微米級的薄膜流體而維持密封端面摩擦副的非接觸狀態(tài),從而避免了端面固體直接接觸而引起的密封副摩擦磨損現(xiàn)象,降低了密封環(huán)因摩擦磨損而出現(xiàn)熱裂、皰疤、崩邊等密封失效行為的可能性;且由于該結(jié)構(gòu)通過端面開設(shè)的動壓槽在旋轉(zhuǎn)條件下的黏性剪切作用,把高壓側(cè)泄漏到低壓側(cè)的被密封介質(zhì)再反輸至高壓側(cè),降低了被密封介質(zhì)由高壓側(cè)向低壓側(cè)的泄漏,大大提高了密封運(yùn)行可靠性、穩(wěn)定性。

        圖7(b)所示為預(yù)設(shè)計(jì)結(jié)構(gòu)參數(shù)下人字槽型密封端面流體溫度分布??擅黠@地看出,密封端面流體溫度基本在120~140 ℃左右,靠近外徑側(cè)槽區(qū)溫度略低于靠近內(nèi)側(cè)處槽區(qū)溫度,這是由于靠近內(nèi)徑槽區(qū)面積較大,流體動壓效應(yīng)強(qiáng),流體與槽壩的摩擦熱高于靠近外徑處溫度。另外雖然該結(jié)構(gòu)密封端面存在一層薄膜流體避免了密封副的直接固體接觸,相對減少了摩擦生熱功率并使密封端面溫度降低,但由于運(yùn)行工況的高轉(zhuǎn)速與密封燃油介質(zhì)具有較大的黏度,故密封端面溫度仍較大于密封內(nèi)外徑處介質(zhì)溫度。

        基于能量守恒,對密封動靜環(huán)及周圍流體進(jìn)行傳熱計(jì)算,得到人字槽型密封動環(huán)、靜環(huán)溫度分布,如圖7(c)、(d)所示??梢园l(fā)現(xiàn),密封動靜環(huán)溫度最大值均位于密封摩擦副接觸端面處,但動環(huán)溫度最值低于靜環(huán)溫度相應(yīng)最值,這是由于動環(huán)材料9Cr18的導(dǎo)熱系數(shù)高于靜環(huán)材料M298的導(dǎo)熱系數(shù),增大了熱量在動環(huán)內(nèi)部的傳遞,且由于動環(huán)隨轉(zhuǎn)軸的運(yùn)動增大了動環(huán)與周圍流體的對流換熱,故雖然密封端面流體黏性剪切熱量流入動環(huán)量大于密封靜環(huán),但動環(huán)內(nèi)大部分區(qū)域溫度仍略低于靜環(huán)溫度。同時可看出,動環(huán)內(nèi)溫度梯度分布明顯大于靜環(huán)內(nèi)溫度梯度分布,這也是由動環(huán)材料9Cr18的導(dǎo)熱系數(shù)高于靜環(huán)材料M298的導(dǎo)熱系數(shù)所決定的。另外,雖然密封外徑側(cè)對流傳熱系數(shù)高于密封內(nèi)徑側(cè),但由于密封外徑側(cè)介質(zhì)溫度大于內(nèi)徑處溫度,故密封動靜環(huán)外徑側(cè)溫度仍均相對高于內(nèi)徑側(cè)溫度。

        圖7(e)、(f)所示為人字槽型密封動環(huán)、靜環(huán)軸向力熱變形??梢园l(fā)現(xiàn),密封動環(huán)最大變形發(fā)生在密封摩擦副接觸端面外徑處,該處邊界條件為液膜壓力與溫度分布;密封靜環(huán)最大變形發(fā)生在靜環(huán)遠(yuǎn)離接觸端面處,該處為軸向約束邊界條件;密封靜環(huán)最小變形發(fā)生在密封摩擦副接觸端面內(nèi)徑處,該處邊界條件為液膜壓力與溫度分布。分析可知,由于密封動環(huán)材料9Cr18的彈性模量遠(yuǎn)大于密封靜環(huán)材料M298的彈性模量,因此密封動環(huán)發(fā)生力熱變形能力弱于密封靜環(huán),故雖然在密封摩擦副接觸端面處具有相同的力熱邊界條件,但靜環(huán)端面最小變形量6.26 μm仍大于動環(huán)端面最大變形量3.72 μm。根據(jù)密封動、靜環(huán)端面變形計(jì)算可知,密封環(huán)變形量與密封端面流體膜厚、端面動壓槽深皆為微米量級,密封環(huán)變形的發(fā)生將通過對端面流體膜厚的影響而進(jìn)一步對密封性能產(chǎn)生影響。因此在模擬分析時采用了熱流固雙向耦合的方法,即通過計(jì)算軟件MatLab與ANSYS的相互調(diào)用,實(shí)現(xiàn)密封端面流場、溫度場、密封環(huán)變形的雙向計(jì)算,提高了該仿真計(jì)算的精確性。

        圖7 人字槽型密封熱流固耦合計(jì)算結(jié)果

        4 試驗(yàn)驗(yàn)證分析

        試驗(yàn)總時長2 h,為穩(wěn)態(tài)試驗(yàn)過程,轉(zhuǎn)速9 210 r/min,介質(zhì)壓力為0.5 MPa。對泄漏量、工裝處溫度、工裝處介質(zhì)溫度及軸套溫度進(jìn)行測量和記錄,每30 min記錄一次數(shù)據(jù)。試驗(yàn)前后密封環(huán)對比如圖8所示。

        圖8 試驗(yàn)前后密封環(huán)端面情況

        從圖8(b)中可以發(fā)現(xiàn):動環(huán)端面磨痕呈周向布置,均勻、無斷續(xù)痕跡;端面有明晰的3個區(qū)域,對應(yīng)外側(cè)壩區(qū)、中間槽區(qū)、內(nèi)徑壩區(qū);對應(yīng)中間槽區(qū)部分磨痕較兩側(cè)為重,符合啟停過程中槽邊刮削所致,但無任何崩邊、掉塊等接觸式密封常見問題。如圖8(b)所示,動環(huán)內(nèi)徑側(cè)劃痕明顯,這與圖7(f)的靜環(huán)變形方向及規(guī)律相對應(yīng),靜環(huán)軸向變形形成沿內(nèi)徑到外徑的收斂形錐度,導(dǎo)致動環(huán)內(nèi)徑處接觸比壓更大,磨痕相對明顯。這一定程度驗(yàn)證理論計(jì)算結(jié)果的準(zhǔn)確性。

        如圖9所示是溫度和泄漏量隨時間的變化。在試驗(yàn)過程中,各處溫度穩(wěn)定,無較大波動。未觀察到有泄漏產(chǎn)生,泄漏量始終為0,可滿足密封使用需求。

        圖9 溫度和泄漏量隨時間的變化

        5 結(jié)論

        (1)以描述粗糙表面的流體動壓潤滑方程、熱效應(yīng)能量方程、粗糙表面微凸體彈塑性接觸力方程、固體傳熱學(xué)方程及固體彈性力學(xué)方程等基本理論為指導(dǎo),建立了適用于航空發(fā)動機(jī)燃油增壓泵用機(jī)械密封的熱流固多物理場耦合理論模型,基于計(jì)算軟件MATLAB與ANSYS的耦合求解原理確定了仿真計(jì)算流程,給出了機(jī)械密封動靜環(huán)的計(jì)算物理模型、計(jì)算工況參數(shù)、密封環(huán)計(jì)算邊界條件,并使用數(shù)字式黏度計(jì)測量、擬合分析出燃油黏溫關(guān)系方程。

        (2)依據(jù)所建立的燃油增壓泵機(jī)械密封熱流固耦合理論模型,開展了人字槽型密封的仿真分析,與平端面接觸式機(jī)械密封進(jìn)行對比,人字槽型密封端面平均溫度降低約50 ℃,有利于密封運(yùn)行穩(wěn)定性的提高。

        (3)對人字槽型機(jī)械密封進(jìn)行了驗(yàn)證試驗(yàn)。密封測試過程中,各處溫度始終保持在穩(wěn)定范圍,且未觀察到有泄漏產(chǎn)生,泄漏量始終為0,滿足航空發(fā)動機(jī)燃油增壓泵對密封性能的要求。

        猜你喜歡
        機(jī)械變形
        談詩的變形
        中華詩詞(2020年1期)2020-09-21 09:24:52
        機(jī)械革命Code01
        電腦報(2020年35期)2020-09-17 13:25:53
        調(diào)試機(jī)械臂
        ikbc R300機(jī)械鍵盤
        電腦報(2019年40期)2019-09-10 07:22:44
        “我”的變形計(jì)
        簡單機(jī)械
        變形巧算
        例談拼圖與整式變形
        會變形的餅
        機(jī)械班長
        日本在线一区二区三区观看 | 中文人妻av久久人妻18| 国产精品一区二区韩国AV | 日韩女优图播一区二区| av中文字幕潮喷人妻系列| 精品亚洲a∨无码一区二区三区| 18禁裸男晨勃露j毛网站| 亚洲色大网站www永久网站| 国产成人精品三级麻豆| 蜜桃视频中文字幕一区二区三区| 日韩av中文字幕波多野九色| 国产午夜av秒播在线观看| 黑人玩弄人妻中文在线| 精品91精品91精品国产片| 免费看片的网站国产亚洲| 寂寞人妻渴望被中出中文字幕| 农村欧美丰满熟妇xxxx| 国产精品亚洲成在人线| 玖玖资源站无码专区| 青青草原综合久久大伊人精品| 熟女少妇精品一区二区| 98bb国产精品视频| 少妇激情一区二区三区久久大香香 | 女人被爽到高潮视频免费国产| 人与嘼av免费| 激情视频在线观看国产中文| 精品亚洲天堂一区二区三区| 后入到高潮免费观看| 国产免费一区二区三区在线观看| 蜜臀av一区二区三区人妻在线| 射进去av一区二区三区| 免费观看a级片| 国产精品午睡沙发系列| 丰满少妇高潮在线观看| 在线视频国产91自拍| 麻豆一区二区99久久久久| 国产剧情无码中文字幕在线观看不卡视频 | 93精91精品国产综合久久香蕉| 中文字幕二区三区在线| 国产大屁股视频免费区| 国精产品一品二品国在线|