• <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
    夾在兩陣風之間
    欧美大码av| 人人妻人人澡欧美一区二区| √禁漫天堂资源中文www| 两个人视频免费观看高清| 久久亚洲精品不卡| av国产免费在线观看| 黑人巨大精品欧美一区二区mp4| 最近最新中文字幕大全电影3| 午夜福利在线观看吧| 欧美性猛交黑人性爽| 亚洲国产精品sss在线观看| 亚洲人成网站在线播放欧美日韩| 成人永久免费在线观看视频| 中文在线观看免费www的网站 | 人妻久久中文字幕网| 色播亚洲综合网| 亚洲人成77777在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲美女视频黄频| 国产三级中文精品| av福利片在线观看| 丁香六月欧美| 在线播放国产精品三级| xxxwww97欧美| 成人手机av| 一个人观看的视频www高清免费观看 | 天天添夜夜摸| 日韩 欧美 亚洲 中文字幕| 国产视频一区二区在线看| 国产激情偷乱视频一区二区| 国产精品 欧美亚洲| 两人在一起打扑克的视频| 十八禁人妻一区二区| 黄色视频,在线免费观看| 十八禁网站免费在线| 啦啦啦免费观看视频1| av中文乱码字幕在线| 国产亚洲精品久久久久久毛片| 亚洲av成人一区二区三| 久久久久精品国产欧美久久久| 午夜老司机福利片| 少妇熟女aⅴ在线视频| 成年女人毛片免费观看观看9| 麻豆av在线久日| 叶爱在线成人免费视频播放| 欧美黑人欧美精品刺激| 日韩成人在线观看一区二区三区| 国产精品久久久久久久电影 | 天堂动漫精品| 久久精品国产亚洲av高清一级| 久久九九热精品免费| 成人国语在线视频| 美女扒开内裤让男人捅视频| 日本三级黄在线观看| 久久久久久久久免费视频了| 精品高清国产在线一区| 久久 成人 亚洲| 香蕉久久夜色| 一本大道久久a久久精品| 一边摸一边做爽爽视频免费| 妹子高潮喷水视频| 男人舔女人的私密视频| 亚洲精品中文字幕在线视频| 女同久久另类99精品国产91| 日本熟妇午夜| 国产激情欧美一区二区| 国内久久婷婷六月综合欲色啪| www国产在线视频色| 中文在线观看免费www的网站 | 午夜影院日韩av| 亚洲欧美日韩无卡精品| 久久人人精品亚洲av| 男女下面进入的视频免费午夜| 国内揄拍国产精品人妻在线| 欧美黑人欧美精品刺激| 亚洲av成人av| 国产午夜福利久久久久久| 免费观看精品视频网站| 最近最新中文字幕大全电影3| 精品高清国产在线一区| 国产精品电影一区二区三区| 日韩国内少妇激情av| 亚洲专区字幕在线| 不卡av一区二区三区| 久久久久国产精品人妻aⅴ院| 超碰成人久久| 老司机深夜福利视频在线观看| 亚洲av五月六月丁香网| 两性午夜刺激爽爽歪歪视频在线观看 | 成人18禁高潮啪啪吃奶动态图| 久久久久久久精品吃奶| 天堂√8在线中文| 亚洲男人天堂网一区| 在线永久观看黄色视频| 每晚都被弄得嗷嗷叫到高潮| 欧美成人免费av一区二区三区| av福利片在线| 亚洲人成电影免费在线| 日本撒尿小便嘘嘘汇集6| 成人三级做爰电影| 国产在线精品亚洲第一网站| 丁香六月欧美| 国产亚洲精品一区二区www| 成人精品一区二区免费| 怎么达到女性高潮| 亚洲av成人一区二区三| 国产精品乱码一区二三区的特点| 妹子高潮喷水视频| 亚洲乱码一区二区免费版| 身体一侧抽搐| 亚洲一区中文字幕在线| 国产伦一二天堂av在线观看| 熟妇人妻久久中文字幕3abv| 黑人操中国人逼视频| 免费看日本二区| 黄色片一级片一级黄色片| a在线观看视频网站| 亚洲熟妇中文字幕五十中出| 日韩 欧美 亚洲 中文字幕| 色哟哟哟哟哟哟| 老司机在亚洲福利影院| 9191精品国产免费久久| www.精华液| 女警被强在线播放| 精品久久久久久久末码| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩无卡精品| 色哟哟哟哟哟哟| 757午夜福利合集在线观看| 精品第一国产精品| 最近在线观看免费完整版| 欧美性猛交╳xxx乱大交人| 中文字幕人成人乱码亚洲影| 久久香蕉精品热| 亚洲色图 男人天堂 中文字幕| 国产精品日韩av在线免费观看| 亚洲av成人不卡在线观看播放网| 妹子高潮喷水视频| 亚洲国产精品合色在线| 淫妇啪啪啪对白视频| 欧美成人免费av一区二区三区| 18禁裸乳无遮挡免费网站照片| 禁无遮挡网站| 欧美日韩乱码在线| 中文字幕人妻丝袜一区二区| 亚洲男人天堂网一区| 国产一区二区在线av高清观看| 亚洲性夜色夜夜综合| 亚洲国产欧洲综合997久久,| АⅤ资源中文在线天堂| 91av网站免费观看| 成人手机av| 非洲黑人性xxxx精品又粗又长| 日韩欧美国产在线观看| 精品第一国产精品| 中出人妻视频一区二区| 日本一二三区视频观看| 国产午夜精品久久久久久| 国产一区二区在线观看日韩 | 亚洲全国av大片| 国产高清视频在线播放一区| 高潮久久久久久久久久久不卡| 国产区一区二久久| 1024视频免费在线观看| 变态另类丝袜制服| av国产免费在线观看| 亚洲精品中文字幕在线视频| 欧美日本视频| 高潮久久久久久久久久久不卡| 久久国产精品人妻蜜桃| 男女那种视频在线观看| bbb黄色大片| 欧美人与性动交α欧美精品济南到| 成人手机av| 黄频高清免费视频| 亚洲国产欧美人成| 两个人免费观看高清视频| 精品少妇一区二区三区视频日本电影| 国产午夜精品论理片| 最新美女视频免费是黄的| 怎么达到女性高潮| 老司机在亚洲福利影院| av中文乱码字幕在线| av天堂在线播放| 白带黄色成豆腐渣| 午夜日韩欧美国产| 正在播放国产对白刺激| 精品欧美国产一区二区三| 又大又爽又粗| 岛国在线免费视频观看| 国内毛片毛片毛片毛片毛片| a级毛片a级免费在线| 亚洲男人天堂网一区| 此物有八面人人有两片| 99久久99久久久精品蜜桃| 精品日产1卡2卡| 嫁个100分男人电影在线观看| 国模一区二区三区四区视频 | 国产真人三级小视频在线观看| 色哟哟哟哟哟哟| 国产精品一区二区精品视频观看| 国产亚洲精品综合一区在线观看 | 亚洲成人免费电影在线观看| 国产乱人伦免费视频| 成人一区二区视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲精品国产精品久久久不卡| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩无卡精品| 日韩中文字幕欧美一区二区| 久久精品亚洲精品国产色婷小说| 99久久国产精品久久久| 日韩高清综合在线| 最近视频中文字幕2019在线8| 亚洲精品一区av在线观看| 亚洲国产精品合色在线| 国产熟女xx| 99久久99久久久精品蜜桃| 午夜激情av网站| 变态另类成人亚洲欧美熟女| 亚洲精品一卡2卡三卡4卡5卡| 成人高潮视频无遮挡免费网站| 在线观看舔阴道视频| 人人妻人人澡欧美一区二区| 全区人妻精品视频| 午夜视频精品福利| 村上凉子中文字幕在线| 亚洲国产精品sss在线观看| 国产一区二区在线观看日韩 | 在线观看www视频免费| 久热爱精品视频在线9| 麻豆成人av在线观看| 国产午夜福利久久久久久| 身体一侧抽搐| ponron亚洲| 此物有八面人人有两片| 国产伦在线观看视频一区| 色噜噜av男人的天堂激情| 国产精品久久久人人做人人爽| 久久这里只有精品中国| 亚洲精品粉嫩美女一区| 老鸭窝网址在线观看| 老鸭窝网址在线观看| 99国产精品一区二区蜜桃av| 色播亚洲综合网| a在线观看视频网站| 操出白浆在线播放| 非洲黑人性xxxx精品又粗又长| 别揉我奶头~嗯~啊~动态视频| 国模一区二区三区四区视频 | 天天躁狠狠躁夜夜躁狠狠躁| 国产成人系列免费观看| 男女下面进入的视频免费午夜| 久久精品国产清高在天天线| 亚洲熟妇中文字幕五十中出| 久久香蕉国产精品| 久久久久久国产a免费观看| 大型av网站在线播放| 午夜免费激情av| 免费搜索国产男女视频| 搡老妇女老女人老熟妇| www.www免费av| 成年人黄色毛片网站| 日韩欧美国产在线观看| 日本成人三级电影网站| 亚洲国产日韩欧美精品在线观看 | 色在线成人网| 色噜噜av男人的天堂激情| 国产欧美日韩一区二区三| 琪琪午夜伦伦电影理论片6080| 最近最新中文字幕大全免费视频| 久久精品国产亚洲av高清一级| 一本久久中文字幕| 免费在线观看成人毛片| 亚洲男人的天堂狠狠| 日日爽夜夜爽网站| 久久久久国产精品人妻aⅴ院| 欧美+亚洲+日韩+国产| 99国产极品粉嫩在线观看| 久久天堂一区二区三区四区| 亚洲18禁久久av| 国产黄a三级三级三级人| 老鸭窝网址在线观看| 亚洲五月天丁香| 国产成人aa在线观看| 性欧美人与动物交配| 亚洲国产中文字幕在线视频| 黄色a级毛片大全视频| 国产亚洲精品av在线| 老司机深夜福利视频在线观看| 国产精品乱码一区二三区的特点| 欧美日韩黄片免| 久久久久久人人人人人| av中文乱码字幕在线| 日韩欧美一区二区三区在线观看| 国产黄a三级三级三级人| 成人午夜高清在线视频| a在线观看视频网站| x7x7x7水蜜桃| 级片在线观看| 老司机福利观看| 久久久久久久久免费视频了| 免费在线观看日本一区| 久久久久免费精品人妻一区二区| 国产精品亚洲一级av第二区| 午夜福利在线在线| 一进一出好大好爽视频| 男女之事视频高清在线观看| 亚洲成av人片免费观看| aaaaa片日本免费| 男女那种视频在线观看| 国产精品一区二区三区四区免费观看 | 国产v大片淫在线免费观看| 久久人妻av系列| 又黄又爽又免费观看的视频| 2021天堂中文幕一二区在线观| 十八禁网站免费在线| 丝袜人妻中文字幕| aaaaa片日本免费| 欧美另类亚洲清纯唯美| 一级毛片高清免费大全| 国产精品一区二区精品视频观看| 午夜福利18| 欧美一区二区国产精品久久精品 | 精品久久久久久久末码| 非洲黑人性xxxx精品又粗又长| 免费在线观看日本一区| 此物有八面人人有两片| 免费在线观看完整版高清| 午夜亚洲福利在线播放| 午夜激情福利司机影院| 国产成人精品久久二区二区免费| 久久久久久人人人人人| 88av欧美| 色精品久久人妻99蜜桃| 亚洲一区高清亚洲精品| 伦理电影免费视频| 9191精品国产免费久久| 国产一级毛片七仙女欲春2| 免费看a级黄色片| 色综合站精品国产| 91国产中文字幕| 日韩欧美在线二视频| 久久久久免费精品人妻一区二区| 女同久久另类99精品国产91| 免费看a级黄色片| 老汉色av国产亚洲站长工具| 国产午夜精品久久久久久| 黄色视频不卡| 国产午夜精品久久久久久| 九九热线精品视视频播放| 久久久久久久午夜电影| 免费电影在线观看免费观看| 国产成人系列免费观看| 97碰自拍视频| 亚洲电影在线观看av| 久热爱精品视频在线9| 国产黄片美女视频| 18禁黄网站禁片免费观看直播| 亚洲精品色激情综合| 一夜夜www| 成人三级做爰电影| 天天一区二区日本电影三级| 国产精品久久久久久久电影 | 黄色 视频免费看| 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区三| 天堂√8在线中文| 老鸭窝网址在线观看| 麻豆成人午夜福利视频| 91麻豆av在线| 久久久久免费精品人妻一区二区| 一区二区三区高清视频在线| 十八禁网站免费在线| 嫁个100分男人电影在线观看| www.精华液| www日本在线高清视频| 国产av在哪里看| 视频区欧美日本亚洲| 很黄的视频免费| 看黄色毛片网站| 人人妻,人人澡人人爽秒播| 麻豆国产97在线/欧美 | 国产麻豆成人av免费视频| 免费观看精品视频网站| 精品国产超薄肉色丝袜足j| 日本精品一区二区三区蜜桃| 黄片大片在线免费观看| 久久久精品大字幕| 精品国产超薄肉色丝袜足j| 中文字幕人妻丝袜一区二区| 国产av又大| 男女那种视频在线观看| 欧美三级亚洲精品| 波多野结衣高清作品| 在线观看免费日韩欧美大片| 99热这里只有是精品50| 亚洲中文字幕日韩| 久久天堂一区二区三区四区| 1024手机看黄色片| 丁香六月欧美| 亚洲专区中文字幕在线| 国产亚洲欧美在线一区二区| 99热6这里只有精品| 男男h啪啪无遮挡| 亚洲成人中文字幕在线播放| 欧美zozozo另类| 日本在线视频免费播放| 亚洲国产高清在线一区二区三| 制服丝袜大香蕉在线| 日日爽夜夜爽网站| 国产又色又爽无遮挡免费看| 很黄的视频免费| 精品少妇一区二区三区视频日本电影| 久久精品综合一区二区三区| 成年女人毛片免费观看观看9| 小说图片视频综合网站| 深夜精品福利| 一进一出抽搐动态| 国产亚洲精品综合一区在线观看 | 少妇裸体淫交视频免费看高清 | 精品久久久久久久久久免费视频| 熟女电影av网| www.999成人在线观看| 可以在线观看的亚洲视频| 欧美日本视频| 最好的美女福利视频网| 亚洲美女黄片视频| 久久亚洲真实| 国产99久久九九免费精品| av天堂在线播放| 制服丝袜大香蕉在线| 久久久久国产一级毛片高清牌| 国产精品98久久久久久宅男小说| 在线免费观看的www视频| 精品欧美国产一区二区三| 久久 成人 亚洲| 久久久久亚洲av毛片大全| 淫秽高清视频在线观看| 久久久久性生活片| 白带黄色成豆腐渣| 天堂动漫精品| 国产精品美女特级片免费视频播放器 | 日韩 欧美 亚洲 中文字幕| 一个人观看的视频www高清免费观看 | 熟女电影av网| 国产高清视频在线观看网站| 18禁黄网站禁片午夜丰满| 亚洲午夜理论影院| 小说图片视频综合网站| 国产一区二区三区视频了| 色综合婷婷激情| 午夜视频精品福利| 亚洲av熟女| 男女那种视频在线观看| 久久欧美精品欧美久久欧美| 国产亚洲精品久久久久5区| 正在播放国产对白刺激| 老熟妇仑乱视频hdxx| 人成视频在线观看免费观看| 男女那种视频在线观看| 午夜福利在线在线| 啪啪无遮挡十八禁网站| xxx96com| 中文字幕人成人乱码亚洲影| 十八禁人妻一区二区| 亚洲国产日韩欧美精品在线观看 | 中亚洲国语对白在线视频| 天堂影院成人在线观看| 老司机靠b影院| 国产一区二区三区在线臀色熟女| √禁漫天堂资源中文www| 国产精品综合久久久久久久免费| 免费观看人在逋| 国产高清videossex| 夜夜看夜夜爽夜夜摸| av在线天堂中文字幕| 国产av一区在线观看免费| 2021天堂中文幕一二区在线观| 久久精品91无色码中文字幕| 欧美绝顶高潮抽搐喷水| 亚洲狠狠婷婷综合久久图片| 欧美大码av| 亚洲 国产 在线| 欧美zozozo另类| 国内精品久久久久精免费| 亚洲自拍偷在线| 亚洲精品久久国产高清桃花| 日日夜夜操网爽| 美女大奶头视频| 桃红色精品国产亚洲av| 香蕉久久夜色| 国产精品,欧美在线| 亚洲欧美激情综合另类| 国产亚洲精品久久久久久毛片| 欧美性猛交╳xxx乱大交人| 久久久久精品国产欧美久久久| 悠悠久久av| 国产高清有码在线观看视频 | 午夜日韩欧美国产| 午夜免费观看网址| 国产成人欧美在线观看| 夜夜看夜夜爽夜夜摸| 久久亚洲真实| 天天一区二区日本电影三级| 黄色视频,在线免费观看| 成人精品一区二区免费| 亚洲无线在线观看| 日日夜夜操网爽| 免费高清视频大片| 亚洲av成人精品一区久久| 欧美午夜高清在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品国产精品久久久不卡| 人妻丰满熟妇av一区二区三区| 精品国产美女av久久久久小说| 法律面前人人平等表现在哪些方面| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 老司机靠b影院| 99精品久久久久人妻精品| 老汉色∧v一级毛片| 国产成人影院久久av| 在线十欧美十亚洲十日本专区| 亚洲性夜色夜夜综合| 日韩国内少妇激情av| 性色av乱码一区二区三区2| 超碰成人久久| 亚洲一码二码三码区别大吗| 曰老女人黄片| 午夜福利视频1000在线观看| av福利片在线观看| 午夜免费激情av| 日本黄大片高清| 久久香蕉激情| 午夜视频精品福利| 国产成人av教育| 村上凉子中文字幕在线| 香蕉国产在线看| 不卡一级毛片| 麻豆成人午夜福利视频| 悠悠久久av| 欧美色视频一区免费| 18禁观看日本| 女警被强在线播放| 精品电影一区二区在线| 午夜精品久久久久久毛片777| 999精品在线视频| 精品熟女少妇八av免费久了| 亚洲 国产 在线| 亚洲成a人片在线一区二区| 亚洲成人精品中文字幕电影| 人妻久久中文字幕网| 亚洲美女视频黄频| 亚洲欧洲精品一区二区精品久久久| 国产精品免费一区二区三区在线| 国内精品久久久久精免费| 久久久久久国产a免费观看| 欧美高清成人免费视频www| 岛国在线免费视频观看| 亚洲全国av大片| 成人18禁在线播放| 久久中文字幕人妻熟女| 亚洲中文字幕一区二区三区有码在线看 | 国产av又大| 精华霜和精华液先用哪个| 99精品在免费线老司机午夜| 99久久无色码亚洲精品果冻| 国产三级中文精品| 国产精品香港三级国产av潘金莲| 日日爽夜夜爽网站| 中亚洲国语对白在线视频| av中文乱码字幕在线| 精品久久久久久久久久免费视频| 日本a在线网址| 两个人的视频大全免费| 法律面前人人平等表现在哪些方面| 视频区欧美日本亚洲| 欧美成人午夜精品| 成年女人毛片免费观看观看9| 亚洲熟妇熟女久久| 婷婷精品国产亚洲av| 亚洲成a人片在线一区二区| 亚洲18禁久久av| 亚洲av五月六月丁香网| 亚洲人成伊人成综合网2020| 88av欧美| 老熟妇乱子伦视频在线观看| 老司机午夜十八禁免费视频| 日本精品一区二区三区蜜桃| 白带黄色成豆腐渣| xxxwww97欧美| АⅤ资源中文在线天堂| 黄色成人免费大全| 性欧美人与动物交配| 九色成人免费人妻av| 久久国产精品影院| 欧美黄色淫秽网站| 亚洲精品国产精品久久久不卡| 麻豆成人av在线观看| 可以免费在线观看a视频的电影网站| 国产午夜精品久久久久久| 黑人巨大精品欧美一区二区mp4| 国产单亲对白刺激| 久久久久久久午夜电影| 亚洲成人免费电影在线观看| 色综合欧美亚洲国产小说| 亚洲av电影在线进入| 不卡一级毛片| 成人18禁在线播放| 一个人免费在线观看电影 | 黄色视频不卡| 亚洲人成网站在线播放欧美日韩| 色综合欧美亚洲国产小说| 中文字幕熟女人妻在线|