土-结构动力相互作用研究综述 (1)
桩—土—结构动力相互作用体系振动台模型试验研究

土-箱基-框架结构动力相互作用体系试验研究

a tr tc n ce sd t en t r I e ido h y tm .Th i l u e as aiae h E n lss ce i i a d i ra e h au a ro ft es se s s n p esmi td l v l td t eF a ay i. i o d Ke r s d n i; x f u d t n;b scfe u n y rw n ees ywo d : y a c b o n a i m o o ai r e c ;d a a dr l e q a
S NG h upn t Z HA S o —ig 。 HAN J sn 。 G i o g WANG l o g HU h—u — Ha d n .Z — Zi i h
(ol eo il n ier g Hu a i。 h n sa Hu a 4 0 8 。C ia C lg f v gnei 。 n Unv C agh , nn 10 2 hn ) e C iE n n
t n o h t cu eb x fu d t n s i s se .wh c dc tdt a h tr cin c a e h yai - l y tm o u o o ih i iae h tt ei e a t h n d t ed n m c h r n n o g
adFnt l n( E)aa s f1'l g a d1 ycmpr gtebs rqec ft tadF n ii Ee tF e me nl io / a es lmo e.B o ai h ai f uni o n E ys 2 r ce n c e s e s e
维普资讯
第3 4卷
第 1期
湖 南 大 学 学 报 (自 然 科 学 版 )
一种三维饱和土--基础--结构动力相互作用分析方法

