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

    基于風(fēng)壓譜和Hermite模型的大跨干煤棚風(fēng)壓場(chǎng)數(shù)值模擬研究

    2018-12-21 12:09:56黃銘楓孫軒濤樓文娟胡德軍
    振動(dòng)與沖擊 2018年23期
    關(guān)鍵詞:角下峰度風(fēng)洞試驗(yàn)

    黃銘楓, 孫軒濤, 馮 鶴, 樓文娟, 胡德軍

    (1. 浙江大學(xué) 結(jié)構(gòu)工程研究所, 杭州 310058; 2. 浙江耀廈控股集團(tuán)有限公司, 杭州 310011)

    風(fēng)荷載是控制結(jié)構(gòu)安全的重要荷載,特別是對(duì)于大跨和高層結(jié)構(gòu),風(fēng)荷載及其動(dòng)力效應(yīng)常常成為設(shè)計(jì)的主導(dǎo)因素。研究表明在屋面結(jié)構(gòu)的來流前緣、后部尾流區(qū)及角區(qū)附近,風(fēng)荷載常常表現(xiàn)出強(qiáng)烈的脈動(dòng)和非高斯特性[1]。與高斯過程相比,非高斯隨機(jī)過程的峰值更大,往往是導(dǎo)致屋面破壞的主要原因,因此對(duì)非高斯風(fēng)荷載的研究具有重要意義。對(duì)于干煤棚等大跨建筑,結(jié)構(gòu)的幾何非線性程度較高,相比于頻域風(fēng)振分析法,時(shí)域法可以考慮結(jié)構(gòu)的非線性因素,能更為精確的獲得結(jié)構(gòu)風(fēng)致動(dòng)力效應(yīng)。然而,大跨屋蓋結(jié)構(gòu)風(fēng)洞測(cè)壓試驗(yàn)得到的風(fēng)壓時(shí)程數(shù)據(jù)往往有限,難以滿足結(jié)構(gòu)動(dòng)力可靠度研究[2]的需要。因此,有必要研究非高斯風(fēng)壓場(chǎng)和風(fēng)荷載的數(shù)值模擬方法[3]。

    屋面脈動(dòng)風(fēng)壓場(chǎng)的特性可由時(shí)域統(tǒng)計(jì)量,如均值、方差、偏度、峰度等,和頻域的脈動(dòng)風(fēng)壓譜來描述。屋面脈動(dòng)風(fēng)壓譜的特性十分復(fù)雜,不僅受到建筑尺寸、屋蓋形狀、來流方向及湍流度影響,還與測(cè)壓系統(tǒng)采樣頻率密切相關(guān)。傳統(tǒng)的準(zhǔn)定常假定認(rèn)為,建筑表面風(fēng)壓的脈動(dòng)與來流風(fēng)速脈動(dòng)一致,因此采用來流風(fēng)譜轉(zhuǎn)換得到建筑表面風(fēng)壓譜。而實(shí)際上,大跨屋蓋表面的脈動(dòng)風(fēng)壓主要受特征湍流的影響[4],作用于屋蓋的風(fēng)場(chǎng)與自然來流特性存在較大差別,不適用準(zhǔn)定常假定[5],因而有必要引入風(fēng)壓譜模型。Kumar[6]對(duì)屋蓋表面的風(fēng)壓譜進(jìn)行研究,得到了風(fēng)壓譜的指數(shù)表達(dá)形式。孫瑛[7]分析了大跨平屋蓋和鞍形屋蓋在不同風(fēng)向角下的風(fēng)壓譜特征,根據(jù)風(fēng)壓譜形狀和峰值對(duì)應(yīng)頻率對(duì)屋蓋表面進(jìn)行了分區(qū),并給出了不同區(qū)域的三參數(shù)風(fēng)壓譜經(jīng)驗(yàn)表達(dá)式。Su等[8]基于115種典型大跨屋蓋的風(fēng)洞試驗(yàn)數(shù)據(jù),對(duì)三參數(shù)風(fēng)壓譜模型進(jìn)行了改進(jìn),并通過風(fēng)致響應(yīng)分析驗(yàn)證了該模型在工程實(shí)踐中的適用性。之后,蘇寧等[9]基于不同矢跨比、高跨比的球面屋蓋模型風(fēng)洞試驗(yàn)數(shù)據(jù),利用廣義回歸神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)了屋蓋表面風(fēng)荷載的特征參數(shù)。

    靜態(tài)轉(zhuǎn)換法(Static Transformation Methods)是非高斯風(fēng)壓模擬的一種常用方法,即將高斯隨機(jī)過程通過一定變換關(guān)系轉(zhuǎn)化成滿足指定功率譜和高階統(tǒng)計(jì)量的非高斯隨機(jī)過程。李錦華等基于Johnson轉(zhuǎn)換和數(shù)字濾波理論,提出了一種能快速生成指定偏度和峰度的非高斯脈動(dòng)風(fēng)壓的方法。李璟等[10]采用三次多項(xiàng)式將非高斯風(fēng)壓場(chǎng)轉(zhuǎn)換為一個(gè)虛擬的多變量高斯過程,并利用諧波合成法生成風(fēng)壓時(shí)程樣本,該方法能快速模擬非高斯風(fēng)壓場(chǎng),但對(duì)功率譜的模擬效果較差。 Winterstein[11]首次提出利用Hermite多項(xiàng)式將一個(gè)高斯過程轉(zhuǎn)換成非高斯過程,但該方法只能用于軟響應(yīng)過程,且對(duì)功率譜的模擬不太理想。 Gruley等[12]基于Hermit多項(xiàng)式的譜修正法提出了平穩(wěn)非高斯風(fēng)壓時(shí)程的模擬方法,但該方法模擬得到的單樣本的偏度和峰度精度較差,需要多次模擬取平均才能得到較好的結(jié)果。之后,Masters等[13]改進(jìn)了Hermit多項(xiàng)式變換函數(shù),使之在模擬非高斯時(shí)程時(shí)能更快的收斂到目標(biāo)值。徐飛等[14]利用實(shí)測(cè)時(shí)程數(shù)據(jù)的分布函數(shù)和Hermit多項(xiàng)式來構(gòu)造幅值調(diào)制信號(hào),并通過與高斯過程相乘調(diào)制來得到非高斯時(shí)程。此方法不僅可以模擬非平穩(wěn)過程,而且對(duì)于峰度較大的非高斯時(shí)程也不會(huì)發(fā)生失真現(xiàn)象,具有較高的精度。Huang等[15]等也基于Hermit多項(xiàng)式變換,利用TPP方法和Deodatis的譜修正方法[16]進(jìn)行了非高斯風(fēng)壓時(shí)程的模擬。羅俊杰等[17]采用非線性無記憶轉(zhuǎn)換法模擬多點(diǎn)非高斯過程,并提出了分解譜密度函數(shù)修正法,以保證功率譜密度函數(shù)矩陣的非負(fù)定性,但該方法只適用于隨機(jī)風(fēng)壓場(chǎng)樣本服從對(duì)數(shù)正態(tài)分布和韋布爾分布的情況。

    上述已有方法雖然都能對(duì)非高斯時(shí)程進(jìn)行模擬,但大多只能模擬單點(diǎn)非高斯時(shí)程,或是只能模擬滿足特定條件的非高斯時(shí)程。本文作者對(duì)Huang等的單點(diǎn)非高斯風(fēng)壓模擬方法進(jìn)行了改進(jìn),包括調(diào)整了原有譜修正方法中的相關(guān)參數(shù),使算法能更快速收斂;另外利用蘇寧等提出的四參數(shù)風(fēng)壓譜模型來構(gòu)造大跨煤棚表面多點(diǎn)風(fēng)壓譜矩陣,提出了一種新的大跨煤棚網(wǎng)架非高斯風(fēng)壓場(chǎng)高效模擬方法。本文方法中只需要少數(shù)測(cè)點(diǎn)數(shù)據(jù)即可歸納得出煤棚表面風(fēng)壓經(jīng)驗(yàn)譜,從而對(duì)整個(gè)大跨煤棚屋面非高斯風(fēng)壓場(chǎng)進(jìn)行高效模擬,風(fēng)壓譜及非高斯特性的模擬精度均可以通過相關(guān)評(píng)價(jià)指標(biāo)進(jìn)行有效控制。以某大跨開敞式干煤棚網(wǎng)架風(fēng)洞試驗(yàn)數(shù)據(jù)為基準(zhǔn),對(duì)多個(gè)代表性測(cè)點(diǎn)脈動(dòng)風(fēng)壓的模擬結(jié)果進(jìn)行了分析,驗(yàn)證了本文非高斯風(fēng)壓場(chǎng)模擬方法的精度和有效性。

    1 非高斯風(fēng)壓場(chǎng)模擬理論

    1.1 四參數(shù)風(fēng)壓譜模型

    各測(cè)點(diǎn)風(fēng)壓系數(shù)按下式計(jì)算

    (1)

    式中:Cpi是測(cè)點(diǎn)i的凈風(fēng)壓系數(shù);Pi為測(cè)點(diǎn)i的風(fēng)壓值;P∞為參考點(diǎn)靜壓力值;V∞為參考點(diǎn)的風(fēng)速。

    風(fēng)壓時(shí)程的平均值和均方根分別按式(2)、式(3)計(jì)算

    (2)

    (3)

    式中:Cpk為測(cè)點(diǎn)風(fēng)壓系數(shù)時(shí)程第k步的值,K為測(cè)點(diǎn)風(fēng)壓時(shí)程數(shù)據(jù)總數(shù)。

    脈動(dòng)風(fēng)壓系數(shù)自功率譜(PSD)采用四參數(shù)模型表示

    (4)

    (5)

    式中:S為無量綱脈動(dòng)風(fēng)壓系數(shù)自功率譜密度,f為頻率(Hz),Scp(f)為脈動(dòng)風(fēng)壓系數(shù)自功率譜密度,F(xiàn)為基于煤棚平均高度h以及該高度處來流平均風(fēng)速Uh的折減頻率,F(xiàn)′為無量綱頻率,Sm、Fm、k1、k2為風(fēng)壓經(jīng)驗(yàn)譜的四個(gè)參數(shù),分別表示無量綱自功率譜的峰值、峰值頻率、低頻段上升斜率和高頻段衰減斜率,四個(gè)參數(shù)的取值與屋蓋形狀、尺寸、測(cè)點(diǎn)位置和風(fēng)向角有關(guān),蘇寧等認(rèn)為一般情況下k1=1。α為待定參數(shù),根據(jù)下式確定

    (6)

    式中:B(x,y)為Beta函數(shù)。

    互功率譜(CPSD)表達(dá)式為

    (7)

    式中:SCpij為i,j兩測(cè)點(diǎn)風(fēng)壓系數(shù)互功率譜,SCpi(f),SCpi(f)分別為i,j兩測(cè)點(diǎn)風(fēng)壓系數(shù)的自功率譜密度函數(shù),Δij為i,j兩點(diǎn)間距,kc為衡量脈動(dòng)風(fēng)壓場(chǎng)相干性的指數(shù)衰減系數(shù),即相干系數(shù)。邵帥[18]根據(jù)多組風(fēng)洞試驗(yàn)數(shù)據(jù),建議開敞式大曲率拱形屋蓋kc為3.5。

    1.2 Hermite變換關(guān)系

    基于Hermite矩的變換式是利用偏度和峰度來轉(zhuǎn)換高斯過程與非高斯過程的關(guān)系式。其中Hermite矩正變換關(guān)系式如下

    x=κ[u+h3(u2-1)+h4(u3-3u)];γ4>3

    (8)

    x為非高斯時(shí)程,u為高斯時(shí)程,γ3為偏度,γ4為峰度。

    Hermite矩逆變換關(guān)系式如下

    u=x-h3(x2-1)-h4(x3-3x);γ4<3

    (9)

    1.3 諧波疊加法

    諧波疊加法是通過一組三角余弦函數(shù)的疊加來模擬穩(wěn)態(tài)高斯過程的標(biāo)準(zhǔn)算法,具有完善的理論和較高的精度,合成公式如下

    (10)

    式中:fj(t)為t時(shí)刻j節(jié)點(diǎn)風(fēng)壓系數(shù)模擬值,Φml為均勻分布于[0,2π]上的獨(dú)立隨機(jī)相位角,ωu為截止圓頻率,頻率增量Δω=ωu/N,N為一個(gè)足夠大的頻率采樣點(diǎn)個(gè)數(shù)。ωml的表達(dá)式如下

    (11)

    H(ω)為互功率譜矩陣S(ω)的Cholesky分解

    S(ω)=H(ω)H*T(ω)

    (12)

    θjm(ω)為Hjm(ω)的幅角,即相位差。

    2 非高斯風(fēng)壓場(chǎng)模擬算法

    本文的非高斯風(fēng)壓場(chǎng)模擬算法流程圖,如圖1所示。詳細(xì)步驟說明如下。

    步驟1通過測(cè)壓風(fēng)洞試驗(yàn),獲取風(fēng)壓場(chǎng)數(shù)據(jù)試驗(yàn)樣本

    P(x,y,t)={p1(t),p2(t),…,pK(t)},減去各測(cè)點(diǎn)平均風(fēng)壓后得到目標(biāo)脈動(dòng)風(fēng)壓場(chǎng)數(shù)據(jù)樣本

    步驟3

    圖1 非高斯風(fēng)壓場(chǎng)模擬算法流程圖

    (13)

    式中:β為修正常數(shù),取值1.1,在保證精度的情況下使之快速收斂。利用修正功率譜模擬該點(diǎn)風(fēng)壓時(shí)程后返回至“(4)”,進(jìn)行下一步迭代計(jì)算

    (7) 重復(fù)上述“(3)”~“(6)”,直至完成所有測(cè)點(diǎn)非高斯風(fēng)壓時(shí)程模擬。

    3 某大跨開敞式干煤棚網(wǎng)架算例分析

    3.1 風(fēng)洞試驗(yàn)

    某大跨開敞式干煤棚網(wǎng)架主體結(jié)構(gòu)縱向長(zhǎng)240 m,跨度120 m,總高度42.757 m。根據(jù)“建筑結(jié)構(gòu)荷載規(guī)范:GB 50009—2012”,該煤棚處于A類地貌,地貌粗糙度指數(shù)α=0.12。百年一遇的基本風(fēng)壓為0.9 kPa,相當(dāng)于離地面10 m高度處的設(shè)計(jì)風(fēng)速U0=37.95 m/s。風(fēng)洞試驗(yàn)在浙江大學(xué)ZD-1邊界層風(fēng)洞中進(jìn)行,模型縮尺比為1∶150。模型表面布置了250對(duì)同步測(cè)量的內(nèi)、外壓測(cè)點(diǎn),沿縱向分為25個(gè)測(cè)區(qū),以編號(hào)A-Y(外表面)和AX-YX(內(nèi)表面)表示。測(cè)壓點(diǎn)布置圖及風(fēng)洞試驗(yàn)?zāi)P?,如圖2和圖3所示。風(fēng)洞試驗(yàn)中脈動(dòng)風(fēng)壓測(cè)量采樣頻率為625 Hz,取30 s采樣時(shí)長(zhǎng)。對(duì)應(yīng)原型的脈動(dòng)風(fēng)壓時(shí)程數(shù)據(jù)時(shí)間步長(zhǎng)為0.070 4 s,總時(shí)長(zhǎng)為1 320 s。

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

    圖3 風(fēng)洞試驗(yàn)?zāi)P?/p>

    圖4為干煤棚表面平均風(fēng)壓系數(shù)及脈動(dòng)風(fēng)壓系數(shù)均方根云圖,從圖4可知,兩者有一定相似之處:在迎風(fēng)側(cè),平均風(fēng)壓系數(shù)與脈動(dòng)風(fēng)壓系數(shù)均方根值都比較大。在0°風(fēng)向角下,煤棚頂部?jī)烧叨急容^小。但在背風(fēng)側(cè)的角點(diǎn)處,脈動(dòng)風(fēng)壓系數(shù)均方根值有明顯的增大。在90°風(fēng)向角下,上游迎風(fēng)側(cè)煤棚兩端脈動(dòng)風(fēng)壓系數(shù)均方根值較大。

    (a) 0°平均風(fēng)壓系數(shù)(b) 0°脈動(dòng)風(fēng)壓系數(shù)

    (c) 90°平均風(fēng)壓系數(shù)(d) 90°脈動(dòng)風(fēng)壓系數(shù)

    圖4 平均風(fēng)壓與脈動(dòng)風(fēng)壓均方根云圖

    Fig.4 Contour of mean and RMS values for fluctuating wind pressure

    3.2 風(fēng)壓譜擬合

    本文將0°風(fēng)向角下煤棚表面分為三個(gè)區(qū)域,分區(qū)圖如圖5(a)所示,z為最小屋蓋平面尺寸的10%和建筑高度的40%中的較小值。I區(qū)為迎風(fēng)側(cè)和頂部?jī)啥藚^(qū)域, II區(qū)為煤棚頂部中心區(qū)域, III區(qū)為背風(fēng)側(cè)。90°風(fēng)向角下煤棚表面分為兩個(gè)區(qū)域,分區(qū)圖如圖5(b)所示,絕大部分區(qū)域?yàn)橛L(fēng)側(cè)I區(qū),靠近背風(fēng)側(cè)小塊區(qū)域?yàn)镮II區(qū)。

    由于屋蓋表面的漩渦作用,各測(cè)點(diǎn)風(fēng)壓自功率譜與測(cè)點(diǎn)位置和風(fēng)向角有關(guān),圖6給出了各測(cè)點(diǎn)無量綱風(fēng)壓功率譜的比較。0°風(fēng)向角下,迎風(fēng)側(cè)測(cè)點(diǎn)受來流湍流影響較大,風(fēng)壓脈動(dòng)能量集中在中低頻處;在頂部的氣流分離區(qū),大尺度分離渦向小尺度渦轉(zhuǎn)變,風(fēng)荷載脈動(dòng)頻率較低,衰減較慢;而在背風(fēng)側(cè)尾部的再附區(qū),風(fēng)壓受小尺度渦主導(dǎo),峰值頻率較高。90°風(fēng)向角下,由于風(fēng)洞試驗(yàn)工況為無煤堆工況,氣流直接從內(nèi)部穿過,風(fēng)壓脈動(dòng)特性變化不大,大部分區(qū)域受來流湍流主導(dǎo), 只有背風(fēng)側(cè)小部分區(qū)域,風(fēng)壓譜主導(dǎo)頻率向高頻移動(dòng)。參考孫瑛所用分區(qū),結(jié)合風(fēng)洞試驗(yàn)所得干煤棚自身風(fēng)壓譜規(guī)律。

    (a) 0°(b) 90°

    圖5 風(fēng)壓譜分區(qū)圖

    Fig.5 Zoning map for wind pressure spectra

    對(duì)各區(qū)域所有測(cè)點(diǎn)無量綱風(fēng)壓譜進(jìn)行擬合,為定量評(píng)價(jià)擬合結(jié)果,定義擬合誤差ε為

    (14)

    各區(qū)域擬合風(fēng)壓譜參數(shù),見表1。從表1可知,0°風(fēng)向角下,k1為前文所述取1,各區(qū)域無量綱風(fēng)壓譜的峰值Sm較為接近,但峰值頻率Fm和衰減斜率k2有很大差別。II區(qū)Fm最小,即頂部風(fēng)壓能量集中在低頻段,k2也最??;III區(qū)Fm最大,k2也最大;I區(qū)介于II區(qū)與III區(qū)之間,風(fēng)壓譜峰值不突出,能量分布較為均勻。90°風(fēng)向角下,k1取大于1的值時(shí)擬合功率譜與實(shí)際譜較為接近;I區(qū)峰值Sm遠(yuǎn)大于III區(qū)峰值,與0°時(shí)各區(qū)峰值較為接近;對(duì)于峰值頻率Fm,I區(qū)小于III區(qū),與0°結(jié)果一致;衰減斜率k2則是I區(qū)大于III區(qū)。兩個(gè)風(fēng)向角下背風(fēng)面風(fēng)壓譜離散性均較大,擬合精度略差于迎風(fēng)面及頂部。

    限于篇幅,本文僅對(duì)0°和90°兩個(gè)代表性風(fēng)向角進(jìn)行分析和模擬。其他風(fēng)向角也有類似規(guī)律,區(qū)別僅在于煤棚表面風(fēng)壓經(jīng)驗(yàn)譜的分區(qū)及對(duì)應(yīng)的四個(gè)譜參數(shù)不同,對(duì)不同風(fēng)向角下煤棚表面進(jìn)行重新分區(qū)和擬合后,可將本方法用于其他風(fēng)向角。

    表1 煤棚表面各區(qū)域風(fēng)壓譜參數(shù)

    (a) 0°

    (b) 90°

    3.3 模擬結(jié)果分析

    3.3.1 測(cè)點(diǎn)脈動(dòng)風(fēng)壓譜模擬結(jié)果分析

    根據(jù)風(fēng)洞試驗(yàn)結(jié)果,在各個(gè)區(qū)域選取脈動(dòng)風(fēng)壓系數(shù)均方根值較大的代表性測(cè)點(diǎn)進(jìn)行風(fēng)壓模擬和結(jié)果分析,0°風(fēng)向角下選取L4、L8和L9,90°風(fēng)向角下選取A7、S7和W7。圖7給出了0°風(fēng)向角下三個(gè)代表性測(cè)點(diǎn)模擬所得的脈動(dòng)風(fēng)壓時(shí)程與風(fēng)洞試驗(yàn)數(shù)據(jù)的比較。圖8為各測(cè)點(diǎn)模擬所得脈動(dòng)風(fēng)壓功率譜與風(fēng)洞試驗(yàn)結(jié)果(目標(biāo)譜)的比較。限于篇幅,時(shí)程比較只給出0°風(fēng)向角下的結(jié)果。

    (a) 測(cè)點(diǎn)L4

    (b) 測(cè)點(diǎn)L8

    (c) 測(cè)點(diǎn)L9

    (b) 90°

    從圖7可知,各測(cè)點(diǎn)模擬所得的風(fēng)壓時(shí)程與風(fēng)洞試驗(yàn)數(shù)據(jù)的波動(dòng)形式較為一致。圖8所示代表性測(cè)點(diǎn)模擬所得和風(fēng)洞試驗(yàn)的風(fēng)壓譜也有較好的可比性。在0°風(fēng)向角下,模擬所得脈動(dòng)風(fēng)壓的功率譜在各個(gè)頻段與風(fēng)洞試驗(yàn)譜曲線均吻合較好,90°風(fēng)向角下,模擬所得脈動(dòng)風(fēng)壓的功率譜在高頻段略大于風(fēng)洞試驗(yàn)譜曲線,考慮到高頻區(qū)域功率譜密度小能量低,模擬結(jié)果已較為理想。

    為進(jìn)一步評(píng)價(jià)脈動(dòng)風(fēng)壓場(chǎng)的模擬結(jié)果,對(duì)不同測(cè)點(diǎn)之間的脈動(dòng)風(fēng)壓相干性進(jìn)行分析。0°風(fēng)向角下選取代表性測(cè)點(diǎn)L8,90°風(fēng)向角下選取代表性測(cè)點(diǎn)B7為考察對(duì)象,分別作出其與順風(fēng)向相鄰測(cè)點(diǎn)與橫風(fēng)向相鄰測(cè)點(diǎn)的互功率譜密度圖,如圖9和圖10所示。從圖中可以看出:0°風(fēng)向角,沿著順風(fēng)向,模擬所得測(cè)點(diǎn)L8與L7脈動(dòng)風(fēng)壓的互功率譜在低頻區(qū)域略小于風(fēng)洞試驗(yàn)譜曲線(目標(biāo)譜)(見圖9(a)),模擬所得測(cè)點(diǎn)L8與下游測(cè)點(diǎn)L9脈動(dòng)風(fēng)壓的互功率譜與風(fēng)洞試驗(yàn)譜曲線吻合程度較好(見圖9(b));沿著橫風(fēng)向,模擬所得測(cè)點(diǎn)L8與K8、M8脈動(dòng)風(fēng)壓的互功率譜在低頻區(qū)域均略小于風(fēng)洞試驗(yàn)譜曲線(見圖10(a)、圖10(b)),在高頻區(qū)域則吻合較好。90°風(fēng)向角下,模擬所得測(cè)點(diǎn)B7與順風(fēng)向測(cè)點(diǎn)A7、C7脈動(dòng)風(fēng)壓的互功率譜與風(fēng)洞試驗(yàn)譜曲線吻合程度較好(見圖9(c)、圖9(d));沿著橫風(fēng)向,所考察測(cè)點(diǎn)B7與橫風(fēng)向測(cè)點(diǎn)B6和B8的模擬所得脈動(dòng)風(fēng)壓的互功率譜在各個(gè)頻段均略小于目標(biāo)風(fēng)洞試驗(yàn)譜曲線(見圖10(c)、圖10(d))。由于本文模擬方法中風(fēng)壓互功率譜其相干系數(shù)取值均為3.5,不能區(qū)分不同風(fēng)向角下沿著不同方向上風(fēng)壓相干性的差異。這樣導(dǎo)致模擬所得各測(cè)點(diǎn)脈動(dòng)風(fēng)壓之間的互功率譜與實(shí)際風(fēng)洞試驗(yàn)的譜曲線有一定差距,值得在將來的工作中研究改進(jìn)。

    3.3.2 測(cè)點(diǎn)脈動(dòng)風(fēng)壓非高斯特性分析

    (a) 測(cè)點(diǎn)L8和L7(0°風(fēng)向角)

    (b) 測(cè)點(diǎn)L8和L9(0°風(fēng)向角)

    (c) 測(cè)點(diǎn)B7和A7(90°風(fēng)向角)

    (d) 測(cè)點(diǎn)B7和C7(90°風(fēng)向角)

    (a) 測(cè)點(diǎn)L8和K8(0°風(fēng)向角)

    (b) 測(cè)點(diǎn)L8和M8(0°風(fēng)向角)

    (c) 測(cè)點(diǎn)B7和B6(90°風(fēng)向角)

    (d) 測(cè)點(diǎn)B7和B8(90°風(fēng)向角)

    為評(píng)價(jià)脈動(dòng)風(fēng)壓場(chǎng)的非高斯特性,對(duì)測(cè)點(diǎn)脈動(dòng)風(fēng)壓的高階統(tǒng)計(jì)量進(jìn)行了分析。圖11為0°風(fēng)向角下脈動(dòng)風(fēng)壓場(chǎng)的偏度及峰度等值線云圖,其中圖11(a)和圖11(c)為根據(jù)風(fēng)洞試驗(yàn)測(cè)壓數(shù)據(jù)獲得的偏度和峰度分布圖,圖11(b)和圖11(d)為基于模擬所得風(fēng)壓場(chǎng)數(shù)據(jù)得出的偏度和峰度分布。從圖11可知,模擬所得高階統(tǒng)計(jì)量分布云圖與試驗(yàn)結(jié)果較為吻合。根據(jù)試驗(yàn)所得各測(cè)點(diǎn)脈動(dòng)風(fēng)壓高階統(tǒng)計(jì)量,在0°和90°風(fēng)向角下各選取4個(gè)具有較強(qiáng)非高斯性的測(cè)點(diǎn)(即,絕對(duì)偏度>0.5,峰度>4)。表2為選取測(cè)點(diǎn)的偏度和峰度模擬值和風(fēng)洞試驗(yàn)數(shù)據(jù)比較。表2的結(jié)果說明各代表性測(cè)點(diǎn)高階統(tǒng)計(jì)量的模擬值與試驗(yàn)值比較吻合。因?yàn)樵诒疚姆椒ㄖ胁捎镁哂薪馕瞿J降腍ermite變換關(guān)系,即式(8),得到相應(yīng)的非高斯時(shí)程數(shù)據(jù),通過反復(fù)的迭代修正,其高階統(tǒng)計(jì)量可與風(fēng)洞試驗(yàn)?zāi)繕?biāo)值非常接近,滿足給定的誤差要求。

    (a) 試驗(yàn)所得偏度分布(b) 模擬所得偏度分布

    (c) 試驗(yàn)所得峰度分布(d) 模擬所得峰度分布

    圖11 0°風(fēng)向角下脈動(dòng)風(fēng)壓場(chǎng)偏度與峰度的分布云圖比較

    為比較本文方法與已有相關(guān)算法的優(yōu)缺點(diǎn),利用李璟等的三次多項(xiàng)式轉(zhuǎn)換法對(duì)90°風(fēng)向角下干煤棚4個(gè)代表性測(cè)點(diǎn)的非高斯脈動(dòng)風(fēng)壓進(jìn)行模擬結(jié)果,如表3所示。對(duì)比表2結(jié)果可知,本文方法模擬所得非高斯風(fēng)壓場(chǎng)在偏度和峰度指標(biāo)上的精度均明顯高于三次多項(xiàng)式轉(zhuǎn)換法。

    表3李璟方法代表性測(cè)點(diǎn)脈動(dòng)風(fēng)壓高階統(tǒng)計(jì)量比較

    Tab.3ComparisonofthesimulatedandmeasuredhigherorderstatisticsforrepresentativemeasurementtapsbyLiJ’smethod

    風(fēng)向角測(cè)點(diǎn)偏度峰度風(fēng)洞試驗(yàn)?zāi)M結(jié)果誤差/%風(fēng)洞試驗(yàn)?zāi)M結(jié)果誤差/%90°D1-0.957-0.9943.874.2204.202-0.43D9-1.090-0.95012.844.6743.949-15.5E1-0.709-0.749-5.644.1694.057-2.69P5-0.735-0.6777.894.0093.7366.81

    4 結(jié) 論

    本文提出了一種基于經(jīng)驗(yàn)風(fēng)壓譜和Hermite模型的大跨干煤棚非高斯風(fēng)壓場(chǎng)模擬算法,并以某干煤棚網(wǎng)架結(jié)構(gòu)風(fēng)洞試驗(yàn)為工程實(shí)例,對(duì)非高斯風(fēng)壓場(chǎng)模擬結(jié)果進(jìn)行了分析和驗(yàn)證。主要結(jié)論如下:

    (1) 大跨干煤棚屋面經(jīng)驗(yàn)風(fēng)壓譜可根據(jù)所處位置分成若干個(gè)區(qū)域,各區(qū)域的無量綱風(fēng)壓譜峰值頻率和形狀有較大不同。本文給出了不同風(fēng)向角下干煤棚屋面各區(qū)域的劃分方法和確定每個(gè)區(qū)域經(jīng)驗(yàn)風(fēng)壓譜的四個(gè)參數(shù)。因此,僅需要個(gè)別點(diǎn)位的風(fēng)壓風(fēng)洞試驗(yàn)數(shù)據(jù)就可以歸納得出基本覆蓋整個(gè)干煤棚屋面的經(jīng)驗(yàn)風(fēng)壓譜,并可以根據(jù)相應(yīng)的經(jīng)驗(yàn)風(fēng)壓譜開展大跨干煤棚非高斯風(fēng)壓場(chǎng)模擬,獲得任意位置的風(fēng)壓時(shí)程數(shù)據(jù)樣本。

    (2) 數(shù)值模擬計(jì)算得到的測(cè)點(diǎn)脈動(dòng)風(fēng)壓時(shí)程在自功率譜密度上與風(fēng)洞試驗(yàn)譜差別不大,但相鄰風(fēng)壓測(cè)點(diǎn)之間的互功率譜密度與風(fēng)洞試驗(yàn)譜曲線還存在一定差距。而且不同風(fēng)向角下的模擬結(jié)果精度也有所不同,這個(gè)主要是由于不同風(fēng)向角下沿著不同方向上屋面風(fēng)壓場(chǎng)相干性的不同而造成的。本文采用的風(fēng)壓互功率譜其相干系數(shù)取值均為3.5,可以在后續(xù)工作中考慮取不同的相關(guān)系數(shù)來反映不同風(fēng)向角下風(fēng)壓場(chǎng)相干性的差別。

    (3) 在衡量脈動(dòng)風(fēng)壓場(chǎng)非高斯特性的高階統(tǒng)計(jì)量上,無論是偏度還是峰度指標(biāo),對(duì)比已有方法,本文方法模擬所得的多點(diǎn)風(fēng)壓時(shí)程數(shù)據(jù)結(jié)果與風(fēng)洞試驗(yàn)?zāi)繕?biāo)值吻合程度都較好,說明本文方法在非高斯風(fēng)壓場(chǎng)模擬方面具有較高的精度。

    (4) 本文算例利用某大跨干煤棚風(fēng)洞試驗(yàn)數(shù)據(jù)歸納得到了干煤棚表面的風(fēng)壓譜分區(qū)和參數(shù)取值,完成了大跨干煤棚非高斯風(fēng)壓場(chǎng)的數(shù)值模擬。但本文方法并不僅僅適用于大跨干煤棚。利用其他大跨建筑的風(fēng)洞試驗(yàn)數(shù)據(jù)或相關(guān)資料可以對(duì)本文采用的干煤棚表面風(fēng)壓經(jīng)驗(yàn)譜分區(qū)和各區(qū)域風(fēng)壓譜參數(shù)進(jìn)行修改,這樣該方法即可推廣應(yīng)用至其他大跨屋蓋建筑非高斯風(fēng)壓場(chǎng)的模擬。

    猜你喜歡
    角下峰度風(fēng)洞試驗(yàn)
    建筑物對(duì)塔機(jī)順風(fēng)向風(fēng)力干擾效應(yīng)研究
    擴(kuò)散峰度成像技術(shù)檢測(cè)急性期癲癇大鼠模型的成像改變
    磁共振擴(kuò)散峰度成像在肝臟病變中的研究進(jìn)展
    兇手是A角
    基于自動(dòng)反相校正和峰度值比較的探地雷達(dá)回波信號(hào)去噪方法
    低風(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)
    磁共振擴(kuò)散峰度成像MK值、FA值在鑒別高級(jí)別膠質(zhì)瘤與轉(zhuǎn)移瘤的價(jià)值分析
    遮擋條件下超高層建筑風(fēng)洞試驗(yàn)研究
    重慶建筑(2014年12期)2014-07-24 14:00:32
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    此物有八面人人有两片| 法律面前人人平等表现在哪些方面| 在线国产一区二区在线| 免费在线观看影片大全网站| 精品乱码久久久久久99久播| 久久亚洲真实| 成人亚洲精品av一区二区| 国产亚洲精品综合一区在线观看 | 大型黄色视频在线免费观看| 国产黄a三级三级三级人| 午夜福利视频1000在线观看| 亚洲男人的天堂狠狠| 无限看片的www在线观看| 男人舔奶头视频| 天堂动漫精品| 亚洲国产中文字幕在线视频| 亚洲在线自拍视频| 成人一区二区视频在线观看| 欧美乱色亚洲激情| 国产亚洲av高清不卡| 男女午夜视频在线观看| 久久狼人影院| 亚洲中文av在线| 男女视频在线观看网站免费 | avwww免费| 岛国视频午夜一区免费看| 欧美日韩乱码在线| 欧美乱妇无乱码| 亚洲人成伊人成综合网2020| 午夜久久久在线观看| 久久国产精品人妻蜜桃| 男男h啪啪无遮挡| 99精品在免费线老司机午夜| 成人一区二区视频在线观看| 黑人操中国人逼视频| 国产精品国产高清国产av| 在线观看66精品国产| 日韩免费av在线播放| 亚洲第一av免费看| 免费av毛片视频| 久久久久久大精品| 欧美zozozo另类| 欧美乱妇无乱码| 精品久久蜜臀av无| 最新美女视频免费是黄的| 亚洲在线自拍视频| 国产乱人伦免费视频| 日韩免费av在线播放| 国产野战对白在线观看| 午夜福利欧美成人| 国产精品二区激情视频| 亚洲精品美女久久久久99蜜臀| 亚洲国产欧美一区二区综合| 九色国产91popny在线| 国产亚洲精品久久久久久毛片| 婷婷六月久久综合丁香| 国产成人精品久久二区二区91| 婷婷丁香在线五月| 久久国产精品影院| 亚洲男人天堂网一区| 午夜福利视频1000在线观看| 黄片播放在线免费| 精品福利观看| 国产一级毛片七仙女欲春2 | 久久中文字幕人妻熟女| 91av网站免费观看| 十八禁人妻一区二区| av片东京热男人的天堂| 久久人妻av系列| 老司机深夜福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 久久人人精品亚洲av| 国产欧美日韩一区二区三| 91成年电影在线观看| 国产精品久久久久久人妻精品电影| 大香蕉久久成人网| 欧美中文日本在线观看视频| 好看av亚洲va欧美ⅴa在| 午夜影院日韩av| 欧美性猛交╳xxx乱大交人| 在线观看免费午夜福利视频| 亚洲自偷自拍图片 自拍| 国产精品久久久av美女十八| 老鸭窝网址在线观看| 色尼玛亚洲综合影院| 欧美国产精品va在线观看不卡| 日韩成人在线观看一区二区三区| 视频在线观看一区二区三区| 国内精品久久久久精免费| 免费在线观看影片大全网站| 黑丝袜美女国产一区| 国产亚洲欧美98| 无限看片的www在线观看| 美女大奶头视频| 日日干狠狠操夜夜爽| 成人亚洲精品av一区二区| 美女国产高潮福利片在线看| 香蕉av资源在线| 给我免费播放毛片高清在线观看| 亚洲欧美日韩高清在线视频| ponron亚洲| 亚洲av日韩精品久久久久久密| 女警被强在线播放| 欧美 亚洲 国产 日韩一| 国产成人系列免费观看| 一边摸一边抽搐一进一小说| 啦啦啦韩国在线观看视频| 亚洲全国av大片| 久久久国产精品麻豆| 国产视频内射| 国产黄片美女视频| 亚洲国产欧洲综合997久久, | 精品人妻1区二区| 久久中文字幕人妻熟女| svipshipincom国产片| 国产精品一区二区精品视频观看| 免费高清视频大片| 在线观看免费午夜福利视频| 香蕉久久夜色| 丰满的人妻完整版| 女警被强在线播放| www日本在线高清视频| 国产av在哪里看| 亚洲精品国产精品久久久不卡| 精品久久久久久久末码| 中文字幕精品亚洲无线码一区 | 国产精品一区二区精品视频观看| 此物有八面人人有两片| 1024香蕉在线观看| 日韩欧美免费精品| 一本一本综合久久| 特大巨黑吊av在线直播 | 久久精品aⅴ一区二区三区四区| 黄色女人牲交| 视频区欧美日本亚洲| 99久久无色码亚洲精品果冻| 久久九九热精品免费| 人人妻人人澡人人看| 欧美成人免费av一区二区三区| 精品第一国产精品| 真人做人爱边吃奶动态| 欧美精品亚洲一区二区| 国产男靠女视频免费网站| 亚洲成人国产一区在线观看| 自线自在国产av| 十八禁网站免费在线| 午夜激情福利司机影院| 欧美激情久久久久久爽电影| 国产av不卡久久| 国产97色在线日韩免费| 在线观看66精品国产| 国产精品一区二区三区四区久久 | 一级黄色大片毛片| 国产三级黄色录像| 人妻久久中文字幕网| 在线观看免费视频日本深夜| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| 99国产精品一区二区三区| 国产av一区二区精品久久| 制服诱惑二区| www日本在线高清视频| 在线永久观看黄色视频| 无人区码免费观看不卡| 又紧又爽又黄一区二区| 90打野战视频偷拍视频| 级片在线观看| 91成人精品电影| 99国产精品一区二区蜜桃av| 我的亚洲天堂| 在线观看一区二区三区| 在线观看免费日韩欧美大片| 9191精品国产免费久久| 日韩欧美一区视频在线观看| 在线观看免费午夜福利视频| 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| 久久久久久亚洲精品国产蜜桃av| 亚洲精品在线美女| 18禁美女被吸乳视频| 九色国产91popny在线| 韩国精品一区二区三区| 久久久久久久精品吃奶| 久久久久久久久中文| 99re在线观看精品视频| 久久久久国内视频| 搡老熟女国产l中国老女人| 国内揄拍国产精品人妻在线 | 免费无遮挡裸体视频| 午夜福利一区二区在线看| 欧美人与性动交α欧美精品济南到| 久久久久免费精品人妻一区二区 | 男人操女人黄网站| 欧美久久黑人一区二区| 亚洲午夜理论影院| 亚洲国产欧美日韩在线播放| 亚洲国产精品sss在线观看| 国产亚洲欧美在线一区二区| 法律面前人人平等表现在哪些方面| 一进一出抽搐gif免费好疼| 欧美日韩亚洲综合一区二区三区_| 美女扒开内裤让男人捅视频| 久久精品人妻少妇| 一区二区三区精品91| 国产亚洲精品av在线| 亚洲片人在线观看| 日韩欧美免费精品| 亚洲成a人片在线一区二区| 成人欧美大片| 亚洲成人免费电影在线观看| 一级a爱片免费观看的视频| 久久中文字幕一级| 久久久精品国产亚洲av高清涩受| 看片在线看免费视频| 成年版毛片免费区| 午夜免费成人在线视频| 免费高清视频大片| 国产v大片淫在线免费观看| 亚洲色图av天堂| 亚洲久久久国产精品| 亚洲欧美激情综合另类| 欧美成人午夜精品| 男女之事视频高清在线观看| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 欧美丝袜亚洲另类 | 国产av一区在线观看免费| 亚洲专区中文字幕在线| 免费高清在线观看日韩| 免费看日本二区| 国产精品1区2区在线观看.| 欧美av亚洲av综合av国产av| 亚洲精品久久国产高清桃花| 国产伦一二天堂av在线观看| 国产精品一区二区三区四区久久 | 亚洲专区中文字幕在线| 午夜免费激情av| 亚洲专区字幕在线| 女性生殖器流出的白浆| 日韩欧美国产在线观看| 亚洲国产精品sss在线观看| 午夜免费激情av| 国产精品亚洲美女久久久| 国产一卡二卡三卡精品| 亚洲一区高清亚洲精品| 一区二区三区高清视频在线| 十八禁人妻一区二区| 国产亚洲精品综合一区在线观看 | 18美女黄网站色大片免费观看| 一级毛片女人18水好多| 国产欧美日韩一区二区三| 国产视频内射| 免费在线观看影片大全网站| 中亚洲国语对白在线视频| 99国产精品99久久久久| 欧美黄色片欧美黄色片| 91麻豆av在线| 日韩精品青青久久久久久| a级毛片在线看网站| 国产成人精品久久二区二区免费| 日韩三级视频一区二区三区| 在线播放国产精品三级| 欧美在线黄色| 午夜老司机福利片| 欧美三级亚洲精品| 99国产精品一区二区蜜桃av| 欧美日本视频| 国产午夜精品久久久久久| 丰满人妻熟妇乱又伦精品不卡| 女人高潮潮喷娇喘18禁视频| 中文字幕人妻熟女乱码| 在线观看免费视频日本深夜| 亚洲av第一区精品v没综合| 亚洲在线自拍视频| 国产蜜桃级精品一区二区三区| 在线观看免费午夜福利视频| 国产亚洲av高清不卡| 午夜亚洲福利在线播放| 久久久久久久久中文| 亚洲av美国av| 成人18禁在线播放| 精品国产一区二区三区四区第35| 97超级碰碰碰精品色视频在线观看| 亚洲熟女毛片儿| 日韩欧美免费精品| 一级毛片高清免费大全| 亚洲国产欧洲综合997久久, | 国产精品二区激情视频| 国产精品永久免费网站| 精华霜和精华液先用哪个| 麻豆成人午夜福利视频| 亚洲aⅴ乱码一区二区在线播放 | 无遮挡黄片免费观看| 午夜免费激情av| 欧美激情久久久久久爽电影| 制服丝袜大香蕉在线| 免费一级毛片在线播放高清视频| 国产91精品成人一区二区三区| 女人被狂操c到高潮| 老汉色∧v一级毛片| 亚洲一码二码三码区别大吗| 国产亚洲精品一区二区www| 久久中文看片网| 亚洲一区高清亚洲精品| 免费在线观看影片大全网站| www.www免费av| 岛国在线观看网站| 性色av乱码一区二区三区2| 啦啦啦观看免费观看视频高清| 男人舔女人下体高潮全视频| 欧美黄色淫秽网站| 波多野结衣高清无吗| 青草久久国产| 国产亚洲av高清不卡| 久久中文字幕一级| 国产亚洲精品av在线| 亚洲中文av在线| 最新美女视频免费是黄的| 久久国产精品男人的天堂亚洲| 国产成人一区二区三区免费视频网站| 999精品在线视频| 色精品久久人妻99蜜桃| 欧美日韩精品网址| 亚洲精品中文字幕一二三四区| 一区二区三区国产精品乱码| 久久久水蜜桃国产精品网| 国产精品久久电影中文字幕| 18禁黄网站禁片免费观看直播| 日韩一卡2卡3卡4卡2021年| 999久久久精品免费观看国产| 69av精品久久久久久| 国产精品亚洲美女久久久| videosex国产| 一级毛片精品| 在线天堂中文资源库| 黄色丝袜av网址大全| tocl精华| 黑人巨大精品欧美一区二区mp4| 母亲3免费完整高清在线观看| 好看av亚洲va欧美ⅴa在| 欧美性猛交黑人性爽| 日日夜夜操网爽| 免费av毛片视频| 国产激情久久老熟女| 亚洲欧美日韩无卡精品| 国产单亲对白刺激| 亚洲美女黄片视频| 国内精品久久久久久久电影| 亚洲av成人不卡在线观看播放网| 免费女性裸体啪啪无遮挡网站| www.精华液| 亚洲电影在线观看av| 成人亚洲精品一区在线观看| 桃色一区二区三区在线观看| 日本熟妇午夜| 国产精品九九99| 国产在线精品亚洲第一网站| 亚洲,欧美精品.| 欧美另类亚洲清纯唯美| 一区福利在线观看| 午夜久久久在线观看| 男女下面进入的视频免费午夜 | 悠悠久久av| 午夜免费激情av| 国产精品99久久99久久久不卡| 制服丝袜大香蕉在线| 在线国产一区二区在线| 久久久久久大精品| 亚洲精品av麻豆狂野| 久久中文看片网| 伦理电影免费视频| 国产私拍福利视频在线观看| 两个人看的免费小视频| 国产一卡二卡三卡精品| 久久午夜综合久久蜜桃| 一区二区日韩欧美中文字幕| svipshipincom国产片| 日韩 欧美 亚洲 中文字幕| 性欧美人与动物交配| 免费看a级黄色片| 国产日本99.免费观看| 色婷婷久久久亚洲欧美| 国产精品久久久av美女十八| 伊人久久大香线蕉亚洲五| or卡值多少钱| 精品日产1卡2卡| 熟女电影av网| 91成年电影在线观看| 精品欧美国产一区二区三| 最近在线观看免费完整版| 9191精品国产免费久久| 免费看a级黄色片| www.www免费av| 亚洲国产日韩欧美精品在线观看 | 亚洲第一av免费看| 久99久视频精品免费| 国产99白浆流出| 亚洲黑人精品在线| 国产精品99久久99久久久不卡| 夜夜爽天天搞| 亚洲欧美精品综合一区二区三区| 夜夜爽天天搞| 少妇 在线观看| av超薄肉色丝袜交足视频| 婷婷精品国产亚洲av| 啪啪无遮挡十八禁网站| 国产黄片美女视频| 色精品久久人妻99蜜桃| 久久精品91蜜桃| 大香蕉久久成人网| 亚洲精品国产区一区二| 操出白浆在线播放| av中文乱码字幕在线| 很黄的视频免费| 午夜福利在线在线| 大香蕉久久成人网| 国产精品亚洲美女久久久| 久久天堂一区二区三区四区| 国产色视频综合| 亚洲精品av麻豆狂野| 国产又爽黄色视频| 成人手机av| 美女国产高潮福利片在线看| 亚洲精品美女久久久久99蜜臀| 欧美中文日本在线观看视频| 国产高清视频在线播放一区| 亚洲av电影在线进入| 久久热在线av| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲成人国产一区在线观看| 非洲黑人性xxxx精品又粗又长| av免费在线观看网站| 欧美激情 高清一区二区三区| 国产极品粉嫩免费观看在线| 欧美激情久久久久久爽电影| 久久久水蜜桃国产精品网| 成人三级做爰电影| 久久久久久人人人人人| 人成视频在线观看免费观看| 精品国产乱子伦一区二区三区| 亚洲国产毛片av蜜桃av| 亚洲色图 男人天堂 中文字幕| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| 日本在线视频免费播放| 久久久国产欧美日韩av| 欧美人与性动交α欧美精品济南到| 少妇 在线观看| 国产黄色小视频在线观看| 色播在线永久视频| 丝袜美腿诱惑在线| 一级a爱视频在线免费观看| 18美女黄网站色大片免费观看| 巨乳人妻的诱惑在线观看| x7x7x7水蜜桃| 日韩 欧美 亚洲 中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 在线永久观看黄色视频| 免费在线观看日本一区| 欧美日韩乱码在线| 一二三四在线观看免费中文在| 桃红色精品国产亚洲av| 91字幕亚洲| 少妇的丰满在线观看| 日韩精品免费视频一区二区三区| 两个人看的免费小视频| 亚洲人成77777在线视频| 熟妇人妻久久中文字幕3abv| 久久精品91无色码中文字幕| 妹子高潮喷水视频| e午夜精品久久久久久久| 精华霜和精华液先用哪个| 午夜成年电影在线免费观看| 国产黄a三级三级三级人| 亚洲国产日韩欧美精品在线观看 | 黑人巨大精品欧美一区二区mp4| tocl精华| 国产亚洲精品第一综合不卡| 欧美国产日韩亚洲一区| 亚洲国产欧洲综合997久久, | 夜夜爽天天搞| 999精品在线视频| 免费在线观看亚洲国产| 日韩精品免费视频一区二区三区| 亚洲人成网站高清观看| 黄色视频,在线免费观看| 久久久久国产精品人妻aⅴ院| 日韩高清综合在线| 天天躁夜夜躁狠狠躁躁| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产清高在天天线| 国产1区2区3区精品| 日日干狠狠操夜夜爽| 国产精品亚洲一级av第二区| 一进一出抽搐动态| 一级a爱片免费观看的视频| 国产又黄又爽又无遮挡在线| videosex国产| 欧美不卡视频在线免费观看 | 免费在线观看日本一区| 久久人妻福利社区极品人妻图片| 少妇 在线观看| 国产麻豆成人av免费视频| 国产精品久久久av美女十八| 一本精品99久久精品77| or卡值多少钱| 国产精品98久久久久久宅男小说| 国产极品粉嫩免费观看在线| 亚洲第一av免费看| 精品国产乱子伦一区二区三区| 国产精品精品国产色婷婷| 老熟妇乱子伦视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 天天躁狠狠躁夜夜躁狠狠躁| 最新美女视频免费是黄的| or卡值多少钱| 午夜激情av网站| 国产av又大| 国产在线观看jvid| 亚洲av电影不卡..在线观看| 久久狼人影院| 免费在线观看黄色视频的| 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| 1024手机看黄色片| 欧美亚洲日本最大视频资源| 久久精品91无色码中文字幕| 色婷婷久久久亚洲欧美| 99国产综合亚洲精品| 色播在线永久视频| 国产三级黄色录像| 日日爽夜夜爽网站| 亚洲精品美女久久久久99蜜臀| 夜夜躁狠狠躁天天躁| 免费电影在线观看免费观看| 精品久久久久久久久久久久久 | 国产亚洲欧美在线一区二区| 免费观看精品视频网站| 深夜精品福利| 91大片在线观看| 他把我摸到了高潮在线观看| 法律面前人人平等表现在哪些方面| 国产一区在线观看成人免费| 久久久国产成人精品二区| 午夜免费成人在线视频| 中文字幕精品亚洲无线码一区 | 国产精品香港三级国产av潘金莲| 可以免费在线观看a视频的电影网站| 男女做爰动态图高潮gif福利片| 欧美一级a爱片免费观看看 | 99热6这里只有精品| 一级a爱片免费观看的视频| 正在播放国产对白刺激| 久久久久久国产a免费观看| 天天一区二区日本电影三级| 麻豆国产av国片精品| 村上凉子中文字幕在线| 两个人看的免费小视频| 久久久精品欧美日韩精品| 午夜激情av网站| 中文字幕人妻熟女乱码| 国产精品香港三级国产av潘金莲| 一级毛片女人18水好多| 中文字幕另类日韩欧美亚洲嫩草| 一级黄色大片毛片| 人人妻人人澡人人看| 久久精品国产亚洲av香蕉五月| 中文字幕高清在线视频| 69av精品久久久久久| 妹子高潮喷水视频| 欧美一级毛片孕妇| 亚洲精品美女久久久久99蜜臀| 久久热在线av| 波多野结衣高清作品| 日日爽夜夜爽网站| 久久久水蜜桃国产精品网| 久久国产亚洲av麻豆专区| 琪琪午夜伦伦电影理论片6080| 麻豆一二三区av精品| 黄色视频不卡| 成人欧美大片| 国产亚洲欧美在线一区二区| 久久 成人 亚洲| 成人免费观看视频高清| 精品卡一卡二卡四卡免费| 国产精品自产拍在线观看55亚洲| 看片在线看免费视频| 久久久久久国产a免费观看| 一级毛片高清免费大全| 国产一区二区三区在线臀色熟女| 亚洲av成人一区二区三| 亚洲av电影在线进入| 亚洲午夜理论影院| 精品国产乱码久久久久久男人| 超碰成人久久| a级毛片在线看网站| 一级a爱视频在线免费观看| 国产成年人精品一区二区| 日本一本二区三区精品| 亚洲成av片中文字幕在线观看| 怎么达到女性高潮| 日本精品一区二区三区蜜桃| 我的亚洲天堂| 美女高潮到喷水免费观看| 人人妻人人看人人澡| 国产精品av久久久久免费| 亚洲av片天天在线观看| 日韩欧美免费精品| 欧美性长视频在线观看| 亚洲一码二码三码区别大吗| av免费在线观看网站| 久久久久九九精品影院| 精品国产超薄肉色丝袜足j| 欧美日韩一级在线毛片| 欧美+亚洲+日韩+国产|