宋壽鵬，唐益群

(同濟大學 地下建筑與工程系，上海，200092)

20世紀90年代，上海市地下水開采得到控制后，本已有效遏制的地面沉降問題又變得日益嚴重，城市建設引起的工程性地面沉降問題凸顯[1-2]。研究表明：上海城市建設對沉降的影響約占地面沉降總量的30%[3]。2009年與1990年相比，上海高層建筑物的幢數增加了近30倍，建筑面積增加了近24倍[4]。建筑規(guī)模及其增長速度直接導致工程性地面沉降同步增長。建筑容積率是指總建筑面積與建筑用地面積的比值，其間接反映了土地的開發(fā)強度[5]。但為了使經濟效益、社會效益與環(huán)境效益相協(xié)調，建筑容積率需控制在一定范圍內，因此，在詳細規(guī)劃中都對各類居住用地的容積率予以明確規(guī)定。而對于建筑容積率，更多的是基于綠地率，卻忽視了其對城市地質環(huán)境的影響。地面沉降效應是建筑容積率對城市地質環(huán)境影響的直接表現，建筑密度越大，建筑容積率越高，地面沉降量越顯著，地面沉降的影響范圍越大[6]，產生的城市環(huán)境問題可能會越多。為此，本文作者采用改進的西原模型描述土的非線性特性，運用FLAC3D軟件[7]建立密集高層建筑群地面沉降模型，模擬筏板-樁-土間的相互作用，以考察不同建筑容積率條件下筏形基礎中樁基受力變化及區(qū)域地面沉降規(guī)律等，并基于沉降控制條件探討最佳建筑容積率的取值范圍。

1 三維數值模型的建立

密集高層建筑群地面沉降模型是以上海軟土工程地質條件作為背景，地層主要由飽和的黏性土、粉性土等組成，具有成層分布的特點，其主要物理力學指標如表 1所示[8]。平均地下水位處于距地面 0.5 m以下。

1.1 地層本構模型

由文獻[9]可知：上海軟土的應力-應變-時間關系符合西原模型，見圖1。圖1中：Ge為Hooker體的剪切模量；Ke為Hooker體的體積模量；Gve為Kelvin體彈簧的剪切模量；Kve為Kelvin體彈簧的體積模量；ηve為Kelvin體黏壺的黏滯系數；ηvp為Bingham體黏壺的黏滯系數；εe為 Hooker體的應變；εve為 Kelvin體的應變；εvp為Bingham體的應變。西原模型可以考慮材料的彈性、黏彈及彈塑等力學特征，適用范圍較廣，理論較成熟[10]。本文采用簡單的Mohr-Coulomb屈服函數替換線性牛頓體。根據文獻[11]可得到三維狀態(tài)下改進的西原模型新偏應力的表達式為：

圖1 西原流變模型三維情形Fig.1 Nishihara viscous-elastic-plastic model in three dimension conditions

可根據上述改進的西原模型的應力增量表達式，利用面向對象的語言標準C++編寫相應的有限差分本構程序代碼，通過UDM接口程序嵌入FLAC3D中。

根據文獻[9]，G和Gve均與加載的應力p有關：

改進的西原模型中共有8個參數，⑦-1和⑦-2層采用退化的Burger模型，其余各層均采用改進的西原模型；除⑦-1層和⑦-2層的p=σ1外，其余各層的。模型計算參數如表2所示[9]。

表1 數值模型各地層的物理力學指標Table 1 Physico-mechanical properties of soils for numerical model

表2 地層模型計算參數[9]Table 2 Calculation parameters of strata model

1.2 三維模型數值模擬

1.2.1 模型的邊界條件

模型側面和底面為位移邊界，側面限制水平移動；底部為固定邊界，限制水平移動和垂直移動。模型上面是地表，取為自由邊界。

1.2.2 模型的初始應力場

垂直地壓為上覆蓋層的重度，水平地壓為垂直地壓乘以側壓力系數0.65[13]，靜水壓力采用水土耦合計算；利用Mohr-Columd彈塑性模型生成初始應力場。

1.2.3 樁筏基礎模擬

在FLAC3D程序中分別采用pile和liner單元模擬樁和筏板。具體參數取值見表3～4[14]。

1.2.4 模擬條件