一种三维饱和土--基础--结构动力相互作用分析方法陈少林;赵宇昕【摘要】地震波入射时土与结构动力相互作用分析是地震工程领域的重要问题。
由于问题的复杂性,以往的研究大多考虑地基土为干土情形。
而实际工程中,土体中经常充满孔隙水,土层往往是含水层或部分含水层。
孔隙水对土层的地震反应影响较大,进而影响支撑于其上的基础和上部结构的响应。
因此,有必要考虑饱和土-基础-结构动力相互作用问题。
基于此,土体采用Biot模型,利用集中质量显式有限元方法并结合多次透射人工边界进行模拟;结构经有限元离散后,采用Newmark隐式时步积分方法进行分析,可通过ANSYS等有限元软件实现;基础假定为刚性,采用显式时步积分进行求解;土体和结构(及基础)可分别采用不同的时间步距;通过FORTRAN程序实现了三维饱和土-基础-结构系统在地震作用下的整体分析。
从饱和多孔介质动力微分方程出发可知,干土是饱和土的特殊情形,通过将流体体积模量及孔隙率置为零,可将饱和土退化到干土,从而将干土统一到饱和土的计算框架中,通过土-结构相互作用算例对此进行了验证,进一步实现了干土与饱和土混合情形时的土-结构动力相互作用分析,使得问题的模拟更贴近实际工程(如考虑地下水位情形)。
通过算例对比了饱和土地基、干土地基、干土覆盖饱和土地基(考虑地下水位)情形时,土-结构相互作用对基础和结构响应的影响,结果表明地下水位对基础和结构响应的影响较大。
%Analysis of soil-structure dynamic interaction subjected to seismic wave is a key problem in earthquake engineering. In general, the soil stratum consists of two-phase saturated porous zones and single-phase viscoelastic zones due to ground water. In most cases of soil-structure dynamic interaction analysis, the soil has been assumed to be a single-phase viscoelastic medium for simplicityand little attention has been paid to the saturated porous soil case, even less to case of the viscoelastic soil layered on saturated soil. In this study, an efficient method for three-dimensional saturated soil-foundation-structure dynamic interaction analysis is proposed. The unbounded saturated soil is modelled by lumped-mass explicit finite element method and transmitting boundary condition the structure is analysed through implicit finite element method, and response of the rigid foundation is calculated through explicit time integration scheme. The different time steps can be chosen for the explicit and implicit integration scheme, which can greatly improve the efficiency. In addition, based on the fact that single-phase soil is a special case of two-phase saturated soil, the dynamic analysis of single-phase soil can be unified into the analysis of saturated soil by setting the bulk modulus of pore fluid and porosity to be zero. Thus, the soil-structure interaction analysis for the viscoelastic soil layered on saturated soil case is realized, which can approximate the ground water table in practice. The effects of ground water on the response of foundation and structure are examined through numerical examples of soil-structure interaction analysis for saturated soil, single-phase viscoelastic soil and the viscoelastic soil layered on saturated soil, and the results show that the ground water influences the structure and foundation responses greatly.【期刊名称】《力学学报》【年(卷),期】2016(048)006【总页数】10页(P1362-1371)【关键词】土-结构动力相互作用;饱和多孔介质;显-隐式积分方法;人工边界条件;地震响应分析【作者】陈少林;赵宇昕【作者单位】南京航空航天大学土木工程系,南京210016;南京航空航天大学土木工程系,南京210016【正文语种】中文【中图分类】TU435土-结构相互作用是地震工程领域的重要问题,以往对该问题的研究集中在地基土为干土的情形,对地基土为饱和土情形的研究相对较少.而实际情形中,地下水位以上可视为干土,水位以下可视为饱和土,对该情形的研究更是鲜有报导.地震波输入下土-结构相互作用分析一般采用子结构法[1-4]、集中参数法[5-8]和直接法[9-12].子结构法和集中参数法将土-结构相互作用问题分解为如下可独立求解的子问题:(1)自由场问题,即结构和基础不存在时土层在地震波入射时的响应问题;(2)基础的输入问题,即结构不存在时,无质量刚性基础在地震波入射时的响应问题;(3)刚性基础的动力刚度或动力柔度矩阵;(4)结构响应问题.下面仅就饱和土层自由场分析和饱和地基动力刚度研究的相关工作做一简单介绍.自 Biot[13-15]建立多孔弹性介质波的传播理论以来,Deresiewicz等[16-18]对于波在饱和土界面上的传播特性进行了系统的研究,分析了平面波在自由表面、分层多孔介质界面中的反射和透射.Jocker等[19]通过推广 Thomson-Haskell传递矩阵法,以得到一种封闭形式的解析表达式,用于计算水平成层多孔介质的反射率和透射率.Rajapakse等[20]、Degrande等[21]、Liang等[22]分别给出了饱和土层和半空间的精确动力刚度矩阵,采用直接刚度法可求得饱和水平层状场地的自由场响应.李伟华等[23]采用波函数展开法,求得饱和弹性均匀半空间的自由场,进一步分析了饱和沉积谷场地对SV波的散射问题.赵宇昕等[24]采用传递矩阵法分析了成层饱和土层自由场响应,并探讨了数值结果中的非因果原因.1986年,Halpern等[25]通过对基础所覆盖的单元域上 Green函数数值积分给出了三维饱和弹性半空间表面上刚性板的动力柔度系数.Bougacha等[26]采用基础空间离散振动模态的有限元方法得到了饱和弹性土层上长条形基础和圆形基础的动力刚度.Japon等[27]应用边界元方法分析了长条形基础的动力刚度,并考虑了渗透力、附加密度、地层深度、基础刚度对基础振动的影响.陈少林等[28-29]提出了一种时域求解基础阻抗函数的方法,通过给基础输入脉冲位移,应用集中质量显示有限元方法结合局部透射人工边界,得到地基施加于基础的反力时程,然后根据阻抗函数定义,应用傅里叶变换求得基础阻抗函数.子结构法以频域内的动力刚度为基础,集中参数法通过拟合动力刚度得到参数模型,理论上都只适用于线性系统,不能考虑土体的非线性.直接法将地基、基础和上部结构一并计算,可以方便考虑土体的非均匀性和非线性,但计算量一般较大,尤其是饱和土情形.梁建文等[30-31]将单相土-结构相互作用的研究拓展到饱和土领域[32],采用间接边界元法分析了二维线性饱和土-结构相互作用问题,分析了P-SV 波斜入射时,孔隙率、土层深度、土层与基岩间的刚度比等参数的影响.本文拟将显-隐式结合的土-结构相互作用分析方法[10],推广到三维饱和土-结构动力相互作用情形,采用集中质量显示有限元结合透射人工边界的方法分析饱和土体[33-35],采用隐式方法对上部结构进行分析,实现三维饱和土-基础-结构动力相互作用分析.根据干土是饱和土的特殊情形这一事实,将两者情形的土-结构相互作用分析统一到饱和土-结构相互作用分析框架,进而可分析考虑地下水位情形的土-结构相互作用分析.1.1 土体内节点运动根据Biot[1]理论,饱和弹性多孔介质运动方程可以表示为其中u和U为固液相位移矢量,D,N,Q,R为非负弹性常数,ρ11=ρ1+ρα,ρ22=ρ2+ρα,ρ12=−ρα. ρ1=(1−φ)ρs,ρ2=φρw,ρs为固相质量密度,ρw为液相质量密度,ρα为固、液两相惯性耦合质量密度,φ为孔隙率.衰减系数为流体动力黏度,k0为渗透率.从地基土中取一有限土体,采用八节点三维实体元对其进行空间离散.将土体划分成两个区域,即人工边界区和内部计算区(如图1),将土体单元节点划分为人工边界点、与基础相连的点和内节点三类.对于内点采用集中质量显式有限单元法.在人工边界上,采用多次透射公式模拟无限域的影响.对于内节点而言,对方程(1)进行有限元空间离散,得到第i个节点的固相运动方程为[33]液相运动方程为其中,和为节点i(全局编号)的固相加速度矢量和液相加速度矢量,为单元e中节点j(单元局部编号)的固相位移和速度矢量,为单元e中节点j的液相位移和速度矢量.Msi和Mwi分别表示集中在i节点上的固相和液相质量阵;是考虑固液两相互相运动所产生的单元阻尼阵的子矩阵,下标k(i)表示节点i(全局编号)在单元中对应的局部节点编号为k,即全局节点编号和单元内的局部节点编号之间的对应关系.表示单元刚度阵的子矩阵,为单元e中分配给节点i的固相载荷矢量和液相载荷矢量.N为包含节点i的单元个数,J为单元e的节点个数.此外,我们可以注意到,上述阻尼矩阵仅考虑固液两相互相运动所产生的阻尼,对固相骨架考虑瑞雷阻尼那么式(2)可修正为对内节点方程(5)和方程(3)采用中心差分与单边结合的积分格式进行时域积分.假定在p+1及其以前时刻节点i的位移矢量已知,时间离散步长用∆t表示,那么p时刻的加速度和速度可以分别表示为其中W=u,U.上式代入式(5)和式(3)可得节点i的位移递推公式1.2 土体边界节点运动对于人工边界上的节点,其运动状态可通过多次透射公式[36]进行确定其中,为边界节点0在p+1时刻的外传波位移(散射位移),N为透射阶数,二项式可表达为为由人工边界节点及内部节点的散射位移构成的列向量行向量Tn为其中上式中记号“Σ”表示对所有满足以下条件的项求和以一阶Modulation Transfer Function(MTF)为例,那么可知其中S=ca∆t/∆x.式(9)分别用于饱和多孔介质的固相位移和液相位移,可得到边界点在 p+1时刻的散射波场位移,加上边界点的自由场位移,即可得到边界点在p+1时刻的总位移.式(11)中的散射位移向量可通过式(7)和式(8)得到的总位移减去自由场位移得到.对于饱和成层土体的自由场计算可参见文献[24].对于上覆干土,下卧饱和土的情形,可通过满足干土和饱和土分界面上的4个应力位移连续条件,得到相应的传递矩阵.具体计算流程可参见图2.1.3 上部结构的运动对于一般情况下的结构或系统,其动力学方程均可写成在p~p+1时步的间隔∆t内,Newmark积分格式对位移、速度、加速度采用如下的假设关系式中,β和γ为积分参数.因此将式(18)代入式(17)可以得到系统p+1时刻的运动平衡方程由式(19)可以得到up+1,进而通过式(18)可得到p+1时刻速度和加速度.1.4 基础的运动刚性基础作为土体和上部结构的连接部分,起到了承上启下的作用,传递结构与土体之间的作用力.其运动可由6个分量描述,即3个平动分量和3个转动分量土和结构作用在基础上的合力使基础产生刚体运动式中,MF为基础的集中质量阵,对角线元素依次为三方向质量MFx,MFy和MFz,三方向绕质心的转动惯量IFx,IFy和IFz;FS和FD分别为结构和土体作用于基础上的广义力矢量.首先考虑p时刻土体对刚性基础作用力饱和土体与基础接触任一节点k处,没有力矩的传递只有3个方向集中力的传递,这三方向集中力分别由每一时刻固液两相本构力合成因此饱和土体对刚性基础的合力可以表示为其中,m表示土体与基础接触点总数;表示节点k对基础作用力;和是土体对基础三方向的集中作用合力;是土体对基础三方向的转动合力矩;∆Xk,∆Yk和∆Zk是节点k对基础质心的相对坐标;矩阵A是一几何转换矩阵,代表土体接触点或结构接触点和基础质心的相对坐标关系.其次考虑p时刻结构对刚性基础作用力类似于上述的方法,对于任一结构与基础接触节点i处,上部结构对基础作用力可表示为可通过ANSYS分析得到各点反力值因此上部结构对刚性基础的合力可以表示为其中,n表示结构与基础接触点总数;和是结构对基础三方向的集中作用合力;和是结构对基础三方向的转动合力矩.得到p时刻作用力FD与FS后,对式(20)采用中心差分格式积分,可得通过刚性基础位移可得到与其接触的土体点或结构点的位移1.5 饱和土-基础-结构相互作用分析流程上面对饱和土-基础-结构相互作用分析的各个部分计算方法做了详细阐述,下面简述其基本分析步骤.假定p及以前时刻土体、基础与结构位移已知,计算p+1时刻各点的位移,基本步骤可以可描述如下:(1)采用集中质量显式有限元的方法,由式(7)和式(8)得到饱和土体内节点固、液相位移;(2)利用多次透射式(9)和自由场位移得到饱和土体人工边界节点的固、液相位移;(3)根据式(29)得到基础整体位移,进而由式(30)得到与刚性基础相连土体节点位移以及结构节点位移;(4)得到与基础相联的结构节点位移后,将该位移作为结构约束,根据式(19)可得到结构各点位移,进一步可得到该时刻结构对基础的作用力;(5)重复执行步骤(1)~(4),即可得到饱和土-基础-结构体系在各时刻的动力反应. 根据上述原理,通过编制相应的FORTRAN程序实现了地震波入射时饱和土-基础-结构相互作用分析.其中结构可通过ANSYS有限元软件进行分析,通过FORTRAN调用ANSYS实现饱和土-基础-结构系统的整体分析.当流体体积模量Kf以及孔隙率φ均置为零时,饱和土退化为干土[29].下面采用如下模型对上述退化进行数值验证.上部结构通过 ANSYS进行分析,采用悬臂梁模型,截面尺寸0.1m×0.1m,梁高2m,密度2500kg/m3,弹性模量3.00GPa,泊松比0.2,阻尼比0.05,剪切波速707m/s,一阶振型周期0.226s.刚性基础尺寸6m×6m×4m,密度2500kg/m3.考虑表1中前两类地基土,输入相同的单位脉冲SV波:时间步距dt=0.0002s,脉冲宽度t0=0.15s,步数n=1013(见图 3).土体划分单元尺寸均为1m×1m×1m,满足波动模拟精度要求.结构隐式分析的时间步距dT=0.002s,即FORTRAN每10步对ANSYS执行一次调用.图4和图5分别给出了基础和结构顶部的位移.实线为干土情形的结果,虚线为饱和土退化情形的结果,两者完全吻合.因此,地基土为干土情形的土-基础-结构相互作用分析可纳入到饱和土-基础-结构相互作用分析框架,进而可以分析考虑地下水位的干土覆盖饱和土情形的土-结构相互作用分析.另外,结构在基础有效输入为零后(因为结构柔且小,不足以带动基础做相应运动)做有阻尼的自由振动,其周期正好为结构一阶振型周期0.226s.实际工程中,地下水位位于地表以下一定深度,若直接考虑为干土或饱和土并不合适.以下就地震波入射时土-基础-结构动力相互作用问题,比较地基土分别为干土情形、饱和土情形、干土覆盖饱和土情形对基础和结构响应的影响.仍采用悬臂梁结构及刚性基础,首先考虑地基土为表1中的3号、4号、5号的情形,分别称为CASE1,CASE2,CASE3.其中,3号为饱和土,4号为具有与3号饱和土相同剪切波速和密度的干土,在工程中,常常对饱和土进行如此简化;5号为干土,具有与3号饱和土相同的骨架剪切波速和固相密度.最后考虑由上覆4m厚的5号干土,下卧14m厚的3号饱和土所组成的地基土,这种场地描述了当地下水位为4m时的情况,称其为CASE4.图6(a)和图6(b)分别为以上4种地基土情形时基础的位移时程和频谱.图7(a)和图7(b)为结构顶点的位移时程和频谱,为了比较,图中同时给出了不考虑土结相互作用情形时的结果.由图可知,就本文算例而言,若通过两种干土模型(CASE2,CASE3)代替饱和土模型(CASE1)进行分析,基础与结构的位移均将偏小.由于CASE3中的干土剪切波速更大,其到时也比CASE1与CASE2的情况要早.CASE4与前3种情况相比,基础与结构的位移都有较大幅度减小,并不是介于干土和饱和土情形的结果之间,表明一定深度的地下水位对基础和结构响应的影响较大.不考虑土结相互作用时,脉冲位移从结构底部输入,到时较其他情形早(见图7(a)),最大位移与饱和土-结构相互作用情形CASE1相当,但要大于其他情形(CASE2~CASE4).自由振动阶段的位移幅值小于CASE1,大于其余情形.另外,由图7可知,由于结构柔且小,在脉冲输入结束后,结构独自做有阻尼的自由振动(此时基础反应很小),此时土-结构相互作用影响较小,与土体情形基本无关,其周期均为结构一阶振型周期0.226s.土体中的孔隙水将影响土体的地震反应,进而影响支撑其上的基础和结构的响应,为此,本文发展了一种饱和土-基础-结构动力相互作用分析的高效方法,主要结论如下:(1)采用集中质量显式有限元方法结合多次透射人工边界模拟半无限饱和土体,利用隐式积分方法分析上部结构,实现了地震波输入时饱和土-基础-结构动力相互作用的高效分析.(2)通过将流体体积模量Kf,孔隙率φ置为零,将干土统一到饱和土计算框架,通过算例验证了这一方案,并实现了地基为干土覆盖饱和土情形的土-结构相互作用分析.(3)对比了饱和土地基、干土地基、干土覆盖饱和土地基(考虑地下水位)情形时,土-结构相互作用对基础和结构响应的影响,算例表明,地下水位对基础和结构响应的影响较大.【相关文献】1 Luco JE.Dynamic interaction of a shear wall with the soil.J Eng Mech DivASCE,1969,95:333-3462 Kobori T,Minai R,Shinozaki Y.Vibration of a rigid circular disc on an elastic half-space subjected to plane waves.Theoretical and Applied Mechanics,1973,(21):108-1193 Guin J,Banerjee P K.Coupled soil-pile-structure interaction analysis under seismic excitation.J Struct Eng Div ASCE,1998,124(4): 434-4444 宋亚新,蒋通,楼梦麟.桩基-非线性框剪结构相互作用体系(下)-相互作用对结构地震反应的影响.地震工程与工程振动, 2000,20(1):48-55(Song Yaxin,Jiang Tong,Lou Menglin.Pile-nonlinearframe-wall interaction systems.Part II:effects of interaction on seismic responses of structures.Journal of Earthquake Engineering and Engineering Vibration,2000,20(1):48-55(in Chinese))5 Wolf JP.Soil-structure interaction analysis in time domain.NJ: Prentice Hall,Englewood Cliffs,1988:11-356 Jean WY,Lin TW,Penzien J.System parameters of soil foundations for time domain dynamic analysis.Earthq.Earthquake Engineering and Structural Dynamics,1990,19(4):541-5537 王有为,王开顺.建筑物-桩-土相互作用地震反应分析的研究.建筑结构学报,1985,6(5):64-73(Wang Youwei,Wang Kaishun. Research on earthquake response analysis of building structure-pilesoil interaction.Journal of Building Structures,1985,6(5):64-73 (in Chinese)) 8 李永梅,孙国富,王松涛等.桩-土-杆系结构动力相互作用.建筑结构学报,2002,23(1):75-81(Li Yongmei,Sun Guofu,Wang Songtao,et al.Dynamic interaction of pile-soil-frame structure. Journal of Building Structures,2002,23(1):75-81(in Chinese))9 刘晶波,吕彦东.结构-地基动力相互作用问题分析的一种直接方法.土木工程学报,1998,31(3):55-64(Liu Jingbo,Lu Yandong. A direct method for analysis of dynamic soil-structure interaction. China Civil Engineering Journal,1998,31(3):55-64(in Chinese))10 陈少林,唐敢,刘启方等.三维土-结构动力相互作用的一种时域直接分析方法.地震工程与工程振动,2010,30(2):24-31(Chen Shaolin,Tang Gan,Liu Qifang,et al.A direct time-domain method for analysis of three-dimensional soil-structure dynamic interaction. Jounal of Earthquake Engineering and Engineering Vibration,2010, 30(2):24-31(in Chinese))11 杨笑梅,郭达文,杨柏坡.三维土-结构动力相互作用体系分析的两步时域显式波动有限元过程.地震工程与工程振动,2011, 31(4):9-17(Yang Xiaomei,Guo Dawen,Yang Baipo.A two-step explicit finit element method in time domain for 3D dynamic SSI analysis.Jounal of Earthquake Engineering and Engineering Vibration,2011,31(4):9-17(in Chinese))12 姜忻良,张海顺.土-结构非线性相互作用混合约束模态实施方法.振动与冲击,2015,34(6):52-56(Jiang Xinliang,Zhang Haishun. Mixed constraint modal method for nonlinear soilstructure interaction.Journal of Vibration and Shock,2015,34(6): 52-56(in Chinese))13 Biot MA.Theory of propagation of elastic waves in a fluid-saturate porous solid.Acoust Soc Am,1956,28:168-19114 Biot MA,Willis DG.The elastic coefficients of the theory of consolidation.Journal of Applied Mechanics,1957,24:594-60115 Biot MA.Mechanics of deformation and acoustic propagation in porous media.Journal of Applied Physics,1962,33(4):1482-149816 Deresiewicz H,Rice JT.The effect of boundaries on wave propagation in a liquid-fille porous solid:III.Reflectio of plane waves at a free plane boundary(general case).Bull Seism Soc Am,1962,52(3): 595-62517 Deresiewicz H,Rice JT.The effect of boundaries on wave propagation in a liquid-fille porous solid:V.Transmission across a plane interface.Bull Seism Soc Am,1964,54(1):409-41618 Deresiewicz H,Levy A.The effect of boundaries on wave propagation in a liquid-filleporous solid:X.Transmission through a stratifie medium.Bull Seism SocAm,1967,57(3):381-39119 Jocker J,Smeulders D,Drijkoningen G,et al.Matrix propagator method for layered porous media:analytical expressions and stability criteria.Geophysics,2004,69(4):1071-108120 Rajapakse RKND,Senjuntichai T.Dynamic response of a multilayered poroelastic medium.Earthquake Engineering and Structure Dynamic,1995,24(5):703-72221 Degrande G,Roeck G,Broeck P,et al.Wave propagation in layered dry,saturated and unsaturated poroelastic media.Int J Solids Structures,1998,35(34-35):4753-477822 Liang JW,You HB.Dynamic stiffness matrix of a poroelastic multilayered site and its Green’s functions.Earthquake Engineering and Engineering Vibration,2004,3(2):273-282 23 李伟华,赵成刚.饱和土沉积谷场地对平面SV波的散射问题的解析解.地球物理学报,2004,47(5):911-919(Li Weihua,Zhao Chenggang.Scattering of Plane SV waves by circular-arc alluvial valleys with saturated soil deposits.Chinese J Geophys,2004,47(5): 911-919(in Chinese))24 赵宇昕,陈少林.关于传递矩阵法分析饱和成层介质响应问题的讨论.力学学报,2016,48(5):1145-1158(Zhao Yuxin,Chen Shaolin.Discussion on the matrix propagator method to analyze the response of saturated layered media.Chinese Journal of Theoretical and Applied Mechanics,2016,48(5):1145-1158(in Chinese))25 Halpern MR,Christiano P.Steady-state harmonic response of a rigid plate bearing on a liquid-saturated poroelastic half-space.Earthquake Engineering and Structural Dynamics,1986,14:439-45426 Bougacha SJ,Roesset M,Tassoulas JL.Dynamic stiffness of foundation on flui fille poroelastic stratum.Journal of Engineering Mechechanics,ASCE,1993,119(8):1649-1662 27 Japon BR,Gallego R,Dominguez J.Dynamic stiffness of foundations on saturated poroelastic soils.Journal of Engineering Mechechanics,ASCE,1997,123(11):1121-112928 陈少林,甄澄.饱和地基上基础动力阻抗函数的一种分析方法.力学学报,2012,44(2):393-400(Chen Shaolin,Zhen Cheng.Dynamic impedance of foundation on saturated poroelastic soil.Chinese Journal of Theoretical and Applied Mechanics,2012,44(2): 393-400(in Chinese))29 陈少林,甄澄.下卧饱和半空间粘弹性土层上基础阻抗函数的分析.振动工程学报,2012,25(4):446-452(Chen Shaolin,Zhen Cheng.Dynamic impedance of foundations on viscoelastic stratum and saturated elastic half-space.Journal of Vibration Engineering, 2012,25(4):446-452(inChinese))30 Liang JW,Fu J,Todorovska MI,et al.Effects of the site dynamic characteristics on soil-structure interaction(I):incident SH waves. Soil Dynamics and EarthquakeEngineering,2013,44:27-3731 Liang JW,Fu J,Todorovska MI,et al.Effects of the site dynamic characteristics on soil-structure interaction(II):incident P and SV waves.Soil Dynamics and Earthquake Engineering,2013,51:58-7632 Liang JW,Fu J,Todorovska MI,et al.In-plane soil-structure interaction in layered,fluid-saturated poroelastic half-space I:Structuralresponse.Soil Dynamics and Earthquake Engineering,2016,81: 84-11133 陈少林,廖振鹏,陈进.两相介质近场波动模拟的一种解耦有限元方法.地球物理学报,2005,48(4):909-917(Chen Shaolin,Liao Zhenpeng,Chen Jin.A decoupling FEM for simulating near-fiel wave motions in two-phase media.Chinese J Geophys,2005,84(4):909-917(in Chinese))34 赵成刚,王进廷,史培新等.流体饱和两相多孔介质动力反应分析的显式有限元法.岩土工程学报,2001,23(2):178-182(Zhao Chenggang,Wang Jinting,Shi Peixin,et al.Dynamic analysis of fluid-saturate porous media by using explicit finit element method Chinese Journal of Geotechnical Engineering,2001,23(2):178-182 (in Chinese))35 王进廷,杜修力,赵成刚.液固两相饱和介质动力分析的一种显式有限元方法.岩石力学与工程学报,2002,21(8):1199-1204 (Wang Jinting,Du Xiuli,Zhao Chenggang.An explicit finit element method for dynamic analysis of fluid-saturate porous media. Chinese Journal of Rock Mechanics and Engineering,2002,21(8): 1199-1204(in Chinese))36 廖振鹏.工程波动理论导论(第2版).北京:科学出版社,2002. 1-298(LiaoZhenpeng.Introduction to Wave Motion Theories in Engineering(2nd Edn).Beijing:Science Press,2002.1-298(in Chinese))。
土-结构相互作用对隔震体系桥梁的影响

