基于ANSYS二次開發(fā)冷卻塔施工全過程風(fēng)致極限承載性能研究

杜凌云, 柯世堂

(南京航空航天大學(xué) 土木工程系，南京　210016)

?

基于ANSYS二次開發(fā)冷卻塔施工全過程風(fēng)致極限承載性能研究

杜凌云, 柯世堂

(南京航空航天大學(xué) 土木工程系，南京210016)

為研究施工全過程中冷卻塔風(fēng)致強(qiáng)度和穩(wěn)定極限承載性能，首先基于UIDL和APDL對(duì)ANSYS進(jìn)行相應(yīng)功能的界面化模塊二次開發(fā)，在此基礎(chǔ)上進(jìn)行智能化有限元建模，通過分別施加規(guī)范與風(fēng)洞試驗(yàn)風(fēng)荷載，進(jìn)行施工全過程中冷卻塔強(qiáng)度及屈曲失穩(wěn)極限承載能力計(jì)算，且考慮對(duì)比了混凝土齡期變化對(duì)屈曲穩(wěn)定性的影響。結(jié)果表明：臨界風(fēng)速在冷卻塔施工初始階段迅速降低后漸趨平穩(wěn)，整體風(fēng)致強(qiáng)度極限承載能力要顯著大于屈曲失穩(wěn)極限承載能力，計(jì)入齡期變化后冷卻塔屈曲模態(tài)與不計(jì)時(shí)相似，但臨界風(fēng)速相對(duì)減??；對(duì)比各施工階段不同工況下的兩種承載能力發(fā)現(xiàn)最小臨界風(fēng)速為100.7 m/s，遠(yuǎn)大于當(dāng)?shù)?0年一遇最大設(shè)計(jì)風(fēng)速29.7 m/s。算例分析表明本文二次開發(fā)的冷卻塔結(jié)構(gòu)智能化分析模塊實(shí)用高效，實(shí)現(xiàn)了冷卻塔施工全過程極限承載能力有限元分析的參數(shù)化和可視化。

大型冷卻塔;二次開發(fā);施工全過程;屈曲失穩(wěn);極限承載性能

大型雙曲冷卻塔[1-3]是火/核電廠重要構(gòu)筑物之一，屬于典型的高聳薄壁結(jié)構(gòu)，其自振頻率低且分布密集，模態(tài)耦合性強(qiáng)，風(fēng)荷載是結(jié)構(gòu)設(shè)計(jì)的主要控制荷載之一[4-5]，在風(fēng)荷載作用下塔筒的極限承載力和穩(wěn)定性能研究具有重要的工程價(jià)值。國(guó)內(nèi)外學(xué)者很早就開展了對(duì)運(yùn)行中冷卻塔穩(wěn)定性及極限承載力的研究工作，文獻(xiàn)[6]采用Abaqus軟件分析了大型雙曲冷卻塔結(jié)構(gòu)的強(qiáng)度和穩(wěn)定性；文獻(xiàn)[7]研究了初始缺陷對(duì)冷卻塔結(jié)構(gòu)極限承載力的影響；文獻(xiàn)[8]探討了開孔對(duì)排煙冷卻塔強(qiáng)度和穩(wěn)定性的影響；文獻(xiàn)[9]借助Staad有限元軟件對(duì)施工過程中冷卻塔的整體穩(wěn)定性能進(jìn)行了分析。已有研究成果為大型冷卻塔的結(jié)構(gòu)設(shè)計(jì)提供了一定的參考依據(jù)，但針對(duì)施工全過程中冷卻塔的穩(wěn)定和強(qiáng)度極限承載性能進(jìn)行系統(tǒng)分析的工作較少，考慮到施工過程中混凝土的強(qiáng)度和彈性模量隨時(shí)間逐步增長(zhǎng)，隨施工高度及淋水面積的逐漸增大，在自重及風(fēng)載作用下通風(fēng)筒的極限承載能力是否滿足要求值得進(jìn)行深入研究，以控制冷卻塔的施工進(jìn)度[10]。

