李永波， 張鴻儒

(1.天津市市政工程設(shè)計研究院，天津 300050； 2.北京交通大學(xué)土木建筑工程學(xué)院，北京 100044)

基于動力BNWF法的凍土-樁相互作用模型研究①

李永波1， 張鴻儒2

(1.天津市市政工程設(shè)計研究院，天津 300050； 2.北京交通大學(xué)土木建筑工程學(xué)院，北京 100044)

基于水平循環(huán)荷載作用下不同負(fù)溫凍土環(huán)境中單樁動力特性模型試驗結(jié)果，在已有分析樁-土-結(jié)構(gòu)相互作用的動力BNWF模型的基礎(chǔ)上，提出改進(jìn)的凍土-樁基動力相互作用非線性反應(yīng)分析模型。在該模型中，利用改進(jìn)的雙向無拉力多段屈服彈簧考慮樁側(cè)凍土的水平非線性力學(xué)特性，同時兼顧樁側(cè)與凍土間的豎向非線性摩擦效應(yīng)、樁尖土的擠壓與分離作用以及遠(yuǎn)場土體阻尼對樁基動力特性的影響。其中樁側(cè)水平多段屈服彈簧參數(shù)根據(jù)凍土非線性p-y關(guān)系獲得，該關(guān)系曲線以三次函數(shù)曲線段及常值函數(shù)段共同模擬，并由室內(nèi)凍土壓縮試驗結(jié)果確定。最后基于改進(jìn)的動力BNWF模型，提取動位移荷載作用下該樁頂力-位移滯回曲線及樁身不同埋深處的彎矩動響應(yīng)數(shù)值分析結(jié)果，并與相應(yīng)的模型試驗結(jié)果對比，二者具有較好的擬合度，表明本文所提出的改進(jìn)模型在分析凍土-樁動力相互作用時有較好的適用性。

凍土-樁相互作用； 模型試驗； 地基非線性； 動力BNWF模型； 雙向無拉力多段屈服彈簧

0 引言

多年凍土在我國分布廣泛，約占21.5%的國土總面積，且主要位于強(qiáng)震高發(fā)區(qū)。對于凍土區(qū)建筑和穿越不穩(wěn)定凍土帶的橋梁結(jié)構(gòu)等重要工程來說，樁基礎(chǔ)以其良好的場地適應(yīng)性及強(qiáng)大的抵抗復(fù)雜荷載能力成為該類工程的主要基礎(chǔ)形式。凍土力學(xué)特性受溫度、含水量等因素影響強(qiáng)烈，強(qiáng)動載作用下凍土中樁基的損傷情況遠(yuǎn)比一般土體復(fù)雜。為研究不同凍土環(huán)境中樁基的動力響應(yīng)特性，國內(nèi)外學(xué)者開展了一些探索試驗：Suleiman等[1-2]研究了季節(jié)凍土中大比例模型柱-樁結(jié)構(gòu)在循環(huán)擬靜力作用下反應(yīng)情況，Wotherspoona等[3]對該試驗進(jìn)行了數(shù)值模擬與分析；Vaziri和Han[4-5]研究了季節(jié)凍土中足尺混凝土樁基在不同激振情況下動力響應(yīng)情況，并將試驗結(jié)果驗證其邊界域理論的正確性；吳志堅等[6]主要對高溫多年凍土區(qū)樁基礎(chǔ)在地震荷載作用下樁側(cè)凍土溫度變化進(jìn)行了研究；李永波等[7]通過模型試驗研究了凍土凍融作用對樁基的水平動力響應(yīng)的影響情況。以上研究均表明樁基礎(chǔ)的動力特性受凍土負(fù)溫及凍融作用影響顯著，在強(qiáng)水平動載作用下樁基礎(chǔ)與周圍土體容易進(jìn)入強(qiáng)非線性狀態(tài)，樁基與土體之間會產(chǎn)生擠密與分離效應(yīng)，因此在對凍土區(qū)樁基礎(chǔ)設(shè)計時要重視土-結(jié)相互作用的影響，并考慮地基土的非線性作用。

