地震波应用

建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定，正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。

频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同，简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。

加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。

持续时间的概念不是指地震波数据中总的时间长度。持时Td的定义可分为两大类，一类是以地震动幅值的绝对值来定义的绝对持时，即指地震地面加速度值大于某值的时间总和，即绝对值｜a(t)｜＞k*g的时间总和，k常取为0.05；另一类为以相对值定义的相对持时，即最先与最后一个k*amax之间的时段长度，k一般取0.3～0.5。不论实际的强震记录还是人工模拟波形，一般持续时间取结构基本周期的5～10倍。

说明:

有效峰值加速度 EPA＝Sa/2.5 （1）

有效峰值速度 EPV＝Sv/2.5 （2）

特征周期Tg = 2π*EPV/EPA （3）

1978年美国ATC－3规范中将阻尼比为5％的加速度反应谱取周期为0.1-0.5秒之间的值平均为Sa，将阻尼比为5%的速度反应谱取周期为0.5-2秒之间的值平均为Sv(或取1s附近的平均速度反应谱)，上面公式中常数2.5为0.05组尼比加速度反应谱的平均放大系数。

—1 —

上述方法使用的是将频段固定的方法来求EPA和EPV，1990年的《中国地震烈度区划图》采用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是：在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱，找出加速度反应谱平台段的起始周期T0和结束周期T1，然后在拟速度反应谱上选定平台段，其起始周期为T1(即加速度反应谱平台段的结束周期T1)，结束周期为T2，将加速度反应谱在T0至T1之间的谱值求平均得Sa，拟速度反应谱在T1至T2之间的谱值求平均得Sv，加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5，按公式(1)、（2）、（3）求得EPA、EPV、Tg。

在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功能(工具>地震波数据生成器，生成后保存为SGS文件)，用户可利用保存的SGS文件(文本格式文件)根据上面所述方法计算Sv、Sa、Tg。通过Tg值可判断该地震波是否适合当地场地和地震设计分组，然后将抗震规范中表5.1.2-2中的EPA值与Sa相比求出调整系数，将其代入到地震波调整系数中。将地震波转换为绝对加速度反应谱和拟速度反应谱时注意周期范围要到6秒(建筑抗震规范规定)。

建筑抗震设计规范5.1.2条中规定，采用时程分析方法时，应按照场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。所谓“在统计意义上相符”指的是，其平均影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在各周期点上相差不大于20%。

在MIDAS程序中，可选取两组实际强震记录生成两个SGS文件(调整Sa后的)，然后将一组人工模拟的加速度时程曲线也保存为SGS文件，将三个SGS文件的数值取平均后与振型分解反应谱法所采用的地震影响系数曲线相比较看是否满足“在统计意义上相符”，由此也可判断选取的地震波是否合适。

—2 —

另外，弹性时程分析时，每条时程曲线计算所得到的结构底部剪力不应小于振型分解反应谱法计算结果的65%，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80%。

建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2中的值即为加速度有效峰值。

表5.1.2-2 时程分析所用地震加速度时程曲线的最大值 (cm/s^2) 地震影响 6度7度 8度 9度

多遇地震 18 35(55) 70(110) 140

罕遇地震- 220(310) 400(510) 620

注: 括号内数值分别用于设计地震加速度为0.15g和0.30g的地区。

另外，表 5.1.2-2中的数值(加速度有效峰值)可使用水平地震影响系数最大值(表5.1.4-1)除2.25得到(5.1.4条文说明)。

表5.1.4-1 水平地震影响系数最大值

地震影响6度7度8度 9度

多遇地震 0.04 0.08(0.12) 0.16(0.24) 0.32

罕遇地震 - 0.50(0.72) 0.90(1.20) 1.40

注: 括号内数值分别用于设计地震加速度为0.15g和0.30g的地区。

举例说明: 7度多遇地震的水平地震影响系数最大值为0.08，重力加速度为9.801m/s^2，所以7度多遇地震区域的加速度有效峰值为0.08*9.801/2.25=0.348 m/s^2(近似为 35 cm/s^2)。

—3 —

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

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

地震波运动学理论