鑒于此，本文首先基于UIDL和APDL對(duì)ANSYS有限元軟件進(jìn)行中文界面智能化二次開發(fā)，以某大型雙曲冷卻塔為例，采用自行研發(fā)的冷卻塔結(jié)構(gòu)分析模塊建立參數(shù)化實(shí)體模型，針對(duì)多種組合工況在不同施工高度下對(duì)冷卻塔施加風(fēng)荷載，分別采用逐級(jí)加載法及屈曲系數(shù)法對(duì)施工全過程中的強(qiáng)度極限承載能力和屈曲失穩(wěn)極限承載能力進(jìn)行分析，并考慮對(duì)比了混凝土齡期變化對(duì)屈曲穩(wěn)定性的影響。最終綜合比選了不同工況下冷卻塔施工全過程的臨界風(fēng)速，為施工狀態(tài)下冷卻塔的抗風(fēng)穩(wěn)定和安全性能提供參考依據(jù)。

1　ANSYS二次開發(fā)

ANSYS作為一種大型通用有限元分析軟件，在界面中直接建立復(fù)雜的幾何模型時(shí)操作繁瑣、建模效率低且不易修改和掌握，但由于ANSYS具有良好的開放性和可定制性，為多種語言的二次開發(fā)提供接口。其中，APDL作為ANSYS的專用解釋型語言之一，用來完成有限元常規(guī)分析操作或通過參數(shù)化變量方式建立分析模型[11]，可實(shí)現(xiàn)智能化的有限元分析過程，滿足用戶對(duì)重復(fù)性設(shè)計(jì)分析的需要，其分析功能強(qiáng)大，但文件方式不直觀、可視化程度差，僅能開發(fā)出簡(jiǎn)單的界面。UIDL是ANSYS專門為用戶提供編寫或改造圖形界面的設(shè)計(jì)語言，允許用戶靈活改變GUI的組件，實(shí)現(xiàn)自己的個(gè)性菜單定制，進(jìn)而完成主菜單系統(tǒng)及菜單項(xiàng)、對(duì)話框及聯(lián)機(jī)幫助系統(tǒng)三種圖形界面的設(shè)計(jì)，是進(jìn)行ANSYS二次開發(fā)強(qiáng)有力的工具[12]，其控制GUI圖形界面的開發(fā)，但不能實(shí)現(xiàn)復(fù)雜問題的建模和分析。故用戶可結(jié)合自身的需要，基于APDL和UIDL兩種設(shè)計(jì)語言針對(duì)特定專業(yè)問題在標(biāo)準(zhǔn)ANSYS版本上進(jìn)行修改編制和功能補(bǔ)充，開發(fā)出具有行業(yè)特點(diǎn)、操作簡(jiǎn)單、易用高效、界面友好的專用有限元程序模塊。

鑒于此，本文結(jié)合APDL和UIDL語言進(jìn)行ANSYS二次開發(fā)，并應(yīng)用于大型雙曲冷卻塔有限元結(jié)構(gòu)分析，在UIDL開發(fā)過程中嵌入APDL程序，對(duì)復(fù)雜的建模、計(jì)算、處理過程進(jìn)行封裝，將有限元建模過程及結(jié)構(gòu)抗風(fēng)分析等功能加入主菜單，開發(fā)出正確、實(shí)用且漢化版的冷卻塔程序模塊，以便于用戶和程序的交互，使非專業(yè)的有限元分析人員也可通過輸入簡(jiǎn)單的參數(shù)，實(shí)現(xiàn)大型雙曲冷卻塔參數(shù)化、可視化的建模及分析。二次開發(fā)主要步驟如下：

(1) 建立MYMENU.GRN菜單控制文件和MYFUNC.GRN對(duì)話框控制文件；

