地震作用下預(yù)應(yīng)力混凝土管樁運動響應(yīng)三維數(shù)值分析

地震作用下預(yù)應(yīng)力混凝土管樁運動響應(yīng)三維數(shù)值分析

陳思奇1，2， 梁發(fā)云1，2， 陳海兵1， 黃茂松1，2

(1.同濟(jì)大學(xué)巖土及地下工程教育部重點實驗室，上海200092； 2.同濟(jì)大學(xué)地下建筑與工程系，上海200092)

摘要：基于黏彈性人工邊界，建立上部結(jié)構(gòu)-樁-土的共同作用三維有限元模型，分析地震作用下預(yù)應(yīng)力混凝土管樁的運動響應(yīng)特性。分別針對預(yù)應(yīng)力混凝土管樁的樁徑、雙層軟硬土剪切波速比值、上覆土層厚度、上部結(jié)構(gòu)荷載等影響因素進(jìn)行數(shù)值計算。參數(shù)分析表明：在地震作用下，樁徑的增大會導(dǎo)致樁身整體彎矩相應(yīng)增加，特別是樁身土層分界面處增大明顯；軟硬土層剪切波速比及上覆土層厚度的增加，引起土層分界面處樁身峰值彎矩增加；固定樁頭條件下，樁頭與樁身軟硬土層分界面處均會產(chǎn)生較大的運動彎矩；上部結(jié)構(gòu)的慣性荷載對固定樁頭的內(nèi)力有著較大影響，對樁身深處段彎矩影響較小。本文研究結(jié)論可為預(yù)應(yīng)力混凝土管樁抗震設(shè)計提供有益的理論參考。

關(guān)鍵詞：預(yù)應(yīng)力管樁； 地震作用； 運動響應(yīng)； 數(shù)值分析； 三維有限元

收稿日期：*2014-08-20

基金項目：國家自然科學(xué)

通訊作者：梁發(fā)云(1976-)，男，安徽肥東人，博士，教授，博士生導(dǎo)師，主要從事樁基礎(chǔ)、深基坑工程及橋梁沖刷等方面的研究與教學(xué)工作.E-mail: fyliang@tongji.edu.cn

中圖分類號:TU47文獻(xiàn)標(biāo)志碼：A

DOI:10.3969/j.issn.1000-0844.2015.01.0055

Characteristics of Bending Moment Distribution in

Prestressed Concrete High Pilings During Earthquakes

CHEN Si-qi1，2， LIANG Fa-yun1，2， CHEN Hai-bing1，2， HUANG Mao-song1，2

(1.KeyLaboratoryofGeotechnicalandUndergroundEngineeringofMinistryofEducation，TongjiUniversity，Shanghai200092，China；

2.DepartmentofGeotechnicalEngineering，TongjiUniversity，Shanghai200092，China)

Abstract：During earthquakes，piles undergo stresses due to both the motion of the superstructure(inertial interaction)and that of the surrounding soil (kinematic interaction).Based on the artificial viscoelastic boundary，a 3D numerical model of seismic performance for superstructure-pile-soil interaction was developed to analyze the kinematic response of PHC pile (prestressed high concrete pipe pile) during earthquakes.The kinematic seismic interaction of single piles embedded in soil was evaluated by focusing on the bending moments induced by the transient motion.Considering the conditions of layered soil，factors influencing the analysis included PHC pile diameter，the ratio of shear wave velocities between two soil layers，the depth of embedded soil，and the inertial loading of the superstructure.The results indicated that pile diameter affects the amplitude of bending moments at the pile head.Specifically，at the interface between layers for a given soil deposit，the bending moment increased as the pile diameter increased，especially at the interface of two soil layers.In soil profiles where the transition between layers was distinct，the bending moments in the pile were significant，especially near the interface of soil layers with highly contrasting stiffness.The bending moment at the interface between two soil layers of soil increased when the ratio of shear wave velocities increased.An increase in the depth of overlying soil in two soil layers would also cause an increase in the pile bending moment at the interface between the two layers.The fixed head of a pile and the interface of two soil layers would increase the bending moment.The superstructure had an important influence on the bending moment of the pile head, but not for the deep section of the pile.The conclusions of this study will lead to recommended reference criteria for the seismic design PHC piles