第二章地震波运动学理论 一、名词解释 1. 地震波运动学：研究在地震波传播过程中的地震波波前的空间位置与其传播时间的关系，即研究波的传播规律，以及这种时空关系与地下地质构造的关系。 2. 地震波动力学：研究地震波在传播过程中波形、振幅、频率、相位等特征的及其变化规律，以及这些变化规律与地下的地层结构，岩石性质及流体性质之间存在的联系。 3. 地震波：是一种在岩层中传播的，频率较低(与天然地震的频率相近)的波，弹性波在 岩层中传播的一种通俗说法。地震波由一个震源激发。 4. 地震子波：爆炸产生的是一个延续时间很短的尖脉冲，这一尖脉冲造成破坏圈、塑性带，最后使离震源较远的介质产生弹性形变，形成地震波，地震波向外传播一定距离后，波形逐渐稳定，成为一个具有2-3个相位（极值）、延续时间60-100毫秒的地震波，称为地震子波。地震子波看作组成一道地震记录的基本元素。 5.波前：振动刚开始与静止时的分界面，即刚要开始振动的那一时刻。 6.射线：是用来描述波的传播路线的一种表示。在一定条件下，认为波及其能量是沿着一条“路径”从波源传到所观测的一点P。这是一条假想的路径，也叫波线。射线总是与波阵面垂直，波动经过每一点都可以设想有这么一条波线。 7. 振动图和波剖面：某点振动随时间的变化的曲线称为振动曲线，也称振动图。地震勘探中，沿测线画出的波形曲线，也称波剖面。 8. 折射波：当入射波大于临界角时，出现滑行波和全反射。在分界面上的滑行波有另一种特性，即会影响第一界面，并激发新的波。在地震勘探中，由滑行波引起的波叫折射波，也叫做首波。入射波以临界角或大于临界角入射高速介质所产生的波 9.滑行波：由透射定律可知，如果V2>V1 ，即sinθ2 > sinθ1 ，θ2 > θ1。当θ1还没到90o时，θ2 到达90o，此时透射波在第二种介质中沿界面滑行，产生的波为滑行波。 10.同相轴和等相位面：同向轴是一组地震道上整齐排列的相位，表示一个新的地震波的到达，由地震记录上系统的相位或振幅变化表示。 11.地震视速度：当波的传播方向与观测方向不一致(夹角θ)时，观测到的速度并不是波前的真速度V，而是视速度Va。即波沿测线方向传播速度。 12 波阻抗：指的是介质(地层)的密度和波的速度的乘积(Zi=ρiVi，i为地层)，在声学中称为声阻抗，在地震学中称波阻抗。波的反射和透射与分界面两边介质的波阻抗有关。只有在Z1≠Z2的条件下，地震波才会发生反射，差别越大，反射也越强。 13.纵波：质点振动方向与波的传播方向一致，传播速度最快。又称压缩波、膨胀波、纵波或P-波。 14.横波：质点振动方向与波的传播方向垂直，速度比纵波慢，也称剪切波、旋转波、横波或S-波，速度小于纵波约0.7倍。横波分为SV和SH波两种形式。 15.体波：波在无穷大均匀介质(固体)中传播时有两种类型的波（纵波和横波），它们在介质的整个立体空间中传播，合称体波。 16共炮点反射道集：在同一炮点激发，不同接收点上接收的反射波记录，称为共炮点道集。在野外的数据采集原始记录中，常以这种记录形式。可分单边放炮和中间放炮。 17.面波：波在自由表面或岩体分界面上传播的一种类型的波。 18.纵测线和非纵测线：激发点与接收点在同一条直线上，这样的测线称为纵测线。用纵测线进行观测得到的时距曲线称为纵时距曲线。激发点不在测线上，用非纵测线进行观测得到的时距曲线称为非纵时距曲线。

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

