资源描述:
《改的时程分析.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、finish/config,nres,/config,fsplit,7500!!!文件分割/prep7*create,ed*dim,wind1,table,4000,19!定义数组,数组名,数组类型,数组下标最大值(改成你自己的点数)*tread,wind1,td,txt,,0!wind为参数名称li带路径的文件名称(改成你自己的txt文件名)*end/input,ed*dim,npa,array,19!(改成你自己的点数)npa(1)=6,11,16,21,26,31,36,41,46,51npa(11)=56,61,66,71,76,81,86,91
2、,96!(这里是19个点的点号,改成你自己的点号和点数)!以上为输入数组,下面为求解finish/prep7nsel,s,node,,1,100,1,0cm,jiedian,node!输出节点编组esel,s,elem,,1,100,1,0cm,danyuan,elem!输出单元编组allsel,allfinish/config,nproc,2!双核运算/soluantype,transallsel,alltrnopt,fullupcoord,1,on!对变形后的结构重新计算TIMINT,on!使用时间积分nlgeom,onsstif,onautots,
3、onsolcontrol,onoutres,all,none!不输出结果信息outres,nsol,5,jiedian!输出节点自由度解!outres,strs,5,danyuan!输出单元节点应力outres,esol,5,danyuan!输出单元节点应力,每2个荷载步输出一次解outres,nload,5,danyuanKBC,0!斜坡荷载DELTIM,0.1!dmprat,0.02!常数阻尼0.02,跟瑞雷阻尼差别不大alpha,0.0015beta,0.0017*do,t,1,4000,1time,t*0.1nsubst,1!指定1个荷载子步kb
4、c,0*do,jj,1,16,1!52层节点循环!(改成你自己的点数)f,npa(jj),fxwind1(t,jj)!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!fx方向(改成你自己的风作用方向)*enddoacel,,,9800!施加自重,这个自重容易造成塔线出现smallequation,进而不收敛(这个是在Z方向施加的重力)allsel,allsolve*enddofinishsave/post26ESOL,3,879,1150,smisc,1!!!!!!!!!!!!!变量号单元号节点号轴力/axlab,x,time/axlab,
5、y,valuxvar,1plvar,3/post26nsol,2,51,u,y,deflection!!!变量号单元号节点位移/axlab,y,valuxvar,1plvar,2/post1etable,fxi,smisc,1/post26numvar,50!!!!!指定变量的数目nsol,3,51,u,ylines,10000!!!!!!!指定每一页出现的数量prvar,3!列表显示变量