(2) 在menulist120.ans添加“當(dāng)前工作路徑MYMENU.GRN”和“當(dāng)前工作路徑MYFUNC.GRN”指向命令；

(3) 在控制文件構(gòu)造塊尾部前添加Cmd_Filename命令即可實(shí)現(xiàn)對(duì)APDL程序的調(diào)用，將菜單控制文件保存為ASIN格式，對(duì)話框控制文件保存為UTF-8格式，可實(shí)現(xiàn)菜單與對(duì)話框的漢化；

(4) 將三者放置于當(dāng)前工作目錄下，啟動(dòng)ANSYS自動(dòng)調(diào)用安裝目錄的uidl文件夾中的menulist120.ans文件，并據(jù)此運(yùn)行其指向的UIDL控制文件，實(shí)現(xiàn)自定義的GUI界面，主菜單顯示及部分參數(shù)對(duì)話框如圖1和2所示。

圖1　模塊主菜單示意圖Fig.1 The diagram of main menu

圖2　支柱參數(shù)對(duì)話框Fig.2 The dialog box of columns’ parameters

2　工程概況

該工程冷卻塔為大型雙曲間冷塔，塔高180 m，喉部相對(duì)標(biāo)高158.4 m，進(jìn)風(fēng)口標(biāo)高27.5 m，塔頂中面直徑103.48 m，喉部中面直徑102 m，底部直徑163.74 m，塔筒采用C40現(xiàn)澆鋼筋混凝土，殼體最大厚度為2.45 m，最小厚度為0.35 m，由40對(duì)矩形截面的X型支柱將風(fēng)筒與環(huán)板基礎(chǔ)連接，支柱采用C45混凝土，截面長(zhǎng)2 m，寬1.2 m，環(huán)板基礎(chǔ)采用C35混凝土，寬12 m，厚2 m。冷卻塔主要結(jié)構(gòu)尺寸如圖3所示。

圖3　冷卻塔立面示意圖(m)Fig.3 The elevation diagram b of cooling tower(m)

塔筒采用SHELL63單元建模，下部與支柱通過節(jié)點(diǎn)自由度耦合連接，支柱和環(huán)基采用BEAM188單元，二者通過剛域連接。采用自行研發(fā)的冷卻塔智能化模塊建模并進(jìn)行加載分析，建模、加載、后處理對(duì)話框分別如圖4～6所示。

圖4　有限元模型建立對(duì)話框Fig.4 The dialog box of finite element model’s establishment

圖5　施加風(fēng)荷載參數(shù)對(duì)話框Fig.5 The dialog box of applying wind load’s parameters

圖6　位移云圖對(duì)話框Fig.6 The dialog box of the displacement contour

3　施工全過程極限承載性能分析

冷卻塔屬于鋼筋混凝土結(jié)構(gòu)，本文采用最廣泛使用的線彈性理論本構(gòu)模型和雙軸荷載下的修正摩爾-庫倫破壞準(zhǔn)則進(jìn)行承載力分析。表面風(fēng)荷載采用風(fēng)洞試驗(yàn)獲得的最不利風(fēng)向角下風(fēng)壓系數(shù)和規(guī)范對(duì)稱風(fēng)壓兩種分布模式，內(nèi)壓系數(shù)取為-0.5。鑒于本文的研究目的，該冷卻塔風(fēng)洞試驗(yàn)和最不利風(fēng)壓系數(shù)分布詳見文獻(xiàn)[13]。風(fēng)振系數(shù)采用風(fēng)洞試驗(yàn)和風(fēng)振動(dòng)力計(jì)算得到的冷卻塔子午向不同高度的平均風(fēng)振系數(shù)，圖7給出了冷卻塔風(fēng)振系數(shù)沿子午向高度分布曲線。本文分四種組合工況進(jìn)行施工全過程冷卻塔極限承載能力計(jì)算：工況1為規(guī)范對(duì)稱風(fēng)壓計(jì)入內(nèi)吸力、工況2為風(fēng)洞試驗(yàn)風(fēng)壓計(jì)入內(nèi)吸力、工況3為規(guī)范對(duì)稱風(fēng)壓不計(jì)內(nèi)吸力、工況4為風(fēng)洞試驗(yàn)風(fēng)壓不計(jì)內(nèi)吸力。

