重力式挡土墙地震响应分析.docxVIP

重力式挡土墙地震响应分析 1 小应变条件下的耦合方法在挡土墙地震模拟中的应用 基于颗粒介质力学理论，从薄视力学的角度研究了砂、土和结构之间的相互作用，成为一个新的焦点。由于计算机性能的限制,用离散颗粒模拟砂土时,颗粒数目过多的情况下不可避免地遇到运算无法执行或计算过于缓慢,因此出现了土体离散-连续耦合方法以减少颗粒数量,只在关心的区域进行细观层次模拟。Cai等在临近隧道开挖的局部土体区域用离散颗粒模拟,其他区域用连续模型模拟,但连续模型是按照静力条件设定的。对于挡土墙和砂土的相互作用,已有基于离散颗粒方法的挡土墙主动和被动土压力模拟。挡土墙的地震响应动力分析中,若土体全部采用离散颗粒模拟,运算时间会过长,因此在分析中尝试用耦合方法进行模拟。耦合方法涉及颗粒材料宏细观参数联系的问题,对二维问题已有的模拟方法有:双轴数值试验模拟静模量和循环双轴试验模拟动力性质。但对于动力耦合分析来说,小应变情况下的颗粒集合体动模量的模拟方法仍不完善。 本文使用离散-连续耦合动力模型,模拟了地震作用下挡墙的动力响应。对靠近结构物的近场土体采用离散元方法进行精细化分析,砂土与地下结构之间满足非连续动力接触关系,振动过程中砂土对结构物的作用直接以接触力的形式作用于结构物上;远离结构的土体采用连续模型模拟。耦合方法的两个核心问题:(1)保证离散-连续土体在边界上的连续性,这里保证边界处速度和力的连续性,本文提出新的边界处耦合力提取方法;(2)模拟中细观区域颗粒的宏观性质需与连续土体模型一致,为完善离散颗粒宏细观参数联系的模拟方法,首次给出颗粒材料的自振柱模拟,获得了小应变下颗粒集合体动模量随应变衰减的关系。对于结构与砂土的相互作用,模拟中采用与离散-连续土体边界耦合一样的方法,但结构-砂土边界处耦合力直接采用颗粒与结构的接触力。 本文使用的耦合方法易于实现,细观土体区域采用颗粒流软件PFC2D模拟,远离结构的土体用有限差分软件FLAC2D模拟,每个循环时步中基于两个软件的数据传输实现宏细观土体的耦合。为避免有限差分算法中结构与土体模量差别过大而导致的时步过小问题,FLAC2D软件只模拟了连续土体区域,对于挡土墙编写了梁单元受碰撞的动力有限元程序嵌入颗粒流软件PFC2D中。 通过模拟文献中重力式挡土墙地震响应的离心机试验,说明了耦合的实现过程。耦合方法从细观层次上模拟了结构受砂土作用的动力过程,并有效减少了颗粒数目,节省运算时间。 2 小应变下颗粒单体的动模量模拟方法 离散-连续耦合方法中,必须保证颗粒集合体的宏观性质与连续体模型中的参数一致,这涉及颗粒集合体宏观性质的模拟。砂这类离散颗粒的动模量随应变的增加而衰减,在室内试验中对于应变小于10-4情况下的动模量一般用共振柱试验获得,而对于应变大于10-4动模量一般用动三轴试验获得。对于颗粒集合体的循环双轴模拟方法文献中已有叙述,但对适用于小应变下颗粒集合体的动模量模拟方法尚未见文献报道。连续土体模型的控制方程中,需代入完整的颗粒集合体动模量与应变的关系曲线,因此建立小应变下颗粒集合体的动模量模拟方法是必须的。 应当指出,用离散元方法直接模拟室内共振柱试验是非常耗时的,原因在于共振柱方法是采用频率扫描的方法,确定圆柱试样的一阶固有频率,再结合试样尺寸和密度,根据弹性波在一维杆中的传播公式推求试样的动模量。对于单一频率的动载荷,用离散元方法模拟试样的振幅响应较为耗时。实际模拟时需对某一带宽的频率进行扫描,还需模拟不同应变下的动模量,由于计算机性能的限制,计算时间过长,难以实现。 1981-1982年间,在共振柱试验的基础上发展了自振柱试验。这种新的试验方法是预先对圆柱试样施加一给定应变,释放应变限制后让试样按其固有特性作自由振动,根据振动的自由衰减曲线获得试样的固有频率,再推求试样的动模量。由于不需要频率扫描,用离散元模拟自振柱试验耗时远小于模拟共振柱试验。 这里给出了自振柱试验的离散元模拟方法,结合动三轴模拟结果,得到完整的动模量随应变变化关系曲线。 2.1 砂级配线的模拟 试验中所用的砂为Nevada-120砂,砂土性质见表1,砂土级配见表2。 在离散元数值模拟中,由于计算机性能的限制,即使采用宏细观耦合方法,仍须采用将实际砂土粒径放大的办法减少颗粒数量。数值模拟试样和离心机试验所采用的Nevada-120砂的级配曲线如图1所示。 模拟时,需使数值试样的宏观性质对应于实际试验中砂的宏观性质。文献虽然是基于干砂的离心机试验,但没有提到干砂试样的三轴试验结果,仅提到饱和试样的三轴试验结果。文献通过对双轴数值试样施加等幅应变,并控制加荷过程中试样面积不变的方法模拟砂土的液化性质。这里采用其加荷过程中试样面积不变的方法,模拟饱和砂土的双轴试验,不断调整颗粒的细观参数,使数值试样在饱和双轴试验中摩擦角与实际饱和砂