陳 蕾,曹邱林

(揚州大學水利科學與工程學院,江蘇 揚州225009)

1 問題的提出

劉集地涵工程的主要結構設施由穿房亭河底的地下涵洞、泄洪口涵洞及南、中、北3座閘共30孔閘門,閘總長325 m,總寬95.8 m。地下涵洞共10孔,每孔寬 3.6 m,高4.0 m,長224m,分12節(jié),底高程 13.5 m。劉集地涵的洞身結構為鋼筋混凝土箱形結構,其優(yōu)點為結構承載力較高。

傳統(tǒng)的箱涵結構計算主要為內力計算及鋼筋混凝土結構計算。內力計算包括各控制截面的彎矩、軸力及剪力的計算。常用的彎矩計算方法有查表法、彎矩分配法及迭代法等[1]。傳統(tǒng)計算方法比較繁瑣,本文采用有限元法對箱涵進行結構計算,能夠簡化繁瑣的計算過程,很簡便的得到內力計算成果,對工程進行結構復核[2]。

2 有限元基本理論方法

對結構進行有限元法分析時,首先進行結構離散化,即選擇單元形態(tài),將計算區(qū)域劃分為有限個單元;然后進行單元分析,即構造單元位移函數(shù),建立單元的等效結點荷載列陣、幾何矩陣、彈性矩陣和單元剛度矩陣等;最后進行整體分析,由單元的等效結點荷載列陣集合成整體結點荷載列陣,單元剛度陣集合為整體剛度矩陣,引進邊界約束條件,解總體剛度方程求出結點位移分量[3]。

由于三角形單元對復雜邊界有較強的適應能力,因此一般將一個二維域離散成有限個三角形單元。典型的3結點三角形單元結點編碼為i,j,m,以逆時針方向編碼為正向,每個結點有2個位移分量,如圖1所示。3結點三角形單元的自由度為6,位移模式為:

形函數(shù)Ni表達式為:

其中ai,bi,ci為與坐標有關的系數(shù),

3 劉集地涵結構計算

劉集地涵洞身共12節(jié),1、2、11、12節(jié)洞身,4～8節(jié)洞身,9、10節(jié)洞身斷面形狀不同?，F(xiàn)取1、2、11、12節(jié)洞身為例,復核其的結構強度。

圖1 3結點三角形單元

模型建立:取邊墩、中墩、縫墩、涵洞底板,涵洞頂板為一計算單元。地基按彈性地基考慮[5],涵洞的底高程為14.4 m,根據(jù)徐州水利局勘探試驗報告,選取第5層土作為地基土,為滿足計算精度要求土層厚度取10 m,選取填土高度最高的最不利截面進行內力計算。

定義材料屬性:邊墩、中墩、縫墩、涵洞底板,涵洞頂板定義為梁桿材料,取單位長度1 m進行定義,截面形狀皆是矩形,涵洞邊墩截面尺寸為1.0 m×0.8 m×5.2 m,涵洞中墩截面尺寸為1.0 m×0.6 m×5.2 m,涵洞中墩截面尺寸為1.0 m×0.6 m×5.2 m,涵洞底板截面尺寸為1.0 m×21.8 m×0.6 m,涵洞頂板截面尺寸為1.0 m×21.8m×0.6 m;混凝土標號為C20;地基定義為線彈性材料,地基土的彈性模量5 000 kPa,泊松比0.35,重度20.4 kN/m3,凝結力 c值 47.0 kPa,摩擦角 22°。

網(wǎng)格劃分:將地基劃分有限三角形單元,單元形狀應盡可能規(guī)則,內角不應相差太多,不能出現(xiàn)鈍角,網(wǎng)格劃分越細,結點越多,計算結果越精確。但網(wǎng)格加密到一定程度后計算精度的提高就不明顯,同時要兼顧機時,綜合考慮對靠近底板的區(qū)域網(wǎng)格加密,遠離底板的區(qū)域劃分稀疏,并且使相鄰分區(qū)單元的邊界節(jié)點重合[6]。

邊界條件:底部取豎向約束狀態(tài),兩側邊界取法向約束狀態(tài)。

施加荷載:涵洞荷載主要有自身重力,涵洞頂?shù)奶钔翂毫?兩側的邊荷載和涵洞內的內水壓力。本次考慮的計算工況為涵洞內滿水時的淹沒壓力流。有限元模型共剖分三角形單元1 895個,結點896個,有限元網(wǎng)格見圖2,得出的內力圖見圖3。

圖2 計算模型

圖3 結構內力圖(5孔)

根據(jù)圖3所示的內力圖,對涵洞結構強度進行復核,結果如表1。

4 結 語

(1)劉集地涵洞身的邊墩、中墩、縫墩、涵洞頂板、涵洞底板的配筋均滿足要求[7]。

(2)有限元分析是近似求解微分方程的一種方法,使用有限元計算軟件進行涵洞結構分析計算簡便,易于工程應用,具有較好的計算精度,能滿足工程計算要求。和傳統(tǒng)的箱涵結構計算方法相比,有限元方法能化繁為簡,化抽象為具體,使分析計算過程清晰明了,并且適于模塊化編程操作,節(jié)省大量的工程量,有限元分析對處理復雜的實際工程具有獨到的優(yōu)越性。

表1 地下涵洞結構配筋計算成果表

[1]熊啟鈞著.涵洞[M].北京:中國水利水電出版社,2006.

[2]朱伯芳著.有限元原理與應用[M].北京:中國水利水電出版社,1998:510-514.

[3]王元漢,李麗娟著.有限元法基礎與程序設計[M].廣州:華南理工大學出版社,2001.

[4]雷曉燕著.有限元法[M].北京:中國鐵道出版社,2000.

[5]陳德亮著.水工建筑物[M].北京:中國水利水電出版社,2005:311-316.

[6]傅永華著.有限元分析基礎[M].武漢:武漢大學出版社,2003.

[7]SL191-2008.水工混凝土結構設計規(guī)范[S].北京:中國水利水電出版社,2008.