陳新佳

(遼寧省水利水電勘測設(shè)計研究院有限責(zé)任公司，遼寧 沈陽 110006)

隨著國內(nèi)城鎮(zhèn)規(guī)模的不斷擴大，引調(diào)水工程項目發(fā)展非常迅速。泵站作為此類工程中的重要組成部分，在保證水資源的合理調(diào)度和管理方面起著不可替代的作用。本工程位于凌源市應(yīng)急供水工程四級泵站中的頭部位置，頭部泵站進水前池的上游為青龍河，河流多年水位變化不大，流量較穩(wěn)定，泵房基礎(chǔ)下部地質(zhì)主要分兩層：上層為砂礫石，下層的全風(fēng)化巖較深，地基承載力驗算后滿足設(shè)計要求。泵站基礎(chǔ)設(shè)計為筏型基礎(chǔ)結(jié)構(gòu)，選擇該類型的主要優(yōu)點有：結(jié)構(gòu)整體性好，適用于各種地基條件，自身重量大，抗浮和抗滑穩(wěn)定性較好[1]。

1 荷載組合及模型建立

1.1 郭氏法計算模型及荷載組合

頭部泵房基礎(chǔ)尺寸為48.89m×12.94m，厚度為1.5m，長度方向設(shè)有沉降縫。模型按照平面變形問題考慮，即地基是半無限空間體，基礎(chǔ)結(jié)構(gòu)在短跨方向各截面形狀和尺寸相同，同時荷載也是均勻分布的。據(jù)此假定建立地基梁模型，主要荷載考慮地上、地下部結(jié)構(gòu)在基礎(chǔ)邊緣產(chǎn)生的軸向力、彎矩，考慮單臺設(shè)備均布荷載，考慮側(cè)向土壓力在基礎(chǔ)邊緣產(chǎn)生的彎矩。根據(jù)設(shè)計報告中地震基本烈度為Ⅵ度，未考慮地震荷載；根據(jù)地質(zhì)報告中地下水埋深較深情況，未考慮泵房下部結(jié)構(gòu)受到地下水帶來的浮力及揚壓力，工況及荷載組合見表1。

對表1中各工況下的荷載組合進行比較后，確定最高水位運行期工況(前池水位最低)為最不利工況，即設(shè)計計算的控制工況。

表1 青龍河頭部泵站各工況荷載組合表

圖1 頭部泵房剖面圖

1.2 有限元法計算模型及荷載組合

頭部泵房的有限元模型取跨度為兩軸線間距12m寬作為計算對象，寬度內(nèi)取整個泵房上、下部結(jié)構(gòu)作為計算對象，地基范圍取為2倍泵房模型高度，平面方向延伸2倍泵房模型跨度[2]；考慮基礎(chǔ)、結(jié)構(gòu)和地基的相互作用，上部房頂結(jié)構(gòu)作為荷載傳遞至立柱上。主泵房共兩層，分為地上及地下部分。其中上部結(jié)構(gòu)(包括屋面、吊車梁及排架)總高9.6m，排架柱中間為填充磚墻；下部結(jié)構(gòu)總高10.2m，為鋼筋混凝土結(jié)構(gòu)，排架柱及下部結(jié)構(gòu)側(cè)墻與基礎(chǔ)底板整體澆筑，剖面如圖1所示。

2 郭氏彈性地基梁法的泵房基礎(chǔ)內(nèi)力

2.1 郭氏法地基反力求解方法

應(yīng)用郭氏法求解地基對基礎(chǔ)的反力時，一般先根據(jù)實際情況假定地基反力分布規(guī)律表達式，基于半無限體模型的郭氏法將梁下任意一點的地基反力p(x)近似地表示為有限項冪級數(shù)，即

p(x)=a0+a1x+a2x2+…+anxn

(1)

式中，ai(i=0，1，2，…，n)—待定的常系數(shù)。

再利用以下兩個方程建立包含ai常數(shù)的線性方程最終求得。