时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[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的乘积即为单自由度体系的地震影响系数。 因此，从原理上讲，时程分析是将地震波的加速度时程曲线作用到结构的基础约束处，得到上部结构的各种地震反应。但是在不同的软件中，其实现方

论地震勘探中几种主要地震波

论地震勘探中的几种主要地震波 论文提要 地震勘探，就是通过人工方法激发地震波，研究地震波在地层中传播的情况，以查明地下地质构造，为寻找油气田或其它勘探目的服务的一种方法。也可以理解为就是利用地震子波从地下地层界面反射回地面时带回来的旅行时间和形状变化的信息，用以推断地下的底层构造和岩性。地震勘探在勘探已有的各种物探方法中，是最有效地方法。在地震勘探中用炸药激发时，一声炮响之后会产生各种各样的地震波。按波在传播过程中质点震动的方向来区分，可以纵波和横波；根据波动所能传播的空间范围而言，地震波又可以分为体波和面波；按照波在传播过程中的传播路径的特点，又可以把地震波分为直达波、反射波、透射波、折射波，等等。地震勘探在石油勘探中除了能产生来自地层界面有用的反射波外，还会产生各种各样的干扰波。因此，我们要更好的了解各种波的产生、特点、用途，等等。下面简单介绍几种地震勘探中产生的地震波。 正文 一、反射波 （一）反射波的形成 1、几何地震学的观点 当炸药在井中爆炸激发地震波时，在雷管引爆几百微妙之内爆炸便完成了，在接近爆炸点的压强是一个延续时间很短的尖脉冲，爆炸脉冲向外传播，压强逐渐减少，地层开始产生弹性形变，形成地震波。地震波继续传播，由于介质对高频的吸收，地震波信号减小。当波入射到两种介质的分界面时（当上层介质波阻抗与下层介质波阻抗不等时，弹性地震波才会发生反射；上层介质波阻抗与下层介质波阻抗差别越大，反射波越强——反射波条件），一部分波回到第一种介质中，这就是所谓的反射波。如图所示 2、物理地震学观点 地震波从震源出发以球面波的方式向下传播，到达反射界面S，S可以就看成有许多

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

Midas地震波的选取方法

地震波的选取方法 建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定，正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。 频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同，简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。 加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。 持续时间的概念不是指地震波数据中总的时间长度。持时T d的定义可分为两大类，一类是以地震动幅值的绝对值来定义的绝对持时，即指地震地面加速度值大于某值的时间总和，即绝对值｜a(t)｜＞k*g的时间总和，k常取为0.05；另一类为以相对值定义的相对持时，即最先与最后一个k*a max之间的时段长度，k一般取0.3～0.5。不论实际的强震记录还是人工模拟波形，一般持续时间取结构基本周期的5～10倍。 说明: 有效峰值加速度 EPA＝Sa/2.5 （1） 有效峰值速度 EPV＝Sv/2.5 （2） 特征周期 Tg = 2π*EPV/EPA（3） 1978年美国ATC－3规范中将阻尼比为5％的加速度反应谱取周期为0.1-0.5秒之间的值平均为Sa，将阻尼比为5%的速度反应谱取周期为0.5-2秒之间的值平均为Sv(或取1s附近的平均速度反应谱)，上面公式中常数2.5为0.05组尼比加速度反应谱的平均放大系数。 上述方法使用的是将频段固定的方法来求EPA和EPV，1990年的《中国地震烈度区划图》采用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是：在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱，找出加速度反应谱平台段的起始周期T0和结束周期T1，然后在拟速度反应谱上选定平台段，其起始周期为T1(即加速度反应谱平台段的结束周期T1)，结束周期为T2，将加速度反应谱在T0至T1之间的谱值求平均得Sa，拟速度反应谱在T1至T2之间的谱值求平均得Sv（注：生成谱的时候一定要用对数谱），加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5，按公式(1)、（2）、（3）求得EPA、EPV、Tg。 在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功能(工具>地震波数据生成器，生成后保存为SGS文件)，用户可利用保存的SGS文件(文本格式文件)根据上面所述方法计算Sv、Sa、Tg＝Sv/Sa。通过Tg值可判断该地震波是否适合当地场地和地震设计分组，然后将抗震规范中表5.1.2-2中的EPA值与Sa相比求出调整系数（即放大系数），将其代入到地震波调整系数中。将地震波转换为绝对加速度反应谱和拟速度反应谱时注意周期范围要到6秒(建筑抗震规范规定)。 建筑抗震设计规范5.1.2条中规定，采用时程分析方法时，应按照场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。所谓“在统计意义上相符”指的是，其平均影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在各周期点上相差不大于20%。 在MIDAS程序中，可选取两组实际强震记录生成两个SGS文件(调整Sa后的)，然后将一组人

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

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

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

地震波传播原理

菲涅尔体和透射波 摘要 在地震成像实验中，通常使用基于波动方程高频渐进解的几何射线理论，因此，通常假设地震波沿着空间中一条连接激发点和接受点的无限窄的线传播，称为射线。事实上，地震记录有非常多的频率成分。地震波频率的带限性就表明波的传播应该扩展到几何射线周围的有限空间。这一空间范围就成为菲涅尔体。在这片教案中，我们讲介绍关于菲涅尔体的物理理论，展示适用于带限地震波的波动方程的解。波动方程的有限频理论通过敏感核函数精确地描述了带限透射波和反射波的旅行时与振幅和地球介质中慢度扰动之间的线性关系。菲涅尔体和有限频敏感核函数可以通过地震波相长干涉的概念联系起来。波动方程的有限频理论引出了一个反直觉的结论-在三维几何射线上的点状速度扰动不会不会造成波长的相位扰动。因此，这说明在射线理论下的菲涅尔体理论是波动方程有限频理论在有限频下的一个特例。最后，我们还澄清了关于菲涅尔体宽度限制成像实验分辨率的误解。 引言 在地震成像技术中，射线理论通常在正演和反演中被用有构建正反演波长算子。射线理论之所以收到欢迎部分是由于计算机速度和内存的限制，因为射线理论具有较高的计算效率并且对于各种地震成像方法的应用也比较容易。而另一方面，地震成像实验清晰的表明，射线理论，由于他对波场传播的近似描述，对于散射效应严重的波场的成像是不完备的。Cerveny 给出了对于地震波射线理论的一个全面的理解。 在地震成像实验中，记录到的透射波和反射波信号都是由一个主要由低频信号组成的宽带震源激发产生的，因为地震波的高频信号在地层中很容易衰减。但是射线理论是基于高频近似的，这表明基于射线理论的成像技术和和测量波场这件之能会存在方法上的冲突。这个围绕射线且对带限地震波的传播起主要影响的空间范围就被叫做菲涅尔体。射线理论在地下构造尺度大于记录波场的第一菲涅尔带的介质中能够取得较好的效果。对于低频反射波（频率成分在10-70Hz之间）和透射波（频率成分在300-800Hz之间），第一菲涅尔体的宽度可以分别达到500m和50m的量级。这个宽度要大于我们在陆地和海洋的反射波地震勘探以及井间和垂直地震剖面中想要成像的地下地质特征。 在这篇教案中，我们将看到如何将地震分辨率扩展到识别体积小于第一菲涅尔带的不均匀体。我们将展示如把射线理论下的旅行时和振幅公式扩展到更精确的、可以应用与带限反射和透射地震信号波场近似理论。波动方程的有限频理论提出了反射和透射地震波的敏感核函数（也称作Frechet核函数）。这些有限频Frechet核函数将速度扰动和旅行时与振幅的扰动线性的联系起来。有限频波长近似被直接应用到各种地震成

地震波数据生成器SGSw

地震波数据生成器 除了程序提供的30多条实测地震波，一些复杂超限工程在做时程分析时往往需要利用当地安评报告的地震波数据生成自己的时程函数，具体的转换过程是被经常提到的一个问题。 相关命令 工具〉地震波数据生成器... 问题解答 midas提供地震波数据生成器这个专门的工具用于生成自己的时程函数，具体操作步骤如下： 1)打开已安装midas软件的文件夹，找到Dbase文件夹，用记事本打开其中任何 一个后缀为dbs的文件；

