• <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
    两人在一起打扑克的视频| 热re99久久国产66热| 国产又爽黄色视频| 怎么达到女性高潮| 午夜精品在线福利| 国产91精品成人一区二区三区| 精品国产美女av久久久久小说| 涩涩av久久男人的天堂| 色综合婷婷激情| 超碰成人久久| 亚洲欧美日韩高清在线视频| 少妇猛男粗大的猛烈进出视频| 日韩免费高清中文字幕av| 亚洲欧美激情综合另类| 欧美一级毛片孕妇| 丰满迷人的少妇在线观看| 男男h啪啪无遮挡| 国产精品久久久久久精品古装| ponron亚洲| 深夜精品福利| 久久香蕉激情| 美女福利国产在线| 十八禁高潮呻吟视频| 女性生殖器流出的白浆| aaaaa片日本免费| 脱女人内裤的视频| 757午夜福利合集在线观看| 国产亚洲精品第一综合不卡| 久久人妻熟女aⅴ| 日本vs欧美在线观看视频| 久久 成人 亚洲| 亚洲精品国产区一区二| 亚洲成av片中文字幕在线观看| 一边摸一边抽搐一进一小说 | 亚洲精品国产区一区二| 国产一区二区激情短视频| 欧美日韩视频精品一区| 免费一级毛片在线播放高清视频 | 下体分泌物呈黄色| 丰满的人妻完整版| 高清毛片免费观看视频网站 | 国产精品久久久久久人妻精品电影| 久久久久久亚洲精品国产蜜桃av| 欧美乱色亚洲激情| 免费黄频网站在线观看国产| 在线看a的网站| 久久久久国产一级毛片高清牌| 成人三级做爰电影| 在线播放国产精品三级| 国产麻豆69| 老司机午夜福利在线观看视频| 18禁黄网站禁片午夜丰满| 欧美黄色淫秽网站| svipshipincom国产片| 亚洲精品久久午夜乱码| 中文字幕精品免费在线观看视频| 亚洲五月色婷婷综合| 日韩制服丝袜自拍偷拍| 午夜福利乱码中文字幕| 久久久水蜜桃国产精品网| 18禁美女被吸乳视频| 后天国语完整版免费观看| 无人区码免费观看不卡| 两个人看的免费小视频| 一二三四在线观看免费中文在| tube8黄色片| 亚洲精品国产精品久久久不卡| 99国产综合亚洲精品| 美女高潮到喷水免费观看| 9热在线视频观看99| netflix在线观看网站| 侵犯人妻中文字幕一二三四区| 欧美亚洲日本最大视频资源| 99久久人妻综合| 精品国产一区二区三区四区第35| 久久精品aⅴ一区二区三区四区| 人妻一区二区av| 久久中文字幕人妻熟女| 视频在线观看一区二区三区| 欧美日韩一级在线毛片| 老司机福利观看| 啦啦啦在线免费观看视频4| 女人爽到高潮嗷嗷叫在线视频| 国产精品一区二区精品视频观看| 国产精品综合久久久久久久免费 | 国产91精品成人一区二区三区| 久久人妻福利社区极品人妻图片| 久久九九热精品免费| 一进一出抽搐动态| 亚洲精品乱久久久久久| 久久影院123| 丝袜美足系列| 变态另类成人亚洲欧美熟女 | 别揉我奶头~嗯~啊~动态视频| 婷婷丁香在线五月| 欧美久久黑人一区二区| 韩国av一区二区三区四区| 女同久久另类99精品国产91| 乱人伦中国视频| 91麻豆精品激情在线观看国产 | 天天添夜夜摸| 国产精品偷伦视频观看了| 女人被狂操c到高潮| 黑人欧美特级aaaaaa片| 国产精品乱码一区二三区的特点 | 久久香蕉激情| 在线观看舔阴道视频| 丝瓜视频免费看黄片| 久久久水蜜桃国产精品网| 欧美日韩成人在线一区二区| 亚洲成国产人片在线观看| 久久草成人影院| 一级片'在线观看视频| 久久人人爽av亚洲精品天堂| 精品高清国产在线一区| 欧美 亚洲 国产 日韩一| 少妇 在线观看| 精品电影一区二区在线| 日韩精品免费视频一区二区三区| 露出奶头的视频| 久久亚洲真实| 夜夜夜夜夜久久久久| 日韩欧美在线二视频 | 天堂动漫精品| 大陆偷拍与自拍| 如日韩欧美国产精品一区二区三区| 精品午夜福利视频在线观看一区| 丝袜美足系列| 伊人久久大香线蕉亚洲五| 亚洲人成77777在线视频| 日韩一卡2卡3卡4卡2021年| 女人爽到高潮嗷嗷叫在线视频| 丁香欧美五月| 青草久久国产| 久久精品亚洲精品国产色婷小说| 成年动漫av网址| 高清毛片免费观看视频网站 | av电影中文网址| 日韩 欧美 亚洲 中文字幕| 午夜福利在线观看吧| 久久国产精品人妻蜜桃| videos熟女内射| а√天堂www在线а√下载 | 欧美日韩乱码在线| 成人18禁高潮啪啪吃奶动态图| 纯流量卡能插随身wifi吗| 国产成人精品在线电影| √禁漫天堂资源中文www| 国产欧美日韩一区二区三| 日本a在线网址| 午夜久久久在线观看| 女同久久另类99精品国产91| 亚洲av成人av| 国产精品亚洲av一区麻豆| 又黄又爽又免费观看的视频| 国产av一区二区精品久久| 国产精品一区二区在线观看99| 极品少妇高潮喷水抽搐| 脱女人内裤的视频| 99久久综合精品五月天人人| 亚洲欧美日韩高清在线视频| 国产乱人伦免费视频| 亚洲av成人不卡在线观看播放网| 91麻豆精品激情在线观看国产 | 久久香蕉国产精品| 久久性视频一级片| 国产亚洲欧美精品永久| 一区二区三区国产精品乱码| 一区二区三区激情视频| 亚洲片人在线观看| 日本a在线网址| 久久这里只有精品19| 看黄色毛片网站| 老司机靠b影院| 久久精品国产清高在天天线| 天堂√8在线中文| 国产淫语在线视频| 免费人成视频x8x8入口观看| 黄色丝袜av网址大全| a级毛片黄视频| 国产在线观看jvid| 五月开心婷婷网| 久9热在线精品视频| 9色porny在线观看| 美女国产高潮福利片在线看| av中文乱码字幕在线| 国产午夜精品久久久久久| 又大又爽又粗| 一级片'在线观看视频| 美国免费a级毛片| 天天躁夜夜躁狠狠躁躁| 午夜免费观看网址| 亚洲一码二码三码区别大吗| 国产麻豆69| 国精品久久久久久国模美| 一二三四在线观看免费中文在| 亚洲第一av免费看| 国产亚洲精品久久久久5区| 18禁黄网站禁片午夜丰满| 国产99久久九九免费精品| 极品少妇高潮喷水抽搐| 久久午夜亚洲精品久久| 手机成人av网站| 欧美老熟妇乱子伦牲交| 很黄的视频免费| 九色亚洲精品在线播放| 黄色视频不卡| 国产片内射在线| 91字幕亚洲| 亚洲av熟女| 91大片在线观看| 日韩制服丝袜自拍偷拍| 国产精华一区二区三区| 亚洲专区字幕在线| 一进一出好大好爽视频| 少妇的丰满在线观看| 亚洲三区欧美一区| cao死你这个sao货| 欧美日韩国产mv在线观看视频| 亚洲视频免费观看视频| 中文字幕色久视频| 18禁观看日本| 麻豆av在线久日| 男男h啪啪无遮挡| 久久久精品免费免费高清| 亚洲一卡2卡3卡4卡5卡精品中文| 香蕉久久夜色| 一进一出好大好爽视频| 久久久久久亚洲精品国产蜜桃av| 免费在线观看影片大全网站| a在线观看视频网站| 欧美乱码精品一区二区三区| 成年人免费黄色播放视频| 亚洲成国产人片在线观看| 欧美成人午夜精品| 露出奶头的视频| 欧洲精品卡2卡3卡4卡5卡区| 乱人伦中国视频| 久久亚洲真实| 无限看片的www在线观看| 亚洲熟妇中文字幕五十中出 | 国产高清videossex| 男人舔女人的私密视频| 亚洲精华国产精华精| 9191精品国产免费久久| 捣出白浆h1v1| 伦理电影免费视频| 啦啦啦免费观看视频1| 如日韩欧美国产精品一区二区三区| 国产精品久久久久成人av| 欧美黄色片欧美黄色片| 国产单亲对白刺激| 18禁观看日本| 成年人黄色毛片网站| 又黄又粗又硬又大视频| 亚洲熟妇熟女久久| 国产精品欧美亚洲77777| 91字幕亚洲| a级毛片在线看网站| 男男h啪啪无遮挡| 亚洲视频免费观看视频| 免费av中文字幕在线| 韩国精品一区二区三区| 国产精品偷伦视频观看了| 热99国产精品久久久久久7| 一二三四社区在线视频社区8| 少妇的丰满在线观看| 999精品在线视频| 亚洲欧美一区二区三区黑人| 真人做人爱边吃奶动态| 一本一本久久a久久精品综合妖精| 大码成人一级视频| 丝瓜视频免费看黄片| 欧美乱妇无乱码| 欧美乱色亚洲激情| videosex国产| 国产精品一区二区免费欧美| 欧美在线黄色| 在线av久久热| 韩国精品一区二区三区| tocl精华| 国产精品免费视频内射| а√天堂www在线а√下载 | 最近最新免费中文字幕在线| avwww免费| 精品国产一区二区久久| 少妇粗大呻吟视频| cao死你这个sao货| 久久久久国产精品人妻aⅴ院 | 久99久视频精品免费| 国产精品亚洲一级av第二区| 亚洲欧美一区二区三区黑人| 韩国av一区二区三区四区| 国产主播在线观看一区二区| 亚洲精品国产一区二区精华液| 自线自在国产av| 最近最新免费中文字幕在线| 国产高清国产精品国产三级| 午夜日韩欧美国产| 欧洲精品卡2卡3卡4卡5卡区| 欧美国产精品va在线观看不卡| 欧美黄色片欧美黄色片| 国产成人av教育| 韩国av一区二区三区四区| 搡老岳熟女国产| 最近最新免费中文字幕在线| 免费观看精品视频网站| 高清在线国产一区| 在线国产一区二区在线| 国产深夜福利视频在线观看| 丝袜美腿诱惑在线| 欧美国产精品一级二级三级| av不卡在线播放| 日韩成人在线观看一区二区三区| 黑丝袜美女国产一区| av网站在线播放免费| 少妇 在线观看| 在线天堂中文资源库| 久久久久视频综合| 国产精品免费视频内射| 中文字幕人妻丝袜一区二区| 国产极品粉嫩免费观看在线| 香蕉久久夜色| 美女午夜性视频免费| 高清毛片免费观看视频网站 | 亚洲国产欧美一区二区综合| 欧美大码av| 热re99久久精品国产66热6| 99re6热这里在线精品视频| 成熟少妇高潮喷水视频| 久久久久久久国产电影| 18禁美女被吸乳视频| 久久久久久久久免费视频了| 另类亚洲欧美激情| 国产精品久久电影中文字幕 | 欧美色视频一区免费| 一边摸一边做爽爽视频免费| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品一区二区www | 欧美不卡视频在线免费观看 | 精品一区二区三区四区五区乱码| 身体一侧抽搐| 久久精品亚洲av国产电影网| 热re99久久精品国产66热6| 午夜福利在线观看吧| 热99国产精品久久久久久7| 深夜精品福利| 精品熟女少妇八av免费久了| 国产成+人综合+亚洲专区| av在线播放免费不卡| 日本vs欧美在线观看视频| 丰满饥渴人妻一区二区三| 亚洲黑人精品在线| 精品久久久精品久久久| 亚洲性夜色夜夜综合| 日日爽夜夜爽网站| 国产一区二区三区综合在线观看| 国产成人免费无遮挡视频| 欧美老熟妇乱子伦牲交| www.精华液| 伊人久久大香线蕉亚洲五| 老鸭窝网址在线观看| 黑人欧美特级aaaaaa片| 国产蜜桃级精品一区二区三区| 久久草成人影院| 国产午夜精品久久久久久一区二区三区 | 噜噜噜噜噜久久久久久91| 国产色婷婷99| 亚洲精品色激情综合| 成人18禁在线播放| 久久99热这里只有精品18| e午夜精品久久久久久久| 国产精品综合久久久久久久免费| 精品一区二区三区视频在线 | 欧美丝袜亚洲另类 | 中文亚洲av片在线观看爽| 亚洲国产精品成人综合色| 精品久久久久久久人妻蜜臀av| 色吧在线观看| 亚洲在线观看片| 国产精品,欧美在线| 精品午夜福利视频在线观看一区| 在线视频色国产色| 亚洲第一欧美日韩一区二区三区| 国产一区二区在线观看日韩 | 最好的美女福利视频网| 国产欧美日韩一区二区精品| 精品一区二区三区人妻视频| 久9热在线精品视频| 国产一级毛片七仙女欲春2| 国产成人影院久久av| av福利片在线观看| 十八禁人妻一区二区| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 国产精品亚洲一级av第二区| 成年人黄色毛片网站| 女同久久另类99精品国产91| 十八禁人妻一区二区| 一区二区三区激情视频| 操出白浆在线播放| 亚洲精品亚洲一区二区| 热99re8久久精品国产| 12—13女人毛片做爰片一| 别揉我奶头~嗯~啊~动态视频| 亚洲成人久久性| 成人av在线播放网站| 日韩欧美 国产精品| 久久性视频一级片| 三级国产精品欧美在线观看| 国产真实乱freesex| 精品人妻偷拍中文字幕| 成人鲁丝片一二三区免费| 亚洲乱码一区二区免费版| 又爽又黄无遮挡网站| 亚洲自拍偷在线| 国产视频内射| 精品久久久久久,| 成人鲁丝片一二三区免费| 中文字幕人成人乱码亚洲影| 亚洲精品乱码久久久v下载方式 | 国产亚洲精品av在线| 欧美一区二区亚洲| 狠狠狠狠99中文字幕| 女人十人毛片免费观看3o分钟| 亚洲欧美精品综合久久99| 悠悠久久av| 免费看美女性在线毛片视频| 真实男女啪啪啪动态图| 日韩欧美三级三区| 国产探花在线观看一区二区| 欧美一区二区国产精品久久精品| 一级毛片高清免费大全| 精品久久久久久久毛片微露脸| 99久久99久久久精品蜜桃| 国产精品自产拍在线观看55亚洲| 欧美乱码精品一区二区三区| 观看美女的网站| 少妇高潮的动态图| 午夜福利视频1000在线观看| 日本免费一区二区三区高清不卡| x7x7x7水蜜桃| 日本与韩国留学比较| 久久精品国产99精品国产亚洲性色| 国产真实伦视频高清在线观看 | 国产一区二区激情短视频| 日本在线视频免费播放| 大型黄色视频在线免费观看| 天天躁日日操中文字幕| 757午夜福利合集在线观看| 91在线精品国自产拍蜜月 | 亚洲av免费高清在线观看| 色哟哟哟哟哟哟| 麻豆成人午夜福利视频| 此物有八面人人有两片| 色视频www国产| 亚洲精品一区av在线观看| 国产成人欧美在线观看| 精品国产三级普通话版| 成人国产综合亚洲| 一个人观看的视频www高清免费观看| 一级黄色大片毛片| 久久久国产成人免费| 亚洲精品亚洲一区二区| 婷婷丁香在线五月| 99精品欧美一区二区三区四区| 99精品在免费线老司机午夜| 国产精品精品国产色婷婷| 欧美日韩瑟瑟在线播放| 国产成人影院久久av| 一区二区三区高清视频在线| 精品无人区乱码1区二区| 丁香欧美五月| 中文字幕人妻熟人妻熟丝袜美 | 亚洲无线观看免费| 老鸭窝网址在线观看| av在线蜜桃| 三级毛片av免费| 午夜福利18| 亚洲精品乱码久久久v下载方式 | 熟妇人妻久久中文字幕3abv| 在线观看美女被高潮喷水网站 | 国产色爽女视频免费观看| 99久久精品热视频| 久久久久久久亚洲中文字幕 | 热99re8久久精品国产| 波多野结衣高清无吗| 少妇的丰满在线观看| 日韩 欧美 亚洲 中文字幕| 精品人妻一区二区三区麻豆 | 人妻丰满熟妇av一区二区三区| 国产亚洲欧美在线一区二区| 在线播放无遮挡| av在线天堂中文字幕| 网址你懂的国产日韩在线| 嫩草影院入口| 午夜福利高清视频| 亚洲精品456在线播放app | 神马国产精品三级电影在线观看| 亚洲成av人片免费观看| 亚洲 国产 在线| 欧美3d第一页| 成人三级黄色视频| 制服人妻中文乱码| 日本精品一区二区三区蜜桃| a级一级毛片免费在线观看| 美女高潮喷水抽搐中文字幕| 中文资源天堂在线| 亚洲av熟女| 在线观看免费视频日本深夜| 亚洲五月婷婷丁香| 亚洲精品一区av在线观看| 法律面前人人平等表现在哪些方面| 国产私拍福利视频在线观看| 啦啦啦韩国在线观看视频| 亚洲无线在线观看| 国产毛片a区久久久久| 99riav亚洲国产免费| 免费在线观看影片大全网站| 狂野欧美白嫩少妇大欣赏| 欧美av亚洲av综合av国产av| 国产v大片淫在线免费观看| 一a级毛片在线观看| 婷婷丁香在线五月| 日韩欧美三级三区| 午夜免费男女啪啪视频观看 | 午夜福利在线观看免费完整高清在 | 国产成人aa在线观看| 亚洲精品成人久久久久久| 国产aⅴ精品一区二区三区波| 岛国视频午夜一区免费看| 韩国av一区二区三区四区| 国模一区二区三区四区视频| 亚洲专区中文字幕在线| 国产精品永久免费网站| 丰满乱子伦码专区| 国产精品99久久99久久久不卡| 免费观看的影片在线观看| 在线观看av片永久免费下载| 中出人妻视频一区二区| 波多野结衣高清作品| 熟妇人妻久久中文字幕3abv| 国产亚洲av嫩草精品影院| 午夜影院日韩av| 亚洲第一电影网av| 琪琪午夜伦伦电影理论片6080| 欧美激情久久久久久爽电影| 1024手机看黄色片| 国产真人三级小视频在线观看| 欧美一级a爱片免费观看看| 日本在线视频免费播放| 丰满人妻一区二区三区视频av | 国产高清视频在线观看网站| 亚洲中文日韩欧美视频| 久久6这里有精品| 亚洲中文字幕一区二区三区有码在线看| 美女 人体艺术 gogo| 色在线成人网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一夜夜www| 九九在线视频观看精品| 午夜精品在线福利| 我的老师免费观看完整版| 日韩亚洲欧美综合| 国产一区二区三区视频了| a在线观看视频网站| 一级黄色大片毛片| 亚洲精华国产精华精| 国产精品99久久99久久久不卡| 俄罗斯特黄特色一大片| 亚洲国产日韩欧美精品在线观看 | 少妇的逼水好多| 国产一区二区激情短视频| 麻豆一二三区av精品| 九九热线精品视视频播放| 午夜免费男女啪啪视频观看 | 噜噜噜噜噜久久久久久91| 午夜免费男女啪啪视频观看 | 男人和女人高潮做爰伦理| 国产精品亚洲一级av第二区| 桃红色精品国产亚洲av| 精品一区二区三区视频在线观看免费| 俺也久久电影网| 精品不卡国产一区二区三区| 美女cb高潮喷水在线观看| 日本五十路高清| 国产高清videossex| 露出奶头的视频| 变态另类成人亚洲欧美熟女| 欧美色欧美亚洲另类二区| 高清日韩中文字幕在线| av在线天堂中文字幕| 在线观看免费视频日本深夜| aaaaa片日本免费| 99久久精品热视频| 韩国av一区二区三区四区| 国产亚洲av嫩草精品影院| 精品久久久久久,| 啦啦啦观看免费观看视频高清| 亚洲aⅴ乱码一区二区在线播放| 三级国产精品欧美在线观看| 91麻豆av在线| 变态另类丝袜制服| 欧美不卡视频在线免费观看| 成人av在线播放网站| 一区福利在线观看| 一级毛片高清免费大全| 国产成人影院久久av| 国产精品久久久久久亚洲av鲁大| 国内久久婷婷六月综合欲色啪| 亚洲国产欧美人成| 国产精品99久久99久久久不卡| 午夜激情福利司机影院| 国产一区二区在线观看日韩 | 在线观看av片永久免费下载|