曲线桥梁如何显示弯矩

chillyer 发表于 2007-8-1 15:17:19 | 显示全部楼层 | 阅读模式
各位老大,我用etable定义mx,my,怎么显示的不对啊,还有扭矩应该如何定义啊。
精彩评论倒序浏览

4689查看6评论

水幽寒 发表于 2007-8-1 22:09:50
你用的是什么单元建立的模型,请贴出部分命令流,或者把问题详细的说清楚
举报 回复
chillyer 发表于 2007-8-2 08:38:44
fini
/clear
/title,three span continous bridge
/prep7
rad=200
h=1.6
a=18.5
b=12.8
b0=0.50!中隔板厚度
a0=3.60!箱室宽度
csys,1
k,,rad-b/2
k,,rad-a0/2-b0/2
k,,rad+a0/2+b0/2
k,,rad+b/2
k,,rad-b/2-0.2,,h
k,,rad-a0/2-b0/2,,h
k,,rad+a0/2+b0/2,,h
k,,rad+b/2+0.2,,h
k,,rad-a/2,,h
k,,rad+a/2,,h

kgen,2,1,10,1,,11.55/rad*180/3.14,,10
kgen,9,11,20,1,,1.25/rad*180/3.14,,10
kgen,3,91,100,1,,11.55/rad*180/3.14,,10
kgen,5,111,120,1,,1.25/rad*180/3.14,,10

et,1,shell63
mp,dens,1,2500
mp,ex,1,3.5e10
mp,prxy,1,0.1667
r,1,0.97
a,1,2,3,4,8,7,6,5
aatt,1,1,1
esize,1
amesh,all
asel,none
r,2,2
a,51,52,53,54,58,57,56,55
aatt,1,2,1
esize,1
amesh,all
asel,none
r,3,0.3
a,101,102,103,104,108,107,106,105
aatt,1,3,1
esize,1
amesh,all
asel,none
r,4,0.97
a,151,152,153,154,158,157,156,155
aatt,1,4,1
esize,1
amesh,all
asel,none

r,5,0.7
*do,i,1,15
a,5+(i-1)*10,5+10*i,1+10*i,1+(i-1)*10
a,8+(i-1)*10,8+10*i,4+10*i,4+(i-1)*10
aatt,1,5,1
esize,1
amesh,all
asel,none
*enddo
r,6,0.22
*do,k,1,3,1
a,k,k+1,k+1+10,k+10
aatt,1,6,1
esize,1.0
amesh,all
asel,none   
*enddo

!底板厚度渐变
*dim,hd,array,5
hd(1)=0.42-(0.42-0.22)/4*4
hd(2)=0.42-(0.42-0.22)/4*3
hd(3)=0.42-(0.42-0.22)/4*2
hd(4)=0.42-(0.42-0.22)/4*1
hd(5)=0.42
*do,i,1,4,1
r,10+i,hd(i),hd(i),hd(i+1),hd(i+1)   !定义渐变的底板厚度
*enddo
*do,i,1,4,1
*do,k,1,3,1
  a,k+i*10,k+1+i*10,k+1+(i+1)*10,k+(i+1)*10
  aatt,1,10+i,1
esize,1.0
amesh,all
asel,none   
*enddo
*enddo
*do,i,5,8,1
*do,k,1,3,1
a,k+1+(i+1)*10,k+(i+1)*10,k+i*10,k+1+i*10
aatt,1,19-i,1
esize,1.0
amesh,all
asel,none   
*enddo
*enddo
*do,i,10,11,1
*do,k,1,3,1
a,k+(i-1)*10,k+1+(i-1)*10,k+1+i*10,k+i*10
aatt,1,6,1
esize,1.0
amesh,all
asel,none   
*enddo
*enddo
*do,i,11,14,1
*do,k,1,3,1
a,k+i*10,k+1+i*10,k+1+(i+1)*10,k+(i+1)*10
aatt,1,i,1
esize,1.0
amesh,all
asel,none   
*enddo
*enddo

r,7,0.5
*do,i,1,15
a,6+(i-1)*10,2+(i-1)*10,2+10*i,6+10*i
a,7+(i-1)*10,3+(i-1)*10,3+10*i,7+10*i
aatt,1,7,1
esize,1.0
amesh,all
asel,none
*enddo

r,8,0.23
*do,i,1,15
a,9+(i-1)*10,5+(i-1)*10,5+10*i,9+10*i
a,5+(i-1)*10,6+(i-1)*10,6+10*i,5+10*i
a,6+(i-1)*10,7+(i-1)*10,7+10*i,6+10*i
a,7+(i-1)*10,8+(i-1)*10,8+10*i,7+10*i
a,8+(i-1)*10,10+(i-1)*10,10+10*i,8+10*i
aatt,1,8,1
esize,1.0
amesh,all
asel,none
*enddo
allsel,all
nsym,y,100000,all
esym,y,100000,all
!支座
nsel,s,,,54
nsel,a,,,61
nsel,a,,,150
nsel,a,,,157
nsel,a,,,100054
nsel,a,,,100061
nsel,a,,,100150
nsel,a,,,100157

d,all,uy,,,,,,uz,rotx,rotz
allsel,all
numcmp,all
nummrg,all

!工况1  一期恒载
/solu
acel,0,0,9.8
solve
举报 回复
chillyer 发表于 2007-8-2 08:40:30
命令流如上,请看有没有问题,shell63建模,etable定义输出mx,my和mxy
输出结果好像不对哦,是我理解错了,还是命令流有错
请老大指正!
举报 回复
jerlos 发表于 2008-1-21 08:50:33
曲线桥要正确输出内力结果是要定义单元坐标的
举报 回复
qlql123 发表于 2009-9-17 19:48:05
楼主搞清楚怎么显示弯矩了吗?有什么收获分享一下
举报 回复
fengzi13 发表于 2009-10-14 14:09:28
实体单元的弯矩显示要采用映射实现。只有梁单元才能直接显示扭矩。
举报 回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00