近年來，圍繞樁-土-結(jié)構(gòu)相互作用問題的研究不斷深入。劉曉強(qiáng)等[8]基于室內(nèi)動三軸試驗對振動前后樁周土體的強(qiáng)度弱化情況進(jìn)行了研究，并提出一個較簡單且符合試驗規(guī)律的天津地區(qū)樁周土的強(qiáng)度弱化公式；El Naggar等[9]分析海岸工程樁基在非線性地震動輸入下的動力響應(yīng)時，提出了一種以簡單非線性彈簧、黏壺及接觸單元考慮了樁-土界面處分離效應(yīng)的簡化動力BNWF模型；Allotey和El Naggar[10]提出了一種分析淺埋及深埋基礎(chǔ)地震響應(yīng)的動力BNWF模型，該模型考慮了樁周土體的加-卸載準(zhǔn)則、變形區(qū)累積發(fā)展過程以及強(qiáng)度的周期衰減和輻射阻尼；Allotey等[11]考慮樁-土的p-y和t-z耦合效應(yīng)，研究了無黏性土中單樁在豎向和水平向地震動作用下動力響應(yīng)。

本文針對凍土環(huán)境中樁基在水平動載作用下的非線性響應(yīng)特性進(jìn)行研究，結(jié)合模型試驗結(jié)果，給出凍土-樁之間連接彈簧的非線性p-y關(guān)系的確定方法，并基于動力BNWF模型理論，提出雙向無拉力多段屈服彈簧模型和改進(jìn)的樁基礎(chǔ)結(jié)構(gòu)非線性地震反應(yīng)分析模型，并進(jìn)行模型試驗驗證。利用改進(jìn)的分析模型有助于對凍土環(huán)境中樁基結(jié)構(gòu)的地震反應(yīng)及變形規(guī)律進(jìn)行研究。

1 凍土-樁動力相互作用模型試驗

為分析多年凍土環(huán)境中樁基在水平動載作用下的動力響應(yīng)特性，對該類場地樁基橋梁結(jié)構(gòu)在地震作用下的損害情況進(jìn)行預(yù)測，開展凍融條件下模型樁基水平動力試驗研究[7]。

1.1 試驗簡介

基于北京交通大學(xué)巖土工程中心自主研制的凍土樁動力相互作用模型試驗系統(tǒng)，以外徑32 mm、壁厚3 mm、長650 mm(其中埋置于土體中長度約為600 mm)的Q235無縫鋼管作用動力試驗?zāi)Ｐ蜆叮瑢鶞胤謩e為-5 ℃、-3 ℃及淺表層解凍凍土環(huán)境中模型樁的水平動力響應(yīng)特性進(jìn)行試驗研究。試驗時動載頻率為2 Hz，每級荷載持續(xù)約30 s。

1.2 試驗結(jié)果

以-3 ℃均溫凍土中模型樁響應(yīng)結(jié)果為例，在幅值設(shè)定值1.5、2.5、3.5 kN 3級正弦動載作用下實測樁頭的位移-荷載關(guān)系及樁身彎矩(圖1)。

主要試驗結(jié)果如下：

(1) 樁基前后凍土因受壓流變產(chǎn)生塑性變形，樁與凍土接觸面間出現(xiàn)較大間隙；土溫越低，相同動載作用下的間隙越小。

圖1 -3 ℃凍土中模型樁動力響應(yīng)圖Fig.1 Dynamic responses of model pile embeded in frozen soils at -3 ℃

(2) 由于樁-土間的分離效應(yīng)，樁基在動載牽引下經(jīng)過位移零值時荷載曲線上會出現(xiàn)平緩變化段，表現(xiàn)在位移-荷載曲線呈現(xiàn)反S型。

(3) 樁基水平剛度隨樁頭承受動載的幅值增大逐漸減小。

