资源描述：
《ansys一个火车过桥动态加载的例子by paradiseboy》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、ANSYS一个火车过桥动态加载的例子byparadiseboy[分享][转帖]ANSYS一个火车过桥动态加载的例子!用对话框提示用户输入基本数据!:reread!multipro,'start',3!*cset,1,3,y_bot,'enterthevalueofy_bot',0!*cset,4,6,y_top,'enterthevalueofy_bot',0!*cset,7,9,ttype,'enterthevalueoftraintype',0!输入火车活载的类型!multipro,'end'y_bot=13.9y_top=14.1ttype=0*if,ttype,eq,0,
2、thenq1=147000!中-活载q2=92000q3=80000L1=7.5L2=30*elseif,ttype,eq,1,thenq1=0q2=0q3=0L1=0L2=0*else*msg,uierrorvalueofttype(1or2),pleaseenteragain*go,:reread*endif!识别桥面过列车的单元，并将单元按节点坐标从小到大排序lsel,s,loc,y,y_bot,y_topesll,s!选中桥面单元elmax=elmiqr(0,14)!最大单元号ndmax=ndinqr(0,14)!最大节点号*dim,eldk,,elmax!存放桥面过车单
3、元*dim,ndx,,elmax!存放过车单元的节点x坐标（两节点坐标中的较小值）*dim,nodes,,ndmax!存放桥面过车的节点j=0*do,i,1,elmax*if,esel(i),gt,0,thenj=j+1eldk(j)=i*endif*enddo!记录桥面单元的标号，保存于数组eldk中eldknum=j!总桥面单元数!对eldk数组排序，按节点坐标从小到大排列*do,i,1,eldknumii=eldk(i)ndi=nelem(ii,1)ndj=nelem(ii,2)i1=nx(ndi)i2=nx(ndj)ndx(i)=i1*if,i2,lt,i1,thennd
4、x(i)=i2*endif*enddo*do,i,1,eldknum-1*do,j,i+1,eldknum*if,ndx(j),lt,ndx(i),thentemp=ndx(i)ndx(i)=ndx(j)ndx(j)=temptemp=eldk(i)eldk(i)=eldk(j)eldk(j)=temp*endif*enddo*enddonsle,s!选中桥面单元的节点j=0*do,i,1,ndmax*if,nsel(i),gt,0,thenj=j+1nodes(j)=i*endif*enddondnum=j!总桥面节点数xmin=nx(nelem(eldk(1),1))temp
5、=nx(nelem(eldk(1),2))*if,temp,lt,xmin,thenxmin=temp*endifxmax=nx(nelem(eldk(eldknum),1))temp=nx(nelem(eldk(eldknum),2))*if,temp,gt,xmax,thenxmax=temp*endifqth=xmax-xmin!桥面过车长度lth1=2*qth+L1+L2+182.5lth=(nint(lth1/5)+8)*5!总长度dim=lth/5+1*dim,disy,table,dim,ndnum!桥面节点的y方向的位移值*dim,disy1,table,dim,
6、ndnum!桥面节点的y方向的位移值(从大到小)*dim,disx,table,dim,ndnum!桥面节点的x方向的位移值*dim,disa,table,dim,ndnum!桥面节点的转角*dim,disymax,,ndnum*dim,diy,table,dim*dim,ffx,table,elmax,dim!单元轴力*dim,fmz,table,elmax,dim!单元弯矩*dim,detaf,,elmax!每个单元轴力变化值*dim,detaf1,,elmax!每个单元轴力变化值（从小到大）*dim,elnum,,elmax!轴力变化最大的单元号*dim,detaf2,,2
7、0!前20个轴力变化最大值*dim,elnum2,,20!前20个轴力变化最大值的单元*dim,fmax,,elmax!每个单元的轴力最大值*dim,fmin,,elmax!每个单元的轴力最小值*dim,fabs,,elmax!每个单元轴力最大值（绝对值）*dim,elnum3,,elmax!轴力最大的单元号*dim,elnum4,,20!前20个轴力最大的单元号*dim,nmax,,20*dim,dfmaxel,table,dim,20!前20各轴力变化最大单元的轴力*dim,
