盛 超,杜 娟

(宿遷學(xué)院建筑工程系)

0 引 言

受自然環(huán)境變化影響的連續(xù)剛構(gòu)箱梁橋,在外界溫度變化作用下,其表面溫度迅速上升或下降,但內(nèi)部大部分區(qū)域仍處于原來的溫度狀態(tài),從而在箱梁中形成較大的溫度梯度,進(jìn)而引起溫度變形,這種溫差作用下產(chǎn)生的變形,受到箱身截面的纖維約束時(shí)將產(chǎn)生溫度應(yīng)力。

由世界上最大的有限元分析軟件公司之一的美國 ANSYS開發(fā)的 ANSYS軟件是融結(jié)構(gòu)、流體、電場、磁場、聲場分析于一體的大型通用有限元分析軟件,其用戶涵蓋了機(jī)械、航空航天、能源、交通運(yùn)輸、土木建筑等眾多領(lǐng)域。ANSYS是這些領(lǐng)域進(jìn)行分析設(shè)計(jì)技術(shù)交流的分析平臺,是個(gè)功能強(qiáng)大的有限元分析程序,ANSYS主要由前置處理(Preprocessing)、解題程序(solution)以及后置處理 (Postprocessing)以及時(shí)間歷程等組成。

1 工程概況

某橋橋梁軸線與水流交角為 90°,橋梁全長 287.54m。主橋上部采用 35m+60m+90m+60m+35m預(yù)應(yīng)力混凝土剛構(gòu) -連續(xù)箱梁體系;主橋主墩采用雙薄壁式墩,主橋邊墩采用板式橋墩;基礎(chǔ)為鋼筋混凝土沉井基礎(chǔ)、剛性擴(kuò)大基礎(chǔ)。箱梁采用單箱單室,主跨墩頂高度為 5.0m,跨中高度2.0m,箱梁梁高按二次拋物線(y=0.001829x2)變化,頂板厚 0.25m,在墩頂范圍內(nèi)為 0.5m,底板厚度在墩頂截面為1.0m,墩頂根部為 0.7m,在跨中為 0.3m,其間按而此拋物線連續(xù)變化,腹板厚度按 0.7～0.55～0.4m變化。全橋總體布置及橫斷面形式如圖1所示。

圖 1 箱梁橫斷面圖 (單位:cm)

2 有限元模型的建立

在 ANSYS仿真分析過程中,有兩種輸入方式:人機(jī)交互輸入和用命令流的方式輸入,對于復(fù)雜的結(jié)構(gòu)用命令流或兩者配合輸入比較方便,另外,建模過程中,可以利用結(jié)構(gòu)的對稱性通過映像的方法也可以起到事半功倍的效果。

根據(jù)該連續(xù)剛構(gòu)箱梁橋的實(shí)際尺寸,用 ANSYS建立全橋模型。

ANSYS具有豐富的單元庫和材料庫,可以通過對各種荷載工況的組合,較為精確的反映出橋梁在各種因素下的綜合特征。針對不同的結(jié)構(gòu),進(jìn)行單元選擇的時(shí)候,應(yīng)盡量選擇適用于分析目的的單元,在 ANSYS的單元庫中,solid90單元為二十節(jié)點(diǎn)六面體單元,該二十節(jié)點(diǎn)單元有協(xié)調(diào)的溫度形函數(shù),每個(gè)節(jié)點(diǎn)具有一個(gè)溫度自由度,本文選用 solid90單元。進(jìn)行單元劃分時(shí),可以使用自由網(wǎng)格劃分,也可以使用映射網(wǎng)格劃分,前者對實(shí)體模型無特殊要求,后者要求實(shí)體模型具有規(guī)則的形狀。對該模型進(jìn)行單元劃分,全橋共分 181472個(gè)單元。

3 荷載的施加與有限元求解

ANSYS具有多場藕合分析功能,能夠?qū)ψ匀唤缌?、流場、熱場、磁場?shí)現(xiàn)全面分析,在 ANSYS中,熱場與結(jié)構(gòu)應(yīng)力場屬于不同的物理場,要在熱場和結(jié)構(gòu)應(yīng)力場之間建立關(guān)系需要進(jìn)行熱應(yīng)力耦合分析,熱應(yīng)力問題實(shí)際上是熱和應(yīng)力兩個(gè)物理場之間的相互作用。

在ANSYS中求解連續(xù)剛構(gòu)橋溫度效應(yīng)時(shí),首先要確定溫度梯度模式,各橋涵設(shè)計(jì)規(guī)范對溫度梯度模式都作了規(guī)定,并且都有一定程度的差異。我國公路橋涵設(shè)計(jì)通用規(guī)范JTGD60-2004規(guī)定:計(jì)算橋梁結(jié)構(gòu)由溫度引起的效應(yīng)時(shí),可采用圖 2所示的豎向溫度梯度曲線。

圖 2 中國公路規(guī)范溫度梯度模式(尺寸單位:mm)