圖7　冷卻塔層平均風(fēng)振系數(shù)沿子午向高度分布曲線Fig.7 The distribution curve of layer average wind vibration coefficient along meridian height for cooling tower

3.1強(qiáng)度極限承載能力

3.1.1逐級(jí)加載法說明

以10 m高度處29.7 m/s的初始設(shè)計(jì)風(fēng)速作為基本風(fēng)速逐級(jí)施加風(fēng)荷載，加載步長(zhǎng)為1～10 m/s。風(fēng)速增大至一定值時(shí)，冷卻塔筒壁局部區(qū)域混凝土產(chǎn)生拉伸破壞退出工作，此時(shí)鋼筋受拉，隨著風(fēng)速進(jìn)一步增大，冷卻塔混凝土筒壁受壓區(qū)接近極限受力狀態(tài)，此時(shí)徑向位移迅速增大[14]。

3.1.2冷卻塔整體強(qiáng)度極限承載能力

為了比較不同風(fēng)壓類型、計(jì)入內(nèi)壓效應(yīng)與否對(duì)冷卻塔極限承載能力的影響，圖8給出了四種工況下冷卻塔塔筒最大徑向位移、位移梯度及出現(xiàn)位置隨基本風(fēng)速變化曲線。

對(duì)比分析得出：① 隨著基本風(fēng)速的遞增，四種工況下塔筒最大徑向位移均迅速增加；② 在風(fēng)速遞增初始階段內(nèi)壓效應(yīng)并不明顯，但隨著風(fēng)速的進(jìn)一步增大，內(nèi)壓效應(yīng)逐漸突出，此時(shí)不計(jì)內(nèi)吸力相比計(jì)入內(nèi)吸力的徑向位移稍??；③ 規(guī)范風(fēng)壓與風(fēng)洞試驗(yàn)最不利風(fēng)壓的徑向位移相差不大，后者比前者略小5.35%左右；④ 冷卻塔的極限承載能力工況三最大，工況四次之，工況一較小，工況二最小，綜合四種工況可知該冷卻塔的臨界風(fēng)速為143 m/s，遠(yuǎn)大于當(dāng)?shù)?0年一遇基本風(fēng)速29.7 m/s；⑤ 隨計(jì)算風(fēng)速不斷增大，按規(guī)范風(fēng)壓計(jì)算塔筒的最大徑向位移出現(xiàn)位置較風(fēng)洞試驗(yàn)計(jì)算的低，且前者隨風(fēng)速呈降低趨勢(shì)，后者反之；當(dāng)采用規(guī)范對(duì)稱風(fēng)壓進(jìn)行加載時(shí)最大位移出現(xiàn)位置計(jì)入內(nèi)吸力相比不計(jì)內(nèi)吸力較低，當(dāng)采用風(fēng)洞試驗(yàn)最不利風(fēng)壓進(jìn)行加載時(shí)，計(jì)入內(nèi)吸力最大位移出現(xiàn)的位置則高于不計(jì)內(nèi)吸力，原因是規(guī)范風(fēng)壓呈對(duì)稱分布，而風(fēng)洞最不利風(fēng)壓考慮了群塔效應(yīng)，改變了表面風(fēng)壓的分布模式，故同等情況下計(jì)入內(nèi)吸力與否對(duì)塔筒徑向位移分布造成很大的差異。

圖8　冷卻塔強(qiáng)度極限承載狀態(tài)(施工模板層數(shù)為120)Fig.8 Cooling tower’s ultimate bearing state(construction template is 120)

3.1.3施工全過程冷卻塔強(qiáng)度極限承載能力

