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

    海上浮式風機在支撐平臺運動影響下的氣動特性研究

    2017-10-11 11:10:27劉格梁胡志強
    海洋工程 2017年1期
    關(guān)鍵詞:葉素浮式葉輪

    劉格梁,胡志強,段 斐

    (上海交通大學 海洋工程國家重點實驗室, 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    海上浮式風機在支撐平臺運動影響下的氣動特性研究

    劉格梁,胡志強,段 斐

    (上海交通大學 海洋工程國家重點實驗室, 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    浮式風機的平臺作為支撐系統(tǒng),其六自由度剛體運動響應(yīng)決定上層風機的運動狀態(tài),也影響葉輪的環(huán)境載荷。葉輪的運動使周圍流場變得復雜,葉片承受著非定常的氣動載荷。研究中考慮葉片結(jié)構(gòu)的運動,將基于定常流場的葉素動量理論用于局部葉素的氣動分析,考慮偏航模型,充分模擬氣流與葉輪的相互作用。在給定平臺各個自由度運動下,計算葉輪氣動參數(shù),分析平臺運動帶來的影響,從入流速度方面探究引起氣動載荷變化的本質(zhì)原因。研究發(fā)現(xiàn),葉輪氣動性能和轉(zhuǎn)動軸的方向有關(guān),若平臺運動改變轉(zhuǎn)動軸方向,會引起氣動載荷以多倍于葉輪旋轉(zhuǎn)的頻率發(fā)生周期變化,平臺轉(zhuǎn)動自由度的運動對葉輪氣動性能影響顯著。

    海上浮式風機;氣動載荷;支撐平臺運動;葉素動量法;非定常分析

    Abstract: This paper addresses a research on the characteristics of aerodynamic loads of Offshore Floating Wind Turbine (OFWT) associated with 6-DOF motions of the supporting platform. The performance of the wind turbine is governed by 6-DOF motion of the platform, which is deemed as the supporting base of an OFWT. Motions of the rotor complicate the surrounding flow field resulting in unsteady aerodynamic loads on the blades. A quasi-steady method is presented and the motion of the blades is combined with the Blade-Element-Momentum method to analyze the local blade load. A skewed model is introduced to simulate the interaction between wind air flow and the rotor. The calculation of aerodynamic load under prescribed single-DOF motion of the platform shows the effect of the motion on the rotor. Essential reasons for changing of aerodynamic loads are investigated from the view of inflow velocity. Results show that the aerodynamic performance of rotor is related to the orientation of rotational axis. If the orientation of the rotation axis is altered by motions of the platform, it will cause aerodynamic loads to change periodically at a frequency, which is times more than that of the rotor rotation. Thus, platform′s rotational DOFs have a significant impact on the aerodynamic behavior of the rotor.

    Keywords: offshore floating wind turbine; aerodynamic loads; supporting platform′s motions; blade-element-momentum method; unsteady analysis

    風機用于將空氣動能轉(zhuǎn)化為電能,為提供能源。一般風機的塔架結(jié)構(gòu)都固定在大地或者海底,但是隨著水深的增加,固定式支撐結(jié)構(gòu)的設(shè)計難度以及成本都急劇增加,浮式平臺支撐的浮式風機成為未來深水風機發(fā)展的必然趨勢。和固定式風機不同,海上浮式風機的平臺在風浪流的作用下,將發(fā)生六自由度運動,風機的氣動載荷加劇了平臺運動響應(yīng);同時,平臺的運動導致了葉輪姿態(tài)和速度的變化,引起葉輪周圍流場的非定常分布,葉輪氣動載荷從而發(fā)生相應(yīng)變化。葉輪不斷地旋轉(zhuǎn),其氣動載荷一般會隨著旋轉(zhuǎn)位置而變化,平臺運動會影響到這種周期性變化的載荷。因此,平臺運動與葉輪載荷是相互耦合的過程。

    這一耦合問題為浮式風機的設(shè)計與分析帶來了挑戰(zhàn),近年來也深受研究人員們關(guān)注。Sebastian等基于自由尾渦方法(FVM)分析了平臺縱搖運動下,葉輪與下風向的尾渦相互作用,認為由平臺縱搖運動引起的葉輪氣動載荷的變化是不可忽視的[1]。Sebastian等認為平臺運動引起了非定常的流場,對比了固定式風機與浮式風機的氣動載荷,研究了不同平臺形式下,對氣動影響較大的平臺運動自由度[2]。Jeon等采用渦格法(VLM)同樣分析了縱搖運動的影響[3]。Farrugia等開展浮式平臺對氣動性能的試驗研究,測量了葉輪的載荷和尾流狀態(tài),提出浮式風機應(yīng)減小縱蕩運動,以避免在葉尖速比較高狀態(tài)下運行[4]。Rockel等通過風洞試驗研究了平臺縱搖運動對尾流的影響,浮體其他自由度的運動也可能使尾流更復雜[5]。Bayati等認識到研究浮體對氣動影響的必要性,設(shè)計了在風洞試驗中能為平臺提供給定單/多自由度的機械裝置[6]。Tran等應(yīng)用CFD方法研究了周期縱搖運動對浮式風機氣動的影響,并對比了其他計算模型,在大幅值縱搖運動的模擬中,CFD方法比傳統(tǒng)氣動模型計算的載荷小,且能更精確計算非定常載荷和平臺隨機運動的影響[7]。Sant等根據(jù)多種模型研究了平臺在不同首搖角度下的葉輪性質(zhì)[8]。Wu等利用CFD方法研究了平臺的縱蕩、縱搖和首搖對整體葉輪、單根葉片的氣動載荷的影響,結(jié)合盤面和尾流的速度分布,對平臺運動帶來的葉輪功率、載荷變化展開分析,認為這些變化將是研究風機控制系統(tǒng)、結(jié)構(gòu)疲勞和強度所面臨的重要挑戰(zhàn)[9]。但是,這些研究僅針對平臺部分自由度運動,缺乏對六自由度運動的全面分析和對比。本文著眼于平臺六自由度運動,從偏航、速度分布等方面揭示了氣動載荷變化的本質(zhì)原因。

    本文研究海上浮式風機支撐平臺運動對葉輪氣動載荷的影響。先根據(jù)運動轉(zhuǎn)化關(guān)系,得到葉素在平臺運動影響下的運動,疊加環(huán)境風速,考慮葉輪附近相對來流的非定常性。再將受葉輪影響的空氣流域分割成局部葉素影響的扇環(huán)流管,應(yīng)用經(jīng)典葉素動量理論分析該區(qū)域,并引入偏航模型研究來流與盤面的夾角對流場的影響。利用非定常氣動模型,分析在給定平臺各自由度運動情況下葉輪氣動載荷的變化規(guī)律。

    1 平臺和葉輪的運動關(guān)系

    海上浮式風機的支撐系統(tǒng)除了塔架,最重要的是支撐浮式平臺。平臺受到水動力載荷、系泊載荷,以及由塔架傳遞的上層結(jié)構(gòu)的環(huán)境載荷和慣性載荷,其運動響應(yīng)復雜。為了方便描述平臺的運動,首先定義分析中所使用的坐標系,并介紹平臺和葉輪運動關(guān)系。如圖 1所示,定義固結(jié)于靜水面的大地坐標系ox0y0z0,平行于大地坐標系且與平臺固結(jié)的參考坐標系ox1y1z1,固結(jié)于平臺并隨之轉(zhuǎn)動的隨體坐標系oxyz,上層結(jié)構(gòu)都在隨體坐標系中表達[10-11]。平臺參考坐標系原點相對于大地坐標系原點的矢徑表示平臺的三個平動位移,分別是沿x0方向的縱蕩、沿y0方向的橫蕩、沿z0方向的垂蕩,文中分別使用自由度編號1-3表示;平臺隨體坐標系相對于參考坐標系坐標軸的轉(zhuǎn)動量是平臺的轉(zhuǎn)動位移,即繞x軸轉(zhuǎn)動的橫搖α、繞y軸轉(zhuǎn)動的縱搖β、繞z軸轉(zhuǎn)動的首搖γ,文中分別使用自由度編號4-6表示,采用卡爾丹角描述。用方向余弦陣表征坐標系之間的轉(zhuǎn)化關(guān)系,方向余弦陣為:

    平臺的轉(zhuǎn)動角速度為:

    固結(jié)于平臺隨體系的某質(zhì)點的矢徑為:

    質(zhì)點的速度:

    一般認為機艙和塔架剛性固定于平臺上,安裝在風機頂部的葉輪隨平臺運動,同時還繞轉(zhuǎn)動軸旋轉(zhuǎn)。葉輪上各葉素的速度表示為:

    其中,rhub為輪轂中心在平臺隨體系中的矢徑,rBE為葉素相對輪轂中心的矢徑,ωR為葉輪旋轉(zhuǎn)角速度。上式右端的前兩項分別是平臺平動和轉(zhuǎn)動運動對葉素速度的貢獻。

    圖1 平臺坐標系以及風機運動的示意Fig. 1 Sketch of platform coordinate system and motions of the wind turbine

    2 氣動載荷計算模型

    2.1經(jīng)典葉素動量法

    葉素動量法是計算風機氣動載荷的經(jīng)典方法[12],是一維動量理論與二維翼型理論的結(jié)合。動量理論假定葉片數(shù)無窮,葉輪沿徑向的流管是均勻的,流管之間相互獨立互不干擾,盤面推力均勻。圖 2中所示的是流管內(nèi)氣流速度和壓強的變化關(guān)系。根據(jù)盤前盤后的運動關(guān)系,利用連續(xù)性方程和伯努利方程計算葉輪作用在流管上的力。翼型理論假設(shè)翼型展向無限長,分析流經(jīng)翼型的二維流場,計算氣動載荷。圖 3所示,vrel為相對入流速度,α為決定載荷的攻角。假定通過上述兩種方法得到的載荷相等,即可實現(xiàn)關(guān)鍵參數(shù)的求解計算。葉素動量法中最為關(guān)鍵的是誘導速度,表征了葉輪對氣流的影響效果,通常用遠前方的來流速度及軸向誘導因子a表示,用周向誘導因子a′考慮葉輪的旋轉(zhuǎn)效應(yīng)。根據(jù)一維和二維理論在動力上的等效關(guān)系,確定誘導速度和盤面附近空氣的流動狀態(tài),計算氣動載荷。

    圖2 一維動量理論Fig. 2 1-D momentum theorem

    圖3 二維翼型理論Fig. 3 2-D airfoil theory

    經(jīng)典葉素動量理論的推力系數(shù)表達如下:

    式(6)對應(yīng)圖2中的動量理論,其中F為普朗特修正因子,為葛勞沃特臨界誘導因子。式(7)對應(yīng)圖3中的翼型理論,σ為葉素實度,Cn為葉素軸向力系數(shù)。

    2.2非定常氣動模型

    經(jīng)典葉素動量法中,假定葉輪固定,來流到入流速度都是在盤面均勻定常分布的。而實際上,由于支撐系統(tǒng)的運動或變形會直接影響葉輪的速度、時變的風速等,任意時刻,盤面不同位置的相對來流速度是不一樣的,誘導速度也呈現(xiàn)出隨時間和空間的變化分布,很多模型都研究這類非定常氣動特性[13]。本文提出使用局部葉素分析法,考慮將上述定常方法擴展至非定常氣動計算,并進行局部分析。

    局部葉素分析方法中,氣動控制體由圓環(huán)柱狀流管,變?yōu)閮H受單根葉片影響的扇環(huán)流管,前述的經(jīng)典葉素動量理論可以方便適用,圖4給出的是三葉片風機扇環(huán)流管的示意圖。

    圖4 受局部葉素影響的扇環(huán)流管示意Fig. 4 Sketch for a sector ring tube subjected to local blade element

    通過運動學關(guān)系得到局部葉素的位置和速度vBE,根據(jù)環(huán)境條件計算該葉素位置對應(yīng)的來流風速vwind,疊加得到遠前方的相對遠前方來流速度v,假定誘導因子a和a′,得到在葉素上的相對入流速度vrel:

    此時得到的氣流速度應(yīng)用于翼型理論,確定局部葉素的誘導速度并計算氣動載荷。

    圖5 葉輪偏航位置示意Fig. 5 Yawed rotor sketch

    2.3偏航模型

    在引入了非定常的局部葉素氣動計算方法后,需要考慮局部偏航,葉輪軸向與入流速度方向不共線的位置稱為偏航位置。相對遠前方來流速度與葉片的局部扇形盤面非正交時,葉素處于偏航狀態(tài),引起誘導因子在空間上分布不均勻。圖5所示,氣流流過葉輪,尾流會傾斜,同樣導致誘導速度隨方位角而變化,葉片指向上風方向比葉片轉(zhuǎn)動半圈后指向下風方向的誘導速度要小。偏航狀態(tài)和相對遠前方來流速度密切相關(guān),考慮針對局部葉素的簡單偏航修正模型[14]:

    其中,γ為偏航角度,是來風速度與葉輪軸向的夾角;χ為尾流傾斜與軸向的夾角,是氣流經(jīng)過葉輪后與葉輪軸向的夾角;θwing為葉片當前的方位角,θ0葉片指向尾流最深處時的方位角,此時葉片所在的直徑為偏航線。圖5中,入流速度在xy平面內(nèi),y軸為偏航線,y軸負方向為θ0的位置。

    2.4浮式風機氣動載荷計算流程

    平臺的位移決定了風輪的位置和姿態(tài),為了準確分析浮式風機的氣動特性,在使用非定常氣動模型時,計入平臺運動的影響,并分析考察其中的動態(tài)規(guī)律。根據(jù)平臺的運動、風機幾何尺寸以及葉輪轉(zhuǎn)動的速度,求出局部葉素相應(yīng)的位置、速度以及軸向方向,并用于氣動載荷的計算。葉輪氣動載荷的計算過程如下:

    1) 根據(jù)式(3)和式(4),利用平臺六自由度位移量及其導數(shù),計算平臺剛體運動速度;

    2) 根據(jù)結(jié)構(gòu)幾何參數(shù),計算葉素的空間位置,根據(jù)式(5)計算葉素的絕對速度;

    3) 根據(jù)風參數(shù),計算葉素的遠前方來流速度,疊加葉素自身速度,由式(8)得到相對遠前方來流速度;

    4) 假設(shè)初始誘導因子,計算攻角,根據(jù)翼型得到動力系數(shù),求解平衡方程,并迭代至收斂;

    5) 根據(jù)葉素方向、相對遠前方來流速度,計算偏航角等,利用偏航模型計算修正后的誘導因子;

    6) 計算速度分布,確定該葉素的軸向、切向和轉(zhuǎn)矩載荷,計算葉素對轉(zhuǎn)軸的局部扭矩;

    7) 根據(jù)以上步驟,計算所有葉素,得到整個葉輪的氣動載荷。

    本研究中,利用Fortran 2003編寫了基于上述理論的計算程序,用于分析浮式風機的性能。

    3 氣動特性分析

    以NREL 5MW上風型水平軸風機[15]為研究對象,外界環(huán)境為均勻定常風,考察在額定風速11.4 m/s、額定轉(zhuǎn)速12.1 r/min情況下的氣動載荷。首先研究了平臺固定時風機自身的載荷特點,然后分別給定平臺各自由度的位移和運動,通過載荷結(jié)果評價平臺運動對氣動性能的影響。

    3.1平臺固定時葉輪的氣動特點

    由于安裝和安全等方面的考慮,通常會設(shè)置葉輪轉(zhuǎn)動軸的傾斜角和葉片結(jié)構(gòu)的錐角,文中稱此類因結(jié)構(gòu)布置而出現(xiàn)的轉(zhuǎn)動角為結(jié)構(gòu)角。圖6中,NREL 5MW風機的葉輪設(shè)置了5°的上傾角和2.5°的錐角。

    圖6 NREL-5MW風機葉輪的結(jié)構(gòu)角示意Fig. 6 Sketch for structural angles of NREL-5MW rotor

    在計算平臺運動的影響之前,首先研究了在風機的平臺固定時,葉輪結(jié)構(gòu)角帶來的影響。結(jié)果表明,如果不存在任何結(jié)構(gòu)角或者僅有葉片錐角,風機的功率和載荷等參數(shù)都是定常值,而對于設(shè)置了轉(zhuǎn)動軸傾角的葉輪,其氣動參數(shù)發(fā)生周期變化的。表1列出了平臺固定時,風機的氣動輸入功率、軸向推力統(tǒng)計值。從表中可知,存在結(jié)構(gòu)角的葉輪結(jié)構(gòu),比沒有結(jié)構(gòu)角對氣體流動的抑制作用小,轉(zhuǎn)化的風能以及受到的氣動載荷也相應(yīng)減小。在結(jié)構(gòu)角中,轉(zhuǎn)動軸傾角引起葉輪偏航,是影響載荷的主要原因。

    表1 結(jié)構(gòu)角對風機功率和軸向推力的影響Tab. 1 Effect of structural angles on power and thrust

    圖7~圖9給出了在平臺固定時,僅有轉(zhuǎn)動軸傾角的葉輪的氣動功率、軸向推力和單根葉片的時歷。傾角使得葉輪處于定常角偏航狀態(tài),單根葉片載荷呈1P變化關(guān)系,3葉片的葉輪整體載荷發(fā)生3P變化。從時歷的結(jié)果可知,尤其是從軸向推力上,葉輪載荷關(guān)于某位置的對稱性。這是由于偏航計算時,來流速度在葉片隨體坐標系翼型剖面內(nèi)的速度分量不變,而翼型剖面的軸向速度,即風速在葉片展向的速度分量對載荷的影響,是通過偏航修正得到的。從式(11)分析可得,葉素的誘導因子是關(guān)于偏航線對稱的,故載荷關(guān)于偏航線對稱。當同時存在轉(zhuǎn)動軸傾角和葉片錐角時,如圖8所示,此時翼型剖面內(nèi)的速度在不同葉片上分量不等,仍然有3P響應(yīng),但誘導因子不再關(guān)于偏航線對稱。葉片錐角的存在減少了捕獲的空氣能量,而葉輪轉(zhuǎn)動軸的傾斜會引起偏航,因此結(jié)構(gòu)角中轉(zhuǎn)軸傾角決定了載荷變化的特點??紤]到結(jié)構(gòu)角對氣動載荷的擾動特點,在研究平臺運動影響的過程中,不設(shè)置葉輪的結(jié)構(gòu)角,以避免出現(xiàn)混淆。

    圖7 僅有轉(zhuǎn)動軸傾角的葉輪氣動載荷時歷Fig. 7 Time history for aerodynamic loads of rotor with shaft tilt only

    圖8 同時存在傾角和錐角的葉輪氣動載荷時歷Fig. 8 Aerodynamic loads of rotor with shaft tilt and cone angle

    圖9 僅存在轉(zhuǎn)動軸傾角時,單根葉片所受推力分布隨時間的變化情況Fig. 9 Time-varying thrust distribution along a single blade only with shaft tilt

    3.2平臺各自由度對氣動載荷的影響分析

    為了探究平臺運動對氣動載荷的影響,根據(jù)表中的計算情況,為平臺各自由度預(yù)設(shè)了位移和運動,對無結(jié)構(gòu)角的葉輪結(jié)構(gòu)進行了計算與分析,參數(shù)見表2。計算編號表示為靜態(tài)定常位移或者正弦運動(S/D)-六自由度編號(1-6)-幅值編號(1-6)-頻率編號(1-6),自由度1-3為線位移,幅值單位為m,自由度4-6為角位移,幅值單位為°,定常位移末位編號為0。如D532表示計算的是無結(jié)構(gòu)角葉輪在平臺縱搖自由度發(fā)生幅值為8°、頻率為0.15 rad/s正弦運動對應(yīng)的氣動情況。

    表2 運動自由度的計算參數(shù)編號Tab. 2 Numbers for computing parameters of motion DOF

    針對各自由度位移為定常值的情況。當平臺線位移為定常值,或發(fā)生常值的橫搖角位移時,葉輪的軸向方向不變,空氣入流速度和葉素的速度計算都未受到平臺影響,葉輪的狀態(tài)與表1中平臺固定無結(jié)構(gòu)角的情況相同,氣動載荷不變。當平臺縱搖和首搖位移固定時,葉輪平面的軸向發(fā)生了變化,與遠前方來流不正交,葉輪處于定常偏航狀態(tài)。圖10給出了S530葉輪的氣動輸入功率和軸向推力的時歷,定常首搖或縱搖同葉輪僅存在轉(zhuǎn)動軸傾角的載荷規(guī)律是一致的,都存在3P變化規(guī)律,而且關(guān)于偏航線對稱的特點更為明顯,都是定常偏航的結(jié)果。圖11的是氣動載荷隨著縱搖角度的變化趨勢,功率和軸向推力都隨著角度增大而減小。引起葉輪偏航的縱搖和首搖既限制了氣動輸入功率和載荷,又引起載荷高頻振蕩,可見固定縱搖和首搖位移同葉輪轉(zhuǎn)軸傾斜引起的偏航現(xiàn)象一樣,對氣動載荷的影響規(guī)律一致。

    圖10 S530(縱搖)的功率和軸向推力的時歷Fig. 10 Time history for power and thrust of S530(pitch)

    圖11 功率和軸向推力隨縱搖角的變化Fig. 11 Change of power and thrust with pitch angle

    平臺各自由度發(fā)生的正弦運動對氣動載荷的影響不盡相同。從結(jié)果中發(fā)現(xiàn),縱蕩、橫搖、縱搖自由度發(fā)生正弦運動時,氣動載荷主要由平臺運動決定,變化周期與平臺運動周期相同。從圖12可知,D133的氣動結(jié)果,縱蕩運動直接改變了相對風速,而葉輪的軸向不變,各葉素處的氣動載荷受縱蕩速度影響。橫搖和縱搖運動影響下的結(jié)果基本一致,圖13中是考慮僅有縱搖自由度運動的結(jié)果。平臺發(fā)生橫搖和縱搖運動,同僅有轉(zhuǎn)動定常角位移帶來的結(jié)果是不一樣的。轉(zhuǎn)動固定角度之后,平臺沒有速度,不影響葉輪上各處速度的分布,相對來流仍然是均勻定常風,此時定常偏航角度對氣動影響明顯。而平臺搖動時,平臺角速度引起各葉素的相對來流速度在空間上分布不均勻,且隨時間不斷變化,需要根據(jù)轉(zhuǎn)子平面轉(zhuǎn)動瞬間的角速度確定葉輪平面各點的速度,與風速疊加后,得到的相對來流隨著平臺的位置和速度而變化。圖14中是當僅有縱搖情況下風機縱剖面的速度分布示意圖。從圖中可知,由于浮式風機系統(tǒng)縱剖面的轉(zhuǎn)動瞬心一般在水線面以下,上層結(jié)構(gòu)隨著平臺發(fā)生的轉(zhuǎn)動時,葉素的速度隨著高度增加而增大。從式(5)分析,由于葉輪至轉(zhuǎn)動瞬心的距離和葉片長度較大,平臺轉(zhuǎn)動引起的葉素速度變化遠大于葉輪轉(zhuǎn)動的影響,盤面上葉素速度的跨度較大。因此與平臺僅有縱蕩運動時的結(jié)果相比,在縱搖運動下,氣動載荷的振蕩幅度更大。

    圖12 D133(縱蕩)的功率、軸向推力時歷Fig. 12 Time history for power and thrust of D133(surge)

    圖13 D533(縱搖)的功率、軸向推力時歷Fig. 13 Time history for power and thrust of D533(pitch)

    相對于上述三個方向的運動,平臺發(fā)生橫蕩、垂蕩和首搖運動時,葉輪的旋轉(zhuǎn)會造成載荷發(fā)生周期振蕩。圖15和圖16是平臺首搖運動的結(jié)果,葉輪關(guān)于首搖軸對稱,葉素的牽連速度線性分布,且葉輪的兩個半面交替進入入流和尾流中,氣動載荷會發(fā)生二倍于首搖頻率的變化。此時葉素的速度以及偏航角度受轉(zhuǎn)動的影響很明顯,載荷3P響應(yīng)的振蕩幅值和平臺轉(zhuǎn)動引起的振蕩幅值比較大,不利于葉輪功率穩(wěn)定輸出。橫蕩和垂蕩運動結(jié)果是一致的,圖17中是平臺僅發(fā)生橫蕩的氣動結(jié)果,在風速和平臺運動速度疊加之后,相對來流等效為偏航角正弦變化的偏航風速,其周期為平臺運動周期的一半。此時葉輪旋轉(zhuǎn)引起載荷的3P振蕩,但與首搖相比,振蕩幅值比較小,對葉輪的影響可忽略不計。

    圖14 縱搖運動時,葉輪縱剖面的速度分布Fig. 14 Velocity distribution in the longitudinal section of rotor when pitching

    圖15 D633(首搖)的的功率時歷Fig. 15 Time history for power of D633(yaw)

    圖16 D633(首搖)的軸向推力和氣動首搖力矩時歷Fig. 16 Time history for thrust and aero yaw moment of D633(yaw)

    圖17 D233(橫蕩)的功率、軸向推力時歷Fig. 17 Time history for power and thrust of D233(sway)

    4 結(jié) 語

    基于葉素動量理論,充分考慮葉輪周圍流場的非定常性,計算風機的氣動載荷。通過對葉輪載荷的分析,以及對平臺運動與氣動載荷之間對應(yīng)關(guān)系的分析,探究平臺六自由度運動對于氣動載荷的影響程度,并得到以下結(jié)論:

    1) 葉輪軸向方向是影響葉輪氣動性能的主要因素,與氣流完全正交的葉輪具有吸收最大風能的能力,氣動載荷也最大;葉輪軸向方向變化會引起偏航,功率和載荷都減小。

    2) 平臺定常的縱搖或首搖位移會引起偏航,氣動載荷減小,并發(fā)生3P振蕩,其他自由度的定常位移影響較小。

    3) 平臺發(fā)生縱蕩、橫搖和縱搖運動時,非定常相對入流的速度主要受平臺運動的影響。

    4) 平臺首搖運動導致氣動載荷發(fā)生3P振蕩,其振幅較大,影響不可忽略。

    因此,為了使浮式風機獲得更大的風能利用效率,應(yīng)增大平臺在縱搖和首搖方向上的剛度,盡可能減小因平臺轉(zhuǎn)動引起的偏航角,確保葉輪能夠正對著來流方向。為了保證浮式風機能夠安全可靠地工作,應(yīng)努力限制平臺縱蕩和三個轉(zhuǎn)動自由度的運動。

    [1] SEBASTIAN T, LACKNER M A. Analysis of the induction and wake evolution of an offshore floating wind turbine[J]. Energies, 2012, 5(4): 968-1000.

    [2] SEBASTIAN T, LACKNER M A. Characterization of the unsteady aerodynamics of offshore floating wind turbines [J]. Wind Energy, 2013, 16(3): 339-352.

    [3] JEON M, LEE S, LEE S. Unsteady aerodynamics of offshore floating wind turbines in platform pitching motion using vortex lattice method [J]. Renewable Energy, 2014, 65: 207-212.

    [4] FARRUGIA R, SANT T, MICALLEF D. Investigating the aerodynamic performance of a model offshore floating wind turbine [J]. Renewable Energy, 2014, 70: 24-30.

    [5] ROCKEL S, CAMP E, SCHMIDT J, et al. Experimental study on influence of pitch motion on the wake of a floating wind turbine model [J]. Energies, 2014, 7(4): 1954-1985.

    [6] BAYATI I, BELLOLI M, FERRARI D, et al. Design of a 6-DOF robotic platform for wind tunnel tests of floating wind turbines [J]. Energy Procedia, 2014, 53: 313-323.

    [7] TRAN T, KIM D, SONG J. Computational fluid dynamic analysis of a floating offshore wind turbine experiencing platform pitching motion [J]. Energies, 2014, 7(8): 5011-5026.

    [8] SANT T, CUSCHIERI K. Numerical modelling of the aerodynamic characteristics of a floating offshore wind turbine under yawed rotor conditions[C]//Proceedings of the 34th International Conference on Ocean, Offshore and Arctic Engineering. 2015: OMAE 2015-41604.

    [9] WU J, DING J H, HE Y P, et al. Study on unsteady aerodynamic performance of floating offshore wind turbine by CFD method[C]//Proceedings of the 25th International Ocean and Polar Engineering Conference. 2015: ISOPE 2015-DCW-13.

    [10] 劉延柱,洪嘉振,楊海興.多剛體系統(tǒng)動力學[M].北京:高等教育出版社, 1989. (LIU Y Z, HONG J Z, YANG H X. Dynamics for multi-rigid-body system[M]. Beijing: Higher Education Press, 1989. (in Chinese))

    [11] NEWMAN J N. Marine hydrodynamics[M]. MIT Press, 1977.

    [12] HANSEN M O L. Aerodynamics of wind turbines[J]. Rotors Loads and Structure James, 2015, 5(2-3):141-167.

    [13] LEISHMAN J G. Challenges in modelling the unsteady aerodynamics of wind turbines [J]. Wind Energy, 2002, 5(2-3): 85-132.

    [14] BURTON T, SHARPE D, JENKINS N, et al. Wind energy handbook [M]. John Wiley & Sons, 2001.

    [15] JONKMAN J M, BUTTERFIELD S, MUSIAL W, et al. Definition of a 5 MW reference wind turbine for offshore system development[J]. National Renewable Energy Laboratory Tech. Rep., 2009: 10.2172/947422.

    Research about influence of supporting platform′s motions on aerodynamics of offshore floating wind turbine

    LIU Geliang, HU Zhiqiang, DUAN Fei

    (State Key Laboratory of Ocean Engineering, Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai Jiao Tong University, Shanghai 200240, China)

    TK83;P751

    A

    10.16483/j.issn.1005-9865.2017.01.005

    1005-9865(2017)01-0042-09

    2016-01-09

    國家自然科學基金資助項目(51239007)

    劉格梁(1991-),男,江蘇靖江人,碩士研究生,研究海上浮式風機。E-mail: rascalpirate@sjtu.edu.cn

    胡志強(1975-),男,博士,副教授。E-mail: zhqhu@sjtu.edu.cn

    猜你喜歡
    葉素浮式葉輪
    硫磺上浮式絡(luò)合鐵脫硫工藝緩解溶液起泡研究
    云南化工(2021年5期)2021-12-21 07:41:52
    1.4317 QT2鋼在高能泵葉輪上的應(yīng)用
    關(guān)于浮式防波堤消能效果及透射系數(shù)的研究
    浮式LNG儲存及再氣化裝置(FSRU)淺析及國內(nèi)應(yīng)用推廣展望
    應(yīng)用石膏型快速精密鑄造技術(shù)制造葉輪
    不能說的秘密
    晚報文萃(2016年7期)2017-01-10 13:11:54
    精心設(shè)計老婆出軌釀滅門慘案
    新傳奇(2016年29期)2016-09-18 08:30:12
    那個不能說的秘密
    中外文摘(2016年12期)2016-08-29 08:14:40
    離心泵葉輪切割方法
    全球首座浮式核電站于今年9月完工
    在线看a的网站| 丁香六月天网| 熟妇人妻不卡中文字幕| 纯流量卡能插随身wifi吗| 色婷婷久久久亚洲欧美| 欧美人与性动交α欧美软件| 超碰成人久久| 亚洲精品av麻豆狂野| 高清不卡的av网站| 亚洲熟女毛片儿| 精品一品国产午夜福利视频| 久久久亚洲精品成人影院| 亚洲精品乱久久久久久| 一本久久精品| 久久精品国产a三级三级三级| 欧美黄色片欧美黄色片| 国产精品成人在线| av一本久久久久| 老司机亚洲免费影院| av在线app专区| 女性生殖器流出的白浆| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久精品精品| 18禁观看日本| 午夜福利免费观看在线| 色婷婷av一区二区三区视频| 亚洲精品国产av成人精品| 日韩制服丝袜自拍偷拍| 色94色欧美一区二区| 日本色播在线视频| 日日摸夜夜添夜夜爱| 亚洲精品国产区一区二| 精品少妇黑人巨大在线播放| 精品国产乱码久久久久久小说| 精品一区二区三区av网在线观看 | 99国产精品免费福利视频| 1024视频免费在线观看| 国产一区二区 视频在线| 国产精品久久久久久精品电影小说| 亚洲成人手机| 亚洲av日韩在线播放| 高清黄色对白视频在线免费看| av网站免费在线观看视频| 免费观看人在逋| 成人亚洲欧美一区二区av| 午夜影院在线不卡| videos熟女内射| 美女午夜性视频免费| 国产午夜精品一二区理论片| 人人妻人人爽人人添夜夜欢视频| 久久久久精品人妻al黑| 国产精品.久久久| 香蕉丝袜av| 夫妻性生交免费视频一级片| 99精国产麻豆久久婷婷| 18在线观看网站| 又大又爽又粗| 免费观看a级毛片全部| 国产片内射在线| 99热全是精品| 九草在线视频观看| 七月丁香在线播放| 五月天丁香电影| 国产色婷婷99| 视频在线观看一区二区三区| 久久ye,这里只有精品| 最近的中文字幕免费完整| 久久久欧美国产精品| 十八禁人妻一区二区| 成人影院久久| 亚洲精品成人av观看孕妇| 国产不卡av网站在线观看| 亚洲av国产av综合av卡| 伦理电影免费视频| 亚洲精品第二区| 两个人免费观看高清视频| 久久免费观看电影| 久久久久久久久免费视频了| 自线自在国产av| 成人免费观看视频高清| 少妇的丰满在线观看| 国产精品久久久久成人av| 丝袜喷水一区| 久久天躁狠狠躁夜夜2o2o | 午夜日韩欧美国产| 天天躁狠狠躁夜夜躁狠狠躁| 男女边摸边吃奶| 青春草亚洲视频在线观看| 国产精品国产三级国产专区5o| 18禁国产床啪视频网站| 成人国产麻豆网| 高清不卡的av网站| 午夜福利视频在线观看免费| a级毛片在线看网站| 中文字幕最新亚洲高清| 亚洲精品国产av蜜桃| 欧美激情极品国产一区二区三区| 日韩大片免费观看网站| 97精品久久久久久久久久精品| 另类亚洲欧美激情| 国产精品成人在线| 高清不卡的av网站| 男人添女人高潮全过程视频| 街头女战士在线观看网站| 天堂8中文在线网| 国产精品久久久久久久久免| 国产熟女欧美一区二区| 丰满乱子伦码专区| 国产精品久久久久久久久免| 国产成人精品久久久久久| 国产男人的电影天堂91| 热re99久久国产66热| 亚洲欧美激情在线| 国产伦人伦偷精品视频| 亚洲一区二区三区欧美精品| 国产 一区精品| 午夜av观看不卡| 成人三级做爰电影| av网站在线播放免费| 久久热在线av| netflix在线观看网站| 精品少妇一区二区三区视频日本电影 | 两性夫妻黄色片| 亚洲av福利一区| 卡戴珊不雅视频在线播放| 亚洲欧美日韩另类电影网站| 国产又爽黄色视频| 亚洲精品视频女| 激情五月婷婷亚洲| 超碰97精品在线观看| 国产一区二区激情短视频 | 欧美激情高清一区二区三区 | 亚洲国产看品久久| 免费高清在线观看视频在线观看| 丝袜美足系列| 成人黄色视频免费在线看| 天美传媒精品一区二区| 欧美成人精品欧美一级黄| 欧美激情极品国产一区二区三区| 亚洲av福利一区| 亚洲欧美激情在线| 老司机影院毛片| 亚洲国产欧美在线一区| 777久久人妻少妇嫩草av网站| 午夜免费男女啪啪视频观看| 人人妻人人澡人人看| 亚洲精品国产av成人精品| avwww免费| 国产精品国产三级专区第一集| 亚洲精品久久午夜乱码| 久久国产精品男人的天堂亚洲| 色精品久久人妻99蜜桃| 精品少妇内射三级| 亚洲人成网站在线观看播放| 菩萨蛮人人尽说江南好唐韦庄| 中文乱码字字幕精品一区二区三区| a级片在线免费高清观看视频| 香蕉国产在线看| 美女福利国产在线| 国产色婷婷99| 午夜老司机福利片| 在线免费观看不下载黄p国产| 国产探花极品一区二区| 亚洲av男天堂| 激情视频va一区二区三区| 精品酒店卫生间| 成人国产av品久久久| 免费高清在线观看日韩| 成年美女黄网站色视频大全免费| 日日啪夜夜爽| 在线观看免费视频网站a站| 宅男免费午夜| 在线看a的网站| 99热网站在线观看| 两个人看的免费小视频| 一本久久精品| 在线观看免费日韩欧美大片| 午夜福利网站1000一区二区三区| 波多野结衣一区麻豆| 一本色道久久久久久精品综合| 又黄又粗又硬又大视频| 巨乳人妻的诱惑在线观看| 午夜福利在线免费观看网站| 午夜av观看不卡| 男女国产视频网站| 在现免费观看毛片| 亚洲少妇的诱惑av| 国产精品二区激情视频| 永久免费av网站大全| 日本vs欧美在线观看视频| 卡戴珊不雅视频在线播放| 18禁观看日本| 欧美日韩成人在线一区二区| 日韩欧美精品免费久久| a级片在线免费高清观看视频| 精品国产一区二区久久| 亚洲欧美激情在线| 精品国产乱码久久久久久男人| 免费av中文字幕在线| 老司机靠b影院| 国产有黄有色有爽视频| 国精品久久久久久国模美| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| 丰满少妇做爰视频| 国产不卡av网站在线观看| 国产一区二区在线观看av| 黄色视频在线播放观看不卡| 亚洲av电影在线进入| 亚洲欧洲日产国产| 国产乱人偷精品视频| 亚洲av国产av综合av卡| 色婷婷久久久亚洲欧美| 亚洲成色77777| 欧美中文综合在线视频| 国产亚洲精品第一综合不卡| 一区在线观看完整版| 色婷婷久久久亚洲欧美| 人人妻人人澡人人爽人人夜夜| 男女边吃奶边做爰视频| 国产野战对白在线观看| 超色免费av| av国产久精品久网站免费入址| 国产精品国产av在线观看| 国产成人91sexporn| 麻豆乱淫一区二区| 午夜福利网站1000一区二区三区| 一区二区三区乱码不卡18| 亚洲视频免费观看视频| 亚洲综合色网址| 国产成人欧美| 色精品久久人妻99蜜桃| 又黄又粗又硬又大视频| 捣出白浆h1v1| 成人漫画全彩无遮挡| 男女午夜视频在线观看| 99热国产这里只有精品6| 久久久久精品性色| 国产av码专区亚洲av| 欧美久久黑人一区二区| 成年动漫av网址| 婷婷色综合大香蕉| 中文字幕av电影在线播放| 亚洲欧美精品综合一区二区三区| www.熟女人妻精品国产| av不卡在线播放| 国产老妇伦熟女老妇高清| 一区福利在线观看| 国产精品av久久久久免费| 女人精品久久久久毛片| 精品国产国语对白av| 亚洲伊人色综图| 国产成人精品久久久久久| 亚洲av福利一区| 97精品久久久久久久久久精品| 午夜福利一区二区在线看| 国产极品粉嫩免费观看在线| 最近的中文字幕免费完整| 在线天堂最新版资源| 国产日韩一区二区三区精品不卡| av视频免费观看在线观看| 亚洲人成电影观看| 久久综合国产亚洲精品| 丝瓜视频免费看黄片| 一区福利在线观看| tube8黄色片| 日韩一本色道免费dvd| 亚洲国产精品一区三区| 麻豆乱淫一区二区| 国产av精品麻豆| 成人影院久久| 老熟女久久久| 久久免费观看电影| 亚洲欧美一区二区三区久久| 一区二区日韩欧美中文字幕| 91国产中文字幕| 咕卡用的链子| 在线观看一区二区三区激情| 在线观看免费日韩欧美大片| 在线观看免费高清a一片| 欧美在线黄色| 日韩伦理黄色片| 老司机亚洲免费影院| 国产在线免费精品| 国产精品三级大全| 悠悠久久av| a级毛片在线看网站| 精品卡一卡二卡四卡免费| 考比视频在线观看| 一区二区av电影网| 最近2019中文字幕mv第一页| 亚洲成av片中文字幕在线观看| 哪个播放器可以免费观看大片| 丰满少妇做爰视频| 日韩欧美一区视频在线观看| 视频在线观看一区二区三区| 欧美日韩视频高清一区二区三区二| 妹子高潮喷水视频| 捣出白浆h1v1| 亚洲欧美成人综合另类久久久| 日韩免费高清中文字幕av| 欧美av亚洲av综合av国产av | 亚洲精品国产色婷婷电影| 亚洲av综合色区一区| 不卡视频在线观看欧美| 久久久精品国产亚洲av高清涩受| 亚洲自偷自拍图片 自拍| 精品视频人人做人人爽| 这个男人来自地球电影免费观看 | 亚洲精品美女久久av网站| 欧美日韩一级在线毛片| 精品少妇内射三级| 老汉色av国产亚洲站长工具| xxxhd国产人妻xxx| 久久精品国产亚洲av涩爱| 激情五月婷婷亚洲| 国产精品久久久人人做人人爽| 波多野结衣av一区二区av| 欧美日韩一级在线毛片| 久久久久久久国产电影| 亚洲av成人精品一二三区| 亚洲精品中文字幕在线视频| 两性夫妻黄色片| 亚洲欧美激情在线| 视频在线观看一区二区三区| 亚洲欧美中文字幕日韩二区| 男女下面插进去视频免费观看| 一本大道久久a久久精品| 老司机影院毛片| 中文字幕高清在线视频| 夫妻性生交免费视频一级片| 欧美乱码精品一区二区三区| 一级毛片黄色毛片免费观看视频| av电影中文网址| 99精国产麻豆久久婷婷| 亚洲精品国产av成人精品| 成人亚洲欧美一区二区av| 国精品久久久久久国模美| 亚洲成av片中文字幕在线观看| 制服人妻中文乱码| 人人妻人人爽人人添夜夜欢视频| 校园人妻丝袜中文字幕| 国产欧美亚洲国产| 热re99久久国产66热| 性高湖久久久久久久久免费观看| 丝袜美足系列| 欧美最新免费一区二区三区| 美国免费a级毛片| 美女高潮到喷水免费观看| 久久久久久久精品精品| 岛国毛片在线播放| 久久国产亚洲av麻豆专区| 黄色一级大片看看| 中文字幕av电影在线播放| 在线天堂中文资源库| 天堂8中文在线网| 欧美黑人欧美精品刺激| 久久97久久精品| 看免费成人av毛片| 国产成人免费无遮挡视频| 国产精品二区激情视频| 人人澡人人妻人| 国产毛片在线视频| 日韩成人av中文字幕在线观看| 国产精品二区激情视频| 日日爽夜夜爽网站| 黄色一级大片看看| 热re99久久国产66热| 午夜福利免费观看在线| 午夜福利影视在线免费观看| 国产福利在线免费观看视频| 久久久久精品性色| 国产精品国产三级专区第一集| 亚洲熟女毛片儿| 国产黄色免费在线视频| 国产精品久久久久久精品古装| 超色免费av| 日韩一区二区三区影片| 青草久久国产| 亚洲成人一二三区av| 久久天躁狠狠躁夜夜2o2o | 亚洲美女搞黄在线观看| 日韩不卡一区二区三区视频在线| 亚洲av在线观看美女高潮| 日本wwww免费看| 韩国av在线不卡| 日韩大码丰满熟妇| 飞空精品影院首页| 精品一区二区免费观看| 91成人精品电影| 国产精品人妻久久久影院| 亚洲国产毛片av蜜桃av| 777米奇影视久久| 美女午夜性视频免费| 亚洲国产中文字幕在线视频| 欧美激情 高清一区二区三区| 在线观看国产h片| 亚洲成人国产一区在线观看 | 国产一级毛片在线| 免费在线观看黄色视频的| 国产精品久久久久久精品古装| 亚洲国产精品一区二区三区在线| 啦啦啦视频在线资源免费观看| 欧美日韩亚洲国产一区二区在线观看 | av不卡在线播放| 免费看不卡的av| 国产麻豆69| 老汉色∧v一级毛片| 久久久久久久久久久久大奶| 午夜日韩欧美国产| 99国产精品免费福利视频| www.精华液| 久久久国产欧美日韩av| 国产欧美日韩综合在线一区二区| 一二三四在线观看免费中文在| 欧美日韩一区二区视频在线观看视频在线| 天美传媒精品一区二区| 两个人免费观看高清视频| 国产一区亚洲一区在线观看| 亚洲av成人不卡在线观看播放网 | 午夜福利乱码中文字幕| 成年女人毛片免费观看观看9 | 亚洲色图综合在线观看| 成人亚洲欧美一区二区av| 不卡av一区二区三区| 天天躁夜夜躁狠狠久久av| 日韩中文字幕欧美一区二区 | 久久韩国三级中文字幕| 成年美女黄网站色视频大全免费| 在线免费观看不下载黄p国产| 亚洲美女黄色视频免费看| 国产爽快片一区二区三区| 亚洲av男天堂| 在线观看免费视频网站a站| 制服人妻中文乱码| 成人免费观看视频高清| 免费日韩欧美在线观看| 超色免费av| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕免费大全7| 日韩精品免费视频一区二区三区| 91老司机精品| 亚洲成人av在线免费| 五月天丁香电影| 亚洲精品久久久久久婷婷小说| 悠悠久久av| 亚洲av男天堂| 成人午夜精彩视频在线观看| 国产99久久九九免费精品| 国产高清不卡午夜福利| 国产片特级美女逼逼视频| 十八禁人妻一区二区| 国产欧美日韩一区二区三区在线| 日本vs欧美在线观看视频| 日韩 欧美 亚洲 中文字幕| 青春草亚洲视频在线观看| 国产一区二区激情短视频 | 亚洲成人手机| 国产女主播在线喷水免费视频网站| 美女国产高潮福利片在线看| 天天操日日干夜夜撸| 国产亚洲一区二区精品| 人人妻人人澡人人看| 国产精品麻豆人妻色哟哟久久| 欧美少妇被猛烈插入视频| 亚洲在久久综合| netflix在线观看网站| 午夜福利网站1000一区二区三区| 啦啦啦啦在线视频资源| 国产成人免费观看mmmm| 日韩欧美一区视频在线观看| 亚洲色图综合在线观看| 亚洲欧美精品自产自拍| 99热全是精品| 国产一区二区在线观看av| 久久狼人影院| 观看美女的网站| 亚洲成色77777| 亚洲一区中文字幕在线| 哪个播放器可以免费观看大片| 国产日韩欧美亚洲二区| 亚洲欧洲精品一区二区精品久久久 | 日韩大码丰满熟妇| av线在线观看网站| 黑人猛操日本美女一级片| 亚洲国产日韩一区二区| 免费日韩欧美在线观看| 国产精品一国产av| 欧美精品高潮呻吟av久久| 激情五月婷婷亚洲| 久久久国产精品麻豆| 亚洲一区中文字幕在线| 制服诱惑二区| 免费黄频网站在线观看国产| 青春草亚洲视频在线观看| 90打野战视频偷拍视频| 下体分泌物呈黄色| 操出白浆在线播放| 日韩视频在线欧美| 91老司机精品| 亚洲第一区二区三区不卡| 青春草国产在线视频| 精品久久久久久电影网| 国产精品国产av在线观看| 美女福利国产在线| 国产精品 国内视频| 中文字幕亚洲精品专区| 美女扒开内裤让男人捅视频| 精品少妇黑人巨大在线播放| 2021少妇久久久久久久久久久| 亚洲精品,欧美精品| 女人高潮潮喷娇喘18禁视频| 中文字幕最新亚洲高清| 人人妻人人澡人人爽人人夜夜| 王馨瑶露胸无遮挡在线观看| 亚洲精品日韩在线中文字幕| 熟妇人妻不卡中文字幕| svipshipincom国产片| 国产精品.久久久| 丝袜脚勾引网站| 亚洲精品久久久久久婷婷小说| 麻豆精品久久久久久蜜桃| 国产成人精品福利久久| 国产极品天堂在线| 一二三四在线观看免费中文在| 一级,二级,三级黄色视频| 最近2019中文字幕mv第一页| 大话2 男鬼变身卡| 高清欧美精品videossex| 精品酒店卫生间| 91精品伊人久久大香线蕉| 久久久久久久久久久免费av| 国产av国产精品国产| 精品少妇黑人巨大在线播放| 国产亚洲最大av| 欧美激情 高清一区二区三区| 国产日韩一区二区三区精品不卡| 久久久久久久久久久免费av| 久热爱精品视频在线9| 久久这里只有精品19| 国产精品国产av在线观看| 99久久99久久久精品蜜桃| 日韩人妻精品一区2区三区| 亚洲美女搞黄在线观看| 久久久久视频综合| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 国产又爽黄色视频| 日韩精品有码人妻一区| 亚洲三区欧美一区| 最新的欧美精品一区二区| 国产亚洲一区二区精品| 大香蕉久久成人网| 激情视频va一区二区三区| 成人手机av| 亚洲精品在线美女| 国产熟女午夜一区二区三区| 老汉色∧v一级毛片| 少妇被粗大猛烈的视频| 美女视频免费永久观看网站| 成人手机av| 亚洲精品在线美女| 亚洲一卡2卡3卡4卡5卡精品中文| 999久久久国产精品视频| 精品一区二区免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av福利一区| 午夜福利乱码中文字幕| 男女高潮啪啪啪动态图| 国产免费又黄又爽又色| 日韩av免费高清视频| av视频免费观看在线观看| 成人18禁高潮啪啪吃奶动态图| 国产欧美亚洲国产| 无限看片的www在线观看| 国产福利在线免费观看视频| 久久久久久人妻| 久久ye,这里只有精品| 香蕉丝袜av| 欧美成人精品欧美一级黄| 多毛熟女@视频| 午夜福利在线免费观看网站| 亚洲精品在线美女| 一边亲一边摸免费视频| 亚洲av电影在线进入| 精品国产超薄肉色丝袜足j| 一级片'在线观看视频| 国产成人精品在线电影| 亚洲精品在线美女| 国产男女内射视频| 99精国产麻豆久久婷婷| 婷婷成人精品国产| 成人国产麻豆网| 国产又爽黄色视频| 又大又黄又爽视频免费| 久久ye,这里只有精品| 国产精品久久久人人做人人爽| 在线观看国产h片| 亚洲男人天堂网一区| 校园人妻丝袜中文字幕| 天天添夜夜摸| 国产成人欧美| 精品一区在线观看国产| 国产伦人伦偷精品视频| 久久人人97超碰香蕉20202| 成人漫画全彩无遮挡| 国产成人精品无人区| 国产乱人偷精品视频| xxxhd国产人妻xxx| 亚洲一区二区三区欧美精品| 国产日韩欧美在线精品| 午夜日韩欧美国产| 成人黄色视频免费在线看| 国产精品久久久人人做人人爽| 超碰成人久久|