用自定义截面梁建立模型.doc

用自定义截面梁建立模型.txt你看得见我打在屏幕上的字，却看不到我掉在键盘上的泪！自己选择45仰视别人，就休怪他人135俯视着看你。先画好梁的截面图形，用plane82划分平面单元（只能用plane82）后，采用secwrite命令写为截面特性文件name.sect，有许多不同截面或者是变截面时可以一次在一个辅助程序里把所有截面特性文件都建立好，再在主程序中用secread读入后建立模型。 ！下面是变截面箱梁，建立数组存放个截面处箱梁控制数据。通过一个循环生成全部截面特性文件。 !/CLEAR,START *dim,H,array,100, *dim,Htop,array,100 *dim,Hbot,array,100 *dim,Hfu,array,100 *dim,name,char,100 ! 各个截面 H(01)= 3 $ Htop(01)=0.28 $ Hbot(01)=0.32 $ Hfu(01)=0.45 $name(01)=A3 ! A3截面 H(02)= 3.082 $ Htop(02)=0.28 $ Hbot(02)=0.328 $ Hfu(02)=0.45 $name(02)=A5 !L5 H(03)= 3.326 $ Htop(03)=0.28 $ Hbot(03)=0.350 $ Hfu(03)=0.45 $name(03)=A7 !L5 H(04)= 3.735 $ Htop(04)=0.28 $ Hbot(04)=0.389 $ Hfu(04)=0.45 $name(04)=A9 !L7 H(05)= 4.226 $ Htop(05)=0.28 $ Hbot(05)=0.434 $ Hfu(05)=0.45 $name(05)=A11 !L9 H(06)= 4.746 $ Htop(06)=0.28 $ Hbot(06)=0.458 $ Hfu(06)=0.45 $name(06)=A13 !L11 H(07)= 5.358 $ Htop(07)=0.28 $ Hbot(07)=0.540 $ Hfu(07)=0.45 $name(07)=A15 !L13 H(08)= 6. $ Htop(08)=0.28 $ Hbot(08)=0.800 $ Hfu(08)=0.45 $name(08)=A17 !L15 ! H(9)= 7.5 $ Htop(9)=0.45 $ Hbot(9)=0.95 $ Hfu(9)=0.65 $name(9)=A0 !0号块 ! /prep7 et,1,plane82 ! *do,i,1,9,1 k,1,-22.5/2,0 k,2,22.5/2,0 k,3,22.5/2,-0.2 k,4,22.5/2-4,-0.45 k,5,22.5/2-5.5,-0.65 k,6,22.5/2-5.5,-H(i) k,7,-22.5/2+5.5,-H(i) k,8,-22.5/2+5.5,-0.65 k,9,-22.5/2+4,-0.45 k,10,-22.5/2,-0.2, a,1,2,3,4,5,6,7,8,9,10 k,11,-22.5/2+5.5+Hfu(i),-0.65, k,12,-3.75,-0.45, k,13,0,-Htop(i) k,14,3.75,-0.45, k,15,22.5/2-5.5-Hfu(i),-0.65 k,16,22.5/2-5.5-Hfu(i),-H(i)+Hbot(i), k,17,-22.5/2+5.5+Hfu(i),-H(i)+Hbot(i), a,11,12,13,14,15,16,17 asba,1,2 esize,1, mshkey,2 amesh,all secwrite,name(i),sect !写入文件 .sect ACLEAR, all ADELE, all , ,1 LDELE, all, , ,1 NUMSTR,NODE,1, NUMSTR,ELEM,1, NUMSTR,KP,1, NUMSTR,LINE,1, NUMSTR,AREA,1, NUMSTR,VOLU,1, *enddo ! A1 k,1,-22.5/2,0 k,2,22.5/2,0 k,3,22.5/2,-0.65 k,4,22.5/2-4,-0.65 k,5,22.5/2-5.5,-0.65 k,6,22.5/2-5.5,-3 k,7,-22.5/2+5.5,-3 k,8,-22.5/2+5.5,-0.65 k,9,-22.5/2+4,-0.65 k,10,-22.5/2,-0.65 a,1,2,3,4,5,6,7,8,9,10 k,11,-22.5/2+5.5+1.6,-0.65, k,12,-3.75,-0.65, k,13,0,-0.65 k,14,3.75,-0.65, k,15,22.5/2-5.5-1.6,-0.65 k,16,22.5/2-5.5-1.6,-3+0.8, k,17,-22.5/2+5.5+1.6,-3+0.8, a,11,12,13,14,15,16,17 asba,1,2 !LESIZE,ALL,1, , , ,1, , ,1, esize,1, mshkey,2 amesh,all secwrite,A1,sect !写入文件 .sect ACLEAR, all ADELE, all , ,1 LDELE, all, , ,1 NUMSTR,NODE,1, NUMSTR,ELEM,1, NUMSTR,KP,1, NUMSTR,LINE,1, NUMSTR,AREA,1, NUMSTR,VOLU,1, finish !11 !调入和验证- /clear /PREP7 k,1 k,2,100 k,3,500,200 l,1,2 et,1,beam189 mp,ex,1,2e5 mp,prxy,1,0.3 lesize,all,20 latt,1,1,3,1 lmesh,all sectype,1,beam,mesh !定义并读入截面 secread,A1,sect, ,mesh secplot,1,1 finish /solu dk,1,all fk,2,fy,-1000 solve /post1 /plopts,info,3 pldisp,1 etable,mi,smisc,2 etable,mj,smisc,15 plls,mi,mj /eof ! 下面是主程序 /clear /PREP7 ET,1,BEAM44 mp,ex,1,3.575e7 $ mp,prxy,1,1/6 mp,dens,1,2.5, ! 建立桥墩的截面特性， SECTYPE, 100, BEAM, RECT, dun1, 0 SECOFFSET, CENT SECDATA,1.5,11.5,0,0,0,0,0,0,0,0 SECTYPE, 101, BEAM, RECT, dun2, 0 SECOFFSET, CENT SECDATA,0.8,11.5,0,0,0,0,0,0,0,0 ! ! 开始建立结构单元，顺桥向为X轴，竖向为Y轴， !主墩1 n,1000,-10,0,0 N,1,60-2.5,-70.4+(60-2.5)*0.03,0, SECNUM, 100 n,5,60-2.5,-70.4+(60-2.5)*0.03+22, n,10,60-2.5,-70.4+(60-2.5)*0.03+22+24, n,15,60-2.5,0+(60-2.5)*0.03, fill,1,5, fill,5,10, fill,10,15, *do,i,1,14,1 e,i,i+1, *enddo NGEN,2,20,all, , ,5,5*0.03, ,1, EGEN,2,20,all, , , , , , , ,5,5*0.03, , SECNUM,101, e,5,25 e,10,30 !主墩2 N,41,166-2.5,-34.24+(166-2.5)*0.03,0, n,45,166-2.5,-34.24+(166-2.5)*0.03+18,0 n,48,166-2.5,(166-2.5)*0.03,0 fill,41,45 fill,45,48 SECNUM, 100 *do,i,41,47,1 e,i,i+1, *enddo NGEN,2,10,41,48,1,5,5*0.03 , ,1, EGEN,2,10,31,37,1 , , , , , ,5,5*0.03, , !A1 sectype,1,beam,mesh, $ SECOFFSET,ORIG, secread,A1,sect,mesh secnum,1 ! n,999,58-2.5,25,0, n,1001,0,0,0 n,1002,4,4*0.03,0 e,1001,1002,999 !L3 sectype,2,beam,mesh, $ SECOFFSET,ORIG, secread,A3,sect,mesh secnum,2 ! n,1003,12,12*0.03,0 e,1002,1003,999 !L5 sectype,3,beam,mesh, $ SECOFFSET,ORIG, secread,A5,sect,mesh secnum,3 ! n,1004,20,20*0.03,0 e,1003,1004,999 !L7 sectype,4,beam,mesh, $ SECOFFSET,ORIG, secread,A7,sect,mesh secnum,4 ! n,1005,28,28*0.03,0 e,1004,1005,999 !L9 sectype,5,beam,mesh, $ SECOFFSET,ORIG, secread,A9,sect,mesh secnum,5 ! n,1006,36,36*0.03,0 e,1005,1006,999 !L11 sectype,6,beam,mesh, $ SECOFFSET,ORIG, secread,A11,sect,mesh secnum,6 ! n,1007,42,42*0.03,0 e,1006,1007,999 !L13 sectype,7,beam,mesh, $ SECOFFSET,ORIG, secread,A13,sect,mesh secnum,7 ! n,1008,48,48*0.03,0 e,1007,1008,999 !L15 sectype,8,beam,mesh, $ SECOFFSET,ORIG, secread,A15,sect,mesh secnum,8 ! n,1009,54,54*0.03,0 e,1008,1009,999 !L17 ! 连接段 sectype,9,beam,mesh, $ SECOFFSET,ORIG, secread,A17,sect,mesh secnum,9 ! e,1009,15,999 !L0 !0号块 sectype,10,beam,mesh, $ SECOFFSET,ORIG, secread,A0,sect,mesh secnum,10 ! e,15,35,999 wpoff,60,0,0 wprot,1.718358001,0,0 ! 旋转3% 的坡度 CSYS,4 n,1499,58-2.5,25,0, NSYM,X,500,1003,1009,1 ESYM, ,500,47,52,1 !L21 secnum,9 ! e,1509,35,999 NGEN,2,1000,1003,1009,1,106, , ,1, EGEN,2,1000,47,52,1, , , , , ,106, , , !L21 e,2009,48,999 !L3 CSYS,0 n,2050,60+53,(60+53)*0.03, secnum,2 e,1503,2050,999 e,2050,2003,999 WPCSYS,-1,0 wpoff,60+53,0,0 wprot,1.718358001,0,0 ! 旋转3% 的坡度 CSYS,4 n,3999,58-2.5,25,0, NSYM,X,3000,1001,1009,1 ESYM, ,3000,45,52,1 csys,0 !L21 secnum,9 ! e,58,4009,999 !L0 secnum,10 ! e,48,58,999 SECNUM,101, e,45,55 /solu D,1,all , D,21,all, D,41,all, D,51,all, d,1001,ux d,1001,uy, d,1001,uz, d,1001,rotx d,4001,ux, d,4001,uy, d,4001,uz, d,4001,rotx, !d,all,uz, ! acel,0,9.8,0, ANTYPE,MODAL ! Cho