• <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
    国产白丝娇喘喷水9色精品| 亚洲av欧美aⅴ国产| 国产精品 国内视频| 午夜av观看不卡| 少妇的丰满在线观看| 国产激情久久老熟女| 亚洲国产精品一区三区| 丝袜美足系列| 九色亚洲精品在线播放| 中文字幕亚洲精品专区| 丁香六月天网| 秋霞伦理黄片| 国产在线视频一区二区| 亚洲第一区二区三区不卡| 免费观看在线日韩| 免费人妻精品一区二区三区视频| 亚洲美女黄色视频免费看| 色婷婷av一区二区三区视频| av免费在线看不卡| 亚洲av电影在线观看一区二区三区| 久久久久久久精品精品| 亚洲成人手机| 97人妻天天添夜夜摸| 亚洲精品国产av蜜桃| 夜夜爽夜夜爽视频| 精品少妇黑人巨大在线播放| 精品国产露脸久久av麻豆| 日韩人妻精品一区2区三区| 最黄视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 美女xxoo啪啪120秒动态图| 国产精品久久久久久av不卡| 下体分泌物呈黄色| 三上悠亚av全集在线观看| 考比视频在线观看| 黄色配什么色好看| 国产精品.久久久| 丝袜脚勾引网站| 九色亚洲精品在线播放| 成年av动漫网址| 成年女人在线观看亚洲视频| 女性生殖器流出的白浆| 美女大奶头黄色视频| 免费在线观看黄色视频的| 一边亲一边摸免费视频| 黄色毛片三级朝国网站| 这个男人来自地球电影免费观看 | 水蜜桃什么品种好| 久久久久久久国产电影| 两个人免费观看高清视频| 少妇高潮的动态图| 男女边摸边吃奶| 欧美日韩视频精品一区| 亚洲国产欧美在线一区| 人人妻人人添人人爽欧美一区卜| 亚洲欧美清纯卡通| 久久精品久久久久久久性| 少妇精品久久久久久久| 伦精品一区二区三区| 午夜福利在线观看免费完整高清在| 香蕉国产在线看| 高清av免费在线| 国语对白做爰xxxⅹ性视频网站| 99热这里只有是精品在线观看| 国产日韩欧美亚洲二区| 欧美人与善性xxx| 大话2 男鬼变身卡| 哪个播放器可以免费观看大片| 一级毛片 在线播放| 一边亲一边摸免费视频| 欧美少妇被猛烈插入视频| 久久午夜福利片| 国产精品一二三区在线看| 日本欧美视频一区| 少妇高潮的动态图| 久久免费观看电影| 亚洲国产av影院在线观看| 9色porny在线观看| 五月天丁香电影| 精品99又大又爽又粗少妇毛片| 精品一品国产午夜福利视频| 少妇高潮的动态图| 久久99精品国语久久久| 国产精品人妻久久久影院| 免费观看性生交大片5| 又黄又粗又硬又大视频| 国产精品久久久久久久久免| 午夜福利乱码中文字幕| 亚洲经典国产精华液单| 国产精品偷伦视频观看了| 国产男人的电影天堂91| 高清不卡的av网站| 天堂俺去俺来也www色官网| 精品国产一区二区三区四区第35| 亚洲成色77777| 男人操女人黄网站| 日本猛色少妇xxxxx猛交久久| 精品少妇久久久久久888优播| 亚洲欧美中文字幕日韩二区| 97在线人人人人妻| 亚洲精品美女久久久久99蜜臀 | xxx大片免费视频| 国产成人免费无遮挡视频| 人妻系列 视频| 欧美97在线视频| 欧美日韩综合久久久久久| 欧美人与性动交α欧美精品济南到 | 精品99又大又爽又粗少妇毛片| 一个人免费看片子| 日韩欧美一区视频在线观看| 亚洲人成77777在线视频| 中文欧美无线码| 国产亚洲欧美精品永久| 嫩草影院入口| 国产无遮挡羞羞视频在线观看| 国产淫语在线视频| 免费观看在线日韩| 少妇人妻 视频| 久久人人爽av亚洲精品天堂| 中文乱码字字幕精品一区二区三区| 自线自在国产av| 捣出白浆h1v1| freevideosex欧美| 欧美精品国产亚洲| 一区二区三区乱码不卡18| 亚洲精品456在线播放app| 蜜桃在线观看..| 成人毛片60女人毛片免费| 最近最新中文字幕免费大全7| 国产成人a∨麻豆精品| 免费看av在线观看网站| 另类亚洲欧美激情| 亚洲国产精品一区二区三区在线| 精品久久久久久电影网| 色视频在线一区二区三区| 51国产日韩欧美| 中文字幕精品免费在线观看视频 | 最黄视频免费看| 天堂8中文在线网| 免费黄色在线免费观看| 日本欧美视频一区| 在线天堂最新版资源| 日韩伦理黄色片| 色视频在线一区二区三区| 少妇熟女欧美另类| 成人午夜精彩视频在线观看| 久久99一区二区三区| 成人国语在线视频| 99久久综合免费| 国产亚洲一区二区精品| 51国产日韩欧美| 亚洲国产看品久久| 在线天堂最新版资源| 国产一区二区在线观看av| 亚洲精品久久午夜乱码| 亚洲欧美一区二区三区国产| 夜夜爽夜夜爽视频| 亚洲欧美一区二区三区黑人 | 制服人妻中文乱码| 午夜激情av网站| 久久这里有精品视频免费| 国产片内射在线| 人妻人人澡人人爽人人| 亚洲国产精品国产精品| 精品酒店卫生间| 免费观看a级毛片全部| 亚洲国产精品一区三区| 免费黄色在线免费观看| 人妻 亚洲 视频| 人成视频在线观看免费观看| 亚洲综合色惰| 国产精品久久久久久av不卡| 女人被躁到高潮嗷嗷叫费观| 不卡视频在线观看欧美| 只有这里有精品99| 美女xxoo啪啪120秒动态图| 青春草国产在线视频| 80岁老熟妇乱子伦牲交| av网站免费在线观看视频| 少妇的逼好多水| 色婷婷久久久亚洲欧美| 国产淫语在线视频| 伊人亚洲综合成人网| 国产男女内射视频| 五月天丁香电影| 精品一区二区三卡| 男女下面插进去视频免费观看 | 在线观看一区二区三区激情| 国产精品久久久久久久电影| 精品一区二区三区四区五区乱码 | 久久精品aⅴ一区二区三区四区 | 全区人妻精品视频| 菩萨蛮人人尽说江南好唐韦庄| 国产免费视频播放在线视频| 免费黄网站久久成人精品| 精品酒店卫生间| 搡女人真爽免费视频火全软件| 国产精品久久久久久av不卡| 韩国精品一区二区三区 | 精品卡一卡二卡四卡免费| 熟妇人妻不卡中文字幕| 国产白丝娇喘喷水9色精品| 亚洲欧美成人综合另类久久久| 久久 成人 亚洲| 免费黄网站久久成人精品| 久久久国产精品麻豆| 国产精品秋霞免费鲁丝片| 美女福利国产在线| 亚洲少妇的诱惑av| 在线天堂最新版资源| 久久久久视频综合| 国产免费一区二区三区四区乱码| 黄色视频在线播放观看不卡| 久久久久久久国产电影| 80岁老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 久久ye,这里只有精品| 三级国产精品片| 欧美激情国产日韩精品一区| 日本色播在线视频| 啦啦啦视频在线资源免费观看| 国产爽快片一区二区三区| 久久女婷五月综合色啪小说| 亚洲国产精品专区欧美| 国语对白做爰xxxⅹ性视频网站| 午夜福利视频在线观看免费| 亚洲欧美一区二区三区黑人 | 日日摸夜夜添夜夜爱| 久久精品国产亚洲av天美| 女性生殖器流出的白浆| 亚洲综合色惰| 免费日韩欧美在线观看| 精品一品国产午夜福利视频| 最新的欧美精品一区二区| 国产日韩一区二区三区精品不卡| 99热6这里只有精品| 久久国产精品男人的天堂亚洲 | 午夜91福利影院| 国产日韩欧美在线精品| 久久久久久久久久久久大奶| 男女高潮啪啪啪动态图| 建设人人有责人人尽责人人享有的| 99国产精品免费福利视频| 26uuu在线亚洲综合色| 久久精品国产综合久久久 | 啦啦啦在线观看免费高清www| 国产精品久久久久久精品电影小说| av视频免费观看在线观看| 一区二区三区精品91| 日本av手机在线免费观看| 欧美成人午夜免费资源| 亚洲第一区二区三区不卡| 夫妻午夜视频| 国产精品一二三区在线看| 国产在线免费精品| 99久久中文字幕三级久久日本| 五月天丁香电影| 18+在线观看网站| 一区在线观看完整版| 国产免费一区二区三区四区乱码| 亚洲欧美色中文字幕在线| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图| 国产成人欧美| 国产高清不卡午夜福利| 777米奇影视久久| 亚洲欧美一区二区三区国产| 欧美+日韩+精品| 亚洲国产av影院在线观看| 成年人午夜在线观看视频| 久久亚洲国产成人精品v| 国产成人精品在线电影| 黑人高潮一二区| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| 日本wwww免费看| 国产精品嫩草影院av在线观看| 久久这里有精品视频免费| 日本av免费视频播放| 亚洲 欧美一区二区三区| 久久精品国产鲁丝片午夜精品| 亚洲一码二码三码区别大吗| 91久久精品国产一区二区三区| 国产成人欧美| 你懂的网址亚洲精品在线观看| 国产精品欧美亚洲77777| 国产精品久久久久久精品古装| 久久久久久久久久人人人人人人| 国产成人精品婷婷| 国产精品一二三区在线看| videosex国产| 2018国产大陆天天弄谢| 少妇熟女欧美另类| 肉色欧美久久久久久久蜜桃| 精品久久国产蜜桃| 欧美 日韩 精品 国产| 一级毛片黄色毛片免费观看视频| 九色成人免费人妻av| 少妇熟女欧美另类| 五月开心婷婷网| 啦啦啦在线观看免费高清www| 国产乱来视频区| 久久精品人人爽人人爽视色| av又黄又爽大尺度在线免费看| 亚洲欧美日韩另类电影网站| 亚洲国产日韩一区二区| 日韩中文字幕视频在线看片| 少妇 在线观看| 欧美激情极品国产一区二区三区 | 久久久久久久亚洲中文字幕| 日产精品乱码卡一卡2卡三| 男男h啪啪无遮挡| 乱人伦中国视频| 亚洲成人av在线免费| 久久这里只有精品19| 丰满少妇做爰视频| h视频一区二区三区| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 两个人免费观看高清视频| 精品国产乱码久久久久久小说| 热99国产精品久久久久久7| 熟女av电影| 日韩制服骚丝袜av| 最近的中文字幕免费完整| 成年人免费黄色播放视频| 免费不卡的大黄色大毛片视频在线观看| 一边摸一边做爽爽视频免费| 女性生殖器流出的白浆| 国产亚洲午夜精品一区二区久久| 性色av一级| 全区人妻精品视频| 欧美精品av麻豆av| 欧美成人精品欧美一级黄| 成年美女黄网站色视频大全免费| 国产成人av激情在线播放| 国产高清三级在线| 欧美精品一区二区免费开放| 99久久精品国产国产毛片| 99re6热这里在线精品视频| 国产精品人妻久久久影院| 另类精品久久| 天堂中文最新版在线下载| 午夜免费鲁丝| 伊人久久国产一区二区| 亚洲精品自拍成人| 精品国产国语对白av| h视频一区二区三区| 亚洲欧洲国产日韩| 最近的中文字幕免费完整| 久久久久久久久久人人人人人人| 精品一区在线观看国产| 亚洲伊人色综图| 18在线观看网站| av国产久精品久网站免费入址| 国产精品人妻久久久久久| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 亚洲一码二码三码区别大吗| 日日爽夜夜爽网站| 成年女人在线观看亚洲视频| 日韩不卡一区二区三区视频在线| 亚洲五月色婷婷综合| 亚洲成人av在线免费| 欧美性感艳星| 午夜影院在线不卡| 99香蕉大伊视频| 熟女人妻精品中文字幕| 国产免费现黄频在线看| 激情五月婷婷亚洲| 久久久a久久爽久久v久久| 午夜av观看不卡| 精品一区二区三区视频在线| 性高湖久久久久久久久免费观看| 精品一区在线观看国产| 乱码一卡2卡4卡精品| 街头女战士在线观看网站| 国产精品嫩草影院av在线观看| 国产成人精品福利久久| 欧美少妇被猛烈插入视频| 乱码一卡2卡4卡精品| 日本色播在线视频| 黄色视频在线播放观看不卡| 国产男女超爽视频在线观看| 久久久久人妻精品一区果冻| 亚洲国产毛片av蜜桃av| 欧美人与性动交α欧美精品济南到 | 99香蕉大伊视频| 少妇 在线观看| 免费观看在线日韩| 91精品伊人久久大香线蕉| 国产伦理片在线播放av一区| 亚洲综合色网址| 国产爽快片一区二区三区| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频 | 91精品伊人久久大香线蕉| 日韩成人伦理影院| 啦啦啦在线观看免费高清www| 久久精品国产a三级三级三级| 国产精品一国产av| 9热在线视频观看99| 在线天堂最新版资源| 亚洲精品av麻豆狂野| 97超碰精品成人国产| 大片电影免费在线观看免费| 精品亚洲乱码少妇综合久久| 国产一区二区三区av在线| av女优亚洲男人天堂| 国产av码专区亚洲av| 久久97久久精品| 国产精品不卡视频一区二区| 久久久国产一区二区| av线在线观看网站| 亚洲精品久久久久久婷婷小说| 人妻人人澡人人爽人人| 国产欧美亚洲国产| 久久影院123| 精品视频人人做人人爽| 99久久中文字幕三级久久日本| 韩国精品一区二区三区 | 麻豆乱淫一区二区| 免费大片黄手机在线观看| 欧美精品国产亚洲| 九九爱精品视频在线观看| 国产成人精品一,二区| 肉色欧美久久久久久久蜜桃| 日日摸夜夜添夜夜爱| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人毛片a级毛片在线播放| 午夜老司机福利剧场| av有码第一页| 视频中文字幕在线观看| 成年女人在线观看亚洲视频| 99久久中文字幕三级久久日本| 韩国精品一区二区三区 | videossex国产| a级片在线免费高清观看视频| 日韩三级伦理在线观看| 日韩 亚洲 欧美在线| 岛国毛片在线播放| 色哟哟·www| 国产在线一区二区三区精| 两性夫妻黄色片 | 久久久精品免费免费高清| 一区二区av电影网| 三级国产精品片| 久久人妻熟女aⅴ| 91在线精品国自产拍蜜月| 成人午夜精彩视频在线观看| 亚洲激情五月婷婷啪啪| 99久久精品国产国产毛片| 人体艺术视频欧美日本| 午夜91福利影院| 全区人妻精品视频| 色5月婷婷丁香| 国产精品女同一区二区软件| 久久亚洲国产成人精品v| 精品第一国产精品| 亚洲,欧美,日韩| 久久热在线av| 蜜桃在线观看..| 国产高清国产精品国产三级| 国产午夜精品一二区理论片| 男女国产视频网站| a级片在线免费高清观看视频| 搡老乐熟女国产| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 久久久久久久久久成人| www.熟女人妻精品国产 | 中文字幕人妻丝袜制服| 男女国产视频网站| 免费黄网站久久成人精品| 免费在线观看黄色视频的| 韩国高清视频一区二区三区| 国产极品天堂在线| 一边摸一边做爽爽视频免费| 亚洲av中文av极速乱| 一二三四在线观看免费中文在 | 精品第一国产精品| 欧美激情 高清一区二区三区| 国产亚洲最大av| 午夜免费男女啪啪视频观看| 国产精品三级大全| 秋霞伦理黄片| 欧美国产精品一级二级三级| 天天影视国产精品| 亚洲综合色网址| 午夜免费观看性视频| 久久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 国产亚洲欧美精品永久| 在线观看免费视频网站a站| 又粗又硬又长又爽又黄的视频| 成人综合一区亚洲| 九色成人免费人妻av| 精品人妻在线不人妻| 色哟哟·www| 少妇的逼水好多| 欧美国产精品一级二级三级| 国产一卡二卡三卡精品| 在线播放国产精品三级| 宅男免费午夜| 91老司机精品| 日韩三级视频一区二区三区| 久久青草综合色| 岛国毛片在线播放| 看黄色毛片网站| 欧美日韩亚洲综合一区二区三区_| 老司机深夜福利视频在线观看| 国产成人精品久久二区二区91| 成年版毛片免费区| 视频在线观看一区二区三区| 不卡av一区二区三区| 国产精品自产拍在线观看55亚洲 | 热99久久久久精品小说推荐| 丰满人妻熟妇乱又伦精品不卡| 国产91精品成人一区二区三区| 国产在线精品亚洲第一网站| 亚洲成av片中文字幕在线观看| 精品一区二区三区av网在线观看| 成年人免费黄色播放视频| 热re99久久国产66热| 久久中文字幕人妻熟女| 不卡av一区二区三区| 欧美精品人与动牲交sv欧美| 在线观看免费视频日本深夜| 国产色视频综合| 午夜精品国产一区二区电影| 无人区码免费观看不卡| 精品亚洲成国产av| 精品久久蜜臀av无| 国产av一区二区精品久久| 久久热在线av| 亚洲精品粉嫩美女一区| 亚洲av成人一区二区三| 久久中文字幕一级| 美女国产高潮福利片在线看| 1024视频免费在线观看| 久久久精品区二区三区| 欧美日韩视频精品一区| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 亚洲成人免费av在线播放| 中出人妻视频一区二区| 黄色女人牲交| 精品少妇久久久久久888优播| 国产亚洲欧美在线一区二区| 9191精品国产免费久久| 大型av网站在线播放| 麻豆乱淫一区二区| 在线观看日韩欧美| 精品亚洲成a人片在线观看| 欧美另类亚洲清纯唯美| 国产蜜桃级精品一区二区三区 | 宅男免费午夜| 99精品在免费线老司机午夜| 色94色欧美一区二区| 亚洲自偷自拍图片 自拍| 亚洲av欧美aⅴ国产| 亚洲精品国产区一区二| 久久香蕉激情| 午夜福利视频在线观看免费| 午夜成年电影在线免费观看| 欧美黑人欧美精品刺激| 精品国产国语对白av| cao死你这个sao货| 国产99白浆流出| 精品久久久久久电影网| 一a级毛片在线观看| 国产在线精品亚洲第一网站| videosex国产| 女性被躁到高潮视频| 91av网站免费观看| 又紧又爽又黄一区二区| 视频在线观看一区二区三区| 激情视频va一区二区三区| 国产精华一区二区三区| 亚洲精品在线美女| 欧美成狂野欧美在线观看| 欧美 日韩 精品 国产| 久久这里只有精品19| 亚洲国产精品sss在线观看 | 精品熟女少妇八av免费久了| 91成人精品电影| 夫妻午夜视频| 日韩三级视频一区二区三区| 国产免费男女视频| 一级,二级,三级黄色视频| 身体一侧抽搐| 很黄的视频免费| 色在线成人网| 日韩欧美在线二视频 | 一级毛片精品| 日本五十路高清| 国产精品久久电影中文字幕 | 久久久精品国产亚洲av高清涩受| 啦啦啦免费观看视频1| 精品久久久久久久久久免费视频 | 亚洲熟妇熟女久久| 欧美乱色亚洲激情| 欧美日韩成人在线一区二区| 嫩草影视91久久| 欧美中文综合在线视频| 国产欧美亚洲国产| 757午夜福利合集在线观看| 成年人黄色毛片网站| 午夜福利视频在线观看免费| 亚洲精品国产色婷婷电影| 丰满的人妻完整版| 亚洲欧洲精品一区二区精品久久久| 91精品三级在线观看| 欧美人与性动交α欧美软件| 久久久久久久精品吃奶|