預(yù)應(yīng)力混凝土收縮徐變空間計(jì)算程序

李君鳳

(同濟(jì)大學(xué)建筑設(shè)計(jì)研究院(集團(tuán))有限公司，上海 200092)

0 引言

大型通用軟件計(jì)算混凝土的收縮、徐變不但要將徐變本構(gòu)關(guān)系用金屬蠕變表示，收縮用溫度變化表示，而且計(jì)算并不準(zhǔn)確，一般很少用大型通用軟件計(jì)算混凝土收縮、徐變。MIDAS/Civil中計(jì)算混凝土的收縮徐變是根據(jù)定義材料的時(shí)間依存性實(shí)現(xiàn)的，應(yīng)該是如今計(jì)算混凝土的收縮、徐變方面最成熟的軟件了，它計(jì)算準(zhǔn)確，可以進(jìn)行平面和空間的計(jì)算。但是它沒有引入我國規(guī)范的收縮徐變模式。橋梁博士Dr．Bridge結(jié)合規(guī)范比較緊密，運(yùn)用方便，但是Dr．Bridge是平面分析軟件，對(duì)于稍稍復(fù)雜的需要用空間進(jìn)行分析的結(jié)構(gòu)無法計(jì)算。

針對(duì)以上問題，本程序添加了混凝土非線性分析模塊。根據(jù)規(guī)范中對(duì)混凝土收縮徐變的規(guī)定，利用空間有限單元逐步計(jì)算混凝土收縮、徐變效應(yīng)及預(yù)應(yīng)力損失效應(yīng)并將其形象的體現(xiàn)在圖形中。

1 程序中徐變、收縮的分析方法

假定在某一時(shí)刻混凝土的彈性模量E保持不變，用齡期調(diào)整的有效模量Eφ(ti，ti－1)代替混凝土的彈性模量E，則應(yīng)力與應(yīng)變?cè)隽康年P(guān)系為:

將其應(yīng)用于空間8節(jié)點(diǎn)等參單元，根據(jù)應(yīng)力與應(yīng)變之間的物理關(guān)系σ=Dε，將混凝土看作各向同性線彈性材料，利用最小位能原理建立有限元方程Ka=P時(shí)，在形成單元?jiǎng)偠染仃嘖e時(shí)，將中的 D 用 Dφ(ti，ti－1)代替即可。

在具體計(jì)算某一時(shí)刻段ti－ti－1時(shí)，根據(jù)有限單元法形成荷載矩陣的原理，對(duì)結(jié)構(gòu)中任一空間8節(jié)點(diǎn)等參單元施加約束，使在第titi－1個(gè)時(shí)間內(nèi)節(jié)點(diǎn)變位增量保持為0，則從式(1)可得到約束產(chǎn)生的單元應(yīng)力增量:

其中，

將單元應(yīng)力增量作為初應(yīng)力，運(yùn)用考慮單元存在初應(yīng)力和初應(yīng)變情況的單元等效節(jié)點(diǎn)荷載列陣一般表達(dá)式(4)，轉(zhuǎn)化為單元節(jié)點(diǎn)等效荷載增量ΔPe*，然后形成整體荷載列陣ΔP*，即ti時(shí)刻在整體荷載列陣ΔP*的作用下，混凝土由收縮徐變引起的整體位移為0，也就是各單元的單元應(yīng)變?yōu)?。那么，將整體荷載列陣ΔP*反號(hào)加在結(jié)構(gòu)上，根據(jù)有限元求解的一般方法，按照式(4)便得到在第ti－ti－1個(gè)時(shí)間內(nèi)由混凝土收縮、徐變?cè)斐傻慕Y(jié)構(gòu)變形δcs，同時(shí)還可以求出在第ti－ti－1個(gè)時(shí)間內(nèi)由混凝土收縮、徐變引起的單元應(yīng)變?cè)隽?/p>

然后將式(1)轉(zhuǎn)化成計(jì)算8節(jié)點(diǎn)有限元的式(6)，計(jì)算出由收縮、徐變引起的應(yīng)力增量。

在計(jì)算某一工況時(shí)，首先計(jì)算在該工況受力狀況下的結(jié)構(gòu)位移和應(yīng)力，然后計(jì)算該工況從開始到結(jié)束時(shí)間段內(nèi)由收縮徐變引起的應(yīng)力和位移增量，并存儲(chǔ)。計(jì)算收縮徐變時(shí)，需將工況從開始到結(jié)束的時(shí)間分成若干計(jì)算時(shí)間段，根據(jù)收縮、徐變所引起應(yīng)變?cè)隽侩S時(shí)間的變化關(guān)系，劃分計(jì)算時(shí)間段時(shí)，按照冪函數(shù)劃分。計(jì)算的具體步驟見圖1。

圖1 工況的計(jì)算過程

