胡 煒,譚信榮,喻 渝

(中鐵二院工程集團(tuán)有限責(zé)任公司,成都 610031)

1 研究背景

目前，隧道襯砌設(shè)計需利用多款軟件協(xié)同作業(yè)完成。具體而言，襯砌內(nèi)力計算需用到ANSYS、Midas-GTS等有限元軟件，內(nèi)力計算完成后，需將計算結(jié)果導(dǎo)出后借助Excel等軟件處理數(shù)據(jù)并完成素混凝土承載力驗算，素混凝土不滿足設(shè)計要求時，還需在Excel中反復(fù)試算配筋[1-5]。多個程序協(xié)同作業(yè)導(dǎo)致軟件之間多次反復(fù)的數(shù)據(jù)相互傳輸，不僅嚴(yán)重降低工作效率，且容易引起數(shù)據(jù)混淆和計算錯誤，中間任一環(huán)節(jié)出了偏差，就可能導(dǎo)致設(shè)計結(jié)果出現(xiàn)較大差異[6-10]。

此外，鐵路隧道設(shè)計方法已由安全系數(shù)法向概率極限狀態(tài)法轉(zhuǎn)軌[11-16]，而目前國內(nèi)缺乏成熟的基于極限狀態(tài)法的襯砌結(jié)構(gòu)設(shè)計軟件。欲實現(xiàn)鐵路隧道襯砌的極限狀態(tài)法設(shè)計，不僅需采用多軟件協(xié)同作業(yè)的傳統(tǒng)方法，還需要重新編制有限元計算命令流及承載力驗算和配筋計算單，進(jìn)一步降低了隧道襯砌設(shè)計效率[17-21]。

因此，開發(fā)一款基于VB6.0開發(fā)的可實現(xiàn)鐵路隧道襯砌極限狀態(tài)設(shè)計的隧道襯砌設(shè)計軟件，該軟件集建模、有限元計算、承載力驗算、配筋設(shè)計、結(jié)果顯示及輸出等功能于一體，避免傳統(tǒng)方法的不足，顯著提升隧道襯砌設(shè)計效率。

2 軟件構(gòu)建思路及功能

如圖1所示，軟件的構(gòu)建思路主要分為6個步驟，具體如下。

圖1 鐵路隧道襯砌極限狀態(tài)法設(shè)計軟件構(gòu)建思路

(1)參數(shù)輸入。輸入?yún)?shù)類型包括模型幾何參數(shù)、襯砌材料參數(shù)、荷載參數(shù)、邊界參數(shù)及系數(shù)等。

(2)有限元計算。有限元計算得到襯砌單元的內(nèi)力及節(jié)點位移等。

(3)素混凝土驗算。對素混凝土襯砌進(jìn)行承載力驗算，判斷素混凝土承載力是否滿足要求。

(4)當(dāng)承載力滿足要求時，可直接顯示計算結(jié)果，包括襯砌內(nèi)力、位移、抗力狀態(tài)等結(jié)果顯示；若不滿足，則需進(jìn)行配筋設(shè)計。

(5)配筋設(shè)計。輸入鋼筋參數(shù)及配筋控制標(biāo)準(zhǔn)后，計算配筋量。

(6)結(jié)果輸出。輸出內(nèi)力計算結(jié)果及配筋設(shè)計結(jié)果等。

同樣從圖1可看出，該軟件具備的功能包括：①建模，包括襯砌模型建立、邊界條件建立及荷載施加；②有限元計算；③承載力驗算；④配筋設(shè)計；⑤結(jié)果顯示及輸出。

3 軟件模塊

3.1 前處理模塊

前處理模塊的主要任務(wù)為完成有限元計算所需參數(shù)的輸入。前處理模塊包括建模、施加邊界、施加荷載3個子模塊。其中，建?？赏ㄟ^Nodefile函數(shù)和Nodefilecad函數(shù)分別導(dǎo)入存有節(jié)點信息和單元信息的文本文檔以及存有隧道模型輪廓的dxf格式文件。