Key words： prestressed high concrete pipe pile； earthquake action； motion response； bending moment distribution； 3D finite element

0引言

預(yù)應(yīng)力混凝土管樁(簡稱PHC管樁)由于具有較高的豎向承載力和抗彎性能，以及單位豎向承載力造價便宜等優(yōu)勢，在基礎(chǔ)工程中得到了廣泛應(yīng)用。在地震作用下，慣性響應(yīng)和運動響應(yīng)是影響其抗震性能的主要方面。然而，工程師在設(shè)計中較為關(guān)注上部結(jié)構(gòu)產(chǎn)生的慣性荷載對樁基礎(chǔ)的影響，對樁土運動響應(yīng)則大多未加以考慮，在我國的相關(guān)設(shè)計規(guī)范中對結(jié)構(gòu)-樁-土的動力相互作用的分析也不夠充分。PHC管樁在歷次地震中的抗震性能表現(xiàn)不盡如人意，例如日本阪神地震的震害調(diào)查表明，PHC管樁的樁頭與土層界面部位的樁身破壞較為普遍[1]。究其原因，不僅是因為PHC管樁的水平抗剪和抗側(cè)移能力較差，而且地震荷載對樁身的動力響應(yīng)也可能造成樁身應(yīng)力在這些位置處發(fā)生突變。所以地震作用下樁土運動響應(yīng)對PHC管樁樁身受力影響需要進(jìn)行深入分析。

目前，樁土運動響應(yīng)研究主要采用動力Winkler地基模型進(jìn)行分析[2-3]，盡管其簡單實用，但畢竟難以合理地描述半無限連續(xù)土體介質(zhì)的阻尼和上部結(jié)構(gòu)的慣性效應(yīng)。為了合理描述土體的非線性以及樁土界面的分離，需要完整的樁-土-上部結(jié)構(gòu)動力相互作用三維有限元分析。Maiorano等[4]利用三維有限元分析法研究分層土不同剪切波速比下固定樁頭樁身彎矩響應(yīng)。Sica等[5]研究了地震荷載下土體塑性指數(shù)并模擬分析了多種不同典型地震波對樁身運動響應(yīng)的影響，提出相關(guān)經(jīng)驗公式。鄭剛等[6]采用ANASYS軟件建立考慮框架-樁-成層土的二維數(shù)值分析模型，針對樁身內(nèi)力的重要影響因素進(jìn)行參數(shù)分析。

為評估預(yù)應(yīng)力管樁在地震作用下的抗彎性能，針對預(yù)應(yīng)力管樁的樁徑、上覆土層厚度、雙層軟硬土動剪切波速比值、上部結(jié)構(gòu)荷載等影響因素，進(jìn)行綜合全面地分析。本文基于黏彈性人工邊界，采用Abaqus有限元軟件，建立上部結(jié)構(gòu)-樁-土的共同作用的三維有限元模型，分析地震荷載下PHC管樁的運動響應(yīng)特性，為PHC管樁受彎性能的抗震設(shè)計提供依據(jù)。

1有限元模型建立

1.1模型基本尺寸

參照文獻(xiàn)[7]的建模方法，建立如圖1所示模型，模型土體長寬均為40 m，樁身采用外徑為600 mm，壁厚130 mm的B型PHC管樁[8]，樁長均為20 m。上部結(jié)構(gòu)采用等間距矩形質(zhì)量塊代替對模型進(jìn)行簡化，以便對其進(jìn)行定量分析。圖1中H為土層總厚度，H1、H2分別為上覆土與下臥土厚度，VS1、VS2分別為兩層土的剪切波速，L為樁身長度。質(zhì)量塊為100 kN，相隔間距為2 m。

圖1　簡化模型示意圖及土層分布 Fig.1　The simplified model and soil layers distribution

1.2本構(gòu)模型和材料參數(shù)

