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

    風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)與單個(gè)葉片模態(tài)參數(shù)對(duì)比研究*

    2022-02-03 02:59:04李志軍馬江波
    關(guān)鍵詞:塔架風(fēng)力機(jī)振型

    李志軍,董 彬,馬江波,李 剛

    (1.西安工業(yè)大學(xué) 建筑工程學(xué)院,西安 710021;2.西安工業(yè)大學(xué) 機(jī)電工程學(xué)院,西安 710021)

    隨著我國(guó)風(fēng)電技術(shù)的快速發(fā)展,到2020年,我國(guó)的風(fēng)電新增裝機(jī)容量在全球排名第一[1]。風(fēng)力機(jī)主要由塔架、機(jī)艙、葉片、輪轂等結(jié)構(gòu)組成。其中葉片是獲取風(fēng)能的核心結(jié)構(gòu)。隨著風(fēng)力機(jī)大型化,葉片的長(zhǎng)度越來(lái)越大,從而呈現(xiàn)剛度低、阻尼弱的特點(diǎn)[2]。葉片應(yīng)具有更加良好的翼型結(jié)構(gòu)、材料和加工工藝,保證葉片具有足夠的結(jié)構(gòu)剛度和韌性,以便提高風(fēng)能利用率以及承受更大的外部載荷。否則,因工作環(huán)境因素產(chǎn)生的葉片振動(dòng)會(huì)導(dǎo)致葉片開裂或者折翼,更嚴(yán)重時(shí)會(huì)導(dǎo)致葉片折斷報(bào)廢,造成倒塌[3]。

    文獻(xiàn)[4]對(duì)單個(gè)葉片進(jìn)行了模態(tài)分析和風(fēng)致振動(dòng)控制分析,仿真結(jié)果有良好的控制效果,但研究中沒(méi)有考慮機(jī)艙、塔架對(duì)葉片振動(dòng)的影響。文獻(xiàn)[5]對(duì)5 MW風(fēng)力機(jī)單個(gè)葉片基于葉片根部固定進(jìn)行了動(dòng)力響應(yīng)分析,通過(guò)加載不同工況的風(fēng)荷載,結(jié)果表明隨著風(fēng)速的增大,葉片根部的應(yīng)力是增大的,但是葉片根部按理想的固定端約束處理會(huì)引起一定的計(jì)算誤差。文獻(xiàn)[6]對(duì)風(fēng)力機(jī)葉片平面內(nèi)的振動(dòng)進(jìn)行了主動(dòng)控制研究,采用歐拉-拉格朗日法計(jì)算了塔架對(duì)葉片的耦合結(jié)構(gòu)的模態(tài)分析,隨后進(jìn)行了振動(dòng)控制,但未考慮塔架對(duì)葉片振動(dòng)的影響。文獻(xiàn)[7]對(duì)風(fēng)力機(jī)風(fēng)輪的氣動(dòng)特性以及風(fēng)力機(jī)組振動(dòng)特性進(jìn)行了研究,考慮了多個(gè)葉片結(jié)構(gòu),但沒(méi)有考慮塔架對(duì)葉輪的影響。文獻(xiàn)[8]同樣研究了葉輪的流固耦合仿真分析,通過(guò)ANSYS WorkBench軟件搭建了葉輪流固耦合模型,通過(guò)加載風(fēng)荷載,對(duì)風(fēng)輪的整體穩(wěn)定性進(jìn)行研究,但并未結(jié)合實(shí)際工況研究塔架、機(jī)艙對(duì)葉片振動(dòng)的影響。文獻(xiàn)[9]對(duì)大型風(fēng)力機(jī)塔架進(jìn)行了風(fēng)-沙工況響應(yīng)分析,主要探究了風(fēng)力機(jī)塔架葉片耦合結(jié)構(gòu)在風(fēng)沙共同作用下的動(dòng)力學(xué)響應(yīng),對(duì)塔架葉片耦合結(jié)構(gòu)的模態(tài)分析缺少研究。文獻(xiàn)[10]研究了風(fēng)力機(jī)塔架的振動(dòng)特性,分別對(duì)風(fēng)力機(jī)停機(jī)和運(yùn)行兩種工況進(jìn)行了仿真分析,結(jié)果表明風(fēng)力機(jī)的運(yùn)行對(duì)塔架的振動(dòng)頻率有一定的影響。同樣是僅對(duì)塔架進(jìn)行了振動(dòng)響應(yīng)分析,沒(méi)有考慮葉片與塔架的耦合作用??紤]到風(fēng)力機(jī)組工作環(huán)境的復(fù)雜程度,傳統(tǒng)的對(duì)單個(gè)葉片進(jìn)行振動(dòng)分析相比實(shí)際工況下的葉片振動(dòng)有較為明顯的差距,忽略了對(duì)于風(fēng)力機(jī)本身模態(tài)的研究。

    為獲得復(fù)雜工況下的葉片風(fēng)致振動(dòng)響應(yīng),需要對(duì)風(fēng)力機(jī)進(jìn)行整體結(jié)構(gòu)建模,以考慮塔架、機(jī)艙、輪轂對(duì)葉片振動(dòng)的影響,文中建立了塔架-葉片-機(jī)艙的整體三維模型和動(dòng)力學(xué)模型。將機(jī)艙底部固定在塔架上端,將葉片固定在輪轂上,然后通過(guò)機(jī)艙與塔架連接,按產(chǎn)品的實(shí)際參數(shù)建模,從而實(shí)現(xiàn)仿真分析過(guò)程中考慮塔架對(duì)葉片振動(dòng)的耦合作用,并對(duì)單個(gè)葉片和塔架-葉片耦合結(jié)構(gòu)兩種不同的模型進(jìn)行模態(tài)分析,為風(fēng)力機(jī)葉片振動(dòng)控制技術(shù)發(fā)展提出依據(jù)。

    1 風(fēng)力機(jī)塔架-葉片耦合模型

    1.1 風(fēng)力機(jī)參數(shù)選取

    文中采用目前技術(shù)成熟的NREL-5 MW渦輪風(fēng)力機(jī)作為研究對(duì)象,其機(jī)組性能和結(jié)構(gòu)參數(shù)[11-12]具體見表1。

    表1 NREL-5 MW風(fēng)力機(jī)參數(shù)Tab.1 Structural parameters of the NREL-5 MW wind turbine

    1.2 單個(gè)葉片三維模型

    通過(guò)NREL-5 MW渦輪風(fēng)力機(jī)的性能參數(shù),計(jì)算可得到葉片各翼型參數(shù)見表2。對(duì)于Cylinder1,DU40_A17系列以及NACA64_A17翼型坐標(biāo),可通過(guò)翼型設(shè)計(jì)軟件來(lái)對(duì)葉片翼型截面的二維坐標(biāo)進(jìn)行查詢,但是二維坐標(biāo)難以完整表述葉片的三維形狀,因此,需要將二維坐標(biāo)轉(zhuǎn)換成三維坐標(biāo)去進(jìn)行三維建模。坐標(biāo)變換法正是通過(guò)已有的每個(gè)葉片的弦長(zhǎng)和扭矩,將二維坐標(biāo)數(shù)據(jù)進(jìn)行坐標(biāo)變換計(jì)算,得到相應(yīng)的三維坐標(biāo),將三維坐標(biāo)保存為.bat文件,導(dǎo)入到Solidworks中,最終生成葉片三維模型,再對(duì)邊緣接縫處進(jìn)行處理,具體操作步驟如下:通過(guò)坐標(biāo)變換法得到葉片的截面三維空間坐標(biāo),將其保存為.bat文件,.bat文件可以直接加載進(jìn)Solidworks軟件中。通過(guò)一系列的草圖繪制,得到了葉片各截面的草圖,通過(guò)曲線擬合建立完整的葉片三維模型。

    以文中采用的5 MW風(fēng)力機(jī)葉片為例,其額定功率設(shè)計(jì)P=5 MW,葉片截面弦長(zhǎng)等參數(shù)見表2,將從profili中獲取的多種翼型二維坐標(biāo)按照相應(yīng)的截面在Solidworks中繪制草圖,通過(guò)曲線自動(dòng)將各截面端部依次連接起來(lái),從而完成翼型圖的繪制,最終完成葉片模型的建立。葉片模型具體如圖1所示。

    表2 NREL-5 MW風(fēng)力機(jī)葉片結(jié)構(gòu)參數(shù)Tab.2 Structural parameters of the NREL-5 MW wind turbine blades

    圖1 風(fēng)力機(jī)葉片三維模型Fig.1 Three-dimensional model of a wind turbine blade

    葉片結(jié)構(gòu)具體其根部截面最大,葉尖處截面最小,中間段逐漸變小,從結(jié)構(gòu)上看,葉片的葉尖附近承載能力較小,振動(dòng)產(chǎn)生的位移比較大。

    1.3 塔架-葉片耦合三維模型

    NREL-5 MW風(fēng)力機(jī)的塔架性能主要取決于安裝轉(zhuǎn)子和機(jī)艙等的支撐結(jié)構(gòu),支撐結(jié)構(gòu)的類型取決于安裝地點(diǎn),其屬性參數(shù)主要有土壤類型、風(fēng)力等級(jí)等參數(shù)來(lái)確定。文中選取的5 MW風(fēng)力機(jī)的塔架、機(jī)艙、輪轂參數(shù)具體見表3。

    表3 機(jī)艙和輪轂參數(shù)Tab.3 Parameters of the nacelle and the hub

    根據(jù)上述風(fēng)力機(jī)各零部件的結(jié)構(gòu)參數(shù),在Solidworks三維建模軟件中建立并裝配風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu),具體如圖2所示。

    圖2 風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)三維模型Fig.2 Three-dimensional model of the wind turbine tower-blade coupled structure

    1.4 動(dòng)力學(xué)模型

    針對(duì)風(fēng)力機(jī)葉片的振動(dòng)問(wèn)題,采用假設(shè)模態(tài)法,基于拉格朗日方程通過(guò)系統(tǒng)的動(dòng)能、勢(shì)能和一般荷載建立塔架-葉片耦合運(yùn)動(dòng)模型。動(dòng)力學(xué)低階矩陣微分方程如式(1)所示[13-14]。

    (1)

    式中:q={q1,in,q2,in,q3,in,qn,in}T為位移矢量;M為質(zhì)量矩陣;C為阻尼矩陣;K為剛度矩陣;Ql為風(fēng)荷載;Qg為重力荷載;Um為控制力。單個(gè)葉片質(zhì)量矩陣[15]為

    (2)

    單個(gè)葉片阻尼矩陣[15]為

    (3)

    單個(gè)葉片剛度矩陣[15]為

    (4)

    式中:y為葉片劃分為n個(gè)單元中的第y個(gè)微單元?;谀B(tài)空間,單個(gè)葉片質(zhì)量可以表示為

    (5)

    式中:μ(x)為葉片距葉根x處的單位質(zhì)量函數(shù)。葉片面內(nèi)振動(dòng)一階、二階模態(tài)質(zhì)量分別為

    (6)

    (7)

    Mtj,h=3m0+Mjc+mtj,h。

    (8)

    塔架前后振動(dòng)耦合結(jié)構(gòu)總等效質(zhì)量為

    Mtjf=3m0+Mjc+mtjf,1。

    (9)

    風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)的質(zhì)量[16-17]為

    (10)

    式中:O為2維零矩陣;M2=diag(m1,21,m2,21)為葉片模態(tài)質(zhì)量矩陣;Mtj=diag(Mtj1,Mtjf1)為塔架模態(tài)質(zhì)量矩陣;Mjφ=diag(Mj11,Mj12)為葉片與塔架耦合模態(tài)質(zhì)量矩陣;其中Mj11=cos(φj)m1,1i,Mj12=m2,11;φj為葉片在t時(shí)刻的偏轉(zhuǎn)角度。

    (11)

    Cjg=diag[cb,cb,cb,c4],

    (12)

    式中,cb=diag(c1,c2),c4=diag(ctj,ctjf);

    (13)

    (14)

    j=1,2,3,

    式中:K1(ωr)=k1,21;K2(ωr)=k2,21;Kio=kio,11;K4=diag(Ktj,Ktjf)

    式中:k1,21為葉片面內(nèi)彈性重力及離心力導(dǎo)致的剛度;k2,21為葉片面外彈性重力及離心力導(dǎo)致的剛度;kio,11為葉片面內(nèi)與面外耦合導(dǎo)致的剛度。

    2 模態(tài)分析

    2.1 有限元分析模型

    1) 導(dǎo)入三維模型。Solidworks建模完成后保存的文件格式為.SLDPRT文件,通用有限元軟件可以直接打開.SLDPRT文件,在菜單欄中選擇要分析的類型,選擇要分析的模型單擊確定,即可載入要分析的三維模型,在Geometry工作空間中單擊Generate進(jìn)行模型生成[16],生成的模型如圖3所示。

    圖3 導(dǎo)入有限元分析軟件中的模型Fig.3 Finite element models

    考慮葉片平面內(nèi)外的模態(tài)分析,可得出耦合結(jié)構(gòu)的第1,2階葉片面內(nèi)面外的振型函數(shù)為

    φ1,i(x)=-0.6923x6+1.3754x5-3.5963x4+

    2.5421x3+0.3549x2

    φ2,i(x)=-2.3855x6+4.8723x5-3.38528x4+

    1.7102x3+0.0623x2

    1.7528x3+1.3759x2

    3.1796x4+0.6879x2

    2) 基本設(shè)置。結(jié)構(gòu)構(gòu)件的材料、密度、彈性模量及泊松比是形成結(jié)構(gòu)有限元模型的基礎(chǔ)物理參數(shù),各構(gòu)件單元類型是影響結(jié)構(gòu)模型動(dòng)力特性的重要參數(shù)。主要構(gòu)件采用的材料及物理參數(shù)見表4。

    表4 結(jié)構(gòu)材料與物理參數(shù)Tab.4 Material and physical parameters of each structural member

    3) 塔架-葉片的網(wǎng)格劃分及施加約束。設(shè)置好各部分結(jié)構(gòu)材料參數(shù)后,對(duì)模型進(jìn)行網(wǎng)格劃分以及施加約束。為保證結(jié)構(gòu)模態(tài)分析的精確性,文中對(duì)葉片劃分網(wǎng)格時(shí)選用自適應(yīng)網(wǎng)絡(luò)來(lái)進(jìn)行劃分。

    風(fēng)力機(jī)葉片根部固定于輪轂上,設(shè)置水平和扭轉(zhuǎn)約束;輪轂連接軸安裝在機(jī)艙端部,設(shè)置水平約束;艙底部固定于塔架頂部,設(shè)置水平約束和扭轉(zhuǎn)約束。

    2.2 數(shù)值模擬及分析

    大型風(fēng)力機(jī)為避免發(fā)生共振,要避免結(jié)構(gòu)主要模態(tài)頻率接近風(fēng)載主要頻率。對(duì)單個(gè)葉片和考慮塔架耦合作用的葉片進(jìn)行模態(tài)分析,計(jì)算出兩種模型的主要模態(tài)頻率和振型圖,風(fēng)力機(jī)葉片振動(dòng)的能量主要集中在前6階振型內(nèi),因此選取前6階模態(tài)頻率和振型進(jìn)行研究,模態(tài)頻率具體數(shù)值見表5,振型如圖4和圖5。

    表5 單個(gè)葉片和塔架-葉片的耦合結(jié)構(gòu)前6階模態(tài)頻率分析Tab.5 The first six order modal frequencies of the single blade structure and the tower-blade coupled structure

    圖4 單個(gè)葉片結(jié)構(gòu)前6階振型圖Fig.4 The first six vibration modes of a single blade

    圖5 風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)前6階振型圖Fig.5 The first six order vibration mode of the tower-blade coupled structure

    由表5可知,對(duì)于單個(gè)葉片和塔架-葉片耦合結(jié)構(gòu)模型,對(duì)比這兩組固有頻率可以得出:塔架-葉片耦合結(jié)構(gòu)相比于單個(gè)葉片的固有頻率在第一階模態(tài)時(shí)只有單個(gè)葉片的44.44%,在第二階模態(tài)時(shí)只有單個(gè)葉片的32.50%,在第三階模態(tài)時(shí)只有單個(gè)葉片的28.00%,在第四階模態(tài)時(shí)只有單個(gè)葉片的13.37%,在第五階模態(tài)時(shí)只有單個(gè)葉片的13.11%,在第6階模態(tài)時(shí)只有單個(gè)葉片的15.31%。

    由圖4對(duì)單個(gè)葉片的前6階振型分析可知,單個(gè)葉片在前6階模態(tài)中主要以揮舞和擺振兩種振型為主,扭轉(zhuǎn)振型次之,僅發(fā)生在第6階模態(tài)中,但從單個(gè)葉片的前6階振動(dòng)形式來(lái)看,主要振動(dòng)形式還是以平面內(nèi)振動(dòng)和平面外振動(dòng)為主。由單個(gè)葉片振型圖可知,單個(gè)葉片的模態(tài)振動(dòng)最大位移集中發(fā)生在葉片的尖部,到葉片尖部頂端距離約是葉片總長(zhǎng)度的20%。

    經(jīng)過(guò)對(duì)風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)前6階振型圖5進(jìn)行分析,耦合結(jié)構(gòu)的葉片的第1,第2階模態(tài)為揮舞振動(dòng),第3,第4階模態(tài)為擺振,第5,第6階為扭轉(zhuǎn)振動(dòng)。由振型圖可知,耦合結(jié)構(gòu)葉片的模態(tài)位移發(fā)生最大處集中在葉片的尖部。

    通過(guò)對(duì)表5、圖4和圖5的對(duì)比分析可知,單個(gè)葉片和風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)中葉片在前6階的模態(tài)振型有差異??紤]塔架等結(jié)構(gòu)對(duì)葉片的影響,葉片出現(xiàn)揮舞振動(dòng)的模態(tài)階數(shù)要?。坏珕蝹€(gè)葉片和風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)的葉片模態(tài)振動(dòng)位移最大處基本一致,這說(shuō)明耦合結(jié)構(gòu)對(duì)葉片結(jié)構(gòu)的薄弱部位不產(chǎn)生影響。

    從模態(tài)頻率來(lái)看,考慮風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)時(shí),風(fēng)力機(jī)葉片的振動(dòng)頻率比只考慮單個(gè)葉片振動(dòng)的固有頻率要低,考慮到風(fēng)譜頻率主要集中在0~1 Hz的低頻區(qū),1~10 Hz也有一定的比例,這表明在實(shí)際工作環(huán)境中,風(fēng)力機(jī)塔基-葉片耦合結(jié)構(gòu)中的葉片更容易與風(fēng)產(chǎn)生共振,因此考慮風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)對(duì)風(fēng)力機(jī)葉片的振動(dòng)是非常有必要的。

    從振型圖來(lái)看,葉片發(fā)生扭轉(zhuǎn)振動(dòng)的模態(tài)階數(shù)提前了,這說(shuō)明風(fēng)力機(jī)塔架-葉片耦合作用會(huì)導(dǎo)致葉片在工作環(huán)境中更容易發(fā)生最不利振動(dòng)形式,導(dǎo)致葉片損壞等情況。進(jìn)而說(shuō)明考慮風(fēng)力機(jī)塔架-葉片耦合結(jié)構(gòu)對(duì)風(fēng)力機(jī)葉片的振動(dòng)是非常有必要的。

    3 結(jié) 論

    1) 仿真結(jié)果表明,單個(gè)葉片的固有頻率在一階時(shí)最小為0.27 Hz,二階為0.4 Hz,三階為0.75 Hz,四階為1.72 Hz,五階為1.83 Hz,在六階時(shí)最大為3.2 Hz;塔架-葉片耦合結(jié)構(gòu)的葉片固有頻率同樣在一階最小為0.12 Hz,二階為0.13 Hz,三階為0.21 Hz,四階為0.23 Hz,五階為0.24 Hz,在六階最大為0.49 Hz。塔架-葉片耦合結(jié)構(gòu)的葉片固有頻率相比單個(gè)葉片的固有頻率在一階時(shí)減小55.55%,在二階時(shí)減小67.5%,在三階時(shí)減小72%,在四階時(shí)減小86.63%,在五階時(shí)減小86.89%,在六階時(shí)減小84.69%。

    2) 兩種模型的振型圖表明,葉片的振動(dòng)位移最大處均在葉片的前端,兩種工況下葉片在前6階模態(tài)振型圖中均以平面內(nèi)的擺振和平面外的揮舞振動(dòng)為主要振動(dòng)形式。其中單個(gè)葉片在兩種工況下第1,2階模態(tài)為揮舞振動(dòng),第3,4階模態(tài)為擺振,第6階模態(tài)為扭轉(zhuǎn)振動(dòng);第5階模態(tài)略有不同。

    3) 綜合考慮,由于風(fēng)荷載的頻率較低,塔架與葉片的耦合會(huì)導(dǎo)致葉片的固有頻率降低,葉片運(yùn)行時(shí)會(huì)更容易產(chǎn)生風(fēng)致振動(dòng)。因此,在對(duì)葉片進(jìn)行振動(dòng)響應(yīng)分析時(shí),需要考慮塔架、機(jī)艙、輪轂對(duì)葉片風(fēng)致振動(dòng)產(chǎn)生的耦合作用,從風(fēng)力機(jī)整體建模分析葉片的振動(dòng)和振動(dòng)控制。

    猜你喜歡
    塔架風(fēng)力機(jī)振型
    長(zhǎng)征六號(hào)甲火箭矗立在塔架旁
    上海航天(2022年5期)2022-12-05 01:55:46
    關(guān)于模態(tài)綜合法的注記
    縱向激勵(lì)下大跨鋼桁拱橋高階振型效應(yīng)分析
    塔腿加過(guò)渡段輸電塔動(dòng)力特性分析
    基于UIOs的風(fēng)力機(jī)傳動(dòng)系統(tǒng)多故障診斷
    門式起重機(jī)塔架系統(tǒng)穩(wěn)定性分析
    雙塔式低塔架自平衡液壓提升裝置與吊裝技術(shù)
    風(fēng)力發(fā)電機(jī)設(shè)備塔架設(shè)計(jì)探析
    結(jié)構(gòu)振型幾何辨識(shí)及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    大型風(fēng)力機(jī)整機(jī)氣動(dòng)彈性響應(yīng)計(jì)算
    免费不卡的大黄色大毛片视频在线观看| 精品一品国产午夜福利视频| 香蕉精品网在线| 亚洲精品成人av观看孕妇| 激情视频va一区二区三区| 天堂俺去俺来也www色官网| 久久精品国产亚洲av天美| 999久久久国产精品视频| 韩国高清视频一区二区三区| 欧美日韩亚洲高清精品| 成人二区视频| 亚洲少妇的诱惑av| 精品一品国产午夜福利视频| 国产在线视频一区二区| 永久网站在线| 亚洲经典国产精华液单| 亚洲国产精品999| 国产毛片在线视频| 秋霞伦理黄片| 久热这里只有精品99| 十八禁高潮呻吟视频| 精品亚洲成a人片在线观看| 久久久久久久亚洲中文字幕| 成人午夜精彩视频在线观看| 精品人妻熟女毛片av久久网站| 老司机影院毛片| 久久久久久久久久久久大奶| 久久免费观看电影| 哪个播放器可以免费观看大片| 日本色播在线视频| 人妻系列 视频| 久久青草综合色| 嫩草影院入口| 免费日韩欧美在线观看| 免费黄色在线免费观看| 久久精品久久久久久久性| 亚洲精品一区蜜桃| 午夜福利乱码中文字幕| 久久久久久久久久久免费av| 亚洲欧美一区二区三区黑人 | 国产日韩一区二区三区精品不卡| videossex国产| 热99国产精品久久久久久7| 一二三四在线观看免费中文在| 免费看av在线观看网站| 国产精品久久久久久久久免| 久久久久久人妻| 国产xxxxx性猛交| 国产成人精品福利久久| 日本-黄色视频高清免费观看| 日产精品乱码卡一卡2卡三| 最新中文字幕久久久久| 国产精品人妻久久久影院| 天堂8中文在线网| 久久国内精品自在自线图片| 人妻系列 视频| 国产一区有黄有色的免费视频| 国产午夜精品一二区理论片| 久久久国产精品麻豆| 一级毛片黄色毛片免费观看视频| 欧美精品人与动牲交sv欧美| 日韩av不卡免费在线播放| 免费观看a级毛片全部| 哪个播放器可以免费观看大片| 欧美bdsm另类| av免费观看日本| 欧美人与性动交α欧美精品济南到 | 国产精品香港三级国产av潘金莲 | 久久久久网色| 国产精品蜜桃在线观看| 久久久国产欧美日韩av| 伦精品一区二区三区| 日本免费在线观看一区| 国产精品无大码| 精品午夜福利在线看| 一区二区三区四区激情视频| 下体分泌物呈黄色| 天天躁夜夜躁狠狠躁躁| 考比视频在线观看| 国产免费现黄频在线看| 国产熟女欧美一区二区| 亚洲精品aⅴ在线观看| 欧美激情高清一区二区三区 | 男女啪啪激烈高潮av片| 欧美亚洲 丝袜 人妻 在线| 又大又黄又爽视频免费| 一级片'在线观看视频| 99久久综合免费| 纯流量卡能插随身wifi吗| 亚洲视频免费观看视频| 久久精品人人爽人人爽视色| 精品国产一区二区三区久久久樱花| 人体艺术视频欧美日本| 美女大奶头黄色视频| 国产 一区精品| 亚洲欧美成人综合另类久久久| 又粗又硬又长又爽又黄的视频| 人人妻人人添人人爽欧美一区卜| 亚洲成人手机| 亚洲欧美中文字幕日韩二区| 亚洲精品国产av蜜桃| 在线观看免费视频网站a站| 国产成人精品在线电影| 成人国产麻豆网| 黄色毛片三级朝国网站| 国产精品久久久久久精品电影小说| 精品人妻一区二区三区麻豆| 黄片小视频在线播放| 国产深夜福利视频在线观看| 亚洲精品国产av成人精品| 亚洲精品成人av观看孕妇| 免费少妇av软件| 黑丝袜美女国产一区| 国产精品成人在线| 高清黄色对白视频在线免费看| 两个人免费观看高清视频| 少妇的逼水好多| 亚洲久久久国产精品| 色哟哟·www| 咕卡用的链子| 国产爽快片一区二区三区| 国产精品香港三级国产av潘金莲 | 欧美 亚洲 国产 日韩一| 国产无遮挡羞羞视频在线观看| 午夜免费男女啪啪视频观看| 亚洲欧美一区二区三区黑人 | 国产精品一二三区在线看| 老女人水多毛片| 国产精品.久久久| 两性夫妻黄色片| 亚洲国产欧美在线一区| 嫩草影院入口| 女的被弄到高潮叫床怎么办| 国产精品 欧美亚洲| 一级毛片电影观看| videossex国产| 国产精品一二三区在线看| 精品国产一区二区三区四区第35| 免费在线观看完整版高清| 久久精品国产亚洲av高清一级| 成人国语在线视频| 亚洲国产最新在线播放| 蜜桃国产av成人99| kizo精华| 最近最新中文字幕大全免费视频 | 欧美黄色片欧美黄色片| 一级,二级,三级黄色视频| 久久精品久久久久久噜噜老黄| 少妇的丰满在线观看| 免费观看性生交大片5| 妹子高潮喷水视频| 久热这里只有精品99| 亚洲欧美一区二区三区黑人 | www.自偷自拍.com| 亚洲在久久综合| 国产女主播在线喷水免费视频网站| 菩萨蛮人人尽说江南好唐韦庄| 精品国产国语对白av| videos熟女内射| 黄频高清免费视频| 看非洲黑人一级黄片| 美女xxoo啪啪120秒动态图| 最近手机中文字幕大全| 亚洲精品第二区| av视频免费观看在线观看| 精品福利永久在线观看| 女人精品久久久久毛片| 久久久久视频综合| 国产精品一区二区在线观看99| 国产成人精品无人区| 最新的欧美精品一区二区| 爱豆传媒免费全集在线观看| 国产av国产精品国产| 亚洲国产精品成人久久小说| 男人舔女人的私密视频| 亚洲第一青青草原| 男女下面插进去视频免费观看| 亚洲美女视频黄频| 黄片播放在线免费| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 精品少妇黑人巨大在线播放| 亚洲av电影在线进入| 巨乳人妻的诱惑在线观看| 美女高潮到喷水免费观看| 哪个播放器可以免费观看大片| 日韩中字成人| 精品第一国产精品| 国产在视频线精品| 丝瓜视频免费看黄片| 亚洲视频免费观看视频| 亚洲av欧美aⅴ国产| 少妇人妻 视频| 青春草亚洲视频在线观看| 少妇人妻精品综合一区二区| 国产淫语在线视频| 超色免费av| 少妇猛男粗大的猛烈进出视频| 免费黄频网站在线观看国产| 天天躁夜夜躁狠狠久久av| 观看美女的网站| 麻豆乱淫一区二区| 日韩在线高清观看一区二区三区| 精品少妇黑人巨大在线播放| 午夜福利网站1000一区二区三区| 国产一区亚洲一区在线观看| 色播在线永久视频| 亚洲国产精品国产精品| 国产白丝娇喘喷水9色精品| 成人毛片60女人毛片免费| av片东京热男人的天堂| 午夜91福利影院| 91成人精品电影| 精品国产露脸久久av麻豆| 一级毛片 在线播放| 久久国产亚洲av麻豆专区| 一级黄片播放器| 亚洲精华国产精华液的使用体验| 伊人亚洲综合成人网| 欧美日韩亚洲国产一区二区在线观看 | 汤姆久久久久久久影院中文字幕| 国产日韩一区二区三区精品不卡| 777久久人妻少妇嫩草av网站| 男女无遮挡免费网站观看| 久久久亚洲精品成人影院| 久久久久久伊人网av| 亚洲欧美中文字幕日韩二区| 国产色婷婷99| 高清黄色对白视频在线免费看| 久久久精品94久久精品| 国产老妇伦熟女老妇高清| 午夜免费男女啪啪视频观看| 男女下面插进去视频免费观看| 亚洲av福利一区| 午夜老司机福利剧场| 成年女人在线观看亚洲视频| 日韩人妻精品一区2区三区| av免费在线看不卡| 欧美 亚洲 国产 日韩一| 亚洲精品一二三| 99热国产这里只有精品6| 大片免费播放器 马上看| av在线老鸭窝| 国产av码专区亚洲av| 成人国语在线视频| 亚洲国产欧美网| 国产野战对白在线观看| 99精国产麻豆久久婷婷| 成年女人毛片免费观看观看9 | 亚洲欧美成人综合另类久久久| 好男人视频免费观看在线| 丝袜喷水一区| 精品国产露脸久久av麻豆| av网站在线播放免费| 老司机影院成人| h视频一区二区三区| 久久99一区二区三区| 精品一区在线观看国产| 欧美人与性动交α欧美软件| 成年女人在线观看亚洲视频| 最近最新中文字幕大全免费视频 | 亚洲精品美女久久久久99蜜臀 | 国产av一区二区精品久久| 久久久久久久久久久久大奶| 国产高清国产精品国产三级| 亚洲成人一二三区av| 天天影视国产精品| 中文天堂在线官网| 在线观看三级黄色| 9191精品国产免费久久| 亚洲国产毛片av蜜桃av| 制服丝袜香蕉在线| 亚洲欧洲精品一区二区精品久久久 | 免费女性裸体啪啪无遮挡网站| 看非洲黑人一级黄片| 少妇猛男粗大的猛烈进出视频| 九九爱精品视频在线观看| 日本-黄色视频高清免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜脚勾引网站| 99re6热这里在线精品视频| 91aial.com中文字幕在线观看| 亚洲 欧美一区二区三区| 老女人水多毛片| 国产极品天堂在线| xxxhd国产人妻xxx| 久久精品国产a三级三级三级| av网站免费在线观看视频| 欧美日韩视频高清一区二区三区二| 欧美亚洲日本最大视频资源| 精品国产露脸久久av麻豆| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩另类电影网站| 制服丝袜香蕉在线| 久久久久久久久久久免费av| av片东京热男人的天堂| 一本大道久久a久久精品| 成年动漫av网址| 成人毛片a级毛片在线播放| 亚洲精品国产色婷婷电影| 亚洲成人一二三区av| 性高湖久久久久久久久免费观看| 亚洲国产av影院在线观看| 色播在线永久视频| 久久久a久久爽久久v久久| 大香蕉久久成人网| 久久99热这里只频精品6学生| 18禁观看日本| 你懂的网址亚洲精品在线观看| av免费在线看不卡| 亚洲欧美中文字幕日韩二区| 国产精品香港三级国产av潘金莲 | 国产乱来视频区| 黄片小视频在线播放| 国产欧美日韩综合在线一区二区| 国产精品熟女久久久久浪| 欧美激情 高清一区二区三区| 少妇被粗大猛烈的视频| 丝袜美足系列| 免费黄频网站在线观看国产| 国产精品 国内视频| 啦啦啦在线免费观看视频4| 亚洲精品国产av蜜桃| 精品国产露脸久久av麻豆| 国产视频首页在线观看| 久久久久久人人人人人| 狠狠婷婷综合久久久久久88av| 大码成人一级视频| 亚洲欧美一区二区三区国产| 中文字幕人妻丝袜一区二区 | 国产国语露脸激情在线看| 国产成人a∨麻豆精品| 久久久久久人妻| 色网站视频免费| 精品一区在线观看国产| 久久久久久久国产电影| 免费高清在线观看日韩| 久久久久久伊人网av| 777久久人妻少妇嫩草av网站| 99国产精品免费福利视频| 精品一区二区免费观看| 国产片特级美女逼逼视频| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 97人妻天天添夜夜摸| 亚洲精品美女久久久久99蜜臀 | 少妇的逼水好多| 十分钟在线观看高清视频www| 亚洲av在线观看美女高潮| 中文欧美无线码| 少妇熟女欧美另类| 视频在线观看一区二区三区| 欧美日韩精品网址| 欧美精品人与动牲交sv欧美| a 毛片基地| 满18在线观看网站| 国产精品熟女久久久久浪| 免费女性裸体啪啪无遮挡网站| 亚洲综合精品二区| 国产熟女欧美一区二区| 亚洲精品国产av成人精品| 国产av码专区亚洲av| 亚洲av电影在线观看一区二区三区| 一个人免费看片子| 嫩草影院入口| 欧美97在线视频| 国产精品成人在线| 久久精品国产综合久久久| 1024香蕉在线观看| 日韩大片免费观看网站| 久久精品久久久久久久性| 国产成人精品久久久久久| 国产av国产精品国产| av视频免费观看在线观看| 精品一区二区三区四区五区乱码 | 美女国产视频在线观看| 秋霞在线观看毛片| 精品国产一区二区三区久久久樱花| 久久影院123| 国产精品麻豆人妻色哟哟久久| 国产成人免费观看mmmm| 国产野战对白在线观看| 少妇的丰满在线观看| 天天影视国产精品| 免费久久久久久久精品成人欧美视频| 日韩欧美一区视频在线观看| 亚洲av福利一区| 韩国精品一区二区三区| 91aial.com中文字幕在线观看| 日韩一区二区三区影片| 国产一区二区三区av在线| 国产成人欧美| 久久精品国产自在天天线| 只有这里有精品99| 久久精品夜色国产| 在线观看免费日韩欧美大片| 婷婷色av中文字幕| 婷婷色综合大香蕉| 日韩制服丝袜自拍偷拍| 国产成人91sexporn| 国产亚洲av片在线观看秒播厂| 桃花免费在线播放| 国产一区有黄有色的免费视频| 中文精品一卡2卡3卡4更新| 在线观看免费视频网站a站| 男女高潮啪啪啪动态图| 久久av网站| 日本爱情动作片www.在线观看| 日本黄色日本黄色录像| 亚洲成人手机| 免费女性裸体啪啪无遮挡网站| 国产老妇伦熟女老妇高清| 美女大奶头黄色视频| 一边摸一边做爽爽视频免费| 街头女战士在线观看网站| 欧美日韩亚洲高清精品| 精品久久久久久电影网| 一级爰片在线观看| 97在线视频观看| 免费看不卡的av| 国产成人av激情在线播放| 日本黄色日本黄色录像| 久久午夜福利片| 精品国产乱码久久久久久男人| 欧美国产精品va在线观看不卡| av网站免费在线观看视频| 欧美成人午夜精品| 欧美日韩精品成人综合77777| 韩国精品一区二区三区| 成年女人毛片免费观看观看9 | 久久人人爽人人片av| 男女高潮啪啪啪动态图| 国产色婷婷99| 成人18禁高潮啪啪吃奶动态图| 人人妻人人爽人人添夜夜欢视频| 嫩草影院入口| 国产福利在线免费观看视频| 亚洲精品第二区| 国产一区二区三区av在线| 18禁观看日本| 精品午夜福利在线看| 赤兔流量卡办理| 午夜免费鲁丝| 国产高清不卡午夜福利| 国产精品久久久久久精品电影小说| 丝袜人妻中文字幕| 欧美 日韩 精品 国产| 国产成人91sexporn| 亚洲国产看品久久| 桃花免费在线播放| 在线亚洲精品国产二区图片欧美| 国产黄频视频在线观看| 亚洲欧美日韩另类电影网站| 超色免费av| 在线亚洲精品国产二区图片欧美| 少妇的丰满在线观看| av视频免费观看在线观看| 人人妻人人添人人爽欧美一区卜| 成人亚洲欧美一区二区av| 亚洲av免费高清在线观看| 男女边摸边吃奶| 亚洲av综合色区一区| 成人国产麻豆网| 亚洲第一av免费看| 中文字幕最新亚洲高清| 看十八女毛片水多多多| 热99国产精品久久久久久7| 这个男人来自地球电影免费观看 | 国产在线免费精品| 亚洲av免费高清在线观看| 这个男人来自地球电影免费观看 | 亚洲av电影在线进入| 国产老妇伦熟女老妇高清| 国产精品偷伦视频观看了| www.自偷自拍.com| 2018国产大陆天天弄谢| 亚洲精品美女久久久久99蜜臀 | 中文字幕av电影在线播放| 中文字幕人妻丝袜制服| 少妇 在线观看| 久久精品国产自在天天线| 欧美精品高潮呻吟av久久| 91久久精品国产一区二区三区| 看免费成人av毛片| 亚洲四区av| 久久99精品国语久久久| 女人高潮潮喷娇喘18禁视频| av一本久久久久| 欧美激情极品国产一区二区三区| 国产精品熟女久久久久浪| 亚洲精品,欧美精品| 韩国精品一区二区三区| √禁漫天堂资源中文www| 久久鲁丝午夜福利片| 亚洲综合精品二区| 欧美日韩成人在线一区二区| 97在线人人人人妻| 欧美精品av麻豆av| 久久国产精品大桥未久av| 日韩 亚洲 欧美在线| 五月伊人婷婷丁香| 巨乳人妻的诱惑在线观看| 黄片播放在线免费| 欧美激情高清一区二区三区 | 亚洲第一青青草原| 午夜日韩欧美国产| 麻豆av在线久日| 伊人久久大香线蕉亚洲五| 中文精品一卡2卡3卡4更新| 欧美人与善性xxx| 丰满少妇做爰视频| 国产成人精品无人区| 丝袜在线中文字幕| 亚洲av福利一区| 精品少妇黑人巨大在线播放| 亚洲欧美清纯卡通| 亚洲美女黄色视频免费看| 国产日韩欧美亚洲二区| 成人18禁高潮啪啪吃奶动态图| 免费高清在线观看日韩| 少妇熟女欧美另类| 日日撸夜夜添| 欧美激情 高清一区二区三区| 久久精品国产综合久久久| 视频在线观看一区二区三区| 伊人久久国产一区二区| 亚洲欧洲精品一区二区精品久久久 | 少妇的逼水好多| 热re99久久国产66热| 边亲边吃奶的免费视频| 天天躁夜夜躁狠狠躁躁| 欧美国产精品一级二级三级| 伦理电影大哥的女人| 亚洲欧美中文字幕日韩二区| 国产欧美日韩综合在线一区二区| 午夜福利乱码中文字幕| 热re99久久国产66热| 男女边摸边吃奶| 亚洲伊人色综图| 一级毛片电影观看| 黑丝袜美女国产一区| 少妇 在线观看| 少妇人妻久久综合中文| 欧美激情 高清一区二区三区| 一本久久精品| 亚洲精品久久久久久婷婷小说| 精品人妻一区二区三区麻豆| 精品第一国产精品| 永久网站在线| 久热久热在线精品观看| 日本猛色少妇xxxxx猛交久久| 伦理电影大哥的女人| 国产精品免费视频内射| 久久久久久伊人网av| 久久av网站| 欧美国产精品va在线观看不卡| 美女主播在线视频| 亚洲精品久久成人aⅴ小说| 国产免费视频播放在线视频| av免费观看日本| 国产色婷婷99| 久久韩国三级中文字幕| 1024视频免费在线观看| 男女午夜视频在线观看| 色视频在线一区二区三区| 国产成人精品久久二区二区91 | 精品一区二区三卡| 一区二区三区四区激情视频| 久久久亚洲精品成人影院| 亚洲少妇的诱惑av| 在线观看美女被高潮喷水网站| 最近最新中文字幕大全免费视频 | 蜜桃国产av成人99| 日韩av免费高清视频| 亚洲成色77777| 亚洲精品成人av观看孕妇| 中国三级夫妇交换| 在线精品无人区一区二区三| 中文乱码字字幕精品一区二区三区| 最新中文字幕久久久久| 午夜免费男女啪啪视频观看| 婷婷色综合www| 一区在线观看完整版| 日韩欧美一区视频在线观看| 午夜老司机福利剧场| 男人操女人黄网站| 在线观看免费视频网站a站| 精品国产国语对白av| 91aial.com中文字幕在线观看| 日韩电影二区| 久久精品国产综合久久久| 午夜av观看不卡| 精品视频人人做人人爽| av国产久精品久网站免费入址| 国产精品成人在线| 成年动漫av网址| 丰满少妇做爰视频| 成人黄色视频免费在线看| 国产av一区二区精品久久| 午夜福利影视在线免费观看| 国产极品天堂在线| 乱人伦中国视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美成人综合另类久久久| 国产精品不卡视频一区二区| 五月天丁香电影| 极品人妻少妇av视频| 久久精品国产亚洲av涩爱| 狂野欧美激情性bbbbbb| 成人国产av品久久久| 久久免费观看电影| 久久久国产欧美日韩av| 黄色毛片三级朝国网站| 大香蕉久久成人网| 在线观看免费视频网站a站|