使得在两者的接触 面上 产生 了相 互作用力 , 进 而产生 了土与基 础
以至结构的相互作 用。
图 1 考虑土一结构相互作用的力学计算模型
2 桥 梁隔震体 系的分析 模型
. 1 桥 梁模 型。 。 土一结构动力相互作用 问题首 先要 确定 的是结 构 的动力 刚 2 度矩阵。对于简化的桥 梁模 型 , 可以将动力学方程 写为 J : 采用隔震支座后 , 上部 桥梁结 构可用 线 弹性模 型加 以描 述 ,
1 土一 结构 动力 相互作 用分 析原理
结构物与支撑 其的地基是一个共 同工作 的整体 系统 , 在地震 作用下 , 其土 与结构 在其 接触面有 着特 定的交 联关 系 , 称 为相 互 作用 。土一结构相互 作用 的实质就是 在地 震作用 下 由于土 体 与建筑物基础的材料差异引起 了两者变形能力 的差异 , 这种差 异
铅 芯橡胶支座并且考虑土一结构相互作用的连续梁桥的三种情况进行对 比分 析 , 结 果表 明土一结构相 互作 用对结构 的动力 响应 产 生不可忽 略的影响。
关键词 : 隔震体系 , 连续梁桥 , 实体模型 , 土一 结构相互作 用
中图分类号 : U 4 4 1 . 3 文献标识码 : A
的不合理。因此 , 在研究 桥梁 的地震 响应规 律时 , 关键 是要 分析
一
+
0
一
…
0
清楚 桥梁周 围场地 的地震动特性 , 即分析 土一结 构相互作 用对桥
梁结构的地震 动影响。
K=
0
一
。
‘
.
.
'
土一 结构体系共 同作用 的问题主要表现在三个方面 J :
收稿 日期 : 2 0 1 3 - 1 2 . 1 1
挤扩支盘灌注桩-土-上部结构体系动力相互作用的研究

