CN111259568A - 一种构造多参数化时域调整曲线的多阻尼人工波拟合方法 - Google Patents
一种构造多参数化时域调整曲线的多阻尼人工波拟合方法 Download PDFInfo
- Publication number
- CN111259568A CN111259568A CN202010119636.3A CN202010119636A CN111259568A CN 111259568 A CN111259568 A CN 111259568A CN 202010119636 A CN202010119636 A CN 202010119636A CN 111259568 A CN111259568 A CN 111259568A
- Authority
- CN
- China
- Prior art keywords
- damping
- standard input
- time
- artificial
- frequency control
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于核电设备工程抗震安全分析与评价技术领域,涉及一种构造多参数化时域调整曲线的多阻尼人工波拟合方法。该窄带时程构造方法,其在时域和频域内均具有窄带特性,且避免了由于傅里叶变换导致的振铃效应,这种新的窄带时程被称为标准输入时程;使用每个频率控制点每个阻尼比单独叠加标准输入时程的方式粗略调整人工波;之后,寻找绝对偏差最大的频率控制点,使用模拟退火算法寻找此频率处,各个阻尼比的标准输入时程的最优权重组合,以同时叠加多个阻尼比的标准输入时程,从而控制时域调整过程中叠加时程对响应谱的影响,进而提高此频率处所有阻尼比反应谱的拟合精度;寻找相对偏差最大的频率控制点,以同样的方式循环调整;最后,在后处理模块调整人工波,使其满足规范要求,并进一步提高拟合精度。本发明将由频率和阻尼比构成的多阻尼人工波拟合问题降维处理,简化为只需关注频率的一维问题,可生成高精度拟合多阻尼核电设计谱的人工地震波。
Description
技术领域
本发明属于核电设备工程抗震安全分析与评价技术领域,涉及一种构造多参数化时域调整曲线的多阻尼人工波拟合方法,具体是指一种基于最优原理构造多参数化时域调整曲线的多阻尼人工波拟合方法。
背景技术
确保核电站抗震安全是确保核电厂建设顺利实施和安全运行的前提与保障。使用数值计算的手段,分析工程结构的地震响应是有效评估其安全性的必要途径,其中选择具有代表性的地震动进行动力分析尤为重要。地震动代表性的一个关键因素是其与设计反应谱的一致性。阻尼比形式与结构形式、材料相关,因而核电设备阻尼比众多、复杂多变。进行多阻尼核电设备抗震分析时,主要有两种方法:
1)使用多条单阻尼人工波分别进行抗震分析,最后再综合评价;
2)直接使用与多阻尼设计谱拟合的一组人工波进行抗震分析。研究证明,使用与多阻尼设计谱拟合的人工波进行动力分析得到的结果具有离散型低,可靠性高等优点,并能有效减少结构抗震分析的计算量。
现有的核电抗震设计规范,如CSA N289.3(2010)、USNRC SPR 3.7.1(2014)、ASCE4-16(2016)等,均允许使用拟合多阻尼设计谱的人工波进行抗震分析并对其提出了具体的要求,如USNRC SRP 3.7.1中规定人工波计算反应谱值低于设计谱值的频率控制点不得超过五个。
目前多阻尼人工波拟合的方法主要为使用校正时程,一次调整一个频率一个阻尼比处的峰值响应,以消除求解的复杂性,避免由于强制同时拟合多个点所造成的奇异现象。然而,这种迭代技巧需要在所有频率控制点和阻尼比处循环迭代调整,且无法控制叠加时程对峰值响应之外的其余响应的影响,因而计算效率和拟合精度随频率控制点个数和阻尼比的增多而显著降低。
发明内容
本发明是基于最优原理构造多参数化时域调整曲线,将人工波时域调整法和模拟退火算法相结合的混合方法。该混合方法将由频率和阻尼比构成的多阻尼人工波拟合问题降维处理,简化为只需关注频率的一维问题。本发明使用模拟退火算法寻找某一特定控制频率处,各个阻尼比的时域调整曲线的最优权重组合,以同时叠加多个阻尼比的时域调整曲线,从而控制时域调整过程中叠加时程对响应谱的影响,进而生成高精度拟合多阻尼设计谱的人工地震波。
本发明包括如下步骤:
第一步,按如下公式构造所有频率控制点和阻尼比处的标准输入时程:
式中,h(ζi,ωj,t)为赵凤新(2010)构造的增量窄带时程;m(t)为截断函数;C(ζi,ωj)为幅值标准化系数。
赵凤新(2010)利用三角函数和傅里叶正逆变换构造具有窄带特性的校正时程,以在时域内调整人工地震波,其构造的增量窄带时程h(ζi,ωj,t)如下所示:
使用截断函数m(t)在一个完整周期处截断增量窄带时程h(ζi,ωj,t),其表达式如下所示:
式中,H(ζi,ωj,tmax(ζi,ωj))为增量窄带时程h(ζi,ωj,t)的响应时程H(ζi,ωj,t)在其峰值时刻tmax(ζi,ωj)处的值,即响应时程H(ζi,ωj,t)的峰值。
第二步,使用预处理模块粗略调整初始地震波。预处理模块的主要作用为当初始波的拟合误差过大时,使用每个阻尼比、每个频率控制点循环迭代调整的策略叠加标准输入时程。
其具体叠加方式如下所示:
幅值调整系数ΔSk(ζi,ωj)的计算方法为:
当时域调整曲线的峰值时间与原时程的峰值时间接近时,才能达到预期的调整效果。因此,在叠加标准输入时程时,需将其沿时间轴平移,从而调整其峰值时间。移动距离Δt由下式计算:
Δt=tmax(ζi,ωj)-t0 (8)
第三步,使用模拟退火算法寻找最优权重:
(1)确定初始温度T0、可接受温度Tf、降温系数a、初始权重α0以及变量x的初始值;
0<αk(ζi,ωj)≤1 (9)
其中,k表示此标准输入时程第k个被接受的权重系数。
(4)由如下概率表达式判断是否接受α*作为一个新解αk+1:
式中,ηj(αk)表示以权重αk叠加标准输入时程后的人工波,在频率控制点ωj处的所有阻尼比的计算反应谱与设计谱的误差的平方和。
(5)进行热传递,即若接受α*作为一个新解αk+1,则温度发生如下变化:
Tk+1=aTk (12)
式中,a为降温系数。
(6)判断温度Tk是否小于可接受温度Tf,若Tk>Tf,则由二分法生成新的变量x,并重复上述过程;若Tk≤Tf,则降温过程结束,此时的权重系数αk(ζi,ωj)为模拟退火算法得到的最优权重,记为αopt(ζi,ωj)。
第四步,使用核心模块中最大绝对偏差调整同时叠加多个阻尼比的标准输入时程。绝对偏差的定义如下式所示:
计算输入的人工波在全部阻尼比和控制频率处的绝对偏差,寻找最大绝对偏差所对应的控制频率。然后,依据步骤三所示的模拟退火算法具体步骤,寻找此频率控制点处所有阻尼比的标准输入时程的最优权重组合,按照下式叠加:
第五步,使用核心模块中最大相对偏差调整同时叠加多个阻尼比的标准输入时程。相对偏差,即拟合误差,其计算方法如下式计算:
计算输入的人工波在全部阻尼比和控制频率处的相对偏差,寻找最大相对偏差所对应的控制频率。然后,依据步骤三所示的模拟退火算法具体步骤,寻找此频率控制点处所有阻尼比的标准输入时程的最优权重组合,按照式(14)同时叠加此频率控制点处的所有阻尼比的标准输入时程,以进一步提高人工波对多阻尼设计谱的拟合精度。
第六步,在后处理模块,将美国核电规范委员会的标准评定大纲(USNRC SRP3.7.1)中给出的人工地震动的计算反应谱包络设计谱的条件:1)计算反应谱值低于设计谱值的频率控制点不超过五个;2)计算反应谱值最小不低于相应频率控制点的设计谱值的90%;以及作为地震动输入的同一组人工地震波的相关系数不得大于0.16,归纳为如下多目标优化模型:
Find:γi,j
Subject to:
ρ≤0.16
之后,按下式调整人工波:
重复上述过程,直至满足多目标优化模型所示的约束条件并且人工波的拟合精度不再提高后,调整过程结束,输出最终结果。
本发明将模拟退火算法和传统的人工波时域调整方法相结合,提出一种基于最优原理构造多参数化时域调整曲线的多阻尼人工波拟合方法。该方法不仅可以将多阻尼人工波拟合问题降维处理,还能有效控制叠加时程对除校正反应谱值之外的其余反应谱值的影响,配合使用在窄带时程基础上发展得到的标准输入时程,可得到满足约束条件的高精度拟合多阻尼目标谱的人工波。
附图说明
图1是本方法生成拟合多阻尼核电设计谱人工地震波的流程图。
图2a是实施例中5%阻尼比15Hz处的标准输入时程。
图2b是实施例中5%阻尼比15Hz处标准输入时程的傅里叶幅值谱。
图2c是实施例中5%阻尼比15Hz处标准输入时程的反应谱。
图3是实施例中作为初始波的归一化的Elcentro波。
图4是实施例中归一化的Elcentro波计算反应谱与RG1.60设计谱对比图。
图5是实施例中调整后的人工波加速度时程
图6是实施例中调整后人工波加速度时程的计算反应谱与RG1.60设计谱对比图
具体实施方式
下面结合具体方案和附图,详细叙述本发明的具体实施例。
本实施例方法以归一化的Elcentro波做为初始波,拟合0.5%、2%、5%、7%和10%五个阻尼比的RG1.60设计谱为例进行说明。包括以下步骤:
(1)按式(1)到(4)所示构造所有频率控制点和阻尼比处的标准输入时程,形成标准输入时程库,便于在后续调整中调用标准输入时程。5%阻尼比15Hz处的标准输入时程及其傅里叶幅值谱和反应谱如图2a、2b和2c所示;
(2)归一化的Elcentro波和其计算反应谱如图3和4所示,计算反应谱与设计谱的最大相对偏误差εmax=80.2%,并且各频率控制点的误差几乎全为负值,故使用预处理模块粗略调整初始波,缩小拟合误差。当最大拟合误差小于30%时,预处理模块结束,进入核心模块的最大绝对偏差调整;
(3)按照式(13)所示公式计算人工地震波在五个阻尼比所有频率控制点的绝对偏差,找出最大绝对偏差所对应的频率控制点;
(4)使用模拟退火算法,计算最大绝对偏差所对应的控制频率处五个阻尼比标准输入时程的最优权重组合,其中初始温度T0为1,可接受温度Tf为0.6,降温系数a为0.99,初始权重α0为0.25,变量x的初始值为0.6。并依据此权重组合同时叠加五个标准输入时程,以同时提高此频率控制点处所有阻尼比计算反应谱与设计谱的拟合精度;
(5)重复步骤(3)和(4)40次后,进入核心模块的最大相对偏差调整;
(6)按照式(15)所示公式计算人工地震波在五个阻尼比所有频率控制点的相对偏差,找出最大相对偏差所对应的频率控制点;
(7)使用模拟退火算法,计算最大相对偏差所对应的控制频率处五个阻尼比标准输入时程的最优权重组合,其中初始温度T0为1,可接受温度Tf为0.6,降温系数a为0.99,初始权重α0为0.25,变量x的初始值为0.6。并依据此权重组合同时叠加五个标准输入时程,以同时提高此频率控制点处所有阻尼比计算反应谱与设计谱的拟合精度;
(8)重复步骤(6)和(7)直至拟合精度不再提高,进入后处理模块;
(9)后处理模块,采用每个频率点每个阻尼比循环迭代调整的方式,调整人工波,使其满足规范USNRC SRP 3.7.1中的约束条件,并进一步提高其拟合精度。若某个频率点的调整不能提高拟合精度,则通过γi,j=0跳过此次调整。当人工波的拟合精度无法提高后,输出最终结果。调整后的人工波如图5所示,其计算反应谱与RG1.60设计谱的最大拟合误差εmax=4.98%,计算反应谱如图6所示。
由上可知,实施例方法合成的拟合多阻尼设计谱人工波在满足设计谱安全性、保守性要求的同时,与通常人工波拟合过程获得的人工波相比,其与目标谱的拟合精度更高。
Claims (1)
1.一种构造多参数化时域调整曲线的多阻尼人工波拟合方法,其特征在于包含以下步骤:
第一步,按如下公式构造所有频率控制点和阻尼比处的标准输入时程:
式中,h(ζi,ωj,t)为赵凤新(2010)构造的增量窄带时程;m(t)为截断函数;C(ζi,ωj)为幅值标准化系数;
赵凤新(2010)利用三角函数和傅里叶正逆变换构造具有窄带特性的校正时程,以在时域内调整人工地震波,其构造的增量窄带时程h(ζi,ωj,t)如下所示:
使用截断函数m(t)在一个完整周期处截断增量窄带时程h(ζi,ωj,t),其表达式如下所示:
式中,H(ζi,ωj,tmax(ζi,ωj))为增量窄带时程h(ζi,ωj,t)的响应时程H(ζi,ωj,t)在其峰值时刻tmax(ζi,ωj)处的值,即响应时程H(ζi,ωj,t)的峰值;
第二步,使用预处理模块粗略调整初始地震波;预处理模块的主要作用为当初始波的拟合误差过大时,使用每个阻尼比、每个频率控制点循环迭代调整的策略叠加标准输入时程;其具体叠加方式如下所示:
幅值调整系数ΔSk(ζi,ωj)的计算方法为:
当时域调整曲线的峰值时间与原时程的峰值时间接近时,才能达到预期的调整效果;因此,在叠加标准输入时程时,需将其沿时间轴平移,从而调整其峰值时间;移动距离Δt由下式计算:
t=tmax(ζi,ωj)-t0 (8)
第三步,使用模拟退火算法寻找最优权重:
(1)确定初始温度T0、可接受温度Tf、降温系数a、初始权重α0以及变量x的初始值;
0<αk(ζi,ωj)≤1 (9)
其中,k表示此标准输入时程第k个被接受的权重系数;
(4)由如下概率表达式判断是否接受α*作为一个新解αk+1:
式中,ηj(αk)表示以权重αk叠加标准输入时程后的人工波,在频率控制点ωj处的所有阻尼比的计算反应谱与设计谱的误差的平方和;
(5)进行热传递,即若接受α*作为一个新解αk+1,则温度发生如下变化:
Tk+1=aTk (12)
式中,a为降温系数;
(6)判断温度Tk是否小于可接受温度Tf,若Tk>Tf,则由二分法生成新的变量x,并重复上述过程;若Tk≤Tf,则降温过程结束,此时的权重系数αk(ζi,ωj)为模拟退火算法得到的最优权重,记为αopt(ζi,ωj);
第四步,使用核心模块中最大绝对偏差调整同时叠加多个阻尼比的标准输入时程;绝对偏差的定义如下式所示:
计算输入的人工波在全部阻尼比和控制频率处的绝对偏差,寻找最大绝对偏差所对应的控制频率;然后,依据步骤三所示的模拟退火算法具体步骤,寻找此频率控制点处所有阻尼比的标准输入时程的最优权重组合,按照下式叠加:
第五步,使用核心模块中最大相对偏差调整同时叠加多个阻尼比的标准输入时程;相对偏差,即拟合误差,其计算方法如下式计算:
计算输入的人工波在全部阻尼比和控制频率处的相对偏差,寻找最大相对偏差所对应的控制频率;然后,依据步骤三所示的模拟退火算法具体步骤,寻找此频率控制点处所有阻尼比的标准输入时程的最优权重组合,按照式(14)同时叠加此频率控制点处的所有阻尼比的标准输入时程,以进一步提高人工波对多阻尼设计谱的拟合精度;
第六步,在后处理模块,将美国核电规范委员会的标准评定大纲(USNRC SRP 3.7.1)中给出的人工地震动的计算反应谱包络设计谱的条件:
1)计算反应谱值低于设计谱值的频率控制点不超过五个;
2)计算反应谱值最小不低于相应频率控制点的设计谱值的90%;以及作为地震动输入的同一组人工地震波的相关系数不得大于0.16,归纳为如下多目标优化模型:
Find:γi,j
Subject to:
ρ≤0.16
之后,按下式调整人工波:
重复上述过程,直至满足多目标优化模型所示的约束条件并且人工波的拟合精度不再提高后,调整过程结束,输出最终结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010119636.3A CN111259568B (zh) | 2020-02-26 | 2020-02-26 | 一种构造多参数化时域调整曲线的多阻尼人工波拟合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010119636.3A CN111259568B (zh) | 2020-02-26 | 2020-02-26 | 一种构造多参数化时域调整曲线的多阻尼人工波拟合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111259568A true CN111259568A (zh) | 2020-06-09 |
CN111259568B CN111259568B (zh) | 2022-09-20 |
Family
ID=70954623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010119636.3A Active CN111259568B (zh) | 2020-02-26 | 2020-02-26 | 一种构造多参数化时域调整曲线的多阻尼人工波拟合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111259568B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112462424A (zh) * | 2020-11-13 | 2021-03-09 | 大连理工大学 | 一种考虑真实波相位特征的人工波生成方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236256A (zh) * | 2008-01-25 | 2008-08-06 | 北京工业大学 | 多阻尼比目标反应谱兼容的人工地震波合成方法 |
-
2020
- 2020-02-26 CN CN202010119636.3A patent/CN111259568B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236256A (zh) * | 2008-01-25 | 2008-08-06 | 北京工业大学 | 多阻尼比目标反应谱兼容的人工地震波合成方法 |
Non-Patent Citations (3)
Title |
---|
何佳等: "用于窄带叠加人工时程模拟的一种窄带构造算法", 《核动力工程》 * |
张伟明等: "考虑增量时程贡献趋向和误差排序的多阻尼目标反应谱拟合", 《地震研究》 * |
赵凤新等: "多阻尼反应谱拟合的时域叠加法", 《核动力工程》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112462424A (zh) * | 2020-11-13 | 2021-03-09 | 大连理工大学 | 一种考虑真实波相位特征的人工波生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111259568B (zh) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108462165B (zh) | 一种新能源接入电力系统的负荷特性评估方法 | |
CN105181121B (zh) | 采用加权迭代等效源法的高精度近场声全息方法 | |
CN109299430A (zh) | 基于两阶段分解与极限学习机的短期风速预测方法 | |
CN111259568B (zh) | 一种构造多参数化时域调整曲线的多阻尼人工波拟合方法 | |
CN110618453A (zh) | 一种基于改进阻尼最小二乘法的波阻抗反演方法 | |
CN103972908A (zh) | 一种基于自适应混沌粒子群算法的多目标无功优化方法 | |
CN106842306A (zh) | 一种全局优化的交错网格有限差分正演模拟方法和装置 | |
CN110598923A (zh) | 基于支持向量回归机优化与误差修正的空调负荷预测方法 | |
CN109100935A (zh) | 大时滞系统的阻尼智慧pi控制方法 | |
CN116933922A (zh) | 基于woa-cnn-lstm光伏发电功率预测的方法 | |
CN112821420B (zh) | 一种基于XGBoost的ASFR模型中动态阻尼因子、多维频率指标的预测方法及系统 | |
CN108566256A (zh) | 一种频谱地图的构建方法 | |
Mihálka et al. | Application of the Cauchy integral formula as a tool of analytic continuation for the resummation of divergent perturbation series | |
CN116300422A (zh) | 一种水电机组控制优化方法及终端 | |
Abbas et al. | State-space realization of LPV input-output models: practical methods for the user | |
CN105096328A (zh) | 一种基于ios-pso的摄像机自标定方法 | |
Guozhen et al. | Tidal current short-term prediction based on support vector regression | |
CN111522226B (zh) | 针对伺服转台的多目标优化高型pid最优控制器设计方法 | |
CN112462424A (zh) | 一种考虑真实波相位特征的人工波生成方法 | |
Bensenouci et al. | Performance analysis and comparative study of LMI-based iterative PID load-frequency controllers of a singlearea power system | |
Maknunah et al. | Estimation of solar radiation per month on horizontal surface using adaptive neuro-fuzzy inference system (case study in Surabaya) | |
CN107145667B (zh) | 一种波前结构函数的快速计算方法 | |
CN117610317B (zh) | 基于深度学习的多bit超表面相位排布优化方法 | |
CN106126818B (zh) | 一种自适应的在线产生子群参数的方法 | |
CN112649869A (zh) | 基于ga-wnn的储层特征参数预测方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |