亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        遮擋型格構(gòu)式塔架風力特性試驗研究

        2016-01-15 02:09:38李正良劉欣鵬晏致濤肖正直俞登科
        振動與沖擊 2015年11期

        李正良,劉欣鵬,晏致濤,肖正直,俞登科

        (1.重慶大學土木工程學院,重慶 400045; 2.重慶大學山地城鎮(zhèn)建設(shè)與新技術(shù)教育部重點實驗室,重慶 400045)

        第一作者李正良男,博士, 教授,1963年生

        遮擋型格構(gòu)式塔架風力特性試驗研究

        李正良1,2,劉欣鵬1,晏致濤1,2,肖正直1,2,俞登科1

        (1.重慶大學土木工程學院,重慶400045; 2.重慶大學山地城鎮(zhèn)建設(shè)與新技術(shù)教育部重點實驗室,重慶400045)

        摘要:基于高頻天平測力技術(shù),對遮擋型格構(gòu)式塔架模型進行不同遮擋距離、不同試驗風速下的風洞試驗,研究其對遮擋型格構(gòu)式塔架順風向、橫風向、扭轉(zhuǎn)向平均風力系數(shù)以及對順風向、橫風向基底彎矩功率譜的影響,并在此基礎(chǔ)上利用最小二乘法對兩個方向的功率譜進行了參數(shù)擬合。結(jié)果表明,遮擋距離是影響平均風力系數(shù)的重要影響因素,而試驗風速僅對橫風向平均風力系數(shù)有影響;與平均風力系數(shù)相反,對于基底彎矩功率譜而言,試驗風速對其影響顯著,而遮擋距離在整個頻段上對功率譜密度影響并不大;特定工況下極值脈動功率譜擬合效果較好,可用于遮擋型格構(gòu)式塔架動力響應的頻域計算中。

        關(guān)鍵詞:遮擋型格構(gòu)式塔架;風洞試驗;風力系數(shù);功率譜;解析模型

        基金項目:國家自然科學基金(51308568)

        收稿日期:2014-03-10修改稿收到日期:2014-06-10

        中圖分類號:TU279.7+44;TU973+.213

        文獻標志碼:A

        DOI:10.13465/j.cnki.jvs.2015.11.018

        Abstract:Based on high-frequency-force-balance technique, a model of an occluded latticed tower was tested in a wind tunnel under conditions of 3 kinds of wind speed, and 6 kinds of occluding distance to investigate mean force coefficients and power spectral densities(PSDs) of base moment. Based on these results, the analytical models of PSD for the base moment in two directions were fitted by applying the least-square method. It was shown that the occluding distance is an important factor to affect mean force coefficients; however, the wind speed is an important factor to affect PSDs of base moment; the analytical models of PSD of base moment agree well with the experimental results obtained with wind tunnel tests, they can be used to calculate wind-induced dynamic responses of occluded latticed towers in frequency domain.

        Tests for wind force on occluded latticed towers

        LIZheng-liang1,2,LIUXin-peng1,YANZhi-tao1,2,XIAOZheng-zhi1,2,YUDeng-ke1(1. School of Civil Engineering, Chongqing University, Chongqing 400045, China;2. Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, Ministry of Education, Chongqing 400045, China)

        Key words:occluded latticed tower; wind tunnel test; wind force coefficient; power spectrum; analytical model

        近年來,格構(gòu)式塔架已被廣泛應用于特高壓輸電系統(tǒng)中,故其受風荷載后的作用機理已得到了廣泛關(guān)注[1-4]。

        目前,國內(nèi)外風工程界對格構(gòu)式塔架順風向風致振動的研究已較為成熟。Holmes[5-7]基于準定常假定,對格構(gòu)式塔架順風向風振響應做了較為充分的研究。Bayar[8]通過剛性模型高頻天平測力試驗, 測定了格構(gòu)式塔架的靜態(tài)風效應。Carril Jr[9]等在得到塔架平均和脈動阻力系數(shù)的基礎(chǔ)上, 研究了風向角、擋風系數(shù)以及紊流度對阻力系數(shù)的影響。然而,大量現(xiàn)場實測和氣動彈性模型風洞試驗表明[10-13],格構(gòu)式塔架橫風向風振響應甚至比順風向還要大。鄒良浩[14]等通過半剛性格構(gòu)式塔架模型的高頻測力天平風洞試驗,推導了消除了模型共振對基底力譜的影響的公式,得到了格構(gòu)式塔架的一階廣義荷載譜。梁樞果[15]等基于高頻底座天平測力風洞試驗,得到了三種典型的格構(gòu)式塔架的順風向、橫風向與扭轉(zhuǎn)向一階振型廣義荷載譜與均方根廣義力系數(shù)。

        圖1 塔架結(jié)構(gòu)形式Fig.1 Structure of lattice tower

        然而隨著我國西電東送項目等的展開,在戈壁荒漠地區(qū)出現(xiàn)了許多結(jié)構(gòu)形式新穎的拉線塔。例如:雙柱懸索拉線塔、拉V塔等。對于這一類拉線塔,往往存在格構(gòu)式塔架相互遮擋的情形。對于此類塔架,本文定義為遮擋型格構(gòu)式塔架。

        雖然如上所述眾多研究者對格構(gòu)式塔架進行了大量研究,然而不難發(fā)現(xiàn),上述研究都僅限于無遮擋型格構(gòu)式塔架,對遮擋型格構(gòu)式塔架的風力特性研究卻尚屬空白。所以研究格構(gòu)式塔架在相互遮擋情況下的動力特性很有必要?;诟哳l底座測力天平的風洞試驗,測得格構(gòu)式塔架在不同遮擋間距、不同風速各工況下的順風向、橫風向、扭轉(zhuǎn)向平均風力系數(shù)以及基底彎矩功率譜,試驗結(jié)果可作為相似結(jié)構(gòu)體系抗風設(shè)計參考。

        1風洞試驗簡介

        1.1試驗模型

        本次風洞試驗模型均由鋁合金材料制成,格構(gòu)式塔架的桿件截面均為L型截面形式,幾何縮尺比為1/30,其結(jié)構(gòu)形式見圖1。由于角鋼寬度和厚度尺寸很小,很難完全滿足幾何縮尺比,所以在制作模型的過程中,盡量滿足幾何相似的情況下,為了增加模型剛度,適當增加L型角鋼厚度。

        圖2 風場模擬結(jié)果Fig.2 Simulation of the wind field

        1.2風場模擬

        本次試驗是在西南交通大學XNJD-3風洞完成,該風洞尺寸為22.5 m×4.5 m×36 m最大風速可達16.5 m/s,最低風速可控制在1 m/s。

        本次風洞試驗是在均勻風場和利用尖塔以及粗糙元模擬的B類風場中進行。試驗中通過安裝尖劈、格柵和多排分布的粗糙元,模擬了縮尺比為1/30的B類地貌紊流風場。在模型放置處測得的風場模擬結(jié)果見圖2。結(jié)果表明風洞試驗模擬的風場環(huán)境符合規(guī)范要求。

        1.3測試方法及數(shù)據(jù)處理

        本次試驗模型由兩座剛性塔架模型組成(見圖3)。

        圖3 遮擋塔架模型Fig.3 Model of occluded tower

        試驗時將被測塔架放置在試驗段中心,通過調(diào)整遮擋距離(遮擋距離分別為:5D、10D、15D、20D、25D、30D,D為模型橫截面寬度),在均勻流場中測得不同遮擋距離下的順風向、橫風向、扭轉(zhuǎn)向平均風力系數(shù)(均勻流場中的試驗風速分別為4.5 m/s、6.1 m/s、7.2 m/s),在B類地表粗糙度紊流場中測得在不同遮擋距離下順風向、橫風向基底彎矩功率譜(遮擋條件同均勻流場),試驗風速分別為VH=4.41 m/s、5.47 m/s、6.82 m/s,其中VH為塔頂平均風速。被擋塔架模型上的氣動力由高頻五分量測力天平測得,天平5個通道的采樣頻率均為256 Hz,采樣時間為200 s。由于塔架整體高度較高,為了消除模型一階振型共振對測得的基底力的影響,采用鄒良浩等提出的方法,對試驗數(shù)據(jù)進行處理。

        2試驗結(jié)果及影響因素分析

        2.1塔架體軸平均風力系數(shù)

        (3)

        (4)

        (5)

        式中:U為試驗風速;ρ為空氣密度;S和D分別為模型迎風面輪廓面積和模型底部寬度;φ(擋風系數(shù))=迎風面擋風面積/迎風面輪廓面積=0.326。

        圖5為在均勻流場中,不同風速、不同遮擋距離下測得格構(gòu)式塔架模型無量綱順風向、橫風向以及扭轉(zhuǎn)向平均風力系數(shù)。

        2.1.1順風向風力系數(shù)

        由圖5(a)可知,風速對順風向平均風力系數(shù)影響不大,在三種風速條件下(4.5 m/s、6.1 m/s、7.2 m/s)測得的基底剪力系數(shù)基本一致,但遮擋距離對其影響比較顯著,從圖中可以看出,三種風速下,順風向平均風力系數(shù)隨遮擋距離的增大先減小后增大,在遮擋距離達到15D時其數(shù)值達到最小值,隨著遮擋距離為20D時達到最大值,之后遮擋距離繼續(xù)增大但順風向風力系數(shù)變化不大。由于雙柱懸索拉線塔結(jié)構(gòu)體系上的特點,主要關(guān)注遮擋距離在10~25D之間平均力系數(shù)變化趨勢,在此范圍內(nèi),風力系數(shù)最大值與最小值分別為2.1(20D)、1.8(15D),由此可知,順風向平均力系數(shù)在15D~20D范圍內(nèi)變化顯著,這種情況正好發(fā)生在雙柱懸索拉線塔的塔中部位。

        圖5 遮擋塔架平均風力系數(shù)Fig.5 Mean force coefficients of occluded lattice tower

        2.1.2橫風向以及扭轉(zhuǎn)向風力系數(shù)

        圖5(b)為格構(gòu)式塔架模型無量綱橫風向平均風力系數(shù)隨風速及遮擋距離的變化規(guī)律,由圖可知,橫風向力系數(shù)隨遮擋距離的變化趨勢與順風向規(guī)律基本相同,隨著遮擋距離的不斷增大,其數(shù)值先減小再增大,分別在15D和30D達到最小和最大值,但當遮擋距離大于20D之后,橫風向力系數(shù)變化比較緩慢。與順風向不同的是,橫風向力系數(shù)不僅與遮擋距離有關(guān),試驗風速也對其有著直接的影響,由圖5(b)可知,橫風向力系數(shù),隨著試驗風速的增大而增大,在三種試驗風速下,其最大值發(fā)生在風速為7.2 m/s,遮擋距離為25D時,其值接近0.1。由圖5(c)可知,扭轉(zhuǎn)向平均力系數(shù)總體上很小,遮擋距離以及風速對其影響并不明顯,故可以忽略不計。

        由上述分析可知,遮擋距離對結(jié)構(gòu)順風向以及橫風向風力系數(shù)影響顯著,且順、橫兩個方向風力系數(shù)隨遮擋距離的變化規(guī)律相似,當遮擋距離為15D時,無論是對順風向力系數(shù)亦或橫向風力系數(shù)遮擋效應最明顯。與遮擋距離不同,試驗風速僅對橫風向風力系數(shù)影響顯著,其風力系數(shù)隨試驗風速的增大而增大,特定工況下變化率高達50%,但就總體而言,橫風向力系數(shù)絕對值較小,對結(jié)構(gòu)的整體影響不大。對于扭轉(zhuǎn)向平均風力系數(shù)而言,其絕對值總體上很小,由圖可知遮擋距離以及試驗風速對其影響并不明顯,可以忽略。

        2.1.3順風向風力系數(shù)與各國規(guī)范的對比

        為了與規(guī)范結(jié)果對比,圖6給出了順風向日本標準(AIJ2004)[16]、英國標準(BS8100:1986)[17]、ASCE標準(ASCE7-02)[18]、澳大利亞/新西蘭標準(AS/NZS1170.2-2002)[19]以及印度規(guī)范(IS:875-1987)[20]規(guī)定的格構(gòu)式塔架阻力系數(shù)。

        圖6 各國規(guī)范阻力系數(shù)Fig.6 Drag coefficient ofdifferent countries

        由圖6可知,就順風向而言,大多數(shù)國家規(guī)范(BS8100:1986、ASCE7-02、AIJ2004、AS/NZS1170.2-2002)規(guī)定的阻力系數(shù)較為接近。在擋風系數(shù)φ≤0.4的范圍內(nèi),印度規(guī)范規(guī)定的順風向阻力系數(shù)比其他國家規(guī)定的阻力系數(shù)都要大10%。本次試驗模型擋風系數(shù)φ=0.326,極值阻力系數(shù)為2.1試驗結(jié)果小于大多數(shù)國家規(guī)范給出的阻力系數(shù)。

        2.2塔架風荷載功率譜特性

        這里重點關(guān)注特定風向角下,遮擋距離以及試驗風速對格構(gòu)式塔架風荷載功率譜的影響,所以本次試驗是在0°風向角下,分別考察遮擋距離(5~30D),試驗風速(4.41 m/s、5.47 m/s、6.82 m/s)對格構(gòu)式塔架順風向、橫風向基底彎矩譜的影響。

        2.2.1遮擋距離對基底彎矩譜的影響

        圖7所示風速分別取4.41 m/s、5.47 m/s、6.82 m/s,遮擋距離(5~30D)對格構(gòu)式塔架順風向、橫風向基底彎矩譜的影響。

        圖7 6種遮擋距離下的基底彎矩功率譜Fig.7 Power spectrum densities of 6 kinds of occluding distances

        對于順風向基底彎矩譜,當試驗風速較低時(VH=4.41 m/s),遮擋距離對彎矩譜影響不大,從圖7可知,6種遮擋條件下順風向基底彎矩譜基本相同,其無量綱廣義氣動力譜峰值在0.16左右,與其對應的折減頻率為0.2,低頻段與高頻段的斜率也基本一致。當試驗風速提高至VH=5.47 m/s時,遮擋距離為20D的順風向基底彎矩譜要明顯大于其他遮擋條件下的基底彎矩譜,其譜峰值及其對應的折減頻率分別為0.08、0.18。當達到最高試驗風速VH=6.82 m/s,遮擋距離為25D的順風向基底彎矩譜在整個頻段上最大,其譜峰值及與之對應的折減頻率分別為0.08、0.13。由此可知,在低風速下,遮擋距離對順風向基底彎矩譜影響并不明顯,當風速提高時,某一遮擋距離會明顯大于其他遮擋距離下的功率譜。這一現(xiàn)象應該是由于遮擋塔體的存在,使得順風向漩渦脫落再附著而引起的。當遮擋距離達到某一臨界值時,遮擋塔架造成的脫落漩渦,重新附著于被擋塔架塔身處,而造成該遮擋距離下阻力譜值急劇增大。

        對于橫風向基底彎矩譜,在相同試驗風速下,遮擋距離對基底彎矩譜的影響并不明顯,即使是在較高試驗風速不同遮擋距離下測得的橫風向基底彎矩譜在整個頻段上也大致相同,由此得出,遮擋距離對橫風向基底彎矩譜的影響并不顯著。

        2.2.2風速對基底彎矩譜的影響

        圖8為遮擋距離相同時(由于篇幅原因僅給出典型遮擋距離),不同試驗風速對順風向、橫風向基底彎矩譜的影響??梢钥闯鲈囼烇L速對基底彎矩譜影響顯著,隨著試驗風速的增加,基底彎矩譜峰值頻率前移,譜峰帶寬逐漸減小,譜峰值隨試驗風速的增大而明顯減小。總體來看,順風向、橫風向基底彎矩譜都體現(xiàn)了這一規(guī)律。

        從總體而言,在相同工況下,橫風向基底彎矩譜均大于順風向基底彎矩譜,出現(xiàn)這一現(xiàn)象的原因從遮擋塔體處脫落的漩渦造成被遮塔體處紊流度增大,紊流中自帶的微小漩渦使被遮塔體脫渦變得復雜,使之逐漸包含了較大的脫落頻率范圍,進而造成被遮塔體橫風向渦激振動加劇。

        2.3塔架基底彎矩功率譜模型

        通過數(shù)據(jù)處理并消除模型的一階振型共振影響后,測得在特定風向角下不同試驗風速、不同遮擋距離下的順風向、橫風向基底彎矩譜。由于結(jié)構(gòu)抗風設(shè)計主要考慮的是極值風響應,故僅給出塔頂平均風速為4.41 m/s,遮擋距離為20D的順風向、橫風向基底彎矩解析模型。

        圖8 3種風速下的基底彎矩功率譜Fig.8 Power spectrum densities of 3 kinds of wind speed

        圖9為順風向和橫風向基底彎矩譜的試驗結(jié)果和擬合的經(jīng)驗公式的對比圖,通過比較可知,擬合結(jié)果非常接近。

        2.3.1順風向基底彎矩功率譜

        格構(gòu)式塔架順風向動力風荷載主要是由順風向風紊流引起的,不同風速、不同遮擋距離下測得的順風向基底彎矩譜較相似,呈單峰形狀,與順風向紊流譜的特征相同。采用式(6)表達順風向無量綱基底彎矩功率譜:

        (6)

        (7)

        2.3.2橫風向基底彎矩功率譜

        與順風向基底彎矩譜相比較,格構(gòu)式塔架橫風向基底彎矩功率譜卓越頻率變化不大,其譜峰位置約在折減頻率0.15~0.2,且比較穩(wěn)定。但是其譜峰值已接近1.8,是順風向基底彎矩譜峰值的10倍左右。橫風向脈動荷載主要以渦旋脫落激勵為主,且在本試驗條件下,未出現(xiàn)橫風向氣流再附著引起的第二尖峰,采用式(8)表達橫風向無量綱基底彎矩功率譜:

        圖9 試驗數(shù)據(jù)與荷載譜模型對比Fig.9 Comparison of spectrum between experiment results and proposed formulae

        (8)

        (9)

        3結(jié)論

        (1)格構(gòu)式塔架順風向平均風力系數(shù)受遮擋距離影響較大,且當遮擋距離在15~20d時,其值變化最敏感,遮擋效應最為顯著。但風速對順風向平均風力系數(shù)影響并不明顯。

        (2)格構(gòu)式塔架橫風向平均風力系數(shù)不僅受遮擋距離的影響,風速對其的影響也較顯著,其值隨風速增大而增大,在特定工況下變化率高達50%。但就整體而言,橫風向風力系數(shù)絕對值較小,故與順風向力系數(shù)相比其對該類塔體的影響較小。

        (3)對于順風向基底彎矩功率譜,遮擋距離對其影響在整個頻段上并不明顯,不同遮擋距離下的功率譜函數(shù)在整個頻段上基本一致。但試驗風速對其影響較顯著,其譜峰值隨試驗風速的增加而顯著減小,頻帶變寬。

        (4)橫風向基底彎矩譜功率譜的變化規(guī)律同順風向基底彎矩譜基本一致,但橫風向譜峰值是順風向譜峰值的10倍左右。

        (5)針對兩個方向基底彎矩功率譜的形成機理以及譜形特征,分別用兩種不同的荷載模型進行參數(shù)擬合,且動力風荷載解析模型與風洞試驗結(jié)果總體吻合較好。

        參考文獻

        [1]于秀雷,梁樞果,鄒良浩. 格構(gòu)式塔架順風向風振響應簡化計算[J]. 華中科技大學學報:城市科學版,2006(S2): 8-10.

        YU Xiu-lei, LIANG Shu-guo, ZOU Liang-hao. Along-wind dynamic responses of lattice towers[J]. Journal of Huazhong University of Science and Technology, 2006(S2): 8-10.

        [2]LIANG Shu-guo, ZOU Liang-hao, ZHAO Lin, et al. The generalized force spectrum of along-wind loads on lattice towers[A]. Proceedings of the Sixth Asia-Pacific Conference on Wind Engineering[C]. Soul, Korea, 2005:943-952.

        [3]李正良,肖正直,韓楓,等. 1000 kV漢江大跨越特高壓輸電塔線體系氣動彈性模型的設(shè)計與風洞試驗[J].電網(wǎng)技術(shù),2008(12):1-5.

        LI Zheng-liang, XIAO Zheng-zhi, HAN Feng, et al. Aeroelastic model design and wind tunnel tests of 1000 kV hanjiang long span transmission line system[J]. Power System Technology, 2008(12):1-5.

        [4]梁樞果,鄒良浩,趙林,等. 格構(gòu)式塔架三維動力風荷載的風洞試驗研究[J]. 空氣動力學學報,2007(03): 311-318.

        LIANG Shu-guo, ZOU Liang-hao, ZHAO Lin, et al.The investigation of 3-D dynamic wind loads on lattice towers by wind tunnel test[J]. Acta Aerodynamica Sinica,2007(03): 311-318.

        [5]Holmes J D. Along-wind responses of lattice tower: part I-derivation of expressions for gust response factors[J]. Engineering Structures, 1994, 16(4): 287.

        [6]Holmes J D. Along-wind responses of lattice tower: part aerodynamics damping and deflections[J]. Engineering Structures, 1996, 18(7): 483.

        [7]Holmes J D. Along-wind responses of lattice tower: part-effective load distributions[J]. Engineering Structures, 1996, 18(7): 489.

        [8]Bayar D C. Drag coefficients of latticed towers[J]. Journal of Structural Engineering, 1986, 112(2):417-430.

        [9]Carril jr C F, Isyumov N, Brasil R M L R F. Experimental study of the wind forces on rectangular latticed communication towers with antennas[J]. Journal of Wind Engineering and Industrial Aerodynamics,2003,(91):1007-1022.

        [10]Glanville M J, Kwok K C S. Dynamic charac-teristics and wind induced responses of a steel frame tower[J]. Journal of Wind Engineering and Industrial Aerodynamics,1995(54/55):133.

        [11]Ballio G, Meberini F, Solari G. A 60 year old 100m high steel tower: limit states under wind actions[J]. Journal of Wind Engineering and Industrial Aerodynamics,1992( 43):2120.

        [12]樓文娟,孫炳楠,唐錦春. 高聳格構(gòu)式結(jié)構(gòu)風振數(shù)值分析及風洞試驗[J]. 振動工程學報, 1996, 9(3):318.

        LOU Wen-juan, SUN Bing-nan, TANG Jin-chun. Wind tunnel test and numerical computation on wind-induced vibration for tall lattice tower[J] Chinese Journal of Vibration Engineering, 1996, 9(3):318.

        [13]梁樞果. 輸電塔-線體系的氣彈模型風洞實驗與動力特性、風振響應分析[R].上海:同濟大學土木工程防災國家重點實驗室,2002.

        [14]鄒良浩,梁樞果.半剛性模型風洞試驗荷載譜的處理方法[J].實驗流體力學,2007(3):76-81.

        ZOU Liang-hao,LIANG Shu-guo. A method to evaluate wind force spectra of semi-rigid model in wind tunnel tests[J]. Journal of Experiments in Fluid Mechanics,2007(3):76-81.

        [15]梁樞果,鄒良浩,趙林,等.格構(gòu)式塔架動力風荷載解析模型[J]. 同濟大學學報:自然科學版,2008(2):166-171.

        LIANG Shu-guo,ZOU Liang-hao,ZHAO Lin,et al. Analytical model of dynamic wind loads on lattice towers[J]. Journal of Tongji University:Natural Science,2008(2):166-171.

        [16]Architectural Institute of Japan. AIJ Recomm-endations for Loads on Building[S].2004.

        [17]British Standard BS8100:Latticed towers and mast,Part1,Code of Practice for Loading[S].1986.

        [18]ASCE 7-02 Minimum Design Loads for Buil-dings and Other Structures[S].

        [19]Australia/New Zealand Standard AS/NZS1170.2-2002. Minimum Design Loads on Structures part2: Wind Loads[S].

        [20]IS.875-1987.Indian Standard Code of Practice for Design Loads(other than Earthquake)for Bui-ldings and Structures. Part 3-Wind Loads[S].

        一区二区三区在线观看日本视频| 久久久精品人妻一区二区三区四 | 天美传媒精品1区2区3区| 无遮挡十八禁在线视频国产制服网站 | 日韩精品免费观看在线| 在线人成视频播放午夜| 久久综合国产乱子伦精品免费| 亚洲人成人网毛片在线播放| 人妻少妇精品视频中文字幕国语| 那有一级内射黄片可以免费看| 制服丝袜中文字幕在线| 亚洲中文字幕无码中字| 亚洲中文字幕av一区二区三区人 | 国产三级精品三级在线观看粤语 | 青青草精品在线视频观看| 少妇饥渴偷公乱a级无码| 中文字幕无码人妻丝袜| 成年女人午夜特黄特色毛片免| 91中文人妻熟女乱又乱| 国产熟妇高潮呻吟喷水| 国产欧美日韩不卡一区二区三区| 日韩av天堂综合网久久| 亚洲熟妇无码av在线播放| 久久久久亚洲av无码专区体验| 亚洲AV日韩Av无码久久| 中文字幕人妻互换av| 亚洲中文字幕无码av| 欧美巨大精品欧美一区二区| 国内露脸中年夫妇交换| 国产精品nv在线观看| 国产精品成人av大片| 精品国产拍国产天天人| 久久久伊人影院| 91桃色在线播放国产| av无码精品一区二区三区| 18禁超污无遮挡无码免费游戏 | 欧美成人网视频| 给我看免费播放的视频在线观看| 国产无遮挡又黄又爽免费网站 | 久久婷婷五月国产色综合| 男女边吃奶边做边爱视频|