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

    水深對軟剛臂單點系泊FPSO動力響應的影響

    2016-10-12 02:33:04劉成義唐友剛
    海洋工程 2016年1期
    關(guān)鍵詞:單點系泊淺水系泊

    劉成義,唐友剛,李 焱

    (1.天津大學 建筑工程學院,天津 300072;2.天津大學 水利工程仿真與安全國家重點實驗室,天津 300072)

    水深對軟剛臂單點系泊FPSO動力響應的影響

    劉成義1,2,唐友剛1,2,李 焱1,2

    (1.天津大學 建筑工程學院,天津 300072;2.天津大學 水利工程仿真與安全國家重點實驗室,天津 300072)

    系泊系統(tǒng)的定位能力是淺水油田作業(yè)的軟剛臂式單點系泊FPSO安全作業(yè)的重要保障,為研究不同水深/吃水比下單點系泊系統(tǒng)的受力性能,針對一艘16萬噸級軟剛臂單點系泊FPSO,在線性三維勢流理論的基礎(chǔ)上,基于多體動力學方法,建立FPSO-系泊腿-軟剛臂的耦合模型,采用Newman近似法和Pinkster近似法分析了FPSO所受二階波浪力,在時域內(nèi)計算了不同水深/吃水比對系泊系統(tǒng)動力響應性能的影響。結(jié)果表明,隨著水深/吃水比的增加,Newman近似法計算得到二階波浪力先增大后減小,引起單點系泊系統(tǒng)載荷先增大后減??;而Pinkster近似法計算得到的二階波浪力逐漸減下,引起單點載荷逐漸減下。在淺水條件下,Pinkster近似法具有較好的適用性,Newman近似法嚴重低估了FPSO所受的二階波浪力;在深水條件下,Newman近似法能滿足工程計算的要求;適用兩種方法的臨界水深/吃水比為1.64。

    多體動力學;軟剛臂系泊系統(tǒng);Newman/Pinkster近似法;二階波浪力;動力響應性能;水深/吃水比

    Abstract:With the increasing operation life of FPSO and mooring systems,there have been some accidents,which have a great impact on safety in production in oil field.More attentions have been paid to the safety of the Soft Yoke Mooring Systems(SYMS).Dynamic response capability of the SYMS under different WD/T(Ratio of Water depth and FPSO Draft) was analyzed based on multi-body dynamic method.The multi-body coupling model including SYMS and FPSO was established and the hydrodynamic parameter of each body was calculated with linear 3-D potential theory.Second-order wave drift force was calculated by Newman approximation and Pinkster approximation,and the forces on the SYMS under different WD/T were analyzed afterward.A comparison of results using different approximation methods was carried out.With the Newman approximation,the second-order wave force increases first and decrease afterwards with the increase of WD/T,leading to the forces on the SYMS changing synchronously; with the Pinkster approximation,the second-order wave force continually decreases with the increase of WD/T,leading to the forces on the SYMS changing synchronously.It is concluded that Pinkster approximation has a good applicability in shallow water while Newman approximation can meet the requirements of engineering calculations in deep water; and the threshold of WD/T is 1.64.

    Keywords:multi-body method; soft yoke mooring system; Newman/Pinkster approximation; second-order wave drift force; dynamic response; ratio of water depth and draft

    對于淺水油田作業(yè)的軟剛臂式單點系泊FPSO,其系泊系統(tǒng)的定位能力是FPSO安全作業(yè)的重要保障。在淺水條件下,水深變化將會極大地影響FPSO所受的二階波浪力,進而對系泊系統(tǒng)的載荷產(chǎn)生影響。

    FPSO所受二階波浪力主要包括平均漂移力、二階慢漂波浪力及二階和頻波浪力,國內(nèi)外學者已經(jīng)對FPSO所受二階波浪力受水深的影響做了大量地研究。

    在二階波浪力的計算方法上,Guilaume等[1]對二階低頻載荷的計算方法進行總結(jié),指出目前計算二階波浪載荷傳遞函數(shù)(QTF)矩陣的主要方法有:全QTF法、忽略自由表面積分的全QTF法、Pinkster近似法以及Newman近似法四種。QTF矩陣主對角線元素為平均漂移力的傳遞函數(shù),主要計算方法為近場法、遠場法兩種;非主對角線元素為二階差頻波浪載荷傳遞函數(shù),完整的數(shù)值計算將耗費大量的時間及計算資源,通常通過Newman近似或Pinkster近似法求解。也有學者對以上計算方法進行了改進,Chen等[2]將遠場法與近場的壓力積分法相結(jié)合,提出了計算平均漂移力的中場法;Cong等[3]基于高階邊界元和自適應的單元積分法,提出了一種新的低頻慢漂力的計算方法;但這些改進算法應用并不廣泛。目前,二階波浪力的計算仍然以Newman近似法和Pinkster近似法為主。

    關(guān)于水深變化對FPSO所受二階波浪力的影響,陳憶寧等[4]在不考慮系泊系統(tǒng)的情況下,以一艘16萬噸級的FPSO為對象,分別采用Newman近似法和全QTF方法計算不同水深下FPSO所受的二階波浪力,結(jié)果表明,Newman近似法在深水條件下具有較好的適用性。Naciri等[5]采用有限元方法對淺水 LNG 船水動力特性進行分析,并與實驗結(jié)果進行對比,認為對于淺水中船形海洋結(jié)構(gòu)物所受二階波浪力,Pinkster 近似法可以滿足計算和研究的要求。李欣等[6]分別采用Pinkster近似法和實驗的方法,研究淺水條件下不同水深對FPSO所受二階波浪力的影響,結(jié)論認為二階波浪力/力矩隨水深的增大而減小。Fonseca和Pessoa等[7-8]用數(shù)值和實驗方法研究了雙色波作用下固定浮體上所受的一階波浪力、平均漂移力以及二階差頻波浪力隨水深變化的規(guī)律,認為隨水深的減小,波浪二階力明顯增大,一階力略有增大,且淺水中二階波浪力的計算較之于深水中更為重要。Xiao等[9]以軟剛臂單點系泊方式定位于淺水油田的16萬噸級FPSO為對象,通過實驗研究不同水深下的低頻縱蕩運動和系泊力,結(jié)論認為:低頻縱蕩運動和系泊力隨著水深的減小大大增加,在水深吃水比較小時,力的增加非常劇烈。

    可見,目前大部分研究將FPSO視作固定的海洋結(jié)構(gòu)物,研究水深變化對其所受二階波浪力的影響,并未考慮FPSO在實際系泊狀態(tài)下的動態(tài)效應;亦未考慮水深變化對軟剛臂單點系泊系統(tǒng)的動力響應性能的影響。為研究水深變化對軟剛臂式單點系泊FPSO系泊系統(tǒng)動力響應性能的影響,本文基于多體動力學的方法,分別采用Newman近似法和Pinkster近似法考慮二階波浪載荷,在時域內(nèi)計算了不同水深/吃水比(WD/T)下單點系泊系統(tǒng)載荷,得到系泊系統(tǒng)載荷隨WD/T的變化規(guī)律,并對比兩種方法的計算結(jié)果,分析產(chǎn)生這種差別的原因;考慮到Newman近似法在深水條件下具有很好的適用性,而Pinkster近似法對淺水條件下的二階波浪載荷具有更高的模擬精度,并進一步給出了適合用Pinkster近似法計算和Newman近似法計算的臨界WD/T值,為工程實際提供參考。

    1 多體耦合計算模型

    1.1FPSO及系泊系統(tǒng)主尺度參數(shù)

    FPSO所在海域水深為17.83 m,其主要參數(shù)如表1所示。軟剛臂系泊系統(tǒng)(soft yoke mooring system,SYMS)由軟剛臂及兩側(cè)系泊腿組成,主尺度參數(shù)以及軟剛臂結(jié)構(gòu)坐標系、系泊腿結(jié)構(gòu)坐標系如圖1所示,其中軟剛臂與壓載艙總重量為1 850 t,單個系泊腿重量為75 t。

    1.2全耦合有限元模型的建立

    對于軟剛臂單點系泊系統(tǒng),傳統(tǒng)的計算方法通常將單點簡化為約束縱蕩運動的非線性剛度彈簧,研究FPSO的運動響應或者單點的受力,這種方法很難模擬系泊機構(gòu)的實際運動特性。對此,許多學者基于多體動力學方法,對軟剛臂式單點系泊FPSO進行了實驗和數(shù)值研究。Wang[10]等以淺水軟剛臂系泊的FPSO為研究對象,采用多體方法模擬水下軟剛臂系泊實際結(jié)構(gòu)和FPSO的耦合作用,對FPSO的運動和動力響應進行研究;Guo等[11]通過時域內(nèi)的數(shù)值模擬和模型試驗,研究了非線性海洋環(huán)境下淺水FPSO運動響應特性,并進行甲板上浪分析;Phadke等[12]基于多體方法,對軟剛臂單點系泊的“蓬勃號”FPSO系泊系統(tǒng)的安裝進行數(shù)值分析,得到較為合理的結(jié)果。這些研究很好地證明了多體方法的適用性,文中對FPSO系泊系統(tǒng)的建模,也是基于多體方法進行的。

    采用ANSYS建立全船及系泊剛架的有限元模型,導入AQWA中進行水動力及時域響應計算??紤]到船艏尾的形狀較為復雜,定義船艏尾的網(wǎng)格密度為2 m,船中的網(wǎng)格密度為3 m。

    FPSO及系泊剛架的有限元模型如圖2所示,模型的單元總數(shù)為8 063。多體耦合模型中包含四部分結(jié)構(gòu):FPSO及船艏系泊剛架通過剛性連接組成第一部分結(jié)構(gòu);兩側(cè)系泊腿為兩個結(jié)構(gòu),上端與船艏系泊剛架通過雙軸鉸連接,釋放橫搖、縱搖兩個自由度的約束;下端與軟剛臂通過三軸鉸連接,釋放3個轉(zhuǎn)動自由度的約束;包含壓載艙的軟剛臂為第四部分結(jié)構(gòu),軟剛臂與塔架系泊點處釋放3個轉(zhuǎn)動自由度的約束。

    表1 FPSO主要參數(shù)Tab.1 Main particulars of FPSO vessel

    圖1 軟剛臂及系泊腿的主要尺度Fig.1 Main dimensions of the yoke and pendent

    圖2 軟剛臂單點系泊FPSO全耦合有限元模型Fig.2 Completely-coupled model of the FPSO and SYMS

    2 計算理論

    2.1多體耦合時域運動方程

    綜合考慮波浪、風、流載荷以及單點系泊系統(tǒng)的回復力的聯(lián)合作用,軟鋼臂式單點系泊FPSO在時域內(nèi)的運動方程為[13]:

    其中,M為FPSO質(zhì)量矩陣,通過FPSO質(zhì)量以及各回轉(zhuǎn)半徑求得;A()為最大計算頻率對應的附連水質(zhì)量矩陣,r(t)為輻射阻尼的脈沖響應函數(shù)矩陣,附連水質(zhì)量、輻射阻尼系數(shù)均采用三維勢流理論應用源匯分布法計算得到;D為FPSO慢漂阻尼矩陣,對于文中軟剛臂式單點系泊FPSO,縱蕩、橫蕩和艏搖自由度低頻漂移阻尼參照BV規(guī)范中的公式進行計算;K為FPSO靜水回復力剛度矩陣,根據(jù)FPSO水動力性能求得;F(1,2)(t)為一階、二階波浪載荷;Fw(t)為風載荷,F(xiàn)c(t)為流載荷,風、流載荷參照OCIMF規(guī)范進行計算;FSPM(t)為單點系泊系統(tǒng)提供的回復力。

    A字型剛臂結(jié)構(gòu)是在系泊單點、兩側(cè)系泊腿拉力下保持平衡,其受力形式較為簡單,此處不列出結(jié)構(gòu)動力方程。

    兩側(cè)系泊腿作為FPSO與A字型剛臂的連接桿,其時域運動方程:

    式中:mp為系泊腿質(zhì)量矩陣,Kp為考慮結(jié)構(gòu)之間相互影響的12×12的全耦合剛度矩陣。

    2.2波浪載荷計算

    選用JONSWAP譜模擬不規(guī)則海浪,根據(jù)譜峰因子取為1.8?;谶x定的波浪譜,通過傅里葉變換可以得到時域內(nèi)的隨機波面升高,進而可通過卷積積分的方式生成一階、二階波浪載荷的時域歷程。在不規(guī)則波浪作用下,作用在結(jié)構(gòu)物上的瞬時波浪力可以寫為:

    基于Newman近似法可得:

    基于Pinkster近似法可得:

    2.3動力響應計算結(jié)果

    動力響應計算方法:在多體系泊分析模塊AQWA-DRIFT中,基于多剛體動力學方法,根據(jù)FPSO、軟剛臂以及系泊腿的時域運動方程,分別建立兩側(cè)系泊腿及軟剛臂的振動方程,在時域內(nèi)通過連接點處的載荷傳遞實現(xiàn)耦合,建立多體的時域耦合方程,進行迭代計算。由于計算中計入了系統(tǒng)內(nèi)包括慢漂阻尼、非線性風流載荷等在內(nèi)的非線性項,采用數(shù)值解法求解非線性多體耦合振動方程組,因此,該方法可模擬實際環(huán)境工況下FPSO的運動響應及單點系統(tǒng)的受力。

    系泊系統(tǒng)的動力響應計算所選取兩種典型的海洋環(huán)境工況如表2所示,分別考慮了風、浪、流同向與不同向兩種工況,計算水深為17.8 m。

    根據(jù)表2給出的工況,采用多體動力學方法進行FPSO及單點系泊系統(tǒng)的耦合動力分析,模擬軟剛臂單點系泊機構(gòu)的具體連接形式,考慮軟剛臂及系泊腿加速度慣性力項對FPSO運動響應及單點系統(tǒng)載荷的影響,分別采用Newman近似法和Pinkster近似法模擬FPSO所受的二階波浪力,計算塔架系泊點處軟剛臂結(jié)構(gòu)坐標系下三個方向的載荷Fx、Fy、Fz以及系泊腿的軸向載荷Pt,在時域內(nèi)統(tǒng)計各載荷的最大值,結(jié)果如表3所示。

    根據(jù)動力響應計算結(jié)果看出,在工況1中,風、浪、流同向,由于系泊系統(tǒng)具有風標效應,單點y方向載荷較小,接近于0,可以忽略不計, Pinkster近似法的計算結(jié)果略大于Newman近似法的計算結(jié)果;原因在于FPSO迎浪條件下受力面積較小,低頻二階力作用不明顯。在工況2中,風、浪、流不同向,F(xiàn)PSO與波浪有夾角,單點各方向載荷皆不為0,且Pinkster近似法的計算結(jié)果遠大于Newman近似法的計算結(jié)果;原因在于,F(xiàn)PSO此時與波浪方向有一定的夾角,二階力作用范圍大,淺水波浪力的二階非線性效應明顯,Newman近似法嚴重低估了波浪二階力的作用,造成其計算結(jié)果遠小于Pinkster近似法的計算結(jié)果。

    表2 海洋環(huán)境條件Tab.2 Environment conditions

    表3 塔架系泊點載荷及系泊腿軸力最大值Tab.3 Maximum dynamic response of the SYMS

    注:單點載荷(Fx、Fy、Fz)皆為軟剛臂結(jié)構(gòu)坐標系下載荷結(jié)果。

    3 不同WD/T動力響應分析

    在水深的敏感性分析中,F(xiàn)PSO吃水(T)為14.5 m,以該海域?qū)嶋H水深(WD)17.8 m為基準,每次增加3 m,選取五個水深作為淺水海況,并額外選取了50 m、75 m、100 m三個水深作為深水海況,所選取的水深為17.8、20.8、23.8、26.8、29.8、50、75、100 m,對應的WD/T依次為1.23、1.43、1.64、1.85、 2.06、3.45、5.17、6.90。分別采用Newman近似法和Pinkster近似法計算軟剛臂單點系泊系統(tǒng)載荷,分析系泊系統(tǒng)載荷隨WD/T變化的規(guī)律,并對比兩種方法的計算結(jié)果。

    此外,由于Newman近似法與Pinkster近似法的根本差別在于二階波浪力的計算方法不同,為分析單點系泊系統(tǒng)載荷隨WD/T的變化規(guī)律產(chǎn)生的原因,文中對兩種近似方法下FPSO所受二階波浪力隨WD/T的變化規(guī)律進行了對比分析。

    3.1工況1系泊系統(tǒng)受力結(jié)果

    對于風、流、浪同向的工況1,軟剛臂單點系泊系統(tǒng)具有風標效應,F(xiàn)PSO僅受到縱蕩方向的載荷,單點系泊系統(tǒng)受力主要包括單點載荷Fx、Fz和系泊腿軸向載荷Pt。根據(jù)2.3節(jié)動力響應分析方法,進行不同WD/T條件下時域計算,統(tǒng)計單點系泊系統(tǒng)載荷最大值,計算結(jié)果如圖3所示。

    根據(jù)計算結(jié)果,在風、浪、流同向條件下,可得出如下結(jié)論:

    1)采用Newman近似法計算時,隨WD/T的增加,單點x方向的載荷Fx先減小后增大,到深水區(qū)域后增大趨勢平緩;單點z方向的載荷Fz和系泊腿軸力Pt均逐漸增加,但變化不大。

    2)采用Pinkster近似法計算時,隨著WD/T的增加,單點x、z方向載荷Fx、Fz和系泊腿軸力Pt呈現(xiàn)相同的規(guī)律,在淺水區(qū)域逐漸減小,深水區(qū)域趨于平緩;且Pt變化明顯,F(xiàn)x、Fz變化不大。

    3)對比兩種方法下系泊系統(tǒng)載荷的計算結(jié)果,可以看出:在淺水區(qū)域,采用Pinkster近似法計算得到的單點系泊系統(tǒng)載荷大于Newman近似法所得到的結(jié)果;深水區(qū)域,Pinkster近似法得到的結(jié)果反而小于Newman近似法;其臨界水深吃水比為1.64,因此,認為在WD/T≥1.64時,Newman近似法可滿足工程計算的要求。

    圖3 工況1系泊系統(tǒng)載荷隨WD/T變化曲線Fig.3 Forces of SYMS vs.WD/T in Condition 1

    圖4 工況1中FPSO縱蕩漂移力隨WD/T變化曲線Fig.4 Wave drift forces of FPSO in surge direction vs.WD/T in Condition 1

    3.2工況1漂移力計算結(jié)果分析

    對于風、浪、流同向的工況1,圖4給出兩種方法下FPSO所受二階波浪漂移力隨WD/T增加而變化的規(guī)律。

    隨著WD/T增加,Newman近似法計算得到的FPSO漂移力先增大后減小,到深水區(qū)域又有微弱的增大趨勢,與單點x方向的載荷呈現(xiàn)相互對應的規(guī)律;Pinkster近似法計算得到的FPSO漂移力先減小后逐漸趨于平穩(wěn),與單點系泊系統(tǒng)載荷的變化趨勢呈現(xiàn)較好的對應規(guī)律;可見,單點系泊系統(tǒng)載荷Fx、Fz、Pt隨WD/T增加而變化的情況直接受到FPSO所受的漂移力的影響。

    3.3工況2系泊系統(tǒng)受力結(jié)果

    對于風、浪、流不同向的工況2,在軟剛臂結(jié)構(gòu)坐標系下,單點系泊系統(tǒng)載荷包括Fx、Fy、Fz以及系泊腿軸向載荷Pt,統(tǒng)計不同 下各系泊系統(tǒng)載荷最大值,結(jié)果如圖5所示。由結(jié)果可知:

    1)采用Newman近似法計算時,隨著WD/T增加,單點系泊系統(tǒng)各載荷括Fx、Fy、Fz都先增大后減小,在深水區(qū)域趨于平緩,在WD/T=2.06(水深29.8 m)處出現(xiàn)峰值。

    2)采用Pinkster近似法計算,WD/T介于1.23~1.64(水深17.8~23.8 m)之間時,隨WD/T增加,系泊系統(tǒng)各項載荷逐漸減?。籛D/T介于1.64~2.05(水深23.8~29.8 m)之間時,隨WD/T增加,系泊系統(tǒng)各項載荷逐漸增大;在WD/T大于2.06后(水深大于29.8 m),隨WD/T增加,系泊系統(tǒng)各項載荷逐漸減小,在深水區(qū)域趨于平緩。

    圖5 工況2系泊系統(tǒng)載荷最大值隨WD/T變化曲線Fig.5 Forces of SYSM vs.WD/T in Condition 2

    3)對比Pinkster近似法和Newman近似法下計算得到的單點載荷結(jié)果,可以看出:①在各個WD/T條件下,Pinkster近似法的計算結(jié)果均大于Newman近似法計算結(jié)果;②在極淺水區(qū)域,WD/T介于1.23~1.64之間時,隨著WD/T的增加,Pinkster近似法計算時系泊系統(tǒng)各項載荷逐漸減小,而Newman近似法計算得到的各項載荷逐漸增加,二者呈現(xiàn)出相反的變化規(guī)律,且Pinkster近似法計算結(jié)果遠大于Newman近似法的計算結(jié)果;③ 超過1.64(水深23.8 m)后,隨著WD/T的增加,兩種方法計算得到的系泊系統(tǒng)載荷都呈現(xiàn)先增大后減小變化規(guī)律,到深水區(qū)域后趨于平緩,都在WD/T(水深29.8 m)處出現(xiàn)峰值; Pinkster近似法計算得到的載荷略大于Newman近似法計算結(jié)果。

    4)在WD/T≥1.64時,Newman近似法和Pinkster近似法計算的結(jié)果隨水深的變化規(guī)律一致,且二者之間差別不大,可以認為Newman近似法滿足工程實際的計算要求。

    3.4工況2漂移力計算結(jié)果分析

    對于風、浪、流不同向的工況2,圖6給出不同方法下WD/T改變時FPSO在橫蕩、縱蕩自由度所受漂移力,圖7給出FPSO所受艏搖自由度漂移力矩。根據(jù)漂移力計算結(jié)果可知:

    1)采用Newman近似法計算時,隨WD/T的增加,F(xiàn)PSO所受的二階波浪力變化規(guī)律為:在縱蕩方向上,漂移力變化較小,從淺水區(qū)域到深水區(qū)域略有增大;橫蕩方向上,漂移力先減小后逐漸趨于平緩;艏搖方向上,漂移力先增大后減小,在深水區(qū)趨于平緩。對比系泊系統(tǒng)載荷與FPSO所受二階波浪力隨WD/T的變化規(guī)律可知,在Newman近似法條件下,系泊系統(tǒng)載荷與艏搖方向漂移力的變化規(guī)律一致。

    2)采用Pinkster近似法計算FPSO二階差頻力時,隨WD/T的增加,F(xiàn)PSO所受的漂移力變化規(guī)律為:在縱蕩方向上,漂移力逐漸減小后趨于平緩,整體變化較?。粰M蕩方向上,當WD/T≤2.06時,漂移力先減小后增大,隨后漂移力逐漸減小并趨于平緩;艏搖方向上,漂移力先減小后趨于平緩,在淺水區(qū)域內(nèi),艏搖漂移力矩接近一階艏搖波浪力矩。可見,在Pinkster近似法條件下,單點系泊系統(tǒng)載荷隨WD/T增加而變化的規(guī)律與FPSO橫蕩、縱蕩、艏搖方向度所受漂移力變化的規(guī)律都有很好的一致性。

    3)對比兩種近似法下漂移力計算結(jié)果,發(fā)現(xiàn)FPSO橫蕩方向漂移力差別不大,在淺水區(qū)域,Pinkster近似法載荷計算結(jié)果較大;在深水區(qū)域,Newman近似法的載荷計算結(jié)果較大。

    綜上所述,在風、浪、流不同向條件下,軟剛臂單點系泊系統(tǒng)載荷的變化規(guī)律主要由FPSO艏搖方向所受的漂移力矩決定。

    圖6 工況2FPSO漂移力隨WD/T變化曲線Fig.6 Wave drift forces of FPSO vs.WD/T in Condition 2

    圖7 工況2FPSO艏搖漂移力隨WD/T變化曲線Fig.7 Wave drift forces of FPSO in yaw direction vs.WD/T in Condition 2

    4 結(jié) 語

    針對一艘淺水田油田作業(yè)、吃水為14.5 m的16萬噸級軟剛臂單點系泊FPSO,基于多體動力學方法,分別采用Newman近似法和Pinkster近似法計算了不同WD/T下FPSO所受的二階波浪力和軟剛臂單點系泊系統(tǒng)載荷,研究各載荷隨WD/T的變化規(guī)律,得出以下結(jié)論:

    1)采用Newman近似法計算時,隨著WD/T的增加,F(xiàn)PSO所受的漂移力先增大后減小,引起單點系泊系統(tǒng)載荷的先增大后減小。

    2)采用Pinkster近似法計算時,隨著WD/T的增加,F(xiàn)PSO所受的漂移力逐漸減小,引起單點系泊系統(tǒng)載荷的逐漸減小。

    3)在淺水區(qū)域,風、浪、流同向時,單點系泊系統(tǒng)載荷隨WD/T的變化規(guī)律由FPSO所受到縱蕩方向的二階波浪力決定;風、浪、流不同向時,系泊系統(tǒng)載荷隨WD/T變化的規(guī)律由艏搖方向的二階波浪力矩決定。

    4)淺水條件下,二階波浪力非線性效應明顯,采用Newman近似法計算僅能獲得定常的低頻漂移力成分,無法獲得其高頻與低頻成分,因此該方法不能滿足實際工程要求。

    5)對于吃水為14.5 m的軟鋼臂單點系泊FPSO,在WD/T≤1.64時,淺水效應明顯,采用Pinkster近似法較為準確的估算FPSO所受的二階力更符合工程實際狀況;在WD/T≥1.64時,Newman近似法可滿足工程計算的要求;這兩種方法計算二階波浪力的臨界WD/T為1.64。

    [1] DE HAUTECLOCQUE G,REZENDE F,WAALS O,et al.Review of approximations to evaluate second-order low-frequency load[C]//Proceedings of the 31st International Conference on Ocean,Offshore and Arctic Engineering.2012:363-371.

    [2] CHEN X B,REZENDE F.Efficient computations of second-order low-frequency wave load[C]//Proceedings of the 28th International Conference on Ocean,Offshore and Arctic Engineering.2009:525-532.

    [3] CONG P W,GOU Y,TENG B.A new approach to low-frequency QTF and its application in predicting slow drift force[J].Ocean Engineering,2012,53:25-37.

    [4] 陳憶寧,陳剛,肖龍飛.淺水 FPSO 二階力與慢漂阻尼計算分析[C]//第十六屆中國海洋 (岸) 工程學術(shù)討論會論文集 (上冊).北京:海洋出版社,2013.(CHEN Y N,CHEN G,XIAO L F.Research of second order wave drift force and damping in shallow water[C]//The 16th China Ocean (Shore) Engineering Symposium Proceedings.Beijing:Ocean Press,2013.(in Chinese))

    [5] NACIRI M,SERGENT E.Diffraction/Radiation of 135 000 m3storage capacity LNG carrier in shallow water:A benchmark study[C]//Proceedings of the 28th International Conference on Ocean,Offshore and Arctic Engineering.2009:637-647.

    [6] 李欣,王磊,楊建民.淺水浮式生產(chǎn)儲油系統(tǒng)二階波浪慢漂力數(shù)值計算[J].上海交通大學學報,2006,40(6):997-999.(LI X,WANG L,YANG J M.The numerical calculation on the second order wave drift force of FPSO in shallow water[J].Journal of Shanghai Jiao Tong University,2006,40(6):997-999.(in Chinese))

    [7] FONSECA N,PESSOA J,MAVRAKOS S,et al.Experimental and numerical investigation of the slowly varying wave exciting drift forces on a restrained body in bi-chromatic waves[J].Ocean Engineering,2011,38(17):2000-2014.

    [8] PESSOA J,FONSECA N.Investigation of depth effects on the wave exciting low frequency drift forces by different approximation methods[J].Applied Ocean Research,2013,42:182-199.

    [9] XIAO Longfei ,YANG Jianmin,LI Xin.Shallow water effects on surge motion and load of soft yoke moored FPSO[J].Chinese Ocean Engineering,2007,21(2):187-196.

    [10] WANG S Q,LI S Y,CHEN X H.Dynamical analysis of a soft yoke moored FPSO in shallow waters[C]//Proceeding of the 6th International Conference on Asian and Pacific Coasts (APAC).2012.

    [11] GUO B,XIAO L F,YANG J M.Analysis on motions and green water of FPSOs in shallow water with non-collinear environments[C]//Proceedings of the 29th International Conference on Ocean,Offshore and Arctic Engineering.2010:419-427.

    [12] PHADKE A C,LIAO S,CHEN D,et al.Assessment of SYMS survival limit for “Peng Bo” FPSO during installation[C]//Proceedings of the 29th International Conference on Ocean,Offshore and Arctic Engineering.2010:541-548.

    [13] Institute N M T R.SIMO-theory manual[G].Norway:2004.

    [14] TENG B,CONG P,GOU Y.Examination on low-frequency QTF of a platform[J].Journal of Hydrodynamics,Ser.B,2010,22(5,Supplement 1):969-974.

    [15] Oil Companies International Marine Forum.Prediction of wind and current loads on VLCCs[M].London:Witherby & Co.LTD.,1994.

    Investigation of depth effects on the dynamic response performance of soft yoke mooring FPSO

    LIU Chengyi1,2,TANG Yougang1,2,LI Yan1,2

    (1.Tianjin University,School of Civil Engineering,Tianjin 300072,China;2.State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin 300072,China)

    P751

    A

    10.16483/j.issn.1005-9865.2016.01.004

    1005-9865(2016)01-0025-08

    2014-09-15

    國家自然科學基金(51279130);國家自然科學基金重點項目(51239008)

    劉成義(1990-),男,湖北隨州人,碩士,主要研究FPSO系泊系統(tǒng)動力響應。E-mail:liuchengyii@163.com

    唐友剛(1952-),男,教授,博士生導師。 E-mail:tangyougang_td@163.com

    猜你喜歡
    單點系泊淺水系泊
    基于CFD方法的系泊狀態(tài)FSO舵載荷計算
    高韌性海洋用三級系泊鏈鋼的研制
    山東冶金(2022年2期)2022-08-08 01:50:50
    新型淺水浮托導管架的應用介紹
    云南化工(2021年10期)2021-12-21 07:33:40
    什么是單點系泊
    石油知識(2019年3期)2019-02-13 03:19:08
    帶阻尼的隨機淺水波方程的隨機吸引子
    基于懸鏈線方程的系泊系統(tǒng)分析
    30萬噸單點的海上責任
    珠江水運(2016年3期)2016-03-23 14:34:47
    (2+1)維廣義淺水波方程的Backlund變換和新精確解的構(gòu)建
    極限海況下單點系泊系統(tǒng)纜索動張力研究
    找不同
    中国国产av一级| 中文乱码字字幕精品一区二区三区 | 男女边吃奶边做爰视频| 久久精品人妻少妇| 我的老师免费观看完整版| 国产老妇伦熟女老妇高清| 一级毛片久久久久久久久女| 国产亚洲91精品色在线| 亚洲精品成人久久久久久| 国产免费福利视频在线观看| 国产午夜精品论理片| 日韩一区二区三区影片| 久久午夜福利片| 深爱激情五月婷婷| 一级片'在线观看视频| eeuss影院久久| 久久99热这里只有精品18| 亚洲欧洲日产国产| 卡戴珊不雅视频在线播放| 日本欧美国产在线视频| 精品99又大又爽又粗少妇毛片| 一二三四中文在线观看免费高清| 99久久精品热视频| 亚洲精品国产成人久久av| 亚洲精品久久午夜乱码| 亚洲国产色片| 亚洲精品乱码久久久v下载方式| 久久精品国产亚洲av天美| 久久精品熟女亚洲av麻豆精品 | 免费av观看视频| 国内精品一区二区在线观看| 欧美xxⅹ黑人| 久久精品久久久久久噜噜老黄| 男人舔奶头视频| 久久久久精品久久久久真实原创| 国产黄色免费在线视频| 美女主播在线视频| av一本久久久久| 好男人视频免费观看在线| 韩国高清视频一区二区三区| 极品少妇高潮喷水抽搐| 国产成人freesex在线| 大香蕉久久网| 十八禁网站网址无遮挡 | 秋霞在线观看毛片| 高清在线视频一区二区三区| 久久久久国产网址| 亚洲av国产av综合av卡| 国产一区二区亚洲精品在线观看| av在线天堂中文字幕| 欧美xxxx黑人xx丫x性爽| 免费看av在线观看网站| 看黄色毛片网站| av天堂中文字幕网| 国产一区二区三区av在线| av免费在线看不卡| 99久久中文字幕三级久久日本| 在线观看一区二区三区| 中国国产av一级| 男人爽女人下面视频在线观看| 国产乱来视频区| 午夜久久久久精精品| 大又大粗又爽又黄少妇毛片口| 成人欧美大片| 成人高潮视频无遮挡免费网站| 熟妇人妻久久中文字幕3abv| 欧美日韩精品成人综合77777| 91久久精品国产一区二区三区| 嫩草影院精品99| 偷拍熟女少妇极品色| 欧美xxxx黑人xx丫x性爽| 婷婷色av中文字幕| 日韩成人伦理影院| 国产黄色免费在线视频| 少妇的逼好多水| av在线播放精品| 水蜜桃什么品种好| 精品一区二区三卡| 三级男女做爰猛烈吃奶摸视频| 欧美激情久久久久久爽电影| 亚洲精品日韩在线中文字幕| 久久久久久国产a免费观看| 久久韩国三级中文字幕| 三级毛片av免费| 国产在线一区二区三区精| 亚洲精品一二三| 亚洲精品aⅴ在线观看| 九九久久精品国产亚洲av麻豆| 国内揄拍国产精品人妻在线| 毛片一级片免费看久久久久| 国产精品一二三区在线看| 国产精品久久久久久久久免| 免费av观看视频| 国产高潮美女av| 高清在线视频一区二区三区| 日本免费在线观看一区| 国产永久视频网站| 蜜臀久久99精品久久宅男| 亚洲精品日韩在线中文字幕| 成人亚洲精品一区在线观看 | 亚洲欧洲国产日韩| 在线免费观看的www视频| 最新中文字幕久久久久| 日韩欧美精品v在线| 寂寞人妻少妇视频99o| 国产成人精品久久久久久| 中文资源天堂在线| 老女人水多毛片| 国产伦在线观看视频一区| 亚洲av二区三区四区| 亚洲乱码一区二区免费版| 国产成人精品婷婷| 卡戴珊不雅视频在线播放| 女人被狂操c到高潮| 精品国产露脸久久av麻豆 | 最近中文字幕2019免费版| 精华霜和精华液先用哪个| 欧美另类一区| 国产精品久久久久久精品电影小说 | 国产免费视频播放在线视频 | 九九爱精品视频在线观看| 一个人看视频在线观看www免费| 免费av观看视频| 精品一区二区三区人妻视频| 免费电影在线观看免费观看| 久久国内精品自在自线图片| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av天美| 欧美高清成人免费视频www| 夫妻性生交免费视频一级片| 在线观看一区二区三区| 最近视频中文字幕2019在线8| 综合色丁香网| 99热6这里只有精品| 蜜桃久久精品国产亚洲av| 免费av观看视频| 国产精品综合久久久久久久免费| 亚洲精品亚洲一区二区| 日本猛色少妇xxxxx猛交久久| 国产精品爽爽va在线观看网站| 国产黄色小视频在线观看| av国产久精品久网站免费入址| 嫩草影院入口| 欧美xxxx性猛交bbbb| 观看免费一级毛片| 成年版毛片免费区| 久久久久精品久久久久真实原创| 国产精品爽爽va在线观看网站| 天美传媒精品一区二区| 国产69精品久久久久777片| 精品久久久久久久人妻蜜臀av| av在线亚洲专区| 国内少妇人妻偷人精品xxx网站| 五月伊人婷婷丁香| 少妇猛男粗大的猛烈进出视频 | 校园人妻丝袜中文字幕| 成年人午夜在线观看视频 | 免费看光身美女| 一区二区三区高清视频在线| 久久精品夜夜夜夜夜久久蜜豆| 黑人高潮一二区| 91在线精品国自产拍蜜月| 麻豆精品久久久久久蜜桃| 美女被艹到高潮喷水动态| 麻豆久久精品国产亚洲av| 亚洲av成人精品一二三区| 大香蕉97超碰在线| 精品一区二区三区人妻视频| 国产午夜精品久久久久久一区二区三区| 日韩av在线免费看完整版不卡| 日韩精品青青久久久久久| 日韩大片免费观看网站| 午夜免费激情av| 狠狠精品人妻久久久久久综合| 亚洲在线自拍视频| 身体一侧抽搐| 有码 亚洲区| 国产淫片久久久久久久久| 成年女人在线观看亚洲视频 | 成人高潮视频无遮挡免费网站| 亚洲激情五月婷婷啪啪| av一本久久久久| 久久人人爽人人片av| 亚洲欧美成人精品一区二区| av播播在线观看一区| 一个人免费在线观看电影| 国产单亲对白刺激| 色网站视频免费| 男女啪啪激烈高潮av片| av国产久精品久网站免费入址| 欧美极品一区二区三区四区| 成年女人看的毛片在线观看| 亚洲欧美成人综合另类久久久| 国产av不卡久久| 激情 狠狠 欧美| 国模一区二区三区四区视频| 久久草成人影院| 欧美+日韩+精品| 成人毛片a级毛片在线播放| 最近最新中文字幕免费大全7| 欧美性猛交╳xxx乱大交人| 成年人午夜在线观看视频 | 舔av片在线| 亚洲va在线va天堂va国产| 成人性生交大片免费视频hd| 亚洲精品456在线播放app| 免费少妇av软件| 中文字幕久久专区| 人妻一区二区av| 午夜福利高清视频| 亚洲成人中文字幕在线播放| 狠狠精品人妻久久久久久综合| 少妇的逼水好多| 国产爱豆传媒在线观看| 成人特级av手机在线观看| 国产精品久久久久久精品电影小说 | 夜夜看夜夜爽夜夜摸| 亚洲精品色激情综合| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区三区四区免费观看| 亚洲国产欧美在线一区| 内地一区二区视频在线| 自拍偷自拍亚洲精品老妇| 亚洲婷婷狠狠爱综合网| 人妻一区二区av| 国产一区亚洲一区在线观看| 免费看av在线观看网站| 一区二区三区高清视频在线| av在线播放精品| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 国产黄a三级三级三级人| 毛片女人毛片| 亚洲,欧美,日韩| 欧美bdsm另类| 美女内射精品一级片tv| 联通29元200g的流量卡| 2022亚洲国产成人精品| 国产综合精华液| 亚洲久久久久久中文字幕| 插逼视频在线观看| 色综合站精品国产| 中文字幕亚洲精品专区| 亚洲电影在线观看av| 国产亚洲91精品色在线| 亚洲一区高清亚洲精品| 插逼视频在线观看| 国产国拍精品亚洲av在线观看| 高清av免费在线| 麻豆精品久久久久久蜜桃| a级一级毛片免费在线观看| 免费高清在线观看视频在线观看| 最新中文字幕久久久久| 国语对白做爰xxxⅹ性视频网站| 婷婷色av中文字幕| 亚洲欧美精品自产自拍| 乱人视频在线观看| 欧美成人精品欧美一级黄| 国产精品一区二区三区四区免费观看| 久久久午夜欧美精品| 丰满乱子伦码专区| 亚洲av电影不卡..在线观看| 国产老妇女一区| www.色视频.com| 人妻少妇偷人精品九色| 国产精品不卡视频一区二区| 国产精品一及| 日本一二三区视频观看| 国产亚洲精品av在线| 婷婷色麻豆天堂久久| 久久99热这里只频精品6学生| 综合色丁香网| 少妇熟女欧美另类| 天天一区二区日本电影三级| 日韩制服骚丝袜av| 久久韩国三级中文字幕| 日本免费在线观看一区| 禁无遮挡网站| 日韩av在线免费看完整版不卡| 天堂中文最新版在线下载 | 99久久精品一区二区三区| 日韩av在线大香蕉| 麻豆国产97在线/欧美| 国国产精品蜜臀av免费| 草草在线视频免费看| 亚洲精品乱码久久久久久按摩| 别揉我奶头 嗯啊视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲最大成人中文| 成年女人在线观看亚洲视频 | 国产精品一区二区三区四区免费观看| 美女高潮的动态| av线在线观看网站| 大又大粗又爽又黄少妇毛片口| 69av精品久久久久久| av福利片在线观看| 中文字幕人妻熟人妻熟丝袜美| 啦啦啦啦在线视频资源| 日本三级黄在线观看| 免费看日本二区| 一级黄片播放器| 中文字幕制服av| 欧美日本视频| 亚洲色图av天堂| 国产免费福利视频在线观看| 亚洲精品久久久久久婷婷小说| 国产有黄有色有爽视频| 国内精品美女久久久久久| 永久网站在线| 欧美xxxx性猛交bbbb| 美女xxoo啪啪120秒动态图| 天美传媒精品一区二区| 精品欧美国产一区二区三| 最后的刺客免费高清国语| 欧美日韩视频高清一区二区三区二| 搡女人真爽免费视频火全软件| 亚洲欧美精品自产自拍| 精品欧美国产一区二区三| 国产亚洲91精品色在线| 国产午夜精品久久久久久一区二区三区| 十八禁网站网址无遮挡 | 国精品久久久久久国模美| 亚洲成人久久爱视频| 一级黄片播放器| 国产精品蜜桃在线观看| 99热这里只有是精品50| 色尼玛亚洲综合影院| 丝瓜视频免费看黄片| 中文字幕制服av| 乱码一卡2卡4卡精品| 午夜福利网站1000一区二区三区| 亚洲色图av天堂| 日韩欧美精品免费久久| 男女那种视频在线观看| 深夜a级毛片| 少妇熟女欧美另类| 亚洲精品乱码久久久久久按摩| 18禁在线播放成人免费| 欧美高清成人免费视频www| 爱豆传媒免费全集在线观看| 久久精品国产鲁丝片午夜精品| 久久久亚洲精品成人影院| 亚洲四区av| 日日啪夜夜爽| 菩萨蛮人人尽说江南好唐韦庄| 97精品久久久久久久久久精品| 亚洲色图av天堂| 国产av国产精品国产| 国产又色又爽无遮挡免| 亚洲av成人av| 黄色欧美视频在线观看| 一边亲一边摸免费视频| 成人欧美大片| 日本猛色少妇xxxxx猛交久久| 观看免费一级毛片| 亚州av有码| 男女那种视频在线观看| 国产成人精品婷婷| 国产伦精品一区二区三区四那| 亚洲国产精品专区欧美| 天堂影院成人在线观看| 亚洲av男天堂| 国产成人精品一,二区| 亚洲成人久久爱视频| 在线a可以看的网站| 两个人视频免费观看高清| 亚洲av成人精品一区久久| 男的添女的下面高潮视频| 免费大片18禁| 国产中年淑女户外野战色| 人妻制服诱惑在线中文字幕| 九九久久精品国产亚洲av麻豆| 亚洲aⅴ乱码一区二区在线播放| 又爽又黄无遮挡网站| 国产成人精品福利久久| 国产成人精品婷婷| 亚洲欧美一区二区三区国产| 少妇熟女aⅴ在线视频| 亚洲av福利一区| 三级国产精品片| av福利片在线观看| 亚洲国产欧美人成| 精品久久久久久成人av| av网站免费在线观看视频 | 最近中文字幕高清免费大全6| 久久久久久久国产电影| 永久免费av网站大全| 舔av片在线| 老师上课跳d突然被开到最大视频| 亚洲精品aⅴ在线观看| 亚洲av不卡在线观看| ponron亚洲| 久久久久精品性色| 日韩av不卡免费在线播放| 亚洲成人久久爱视频| 在线观看av片永久免费下载| 国产精品美女特级片免费视频播放器| av网站免费在线观看视频 | 人妻制服诱惑在线中文字幕| 一个人免费在线观看电影| 久久精品国产亚洲av涩爱| 欧美一区二区亚洲| videos熟女内射| 一级毛片 在线播放| 亚洲av成人精品一区久久| 亚洲国产欧美在线一区| 赤兔流量卡办理| av在线天堂中文字幕| 久久精品国产亚洲av涩爱| 久久97久久精品| 大香蕉久久网| 亚洲欧美清纯卡通| ponron亚洲| 麻豆乱淫一区二区| 亚洲熟女精品中文字幕| 午夜福利成人在线免费观看| 高清日韩中文字幕在线| 久久久a久久爽久久v久久| 欧美精品国产亚洲| 免费av观看视频| 观看美女的网站| 日韩,欧美,国产一区二区三区| 国产综合精华液| 欧美人与善性xxx| 免费看美女性在线毛片视频| 天天一区二区日本电影三级| 欧美日韩国产mv在线观看视频 | 天堂av国产一区二区熟女人妻| 亚洲国产av新网站| 亚洲国产成人一精品久久久| 色吧在线观看| 欧美变态另类bdsm刘玥| 免费看美女性在线毛片视频| 国产精品一区二区三区四区免费观看| 亚洲精品影视一区二区三区av| 欧美激情久久久久久爽电影| 一区二区三区免费毛片| 两个人的视频大全免费| 非洲黑人性xxxx精品又粗又长| 国产极品天堂在线| 91精品国产九色| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 一级片'在线观看视频| 国产成人午夜福利电影在线观看| 我要看日韩黄色一级片| 国产一区二区三区综合在线观看 | 九九久久精品国产亚洲av麻豆| 99久久人妻综合| av又黄又爽大尺度在线免费看| 99热这里只有是精品50| 在线观看一区二区三区| 亚洲av在线观看美女高潮| 国产成年人精品一区二区| 80岁老熟妇乱子伦牲交| 啦啦啦中文免费视频观看日本| 大香蕉久久网| 亚洲最大成人手机在线| 日产精品乱码卡一卡2卡三| 欧美日韩亚洲高清精品| 国产精品国产三级国产av玫瑰| 国产成人精品一,二区| 成人特级av手机在线观看| 一个人看视频在线观看www免费| 久久97久久精品| 国产黄片美女视频| 一区二区三区乱码不卡18| 精品一区二区三区人妻视频| 天堂俺去俺来也www色官网 | 久久久久精品性色| 欧美日韩在线观看h| 熟妇人妻不卡中文字幕| 天堂√8在线中文| 小蜜桃在线观看免费完整版高清| 舔av片在线| 国产毛片a区久久久久| 日本一二三区视频观看| 91午夜精品亚洲一区二区三区| 99热全是精品| 日本wwww免费看| 联通29元200g的流量卡| 一级a做视频免费观看| 搞女人的毛片| 久久6这里有精品| 国产精品一区二区三区四区免费观看| 亚洲欧美日韩卡通动漫| 寂寞人妻少妇视频99o| 亚洲精品第二区| 免费观看a级毛片全部| 国产亚洲av片在线观看秒播厂 | 狠狠精品人妻久久久久久综合| 菩萨蛮人人尽说江南好唐韦庄| av在线老鸭窝| 一级二级三级毛片免费看| 国产精品1区2区在线观看.| 99久久精品国产国产毛片| 婷婷色麻豆天堂久久| 熟妇人妻不卡中文字幕| 床上黄色一级片| 国产成人freesex在线| 春色校园在线视频观看| 国产黄色免费在线视频| 99热这里只有精品一区| 国产成人福利小说| 亚洲在久久综合| 美女内射精品一级片tv| 男女边摸边吃奶| 晚上一个人看的免费电影| 高清日韩中文字幕在线| 国产亚洲精品av在线| 国产精品久久视频播放| 国产爱豆传媒在线观看| 亚洲久久久久久中文字幕| 国产国拍精品亚洲av在线观看| av专区在线播放| 日本免费a在线| 在线天堂最新版资源| 国产乱来视频区| 精品一区二区免费观看| 久久久精品免费免费高清| 69av精品久久久久久| 国产免费福利视频在线观看| 蜜桃亚洲精品一区二区三区| 国产 亚洲一区二区三区 | 亚洲人成网站高清观看| 久久久精品欧美日韩精品| 国产伦在线观看视频一区| 欧美xxxx黑人xx丫x性爽| 亚洲精品自拍成人| 亚洲不卡免费看| 国国产精品蜜臀av免费| 大香蕉久久网| 亚洲成人久久爱视频| 日韩精品有码人妻一区| 别揉我奶头 嗯啊视频| 男女边吃奶边做爰视频| 成人亚洲精品一区在线观看 | 免费看不卡的av| av免费在线看不卡| 女人被狂操c到高潮| 成人毛片a级毛片在线播放| 一本久久精品| 精品不卡国产一区二区三区| 一个人观看的视频www高清免费观看| 高清视频免费观看一区二区 | 中文在线观看免费www的网站| 亚洲精品日韩av片在线观看| 国产亚洲一区二区精品| 日本与韩国留学比较| 97精品久久久久久久久久精品| 免费观看在线日韩| 国内精品宾馆在线| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产鲁丝片午夜精品| 久久99蜜桃精品久久| 美女高潮的动态| 五月伊人婷婷丁香| 欧美高清性xxxxhd video| 又大又黄又爽视频免费| 久久草成人影院| 亚洲av免费在线观看| 亚洲精品一二三| 少妇猛男粗大的猛烈进出视频 | 黄片wwwwww| 午夜免费观看性视频| 精品久久久久久久末码| 蜜臀久久99精品久久宅男| 丝袜喷水一区| 91在线精品国自产拍蜜月| 国产精品综合久久久久久久免费| 亚洲精品中文字幕在线视频 | 欧美3d第一页| 五月玫瑰六月丁香| 欧美激情久久久久久爽电影| 国产麻豆成人av免费视频| 免费av毛片视频| 亚洲,欧美,日韩| 中文乱码字字幕精品一区二区三区 | 久久久久九九精品影院| 亚洲av中文av极速乱| 国产 一区 欧美 日韩| 免费大片18禁| 国产女主播在线喷水免费视频网站 | 久久久久久久久久久免费av| 亚洲精品成人av观看孕妇| 亚洲一级一片aⅴ在线观看| 99re6热这里在线精品视频| 国产精品嫩草影院av在线观看| 国产精品国产三级专区第一集| 久久精品国产鲁丝片午夜精品| av一本久久久久| 欧美日韩一区二区视频在线观看视频在线 | 国产免费视频播放在线视频 | 如何舔出高潮| 久久久精品94久久精品| 国产精品久久久久久久电影| 国产精品蜜桃在线观看| 床上黄色一级片| av天堂中文字幕网| 青春草国产在线视频| 又粗又硬又长又爽又黄的视频| 黄色配什么色好看| 亚洲精品成人av观看孕妇| 亚洲av免费在线观看| 国产精品.久久久| 在线免费观看不下载黄p国产| 久久精品夜夜夜夜夜久久蜜豆| 天堂av国产一区二区熟女人妻| 六月丁香七月| 久久99精品国语久久久| 久久久a久久爽久久v久久| 女人十人毛片免费观看3o分钟| 日本黄大片高清| 亚洲精品乱码久久久v下载方式| 国产精品爽爽va在线观看网站|