结构动力学.docx

《结构动力学.docx》由会员分享，可在线阅读，更多相关《结构动力学.docx（11页珍藏版）》请在第一文库网上搜索。

1、斜拉桥地震响应分析摘要：斜拉桥在地震波荷载作用下有极其复杂的振动响应，本文采用ANSYS有限元软件对某斜拉桥在centro波作用下动力响应进行了分析。得出结论：ANSYS有限元软件能为复杂大跨度结构的抗震性能分析提供高效、可靠的计算平台；对于复杂结构或异性结构,谱分析的结果未必偏于安全,这时采用地震波瞬态分析更精确。因此，应用ANSYS有限元软件分析斜拉桥的动力响应有较好的效果，并且centro波可以作为结构动荷载的近似标准波使用。关键词：斜拉桥；动力分析；CentrO波；ANSYS有限元一、概述对于桥梁而言，地震所带来的破坏，无论从数量上，还是从程度上，都大大超过其他自然灾害的破坏。严重的桥

2、梁灾害不仅直接影响交通，而且经常引发次生灾害，从而加剧地震灾害的严重性。为了减轻地震所造成的损失，既要对桥梁做好抗震加固工作，更需在桥梁设计上采取措施以满足抗震要求。因此，对桥梁的地震响应进行相应的分析是有必要的。1 .地震作用理论(1)直接动力分析理论1900年，日本大森房吉教授提出了静力理论。静力理论不考虑建筑物的动力特性。假设结构物为绝对刚性，地震时建筑物的运动与地面运动完全一致，建筑物的最大加速度等于地面运动的最大加速度。建筑物所受的最大地震荷载F等于其质量m与地面最大加速度anm的乘积，即：由于这种方法比较简单，且用这种方法设计的建筑物大多经受了一般地震的考验，所以，它稍作修改后至今

3、仍被某些国家的抗震设计规范所采用。但是这种方法完全忽略了结构本身动力特性的影响。因为只有当结构的基本固有周期比地面运动周期小的多时，结构在地震时才有可能不产生变形而被视为刚体。所以静力理论只适合于低矮的、刚性较大的建筑，如路基、挡土墙和重力式桥台等。(2)直接动力分析理论随着动力学的发展，美国学者在20世纪40年代提出了计算地震力的反应谱理论，也称动力法。它考虑了地震时地面的运动特性与结构物自身的动力特性，是当前工程设计应用最为广泛的抗震设计方法。反应谱理论是以单质点体系在实际地震作用下的反应为基础来分析结构反应的方法。质量为m的单质点体系在水平方向地面运动分量作用下质点绝对加速度为a(t),

4、则质点所受的水平地震作用为：F(t)=ma(t)上式表示在地震过程中，质点水平地震作用的大小与方向随时间t变化。抗震设计中通常只需地震作用的最大值，其值可表示为：F=mama,=(W/g)manax=W式中，W=mg为质点重量，g为重力加速度，=anJg称为水平地震影响系数。我国学者根据国内外数百条地震记录的反应谱进行统计分析后，建立了地震响应系数与结构体系自振周期T的关系曲线(T),相关参数可查相应规范。一般认为竖向地震影响系数4(丁)的曲线形状与水平地震响应系数大体相同，可直接使用(T)曲线，数值上一般取%=(1223),竖向地震作用可由下式计算。F=a1W(3)直接动力分析理论尽管反应谱

5、分析方法在结构抗震设计中得到广泛应用，将实际地震的加速度时程记录输入结构计算模型，直接分析结构的地震反应已成为可能。可直接获得地震过程中结构节点各时刻的位移、速度和加速度，从而计算各时刻竖向地震作用和构件的地震内力。这种方法也称时程分析法。2 .ANSYS分析方法(1)谱分析技术谱分析是一种将模态分析结果与一个已知的谱联系起来，然后计算模型的位移和应力的分析技术。主要用于确定结构在随机载荷或随时间变化载荷作用下的动力响应，如地震、风载等。谱分析必须要已知结构的振型和固有频率，因此需先进行模态分析。在扩展模态时，只需扩展到对最后进行谱分析有影响的模态即可。(2)瞬态动力学分析瞬态动力学分析，也称

6、时间历程分析，可以用来分析结构在承受任意的、随时间变化的载荷作用时的动力响应。瞬态分析必须指定系统的杨氏模量，密度等。二、算例分析1.地震荷载(1)地震谱由于为近震，场地为In类，7；=0.4,tzm=0.08,所以根据前面反应谱理论可得以下的地震响应系数方程：=0.44T0.036,0T0.1a=0.08,0.1T0.4a=(0.4/T)0-9*0.08,0.4T3由上式可得地震谱如表1所示。表1地震频谱值自振周期频率谱值自振周期频率谱值0.02500.04480.33.3330.080.04250.05360.362.7780.080.0616.6670.06240.42.50.080.0

7、812.50.07120.81.250.0428710.1100.081.20.8330.0297630.128.3330.081.60.6250.0229740.166.250.0820.50.0187940.250.082.40.4160.015950.244.1670.082.80.3570.0138840.283.5710.0830.3330.013(2)地震波地震波采用Centro波输入，考虑纵波和横波的双向作用情况。限于篇幅，此处不再列举地震波的加速度参数。2 .斜拉桥技术参数(1)材料性能主梁、索塔:E=35e10,p=2500/,=0.17刚性鱼刺横梁和主塔连接横梁：E=1.0

