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

    某高層建筑實測風(fēng)場和風(fēng)壓的相關(guān)性研究

    2015-06-22 14:07:56李正農(nóng)李紅益羅疊峰潘月月
    實驗流體力學(xué) 2015年4期
    關(guān)鍵詞:方根值時距風(fēng)場

    李正農(nóng), 李紅益, 羅疊峰, 潘月月

    (湖南大學(xué) 建筑安全與節(jié)能教育部重點實驗室, 長沙 410082)

    某高層建筑實測風(fēng)場和風(fēng)壓的相關(guān)性研究

    李正農(nóng), 李紅益*, 羅疊峰, 潘月月

    (湖南大學(xué) 建筑安全與節(jié)能教育部重點實驗室, 長沙 410082)

    通過對臺風(fēng)“天兔”登陸時廈門沿海某高層建筑的風(fēng)場及建筑物迎風(fēng)墻面風(fēng)壓的現(xiàn)場實測數(shù)據(jù)進行分析,研究了高層建筑風(fēng)場和其迎風(fēng)墻面風(fēng)壓的狀況及它們的相關(guān)性。對臺風(fēng)登陸前后的風(fēng)場和風(fēng)壓數(shù)據(jù)進行的分析表明:城市上空強風(fēng)“天兔”的湍流度不是很大,且隨著風(fēng)速的增加變化比較平穩(wěn);迎風(fēng)墻面各測點的風(fēng)壓系數(shù)呈現(xiàn)出中間大、兩頭小的特征,但差別不是很大。在不同平均風(fēng)時距和不同來流風(fēng)向角條件下,對迎風(fēng)墻面風(fēng)場和風(fēng)壓的相關(guān)性進行了分析。結(jié)果表明:當(dāng)基本時距取為3s時風(fēng)場同風(fēng)壓的相關(guān)性較高,當(dāng)來流風(fēng)向垂直于結(jié)構(gòu)迎風(fēng)墻面時,風(fēng)場同風(fēng)壓場的相關(guān)性比較大。

    高層建筑;現(xiàn)場實測;風(fēng)場;風(fēng)壓;相關(guān)性

    0 引 言

    由于現(xiàn)代高層建筑呈現(xiàn)出質(zhì)量輕、柔性大、阻尼低等特點,其對風(fēng)的敏感性也隨著這一特征的放大而增加,高層建筑(尤其是沿海多臺風(fēng)地區(qū)的高層建筑)的風(fēng)荷載研究變得越來越重要和必不可少。 現(xiàn)場實測作為結(jié)構(gòu)風(fēng)荷載研究最直接、最有效的方法之一,近年來在國內(nèi)風(fēng)工程領(lǐng)域逐步興起,并取得了一些實質(zhì)性的突破和發(fā)展。但由于受臺風(fēng)實測全尺度測量條件、露天環(huán)境惡劣、實測成本高等因素的制約,結(jié)構(gòu)抗風(fēng)的實測研究相對來說還比較少,特別是高層建筑結(jié)構(gòu)抗風(fēng)方面的實測數(shù)據(jù)更是缺乏,開展這方面的現(xiàn)場實測, 對高層建筑抗風(fēng)具有十分重要的意義[1-2]。

    研究分析各國的結(jié)構(gòu)設(shè)計荷載規(guī)范可知[3-4],在風(fēng)荷載計算時,首先需要確定的是基本風(fēng)速和基本風(fēng)壓,而這又與平均風(fēng)速的確定直接相關(guān)。平均風(fēng)速與所取某一中心時刻附近的時距有關(guān),隨著所取平均風(fēng)時距的縮短,對應(yīng)于這一時距的最大平均風(fēng)速將增大,因為在較小的時距內(nèi)能集中反映較大波峰的影響[5],而較小的波峰未能得以反應(yīng)。各國荷載規(guī)范對平均風(fēng)時距的規(guī)定不盡相同,國際上許多國家(包括我國) 將平均風(fēng)速計算時距值取為 10min,但也有國家取為1h(如加拿大等),甚至有的國家取為3~5s時距的瞬時風(fēng)速(如美國規(guī)范取為3s)。英國規(guī)范規(guī)定對所有圍護構(gòu)件、玻璃及屋面,都采用3s陣風(fēng)風(fēng)速。對于豎向和水平最大尺寸大于50m的房屋或結(jié)構(gòu)物,采用15s的平均風(fēng)速。顯然在風(fēng)荷載計算時,取不同的時距得到的結(jié)果是不同的。對于超高層建筑而言,究竟時距長度取何值時最為適宜,在國內(nèi)外尚沒有定論。本文通過對廈門觀音山商業(yè)運用中心11號樓在臺風(fēng)“天兔”登陸時的風(fēng)場及其迎風(fēng)墻面風(fēng)壓的實測數(shù)據(jù)進行了相關(guān)性分析, 研究了時距和來流風(fēng)向角對風(fēng)場和風(fēng)壓相關(guān)性大小的影響, 以供今后在高層建筑抗風(fēng)研究中,在基本時距選取時作一個參考。

    1 觀測位置和儀器

    為了獲取在臺風(fēng)“天兔”作用下超高層建筑的風(fēng)場特性及建筑迎風(fēng)墻面風(fēng)壓場特征,在位于廈門觀音山商業(yè)運營中心的11號樓頂部安裝了一個風(fēng)速風(fēng)向儀,風(fēng)速風(fēng)向儀安裝位置如圖1所示,其高出屋面柵欄6.4m,超出柵欄4.2m,比前方屋面高出約15m。并在該棟樓的第17層(避難層)幕墻外表面布置了18個風(fēng)壓測點。該樓距離海邊大概450m左右,地貌類別為A類。該建筑主體為36層,地上部分的高度為149.8m?,F(xiàn)場實測中采用RM Young 05103v螺旋槳式風(fēng)速風(fēng)向儀,在實測時,風(fēng)速儀零度風(fēng)向角對應(yīng)的地理方向為正北方向,角度按順時針方向遞增(圖2),此處定義結(jié)構(gòu)主軸X向為建筑物的長邊方向,結(jié)構(gòu)主軸Y向為建筑物的短邊方向;在該樓的第17層(避難層,離地面的高度約為70m)開展風(fēng)壓實測,風(fēng)壓傳感器采用CY2000FA1P型風(fēng)壓傳感器。風(fēng)壓傳感器粘貼在東北和東南的外墻面上,測點布置如圖2所示,共布置了18個測點(本文只采用了1~15測點的數(shù)據(jù))。當(dāng)臺風(fēng)“天兔”2013年9月20日16 時32分登陸福建沿海地區(qū)時,對該樓的樓頂風(fēng)場和避難層的風(fēng)壓情況進行了現(xiàn)場實測,風(fēng)速和風(fēng)壓的實測頻率分別為20Hz和25.6Hz,同步連續(xù)記錄了從2013年9月20日16時至9月23日9時長達(dá)65h的風(fēng)場及風(fēng)壓數(shù)據(jù)。

    圖1 實測現(xiàn)場信息圖Fig.1 The scene graph of measured information

    2 數(shù)據(jù)的選取

    論文在數(shù)據(jù)選取時,主要選取風(fēng)向角與迎風(fēng)墻面基本垂直范圍內(nèi)的風(fēng)場數(shù)據(jù)。同時,文獻[10]中指出風(fēng)壓傳感器所測到的風(fēng)壓理論上為垂直于墻面的風(fēng)壓,平行于墻面的部分可以忽略。因此,本文提取風(fēng)速垂直于墻面的分量進行分析。并且,在風(fēng)速的分解中發(fā)現(xiàn),所取風(fēng)向范圍內(nèi)的橫風(fēng)向風(fēng)速很小,基本可不做考慮。在后續(xù)的風(fēng)場參數(shù)的定義中均以垂直于墻面的風(fēng)速分量作為基本參數(shù)進行定義與計算分析,以保持風(fēng)場參數(shù)方向與所測風(fēng)壓一致。

    根據(jù)上文的分析,本文在進行數(shù)據(jù)選取時,選取了來流風(fēng)向角基本與建筑物迎風(fēng)墻面垂直,風(fēng)速變化比較平穩(wěn)的一段4h內(nèi)風(fēng)場時程作計算分析,風(fēng)壓取時間上同風(fēng)場同步的實測風(fēng)壓數(shù)據(jù)?;谝陨峡紤],本文對實測風(fēng)速進行如圖2所示的Y和X兩個方向進行實時分解,得到Y(jié)方向的風(fēng)速時程,在此基礎(chǔ)上再求Y方向的平均風(fēng)速和湍流度等風(fēng)場特性值。

    圖2 風(fēng)壓傳感器的平面布置圖Fig.2 Plan view of wind pressure sensor location

    圖3 實測流程線路圖Fig.3 Flowchart of measurement

    3 實測風(fēng)場特性

    3.1 平均風(fēng)速和風(fēng)向

    實測的風(fēng)速、風(fēng)向記錄分為2個時間序列,即水平風(fēng)速u(t)和風(fēng)向φ(t)。其隨時間變化的時程圖如圖5所示。風(fēng)速可根據(jù)以下公式分解為2個坐標(biāo)軸方向的分量ux(t)(建筑物長邊方向)和uy(t)(建筑物短邊方向),如圖4所示:

    (1)

    (2)

    計算分析取T為基本時距,Y方向的平均風(fēng)速U和水平風(fēng)向角可表示為:

    (3)

    (4)

    圖4 風(fēng)速、風(fēng)向示意圖Fig.4 Coordinate of wind speed and wind direction

    圖5 水平風(fēng)速風(fēng)向及Y方向風(fēng)速時程

    Fig.5 Time series of the horizontal wind speed and direction, wind speed inYdirection

    實測風(fēng)速儀零度風(fēng)向角對應(yīng)于地理上的正北方向,而上述Y方向與地理正北方向之間的夾角為39°(風(fēng)向角沿著順時針方向增大)。本文所有提及的風(fēng)場特征值(除風(fēng)向角外,其它的風(fēng)場特征值),如無特別說明,均指實際量在Y方向的分量。

    經(jīng)公式(2)實時分解后Y方向的風(fēng)速時程如圖5(c)所示,由于所取風(fēng)場數(shù)據(jù)段的風(fēng)向角基本與建筑迎風(fēng)墻面垂直(Y方向),因此經(jīng)公式(2)分解后Y方向的風(fēng)速時程同水平風(fēng)速時程相比,其值變化不大,變化趨勢基本一致,由于所取水平風(fēng)速可假定是平穩(wěn)隨機過程[6],因此Y方向的風(fēng)速時程也可以認(rèn)為是平穩(wěn)隨機過程。

    在本次臺風(fēng)實測的時程數(shù)據(jù)中選取4h的實測數(shù)據(jù)進行分析,以10min為基本時距的Y方向平均風(fēng)速和風(fēng)向圖如圖6所示,本次觀測的Y方向分量風(fēng)的總體平均風(fēng)速為11.08m/s,最大10min平均風(fēng)速為12.46m/s,順風(fēng)向的最大瞬時風(fēng)速為16.77m/s。實測的平均風(fēng)向角先變大,然后趨于平穩(wěn),平均風(fēng)向處于東北偏東方向。由圖及以上分析可知風(fēng)速和風(fēng)向變化較為平穩(wěn),表明臺風(fēng)登陸后,陸地上的風(fēng)速和風(fēng)向趨于穩(wěn)定。此外,同水平方向平均風(fēng)速相比,Y方向的平均風(fēng)速略小,但其變化趨勢基本一致。

    圖6 10min Y方向平均風(fēng)速和平均風(fēng)向時程,水平平均風(fēng)速時程

    Fig.6 Time series of 10min mean wind inYdirection and mean wind direction,the level of average wind speed

    3.2Y方向分量風(fēng)的脈動風(fēng)速均方根值和湍流度

    在結(jié)構(gòu)抗風(fēng)中,風(fēng)的湍流特征對于柔性結(jié)構(gòu)的共振和剛性結(jié)構(gòu)的振動有著重要的作用。湍流度是結(jié)構(gòu)抗風(fēng)設(shè)計中基本的湍流特征信息之一。Y方向風(fēng)的脈動風(fēng)速、湍流度、脈動風(fēng)速均方根值定義公式如下:

    (5)

    (6)

    (7)

    風(fēng)場中脈動風(fēng)速的變化過程可以看作是一個隨機過程[6-8],進而可利用隨機過程理論建立脈動風(fēng)速模型,在滿足工程計算精度要求的前提下,可以假定脈動風(fēng)速是零均值的平穩(wěn)隨機過程[9]。理論上,假定脈動風(fēng)速和脈動風(fēng)壓是平穩(wěn)隨機過程(對本文公式(5)等號兩邊求平均,可知等號右邊的值為0)。但是,通過現(xiàn)場實測得到的脈動風(fēng)速和脈動風(fēng)壓的均值不是0,而是接近于0的值,且隨著基本時距的增大而趨向于0。因此,本文采用脈動風(fēng)壓的均方根,并分析基本時距對兩者之間的相關(guān)性的影響規(guī)律。以10min作為基本時距,選擇垂直于建筑物迎風(fēng)墻面方向的脈動風(fēng)速均方根值做為脈動風(fēng)速特征值,其隨時間的變化歷程如圖7(a)所示。從圖中可以看出,在所取的4h內(nèi)的脈動風(fēng)速均方根值歷程中,在小范圍內(nèi)有較大波動,但波動范圍不大,總體變化比較平穩(wěn)。

    圖7(b)描述了Y方向?qū)崪y風(fēng)速的10min平均湍流度,湍流度的變化范圍為7.8%~12.3%,平均湍流度為9.3%,可以看出10min平均湍流度的變化平緩,波動不是很大。

    4 實測風(fēng)壓特性

    4.1 高空實測墻面的風(fēng)壓特征

    實測風(fēng)場是隨機的,因此風(fēng)速和風(fēng)向是隨機過程,導(dǎo)致建筑物表面所測得的風(fēng)壓顯然也是隨機過程[9]。由于來流與建筑物呈一交角,因此作用于建筑物表面的風(fēng)壓可分為沿法線方向和切線方向的2個分量,由于空氣的粘性很小,抗剪能力較差,因此起主要作用的是沿結(jié)構(gòu)法線方向的風(fēng)壓[10],即垂直于結(jié)構(gòu)表面的風(fēng)壓。同風(fēng)速一樣,風(fēng)壓也可以提取一系列代表風(fēng)壓的特征值,本文主要計算了建筑物表面實測的平均風(fēng)壓、脈動風(fēng)壓均方根值。平均風(fēng)壓的計算公式如下:

    圖7 脈動風(fēng)速均方根值和湍流度時程Fig.7 Variation of fluctuating wind speed and turbulence intensity

    (8)

    T為計算時所取基本時距,本節(jié)T的取值為10min,其它時距下也可以采用該式求解。

    由于風(fēng)場的隨機性,由前文可知,脈動風(fēng)速為一定時距上均值為零的平穩(wěn)隨機過程,因此脈動風(fēng)壓也可以認(rèn)為是一定時距上均值為零的平穩(wěn)隨機過程,因此,同3.2節(jié)中的脈動風(fēng)速一樣,將脈動風(fēng)壓均方根值作為風(fēng)壓脈動的特征值,用于后文的相關(guān)性分析,其定義如下:

    (9)

    T為計算所取的基本時距,同前面風(fēng)場特性的定義一致。 以10min為基本時距,計算了迎風(fēng)墻面測點7的平均風(fēng)壓和脈動風(fēng)壓均方根值與風(fēng)場同步4h內(nèi)時間變化歷程,如圖8所示。

    從圖8和圖6(b)、圖7(a)對比可以看出,風(fēng)壓的變化趨勢同風(fēng)場特性的變化趨勢相一致,說明實測風(fēng)場數(shù)據(jù)與實測的風(fēng)壓數(shù)據(jù)呈現(xiàn)某種程度的相關(guān)性,對于迎風(fēng)墻面的其它測點采用同樣的方法分析,得到類似的結(jié)論。對于其它時距條件,也可以采用同樣的方法進行分析與計算。

    圖8 10min平均風(fēng)壓和脈動風(fēng)壓均方根值時程

    Fig.8 Time series of 10min mean wind pressure and the RMS value of fluctuating wind pressure

    4.2 風(fēng)壓系數(shù)

    根據(jù)文獻[11],在風(fēng)壓的全尺度實測中,平均風(fēng)壓系數(shù)定義如下[11]:

    (10)

    各測點的風(fēng)壓系數(shù)值如表1所示。平均風(fēng)壓系數(shù)呈現(xiàn)中間大、兩頭小的規(guī)律,且均為小于1的正值,但靠近東邊的測點與靠近西邊的對稱測點相比,平均風(fēng)壓系數(shù)要大,這可能與風(fēng)是從東北偏東方向吹來有關(guān)。各測點均方根風(fēng)壓系數(shù)有類似規(guī)律。

    表1 各測點風(fēng)壓參數(shù)值Table 1 Pressure coefficients of each test point

    圖9為各測點在10min時距下的平均風(fēng)壓系數(shù)隨風(fēng)向角變化趨勢。由圖可知,迎風(fēng)墻面各測點隨著風(fēng)向角的變化趨勢基本趨于一致,從中間測點往兩邊測點呈現(xiàn)遞減的規(guī)律,這可能與轉(zhuǎn)角處的脈動較大有關(guān)。在不同風(fēng)向角下,靠近北面的各測點平均風(fēng)壓系數(shù)隨風(fēng)向角的增大而減小,而靠近南面的各測點的平均風(fēng)壓系數(shù)隨風(fēng)向角的變化趨勢剛好相反,中間測點隨著風(fēng)向角呈現(xiàn)先增大后減小的趨勢。這說明同來流風(fēng)向靠近的測點風(fēng)壓更大,而遠(yuǎn)端的測點相對要小。中間測點始終是最大的,這可能與靠邊測點的脈動較大有關(guān)。

    圖9 不同風(fēng)向角下各測點的平均風(fēng)壓系數(shù)

    Fig.9 Mean pressure coefficient distributions of every test point in different wind directions

    5 相關(guān)性分析

    相關(guān)性分析是指對兩個或多個具備相關(guān)性的變量元素進行分析,從而衡量兩個變量因素相關(guān)程度的一種分析手段。參與分析的各元素之間需要存在一定的聯(lián)系或者概率才可以進行相關(guān)性分析[12]。由前文可知,風(fēng)場和風(fēng)壓顯然滿足相關(guān)性分析的條件。

    為了探討基本時距和風(fēng)向角對風(fēng)場和風(fēng)壓相關(guān)性的影響, 本文將根據(jù)前述實測所獲得的風(fēng)場數(shù)據(jù)和風(fēng)壓場數(shù)據(jù),通過相關(guān)性分析的方法研究在不同的時距條件下兩者之間的相關(guān)性, 進而通過比較在不同時距條件下的相關(guān)性,確定相關(guān)性較大時所對應(yīng)的時距。另外對比了10min基本時距時,不同來流風(fēng)向角下二者的相關(guān)性,以期得到相關(guān)性最大時來流的角度。具體的計算方法和公式介紹如下。

    5.1 相關(guān)性及其計算[13]

    對于一組復(fù)雜數(shù)據(jù)可以使用一定的數(shù)學(xué)手段產(chǎn)生一個相對簡單的類結(jié)構(gòu),進行“相關(guān)性”或“相似性”分析,最常用的是使用相似系數(shù)Cxy來表示指標(biāo)x和指標(biāo)y之間的相似關(guān)系。Cxy的取值在-1~1之間,其絕對值越接近1,表示指標(biāo)x和指標(biāo)y之間的關(guān)系越密切,其絕對值越接近0,表示指標(biāo)x和指標(biāo)y之間的關(guān)系越疏遠(yuǎn)。對于兩條曲線形狀變化規(guī)律之間的間隔尺度,常用的相似系數(shù)有夾角余弦和相關(guān)系數(shù)。

    (1)夾角余弦;這是受相似形的啟發(fā)而來,對于兩條曲線,如果長度不一,但是形狀相似,當(dāng)長度不是主要矛盾時,可以定義一種相似系數(shù)來表示兩條曲線之間比較密切的相似關(guān)系從而加以描述。通常用夾角余弦來描述這種關(guān)系,其定義為

    (11)

    式中,vkx,vky表示所要研究的兩個指標(biāo)向量,在本文研究分析中,vkx分別對應(yīng)于平均風(fēng)速、脈動風(fēng)速均方根值和湍流度;vky分別對應(yīng)于平均風(fēng)壓、脈動風(fēng)壓均方根值,k表示指標(biāo)向量的第k個元素,n表示指標(biāo)向量中元素的個數(shù)。

    (2)相關(guān)系數(shù);相關(guān)系數(shù)Cxy就是求將數(shù)據(jù)標(biāo)準(zhǔn)化 (將指標(biāo)向量的每個元素減去指標(biāo)向量元素平均值) 后的夾角余弦。為了簡便計,在計算分析時,相關(guān)性系數(shù)取計算各指標(biāo)向量之間的夾角余弦。

    5.2 相關(guān)性分析結(jié)果

    在研究不同時距和不同來流風(fēng)向角下風(fēng)場與建筑物迎風(fēng)墻面風(fēng)壓場的相關(guān)性時,集中計算了在不同平均風(fēng)時距和不同風(fēng)向角下平均風(fēng)速同平均風(fēng)壓,脈動風(fēng)速均方根值、湍流度分別同脈動風(fēng)壓均方根值的相關(guān)性。從數(shù)理統(tǒng)計的角度而言,參與計算的數(shù)據(jù)個數(shù)的多少對于相關(guān)性系數(shù)的值有一定的影響。

    5.2.1 不同的基本時距條件下風(fēng)場和風(fēng)壓的相關(guān)性分析結(jié)果

    在總的采樣時間長度相同,不同基本時距條件下劃分的數(shù)據(jù)段個數(shù)顯然是不同的,平均風(fēng)時距越大劃分的數(shù)據(jù)段個數(shù)越少,而采用上述相關(guān)性分析方法,數(shù)據(jù)段個數(shù)的多少會直接影響到計算得到的相關(guān)系數(shù)值。因此,為了體現(xiàn)最后結(jié)果的一般性,本文分別取等時間長度(取計算總時長為1h)和等數(shù)據(jù)段個數(shù)(在不同時距條件取相同的數(shù)據(jù)段個數(shù)(60個))兩種情況進行分析。本文選取了1s,2s,3s,10s,0.5min,1min,3min,10min共8個時距下的相關(guān)性分析。

    (1) 等時間長度時的相關(guān)度

    取計算的時間長度為1h,在此條件下,計算得到了在不同時距條件下,平均風(fēng)速同平均風(fēng)壓以及脈動風(fēng)速均方根值和湍流度分別同脈動分壓均方根值的的相關(guān)性。具體計算結(jié)果如圖10所示。

    由圖10可知,對于脈動風(fēng)速均方根值同脈動風(fēng)壓均方根值以及湍流度同脈動風(fēng)壓均方根值的相關(guān)性變化趨勢基本一致,由脈動風(fēng)速均方根值同湍流度的關(guān)系可知,此時平均風(fēng)速對其相關(guān)性的影響較小。具體來說,當(dāng)時距小于3s時,相關(guān)性呈現(xiàn)陡降的趨勢,當(dāng)時距取為1s時,基本上沒有相關(guān)性。這可能與脈動風(fēng)速的周期正好處于這一時距段,因而得到的數(shù)據(jù)隨機性大有關(guān)。當(dāng)時距處于3s~0.5min的區(qū)間時,相關(guān)性呈現(xiàn)穩(wěn)步遞增的趨勢,變化相對比較平穩(wěn)。當(dāng)時距大于0.5min時,相關(guān)性系數(shù)基本呈一常數(shù)。分析原因可知,隨著時距的增大,脈動的瞬時性減弱,而呈現(xiàn)出長時距上的整體平均,偏離脈動風(fēng)自身的周期越大,這種現(xiàn)象越為明顯。平均風(fēng)速同平均風(fēng)壓的相關(guān)性在時距大于等于3s這一時距區(qū)間內(nèi)基本呈現(xiàn)和上述脈動風(fēng)速均方根值和脈動風(fēng)壓均方根值之間相同的變化趨勢。當(dāng)時距小于3s時,雖然相關(guān)性系數(shù)有所下降,但下降得不多,且在這一低時距區(qū)間內(nèi),還呈現(xiàn)出一定的相關(guān)性。這可能與平均風(fēng)的計算分析與脈動風(fēng)的周期計算相關(guān)不大有關(guān)系,再有,當(dāng)時距小時參與相關(guān)性計算的數(shù)據(jù)個數(shù)比時距大時要多,更能詳盡描述兩條曲線的變化趨勢,因此,在較小的時距段均呈現(xiàn)較低或幾乎沒有相關(guān)性。此外,為了分析迎風(fēng)墻面各個測點的相關(guān)性,分別對迎風(fēng)墻面中間測點7、邊測點1和邊測點15進行了分析,由圖可知3組數(shù)據(jù)下的相關(guān)性曲線走勢基本相同,對其它的測點進行同樣的分析也有類似的結(jié)果。這說明,迎風(fēng)墻面各個測點的風(fēng)壓同風(fēng)場有一致的相關(guān)性。

    圖10 風(fēng)場和風(fēng)壓的相關(guān)性

    Fig.10 The correlation between the measured wind field and wind pressure

    (2) 等數(shù)據(jù)個數(shù)時相關(guān)性結(jié)果

    在相關(guān)性系數(shù)計算的數(shù)學(xué)模型中,參與計算的數(shù)據(jù)個數(shù)對相關(guān)性系數(shù)值有影響。當(dāng)參與分析的數(shù)據(jù)個數(shù)比較少時,不能詳盡描述2個隨機變量之間或2條曲線變化趨勢的相似程度,這往往會導(dǎo)致計算出來的相關(guān)性系數(shù)比較大。為了消除這種影響,對不同時距條件下的數(shù)據(jù)均取60個進行分析。

    由圖11可知,取相同數(shù)據(jù)個數(shù)分析時,在時距大的區(qū)間內(nèi)相關(guān)性有所降低,而時距小的區(qū)間內(nèi)相關(guān)性有所增加。具體來說,對平均風(fēng)速和平均風(fēng)壓之間的相關(guān)性,當(dāng)時距小于等于2s時,基本為一常數(shù),同圖10比較可知,在這一時距區(qū)間內(nèi)相關(guān)性增大,說明參與分析的數(shù)據(jù)個數(shù)越少,相關(guān)性會越大。3s~2s的區(qū)間存在陡降趨勢,當(dāng)大于3s時,呈現(xiàn)遞減趨勢,但是變化不大。在這一時距區(qū)間內(nèi)有所減小,說明時距的影響對平均風(fēng)成分的影響是有限的。對脈動風(fēng)速均方根值和湍流度同脈動風(fēng)壓均方根值的相關(guān)性,當(dāng)時距小于等于2s時,相關(guān)性基本呈一常數(shù),3s~2s的區(qū)間存在陡降趨勢。說明所取時距小于脈動風(fēng)速的周期時,由于脈動風(fēng)場的瞬時性,在此時距區(qū)間內(nèi)的相關(guān)性較低,由于參與分析的數(shù)據(jù)個數(shù)相對圖10的情況來說很少,因此分析時還是有很小的相關(guān)性。3s~0.5min的區(qū)間內(nèi),基本呈一常數(shù),當(dāng)時距大于0.5min時,相關(guān)性呈下降趨勢,且下降比較多。這可能與所取時距過長,脈動成分表現(xiàn)出被中和抵消的情況有關(guān),但同圖10的情況相比,相關(guān)性下降較多,這與在大的時距區(qū)間參與分析數(shù)據(jù)的個數(shù)比圖10的情況要多有關(guān)。綜合以上分析可知,基本時距取3s時,風(fēng)場和風(fēng)壓之間各參數(shù)的相關(guān)性較大。

    5.2.2 不同來流風(fēng)向角下相關(guān)性結(jié)果

    取來流風(fēng)向角穩(wěn)定在42°~65°時,1h的風(fēng)場和風(fēng)壓數(shù)據(jù),集中研究了當(dāng)基本時距取10min時,在不同風(fēng)向角下平均風(fēng)速同平均風(fēng)壓,脈動風(fēng)速均方根值和湍流度分別同脈動風(fēng)壓均方根值的相關(guān)性。當(dāng)來流風(fēng)向同建筑物表面呈一夾角時,隨著相交角度的變化,風(fēng)壓隨著風(fēng)速變化的影響程度顯然是不一樣的。在不同風(fēng)向角下,風(fēng)場同風(fēng)壓的相關(guān)性如圖12所示。

    圖11 風(fēng)場和風(fēng)壓的相關(guān)性

    Fig.11 The correlation between the measured wind field and wind pressure

    由圖12可知,隨著來流角度的增大,相關(guān)性系數(shù)呈現(xiàn)先增大后減小的規(guī)律,表明來流風(fēng)向角對相關(guān)性影響較大,由于規(guī)定正北方向為來流風(fēng)向的零度角,當(dāng)來流風(fēng)向為39°時,來流風(fēng)速恰好與建筑物迎風(fēng)墻面垂直,由圖可知當(dāng)風(fēng)向角為39.4°,各測點風(fēng)場同風(fēng)壓之間各個對應(yīng)參數(shù)之間的相關(guān)性最大。說明來流風(fēng)向與迎風(fēng)墻面垂直時,相關(guān)性最大。對于其它風(fēng)向角下的相關(guān)性,當(dāng)風(fēng)向角從39.4°往兩邊變化基本呈對稱遞減的特征。因為建筑物表面粘貼的風(fēng)壓傳感器所能感應(yīng)到的風(fēng)是垂直于建筑物表面方向的風(fēng)速分量,因此由分析可知,當(dāng)風(fēng)速垂直吹向風(fēng)壓傳感器時,測到的風(fēng)壓更能準(zhǔn)確地反映風(fēng)速大小的變化,兩者的變化趨勢自然會呈現(xiàn)高度的相似,因此實測相關(guān)性分析得到的結(jié)果同我們理論分析的結(jié)果相同。

    圖12 風(fēng)場和風(fēng)壓的相關(guān)性

    Fig.12 The correlation between the measured wind field and wind pressure

    6 結(jié) 論

    通過上述對風(fēng)場和建筑物迎風(fēng)墻面實測風(fēng)壓及風(fēng)場參數(shù)和風(fēng)壓場參數(shù)的相似性分析得到的結(jié)果,得出以下結(jié)論:

    (1) 通過對風(fēng)場特性及風(fēng)壓特性實測數(shù)據(jù)的分析,強風(fēng)“天兔”的風(fēng)速和風(fēng)向變化比較平穩(wěn),湍流度不是很大,脈動風(fēng)速均方根值比較小,且隨著風(fēng)速的增加變化比較平穩(wěn);

    (2) 實測風(fēng)壓場與Y方向風(fēng)速的變化趨勢類似,對于迎風(fēng)墻面各測點風(fēng)壓系數(shù)呈現(xiàn)出中間大、兩頭小的特征,但差別不是很大;

    (3) 通過對風(fēng)場參數(shù)與風(fēng)壓參數(shù)的相關(guān)性分析,當(dāng)取等時長1h進行分析時,由于各個時距下,參與分析的數(shù)據(jù)個數(shù)不同,在時距大于等于0.5min時,相關(guān)性最高且基本呈一常數(shù),當(dāng)時距小于3s時,相關(guān)性較低,尤其是脈動風(fēng)壓均方根值同湍流度及脈動風(fēng)速均方根值的相關(guān)度非常小,下降比較快。當(dāng)基本時距為3s時,呈現(xiàn)比較高的相關(guān)性,3s~0.5min的區(qū)間內(nèi)穩(wěn)步上升。

    (4) 當(dāng)取相同的數(shù)據(jù)個數(shù)進行分析時,時距為3s時的相關(guān)性最高,3s~2s的區(qū)間內(nèi)呈現(xiàn)陡降趨勢,小于等于2s時相關(guān)性基本呈一常數(shù),且數(shù)值較低。當(dāng)時距大于等于3s時,隨著時距的增加,平均風(fēng)速同平均風(fēng)壓之間的相關(guān)性越來越小,但減小比較緩慢。脈動風(fēng)壓均方根值同脈動風(fēng)速均方根值以及湍流度的相關(guān)度,在3s~0.5min的區(qū)間內(nèi)基本不變,當(dāng)時距大于0.5min時,相關(guān)性值越來越小,且下降比較多。

    (5) 當(dāng)基本時距取3s時,風(fēng)場同風(fēng)壓場的相關(guān)性較高;當(dāng)來流風(fēng)向垂直于結(jié)構(gòu)迎風(fēng)墻面時,風(fēng)場同風(fēng)壓場的相關(guān)系數(shù)值較高,以上結(jié)果是在特定臺風(fēng)下的研究結(jié)果,結(jié)論是否具有普遍性,還有待進一步探索。

    [1] Li Q S, Fang J Q, Jeary A P, et al. Full-scale measurement of wind effects on tall buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1998, 74-76: 741-750.

    [2] 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 rein-forced concrete super-tall building[J]. Engineering Structures, 2004, 26: 1779-1794.

    [3] 張相庭. 國內(nèi)外風(fēng)載規(guī)范的評估和展望[J]. 同濟大學(xué)學(xué)報: 自然科學(xué)版, 2002, 30(5): 539-543.

    Zhang Xiangting. Evaluation and prosper for wind loading codes at home and abroad[J]. Journal Tongji University: Natural Science Edition, 2002, 30(5): 539-543.

    [4] GB50009-2001. 建筑結(jié)構(gòu)荷載規(guī)范[S]. 北京: 中國建筑工業(yè)出版社, 2006.

    GB50009-2001. Load code for the design building structure[S]. Beijing: China Architecture & Building Press, 2006.

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

    Huang Bencai, Wang Congjun. Wind resistant analysis of the principle and application[M]. 2nd ed. Shanghai: Tongji University, 2008: 28.

    [6] 陳嚴(yán), 張錦源, 王楠, 等. 風(fēng)力機風(fēng)場模型的研究及紊流風(fēng)場的Matlab數(shù)值模擬[J]. 太陽能學(xué)報, 2006, 27(9): 955-960.

    Chen Yan, Zhang Jinyuan, Wang Nan, et al. Study on turbine wind-field model and numerical simulation in turbulent flow[J]. Acta Energiae Solaris Sinica, 2006, 27(9): 955-960.

    [7] 谷興凱, 范高鋒, 王曉蓉, 等. 風(fēng)電功率預(yù)測技術(shù)綜述[J]. 電網(wǎng)技術(shù), 2007, 31(S2): 335-338.

    Gu Xingkai, Fan Gaofeng, Wang Xiaorong, et al. Summarization of wind power technology[J]. Power System Technology, 2007, 31(S2): 335-338.

    [8] 婁素華, 李志恒, 高蘇杰, 等. 風(fēng)電場模型及其對電力系統(tǒng)的影響[J]. 電網(wǎng)技術(shù), 2007, 31(S2): 330-334.

    Lou Suhua, Li Zhiheng, Gao Sujie, et al. Wind farms models and its impacts on wind farms integration into power system[J]. Power System Technology, 2007, 31(S2): 330-334.

    [9] 舒新玲, 周岱. 風(fēng)速時程AR模型及其快速實現(xiàn)[J]. 空間結(jié)構(gòu), 2004, 9(4): 27-32.

    Shu Xinling, Zhou Dai. AR model of wind speed time series and its rapid implementation[J]. Spatial Structures, 2004, 9(4): 27-32.

    [10] 楊慶山, 黃韜穎. 中美澳三國風(fēng)荷載規(guī)范簡介及其關(guān)于設(shè)計風(fēng)速規(guī)定之比較[C]//第十三屆全國結(jié)構(gòu)風(fēng)工程學(xué)術(shù)會議論文集 (下冊), 2007.

    Yang Qingshan, Huang Taoying. Brief introduction and comparative study on the design wind speed regulation of wind loads among the codes of China, America and Australia[C]// Proceedings of the Thirteenth National Academic Conference on Wind Engineering Structure, 2007.

    [11] Simiu E, Scanlan R H. Wind effects on structures[M]. Wiley, 1996.

    [12] 白麗, 張向陽. 廣州大學(xué)生網(wǎng)上創(chuàng)業(yè)與高校創(chuàng)業(yè)教育相關(guān)性探討[J]. 繼續(xù)教育研究, 2012, 2012(6): 123-125.

    Bai Li, Zhang Xiangyang. Correlation discussion between Guangzhou college students′ internet entrepreneurship and entrepreneurship education in colleges and universities[J]. Continue Education Research, 2012, 2012(6): 123-125.

    [13] 李正農(nóng), 宋克, 李秋勝, 等. 廣州中信廣場臺風(fēng)特性與結(jié)構(gòu)響應(yīng)的相關(guān)性分析[J]. 實驗流體力學(xué), 2009, 23(04): 21-27.

    Li Zhengnong, Song Ke, Li Qiusheng, et al. Correlation analysis of wind characteristic and wind-induced response of CITIC Plaza[J]. Journal of Experiments in Fluid Mechanics, 2009, 23(4): 21-27.

    (編輯:張巧蕓)

    Study of the correlation between the measured wind field and wind pressure on a high-rise building

    Li Zhengnong, Li Hongyi*, Luo Diefeng, Pan Yueyue

    (Key Laboratory of Building Safety and Efficiency of the Ministry of Education, Hunan University, Changsha 410082, China)

    Based on the field measurement and wind pressure measurement on the surface of a super-tall building on the coast of the Xiamen city during the passage of typhoon Usagi in September 2013, the wind field was measured by a propeller anemometer placed at the roof of the test building, and the wind pressure was measured by 18 pressure sensors stuck on the glass wall on the refuge floor at the height of about 70m. This paper presents selected measurement data (such as the mean wind direction, mean wind speed, turbulence intensity, mean pressure, and soon). Through in-depth analysis it is shown that the 10min wind direction and wind speed are relatively stable,the filed turbulence intensity of typhoon Usagi is normal, and with the increase of wind velocity the turbulence intersity varies steadily. The wind pressure coefficient shows the following characteristics: the wind pressure coefficient of each measuring point on the windward side shows a relatively large value in the middlepoint and smaller values from the middle to each end point, but the gap is not great. Meanwhile, the correlation between the measured wind field(the mean and RMS wind velocity,the turbulence intensity) and wind pressure (mean and RMS wind pressure)on the windward side is analyzed from two aspects. One is to investigate the correlation in 8 selected different time intervals: 1s, 2s, 3s, 10s, 30s, 1min, 3min, 10min, and the results of the research clearly reveal that: when the average time interval is taken as 3s, there is a relatively high correlation between the measured wind field and wind pressure. The other one is to evaluate the correlation under different wind directions: 30.8°, 33.2°, 36.7°, 39.4°, 42.1°, 45.3°, 47.6°, 49.5°, and as the wind direction is roughly perpendicular to the windward wall (when the wind direction is 39.4°), the correlation between the wind field and the wind pressure is relatively large.

    high-rise building;field measurement;wind field;wind pressure;correlation

    1672-9897(2015)04-0032-09

    10.11729/syltlx20140125

    2014-11-03;

    2015-04-02

    國家自然科學(xué)基金資助(91215302,51178180,51278190)

    LiZN,LiHY,LuoDF,etal.Studyofthecorrelationbetweenthemeasuredwindfieldandwindpressureonahigh-risebuilding.JournalofExperimentsinFluidMechanics, 2015, 29(4): 32-40. 李正農(nóng), 李紅益, 羅疊峰, 等. 某高層建筑實測風(fēng)場和風(fēng)壓的相關(guān)性研究. 實驗流體力學(xué), 2015, 29(4): 32-40.

    P425.6+5

    A

    李正農(nóng)(1962-),男,湖北武漢人,教授,博士生導(dǎo)師。研究方向:結(jié)構(gòu)抗震抗風(fēng)。通信地址:湖南省長沙市岳麓區(qū)湖南大學(xué)土木工程學(xué)院(410082)。E-mail:zhn88@263.net

    *通信作者 E-mail: 799425261@qq.com

    猜你喜歡
    方根值時距風(fēng)場
    磁流變彈性體減振單元動力學(xué)分析
    趣味性進度提示設(shè)計對時距知覺的影響
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    客車平順性仿真及優(yōu)化
    學(xué)習(xí)障礙學(xué)生與普通學(xué)生短時距估計的比較
    回溯式時距估計的計時機制*
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場中無人機的飛行研究
    標(biāo)量計時模型的影響因素及發(fā)展*
    Studyof velocityfluctuations in the plenum of a 3/4open jet automotive wind tunnel
    99久久精品国产国产毛片| 亚洲av在线观看美女高潮| 精品亚洲成国产av| 亚洲精品一二三| 免费看不卡的av| 久久久久久久久久成人| 一个人免费看片子| av有码第一页| 国产黄频视频在线观看| 久久人人爽人人爽人人片va| 啦啦啦视频在线资源免费观看| 色婷婷久久久亚洲欧美| 国产一级毛片在线| 亚洲精品国产av蜜桃| 国产成人一区二区在线| 99视频精品全部免费 在线| 蜜臀久久99精品久久宅男| 日韩制服骚丝袜av| 看免费av毛片| 51国产日韩欧美| 亚洲精品成人av观看孕妇| 久久99热6这里只有精品| av线在线观看网站| 国产成人午夜福利电影在线观看| 国产有黄有色有爽视频| 中文字幕人妻熟女乱码| 国产高清国产精品国产三级| 99久国产av精品国产电影| 超碰97精品在线观看| 老女人水多毛片| 一级爰片在线观看| 国产毛片在线视频| 免费人妻精品一区二区三区视频| 少妇的逼好多水| 黄色毛片三级朝国网站| 熟女人妻精品中文字幕| 丝袜喷水一区| 国产成人精品在线电影| 国产成人精品福利久久| 母亲3免费完整高清在线观看 | 大香蕉久久网| 中文字幕另类日韩欧美亚洲嫩草| 日产精品乱码卡一卡2卡三| 人人妻人人澡人人看| 一级a做视频免费观看| 国产1区2区3区精品| 国产精品久久久久成人av| 国产精品嫩草影院av在线观看| 免费看av在线观看网站| 日韩精品免费视频一区二区三区 | 免费人妻精品一区二区三区视频| 一级片'在线观看视频| 2021少妇久久久久久久久久久| 国产精品国产三级国产专区5o| 国产色爽女视频免费观看| 久久久久人妻精品一区果冻| 亚洲精品国产av成人精品| 国产乱来视频区| 午夜福利视频在线观看免费| 九草在线视频观看| 国产福利在线免费观看视频| 午夜激情久久久久久久| 日韩免费高清中文字幕av| 国产69精品久久久久777片| 极品少妇高潮喷水抽搐| 国产白丝娇喘喷水9色精品| 91午夜精品亚洲一区二区三区| 国产一区二区三区av在线| 国产综合精华液| 婷婷色综合www| 久久久久久人人人人人| 久热这里只有精品99| 免费高清在线观看日韩| 国产一区二区三区av在线| 多毛熟女@视频| 精品少妇内射三级| 欧美 日韩 精品 国产| 日韩视频在线欧美| 亚洲激情五月婷婷啪啪| 亚洲精品色激情综合| 亚洲色图综合在线观看| 欧美 日韩 精品 国产| 免费看不卡的av| 亚洲欧美精品自产自拍| 中文天堂在线官网| 成人黄色视频免费在线看| 欧美精品国产亚洲| 午夜福利视频精品| 亚洲国产成人一精品久久久| 80岁老熟妇乱子伦牲交| 视频在线观看一区二区三区| 多毛熟女@视频| 婷婷色综合www| 国产69精品久久久久777片| 亚洲国产看品久久| 国产在线免费精品| 国产 一区精品| 国国产精品蜜臀av免费| 午夜福利影视在线免费观看| 国产乱人偷精品视频| 久久鲁丝午夜福利片| 22中文网久久字幕| 视频中文字幕在线观看| 婷婷色av中文字幕| 18禁国产床啪视频网站| 亚洲精品自拍成人| 欧美少妇被猛烈插入视频| 中国国产av一级| 五月玫瑰六月丁香| 欧美人与性动交α欧美软件 | 午夜日本视频在线| 久久久久视频综合| 亚洲精品成人av观看孕妇| 日韩电影二区| 纵有疾风起免费观看全集完整版| 国产激情久久老熟女| 成人无遮挡网站| 久久久国产欧美日韩av| 国产熟女欧美一区二区| 国产欧美另类精品又又久久亚洲欧美| 在现免费观看毛片| 国产高清三级在线| 国产av精品麻豆| 亚洲精品久久久久久婷婷小说| 啦啦啦在线观看免费高清www| 日本爱情动作片www.在线观看| 韩国av在线不卡| 日韩精品免费视频一区二区三区 | 97在线视频观看| 午夜福利影视在线免费观看| 亚洲美女搞黄在线观看| 亚洲av电影在线进入| 亚洲欧洲国产日韩| 最后的刺客免费高清国语| 国产精品一二三区在线看| 人妻一区二区av| 亚洲精品自拍成人| 亚洲国产精品一区三区| 国产一级毛片在线| 国产白丝娇喘喷水9色精品| av福利片在线| 国产精品久久久久久av不卡| 三上悠亚av全集在线观看| 精品少妇黑人巨大在线播放| 哪个播放器可以免费观看大片| 久久热在线av| a级毛色黄片| 成年人免费黄色播放视频| 香蕉国产在线看| 在线观看国产h片| 国产 精品1| 国产又爽黄色视频| 母亲3免费完整高清在线观看 | av一本久久久久| 一二三四在线观看免费中文在 | 丝瓜视频免费看黄片| 午夜免费观看性视频| 国产熟女欧美一区二区| 日日撸夜夜添| 黄色视频在线播放观看不卡| 久久久久久久亚洲中文字幕| 亚洲第一区二区三区不卡| 日本wwww免费看| 中文字幕精品免费在线观看视频 | 街头女战士在线观看网站| 99热网站在线观看| 精品久久久精品久久久| 97在线人人人人妻| 99国产综合亚洲精品| 国产男女超爽视频在线观看| 亚洲av男天堂| 狂野欧美激情性bbbbbb| 99re6热这里在线精品视频| 国产精品无大码| 制服人妻中文乱码| 极品少妇高潮喷水抽搐| 你懂的网址亚洲精品在线观看| 成人国产麻豆网| 高清不卡的av网站| 日本91视频免费播放| 亚洲成人av在线免费| 精品一区二区三区四区五区乱码 | 国产精品免费大片| 一边摸一边做爽爽视频免费| 午夜免费观看性视频| 国产一区二区三区av在线| 久久久久久久精品精品| 涩涩av久久男人的天堂| 一级片免费观看大全| 日本av手机在线免费观看| 香蕉精品网在线| 国产在线视频一区二区| 久久精品久久久久久噜噜老黄| 一区二区三区四区激情视频| 一个人免费看片子| 满18在线观看网站| 日本av手机在线免费观看| 在线亚洲精品国产二区图片欧美| 中文字幕另类日韩欧美亚洲嫩草| 97在线视频观看| 18禁观看日本| 国产av码专区亚洲av| 中国三级夫妇交换| 9热在线视频观看99| 久久婷婷青草| 成人手机av| 国产精品一区二区在线不卡| 丝袜喷水一区| 美女大奶头黄色视频| 激情五月婷婷亚洲| 国产精品蜜桃在线观看| 永久免费av网站大全| 夫妻午夜视频| 飞空精品影院首页| 免费av不卡在线播放| 色婷婷av一区二区三区视频| 一本久久精品| 五月开心婷婷网| 亚洲激情五月婷婷啪啪| 99久久中文字幕三级久久日本| 一本久久精品| 欧美97在线视频| 亚洲国产精品一区二区三区在线| 欧美精品一区二区免费开放| 久久青草综合色| 国产无遮挡羞羞视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲精品久久久久久婷婷小说| 成人18禁高潮啪啪吃奶动态图| 免费看av在线观看网站| 99精国产麻豆久久婷婷| 国产男人的电影天堂91| 国产在视频线精品| 亚洲中文av在线| av片东京热男人的天堂| 曰老女人黄片| 国产av码专区亚洲av| 少妇的逼好多水| 秋霞在线观看毛片| 伦精品一区二区三区| 看非洲黑人一级黄片| 亚洲欧美日韩卡通动漫| 天天影视国产精品| 日本与韩国留学比较| 久久99热这里只频精品6学生| 成年av动漫网址| 精品国产露脸久久av麻豆| 免费不卡的大黄色大毛片视频在线观看| videossex国产| 在线观看一区二区三区激情| 久久99一区二区三区| 国国产精品蜜臀av免费| 久久久久久久大尺度免费视频| 人人妻人人澡人人爽人人夜夜| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久精品电影小说| 热99久久久久精品小说推荐| 亚洲av成人精品一二三区| 亚洲精品国产av成人精品| 国产xxxxx性猛交| 永久免费av网站大全| 精品一区在线观看国产| 亚洲精品中文字幕在线视频| 亚洲,欧美精品.| 永久免费av网站大全| av免费观看日本| 日韩成人伦理影院| 日本黄大片高清| 日韩大片免费观看网站| 黑人高潮一二区| 久久久久久久精品精品| 黄色毛片三级朝国网站| 妹子高潮喷水视频| 天天躁夜夜躁狠狠躁躁| 国产不卡av网站在线观看| 看十八女毛片水多多多| 熟女电影av网| 黄片无遮挡物在线观看| 国产高清国产精品国产三级| 99久久综合免费| 人妻少妇偷人精品九色| 国产一级毛片在线| 欧美变态另类bdsm刘玥| 亚洲成av片中文字幕在线观看 | 久久99一区二区三区| 一区在线观看完整版| 激情五月婷婷亚洲| www日本在线高清视频| 赤兔流量卡办理| 五月伊人婷婷丁香| 99久久综合免费| 午夜福利乱码中文字幕| 欧美人与性动交α欧美精品济南到 | 久久久久久久亚洲中文字幕| 蜜臀久久99精品久久宅男| 免费观看无遮挡的男女| 午夜影院在线不卡| 亚洲国产日韩一区二区| 亚洲av男天堂| 侵犯人妻中文字幕一二三四区| 大香蕉久久网| 国产一区亚洲一区在线观看| 国产精品久久久久久久久免| 亚洲国产精品一区二区三区在线| 人人澡人人妻人| 丝瓜视频免费看黄片| 母亲3免费完整高清在线观看 | 欧美日韩一区二区视频在线观看视频在线| 丰满乱子伦码专区| 熟妇人妻不卡中文字幕| 午夜日本视频在线| 18禁在线无遮挡免费观看视频| 国产色爽女视频免费观看| 香蕉丝袜av| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边做爽爽视频免费| 90打野战视频偷拍视频| 一级片'在线观看视频| 久久99热6这里只有精品| 免费高清在线观看日韩| 夫妻午夜视频| 久久精品国产亚洲av天美| 精品福利永久在线观看| 一区二区三区四区激情视频| 少妇精品久久久久久久| 日韩人妻精品一区2区三区| 亚洲国产精品一区三区| 久久国产精品男人的天堂亚洲 | 曰老女人黄片| 国产在线免费精品| 一级a做视频免费观看| 中文字幕亚洲精品专区| 日韩人妻精品一区2区三区| 日日摸夜夜添夜夜爱| 又黄又爽又刺激的免费视频.| 汤姆久久久久久久影院中文字幕| 国产精品久久久久久av不卡| 热99久久久久精品小说推荐| 久久这里只有精品19| av在线播放精品| 少妇精品久久久久久久| 汤姆久久久久久久影院中文字幕| av一本久久久久| 热re99久久精品国产66热6| 日韩,欧美,国产一区二区三区| 国产成人91sexporn| 亚洲人与动物交配视频| 大香蕉久久网| 一级毛片电影观看| 亚洲精品一区蜜桃| 亚洲美女黄色视频免费看| 日韩欧美精品免费久久| 欧美激情 高清一区二区三区| 一级毛片 在线播放| 国产乱来视频区| 在线天堂最新版资源| av天堂久久9| 在线免费观看不下载黄p国产| 人人妻人人澡人人爽人人夜夜| 午夜久久久在线观看| 18禁观看日本| 日本与韩国留学比较| 久久久精品区二区三区| 街头女战士在线观看网站| 韩国高清视频一区二区三区| 亚洲天堂av无毛| 久热久热在线精品观看| 观看av在线不卡| 成人二区视频| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 少妇人妻 视频| 18禁观看日本| www.色视频.com| 高清毛片免费看| 国产欧美日韩综合在线一区二区| 九九爱精品视频在线观看| 亚洲美女视频黄频| 国产乱人偷精品视频| 日日撸夜夜添| 久久国产精品男人的天堂亚洲 | 国产淫语在线视频| 亚洲欧洲日产国产| 日韩在线高清观看一区二区三区| 亚洲情色 制服丝袜| 99re6热这里在线精品视频| 天美传媒精品一区二区| 黄色一级大片看看| 免费日韩欧美在线观看| 午夜老司机福利剧场| 久久国产精品大桥未久av| 热re99久久国产66热| 亚洲人成77777在线视频| 成人无遮挡网站| 亚洲av电影在线进入| 咕卡用的链子| 最近最新中文字幕大全免费视频 | 精品熟女少妇av免费看| 老司机影院成人| 成人毛片a级毛片在线播放| 人妻少妇偷人精品九色| 只有这里有精品99| 人妻系列 视频| 精品亚洲成a人片在线观看| 18禁裸乳无遮挡动漫免费视频| 婷婷色综合www| 狠狠婷婷综合久久久久久88av| 最近中文字幕高清免费大全6| 久久精品aⅴ一区二区三区四区 | 日韩av不卡免费在线播放| 一区二区三区乱码不卡18| 日本黄色日本黄色录像| 亚洲国产精品一区三区| 午夜免费男女啪啪视频观看| 日韩大片免费观看网站| 99久久人妻综合| 久久这里只有精品19| 国产成人欧美| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 最近中文字幕2019免费版| 免费大片18禁| 免费少妇av软件| 精品久久久久久电影网| 香蕉精品网在线| 亚洲精品自拍成人| 欧美精品亚洲一区二区| 国产成人精品无人区| 少妇猛男粗大的猛烈进出视频| 亚洲经典国产精华液单| 观看美女的网站| 久久精品aⅴ一区二区三区四区 | 成人免费观看视频高清| 91久久精品国产一区二区三区| 国产色爽女视频免费观看| 久久这里只有精品19| 女人被躁到高潮嗷嗷叫费观| 亚洲熟女精品中文字幕| 啦啦啦啦在线视频资源| 欧美最新免费一区二区三区| 国产精品嫩草影院av在线观看| 大陆偷拍与自拍| 18禁在线无遮挡免费观看视频| 欧美人与善性xxx| 91精品伊人久久大香线蕉| 精品人妻熟女毛片av久久网站| av一本久久久久| 久久精品熟女亚洲av麻豆精品| 桃花免费在线播放| 新久久久久国产一级毛片| av福利片在线| 国产一级毛片在线| 日本与韩国留学比较| a级毛色黄片| 日韩av在线免费看完整版不卡| 久久久久网色| 日本黄大片高清| 中文字幕精品免费在线观看视频 | 亚洲婷婷狠狠爱综合网| 看十八女毛片水多多多| 91精品伊人久久大香线蕉| 十八禁网站网址无遮挡| 黑人欧美特级aaaaaa片| 亚洲精华国产精华液的使用体验| 有码 亚洲区| 久久精品国产亚洲av涩爱| 97在线视频观看| 亚洲丝袜综合中文字幕| 大码成人一级视频| av在线老鸭窝| 99热网站在线观看| 久久久久国产精品人妻一区二区| 久久免费观看电影| 一区二区三区乱码不卡18| 免费高清在线观看日韩| av国产久精品久网站免费入址| 国产在线免费精品| 亚洲欧美日韩卡通动漫| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 亚洲性久久影院| 久久热在线av| 日韩成人伦理影院| 久久久久人妻精品一区果冻| 日韩成人伦理影院| 久久人人爽av亚洲精品天堂| 欧美激情 高清一区二区三区| 国产片内射在线| 80岁老熟妇乱子伦牲交| www.熟女人妻精品国产 | 汤姆久久久久久久影院中文字幕| 亚洲五月色婷婷综合| 久久免费观看电影| 人妻少妇偷人精品九色| 男的添女的下面高潮视频| 五月天丁香电影| 插逼视频在线观看| 9191精品国产免费久久| 国产免费现黄频在线看| 卡戴珊不雅视频在线播放| 日本黄色日本黄色录像| 亚洲精品久久午夜乱码| 在线精品无人区一区二区三| 丰满乱子伦码专区| 久久精品aⅴ一区二区三区四区 | a级片在线免费高清观看视频| 多毛熟女@视频| 成年动漫av网址| 啦啦啦在线观看免费高清www| 国产又爽黄色视频| 国产精品人妻久久久影院| 国产一区有黄有色的免费视频| 国产成人免费观看mmmm| 国产乱人偷精品视频| av片东京热男人的天堂| 日韩制服丝袜自拍偷拍| 男人操女人黄网站| 热99国产精品久久久久久7| 国产福利在线免费观看视频| 在线 av 中文字幕| 成人亚洲欧美一区二区av| 国语对白做爰xxxⅹ性视频网站| 欧美最新免费一区二区三区| 美女内射精品一级片tv| 欧美日韩综合久久久久久| 卡戴珊不雅视频在线播放| www.av在线官网国产| 亚洲欧美日韩卡通动漫| 国产免费现黄频在线看| 国产免费一区二区三区四区乱码| 国产av精品麻豆| 欧美人与性动交α欧美软件 | 国产成人精品婷婷| 国产又爽黄色视频| 2018国产大陆天天弄谢| 在线观看免费视频网站a站| 99热6这里只有精品| 少妇人妻 视频| 精品人妻一区二区三区麻豆| 十分钟在线观看高清视频www| 国产亚洲精品第一综合不卡 | 国产成人免费无遮挡视频| 国产精品久久久久久久电影| 欧美丝袜亚洲另类| videosex国产| 咕卡用的链子| 久久99一区二区三区| 国产成人欧美| 九色成人免费人妻av| 欧美日本中文国产一区发布| a级毛片黄视频| 久久精品久久久久久噜噜老黄| 日韩电影二区| 欧美国产精品一级二级三级| 成人亚洲精品一区在线观看| 日韩 亚洲 欧美在线| 国产国语露脸激情在线看| 人妻 亚洲 视频| 久久 成人 亚洲| 天天操日日干夜夜撸| 国产爽快片一区二区三区| 免费av中文字幕在线| 亚洲人成77777在线视频| 日韩一区二区三区影片| 国产成人免费观看mmmm| 欧美激情国产日韩精品一区| 成年人午夜在线观看视频| av福利片在线| 少妇的逼水好多| 好男人视频免费观看在线| freevideosex欧美| 精品亚洲成a人片在线观看| 黄片播放在线免费| 亚洲一码二码三码区别大吗| 自线自在国产av| 国产片内射在线| 香蕉丝袜av| 在线观看www视频免费| av片东京热男人的天堂| 国产毛片在线视频| www日本在线高清视频| 日本色播在线视频| 精品国产一区二区三区四区第35| 久久 成人 亚洲| 2022亚洲国产成人精品| 丁香六月天网| 国产成人av激情在线播放| 人成视频在线观看免费观看| 久久久久久久国产电影| 国产亚洲午夜精品一区二区久久| 国产一区二区三区综合在线观看 | 久久免费观看电影| 新久久久久国产一级毛片| 久久久久久久精品精品| 天天躁夜夜躁狠狠久久av| 中文字幕免费在线视频6| 欧美少妇被猛烈插入视频| 你懂的网址亚洲精品在线观看| 欧美成人精品欧美一级黄| 久久久精品94久久精品| 久久女婷五月综合色啪小说| 秋霞在线观看毛片| 大码成人一级视频| www.av在线官网国产| 久久精品夜色国产| 亚洲精品乱久久久久久| 日本黄色日本黄色录像| 最近的中文字幕免费完整| 国产成人91sexporn| 久久久久精品人妻al黑| 熟女人妻精品中文字幕| 亚洲av.av天堂| 免费大片黄手机在线观看| 在线免费观看不下载黄p国产| 日韩三级伦理在线观看|