2)将安评报告的实测地震波数据完全按上述dbs文件的格式输入后另存，修改后 缀txt为dbs； 3)打开地震波数据生成器，执行菜单操作Generete-Earthquake Record；

4)点击Import，导入第2）步中生成的dbs文件，同时可修改地震波三要素中的 有效峰值和持时，保存为一个sgs文件； 5)midas软件中添加时程函数时，导入第4）步生成的sgs文件即可。 相关知识 时程分析往往作为多遇地震的补充计算手段，规范中要求每条时程曲线计算底部剪力结

果不应小于振型分解反应谱法相应结果的65% ，多条时程曲线计算所得底部剪力结果平均值不应小于振型分解反应谱法计算结果的80%。所以选择合适的波很重要，地震波数据生成器还提供时程函数到反应谱的转换，可以和反应谱分析中地震影响系数曲线进行大致的比较，对结果的正确性给予一定的保证。 具体操作步骤如下： 1）同上。 2）同上。 3）打开地震波数据生成器，执行菜单操作Generete-Earthquake Response Spectra；

4）点击Import，导入第2）步中生成的dbs文件，可选择生成多种形式的反应谱，如绝对加速度、相对速度、相对位移等，保存为sgs文件； 5）和时程函数一样，也可以在定义反应谱函数的时候导入第4）步生成的sgs文件。

