part1 tm0817.docx
《part1 tm0817.docx》由会员分享，可在线阅读，更多相关《part1 tm0817.docx（7页珍藏版）》请在冰豆网上搜索。
part1tm0817
finish
/clear
/prep7
/nerr,0
et,1,82!
定义截面网分单元
*dim,name,char,20!
定义截面名称数组
name
（1）='n01'!
主梁截面号
name
（2）='n02'$name（3）='n03'$name（4）='n04'
name（5）='n05'$name（6）='n06'$name（7）='n07'
name（8）='n08'$name（9）='n09'$name（10）='n10'
name（11）='n11'!
桥墩截面号
!
~~~~~~~~~~~~定义箱梁截面数据~~~~~~~~~~~~~~~~~~~~~~~
*dim,aa,array,10!
截面顶板半宽-翼缘板外挑净宽2.5m
aa
（1）=3.25,3.25,3.25,3.25,3.25,3.25,3.25,3.25,3.25,3.25
*dim,bb,array,10!
截面底板半宽
bb
（1）=3.25,3.25,3.25,3.25,3.25,3.25,3.25,3.25,3.25,3.25
*dim,ttop,array,10!
截面顶板厚度
ttop
（1）=0.55,0.498,0.28,0.28,0.28,0.28,0.28,0.28,0.28,0.55
*dim,tf,array,10!
截面腹板厚度
tf
（1）=0.8,0.771,0.55,0.55,0.55,0.55,0.55,0.55,0.55,0.8
*dim,tbot,array,10!
截面底板厚度
tbot
（1）=1,0.6,0.498,0.431,0.375,0.325,0.289,0.267,0.26,0.6
*dim,height,array,10!
截面箱梁高度
height
（1）=3.5,3.5,3.049,2.754,2.508,2.286,2.127,2.032,2,2
*dim,k8yy,arry,10!
定义截面关键点8的y坐标
k8yy
（1）=-0.85,-0.798,-0.58,-0.58,-0.58,-0.58,-0.58,-0.58,-0.58,-0.58
*dim,k9xx,arry,10!
定义截面关键点9的x坐标
k9xx
（1）=-1.25,-1.25,-1.25,-1.25,-1.25,-1.25,-1.25,-1.25,-1.25,-1.25
!
~~~~~~~~~~~~定义箱梁截面~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*do,i,1,10,1
*set,a,aa（i）!
顶板半宽-翼缘板外挑净宽
*set,b,bb（i）!
底板半宽
*set,tt,ttop（i）!
顶板厚度
*set,ft,tf（i）!
腹板厚度
*set,tbt,tbot（i）!
底板厚度
*set,h,height（i）!
梁高
*set,k8y,k8yy（i）!
K8纵坐标
*set,k9x,k9xx（i）!
K9横坐标
k,1,0,0$k,2,-1.25,0$k,3,-a+ft,0
k,4,-a,0$k,5,-5.75,0$k,6,-5.75,-0.2
k,7,-a,-0.55$k,8,-a+ft,k8y$k,9,k9x,-tt
k,10,0,-tt$k,11,0,-h+tbt$k,12,-a+ft+0.4,-h+tbt
k,13,-a+ft,-h+tbt+0.2$k,14,-a,-h+tbt+0.2
k,15,-a,-h$k,16,-a+ft,-h$k,17,-a+ft+0.4,-h
k,18,0,-h
l,1,10,1,1$l,1,2,4,1$l,2,3,2,1$l,3,4,2,1
l,4,5,4,1$l,4,7,1,1$l,5,6,1,1$l,6,7,4,1
l,8,9,2,1$l,8,3,1,1$l,8,7,2,1$l,7,14,4,1
l,9,2,1,1$l,9,10,4,1$l,8,13,4,1$l,14,13,2,1
l,14,15,2,1$l,13,16,2,1$l,15,16,2,1$l,13,12,1,1
l,16,17,1,1$l,12,17,2,1$l,12,11,4,1$l,17,18,4,1
l,11,18,2,1
a,1,2,9,10$a,2,3,8,9$a,3,4,7,8$a,4,5,6,7
a,7,14,13,8$a,14,13,16,15$a,12,13,16,17
a,11,12,17,18
mshape,0,2d$mshkey,0$amesh,all
arsym,x,all!
对称生成另一半截面
nummrg,all,,,,low$numcmp,all
secwrite,name（i）,,,1!
写入截面
sectype,i,beam,mesh,name（i）
secoffset,user,0,0!
截面偏移
secread,name（i）,,,mesh$secplot,i,
allsel
aclear,all$adele,all,,,1
*enddo
!
~~~~~~~~~~~建立桥墩截面（6.5×1.8m矩形,截面号11）~~~~~~~~~
rectng,-3.25,3.25,-0.9,0.9
allsel
lesize,all,,,4
mshape,0,2d$mshkey,1$amesh,all
secwrite,n11,,,1
sectype,11,beam,mesh,n11
secread,n11,,,mesh$secplot,11
allsel
aclear,all$adele,all,,,1
etdele,1!
删除辅助单元类型
!
~~~~~~~~~~~~定义所有材料特性~~~~~~~~~~~~~~~~~~~~~~~
et,1,beam188!
主梁C50
keyopt,1,3,2
MP,EX,1,3.45e10$mp,dens,1,2600
MP,PRXY,1,0.2
et,2,beam188!
桥墩C40
keyopt,2,3,2$MP,EX,2,3.25e10
mp,dens,2,2600$MP,PRXY,2,0.2
et,3,mpc184!
定义MPC-184刚性梁
keyopt,3,1,1$keyopt,3,2,0
!
~~~~~~~~~~建立变截面~~~~~~~~~~~~~~~~~~~~~~~~~~~~
SECTYPE,21,TAPER!
左边跨现浇段
SECDATE,10,-65$SECDATE,10,-64.5
SECTYPE,22,TAPER$SECDATE,10,-64.5$SECDATE,10,-64
SECTYPE,23,TAPER$SECDATE,10,-64$SECDATE,9,-62
SECTYPE,24,TAPER$SECDATE,9,-62$SECDATE,9,-61
SECTYPE,25,TAPER!
左边跨合龙段
SECDATE,9,-61$SECDATE,9,-59
*do,i,1,4,1!
左边跨悬浇段3、4、5、6
SECTYPE,25+i,TAPER
SECDATE,10-i,-59+（i-1）*4$SECDATE,9-i,-55+（i-1）*4
*enddo
*do,i,1,2,1!
左边跨悬浇段2、1
SECTYPE,29+i,TAPER
SECDATE,6-i,-43+（i-1）*3.5$SECDATE,5-i,-39.5+（i-1）*3.5
*enddo
SECTYPE,32,TAPER!
左墩0号段
SECDATE,3,-36$SECDATE,2,-31.5
SECTYPE,33,TAPER$SECDATE,2,-31.5$SECDATE,1,-30.9
*do,i,1,2,1
SECTYPE,33+i,TAPER$SECDATE,1,-30.9+（i-1）*0.9
SECDATE,1,-30+（i-1）*0.9$*enddo
SECTYPE,36,TAPER$SECDATE,1,-29.1$SECDATE,2,-28.5
SECTYPE,37,TAPER$SECDATE,2,-28.5$SECDATE,3,-24
*do,i,1,2,1!
中跨左悬浇段1、2
SECTYPE,37+i,TAPER$SECDATE,2+i,-24+（i-1）*3.5
SECDATE,3+i,-20.5+（i-1）*3.5$*enddo
*do,i,1,4,1!
中跨左悬浇段3、4、5、6
SECTYPE,39+i,TAPER$SECDATE,4+i,-17+（i-1）*4
SECDATE,5+i,-13+（i-1）*4$*enddo
SECTYPE,44,TAPER!
中跨合拢段
SECDATE,9,-1$SECDATE,9,0
SECTYPE,45,TAPER!
中跨合拢段
SECDATE,9,0$SECDATE,9,1
*do,i,1,4,1!
中跨右悬浇段3、4、5、6
SECTYPE,45+i,TAPER$SECDATE,10-i,1+（i-1）*4
SECDATE,9-i,5+（i-1）*4$*enddo
*do,i,1,2,1!
中跨右悬浇段1、2
SECTYPE,49+i,TAPER$SECDATE,6-i,17+（i-1）*3.5
SECDATE,5-i,20.5+（i-1）*3.5$*enddo
SECTYPE,52,TAPER!
0号段
SECDATE,3,24$SECDATE,2,28.5
SECTYPE,53,TAPER!
0号段
SECDATE,2,28.5$SECDATE,1,29.1
sectype,54,taper$secdata,1,29.1$secdata,1,30
sectype,55,taper$secdata,1,30$secdata,1,30.9
sectype,56,taper$secdata,1,30.9$secdata,2,31.5
sectype,57,taper$secdata,2,31.5$secdata,3,36
sectype,58,taper!
右边跨悬浇段1、2
secdata,3,36$secdata,4,39.5
sectype,59,taper$secdata,4,39.5$secdata,5,43
*do,i,1,4,1!
右边跨悬浇段3、4、5、6
sectype,59+i,taper$secdata,4+i,43+（i-1）*4
secdata,5+i,47+（i-1）*4$*enddo
sectype,64,taper!
右边跨合龙段
secdata,9,59$secdata,9,61
sectype,65,taper!
右边跨现浇段
secdata,9,61$secdata,9,62
sectype,66,taper$secdata,9,62$secdata,10,64
sectype,67,taper$secdata,10,64$secdata,10,64.5
sectype,68,taper$secdata,10,64.5$secdata,10,65
!
~~~~~~~~~~~~~~~~建立主梁及桥墩~~~~~~~~~~~~~~~~~~~~~~~
k,1,-65,0!
桥梁始点
k,2,-64.5!
左跨边支点
k,3,-64,0$k,4,-62,0$k,5,-61,0$k,6,-59,0
*do,i,1,4$k,6+i,-59+i*4$*enddo
k,11,-39.5,0$k,12,-36.0,0$k,13,-31.5,0
k,14,-30.9,0
k,15,-30.0,0!
左墩中点
k,16,-29.1,0$k,17,-28.5,0$k,18,-24,0
k,19,-20.5,0$k,20,-17,0
*do,i,1,4$k,20+i,-17+i*4$*enddo
k,25,0,0!
桥中点
k,26,1,0
*do,i,1,4$k,26+i,1+i*4$*enddo
k,31,20.5,0$k,32,24,0$k,33,28.5
k,34,29.1
k,35,30!
右墩中点
k,36,30.9$k,37,31.5$k,38,36$k,39,39.5$k,40,43
*do,i,1,4$k,40+i,43+i*4$*enddo
k,45,61$k,46,62$k,47,64
k,48,64.5!
右跨边支点
k,49,65!
桥梁终点
!
指定线的单元属性
*do,i,1,48$l,i,i+1$*enddo
*do,i,1,48
allsel
lsel,s,line,,i$latt,1,,1,,,,20+i$*enddo
allsel
lesize,all,,,1
!
建立桥墩
k,50,-30,,-3.5$k,51,-30,,-19.5
l,50,51,8,1!
左墩
k,52,30,,-3.5$k,53,30,,-19.5
l,52,53,8,1!
右墩
lsel,s,line,,49,50,1
latt,2,,2,,,,11!
赋予材料属性2，单元类型2，截面11
!
墩梁固结单元
l,15,50,1,1!
左墩处
l,35,52,1,1!
右墩处
lsel,s,line,,51,52,1
latt,,,3!
指定墩梁固结处属性为mpc-184刚性梁
!
全桥划分单元
lsel,all$LMESH,all
allsel
eplot
/VUP,1,Z$/ESHAPE,1.0!
显示单元形状
/VIEW,1,,-1$/replot
nummrg,all$numcmp,all
save!
保存
升级会员 