CN109613611A - 用于结构抗震时程分析的输入地震波的确定方法及系统 - Google Patents
用于结构抗震时程分析的输入地震波的确定方法及系统 Download PDFInfo
- Publication number
- CN109613611A CN109613611A CN201910066254.6A CN201910066254A CN109613611A CN 109613611 A CN109613611 A CN 109613611A CN 201910066254 A CN201910066254 A CN 201910066254A CN 109613611 A CN109613611 A CN 109613611A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- newmark
- tri
- response spectrum
- quasi
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 47
- 230000006641 stabilisation Effects 0.000 title claims abstract description 36
- 238000011105 stabilization Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 208
- 238000012216 screening Methods 0.000 claims abstract description 14
- 230000003321 amplification Effects 0.000 claims description 41
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 41
- 230000001133 acceleration Effects 0.000 claims description 34
- 238000006073 displacement reaction Methods 0.000 claims description 25
- 230000000737 periodic effect Effects 0.000 claims description 17
- 230000003595 spectral effect Effects 0.000 claims description 16
- 238000012502 risk assessment Methods 0.000 claims description 11
- 238000009499 grossing Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 2
- UOORRWUZONOOLO-OWOJBTEDSA-N (E)-1,3-dichloropropene Chemical compound ClC\C=C\Cl UOORRWUZONOOLO-OWOJBTEDSA-N 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 208000012788 shakes Diseases 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种用于结构抗震时程分析的输入地震波的确定方法及系统。包括:获取多条备选地震波数据;获取待分析场地的Newmark三联谱;将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式;将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数;筛选设定缩放系数范围内的拟速度反应谱,按照匹配误差的大小进行排序,将匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。本发明可以提高结构抗震时程分析的准确性,并可有效降低结构时程分析反应结果的离散性。
Description
技术领域
本发明涉及结构抗震分析领域,特别是涉及一种用于结构抗震时程分析的输入地震波的确定方法及系统。
背景技术
时程分析法在抗震设计中称为“动态设计”,是指由结构基本运动方程输入地面加速度记录进行积分求解,以求得整个时间历程的地震反应的方法。此法输入对应于工程场地的若干条地震加速度记录或人工加速度时程曲线,通过积分运算求得在地面加速度随时间变化期间结构的内力和变形状态随时间变化的全过程,并以此进行结构构件的截面抗震承载力验算和变形验算。在确定输入地震波的过程中,需要采用目标谱进行选择。
目前国内外选波研究采用的目标谱有规范设计谱、一致概率谱(UHS)、条件均值谱(CMS或CS)以及其它一些类型的目标谱等。采用目标谱不同,选波过程及结果也有所差异。规范设计谱(加速度谱或影响系数谱或放大系数谱)是目前各国学者和工程师最常选用的目标谱,但统计依据简单且地震地质背景信息不清晰,使其很难保证结构反应预估的准确性。一致概率反应谱(UHS)是基于概率地震危险性分析(PSHA)建立的设计谱,能够考虑地震地质背景,广泛应用于美国核工业工程的设计规范。UHS本质的局限在于认为各个周期点具有相同的超越概率,而这并不符合一次地震事件的概率分布特征,由于其过于保守,很多国家规范并不采用。条件均值谱(CMS)由Baker(2011)采用Sa(T1)与谱形系数Epsilon两个向量地震动强度指标,考虑了不同周期点地震反应的相关性并引入条件分布,构建了CMS不仅避免了UHS的局限性,谱形系数Epsilon的采用也使以CMS为目标谱的选波方法在多自由度体系非线性反应分析中显现出明显的优势。此外,Baker基于CMS还提出了能够考虑反应谱分布的条件分布谱(CS)。条件均值谱的建立受衰减关系、震级M、距离R及各周期谱值间相关性(即相关系数ρ)等因素的影响。我国规范在地震动参数区划图制定中已建立了衰减关系,但与PSHA相关的地震危险性分析信息尚未很好地提供(无法获得不同震源地震风险贡献率),将CMS与我国规范相结合用于选波工作,仍有多处细节需考虑。
上述目标谱多为加速度反应谱,主要是反映地震峰值加速度PGA或加速度反应谱Sa的地震衰减或统计特征,与中长周期结构反应的相关性不够密切,仍具有一定的局限性。虽然也有学者采用非线性位移目标谱分析长周期结构,但由于未与PSHA相关联,且参考依据不统一,应用及发展均受到限制。
发明内容
本发明的目的是提供一种用于结构抗震时程分析的输入地震波的确定方法及系统,可应用于高层、大跨桥梁等重要及复杂工程结构在进行弹性或弹塑性时程分析时进行输入地震波的确定,以提高结构抗震时程分析的准确性,并可有效降低结构时程分析反应结果的离散性。
为实现上述目的,本发明提供了如下方案:
一种用于结构抗震时程分析的输入地震波的确定方法,包括:
获取多条备选地震波数据;
获取待分析场地的Newmark三联谱;
将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式;
将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数;
筛选设定缩放系数范围内的拟速度反应谱,得到第一拟速度反应谱集合;所述第一速度反应谱集合中的所有拟速度反应谱的缩放系数均属于所述设定缩放系数范围内;
将所述第一拟速度反应谱集合中的元素按照匹配误差的大小进行排序,得到第一序列;
将所述第一序列中匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。
可选的,所述获取多条备选地震波数据,具体包括:
获取地震信息数据库;
根据待分析场地的参数,从所述地震信息数据库中筛选多条备选地震波数据。
可选的,所述获取待分析场地的Newmark三联谱,具体包括:
获取待分析场地的地震危险性分析结果;
确定设定超越概率的地震风险水平对应的加速度峰值、速度峰值和位移峰值;
获取所述加速度峰值对应的加速度反应谱的放大系数、所述速度峰值对应的速度反应谱的放大系数和所述位移峰值对应的位移反应谱的放大系数;
根据所述加速度反应谱的放大系数、所述速度反应谱的放大系数和所述位移反应谱的放大系数,构建平滑的拟速度反应谱形式的Newmark三联谱,得到所述待分析场地的Newmark三联谱。
可选的,所述将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式,具体包括:
确定Newmark三联谱的各拐点周期;所述Newmark三联谱的拐点包括Ta、Tb、Tc、Td和Te,其中,Ta=1/33s,Tb=1/8s,Te=10s,Tf=33s,Tc和Td为变量,Tc大于Tb且小于Td,Td小于Te;
根据所述Newmark三联谱的拐点绘制Newmark三联谱曲线,每条备选地震波的拟速度反应谱在对应的Newmark三联谱曲线上方的区域和下方的区域的面积相等。
可选的,所述将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数,具体包括:
采用公式获取每条拟速度反应谱与所述目标谱的匹配误差SSE;
采用公式获取每条拟速度反应谱与所述目标谱的缩放系数SF;
其中,PSv(Ti)为拟速度反应谱在Ti周期点对应的拟速度谱值,为目标谱在Ti周期点对应的拟速度谱值;n为整个匹配周期段内间隔0.05s确定的周期点总数。
本发明还提供一种用于结构抗震时程分析的输入地震波的确定系统,包括:
备选地震波数据获取模块,用于获取多条备选地震波数据;
Newmark三联谱获取模块,用于获取待分析场地的Newmark三联谱;
平滑化模块,用于将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式;
匹配参数获取模块,用于将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数;
筛选模块,用于筛选设定缩放系数范围内的拟速度反应谱,得到第一拟速度反应谱集合;所述第一速度反应谱集合中的所有拟速度反应谱的缩放系数均属于所述设定缩放系数范围内;
排序模块,用于将所述第一拟速度反应谱集合中的元素按照匹配误差的大小进行排序,得到第一序列;
输入地震波确定模块,用于将所述第一序列中匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。
可选的,所述备选地震波数据获取模块具体包括:
地震信息数据库获取单元,用于获取地震信息数据库;
筛选单元,用于根据待分析场地的参数,从所述地震信息数据库中筛选多条备选地震波数据。
可选的,所述Newmark三联谱获取模块具体包括:
地震危险性分析结果获取单元,用于获取待分析场地的地震危险性分析结果;
地震风险水平参数确定单元,用于确定设定超越概率的地震风险水平对应的加速度峰值、速度峰值和位移峰值;
放大系数获取单元,用于获取所述加速度峰值对应的加速度反应谱的放大系数、所述速度峰值对应的速度反应谱的放大系数和所述位移峰值对应的位移反应谱的放大系数;
Newmark三联谱构建单元,用于根据所述加速度反应谱的放大系数、所述速度反应谱的放大系数和所述位移反应谱的放大系数,构建平滑的拟速度反应谱形式的Newmark三联谱,得到所述待分析场地的Newmark三联谱。
可选的,所述平滑化模块具体包括:
各拐点周期确定单元,用于确定Newmark三联谱的各拐点周期;所述Newmark三联谱的拐点包括Ta、Tb、Tc、Td和Te,其中,Ta=1/33s,Tb=1/8s,Te=10s,Tf=33s,Tc和Td为变量,Tc大于Tb且小于Td,Td小于Te;
Newmark三联谱曲线绘制单元,用于根据所述Newmark三联谱的拐点绘制Newmark三联谱曲线,每条备选地震波的拟速度反应谱在对应的Newmark三联谱曲线上方的区域和下方的区域的面积相等。
可选的,所述匹配参数获取模块具体包括:
匹配误差获取单元,用于采用公式获取每条拟速度反应谱与所述目标谱的匹配误差SSE;
缩放系数获取单元,用于采用公式获取每条拟速度反应谱与所述目标谱的缩放系数SF;
其中,PSv(Ti)为拟速度反应谱在Ti周期点对应的拟速度谱值,为目标谱在Ti周期点对应的拟速度谱值;n为整个匹配周期段内间隔0.05s确定的周期点总数。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明采用Newmark三联谱作为目标谱,Newmark三联谱的短、中和长周期段分别基于地震动强度指标PGA、PGV和PGD构建,与不同周期结构反应具有“天然”的良好相关性,将其作为目标谱,可为结构时程分析选波提供新的解决途径。传统的选波过程所用目标谱多为加速度反应谱,主要是反映PGA或Sa的地震衰减或统计特征,与中长周期结构反应的相关性不够密切,仍具有一定的局限性。本发明与传统的采用加速度反应谱为目标谱的选波方法相比,可有效降低结构反应离散性,当结构高度周期较长、结构非线性程度较大时,这一优势则更为明显。时程反应结果较低的离散性对工程师意味着选择和处理上的方便性,也一定程度上提升了抗震设计结果的可信性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明用于结构抗震时程分析的输入地震波的确定方法的流程示意图;
图2为本发明用于结构抗震时程分析的输入地震波的确定方法中将拟速度反应谱平滑化成Newmark三联谱的示意图;
图3为本发明用于结构抗震时程分析的输入地震波的确定方法中匹配误差计算示意图;
图4为本发明用于结构抗震时程分析的输入地震波的确定系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明用于结构抗震时程分析的输入地震波的确定方法的流程示意图。如图1所示,所述确定方法包括以下步骤:
步骤100:获取多条备选地震波数据。备选地震波数据为从地震信息数据库中,根据待分析场地的参数筛选的多条备选地震波数据。例如,可从美国的太平洋地震工程中心(PEER)的强震记录数据库(NGA)、日本的强震观测计划(K-NET和KiK-NET)等公开数据库获得,这些数据库提供免费自由下载服务。待分析场地的参数可以为一定的地震信息,如限定震级等级大于6、限定震中距或断层距小于60km等,也可以为限定场地条件等,根据待分析场地的参数设定备选地震波的初选条件,即可得到筛选后的多条备选地震波数据。大多数公开数据库都支持上述操作。
步骤200:获取待分析场地的Newmark三联谱。具体过程如下:
(1)获取待分析场地的地震危险性分析结果。
(2)根据地震危险性分析结果,确定设定超越概率的地震风险水平对应的加速度峰值PGA、速度峰值PGV和位移峰值PGD;此处,设定超越概率根据实际情况而定,例如设计规范规定为50年超越概率10%,此处设定超越概率即为50年超越概率10%。
(3)获取所述加速度峰值对应的加速度反应谱的放大系数、所述速度峰值对应的速度反应谱的放大系数和所述位移峰值对应的位移反应谱的放大系数。
(4)根据所述加速度反应谱的放大系数、所述速度反应谱的放大系数和所述位移反应谱的放大系数,构建平滑的拟速度反应谱形式的Newmark三联谱,得到所述待分析场地的Newmark三联谱,即为目标谱。
步骤300:将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式。首先,确定Newmark三联谱的各拐点周期。如图2所示,图2为本发明用于结构抗震时程分析的输入地震波的确定方法中将拟速度反应谱平滑化成Newmark三联谱的示意图。Newmark三联谱中各拐点周期的确定,可取如下值Ta=1/33s,Tb=1/8s,Te=10s,Tf=33s,Tc大于Tb且小于Td,Td小于Te;Tc和Td对于各条备选地震波是变量,它们确定的原则是:保证备选波的拟速度反应谱(PSv谱)在平滑化的三联谱曲线上方的区域和下方的区域的面积相等,即图2中折线上方阴影部分和下方阴影部分的面积相等。然后,根据所述Newmark三联谱的拐点绘制Newmark三联谱曲线,每条备选地震波的拟速度反应谱在对应的Newmark三联谱曲线上方的区域和下方的区域的面积相等。
步骤400:将待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与目标谱的匹配误差和缩放系数。
具体的,采用公式获取每条拟速度反应谱与所述目标谱的匹配误差SSE。得到的结果如图3示意,图3为本发明用于结构抗震时程分析的输入地震波的确定方法中匹配误差计算示意图;
采用公式获取每条拟速度反应谱与所述目标谱的缩放系数SF;其中,PSv(Ti)为拟速度反应谱在Ti周期点对应的拟速度谱值,为目标谱在Ti周期点对应的拟速度谱值;n为整个匹配周期段内间隔0.05s确定的周期点总数。匹配周期段的下限为Tm,其中m为保证振型参与质量不小于总质量的90%需考虑的振型数,匹配周期段上限的取值可参考ASCE7-10或参考欧洲规范。
步骤500:筛选设定缩放系数范围内的拟速度反应谱,得到第一拟速度反应谱集合。第一速度反应谱集合中的所有拟速度反应谱的缩放系数均属于设定缩放系数范围内。
步骤600:将第一拟速度反应谱集合中的元素按照匹配误差的大小进行排序,得到第一序列。
步骤700:将第一序列中匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。
设定缩放系数范围和N的大小根据实际情况而定,例如,可依据匹配误差SSE较小且缩放系数SF接近1为原则将备选地震波对应的拟速度反应谱进行较简单的排序,如设定缩放系数范围为SF不大于3,不小于0.5,然后在设定缩放系数范围内,按照SSE从小到大排列,取一定数量地震波(例如,通常设计规范要求7~10条),作为输入输入地震波,进行结构时程分析。
本发明针对基于目标谱进行时程分析选波研究,提出了一种新的目标谱,即Newmark三联谱。因Newmark三联谱的短、中、长周期段(频率敏感范围)分别与地震动的PGA、PGV和PGD相关,其谱值为各个区段的放大系数乘以对应的地震动幅值参数,具有与不同周期结构反应“天然”的良好相关性,同时还可突破选波研究中大多基于加速度目标谱的局限。本方法可应用于高层、大跨桥梁等重要及复杂工程结构在进行弹性或弹塑性时程分析时所需输入地震波的选择及(线性)调整,能够保证结构时程反应估计的准确性,并可有效降低结构时程分析反应结果的离散性,且对于较长周期结构以及较高地震危险性水平条件下,这一优势更为凸显。
图4为本发明用于结构抗震时程分析的输入地震波的确定系统的结构示意图。如图4所示,所述确定系统包括:
备选地震波数据获取模块401,用于获取多条备选地震波数据;
Newmark三联谱获取模块402,用于获取待分析场地的Newmark三联谱;
平滑化模块403,用于将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式;
匹配参数获取模块404,用于将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数;
筛选模块405,用于筛选设定缩放系数范围内的拟速度反应谱,得到第一拟速度反应谱集合;所述第一速度反应谱集合中的所有拟速度反应谱的缩放系数均属于所述设定缩放系数范围内;
排序模块406,用于将所述第一拟速度反应谱集合中的元素按照匹配误差的大小进行排序,得到第一序列;
输入地震波确定模块407,用于将所述第一序列中匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。
所述备选地震波数据获取模块401具体包括:
地震信息数据库获取单元,用于获取地震信息数据库;
筛选单元,用于根据待分析场地的参数,从所述地震信息数据库中筛选多条备选地震波数据。
所述Newmark三联谱获取模块402具体包括:
地震危险性分析结果获取单元,用于获取待分析场地的地震危险性分析结果;
地震风险水平参数确定单元,用于确定设定超越概率的地震风险水平对应的加速度峰值、速度峰值和位移峰值;
放大系数获取单元,用于获取所述加速度峰值对应的加速度反应谱的放大系数、所述速度峰值对应的速度反应谱的放大系数和所述位移峰值对应的位移反应谱的放大系数;
Newmark三联谱构建单元,用于根据所述加速度反应谱的放大系数、所述速度反应谱的放大系数和所述位移反应谱的放大系数,构建平滑的拟速度反应谱形式的Newmark三联谱,得到所述待分析场地的Newmark三联谱。
所述平滑化模块403具体包括:
各拐点周期确定单元,用于确定Newmark三联谱的各拐点周期;所述Newmark三联谱的拐点包括Ta、Tb、Tc、Td和Te,其中,Ta=1/33s,Tb=1/8s,Te=10s,Tf=33s,Tc和Td为变量,Tc大于Tb且小于Td,Td小于Te;
Newmark三联谱曲线绘制单元,用于根据所述Newmark三联谱的拐点绘制Newmark三联谱曲线,每条备选地震波的拟速度反应谱在对应的Newmark三联谱曲线上方的区域和下方的区域的面积相等。
所述匹配参数获取模块404具体包括:
匹配误差获取单元,用于采用公式获取每条拟速度反应谱与所述目标谱的匹配误差SSE;
缩放系数获取单元,用于采用公式获取每条拟速度反应谱与所述目标谱的缩放系数SF;
其中,PSv(Ti)为拟速度反应谱在Ti周期点对应的拟速度谱值,为目标谱在Ti周期点对应的拟速度谱值;n为整个匹配周期段内间隔0.05s确定的周期点总数。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种用于结构抗震时程分析的输入地震波的确定方法,其特征在于,包括:
获取多条备选地震波数据;
获取待分析场地的Newmark三联谱;
将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式;
将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数;
筛选设定缩放系数范围内的拟速度反应谱,得到第一拟速度反应谱集合;所述第一速度反应谱集合中的所有拟速度反应谱的缩放系数均属于所述设定缩放系数范围内;
将所述第一拟速度反应谱集合中的元素按照匹配误差的大小进行排序,得到第一序列;
将所述第一序列中匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。
2.根据权利要求1所述的用于结构抗震时程分析的输入地震波的确定方法,其特征在于,所述获取多条备选地震波数据,具体包括:
获取地震信息数据库;
根据待分析场地的参数,从所述地震信息数据库中筛选多条备选地震波数据。
3.根据权利要求1所述的用于结构抗震时程分析的输入地震波的确定方法,其特征在于,所述获取待分析场地的Newmark三联谱,具体包括:
获取待分析场地的地震危险性分析结果;
确定设定超越概率的地震风险水平对应的加速度峰值、速度峰值和位移峰值;
获取所述加速度峰值对应的加速度反应谱的放大系数、所述速度峰值对应的速度反应谱的放大系数和所述位移峰值对应的位移反应谱的放大系数;
根据所述加速度反应谱的放大系数、所述速度反应谱的放大系数和所述位移反应谱的放大系数,构建平滑的拟速度反应谱形式的Newmark三联谱,得到所述待分析场地的Newmark三联谱。
4.根据权利要求1所述的用于结构抗震时程分析的输入地震波的确定方法,其特征在于,所述将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式,具体包括:
确定Newmark三联谱的各拐点周期;所述Newmark三联谱的拐点包括Ta、Tb、Tc、Td和Te,其中,Ta=1/33s,Tb=1/8s,Te=10s,Tf=33s,Tc和Td为变量,Tc大于Tb且小于Td,Td小于Te;
根据所述Newmark三联谱的拐点绘制Newmark三联谱曲线,每条备选地震波的拟速度反应谱在对应的Newmark三联谱曲线上方的区域和下方的区域的面积相等。
5.根据权利要求1所述的用于结构抗震时程分析的输入地震波的确定方法,其特征在于,将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数,具体包括:
采用公式获取每条拟速度反应谱与所述目标谱的匹配误差SSE;
采用公式获取每条拟速度反应谱与所述目标谱的缩放系数SF;
其中,PSv(Ti)为拟速度反应谱在Ti周期点对应的拟速度谱值,为目标谱在Ti周期点对应的拟速度谱值;n为整个匹配周期段内间隔0.05s确定的周期点总数。
6.一种用于结构抗震时程分析的输入地震波的确定系统,其特征在于,包括:
备选地震波数据获取模块,用于获取多条备选地震波数据;
Newmark三联谱获取模块,用于获取待分析场地的Newmark三联谱;
平滑化模块,用于将每条备选地震波的拟速度反应谱平滑化成Newmark三联谱形式;
匹配参数获取模块,用于将所述待分析场地的Newmark三联谱作为目标谱,获取每条拟速度反应谱与所述目标谱的匹配误差和缩放系数;
筛选模块,用于筛选设定缩放系数范围内的拟速度反应谱,得到第一拟速度反应谱集合;所述第一速度反应谱集合中的所有拟速度反应谱的缩放系数均属于所述设定缩放系数范围内;
排序模块,用于将所述第一拟速度反应谱集合中的元素按照匹配误差的大小进行排序,得到第一序列;
输入地震波确定模块,用于将所述第一序列中匹配误差最小的N条拟速度反应谱对应的备选地震波,确定为结构抗震时程分析的输入地震波。
7.根据权利要求6所述的用于结构抗震时程分析的输入地震波的确定系统,其特征在于,所述备选地震波数据获取模块具体包括:
地震信息数据库获取单元,用于获取地震信息数据库;
筛选单元,用于根据待分析场地的参数,从所述地震信息数据库中筛选多条备选地震波数据。
8.根据权利要求6所述的用于结构抗震时程分析的输入地震波的确定系统,其特征在于,所述Newmark三联谱获取模块具体包括:
地震危险性分析结果获取单元,用于获取待分析场地的地震危险性分析结果;
地震风险水平参数确定单元,用于确定设定超越概率的地震风险水平对应的加速度峰值、速度峰值和位移峰值;
放大系数获取单元,用于获取所述加速度峰值对应的加速度反应谱的放大系数、所述速度峰值对应的速度反应谱的放大系数和所述位移峰值对应的位移反应谱的放大系数;
Newmark三联谱构建单元,用于根据所述加速度反应谱的放大系数、所述速度反应谱的放大系数和所述位移反应谱的放大系数,构建平滑的拟速度反应谱形式的Newmark三联谱,得到所述待分析场地的Newmark三联谱。
9.根据权利要求6所述的用于结构抗震时程分析的输入地震波的确定系统,其特征在于,所述平滑化模块具体包括:
各拐点周期确定单元,用于确定Newmark三联谱的各拐点周期;所述Newmark三联谱的拐点包括Ta、Tb、Tc、Td和Te,其中,Ta=1/33s,Tb=1/8s,Te=10s,Tf=33s,Tc和Td为变量,Tc大于Tb且小于Td,Td小于Te;
Newmark三联谱曲线绘制单元,用于根据所述Newmark三联谱的拐点绘制Newmark三联谱曲线,每条备选地震波的拟速度反应谱在对应的Newmark三联谱曲线上方的区域和下方的区域的面积相等。
10.根据权利要求6所述的用于结构抗震时程分析的输入地震波的确定系统,其特征在于,所述匹配参数获取模块具体包括:
匹配误差获取单元,用于采用公式获取每条拟速度反应谱与所述目标谱的匹配误差SSE;
缩放系数获取单元,用于采用公式获取每条拟速度反应谱与所述目标谱的缩放系数SF;
其中,PSv(Ti)为拟速度反应谱在Ti周期点对应的拟速度谱值,为目标谱在Ti周期点对应的拟速度谱值;n为整个匹配周期段内间隔0.05s确定的周期点总数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910066254.6A CN109613611B (zh) | 2019-01-24 | 2019-01-24 | 用于结构抗震时程分析的输入地震波的确定方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910066254.6A CN109613611B (zh) | 2019-01-24 | 2019-01-24 | 用于结构抗震时程分析的输入地震波的确定方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109613611A true CN109613611A (zh) | 2019-04-12 |
CN109613611B CN109613611B (zh) | 2020-07-03 |
Family
ID=66019271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910066254.6A Active CN109613611B (zh) | 2019-01-24 | 2019-01-24 | 用于结构抗震时程分析的输入地震波的确定方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109613611B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110134682A (zh) * | 2019-04-16 | 2019-08-16 | 同济大学 | 基于随机地震动数据库构建的数据交互方法及数据库装置 |
CN112613096A (zh) * | 2020-12-15 | 2021-04-06 | 应急管理部国家自然灾害防治研究院 | 一种面向强震前后不同阶段的地质灾害评价方法 |
CN112699438A (zh) * | 2020-12-25 | 2021-04-23 | 中国地震局工程力学研究所 | 基于破坏强度的全球地震动定量排序的输入地震动选取方法 |
CN112989475A (zh) * | 2021-04-06 | 2021-06-18 | 哈尔滨工业大学(威海) | 一种适用于中国抗震设计规范的位移设计谱方法 |
CN115267883A (zh) * | 2022-06-23 | 2022-11-01 | 华南理工大学 | 地震响应预测模型训练及预测方法、系统、设备及介质 |
CN116168115A (zh) * | 2023-04-20 | 2023-05-26 | 山东省地震工程研究院 | 一种基岩反应谱和超越概率曲线自动化绘图方法和系统 |
CN116776208A (zh) * | 2023-08-25 | 2023-09-19 | 华侨大学 | 地震波分类模型的训练方法、地震选波方法、设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090259405A1 (en) * | 2008-04-15 | 2009-10-15 | Battelle Energy Alliance, Llc | Methods, systems, and computer-readable media for generating seismic event time histories |
CN103364829A (zh) * | 2013-07-17 | 2013-10-23 | 大连海事大学 | 复杂结构抗震时程分析输入地震波的选择方法 |
CN104200128A (zh) * | 2014-09-26 | 2014-12-10 | 广西交通科学研究院 | 一种结构地震反应分析中地震波的选取及调整方法 |
CN106842319A (zh) * | 2017-01-09 | 2017-06-13 | 青岛理工大学 | 一种时程分析地震动备选库建立方法 |
-
2019
- 2019-01-24 CN CN201910066254.6A patent/CN109613611B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090259405A1 (en) * | 2008-04-15 | 2009-10-15 | Battelle Energy Alliance, Llc | Methods, systems, and computer-readable media for generating seismic event time histories |
CN103364829A (zh) * | 2013-07-17 | 2013-10-23 | 大连海事大学 | 复杂结构抗震时程分析输入地震波的选择方法 |
CN104200128A (zh) * | 2014-09-26 | 2014-12-10 | 广西交通科学研究院 | 一种结构地震反应分析中地震波的选取及调整方法 |
CN106842319A (zh) * | 2017-01-09 | 2017-06-13 | 青岛理工大学 | 一种时程分析地震动备选库建立方法 |
Non-Patent Citations (1)
Title |
---|
张锐 等: "时程分析考虑高阶振型影响的多频段地震波选择方法研究", 《工程力学》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110134682A (zh) * | 2019-04-16 | 2019-08-16 | 同济大学 | 基于随机地震动数据库构建的数据交互方法及数据库装置 |
CN112613096A (zh) * | 2020-12-15 | 2021-04-06 | 应急管理部国家自然灾害防治研究院 | 一种面向强震前后不同阶段的地质灾害评价方法 |
CN112613096B (zh) * | 2020-12-15 | 2024-02-23 | 应急管理部国家自然灾害防治研究院 | 一种面向强震前后不同阶段的地质灾害评价方法 |
CN112699438A (zh) * | 2020-12-25 | 2021-04-23 | 中国地震局工程力学研究所 | 基于破坏强度的全球地震动定量排序的输入地震动选取方法 |
CN112989475A (zh) * | 2021-04-06 | 2021-06-18 | 哈尔滨工业大学(威海) | 一种适用于中国抗震设计规范的位移设计谱方法 |
CN115267883A (zh) * | 2022-06-23 | 2022-11-01 | 华南理工大学 | 地震响应预测模型训练及预测方法、系统、设备及介质 |
CN115267883B (zh) * | 2022-06-23 | 2023-08-29 | 华南理工大学 | 地震响应预测模型训练及预测方法、系统、设备及介质 |
CN116168115A (zh) * | 2023-04-20 | 2023-05-26 | 山东省地震工程研究院 | 一种基岩反应谱和超越概率曲线自动化绘图方法和系统 |
CN116776208A (zh) * | 2023-08-25 | 2023-09-19 | 华侨大学 | 地震波分类模型的训练方法、地震选波方法、设备及介质 |
CN116776208B (zh) * | 2023-08-25 | 2023-10-20 | 华侨大学 | 地震波分类模型的训练方法、地震选波方法、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109613611B (zh) | 2020-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109613611A (zh) | 用于结构抗震时程分析的输入地震波的确定方法及系统 | |
US20230260599A1 (en) | Method of quantifying soil carbon | |
Juang et al. | Probabilistic framework for liquefaction potential by shear wave velocity | |
Li et al. | Bayesian updating of fragility functions using hybrid simulation | |
CN106404620A (zh) | 地统计插值与卫星遥感联合反演地面pm2.5的方法及系统 | |
CN106649954A (zh) | 一种基于扩展pbee2理论框架的地震易损性分析方法 | |
CN105974495A (zh) | 利用分类拟合法预判目标区域未来平均云量的方法 | |
Miano et al. | Considering structural modeling uncertainties using Bayesian cloud analysis | |
Pomeroy et al. | Putting prediction in ungauged basins into practice | |
CN111914457A (zh) | 一种输电线路塔基边坡稳定性判别方法、装置和存储介质 | |
CN110096805A (zh) | 一种有限观测数据下基于改进自助法的结构参数不确定性量化及传递方法 | |
Mu et al. | Seismic attenuation relationship with homogeneous and heterogeneous prediction-error variance models | |
CN104021288B (zh) | 用于导管架平台频谱疲劳分析的基本波确定方法 | |
CN107957982A (zh) | 震后次生地质灾害易发性快速评估方法和系统 | |
CN117235510A (zh) | 节理粗糙度预测方法与节理粗糙度预测模型训练方法 | |
CN115630577A (zh) | 一种降尺度降雨量预测方法及系统 | |
NAZRI et al. | Fragility curves | |
CN112214825B (zh) | 一种复杂结构抗震易损性输入地震波的选择方法 | |
Ryu et al. | A probabilistic method for the magnitude estimation of a historical damaging earthquake using structural fragility functions | |
Carrao et al. | Sample design and analysis for thematic map accuracy assessment: an approach based on domain estimation for the validation of land cover products | |
Minas et al. | Spectral shape proxies and simplified fragility analysis of mid-rise reinforced concrete buildings | |
Cruz | Evaluation of damping ratios inferred from the seismic response of buildings | |
Holmes | Seismic evaluation of existing buildings: state of the practice | |
CN105528657A (zh) | 基于北斗和向量机的建筑物震害预测方法 | |
CN109782349B (zh) | 一种用于结构抗震时程分析的选波方法及系统 |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230605 Address after: 053000 No. 989, Hengzao Road, Binhu New District, Hengshui City, Hebei Province Patentee after: HENGSHUI RUBBER GENERAL PLANT Co.,Ltd. Address before: 300000 5340 Xiping Road, Shuang Kou Town, Beichen District, Tianjin. Patentee before: Hebei University of Technology |