節(jié)點信息格式為：n，kn，x，y。其中，kn為節(jié)點編號，x為節(jié)點橫坐標(biāo)，y為節(jié)點縱坐標(biāo)。

單元信息格式為：e，ke，n1，n2，h。其中，ke為單元編號，n1為單元起點節(jié)點號，n2為單元終點節(jié)點號，h為單元厚度。

模型導(dǎo)入后，對襯砌的彈性模量、強度、泊松比等參數(shù)進(jìn)行賦值，完成模型建立。

根據(jù)二次襯砌的受力情況，邊界類型分為彈性支承及固定約束。在前處理模塊中，完成彈性支承邊界彈性反力系數(shù)的賦值，有固定約束時確定固定約束的節(jié)點及約束方向。另外，需根據(jù)用戶選擇的荷載類型和輸入的荷載參數(shù)，通過相應(yīng)的荷載計算函數(shù)分別計算荷載大小，為有限元計算模塊中的節(jié)點荷載計算奠定基礎(chǔ)。由于目前極限狀態(tài)法轉(zhuǎn)軌僅考慮了鐵路隧道襯砌及明洞結(jié)構(gòu)的自重及深埋、淺埋、偏壓圍巖荷載，Q/CR 9129—2018《鐵路隧道設(shè)計規(guī)范(極限狀態(tài)法)》中的襯砌目標(biāo)可靠指標(biāo)及分項系數(shù)也是僅在考慮了襯砌自重及圍巖荷載的基礎(chǔ)上提出的，而襯砌承受水壓力、地震荷載等特殊荷載時的極限狀態(tài)設(shè)計方法還有待進(jìn)一步研究。因此，本軟件在施加水壓荷載、地震荷載及自定義荷載時，由于分項系數(shù)暫不明確，用戶需自定義此類荷載的分項系數(shù)，而施加深埋、淺埋、偏壓及明洞荷載此類一般荷載時，軟件可自動按照Q/CR 9129—2018《鐵路隧道設(shè)計規(guī)范(極限狀態(tài)法)》的規(guī)定計算荷載大小并確定分項系數(shù)。

3.2 有限元計算模塊

二襯的有限元計算通過main函數(shù)完成。main函數(shù)包括幾個子函數(shù)，首先，通過函數(shù)input()將襯砌節(jié)點、單元及材料參數(shù)導(dǎo)入主函數(shù)；然后，通過stiff(ie)函數(shù)組建單元剛度矩陣，并通過locat(ie)函數(shù)形成單元定位向量后通過wstiff()函數(shù)組建襯砌總剛度矩陣；最后，通過fload()函數(shù)和bound()函數(shù)分別引入節(jié)點荷載和邊界條件后，通過gauss()函數(shù)求解有限元方程組。有限元計算實現(xiàn)方法如圖2所示。

圖2 有限元計算實現(xiàn)方法

3.3 配筋設(shè)計模塊

傳統(tǒng)配筋方法是通過假定襯砌截面鋼筋面積來試算最大裂縫寬度和承載力，將試算值與限制值對比，逐步增大鋼筋面積，最終滿足承載力及裂縫限制要求。當(dāng)最大裂縫寬度控制襯砌截面鋼筋面積時，采用傳統(tǒng)方法計算可能需要多次反復(fù)計算得到鋼筋面積，計算過程較為繁瑣。而按照圖3中的流程，將最大裂縫寬度限值作為已知參數(shù)來計算鋼筋混凝土襯砌截面鋼筋面積，僅需進(jìn)行一個流程的計算，便可得到同時滿足承載力和最大裂縫寬度要求的截面設(shè)計鋼筋面積，有效避免了反復(fù)調(diào)整鋼筋面積試算最大裂縫寬度的過程，提高了鋼筋混凝土襯砌截面鋼筋面積設(shè)計效率。

