地震作用下山区高墩桥梁随机地震响应分析.docxVIP

地震作用下山区高墩桥梁随机地震响应分析 0 主梁和斜拉桥的地震动响应 山区的高速铁路桥梁具有墩高大、墩高不规则、墩高差异大等特点。同时，山区地形的变化复杂，严重影响了山区高墩桥的抗疲劳动性能。因此，考虑到振动的空间效应（位置效应、相干效应和行波效应）对高提桥的抗强震力作用非常重要。在考虑地震动空间性的基础上,进一步考虑墩高变化对墩内力及桥梁上部结构地震响应影响是必要的。文献引入过滤白噪声随机地震动模型分别对大跨度拱桥和斜拉桥进行反应谱和功率谱分析,考虑场地效应、相干效应及行波效应对地震响应的影响;文献研究了地震作用下多跨简支梁桥碰撞参数与空间变化的地震动关系,建立三维有限元模型分析了梁体间初始间隙随地震动空间变化对碰撞的放大作用;文献的研究表明忽略地震动空间变化效应会低估或高估结构的地震动响应,结构地震响应不仅与结构本身动力特性有关而且与地震动空间变化关系密切;文献对南京大胜关桥进行了非一致地震响应分析,考虑了不同地震波速对结构响应的影响;文献对基于随机振动理论,详细研究了场地效应对高墩桥梁的影响,同时,也研究了平稳虚拟激励法在有限元软件中的实现,提高了计算效率; 文献考虑了局部场地效应对大跨度斜拉桥的地震动随机响应的影响,与一致激励相比,局部场地效应对结构的内力和位移影响显著,场地效应对主梁、索塔内力和位移的影响程度及规律也不尽相同,需区别对待分析。 基于以上研究,国内外对山区高墩桥梁研究不多,用抗震理论发展必然趋势的随机振动理论研究高墩桥梁抗震性能就更少,故本文基于随机振动理论,建立了三维有限元模型,对山区高墩桥梁进行多点激励随机振动分析,研究了强地震作用下高墩桥梁响应与地震动空间性和墩高之间的关系。 1 相对动力响应方程 地震作用下,结构系统的运动方程可分解为关于非支座节点和支座节点的分块矩阵形式 [ΜffΜfrΜrfΜrr]{¨uf¨ur}+[CffCfrCrfCrr]{˙uf˙ur}+[ΚffΚfrΚrfΚrr]{ufur}={0F}?(1)[MffMrfMfrMrr]{u¨fu¨r}+[CffCrfCfrCrr]{u˙fu˙r}+[KffKrfKfrKrr]{ufur}={0F}?(1) 式中,下标f表示非支座节点,r表示支座节点;M、C、K分别为质量矩阵、阻尼矩阵和刚度矩阵;uf为非支座节点自由度;ur为支座节点自由度;F为支座激励力。将非支座节点位移分解为拟静力位移和动力位移,即: uf=us+ud。(2)uf=us+ud。(2) 将(2)代入(1)并忽略惯性力和阻尼力项,由静力计算得到拟静力位移响应: us=-Κ-1ffΚfrur=Aur?(3)us=?K?1ffKfrur=Aur?(3) 式中A表示支座单位位移所引起的结构的拟静力位移。将(2)和(3)代入(1)得相对动力位移响应方程为: Μff?¨ud+Cff˙ud+Κffud?-(ΜffA+Μfr)¨ur。(4) 运用模态叠加原理对上述方程解耦,可设: ud(t)=?y(t)?(5) 将(5)代入(4)被解耦后方程为: ¨yj+2ξjωj˙yj+ω2jyj=Gj?(j=1,2,3,?,n)?(6) 式中,n为模态数;yj为广义坐标;ωj和ξj分别为第j阶频率和模态阻尼比;Gj为模态荷载。 Gj=ΓΤj¨ur。(7) 对应支座激励的模态参与系数为: Γj=-(ΜffA+Μfr)Τφj。(8) 相对动力和拟静力功率谱密度函数分别为: Sdi(ω)=n∑j=1n∑k=1φijφik[r1∑?=1r2∑m=1γ?jγmkΗ*j(ω)Ηk(ω)ˉS?m(ω)]+n∑j=1n∑k=1φijφik[r2∑?=1r2∑m=1Γ?jΓmkΗ*j(ω)Ηk(ω)?S?m(ω)],(9)Ssi(ω)=r1∑?=1r2∑m=1Ai?Aim[1ω4?S?m(ω)]。(10) 拟静力项和相对动力项耦合部分为: Ssdi=n∑j=1r1∑?=1r2∑m=1?ijAie[-1ω2ΓmjΗj(ω)?S?m(ω)]?(11) 式中,n为所选模态数;r1和r2分别为非支座节点数和支座节点数。结构频率响应函数为: Ηj(ω)=1ω2j-ω2+i(2ξjωjω)?(12) 式中,ω为荷载频率;ωj为结构第j阶频率。 结构反应均方值表示为: σ2fi=∫∞0Sdi(ω)dω+∫∞0Ssi(ω)dω+ 2|∫∞0Ssdi(ω)dω|Re=σ2di+σ2si+2Cv(usi,udi), (13) 式中,|·|表示取实部;σ2di表示第i阶相对动力位移响应方差;σ2si表示第i阶拟静力位移响应方差;Cv(usi,udi)表示拟静力和相对动力项协方差。多点激励下,功率谱矩阵为式(14): [S(ω)]=[S11(ω)C12(ω)+iQ12(ω)C13(ω)+iQ13(ω)?C1k(ω)+iQ1