宋力，吳明晗

（1.江蘇省水利科學(xué)研究院材料結(jié)構(gòu)研究室，江蘇 揚州 225002；2.江蘇省水利工程科技咨詢有限公司，江蘇 南京 210000）

1 研究背景

閉口框架結(jié)構(gòu)是下穿立交工程中常用的結(jié)構(gòu)形式，其結(jié)構(gòu)內(nèi)力計算因邊界約束條件定義的不同會出現(xiàn)很大差別。在很多類似結(jié)構(gòu)的設(shè)計中，常將底板厚度設(shè)計較厚，而頂板較薄，其設(shè)計的理論依據(jù)正是底板的內(nèi)力較頂板大。本文針對這一問題采用基于位移法的反力自平衡法和彈性地基梁法分別進(jìn)行計算，通過對比分析頂、底板的內(nèi)力分布情況，為結(jié)構(gòu)設(shè)計提供依據(jù)。

2 力學(xué)模型的建立

本文選取的閉口框架結(jié)構(gòu)頂、底板計算長度15.9m，側(cè)墻計算長度6m，頂板覆土2.8m；頂板厚1.2m，底板厚1.3m，側(cè)墻厚1m。結(jié)構(gòu)荷載包括：結(jié)構(gòu)自重G1、覆土重量G2、側(cè)墻土壓力F1?？蚣芙Y(jié)構(gòu)容重25kN/m3，覆土容重19kN/m3。

3 內(nèi)力分析

目前，在閉口框架結(jié)構(gòu)內(nèi)力計算過程中，常用的方法有兩種：一種為基于位移法的反力自平衡法；另一種為彈性地基梁法。

文章以有限元軟件MIDAS為基礎(chǔ)，沿結(jié)構(gòu)縱向單位長度（1m）建立計算模型。

3.1 反力自平衡法計算分析

梁單元內(nèi)力計算理論基礎(chǔ)為梁單元的平衡方程：

圖1 閉口框架結(jié)構(gòu)簡化模型

圖2 基于反力自平衡法計算彎矩圖

圖3 基于彈性地基梁法計算彎矩圖

上式作為梁單元有限元計算程序的理論基礎(chǔ)，對梁單元離散化，在外荷載q、P 的作用下，求解各節(jié)點位移w，進(jìn)一步求得梁單元內(nèi)力。

反力自平衡法是在有限元計算的基礎(chǔ)上，假定對閉口框架結(jié)構(gòu)施加初始邊界條件（鉸支座及鏈桿），在計算過程中消除鉸支座及鏈桿的支座反力，使結(jié)構(gòu)自平衡，其基本思路如下：

①結(jié)構(gòu)自重G1單獨作用于結(jié)構(gòu)，計算豎向支反力RY；

②在鉸支座及鏈桿支座節(jié)點處RY之和折算為底板線荷載，反作用于底板，模擬基底反力，重新計算，使支反力RY=0；

③施加頂板覆土重量G2（按線荷載計），同時底板施加大小相等、方向相反的基底反力；

④施加側(cè)墻土壓力。

運行計算程序，計算結(jié)構(gòu)內(nèi)力。

3.2 基于彈性地基梁法計算分析

建立彈性地基梁單元的偏微分方程：

式中k 為地基彈性支承剛度，其值與地基土水平抗力的比例系數(shù)m 相關(guān)。

文章針對地基土水平抗力的比例系數(shù)m 分別選取1000kN/m4、10000kN/m4、20000kN/m4、30000kN/m4進(jìn)行比較。底板單元長度a=0.53m，縱向計算長度為s=1m，底板底部埋深h=10m，底板底部由土體提供的彈性支承剛度k 由下式計算：

四種m值對應(yīng)的彈性支承剛度k 分別為5300kN/m、53000kN/m、106000kN/m。由于土體不提供拉力約束，因此本文采用只受壓彈簧單元定義土體對結(jié)構(gòu)底板的約束。

3.3 計算成果對比分析

彎矩計算成果匯總表(單位：kN·m)

從上表可以看出，反力自平衡法計算結(jié)果與m 值無關(guān)，即假定地基對結(jié)構(gòu)不產(chǎn)生約束（m=0），地基僅對結(jié)構(gòu)產(chǎn)生反作用（地基反力）；彈性地基梁法計算結(jié)果表明，m 值的選取對內(nèi)力分布及框架結(jié)構(gòu)的剛度分配產(chǎn)生較大影響，m 趨近于0時，內(nèi)力計算結(jié)果趨近于反力平衡法計算結(jié)果。由此可見，采用彈性地基梁法進(jìn)行計算時，m 值的合理選用十分重要，反力平衡法作為m=0時一個特例，其內(nèi)力計算結(jié)果是趨于保守的。

4 結(jié) 論

反力自平衡法將各種豎向荷載對應(yīng)的基底反力反作用于底板，實現(xiàn)結(jié)構(gòu)的自平衡（支反力為零），最初定義的邊界條件（鉸支座及鏈桿支座）對最終計算結(jié)果沒有影響，從而回避了因邊界條件與實際約束不符所帶來的影響，這是該方法的優(yōu)點；但其缺點是忽略了地基與結(jié)構(gòu)的協(xié)調(diào)變形對結(jié)構(gòu)剛度分配的有利作用，內(nèi)力結(jié)果過于保守，這在底板上的內(nèi)力分布表現(xiàn)的尤為明顯。

M法在底板、土體之間建立只受壓的彈簧單元，較好地反映了地基與底板之間的協(xié)調(diào)變形，邊界條件更為合理。從表中可以看出，在m 值較大時底板內(nèi)力較頂板內(nèi)力小很多，因此底板厚度并不需要較頂板厚，甚至可適當(dāng)降低，對于節(jié)約工程成本具有較好的經(jīng)濟(jì)意義。但缺點是m 值作為M法的一個關(guān)鍵指標(biāo)，不易取得精確值，實際應(yīng)用時可根據(jù)工程所處土層的狀態(tài)按下限取值。

[1]I·M·Smith,D·V·Griffiths.有限元方法編程（第三版）[M].王崧,等譯.北京:電子工業(yè)出版社,2003.

[2]王勖成.有限單元法[M].北京:清華大學(xué)出版社,2003.

[3]朱合華.地下建筑結(jié)構(gòu)[M].北京:中國建筑工業(yè)出版社,2006.

[4]劉建航,侯學(xué)淵.基坑工程手冊[M].北京:中國建筑工業(yè)出版社,2003.