新输入地震波

[结构分析] 地震波输入的问题 三向输入简化后的单向输入 首先，将三个方向的地震加速度放到一个文本文件里，如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

第一章地震波运动学

第一章地震波运动学 1.斯奈尔定律与费马原理的关系： 作出各种不同入射角的射线路径(从S 点到D 点)，并计算其相应的旅行时间，作出θ～t(单程)图，从图中找出费马路径，即Tmin 由l 对应的θ；再根据给出的两种介质的速度值，验证这一路径是否符合斯奈尔定律。 2.依据惠更斯原理用做图法证明折射波的出射角等于临界角θ。 3.在0点放炮，在离O 点200米处布置一个排列，有14道，道间距为10米，放一炮后得到的地震记录的一部分如图3—2所示，在该记录上看到的是一个直达波的一组振动图。 请分析这张记录，回答下列问题： (1)读出直达波的到达时间，画出直达波的时距曲线，并根据时距曲线的斜率求出直达波的速度。 (2)根据这张记录，试画出下列各时刻的波剖面，t i =0．1l ；0．13；O ．16；0.17；0．20秒，作图时用一张15×25平方厘米的方格纸，距离x 的比例尺：l 毫米=2米，振幅的比例尺与地震记录上振幅的比例尺相同。 (3)从哪个时刻的波剖面上可以读出这个波的视波长数值来，棍波长等于多少?根据视波长和视周期的公式，从地震记录上得到有关数值，再用公式计算出视波长值，把计算出的值与从波部面上读出的值比较一下。 (4)这个波的波剖面长度是多少?振动图的延续时间是多少? (5)把t=O.16秒时刻的那个完整的波剖面图形与地震记录上的振动图比较一下，能否看出它们之间有什么关系?为什么会有这种关系? 4．已知波速V=1000m ／s ，利用虚爆炸点做下列各图 a)已知反射界面的位置定时距曲线的形状和长度 b)已知时距曲线上t O =1．000秒，极小点坐标t m =0．865秒如图2—5，求反射界面的位置及产状。 5．关于正常时差、倾角时差的计算。 (1)水平界面，均匀覆盖介质，V=2500米／秒，h=1250米，计算炮检距x=0米，100米，200米，……1000米的反射波旅行时t 平。 t = 平 计算各x 值的正常时差：0n t t t ?=-平 (2)倾斜界面，φ=10O ，激发点O 处的界面法浅深度h o =1250米，均匀覆盖层波速V=2500米／秒，计算x=0米，士100米，士200米，……士1000米的反射波旅行时t 斜 t = 斜注意：本题设界面上倾方向与x 的负方向一致，取正号，但x 本身有正负号。 (3)用公式2sin d t V χ? ?= ，计算x=100米，200米，……1000米的倾角时差，并回答此公式计算出的是哪两点间的倾角时差? (4)用公式/ 0d n t t t t ?=-?-斜，计算x=100米，200米，……1000米的倾角时差?并回 答这是相对于哪一点的倾角时差?/ d t ?与d t ?有什么关系?哪一个精确?为什么?

