彭文件 盧良青

(貴州省交通規(guī)劃勘察設(shè)計(jì)研究院股份有限公司，貴州貴陽(yáng) 550081)

基于有限元法的過(guò)渡段沉降變形預(yù)測(cè)

彭文件 盧良青

(貴州省交通規(guī)劃勘察設(shè)計(jì)研究院股份有限公司，貴州貴陽(yáng) 550081)

依托于“新建鐵路成都至重慶客運(yùn)專線Ⅳ標(biāo)段四分部線下工程沉降與變形監(jiān)測(cè)”項(xiàng)目，借助ABAQUS 6.10有限元軟件進(jìn)行沉降變形模擬計(jì)算，探討基于ABAQUS的有限元法在客運(yùn)專線線下工程過(guò)渡段沉降變形預(yù)測(cè)方面的可行性。

客運(yùn)專線 沉降預(yù)測(cè) 有限元法 過(guò)渡段

過(guò)渡段建設(shè)是否合理關(guān)系整個(gè)路段的舒適性，處理不合理會(huì)出現(xiàn)“跳車”的現(xiàn)象，輕則影響旅客的舒適度，重則出現(xiàn)預(yù)想不到的災(zāi)害。

現(xiàn)階段對(duì)過(guò)渡段的預(yù)測(cè)常采用的方法是基于前期的監(jiān)測(cè)數(shù)據(jù)對(duì)后期的沉降變形作出預(yù)測(cè)，但因施工現(xiàn)場(chǎng)的復(fù)雜性，觀測(cè)元件時(shí)常被破壞，無(wú)法采集到沉降變形數(shù)據(jù)或采集到的數(shù)據(jù)貧乏，不能有效地對(duì)后期沉降變形作出預(yù)測(cè)。朱訓(xùn)國(guó)等人提出利用ABAQUS模擬NATM隧道施工過(guò)程[1]，孫吉主等人提出ABAQUS在軟基固結(jié)過(guò)程分析中的應(yīng)用[2]，劉學(xué)提出采用ABAQUS的隧道穩(wěn)定性分析[3]，王丙興提出基于ABAQUS的公路軟基沉降預(yù)測(cè)非概率可靠度研究[4]，夏力農(nóng)等人提出一種直接加固軟弱下臥土層地基方法的三維有限元分析[5]，前人開(kāi)始借助于ABAQUS有限元軟件來(lái)解決土木方面的問(wèn)題，而且取得的效果明顯。本文嘗試借助于ABAQUS有限元軟件來(lái)解決過(guò)渡段沉降變形預(yù)測(cè)問(wèn)題。

1 ABAQUS基本原理

ABAQUS是一套大型有限元通用軟件，具有強(qiáng)大的模擬功能，在它的內(nèi)部有多種模型，而各種模型都有著不同的用途。根據(jù)所模擬結(jié)構(gòu)物的特性選取不同的模型，選用最為經(jīng)典的Mohr-Coulomb(摩爾庫(kù)侖)模型。

Mohr-Coulomb塑性模型一般在顆粒結(jié)構(gòu)的材料以及單調(diào)荷載施壓的情況下使用，如土體結(jié)構(gòu)，該結(jié)構(gòu)與速率變化完全無(wú)關(guān)[6]。

1.1 Mohr-Coulomb模型的公式和參數(shù)

彈性階段中的Mohr-Coulomb模型，需要是線性、各向同性的，其屈服函數(shù)為

(1)

式中,Rmc(Θ，Φ)為π平面上屈服面形狀的一個(gè)變量。

(2)

式中,C是材料的粘聚力；

φ為q-p應(yīng)力面上Mohr-Coulomb屈服面的傾斜角，稱之為材料的摩擦角(見(jiàn)圖1)。

圖1 Mohr-Coulomb模型中的屈服面

在Mohr-Coulomb模型里，首先假設(shè)材料的硬化具有各向同性，并且硬化通過(guò)它的粘聚力系數(shù)來(lái)決定。同樣，C(黏聚力系數(shù))由塑性應(yīng)變、溫度或場(chǎng)變量的函數(shù)來(lái)表示。

