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

    基于改進(jìn)大渦模擬模型的風(fēng)荷載驗(yàn)證

    2016-10-20 05:03:07婁伶伶吳玖榮
    關(guān)鍵詞:大渦風(fēng)洞試驗(yàn)風(fēng)壓

    婁伶伶, 吳玖榮

    (廣州大學(xué) 廣州大學(xué)-淡江大學(xué) 工程結(jié)構(gòu)災(zāi)害與控制聯(lián)合研究中心, 廣東 廣州 510006)

    ?

    基于改進(jìn)大渦模擬模型的風(fēng)荷載驗(yàn)證

    婁伶伶, 吳玖榮*

    (廣州大學(xué) 廣州大學(xué)-淡江大學(xué) 工程結(jié)構(gòu)災(zāi)害與控制聯(lián)合研究中心, 廣東 廣州510006)

    為了獲得與近地大氣邊界層實(shí)際脈動(dòng)風(fēng)場(chǎng)特征相一致的三維脈動(dòng)風(fēng)速,文章首先應(yīng)用DSRFG (Discretizing and Synthesizing of Random Flow Generation)方法,模擬生成計(jì)算流體動(dòng)力學(xué)中大渦模擬算法所需的模擬入口邊界脈動(dòng)風(fēng)速時(shí)程,然后采用改進(jìn)的動(dòng)態(tài)一方程模型作為大渦模擬的亞格子模型,計(jì)算長(zhǎng)寬高為1∶1∶3的高層單體矩形建筑在0度風(fēng)向角下的風(fēng)荷載分布及特性,并與日本TPU風(fēng)洞試驗(yàn)進(jìn)行了對(duì)比.結(jié)果表明:①當(dāng)高層建筑長(zhǎng)寬比(B/D)=1時(shí),扭矩的能量主要來自渦旋脫落而不是入口湍流的影響;②DSRFG結(jié)合改進(jìn)的動(dòng)態(tài)一方程模型能夠較為準(zhǔn)確地模擬得到高層建筑表面的平均及脈動(dòng)風(fēng)荷載特征.此方法有望應(yīng)用到實(shí)際高層建筑風(fēng)荷載評(píng)估中.

    高層建筑; DSRFG; 大渦模擬; 亞格子模型; 風(fēng)荷載; 風(fēng)洞試驗(yàn)

    高層建筑具有高、細(xì)、柔的特征,其自振頻率較低且阻尼較小,過大的風(fēng)速會(huì)使得高層建筑外圍結(jié)構(gòu)甚至主體結(jié)構(gòu)損傷.由于能夠獲得風(fēng)場(chǎng)作用下風(fēng)荷載的脈動(dòng)時(shí)程及相關(guān)統(tǒng)計(jì)信息,大渦模擬(Large Eddy Simulation, LES)已經(jīng)成為計(jì)算流體動(dòng)力學(xué)(Computational Fluid Dynamics, CFD)對(duì)高層建筑風(fēng)荷載數(shù)值模擬研究的主要手段之一.能否較好實(shí)現(xiàn)LES在風(fēng)工程中的應(yīng)用,TAMURA[1]總結(jié)了3個(gè)關(guān)鍵點(diǎn):①入口湍流風(fēng)速的準(zhǔn)確模擬;②亞格子湍流模型的選用及處理;③計(jì)算流體動(dòng)力學(xué)方法中相關(guān)數(shù)值離散化算法的處理.

    近地邊界層風(fēng)為具有高雷諾數(shù)的非定常流,且表現(xiàn)出較強(qiáng)的不均勻性和各向異性.因此,為獲得與近地邊界層實(shí)際脈動(dòng)風(fēng)場(chǎng)特征相一致的三維脈動(dòng)風(fēng)速,需要較為準(zhǔn)確地模擬生成滿足相空間相關(guān)性和湍流特性的入口邊界脈動(dòng)風(fēng)速時(shí)程.目前CFD中主要采用2種入口邊界脈動(dòng)風(fēng)速生成方法:①在輔助的流體計(jì)算域內(nèi)通過recycling方法生成湍流風(fēng);②通過人工數(shù)值模型生成入口邊界風(fēng).NOZAWA等[2]用第一種方法生成了入口邊界.DSRFG[3]屬于第二種方法,DSRFG方法可以生成滿足特定頻譜和空間相關(guān)性的脈動(dòng)湍流場(chǎng),且生成的脈動(dòng)速度滿足連續(xù)和自保持條件.

    大渦模擬算法中采用的亞格子模型一般有如下幾種方式:代數(shù)渦粘模型[4-5]、一方程模型[6]、兩方程模型[7].然而上述亞格子模型在計(jì)算風(fēng)工程實(shí)際應(yīng)用中存在2個(gè)有待改進(jìn)之處:①大多數(shù)亞格子模型適用于高精度數(shù)值格式和結(jié)構(gòu)網(wǎng)格[8],但是實(shí)際工程中結(jié)構(gòu)比較復(fù)雜,一般采用精度相對(duì)較低的數(shù)值分析算法及無結(jié)構(gòu)網(wǎng)格;②在高雷諾數(shù)流和網(wǎng)格相對(duì)粗糙情況時(shí),代數(shù)渦粘模型的局部平衡假定與其表現(xiàn)出來的不平衡特征相矛盾.因此,在結(jié)合上述亞格子模型各自優(yōu)缺點(diǎn)的基礎(chǔ)上,有必要建立一種與實(shí)際工程CFD風(fēng)荷載模擬需求情況較為接近的改進(jìn)亞格子模型.HUANG[9]等采用結(jié)合動(dòng)態(tài)一方程模型[6]和WALE模型[5]各自優(yōu)點(diǎn)的改進(jìn)亞格子模型,進(jìn)行大渦模擬高層建筑風(fēng)荷載,其特點(diǎn)有:①可用于高雷諾數(shù)流和粗糙網(wǎng)格的計(jì)算;②確定動(dòng)態(tài)參數(shù)時(shí)不需要試驗(yàn)過濾,可用于低階數(shù)值離散和無結(jié)構(gòu)化網(wǎng)格;③用不同的機(jī)制處理亞格子能量的生成及GS(grid-scale)能量耗散.

    本文采用DSRFG方法,模擬與高層建筑所處的近地大氣邊界層風(fēng)場(chǎng)特征較為一致的入口邊界風(fēng)速時(shí)程,并采用改進(jìn)的動(dòng)態(tài)一方程模型[9]相對(duì)應(yīng)的大渦模擬算法,計(jì)算了長(zhǎng)寬高比為1∶1∶3的高層單體矩形建筑在0度風(fēng)向角下的風(fēng)荷載分布及特性,并與日本TPU風(fēng)洞試驗(yàn)[10]結(jié)果進(jìn)行了對(duì)比研究,分析LES獲得的平均及脈動(dòng)風(fēng)荷載特性,驗(yàn)證了上述方法的有效性和準(zhǔn)確性.

    1 本文采用的大渦模擬相關(guān)算法概述

    本文的大渦模擬CFD分析是在含有48個(gè)CPU的多處理器工作站進(jìn)行的,該工作站可用于大規(guī)模的并行計(jì)算任務(wù),采用計(jì)算軟件為Flunet 6.3.

    1.1計(jì)算模型與網(wǎng)格劃分

    采用的高層建筑模型[10]尺寸為40 m×40 m×120 m,幾何縮尺比1∶400.CFD計(jì)算域尺寸:X、Y、Z3個(gè)方向?yàn)? 600 m×1 000 m×1 000 m.CFD計(jì)算時(shí)縮小模型風(fēng)洞阻塞比不應(yīng)超過3%[11],本文模型的計(jì)算阻塞比為3%,滿足模擬要求.

    CFD計(jì)算時(shí)采用六面體結(jié)構(gòu)化網(wǎng)格,圖1為計(jì)算域網(wǎng)格的劃分情況.為保證模擬精度并節(jié)省計(jì)算時(shí)間,對(duì)建筑物表面和模型附近的網(wǎng)格進(jìn)行加密(圖2),并向外逐漸增大網(wǎng)格密度,劃分出約70萬個(gè)網(wǎng)格單元.

    圖1 CFD計(jì)算區(qū)域示意圖

    圖2 模擬建筑物周邊網(wǎng)格局部放大示意圖

    1.2CFD中入口邊界脈動(dòng)風(fēng)速時(shí)程模擬

    風(fēng)工程中大多數(shù)實(shí)測(cè)的近地大氣邊界層脈動(dòng)風(fēng)速譜都服從馮·卡門(von Karman)譜[12],為此本文采用DSRFG方法[3,13],以模擬生成滿足目標(biāo)譜(馮·卡門譜)的CFD入口邊界風(fēng)速時(shí)程.其主要步驟如下:

    (1)對(duì)于具有任意形式的3D脈動(dòng)風(fēng)速能量譜E(k),均可由一系列離散的E(km)構(gòu)成:

    (1)

    對(duì)于任一方向(此處假設(shè)為x方向)具有離散風(fēng)速能量譜的脈動(dòng)風(fēng)速,其脈動(dòng)速度場(chǎng)可由下式模擬:

    (2)

    (2)脈動(dòng)風(fēng)速在頻域的合成:

    (3)

    (3.1)

    (3.2)

    ωm,n∈N(0,2πfm),fm=kmUavg

    (3.3)

    其中,ζ,ξ∈N(0,1),fm為頻率,Uavg為平均風(fēng)速,Ls為湍流積分尺度.a滿足(0~1)之間的均勻分布.kmax=500,N=50時(shí)模擬結(jié)果與目標(biāo)譜較為接近[12].因此,本文中kmax=500,N=50.

    本文采用B類地貌,α=0.15.入口邊界處平均風(fēng)速度剖面采用如下的指數(shù)率形式:

    (4)

    其中,Z0=10 m,U0代表10 m高度處的平均風(fēng)速,本文取15 m·s-1.

    根據(jù)日本《AIJ建筑荷載建議》[14],湍流強(qiáng)度表達(dá)式如下:

    (5)

    其中,zb為近地高度尺度,取7.5 m,zG代表梯度風(fēng)高度取350 m.圖3、圖4分別為入口處平均風(fēng)速、X方向湍流強(qiáng)度隨高度變化示意圖.

    圖3 平均風(fēng)剖面

    圖4 X向湍流強(qiáng)度

    圖5給出了DSRFG法模擬的風(fēng)場(chǎng)中的脈動(dòng)風(fēng)速功率譜與馮卡門譜的對(duì)比圖,可看出2者吻合較好,說明入口湍流脈動(dòng)風(fēng)可以較為真實(shí)地模擬實(shí)際大氣湍流邊界層情況[3,13].

    圖5 DSRFG法生成的脈動(dòng)風(fēng)速功率譜與馮·卡門譜對(duì)比

    1.3大渦模擬中改進(jìn)亞格子模型概述

    基于KAJISHIMA[6]的一方程模型,ksgs可以從方程(6)求得[9]:

    (6)

    (7)

    最后一項(xiàng)εω是額外的耗散項(xiàng),可以解釋近壁面低雷諾數(shù)的影響.

    根冠比是指植物地下部分與地上部分的鮮重或干重的比值,它的大小能反映出植株地上、地下部生長(zhǎng)發(fā)育情況以及土壤的營(yíng)養(yǎng)供應(yīng)狀況,該指標(biāo)高則根系機(jī)能活性強(qiáng),低則弱[22]。錦紫蘇在不同栽培基質(zhì)下,除處理⑥外,其他處理根冠比均高于對(duì)照組,說明不同配比基質(zhì)的城市堆肥污泥有利于植物根系發(fā)達(dá),增加植物的根系機(jī)能活性,從而有利于培養(yǎng)壯株,有利于作物的地上部分生長(zhǎng)。試驗(yàn)證明,城市污泥堆肥對(duì)于錦紫蘇生長(zhǎng)具有好的適應(yīng)性,并為污泥堆肥直接用于觀葉類草本植物的種植提供了理論基礎(chǔ),在一定程度上促進(jìn)了城市污泥堆肥的資源化利用。

    (8)

    本文所采用的改進(jìn)亞格子模型具有如下優(yōu)點(diǎn):不需要試驗(yàn)過濾,可用于無結(jié)構(gòu)網(wǎng)格計(jì)算,可用于高雷諾數(shù)流的模擬,采用Fluent UDF編程加以完成.

    1.4CFD求解設(shè)置及邊界設(shè)置

    CFD中有限體積法對(duì)應(yīng)的離散方程組的求解,一般采用收斂性較好的SIMPLEC算法[15].對(duì)一定偏斜的網(wǎng)格,Skewness Correction可減小收斂困難,默認(rèn)為0.LES計(jì)算中,動(dòng)量方程對(duì)流項(xiàng)的離散采用Bounded Central Differencing算法.時(shí)間項(xiàng)的離散采用二階隱式方案.各項(xiàng)松弛因子均采用默認(rèn)值.計(jì)算時(shí)間步長(zhǎng)0.001 s,CFD模擬計(jì)算采用10 000步.

    流場(chǎng)出口采用完全發(fā)展的出口邊界條件,計(jì)算域頂部和兩側(cè)采用自由滑移壁面,計(jì)算域底面和高層建筑表面均采用不可滑移的壁面.

    2 CFD數(shù)值模擬結(jié)果

    日本TPU風(fēng)洞試驗(yàn)包含了具有多種不同長(zhǎng)寬高比值的矩形單體建筑模型,位于各類常用的地貌條件下的多測(cè)點(diǎn)高頻測(cè)壓風(fēng)洞試驗(yàn)結(jié)果,為高層建筑結(jié)構(gòu)抗風(fēng)設(shè)計(jì)提供了在多風(fēng)向角條件下的結(jié)構(gòu)物表面詳細(xì)的風(fēng)荷載時(shí)程風(fēng)洞試驗(yàn)數(shù)據(jù)[10].本文CFD模擬結(jié)果將與TPU數(shù)據(jù)庫(kù)中長(zhǎng)寬高之比為1∶1∶3的矩形單體高層建筑模型,在0度風(fēng)向角下的風(fēng)洞試驗(yàn)數(shù)據(jù)結(jié)果進(jìn)行對(duì)比分析,以分析本文所采用大渦模擬CFD計(jì)算算法的有效性.

    2.1平均風(fēng)壓系數(shù)

    以本文所分析的高層建筑模型在Z=0.25 m、Z=0.15 m、Z=0.05 m 3個(gè)高度的部分測(cè)點(diǎn)為研究對(duì)象,各層布置的測(cè)點(diǎn)數(shù)與風(fēng)洞試驗(yàn)相同[10],為20個(gè),測(cè)點(diǎn)布置見圖6.

    圖6 測(cè)點(diǎn)布置圖

    (9)

    圖7為大渦模擬與風(fēng)洞試驗(yàn)在模型3個(gè)高度處不同測(cè)點(diǎn)平均風(fēng)壓系數(shù)的對(duì)比圖.1~5號(hào)測(cè)點(diǎn)在高層建筑迎風(fēng)面,6~10號(hào)測(cè)點(diǎn)在右側(cè)面,11~15號(hào)測(cè)點(diǎn)在背風(fēng)面,16~20號(hào)測(cè)點(diǎn)在左側(cè)面.

    圖7風(fēng)洞試驗(yàn)與數(shù)值模擬各測(cè)點(diǎn)的平均風(fēng)壓系數(shù)對(duì)比

    Fig.7Comparison of mean wind pressure coefficient by CFD and wind tunnel

    通過對(duì)比分析可知:①在模型Z=0.25 m高度,迎風(fēng)面為正壓,呈現(xiàn)出中間大、邊緣小(最大0.89,最小0.49)的變化趨勢(shì).5號(hào)測(cè)點(diǎn)與6號(hào)測(cè)點(diǎn)之間風(fēng)壓系數(shù)迅速降低,說明來流遇到建筑物的阻擋在模型兩側(cè)分離較強(qiáng),邊緣處負(fù)壓梯度較大.側(cè)風(fēng)面與背風(fēng)面均為負(fù)壓,且各測(cè)點(diǎn)所得的模擬結(jié)果精度比較高,特別是10號(hào)測(cè)點(diǎn)至17號(hào)測(cè)點(diǎn),誤差均在10%以內(nèi).迎風(fēng)面誤差較大,可能是迎風(fēng)面與湍流入口距離不足或邊界條件設(shè)置的影響所導(dǎo)致,但是整體趨勢(shì)吻合較好;②當(dāng)來流平行于結(jié)構(gòu)對(duì)稱軸時(shí),CFD數(shù)值模擬的平均風(fēng)壓基本為對(duì)稱分布,而風(fēng)洞試驗(yàn)數(shù)據(jù)對(duì)稱性較差;③Z=0.15 m、Z=0.05 m時(shí)與Z=0.25 m時(shí)平均風(fēng)壓系數(shù)的分布類似.

    2.2脈動(dòng)風(fēng)壓均方根系數(shù)

    脈動(dòng)風(fēng)壓系數(shù)反映脈動(dòng)能量的大小是脈動(dòng)風(fēng)荷載的重要特征.脈動(dòng)風(fēng)壓均方根系數(shù)如下:

    (10)

    圖8可見:①脈動(dòng)風(fēng)壓均方根系數(shù)與平均風(fēng)壓系數(shù)具有相似的分布趨勢(shì);②來流在迎風(fēng)面邊緣分離再附著,導(dǎo)致迎風(fēng)面脈動(dòng)風(fēng)壓值較大,變化梯度也相對(duì)較大,在其他部位變化相對(duì)平緩;③隨著高度的變化,建筑背風(fēng)面脈動(dòng)風(fēng)壓均方根系數(shù)差值在10%以內(nèi),說明背風(fēng)面渦旋的作用比較均勻;④脈動(dòng)風(fēng)壓系數(shù)隨著高度的降低而逐步變小,說明建筑低處含有的脈動(dòng)能量相對(duì)較小.

    圖8 數(shù)值模擬各測(cè)點(diǎn)的脈動(dòng)風(fēng)壓均方根系數(shù)

    Fig.8Root mean square coefficient of fluctuating wind pressure by CFD

    2.3層脈動(dòng)風(fēng)荷載功率譜

    高層建筑風(fēng)荷載分為2個(gè)主軸方向上的順風(fēng)向、橫風(fēng)向荷載以及繞Z軸的扭矩.沿高度將建筑分為6層,各測(cè)點(diǎn)與風(fēng)洞試驗(yàn)測(cè)點(diǎn)布置相對(duì)應(yīng).將CFD計(jì)算時(shí)各層測(cè)點(diǎn)每個(gè)時(shí)間步上的脈動(dòng)風(fēng)壓,進(jìn)行積分可得到各層的脈動(dòng)風(fēng)荷載時(shí)程數(shù)據(jù),然后從頻域的角度表達(dá)層風(fēng)荷載的特性.限于篇幅,本文只取第2、4、6層進(jìn)行對(duì)比分析.

    2.3.1層順風(fēng)向力功率譜

    ①圖9a、b、c分別為第2、4、6層的順風(fēng)荷載歸一化功率譜,低頻段(無量綱頻率0.04~0.40)能量分布較均勻,隨后譜密度開始下降,但是衰減速度比風(fēng)洞試驗(yàn)結(jié)果慢,可能是因?yàn)镃FD模擬中入口處與模型迎風(fēng)面的距離比風(fēng)洞試驗(yàn)?。虎谀芰恐饕性诘皖l段,且與圖5所示的入口脈動(dòng)風(fēng)速功率譜形狀相似,說明入口風(fēng)湍流是產(chǎn)生順風(fēng)向風(fēng)荷載的主要原因;③圖9d為前述3層數(shù)值模擬功率譜對(duì)比,在低頻段第6層功率譜值最小,第2層最大,因此,低頻段較高層的功率譜能量比低層小,可能是脈動(dòng)風(fēng)在建筑迎風(fēng)面下沉,在地面形成回流,湍流較大導(dǎo)致的.高頻段不同層上沒有明顯的變化.

    2.3.2層橫風(fēng)向力功率譜

    ①圖10a、b、c分別為第2、4、6層的橫風(fēng)向荷載歸一化功率譜,低頻段(無量綱頻率0.04~0.10)能量高于風(fēng)洞試驗(yàn)結(jié)果.高頻段與風(fēng)洞試驗(yàn)結(jié)果較為接近;②不同層功率譜曲線均在無量綱頻率0.1處出現(xiàn)一個(gè)明顯的尖峰,其主要是渦旋脫離造成的;③由圖10d知,低頻段第2層峰值最大,第6層峰值最小,高頻段3者數(shù)值基本接近.

    2.3.2層扭轉(zhuǎn)風(fēng)荷載功率譜

    ①圖11a、b、c分別為第2、4、6層的扭轉(zhuǎn)風(fēng)荷載歸一化功率譜,當(dāng)頻率小于無量綱頻率1時(shí)模擬值結(jié)果大于風(fēng)洞試驗(yàn)值,大于1時(shí)模擬值小于風(fēng)洞實(shí)驗(yàn)值,在無量綱頻率0.10附近,出現(xiàn)一個(gè)明顯的尖峰,與橫風(fēng)向功率譜峰值對(duì)應(yīng)的頻率接近,渦旋脫落現(xiàn)象明顯;②扭轉(zhuǎn)譜與風(fēng)洞試驗(yàn)結(jié)果的誤差,與層順風(fēng)向及層橫風(fēng)向譜和相應(yīng)風(fēng)洞試驗(yàn)結(jié)果對(duì)比來看,說明本文CFD模擬結(jié)果對(duì)層扭轉(zhuǎn)譜的模擬精度要相對(duì)低一些.但CFD數(shù)值模擬結(jié)果所反映的扭轉(zhuǎn)譜的渦旋脫落現(xiàn)象,與風(fēng)洞試驗(yàn)結(jié)果基本相似.

    圖9 層順風(fēng)向荷載功率譜

    圖10 層橫風(fēng)向荷載功率譜

    圖11 層扭轉(zhuǎn)風(fēng)荷載功率譜

    3 結(jié) 論

    本文應(yīng)用DSRFG方法,模擬生成CFD大渦模擬算法所需的入口邊界脈動(dòng)風(fēng)速時(shí)程,采用改進(jìn)的動(dòng)態(tài)一方程模型作為大渦模擬的亞格子模型,計(jì)算了長(zhǎng)寬高為1∶1∶3的高層單體矩形建筑在0度風(fēng)向角下的風(fēng)荷載分布及特性,并與日本TPU風(fēng)洞試驗(yàn)進(jìn)行了對(duì)比分析,得到以下結(jié)論:

    (1)DSRFG可以生成大渦模擬算法中與近地大氣邊界層實(shí)際脈動(dòng)風(fēng)場(chǎng)特征相一致的入口邊界脈動(dòng)風(fēng)速場(chǎng).DSRFG方法結(jié)合改進(jìn)的動(dòng)態(tài)一方程亞格子模型,運(yùn)用于CFD大渦模擬算法中,可以較好地預(yù)測(cè)高層建筑表面的風(fēng)荷載,數(shù)值模擬與風(fēng)洞試驗(yàn)數(shù)據(jù)結(jié)果基本一致.

    (2)LES模擬得到的高層建筑表面平均風(fēng)壓系數(shù)與風(fēng)洞試驗(yàn)數(shù)據(jù)吻合較好,反映了典型鈍體繞流的特征;來流平行于結(jié)構(gòu)對(duì)稱軸時(shí),風(fēng)洞試驗(yàn)的平均風(fēng)壓對(duì)稱性很差,而數(shù)值模擬的對(duì)稱性很好,體現(xiàn)了數(shù)值模擬的優(yōu)勢(shì).

    (3)LES的脈動(dòng)風(fēng)壓均方根系數(shù)與平均風(fēng)壓系數(shù)具有相似的分布趨勢(shì).隨著高度的增加,迎風(fēng)面及兩側(cè)面脈動(dòng)能量有逐步增大的趨勢(shì),背風(fēng)面渦旋作用較均勻,沒有明顯變化.

    (4)LES對(duì)層風(fēng)荷載歸一化功率譜的模擬表明,入口湍流風(fēng)是產(chǎn)生順風(fēng)向荷載的主要原因;渦旋脫落是造成橫風(fēng)向荷載的主要原因,功率譜曲線在無量綱頻率0.1處出現(xiàn)明顯的尖峰;扭矩一般由入口風(fēng)湍流與渦旋脫落相互作用造成,但對(duì)于本文的模型B/D=1,扭矩主要來自渦旋脫落的影響.

    [1]TAMURA T. Towards practical use of LES in wind engineering[J]. J Wind Eng Indust Aerod, 2008, 96:1451-1471.

    [2]NOZAWA K, TAMURA T. Large eddy simulation of wind flows over large roughness elements[C]∥Proceedings of the 4th European and African Conference on Wind Engineering (EACWE4),2005, Prague, Cezch Republic.

    [3]HUANG S H, LI Q S, WU J R. A general inflow turbulence generator for large eddy simulation[J]. J Wind Eng Indust Aerod, 2010, 98(10/11):600-617.

    [4]GERMANO M, PIOMELLI U, MOIN P, et al. A dynamic subgrid-scale eddy viscosity model[J]. Phys Fluid, 1991, 3(7):1760-1765.

    [5]NICOUD F, DUCROS F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor[J]. Flow Turbul Combust, 1999, 62(3):183-200.

    [6]KAJISHIMA T, NOMACHI T, KAJISHIMA T, et al. One-equation subgrid scale model using dynamic procedure for the energy production[J]. J Appl Mech, 2005, 73(3):368-373.

    [7]GALLERANO F, PASERO E, CANNATA G. A dynamic two-equation sub-grid scale model[J]. Contin Mech Thermod, 2005, 17(2): 101-123.

    [8]TSUBOKURA M. Proper representation of the subgrid-scale eddy viscosity for the dynamic procedure in large eddy simulation using finite difference method[J]. Phys Fluid, 2001, 13(2):500-504.

    [9]HUANG S, LI Q S. A new dynamic one-equation subgrid-scale model for large eddy simulations[J]. Intern J Numer Method Eng, 2009, 81(7):835-865.

    [10]林寧,梁波,田村幸雄. 高層建筑層風(fēng)力特性實(shí)驗(yàn)研究[J]. 振動(dòng)工程學(xué)報(bào),2003(12): 409-414.

    LIN N, LIANG B, YUKIO T. Experimental investigation on local wind force of high-r ise building[J]. J Vibr Eng,2003(12):409-414.

    [11]黃本才. 結(jié)構(gòu)抗風(fēng)分析原理及應(yīng)用[M].上海:同濟(jì)大學(xué)出版社, 2008.

    HUANG B C. Principle and application of structural wind resistance analysis[M]. Shanghai: Tongji University Press, 2008.

    [12]LI Q S, WU J R, LIANG S G, et al. Full-scale measurements and numerical evaluation of wind-induced vibration of a 63-story reinforced concrete tall building[J]. Eng Struct, 2004, 26(12):1779-1794.

    [13]熊濤, 吳玖榮. 大渦模擬中入流邊界脈動(dòng)風(fēng)速模擬相關(guān)參數(shù)取值分析[J]. 廣州大學(xué)學(xué)報(bào):自然科學(xué)版, 2014, 13(3):60-67.

    XIONG T, WU J R. Parameter analysis for modeling the inlet fluctuated wind velocity in large eddy simulation[J]. J Guangzhou Univ:Nat Sci Edi,2014,13(3):60-67.

    [14]Architectural Institute of Japan. Recommendations for loads on buildings[M]. (English Version). Tokyo: AIJ Press, 2004.

    [15]VAN DOORMAAL J P, RAITHBY G D. Enhancements of the simple method for predicting incompressible fluid flows[J]. Numer Heat Transf Appl, 1984, 7(2):147-163.

    【責(zé)任編輯: 周全】

    Verification and evaluation on wind loading of large eddy simulation by the improved sub-grid scale model

    LOU Ling-ling, WU Jiu-rong

    (Guangzhou University-Tamkang University Joint Research Center for Engineering Structure Disaster Prevention and Control, Guangzhou University, Guangzhou 510006, China)

    In order to obtain the time history of 3-dimensional fluctuating wind speed with the same statistical characteristics as the wind of near ground atmosphere boundary layer, an inflow turbulence generation method called DSRFG (discretizing and synthesizing of random flow generation) is applied to generate the time record of inflow turbulence wind speed for Large Eddy Simulation in CFD, meanwhile an improved dynamic one-equation model is adopted as the sub-grid scale model for Large Eddy Simulation method. These two methods are combined for the CFD simulation on the wind load distribution and characteristics of an isolated rectangular tall building with 1∶1∶3 ratio (the ratio between length, width and height) in zero degree wind direction. The CFD numerical simulation results are further compared with those experimental wind tunnel test results on the same building from TPU wind tunnel laboratory in Japan. All these comparisons indicate that: the fluctuating kinetic energy of wind-induced torsional moment mainly results from the influence of vortex shedding other than the inlet turbulence for the test building; the combination of DSRFG with the improved dynamic one-equation model in the large eddy simulation method in CFD can accurately obtain the value and their characteristics of mean and fluctuating wind loads acting on the isolated rectangular tall building, which jointly testify the feasibility of its application to the design wind load evaluation for actual high rise buildings.

    tall building; DSRFG; large eddy simulation; sub-grid scale model; wind load; wind tunnel test

    2016-02-22;

    2016-04-12

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51378134,51578169);廣東省高層次人才資助項(xiàng)目

    婁伶伶(1989-),女,碩士研究生.E-mail:loulinghcool@163.com

    . E-mail: jrwuce@qq.com

    1671- 4229(2016)04-0050-07

    TU 312.1

    A

    猜你喜歡
    大渦風(fēng)洞試驗(yàn)風(fēng)壓
    基于壁面射流的下?lián)舯┝鞣欠€(wěn)態(tài)風(fēng)場(chǎng)大渦模擬
    軸流風(fēng)機(jī)葉尖泄漏流動(dòng)的大渦模擬
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    滾轉(zhuǎn)機(jī)動(dòng)載荷減緩風(fēng)洞試驗(yàn)
    低風(fēng)壓導(dǎo)線的發(fā)展和測(cè)試
    電線電纜(2016年5期)2016-02-27 09:02:06
    基于大渦模擬的旋風(fēng)分離器錐體結(jié)構(gòu)影響研究
    遮擋條件下超高層建筑風(fēng)洞試驗(yàn)研究
    重慶建筑(2014年12期)2014-07-24 14:00:32
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    高層建筑風(fēng)荷載有哪些計(jì)算要求
    重慶建筑(2014年10期)2014-03-27 02:20:07
    二1-21061機(jī)巷低抽巷全風(fēng)壓通風(fēng)應(yīng)用方案研究
    河南科技(2014年10期)2014-02-27 14:09:23
    国产区一区二久久| 亚洲精品粉嫩美女一区| 动漫黄色视频在线观看| 色综合婷婷激情| 黄色丝袜av网址大全| 真人一进一出gif抽搐免费| 日韩三级视频一区二区三区| 精品少妇一区二区三区视频日本电影| 欧美另类亚洲清纯唯美| 欧美黑人精品巨大| 国产精品1区2区在线观看.| 亚洲第一欧美日韩一区二区三区| 麻豆国产av国片精品| 国产精品av久久久久免费| 草草在线视频免费看| 韩国av一区二区三区四区| av片东京热男人的天堂| 一级片免费观看大全| 哪里可以看免费的av片| 国产免费av片在线观看野外av| 手机成人av网站| 欧美中文综合在线视频| 精品熟女少妇八av免费久了| 成在线人永久免费视频| 国内揄拍国产精品人妻在线 | 亚洲精品av麻豆狂野| 三级毛片av免费| a级毛片在线看网站| 国产精品久久久久久人妻精品电影| 免费在线观看完整版高清| 久久人人精品亚洲av| 999久久久精品免费观看国产| 亚洲一区二区三区不卡视频| 中文字幕最新亚洲高清| 操出白浆在线播放| 中国美女看黄片| а√天堂www在线а√下载| 国产精品爽爽va在线观看网站 | 一a级毛片在线观看| 国内少妇人妻偷人精品xxx网站 | 中亚洲国语对白在线视频| 欧美成人一区二区免费高清观看 | 国产乱人伦免费视频| 成年版毛片免费区| 手机成人av网站| 亚洲国产欧美日韩在线播放| 欧美 亚洲 国产 日韩一| 窝窝影院91人妻| 中文在线观看免费www的网站 | 丝袜人妻中文字幕| 免费看十八禁软件| 免费看十八禁软件| 国产黄色小视频在线观看| 欧美性猛交黑人性爽| 亚洲片人在线观看| 久久亚洲精品不卡| 男人舔女人的私密视频| 国产精品乱码一区二三区的特点| 精品一区二区三区av网在线观看| 久久久久久久久免费视频了| 好男人在线观看高清免费视频 | 午夜免费鲁丝| 国产av一区在线观看免费| 日本撒尿小便嘘嘘汇集6| 国产伦人伦偷精品视频| 长腿黑丝高跟| 午夜福利18| 一进一出抽搐gif免费好疼| 亚洲国产精品合色在线| 夜夜躁狠狠躁天天躁| 午夜福利成人在线免费观看| 1024手机看黄色片| 18美女黄网站色大片免费观看| 久久中文字幕一级| 成人亚洲精品av一区二区| 亚洲九九香蕉| 日韩欧美三级三区| 国产成年人精品一区二区| 亚洲真实伦在线观看| 午夜福利高清视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲第一欧美日韩一区二区三区| 国内精品久久久久久久电影| 欧美日韩福利视频一区二区| 久久香蕉激情| 精品熟女少妇八av免费久了| 母亲3免费完整高清在线观看| 国产精品爽爽va在线观看网站 | 午夜福利免费观看在线| 长腿黑丝高跟| 精品不卡国产一区二区三区| 一级毛片女人18水好多| 午夜免费鲁丝| 免费在线观看黄色视频的| 一本一本综合久久| 麻豆成人午夜福利视频| 一本久久中文字幕| 午夜福利在线观看吧| av视频在线观看入口| 久久人妻福利社区极品人妻图片| 免费女性裸体啪啪无遮挡网站| 亚洲av日韩精品久久久久久密| 午夜亚洲福利在线播放| 日韩av在线大香蕉| 男女视频在线观看网站免费 | 国产精品综合久久久久久久免费| 亚洲精品一卡2卡三卡4卡5卡| 老熟妇乱子伦视频在线观看| 国产免费男女视频| 99久久99久久久精品蜜桃| 国产人伦9x9x在线观看| 丰满人妻熟妇乱又伦精品不卡| 狂野欧美激情性xxxx| 99热6这里只有精品| 在线av久久热| 十分钟在线观看高清视频www| 国产97色在线日韩免费| 日本在线视频免费播放| 色av中文字幕| svipshipincom国产片| 久久久久精品国产欧美久久久| 国产精品 欧美亚洲| 国产精品 欧美亚洲| 午夜成年电影在线免费观看| 国产激情久久老熟女| 淫妇啪啪啪对白视频| 亚洲av五月六月丁香网| 国产成人欧美| 欧美成人午夜精品| 免费无遮挡裸体视频| 久久精品亚洲精品国产色婷小说| 国产精品99久久99久久久不卡| 欧美亚洲日本最大视频资源| 久久精品91蜜桃| 欧美又色又爽又黄视频| 丝袜在线中文字幕| 精品国产亚洲在线| 国产一级毛片七仙女欲春2 | 中文字幕人妻丝袜一区二区| 成人三级黄色视频| 免费看a级黄色片| 久久国产精品男人的天堂亚洲| 亚洲,欧美精品.| 色综合站精品国产| 亚洲激情在线av| 美女高潮喷水抽搐中文字幕| 国产av不卡久久| 国产亚洲精品久久久久5区| 国产精品 国内视频| 好男人在线观看高清免费视频 | 少妇熟女aⅴ在线视频| 亚洲 国产 在线| or卡值多少钱| 99精品久久久久人妻精品| 成人18禁在线播放| 午夜免费鲁丝| 久久精品亚洲精品国产色婷小说| www日本在线高清视频| 成年版毛片免费区| 免费在线观看影片大全网站| 999久久久国产精品视频| 亚洲自拍偷在线| 国产高清激情床上av| 99精品欧美一区二区三区四区| АⅤ资源中文在线天堂| 俺也久久电影网| 丰满人妻熟妇乱又伦精品不卡| 高潮久久久久久久久久久不卡| 看免费av毛片| 午夜福利免费观看在线| 亚洲精品国产精品久久久不卡| 色精品久久人妻99蜜桃| 免费一级毛片在线播放高清视频| 国产在线观看jvid| 麻豆成人午夜福利视频| 国产精品免费视频内射| 女性生殖器流出的白浆| 欧美中文综合在线视频| 12—13女人毛片做爰片一| 极品教师在线免费播放| 欧美久久黑人一区二区| 色在线成人网| 别揉我奶头~嗯~啊~动态视频| 美国免费a级毛片| 日韩欧美一区视频在线观看| 日韩一卡2卡3卡4卡2021年| 国产精品国产高清国产av| 欧美激情 高清一区二区三区| 又大又爽又粗| 欧美 亚洲 国产 日韩一| 国产一卡二卡三卡精品| 国产单亲对白刺激| 每晚都被弄得嗷嗷叫到高潮| 一级片免费观看大全| 91老司机精品| 成人三级黄色视频| bbb黄色大片| 成年版毛片免费区| 久久久精品欧美日韩精品| 美女扒开内裤让男人捅视频| 一进一出抽搐gif免费好疼| 亚洲七黄色美女视频| 国产精品一区二区免费欧美| 国产不卡一卡二| 1024手机看黄色片| 国语自产精品视频在线第100页| 88av欧美| 亚洲avbb在线观看| 亚洲精品一区av在线观看| 国产成人影院久久av| 国产成人精品久久二区二区免费| 国产精品免费一区二区三区在线| 欧美一级毛片孕妇| 99国产极品粉嫩在线观看| 中亚洲国语对白在线视频| 热re99久久国产66热| 日韩精品免费视频一区二区三区| 一本久久中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久人人做人人爽| 国产成人一区二区三区免费视频网站| 国产亚洲精品久久久久久毛片| 国产精品精品国产色婷婷| 成人18禁高潮啪啪吃奶动态图| 久久人人精品亚洲av| a在线观看视频网站| 一本综合久久免费| 久久久久久大精品| 亚洲 欧美一区二区三区| avwww免费| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产精品影院| 男女之事视频高清在线观看| 亚洲av成人一区二区三| 一进一出好大好爽视频| 久久草成人影院| 一进一出抽搐gif免费好疼| 亚洲一码二码三码区别大吗| 亚洲午夜理论影院| 国产午夜精品久久久久久| 婷婷丁香在线五月| 人人妻,人人澡人人爽秒播| 黄片大片在线免费观看| 黄色成人免费大全| 国产成人精品无人区| 脱女人内裤的视频| 宅男免费午夜| 国产精品久久电影中文字幕| 丁香欧美五月| 别揉我奶头~嗯~啊~动态视频| 欧美不卡视频在线免费观看 | 国产欧美日韩一区二区精品| 中文亚洲av片在线观看爽| 国产黄色小视频在线观看| 窝窝影院91人妻| 丰满的人妻完整版| 国产成人啪精品午夜网站| 日韩免费av在线播放| 久久精品91无色码中文字幕| 亚洲一区二区三区不卡视频| 黄色视频不卡| 91老司机精品| 国产精品免费一区二区三区在线| 色尼玛亚洲综合影院| 精品久久久久久,| 一本一本综合久久| 国产黄片美女视频| 精品国产美女av久久久久小说| 在线观看免费日韩欧美大片| 91老司机精品| 婷婷精品国产亚洲av| 美女午夜性视频免费| 久久久久久久精品吃奶| 婷婷丁香在线五月| 一区二区三区激情视频| 黄色视频不卡| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 国产黄色小视频在线观看| 黑人欧美特级aaaaaa片| 97碰自拍视频| cao死你这个sao货| 黄片播放在线免费| 999久久久国产精品视频| 免费女性裸体啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| 国产av在哪里看| 给我免费播放毛片高清在线观看| 亚洲成人免费电影在线观看| 久久精品91无色码中文字幕| 国语自产精品视频在线第100页| 国内毛片毛片毛片毛片毛片| 久久久久免费精品人妻一区二区 | 精品欧美一区二区三区在线| 国产亚洲精品第一综合不卡| 三级毛片av免费| 久久精品91蜜桃| 亚洲五月天丁香| 少妇熟女aⅴ在线视频| 黑丝袜美女国产一区| 精品国产乱子伦一区二区三区| 在线观看一区二区三区| 老司机午夜福利在线观看视频| 国产精品久久久久久亚洲av鲁大| 亚洲人成电影免费在线| 国产爱豆传媒在线观看 | 亚洲欧洲精品一区二区精品久久久| 久久久国产欧美日韩av| 18禁黄网站禁片午夜丰满| 夜夜看夜夜爽夜夜摸| 久久久久久久精品吃奶| 曰老女人黄片| cao死你这个sao货| 老熟妇乱子伦视频在线观看| 叶爱在线成人免费视频播放| 一本综合久久免费| 极品教师在线免费播放| 国产精品亚洲美女久久久| 国产精品爽爽va在线观看网站 | 1024香蕉在线观看| 最新美女视频免费是黄的| 两个人免费观看高清视频| www国产在线视频色| 国产精品久久久久久亚洲av鲁大| 非洲黑人性xxxx精品又粗又长| 不卡一级毛片| 欧美成狂野欧美在线观看| 国产亚洲精品久久久久久毛片| 亚洲五月天丁香| 成人国产综合亚洲| 亚洲aⅴ乱码一区二区在线播放 | 亚洲片人在线观看| 国产单亲对白刺激| 亚洲第一av免费看| 桃色一区二区三区在线观看| 天天躁夜夜躁狠狠躁躁| 一进一出抽搐动态| 男人舔女人的私密视频| 国产国语露脸激情在线看| 老熟妇仑乱视频hdxx| 国产色视频综合| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一小说| 女同久久另类99精品国产91| 男人的好看免费观看在线视频 | 欧美黄色片欧美黄色片| 日韩欧美在线二视频| 午夜福利欧美成人| 日本三级黄在线观看| 日韩 欧美 亚洲 中文字幕| 中出人妻视频一区二区| 国产又爽黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 一区福利在线观看| 久久久久精品国产欧美久久久| av在线播放免费不卡| 精品日产1卡2卡| 午夜福利18| 首页视频小说图片口味搜索| 欧美在线黄色| 国产久久久一区二区三区| 久久狼人影院| 午夜成年电影在线免费观看| 亚洲激情在线av| 午夜久久久久精精品| 午夜免费观看网址| 我的亚洲天堂| 琪琪午夜伦伦电影理论片6080| 国产1区2区3区精品| 哪里可以看免费的av片| 亚洲av中文字字幕乱码综合 | 美女国产高潮福利片在线看| 国产精品久久视频播放| 亚洲五月色婷婷综合| 日韩高清综合在线| 亚洲精品在线美女| 免费在线观看黄色视频的| 男人舔女人的私密视频| 成年免费大片在线观看| 亚洲欧美日韩无卡精品| 午夜久久久在线观看| 男女下面进入的视频免费午夜 | 最好的美女福利视频网| 夜夜夜夜夜久久久久| 亚洲真实伦在线观看| 国产97色在线日韩免费| 国产又色又爽无遮挡免费看| 校园春色视频在线观看| 日韩欧美一区视频在线观看| 久久热在线av| 老汉色av国产亚洲站长工具| 久久 成人 亚洲| 黄色视频不卡| 香蕉久久夜色| 满18在线观看网站| 国产精品一区二区精品视频观看| 久久精品成人免费网站| 欧美三级亚洲精品| 国产真实乱freesex| 久久精品人妻少妇| e午夜精品久久久久久久| 日韩成人在线观看一区二区三区| 国产激情久久老熟女| 一区二区日韩欧美中文字幕| 日本一本二区三区精品| 国产一区二区激情短视频| 国产亚洲精品综合一区在线观看 | 亚洲免费av在线视频| 欧美日本亚洲视频在线播放| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| 50天的宝宝边吃奶边哭怎么回事| 欧美性猛交黑人性爽| 久久久久免费精品人妻一区二区 | 人人妻人人看人人澡| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品合色在线| 午夜免费激情av| 黄色 视频免费看| 热99re8久久精品国产| 村上凉子中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 国产蜜桃级精品一区二区三区| www.熟女人妻精品国产| 亚洲成av片中文字幕在线观看| 亚洲精品av麻豆狂野| 日韩视频一区二区在线观看| 欧美 亚洲 国产 日韩一| 免费搜索国产男女视频| 亚洲中文字幕日韩| 欧美日韩一级在线毛片| 国产亚洲精品第一综合不卡| 亚洲国产精品久久男人天堂| 亚洲久久久国产精品| 中文字幕av电影在线播放| 欧美性猛交╳xxx乱大交人| 国产精品亚洲美女久久久| cao死你这个sao货| 久久精品国产亚洲av香蕉五月| 国产精品久久电影中文字幕| 深夜精品福利| 精品高清国产在线一区| 欧美成人性av电影在线观看| 久久国产精品影院| 亚洲人成电影免费在线| 亚洲欧美日韩无卡精品| 中文字幕最新亚洲高清| 久久久水蜜桃国产精品网| 亚洲午夜理论影院| 夜夜爽天天搞| 亚洲国产欧美日韩在线播放| 欧美黄色淫秽网站| 婷婷六月久久综合丁香| 欧美一级a爱片免费观看看 | 国产成人欧美在线观看| 免费av毛片视频| 我的亚洲天堂| av在线天堂中文字幕| 手机成人av网站| 久久性视频一级片| 好男人电影高清在线观看| 桃色一区二区三区在线观看| 黄色女人牲交| 啦啦啦观看免费观看视频高清| 女生性感内裤真人,穿戴方法视频| 女人爽到高潮嗷嗷叫在线视频| 最近最新中文字幕大全免费视频| 久久精品91蜜桃| 久久久精品国产亚洲av高清涩受| 啦啦啦免费观看视频1| www.精华液| 久久久久久久久久黄片| 久久精品国产综合久久久| 免费人成视频x8x8入口观看| 久久精品91蜜桃| 99热6这里只有精品| 日韩 欧美 亚洲 中文字幕| 九色国产91popny在线| 久久久国产精品麻豆| 女性被躁到高潮视频| 亚洲精品久久国产高清桃花| 欧美一区二区精品小视频在线| 日本撒尿小便嘘嘘汇集6| 一区福利在线观看| 亚洲 欧美一区二区三区| 欧美日韩乱码在线| 黄网站色视频无遮挡免费观看| 午夜激情福利司机影院| 亚洲国产看品久久| 黄色a级毛片大全视频| 欧美性猛交╳xxx乱大交人| 欧美精品亚洲一区二区| 国产极品粉嫩免费观看在线| 一a级毛片在线观看| 午夜免费鲁丝| 一夜夜www| 国产精品亚洲一级av第二区| 精品一区二区三区四区五区乱码| 久久狼人影院| 国产精品电影一区二区三区| 两个人视频免费观看高清| 国产亚洲精品久久久久久毛片| 人人妻人人澡人人看| 日韩成人在线观看一区二区三区| 久久热在线av| 日韩欧美免费精品| 日本黄色视频三级网站网址| 久久久久久亚洲精品国产蜜桃av| 久久中文看片网| 久久精品国产综合久久久| 精品久久久久久久末码| 此物有八面人人有两片| 午夜久久久在线观看| 国产午夜精品久久久久久| 女性被躁到高潮视频| 真人一进一出gif抽搐免费| 最好的美女福利视频网| 国产又黄又爽又无遮挡在线| 欧美性猛交╳xxx乱大交人| 精品不卡国产一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产蜜桃级精品一区二区三区| 真人做人爱边吃奶动态| a级毛片在线看网站| 大型黄色视频在线免费观看| 熟女少妇亚洲综合色aaa.| www.自偷自拍.com| 在线播放国产精品三级| 成人亚洲精品一区在线观看| 亚洲专区国产一区二区| 很黄的视频免费| 国产高清激情床上av| 麻豆一二三区av精品| 一个人观看的视频www高清免费观看 | 亚洲性夜色夜夜综合| av免费在线观看网站| 中文字幕人妻熟女乱码| 欧美在线黄色| 日日摸夜夜添夜夜添小说| 99久久无色码亚洲精品果冻| 国产高清视频在线播放一区| 国产激情久久老熟女| 琪琪午夜伦伦电影理论片6080| 午夜激情av网站| 日本一区二区免费在线视频| 国产av一区二区精品久久| 精品高清国产在线一区| 又大又爽又粗| 波多野结衣av一区二区av| 九色国产91popny在线| a级毛片a级免费在线| 一二三四在线观看免费中文在| 非洲黑人性xxxx精品又粗又长| 色综合亚洲欧美另类图片| 黑人操中国人逼视频| 嫁个100分男人电影在线观看| av天堂在线播放| 国产片内射在线| 久久中文字幕人妻熟女| 国产一区二区三区视频了| 成年女人毛片免费观看观看9| 欧美黄色片欧美黄色片| 色婷婷久久久亚洲欧美| 亚洲人成电影免费在线| 国产av一区二区精品久久| 久久精品国产综合久久久| 亚洲五月婷婷丁香| 亚洲自拍偷在线| 日韩大尺度精品在线看网址| 满18在线观看网站| 好男人电影高清在线观看| 国产人伦9x9x在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产成人欧美| 国产精品 欧美亚洲| 亚洲性夜色夜夜综合| 亚洲中文字幕一区二区三区有码在线看 | 久久人妻福利社区极品人妻图片| 亚洲国产高清在线一区二区三 | 午夜福利在线在线| 国产精品一区二区免费欧美| 成人午夜高清在线视频 | 国产成人啪精品午夜网站| 亚洲av电影在线进入| 精品一区二区三区视频在线观看免费| 成年女人毛片免费观看观看9| 啦啦啦 在线观看视频| 老司机在亚洲福利影院| 国产视频内射| 欧美黑人欧美精品刺激| 成年版毛片免费区| av片东京热男人的天堂| 岛国视频午夜一区免费看| 欧美激情久久久久久爽电影| 非洲黑人性xxxx精品又粗又长| 国产精品,欧美在线| 搞女人的毛片| 国产麻豆成人av免费视频| 国产精品,欧美在线| 又黄又爽又免费观看的视频| 亚洲人成伊人成综合网2020| 日本精品一区二区三区蜜桃| 757午夜福利合集在线观看| 久久久精品国产亚洲av高清涩受| 久久久久久久久久黄片| 成年版毛片免费区| 成年免费大片在线观看| 在线天堂中文资源库| 国产精品免费视频内射| 狂野欧美激情性xxxx| 视频区欧美日本亚洲| 国产午夜福利久久久久久| 法律面前人人平等表现在哪些方面| 俄罗斯特黄特色一大片| 久久久久久人人人人人|