一种匹配多目标的频率非平稳地震动拟合方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN111551994A(43)申请公布日2020.08.18(21)申请号202010393272.8(22)申请日2020.05.11(71)申请人中国地震局地球物理研究所地址100081北京市海淀区民族大学南路5号(72)发明人俞瑞芳王少卿(74)专利代理机构北京煦润律师事务所11522代理人惠磊(51)Int.Cl.G01V1/36(2006.01)G01V1/30(2006.01)权利要求书2页说明书5页附图2页(54)发明名称一种匹配多目标的频率非平稳地震动拟合方法(57)摘要本申请涉及一种匹配多目标的频率非平稳地震动拟合方法，包括步骤：根据目标地震动拟合参数计算演变功率谱时频包线函数；选定初始相位参数，合成初始地震动并计算小波时频谱；将目标加速度反应谱离散化处理，对全部控制点进行多轮调整并进行基线校正，直至地震反应谱和峰值加速度都满足目标限值；记录地震反应谱调整轮数，并计算满足反应谱和峰值加速度限值要求的地震动小波谱与初始地震动小波谱之间的二维相关系数；若不满足限值，则重新生成新的初始相位参数，并重复上述步骤。该方法模拟的地震动不仅能够同时匹配多个目标参数的要求，而且还能够表征特定的目标强度和频率非平稳特性。CN111551994ACN111551994A权利要求书1/2页1.一种匹配多目标的频率非平稳地震动拟合方法，其特征在于：包括以下步骤：(1)根据目标地震动拟合参数计算演变功率谱时频包线函数；(2)选定初始相位参数，合成初始地震动并计算小波时频谱；(3)将目标加速度反应谱离散化处理，对全部控制点进行多轮调整并进行基线校正，直至地震反应谱和峰值加速度都满足目标限值；(4)记录地震反应谱调整轮数，并计算满足反应谱和峰值加速度限值要求的地震动小波谱与初始地震动小波谱之间的二维相关系数；(5)若地震反应谱调整轮数和上述二维相关系数不满足限值，则重新生成新的初始相位参数，并重复步骤(2)-(4)，直至地震动反应谱调整轮数和上述二维相关系数满足限值要求，从而获得满足所有目标的地震加速度时程。2.根据权利要求1所述的频率非平稳地震动拟合方法，其特征在于：在步骤(1)中，目标地震动拟合参数包括地震动主频率变化、地震动强度包线、目标加速度反应谱和峰值加速度，上述目标地震动拟合参数根据场地条件、地震地质条件分析得到。3.根据权利要求1或2所述的频率非平稳地震动拟合方法，其特征在于：在步骤(1)中，计算演变功率谱时频包线函数的方法包括：利用给定地震动主频率变化函数计算频率成分随时间的分布函数，结合地震动强度包线得到时频谱等效平均强度包线，从而得到地震动的时频联合分布和地震动的演变功率谱时频包线函数。4.根据权利要求3所述的频率非平稳地震动拟合方法，其特征在于：在步骤(1)中，给定地震动主频率变化Fp(t)，用下式的模型描述频率成分随时间t的分布：结合地震动强度包线E(t)，可得时频谱等效平均强度包线2式中ω为自振圆频率；Ck＝max[E(t)L(t,ωk)]；S(ω)为目标地震动的近似功率谱，由目标加速度反应谱计算得到：式中：ξ为结构阻尼比，T为结构自振周期，P为超越概率；由此可得地震动的时频联合分布为：地震动的演变功率谱时频包线函数为：5.根据权利要求1-4中任一项所述的频率非平稳地震动拟合方法，其特征在于：在步骤(3)中，具体步骤包括：2CN111551994A权利要求书2/2页将目标加速度反应谱离散为M个单自由度体系；对全部M个控制点进行第i轮调整；将第i轮调整后的时程峰值调整为目标加速度峰值，并进行基线校正，然后作为第i+1轮调整的初始地震动进行第i+1轮调整，直至地震动反应谱和峰值加速度都满足目标限值。6.根据权利要求4所述的频率非平稳地震动拟合方法，其特征在于：在步骤(3)中，具体步骤包括：(1)将目标加速度反应谱离散为圆频率为ω1,ω2,…,ωM的M个单自由度体系；(2)对全部M个控制点进行第i轮调整：第i轮反应谱调整中，若对第k个圆频率ωk调整前的加速度时程为则调整后的加速度时程为：式中，为第i轮调整中对ωk进行调整的增量时程，采用下式计算：上式中，tm为单自由度体系最大反应发生的时刻；B(t,ωk)为演变功率谱时频包线函数；h(tm-t)为单位脉冲响应函数，表示为：其中：为调整时程的初始相位，保证叠加的单位脉冲函数在tm处取得峰值；第i轮调整中，若此次调整前地震动的反应谱在ωk与目标反应谱的差值为则公式(8)中的增量时程振幅调整系数R由下式计算其中：(3)将第i轮调整后的时程峰值调整为目标加速度峰值，并进行基线校正，然后作为第i+1轮调整的初始地震动进行第i+1轮调整，直至地震反应谱和峰值加速度都满足目标限值。7.一种工程结构抗震设计的方