經(jīng)典Mohr-Coulomb模型的屈服面具有明顯尖角，在計(jì)算時(shí)若使用相關(guān)聯(lián)的流動(dòng)法則，也就是塑性勢(shì)面與屈服面兩者一樣時(shí)，在屈服面的尖角處會(huì)產(chǎn)生不唯一的塑性流動(dòng)方向的情況。由于這一情況的存在，將使其整個(gè)模型的數(shù)值計(jì)算繁瑣、收斂緩慢。為了避免出現(xiàn)上述問(wèn)題，ABAQUS在進(jìn)行模擬計(jì)算時(shí)需選取如圖2所示連續(xù)光滑的橢圓函數(shù)模擬塑性勢(shì)面。

圖2 Mohr-Coulomb模型中的塑性勢(shì)面

Mohr-Coulomb模型通過(guò)控制粘聚力C的大小來(lái)模擬屈服面大小的變化，即材料的硬化或軟化。用戶需要指定出粘聚力C與等效塑性應(yīng)變之間的變化關(guān)系。

1.2 注意事項(xiàng)

①在ABAQUS 中，Mohr-Coulomb模型使用時(shí)需選用的求解方式為非對(duì)稱性求解，主要原因來(lái)自Mohr-Coulomb模型所采用塑性流動(dòng)法則具有非關(guān)聯(lián)性。

②Mohr-Coulomb模型在使用之前需進(jìn)行一定的標(biāo)定(即Calibrated)，典型的標(biāo)定方法是通過(guò)從不相同的三維實(shí)驗(yàn)中所得到極限應(yīng)力狀態(tài)，再在子午線平面上將不同的狀態(tài)定位出來(lái)，從而估算出摩擦角φ與粘聚力C。由此可見(jiàn)膨脹系數(shù)ψ的確定需慎重，才能保證標(biāo)定過(guò)程和塑性變形時(shí)體積的改變具有一致性。若塑性變形趨向于硬化，那么在進(jìn)行標(biāo)定的時(shí)候，最少需要一組三維實(shí)驗(yàn)作為硬化參數(shù)。

③在ABAQUS有限元模擬軟件之中，所選用Mohr-Coulomb模型的塑性流動(dòng)勢(shì)面是光滑的，在這點(diǎn)上它與經(jīng)典的相關(guān)聯(lián)Mohr-Coulomb準(zhǔn)則所采用的流動(dòng)勢(shì)面截然不同。因此，在ABAQUS有限元模擬軟件中采用的Mohr-Coulomb模型和經(jīng)典的Mohr-Coulomb模型兩者之間的塑性行為有一定的差異[7]。

2 模型的建立及材料參數(shù)確定

2.1 模型的建立

選取DK231+715～ DK231+804區(qū)間段，此段包含路基(DK231+715～ DK231+755，長(zhǎng)20 m)、路橋過(guò)渡段(DK231+755～ DK231+775，長(zhǎng)20 m)、橋墩臺(tái)(中心里程:DK231+780，沿里程方向長(zhǎng)9 m)作為建模對(duì)象，見(jiàn)圖3。即同一模型能共同研究路基、過(guò)渡段、橋墩三種不同結(jié)構(gòu)物沉降變形，避免了分開(kāi)建模未考慮各個(gè)不同結(jié)構(gòu)物相互之間影響的缺陷。本模型為了更好反應(yīng)實(shí)際情況，盡量在與現(xiàn)場(chǎng)工況相符合的情況下進(jìn)行簡(jiǎn)化。考慮到模型邊界的影響，將左、右部分尺寸設(shè)置為本次模型需考慮計(jì)算結(jié)構(gòu)物沿伸20 m；底部橋墩樁基底沿伸10 m(因橋墩為群樁基礎(chǔ)且為已經(jīng)伸入基巖中，故只延伸10 m)；因樁基在橫斷面上的間距為2 m且樁基直徑大小為1 m，為了簡(jiǎn)化模型故取值一半，即本次模型長(zhǎng)度、深度、厚度分別取89 m、34 m、1.5 m。

圖3 模型裝配

2.2 材料參數(shù)

