申顯柱，李水生，鮑 偉，張 凱

(中國電建集團(tuán)貴陽勘測設(shè)計研究院有限公司，貴州 貴陽 550081)

在水利水電工程中，月牙肋岔管具有受力明確合理、設(shè)計方便、水流流態(tài)好、水頭損失小、結(jié)構(gòu)可靠、制作安裝容易等優(yōu)點(diǎn)，故被廣泛采用。月牙肋岔管傳統(tǒng)的設(shè)計流程：①結(jié)合主支管參數(shù)，依據(jù)規(guī)范及設(shè)計經(jīng)驗(yàn)采用AutoCAD初擬月牙肋岔管體型參數(shù)；②采用Excel表格按照結(jié)構(gòu)力學(xué)的方法進(jìn)行結(jié)構(gòu)計算驗(yàn)證；③根據(jù)計算結(jié)果調(diào)整體型后再進(jìn)行結(jié)構(gòu)計算，直至體型設(shè)計及結(jié)構(gòu)計算均滿足要求；④對于重要的工程尚需有限元法進(jìn)行計算分析；⑤最后采用AutoCAD繪制二維圖和采用Word編寫結(jié)構(gòu)計算書。

實(shí)踐表明，因?yàn)樵卵览卟砉芙Y(jié)構(gòu)較復(fù)雜，傳統(tǒng)設(shè)計過程以人工為主，設(shè)計過程時間長、易出錯、步驟多、內(nèi)容繁瑣，所以體型參數(shù)設(shè)計、結(jié)構(gòu)計算分析和二維施工圖繪制一直是設(shè)計工作中最耗工時的。雖然結(jié)構(gòu)設(shè)計、有限元計算等有相關(guān)的計算軟件，但是整體設(shè)計的自動化程度不高，反復(fù)計算和調(diào)整體型的工作量很大；工程制圖也是使用AutoCAD人工繪制，然后人工編寫結(jié)構(gòu)計算書。

本文提供了一種基于Catia-V5三維軟件設(shè)計平臺的埋藏式月牙肋岔管體型設(shè)計、結(jié)構(gòu)計算、編寫計算書、工程制圖等全過程輔助設(shè)計系統(tǒng)，能快速、自動化和規(guī)范化的完成月牙肋岔管設(shè)計工作，從而提高了生產(chǎn)效率和設(shè)計質(zhì)量，且降低了月牙肋岔管的設(shè)計難度。月牙肋岔管結(jié)構(gòu)計算編程及自動編寫計算書功能主要依托Visual Basic6.0軟件，參數(shù)化三維模板制作及二維圖繪制功能主要依托Catia-V5軟件，圍巖和岔管聯(lián)合承載有限元計算方法功能主要依托Ansys和Hypermesh等軟件進(jìn)行開發(fā)。月牙肋岔管結(jié)構(gòu)型式有對稱Y形岔管和非對稱Y形岔管，因篇幅有限，僅針對對稱Y形岔管并考慮鋼管與回填混凝土之間脫空情況進(jìn)行闡述。月牙肋岔管設(shè)計技術(shù)路線流程如圖1所示。

1 結(jié)構(gòu)計算編程及自動編寫計算書

1.1 編程依據(jù)及計算內(nèi)容

按結(jié)構(gòu)力學(xué)法，以SL 281—2003《水電站壓力銅管設(shè)計規(guī)范》[1]為依據(jù)，利用VB6.0語言進(jìn)行編程，實(shí)現(xiàn)月牙肋岔管選材、體型設(shè)計、結(jié)構(gòu)計算及自動編寫標(biāo)準(zhǔn)化計算書的功能，并包含岔管結(jié)構(gòu)計算程序結(jié)果與Catia-V5三維模板的鏈接功能。

1.2 界面設(shè)計及成果展示

(1)首先必須按照預(yù)先給定的用戶名及密碼登錄軟件界面，如圖2所示，確保了軟件的安全性。

(2)然后進(jìn)入到具體的岔管幾何參數(shù)輸入界面，依次點(diǎn)擊左邊計算樹中的各模塊，如圖3所示，體型參數(shù)輸入→物理參數(shù)輸入→工況荷載輸入→應(yīng)力計算參數(shù)→結(jié)構(gòu)應(yīng)力計算。

圖3 計算參數(shù)輸入界面