2 收縮徐變與預(yù)應(yīng)力損失

本程序中的預(yù)應(yīng)力是用等效荷載來考慮的，而等效荷載列陣在收縮、徐變引起的結(jié)構(gòu)變形中并沒有變化，但實(shí)際結(jié)構(gòu)中則產(chǎn)生預(yù)應(yīng)力損失。

由混凝土收縮、徐變引起的構(gòu)件預(yù)應(yīng)力鋼筋的預(yù)應(yīng)力損失屬于與結(jié)構(gòu)變形有關(guān)的預(yù)應(yīng)力損失。采用空間8節(jié)點(diǎn)等參單元，在計(jì)算單元?jiǎng)偠葧r(shí)，將預(yù)應(yīng)力提供的剛度加進(jìn)去，預(yù)應(yīng)力與結(jié)構(gòu)物一同參加有限元計(jì)算，這對(duì)于由混凝土彈性壓縮引起的預(yù)應(yīng)力損失是適合的。根據(jù)圖1可知，在鎖定結(jié)構(gòu)，即使單元節(jié)點(diǎn)的應(yīng)變?yōu)?，得到在收縮徐變作用下由約束產(chǎn)生的單元應(yīng)力增量并形成等效荷載列陣 ΔP*時(shí)，和預(yù)應(yīng)力鋼筋并沒有任何的關(guān)系。在式(3)的計(jì)算中，只需考慮單元中混凝土的彈性模量，預(yù)應(yīng)力鋼筋引起的增量不予考慮。但是在釋放約束后將ΔP*反號(hào)作用于結(jié)構(gòu)，求解結(jié)構(gòu)由收縮、徐變引起的位移增量δcs，單元應(yīng)變?cè)隽浚瑔卧獞?yīng)力增量及結(jié)構(gòu)約束力增量時(shí)，預(yù)應(yīng)力和混凝土則是同時(shí)工作，應(yīng)該考慮預(yù)應(yīng)力鋼筋的影響，在計(jì)算剛度矩陣時(shí)需要將預(yù)應(yīng)力鋼筋引起的剛度的增量加上去。這樣，由混凝土收縮、徐變引起的結(jié)構(gòu)變形導(dǎo)致的預(yù)應(yīng)力損失在有限元的計(jì)算中自動(dòng)扣除，無需額外計(jì)算。

3 工程實(shí)例

以一次落架的預(yù)應(yīng)力混凝土連續(xù)箱梁橋?yàn)槔涑叽缛鐖D2所示。

圖2 三跨連續(xù)梁立面圖（單位：m）

在建模時(shí)各工況的起始時(shí)間如表1所示。計(jì)算結(jié)果如圖3～圖5所示。

圖3 最后階段縱向應(yīng)力圖

圖4 收縮、徐變縱向次內(nèi)力應(yīng)力圖

圖6 最后階段位移圖

表1 連續(xù)梁各工況表

對(duì)于一次落架的連續(xù)梁，如果按照梁?jiǎn)卧?jì)算，混凝土徐變只導(dǎo)致結(jié)構(gòu)變形的增加并不引起次內(nèi)力。用本程序計(jì)算，由于在計(jì)算中計(jì)入了預(yù)應(yīng)力損失的影響，故計(jì)算結(jié)果中由收縮、徐變引起的次內(nèi)力并不為0。

4 結(jié)語

本文針對(duì)目前橋梁結(jié)構(gòu)非線性分析的現(xiàn)狀，按有效彈性模量的計(jì)算方法，編制了相應(yīng)的橋梁結(jié)構(gòu)分析程序。利用工程實(shí)例，驗(yàn)證了本程序計(jì)算的合理性和實(shí)用性。本文所述的自編程序，有著廣泛的工程應(yīng)用前景，但求解以及本程序與相關(guān)數(shù)據(jù)庫的連接接口等有待進(jìn)一步的深入研究。

［1] 中國公路建設(shè)市場(chǎng)編輯部．橋梁結(jié)構(gòu)分析程序及軟件的發(fā)展［J］．中國公路建設(shè)市場(chǎng)，2004(5):9-11．

［2] 陳明憲，彭建新，顏東煌，等．按齡期調(diào)整的有效彈性模量法分析混凝土收縮徐變［J］．長沙交通學(xué)院學(xué)報(bào)，2005，20(3):92-95．

［3] 李 毅．基于空間等參元并考慮預(yù)應(yīng)力效應(yīng)的橋梁分析程序開發(fā)［D］．上海:同濟(jì)大學(xué)碩士學(xué)位論文，2008．

［4] 李君鳳，吳 迅，李 毅．混凝土收縮徐變空間模擬［J］．山西建筑，2010，36(6):1-3．

［5] JTG D62-2004，公路鋼筋混凝土及預(yù)應(yīng)力混凝土橋涵設(shè)計(jì)規(guī)范［S］．