• <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日韩精品久久久久久密| 久久久国产成人免费| 精品国产美女av久久久久小说| 999久久久精品免费观看国产| 99国产极品粉嫩在线观看| 在线永久观看黄色视频| 日韩三级视频一区二区三区| 制服诱惑二区| 亚洲在线自拍视频| 一边摸一边做爽爽视频免费| 精品国产美女av久久久久小说| 成在线人永久免费视频| 在线观看免费日韩欧美大片| 国产无遮挡羞羞视频在线观看| 亚洲专区中文字幕在线| 日韩三级视频一区二区三区| 久久久久国内视频| 国产精品欧美亚洲77777| 亚洲 欧美一区二区三区| 精品福利观看| 久久亚洲精品不卡| 免费观看人在逋| 亚洲人成电影观看| 成人18禁高潮啪啪吃奶动态图| 久久精品亚洲精品国产色婷小说| 久久久久久人人人人人| 美女高潮喷水抽搐中文字幕| 国产精品成人在线| 国产aⅴ精品一区二区三区波| av电影中文网址| 精品电影一区二区在线| 女人久久www免费人成看片| 18在线观看网站| 老司机午夜福利在线观看视频| 精品一区二区三区av网在线观看| 亚洲专区国产一区二区| 国产精品国产av在线观看| 男女之事视频高清在线观看| 色综合欧美亚洲国产小说| 咕卡用的链子| 国产男靠女视频免费网站| 捣出白浆h1v1| 在线观看免费日韩欧美大片| 午夜免费鲁丝| 国产精品综合久久久久久久免费 | 成年女人毛片免费观看观看9 | 亚洲精品久久成人aⅴ小说| 亚洲av电影在线进入| 制服诱惑二区| 久久人人爽av亚洲精品天堂| 在线观看www视频免费| 两个人免费观看高清视频| 满18在线观看网站| 国产男女内射视频| av不卡在线播放| 99国产极品粉嫩在线观看| 亚洲中文av在线| 久久精品亚洲av国产电影网| 91老司机精品| 国产视频一区二区在线看| 亚洲av成人不卡在线观看播放网| 真人做人爱边吃奶动态| 老汉色∧v一级毛片| 欧美人与性动交α欧美软件| 变态另类成人亚洲欧美熟女 | 一个人免费在线观看的高清视频| 欧美日韩中文字幕国产精品一区二区三区 | tube8黄色片| 午夜福利一区二区在线看| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产精品久久久不卡| 国产精品自产拍在线观看55亚洲 | a级毛片黄视频| 欧美日韩亚洲高清精品| 欧美激情极品国产一区二区三区| 在线观看免费午夜福利视频| 9191精品国产免费久久| 欧美av亚洲av综合av国产av| 夜夜爽天天搞| 老司机福利观看| 免费日韩欧美在线观看| 黄色女人牲交| 一级片免费观看大全| 成熟少妇高潮喷水视频| 亚洲五月色婷婷综合| 午夜福利一区二区在线看| 天天操日日干夜夜撸| 国产单亲对白刺激| 国产精品99久久99久久久不卡| 免费看a级黄色片| 女人久久www免费人成看片| 久久亚洲真实| 在线观看日韩欧美| 最近最新免费中文字幕在线| 日本五十路高清| 建设人人有责人人尽责人人享有的| 精品卡一卡二卡四卡免费| 午夜免费鲁丝| 国产成人影院久久av| www.999成人在线观看| 在线观看66精品国产| 国产成人免费无遮挡视频| 免费av中文字幕在线| 老司机福利观看| 人人妻人人添人人爽欧美一区卜| av不卡在线播放| 王馨瑶露胸无遮挡在线观看| 欧美激情久久久久久爽电影 | 电影成人av| 后天国语完整版免费观看| 免费不卡黄色视频| 免费av中文字幕在线| 真人做人爱边吃奶动态| 午夜激情av网站| 18禁观看日本| 在线观看免费视频网站a站| a级片在线免费高清观看视频| 亚洲av成人av| 亚洲一码二码三码区别大吗| 在线观看免费日韩欧美大片| 十分钟在线观看高清视频www| 日韩欧美在线二视频 | 亚洲第一青青草原| 国产精品av久久久久免费| 日韩大码丰满熟妇| 精品久久久久久电影网| 亚洲国产中文字幕在线视频| 99精品在免费线老司机午夜| 91九色精品人成在线观看| 久久香蕉精品热| 日韩制服丝袜自拍偷拍| 国产精品香港三级国产av潘金莲| 亚洲av美国av| 性色av乱码一区二区三区2| 亚洲色图 男人天堂 中文字幕| 午夜老司机福利片| 国产精品 欧美亚洲| 国产成人精品久久二区二区免费| 99国产精品一区二区三区| 黄色 视频免费看| 国产男靠女视频免费网站| 国产一卡二卡三卡精品| 亚洲九九香蕉| 午夜精品在线福利| 中出人妻视频一区二区| 91大片在线观看| 热re99久久国产66热| 午夜激情av网站| 大型黄色视频在线免费观看| 日韩欧美一区视频在线观看| 国产三级黄色录像| av超薄肉色丝袜交足视频| 老司机午夜十八禁免费视频| 欧洲精品卡2卡3卡4卡5卡区| av一本久久久久| 成人黄色视频免费在线看| 久久99一区二区三区| 亚洲七黄色美女视频| 日韩人妻精品一区2区三区| 色精品久久人妻99蜜桃| 精品无人区乱码1区二区| www.精华液| 热re99久久国产66热| 亚洲熟女毛片儿| 一级a爱片免费观看的视频| 91成年电影在线观看| 成年人黄色毛片网站| 电影成人av| 亚洲成人国产一区在线观看| 久久人妻福利社区极品人妻图片| 亚洲av日韩在线播放| 国产国语露脸激情在线看| 久久人妻av系列| 国产精品国产高清国产av | 91成年电影在线观看| 在线观看午夜福利视频| 黄色丝袜av网址大全| 国产精品亚洲一级av第二区| 自线自在国产av| 国产男女超爽视频在线观看| aaaaa片日本免费| 精品国内亚洲2022精品成人 | 精品亚洲成国产av| 91成年电影在线观看| 男人的好看免费观看在线视频 | 大香蕉久久网| tocl精华| xxxhd国产人妻xxx| 亚洲 国产 在线| 成人亚洲精品一区在线观看| 一级作爱视频免费观看| 国产精品一区二区在线不卡| 一边摸一边做爽爽视频免费| 日韩欧美一区视频在线观看| 欧美日韩亚洲综合一区二区三区_| 国产精华一区二区三区| 国产片内射在线| 国产乱人伦免费视频| 免费黄频网站在线观看国产| 免费女性裸体啪啪无遮挡网站| 精品久久久久久电影网| 9191精品国产免费久久| 极品少妇高潮喷水抽搐| 亚洲精品粉嫩美女一区| 国产精华一区二区三区| 黄色怎么调成土黄色| 999精品在线视频| 天堂√8在线中文| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜免费鲁丝| 日韩欧美国产一区二区入口| 91av网站免费观看| 成熟少妇高潮喷水视频| 亚洲欧美日韩另类电影网站| 99精品欧美一区二区三区四区| 麻豆av在线久日| 国产aⅴ精品一区二区三区波| 久久精品亚洲熟妇少妇任你| 亚洲欧美色中文字幕在线| 男女免费视频国产| 成熟少妇高潮喷水视频| 下体分泌物呈黄色| 国产精品久久久人人做人人爽| 亚洲av欧美aⅴ国产| 成人影院久久| 亚洲精品国产区一区二| av有码第一页| x7x7x7水蜜桃| 午夜日韩欧美国产| 精品亚洲成国产av| av网站免费在线观看视频| 国产男女内射视频| 国产精品99久久99久久久不卡| 亚洲国产看品久久| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 中文字幕制服av| 国产av精品麻豆| 一级片'在线观看视频| 欧美日韩精品网址| 久久中文看片网| 亚洲熟女精品中文字幕| 色综合婷婷激情| 亚洲少妇的诱惑av| 国产精品久久久久久人妻精品电影| 两人在一起打扑克的视频| 精品高清国产在线一区| 精品国产超薄肉色丝袜足j| videos熟女内射| 亚洲av成人av| 日韩人妻精品一区2区三区| 夜夜躁狠狠躁天天躁| 老熟女久久久| 午夜福利在线免费观看网站| 女人精品久久久久毛片| 老熟妇乱子伦视频在线观看| 久热这里只有精品99| 国产成人影院久久av| 最近最新中文字幕大全电影3 | 亚洲av成人一区二区三| 18禁观看日本| 午夜亚洲福利在线播放| 国产精品久久电影中文字幕 | 黄色女人牲交| 男人操女人黄网站| 怎么达到女性高潮| 老司机靠b影院| 国产免费av片在线观看野外av| 视频在线观看一区二区三区| 日日摸夜夜添夜夜添小说| 人人妻,人人澡人人爽秒播| 日本精品一区二区三区蜜桃| 免费观看a级毛片全部| 精品国内亚洲2022精品成人 | 国产一卡二卡三卡精品| 精品国产乱码久久久久久男人| 国产蜜桃级精品一区二区三区 | 色在线成人网| 日韩熟女老妇一区二区性免费视频| 亚洲黑人精品在线| 中出人妻视频一区二区| av电影中文网址| 亚洲av成人不卡在线观看播放网| 一二三四社区在线视频社区8| a在线观看视频网站| 黄片播放在线免费| 美女视频免费永久观看网站| 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| 亚洲精品久久午夜乱码| 免费在线观看影片大全网站| 啦啦啦 在线观看视频| 国产成人一区二区三区免费视频网站| 啦啦啦视频在线资源免费观看| 免费看a级黄色片| 欧美乱妇无乱码| 国产精品国产高清国产av | 精品卡一卡二卡四卡免费| 国产精品久久久人人做人人爽| 视频区图区小说| 国产精品亚洲一级av第二区| а√天堂www在线а√下载 | 一级片免费观看大全| 亚洲午夜理论影院| 欧美激情极品国产一区二区三区| 亚洲精品国产色婷婷电影| 亚洲美女黄片视频| 免费在线观看亚洲国产| 天天躁夜夜躁狠狠躁躁| 大香蕉久久成人网| 老司机在亚洲福利影院| а√天堂www在线а√下载 | 亚洲国产欧美一区二区综合| 99国产精品一区二区三区| 国产一区二区激情短视频| 国产高清国产精品国产三级| 日韩 欧美 亚洲 中文字幕| 国产成人一区二区三区免费视频网站| 日韩中文字幕欧美一区二区| 天堂中文最新版在线下载| 中文字幕av电影在线播放| 国产蜜桃级精品一区二区三区 | 国产高清激情床上av| 日韩有码中文字幕| 免费日韩欧美在线观看| 精品一区二区三卡| 欧美 亚洲 国产 日韩一| 亚洲欧美激情在线| 男女床上黄色一级片免费看| 在线观看免费日韩欧美大片| 夫妻午夜视频| 国产精品99久久99久久久不卡| 精品国产一区二区三区四区第35| 亚洲av熟女| 欧美国产精品一级二级三级| 久久久久久免费高清国产稀缺| 天堂√8在线中文| 女人精品久久久久毛片| 曰老女人黄片| 99re6热这里在线精品视频| 日本一区二区免费在线视频| 国产成+人综合+亚洲专区| av在线播放免费不卡| 两性夫妻黄色片| 亚洲国产欧美网| 日韩中文字幕欧美一区二区| 19禁男女啪啪无遮挡网站| 女警被强在线播放| 老司机亚洲免费影院| 人妻一区二区av| 日本vs欧美在线观看视频| a级毛片黄视频| 久久中文看片网| 天天操日日干夜夜撸| av中文乱码字幕在线| 美女高潮到喷水免费观看| 后天国语完整版免费观看| 亚洲午夜理论影院| 精品视频人人做人人爽| 中文欧美无线码| 免费一级毛片在线播放高清视频 | 欧美日韩成人在线一区二区| 男男h啪啪无遮挡| 亚洲国产毛片av蜜桃av| 中国美女看黄片| 窝窝影院91人妻| 91成人精品电影| 十八禁人妻一区二区| 男男h啪啪无遮挡| 国产av一区二区精品久久| 国产成人欧美在线观看 | 少妇的丰满在线观看| 久久精品国产a三级三级三级| 成人手机av| 一级黄色大片毛片| 国产黄色免费在线视频| 一级a爱视频在线免费观看| 精品免费久久久久久久清纯 | 欧美在线一区亚洲| 在线十欧美十亚洲十日本专区| 久久久国产精品麻豆| 国产99久久九九免费精品| 99re6热这里在线精品视频| 久久影院123| av视频免费观看在线观看| 婷婷丁香在线五月| 免费少妇av软件| 久久久久久亚洲精品国产蜜桃av| 韩国精品一区二区三区| 欧美日韩一级在线毛片| 亚洲全国av大片| 99精品久久久久人妻精品| 久久婷婷成人综合色麻豆| 欧美乱妇无乱码| 动漫黄色视频在线观看| 男人操女人黄网站| 中文字幕制服av| 久久香蕉激情| 国产欧美日韩一区二区三区在线| 一区在线观看完整版| 欧美日韩黄片免| 欧美日韩成人在线一区二区| 久久久久久免费高清国产稀缺| 中国美女看黄片| 在线观看66精品国产| 色精品久久人妻99蜜桃| 国产精品亚洲一级av第二区| 91国产中文字幕| 国产精品亚洲一级av第二区| 国产精品成人在线| 老汉色∧v一级毛片| 国产av又大| 美女视频免费永久观看网站| 99精品久久久久人妻精品| 午夜视频精品福利| 黄色片一级片一级黄色片| 成人av一区二区三区在线看| 99久久国产精品久久久| 极品少妇高潮喷水抽搐| 少妇的丰满在线观看| 亚洲国产毛片av蜜桃av| 高清av免费在线| 国产淫语在线视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产看品久久| 国产精品 欧美亚洲| 两个人免费观看高清视频| 男女高潮啪啪啪动态图| 国产91精品成人一区二区三区| 国产激情久久老熟女| 国产精品香港三级国产av潘金莲| 成熟少妇高潮喷水视频| 国产av一区二区精品久久| www.熟女人妻精品国产| 国产真人三级小视频在线观看| 免费看十八禁软件| 女警被强在线播放| 可以免费在线观看a视频的电影网站| a级毛片黄视频| 黄片小视频在线播放| 一进一出抽搐动态| 亚洲va日本ⅴa欧美va伊人久久| 久久久精品免费免费高清| bbb黄色大片| 亚洲欧美精品综合一区二区三区| www.999成人在线观看| 亚洲一码二码三码区别大吗| 一本综合久久免费| 亚洲国产精品一区二区三区在线| 成人永久免费在线观看视频| 国产伦人伦偷精品视频| 久久久久国产精品人妻aⅴ院 | 免费观看精品视频网站| 国产精品.久久久| 王馨瑶露胸无遮挡在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美 亚洲 国产 日韩一| 亚洲av熟女| 一二三四在线观看免费中文在| a级毛片在线看网站| 国产视频一区二区在线看| 中文欧美无线码| 欧美日韩亚洲高清精品| 99久久国产精品久久久| tube8黄色片| 欧美另类亚洲清纯唯美| a级片在线免费高清观看视频| 制服人妻中文乱码| 在线天堂中文资源库| 亚洲成人手机| 精品免费久久久久久久清纯 | 中文欧美无线码| 国产精品综合久久久久久久免费 | 日韩三级视频一区二区三区| 9191精品国产免费久久| 在线看a的网站| 国产精品1区2区在线观看. | av不卡在线播放| 国产成人精品久久二区二区91| 女性被躁到高潮视频| 亚洲五月天丁香| 国产精品一区二区精品视频观看| 久久精品成人免费网站| 黄色 视频免费看| 制服诱惑二区| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 久久精品国产亚洲av香蕉五月 | 欧美激情极品国产一区二区三区| 少妇猛男粗大的猛烈进出视频| 国产免费现黄频在线看| 色婷婷久久久亚洲欧美| 免费观看人在逋| 一区二区三区国产精品乱码| 精品国产一区二区三区四区第35| 好男人电影高清在线观看| 99精品在免费线老司机午夜| 日本五十路高清| 欧美色视频一区免费| 9热在线视频观看99| 久久国产精品大桥未久av| 精品久久久久久久久久免费视频 | 亚洲色图综合在线观看| 欧美午夜高清在线| 99久久综合精品五月天人人| 黑人巨大精品欧美一区二区蜜桃| 超色免费av| 在线观看免费视频日本深夜| 人成视频在线观看免费观看| 成人永久免费在线观看视频| 亚洲人成77777在线视频| 精品国内亚洲2022精品成人 | 亚洲一区二区三区欧美精品| 美女国产高潮福利片在线看| 久久人妻福利社区极品人妻图片| 精品久久久久久久久久免费视频 | 午夜视频精品福利| 久久香蕉激情| 成人av一区二区三区在线看| 成人黄色视频免费在线看| 久久国产精品影院| cao死你这个sao货| 亚洲欧美一区二区三区久久| 亚洲男人天堂网一区| 亚洲avbb在线观看| 在线视频色国产色| 无人区码免费观看不卡| 亚洲午夜精品一区,二区,三区| 中文字幕高清在线视频| 精品国产美女av久久久久小说| 无人区码免费观看不卡| 国产一区二区三区视频了| 久久精品成人免费网站| 亚洲欧美精品综合一区二区三区| 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 99re6热这里在线精品视频| 久久影院123| 97超视频在线观看视频| 制服人妻中文乱码| 国产私拍福利视频在线观看| 黄色丝袜av网址大全| 搡女人真爽免费视频火全软件 | 久久精品91无色码中文字幕| 欧美激情久久久久久爽电影| 国产高清有码在线观看视频| www日本在线高清视频| 18禁美女被吸乳视频| 日韩国内少妇激情av| 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 亚洲国产欧美人成| 床上黄色一级片| 久久久国产成人精品二区| 99国产极品粉嫩在线观看| 久久久国产成人免费| 国产黄片美女视频| 嫩草影院精品99| 中文字幕av在线有码专区| 国产一区二区激情短视频| 99精品在免费线老司机午夜| 亚洲熟妇熟女久久| 国产av麻豆久久久久久久| 久久香蕉国产精品| 国产av麻豆久久久久久久| 成年女人永久免费观看视频| a级一级毛片免费在线观看| 久久香蕉国产精品| 精品久久久久久久人妻蜜臀av| 国产高清视频在线观看网站| 国产在线精品亚洲第一网站| 男女床上黄色一级片免费看| 亚洲在线观看片| 亚洲熟妇熟女久久| 窝窝影院91人妻| a级毛片a级免费在线| 亚洲国产欧美网| 国产精品精品国产色婷婷| 欧美乱码精品一区二区三区| 极品教师在线免费播放| 又黄又粗又硬又大视频| 亚洲欧美一区二区三区黑人| 又黄又爽又免费观看的视频| 老鸭窝网址在线观看| 欧美最新免费一区二区三区 | 欧美区成人在线视频| 超碰av人人做人人爽久久 | 日韩人妻高清精品专区| 在线观看日韩欧美| 国产高潮美女av| 无人区码免费观看不卡| 婷婷精品国产亚洲av在线| 日本黄色片子视频| 18禁黄网站禁片免费观看直播| 欧美中文日本在线观看视频| 91av网一区二区| 亚洲不卡免费看| 香蕉av资源在线| 久久久久久久亚洲中文字幕 | 日韩欧美三级三区| 99热这里只有是精品50| 97超视频在线观看视频| 精品久久久久久久久久久久久| eeuss影院久久| 国产成人啪精品午夜网站| 久久久久久久久大av| 国产高清三级在线| 国产成人福利小说| 一区二区三区国产精品乱码| 久久久精品欧美日韩精品| 欧美绝顶高潮抽搐喷水| 欧美在线一区亚洲| 亚洲成人免费电影在线观看| 久9热在线精品视频| 欧美+亚洲+日韩+国产|