马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- /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
复制代码
|