悬索桥受力分析实例【转】
/filname,cable,1
keyw,pr_struc,1
/prep7
et,1,link10
et,2,beam3
r,1,1,1e-5
r,2,0.02,1e-5
r,3,0.5,1/24,1
mp,ex,1,25e10
mp,prxy,1,0.1
mp,dens,1,1000
mp,ex,2,30e10
mp,prxy,2,0.1
mp,dens,2,1000
!生成节点
n,1
n,2,10,0.200026668
n,3,20,0.800426758
n,4,30,1.802161037
n,5,40,3.206832495
n,6,50,5.016688905
n,7,60,7.234626423
n,8,70,9.864194224
n,9,80,12.90960018
n,10,90,16.37571758
n,11,100,20.26809296
n,12,10
n,13,20
n,14,30
n,15,40
n,16,50
n,17,60
n,18,70
n,19,80
n,20,90
n,21,100
!生成单元
!生成索单元
*do,i,1,10,1
e,i,i+1
*enddo
!生成吊杆单元
real,2
*do,i,2,10,1
e,i,i+10
*enddo
!生成梁单元
type,2
mat,2
real,3
e,1,12
*do,i,12,20,1
e,i,i+1
*enddo
!生成全桥节点
nsym,x,21,all
esym,,21,all
!合并、压缩节点和单元
asell
nummrg,all
numcmp,all
/replot
!加载与求解
!由于link10是空间单元，所以如果有节点与此单元相连，那么应该对那些节点给与适当的约束，这里应该给与z向约束。悬索为铰约束，大梁为简支
!大梁约束
d,41,,,,,,ux,uy
d,21,uy
!悬索端约束
d,11,,,,,,ux,uy
d,31,,,,,,ux,uy
!link10单元节点设置z向约束
nsel,s,node,,21,42,21
nsel,invert
d,all,uz
allsel
!设定重力场
acel,,10
!运用单元生死模拟悬索的找形和加劲梁的施工过程
!悬索的找形
!牛顿拉普森选项
nropt,full
!杀死除悬索外的其他单元
esel,s,real,,2
esel,a,real,,3
ekill,all
allsel
!设置求解选项
/solu
nlgeom,on
NSUBST,10,20,5
time,1
!设置应力刚化
sstif,on
solve
save,,,,all
!画y向节点位移图，检查悬索是否已达到零位移。（只选活单元）
!esel,s,live
!UPGEOM,1,LAST,LAST,'cable','rst',' '!更新几何形状（使其位移在同侧且较小即可）
!加劲梁的施工模拟
!施工过程为从跨中向两端对称进行，每一个工期左右各架设一段加劲梁（beam3单元）。由于在施工过程中加劲梁各单元之间是铰接的，所以可以把它当作集中荷载作用在悬索上，而这里如直接用激活单元的办法来模拟是不合适的，每段加劲梁的自重W=50000，两端节点各承受25000.只有当加劲梁已经完成安装时，才激活全部单元，从而达到体系转化目的
!中间段梁的架设
f,1,fy,-50000
f,2,fy,-25000
f,22,fy,-25000
!求解
allsel
time,2
solve
!下一阶段
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,2,fy,-25000
f,3,fy,-25000
f,22,fy,-25000
f,23,fy,-25000
time,3
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,3,fy,-25000
f,4,fy,-25000
f,24,fy,-25000
f,23,fy,-25000
time,4
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,5,fy,-25000
f,4,fy,-25000
f,24,fy,-25000
f,25,fy,-25000
time,5
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,5,fy,-25000
f,6,fy,-25000
f,26,fy,-25000
f,25,fy,-25000
time,6
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,7,fy,-25000
f,6,fy,-25000
f,26,fy,-25000
f,27,fy,-25000
time,7
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,7,fy,-25000
f,8,fy,-25000
f,28,fy,-25000
f,27,fy,-25000
time,8
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,9,fy,-25000
f,8,fy,-25000
f,28,fy,-25000
f,29,fy,-25000
time,9
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,9,fy,-25000
f,10,fy,-25000
f,30,fy,-25000
f,29,fy,-25000
time,10
solve
dofsel,s,,fy
fcum,add,1,
dofsel,all
f,10,fy,-25000
f,30,fy,-25000
time,11
solve
!悬索桥结构的体系转换
!施工完成后，梁将参与结构的受力，此时应删除由于梁的作用而加在索上的集中力，用梁本身的重量来代替，因此需要激活所有吊杆和梁单元
!删除所有悬索上的集中力，包括铰支节点上的集中力
nsel,s,node,,1,10,1
nsel,a,node,,22,30,1 !有问题，计算不收敛，可能需要重启动
fdele,all,fy
esel,s,live
esel,invert
ealive,all
nsel,all
esel,all
time,12
solve
有没有桁架式起重机的资料啊？模型及建模过程 本帖最后由 Chelsea 于 2011-3-15 17:30 编辑
kakakey 发表于 2011-3-15 17:02 static/image/common/back.gif
有没有桁架式起重机的资料啊？模型及建模过程
参考http://www.chinavib.com/thread-100533-1-1.html
http://www.chinavib.com/forum-redirect-goto-findpost-ptid-8989-pid-60157-fromuid-78822.html {:{39}:}这个东西很需要，下载下来了。准备好好研究下。悬索斜拉都有非线性问题，看看这个能不能研究下。
页:
[1]