圖3 配筋計算流程

圖3中各參數(shù)計算公式如式(1)～式(5)所示。

(1)

ρte=As/Ate

(2)

ψ′=1.1-65ftk/σsq

(3)

(4)

(5)

式中，ωmax為設(shè)計最大裂縫寬度；As為受拉區(qū)縱向鋼筋截面面積；αcr為構(gòu)件受力特征系數(shù)；cs為最外層縱向受拉鋼筋外邊緣到受拉底邊的距離；deq為受拉縱向鋼筋等效直徑；Nq為按荷載標(biāo)準(zhǔn)或準(zhǔn)永久組合計算的軸向力值；Es為鋼筋的彈性模量；ρte為按有效受拉混凝土截面面積計算的縱向受拉鋼筋配筋率；ftk為混凝土軸心抗拉強度標(biāo)準(zhǔn)值；σsq為按荷載標(biāo)準(zhǔn)或準(zhǔn)永久組合計算的鋼筋混凝土構(gòu)件縱向受拉鋼筋應(yīng)力；ψ為縱向受拉鋼筋應(yīng)變不均勻系數(shù)；Ate為有效受拉混凝土截面面積；e為軸向壓力作用點至縱向受拉鋼筋合力點的距離；z為縱向受拉鋼筋合力點至截面受壓區(qū)合力點的距離。

式中各參數(shù)的取值方法參考Q/CR 9129—2018《鐵路隧道設(shè)計規(guī)范(極限狀態(tài)法)》。

3.4 后處理模塊

后處理模塊主要包括結(jié)果顯示模塊和結(jié)果輸出模塊。結(jié)果顯示模塊通過繪圖函數(shù)draw()完成，其中，draw()函數(shù)又包括6個子函數(shù)，分別負(fù)責(zé)襯砌彎矩、軸力、剪力、位移、抗力狀態(tài)、彈簧軸力的顯示；結(jié)果輸出模塊主要包含函數(shù)output()，通過將計算參數(shù)及計算結(jié)果寫入xlsx、doc和dwg格式的模板中，完成素混凝土結(jié)算結(jié)果、鋼筋混凝土計算結(jié)果、計算報告和鋼筋圖的輸出。

4 應(yīng)用案例

以時速250 km的IV級深埋雙線隧道為例，介紹軟件襯砌極限狀態(tài)法設(shè)計。首先，通過導(dǎo)入文本文檔或dxf文件建立襯砌模型，并輸入襯砌材料參數(shù)，所建模型如圖4所示。

圖4 襯砌模型

襯砌建立完成后，選擇荷載類型為深埋圍巖荷載，如圖5所示。定義深埋圍巖荷載時，首先選擇圍巖級別，軟件會根據(jù)圍巖級別規(guī)范自動輸入各參數(shù)的規(guī)范建議值，依次點擊“計算荷載”“施加荷載”后即可完成荷載的定義。邊界條件為沿襯砌全環(huán)設(shè)置彈性支承，計算模型如圖6所示。

圖5 深埋荷載計算及施加

圖6 計算模型

計算模型定義完成后，定義荷載分項系數(shù)及結(jié)構(gòu)抗裂和抗壓調(diào)整系數(shù)，如圖7所示。設(shè)置完成后即可進(jìn)行有限元計算。

圖7 系數(shù)設(shè)置

計算完成后，軟件可顯示襯砌的彎矩、軸力、剪力、位移、接觸狀態(tài)、彈簧軸力等分布圖，其中，襯砌彎矩及軸力的分布云圖分別如圖8、圖9所示。同時軟件對襯砌進(jìn)行素混凝土承載力驗算，驗算結(jié)果如圖10所示，驗算不通過需進(jìn)行配筋設(shè)計，如圖11所示。配筋設(shè)計完成后，輸出計算結(jié)果，輸出結(jié)果可包括素混凝土驗算結(jié)果、鋼筋混凝土驗算結(jié)果、配筋結(jié)果及計算報告等。

