时程分析地震波选取问题

地震响应的反应谱法与时程分析比较 (1)

发电厂房墙体地震响应的反应谱法与时程分析比较 1问题描述 发电厂房墙体的基本模型如图1所示： 图1 发电厂墙体几何模型 基本要求：依据class 9_10.pdf的最后一页的作业建立ansys模型,考虑两个水平向地震波的共同作用(地震载荷按RG1.60标准谱缩放,谱值如下)，主要计算底部跨中单宽上的剪力与弯矩最大值,及顶部水平位移。要求详细的ansys反应谱法命令流与手算验证过程。以时程法结果进行比较。分析不同阻尼值(0.02,0.05,0.10)的影响。 RG1.60标准谱 (1g=9.81m/s2) (设计地震动值为0.1g) 频率谱值(g) 33 0.1 9 0.261 2.5 0.313 0.25 0.047 与RG1.60标准谱对应的两条人工波见文件rg160x.txt与rg160y.txt 2数值分析框图思路与理论简介 2.1理论简介 该问题主要牵涉到结构动力分析当中的时程分析和谱分析。时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。谱分析是模态分析的扩展，是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。 2.2 分析框架： 时程分析：在X和Z两个水平方向地震波作用下，提取底部跨中单宽上的剪力、弯矩值和顶部水平位移，并求出最大响应。 谱分析：先做模态分析，再求谱解，由于X和Z两个方向的单点谱激励，因此需进行两次谱分析，分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。 3有限元模型与荷载说明 3.1 有限元模型 考虑结构的几何特性建立有限元模型，首先建立平面几何模型，并将模型进行合理的切割，采用plane42单元，使用映射划分网格的方法生产平面单元（XOY平面）。然后，采用solid45

在 ansys 中如何 施加 地震波

三向输入简化后的单向输入 首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据： -0.227109E-02 -0.209046E+00 0.467072E+01 -0.413893E-02 -0.168195E+00 0.261523E+01 -0.574753E-02 -0.157890E+00 0.809014E-01 -0.731227E-02 -0.152996E+00 0.119975E+01 -0.876865E-02 -0.138102E+00 0.130902E+01 -0.101067E-01 -0.131582E+00 0.143611E+00 ....................... 然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据： 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01 ....................... 编写如下的命令流文件，并命名为acce.inp *dim,ACCEXYZ,TABLE,2000,3 !01行 *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行 *vread,ACCEXYZ(1,0),time,txt !04行 (e16.6) !05行 ACCEXYZ(0,1)=1 !06行 ACCEXYZ(0,2)=2 !07行，同上 ACCEXYZ(0,3)=3 !08行，同上 finish /SOLU ANTYPE,trans btime=0.01 !定义计算起始时间 etime=15.00 !定义计算结束时间 dtime=0.01 !定义计算时间步长 *DO,itime,btime,etime,dtime time,itime AUTOTS,0 NSUBST,1, , ,1 KBC,1 acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度 SOLVE

时程分析中地震波输入位置的讨论

时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[1]。目前有限元软件可以实现结构的时程分析，但是在不同的软件中，其实现方式不同，主要区别在地震波的输入位置不同。本文通过有限元软件ABAQUS采用不同的地震波输入位置对同一结构进行时程分析分析，对比结构相同位置的时程位移曲线，结果表明结构在采用不同地震波输入位置的时程分析中，结构的地震响应基本一致。 关键词：时程分析、有限元软件、钢筋混凝土剪力墙 Abstract: The time history analysis method to analyze the available structure through direct power to the relationship between the corresponding changes over time, truly reflect the structure of earthquake corresponding to the whole process of change over time, is an important method of seismic analysis [1]. Finite element software can be time-history analysis of the structure, but in different software in different ways, the main difference between the different positions in the seismic wave input. In this paper the finite element software ABAQUS using different seismic wave input location on the same structure, process analysis analysis, contrast structure the same location of when the process displacement curve, the results show that the structure using different seismic waves enter the position time history analysis, the seismic response basically the same. Keywords: time history analysis, finite element software, reinforced concrete shear walls 一、引言 在时程分析等动力学问题中，地震力以加速度形式从基础固定处输入。由于结构的刚度不是无限大，在结构上的加速度反应与基础输入的加速度并不相同。在很多时候，结构的加速度比基础输入的加速度更大，即对输入的加速度有一个动力放大效应。在单自由度弹性体系中，体系最大绝对加速度与地面运动最大加速度的比值，即称为动力系数[2] （1） 动力系数与结构的动力学特性和输入的地震波的频率特性有关。它与地震系数k的乘积即为单自由度体系的地震影响系数。 因此，从原理上讲，时程分析是将地震波的加速度时程曲线作用到结构的基础约束处，得到上部结构的各种地震反应。但是在不同的软件中，其实现方