(4) 同一凍土中模型樁基在不同動載幅值作用下樁身彎矩分布規(guī)律基本相同。隨著土體溫度升高，樁身最大彎矩點埋深逐漸增大，且最大彎矩呈現(xiàn)出增大趨勢。

2 樁-土-結(jié)構(gòu)相互作用分析模型

2.1 分析模型

為描述樁基結(jié)構(gòu)在水平向地震動作用下的響應(yīng)特性，針對樁周土體的力學(xué)狀態(tài)及其對樁基動力阻抗行為的影響，學(xué)者提出了多種力學(xué)模型。較常見的有Penzien模型及改進(jìn)Penzien模型[12](圖2)、橋梁樁基平面桿系結(jié)構(gòu)模型，以及大量簡化的動力BNWF模型，如Kagawa模型[13]、El Naggar和Novak模型[14]與Boulanger模型[15](圖3)等。

圖2 改進(jìn)的Penzien模型Fig.2 Improved Penzien model

圖3 Boulanger分析模型Fig.3 Boulanger analytical model

以上各分析模型均有其優(yōu)缺點及適用范圍，其中改進(jìn)的Penzien模型能夠考慮軸力變化的影響，同時直接考慮多排樁，避免了設(shè)置基礎(chǔ)轉(zhuǎn)動彈簧，并把相應(yīng)估計范圍內(nèi)等效自由場地的質(zhì)量以質(zhì)點的形式模擬。Boulager模型考慮的因素較為全面：以拖拉和分離彈簧考慮樁-土間的分離閉合效應(yīng)，以塑性彈簧考慮土體的屈服效應(yīng)，以彈簧-阻尼并聯(lián)單元模擬遠(yuǎn)場土體。

基于Boulanger模型的思想，結(jié)合前述試驗結(jié)果，本文提出一種可以考慮動載作用下樁基與其周圍凍土分離及摩擦效應(yīng)的樁-土相互作用計算模型，如圖4所示。

圖4 本文樁-土動力相互作用計算模型Fig.4 The calculation model for frozon soil-pile dynamic interaction proposed in this study

該模型具有以下特點：

(1) 考慮樁-土水平向接觸作用以及樁、土材料非線性特性；凍土水平彈簧為雙向無拉力多段屈服彈簧單元，可以考慮樁-土間的分離與擠壓；

(2) 考慮樁基與周圍土體的豎向摩擦作用，按照雙線性屈服準(zhǔn)則的豎向摩擦非線性彈簧單元進(jìn)行模擬；

(3) 考慮樁底與土體的擠壓作用，并以無拉力彈簧模擬；

(4) 考慮遠(yuǎn)場土體阻尼對樁基動力特性的影響。

3 樁-土彈簧參數(shù)的確定

3.1 凍土非線性p-y曲線

樁基受水平動載作用時，土體的非線性p-y關(guān)系是樁-土相互作用研究中最受關(guān)注的部分。土體凍結(jié)時的剛度約為其非凍結(jié)時的兩個數(shù)量級倍，當(dāng)有足夠量冰晶充填于土顆粒間空隙時，凍土的力學(xué)特性與弱巖相似。因此，本文基于模型試驗數(shù)據(jù)，采用修正的弱巖p-y曲線建立文中凍土的非線性p-y關(guān)系。

圖5 凍土p-y曲線示意圖Fig.5 Sketch of the p-y curve for frozen soil used in this study

如圖5所示，本文所采用的凍土非線性p-y關(guān)系曲線由兩部分線段組成：第一部分為三次函數(shù)曲線段，另一部分則為常值函數(shù)段。兩條線段的函數(shù)方程表達(dá)式為

