摘要：为研究适用于跨断层桥梁的地震动输入分析模型，以1座跨断层桥梁为例，分别计算了非一致激励位移输入模型和大质量法加速度输入模型的地震响应，探讨了2种输入模型的适用性及模型误差产生原因。结果表明：位移输入模型可考虑永久地面位移对结构地震响应的影响，能够真实反映跨断层桥梁在地面运动结束之后桥墩具有的残余内力和变形，适用于跨断层桥梁非一致激励地震响应分析，而大质量法加速度输入模型对此无法考虑，可能导致不合理也是不安全的计算结果；在运用通用有限元计算软件进行加速度输入模型多点激励分析时，应首先核查与拟输入加速度时程对应的位移时程是否存在永久地面位移，若存在永久地面位移，则建议改为位移输入模型进行计算，否则将可能得到失真的计算结果。
关键词：多点激励；跨断层桥梁；位移输入模型；加速度输入模型
中图分类号：U442.55 文献标志码：A
0 引 言
跨越活动断层的桥梁结构地震时具有较大的破坏风险，为此许多国家和地区出台了相关规范和条例以禁止在活动断层之上新建桥梁。然而对于跨越峡谷、河流等障碍物的桥梁结构，往往因各种客观条件限制，不得不采用桥梁形式跨越活动断层。
目前针对跨断层桥梁的研究工作尚处在初始阶段，仅有少量学者对此进行了研究。Park等[12]对土耳其杜杰地震中发生严重破坏的跨断层桥梁博卢1号高架桥的抗震性能进行了分析，认为跨断层效应是影响减隔震桥梁地震响应的重要因素。Goel等[34]研究了跨断层常规中小跨直线桥和曲线桥地震响应需求的简化计算方法。Saiidi等[5]以两跨刚构桥为研究对象，首次开展了跨断层桥梁的振动台试验，推进了该领域的研究工作。杨怀宇等[6]以跨断层简支桥梁为例，采用多点激励时程分析方法研究了跨断层效应对结构地震响应的影响。惠迎新等[79]从震害、地震动输入、分析模型及地震响应等方面对跨断层桥梁的抗震问题进行了较为系统的研究。
跨断层桥梁由于发震断层地表相对错动导致断层两侧支撑具有不同的甚至完全相反的地面运动特征，应采用多点激励输入模型计算其在地震作用下的结构地震响应。在进行结构多点激励分析时，基于位移输入法的位移输入模型和大质量法加速度输入模型在各类大跨桥梁、房屋建筑结构中广泛应用[1011]，以往很多文献对这2种模型在多点激励输入中的适用性进行了探讨、论证[1213]，认为2种输入模型的运动方程均建立在绝对坐标系下，计算结果中包含了结构的拟静力反应和动力反应，既适用于线性结构也适用于非线性结构的地震响应分析。在既有考虑行波效应、部分相干效应以及局部场地效应的多点激励研究中，采用的输入地震动时程均为未包含永久地面位移的简单地震波，而对于以长周期、永久地面位移为主要特征的跨断层桥梁地震动激励，以往研究中采用的位移输入模型和大质量法加速度模型是否同样适用有待于进一步研究。
为此，本文以1座跨断层桥梁为例，结合跨断层桥梁的地面运动特点，以位移输入模型和大质量法加速度输入模型作为地震动输入，对比分析不同模型和不同地震动激励下的结构响应，探讨适用于跨断层桥梁的地震动输入模型。
3 模型建立及地震动输入
3.1 模型的建立
以1座跨越走滑断层的独塔斜拉桥为研究对象，活动断层穿过该桥的第1跨，断层走向与桥位基本垂直，桥梁布置及断层走向如图1所示，其中Fp为平行于断层方向的地震动分量。采用SAP2000
结构分析程序建立该桥计算模型；由于本文仅对地震动输入模型的适用性进行定性研究，计算单元均采用线弹性单元，其中主梁和桥墩采用三维线弹性单元模拟，斜拉索采用桁架单元模拟，桩土相互作用采用承台底加6个自由度的弹簧模拟。时程分析采用Newmarkβ直接积分法，其中参数α=0.5，参数β=0.25。为避免采用Rayleigh阻尼对2种输入模型计算结果可能产生的误差，时程分析中采用常数阻尼，阻尼比为5%。
3.2 地震动输入
地震时走滑断层在平行于断层方向表现为滑冲效应，即由于断层滑动突然升高或降低形成台阶，也称之为永久地面位移。滑冲效应使得断层两侧支承在平行于断层方向呈典型多点激励运动特征[17]。因此，本文分析中以该方向为例，进行多点激励地震响应分析。对于断层中的某点而言，可认为在平行于断层方向断层两侧桥梁支承地震动输入方向相反，数值相等[18]。图2为断层两侧桥墩地震动输入方向示意。
为对比分析多点激励位移输入模型和大质量法加速度输入模型在跨断层桥梁中的适用性，计算中以未包含永久地面位移的常规地震动El Centro波和包含永久地面位移的集集地震TCU049波作为地震动激励，分别采用上述2种输入模型进行结构地震响应分析。
El Centro波和集集地震TCU049波对应2种输入模型的加速度时程和位移时程分别如图3，4所示，其中g为重力加速度。图3（b）和图4（b）所示位移时程分别由图3（a）和图4（a）所示加速度时程经基线修正后积分得到。对于包含永久地面位移的地震动时程，采用改进基线修正方法[8]以满足速度时程末尾部分速度为0，位移时程末尾部分与时间轴基本平行的基线修正准则[12，19]。
4 计算结果及分析
限于篇幅，本文未将各墩墩底内力及加速度、位移时程比较全部罗列，仅列出了1#墩和3#墩的墩底弯矩响应时程和墩顶位移时程比较。由于对位移输入模型直接求解所得计算结果未经过任何近似处理，因此将其作为精确解，相对误差表示加速度输入模型计算值相对于位移输入模型的差值。
4.1 El Centro波作用下的结构地震响应
表1为2种输入模型下各墩墩底内力及墩（塔）顶加速度、位移最大值比较。图5，6分别为1#墩和3#墩的墩底弯矩响应和墩顶位移时程比较。
由表1可知：2种输入模型作用下各墩墩底内力及墩（塔）顶加速度、位移最大值较为接近，墩底内力最大相对误差出现在墩底轴力响应，相对误差为墩底内力和墩顶加速度、位移最大值比较。由表2可知，与常规地震动激励的计算结果不同，2种输入模型各墩墩底内力和墩顶加速度、位移最大值差别较大，除2#墩个别计算响应参数外，位移输入模型计算响应最大值总体大于加速度输入模型，且相差较大，最大相对误差达23.9%。 的墩底弯矩响应和墩顶位移时程比较。从结构响应看，断层错动之前，2种输入模型作用下1#墩和3#墩的墩底弯矩响应和墩顶位移时程的大小和曲线形状基本一致。断层错动后，位移输入模型的结构响应向不同方向偏离原振动平衡位置，随后沿新的平衡位置往复振荡，地震动结束后存在残余弯矩和位移。加速度输入模型在断层错动后弯矩响应和位移响应仍沿原振动平衡位置振荡，地震动结束后变形和内力值趋于0。
4.3 2种输入模型的适用性分析
由以上不同特性地震动输入下的计算结果可得到以下结论：
（1）在以常规地震动（未包含永久地面位移）作为输入对跨断层桥梁进行结构地震响应分析时，位移输入模型和改进大质量法加速度输入模型所得到的结构响应规律、趋势基本一致，数值非常接近，最大相对误差仅为3.2%，满足工程计算精度小于5%的要求。显然，2种方法均适用于常规地震动作用下的多点激励地震响应分析。
（2）在以包含永久位移地震动作为输入对跨断层桥梁进行结构地震响应分析时，位移输入模型和改进大质量法加速度输入模型所得到的结构响应规律、趋势及响应峰值存在显著差异。位移输入模型可考虑断层错动导致的永久地面位移对结构地震响应的影响，能够真实反映跨断层桥梁在地面运动结束之后桥墩具有的残余内力和变形，符合跨断层桥梁震害特征[7]，适用于跨断层桥梁的多点激励响应分析。加速度输入模型对此无法考虑，可能导致不合理也是不安全的计算结果。
4.4 大质量法加速度输入模型误差原因分析
由大质量法求解结构动力响应平衡方程的推导过程可知，尽管直观上大质量法的地震动输入形式为基底加速度（g），但其本质仍为通过基底位移（ug）求解结构地震响应[式（11）]。若将大质量法输入模型动力平衡方程式（11）中的约等号用等号代替，则与位移输入模型结构动力平衡方程式（3）一致。采用大质量法加速度输入模型计算地震响应时，gug的计算过程是导致地震响应计算出现误差的关键环节，而非大质量法本身。
图9为不含永久地面位移的El Centro波分别在2种输入模型下计算所得到的1#墩和3#墩墩底位移时程。由图9可知，2种输入模型的墩底位移时程基本相同，故2种输入模型计算所得的内力响应相同。对于未包含永久地面位移的常规地震动，大质量法加速度输入模型在gug过程中，计算结构内力响应的墩底位移时程实质上是由输入墩底的加速度时程经数值积分直接得到的，这与文献[12]的研究结论一致。
然而，对于包含有永久地面位移的地震动，这一结论并不一定适用。图10为包含永久地面位移的TCU049波分别在2种输入模型下计算所得到的墩和3#墩墩底位移时程。由图10可知，2种输入模型的墩底位移时程波形存在相似之处，但在断层错动后2种输入模型的变化趋势不同，大质量法基底输入在由加速度时程转化为位移时程的过程中滤去了永久地面位移部分，致使地震动结束后位移时程趋于0，故2种输入模型计算所得的内力响应存在差异。
综合上述分析可得如下结论：采用大质量法加速度输入模型对包含永久地面位移的地震动激励进行计算时，结构响应误差并不是大质量法本身引起的，而是由于SAP2000有限元分析软件内置计算程序未考虑包含永久地面位移的情况，在gug数值积分过程中滤去了永久地面位移部分，致使地震动结束后位移时程趋于0。
基于上述分析，在运用通用有限元计算软件进行加速度输入模型多点激励分析时，应首先对拟输入的加速度时程进行数值积分，确认对应的位移时程是否存在永久地面位移。若存在永久地面位移，则建议改为位移输入模型进行计算，否则将可能得到失真的计算结果。5 结 语
（1）在对跨断层桥梁进行非一致激励分析时，位移输入模型可考虑断层错动导致的永久地面位移对结构地震响应的影响，能够真实反映跨断层桥梁在地面运动结束之后桥墩具有的残余内力和变形，符合实际震害特征，适用于跨断层桥梁的非一致激励响应分析。加速度输入模型对此无法考虑，可能导致不合理也是不安全的计算结果。
（2） 在运用通用有限元计算软件进行加速度输入模型多点激励分析时，应首先对拟输入的加速度时程进行数值积分，确认对应的位移时程是否存在永久地面位移。若存在永久地面位移，则建议改为位移输入模型进行计算，否则将可能得到失真的计算结果。
参考文献：
References：
[1] PARK S W，GHASEMI H，SHEN J，et al.Simulation of the Seismic Performance of the Bolu Viaduct Subjected to Nearfault Ground Motions[J].Earthquake Engineering & Structural Dynamics，2004，33（13）：12491270. [2]UCAK A，MAVROEIDIS G P，TSOPELAS P.The Response of the Seismically Isolated Bolu Viaduct Subjected to Fault Crossing[C]//LESHKO B J，MCHUGH J.Structures Congress 2013：Bridging Your Passion with Your Profession.Reston：ASCE，2013：817826. [3]GOEL R K，CHOPRA A K.Linear Analysis of Ordinary Bridges Crossing Faultrupture Zones[J].Journal of Bridge Engineering，2009，14（3）：203215. [4]GOEL R，QU B，TURES J，et al.Validation of Fault Ruptureresponse Spectrum Analysis Method for Curved Bridges Crossing Strikeslip Fault Rupture Zones[J].Journal of Bridge Engineering，2014，19（5）：644651. [5]SAIIDI M S，VOSOOGHI A，CHOI H，et al.Shake Table Studies and Analysis of a Twospan RC Bridge Model Subjected to a Fault Rupture[J].Journal of Bridge Engineering，2014，19（8）：182190. [6]杨怀宇，李建中.断层地震动对隔震桥梁地震响应的影响[J].同济大学学报：自然科学版，2015，43（8）：11441152.YANG Huaiyu，LI Jianzhong.Response Analysis of Seismic Isolated Bridge Under Influence of Faultcrossing Groundmotions[J].Journal of Tongji University：Natural Science，2015，43（8）：11441152.
[7]惠迎新，王克海，李 冲.跨断层地表破裂带桥梁震害研究及抗震概念设计[J].公路交通科技，2014，31（10）：5157.HUI Yingxin，WANG Kehai，LI Chong.Study of Seismic Damage and Seismic Conceptual Design of Bridges Across Fault Surface Rupture Zones[J].Journal of Highway and Transportation Research and Development，2014，31（10）：5157.
[8]惠迎新，王克海.基于多点激励位移输入模型的跨断层桥梁地震动输入方法[J].东南大学学报：自然科学版，2015，45（3）：557562.HUI Yingxin，WANG Kehai.Earthquake Motion Input Method for Bridges Crossing Fault Based on Multisupport Excitation Displacement Input Model[J].Journal of Southeast University：Natural Science Edition，2015，45（3）：557562.
[9]惠迎新，王克海，吴 刚，等.跨断层桥梁地震响应分析及合理跨越角度研究[J].振动与冲击，2015，34（13）：611，17.HUI Yingxin，WANG Kehai，WU Gang，et al.Seismic Responses of Bridges Crossing Faults and Their Best Crossing Angles[J].Journal of Vibration and Shock，2015，34（13）：611，17.
[10]柳国环，李宏男，田 利.九江长江大桥在多点多维地震激励下的反应分析[J].振动与冲击，2009，28（9）：204209，224.LIU Guohuan，LI Hongnan，TIAN Li.Response Analysis of Jiujiang Yangtze River Highway Bridge Under Spatially Variable Earthquake Ground Motions[J].Journal of Vibration and Shock，2009，28（9）：204209，224.
[11]杨庆山，刘文华，田玉基.国家体育场在多点激励作用下的地震反应分析[J].土木工程学报，2008，41（2）：3541.YANG Qingshan，LIU Wenhua，TIAN Yuji.Response Analysis of National Stadium Under Specially Variable Earthquake Ground Motions[J].China Civil Engineering Journal，2008，41（2）：3541.
[12]罗 超，楼梦麟，桂国庆.大跨度结构多点地震反应计算方法的比较[J].同济大学学报：自然科学版，2015，43（1）：815.LUO Chao，LOU Menglin，GUI Guoqing.Comparison for Calculation Methods of Longspan Structure Under Multisupport Seismic Excitation[J].Journal of Tongji University：Natural Science，2015，43（1）：815.
[13]周国良，李小军，喻 x，等.结构地震作用的基底激励模型及其适用性[J].建筑结构学报，2010，31（增2）：8288.ZHOU Guoliang，LI Xiaojun，YU Tian，et al.Applicability Research on Base Excitation Models Used in Structural Seismic Response Analysis[J].Journal of Building Structures，2010，31（S2）：8288.
[14]CHOPRA A K.Dynamic of Structures：Theory and Applications to Earthquake Engineering[M].3rd ed.Upper Saddle River：Prentice Hall，2007. [15]田玉基，杨庆山.地震地面运动作用下结构反应的分析模型[J].工程力学，2005，22（5）：170174.TIAN Yuji，YANG Qingshan.Analysis Models and Methods for Structural Seismic Responses[J].Engineering Mechanics，2005，22（5）：170174.
[16]LEGER P，IDE I M，PAULTRE P.Multiplesupport Seismic Analysis of Large Structures[J].Computers & Structures，1990，36（6）：11531158. [17]SOMERVILLE P G，SMITH N F，GRAVES R W，et al.Modification of Empirical Strong Ground Motion Attenuation Relations to Include the Amplitude and Duration Effects of Rupture Directivity[J].Seismological Research Letters，1997，68（1）：199222. [18]PARK S W，GHASEMI H，SHEN J，et al.Simulation of the Seismic Performance of the Bolu Viaduct Subjected to Nearfault Ground Motions[J].Earthquake Engineering & Structural Dynamics，2004，33（13）：12491270. [19]于海英，江汶乡，解全才，等.近场数字强震仪记录误差分析与零线校正方法[J].地震工程与工程振动，2009，29（6）：112.YU Haiying，JIANG Wenxiang，XIE Quancai，et al.Baseline Correction of Digital Strongmotion Records in Nearfield[J].Journal of Earthquake Engineering and Engineering Vibration，2009，29（6）：112.