材料參數(shù)設(shè)置的合理性直接影響著模型計(jì)算結(jié)果的準(zhǔn)確度，在有限元模型的計(jì)算中，材料參數(shù)的設(shè)置是否成功直接決定模型計(jì)算的成敗。本模型材料參數(shù)來(lái)源于設(shè)計(jì)圖紙、地勘資料與其他文獻(xiàn)。根據(jù)設(shè)計(jì)圖紙，橋墩樁基和承臺(tái)采用混凝土標(biāo)號(hào)為C40，而臺(tái)身則采用混凝土標(biāo)號(hào)為C35，同時(shí)各部件配筋率與配筋方式都有所不同，過(guò)渡段與路基填料也有一定的差異。根據(jù)地勘資料對(duì)土體進(jìn)行分層，查閱模型涉及到的設(shè)計(jì)圖紙與其它文獻(xiàn)，得到各部件的材料參數(shù)(見(jiàn)表1)。

表1 材料參數(shù)

2.3 分析步

本模型分成46個(gè)分析步進(jìn)行分析，分析步的大小、順序嚴(yán)格根據(jù)施工現(xiàn)場(chǎng)的施工工況設(shè)置?，F(xiàn)場(chǎng)于2012年7月5日開(kāi)始施工，2013年1月23日架梁完成，整個(gè)施工過(guò)程總計(jì)203天，且《成渝客運(yùn)專線線下工程沉降變形監(jiān)測(cè)方案》要求路基、過(guò)渡段施工完成后監(jiān)測(cè)6個(gè)月方可參加初評(píng)，故該模型總步長(zhǎng)為383，以便于分析每步對(duì)應(yīng)工程屬性的沉降變形情況。具體分析步見(jiàn)表2。

表2 分析步屬?zèng)r

2.4 荷載與網(wǎng)格

(1)邊界條件設(shè)置

為了保證底部不出現(xiàn)任何位移，對(duì)底部邊界節(jié)點(diǎn)x、y、z三個(gè)方向的位移進(jìn)行約束，相當(dāng)于固定支座。根據(jù)每個(gè)分析步的需要，從模型z方向向上分別對(duì)x兩側(cè)的邊界進(jìn)行y方向位移約束及y兩側(cè)的邊界進(jìn)行y方向位移約束。

(2)模型荷載加載方式

第一步，模型原土體在多年自重的影響下早已自然沉降完成，故需要考慮地應(yīng)力平衡，即在模擬施工前將初始應(yīng)力計(jì)算過(guò)程中產(chǎn)生的質(zhì)心6個(gè)方向的應(yīng)力導(dǎo)入模型下一步計(jì)算，相當(dāng)于對(duì)施工前原土體進(jìn)行穩(wěn)定處理。本模型中路基與過(guò)渡段都進(jìn)行分層填筑，每層填筑0.3 m，過(guò)渡段填筑高度為6.0 m(分20層填筑)，路基填筑高度為1.8 m(分6層填筑)，荷載結(jié)合分析步所屬的工況進(jìn)行逐步考慮，采用分級(jí)加載的方法激活土體單元，進(jìn)行施工過(guò)程模擬。待整個(gè)施工完成后，再利用該模型整體進(jìn)行地應(yīng)力平衡計(jì)算，借助ABAQUS模擬客運(yùn)專線線下工程工后沉降，分析工后沉降及不均勻沉降是否滿足工程設(shè)計(jì)要求。

(3)相互作用與網(wǎng)格

樁基周圍與土體、填筑層與橋墩臺(tái)后臺(tái)之間的法向采用“硬”接觸，允許接觸后分離，切向采用摩擦系數(shù)為0.53的摩擦接觸； 樁基底與土體、填筑層與橋墩臺(tái)后臺(tái)之間法向采用“硬”接觸，允許接觸后分離，切向采用無(wú)摩擦。因橋墩臺(tái)結(jié)構(gòu)物較多，原土體巖性分層也較為復(fù)雜，且原土體不是本模型重點(diǎn)考慮對(duì)象，故土體采用C3D4( 四結(jié)點(diǎn)線性四面體單元)，而過(guò)渡段與路基填筑部分采用C3D8R( 八結(jié)點(diǎn)線性六面體單元，沙漏控制，減縮積分)，便于更加準(zhǔn)確研究整個(gè)填筑區(qū)間段的沉降變形情況。整個(gè)模型節(jié)點(diǎn)總數(shù)為46 917個(gè)，如圖4所示。

圖4 模型網(wǎng)格

(4)作用荷載

荷載作用是否恰當(dāng)，直接關(guān)系到模型計(jì)算結(jié)果是否符合實(shí)際工況。荷載作用的時(shí)間與大小如表3所示。

