桩土软支撑）结构非线性动力相互作用分析

熊辉　杨佳

摘要：基于群桩土软支撑模型，建立了土结构相互作用有限元简化分析模型.针对不同地震激励，不同桩土条件下对模型进行了动力非线性时程分析.结果表明：在某些地震动和土基础条件下，上部结构非线性作用动力反应结果可能大于固基假定情形，且下部桩土软支撑会对上部结构柔弱层位置产生影响，分析结果与现有的试验成果具有良好的一致性.以简化非线性弹阻单元的形式等代桩土实体单元接触面，能有效和快速地进行复杂的SSI动力时程分析及抗震评估.

关键词：群桩土软支撑模型；土结构相互作用；非线性时程；动力反应

中图分类号：TU473 文献标识码：A

Abstract：Based on the pilesoil soft support model， a simplified interaction finite element analysis model of soil structure was established， and the nonlinear dynamic time history analysis was conducted for this model under different seismic excitations and different conditions of the pilesoil. The result shows that， under certain ground motion and soilfoundation conditions， the nonlinear effect of the upper structure can be greater than fixed base assumption， and pilegroupsoil soft support model also influences the weak layer position on the upper structure. The analysis results and the existing test results are in good agreement. By using structural units to solve complex computational problems in geotechnical engineering， designers can carry out dynamic time history analysis and seismic evaluation for the superstructure quickly， accurately and efficiently.

Key words：pilegroupssoil soft support model； soilstructure interaction； nonlinear time history； dynamic response

传统抗震设计的刚性地基假定分析不考虑桩土协同作用，但对于现在普遍采用桩基础的高层建筑，尤其在软土地基条件下，由于桩土上部结构体系的相互作用，长周期的地震动放大效应可能会对周期较长的高层建筑产生严重的影响.同时地震动作用下，桩土交界面上发生的脱开再闭合现象，如何快速而较为准确的进行桩土结构动力相互作用分析成为了当前研究的难点和热点.任重翠等[1]采用双向线性弹簧模拟桩土相互作用，并且考虑了结构的二阶效应及材料弹塑性等因素，研究了独柱高架车站考虑桩土相互作用的抗震性能；张亚旭等[2]利用ABAQUS建立了土桩框架结构非线性相互作用的有限元精细模型，采用接触面对法进行桩土界面的模拟，认为桩与桩周土的碰撞可能会导致桩基破坏.本文所建立的平面桩筏高层结构相互作用体系在考虑桩土桩相互作用的基础上，将桩周土分为近域和远域，在近域场提出了软支撑模型，即采用带有滑动元件的非线性弹阻单元来模拟桩土接触界面的非连续变形现象，从而反映了桩土基础的柔性参与以及桩土动力作用.

1分析模型

1.1整体有限元模型

基于桩基结构的复杂性，尤其在涉及桩土接触非线性问题时，若按照岩土实体单元划分，即便是较简单的问题桩土结构体系相互作用的地震反应计算量也是非常大的[3]，因此桩土结构体系模型规模和复杂程度不宜太高，故本文采用如图1所示的二维平面整体有限元模型来模拟桩土剪切型框架（实际上，上部结构可根据需要选择不同体型）的动力特性：框架为二结点梁元，其梁的刚度值远较柱刚度为大，以满足剪切型结构特点；对于桩身同样采用梁单元进行模拟，其范围内的土层采取“人工划分”.

桩间土由质量弹簧单元模拟，combin40非线性单元模拟桩周近域土抗力，远域土为串联式弹阻质量单元，且由结构自振分析结果截取得到相应的参考频率阻抗作为桩间土、桩周近、远域土弹簧阻尼单元的输入特性参数.该过程利用通用有限元Ansys软件APDL过程来实现.

1.2上、下部动力共同作用基本运动方程

依达郎贝尔定理，根据结构以及联系结构与桩基的承台基础发生平动及摇摆反应时的瞬时平衡状态，可列出仅考虑垂直入射地震SH波且上部结构为纯剪切形结构时上、下部相互作用平衡方程的矩阵形式：

2桩土软支撑模型

本文桩采用梁单元模拟，质量单元考虑土体的参振特征，引入具有桩桩、桩土、土土、筏土阻抗特性的弹阻单元来描述不同群桩的布置、土层状况因素对体系反应的参与作用，如图2所示.考虑到动力分析时桩周弱化土域的变形非线性行为、桩土界面在大位移条件下的相对滑动和相对分离等非连续现象以及其计算区域的滞回阻尼效应，本文在[4]的基础上给出了一种具有一致质量的双质点耦联弹簧形式的近域非线性文克尔模型，即软支撑弹阻单元，由非线性弹簧、阻尼器以及滑动块组成.该单元组集刚度及质量可表示为：

（a）桩土相互作用示意图

（b）桩间弹簧质量阻尼单元

