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

    考慮到臺風(fēng)海況浮式風(fēng)力機半潛式平臺風(fēng)浪流載荷動態(tài)響應(yīng)及系泊性能分析

    2017-12-19 02:54:36黃致謙丁勤衛(wèi)周紅杰陳福東
    動力工程學(xué)報 2017年12期
    關(guān)鍵詞:潛式海況浮式

    黃致謙, 李 春,2, 丁勤衛(wèi), 周紅杰, 陳福東

    (1.上海理工大學(xué) 能源與動力工程學(xué)院,上海 200093; 2.上海市動力工程多相流動與傳熱重點實驗室,上海 200093)

    考慮到臺風(fēng)海況浮式風(fēng)力機半潛式平臺風(fēng)浪流載荷動態(tài)響應(yīng)及系泊性能分析

    黃致謙1, 李 春1,2, 丁勤衛(wèi)1, 周紅杰1, 陳福東1

    (1.上海理工大學(xué) 能源與動力工程學(xué)院,上海 200093; 2.上海市動力工程多相流動與傳熱重點實驗室,上海 200093)

    為研究浮式風(fēng)力機半潛式平臺的動態(tài)響應(yīng)及系泊性能,建立了基于半潛式平臺的NREL 5 MW浮式風(fēng)力機模型.結(jié)合有限元方法,并利用輻射和衍射理論,在考慮風(fēng)、浪、流載荷的聯(lián)合作用下,運用Aqwa軟件對極端海況下半潛式平臺的動態(tài)響應(yīng)及系泊性能進行數(shù)值仿真,得到了幅值響應(yīng)算子(RAO)、附加質(zhì)量、輻射阻尼隨波浪頻率的變化趨勢以及極端海況下半潛式平臺的動態(tài)響應(yīng)和系泊線的張力響應(yīng).結(jié)果表明:半潛式平臺在高頻區(qū)域(大于1.6 rad/s)具有良好的運動性能,動態(tài)響應(yīng)?。淮故幏较騽討B(tài)響應(yīng)不受波浪方向變化的影響,且波浪方向為0°時,縱蕩和縱搖動態(tài)響應(yīng)最為劇烈;臺風(fēng)海況下半潛式平臺的動態(tài)響應(yīng)和系泊線的張力響應(yīng)均大于無臺風(fēng)海況;隨著海況惡劣程度加深,半潛式平臺的動態(tài)響應(yīng)和系泊線的張力響應(yīng)均明顯增大.

    浮式風(fēng)力機; 半潛式平臺; 臺風(fēng)海況; 系泊系統(tǒng); 動態(tài)響應(yīng)

    化石能源支撐著世界經(jīng)濟發(fā)展,但現(xiàn)今化石能源日益枯竭,各國經(jīng)濟發(fā)展均受到能源緊缺等問題的限制,人們越來越多地開始關(guān)注能源問題[1-2].風(fēng)能是太陽能的一種轉(zhuǎn)換形式,是取之不盡、用之不竭的可再生能源,具有蘊藏量大、無污染、可再生、分布廣等特點[3].與陸上風(fēng)能相比,海上風(fēng)能具有充裕、平穩(wěn)、海面摩擦系數(shù)小、風(fēng)能密度大等優(yōu)點[4].相比陸上風(fēng)電場,海上風(fēng)電場的建設(shè)不存在噪音污染問題,且不會占用寶貴的陸地資源.根據(jù)支撐風(fēng)力機的塔架形式,海上風(fēng)力機分為固定式和漂浮式.隨著大型現(xiàn)代風(fēng)力機技術(shù)的成熟和海上風(fēng)力發(fā)電的發(fā)展,利用浮式平臺技術(shù)獲取海上優(yōu)質(zhì)風(fēng)資源已成為公認(rèn)的發(fā)展方向和研究熱點[5].半潛式平臺因總投資較小、甲板空間及甲板載荷可變范圍大、抗風(fēng)浪能力強、適應(yīng)水深范圍廣和作業(yè)功能全面等優(yōu)勢,具有較好的應(yīng)用前景[6].相比于固定式風(fēng)力機,浮式風(fēng)力機增加了浮式基礎(chǔ)和錨泊系統(tǒng),還承受復(fù)雜的風(fēng)、浪、流載荷影響,導(dǎo)致外界載荷條件更加復(fù)雜,其搖蕩特性對風(fēng)電機組的發(fā)電性能也有較大影響.因此,研究在風(fēng)、浪、流載荷作用下浮式平臺的動態(tài)響應(yīng)非常必要[7].

    1972年,Heronemus[8]率先提出具有浮式基礎(chǔ)的海上漂浮式風(fēng)力機(以下簡稱浮式風(fēng)力機)概念后,國內(nèi)外學(xué)者開始對半潛式平臺等進行探索.2000年,S?ylemez等[9]建立半潛式平臺模型,將莫里森方程法和格林函數(shù)積分方程法2種實用方法計算所得半潛式平臺波浪力和動態(tài)響應(yīng)進行比較,并將其與實驗結(jié)果對比,發(fā)現(xiàn)在主要頻域范圍內(nèi)2種方法吻合良好.2006年, Zambrano等[10]運用Wamit軟件計算了作用在半潛式平臺上的波浪力,并分析了在墨西哥風(fēng)暴海況下半潛式平臺6個自由度的動態(tài)響應(yīng).2008年,Manabe等[11]提出一種基于半潛式船體(semi-sub hull)結(jié)構(gòu)的浮式風(fēng)力機系統(tǒng),分析了浮體結(jié)構(gòu)的水動力響應(yīng),并優(yōu)化了該浮體結(jié)構(gòu).2012年,Lefebvre等[12]對蘇格蘭北海的浮式風(fēng)力機進行調(diào)查,提出了一種改進的半潛式平臺浮式風(fēng)力機平臺概念,并利用Sesam軟件對其整體穩(wěn)定性、經(jīng)濟性和技術(shù)性進行了分析和比較.Tang等[13]設(shè)計了60 m水深的600 kW半潛式平臺海上水平軸風(fēng)力機,對包括系泊系統(tǒng)、浮式基礎(chǔ)、塔架及風(fēng)力機在內(nèi)的結(jié)構(gòu)建立了精確的數(shù)學(xué)模型,并分析了在波浪載荷影響下其結(jié)構(gòu)的頻域動態(tài)響應(yīng).張亮等[14]針對其課題組設(shè)計的一種半潛式平臺進行了穩(wěn)定性分析,采用數(shù)值模擬方法分別計算了該半潛式平臺在不同工況下的力矩曲線.肖宇維等[15]通過Ansys/Aqwa軟件計算了規(guī)則波作用下不同形式輕型張力腿半潛式平臺6個自由度的動態(tài)響應(yīng),并在不規(guī)則波作用下進行了數(shù)值分析,得到時域分析結(jié)果.上述研究大多僅針對浮式平臺的動態(tài)響應(yīng)或系泊性能,較少有人綜合考慮浮式平臺的動態(tài)響應(yīng)和系泊性能,且在對浮式平臺進行動態(tài)響應(yīng)或系泊性能研究時,大多采用單一載荷作用.

    目前主要根據(jù)海洋石油平臺的相關(guān)規(guī)范和經(jīng)驗進行探索,尚缺乏有關(guān)浮式風(fēng)力機平臺所受載荷及穩(wěn)定性的設(shè)計標(biāo)準(zhǔn)[16].平臺結(jié)構(gòu)穩(wěn)定是浮式風(fēng)力機安全運行的根本保障[4],而動態(tài)響應(yīng)特性是浮式平臺相關(guān)研究中的重點和難點.動態(tài)響應(yīng)特性不僅是結(jié)構(gòu)形式創(chuàng)新及設(shè)計優(yōu)化的重要內(nèi)容,同時也為平臺的施工、安裝及檢測提供了支持[1].因此,筆者研究了浮式風(fēng)力機半潛式平臺的動態(tài)響應(yīng)及系泊性能,以期為浮式風(fēng)力機平臺的設(shè)計提供參考.

    1 浮式風(fēng)力機參數(shù)及系泊系統(tǒng)

    1.1 半潛式平臺及風(fēng)力機參數(shù)

    以半潛式平臺為研究對象,半潛式平臺選為Semi-submersible[17],浮式風(fēng)力機選為NREL 5 MW[18].半潛式平臺的詳細(xì)參數(shù)見表1,浮式風(fēng)力機的詳細(xì)參數(shù)見表2.

    基于NREL 5 MW風(fēng)力機參數(shù)和半潛式平臺參數(shù)建立的浮式風(fēng)力機模型如圖1所示.模型的網(wǎng)格分布如圖2所示.計算要求1個波長至少要覆蓋7個最大網(wǎng)格單元,但網(wǎng)格越細(xì),可計算的波浪頻率越大,計算耗時越多.故劃分網(wǎng)格時,需根據(jù)計算的最大波浪頻率來設(shè)定網(wǎng)格的控制尺寸.

    表1 半潛式平臺參數(shù)

    表2 NREL 5 MW風(fēng)力機參數(shù)

    圖1 浮式風(fēng)力機模型

    1.2 系泊系統(tǒng)

    根據(jù)API RP 2SK規(guī)范,懸鏈線式系泊系統(tǒng)采用“錨鏈(上段)-鋼纜(中段)-錨鏈(下段)”布置形式,具有重量小、經(jīng)濟性好等優(yōu)點,具體設(shè)計參數(shù)如表3所示.

    圖2 網(wǎng)格分布

    表3 系泊線參數(shù)

    設(shè)計系泊系統(tǒng)時,系泊線為6根,布置方式為3×2對稱布置,具體系泊線布置形式見圖3和圖4.

    圖3 系泊系統(tǒng)示意圖

    2 計算流程、環(huán)境載荷及參數(shù)

    2.1 計算流程

    應(yīng)用Aqwa軟件對半潛式平臺的動態(tài)響應(yīng)及極端海況下的系泊性能進行數(shù)值模擬和分析,計算時采用Aqwa-Line模塊和Aqwa-Drift模塊.

    圖4 局部系泊線水下示意圖

    Aqwa-Line模塊中,可通過輻射和衍射理論計算浮體結(jié)構(gòu)周圍的波浪力,也可獲得附加質(zhì)量、輻射阻尼及頻域動態(tài)響應(yīng).

    Aqwa-Drift模塊中,可計算在風(fēng)、浪、流載荷的聯(lián)合作用下浮體結(jié)構(gòu)的動態(tài)響應(yīng)時間歷程和系泊系統(tǒng)的張力響應(yīng)時間歷程.

    Aqwa-Drift模塊可使用Aqwa-Line模塊中的附加質(zhì)量力、輻射阻尼力及各波長與波浪方向漂移力等數(shù)據(jù),可得到風(fēng)、浪、流載荷聯(lián)合作用下的浮體結(jié)構(gòu)運動及系泊系統(tǒng)張力響應(yīng)的時間歷程,計算流程圖如圖5所示.

    圖5 計算流程圖

    2.2 風(fēng)載荷

    在船舶與海洋平臺等工程設(shè)計過程中,風(fēng)載荷是必須考慮的基本載荷之一.風(fēng)載荷作用在半潛式平臺上,對其產(chǎn)生水平力和力矩,會影響漂浮穩(wěn)定性及半潛式平臺的結(jié)構(gòu)強度.風(fēng)速及其方向通常隨時間變化,但在較長時間段內(nèi),風(fēng)速可作固定值(一般取海平面以上平均風(fēng)速為10 m/h).假定風(fēng)速不變,且沿y軸負(fù)方向,垂直風(fēng)力機方向入射.由于設(shè)定為惡劣海況,風(fēng)速已超過風(fēng)力機的切出風(fēng)速,風(fēng)力機處于停機狀態(tài),故不考慮風(fēng)力機葉片旋轉(zhuǎn)產(chǎn)生的載荷.作用在半潛式平臺上的風(fēng)載荷按下式計算:

    (1)

    Fwi=CwiV2

    (2)

    式中:Cwi為i自由度的風(fēng)力系數(shù);CS為受風(fēng)構(gòu)件形狀系數(shù);Ch為風(fēng)壓沿高度變化的系數(shù);Sj為受風(fēng)結(jié)構(gòu)在風(fēng)向上的投影面積;Fwi為i自由度的風(fēng)力;V為平均風(fēng)速.

    2.3 浪載荷

    與陸上風(fēng)力機不同,作為海洋能源開發(fā)的重要裝備,半潛式平臺風(fēng)力機長期在惡劣海況中作業(yè)并承受多種載荷作用,且受浪載荷的影響,這也是工程設(shè)計中關(guān)心的重要問題之一.

    采用線性勢流理論計算浪載荷及運動,流體為不可壓縮、無黏、無旋.整個流域為由浮體濕表面SB、波浪自由表面SF和海底表面SD構(gòu)成的半無限大空間,如圖6所示.

    圖6 波浪對物體作用示意圖

    總速度勢包括6個自由度上的入射波浪勢、繞射波浪勢和輻射波浪勢,可定義為:

    (3)

    式中:φ為總速度勢;φI為入射波浪勢;φD為繞射波浪勢;φR為輻射波浪勢.

    流體的總速度勢滿足拉普拉斯方程及無窮遠(yuǎn)的輻射條件(Sommerfeld條件):

    ▽2φ=0

    (4)

    為使方程有唯一確定解,須滿足定解條件,包括邊界條件和初始條件:

    (5)

    (6)

    (7)

    式中:g為重力加速度,m/s2;n為浮體濕表面外法向量;vn為浮體沿向量n方向的速度,m/s.

    浮體濕表面的水動壓力為:

    (8)

    波浪力Fw和波浪力矩Mw分別為:

    Fw=?SB-pndSB

    (9)

    Mw=?SB-p(r×n)dSB

    (10)

    式中:r為波浪力作用點到重心的距離;ρ為海水密度.

    2.4 流載荷

    海流是海洋環(huán)境中各種流動的集合,是主要動力現(xiàn)象之一.設(shè)計水深為320 m,屬于深海領(lǐng)域,故海流可視為剪切流.在水平方向上,海流近似為恒速;在垂直方向,海流流速呈梯度遞減,海底流速幾乎為0.相對風(fēng)載荷和浪載荷,流載荷較為穩(wěn)定.計算流載荷時,可認(rèn)為海流對結(jié)構(gòu)物的作用力僅為拖拽力Fc:

    Fc=CSCD(Ac+Af)

    (11)

    式中:CD為阻力系數(shù);CS為半潛式平臺的海流系數(shù);Ac為海面下圓柱構(gòu)件的投影面積;Af為海面下平面構(gòu)件的投影面積.

    2.5 運動方程

    整個浮式風(fēng)力機系統(tǒng)在波浪環(huán)境中的運動方程為:

    Fw+Fc+Fi

    (12)

    式中:M為浮式風(fēng)力機系統(tǒng)的質(zhì)量矩陣;C為浮式風(fēng)力機系統(tǒng)的阻尼矩陣;D1和D2分別為線性阻尼矩陣和平方阻尼矩陣;F為矢量函數(shù);K(x)為以位移為函數(shù)的靜水力剛度矩陣;x為位移矢量;Fw為浪載荷;Fi為風(fēng)載荷.

    2.6 環(huán)境參數(shù)

    通常情況下,當(dāng)風(fēng)、浪、流載荷以相同方向入射半潛式平臺時,將對半潛式平臺的系泊系統(tǒng)產(chǎn)生最大影響.因此, Aqwa-Drift模塊上設(shè)置風(fēng)、浪、流載荷以相同方向入射半潛式平臺,即沿y軸負(fù)方向、垂直于半潛式平臺入射.時域分析選取南海有、無臺風(fēng)2種海況,具體參數(shù)如表4所示.

    表4 不同海況參數(shù)

    3 模擬結(jié)果及分析

    3.1 半潛式平臺頻域動態(tài)響應(yīng)

    3.1.1RAO頻域動態(tài)響應(yīng)

    浮式風(fēng)力機半潛式平臺有6個自由度,分別為沿x軸、y軸和z軸的平動及繞各軸的轉(zhuǎn)動.其中,由長度單位表示的平動包括縱蕩、橫蕩和垂蕩;由角度單位表示的轉(zhuǎn)動包括橫搖、縱搖和首搖.浮式風(fēng)力機半潛式平臺6個自由度的運動如圖7所示.

    圖7 半潛式平臺6個自由度的運動示意圖

    受規(guī)則波作用,浮體結(jié)構(gòu)產(chǎn)生6個自由度波頻動態(tài)響應(yīng)的穩(wěn)態(tài)幅值,稱為幅值響應(yīng)算子(RAO).RAO可反映不同海況對浮體結(jié)構(gòu)的影響,表現(xiàn)為浮體結(jié)構(gòu)受規(guī)則波影響下隨頻率變化的一階響應(yīng)幅值:

    RAO=ηi/ξ

    (13)

    式中:ηi為半潛式平臺運動第i個自由度的值;ξ為某一頻率波浪高度的幅值.

    通常認(rèn)為海洋中波浪高度變化是一個具有零均值、各態(tài)歷經(jīng)的高斯隨機過程,浮體結(jié)構(gòu)受任意波浪組分作用的響應(yīng)是該組分波波幅的線性函數(shù),且與對其他波浪組分的響應(yīng)無關(guān),故可根據(jù)所有構(gòu)成波浪的各組分單元波響應(yīng)之和得到浮體結(jié)構(gòu)對隨機波浪的響應(yīng).因此,可用半潛式平臺位置RAO來反映在外界載荷影響下半潛式平臺在6個自由度的運動情況.

    以x軸負(fù)方向為0°,逆時針為正,不同波浪方向下(以下簡稱浪向)半潛式平臺縱蕩、垂蕩和縱搖方向RAO隨波浪頻率的變化如圖8所示,其中縱軸表示對應(yīng)單位波幅的位移或偏轉(zhuǎn)角.

    由圖8(a)可知,在縱蕩方向上,浪向為0°和45°時兩者雖峰值不同,但RAO均隨波浪頻率的增大而減小,尤其當(dāng)01.6 rad/s時,RAO逐漸趨于0;浪向為90°時,RAO很小,可以忽略.由圖8(b)可知,在垂蕩方向上,在不同浪向上各RAO曲線幾乎重合,均隨波浪頻率的增大先小范圍上升、后急劇下降,最后趨于平穩(wěn),說明垂蕩方向上RAO不受浪向變化的影響.由圖8(c)可知,縱搖方向與縱蕩方向類似,浪向為90°時,RAO很小,可以忽略;浪向為0°和45°時,兩者變化幅度大體相同,但峰值不同.總體上看,在波浪低頻區(qū)域,半潛式平臺的動態(tài)響應(yīng)特別劇烈;在波浪高頻區(qū)域,當(dāng)f>1.6 rad/s時,半潛式平臺運動性能表現(xiàn)良好,響應(yīng)幅值逐漸趨于0.

    (a)縱蕩

    (b)垂蕩

    (c)縱搖

    3.1.2 水動力頻域響應(yīng)

    在非定常運動過程中,浮式結(jié)構(gòu)物與周圍的水相互作用,可使浮式結(jié)構(gòu)物獲得加速度及速度.水對物體產(chǎn)生的反作用力包括與加速度成正比的附加質(zhì)量和與速度成正比的輻射阻尼,附加質(zhì)量和輻射阻尼是重要的水動力學(xué)參數(shù).

    圖9為各自由度的附加質(zhì)量和輻射阻尼.由圖9(a)可知,垂蕩方向上的附加質(zhì)量大于橫蕩和縱蕩方向上,且其變化最后趨于平穩(wěn).由于半潛式平臺為對稱結(jié)構(gòu),橫蕩和縱蕩方向的附加質(zhì)量曲線重合.在常見的波浪周期(0.8 rad/s1.65 rad/s時附加質(zhì)量趨于穩(wěn)定;首搖方向上,0.03 rad/s1.05 rad/s時附加質(zhì)量小于橫搖和縱搖方向.由圖9(c)可知,縱蕩和垂蕩方向的輻射阻尼完全一致,均隨波浪周期的增大呈先增大后減小的趨勢.與縱蕩和垂蕩方向相比,橫蕩方向的輻射阻尼較小,可忽略.由圖9(d)可知,首搖方向上輻射阻尼的變化趨勢與圖9(c)中縱蕩和垂蕩方向一致,呈先增大后減小的趨勢;橫搖和縱搖方向上輻射阻尼曲線重合,且數(shù)值遠(yuǎn)小于首搖方向,可忽略.

    (a)縱蕩、橫蕩和垂蕩方向的附加質(zhì)量

    (b)縱搖、橫搖和首搖方向的附加質(zhì)量

    (c)縱蕩、橫蕩和垂蕩方向的輻射阻尼

    (d)縱搖、橫搖和首搖方向的輻射阻尼

    3.2 半潛式平臺時域動態(tài)響應(yīng)

    3.2.1RAOs時域動態(tài)響應(yīng)

    描述浮式風(fēng)力機半潛式平臺等浮體結(jié)構(gòu)最直接的參數(shù)是位移量和轉(zhuǎn)動偏轉(zhuǎn)角.在南海有無臺風(fēng)2種極限海況下,風(fēng)、浪、流載荷聯(lián)合作用時半潛式平臺的位移量和轉(zhuǎn)動偏轉(zhuǎn)角如圖10所示.

    (a)橫蕩

    (b)垂蕩

    (c)橫搖

    由圖10可知,在臺風(fēng)海況下半潛式平臺橫蕩、垂蕩和橫搖3個方向的動態(tài)響應(yīng)均比無臺風(fēng)海況下的動態(tài)響應(yīng)劇烈,說明隨著海況惡劣程度的加深,半潛式平臺的動態(tài)響應(yīng)增大.在同一海況下,由于風(fēng)、浪、流載荷沿橫蕩方向入射,因此橫蕩運動的響應(yīng)最大,垂蕩運動次之,橫搖運動最不明顯.

    3.2.2 系泊線的張力響應(yīng)分析

    半潛式平臺系泊系統(tǒng)受上述2種海況的風(fēng)、浪、流載荷作用時,通過Aqwa-Drift計算可得到系泊線的張力響應(yīng).由于相鄰2根系泊線非常接近,受到的載荷作用也極其相似,故僅需分析1號、3號和5號系泊線的張力響應(yīng),其系泊線位置如圖11所示.

    圖11 系泊線位置

    圖12給出了在2種極限海況下1號、3號和5號系泊線張力隨時間的變化曲線.由圖12可知,2種海況下所有系泊線的張力響應(yīng)均達到最大值,且臺風(fēng)海況下系泊線張力響應(yīng)幅度明顯大于無臺風(fēng)海況.相同海況下,由于5號系泊線所處方向為風(fēng)、浪、流載荷正對方向(沿y軸負(fù)方向且垂直于風(fēng)力機),故其張力響應(yīng)最大,所受張力響應(yīng)極值也最大,臺風(fēng)海況時為1 589.31 kN,無臺風(fēng)海況時為581.23 kN;由于1號和3號系泊線處于沿y軸負(fù)方向?qū)ΨQ位置,故其張力響應(yīng)時域曲線基本相同,張力響應(yīng)極值也相近,臺風(fēng)海況時分別為454.83 kN、442.52 kN,無臺風(fēng)海況時分別為391.03 kN、399.92 kN.

    (a)1號系泊線

    (b)3號系泊線

    (c)5號系泊線

    4 結(jié) 論

    (1)在波浪低頻區(qū)域,半潛式平臺的動態(tài)響應(yīng)特別劇烈;在波浪高頻區(qū)域,特別在f>1.6 rad/s區(qū)域,半潛式平臺運動性能表現(xiàn)良好,響應(yīng)幅值逐漸趨于0.

    (2)垂蕩方向的動態(tài)響應(yīng)不受浪向的影響;當(dāng)浪向為0°時,縱蕩和縱搖方向半潛式平臺的動態(tài)響應(yīng)程度最劇烈.

    (3)不同海況下,各自由度的動態(tài)響應(yīng)不同.隨著海況惡劣程度的加深,半潛式平臺的動態(tài)響應(yīng)幅度也增大.

    (4)隨著海況惡劣程度的加劇,各系泊線的張力響應(yīng)顯著增大.

    (5)由于5號系泊線位于風(fēng)、浪、流載荷的入射方向,且半潛式平臺結(jié)構(gòu)對稱,因此相同海況下5號系泊線張力響應(yīng)最大,1號和3號系泊線張力響應(yīng)曲線變化趨勢一致.

    [1] 丁勤衛(wèi), 李春, 周國龍, 等. 基于AQWA的漂浮式風(fēng)力機駁船式平臺的動態(tài)響應(yīng)[J].水資源與水工程學(xué)報, 2015, 26(2): 150-155.

    DING Qinwei, LI Chun, ZHOU Guolong, et al. Dynamic response of the floating wind turbine on barge platform based on AQWA[J].JournalofWaterResourcesandWaterEngineering, 2015, 26(2): 150-155.

    [2] 楊陽, 李春, 繆維跑, 等. 湍流風(fēng)場與地震激勵聯(lián)合作用下的風(fēng)力機結(jié)構(gòu)動力學(xué)響應(yīng)[J].振動與沖擊, 2015, 34(21): 136-143.

    YANG Yang, LI Chun, MIAO Weipao, et al. Response of structural dynamic characteristics of wind turbine operating in turbulent wind combined with seismic motion[J].JournalofVibrationandShock, 2015, 34(21): 136-143.

    [3] 丁勤衛(wèi), 李春, 周國龍, 等. 陸海風(fēng)力機動態(tài)響應(yīng)對比[J].動力工程學(xué)報, 2016, 36(1): 65-73.

    DING Qinwei, LI Chun, ZHOU Guolong, et al. Comparison of dynamic response between stationary and floating wind turbine[J].JournalofChineseSocietyofPowerEngineering, 2016, 36(1): 65-73.

    [4] 于華明, 劉容子, 鮑獻文. 海洋可再生能源發(fā)展現(xiàn)狀與展望[M]. 青島: 中國海洋大學(xué)出版社, 2012.

    [5] 李春, 葉舟, 高偉, 等. 現(xiàn)代大型風(fēng)力機設(shè)計原理[M]. 上海: 上??茖W(xué)技術(shù)出版社, 2013.

    [6] 陳鵬, 馬駿, 黃進浩, 等. 基于AQWA的半潛式平臺水動力分析及系泊性能計算分析[J].船海工程, 2013, 42(3): 44-47.

    CHEN Peng, MA Jun, HUANG Jinhao, et al. Hydrodynamic analysis and mooring system calculation for semi-submersible platform[J].Ship&OceanEngineering, 2013, 42(3): 44-47.

    [7] 王志新. 海上風(fēng)力發(fā)電技術(shù)[M]. 北京: 機械工業(yè)出版社, 2013.

    [8] HERONEMUS W E. Pollution-free energy from the offshore winds[M].Washington D C, USA:Marine Technology Society, 1972.

    [9] S?YLEMEZ M, ATLAR M. Comparative study of two practical methods for estimating the hydrodynamic loads and motions of a semi-submersible[J].JournalofOffshoreMechanicsandArcticEngineering, 2000, 122(1): 57-63.

    [10] ZAMBRANO T, MACCREADY T, KICENIUK T, et al. Dynamic modeling of deepwater offshore wind turbine structures in gulf of mexico storm condition[C]//Proceedingsof25thInternationalConferenceonOffshoreMechanicsandArcticEngineering. Hamburg, Germany: ASME, 2006.

    [11] MANABE H, UEHIRO T, UTIYAMA M, et al. Development of the floating structure for the sailing-type offshore wind farm[C]//OCEANS2008-MTS/IEEEKobeTechno-Ocean. Kobe, Japan: IEEE, 2008.

    [12] LEFEBVRE S, COLLU M. Preliminary design of a floating support structure for a 5 MW offshore wind turbine[J].OceanEngineering, 2012, 40(2): 15-26.

    [13] TANG Y G, HU J, LIU L Q. Study on the dynamic response for floating foundation of offshore wind turbine[C]//TheASME201130thInternationalConferenceonOcean,OffshoreandArcticEngineering.Rotterdam, Netherlands: ASME, 2011.

    [14] 張亮, 鄧慧靜. 浮式風(fēng)機半潛平臺穩(wěn)性數(shù)值分析[J].應(yīng)用科技, 2011, 38(10): 13-17.

    ZHANG Liang, DENG Huijing. Numerical analysis on stability of the semi-submersible platform of floating wind turbines[J].AppliedScienceandTechnology, 2011, 38(10): 13-17.

    [15] 肖宇維, 孫樹民. 基于AQWA的輕型張力腿平臺型式研究[J].科學(xué)技術(shù)與工程, 2012, 12(33): 8987-8989.

    XIAO Weiyu, SUN Shumin. A study of light tension leg platform types based on AQWA[J].ScienceTechnologyandEngineering, 2012, 12(33): 8987-8989.

    [16] 成欣, 葉舟, 周國龍, 等. 漂浮式風(fēng)力機平臺在不同水深下的響應(yīng)特性分析[J].水資源與水工程學(xué)報, 2014, 25(4): 97-103.

    CHENG Xin, YE Zhou, ZHOU Guolong, et al. Analysis on response performance for floating wind turbine platforms under different water depth[J].JournalofWaterResourcesandWaterEngineering, 2014, 25(4): 97-103.

    [17] ROBERTSON A, JONKMAN J, MASCIOLA M, et al. Definition of the semisubmersible floating system for phase II of OC4[R]. National Renewable Energy Laboratory, Washington D C, USA: 2014.

    [18] ROBERTSON A N, JONKMAN J M. Loads analysis of several offshore floating wind turbine concepts[C]//TheTwenty-firstInternationalOffshoreandPolarEngineeringConference. Hawaii, USA: International Society of Offshore and Polar Engineers, 2011.

    DynamicResponseandMooringPerformanceofaSemi-submersiblePlatformforFloatingWindTurbineConsideringTyphoonSeaStates

    HUANGZhiqian1,LIChun1,2,DINGQinwei1,ZHOUHongjie1,CHENFudong1
    (1. School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China; 2. Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering, Shanghai 200093, China)

    To research the dynamic response and mooring performance of a semi-submersible platform for floating wind turbine, a NREL 5 MW wind turbine model was established based on the semi-submersible platform, with which, numerical simulation was conducted on dynamic response and mooring performance of the platform under extreme sea conditions using Aqwa software by finite element method considering the combined action of random wave, wind and current load with the radiation/diffraction theory, thus obtaining the response amplitude operator (RAO), the additional mass and the rate of radiation damping changing with the wave frequency, as well as the dynamic response of the platform and the tension response of the mooring line under extreme sea conditions. Results show that the platform would have good motion performance with little dynamic response, as the wave frequency is high (above 1.6 rad/s); changes of wave directions have little effect on the heave response; severest responses of surge and pitch would happen in the wave direction of 0 degree; under typhoon sea states, both the peak of dynamic response of platform and the peak of tension response amplitude are above those without typhoon; the dynamic response of platform and the tension response of mooring line increase with worsening sea states.

    floating wind turbine; semi-submersible platform; typhoon sea state; mooring system; dynamic response

    2016-12-05

    2016-12-15

    國家自然科學(xué)基金資助項目(51676131,51176129);上海市科委資助項目(13DZ2260900)

    黃致謙(1993-),男,海南屯昌人,碩士,主要從事漂浮式風(fēng)力機動態(tài)響應(yīng)及結(jié)構(gòu)動力學(xué)響應(yīng)方面的研究.

    李 春(通信作者),男,教授,博導(dǎo),電話(Tel.):13301777337;E-mail:lichunusst@163.com.

    1674-7607(2017)12-1015-08

    TK83

    A

    480.60

    猜你喜歡
    潛式海況浮式
    硫磺上浮式絡(luò)合鐵脫硫工藝緩解溶液起泡研究
    云南化工(2021年5期)2021-12-21 07:41:52
    半潛式平臺表層卡套管處理與認(rèn)識
    海洋石油(2021年3期)2021-11-05 07:43:10
    典型海況下艦載發(fā)射箱結(jié)構(gòu)強度仿真分析
    陵水17-2半潛式生產(chǎn)儲油平臺落物風(fēng)險分析
    關(guān)于浮式防波堤消能效果及透射系數(shù)的研究
    浮式LNG儲存及再氣化裝置(FSRU)淺析及國內(nèi)應(yīng)用推廣展望
    全球首座浮式核電站于今年9月完工
    半潛式鉆井平臺水密滑動門結(jié)構(gòu)疲勞強度評估
    船海工程(2015年4期)2016-01-05 15:53:40
    惡劣海況下海洋石油116內(nèi)轉(zhuǎn)塔式FPSO裝配載優(yōu)化
    船海工程(2015年4期)2016-01-05 15:53:38
    極端海況下軟質(zhì)海崖大規(guī)模蝕退計算模型研究
    亚洲欧美精品综合久久99| 欧美激情在线99| 国产精品亚洲美女久久久| 一边摸一边抽搐一进一小说| 亚洲精品在线观看二区| 国产精品野战在线观看| 精品一区二区三区视频在线 | 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧美人成| 老司机在亚洲福利影院| 精品一区二区三区四区五区乱码| 日本 av在线| 1000部很黄的大片| 看免费av毛片| 欧美黑人巨大hd| 国产精品影院久久| 99热这里只有是精品50| 亚洲九九香蕉| 久久精品国产亚洲av香蕉五月| 不卡av一区二区三区| 午夜激情欧美在线| 亚洲熟妇中文字幕五十中出| 国产精品98久久久久久宅男小说| 国内精品久久久久久久电影| 噜噜噜噜噜久久久久久91| 巨乳人妻的诱惑在线观看| 国产三级中文精品| 成熟少妇高潮喷水视频| 国产一区二区三区在线臀色熟女| 国产高潮美女av| 三级男女做爰猛烈吃奶摸视频| 少妇的逼水好多| 国产高清三级在线| svipshipincom国产片| www国产在线视频色| 欧美3d第一页| 亚洲一区高清亚洲精品| 狂野欧美白嫩少妇大欣赏| 成人无遮挡网站| 国产精品日韩av在线免费观看| 99久国产av精品| 国产69精品久久久久777片 | 国产三级中文精品| 一区福利在线观看| 亚洲熟妇熟女久久| 黄色视频,在线免费观看| 国产蜜桃级精品一区二区三区| 午夜免费激情av| 波多野结衣高清无吗| 午夜两性在线视频| 最新中文字幕久久久久 | 啪啪无遮挡十八禁网站| 一进一出抽搐动态| avwww免费| 啦啦啦韩国在线观看视频| 久久天堂一区二区三区四区| 久久香蕉精品热| 亚洲国产欧美网| 精品99又大又爽又粗少妇毛片 | 日日干狠狠操夜夜爽| 国产成人一区二区三区免费视频网站| 成人欧美大片| 国语自产精品视频在线第100页| 人人妻,人人澡人人爽秒播| 久久精品国产清高在天天线| 禁无遮挡网站| 精品久久久久久久久久免费视频| 不卡av一区二区三区| 欧美+亚洲+日韩+国产| 九色成人免费人妻av| 男插女下体视频免费在线播放| 精品久久蜜臀av无| 成年版毛片免费区| 少妇人妻一区二区三区视频| 欧美激情在线99| 97超级碰碰碰精品色视频在线观看| 日本一本二区三区精品| 国产精品综合久久久久久久免费| 国产黄a三级三级三级人| 给我免费播放毛片高清在线观看| 亚洲欧美日韩高清在线视频| 日本黄色视频三级网站网址| 天堂动漫精品| 国产麻豆成人av免费视频| 男女那种视频在线观看| av在线蜜桃| 宅男免费午夜| 欧美色视频一区免费| 日本一本二区三区精品| 99久久久亚洲精品蜜臀av| 国产成年人精品一区二区| 一边摸一边抽搐一进一小说| 精品一区二区三区av网在线观看| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女| 国产单亲对白刺激| 婷婷精品国产亚洲av| 午夜福利高清视频| 久久国产精品人妻蜜桃| 亚洲专区字幕在线| 黄色片一级片一级黄色片| 美女高潮喷水抽搐中文字幕| 国产精品亚洲一级av第二区| 在线免费观看不下载黄p国产 | 又粗又爽又猛毛片免费看| 亚洲人与动物交配视频| 午夜亚洲福利在线播放| 黄片小视频在线播放| 波多野结衣高清无吗| 亚洲人成网站在线播放欧美日韩| 一区二区三区国产精品乱码| 亚洲欧美日韩高清在线视频| 一本精品99久久精品77| 国内久久婷婷六月综合欲色啪| 欧美成人免费av一区二区三区| 午夜激情欧美在线| 久久午夜亚洲精品久久| 日本五十路高清| 免费观看精品视频网站| 99热精品在线国产| 精品日产1卡2卡| 国产高清videossex| 国产精品电影一区二区三区| 亚洲国产精品合色在线| 欧美乱码精品一区二区三区| 婷婷六月久久综合丁香| 久久精品夜夜夜夜夜久久蜜豆| 天堂动漫精品| 国产精品日韩av在线免费观看| av在线天堂中文字幕| 亚洲欧美日韩卡通动漫| 中亚洲国语对白在线视频| 日本a在线网址| 国产三级黄色录像| 国产精品久久久久久精品电影| 婷婷亚洲欧美| 国产精品99久久久久久久久| 国产精品电影一区二区三区| 国产精品一区二区免费欧美| 亚洲欧美日韩高清在线视频| 亚洲精品乱码久久久v下载方式 | 黄片大片在线免费观看| cao死你这个sao货| 嫩草影院入口| 男女视频在线观看网站免费| 这个男人来自地球电影免费观看| 成人一区二区视频在线观看| 精品久久久久久,| 亚洲精品中文字幕一二三四区| 欧美极品一区二区三区四区| 人妻夜夜爽99麻豆av| 国产一区二区激情短视频| 很黄的视频免费| 国产伦在线观看视频一区| 国产成人系列免费观看| 男人的好看免费观看在线视频| 免费观看精品视频网站| 一级a爱片免费观看的视频| 久久亚洲真实| 国产av麻豆久久久久久久| 此物有八面人人有两片| 免费人成视频x8x8入口观看| 午夜精品一区二区三区免费看| 好男人电影高清在线观看| 99国产综合亚洲精品| 狂野欧美激情性xxxx| 免费无遮挡裸体视频| 一进一出好大好爽视频| 三级毛片av免费| 亚洲专区字幕在线| 日本精品一区二区三区蜜桃| 日本黄色片子视频| 亚洲国产高清在线一区二区三| 天堂√8在线中文| 中文字幕精品亚洲无线码一区| 国产免费男女视频| 亚洲,欧美精品.| 亚洲人与动物交配视频| 亚洲欧美精品综合一区二区三区| 九九热线精品视视频播放| 中出人妻视频一区二区| 国产99白浆流出| 国产精品九九99| 久久久国产成人精品二区| 亚洲av五月六月丁香网| 19禁男女啪啪无遮挡网站| 日本撒尿小便嘘嘘汇集6| 亚洲最大成人中文| 九色成人免费人妻av| 性欧美人与动物交配| 激情在线观看视频在线高清| 哪里可以看免费的av片| 欧美日韩亚洲国产一区二区在线观看| 久久久久久九九精品二区国产| av在线天堂中文字幕| 日本三级黄在线观看| 国产一区二区三区视频了| a级毛片a级免费在线| 欧美另类亚洲清纯唯美| 91字幕亚洲| 亚洲欧美一区二区三区黑人| 欧美一级a爱片免费观看看| 人人妻人人澡欧美一区二区| 俺也久久电影网| 麻豆一二三区av精品| 精华霜和精华液先用哪个| 桃红色精品国产亚洲av| 久久精品亚洲精品国产色婷小说| 亚洲av成人av| 法律面前人人平等表现在哪些方面| 亚洲一区二区三区色噜噜| 国产一区二区三区在线臀色熟女| 国产aⅴ精品一区二区三区波| 国产乱人视频| 免费看十八禁软件| 日韩国内少妇激情av| 久久中文看片网| 美女 人体艺术 gogo| 国产精品国产高清国产av| 深夜精品福利| 欧美成人免费av一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 99国产综合亚洲精品| 亚洲第一电影网av| 日本与韩国留学比较| 午夜福利在线观看吧| 日韩欧美在线乱码| 国产成年人精品一区二区| 久久久久国产一级毛片高清牌| 成人特级黄色片久久久久久久| 精品久久久久久久毛片微露脸| 色尼玛亚洲综合影院| 久久精品人妻少妇| 亚洲五月婷婷丁香| 国产精品一区二区三区四区免费观看 | 久久久水蜜桃国产精品网| 99久国产av精品| 色哟哟哟哟哟哟| 搡老妇女老女人老熟妇| 天堂av国产一区二区熟女人妻| 午夜福利欧美成人| 国产日本99.免费观看| 国产亚洲精品久久久久久毛片| 久久久精品欧美日韩精品| 日本在线视频免费播放| 亚洲国产精品999在线| 男插女下体视频免费在线播放| 国产免费男女视频| 欧美性猛交黑人性爽| www国产在线视频色| 国内精品美女久久久久久| 国产成人精品无人区| 亚洲人成网站高清观看| 一a级毛片在线观看| 在线观看免费午夜福利视频| av天堂中文字幕网| 久久久久久大精品| 人妻丰满熟妇av一区二区三区| 国产真人三级小视频在线观看| 99精品久久久久人妻精品| 国产精品国产高清国产av| 九九热线精品视视频播放| 国产精品久久久久久久电影 | 99热6这里只有精品| 国产精品亚洲一级av第二区| 亚洲av电影在线进入| 国产高清三级在线| 久久99热这里只有精品18| 国产97色在线日韩免费| 亚洲精品在线观看二区| 久久午夜亚洲精品久久| 一二三四在线观看免费中文在| 国产精品九九99| 欧美大码av| 国产成人系列免费观看| 夜夜夜夜夜久久久久| 三级国产精品欧美在线观看 | 国产成+人综合+亚洲专区| 亚洲熟妇中文字幕五十中出| 亚洲国产中文字幕在线视频| 久久久国产欧美日韩av| 国产日本99.免费观看| 免费在线观看影片大全网站| 欧美成狂野欧美在线观看| 国产精品国产高清国产av| 一个人免费在线观看的高清视频| 91老司机精品| 久久婷婷人人爽人人干人人爱| av片东京热男人的天堂| 伊人久久大香线蕉亚洲五| 黄片大片在线免费观看| 亚洲熟女毛片儿| av视频在线观看入口| 欧美日本视频| 亚洲国产欧洲综合997久久,| 精品久久久久久久久久久久久| 婷婷精品国产亚洲av在线| 禁无遮挡网站| 午夜两性在线视频| 久久久国产精品麻豆| 午夜久久久久精精品| 日韩欧美精品v在线| 最近在线观看免费完整版| 国产精品98久久久久久宅男小说| 男女之事视频高清在线观看| 国产亚洲精品一区二区www| 亚洲18禁久久av| 一卡2卡三卡四卡精品乱码亚洲| 男人舔奶头视频| 日本一二三区视频观看| 黄色日韩在线| 国产成人精品无人区| 色综合婷婷激情| 亚洲av免费在线观看| 色综合欧美亚洲国产小说| av中文乱码字幕在线| 精品欧美国产一区二区三| 久久久久国产精品人妻aⅴ院| 国产亚洲欧美在线一区二区| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 欧美一级毛片孕妇| 男人舔奶头视频| 亚洲国产欧美一区二区综合| 老司机在亚洲福利影院| 91av网站免费观看| 久久久久久久久中文| 亚洲熟妇中文字幕五十中出| 日韩欧美国产在线观看| 男女做爰动态图高潮gif福利片| 日本黄色片子视频| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 欧美三级亚洲精品| 国产精华一区二区三区| 国产精品1区2区在线观看.| 免费观看人在逋| 麻豆一二三区av精品| 波多野结衣高清作品| 美女扒开内裤让男人捅视频| 97人妻精品一区二区三区麻豆| 老司机午夜十八禁免费视频| 亚洲精品在线观看二区| 精品熟女少妇八av免费久了| 国产午夜精品论理片| 久久婷婷人人爽人人干人人爱| 国产欧美日韩一区二区三| 午夜免费观看网址| 欧美日韩福利视频一区二区| 欧美绝顶高潮抽搐喷水| 18禁国产床啪视频网站| 午夜激情欧美在线| 熟妇人妻久久中文字幕3abv| 久久精品人妻少妇| 免费看美女性在线毛片视频| 蜜桃久久精品国产亚洲av| 黑人巨大精品欧美一区二区mp4| 波多野结衣高清作品| 麻豆国产av国片精品| 国产亚洲精品综合一区在线观看| 国内揄拍国产精品人妻在线| 免费看十八禁软件| 成人一区二区视频在线观看| 亚洲精品在线观看二区| www.www免费av| 国产亚洲精品综合一区在线观看| 久久精品国产清高在天天线| www国产在线视频色| 人妻夜夜爽99麻豆av| 毛片女人毛片| 悠悠久久av| av女优亚洲男人天堂 | 日本一二三区视频观看| 国内精品美女久久久久久| 精品欧美国产一区二区三| 不卡一级毛片| 久久性视频一级片| 老熟妇仑乱视频hdxx| 亚洲色图 男人天堂 中文字幕| 免费看光身美女| 亚洲18禁久久av| 人妻丰满熟妇av一区二区三区| 一区福利在线观看| 老司机午夜十八禁免费视频| 日韩精品青青久久久久久| 18禁黄网站禁片午夜丰满| 伊人久久大香线蕉亚洲五| 三级毛片av免费| 婷婷精品国产亚洲av在线| 国产视频内射| 亚洲国产精品999在线| 91麻豆av在线| 午夜福利欧美成人| 国内久久婷婷六月综合欲色啪| 丰满的人妻完整版| 男女那种视频在线观看| 欧美黄色片欧美黄色片| 少妇的丰满在线观看| 国产伦精品一区二区三区四那| 精品久久久久久久久久免费视频| 韩国av一区二区三区四区| 黄色女人牲交| 美女高潮的动态| 小说图片视频综合网站| 天天添夜夜摸| 欧美丝袜亚洲另类 | 久久久精品大字幕| 啦啦啦免费观看视频1| 香蕉国产在线看| 黄色成人免费大全| 一个人免费在线观看电影 | 长腿黑丝高跟| 高清在线国产一区| 日韩欧美精品v在线| 18禁黄网站禁片午夜丰满| 国产黄片美女视频| 天天添夜夜摸| 婷婷丁香在线五月| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 国产精品永久免费网站| 精品电影一区二区在线| 这个男人来自地球电影免费观看| 午夜亚洲福利在线播放| 日韩大尺度精品在线看网址| 色噜噜av男人的天堂激情| 三级毛片av免费| av天堂中文字幕网| 午夜亚洲福利在线播放| 一级毛片高清免费大全| 国产私拍福利视频在线观看| 少妇的丰满在线观看| 男人舔奶头视频| tocl精华| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 少妇的逼水好多| 天堂√8在线中文| 制服人妻中文乱码| 男人舔女人下体高潮全视频| 日本黄色片子视频| 国产av在哪里看| 18禁黄网站禁片午夜丰满| 三级国产精品欧美在线观看 | 免费无遮挡裸体视频| 制服丝袜大香蕉在线| 人人妻人人看人人澡| 国产高清视频在线观看网站| 成熟少妇高潮喷水视频| 免费在线观看成人毛片| 91在线精品国自产拍蜜月 | 非洲黑人性xxxx精品又粗又长| 女同久久另类99精品国产91| 久久精品91无色码中文字幕| 老司机在亚洲福利影院| 久久天躁狠狠躁夜夜2o2o| 三级男女做爰猛烈吃奶摸视频| 成人午夜高清在线视频| 亚洲一区二区三区不卡视频| 国产精品日韩av在线免费观看| 国产又黄又爽又无遮挡在线| 男人舔女人下体高潮全视频| 亚洲av电影不卡..在线观看| 午夜影院日韩av| 夜夜夜夜夜久久久久| 欧美一区二区精品小视频在线| 久久久国产成人免费| 999久久久国产精品视频| 国产欧美日韩一区二区精品| 国产高潮美女av| tocl精华| 国产99白浆流出| 又大又爽又粗| 国产精品久久久久久久电影 | 在线免费观看不下载黄p国产 | 亚洲精品国产精品久久久不卡| 久久久久久久久免费视频了| 欧美黑人欧美精品刺激| 国产精品一区二区免费欧美| 亚洲熟妇中文字幕五十中出| 久久久精品大字幕| 久久天堂一区二区三区四区| 国产激情久久老熟女| 99re在线观看精品视频| 两人在一起打扑克的视频| 悠悠久久av| 精品欧美国产一区二区三| 久久久久性生活片| 韩国av一区二区三区四区| 天堂动漫精品| xxx96com| 日本 av在线| 热99re8久久精品国产| 12—13女人毛片做爰片一| 后天国语完整版免费观看| 亚洲精品在线美女| 国产亚洲精品久久久com| 国产视频一区二区在线看| 在线观看午夜福利视频| 一二三四在线观看免费中文在| 亚洲精华国产精华精| 亚洲欧美精品综合一区二区三区| 啪啪无遮挡十八禁网站| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站 | 午夜福利18| av片东京热男人的天堂| 国产男靠女视频免费网站| 欧美日韩福利视频一区二区| ponron亚洲| 一进一出抽搐动态| 又粗又爽又猛毛片免费看| 午夜激情欧美在线| 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 国产野战对白在线观看| 国产麻豆成人av免费视频| 欧美日本亚洲视频在线播放| 久久午夜综合久久蜜桃| 俄罗斯特黄特色一大片| 国产精品久久久久久久电影 | 中文资源天堂在线| 麻豆av在线久日| 国产高潮美女av| 成熟少妇高潮喷水视频| 国产成人福利小说| 国内精品一区二区在线观看| 一本一本综合久久| 国产精品久久久av美女十八| 真实男女啪啪啪动态图| 午夜亚洲福利在线播放| 999精品在线视频| 老司机午夜福利在线观看视频| 免费电影在线观看免费观看| 国产亚洲精品综合一区在线观看| 久久久久国内视频| 婷婷精品国产亚洲av在线| 国产高清videossex| or卡值多少钱| 久久国产精品影院| 两人在一起打扑克的视频| 在线永久观看黄色视频| 真人做人爱边吃奶动态| 亚洲aⅴ乱码一区二区在线播放| 国内精品美女久久久久久| 91麻豆av在线| 天堂√8在线中文| 亚洲午夜精品一区,二区,三区| 欧美色视频一区免费| 啦啦啦免费观看视频1| 天堂网av新在线| 在线观看美女被高潮喷水网站 | 制服人妻中文乱码| 午夜福利免费观看在线| 免费无遮挡裸体视频| 国产成人av教育| 丰满的人妻完整版| 国产午夜精品论理片| 亚洲成av人片免费观看| 亚洲国产精品999在线| 欧美最黄视频在线播放免费| 90打野战视频偷拍视频| 91九色精品人成在线观看| 亚洲国产欧美人成| 一级作爱视频免费观看| 成人av一区二区三区在线看| 99久久精品热视频| 制服人妻中文乱码| 欧美性猛交黑人性爽| 国内精品一区二区在线观看| 在线观看午夜福利视频| 中文亚洲av片在线观看爽| 可以在线观看的亚洲视频| 日本五十路高清| 精品久久久久久成人av| 色综合亚洲欧美另类图片| 国产高清三级在线| 女生性感内裤真人,穿戴方法视频| 一区二区三区国产精品乱码| 国内毛片毛片毛片毛片毛片| 国产免费男女视频| 日韩大尺度精品在线看网址| 国产精品九九99| 国产午夜精品论理片| 亚洲五月天丁香| 99久久久亚洲精品蜜臀av| 日本黄色片子视频| 香蕉av资源在线| 色综合婷婷激情| 欧美在线一区亚洲| 成年女人看的毛片在线观看| 两性夫妻黄色片| 国产野战对白在线观看| 亚洲欧美精品综合久久99| 精品国产美女av久久久久小说| 免费在线观看日本一区| 琪琪午夜伦伦电影理论片6080| 日本熟妇午夜| 亚洲色图 男人天堂 中文字幕| 午夜影院日韩av| 午夜久久久久精精品| 亚洲熟妇熟女久久| 欧美绝顶高潮抽搐喷水| 午夜福利视频1000在线观看| 日韩人妻高清精品专区| 91麻豆av在线| 中文字幕人妻丝袜一区二区| 久久久久久久久中文| 亚洲av第一区精品v没综合| 又紧又爽又黄一区二区| 精品熟女少妇八av免费久了| 美女免费视频网站| 国产在线精品亚洲第一网站| 极品教师在线免费播放| av福利片在线观看|