密集高層建筑群地面沉降模型見圖 2，模型中的最大建筑基地面積為180 m×180 m；4棟建筑物的布置方式采用上海住宅區(qū)常用的平行式，其樁筏基礎尺寸完全相同(樁長為56 m，樁徑為0.4 m，樁間距為2 m，筏板長×寬為12 m×12 m)，建筑間距為20 m。樁基礎的尺寸滿足單樁豎向承載力要求；由文獻[15]可知：筏板的尺寸只對建筑物的豎向沉降量有影響，而與地面沉降影響范圍無關，故可依據實際計算要求選擇筏板尺寸。

在沉降模型中，建筑物底層面積和模型占地面積一定的條件下，建筑容積率僅與建筑面積有關。建筑面積對建筑基礎的影響僅體現于基礎附加應力方面，以建筑模型增重的方式來實現[16]。模型中建筑物層數分別設定為18層、24層、32層，其對應的最小建筑容積率分別為0.08，0.11和0.14。的位置有很大關系。

表3 pile單元參數Table 3 Parameters of pile element

圖2 模型剖面位置及典型樁位置圖Fig.2 Location of model profiles and typical piles

表4 liner單元參數[14]Table 4 Parameters of liner element

2 數值模擬結果分析

2.1 樁基受力分析

利用上述參數對不同建筑容積率下密集高層建筑群的工作性狀進行模擬，得到圖2中1號建筑物的6根典型樁的樁身軸力分布圖及樁側阻力圖，分別見圖3和圖4。

從圖3和圖4可見：(1) 根樁的樁身軸力及樁側阻力均隨建筑層數的增加而增大，樁身軸力及樁側阻力的增加最終導致樁筏基礎沉降量的增大；(2) 在同一建筑層數下，樁側阻力由大至小分別為 3號樁、2號樁、4號樁、1號樁、5號樁和6號樁。但考慮2號樁與6號樁、3號樁與5號樁在1號建筑物中位置的對稱關系，建筑物間存在明顯的應力疊加效應才造成了2組樁樁側阻力的差異，這種差異性與各樁在模型

2.2 區(qū)域地面沉降規(guī)律

為了更好地描述不同建筑層數下密集建筑群地面沉降模型中的區(qū)域地面沉降規(guī)律(見圖 5)，選取 2個點、3個剖面隨不同建筑層數的變化特征來討論。所選擇的2個點為模型中心A點、兩建筑物中心B點，3個剖面為過A和B這2個點的Ⅰ-Ⅰ剖面、過兩建筑物及B點的Ⅱ-Ⅱ剖面、過兩建筑物及A點的Ⅲ-Ⅲ剖面。2個點與3個剖面的具體位置如圖2所示。由模型的模擬結果即得出A點和B點沉降量隨時間的變化曲線，見圖6。

從圖6可以看出：(1) A和B 2點的沉降量均隨時間的增加而增大，第9年末A點的沉降均穩(wěn)定，第11年末B點的沉降全部穩(wěn)定；(2) A和B 2點的最大沉降量均隨建筑層數的增加而增大；(3) A和B 2點的3條沉降曲線趨勢基本一致，緩慢沉降、加速沉降、沉降穩(wěn)定的時間段基本相同。

圖3 模型典型樁的樁身軸力分布圖Fig.3 Press about typical piles of model

圖4 模型典型樁的樁側阻力圖Fig.4 Skin friction about typical piles of model

圖5 區(qū)域地面沉降變化圖Fig.5 Variations of regional land subsidence

圖6 模型典型點沉降量隨時間的變化曲線Fig.6 Curves of settlement with time about typical points of model

11年末不同建筑層數下各剖面沉降量曲線見圖7，11年末密集建筑群與單棟建筑下Ⅱ-Ⅱ剖面沉降曲線見圖8。從圖7可見：(1) 建筑層數的改變對4棟建筑物的中間區(qū)域影響較大，而對外部區(qū)域的影響較小，在Ⅰ-Ⅰ剖面上沉降量最大的點出現在A和B 2點之間(A點和B點位置見圖2)；(2) Ⅱ-Ⅱ剖面的沉降量以建筑物筏板上某點最大，Ⅱ-Ⅱ剖面上各點的沉降量都隨建筑間距的增加而減小；(3) Ⅲ-Ⅲ剖面的沉降量也以建筑物筏板上某點最大；Ⅲ-Ⅲ剖面上各點的沉降量也都隨建筑間距的增加而減小。

