火车过桥动态加载实例(部分)

发布日期:[08-12-04 14:51:08] 浏览人次:[]

变化最大的单元
elnum2 = elnum
*enddo
*do,i,1,20 ! 前20个内力变化最大的单元
numel = elnum2
dfmaxel(0,i) = i
dmmaxel(0,i) = i
*do,j,1,dim
kk = (j-1)*5
dfmaxel(j,0) = kk
dfmaxel(j,i) = ffx(numel,kk) ! 轴力变化最大单元的轴力
dmmaxel(j,0) = kk
dmmaxel(j,i) = fmz(numel,kk) ! 相应的弯矩
*enddo
*enddo
! 前20个内力最大的单元
*do,i,1,elmax
fabs =abs(fmax )
*if,fabs ,lt,abs(fmin ),then
fabs =abs(fmin )
*endif
*enddo
*do,i,1,elmax
elnum3 =i
*enddo
*do,i,1,elmax-1
*do,j,i+1,elmax
*if,fabs(j),gt,fabs ,then
temp = fabs
fabs = fabs(j)
fabs(j) =temp
temp = elnum3
elnum3 = elnum3(j)
elnum3(j) = temp
*endif
*enddo
*enddo
*do,i,1,20
elnum4 = elnum3
nmax = fabs
*enddo
*do,i,1,20 ! 前20个内力最大的单元
numel = elnum4
fmaxel(0,i) = i
mmaxel(0,i) = i
*do,j,1,dim
kk = (j-1)*5
fmaxel(j,0) = kk
fmaxel(j,i) = ffx(numel,kk) ! 轴力最大单元的轴力
mmaxel(j,0) = kk
mmaxel(j,i) = fmz(numel,kk) ! 相应的弯矩
*enddo
*enddo
save,,,,all
!/title,the 91 node displament vs train position ! 显示图形的名称
!*vplot,diy(1,0),diy(1,1)
!/ui,copy,save,bmp,graph
!/title
!/replot
*do,i,1,30
cc = nodes
/title, y displacement of node %cc% VS. train position
*vplot,disy1(1,0),disy1(1,i)
/ui,copy,save,bmp,graph
*enddo
/title
/sys,rename njq*.bmp aaa*.bmp
*do,i,31,ndnum
cc = nodes
/title, y displacement of node %cc% VS. train position
*vplot,disy1(1,0),disy1(1,i)
/ui,copy,save,bmp,graph
*enddo
/title
/sys,rename njq*.bmp aa*.bmp
!*do,i,1,20
!cc = elnum2
!/title,axis of element %cc% VS. train position
! *vplot,dfmaxel(1,0),dfmaxel(1,i)
!/ui,copy,save,bmp,graph
!/title,moment of element %cc% VS. train position
! *vplot,dmmaxel(1,0),dmmaxel(1,i)
!/ui,copy,save,bmp,graph
!*enddo
!/title
*do,i,1,20
cc = elnum4
/title,axis of element %cc% vs train position
*vplot,fmaxel(1,0),fmaxel(1,i)
/ui,copy,save,bmp,graph
/title,moment of element %cc% vs train position
*vplot,mmaxel(1,0),mmaxel(1,i)
/ui,copy,save,bmp,graph
*enddo
/title
*cfopen,result,txt,,,
*vwrite
('' '')
*vwrite
(''**********************南京桥计算结果*****************************'')
*vwrite
('' '')
*vwrite
('' ****节点最大位移**** '')
ndnum1 = 5*nint(ndnum/5)
*do,i,1,ndnum1,5
a=nodes
b=nodes(i+1)
c=nodes(i+2)
d=nodes(i+3)
e=nodes(i+4)
*vwrite,a,b,c,d,e
(''荷载位置'',5(f8.0,5x))
*vwrite,disy1(1,0),disy1(1,i),disy1(1,i+1),disy1(1,i+2),disy1(1,i+3),disy1(
1,i+4)
(f5.0,,5(3x,f10.5))
*enddo
*vwrite
('' ****前20个轴力变化最大的单元和相应的轴力值**** '')
*vwrite,elnum2(1),detaf2(1)
(f15.0,'' '', f15.1)
*vwrite
('' '')
*do,i,1,18,3
*vwrite
('' '')
a= elnum2
b= elnum2(i+1)
c= elnum2(i+2)
*vwrite,a,b,c
(''荷载位置'',9x,''单元'',f4.0,14x,''单元'',f4.0,14x,''单元'',f4.0)
*vwrite,dfmaxel(1,0),dfmaxel(1,i),dmmaxel(1,i),dfmaxel(1,i+1),dmmaxel(1,i+1)
,dfmaxel(1,i+2),dmmaxel(1,i+2)
(f5.0,3(f11.1,'' '',f11
|<< << < 1 2 3 4 5 6 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:火车  过桥  动态  加载  实例  部分
相关资讯
热点资讯
推荐资讯

关于我们 | 站点导航 | 使用帮助 | 友情链接 | 广告服务 | 免责声明 | 新手上路
设为首页 | 加入收藏 | 在线留言 | 马棚网QQ群:{92562572}{102901272}{333259257} | 交流QQ: 客户服务 客户服务 客户服务