手册-时程_地震波

时程反应谱换算时程是描述地震动运动的时间历程。

地震动是地震引起的地面振动，其运动特征可以通过时程来描述。

时程可以用于评估结构物的地震响应，是地震工程中重要的参数之一。

时程可以通过观测实际地震事件得到，也可以通过合成地震动进行模拟获得。

观测到的地震动通常以地震波形的形式进行记录，地震波形是时域的表达方式。

地震波形的单位是加速度、速度或位移，时间单位是秒。

地震波形可以反映地面在某一特定位置上的振动情况。

在地震工程中，时程是用来评估结构物抗震性能的重要参数之一。

通过将地震波形输入结构物的动力分析模型，可以计算结构物的地震响应，如位移、加速度、应力等。

对于某一特定地震事件来说，不同位置上的结构物受到的地震动可能会有所不同，因此需要根据不同位置的地震动记录得到相应的时程数据。

时程分析可以有线性和非线性两种模型。

线性时程分析假设结构物的反应是线性的，即结构物的刚度和阻尼在地震作用下不发生变化。

非线性时程分析考虑结构物在地震作用下发生的非线性行为，如屈曲、破坏等。

根据结构物模型的复杂程度和分析的目的，可以选择使用线性时程分析或非线性时程分析方法。

反应谱是描述结构物地震响应的频域表达方式。

反应谱可以通过将时程信号进行傅里叶变换得到。

在反应谱中，横坐标代表结构物的周期，纵坐标代表结构物的地震反应。

地震反应可以是加速度、速度或位移的峰值值或响应谱。

反应谱可以用来研究结构物在不同周期下的地震反应情况，从而评估结构物的抗震性能。

反应谱分析是地震工程中常用的一种分析方法。

通过反应谱分析，可以得到结构物在不同周期下的地震反应谱，进而评估结构物的抗震性能。

通常，设计地震动会以一种反应谱的形式给定，结构物的设计需满足给定反应谱的限制条件。

反应谱和时程之间存在一种转换关系。

时程可以通过反傅里叶变换获得频域信息，从而得到反应谱。

反应谱可以通过傅里叶变换得到时程信息，从而得到结构物在地震作用下的运动情况。

反应谱和时程相互转换的过程中，会有信息的丢失，因此需要根据分析的目的和应用的要求来选择使用时程还是反应谱来描述地震动。

