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

        ?

        頂部定常吸氣對高層建筑模型氣動力的影響

        2017-12-12 02:50:18王漢封李石清
        關(guān)鍵詞:方柱柱體氣動力

        彭 思, 王漢封,2, 李石清

        (1.中南大學(xué) 土木工程學(xué)院,長沙 410075; 2.高速鐵路建造技術(shù)國家工程實(shí)驗(yàn)室(中南大學(xué)), 長沙410075)

        頂部定常吸氣對高層建筑模型氣動力的影響

        彭 思1, 王漢封1,2, 李石清1

        (1.中南大學(xué) 土木工程學(xué)院,長沙 410075; 2.高速鐵路建造技術(shù)國家工程實(shí)驗(yàn)室(中南大學(xué)), 長沙410075)

        為減少高層建筑受到的風(fēng)荷載,提高高層建筑抗風(fēng)性能,提出一種頂部吸氣的主動控制方法。在一高寬比為H/d=5的正方形截面柱體的高層建筑模型上,通過在其頂部前邊緣開設(shè)狹縫進(jìn)行吸氣來實(shí)現(xiàn)流動控制。通過風(fēng)洞試驗(yàn)研究了頂部吸氣對氣動力與頂部分離流的影響. 利用流動可視化與流場測試結(jié)果揭示了頂部狹縫吸氣的控制機(jī)理. 實(shí)驗(yàn)對比了不同吸氣系數(shù)Q(=U/U∞,U為狹縫吸氣速度) 對氣動力控制效果的影響, 并對Q=0, 1和3的三種工況下風(fēng)壓分布與頂部剪切流進(jìn)行了詳細(xì)對比. 研究發(fā)現(xiàn)狹縫吸氣改變了頂部分離流特性, 并對模型所有高度上的氣動力均有顯著影響.Q=1時控制效果最佳, 脈動阻力與脈動升力分別減小17.8%和45.5%. 此時頂面分離流被削弱并伴隨有再附現(xiàn)象且湍流強(qiáng)度較大, 最有利于頂部剪切流與尾流間的動量交換, 從而削弱柱體展向渦脫落與脈動氣動力.

        空氣動力學(xué);有限長方柱體;定常吸氣;流動分離;流動控制

        針對超高層建筑的氣動優(yōu)化措施大體可以分為被動控制和主動控制[1]. 被動控制最常見的是結(jié)構(gòu)物的外形優(yōu)化, 如:采用倒角、圓角、凹角或角部擾流片, 柱體貫穿開洞, 沿高度改變建筑橫截面面積, 增加棱柱體頂面棱邊和側(cè)面螺旋線的數(shù)量[2-7]. 相比于被動控制, 主動控制方法以其可控性更能適用于復(fù)雜的風(fēng)環(huán)境. 其中吹、吸氣方法是常見的主動控制措施之一. 文獻(xiàn)[8]在高寬比H/d為6的高層建筑模型不同高度上開設(shè)吸氣孔, 研究了不同風(fēng)向角下模型側(cè)面和背面定常吸氣的控制效果. 實(shí)驗(yàn)發(fā)現(xiàn)側(cè)面吸氣比背面吸氣效果更好, 且開孔距流動分離點(diǎn)越近, 控制效果越好. 文獻(xiàn)[9-11]通過數(shù)值模擬研究了高層建筑側(cè)面全高和分段的吸、吹氣控制, 發(fā)現(xiàn)氣流角度、開孔位置和開孔寬度對控制效果都有影響, 但主要受無量綱的吸氣量系數(shù)控制. 從能耗和折減效率的角度分析, 發(fā)現(xiàn)吸、吹氣有不同的適用范圍, 且分段控制效果不如全高控制. 文獻(xiàn)[12-13]進(jìn)一步探討了被動控制和主動控制相結(jié)合的方法, 研究了方形、Y形截面形式、錐臺形式和帶切角或凹角形式的高層建筑模型在0°至45°風(fēng)向角下通過兩側(cè)吸氣的控制效果. 實(shí)驗(yàn)發(fā)現(xiàn), 除錐臺模型外, 其他各吸氣模型對橫風(fēng)向脈動力都有很好的抑制效果.

        有研究發(fā)現(xiàn), 有限長柱體頂部形成下掃流對其氣動力有很大影響[14-16], 但是能否通過控制下掃流實(shí)現(xiàn)對高層建筑氣動力的控制仍有待進(jìn)一步驗(yàn)證. 本文通過風(fēng)洞試驗(yàn), 在正方形截面的高層建筑簡化模型頂部施加定常吸氣, 系統(tǒng)地研究了吸氣對氣動力的影響規(guī)律, 并通過流動可視化、流場測量等方法對頂部吸氣控制的機(jī)理進(jìn)行了探討.

        1 實(shí)驗(yàn)介紹

        所有試驗(yàn)在一直流式風(fēng)洞內(nèi)進(jìn)行, 其實(shí)驗(yàn)段截面尺寸為450 mm×450 mm, 風(fēng)速范圍0~40 m/s, 湍流度小于0.5%. 被測高層建筑簡化為一正方形截面棱柱, 寬d=40 mm, 高寬比H/d=5. 模型由韌性類ABS樹脂經(jīng)3D打印制成, 剛度及表面光滑度都能滿足實(shí)驗(yàn)要求. 安裝模型前, 對風(fēng)洞壁面邊界層進(jìn)行標(biāo)定, 其厚度約為0.4d, 即模型0.4d以上高度都處于均勻來流中. 圖1(a)給出了模型示意圖與坐標(biāo)系的定義. 坐標(biāo)系原點(diǎn)為方柱底面中心. 模型頂面迎風(fēng)側(cè)距邊緣1 mm處開設(shè)有一寬1 mm長36 mm的吸氣細(xì)縫. 模型底部的出氣管與放置在風(fēng)洞外的抽氣泵相連, 并可通過流量計監(jiān)測吸氣流量, 見圖1(b).在z*= 1, 2, 3, 4和4.5等(z*=h/d,h為測點(diǎn)高度,d為方柱寬度) 不同高度上, 布置五層測壓孔. 每層各個面的測壓點(diǎn)為5個, 沿柱體表面均勻布置.測壓試驗(yàn)使用美國Scannivalve公司ZOC33/64PxX2高頻壓力掃描閥進(jìn)行測量, 測壓試驗(yàn)采樣頻率為625 Hz, 采樣時間32 s, 測壓系統(tǒng)見圖1(b). 實(shí)驗(yàn)中自由來流風(fēng)速U∞= 9.8 m/s, 基于d與U∞的Re=2.74×104. 通過流量計監(jiān)測吸氣的體積流量, 并依據(jù)吸氣狹縫面積, 估算狹縫吸氣速度U, 定義吸氣系數(shù)Q=U/U∞.試驗(yàn)研究了Q=0~4范圍內(nèi)頂部狹縫吸氣對模型氣動力與頂部剪切流的影響規(guī)律.

        (1)

        (2)

        式中:pi表示各測點(diǎn)風(fēng)壓,p為風(fēng)洞靜壓,ρ為空氣密度,為測壓點(diǎn)脈動風(fēng)壓的均方根值(RMS值). 根據(jù)測點(diǎn)代表的面積, 可計算得到模型總體所受阻力系數(shù)Cd,whole, 脈動阻力系數(shù)和脈動升力系數(shù).

        圖1 實(shí)驗(yàn)裝置

        2 結(jié)果分析與討論

        2.1 氣動力

        圖4給出了模型各高度上局部的升力系數(shù)自功率譜. 為方便對比, 所有工況下各層的功率譜均以無控制時z*=1的譜峰值強(qiáng)度作為參考值,按式(3)無量綱化為dB. 頻率則以來流速度U∞和方柱寬度d按式(4)無量綱.

        (3)

        (4)

        式中:pQ=i,z*=j為Q=i,z*=j時功率譜峰值強(qiáng)度,pQ=0,z*=1為Q=0,z*=1時功率譜峰值強(qiáng)度,f為頻率.

        圖3 Q=0、1和3時方柱氣動力沿展向的分布

        Fig.3 Distribution of the aerodynamic forces at different spanwise positions withQ=0, 1 and 3

        圖4 方柱在Q=0、1和3時的各層升力系數(shù)自功率譜

        Fig.4 Power spectral density of the lift at different spanwise positions withQ=0,1 and 3

        無論吸氣與否, 模型所有高度上的無量綱渦脫頻率都為0.116, 這表明頂部吸氣雖然對脈動升力等有明顯削弱作用, 但其不會改變模型渦脫落頻率.Q=0時, 模型所有高度上的升力系數(shù)譜都具有明顯峰值, 說明在不吸氣時, 所有高度上展向渦都

        具有強(qiáng)烈的周期性. 此外, 由于模型頂部所形成的下掃流具有削弱展向渦的特性, 因此, 越接近模型頂部功率譜峰值越弱. 當(dāng)Q=1時, 所有高度上功率譜峰值相對于Q=0時都明顯減弱了, 且越接近頂部削弱得越顯著. 而當(dāng)Q=3時, 功率譜峰值又有所增強(qiáng).

        為進(jìn)一步揭示頂部吸氣對模型氣動力的影響機(jī)理, 圖5給出了模型總體Cd,whole和Cl,whole的時程曲線以及Cl,whole對應(yīng)的時頻分析結(jié)果. 由圖 5可得到:首先,Q=0時,Cl,whole大幅波動的情形在時間上具有明顯的間歇性, 時頻譜中峰值是不連續(xù)的. 這與文獻(xiàn)[18-20]所觀察到的現(xiàn)象一致, 即有限長方柱尾流中會隨機(jī)出現(xiàn)兩種典型狀態(tài): 一種是方柱兩側(cè)展向渦有規(guī)律地交替脫落, 此時Cl,whole大幅振動;另一種是展向渦準(zhǔn)對稱地脫落, 對應(yīng)Cl,whole振幅很小且無明顯周期性. 顯然, 前者所對應(yīng)的升力脈動值要顯著大于后者. 其次,Q=1時,Cd,whole和Cl,whole時程曲線表現(xiàn)得較為平穩(wěn), 這說明類似卡門渦街的交替渦脫模態(tài)在很大程度上被抑制了. 而當(dāng)Q=3時,Cl,whole的周期性波動再次變強(qiáng).

        圖5 模型總體氣動力特性

        2.2 頂部剪切流

        為進(jìn)一步探討頂部剪切流的特性, 進(jìn)行了流動顯示實(shí)驗(yàn). 煙線布置見圖1(a), 發(fā)煙材料為石蠟. 實(shí)驗(yàn)風(fēng)速為1.12 m/s, 并以此調(diào)節(jié)吸氣速度, 研究Q=0、1和3工況下方柱頂部流動特性. 由于方柱沒有雷諾數(shù)效應(yīng), 故該風(fēng)速下的實(shí)驗(yàn)結(jié)果仍能為測壓實(shí)驗(yàn)現(xiàn)象的解釋提供參考. 圖6給出了Q=0、1和3時, 模型頂部剪切流的可視化結(jié)果. 由圖 6可以看出,Q=0時, 剪切流在模型頂部邊緣發(fā)生分離, 分離流直接掠過模型頂部, 在其下游形成下掃流, 并與模型兩側(cè)形成的展向渦發(fā)生作用. 這一下掃流具有將展向渦向兩側(cè)分開的趨勢并能削弱展向渦的強(qiáng)度, 降低其脫落頻率(相對于二維柱體)[21-22]. 當(dāng)頂部前邊緣吸氣系數(shù)Q=1時, 頂部分離流被明顯抑制并向模型頂部靠近. 在狹縫吸氣的作用下, 分離流在頂部的下游部分發(fā)生了再附著. 隨著向下游的發(fā)展, 流動再次在模型頂部后邊緣發(fā)生分離, 并伴隨有較為顯著的不穩(wěn)定性. 從流動可視化圖中可以看出, 該剪切流與模型展向渦在模型后部發(fā)生了較為強(qiáng)烈的相互作用, 這可能是此時脈動升力達(dá)到最小的原因. 當(dāng)Q增大到3時, 模型前邊緣剪切流出現(xiàn)了另一種完全不同的形態(tài). 此時, 頂部前邊緣的流動分離被完全抑制, 流動緊貼頂部表面向下游運(yùn)動, 在頂面尾部發(fā)生分離, 且剪切流沒出現(xiàn)如Q=1時的不穩(wěn)定擺動.

        圖6 Q=0、1和3時方柱頂部流動顯示

        圖7 方柱頂部剪切流測量結(jié)果

        圖 7(b)給出了上述測點(diǎn)位置上的流向脈動速度Iuu.Q=0時, 頂部邊緣形成的分離流的剪切層內(nèi)Iuu最強(qiáng), 但其距離頂面較遠(yuǎn). 前緣吸氣在Q=1時不僅顯著縮小了分離區(qū), 而且仍保持了頂面附近較強(qiáng)的Iuu. 而當(dāng)Q增大到3時, 分離流被完全抑制后, 來自柱體前部的自由流緊貼著頂面流過, 此時頂面附近Iuu與自由流相同, 其強(qiáng)度幾乎為零. 由頂面剪切流的流動可視化與測量結(jié)果可知,Q=0時, 頂面前緣分離流直接掠過模型頂部, 并會在頂面上形成分離泡. 而當(dāng)Q=1時, 頂面上的分離流會發(fā)生再附著并明顯抑制分離泡, 且此時頂面流動具有較高的湍流度, 這有利于模型尾流與外部高速流體進(jìn)行動量交換, 并削弱展向漩渦脫落. 而當(dāng)Q增大到3時, 盡管頂部分離流被完全抑制, 但頂面流動湍流度幾乎為零, 這可能反而削弱了尾流與外部高速流體的動量交換, 所以控制效果反而較Q=1時有所減弱.

        3 結(jié) 論

        采用風(fēng)洞試驗(yàn), 研究了H/d=5的正方形截面高層建筑簡化模型頂部狹縫吸氣對其氣動力與頂部剪切流的控制效果, 得到如下結(jié)論:

        1)頂部狹縫吸氣能夠略微減小模型時均阻力, 并顯著抑制模型脈動阻力與脈動升力. 就本文吸氣設(shè)置情況而言, 當(dāng)吸氣系數(shù)Q=1時, 氣動力控制效果最為顯著, 相對未控制工況, 時均阻力減小3.60%, 脈動阻力與脈動升力分別減小17.85%和45.57%.

        2)頂部狹縫吸氣不僅對頂部附近氣動力有影響, 而且對模型所有高度上的氣動力均有顯著控制效果. 頂部吸氣可削弱模型展向旋渦脫落, 但并未改變其渦脫落頻率.

        3)狹縫吸氣對模型頂部分離流有顯著的加速作用, 并顯著改變了分離流特性.Q=1時, 模型頂部流動分離被明顯削弱, 并在頂面后部發(fā)生再附著;而當(dāng)Q=3時, 流動分離被完全抑制, 且頂面上的邊界層也基本消失.

        [1] TAMURA Y, KAREEM A. Advancedstructural wind engineering[M].[S.l.]: Springer Japan, 2013: 377-383. DOI:10.1007/978-4-431-54337-4.

        [2] KIM Y C, BANDI E K, YOSHIDA A, et al. Response characteristics of super-tall buildings-Effects of number of sides and helical angle[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2015, 145: 252-262. DOI:10.1016/j.jweia.2015.07.001.

        [3] TANAKA H, TAMURA Y, OHTAKE K, et al. Experimental investigation of aerodynamic forces and wind pressures acting on tall buildings with various unconventional configurations[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2012, 107/108: 179-191. DOI:10.1016/j.jweia.2012.04.014.

        [4] KAREEM A, KIJEWSKI T, TAMURA Y. Mitigation of motions of tall buildings with specific examples of recent applications[J]. Wind & Structures And International Journal, 1999, 2(3): 201-251.

        [5] KWOK K C S, BAILEY P A. Aerodynamicdevices for tall buildings and structures[J]. Journal of Engineering Mechanics, 1987, 113(3): 349-365.

        [6] DUTTON R, ISYUMOV N. Reduction of tall building motion by aerodynamic treatments[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1990, 36(1-3): 739-744.

        [7] GU Ming, QUAN Yong. Across-wind loads of typical tall buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2004, 92(13):1147-1165. DOI:10.1016/j.jweia.2004.06.004.

        [8] ZHANG Hongfu, XIN Dabo, OU Jinping. Steady suction for controlling across-wind loading of high-rise buildings[J]. The Structural Design of Tall and Special Buildings, 2016, 25(15): 785-800. DOI:10.1002/tal.1283.

        [9] ZHENG Chaorong, ZHANG Yaochun. Numerical investigation on the drag reduction properties of a suction controlled high-rise building[J]. Journal of Zhejiang University-Science A, 2010, 11(7): 477-487. DOI:10.1631/jzus.A0900593.

        [10]鄭朝榮. 高層建筑風(fēng)荷載吸/吹氣控制的數(shù)值模擬研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2010.

        ZHENG Chaorong. Numerical investigation of wind loads on high-rise buildings controlled by suction/blowing[D]. Harbin: Harbin Institute of Technology, 2010.

        [11]鄭朝榮,張耀春. 分段吸氣高層建筑減阻性能的數(shù)值研究[J]. 力學(xué)學(xué)報,2011, 43(2): 372-380.

        ZHENG Chaorong, ZHANG Yaohchun. Numerical investigation of drag-reducing property for a high-rise building under subsection suction[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(2):372-380.

        [12]張繼同. 主動吸氣與被動氣動控制下超高層建筑風(fēng)荷載特性研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2015.

        ZHANG Jitong. Characteristics of wind loads on super high-rise buildings under active suction control and passive aerodynamic control[D].Harbin: Harbin Institute of Technology, 2015.

        [13]鄭朝榮,張繼同,張智棟. 凹角與吸氣控制下高層建筑平均風(fēng)荷載特性試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報,2016(10): 125-131.

        ZHENG Chaorong, ZHANG Jitong, ZHANG Zhidong. Experimental investigation on characteristics of mean wind loads of high-rise buildings controlled by corner recession and air-suction[J]. Journal of Building Structures, 2016(10): 125-131.

        [14]WANG Hanfeng, ZHOU Yu. The finite-length square cylinder near wake[J]. Journal of Fluid Mechanics, 2009, 638: 453-490. DOI:10.1017/S0022112009990693.

        [15]WANG Hanfeng, ZHOU Yu, CHAN C K, et al. Effect of initial conditions on interaction between a boundary layer and a wall-mounted finite-length-cylinder wake[J]. Physics of Fluids, 2006, 18(6): 065106. DOI:10.1063/1.2212329.

        [16]UFFINGER T, ALI I, BECKER S. Experimental and numerical investigations of the flow around three different wall-mounted cylinder geometries of finite length[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013, 119: 13-27. DOI: 10.1016/j.jweia.2013.05.006.

        [17]CHEN Wenli, LI Hui, HU Hui. An experimental study on a suction flow control method to reduce the unsteadiness of the wind loads acting on a circular cylinder[J]. Experiments in Fluids, 2014, 55:1707. DOI:10.1007/s00348-014-1707-7.

        [18]SATTARI P, BOURGEOIS J A, MARTINUZZI R J. On the vortex dynamics in the wake of a finite surface-mounted square cylinder[J]. Experiments in Fluids, 2012, 52(5): 1149-1167. DOI:10.1007/s00348-011-1244-6.

        [19]王漢封,徐萌霞. 有限長正方形棱柱繞流的雙穩(wěn)態(tài)現(xiàn)象[J]. 力學(xué)與實(shí)踐, 2013, 35(3): 46-50.

        WANG Hanfeng, XU Mengxia. Bistable phenomenon of the flow around a finite-length square prism[J]. Mechanics in Engineering, 2013, 35(3): 46-50.

        [20]王漢封,趙東偉,鄒超,等. 高寬比為5的懸臂正方形棱柱氣動力特性[J]. 實(shí)驗(yàn)流體力學(xué),2015, 29(6): 8-15.

        WANG Hanfeng, ZHAO Dongwei,ZOU Chao, et al. Aerodynamic forces of a cantilevered square prism with aspect ratio of 5[J]. Journal of Experiments in Fluid Mechanics, 2015, 29(6): 8-15.

        [21]SAKAMOTO H, ARIE M. Vortex shedding from a rectangular prism and a circular cylinder placed vertically in a turbulent boundary layer[J]. Journal of Fluid Mechanics, 1983, 126(126): 147-165.

        [22]PARK C, LEE S. Free end effects on the near wake flow structure behind a finite circular cylinder[J].Journal of Wind Engineering and Industrial Aerodynamics, 2000, 88(2): 231-246. DOI:10.1016/S0167-6105(00)00051-9.

        (編輯趙麗瑩)

        Thefree-endsteadysuctioneffectsontheaerodynamicforcesofhigh-risebuildingmodel

        PENG Si1, WANG Hanfeng1,2, LI Shiqing1

        (1.School of Civil Engineering, Central South University, Changsha 410075, China; 2.National Engineering Laboratory for High-speed Railway Construction (Central South University), Changsha 410075, China)

        To reduce wind loads on high-rise buildings, a new active aerodynamic control method named steady suction on the top of the building was introduced. The effects of free-end leading edge steady suction on the aerodynamic forces and free-end separation of a finite-length square cylinder were experimentally investigated in a wind tunnel. Flow visualization and flow field measurement were utilized to reveal the control mechanism. The dependence of aerodynamic forces on the suction ratioQ(=U/U∞) was studied. The pressure distribution and the shear flow above the free end were compared in detail forQ=0, 1 and 3. The results show that the slot suction changes the flow separation on the free end significantly and also the aerodynamic forces on the whole cylinder span. The best control result appears atQ=1, with the fluctuation drag and lift reduced by 17.8% and 45.5%, respectively. At this suction ratio, the flow separation at the leading edge is weakened with reattachment and strong turbulence fluctuation presence on the free end, which results in strong momentum transport between the free-end shear flow and the wake and suppresses the vortex shedding and aerodynamic forces efficiently.

        aerodynamic forces; finite-length square cylinder; steady suction; flow separation; flow control

        10.11918/j.issn.0367-6234.201703123

        O355

        A

        0367-6234(2017)12-0159-06

        2017-03-24

        國家自然科學(xué)基金(11472312)

        彭 思(1990—),男,碩士研究生;

        王漢封(1976—),男,教授,博士生導(dǎo)師

        王漢封,wanghf@csu.edu.cn

        猜你喜歡
        方柱柱體氣動力
        上游切角倒角小間距比串列方柱大渦模擬研究
        振動與沖擊(2022年7期)2022-05-04 05:18:38
        串列多方柱氣動特性的試驗(yàn)研究
        振動與沖擊(2021年9期)2021-05-17 05:32:12
        不同倒角半徑四柱體繞流數(shù)值模擬及水動力特性分析
        海洋工程(2021年1期)2021-02-02 02:48:12
        飛行載荷外部氣動力的二次規(guī)劃等效映射方法
        基于多介質(zhì)ALE算法的柱體高速垂直入水仿真
        談擬柱體的體積
        外注式單體液壓支柱頂蓋與活柱體連接結(jié)構(gòu)的改進(jìn)
        2017年中考數(shù)學(xué)模擬試題(十)
        側(cè)風(fēng)對拍動翅氣動力的影響
        方柱繞流中聚乙烯熔體流變行為的數(shù)值模擬
        中國塑料(2015年10期)2015-10-14 01:13:20
        国产日产桃色精品久久久| 丰满少妇被猛烈进入无码| 日韩欧美第一区二区三区| 扒下语文老师的丝袜美腿| 国产精品女同av在线观看| 欧美精品色婷婷五月综合| 国产精品18久久久| 无套内谢孕妇毛片免费看看| 噜噜噜色97| 麻豆视频黄片在线免费观看| 77777亚洲午夜久久多喷| 久久精品无码一区二区三区免费| 国产欧美日韩a片免费软件| 2022AV一区在线| 精品一区二区三区a桃蜜| 亚洲 日韩 激情 无码 中出| 97人人模人人爽人人喊电影 | 亚洲AV无码成人精品区网页| av在线免费播放网站| 亚洲av熟女中文字幕| 久久综合狠狠综合久久综合88| 欧美疯狂做受xxxx高潮小说| 国内精品久久久久久久久蜜桃| 亚洲精品国产av成人网| 91超精品碰国产在线观看| 中国少妇内射xxxx狠干| 日产精品久久久久久久| 日韩av在线不卡观看| 国产一区二区av免费观看| 欧美成人看片一区二区三区尤物 | 在线观看中文字幕一区二区三区| 日本一区二区三级在线| 国产综合色在线视频区| 国产乱人伦av在线无码| 色二av手机版在线| 久久伊人精品中文字幕有| 成品人视频ww入口| 欧美激情内射喷水高潮| 99在线无码精品秘 入口九色| 黄片国产一区二区三区| 亚洲情综合五月天|