本文主要對根據(jù)雙目視覺的原理，通過雙目攝像機(jī)獲取場景的信息，得到距離信息，進(jìn)而實現(xiàn)避障的方法進(jìn)行研究，利用opencv實現(xiàn)對攝像機(jī)標(biāo)定、特征點檢測與匹配、距離信息提取等方面進(jìn)行探究，設(shè)計了一套可行的方案。但在匹配精度和測距精度上仍有一些地方需要改進(jìn)。通過實現(xiàn)雙目測距，可以實現(xiàn)測速，進(jìn)而可以實現(xiàn)無人機(jī)導(dǎo)航勘測等功能。

(1)

式中：pu為凍土的極限抗阻力；ym為凍土體所受荷載壓力為pu值一半時的樁身變形，二者可由式(2)及式(3)求解獲得。

(2)

(3)

式中：qu為凍土的壓縮強(qiáng)度，可根據(jù)文獻(xiàn)[16]中凍土壓縮試驗結(jié)果獲得，如式(4)所示；d為樁基直徑或?qū)挾?；xfs為地表以下凍土深度；km為凍土極限強(qiáng)度50%時相應(yīng)的應(yīng)變，本文中取km=1/25。

(4)

式中：T為凍土負(fù)溫。

根據(jù)凍土材料的非線性p-y關(guān)系，在本文所提出的凍土-樁動力相互作用計算模型中，凍土水平彈簧設(shè)為雙向無拉力三段屈服彈簧單元。該彈簧模型卸載斜率與初始加載斜率相同；循環(huán)動力加載過程中，彈簧受力與變形過程如圖6所示。

圖6 雙向無拉力多段屈服彈簧模型Fig.6 Model of the bidirectional tensionless multi-section yield spring

雙向無拉力多段屈服非線性彈簧可以考慮土體分層的影響，在樁基的不同深度處根據(jù)土介情況施加如上所述的不同剛度、屈服力的無拉力多段屈服彈簧即可。

3.2 遠(yuǎn)場彈簧剛度與阻尼

連接遠(yuǎn)場土體與樁周土體的水平簧剛度KHi可以根據(jù)Mindlin介質(zhì)中位移應(yīng)力分層計算歸結(jié)到各層界面上得到[17]，計算公式為：

(5)

式中：Ei為第i層土的彈性模量；hi為第i層土的厚度；zi為第i層土的中心深度；B為樁半徑。

遠(yuǎn)場土體水平阻尼系數(shù)CH的確定可根據(jù)Lysmer等[17]提出的用黏性阻尼器模擬波動能量逸散的方法，計算公式為：

(6)

式中：B為樁半徑；hi為第i層土的厚度；cP、cS分別為場地土體縱波波速及剪切波波速，可由式(7)及式(8)求解；ρ為土體密度；i=2，3，…。

(7)

(8)

式中：λ為拉梅常數(shù)；G為土體剪切模量。

3.3 樁周及樁底非線性彈簧

樁基受載過程中，若樁尖的豎直地基抗力達(dá)到樁尖支承力、樁周的豎向地基抗力達(dá)到計算范圍內(nèi)的最大樁周摩阻力，兩者接觸土體均進(jìn)入塑性。其中樁尖豎直彈簧只承壓而無拉力特性。地基土彈簧的抗力與位移關(guān)系的骨架曲線如圖7所示，其值求解公式及方法參見文獻(xiàn)[18]。

圖7 樁周及樁底彈簧骨架曲線Fig.7 Skelton curve of soil springs in pile bottom and around the pile

4 數(shù)值算例

根據(jù)本文提出的凍土-樁動力相互作用計算模型，結(jié)合文獻(xiàn)[16]凍土動力特性試驗獲得的-3 ℃時凍土的彈性模量等，可確定樁側(cè)水平地基土雙向無接力多段屈服非線性彈簧、樁周土豎直非線性彈簧、樁底彈簧以及遠(yuǎn)場彈簧及阻尼等參數(shù)?？紤]到模型試驗中樁頭受載主要由近似按正弦循環(huán)動位移控制，因此在數(shù)值模型分析時，亦在模型樁頭施加一按正弦規(guī)律變化的動位移，動載頻率同模型試驗，為2 Hz。