建立不同施工高度冷卻塔模型，按上述加載方法計(jì)算施工全過程冷卻塔的臨界風(fēng)速并進(jìn)行對(duì)比分析，通過設(shè)置施工模板數(shù)為20、40、60、80、100、120層實(shí)現(xiàn)施工高度的遞增。圖9給出四種工況下施工全過程中冷卻塔的極限承載狀態(tài)。

圖9　施工全過程中冷卻塔強(qiáng)度極限承載狀態(tài)Fig.9 Strength ultimate bearing state in construction process

圖中結(jié)果表明：① 隨著施工模板數(shù)遞增，冷卻塔的臨界風(fēng)速迅速降低，施工模板數(shù)達(dá)到40層后，臨界風(fēng)速下降速率減緩； ② 初期施工階段工況四極限承載能力最大，后期施工階段工況四與工況三極限承載能力相當(dāng)；③ 綜合四種工況和不同施工高度下冷卻塔的臨界風(fēng)速，可知該冷卻塔在整個(gè)施工過程中的最小強(qiáng)度極限臨界風(fēng)速為143 m/s左右，遠(yuǎn)大于當(dāng)?shù)鼗驹O(shè)計(jì)風(fēng)速。

表1給出了工況二施工全過程中冷卻塔的極限承載狀態(tài)，對(duì)比發(fā)現(xiàn)： ① 不同施工階段最大徑向位移出現(xiàn)位置沿環(huán)向基本不變，沿子午向逐漸上升；② 臨界風(fēng)速隨模板數(shù)遞增逐漸減小，徑向最大位移逐漸增大，最大值為2.272 m，當(dāng)基本風(fēng)速超過臨界風(fēng)速，冷卻塔將達(dá)到極限承載狀態(tài)，繼而發(fā)生倒塌。

表1　工況二施工全過程中冷卻塔極限承載狀態(tài)

3.2屈曲失穩(wěn)極限承載能力

3.2.1屈曲系數(shù)法說明

建立不同施工高度冷卻塔模型，對(duì)其施加不同組合工況荷載然后進(jìn)行屈曲分析，獲取該工況下冷卻塔整體失穩(wěn)的屈曲系數(shù)，進(jìn)而獲得對(duì)應(yīng)施工高度下冷卻塔的極限承載能力即臨界風(fēng)速。本節(jié)將不考慮混凝土齡期變化和考慮混凝土齡期變化兩者進(jìn)行對(duì)比研究?；炷翉椥阅Ａ侩S齡期變化關(guān)系為：

(1)

式中：Ec(t)和Ec分別為齡期t天和28天時(shí)塔筒C40的彈性模量；βt為系數(shù)，s取決于水泥種類，普通和快硬水泥取0.25；t為混凝土齡期(d)。

3.2.2自重+表面風(fēng)荷載作用

不考慮內(nèi)吸力在重力和表面風(fēng)荷載作用下計(jì)算不同高度冷卻塔的屈曲失穩(wěn)極限承載能力，屈曲荷載組合為G+λW(G表示自重，W表示風(fēng)荷載，λ表示屈曲系數(shù))。表2給出了施工全過程中按風(fēng)洞試驗(yàn)和規(guī)范風(fēng)荷載作用下冷卻塔屈曲失穩(wěn)極限承載能力。

對(duì)比分析發(fā)現(xiàn)：① 隨著施工高度的上升，臨界風(fēng)速逐漸減??；② 當(dāng)該冷卻塔施工成型時(shí)最小臨界風(fēng)速為110.70 m/s遠(yuǎn)大于當(dāng)?shù)?0年一遇最大基本風(fēng)速29.7 m/s；③ 對(duì)比兩種風(fēng)荷載作用結(jié)果可知不同施工狀態(tài)下二者臨界風(fēng)速相近，從而驗(yàn)證了規(guī)范[15]中對(duì)雙曲冷卻塔平均風(fēng)壓分布系數(shù)的規(guī)定的合理性與正確性；④ 對(duì)比可發(fā)現(xiàn)計(jì)入齡期變化后屈曲模態(tài)與不計(jì)時(shí)相似，但由于上部混凝土相對(duì)較薄弱，故屈曲系數(shù)低于不計(jì)入齡期變化的冷卻塔的屈曲系數(shù)，且隨施工模板數(shù)的增加，屈曲系數(shù)降低量減少。

