基于ANSYS的巨型框架结构体系地震响应分析

  • 格式:pdf
  • 大小:446.48 KB
  • 文档页数:4

第34卷第4期 2011年4月合肥工业大学学报(自然科学版)J OU RNAL OF H EFEI UN IV ERSIT Y OF TECHNOLO GYVol.34No.4 Apr.2011 收稿日期:2010203210;修回日期:2011203214基金项目:国家自然科学基金重大研究计划资助项目(90715016)作者简介:花能新(1986-),男,安徽肥西人,合肥工业大学硕士生;叶献国(1954-),男,安徽无为人,博士,合肥工业大学教授,博士生导师.Doi :10.3969/j.issn.100325060.2011.04.016基于AN S YS 的巨型框架结构体系地震响应分析花能新, 叶献国(合肥工业大学土木与水利工程学院,安徽合肥 230009)摘 要:文章采用有限元软件ANSYS 自带的A PDL 语言进行参数化三维结构建模,对巨型框架结构进行抗震性能计算分析,包括模态分析、谱分析及时程分析,讨论了不同主次框架刚度比下巨型框架结构的动力特性及地震响应。

结果表明:钢筋混凝土巨型框架在地震作用下的侧移曲线为弯剪型;主次框架刚度比对巨型框架弹性地震反应影响较大;当结构高阶频率对应周期与地震波卓越周期相近时有可能发生类共振现象。

关键词:巨型框架;A PDL 语言;模态分析;谱分析;动力时程分析中图分类号:TU31113 文献标识码:A 文章编号:100325060(2011)0420545204Analysis of seismic response of RC mega 2frame structures by ANSYSHUA Neng 2xin , YE Xian 2guo(School of Civil and Hydraulic Engineering ,Hefei University of Technology ,Hefei 230009,China )Abstract :U sing t he ANS YS paramet ric design language (A PDL ),a t hree 2dimensional parametric finite element model is set up in t his paper.The aseismatic performance of t he mega 2f rame st ruct ure is dis 2cussed including t he modal analysis ,spect ral analysis ,and time 2history analysis.And t he dynamic characteristic and seismic response of t he mega 2frame st ruct ure wit h different stiff ness ratio between t he major and minor frames are st udied.The result s show t hat t he displacement curve of reinforced concrete (RC )mega 2f rame struct ures is a type of flexure 2shear.The stiff ness ratio between t he major and minor f rames has a notable influence on t he elastic seismic response of t he st ruct ure.The higher order frequency of t he struct ure has certain influence on t he st ruct ure ,and when it s period is close to t he eart hquake wave p redominant one ,a p seudo 2resonance p henomenon may occur.K ey w ords :mega 2f rame st ruct ure ;ANS YS parametric design language (A PDL );modal analysis ;spec 2t ral analysis ;time 2history analysis 钢筋混凝土巨型框架结构由主次2级结构组成,它能够充分发挥材料和结构性能,开展巨型框架结构体系的研究,有重要的理论价值和实际意义[1,2]。

本文采用有限元软件ANS YS 自带的A PDL 语言进行三维结构建模,对巨型框架结构进行了抗震性能的计算分析,包括模态分析、谱分析和时程分析。

考虑到影响结构弹性地震反应的因素有很多,如场地条件、地震烈度、地震波的频谱特性和结构本身的特性如刚度、质量、阻尼、周期等,本文对巨型框架结构在不同主框架和次框架刚度比下进行了弹性地震响应动力时程分析。

1 巨型框架结构有限元模型模型共55层,总高度201m ,高宽比为5158。

普通框架层高316m ,巨型层层高412m ,第11、22、33、44及55层为巨型层。

模型平面尺寸为(319m ×10)×(316m ×10),四角设有718m ×712m 的筒体巨型柱,各筒四角设有方钢管柱。

筒体墙厚度由高到低从40cm递增至60cm,混凝土标号从C30递增至C60。

巨型层上下楼板处设交叉钢骨梁,节点处设方钢管柱,并设柱间双向交叉支撑。

巨型层柱为方钢管混凝土柱(40mm×40mm×500mm×500mm),梁为工字形钢骨梁(300mm×650mm2 40mm×400mm×150mm×40mm)。

巨型层混凝土标号均为C30(包括钢管支撑、钢管梁及钢管柱,不包括角筒)。

普通框架层x、y向主梁截面尺寸分别为300mm×650mm、300mm×600mm。

巨型层板厚20cm,普通层板厚10cm,混凝土标号均为C30。

巨型框架7度设防,处于Ⅱ类场地,第1设计分组[2]。

