高架支梁非线性碰撞地震反应分析.docVIP

高架简支梁桥非线性碰撞地震反应分析 郭维 铁道部科学研究院 博士 [提要] 本文采用MSC.Nastran非线性单元库中的gap单元（间隙元）来模拟高架简支梁桥在伸缩缝处的碰撞作用，并就非线性碰撞单元的力学特性和求解非线性过程中需注意的问题进行了分析和探讨，最后研究了伸缩缝宽度、单（双）边碰撞和邻梁摩擦对碰撞地震反应的影响，结果表明，NASTRAN的非线性迭代功能及gap元的引入对处理类似的桥梁碰撞问题具有很大的优越性，计算结果令人满意。 1. 研究背景 在我国历史上记录着许多由于碰撞引起的桥梁震害现象。例如唐山大地震中滦河桥的严重落梁破坏就是一定程度的碰撞作用引起的，这种碰撞行为一般造成主梁梁端开裂、桥台胸墙局部混凝土脱落、伸缩缝挤压等结构轻微的破坏，但在某些情况下碰撞却会引起桥梁结构发生严重的落梁破坏。在国外，1971年圣·费尔南多地震中曾观察到桥台与桥面的碰撞现象；在1994年Northridge地震中，在洲际5号道和洲路14号道交叉道的伸缩铰处和附近一些桥梁的耸立桥台处产生了严重的碰撞破坏，这引起了有关专家的高度重视[1]；1989年Loma Prieta地震中，一段低矮的路面与支承高架桥南段桥面的高柱发生的碰撞导致路面和柱的两侧遭受严重损坏。 为从本质上认识桥梁碰撞的发生机制和作用原理，本文将利用MSC.Nastran这一国际著名的有限元分析程序对此加以探讨。 2．非线性碰撞单元的力学特性模拟 2.1 碰撞单元力学特性[2] MSC.Nastran是基于罚的方法来求解点对点的接触问题，这种方法本身具有的参数敏感性决定了在使用gap元时将出现难以收敛或发散的问题，在MSC.Nastran V70.7中，gap元的求解能力较前几个版本有了很大提高，它基于子增量法，具有适时更新刚度、二分增量步长、在动力和非线性环境下自动调整罚刚度值等优点，而且还引入了动摩擦的影响，从而使数值稳定性和解的精度得到了控制。 gap元的单元坐标系如（图1）所示，gap元的轴向力和侧向力—位移曲线如（图2）和（图3）所示。gap元的内力有三个，沿单元坐标系的位移(U，V，W)方向分别是轴向压力（Fx）和其它两个方向的摩擦力（Fy，Fz）。下面就根据gap元所处的状态来计算这些内力，其中： U0——初始间隙；F0——预加载； Kb——张开刚度，缺省Kb=10-14Ka；Kt——横向剪切刚度； Ka——闭合刚度，当UA-UB(U0时起作用； (s——静摩擦系数；(k——动摩擦系数； UA，UB——节点A，B沿gap元轴向的位移； VA，VB，WA，WB——节点A，B沿垂直于gap元轴向方向的位移； 当gap元张开时，两点无接触，无剪切刚度。 Fx=Kb·U(0 且Fy=Fz=0，Kb是一个任意的张开刚度。 当gap元闭合且滑动时，无摩擦，称为“sliding”状态。 Fx=Ka·U0 且Fy=Fz=0，Ka是闭合刚度罚值。 当gap元闭合且粘附着时，有静摩擦，称为“sticking”状态。 Fy2+Fz2(((s·Fx)2 其中Fx=Ka·U，Fy=Kt·V，Fz=Kt·W，注意到“sticking”状态的出现是以静摩擦系数(s和横向剪切刚度Kt为标志的。 当gap元闭合且滑动时，有动摩擦，称为“slipping ”状态。 Fy2+Fz2((k·Fx)2 其中 Fx=Ka·U，，，“slipping”状态的出现是以动摩擦系数(k为标志的。 gap元的单刚形式如下： ， 依据gap元的不同状态，的形式有所不同。 （图1） gap单元坐标系 （图2） gap元轴向力—变形关系 （图3） gap元侧向力—变形关系 （图4） 高架桥动力分析模型 当gap元张开时，，当gap元闭合处于“sliding”状态时， 当gap元闭合处于“sticking”状态时， 当gap元闭合处于“slipping”状态时， 2.2 采用gap单元求解非线性应注意的问题 (一)、gap 元属于非线性单元，它的轴向和侧向力—位移关系为非线性，当其处于不同状态时，刚度发生突变，可能引起发散或产生振荡解，因此必须采用适应性的刚度更新准则来适时地调整刚度矩阵以确保取得收敛解。 (二)、在用gap元求解非线性问题时，结果的非线性力和非线性位移很大程度上依赖于参数的选取，尤其是Ka，Kt，即自身刚度的选择至关重要。Ka的取值过大将降低解的收敛甚至引起发散，而较小的Ka又会导致数值模拟的不精确。因此，gap元的刚度必须为有限，且大小必须小心的加以选择以避免有限元模拟的失真。经验表明，采用Ka=108，Kb=10-4来进行试算是可行的。缺省的横向剪切刚度Kt=(s·Ka，推荐Kt