基于规范反应谱编制的matlab 选波程序
一、选波程序简介
本程序采用Newmark-beta 法求解单自由度弹性体系运动方程，地震波库来源于PEER 强震数据库（旧版），目前共收集约3600多条地震动数据。
从实际工程应用出发，编制该matlab 程序主要目的是依据规范反应谱来选波，力求能够从强震数据库中选出满足规范要求的地震动。在《抗规》第5.1.2条中提到，进行时程分析所使用的地震动主要满足以下两点：（1）与规范应谱所采用的地震影响系数曲线在统计意义上相符；（2）结构基底剪力要求。而本程序主要是满足第（1）点要求。
因此，由本程序选出的地震动并不一定能满足规范对结构基底剪力的要求。但通过本程序能够初步筛选出一批符合规范在统计意义上要求的地震动，减少盲目选波耗费时间。
如有需要选波，可通过邮件提供相关参数（详第三点），本人会尽快回复相应筛选结果。邮箱：185541812@https://www.doczj.com/doc/e4918647.html, 。同时也希望大家能反馈选波是否有效，以便进一步改善程序。
另注：本程序仅提供选波数据，对这些地震波应用于实际工程分析的结果并不负责，烦请各位自行判断结果合理性和正确性。 二、选波原则
目前主要提供以下两种选波方式：
（1）根据《抗规》要求，在结构主要周期点上，控制规范反应谱与地震动反应谱值相差不超过一定误差范围（如10%~20%的差异），如图1所示；
（2）根据杨溥和李英民提出的双频段控制方式选波。具体详《结构时程分析法输入地震波的选择控制指标》一文，如图2所示。
图1 选波方式一（主周期控制） 图2 选波方式二（双频段控制）
三、选波需提供输入参数
主要提供以下参数：①结构第一、二周期T1，T2；②结构阻尼比ζ；③场地特征周期Tg；
④水平地震影响系数最大值αmax；⑤时程分析地震动的峰值加速度；⑥规范反应谱与地震动反应谱在结构第一、二周期点上的误差控制（建议取10%~20%）。
四、输出结果
主要提供以下三个数据文件（以excel文件给出）：
（1）筛选出的地震动反应谱数据（图3）
蓝色框中数据依次是地震动编号、在T1和T2周期点上规范反应谱与地震动反应谱的误差；绿色框中数据是该地震动编号对应的反应谱（即水平地震影响系数α），该列数据是根据地震动的峰值加速度进行调幅后的结果；红色框中数据是该地震动反应谱，与水平地震影响系数α值相对应的周期（s）。
图3 程序选出的地震动反应谱数据
（2）筛选出地震动的原加速度时程（图4）
蓝色框中数据是地震动编号；绿色框中数据是该地震动编号对应的加速度，单位为cm/s2，该加速度数据是原波数据，并未调幅；红色框中数据是与加速度对应的时间，单位为s秒。
图4 程序选出的地震动反应谱数据
（3）地震动编号对应信息（图5）
图5 地震动编号对应信息
五、反应谱计算验证
通过与seismosignal 计算对比，验证本程序求解单自由度弹性体系运动方程的正确性。
图6 加速度谱对比 图7 速度谱对比
1
2
3
4
5
6
50
100150200250300
350加速度 (c m /s /s )
周期 (s)
0123456
50
100
150
200
速度 (c m /s )
周期 (s)
图8 位移谱对比
0123456
50
100
150
200
位移 (c m )
周期 (s)