地震反应谱分析实例

结构地震反应谱分析实例 在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑X与Y向地震效应作用。已知地震影响系数a与周期T的关系： a(T)= 0.4853*(0.4444+2.2222*T) 0

!进行模态求解 ANTYPE,MODAL MODOPT,LANB,30 SOLVE FINISH !进行谱分析 /SOLU ANTYPE,SPECTR SPOPT,SPRS,30,YES SVTYP,2 !加速度反应谱 SED,1,1 !X与Y向 FREQ,0.2500,0.2632,0.2778,0.2941,0.3125,0.3333,0.3571,0.3846,0.4167 FREQ,0.4545,0.5000,0.5556,0.6250,0.7143,0.8333,1.1111,2.0000,10.0000 FREQ,25.0000,1000.0000 SV,0.05,0.0797,0.0861,0.0934,0.1018,0.1114,0.1228,0.1362,0.1522,0.1716 SV,0.05,0.1955,0.2255,0.2642,0.3152,0.3851,0.4853,0.4853,0.4853,0.4853 SV,0.05,0.2588,0.2167 SOLVE FINISH !进行模态求解(模态扩展) /SOLU ANTYPE,MODAL EXPASS,ON MXPAND,30,,,YES,0.005 SOLVE FINISH !进行谱分析(合并模态) /SOLU ANTYPE,SPECTR SRSS,0.15,disp SOLVE FINISH /POST1 SET,LIST !结果1 /INP,,mcom

时程分析中地震波选取浅析

时程分析中地震波选取浅析 通过介绍时程分析法中输入地震波的选择原则、地震动幅值和频率特性等一系列问题，使初学者对输入地震波的选择有初步认识和了解，为以后更深层次的研究打下基础。 标签：时程分析法；地震波选择 1、引言 随着社会、经济和科技的不断发展以及人口数量的迅速膨胀，高层、超高层以及复杂形状的建筑的数量定会快速增长。抗震设计规范规定，对于此类重要、复杂并超过规定高度的建筑，其抗震设计中的地震作用计算都要通过时程分析法进行补充验证。而在时程分析法的计算过程中最重要，最影响地震作用计算结果的莫过于地震波的选取。所以，本文将从地震波选取原则、地震动幅值、频谱特性、持续时间、地震波数量、地震波转动分量等多个方面对地震波的选取进行浅析。 2、地震波的选取原则 时程分析中的地震波如何选取的问题，一直是时程分析法中的一个难点。在选择地震波输入时，要满足两点要求： 1）首先要使选择输入的地震波的某些参数和建筑物所在地的条件相一致。参数主要包括：场地的土壤类别、地震烈度、地震强度参数、卓越周期和反应谱等。 2）其次还要满足地震活动三要素的要求。即频谱特性、地震加速度时程曲线持续时间和幅值，选取的地震波中的这三者，要满足相关规定。相关规定要求：选用数字化的地震波应按照建筑场地类别和设计地震分组进行选取，选用不少于两组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱分析法所采用的地震影响系数曲线在统计意义上相符。在统计意义上相符是指：其平均地震影响曲线与振型分解反应谱法所用到的地震影响系数曲线相比，在各个周期点上相差不大于20%。弹性时程分析时，每条时程曲线计算所得的结构底部剪力不应小于阵型分解反应谱法计算结果的65%。多条时程曲线计算结果的结构底部剪力平均值不应小于振型分解反应谱计算结果的80%[1]。 3、地震动幅值 地震动幅值有两种意义，即可以指地震加速度、位移和速度中的任何一种的最大值，又可以指在某种意义下的等代值。在一定程度上，地震波的峰值能够反应并代表地震波的强度，所以，建筑物所在地的设防烈度所要求的多遇地震或罕

ANSYS地震反应谱SRSS分析共24页