對照模型試驗鋼管樁各參數(shù)，數(shù)值分析中樁長650 mm、樁徑32 mm、壁厚3 mm，其中埋置于土體中的長度為600 mm。鋼管樁按隨動雙線性屈服材料模擬，屈服前彈性模型為206 GPa，屈服應(yīng)力為235 MPa；屈服后切線模量為20.6 GPa。

4.1 樁頭位移-荷載關(guān)系

提取幅值分別為2 mm、3 mm及4 mm的動位移作用下模型樁樁頭位移-荷載關(guān)系曲線，如圖8所示。

從圖8中可以看出，在不同等級的荷載作用下，數(shù)值分析所得到樁頭的位移-荷載滯回曲線與模型試驗所得形狀極為近似，且數(shù)值大小相近，較好地模擬了試驗結(jié)果，證明本章所提出的基于動力BNWF法的樁-土動力相互作用分析模型是可行的。

圖8 不同動載下樁頭滯回曲線Fig.8 Displacement-force responses of the pile head under different dynamic loadings

數(shù)值模擬與模型試驗所得樁頭滯回曲線形狀均呈反S形，即曲線在位移零點左右明顯收縮，而在相對較大位移處漸變?yōu)樨S滿，間接說明了在位移零點處樁土分離后，由鋼管樁承擔(dān)較多動荷載。隨著動荷載增大，曲線在位移零點的飽滿度增大，表明此時鋼管樁已發(fā)生一定程度的屈服。

圖8中不同峰值動位移作用下樁頭滯回曲線與試驗所得曲線的差異表現(xiàn)在：在荷載作用下樁頭的位移及荷載大小基本上呈對稱分布，而試驗所得為模型樁受推一側(cè)的荷載明顯較受拉側(cè)大，這是由試驗條件局限導(dǎo)致的。

4.2 樁身動力響應(yīng)分析

水平動載作用下樁身各點的位移、彎矩響應(yīng)情況是樁-土動力相互作用研究中比較關(guān)心的問題。提取峰值為3 mm的正弦動載作用2 s時間內(nèi)數(shù)值，分析模型中樁身距土表不同深度處的位移和彎矩時程，如圖9所示。從圖中可以看出：

圖9 動載作用下樁身位移和彎矩響應(yīng)情況圖Fig.9 Responses of displacement and bending moment of pile shaft under dynamic loadings

(1) 樁基在正弦動載作用下，其不同埋深點的位移、彎矩時程響應(yīng)也近似按正弦規(guī)律變化，但相對于樁頭，其響應(yīng)值有不同程度的滯后性。

(2) 動載作用下樁身的振動響應(yīng)是一個復(fù)雜的過程，總體上隨著樁身埋置深度的加大，位移響應(yīng)值不斷減小，而樁身彎矩值則隨埋深先增大后減小。從樁身位移響應(yīng)來看，埋深18 cm以上樁身有明顯的位移響應(yīng)，之下的部分位移很小，這與樁基震害調(diào)查發(fā)現(xiàn)樁基破壞位移主要集中于樁的上段位置相一致。

5 結(jié)語

本文針對凍土環(huán)境中樁基礎(chǔ)在水平動載作用下的響應(yīng)特性進(jìn)行試驗及理論分析模型研究，所得主要結(jié)論如下：

(1) 基于研究樁-土-結(jié)構(gòu)相互作用的動力BNWF模型，提出改進(jìn)的凍土地基中樁基礎(chǔ)動力分析模型。該模型可綜合考慮樁側(cè)凍土的水平非線性作用、樁側(cè)與土體間的非線性摩擦作用、樁尖土的非線性擠壓與提離作用以及遠(yuǎn)場土體阻尼對樁基動力特性的影響。

(2) 提出模擬樁側(cè)凍土水平非線性作用的雙向無拉力三段屈服彈簧單元及其參數(shù)的確定方法，可較好地分析樁-土接觸分離及凍土材料的非線性對樁基動力特性的影響。