8、d6,p=0,=0斜拉索：E=1.611,p-1600,=0.3(2)截面特性主梁:。=16,b=1.6,A=25.6,v=2000/=20/=21.8上索塔:。=4.7,b=3.4,A=16,4=39.7/=77.7,4=61.6中索塔:=9,b=6,A=54,1,=200/=450/=61.6下索塔:。=8,力=5,A=40,=83.3/=213.34=333.3刚性鱼刺横梁和索塔横梁:a=1*=1,A=I,/y=1/12,/二=1/12,4=1/3主梁:A=0.012(3)结构尺寸本例双塔斜拉桥，主跨360m,边跨174m,塔高162m。桥宽28m,主塔连接横梁长30m。塔的倒Y分叉点距

9、桥面60m,塔底距桥面30m,塔底横桥向距离20mo斜索在主梁每6m处布置一根，索塔处也布置横梁且被斜拉索吊住。主塔从塔顶往下每隔18m设置一个斜索张拉集中点，共分4个张拉点。上面3个集中点每个单侧可以张拉7条索，第4个点张拉8条斜拉索，总计116条索。(4)边界条件左桥端仅给予竖向和横向的平移自由度，右桥端仅给予横向的平移自由度约束，索塔底部完全约束，刚横梁在索塔处仅给予横向和竖向约束，索单元和梁单元给予完全较约束。3 .斜拉桥ANSYS模型的建立采用空间鱼刺模型。所有截面几何特性以及质量集中于主梁上，横梁相对为刚性构件，只起到传递力的作用。主梁处于全漂浮状态，主梁在索塔处以及两个桥端都释放

10、纵向约束。中、下塔柱使用一个梁单元模拟，上塔柱使用四个梁单元模拟，索塔横梁采用两个刚性单元模拟。地震载荷取横桥向和纵桥向双向作用的情况，采用瞬态动力分析的方法求解，不考虑阻尼地震波，采用桥梁地震反应分析数据计算。瞬态动力分析需要先分步写激励数据文件步，如果有N个载荷子步，那么就需要写N个载荷文件。这里每个子步时间间隔为0.1秒，沿桥纵向共116个单元,每个单元117个节点，每个单元长6m,总计将计算117步。由于该斜拉桥结构较为复杂，采用APD1语言创建节点和单元，直接建立模型，图1即为所得的斜拉桥模型。ANSYSE-NJU1102011U14:55:09ROT图1.斜拉桥模型4 .地震响应分

11、析（1）地震谱响应分析启用1ANB求解器，提取模态为20阶，扩展前10阶。频率起始值和终止值分别为0和100。模态计算结果得到前10阶频率如表2所示。这里提取适当的几阶振型以说明问题，第一阶为主跨横向弯曲，第二阶为全桥纵向平移，以下几阶基本是桥面在水平面内的弯曲，到了第六阶就开始变成竖向挠曲振动了。这里列出了1、2、5、6阶振型图，分别如图2图9所示。一般来说纵飘发生在第一阶频率处，然而这里主跨横弯变成了第一阶，所以横向弯曲刚度是不够的。总的来说纵飘处于横向和竖向运动的前面，横向运动处于竖向运动的前面，竖向挠曲刚度是三者中最小的。表2.前十阶频率SETTIME/FREQ1OADSTEPSUBS

12、TEP10.89634E-011120.270671230.272911340.317111450.501081560.604121670.657001780.659831890.7120119100.81878110NSYSJU111201108:26:441NODA1SO1UTIONSTEP=ISUB-1FREQ=.089642UY(AVG)RSYS=ODMX=.309E-03-.309E-03-.172E-03-.343E-04.103E-03-.240E-03-.103E-03.343E-04.240E-03.172E-03.309E-03图2.第一阶频率振型图（俯视）STEP-ISU

13、B=1FREQ=.089642UY(AVG)RSYS=ODMX=.309E-03SMN=-.104E-05SMX=.309E-03JU111201108:30:30-.104E-05.679E-04.137E-03.206E-03.275E-03.334E-04.102E-03.171E-03.240E-03.309E-03ANSYSJU111201108:31:53图3.第一阶频率振型图(等侧)NODA1SO1UTIONSTEP=ISUB2FREQ=.270711UY(AVG)RSYS-ODMX=.424E-03SMKSMX.330E-03.235E-03.424E-03-.424E-03-

14、.235E-03-.471E-04.141E-03-.330E-03-.141E-03.471E-04图4.第二阶频率振型图（侧视）JU111201108:35:41FREQ=.270711NSYSJU111201108:47:50-.240E-03-.114E-03.124E-04.139E-03-.177E-03-.508E-04.757E-04.265E-03.202E-03.328E-03图6.第五阶频率振型图（俯视）STEP1SUB=2-.827E-O4.299E-04.142E-03.255E-03.367E-03-.264E-O4.861E-04.199E-03.311E-03.424E-03图5.第二阶频率振型图(等侧)1NODA1SO1UTIONSTEP=ISUB=5FREQ=.501139UY(AVG)RSYS=ODMX=.329E-03SMN=-.240E-03SMX-.328E-03JU111201108:37:08STEP=ISUB=5FREQ=.501139UY(AVG)-.328E-03-.182E-03-.365E-04.109E-03-.255