ANSYS地震反应谱SRSS分析 我在ANSYS中作地震分解反应谱分析，一次X方向，一次Y 方向，他们要求是独立互不干扰的，可是采用直进行一次模态分析的话，他生成的*.mcom文件好像是包含了前面的计算 结果，命令流如下： !进入PREP7并建模 /PREP7 B=15 !基本尺寸 A1=1000 !第一个面积 A2=1000 !第二个面积 A3=1000 !第三个面积 ET,1,beam4 !二维杆单元 R,1,0.25,0.0052,0.0052,0.5,0.5 !以参数形式的实参 MP,EX,1,2.0E11 !杨氏模量 mp,PRXY,1,,0.3 mp,dens,1,7.8e3 N,1,-B,0,0 !定义结点 N,2,0,0,0 N,3,-B,0,b

N,4,0,0,b N,5,-B,0,2*b N,6,0,0,2*b N,7,-B,0,3*b N,8,0,0,3*b E,1,3 !定义单元 E,2,4 E,3,5 E,4,6 E,3,4 E,5,6 e,5,7 e,6,8 e,7,8 D,1,ALL,0,,2 FINISH ! !进入求解器，定义载荷和求解 /SOLU D,1,ALL,0,,2 !结点UX=UY=0

sfbeam,1,1,PRES,100000, sfbeam,3,1,PRES,100000, sfbeam,7,1,PRES,100000, SOLVE FINISH allsel NMODE=10 /SOL !* ANTYPE,2 !* MSAVE,0 !* MODOPT,LANB,NMODE EQSLV,SPAR MXPAND,NMODE , , ,1 LUMPM,0 PSTRES,0 !* MODOPT,LANB,NMODE ,0,0, ,OFF

反应谱与时程理论对比

反应谱是在给定的地震加速度作用期间内，单质点体系的最大位移反应、速度反应和加速度反应随质点自振周期变化的曲线。用作计算在地震作用下结构的内力和变形。更直观的定义为：一组具有相同阻尼、不同自振周期的单质点体系，在某一地震动时程作用下的最大反应，为该地震动的反应谱。 反应谱理论考虑了结构动力特性与地震动特性之间的动力关系，通过反应谱来计算由结构动力特性(自振周期、振型和阻尼)所产生的共振效应，但其计算公式仍保留了早期静 力理论的形式。地震时结构所受的最大水平基底剪力，即总水平地震作用为: FEK= αG 其中α为地震影响系数，即单质点弹性体系在地震时最大反应加速度。另一方面地震影响系数也可视为作用在质点上的地震作用与结构重力荷载代表值之比。 目前，反应谱分析法比较成熟，一些主要国家的抗震规范均将它作为基本设计方法。不过，它主要适合用于规则结构。对于不规则结构以及高层建筑，各国规范多要求采用时程分析法进行补充计算。 地震作用反应谱分析本质上是一种拟动力分析，它首先使用动力法计算质点地震响应，并使用统计的方法形成反应谱曲线，然后使用静力法进行结构分析。但它并不是结构真实的动力响应分析，只是对于结构动力响应最大值进行估算的近似方法，在线弹性范围内，反应谱分析法被认为是高效而且合理的方法。反应谱分为加速度反应谱、速度反应谱和位移反应谱。基于不同周期结构相应峰值的大小，我们可以绘制结构速度及加速度的反应谱曲线。一般情况下，随着周期的延长，位移反应谱为上升曲线，速度反应谱为平直曲线，加速度反应谱为下降曲线，目前结构设计主要依据加速度反应谱。 加速度反应谱在短周期部分为快速上升曲线，并且在结构周期与场地特征周期接近时出现峰值，后面更大范围为逐渐下降阶段。峰值出现的时间与对应的结构周期和场地特征周期有关。一般来说结构自振周期的延长，地震作用将减小。当结构自振周期接近场地特征周期时，地震作用最大。 反应谱分析方法需要先求解一个方向地震作用响应，再基于三个正交方向的分量考虑结构总响应，即基于振型组合求解一个方向的地震响应，再基于方向组合求解结构总响应。 振型组合方法有SRSS法，CQC法。 1.SRSS法 SRSS法是平方和平方根法，这种方法假定所有最大模态值在统计上都是相互独立的，通过求各参与阵型的平方和平方根来进行组合。该法不考虑各振型间的藕联作用，实际上结构模态都是相互关联的，不可避免的存在藕联效应，对那些相邻周期几乎相等的结构，或者不规则结构不适用此法。《抗规》GB50011-2010规定的SRSS法为如下所示：

地震波使用说明