為簡化計算，樁身、土體、承臺、鋼筋均采用線彈性本構(gòu)模型，并且假定土層間呈均質(zhì)水平分布且各向同性。為減小計算量，取半結(jié)構(gòu)進(jìn)行計算。上部結(jié)構(gòu)-樁-土模型劃分為39 340個單元，土體、承臺、樁身均采用八結(jié)點線性六面體單元、鋼筋采用三結(jié)點二次三維桁架單元模擬，網(wǎng)格劃分如圖2所示，樁身以及樁周土劃分較密，樁底土與上部結(jié)構(gòu)網(wǎng)格劃分較粗。

本文考慮在地震作用下，樁土接觸面可能發(fā)生脫開現(xiàn)象，需要在樁土相互作用中設(shè)置接觸條件，接觸面切向采用罰函數(shù)，接觸面摩擦系數(shù)為0.364，即接觸面內(nèi)摩擦角tan20°；法向采用硬接觸，樁底采用“Tie”與樁底土約束，樁身內(nèi)部鋼筋與樁體相互作用采用內(nèi)置區(qū)域接觸，對稱面采用垂直于立面固定約束。為了方便比較，從樁頭至樁底依次設(shè)置21個等距參考面。

圖2　三維有限元網(wǎng)格 Fig.2　Grid graph of the 3-D finite element model

系統(tǒng)阻尼比ζ采用文獻(xiàn)[3]中的建議值0.05，土體材料阻尼選用Rayleigh阻尼即：

其中α，β值有如下關(guān)系：

通過線性攝動中的頻率提取步進(jìn)行模態(tài)分析，并根據(jù)Amorosi[9]建議方法選擇i、j得到ωi、ωj；帶入上式求得α，β。

B型PHC管樁和承臺的材料參數(shù)如表1所示，土體材料的動力參數(shù)如表2所示。

表 1　樁材料參數(shù)

表 2　土體材料動力參數(shù)

1.3地震波的輸入與邊界條件

地震波輸入采用在模型底面即基巖處施加等效節(jié)點加速度邊界條件的方法實現(xiàn)。本文所選用的地震波為《建筑抗震設(shè)計規(guī)程DGJ08-9-2003》[10]中加速度峰值為35 cm/s2的上海人工波1(以下簡稱SHWN1)，并根據(jù)規(guī)范中規(guī)定，抗震設(shè)防烈度Ⅷ時設(shè)計地震加速度為0.20 g，將SHWN1波乘以常數(shù)加以放大至加速度峰值1.96 m/s2，放大后的加速時程曲線圖與5%加速度反應(yīng)譜如圖3所示。

圖3　SHWN1波加速時程曲線和5%加速度反應(yīng)譜圖 Fig.3　Acceleration time history curve and 5% acceleration response pectrum under SHWN1 wave

在分析土體地震受荷響應(yīng)問題時，真實土體為半無限空間，由于人為截取的地基空間限制，忽略了遠(yuǎn)場剛度和阻尼對系統(tǒng)振動的影響。彈性波會因此在有限域內(nèi)發(fā)生反射形成“盒子效應(yīng)”導(dǎo)致模擬失真。

本文模型邊界條件采用黏彈性人工邊界，即在土體邊界采用相互平行的彈簧和阻尼器組成的單元，防止邊界外發(fā)生地震波反射。單元系數(shù)采用Novak[11]建議取值方法，即：

S1、S2為封閉解得到的無量綱參數(shù)，并給出了剛度阻尼系數(shù)與無量綱頻率的關(guān)系，可得到：

其中Gs為土體剪切剛度；r0為樁軸至邊界距離；ω為荷載主頻率，可將地震波進(jìn)行Fourier變換，由Fourier振幅譜得到荷載主頻率。

2算例驗證

算例1與文獻(xiàn)[4]中有限元模型一致的條件進(jìn)行驗證，樁與上覆土層剛度比為460，阻尼比為5%，上覆土層厚度為15 m，下臥土層厚度15 m。算例2與文獻(xiàn)[5]進(jìn)行驗證，樁與上覆土層剛度比為470，阻尼比為20%，上覆土層厚度為15 m，下臥土層厚度20 m。算例1與算例2均采用地震波Tolmezzo(A-TMZ000)，加速度峰值為0.357 g，上覆土剪切波速為100 m/s，下臥土剪切波速為400 m/s，樁長20 m。加速度時程曲線與5%加速度反應(yīng)譜如圖4所示。

