用MIDAS模擬樁土相互作用

1、用MIDAS模擬樁-土相互作用（“m法”確定土彈簧剛度）北京邁達斯技術有限公司2009年05月1、引言 6070年代，美國學者J.penzien等在解決泥沼地上大橋動力分析時提出了集中質量法，目前已在國內外得到了廣泛的應用。集中質量法將橋梁上部結構多質點體系和樁一土體系的質量聯(lián)合作為一個整體，來建立整體耦聯(lián)的地震振動微分方程組進行求解。該模型假定樁側土是Winkler連續(xù)介質。以半空間的Mindlin靜力基本解為基礎，將樁-土體系的質量按一定的厚度簡化并集中為一系列質點，離散成一理想化的參數系統(tǒng)。并用彈簧和阻尼器模擬土介質的動力性質，形成一個包括地下部分的多質點體系。土彈簧剛度的確定，除考慮使

2、用較為精確的有限元或邊界元方法外，較為簡便的方法是采用Penzien模型中提供的土彈簧計算方法或參照現行規(guī)范中土彈簧的計算方法。我國公路橋涵地基與基礎設計規(guī)范(JTG D63-2007)用的“m法”計算方法和參數選取方面比Penzien的方法要簡單和方便，且為國內廣大工程師所熟. “m法”的基本原理是將樁作為彈性地基梁，按Winkler假定（梁身任一點的土抗力和該點的位移成正比）求解。但是，由于樁-土相互作用的實驗數據不足，土的物性取值有時亦缺乏合理性，在確定土彈簧的剛度時，仍有不少問題未能很好解決。特別是，“m法”中m的取值對彈簧剛度的計算結果影響很大，且不能反映地震波的頻率特性和強度帶來的

3、影響。本次介紹的土彈簧的模擬是采用規(guī)范中的“m法”確定土的地基系數C（ m的取值根據土的物性而定）,再由其算出土彈簧的水平剛度。 2、采用“m”法，確定土彈簧剛度樁在橫軸向荷載作用下樁身的內力和位移計算，國內外學者提出了許多方法。目前較為普遍的是樁側土采用文克爾假定，通過求解撓曲微分方程，再結合力的平衡條件，求出樁各部位的內力和位移，該方法稱為彈性地基梁法。以文克爾假定為基礎的彈性地基梁法從土力學觀點看是不夠嚴密的，但其基本概念明確，方法簡單，所得結果一般較安全，在國內外工程界得到廣泛應用。我國公路、鐵路在樁基礎的設計中常用的“m”法、就屬此種方法，本文將主要根據“m”法來確定土彈簧

4、剛度。2.1地基系數的概念及確定方法 （1）概念地基系數C表示單位面積土在彈性限度內產生單位變形時所需施加的力，單位為kN/m3或MN/m3。（2）確定方法地基系數大小與地基土的類別、物理力學性質有關。地基系數C值是通過對試樁在不同類別土質及不同深度進行實測及后反算得到。大量的試驗表明，地基系數C值不僅與土的類別及其性質有關，而且也隨著深度而變化。由于實測的客觀條件和分析方法不盡相同等原因，所采用的C值隨深度的分布規(guī)律也各有不同。常采用的地基系數分布規(guī)律有圖下所示的幾種形式，因此也就產生了與之相應的基樁內力和位移的計算方法。圖1 地基系數變化規(guī)律現將樁的幾種有代表性的彈性地基梁計算方法概括在表

5、下中。 樁的幾種典型的彈性地基梁法 表4-1計算方法圖 號地基系數隨深度分布地基系數C表達式說 明m法4-1a）與深度成正比C=mZm為地基土比例系數K法4-1b）樁身第一撓曲零點以上拋物線變化，以下不隨深度變化C=KK為常數C值法4-1c）與深度呈拋物線變化C=cZ0.5c為地基土比例系數張有齡法4-1d）沿深度均勻分布C=K0K0為常數上述的四種方法各自假定的地基系數隨深度分布規(guī)律不同，其計算結果是有差異的。實驗資料分析表明，宜根據土質特性來選擇恰當的計算方法。2.2樁的計算寬度1定義計算樁的內力與位移時不直接采用樁的設計寬度（直徑），而是換算成實際工作條件下相當于矩形截面樁的寬度b1，b

