本帖最后由 印记央央 于 2012-11-13 16:54 编辑
我用命令流来说明吧
不知道为何我纵向的划分对了 竖向的就是不行,明明命令流都差不多的
Finish /cle /filn,watermelon,1 !柱尺寸 Lp=6250 Hp=700 Bp=240 Tpw=16 Tpf=10 !次梁尺寸 lg=7600 hg=450 bg=160 tgw=8 tgf=6 !螺栓尺寸 Bolt0=23.5 Bolt=22 Nut=35 Nut_h=14 S1=80 !线距 Sa=50 !水平方向端距 Sb=50 !竖直方向端距 Sc=55 !加劲肋离次梁翼缘的距离 /prep7 et,1,solid65,,,,,,,1 !定义混凝土板单元为soild65 et,2,solid45 et,3,combin39 et,4,solid95 !过渡金字塔单元 et,5,solid92 !螺栓单元 keyopt,4,11,1 keyopt,3,2,0 keyopt,3,3,3 !沿Z轴的平动自由度 keyopt,1,1,0 !考虑大变形 keyopt,1,7,1 !考虑应力松弛,为半脆性开裂,有助于计算收敛 mp,ex,1,30400 !定义混凝土材料的弹性模量,泊松比 mp,prxy,1,0.2 fc=32$ft=2.64 tb,concr,1$tbdata,,0.5,0.95,ft,-1 !关闭压碎 tb,miso,1,,11 !定义混凝土应力应变曲线,采用MISO模型 tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64$tbpt,,0.0010,fc*0.75$tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99 tbpt,,0.002,fc$tbpt,,0.0033,fc mp,ex,2,2.0e5 !定义钢筋的弹性模量,泊松比 mp,prxy,2,0.3 tb,biso,2 !钢筋的应力应变曲线,用BISO模型 tbdata,,290,0 mp,ex,3,2.06e5 !定义材料三,钢梁的模量,泊松比 mp,prxy,3,0.3 tb,bkin,3 tbdata,1,310 mp,ex,4,2.06e5 mp,prxy,4,0.3 tb,bkin,3 tbdata,1,400 !假设的 !tb,mkin,3 !tbtemp,,strain !tbdata,,0.001,0.0015,0.025 !tbtemp,0.0 !tbdata,,206,310,420 !tb,mkin,4 !tbtemp,,strain !tbdata,,0,003,0.006,0.008 !tbtemp,0.0 !tbdata,,618,940,1040 tbplot pi=acos(-1) v1=0.25*pi*6*6*6/(400*50) !上部区纵筋的体积配筋率 v2=0.25*pi*6*6*3.5/(400*50) !下部区纵筋的体积配筋率 v4=0.25*pi*6*6*24.5/(1920*50) !水平箍筋的体积配筋率 r,2,2,v1,0,90,2,v4,0,0 !定义拉区混凝土单元的实常数 r,3,2,v2,0,90,2,v4,0,0 !定义压区混凝土单元的实常数 r,1 r,4 r,5,0,0,0.19,20099,0.5,29291 rmore,1.0,38012,1.5,43978,2,48550 rmore,2.5,52244,3,55323 !创建局部坐标系 local,11,0,,,lp/2 local,12,0,tpf/2+(bp-tpf)/4,,lp/2+tgf/2 local,13,0,bp/2+(2*sa+10)/2,(hp-hg)/2, lp/2+tgf/2 local,14,0,bp/2+10, (hp-hg)/2, lp/2 local,15,0,bp/2+10+sa,hp/2-sc-sb, lp/2+tgf/2+8 local,16,0,lg/2,hp/2,lp/2,,-90 !创建主梁 Csys,0 K,1,0,-hp/2$k,2,bp/2,-hp/2$k,3,bp/2,-hp/2+tpw$k,4,tpf/2,-hp/2+tpw K,5,tpf/2,hp/2-tpw$k,6,bp/2,hp/2-tpw$k,7,bp/2,hp/2$k,8,0,hp/2 A,1,2,3,4,5,6,7,8 Vext,1,,,,,lp Allsel Wpcsys,,12 Blc5,,, (bp-tpf)/2,hp-2*tpw,8 Wpcsys,,13 Blc5,,,2*sa+10,hg-2*sc,8 vadd,2,3 !创建次梁 Wpcsyc,-1 Csys,14 K,33,0,-hg/2 K,34,0,-hg/2,bg/2 K,35,0,-hg/2+tgw,bg/2 K,36,0,-hg/2+tgw,tgf/2 K,37,0,hg/2-tgw,tgf/2 K,38,0,hg/2-tgw,bg/2 K,39,0, hg/2,bg/2 K,40,0, hg/2 A,33,34,35,36,37,38,39,40 Vext,11,,,lg-bp/2-10 Vsymm,z,2 Vadd,2,3 allsel Cm,v1,volu !创建高强螺栓 Wpcsys,,15 Cyl4,0,0,bolt0/2,,,,-(8+tgf)$vgen,2,2,,,,-s1$vgen,2,3,,,,-s1$vgen,2,6,,,,-s1 Vsel,u,,,v1 cm,v2,volu allsel vsbv,v1,v2 cm,v3,volu wpcsys,,15 cylind,bolt/2,nut/2,0,nut_h vgen,2,2,,,,,-(tgf+8+nut_h) wpoff,,,nut_h cylind,0,bolt/2,-(2*nut_h+8+tgf) vsel,u,,,v3 cm,b1,volu vgen,2,b1,,,,-s1 cm,b2,volu vgen,2,b2,,,,-2*s1 allsel vsel,u,,,v3 vglue,all cm,bolt,volu allsel wpcsys,,16 Blc5,,,lg,lp,100 !**********************网格划分*************************! !划分高强螺栓 Allsel$cmsel,s,bolt Vatt,4,6,5,0 Mshape,1,3d Mshkey,0 Vmesh,all Allsel Vsel,s,,,9 Aslv$lsla Wpcsys,14$wpoff,2*sa$wpro,,,90$vsbw,all Wpoff,,,70$vsbw,all Wpoff,-tgf/2$wpro,,,90$vsbw,all Wpoff,,,tgf/2$vsbw,all Wpoff,,,tgf/2$vsbw,all Wpoff,,hg/2-tgw$wpro,,90$vsbw,all Wpoff,,,sc-tgw+sb+s1/2$vsbw,all Wpoff,,,s1$vsbw,all Wpoff,,,s1$vsbw,all Wpoff,,,s1/2+sb+sc-tgw$vsbw,all Cm,liang,volu Cm,liang_line,line Csys,14 Vsel,r,loc,x,2*sa+70,7470$aslv$lsla Cm,liang_ling1,line Lsel,r,length,,tgw$lsel,a,length,,tgf/2$lesize,all,,,1 Cmsel,s,liang_line$Lsel,r,length,,7300$lesize,all,,,147 Cmsel,s,liang_line$Lsel,r,length,,s1$lesize,all,,,2 Cmsel,s,liang_line$Lsel,r,length,, s1/2+sb+sc-tgw $lesize,all,,,3 !这两条命令就是将图中最右边的紫色部分所指的2段和3段 Cmsel,s,liang_line$lsel,r,length,,(bg-tgf)/2$lesize,all,,,3 Vatt,3,4,2,0 !solid45 Mshape,0,3D Mshkey,1 Vmesh,all Allsel$cmsel,s,liang Vsel,r,loc,x,0,2*sa+70$vsel,r,,,20,22,1$vsel,a,,,25,28,1$vsel,a,,,30,31$vsel,a,,,7,15,8 vsel,a,,,44,45$vsel,a,,,52,56,4$aslv$lsla Lsel,r,,,67,97,30$Lsel,a,,,156,254,98$Lsel,a,,,304,443,139$lesize,all,,,3 Lsel,a,,,367,374,7$ Lsel,a,,,353,447,94$Lsel,a,,,410,450,40$lesize,all,,,5 !这两条命令就是梁端上下翼缘划分,即上面图中最上面大红色标出的5段和3段 Vatt,3,4,2,0 !solid45 Mshape,0,3D Mshkey,1 Vmesh,all Vsel,s,,,48,49$vsel,a,,,9,12,3$vsel,a,,,54,55$vsel,a,,,40,41 Vsel,a,,,60,61$vsel,a,,,35,36$vsel,a,,,68,69$vsel,a,,,64,65 aslv$lsla$cm,liang_line2,line Lsel,r,radius,,bolt0/2$cm,r1,line$lsel,r,,,75$lsel,a,,,81,83,1$lsel,a,,,105,108,1 Lsel,a,,,115,118,1$lsel,a,,,125,128,1$lesize,all,,,4$cm,r11,line Cmsel,s,r1$cmsel,u,r11$lesize,all,,,8 Cmsel,s,liang_line2$lsel,r,length,,s1$lesize,all,,,8 Cmsel,s,liang_line2$lsel,r,length,,s1/2+sb+sc-tgw$lesize,all,,,10 !这两条命令就是图中左边紫色的划分,但是未成功 Cmsel,s,liang_line2$ lsel,r,length,,2*sa$lesize,all,,,10 Cmsel,s,liang_line2$ lsel,r,length,,70$lesize,all,,,7 !这两条命令就是图中中间红色的划分,已经成功了 Vatt,3,4,4,0 !solid95 Mshape,1,3D Mshkey,0 Vmesh,all Tchg,solid95,solid92
|