對混凝土結(jié)構(gòu),當(dāng)梁高 H小于 400mm時(shí),A=H-100(mm);梁高H等于或大于 400mm時(shí),A=300mm。對帶混凝土橋面板的鋼結(jié)構(gòu),A=300mm,圖中t為混凝土橋面板的厚度(mm)。混凝土上部結(jié)構(gòu)和帶混凝土橋面板的鋼結(jié)構(gòu)的豎向日照反溫差為正溫差乘以-0.5。

表 1 豎向日照正溫差計(jì)算的溫度基數(shù)

在ANSYS進(jìn)行熱—結(jié)構(gòu)藕合分析時(shí),實(shí)際上是把溫差作為一種荷載施加在有限元模型上,然后施加邊界條件進(jìn)行運(yùn)算求解。

ANSYS中藕合場的分析最終可歸結(jié)為兩種不同的方法即直接藕合方法和序貫藕合方法:直接藕合解法,通過計(jì)算包含所有必需項(xiàng)的單元矩陣或單元荷載向量來實(shí)現(xiàn),利用包含所有必需自由度的藕合單元類型通過一次求解得出耦合場分析結(jié)果;序貫藕合解法按照順序進(jìn)行兩次或更多次的相關(guān)場分析,通過把第一次場分析的結(jié)果作為第二次場分析的載荷來實(shí)現(xiàn)兩種場的藕合。序貫藕合方法簡單而且適應(yīng)性強(qiáng),對熱—結(jié)構(gòu)模型限制較少,故本文運(yùn)用序貫熱—應(yīng)力藕合法對連續(xù)剛構(gòu)箱梁進(jìn)行熱—結(jié)構(gòu)藕合分析。

序貫熱—應(yīng)力藕合分析中需要將熱單元轉(zhuǎn)換為相應(yīng)的結(jié)構(gòu)單元,與熱單元 solid90相對應(yīng)的結(jié)構(gòu)單元是 solid95,Solid95是Solid90的等效結(jié)構(gòu)單元,即 3-D結(jié)構(gòu)實(shí)體二十節(jié)點(diǎn)六面體單元,該二十節(jié)點(diǎn)單元有協(xié)調(diào)的彎曲邊界形函數(shù),每個(gè)節(jié)點(diǎn)具有三個(gè)自由度(X,Y,Z三個(gè)方向的位移)。本文用熱單元solid90及與其相對應(yīng)的結(jié)構(gòu)單元solid95進(jìn)行序貫熱—應(yīng)力藕合分析。

施加溫度荷載進(jìn)行熱分析后,將溫度作為荷載施加在有限元模型上,按公路規(guī)范新規(guī)范規(guī)定的溫度模式施加溫度荷載后的有限元模型。

4 后處理

在 ANSYS中,檢查分析結(jié)果可以使用兩種后處理器:POST1通用后處理器,可用來檢查整個(gè)模型在某一荷載步和子步的結(jié)果;POST26時(shí)間歷程后處理器,可以檢查模型的指定點(diǎn)的特定結(jié)果相對于時(shí)間、頻率或其他結(jié)果項(xiàng)的變化。選用通用后處理器處理仿真計(jì)算數(shù)據(jù)。

使用通用后處理器可以把仿真計(jì)算結(jié)果用應(yīng)力等值線的方式顯示出來,通過有限元分析,求得的連續(xù)剛構(gòu)橋跨中斷面的應(yīng)力云。

使用通用后處理器還可以顯示模型仿真分析的結(jié)構(gòu)變形,通過有限元分析,求得的連續(xù)剛構(gòu)橋的變形。

5 結(jié) 語

用ANSYS有限元分析軟件求解連續(xù)剛構(gòu)箱梁橋的溫度效應(yīng),涉及到有限元模型的建立、單元類型的選擇、熱—應(yīng)力耦合分析類型選擇以及后處理等問題。我們應(yīng)該根據(jù)具體情況尋求適宜的方法,滿足我們學(xué)習(xí)研究的需要。

[1]蔣國富.大跨徑橋梁高墩日照溫度效應(yīng)的研究[D].長安大學(xué),2005.

[2]康為江.鋼筋混凝土箱梁日照溫度效應(yīng)研究[D].湖南大學(xué)碩士論文,2000.

[3]柯尊禮.大跨度PC連續(xù)箱梁橋的溫度場及其效應(yīng)分析[D].武漢理工大學(xué)碩士論文,2004.

[4]王琳,項(xiàng)貽強(qiáng),王勁豐,王建江.各國規(guī)范關(guān)于混凝土箱梁橋溫度應(yīng)力計(jì)算的分析與比較[J].公路,2004,(6):76～78.

[5]中華人民共和國行業(yè)標(biāo)準(zhǔn).JTGD60-2004公路橋涵設(shè)計(jì)通用規(guī)范[M].北京:人民交通出版社,2004.

[6]王富恥,張朝暉.ANSYS10.0有限元分析理論與工程應(yīng)用[M].北京:電子工業(yè)出版社,2006.

[7]陳精一,蔡國忠.電腦輔助工程分析 ANSYS使用指南[M].北京:中國鐵道出版社,2001.

[8]王建江,胡仁喜,劉英林等編著.ANSYS11.0結(jié)構(gòu)與熱力學(xué)有限元分析實(shí)例指導(dǎo)教程[M].北京:機(jī)械工業(yè)出版社,2008.