(3)最后點(diǎn)擊“結(jié)果輸出”可以進(jìn)入到計算結(jié)果輸出界面，本軟件提供了兩種結(jié)果輸出方式，一種是以Word算稿的型式進(jìn)行輸出，如圖4所示，另外一種是以三維軟件Catia-V5的型式輸出三維空間結(jié)構(gòu)，如圖5所示。

圖4 Word格式文件輸出界面

圖5 Catia三維格式文件輸出界面

綜上所述，月牙肋岔管程序計算軟件可以方便地實(shí)現(xiàn)結(jié)構(gòu)計算，并且還可以生成相應(yīng)的標(biāo)準(zhǔn)化計算書和Catia-V5的三維初步模型，和傳統(tǒng)的Excel表格計算和人工編寫Word計算書相比較，不僅降低了表格計算的錯誤率，而且減少了計算書編寫中的許多重復(fù)工作，從而有效地提高了工作效率。

2 參數(shù)化三維模板制作及二維圖繪制

2.1 主要研究目的及內(nèi)容

2.1.1研究目的

基于Catia-V5創(chuàng)建滿足設(shè)計要求的參數(shù)化模板，并與牙肋岔管結(jié)構(gòu)計算編程模塊成果自動鏈接和更新，再通過局部調(diào)整其他參數(shù)，最后形成三維效果圖展示、繪制二維圖和計算工程量，同時也將Catia-V5結(jié)構(gòu)模型導(dǎo)入Ansys中完成網(wǎng)格劃分工作。

2.1.2模板制作內(nèi)容

月牙肋岔管三維模板制作以約束造型為核心，以尺寸驅(qū)動為特征，實(shí)現(xiàn)參數(shù)化、標(biāo)準(zhǔn)化及模板化設(shè)計，包括體形設(shè)計和管節(jié)展開。模板文件分別為有限元計算模型和管節(jié)展開模型。

2.1.3模板制作思路

(1)利用Catia-V5的參數(shù)化設(shè)計和知識工程功能，建模時將相關(guān)的設(shè)計參數(shù)提取出來，放在模型樹中，定義參數(shù)名稱以便參數(shù)化設(shè)計，只要更改模型樹中的參數(shù)，就可以很容易地更新模型。(2)將規(guī)范中的限制性條件寫成規(guī)則，由軟件自動判斷岔管的體型參數(shù)是否滿足規(guī)范要求。(3)利用Catia-V5工程制圖模塊快速將岔管管節(jié)展開圖投影生成平面線框圖，以方便岔管施工圖繪制。

2.2 三維模板制作及成果展示

月牙肋岔管模板建模過程：規(guī)劃岔管的結(jié)構(gòu)分解→參數(shù)規(guī)劃→草圖繪制→三錐體中間面建?！甙逯虚g面建模→生成實(shí)體模型→有限元模型→出圖模型→知識工程。

主要成果有三錐體的曲面(鋼管中間面)如圖6所示，三維模型如圖7所示，二維圖如圖8所示。

圖6 對稱月牙肋岔管三錐體曲面

圖7 對稱月牙肋岔管三維模型圖

圖8 對稱月牙肋岔管二維圖

綜上所述，采用傳統(tǒng)設(shè)計方法進(jìn)行月牙肋岔管設(shè)計時，體型設(shè)計幾乎都采用數(shù)解法，即通過已經(jīng)定型的公式來計算岔管各點(diǎn)坐標(biāo)、參數(shù)等，從而最后繪出岔管幾何體型。但利用Catia-V5軟件的參數(shù)化建模技術(shù)進(jìn)行月牙肋岔管設(shè)計，具有準(zhǔn)確、快速、直觀、修改方便等優(yōu)點(diǎn)，減少了設(shè)計者工作量和人為計算誤差，提高了設(shè)計效率和質(zhì)量。

3 圍巖和岔管聯(lián)合承載有限元計算方法

3.1 研究目的及計算原則

3.1.1研究目的

由Catia-V5參數(shù)化模板生成有限元模型后，利用Ansys等軟件進(jìn)行有限元分析計算，掌握圍巖和岔管聯(lián)合承載有限元計算理論及方法。

3.1.2計算原則