（c）近域土的非线性弹阻单元

（d）线性远域土弹阻模拟

3.1时程曲线

为了考查上下部动力共同作用中的刚体（运动）相互作用效应，本文通过SFD波输入来反映承台与邻近自由场运动.图3，图4分别为SFD波激励下固基假定、线性相互作用及考虑桩土接触非线性时的结构顶层（A点）、桩头与承台连接点（B点）以及桩身（C点）的位移和加速度时程曲线.从反应结果来看，线性情形下B和C两点的位移反应几乎是重合的；对于桩土间发生了非线性滑移的情形，B和C的位移和加速度峰值较线性情形高出2～3倍，刚体相互作用明显.同时可知：1）由于桩土软支撑对结构动力反应的柔性参与作用，考虑相互作用的体系顶点位移幅值大于固基假定的情形，这与相互作用理论的普遍认识一致；2）非线性与线性位移反应差异不大，说明软支撑弹阻单元非线性对体系的顶点位移反应贡献较小.

3.2各层加速度及位移峰值

各层的峰值位移均随楼层增大而增大呈剪切型变化（见表1）.对本模型而言，刚性假定情况下位移峰值均小于考虑相互作用分析情况，且非线性分析结果大于线性分析结果.由表1可知，刚性假定情况下各层加速度峰值总体大于考虑相互作用分析情况，但在考虑桩土接触非线性情况时，PAK波作用下出现了加速度峰值大于刚性假定的情况（工况2顶层加速度峰值增大了6.7%）；上部结构各层加速度峰值总体呈现上大下小的“剪切型”趋势，但在SFD波输入的工况1条件下出现了中部小而上部和底部大的“弯剪型”变化.原因可能是在地震作用下，桩土接触面会发生接触分开接触的过程，这影响了结构的振动特性，不同激励及不同工况下土层对地震动可能出现放大效果，但主要还是体现了软支撑的减震作用.

3.3层间位移、层间剪力

以考虑相互作用结构的内力和变形与固基假定情况的比值K作为分析指标，考察在不同地震作用，不同工况条件下，考虑相互作用线性分析与非线性分析中框架结构的层间位移、剪力的变化情况，结果如图5，图6所示.层间位移K值表现为SFD波作用下非线性分析对工况1影响较大，且非线性结果大于线性结果，Pak波作用下上部结构非线性分析结果略小于线性分析结果；对于工况2条件下，PAK波作用下非线性分析影响不大，而SFD波输入下部结构非线性结果小于线性分析结果.

3.4软支撑模型柔性贡献分析

为了进一步反映桩土（软支撑）结构非线性 SSI效应对上部结构的柔性贡献，本文还取桩土组合采用R0=0.15 m， Gs1=1 000 kPa的情形作为工况3进行了分析，模态分析得到PGSS体系前3周期与固基周期比值Ti/T （i=1～3）分别为2.564，2.096，1.705，对比上述工况1和工况2的结果可知与文献[9-11]中尚守平等进行的上部结构与地基相对刚度比对土结构体系基频影响试验结论是一致的：基频折减率和上部结构与地基相对刚度比有关， 相对刚度比越大基频折减率越大[12].

通过时程分析可知工况3情形下结构顶部在非线性分析、线性分析下最大位移分别为PAK波作用下7.06 cm，6.84 cm；SFD波作用下5.98 cm，6.13 cm.同样结构顶层加速度分别为PAK波作用下2.03 m/s2，2.46 m/s2（0.825∶1.0）；SFD波作用下1.23 m/s2，1.51 m/s2（0.814∶1.0）.对比工况1，工况2可知下部越软，桩土（软支撑）结构非线性 SSI效应对上部结构的柔性贡献越大，符合吕西林等[3]分层土基础高层框架结构相互作用体系振动台模型试验分析结果：土体软化，非线性发展加强，土体传递振动能力减弱，上部结构加速度反应小.

4结论

通过对考虑桩土接触非线性的桩土（软支撑）框架结构相互作用体系的地震反应分析，得到一些结论：

1）利用各种具有软支撑阻抗特征的弹阻单元作为桩桩、桩周土以及土筏基础之间的一种反馈力边界来模拟SSI体系动力反应较岩土实体划分单元计算效率成倍提高；

2）考虑桩土结构体系相互作用数值分析的上部结构效应明显不同于固基假定的结果.SFD波作用激励下，非线性情形下刚体相互作用效应显著于线性情形，且在较软土及较小桩径组合工况中，考虑软支撑弹阻效应可改变上部结构的柔弱层部位，使结构各层加速度峰值变化趋势由“剪切型”转为“弯剪型”；

3）地震动作用过程中桩与土界面处出现的相对滑移、脱开及闭合等非线性行为，使得藕联系统的振动特性发生了根本的变化，不同的地震波作用，不同的场地条件，结构的反应无一定的规律，系统变化后的频率若和地震波主频接近，则土层对地震动起到放大效果，反之则会起到隔震作用，需要“一对一”的工程验算分析；

