爆破地震波作用下建筑物的动力响应分析方法
在 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.1 机构运动分析的任务、目的和方法 曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构，它可以用来实现转动和摆动之间运动形式的转换或传递动力。 对四杆机构进行运动分析的意义是：在机构尺寸参数已知的情况下，假定主动件（曲柄）做匀速转动，撇开力的作用，仅从运动几何关系上分析从动件（连杆、摇杆）的角位移、角速度、角加速度等运动参数的变化情况。还可以根据机构闭环矢量方程计算从动件的位移偏差。上述这些内容，无论是设计新的机械，还是为了了解现有机械的运动性能，都是十分必要的，而且它还是研究机械运动性能和动力性能提供必要的依据。 机构运动分析的方法很多，主要有图解法和解析法。当需要简捷直观地了解机构的某个或某几个位置的运动特性时，采用图解法比较方便，而且精度也能满足实际问题的要求。而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时，采用解析法并借助计算机，不仅可获得很高的计算精度及一系列位置的分析结果，并能绘制机构相应的运动线图，同时还可以把机构分析和机构综合问题联系起来，以便于机构的优化设计。 1.2 机构的工作原理 在平面四杆机构中，其具有曲柄的条件为： a.各杆的长度应满足杆长条件，即： 最短杆长度+最长杆长度≤其余两杆长度之和。 b.组成该周转副的两杆中必有一杆为最短杆，且其最短杆为连架杆或机架（当最短杆为连架杆时，四杆机构为曲柄摇杆机构；当最短杆为机架时，则为双曲柄机构）。 第一组(2代一套)四杆机构L1=125.36mm,L2=73.4mm,L3=103.4mm,L4=103.52mm 最短杆长度+最长杆长度(125.36+73.4) ≤其余两杆长度之和(103.4+103.52) 最短杆为连架杆，四杆机构为曲柄摇杆机构 第二组(2代二套)四杆机构L1=125.36mm,L2=50.1mm,L3=109.8mm,L4=72.85mm 最短杆长度+最长杆长度(125.36+50.1) ≤其余两杆长度之和(109.8+72.85) 最短杆为连架杆，四杆机构为曲柄摇杆机构 第三组(3代)四杆机构L1=163.2mm,L2=61.6mm,L3=150mm,L4=90mm 最短杆长度+最长杆长度(163.2+61.6) ≤其余两杆长度之和(150+90) 最短杆为连架杆，四杆机构为曲柄摇杆机构 在如下图1所示的曲柄摇杆机构中，构件AB为曲柄，则B点应能通过曲柄与连杆两次共线的位置。 1.3 机构的数学模型的建立 图1机构结构简图 在用矢量法建立机构的位置方程时，需将构件用矢量来表示，并作出机构的封闭矢量多边形。如图1所示，先建立一直角坐标系。设各构件的长度分别为L1 、L2 、L3 、L4 ， 其方位角为、、、。以各杆矢量组成一个封闭矢量多边形，即ABCDA。其个矢量之和必等于零。即：
地震波运动学理论
第二章地震波运动学理论 一、名词解释 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.纵测线和非纵测线：激发点与接收点在同一条直线上，这样的测线称为纵测线。用纵测线进行观测得到的时距曲线称为纵时距曲线。激发点不在测线上，用非纵测线进行观测得到的时距曲线称为非纵时距曲线。
地震波的概念种类特点及地表影响
关于地震波 摘要：地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面，将地球内部分为地壳、地幔和地核三个圈层。 关键词：地震波辐射地球内部 一：背景 ①2008年5月12日14时28分04秒，四川汶川、北川，8级强震猝然袭来，大地颤抖， 山河移位，满目疮痍，生离死别……西南处，国有殇。这是新中国成立以来破坏性最强、波及范围最大的一次地震。此次地震重创约50万平方公里的中国大地！为表达全国各族人民对四川汶川大地震遇难同胞的深切哀悼，国务院决定，2008年5月19日至21日为全国哀悼日。自2009年起，每年5月12日为全国防灾减灾日。 ②1976年7月28日北京时间03时42分53.8秒，在中国河北省唐山、丰南一带（东经118.2°， 北纬39.6°）发生了强度里氏7.8级（矩震级7.5级），震中烈度Ⅺ度，震源深度23千米的地震。地震持续约12秒。有感范围广达14个省、市、自治区，其中北京市和天津市受到严重波及。强震产生的能量相当于400颗广岛原子弹爆炸。整个唐山市顷刻间夷为平地，全市交通、通讯、供水、供电中断。唐山地震没有小规模前震，而且发生于凌晨人们熟睡之时，使得绝大部分人毫无防备，造成24.2万人死亡，重伤16.4万人，名列20世纪世界地震史死亡人数第一。 ③邢台地震由两个大地震组成：1966年3月8日5时29分14秒，河北省邢台专区隆尧县 （北纬37度21分，东经114度55分）发生震级为6.8级的大地震，震中烈度9度强； 1966年3月22日16时19分46秒，河北省邢台专区宁晋县（北纬37度32分，东经115度03分）发生震级为7.2级的大地震，震中烈度10度。两次地震共死亡8064人，伤38000人，经济损失10亿元。这是一次久旱之后的大震。
时程分析中地震波输入位置的讨论
时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[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，3或4），选择时程曲线。在定义时程工况时，对于多遇或罕遇地震，按比例调整时程曲线的最大值。中国抗震规范规定，作为抗震计算中底部剪力法和振型分解反应谱法的补充方法，对于特别不规则，特别重要的和较高的结构应采用时程分析法进行多遇地震下的补充计算。 可取多条时程曲线的计算结果的平均值与振型分解反应谱法计算结果的较大值。 采用时程分析法时，应咱建筑场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。 其加速度时程最大值可按规范中对于多遇和罕遇地震在不同烈度下的值。 弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65％，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80％ 。 可使用弹塑性时程分析法计算罕遇地震下结构的变形。 时程分析是一个承受随时间变化的指定荷载结构的逐步动态反应分析，可以是线性或非线性的。 此章对时程分析进行一般的描述，特别是线性时程分析。 定义时程函数 用户可使用“从文件中添加函数”，导入已定义的文本文件，即实测的时程曲线；也可使用程序内置的时程函数。
时程函数定义对话框 时程函数定义对话框中的条目解释如下： ?函数名 通过在编辑框中直接键入以指定或修改时程函数的名称。 ?函数文件 1.在函数文件域点击浏览按钮以调出一个对话框，在此可找出包含时程函数的 文本文件名。注意文件名显示在文件名框中 2.在 "要跳过的标题行" 编辑框中输入一个希望ETABS在文本文件中跳过的 行数。 3.在 "每行要跳过的前缀字符" 编辑框中输入一个希望ETABS在文本文件中 每行要跳过的字符数。 4.在 "每行的点数" 编辑框中输入一个数告诉ETABS文本文件每行的绘图点 数。
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 ………………
动力学分析方法
1 动力学分析方法 结构动力学的研究方法可分为分析方法（结构动力分析）和试验方法（结构动力试验）两大类。[7-10] 分析方法的主要任务是建模（modeling），建模的过程是对问题的去粗取精、去伪存真的过程。在结构动力学中，着重研究力学模型（物理模型）和数学模型。建模方法很多，一般可分为正问题建模方法和反问题建模方法。正问题建模方法所建立的模型称为分析模型（或机理模型）。因为在正问题中，对所研究的结构（系统）有足够的了解，这种系统成为白箱系统。我们可以把一个实际系统分为若干个元素或元件（element），对每个元素或元件直接应用力学原理建立方程（如平衡方程、本构方程、汉密尔顿原理等），再考虑几何约束条件综合建立系统的数学模型。如果所取的元素是一无限小的单元，则建立的是连续模型；如果是有限的单元或元件，则建立的是离散模型。这是传统的建模方法，也称为理论建模方法。反问题建模方法适用于对系统了解（称黑箱系统——black box system）或不完全了解（称灰箱系统——grey box system）的情况，它必须对系统进行动力学实验，利用系统的输入（载荷）和输出（响应——response）数据，然后根据一定的准则建立系统的数学模型，这种方法称为试验建模方法，所建立的模型称为统计模型。 在动力平衡方程中，为了方便起见一般将惯性力一项隔离出来，单独列出，因此通常表达式为： +P M (2) u I - = 其中M为质量矩阵，通常是一个不随时间改变的产量；I和P是与位移和速度有关的向量，而与对时间的更高阶导数无关。因此系统是一个关于时间二级导数的平衡系统，而阻尼和耗能的影响将在I和P中体现。可以定义： + = (3) I Ku C u 如果其中的刚度矩阵K和阻尼矩阵C为常数，系统的求解将是一个线性的问题；否则将需要求解非线性系统。可见线性动力问题的前提是假设I是与节点位移和速度是线性相关的。 将公式(2)代入(1)中，则有 (4) + M= + u P Ku C u
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
模拟地震波传播可视化
模拟地震波传播的可视化研究 摘要：实验中选取了与地壳平均波速相近的光学玻璃作样品，利用动态光弹的成像系统，来观测波在光学玻璃及波从光学玻璃透射到水中的传播过程，并记录0～50μs内的波的传播过程，以此来模拟地震波在地壳中反射、透射等传播行为。 abstract： the experiments selected optical glass similar with average velocity， and used imaging system of dynamic photoelasticity to observe the communication process of light in optical glass and light refraction from optical glass to water， and record the wave transmission during 0～50μs，for simulating reflection and transmission of seismic wave in crust. 关键词：地震波；动态光弹；反射；透射 key words： seismic waves；dynamic photoelastic；reflection；transmission 中图分类号：p315.3+1 文献标识码：a 文章编号：1006-4311（2013）04-0297-02 0 引言 在地震勘探中，通常是通过检波器来记录地下地震波带来的信息，根据相应的数学和物理模型进行复杂的计算机处理以获得地下的构造情况，虽然地震勘探的相关理论有很大的发展，但是由于理论结果难以获得，并且对于复杂形状的结构，解析方法变得相当繁
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核函数将速度扰动和旅行时与振幅的扰动线性的联系起来。有限频波长近似被直接应用到各种地震成
新输入地震波
[结构分析] 地震波输入的问题 三向输入简化后的单向输入 首先，将三个方向的地震加速度放到一个文本文件里，如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
空间机构动力学分析方法的研究
空间机构动力学分析方法的研究 RESEARCH OF SPATIAL MECHANISM DYNAMICS ANALYSIS METHODS 袁清珂 1 刘大慧1 惠延波2 张明天1 成思源 1 (1.广东工业大学机电工程学院,广州510006) (2.河南工业大学机电工程学院,郑州450007) YU AN QingKe 1 LIU DaHui 1 HUI YanBo 2 ZHANG MingTian 1 CHENG SiYuan 1 (1.College o f Mechanical &Electrical Engineering ,Guangdong University o f Technology ,Guangzhou 510006,China ) (2.College o f Mechanical Engineering ,Henan University o f Technology ,Zhengzhou 450007,China ) 摘要 在空间机构中约束运动副处的构件上建立笛卡尔直角坐标系,开发描述空间机构结构形态的符号体系,讨论杆件形状矩阵和约束运动矩阵,在运动链上连续使用变换矩阵,建立空间机构的运动方程,分析空间机构二种类型的动力学建模与分析方法,在此基础上开发空间机构动力学通用分析软件,并给出应用实例。 关键词 空间机构 机械学 运动学 动力学 数值方法 软件工程中图分类号 TH112 TH113 Abstract Cartesian coordinate sys tems were built on the two links at a cons train t kinematics pair respectively,so a notation set to describing spatial mechanis ms was established.Link shape matrices and constraint motion matrices were discussed.By in series usin g these transform matrices on a kinematics chain,the motion equations of spatial mechanisms were set up.The modeling and analyzin g methods for two kinds of spatial mechanis m dynamics problems were explored.Based on the above research,the general analysis software spatial mechanisms dynamics were developed,an example was given. Key words Spatial mechanism ;Mechanics ;Kinem atics ;Dynamics ;Numerical methods ;Software eng ineering Correspon ding author :YUAN Qing Ke ,E mail :qkyuan @gdut .edu .cn The project supported by the National Natural Science Foundation of Chi na (No.50805025),and Guangdong Technological Plans Projects (No.2008B010400011),and Guangzhou T echnological Plans Projects (No.2008Z1 D371). Manuscript received 20080821,in revi sed form 20090307. 引言 研制开发通用机构计算机自动分析软件,首先遇到的问题是如何以一种适当的计算机能够理解的 方式来描述机构的结构形态,使计算机能够自动识别机构、自动建立机构的运动方程、自动求解运动方程,并以数字和图形的方式输出结果。目前,常用的方法是基于机构分组的方法,通过数据文件表达机构的结构形态,这种方法存在描述机构范围有限和用户使用不便等不足。要实现机构分析软件的真正通用化和自动化,必须建立描述机构的通用方法和语言。通过通用方法和专用语言描述各种机构,计算机能识别这种描述,并且能通过这种描述自动识别机构的结构形态和运动链,自动建立机构运动方程,自动进行求解并输出结果。 本文在Denavit 和Hartenberg 提出的用于描述低副机构的描述方法(Denavit Hartenberg,D H )[1 2] 的基础 上,结合有关分析方法 [3 9] ,提出一种新的空间机构运 动建模与分析方法。在空间机构中约束运动副处的构件上建立笛卡尔直角坐标系,开发描述空间机构结构形态的符号体系,讨论杆件形状矩阵和约束运动矩阵,在运动链上连续变换矩阵,建立空间机构的运动方程,分析空间机构两种类型的动力学建模与分析方法,在此基础上,开发通用分析软件,并给出应用实例。 1 运动方程的建立 1.1 坐标系的建立 在空间机构中每一运动副处,分别在构成该运动副的两个构件上,根据运动副的性质和特征,按照不同的规律建立固结于构件上的直角坐标系, 在机构运动 Journal of Mechanical Strength 2011,33(1):040 044 袁清珂,男,1963年1月生,山东青岛人,汉族。广东工业大学教授,从事知识工程与智能设计、机电控制、多体动力学与计算机仿真、企业 信息化、电子商务与网络化制造的研究。 20080821收到初稿,20090307收到修改稿。国家自然科学基金(50805025)、广东省教育部产学研结合项目(2009B090300340)、广东省科技计划 (2008B010400011)、广州市科技计划(2008Z1 D371)资助。
时程分析时地震波的选取及地震波的反应谱化
时程分析时地震波的选取及地震波的反应谱化 发表时间：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] 积分法是逐步法中的一般性的方法之一。积分法对每一时间步，从初始最终条件应用积分向前进一步，速度的变化依赖于加速度历程的积分，而位移的变化依赖于相应的速度积分。 用积分法进行分析时，首先需要假设在时间步内加速度是如何变化的。根据假设的不同，积分法可分为“基于常平均加速度”法（如