基于MATLAB的地震反应谱计算方法的比较

地震响应的反应谱法与时程分析比较 (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

地震反应谱分析实例

结构地震反应谱分析实例 在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑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

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法为如下所示：

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

反应谱理论与人工模拟地震波技术简介

第３３卷第２６期?１０６?２００７年９月山西建筑 ＳＨＡＮＸＩＡＲＣＨＩＴＥＣＴＵＲＥ Ｖｄ３３Ｎｏ．２６ Ｓｅｐ．２００７ 文章编号：１００９—６８２５｛２００７）２６—０１０６—０３ 反应谱理论与人工模拟地震波技术简介 邱玉国王玉富 摘要：介绍了反应谱理论的发展历程和国内外研究现状，分析了研究问题的思路，指出了利用反应谱理论来解决实际工程时遇到的问题，并简单介绍了国外对人工模拟地震波技术的应用和研究，为抗震理论提供了参考依据。 关键词：反应谱理论，地震波，随机振动，非弹性地震波 中图分类号：ＴＵ３５２文献标识码：Ａ １概述 反应谱理论是建筑结构抗震设计的重要理论基础之一。从２０世纪５０年代开始，反应谱理论逐渐成为结构抗震设计的重要方法，经过５０多年的发展，目前这种方法已经为世界上大多数国家的设计规范所采用。但是，由于地震产生机理和作用效果的复杂性，采用反应谱理论进行分析和设计与工程实践还存在很多与实际不相符合之处。此外，对于反应地震重要特性的时间问题，反应谱法也无能为力。 人工模拟地震波技术是近年来才发展起来的一项新的结构抗震设计的技术手段，目前主要用于计算机模拟和特别重要结构模型的振动台试验。它能够通过模拟地震波的特性来用于对结构进行时程分析，是～种新兴的、具有革命性意义的试验手段。 图２数值模拟结果２．３计算结果分析 通过数值模拟和试验得到瓦斯管承载力等数值如表２所示。 表２数值模拟和试验结果 Ｉ研究方法承载力仆但ａ最大应变／％最大剪应力／ＳＰａＩ数值模拟７．１４Ｏ．０８４２１６０室内试验６．６２０．０９６４ ３结语 通过对丁集煤矿瓦斯管材质和整体抗外压的试验研究以及数值模拟分析，可以获得如下重要结论： １）通过对管材材质的试验研究表明：工作管材质采用Ｑ３４５，尺寸为柘３０ｒｆｌｌＴｌ×１４ｉｎｌｎ，能够满足强度和稳定性要求。 ２）瓦斯管整体抗外压试验结果表明：工作管抗外压承载力为６，６２ＭＰａ；通过大变形有限元数值计算，采用变形稳定性控制其承载力，结果为７．１４ＭＰａ，两者数值十分接近，说明用文中方法模拟大直径瓦斯管的承载力是可行的。 参考文献： ［１］李正来．瓦斯抽排钻孔定向技术的改进［Ｊ］．安徽科技，２００６（３）：４９—５０． ［２］汪东生．瓦斯抽排技术治理本煤层采空区瓦斯涌出的实践［Ｊ］．煤矿安全，２００６（１）：１３—１５． ［３］张敦伍，任胜杰．瓦斯抽排钻孔防偏斜实践［Ｊ］．矿业安全与环保，２００５（８）：６７—６８． ［４］刘克功，范再良，赵新华．采空区瓦斯抽排法治理综放面瓦斯超限［Ｊ］．煤，１９９８（２）：４８—５０． Ｓｔｕｄｙｉｎｇｏｎｒａｄｉａｌｓｔａｂｉｌｉｔｙｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｌａｒｇｅｐｉｐｅｉｎｍｉｎｅ ＴＯＮＧＷｅｎ－ｌｉｎ Ａｂｓｔｒａｃｔ：ＴｈｅｅｘｐｅｒｉｍｅｎｔａｌａｎｄｖａｌｕｅｓｉｍｕｌａｔｉｏｎｍｅｔｈｏｄｓｈａｖｅｓｔｕｄｉｅｄｔｈｅＤｉｎｇｉｉｃｏａｌｍｉｎｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅｕｎｄｅｒｍｅｃｈａｎｉｃｓｃｈａｒａｃｔｅｒ—ｉｓｔｉｅ．Ｒｅｓｕｌｔｉｎｄｉｃａｔｅｄ：ｔｈｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅｉｓｐｒｅｓｅｎｔｅｄｓｔａｂｉｌｉｔｙｆａｉｌｕｒｅｍｏｄｅｌｉｎｅｎｃｉｒｃｌｅｓｐｒｅｓｓｅｓｓｈａｐｅ，ｉｔｓｓａｆｅｔｙｆａｃｔｏｒｒｅａｃｈｅｓ３．０，ｉｔｉｓｄｅｓｉｇｎｔｈｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅａｎｄｔｈｅｃｏｎｓｔｒｕｃｔｐｍｖｉｄｅｓｔｈｅｒｅｆｅｒｅｎｃｅ． Ｋｅｙｗｏｒｄｓ：ｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅ，ｅｘｐｅｒｉｍｅｎｔａｌｉｎｌａｂ，ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ，ｓｔａｂｉｌｉｔｙｆａｉｌｕｒｅｍｏｄｅｌ 收稿日期：２００７．０４．０６ 作者简介：邱玉国（１９７３。），男，工程师，辽宁工程技术大学软件学院，辽宁阜新１２３０００ 王玉富（１９７０．），男，工程师，中铁十九局集团第三工程有限公司，辽宁辽阳１１１０００

人工地震波生成程序简介

姓名：郭勇 学号：022******* 人工地震波生成程序简介 一、程序设计内容及方法 1、程序内容 本程序根据特征周期、水平地震波影响系数最大值和地震波幅值等初始条件生成人工地震波，为结构动力分析的时程分析法提供地震波来源。 2、程序设计方法 (1) 理论依据 本程序采用三角级数法生成人工地震波。 对于给定的功率谱密度函数，按照下面的公式可以方便的生成以为功率谱密度函数、均值为零的高斯平稳过程。 (1) 式中： (2) 为内均匀分布的随机相角；，分别为正域内的上、下限值，即认为的有效功率在范围内，而范围外的值可视为零。 为了反映地面运动的非平稳性，采用包络函数乘以平稳过程， (3) (3)式即为人工地震波模型。 可根据下式确定： (4) 式中：为衰减系数，通常取值范围为0.1~1.0，本程序取0.15；，和根据不同实际情况取值，为地震波持时，本程序取，分别为4s，15s，和均为40s。 本程序采用《建筑抗震设计规范》(GB50011-2001)中的反应谱作为目标谱，通过Kaul 提出的平稳过程反应谱与功率谱的近似关系 (5) 式中：为规范反应谱；为阻尼比；为地震动持时；为反应不超过反应谱值的概率，本程序取0.85。通过(3)式和(5)式即可生成人工地震波。 (2) 程序实现方法 首先建立基于对话框的应用程序框架，添加的主要控件为3个编辑框和4个按钮。3个编辑框分别作为程序中的特征周期(对应成员变量为m_dTg)、水平地震影响系数最大值(对应成员变量为m_dAmax)和地震波幅值(对应成员变量为m_pd)3个数据的交互输入处；4个按钮分别为"生成地震波"、"输出地震波"、"输入地震波"和"退出"。 添加的成员函数有：Wavegener()(生成地震波)、Wavedrawing()(绘制地震波加速度时程曲线)、OnSTART()(对应"生成地震波"按钮，实现生成地震波的功能)、OnOutput()(对应"输出地震波"按钮，实现输出数字化的地震波记录的功能)和OnInput(对应"输入地震波"按钮，实现输入数字化的地震波记录并绘制其加速度时程曲线的功能)。 几点说明： a 生成随机相角的程序如下： srand((unsigned)time( NULL ));

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

三设计地震动反应谱确定的规范方法 设计地震动是通过对地震环境和场地环境的分析判断和分类方法确定。工程勘察单位至少提供： 设计基本地震加速度和设计特征周期 场地环境：覆盖层厚度、剪切波速、土层钻孔资料 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

反应谱生成人工地震波

反应谱生成人工地震波 一、软件SIMQKE_GR使用说明 1.先安装程序 2.使用方法 双击，打开程序，可以得到如图1界面。 图1 程序开始界面 如图1所示，由于程序本身提供的反应谱是适用于欧洲规范的，不适合于我国的规范反应谱，因此不能通过调整参数来获得符合我国规范的反应谱。可以采用导入的方法来输入反应谱。 3.点击菜单栏“file”—“Import spectra data”，出现打开对话框，如图2所示， 要求打开一个已经存在的反应谱文件（如 1.srf）。

图2 导入反应谱文件对话框 4.文件格式如下所示（红字部分不能修改，注意反应谱单位为g），下面部分 可以替换。 response spectrum time(s) acc(g) 0 0.1215 0.01 0.13635 0.02 0.1512 0.03 0.16605 0.04 0.1809 0.05 0.19575 0.06 0.2106 0.07 0.22545 0.08 0.2403 0.09 0.25515 0.1 0.27 0.15 0.27 0.2 0.27 0.25 0.27 0.3 0.27 0.35 0.27 0.4 0.27 0.45 0.27

0.5 0.243 0.6 0.2025 0.7 0.173571429 0.8 0.151875 0.9 0.135 1 0.1215 1.1 0.110454545 1.2 0.10125 1.3 0.093461538 1.4 0.086785714 1.5 0.081 1.6 0.0759375 1.7 0.071470588 1.8 0.0675 1.9 0.063947368 2 0.06075 2.1 0.057857143 2.2 0.055227273 2.3 0.052826087 2.4 0.050625 2.5 0.0486 2.6 0.046730769 2.7 0.045 2.8 0.043392857 2.9 0.041896552 3 0.0405 3.1 0.039193548 3.2 0.03796875 3.3 0.036818182 3.4 0.035735294 3.5 0.034714286 3.6 0.03375 3.7 0.032837838 3.8 0.031973684 3.9 0.031153846 4 0.030375 4.1 0.029634146 4.2 0.028928571 4.3 0.028255814 4.4 0.027613636 4.5 0.027 4.6 0.026413043 4.7 0.025851064 4.8 0.0253125

地震反应谱的绘制

地震时程曲线与反应谱的绘制 ①地震反应谱的意义 地震反应谱表示的是在一定的地震动下结构的最大反应，是结构进行抗震分析与设计的重要工具。 由于同一结构在遭遇不同的地震作用时的反应并不相同，单独一个地震记录的反应谱不能用于结构设计。但是地震记录的反应谱又有一定的相似性，我们可以将具有普遍特性记录的反应谱进行平均和平滑处理，以用于抗震设计。现在，地震反应谱不但是工程抗震学中最重要的概念之一，还是整个地震工程学中最重要的概念之一。 ②地震反应谱的计算方法 反应谱的计算方法涉及到时域分析方法和频域分析方法。 时域分析方法中的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')

人工地震动生成程序

clear clc close all hidden %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% fni=input('生成人工地震波-输入数据文件名(20041012):','s'); fid=fopen(fni,'r'); fs=fscanf(fid,'%f',1);%采样频率 tu=fscanf(fid,'%f',1);%上升时间长度 %上升时间包络线线形(1-直线、2-抛物线、3-指数曲线) iu=fscanf(fid,'%f',1); %上升时间包络线线形参数（只有指数曲线需要具体参数，其均为1） cu=fscanf(fid,'%f',1); ta=fscanf(fid,'%f',1);%持时时间长度 td=fscanf(fid,'%f',1);%下降时间长度 %下降时间包络线线形(1-直线、2-抛物线、3-指数曲线) id= fscanf(fid,'%f',1); %下降时间包络线线形（只有抛物线，指数曲线需要具体参数，其余为1） cd=fscanf(fid,'%f',1); dp=fscanf(fid,'%f',1);%阴尼比值 p=fscanf(fid,'%f',1);%概率系数（一般可取P=0.85） nn=fscanf(fid,'%f',1);%迭代次数 fno=fscanf(fid,'%f',1);%输出数据文件名

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %对目标反应谱取值 x=fscanf(fid,'%f',[2,inf]);%反应谱频率和幅值数据 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% tatus=fclose(fid); %计算生成地震波的数据长度 tl=tu+ta+td; %计算生成地震波的数据长度 nt=round(fs*tl+1); %大于并最接近nt的2的幂次方为FFT长度 nfft=2^nestpow2(nt) %计算频率间隔（Hz） df=fs/nfft %定义反应谱的离散频率向量 f=0:df:(nfft/2-1)*df %计算时间间隔（s） dt=1/fs; %定义的离散时间向量 t=0:dt:(nt-1)*dt %生成0到2PI的随机数为随机相位 g=rand(1,nfft/s)*2*pi; %建立时间包络线 %建立与地震波长度相同元素为1的向量 en=ones(1,nt); %上升时间阶段 %确定上升时间段的长度 l=round(tu*fs)+1 %产生上升时间段的包络线数组元素 switch iu case 1 %直线 en(1:l)=linspace(0,1,1);% y = linspace(a,b,n) generates a row vector y of n points linearly

上海地震波-三向输入选取(说明)m

上海地区抗震设计输入地震时程说明 （共8页） 同济大学房结构工程与防灾研究所 二〇一二年六月

目录 1 天然地震时程选取原则 (3) 2 峰值调整 (3) 3 频谱特性 (3) 4 地震动持时 (3) 5 人造地震动生成的方法 (3) 6 目标反应谱的确定 (4) 7 所选地震时程的基本信息 (4) 8 地震时程反应谱与规范反应谱对比 (5)

上海地区抗震设计输入地震时程说明 1 天然地震时程选取原则 天然地震动具有很强的随机性，随着输入地震波的不同结构的地震响应也会有很大的差异，故要保证时程分析结果的合理性，在选择地震波时必须遵循一定的原则。一般而言，选择输入地震波时应以地震波的三要素（峰值、频谱特性、地震动持时）为主要考虑因素。 2 峰值调整 地震波的峰值一定程度上反应了地震波的强度，因此要求输入结构的地震波峰值应与设防烈度要求的多遇地震或罕遇地震的峰值相当（峰值相当并非峰值相等，而是在峰值相近的情况下所选用地震波的反应谱与规范反应谱基本相符）。 3 频谱特性 频谱是地面运动的频率成分及各频率的影响程度。它与地震传播距离、区域、介质及结构所在的场地土性质有密切关系。一般来说，在震中附近或岩石等坚硬场地土中，地震波中的短周期成分较多，在震中距较远或软弱场地土中，地震波的长期成分较多。输入地震波的卓越周期应尽可能与拟建场地的特征周期一致，且在一定的周期段内与规范反应谱尽量接近。对于天然地震记录而言，3个方向地震波同时都与规范反应谱很接近的条件是很难满足的，但应保证至少一个水平向地震波反应谱与规范反应谱基本吻合。 4 地震动持时 地震持时也是结构破坏和倒塌的重要因素，工程实践中确定地震动持续时间的原则是：1）地震记录最强烈部分应包含在所选持续时间内，2）若对结构进行弹塑性地震反应分析（考虑累计损伤效应），持续时间可取长些。另外，在截取地震波时尚需注意尽量在速度/位移零点处截断以尽量避免加速度积分时速度或位移的``漂移''现象。 5 人造地震动生成的方法 工程中较为实用的人造地震动的生成方法主要有两种：一是将不同频率具有随机相位的三角波进行叠加并根据目标反应谱或功率谱进行迭代修正；二是选择满足场地条件等要求的天然地震记录，保留其相位等随机特征，然后修正其不同频段的幅值以逼近目标反应谱或功率谱。由于三个方向地震动间的随机相关性关系很复杂，采用三角波叠加的方法生成地震波时各方向地震波的相关性难以确定，故本文采

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

正确选取地震波

地震波的选取方法（MIDAS (2009-05-16 22:51:32) 转载▼ 分类：结构专业 标签： 杂谈 建筑抗震设计规范(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组尼比加速度反应谱的平均放大系数。 上述方法使用的是将频段固定的方法来求EPA和EPV，1990年的《中国地震烈度区划图》采用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是：在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱，找出加速度反应谱平台段的起始周期T0和结束周期T1，然后在拟速度反应谱上选定平台段，其起始周期为T1(即加速度反应谱平台段的结束周期T1)，结束周期为T2，将加速度反应谱在T0至T1之间的谱值求平均得Sa，拟速度反应谱在T1至T2之间的谱值求平均得Sv，加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5，按公式(1)、（2）、（3）求得EPA、EPV、Tg。 在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功

地震反应谱的特性

地震反应谱的特性 崔济东(JiDong Cui) (华南理工大学土木与交通学院，广东广州，510640) 1反应谱的基本概念（Introduction to Response Spectra） 地震动反应谱：单自由度弹性系统对于某个实际地震加速度的最大反应（可以是加速度、速度和位移）和体系的自振特征（自振周期或频率和阻尼比）之间的关系。前一篇博文《Earthquake Response Spectra地震反应谱》介绍了反应谱和伪反应谱的基本概念，并编制了相应的反应谱计算程序——SPECTR。本文利用该软件，通过几个实测地震记录的反应谱分析，总结地震反应的一般谱特性。 2本文用到的地震加速度记录（Acceleration Time History Records） 2.11999年台湾集集地震记录的加速度记录： （1）加速度记录信息： The Chi-Chi (Taiwan) earthquake of September 20, 1999. Source: PEER Strong Motion database Recording station: TCU045 Frequency range: 0.02-50.0 Hz Maximum Absolute Acceleration: 0.361g （2）加速度时程与相应的速度和位移

图2-1 ChiChi地震加速度时程2.21994年美国北岭地震记录的加速度时程： （1）加速度记录信息： The Northridge (USA) earthquake of January 17, 1994. Source: PEER Strong Motion Database Recording station: 090 CDMG STATION 24278 Frequency range: 0.12-23.0 Hz Maximum Absolute Acceleration: 0.5683g （2）加速度时程与相应的速度和位移 作者：崔济东（1988- ），男，结构工程专业，博士研究生。

地震波描述

1.1设计加速度过程线 依据GB18306-2001《中国地震动参数区划图》，50年超越概率为10%时，工程区地震动峰值加速度为0.15 g，地震动反应谱特征周期为0.45 s，相应地震基本烈度为7度。场地土属中软场地土，场地类别为Ⅱ类。根据（DL5073-2000）《水工建筑物抗震设计规范》的规定，本工程壅水建筑物抗震设防类别为乙类，设计烈度按7度取。 参考工程地质报告，本课题选取美国Taft地震波、人工地震波与实测地震波共三条地震波进行分析。 Taft地震波，1952年7月21日发生于美国的加利弗里亚州地震（California Earthquake，震级7.4级），是位于加州Kern County林肯学校的No.1095地震台测得的地震记录，该记录地距震中约43.5 km。地震仪设于学校附近一隧洞混凝土地板上，测得完整的三向地震波，记录长达54 s，最大地震加速度175.9 cm/s2，最大速度17.7 cm/s，最大位移9.15 cm。Taft地震波由于记录完整、数据可靠，在国际地震工程界被广泛引用。本报告中将其峰值加速度调整至0.15 g得到设计地震加速度过程线进行动力反应分析，通过SHAKE91程序反演后，坝基水平向基岩地震波峰值为0.12 g，竖直向基岩地震波峰值为0.08 g。横河向、顺河向和竖直向输入加速度之比为3:3:2。计算地震时长20 s，时间步长为0.02 s，各方向地震波时程如图1.1-1至图1.1-3所示。 人工地震波，是根据《水工建筑物抗震设计规范》选取规范标准反应谱为目标谱生成。人工波生成时，迭代误差取为5%，其中特征周期T g按照基岩场地取0.3 s，反应谱最大值的代表值βmax取为2，设计加速度代表值为0.15 g。由此得到设计地震加速度过程线进行动力反应分析，通过SHAKE91程序反演后，坝基水平向基岩地震波峰值为0.12 g，竖直向基岩地震波峰值为0.08 g。横河向、顺河向和竖直向输入加速度之比为3:3:2。计算地震时长20 s，时间步长为0.02 s，各方向地震时程如图1.1-4至图1.1-6所示。

结构地震反应谱分析实例

在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑X 与Y向地震效应作用。已知地震影响系数a与周期T的关系： a(T)= 0.4853*(0.4444+2.2222*T) 0<T<=0.04 秒 0.4853*(0.10/T)^(-0.686) 0.04<T<=0.1 秒 0.4853 0.1<T<=1.2 秒 0.4853*(1.2/T)^1.5 1.2<T<=4 秒 以下是命令流程序 ---------------------------------------------------------------------------------------------------- /filname,SPEC,1 /PREP7 !定义单元类型及材料特性 ET,1,45 MP,EX,1,2.8E10 MP,DENS,1,2.4E3 MP,NUXY,1,0.18 !建立模型 BLOCK,0,1,0,1,0,5 !网格剖分 ESIZE,0.5 VMESH,all /VIEW,,-0.3,-1,1 EPLOT FINISH /SOLU !施加底部约束 ASEL,,LOC,Z,0 DA,ALL,ALL ALLSEL !施加自重荷载 ACEL,0,0,10 !进行模态求解

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 lcwrite,11