圖8 襯砌彎矩云圖(單位：N·m)

圖9 襯砌軸力云圖(單位：N)

圖10 素混凝土驗算結(jié)果

圖11 配筋設(shè)計

目前，該軟件在浦梅線等多條線路的極限狀態(tài)試設(shè)計項目中廣泛使用。大量應(yīng)用實踐表明，該軟件的計算結(jié)果與Ansys軟件計算結(jié)果基本相同，襯砌內(nèi)力的最大偏差在2%以內(nèi)，且完成一次計算僅需5 min左右，而采用多軟件協(xié)同工作的傳統(tǒng)方法完成一次計算一般需20 min左右，由此可見，該軟件在計算效率上有明顯優(yōu)勢。

5 應(yīng)用成果

在浦梅線的筆架山隧道、岐山隧道等工點的試設(shè)計項目中，采用該軟件分別根據(jù)基于極限狀態(tài)法的Q/CR 9129—2018《鐵路隧道設(shè)計規(guī)范(極限狀態(tài)法)》和基于安全系數(shù)法的TB 1003—2016《鐵路隧道設(shè)計規(guī)范》對鋼筋混凝土襯砌、素混凝土襯砌及明洞結(jié)構(gòu)進(jìn)行了平行設(shè)計。從設(shè)計結(jié)果上來看，兩種方法的設(shè)計成果基本相同，少數(shù)工況存在較小差別，部分設(shè)計結(jié)果如表1所示。

表1 極限狀態(tài)法和安全系數(shù)法設(shè)計結(jié)果對比

表2是筆架山隧道和岐山隧道分別采用兩種方法時的經(jīng)濟(jì)性對比。從經(jīng)濟(jì)性分析結(jié)果來看，筆架山隧道和岐山隧道采用極限狀態(tài)法與安全系數(shù)法試設(shè)計工程投資相當(dāng)，筆架山隧道極限狀態(tài)法較安全系數(shù)法工程投資減少15.410萬元，節(jié)約0.73%；岐山隧道極限狀態(tài)法較安全系數(shù)法工程投資減少14.91萬元，節(jié)約1.6%。

表2 筆架山隧道投資對比

6 結(jié)論

介紹基于VB6.0開發(fā)的鐵路隧道襯砌極限狀態(tài)法設(shè)計軟件的構(gòu)建思路、軟件功能及軟件模塊，并將軟件應(yīng)用于鐵路隧道極限狀態(tài)試設(shè)計中。該軟件集建模、有限元計算、承載力驗算、配筋設(shè)計、結(jié)果顯示及輸出等功能于一體，避免了多軟件共同作業(yè)的繁瑣操作。

目前，軟件已在鐵路隧道極限狀態(tài)法試設(shè)計中廣泛應(yīng)用，并在試設(shè)計中采用該軟件對鐵路隧道襯砌進(jìn)行了基于極限狀態(tài)法和安全系數(shù)的平行設(shè)計，得到以下結(jié)論。

(1)該軟件將隧道襯砌及明洞的單次設(shè)計時間由傳統(tǒng)方法所需的20 min縮短為5 min左右，大幅度提高了設(shè)計效率。

(2)Q/CR 9129—2018《鐵路隧道設(shè)計規(guī)范(極限狀態(tài)法)》和TB 1003-2016《鐵路隧道設(shè)計規(guī)范》對素混凝土襯砌、鋼筋混凝土襯砌及明洞結(jié)構(gòu)的設(shè)計結(jié)果基本相同，工程投資大致相當(dāng)。

該軟件不僅在鐵路隧道設(shè)計方法轉(zhuǎn)軌工作中發(fā)揮了重要作用，同時也為鐵路隧道極限狀態(tài)設(shè)計體系在實際工程中不斷修正和完善建立了一個可持續(xù)平臺。