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

    基于逾滲理論的液體火箭發(fā)動(dòng)機(jī)管路接頭密封性能研究

    2023-01-07 05:41:50金輝趙有磊尚現(xiàn)偉金平蔡國(guó)飆
    強(qiáng)度與環(huán)境 2022年6期
    關(guān)鍵詞:管路尺度介質(zhì)

    金輝 趙有磊 尚現(xiàn)偉 金平 蔡國(guó)飆

    (1 北京航空航天大學(xué) 宇航學(xué)院,北京 102206;2 北京航天動(dòng)力研究所,北京 100076)

    0 引言

    液體火箭發(fā)動(dòng)機(jī)管路接頭的密封問(wèn)題屬于靜密封的研究范疇,靜密封的定義為密封界面不存在明顯的相對(duì)滑動(dòng)運(yùn)動(dòng)[1]。密封部位發(fā)生泄漏將給火箭帶來(lái)巨大安全隱患,2010年發(fā)現(xiàn)號(hào)航天飛機(jī)的氫燃料排放管路發(fā)生泄漏導(dǎo)致發(fā)射推遲[2]。發(fā)動(dòng)機(jī)中的渦輪泵端面密封、閥門、管路接頭和法蘭等零組件中采用了靜密封方式來(lái)阻止介質(zhì)發(fā)生泄漏。在預(yù)緊力的作用下,如圖1所示的待密封兩表面之間發(fā)生接觸,從而形成靜密封區(qū)域。液體火箭發(fā)動(dòng)機(jī)的靜密封結(jié)構(gòu)面臨著諸多挑戰(zhàn):密封部位處于多場(chǎng)耦合作用狀態(tài),介質(zhì)溫差大、壓力高和結(jié)構(gòu)振動(dòng)等[3],這些復(fù)雜的載荷條件加劇了發(fā)動(dòng)機(jī)靜密封失效的風(fēng)險(xiǎn)。

    圖1 液體火箭發(fā)動(dòng)機(jī)的靜密封結(jié)構(gòu)[4],1-法蘭端面,2-墊片,3-接頭體,4-螺母,5-接頭,6-活門,7-靜環(huán), 8-墊圈,9-彈簧10-殼體,11-活門 Fig.1 Liquid rocket engine turbo pump components and valves

    大量研究對(duì)靜密封結(jié)構(gòu)的接觸過(guò)程開(kāi)展數(shù)值計(jì)算和試驗(yàn),獲得密封面的接觸壓力、接觸面積和接觸長(zhǎng)度,并且以此作為密封性能判斷準(zhǔn)則。韓沖[5]對(duì)火箭發(fā)動(dòng)機(jī)60°球形管路接頭的靜密封進(jìn)行研究,分析了密封面寬度、平均接觸壓力與擰緊力矩之間的關(guān)系。進(jìn)一步的,楊卓然和Van Zyl G[6-8]等人給出了接觸壓力和接觸長(zhǎng)度下限,認(rèn)為高于下限值才能滿足形成可靠的密封。上述的研究都忽略了靜密封接觸部位的表面形貌,實(shí)際上從微觀尺度來(lái)看,接觸面是由許多粗糙峰和凹陷谷構(gòu)成的不規(guī)則粗糙面[9],從微觀角度進(jìn)行靜密封的接觸分析有助于進(jìn)一步提高對(duì)靜密封形成過(guò)程的理解。閆洋洋[10]采用多尺度有限元模型分析了卡套式接頭在擰緊過(guò)程中的密封狀態(tài)和特性變化。Patel H[11]假設(shè)法蘭墊片的表面由高度不規(guī)則的球形微凸體構(gòu)成,認(rèn)為密封性能取決于預(yù)緊力、流體壓力和彎矩等因素。Liua M[12]將靜密封兩表面等效為剛性表面與分形表面,并以此為基礎(chǔ)分析了泄漏機(jī)理。在獲得各影響因素對(duì)靜密封性能影響規(guī)律的基礎(chǔ)上,研究者們希望從形成密封的原理出發(fā),得到滿足靜密封性能要求的條件。Broadbent S R[13]提出逾滲理論來(lái)分析流體介質(zhì)的逾滲特性,為靜密封性能的評(píng)價(jià)提供理論基礎(chǔ)。所謂逾滲是指體系以外的一種介質(zhì)通過(guò)一定的路徑進(jìn)入體系內(nèi)的過(guò)程。該理論與流體自由擴(kuò)散的區(qū)別在于:認(rèn)為流體的運(yùn)動(dòng)不具有隨機(jī)性,而逾滲結(jié)構(gòu)具有隨機(jī)性,結(jié)構(gòu)的形狀和分布決定了流體的運(yùn)動(dòng)方向。Bottiglione F和張彤[14-15]等人認(rèn)為密封界面并不是光滑的,當(dāng)粗糙面的接觸面積與光滑面的總面積之比(即接觸面積比)小于一定比例時(shí)就會(huì)出現(xiàn)臨界泄漏通道。進(jìn)一步的,史建成[16]利用種子單元擴(kuò)張的方法,計(jì)算得到出現(xiàn)臨界泄漏通道的接觸面積比值收斂于0.41。Ying Cui[17]建立粗糙表面的有限元接觸模型,計(jì)算得到整個(gè)密封面無(wú)泄漏通道對(duì)應(yīng)的接觸壓力,并將接觸壓力作為密封性能的評(píng)價(jià)指標(biāo)。

    液體火箭發(fā)動(dòng)機(jī)靜密封部位充斥著工況惡劣的流體介質(zhì),本文選取發(fā)動(dòng)機(jī)中應(yīng)用靜密封形式的管路接頭作為研究對(duì)象,開(kāi)展微觀尺度和宏觀尺度有限元接觸分析。首先通過(guò)建立靜密封微尺度接觸模型,得到接觸面積比和平均接觸壓力的關(guān)系,然后基于逾滲理論得到微尺度臨界平均接觸壓力(簡(jiǎn)稱臨界接觸壓力)這一性能評(píng)價(jià)指標(biāo),并將應(yīng)用該指標(biāo)分析得到的結(jié)論和試驗(yàn)進(jìn)行對(duì)比。最后將該指標(biāo)應(yīng)用于宏觀尺度失效分析的評(píng)價(jià)當(dāng)中,分析預(yù)緊力、介質(zhì)壓力和材料搭配對(duì)管路接頭靜密封性能的影響,為發(fā)動(dòng)機(jī)管路接頭的裝配與使用提供理論指導(dǎo)。

    1 基于逾滲理論的靜密封性能分析方法

    1.1 靜密封表面的形貌模擬

    粗糙面的三維形貌是無(wú)序的,可以通過(guò)統(tǒng)計(jì)學(xué)方法來(lái)描述。提取表面的統(tǒng)計(jì)學(xué)信息,如形狀誤差、粗糙度特征和波紋度等[18],并通過(guò)統(tǒng)計(jì)學(xué)方法擬合出表面形貌。本文參考相關(guān)文獻(xiàn)[9][18]認(rèn)為管路接頭的加工表面具有各向同性Gauss分布表面特征。構(gòu)建Gauss粗糙表面的關(guān)鍵步驟是建立Gauss數(shù)組內(nèi)各相鄰元素的空間相關(guān)性。Gauss數(shù)組內(nèi)每一點(diǎn)高度的概率服從Gauss分布,但這些點(diǎn)在空間上沒(méi)有相關(guān)性,需要依靠構(gòu)造自相關(guān)函數(shù)ACF來(lái)組成連續(xù)的表面,自相關(guān)函數(shù)描述了隨機(jī)數(shù)組內(nèi)每一點(diǎn)的空間分布。將尺寸為L(zhǎng)×L的正方形區(qū)域沿著x、y軸劃分為N×N個(gè)網(wǎng)格單元,所有單元節(jié)點(diǎn)的高度組成服從Gauss分布的二維數(shù)組z0,經(jīng)過(guò)式(1)的變換后數(shù)組z0成為具有連續(xù)性的數(shù)組z。

    其中z0=θ·randn(N,N),θ為數(shù)組z0的均方根,randn(N,N)為相關(guān)性等于0的二維Gauss數(shù)組。自相關(guān)函數(shù)ACF的表達(dá)式如(2)所示:

    其中x和y為網(wǎng)格單元節(jié)點(diǎn)的坐標(biāo)。cl為相關(guān)長(zhǎng)度,代表兩點(diǎn)相關(guān)性。FFT和FFT-1分別為快速傅里葉變換及其逆變換:

    由于FFT對(duì)有限數(shù)組的變換存在截?cái)嗾`差,需要使用偏態(tài)和峰度對(duì)求得的隨機(jī)粗糙表面進(jìn)行驗(yàn)證。偏態(tài)Sk是描述統(tǒng)計(jì)分布對(duì)稱性的指標(biāo),Gauss分布數(shù)組的偏態(tài)Sk=0,這意味著表面輪廓的粗糙峰和凹陷谷均勻分布。峰度Ku是對(duì)高度分布曲線振蕩程度的度量,Gauss分布數(shù)組的峰度Ku=3[23]。

    1.2 逾滲理論

    在無(wú)序系統(tǒng)中,構(gòu)成系統(tǒng)的各元素之間關(guān)聯(lián)程度變化,引起系統(tǒng)的某一特性發(fā)生階躍,逾滲描述的正是這一階躍現(xiàn)象[21]。由此可以看出適用逾滲理論處理的系統(tǒng)具有兩個(gè)特點(diǎn)(1)系統(tǒng)具有隨機(jī)性(2)系統(tǒng)的元素之間存在某種聯(lián)系。逾滲理論給復(fù)雜、無(wú)序空間內(nèi)的連通性相關(guān)問(wèn)題的研究提供了新的思路。逾滲理論認(rèn)為固體域的結(jié)構(gòu)形狀決定了流體的流動(dòng)方向,但實(shí)際上流體會(huì)導(dǎo)致固體域發(fā)生微小形變,從而導(dǎo)致流體的流動(dòng)區(qū)域發(fā)生改變。當(dāng)忽略流體域?qū)腆w域的作用,逾滲理論可以處理無(wú)序系統(tǒng)的連通性問(wèn)題。

    微尺度下靜密封的兩個(gè)接觸表面由如圖2a)所示的許多粗糙峰和凹陷谷構(gòu)成,不平整的兩表面存在如圖2b)所示的黑色接觸部位和白色未接觸部位,未接觸的部位形成空隙,連續(xù)的空隙構(gòu)成了流體介質(zhì)的泄漏通道,未接觸部位的數(shù)量越多意味著密封部位發(fā)生泄漏的可能性越大,尋找泄漏通道為靜密封的性能分析提供了思路。從微尺度觀察密封接觸面的每一區(qū)域,區(qū)域內(nèi)是否發(fā)生接觸存在隨機(jī)性,因此靜密封滿足逾滲理論的假設(shè)條件,逐一分析粗糙表面的未接觸點(diǎn)就能夠判斷是否存在連通的泄漏通道。

    圖2 微尺度下靜密封接觸部位示意圖 Fig.2 Schematic diagram of the static seal contact part at the microscale

    將密封界面離散成N×N個(gè)如圖3所示的正方形網(wǎng)格單元,每一個(gè)網(wǎng)格單元只有兩種狀態(tài),分別是黑色代表的“占有狀態(tài)”和綠色代表的“空狀態(tài)”,并且每個(gè)網(wǎng)格單元被占有的概率均為γ,則為空狀態(tài)的概率為1-γ。參考史建成[16]對(duì)靜密封界面狀態(tài)演變特性的研究,介質(zhì)以四連通的方式流動(dòng)。假設(shè)粗糙界面的微尺度接觸面積為A,宏觀尺度接觸面積為An,接觸面積比α=A/An,其中α∈[0,1]。當(dāng)α≈0時(shí),存在貫穿入口至出口的綠色網(wǎng)格單元簇,這意味著界面發(fā)生泄漏,而隨著α的增大,介質(zhì)能夠流過(guò)的綠色網(wǎng)格單元數(shù)量不斷減少,因此存在臨界接觸面積αc,當(dāng)α<αc時(shí)界面不存在貫穿的流動(dòng)通道。文獻(xiàn)[22]中分析了當(dāng)網(wǎng)格數(shù)量趨近于無(wú)限大時(shí),界面不存在泄漏通道的臨界接觸面積比收斂于αc≈0.41。

    圖3 逾滲理論示意圖 Fig.3 Schematic diagram of percolation theory

    1.3 靜密封性能評(píng)價(jià)指標(biāo)

    評(píng)價(jià)靜密封性能的指標(biāo)來(lái)源于接觸面的信息,包括接觸面的應(yīng)力和長(zhǎng)度。接觸壓力起到構(gòu)建密封區(qū)域的作用,兩表面從初始接觸到發(fā)生擠壓變形的過(guò)程中接觸壓力不斷升高,一般認(rèn)為接觸壓力達(dá)到臨界值即可滿足密封的要求。如圖4所示為計(jì)算靜密封性能評(píng)價(jià)指標(biāo)的框架圖,在宏觀尺度接觸部位內(nèi)選取局部區(qū)域作為微尺度的研究對(duì)象,找到滿足逾滲理論要求的臨界接觸面積比,然后得到面積比對(duì)應(yīng)的微尺度臨界平均接觸壓力。由于微觀尺度模型的面積遠(yuǎn)小于宏觀尺度接觸面的面積,可將微觀區(qū)域視為宏觀尺度接觸面上的某一點(diǎn),因此微尺度平均接觸壓力等于宏觀尺度接觸面上任意一點(diǎn)的接觸壓力P,當(dāng)宏觀尺度接觸面上某一點(diǎn)的接觸壓力時(shí),流體介質(zhì)無(wú)法通過(guò)該點(diǎn),考慮到管路接頭具有中心軸對(duì)稱的特征,則該點(diǎn)經(jīng)過(guò)中心軸旋轉(zhuǎn)之后構(gòu)成了一條阻礙介質(zhì)泄漏的密封線。

    圖4 靜密封性能評(píng)價(jià)指標(biāo)的框架圖 Fig.4 Frame diagram of static sealing performance evaluation index

    2 靜密封有限元接觸模型

    2.1 微觀尺度接觸模型

    以液體火箭發(fā)動(dòng)機(jī)中常用37°DN6球形管路接頭作為靜密封性能分析的研究對(duì)象,選取宏觀尺度兩接觸面的部分區(qū)域建立微尺度接觸模型,計(jì)算微尺度平均接觸壓力。管路接頭的表面粗糙度Ra=1.6um,各輸入?yún)?shù)如表1所示。

    計(jì)算得到粗糙表面的偏態(tài)Skz=-0.0052,峰度Kuz=2.9985,與Gauss分布特征表面的理論值相比誤差都在1%以內(nèi),因此認(rèn)為生成的隨機(jī)表面服從Gauss分布。將隨機(jī)粗糙表面的高度坐標(biāo)z導(dǎo)入軟件Comsol中,通過(guò)插值擬合出連續(xù)的表面,再經(jīng)過(guò)建模軟件的處理得到如圖5所示上表面粗糙的長(zhǎng)方體。

    表1 Gauss分布粗糙面輸入?yún)?shù) Table1 Gauss distribution rough surface input parameters

    考慮到接觸模型仿真的計(jì)算規(guī)模,從尺寸為250μm×250μm的粗糙體中截取尺寸為100μm×100μm的區(qū)域作為仿真對(duì)象。在軟件ANSYS中對(duì)如圖6所示兩個(gè)表面粗糙的長(zhǎng)方體進(jìn)行有限元接觸仿真,上長(zhǎng)方體(上體)和下長(zhǎng)方體(下體)的材料均為雙線性各向同性硬化本構(gòu)的不銹鋼,材料的泊松比υ=0.31,彈性模量E=1.93GPa,屈服極限σs=210MPa,切線模量G=180MPa。采用四面體網(wǎng)格劃分,非接觸區(qū)域的網(wǎng)格單元類型為Solid187,上體的粗糙面作為接觸面,單元尺寸為0.7μm,網(wǎng)格單元類型為Conta174。下體的粗糙面作為目標(biāo)面,單元尺寸為0.9μm,網(wǎng)格單元類型為Targe170,仿真模型的網(wǎng)格單元總數(shù)為616997。固定下體的底面,為了使得計(jì)算過(guò)程更易于收斂,在上體的上平面施加強(qiáng)制位移,上體和下體之間為無(wú)摩擦接觸,采用無(wú)摩擦支撐分別約束上、下體的其它四個(gè)側(cè)面。

    圖5 隨機(jī)粗糙表面模擬結(jié)果圖 Fig.5 Random rough surface simulation result

    2.2 宏觀尺度接觸模型

    如圖7所示為管路接頭的結(jié)構(gòu)示意圖,由球形接頭、接頭體和螺母三部分構(gòu)成,其中螺母對(duì)球形接頭的后端面施加軸向預(yù)緊力,使得接頭體和球形接頭靠近、接觸直至相互擠壓。由此可知管路接頭的靜密封區(qū)域是球形接頭和接頭體的擠壓區(qū)域,因此在仿真的過(guò)程中忽略螺母具體結(jié)構(gòu),僅保留其預(yù)緊功能。如圖7所示為管路接頭的主要結(jié)構(gòu)尺寸實(shí)體圖,包括通徑d,接頭體的錐面角β和球頭半徑R。

    圖6 微尺度接觸模型及其邊界條件 Fig.6 Microscale contact model and its boundary conditions

    圖7 管路接頭關(guān)鍵尺寸示意圖 Fig.7 Schematic diagram of key dimensions of pipe joint

    管路接頭在發(fā)動(dòng)機(jī)內(nèi)工作時(shí)受到的載荷是軸對(duì)稱的,因此可將三維模型簡(jiǎn)化為二維軸對(duì)稱模型。如圖8所示為管路接頭宏觀尺度接觸模型,采用平面單元進(jìn)行網(wǎng)格劃分,非接觸區(qū)域的單元類型為Plane182,球形接頭的球面是接觸面,其單元類型為Conta172,接頭體的錐面為目標(biāo)面,其單元類型為Trage169。本文主要關(guān)注接觸部位的信息,因此對(duì)接觸部位的網(wǎng)格加密處理,球形接頭和接頭體的加密尺寸分別為15μm、30μm。固定約束接頭體的下端面,在球形接頭的左端面施加軸向預(yù)緊力,接觸面之間為摩擦接觸,摩擦系數(shù)取為0.15,并且在管路接頭的內(nèi)部分散著流體壓力載荷。

    3 結(jié)果與討論

    3.1 Gauss粗糙面參數(shù)對(duì)幾何形貌的影響

    圖9所示為0.25mm×0.25mm區(qū)域內(nèi),不同粗糙度Ra和自相關(guān)長(zhǎng)度cl對(duì)應(yīng)的粗糙面高度分布云圖,高度值越大的區(qū)域顏色越深。對(duì)比圖9(a)-圖9(c),當(dāng)表面粗糙度相等時(shí),隨著自相關(guān)長(zhǎng)度的增加,不同顏色區(qū)域的面積不斷升高,這說(shuō)明自相關(guān)長(zhǎng)度能夠增加粗糙表面的連續(xù)性。對(duì)比圖9(b)與圖9(d),當(dāng)表面的自相關(guān)長(zhǎng)度相等時(shí),可以看出隨著表面粗糙度的增加,表面擁有更大的高度幅值和高度差,這意味著粗糙度較大的表面不平坦。

    圖8 宏觀尺度接觸模型及其邊界條件 Fig.8 Macroscale contact model and its boundary conditions

    圖9 不同粗糙度和自相關(guān)長(zhǎng)度的表面高度分布云圖 Fig.9 Surface height distribution cloud images with different roughness and autocorrelation length

    3.2 密封性能評(píng)價(jià)指標(biāo)的取值

    圖10 平均接觸壓力與接觸面積比的關(guān)系曲線 Fig.10 The relationship between the average contact stress and the contact area ratio

    圖11所示為兩種不同接觸面積比對(duì)應(yīng)的密封面接觸情況,隨著接觸面積比α的增大,圖中接觸區(qū)域的面積也隨之增大。如圖11(a)所示,α=0.3的區(qū)域內(nèi)存在大量未接觸的空白區(qū)域,這些連通的空白區(qū)域會(huì)構(gòu)成泄漏通道,而圖11(b)中α=0.41的區(qū)域內(nèi)無(wú)明顯的泄漏通道。綜上可知,由不銹鋼材料加工而成的靜密封部位,只有當(dāng)宏觀尺度接觸區(qū)域的接觸壓力P≥520MPa時(shí)才能形成良好的密封效果。

    對(duì)37°DN6球形管路接頭開(kāi)展氣密性試驗(yàn),將試驗(yàn)件浸入水中,用空氣給試驗(yàn)件緩慢加壓到45MPa高壓,保持5min,試驗(yàn)過(guò)程中觀察接頭連接處是存在氣泡,以及氣泡的數(shù)量。7組試驗(yàn)的結(jié)果表明,當(dāng)預(yù)緊扭矩為30~35N·m之間時(shí),管路接頭無(wú)氣泡冒出。參考文獻(xiàn)[5]的方法計(jì)算得到30N·m扭矩對(duì)應(yīng)的軸向預(yù)緊力為9kN左右,該預(yù)緊條件下接觸面的最大接觸壓力P=623MPa,大于微尺度臨界平均接觸壓力,利用本文提出的密封性能評(píng)價(jià)方法得到的結(jié)論和試驗(yàn)一致。

    3.3 預(yù)緊力對(duì)密封性能的影響

    預(yù)緊力F對(duì)靜密封區(qū)域的形成起著決定性的作用。圖12為預(yù)緊過(guò)程中接觸長(zhǎng)度和接觸區(qū)域Von Mises等效應(yīng)力的變化規(guī)律。管接頭預(yù)緊過(guò)程中接觸長(zhǎng)度隨預(yù)緊力線性增加,這表明預(yù)緊力能夠增強(qiáng)接頭的密封性能。在預(yù)緊初始階段F=0.5kN 時(shí),從接觸部位的等效應(yīng)力云圖可以看出,接觸部位的最大等效應(yīng)力σmax=220MPa ,超過(guò)材料的屈服極限,這是因?yàn)榍蛎婧湾F面的初始接觸類型是點(diǎn)、線接觸,接觸部位存在應(yīng)力集中的現(xiàn)象。隨著預(yù)緊力的增大,接觸部位變成線、線接觸,高應(yīng)力區(qū)面積不斷增加,但高應(yīng)力區(qū)內(nèi)的球面和錐面將發(fā)生不可逆轉(zhuǎn)的塑性變形,塑性變形削弱了管路接頭的回彈能力,因此需要找到密封要求對(duì)應(yīng)的最小預(yù)緊力,減小預(yù)緊力對(duì)密封結(jié)構(gòu)的破壞。如圖13所示為預(yù)緊過(guò)程中的接觸壓力分布情況,隨著預(yù)緊力的增大,接觸壓力的幅值和分布長(zhǎng)度不斷增加,并且高接觸壓力區(qū)主要集中在接觸部位,當(dāng)F=3kN時(shí)接觸壓力P>520MPa,這意味著能夠形成密封區(qū)域。

    圖11 不同接觸面積比的密封面接觸情況 Fig.11 The contact situation of the sealing surface with different contact area ratio

    圖12 預(yù)緊過(guò)程中接觸長(zhǎng)度和應(yīng)力面積Fig.12 Contact length and stress area during assembly process

    圖13 預(yù)緊過(guò)程中接觸壓力Fig.13 Contact pressure during assembly process

    3.4 介質(zhì)壓力對(duì)密封性能的影響

    液體火箭發(fā)動(dòng)機(jī)的流體介質(zhì)處于超高壓狀態(tài),增加了發(fā)動(dòng)機(jī)靜密封的泄漏風(fēng)險(xiǎn),分析介質(zhì)壓力pf對(duì)管路接頭密封部位的影響。在上述預(yù)緊力對(duì)管路接頭密封性能的分析過(guò)程中,認(rèn)為無(wú)介質(zhì)壓力情況下預(yù)緊力F=3kN即可滿足管路接頭的密封要求。選取F=3kN,介質(zhì)壓力pf∈ [0,40]MPa 的工況分析介質(zhì)壓力對(duì)管路接頭密封性能的影響。如圖14所示接觸部位的長(zhǎng)度和應(yīng)力隨介質(zhì)壓力的變化情況,當(dāng)pf∈ [0,15]MPa 時(shí),接觸長(zhǎng)度保持不變,高應(yīng)力區(qū)面積減小,高應(yīng)力區(qū)的形狀也由無(wú)介質(zhì)壓力時(shí)的橢圓形變得不規(guī)則,兩者的最大等效應(yīng)力都位于靠近流體介質(zhì)的一側(cè),并且球面的高應(yīng)力區(qū)面積比錐面的減少更多。當(dāng)pf∈[20,30]MPa 時(shí)接觸長(zhǎng)度減少3%,高應(yīng)力區(qū)的形狀不變但面積進(jìn)一步減小,此時(shí)兩接觸面發(fā)生微小滑移運(yùn)動(dòng)。介質(zhì)壓力pf>30MPa 的接觸長(zhǎng)度急劇減小,pf=30MPa 對(duì)應(yīng)的接觸長(zhǎng)度減小16%,這意味著兩接觸表面發(fā)生的顯著相對(duì)運(yùn)動(dòng)。

    圖14 不同介質(zhì)壓力的接觸長(zhǎng)度和應(yīng)力面積Fig.14 Contact length and stress area for different media pressures

    如圖15所示為不同介質(zhì)壓力的接觸壓力分布變化情況,隨著介質(zhì)壓力的升高,密封區(qū)域的接觸壓力幅值和分布范圍都呈現(xiàn)下降的趨勢(shì),這是因?yàn)榻橘|(zhì)壓力抵消了一部分預(yù)緊力的作用,導(dǎo)致密封部位的性能降低,管路接頭在預(yù)緊力F= 3kN的條件下可以承受pf∈[0,10]MPa 的壓力,當(dāng)介質(zhì)壓力過(guò)大時(shí)只能通過(guò)增加預(yù)緊力的方式來(lái)提高管路接頭的密封性能。

    圖15 不同介質(zhì)壓力的接觸壓力Fig.15 Contact pressure for different media pressures

    3.5 材料搭配的影響

    球形接頭的材料為不銹鋼,接管嘴的材料分別為GH4169、鈦合金、不銹鋼和鋁合金,比較不同材料對(duì)管路接頭密封性能的影響,如表2所示為相關(guān)的材料參數(shù),材料本構(gòu)均為雙線性各向同性硬化。如圖16所示為不同材料在預(yù)緊力作用下的接觸長(zhǎng)度變化,當(dāng)預(yù)緊力F<3kN時(shí),鋁合金的接觸長(zhǎng)度明顯高于其它材料的接觸長(zhǎng)度,這是因?yàn)殇X合金的彈性模量在所有材料中最小,接觸面在擠壓的過(guò)程中容易發(fā)生變形,因此發(fā)生接觸的區(qū)域更多。當(dāng)F=8kN時(shí),鈦合金、不銹鋼和鋁合金的接觸長(zhǎng)度相等,相比之下,GH4169的接觸長(zhǎng)度比其他三種材料少5.7%。保持預(yù)緊力F=7kN不變,如圖17所示為預(yù)緊力相等的條件下,四種材料對(duì)應(yīng)的接觸壓力分布。屈服極限更大的GH4169和鈦合金,它們的接觸壓力幅值高于屈服極限較低的鋁合金與不銹鋼,并且接管嘴和球形接頭材料相同時(shí)的接觸壓力最低。因此為了提高管路接頭的密封性能,推薦使用不同材料搭配的接管嘴和球形接頭。

    表2 材料參數(shù) Table2 Material parameters

    圖16 不同材料的接觸長(zhǎng)度Fig.16 Contact length of different materials

    圖17 不同材料的接觸壓力Fig.17 Contact pressure of different materials

    4 結(jié)論

    1)具有Gauss分布特征的隨機(jī)粗糙表面,隨著自相關(guān)函數(shù)的增加,表面的連續(xù)性增強(qiáng)。粗糙度越大的隨機(jī)表面具有更大的高度差。建立表面粗糙度Ra=1.6um,自相關(guān)長(zhǎng)度cl=3um的兩不銹鋼粗糙表面有限元接觸模型,得到在臨界接觸面積αc=0.41對(duì)應(yīng)的臨界接觸壓力=520MPa,并將作為靜密封性能的評(píng)價(jià)標(biāo)準(zhǔn),應(yīng)用該指標(biāo)判斷密封性能和45MPa高壓氣密試驗(yàn)得到的結(jié)果相吻合。

    2)選取液體火箭發(fā)動(dòng)機(jī)DN6不銹鋼球形管路接頭的密封區(qū)域作為研究對(duì)象,建立宏觀尺度有限元接觸模型。在無(wú)壓力載荷的條件下,隨著預(yù)緊力的增大,接觸部位的高應(yīng)力區(qū)面積不斷升高,接觸壓力的分布范圍和接觸長(zhǎng)度也在增加,當(dāng)預(yù)緊力F=3kN時(shí)能夠滿足密封性能評(píng)價(jià)標(biāo)準(zhǔn)。介質(zhì)壓力對(duì)預(yù)緊力起到抵消的作用,并且改變了接觸部位的應(yīng)力分布形狀,使接觸壓力和接觸長(zhǎng)度下降,預(yù)緊力F=3kN的接頭可承受[0,10]MPa范圍內(nèi)的介質(zhì)壓力。不同屈服強(qiáng)度的材料搭配可以提高管路接頭的密封性能。

    猜你喜歡
    管路尺度介質(zhì)
    基于水質(zhì)變化的供熱采暖管路設(shè)計(jì)
    信息交流介質(zhì)的演化與選擇偏好
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    液壓管路系統(tǒng)隨機(jī)振動(dòng)下疲勞分析
    淬火冷卻介質(zhì)在航空工業(yè)的應(yīng)用
    硅鋼軋制過(guò)程中乳化液流量控制解耦研究及應(yīng)用
    山西冶金(2019年2期)2019-05-31 11:30:04
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    9
    美航天服漏水或因管路堵塞
    太空探索(2014年4期)2014-07-19 10:08:58
    考慮中間介質(zhì)換熱的廠際熱聯(lián)合
    久久久久久久午夜电影| 亚洲人成网站高清观看| 尾随美女入室| 午夜日本视频在线| 精品一区在线观看国产| 亚洲人与动物交配视频| a级一级毛片免费在线观看| 色网站视频免费| eeuss影院久久| 欧美性感艳星| 韩国高清视频一区二区三区| 免费观看精品视频网站| 欧美精品一区二区大全| 国国产精品蜜臀av免费| 99热6这里只有精品| 国产探花在线观看一区二区| 国产熟女欧美一区二区| 最近最新中文字幕免费大全7| 国产精品久久久久久av不卡| 狠狠精品人妻久久久久久综合| 尾随美女入室| 亚洲熟女精品中文字幕| 欧美zozozo另类| 一级黄片播放器| 国产欧美日韩精品一区二区| 少妇的逼好多水| 国产黄色小视频在线观看| a级毛色黄片| 少妇的逼水好多| 男的添女的下面高潮视频| 亚洲精品色激情综合| 欧美日韩国产mv在线观看视频 | 欧美一区二区亚洲| 美女高潮的动态| 91精品一卡2卡3卡4卡| 22中文网久久字幕| 一级毛片电影观看| 三级国产精品片| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 青春草亚洲视频在线观看| 国产成人aa在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲国产高清在线一区二区三| 色综合亚洲欧美另类图片| 日韩av在线大香蕉| 午夜精品国产一区二区电影 | 插逼视频在线观看| 免费av毛片视频| 亚洲国产最新在线播放| 国产成人a区在线观看| 国产av码专区亚洲av| 色吧在线观看| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验| 国产精品国产三级专区第一集| 久久精品综合一区二区三区| 在线免费观看不下载黄p国产| 秋霞伦理黄片| 久久99精品国语久久久| 69av精品久久久久久| 久久草成人影院| 国产精品99久久久久久久久| 日韩伦理黄色片| 国精品久久久久久国模美| 91精品国产九色| 伊人久久精品亚洲午夜| 街头女战士在线观看网站| 国产一区亚洲一区在线观看| 亚洲自偷自拍三级| 欧美极品一区二区三区四区| 午夜精品国产一区二区电影 | 久久久久九九精品影院| 亚洲av不卡在线观看| 国产在线一区二区三区精| 极品少妇高潮喷水抽搐| 成人一区二区视频在线观看| 一个人看的www免费观看视频| 亚洲精品自拍成人| 日日撸夜夜添| 欧美性感艳星| 国产老妇伦熟女老妇高清| 麻豆av噜噜一区二区三区| 直男gayav资源| 日本一本二区三区精品| 欧美日本视频| 色网站视频免费| freevideosex欧美| 国产成人91sexporn| 舔av片在线| 欧美成人一区二区免费高清观看| 久久精品夜色国产| 伊人久久国产一区二区| 日日摸夜夜添夜夜添av毛片| 寂寞人妻少妇视频99o| 亚洲自拍偷在线| 亚洲精品乱久久久久久| 一区二区三区四区激情视频| 国产人妻一区二区三区在| 中文字幕免费在线视频6| 国产午夜精品一二区理论片| 国产午夜福利久久久久久| 亚洲av中文av极速乱| 国产精品一及| 淫秽高清视频在线观看| 大片免费播放器 马上看| 熟女人妻精品中文字幕| 特大巨黑吊av在线直播| 成人综合一区亚洲| 欧美激情国产日韩精品一区| 插逼视频在线观看| av免费观看日本| 国产成人精品婷婷| 最后的刺客免费高清国语| 麻豆乱淫一区二区| 久久精品国产亚洲av天美| 亚洲丝袜综合中文字幕| 国产精品1区2区在线观看.| 国产精品国产三级专区第一集| 在线免费观看的www视频| or卡值多少钱| 性色avwww在线观看| 国产69精品久久久久777片| 国产色爽女视频免费观看| 一个人看视频在线观看www免费| 男的添女的下面高潮视频| 菩萨蛮人人尽说江南好唐韦庄| 国国产精品蜜臀av免费| 午夜老司机福利剧场| 色视频www国产| av在线老鸭窝| 国内精品一区二区在线观看| 亚洲欧美一区二区三区黑人 | 亚洲成人av在线免费| 亚洲乱码一区二区免费版| 亚洲国产欧美在线一区| 嫩草影院新地址| 又爽又黄a免费视频| 亚洲美女搞黄在线观看| 国产精品人妻久久久久久| 国产有黄有色有爽视频| 日韩人妻高清精品专区| 国产伦精品一区二区三区视频9| 午夜激情久久久久久久| 91在线精品国自产拍蜜月| 国产视频内射| 免费人成在线观看视频色| 欧美日韩综合久久久久久| 日本与韩国留学比较| 国产av不卡久久| 午夜老司机福利剧场| 国产91av在线免费观看| 亚洲精华国产精华液的使用体验| 纵有疾风起免费观看全集完整版 | 国产av国产精品国产| 国产美女午夜福利| 最近最新中文字幕大全电影3| 九草在线视频观看| 久久99蜜桃精品久久| 一边亲一边摸免费视频| 久久6这里有精品| 国产亚洲一区二区精品| 一级爰片在线观看| www.av在线官网国产| 蜜臀久久99精品久久宅男| 2018国产大陆天天弄谢| 嫩草影院入口| 亚洲精品中文字幕在线视频 | 国产色婷婷99| 亚洲国产高清在线一区二区三| 菩萨蛮人人尽说江南好唐韦庄| 97热精品久久久久久| 老女人水多毛片| 三级经典国产精品| 汤姆久久久久久久影院中文字幕 | 能在线免费观看的黄片| 亚洲在线自拍视频| 欧美区成人在线视频| 少妇丰满av| 欧美zozozo另类| 成人二区视频| 国产伦在线观看视频一区| 边亲边吃奶的免费视频| 丰满人妻一区二区三区视频av| 亚洲欧美日韩卡通动漫| 成年免费大片在线观看| 看黄色毛片网站| 91久久精品电影网| 在线免费观看的www视频| 少妇猛男粗大的猛烈进出视频 | 国产久久久一区二区三区| 日韩国内少妇激情av| 22中文网久久字幕| 国产一区亚洲一区在线观看| 一区二区三区四区激情视频| 久久久久精品性色| 成人综合一区亚洲| 两个人的视频大全免费| 高清日韩中文字幕在线| 一级毛片我不卡| 直男gayav资源| 夜夜看夜夜爽夜夜摸| 亚洲av成人精品一二三区| 嫩草影院入口| 日日啪夜夜爽| 国产白丝娇喘喷水9色精品| 黄片wwwwww| 日韩欧美精品免费久久| 丝瓜视频免费看黄片| 少妇熟女aⅴ在线视频| 最新中文字幕久久久久| 女人十人毛片免费观看3o分钟| 嘟嘟电影网在线观看| 白带黄色成豆腐渣| 国产精品av视频在线免费观看| 午夜老司机福利剧场| 国产大屁股一区二区在线视频| 国产老妇女一区| 国产精品一及| 欧美成人一区二区免费高清观看| 日韩欧美国产在线观看| 精品国产一区二区三区久久久樱花 | 哪个播放器可以免费观看大片| 色哟哟·www| 午夜爱爱视频在线播放| 免费av观看视频| 亚洲av中文av极速乱| 亚洲av电影在线观看一区二区三区 | 卡戴珊不雅视频在线播放| 久久精品国产亚洲av天美| 国产高清不卡午夜福利| 午夜激情久久久久久久| 51国产日韩欧美| 免费少妇av软件| 美女脱内裤让男人舔精品视频| 哪个播放器可以免费观看大片| 欧美激情在线99| 狠狠精品人妻久久久久久综合| 国产精品.久久久| 少妇丰满av| 美女cb高潮喷水在线观看| 国产高清三级在线| 日本黄大片高清| 色综合站精品国产| 超碰97精品在线观看| 九草在线视频观看| 搡女人真爽免费视频火全软件| 午夜福利视频精品| av线在线观看网站| 日本免费在线观看一区| 丰满人妻一区二区三区视频av| av女优亚洲男人天堂| a级一级毛片免费在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品色激情综合| 大香蕉久久网| 街头女战士在线观看网站| 欧美日韩视频高清一区二区三区二| 亚洲人与动物交配视频| 欧美+日韩+精品| 男女下面进入的视频免费午夜| 亚洲va在线va天堂va国产| 免费黄频网站在线观看国产| 精品一区二区三区人妻视频| 亚洲欧洲日产国产| 成人鲁丝片一二三区免费| 黄色配什么色好看| 婷婷六月久久综合丁香| 国产老妇伦熟女老妇高清| av专区在线播放| av国产免费在线观看| 亚洲高清免费不卡视频| 中文字幕人妻熟人妻熟丝袜美| 人体艺术视频欧美日本| 国产 亚洲一区二区三区 | 午夜福利在线观看吧| 欧美成人午夜免费资源| 国产精品1区2区在线观看.| 亚洲av男天堂| 婷婷色综合www| 久久这里有精品视频免费| 汤姆久久久久久久影院中文字幕 | 久热久热在线精品观看| 亚洲成人av在线免费| 国产精品国产三级国产av玫瑰| 亚洲国产色片| 99久久精品一区二区三区| 国产精品福利在线免费观看| 91在线精品国自产拍蜜月| 日日啪夜夜爽| 国内少妇人妻偷人精品xxx网站| 嫩草影院精品99| 国产精品一二三区在线看| 国产男女超爽视频在线观看| 久久人人爽人人爽人人片va| 亚洲精品成人av观看孕妇| 亚洲精品一二三| 免费播放大片免费观看视频在线观看| 简卡轻食公司| 欧美不卡视频在线免费观看| 亚洲美女视频黄频| 亚洲国产色片| 国内精品一区二区在线观看| 国产成人91sexporn| 三级国产精品片| 免费观看av网站的网址| 成人漫画全彩无遮挡| 女的被弄到高潮叫床怎么办| 夜夜看夜夜爽夜夜摸| 一级毛片电影观看| 欧美极品一区二区三区四区| 肉色欧美久久久久久久蜜桃 | 欧美日韩视频高清一区二区三区二| 日韩欧美国产在线观看| 国产精品综合久久久久久久免费| 天美传媒精品一区二区| 蜜桃久久精品国产亚洲av| 日韩欧美 国产精品| 少妇裸体淫交视频免费看高清| 国产乱人偷精品视频| 日韩av不卡免费在线播放| 视频中文字幕在线观看| 国产在线一区二区三区精| 纵有疾风起免费观看全集完整版 | 国产91av在线免费观看| 亚洲欧美精品自产自拍| 视频中文字幕在线观看| 老女人水多毛片| 精品人妻偷拍中文字幕| 校园人妻丝袜中文字幕| 高清在线视频一区二区三区| 亚洲av国产av综合av卡| 久久综合国产亚洲精品| 亚洲人成网站高清观看| 久久国产乱子免费精品| 18禁在线无遮挡免费观看视频| 黄色欧美视频在线观看| 国产精品无大码| 国语对白做爰xxxⅹ性视频网站| 精品久久久久久久人妻蜜臀av| 欧美bdsm另类| 久久久色成人| 精品久久久久久久人妻蜜臀av| 成人国产麻豆网| 简卡轻食公司| 免费大片黄手机在线观看| 床上黄色一级片| 欧美3d第一页| 免费观看av网站的网址| 视频中文字幕在线观看| 99热这里只有是精品50| 男人舔奶头视频| 一二三四中文在线观看免费高清| 亚洲在线观看片| 久久精品国产自在天天线| 2021少妇久久久久久久久久久| 国产探花极品一区二区| 国产午夜精品久久久久久一区二区三区| 午夜福利在线观看吧| 国产又色又爽无遮挡免| 久久久久久伊人网av| 国产亚洲一区二区精品| 人体艺术视频欧美日本| 国产精品精品国产色婷婷| 亚洲av在线观看美女高潮| 日日啪夜夜撸| 禁无遮挡网站| 岛国毛片在线播放| 亚洲国产欧美在线一区| 亚洲成人一二三区av| 在线天堂最新版资源| 国产综合精华液| 啦啦啦韩国在线观看视频| 精品国产三级普通话版| 99热全是精品| 欧美日韩一区二区视频在线观看视频在线 | 极品少妇高潮喷水抽搐| 国产爱豆传媒在线观看| 69av精品久久久久久| 国产在线男女| 精品久久久久久电影网| 麻豆精品久久久久久蜜桃| 夜夜看夜夜爽夜夜摸| 日本av手机在线免费观看| 国产色爽女视频免费观看| 免费无遮挡裸体视频| 亚洲自拍偷在线| 精品欧美国产一区二区三| 看免费成人av毛片| 2018国产大陆天天弄谢| 国产精品蜜桃在线观看| 免费大片18禁| 最近最新中文字幕大全电影3| 免费高清在线观看视频在线观看| 久久97久久精品| 欧美xxxx性猛交bbbb| 亚洲精品第二区| 成年免费大片在线观看| 蜜桃亚洲精品一区二区三区| 日本一二三区视频观看| 国产免费视频播放在线视频 | 97超碰精品成人国产| 久久鲁丝午夜福利片| 精品少妇黑人巨大在线播放| 精品久久久久久久久亚洲| 亚洲国产最新在线播放| 看黄色毛片网站| 97人妻精品一区二区三区麻豆| 男人爽女人下面视频在线观看| 波多野结衣巨乳人妻| 久久精品国产自在天天线| 99热全是精品| 永久网站在线| 亚洲精品影视一区二区三区av| 国产精品不卡视频一区二区| 九色成人免费人妻av| 国产91av在线免费观看| 神马国产精品三级电影在线观看| 丰满少妇做爰视频| 中文资源天堂在线| 黑人高潮一二区| 免费看日本二区| 亚洲av国产av综合av卡| videossex国产| 亚洲最大成人中文| 床上黄色一级片| 国产91av在线免费观看| 国内少妇人妻偷人精品xxx网站| videossex国产| 国产探花极品一区二区| 精品久久久久久久久av| 美女高潮的动态| 久久久亚洲精品成人影院| 精品久久久精品久久久| 久久亚洲国产成人精品v| 一级片'在线观看视频| 中文字幕制服av| 大香蕉久久成人网| 亚洲欧美中文字幕日韩二区| 女性生殖器流出的白浆| 两个人看的免费小视频| 日韩电影二区| 99热网站在线观看| 麻豆乱淫一区二区| 欧美精品一区二区大全| 人妻人人澡人人爽人人| 伦理电影大哥的女人| 久久久久久久久免费视频了| av视频免费观看在线观看| 女人精品久久久久毛片| 自线自在国产av| 国产精品成人在线| 9热在线视频观看99| 国产av码专区亚洲av| 卡戴珊不雅视频在线播放| 欧美+日韩+精品| 搡老乐熟女国产| 欧美在线黄色| 日韩一卡2卡3卡4卡2021年| 激情视频va一区二区三区| 这个男人来自地球电影免费观看 | 亚洲欧美中文字幕日韩二区| 国产男女超爽视频在线观看| 日本爱情动作片www.在线观看| 亚洲第一青青草原| 精品一区在线观看国产| 久久精品熟女亚洲av麻豆精品| 日本vs欧美在线观看视频| 18+在线观看网站| 青草久久国产| 久久久久国产精品人妻一区二区| 一级毛片黄色毛片免费观看视频| 日本爱情动作片www.在线观看| 蜜桃在线观看..| 亚洲综合色网址| 在线观看一区二区三区激情| 菩萨蛮人人尽说江南好唐韦庄| 丝袜在线中文字幕| 欧美亚洲 丝袜 人妻 在线| 国产精品99久久99久久久不卡 | 亚洲av欧美aⅴ国产| 超碰成人久久| 美女脱内裤让男人舔精品视频| 大片免费播放器 马上看| 美女高潮到喷水免费观看| 天美传媒精品一区二区| 免费在线观看完整版高清| 亚洲欧洲日产国产| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆 | 国产精品香港三级国产av潘金莲 | 亚洲精品美女久久久久99蜜臀 | 精品国产露脸久久av麻豆| 日日摸夜夜添夜夜爱| 在线精品无人区一区二区三| av一本久久久久| 国产av一区二区精品久久| 国产激情久久老熟女| av视频免费观看在线观看| 国产男女内射视频| 天天影视国产精品| 18禁裸乳无遮挡动漫免费视频| 九九爱精品视频在线观看| 日韩中字成人| 中文字幕制服av| 国产熟女午夜一区二区三区| 热re99久久精品国产66热6| 美女大奶头黄色视频| 亚洲一区中文字幕在线| 大话2 男鬼变身卡| www日本在线高清视频| 国产熟女欧美一区二区| 午夜91福利影院| 亚洲成色77777| 美女福利国产在线| 亚洲中文av在线| 午夜福利视频精品| 亚洲中文av在线| 卡戴珊不雅视频在线播放| 岛国毛片在线播放| 亚洲 欧美一区二区三区| 免费大片黄手机在线观看| 久久97久久精品| 国产成人精品福利久久| 日本爱情动作片www.在线观看| 亚洲成av片中文字幕在线观看 | 成人二区视频| 男人操女人黄网站| 国产成人欧美| 欧美日韩精品成人综合77777| 国产成人欧美| 69精品国产乱码久久久| 人妻系列 视频| 街头女战士在线观看网站| 热re99久久精品国产66热6| 王馨瑶露胸无遮挡在线观看| 亚洲久久久国产精品| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品999| 国产一区二区 视频在线| 久久久久精品人妻al黑| 观看美女的网站| 精品少妇黑人巨大在线播放| 日本av手机在线免费观看| 九草在线视频观看| 在线观看免费高清a一片| 国产成人精品在线电影| 丝袜美足系列| 国产高清不卡午夜福利| 亚洲欧美一区二区三区久久| 中文字幕人妻丝袜制服| 色婷婷久久久亚洲欧美| 亚洲国产精品一区二区三区在线| 免费在线观看视频国产中文字幕亚洲 | av电影中文网址| 99国产综合亚洲精品| 看免费av毛片| 久久久欧美国产精品| 青青草视频在线视频观看| 免费大片黄手机在线观看| 91精品国产国语对白视频| 久久久久久久久久久免费av| 精品卡一卡二卡四卡免费| 免费观看av网站的网址| 曰老女人黄片| 亚洲精品日韩在线中文字幕| 韩国高清视频一区二区三区| 一区在线观看完整版| 成年av动漫网址| 亚洲伊人久久精品综合| 国产日韩欧美亚洲二区| 啦啦啦在线观看免费高清www| 一区二区日韩欧美中文字幕| av一本久久久久| 亚洲,欧美精品.| 亚洲 欧美一区二区三区| 女性生殖器流出的白浆| 最近手机中文字幕大全| 精品福利永久在线观看| 亚洲精品第二区| 亚洲熟女精品中文字幕| 99久久人妻综合| 国产成人精品婷婷| 欧美老熟妇乱子伦牲交| 91在线精品国自产拍蜜月| 成人亚洲欧美一区二区av| 热99国产精品久久久久久7| 久久久亚洲精品成人影院| 妹子高潮喷水视频| 搡女人真爽免费视频火全软件| 亚洲国产毛片av蜜桃av| 国产亚洲午夜精品一区二区久久| 免费在线观看完整版高清| 91精品国产国语对白视频| 叶爱在线成人免费视频播放| 热99国产精品久久久久久7| 亚洲av国产av综合av卡| 美女视频免费永久观看网站| 婷婷色综合大香蕉| 国产爽快片一区二区三区| 18禁裸乳无遮挡动漫免费视频| 久久精品国产鲁丝片午夜精品| 在线观看三级黄色| 国产伦理片在线播放av一区| 美女福利国产在线| 免费观看在线日韩| 欧美成人午夜精品| 视频区图区小说| h视频一区二区三区| 久久人人爽av亚洲精品天堂| 久久国产精品大桥未久av| 午夜激情av网站| 日韩中字成人| 国产精品久久久av美女十八| 成人影院久久|