在 ansys 中如何 施加 地震波
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。
ANSYS地震时程分析
在ANSYS里做地震分析时,需要对结构施加地震惯性荷载,地震惯性力是通过加速度的方式输入进结构的,然后与结构的质量一起形成动力计算时的惯性荷载,下面说一下在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-010.200000E-010.300000E-010.400000E-010.500000E-010.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/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
关于ANSYS地震响应分析的一些讨论剖析
地震响应分析1模态组合就是根据模态分析中的几阶振型(也可以少于这几阶,看你要求的精度)进行组合(类似于结构最不利组合),从而求出地震响应的最大值。
2组合各振型反应的最大值,求得结构地震响应的最大值。
这个问题在论坛上已经有很多人问过,也有各种各样的回答,但是至今没有令人满意的解答。
我自己试过很多种方法,加上论坛上其他人提到的方法,大致归类如下:1.先做静力恒载工况分析,打开预应力pstres开关;然后转到时程分析。
结果:恒载对后面的时程计算不起作用,时程计算依然从0开始。
2.直接在antype,trans中考虑恒载:先把timint,off加acel,,9.81,打开应力刚化,sstif,on,lswrite,1,然后timint,on开始时程计算。
结果:恒载9.81起作用了,但结果是错的,它被积分了。
3.不用什么前处理,直接把9.81加在地震波上acel,9.81+ac(i)。
结果,同2,9.81带入了积分,这个9.81相当于阶跃荷载,而不是产生恒载。
4.ansys帮助中施加初始加速度的方法(篇幅限制请自己看帮助)。
结果,同2、3,9.81还是带进时间积分。
5.这种是我受到别人的启发,通过结构受ramp荷载的特点施加的,可以近似的解决问题。
即1)求出结构的自振一阶频率w2)令tr=1/w3 定义ramp荷载为从0到tr加到9.81,然后在整个时间积分中保持不变4)antype,trans中分几个荷载步将荷载从0加到9.815 在随后的荷载步中acel,,9.81+ac(i)这种做法虽然也是将9.81++加到地震波中,但是因为满足TR的要求,所以这个动力效应被削弱到了静力效应,它作用在结构上就像静载一样。
对于单自由度结构理论上跟静载是完全一样的,但是多自由度会子静力效应上下很小的范围内波动,所以可以认为相当于静载的作用,这样我们就可以达到考虑恒载的目的了。
第5种是我至今为止考虑恒载的做法,我也很想知道还有没有更简单精确的方法,或者在前4种方法中就有只是我使用不正确,希望大家能一起来讨论,彻底解决这个问题。
ANSYS地震分析算例
02 地震分析算例 (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。
ANSYS Example02地震分析算例 (ANSYS)
02地震分析算例(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-> MaterialModels菜单,在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic属性,输入材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。
(6)继续给Material Model Number 1添加Density属性,输入密度为7800。
基于ANSYS的框架结构地震分析教程(静力分析+模态分析+反应谱分析+LS-DYNA时程分析)
!---!---! !
QQ: 1549221758 Website: phipsi.top
-------
!****************************************************************
!**************************************************************** !----------------------------初始化-----------------------------!**************************************************************** finish /clear /FILNAME,ANSYS_Seismic_analysis_Static /Titie,Seismic analysis of the frame structure /PREP7 /DSCALE,ALL,10.0 /eshape,1 /UIS,MSGPOP,3 !进入前处理器 !后处理变形缩放系数 10 倍 !显示单元真实形状 !忽略警告信息
ANSYS软件在结构地震反应分析中所用的方法
ANSYS软件在结构地震反应分析中所用的方法徐旻洋 1110109132工程体系多自由度运动方程为:[M]{x’’}+[C]{x’}+[K]{x}={F(t)}(1)式中,[M]表示质量矩阵;[C]为阻尼矩阵;[K]为刚度矩阵;{x}为结构体系的位移向量;{F(t)}表示t时刻的载荷向量。
典型的无阻尼模态分析求解的基本方程就是上式(1)的特征值问题:[K]{Φi}=ωi2[M]{Φi} (2)式中,{Φi}是第i阶模态的振型向量(特征向量),ωi是第i阶模态的固有频率。
ANSYS软件可以求解式(2),计算结构的固有频率ωi,然后计算相应的振型向量{Φi},即模态分析。
当式(1)右边{F(t)}是一个已知的谱(如地震反应谱)时,可以用ANSYS软件进行谱分析。
当{F(t)}是任意的随时间变化载荷时,ANSYS软件可进行瞬态动力分析。
ANSYS结构振型分解反应谱分析有如下内容:1)首先要定义好加速度反应谱。
这里需要注意的是,规范上给的是地震影响系数谱曲线,这个曲线的函数值是以地面加速度为单位的。
而我们在用这个软件算的时候就需要给出绝对的加速度值,这个绝对加速度值当然就是要在地震影响系数的基础上再乘上一个地面加速度。
而地面加速度也并不一定是9.8,这与我们使用的单位制有关,如果是N/M/S,就应该是9.8,如果是N/MM/S就应该是9800。
2)求振型。
一定要是相对质量矩阵进行归一化,使用modopt命令默认的方法就可以了。
这个式子是求振型参与系数的,显然这个式子里面不是完整的求振型参与系数的式子,它少了分母,但是,由于对振型相对质量矩阵进行了归一化,这个分母就等于1了,这就是为什么必须要对振型相对质量进行归一化的原因了。
在这一步中,可以这样理解,程序只进行了一次特征值求解,即只求出了周期和振型。
如果需要看某个振型的“内力/应力/反力”,就需要对其进行模态扩展。
3)求谱解。
其实在这一步中,程序只做了一件事,那就是求模态系数。
利用ANSYS程序的二次开发功能进行面板堆石坝地震反应分析
1 工 程 背 景
九甸 峡水 利枢纽 工程 混凝 土 面板堆 石坝最 大 高 度 1 65 3 . m,电站 总装 机容 量为 3 0 0 MW 。
该工 程坝 址 区地震 基本烈 度 经鉴定 ,并 由国家
图 I S I OL I M5单 元 图
压 的应 力 应 变 关 系 ,直 接 用 于 堆 石 坝 地 震 反 应 的 非 线 性 瞬 态 分 析 。 【 关键 词】 面 板堆 石 坝 地 震 波 A Y NS S程 序 AP DL语 言 非 线 性 应 力 应 变 曲 线 弹 塑 性 模 型 围压 【 图分类号】 中 T 2 【 V3 3 文献 标 识 码 】 A 【 章 编 号 】 1 7 — 2 6 (0 6 4 02 一O 文 6 2 4 9 2 0 )0 — 0 8 3
・ 2 ・ 8
维普资讯
设计与施工
表 1 材 料 密 度 表
材 料 位置 垫层 料 过 渡料 主 堆 石料
水 利 规 划 与设 计
型中 ,用虚 拟温度 值作 为 围压值 。
20 年第 4 06 期
干 密度 ( / m。 gc ) 22 .8 2 2 .5 2 2 .O
・ — ’
P , 7 ’ f 、
jt ‘ 。
T, 。 。 0O器● , { l O
霎
:
二
— -o 酬 .-
—
’
/ ● .. _ .
×
r.
、
| /
SI £
I
『 ,
『 /
p- ’ ●
,
O 轴 向成 蹙 £ ( 1 )
ANSYS地震分析算例
ANSYS地震分析算例算例描述:在一个城市里,设计师计划建造一座高楼大厦。
然而,由于该地区位于地震活动带,设计师决定对该高楼进行地震分析,以确保其在地震时的稳定性。
假设该高楼大厦的结构为钢筋混凝土框架结构。
设计师希望通过地震分析来确定该结构的最大位移、最大应力和最大应变等参数,并评估结构的稳定性。
地震分析步骤:1.确定模型和边界条件:-在ANSYS中创建高楼大厦的模型。
可以使用ANSYS提供的建模工具,通过绘制位于平面上的轮廓线进行模型创建。
-定义地震加载条件:确定地震波的性质和加载方向,并将其应用于模型的适当位置。
-设置边界条件:确定结构的支撑方式,如固定支座或固定边界条件。
2.材料和单元属性设置:-模型中的材料属性:定义使用的材料的弹性模量、泊松比和密度等特性。
-定义单元属性:选择适当的单元类型,并设置单元的尺寸和属性值。
3.进行静态分析:-应用静态载荷:为了模拟自重和其他永久载荷,将这些载荷应用到模型中。
-执行静态分析:运行ANSYS分析以计算静态应力和变形。
4.开展地震分析:-应用地震波加载:将地震波加载应用于模型的适当位置。
-执行动态分析:运行ANSYS分析以计算结构在地震荷载下的动态响应。
5.结果分析:-输出结果:分析完毕后,ANSYS将提供关于位移、应力、应变和变形等参数的结果。
-结果评估:根据结果评估结构的稳定性和安全性。
可以根据设计准则或标准来判断结构是否合格。
总结:地震分析是建筑和结构设计中至关重要的一步。
ANSYS提供了强大的分析工具,可以帮助工程师和建筑师评估建筑结构在地震荷载下的响应。
通过ANSYS地震分析,设计者可以确定结构的稳定性和安全性,并采取必要的措施来增加其抗震能力。
二维黏弹性人工边界单元及地震波输入在ANSYS中的实现
水利与建筑工程学报
JournalofWaterResourcesandArchitecturalEngineering
DOI:10.3969/j.issn.1672-1144.2019.05.010
Vol.17No.5 Oct.,2019
二维黏弹性人工边界单元及地震波输入 在 ANSYS中的实现
对无限域或半无限域进行地震响应分析时,往 往采用有限元法。该方法通常采取的处理措施是在 无限域或半无限域空间中人为地截取一定范围的有 限域,进而可以将无限区域转化为有限区域,然后在 有限域的基础上,对其进行离散法处理,进而无限域
或半无限域最终被转化为有限个单元。同时为了达 到近似处理的效果,通常在截取的有限域的边界上 人为的施加约束,但是这一处理的缺点是在处理过 程中,不仅没有考虑远场无限地基的辐射阻尼效应, 同时半空间无限地基弹卷
截取的有限域的范围存在不确定性,随着有限域截 取范围的增加,自由度也随之增加。基于此离散范 围也更加宽广,进而网格划分的单元数及结点数成 倍增加,因此,在地震动的动力计算中工作量也进一 步被增加,对计算机的要求也更高,进而,在动力学 中非线性问题的求解难度进一步增大。目前,针对 该关键性问题处理的诸多措施中,较好的处理方法 是在截取的有限域边界处设置局部人工边界。目前 在工程中较常用的局部边界主要有黏性边界[1]、透 射边界[2]、黏弹性边界[3]等。黏性边界的优点是物 理概念清晰,在实际动力计算程序中实现起来相对 比较容易,对于 1阶动力响应精确度较高。基于此, 黏性边界单元被嵌入到了很多大型商用软件中以便 于工程的实际应用,如 LSDYNA、FIAC、ABAQUS[4] 等。它的缺点是在施加人工边界过程中只着重考虑 了人工边界辐射阻尼效应,忽略了人工边界处半无 限空间介质的弹性恢复性能。因此在实际工程应用 中容易发生整体漂移的现象,进而使计算结果精度 不高。透射边界的优点是对于 2阶动力响应计算结 果存在较高的精度,其缺点是动力计算分析程序的 代码编写较复杂、在实际工程应用中可能存在高频 失稳的现象。黏弹性边界不仅能有效的模拟半无限 空间介质的辐射阻尼效应,同时还能较好的模拟半 无限地基的弹性恢复性能。此外还具有良好的频率 稳定性和较高的精度。因此在实际工程中得到了广 泛的应用[5-10]。
ANSYS软件在高层建筑地震反应分析中的应用
ANSYS软件在高层建筑地震反应分析中的应用一、概要随着全球地震灾害的频发,高层建筑在地震作用下的安全性受到了广泛关注。
地震反应分析是评估高层建筑在地震作用下结构响应及损坏情况的关键技术手段。
ANSYS软件作为一款广泛应用于工程领域的有限元分析软件,在高层建筑地震反应分析中发挥着重要作用。
1. 抗震设计的重要性随着全球地震活动的增加,高层建筑在其生命周期中的抗震性能显得尤为重要。
高层建筑由于其特殊的结构形式和高度,往往成为地震作用下的关键受力构件。
如果在地震作用下发生倒塌,将造成巨大的人员伤亡和财产损失。
抗震设计的核心目标是确保建筑物在可能发生的地震中能够保持足够的稳定性和完整性,从而保护人员安全并减少财产损失。
ANSYS 软件作为一种强大的有限元分析工具,在高层建筑地震反应分析中发挥着至关重要的作用。
通过ANSYS,工程师们可以模拟和分析建筑物在地震作用下的动态行为,包括应力和变形分布、结构的失效模式以及能量耗散等。
这有助于设计师在建筑设计阶段就识别出潜在的薄弱环节,并采取相应的加固措施来提高建筑的抗震性能。
ANSYS还可以用于验证设计的合理性,通过与其他软件或实验结果的对比,确保建筑物在实际地震中的表现符合预期。
这对于保证高层建筑在地震中的安全性至关重要。
抗震设计是高层建筑安全性的重要保障。
ANSYS软件的应用使得这一过程更加高效、准确,为设计师提供了强有力的工具来应对地震带来的挑战。
2. 高层建筑地震反应分析的挑战随着城市化的加速和土地资源的紧张,高层建筑越来越多。
高层建筑在地震作用下的地震反应分析是一个复杂且具有挑战性的问题。
在地震作用下,高层建筑会受到水平、竖向和扭转等多种振动模态的影响,使得地震反应分析变得非常复杂。
高层建筑结构的耦合作用使得地震反应分析更加困难。
高层建筑结构中,各构件之间存在复杂的相互作用,如梁柱、梁墙、墙柱等。
这些相互作用使得地震力在建筑物内的传递变得复杂,难以准确模拟实际的地震反应。
ANSYS 地震分析算例
02 地震分析算例 (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。
在ANSYS里施加地震惯性力的方法
[分享]在ANSYS里施加地震惯性力的方法[b][size=16px]在ANSYS里施加地震惯性力的方法[/size][/b]在ANSYS里做地震分析时,需要对结构施加地震惯性荷载,地震惯性力是通过加速度的方式输入进结构的,然后与结构的质量一起形成动力计算时的惯性荷载,下面说一下在ANSYS 里施加地震惯性力的方法。
首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据:[font=宋体] -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[/font] .......................然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据:[font=宋体] 0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01[/font].......................编写如下的命令流文件,并命名为acce.inp[font=宋体]*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行[/font][font=宋体](3e16.6) !03行[/font][font=宋体]*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行[/font][font=宋体]ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行,同上ACCEXYZ(0,3)=3 !08行,同上[/font]finish[font=宋体]/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itime[/font][font=宋体][font=宋体]AUTOTS,0NSUBST,1, , ,1KBC,1[/font]acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO[/font]最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
地震分析算例_ANSYS
地震分析算例_ANSYS地震分析是指通过数值模拟和分析地震过程及其对结构物的影响,以评估结构物在地震中的性能和安全性。
在这个算例中,我们将使用ANSYS 软件进行地震分析,分析一个简单的二维框架结构在地震中的响应。
下面是详细的步骤和算例设置:1.几何建模:我们首先在ANSYS中进行几何建模,绘制一个二维的框架结构。
框架由4个节点组成,其中1号和4号节点是固定支座,2号和3号节点是自由节点。
我们可以设置框架的长度、宽度、高度等参数。
2.材料属性:我们需要为框架结构定义材料属性,包括弹性模量和泊松比等。
这些参数可以根据实际的材料特性进行设置。
3.边界条件:我们将1号和4号节点设置为固定支座,以防止结构物在地震中发生位移。
4.地震负荷:我们需要定义地震负荷,即地震的加速度记录。
这些加速度记录可以根据地震现场的实测数据来确定。
在ANSYS中,可以将地震加速度记录分为不同的时程,并将其作为负荷应用在结构上。
5.模型分析:在所有参数设置好后,我们可以进行模型分析。
在ANSYS中,可以选择静力分析或动力分析进行地震分析。
如果选择静力分析,将根据结构物的初始状态和地震负荷计算结构物的响应。
如果选择动力分析,则可以考虑结构物的动态特性和阻尼效应。
6.结果评估:结果评估可以包括结构物的最大位移、最大应力、最大应变等信息,以及结构物的破坏模式和安全性评估等。
在ANSYS中,可以通过可视化和结果输出等方式来进行结果评估。
总结:在这个地震分析算例中,我们使用ANSYS软件对一个二维框架结构进行了地震响应的模拟和分析。
通过设置几何模型、材料属性、边界条件和地震负荷等参数,进行模型分析并评估结构物的性能和安全性。
通过这个算例,我们可以更好地理解地震分析的过程和方法,并为实际工程项目提供参考和指导。
ANSYS时程分析-考虑地震作用的建筑物加速度瞬态分析
! 【算例】考虑地震作用的建筑物加速度瞬态分析C*** 采用等效实体单元进行建筑物的地震加速度瞬态分析c*** 进行中央电视台新址主楼的等效几何建模*** begin ***FINISH/CLEARc*** 设置参数L=150W=150H=235PI=3.1415926 ! 底座长度! 底座宽度! 整个建筑的高度A=PI*6/180L_HL=L/2T_HL=L/2-TAN(A)*H LO_H=40TO_H=75WW=40/PREP7! 内倾角! 底座长度的1/2 ,作为角点坐标! 顶部长度的1/2 ,作为焦点的坐标! 底座的高度! 顶部悬空部分的高度! “腿”部的宽度/TITLE, EX 9.2(12) by Zeng P, Lei L P, Fang GET,1,SOLID45MP,EX,1,1e10MP,PRXY,1,0.23 MP,DENS,1,1e3 MP,DAMP,1,0.05 ! 几何建模K,,-L_HL,-L_HL, K,,L_HL,-L_HL, K,,L_HL,L_HL, K,,-L_HL,L_HL, K,,-T_HL,-T_HL,H,K,,T_HL,-T_HL,H, K,,T_HL,T_HL,H, K,,-T_HL,T_HL,H,V,1,2,3,4,5,6,7,8 ! 定义单元和材料参数! 由于采用了的简化实体模型,材料参数也进行了等效! 定义8 个角点! 根据8 个角点生成方锥体V1VGEN,2,1, , ,WW,WW,LO_H,, ! V2 (THE COPY OF V1) VGEN,2,1, , ,-WW,-WW,-LO_H,, ! V3 (THE COPY OF V1)VSBV,1,2 VSBV,4,3! V4 (V4=V1-V2) 切掉一个角! NEW_V1 (NEW_V1=V4-V3) 切掉另一个角c*** 为了实现六面体网格划分,需要对实体进行分块处理CSKP,11,0,14,34,10,1,1, ! 通过 3 个点建立局部坐标系WPCSYS,-1,11 VSBW,ALL! 基于局部坐标系建立新的工作面( workplane ) ! 用当前的工作面剖分实体CSKP,12,0,18,30,34,1,1, ! 以下的命令流与前三条的作用相同WPCSYS,-1,12 $VSBW,ALL $CSKP,13,0,27,26,25,1,1,WPCSYS,-1,13 $VSBW,ALL $CSKP,14,0,32,33,34,1,1,WPCSYS,-1,14 $VSBW,ALL $CSKP,15,0,15,31,32,1,1,WPCSYS,-1,15 $VSBW,ALL $CSKP,16,0,29,28,25,1,1,WPCSYS,-1,16 $VSBW,ALLESIZE,L/10! 定义单元尺度 VMESH,ALL ! 对体划分网格CSDELE,11,16,1, ! 删除前面建立的局部坐标系ASEL,S,LOC,Z,,, ! 选择 z=0 处的面,DA,ALL,ALL,, ! 对所选择的面施加完全的约束,即将地面固定 ALLSELFINISHc*** 进行中央电视台新址主楼的等效几何建模*** end *** c*** 进行瞬态分析/SOLU*DIM,ACC,ARRAY,6,4, ! 定义数组,在后面的语句中分别为时间和三个方! 向的加速度赋值(天津波简化数据)ACC(1,1)=0.1 $ACC(1,2)= 0.03599 $ACC(1,3)= 0.24653 $ACC(1,4)= -0.24363ACC(2,1)=0.7 $ACC(2,2)= 0.0293 $ACC(2,3)= 0.14177 $ACC(2,4)= -0.31038ACC(3,1)=1.2 $ACC(3,2)= 0.03695 $ACC(3,3)= 0.71123 $ACC(3,4)= -0.48752ACC(4,1)=1.7 $ACC(4,2)= -0.00725 $ACC(4,3)= -0.2167 $ACC(4,4)= -0.10452ACC(5,1)=2.4 $ACC(5,2)= 0.00734 $ACC(5,3)= -0.81282 $ACC(5,4)= 0.66357ACC(6,1)=3.0 $ACC(6,2)= 0.01759 $ACC(6,3)= -0.00823 $ACC(6,4)= -0.0361ANTYPE,TRANS ! 设置瞬态分析方式TIME,0.001 ! 设置初始加载步,设一个很小的数,但不能为零ACEL,0,0,0 !三个方向都施加零加速度,作为初始计算 SOLVE*DO,TT,1,6,1 !根据不同的时间步,施加加速度给整个结构 TIME,ACC(TT,1)ACEL,ACC(TT,2),ACC(TT,3),ACC(TT,4) ! 三个方向都施加加速度SOLVE *ENDDOc*** 进入时间历程后处理器 /POST26 ! 进入时间 -历程后处理 /AXLAB,x,Time(s) ! /AXLAB,y,Displacement (m) ! NSOL,2,node(-T_HL,-T_HL,H),u,X NSOL,3,node(-T_HL,-T_HL,H),u,Y NSOL,4,node(-T_HL,-T_HL,H),u,Z PLVAR,2,3,4 ! 输出该点位移随时间的变化关系曲线/AXLAB,y,von Mises Equivalent Stress (Pa) ! 定义 y 轴标题ANSOL,5,node(-L_HL+LO_H,L_HL-LO_H*tan(A),LO_H),s,eqv !定义" 腿部 " 角点(受力最!大)等效应力为 y 轴的变量plvar,5 ! 输出该点等效应力随时间的变化关系曲线FINISH!%%%%% %%%% end %%%%%%%%%%%定义 x 轴标识定义 y 轴标识!将建筑物悬空部最高点的位移设置为 2 号变量欢迎您的下载,资料仅供参考!致力为企业和个人提供合同协议,策划案计划书,学习资料等等打造全网一站式需求。
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。
在ansys中如何施加地震波
三向输入简化后的单向输入首先,将三个方向的地震加速度放到一个文本文件里,如accex yz.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-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件,并命名为ac ce.inp*dim,ACCEXY Z,TABLE,2000,3 !01行*vread,A CCEXY Z(1,1),accexy z,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,A CCEXY Z(1,0),time,txt !04行(e16.6) !05行ACCEXY Z(0,1)=1 !06行XYZ(0,2)=2 !07行,同上ACCEXY Z(0,3)=3 !08行,同上sh/SOLU9ANTYPE,transbtime=0.01 !定义计算起始时间e=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtime,itimeTS,0NSUBST,1, , ,1KBC,1acel,ACCEXY Z(itime,1),ACCEXY Z(itime,2),ACCEXY Z(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
在 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-010.200000E-01 0.300000E-01 0.400000E-010.500000E-010.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/SOLUANTYPE,transbtime=0.01!定义计算起始时间etime=15.00!定义计算结束时间dtime=0.01!定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ<itime,1>,ACCEXYZ<itime,2>,ACCEXYZ<itime,3>!施加三个方向的地震加速度SOLVE*ENDDO最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
三向输入简化后的单向输入首先,将三个方向的地震加速度放到一个文本文件里,如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-010.200000E-010.300000E-010.400000E-010.500000E-010.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/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。
说明和讨论:1、命令流中各行说明:01行:定义2000行,3列的数组,(行数根据数据文件里加速度点数来定)ACCEXYZ用来存放三个方向的加速度值。
02行:从数据文件accexyz.txt里读加速度值到数组ACCEXYZ,2000为行数,可根据情况修改。
其中的JIK,3,2000非常重要,它决定着将从加速度数据文件中数据输入到ACCEXYZ数组时的存放格式,这里用到按JIK方式变化,即读进来的数据依次放入ACCEXYZ(I,J)(J从1到3,I从1到2000,K默认从1到1)。
根据这行命令下面(3e16.6)的格式,每次从数据文件里读一行三个数据,分别放入ACCEXYZ(1,1)、ACCEXYZ(1,2)、ACCEXYZ(1,3),接着再读下一行的三个数据,分别放入ACCEXYZ(2,1)、ACCEXYZ(2,2)、ACCEXYZ(2,3),依次类推。
当然依据数据文件的格式,也可以选用IJK, IKJ, JIK, JKI, KIJ, KJI等格式,其中IJK为默认。
03行:读数据的格式,每行三个数值。
由于数据文件中的数据是用科学记数年法表示的,因此,这里也用相应FORTRAN的科学记数法的格式。
如果数据文件里的数值是如“0.2876 1.2333 2.9938”这样的格式,此行的数据格式也就修改为“(3f10.4)”这样的格式。
04行:从数据文件time.txt里读时间值到数组ACCEXYZ的第零列。
06行:将数组ACCEXYZ的第零行赋值,如果不对行插值的话也可以不赋值对于地震波的输入,可以把荷载记录做成文件,利用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),ACETT,,,1(e17.6)ACEX(0,1)=1*end/input,ff*creat,ff*vread,aceY(1,1),acey,txt,,1(e16.6)*vread,aceY(1,0),ACETT,,,1(e17.6)ACEY(0,1)=1*end/input,ff*creat,ff*vread,aceZ(1,1),acez,txt,,1(e16.6)*vread,aceZ(1,0),ACETT,,,1(e17.6)ACEZ(0,1)=1*end/input,ff!地震波时程记录分成了3个文件,每个文件是一列。
分别记录x,y,z方向的加速度。
acett 是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件,这样程序就简单多了。
大家可以试看看修改一下。
下面是计算部分语句:/SOLUANTYPE,trans!求解其自己选了TM_START=0.01TM_END=15.00TM_INCR=0.01*DO,TM,TM_START,TM_END,TM_INCRTIME,tmalpha,BETAD,ACEL,acex(tm),acey(tm),acez(tm)SOLVE*ENDDOfini!借助于ANSYS动力分析基本功能和APDL语言以及有关的地震资料编写的SEISMIC程序NT=100DT=0.02*DIM,ac,,NT/INPUT,FF,txt/SOLUNSUBST,1,,,1OUTRES,ALL,1ANTYPE,TRANS*DO,I,1.NTACEL,0.01*ac(I)*1.3,0,0 !ac(I)的值由地震资料获得TIME,I*DTSOLVE*ENDDO地震分析算例(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(8) 接着建立梁单元的几何属性,和上一个例子一样,采用Sections建模,进入ANSYS 主菜单Preprocessor->Sections->Beam->Common Sections,选择Sub-Type为工字型,截面尺寸W1=0.2,W2=0.2, W3=0.5, t1=0.01, t2=0.01, t3=0.008(9) 通过实参数输入集中质量单元的质量和转动惯量,在ANSYS主菜单中选择Preprocessor->Real Constants->Add/Edit/Delete菜单,在Real Constants窗口中选择Add,在Element type for Real Constants选择Mass 21,在Real Constant for Number 1窗口中输入1.6E2, 1.6E2,如图所示。
即该质量单元在X和Y方向的质量都为160,由于本例子模型为平面问题,所以不必考虑Z方向的质量,同样也不考虑单元的转动惯量。
(10) 继续添加第二类集中质量,过程和上面一样,但是输入的质量数值为1.2E2, 1.2E2(11) 完成以上工作就完成了模型的基本数据准备,下面开始建立物理模型。
(12) 在ANSYS主菜单中选择Preprocessor->Modeling->Create->Keypoints->In Active CS,依次输入关键点编号和坐标:需要说明的是,关键点6为后面建立梁单元所需的截面方向控制点,在上一个例子中已经做过介绍。
(13) 完成关键点输入后下面建立直线模型。
在ANSYS主菜单中选择Preprocessor->Modeling->Create->Lines->Lines->Straight Line,依次连接关键点1~5。
(14) 下面给建立完的几何模型赋予材料属性,在ANSYS菜单中选择Meshing->Mesh Attributes->Picked Lines,选中所有的直线,进入Line Attributes窗口,选择相关选项材料属性,实参数,单元类型和截面类型都为1,点击OK后输入关键点6作为截面方向控制点。