地震波使用说明 此目录下提供了四类场地土的地震波时程曲线和上海人工波。 按照场地土类型（1，2，3或4），选择时程曲线。在定义时程工况时，对于多遇或罕遇地震，按比例调整时程曲线的最大值。中国抗震规范规定，作为抗震计算中底部剪力法和振型分解反应谱法的补充方法，对于特别不规则，特别重要的和较高的结构应采用时程分析法进行多遇地震下的补充计算。 可取多条时程曲线的计算结果的平均值与振型分解反应谱法计算结果的较大值。 采用时程分析法时，应咱建筑场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。 其加速度时程最大值可按规范中对于多遇和罕遇地震在不同烈度下的值。 弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65％，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80％ 。 可使用弹塑性时程分析法计算罕遇地震下结构的变形。 时程分析是一个承受随时间变化的指定荷载结构的逐步动态反应分析，可以是线性或非线性的。 此章对时程分析进行一般的描述，特别是线性时程分析。 定义时程函数 用户可使用“从文件中添加函数”，导入已定义的文本文件，即实测的时程曲线；也可使用程序内置的时程函数。

时程函数定义对话框 时程函数定义对话框中的条目解释如下： ?函数名 通过在编辑框中直接键入以指定或修改时程函数的名称。 ?函数文件 1.在函数文件域点击浏览按钮以调出一个对话框，在此可找出包含时程函数的 文本文件名。注意文件名显示在文件名框中 2.在 "要跳过的标题行" 编辑框中输入一个希望ETABS在文本文件中跳过的 行数。 3.在 "每行要跳过的前缀字符" 编辑框中输入一个希望ETABS在文本文件中 每行要跳过的字符数。 4.在 "每行的点数" 编辑框中输入一个数告诉ETABS文本文件每行的绘图点 数。

ANSYS地震分析实例

ANSYS地震分析实例 土木工程中除了常见的静力分析以外，动力分析，特别是结构在地震荷载作用下的受力分析，也是土木工程中经常碰到的题目。结构的地震分析根据现行抗震规范要求，一般分为以下两类：基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程分析。 本算例将以一个4质点的弹簧－质点体系来说明如何使用有限元软件进行地震分析。更复杂结构的分析其基本过程也与之类似。 关键知识点： (a) 模态分析 (b) 谱分析 (c) 地震反应谱输进 (d) 地震时程输进 (e) 时程动力分析 (1) 在ANSYS窗口顶部静态菜单，进进Parameters菜单，选择Scalar Parameters选项，在输进窗口中填进DAMPRATIO=0.02，即所有振型的阻尼比为2% (2) ANSYS主菜单Preprocessor->Element type->Add/Edit/Delete，添加Beam 188单元 (3) 在Element Types窗口中，选择Beam 188单元，选择Options，进进Beam 188的选项窗口，将第7个和第8个选项，Stress/Strain (Sect Points) K7, Stress/Strain (Sect Nods) K8，从None 改为Max and Min Only。即要求Beam 188单元输出积分点和节点上的最大、最小应力和应变 (4) 在Element Types 窗口中，继续添加Mass 21集中质量单元 (5) 下面输进材料参数，进进ANSYS主菜单Preprocessor->Material Props-> Material Models菜单，在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic 属性，输进材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。 (6) 继续给Material Model Number 1添加Density属性，输进密度为7800。 (7) 继续给Material Model Number 1添加Damping属性，采用参数化建模，输进阻尼类型为Constant，数值为DAMPRATIO

汶川地震波加速度时程

FILTER POINTS: NPTS= 25000, DT= .0200 SEC 单位cm/s2 -85.18660303 -26.81310797 29.69681739 17.66375925 10.74725598 5.89307633 10.53329077 8.38328838 -6.60304626 -7.29547636 11.29502428 13.67286273 -5.92602216 -23.83117937 -24.01378573 -10.41018086 -1.73061627 -1.68681541 -0.25120385 5.69817713 22.27844984 29.26687041 14.37719476 -0.78176094 14.85243403 43.59400018 29.38438838 3.43393402 5.49481755 6.23200112 -21.87622974 -52.48408787 -47.05050653 -6.35843868 26.04376039 16.94997018 8.10590549 14.56990612 8.61432565 -0.67444454 1.37721874 -3.94592499