(3) 與水平動載作用下模型樁基試驗結(jié)果對比，表明本文提出的凍土-樁計算模型能有效地考慮樁-土間的擠壓與提離效應(yīng)時樁基的動力特性，具有較好的適用性。

References)

[1] Suleiman M T，Sritharan S，White D J.Cyclic Lateral Load Response of Bridge Column-foundation-soil Systems in Freezing Conditions[J].Journal of Structural Engineering，2006，132(11)：1745-1754.

[2] Sritharan S，Suleiman M T，White D J.Effects of Seasonal Freezing on Bridge Column-foundation-soil Interaction and Their Implications[J].Earthquake Spectra，2007， 23(1)：199-222.

[3] Wotherspoona L M，Sritharan S，Pender M J.Modelling the Response of Cyclically Loaded Bridge Columns Embedded in Warm and Seasonally Frozen Soils[J].Engineering Structures，2010，32(4)：933-943.

[4] Vaziri H，Han Y C.Full-scale Field Studies of the Dynamic Response of Piles Embedded in Partially Frozen Soils[J].Canadian Geotechnical Journal，1991，28(5)：708-718.

[5] Han Y C，Vaziri H H，Hubble D.Dynamic Response of Pile Foundations in Frozen Soils[C]//Canadian Geo-technical Conference.1991，40：1-7.

[6] 吳志堅，王平，霍元坤，等.多年凍土區(qū)橋梁樁基礎(chǔ)地震響應(yīng)的模型振動試驗研究[J].西北地震學(xué)報，2009，31(4)：319-326.WU Zhi-jian，WANG Ping，HUO Yuan-kun，et al.Study on Shaking TableTest for Seismic Response of Pile Foundation of Bridges at the Permafrost Regions[J].Northwest Seismological Journal，2009，31(4)：319-326.(in Chinese)

[7] 李永波，張鴻儒，全克江，等.凍融條件下模型樁基水平動力試驗研究[J].巖土力學(xué)，2012，33(2)：433-438.LI Yong-bo，ZHANG Hong-ru，QUAN Ke-jiang，et al.Experimental Study of Model Pile Foundations Under Lateral Dynamic Load in Frozen and Thawed Soils[J].Rock and Soil Mechanics，2012，33(2)：433-438.(in Chinese)

[8] 劉曉強(qiáng)，趙楠，夏開文，等.循環(huán)荷載作用下樁周土體強(qiáng)度弱化特性試驗研究[J].地震工程學(xué)報，2014， 36(4)：987-990.LIU Xiao-qiang，ZHAO Nan，XIA Kai-wen，et al.Test and Study on Strength Weakening Behavior of Soil Around a Pile Under Cyclic Loading[J].China Earthquak Engineering Journal，2014，36(4)：987-990.(in Chinese)

[9] El Naggar M H，Shayanfar M A，Kimiaei M，et al.Simplified BNWF Model for Nonlinear Seismic Response Analysis of Offshore Piles With Nonlinear Input Ground Motion Analysis[J].Canadian Geotechnical Journal，2005，42(2)：365-380.

[10] Allotey N，El Naggar M H.Generalized Dynamic Winkler Model for Nonlinear Soil-structure Interaction Analysis[J].Canadian Geotechnical Journal，2008，45(4)：560-573.

[11] Allotery N，F(xiàn)oschi R.Coupledp-yt-zAnalysis of Single Piles in Cohesionless Soil Under Vertical and /or Horizontal Ground Motion[J].Journal of Earthquake Engineering，2005，9(6)：755-775.

[12] 孫利民，張晨南，范立礎(chǔ).橋梁樁土相互作用的集中質(zhì)量模型及參數(shù)確定[J].同濟(jì)大學(xué)學(xué)報，2002，30(4)：409-415.SUN Li-min，ZHANG Chen-nan，PAN Long，et al.Lumped-mass Model and Its Parameters for Dynamic Analysis of Bridge Pier-pile-soil System[J].Journal of Tongji University，2002，30(4)：409-415.(in Chinese)

