• <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级片在线免费高清观看视频| 国产欧美日韩一区二区三区在线 | 久久国产精品男人的天堂亚洲 | 欧美成人午夜免费资源| 美女大奶头黄色视频| 91午夜精品亚洲一区二区三区| 2022亚洲国产成人精品| 亚洲精品视频女| 最新的欧美精品一区二区| 久久99蜜桃精品久久| 国产熟女午夜一区二区三区 | 人人妻人人添人人爽欧美一区卜| 99九九线精品视频在线观看视频| 丝瓜视频免费看黄片| 嫩草影院入口| 夫妻午夜视频| 91精品伊人久久大香线蕉| 午夜91福利影院| 黄色视频在线播放观看不卡| 国产欧美日韩一区二区三区在线 | 国产有黄有色有爽视频| .国产精品久久| 秋霞伦理黄片| 日日摸夜夜添夜夜添av毛片| 国产免费一区二区三区四区乱码| 国产成人一区二区在线| 男女无遮挡免费网站观看| av视频免费观看在线观看| 国产又色又爽无遮挡免| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av天美| 欧美成人精品欧美一级黄| 日韩强制内射视频| 国产高清不卡午夜福利| 九九在线视频观看精品| 国产成人精品在线电影| 99视频精品全部免费 在线| 91久久精品电影网| 天天躁夜夜躁狠狠久久av| av网站免费在线观看视频| 看非洲黑人一级黄片| 最新中文字幕久久久久| 国产亚洲午夜精品一区二区久久| 日韩强制内射视频| 国产黄片视频在线免费观看| 亚洲高清免费不卡视频| 男男h啪啪无遮挡| 久久青草综合色| 免费av中文字幕在线| 欧美少妇被猛烈插入视频| 午夜免费鲁丝| 精品少妇内射三级| 亚洲高清免费不卡视频| 亚洲国产精品成人久久小说| 一边摸一边做爽爽视频免费| 婷婷成人精品国产| 精品人妻一区二区三区麻豆| 日日啪夜夜爽| 免费黄网站久久成人精品| 国产欧美日韩一区二区三区在线 | 国产黄色免费在线视频| 日本免费在线观看一区| 国产在线一区二区三区精| 精品人妻在线不人妻| 黑丝袜美女国产一区| 亚洲av中文av极速乱| 久久久久久久久大av| 爱豆传媒免费全集在线观看| 女人精品久久久久毛片| 国产精品一区二区在线不卡| 午夜老司机福利剧场| 久久精品国产亚洲网站| 亚洲国产欧美日韩在线播放| 亚洲少妇的诱惑av| 男人添女人高潮全过程视频| 午夜视频国产福利| 少妇的逼水好多| 国产精品.久久久| 老司机影院毛片| 日本黄大片高清| 观看av在线不卡| 精品午夜福利在线看| 青春草亚洲视频在线观看| av在线播放精品| 中文字幕精品免费在线观看视频 | 久久久久视频综合| 欧美精品亚洲一区二区| 狂野欧美激情性bbbbbb| 交换朋友夫妻互换小说| 日日摸夜夜添夜夜爱| 亚洲av免费高清在线观看| 国产乱来视频区| 亚洲av不卡在线观看| 国产成人精品一,二区| 亚洲久久久国产精品| 国产av精品麻豆| 国产无遮挡羞羞视频在线观看| 校园人妻丝袜中文字幕| 国语对白做爰xxxⅹ性视频网站| 少妇 在线观看| 超碰97精品在线观看| 亚洲三级黄色毛片| 黄色欧美视频在线观看| 热99久久久久精品小说推荐| 在线播放无遮挡| 欧美日本中文国产一区发布| 国产欧美亚洲国产| 大香蕉久久网| 国产永久视频网站| 日韩精品有码人妻一区| 五月玫瑰六月丁香| 国产综合精华液| 女性被躁到高潮视频| 日韩伦理黄色片| videosex国产| 母亲3免费完整高清在线观看 | 一区二区三区乱码不卡18| 自线自在国产av| 黄色配什么色好看| 欧美日韩视频精品一区| 久久午夜综合久久蜜桃| 九色亚洲精品在线播放| 建设人人有责人人尽责人人享有的| av黄色大香蕉| 老熟女久久久| 精品久久久久久久久av| 日韩欧美精品免费久久| 三上悠亚av全集在线观看| av不卡在线播放| 丝袜美足系列| 国产精品久久久久成人av| 在线 av 中文字幕| 九色亚洲精品在线播放| 纯流量卡能插随身wifi吗| 国产午夜精品一二区理论片| 成人亚洲精品一区在线观看| 免费高清在线观看视频在线观看| 日韩三级伦理在线观看| 精品久久久精品久久久| 在线观看www视频免费| 国产极品粉嫩免费观看在线 | 精品午夜福利在线看| 久久青草综合色| 亚洲国产最新在线播放| 一边摸一边做爽爽视频免费| 熟妇人妻不卡中文字幕| 人妻一区二区av| 精品熟女少妇av免费看| 国产黄色免费在线视频| 菩萨蛮人人尽说江南好唐韦庄| 观看av在线不卡| 91精品一卡2卡3卡4卡| 欧美激情极品国产一区二区三区 | 日日摸夜夜添夜夜添av毛片| 亚洲久久久国产精品| 久久久久久久久久久久大奶| 99精国产麻豆久久婷婷| 午夜日本视频在线| 国产成人精品婷婷| 亚洲欧洲日产国产| 国产无遮挡羞羞视频在线观看| 午夜福利视频在线观看免费| 久久久国产精品麻豆| 国产有黄有色有爽视频| av线在线观看网站| 欧美精品高潮呻吟av久久| 国产精品.久久久| 精品国产一区二区久久| 91在线精品国自产拍蜜月| 亚洲国产精品成人久久小说| 99久久人妻综合| 亚洲无线观看免费| av网站免费在线观看视频| 亚洲美女视频黄频| 美女福利国产在线| 麻豆乱淫一区二区| 国产午夜精品一二区理论片| 午夜福利在线观看免费完整高清在| 人人妻人人澡人人看| 国产国拍精品亚洲av在线观看| 亚洲婷婷狠狠爱综合网| 自线自在国产av| 国产在线视频一区二区| 人成视频在线观看免费观看| 欧美日韩精品成人综合77777| 色哟哟·www| 免费观看的影片在线观看| 亚洲av在线观看美女高潮| 高清视频免费观看一区二区| 精品熟女少妇av免费看| 草草在线视频免费看| 国产极品粉嫩免费观看在线 | 国语对白做爰xxxⅹ性视频网站| 十八禁高潮呻吟视频| 麻豆成人av视频| 亚洲av欧美aⅴ国产| 日韩大片免费观看网站| 91久久精品国产一区二区成人| 乱人伦中国视频| 久久鲁丝午夜福利片| a级毛片在线看网站| 久久人人爽人人爽人人片va| 精品卡一卡二卡四卡免费| 亚洲国产精品成人久久小说| 乱人伦中国视频| 赤兔流量卡办理| 国产高清不卡午夜福利| 全区人妻精品视频| a级毛片在线看网站| 国产精品三级大全| 久久久久网色| 亚洲精品456在线播放app| 人妻 亚洲 视频| 母亲3免费完整高清在线观看 | 麻豆乱淫一区二区| 涩涩av久久男人的天堂| 久久精品人人爽人人爽视色| 街头女战士在线观看网站| 日本猛色少妇xxxxx猛交久久| av不卡在线播放| 在线观看一区二区三区激情| 天堂中文最新版在线下载| 国产综合精华液| 国产成人精品一,二区| 日本色播在线视频| 久久精品国产亚洲网站| 男的添女的下面高潮视频| 黄色毛片三级朝国网站| 色5月婷婷丁香| 亚洲精品乱久久久久久| 国产一级毛片在线| 国产极品天堂在线| 色婷婷久久久亚洲欧美| 插逼视频在线观看| 中文字幕av电影在线播放| 久热这里只有精品99| 国产国语露脸激情在线看| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人 | 狂野欧美激情性xxxx在线观看| 精品亚洲乱码少妇综合久久| 看非洲黑人一级黄片| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 色视频在线一区二区三区| 又黄又爽又刺激的免费视频.| 欧美日韩视频精品一区| 国产日韩欧美视频二区| 永久网站在线| 精品久久久精品久久久| 久久久久久久久久久免费av| 91在线精品国自产拍蜜月| 久热这里只有精品99| 亚洲综合色惰| 高清毛片免费看| 2022亚洲国产成人精品| 午夜福利影视在线免费观看| 国产在线免费精品| 成人漫画全彩无遮挡| 国产免费又黄又爽又色| 最新中文字幕久久久久| av电影中文网址| 国产成人freesex在线| 精品久久国产蜜桃| 又黄又爽又刺激的免费视频.| 久久精品国产自在天天线| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| 91精品国产国语对白视频| 日韩一区二区视频免费看| 少妇人妻 视频| 久久久精品区二区三区| 精品熟女少妇av免费看| 午夜福利影视在线免费观看| 免费高清在线观看日韩| 男女高潮啪啪啪动态图| 美女中出高潮动态图| 国产伦精品一区二区三区视频9| 久久精品国产亚洲网站| 视频在线观看一区二区三区| 美女国产高潮福利片在线看| 免费高清在线观看日韩| 国产精品久久久久久av不卡| 国产一区二区三区综合在线观看 | 三级国产精品片| 男女免费视频国产| 美女cb高潮喷水在线观看| 精品酒店卫生间| 久久影院123| 亚洲欧美色中文字幕在线| 又大又黄又爽视频免费| 18禁动态无遮挡网站| 亚洲精品亚洲一区二区| 狂野欧美白嫩少妇大欣赏| 久久人妻熟女aⅴ| 美女内射精品一级片tv| 亚洲av二区三区四区| 视频中文字幕在线观看| 国产免费一区二区三区四区乱码| 国产在线视频一区二区| 黄色怎么调成土黄色| 18+在线观看网站| 国产亚洲欧美精品永久| 精品一区二区三卡| 精品久久蜜臀av无| tube8黄色片| 亚洲国产av新网站| 黑人巨大精品欧美一区二区蜜桃 | av国产精品久久久久影院| 国产又色又爽无遮挡免| 伦理电影免费视频| 国产成人午夜福利电影在线观看| 男女国产视频网站| 国产精品人妻久久久久久| 亚洲av男天堂| 日本91视频免费播放| 99热国产这里只有精品6| 久久久精品区二区三区| 国产午夜精品久久久久久一区二区三区| 简卡轻食公司| 久久热精品热| 简卡轻食公司| 久久久亚洲精品成人影院| 在现免费观看毛片| 永久免费av网站大全| 99久久精品国产国产毛片| 黄色一级大片看看| 高清不卡的av网站| 大香蕉久久网| 亚洲成人手机| 国产精品一区二区在线观看99| 免费观看性生交大片5| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久久久久久人人人人人人| 国产免费视频播放在线视频| 少妇高潮的动态图| 国产成人a∨麻豆精品| 看免费成人av毛片| 卡戴珊不雅视频在线播放| 99国产精品免费福利视频| 久久影院123| 国产又色又爽无遮挡免| 制服人妻中文乱码| 熟女av电影| 国产高清国产精品国产三级| 熟女av电影| 成人影院久久| 国产伦理片在线播放av一区| 综合色丁香网| 国产在线免费精品| 亚洲高清免费不卡视频| 插逼视频在线观看| 免费日韩欧美在线观看| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 久久人妻熟女aⅴ| 一本大道久久a久久精品| 欧美日韩综合久久久久久| 久久午夜福利片| 大香蕉久久网| 精品国产国语对白av| 色94色欧美一区二区| 男男h啪啪无遮挡| 久久精品人人爽人人爽视色| 一区二区三区乱码不卡18| 欧美日韩综合久久久久久| 久久午夜福利片| 国产高清不卡午夜福利| 国产精品久久久久成人av| 免费观看在线日韩| 欧美日韩av久久| 欧美精品一区二区大全| 最近2019中文字幕mv第一页| 汤姆久久久久久久影院中文字幕| 欧美丝袜亚洲另类| 欧美日韩精品成人综合77777| 日韩欧美一区视频在线观看| 一级毛片aaaaaa免费看小| 成年人午夜在线观看视频| 热99国产精品久久久久久7| 亚洲av在线观看美女高潮| 夜夜看夜夜爽夜夜摸| 久久久久人妻精品一区果冻| 精品亚洲成国产av| 精品人妻偷拍中文字幕| 婷婷色综合大香蕉| 亚洲成人手机| 欧美激情 高清一区二区三区| 69精品国产乱码久久久| 免费日韩欧美在线观看| 免费看光身美女| 激情五月婷婷亚洲| 99热这里只有精品一区| 大话2 男鬼变身卡| 欧美日韩视频高清一区二区三区二| 欧美成人午夜免费资源| av在线观看视频网站免费| 精品久久久精品久久久| 国产成人一区二区在线| 女人精品久久久久毛片| 国产无遮挡羞羞视频在线观看| 人体艺术视频欧美日本| 亚洲精品国产色婷婷电影| 亚洲精品国产av蜜桃| 日本-黄色视频高清免费观看| 亚洲欧洲日产国产| 新久久久久国产一级毛片| 免费观看a级毛片全部| 最近的中文字幕免费完整| 日本av手机在线免费观看| 最近2019中文字幕mv第一页| .国产精品久久| 成人影院久久| 乱码一卡2卡4卡精品| 最黄视频免费看| 99热6这里只有精品| 国产精品久久久久久久电影| 黄片播放在线免费| 久久婷婷青草| 亚洲精品色激情综合| 欧美成人午夜免费资源| 一区在线观看完整版| 韩国高清视频一区二区三区| av网站免费在线观看视频| 亚洲四区av| 五月开心婷婷网| 9色porny在线观看| 久久久精品免费免费高清| 少妇人妻 视频| 香蕉精品网在线| 亚洲精品亚洲一区二区| av电影中文网址| .国产精品久久| 精品少妇内射三级| 黄色配什么色好看| 亚洲情色 制服丝袜| 精品一区二区三卡| 成人国语在线视频| 久久精品国产a三级三级三级| 欧美日韩成人在线一区二区| 一本大道久久a久久精品| 少妇高潮的动态图| 在线观看国产h片| 九九爱精品视频在线观看| 国产永久视频网站| 久久精品久久精品一区二区三区| 亚洲无线观看免费| 免费久久久久久久精品成人欧美视频 | 18+在线观看网站| 777米奇影视久久| 亚洲精品一区蜜桃| 成人国产麻豆网| 午夜激情福利司机影院| 边亲边吃奶的免费视频| 国内精品宾馆在线| 久久人人爽人人片av| 国产老妇伦熟女老妇高清| 色5月婷婷丁香| 国产精品人妻久久久影院| 欧美日韩av久久| 十八禁网站网址无遮挡| 日本猛色少妇xxxxx猛交久久| 欧美成人精品欧美一级黄| 精品人妻偷拍中文字幕| 美女cb高潮喷水在线观看| 黑人欧美特级aaaaaa片| 又大又黄又爽视频免费| 精品一区二区三区视频在线| 国产熟女欧美一区二区| 在线观看美女被高潮喷水网站| 亚洲精华国产精华液的使用体验| 久久久久视频综合| 午夜激情福利司机影院| 免费观看在线日韩| 卡戴珊不雅视频在线播放| 亚洲性久久影院| 日韩三级伦理在线观看| 精品人妻一区二区三区麻豆| 2022亚洲国产成人精品| 校园人妻丝袜中文字幕| a级毛片黄视频| 全区人妻精品视频| 午夜日本视频在线| av卡一久久| 午夜福利视频精品| 欧美激情 高清一区二区三区| 国产视频首页在线观看| 丰满饥渴人妻一区二区三| 韩国av在线不卡| av.在线天堂| 观看美女的网站| 免费看不卡的av| 中文精品一卡2卡3卡4更新| 免费看av在线观看网站| 我的老师免费观看完整版| 一区二区三区免费毛片| 十八禁网站网址无遮挡| 激情五月婷婷亚洲| 午夜影院在线不卡| 日韩大片免费观看网站| 婷婷成人精品国产| 欧美亚洲 丝袜 人妻 在线| 丝袜脚勾引网站| 国产av码专区亚洲av| 国产成人a∨麻豆精品| 亚洲精品美女久久av网站| 久久影院123| 精品人妻在线不人妻| av又黄又爽大尺度在线免费看| 欧美成人精品欧美一级黄| av福利片在线| 久久久久久久国产电影| 成人国产麻豆网| 男人添女人高潮全过程视频| 黄色怎么调成土黄色| 国产av码专区亚洲av| 久久精品国产亚洲网站| 久久精品国产自在天天线| 免费少妇av软件| 久久综合国产亚洲精品| 老司机影院毛片| 精品久久久精品久久久| 伦理电影大哥的女人| 国精品久久久久久国模美| 欧美激情国产日韩精品一区| 观看美女的网站| 麻豆成人av视频| 丰满乱子伦码专区| 男女无遮挡免费网站观看| 一边亲一边摸免费视频| 午夜影院在线不卡| 在线观看三级黄色| 精品一区二区三卡| 大片电影免费在线观看免费| 午夜影院在线不卡| 久久久国产一区二区| 久久国产精品男人的天堂亚洲 | 亚洲精品亚洲一区二区| 久久久久久久久大av| 亚洲精品久久午夜乱码| 欧美+日韩+精品| 亚洲精品久久午夜乱码| 一个人看视频在线观看www免费| 国产精品女同一区二区软件| 日韩电影二区| 国产男人的电影天堂91| 一级毛片aaaaaa免费看小| 国产片内射在线| 久久精品国产亚洲av涩爱| 美女福利国产在线| 一级毛片黄色毛片免费观看视频| 国模一区二区三区四区视频| 国产av码专区亚洲av| 九色亚洲精品在线播放| 亚洲美女黄色视频免费看| 中文字幕亚洲精品专区| 亚洲国产av影院在线观看| 国产精品麻豆人妻色哟哟久久| 2018国产大陆天天弄谢| 日韩在线高清观看一区二区三区| 久久99热这里只频精品6学生| www.av在线官网国产| 免费大片18禁| 十八禁网站网址无遮挡| 日日摸夜夜添夜夜爱| 夫妻性生交免费视频一级片| 丰满少妇做爰视频| 国产免费一级a男人的天堂| 亚洲国产精品成人久久小说| 男女边摸边吃奶| 91精品国产九色| 国产爽快片一区二区三区| 国产熟女午夜一区二区三区 | 婷婷色麻豆天堂久久| 欧美成人精品欧美一级黄| 午夜av观看不卡| 久久精品国产a三级三级三级| 中文精品一卡2卡3卡4更新| 色5月婷婷丁香| 亚洲婷婷狠狠爱综合网| 国产永久视频网站| 国产69精品久久久久777片| 欧美bdsm另类| 天天影视国产精品| 久久久久人妻精品一区果冻| 国产精品.久久久| 在线 av 中文字幕| 国产午夜精品久久久久久一区二区三区| 狠狠精品人妻久久久久久综合| 最近最新中文字幕免费大全7| h视频一区二区三区| 国产精品女同一区二区软件| 黑丝袜美女国产一区| a 毛片基地| 哪个播放器可以免费观看大片| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲欧美精品永久| 久久亚洲国产成人精品v| 亚洲成人手机| 高清黄色对白视频在线免费看| 免费看av在线观看网站| 亚洲综合色惰| 亚洲精品久久成人aⅴ小说 | 蜜桃国产av成人99| 精品久久久久久电影网| 亚洲三级黄色毛片| 国产日韩一区二区三区精品不卡 | 亚洲国产av新网站| 一级毛片aaaaaa免费看小| 精品少妇久久久久久888优播| 免费av不卡在线播放| 亚洲国产色片| 18在线观看网站| 午夜精品国产一区二区电影|