一种构造多参数化时域调整曲线的多阻尼人工波拟合方法技术

本发明专利技术属于核电设备工程抗震安全分析与评价技术领域，涉及一种构造多参数化时域调整曲线的多阻尼人工波拟合方法。该窄带时程构造方法，其在时域和频域内均具有窄带特性，且避免了由于傅里叶变换导致的振铃效应，这种新的窄带时程被称为标准输入时程；使用每个频率控制点每个阻尼比单独叠加标准输入时程的方式粗略调整人工波；之后，寻找绝对偏差最大的频率控制点，使用模拟退火算法寻找此频率处，各个阻尼比的标准输入时程的最优权重组合，以同时叠加多个阻尼比的标准输入时程，从而控制时域调整过程中叠加时程对响应谱的影响，进而提高此频率处所有阻尼比反应谱的拟合精度；寻找相对偏差最大的频率控制点，以同样的方式循环调整；最后，在后处理模块调整人工波，使其满足规范要求，并进一步提高拟合精度。本发明专利技术将由频率和阻尼比构成的多阻尼人工波拟合问题降维处理，简化为只需关注频率的一维问题，可生成高精度拟合多阻尼核电设计谱的人工地震波。

A multi damping artificial wave fitting method for constructing multi parametric time domain adjustment curve

【技术实现步骤摘要】
一种构造多参数化时域调整曲线的多阻尼人工波拟合方法
本专利技术属于核电设备工程抗震安全分析与评价
，涉及一种构造多参数化时域调整曲线的多阻尼人工波拟合方法，具体是指一种基于最优原理构造多参数化时域调整曲线的多阻尼人工波拟合方法。
技术介绍
确保核电站抗震安全是确保核电厂建设顺利实施和安全运行的前提与保障。使用数值计算的手段，分析工程结构的地震响应是有效评估其安全性的必要途径，其中选择具有代表性的地震动进行动力分析尤为重要。地震动代表性的一个关键因素是其与设计反应谱的一致性。阻尼比形式与结构形式、材料相关，因而核电设备阻尼比众多、复杂多变。进行多阻尼核电设备抗震分析时，主要有两种方法：1)使用多条单阻尼人工波分别进行抗震分析，最后再综合评价；2)直接使用与多阻尼设计谱拟合的一组人工波进行抗震分析。研究证明，使用与多阻尼设计谱拟合的人工波进行动力分析得到的结果具有离散型低，可靠性高等优点，并能有效减少结构抗震分析的计算量。现有的核电抗震设计规范，如CSAN289.3(2010)、USNRCSPR3.7.1(2014)、本文档来自技高网...

【技术保护点】
1.一种构造多参数化时域调整曲线的多阻尼人工波拟合方法，其特征在于包含以下步骤：/n第一步，按如下公式构造所有频率控制点和阻尼比处的标准输入时程：/n

【技术特征摘要】
1.一种构造多参数化时域调整曲线的多阻尼人工波拟合方法，其特征在于包含以下步骤：
第一步，按如下公式构造所有频率控制点和阻尼比处的标准输入时程：



式中，h(ζi,ωj,t)为赵凤新(2010)构造的增量窄带时程；m(t)为截断函数；C(ζi,ωj)为幅值标准化系数；
赵凤新(2010)利用三角函数和傅里叶正逆变换构造具有窄带特性的校正时程，以在时域内调整人工地震波，其构造的增量窄带时程h(ζi,ωj,t)如下所示：



式中，t0为所有带宽中最小带宽所对应的周期，即t0＝2π/min(ωb1,ωb2,...,ωbM)；显然，函数h(ζi,ωj,t)的峰值时间为t0，且标准输入时程的峰值时间也为t0；
使用截断函数m(t)在一个完整周期处截断增量窄带时程h(ζi,ωj,t)，其表达式如下所示：



幅值标准化系数C(ζi,ωj)的作用为调整标准输入时程的幅值，从而使其响应时程的峰值为1；幅值标准化系数C(ζi,ωj)的计算方法如下：



式中，H(ζi,ωj,tmax(ζi,ωj))为增量窄带时程h(ζi,ωj,t)的响应时程H(ζi,ωj,t)在其峰值时刻tmax(ζi,ωj)处的值，即响应时程H(ζi,ωj,t)的峰值；
第二步，使用预处理模块粗略调整初始地震波；预处理模块的主要作用为当初始波的拟合误差过大时，使用每个阻尼比、每个频率控制点循环迭代调整的策略叠加标准输入时程；其具体叠加方式如下所示：



式中，为第k次迭代调整后的人工波；Sk(ζi,ωj)为第k次迭代调整中标准输入时程的幅值调整系数；为移动后的标准输入时程；
幅值调整系数ΔSk(ζi,ωj)的计算方法为：



式中，ST(ζi,ωj)为在阻尼比ζi、频率控制点ωj处的设计谱值；为在阻尼比ζi、频率控制点ωj处时程的计算反应谱值；为符号函数，其含义如下所示：




为在阻尼比ζi、频率控制点ωj处的响应时程，tmax(ζi,ωj)为响应时程的峰值时间，故为在阻尼比ζi、频率控制点ωj处的响应时程的峰值；
当时域调整曲线的峰值时间与原时程的峰值时间接近时，才能达到预期的调整效果；因此，在叠加标准输入时程时，需将其沿时间轴平移，从而调整其峰值时间；移动距离Δt由下式计算：
t＝tmax(ζi,ωj)-t0(8)
式中，t0为式(2.5)中的标准输入时程的峰值时间；故移动后的标准输入时程的峰值时间为tmax(ζi,ωj)；
第三步，使用模拟退火算法寻找最优权重：
(1)确定初始温度T0、可接受温度Tf、降温系数a、初始权重α0以及变量x的初始值；
(2)确定值域空间，即标准输入时程的权重系数αk(ζi,ωj)的取值范围；其中：
0＜αk(...

【专利技术属性】
技术研发人员：李建波，程峰，林皋，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21