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

    大型冷卻塔結(jié)構(gòu)風(fēng)致氣動(dòng)和氣彈效應(yīng)綜合研究

    2017-07-07 13:45:00展艷艷王志男梁譽(yù)文劉曉鵬程霄翔張軍鋒柯世堂王小松曹曙陽(yáng)葛耀君
    實(shí)驗(yàn)流體力學(xué) 2017年3期
    關(guān)鍵詞:趙林冷卻塔風(fēng)壓

    趙 林,展艷艷,王志男,梁譽(yù)文,劉曉鵬,程霄翔,張軍鋒,柯世堂,王小松,曹曙陽(yáng),葛耀君

    (同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200092)

    大型冷卻塔結(jié)構(gòu)風(fēng)致氣動(dòng)和氣彈效應(yīng)綜合研究

    趙 林,展艷艷*,王志男,梁譽(yù)文,劉曉鵬,程霄翔,張軍鋒,柯世堂,王小松,曹曙陽(yáng),葛耀君

    (同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200092)

    大型冷卻塔作為空間薄壁高聳結(jié)構(gòu),結(jié)構(gòu)風(fēng)效應(yīng)突出,在設(shè)計(jì)中風(fēng)荷載屬重要的控制因素。本文系統(tǒng)地介紹了同濟(jì)大學(xué)結(jié)構(gòu)風(fēng)工程研究團(tuán)隊(duì)在冷卻塔結(jié)構(gòu)抗風(fēng)方向取得的研究進(jìn)展,涉及風(fēng)荷載特征與統(tǒng)計(jì)、結(jié)構(gòu)風(fēng)效應(yīng)試驗(yàn)與分析、整體結(jié)構(gòu)多目標(biāo)優(yōu)化、災(zāi)害氣候條件荷載與效應(yīng)分析和全過(guò)程集成軟件平臺(tái)開(kāi)發(fā)等5個(gè)方面,在冷卻塔結(jié)構(gòu)抗風(fēng)問(wèn)題中的風(fēng)荷載、結(jié)構(gòu)響應(yīng)和試驗(yàn)?zāi)P驮O(shè)計(jì)等方面得出相應(yīng)結(jié)論。

    大型冷卻塔;風(fēng)荷載;結(jié)構(gòu)風(fēng)效應(yīng);多目標(biāo)優(yōu)化;災(zāi)害氣候;全過(guò)程分析

    0 引 言

    大型冷卻塔為空間薄壁高聳殼體結(jié)構(gòu),對(duì)風(fēng)荷載作用敏感,在設(shè)計(jì)中風(fēng)荷載是一個(gè)重要的控制因素,結(jié)構(gòu)風(fēng)荷載的準(zhǔn)確評(píng)價(jià)是其抗風(fēng)研究的基礎(chǔ)。結(jié)合風(fēng)洞試驗(yàn)和現(xiàn)場(chǎng)實(shí)測(cè)手段,研究了塔筒表面內(nèi)外壓分布規(guī)律[1-6],發(fā)現(xiàn)了超高雷諾數(shù)條件結(jié)構(gòu)表面動(dòng)態(tài)繞流發(fā)展規(guī)律[7-10],建立了基于計(jì)算流體力學(xué)方法的冷卻塔結(jié)構(gòu)風(fēng)荷載計(jì)算方法,揭示了復(fù)雜群塔干擾條件風(fēng)荷載分布模式和結(jié)構(gòu)靜、動(dòng)力荷載響應(yīng)規(guī)律[11],并對(duì)原型冷卻塔開(kāi)展了長(zhǎng)期動(dòng)態(tài)風(fēng)壓實(shí)測(cè)以驗(yàn)證試驗(yàn)和計(jì)算結(jié)果[12-14]。

    大型冷卻塔振型復(fù)雜,風(fēng)振效應(yīng)突出[15-17]。為從風(fēng)洞試驗(yàn)的角度研究結(jié)構(gòu)風(fēng)效應(yīng),提出了等效梁格氣動(dòng)彈性模型設(shè)計(jì)理論和方法[18]并開(kāi)展同步測(cè)壓和測(cè)振試驗(yàn)[19],用一致耦合分析法分析冷卻塔風(fēng)致響應(yīng)[20-21];建立了計(jì)入結(jié)構(gòu)荷載行為[22]的多目標(biāo)一致原則的等效風(fēng)荷載分析方法[2,23-25];考慮不同評(píng)價(jià)指標(biāo)[26-29],研究雙塔[17,30-31]、四塔[27]、六塔[26]和八塔[28,32]群塔組合干擾效應(yīng)[33-34];開(kāi)展了考慮非線性效應(yīng)的冷卻塔結(jié)構(gòu)屈曲穩(wěn)定分析。

    基于大型冷卻塔全壽命經(jīng)濟(jì)性,提出了考慮各種荷載模式下[35-36]結(jié)構(gòu)穩(wěn)定和強(qiáng)度安全等綜合效應(yīng)的塔筒選型多目標(biāo)優(yōu)化方法[37-39],優(yōu)化過(guò)程同時(shí)考慮復(fù)雜荷載組合作用下的結(jié)構(gòu)安全性和工程經(jīng)濟(jì)性。

    傳統(tǒng)冷卻塔研究工作中的風(fēng)洞試驗(yàn)和計(jì)算分析完全基于良態(tài)氣候模式(季風(fēng)等),大型冷卻塔對(duì)臺(tái)風(fēng)、龍卷風(fēng)等特殊極端風(fēng)作用非常敏感,此類風(fēng)環(huán)境特征不同于良態(tài)氣候模式。為此,開(kāi)展了臺(tái)風(fēng)和龍卷風(fēng)2類罕遇災(zāi)害模式作用下的理論、試驗(yàn)和實(shí)測(cè)研究[40-42],在臺(tái)風(fēng)脈動(dòng)特性[43-44]、極值風(fēng)特性隨機(jī)模擬[42]等方面取得進(jìn)展;分別利用多風(fēng)扇主動(dòng)控制風(fēng)洞[45-48]和自主開(kāi)發(fā)的龍卷風(fēng)模擬裝置[49]對(duì)臺(tái)風(fēng)和龍卷風(fēng)風(fēng)場(chǎng)特性進(jìn)行試驗(yàn)研究工作,保證了災(zāi)害氣候條件下冷卻塔環(huán)向風(fēng)壓分布、風(fēng)振系數(shù)及內(nèi)吸力影響等因素相關(guān)技術(shù)參數(shù)取值正確,滿足在我國(guó)沿海風(fēng)災(zāi)多發(fā)區(qū)建設(shè)大型冷卻塔的需求。

    基于冷卻塔設(shè)計(jì)預(yù)研、施工建造和整體運(yùn)營(yíng)等全壽命各個(gè)環(huán)節(jié),研發(fā)了涵蓋結(jié)構(gòu)建模、內(nèi)力計(jì)算、雙向配筋、整體優(yōu)化和自動(dòng)繪圖[50]等冷卻塔設(shè)計(jì)的軟件集成平臺(tái),并可針對(duì)多種災(zāi)害氣候模式(臺(tái)風(fēng)和龍卷風(fēng))分析風(fēng)效應(yīng)規(guī)律和抗風(fēng)穩(wěn)定性。

    1 風(fēng)荷載特征與統(tǒng)計(jì)

    1.1 動(dòng)態(tài)風(fēng)壓現(xiàn)場(chǎng)實(shí)測(cè)

    冷卻塔風(fēng)洞試驗(yàn)?zāi)P蜔o(wú)法完整模擬雷諾數(shù)效應(yīng),難于準(zhǔn)確再現(xiàn)冷卻塔表面動(dòng)態(tài)風(fēng)荷載與來(lái)流條件、塔群組合等參數(shù)間的關(guān)系,全尺寸現(xiàn)場(chǎng)實(shí)測(cè)是對(duì)風(fēng)洞試驗(yàn)的有力補(bǔ)充和驗(yàn)證。

    研發(fā)了全天候動(dòng)態(tài)風(fēng)壓采集設(shè)備,對(duì)徐州某電廠冷卻塔(高約167m)進(jìn)行了塔筒表面動(dòng)、靜態(tài)風(fēng)壓長(zhǎng)期現(xiàn)場(chǎng)觀測(cè)[12-13]。實(shí)測(cè)現(xiàn)場(chǎng)地處廣闊平原(可近似看作B類地貌),在正北與正東方向無(wú)任何高大建筑,在正西方向有少部分高大建筑(見(jiàn)圖1),考慮風(fēng)速風(fēng)向現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)并結(jié)合冷卻塔剛性模型測(cè)壓風(fēng)洞試驗(yàn)對(duì)比,規(guī)避了特定風(fēng)向臨近建筑物的干擾影響[13]。為確定風(fēng)壓測(cè)量時(shí)的來(lái)流狀態(tài),在測(cè)量塔正東向、距地面高度約10m處設(shè)置超聲波風(fēng)速儀。如圖2所示,在冷卻塔130m高度處沿塔筒外表面環(huán)向均勻布置36個(gè)壓阻式風(fēng)壓傳感器,在90m高度處外表面環(huán)向均勻布置18個(gè)壓阻式風(fēng)壓傳感器。測(cè)試期間記錄到了6次風(fēng)速較大(130m高度處風(fēng)速介于11~17m/s)且不受周邊高大建筑物干擾的風(fēng)壓時(shí)程。圖3給出了6次不同風(fēng)速下130m高度處沿環(huán)向表面脈動(dòng)壓力分布曲線及相應(yīng)的最小二乘平均值擬合曲線,同時(shí)與國(guó)內(nèi)外歷次冷卻塔脈動(dòng)風(fēng)壓的實(shí)測(cè)結(jié)果進(jìn)行對(duì)比。環(huán)向脈動(dòng)壓力分布曲線可以分為3個(gè)區(qū)域:迎風(fēng)區(qū)(0~40°)、側(cè)風(fēng)區(qū)(40°~120°)和背風(fēng)區(qū)(120°~180°)。迎風(fēng)區(qū)在駐點(diǎn)處的脈動(dòng)風(fēng)壓最大,在40°附近達(dá)到最小值,此點(diǎn)為平均壓力系數(shù)為0的角度;環(huán)向脈動(dòng)壓力的最大值出現(xiàn)在側(cè)風(fēng)區(qū)70°~90°的區(qū)域內(nèi),該點(diǎn)同時(shí)為平均壓力分布中最小負(fù)壓出現(xiàn)的角度;110°~120°為側(cè)風(fēng)區(qū)與尾流區(qū)的過(guò)渡區(qū)域,脈動(dòng)風(fēng)壓急劇減小,來(lái)流在此區(qū)域脫離;背風(fēng)區(qū)的風(fēng)壓脈動(dòng)平穩(wěn)且數(shù)值較小。

    1.2 塔筒表面風(fēng)荷載非高斯特性

    為了說(shuō)明不同斷面不同區(qū)域測(cè)點(diǎn)間風(fēng)壓的空間相關(guān)特性,圖4[25]給出了單體冷卻塔塔底、喉部和塔頂3個(gè)典型斷面迎風(fēng)區(qū)、負(fù)壓極值區(qū)和分離區(qū)測(cè)點(diǎn)風(fēng)壓與環(huán)向所有測(cè)點(diǎn)的相關(guān)性結(jié)果。由于塔底斷面存在明顯的三維端部效應(yīng),測(cè)點(diǎn)間的相關(guān)性較弱,主要呈現(xiàn)的是高斯分布特性;在中部和塔頂區(qū)域,迎風(fēng)面互相關(guān)性衰減極快,測(cè)點(diǎn)表現(xiàn)出高斯分布特性,在負(fù)壓極值至分離區(qū)域,測(cè)點(diǎn)間的相關(guān)性較強(qiáng),并存在明顯的大偏斜和高峰態(tài)現(xiàn)象,風(fēng)壓信號(hào)表現(xiàn)出非高斯特性[56],進(jìn)入背風(fēng)區(qū)其測(cè)點(diǎn)的相關(guān)性衰減,高斯特性又逐漸明顯。

    圖5給出了冷卻塔中上部斷面的環(huán)向高斯與非高斯區(qū)域的劃分圖,簡(jiǎn)單直接地反應(yīng)不同區(qū)域的風(fēng)壓分布特性[25]。圖中冷卻塔表面劃分高斯及非高斯區(qū)域的標(biāo)準(zhǔn)為:偏斜值|Cpisk|>0.2且峰態(tài)值|Cpiku|>3.5的風(fēng)壓信號(hào)為非高斯分布。

    1.3 統(tǒng)一脈動(dòng)風(fēng)壓公式

    根據(jù)脈動(dòng)風(fēng)壓現(xiàn)場(chǎng)實(shí)測(cè)和風(fēng)洞試驗(yàn)數(shù)據(jù)所表明的脈動(dòng)風(fēng)壓與紊流度之間的相關(guān)性,建立了不同紊流度條件下統(tǒng)一的脈動(dòng)風(fēng)壓環(huán)向分布模式,提出了考慮迎風(fēng)點(diǎn)脈動(dòng)風(fēng)壓系數(shù)修正擬合的統(tǒng)一脈動(dòng)風(fēng)壓系數(shù)均方值環(huán)向分布公式:

    式中:αk(Pu)為與迎風(fēng)點(diǎn)脈動(dòng)風(fēng)壓系數(shù)相關(guān)的八項(xiàng)式參數(shù);k=0~7;Pu是迎風(fēng)點(diǎn)脈動(dòng)風(fēng)壓系數(shù);Iu是來(lái)流順風(fēng)向紊流度。擬合參數(shù)如表1所示。

    表1 脈動(dòng)風(fēng)壓分布統(tǒng)一擬合式及擬合參數(shù)Table 1 Uniform fitting formula and fitting parameters of fluctuating wind pressure distribution

    為了進(jìn)行直觀的對(duì)比,將現(xiàn)有現(xiàn)場(chǎng)實(shí)測(cè)曲線按照迎風(fēng)點(diǎn)脈動(dòng)風(fēng)壓系數(shù)數(shù)值大小分為3組:第1組:徐州現(xiàn)場(chǎng)實(shí)測(cè),迎風(fēng)點(diǎn)脈動(dòng)系數(shù)約0.115,來(lái)流紊流度約7.32%;第2組:Ruscheweyh[51]/Davenport[52]/Sageau[53]現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果[57],迎風(fēng)點(diǎn)脈動(dòng)系數(shù)約0.2,來(lái)流紊流度約12.52%;第3組:周良茂現(xiàn)場(chǎng)實(shí)測(cè),迎風(fēng)點(diǎn)脈動(dòng)系數(shù)約0.245,來(lái)流紊流度約15.2%。分別采用迎風(fēng)點(diǎn)脈動(dòng)風(fēng)壓系數(shù),利用修正脈動(dòng)風(fēng)壓公式計(jì)算得到的結(jié)果,與上述3組現(xiàn)場(chǎng)實(shí)測(cè)脈動(dòng)風(fēng)壓系數(shù)進(jìn)行對(duì)比。圖6為徐州現(xiàn)場(chǎng)實(shí)測(cè)值與統(tǒng)一脈動(dòng)風(fēng)壓公式計(jì)算結(jié)果的對(duì)比,可以看出,徐州實(shí)測(cè)脈動(dòng)風(fēng)壓曲線與統(tǒng)一脈動(dòng)風(fēng)壓分布類似。圖7為已有現(xiàn)場(chǎng)實(shí)測(cè)值與統(tǒng)一脈動(dòng)風(fēng)壓系數(shù)的對(duì)比,Ruscheweyh和Davenport實(shí)測(cè)與統(tǒng)一脈動(dòng)風(fēng)壓系數(shù)在各個(gè)區(qū)域較為接近,而Sageau實(shí)測(cè)值在各個(gè)區(qū)域均大于其余測(cè)量值;周良茂實(shí)測(cè)脈動(dòng)風(fēng)壓曲線與統(tǒng)一脈動(dòng)風(fēng)壓系數(shù)較為接近,除了90°附近略有差異。綜上所述,除了Sageau實(shí)測(cè)脈動(dòng)風(fēng)壓系數(shù)與統(tǒng)一脈動(dòng)風(fēng)壓系數(shù)差別較大外,徐州冷卻塔脈動(dòng)風(fēng)壓實(shí)測(cè)結(jié)果與Davenport[52]、Ruscheweyh[51]、周良茂[54]實(shí)測(cè)曲線當(dāng)考慮來(lái)流紊流度效應(yīng)后與統(tǒng)一脈動(dòng)風(fēng)壓擬合結(jié)果在數(shù)值上差別不大,故可以定性地說(shuō)明來(lái)流紊流度的不同是造成實(shí)測(cè)脈動(dòng)風(fēng)壓系數(shù)產(chǎn)生差異的主要原因。結(jié)合本文推薦的脈動(dòng)風(fēng)壓曲線統(tǒng)一擬合式(即公式(1)和(2))解釋了圖3中所示歷史上多次脈動(dòng)風(fēng)壓分布存在差異的原因。

    2 結(jié)構(gòu)風(fēng)效應(yīng)試驗(yàn)與分析

    2.1 等效梁格氣彈模型設(shè)計(jì)

    氣彈模型風(fēng)洞試驗(yàn)是研究冷卻塔風(fēng)振響應(yīng)的有效手段,氣彈模型等效梁格設(shè)計(jì)方法針對(duì)傳統(tǒng)“連續(xù)介質(zhì)氣彈模型”在物理參數(shù)和氣動(dòng)力參數(shù)相似比模擬等方面存在的不足做出相應(yīng)的改進(jìn)(見(jiàn)表2),在模型頻率、振型、氣動(dòng)力參數(shù)模擬、風(fēng)振系數(shù)分布和試驗(yàn)易操作性方面均取得良好的效果[18]。等效梁格方法采用空間縱橫垂直交叉桁梁網(wǎng)狀結(jié)構(gòu)近似模擬連續(xù)殼筒體結(jié)構(gòu)動(dòng)力特性,解決了連續(xù)介質(zhì)模型彎扭剛度與軸向剛度不協(xié)調(diào)的問(wèn)題;按照幾何相似比的要求,采用具有可張拉性能的彈性、輕質(zhì)薄膜整體張貼在鋼骨架外表面模擬實(shí)際冷卻塔結(jié)構(gòu)的外形;為實(shí)現(xiàn)質(zhì)量系統(tǒng)的模擬,用銅鉛塊為配重補(bǔ)充不足部分的質(zhì)量;模型環(huán)向和子午向構(gòu)件采用A3鋼線切割加工,考慮外置彈性蒙皮附加阻尼比效應(yīng)后中,總阻尼比仍可以控制在3.5%左右,基本滿足鋼筋混凝土冷卻塔5%的阻尼比要求。根據(jù)等效梁格方法設(shè)計(jì)的冷卻塔氣彈模型如圖8所示。由表3可見(jiàn),等效梁格冷卻塔模型與原型結(jié)構(gòu)振型具有較好的相似性且表面壓力分布與規(guī)范值較好吻合。

    表2 冷卻塔氣彈模型等效梁格設(shè)計(jì)方法與傳統(tǒng)方法比較Table 2 Comparison between equivalent beam-net method and traditional method

    表3 冷卻塔氣彈模型設(shè)計(jì)參數(shù)與實(shí)測(cè)值Table 3 Design parameters and measurement values of cooling tower model

    2.2 共振響應(yīng)

    冷卻塔結(jié)構(gòu)的脈動(dòng)風(fēng)振響應(yīng)以共振分量為主[16],背景和交叉項(xiàng)分量的貢獻(xiàn)比較接近,一般均在10%以下,最大可達(dá)到20%左右。圖9給出了超大塔結(jié)構(gòu)總的脈動(dòng)風(fēng)振響應(yīng)三維分布圖。由圖可見(jiàn),隨著子午向高度的增大,總脈動(dòng)風(fēng)振響應(yīng)逐漸變大,在80m高度處出現(xiàn)第1個(gè)峰值,然后再減小到140m高度后又逐步增大至塔頂,其在子午向上出現(xiàn)2個(gè)明顯的峰值;在環(huán)向斷面上其響應(yīng)也出現(xiàn)多個(gè)峰值,并且沒(méi)有對(duì)稱特性;結(jié)構(gòu)的最大脈動(dòng)風(fēng)振響應(yīng)數(shù)值達(dá)到6.6mm。

    2.3 群塔干擾評(píng)價(jià)準(zhǔn)則

    為量化風(fēng)荷載條件冷卻塔群塔干擾效應(yīng),形成了眾多干擾效應(yīng)定義準(zhǔn)則,派生出數(shù)值結(jié)果差異性較大的荷載比例放大因子——群塔比例系數(shù),認(rèn)識(shí)和評(píng)判準(zhǔn)則存在明顯的差異。為了系統(tǒng)地研究冷卻塔復(fù)雜群塔條件荷載干擾效應(yīng),本項(xiàng)研究以某超大型冷卻塔為例,基于風(fēng)洞試驗(yàn)、結(jié)構(gòu)有限元分析和結(jié)構(gòu)設(shè)計(jì)配筋方案,系統(tǒng)分析了六塔典型布置條件的群塔干擾效應(yīng),詳細(xì)比較了不同塔距條件矩陣、菱形布置方案的基于荷載、受力和設(shè)計(jì)配筋3個(gè)比較準(zhǔn)則層面的群塔干擾效應(yīng)。圖10定義和比較了25種指標(biāo)下的群塔比例系數(shù),可見(jiàn)多種群塔比例系數(shù)準(zhǔn)則在數(shù)值大小和離散性方面存在明顯差異,在描述風(fēng)荷載最不利來(lái)流角度方面具有較好的一致性,其中配筋層面的群塔比例系數(shù)離散性顯著小于荷載層面結(jié)果[58]。

    根據(jù)配筋包絡(luò),采取分項(xiàng)群塔比例系數(shù)放大規(guī)范風(fēng)壓(見(jiàn)圖11),配筋曲線與按全局統(tǒng)一的群塔比例系數(shù)的配筋曲線的對(duì)比如圖12所示。由圖可知,按照分項(xiàng)群塔比例系數(shù)計(jì)算的塔筒配筋曲線在實(shí)現(xiàn)配筋包絡(luò)的同時(shí)減少其它模板位置的配筋富余[59]。推薦多風(fēng)向角多設(shè)計(jì)荷載組合獲得的塔筒配筋包絡(luò)作為冷卻塔設(shè)計(jì)參考依據(jù),基于配筋包絡(luò)比選的在塔筒高度范圍內(nèi)變化的分項(xiàng)干擾系數(shù)作為工程應(yīng)用群塔比例系數(shù)。

    2.4 多目標(biāo)等效風(fēng)荷載

    多目標(biāo)等效靜風(fēng)荷載的實(shí)質(zhì)是通過(guò)一定的數(shù)值優(yōu)化方法找到一個(gè)“誤差最小”的等效靜風(fēng)荷載分布模式,同時(shí)接近多個(gè)響應(yīng)目標(biāo)的極值。以喉部斷面環(huán)向18個(gè)節(jié)點(diǎn)(平均間隔20°取點(diǎn))位移響應(yīng)極值作為等效目標(biāo),得出了對(duì)應(yīng)的多目標(biāo)等效靜風(fēng)荷載三維分布圖和等值線圖,如圖13和14所示[25]。由圖可知:多目標(biāo)等效靜風(fēng)荷載沒(méi)有具體的物理含義,在數(shù)值上沒(méi)有出現(xiàn)分布很不合理的區(qū)域或現(xiàn)象,均在-15~18k N之間;在頂部斷面出現(xiàn)的峰值數(shù)量減小,并且數(shù)值也更加合理,與實(shí)際的等效靜風(fēng)荷載分布特征接近。

    為對(duì)比多目標(biāo)等效靜風(fēng)荷載和規(guī)范給出的單塔表面風(fēng)荷載分布特征的差別,圖15給出了冷卻塔喉部環(huán)向斷面上以典型結(jié)點(diǎn)位移響應(yīng)極值作為等效目標(biāo)得到的等效風(fēng)荷載和規(guī)范分布曲線圖,并給出了各個(gè)典型斷面上隨環(huán)向角度變化的等效風(fēng)荷載擬合多項(xiàng)式。由圖可知,考慮多目標(biāo)的等效風(fēng)荷載曲線在分布模式和數(shù)值大小上均與規(guī)范數(shù)值有較大的差別;在不同斷面分布模式不同,不存在絕對(duì)的正、負(fù)壓區(qū)域,并且其峰值也不一定出現(xiàn)在迎風(fēng)面,因此分?jǐn)嗝姹硎靖鼮楹侠怼?/p>

    為驗(yàn)證多目標(biāo)等效靜風(fēng)荷載的精度,將其作為靜力荷載加載到大型冷卻塔結(jié)構(gòu)上,求解相應(yīng)的風(fēng)振響應(yīng),并與結(jié)構(gòu)實(shí)際響應(yīng)極值進(jìn)行對(duì)比,對(duì)比結(jié)果如圖16所示[25]。由圖可見(jiàn),采用多目標(biāo)等效靜風(fēng)荷載計(jì)算出來(lái)的結(jié)果與實(shí)際響應(yīng)極值相差很小,誤差控制滿足要求。

    3 整體結(jié)構(gòu)多目標(biāo)優(yōu)化

    塔筒選型多目標(biāo)優(yōu)化方法考慮各種荷載組合下結(jié)構(gòu)穩(wěn)定、強(qiáng)度安全和總體造價(jià)等的綜合效應(yīng),結(jié)合響應(yīng)面法與梯度搜索法,得到了適用于多種風(fēng)荷載分布模式的塔型。結(jié)構(gòu)優(yōu)化選型過(guò)程分為2個(gè)階段[38]:第一階段確定在不同風(fēng)荷載分布模式下關(guān)于荷載加載方向角的最不利工況,優(yōu)化目標(biāo)為PT→PTmax;第二階段在對(duì)應(yīng)風(fēng)荷載分布模式下對(duì)典型結(jié)構(gòu)參量取值構(gòu)成的可行域內(nèi)進(jìn)行整體結(jié)構(gòu)優(yōu)化選出經(jīng)濟(jì)結(jié)構(gòu)尺寸,優(yōu)化目標(biāo)為PT→PTmin,待優(yōu)化參數(shù)如表4和圖17所示。

    表4 主要優(yōu)化變量Table 4 Main optimization variables

    選取規(guī)范對(duì)稱風(fēng)壓、試驗(yàn)迎風(fēng)非對(duì)稱風(fēng)壓、試驗(yàn)側(cè)風(fēng)非對(duì)稱風(fēng)壓和試驗(yàn)對(duì)稱放大風(fēng)壓4種風(fēng)荷載分布模式,在4種風(fēng)荷載分布模式下分別進(jìn)行冷卻塔結(jié)構(gòu)優(yōu)化,得到各荷載模式對(duì)應(yīng)的優(yōu)化塔型子午線型。優(yōu)化塔型子午線型的半徑及壁厚隨高度的變化如圖18所示[38]。

    從塔筒配筋的角度考察優(yōu)化塔型1的強(qiáng)度安全性,圖19[38]給出其在4種風(fēng)荷載分布模式下的塔筒理論配筋量。由圖可知,在規(guī)范對(duì)稱風(fēng)壓下,優(yōu)化塔型1的子午向外側(cè)、子午向內(nèi)側(cè)理論配筋量在絕大多數(shù)子午向模板都大于其他風(fēng)荷載分布模式的相應(yīng)值,環(huán)向外側(cè)、環(huán)向內(nèi)側(cè)理論配筋量至少70%的子午向模板理論配筋量都不小于其他風(fēng)荷載分布模式的相應(yīng)值,其余曲線交叉情況大多出現(xiàn)在第80階及以上或第5階及以下的子午向模板。在結(jié)構(gòu)設(shè)計(jì)中可通過(guò)實(shí)際局部配筋加強(qiáng)保證每一層子午向模板的強(qiáng)度安全。綜上所述,優(yōu)化后的塔筒子午向線形1適用于多種風(fēng)荷載分布模式。

    在規(guī)范對(duì)稱風(fēng)壓下,優(yōu)化塔型1的穩(wěn)定系數(shù)相較于初始塔型有所提高,鋼筋造價(jià)比降低5.4%,總造價(jià)降低19.0%。具體優(yōu)化效果隨優(yōu)化迭代步數(shù)的變化趨勢(shì)如圖20所示。

    表5 優(yōu)化塔型1與初始塔型性能對(duì)比Table 5 Comparison between optimized and initial tower

    4 災(zāi)害氣候條件荷載與效應(yīng)分析

    4.1 強(qiáng)臺(tái)風(fēng)環(huán)境主動(dòng)風(fēng)洞模擬

    由于氣象站在進(jìn)行風(fēng)速統(tǒng)計(jì)時(shí)無(wú)法區(qū)分極值風(fēng)速是良態(tài)風(fēng)引起還是臺(tái)風(fēng)引起,當(dāng)某地區(qū)良態(tài)風(fēng)的發(fā)生權(quán)重大于臺(tái)風(fēng)時(shí),就會(huì)降低規(guī)范所統(tǒng)計(jì)的基本風(fēng)速,若在該地區(qū)僅采用規(guī)范風(fēng)荷載進(jìn)行結(jié)構(gòu)抗風(fēng)設(shè)計(jì)而不考慮臺(tái)風(fēng)氣候的影響可能會(huì)導(dǎo)致結(jié)構(gòu)潛在的風(fēng)險(xiǎn)。

    為了模擬臺(tái)風(fēng)場(chǎng)極值風(fēng)環(huán)境特征,選取了上海地區(qū)臺(tái)風(fēng)極值風(fēng)環(huán)境作為主動(dòng)風(fēng)洞試驗(yàn)?zāi)M對(duì)象,本次試驗(yàn)利用了日本宮崎大學(xué)的三維多風(fēng)扇主動(dòng)控制風(fēng)洞(見(jiàn)圖21)。臺(tái)風(fēng)氣候下B類工程地貌的平均風(fēng)、紊流度、紊流積分尺度剖面模擬如圖22所示。

    以上海地區(qū)登陸的臺(tái)風(fēng)作為例,收集和整理了1949~2005年間該地臺(tái)風(fēng)年鑒統(tǒng)計(jì)資料和上海崇明島侯家鎮(zhèn)氣象站的逐時(shí)(連續(xù))氣象觀測(cè)記錄[41],通過(guò)敏感度分析方法選出了對(duì)臺(tái)風(fēng)極值風(fēng)環(huán)境貢獻(xiàn)程度較大的參數(shù),采用Monte-Carlo臺(tái)風(fēng)風(fēng)場(chǎng)隨機(jī)模型以越界峰值法和廣義Pareto分布探討了幾類工程場(chǎng)地目標(biāo)重現(xiàn)期內(nèi)極值風(fēng)速預(yù)測(cè)過(guò)程,給出了3類典型工程場(chǎng)地臺(tái)風(fēng)風(fēng)環(huán)境的梯度風(fēng)高度、平均風(fēng)剖面、極值風(fēng)速和陣風(fēng)因子取值(見(jiàn)表6)??梢钥闯?臺(tái)風(fēng)氣候A類、B類場(chǎng)地下,梯度風(fēng)高度普遍低于規(guī)范值,平均風(fēng)速剖面也較規(guī)范值陡;同時(shí)臺(tái)風(fēng)氣候下的10min平均風(fēng)速及瞬時(shí)風(fēng)速一般低于規(guī)范值。

    表6 3種典型工程場(chǎng)地臺(tái)風(fēng)與良態(tài)氣候模式風(fēng)環(huán)境參數(shù)比較(百年設(shè)計(jì)風(fēng)速)Table 6 Parameter comparison between typhoon and normal wind at 3 typical engineering sites

    基于某電廠大型冷卻塔現(xiàn)場(chǎng)實(shí)測(cè)和風(fēng)洞試驗(yàn),發(fā)現(xiàn)臺(tái)風(fēng)氣候下的紊流度Iu(TC)與規(guī)范風(fēng)荷載下的紊流度Iu(Code)存在如公式(3)所示的關(guān)系:

    式中:α值對(duì)應(yīng)A、B、C 3類工程場(chǎng)地分別取1.60、1.48和1.36;3類工程場(chǎng)地下的Iu(Code)分別 為IuA(z)=12% × (z/10)-0.12、IuB(z)=14% ×(z/10)-0.15、IuC(z)=23%×(z/10)-0.22。圖23給出了臺(tái)風(fēng)氣候與規(guī)范荷載下3類工程場(chǎng)地的紊流度剖面。圖24(a)和(b)分別給出了冷卻塔喉部區(qū)域與迎風(fēng)駐點(diǎn)的相關(guān)系數(shù)及與風(fēng)壓最小點(diǎn)的相關(guān)系數(shù),可以看出:臺(tái)風(fēng)氣候下紊流度高于規(guī)范風(fēng)荷載的紊流度,因此相關(guān)系數(shù)也較大,即紊流度越大,相關(guān)系數(shù)越強(qiáng)。

    4.2 龍卷風(fēng)條件

    同濟(jì)大學(xué)設(shè)計(jì)開(kāi)發(fā)的龍卷風(fēng)模擬器[49]可以通過(guò)改變導(dǎo)流板角度得到不同渦流比(不同尺度、不同流動(dòng)形態(tài))的龍卷風(fēng),通過(guò)改變風(fēng)機(jī)轉(zhuǎn)速得到不同轉(zhuǎn)速的龍卷風(fēng),同時(shí)可模擬龍卷風(fēng)水平移動(dòng)。龍卷風(fēng)風(fēng)洞及示意圖如圖25和26所示。由龍卷風(fēng)風(fēng)洞模擬器得到的龍卷風(fēng)三維風(fēng)速和風(fēng)壓場(chǎng)[61]如圖27所示。

    改變龍卷風(fēng)渦核與冷卻塔塔筒中心的相對(duì)位置關(guān)系,測(cè)量塔筒內(nèi)外表面風(fēng)壓分布。如圖28[62]所示,試驗(yàn)中龍卷風(fēng)和冷卻塔模型的相對(duì)位置關(guān)系通過(guò)r/rc量化,其中r代表冷卻塔中心與龍卷風(fēng)中心之間的距離,rc代表冷卻塔塔筒高度處的龍卷風(fēng)半徑;沿冷卻塔高度布置3層測(cè)壓點(diǎn),環(huán)向每層均勻布置12個(gè),測(cè)壓點(diǎn)的布置和編號(hào)如圖29所示。定義渦流比,變化范圍為1~12。

    以上部測(cè)點(diǎn)的平均風(fēng)壓系數(shù)分布為例,給出不同渦流比S、不同相對(duì)位置條件下的試驗(yàn)結(jié)果[62]如圖30、31和32所示。圖中橫坐標(biāo)為測(cè)壓點(diǎn)編號(hào),縱坐標(biāo)為平均壓力系數(shù)。由圖可見(jiàn),在龍卷風(fēng)特異氣流影響下冷卻塔內(nèi)、外表面均呈現(xiàn)負(fù)壓,與規(guī)范風(fēng)荷載與臺(tái)風(fēng)氣候下的冷卻塔表面風(fēng)荷載分布存在本質(zhì)區(qū)別;冷卻塔距離龍卷風(fēng)越近,塔筒內(nèi)、外表面吸力越大,反之,塔筒內(nèi)、外表面吸力越小。

    5 全過(guò)程集成軟件平臺(tái)開(kāi)發(fā)

    目前國(guó)內(nèi)采用的冷卻塔計(jì)算分析及設(shè)計(jì)程序難于滿足超出規(guī)范要求的超大型冷卻塔結(jié)構(gòu)抗風(fēng)和抗震設(shè)計(jì)要求,缺少靈活的風(fēng)洞試驗(yàn)與現(xiàn)場(chǎng)實(shí)測(cè)風(fēng)荷載與國(guó)外相關(guān)規(guī)范規(guī)定的荷載輸入條件,不具備整體結(jié)構(gòu)考慮多種荷載組合條件全局優(yōu)化設(shè)計(jì)功能,分析結(jié)果未與通用商業(yè)軟件進(jìn)行系統(tǒng)的校核,不方便設(shè)計(jì)人員完成冷卻塔結(jié)構(gòu)設(shè)計(jì)工作。為適應(yīng)我國(guó)新形式下超大型冷卻塔建設(shè)工作的順利展開(kāi),基于設(shè)計(jì)預(yù)研、施工建造和整體運(yùn)營(yíng)等全壽命各個(gè)環(huán)節(jié),研發(fā)了涵蓋結(jié)構(gòu)建模、內(nèi)力計(jì)算、雙向配筋、整體優(yōu)化和自動(dòng)繪圖等冷卻塔設(shè)計(jì)的軟件集成平臺(tái)——同濟(jì)風(fēng)向標(biāo)Wind Lock軟件,并可針對(duì)多種災(zāi)害氣候模式(臺(tái)風(fēng)和龍卷風(fēng))分析風(fēng)效應(yīng)規(guī)律和抗風(fēng)穩(wěn)定性。

    Wind Lock功能模塊主要包括:(1)良態(tài)與臺(tái)風(fēng)氣候極值風(fēng)環(huán)境模擬與預(yù)測(cè);(2)大跨空間結(jié)構(gòu)(冷卻塔)動(dòng)力及等效風(fēng)荷載分析;(3)大跨空間結(jié)構(gòu)(冷卻塔)多種荷載組合與設(shè)計(jì)分析;(4)輔助模塊 (氣動(dòng)力參數(shù)數(shù)據(jù)庫(kù)、復(fù)雜群樁特性分析等)。具體分項(xiàng)功能實(shí)現(xiàn)示例如圖33所示。

    6 結(jié) 語(yǔ)

    本文在概括介紹同濟(jì)大學(xué)結(jié)構(gòu)風(fēng)工程研究團(tuán)隊(duì)在冷卻塔結(jié)構(gòu)抗風(fēng)方向取得的研究進(jìn)展的基礎(chǔ)上,重點(diǎn)說(shuō)明了在風(fēng)荷載特征與統(tǒng)計(jì)、結(jié)構(gòu)風(fēng)效應(yīng)試驗(yàn)與分析、整體結(jié)構(gòu)多目標(biāo)優(yōu)化、災(zāi)害氣候條件荷載與效應(yīng)分析和全過(guò)程集成軟件平臺(tái)開(kāi)發(fā)等5個(gè)方面取得的代表性研究成果。主要結(jié)論如下:

    (1)進(jìn)行了長(zhǎng)期的原形冷卻塔動(dòng)態(tài)風(fēng)壓現(xiàn)場(chǎng)實(shí)測(cè),得到了環(huán)向脈動(dòng)風(fēng)壓分布曲線;發(fā)現(xiàn)了塔筒表面風(fēng)荷載的非高斯分布特性,給出了塔筒環(huán)向非高斯區(qū)域劃分;建立了考慮紊流度影響的統(tǒng)一脈動(dòng)風(fēng)壓擬合公式,對(duì)歷次現(xiàn)場(chǎng)實(shí)測(cè)值均有較好的擬合效果。

    (2)提出并應(yīng)用等效梁格氣彈模型研究冷卻塔結(jié)構(gòu)動(dòng)力響應(yīng),量化了沿塔筒高度脈動(dòng)響應(yīng)分布;比較不同指標(biāo)對(duì)群塔干擾效應(yīng)的評(píng)價(jià)異同;研究了多目標(biāo)等效靜風(fēng)荷載的分布規(guī)律及其合理性。

    (3)提出針對(duì)最不利工況確定、冷卻塔結(jié)構(gòu)參數(shù)選取和子午線性改進(jìn)方面的整體結(jié)構(gòu)多目標(biāo)優(yōu)化方法。

    (4)研究了臺(tái)風(fēng)、龍卷風(fēng)等災(zāi)害氣候條件下冷卻塔結(jié)構(gòu)的風(fēng)效應(yīng),對(duì)比良態(tài)風(fēng)條件及災(zāi)害氣候條件下冷卻塔內(nèi)、外表面風(fēng)荷載分布差異。

    (5)自主開(kāi)發(fā)了全過(guò)程集成軟件平臺(tái)Wind-Lock,實(shí)現(xiàn)大型冷卻塔結(jié)構(gòu)參數(shù)化建模、動(dòng)靜力分析、整體結(jié)構(gòu)優(yōu)化分析等。

    [1] 李鵬飛,趙林,葛耀君,等.超大型冷卻塔風(fēng)荷載特性風(fēng)洞試驗(yàn)研究[J].工程力學(xué),2008,25(6):60-67.Li P F,Zhao L,Ge Y J,et al.Wind tunnel investigation on wind load characteristics for super large cooling towers[J].Engineering Mechanics,2008,25(6):60-67.

    [2] 趙林,李鵬飛,葛耀君.等效靜風(fēng)荷載下超大型冷卻塔受力性能分析[J].工程力學(xué),2008,25(7):79-86.Zhao L,Li P F,Ge Y J.Numerical investigation on performance of super large cooling towers under equivalent static wind load[J].Engineering Mechanics,2008,25(7):79-86.

    [3]Zhao L,Ge Y J.Wind loading characteristics of super-large cooling towers[J].Wind and Structures,2010,13(3):257-273.

    [4]Cao S,Wang J,Cao J,et al.Experimental study of wind pressures acting on a cooling tower exposed to stationary tornadolike vortices[J].Journal of Wind Engineering and Industrial Aerodynamics,2015,145:75-86.

    [5]Cheng X X,Zhao L,Ge Y J.Multiple loading effects on windinduced static performance of super-large cooling towers[J].International Journal of Structural Stability and Dynamics,2013,13(8):1-13.

    [6] 李鵬飛.脈動(dòng)風(fēng)特性及其對(duì)橋梁主梁斷面的抖振作用研究[D].上海:同濟(jì)大學(xué),2007.Li P F.Research on turbulent wind characteristics and its effects on buffeting responses of bridge girder sections[D].Shanghai:Tongji University,2007.

    [7] 操金鑫,趙林,葛耀君,等.雙曲線圓截面建筑結(jié)構(gòu)雷諾數(shù)效應(yīng)模擬實(shí)踐[J].實(shí)驗(yàn)流體力學(xué),2009,23(4):46-50,55.Cao J X,Zhao L,Ge Y J,et al.Practice on simulation of Reynolds number effects for the hyperbolic circular section structures[J].Journal of Experiments in Fluid Mechanics,2009,23(4):46-50,55.

    [8] 董銳,趙林,葛耀君,等.雙曲圓截面冷卻塔壁面粗糙度對(duì)其繞流動(dòng)態(tài)特性影響[J].空氣動(dòng)力學(xué)學(xué)報(bào),2013,31(2):250-259.Dong R,Zhao L,Ge Y J,et al.Investigation of surface roughness and its influence to flow dynamic characteristics of hyperbolic cooling tower[J].Acta Aerodynamica Sinica,2013,31(2):250-259.

    [9] 劉天成,趙林,丁志斌.圓形截面冷卻塔不同表面粗糙度時(shí)繞流特性的試驗(yàn)研究[J].工業(yè)建筑,2006,36:301-304.Liu T C,Zhao L,Ding Z B.Test research of flow feature for hyperbolic circular section cooling tower with different superfacial roughness[J].Industrial Construction,2006,36:301-304.

    [10]Liu X P,Zhao L,Ge Y J.Extreme value distribution of surface aerodynamic pressure of hyperbolic cooling tower[C].The 7thinternational colloquium on bluff body areodynamics and applications,Shanghai,2012.

    [11]Zhang J F,Ge Y J,Zhao L.Influence of latitude wind pressure distribution on the responses of hyperbolodial cooling tower shell[J].Wind and Structures,2013,16(6):579-601.

    [12]Zhao L,Ge Y J,Ahsan K.Fluctuating Wind Pressure Distribution around Full-scale Cooling Tower[J].Journal of Wind Engineering&Industrial Aerodynamics,2017,165:34-45.

    [13]劉曉鵬.大型冷卻塔動(dòng)態(tài)風(fēng)壓現(xiàn)場(chǎng)實(shí)測(cè)及其效應(yīng)綜合評(píng)價(jià)[D].上海:同濟(jì)大學(xué),2013.Liu X P.On-spot measurement of stochastic wind pressure and its evaluation of multiple effects[D].Shanghai:Tongji University,2013.

    [14]Zhao L,Ge Y J.Dynamic wind pressure distribution around full size cooling tower[C].Proceedings of the International Association for Shell and Spatial Structures(IASS)Symposium 2015,2015.

    [15]Ke S T,Ge Y J.The influence of self-excited forces on wind loads and wind effects for super-large cooling towers[J].Journal of Wind Engineering and Industrial Aerodynamics,2014,132:125-135.

    [16]Ke S T,Ge Y J,Zhao L,et al.A new methodology for analysis of equivalent static wind loads on super-large cooling towers[J].Journal of Wind Engineering and Industrial Aerodynamics,2012,111:30-39.

    [17]許林汕,趙林,葛耀君.超大型冷卻塔隨機(jī)風(fēng)振響應(yīng)分析[J].振動(dòng)與沖擊,2009,28(4):180-184.Xu L S,Zhao L,Ge Y J.Wind-excited stochastic responses of super large cooling towers[J].Journal of Vibration and Shock,2009,28(4):180-184.

    [18]趙林,葛耀君,曹豐產(chǎn).雙曲薄殼冷卻塔氣彈模型的等效梁格方法和實(shí)驗(yàn)研究[J].振動(dòng)工程學(xué)報(bào),2008,21(1):31-37.Zhao L,Ge Y J,Cao F C.Equivalent beam-net design theory of aero-elastic model about hyperbolic thin-shell cooling towers and its experimental investigation[J].Journal of Vibration Engineering,2008,21(1):31-37.

    [19]柯世堂,侯憲安,趙林,等.超大型冷卻塔風(fēng)荷載和風(fēng)振響應(yīng)參數(shù)分析:自激力效應(yīng)[J].土木工程學(xué)報(bào),2012,45(12):45-53.Ke S T,Hou X A,Zhao L,et al.Parameter analysis of wind loads and wind induced responses forsuper-large cooling towers:self-excited force effect[J].China Civil Engineering Journal,2012,45(12):45-53.

    [20]Zhao L,Ke S T,Ge Y J,et al.Wind-induced performance of super large cooling tower considering coupling effects between resonant and background components[C].International Symposium on Industrial Chimneys and Cooling Towers.Prague,2014.

    [21]葛耀君,趙林.大型冷卻塔風(fēng)荷載共振效應(yīng)的耦合分析[C].全國(guó)結(jié)構(gòu)振動(dòng)與動(dòng)力學(xué)學(xué)術(shù)研討會(huì)暨第四屆結(jié)構(gòu)動(dòng)力學(xué)專業(yè)委員會(huì)會(huì)議,中國(guó) 蘇州,2011.

    [22]Ding Z B,Tamura Y,Yoshida A.Contributions to member stresses due to overall wind-induced behaviors of thin-walled cylindrical shell[J].Journal of Wind Engineering and Industrial Aerodynamics,2012,107:192-201.

    [23]柯世堂,陳少林,趙林,等.超大型冷卻塔等效靜風(fēng)荷載精細(xì)化計(jì)算及應(yīng)用[J].振動(dòng)測(cè)試與診斷,2013,33(5):824-830.Ke S T,Chen S L,Zhao L,et al.Fine calculation and application of equivalent static wind load of super large cooling tower[J].Journal of Vibration,Measurement&Diagnosis,2013,33(5):824-830.

    [24]Ke S T,Ge Y J,Zhao L,et al.Wind-induced responses of super-large cooling towers[J].Journal of Central South University,2013,20(11):3216-3228.

    [25]柯世堂.大型冷卻塔結(jié)構(gòu)風(fēng)效應(yīng)和等效靜風(fēng)荷載研究[D].上海:同濟(jì)大學(xué),2011.Ke S T.Investigation of structrual wind effects and equivalent static wind loads of large cooling towers[D].Shanghai:Tongji University,2011.

    [26]程霄翔.大型冷卻塔表面風(fēng)壓及其雷諾數(shù)效應(yīng)的風(fēng)洞試驗(yàn)和現(xiàn)場(chǎng)實(shí)測(cè)[D].上海:同濟(jì)大學(xué),2015.Cheng X X.Study of wind pressures and their Reynold’s number effects on large cooling towers using wind tunnel tests and field measurements[D].Shanghai:Tongji University,2015.

    [27]Zhao L,Chen X,Ke S T,et al.Aerodynamic and aero-elastic performances of super-large cooling towers[J].Wind and Structures,2014,19(4):443-465.

    [28]程霄翔,趙林,葛耀君.典型矩形八塔超大型冷卻塔塔群風(fēng)致干擾效應(yīng)試驗(yàn)[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,44(1):372-380.Cheng X X,Zhao L,Ge Y J.Wind tunnel investigation on interference effect of eight grouped super large cooling towers with rectangular arrangement[J].Journal of Central South U-niversity(Science and Technology),2013,44(1):372-380.

    [29]張軍鋒,葛耀君,趙林.群塔布置對(duì)冷卻塔整體風(fēng)荷載和風(fēng)致響應(yīng)的不同干擾效應(yīng)[J].工程力學(xué),2016,(8):15-23,44.Zhang J F,Ge Y J,Zhao L.Interference effects on global wind loads and wind induced responses for group hyperboloidal cooling towers[J].Engineering Mechanics,2016,(8):15-23,44.

    [30]趙林,葛耀君,許林汕,等.超大型冷卻塔風(fēng)致干擾效應(yīng)試驗(yàn)研究[J].工程力學(xué),2009,(1):149-154.Zhao L,Ge Y J,Xu L S.et al.Wind tunnel investigation on wind-induced interference effects for super large cooling towers[J].Engineering Mechanics,2009,(1):149-154.

    [31]張軍鋒,趙林,柯世堂,等.大型雙冷卻塔雙塔組合表面風(fēng)壓干擾效應(yīng)試驗(yàn)研究[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào),2011,43(4):81-87.Zhang J F,Zhao L,Ke S T,et al.Wind tunnel investigation on wind-pressure interference effects for two large hyperbolodial cooling towers[J].Journal of Harbin Institute of Technology,2011,43(4):81-87.

    [32]Zhao L,Zhan Y Y,Ge Y J.Unified wind-induced interference criteria based on the reinforcement envelope of cooling towers[C].International Symposium on Industrial Chimneys and Cooling Towers,Rotterdam,2016.

    [33]張軍鋒.大型冷卻塔結(jié)構(gòu)特性與風(fēng)致干擾效應(yīng)研究[D].上海:同濟(jì)大學(xué),2012.Zhang J F.Structure behavior and wind-induced interference effects of large cooling towers[D].Shanghai:Tongji University,2012.

    [34]Zhao L,Chen X,Ge Y J.Investigations of adverse wind loads on a large cooling tower for the six-tower combination[J].Applied Thermal Engineering,2016,105:988-999.

    [35]張軍鋒,葛耀君,柯世堂,等.中美日三國(guó)規(guī)范高層結(jié)構(gòu)風(fēng)荷載標(biāo)準(zhǔn)值對(duì)比[J].湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,38(10):18-25.Zhang J F,Ge Y J,Ke S T,et al.Comparative study on the nominal value of wind loads on tall buildings among the codes of China,America and Japan[J].Journal of Hunan University(Natural Sciences),2011,38(10):18-25.

    [36]張軍鋒,葛耀君,趙林.雙曲冷卻塔結(jié)構(gòu)規(guī)范風(fēng)荷載標(biāo)準(zhǔn)值對(duì)比[J].建筑科學(xué)與工程學(xué)報(bào),2011,28(2):52-61.Zhang J F,Ge Y J,Zhao L.Comparisons of nominal values of wind loads on hyperboloidal cooling towers[J].Journal of Architecture and Civil Engineering,2011,28(2):52-61.

    [37]張軍鋒,葛耀君,趙林.雙曲冷卻塔下部子午線形對(duì)結(jié)構(gòu)性能的影響[J].建筑結(jié)構(gòu),2012,(2):118-121.Zhang J F,Ge Y J,Zhao L.Effect of bottom meridian curve on the mechanics behavior of hyperboloidal cooling tower[J].Building Structure,2012,(2):118-121.

    [38]梁譽(yù)文.考慮多種風(fēng)荷載分布模式的冷卻塔結(jié)構(gòu)優(yōu)化選型[D].上海:同濟(jì)大學(xué),2016.Liang Y W.Structural optimization of large cooling towers considering various distribution patterns of wind loadings[D].Shanghai:Tongji University,2016.

    [39]Zhao L,Liang Y W,Ge Y J,et al.Multi-object optimazation analysis of structural design for large cooling towers[J].Heat Transfer engineering,2017,38(11-12):1135-1145.

    [40]趙林,葛耀君,項(xiàng)海帆.臺(tái)風(fēng)風(fēng)場(chǎng)隨機(jī)參數(shù)敏感性分析[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,33(6):727-731.Zhao L,Ge Y J,Xiang H F.Stochastic parameter sensitivity analysis of typhoon wind field[J].Journal of Tongji University(Natural Science),2005,33(6):727-731.

    [41]趙林,朱樂(lè)東,葛耀君.上海地區(qū)臺(tái)風(fēng)風(fēng)特性Monte-Carlo隨機(jī)模擬研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2009,(1):25-31.Zhao L,Zhu L D,Ge Y J.Monte-Carlo simulation about typhoon extreme value wind characteristics in Shanghai region[J].Acta Aerodynamica Sinica,2009,(1):25-31.

    [42]趙林,葛耀君,宋麗莉,等.廣州地區(qū)臺(tái)風(fēng)極值風(fēng)特性蒙特卡羅隨機(jī)模擬[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,08:1034-1038,1068.Zhao L,Ge Y J,Song L L,et al.Monte-Carlo simulation analysis of typhoon extreme value wind characteristics in Guangzhou[J].Journal of Tongji University(Natural Science),2007,08:1034-1038,1068.

    [43]武占科,趙林,朱樂(lè)東,等.“羅莎”(0716)臺(tái)風(fēng)高空實(shí)測(cè)脈動(dòng)風(fēng)特性分析[J].空氣動(dòng)力學(xué)學(xué)報(bào),2010,28(3):291-296.Wu Z K,Zhao L,Zhu L D,et al.High-altitude observation about turbulence characteristics for“Krosa”(0716)strong typhoon[J].Acta Aerodynamica Sinica,2010,28(3:)291-296.

    [44]潘晶晶,趙林,冀春曉,等.東南沿海登陸臺(tái)風(fēng)近地脈動(dòng)特性分析[J].建筑結(jié)構(gòu)學(xué)報(bào),2016,(1):85-90.Pan J J,Zhao L,Ji C X,et al.Fluctuating wind characteristics near flat ground during typhoon landfall in southeast coastal areas of China[J].Journal of Building Structures,2016,(1):85-90.

    [45]潘韜,趙林,曹曙陽(yáng),等.多風(fēng)扇主動(dòng)控制風(fēng)洞類平板斷面抖振力識(shí)別研究[J].振動(dòng)與沖擊,2010,29(6):178-183.Pan T,Zhao L,Cao S Y,et al.Buffeting force analysis of thin plate section in multiple fans active control wind tunnel[J].Journal of Vibration and Shock,2010,29(6):178-183.

    [46]許坤,趙林,葛耀君,等.主動(dòng)來(lái)流閉口箱梁氣動(dòng)力時(shí)頻相關(guān)性[J].同濟(jì)大學(xué)學(xué)報(bào)自然科學(xué)版,2014,42(3):338-345.Xu K,Zhao L,Ge Y J,et al.Time-frequency coherences of aerodynamic loads on a streamlined bridge girder in an actively controlled wind tunnel[J].Journal of Tongji University(Natural Science),2014,42(3):338-345.

    [47]潘韜,趙林,曹曙陽(yáng),等.主動(dòng)來(lái)流條件類平板斷面氣動(dòng)力荷載效應(yīng)分析[J].實(shí)驗(yàn)流體力學(xué),2010,24(6):32-37.Pan T,Zhao L,Cao S Y,et al.Analysis of aerodynamic load effects on thin plat section under active control flow condition[J].Journal of Experiments in Fluid Mechanics,2010,24(6):32-37.

    [48]趙楊,曹曙陽(yáng),Yukio Tamura,等.雷暴沖擊風(fēng)模擬及其荷載的風(fēng)洞試驗(yàn)研究[J].振動(dòng)與沖擊,2009,(4):1-3,9.Zhao Y,Cao S Y,Yukio T,et al.Simulation of downburst and its loads with wind tunnel test[J].Journal of Vibration and Shock,2009,(4):1-3,9.

    [49]王錦,周強(qiáng),曹曙陽(yáng),等.龍卷風(fēng)風(fēng)場(chǎng)的試驗(yàn)?zāi)M[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,11:1654-1659.Wang J,Zhou Q,Cao S Y,et al.Physical study on tornadolike flow based on tornado vortex simulator[J].Journal of Tongji University(Natural Science),2014,11:1654-1659.

    [50]鄭鴻飛,王小松,毛國(guó)輝.CAD二次開(kāi)發(fā)在冷卻塔軟件開(kāi)發(fā)中的應(yīng)用[J].計(jì)算機(jī)技術(shù)與發(fā)展,2013,23(10):188-192.Zheng H F,Wang X S,Mao G H.Application of CAD secondary development in cooling tower software development[J].Computer Technology and Development,2013,23(10):188-192.

    [51]Ruscheweyh H.Wind loadings on hyperbolic natural draft cooling towers[J].Journal of Wind Engineering and Industrial Aerodynamics,1976,1:335-340.

    [52]Davenport A G,Isyumov N.The dynamic and static action of wind on hyperbolic cooling towers[R].Rresearch report No.BLWT1-66,Univ of Western Ontario,London,Ontario,Canada,1966.

    [53]Sageau J F.Caractérisation des champs de pression moyens et fluctuantsàla surface des grands aérorefrigérants,Electricitéde France[C].Direction desétudes et Recherches,6 quai Watier,Chatou,France,1979.

    [54]周良茂,李培華.兩個(gè)鄰近全尺寸雙曲型冷卻塔風(fēng)壓分布的測(cè)量[J].氣動(dòng)實(shí)驗(yàn)與測(cè)量控制,1992,6(3):37-44.Zhou L M,Li P H.A measurement of the wind pressure distribution on two neighboring full-scale hyperbolic cooling towers[J].Aerodynamic Experiment and Measurement&Control,1992,6(3):37-44.

    [55]趙林,劉曉鵬,高玲,等.大型冷卻塔表面脈動(dòng)風(fēng)壓原型實(shí)測(cè)與分布準(zhǔn)則[J].土木工程學(xué)報(bào),2017,50(1):1-11.Zhao L,Liu X P,Gao L,et al.Full scale measurement and distribution rules of surface fluctuating wind pressure of a large cooling tower[J].China Civil Engineering Journal,2017,50(1):1-11.

    [56]柯世堂,趙林,葛耀君.大型雙曲冷卻塔表面脈動(dòng)風(fēng)壓隨機(jī)特性——風(fēng)壓極值探討[J].實(shí)驗(yàn)流體力學(xué),2010,24(4):7-12.Ke S T,Zhao L,Ge Y J.Features of fluctuating wind pressure on large hyperbolic cooling tower:discussions on extreme wind pressure[J].Journal of Experiments in Fluid Mechanics,2010,24(4):7-12.

    [57]Simiu E,Scanlan R H,Simiu E,et al.Wind effects on structures:funddementals and applications to design[M].New York:John Wiley,1996.

    [58]Zhao L,Zhan Y Y,Liang Y W,et al.Comparison of wind-induced interference criteria and effects under typical towersgroup arrangements[C]//Spatial Structures in the 21st Century,Proceedings of the IASS Annual Symposium 2016,Tokyo,Japan,2016.

    [59]Zhao L,Zhan Y Y,Ge Y J.Unified wind-induced interference criteria based on the reinforcement envelope of cooling towers[C].International Symposium on Industrial Chimneys and Cooling Towers,Rotterdam,Holland,2016.

    [60]Nishi A,Kikugawa H,Matsuda Y,et al.Turbulence control in multiple-fan wind tunnels[J].Journal of Wind Engineering&Industrial Aerodynamics,1997,67-68(4):861-872.

    [61]Wang J,Cao S,Pang W,et al.Experimental study on effects of ground roughness on flow characteristics of tornado-like vortices[J].Boundary-Layer Meteorology,2016:1-21.

    [62]Cao S,Wang J,Cao J,et al.Experimental study of wind pressures acting on a cooling tower exposed to stationary tornadolike vortices[J].Journal of Wind Engineering&Industrial Aerodynamics,2015,145:75-86.

    Comprehensive investigation into wind-induced aerodynamic and aeroelastic effects of large cooling towers

    Zhao Lin,Zhan Yanyan*,Wang Zhinan,Liang Yuwen,Liu Xiaopeng,Cheng Xiaoxiang,Zhang Junfeng,Ke Shitang,Wang Xiaosong,Cao Shuyang,Ge Yaojun
    (State Key Laboratory of Disaster Reduction in Civil Engineering,Tongji University,Shanghai 200092,China)

    As typical thin-wall flexible structures,large cooling towers are of critical sensitivity to wind loads during its life-cycle.It's introduced systematically the research development of large cooling towers made by Tongji wind-resistance group during the past decades,including the study on statistic characteristics of the wind load from the full-scale observation,experiment and analysis of wind-induced structural responses,multi-objective structural optimization,the study on wind-induced loads and effects under disaster climates,the development of the self-developed whole process integration software platform,and so on.Some key problems in engineering application and fundamental investigation discovered during the development of large cooling towers in China have been resolved successfully.

    large cooling towers;wind loads;wind-induced structural effects;multi-objective optimization;disaster climate;full-process analysis

    TU973.2+13

    :A

    (編輯:李金勇)

    2016-12-21;

    :2017-03-13

    國(guó)家自然科學(xué)基金項(xiàng)目(51678451和50978203);科技部重大科技項(xiàng)目(核電超大型冷卻塔結(jié)構(gòu)研究及技術(shù)支持,2009ZX06004-010-HYJY);國(guó)家自然科學(xué)基金委優(yōu)秀青年基金項(xiàng)目(51222809);教育部“新世紀(jì)優(yōu)秀人才支持計(jì)劃”(NCET-13-0429)聯(lián)合資助項(xiàng)目

    *通信作者E-mail:Zhan-yanyan@126.com

    Zhao L,Zhan Y Y,Wang Z N,et al.Comprehensive investigation into wind-induced aerodynamic and aeroelastic effects of large cooling towers.Journal of Experiments in Fluid Mechanics,2017,31(3):1-15.趙 林,展艷艷,王志男,等.大型冷卻塔結(jié)構(gòu)風(fēng)致氣動(dòng)和氣彈效應(yīng)綜合研究.實(shí)驗(yàn)流體力學(xué),2017,31(3):1-15.

    1672-9897(2017)03-0001-15

    10.11729/syltlx20160201

    趙 林(1974-),黑龍江牡丹江人,教授,博導(dǎo)。研究方向:橋梁和建筑結(jié)構(gòu)風(fēng)工程。通信地址:上海市四平路1239號(hào)同濟(jì)大學(xué)風(fēng)工程館305室(200092)。E-mail:zhaolin@#edu.cn

    猜你喜歡
    趙林冷卻塔風(fēng)壓
    巧妙的接頭
    “蔥油餅大姐”趙林的別樣人生
    愉快的一天
    昆崳(2018年2期)2018-10-09 07:42:04
    鋼結(jié)構(gòu)冷卻塔風(fēng)荷載的模擬和計(jì)算
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    空氣冷卻塔除沫器改造與管理
    低風(fēng)壓導(dǎo)線的發(fā)展和測(cè)試
    電線電纜(2016年5期)2016-02-27 09:02:06
    西門(mén)子S7-200PLC在火電廠冷卻塔風(fēng)機(jī)控制中的應(yīng)用
    西藏科技(2015年3期)2015-09-26 12:11:11
    “蔥油餅大姐”的小說(shuō)夢(mèng)
    高層建筑風(fēng)荷載有哪些計(jì)算要求
    重慶建筑(2014年10期)2014-03-27 02:20:07
    国产精品九九99| 法律面前人人平等表现在哪些方面| 人人妻人人添人人爽欧美一区卜| 欧美精品亚洲一区二区| 亚洲伊人色综图| 亚洲第一青青草原| 色精品久久人妻99蜜桃| 大码成人一级视频| 咕卡用的链子| 村上凉子中文字幕在线| 大香蕉久久网| 99精品久久久久人妻精品| 久久久久精品国产欧美久久久| 一本大道久久a久久精品| 久久久久久免费高清国产稀缺| av欧美777| 日本黄色日本黄色录像| 久久久久国内视频| 伦理电影免费视频| aaaaa片日本免费| 午夜免费观看网址| 人妻久久中文字幕网| 夜夜夜夜夜久久久久| 18禁黄网站禁片午夜丰满| 99精品久久久久人妻精品| 久9热在线精品视频| 精品熟女少妇八av免费久了| 熟女少妇亚洲综合色aaa.| 91在线观看av| 国产精品免费视频内射| 美女国产高潮福利片在线看| 国产成人免费观看mmmm| 久久久国产一区二区| 免费av中文字幕在线| 中文字幕最新亚洲高清| 亚洲 欧美一区二区三区| 丝袜美足系列| 日韩免费高清中文字幕av| 久久精品亚洲熟妇少妇任你| 久久国产精品影院| 乱人伦中国视频| 日韩欧美三级三区| 性色av乱码一区二区三区2| 成年人黄色毛片网站| 18禁黄网站禁片午夜丰满| 国产精品免费大片| xxx96com| 黄色毛片三级朝国网站| 欧美激情久久久久久爽电影 | 亚洲精品一二三| 大香蕉久久网| 亚洲七黄色美女视频| 久9热在线精品视频| 国产单亲对白刺激| 50天的宝宝边吃奶边哭怎么回事| 看片在线看免费视频| 变态另类成人亚洲欧美熟女 | 亚洲成国产人片在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲av片天天在线观看| 亚洲专区国产一区二区| av线在线观看网站| 亚洲国产毛片av蜜桃av| 亚洲,欧美精品.| 极品教师在线免费播放| 这个男人来自地球电影免费观看| av福利片在线| 午夜久久久在线观看| 在线观看免费高清a一片| 国产亚洲精品久久久久久毛片 | 操美女的视频在线观看| 国产精品永久免费网站| 看片在线看免费视频| av免费在线观看网站| 悠悠久久av| 欧美av亚洲av综合av国产av| 久久中文看片网| 91精品三级在线观看| 精品久久久久久,| 国产淫语在线视频| 国产欧美日韩精品亚洲av| 极品教师在线免费播放| 久久久久国产一级毛片高清牌| 韩国av一区二区三区四区| 丝瓜视频免费看黄片| 国产亚洲欧美在线一区二区| 在线av久久热| 黄色成人免费大全| 露出奶头的视频| 香蕉国产在线看| 黄色怎么调成土黄色| 动漫黄色视频在线观看| av有码第一页| 18禁裸乳无遮挡动漫免费视频| netflix在线观看网站| 丝瓜视频免费看黄片| 久久精品人人爽人人爽视色| 国产成人av教育| 亚洲av成人一区二区三| 精品国产一区二区三区四区第35| 亚洲熟妇熟女久久| 一级作爱视频免费观看| 啦啦啦在线免费观看视频4| 每晚都被弄得嗷嗷叫到高潮| 大香蕉久久网| 精品欧美一区二区三区在线| 中文亚洲av片在线观看爽 | 淫妇啪啪啪对白视频| 欧美 亚洲 国产 日韩一| 欧美激情久久久久久爽电影 | 亚洲黑人精品在线| 99riav亚洲国产免费| 首页视频小说图片口味搜索| xxx96com| 午夜激情av网站| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲精品第一综合不卡| 波多野结衣av一区二区av| 一区福利在线观看| 精品乱码久久久久久99久播| 最近最新中文字幕大全免费视频| 国产人伦9x9x在线观看| ponron亚洲| 亚洲视频免费观看视频| 国产精品自产拍在线观看55亚洲 | 在线观看免费高清a一片| 美国免费a级毛片| 亚洲va日本ⅴa欧美va伊人久久| 9热在线视频观看99| 午夜两性在线视频| av中文乱码字幕在线| 免费看a级黄色片| 亚洲少妇的诱惑av| av视频免费观看在线观看| 欧美老熟妇乱子伦牲交| 亚洲av第一区精品v没综合| 麻豆国产av国片精品| 午夜精品在线福利| 1024香蕉在线观看| 女性生殖器流出的白浆| 午夜福利在线观看吧| 这个男人来自地球电影免费观看| 精品国产一区二区久久| bbb黄色大片| 777米奇影视久久| 窝窝影院91人妻| 天堂俺去俺来也www色官网| 91成人精品电影| 欧美午夜高清在线| 女性生殖器流出的白浆| 亚洲国产精品一区二区三区在线| 母亲3免费完整高清在线观看| 男人操女人黄网站| 十分钟在线观看高清视频www| 午夜成年电影在线免费观看| 午夜免费鲁丝| 又黄又粗又硬又大视频| 老汉色∧v一级毛片| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av | 黄色 视频免费看| 国产精品.久久久| 欧美黄色淫秽网站| 成年人午夜在线观看视频| 亚洲成人免费av在线播放| 女性生殖器流出的白浆| 69精品国产乱码久久久| 波多野结衣一区麻豆| 国产精品久久电影中文字幕 | 12—13女人毛片做爰片一| 黄色怎么调成土黄色| 超碰成人久久| 日本黄色视频三级网站网址 | 这个男人来自地球电影免费观看| 欧美日韩福利视频一区二区| 99久久综合精品五月天人人| 亚洲自偷自拍图片 自拍| 韩国精品一区二区三区| 亚洲五月天丁香| 免费观看精品视频网站| 国产又色又爽无遮挡免费看| av福利片在线| 一二三四在线观看免费中文在| 深夜精品福利| 国产亚洲欧美98| 欧美不卡视频在线免费观看 | 久久精品91无色码中文字幕| 精品高清国产在线一区| 欧美午夜高清在线| 视频区图区小说| 在线十欧美十亚洲十日本专区| 久久国产精品男人的天堂亚洲| 色精品久久人妻99蜜桃| 美国免费a级毛片| 国产高清激情床上av| 国产色视频综合| 在线观看舔阴道视频| 久久亚洲精品不卡| 美女 人体艺术 gogo| 老司机亚洲免费影院| 国产不卡av网站在线观看| 国产成人av教育| 在线十欧美十亚洲十日本专区| 精品国内亚洲2022精品成人 | 女人被狂操c到高潮| 欧美日韩国产mv在线观看视频| 一级毛片女人18水好多| 男女高潮啪啪啪动态图| 老司机亚洲免费影院| 精品视频人人做人人爽| 精品一区二区三卡| a在线观看视频网站| 免费av中文字幕在线| 又紧又爽又黄一区二区| av视频免费观看在线观看| 久久天躁狠狠躁夜夜2o2o| 正在播放国产对白刺激| 女人被躁到高潮嗷嗷叫费观| 亚洲熟妇熟女久久| 制服诱惑二区| 亚洲精品久久午夜乱码| 国产精品1区2区在线观看. | 亚洲精品粉嫩美女一区| 午夜精品久久久久久毛片777| 制服人妻中文乱码| 在线观看66精品国产| 女人被躁到高潮嗷嗷叫费观| 高清视频免费观看一区二区| 三级毛片av免费| 午夜日韩欧美国产| 在线观看一区二区三区激情| 日韩欧美免费精品| 热99re8久久精品国产| 夫妻午夜视频| 黑人欧美特级aaaaaa片| 变态另类成人亚洲欧美熟女 | 制服人妻中文乱码| 亚洲精品成人av观看孕妇| 91成年电影在线观看| 免费久久久久久久精品成人欧美视频| 妹子高潮喷水视频| 久久人妻熟女aⅴ| 黄色怎么调成土黄色| 国产又爽黄色视频| 色婷婷av一区二区三区视频| 不卡一级毛片| cao死你这个sao货| 久久精品国产综合久久久| 在线永久观看黄色视频| av片东京热男人的天堂| 精品久久久久久久久久免费视频 | 村上凉子中文字幕在线| 欧美精品一区二区免费开放| 欧美不卡视频在线免费观看 | 一区二区三区精品91| 成人精品一区二区免费| 亚洲精品一二三| 18禁裸乳无遮挡动漫免费视频| 在线观看www视频免费| 天堂中文最新版在线下载| 人妻丰满熟妇av一区二区三区 | 人人妻人人澡人人看| 精品福利观看| 狠狠婷婷综合久久久久久88av| 9色porny在线观看| 日本欧美视频一区| 夜夜爽天天搞| 777米奇影视久久| 黄色片一级片一级黄色片| 男人的好看免费观看在线视频 | 天堂俺去俺来也www色官网| tube8黄色片| 国产成人精品久久二区二区91| 久久ye,这里只有精品| www.自偷自拍.com| 51午夜福利影视在线观看| 91老司机精品| av中文乱码字幕在线| 精品国产乱码久久久久久男人| 久久草成人影院| 久久久久国产精品人妻aⅴ院 | 精品国产一区二区久久| 午夜成年电影在线免费观看| 欧美亚洲日本最大视频资源| 欧美日本中文国产一区发布| 国产欧美日韩综合在线一区二区| 在线国产一区二区在线| 99re在线观看精品视频| 久久亚洲真实| 两个人看的免费小视频| 一进一出抽搐gif免费好疼 | 亚洲欧美色中文字幕在线| 午夜亚洲福利在线播放| 99热网站在线观看| 9191精品国产免费久久| 国产精品免费一区二区三区在线 | 欧美性长视频在线观看| 岛国在线观看网站| 国产1区2区3区精品| 老熟妇仑乱视频hdxx| 国产蜜桃级精品一区二区三区 | 亚洲欧洲精品一区二区精品久久久| 老熟妇乱子伦视频在线观看| 巨乳人妻的诱惑在线观看| 久久人妻熟女aⅴ| 高清av免费在线| 亚洲 国产 在线| 交换朋友夫妻互换小说| 黑人欧美特级aaaaaa片| 十八禁高潮呻吟视频| 成年人黄色毛片网站| 黄色成人免费大全| 久久草成人影院| 欧美在线一区亚洲| 高清av免费在线| 亚洲熟女毛片儿| 国产免费现黄频在线看| 最新在线观看一区二区三区| 成年版毛片免费区| 日韩精品免费视频一区二区三区| 精品国产一区二区三区久久久樱花| 精品久久久久久电影网| 久久精品91无色码中文字幕| 国产麻豆69| 亚洲精品在线美女| 欧美日韩成人在线一区二区| 亚洲五月色婷婷综合| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 国产不卡av网站在线观看| 成熟少妇高潮喷水视频| 亚洲少妇的诱惑av| 99久久精品国产亚洲精品| 午夜视频精品福利| 如日韩欧美国产精品一区二区三区| 一二三四在线观看免费中文在| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 在线天堂中文资源库| 久久人妻av系列| 99re在线观看精品视频| 色综合欧美亚洲国产小说| 国产在线一区二区三区精| 69精品国产乱码久久久| 欧美 日韩 精品 国产| 久久精品亚洲av国产电影网| 国产亚洲一区二区精品| 国产色视频综合| 亚洲精品成人av观看孕妇| 亚洲aⅴ乱码一区二区在线播放 | 亚洲欧美激情在线| 日韩人妻精品一区2区三区| 99国产精品一区二区三区| 99精国产麻豆久久婷婷| 国产亚洲精品一区二区www | 又紧又爽又黄一区二区| 视频在线观看一区二区三区| 麻豆成人av在线观看| 亚洲av熟女| 天天躁日日躁夜夜躁夜夜| 久久久久久久午夜电影 | 欧美日韩黄片免| 久久精品成人免费网站| 少妇 在线观看| 男女下面插进去视频免费观看| 天天躁日日躁夜夜躁夜夜| 人人妻人人澡人人爽人人夜夜| 精品少妇久久久久久888优播| 亚洲成人国产一区在线观看| 国产片内射在线| 久久狼人影院| 亚洲精品国产一区二区精华液| 国产一区有黄有色的免费视频| 很黄的视频免费| 嫩草影视91久久| 一本大道久久a久久精品| 女性生殖器流出的白浆| 王馨瑶露胸无遮挡在线观看| 日韩熟女老妇一区二区性免费视频| xxxhd国产人妻xxx| 亚洲精品国产精品久久久不卡| 黑人欧美特级aaaaaa片| 亚洲熟妇中文字幕五十中出 | 国产无遮挡羞羞视频在线观看| 99久久精品国产亚洲精品| 欧美乱色亚洲激情| 欧美日韩亚洲综合一区二区三区_| 一级毛片高清免费大全| 黑丝袜美女国产一区| 9191精品国产免费久久| 国产精品av久久久久免费| 在线观看舔阴道视频| 大陆偷拍与自拍| av电影中文网址| 老司机影院毛片| 久久 成人 亚洲| 成年人免费黄色播放视频| 午夜免费鲁丝| 国产精品.久久久| 99国产精品99久久久久| 视频区欧美日本亚洲| 91精品国产国语对白视频| 在线天堂中文资源库| 一级毛片精品| 亚洲av日韩精品久久久久久密| 日日摸夜夜添夜夜添小说| 女性生殖器流出的白浆| 女人爽到高潮嗷嗷叫在线视频| 久久中文字幕一级| 又黄又爽又免费观看的视频| 亚洲av欧美aⅴ国产| 国产精品1区2区在线观看. | 美女国产高潮福利片在线看| 久久精品91无色码中文字幕| 亚洲精品成人av观看孕妇| a在线观看视频网站| 一区二区三区激情视频| 久久九九热精品免费| 久久久精品区二区三区| 欧美日韩福利视频一区二区| av线在线观看网站| 夜夜夜夜夜久久久久| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| av在线播放免费不卡| 人人妻人人澡人人看| 国产亚洲一区二区精品| 亚洲精品成人av观看孕妇| 国产野战对白在线观看| 超色免费av| 欧美日韩国产mv在线观看视频| 成年动漫av网址| 免费观看a级毛片全部| 成年人午夜在线观看视频| 成年人免费黄色播放视频| 狂野欧美激情性xxxx| 欧美日韩黄片免| 1024视频免费在线观看| 日日夜夜操网爽| 极品少妇高潮喷水抽搐| 成人亚洲精品一区在线观看| 老司机午夜十八禁免费视频| av不卡在线播放| 欧美黑人欧美精品刺激| 18禁裸乳无遮挡动漫免费视频| 中文亚洲av片在线观看爽 | 国产高清国产精品国产三级| av天堂在线播放| 久久精品亚洲av国产电影网| 757午夜福利合集在线观看| 视频区欧美日本亚洲| 两个人看的免费小视频| 亚洲精品久久成人aⅴ小说| 欧美丝袜亚洲另类 | 热99久久久久精品小说推荐| 成人黄色视频免费在线看| 九色亚洲精品在线播放| 日韩欧美三级三区| 99热国产这里只有精品6| av天堂在线播放| 大码成人一级视频| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 国产三级黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一区有黄有色的免费视频| 大片电影免费在线观看免费| 日韩欧美三级三区| 亚洲成人手机| 亚洲欧美一区二区三区黑人| 国产精品亚洲一级av第二区| 操出白浆在线播放| 在线视频色国产色| 成人18禁在线播放| 国产午夜精品久久久久久| 搡老熟女国产l中国老女人| 久久99一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩成人在线一区二区| 日韩制服丝袜自拍偷拍| 亚洲综合色网址| 国产成人系列免费观看| 亚洲欧美激情综合另类| 少妇粗大呻吟视频| 美女高潮喷水抽搐中文字幕| 亚洲人成电影观看| a级毛片黄视频| 亚洲黑人精品在线| 韩国av一区二区三区四区| 亚洲精品一二三| 成人18禁高潮啪啪吃奶动态图| 一区二区日韩欧美中文字幕| 美女扒开内裤让男人捅视频| 国产又爽黄色视频| 亚洲精品乱久久久久久| 嫩草影视91久久| 超色免费av| 国产精品1区2区在线观看. | 国产熟女午夜一区二区三区| 亚洲精品成人av观看孕妇| 最近最新免费中文字幕在线| 国产免费男女视频| 免费高清在线观看日韩| 免费在线观看日本一区| 国产一区二区三区综合在线观看| 丰满迷人的少妇在线观看| 国产不卡av网站在线观看| av一本久久久久| 另类亚洲欧美激情| 成熟少妇高潮喷水视频| av欧美777| 亚洲国产中文字幕在线视频| 日韩欧美一区视频在线观看| 黑人猛操日本美女一级片| 精品视频人人做人人爽| √禁漫天堂资源中文www| 亚洲片人在线观看| 黄色丝袜av网址大全| 亚洲片人在线观看| 久久久久久人人人人人| 免费人成视频x8x8入口观看| 免费在线观看日本一区| 美女午夜性视频免费| 免费高清在线观看日韩| 19禁男女啪啪无遮挡网站| 国产单亲对白刺激| 精品一区二区三区视频在线观看免费 | 人人妻人人澡人人爽人人夜夜| 校园春色视频在线观看| 亚洲精品中文字幕在线视频| 亚洲av日韩精品久久久久久密| 可以免费在线观看a视频的电影网站| 国产亚洲av高清不卡| 视频在线观看一区二区三区| 久久香蕉精品热| 青草久久国产| 12—13女人毛片做爰片一| 精品熟女少妇八av免费久了| 国产亚洲精品第一综合不卡| 两个人看的免费小视频| 狠狠婷婷综合久久久久久88av| 亚洲中文字幕日韩| 成人国语在线视频| 一本一本久久a久久精品综合妖精| 黑人操中国人逼视频| 久久精品91无色码中文字幕| 亚洲精品在线观看二区| 中文字幕最新亚洲高清| 90打野战视频偷拍视频| 在线观看免费视频日本深夜| 国产成人免费观看mmmm| 一本一本久久a久久精品综合妖精| 国产成+人综合+亚洲专区| 中文欧美无线码| 亚洲欧美色中文字幕在线| 午夜免费鲁丝| 99国产精品一区二区蜜桃av | 久久这里只有精品19| 热99久久久久精品小说推荐| 亚洲欧美一区二区三区久久| 国产精品偷伦视频观看了| 韩国精品一区二区三区| 搡老乐熟女国产| 精品久久久久久久毛片微露脸| 又紧又爽又黄一区二区| 久久久久久亚洲精品国产蜜桃av| 深夜精品福利| 欧美激情久久久久久爽电影 | 黄片大片在线免费观看| 男女午夜视频在线观看| 黄色a级毛片大全视频| 黄频高清免费视频| 少妇粗大呻吟视频| 法律面前人人平等表现在哪些方面| 国产免费现黄频在线看| 少妇 在线观看| 久久香蕉激情| 男女午夜视频在线观看| 精品久久久久久电影网| xxxhd国产人妻xxx| 美女高潮到喷水免费观看| 精品久久久精品久久久| 久久亚洲真实| 国产精品久久久久成人av| 国产精品电影一区二区三区 | 国产一区二区三区在线臀色熟女 | 亚洲欧美激情在线| 69精品国产乱码久久久| 婷婷精品国产亚洲av在线 | 在线观看一区二区三区激情| 夫妻午夜视频| 亚洲熟妇熟女久久| av电影中文网址| 丝袜人妻中文字幕| 国产高清视频在线播放一区| 可以免费在线观看a视频的电影网站| 久久久精品免费免费高清| 国产精品电影一区二区三区 | 无限看片的www在线观看| 亚洲av欧美aⅴ国产| 久久婷婷成人综合色麻豆| 夜夜夜夜夜久久久久| 国产高清激情床上av| 国产99久久九九免费精品| 天堂√8在线中文| 又黄又粗又硬又大视频| 别揉我奶头~嗯~啊~动态视频| 99久久99久久久精品蜜桃| 亚洲熟妇熟女久久| 精品少妇久久久久久888优播| 中文字幕人妻丝袜一区二区| 最近最新中文字幕大全电影3 | 免费观看精品视频网站| 国产精品国产av在线观看| 精品久久蜜臀av无|