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

    近海岸強風風場特性現(xiàn)場實測研究

    2017-04-28 04:52:49胡尚瑜聶功恒李秋勝
    空氣動力學(xué)學(xué)報 2017年2期
    關(guān)鍵詞:季風陣風風場

    胡尚瑜, 聶功恒, 李秋勝, 張 敏

    (1. 桂林理工大學(xué) 廣西建筑新能源與節(jié)能重點實驗室, 廣西 桂林 541004;2. 香港城市大學(xué) 土木及建筑工程系, 中國 香港 999077)

    ?

    近海岸強風風場特性現(xiàn)場實測研究

    胡尚瑜1,*, 聶功恒1, 李秋勝2, 張 敏1

    (1. 桂林理工大學(xué) 廣西建筑新能源與節(jié)能重點實驗室, 廣西 桂林 541004;2. 香港城市大學(xué) 土木及建筑工程系, 中國 香港 999077)

    基于近海岸100 m測風塔獲取的近海面地貌條件5個不同高度風速數(shù)據(jù)。以10 m標準高度平均風速大于10 m/s為強風標準,選取臺風和季風風場強風樣本,應(yīng)用Bootstrap統(tǒng)計分析法,對強風條件下近地層平均風剖面、湍流強度剖面、湍流特征參數(shù)進行分析,探討了臺風和季風條件下這些參數(shù)的差異,并將實測結(jié)果與風荷載規(guī)范規(guī)定值進行比較,檢驗風載荷規(guī)范的適宜性。分析結(jié)果表明 :近地邊界層的平均風剖面符合指數(shù)律或?qū)?shù)律,臺風和季風風剖面指數(shù)分別約為0.12和0.06;近地平均湍流度剖面符合指數(shù)律,臺風風場各高度平均湍流度與現(xiàn)行規(guī)范相接近;臺風風場各高度陣風因子均值較風荷載規(guī)范推算值要大,但季風風場陣風因子均值與規(guī)范推算值相接近。

    風荷載;平均風剖面;湍流度;陣風因子;Bootstrap統(tǒng)計

    0 引 言

    近年來近海岸建筑和高聳結(jié)構(gòu)及超大跨度橋梁等工程不斷興建,近海岸工程的結(jié)構(gòu)抗風設(shè)計和安全可靠性顯得尤為重要。在強/臺風極端條件下,近地風場特性和結(jié)構(gòu)極值風荷載的現(xiàn)場實測研究,對近海岸工程結(jié)構(gòu)抗風設(shè)計和結(jié)構(gòu)安全評估有著重要的意義。近年來國內(nèi)外學(xué)者[1-10]通過移動測風塔和近海岸觀測塔等各種測試平臺,獲取了大量的近地臺風風場數(shù)據(jù)?;诂F(xiàn)場實測數(shù)據(jù),統(tǒng)計分析了近地平均風剖面指數(shù)、湍流度、陣風因子等物理量。由于受臺風自身湍流結(jié)構(gòu)及隨機性,不同臺風區(qū)域和不同強風標準及統(tǒng)計方法等因素影響,近地臺風風場湍流特性參數(shù)取值存在較大差異。

    本文基于近海岸100 m測風塔在2009~2011年期間測得的10 m、50 m、65m、80 m、100 m 5個不同高度實測風速數(shù)據(jù),以10 m高度平均風速大于10 m/s為強風標準,選取有效分析風速樣本,應(yīng)用Bootstrap統(tǒng)計方法,對近海岸海面來流地貌條件,臺風和季風風場的平均風剖面、湍流強度剖面及相關(guān)湍流特征參數(shù)進行分析,并探討了臺風和季風風場近地層的平均風剖面指數(shù)、湍流度和陣風因子等參數(shù)的差異,同時將實測結(jié)果與風荷載規(guī)范規(guī)定值進行比較,檢驗風荷載規(guī)范的適宜性。本文研究結(jié)果可為近海岸工程結(jié)構(gòu)抗風設(shè)計風剖面參數(shù)取值提供參考。

    1 近地邊界層風特性

    1.1 平均風剖面

    在強風作用下,近地中性穩(wěn)定大氣邊界層的湍流主要由地面粗糙度控制,平均風速剖面變化規(guī)律常認為符合指數(shù)律或?qū)?shù)律[11-12]。

    指數(shù)律模型:

    式中:zR為參考高度;α為冪指數(shù),其值與地貌狀況相關(guān)。

    對數(shù)律模型:

    式中:k為卡曼常數(shù),取0.4;u*為摩擦速度;z0為地面粗糙度長度。

    1.1.1 指數(shù)律模型擬合

    依據(jù)指數(shù)律模型式可得:

    各高度擬合計算風速與實測風速離差平方和計算式如下:

    式中:U(z)為實測值,U(zi)為擬合計算風速。

    各高度擬合計算風速與實測風速離差平方和均方根誤差計算式如下:

    擬合風剖面的相對誤差大小定義為離差平方和均方根誤差與參考風速大小的比值如下:

    1.1.2 對數(shù)律模型擬合

    依據(jù)對數(shù)律模型計算式可得:

    應(yīng)用最小二乘法對各高度實測平均風速進行線性擬合,當擬合計算風速與實測風速離差平方和趨于零,可得地面粗糙度z0的計算公式[13]如下:

    摩擦速度定義為:

    式中:τ0為地表的剪切應(yīng)力,ρ為空氣密度。本文摩擦速度依據(jù)地面粗糙度z0擬合值和不同高度的風速代入公式(2),運用最小二乘法進行線性擬合,可得其對應(yīng)的摩擦速度u*擬合值。

    1.2 湍流強度

    湍流度反映了風的脈動強度,是確定脈動風荷載的關(guān)鍵參數(shù)之一。湍流度定義為平均時距的脈動風速標準差與水平方向平均風速的比值:

    式中σu為1s采樣頻率脈動風速的標準差,文中的湍流度主要反映大尺度湍流脈動的特征。依據(jù)文獻[14]比較超聲測風儀和NRG杯式測風儀實測湍流強度分析,認為兩個風速儀獲得的湍流強度具有相同趨勢,超聲風速儀的湍流強度略大于NRG杯式測風儀,因而本文對NRG杯式測風儀實測各高度的湍流強度乘以1.1進行了相應(yīng)的修訂。

    GB50009-2012建筑結(jié)構(gòu)荷載規(guī)范[15]給出了各高度平均湍流強度剖面計算式:

    式中α為風剖面指數(shù),其值與地貌狀況相關(guān)。對應(yīng)于A類地貌,規(guī)范給定值為0.12。I10為10m高度名義湍流度,對應(yīng)A類地貌取值為0.12。

    1.3 陣風因子

    風的脈動強度也可用陣風因子來表示。陣風因子定義為陣風持續(xù)期的平均風速最大值與平均時距內(nèi)平均風速的比值:

    本文取陣風持續(xù)期tg=3 s,平均時距T=10 min。

    陣風峰值因子定義為:

    式中:Umax,tg為陣風時距內(nèi)最大平均風速。

    陣風因子和湍流度在物理上都表示風速的脈動強度,相互關(guān)系可表示為:

    依據(jù)式湍流強度剖面計算式(12)和式(15)可得各高度的陣風因子剖面計算式:

    1.4 Boostrap法分析思路

    以選取梯度風各高度范圍內(nèi)的強風風速樣本,采用指數(shù)律模型擬合得到各樣本的風剖面指數(shù)計算序列(αs1,αs2,…,αsn)為例。運用Bootstrap方法,統(tǒng)計分析風剖面指數(shù)置信區(qū)間的分析思路如下:

    4) 按步驟3)連續(xù)重復(fù)B次,即得到B個Bootstrap再生樣本,用矩陣表示為:

    計算再生樣本的統(tǒng)計量:

    其中:y″i=y″i1,y″i2,…,y″iB,i=1,2。

    5) 根據(jù)再生樣本y″i,可以得到μ和σ2的估計:

    2 現(xiàn)場實測及數(shù)據(jù)選取

    2.1 實測系統(tǒng)及場地條件

    測風塔位于海南文昌木蘭頭,東經(jīng)110°44′,北緯20°09′,在離地面10 m、50 m、65 m、80 m、100 m高度處布設(shè)安裝了NRG-Symphonie型杯式風速儀進行風速觀測,同時在50 m、65 m、100 m高度處安裝了風向儀進行風向角的觀測。測試系統(tǒng)采樣頻率為1 Hz,每10 min輸出1組平均風速風向、極大風速風向和風速、風向標準差等。觀測塔觀測地貌狀況:東北方位來流地貌條件為海面,距離海岸線約為50 m;西南方位為平坦開闊場地,下墊面為近地有2~3 m高的馬尾松防護林。

    2.2 強風樣本選取

    2009年7月至2011年7月兩年期間,近海岸梯度風觀測塔獲取的近地臺風風場及季風風場有效觀測樣本的極值風速統(tǒng)計如表1所示,臺風“彩虹” 實測3 s最大平均風速為26.2 m/s,10 min最大平均風速為20.4 m/s。臺風“凱薩娜” 實測3 s最大平均風速為26.2 m/s, 10 min最大平均風速為20.3 m/s。臺風“芭瑪”實測3 s最大平均風速為20.1 m/s,10 min最大平均風速為14.7 m/s。

    表1 觀測樣本資料 Table 1 Observation of wind records used in study

    本文以測風塔10 m高度處10 min平均風速大于10 m/s為強風標準,選取臺風和季風風場的強風樣本數(shù)據(jù),對近地邊界層風特性進行分析。10 m高度處臺風和季風風場強風風速、風向角樣本分布如圖1(a、b)所示:臺風和季風風場強風樣本主要分布在30°~60°范圍內(nèi),因而本文選取風向角分布在30°~60°范圍內(nèi)的強風樣本,該方位來流地貌為近海面,考慮來流經(jīng)歷由海面地貌到近海岸地貌條件,存在內(nèi)部邊界層高度需確定問題。參考文獻[17]給出的內(nèi)部邊界高度計算公式,根據(jù)測風塔距其海岸線大約為50 m,假設(shè)海面地貌粗糙度長度由0.002 m過渡到0.02 m地貌條件時,其內(nèi)部邊界層高度值約為7.1 m。因而本文實測各高度風速實測大小反映來流為海面地貌條件下的邊界層風速,旨在分析近海岸海面地貌狀況下臺風風場和季風風場特性。各觀測層的高度應(yīng)從海平面開始,根據(jù)實際量測測風塔相對海面高度差約為4 m。同時考慮到強風條件下受海浪的影響,測風塔平面位置相對海面高度差近似取為3 m,即各觀測層的高度分別為13 m、53 m、68 m、83 m、103 m。

    (a) 臺風樣本

    (b) 季風樣本

    3 風特性分析

    3.1 平均風特性

    選取來流為近海面地貌條件下,實測臺風“彩虹”、“凱薩娜”、“芭瑪”及季風風場各觀測高度的平均風速數(shù)據(jù),應(yīng)用對數(shù)律模型和指數(shù)律模型對風剖面進行最小二乘線性擬合。考慮當10 m高度風速達到17.2 m/s時,風力等級將達到八級大風級。因而以風速17.2 m/s為分界線對臺風“凱薩娜”風速劃分區(qū)間,區(qū)間[10,17.2)為“凱薩娜”(I),區(qū)間[17.2,20.3)為“凱薩娜”(II)。指數(shù)律、對數(shù)律風剖面參數(shù)擬合值如表2所示:臺風“彩虹”風剖面指數(shù)α均值為0.122,應(yīng)用Boostrap法計算得到的風剖面指數(shù)α均值在95%置信度條件下的置信區(qū)間為[0.115,0.130];臺風 “凱薩娜”(I) 平均風速小于17.2 m/s樣本的風剖面指數(shù)α均值為0.128;臺風 “凱薩娜”(II) 平均風速大于等于17.2 m/s大風樣本的風剖面指數(shù)α均值為0.119;“芭瑪”風場風剖面α指數(shù)均值范圍為0.118,各臺風風場與規(guī)范GB50009-2012中 A類地貌條件風剖面指數(shù)α規(guī)定值0.120相近。季風風場風剖面指數(shù)α均值為0.062,其95%置信區(qū)間為[0.060,0.063]。季風風場風剖面指數(shù)α均值小于規(guī)范值和臺風風場值。臺風風場地面粗糙度長度和摩擦速度與平均風速大小相關(guān),與季風風場相比,地面粗糙度長度和摩擦速度相對增大。近地100 m 高度范圍內(nèi)的指數(shù)律和對數(shù)律剖面以及相關(guān)擬合曲線如圖2所示。由圖2(a)可知:“彩虹”臺風近地層平均風剖面各高度平均風速變化規(guī)律符合指數(shù)律和對數(shù)律模型。同理由圖2(b-e)可知,指數(shù)律和對數(shù)律模型對臺風“凱薩娜(I)、(II)”和 “芭瑪”及季風風場近地層平均風剖面各高度風速擬合基本較好,各臺風指數(shù)律風剖面擬合相對誤差值與參考風速變化關(guān)系如圖2(f)所示,整體上相對誤差約為2%左右,表明強風條件下指數(shù)律和對數(shù)律模型都能較好的反映近地100 m高度范圍內(nèi)的平均風速變化規(guī)律。值得說明的是,在50~70 m 高度范圍,近地層存在一段風速變化大小變化不明顯區(qū)間,這與文獻[18]在近海平坦地貌條件,海面來風時,在50~70 m高度范圍,風速大小變化不明顯,具有相似的現(xiàn)象,在高度60 m附近范圍存在強風區(qū)段。其主要原因仍需要通過采用超聲風速儀等資料分析其區(qū)間梯度風變化規(guī)律及差異產(chǎn)生的原因。

    表2 平均風剖面參數(shù)Table 2 Parameters of mean wind speed

    (a) 臺風“彩虹”

    (b) 臺風“凱薩娜”(I)

    (c) 臺風“凱薩娜”(II)

    (d) 臺風“芭瑪”

    (e) 季風風場

    (f) 相對誤差

    3.2 湍流特性

    3.2.1 湍流強度及剖面

    選取來流為近海面地貌條件下,實測臺風“彩虹”、“凱薩娜”、“芭瑪”及季風風場各觀測高度的平均風速及風速標準差的強風樣本,統(tǒng)計分析了近海岸地貌下,臺風和季風風場近地邊界層湍流度分布規(guī)律。如表3所示:“彩虹”、“凱薩娜”樣本(I)、(II)和“芭瑪”臺風風場13 m高度湍流強度均值分別為0.119、0.135、0.123和0.125,均值范圍為0.119~0.135,略大于現(xiàn)行規(guī)范A類地貌規(guī)范計算值0.116。季風風場13 m高度湍流強度均值為0.104,略小于現(xiàn)行規(guī)范A類地貌規(guī)范計算值。

    如圖3所示,臺風風場平均湍流度沿高度變化規(guī)律符合指數(shù)律分布,各臺風湍流強度平均剖面與現(xiàn)行規(guī)范給定值相接近。應(yīng)用公式(12)和實測風剖面指數(shù)z0擬合值,采用最小二乘法線性擬合,可得“彩虹”、“凱薩娜”樣本(I)、(II)和 “芭瑪”臺風風場13m高度湍流度擬合值。同時考慮到13m高度的湍流實測值可能受到近地面干擾影響,在湍流特性參數(shù)擬合計算時未采用13m高度湍流實測值。擬合所得各臺風13m高度湍流度擬合值分別為0.116、0.114、0.108和0.109,與規(guī)范給定計算值0.116較為接近,小于實測值。各高度的規(guī)范推薦值能較好反映實測各臺風風場平均湍流度剖面均值分布特性。從圖3(e)可以看出,季風實測各高度湍流度明顯小于規(guī)范推薦值。同理對季風風場湍流剖面進行擬合,可得季風風場13 m高度湍流度擬合值為0.078,小于13 m高度實測湍流度均值0.104和規(guī)范計算值0.116。近地13 m高度的湍流度的擬合值與實測值差別的主要原因是:通常湍流參數(shù)的內(nèi)邊界層高度是平均風速的2倍,相對近地13 m高以下大氣湍流的更易受近地植被擾動作用影響。

    表3 各高度的湍流強度與陣風因子統(tǒng)計Table 3 Parameter of wind turbulence intensity and gust factors statistics

    (a) 臺風“彩虹”

    (c) 臺風“凱薩娜”(II)

    (b) 臺風“凱薩娜”(I)

    (d) 臺風“芭瑪”

    (e) 季風風場

    3.2.2 陣風因子及剖面

    基于13 m、53 m、68 m、83 m、103 m 5個不同高度處相對應(yīng)的平均風速及3 s最大陣風風速,統(tǒng)計分析臺風和季風風場近地邊界層陣風因子及剖面。如表3所示:臺風 “彩虹”、“凱薩娜(I)”、“凱薩娜(II)”和 “芭瑪”在13 m高度處的陣風因子均值分別為1.33、1.38、1.36、1.37,略大于現(xiàn)行規(guī)范A類地貌規(guī)范計算值1.29。季風風場13 m高度陣風因子均值為1.30,與現(xiàn)行規(guī)范A類地貌規(guī)范計算值接近。如圖4所示,實測臺風和季風風場陣風因子沿高度變化規(guī)律同樣符合指數(shù)律分布,臺風風場陣風因子剖面大于現(xiàn)行規(guī)范給定值,季風陣風因子與規(guī)范推算值相接近。依據(jù)陣風因子計算式(16),采用實測風剖面指數(shù)α擬合值和各臺風風場湍流強度剖面擬合13 m高度湍流度擬合值即I13取0.116、0.114、0.108和0.109,對各臺風風場陣風因子剖面進行最小二乘線性擬合。擬合“彩虹”、“凱薩娜”(I)、(II)和 “芭瑪”臺風風場的峰值因子g值分別為2.74、2.68、2.73和2.84,略大于規(guī)范給定值2.5。同理對季風風場陣風因子剖面進行擬合得季風峰值因子g值為3.07,大于規(guī)范給定值。相對而言,近地13m高度的陣風因子的擬合值與實測值差別的主要原因是受近地面干擾作用影響產(chǎn)生一定大尺度陣風湍流。

    (a) 臺風“彩虹”

    (b) 臺風“凱薩娜”(I)

    (c) 臺風“凱薩娜”(II)

    (d) 臺風“芭瑪”

    (e) 季風

    3.3 風特性參數(shù)與平均風速

    選取臺風“凱薩娜”(I)、(II)風場樣本統(tǒng)計的風特性參數(shù),探討風特性參數(shù)與平均風速的相關(guān)性。臺風風場風剖面指數(shù)與平均風速的關(guān)系如圖5所示:臺風風剖面指數(shù)α隨平均風速增大有減小趨勢。類似對數(shù)律擬合得臺風風場摩擦速度和地面粗糙度長度與平均風速大小相關(guān),如圖6、圖7所示,臺風風場地面粗糙度長度隨平均風速增大而減少,而摩擦速度隨平均風速增大而增大。

    臺風“凱薩娜”(I)、(II)風場樣本的相對海面高度為13 m處的湍流強度與平均風速的關(guān)系如圖8所均風速增大而減少,而在平均風速17.2~22 m/s區(qū)間范圍相對趨于平穩(wěn)。

    圖5 風剖面指數(shù)與平均風速的關(guān)系Fig.5 Variation of power law index with mean wind speed

    圖6 摩擦速度與平均風速的關(guān)系Fig.6 Variation of friction with mean wind speed

    示,在平均風速10~22 m/s范圍內(nèi),臺風的湍流強度隨平均風速增大而減小,隨平均風速的進一步增大而相對趨于穩(wěn)定。相對海面高度為13 m處的陣風因子與平均風速的關(guān)系如圖9所示,臺風風場在平均風速10~22 m/s范圍內(nèi),陣風因子隨平均風速的變化相對穩(wěn)定。綜上所述,近海岸地貌強風條件下,臺風風場風剖面指數(shù)、湍流強度和陣風因子大于季風風場實測值。各臺風風場湍流參數(shù)隨著平均風速增大有減小的趨勢,隨著風速的繼續(xù)增大而趨于穩(wěn)定。與文獻[6-7]具有相似的結(jié)論,各臺風風場湍流參數(shù)隨著平均風速增大而減少,而在平均風速17.2~22 m/s區(qū)間范圍相對趨于平穩(wěn)。

    圖7 地貌粗糙度與平均風速的關(guān)系Fig.7 Variation of roughness with mean wind speed

    圖9 陣風因子與平均風速的關(guān)系Fig.9 Variation of gust factor with mean wind speed

    4 結(jié) 論

    基于近地層風場風特性的實測強風風速樣本,對近海岸地貌強風條件下臺風和季風風場參數(shù)如風剖面指數(shù)、湍流度和陣風因子等進行了統(tǒng)計分析,得到如下結(jié)論:

    1) 強風風場近地平均風速剖面符合指數(shù)律或?qū)?shù)律,臺風風場風剖面指數(shù)α均值約為0.12,與規(guī)范規(guī)定值接近。季風風場風剖面指數(shù)α均值為0.06,小于規(guī)范值。臺風風場地面粗糙度長度及摩擦速度較季風風場值偏大達70%;

    2) 臺風風場平均湍流度剖面符合指數(shù)律分布,臺風湍流強度平均剖面與規(guī)范A類地貌湍流強度推算值相一致。季風風場各高度實測湍流度明顯小于規(guī)范給定值;

    3) 臺風風場陣風因子剖面變化規(guī)律同樣符合相應(yīng)的平均風速剖面α指數(shù)律分布。臺風風場陣風因子剖面大于現(xiàn)行規(guī)范給定值,季風風場陣風因子與規(guī)范給定值接近;

    4) 臺風風場平均風剖面α指數(shù)、湍流度和陣風因子大于季風風場實測值。各臺風風場湍流參數(shù)隨著平均風速增大而減少,而在平均風速17.2~22 m/s區(qū)間范圍相對趨于平穩(wěn)。

    討論:本文依據(jù)近海岸觀測塔獲取的臺風和季風強風特性的實測數(shù)據(jù)進行統(tǒng)計分析,得出了近海岸海面地貌條件下近地平均風速剖面、湍流度剖面及陣風因子剖面的變化規(guī)律與現(xiàn)行規(guī)范A類地貌條件推薦值具有較好的一致性,但具體量值存在的差異,風速剖面參數(shù)化仍需要更多的有效實測強風樣本確定和驗證。

    致謝:感謝中國華能海南文昌風電廠提供的實測風速數(shù)據(jù)。

    [1]Schroeder J L, Edwards B P, Giammanco I M. Observed tropical cyclone wind flow characteristics[J]. Wind and Structures, 2009, 12(4): 349-381.

    [2]Masters F J, Tieleman H W, Balderrama J A. Surface wind measurements in three Gulf Coast hurricanes of 2005[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2010, 98(10): 533-547.

    [3]Fujimura M, Maeda J. Cross-correlation of fluctuating components of wind speed based on strong wind measurement[C]//Proceedings of 7th Asia-Pacific Conference on Wind Engineering, T1-D, Taipei. 2009.

    [4]Cao S, Tamura Y, Kikuchi N, et al. Wind characteristics of a strong typhoon[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2009, 97(1): 11-21.

    [5]Zheng Q M, Luo Y L. Investigations on the terrain classification and characteristics of atmospheric turbulent boundary layer[J]. Journal of Architecture, 2009(68_S): 103-120. 鄭啟明, 羅元隆. 風力規(guī)范之標準地況分類與紊流邊界層特性研究[J]. 建筑學(xué)報, 2009(68_S): 103-120.

    [6]Hu S Y, Song L L, L Q S. Monitoring of typhoons in surface boundary layer and analysis of wind turbulence characteristics[J]. Journal of Building structures, 2011, 32(4): 1-8. 胡尚瑜, 宋麗莉, 李秋勝. 近地邊界層臺風觀測及湍流特征參數(shù)分析[J]. 建筑結(jié)構(gòu)學(xué)報, 2011, 32(4): 1-8.

    [7]Xiao Y Q, Li L X, Song L L, et al. Study on wind characteristics of typhoonhagupit based on offshore sea surface measurements[J]. Acta Aerodynamica Sinica, 2012, 30(3): 380-389. 肖儀清, 李利孝, 宋麗莉, 等. 基于近海海面觀測的臺風黑格比風特性研究[J]. 空氣動力學(xué)學(xué)報, 2012, 30(3): 380-389.

    [8]Song L l, Mao H Q, Huang H H, et al. Analysis on boundary layer turbulent features of landfalling typhoon[J]. Acta Meteorologica Sinica, 2005, 63(6): 915-921. 宋麗莉, 毛慧琴, 黃浩輝, 等. 登陸臺風近地層湍流特征觀測分析[J]. 氣象學(xué)報, 2005, 63(6): 915-921.

    [9]Song L L, Pang J B, Jiang C L, et al. Field measurement and analysis of turbulence coherence for typhoon nuri at macao friendship bridge[J]. Science China Technological Sciences, 2010, 53(10): 2647-2657.

    [10]Li Q S, Dai Y M, Li Z N, et al. Surface layer wind field characteristics during a severe typhoon ‘Hagupit’ landfalling[J]. Journal of Building structures, 2010, (4): 54-61. 李秋勝, 戴益民, 李正農(nóng), 等. 強臺風“黑格比”登陸過程中近地風場特性[J]. 建筑結(jié)構(gòu)學(xué)報, 2010, (4): 54-61.

    [11]Powell M D, Vickery P J, Reinhold T A. Reduced drag coefficient for high wind speeds in tropical cyclones[J]. Nature, 2003, 422: 279-283.

    [12]Choi E C C. Characteristics of typhoons over the south China sea[J]. Journal of wind engineering and industrial aerodynamics, 1978, 3: 353-365.

    [13]Archer C L, Jacobson M Z. Spatial and temporal distributions of US winds and wind power at 80 m derived from measurements[J]. Journal of Geophysical Research: Atmospheres(1984-2012), 2003, 108(D9).

    [14]Song L, Li Q S, Chen W, et al. Wind characteristics of a strong typhoon in marine sur-face boundary layer[J]. Wind and Structures, 2012, 15(1): 1-15.

    [15]Load code for the design of building structures. GB50009-2012[S]. Beijing: China Architecture snd Building Press, 2012: 55. 建筑結(jié)構(gòu)荷載規(guī)范. GB50009-2012[S]. 北京: 中國建筑工業(yè)出版社, 2012: 55.

    [16]Bradley E. Bootstrap methods: another look at the Jackknife[J]. The Annals of Statistics, 1979, 7(1): 1-20.

    [17]Savelyev S A, Taylor P A. Internal boundary layers: I. Height formulae for neutral and diabatic flows[J]. Boundary-Layer Meteorology, 2005, 115(1): 1-25.

    [18]Li L X, Xiao Y Q, Song L L, et al. Study on wind profile of typhoon hagupit using wind observed tower and wind profile radar measurements[J]. Engineering Mechanics, 2012, 29(9): 284-293. 李利孝, 肖儀清, 宋麗莉, 等. 基于風觀測塔和風廓線雷達實測的強臺風黑格比風剖面研究[J]. 工程力學(xué), 2012, 29(9): 284-293.

    Field measurement study on wind characteristics of strong windstorm in coastal terrain

    Hu Shangyu1,*, Nie Gongheng1, Li Qiusheng2, Zhang Min1

    (1.GuangxiKeyLaboratoryofNewEnergyandBuildingEnergySaving,GuilinUniversityofTechnology,Guilin541004,China; 2.DepartmentofBuildingandConstruction,CityUniversityofHongKong,HongKong999077)

    Mean wind speed and wind turbulence characteristics of strong winds in surface boundary layer were studied by bootstrap method based on strong wind data sample in accordance with 10m height and mean wind speed higher than 10 m/s. The strong wind data were recorded from 100 m meteorological tower at sites near seashore with strong wind conditions and coastal terrain. The major objective of the paper is to further understand differences between near ground typhoon-generated and monsoon wind characteristics. Furthermore, the estimated wind profile and turbulence profiles were compared to those stipulated by wind loads Standard to verify the suitability of the wind load specification.The results reveal that in the near-surface range vertical distribution of mean wind speed can be well described by a logarithmic law and a power law. The mean values of exponent of the power-law profile are 0.12 and 0.06 for the coastal terrain in the typhoons and monsoon wind climates, respectively. The variation of the mean longitudinal turbulence intensity with varing height approximately follows a power law. The mean value of longitudinal turbulence intensity at five heights was fitting the results of the turbulence intensity profile. The profile is which was approximately same as that stipulated in wind loads Standard for coastal terrain exposure. Meanwhile, the measured gust factor of typhoons is larger than that stipulated in wind loads Standard for the coastal terrain exposure, but for the monsoon, it is close to the code value.

    wind loads; mean wind speed profile; turbulence intensity; gust factor; bootstrap statistics

    0258-1825(2017)02-0242-09

    2015-10-21;

    2015-12-08

    國家自然科學(xué)基金(51308140,51668015); 廣西自然科學(xué)基金(2015GXNSFAA139251);廣西建筑新能源與節(jié)能重點實驗室(16-J-21-7);廣西高校科研項目(2013YB103)

    胡尚瑜*(1981-),男,湖南永州人,工學(xué)博士,副教授, 研究方向: 結(jié)構(gòu)風工程. E-mail:hushangyu@glut.edu.cn

    胡尚瑜, 聶功恒, 李秋勝. 近海岸強風風場特性現(xiàn)場實測研究[J]. 空氣動力學(xué)學(xué)報, 2017, 35(2): 242-250.

    10.7638/kqdlxxb-2015.0009 Hu S Y, Nie G H, Li Q S. Field measurement studyon wind characteristics of strong windstorm in coastal terrain[J]. Acta Aerodynamica Sinica, 2017, 35(2): 242-250.

    TU312.1

    A doi: 10.7638/kqdlxxb-2015.0009

    猜你喜歡
    季風陣風風場
    基于FLUENT的下?lián)舯┝魅S風場建模
    陣風戰(zhàn)斗機
    軍事文摘(2021年19期)2021-10-10 13:29:36
    法國陣風戰(zhàn)斗機
    軍事文摘(2021年17期)2021-09-24 01:11:46
    綠水青山圖——海洋季風的贊歌
    陣風勁吹
    航空世界(2018年12期)2018-07-16 08:34:50
    “最美風場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    戶撒刀
    側(cè)向風場中無人機的飛行研究
    萬馬奔騰
    中華奇石(2015年9期)2015-07-09 18:32:58
    夾在兩陣風之間
    亚洲精品国产精品久久久不卡| 久久精品国产99精品国产亚洲性色| www.999成人在线观看| 亚洲18禁久久av| 波多野结衣高清作品| 精品久久久久久久毛片微露脸| 欧美性猛交╳xxx乱大交人| 狠狠狠狠99中文字幕| 丝袜人妻中文字幕| 中文字幕精品亚洲无线码一区| 精品久久蜜臀av无| 精品久久久久久久人妻蜜臀av| 特级一级黄色大片| 国产亚洲精品av在线| 99久久99久久久精品蜜桃| 成年版毛片免费区| 动漫黄色视频在线观看| 两个人的视频大全免费| 日韩精品免费视频一区二区三区| 琪琪午夜伦伦电影理论片6080| 亚洲一码二码三码区别大吗| 天天添夜夜摸| 狂野欧美白嫩少妇大欣赏| 成人三级做爰电影| 国产99久久九九免费精品| 亚洲国产日韩欧美精品在线观看 | 淫秽高清视频在线观看| 母亲3免费完整高清在线观看| 久久久久国内视频| 久久久国产成人精品二区| 欧美国产日韩亚洲一区| 中文字幕av在线有码专区| 欧美日韩精品网址| 国产精品98久久久久久宅男小说| 一区福利在线观看| 一进一出抽搐gif免费好疼| 欧美绝顶高潮抽搐喷水| 日本在线视频免费播放| 一进一出抽搐动态| 人成视频在线观看免费观看| a级毛片a级免费在线| 国产又色又爽无遮挡免费看| 精品福利观看| 亚洲av中文字字幕乱码综合| 色尼玛亚洲综合影院| 国产成人欧美在线观看| 一级毛片高清免费大全| 小说图片视频综合网站| 午夜福利在线在线| 亚洲成人国产一区在线观看| 国产精品精品国产色婷婷| 亚洲精品美女久久av网站| 宅男免费午夜| 日韩av在线大香蕉| 99riav亚洲国产免费| 岛国在线免费视频观看| 性色av乱码一区二区三区2| 国产黄色小视频在线观看| 亚洲熟妇中文字幕五十中出| 五月伊人婷婷丁香| 欧美av亚洲av综合av国产av| av有码第一页| 91国产中文字幕| 校园春色视频在线观看| 国产区一区二久久| 丁香六月欧美| 婷婷精品国产亚洲av在线| 美女黄网站色视频| 日韩av在线大香蕉| 老司机午夜福利在线观看视频| 久久人妻av系列| 国产亚洲av高清不卡| 午夜影院日韩av| 99国产极品粉嫩在线观看| 亚洲成人精品中文字幕电影| 日本成人三级电影网站| 精品免费久久久久久久清纯| 亚洲一区二区三区色噜噜| 一边摸一边做爽爽视频免费| 亚洲成人免费电影在线观看| 全区人妻精品视频| 日韩欧美 国产精品| 嫩草影视91久久| 男女午夜视频在线观看| 日韩欧美精品v在线| av福利片在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 88av欧美| 性色av乱码一区二区三区2| 国产亚洲精品第一综合不卡| 久久香蕉激情| 丝袜人妻中文字幕| 麻豆成人午夜福利视频| 亚洲欧美一区二区三区黑人| 免费人成视频x8x8入口观看| 欧美性猛交黑人性爽| 18禁黄网站禁片免费观看直播| 国产亚洲av高清不卡| 丁香六月欧美| 啦啦啦观看免费观看视频高清| 两个人看的免费小视频| 啦啦啦韩国在线观看视频| 九色国产91popny在线| 日韩国内少妇激情av| 日韩国内少妇激情av| 国产午夜福利久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产av不卡久久| 久久精品91无色码中文字幕| 国产精品亚洲av一区麻豆| 中文字幕久久专区| 曰老女人黄片| 女同久久另类99精品国产91| 可以在线观看毛片的网站| 中出人妻视频一区二区| 国产熟女午夜一区二区三区| 亚洲avbb在线观看| 色av中文字幕| 亚洲av中文字字幕乱码综合| 岛国在线观看网站| 999久久久精品免费观看国产| 美女午夜性视频免费| 女生性感内裤真人,穿戴方法视频| 人成视频在线观看免费观看| 欧美3d第一页| 欧美日本视频| 国产三级中文精品| 久久久久久久午夜电影| 禁无遮挡网站| 久久久久精品国产欧美久久久| 久久久久免费精品人妻一区二区| 成人av一区二区三区在线看| 一二三四社区在线视频社区8| 亚洲一码二码三码区别大吗| 久久久久性生活片| 国产精品一区二区免费欧美| 国产视频一区二区在线看| 男女午夜视频在线观看| 叶爱在线成人免费视频播放| 午夜老司机福利片| 国产精品一区二区免费欧美| 国产精品久久久久久久电影 | 一个人观看的视频www高清免费观看 | 国产熟女午夜一区二区三区| 久久国产精品人妻蜜桃| 精品少妇一区二区三区视频日本电影| 国产免费av片在线观看野外av| 叶爱在线成人免费视频播放| 夜夜夜夜夜久久久久| 窝窝影院91人妻| 国产亚洲精品久久久久久毛片| 午夜福利18| 麻豆国产97在线/欧美 | 美女扒开内裤让男人捅视频| 国产熟女xx| 老汉色∧v一级毛片| 少妇被粗大的猛进出69影院| 一级片免费观看大全| 夜夜夜夜夜久久久久| 看免费av毛片| 麻豆成人av在线观看| 欧美绝顶高潮抽搐喷水| 国产一区在线观看成人免费| 在线观看美女被高潮喷水网站 | 日韩欧美三级三区| 黄片大片在线免费观看| √禁漫天堂资源中文www| 757午夜福利合集在线观看| 亚洲精品一区av在线观看| 哪里可以看免费的av片| 亚洲精品久久国产高清桃花| 法律面前人人平等表现在哪些方面| 国产1区2区3区精品| 国产视频一区二区在线看| 999久久久国产精品视频| 波多野结衣高清无吗| 久久精品91蜜桃| 久久久久久免费高清国产稀缺| 亚洲国产精品999在线| 国产精品自产拍在线观看55亚洲| 热99re8久久精品国产| 天堂影院成人在线观看| 久久久久久免费高清国产稀缺| 久久这里只有精品中国| 麻豆av在线久日| 波多野结衣巨乳人妻| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 成人18禁在线播放| 99riav亚洲国产免费| 国产欧美日韩一区二区精品| 欧美成人免费av一区二区三区| 久99久视频精品免费| 久久天躁狠狠躁夜夜2o2o| 亚洲人成77777在线视频| 香蕉丝袜av| 国产精品爽爽va在线观看网站| 舔av片在线| videosex国产| 久久久久久久午夜电影| 欧美性猛交╳xxx乱大交人| 亚洲黑人精品在线| 国产亚洲av嫩草精品影院| 可以在线观看的亚洲视频| 狂野欧美激情性xxxx| 性欧美人与动物交配| 久久热在线av| 一级毛片精品| 亚洲人成伊人成综合网2020| 男女之事视频高清在线观看| 国内毛片毛片毛片毛片毛片| a级毛片a级免费在线| 日韩三级视频一区二区三区| 在线看三级毛片| 成人精品一区二区免费| av中文乱码字幕在线| 欧美午夜高清在线| 久久精品91无色码中文字幕| 欧美另类亚洲清纯唯美| 久久久久久久久免费视频了| 一区二区三区激情视频| 国产主播在线观看一区二区| 日本 欧美在线| 国产精品美女特级片免费视频播放器 | 国产激情欧美一区二区| 欧美一区二区精品小视频在线| 日韩欧美精品v在线| 免费观看人在逋| 变态另类丝袜制服| 国产成人欧美在线观看| 曰老女人黄片| 国产91精品成人一区二区三区| a在线观看视频网站| 国产精品久久电影中文字幕| 国产黄片美女视频| 国产私拍福利视频在线观看| 国产精品av久久久久免费| 国产成人精品无人区| 亚洲五月婷婷丁香| 精品无人区乱码1区二区| 精品午夜福利视频在线观看一区| 欧美大码av| 国产99久久九九免费精品| 亚洲av成人一区二区三| 夜夜躁狠狠躁天天躁| 黑人操中国人逼视频| 一本大道久久a久久精品| 国产乱人伦免费视频| 欧美三级亚洲精品| 国产男靠女视频免费网站| 在线永久观看黄色视频| 长腿黑丝高跟| 免费看美女性在线毛片视频| 一级片免费观看大全| 一本一本综合久久| 两人在一起打扑克的视频| 日本一区二区免费在线视频| 女人被狂操c到高潮| 级片在线观看| 91老司机精品| svipshipincom国产片| 亚洲第一欧美日韩一区二区三区| 看片在线看免费视频| 每晚都被弄得嗷嗷叫到高潮| 国产午夜精品久久久久久| 在线观看午夜福利视频| 国产成人av激情在线播放| 久久精品91蜜桃| АⅤ资源中文在线天堂| 亚洲人成电影免费在线| 午夜日韩欧美国产| 欧美在线黄色| 男女视频在线观看网站免费 | 国产精品永久免费网站| 真人做人爱边吃奶动态| 在线视频色国产色| 欧美zozozo另类| 色av中文字幕| 国产男靠女视频免费网站| 亚洲精品美女久久av网站| 黄片大片在线免费观看| 精品不卡国产一区二区三区| 久久亚洲真实| 亚洲欧美日韩东京热| 90打野战视频偷拍视频| 亚洲av中文字字幕乱码综合| 黄色视频不卡| 国产成人av激情在线播放| 日韩精品免费视频一区二区三区| 国产精华一区二区三区| 韩国av一区二区三区四区| 夜夜躁狠狠躁天天躁| 在线观看舔阴道视频| 久久精品国产亚洲av香蕉五月| tocl精华| 欧美日本视频| 国产精品一及| 中文字幕精品亚洲无线码一区| 免费无遮挡裸体视频| 1024香蕉在线观看| 丝袜人妻中文字幕| 精品久久久久久久久久免费视频| 久久婷婷成人综合色麻豆| 在线观看美女被高潮喷水网站 | 91大片在线观看| 欧美三级亚洲精品| av有码第一页| www国产在线视频色| 国产精品亚洲一级av第二区| 一个人免费在线观看的高清视频| 国产片内射在线| 一进一出抽搐gif免费好疼| 亚洲精品中文字幕一二三四区| 夜夜躁狠狠躁天天躁| av国产免费在线观看| 国产欧美日韩一区二区精品| 国产人伦9x9x在线观看| 曰老女人黄片| 国产av一区二区精品久久| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 亚洲精品美女久久av网站| 蜜桃久久精品国产亚洲av| 婷婷精品国产亚洲av在线| 色播亚洲综合网| 日韩大尺度精品在线看网址| 国产爱豆传媒在线观看 | 无遮挡黄片免费观看| 国产成年人精品一区二区| 99久久国产精品久久久| 动漫黄色视频在线观看| 亚洲九九香蕉| 久久久精品大字幕| 亚洲片人在线观看| 操出白浆在线播放| 久热爱精品视频在线9| 亚洲精品美女久久久久99蜜臀| 亚洲成av人片免费观看| 欧美不卡视频在线免费观看 | 国产av一区在线观看免费| 日韩中文字幕欧美一区二区| 亚洲国产欧美人成| e午夜精品久久久久久久| 精品不卡国产一区二区三区| 亚洲av电影在线进入| 天堂影院成人在线观看| www日本黄色视频网| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 国产欧美日韩一区二区三| 99久久国产精品久久久| av福利片在线观看| 99在线人妻在线中文字幕| 91麻豆av在线| 国产亚洲精品久久久久久毛片| 哪里可以看免费的av片| 亚洲国产中文字幕在线视频| 亚洲自拍偷在线| 亚洲欧美日韩高清专用| 欧美精品啪啪一区二区三区| 国产成人系列免费观看| 两性夫妻黄色片| 日韩国内少妇激情av| 国产高清视频在线观看网站| 老司机午夜十八禁免费视频| 精品福利观看| 久久久久久久午夜电影| 高清在线国产一区| 午夜免费成人在线视频| 中亚洲国语对白在线视频| 看黄色毛片网站| 婷婷六月久久综合丁香| 黄色片一级片一级黄色片| 激情在线观看视频在线高清| 久久久久久国产a免费观看| 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| 久久午夜亚洲精品久久| 精品国内亚洲2022精品成人| 在线永久观看黄色视频| 中出人妻视频一区二区| 亚洲全国av大片| av福利片在线| 亚洲真实伦在线观看| 国产精品久久视频播放| 久久香蕉精品热| 丝袜人妻中文字幕| 三级国产精品欧美在线观看 | 欧美久久黑人一区二区| 国产av不卡久久| 欧美最黄视频在线播放免费| 国产三级黄色录像| 国产亚洲精品一区二区www| 免费观看人在逋| 国产精品乱码一区二三区的特点| 国产成人精品久久二区二区91| 曰老女人黄片| 久久久久九九精品影院| 亚洲中文字幕一区二区三区有码在线看 | 男女视频在线观看网站免费 | 无人区码免费观看不卡| 欧美又色又爽又黄视频| 成人国语在线视频| 国产又色又爽无遮挡免费看| avwww免费| 色尼玛亚洲综合影院| 亚洲激情在线av| 日本一本二区三区精品| av有码第一页| 又黄又粗又硬又大视频| 亚洲aⅴ乱码一区二区在线播放 | 可以在线观看的亚洲视频| 97碰自拍视频| 中出人妻视频一区二区| 国产精品久久久久久精品电影| 最近最新中文字幕大全电影3| 国产精品一区二区三区四区久久| 久久 成人 亚洲| 舔av片在线| 老司机深夜福利视频在线观看| 两个人视频免费观看高清| 国产乱人伦免费视频| 十八禁人妻一区二区| 动漫黄色视频在线观看| 中文字幕精品亚洲无线码一区| 国产三级在线视频| 午夜福利成人在线免费观看| 国产精品免费一区二区三区在线| 日本在线视频免费播放| www.999成人在线观看| 亚洲av中文字字幕乱码综合| 色尼玛亚洲综合影院| 搡老岳熟女国产| 免费在线观看成人毛片| 老汉色∧v一级毛片| 高清毛片免费观看视频网站| 中文资源天堂在线| 国产精品免费一区二区三区在线| 国产av一区在线观看免费| 男女下面进入的视频免费午夜| 色av中文字幕| 在线观看一区二区三区| 国产97色在线日韩免费| 国产精品一区二区精品视频观看| 制服丝袜大香蕉在线| 日日摸夜夜添夜夜添小说| 免费av毛片视频| 可以免费在线观看a视频的电影网站| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美精品v在线| 国产麻豆成人av免费视频| 全区人妻精品视频| 亚洲国产欧美一区二区综合| 日韩国内少妇激情av| 欧美性猛交黑人性爽| 亚洲av成人精品一区久久| 婷婷六月久久综合丁香| 老司机午夜福利在线观看视频| 免费看a级黄色片| 免费在线观看成人毛片| 国产精品香港三级国产av潘金莲| 国产aⅴ精品一区二区三区波| 黄色视频不卡| 99国产精品一区二区三区| 夜夜爽天天搞| 一区二区三区高清视频在线| 两人在一起打扑克的视频| 欧美av亚洲av综合av国产av| 亚洲国产中文字幕在线视频| 久久精品影院6| 丁香六月欧美| 黄色视频不卡| 淫妇啪啪啪对白视频| 两个人免费观看高清视频| 免费在线观看亚洲国产| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 国产精品1区2区在线观看.| 久久久久九九精品影院| 久久久久久人人人人人| 女人被狂操c到高潮| 国产精品久久久久久精品电影| 成人一区二区视频在线观看| 18禁美女被吸乳视频| 国产视频一区二区在线看| 一进一出抽搐动态| 男女下面进入的视频免费午夜| 91老司机精品| 少妇裸体淫交视频免费看高清 | 久久精品综合一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 这个男人来自地球电影免费观看| 午夜福利高清视频| 午夜福利成人在线免费观看| 亚洲精品美女久久久久99蜜臀| 好看av亚洲va欧美ⅴa在| 色精品久久人妻99蜜桃| 亚洲全国av大片| 人妻久久中文字幕网| 亚洲,欧美精品.| 9191精品国产免费久久| 国模一区二区三区四区视频 | 精品久久久久久久毛片微露脸| 激情在线观看视频在线高清| 欧美人与性动交α欧美精品济南到| 亚洲真实伦在线观看| 国内久久婷婷六月综合欲色啪| 黄色视频不卡| 婷婷亚洲欧美| 亚洲人成电影免费在线| 男人舔女人下体高潮全视频| 看免费av毛片| 成人午夜高清在线视频| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| 久久天堂一区二区三区四区| 黑人操中国人逼视频| 免费看十八禁软件| 午夜免费成人在线视频| 天堂影院成人在线观看| 亚洲性夜色夜夜综合| 精品不卡国产一区二区三区| 午夜精品久久久久久毛片777| 在线视频色国产色| 久久久精品大字幕| 51午夜福利影视在线观看| 全区人妻精品视频| 亚洲免费av在线视频| 18禁国产床啪视频网站| 好看av亚洲va欧美ⅴa在| 18禁裸乳无遮挡免费网站照片| 午夜福利在线在线| 久久午夜综合久久蜜桃| 久久中文字幕一级| 精品国产超薄肉色丝袜足j| 亚洲一区高清亚洲精品| 99精品欧美一区二区三区四区| 一本大道久久a久久精品| 毛片女人毛片| 无遮挡黄片免费观看| 久久久久久久久免费视频了| 亚洲av第一区精品v没综合| 999精品在线视频| 欧美一级毛片孕妇| 19禁男女啪啪无遮挡网站| 婷婷亚洲欧美| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月| 人人妻人人看人人澡| 一边摸一边抽搐一进一小说| 精品国产亚洲在线| 国产熟女午夜一区二区三区| 日韩欧美国产一区二区入口| 国产单亲对白刺激| 久久婷婷人人爽人人干人人爱| 最新美女视频免费是黄的| 变态另类丝袜制服| 欧美日韩精品网址| 很黄的视频免费| 国产蜜桃级精品一区二区三区| 亚洲av成人一区二区三| 黄片大片在线免费观看| 亚洲国产精品sss在线观看| 日韩大尺度精品在线看网址| 亚洲aⅴ乱码一区二区在线播放 | 人妻久久中文字幕网| 麻豆国产97在线/欧美 | 久久精品国产亚洲av高清一级| 12—13女人毛片做爰片一| 国产激情久久老熟女| 91成年电影在线观看| 亚洲专区国产一区二区| 亚洲人成伊人成综合网2020| 99热只有精品国产| 特级一级黄色大片| 午夜精品久久久久久毛片777| 久久精品夜夜夜夜夜久久蜜豆 | 深夜精品福利| 欧美性长视频在线观看| 国产成人aa在线观看| 国产精品久久久久久精品电影| 国产亚洲精品一区二区www| 一本一本综合久久| 久久精品夜夜夜夜夜久久蜜豆 | 99在线人妻在线中文字幕| 午夜激情av网站| 国产视频内射| 免费高清视频大片| 亚洲精品一卡2卡三卡4卡5卡| e午夜精品久久久久久久| 国产高清视频在线观看网站| 日韩国内少妇激情av| 免费在线观看黄色视频的| 国产精品av视频在线免费观看| 国产熟女xx| 久9热在线精品视频| 男女做爰动态图高潮gif福利片| 搡老熟女国产l中国老女人| 好男人在线观看高清免费视频| 欧美黑人精品巨大| 亚洲午夜理论影院| 午夜福利视频1000在线观看| 亚洲国产看品久久| 黑人操中国人逼视频| 后天国语完整版免费观看| 亚洲美女视频黄频| 国产伦一二天堂av在线观看| 舔av片在线| 老汉色∧v一级毛片| 欧美国产日韩亚洲一区| 免费看美女性在线毛片视频| 国产高清视频在线播放一区| 亚洲欧洲精品一区二区精品久久久| 精品免费久久久久久久清纯| 欧美+亚洲+日韩+国产|