时程分析计算精辟解读（值得收藏）时程分析法是20世纪60年代逐步发展起来的抗震分析方法.用以进行超高层建筑的抗震分析和工程抗震研究等.至80年代,已成为多数国家抗震设计规范或规程的分析方法之一.“时程分析法”是由结构基本运动方程输入地震加速度记录进行积分,求得整个时间历程内结构地震作用效应的一种结构动力计算方法,也为国际通用的动力分析方法.“时程分析法”常作为计算高层或超高层的一种（补充计算）方法,也就是说满足了规范要求的时候是可以不用它计算结构的.规范规定：对于特别不规则的建筑、甲类建筑及超过一定高度的高层建筑,宜采用时程分析法进行补充计算.所以有较多设计人员对应用时程分析法进行抗震设计感到生疏.近年来,随着高层建筑和复杂结构的发展,时程分析在工程中的应用也越来越广泛了.1输入地震动准则输入地震动准则即为结构时程分析选择输入地震加速度记录时程（简称地震波）的基本要求,包括：地震环境（场地类别和地震分组）、数量、持续时间、检验方法等.地震波的合理选择是时程分析结果能否既反映结构最大可能遭遇的地震作用,又满足工程抗震设计基于安全和功能要求的基础.在这里不提“真实”地反映地震作用,也不提计算结果的“精确”性,正是基于对结构可能遭遇地震的极大不确定性和计算中结构建模的近似性.在工程实际应用中经常出现对同一个建筑结构进行时程分析时,由于输入地震波的不同,造成计算结果的数倍乃至数十倍之差,使工程师无所适从.《建筑抗震设计规范》（GB50011—2010）（简称2010规范）5.1.2-3条要求“采用时程分析法时,应按建筑场地类别和设计地震分组选用实际强震记录和人工模拟的加速度时程曲线,其中实际强震记录的数量不应少于总数的2/3,多组时程曲线的平均地震影响系数曲线（即反应谱）应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符”.1.1“选波”要求1.1.1地震环境要求2010规范在构建设计反应谱时,按不同场地类别和震级、震中距从全球强震加速度记录数据库中挑选了数百条地面加速度记录,求出每条记录的反应谱.同时收集这些记录台站的地质剖面和地震震级、震中距等参数,按照2010规范的场地类别划分标准,场地分成Ⅰ~Ⅳ类和远、中、近震分组,共计12组,再经平滑处理得到2010规范5.1.5条的地震影响系数曲线,即设计反应谱.时程分析法输入地震波的选择应遵循上述构建设计反应谱的原则,考虑建设场地与记录台站场地的地震环境.1.1.2数量要求对于高度不是太高、体型比较规则的高层建筑,取2+1,即选用不少于2条天然地震波和1条拟合目标谱的人工地震波,计算结果宜取包络值.对于超高、大跨、体型复杂的建筑结构,取5+2,即不少于7组地震波,其中,天然地震波数量不少于总数的2/3,计算结果取平均值.1.1.3持续时间要求为了充分地激励建筑结构,一般要求输入的地震动有效持续时间为结构基本周期的5倍左右.对于结构动力时程分析,只有加速度记录的强震部分的时长,即有效持续时间才有意义.最常用的有效持续时间定义是：取记录最大峰值的10%~15%作为起始峰值和结束峰值,在此之间的时间段为有效持续时间.图1表示编号为US185地震加速度记录的波形,用于7度小震下结构时程分析,最大加速度峰值是35gal,取首、尾两个峰值为3.5gal之间的时间长度为有效持续时间,大约为30s,可用于基本周期小于6s的结构.ps:持续时间不是指整个时程的记录时间图1加速度记录有效持续时间的定义1.1.4统计特征要求规范规定,时程分析所采用的地震波的平均反应谱与振型分解反应谱法所采用的反应谱应“在统计意义上相符”.如前所述,天然地震波具有千变万化的特征,不同结构的动力特性也千差万别.对同一个结构,输入不同的地震波进行时程分析会得到完全不同的结果.所以,遵循“在统计意义上相符”的原则选择天然地震波时,只要求所选的天然地震加速度记录的反应谱值在对应于结构主要周期点（而不是每个周期点）上与规范反应谱相差不大于20%.这个要求只是一种参考,便于数据库管理员在数据库中挑选合适的记录.一般情况下,照此要求选择的地震波可以满足时程分析要求.但是,不宜将此作为检验地震波的标准,检验标准仍然是规范规定的结构底部剪力.为什么既要求有天然地震波,又要求有人工地震波作为输入？原因是,所谓人工地震波,是应用数学方法,将足够多的具有不同周期的正弦波叠加组合形成一个平稳或非平稳的随机时间历程,对叠加组合过程不断进行迭代修正,使它的反应谱逐步逼近规范的设计反应谱.当拟合精度达到在各个周期点上的反应谱值与规范反应谱值相差小于10%,即认为“在统计意义上相符”了.这样合成的人工地震波具有足够多的周期分量,可以均匀地“激发”结构的各个振型响应.但是,由于人工地震波是“拟合”设计反应谱的加速度时间过程,不具备天然地震波的完全非平稳随机过程特性,特别是缺少强烈变化的短周期成分.因此它只能在设计反应谱的“框架”内激励结构,无法“激发”结构的高阶振型响应,所以时程分析要求以天然地震波为主,同时辅以人工地震波作为地震动输入.ps:人工波对低阶振型激发较好,而对于高阶振型的激发不够（如肖总所说）,因此对于高阶振型部分,必须仰仗天然波来激发.本人理解,作者建议采用EPA,就是为了保证天然波对于高阶振型的激发.弹性时程分析与振型分解反应谱分析的关系,实质上是事物的特殊性与一般性的关系,多条地震波时程分析结果的平均值近似于反应谱法计算结果,输入的地震波数量越多,这种近似性越好.ps:现在很多软件能够根据规范相关要求,自动选波,比如YJK弹性时程分析时就可以做相关的自动选波.自动筛选最优地震波组合这块就给设计师在筛选地震波这块提供了相当大的便利.选择框中列出了程序根据特征周期归类后的波库中天然波和人工波,用户可从中选择参与筛选的备选地震波到中间列表框.如全选,筛选出的地震波组合可能多一些,但计算时间稍长.可根据规范在对话框下部设定地震波组合的人工波数,天然波数.按照规范要求,实际强震记录的数量不应少于总数的2/3.若选用不少于二组实际记录和一组人工模拟的加速度时程曲线作为输入,计算的平均地震效应值不小于大样本容量平均值的保证率的85%以上.YJK计算程序即根据设置好的限定条件计算每条地震波的基底剪力与结构周期点上所对应的反应谱值.最终满足要求的所有组合结果将在该按钮下方的列表框中按最优至次优的顺序显示.列表中的组合可以通过选择地震波组合按钮选择,选中的地震波组合包含的地震波将在下方列表框中显示.如下图所示图1自动筛选最优地震波组合对话框根据《建筑抗震设计规范》GB50011-2010中的规定,程序遵循的地震波组合筛选原则如下：1：单条地震波满足限制条件每条地震波输入的计算结果不会小于65%,不大于135%.2：多条地震波组合满足限制条件（1）“在统计意义上相符”,即多组时程波的平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比,在对应于结构主要振型的周期点上相差不大于20%,即：>80%并<120%（2）多条地震波计算结果在结构主方向的平均底部剪力一般不小于振型分解反应谱计算结果的80%,不大于120%.（3）按照平均底部剪力与振型分解反应谱法计算的底部剪力偏差最小的原则对已经满足上述限制的组合再进行排序,默认选出偏差最小的组合作为最有组合.在搜索过程中,当程序提示未搜索到符合要求的地震组合时,用户可根据抗震规范规定适当增加相邻特征周期的可选地震波或者放宽主次方向地震峰值加速度值以满足以上的限制条件.点击“查看计算结果文本”,程序将打开结果文件,内容包括了地震波在筛选地震波组合时计算的地震波基底剪力,周期点谱值及地震波组合计算的统计结果.用户可根据该计算统计结果适当改变地震波组合方案.图2筛选方案排序示例对于未筛选出满足要求的地震波组合工程,用户可尝试从以下几方面检查参数设置或者进行适度调整.（1）主次波峰值加速度对应地震烈度是否与前处理中地震烈度设置一致.（2）前处理中周期折减系数是否过小.（3）根据规范相关阐述,在选取不到恰当地震波组合情况下,可选取相邻特征周期地震波或增大减小地震波峰值加速度以满足剪力即谱值要求.1.2天然地震波加速度值的调整如前所述,结构时程分析法补充计算被用于校核振型分解反应谱法的计算结果.反应谱法以反应谱作为输入地震动,时程分析以加速度时程（地震波）作为输入,需要对它的加速度值进行调整.2010规范以中国地震动参数区划图定义的地面峰值加速度GPA为设防地震（中震）基本地震峰值加速度,如表3.2.2所示；表5.1.2-2分别给出多遇地震（小震）和罕遇地震（大震）加速度峰值,与之相对应的规范设计反应谱是基于大量的天然地震加速度记录,并经平滑处理和统计平均后构建的,是地震动的预期均值.对每一条天然地震波加速度时程进行调整的步骤是：根据规范给定的加速度峰值GPA,按比例调整后求得其加速度反应谱,经平滑处理得到归一化的反应谱.运用式（1）求得有效峰值加速度EPA,以其为基准对地震波加速度时程进行再调整,得到结构时程分析所需要的加速度时程.需要指出的是,有效峰值加速度EPA不等于地面峰值加速度GPA,当地震波的短周期成分显著时,GPA大于EPA.如前所述,人工地震波是采用拟合规范反应谱的数值合成方法得到的加速度时程,按GPA比例调整后即可作为时程分析的输入地震动.美国地震危险区划图定义,有效峰值加速度EPA、加速度反应谱最大值Sa（对应于中国规范的地震影响系数）和放大系数β存在如下关系：式中：Sa（0.2）为周期0.2s处的谱加速度值；β为动力放大系数,取2.50（中国规范取2.25）.下面以位于7度区III类场地（Tg=0.70s）的设计地震分组为第三组的某一高层建筑为例,大震作用下结构弹塑性时程分析选用7组输入地震波,其归一化的加速度反应谱及其平均、平滑处理后的结果如图2所示.图2地震加速度反应谱表1和图3是每条地震波调整前后的地震动参数与规范的对比（大震作用GPA 取2.20m/s2）.可以看到,尽管各条地震波的三个参数差别较大,但经平滑平均后接近于规范反应谱,且EPA＜GPA.由此也可证明,2010规范对地震波数量的要求是必要而且合理的.ps:上表的平均值为平均谱所得的值.如amax,并不是每个波amax的平均,而是由平均谱求得的,因为每个谱的极值点不会都在同一个周期,故平均谱的amax比每个波amax的平均值小.图3地震动参数对比1.3检验要求《建筑抗震设计规范》（GB50011—2001）（简称2001规范）和2010规范提出：弹性时程分析时,每条时程曲线计算所得底部剪力不应小于振型分解反应谱法计算结果的65%,多条时程曲线计算所得底部剪力的平均值不应小于振型分解反应谱计算结果的80%.具体操作时,当采用一组（单向或两向水平）地震波输入进行时程分析,结构主方向基底总剪力为同方向反应谱CQC计算结果的65%～130%,多组地震波输入的平均值为80%～120%.不要求结构主、次两个方向的基底剪力同时满足这个要求.需要说明的是,对结构可以按第一、二主振型认定主、次方向,而一组地震记录的两个水平方向无法区分主、次方向.ps:X向为主时仅要求X向满足,Y向为主是仅要求Y向满足.1.4选波实例下面以两组天然地震波和一组人工合成地震波为例说明选波过程及效果.（1）图4为所选择的一组3分量天然地震波时程及反应谱,其中编号US2570和US2571为两向水平分量,US2569为竖向分量,需要按小震作用所对应的最大加速度峰值进行调整,除有特殊要求外,通常取两向水平峰值与竖向峰值之比为1.00：0.85：0.65.从波形和反应谱可以看到,竖向分量的短周期成分十分显著,水平分量在短周期部分的波动也很显著,各向分量的反应谱曲线相差明显.图4第一组天然地震波和反应谱（2）图5为另一组3分量天然地震波时程及反应谱,其中编号US184和US185为两向水平分量,US186为竖向分量.同样可以看到,竖向和水平分量在短周期部分的波动很明显,但是两个水平分量的反应谱曲线比较一致.两组地震波反应谱的明显差异反映了天然地震波特征的不确定性,用于结构时程分析时,很难做到两向水平输入的地震波均能满足规范要求,一般只要求结构主方向的底部总剪力满足规范要求即可.图5第二组天然地震波和反应谱（3）图6为三条人工地震波及反应谱.图6三条人工地震波及反应谱2时程分析输出结果解读结构时程分析一般要求进行小震作用下弹性和大震作用下弹塑性计算.对计算结果的解读可以判断结构的动力响应和损伤情况.2.1小震作用下的计算结果（1）楼层水平地震剪力分布：对于高层建筑,通常可由此判断结构是否存在高阶振型响应并发现薄弱楼层.图7为某幢高层建筑结构小震弹性时程分析得到的楼层剪力分布,可见结构存在高阶振型响应,应对结构上部相关楼层地震剪力加以调整放大.图7楼层地震剪力分布（2）弹性层间位移角分布：如图8所示,上部结构部分楼层的层间位移角大于规范限值.从图7和图8可以看到,输入3组地震波进行时程分析,结构高阶振型响应明显,上部楼层剪力和位移均放大了,应对反应谱法结果进行调整,采用包络设计.图8弹性层间位移角分布2.2大震作用下的计算结果（1）层间位移角分布：按照规范要求进行大震作用下结构的时程分析,主要是弹塑性变形计算,力的计算并不重要.计算结果通常给出弹性和弹塑性层间位移角分布的对比,如图9所示.X向最大层间位移角为1/178,Y向为1/138,均满足规范限值1/100.一般情况下,最大弹性位移角大于弹塑性位移角.图9弹塑性层间位移角分布（弹塑性/弹性）（2）结构顶点位移时程曲线：从结构顶点位移时程曲线除了可以看出位移是否满足规范限值外,更重要的是可以判断结构整体刚度退化程度,并推测结构的塑性损伤程度.如图10所示,弹塑性位移时程曲线表明,结构的周期逐步变长,说明有部分构件累积损伤,导致结构整体刚度退化.图10结构顶点位移时程曲线对比（弹性/弹塑性）（3）构件损伤：通常要求给出主要抗侧力构件,如剪力墙、框架柱、支撑、环带桁架、伸臂桁架等,以及耗能构件,如连梁、框架梁等的损伤,以应力比、应变、损伤因子等表示.图11表示某高层建筑核心筒剪力墙受压、受拉和框架柱的损伤.图11核心筒剪力墙和框架柱损伤（4）能量分布：有的软件可以提供在地震作用下结构的能量分布情况.如图12所示,从上至下的区域分别表示结构动能、弹性应变能、与质量M相关的粘滞阻尼耗能、与刚度K相关的粘滞阻尼耗能、塑性耗能.其中,塑性耗能属于不可恢复的能量耗散,所占比例越大,表明结构整体破坏越严重.图12结构能量分布。