表2　風(fēng)與自重作用下施工全過程冷卻塔屈曲極限承載能力

3.2.3自重+表面風(fēng)荷載+內(nèi)吸力作用

在重力、表面風(fēng)荷載及內(nèi)吸力共同作用下計(jì)算了不同高度冷卻塔的屈曲失穩(wěn)極限承載能力，圖10給出冷卻塔的臨界風(fēng)速隨施工模板層數(shù)遞增變化曲線。

對(duì)比研究發(fā)現(xiàn)：① 四種工況下冷卻塔的臨界風(fēng)速均隨施工高度的增加而減小，計(jì)入齡期變化時(shí)施加風(fēng)洞風(fēng)荷載且計(jì)入內(nèi)吸力的屈曲失穩(wěn)臨界風(fēng)速最小，為100.7 m/s，遠(yuǎn)大于當(dāng)?shù)鼗撅L(fēng)速29.7 m/s；② 規(guī)范風(fēng)荷載和風(fēng)洞試驗(yàn)風(fēng)荷載作用下的屈曲失穩(wěn)臨界風(fēng)速極為接近，計(jì)入混凝土彈性模量隨施工齡期變化時(shí)的臨界風(fēng)速小于不計(jì)齡期變化的臨界風(fēng)速；③ 在某一特定施工高度下，內(nèi)吸力的計(jì)入使冷卻塔的計(jì)算臨界風(fēng)速減小，因?yàn)閮?nèi)吸力的存在對(duì)冷卻塔的穩(wěn)定性起著不利的影響，使塔體在較低的風(fēng)速下更易發(fā)生屈曲失穩(wěn)達(dá)到極限承載狀態(tài)。

圖10　冷卻塔臨界風(fēng)速隨施工高度變化圖Fig.10 Critical wind velocity of different construction height for cooling tower

4　結(jié)　論

針對(duì)冷卻塔極限承載力問題，進(jìn)行了基于APDL和UIDL語言的ANSYS二次開發(fā)，采用逐級(jí)加載法和屈曲系數(shù)法對(duì)某大型冷卻塔施工全過程中強(qiáng)度和穩(wěn)定極限承載性能進(jìn)行了對(duì)比研究，主要結(jié)論如下：

(1) 實(shí)例計(jì)算表明本文研發(fā)的冷卻塔結(jié)構(gòu)分析智能化模塊能夠很好地應(yīng)用于實(shí)際工程，實(shí)現(xiàn)了冷卻塔有限元分析的參數(shù)化和可視化，大大提高了用戶的可接受程度和分析效率；

(2) 隨著施工高度的逐漸增大，冷卻塔極限承載能力迅速降低，按規(guī)范風(fēng)壓計(jì)算的極限承載能力相比按風(fēng)洞試驗(yàn)最不利風(fēng)壓計(jì)算的結(jié)果較?。粌?nèi)壓效應(yīng)對(duì)極限承載性能的影響并不突出，冷卻塔計(jì)入內(nèi)壓時(shí)的計(jì)算結(jié)果比不計(jì)內(nèi)壓時(shí)略?。挥?jì)入混凝土彈性模量隨施工齡期變化時(shí)的屈曲失穩(wěn)臨界風(fēng)速小于不計(jì)齡期變化的臨界風(fēng)速，但二者屈曲模態(tài)相似；

