结构动力学分析.doc

MIDAS/GEN六层框架结构的动力分析工程概况建筑地点：北京市建筑类型：六层综合办公楼，框架填充墙结构。地质条件：根据设计任务说明地震设防烈度为8度。柱网与层高：本办公楼采用柱距为6.0m的内廊式小柱网，边跨为6.0m，中间跨为2.7m，层高取首层为4.5m,其余为3.3m，如下图所示：框架结构的计算简图： 梁、柱、板截面尺寸的初步确定：1、梁截面高度一般取梁跨度的1/12至1/8。本方案取1/106000=600mm，截面宽度取6001/2=250mm，可得梁的截面初步定为bh=250*600。楼板取120mm，楼梯板及休息平台板为100mm,平台梁250400。2、框架柱的截面尺寸梁截面尺寸（mm）混凝土等级横梁（bh）纵梁（bh）AB跨、CD跨BC跨C30250600250400250600柱截面尺寸（mm）层次混凝土等级bh1-6C30600600结构动力学分析用来求解随时间变化的载荷对结构或部件的影响。与静力分析不同，动力分析要考虑随时间变化的力载荷以及它对阻尼和惯性的影响。MIDAS/GEN可进行的结构动力学分析类型包括:瞬态动力学分析、模态分析、屈曲分析、动力非线性分析等。本文将以一个六层框架结构为例对结构进行模态分析和谱分析。一模态分析模态分析是用于确定设计中的结构或机器部件的振动特性。它也是其他更详细动力学分析的起点，例如瞬态动力学分析和谱分析等，可以通过模态分析确定结构部件的频率响应和模态。一般对于动力加载条件下的结构设计而言，频率响应和模态是非常重要的参数，即使在谱分析及瞬态分析中也是需要的。1.1动力学求解方法 MIDAS目前提供了三种特征值分析方法，它们是子空间法、分块Lanczos算法、多重Ritz向量法。本文采用子空间法进行计算求解。子空间法使用迭代技术，求出结构的前r阶振型，它内部使用广义Jacobi迭代算法。由于该方法采用了完整的质量和刚度矩阵，因此精度很高，但计算速度较慢，特别适用于大型对称特征值求解问题。分块Lanczos法特征值求解器采用Lanczos算法，Lanczos算法是用一组向量来实现递归计算。这种方法和子空间法一样精确，但速度较快。多重Ritz向量法以变分原理为基础，直接迭加法求出的是和激发荷载向量直接相关的振型，其收敛具有严格的理论基础，在物理和、力学的微分方程中占有很重要的位置，得到广泛的应用。1.2本工程模态分析结果1.2.1自振周期与振型: 使用MIDAS/GEN中的模态分析计算结构的自振周期和振型。模态分析所使用的方法是子空间迭代法。高层建筑结构振型多，分布规律很难掌握，扭转振动会对结构产生教大影响，因此不能简单的取前几阶进行计算。规范中规定对于高层结构一般取35阶振型。为使高层建筑的分析精度有所改进，其组合的振型个数适当增加。考虑到MIDAS/GEN软件的强大快速的数据处理能力和精度的要求，本文取30阶振型。从国内高层建筑结构设计经验来看，建议基本自振周期按以下的几个公式估计，其中N为地面以上建筑物结构层数。经验公式表达简单，使用方便，但比较粗糙，而且只有基本周期，但经常用于对理论计算值的计算与评价。 框架:T1=(0.080.1) N 框架一剪力墙:T1=(0.060.09) N钢结构:T1=0.1N本工程得经过MIDAS/GEN的计算得到固有周期、固有频率、振型参与质量等的数值结果;X方向振型参与达到总质量的95.57%, Y方向振型参与达到总质量的94.83%，经过整理取前十阶列表可得到表1表1 结构的自阵周期、频率、运动形态模态周期（s）频率（Hz）运动形态11.01490.9853Y方向平动20.97681.0236Z方向转动30.88261.133X方向平动40.30533.2757Y方向平动50.29363.406Z方向转动60.27163.6824X方向平动 由表1可知道结构以水平振动为主。第一阶振型为X方向的平动，第二阶振型为Y方向平动，第三阶振型为绕Z轴的转动。以平动为主的第一自振周期为T1= 1.0149s，以扭转为主的第一自振周期为T2=0.9768 s，其比值T2/T1=0.96满足高层建筑混凝土结构技术规程规定的0.90 （ A级高度）的严格限值，可见结构的平面刚度比较均匀。1.2.2前三阶模态图形:图1 第一阵型图 Y方向平动图2 第一阵型图 Y方向平动（侧面图）图3 第二阵型图 X方向平动图4 第二阵型图 X方向平动（侧面图）图3 第二阵型图 Z方向转动图4 第二阵型图 Y方向平动（平面图）由计算结果看出，结构的第一周期为0.08N 0.1N，即0.480.6，而计算结果为1.01s，超出了这个经验公式的区间内，在误差范围内属合理的周期值。经计算，结构的前三阶振型如图1图4所示。图中的各阶振型是从正视、侧视和俯视这三个角度观察的，以便一目了然。如果不给出别的视图，就会感觉到该振型在视角观察上不明晰。表2为计算的各个振型的参与系数模态号TRAN-X(%)TRAN-Y(%)ROTN-Z(%)1084.32750.1451200.144884.4453385.9391004010.34030.0164500.016510.28269.634100由前三阶振型图可以看出，该结构的动力特征具有以下特点:1).第一振型以纵向平动为主2).第三振型表现为横向平动3).第二振型是扭转的振动，这主要是因为该结构纵横方向刚度不同引起的。由前三个振型图可以看出，两个水平方向的振型参与系数不在同一个数量级，即不存在明显的扭转效应，可见本结构具有良好的抗扭刚度。二结构的在地震作用下的响应分析 工程结构的抗震设计，除了按建筑物重要性分类确定设防烈度、选择有利于建筑物的地段和场地类型、确定有利抗震的结构方案外，还要由分析计算的方法，求出预期地震作用下结构所产生的内力和变形，这就是结构的地震反应分析。目前，用于结构地震反应分析的方法，大多属惯性力作用一下的振动解法。高层建筑结构采用的方法主要有振型分解反应谱法和时程分析法。2.1谱分析基本理论谱分析是一种将模态分析的结果与一个己知的谱联系起来计算模型的位移和应力的分析技术。谱分析替代时间一历程分析，主要用于确定结构对随机载荷或随时间变化载荷的动力响应情况振型分解反应谱理论的基本假设：（1）是结构的地震反应是线弹性的，可采用叠加原理进行组合; （2）结构的基础是刚性的，所有支承处地震动完全相同，在MIDAS软件中即为单点响应谱分析; （3）结构最不利地震反应为其最大地震反; （4）地震动过程是平稳随机过程。 本结构为钢筋混凝土结构，其阻尼比取0.05。地震影响系数的阻尼调整系数取1.0，本工程为7度(0.10g)，设计地震分组为一组，类场地，于是得到设计反应谱如图5所表示.图5加速度反应谱2.2反应谱分析结果 本工程采用单点响应谱进行了结构的抗震分析，即对基底各节点分别沿X,Y方向输入地震加速度反应谱曲线来计算结构的地震响应。本工程考虑各振型间的祸合作用，所以选用CQC法。对结构采用振型分解反应谱法计算地震作用下的结构位移和内力。结构的底部为固结，下面为MIDAS/GEN计算的结果。2.3位移分析: 由结构的变形图可看出结构的水平侧移随楼层的变化比较均匀，也说明结构刚度均匀。具体图形如下面图形所示。(1) 结构的整体变形和位移图6 X一反应谱作用下整体位移图7 Y一反应谱作用下整体位移(2) 结构的位移沿层高变化表3 X一反应谱作用下楼层位移层标高(m)层高度(m)最大位移(m)平均位移(m)屋顶16.500.01250.01256F13.23.30.01150.01155F9.93.30.00990.00994F6.63.30.00790.00793F3.33.30.00550.00552F03.30.00290.00291F-4.54.500表4 Y一反应谱作用下楼层位移层标高(m)层高度(m)最大位移(m)平均位移(m)屋顶16.500.00020.00016F13.23.30.00020.00015F9.93.30.00020.00014F6.63.30.00020.00013F3.33.30.00010.00012F03.30.000101F-4.54.500 由结构在X, Y反应谱作用下的整体位移沿层高变化图形可知，X方向屋顶最大位移为12.5毫米，Y方向得到的位移很小，可以忽略不记。2.4层间位移角分析图5 X向反应潜层一层间位移角层层高度(m)层间位移角限值全部竖向单元的最大层间位移质心计算的层间位移节点层间位移(m)层间位移角验算层间位移(m)层间位移角(最大/当前方法)层间位移角验算6F3.31/5503010.00121/2753OK0.001211/2753OK5F3.31/5502410.00171/1894OK0.001711/1894OK4F3.31/5502110.00211/1541OK0.002111/1541OK3F3.31/5501210.00251/1344OK0.002511/1344OK2F3.31/55010.00271/1232OK0.002711/1232OK1F4.51/550610.00291/1571OK0.002911/1571OK图6 Y向反应潜层一层间位移角层层高度(m)层间位移角限值全部竖向单元的最大层间位移质心计算的层间位移节点层间位移(m)层间位移角验算层间位移(m)层间位移角(最大/当前方法)层间位移角验算6F3.31/55030101/123492OK000OK5F3.31/55024101/89932OK000OK4F3.31/55018101/74693OK000OK3F3.31/55012101/66226OK000OK2F3.31/55010.00011/62998OK000OK1F4.51/550610.00011/89786OK000OK 反应谱分别从X, Y向输入，得到各点的位移，经过整理，得到各个各层位移的表格。X方向最大层间位移角为6F的1/2753,Y方向很小的，可以忽略不记。小于规范规定的1 /550，完全满足抗震变形验算的要求。2.5最大剪力分析 从图中可以看出在X向反应谱作用下整个楼层的剪力最大值为200.7KN，在Y向反应谱作用下整个楼层的剪力最大值为167.9KN，发生的位置在楼层的最底层1F。通过框架和抗震墙之间的协同工作分析，进行地震内力分配。 图7 X向反应谱层一层剪力图8 Y向反应谱层一层剪力2.6弯矩分析各个楼层在反应谱X方向、Y方向作用下楼层的倾覆弯矩如下。荷载工况层标高(m)层高度(m)角度1(deg)竖向构件的倾覆弯矩(kN*m)角度2(deg)竖向构件的倾覆弯矩(kN*m)框架墙单元框架墙单元Value比值Value比值Value比值Value比值静力荷载工况结果角度: 0 度输入角度后请按 适用 键。0rx(RS)6F13.23.304608.62100900000rx(RS)5F9.93.3011539.94100900000rx(RS)4F6.63.3019738.5100900000rx(RS)3F3.33.3028692.85100900000rx(RS)2F03.3038461.7100900000rx(RS)1F-4.54.5053124.2100900000表7 X向反应谱层一倾覆弯矩表8 Y向反应谱层一倾覆弯荷载工况层标高(m)层高度(m)角度1(deg)竖向构件的倾覆弯矩(kN*m)角度2(deg)竖向构件的倾覆弯矩(kN*m)框架墙单元框架墙单元Value比值Value比值Value比值Value比值静力荷载工况结果角度: 0 度输入角度后请按 适用 键。0ry(RS)6F13.23.39001001800000ry(RS)5F9.93.39001001800000ry(RS)4F6.63.39001001800000ry(RS)3F3.33.39001001800000ry(RS)2F03.39001001800000ry(RS)1F-4.54.590010018000002.7反应谱方法的优缺点反应谱方法采用动力方法计算地震反应，考虑了地面运动的强弱、场地土性质及结构动力特性对地震惯性力的影响，能够在一定程度上代表地震对房屋的作用。反应谱法用于设计比较方便，求出等效地震荷载后按照静力方法进行计算，所得到的内力能够代表在地震作用下的不利内力，可以根据它们设计截面。 但是，反应谱方法也存在一些缺陷:(1)反应谱主要依据的是加速度反应谱，等效地震荷载的大小与地面加速度的峰值有密切关系，但未反映地面运动中速度、位移及持续时间等参数的影响。研究表明，地面运动中速度对结构反应影响很大，在相同加速度峰值下，速度值越大，结构反应也越强烈。持续时间对地震反应也有很大影响，特别是当结构进入弹塑性状态后，长持时地震的危害更大，但反应谱方法中未能考虑这些不利因素。(2)反应谱计算是建立在弹性动力分析基础上，它未考虑结构弹塑性性能的影响。我国抗震规范规定，在小震下用等效静力地震荷载计算弹性状态下的内力和位移，是比较合理的。如果用反应谱方法计算弹塑性状态下的结构往往得不到合理的结果，而通常在设防烈度地震作用下，结构大多是进入弹塑性状态的。(3)高层建筑都是多质点体系，而反应谱曲线是从单质点体系得到的。虽然用