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

    大型FLNG液艙晃蕩和船體耦合運動研究現(xiàn)狀及展望*

    2017-06-21 15:12:33韓旭亮朱小松王俊榮
    中國海上油氣 2017年1期
    關(guān)鍵詞:液艙模型試驗船體

    韓旭亮 謝 彬 朱小松 王俊榮

    (1.中海油研究總院 北京 100028;2.中國石油大學(xué)(北京) 北京 102249)

    大型FLNG液艙晃蕩和船體耦合運動研究現(xiàn)狀及展望*

    韓旭亮1,2謝 彬1朱小松1王俊榮1

    (1.中海油研究總院 北京 100028;2.中國石油大學(xué)(北京) 北京 102249)

    對FLNG應(yīng)用背景進行了介紹,指出了液艙晃蕩和船體耦合運動是FLNG面臨的關(guān)鍵問題之一,并分別從圍護系統(tǒng)選型、模型試驗技術(shù)和數(shù)值模擬技術(shù)方面對液艙晃蕩和船體耦合運動研究現(xiàn)狀進行了總結(jié)分析,探討了試驗技術(shù)研究的特點和數(shù)值模擬技術(shù)的應(yīng)用范圍,并對大型FLNG液艙晃蕩和船體耦合運動下一步的研究方向提出了建議。

    FLNG;液艙晃蕩;耦合運動;圍護系統(tǒng);模型試驗;數(shù)值模擬

    浮式液化天然氣生產(chǎn)儲卸裝置FLNG是一種用于海上天然氣田開發(fā)的浮式生產(chǎn)裝置,通過系泊系統(tǒng)定位于海上,具有開采、處理、液化、儲存和裝卸天然氣的功能,并通過與液化天然氣船搭配使用,實現(xiàn)海上天然氣田的開采和運輸。利用FLNG進行深水氣田開發(fā)結(jié)束了采用管道運輸上岸的單一模式,避免了深水鋪管的難題,既節(jié)約了成本,也為邊際氣田開發(fā)提供了靈活配置、經(jīng)濟有效的方案。隨著天然氣供應(yīng)緊張和相關(guān)應(yīng)用技術(shù)逐步成熟,F(xiàn)LNG概念的工程化已被眾多能源公司所接受,但是目前世界上尚未有FLNG正式投產(chǎn)使用。2011年5月荷蘭Shell公司與三星重工及Technip公司簽訂建造全球第一艘FLNG合同,該Prelude FLNG預(yù)計2017年投產(chǎn),用于澳大利亞西北天然氣田。馬來西亞國家石油公司分別與大宇造船海洋和三星重工簽訂建造PFLNG1和PFLNG2合同,計劃2016年底和2018年分別用于Kanowit和Rotan氣田作業(yè)。國內(nèi)南通惠生重工也在為國外Exmar公司建造年產(chǎn)50萬t的駁船型FLNG。中海油研究總院與多家科研院所及高校合作,從“十一五”期間開始以南海某深水氣田為目標(biāo),針對FLNG裝置開展系統(tǒng)研究,并取得階段性研究成果。同時,我國工業(yè)和信息化部在2015年10月30日正式發(fā)布的《中國制造2025》重點領(lǐng)域技術(shù)線路藍(lán)圖中已明確提出將FLNG等新型海洋油氣資源開發(fā)裝備列為重點產(chǎn)品。

    雖然FLNG船仍處于設(shè)計、建造階段,但是潛在風(fēng)險值得分析,須借鑒LNG運輸船的典型海損事故對其特有的風(fēng)險因素進行分析,即液貨系統(tǒng)故障分析。液艙內(nèi)低溫液體在船體運動過程中對艙壁結(jié)構(gòu)的沖擊作用是導(dǎo)致液體泄漏的重要原因,比較典型的LNG運輸船事故有:由于低溫液體劇烈沖擊液艙,Polar Alaska、Arctic Tokyo和Catalunya Spirit均在液艙的內(nèi)薄膜或角點多處發(fā)生變形或受損,EI Paso公司12.9 萬m3的LNG運輸船在外界波浪運動頻率和艙內(nèi)低溫液體晃蕩頻率接近時加重了艙內(nèi)低溫液體晃動,造成了艙壁破損液體泄漏。大型FLNG船不能自航,長期停泊外海,無法避免惡劣海況作用,連續(xù)生產(chǎn)作業(yè)時內(nèi)部液艙會經(jīng)歷各種裝載情況,無法按照LNG運輸船的指導(dǎo)規(guī)定使艙內(nèi)液體低于10%或者高于70%,因此,F(xiàn)LNG船的液艙會比LNG運輸船更易受到晃蕩的影響。FLNG船體運動激勵液艙內(nèi)低溫液體產(chǎn)生晃蕩,同時液艙內(nèi)低溫液體的晃蕩產(chǎn)生的力矩會改變FLNG船體的運動性能,兩者相互作用存在耦合關(guān)系。不同裝載工況的艙內(nèi)低溫液體晃蕩可能與FLNG船體運動發(fā)生共振響應(yīng),從而引起船體劇烈運動并造成較大傷害,降低了船體穩(wěn)定性,進而引起上部甲板液化裝置的失效,極端情況下將引起船體事故。由此可見,液艙晃蕩和船體水動力耦合運動的整體水動力特性是保證FLNG裝置安全運行的重要問題。

    本文圍繞大型FLNG液艙晃蕩和船體耦合運動這一具體并迫切需要解決的問題,借助于試驗技術(shù)和數(shù)值模擬技術(shù),對液艙晃蕩和船體耦合運動問題進行總結(jié)分析,以期為我國FLNG船的設(shè)計提供理論依據(jù)和技術(shù)支撐。

    1 圍護系統(tǒng)選型

    圍護系統(tǒng)的主要作用是將超低溫的LNG與船體結(jié)構(gòu)隔離并實現(xiàn)保冷存儲,通過絕熱材料保持所裝載的LNG處于液態(tài),避免并保護船體結(jié)構(gòu)直接承受低溫貨物的影響。目前LNG圍護系統(tǒng)的主要型式為SPB棱柱型液艙、MOSS獨立球型液艙和GTT薄膜型液艙[1]。FLNG裝置的圍護系統(tǒng)除達(dá)到常規(guī)LNG運輸船的要求外,還應(yīng)當(dāng)滿足的條件有:FLNG裝置需要平整、寬闊、大的甲板面積提供足夠的空間來布置生產(chǎn)模塊、管系、設(shè)備等,還要充分考慮每組模塊支墩的數(shù)量、支撐位置、船體變形量與結(jié)構(gòu)加強等;貨艙裝載率不受限制的工況下能抵御惡劣海況下的液艙晃蕩沖擊載荷;可靠的長壽命設(shè)計,可以有效降低停產(chǎn)風(fēng)險;易于作業(yè)現(xiàn)場進行檢驗和維護。目前對于LNG產(chǎn)能大于100萬t/a的FLNG,滿足上述特點且技術(shù)成熟的圍護系統(tǒng)主要有GTT薄膜型液艙和SPB棱型液艙[2]。為了使艙內(nèi)裝載高度不受限制,并降低FLNG裝置圍護系統(tǒng)受到的液艙晃蕩沖擊載荷,GTT薄膜型液艙采用雙排艙設(shè)計。2種FLNG圍護系統(tǒng)的橫剖面如圖1所示。

    從圖1可以看出,GTT薄膜型貨艙的縱向和橫向設(shè)有雙層隔離艙壁,圍護系統(tǒng)區(qū)域的船體強度較好,且中縱隔離艙壁可為上部模塊支墩提供較好的支撐;而SPB型貨艙中間無縱艙壁,橫向為單層,圍護系統(tǒng)區(qū)域的船體強度偏弱,且上部模塊支墩跨距較大,自身的設(shè)計強度要求較高。GTT薄膜型液艙受艙內(nèi)液體晃蕩影響比較嚴(yán)重,船體運動可能會導(dǎo)致艙內(nèi)劇烈晃蕩,從而在艙壁產(chǎn)生非常大的沖擊載荷,但雙排艙設(shè)計可以有效減弱液艙晃蕩沖擊載荷對圍護系統(tǒng)的影響;而SPB型圍護系統(tǒng)不存在晃蕩問題,不需限制載液率。GTT薄膜型圍護系統(tǒng)面臨加熱系統(tǒng)的限制,而SPB型圍護系統(tǒng)過度依賴IHI,絕緣層的安裝存在風(fēng)險??梢姡珿TT薄膜型圍護系統(tǒng)和SPB型圍護系統(tǒng)都滿足FLNG裝置的要求并各有優(yōu)缺點(表1)。結(jié)合建造成本和國內(nèi)自主建造能力等方面綜合考慮,GTT薄膜型圍護系統(tǒng)具有較明顯的優(yōu)勢。因此,建議FLNG裝置采用GTT薄膜型圍護系統(tǒng)。

    圖1 FLNG圍護系統(tǒng)的橫剖面Fig.1 Transverse section of containment system for FLNG

    表1 GTT薄膜型和SPB型貨物圍護系統(tǒng)對比Table 1 Comparison of containment system between GTT and SPB

    2 模型試驗技術(shù)研究

    2.1 液艙晃蕩模型試驗研究

    由于晃蕩現(xiàn)象具有較強的非線性,如波浪破碎、波面翻轉(zhuǎn)和氣泡效應(yīng)等,而這些特殊的物理現(xiàn)象對晃蕩沖擊載荷研究又至關(guān)重要,但目前數(shù)值模擬技術(shù)水平有限。因此,借助模型試驗可以真實復(fù)原液艙晃蕩復(fù)雜水動力現(xiàn)象,并對晃蕩物理現(xiàn)象和艙內(nèi)沖擊載荷進行深入研究。液艙晃蕩模型試驗主要通過大噸位運動平臺來模擬船體運動。早期的模型試驗主要集中于矩形或者球型單自由度規(guī)則激勵的晃蕩試驗研究。隨著試驗技術(shù)的不斷提高,在室內(nèi)利用六自由度運動平臺模擬船體運動變成了現(xiàn)實。表2給出了目前國內(nèi)外擁有先進的多自由度運動液艙晃蕩試驗平臺的對比情況,利用這些平臺都可以精確實現(xiàn)規(guī)則和不規(guī)則的單自由度和多自由度耦合運動的模型試驗。從2012年開始,多所研究機構(gòu)參與了一系列矩形標(biāo)準(zhǔn)試驗[3-4],旨在對比分析相同運動激勵在不同運動平臺下的復(fù)現(xiàn)情況,以及六自由度運動平臺實現(xiàn)方式和精度對晃蕩載荷的影響。

    表2 多自由度運動液艙晃蕩試驗平臺對比Table 2 Comparison of multi DOF motion sloshing test platform

    如何才能使液艙晃蕩模型試驗和原型具有較好的相似性,是目前許多學(xué)者研究的熱點。Bass等[5]通過模型試驗分別討論了流體黏性、液體壓縮率、氣體壓縮率、氣穴效應(yīng)和表面張力等參數(shù)對載荷相似性的影響,結(jié)果表明影響相似處理的主要影響因素是氣體和液體壓縮率,而不是流體黏性和表面張力。Yung等[6]通過典型矩形艙晃蕩試驗給出了沖擊載荷、生成時間、沖擊速度和加速度、力和力矩等模型和實物的相似準(zhǔn)則。Kimmoun等[7]在試驗中發(fā)現(xiàn)壓力變換系數(shù)是λ0.5,而不是Yung等[6]認(rèn)為的λ。多次試驗表明,液艙晃蕩模型試驗除了受Froude系數(shù)主導(dǎo)外,還受到其他參數(shù)影響,諸如時間和空間對艙內(nèi)液體晃蕩壓力的影響、艙內(nèi)氣液密度比的影響以及艙內(nèi)氣液交互作用的影響等都會導(dǎo)致比例尺還原系數(shù)不夠明確,因此模型試驗與全尺度試驗結(jié)果間的差異還有待今后進一步研究。

    由于液艙晃蕩小尺度模型試驗得到的晃蕩載荷存在縮尺比效應(yīng),近年來也有相關(guān)學(xué)者進行了全尺度原型試驗。這類試驗雖然尺度上與實際液艙相同,且能直觀地展現(xiàn)液艙內(nèi)晃蕩的劇烈程度和得到晃蕩產(chǎn)生的沖擊載荷,但是實船中液艙晃蕩是長期時歷過程而不是單次的砰擊作用,因此全尺度原型試驗對液艙運動及真實晃蕩過程還不能完全等效模擬,僅能模擬波浪對艙壁的單次砰擊過程并測量砰擊壓力。由Marin和GTT等單位聯(lián)合進行了全尺度原型試驗Sloshed[8],如圖2所示。該試驗詳細(xì)研究了晃蕩壓力、氣穴、波浪飛濺等特性,并對晃蕩載荷與彈性結(jié)構(gòu)的流固耦合問題也進行了初步研究。通過這些試驗可以發(fā)現(xiàn)局部晃蕩壓強與液體拍擊壁面的形式有關(guān),但是它們對墻體的載荷總體相當(dāng)。

    圖2 液艙晃蕩全尺度原型試驗[8]Fig.2 Full scale prototype test of liquid tank sloshing

    2.2 液艙晃蕩和船體耦合運動試驗研究

    與常規(guī)船舶不同,F(xiàn)LNG船帶有載液艙室,在模型試驗中無法保證所有的相似準(zhǔn)則數(shù)都相同,通常保證Froude數(shù)相等。許多學(xué)者通過模型試驗研究液艙晃蕩和船體耦合運動問題,并對開發(fā)的數(shù)值計算方法進行了計算驗證。Molin等[9]對帶有2個部分充液的矩形液艙的駁船在不規(guī)則波作用下的橫搖運動進行了模型試驗(圖3),重點研究了液艙大幅晃蕩運動和不同艙頂(平板和倒角)距艙內(nèi)液面高度的非線性影響,結(jié)果表明:隨著海況增加,液艙晃蕩非線性效應(yīng)顯著;橫搖運動第2個峰值逐漸減小至變平緩,并向低頻方向移動,但是艙頂和艙內(nèi)液面之間的氣穴并不因為艙頂撞擊數(shù)量和強度而有所不同,這說明第2個峰值的降低可能與內(nèi)部耗散關(guān)系不大,更主要是由于艙內(nèi)液面的非線性特性引起的。但是,這些結(jié)論是否適用于其他自由度、液艙形狀、浮體形狀還需進一步研究。Nam等[10]進行了雙液艙FLNG船耦合運動的模型試驗,研究了不同裝載下液艙晃蕩對船體整體運動響應(yīng)及艙內(nèi)流體壓強的影響(圖4),發(fā)現(xiàn)液艙內(nèi)自由水面非線性對船舶運動響應(yīng)影響較大,且運動幅值和波浪幅值并非呈線性關(guān)系;而艙內(nèi)液體壓力與船型、液艙形狀、波浪頻率和幅值等有關(guān)。中海油研究總院和上海交通大學(xué)合作[11]進行了多液艙FLNG船耦合運動的模型試驗(圖5),研究了固體和液體多種載況下液艙晃蕩和船體耦合運動的響應(yīng)情況,結(jié)果表明液艙晃蕩對FLNG船的橫蕩和橫搖運動影響嚴(yán)重,且低液位裝載比高液位裝載在相同外界激勵下更易引起較嚴(yán)重的晃蕩響應(yīng)。該試驗也證實了液艙內(nèi)液面運動幅值和入射波幅呈非線性關(guān)系,并提出了在液艙內(nèi)安裝浪高儀測量艙內(nèi)液面變化的試驗技術(shù)。

    可見,液艙晃蕩和船體耦合運動水池模型試驗的振動方式由單自由度逐漸趨向于多自由度聯(lián)合耦合激勵,集中研究其非線性問題,但多數(shù)研究對象為液艙晃蕩和剛性船體的耦合運動,忽略了液艙晃蕩對船體結(jié)構(gòu)強度、振動破壞形式等方面的影響。目前試驗研究僅局限于模型試驗,經(jīng)濟成本過高、不受尺度效應(yīng)影響的實船試驗還未開展。

    圖3 帶液艙駁船橫搖運動的模型試驗[9]Fig.3 Model test for roll motion of liquid tank barge

    圖4 雙液艙FLNG船耦合運動的模型試驗[10]Fig.4 Mdoel test of coupling motion of FLNG with double tanks

    圖5 多液艙FLNG船耦合運動的模型試驗[11]Fig.5 Mdoel test of coupling motion of FLNG with multi-tanks

    3 數(shù)值模擬技術(shù)研究

    3.1 液艙晃蕩數(shù)值模擬研究

    液艙晃蕩的數(shù)值模擬方法大體可分為勢流理論方法、CFD (Computational Fluid Dynamics)有網(wǎng)格方法和CFD無網(wǎng)格方法等3種,各方法優(yōu)缺點見表3。勢流理論方法包括基于勢流理論的線性解析方法[12]、非線性復(fù)合模態(tài)方法[13]、非線性邊界元方法[14]和非線性有限元方法[15]等,是在流域邊界上配置一組奇點,通過配置點上的邊界條件求解相應(yīng)奇點上的未知量封閉方程組,求出流域各點處的速度勢和速度。與其他方法相比,勢流理論方法的最大特點是計算速度快,求解精度高,但是難以處理流體黏性效應(yīng),盡管引入了人工黏性阻尼[16],仍無法模擬艙內(nèi)波面發(fā)生翻轉(zhuǎn)和破碎等強非線性現(xiàn)象的晃蕩問題。CFD有網(wǎng)格方法是基于求解不可壓縮Navier-Stokes方程的MAC (Marker-and-cell)[17]、ALE (Arbitrary Lagrangian Eulerian)[18]、VOF (Vloume of Fluid)[19]和Level Set[20]方法等,基于歐拉網(wǎng)格將流場離散可以模擬液艙的劇烈晃蕩,但是生成網(wǎng)格困難且耗時,且液艙晃蕩的非線性現(xiàn)象模擬受網(wǎng)格影響較大,易因幾何大變形而失敗。CFD無網(wǎng)格方法包括無網(wǎng)格的SPH (Smoothed Particle Hydrodynamics)[21]、MLPG (Meshless Local Petrov-Galerkin)[22]和MPS (Moving Particle Semi-implicit)[23]方法等,與前2種方法相比,該方法不依賴網(wǎng)格,適合處理艙內(nèi)液面發(fā)生翻轉(zhuǎn)和破碎等大變形的強非線性問題,但是在計算效率、流場壓力預(yù)報精度[24]和未確定參數(shù)等方面尚需進一步研究。可見,數(shù)值模擬技術(shù)可以較好模擬液艙晃蕩問題,但還不能完全模擬該現(xiàn)象,其精度還有待進一步提高,需要耗費較多的計算時間才能獲取較好的數(shù)值分析結(jié)果。

    表3 液艙晃蕩數(shù)值模擬方法優(yōu)缺點對比Table 3 Advantages and disadvantages comparison of numerical simulation methods for sloshing

    3.2 液艙晃蕩和船體耦合運動數(shù)值模擬研究

    起初有關(guān)液艙晃蕩和船體耦合運動的研究較多是基于頻域勢流理論進行數(shù)值模擬[25-26],頻域理論的前提條件是液艙晃蕩的波動頻率與外界波浪的激勵頻率相同。但是研究發(fā)現(xiàn),液艙晃蕩并非是外界單頻激勵引起的波動,還包括液艙自身固有頻率的振蕩,所導(dǎo)致的船體耦合運動并不滿足頻域理論的單頻假設(shè),很可能存在誤差[27]。雖然艙內(nèi)液體會發(fā)生劇烈晃蕩,但液艙穩(wěn)定后的運動仍將保持線性運動特征,其運動頻率與外界波浪的激勵頻率相同,液艙內(nèi)流體晃蕩所產(chǎn)生的高頻諧波作用力將被系統(tǒng)過濾掉。因此,頻域理論適用于液艙晃蕩和船體耦合運動的穩(wěn)態(tài)響應(yīng)分析[28-29]。

    外流場船體運動引發(fā)艙內(nèi)液體晃蕩,而艙內(nèi)液體晃蕩反過來也會影響船體水動力,兩者相互作用嚴(yán)重影響整體水動力性能。非線性液艙晃蕩的模擬已經(jīng)相當(dāng)耗時,再耦合非線性船體運動后必將大大增加計算負(fù)擔(dān),所以強非線性劇烈晃蕩和船體耦合運動的時域數(shù)值模擬是一個挑戰(zhàn),需要選擇合適的非線性數(shù)值模型進行簡化。近年來,許多學(xué)者從時域理論角度對液艙晃蕩和船體耦合運動非線性問題進行了不同數(shù)值模擬方法的研究。表4給出了當(dāng)前液艙晃蕩和船體耦合運動數(shù)值模擬方法的研究發(fā)展方向。

    表4 液艙晃蕩和船體耦合運動數(shù)值模擬方法的研究方向Table 4 Research direction for numerical simulation of the liquid sloshing and hull coupling Motion

    從表4可以看出,液艙晃蕩和船體耦合運動數(shù)值模擬方法的發(fā)展方向大致分為4類。第1類,基于勢流理論采用脈沖響應(yīng)函數(shù)(Impulse response function,IRF)的間接時域方法計算船體運動,對于艙內(nèi)液體基于勢流理論采用非線性方法模擬艙內(nèi)液體未發(fā)生翻轉(zhuǎn)或破碎情況的晃蕩問題[30-32],相比CFD方法,其計算效率較高,具有工程使用價值。第2類,基于勢流理論采用間接時域方法計算船體運動,對于艙內(nèi)液體基于CFD采用Fluent[33-34]、OpenFOAM[35-36]或者SPH[37]方法模擬液艙非線性晃蕩。前2類中都涉及到IRF方法,即在頻域中計算船體本身的水動力系數(shù)、穩(wěn)態(tài)水動力及平均漂移力,通過傅里葉變換將頻域水動力系數(shù)轉(zhuǎn)為時域脈沖響應(yīng)函數(shù),進而在時域內(nèi)建立液艙晃蕩和船體耦合運動方程。該方法涉及到較多人工參與和高頻計算不準(zhǔn)確問題,但其計算效率遠(yuǎn)高于直接時域方法,仍然為分析液艙晃蕩和船體耦合運動的實用方法。前2類都集中于研究液艙晃蕩的非線性問題,未考慮船體外流場非線性的影響。實際上,在外界激勵頻率接近液艙晃蕩固有頻率時,即使外界入射波高較小,也會由于共振導(dǎo)致液艙晃蕩和船體產(chǎn)生大幅運動。第3類,基于勢流理論采用時域格林函數(shù)方法,考慮船體濕表面瞬時變化計算外流場船體運動,對于艙內(nèi)液體基于勢流理論采用非線性方法模擬晃蕩[38]。研究發(fā)現(xiàn),液艙晃蕩非線性對耦合運動的影響最大,船體濕表面變化非線性影響次之,且非線性水動力對耦合運動穩(wěn)態(tài)響應(yīng)幅值有較大影響。第4類,在內(nèi)、外流場都是基于CFD方法模擬液艙晃蕩和船體耦合運動[39],但是計算代價較高,效率較低。

    4 研究展望

    由于液艙晃蕩和船體水動力耦合運動的整體水動力特性是保證FLNG正常安全生產(chǎn)運行的關(guān)鍵問題,因此,今后有必要對以下幾方面進行深入研究,以期為研究FLNG液艙晃蕩和船體耦合運動的水動力機理提供更可靠的依據(jù)。

    1) 確定最佳的尺度效應(yīng)。通過對目前液艙晃蕩模型試驗結(jié)果的研究發(fā)現(xiàn),液艙晃蕩模型試驗都會滿足Froude相似,但是比例尺還原系數(shù)不是很明確,還受到諸如時間和空間對艙內(nèi)晃蕩壓力的影響、氣液密度比的影響以及氣液交互作用的影響等,無法獲得實際液艙所受的載荷。因此,還需要尋找并開展遵循不同相似原則的模型試驗,并對試驗結(jié)果進行比較分析,以確定最佳的尺度效應(yīng),實現(xiàn)模型試驗和實尺度試驗的最小差異化。

    2) 建立滿足真實條件的試驗。試驗?zāi)軌蛘鎸崗?fù)原展現(xiàn)艙內(nèi)液體晃蕩復(fù)雜水動力物理現(xiàn)象,得到艙內(nèi)液體晃蕩產(chǎn)生的沖擊載荷。但是一般的試驗很難滿足真實的FLNG液艙晃蕩的要求,只是在尺度上與實際液艙相同,而激勵形式和艙內(nèi)液體運動都與實際情況存在差異,且試驗中未能采用真實的LNG進行試驗,往往不能滿足低溫條件。因此,反映FLNG液艙晃蕩真實條件的試驗還有待今后進一步驗證研究。

    3) 開展液艙晃蕩和船體水彈性響應(yīng)研究。液艙晃蕩和船體耦合運動的振動方式趨于多自由度聯(lián)合耦合激勵,通常仍假設(shè)研究對象為剛體,然而大型FLNG屬于大尺度海洋結(jié)構(gòu)物,存在相關(guān)水彈性問題。船體結(jié)構(gòu)在外界波浪激勵作用下會產(chǎn)生變形,需要考慮液艙晃蕩和彈性大尺度船體相互之間影響;而且艙內(nèi)液體晃蕩沖擊壓力也會使得液艙結(jié)構(gòu)發(fā)生變形,影響艙內(nèi)液體的運動情況。這方面工作難度較大,尚處于空白。

    4) 建立液艙晃蕩和船體耦合運動的理論模型。從液艙晃蕩和船體耦合運動數(shù)值模擬方法的4類研究方向可以看出,數(shù)值模擬方法在該問題中已得到比較廣泛的應(yīng)用,可以較好地模擬和處理液艙晃蕩和船體耦合運動問題。但是目前理論模型主要是基于部分非線性理論,數(shù)值模擬結(jié)果與試驗結(jié)果存在差異,還不能完全真實地模擬求解該問題,數(shù)值分析方法的精度還有待進一步提高。因此,有必要建立液艙劇烈晃蕩和船體大幅運動的強非線性理論模型,使得數(shù)值模擬更接近實際問題。

    5) 提高液艙晃蕩和船體耦合運動的數(shù)值模擬效率。提高計算效率始終是數(shù)值模擬中的重要問題。液艙晃蕩和船體耦合運動的機理和特性非常復(fù)雜,液艙晃蕩的模擬已然相當(dāng)耗時,再將它和船體水動力相互作用,計算量勢必大大增加。因此,還需要簡化或者改進相關(guān)理論模型,提高液艙晃蕩和船體耦合運動的數(shù)值模擬效率。

    [1] 朱小松,謝彬,喻西崇.LNG/LPG液艙晃蕩研究進展綜述[J].中國造船,2013,54(1):229-236.ZHU Xiaosong,XIE Bin,YU Xichong.Research progress of liquid sloshing in LNG/LPG tanks[J].Ship Building of China,2013,54(1):229-236.

    [2] 宋吉衛(wèi).FLNG總體設(shè)計關(guān)鍵技術(shù)研究[J].中國造船,2015,56(2):81-86.SONG Jiwei.Key technology research in FLNG general design[J].Ship Building of China,2015,56(2):81-86.

    [3] LOYSEL T,CHOLLET S,GERVAISE E,et al.Results of the first sloshing model test benchmark[R]∥Proceedings of the 22th International Offshore and Polar Engineering Conference.Rhodes,Greece,2012:398-408.

    [4] 衛(wèi)志軍,張文首,王安良,等.基于SPH方法的二維矩形艙液體晃蕩數(shù)值研究[J].大連理工大學(xué)學(xué)報,2014,54(6):597-603.WEI Zhijun,ZHANG Wenshou,WANG Anliang,et al.Numerical investigation of liquid sloshing in a 2D rectangular tank based on SPH method[J].Journal of Dalian University of Technology,2014,54(6):597-603.

    [5] BASS R L,TRUDELL R W,YOSHIMURA N,et al.Modeling criteria for scaled LNG sloshing experiments[J].Journal of Fluids Engineering,1985,107(2):272-280.

    [6] YUNG T W,DING Z,HE H.LNG sloshing:Characteristics and scaling laws[J].International Journal of Offshore and Polar Engineering,2009,19(4):264-270.

    [7] KIMMOUN O,RATOUIS A,BROSSET L.Sloshing and scaling:experimental study in a wave canal at two different scales[R]∥Proceedings of the 20th International Offshore and Polar Engineering Conference.Beijing,China,2010:33-43.

    [8] KAMINSKI M L,BOGAERT H.Full scale sloshing impact tests[R]∥Proceedings of the 19th International Offshore and Polar Engineering Conference.Osaka,Japan,2009:125-134.

    [9] MOLIN B,REMY F.Effect of roof impacts on coupling between wave response and sloshing in tanks of LNG-carriers[C]∥Proceeding of the ASME 27th International Conference on Offshore Mechanics and Arctic Engineering,2008.

    [10] NAM B W,KIM Y,KIM D W,et al.Experimental and numerical studies on ship motion responses coupled with sloshing in waves[J].Journal of Ship Research,2009,53(2):68-82.

    [11] ZHAO Wenhu,YANG Jianmin,HU Zhiqiang,et al.Experimental investigation of effects of inner-tank sloshing on hydrodynamics of an FLNG system[J].Journal of Hydrodynamics.Ser.B,2012,24(1):107-115.

    [12] FALTINSEN O M,TIMOKHA A N.Sloshing[M].Cambridge:Cambridge University Press,2009.

    [13] FALTINSEN O M,TIMOKHA A N.Sloshing with marine and land-based applications[C]∥ International Conference on Hydrodynamics,Shanghai,China,2010.

    [14] KIM B,SHIN Y.Coupled seakeeping with liquid sloshing in ship tanks[C]∥Proceedings of the ASME 27th International Conference on Offshore Mechanics and Arctic Engineering,2008:247-257.

    [15] WU G X,MA Q W,EATOCK T R.Numerical simulation of sloshing waves in a 3D tank based on a finite element method[J].Applied Ocean Research,1998,20(6):335-337.

    [16] HUANG Shuo,DUAN Wenyang,MA Qingwei.An approximation to energy dissipation in time domain simulation of sloshing waves based on linear potential theory[J].China Ocean Engineering,2011,25(2):189-200.

    [18] 丁峻宏,金允龍,王惠.基于ALE的礦粉貨物液化晃蕩問題并行數(shù)值模擬[J].船舶力學(xué),2015,19(8):927-933.DING Junhong,JIN Yunlong,WANG Hui.ALE-based parallel numerical simulation for sloshing problem of liquefied ore fines cargo[J].Journal of Ship Mechanics,2015,19(8):927-933.

    [19] ZHU R Q,WANG Z D,FANG Z Y.Prediction of pressure induced by large amplitude sloshing in a tank[J].Journal of Ship Mechanics,2004,8(6):71-78.

    [20] ZHANG Yuxin,WANG Decheng,HINO,et al.Comparative study of MPS method and level-set method for sloshing flows[J].Journal of Hydrodynamics,2014,26(4):577-585.

    [21] ZHENG Xing,DUAN Wenyang,MA Qingwei.Comparison of improved meshless interpolation schemes for SPH method and accuracy analysis[J].Journal of Marine Science and Application,2010,9(3):223-230.

    [22] ATLURI S N,ZHU T.A new meshless local Petrov-Ganlerkin (MLPG) approach in computational mechanics [J].Computer Mechanics,1998,22:117-127.

    [23] KOSHIZUKA S,OKA Y.Moving particle semi-implicit for fragmentation of incompressible fluid[J].Nuclear Science and Engineering,1996,123:421-434.

    [24] ZHENG X,MA Q W,DUAN W Y.Incompressible SPH method based on Rankine source solution for violent water wave simulation [J].Journal of Computational Physics,2014,276:291-314.

    [25] NEWMAN J.Effects on vessels with internal tanks[C]∥ Proceedings of the 20th Interactional Workshop on Water Waves and Floating Bodies,Spitsbergen,Norway,2005.

    [26] BOOKI K,YUNG S.Coupled seakeeping with liquid sloshing in ship tanks[C]∥Proceedings of the ASME 27th Interactional Conference on Offshore Mechanics and Arctic Engineering,2008:165-175.

    [27] FALTINSEN O M.A nonlinear theory of sloshing in rectangular tanks[J].Journal of Ship Research,1974,18(4):224-241.

    [28] 徐海霞.液艙晃蕩對船體運動影響的二維與三維頻域線性分析[D].哈爾濱:哈爾濱工程大學(xué),2009.XU Haixia.2-D and 3-D analysis of sloshing effect on ship motions in frequency domain[D].Harbin:Harbin Engineering University,2009.

    [29] 駱陽,朱仁慶,劉永濤.FPSO與運輸船旁靠時液艙晃蕩與船舶運動耦合效應(yīng)分析[J].江蘇科技大學(xué)學(xué)報(自然科學(xué)版),2015,29(4):307-316.LUO Yang,ZHU Renqing,LIU Yongtao.Analysis of coupling effects of sloshing and ship motion for FPSO and transport ship in side-by-side case[J].Journal of Jiangsu University of Science and Technology(Natural Science Edition),2015,29(4):307-316.

    [30] HUANG S,DUAN W Y,ZHU X.Time-domain simulation of tank sloshing pressure and experimental validation[J].Journal of Hydrodynamics.Ser.B,2010,22(5):556-563.

    [31] ZHAO W H,YANG J M,HU Z Q,et al.Hydrodynamics of a 2D vessel including internal sloshing flows[J].Ocean Engineering,2014,84:45-53.

    [32] GOU Y,KIM Y,KIM T Y.A numerical study on coupling between ship motions and sloshing in frequency and time domain[C]∥Proceedings of the 21th International Offshore and Polar Engineering Conference.Hawaii,USA,2011:158-164.

    [33] 文攀.船舶與船舶液艙晃蕩耦合運動分析[D].武漢:華中科技大學(xué),2013.WEN Pan.Analysis on coupling effects of ship motion and tank sloshing[D].Wuhan:Huazhong University of Science and Technology,2013.

    [34] 操戈,李旭,張詠鷗,等.FPSO液艙晃蕩與船舶時域耦合運動數(shù)值模擬[J].中國艦船研究,2015,10(1):88-96.CAO Ge,LI Xu,ZHANG Yongou,et al.Numerical simulation of ship motion coupled with tank sloshing for FPSO[J].Chinese Journal of Ship Research,2015,10(1):88-96.

    [35] 李裕龍,朱仁傳,繆國平,等.基于OpenFOAM的船舶與液艙流體晃蕩在波浪中時域耦合運動的數(shù)值模擬[J].船舶力學(xué),2012,16(7):750-758.LI Yulong,ZHU Renchuan,MIAO Guoping,et al.Simulation of ship motions coupled with tank sloshing in time domain based on OpenFOAM[J].Journal of Ship Mechanics,2012,16(7):750-758.

    [36] 姜勝超.波浪作用下船舶運動與液艙內(nèi)流體晃蕩的耦合數(shù)值分析[D].大連:大連理工大學(xué),2013.JIANG Shengchao.Numerical simulation on coupled effect between ship motion and liquid sloshing under wave action[D].Dalian:Dalian University of Technology,2013.

    [37] 劉玉龍.波浪中浮箱運動與三維液艙晃蕩的耦合研究[D].大連:大連理工大學(xué),2014.LIU Yulong.Study on coupling effects of floating box motion and 3D sloshing in waves[D].Dalian:Dalian University of Technology,2014.

    [38] 黃碩,段文洋,游亞戈,等.液艙晃蕩與船體非線性時域耦合運動計算[J].哈爾濱工程大學(xué)學(xué)報,2014,35(9):1045-1052.HUANG Shuo,DUAN Wenyang,YOU Yage,et al.Nonlinear time domain simulation of sloshing and coupled ship motion[J].Journal of Harbin Engineering University,2014,35(9):1045-1052.

    [39] MILOVAN P,ZORN T.Simulation of sloshing loads on moving tanks[C]∥Proceedings of the 24th International Conference on Offshore Mechanics and Arctic Engineering.Halkidiki,Greece,2005.

    (編輯:呂歡歡)

    Status and outlook of research on the coupling motion between liquid sloshing and the hull body of large FLNG

    HAN Xuliang1,2XIE Bin1ZHU Xiaosong1WANG Junrong1

    (1.CNOOCResearchInstitute,Beijing100028,China; 2.ChinaUniversityofPetroleum,Beijing102249,China)

    Upon the introduction to FLNG application background, it is pointed out that the coupling motion between liquid sloshing and the hull body is one of the most serious problems for FLNG.The research status of the coupling motion between liquid sloshing and the hull body is summarized and analyzed, mainly concerning the selection of containment systems, model tests and numerical simulation technology.The characteristics of the experimental technology and the application of the numerical simulation technology are discussed.Finally, some suggestions are put forward for the directions of future research on the coupling motion.

    FLNG; liquid sloshing; coupling motion; containment system; model tests; numerical simulation

    *國家自然科學(xué)基金“南海惡劣海況下FLNG船體運動與液艙晃蕩耦合分析方法研究(編號:51609267)”、“十二五”國家科技重大專項“大型FLNG/FLPG、FDPSO關(guān)鍵技術(shù)(編號:2011ZX05026-006)”、中海油研究總院自立課題“深水油氣田開發(fā)工程技術(shù)協(xié)同創(chuàng)新(編號:2016-KJZC-013)”部分研究成果。

    韓旭亮,男,博士后,2015年畢業(yè)于哈爾濱工程大學(xué)并獲博士學(xué)位,現(xiàn)主要從事深水浮式結(jié)構(gòu)物總體性能研究工作。地址:北京市朝陽區(qū)太陽宮南街6號院 (郵編:100028)。E-mail:hanxl9@cnooc.com.cn。

    1673-1506(2017)01-0116-08

    10.11935/j.issn.1673-1506.2017.01.018

    韓旭亮,謝彬,朱小松,等.大型FLNG液艙晃蕩和船體耦合運動研究現(xiàn)狀及展望[J].中國海上油氣,2017,29(1):116-123.

    HAN Xuliang,XIE Bin,ZHU Xiaosong,et al.Status and outlook of research on the coupling motion between liquid sloshing and the hull body of large FLNG[J].China Offshore Oil and Gas,2017,29(1):116-123.

    TE835

    A

    2016-05-18 改回日期:2016-07-20

    猜你喜歡
    液艙模型試驗船體
    B型LNG液艙支座縱骨趾端處表面裂紋擴展計算
    船體行駛過程中的壓力監(jiān)測方法
    基于CFD的大型船舶液艙晃蕩研究
    反推力裝置模型試驗臺的研制及驗證
    考慮晃蕩效應(yīng)的獨立B型LNG液艙結(jié)構(gòu)多目標(biāo)優(yōu)化
    海洋工程(2016年2期)2016-10-12 05:08:07
    臺階式短加筋土擋墻行為特征的離心模型試驗
    FPSO與運輸船旁靠時液艙晃蕩與船舶運動耦合效應(yīng)分析
    巨厚堅硬巖漿巖不同配比的模型試驗研究
    焊接殘余應(yīng)力對船體結(jié)構(gòu)疲勞強度的影響分析
    焊接(2015年9期)2015-07-18 11:03:51
    電滲—堆載聯(lián)合氣壓劈烈的室內(nèi)模型試驗
    一级毛片精品| 亚洲五月婷婷丁香| 欧美少妇被猛烈插入视频| 啦啦啦视频在线资源免费观看| 国产精品1区2区在线观看. | 国产成人啪精品午夜网站| 黄色毛片三级朝国网站| 亚洲自偷自拍图片 自拍| 亚洲天堂av无毛| 午夜免费观看性视频| 免费少妇av软件| 日韩欧美国产一区二区入口| 欧美日韩精品网址| 丝袜脚勾引网站| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三区在线| 丝瓜视频免费看黄片| 日韩欧美国产一区二区入口| 国产成人系列免费观看| 久久精品熟女亚洲av麻豆精品| 日韩精品免费视频一区二区三区| 亚洲av美国av| 欧美黑人欧美精品刺激| 一个人免费在线观看的高清视频 | 国产成人精品久久二区二区91| 男人操女人黄网站| 日韩有码中文字幕| 午夜激情久久久久久久| 亚洲精品乱久久久久久| 久久综合国产亚洲精品| 在线观看免费高清a一片| av一本久久久久| www.熟女人妻精品国产| 男女高潮啪啪啪动态图| 精品国产一区二区久久| 国产av又大| 狠狠狠狠99中文字幕| 久久精品国产亚洲av香蕉五月 | 菩萨蛮人人尽说江南好唐韦庄| 丝袜人妻中文字幕| 亚洲人成77777在线视频| 人人妻人人爽人人添夜夜欢视频| 久久久精品免费免费高清| 久久精品aⅴ一区二区三区四区| 宅男免费午夜| 国产在线视频一区二区| 日本91视频免费播放| 精品国产超薄肉色丝袜足j| 亚洲精品国产av成人精品| 国产高清国产精品国产三级| 亚洲综合色网址| bbb黄色大片| tocl精华| 亚洲精品乱久久久久久| 亚洲情色 制服丝袜| 99九九在线精品视频| 日韩中文字幕欧美一区二区| 黄色 视频免费看| 蜜桃在线观看..| 久久精品成人免费网站| 高潮久久久久久久久久久不卡| 激情视频va一区二区三区| 91麻豆精品激情在线观看国产 | 啦啦啦啦在线视频资源| 亚洲av美国av| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产综合久久久| 99久久人妻综合| 亚洲av电影在线进入| 国产精品av久久久久免费| 久久国产精品大桥未久av| 国产一区二区三区综合在线观看| 两个人免费观看高清视频| 久久精品国产a三级三级三级| 国产亚洲一区二区精品| 一本久久精品| 亚洲精品av麻豆狂野| 王馨瑶露胸无遮挡在线观看| 亚洲自偷自拍图片 自拍| 欧美日韩黄片免| 国产又爽黄色视频| 91精品三级在线观看| 欧美在线一区亚洲| 欧美另类亚洲清纯唯美| 亚洲专区国产一区二区| 国产精品九九99| 亚洲国产av新网站| 丝袜脚勾引网站| 一区二区日韩欧美中文字幕| 19禁男女啪啪无遮挡网站| 国产男女超爽视频在线观看| 秋霞在线观看毛片| 女性被躁到高潮视频| 91精品国产国语对白视频| 精品卡一卡二卡四卡免费| 色播在线永久视频| 我要看黄色一级片免费的| 下体分泌物呈黄色| 久久久久久人人人人人| 丁香六月天网| 一区二区日韩欧美中文字幕| 欧美午夜高清在线| 国产免费现黄频在线看| 国产一区二区激情短视频 | 爱豆传媒免费全集在线观看| 黄片播放在线免费| 婷婷丁香在线五月| 精品国产一区二区久久| 久久99一区二区三区| 亚洲欧美激情在线| 成人黄色视频免费在线看| 精品国产一区二区三区四区第35| 9色porny在线观看| 99久久精品国产亚洲精品| 亚洲九九香蕉| 精品亚洲成a人片在线观看| 狂野欧美激情性bbbbbb| 亚洲国产精品一区三区| www.精华液| 波多野结衣一区麻豆| 青草久久国产| 十八禁人妻一区二区| 亚洲成人国产一区在线观看| 亚洲精品一区蜜桃| 岛国在线观看网站| 一本综合久久免费| 脱女人内裤的视频| 久久国产精品人妻蜜桃| 国产精品免费大片| 麻豆国产av国片精品| 一级,二级,三级黄色视频| a在线观看视频网站| 久久久国产精品麻豆| a级片在线免费高清观看视频| 久热这里只有精品99| 欧美精品高潮呻吟av久久| 精品人妻一区二区三区麻豆| 天天躁夜夜躁狠狠躁躁| 老汉色∧v一级毛片| 国产精品久久久久久精品电影小说| 大片免费播放器 马上看| 成人18禁高潮啪啪吃奶动态图| 亚洲人成电影免费在线| 免费在线观看黄色视频的| 午夜精品久久久久久毛片777| 久久精品成人免费网站| a级毛片在线看网站| 亚洲精品美女久久久久99蜜臀| 精品人妻在线不人妻| 久久精品熟女亚洲av麻豆精品| 欧美97在线视频| 在线观看一区二区三区激情| 亚洲全国av大片| 丝袜脚勾引网站| 人妻人人澡人人爽人人| 午夜激情久久久久久久| 嫁个100分男人电影在线观看| 男女之事视频高清在线观看| 一二三四社区在线视频社区8| 日韩电影二区| 丁香六月天网| 国产欧美日韩一区二区三区在线| 国产视频一区二区在线看| 9热在线视频观看99| 国产成+人综合+亚洲专区| 天堂8中文在线网| 精品一区在线观看国产| 18禁裸乳无遮挡动漫免费视频| 99热国产这里只有精品6| 国产一区二区三区av在线| a在线观看视频网站| 91av网站免费观看| 热99久久久久精品小说推荐| 波多野结衣av一区二区av| 韩国精品一区二区三区| 波多野结衣av一区二区av| avwww免费| 国产精品久久久av美女十八| 亚洲少妇的诱惑av| 日本五十路高清| 欧美黑人精品巨大| 久热这里只有精品99| a级片在线免费高清观看视频| 国产精品久久久av美女十八| 考比视频在线观看| 日韩欧美一区视频在线观看| 精品人妻熟女毛片av久久网站| 啦啦啦 在线观看视频| 男女之事视频高清在线观看| 黄片大片在线免费观看| 欧美日韩成人在线一区二区| 亚洲综合色网址| 亚洲 欧美一区二区三区| 国产亚洲精品一区二区www | 精品熟女少妇八av免费久了| 国产日韩欧美视频二区| 精品亚洲成国产av| 美女午夜性视频免费| 国产免费av片在线观看野外av| 久久久久精品人妻al黑| 亚洲黑人精品在线| 超色免费av| 亚洲精品自拍成人| 69精品国产乱码久久久| 日本av免费视频播放| 在线av久久热| www.熟女人妻精品国产| 久久ye,这里只有精品| 美国免费a级毛片| 丝袜在线中文字幕| 久久精品熟女亚洲av麻豆精品| 久久亚洲国产成人精品v| 成年美女黄网站色视频大全免费| 一级片免费观看大全| 99久久精品国产亚洲精品| 精品第一国产精品| 亚洲专区字幕在线| 亚洲欧美精品综合一区二区三区| 亚洲专区中文字幕在线| 精品视频人人做人人爽| 国产精品九九99| 久久天躁狠狠躁夜夜2o2o| 亚洲免费av在线视频| av天堂久久9| 男男h啪啪无遮挡| 国产精品久久久久久精品电影小说| 狂野欧美激情性bbbbbb| 99国产极品粉嫩在线观看| 亚洲国产毛片av蜜桃av| 国产免费av片在线观看野外av| 国产精品一区二区在线观看99| 99香蕉大伊视频| 菩萨蛮人人尽说江南好唐韦庄| 久久综合国产亚洲精品| 一边摸一边抽搐一进一出视频| 飞空精品影院首页| 亚洲成人免费av在线播放| 黄频高清免费视频| 日韩大片免费观看网站| 欧美人与性动交α欧美软件| 女人爽到高潮嗷嗷叫在线视频| 亚洲全国av大片| 国产欧美日韩精品亚洲av| 亚洲欧美精品自产自拍| 成年人黄色毛片网站| 欧美精品av麻豆av| 九色亚洲精品在线播放| 午夜精品国产一区二区电影| 久久av网站| 嫩草影视91久久| 午夜福利视频在线观看免费| 黄色视频不卡| 后天国语完整版免费观看| 精品第一国产精品| 亚洲欧美日韩另类电影网站| 国产又爽黄色视频| 成人免费观看视频高清| 亚洲av国产av综合av卡| 国产一区二区在线观看av| 9色porny在线观看| 久久人人爽人人片av| 丝袜喷水一区| 精品少妇黑人巨大在线播放| 永久免费av网站大全| 色婷婷av一区二区三区视频| 国产成人免费无遮挡视频| 亚洲精品在线美女| 亚洲,欧美精品.| av天堂久久9| 亚洲伊人色综图| 午夜久久久在线观看| 亚洲精品日韩在线中文字幕| 亚洲欧美激情在线| 亚洲国产毛片av蜜桃av| 国产精品av久久久久免费| 亚洲av男天堂| 王馨瑶露胸无遮挡在线观看| 丝瓜视频免费看黄片| 免费观看人在逋| 久久久久久亚洲精品国产蜜桃av| 天天躁日日躁夜夜躁夜夜| 精品亚洲乱码少妇综合久久| 少妇裸体淫交视频免费看高清 | 精品国产超薄肉色丝袜足j| 午夜免费观看性视频| 国产精品自产拍在线观看55亚洲 | 91字幕亚洲| 宅男免费午夜| 国产一区二区三区av在线| 亚洲av电影在线观看一区二区三区| 亚洲欧美精品综合一区二区三区| 精品一区二区三卡| 女性生殖器流出的白浆| 99精品久久久久人妻精品| 人人妻人人添人人爽欧美一区卜| 欧美亚洲日本最大视频资源| 日本黄色日本黄色录像| 欧美日韩福利视频一区二区| 国产有黄有色有爽视频| 欧美日韩视频精品一区| 久久久久久久久免费视频了| 女人被躁到高潮嗷嗷叫费观| 亚洲精品粉嫩美女一区| 久久久国产成人免费| 一级毛片精品| 丝袜在线中文字幕| 人人妻人人澡人人爽人人夜夜| 国产麻豆69| 亚洲 欧美一区二区三区| 日韩制服丝袜自拍偷拍| 国产av精品麻豆| 夜夜夜夜夜久久久久| 美女国产高潮福利片在线看| 国产日韩欧美在线精品| 国产精品av久久久久免费| 欧美黄色淫秽网站| a级毛片黄视频| 久久毛片免费看一区二区三区| a在线观看视频网站| 看免费av毛片| 女人高潮潮喷娇喘18禁视频| 日本a在线网址| 国产精品免费大片| 建设人人有责人人尽责人人享有的| 国产欧美亚洲国产| 国产精品亚洲av一区麻豆| 欧美日韩av久久| 国产精品二区激情视频| 18禁黄网站禁片午夜丰满| 91九色精品人成在线观看| 在线观看一区二区三区激情| 国产福利在线免费观看视频| 精品一区在线观看国产| 一区二区三区激情视频| 日韩,欧美,国产一区二区三区| 操出白浆在线播放| 国产免费一区二区三区四区乱码| 国产精品麻豆人妻色哟哟久久| 久久久久国内视频| 人人妻人人添人人爽欧美一区卜| 好男人电影高清在线观看| 欧美激情高清一区二区三区| 大片电影免费在线观看免费| 国产精品1区2区在线观看. | 欧美成狂野欧美在线观看| 国产1区2区3区精品| av在线老鸭窝| 一区二区三区精品91| 一级毛片女人18水好多| 亚洲av成人不卡在线观看播放网 | 久久久久网色| 久久久久精品人妻al黑| 十分钟在线观看高清视频www| 777久久人妻少妇嫩草av网站| 高清视频免费观看一区二区| 热99re8久久精品国产| 国产精品熟女久久久久浪| 欧美午夜高清在线| 国产欧美日韩一区二区精品| 天天躁日日躁夜夜躁夜夜| 日韩欧美一区二区三区在线观看 | 欧美日韩国产mv在线观看视频| 91字幕亚洲| 叶爱在线成人免费视频播放| 中文精品一卡2卡3卡4更新| 日韩 亚洲 欧美在线| 日韩熟女老妇一区二区性免费视频| 两人在一起打扑克的视频| 下体分泌物呈黄色| 久久av网站| 国产一区二区在线观看av| 欧美精品亚洲一区二区| 日本精品一区二区三区蜜桃| 操出白浆在线播放| 在线十欧美十亚洲十日本专区| 亚洲久久久国产精品| 国产区一区二久久| 欧美日韩av久久| 精品久久久久久久毛片微露脸 | 搡老熟女国产l中国老女人| 人妻 亚洲 视频| 日本一区二区免费在线视频| 免费日韩欧美在线观看| 宅男免费午夜| 亚洲伊人久久精品综合| 亚洲成人免费电影在线观看| 女人久久www免费人成看片| 欧美另类亚洲清纯唯美| 国产成人欧美| 国产免费福利视频在线观看| 黄色怎么调成土黄色| 亚洲成人免费av在线播放| 又大又爽又粗| 午夜日韩欧美国产| 在线永久观看黄色视频| 日本wwww免费看| 久热爱精品视频在线9| 久久国产精品大桥未久av| 法律面前人人平等表现在哪些方面 | 亚洲五月色婷婷综合| 成在线人永久免费视频| 夫妻午夜视频| av在线app专区| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| 亚洲成av片中文字幕在线观看| av片东京热男人的天堂| 丰满饥渴人妻一区二区三| 亚洲黑人精品在线| 在线观看舔阴道视频| 亚洲成人国产一区在线观看| 爱豆传媒免费全集在线观看| 久久久精品区二区三区| 在线看a的网站| 在线观看免费高清a一片| avwww免费| 99久久国产精品久久久| 亚洲精品第二区| 久久人妻熟女aⅴ| 91精品国产国语对白视频| 日韩欧美国产一区二区入口| 桃红色精品国产亚洲av| 国产老妇伦熟女老妇高清| 国产免费现黄频在线看| 国产高清国产精品国产三级| 国产日韩欧美亚洲二区| 嫁个100分男人电影在线观看| 欧美国产精品va在线观看不卡| 欧美激情极品国产一区二区三区| 亚洲精品国产色婷婷电影| 日本91视频免费播放| 欧美日韩一级在线毛片| 97精品久久久久久久久久精品| 中文字幕另类日韩欧美亚洲嫩草| xxxhd国产人妻xxx| 久久毛片免费看一区二区三区| 爱豆传媒免费全集在线观看| 久久久久久久大尺度免费视频| av超薄肉色丝袜交足视频| 国产精品久久久久成人av| netflix在线观看网站| kizo精华| 一级片'在线观看视频| 日本91视频免费播放| 久久久国产一区二区| kizo精华| 男男h啪啪无遮挡| 狠狠婷婷综合久久久久久88av| 手机成人av网站| 嫁个100分男人电影在线观看| 12—13女人毛片做爰片一| 黄色片一级片一级黄色片| av天堂在线播放| 午夜激情久久久久久久| e午夜精品久久久久久久| 99国产极品粉嫩在线观看| 一本—道久久a久久精品蜜桃钙片| 日韩精品免费视频一区二区三区| 青青草视频在线视频观看| 视频区欧美日本亚洲| 人妻人人澡人人爽人人| 狠狠狠狠99中文字幕| 一边摸一边做爽爽视频免费| 视频区图区小说| 国产在线观看jvid| 中国国产av一级| 女人久久www免费人成看片| 国产精品久久久久成人av| 亚洲精品自拍成人| 男人舔女人的私密视频| 午夜激情久久久久久久| 亚洲久久久国产精品| netflix在线观看网站| 99国产精品一区二区蜜桃av | 老司机深夜福利视频在线观看 | 欧美+亚洲+日韩+国产| 国产97色在线日韩免费| 他把我摸到了高潮在线观看 | 久久人人97超碰香蕉20202| 可以免费在线观看a视频的电影网站| 国产精品香港三级国产av潘金莲| 成人18禁高潮啪啪吃奶动态图| 国产亚洲欧美在线一区二区| 成年人午夜在线观看视频| 亚洲精品久久久久久婷婷小说| 欧美 日韩 精品 国产| 操美女的视频在线观看| 在线永久观看黄色视频| 精品高清国产在线一区| 满18在线观看网站| 日韩中文字幕欧美一区二区| 色婷婷av一区二区三区视频| 在线永久观看黄色视频| 国产成人精品久久二区二区免费| 欧美精品人与动牲交sv欧美| 老熟妇仑乱视频hdxx| 成人国产av品久久久| 视频区欧美日本亚洲| 国产老妇伦熟女老妇高清| 欧美精品一区二区免费开放| 国产亚洲av片在线观看秒播厂| 伦理电影免费视频| 国产亚洲欧美在线一区二区| 美女国产高潮福利片在线看| 69精品国产乱码久久久| 日韩视频在线欧美| 脱女人内裤的视频| 满18在线观看网站| 成人av一区二区三区在线看 | 久久毛片免费看一区二区三区| 久久 成人 亚洲| 日本a在线网址| 999久久久精品免费观看国产| av片东京热男人的天堂| 1024香蕉在线观看| 亚洲成人手机| 伊人亚洲综合成人网| 老熟妇仑乱视频hdxx| 中文精品一卡2卡3卡4更新| 99国产精品免费福利视频| 91九色精品人成在线观看| 纵有疾风起免费观看全集完整版| 黄色视频在线播放观看不卡| 天天操日日干夜夜撸| av免费在线观看网站| 亚洲精品一二三| 久久国产精品大桥未久av| 黑人操中国人逼视频| 亚洲av电影在线进入| 热99国产精品久久久久久7| 母亲3免费完整高清在线观看| 侵犯人妻中文字幕一二三四区| 免费人妻精品一区二区三区视频| 亚洲专区国产一区二区| 1024视频免费在线观看| 九色亚洲精品在线播放| 国产精品成人在线| 国产高清视频在线播放一区 | 国产日韩欧美视频二区| 波多野结衣一区麻豆| 自线自在国产av| 国产精品亚洲av一区麻豆| 久久久精品免费免费高清| 精品第一国产精品| 十八禁网站网址无遮挡| 精品国产国语对白av| 男女之事视频高清在线观看| 欧美日韩av久久| 成人黄色视频免费在线看| av天堂在线播放| 老司机午夜福利在线观看视频 | 99久久人妻综合| 国产一区二区在线观看av| 国产亚洲av片在线观看秒播厂| 欧美国产精品一级二级三级| 亚洲久久久国产精品| 免费一级毛片在线播放高清视频 | 午夜免费观看性视频| 亚洲国产精品一区三区| 亚洲精品国产精品久久久不卡| 一本大道久久a久久精品| 伊人亚洲综合成人网| 亚洲精品在线美女| 一进一出抽搐动态| 亚洲成人免费av在线播放| 国产精品 欧美亚洲| 久久久久精品人妻al黑| 久久精品国产综合久久久| 亚洲七黄色美女视频| 国产99久久九九免费精品| 91成年电影在线观看| 午夜福利影视在线免费观看| 老鸭窝网址在线观看| 老司机靠b影院| 亚洲国产av新网站| 日韩电影二区| 人人妻人人添人人爽欧美一区卜| 久久女婷五月综合色啪小说| 下体分泌物呈黄色| 久久天躁狠狠躁夜夜2o2o| 老熟妇仑乱视频hdxx| 亚洲第一欧美日韩一区二区三区 | 国产高清视频在线播放一区 | 男女无遮挡免费网站观看| 青青草视频在线视频观看| 天天添夜夜摸| 欧美国产精品va在线观看不卡| 十八禁人妻一区二区| 国产欧美日韩一区二区三区在线| 国产精品亚洲av一区麻豆| 久久天堂一区二区三区四区| 国产成人免费无遮挡视频| 大型av网站在线播放| 国产精品国产三级国产专区5o| 别揉我奶头~嗯~啊~动态视频 | 亚洲综合色网址| 国产欧美日韩精品亚洲av| 高清欧美精品videossex| 最新的欧美精品一区二区| 亚洲专区中文字幕在线| 欧美日韩亚洲国产一区二区在线观看 | 亚洲一码二码三码区别大吗| 啦啦啦啦在线视频资源| 亚洲国产av新网站| 精品视频人人做人人爽| 黄色怎么调成土黄色| 国产有黄有色有爽视频| 777久久人妻少妇嫩草av网站| 一级片免费观看大全| 亚洲专区中文字幕在线| 亚洲综合色网址| 日本黄色日本黄色录像| 免费日韩欧美在线观看| 最近最新免费中文字幕在线| 国产高清视频在线播放一区 |