表3 荷載作用的大小與步長(zhǎng)

3 數(shù)據(jù)提取與分析

建模完成后提交計(jì)算，待計(jì)算完成后，提取和分析其計(jì)算結(jié)果。在ABAQUS中其計(jì)算結(jié)果先以最直觀的云圖表示出來(lái)，該模型Y方向的位移量云圖如圖5。根據(jù)建模設(shè)置因素，其該模型的總個(gè)沉降變化趨勢(shì)即為Y方向的位移量，也就是說(shuō)研究該模型Y方向的位移變化量與變化趨勢(shì)即為預(yù)測(cè)的沉降量與發(fā)展趨勢(shì)。

3.1 數(shù)據(jù)的提取

根據(jù)成渝客運(yùn)專線線下工程沉降變形監(jiān)測(cè)設(shè)計(jì)要求，該模型包括7個(gè)斷面，分別為里程231 746、231 755兩個(gè)路基斷面，里程231 757、231 766、231 774三個(gè)過(guò)渡段斷面及里程231 775、231 780兩個(gè)橋墩臺(tái)斷面。

在模型橋墩臺(tái)提取橋臺(tái)身的沉降變形數(shù)據(jù)，而在路基與過(guò)渡段中沉降板埋設(shè)在基底，沉降板監(jiān)測(cè)隨著施工開(kāi)始進(jìn)行監(jiān)測(cè)，伴隨著整個(gè)施工，能反映整個(gè)沉降變形過(guò)程，故路基與過(guò)渡段中五個(gè)斷面提取沉降板的沉降變形數(shù)據(jù)與監(jiān)測(cè)變形數(shù)據(jù)進(jìn)行對(duì)比，探討ABAQUS模擬沉降變形預(yù)測(cè)的可行性。

圖5 沉降量云圖

3.2 數(shù)據(jù)的分析

(1)沉降變形趨勢(shì)

從ABAQUS模擬計(jì)算結(jié)果中逐個(gè)提取這7個(gè)斷面的每個(gè)分析步沉降預(yù)測(cè)值，并與監(jiān)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，比較每個(gè)施工時(shí)間段的沉降變化量與變化趨勢(shì)。ABAQUS模擬的7個(gè)斷面地沉降變形趨勢(shì)與監(jiān)測(cè)數(shù)據(jù)趨勢(shì)對(duì)比如圖6～圖12。

圖6 231746L1點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

圖7 231755L1點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

圖8 231757L1點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

圖9 231766L1點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

圖10 231774L1點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

圖11 231775T1點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

圖12 231780T3點(diǎn)預(yù)測(cè)值與監(jiān)測(cè)值對(duì)比

ABAQUS模擬客運(yùn)專線線下工程沉降變形的沉降趨勢(shì)與所監(jiān)測(cè)得到的沉降趨勢(shì)基本吻合。

ABAQUS模擬沉降量在前期比監(jiān)測(cè)沉降量略偏小，而在后期偏大，說(shuō)明ABAQUS模擬客運(yùn)專線線下工程沉降變形計(jì)算出來(lái)的結(jié)果偏于保守。

不管是在模擬曲線圖還是監(jiān)測(cè)曲線圖，在其中間部分出現(xiàn)突然沉降量，主要原因是運(yùn)梁車通過(guò)或架梁時(shí)期，而且可以明顯看出這部分沉降量值占整個(gè)沉降量的比重很大。

總體上來(lái)看，在缺少監(jiān)測(cè)數(shù)據(jù)或監(jiān)測(cè)數(shù)據(jù)貧時(shí)，可以借助ABAUQS來(lái)模擬客運(yùn)專線線下工程的沉降變形。

(2)沉降總量

ABAQUS模擬的這7個(gè)斷面的預(yù)測(cè)總沉降量與監(jiān)測(cè)總沉降量對(duì)比見(jiàn)表4。

表4 預(yù)測(cè)與監(jiān)測(cè)的總沉降量對(duì)比

從表4可知：