《地震波理论》复习最终版

《地震波理论》复习内容 一、弹性理论基础 1. 柯西公式的意义； 因此弹性体内一点的应力状态可以完全由作用于垂直坐标轴方向的三个截面上的应力向量或其分量所确定。 2. 应力与应变的关系； （为单位函数） 3. 杨氏模量E（纵向应力与纵向应变的比例常数就是材料的弹性模量E, 也叫杨氏模量） 泊松比ν（横向应变与纵向应变之比值称为泊松比，也叫横向变性系数,它是反映材料横向变形的弹性常数）； 4. 拉梅常数λ、μ； 为引入均匀各向同性介质中应力与应变关系，引入λ、μ，μ表示剪切模量。 5. 运动的应力方程和位移方程； 运动应力方程：

运动位移方程： 6. 介质受应力作用产生位移由哪几部分组成； 由式上式可以看出处于应力应变状态下的物体其质点位移由三部分组成: ①平动: u，v，w，这是和参考点M一起作同样的运动，它不使物体形 状改变; ②弹性应变: eij，i,j=x,y,z 这是一种使物体形状和体积发生改变的 运动，称为弹性应变.应变有九个分量，考虑到它的对称性，只有其中六个分量独立的。exx，eyy，ezz称为正应变，exy，eyz，ezx称为切应变; ③旋转: ωx，ωy，ωz这是质点围绕参考点M的旋转运动，不使物体形 状和体积发生改变，不属弹性应变范畴. 7. 导出拉梅方程的前提条件； 在对空间求导时，只有λ、μ不随空间变化，即在均匀介质中才能导出拉梅方程。 8. 能流密度。 表示在单位时间内通过与它垂直的单位截面积的机械能。 二、弹性动力学中的基本波 1. 由拉梅方程导出纵波、横波方程；

拉梅方程 对上式进行散度运算，得到： 对上式进行旋度运算，得到： 2. 平面波、不均匀平面波； 平面波：等相位为平面，且与波的传播方向垂直的波动。 不均匀平面波：平面波传播的方向余弦为l 、m 、n 是复数，这样的波为 不均匀平面波。 3. 在什么情况下才能称为平面波； 离震源较远时可以将在局部等相位内，将点震源产生的球面波看成一个 平面。 4. 沿着x 方向传播的平面波的表示方法； 5. xoz 平面内的波剖面（图2-3）；P 39 若该方程波函数为： ，则波剖面如下图所示 1(,,)()lx nz C f x z t f t + =-''''exp()(,,)exp()exp()Ksh z Kc f x z t A j x j t h ωθθ=--12(x,)()()x x C C f t f t f t =-+ +''' j θθθ=+

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

