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

地震响应分析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）令tr1/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,2lstep=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,transtrnopt,fullTIMINT,OFF !*先关闭时间积分效应 SSTIF,OnKBC,1TIME,1e-3!*设一个极短的积分时间 acel,8.7 NSUBST,2!有时候子步数要增大 solvetimint,onkbc,0 !TINTP, ,0.25,0.5,-1,0.5,-1,ki=0.05kj=0.05oi=2*3.14159*0.09oj=2*3.14159*0.16ma=2*(ki*oj-kj*oi)*oi*oj/(oj*oj-oi*oi)mb=2*(kj*oj-ki*oi)/(oj*oj-oi*oi)ALPHAD,ma,BETAD,mb, allsel,all*do,tm,1,1500,1time,(tm)*2/100NSUBST,1, , ,1 AUTOTS,1acel,aei(tm)/314.7 !EI波，只输入X向加速度solve*enddofini*关于动力学分析模拟自重对结构的影响，这个问题，现在偶已经解决这个问题啦，几行命令流就可以了：/PREP7ACEL,0,9800,0, !*先把重力加速度在预处理设好，以后就不在设了FINISH/SOLUANTYPE, TRANSTRNOPT,FULLTIMINT,OFF !*先关闭时间积分效应TIME,1E-8 !*设一个极短的积分时间NSUBST,2KBC,1LSWR,1 !*把这个写入第一步!*TIMINT,ON !*然后再时间积分效应开关，以后就正常写载荷步了.虽然，这只是个很简单的问题，但如果不注意，会使计算的动力响应值受重力影响而偏大，我当时就很郁闷，怎么也不能把重力的影响去掉，现在终于搞定了！现与大家分享偶的一点心得，同时，也感谢各位大侠的关注：-） 模态叠加法（时间可能取的太短了）求地震反应! 步骤1：进行模态分析/SOLU! 1:MODE ANALYSISAnType,ModalModOpt,Lanb,20! EXTRACT 20 MODESMxpand,20,1SOLVEFINISH! 时程分析命令流! 采用模态叠加法*set,st,0.01*set,endpoint,100*dim,xforce,array,endpoint*vread,xforce(1),Seismic_x,txt(7F10.5) *VOPER,xforce(1),xforce(1),MULT,9.8, , ,/SOLU !* ANTYPE,4TRNOPT,MSUP,20,1 ! ! ! 施加位移约束NSEL,S,LOC,Z,0,0! 选取底层所有节点D,ALL, , , , , ,ALL, , , , ,! 约束所有自由度 ALLSEL,ALL!* *dim,shijian,array,endpoint*vfill,shijian(1),ramp,st,stNt_star=1Nt_end=100Nt_incr=1*DO,Nt,Nt_star,Nt_end,Nt_incrTIME,shijian(Nt)ACEL,xforce(Nt),NSEL,ALL!LSWRITE,NtSOLVE*ENDDOFINISH!模态扩展/SOLU EXPASS,ONNUMEXP,all,0,1OUTRES,all,allSOLVEFINISH*谱分析我想，作谱分析之前应该先弄明白谱分析的工程意义。实际上，谱分析是一种线性化的动力分析思路，简单点说，就是在获得