6、1稱為樁的計算寬度。2采用計算寬度的原因為了將空間受力簡化為平面受力，并綜合考慮樁的截面形狀及多排樁樁間的相互遮蔽作用。3計算方法樁的計算寬度b1請按照現行橋涵地基與基礎設計規(guī)范(JTG D63-2007)P83P84頁要求確定。以上的計算方法比較復雜，理論和實踐的根據也是不夠的，因此國內有些規(guī)范建議簡化計算。圓形樁：當d1m時，b1=0.9（1.5d+0.5）；當d1m時，b1=0.9（d+1）。方形樁：當邊寬b1m時，b1=1.5b+0.5；當邊寬1m時，b1=b+1。而國外有些規(guī)范更為簡單：柱樁及樁身尺寸直徑0.8m以下的灌注樁，b1=d+1（m）；其余類型及截面尺寸的樁，b1=1.5d

7、+0.5（m）。2.3剛性樁與彈性樁為計算方便起見，按照樁與土的相對剛度，將樁分為剛性樁和彈性樁。1彈性樁當樁的入土深度時，這時樁的相對剛度小，必須考慮樁的實際剛度，按彈性樁來計算。其中稱為樁的變形系數，2剛性樁當樁的入土深度時，則樁的相對剛度較大，計算時認為屬剛性樁，2.43土的地基比例系數“m”地基土水平抗力系數的比例系數m值宜通過樁的水平靜載試驗確定。但由于試驗費用、時間等原因，某些建筑物不一定進行樁的水平靜載試驗，可采用規(guī)范提供的經驗值如下表所示。 在應用上表時應注意以下事項圖4-5 比例系數m的換算1由于樁的水平荷載與位移關系是非線性的，即m值隨荷載與位移增大而有所減小，因此，m值的

8、確定要與樁的實際荷載相適應。一般結構在地面處最大位移不超過10mm，對位移敏感的結構、橋梁工程為6mm。位移較大時，應適當降低表列m值。2當基樁側面由幾種土層組成時，從地面或局部沖刷線起，應求得主要影響深度hm=2（d+1）米范圍內的平均m值作為整個深度內的m值（見圖4-5），對于剛性樁，hm采用整個深度h。當基礎側面地面或局部沖刷線以下深度內有兩層土時，應將兩層土的比例系數按式（P.0.2-3）換算成一個m值，作為整個深度的m值。當hm（hm>2.5m時）深度內存在兩層不同土時： (4-5)當hm深度內存在三層不同土時： （4-6）3承臺側面地基土水平抗力系數Cn Cn=m·

9、hn （4-7）式中：m承臺埋深范圍內地基土的水平抗力系數，MN/m4； hn承臺埋深，m。4地基土豎向抗力系數C0、Cb和地基土豎向抗力系數的比例系數m0（1）樁底面地基土豎向抗力系數C0 C0=m0h （4-8）式中：m0樁底面地基土豎向抗力系數的比例系數，kN/m4，近似取m0=m； h樁的入土深度(m)，當h小于10m時，按10m計算。（2）承臺底地基土豎向抗力系數Cb Cb=m0hn （4-9）式中：hn承臺埋深(m)，當hn小于1m時，按1m計算。 巖石地基豎向抗力系數C0 表 3-17單軸極限抗壓強度標準值RC（MPa）C0（MN/m3）12530015000注：當RC為表列數值

10、的中間值時，C0采用插入法確定。3、單樁計算1）采用正方形樁寬為1m的，樁長為15m。樁頂假設與地面齊平。根據規(guī)范，各種參數計算如下表。地基土的比例系數m(N/m4)基礎寬度b（m）形狀換算系數k空間工作系數k0=1+1/b基礎的計算寬度b1基礎彈模E(N/m2)2000000011223E+10基礎抗彎I(m4)EI(KNm2)土中基礎變形系數(1/m),=(mb1)/(EI).2基礎入土深度h（m）h<2.5時按照剛性基礎計算，否則彈性基礎0.0833333333E+060.437156.5601722）采用正方形樁寬為1.5m的，樁長為15m。樁頂假設與地面齊平。根據規(guī)范，各種參數

11、計算如下表。地基土的比例系數m(N/m4)基礎寬度b（m）形狀換算系數k空間工作系數k0=1+1/b基礎的計算寬度b1基礎彈模E(N/m2)200000001.511.66672.53E+10基礎抗彎I(m4)EI(KNm2)土中基礎變形系數(1/m),=(mb1)/(EI).2基礎入土深度h（m）h<2.5時按照剛性基礎計算，否則彈性基礎0.4218751E+070.331154.9593443）在上述參數條件下，樁基土彈簧剛度計算見表“M值法計算地基彈性剛度”。根據表中數據在civil中建立樁基模型，進行計算，所得結果與橋博比較，得出下列表格。單樁,樁頂在地面，樁頂受到100KN水平力樁頂位移（mm）樁基最大彎矩（KN/m)midas Civil橋博midas Ci