地震动时程曲线时域校正的一种方法胡汛训;马雪峰【摘要】在做土-结动力相互作用分析时,需要在地基边界上输入速度及位移时程曲线.速度及位移时程往往由加速度时程简单积分得到.由于加速度时程中存在一些长周期分量,这些长周期分量对加速度时程的影响甚微,却会造成积分后的位移时程出现漂移现象,从而影响边界输入数据的真实.为此,通过探索表明,采用最小二乘法编制了位移时程漂移时域校正程序,得到了满意的结果.【期刊名称】《水电与新能源》【年(卷),期】2017(000)007【总页数】3页(P17-19)【关键词】地震动;基线漂移;最小二乘法【作者】胡汛训;马雪峰【作者单位】湖北省水利水电规划勘测设计院,湖北武汉 430072;湖北省水利水电规划勘测设计院,湖北武汉 430072【正文语种】中文【中图分类】P315随着工程抗震研究手段的进步，学者们在进行桁架、桥梁、大坝等一些大跨度空间结构抗震分析时发现，土-结动力相互作用对结构的影响所占比重较大，是不可忽略的部分。

在做土-结动力相互作用工程抗震分析时，因为边界处理的需要，地基边界上输入的是地震动的速度和位移数值，而不是加速度数值。

理论上，对加速度逐次进行积分运算，即可得到速度和位移，但是，按次积分计算得到的速度和位移时程往往会出现基线漂移的现象，即后期的振动偏离水平零轴。