结构 的安全 又可降低 造价 。目前桩基工程 中的一种新 型桩 ,
挤扩支盘灌 注桩 具有 良好 的抗压 抗拔性 能及 抗震性 能 ,在 工程界颇受青 睐。
Y i ,Q A e—l g EQ a n IND i n
( c o l f ii a d I ia o n ie r g S h o o vl n r g t nE gn ei ,Hee U i r t f e h ooyHe e,2 0 0 C r i n fi nv s y o T c n lg i i 3 0 9,C ia ei f hn )
S u n d na i n e a to fs i s r c ur y t m b v t dy o y m c i t r c i n o o l t u t e s s e a o e c nc e e po i ie o qu e e x n d s ppo tpa o r t urng p l fs e z d e pa de u r n
0 引 言
根 据 目前 已有 的 震 害 资 料 ,如 果 是 软 弱 的 地 基 土 ,柔 性 结构 破 坏 严 重 ,刚 性结 构 表 现 较 好 ,如 果 是 坚 硬 的地 基 , 柔 性 结 构 表 现 较 好 ,刚 性 结 构 的 表 现 有 的好 有 的 差 。 在 坚
为非挤土灌注桩 ( 普通灌 注 桩 ) 、部分挤 土灌 注桩 ( 挤扩 支 盘桩 ) 和挤 土灌 注桩 ( 沉管灌注桩 ) 三种 。 挤扩灌 注 桩 由 主桩 、底 盘 、中盘 、顶盘 及 分 支组 成 。
究成 果 ,挤 扩 支盘 灌 注桩 中的 支盘相 当 于阻尼 器 ,可以起 到减 震作 用 。考虑相 互作 用 的体 系,结
地震作用下桩-土-结构相互作用的分析