圖4　A-TMZ000波加速時程曲線圖與5%加速度反應(yīng)譜 Fig.4　Accleration time history curve and 5% acceleration response spectrum under A-TMZ000 wave

圖5表示A-TMZ000地震波作用于固定樁頭形式且無上部慣性荷載情況下，其樁身彎矩峰值分別與文獻(xiàn)[4]、[5]中模擬結(jié)果對比。

從圖5可看出：本文所得樁身最大彎矩點出現(xiàn)在土層分界面處，且在此處樁身彎矩變化明顯，樁身上段彎矩變化不大，樁底部彎矩逐漸趨于零，與文獻(xiàn)[4]、[5]中的樁身彎矩變化趨勢基本一致。說明本文模型反應(yīng)樁土相互作用的運動響應(yīng)是可行的。

圖5　樁身彎矩峰值對比圖 Fig.5　Comparison of the peak bending moments at the pile body

3參數(shù)分析

考慮上部結(jié)構(gòu)-樁-土相互作用的PHC管樁在地震荷載下運動響應(yīng)受到多方因素影響。圖6為放大后的SHWN1波作用于圖1模型下得到的樁身峰值彎矩、剪力圖。

圖6　彎矩和剪力樁身峰值圖 Fig.6　Peak bending moments and shearing forces at the pile body

圖6中可看出：在地震荷載作用下，樁頂與土層分界面處的樁身彎矩峰值明顯大于其他部位，各點剪力峰值曲線也有相同的特征，符合文獻(xiàn)[6]得到的地震作用下分層土樁身受力分布。樁頭與土層分界面處的樁身受力狀況是抗震設(shè)計中重點關(guān)注的位置，本文將通過控制變量法，分別考慮樁徑、相鄰?fù)翆蛹羟胁ㄋ俦取⑸细餐翆雍穸取⑸喜拷Y(jié)構(gòu)荷載對地震作用下樁身運動響應(yīng)彎矩的影響規(guī)律。

3.1樁徑對樁身運動響應(yīng)的影響

首先以圖1模型為基礎(chǔ)，考察相同剪切波速比、上覆土層厚度和上部結(jié)構(gòu)荷載，但樁徑不同情況下的樁土運動響應(yīng)，并通過比較0.2、0.4、0.6、0、8 和1.0 m樁徑情況下的樁身彎矩進(jìn)行分析，彎矩峰值包絡(luò)圖如圖7所示。

從圖7中可看出：在剪切波速比、上覆土層厚度和上部結(jié)構(gòu)荷載不變的情況下，隨著樁徑的增大，除樁頭部分變化較小，其余樁段整體彎矩有明顯增大。相較于樁徑為0.2 m的樁身彎矩峰值，0.4、0.6、0.8和1.0 m樁土分界面處最大彎矩峰值增大倍數(shù)分別為3.2、8.9、20.1和41.2倍，當(dāng)樁徑為0.2 m和0.4 m時，樁身最大彎矩出現(xiàn)在樁頭位置；當(dāng)樁徑為0.6、0.8和1.0 m時，樁身最大彎矩出現(xiàn)在土層分界面處。

圖7　不同樁徑對樁身彎矩影響 Fig.7　Effects of different pile’s diameters on bending moments at the pile body

3.2分層土剪切波速比對樁身運動響應(yīng)的影響

本節(jié)考察樁徑、上覆土層厚度和上部結(jié)構(gòu)荷載相同，但土層剪切波速比不同的情況下，通過改變下臥土層的剪切波速從而改變分層土剪切波速比，討論相鄰分層土剪切波速比對樁身運動響應(yīng)的影響，樁身彎矩峰值包絡(luò)圖如圖8所示。

圖8　動剪切波速比對樁身彎矩影響 Fig.8　Effects of different V 2/V 1 on bending moments

從圖8中可看出：相對于雙層土剪切波速比值為1.5情況下，比值為2、2.5和3的峰值彎矩最大值增幅分別為96.4%、178.4%及274.0%。分層土情況下土層分界面處樁身彎矩剪力明顯增大，土層間剪切波速比對樁土分界面處樁身影響也較大。隨著剪切波速比值的增大，樁土分界面處彎矩峰值也隨之增大，而剪切波速比值對樁土分界面處以外的樁身區(qū)段影響較小。