(1)基礎(chǔ)梁在荷載和反力共同作用下，撓度微分方程為

(2)

式中，E，I—基礎(chǔ)梁的彈性模量、截面慣性矩；q(x)—梁上荷載。

平面問題中，半無限大彈性體由彈性理論導(dǎo)出的計算公式為

隰縣作為一個國家扶貧開發(fā)工作重點縣，脫貧攻堅進入沖刺階段，農(nóng)村安全飲水工程也不斷加碼。飲水安全工程不僅為老百姓送來了健康，也讓不少農(nóng)民借“水”生財。

(3)

式中，P—在半無限大彈性體邊界上的集中力；s—任意點至P的距離；r—要求其沉降的點至P的距離；E0，u0—地基的壓縮模量及泊松比。

根據(jù)基礎(chǔ)梁與地基之間不會出現(xiàn)拉力，在荷載作用下，梁與地基緊密接觸，即變形連續(xù)假定，得到變形連續(xù)方程y(x)=W(x)，求得各個常數(shù)，再代入式(1)求得地基反力p(x)，最終利用截面法求得梁的內(nèi)力。當(dāng)郭氏法中p(x)所取級數(shù)項數(shù)越多，結(jié)果的準(zhǔn)確性越高[3]。

2.2 郭氏法查表計算結(jié)果

查表前，需對地基梁按照柔性指數(shù)t對梁進行分類，計算公式如下：

(4)

式中，l—梁的一半長，h—梁的厚度。

根據(jù)柔性指數(shù)取值范圍可以將基礎(chǔ)梁分為3類：剛性梁(t<1)、有限長梁(集中荷載1≤t≤10，均布荷載1≤t≤50)、無限長梁(集中荷載t﹥10，均布荷載t>50)。其中前兩類按照短梁查表計算，第三類按照長梁查表計算。地基壓實模量E為4.47MPa；梁彈性模量為2.8×106t/m。

本文根據(jù)公式計算后得到底板地基梁柔性指數(shù)t為0.12<1，應(yīng)按照短梁查表。采用1.1節(jié)中最不利荷載組合計算得到的結(jié)果如下：

屋面及吊車荷載：彎矩(M)為±107.5kN·m，向下軸力(↓)為50kN；

上部磚墻結(jié)構(gòu)荷載：向下軸力(↓)為78kN；

下部混凝土結(jié)構(gòu)荷載：向下軸力(↓)為175kN；

單臺設(shè)備荷載：向下軸力(↓)為58kN；

下部混凝土承受側(cè)向土壓力荷載：彎矩(→)為570kN·m；

下部混凝土承受前池最低水位水壓力荷載：彎矩(←)87kN·m。

通過查表得到地基梁分段彎矩控制值(→)為751kN·m，剪力值(↓)為263kN；再由材料力學(xué)應(yīng)力求解公式得到地基梁分段最大應(yīng)力σmax為1.422MPa，在地基梁端部與前池連接位置。

3 基于有限單元法的泵房整體三維模擬

3.1 有限單元法求解過程分析

采用SolidWorks進行三維實體建模，并通過與有限元分析軟件提供的通用格式對泵房進行受力模擬，其中混凝土采用SOLID 65單元模擬，軟件模擬的8節(jié)點3D實體單元能實現(xiàn)開裂、壓碎及溫度升降應(yīng)力、應(yīng)力釋放等特性?；炷帘緲?gòu)關(guān)系模型采用軟件默認的理想彈塑性本構(gòu)模型，同時采用非線性拉應(yīng)力準(zhǔn)則及混凝土多參數(shù)破壞準(zhǔn)則[4]。

主泵房采用上、下部結(jié)構(gòu)與泵房底板整體建模方式，軟件計算分析時選取其中一跨12m進行模擬(同郭氏法)，整個有限元模型8節(jié)點六面體等參單元共有27256個，節(jié)點36439個。三維有限元模型網(wǎng)格如圖2所示。