① 不管是從監(jiān)測(cè)數(shù)據(jù)還是預(yù)測(cè)數(shù)據(jù)來(lái)看，這7個(gè)斷面的總沉降量大小趨勢(shì)是從小里程到大里程先逐漸增加，再逐漸減小，即路基到過(guò)渡段總沉降量先是一個(gè)逐漸增加的趨勢(shì)，再到橋墩臺(tái)逐漸減小。其主要原因與填筑高度有直接關(guān)系，雖然橋墩臺(tái)高度大于過(guò)渡段的填筑高度，但橋墩臺(tái)沉降總量明顯比過(guò)渡段小，可見(jiàn)群樁基礎(chǔ)對(duì)減小沉降變形有顯著作用。

② ABAUQUS在模擬客運(yùn)專線線下工程沉降變形時(shí)，準(zhǔn)確度到達(dá)80%以上，其預(yù)測(cè)結(jié)果可以應(yīng)用于實(shí)際工程中。

(3)工后沉降

在ABAQUS模擬中，待全部施工完成后再考慮地應(yīng)力平衡，其計(jì)算結(jié)果中的沉降量為工后沉降，計(jì)算結(jié)果云圖如圖13。

圖13 工后沉降云圖

從圖13可以看出：

工后沉降最大沉降量出現(xiàn)在橋墩臺(tái)231775T1位置，其工后沉降為2.06 mm，橋墩臺(tái)工后沉降設(shè)計(jì)要求小于5 mm，故該工后沉降沉降量滿足沉降設(shè)計(jì)要求。

在過(guò)渡段與橋墩臺(tái)接觸位置出現(xiàn)不均勻沉降，不均勻沉降值在2 mm左右，其主要原因是建立模型時(shí)，未考慮在此處設(shè)計(jì)搭板，可在此處設(shè)計(jì)搭板來(lái)解決不均勻沉降。

4 結(jié)論

(1)ABAQUS模擬客運(yùn)專線線下工程過(guò)渡段沉降變形其變化趨勢(shì)與實(shí)際沉降變化趨勢(shì)基本吻合，其總沉降量的準(zhǔn)確度都在80%以上，且偏于保守。

(2)利用ABAQUS模擬客運(yùn)專線線下工程過(guò)渡段沉降變形，不僅不依賴于前期的監(jiān)測(cè)數(shù)據(jù)，而且可以提前預(yù)測(cè)出其沉降變形變化趨勢(shì)和沉降總量，可為前期的設(shè)計(jì)服務(wù)，優(yōu)化設(shè)計(jì)方案。

(3)在監(jiān)測(cè)數(shù)據(jù)采集時(shí)，常因施工現(xiàn)場(chǎng)的復(fù)雜性，出現(xiàn)觀測(cè)標(biāo)志被破壞的情況，從而引起沒(méi)有監(jiān)測(cè)數(shù)據(jù)及監(jiān)測(cè)數(shù)據(jù)貧，該情況也可采用ABAQUS模擬預(yù)測(cè)客運(yùn)專線線下工程沉降變形。

[1] 朱訓(xùn)國(guó)，楊慶，欒茂田.利用ABAQUS模擬NATM隧道施工過(guò)程[J].巖土力學(xué)，2006，10(S1):283-289

[2] 孫吉主，高暉.ABAQUS在軟基固結(jié)過(guò)程分析中的應(yīng)用[J].巖土力學(xué)，2007，28(8):919-922

[3] 劉學(xué).采用ABAQUS的隧道穩(wěn)定性分析[J].山西建筑，2009，3(9):312-313

[4] 王丙興.基于ABAQUS的公路軟基沉降預(yù)測(cè)非概率可靠度研究[J]. 公路交通技術(shù)，2012，8(4):1-3,9

[5] 夏力農(nóng)，廖常斌，王倩，等.一種直接加固軟弱下臥土層地基方法的三維有限元分析[J].湘潭大學(xué)自然科學(xué)學(xué)報(bào)，2012，9(3):28-33

[6] 費(fèi)康，張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京:中國(guó)水利水電出版社，2010

[7] 石亦平，周玉蓉.BAQUS有限元分析實(shí)例詳解[M].北京:機(jī)械工業(yè)出版社，2006

Settlementdeformationpredictionoftransitionsectionbasedonfiniteelementmethod

PENG Wen-jian LU Liang-qing

2014-09-25

彭文件(1987—)，男，2014年畢業(yè)于西南交通大學(xué)道路與鐵道專業(yè)，工學(xué)碩士。

1672-7479(2014)06-0037-05

U213.1+4

： A