基于ANSYS的重力坝抗震性能分析
- 格式:doc
- 大小:268.00 KB
- 文档页数:8
ANSYS_地震分析算例地震是指地球上因地壳运动而产生的震动现象。
在地震发生后,建筑物的结构稳定性和抗震性能至关重要,因为地震可以对建筑物造成严重破坏。
因此,在建设和设计建筑物时,地震分析变得非常重要。
在此我将介绍一种用ANSYS进行地震分析的算例。
在地震分析中,我们首先需要建立一个合适的模型。
在这个算例中,我们将使用ANSYS提供的有限元分析方法。
首先,我们需要创建一个建筑物的三维模型。
在建筑物的模型中,我们需要包括所有的结构细节,例如建筑物的基础、柱子、梁和地板等。
我们可以使用ANSYS的几何建模工具来创建这个模型。
接下来,我们需要为建筑物定义材料特性。
建筑物的材料特性会对地震分析的结果产生重要影响。
例如,不同种类的混凝土、钢铁和木材等材料在地震作用下的响应是不同的。
我们需要使用ANSYS的材料库来定义这些材料的特性。
完成模型和材料定义后,我们需要定义地震荷载。
地震荷载是指地震发生时作用在建筑结构上的力量。
地震荷载可以根据地震的震级和地震波的性质来确定。
我们可以使用ANSYS的预处理工具来定义这些地震荷载。
接下来,我们需要定义边界条件。
边界条件是指建筑物与外部环境之间的约束关系。
例如,建筑物的基础是固定的,地震发生时不能移动。
我们需要使用ANSYS的加载工具来定义这些边界条件。
完成了模型、材料、地震荷载和边界条件的定义后,我们可以进行地震分析。
地震分析是指通过模拟地震发生时结构的动力响应来评估建筑物的抗震性能。
在ANSYS中,我们可以使用动力分析工具来进行这个分析。
在地震分析过程中,我们可以观察到各个部位的应力和位移等响应。
这些响应可以帮助我们评估建筑物的破坏机制和结构的安全性能。
例如,我们可以观察到柱子是否出现弯曲、梁是否发生裂缝等。
根据地震分析的结果,我们可以对建筑物的设计和结构进行优化。
例如,我们可以调整柱子和梁的尺寸、材质和布置方式,以提高建筑物的抗震能力。
综上所述,通过ANSYS进行地震分析可以帮助我们评估和优化建筑物的抗震性能。
基于ANSYS的高层钢结构抗震及稳定性分析共3篇基于ANSYS的高层钢结构抗震及稳定性分析1基于ANSYS的高层钢结构抗震及稳定性分析随着城市化进程的不断加快,建筑高度和层数不断增加,高层建筑的结构安全问题越来越受到人们的关注。
而地震是高层建筑结构安全的关键因素之一,抗震设计成为高层建筑结构设计的重点之一。
而对于钢结构而言,钢材的高强度、可塑性好、适应性强等特点,使得钢结构成为高层建筑结构的重要选择。
本文将以基于ANSYS的高层钢结构为对象,探讨其抗震及稳定性分析。
1. 建立高层钢结构有限元模型在进行高层钢结构的抗震及稳定性分析前,需要先通过ANSYS 等有限元软件建立高层钢结构的有限元模型。
建立模型需要考虑高层钢结构的结构特点和工程实际情况,确定结构参数、节点分布及约束情况。
2. 高层钢结构抗震分析地震对高层建筑结构的影响主要体现在地震作用下建筑结构内部产生的地震应力和滞回曲线等。
因此,在进行高层钢结构的抗震分析时,需要考虑其受到的地震作用,分析结构内力和变形等参数。
首先,需要进行地震作用下钢结构模型的动力特性分析。
在这一步中,可以使用ANSYS中的模态分析功能,以得到结构在不同模态下的自然频率和振型。
其次,根据钢结构在地震作用下的动力特性,进行地震反应谱法抗力设计。
地震反应谱是描述结构在不同频率下受到地震作用时的反应的一种方法,可以分析结构受到的地震作用下的最大位移、加速度和力等参数。
对于高层钢结构,可以通过ANSYS中的响应谱分析功能进行计算。
最后,通过引入钢结构弹塑性性能纳入分析中,能够更加精准地分析高层钢结构在地震作用下的受力性能。
3. 高层钢结构稳定性分析高层钢结构的稳定性是结构设计或构件设计中必须考虑的重要问题。
高层钢结构结构体系复杂,其极限状态的稳定性较低。
在进行高层钢结构的稳定性分析时,需对结构进行屈曲分析,以了解梁和柱在地震作用下的稳定性。
在进行屈曲分析时,需要先得到高层钢结构构件的稳定系数。
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的重力坝三维静动态结构分析目录1 引言 (1)2 工程概况 (1)3 基本资料 (1)3.1 反应谱 (1)3.2 材料参数 (2)3.3 规范要求 (2)4 分析简介 (4)4.1 分析模型 (4)4.2 边界条件 (6)4.3 荷载工况 (6)5 计算成果 (7)5.1 工况一 (7)5.2 工况二 (8)5.3 工况三 (10)5.4 工况四 (11)5.5 工况五 (12)5.6 工况六 (14)5.7 结果总结及分析 (15)6 结论及建议 (17)7 分析命令流 (17)1 引言重力坝是我国高坝中的主要坝型,在防洪、发电、灌溉、城镇供水、航运、养殖和旅游等方面发挥了巨大的作用,取得了显著的经济效益和社会效益。
众所周知,重力坝主要依靠其自身的重力来维持稳定,其坝体体积大,稳定性好。
但由于各种原因,仍有可能失事。
因此,重力坝的应力应变状态和坝基稳定性一直都是设计和施工十分重视的问题。
此外,大坝多建于地震频发的地区,因而对重力坝进行地震荷载作用下的安全评估也十分必要。
本次作业采用有限元方法,运用大型通用有限元分析软件ANSYS,对简化的三维重力坝的线弹性模型在静动力工况下进行有限元计算,并对结果加以分析,最后给出安全评价结论及建议。
2 工程概况某水电站是以发电为主,兼有防洪,航运等综合效益的水电枢纽工程。
该工程枢纽总体布置采用砼重力坝挡水,大坝基本坝剖面为上游坝坡铅直,下游坝坡为1:0.75。
坝顶总长270m,坝高180m,坝顶宽18m,坝底宽139.5m,正常蓄水位170m。
重力坝坝低至坝高100m之间使用坝体混凝土Ⅱ,坝高100m至坝顶之间使用坝体混凝土Ⅰ。
上游正常蓄水位为170m,下游无水。
3 基本资料3.1 反应谱谱分析是一种将模态分析的结果与一个已知的谱联系起来计算结构的位移和应力的分析技术。
在土木工程动力响应分析中,谱分析代替时间-历程分析,特别是抗震分析,主要用来确定结构对随机荷载或随时间变化荷载的动力响应。
在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对于坝体的研究分析报告坝体及相关建筑在使用过程中,会承受如重力、净水压力、淤泥荷载、浪压力、扬压力等各种作用,而我们在设计、建造这个建筑之前,要分析其产生的应力、应变进而选取材料和校核材料的安全性。
为分析所需,基于Ansys软件建立相应的模型,并施加荷载和作用,在三种工况下校核结构的安全性。
一:分析对象1:坝体的几何参数:2:基岩的几何参数:二:作用及荷载(1) 约束基岩左右两端受x 方向的位移约束,基岩下端受x 、y 两个方向的位移约束。
(2)静水压力正常蓄水位高程91.75m ,防洪高水位97m ,校核洪水位101m 。
对应下游水位分别为15m ,20m 和25m 。
(3)泥沙荷载坝前泥沙淤积高程: 25m 。
坝前泥沙浮容重:6.0kN/m 3,淤沙内摩擦角:12°。
坝面上单位宽度上的泥沙压力为221(45)22s sk sb s p h tg ϕγ=︒- 式中: sk p ——淤沙压力标准值(KN/m );sb γ——泥沙的浮容重,取6kN/m 3;s h ——泥沙淤积深度(m );s ϕ——淤沙的内摩擦角,12°。
(4)浪压力50年一遇计算风速21m/s ,多年平均最大风速14m/s ,有效吹程1km 。
(5)扬压力取渗透压力强度系数α=0.25,帷幕中心线坐标X=10m 。
三:选用单元及划分网格1) 单元选择:Solid –Quad 4node 422) 材料参数:坝体和基岩分别设置,见上图。
3) 划分网格:坝体部分-外围线按1m 每格划分,整体按自由网格划分。
基岩部分-靠近坝体网格密集,坝基面水平线上基岩外围线按20份、4的比率划分;垂直地表的按20份、0.25的比率划分;基底均分。
整体基岩自由网格划分。
四:三种工况的具体Ansys设置1)正常蓄水位(其中括号内为承载能力极限状态时的分项系数)上游水位高程为91.75m,下游水位高程为15m。
(1)上下游静水压力(分项系数为1.0)gradient 斜率为-9810 沿y轴方向,分别取91.75m和15m在各自位置。
ANSYS地震力分析之时程分析我们知道反应谱分析法是现代抗震设计的基本理论,它能反映出结构在不同自振周期和阻尼比下的最大反应,也就是在给定地震加速度时间过程下,粘滞阻尼体系的最大反应相对于自振周期和阻尼系数的函数关系。
但是反应谱分析中要用来地震参数,如地震动反应谱值受人为因素影响太大。
而且,反应谱分析仅能给出结构各振型反应的最大值,而丢失了与最大值有关且对振型组合又非常重要的信息,如结构位移(或应力、弯矩)最大值发生的时间及其正负号,使得各振型最大值的组合陷入困境。
所以,反应谱方法在实际分析中依然存在较大的误差。
国外相关试验研究也做出了总结,认为:1)模态分析可以获得比较准确的自然频率和振型;2)对于应力应变只能定性的分析应力危险区和应变过大区域,不能定量的分析;但是可以为产品的初期设计提供改善依据和参考;利用ANSYS计算地震力除了反应谱分析法之外还有时间历程响应分析法。
时间历程响应理论的分析方法是通过输入对应于工程场地的若干条地震加速度记录或人工加速度时程曲线,通过积分运算求得在地面加速度随时间变化期间结构的内力和变形状态随时间变化的全过程,并以此进行结构构件的界面抗震承载力验算和变形验算。
而且时间历程方法的另一个特点就是能够进行非线性的动力学性能分析,虽然在分析计算中将消耗大量的时间,但弥补了反应谱理论的缺陷和不足。
在框架地震力分析中,步骤分为以下几个步骤:1)建立模型2)加边界条件3)静力分析4)模态分析5)地震时程分析一、静力分析①施加重力加速度acel,0.-9.8②施加均布于梁上的压力sfbeam,,2,pres,150000,150000③结果后处理etable,saxl,ls,1建立轴力单元结果表格etable,epelaxl,lepel,1建立轴应变单元结果表格二、模态分析antype,modalmodopt,lanb,20,0,100,,offmxpand,20,,,yesmodopt命令,用来定义模态分析选项,参数为MODOPT, Method, NMODE, FREQB, FREQE, Cpxmod/PRMODE, Nrmkey, --, BlockSizelanb是模态提取方法选项,20表示模态扩展数量,0表示自己感兴趣的频率范围,在这里表示频率范围为0~100Hz,freqb开始频率其默认值,程序会自动计算,freqe其默认值为1e8,Nrmkey表示是模态振型是否进行正则化处理,其默认值为对质量矩阵的正则化,mxpand命令,用来定义模态扩展数MXPANDMXPAND, NMODE, FREQB, FREQE, Elcalc, SIGNIF, MSUPkey Nmode的值设为-1时,表示不扩展,如果是0则扩展所有Elcalc表示是否计算单元结果和反力,默认是不计算,Msupkey表示是否将单元结果写入模态分析结果中,当elcalc设置为yes 后,msupkey默认值为yes,会对单元结果写入模态分析结果中,用于后面的模态叠加三、地震时程分析因为时程分析主要了解结构在任意载荷作用下下的动力响应,所以还要考虑阻尼效应antype,transalphad,betad,timint,off !关闭时间积分nlgeom,on !大变形选项打开time,1e-6kbc,1acel,0,-9.8solve由于框架自身重力已经对结构产生了影响,设置一个分析时间为1e-6的瞬态分析,为下面的地震力分析,得到一个框架的预应力和预变形。
基于ANSYS的重力坝抗震性能分析
【摘要】建立一个120m重力坝模型,利用ANSYS分析软件,分析此重力坝挡水坝段在静,动力作用下应力变化规律,并对坝体的抗震安全性能进行评估,为
类似工程设计、施工提供理论依据。
【关键词】重力坝;ANSYS;反应谱;地震
重力坝是世界上最早出现的一种坝型之一。
依据其相对安全可靠,耐久性好,对
不同的地形和地质条件适应性强等特点,重力坝在各个国家都很流行。
由于重力
坝大多都建在高烈度或地震多发地区,一旦失事,损失不可估量,因此在大坝时
对其进行抗震安全分析十分必要。
ANALYSIS OF SEISMIC PERFORMANCE OF GRAVITY DAM
BASED ON ANSYS
【Abstract】Establish a 120m gravity dam model and using ANSYS analysis software, analysis of the gravity dam retaining dam in static and dynamic effect of the stress change rules, and on the dam seismic safety performance assessment, to provide a theoretical basis for the design and construction of similar projects.
【Keywords】gravity dam;Ansys;response spectrum;earthquake 1 有限元模型
1.1 计算基本假定
(1)假定库水为不可压缩流体,库水对坝体的动力相互作用以坝面附加质量的
形式计入;
(2)坝体材料假定为线弹性,并假定不同部位材料有不同的弹性常数;
(3)采用无质量地基方案,近似考虑坝体结构和地基间的动力相互作用;
(4)地基为均匀弹性体,并于坝体紧密联系在一起。
1.2 有限元模型
设计某重力坝,坝高120m,坝底宽为76m,坝顶宽为10m,上游坝面坡度和
下游坝面坡度如图1所示。
因为设计此重力坝结构比较简单,垂直于长度方向的断面结构受力分布情况也基本相同,且大坝的纵向长度远远大于其横断面,因此大坝的抗震性能分析选用单位断面进行平面应变分析是可行的。
图1 重力坝断面结构图
大坝抗震性能分析的计算条件如下:
(1)重力坝基础按嵌入到基岩中考虑,假定地基为刚性。
(2)重力坝计算采用材料参数如下:弹性模量E=35GPa,泊松比ν=0.2,容重γ=25KN/m³。
(3)计算分析大坝水位为120m。
(4)水的密度为1000Kg/m³。
(5)重力坝抗震设防地震烈度为8度,水平方向地震加速度值为0.2g。
2 有限元计算
2.1 计算方法分析
大坝的抗震性能分析主要有三种方法: 拟静力法、时程分析法和反应谱方法。
拟静力法是一种把地震的影响用一种折算的静荷载来表示,求出这种地震荷载后,按照常规的静力法进行坝体的各项应力、位移的抗震分析方法。
它是假定地震时与地面加速度相同的加速度作用在坝体各部位,求出地震时的惯性力,然后根据惯性力来评价大坝的安全性。
时程分析法是将地震动记录或人工波作用在结构上,直接对结构运动方程进行积分,求得结构任意时刻地震反应的分析方法,所以动态时程分析方法也称为直接积分法。
反应谱分析法是以单质点弹性体系在实际地震过程中的反应为基础,来进行结构反应的分析,它通过反应谱巧妙地将动力问题静力化,使得复杂的结构地震反应计算变得简单易行。
按照这一理论,应用地震谱曲线,就可以按照实际地面运动来计算建筑物的反应。
反应谱是单质点弹性体系对于实际地面运动的最大反应和体系自振周期的函数关系。
对于复杂的结构可以简化为若干振型的叠加,每个振型又可转化为一个单质点来考虑。
使用已经确定的设计反应谱计算重力坝在地震作用下的反应,就归结于寻求坝体的自振特性。
地震产生的破坏,与受力大小和受频谱最大振动的持续时间都有关系。
在进行反应谱分析计算之前,首先要计算大坝的自振特性。
模态分析是用于确定结构的振动特性,即结构的固有频率和振型。
它们是结构受动态荷载设计中的重要参数,也是更详细的动力分析的基础。
模态分析计算中采用子空间迭代法提取模态。
水深h 处的地震动水压力的作用按下面公式转化为相应的坝面附加质量。
h H a h P w h 0w 8
7)(ρ= (1) 根据图2所示的大坝设计的反应谱曲线图,可得大坝反应谱曲线方程:
⎪⎪⎪⎩⎪⎪⎪⎨⎧<*≤<≤<-+=T T T
T T T T T g g .............)(1.0..........................1.00).......1(101max 9.00max max ββββ (2)
图2 大坝设计反应谱曲线图
本次重力坝抗震性能分析中,m ax 取值为2,g T 取值为0.2,特征周期0T 取值为0.2S 。
2.2 静力分析
该重力坝静力分析应力应变云图见图3。
由图3可知,在水荷载的作用下,坝体水平方向最大位移发生在坝顶,其值为12.402mm ;坝体垂直方向最大位移发生在坝腹出,其值为0.439mm 。
(a )坝体X 方向位移云图 (b )坝体Y 方向位移云图
(c)X方向应力云图(d)Y方向应力云图(e)第一主应力云图
图3 大坝静力分析应力应变云图
2.3动力分析
绘制坝体振型图,如图4所示。
(a)第1振型(b)第3振型(c)第5振型(d)第7振型
(e)第9振型(f)第11振型(g)第15振型(h)第17振型
图4 坝体振型图
调出模态分析各阶频率,并由式(2)求出反应谱值,如表1所列。
表1 大坝动力计算18阶振动频率及反应谱值 Hz
振型振动频率1/s 振动周期T 反应谱值
1 3.4531 0.2896 1.433
2 7.9984 0.1250 2
3 10.889 0.0918 1.918
4 14.203 0.0704 1.704
5 21.64
6 0.0462 1.462
6 24.488 0.0408 1.408
7 30.586 0.0327 1.327
8 33.863 0.0295 1.295
9 35.599 0.0281 1.281
10 38.801 0.0258 1.258
11 40.335 0.0248 1.248
12 42.531 0.0235 1.235
13 47.751 0.0209 1.209
14 49.194 0.0203 1.203
15 52.122 0.0192 1.192
16 55.161 0.0181 1.181
17 55.505 0.0180 1.180
18 58.095 0.0172 1.172
因为响应谱分析是在频域内进行的,对于结构动力特性依赖于频率的变化而变化,因此在模态分析后要进行模态合并求解,才能得到坝体结构真实的总体效应。
合并模态后,得到坝体在各阶频率的真实位移云图和应力云图,如图5所示。
(a)(b)(c)(d)
(e)(f)(g)(h)
图5 坝体位移、应力云图
(a)第2阶X方向位移云图(b)第 2阶Y方向位移云图(c)第7阶X方向位移云图(d)第 7阶Y方向位移云图(e)第2阶X方向应力云图(f)第2阶Y 方向应力云图(g)第7阶X方向应力云图(h)第7阶Y方向应力云图
进行模态求解时并没有同时进行模态扩展,而是在获得谱解后又单独进行了扩展模态求解。
这样,有明显意义的模态为1、2、3、4、5和7阶模态。
图5仅提供了2、7阶位移、应力云图。
由图4、图5可以看出,在地震作用下,坝顶出现水平向最大位移,最大位移为1.394mm,主要由第2振型的贡献产生。
大坝第7阶X方向最大应力出现在下游约80m处,最大应力值为0.369Mpa,Y方向的最大应力出现在上游约100m 处,最大应力值为1.93Mpa。
总之,重力坝地震反应谱分析获得的各阶振型的最大、最小的应力和位移值的具体位置,可以通过ANSYS查询到,出现应力集中或超出规范规定的部位,在设计和施工时应给予足够的重视。
3 结论
通过对该重力坝的静、动力有限元分析,可以了解到大坝在静力作用下的应力分布情况和大坝在地震荷载作用下的动力响应特征,从而评价大坝在静、动力荷载作用下的安全性能。
静力分析表明:重力坝水平方向的最大位移出现在坝顶,在设计中强度验算应满足相关规范要求;动力分析表明:在地震作用下重力坝下游折坡处和80m高程处出现最大应力,在这些部位应采取一些结构措施或工程加固措施。
参考文献
[1]隋红军,上官子昌. 基于ANSYS重力坝抗震性能分析[J]. 建筑技术开发,2015,06:63-66.
[2]赵廷红,王新峰,梁俊国,梁俊青. 基于Ansys的龙滩碾压混凝土重力坝抗震动力分析[J]. 甘肃科学学报,2013,03:84-87.
[3] 李守华. 重力坝抗震性能的有限元分析[J]. 城市道桥与防洪. 2012(12)。