在结构有限元建模中,普通层框架梁、柱和巨型层梁柱以及斜撑均采用平面梁单元(Beam188)进行模拟,筒体墙以及楼板采用壳单元(Shell63)模拟。

Beam188梁单元最大的特点是支持梁截面形状显示,可以考虑剪切变形和翘曲,同时也支持大转动和大应变等非线性行为,而且可以直接显示梁截面上的应力和变形。

Shell63单元可以同时承受壳面内、外的荷载,具有壳单元算法和膜单元算法,而且还考虑了应力刚化效应以及大变形效应。

壳单元具有平面内和平面外刚度。

梁单元与壳单元直接连接,连接方式默认为刚接[3]。

有限元模型如图1所示。

整个结构有限元计算模型共包括板壳单元24903个,梁柱单元32376个;共计单元59307个,共计节点65130个。

图1 巨型框架有限元模型2 模态分析及谱分析对与地面接触的框架柱和筒体施加x、y、z 方向的零位移约束。

本例采用Subspace法提取前30阶模态,模态分析采用扩展模态法[4]。

主要计算结果见表1所列。

表1 结构主振型频率阶次固有频率/Hz周期/s振 型1012358741239624结构在y方向上平动2012435241106439结构在x方向上平动3015166711935471结构绕z轴转动5017729911293678结构在y方向上平动加强,振型出现S型7017730411293594结构在z方向上振动,项部出现变形 以平动为主的第1振型周期T p=4124s,以扭动为主的第1振型周期T n=1194s,其比值T n/T p=01457。

结构抗扭刚度满足规范要求,扭转效应对结构的影响不明显,另外,5阶以后的高阶振型差别不明显。

谱分析是一种将模态分析的结果与一个已知的谱联系起来,计算模型的位移和应力的分析技术。

本文主要用振型分解反应谱法来确定结构对地震荷载的动力响应[4,5]。

巨型框架所处Ⅱ类场地,抗震设防烈度为7度,在模态分析的基础上,对模型施加x向的7度多遇加速度响应谱,采用SRSS模态合并方法,综合前10阶振型,通过选择合并的输出方式,得到楼层的最大位移响应值为4146cm,满足规范要求。

整个模型的位移与角位移云图如图2所示。

图2 反应谱分析结构位移与角位移云图645 合肥工业大学学报(自然科学版)第34卷 3 地震响应分析对巨型框架进行动力时程分析,地震波选用El2Cent ro(N2S)波和Taft(S69E)波,峰值加速度分别为34117cm/s2和17519cm/s2,时间间隔0102s。

依据文献[6],在进行时程分析时,将加速度峰值调整为35cm/s2,持时30s,相当于7度多遇地震[7]。

结构的阻尼比取0105。

地震波为x方向单向输入,可以得到巨型框架结构在2种地震波作用下楼层绝对位移曲线以及层间位移角曲线,如图3、图4所示。

从图3、图4可以看出,楼层最大水平位移沿高度不断增加,其侧移曲线呈弯剪型,结构的最大层间位移角基本均出现在第9层,这与楼层的刚度突变有关,但都远小于弹性层间位移角限值。

也就是说,第9层为整个结构的薄弱层,在设计时应当加以注意。

参考文献[2]的相关定义,令:β=主框架刚度次框架刚度=∑i a+∑i b ∑i c+∑i d,其中,Σi a、Σi b分别为主框架梁、柱的线刚度和;Σi c、Σi d分别为次框架梁、柱线刚度和。

梁、柱线刚度和分别以小层计算[2]。

该模型的主次框架线刚度比β约为797159。

保持次框架梁、柱截面不变,通过变化主框架梁、柱截面而改变β,共得到6种方案,对6种方案分别进行模态分析,根据所得各阶模态振型,得出每种方案的自振周期,见表2所列。

图3 楼层绝对位移曲线图4 楼层层间位移角变化曲线选取上述2种地震波,对不同方案结构分别进行地震响应计算和分析。

各方案楼层顶点最大位移与加速度反应见表3所列。

表2 不同方案的β值及结构自振周期参数方案1方案2方案3方案4方案5方案6β16210827211647815579715913811361594132周期/s510141724133412441033199表3 不同方案楼层顶点最大位移和最大加速度比较方案El2Centro波s/cm a/(m・s-2)Taft波s/cm a/(m・s-2)方案14186013741510143方案23199013851340148方案33164014641970154方案43158014641420151方案53145014831410155方案63147014931220153 从表3可以看出,同一种结构方案在不同的地震波作用下其响应不同。