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

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

! 用对话框提示用户输入基本数据
!:reread
!multipro,''start'',3
! *cset,1,3,y_bot,''enter the value of y_bot'',0
! *cset,4,6,y_top,''enter the value of y_bot'',0
! *cset,7,9,ttype,''enter the value of train type'',0 !输入火车活载的类型

!multipro,''end''
y_bot=13.9
y_top=14.1
ttype = 0
*if,ttype,eq,0, then
q1 = 147000 !中-活载
q2 = 92000
q3 = 80000
L1 = 7.5
L2 = 30
*elseif,ttype,eq,1,then
q1 = 0
q2 = 0
q3 = 0
L1 = 0
L2 = 0
*else
*msg,ui
error value of ttype(1 or 2),please enter again
*go,: reread
*endif
! 识别桥面过列车的单元,并将单元按节点坐标从小到大排序
lsel,s,loc,y,y_bot,y_top
esll,s ! 选中桥面单元
elmax = elmiqr(0,14) ! 最大单元号
ndmax = ndinqr(0,14) ! 最大节点号
*dim,eldk,,elmax ! 存放桥面过车单元
*dim,ndx,,elmax ! 存放过车单元的节点x坐标(两节点坐标中的较小值)
*dim,nodes,,ndmax ! 存放桥面过车的节点
j=0
*do,i,1,elmax
*if,esel ,gt,0, then
j = j+1
eldk(j) = i
*endif
*enddo ! 记录桥面单元的标号,保存于数组eldk中
eldknum = j ! 总桥面单元数
! 对eldk数组排序,按节点坐标从小到大排列
*do,i,1,eldknum
ii = eldk
ndi = nelem(ii,1)
ndj = nelem(ii,2)
i1 = nx(ndi)
i2 = nx(ndj)
ndx = i1
*if,i2,lt,i1,then
ndx = i2
*endif
*enddo
*do,i,1,eldknum-1
*do,j,i+1,eldknum
*if,ndx(j),lt,ndx ,then
temp = ndx
ndx = ndx(j)
ndx(j) = temp
temp = eldk
eldk = eldk(j)
eldk(j)= temp
*endif
*enddo
*enddo
nsle,s ! 选中桥面单元的节点
j=0
*do,i,1,ndmax
*if,nsel ,gt,0,then
j=j+1
nodes(j)=i
*endif
*enddo
ndnum = j ! 总桥面节点数
xmin = nx(nelem(eldk(1),1))
temp = nx(nelem(eldk(1),2))
*if,temp,lt,xmin,then
xmin = temp
*endif
xmax = nx(nelem(eldk(eldknum),1))
temp = nx(nelem(eldk(eldknum),2))
*if,temp,gt,xmax,then
xmax = temp
*endif
qth = xmax -xmin ! 桥面过车长度
lth1 = 2*qth + L1+L2 +182.5
lth = (nint(lth1/5)+8)*5 ! 总长度
dim = lth/5 +1
*dim,disy,table,dim,ndnum ! 桥面节点的y方向的位移值
*dim,disy1,table,dim,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,,20 ! 前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,dmmaxel,table,dim,20 ! 前20各轴力变化最大单元的弯矩
*dim,fmaxel,table,dim,20 ! 前20各轴力最大单元的轴力
*dim,mmaxel,table,dim,20 ! 前20各轴力最大单元的弯矩
eplot ! plot element
allsel
time = 0
*do,kk,0,lth,5,
tim
|<< << < 1 2 3 4 5 6 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:火车  过桥  动态  加载  实例  部分
相关资讯
热点资讯
推荐资讯

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