(3) 冷卻塔屈曲失穩(wěn)臨界風(fēng)速小于強(qiáng)度破壞臨界風(fēng)速，故在施工過程中，若風(fēng)速能夠達(dá)到最小臨界風(fēng)速，冷卻塔首先進(jìn)入屈曲失穩(wěn)極限狀態(tài)，風(fēng)速進(jìn)一步增大筒壁受壓區(qū)接近極限受力狀態(tài)，冷卻塔達(dá)到強(qiáng)度極限承載能力；

(4) 針對(duì)本文實(shí)際工程，在不考慮施工荷載的前提下，綜合計(jì)入內(nèi)吸力與否兩種狀態(tài)的各施工階段四種工況下的強(qiáng)度及失穩(wěn)承載能力計(jì)算結(jié)果，大型冷卻塔的最小臨界風(fēng)速為100.7 m/s左右，遠(yuǎn)大于當(dāng)?shù)?0年一遇最大風(fēng)速29.7 m/s。

[1] 沈國(guó)輝, 余關(guān)鵬, 孫炳楠,等. 大型冷卻塔雙塔干擾的風(fēng)洞試驗(yàn)研究[J]. 振動(dòng)與沖擊, 2011, 30(3): 109-114.

SHEN Guohui, YU Guanpeng, SUN Bingnan, et al.Study of interference effects on two large hyperbolic cooling towers using wind tunnel test[J]. Journal of Vibration and Shock, 2011, 30(3): 109-114.

[2] 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(3): 30-39.

[3] 鄒云峰, 李壽英, 牛華偉, 等. 雙曲冷卻塔等效靜力風(fēng)荷載中國(guó)規(guī)范適應(yīng)性研究[J]. 振動(dòng)與沖擊, 2013, 32(11): 100-105.

ZOU Yunfeng, LI Shouying, NIU Huawei, et al. Adaptability study of China’s code for equivalent static wind load of hyperbolic cooling towers[J]. Journal of Vibration and Shock, 2013, 32(11): 100-105.

[4] 柯世堂, 侯憲安, 趙林,等. 超大型冷卻塔風(fēng)荷載和風(fēng)振響應(yīng)參數(shù)分析：自激力效應(yīng)[J]. 土木工程學(xué)報(bào),2012,45(12): 45-53.

KE Shitang, HOU Xianan, ZHAO Lin, et al.Parameter analysis of wind loads and wind induced responses for super-large cooling towers: self-excited force effect[J]. China Civil Engineering Journal,2012, 45(12): 45-53.

[5] 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.

[6] 秦文科, 陳明祥, 周劍波. 基于有限元的雙曲線型冷卻塔結(jié)構(gòu)強(qiáng)度及穩(wěn)定性分析[J]. 武漢大學(xué)學(xué)報(bào),2012,45(7):209-212.

QIN Wenke, CHEN Mingxiang, ZHOU Jianbo. Structual strength and stability analyses of a hyperbolic cooling tower based on finite elements[J]. Engineering Journal of Wuhan University,2012,45(7):209-212.

[7] 趙海, 白國(guó)良, 徐亞洲, 等. 初始缺陷對(duì)超大型雙曲冷卻塔極限承載力的影響[J]. 西安建筑科技大學(xué)學(xué)報(bào), 2011, 43(6): 777-783.

ZHAO Hai, BAI Guoliang, XU Yazhou, et al. The influence of initial imperfections on the ultimate bearing capacity of super-large hyperbolic cooling towers[J]. Xi’an Univ of Arch & Tech, 2011, 43(6): 777-783.

[8] 邵亞會(huì), 柯世堂, 葛耀君, 等. 超大型排煙冷卻塔強(qiáng)度及穩(wěn)定性能分析[J]. 工業(yè)建筑, 2014, 44(3): 80-84.

SHAO Yahui, KE Shitang, GE Yaojun, et al. The strength and stability analysis of super large cooling towers with gas flue[J]. Industrial Building, 2014, 44(3): 80-84.

[9] 李輝, 李龍華, 彭旭軍. 雙曲冷卻塔施工過程塔筒整體穩(wěn)定性分析[J]. 施工技術(shù), 2011, 40(6): 422-424.