这种现象可以从图1中看到，图1中加速度和速度曲线围绕坐标轴上下振动，最终归零，没有发生漂移现象。

但是位移曲线随着时间的推移，振动的平衡位置逐渐偏离坐标轴，产生了“残余变形”的现象。

产生这种现象的原因就是加速度波中存在低频分量，经过两次积分后，这种不适宜分量得到放大，造成了积分后的位移在最终时刻非零，即产生的基线漂移，从而造成了抗震分析中位移输入数据的失真，影响了结构地震动分析。

地震波时程校正一般可以分为两类方式，其一是采用了数字滤波校正，将地震波加速度时程中不适宜的波频剔除；二是通过直接修改地震波加速度时程的初期数值，强制积分后的位移时程终点时刻归零。

第二章 地震波的时距曲线在地震勘探工作中，每激发一次人工地震，都要在多个检波点接收地震信号。

炮点和检波点都沿一条直测线布置，炮点到任意检波点的距离称炮检距x ，相邻检波点的距离叫道间距Δx ，来自同一界面的地震波沿不同路径先后到达各检波点，从而形成一张如图所示的地震记录。

图中横坐标表示地震波旅行时间t ，纵坐标表示炮点到任意检波点的距离称炮检距x ，每一条波动曲线是一道地震记录，它反映出一个检波点的振动过程。