以有限元計算軟件Ansys作為主要平臺，輔以Hypermesh軟件進(jìn)行相關(guān)處理?？紤]圍巖分擔(dān)內(nèi)水壓力，考慮鋼管與回填混凝土之間脫空。計算過程中岔管管壁基于薄壁理論假設(shè)采用殼單元進(jìn)行模擬，圍巖作用采用彈性抗力系數(shù)來簡化計算，圍巖與管壁之間相互作用采用接觸單元進(jìn)行模擬。

3.2 圍巖和岔管聯(lián)合承載有限元計算方法

3.2.1埋藏式月牙肋岔管的模擬方法

當(dāng)需要考慮鋼襯與回填混凝土之間脫空時，裂縫的寬度和管道直徑在尺寸上的巨大量級差異給有限元計算帶來了困難。針對地下埋管計算中山體網(wǎng)格劃分工作量大的問題，經(jīng)過研究使用殼單元和點(diǎn)對點(diǎn)接觸單元的方法來模擬地下埋管的方法。這種簡化方式避免了對山體的網(wǎng)格劃分，提高了計算效率。原理簡圖如圖9所示，接觸單元的受力特性曲線如圖10所示。

圖9 混凝土與岔管接觸單元示意圖

圖10 接觸單元應(yīng)力應(yīng)變關(guān)系

第一步：利用前處理軟件Hypermesh劃分岔管本身網(wǎng)格，將其導(dǎo)入Ansys軟件。第二步：施加邊界條件(最外邊使用全約束邊界)。第三步：利用提供的APDL代碼施加圍巖邊界條件，并進(jìn)行計算。

3.2.2邊界條件及有限元網(wǎng)格劃分相關(guān)要求

計算范圍選取原則：模型的截斷邊界要求距離岔管不少于一倍的公切球直徑，靠近岔管(管徑范圍內(nèi))管節(jié)要求模擬其體型。遠(yuǎn)離岔管核心區(qū)域的鋼管端部采用全約束邊界條件。

網(wǎng)格劃分：岔管的有限元網(wǎng)格劃分通常分為兩個部分，一部分為月牙肋的網(wǎng)格，另一部分為岔管主體的網(wǎng)格劃分。鋼管相對較薄，在有限元計算中使用殼單元來模擬；肋板相對較厚，且肋板對平面外應(yīng)力有一定要求，在計算過程中使用實(shí)體單元進(jìn)行模擬。由于模型中同時存在殼單元和實(shí)體單元，在設(shè)置單元大小時候應(yīng)滿足殼單元的薄殼假設(shè)，盡量避免剪切自鎖現(xiàn)象。肋板基本不承平面外受彎矩作用，并且考慮到網(wǎng)格質(zhì)量問題，在網(wǎng)格劃分過程中為避免肋板單元長寬比失衡，在厚度方向僅布置一層單元。管壁網(wǎng)格示意圖如圖11所示。

圖11 管壁網(wǎng)格示意圖

圍巖約束條件：圍巖約束條件采用圍巖彈性抗力系數(shù)來模擬。每個節(jié)點(diǎn)上對應(yīng)的接觸剛度為k=f×A，f為圍巖彈性抗力系數(shù)；A為節(jié)點(diǎn)的特征面積。由于有限元網(wǎng)格在離散過程中不能保證每個節(jié)點(diǎn)代表的面積完全相同，且各個接觸單元方向不同，圍巖約束條件手動施加實(shí)現(xiàn)難度較大。實(shí)施過程中使用Ansys的內(nèi)部APDL編程語言進(jìn)行編程實(shí)現(xiàn)圍巖約束條件的施加。

綜上所述，隨著國內(nèi)水電事業(yè)的發(fā)展，壓力鋼管HD值越來越大，如何減小鋼板厚度，降低制造安裝難度，已經(jīng)成為了比較突出的問題。除了優(yōu)化岔管體型及提高鋼材強(qiáng)度之外，對于地下埋管來說，合理利用圍巖的承載能力也是一個重要的措施，本節(jié)內(nèi)容針對埋藏式月牙肋岔管的有限元計算分析方法進(jìn)行了說明，以期為設(shè)計人員提供參考。

4 工程實(shí)例

4.1 結(jié)構(gòu)計算

(1)基本情況

