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

    柔性壓力管道甩擊特性數(shù)值模擬

    2022-09-14 08:37:14王曉川李月琴趙文勝向美景
    重慶大學(xué)學(xué)報(bào) 2022年8期
    關(guān)鍵詞:柔性流速流體

    王曉川,劉 錚,李月琴,王 滿(mǎn),趙文勝,向美景

    (1.武漢大學(xué) a. 動(dòng)力與機(jī)械學(xué)院;b. 水射流理論與新技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430072;2.煉焦煤資源開(kāi)發(fā)及綜合利用國(guó)家重點(diǎn)實(shí)驗(yàn)室,河南 平頂山 467000;3.中國(guó)平煤神馬能源化工集團(tuán)有限責(zé)任公司,河南 平頂山 467000)

    管道作為一種重要的工業(yè)物料傳輸工具,承擔(dān)著能源輸送、動(dòng)力傳輸、物質(zhì)循環(huán)的重大任務(wù),在國(guó)民生產(chǎn)中發(fā)揮著重要作用,被譽(yù)為工業(yè)生產(chǎn)的動(dòng)脈。柔性壓力管道作為一種重要的流體輸送工具,以適用性強(qiáng)、安裝快捷等特點(diǎn)廣泛應(yīng)用于臨時(shí)或機(jī)動(dòng)工程中,如工程機(jī)械、煤礦工業(yè)、冶金鍛壓、各種機(jī)床以及各工業(yè)部門(mén)機(jī)械化、自動(dòng)化液壓系統(tǒng)中輸送具有一定壓力和溫度的石油基液體、氣體等。當(dāng)柔性壓力管道脫扣斷裂時(shí),管內(nèi)高壓流體噴射而出,管道受到流體的噴射力作用發(fā)生變形,管道的變形也引起流場(chǎng)變化,這種流體和固體間的耦合作用,使管道持續(xù)發(fā)生甩動(dòng),對(duì)周?chē)嘶蛭锏陌踩斐蓸O大的威脅。柔性壓力管道工作壓力高,甩擊事故時(shí)常發(fā)生。當(dāng)應(yīng)用于高壓射流技術(shù)時(shí),其工作壓力普遍在20 MPa以上,在采礦、非常規(guī)油氣開(kāi)采等特定領(lǐng)域中,超高壓水射流工作壓力已達(dá)400 MPa以上。由于管道自身韌性小,管道振甩波及范圍大,在同一尺度下的甩擊破壞性強(qiáng),防控難度更大。柔性壓力管道在管道破斷泄漏過(guò)程中涉及的耦合參數(shù)較多,加上柔性壓力管道在鋪設(shè)過(guò)程中比較隨意,延伸過(guò)程中管道起伏不定情況更為復(fù)雜。目前有關(guān)柔性壓力管道甩擊運(yùn)動(dòng)的特點(diǎn)和防護(hù)理論成果較少,隨著柔性壓力管道的應(yīng)用范圍加大及工作壓力等級(jí)的提高,這項(xiàng)工作變得尤為迫切,特別是礦業(yè)和石油行業(yè)臨時(shí)性使用的柔性壓力管道數(shù)量極多,因此很有必要對(duì)柔性壓力管道的甩擊運(yùn)動(dòng)特性進(jìn)行系統(tǒng)研究,提高我國(guó)的壓力管道安全管理水平。

    柔性壓力管道的甩擊運(yùn)動(dòng)是一種高度非線性化的運(yùn)動(dòng),本質(zhì)上是一個(gè)大變形瞬態(tài)動(dòng)力學(xué)過(guò)程。幾何非線性、材料非線性和接觸非線性都對(duì)管道的甩擊運(yùn)動(dòng)有非常顯著的影響,在分析管道甩擊運(yùn)動(dòng)的過(guò)程中,很多理論都難以同時(shí)考慮到以上幾種情況。隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展,利用數(shù)值模擬的方法對(duì)柔性管道的甩擊行為研究得以實(shí)現(xiàn),給該問(wèn)題帶來(lái)了新的突破方向。在20世紀(jì)70年代和80年代,意大利Dini、Lazzeri等先后發(fā)展了FRUSTRA、PAULA82等有限元軟件,美國(guó)通用電氣也開(kāi)發(fā)了三維非線性動(dòng)力分析程序COMET以及ADINA、ABAQUS等大型通用有限元分析軟件,這些軟件對(duì)簡(jiǎn)化過(guò)后的數(shù)值模型進(jìn)行分析求解。該階段主要采用梁?jiǎn)卧M管道,用彈簧單元模擬限制件,在理論上尚存在諸多不完善的地方。Dini和Lazzeri[1-3]隨后改進(jìn)了梁?jiǎn)卧P?,通過(guò)建立管道扁平化下“彎矩-曲率”的非線性關(guān)系,使其更符合管道甩擊過(guò)程中的變形。Ueda等[4]則在對(duì)管道彎頭碰撞區(qū)進(jìn)行靜力彈塑性分析嘗試使用殼單元模型,取得了不錯(cuò)的效果。Ma等[5]使用ADINA分析程序?qū)艿琅c限制件之間碰撞接觸進(jìn)行了數(shù)值模擬計(jì)算。1997年,Olson等[6]對(duì)直管與管內(nèi)流體的耦合振動(dòng)變形的有限元法進(jìn)行了綜述,系統(tǒng)地總結(jié)了計(jì)算壓力管道流固耦合問(wèn)題的基本思想。Potapov等[7]利用一維、三維混合管道模型對(duì)甩擊運(yùn)動(dòng)的變形和斷裂情況進(jìn)行了模擬計(jì)算,取得了比簡(jiǎn)化的管道模型更好的結(jié)果。劉德輔等[8]利用有限單元法推導(dǎo)在外界環(huán)境載荷條件下管道的側(cè)向振動(dòng)微分方程,研究隨機(jī)波浪對(duì)管道振動(dòng)的影響。張艷萍等[9]采用結(jié)構(gòu)有限元方法和CFD方法實(shí)現(xiàn)管道與流體間的雙向耦合,重點(diǎn)研究了變化流速下管道結(jié)構(gòu)和流體的瞬態(tài)響應(yīng)。Reid等[10]利用小直徑懸臂管的甩擊過(guò)程來(lái)模擬研究管道的甩擊行為,得到了管道甩擊運(yùn)動(dòng)的影響范圍,并進(jìn)行了甩擊測(cè)試實(shí)驗(yàn)對(duì)結(jié)果進(jìn)行驗(yàn)證。Vongmongkol等[11]和Pieters等[12]研究了管道與防甩擊限制件的碰撞問(wèn)題,初步得到了沖擊力和管道變形等參數(shù)的簡(jiǎn)化計(jì)算方法。Zhao等[13]采用有限元法,結(jié)合RANS方程研究了在橫向流作用下管道排列方式對(duì)甩擊振動(dòng)的影響。王春霖等[14]利用有限元分析軟件LS-DYNA模擬了管道在環(huán)向斷裂條件下的甩擊運(yùn)動(dòng)行為,但他的分析沒(méi)有考慮管道的大變形效應(yīng),其結(jié)果可靠性較差。He等[15]基于流體彈性模型,建立了管道非線性方程,推動(dòng)了管道振動(dòng)鎖定區(qū)域研究計(jì)算方法發(fā)展。Li等[16]采用變分迭代法開(kāi)展了對(duì)兩端彈性支承限制下輸流管道的流固耦合振動(dòng)理論研究,并進(jìn)行了實(shí)驗(yàn)加以驗(yàn)證。劉鋒等[17]采用不同工況下的彎矩-曲率關(guān)系對(duì)管道瞬態(tài)甩擊運(yùn)動(dòng)進(jìn)行了數(shù)值模擬分析。

    通過(guò)以上分析,發(fā)現(xiàn)目前學(xué)者對(duì)流體噴射力影響下管道的甩擊運(yùn)動(dòng)和變形特性進(jìn)行的數(shù)值模擬研究取得了一系列較好的研究成果,但總的來(lái)說(shuō)仍存在不少問(wèn)題。其中,最關(guān)鍵的是大多數(shù)的研究都以懸臂直管作為研究對(duì)象,研究思路都是著眼于固體,流體沖擊力的計(jì)算過(guò)于簡(jiǎn)化,多數(shù)情況下只是簡(jiǎn)單將其作為變化荷載加載在固體上,同時(shí)對(duì)彎管這一類(lèi)更復(fù)雜的情景考慮較少。然而,正如前文所言,管道的甩擊運(yùn)動(dòng)是一種典型的流固耦合問(wèn)題,流體的各項(xiàng)參數(shù)會(huì)受到管道變形的影響,僅從固體出發(fā)無(wú)法真正得到關(guān)于整個(gè)管道的運(yùn)動(dòng),流體對(duì)柔性壓力管道的作用不能簡(jiǎn)單地簡(jiǎn)化為動(dòng)荷載或沖擊荷載。柔性壓力管道由于彈性較大,受流固耦合的影響更加明顯。因此,需要進(jìn)一步從整體出發(fā),將流體和固體統(tǒng)一考慮,求解整個(gè)流固耦合方程,得到柔性壓力管道甩擊運(yùn)動(dòng)的特點(diǎn),并對(duì)影響因素進(jìn)行分析。筆者擬在現(xiàn)有剛性管道失穩(wěn)振動(dòng)分析基礎(chǔ)上針對(duì)柔性管道運(yùn)行過(guò)程中發(fā)生的脫扣斷裂大變形振甩行為進(jìn)行數(shù)值模擬研究,系統(tǒng)分析柔性管道流固耦合作用下斷裂失穩(wěn)運(yùn)動(dòng)規(guī)律及影響因素,為柔性壓力管道安全運(yùn)行及管理提供理論支撐。

    1 柔性壓力管道有限元計(jì)算模型

    流固耦合計(jì)算涉及到流體域和固體結(jié)構(gòu)的數(shù)值計(jì)算,柔性壓力管道彈性大,在流體的作用下會(huì)發(fā)生大變形運(yùn)動(dòng),因此不能忽略管道的形變,必須進(jìn)行雙向流固耦合計(jì)算。ANSYS Workbench平臺(tái)可以很好地解決不同求解器之間的數(shù)據(jù)傳遞,并在耦合求解器中進(jìn)行數(shù)據(jù)交換,對(duì)流固交界面進(jìn)行求解計(jì)算,實(shí)現(xiàn)雙向流固耦合。這里擬采用ANSYS Workbench 平臺(tái)進(jìn)行柔性壓力管道甩擊運(yùn)動(dòng)的模擬計(jì)算。

    1.1 建立幾何模型

    圖1 柔性壓力管道模型Fig. 1 Geometric model of flexible pressure pipeline

    柔性管道由于結(jié)構(gòu)剛度低、可變性強(qiáng)、自由度高等特點(diǎn)多數(shù)應(yīng)用在臨時(shí)工程中,在延展過(guò)程中鋪設(shè)比較隨意,管道起伏不定,存在許多彎曲段,管道的流固耦合失穩(wěn)表現(xiàn)得非常強(qiáng)烈。直角彎是柔性管道布設(shè)中常見(jiàn)的一種彎曲類(lèi)型,選取這種彎曲形狀的柔性壓力管道作為研究對(duì)象,分析甩擊過(guò)程中的形變、應(yīng)力分布和能量分布等問(wèn)題。圖1中管道內(nèi)部為流體域。幾何尺寸參數(shù)如表1所示,包括管道長(zhǎng)度L,管道彎曲半徑R,管道彎曲角度α,管道外徑D,管道內(nèi)徑d以及管壁厚度δ。

    表1 管道幾何參數(shù)

    關(guān)于柔性壓力管道材料物理參數(shù),采用的柔性壓力管道由兩層鋼絲包裹復(fù)合而成,其壁面能承受的最大工作壓力為25 MPa。管道內(nèi)流體為水。具體參數(shù)見(jiàn)表2。

    表2 材料物理參數(shù)

    1.2 計(jì)算網(wǎng)格的劃分

    管道的流固耦合計(jì)算流體域時(shí)需抑制固體管道的幾何體,計(jì)算固體域時(shí)也同樣如此,考慮到柔性壓力管道甩擊運(yùn)動(dòng)幅度大,本次計(jì)算采用適應(yīng)性更好的四面體網(wǎng)格。

    計(jì)算中流體域使用映射面網(wǎng)格方式進(jìn)行劃分,并在流體邊界處進(jìn)行加密,邊界層網(wǎng)格共5層,增長(zhǎng)率為1.3;固體域使用sweep方式進(jìn)行網(wǎng)格劃分。為了提高計(jì)算效率獲得可靠的模擬結(jié)果,計(jì)算過(guò)程中對(duì)網(wǎng)格數(shù)量進(jìn)行無(wú)關(guān)性驗(yàn)證,最終確定流體域和固體域網(wǎng)格尺寸均為3 mm,網(wǎng)格數(shù)分別為112 860和17 930。流體域和固體域網(wǎng)格劃分如圖2所示。

    圖2 網(wǎng)格劃分示意圖Fig. 2 Schematic diagram of meshing

    1.3 模型驗(yàn)證

    管道甩擊是典型的流固耦合問(wèn)題,實(shí)際上也可以看成管道的一種振動(dòng)行為,當(dāng)振動(dòng)幅度達(dá)到一定程度時(shí),管道振動(dòng)發(fā)展成為甩擊。因此,在研究管道甩擊特性的數(shù)值模擬研究過(guò)程中,可以參考輸流管道流固耦合有限元模型,利用Workbench實(shí)現(xiàn)雙向流固耦合仿真分析。陳雄[18]研究90°柔性液壓彎管振動(dòng)特性時(shí)發(fā)現(xiàn),隨著流速的增加,管道最大位移和最大幅值增加,且在恒定流速、脈動(dòng)流速和脈動(dòng)壓力下,管道速度和位移的變化趨勢(shì)大體相同。劉發(fā)明等[19]通過(guò)Workbench中的ANSYS+CFX建立流固耦合模型,研究了脈動(dòng)入口壓力下管道的自由甩動(dòng)情況。如圖3所示,(a)圖為劉發(fā)明等[19]所繪管道自由端位移變化圖,(b)圖為本文仿真時(shí)所繪位移變化圖。將兩者位移變化情況進(jìn)行對(duì)比,可以發(fā)現(xiàn)具有相同的變化趨勢(shì),位移呈現(xiàn)周期性波動(dòng),且位移峰值隨著時(shí)間增加逐漸減小。因此,通過(guò)Workbench進(jìn)行雙向流固耦合仿真計(jì)算能較好地對(duì)管道甩擊特性進(jìn)行研究。

    圖3 管道甩擊位移對(duì)比圖Fig. 3 Comparison of pipe displacement

    2 柔性壓力管道斷裂失穩(wěn)甩擊特性分析

    2.1 變形位移和甩擊速度分析

    管道的最大位移和最大速度可衡量管道甩擊的劇烈程度,在本次計(jì)算中流體流速為v=20 m/s,管道彎曲半徑為1 200 mm,彎曲角度為90°,圖4是管道發(fā)生甩擊的整個(gè)過(guò)程中管道最大總位移和最大總速度變化圖像。從圖4中可以看出在柔性壓力管道發(fā)生甩擊運(yùn)動(dòng)時(shí),管道的最大位移和速度呈周期性變化。當(dāng)流體的流速為20 m/s,管道內(nèi)流體為水的條件下,整個(gè)甩擊過(guò)程中最大變形位移達(dá)到62.0 mm,最大甩擊速度則達(dá)到了1 452.4 mm/s。

    圖4 柔性壓力管道甩擊變形位移和速度變化曲線Fig. 4 Deformation displacement and velocity change curves of flexible pressure pipeline

    隨著時(shí)間的推移,最大位移和最大甩擊速度的峰值都在逐漸衰減。從能量轉(zhuǎn)化的角度來(lái)看,柔性壓力管道發(fā)生變形時(shí),部分能量被管道材料吸收,如果材料發(fā)生了塑性變形,則這部分能量將會(huì)損失,導(dǎo)致管道變形的位移和速度降低。因此,經(jīng)過(guò)一段時(shí)間后,管道的狀態(tài)可能會(huì)到達(dá)穩(wěn)定的變化趨勢(shì)。管道變形位移波動(dòng)的平衡點(diǎn)約為25 mm,表明管道始終存在變形,最終在新的位置重新建立了平衡。但在實(shí)際工程中,這種平衡是很不穩(wěn)定的,一旦流體、管道等任何狀態(tài)發(fā)生變化,管道會(huì)馬上發(fā)生二次甩擊。圖5是管道總體最大變形位移和各個(gè)坐標(biāo)軸方向上位移分量隨時(shí)間的變化情況。由圖5可知,在y軸方向,管道的正向位移始終為0,管道在y軸方向上沒(méi)有發(fā)生扭曲。相反,x軸方向和z軸方向上同時(shí)存在正負(fù)2個(gè)方向上的位移,即管道在這2個(gè)方向上發(fā)生了扭曲。

    圖5 柔性壓力管道總變形位移、甩擊速度和各分量變化曲線Fig. 5 Total deformation displacement, slamming speed and various components change curve

    對(duì)圖4中的柔性管道甩擊變形位移變化曲線做快速傅里葉變換,得到甩擊過(guò)程中振動(dòng)信號(hào)的幅頻曲線如圖6所示。從圖6中可以看出,在甩擊過(guò)程中振動(dòng)信號(hào)主要分布在低頻段,主導(dǎo)頻率為6.65 Hz。對(duì)該系統(tǒng)進(jìn)行模態(tài)分析,得到其前六階固有頻率如表3所示。由此可見(jiàn),在柔性管道甩擊過(guò)程中,被激發(fā)的主要是第二階模態(tài)。

    圖6 柔性壓力管道甩擊振動(dòng)幅頻曲線Fig. 6 Amplitude-frequency curve of flexible pressure pipeline

    表3 系統(tǒng)前六階固有頻率

    2.2 管道受力分析

    等效應(yīng)力是材料力學(xué)理論中將主應(yīng)力的某一綜合值與材料單向拉伸軸向拉壓許用應(yīng)力比較綜合值。等效應(yīng)力越大,管道離發(fā)生失效斷裂就越近。柔性壓力管道的等效應(yīng)力分布情況如圖7所示,管道在y軸負(fù)方向上變形位移最大,在相反的方向上發(fā)生了拉伸變形,等效應(yīng)力分布集中且數(shù)值非常大,由圖8可知在管道內(nèi)水的流速為20 m/s時(shí),等效應(yīng)力最大可達(dá)50.4 MPa,而在管道受到壓縮變形的一側(cè),等效應(yīng)力值很小。

    圖7 柔性壓力管道等效應(yīng)力分布圖Fig. 7 Equivalent stress distribution

    圖8 柔性壓力管道甩擊過(guò)程最大應(yīng)力變化曲線Fig. 8 Maximum stress curve

    管道受到的壓縮變形可以用應(yīng)變能來(lái)描述,應(yīng)變能越大,材料就越容易達(dá)到屈服極限。圖9是柔性壓力管道處于最大變形位移時(shí)的應(yīng)變能分布情況,除了前面發(fā)生拉伸變形的部分,管道靠近變形方向的內(nèi)側(cè)的應(yīng)變能也非常明顯。這表明盡管柔性壓力管道靠近固定約束這一側(cè)變形位移較小,但管道材料內(nèi)部反而受到了強(qiáng)烈的拉伸和壓縮作用。如果管道的變形足夠大,柔性壓力管道的這些位置可能因應(yīng)力過(guò)大發(fā)生二次斷裂,造成更大的人員和財(cái)產(chǎn)損失。

    圖9 柔性壓力管道甩擊運(yùn)動(dòng)應(yīng)變能分布Fig. 9 Strain energy distribution

    3 影響因素分析

    3.1 流體速度

    管道中流體流速由實(shí)際應(yīng)用中所需的工況決定,其大小對(duì)柔性壓力管道的甩擊行為具有明顯影響。為了研究不同流速下管道的甩擊行為,分別對(duì)流體流速為20,30,40 m/s時(shí)柔性壓力管道甩擊運(yùn)動(dòng)進(jìn)行計(jì)算分析。從圖10、11可以發(fā)現(xiàn),柔性壓力管道的最大變形位移、甩擊速度、等效應(yīng)力都隨時(shí)間不斷衰減,衰減速度基本相同。當(dāng)v=20 m/s時(shí)第七個(gè)變形位移峰值約為第一個(gè)峰值的78.9%,而當(dāng)v=40 m/s時(shí),這個(gè)比例為78.1%,說(shuō)明流速增大,管道吸收的能量增加,管道的失效速度加快。

    圖10 不同速度下的最大變形位移和甩擊速度變化曲線Fig. 10 Curves of maximum deformation displacement and slamming speed at different speeds

    圖11 不同流速下的最大應(yīng)力變化曲線Fig. 11 Variation curve of maximum stress at different flow rates

    表4是在不同的流體流速下,柔性壓力管道的甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)。從表4中可以看到,隨著流體速度的增加,甩擊變形位移、甩擊速度、等效應(yīng)力急劇增加,應(yīng)變能的增加速度甚至呈現(xiàn)高次方變化。

    表4 柔性壓力管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)最大值

    對(duì)不同流速下管道最大變形位移變化曲線做快速傅里葉變換,得到不同流速下甩擊過(guò)程中振動(dòng)信號(hào)的幅頻曲線如圖12所示。隨著流速的增加,振動(dòng)信號(hào)的幅值明顯增加,但振動(dòng)信號(hào)的頻率分布及主導(dǎo)頻率的變化則相對(duì)較小。對(duì)不同流速下的該系統(tǒng)進(jìn)行模態(tài)分析,得到不同流速下的前六階固有頻率如圖13所示。從圖13中可以看出隨著流速的增大,系統(tǒng)前六階固有頻率變化較小,在一定情況下可以忽略不考慮。在流速為30 m/s和40 m/s時(shí),柔性壓力管道甩擊振動(dòng)的主導(dǎo)頻率均靠近系統(tǒng)第二階固有頻率??梢?jiàn),管道甩擊過(guò)程中被激發(fā)的仍為第二階模態(tài)。

    圖12 不同流速下柔性壓力管道甩擊振動(dòng)幅頻曲線Fig. 12 Amplitude-frequency curve of flexible pressure pipeline under different flow velocities

    圖13 不同流速系統(tǒng)前六階固有頻率Fig. 13 Natural frequencies of system under different flow velocities

    3.2 流體物性參數(shù)

    柔性壓力管道廣泛應(yīng)用于工程實(shí)踐中,輸送介質(zhì)也是多種多樣,如液壓管道中流體介質(zhì)是液壓油,礦井中的管道很多是水和磨料的混合物。這些流體性質(zhì)差異較大,為了研究不同流體的性質(zhì)對(duì)管道甩擊的影響,對(duì)液壓油、水和含磨料的混合水3種流體進(jìn)行柔性壓力管道的甩擊特性分析。其中混合水由體積分?jǐn)?shù)為8%的石榴石作為磨料混合而成。在計(jì)算中,不考慮管道在徑向發(fā)生的變形,即不改變流體域的體積大小,所以不用考慮流體的壓縮性。3種流體的主要性質(zhì)如表5所示。

    表5 3種流體的基本性質(zhì)

    圖14、15分別是流體為油和含磨料水時(shí),流速v=20 m/s的條件下柔性壓力管道甩擊過(guò)程中的變形位移、甩擊速度、等效應(yīng)力的變化曲線。流體的性質(zhì)變化對(duì)柔性壓力管道甩擊運(yùn)動(dòng)的總體趨勢(shì)沒(méi)有影響,整體來(lái)看都是隨著時(shí)間的增加甩擊的強(qiáng)度在逐漸減弱。

    圖14 不同流體作用下的最大變形位移和甩擊速度變化曲線Fig. 14 Variation curve of maximum deformation displacement and slamming velocity under different fluids

    圖15 不同流體作用下的最大應(yīng)力變化曲線Fig. 15 Variation curve of maximum stress under different fluids

    表6分別是3種不同流體下柔性壓力管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)的最大值。含磨料的水和純水相比,密度更大,能量更高,因此在各項(xiàng)甩擊運(yùn)動(dòng)參數(shù)上都相對(duì)較大,二者的數(shù)值與流體的密度大致成正比。即流體密度越大,則柔性壓力管道的甩擊運(yùn)動(dòng)越劇烈,管道材料承受的應(yīng)力等也越大,危險(xiǎn)性也越高。

    液壓油和水相比,密度更小,但黏度系數(shù)更大,從而導(dǎo)致結(jié)果上的差異。根據(jù)水和含磨料水的比較結(jié)果可知,如果液壓油的黏度系數(shù)和水一樣,那么柔性壓力管道的甩擊應(yīng)該更輕微。由于液壓油具有較大的黏度系數(shù),使其在密度更小的條件下,柔性壓力管道的甩擊劇烈程度遠(yuǎn)超流體為水時(shí)的情況。流體黏度越大,其和柔性壓力管道的耦合效應(yīng)就越強(qiáng)烈,可以將更多的能量傳遞到管道,引起更加劇烈的甩擊運(yùn)動(dòng)。

    表6 不同流體下柔性壓力管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)的最大值

    3.3 管道彎曲長(zhǎng)度

    柔性壓力管道的在鋪設(shè)過(guò)程中,彎曲半徑會(huì)隨實(shí)際情況而改變。為了研究不同的彎曲長(zhǎng)度對(duì)柔性壓力管道甩擊運(yùn)動(dòng)的影響,建立了3種不同長(zhǎng)度的管道模型,長(zhǎng)度分別為600π、900π和1 200π,如圖16所示。

    圖16 3種不同彎曲長(zhǎng)度管道模型Fig. 16 Three kinds of pipe models with different bending lengths

    圖17和圖18是管道長(zhǎng)度為900π和1 200π時(shí)管道甩擊的各項(xiàng)參數(shù)變化曲線。在相同的流速下,管道越長(zhǎng),甩擊運(yùn)動(dòng)的變化周期也就長(zhǎng),管道甩擊的頻率越慢。長(zhǎng)度為1 200π時(shí),其周期約為0.45 s,長(zhǎng)度為600π時(shí),周期只有約0.14 s。

    圖17 不同彎曲長(zhǎng)度下的最大變形位移和甩擊速度變化曲線Fig. 17 Variation curve of maximum deformation displacement and flicking speed under different bending lengths

    圖19是由模態(tài)分析得到的不同彎曲長(zhǎng)度系統(tǒng)的前六階固有頻率,從中可以發(fā)現(xiàn),隨著彎曲長(zhǎng)度的增大,系統(tǒng)前六階固有頻率均顯著減?。浑A數(shù)越高,固有頻率減小幅度越大。因此,彎曲長(zhǎng)度對(duì)系統(tǒng)固有頻率的顯著影響導(dǎo)致了柔性壓力管道振幅的顯著增大與其變化周期的顯著延長(zhǎng)。同時(shí),對(duì)不同彎曲長(zhǎng)度的管道最大變形位移變化曲線做快速傅里葉變換,得到不同彎曲長(zhǎng)度管道甩擊過(guò)程中振動(dòng)信號(hào)的幅頻曲線如圖20所示。隨著彎曲長(zhǎng)度的增加,振動(dòng)信號(hào)的幅值顯著增加,振動(dòng)信號(hào)的頻率分布及主導(dǎo)頻率的變化也發(fā)生顯著變化。此時(shí),起主導(dǎo)作用的仍為二階模態(tài)。

    圖18 不同彎曲長(zhǎng)度下的最大應(yīng)力變化曲線Fig. 18 Maximum stress variation curve under different bending lengths

    圖19 不同彎曲長(zhǎng)度系統(tǒng)前六階固有頻率Fig. 19 Natural frequencies of system of different lengths

    圖20 不同彎曲長(zhǎng)度柔性壓力管道甩擊振動(dòng)幅頻曲線Fig. 20 Amplitude-frequency curve of flexible pressure pipeline of different lengths

    表7是3種不同彎曲長(zhǎng)度下甩擊行為參數(shù)的最大值。從表7中可以看出,隨著柔性壓力管道彎曲段長(zhǎng)度的增加,其甩擊行為變得更加劇烈。管道的彎曲長(zhǎng)度增加1倍,管道的變形位移增加了6倍,甩擊速度則增加了1.5倍。彎曲長(zhǎng)度從600π增加到900π的過(guò)程中,柔性壓力管道上的等效應(yīng)力幾乎保持不變,應(yīng)變能略有增加。當(dāng)管道彎曲長(zhǎng)度從900π增加到1 200π的時(shí)候,管道的等效應(yīng)力和應(yīng)變能出現(xiàn)了劇烈的變化,分別為原來(lái)的3倍和6倍。等效應(yīng)力在彎曲長(zhǎng)度為900π時(shí),管道的變形位移圖像出現(xiàn)了次峰,圖21是t=1 s時(shí)的等效應(yīng)力分布情況。

    表7 不同彎曲長(zhǎng)度下,柔性壓力管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)的最大值

    圖21 t=1 s,管道彎曲段長(zhǎng)900π時(shí),管道等效應(yīng)力分布Fig. 21 t=1 s, the equivalent stress distribution of the pipeline when the bending section of the pipeline is 900π long

    圖21中紅色圓圈出的應(yīng)力集中是導(dǎo)致彎曲長(zhǎng)度為900π時(shí)管道等效應(yīng)力曲線出現(xiàn)次峰的原因。也因此,管道發(fā)生次級(jí)變形,吸收了流體的能量,管道的最大等效應(yīng)力變化也較小。這種現(xiàn)象可能與流體初始速度的方向和約束面存在夾角有關(guān)。

    3.4 管道彎曲半徑

    柔性壓力管道的在鋪設(shè)過(guò)程中,障礙物不同,彎曲半徑也不同。對(duì)彎曲半徑分別為900,1 200,1 500 mm,總長(zhǎng)均為2 000 mm的管道進(jìn)行分析,圖22和圖23是不同彎曲半徑的柔性壓力管道甩擊運(yùn)動(dòng)參數(shù)時(shí)間變化曲線,,隨著曲率半徑的增大,柔性壓力管道的最大變形位移和最大甩擊速度都逐漸減小,等效應(yīng)力也隨曲率半徑的增大而減小,但減小的幅度相比之下要小得多。表8是不同彎曲半徑下管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)的最大值,隨著曲率半徑增大,柔性壓力管道的位移、速度、等效應(yīng)力和應(yīng)變能的最大值都減小。

    圖22 不同彎曲半徑下的最大變形位移和甩擊速度變化曲線Fig. 22 Variation curve of maximum deformation displacement and flicking speed under different bending radii

    圖23 不同彎曲半徑下的最大應(yīng)力變化曲線Fig. 23 Maximum stress variation curve under different bending radii

    表8 不同彎曲半徑下,柔性壓力管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)的最大值

    管道的彎曲半徑增大減小了管道的甩擊變形位移和甩擊速度,減小的幅度較為明顯。同時(shí),增大管道的彎曲半徑還可以減小管道固定端受到的等效應(yīng)力,降低管道的應(yīng)變能,增強(qiáng)管道在甩擊運(yùn)動(dòng)過(guò)程中的安全性。

    圖24是由模態(tài)分析得到的不同彎曲長(zhǎng)度系統(tǒng)的前六階固有頻率,從中可以發(fā)現(xiàn),隨著彎曲半徑的增大,系統(tǒng)前兩階固有頻率減小,三到六階固有頻率顯著增大;階數(shù)越高,固有頻率變化幅度越大。同時(shí),對(duì)不同彎曲半徑的管道最大變形位移變化曲線做快速傅里葉變換,得到不同彎曲半徑管道甩擊過(guò)程中振動(dòng)信號(hào)的幅頻曲線如圖25所示。隨著彎曲半徑的增加,振動(dòng)信號(hào)的幅值與主導(dǎo)頻率發(fā)生變化。此時(shí),起主導(dǎo)作用的仍為二階模態(tài)。

    圖24 不同彎曲半徑系統(tǒng)前六階固有頻率Fig. 24 Natural frequencies of system of different curvature radiuses

    圖25 不同彎曲半徑柔性壓力管道甩擊振動(dòng)幅頻曲線Fig. 25 Amplitude-frequency curve of flexible pressure pipeline of different curvature radiuses

    3.5 約束位置

    柔性壓力管道在工作時(shí),一般每隔一段距離就設(shè)置一個(gè)安全扣進(jìn)行約束,安全扣和彎曲管道的位置關(guān)系對(duì)管道另一端脫扣后的甩擊行為有很大的影響。管道的幾何模型如圖26所示,對(duì)d=0,300,600 mm的情況進(jìn)行管道甩擊分析。圖27和28是不同約束位置下管道的甩擊運(yùn)動(dòng)參數(shù)變化曲線。

    圖26 不同約束位置幾何模型Fig. 26 Geometric models with different constrained positions

    圖27 不同約束位置下變形位移和甩擊速度變化曲線Fig. 27 Curves of deformation displacement and slamming speed under different constrained positions

    圖28 不同約束位置下的最大應(yīng)力變化曲線Fig. 28 Maximum stress change curve under different constrained positions

    隨著約束點(diǎn)離彎曲段的距離不斷增大,管道甩擊的整體最大位移不斷減小,最大甩擊速度也不斷減小,但管道靠近固定端的等效應(yīng)力變化不大。從表9可以看出,隨著約束點(diǎn)離彎曲段的距離不斷增大,柔性壓力管道的位移、速度和應(yīng)變能的最大值都減小。從圖29、30可以發(fā)現(xiàn),隨著約束點(diǎn)距離管道彎曲段越來(lái)越遠(yuǎn),彎曲段的變形由于沒(méi)有約束的存在會(huì)繼續(xù)向直管段傳遞,能量也傳遞到直管部分,使得更多的管道變形,每個(gè)單位的應(yīng)力能反而下降。約束點(diǎn)距離管道彎曲段越遠(yuǎn),這種效應(yīng)越顯著,對(duì)防止管道發(fā)生二次斷裂更有利。但約束點(diǎn)越遠(yuǎn),則管道甩擊運(yùn)動(dòng)的范圍會(huì)越大,對(duì)人員和設(shè)備的危險(xiǎn)性反而提升。因此,需要綜合考慮對(duì)管道的約束。

    表9 不同約束點(diǎn)下柔性壓力管道甩擊運(yùn)動(dòng)各項(xiàng)參數(shù)的最大值

    圖29 d=300 mm,t=1 s時(shí),管道變形位移和等效應(yīng)力分布圖Fig. 29 d= 300 mm, t = 1 s, the distribution of deformation and equivalent stress of pipeline

    4 結(jié) 論

    1)管道發(fā)生脫扣甩擊時(shí),運(yùn)動(dòng)具有很強(qiáng)的非線性特征,最大變形位移和甩擊速度的空間分布沒(méi)有明顯規(guī)律。管道靠近自由端的變形位移和甩擊速度最大,對(duì)人員和設(shè)備危險(xiǎn)性高,而固定端的等效應(yīng)力最大,容易產(chǎn)生疲勞破壞造成二次危害。

    2)流體流速增加,柔性壓力管道甩擊變形位移、甩擊速度、等效應(yīng)力也迅速增大,應(yīng)變能受流速的影響更大,與流速呈高次方增加的關(guān)系。

    3)流體對(duì)柔性壓力管道甩擊行為的影響主要體現(xiàn)在流體的密度和黏度系數(shù)兩個(gè)方面。隨著流體密度增加,柔性壓力管道甩擊運(yùn)動(dòng)的主要參數(shù)都線性增加。而流體黏度對(duì)柔性壓力管道的甩擊影響較為復(fù)雜,它會(huì)影響管道和流體的耦合效應(yīng),進(jìn)而影響管道的甩擊行為。隨著流體黏度的增加,管道的甩擊行為更加劇烈。

    4)柔性壓力管道的長(zhǎng)度越長(zhǎng),甩擊就越劇烈,管道的疲勞損壞風(fēng)險(xiǎn)大大增加,其運(yùn)動(dòng)周期也延長(zhǎng)。

    5)管道的彎曲半徑增大,管道的甩擊變形位移和甩擊速度明顯減小。增大彎曲半徑,還可減小管道固定端的等效應(yīng)力,降低管道的應(yīng)變能,改善柔性壓力管道的受力情況,延長(zhǎng)管道使用壽命;此外,利于延長(zhǎng)甩擊運(yùn)動(dòng)的周期,極大地提高柔性壓力管道的總體安全性。

    6)改變管道安全扣的約束位置的本質(zhì)是在調(diào)節(jié)管道運(yùn)行的長(zhǎng)度和彎曲半徑兩個(gè)參數(shù),它們的變化規(guī)律共同影響了約束點(diǎn)對(duì)甩擊運(yùn)動(dòng)的影響。

    7)通過(guò)頻譜分析和模態(tài)分析發(fā)現(xiàn),管道振動(dòng)都是由二階模態(tài)主導(dǎo)。

    猜你喜歡
    柔性流速流體
    “流體壓強(qiáng)與流速的關(guān)系”知識(shí)鞏固
    一種柔性拋光打磨頭設(shè)計(jì)
    『流體壓強(qiáng)與流速的關(guān)系』知識(shí)鞏固
    流體壓強(qiáng)知多少
    灌注式半柔性路面研究進(jìn)展(1)——半柔性混合料組成設(shè)計(jì)
    石油瀝青(2021年5期)2021-12-02 03:21:18
    高校學(xué)生管理工作中柔性管理模式應(yīng)用探索
    山雨欲來(lái)風(fēng)滿(mǎn)樓之流體壓強(qiáng)與流速
    愛(ài)虛張聲勢(shì)的水
    等效流體體積模量直接反演的流體識(shí)別方法
    柔性的思維
    棋藝(2001年7期)2001-06-14 03:21:00
    午夜91福利影院| 久久狼人影院| 91老司机精品| 狠狠精品人妻久久久久久综合| 高清在线国产一区| 制服人妻中文乱码| 人人澡人人妻人| 精品久久久精品久久久| 男人爽女人下面视频在线观看| 操出白浆在线播放| 十八禁高潮呻吟视频| 久久久久久久大尺度免费视频| 午夜精品国产一区二区电影| 精品高清国产在线一区| 日本wwww免费看| 中文字幕另类日韩欧美亚洲嫩草| 国产av一区二区精品久久| www.自偷自拍.com| av在线app专区| 亚洲国产成人一精品久久久| 午夜视频精品福利| 免费在线观看影片大全网站| 日本五十路高清| 亚洲精品一区蜜桃| 国产深夜福利视频在线观看| 午夜激情av网站| 曰老女人黄片| 少妇 在线观看| 欧美性长视频在线观看| 亚洲欧美精品综合一区二区三区| 国产精品99久久99久久久不卡| 成人国产一区最新在线观看| 久久久久精品人妻al黑| 免费观看a级毛片全部| 久久久久视频综合| 亚洲精品久久午夜乱码| 亚洲精品第二区| 久久天堂一区二区三区四区| 男女床上黄色一级片免费看| 中文字幕人妻熟女乱码| 高清av免费在线| 免费高清在线观看视频在线观看| 亚洲五月婷婷丁香| 1024视频免费在线观看| 亚洲一码二码三码区别大吗| 国产av一区二区精品久久| 久久久久久久国产电影| 男人爽女人下面视频在线观看| av免费在线观看网站| 在线 av 中文字幕| 9热在线视频观看99| 亚洲欧美日韩另类电影网站| 精品国产国语对白av| 又黄又粗又硬又大视频| 真人做人爱边吃奶动态| 国产成人欧美在线观看 | 在线观看免费午夜福利视频| 日韩,欧美,国产一区二区三区| 人人妻人人澡人人爽人人夜夜| 男女下面插进去视频免费观看| 美女主播在线视频| 免费不卡黄色视频| 久久综合国产亚洲精品| 国产成人精品久久二区二区免费| 曰老女人黄片| videos熟女内射| 91麻豆精品激情在线观看国产 | 纵有疾风起免费观看全集完整版| 在线观看免费午夜福利视频| 色综合欧美亚洲国产小说| 亚洲 欧美一区二区三区| 午夜激情av网站| 黑人巨大精品欧美一区二区蜜桃| 人人妻,人人澡人人爽秒播| 国产深夜福利视频在线观看| av天堂在线播放| 丰满迷人的少妇在线观看| 久久久国产一区二区| 亚洲av片天天在线观看| 丝袜脚勾引网站| 欧美日韩av久久| 建设人人有责人人尽责人人享有的| 精品亚洲乱码少妇综合久久| 欧美日韩亚洲高清精品| av又黄又爽大尺度在线免费看| 成人国语在线视频| 精品一区在线观看国产| 一级毛片精品| 亚洲自偷自拍图片 自拍| 日日摸夜夜添夜夜添小说| 色视频在线一区二区三区| 黑人欧美特级aaaaaa片| 日日爽夜夜爽网站| 亚洲情色 制服丝袜| 亚洲精品av麻豆狂野| 午夜日韩欧美国产| 三级毛片av免费| 欧美日韩亚洲国产一区二区在线观看 | 免费女性裸体啪啪无遮挡网站| 亚洲一区二区三区欧美精品| 久久久久久久久久久久大奶| xxxhd国产人妻xxx| 成人国语在线视频| 国产成人欧美在线观看 | 国产精品香港三级国产av潘金莲| 黄色视频不卡| 国产成人精品无人区| 国产欧美日韩综合在线一区二区| 看免费av毛片| 国产无遮挡羞羞视频在线观看| 十八禁网站网址无遮挡| 欧美日韩福利视频一区二区| 亚洲国产欧美网| 三上悠亚av全集在线观看| 色94色欧美一区二区| 青草久久国产| 国精品久久久久久国模美| 嫩草影视91久久| 午夜福利视频在线观看免费| 丝袜脚勾引网站| 十八禁高潮呻吟视频| 性色av一级| 欧美日韩av久久| 美女主播在线视频| 亚洲精品久久久久久婷婷小说| 中文字幕人妻熟女乱码| 麻豆国产av国片精品| 午夜福利在线观看吧| 一进一出抽搐动态| 精品久久蜜臀av无| 少妇的丰满在线观看| 宅男免费午夜| 91麻豆精品激情在线观看国产 | 欧美精品人与动牲交sv欧美| 黄色片一级片一级黄色片| 国产日韩一区二区三区精品不卡| 一本久久精品| 一个人免费看片子| 国产男女超爽视频在线观看| 成人国语在线视频| 国产精品一区二区精品视频观看| 热99久久久久精品小说推荐| 91国产中文字幕| 成人国语在线视频| 中国国产av一级| 国产精品一区二区在线观看99| 9191精品国产免费久久| 日日摸夜夜添夜夜添小说| 99久久99久久久精品蜜桃| av免费在线观看网站| av天堂在线播放| 一区二区日韩欧美中文字幕| 十八禁网站免费在线| 亚洲精品粉嫩美女一区| 亚洲中文av在线| 精品久久久久久久毛片微露脸 | 97精品久久久久久久久久精品| 女性被躁到高潮视频| 999精品在线视频| 97在线人人人人妻| 他把我摸到了高潮在线观看 | 18禁观看日本| 久久精品久久久久久噜噜老黄| 美女主播在线视频| 精品一区二区三区四区五区乱码| 国产麻豆69| 欧美人与性动交α欧美精品济南到| 欧美日韩亚洲国产一区二区在线观看 | 成年美女黄网站色视频大全免费| 性高湖久久久久久久久免费观看| 久久99一区二区三区| 高清在线国产一区| 777久久人妻少妇嫩草av网站| 国产亚洲av高清不卡| 啦啦啦啦在线视频资源| 免费人妻精品一区二区三区视频| av超薄肉色丝袜交足视频| 国产色视频综合| 日本av免费视频播放| 自拍欧美九色日韩亚洲蝌蚪91| 黑人操中国人逼视频| 老熟妇仑乱视频hdxx| 黄片小视频在线播放| 一级片免费观看大全| 丰满人妻熟妇乱又伦精品不卡| 婷婷色av中文字幕| 在线看a的网站| 两个人看的免费小视频| 人人妻,人人澡人人爽秒播| 交换朋友夫妻互换小说| 国产一区二区 视频在线| 美女视频免费永久观看网站| 久久99一区二区三区| 伦理电影免费视频| 亚洲 国产 在线| 天堂中文最新版在线下载| 午夜精品久久久久久毛片777| 99久久综合免费| 丁香六月天网| 久久国产精品大桥未久av| 免费久久久久久久精品成人欧美视频| 水蜜桃什么品种好| 巨乳人妻的诱惑在线观看| 久久久久国内视频| 精品国产国语对白av| 男人舔女人的私密视频| 国产成人精品无人区| 黄色a级毛片大全视频| 国产视频一区二区在线看| 欧美黄色片欧美黄色片| 极品人妻少妇av视频| 在线观看免费高清a一片| 啦啦啦免费观看视频1| 日本av免费视频播放| 国产亚洲精品一区二区www | www日本在线高清视频| 91麻豆精品激情在线观看国产 | 久热这里只有精品99| 国产一卡二卡三卡精品| 肉色欧美久久久久久久蜜桃| 高清av免费在线| 在线 av 中文字幕| 国产欧美亚洲国产| 丰满迷人的少妇在线观看| av天堂在线播放| cao死你这个sao货| 91老司机精品| 乱人伦中国视频| 一级,二级,三级黄色视频| 日韩中文字幕视频在线看片| 国产一级毛片在线| 免费女性裸体啪啪无遮挡网站| 色视频在线一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看人在逋| 人人澡人人妻人| 国产一区有黄有色的免费视频| 欧美精品一区二区免费开放| 一边摸一边抽搐一进一出视频| 成年人免费黄色播放视频| 中文精品一卡2卡3卡4更新| 欧美亚洲日本最大视频资源| 久久天躁狠狠躁夜夜2o2o| av天堂久久9| 精品亚洲乱码少妇综合久久| 国产色视频综合| 中文字幕精品免费在线观看视频| 岛国毛片在线播放| 成人三级做爰电影| 亚洲欧美日韩另类电影网站| 国产福利在线免费观看视频| 亚洲色图综合在线观看| 在线永久观看黄色视频| 久久久精品国产亚洲av高清涩受| 欧美 日韩 精品 国产| 性色av乱码一区二区三区2| 91老司机精品| 久久久欧美国产精品| 精品久久久久久电影网| 亚洲熟女毛片儿| 精品国产乱码久久久久久男人| 淫妇啪啪啪对白视频 | 国产麻豆69| 下体分泌物呈黄色| 国产精品一区二区在线观看99| 各种免费的搞黄视频| 蜜桃国产av成人99| 久久久久国内视频| 久久99热这里只频精品6学生| 少妇粗大呻吟视频| 久久女婷五月综合色啪小说| 最近中文字幕2019免费版| 一个人免费在线观看的高清视频 | 欧美黄色淫秽网站| 亚洲熟女毛片儿| 黄色视频,在线免费观看| 一本色道久久久久久精品综合| 亚洲欧美日韩另类电影网站| 国产精品九九99| 久久这里只有精品19| 精品高清国产在线一区| 中文字幕制服av| 母亲3免费完整高清在线观看| 亚洲精品国产一区二区精华液| 国产精品免费视频内射| 日韩三级视频一区二区三区| 亚洲成国产人片在线观看| 中文字幕精品免费在线观看视频| 午夜福利乱码中文字幕| 日韩有码中文字幕| 99国产精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲人成电影免费在线| 色婷婷久久久亚洲欧美| 99精品欧美一区二区三区四区| 一本色道久久久久久精品综合| 精品福利观看| 国产精品九九99| 最新在线观看一区二区三区| 国产视频一区二区在线看| 亚洲精品粉嫩美女一区| 首页视频小说图片口味搜索| 免费在线观看视频国产中文字幕亚洲 | 国精品久久久久久国模美| 老司机亚洲免费影院| 两个人看的免费小视频| 亚洲欧美日韩另类电影网站| 男女床上黄色一级片免费看| 嫁个100分男人电影在线观看| 欧美另类一区| 一边摸一边做爽爽视频免费| 捣出白浆h1v1| 日本av免费视频播放| 亚洲av片天天在线观看| 精品一区二区三卡| 超色免费av| 久久精品人人爽人人爽视色| 性少妇av在线| a级毛片黄视频| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久综合国产亚洲精品| 亚洲国产av影院在线观看| 亚洲色图综合在线观看| 欧美日韩亚洲高清精品| 一本一本久久a久久精品综合妖精| 两性夫妻黄色片| 制服人妻中文乱码| 人人妻,人人澡人人爽秒播| 51午夜福利影视在线观看| 正在播放国产对白刺激| 午夜视频精品福利| 亚洲av美国av| 精品熟女少妇八av免费久了| 我的亚洲天堂| 日韩 亚洲 欧美在线| 亚洲第一av免费看| 亚洲av日韩在线播放| 99久久精品国产亚洲精品| 亚洲精品日韩在线中文字幕| 性少妇av在线| 欧美激情久久久久久爽电影 | 久久国产亚洲av麻豆专区| 国产又色又爽无遮挡免| 久久久久精品国产欧美久久久 | av网站免费在线观看视频| √禁漫天堂资源中文www| 国产在线观看jvid| 一区在线观看完整版| 久久99热这里只频精品6学生| 飞空精品影院首页| 日本精品一区二区三区蜜桃| 精品国产乱码久久久久久男人| 考比视频在线观看| 大型av网站在线播放| av网站在线播放免费| 岛国在线观看网站| 日韩欧美一区二区三区在线观看 | 天堂中文最新版在线下载| 男女下面插进去视频免费观看| 1024视频免费在线观看| 欧美大码av| 丰满饥渴人妻一区二区三| avwww免费| 51午夜福利影视在线观看| 亚洲欧洲精品一区二区精品久久久| av免费在线观看网站| 久久久久精品人妻al黑| 国产97色在线日韩免费| 另类精品久久| 亚洲一区二区三区欧美精品| 99热全是精品| 汤姆久久久久久久影院中文字幕| 欧美性长视频在线观看| 国产福利在线免费观看视频| 黄色视频在线播放观看不卡| 久久精品成人免费网站| 色婷婷久久久亚洲欧美| 最近最新免费中文字幕在线| 在线 av 中文字幕| 国产伦人伦偷精品视频| 亚洲va日本ⅴa欧美va伊人久久 | 两人在一起打扑克的视频| 国产精品久久久久久精品电影小说| 男女之事视频高清在线观看| 窝窝影院91人妻| kizo精华| 黄片播放在线免费| 国产欧美日韩一区二区三 | 亚洲精品中文字幕在线视频| 久久精品亚洲av国产电影网| 久久久久久免费高清国产稀缺| 一区二区三区精品91| 日韩一区二区三区影片| 99精品欧美一区二区三区四区| 国产精品免费大片| 久久久久久久大尺度免费视频| 老鸭窝网址在线观看| 国产精品久久久久久精品电影小说| 天天操日日干夜夜撸| 黑人猛操日本美女一级片| 午夜福利视频在线观看免费| 制服诱惑二区| 91九色精品人成在线观看| a级片在线免费高清观看视频| 麻豆av在线久日| 国产精品免费大片| 久久久国产成人免费| 亚洲成人手机| 午夜福利视频在线观看免费| 欧美黄色片欧美黄色片| 日韩有码中文字幕| 欧美激情久久久久久爽电影 | 欧美日韩av久久| 女性生殖器流出的白浆| 中文精品一卡2卡3卡4更新| 久久久久久人人人人人| 51午夜福利影视在线观看| 国产免费一区二区三区四区乱码| 欧美日韩亚洲综合一区二区三区_| 亚洲成人手机| 香蕉丝袜av| a级毛片黄视频| 天天操日日干夜夜撸| 精品少妇一区二区三区视频日本电影| 日韩人妻精品一区2区三区| av免费在线观看网站| 永久免费av网站大全| 久久久久国产精品人妻一区二区| 国产黄频视频在线观看| 午夜两性在线视频| 欧美日韩一级在线毛片| 国产高清视频在线播放一区 | 99国产极品粉嫩在线观看| 性少妇av在线| h视频一区二区三区| 不卡av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 女人精品久久久久毛片| 天天躁狠狠躁夜夜躁狠狠躁| 正在播放国产对白刺激| 欧美激情 高清一区二区三区| 男女免费视频国产| 精品国产一区二区久久| av在线播放精品| 一区二区日韩欧美中文字幕| 亚洲五月婷婷丁香| 国产精品熟女久久久久浪| 国产精品九九99| 成人国产av品久久久| 50天的宝宝边吃奶边哭怎么回事| 国产又色又爽无遮挡免| 久久久久精品国产欧美久久久 | 亚洲av男天堂| 老熟妇乱子伦视频在线观看 | 美女扒开内裤让男人捅视频| 久久国产精品男人的天堂亚洲| 日本撒尿小便嘘嘘汇集6| 欧美日韩亚洲高清精品| 亚洲精品久久久久久婷婷小说| 日韩 亚洲 欧美在线| 亚洲精品中文字幕一二三四区 | 亚洲一区中文字幕在线| 18禁国产床啪视频网站| 性少妇av在线| 国产深夜福利视频在线观看| 亚洲人成电影免费在线| 老汉色∧v一级毛片| 国产精品偷伦视频观看了| 亚洲欧美精品综合一区二区三区| 亚洲欧美精品自产自拍| 一本综合久久免费| 日本黄色日本黄色录像| 青青草视频在线视频观看| 视频区欧美日本亚洲| 中文字幕精品免费在线观看视频| 亚洲精品久久久久久婷婷小说| 三上悠亚av全集在线观看| 18禁裸乳无遮挡动漫免费视频| 久久精品aⅴ一区二区三区四区| 人人澡人人妻人| 久久国产精品大桥未久av| 手机成人av网站| tocl精华| 亚洲欧美激情在线| 精品久久久久久电影网| 久久久精品免费免费高清| 久久久久久久精品精品| 狠狠婷婷综合久久久久久88av| 侵犯人妻中文字幕一二三四区| 国内毛片毛片毛片毛片毛片| 国产视频一区二区在线看| 国产国语露脸激情在线看| 亚洲第一av免费看| 精品一区二区三区av网在线观看 | 男女边摸边吃奶| 在线亚洲精品国产二区图片欧美| 成年女人毛片免费观看观看9 | 性色av一级| 日本vs欧美在线观看视频| 国产野战对白在线观看| 91大片在线观看| 久久久水蜜桃国产精品网| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美另类亚洲清纯唯美| 亚洲成av片中文字幕在线观看| 欧美日韩国产mv在线观看视频| 久久免费观看电影| 久久人人97超碰香蕉20202| av片东京热男人的天堂| 他把我摸到了高潮在线观看 | av天堂在线播放| 欧美另类一区| 国产成人精品在线电影| 久热爱精品视频在线9| 99久久国产精品久久久| 久久亚洲精品不卡| 一级毛片精品| 99国产精品一区二区蜜桃av | 一二三四在线观看免费中文在| 亚洲人成电影免费在线| 一区福利在线观看| 真人做人爱边吃奶动态| 日韩人妻精品一区2区三区| 超碰97精品在线观看| 亚洲av男天堂| 国产伦理片在线播放av一区| 国产成+人综合+亚洲专区| 黄片大片在线免费观看| 狠狠精品人妻久久久久久综合| 女人精品久久久久毛片| 人妻一区二区av| 免费在线观看日本一区| 两个人免费观看高清视频| 777米奇影视久久| 成年动漫av网址| 精品久久蜜臀av无| 我的亚洲天堂| kizo精华| av国产精品久久久久影院| 国产成人精品久久二区二区91| 少妇粗大呻吟视频| 国产成人啪精品午夜网站| 性高湖久久久久久久久免费观看| 曰老女人黄片| a在线观看视频网站| 丁香六月天网| 黄片播放在线免费| 狠狠狠狠99中文字幕| 在线观看舔阴道视频| 电影成人av| 三上悠亚av全集在线观看| 国产免费福利视频在线观看| 蜜桃国产av成人99| 2018国产大陆天天弄谢| 韩国高清视频一区二区三区| 久久久精品国产亚洲av高清涩受| 久久天堂一区二区三区四区| 99九九在线精品视频| 久久精品国产亚洲av高清一级| 97精品久久久久久久久久精品| 日韩制服丝袜自拍偷拍| 制服诱惑二区| 成人av一区二区三区在线看 | 亚洲精品av麻豆狂野| a级片在线免费高清观看视频| 欧美久久黑人一区二区| 巨乳人妻的诱惑在线观看| 成人国产av品久久久| 国产一区二区三区av在线| 999精品在线视频| 亚洲一码二码三码区别大吗| 搡老熟女国产l中国老女人| 亚洲少妇的诱惑av| 久久久久精品人妻al黑| 人人澡人人妻人| 精品一品国产午夜福利视频| 成年人黄色毛片网站| a级毛片黄视频| 无限看片的www在线观看| 热re99久久精品国产66热6| 亚洲国产毛片av蜜桃av| 性少妇av在线| 亚洲少妇的诱惑av| 操出白浆在线播放| 国产在线一区二区三区精| 国产精品成人在线| 三上悠亚av全集在线观看| 交换朋友夫妻互换小说| www日本在线高清视频| 国产男女超爽视频在线观看| 中文字幕av电影在线播放| 国产亚洲精品第一综合不卡| 十分钟在线观看高清视频www| 国产在线观看jvid| 香蕉国产在线看| 999久久久精品免费观看国产| 电影成人av| 亚洲欧美激情在线| 51午夜福利影视在线观看| 成年人午夜在线观看视频| 黑丝袜美女国产一区| 2018国产大陆天天弄谢| 国产一卡二卡三卡精品| 国产亚洲午夜精品一区二区久久| 人妻一区二区av| 高清av免费在线| 欧美精品高潮呻吟av久久| 叶爱在线成人免费视频播放| 国产精品久久久久久精品电影小说| 人成视频在线观看免费观看| 大片免费播放器 马上看| 免费观看av网站的网址| av电影中文网址| 亚洲欧美一区二区三区黑人|