圖2 三維有限元模型網(wǎng)格劃分圖

泵房荷載包括傳至牛腿柱的屋面荷載、吊車荷載采用集中力的加載方式；底板上設(shè)備荷載采用面荷載的加載方式；兩側(cè)混凝土墻體上的土壓力及前池側(cè)向水壓力采用應(yīng)力梯度以面荷載的方式施加(工況及主要荷載組合同表1)。由于底板的剛度較大，故可以看作埋置于土中的剛性基礎(chǔ)，將地基底部所有節(jié)點進行Y方向的自由度約束，而縱向墻體邊界上的節(jié)點進行X和Y兩個方向的自由度約束[5]，依據(jù)此約束計算得到的泵房進水側(cè)在控制工況下(最高揚程)的最大應(yīng)力。

3.2 泵房整體穩(wěn)定有限元計算結(jié)果

根據(jù)有限元原理中協(xié)調(diào)關(guān)系與本構(gòu)關(guān)系原理及平衡方程，應(yīng)用軟件對模型進行分析，計算參數(shù)、工況及荷載組合在1.2節(jié)中已經(jīng)介紹，求解采用最不利工況下節(jié)點的位移、應(yīng)變和應(yīng)力，以便同郭氏彈性地基梁法進行對比。在此基礎(chǔ)上，應(yīng)用有限元法三維整體建模，不僅可以得到基礎(chǔ)內(nèi)力分布圖，亦可得到泵房上、下部結(jié)構(gòu)的整體內(nèi)分布圖，軟件計算得到的模型整體結(jié)構(gòu)應(yīng)力云圖如圖3—5所示，X、Y方向應(yīng)力云圖如圖4—5所示。

圖3 第一主應(yīng)力云圖

圖4 X方向應(yīng)力云圖

圖5 Y方向應(yīng)力云圖

由軟件分析結(jié)果可以看出，泵房底板混凝土最大拉應(yīng)力位于底板邊緣下側(cè)與前池連接位置，拉應(yīng)力值為1.398MPa，此值與郭氏彈性地基梁法計算結(jié)果較接近，說明模型假定的邊界條件符合實際情況。

此外，圖4—5中泵房底板上部電機和水泵壓應(yīng)力值為0.567MPa，分析時還需考慮事故工況水錘壓力，故推測在設(shè)備底部中心底板所受應(yīng)力較大(具體值在打壓試驗后通過相關(guān)數(shù)據(jù)進行模擬)，應(yīng)加強配筋；分析結(jié)果中排架中部分牛腿以下最大拉應(yīng)力為1.1～1.5MPa，位置在柱與梁、墻相交處，此處也需加強配筋。而混凝土墻體內(nèi)的下柱段應(yīng)力較小，配筋時可與泵房水下墻一并考慮，亦符合以往的工程設(shè)計經(jīng)驗。

4 結(jié)論

(1)郭氏法在計算假定半無限彈性地基模型時，反映了地基的連續(xù)整體性，計算模型相對簡單，可通過編程提高方程冪級數(shù)，從而提高計算精度，對于厚土層或整體巖石基礎(chǔ)的求解較適用[6]。但模型無法反映地基土非彈性性質(zhì)和分層復(fù)雜性，因此在應(yīng)用時需注意所選場地的地質(zhì)復(fù)雜程度。

(2)本文應(yīng)用郭氏法對比有限元軟件進行模擬分析，得到結(jié)果基本相同，復(fù)核了有限元法在建模及參數(shù)、邊界條件設(shè)置時的正確性。分析有限元計算最大應(yīng)力值略小于郭氏法的最大應(yīng)力值，這是建模時將泵房底板與排架柱、水下墻整體考慮的結(jié)果。

本文將兩種方法應(yīng)用于具體案例，計算彎矩、剪力值符合工程實際，在泵房、水閘等筏型基礎(chǔ)的結(jié)構(gòu)設(shè)計[7- 10]、配筋計算中具有一定參考價值。