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

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

1*c2+q2*c3)/(elth2)
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
c4 = xlth3 - tpos3
c5 = elth3 - c4
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q2*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*elseif,tpos1,le,(qth+L1+L2),then !煤水车头出桥
c4 = xlth3 - tpos3
c5 = elth3 - c4
dtele1 = eldknum - tele3
*if,dtele1,eq,0,then
qq = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq,qq,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*else
*do,j,tele3+1,eldknum
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq,qq,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*elseif,tpos1,le,(qth+L1+L2+182.5),then ! 全桥均布列车荷载(182.5为任意取的
常数)
*do,j,1,eldknum
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*elseif,tpos1,le,lth1,then ! 均布列车荷载出桥
c1 = xlth4 - tpos4
c2 = elth4 - c1
*do,j,tele4+1,eldknum
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
qq = q3*c1/elth4
sfbeam,eldk(tele4),1,pres,qq,qq,,,,,
*endif
/psf,pres,norm,1
/replot
solve
finish
/post1
diy(time,0) = tpos1
diy(time) = uy(91) ! 得到中间节点的y坐标
*vplot,diy(1,0),diy(1,1) ! 动态显示中间节点位移结果
disx(time,0) = tpos1
disy(time,0) = tpos1
disa(time,0) = tpos1
*do,j,1,ndnum
disx(0,j) = nodes(j)
disy(0,j) = nodes(j)
disa(0,j) = nodes(j)
disx(time,j) = ux(nodes(j))
disy(time,j) = uy(nodes(j))
disa(time,j) = rotz(nodes(j))
*enddo
etable,forcex,smisc,1 ! 定义单元轴力表
etable,momentz,smisc,6 ! 定义单元弯矩表
*vget,ffx(1,time),elem,1,etab,forcex,,,,
*vget,fmz(1,time),elem,1,etab,momentz,,,,
ffx(0,time) = tpos1
fmz(0,time) = tpos1
*if,time,eq,1,then
*do,j,1,elmax
ffx(j,0) = j
fmz(j,0) = j
*enddo
*endif
*enddo
save,,,,all
! 结果处理
*do,i,1,ndnum
nd = nodes
disymax = abs(disy(1,nd))
*do,j,2,dim
*if,abs(disy(j,nd)),gt,disymax ,then
disymax = abs(disy(j,nd))
*endif
*enddo
*enddo
*do,i,1,ndnum-1
*do,j,i+1,ndnum
*if,disymax(j),gt,disymax ,then
temp = disymax
disymax = disymax(j)
disymax(j) = temp
temp = nodes
nodes = nodes(j)
nodes(j) = temp
*endif
*enddo
*enddo
*do,i,1,ndnum ! 前5个位移最大的节点
numnd = nodes
disy1(0,i) = numnd
*do,j,1,dim
kk = (j-1)*5
disy1(j,0) = kk
disy1(j,i) = disy(kk,numnd)
*enddo
*enddo
*do,i,1,elmax
fmin = ffx(i,1)
fmax = ffx(i,1)
*do,j,2,dim
*if,ffx(i,j),gt,fmax ,then
fmax = ffx(i,j)
*endif
*if,ffx(i,j),lt,fmin ,then
fmin = ffx(i,j)
*endif
*enddo
detaf = fmax - fmin ! 每个单元的内力变化值
detaf1 = detaf
*enddo
*do,i,1,elmax
elnum = i
*enddo
*do,i,1,elmax-1
*do,j,i+1,elmax
*if,detaf1(j),gt,detaf1 ,then
temp = detaf1
detaf1 = detaf1(j)
detaf1(j) = temp
temp = elnum
elnum = elnum(j)
elnum(j)=temp
*endif
*enddo
*enddo
*do,i,1,20
detaf2 = detaf1 !前20个内力
|<< << < 1 2 3 4 5 6 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:火车  过桥  动态  加载  实例  部分
相关资讯
热点资讯
推荐资讯

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