14.4976265 15.78556387 8.14483834 -6.35078639 -12.80601727 -7.48203698 -38.05952669 -47.252601 -22.58282251 4.85801197 24.94937343 15.78910797 -7.89294083 -27.34878362 -7.67345164 21.73423697 33.04934899 5.07072495 5.38605227 29.29248964 32.47681472 63.38543891 78.60146871 41.15192422 16.98456088 24.67915562 29.25445807 -12.85852194 -52.13522708 -27.93999245 9.73808504 20.14985389 -1.49063292 -17.18407036 -9.30291165 -0.51966159 7.96338525 -7.54802444 -11.34591285 19.42029186 0.6700927 -5.71160065 -26.02061476

abaqus如何施加地震波

施加地震波： 1 *amplitude,name=amp,input=seismicdata.dat 输入地震波 2 *boundary,type=acceleration,amplitude=amp施加荷载 方法：module选load，在tools-----amplitude-----creat默认的continue在Edit A mplitude里面输入时间和加速度，点OK。点creat boundary condition，涌现对 话框creat boundary condition，选择acceleration/angular acceleration,continu e---选择要施加的边界---done----涌现对话框edit bondary condition对话框，在 amplitude里选择你所定义的时间和加速度。点ok就完工了。 在网上查了些方法： module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude 里面输入时间和加速度，点OK。点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界---done----出现对话框edit bondary condition对话框，在amplitude里选择你所定义的时间和加速度。点ok就完工了。 这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。 首先在CAE里建好模型，定义两个分析步。 第一个分析步是加自重，采用线性加载的方式。 (a) 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous（瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。 (b)新建加载方式：创建一个新的Amplitude，Type=smooth tpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smooth tpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。 第二个分析步就是加地震波。 输入地震波有两种方法： 1、在如下位置加入下面加黑的字体部分。格式如下：时间,地震波,时间,地震波,时间, 地震波,时间, 地震波…………每行8个数据（我下到的地震波文件是不带时间的，自己用C++处理了一下）。%%%%%%%%%%%%%%%%%%%%%% *End Assembly *Amplitude, name=Amp-1 0.005, -7.5e-08, 0.01, -3.55e-07, 0.015, -7.03e-07, 0.02, -4.53e-07 0.025, 1.82e-06, 0.03, 7.01e-06, 0.035, 1.5e-05, 0.04, 2.49e-05 0.045, 3.54e-05, 0.05, 4.5e-05, 0.055, 5.2e-05, 0.06, 5.5e-05 ………………

ANSYS地震波的输入

对于地震波的输入，可以把荷载记录作成文件，利用apdl的读取功能读入数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。 Fini /config,nres,1000 *dim,aceX,TABLE,3000,1 *dim,aceY,TABLE,3000,1 *dim,aceZ,TABLE,3000,1 *creat,ff *vread,aceX(1,1),acex,txt,,1 (e16.6) *vread,aceX(1,0),acexTT,txt,,,1 (e16.6) ACEX(0,1)=1 *end /input,ff *creat,ff *vread,aceY(1,1),txt,,1 (e16.6) *vread,aceY(1,1),ACETT,,,1 (e17.6) ACEY(0,1)=1 *end /input,ff *creat,ff *vread,aceZ(1,1),txt,,1 (e16.6) *vread,aceZ(1,0),ACETT,,,1 (e17.6) ACEZ(0,1)=1 *end /input,ff !地震波时程记录分成了3个文件，每个文件是一列。分别记录x，y，z方向的加速度。Accett是时间记录。 这样就可以把加速度记录读取到ansys数据库中作为数组。 也可以把加速度记录作成一个文件，这样程序就简单多了。 下面是计算部分语句： /SOLU ANTYPE,trans !求解其自己选了 TM_START=0.01 TM_END=15.00 TM_INCR=0.01 *DO,TM,TM_START,TM_END,TM_INCR TIME,tm

结构抗震设计时程分地震波的选择

