蘇志軍 李競 蘇志營
1 中化地質(zhì)礦山總局,北京 100013
2 中化地質(zhì)礦山總局河南地質(zhì)局,河南 鄭州 450002
3 鄭州中蘇巖土工程有限公司,河南 鄭州 450002
4 河南亞星建筑安裝工程有限公司,河南 鄭州 450012
鄭州地區(qū)空間規(guī)劃與工程建設(shè)日趨緊密,具有能夠承載多層建筑物且地層分布范圍大、特征明顯的巖土體(特征土),研究其特性合理化利用越來越為土木工程師所關(guān)注。持續(xù)發(fā)展的基礎(chǔ)建設(shè)使得工程條件越來越復(fù)雜,通過顆粒介質(zhì)材料承受并傳遞上部荷載,多數(shù)工程幾何形態(tài)多樣、力學(xué)環(huán)境復(fù)雜、地質(zhì)條件復(fù)雜,使得當(dāng)前在結(jié)構(gòu)力學(xué)、彈性力學(xué)、巖土力學(xué)等傳統(tǒng)方法中難以獲得解析解。而隨著計算機(jī)技術(shù)的發(fā)展,數(shù)值分析(數(shù)值模擬)方法已逐漸成為土木工程求解科學(xué)問題不可或缺的分析手段。
因此,采用數(shù)值模擬的方法以鄭州地區(qū)奧蘭花園住宅小區(qū)場地淺部特征土為主要對象,將原位測試、室內(nèi)試驗相結(jié)合來開展本次研究。
特征土是指土體內(nèi)存在可辨識包含物或特征明顯的巖土體,其物理力學(xué)性質(zhì)具規(guī)律性,在較大范圍內(nèi)都有分布;能夠承擔(dān)一定的上部荷載;測試特征曲線明顯;特征土層具有科研研究和經(jīng)濟(jì)利用價值等。鄭州淺部特征土一般是指層底埋深10~15m以上的具上述特征的典型巖土層。
奧蘭花園住宅小區(qū)場地位于鄭州市鄭東新區(qū),東風(fēng)東路與鄭汴路交匯處北約300m路西,北臨康寧街,南側(cè)為福祿街,東臨東風(fēng)東路,西傍康平路。場地地貌單元區(qū)域上屬于黃河泛濫沖積平原,場地內(nèi)勘探孔揭露85.0m范圍內(nèi),5.9~10.3m以上為新近紀(jì)堆積粉土、粉質(zhì)粘土;約10.0~30.0m為第四紀(jì)全新世沖積形成的地層,以粉土、細(xì)砂層為主;30.0~85.0m為第四紀(jì)晚更新世沖積形成的地層,以粉土、粉質(zhì)粘土為主[1]。
奧蘭花園分布特征土上部(層底埋深5.5~9.8m):地層呈灰-灰(褐)黃色,濕,中密。搖振反應(yīng)中等,干強(qiáng)度低,韌性低,無光澤。有粘性,質(zhì)軟,偶見少量蝸牛殼碎片。局部有粉質(zhì)粘土夾層,軟塑。
其下部(層底埋深6.4~12.6m):地層呈淺灰-灰(褐)黃色,濕,密實。搖振反應(yīng)迅速,干強(qiáng)度低,韌性低,無光澤。本層稍有砂感,局部發(fā)育粉砂薄層。
鄭州地區(qū)四至六層、高層低配裙樓及荷載較大的廠房等建(構(gòu))筑物大多直接或間接以地基處理形式采用該特征土層為主要持力層。
根據(jù)鄭州市城市人防工程地質(zhì)勘查資料[2]及東西南北不同區(qū)域結(jié)合項目勘探情況,搜集整理鄭州城區(qū)特征土分布情況見表1。
表1 鄭州城區(qū)特征土分布情況一覽表Table 1 Distribution table of characteristic soils in Zhengzhou urban area
土是由巖石經(jīng)過風(fēng)化后產(chǎn)生的松散物集合體,主要的特點是不連續(xù)性和三相性(土固體顆粒-固相,土中水-液相,土中氣-氣相),這使其變形和強(qiáng)度等力學(xué)性質(zhì)不同于其它連續(xù)介質(zhì),變得極其復(fù)雜。目前,采用較多的數(shù)值模擬方法是離散元法和有限差分法、有限元法等。
1.2.1 離散元法
離散元法克服了傳統(tǒng)連續(xù)介質(zhì)力學(xué)模型的宏觀連續(xù)性假設(shè),從細(xì)觀層面上對土的工程特性進(jìn)行數(shù)值模擬,并通過細(xì)觀參數(shù)的研究來分析宏觀力學(xué)行為,尤其適用于粉土、砂性土等具離散介質(zhì)的力學(xué)分析。由于其研究基本構(gòu)成為顆粒,因此從本質(zhì)上適合研究固體(松散或黏結(jié))介質(zhì)的力學(xué)特性[3]。
1.2.2 有限差分法與有限元法
有限差分法和有限元法是以連續(xù)介質(zhì)為基礎(chǔ)考慮,它們考慮了土體的非線性彈塑性關(guān)系,且都產(chǎn)生一組待解方程組。盡管這些方程是通過不同方式推導(dǎo)出來的,但兩者產(chǎn)生的方程式一致的。有限單元程序需要將單元矩陣組合成大型整體剛度矩陣,有限差分則無需如此。它們的缺點是不能體現(xiàn)顆粒間的復(fù)雜相互作用,不能刻畫類似粉土這類特征土的流動變形特征。
因此,側(cè)重研究與實際應(yīng)用,本文主要選擇離散元法進(jìn)行數(shù)值模擬。
奧蘭花園項目對原狀土樣進(jìn)行了常規(guī)項目的分析(表2),用四聯(lián)電動剪力儀做直剪試驗;靜三軸做不固結(jié)不排水UU剪切試驗(表3);同時進(jìn)行了高壓固結(jié)試驗測試(表4)。
表2 特征土物理性質(zhì)統(tǒng)計表Table 2 Statistical table of physical properties of characteristic soil
表3 特征土力學(xué)性質(zhì)統(tǒng)計表Table 3 Statistical table of mechanical properties of characteristic soils
表4 特征土高壓固結(jié)成果統(tǒng)計表Table 4 Statistical table of high pressure consolidation results of characteristic soils
2.2.1 靜力觸探
靜力觸探采用LT-20A靜力觸探工程車施工,雙橋探頭,配備LMC-310C型自動采集、記錄、數(shù)據(jù)處理打印微機(jī)系統(tǒng)[4]。
2.2.2 標(biāo)準(zhǔn)貫入
標(biāo)準(zhǔn)貫入試驗采用標(biāo)準(zhǔn)貫入儀作現(xiàn)場測試,嚴(yán)格按規(guī)范執(zhí)行,試驗時清孔干凈,扶正鉆桿,詳細(xì)記錄貫入度。采用導(dǎo)向桿變徑自動脫鉤式落錘裝置(錘重為63.5kg,落距為76cm)配合鉆機(jī)進(jìn)行測試。
2.2.3 測試結(jié)果統(tǒng)計
特征土進(jìn)行原位測試數(shù)理統(tǒng)計(表5)。
表5 特征土原位測試統(tǒng)計表Table 5 Statistical table of in-situ test of characteristic soil
從表中可以看出,特征土上部略軟下部偏硬,數(shù)理統(tǒng)計符合規(guī)范要求,經(jīng)測算地基承載力上下相差不大,擬建建筑物存在負(fù)荷小的多層建筑物時可考慮上部承載;存在負(fù)荷大的多層建筑物時可考慮下部為主要持力層。
直剪試驗(直接剪切試驗)是將制備的代表性試樣放入剪切盒內(nèi),將上盒固定,施加垂直壓力,逐級施加水平剪力使下盒沿水平方向滑動,直至試樣剪切破壞。因而對于直剪試驗的顆粒流模型,可僅建立直剪儀剪切盒與盒內(nèi)的土樣。通過定義,使剪切盒發(fā)生一定量的位移并使土樣產(chǎn)生剪切破壞[5]。
實際顆粒流模擬時主要分為試樣生成,圍壓施加,加載三個步驟實施。
直剪試驗測試相對簡單,試驗不能嚴(yán)格控制排水條件;剪切面不是沿最薄弱面產(chǎn)生破壞,土樣的剪切破壞始于邊緣,應(yīng)力集中發(fā)生在邊緣(圖1)。
圖1 PFC顆粒流直剪試驗?zāi)MFig.1 Direct shear test simulation of PFC particle flow
縱觀剪切模擬,細(xì)粒土存在明顯的剪切帶,這在直剪試驗過程中不易直接觀察到,但試樣破壞后模擬與試驗形態(tài)呈現(xiàn)較好的一致性。破壞形式為受限強(qiáng)制性剪切破壞。
三軸試驗(三軸壓縮試驗)是將制備的代表性試樣放置于充滿液體的壓力室內(nèi),試樣用薄橡皮膜包裹,使土樣與膜外液體完全隔開。在給定的三軸壓力室周圍壓力作用下,不斷加大軸向附加壓力,直至試樣破壞。當(dāng)周圍壓力(σ1為大主應(yīng)力,σ2為中主應(yīng)力,σ3為小主應(yīng)力,σ2、σ3為圍壓)為不同方向大?。é?≠σ3)時,為真三軸試驗;當(dāng)三軸不同方向周圍壓力均為同一常數(shù)(σ2=σ3=k)時,為靜三軸試驗。
Robert D. Holtz指出:“莫爾早在1900年就假設(shè)了一個實際材料的破壞準(zhǔn)則:當(dāng)破壞面的剪應(yīng)力達(dá)到了該面上的某個單值函數(shù)時,材料就破壞了[6]?!毕鄬τ讦?=σ3=k的情況,中主應(yīng)力的增加會提高土的強(qiáng)度已經(jīng)成為人們的共識。
實驗表明,主應(yīng)力方向?qū)辜魪?qiáng)度的影響大于中主應(yīng)力的影響,從實用的角度出發(fā),忽略中主應(yīng)力的影響是偏于安全的[7]。綜上所述,本次PFC模擬主要是靜三軸試驗(圖2),不同粒徑下剪切帶分布及取向不同,按試樣級配顆粒組成模擬模型進(jìn)行模擬計算,模擬破壞形態(tài)與試驗取得較好的一致性。
圖2 PFC顆粒流三軸試驗?zāi)MFig.2 Triaxial test simulation of PFC particle flow
根據(jù)粗粒料真三軸試驗結(jié)果結(jié)合顆粒流數(shù)值模擬,河海大學(xué)巖土力學(xué)與堤壩工程教育部重點試驗室得出真三軸建立的各向同性強(qiáng)度準(zhǔn)則[8],結(jié)合鄭州地區(qū)特征土不同區(qū)域取土測試結(jié)果及數(shù)值模擬,建立靜三軸特征土破壞方程式(1):
式中:
φ0為靜三軸壓縮條件下內(nèi)摩擦角
根據(jù)特征土分布情況,采用集管理、專業(yè)及系統(tǒng)分析軟件為一體的信息化模式,以提高試驗?zāi)M與實際工程相關(guān)參數(shù)的準(zhǔn)確度[9]。取特征土土層上下部視為一體考慮,采樣六組(S-1—S-6)進(jìn)行靜三軸不固結(jié)不排水試驗(UU)——將制備好的代表性試樣放入三軸儀壓力室內(nèi),在排水閥門關(guān)閉的情況下,先向土樣施加周圍壓力,隨即施加軸向應(yīng)力進(jìn)行剪切,直至破壞的試驗,試驗成果整理成不固結(jié)不排水強(qiáng)度包線圖(圖3),將試驗結(jié)果與破壞方程反演計算結(jié)果進(jìn)行對比。
依照摩爾-庫倫理論破壞準(zhǔn)則,可以采用最大主應(yīng)力比較法、最小主應(yīng)力比較法以及內(nèi)摩擦角比較法來判別土體單元是否發(fā)生破壞[10-11]。建立的特征土破壞方程反演參數(shù),可以方便進(jìn)行這種判斷。
圖3 三軸(不固結(jié)不排水)試驗強(qiáng)度包線Fig.3 Triaxial (unconsolidated and undrained) test strength enveloping line
表6 試驗值與反演值的比較Table6 Comparison of test value and inversion value
從表6中可以看出,由于土質(zhì)不均勻性,代表性試樣測試結(jié)果存在一定離散性,而采用破壞方程反演計算的結(jié)果一致性比較好。這為技術(shù)人員后期分析提供了便利,也為采樣困難、地層在采取一定質(zhì)量代表性樣品測試后進(jìn)行反演對比選取適宜的評價參數(shù)提供了新的思路。同時,也避免了技術(shù)人員在均值、最小值、修正值、經(jīng)驗值之間模糊取值,為參數(shù)確定提供了有利的支撐。
(1)結(jié)合鄭州城區(qū)地基土特征,提出特征土的概念,以奧蘭花園場地為例,對鄭州淺部具經(jīng)濟(jì)價值的特征土進(jìn)行分析,研究鄭州淺部特征土分布情況,采取現(xiàn)場原位測試與代表性試樣采取進(jìn)行室內(nèi)物理力學(xué)性質(zhì)測試相結(jié)合的方法,得出其上部略軟下部偏硬的一般規(guī)律。由于特征土上下部地基承載力相差不大,可根據(jù)擬建多層建筑物荷載大小優(yōu)選可利用部位為主要持力層,以實現(xiàn)其經(jīng)濟(jì)價值。
(2)比較離散元法與有限元法,選用顆粒流離散元法對特征土進(jìn)行了直剪試驗和三軸試驗?zāi)M,與試驗符合一致,考慮直剪試驗存在受限強(qiáng)制性破壞,與實際應(yīng)用存在差距。因此,采用靜三軸模擬進(jìn)行特征土破壞方程構(gòu)建。
(3)在前人研究基礎(chǔ)上,得出特征土靜三軸破壞方程,并通過采樣測試驗證。采用破壞方程反演計算的結(jié)果一致性較好。這對特征土深入研究分析和充分利用提供了便利,也為采樣困難地層在采取有質(zhì)量保證的代表性樣品測試后,進(jìn)行反演比選確定適宜的評價參數(shù)提供了支撐。