3.3上覆軟土層厚度對樁身運動響應(yīng)的影響

考察樁徑、土層剪切波速比、上部結(jié)構(gòu)相同，但上覆土層厚度不同情況下，即以圖1為基礎(chǔ)，改變H1值進(jìn)行對比分析。圖9表示上覆土層厚度分別為5、10、15和18 m情況下樁身峰值彎矩包絡(luò)圖。

圖9　上覆土層厚度對樁身彎矩影響 Fig.9　Effects of the thickness of overlaying soil on bending moments

從圖9中可看出：樁徑、土層剪切波速比和上部結(jié)構(gòu)荷載相同的情況下，上覆土層厚度H1對樁頭彎矩影響不大，上覆土層厚度對樁土分界面處彎矩峰值有一定影響，樁身峰值最大值出現(xiàn)在上覆土層厚度約為3/4樁長時。樁身彎矩峰值點位置隨上覆土厚度改變而相應(yīng)改變，樁身彎矩變化主要發(fā)生在軟硬土層分界面處。當(dāng)上覆軟土層厚度小于樁長的3/4時，彎矩峰值隨上覆土層厚度增加而增加，當(dāng)上覆土層厚度大于約3/4樁身長度時，彎矩峰值反而隨上覆土厚度增加而略微減小；相較上覆軟土層厚度為5 m的峰值彎矩，上覆土層厚度為10、15和18 m時的峰值彎矩增幅分別為12.5%、35.1%和25.8%。

3.4部結(jié)構(gòu)對于樁身運動響應(yīng)的影響

不考慮上部結(jié)構(gòu)自身受地震荷載內(nèi)力響應(yīng)，僅考慮上部結(jié)構(gòu)所產(chǎn)生的慣性荷載對樁身的影響，并考察相同樁徑、剪切波速比和上覆土層厚度，但上部結(jié)構(gòu)層數(shù)不同，通過改變上部結(jié)構(gòu)層數(shù)，分別計算1～5層上部結(jié)構(gòu)情況下的樁身彎矩峰值包絡(luò)圖。從計算結(jié)果圖10中可知：不同上部結(jié)構(gòu)荷載下，樁頭與樁身土層分界面處均可能成為彎矩最大值點。樁頭彎矩受上部慣性荷載影響明顯，并隨著上部結(jié)構(gòu)質(zhì)量增加而增加，其影響深度也略微隨之增加。相對上部結(jié)構(gòu)數(shù)為1的情況，上部結(jié)構(gòu)數(shù)為2、3、4和5的彎矩峰值增幅分別為81.1%、188.6%、356.1%和634.1%，而上部結(jié)構(gòu)數(shù)對樁身深處彎矩影響較小。

圖10　上部結(jié)構(gòu)荷載對樁身彎矩影響 Fig.10　Effects of the load of superstructure on bending moments

4結(jié)語

本文采用Abaqus軟件建立了三維有限元分析模型，考慮地震作用下，樁-土-上部結(jié)構(gòu)共同作用下預(yù)應(yīng)力混凝土管樁的運動響應(yīng)，并得到以下結(jié)論：

(1) 固定樁頭彎矩較大之外，受分層土影響并不明顯；軟硬土層分界面處樁身彎矩峰值也明顯增大，隨著相鄰?fù)馏w剪切波速差值增大而增大；上覆土層厚度變化對固定樁頭的彎矩影響不大，但會對樁土分界面處彎矩峰值產(chǎn)生一定影響，工程中在土層分界面處進(jìn)行接樁需謹(jǐn)慎處理。

(2) 隨著樁徑的改變，除樁頭部分，其余樁身各處彎矩都受到較大的影響，樁土分界面處彎矩峰值隨著樁徑增加而顯著增大。設(shè)計時應(yīng)考慮樁徑對樁土運動響應(yīng)的影響。

