关于ANSYS地震响应分析的一些讨论.docVIP

地震响应分析 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）求出结构的自振一阶频率w ???2）令tr＝1/w ???3) 定义ramp荷载为从0到tr加到9.81，然后在整个时间积分中保持不变 ???4）antype,trans中分几个荷载步将荷载从0加到9.81 ???5) 在随后的荷载步中acel，，9.81+ac（i） ???这种做法虽然也是将9.81++加到地震波中，但是因为满足TR的要求，所以这个动力效应被削弱到了静力效应，它作用在结构上就像静载一样。对于单自由度结构理论上跟静载是完全一样的，但是多自由度会子静力效应上下很小的范围内波动，所以可以认为相当于静载的作用，这样我们就可以达到考虑恒载的目的了。 ????第5种是我至今为止考虑恒载的做法，我也很想知道还有没有更简单精确的方法，或者在前4种方法中就有只是我使用不正确，希望大家能一起来讨论，彻底解决这个问题。谢谢！ SOLU ANTYPE, TRANS TRNOPT,FULL TIMINT,OFF !*先关闭时间积分效应 TIME,1E-8 !*设一个极短的积分时间 acel,,9.8 NSUBST,2? ?!有时候子步数要增大 KBC,1 LSWR,1 !*把这个写入第一步 TIMINT,ON !*然后再时间积分效应开关，以后就正常写载荷步了 ??这种方法应该是对的,ANSYS帮助文件中也有提到, 可是,有一个问题:由于是阶跃荷载,就会产生动力效应,整个结构的变形大于实际的情况吧?这样与实际结构在重力下受到的变形就不一样了! ??我想,之所以用KBC,1的原因是产生阶跃荷载,第一个子步产生变形,然后后几个子步使其稳定,不产生初始速度吧 ??自己作了一个例子, 终于彻底搞明白啦, KBC,1是产生阶跃荷载,不过,关键就在于此时关掉了时间积分效应,TIMINT,OFF??那么,此时就相当于是静力分析了,当然就不会产生动力效应了,整个结构的变形就是实际变形.??另外KBC,1的作用是不产生初始速度. ??总结下:非线性时程分析时,结构在地震波来之前,已经存在初始位移和初始内力,但是没有初始速度和加速度.因此,一定要注意KBC,1,同时TIMINT,OFF ,设置2个以上的子步,很小的时间. ??一楼的总结有很多不妥的地方, 你那个方法不知道正确不,而且太麻烦了. tmh_cn在第12楼总结的方法是对的,其他人说的方法都有问题. 另外,不用打开应力刚化 .fini /config,nres,20000 !/config,nproc,2 lstep=1500 !el *DIM,aei,ARRAY,lstep,1,1, , , *CREATE,ansuitmp *VREAD,aei(1),el,txt,,1, (e16.6) *END /INPUT,ansuitmp csys,0 /solu antype,trans trnopt,full TIMINT,OFF !*先关闭时间积分效应 SSTIF,On KBC,1 TIME,1e-3!*设一个极短的积分时间 acel,,8.7 NSUBST,2!有时候子步数要增大 solve timint,on kbc,0 !TINTP, ,0.25,0.5,-1,0.5,-1, ki=0.05 kj=0.05 oi=2*3.14159*0.09 oj=