從圖8可見：密集建筑群條件下與單棟建筑條件下Ⅱ-Ⅱ剖面的沉降曲線相比，影響范圍重合于距筏板中心點48，50和52 m處。此現象說明密集建筑群區(qū)的應力疊加效應明顯，應力疊加效應的范圍也隨附加應力即建筑層數的增加而增大。

圖7 11年末不同建筑層數下各剖面沉降量曲線Fig.7 Settlement curves of profiles under different layers in end of 11 years

2.3 最佳建筑容積率的確定

建筑層數的改變不論是對剖面上各點的沉降量還是對密集建筑群影響范圍都有明顯影響。從城市地質環(huán)境保護的角度看，需盡可能減小密集建筑群間的應力疊加效應，以有效控制沉降量。這就要使密集建筑群的建筑用地面積盡可能滿足其造成地面沉降的最大影響范圍，而建筑容積率對其影響范圍起決定作用，因此，最佳建筑容積率的確定對城市地質環(huán)境保護具有重要意義。

圖8 11年末密集建筑群與單棟建筑下Ⅱ—Ⅱ剖面沉降曲線對比Fig.8 Settlement curves of Ⅱ—Ⅱ profile between intensive buildings and single building in the end of 11 years

11年末不同建筑層數下Ⅱ-Ⅱ剖面相鄰2點沉降差與距離的關系曲線見圖 9。參照《工程測量規(guī)范》(GB 50026—2007)，根據相鄰2點沉降量的差值可近似得到密集建筑群區(qū)的最大影響范圍，從而得到不同建筑層數下Ⅱ-Ⅱ剖面相鄰 2點沉降差與距離的關系曲線，并選取0.5 mm作為沉降差界限，并以此從圖9中估算不同建筑層數下密集建筑群區(qū)的最大影響范圍。

圖9 11年末不同建筑層數下Ⅱ-Ⅱ剖面相鄰2點沉降差與距離的關系曲線Fig.9 Curves of adjacent points between differential settlement and distance under Ⅱ-Ⅱ profile of different layers in the end of 11 years

依據圖9得到的影響范圍也就是建筑用地面積，利用公式計算不同建筑層數對應的建筑容積率，得到散點如圖10所示，經擬合得到散點趨勢線便可取得密集建筑群的最佳容積率為0.7～1.2。

圖10 11年末建筑層數與建筑容積率的關系曲線Fig.10 Curves between layers and FAR in the end of 11 years

3 結論

(1) 樁筏基礎中樁身軸力及樁側阻力均隨建筑層數的增加而增大，樁身軸力及樁側阻力的增加最終導致樁筏基礎沉降量增大；與單棟建筑物情況相比，各樁的力學性質均受其他建筑物的影響，各樁位置的差異造成了各樁身應力疊加效應和樁側摩阻力的差異。

(2) 建筑層數的改變對密集建筑群的中間區(qū)域影響較大，最大沉降處仍為建筑物自身；密集建筑群的最大影響范圍隨建筑層數的增加而增大。

(3) 在控制地面沉降的條件下，以0.5 mm作為沉降差界限估算密集建筑群的最大影響范圍，并據此得到最佳建筑容積率的取值范圍為0.7～1.2。

[1] TANG Yi-qun, CUI Zhen-dong. Model test study of land subsidence caused by the high-rise building group[J]. Bull Eng Geol Environ, 2008, 67(2): 173-179.

[2] TANG Yi-qun, CUI Zhen-dong. Application of grey theorybased model to prediction of land subsidence due to engineering environment in Shanghai[J]. Environ Geol, 2008, 55(3):583-593.

[3] 介玉新, 高燕, 李廣信. 建筑物間距的統(tǒng)計規(guī)律及對附加應力的影響[J]. 工業(yè)建筑, 2010, 40(S): 62-65.JIE Yu-xin, GAO Yan, LI Guang-xin. Statistical pattern of building spacing and its influence on additional stress[J].Industrial Construction, 2010, 40(S): 62-65.

[4] 上海市統(tǒng)計局. 上海統(tǒng)計年鑒[M]. 北京: 中國統(tǒng)計出版社,2010: 172-173.Shanghai Statistic Bureau. Statistical yearbook of Shanghai[M].Beijing: China Statistic Press, 2010: 172-173.

[5] 趙延軍, 王曉鳴. 開發(fā)項目最佳容積率研究[J]. 長安大學學報: 社會科學版, 2008, 10(2): 92-95.ZHAO Yan-jun, WANG Xiao-ming. Best floor-area ratio of development project[J]. Journal of Changan University: Social Science Edition, 2008, 10(2): 92-95.