(3) 上部結(jié)構(gòu)慣性荷載作用下，樁頭處會產(chǎn)生較大的彎矩，其彎矩隨上部結(jié)構(gòu)增大而增大。因此，預(yù)應(yīng)力混凝土管樁抗震設(shè)計應(yīng)關(guān)注慣性荷載對受承臺約束的樁頂承載特性的影響。

參考文獻(xiàn)(References)

[1]劉慧珊.樁基震害及原因分析——日本阪神大地震的啟示[J].工程抗震，1999(1)：38-43.

LIU Hui-shan.Damage of Pile Foundations and the Analysis of its Results——the Revelations of Hyogoken-Hanbu Earthquake in Japan[J].Earthquake Resistant Engineering，1999(1)： 38-43.(in Chinese)

[2]Mylonakis G.Simplified Model for Seismic Pile Bending at Soil Layer Interfaces[J].Soils and Foundations，2001，41(4)：47-58.

[3]Dezi F，Carbonari S，Leoni G.Kinematic Bending Moments in Pile Foundations[J].Soil Dynamics and Earthquake Engineering，2010，30(3)：119-132.

[4]Maiorano R M S，Sanctis L D，Aversa S，et al.Kinematic Response Analysis of Piled Foundations Under Seismic Excitation[J].Canadian Geotechnical Journal，2009，46(5)：571-584.

[5]Sica S，Mylonakis G，Simonelli A L.Transient Kinematic Pile Bending in Two-layer Soil[J].Soil Dynamics and Earthquake Engineering，2011，31(7)：891-905.

[6]鄭剛，張楠，巴振寧，等.成層軟土地基PHC管樁抗震性能研究[J].巖土工程學(xué)報，2013，35(增刊)：506-510.

ZHENG Gang，ZHANG Nan，BA Zhen-ning，et al.Seismic Performance of PHC Pipe Piles in Layered Soft Soils[J].Chinese Journal of Geotechnical Engineering，2013，35(S2)：506-510.(in Chinese)

[7]Comodromos E M，Papadopolou M C.Explicit Extension of thep-yMethod to Pile Groups in Cohesive Soils[J].Computers and Geotechnics，2013，47：28-41.

[8]中國建筑標(biāo)準(zhǔn)設(shè)計研究院.10G409預(yù)應(yīng)力混凝土管樁[S].北京：中國計劃出版社，2010.

China Architecture Design and Research Group.10G409 Prestressed Concrete Pipe Pile.Beijing：China Planning Press，2010.(in Chinese)

[9]Amorosi A，Boldini D，Elia G.Parametric Study on Seismic Ground Response by Finite Element Modeling[J].Computers and Geotechnics，2010，37(4)：515-528.

[10]同濟(jì)大學(xué).DGJ09-9-2003建筑抗震設(shè)計規(guī)程[S].上海：上海建設(shè)和管理委員會，2003.

Tongji University.DGJ09-9-2003 Code for Seismic Design of Buildings[S].Shanghai：Shanghai Construction and Management Committee，2003.(in Chinese)

[11]Novak M，Mitwally H.Transmitting Boundary for Axisymmetrical Dilation Problems[J].Journal of Engineering Mechanic，1988，114(1)：No.1，181-187.

[12]Dezi F，Carbonari S，Leoni G. A Model for the 3D Kinematic Interaction Analysis of Pile Groups in Layered Soils[J].Earthquake Engineering and Structural Dynamics，2009，38(11)：1281-1305.

[13]De Sanctis L，Maiorano R M S，Aversa S.A Method for Assessing Kinematic Bending Moments at the Pile Head[J].Earthquake Engineering and Structure Dynamics，2010，39(10)：1133-1154.

[14]Sica S，Mylonakis G，Simonelli A L.Strain Effect on Kinematic Pile Bending in Layered Soil[J].Soil Dynamics and Earthquake Engineering，2013，48：213-242.

[15]Di Laora R，Mylonakis G，Mandolini A.Pile-head Kinematic Bending in Layered Soil[J].Earthquake Engineering and Structural Dynamics，2013，42(3)：319-337.

[16]Anoyatis G，Laora R D，Mylonakis G. Axial Kinematic Response of End-bearing Piles to P Wave[J].International Journal for Number and Analytical Methods in Geomechanics，2013，37(17)：2877-2896.