2 1 程 概 况 .工
0 . 4
J
九
.
兰 ”
…
.
M 《
』 / I 1
- s
。
一 0 2
t t s
某 六 层 框 架 结 构 , 高 4米 , 于上 海 南 汇 区 , 区 位 于 长 江 人 海 层 位 该 口, 冲积 平 原 , 质 条件 较 差 。 构 平 面 图如 图 2所 示 , 主要 承 重 构 属 地 结 其 件 的截 面尺 寸 及 材 料 力 学性 能 参 数 如 下 : 框 架 柱 :0 m 4 0 4 0 mx 0 mm混 凝 土 柱 外 环 梁 :0 m 4 0 3 0 mX 0 mm混 凝 土 梁 楼 面粱 :0 5 0丁 字 钢 , H 5 0 m, B 2 0 1x0 高 = 0 m 宽 = 0 mm, 缘 1 1m 翼 1 6 m, = 腹板 t 1rl 2 0ll = l l 楼 面 板 :O m 混 凝 土楼 面板 lO m 材料特性 : 土 C 0 弹 模 E 30 1 / 泊 松 比 V 02 密度 P 昆凝 3, = .x 0 Nm , = ., =
Wike 地基上的梁模型 、 nlr 波动场模型 以及有限元模 型。 3种模型 中由于 动 力 Wik r 基 上 的梁 模 型 力 学 概 念 清 楚 , 单 实 用 , 于被 T 程 人 nl 地 e 简 易 员接 受 而 被 广 泛 运 用 。 本文 拟 采 用 这 一 模 型 , 桩 离 散 为 梁 单 元 , 节 将 在 点处 考虑 水 平 位 移 和 转 角 , 桩 体 的 质 量 集 中在 节 点 处 。土 体 由离 散 的 将 弹簧 一阻尼器系统代替 , 简化模型见图 1 。
9度
有限元模拟土-箱基-框架结构大比例模型动力相互作用研究

