《用自定义截面梁建立模型》由会员分享，可在线阅读，更多相关《用自定义截面梁建立模型（7页珍藏版）》请在金锄头文库上搜索。
1、用自定义截面梁建立模型.txt 你看得见我打在屏幕上的字，却看不到我掉在键盘上的泪！ 自己选择 45仰视别人，就休怪他人 135俯视着看你。先画好梁的截面图形，用 plane82 划分平面单元（只能用 plane82）后，采用 secwrite 命令写为截面特性文件 name.sect，有 许多不同截面或者是变截面时可以一次在一个辅助程序里把所有截面特性文件都建立好，再 在主程序中用 secread 读入后建立模型。 ！下面是变截面箱梁，建立数组存放个截面处箱梁控制数据。通过一个循环生成全部截面特 性文件。 !/CLEAR,START *dim,H,array,100, *dim,Htop,a
2、rray,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)=
3、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)=
4、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、.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)+
6、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
7、/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+
8、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
9、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 ! 下面是主程序
10、 /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
11、,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,
12、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, $ SECOFF
13、SET,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
14、,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, secr
15、ead,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, secre
