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

    半潛式平臺(tái)極限強(qiáng)度可靠性研究

    2011-02-27 09:06:44祁恩榮龐建華吳東偉
    船舶力學(xué) 2011年4期
    關(guān)鍵詞:方法

    祁恩榮,龐建華,吳東偉

    (1中國(guó)船舶科學(xué)研究中心,江蘇 無錫 214082;2江蘇科技大學(xué),江蘇 鎮(zhèn)江 212003)

    1 前 言

    21世紀(jì)是海洋的世紀(jì),海洋成為國(guó)家綜合國(guó)力競(jìng)爭(zhēng)的重要制高點(diǎn),經(jīng)略海洋、興海強(qiáng)國(guó)已經(jīng)成為我國(guó)國(guó)家發(fā)展戰(zhàn)略的重要內(nèi)容。一面是日益枯竭的陸上油氣資源,一面是海洋這座尚未充分挖掘的能源寶庫,在能源需求不斷增加的未來,向海洋進(jìn)軍是當(dāng)前世界油氣開發(fā)的必然選擇。半潛式平臺(tái)長(zhǎng)期在深海作業(yè),不可避免地會(huì)遭遇惡劣海況,由風(fēng)、浪和流引起的極值環(huán)境載荷可能超過平臺(tái)結(jié)構(gòu)最大承載能力,從而造成人員、平臺(tái)和海洋環(huán)境的嚴(yán)重?fù)p失。極限強(qiáng)度是半潛式平臺(tái)海洋環(huán)境適應(yīng)能力的顯示指標(biāo),為了獲得安全和經(jīng)濟(jì)的結(jié)構(gòu)設(shè)計(jì),需要精確評(píng)估半潛式平臺(tái)極限承載能力。

    長(zhǎng)期以來船體極限強(qiáng)度研究受到各國(guó)船舶界的廣泛重視,國(guó)內(nèi)外都進(jìn)行了大量的研究工作。目前船體極限強(qiáng)度評(píng)估方法可以分為三類[1]:

    ·實(shí)船事故調(diào)查和模型試驗(yàn);

    · 直接方法,如線彈性方法、經(jīng)驗(yàn)公式和解析方法(AM),以及

    · 逐步破壞法,如簡(jiǎn)化方法(SM)、理想結(jié)構(gòu)單元法(ISUM)和非線性有限元法(FEM)。

    其中,解析方法和簡(jiǎn)化方法是面向設(shè)計(jì)的方法。國(guó)際船舶結(jié)構(gòu)大會(huì)(ISSC)對(duì)船體及其構(gòu)件極限強(qiáng)度分析方法進(jìn)行了多次比較研究[2-3]。國(guó)際船級(jí)社協(xié)會(huì)(IACS)于2006年相繼推出正式的雙殼油船和散貨船共同規(guī)范,在這些規(guī)范中明確提出了船體極限強(qiáng)度的校核要求,采用一步法和簡(jiǎn)化方法計(jì)算船體極限強(qiáng)度,也可選擇有限元方法計(jì)算船體極限強(qiáng)度[4-5]。祁恩榮等[6-7]采用解析方法、簡(jiǎn)化方法、理想結(jié)構(gòu)單元法和非線性有限元法對(duì)雙殼油船和液化天然氣船極限強(qiáng)度進(jìn)行了比較研究。然而各類方法在半潛式平臺(tái)極限強(qiáng)度計(jì)算方面應(yīng)用得較少[8],也缺乏比較分析和模型試驗(yàn)的驗(yàn)證。

    由于承載能力和環(huán)境載荷具有明顯的隨機(jī)性,應(yīng)基于極限強(qiáng)度,采用可靠性方法評(píng)價(jià)半潛式平臺(tái)結(jié)構(gòu)安全性。本文首先闡述了簡(jiǎn)化方法和解析方法的基本步驟,并對(duì)某目標(biāo)平臺(tái)總縱極限強(qiáng)度進(jìn)行了比較研究,然后基于三維勢(shì)流理論和Morison方程對(duì)生存海況下的目標(biāo)平臺(tái)垂向波浪彎矩進(jìn)行了預(yù)報(bào),最后基于極限強(qiáng)度建立了半潛式平臺(tái)可靠性計(jì)算模型,并對(duì)目標(biāo)平臺(tái)結(jié)構(gòu)可靠性和目標(biāo)可靠度進(jìn)行了分析。

    2 平臺(tái)極限強(qiáng)度比較研究

    解析方法和簡(jiǎn)化方法是面向設(shè)計(jì)的方法,并得到了船體極限強(qiáng)度比較分析和模型試驗(yàn)的驗(yàn)證。解析方法將平臺(tái)橫剖面離散化為加筋板格,計(jì)算加筋板格屈曲極限強(qiáng)度;假設(shè)橫剖面極限狀態(tài)應(yīng)力分布,確定中和軸位置,計(jì)算平臺(tái)極限強(qiáng)度。而簡(jiǎn)化方法將平臺(tái)橫剖面離散化為加筋板單元和硬角單元,推導(dǎo)加筋板單元平均應(yīng)力—應(yīng)變關(guān)系;基于平斷面假設(shè),忽略各單元間的相互影響,計(jì)算平臺(tái)中拱和中垂彎矩與曲率的變化歷程。

    解析方法的基本步驟是[9]:

    ·將剖面離散化為加筋板格,利用彈塑性方法計(jì)算其屈曲極限強(qiáng)度;

    ·極限狀態(tài)時(shí)剖面拉伸邊緣屈服,壓縮邊緣屈曲,而在剖面中和軸附近保持線彈性狀態(tài);

    ·剖面彈性區(qū)域由完全屈服和屈曲應(yīng)力分布模型中拉伸力心和壓縮力心在垂直于中和軸方向的位置確定;

    ·極限狀態(tài)剖面中和軸的位置和方向由平衡條件確定;

    ·極限彎矩可表示為彈塑性應(yīng)力分布模型中拉伸力與拉伸力心和壓縮力心之間的距離的乘積。

    解析方法的一個(gè)關(guān)鍵問題是加筋板格極限強(qiáng)度的計(jì)算。彈塑性方法[9]將加筋板格的失效模式分為四類:(1)模式A—縱、橫兩向加強(qiáng)筋與板同時(shí)發(fā)生屈曲失效;(2)模式B—兩橫向加強(qiáng)筋間板與縱向加強(qiáng)筋同時(shí)發(fā)生屈曲失效;(3)模式C—筋間板屈曲引起的縱向加強(qiáng)筋屈服或屈曲失效;(4)模式D—筋間板屈曲引起的縱向加強(qiáng)筋扭曲失效。

    簡(jiǎn)化方法的基本步驟是[6]:

    ·將橫剖面離散化為加筋板單元和硬角單元,單元之間的相互影響忽略不計(jì);

    ·利用平均應(yīng)力—應(yīng)變關(guān)系計(jì)算各單元的軸向剛度;

    ·利用單元軸向剛度計(jì)算橫剖面的抗彎剛度;

    ·基于平斷面假設(shè),逐步增加垂向和水平曲率,從而產(chǎn)生對(duì)瞬時(shí)中和軸的彎曲;

    ·計(jì)算相應(yīng)的彎矩增量,以及各單元的應(yīng)力應(yīng)變?cè)隽浚?/p>

    ·累加各增量,從而得到結(jié)構(gòu)逐步破壞的彎矩—曲率變化歷程。

    簡(jiǎn)化方法的關(guān)鍵是確定各單元的平均應(yīng)力—應(yīng)變關(guān)系。為了可靠和快速地計(jì)算加筋板單元平均應(yīng)力—應(yīng)變關(guān)系,基于梁柱理論的簡(jiǎn)化方法是值得關(guān)注的方法,該方法利用彈塑性梁柱理論推導(dǎo)加筋板單元平均應(yīng)力—應(yīng)變關(guān)系[6]。對(duì)于壓縮單元,極限強(qiáng)度是由加強(qiáng)筋端部壓縮破壞引發(fā)的失效模式和由帶板壓縮破壞引發(fā)的失效模式的軸向壓縮應(yīng)力的較小者,平均應(yīng)力—應(yīng)變關(guān)系分為三個(gè)區(qū)域:穩(wěn)定區(qū)、非卸載區(qū)和卸載區(qū)。

    本文選用的目標(biāo)平臺(tái)是NDB半潛式平臺(tái),其結(jié)構(gòu)原型為1986年俄羅斯建造的“shelf-6”,屬于第 2~3代半潛平臺(tái),經(jīng)美國(guó)NOBEL公司改建后,目前具備3 000m深作業(yè)能力。目標(biāo)平臺(tái)總長(zhǎng)111.6m,總寬66.4m,橫剖面如圖1所示。目標(biāo)平臺(tái)由三部分構(gòu)成:(1)上甲板、中間甲板和主甲板;(2)立柱及外側(cè)凸起結(jié)構(gòu);(3)浮體及內(nèi)側(cè)凸起結(jié)構(gòu)。上甲板距基線高30.5m,中間甲板距基線高27.5m,上甲板和中間甲板寬49.2m,板厚6mm,縱桁T8×400/12×200mm;主甲板距基線高 24.5m,寬 64.2m,板厚 8mm,箱型縱桁 U2×8×1 800/8×1 600mm,箱型縱桁上縱骨L150×90×9mm;上甲板、中間甲板和主甲板縱骨L100×75×6mm,縱骨間距600mm,構(gòu)件跨距2.4m;縱艙壁為橫骨架式,板厚8mm。立柱直徑9.8m,高18.5m,外側(cè)凸起結(jié)構(gòu)寬2.6m,立柱和外側(cè)凸起結(jié)構(gòu)為橫骨架式,板厚8mm,縱向梁T10×340/12×200mm。浮體和內(nèi)側(cè)凸起結(jié)構(gòu)型深6m,浮體寬15m,內(nèi)側(cè)凸起結(jié)構(gòu)寬2.7m,板厚12mm,縱骨L150×90×12mm,縱骨間距600mm。材料屈服強(qiáng)度355MPa,彈性模量 206GPa。

    圖1 目標(biāo)平臺(tái)橫剖面Fig.1 Cross section of target platform

    圖2 目標(biāo)平臺(tái)總縱彎矩—曲率關(guān)系Fig.2 Longitudinal moment to curvature relationship of target platform

    采用解析方法和簡(jiǎn)化方法對(duì)目標(biāo)平臺(tái)總縱極限強(qiáng)度進(jìn)行了比較研究。簡(jiǎn)化方法將平臺(tái)共劃分519個(gè)單元,圖2給出了簡(jiǎn)化方法得到的平臺(tái)總縱彎矩—曲率關(guān)系曲線,由簡(jiǎn)化方法得到的平臺(tái)中拱極限彎矩為11 728MNm,中垂極限彎矩為-9 937MNm,平臺(tái)中垂和中拱極限彎矩絕對(duì)值比值為0.847。解析方法共劃分92個(gè)單元,由解析方法得到的平臺(tái)中拱極限彎矩為12 221MNm,中垂極限彎矩為-9 848MNm,中垂和中拱極限彎矩絕對(duì)值比值為0.806。解析方法與簡(jiǎn)化方法的計(jì)算結(jié)果比較分析如表1所示,兩種方法的中垂極限彎矩非常吻合,解析方法的中拱極限彎矩略大。由于解析方法和簡(jiǎn)化方法的計(jì)算模型和計(jì)算方法不同,兩種方法的比較分析驗(yàn)證了目標(biāo)平臺(tái)極限強(qiáng)度計(jì)算結(jié)果的可靠性。

    表1 目標(biāo)平臺(tái)總縱極限強(qiáng)度比較分析Tab.1 Comparative study on ultimate strength of target platform

    3 平臺(tái)波浪載荷預(yù)報(bào)

    半潛式平臺(tái)有遷航、作業(yè)和自存三種工況,靜水彎矩的變化主要取決于平臺(tái)自身重量及其縱向分布。平臺(tái)總重量包括空船重量、可移動(dòng)重量(包括可變甲板載荷和消耗品)和壓載,其中可移動(dòng)重量和壓載隨不同的載荷狀態(tài)發(fā)生變化。只要知道了平臺(tái)的重量分布,就可以計(jì)算平臺(tái)的靜水彎矩。而半潛式平臺(tái)的波浪載荷是最為復(fù)雜和關(guān)鍵,合理地預(yù)報(bào)波浪載荷是指導(dǎo)海洋工程設(shè)計(jì)和進(jìn)行結(jié)構(gòu)強(qiáng)度安全性評(píng)估的基礎(chǔ)。

    平臺(tái)波浪載荷的計(jì)算主要基于三維線性勢(shì)流理論和Morison方程,對(duì)于平臺(tái)上的較大尺寸構(gòu)件(浮筒、立柱),由于波浪在結(jié)構(gòu)周圍的繞射效應(yīng)以及浮體運(yùn)動(dòng)產(chǎn)生的輻射效應(yīng)不能忽略,繞射理論是預(yù)報(bào)平臺(tái)主體載荷比較合適的方法;而對(duì)于小尺寸構(gòu)件(如撐桿、系泊系統(tǒng)等),結(jié)構(gòu)尺寸相對(duì)于入射波長(zhǎng)是個(gè)小量,它的存在對(duì)流體運(yùn)動(dòng)的影響可以忽略,采用半經(jīng)驗(yàn)的Morison方程來計(jì)算環(huán)境載荷方便快捷。在頻域內(nèi)將Morison力進(jìn)行線性化之后與平臺(tái)主體波浪載荷疊加,可通過迭代計(jì)算使運(yùn)動(dòng)達(dá)到穩(wěn)態(tài)。這時(shí)對(duì)浮體表面壓力積分,并計(jì)入平臺(tái)加速度引起的慣性力,便可獲得平臺(tái)剖面的特征載荷響應(yīng)。

    圖3 目標(biāo)平臺(tái)垂向波浪彎矩長(zhǎng)期預(yù)報(bào)Fig.3 Long-term value of vertical wave moment of target platform

    在極值波浪載荷預(yù)報(bào)中僅關(guān)心平臺(tái)生存工況下的載荷響應(yīng),生存工況下目標(biāo)平臺(tái)吃水12.19m,排水量 35 744t,重心位置(0,0,17.58m),質(zhì)量矩(1.96×1010,3.36×1010,4.37×1010kg·m2)。選擇 DNV Sesam軟件計(jì)算目標(biāo)平臺(tái)波浪載荷,應(yīng)用GeinE模塊建立平臺(tái)水動(dòng)力模型(包括平臺(tái)主體濕表面模型和桿件的Morison模型)以及結(jié)構(gòu)質(zhì)量模型。將建好的各模型按步驟順序依次導(dǎo)入Wadam模塊中進(jìn)行水動(dòng)力運(yùn)動(dòng)和波浪載荷傳遞函數(shù)(RAO)的計(jì)算,計(jì)算結(jié)束后,應(yīng)用Postresp后處理系統(tǒng)進(jìn)行載荷的短期和長(zhǎng)期預(yù)報(bào),并輸出計(jì)算結(jié)果。長(zhǎng)期預(yù)報(bào)分別采用中國(guó)南海,IACS推薦的北大西洋(以下稱NA-1C)長(zhǎng)期海況資料,DNV推薦的全球平均長(zhǎng)期海況資料(以下稱DNV-WW)共三種長(zhǎng)期海況資料情況。長(zhǎng)期超越概率從10-2~10-10變化,垂向波浪彎矩長(zhǎng)期預(yù)報(bào)結(jié)果如圖3和表2。經(jīng)計(jì)算初步推斷,百年一遇極值大約對(duì)應(yīng)10-8.60的超越概率[10]。

    表2 目標(biāo)平臺(tái)垂向波浪彎矩長(zhǎng)期預(yù)報(bào)(MNm)Tab.2 Long-term value of vertical wave moment of target platform(MNm)

    4 平臺(tái)結(jié)構(gòu)可靠性分析

    由于承載能力和環(huán)境載荷具有明顯的隨機(jī)性,直接以平臺(tái)極限強(qiáng)度Mu為破壞準(zhǔn)則,確立的平臺(tái)極限狀態(tài)方程為

    式中,Mu,Ms和Mw分別為平臺(tái)極限強(qiáng)度、靜水彎矩和波浪彎矩隨機(jī)變量;xu,xs和xw分別為表示平臺(tái)極限強(qiáng)度、靜水彎矩和波浪彎矩分析模型不確定性的隨機(jī)變量。

    平臺(tái)極限強(qiáng)度可表示為

    式中,ti,σ0i和Ei分別為第i個(gè)構(gòu)件的板厚、材料屈服強(qiáng)度和彈性模量。板厚、材料屈服強(qiáng)度和彈性模量的概率分布和變異系數(shù)(COV)如表3所示[11],材料屈服強(qiáng)度的均值取規(guī)范最小值,這樣可能有15%的保守;板厚和彈性模量的均值取名義值。注意到(2)式的平臺(tái)極限強(qiáng)度實(shí)際上是隨機(jī)變量的隱函數(shù)。雖然可靠性計(jì)算的數(shù)值模擬法可以處理帶有隱函數(shù)的極限狀態(tài)方程,但重分析次數(shù)可能過多。崔維成[12]給出改進(jìn)的Rosenbluthe法,可以較方便地求解隨機(jī)函數(shù)的統(tǒng)計(jì)特征,并能達(dá)到二階精度。平臺(tái)極限強(qiáng)度為正態(tài)分布,為了計(jì)及平臺(tái)極限強(qiáng)度模型的不確定性,引入正態(tài)變量xu,如表4-5所示。

    表3 與強(qiáng)度不確定性有關(guān)的隨機(jī)變量Tab.3 Random variable related to strength uncertainty

    表4 強(qiáng)度和載荷模型不確定性變量Tab.4 Uncertainty variable of strength and load model

    表5 平臺(tái)強(qiáng)度和載荷隨機(jī)變量Tab.5 Random variable of platform strength and load

    平臺(tái)重量分布是隨時(shí)間變化的,通常以正態(tài)分布隨機(jī)變量描述平臺(tái)靜水彎矩。此外,由于計(jì)算方法和輸入數(shù)據(jù)的不確定性,理論預(yù)報(bào)結(jié)果與實(shí)際靜水彎矩總有一定的差異,可引入正態(tài)分布變量xs來計(jì)及靜水彎矩的模型不確定性,如表4-5所示。

    在波浪彎矩的線性長(zhǎng)期預(yù)報(bào)方法中引用了一些基本物理假定:波浪短期過程為窄帶分布;波浪周期為相互統(tǒng)計(jì)獨(dú)立。在以上假定條件下,應(yīng)用線性切片理論計(jì)算傳遞函數(shù),從而獲得波浪長(zhǎng)期分布以及極值分布,就沒有多大的困難。如果假定波浪周期相互獨(dú)立,則波浪彎矩極值分布可用泊松分布來模擬,所獲得的分布快速收斂于Gumbel(極值I型)分布。由于在波浪彎矩總不確定性中模型不確定性占較大部分,應(yīng)引入正態(tài)分布變量xw來計(jì)及波浪彎矩分析模型的不確定性,如表4-5所示。

    由于數(shù)值模擬法對(duì)積分維數(shù)、極限狀態(tài)函數(shù)和概率密度函數(shù)沒有任何限制,因此數(shù)值模擬法在現(xiàn)代工程結(jié)構(gòu)中得到廣泛應(yīng)用。注意到拉丁超立方抽樣具有用較少的模擬次數(shù)獲得期望的失效概率的優(yōu)點(diǎn),但它不能減小估計(jì)值的方差,文獻(xiàn)[13]將拉丁超立方抽樣法與條件期望和對(duì)偶變數(shù)方差減縮技術(shù)有機(jī)地結(jié)合,給出結(jié)構(gòu)可靠性分析的組合方法,比較研究結(jié)果表明組合方法既可得到穩(wěn)定的失效概率估計(jì)值,也可以獲得較大的方差減縮,并能針對(duì)一系列不同的極限狀態(tài)函數(shù)。

    本文采用組合方法得到的目標(biāo)平臺(tái)中垂?fàn)顟B(tài)失效概率Pf=1.064×10-3,安全指數(shù)β=3.072。基于目標(biāo)平臺(tái)可靠性分析結(jié)果和相關(guān)文獻(xiàn)的總結(jié),本文推薦半潛式平臺(tái)極限強(qiáng)度目標(biāo)可靠度取β0=3.09,對(duì)應(yīng)的失效概率Pf=1.0×10-3?;诎霛撌狡脚_(tái)極限強(qiáng)度目標(biāo)可靠度,可進(jìn)行半潛式平臺(tái)極限強(qiáng)度可靠性設(shè)計(jì),從而獲得滿足目標(biāo)可靠度要求的平臺(tái)結(jié)構(gòu)設(shè)計(jì)。

    5 結(jié) 論

    本文采用解析方法和簡(jiǎn)化方法對(duì)半潛式平臺(tái)總縱極限強(qiáng)度進(jìn)行了比較研究,基于三維勢(shì)流理論和Morison方程對(duì)生存海況下的半潛式平臺(tái)垂向波浪彎矩進(jìn)行了預(yù)報(bào),基于極限強(qiáng)度建立了半潛式平臺(tái)可靠性計(jì)算模型,并對(duì)目標(biāo)平臺(tái)結(jié)構(gòu)可靠性和目標(biāo)可靠度進(jìn)行了分析。研究結(jié)果表明,簡(jiǎn)化方法和解析方法可用于半潛式平臺(tái)極限強(qiáng)度計(jì)算,半潛式平臺(tái)中垂?fàn)顟B(tài)結(jié)構(gòu)安全性值得關(guān)注,可靠性分析結(jié)果可用于指導(dǎo)半潛式平臺(tái)結(jié)構(gòu)設(shè)計(jì)。

    [1]Qi E R,Cui W C.Designed-oriented methods of ultimate hull girder strength[C].OMAE,2006,92506.

    [2]ISSC.Report of special task committee VI.2(Ultimate hull girder strength)[C]//Proc of 14th ISSC,2000.Nagasaki,2000,2:321-391.

    [3]ISSC.Report of committee III.1(Ultimate strength)[C]//Proc of 16th ISSC,2006.Southampton,2006,1:369-458.

    [4]IACS.Common structural rules for double hull oil tankers[S].2006.

    [5]IACS.Common structural rules for bulk carriers[S].2006.

    [6]Qi E R,Cui W C,Wan Z Q.Comparative study of ultimate hull girder strength of large double hull tankers[J].Marine Structures,2005,18(3):227-249.

    [7]祁恩榮,張曉杰,滕 蓓,陳小平,蔣彩霞.大型液化天然氣船船體極限強(qiáng)度研究[J].船舶力學(xué),2010,14(1-2):66-73.

    [8]張劍波,曾???肖 熙.半潛式平臺(tái)的極限強(qiáng)度分析研究[J].中國(guó)海洋平臺(tái),2005,20(3):19-22.

    [9]Qi E R,Cui W C.Analytical method of ultimate strength of intact and damaged ship hulls[J].Ships and Offshore Structures,2006,1(2):1-12.

    [10]吳東偉,顧學(xué)康,祁恩榮.半潛式平臺(tái)全壽命期極值載荷的確定[J].艦船科學(xué)技術(shù),2011,33(259)(待發(fā)表).

    [11]Mansour A E,Hovem L.Probability-based ship structural safety analysis[J].Journal of Ship Research,1994,38:329-339.

    [12]崔維成,徐向東,邱 強(qiáng).一種快速計(jì)算隨機(jī)變量函數(shù)均值與標(biāo)準(zhǔn)差的新方法[J].船舶力學(xué),1997,1(6):50-60.

    [13]Qi E R,Cui W C.Latin hypercube sampling in ultimate strength reliability of ship hull girder[J].Journal of Ship Mechanics,2002,6(3):41-50.

    猜你喜歡
    方法
    中醫(yī)特有的急救方法
    中老年保健(2021年9期)2021-08-24 03:52:04
    高中數(shù)學(xué)教學(xué)改革的方法
    化學(xué)反應(yīng)多變幻 “虛擬”方法幫大忙
    變快的方法
    兒童繪本(2020年5期)2020-04-07 17:46:30
    學(xué)習(xí)方法
    可能是方法不對(duì)
    用對(duì)方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    最有效的簡(jiǎn)單方法
    山東青年(2016年1期)2016-02-28 14:25:23
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    亚洲真实伦在线观看| 国产精品久久视频播放| 久久精品国产亚洲网站| 中文字幕av在线有码专区| 欧美最新免费一区二区三区| 男人舔女人下体高潮全视频| 一夜夜www| 国产成人aa在线观看| 男女那种视频在线观看| 成人特级黄色片久久久久久久| 最新中文字幕久久久久| 免费搜索国产男女视频| 中文字幕免费在线视频6| 国产精品亚洲美女久久久| 欧美+亚洲+日韩+国产| 长腿黑丝高跟| 男女视频在线观看网站免费| 男人和女人高潮做爰伦理| 三级毛片av免费| 韩国av在线不卡| 国产激情偷乱视频一区二区| 免费在线观看影片大全网站| 精品午夜福利在线看| 欧美xxxx黑人xx丫x性爽| 免费在线观看成人毛片| 18+在线观看网站| 亚洲成人av在线免费| 久久精品国产自在天天线| 嫩草影院精品99| 麻豆国产av国片精品| 国产高清三级在线| 日本黄色视频三级网站网址| 在线免费观看的www视频| 久久午夜福利片| 亚洲久久久久久中文字幕| 久久鲁丝午夜福利片| 亚洲精品日韩av片在线观看| 综合色丁香网| 欧美中文日本在线观看视频| 成人高潮视频无遮挡免费网站| 2021天堂中文幕一二区在线观| 97在线视频观看| 1024手机看黄色片| 麻豆一二三区av精品| 一进一出抽搐gif免费好疼| 黑人高潮一二区| 高清毛片免费观看视频网站| 丝袜喷水一区| 国产女主播在线喷水免费视频网站 | 亚洲国产精品成人久久小说 | 亚洲成人久久性| 一进一出好大好爽视频| 国产成人影院久久av| 别揉我奶头 嗯啊视频| 伦理电影大哥的女人| 久久精品国产自在天天线| 亚洲av熟女| 成人二区视频| 欧美人与善性xxx| 三级国产精品欧美在线观看| 日日啪夜夜撸| 久久精品国产自在天天线| 变态另类丝袜制服| 99在线视频只有这里精品首页| 国产成人a区在线观看| 全区人妻精品视频| 日本精品一区二区三区蜜桃| 在线免费十八禁| 一级毛片aaaaaa免费看小| 99久久久亚洲精品蜜臀av| 亚洲内射少妇av| 国产精品一区www在线观看| 一级毛片aaaaaa免费看小| 91久久精品电影网| 久久人人爽人人爽人人片va| 亚洲欧美清纯卡通| 欧美色欧美亚洲另类二区| 免费人成在线观看视频色| 中国美白少妇内射xxxbb| 久久热精品热| av在线天堂中文字幕| 女人被狂操c到高潮| 特级一级黄色大片| 热99在线观看视频| 欧美色视频一区免费| 啦啦啦韩国在线观看视频| 午夜福利高清视频| 亚洲精品国产成人久久av| 一级黄片播放器| 尾随美女入室| 久久亚洲精品不卡| 男人和女人高潮做爰伦理| 成人欧美大片| 国产午夜精品久久久久久一区二区三区 | 国产真实乱freesex| 亚洲av中文字字幕乱码综合| 欧美日韩在线观看h| av在线播放精品| 成年av动漫网址| av福利片在线观看| 少妇的逼好多水| 秋霞在线观看毛片| 欧美+日韩+精品| 国产精品福利在线免费观看| 天天躁日日操中文字幕| 成人综合一区亚洲| 日韩av不卡免费在线播放| 我要看日韩黄色一级片| 午夜精品国产一区二区电影 | 97热精品久久久久久| 国产aⅴ精品一区二区三区波| 亚洲精品一卡2卡三卡4卡5卡| 日韩中字成人| 99riav亚洲国产免费| 一本精品99久久精品77| 欧美一区二区精品小视频在线| 亚洲无线在线观看| 欧美丝袜亚洲另类| 国产精品一区二区三区四区免费观看 | 欧美色视频一区免费| 看黄色毛片网站| 午夜福利18| 成人午夜高清在线视频| 午夜久久久久精精品| 精品一区二区三区人妻视频| 成人特级av手机在线观看| 成年免费大片在线观看| 欧美日韩综合久久久久久| 亚洲av免费在线观看| av女优亚洲男人天堂| 人人妻人人看人人澡| 亚洲美女黄片视频| 欧美在线一区亚洲| 日本一本二区三区精品| 国产精品久久久久久亚洲av鲁大| 岛国在线免费视频观看| 国产又黄又爽又无遮挡在线| 久久精品人妻少妇| 夜夜爽天天搞| 国产伦精品一区二区三区四那| 日韩av不卡免费在线播放| 免费观看人在逋| 国产精品永久免费网站| 欧美又色又爽又黄视频| 天天躁日日操中文字幕| 精品久久久久久久人妻蜜臀av| 欧美日韩在线观看h| 免费人成视频x8x8入口观看| 午夜精品一区二区三区免费看| 又粗又爽又猛毛片免费看| 欧美日韩乱码在线| 亚洲在线观看片| 此物有八面人人有两片| 午夜福利在线观看吧| 精品一区二区三区av网在线观看| 国产av在哪里看| 岛国在线免费视频观看| 亚洲在线自拍视频| 国产真实伦视频高清在线观看| 久久天躁狠狠躁夜夜2o2o| 真实男女啪啪啪动态图| 村上凉子中文字幕在线| 男人的好看免费观看在线视频| 成人欧美大片| 国产伦精品一区二区三区视频9| 国产亚洲av嫩草精品影院| 麻豆一二三区av精品| 18禁黄网站禁片免费观看直播| 午夜福利在线观看吧| 成人一区二区视频在线观看| 一区福利在线观看| av女优亚洲男人天堂| 国产精品亚洲美女久久久| 婷婷六月久久综合丁香| 国产黄色视频一区二区在线观看 | 中文字幕久久专区| 3wmmmm亚洲av在线观看| 日韩 亚洲 欧美在线| 一个人看视频在线观看www免费| 香蕉av资源在线| 国产成人aa在线观看| 人妻制服诱惑在线中文字幕| 久久久久久久久久黄片| 一级av片app| 九九在线视频观看精品| 亚洲av一区综合| 日韩欧美国产在线观看| 嫩草影视91久久| 69人妻影院| 十八禁国产超污无遮挡网站| 久久国产乱子免费精品| 一本一本综合久久| 人妻制服诱惑在线中文字幕| 色哟哟·www| 人妻久久中文字幕网| 精品一区二区三区视频在线观看免费| 舔av片在线| 直男gayav资源| 日韩欧美 国产精品| 国产精品人妻久久久影院| 亚洲精品国产av成人精品 | 国产黄色小视频在线观看| 中文在线观看免费www的网站| 男人舔奶头视频| 男女视频在线观看网站免费| 久久亚洲国产成人精品v| 国产精品av视频在线免费观看| 给我免费播放毛片高清在线观看| 久久久久国产网址| 露出奶头的视频| 精品久久久久久久久久久久久| 在线观看av片永久免费下载| 男插女下体视频免费在线播放| 国产成人一区二区在线| 国产精品三级大全| 又黄又爽又刺激的免费视频.| av中文乱码字幕在线| 欧美成人一区二区免费高清观看| 伦精品一区二区三区| 色吧在线观看| 亚洲av不卡在线观看| 精品一区二区三区视频在线观看免费| 少妇猛男粗大的猛烈进出视频 | 此物有八面人人有两片| 久久九九热精品免费| 欧美日韩国产亚洲二区| 国产黄a三级三级三级人| 校园春色视频在线观看| 国产一级毛片七仙女欲春2| 免费大片18禁| 精品久久久久久久人妻蜜臀av| 99久久精品热视频| 精品午夜福利视频在线观看一区| av中文乱码字幕在线| 男插女下体视频免费在线播放| 中国美白少妇内射xxxbb| 97超视频在线观看视频| 色av中文字幕| 午夜福利18| 99久国产av精品国产电影| 成人高潮视频无遮挡免费网站| 精品久久久久久久久av| 国产在线精品亚洲第一网站| 精品一区二区免费观看| 日本撒尿小便嘘嘘汇集6| 你懂的网址亚洲精品在线观看 | 综合色av麻豆| 精品一区二区三区av网在线观看| 国产精品久久久久久精品电影| 成年免费大片在线观看| 日韩av不卡免费在线播放| 亚洲精品久久国产高清桃花| 亚洲最大成人中文| 18禁在线播放成人免费| 真人做人爱边吃奶动态| 亚洲18禁久久av| 校园春色视频在线观看| 亚洲成人精品中文字幕电影| 亚洲自偷自拍三级| 蜜桃亚洲精品一区二区三区| 日本免费一区二区三区高清不卡| 久久国内精品自在自线图片| 亚洲欧美日韩东京热| 欧美一区二区国产精品久久精品| 国产成人a区在线观看| avwww免费| 最近2019中文字幕mv第一页| 欧美日韩综合久久久久久| 日韩成人av中文字幕在线观看 | 在线观看美女被高潮喷水网站| 亚洲欧美日韩高清在线视频| 99热网站在线观看| 美女高潮的动态| 国产老妇女一区| 精品久久久久久久久av| 日韩亚洲欧美综合| 欧美性猛交╳xxx乱大交人| 直男gayav资源| 亚洲人成网站在线观看播放| 精品人妻视频免费看| 国产精品精品国产色婷婷| 欧美3d第一页| 国产午夜精品论理片| 男人舔女人下体高潮全视频| 俺也久久电影网| 成人精品一区二区免费| 97碰自拍视频| 18禁在线播放成人免费| 欧美另类亚洲清纯唯美| 在线a可以看的网站| 日韩欧美一区二区三区在线观看| 性插视频无遮挡在线免费观看| 尤物成人国产欧美一区二区三区| 国产视频一区二区在线看| 国产精品,欧美在线| 亚洲欧美日韩高清专用| 午夜福利成人在线免费观看| 国产探花极品一区二区| 美女免费视频网站| 岛国在线免费视频观看| 免费看av在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清激情床上av| 一边摸一边抽搐一进一小说| 九九爱精品视频在线观看| 欧美最新免费一区二区三区| 热99在线观看视频| 国产精品久久视频播放| 国产综合懂色| 永久网站在线| 精品久久国产蜜桃| 成人性生交大片免费视频hd| 淫妇啪啪啪对白视频| 免费观看在线日韩| 亚洲一级一片aⅴ在线观看| 男插女下体视频免费在线播放| 国产精品一区二区免费欧美| 亚洲精品一卡2卡三卡4卡5卡| 国产蜜桃级精品一区二区三区| 插阴视频在线观看视频| 观看免费一级毛片| 在线天堂最新版资源| 少妇的逼好多水| 国产中年淑女户外野战色| 麻豆一二三区av精品| 亚洲乱码一区二区免费版| 美女xxoo啪啪120秒动态图| 麻豆av噜噜一区二区三区| 特级一级黄色大片| 最新中文字幕久久久久| 欧美成人一区二区免费高清观看| 亚洲精华国产精华液的使用体验 | 亚洲人成网站在线播放欧美日韩| 菩萨蛮人人尽说江南好唐韦庄 | 国产午夜精品论理片| 欧美bdsm另类| 日韩欧美精品免费久久| 老女人水多毛片| 国产成人精品久久久久久| 欧美精品国产亚洲| 欧美xxxx性猛交bbbb| 国内精品一区二区在线观看| 日日撸夜夜添| 久久人人爽人人片av| 日本熟妇午夜| 香蕉av资源在线| 能在线免费观看的黄片| 色5月婷婷丁香| avwww免费| 亚洲欧美日韩高清专用| 小蜜桃在线观看免费完整版高清| 色在线成人网| 夜夜爽天天搞| 久久人妻av系列| 婷婷精品国产亚洲av在线| 日本爱情动作片www.在线观看 | 国产精品久久久久久精品电影| 亚洲国产精品成人久久小说 | 国产精品乱码一区二三区的特点| av免费在线看不卡| 日韩 亚洲 欧美在线| 亚洲色图av天堂| 97超视频在线观看视频| 最近2019中文字幕mv第一页| 久久鲁丝午夜福利片| 国国产精品蜜臀av免费| 成人鲁丝片一二三区免费| av福利片在线观看| 亚洲精品粉嫩美女一区| 国产精品久久久久久久久免| 亚洲中文字幕日韩| a级毛片a级免费在线| 特大巨黑吊av在线直播| 99久久久亚洲精品蜜臀av| 欧美日韩乱码在线| 国产91av在线免费观看| 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站| 亚洲不卡免费看| 我的老师免费观看完整版| 日韩欧美免费精品| 久久久久久久久久久丰满| 一进一出抽搐gif免费好疼| 一区二区三区四区激情视频 | 成人欧美大片| 国产精品无大码| 午夜老司机福利剧场| 亚洲在线观看片| 深夜a级毛片| 一本久久中文字幕| 尾随美女入室| 97超视频在线观看视频| 欧美色欧美亚洲另类二区| 婷婷精品国产亚洲av| 国产精品三级大全| 精品日产1卡2卡| 日韩欧美 国产精品| а√天堂www在线а√下载| 久久久久久国产a免费观看| 国产av一区在线观看免费| 12—13女人毛片做爰片一| а√天堂www在线а√下载| 国产色婷婷99| 日韩精品青青久久久久久| 91久久精品国产一区二区成人| 老女人水多毛片| 午夜影院日韩av| 春色校园在线视频观看| 在线国产一区二区在线| 搡女人真爽免费视频火全软件 | 免费一级毛片在线播放高清视频| 深夜a级毛片| 午夜福利视频1000在线观看| 久久久久久久久久成人| 国产成年人精品一区二区| 在线免费观看不下载黄p国产| 熟妇人妻久久中文字幕3abv| 美女黄网站色视频| 久久久国产成人免费| 亚洲精品一卡2卡三卡4卡5卡| 国国产精品蜜臀av免费| 一级av片app| 亚洲成人久久性| 午夜亚洲福利在线播放| 成人一区二区视频在线观看| 久久精品久久久久久噜噜老黄 | 一区二区三区高清视频在线| 性插视频无遮挡在线免费观看| 亚洲精华国产精华液的使用体验 | 色噜噜av男人的天堂激情| 最后的刺客免费高清国语| 中国美白少妇内射xxxbb| 精品午夜福利在线看| 国产视频内射| 女同久久另类99精品国产91| 狂野欧美激情性xxxx在线观看| 99热这里只有精品一区| 国产午夜福利久久久久久| 国产黄色小视频在线观看| 亚洲av中文av极速乱| 一a级毛片在线观看| 日韩 亚洲 欧美在线| 国产黄片美女视频| 乱人视频在线观看| 黄色欧美视频在线观看| 午夜影院日韩av| 给我免费播放毛片高清在线观看| 久久久久国产网址| 久久久久久久午夜电影| 精品99又大又爽又粗少妇毛片| 两个人的视频大全免费| 搡老妇女老女人老熟妇| 精品久久久噜噜| 久久精品久久久久久噜噜老黄 | 22中文网久久字幕| 99久久无色码亚洲精品果冻| 99热精品在线国产| 国内揄拍国产精品人妻在线| 最近的中文字幕免费完整| 插逼视频在线观看| 最近视频中文字幕2019在线8| 亚洲欧美精品综合久久99| 永久网站在线| 日日摸夜夜添夜夜爱| 国产激情偷乱视频一区二区| 男女视频在线观看网站免费| 韩国av在线不卡| 全区人妻精品视频| 午夜亚洲福利在线播放| 又黄又爽又刺激的免费视频.| 欧美bdsm另类| 国产人妻一区二区三区在| 亚洲色图av天堂| 欧美激情在线99| 国内精品一区二区在线观看| 少妇的逼好多水| 久久久久精品国产欧美久久久| 97碰自拍视频| 一本久久中文字幕| 国内精品一区二区在线观看| 亚洲无线在线观看| 成人美女网站在线观看视频| 婷婷色综合大香蕉| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av| 69人妻影院| 久久久欧美国产精品| 91麻豆精品激情在线观看国产| 在线看三级毛片| 欧美日韩综合久久久久久| 一个人免费在线观看电影| 免费一级毛片在线播放高清视频| 午夜a级毛片| 桃色一区二区三区在线观看| 欧美区成人在线视频| а√天堂www在线а√下载| 久久久a久久爽久久v久久| 久久这里只有精品中国| 欧美色欧美亚洲另类二区| 1000部很黄的大片| 午夜福利在线观看吧| 久久鲁丝午夜福利片| 欧美一区二区亚洲| 久久久久久久久中文| 天堂√8在线中文| 国产色爽女视频免费观看| 久久精品国产清高在天天线| 麻豆av噜噜一区二区三区| 此物有八面人人有两片| 在现免费观看毛片| 国产精品国产高清国产av| 精品久久国产蜜桃| 少妇猛男粗大的猛烈进出视频 | 日韩欧美免费精品| 人人妻人人看人人澡| 在线免费观看的www视频| 亚洲av电影不卡..在线观看| av天堂在线播放| 亚洲最大成人手机在线| 中文亚洲av片在线观看爽| 日日干狠狠操夜夜爽| 91狼人影院| 色哟哟·www| 天堂网av新在线| 色吧在线观看| 亚洲美女搞黄在线观看 | 一进一出抽搐动态| 寂寞人妻少妇视频99o| 最近手机中文字幕大全| 婷婷精品国产亚洲av| 少妇猛男粗大的猛烈进出视频 | 亚洲国产色片| 免费观看在线日韩| 18禁裸乳无遮挡免费网站照片| 三级经典国产精品| 免费看美女性在线毛片视频| 国产亚洲精品综合一区在线观看| 一级毛片电影观看 | 日本色播在线视频| 色视频www国产| 国产精品福利在线免费观看| 三级毛片av免费| 1000部很黄的大片| 麻豆国产av国片精品| 永久网站在线| 日韩高清综合在线| 亚洲精品一卡2卡三卡4卡5卡| 欧美性感艳星| 18+在线观看网站| 久久久久久久午夜电影| 天堂av国产一区二区熟女人妻| 床上黄色一级片| 国产欧美日韩精品一区二区| 免费高清视频大片| 国产精品电影一区二区三区| 亚洲av成人av| 成人美女网站在线观看视频| 91麻豆精品激情在线观看国产| 亚洲天堂国产精品一区在线| 亚洲欧美成人综合另类久久久 | 人人妻人人澡欧美一区二区| aaaaa片日本免费| 此物有八面人人有两片| 亚洲人与动物交配视频| 久久久午夜欧美精品| 亚洲成a人片在线一区二区| 自拍偷自拍亚洲精品老妇| 午夜激情欧美在线| 少妇人妻一区二区三区视频| 成年av动漫网址| 国产精品一及| 一本精品99久久精品77| 亚洲精品粉嫩美女一区| 国产亚洲精品av在线| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩一区二区视频在线观看视频在线 | 晚上一个人看的免费电影| 神马国产精品三级电影在线观看| 熟女人妻精品中文字幕| 欧美又色又爽又黄视频| 久久人人爽人人片av| 神马国产精品三级电影在线观看| 美女 人体艺术 gogo| 午夜影院日韩av| 性欧美人与动物交配| 在线观看一区二区三区| 午夜精品国产一区二区电影 | 精品久久久久久久久久久久久| 国产精品一及| 欧美一级a爱片免费观看看| 国产男人的电影天堂91| 国产精品人妻久久久久久| 国产高清有码在线观看视频| 在线观看免费视频日本深夜| 久久午夜亚洲精品久久| 日韩欧美 国产精品| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 国产v大片淫在线免费观看| 国产精品福利在线免费观看| 欧美成人精品欧美一级黄| 婷婷六月久久综合丁香| 观看免费一级毛片| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 综合色av麻豆| 伦精品一区二区三区| 国产精品女同一区二区软件| 久久久久国产精品人妻aⅴ院| 成人午夜高清在线视频| 大香蕉久久网| 性色avwww在线观看| 国产精品一区二区性色av| 久久精品国产亚洲网站| 国内揄拍国产精品人妻在线|