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

    高層建筑形狀及布局對(duì)城市街區(qū)行人風(fēng)環(huán)境影響研究

    2021-09-13 13:30:44閆渤文,魏民,鄢喬,程勇,舒臻孺,李秋勝,周緒紅
    關(guān)鍵詞:風(fēng)洞

    閆渤文,魏民,鄢喬,程勇,舒臻孺,李秋勝,周緒紅

    摘? ?要:基于風(fēng)洞試驗(yàn)和計(jì)算流體動(dòng)力學(xué)方法(Computational Fluids Dynamics,CFD)研究高層建筑形狀及布局對(duì)城市街區(qū)行人風(fēng)環(huán)境的影響. 采用最大風(fēng)速比和歸一化加速面積比,定量研究五種高層建筑形狀及四類建筑布局對(duì)城市街區(qū)行人風(fēng)環(huán)境的影響,確定全風(fēng)向下的最優(yōu)建筑形狀以及布局,結(jié)合CFD數(shù)值模擬獲得的全域流場(chǎng)信息,揭示建筑形狀和布局對(duì)于城市街區(qū)行人風(fēng)環(huán)境的影響機(jī)理. 結(jié)果表明:在保持建筑高度和街區(qū)容積率一致的情況下,高層建筑群周邊最大風(fēng)速比不會(huì)隨著建筑形狀和建筑布局的改變而發(fā)生明顯變化. 但建筑形狀和建筑布局會(huì)改變建筑群周邊高風(fēng)速區(qū)域的面積大小,全風(fēng)向下的最優(yōu)建筑形狀和布局分別是Y字形和錯(cuò)列式布局,而最不利形狀和布局分別是H字形和圍合式布局. 不同布局下的方形、H字形及X字形高層建筑群的最不利風(fēng)向均位于斜風(fēng)向,而十字形及Y字形則為正風(fēng)向. 高層建筑群在行人高度處的風(fēng)加速現(xiàn)象主要是由狹管效應(yīng)和角部分離效應(yīng)造成的.

    關(guān)鍵詞:城市街區(qū);行人風(fēng)環(huán)境;建筑形狀;建筑布局;風(fēng)洞;CFD

    中圖分類號(hào):TU312.1? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼:A

    Study on Effects of Tall Building Shape and Layout on

    Pedestrian-level Wind Environment in the Urban Area

    YAN Bowen1,WEI Min1,YAN Qiao1,2,CHENG Yong1,SHU Zhenru3,LI Qiusheng4,ZHOU Xuhong1,5

    (1. Chongqing Key Laboratory of Wind Engineering and Wind Engergy Utilization,

    Chongqing University,Chongqing 400045,China;

    2. Hubei United Investment Group Co,Ltd,Wuhan 430061,China;

    3. School of Civil Engineering,Central South University,Changsha 410075,China;

    4. Department of Architecture and Civil Engineering,City University of Hong Kong,Hong Kong 999077,China;

    5. College of Civil Engineering,Hunan University,Changsha 410082,China)

    Abstract:This paper investigates the influences of tall building shape and layout on the pedestrian-level wind (PLW) environment in the urban area by combing the wind tunnel test and Computational Fluids Dynamics (CFD) simulations. The maximum wind speed-up ratio and integrated normalized speed-up area ratio were used to quantify the effects of five different building shapes and four building layouts on the PLW environment. The favored building shape and layout with omnidirectional equal-probability distribution were determined, and the underlying mechanism of building shape and layout that influences the PLW environment in the urban area were elucidated according to the whole flow-domain information attained by CFD simulations. The PLW of groups of Y-shaped tall buildings in the staggered layout was most favored, while the most unfavorable scenario was H-shaped and enclosed layout; In addition, the most unfavorable wind directions of square, H-shaped and X-shaped tall buildings were in the oblique direction, while these of the cross- and Y-shaped buildings were in the normal wind direction. The wind speed-up phenomena of groups of tall buildings at the pedestrian-level were mainly attributed to the flow separation at building corners and channeling effects.

    Key words:urban area;pedestrian-level wind(PLW) environment;building shape;building layout;wind tunnel;CFD(Computational Fluids Dynamics)

    高層建筑物的存在會(huì)顯著改變城市區(qū)域的空氣流動(dòng),形成包括撞擊、繞流、分離及尾流區(qū)域等多種復(fù)雜的鈍體空氣動(dòng)力學(xué)現(xiàn)象[1-3]. 與此同時(shí),隨著城市化進(jìn)程的不斷加快,高密度城市區(qū)域建筑物間的氣動(dòng)干擾現(xiàn)象愈加突出,導(dǎo)致行人高度處出現(xiàn)顯著的高風(fēng)速區(qū)域,進(jìn)一步加劇了建筑群行人風(fēng)環(huán)境的惡化,影響行人舒適度,更有甚者會(huì)危及行人安全,帶來行人風(fēng)環(huán)境問題[4]. 因此,有必要對(duì)高層建筑物周邊區(qū)域的行人風(fēng)環(huán)境進(jìn)行深入研究,提出改善高密度城市區(qū)域行人風(fēng)環(huán)境的方法,進(jìn)一步落實(shí)當(dāng)前我國(guó)發(fā)展可持續(xù)的綠色宜居城市的戰(zhàn)略[5].

    國(guó)內(nèi)外學(xué)者對(duì)建筑物行人風(fēng)環(huán)境做了廣泛研究,并逐漸形成了行人風(fēng)環(huán)境的評(píng)估流程. 主要包含4方面:當(dāng)?shù)仫L(fēng)氣候、周邊地形影響、建筑氣動(dòng)特性以及行人風(fēng)環(huán)境評(píng)估標(biāo)準(zhǔn)[6]. 其中,建筑氣動(dòng)特性對(duì)行人風(fēng)環(huán)境的影響是當(dāng)前的研究熱點(diǎn)[7]. 學(xué)者們對(duì)此展開了廣泛的試驗(yàn)研究. Tsang等[8]對(duì)4個(gè)并列布置的矩形高層建筑開展了風(fēng)洞試驗(yàn),研究了長(zhǎng)寬比、建筑間距對(duì)高層建筑周邊行人高度處風(fēng)速的影響. Xu等[9]構(gòu)建了40種不同形狀的單體高層建筑,研究了建筑形狀對(duì)行人風(fēng)環(huán)境的作用機(jī)理.受到測(cè)點(diǎn)數(shù)量的限制,風(fēng)洞試驗(yàn)中無法提供建筑群周圍及其全域詳細(xì)的三維流場(chǎng). 近些年,隨著計(jì)算流體動(dòng)力學(xué)方法(Computational Fluids Dynamics,CFD)的快速發(fā)展以及高性能計(jì)算資源的增長(zhǎng),雷諾平均湍流模型(Reynolds-Averaged Navier-Stokes,RANS)、大渦模擬方法(Large-eddy Simulation,LES)以及分離渦模擬(Detached-eddy Simulation,DES)已在城市行人風(fēng)環(huán)境舒適性評(píng)估和研究中得到應(yīng)用. 雖然LES和DES模型能在尾流和分離區(qū)獲得比RANS模型更準(zhǔn)確的瞬態(tài)流場(chǎng)結(jié)果,但由于RANS模型計(jì)算效率高,且在高風(fēng)速區(qū)能提供可靠的平均風(fēng)速模擬結(jié)果,仍被廣泛應(yīng)用于行人風(fēng)環(huán)境的研究和工程實(shí)踐[10-12]. Iqbal和Chan[13]結(jié)合風(fēng)洞試驗(yàn)和CFD數(shù)值模擬研究了建筑間距和風(fēng)向?qū)κ中胃邔咏ㄖ盒腥孙L(fēng)環(huán)境的影響. Van等[14]結(jié)合風(fēng)洞試驗(yàn)驗(yàn)證了多種RANS湍流模型,并采用RANS模型研究了雨棚、露臺(tái)以及透風(fēng)樓層等氣動(dòng)外形修正措施對(duì)單棟高層建筑行人風(fēng)環(huán)境的影響. 但目前行人風(fēng)環(huán)境研究的主要研究對(duì)象還只是單棟或者2~3棟建筑,并且現(xiàn)有研究中建筑物外形多以矩形或者方形為主,缺乏建筑布局以及外形對(duì)高層建筑群風(fēng)環(huán)境的影響研究.

    基于以上研究現(xiàn)狀,本文結(jié)合風(fēng)洞試驗(yàn)和CFD數(shù)值模擬,定量分析了等概率全風(fēng)向下建筑外形及布局對(duì)高層建筑群風(fēng)環(huán)境的影響,明確了最優(yōu)建筑外形和布局,并進(jìn)一步基于CFD全域流場(chǎng)結(jié)果揭示了建筑外形和布局對(duì)高層建筑群行人風(fēng)環(huán)境的影響機(jī)理.

    1? ?研究方法

    1.1? ?風(fēng)洞試驗(yàn)

    本次風(fēng)洞試驗(yàn)在香港城市大學(xué)邊界層風(fēng)洞中進(jìn)行(見圖1(a)),試驗(yàn)段尺寸為11.0 m×4.0 m×2.0 m(長(zhǎng)×寬×高). 試驗(yàn)來流依據(jù)日本規(guī)范所規(guī)定的Ⅲ類地貌[15](見圖1(b)). 試驗(yàn)所采用的模型均由PVC材料制成,幾何縮尺比為1 ∶ 400;每個(gè)高層建筑群模型由8個(gè)相同的單棟高層建筑模型構(gòu)成,為了保證容積率相同,每一個(gè)單體建筑具有相同的高度H(210 mm)和平面投影面積(見圖2),模型的阻塞比小于3%. 為了研究建筑形狀及布局對(duì)行人風(fēng)環(huán)境的影響,本文參考了城市街區(qū)常見的5種建筑形狀和4種建筑布局,共20個(gè)工況(見表1),風(fēng)向角間隔30°. 結(jié)合本文采用的建筑布局和建筑外形,考慮實(shí)際試驗(yàn)方案的可行性,在街道、轉(zhuǎn)角以及建筑模型周邊布置了26個(gè)測(cè)點(diǎn),測(cè)點(diǎn)高度為距地面2.0 m(縮尺后為0.5 cm)的行人高度處(見圖3). 風(fēng)速測(cè)量采用kanomax風(fēng)速探頭,精度為±0.1 m/s,采樣頻率選用625 Hz.

    1.2? ?CFD數(shù)值模擬

    CFD數(shù)值模擬中,所有模型的尺寸均與風(fēng)洞試驗(yàn)保持一致. 計(jì)算域尺寸為15H(長(zhǎng))×10H(寬)×6H(高)(如圖4所示),CFD數(shù)值模擬的阻塞率小于3%,滿足日本風(fēng)規(guī)范規(guī)定的數(shù)值模擬計(jì)算域大小的要求[15],無需對(duì)結(jié)果進(jìn)行修正[16-17]. 邊界條件設(shè)定見表2,入口設(shè)置為速度入口(Velocity-inlet),出口設(shè)置壓力出口(Pressure-outlet),兩側(cè)邊界及頂部邊界均采用對(duì)稱性邊界條件(Symmetry),建筑壁面和地面采用無滑移壁面(No-slip wall). 網(wǎng)格劃分采用結(jié)構(gòu)化網(wǎng)格,在建筑物壁面處對(duì)網(wǎng)格加密,首層網(wǎng)格高度為0.000 2 m,建筑壁面首層網(wǎng)格y+為30左右,網(wǎng)格增長(zhǎng)率為1.1,所有工況網(wǎng)格總量為7×106~9×106.

    由于RANS湍流模型會(huì)顯著影響數(shù)值模擬結(jié)果的精度,為了驗(yàn)證不同RANS湍流模型,本文根據(jù)Xu等人在東京工藝大學(xué)(Tokyo Polytechnic University,TPU)邊界層風(fēng)洞中開展的方形單體建筑行人風(fēng)環(huán)境試驗(yàn)結(jié)果[9],分別從定性和定量角度驗(yàn)證了不同RANS湍流模型的模擬精度,包括標(biāo)準(zhǔn)k-ε模型、Realizable k-ε模型、RNG k-ε模型、標(biāo)準(zhǔn)k-ω模型以及SST k-ω模型. 其中標(biāo)準(zhǔn)k-ε模型中的湍流參數(shù)依據(jù)文獻(xiàn)[12]進(jìn)行了修正,驗(yàn)證湍流模型參數(shù)對(duì)高風(fēng)速區(qū)模擬準(zhǔn)確性的影響.

    本文采用基于有限體積法的ANSYS/Fluent 15.0 CFD數(shù)值模擬平臺(tái),相應(yīng)的湍動(dòng)能、湍流耗散率及平均風(fēng)速剖面等入口邊界條件通過編寫自定義函數(shù)(User-Define Function,UDF)實(shí)現(xiàn);求解器為基于壓力求解的不可壓縮流穩(wěn)態(tài)算法,速度-壓力耦合方式為SIMPLEC,動(dòng)量方程和湍流模型輸運(yùn)方程的非線性對(duì)流項(xiàng)離散格式為二階迎風(fēng)格式(Second order Upwind Scheme,SUS);模擬收斂準(zhǔn)則為所有變量的殘差變化穩(wěn)定,最終觀察到k及ε殘差達(dá)到10-6以下,連續(xù)方程殘差達(dá)到10-4以下,且在關(guān)鍵監(jiān)測(cè)點(diǎn)風(fēng)速值達(dá)到平穩(wěn).

    1.3? ?來流自保持性

    行人風(fēng)環(huán)境CFD數(shù)值模擬研究中,由于行人高度距離地面較近,來流特性易受到地面粗糙度的影響而無法保證風(fēng)場(chǎng)特性在來流方向上的一致,影響數(shù)值模擬結(jié)果的準(zhǔn)確性. 因此,基于RANS模型的CFD數(shù)值模擬的關(guān)鍵問題之一是驗(yàn)證來流的自保持性,即流體經(jīng)過地表面時(shí),保證流場(chǎng)沿來流方向的流動(dòng)特性保持一致[19]. 本文通過修正壁面函數(shù)及粗糙度參數(shù)的方法來實(shí)現(xiàn)風(fēng)速的自保持性[20].? 如圖5所示,在空風(fēng)場(chǎng)的校驗(yàn)結(jié)果中,入口和出口的風(fēng)速剖面較為一致,具有較好的自保持性.

    1.4? ?網(wǎng)格無關(guān)性驗(yàn)證

    本節(jié)以工況S-E-180為例(網(wǎng)格總數(shù)為7×106),建立了兩套不同尺寸的網(wǎng)格(稀疏網(wǎng)格首層網(wǎng)格尺寸為0.000 4 m,增長(zhǎng)率為1.1,網(wǎng)格總數(shù)400萬(wàn)左右;加密網(wǎng)格首層網(wǎng)格尺寸為0.000 1 m,增長(zhǎng)率為1.1,網(wǎng)格總數(shù)1 100萬(wàn)左右);3套不同網(wǎng)格對(duì)比見圖6,同時(shí)在圖7中,給出了采用不同網(wǎng)格尺寸得到的26個(gè)測(cè)點(diǎn)模擬結(jié)果.

    圖7結(jié)果表明,稀疏網(wǎng)格與基本網(wǎng)格之間的計(jì)算結(jié)果存在顯著差異,而加密網(wǎng)格與基礎(chǔ)網(wǎng)格之間的差異卻很小. 說明采用基礎(chǔ)網(wǎng)格進(jìn)行模擬不僅可以保證計(jì)算精度,而且可以盡可能小地消耗計(jì)算資源,因此本文所有計(jì)算工況網(wǎng)格均按照基礎(chǔ)網(wǎng)格方式劃分.

    1.5? ?行人風(fēng)環(huán)境評(píng)估方法

    在對(duì)建筑物周邊的行人風(fēng)環(huán)境進(jìn)行評(píng)估時(shí),常常采用風(fēng)速比Ri進(jìn)行分析[21],其相關(guān)定義如下:

    Ri = ■? ? ? ? ?(1)

    式中:Vi為建筑物周邊測(cè)點(diǎn)在行人高度處(本文為離地面高度2 m處)的風(fēng)速;V0為無建筑物時(shí),入口行人高度處風(fēng)速.

    一方面,從定量的角度,本文通過采用誤差度量指標(biāo)q進(jìn)行分析[22],其定義如下:

    q=■■ni,

    ni =1,■≤q或Ri exp-Ri CFD≤q;0,其他.

    (2)

    式中:N為測(cè)點(diǎn)總數(shù);Ri exp、Ri CFD分別代表試驗(yàn)結(jié)果及其對(duì)應(yīng)的數(shù)值模擬結(jié)果;q值最大計(jì)算誤差為20%.

    另一方面,為了定量分析建筑形狀和布局對(duì)高層建筑群行人風(fēng)環(huán)境的影響,本文采用最大風(fēng)速比Rmax和歸一化加速面積比A*? ? ? ?R,avg對(duì)建筑群風(fēng)環(huán)境進(jìn)行量化評(píng)估[9],其相關(guān)定義如下:

    Rmax,θ = max(Ri,θ),A*? ?R,θ = ■? ? ? ? ?(3)

    Rmax = max(Rmax,θ),Rmin = min(Rmax,θ),

    A*? ? ? ?R,avg = ■■A*? ? ?R,θ? ? ? ? ?(4)

    其中AT為評(píng)估區(qū)域的面積,本文為900×900 mm2,如圖8所示.

    由前文定義可知,AR,θ中R值的大小應(yīng)該為一個(gè)大于1.0的數(shù)值. 關(guān)于R的具體取值,在缺乏氣象統(tǒng)計(jì)資料時(shí),要想滿足風(fēng)環(huán)境舒適性,主導(dǎo)風(fēng)向下的風(fēng)速比不宜大于1.2[23]. 因此,本文采用Ri = 1.2,即計(jì)算A1.2,θ的歸一化加速面積.

    2? ?結(jié)果與討論

    2.1? ?CFD結(jié)果驗(yàn)證

    首先,基于TPU風(fēng)洞試驗(yàn)結(jié)果,驗(yàn)證了5種不同的RANS湍流模型,并將其模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比,發(fā)現(xiàn):對(duì)于建筑物兩側(cè)加速區(qū)域的模擬,使用不同的RANS湍流模型,其模擬結(jié)果有所不同,其中RNG k-ε模型以及修正湍流參數(shù)后的標(biāo)準(zhǔn)k-ε模型模擬的效果最好. 限于篇幅,此處不再給出相關(guān)結(jié)果. 這與現(xiàn)有研究中RANS模型對(duì)于高風(fēng)速區(qū)域(Ri > 1.0)的模擬較為準(zhǔn)確的結(jié)論基本一致[10-13]. 考慮到標(biāo)準(zhǔn)k-ε模型計(jì)算效率更高,本文在后續(xù)研究中主要采用經(jīng)過修正湍流參數(shù)的標(biāo)準(zhǔn)k-ε模型開展數(shù)值模擬研究.

    表3為本文開展240個(gè)工況在全部風(fēng)向角下高風(fēng)速區(qū)(Ri > 1.0)的試驗(yàn)與模擬誤差對(duì)比結(jié)果,并與當(dāng)前采用RANS模型開展行人風(fēng)環(huán)境研究的數(shù)值模擬驗(yàn)證結(jié)果進(jìn)行了對(duì)比.

    從表3可以看出:高風(fēng)速區(qū),所有測(cè)點(diǎn)的最大誤差值為20.22%(工況Y-E-150);最小誤差為0.087%(工況H-CO-0);與現(xiàn)有研究相比(鄭朝榮等[24]的最大誤差為22.5%,最小誤差為1.21%,Iqbal等[13]的最大誤差為25.1%,最小誤差為1.59%),本文的模擬精度更高,具有較高的可信度.

    為了進(jìn)一步說明經(jīng)過修正湍流參數(shù)的標(biāo)準(zhǔn)k-ε模型模擬結(jié)果的可信度,采用誤差度量指標(biāo)q進(jìn)行分析[22],q值根據(jù)前文計(jì)算最大誤差選擇20%.

    圖9統(tǒng)計(jì)了所有工況的q值,發(fā)現(xiàn)所有計(jì)算工況的q值均在0.85~0.92之間,表明本文數(shù)值模擬計(jì)算的整體最大誤差小于15%,數(shù)據(jù)置信度在0.85以上. 數(shù)值模擬在高風(fēng)速區(qū)域可以得到與試驗(yàn)較為一致的結(jié)果,但在低風(fēng)速區(qū)RANS模型結(jié)果準(zhǔn)確性則較差,主要是因?yàn)楸疚乃捎玫姆€(wěn)態(tài)雷諾平均模型在建筑物尾流漩渦脫落區(qū)域(低風(fēng)速區(qū))低估了尾流中的湍動(dòng)能,導(dǎo)致數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果相差較大. 現(xiàn)有研究表明:采用RANS模型對(duì)行人風(fēng)環(huán)境進(jìn)行研究時(shí),數(shù)值模擬與試驗(yàn)結(jié)果在高風(fēng)速區(qū)(Ri > 1.0)高度吻合,其誤差控制在0.2左右,低風(fēng)速區(qū)則準(zhǔn)確性較差,與本文模擬結(jié)果相同[5].

    2.2? ?建筑形狀及布局對(duì)行人風(fēng)環(huán)境的影響評(píng)估

    2.2.1? ?建筑形狀對(duì)行人風(fēng)環(huán)境影響

    為了研究建筑形狀對(duì)行人風(fēng)環(huán)境的影響,圖10、圖11分別給出了不同建筑形狀在不同風(fēng)向角下Rmax,θ和A*? ? ? ?1.2,θ的分布情況.

    由于錯(cuò)列式分布的Rmax,θ和A*? ? ? ?1.2,θ分布趨勢(shì)與圍合式分布趨勢(shì)相同,角開口式與中部開口式布局的 Rmax,θ和A*? ? ? ?1.2,θ趨勢(shì)相同,限于篇幅,結(jié)果未給出,可參考圖10和圖11分布情況。

    從圖10 Rmax,θ的分布趨勢(shì)來看:在圍合式布局和角開口式布局下,方形、H字形及X字形建筑群的Rmax出現(xiàn)在斜風(fēng)向下,Rmin則出現(xiàn)在正風(fēng)向下;而十字形及Y字形建筑群則與此相反,Rmax出現(xiàn)在正風(fēng)向下,Rmin則出現(xiàn)于斜風(fēng)向. 對(duì)于角部開口式布局,可以得到相同的結(jié)論,表明這四種布局不影響不同建筑形狀的最有利以及最不利風(fēng)向角. 進(jìn)一步分析比較A*? ? ? ?1.2,θ,其分布規(guī)律與Rmax,θ的分布規(guī)律完全相同,即不同布局下,方形、H字形及X字形的建筑群的最有利風(fēng)向?yàn)檎L(fēng)向,最不利風(fēng)向?yàn)樾憋L(fēng)向;十字形及Y字形建筑群的最有利風(fēng)向?yàn)樾憋L(fēng)向,最不利風(fēng)向?yàn)檎L(fēng)向.

    圖12給出了不同建筑形狀在圍合式布局下,最不利風(fēng)向角下所對(duì)應(yīng)的行人高度處的風(fēng)速比云圖. 通過圖12可以看出:風(fēng)場(chǎng)中的高風(fēng)速區(qū)域集中在建筑群的側(cè)面及建筑物之間形成的天然廊道中. 進(jìn)一步分析,對(duì)于建筑群側(cè)面高風(fēng)速區(qū)域的形成,主要是由于建筑物迎風(fēng)前緣對(duì)風(fēng)場(chǎng)剪切作用和氣流下洗作用聯(lián)合造成的[24];而建筑物通道高風(fēng)速區(qū)域,則是由于建筑物之間形成的通風(fēng)廊道,產(chǎn)生狹管效應(yīng)現(xiàn)象,使得廊道區(qū)域的風(fēng)速明顯增大[25].

    圖13為不同的計(jì)算工況在全風(fēng)向下的最大風(fēng)速比Rmax與歸一化平均加速面積比A*? ? ? ? 1.2,avg的分布情況. 從Rmax和A*? ? ? ? 1.2,avg值的趨勢(shì)來看,對(duì)于這20種不同的計(jì)算工況,Rmax最大值為1.532,最小值為1.468,均位于1.5附近,相互之間差值僅為4.18%;另一方面比較A*? ? ? ? 1.2,avg的大小,明顯發(fā)現(xiàn)不同工況之間A*? ? ? ? 1.2,avg差異很大,最大值為4.02,最小值則為1.67,相互之間的差值超過了140%. 因此,在保持高度和容積率一定的情況下,建筑群周邊最大風(fēng)速比不會(huì)隨著建筑形狀和建筑布局的改變而發(fā)生明顯變化;但建筑形狀和建筑布局會(huì)改變建筑群周邊高風(fēng)速區(qū)域的面積大小,而高風(fēng)速區(qū)域正是本文研究的重點(diǎn). 從Rmax和A*? ? ? ? 1.2,avg值的大小來看,圍合式布局中,方形圍合式布局建筑群的A*? ? ? ? 1.2,avg值最小,因此可認(rèn)為其行人風(fēng)環(huán)境

    整體上最優(yōu);相同的,角開口式及中部開口式布局下,十字形建筑群的整體行人風(fēng)環(huán)境最優(yōu);Y字形建筑群在錯(cuò)列式布局下的行人風(fēng)環(huán)境最優(yōu). 與之前研究結(jié)論基本一致[26-27].

    2.2.2? ?建筑布局對(duì)行人風(fēng)環(huán)境影響

    由于H字形及X字形的Rmax,θ和A*? ? ? ?1.2,θ分布趨勢(shì)與方形完全相同,Y字形則和十字形的Rmax,θ和A*? ? ? ?1.2,θ分布趨勢(shì)相同,因此,本節(jié)僅給出了方形和十字形建筑群在不同布局下的Rmax,θ和A*? ? ? ?1.2,θ分布情況.

    從分布趨勢(shì)來看(如圖14和圖15所示),方形建筑群在四種不同布局下,Rmax對(duì)應(yīng)的風(fēng)向角均為斜方向,而Rmin對(duì)應(yīng)的風(fēng)向角均與建筑物迎風(fēng)面垂直;十字形建筑群在四種布局下的結(jié)論與方形建筑群相反. 這說明在不同布局下,方形、H字形及X字形的建筑群的最有利風(fēng)向角與建筑物側(cè)表面垂直,最不利風(fēng)向角則為斜風(fēng)向;而十字形和Y字形建筑群在不同布局下的結(jié)論與前者截然相反.

    圖16和圖17以方形及十字形建筑群為例,分析其作用機(jī)理,發(fā)現(xiàn)造成建筑群周邊區(qū)域風(fēng)速增大的兩個(gè)主要原因?yàn)榻遣糠蛛x效應(yīng)及狹管效應(yīng). 由于方形建筑前角的分離剪切作用較大,導(dǎo)致其角部分離效應(yīng)更加明顯,這就使得最不利風(fēng)向角出現(xiàn)于斜風(fēng)向下.

    通過計(jì)算不同工況下Rmax和A*? ? ? ? ? 1.2,avg值,從整體上對(duì)不同工況的行人風(fēng)環(huán)境進(jìn)行評(píng)估. 方形、H字形及X字形建筑群均在錯(cuò)列式布局下具有最優(yōu)的行人風(fēng)環(huán)境效應(yīng);十字形及Y字形建筑群則在中部開口式下具有最優(yōu)的行人風(fēng)環(huán)境效應(yīng). 綜合考慮建筑布局和建筑形狀,風(fēng)環(huán)境最好的是Y字形錯(cuò)列式布局建筑群,風(fēng)環(huán)境最差的是H字形圍合式布局建筑群.

    3? ?結(jié)? ?論

    通過風(fēng)洞試驗(yàn)與數(shù)值模擬方法,深入研究了建筑形狀和建筑布局對(duì)高層建筑群行人風(fēng)環(huán)境的影響,得出如下結(jié)論:

    1)經(jīng)過風(fēng)洞試驗(yàn)的驗(yàn)證,表明采用RANS模型研究高層建筑群周邊的行人風(fēng)環(huán)境是可行的,對(duì)于高風(fēng)速區(qū)域的模擬,通過修正湍流參數(shù),可以采用標(biāo)準(zhǔn)k-ε模型獲得與風(fēng)洞試驗(yàn)較為一致的模擬結(jié)果.

    2)造成建筑群周邊區(qū)域風(fēng)速增大的主要原因是建筑群中出現(xiàn)角部分離效應(yīng)和狹管效應(yīng). 建筑群周邊最大風(fēng)速比Rmax的大小不受建筑形狀以及建筑布局的影響,所有工況下Rmax均位于1.5附近.

    3)建筑形狀對(duì)建筑群的最不利風(fēng)向角影響較顯著,5種不同建筑形狀下,方形、H字形及X字形建筑群的最不利風(fēng)向角均出現(xiàn)在斜風(fēng)向角度,十字形和Y字形建筑群的最不利風(fēng)向?yàn)檎L(fēng)向;建筑布局對(duì)于建筑群的最不利風(fēng)向角度影響很小,4種不同布局下,建筑群最不利風(fēng)向角度未發(fā)生變化.

    4)綜合考慮建筑布局和建筑形狀的影響,行人風(fēng)環(huán)境最好的是Y字形錯(cuò)列式布局建筑群,最差的是H字形圍合式布局建筑群.

    5)在保持高度和容積率一定的情況下,建筑群周邊最大風(fēng)速比不會(huì)隨著建筑形狀和建筑布局的改變而發(fā)生明顯變化.

    參考文獻(xiàn)

    [1]? ? KUBOTA T,MIURA M,DE TOMINAGA Y,et al.? Wind tunnel tests on the relationship between building density and pedestrian-level wind velocity:Development of guidelines for realizing acceptable wind environment in residential neighborhoods[J]. Building and Environment,2008,43(10):1699—1708.

    [2]? ? 方智遠(yuǎn),汪之松,李正良. 雷暴沖擊風(fēng)作用下高層建筑風(fēng)荷載頻域特性[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2020,47(1):100—107.

    FANG Z Y,WANG Z S,LI Z L. Frequency domain characteristics of wind loads on high-rise buildings under thunderstorm downburst[J]. Journal of Hunan University (Natural Sciences),2020,47(1):100—107. (In Chinese)

    [3]? ? 汪夢(mèng)甫,梁曉婷. 帶伸臂桁架減振層高層結(jié)構(gòu)抗風(fēng)效果分析[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2018,45(3):1—7.

    WANG M F,LIANG X T. Analysis on wind resistance effect of high-rise building with damped outrigger storeys[J]. Journal of Hunan University (Natural Sciences),2018,45(3):1—7. (In Chinese)

    [4]? ? JUAN Y H,YANG A S,WEN C Y,et al. Optimization procedures for enhancement of city breathability using arcade design in a realistic high-rise urban area[J]. Building and Environment,2017,121:247—261.

    [5]? ? DU Y X,MAK C M. Improving pedestrian level low wind velocity environment in high-density cities:a general framework and case study[J]. Sustainable Cities and Society,2018,42:314—324.

    [6]? ? BLOCKEN B,JANSSEN W D,VAN HOOFF T. CFD simulation for pedestrian wind comfort and wind safety in urban areas:General decision framework and case study for the Eindhoven University campus[J]. Environmental Modelling & Software,2012,30:15—34.

    [7]? ? BLOCKEN B,STATHOPOULOS T,VAN BEECK J P A J. Pedestrian-level wind conditions around buildings:Review of wind-tunnel and CFD techniques and their accuracy for wind comfort assessment[J]. Building and Environment,2016,100:50—81.

    [8]? ? TSANG C W,KWOK K C S,HITCHCOCK P A. Wind tunnel study of pedestrian level wind environment around tall buildings:Effects of building dimensions,separation and podium[J]. Building and Environment,2012,49:167—181.

    [9]? ? XU X D,YANG Q S,YOSHIDA A,et al. Characteristics of pedestrian-level wind around super-tall buildings with various configurations[J]. Journal of Wind Engineering and Industrial Aerodynamics,2017,166:61—73.

    [10]? YOSHIE R,MOCHIDA A,TOMINAGA Y,et al. Cooperative project for CFD prediction of pedestrian wind environment in the Architectural Institute of Japan[J]. Journal of Wind Engineering and Industrial Aerodynamics,2007,95(9/11):1551—1578.

    [11]? MITTAL H,SHARMA A,GAIROLA A. A review on the study of urban wind at the pedestrian level around buildings[J]. Journal of Building Engineering,2018,18:154—163.

    [12]? BLOCKEN B. LES over RANS in building simulation for outdoor and indoor applications:a foregone conclusion?[J]. Building Simulation,2018,11(5):821—870.

    [13]? IQBAL Q M Y Z,CHAN A L S. Pedestrian level wind environment assessment around group of high-rise cross-shaped buildings:Effect of building shape,separation and orientation[J]. Building and Environment,2016,101:45—63.

    [14]? VAN DRUENEN T,VAN HOOFF T,MONTAZERI H,et al. CFD evaluation of building geometry modifications to reduce pedestrian-level wind speed[J]. Building and Environment,2019,163:106293.

    [15]? Recommendations for loads on buildings:AIJ—2004 [S]. Tykyo:Architectural Institute of Japan,2004:10—12.

    [16]? FRANKE J,HIRSCH C,JENSEN A G,et al. Recommendations on the use of CFD in wind engineering [C]//Proceedings of the International Conference on Urban Wind Engineering and Building Aerodynamics. Belgium,DC:Von Karman Institute,2004:568—576.

    [17]? FRANKE J,HELLSTEN A,SCHL?譈NZEN H,et al.? Best practice guideline for the CFD simulation of flows in the urban environment[M]. Brussels:Meteorological Inst. COST Office,2007:16—18.

    [18]? YANG Y,GU M,CHEN S Q,et al. New inflow boundary conditions for modelling the neutral equilibrium atmospheric boundary layer in computational wind engineering[J]. Journal of Wind Engineering and Industrial Aerodynamics,2009,97(2):88—95.

    [19]? BLOCKEN B,STATHOPOULOS T,CARMELIET J. CFD simulation of the atmospheric boundary layer:wall function problems[J]. Atmospheric Environment,2007,41(2):238—252.

    [20]? 羅凱文,楊易,謝壯寧. 基于k-ε模型模擬平衡態(tài)大氣邊界層的比較研究[J]. 工程力學(xué),2018,35(2):21—29.

    LUO K W,YANG Y,XIE Z N. A comparative study on the simulation of neutral atmospheric boundary layer based on the k-ε turbulence model[J]. Engineering Mechanics,2018,35(2):21—29. (In Chinese)

    [21] STATHOPOULOS T,WU H,B?魪DARD C. Wind environment around buildings:a knowledge-based approach[J]. Journal of Wind Engineering and Industrial Aerodynamics,1992,44(1/3):2377—2388.

    [22]? LIU Z Q,ISHIHARA T,TANAKA T,et al. LES study of turbulent flow fields over a smooth 3-D hill and a smooth 2-D ridge[J]. Journal of Wind Engineering and Industrial Aerodynamics,2016,153:1—12.

    [23]? 建筑工程風(fēng)洞試驗(yàn)方法標(biāo)準(zhǔn):JGJ/T 338—2014[S]. 北京:中國(guó)建筑工業(yè)出版社,2014:13—15

    Standard for wind tunnel test method of building engineering:JGJ/T 338—2014[S]. Beijing:China Architecture & Building? Press,2014:13—15. (In Chinese)

    [24]? 鄭朝榮,陳勇,金釗,等. 基于超越閾值概率的某千米級(jí)摩天大樓室外平臺(tái)行人風(fēng)環(huán)境評(píng)估[J]. 建筑結(jié)構(gòu)學(xué)報(bào),2018,39(2):122—129.

    ZHENG C R,CHEN Y,JIN Z,et al. Pedestrian wind environment assessment on the outdoor platform of a certain kilometer-level skyscraper based on the probability of exceeding the threshold [J]. Journal of Building Structures,2018,39(2):122—129. (In Chinese)

    [25]? BERNABEI R,BELLI P,MONTECCHIA F,et al. Possible implications of the channeling effect in NaI(Tl) crystals [J]. The European Physical Journal C,2008,53(2):205—213.

    [26]? 應(yīng)小宇,朱煒,外尾一則. 高層建筑群平面布局類型對(duì)室外風(fēng)環(huán)境影響的對(duì)比研究[J]. 地理科學(xué),2013,33(9):1097—1103.

    YING X Y,ZHU W,HOKAO K. Comparative study on the influence of layout types of high-rise buildings on outdoor wind environment[J]. Geographical Sciences,2013,33(9):1097—1103. (In Chinese)

    [27]? 馬劍,陳水福. 平面布局對(duì)高層建筑群風(fēng)環(huán)境影響的數(shù)值研究[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版),2007,41(9):1477—1481.

    MA J,CHEN S F. Numerical study of plan arrangement effect on wind environment around tall buildings[J].? Journal of Zhejiang University(Engineering Science),2007,41(9):1477—1481. (In Chinese)

    猜你喜歡
    風(fēng)洞
    綜合訓(xùn)練風(fēng)洞為科技奧運(yùn)助力
    斑頭雁進(jìn)風(fēng)洞
    小水滴在風(fēng)洞氣流中的跟隨性
    黃風(fēng)洞貂鼠精
    基于NI cRIO平臺(tái)的脈沖燃燒風(fēng)洞控制系統(tǒng)設(shè)計(jì)
    飛行器風(fēng)洞模型的快速制造技術(shù)
    19禁男女啪啪无遮挡网站| 美女高潮喷水抽搐中文字幕| 精品免费久久久久久久清纯| 国产成人av教育| 亚洲天堂国产精品一区在线| 侵犯人妻中文字幕一二三四区| www.自偷自拍.com| 国产一区二区激情短视频| 一区二区三区精品91| 99久久综合精品五月天人人| 日韩精品中文字幕看吧| 欧美av亚洲av综合av国产av| 久久久久久大精品| 国内精品久久久久久久电影| 一区二区三区国产精品乱码| 亚洲中文字幕一区二区三区有码在线看 | 久久天堂一区二区三区四区| 久久久久国内视频| 久久午夜综合久久蜜桃| 日韩中文字幕欧美一区二区| 亚洲久久久国产精品| 免费看十八禁软件| 最近最新中文字幕大全电影3 | 国内揄拍国产精品人妻在线 | 久久草成人影院| 久久这里只有精品19| 俺也久久电影网| 婷婷亚洲欧美| 天堂影院成人在线观看| 自线自在国产av| 一区二区三区高清视频在线| 999久久久精品免费观看国产| 两性夫妻黄色片| 免费在线观看黄色视频的| 成年女人毛片免费观看观看9| 日本三级黄在线观看| 高潮久久久久久久久久久不卡| 99热这里只有精品一区 | 中文字幕av电影在线播放| 看片在线看免费视频| 9191精品国产免费久久| 免费在线观看视频国产中文字幕亚洲| 久久精品国产综合久久久| 亚洲国产精品合色在线| 在线观看免费视频日本深夜| 亚洲中文字幕日韩| 999久久久国产精品视频| 亚洲自拍偷在线| 99re在线观看精品视频| 日韩欧美一区视频在线观看| 精华霜和精华液先用哪个| 草草在线视频免费看| 在线永久观看黄色视频| 日本免费a在线| 桃色一区二区三区在线观看| 成人永久免费在线观看视频| 日本五十路高清| 一个人免费在线观看的高清视频| 可以免费在线观看a视频的电影网站| 亚洲一区二区三区不卡视频| 午夜福利欧美成人| 99riav亚洲国产免费| 国产久久久一区二区三区| 老司机在亚洲福利影院| 亚洲成人国产一区在线观看| 日韩高清综合在线| 黑人欧美特级aaaaaa片| 一个人免费在线观看的高清视频| 首页视频小说图片口味搜索| 亚洲va日本ⅴa欧美va伊人久久| 久久国产乱子伦精品免费另类| av在线天堂中文字幕| 国产区一区二久久| 婷婷亚洲欧美| 欧美激情 高清一区二区三区| 丝袜人妻中文字幕| 免费在线观看日本一区| 最好的美女福利视频网| 欧美在线一区亚洲| www.熟女人妻精品国产| 精品国产超薄肉色丝袜足j| 亚洲精品色激情综合| 给我免费播放毛片高清在线观看| 国产精品野战在线观看| 村上凉子中文字幕在线| 99国产精品一区二区蜜桃av| 国产aⅴ精品一区二区三区波| 欧美性长视频在线观看| 村上凉子中文字幕在线| 欧美黑人巨大hd| 欧美中文综合在线视频| 成人特级黄色片久久久久久久| 国产精品免费一区二区三区在线| 久久九九热精品免费| 欧美久久黑人一区二区| 好男人电影高清在线观看| 色综合站精品国产| 最近在线观看免费完整版| 怎么达到女性高潮| 又紧又爽又黄一区二区| 久久人人精品亚洲av| 欧美日本视频| 桃色一区二区三区在线观看| 人人妻人人看人人澡| 亚洲午夜精品一区,二区,三区| 久久久久久免费高清国产稀缺| 两个人免费观看高清视频| 成人午夜高清在线视频 | 一级a爱视频在线免费观看| 99久久久亚洲精品蜜臀av| 国产伦一二天堂av在线观看| 中文字幕精品免费在线观看视频| 欧美三级亚洲精品| 美女午夜性视频免费| 日韩欧美一区视频在线观看| 亚洲国产欧洲综合997久久, | 青草久久国产| 别揉我奶头~嗯~啊~动态视频| 成人18禁在线播放| 亚洲精品国产一区二区精华液| 色老头精品视频在线观看| 男女午夜视频在线观看| 一进一出好大好爽视频| 国产高清有码在线观看视频 | 久久精品成人免费网站| 国产亚洲欧美精品永久| 中文在线观看免费www的网站 | 久热爱精品视频在线9| 亚洲三区欧美一区| 麻豆久久精品国产亚洲av| 亚洲成人国产一区在线观看| 99国产精品一区二区蜜桃av| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品久久男人天堂| 日本熟妇午夜| 欧美不卡视频在线免费观看 | 免费av毛片视频| 精品国产一区二区三区四区第35| 在线免费观看的www视频| 中国美女看黄片| 人妻久久中文字幕网| 日韩大尺度精品在线看网址| 夜夜看夜夜爽夜夜摸| 精品久久久久久,| 男女床上黄色一级片免费看| 国产亚洲精品久久久久5区| 麻豆国产av国片精品| 午夜激情福利司机影院| 国产成人av激情在线播放| 婷婷精品国产亚洲av| 久久久国产欧美日韩av| 国产精品一区二区三区四区久久 | 1024手机看黄色片| 一进一出好大好爽视频| 国产亚洲精品久久久久5区| 国产视频内射| 国产精品野战在线观看| 欧美亚洲日本最大视频资源| 国产精品,欧美在线| 制服丝袜大香蕉在线| 免费av毛片视频| 国产蜜桃级精品一区二区三区| 国产亚洲精品久久久久5区| 在线观看免费午夜福利视频| 成人午夜高清在线视频 | 99国产精品一区二区蜜桃av| 久久久久国内视频| 别揉我奶头~嗯~啊~动态视频| www国产在线视频色| 一级毛片女人18水好多| 国产视频一区二区在线看| 97碰自拍视频| 国产精品九九99| 又大又爽又粗| 性欧美人与动物交配| 少妇的丰满在线观看| 精品国产乱子伦一区二区三区| 看免费av毛片| 免费在线观看日本一区| 99精品在免费线老司机午夜| 午夜日韩欧美国产| 亚洲熟女毛片儿| 久久国产乱子伦精品免费另类| 亚洲精品国产精品久久久不卡| 黄片大片在线免费观看| 成熟少妇高潮喷水视频| 中文字幕av电影在线播放| 精品少妇一区二区三区视频日本电影| 欧美+亚洲+日韩+国产| 国产欧美日韩一区二区三| 欧美性猛交╳xxx乱大交人| 国产熟女午夜一区二区三区| 国产激情久久老熟女| 国产亚洲精品一区二区www| 黄色 视频免费看| 久久婷婷人人爽人人干人人爱| 免费女性裸体啪啪无遮挡网站| 热99re8久久精品国产| www.熟女人妻精品国产| 午夜日韩欧美国产| 亚洲精品一区av在线观看| 久久人妻福利社区极品人妻图片| xxxwww97欧美| 欧美乱色亚洲激情| 丁香欧美五月| 天堂√8在线中文| av天堂在线播放| 成人亚洲精品一区在线观看| 日韩有码中文字幕| 亚洲精品av麻豆狂野| 999久久久国产精品视频| 男女床上黄色一级片免费看| tocl精华| 黄片播放在线免费| 亚洲成人免费电影在线观看| 宅男免费午夜| 精品国产国语对白av| 美女免费视频网站| 久久99热这里只有精品18| 极品教师在线免费播放| 久久这里只有精品19| 搡老熟女国产l中国老女人| 亚洲av成人不卡在线观看播放网| 亚洲第一青青草原| 久久久久久久精品吃奶| 可以免费在线观看a视频的电影网站| 中文资源天堂在线| 男女午夜视频在线观看| 国产一区二区激情短视频| 国产熟女午夜一区二区三区| 夜夜躁狠狠躁天天躁| 国产99白浆流出| 亚洲精品美女久久久久99蜜臀| 这个男人来自地球电影免费观看| 熟女少妇亚洲综合色aaa.| 又黄又粗又硬又大视频| 村上凉子中文字幕在线| 黄色视频,在线免费观看| 热99re8久久精品国产| 国产在线精品亚洲第一网站| 999久久久精品免费观看国产| 免费在线观看完整版高清| 最近最新免费中文字幕在线| 最近最新中文字幕大全电影3 | 麻豆久久精品国产亚洲av| 免费av毛片视频| cao死你这个sao货| 亚洲一区高清亚洲精品| 美女国产高潮福利片在线看| 日韩欧美免费精品| 亚洲中文日韩欧美视频| 在线免费观看的www视频| 中出人妻视频一区二区| 精品人妻1区二区| 国产精华一区二区三区| 国产精品1区2区在线观看.| 日韩视频一区二区在线观看| 脱女人内裤的视频| 黄色女人牲交| 99热这里只有精品一区 | 久久精品国产综合久久久| 欧美激情 高清一区二区三区| www国产在线视频色| 日韩免费av在线播放| a级毛片在线看网站| 91在线观看av| 亚洲av成人av| 人人澡人人妻人| avwww免费| 精品一区二区三区四区五区乱码| 麻豆久久精品国产亚洲av| 日本 欧美在线| 国产伦一二天堂av在线观看| 日韩视频一区二区在线观看| 精品久久久久久久人妻蜜臀av| 在线永久观看黄色视频| 久久精品国产99精品国产亚洲性色| 精品国产亚洲在线| 这个男人来自地球电影免费观看| 午夜激情av网站| 国内毛片毛片毛片毛片毛片| 麻豆一二三区av精品| 欧美激情 高清一区二区三区| 成人三级做爰电影| 波多野结衣高清无吗| 亚洲五月婷婷丁香| 亚洲五月色婷婷综合| 亚洲欧美一区二区三区黑人| 欧美乱色亚洲激情| 99国产综合亚洲精品| 久久青草综合色| 亚洲国产看品久久| 老司机午夜福利在线观看视频| www.自偷自拍.com| 一进一出抽搐动态| 日韩中文字幕欧美一区二区| 国产精品久久视频播放| 婷婷精品国产亚洲av| 男女之事视频高清在线观看| 韩国av一区二区三区四区| 在线天堂中文资源库| 亚洲欧洲精品一区二区精品久久久| 亚洲国产毛片av蜜桃av| 成年免费大片在线观看| 神马国产精品三级电影在线观看 | 精品国产乱子伦一区二区三区| 少妇被粗大的猛进出69影院| 欧美另类亚洲清纯唯美| 国产精品影院久久| 国内精品久久久久久久电影| 91字幕亚洲| 欧美色视频一区免费| 美女国产高潮福利片在线看| 亚洲精品美女久久久久99蜜臀| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 日韩三级视频一区二区三区| 午夜福利视频1000在线观看| 一本精品99久久精品77| 一区二区三区高清视频在线| 国产一区二区三区在线臀色熟女| x7x7x7水蜜桃| netflix在线观看网站| 中文字幕人成人乱码亚洲影| 欧美成人一区二区免费高清观看 | 大型av网站在线播放| 天堂影院成人在线观看| 91字幕亚洲| 最近最新中文字幕大全免费视频| 国产av一区在线观看免费| ponron亚洲| 亚洲人成77777在线视频| 色av中文字幕| 国产av在哪里看| 久9热在线精品视频| 一区二区三区激情视频| 高清毛片免费观看视频网站| 一级a爱视频在线免费观看| 最近最新中文字幕大全免费视频| 又黄又爽又免费观看的视频| 免费人成视频x8x8入口观看| 欧美激情久久久久久爽电影| 成人午夜高清在线视频 | 一本久久中文字幕| 亚洲专区国产一区二区| 日韩大码丰满熟妇| 国产亚洲精品久久久久久毛片| 免费人成视频x8x8入口观看| 国产精品一区二区精品视频观看| 中亚洲国语对白在线视频| 成人国产综合亚洲| 亚洲国产精品成人综合色| 色婷婷久久久亚洲欧美| 欧美日韩乱码在线| 超碰成人久久| 欧美激情高清一区二区三区| 欧美国产精品va在线观看不卡| 久久久精品欧美日韩精品| 麻豆一二三区av精品| 国产午夜精品久久久久久| 脱女人内裤的视频| 成人免费观看视频高清| 热re99久久国产66热| 黄色视频不卡| 亚洲国产欧美网| 一个人免费在线观看的高清视频| 久久久久免费精品人妻一区二区 | 国内精品久久久久精免费| 免费在线观看完整版高清| 中文字幕最新亚洲高清| 99精品欧美一区二区三区四区| 中文在线观看免费www的网站 | 桃红色精品国产亚洲av| 不卡一级毛片| 欧美日韩福利视频一区二区| 日韩精品青青久久久久久| 午夜激情福利司机影院| 91成人精品电影| 午夜精品久久久久久毛片777| 欧美亚洲日本最大视频资源| 麻豆成人午夜福利视频| 亚洲av成人不卡在线观看播放网| 久久久精品国产亚洲av高清涩受| 欧美日本亚洲视频在线播放| 性欧美人与动物交配| 亚洲国产毛片av蜜桃av| 亚洲男人的天堂狠狠| 欧美黑人精品巨大| 亚洲欧洲精品一区二区精品久久久| 免费av毛片视频| 午夜久久久久精精品| 国产91精品成人一区二区三区| 国产人伦9x9x在线观看| 亚洲久久久国产精品| 我的亚洲天堂| 国产午夜精品久久久久久| 久久国产亚洲av麻豆专区| 久久久国产精品麻豆| 99久久99久久久精品蜜桃| 亚洲av电影在线进入| 中文资源天堂在线| 满18在线观看网站| 黄色成人免费大全| 禁无遮挡网站| ponron亚洲| 国产成人一区二区三区免费视频网站| 村上凉子中文字幕在线| 国产亚洲欧美精品永久| 1024香蕉在线观看| 欧美一级a爱片免费观看看 | 99久久精品国产亚洲精品| 2021天堂中文幕一二区在线观 | 国产私拍福利视频在线观看| 一级毛片女人18水好多| 精品久久久久久久毛片微露脸| 国产成人精品久久二区二区免费| 国产精品av久久久久免费| 国产伦在线观看视频一区| 无遮挡黄片免费观看| 日韩视频一区二区在线观看| 中文字幕久久专区| 国产蜜桃级精品一区二区三区| 亚洲免费av在线视频| 999久久久国产精品视频| 免费人成视频x8x8入口观看| 美国免费a级毛片| 日韩三级视频一区二区三区| 日本一本二区三区精品| 欧美黑人精品巨大| 少妇的丰满在线观看| 久久久久亚洲av毛片大全| 欧美色视频一区免费| 亚洲欧美精品综合一区二区三区| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 久久中文看片网| 精品午夜福利视频在线观看一区| 亚洲国产高清在线一区二区三 | 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品第一综合不卡| 午夜激情福利司机影院| 丝袜在线中文字幕| 在线观看免费日韩欧美大片| 亚洲欧美精品综合一区二区三区| 中国美女看黄片| 色哟哟哟哟哟哟| 无人区码免费观看不卡| 国产成人系列免费观看| 国产av一区在线观看免费| 国产精品一区二区三区四区久久 | 亚洲精品色激情综合| 久久国产亚洲av麻豆专区| 香蕉丝袜av| 国产乱人伦免费视频| 嫁个100分男人电影在线观看| 丝袜在线中文字幕| 夜夜夜夜夜久久久久| 国产又爽黄色视频| 精品第一国产精品| a级毛片a级免费在线| 欧美绝顶高潮抽搐喷水| 中文字幕另类日韩欧美亚洲嫩草| 国产乱人伦免费视频| 成人手机av| www.自偷自拍.com| 成人亚洲精品av一区二区| 国产91精品成人一区二区三区| 99精品在免费线老司机午夜| 免费看美女性在线毛片视频| 国产精品一区二区精品视频观看| 香蕉国产在线看| 精品一区二区三区视频在线观看免费| 手机成人av网站| 狠狠狠狠99中文字幕| 国产亚洲精品综合一区在线观看 | 午夜福利免费观看在线| 高潮久久久久久久久久久不卡| 国产一卡二卡三卡精品| 亚洲精品国产一区二区精华液| 欧美av亚洲av综合av国产av| 午夜激情福利司机影院| 18禁观看日本| 久久久国产成人免费| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲最大成人中文| 99国产综合亚洲精品| 在线国产一区二区在线| 亚洲中文av在线| 看免费av毛片| 久久精品影院6| 欧美一级a爱片免费观看看 | 90打野战视频偷拍视频| 国产精品野战在线观看| 日韩三级视频一区二区三区| 人妻丰满熟妇av一区二区三区| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 免费av毛片视频| 高清在线国产一区| 他把我摸到了高潮在线观看| 精品久久久久久久久久久久久 | 欧美乱妇无乱码| 婷婷亚洲欧美| 亚洲一码二码三码区别大吗| 亚洲自偷自拍图片 自拍| 久久中文字幕人妻熟女| 国产亚洲欧美在线一区二区| 成人18禁高潮啪啪吃奶动态图| 他把我摸到了高潮在线观看| 国产片内射在线| 国产免费av片在线观看野外av| 亚洲国产精品久久男人天堂| 国产精品久久久人人做人人爽| 欧美中文综合在线视频| 国产主播在线观看一区二区| 日本三级黄在线观看| 久久99热这里只有精品18| 亚洲avbb在线观看| 久久精品aⅴ一区二区三区四区| 国产精品久久电影中文字幕| 国产高清激情床上av| 亚洲国产精品久久男人天堂| 欧美午夜高清在线| 亚洲人成网站高清观看| 国产欧美日韩精品亚洲av| 国产精品野战在线观看| 91成年电影在线观看| 又黄又粗又硬又大视频| 欧美zozozo另类| 一二三四在线观看免费中文在| www.www免费av| 久久欧美精品欧美久久欧美| 国产成人啪精品午夜网站| 亚洲熟妇熟女久久| 在线十欧美十亚洲十日本专区| 最近最新免费中文字幕在线| 男人舔奶头视频| 国产av一区在线观看免费| 午夜福利免费观看在线| 侵犯人妻中文字幕一二三四区| 非洲黑人性xxxx精品又粗又长| 国产男靠女视频免费网站| 欧美中文日本在线观看视频| 亚洲国产看品久久| 久久午夜综合久久蜜桃| 亚洲av美国av| 可以在线观看毛片的网站| 91成年电影在线观看| 又黄又粗又硬又大视频| 69av精品久久久久久| 岛国在线观看网站| 男女视频在线观看网站免费 | 亚洲三区欧美一区| 欧美成狂野欧美在线观看| 12—13女人毛片做爰片一| 99久久国产精品久久久| 免费看a级黄色片| 久久久久国产一级毛片高清牌| 国产不卡一卡二| 美女扒开内裤让男人捅视频| 久久精品国产清高在天天线| 我的亚洲天堂| 欧美性猛交黑人性爽| 国产熟女xx| 51午夜福利影视在线观看| 老熟妇乱子伦视频在线观看| 欧美成人一区二区免费高清观看 | 特大巨黑吊av在线直播 | 免费在线观看成人毛片| 国产又爽黄色视频| 正在播放国产对白刺激| 亚洲av成人av| 中文资源天堂在线| 亚洲aⅴ乱码一区二区在线播放 | 老汉色∧v一级毛片| 久久这里只有精品19| 国产真人三级小视频在线观看| 国产片内射在线| 久久精品影院6| 搡老熟女国产l中国老女人| 99国产精品一区二区蜜桃av| 亚洲人成网站在线播放欧美日韩| 国产熟女午夜一区二区三区| 亚洲精品一区av在线观看| 精品久久久久久成人av| 岛国视频午夜一区免费看| 在线国产一区二区在线| 亚洲av熟女| 亚洲一区二区三区色噜噜| 亚洲免费av在线视频| 精品国产美女av久久久久小说| 免费观看精品视频网站| 亚洲人成77777在线视频| 亚洲一区中文字幕在线| 久久性视频一级片| 亚洲国产看品久久| 色综合婷婷激情| 国产色视频综合| 一区福利在线观看| 国产亚洲精品综合一区在线观看 | 午夜免费鲁丝| 老鸭窝网址在线观看| 亚洲av成人不卡在线观看播放网| 久久热在线av| 日日干狠狠操夜夜爽| 大型黄色视频在线免费观看| 最新美女视频免费是黄的| 日韩视频一区二区在线观看| 国产99白浆流出| 日韩欧美在线二视频| 免费看美女性在线毛片视频| 老熟妇仑乱视频hdxx| 丰满人妻熟妇乱又伦精品不卡| 国内毛片毛片毛片毛片毛片|