来自同一界面的反射波（或折射波）以一定的视速度规律依次到达个检波点，在地震记录中表现为振动极值的规则排列，各道地震记录波按一定规则排列，形成同相轴（它是相同相位点的连线形成的图形）。

同相轴反映出地震波的旅行时间t 与炮检距x 的函数关系。

将它表示在t-x 直角坐标系中，称为地震波的时距曲线。

不同种类的地震波，其时距曲线的形状不同。

如图中的直达波、反射波、折射波、地滚波、声波等都有自己特有的形状。

每一类特定的时距曲线，其曲线参数与地下介质的纵波速度v 及地震界面的产状有着直接的关系。

第一节 反射波的时距曲线一、 两层介质的直达波和反射波时距曲线（一）直达波的时距曲线从震源出发，不经过反射或折射而直线前进到各检波点的地震波成为直达波。

当震源深度为零时，直达波沿测线传播，旅行时间t 与炮检距x 的函数关系为)1.1.2(1v x t ±= 是两条经过原点的、斜率为1/v 1的两条直线。

如图2.1-1，根据直达波时距曲线的斜率，可以求取界面上层介质的波速v 1。

图2.1-1 直达波与水平界面反射波时距曲线(二)水平界面的反射波时距曲线和正常时差由图2.1-1，若界面埋深为h, 炮点0为激发点，到达界面R 点后反射到地面的s 点，设s 点的炮检距为x ，为计算方便，做炮点0关于界面的镜像点0*，称为虚震源，根据图2.1-1的几何关系，反射波旅行时间t 与炮检距x 的函数关系为)2.1.2(4102211*x h v v RS t +== 将反射波在炮点的反射时间称为反射回声时间，102v h t = 则（2.1.1）式可改写为)2.1.2()(2122022120′+=+=v x t t v x t t 或 式（2.1.2）就是水平界面反射波的时距曲线，可化简为以下的标准双曲线方程)2.1.2(1422202′′=−h x t t综上所述：1．反射波时距曲线在x-t 坐标系是双曲线，其极小点在炮点正上方；2．在x 2-t 2坐标系，反射波时距曲线是直线，直线的斜率为1/v 12, 利用直线的斜率可求界面上方介质的速度；3．反射波时距曲线以直达波时距曲线为其渐近线。