基于ansys的铁路悬索桥地震响应分析

基于ansys的铁路悬索桥地震响应分析

大桥的抗疲劳性能得到了高度重视，并进行了大量研究。对其分析过程中关键为确定地震动输入。地震动本质为空间变化,由于桥梁跨度大、基底延伸长,地基各点振幅、相位差别较大,考虑地震动多点非一致输入极为重要。Abdel-Ghaffar等利用实际地震记录分别在时、频域内对金门大桥进行地震反应分析。结果表明,一致输入反应分析并非代表最不利情况,行波效应对悬索桥的地震反应影响显著。Nakamura等用多点激励反应谱方法对金门大桥进行分析,认为受其柔性影响,以动力反应分量为主,行波效应及多点激励对动力反应分量影响显著,且多点激励会致桥塔有小响应。Dumanoglu等对三座欧洲大跨度箱梁悬索桥行波效应比较分析发现,在悬索桥的某些部位行波效应会增大地震反应。Harichandran等认为对悬索桥进行一致输入分析会高估某些部位反应,低估其它部位反应。胡世德等通过对江阴长江大桥地震反应分析,认为行波效应会使结构反应增大,而多点激振影响较小。郭发忠通过对某三跨悬索桥的地震响应分析,认为由于多支承激励引起的柔性结构地震反应中拟静力部分占比较小,总体上行波效应对结构内力反应影响不大。有研究只考虑行波效应或直接生成人工地震波,虽用两种方案进行计算,但均未能揭示两种非一致激励考虑方式对桥梁响应影响的根本区别。只考虑行波效应会高估各点间相关性,直接生成人工地震波而不进行谱拟合则会致生成的人工地震波反应谱值与实际地震波不符,如长周期段谱值偏小等,且各点间地震波反应谱值差异亦会造成桥梁地震响应变化。本文以某三跨铁路悬索桥为例,基于大型通用有限元软件ANSYS采用行波输入与多点激励两种非一致地震激励对铁路悬索桥进行地震响应计算分析,获得该桥的响应变化规律。对归一化后桥梁地震响应时程进行傅里叶变换并分析谱值,获得桥梁各响应控制振型变化,给出非一致激励对桥梁地震响应影响机理,揭示两种非一致激励考虑方式对桥梁响应产生影响的根本区别。1地震波加速度分析各点地震动加速度可表示为式中:F(t)为强度包络函数;Anm(ωk),θnm(ωk)为考虑第n点与第m点相关的第k个频率分量幅值及相位角;φmk为随机相位角,在(0,2π)区间均匀分布,且相互独立。频率为ωk的地震动分量功率谱矩阵可表示为自由场两点地震动互功率谱可表示为式中:ρmn(dmn,ωk)为两点间功率谱相干函数;dmn为两点间距离;反映行波作用对多点激励影响。对相干函数,本文采用适用性较强的屈铁军-王君杰-王前信(QWW)模型。通过Cholesky分解为S(iωk)=L(iωk)L*T(iωk),由此,Anm(ωk),θnm(ωk)可据求得。将各点随机相位角项φmk设为一致,即与各点位置无关φk,便得局部场地收敛的多点地震波。进行反应谱拟合时,本文采用考虑相位谱的人造地震动反应谱拟合方法,可得合理的加速度时程曲线,但会在模拟所得地震地面运动中引入长周期分量,造成直接积分的速度、位移时程曲线严重失真。本文通过时域内基线调整及频域内过滤掉不合理长周期分量方法达到零线调整目的。本文加速度时程与某组地震波位移时程对比见图1。由图1看出,考虑多点相关性多点激励输入的各支承点位移时程较接近,无明显延后现象。2中央扣节点本文以拟建某特大铁路悬索桥为计算分析模型,用ANSYS有限元分析软件建模。桥梁为地锚式悬索桥,跨度98m+660m+98m,中跨设置吊杆,在跨中设置三对柔性中央扣,边跨无吊杆,主塔为混凝土,高196m、156m,主梁用上承式钢桁梁,桥面宽24m,桁高12m。主桁与桥塔横梁相交处只设竖向、横桥向支座,顺桥向无约束,主梁两端只约束竖向、横桥向位移。悬索桥立面图及各支承点编号见图2。为简化计算、明确主梁内力变化,桥梁有限元模型用鱼骨式建模,见图3。桥梁前10阶振型及描述见表1。桥梁前10阶振型见图4。设置柔性中央扣后桥梁纵飘振型出现被延后,桥梁1阶振型为主梁对称横弯。桥梁周期较长,表明该桥属柔性结构。应对桥梁内力、位移进行分析研究。3悬索桥承载体地震波的计算本文将系统对大跨度悬索桥在两种不同非一致激励方式下地震响应进行分析。为反映不同输入对桥梁地震响应统计意义的变化规律,生成5组多点地震波进行计算,每组地震波含3种波速,分别为500m/s、1000m/s及2000m/s,每种波速下均含悬索桥的6个支承点地震波。计算一致激励时用多点输入的第1点地震波计算,行波输入则对第1点地震波做延时处理,以保证几种输入方式的可比性。采用常遇地震顺桥向加载,按50年超越概率63%计算,加速度峰值为0.79m/s2,计算中考虑几何非线性影响,并设地震波从左向右传播。计算结果均为用5组地震波分别计算平均值,且均为与自重荷载组合结果。分析主梁位移时为相对主梁左端相对位移,分析桥塔位移时则为相对各塔底相对位移。3.1吊杆应力分析5组地震波作用、不同输入方式下主缆及吊杆应力平均值见图5。由表2看出,随波速变化不同非一致激励方式下地震响应均有不同。行波输入下随波速减小主缆及吊杆最大应力均有增大,但幅度较小,主缆应力最大增大4.82%,吊杆最大应力最大增加3.12%;多点激励作用下主缆及吊杆应力与一致输入基本一致,幅值变化约1%。考虑本文加载地震动峰值较小,地震造成的响应与自重下内力相比较小,按罕遇地震加载时各响应浮动百分比应增大。多种工况作用下吊杆最小应力为0值时出现于跨中正中部位,表明地震作用下该位置吊杆可能出现松弛,设计时应采取必要构造措施。不同输入方式的吊杆应力随吊杆位置分布见图5。由图5看出,两种非一致输入方式的吊杆应力分布与一致激励基本一致,仅应力幅值有区别。3.2梁端顺桥与顺桥向位移不同输入下的主梁最大地震响应见表3。由表3看出,行波输入下,随波速的减小主梁最大竖向位移、弯矩及剪力均呈增长趋势,波速500m/s时增幅分别为403.03%、10.88%及4.37%;而梁端顺桥向位移及主梁最大轴力变化规律不明显,除波速1000m/s有所减小外基本与一致激励持平,且梁端顺桥向位移有减小趋势;多点激励作用下仅主梁竖向位移有所增长,波速1000m/s时增长100%,但总体随波速变化规律不明显;而梁端位移、主梁最大轴力、最大竖向弯矩及剪力亦与一致激励基本持平。为进一步研究两种不同非一致输入下主梁的响应变化规律,主梁竖向位移、主梁轴力、主梁竖向弯矩及主梁竖向剪力沿梁长包络见图6~图9。由4图看出,不同输入下主梁响应发生变化区域均为中跨,两侧边跨响应基本无影响。行波输入下主梁跨中竖向位移显著增大,包络形状由一致激励的中跨1/4跨、3/4跨双峰变为1/2跨单峰,跨中弯矩明显增大,剪力变化不明显。多点激励输入下主梁跨中竖向位移有一定程度增长,但包络仍为双峰形式,主梁竖向弯矩、剪力包络无明显变化。两种输入下主梁轴力分布均与一致激励相同,区别仅在于峰值大小。因此,两种输入方式对主梁响应有较大不同,考虑非一致激励影响时应谨慎选则。3.3基本与一致激励并存不同输入方式下两侧桥塔最大顺桥向地震响应见表4。由表4看出,行波输入下桥塔最大顺桥向位移与塔底轴力均有所增加,其中矮塔顺桥向位移增加107.28%,塔底轴力增加2.87%,而塔底弯矩及剪力幅值则略有减小,最大减幅分别达27.21%及26.48%,但随波速增大各项响应均越接近一致激励结果;多点激励作用下各项地震响应均基本与一致激励持平,无明显变化规律。图11~图14为桥塔顺桥向位移、轴力、顺桥向弯矩及顺桥向剪力包络,由于包络线基本对称,仅给出一侧。由4图看出,不同输入下,桥塔顺桥向位移包络发生明显变化:一致输入时桥塔最大顺桥向位移发生于塔梁固结处与塔顶间桥塔中部,考虑非一致激励后塔顶顺桥向位移均不同程度增大,矮塔更明显。桥塔轴力包络无明显区别,仅幅值有所区别。桥塔顺桥向弯矩(剪力)包络在波速500m/s时变化较大,并仍以矮塔明显。此矮塔支承点输入地震波延时较长,且矮塔自振周期较高塔短的缘故。对桥塔各项地震响应傅里叶幅值谱进行分析。矮塔塔底顺桥向弯矩幅值谱见图15。由图15看出,该结果与主梁分析结果一致,桥塔地震响应幅值谱仍为行波输入下0.253Hz处幅值增加,而多点激励作用下则无明显变化。结果表明,行波输入对大跨度悬索桥影响主要表现为对一阶竖向振型激发,且能使该振型成为部分响应的控制振型,多点激励对主要振型无显著影响。其原因为行波输入下各点间相对位移较大,拟静力位移导致桥梁一阶振型激发,多点激励因各点间位移时程相差较小,故无此现象产生。4行波输入与多点激励作用下桥梁地震响应本文结合某大跨度铁路悬索桥,基于大型有限元通用软件ANSYS,进行两种非一致激励下地震动响应分析,结论如下:(1)行波输入下桥梁主缆及吊杆应力随波速减小而增加,而多点激励作用下,主缆及吊杆应力与一致输入基本相同。两种输入方式下,吊杆应力分布与一致激励无明显区别。(2)两种非一致激励作用时主梁竖向位移较一致激励均有所增大,行波输入时增大更显著,且随波速减小增幅增大,并伴随主梁竖向弯矩、位移包络形状改变;多点输入作用时主梁竖向位移随波速变化规律不明显,竖向弯矩、位移包络无明显改变。两种非一致激励作用时主梁竖向剪力峰值及包络变化不明显。(3)行波输入时桥塔最大顺桥向位移与塔底轴力均有所增大。而塔底弯矩及剪力幅值则略有减小,但波速越快各项响应越接近一致激励结果;多点激励作用下各项地震响应均与一致激励基本持平,无明显变化规律。两种非一致激励作用时塔顶位移均显著增大。(4)行波输入对大跨度悬索桥影响主要表现为对一阶竖向振型激发,而多点激励对主要控制振型无显著影响。各点间相对位移较大,拟静力位移致桥梁一阶振型激发,多点激励由于各点间位移时程相差较小。(5)不同输入下悬索桥地震响应的主要控制振型未发生改变时其响应幅值变化不会太大。进行悬索桥地震响应分析前应先分析各响应主要控制振型,以便选择合理的地震动输入方式。(6)本文计算结果仅考虑常遇地震,所得桥梁地震响应变化幅度仅作规律研究。罕遇