LI Hui, LI Longhua, PENG Xujun.Overall stability analysis of hyperbolic cooling tower during construction[J]. Construction Technolocy, 2011, 40(6): 422-424.

[10] 柯世堂, 趙林, 張軍峰, 等. 電廠超大型排煙塔風(fēng)洞試驗(yàn)與穩(wěn)定性分析[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2011, 43(2): 114-118.

KE Shitang, ZHAO Lin, ZHANG Junfeng, et al. Wind tunnel test and stability performance analysis of super large cooling towers with gas flue in power plant[J]. Journal of Harbin Industitute of Technology, 2011, 43(2): 114-118.

[11] 龔曙光, 謝桂蘭, 黃云清. ANSYS參數(shù)化編程與命令手冊(cè)[M]. 北京: 機(jī)械工業(yè)出版社, 2009.

[12] 程銘, 吳志軍, 馮平法. 基于UIDL的ANSYS用戶菜單定制技術(shù)[J]. 機(jī)械制造, 2006, 44(10): 44-46.

CHENG Ming, WU Zhijun, FENG Pingfa. ANSYS user menu customization technology based on UIDL[J]. Machinery Manufacture, 2006, 44(10): 44-46.

[13] 柯世堂. 內(nèi)蒙古京能盛樂2×350 MW工程大型冷卻塔群風(fēng)洞試驗(yàn)與風(fēng)振分析總研究報(bào)告[R]. 南京航空航天大學(xué)結(jié)構(gòu)風(fēng)工程研究中心, 2014.

[14] 趙林, 李鵬飛, 葛耀君. 等效靜風(fēng)荷載下超大型冷卻塔受力性能分析[J]. 工程力學(xué), 2008, 25(7): 79-86.

ZHAO Lin, LI Pengfei, GE Yaojun. Numerical investigation on performance of super large cooling towers under equivalent static wind load[J]. Engineering Machannics,2008, 25(7):79-86.

[15] DL/T 5339—2006.火力發(fā)電廠水工設(shè)計(jì)規(guī)范[S]. 北京: 中國(guó)電力出版社, 2006.

Wind-induced limit bearing capacity in the whole construction process for a cooling tower based on ANSYS secondary development

DU Lingyun, KE Shitang

(Department of Civil Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

In order to investigate wind-induced limit bearing capacity of strength and stability in the whole construction process for a cooling tower. ANSYS secondary development of corresponding function interface module by using UIDL and APDL was conducted, based on which an intelligent finite element model was built. Wind tunnel and standard wind load were imposed on it. Wind-induced limit bearing capacity of strength and bulking instability in the whole construction process was calculated. The effects of concrete age change on buckling stability were then compared. It is founded that the critical wind speed of the cooling tower decreases rapidly in the initial stage of construction, and then gradually becomes smooth. Wind-induced limit bearing capacity of strength is larger than the wind-induced limit bearing capacity of bulking instability. With including concrete age change, the buckling mode is similar to the result without considering concrete age change, but the critical wind speed decreases. The overall minimum critical wind speed is 100.7m/s, far more than the maximum, wind speed of fifty years in the local(29.7 m/s). The results demonstrate that the cooling tower structure intelligent analysis module of secondary development is practical and efficient. Parametric and visual finite-element analysis of limit bearing capacity in the whole construction process for the cooling tower has been realized.

large cooling tower; secondary development; whole construction process; buckling instability; ultimate bearing capacity

國(guó)家自然科學(xué)基金(51208254);中國(guó)博士后科學(xué)基金(2013M530255;1202006B)

2015-04-07修改稿收到日期：2015-09-06

杜凌云 女，碩士生，1993年生

柯世堂 男，博士，副教授，國(guó)家一級(jí)注冊(cè)結(jié)構(gòu)師，1982年生E-mail:keshitang@163.com

TU279.7+41

A

10.13465/j.cnki.jvs.2016.16.027