4）针对桩土高层框架进行的地震动作用分析结果与现有的试验成果有良好的一致性，说明土体软支撑的PGSS体系非线性动力分析方法具有相当的理论和工程实践意义.

参考文献

[1]任重翠，徐自国，肖从真，等.独柱高架车站考虑桩土相互作用的抗震性能分析[J].建筑结构学报，2014，35（5）：25-32.

REN Chongcui， XU Ziguo， XIAO Congzhen，et al.Seismic performance analysis of singlecolumn elevated station considering pilesoil interaction[J]. Journal of Building Structures， 2014，35（5）：25-32.（In Chinese）

[2]张亚旭，王修信，庄海洋.土桩框架结构非线性相互作用的精细数值模型及其验证[J].防灾减灾工程学报，2010，30（5）：558-566.

ZHANG Yaxun， WANG Xiuxin， ZHUANG Haiyang.Fine numerical modeling of nonlinear soilpileframe structure interaction system and its verification[J]. Journal of Disaster Prevention and Mitigation Engineering， 2010，30（5）：558-566.（In Chinese）

[3]陈跃庆，吕西林，李培振.分层土基础高层框架结构相互作用体系振动台模型试验研究[J].地震工程与工程振动，2001，21（3）：104-112.

CHEN Yueqing， LV Xilin， LI Peizhen. Shaking table testing for layered soilfoundationstructure interaction system[J]. Eathquake Engineering and Engineering Vibration， 2001，21（3）：104-112.（In Chinese）

[4]熊辉，邹银生.层状土中考虑频域内轴向力分担的群桩水平动力阻抗[J].计算力学学报，2004，23（4）：757-763.

XIONG Hui， ZOU Yinsheng. Horizontal vibration impedance of pilegroups considering allocated frequencydependent axial forces in layered soilfield[J]. Journal of Computational Mechanics， 2004， 23（4）：757-763.（In Chinese）

[5]WANG S， KUTTER B L，EERI M， et al.Nonlinear seismic soilpile structure interaction[J]. Earthquake Spectra， 1998，14（2）：377-396.

[6]NOGAMI T， OTANI J， KONAGAI K. Nonlinear soilpile interaction model for dynamic lateral motion[J]. Journal of Geotechnical Engineering， ASCE，1992，118（1）：89-106.

[7]尚守平，贺志文，王海东，等.上部结构与地基相对刚度比对土结构体系基频影响试验研究[J].地震工程与工程振动，2008，28（5）：94-101.

SHANG Shouping， HE Zhiwen， WANG Haidong， et al. Experimental investigation on the effect of the relative stiffness ratio between superstructure and ground soil on the fundamental frequency of soilstructure system[J].Journal of Earthquake Engineering and Engineering Vibration， 2008，28（5）：94-101.（In Chinese）

[8]熊辉.群桩土结构（PGSS）相互作用体系的动力有限元方法研究[D].上海：同济大学结构工程与防灾研究所，2005：91-134.

XIONG Hui. Dynamic finiteelement method research on pilegroupedsoilstructure（PGSS）interaction system [D]. Shanghai：Research Institute of Structural Engineering and Disaster Reduction，Tongji University，2005：91-134.（In Chinese）

[9]尚守平，陈婉若，鲁华伟，等.土筏基钢框架结构动力相互作用试验研究[J].湖南大学学报：自然科学版，2013，40（3）：1-5.

SHANG Shouping， CHEN Wanruo， LU Huawei，et al. Experimental investigation of the dynamic soilinteraction about steel frameraft foundation model[J].Journal of Hunan University： Natural Sciences， 2013，40（3）：1-5.（In Chinese）

[10]尹华伟，易伟建.简化桩模型在土结构相互作用分析中的应用[J].湖南大学学报：自然科学版，2007，34（2）：11-15.

YIN Huawei， YI Weijian. Application of simplified pile model in soil structure interaction analysis[J].Journal of Hunan University： Natural Sciences， 2007，34（2）：11-15.（In Chinese）

[11]尚守平，朱志辉，涂文戈. 土箱形基础结构动力相互作用的模态试验分析[J].湖南大学学报：自然科学版， 2004， 31（5）：71-76.

SHANG Shouping， ZHU Zhihui， TU Wenge. Modal experimental analysis for soilbox foundationstructure dynamic interaction system [J]. Journal of Hunan University：Natural Sciences， 2004， 31（5）： 71-76.（In Chinese）

[12]周峰，林树枝.实现桩土共同作用的机理及若干方法[J].建筑结构，2012，42（3）：140-143.

ZHOU Feng， LIN Shuzhi. Working mechanism and some realization methods of pilesoil interaction[J].Building Structure， 2012，42（3）：140-143.（In Chinese）