大跨斜拉桥多点激励随机地震响应研究.docVIP

大跨斜拉桥多点激励随机地震响应研究 　　摘要： 虚拟激励法可以高效地计算结构随机响应，但目前并没有整合到既有有限元软件中，阻碍了其推广应用，在很多情况下尚需针对具体模型进行编程计算。将该方法整合到既有的有限元软件中，对虚拟激励法的推广使用具有重要意义。大跨柔性结构的初始平衡对其动力响应具有重要影响。在既有研究的基础上，基于Abaqus/Python平台实现了大跨柔性结构需考虑初始平衡态影响的多点激励随机响应计算流程。以目前跨度最大的双层公路斜拉桥——上海闵浦大桥为例，研究大跨斜拉桥多点激励随机地震响应的特点，并分析了视波速对随机地震响应的影响。关键词： 大跨斜拉桥； 随机地震响应； 虚拟激励法； 多点激励 　　中图分类号： U441+.3文献标识码： A文章编号：2013 　　引言 　　文献研究表明地震多点激励对大跨度桥梁具有不可忽视的影响，是大跨桥梁抗震研究的重要内容之一，其分析方法包括确定性动力分析法（时域、频域、多点反应谱法）和随机振动分析法。尽管确定性分析方法具有诸多优点，但其分析结果对地震输入具有一定的依赖性，因而采用随机振动分析法研究大桥梁具有重要意义。Allam（1999）根据简化平面斜拉桥主梁连续体振动方程，将弹性拉索简化为竖向弹性支承，研究了多点随机激励各种参数对地震响应的影响[1]。Dumanogluid（2003）以Jindo斜拉桥为例，建立2D有限元模型，综合研究上述3种效应对大跨斜拉桥的随机地震影响[2]。DebChaudhury（1982，1988）给出了Mdof系统基于模态叠加的多点地震激励响应均方根求解方法，该方法利用绝对坐标下的结构动力方程构造状态方程，并将输入功率谱模型的滤波方程整合到上述状态方程中，将其转换为一阶常微分方程组，该方法具有一定的普遍性[3，4]；但其计算公式甚为复杂，针对具体问题尚需大量编程工作。相对于将响应分解为动力部分（Dynamic）和拟静力部分（Static）的处理方式[5]，上述方法具有一定的新颖性，并且可应用于非平稳随机响应分析中；对于简单结构具有一定的优势。林家浩提出的随机振动虚拟激励法（Pseudoexcitation Method，abbr. PEM）是一种有效的分析工具[6]。 　　虚拟激励法通过构造虚拟激励，将原来的随机振动分析问题，转化为计算结构在一系列虚拟激励作用下的确定性动力学计算问题，极大地降低了计算量，系一种全新的随机振动求解方法体系[7]。虚拟激励法在随机响应计算方面具有高效性，但目前并未整合到既有商业有限元软件中，一般需编程计算，给虚拟激励法的推广应用带来了很大的障碍，因而将其整合到既有商业有限元软件中具有重要意义。江洋（2010）基于Ansys平台上实现了多点激励随机振动计算的模态叠加的方法，其计算主要包括3个方面：（a）构造力荷载虚拟激励（包括激励功率谱复分解、耦合刚度矩阵和质量矩阵提取）；（b）进行谐分析（Harmonic）分析获取虚拟响应；（c）编程计算结构响应的功率谱、方差等物理量[8]。李永华亦给出了类似的流程及算例[9]。文献[8，9]所提出的分析流程对将虚拟激励法用于大型复杂模型的随机激励响应分析具有重要意义，但其采用Ansys平台时忽略了结构初始平衡态的影响。对于大跨柔性结构，由于几何非线性的影响，其动力响应分析（地震、风等荷载作用下的时域或随机响应分析）一般需以结构的静力初始平衡态为基础上进行相应分析，即以结构静力作用下变形后的构型为初始状态，对其切线刚度矩阵与质量矩阵进行后继的分析。文献[8，9]所给研究算例中均未体现重力的作用，即没有考虑初始平衡态对结构随机响应的影响。在荷载（如恒载、预应力等）对结构频率影响较小的情况下，上述流程是适用的。 　　本文在既有研究的基础上，基于Abaqus/Python平台实现需考虑初始平衡态影响的大跨柔性结构多点随机激励响应分析方法，给出相应分析流程；以目前跨度第一的双层公路斜拉桥——上海闵浦大桥为工程背景，分析大跨斜拉桥多点随机激励地震响应的特点。 　　4结束语 　　本文基于Abaqus/Python平台实现了大跨柔性结构需考虑初始平衡态影响的多点激励随机响应计算流程。利用上述方法，以目前跨度最大的双层公路斜拉桥——上海闵浦大桥为例，研究了大跨斜拉桥在考虑行波效应、多点激励下的随机地震响应。研究表明： 　　（1）视波速对主梁位移随机响应RMS的影响不一，除主梁扭转外，视波速越小、主梁位移随机响应RMS越大，因而选择合适的视波速尤为重要。不同视波速下主梁扭转RMS沿主梁分布趋势有所不同。主梁横向位移、竖向位移以及竖向转角RMS均随视波速的增加而趋向于一致随机地震响应并且由于边跨刚度较大，视波速对主梁边跨横向位移、竖向位移以及扭转位移