时程分析时地震波的选取及地震波的反应谱化 发表时间：2017-12-29T15:40:37.810Z 来源：《防护工程》2017年第22期作者：金林飞 [导读] 目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应谱法。 常州市规划设计院江苏常州 213002 摘要：目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应谱法。时程分析法作为补充计算方法，在不规则、重要或较高建筑中采用。进行时程分析时，首先面临正确选择输入的地震加速度时程曲线的问题。时程曲线的选择是否满足规范的要求，则需要首先将时程曲线进行单自由度反应计算，得到其反应谱曲线，并按规范要求和规范反应谱进行对比和取舍。本文通过介绍常用的数值计算方法及计算步骤，实现将地震加速度时程曲线计算转化成反应谱曲线，从而为特定工程在时程分析时地震波的选取提供帮助。 关键词：时程分析，地震波，反应谱，动力计算 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、时程曲线要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间均要符合规定。其中频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定；加速度的有效峰值按《抗规》表5.1.2-2中所列地震加速度最大值采用；输入的地震加速度时程曲线的有效持续时间，一般从首次达到该时程曲线最大峰值的10%那一点算起，到最后一点达到最大峰值的10%为止。有效持续时间一般为结构的基本周期的（5~10）倍，即结构顶点的位移可按基本周期往复（5~10）次。 根据以上《抗规》对地震波选取的要求，在工程中选取地震波时，首先需将地震波进行单自由度动力求解，得出相应的反应谱曲线，并与规范反应谱所用的地震影响系数曲线进行对比。对比频谱特性，查看是否在统计意义上相符，若相符则可将此地震波用于结构计算进行底部剪力验证，否则则需要更换地震波重新分析。 3 求解地震波单自由度体系反应的数值计算方法 对一般动力荷载反应的求解方法可分为叠加法和逐步法两类。 叠加法包括时域分析和频域分析，总反应计算采用独立反应贡献的组合。时域方法中，荷载被考虑为短暂持续时间的脉冲序列，由每个脉冲自由振动反应的独立贡献得到后续时间的总反应；频遇方法中，假设荷载为周期的，并用Fourier变换为离散的谐振分量Pn。再由这些荷载分量乘以结构的频率反应系数Hn，得到与其相应的结构谐振反应分量Vn。最后，由组合谐振反应分量（Fourier逆变换）获得结构的总反应。由于结果都使用了叠加，因此不适用于非线性反应分析。[2] 逐步法有很多种，但所有方法都是将荷载和反应历程分成一系列时间间隔或“步”。在每步期间均以此步开始时存在的初始条件（位移、速度和加速度）和该步期间的荷载历程来计算反应。因此每步反应是一个独立的分析问题。[2] 积分法是逐步法中的一般性的方法之一。积分法对每一时间步，从初始最终条件应用积分向前进一步，速度的变化依赖于加速度历程的积分，而位移的变化依赖于相应的速度积分。 用积分法进行分析时，首先需要假设在时间步内加速度是如何变化的。根据假设的不同，积分法可分为“基于常平均加速度”法（如

怎样在Sap2000中输入地震波

怎样在Sap2000中输入地震波 1、地震波文件类型： （1）、有些地震波数据是一列的，这种情况最好办。直接在定义fuction-time history时定义numbers per line为一即可。但要清楚时间步长。系统自定义的是0.005秒！ （2）、有些地震波是两列的，第一列是时间，第二列是加速度值。这时在定义fuction-time history中有一项time and fuction value一定要选中，而number of point per line是填1而不是2。这样程序先读时间，然后读加速度值。 （3）、还有一些地震波干脆是是5-8列的（8列的在伯克利很常见），其中的每个数据都是加速度值。加速度的读取顺序应该是沿着“行”读，读完本行的数据再读下一行。在sap中，只需设定numbers per line=列数即可。 （4）、地震波的数据为六列，第一、三、五列是时间，二、四、六列是加速度，正确的读数顺序是先横着读完一行再读下一行。导入这种地震波时候这样设置：the number of points per line =3，并选择time and function values。 （5）、从PEER数据库下载的地震波，有的全属于一行，故numbers per line=地震波总点数(NPTS)。 （2）、重力加速度单位gal： 描述地重力场强度的单位，常用于地震工程学中，用来描述地震加速度。gal，称为“伽”或者“盖”，是为纪念第一个重力测量者意大利科学家伽利略（1564~1642）而命名的。重力场的量纲是厘米每二次方秒，规定1 cm/s^2=1 gal，1 g（重力加速度）= 9.8 m/s2 = 980 gal。我国各类抗震规范中，也用gal来描述地震加速度峰值。 （3）、定义时程函数时，单位无所谓，只要你的系数对应好就可以。sap输入的地震函数本身是没有单位的，它的单位随着你sap的右下角的单位走的。所以才需要将这个单位和原始波单位对应。 3、时程分析输出步长与输出点数： 如果是采用sap2000中自带的地震波时，在导入地震波后，通过点击显示文件（View File）可以显示该地震波的总时间t和时间步长s（输出步长大小），通过这两个数据可以计算输出时段数( )。对于同一个工程，地震波作用的总时间t可以不同，但输出时段大小s应该一致。否则，两次时程分析输出的是不同时间点处的结果，其值当然是不同的。一般来说，输出时段大小一般取0.01或0.02。输入的是同一个地震波,只不过输出的时段数n和时段大小s不