結(jié)合石埡子水電站[2]埋藏式月牙肋岔管進(jìn)行實(shí)例驗(yàn)證，鋼管由直徑為5.7m的主管經(jīng)對稱Y形月牙肋岔管后分岔為兩條直徑3.9m的支管，管外回填C20素混凝土。原設(shè)計過程中，岔管處鋼管承擔(dān)140.0m內(nèi)壓，綜合考慮各方面因素，假設(shè)圍巖分擔(dān)35.36m內(nèi)壓，圍巖分擔(dān)率20.16%。鋼材采用Q345R，外水壓力30m。

(2)計算結(jié)果

岔管結(jié)構(gòu)分岔角65°，管壁計算厚度30mm，肋板厚度68mm，主管半徑2866mm，支管半徑1966mm。公切球半徑3150mm，支管腰線轉(zhuǎn)角15°，主管腰線轉(zhuǎn)角10°。加勁環(huán)尺寸采用130×22mm，間距為1000mm。

采用本文結(jié)構(gòu)計算編程軟件，計算成果及結(jié)論：管壁膜應(yīng)力151.02MPa，主錐與支岔錐連接處的邊緣應(yīng)力207.69MPa，主錐與主管連接處的邊緣應(yīng)力239.03MPa，支岔錐與支管連接處的邊緣應(yīng)力203.71MPa，肋板正應(yīng)力96.87MPa，剪應(yīng)力64.24MPa，第一主應(yīng)力128.88MPa，第二主應(yīng)力-32.02MPa，強(qiáng)度及應(yīng)力校核滿足要求。加勁環(huán)管壁間臨界外壓為6.96MPa，加勁環(huán)的穩(wěn)定臨界外壓為1.88MPa，鋼管外壓穩(wěn)定計算和抗外壓強(qiáng)度計算均滿足設(shè)計要求。計算成果與原設(shè)計過程應(yīng)力值相差約5%左右，精度滿足要求。

4.2 三維設(shè)計及二維圖

輸入?yún)?shù)生成三維、二維模型如圖12—13所示。

圖12 三維模型展示圖

圖13 二維圖

通過對岔管的結(jié)構(gòu)平面圖、肋板立面圖、主岔管展開圖、支岔管展開圖等施工圖與模型投圖對比，吻合度很好。肋板立面圖施工圖與模型投圖對比如圖14所示，圖片中紅色線條為實(shí)際施工圖，黑色線條為三維模型投圖，完全相同。

圖14 肋板施工圖與模型投圖對比

4.3 有限元計算

結(jié)合石埡子水電站岔管實(shí)際情況，驗(yàn)證岔管應(yīng)力以及圍巖分擔(dān)等情況的計算精度。岔管位于II類圍巖中，彈性抗力系數(shù)75MPa/cm，回填混凝土與鋼襯的脫空按1.14mm(縫隙計算中的最終總縫隙取值)考慮，內(nèi)水壓力175.36m。有限元計算應(yīng)力如圖15所示。

圖15 Ansys平臺計算結(jié)果

計算結(jié)果分析，①特征點(diǎn)最大應(yīng)力129.28MPa，肋板最大應(yīng)力94MPa，在考慮1.14mm的脫空后岔管應(yīng)力為明管應(yīng)力的60%左右。②原算稿中圍巖分擔(dān)比例明顯要低于有限元計算結(jié)果，因此可以通過有限元計算分析為岔管設(shè)計過程中圍巖分擔(dān)比例的選擇提供依據(jù)，結(jié)合實(shí)際情況確定合理的圍巖分擔(dān)比，降低岔管現(xiàn)場制作安裝的難度和節(jié)省投資。③在考慮脫空情況下，石埡子水電站岔管的最大應(yīng)力為130MPa，計算結(jié)果與傳統(tǒng)算稿基本一致，計算精度滿足工程需求。

5 結(jié)語

埋藏式月牙肋岔管輔助設(shè)計系統(tǒng)是基于Catia-V5軟件平臺及應(yīng)用較成熟的VB6.0和Ansys等商業(yè)軟件開發(fā)的，其研究成果內(nèi)容豐富、齊全，各功能模塊可視化程度較高，簡單明了，易懂易學(xué)，計人員輸入初始參數(shù)后，計算機(jī)就能自動完成設(shè)計工作。這個設(shè)計系統(tǒng)的開發(fā)不僅極大地提高了埋藏式月牙肋岔管的生產(chǎn)效率和產(chǎn)品質(zhì)量，而且降低了岔管設(shè)計技術(shù)難度，具有推廣應(yīng)用價值，供讀者借鑒。