（1）设计用地震记录的选择和调整 用规范的确定性方法和地震危险性分析方法所确定的设计地震动参数，是选择天然地震加速度记录的依据。 (一)实际地震记录的选择方法 选择地震记录应考虑地震动三要素，即强度（峰值）、频谱和持续时间。对某一建筑的抗震设计，最好是选用该建筑所在场地曾经记录 到的地震加速度时间过程。但是，这种机会极少。为此，人们只能从现有的国内外常用的地震记录中去选择，尽可能挑选那些在震级、震中距和场地条件等方面都比较接近设计地震动参数的记录。他的文章给出了相应的地震数据的记录目录。 （二）实际地震记录的调整 1.强度调整。将地震记录的加速度值按适当的比例放大或缩小，使其峰值加速度等于事先所确定的设计地震加速度峰值。即令 其中a(为记录的加速度值为调整后的加速度值;A众为设计地震加速度峰值;。为记录的加速度峰值。这种调整只是针对原记录的强度进行的，基本上保留了实际地震记录的特征。也就是所说的（强度修正。将地震波的加速度峰值及所有的离散点都按比例放大或缩小以满足场地的烈度要求）

2.频率调整考虑到场地条件对地震地面运动的影响，原则上所选择的实际地震记录的富氏谱或功率谱的卓越周期乃至形状，应尽量与场地土相应的谱的特性一致。如果不一致，可以调整实际地震记录的时间步长，即将记录的时间轴“拉长”或“缩短”，以改变其卓越周期而加速度值不变也可以用数字滤波的方法滤去某些频率成分，改变谱的形状。另外，为了在计算中得到结构的最大反应，也可以根据建筑结构基本自振周期，调整实际地震记录的卓越周期，使二者接近。这种调整的结果，改变了实际地震记录的频率结构，从物理意义上分析是不合理的。 另外，在测定场地土和建筑结构的卓越周期时，运用不同的测试仪器和测试技术，往往得到不同的结果。即使是对同一个测试结果，在频谱上确定卓越周期时，不同的分析方法也会导致不同的结果。有的选取谱的第一个峰值所对应的周期作为卓越周期，有的选最大峰值时的，也有的取某一段周期等，很不一致。对如何确定地震加速度记录的卓越周期，也是各行其是，有的利用加速度反应谱，有的用伪速度谱，有的用富氏谱，结果当然是不一样的。上述各种作法在工程中引起了一些混乱。 王亚勇认为，用脉动测试方法测定场地土和结构的卓越周期及自振周期时，应采用速度摆型或加速度摆型的地震仪测定地运动和结构振动，然后计算其富氏谱或功率谱，以谱的最大峰值所对应的周期作为卓越周期和自振周期比较合适。反应而相应地根据记录的位移谱或速度谱。 这也就是所谓的滤波修正。可按要求设计滤波器，对地震波进行时域或频域的滤波修正。这样修正的地震资料不仅卓越周期满足要求，功率谱的形状和面积也可控制。卓越周期修正。将地震波的离散步长按人为比例改变，

三 设计地震动反应谱确定的规范方法

三设计地震动反应谱确定的规范方法 设计地震动是通过对地震环境和场地环境的分析判断和分类方法确定。工程勘察单位至少提供： 设计基本地震加速度和设计特征周期 场地环境：覆盖层厚度、剪切波速、土层钻孔资料 1.设计基本地震加速度和设计特征周期 根据场地在中国地震动参数区划图上的位置判断确定。

土层剪切波速的测量应符合下列要求: 1 在场地初步勘察阶段对大面积的同一地质单元测量土层剪切波速的钻孔数量不宜少于3。 2 在场地详细勘察阶段对单幢建筑测量土层剪切波速的钻孔数量不宜少于2 个数据变化较大时可适量增加对小区中处于同一地质单元的密集高层建筑群测量土层剪切波速的钻孔数量可适量减少但每幢高层建筑下不得少于一个。 3 对丁类建筑及层数不超过10 层且高度不超过30m 的丙类建筑当无实测剪切波速时可根据岩土名称和性状按表 4.1.3 划分土的类型再利用当地经验在下表的剪切波速范围内估计各土层的剪切波速.

建筑场地覆盖层厚度的确定应符合下列要求: 1 一般情况下应按地面至剪切波速大于500m/s 的土层顶面的距离确定（且其下卧层沿途的剪切波速均不小于500m/s）。 2 当地面5m 以下存在剪切波速大于（其上部各土层）相邻上层土剪切波速2.5 倍的土层且其下卧岩土的剪切波速均不小于400m/s 时可按地面至该土层顶面的距离确定 3 剪切波速大于500m/s 的孤石、透镜体应视同周围土层 4.土层中的火山岩硬夹层应视为刚体其厚度应从覆盖土层中扣除

例题：某类建筑场地位于7度烈度区，设计地震分组为第一组，设计基本地震加速度为0.1g，建筑结构自振周期T=1.4s，阻尼比为0.08，该场地在建筑多遇地震条件下地震影响系数a为多少。 同一个场地上甲乙两座建筑物的结构自震周期分别为T甲=0.25sT乙=0.60s，一建筑场地类别为Ⅱ类，设计地震分组为第一组，若两座建筑的阻尼比都取0.05，问在抗震验算时甲、乙两座建筑的地震影响系数之比最接近下列那个选项。 A 1.6 B 1.2 C 0.6 D 条件不足无法计算 例题：吉林省松原市某民用建筑场地地质资料如下： （1）0-5m粉土，=150 =180m/s (2) 5-12m中砂土=200 =240m/s （3）12-24m粗砂土=230 =310m/s (4) 24-45m硬塑粘土=260 =300m/s （5）45-60m泥岩=500 =520m/s 建筑物采用浅基础，埋深2m，地下水位2.0m，阻尼比为0.05，自震周期为1.8s该建筑进行抗震设计时 （1）进行第一阶段设计时，地震影响系数应取多少 （2）进行第二阶段设计时，地震影响系数应取多少 例题：吉林省松原市某民用建筑场地地质资料如下： （1）0-5m粉土，=150 =180m/s (2) 5-12m中砂土=200 =240m/s

地震反应谱的绘制

地震时程曲线与反应谱的绘制 ①地震反应谱的意义 地震反应谱表示的是在一定的地震动下结构的最大反应，是结构进行抗震分析与设计的重要工具。 由于同一结构在遭遇不同的地震作用时的反应并不相同，单独一个地震记录的反应谱不能用于结构设计。但是地震记录的反应谱又有一定的相似性，我们可以将具有普遍特性记录的反应谱进行平均和平滑处理，以用于抗震设计。现在，地震反应谱不但是工程抗震学中最重要的概念之一，还是整个地震工程学中最重要的概念之一。 ②地震反应谱的计算方法 反应谱的计算方法涉及到时域分析方法和频域分析方法。 时域分析方法中的Duhamel 积分，是现在公认精度最高的方法。 绝对加速度反应谱公式如下：（推导略） 但由于实际结构系统的阻尼比ξ通常都小于0.1，所以有阻尼系统和无阻尼系统的自振 周期ω近似相等即由ωζω21-=d （精确度≥99.5%）简化成ωω=d ，实际计算中通常按无阻尼系统的自振周期确定。 从而上式可以简化为 ()()()max 00max sin )(?-==--t t a d t e x t a S ττωτωτζω ③用matlab 画地震时程曲线与绝对加速度反应谱： 所需准备软件： excel ，notepad2，matlab 以NINGHE 地震波为例 Code ： %NINGHE 地震波时程曲线 % 加载前用excel 和notepad 对数据进行规整

load NINGHE.txt; % 数据放在安装文件的work目录下 NUMERIC=transpose(NINGHE); % matlab read the data by column, ni=reshape(NUMERIC,numel(NUMERIC),1);% make the date one column t_ni=0:0.002:(length(ni)-1)*0.002; % determine the time plot(t_ni,ni); ylabel('Acceleration'); xlabel('time'); title('NINGHE') %NINGHE绝对加速度反应谱 load NINGHE.txt; NUMERIC=transpose(NINGHE); ni=reshape(NUMERIC,numel(NUMERIC),1);%make the date one column d=0;%d is damping ratio for k=1:600; t(k)=0.01*k;%规范的加速度反应谱只关心前6秒的值 w=6.283185/t(k); t_ni=0:0.02:(length(ni)-1)*0.02; Hw=exp(-1*d*w*t_ni).*sin(w*t_ni); y1=conv(ni,Hw).*(0.02*w);y1=max(abs(y1));%卷积积分 c(k)=y1*10; end;plot(t,c,'black')

时程分析时地震波的选取及地震波的反应谱化

时程分析时地震波的选取及地震波的反应谱化 摘要：目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应 谱法。时程分析法作为补充计算方法，在不规则、重要或较高建筑中采用。进行 时程分析时，首先面临正确选择输入的地震加速度时程曲线的问题。时程曲线的 选择是否满足规范的要求，则需要首先将时程曲线进行单自由度反应计算，得到 其反应谱曲线，并按规范要求和规范反应谱进行对比和取舍。本文通过介绍常用 的数值计算方法及计算步骤，实现将地震加速度时程曲线计算转化成反应谱曲线，从而为特定工程在时程分析时地震波的选取提供帮助。 关键词：时程分析，地震波，反应谱，动力计算 1 地震反应分析方法的发展过程 结构的地震反应取决于地震动和结构特性。因此，地震反应分析的水平也是随着人们对 这两个方面认识的深入而提高的。结构地震反应分析的发展可以分为静力法、反应谱法、动 力分析法这三个阶段。在动力分析法阶段中又可分为弹性和非弹性（或非线性）两个阶段。[1] 目前，在我国和其他许多国家的抗震设计规范中，广泛采用反应谱法确定地震作用，其 中以加速度反应谱应用得最多。反应谱是指：单自由度弹性体系在给定的地震作用下，某个 最大反应量（如加速度、速度、位移等）与体系自振周期的关系曲线。反应谱理论是指：结 构物可以简化为多自由度体系，多自由度体系的地震反应可以按振型分解为多个单自由度体 系反应的组合，每个单自由度体系的最大反应可以从反应谱求得。其优点是物理概念清晰， 计算方法较为简单，参数易于确定。 反应谱理论包括如下三个基本假定：1、结构物的地震反应是弹性的，可以采用叠加原理 来进行振型组合；2、现有反应谱假定结构的所有支座处地震动完全相同；3、结构物最不利 的地震反应为其最大地震反应，而与其他动力反应参数，如最大值附近的次数、概率、持时 等无关。[1] 时程分析法是对结构物的运动微分方程直接进行逐步积分求解的一种动力分析方法。由 于此法是对运动方程直接求解，又称直接动力分析法。可直接计算地震期间结构的位移、速 度和加速度时程反应，从而描述结构在强地震作用下弹性和非弹性阶段的内力变化，以及结 构构件逐步开裂、屈服、破坏甚至倒塌全过程。 根据我国《建筑抗震设计规范》(GB5011-2010)(以下简称《抗规》)第5.1.2-3条要求，特 别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑，应采用时程分析法进行多 遇地震下的补充计算。此外《高层建筑混凝土结构技术规程》(JGJ3-2010) (以下简称《高规》)第4.3.4条也有相关要求。 2 时程分析时地震波的选取要求 在进行时程分析时，首先面临地震波选取的问题。所选的地震波需要符合场地条件、设 防类别、震中距远近等因素。《抗规》对于地震波的选取主要有以下几点要求： 1、当取三组加速度时程曲线输入时，计算结果宜取时程法的包络值和振型分解反应谱法 的较大值；当取七组及七组以上的时程曲线时，计算结果可取时程法的平均值和振型分解反 应谱法的较大值（其中实际强震记录的数量不应少于总数的2/3）。 2、弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计 算结果的65%，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计 算结果的80%。 3、多组时程曲线的平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数 曲线在统计意义上相符。根据规范条文说明，所谓“统计意义上相符”指的是，多组时程波的 平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在对应于结构主 要振型的周期点上相差不大于20%。但计算结果也不能太大，每条地震波输入计算不大于135%，平均不大于120%。 4、时程曲线要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间均要符合规

ABAQUS地震反应谱分析

ABAQUS反应谱法计算地震反应的简单实例 Fan.hj 2010年4月4日 清明小长假，琢磨了下ABAQUS如何进行地震反应谱计算。现通过一小算例说明。 问题描述： （本例的问题引用《有限元法及其应用》一书中陆新征博士ANSYS算例的问题） 悬臂柱高12m，工字型截面（图1），密度7800kg/m3，EX=2.1e11Pa，泊松比0.3，所有振型的阻尼比为2%，在3m高处有一集中质量160kg，在6m、9m、12m处分别有120kg的集中质量。反应谱按7度多遇地震，取地震影响系数为0.08，第一组，III类场地，卓越周期Tg=0.45s。 图1 计算对象 几点说明： ●本例建模过程使用CAE； ●添加反应谱必须在inp中加关键词实现，CAE不支持反应谱； ●*Spectrum不可以在keyword editor中添加，keyword editor不支持此关键词读入； ●ABAQUS的反应谱法计算过程以及后处理要比ANSYS方便的多。 操作过程为： （1）打开ABAQUS/CAE，点击create model database。 （2）进入Part模块，点击create part，命名为column，3D、deformation、wire。OK （3）Create lines：connected，分别输入0，0；0，3；0，6；0，9；0，12。OK。退出sketch。（4）进入property模块，create material，name：steel，general-->>density，mass density：7800，mechanical-->>elasticity-->>elastic，young‘s modulus：2.1e11，poisson’s ratio： 0.3.OK