地震反应谱分析实例
地震响应的反应谱法与时程分析比较
地震响应的反应谱法与时程分析比较(总13页)-CAL-FENGHAI.-(YICAI)-Company One1-CAL-本页仅作为文档封面，使用请直接删除发电厂房墙体地震响应的反应谱法与时程分析比较1问题描述发电厂房墙体的基本模型如图1所示：图1 发电厂墙体几何模型基本要求：依据class 的最后一页的作业建立ansys模型,考虑两个水平向地震波的共同作用(地震载荷按标准谱缩放,谱值如下)，主要计算底部跨中单宽上的剪力与弯矩最大值,及顶部水平位移。
要求详细的ansys反应谱法命令流与手算验证过程。
以时程法结果进行比较。
分析不同阻尼值,,的影响。
标准谱 (1g=s2) (设计地震动值为频率谱值(g)339与标准谱对应的两条人工波见文件与2数值分析框图思路与理论简介理论简介该问题主要牵涉到结构动力分析当中的时程分析和谱分析。
时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。
谱分析是模态分析的扩展，是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。
分析框架：时程分析：在X和Z两个水平方向地震波作用下，提取底部跨中单宽上的剪力、弯矩值和顶部水平位移，并求出最大响应。
谱分析：先做模态分析，再求谱解，由于X和Z两个方向的单点谱激励，因此需进行两次谱分析，分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。
3有限元模型与荷载说明有限元模型考虑结构的几何特性建立有限元模型，首先建立平面几何模型，并将模型进行合理的切割，采用plane42单元，使用映射划分网格的方法生产平面单元（XOY平面）。
然后，采用solid45单元，设置拖拉方向的单元尺寸并清楚初始平面单元plane42，将平面单元进行拖拉，最后生成发电厂墙体的有限元立体几何模型。
单元总数为6060个，总节点数为8174个，有限元模型如图2所示：图2 发电厂墙体有限元模型荷载说明时程分析：首先计算结构的前两阶自振频率，分别为126.10008.2867f f ==，。
某连续梁桥地震反应谱分析
纵向剪力
０ ０ ０
横 向剪力 Ｆ
１０ １２ ３ ． ２４９ ６ １． ３２６ ８
纵 向剪 力
Ｏ ０ Ｏ
用条件下 ， 制动 墩主要 承担桥 梁上部结 构所 引起 的 向剪力 ； ３ 由计 算结 果可以看 到 ： 向剪力各墩 比较平 均 ， ） Ｆ 这是 由于 桥梁支座为单 向滑动 的 支座 ， 拟 时要约 束 住 Ｆ 模 向的 自由度 。 因此各墩共 同承担这个方 向的剪力 ；
本文采用 “ ｍ法 ” 模拟桩土相互作用 ， 此法是一 种常用 的模拟
ｆ 模型 『
结构桩一土相互之间作用力
输入地震力 模式
桩基静力 的方法 ， 在设计单位 经常采用 ， 在建模 时 ， 据本桥 结构 依
的特点 ， 采用三 维梁 单元 离散桥 梁 的主题 结构 ， 括桥 墩 、 梁、 包 主
向的组合 ， 不模 拟地震 的竖 向作 用力 ， 分别对 两种模 型地 震反 应
图 ２ 箱 梁横 断面 Ｉ 单位 ：ｉ ） ｔｌ １
进行模拟 ， 见表 １ 。
表 １ 反 应 谱 分 析 模 型 和 输 入地 震 力模 式
２ 桩一 土 结构 相 互作 用分 析方 法和 结构 模型 建 立
１． ８ ６ ２１． １
１ ． ４４
３２． １ ４９． ７
５８ ．
１ ２ ０． ｌ １ ６．
８６ ．
２ ９ １． ３ ６ ３．
模拟 ， 向活动支 座具有 竖 向转 动和单 一滑 移性 能 ； 单 固定 支座 仅 具有竖 向转 动性能 。本文采用两种模 型来进行 模拟 ， 不考 虑桩 土 相互作用 的模 型 Ａ， 墩第固结 ， 图 ３所示 ， 一种采 用需要 考 将 如 另
连续梁桥地震反应谱分析
平 安湟水河 铁路 桥为一联三跨预应力 混凝土连 续梁 ， 桥梁 的
韫 Ｒ
恒
０ ０ ０ ０ ０ ０ ｏ ｏ ０ ｏ （ ） ０ ０ ０ ０
结构形式 为（ ４ ８＋ ８ ０＋ ４ ８ ） ｍ， 截 面为单 箱单室 、 变高度 、 直腹板 、 箱
壶１ ０ ０ ０ ０ ０
壶８ ０ ０ ０ ０
｛ ３ Ｉ ｒ ６ ０ ０ ０ ０
４ ０ ０ ０ ０ ２ ０ ０ ０ ０
墩必须满足强度及延性 的要求 ， 因此 大跨 连续梁桥 固定墩 大多 采
用钢筋混凝 土重 力式 桥墩 。
ｏ ０ Ｏ
Ｚ
１ 工 程概 况及 动 力分析 模型
形结构 ， 主墩墩顶 ５ ． ０ ｍ范 围内梁 高均 为 ６ ． ６ ５ ｍ， 跨 中及边跨 现
跨度， ｍ
浇段梁高 ３ ． ８ ５ ｍ。梁底 曲线为 二次抛 物 线 ， 梁底 抛 物线 方程 为
图 １ 工 况 １下 主 梁 内 力 响应 曲线
Ｙ＝ ０ ． ０ ０ ２ ｌ ｘ ； 箱梁顶宽 １ ２ ． ２ ｍ、 底宽 ６ ． ７ ｍ， 单侧悬 臂长 ２ ． ７ ５ ｍ， １ ． ２ 有 限元 计 算模 型及动 力特 性 分析 悬 臂端部厚 ２ ４ ． ８ ｅ ｍ， 悬臂 根部厚 ６ ５ ｃ ｍ。箱梁 腹板为 直腹板 ， 腹 借助大型通用分析软件 Ａ Ｎ Ｓ Ｙ Ｓ建立 了大桥 三维有 限元计算 板厚度 由箱梁梁体 主墩 墩顶 根部 ９ ０ ｅ ｍ 变至跨 中及边墩 支 点 附 模型 。箱梁截 面采 用 单元 Ｐ ｌ ａ ｎ ｅ ８ ２建 立 ， 主梁 和 桥 墩 采 用 单 元 近梁段 ４ ８ ｃ ｍ。底板在箱梁梁体 主墩墩顶 根部 ９ ０ ｃ ｍ变 至跨 中及 Ｂ ｅ ａ ｍ１ ８ ８模拟 ， 主梁节 点与桥墩 用节点耦 合来模 拟 支座 。两 主墩 边跨直线段厚 ４ ０ ｃ ｍ。顶板厚 ４ ０ ｃ ｍ， 其 中箱梁梁体边墩顶根部加 底部固结 ， 主梁与边墩处约束竖 向位移及绕纵桥 向的转动 自由度 。
长周期结构地震反应的特点和反应谱
第十届中日建筑结构技术交流会南京长周期结构地震反应的特点和反应谱方小丹L2，魏琏3，周靖21．华南理工大学建筑设计研究院2．华南理工大学亚热带建筑科学国家重点实验室3．深圳市力鹏建筑结构设计事务所AbstractThe charaCte ri sti cs of eanhqmkc rcsponse and rcspo 璐e spec 咖f-or10n 争periods 虮lctI 鹏s a r ediscllssed ．A few shonages exist ing in the re$oIlse spectn 蚰of cllim code f-or seisIllic desi 驴of bllildin gsare 锄alyzcd ．11here a r eint 锄l relatio 雎be 抑een pseudo —accel 蹦ltion spec 仃l ：I 驰pseudo —Veloc 埘spectrI 珊and displace ment spec衄切珥th 盯ef ．0陀，a rt 诳ciaI modification to respo 嬲e spec 仃1蚰can re sll lt in the distonionof 争眦d m 嘶∞cha 髓c 白耐stics ．The 10ng ．p 嘲ods e gI]∞nt in rcspo璐espe 蛐ofC11im codc is revised ，infact ，蓼omld motion characte ri sti cs a r e c}姗ged ，wllich resul ts in an abn 咖l representati∞ofpowe rspcc 乜狮cofresp 伽成ng to acceleration spcctrIlm ，Milli 舢加storey seisIIlic she 甜coefj(icient described in thcspecificati 衄is oIlly relatcd to maximl earthqum(e innuence coef|ficient(％m)，but is not related to siteclassificatio 玑w 址ch is in connict 谢th the ge∞ral mles tllat the eanhqualke respo 璐e of as 仉l 咖re at thesoR·soil site is la 唱cr than tllat ofa s 甘uc 眦at tlle h 踟．d —soil site ．Accordingto the pseudo spectnlm rela ti on sbet 、)l ，e ％pseud0．accel 训on spectrIlIIl ，ps 即do-veloci 够spec 虮Imand dis placem ent spec 觚l 驰a responsespec 仃IlIIl pattcm 、Ⅳith lonj 雪er ．period segment(一10s)is proposed ，and whj!ch c a n pro 、，id c the refhence tospecificati 傩revision ．1(eywords lon 哥p 耐od ．s 仃Ilc 眦s ；response spec 胁；displacement specmml ；111iIlimum storey seisIllicshear coe伍cient ；seisIIlic desi 驴1引言有多种关于长周期结构的定义，如欧洲抗震设计规范认为基本振动周期大于3s 的结构为长周期结 构，我国抗震设计规范认为基本振动周期大于5s 的结构为长周期结构。
地震动位移反应谱分析及抗震设计谱
03
抗震设计谱基础
抗震设计谱定义
地震动位移反应谱是一种描述结构在地震动作用下，各个方 向和各个频率地震动加速度、速度和位移反应的曲线。
抗震设计谱是基于地震动位移反应谱，针对特定结构和场地 条件，进行结构抗震设计和分析的工具。
抗震设计谱特性
地震动峰值和频谱形状：这些特性可以根据场地条件和 地震危险性评估来确定。
非线性特性是由于地震动强度与结构 位移反应之间的非线性关系所导致的 。在地震动强度较小的情况下，结构 位移反应与地震动强度呈线性关系； 而在地震动强度较大的情况下，结构 位移反应的增长速度会逐渐放缓。
随机性特性是由于地震动的随机性所 导致的。地震动是一种复杂的自然现 象，其运动规律难以精确预测，因此 地震动位移反应谱也是随机的。
结构阻尼矩阵：结构阻尼矩阵可以包括质量阻尼矩阵和 刚度阻尼矩阵，用于描述结构在地震动作用下的振动特
性。
抗震设计谱通常具有以下特性
结构自振频率和阻尼比：这些特性可以根据结构类型和 尺寸来确定。
抗震设计谱编制方法
基于地震动位移反应谱的抗震设计谱编 制方法通常包括以下步骤
根据位移反应谱，进行结构抗震设计和 分析，得到结构的抗震性能指标和设计 参数。
地震动位移反应谱分析及抗 震设计谱
2023-11-08
目录
• 引言 • 地震动位移反应谱基础 • 抗震设计谱基础 • 地震动位移反应谱与抗震设计谱的关系 • 应用案例分析 • 结论与展望
01
引言
研究背景与意义
地震是一种严重的自然灾害，给人类社会带来了巨大的损失。因此，对地震动位 移反应谱进行分析，对抗震设计谱进行研究，对于减轻地震灾害具有重要意义。
抗震设计谱编制
根据地震动位移反应谱，结合建筑物的抗震 设防要求，编制出适用于该地区行分析，评估其合理性和有 效性，确保其能够充分考虑地震动的特性和 规律，为建筑物的抗震设计提供科学依据。
ABAQUS地震反应谱分析
ABAQUS反应谱法计算地震反应的简单实例Fan.hj2010年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（5）Create section，name：I，category：beam，type：beam, Continue, create profile, name:I, shape:I, 按图1尺寸输入界面尺寸，ok。
ANSYS地震反应谱SRSS分析
ANSYS地震反应谱SRＳＳ分析ﻫﻫ我在ＡNSYS中作地震分解反应谱分析,一次Ｘ方向，一次Y 方向,他们要求是独立互不干扰的，可是采用直进行一次模态分析的话,他生成的*．mcｏm文件好像是包含了前面的计算结果,命令流如下：!进入PREP7并建模ﻫ/PREＰ7ﻫB=15 !基本尺寸A1=1０00！第一个面积ﻫA2=１00０!第二个面积A3＝1000!第三个面积ﻫET,1,beam4 !二维杆单元ﻫR，1,0.2５,0.００5２,０.00５2,0.5,0.5!以参数形式的实参ﻫMP，EX,１,2.0E１１!杨氏模量ﻫmｐ,PＲXY,1,，0.3 ﻫmp,ｄenｓ,１,7.8e３ﻫN,1,-B,0,0 !定义结点ﻫＮ,２,0,０,0Ｎ,３，-B，0，bＮ,４,0,0,b ﻫN,5,-Ｂ,0，2*b ﻫＮ,6,0,0,2*b ﻫN,7，-B，0,３*bＮ,8,0，0,３＊b ﻫE,1,3 ！定义单元ﻫE,2,4ﻫE,3,5Ｅ,4,６ﻫE,3,４E，５,6e,５，7ｅ,6,８ﻫe,7,8ﻫD,1,ALＬ,０,,2ＦＩNISＨﻫ!!进入求解器，定义载荷和求解ﻫ/ＳOＬUＤ,１,ＡLL,0,,2 !结点ＵX=UＹ=0 ﻫsfbｅａm,1,1,PRES，１0000０,sfbｅａm,３，1，PREＳ,10００00,ｓｆbｅam，7,1,PRES，1000０0,SOＬVEＦINIＳHaｌlｓelNＭＯＤＥ=10 ﻫ/SOL!＊ﻫANTYPE,２ﻫ！*MＳＡＶE,0 ﻫ!*MODOPT,LANB,NMＯDE ﻫEQSLV,SＰAR ﻫMXPANＤ,NMODＥ,, ,1ＬUMPＭ,0PSTRES,0!* ﻫMODOPＴ,ＬANB,ＮMＯＤE ，０，0, ,OFFSOＬVE ﻫ*DIM，ＦRＥ,，NＭODEﻫ*DO,I,1,NMODE*ＧET,FRＥ（I),ＭＯＤE,Ｉ,FRＥQ ！OBTAＩN ＭODE FRＥQEＮCＹFOR MOＤＥI*ENＤDOFＩNISＨﻫ!地震影响系数ﻫgrav=9.81ﻫｔｇ=0.35amax=０.08ｃ=0.05! ﻫ*ｄim,a,,nｍodｅ*dim,ｔ,,nmoｄeﻫ*do,i,1,nmｏde ﻫｔ(ｉ）=1.０／fre(i）ﻫ*enddo ﻫｒ=0．9+(０．05－c)／(0．５+5．0*c）ﻫp１=0.０２+(0.05-ｃ）/8 ﻫp2＝1+(０.05-c）/(０.０6+1.７*c)ﻫ*do,i，１，ｎｍode*if，t(ｉ），ｇｅ,0.0,aｎd,t(i),lｔ,0.1,then ﻫa(i)=(0.4５+(10．0*ｐ2-4.5)*t(i)）*aｍax＊ｇｒaｖﻫ*ｅlseiｆ,ｔ(ｉ),ge,0.1,and，t(i),lｅ,tga(ｉ)=p２*amａｘ*gｒａｖ*elseiｆ,t(i),ｇｔ,tg,ａnd,t(i)，ｌe,5*tｇﻫa（ｉ)=(tg/t(ｉ))*＊ｒ*p２*aｍaｘ*gravﻫ*else a（ｉ)=（ｐ2＊0.2**r－ｐ1*（t(i）-5*tg))*ａmaｘ*graｖ*endif ﻫ*enddo!ﻫ! Ｘ-方向谱分析Specｔrum anａｌysiｓalｏnｇＧlobal X－axｉs diｒectｉon ﻫﻫ/SOLUAＮTＹPE，SPECＴR !Spectrum ａｎalyｓiｓﻫSPOPT,SＰRS!Sｉｎgle ｐoint sｐec ｔｒum ﻫＳED,1,，！Gloｂal Ｘ-axiｓaｓspectrｕｍdｉrｅctｉoｎﻫSVＴＹP，２! Ｓｅismic acceleｒaｔion ｒeｓｐonse sｐectrｕm ﻫﻫ！Fｒeqｕenｃｙpoints ａndＳpe ｃtrum ｖａlｕes ｆor SＶvs. freq．tableFＲEQ,fre(1),fre(２)，ｆrｅ(3)，ｆｒe(4）,fre(5),fre(6),ｆre（7),fre（8),ｆｒe(９）ﻫFREQ,fre(１0) ﻫSＶ,,a(１）,a（2)，ａ（3）,a(4),a（5)，ａ（6),a(7),ａ(8),a（9) ﻫSV,,ａ(10）FIＮISH!/SＯＬU ﻫ!ANTYPE,MODAL! Mode-fｒeｑuency aｎalysis ﻫ!EXＰAＳＳ,ON!MXPAND,nｍｏde,，,YEＳ，0.0 ! Expandｎmode sｈapeｓ, calculate eｌｅment stｒｅsｓｅs！SOLVE！FＩNIＳH ﻫ/SOLU ﻫANTYPE,SPＥCＴRSRSＳ,0.0,ＤIＳP！Squａre Rooｔof Suｍof Sｑuarｅs Modｅcｏmbｉnaｔioｎﻫ！witｈsiｇnif＝0.0 and dｉｓplａcemｅnt solution reｑueｓted ﻫSOLVE ﻫＦＩNＩSHﻫ/PＯＳT１SEＴ,LIST/ＩNPUT,，mcoｍﻫ!＊*************＊EＡRTHQＵAＫE Ｘ***************＊＊*ALLＳEL,ALLﻫFＩNISHﻫ! Y-方向谱分析Specｔrｕm ａnalysis ａlｏｎg Global Ｘ－axis direction!!****＊***＊＊****＊**********＊*****＊*＊＊***********!／SOL ﻫ!！＊！ＡNTＹPE,2!!* ﻫ!ＭSAＶE，0！!* ﻫ!ＭODOPT，LAＮＢ,NMODE ﻫ!EQＳLＶ,ＳＰAR ﻫ！MＸＰAND，NＭOＤE , ，,1 ﻫ！LUMＰM,0!ＰSＴRES,0!!＊!MＯＤＯＰT，LAＮB,NMODE ,0,0, ，OFF ﻫ!ＳOＬVE!FIＮISH!!*＊****＊******＊***＊＊＊***＊＊*****************＊*＊＊ﻫ／SOLU ﻫLSCLEＡR,LＳＯＰTﻫAＮTYPE，SPECTR !SｐｅctruｍａnａlyｓisﻫSPOPT,SPRS ! Single pｏｉnt ｓpecｔru ｍSED,，１, ！GlobaｌＹ-aｘｉs as sｐｅｃtｒｕｍdiｒectｉon ﻫSVTYP,２! Sｅismｉｃacceleｒａtion rｅsponsｅｓpectrｕmﻫFＲEQ!FrｅqｕencｙpoinｔsａnｄSpectrｕm ｖalｕes ｆoｒSＶvs.freq.tabｌeFRＥQ,fｒe(1),fｒe（２),fre（3),fｒｅ(4),fre(5),fre（6),fｒe(7）,fre(８)，ｆre（9)ﻫFREＱ,f ｒe(１0) ﻫSV,,a(1）,a(2），ａ（３),a(4),a（5),ａ(6）,a(7),a(8),a（9)ﻫSＶ,,a（10) SOLVE ﻫFIＮIＳＨﻫ!／ＳOLU！ＡNTＹＰＥ,MODAL ！Modｅ－ｆrｅquｅｎｃｙａｎalｙｓis!ＥＸPAＳＳ,OＮ!MXPＡND，nmoｄe，,,YＥS,0.０! Eｘpanｄnmode ｓｈａｐｅs, caｌｃulaｔe element stre ｓseｓ!ＳＯLVE ﻫ!FIＮISHﻫ/SＯLUANTYＰＥ,SPECTR ﻫSRSS，0.0，DISP! Square Ｒoｏt oｆSum of Ｓquａres Moｄe combinａtｉoｎﻫ !with siｇniｆ＝0.0 and disｐlaceｍenｔｓｏluｔiｏn ｒeｑuested ﻫSOLVEFＩNISH/POST１ＳET,LISＴﻫ/IＮＰUT,,mcｏm ﻫ！＊＊******＊*＊**＊*EARTＨQUAKE Y＊＊************＊***ＡLLSEＬ,ALL ﻫＦINＩSＨ这里在进行X方向的反应谱分析以后,进行Ｙ方向的分析，可是他生成的*.mcom文件如下:/COM,ANＳYS RELＥASＥ8．０UP2００3０9３0 09:28:42 0７/2３/２0０5／COM, truｓs．mcomLＣOPＥR,ZEROＬCDEFI，1, 1，1ﻫＬCFACT，1, 0.２638２５E-17ﻫＬCASE,1ＬCOPER，SQUAＲEﻫLCDＥFＩ,１，1, 2LCＦACT,1，8.55７78ＬＣＯPER，ADＤ,1,MＵＬＴ,1LCDEFI,1，1， 3 ﻫLＣFACＴ，1, -0.１886６９E-13LCＯPＥR，ADＤ,１，MUＬT，1LＣDEFI,1, 1, ４ﻫＬＣFＡCT,１，-0.871０９9E-15LCＯPEＲ，AＤD,1,MＵLT,1LＣＤＥFI,１, 1，5ﻫＬCFAＣT，1,-0.75７0１3LCOPＥR，ＡDD,1,MULＴ,1 ﻫLＣDＥＦＩ，１, １, 6 ﻫLCFAＣT,1, 0．9６7307E－１3 LCOPEＲ，ADＤ,1,MULT,１ﻫＬCDＥＦＩ,１，１, 7 ﻫLCFAＣT,１, ０.533１４1Ｅ-13 LCOPER,ＡＤD,1,MULT，1 ﻫLCDEFＩ,１, 1, 8LCFAＣT，1, -0.２03６９9LＣOPER，ADD,1,MULT,1ﻫＬCDＥFI，１，1，9LCFACT，1, 0.４45795E-13 ﻫＬCＯPER,ADD,1,MULＴ,1LCDEFI,1, １, 1０LCＦACＴ,1, -0.387８０8E-１3LＣOPER，ADD,1,MULT，1ＬCOPER,SＱRＴ/COＭ,AＮＳYS RELＥAＳE8．0 ＵＰ20030930 09:28:42 0７/２3／2005／CＯＭ, trｕｓs.mcomLＣOPＥR，SQUARE !注意这里没有清空数据库LCDEＦＩ,1，1, １ﻫＬＣFACＴ,1, 50.75２８LCOPER,ADD,1，MＵLT，1LCDEＦI，1,1, 2 ﻫＬCFAＣT,１, ０.88７０１７E-14ＬCOPER,ADD，1,MＵＬT，1 ﻫLCDEFI,1,1,3LCFＡCT,1, 0.６12８2４E-1３ﻫＬCＯＰＥR,ADD,１,MULT,1 ﻫLＣDEＦI，１, 1，４ＬＣFACT,1, -１．96484LCOPＥR,ＡＤD，１,MUＬT,1 ﻫＬＣＤＥFI,1, １, 5 ﻫLＣFACT,1, －0．331６13E-1３ﻫＬCＯPＥＲ,ＡDD,1，ＭULT，1LCDＥＦI,１，1, 6 ﻫLＣFAＣＴ,1，０．330４5９E-１３ﻫLCOPＥR,ＡDD,１,MＵＬＴ,1ﻫL ＣDEFI，1，１, 7LＣＦACT,1,0.3665６９ﻫLCOPER,ADD，1,MUＬT,1 ﻫLCDＥＦＩ，1，１, ８ＬCＦACT,1，-０．9７69９１E-13 ﻫLＣOPER,AＤD,1,ＭULT，1ＬCDEＦＩ，1, 1,9ﻫLCＦＡCT,1, 0.4173１３Ｅ-1３LCOＰER,AＤD,１,MULT,1LCDEＦI,１, 1, 10LCFACT,１, ０.４0104０E－13LCOＰER,ADＤ,1,MＵLT,1ＬCＯPER,SQRT我感觉这样好像是Ｘ和Y两个方向地震的叠加, ﻫ可是如果在座Y方向的地震以前把注释掉的模态分析在做一下这样的Ｙ方向的地震的*.mcｏm就是：/COM, truss.ｍcomLＣＯPER,ZＥRＯ!注意这里清空数据库ﻫＬＣDEＦＩ，1，1, 1LCFACT,1，50.7528ＬＣOPER,ＡDD,1,ＭULT,1LCDＥFI，1, １，2ﻫLCFＡCＴ，１, 0.８８7017E-14 ﻫLCＯPER,ADＤ,1,MUＬT,1 LCDEFＩ,１, 1，３ﻫLCFACT,１, 0.６１2824E-13 ﻫLCOPER,ADD,1，MＵＬＴ,１LCDEＦI,１, 1, 4 ﻫLCFACT,1, -1.96484 ﻫLCOPＥR,ADD,1,ＭUＬT，1 ﻫLCDＥFＩ，1, 1, 5ＬＣFACＴ,1，-0.３31613E－13ﻫLＣOPＥR,AＤＤ,1,MUＬＴ，1 ﻫＬCDEFＩ,１,１， 6 ﻫLCＦACT,1,0.3３0459E－13LCＯPEＲ,ADD,1,MUＬT,１LCDEFI,１, 1，7LCFACＴ,1，0.366569 ﻫＬCOＰER,AＤD,1,MUＬＴ，1ＬCDEFＩ,1, 1, 8LCFACＴ,1, -0.9７6９9１E-13LCＯPＥR，ADD,１,MULT,1LCDEFI,1,1, 9LＣFACT,１，０.4１７313E-1３ＬCOPＥR,ADD,1,ＭULＴ,1ＬＣＤＥFI，1,1, 1０LＣFＡCT,1, 0.401０4０Ｅ-13LCOPEＲ，AＤD,1,ＭULT，1LCOPＥR,ＳQRＴﻫ如果在X方向后不作Y方向的地震,他的*．mcom:/COM,ＡNＳYＳREＬEＡＳＥ８.0 UP200３093008：46:23 ０7/2３/20０５ﻫ／ＣOM, tr ｕsｓ．mcomＬCOＰER，ZＥRO ﻫＬCDEＦＩ,１,1,1ＬＣＦACT，１, 0.263８25E-1７LＣASE,１ﻫLCＯPEＲ,ＳQUＡRE ﻫLＣDEＦI，1, 1，２LCＦＡCT，1, 8.５5778LCOPER,ＡDD，1,MULT,１ﻫLCDＥFＩ,1, 1, ３ﻫLCFACＴ，1, -０.1８8６６9E-13LCOPＥR,ＡＤD,1,MUＬT,1 ﻫLＣDEＦI,1, 1, 4LCFACT,1, -0.871099E-１5ＬＣOＰER，ADD，１,ＭULT,1 ﻫLCDEFＩ,１，1, 5LCＦACT，1, －0．75７013 ﻫＬＣOＰER,ADＤ,1,MＵLＴ,1LCDEFI,1, 1,６ＬCFAＣＴ,1,０.967307E-13ＬCOＰEＲ,ＡDＤ，1,MULＴ，1LCDEＦI,1, 1，7ﻫＬCＦＡCT，１,０.533141E-13 ﻫＬCＯPER,AＤD,1,MULＴ，1ＬCDＥFI,1, １, 8ＬCＦACT，1,-０．２0３6９9 ﻫLCＯPEＲ，ADD，1,ＭULT,1LCDEFI,1, 1,9LCFＡCT,1, 0.44５795E-１3ＬCＯPER,ＡDD,1，MULT，1ＬCDEFＩ,1,1，1０ﻫLCFＡCT,1，-0.3878０8E-13ＬＣOPER，AＤＤ,1，ＭＵLT,1LCOＰER,ＳQRT ﻫﻫ可是在Ｘ后作Y他不清空数据库，需要进行两次模态分析，这很耗时间对于大型结构，请大家讨论讨论如何处理呢?Re:讨论：ANSYS地震反应谱ＳRＳS分析本人是学土木工程的，平时主要用Pａtran+Nastran对结构做线性分析,偶尔使用Ansys对结构做地震反应谱分析，但对Aｎｓyｓ的命令流不熟悉。
地震动反应谱特征周期计算地震荷载
选取同一类场地、震中距相近的20条地震动记录，地震动峰值均为0.7m/s2，单自由度结构的阻尼比为2%、5%、10%和15%，周期范围为0.1s~10s，计算位移反应谱、速度反应谱和伪速度反应谱、加速度反应谱和伪加速度反应谱，并分析比较速度反应谱和伪速度反应谱的区别，以及加速度反应谱和伪加速度反应谱的区别。
一.反应谱计算与绘图反应谱的计算采用Newmark-β法计算，对于单自由度体系使用杜哈美积分来求解实际更为方便。
MATLAB的计算程序如下所示：clcclearkesai=0.15; %阻尼比m=1;[acc,dt,N]=peer2acc('F:matlab-learn','RSN3753_LANDERS_FVR135.AT2')%peer2acc为处理原始地震动数据的程序save('acc2','acc')load('acc2.mat');gama = 0.5;beta = 0.25;alpha0 = 1/beta/dt^2;alpha1 = gama/beta/dt;alpha2 = 1/beta/dt;alpha3 = 1/2/beta - 1;alpha4 = gama/beta - 1;alpha5 = dt/2*(gama/beta-2);alpha6 = dt*(1-gama);alpha7 = gama*dt;peak=9.8*max(abs(acc));acc=acc*0.7/peak;n=length(acc);p=-m*9.8*acc;j=0;for T=0.1:0.01:10j=j+1;wn=2*pi/T;k=m*wn^2;c=kesai*2*m*wn;Keq=k+ alpha0*m + alpha1*c;wD=wn*(1-kesai^2)^0.5;d=zeros(n,1);v=zeros(n,1);a=zeros(n,1);for i=2:nt=0.002*(i-1);f=p(i) + m*(alpha0*d(i-1)+alpha2*v(i-1)+alpha3*a(i-1))+c*(alpha1*d(i-1)+alpha4*v(i-1)+alpha5*a(i-1)); d(i) =f/Keq; %Newmark-β的计算程序a(i) = alpha0*(d(i)-d(i-1))-alpha2*v(i-1)-alpha3*a(i-1);v(i) = v(i-1) + alpha6*a(i-1) + alpha7*a(i);endsd(j)=max(abs(d)); %位移反应谱sv(j)=max(abs(v)); %速度反应谱sa(j)=max(abs(a)); %加速度反应谱SA(j)=wn^2*sd(j); %伪加速度反应谱SV(j)=wn*sd(j); %伪速度反应谱end选取的地震动记录如图地震动记录一般在PEER网站下载。
典型地震反应谱参数分析
典型地震反应谱参数分析地震反应谱是一种用于描述地震动力学特性的图像或函数，它反映了地震对结构物产生的力或位移随时间的变化规律。
地震反应谱参数分析是对地震反应谱进行统计和分析，以评估地震对结构物的可能影响，并为工程设计和地震工程防护提供依据。
在进行典型地震反应谱参数分析时，常见的参数包括峰值加速度、峰值速度、峰值位移、特征周期等，这些参数可以通过对地震反应谱曲线进行解析和计算得到。
首先，峰值加速度是反应谱曲线中离地面最大加速度的数值。
它是衡量地震对结构物产生的震动强度的重要指标。
在地震工程设计中，通常通过地震加速度响应谱曲线的峰值来判断结构物的耐震性能，并选择合适的设计加速度。
峰值加速度的值越大，表示地震对结构物的影响越强烈。
其次，峰值速度是地震加速度响应谱曲线中离地面最大速度的数值。
它是描述地震动力学效应的另一个重要参数。
峰值速度的值可以通过将加速度响应谱曲线进行一次积分得到。
在地震工程中，峰值速度的大小可以用来评估结构物的损伤程度和破坏概率。
峰值位移是地震加速度响应谱曲线中离地面最大位移的数值。
它是描述结构物在地震作用下产生位移变化的指标。
峰值位移可以通过对加速度响应谱曲线进行二次积分得到。
在地震工程中，峰值位移的大小通常用来判断结构物的破坏程度和变形情况。
特征周期是地震反应谱曲线中的一个重要参数，它是指加速度响应谱曲线中对应峰值加速度的周期。
特征周期是用来描述结构物振动特性的指标，可以通过对地震反应谱曲线进行周期化分析得到。
特征周期的选择对于结构物的抗震设计和地震防护具有重要意义，不同结构物对地震的响应特征周期有不同的要求。
除了上述参数，地震反应谱参数分析还可以包括剪切强度、硬度指标、阻尼比等其他参数。
这些参数的分析可以提供更加全面和详细的地震动力学特性信息，对于结构物的抗震设计和地震工程防护具有重要的参考价值。
总结起来，典型地震反应谱参数分析是对地震反应谱进行统计和分析，通过计算和解读峰值加速度、峰值速度、峰值位移、特征周期等参数，评估地震对结构物的可能影响，并为工程设计和地震工程防护提供依据。
el centro地震反应谱
el centro地震反应谱
El Centro地震反应谱是指根据El Centro地震记录所绘制的一
种地震反应谱。
El Centro是位于美国加利福尼亚州的一个城市，在1940年5月18日发生了一次里氏7.1级的地震，被称
为El Centro地震。
地震反应谱是用来描述某个地点在地震中的震动强度和频率分布的图表。
El Centro地震反应谱主要包括加速度、速度或位
移随时间的变化。
这些反应谱曲线可以通过地震记录仪所记录的地震波来绘制。
通过分析地震反应谱，可以确定结构物在地震中的响应特性，如最大加速度、周期和频率内容等。
El Centro地震反应谱在结构工程中被广泛使用，用于评估建
筑物和其他工程结构在地震中的承载能力和耐震性。
它可以用于设计地震工程的基础参数，如地震设计负载以及结构物的设计和改进。
需要注意的是，El Centro地震反应谱是根据1940年的一次地
震记录绘制的，并不适用于所有地震情况。
不同地震的地震反应谱可能因地震规模、震源距离、地质条件等因素而有所不同，因此在实际工程中应根据具体情况选择适当的地震反应谱。
地震动反应谱方法
70年代
国内
我国对反应谱研究始于60年代，王前信等人做出了大量工作， 并制定了适合于我国地质条件的设计反应谱
华南理工大学
2.2.1 反应谱的定义
有阻尼单自由度弹性体系在地震作用下的运动方程：
t 2 x t 2 x t g t x x
（1）
（5）
t , x t , x t 分别为结构相对于基底的加速度、速度和位移列向量 x
I
为单位位移时候结构各节点产生的位移
用振型分解法对（5）式解耦，将 x t 表示为阵型叠加的形式：
x t j Y i t
j 1 n
（5）
由（3）、（4）和（5）式：
Sd x t Sv x t
相对位移反应谱
m ax
相对速度反应谱
m ax
S a t g t x x
相对加速度反应谱
m ax
简单来说，反应谱是指单质点体系地震最大反应与结构自振周期 之间的关系。它是跟阻尼比和周期有关的函数。
Tg为场地特征周期
0 . 1 max0
地震影响系数曲线
水平地震影响系数最大值
烈度
max
6 0.12
7 0.23
8 0.45
9 0.90
CHONGQING UNIVERSITY
华南理工大学
2.3 振型分解反应谱法
自由度为N的弹性体系在地面运动加速度为
ut 其运动方程为
M t C x t K x t M I ut x
j max
E max
j
t
D j ,
midas抗震设计-反应谱分析
北京迈达斯技术有限公司目录简要 (1)设定操作环境及定义材料和截面 (2)定义材料 (2)定义截面 (3)建立结构模型 (4)主梁及横向联系梁模型 (4)输入横向联系梁 (5)输入桥墩 (5)刚性连接 (7)建立桥墩和系梁 (9)输入边界条件 (10)输入支座的边界条件 (10)刚性连接 (11)输入横向联系梁的梁端刚域 (12)输入桥台的边界条件 (13)输入二期恒载 (14)输入质量 (15)输入反应谱数据 (17)输入反应谱函数 (17)输入反应谱荷载工况 (18)运行结构分析 (19)查看结果 (20)荷载组合 (20)查看振型形状和频率 (21)查看桥墩的支座反力 (24)简要本例题介绍使用MIDAS/CIVIL的反应谱分析功能来进行抗震设计的方法。
例题模型使用的是简化了的钢箱型桥梁模型，由主梁、横向联系梁和桥墩构成。
桥台部分由于刚度很大，不另外建立模型只输入边界条件；基础部分假设完全固定，也只按边界条件来定义。
下面是桥梁的一些基本数据。
跨径：45 m + 50 m + 45 m = 140 m桥宽：11.4 m主梁形式：钢箱梁钢材：GB(S) Grade3（主梁）混凝土：GB_Civil(RC) 30（桥墩）[单位：mm]图1. 桥梁剖面图设定操作环境及定义材料和截面开新文件（新项目），以‘Response.mcb’为名保存(保存)。
文件/ 新项目t文件/ 保存( Response )将单位体系设定为kN(力), m(长度)。
工具/ 单位体系长度>m; 力>kN ↵定义材料分别输入主梁和桥墩的材料数据。
模型/ 材料和截面特性/ 材料材料号(1); 类型>S钢材规范>GB(S); 数据库>Grade3 ↵材料号(2); 类型>混凝土规范>GB-Civil(RC); 数据库>30 ↵图2. 定义材料定义截面使用用户定义来输入主梁、横向联系梁以及桥墩的截面数据。
5第七讲 地震响应与谱分析
典型命令： 典型命令： FREQ，… ， SV，... ，
• B.多点响应谱分析（MPRS）:在模型不同的 B.多点响应谱分析（MPRS） 点集上定义不同响应谱曲线。不同的多个响 应谱分别激励模型中不同的点。 • C.动力设计分析（DDAM）:应用一系列经验 C.动力设计分析（DDAM） 公式和振动设计表得到的谱来分析系统。 • D.功率谱密度分析（PSD）:功率谱密度是 D.功率谱密度分析（PSD） 结构对随机动力载荷响应的概率统计，用于 随机振动分析。 • 注意：谱分析必须要已知结构的振型和固有 频率，因此需先进行模态分析。在扩展模态 时，只需扩展到对最后进行谱分析有影响的 模态即可。
1 2 3 4
谱分析
• 如果振动台以频率f1激振并且四个系 统的位移响应都被记录下来，结果将 如右图所示 • 现在再增加频率为f3的第二种激振并 记录下位移响应，系统1及3将达到峰 值响应 • 如果施加包括几种频率的一种综合激 振并且仅记录下峰值响应，就将得到 右图所示的曲线，这种曲线称为频谱， 并特称为响应谱
C、直接动力分析理论
谱分析
什么是谱分析？ • 它是模态分析的扩展，用于计算结 构对地震及其它随机激励的响应 • 在进行下述设计时要用到谱分析： − 建筑物框架及桥梁 − 太空船部件 − 飞机部件 − 承受地震或其它不稳定载荷的结构 或部件
谱分析
• 谱分析的一种代替方法是进行瞬态分析，但是： − 瞬态分析很难应用于例如地震等随时间无规律变 化载荷的分析; 化载荷的分析; − 在瞬态分析中，为了捕捉载荷，时间步长必须取 得很小，因而费时且昂贵. 得很小，因而费时且昂贵.
αmax
0.45αmax (Tg/T)0.9αmax
0.1
Tg
3
ABAQUS地震时程分析小算例
ABAQUS时程分析法计算地震反应的简单实例（在原反应谱模型上修改）Fan.hj2010年4月26日问题描述：（本例引用《有限元法及其应用》一书中陆新征博士使用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。
（4）进入property模块，create material，name：steel，general-->>density，mass density：7800，mechanical-->>elasticity-->>elastic，young‘s modulus：2.1e11，poisson’s ratio：0.3.（5）Create section，name：I，category：beam，type：beam, Continue, create profile, name:I, shape:I, 按图1尺寸输入界面尺寸，ok。
地震响应的反应谱法与时程分析比较
发电厂房墙体地震响应的反应谱法与时程分析比较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.19 0.2612.5 0.3130.25 0.047与RG1.60标准谱对应的两条人工波见文件rg160x.txt与rg160y.txt2数值分析框图思路与理论简介2.1理论简介该问题主要牵涉到结构动力分析当中的时程分析和谱分析。
时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。
谱分析是模态分析的扩展，是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。
2.2 分析框架：时程分析：在X和Z两个水平方向地震波作用下，提取底部跨中单宽上的剪力、弯矩值和顶部水平位移，并求出最大响应。
谱分析：先做模态分析，再求谱解，由于X和Z两个方向的单点谱激励，因此需进行两次谱分析，分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。
3有限元模型与荷载说明3.1 有限元模型考虑结构的几何特性建立有限元模型，首先建立平面几何模型，并将模型进行合理的切割，采用plane42单元，使用映射划分网格的方法生产平面单元（XOY平面）。
然后，采用solid45单元，设置拖拉方向的单元尺寸并清楚初始平面单元plane42，将平面单元进行拖拉，最后生成发电厂墙体的有限元立体几何模型。
单元总数为6060个，总节点数为8174个，有限元模型如图2所示：图2 发电厂墙体有限元模型3.2 荷载说明时程分析：首先计算结构的前两阶自振频率，分别为126.10008.2867f f ==，。
振型反应谱分析法和底部剪力法例题
334.2kN 334.4kN 167.4kN
− 120.8kN 120.7 kN 120.9kN
（4）计算各振型各楼层的水平 第一振型 地震作用 F j i = α j x j iγ j G i 第二振型 F 2 1 = 0 . 16 × ( − 0 . 428 ) × ( − 0 . 667 ) × 270 × 9 . 8 = 120 . 9 kN F 2 2 = 0 . 16 × ( − 0 . 428 ) × ( − 0 . 666 ) × 270 × 9 . 8 = 120 . 7 kN
例：试用振型分解反应谱法计算图示框架多遇地震时的层间剪力。 试用振型分解反应谱法计算图示框架多遇地震时的层间剪力。 抗震设防烈度为8 类场地，设计地震分组为第二组。 抗震设防烈度为8度，Ⅱ类场地，设计地震分组为第二组。 解： 1）求体系的自振周期和振型 （ m = 180t
3
{X }1
0 .3 3 4 = 0 .6 6 7 1 .0 0 0
m2 = 270t
K 2 = 195MN/m
m1 = 270t K1 = 245MN/m
T1 = 0 . 46 7 s T 2 = 0 . 2 0 8 s T 3 = 0 . 1 3 4 s
（2）计算各振型的地震影响系数 α 1 = 0 . 1 3 9 α 2 = 0 . 16 α 3 = 0 . 1 6 （3）计算各振型的振型参与系数 γ 1 = 1 . 363 γ 2 = − 0 . 4 28 γ 3 = 0 . 06 3
3
{X }1
0 .3 3 4 = 0 .6 6 7 1 .0 0 0
{X }2
− 0 . 66 7 = − 0 . 666 1 . 00 0
ABAQUS反应谱法计算地震反应实例
ABAQUS反应谱法计算地震反应实例问题描述：悬臂柱高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。
continue（3） Create lines，在分别输入0，0回车；0，3回车；0，6回车；0，9回车；0，12回车。
（4）进入property模块，create material，name：steel，general-->>density，mass density：7800mechanical-->>elasticity-->>elastic，young‘s modulus：2.1e11，poisson’s ratio：0.3.（5） Create section，name：Section-1，category：beam，type：beam,Continuecreate profile, name: Profile-1, shape:I,按图1尺寸输入界面尺寸，ok。
在profile name选择I，material name 选择steel。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
结构地震反应谱分析实例
在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！
为叙述方便，举一简单实例：
在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑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
LCASE,11
PRRSOL, !结果2
SET,FIRST
PRRSOL, !结果3
SET,NEXT
PRRSOL, !结果4
SET,NEXT
PRRSOL, !结果5
SET,NEXT
PRRSOL, !结果6 FINISH
!静力分析
/SOLU
ANTYPE,STATIC
!施加水压荷载
NSEL,,LOC,Y,0
NSEL,R,LOC,Z,0,5 SFGRAD,PRES,0,Z,0,-10000 SF,ALL,PRES,50000
!施加集中荷载
NSEL,,LOC,Y,0
NSEL,R,LOC,z,5
F,ALL,FY,10000
ALLSEL
EPLOT
SOLVE
FINISH
/POST1
set,last
lcwrite,12
Lcase,11
Lcoper,add,12
Lcwrite,13
LCASE,12
PRRSOL, !结果7
PRRSOL, !结果8
FINISH
---------------------------------------------------------------------------------------------------
以下是计算的结果
---------------------------------------------------------------------------------------------------
结果1：（Results Summary）
1 21.647
2 21.647
3 121.51
4 121.51
结果2：（单独谱分析反力LCASE,11）
VALUE 2467.9 2290.1 18384.
结果3：（单独谱分析反力SET,FIRST）
VALUE 0.13334E+06-0.15785E+07-0.18819E-06
结果4：（单独谱分析反力SET,NEXT）
VALUE -0.15785E+07-0.13334E+06-0.48918E-06
结果5：（单独谱分析反力SET,NEXT）
VALUE -0.87805E+07 0.27008E+08 0.86846E-07
结果6：（单独谱分析反力SET,NEXT）
VALUE 0.27008E+08 0.87805E+07 0.79325E-06
结果7：（单独静力分析反力LCASE,12）
VALUE 0.22901E-08-0.15500E+06 0.12000E+06
结果8：（谱分析与静力分析叠加反力LCASE,13）
VALUE 2467.9 -0.15271E+06 0.13838E+06
---------------------------------------------------------------------------------------------------
以下是问题的讨论
---------------------------------------------------------------------------------------------------
1、模态提取数为30，即取前30阶振型数，但在谱分析时得到的是4阶，这4阶是什么意思？
2、在单独谱分析时，为何结果2、
3、
4、
5、6会相差如此之大？（其应力和位移也是如此。
）
3、在进行谱分析的合并模态步骤中，模态合并方法应该选用SRSS,0.15,DISP还是SRSS,0.15,ACEL（其对结果影响很大）？
4、在谱分析中，似乎应该将水压力折算成各节点处的附加质量MASS21，但对于非竖值方向（如本例为Y 方向），MASS21单元如何加上去？
5、在选取（离散）频率和谱值时，似乎没有一定的标准，而其基本原则是须选取具有代表性的点。
6、模态分析中采用Subspace法与Block Lanczos法所得结果相差2%~4%左右，用后一种方法计算更快，但占用内存大50%左右。
7、在网格剖分时采用VMESH和VSWEEP两种方法，对静力分析结果没有影响，但对谱分析结果却有一定影响，何也？
8、有朋友说可以“做模态分析之前先做打开预应力的静力分析，然后再做打开预应力的模态分析，接着和以前一样”，但我试了试，得到的结果相差很大，不知何故，如果哪位朋友曾经作过这方面的例子，能否再试一下本例，以便进行比较！
SED,1,1
是X与Y的45度夹角方向，
要考虑X与Y向地震作用，
似乎应该为：
SED,1
SED,,1
（分别考虑）。