[13] Kagawa T，Kraft J，Leland M.Seismicp-yResponses of Flexible Piles[J].Journal of the Geotechnical Engineering Division，1980，106(8)：899-918.

[14] El Naggar M，Novak M.Effect of Foundation Nonlinearity on Modal Properties of Offshore Towers[J].Journal of Geotechnical Engineering，1995，121(9)：660-668.

[15] Boulanger R W，Curras C J，Kutter B L，et al.Seismic Soil-pile-structure Interaction Experiments and Analyses[J].Journal of Geotechnical and Geoenvironmental Engineering，1999，125(9)：750-759.

[16] 李永波.考慮凍土-樁動力相互作用的長大橋梁地震響應(yīng)分析[D].北京：北京交通大學(xué)，2013.LI Yong-bo.Seismic Response Analysis of Grand Bridge Including Dynamic Frozen Soil-pile Interaction[D].Beijing: Beijing Jiaotong University，2013.(in Chinese)

[17] 伍小平.砂土-樁-結(jié)構(gòu)相互作用振動臺試驗研究[D].上海：同濟(jì)大學(xué)，2002.WU Xiao-ping.Shaking TableTest of Sand-pile-structure Interaction System[D].Shanghai：Tongji University，2002.(in Chinese)

[18] 日本鐵道綜合技術(shù)研究所.鐵道構(gòu)造物等設(shè)計標(biāo)準(zhǔn)及解說V抗震設(shè)計[S].2000.Japan Railway Technical Research Institute.Seismic Design standards and Explanations for Railway Structures[S].2000.(in Chinese)

Study on Frozen Soil Pile Interaction Model Based on Dynamic Beam on Nonlinear Winkler Foundation Method

LI Yong-bo1， ZHANG Hong-ru2

(1.TianjinMunicipalEngineeringDesign&ResearchInstitute，Tianjin300050，China；2.SchoolofCivilEngineering，BeijingJiaotongUniversity，Beijing100044，China)

Permafrost is distributed widely in our country.Although pile foundation is the main foundation form for large bridges in frozen zones，few studies have been conducted regarding the dynamic characteristics of pile foundation construction in permafrost.An improved nonlinear analytical model for frozen soil-pile dynamic interaction was developed on the basis of model test results of the dynamic characteristics of pile foundations in various subzero temperature soils under lateral dynamic loads and an analytical model of pile-soil-structure interaction with the dynamic beam on nonlinear Winkler foundation model.This improved analytical model simulates the nonlinear lateral pressure effect of the frozen soil around the pile with an improved bidirectional compression-only multi-yield spring.Vertical friction effects between the frozen soil and pile, extrusion and separation effects of frozen soil under pile tips, and damping effects of the far-field soil on the dynamic characteristics of pile foundations are also considered in this analytical model.The parameters of the bidirectional tensionless multi-section yield spring are determined by the nonlinearp-yrelationship of frozen soil.Thep-yrelationship whose parameters were based on the results of the indoor frozen soil compression test was simulated by the combination of a cubic function and a constant function.Displacement-force response of the pile head and pile shaft responses to bending moments at different depths under dynamic forces were very similar to the model test results.Results indicate that the improved analytical model is potentially helpful during the analysis of frozen soil-pile dynamic interactions.

frozen soil pile interaction； model test； foundation nonlinearity； dynamic beam on nonlinear Winkler foundation (BNWF) model； bidirectional tensionless multi-section yield spring

2014-08-20

國家自然科學(xué)基金(90715013)

李永波(1985-)，男，山東東明人，博士，主要從事凍土-樁動力相互作用、長大橋梁及地下空間結(jié)構(gòu)抗震等方面的研究工作.E-mail: liyb2011@foxmail.com

P642.14；TU473

A

1000-0844(2015)02-0453-007

10.3969/j.issn.1000-0844.2015.02.0453