[6] 嚴學新, 龔士良, 曾正強, 等. 上海城區(qū)建筑密度與地面沉降關系分析[J]. 水文地質工程地質, 2002(6): 21-25.YAN Xue-xin, GONG Shi-liang, ZENG Zheng-qiang, et al.Relationship between building density and land subsidence in Shanghai urban zone[J]. Hydrogeology and Engineering Geology, 2002(6): 21-25.

[7] Itasca Consulting Group Inc. FLAC3D(Version 2.1) users manual[R]. Minnesota: Itasca Consulting Group Inc, 2002: 1-4.

[8] 龔劍, 趙錫宏. 對 101層上海環(huán)球金融中心樁筏基礎性狀的預測[J]. 巖土力學, 2007, 28(8): 1695-1699.GONG Jian, ZHAO Xi-hong. Prediction of behavior of piled raft foundation for Shanghai World Financial Center of 101-storey[J].Rock and Soil Mechanics, 2007, 28(8): 1695-1699.

[9] 詹美禮, 錢家歡, 陳緒祿. 軟土流變特性試驗及流變模型[J].巖土工程學報, 1993, 15(3): 54-62.ZHAN Mei-li, QIAN Jia-huan, CHEN Xun-lu. Tests on rheological behavior of soft soil and rheologic model[J]. Chinese Journal of Geotechnical Engineering, 1993, 15(3): 54-62.

[10] 孫鈞. 巖土材料流變及其工程應用[M]. 北京: 中國建筑工業(yè)出版社, 1999: 459-467.SUN Jun. Rheology of geo-materials and its engineering application[M]. Beijing: China Architecture and Building Press,1999: 459-467.

[11] 褚衛(wèi)江, 徐衛(wèi)亞, 楊圣奇, 等. 基于 FLAC3D巖石黏彈塑性流變模型的二次開發(fā)研究[J]. 巖土力學, 2006, 27(11):2005-2010.CHU Wei-jiang, XU Wei-ya, YANG Sheng-qi, et al. Secondary development of a viscoelasto-plastic rheological constitutive model of rock based on FLAC3D[J]. Rock and Soil Mechanics,2006, 27(11): 2005-2010.

[12] 張云, 薛禹群, 吳吉春, 等. 上海砂土蠕變變形特征的試驗研究[J]. 巖土力學, 2009, 30(5): 1226-1236.ZHANG Yun, XUE Yu-qun, WU Ji-chun, et al. Experimental research on creep of Shanghai sands[J]. Rock and Soil Mechanics, 2009, 30(5): 1226-1236.

[13] 姜洪偉, 趙錫宏, 張保良. 各向異性條件下軟土深基坑抗隆起穩(wěn)定性分析[J]. 巖土工程學報, 1997, 19(1): 1-7.JIANG Hong-wei, ZHAO Xi-hong, ZHANG Bao-liang. Analysis of heave-resistant stability for deep braced excavation in soft clay under anisotropic condition[J]. Chinese Journal of Geotechnical Engineering, 1997, 19(1): 1-7.

[14] 杜曉霞. 上部結構-樁筏基礎-地基-共同作用的非線性數值模擬[D]. 南京: 南京理工大學理學院, 2005: 40-41.DU Xiao-xia. Experimental study on small scale tri-functional water-source (ground-source) heat pump central air-conditioning unit[D]. Nanjing: Nanjing University of Science and Technology.School of Sciences, 2005: 40-41.

[15] 樓曉明, 劉建航, 洪毓康. 兩個相鄰高層建筑樁箱基礎的共同作用分析[J]. 土木工程學報, 2002, 35(3): 55-60.LOU Xiao-ming, LIU Jian-hang, HONG Yu-kang. Study on interaction of two adjacent pile-box foundations[J]. China Civil Engineering Journal, 2002, 35(3): 55-60.

[16] 唐益群, 宋壽鵬, 陳斌, 等. 不同建筑容積率下密集建筑群區(qū)地面沉降規(guī)律研究[J]. 巖石力學與工程學報, 2010, 29(S1):3425-3431.TANG Yi-qun, SONG Shou-peng, CHEN Bin, et al. The study of dense buildings group on ground subsidence rule in soft soil area[J]. Chinese Journal of Rock Mechanics and Engineering,2010, 29(S1): 3425-3431.