小震弹性时程计算结果比反应谱法大_地震波转反应谱 附源码

建筑抗震设计规范(GB50011-2010，以下简称抗规)中5.1.2规定，对特别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围内的高层建筑，应采用时程分析法进行多遇地震下的补充计算；

相应的条文说明中指出：正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间均要符合规定。其中频谱特性的要求中，要求时程分析所选的地震波与抗规的反应谱在统计意义上相符，即多组时程波的平均地震响应系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在对应于结构主要振型的周期点上相差不大于20%。

非解构之前的文章中，已经对地震波的相关概念进行介绍(说说地震波的那些事儿，一文读懂地震波，值得收藏！，Python爬虫实战——从太平洋地震工程研究中心数据库自动下载地震波时程，地震来了到底该跑还是就近躲避？)。那么，要如何生成地震相应系数曲线呢？

1、单条地震波反应谱的生成

反应谱是某个反应量的峰值作为体系的固有振动周期Tn，或者像圆频率ωn(或者循环评率fn)那样的相关参数的函数图形。作为一种实用的工具来描述地面运动及其结构效应。它提供了一种方便的方法来概括所有可能的线性单自由度体系对地面运动的某个特定分量的峰值反应。反应谱是在1932年由M.A.Biot引入的，40年代由Housner应用到抗震中，50年代末刘恢先教授在我国第一部抗震规范编制时引入了反应谱理论，开创了反应谱在我国抗震规范中应用的先河。

下面我们就用一个实际的例子来解释到底什么是反应谱，它是如何生成的。我们用ELCENTRO波作为我们的案例对象，说明ELCENTRO波的反应谱的计算过程。以下是ELCENTRO波的时程曲线(以下求解过程假设阻尼比均为5%)：

回到最基本的求解单自由度体系时程反应的结构动力学方程：

其中m为质量，c为阻尼，u为反应位移，为反应速度，为反应加速度，P(t)为地震波时程。

首先，我们需要求解每条地震波作用在任意单自由度体系时，单自由度体系的反应时程。根基结构动力学的理论，对于任意激励的单自由度体系，其时程反应可以用杜哈梅积分(Duhamel)进行求解：

下图所示为小编利用杜哈梅积分求解将ELCENTRO波作用在不同的单自由度体系上的时程反应：

可以清晰地发现，ELCENTRO波作用在不同周期的单自由度体系时，每条时程的反应的峰值是不同的，其产生峰值的时间点也各不相同 。

假设反应谱的周期范围为0-6s。只要我们按照一定的将这个周期范围按一定的时间步长进行划分，分别求解每个周期点对应的单自由度体系在ELCENTRO波作用下的反应。提取每条反应时程的峰值，以周期点为横轴，以每个周期点对应的时程反应峰值为纵