中图分类号 : U 2 T 58
文献标识码 : A
文章编号 :62 09 20 )1 05 5 17 —7 2 (06 0 —07 —0
F M t d n d n mi it rcin o ol— o o n ain—f me sr cue E su y o y a c ne a t fs i b xf u d t o — o —r a t t r u
h ya c h at sc r l s di t re i ai t F T a t sp s t e T E f i Tednmi ca c r ts eea ot t ef cdv rtnts ( v ) th u rrc r. eF M (i t r e i w s ee nh o i b o e e e t u h u ne
维普资讯
第 3卷 第 1 期 20 06年 2月
铁道科学与工程学报
J OURNALOF P I AY SCINCE AND GI ERI AlW E EN NE NG
V I3 No 1 o . .
Fb e
有 限元模 拟 土. 基 . 架 结构 箱 框 大 比例 模 型 动 力相 互 作 用研 究
M r建模计算与试验对 比分析 。研 究结果表明 : 土 一结构相 互作 用影 响后 , a c 考虑 结构的 自 频率减小 , 延长 ; 堆载 振 周期 通过
改变结构楼层质量分布后 , 基频减小 , 结构 的动 力特性 变化规律 满足动 力相似 关 系。有 限元 计算 与现 场试验 吻合较好 , 为 研 究如何利 用土 一结构共 同作 用指 导结构设计提供 了参 考。
土动力学(第1章)

从总应力法,到有效应力法; 从单独土的动力特性研究,到考虑土和上部结构的相
互作用; 从简单应力状态,到复杂应力状态; 从室内试验,到现场试验和模型试验。
2021/4/9
14
本章内容结束, 谢谢大家!
2021/4/9
15
感谢您的阅读收藏,谢谢!
2021/4/9
16
2021/4/9
12
土动力学的研究途径
土动力学问题的研究,必须建立在土力学、 地震工程学、结构动力学、土工抗震学等一 系列学科的基础上,并且充分运用现有室内 外试验量测技术,广泛积累原型观测和工程 实录资料,作出理论联系实际的分析,建立 科学的理论方法。目前解决这一问题的途径 有:
一是,建立起一定理论下的计算方法,引入一些表征动力作用和土动力 特性的相应指标,作为动力设计的基础,并在不断地实践中加以检验和 修正;
土动力学
Soil dynamics
骆亚生
二零零八年四月
2021/4/9
1
第一章 绪论
土动力学的定义 动荷载的类型及其对土体的作用特点 动荷载作用对工程建筑的影响 土动力学的任务、内容与研究途径 土动力学的发展历史与趋势
2021/4/9
2
什么是土动力学?
土动力学是土力学的一个分支,是研究各种 (由地震、爆炸、海浪以及动力机械等所引起的) 动荷载作用下土的变形强度特性及土体变形稳定 性的一门科学。
①工程建筑中的各种动荷作用及其特点问题; ②动荷所引起的振动和波动及其在土中运动规律问题; ③土的动应力~动应变关系问题; ④土的动强度和动变形问题; ⑤土的振动液化问题; ⑥土动力特性测试设备与测试技术问题; ⑦动荷条件下的地基承载力,土坡稳定及挡土墙土压力问题; ⑧土与结构物的相互作用问题(包括动力机器基础问题)。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
7 & 前言
地震时土体与结构的相互作用是一个普遍存在的问题9 土 > 结构物的动力相互作用问题, 是一个涉及到 土动力学、 结构动力学、 非线性振动理论、 地震工程学、 岩土及结构抗震工程学、 计算力学及计算机技术等众 多学科的交叉性研究课题, 也是一个涉及到非线性、 大变形、 接触面、 局部不连续等当代力学领域众多理论与 技术热点的前沿性研究课题。随着科学计算技术的迅猛发展和实验手段的不断改进, 重大和复杂体系工程 的不断建造, 促进了土与结构动力相互作用的深入研究, 几十年来一直引起国内外的广泛重视和研究。"@$B 年日本新泻地震、 "@?$ 年我国唐山地震、 "@]= 年墨西哥地震和等许多实践课题促进了这门学科的迅速发展,
!" #$%&$’%( #" &%)%*&+, #- ./"*0’+ )#’12)3&4+34&% ’"3%&*+3’#"
D-34 EFG" *H /;7A8G;IJ" ’-34 *7IJ!
( "9 +KL;<MNKIM 7O 4K7MKPFIGP;8 KIJGIKK<GIJ,.7IJQG 1IGRK<SGMT, UF;IJF;G, !###@! , VFGI;; !9 ’HT;IJ V7ISM<HPMG7I :KPF;IGSN V79 *MW, 4H;IJXF7H, ="#B=# , VFGI;)
由于现代数值计算理论和计算机技术的发展, 以及一些重大工程的相继修建, 从理论和实践中都将极大 地推动了土 , 结构动力相互作用问题研究的迅速发展。
! Y 土 , 结构动力相互作用的研究方法
Y Y 土体 , 结构动力相互作用的研究方法可总结为以下几种: ( 包括现场震害调查和足尺试 ! 原型观测 验) ; ( 主要是动力离心模型试验, 实验室内小比例模型试验) ; "模型试验 # 数理模型与计算方法。由于地 万 方数据 震载荷具有随机性、 难以预测性以及土体与结构物及其系统本身的复杂性, 目前还没有哪一种手段能够独立
[ "] "@@= 年日本神户大地震、 "@@@ 年土耳其地震和中国台湾地震 等使土动力学和土与结构动力相互作用的研
究达到了一个新的高潮, 取得了丰硕的成果。
8 & 研究历史与现状
[ !] & & 土 > 结构动力相互作用的研究最早可以追溯到 "@#B 年 *;NY 对弹性地基振动问题的分析。"@%$ 年
系数即动力阻抗来表示力和位移的关系。 研究者采用不同的吸收边界, 利用有限元计算各种埋置基础的阻抗函数。(>$*"& 和 D))$ 在 -.@K 年提 出了一种能够传播洛夫波和瑞利波的吸收边界, 这一边界应用于坐落在坚硬基岩及水平分层介质上的基础
[ -/ ] [ -9 ] 的频域动力分析 。J)1$"2 将这一方法扩展到轴对称问题并得出不同基础形式的阻抗函数 。-.@9 年 [ -@ ] (16; 利用积分的方法得出粘弹性分层介质上刚性基础的阻抗函数 。 M)$51H8)、 J)*"$I)&) 和 !); 把 (>$N [ -K ] [ -A ] *"& 的研究扩展到研究平面、 轴对称及三维基础问题, 同时考虑了材料的非线性及非均质性 。
=)%"&V"2 利用格林函数形成成层土的刚度矩阵, 利用有限差分能量法形成板的刚度矩阵, 给出了圆形基础的
[ K. ] 动力反应 。D;2? 和 4;%5 于 -..3 年提出了微分有限元核方法, 并用于求解时域内无限介质的动态刚度矩 [ ’F ] [ ’- ] 阵 。-..@ 年 B)% 提出非线性土中桩的动力反应 。-..A 年, ="&%)2 和 X;1$"?? 提出了非线性土中结构共 [ ’K ] 同作用的时域频域混合法 。职洪涛、 俞载道、 曹国敖应用参数分布法, 对基础产生提离和滑移的土 , 结构 [ ’’ ] 。 动力相互作用实用的方法得到了发展, 并逐渐向三维分析过渡, 而非线性分析逐渐 成为研究的主流方向。J)&)+)2#$ 和 ="$P;$ 利用 =ST ( 边界元) 在时域中分析了三维刚性地面上明置基础及
[ KK , K’ ] 埋置基础的阻抗函数 。-.A/ 年 :+)$6)2 和 M;*#%51"O 利用 UST ( 有限元), =ST 求得了柔性基础的动 [ K3 ] [ K/ ] 力反应 , 并在 -.A9 年求得基础置于分层非均质土体上的阻抗函数 。 :7*)< 和 =)%"&V"2 ( 2.AA ) 利用频 [ K9 ] 域中的二次单元求得基础的柔度函数, 并利用边界元求得柔性埋置基础的动力反应 。 J)&)+)2#$ 和 ="$N [ K@ ] [ KA ] 以及 W)#8)%)&;$ 和 J)&)+)2#$ 在时域中利用 =ST , UST 分析得出了埋置基础的动力反应。:7*)< 和 P;$
[ ’] !"#$$%"& 通过对 ()*+ 解的积分, 研究了刚性圆形基础板在竖向荷载作用的简化边界条件下的振动问题
( 即通常所称的基础振动问题的 !"#$$%"& 解) , 标志着土 , 结构物动力相互作用研究的开始。此后这一领域
[ 3] 得到了众多研究者的关注。-./’ 年, 01#%2)% 求得刚性基础下的接触压力问题 。 41%5 对 !"#$$%"& 的工作 [ /] 做了有益的补充, 他假定半空间表面上圆形刚性基础接触面内有三种反力分布 ( 均匀分布、 抛物线分布及
刚性基础的压力分布) , 给出了不同泊松比下的解, 并得出了半空间基础上条形基础反力的数学表达式。这 一理论可以应用于计算条形基础的柔度函数。!#67)&8 ( -.9/ ) 通过计算对比后认为这一方法具有相当的精
[ 9] 度 。:&%;2< 等 ( -.// ) 及 =>6&;?8 ( -./9 ) 得出 圆 形 基 础 在 摇 摆 振 动 下 的 解, 第一次得出基础下的应力 A] 分布与频率有关的结论[ @ , 。-.99 年 (>$*"& 和 !#67)&8 提出的解决土 , 结构动力相互作用的集总参数
土与结构相互作用问题涉及到半空间无限域的问题, 这就导致边界元法的出现。因为其能够包含模型 中的无限域问题, 所以边界元法大量应用于土 , 结构相互作用问题中。-.@A 年 M;*#%51"O 第一次利用频域
[ -. ] 中边界元法得出明置或埋置于粘弹性半空间中的矩形基础的阻抗函数 。-.AF 年, E"P"&#$ 和王贻荪等人 [ KF ] 对基础振动半空间理论作了进一步分析, 求得了动 =;1$$#%"$Q 问题的精确解 。-.A- 年, R5167# 和 (16; 得 [ K- ] 出了弹性半空间上矩形基础的动力反应 。
间, 上部结构理想化为带刚性底板的单自由度刚架, 其刚性底板搁置在地基土表面。这一力学模型的提出, 标志着土 , 结构动力相互作用研究工作的深化。 G7;H&), E"&1*)2$I)*# 在分析大坝与基础在地震时的相互
[ -- ] 作用时提出了子结构法 , 使数值计算能够在动力体系中得以有效应用。 (>$*"& 和 J172"*>"& 提出的吸收 [ -K ] [ -’ ] 边界包含黏性阻尼和瑞利波阻尼, 并且得出了条形基础及对称基础的柔度函数 。 (16; 和 D"$8*)% 以 [ -3 ] 及 L"2"8$;$,D"# 得到了复合边界条件的解, 前者是用柔度来表示力和位移的关系。后者则用刚度和阻尼
"期
方< 志等: 土 ! 结构动力相互作用研究综述
+#
地对土体 ! 结构物动力相互作用问题进行全面而真实的解释和模拟。只有将这些不同的手段有机结合起 来, 才可能真正揭示动力相互作用问题的机理, 并对其行为和过程进行较为准确的描述。具体来说, 通过模 型试验和原型观测结果部分的或定性的再现实际现象、 解释物理机制、 推断变化过程、 总结特性规律和分析 灾变后果, 在此基础上建立合理的能够反映实际动力相互作用规律的数理分析模型, 发展相应的解析或数值 分析方法, 再通过模型试验和原型观测结果加以验证。这是研究和解决动力相互作用这一复杂问题的较为 合理的途径。例如从 "#$$ 年开始, 有关研究人员对 %&’&() 图书馆进行激振试验与观测, 对建筑结构在动力 作用下土 ! 结构相互作用有了很大程度的认识。"#*+ 年在台湾罗东建成两座核电厂安全壳模型, 研究土 ! 结构动力相互作用对核电厂结构地震反应的影响, 取得了大量的观测数据及试验数据。"#*, 年日本动力工 程试验中心在福岛核电厂进行了一系列的大比例尺模型试验。"#** 年小掘铎二等对一座五层的钢框架房
( "9 同济大学 地下建筑与工程系,上海 !###@! ; !9 广州五羊建设机械有限公司, 广东 广州 ="#B=# )
摘要: 土与结构动力相互作用是当代力学领域的前沿性研究课题, 具有很强的实践性。对土与结构动 力相互作用的研究历史与现状进行了介绍, 简要综述了当前土与结构动力相互作用的研究方法, 重点 介绍了目前关于土与结构动力相互作用问题中从无限域转化成有限域的人工边界研究进展问题, 并 对该领域今后的研究工作提出了建议。 关键词: 土 > 结构动力相互作用; 土动力学; 地震工程; 现状; 人工边界; 综述 中图分类号: C%"= & & & 文献标识码: -