空间网架结构优化设计命令流.doc

*CREATE,SCRATCH /filname, RACE /prep7 et,1,link8 A=0.00005 !定义单元属性 r=9 pi=3.1415926535 R,1,A !定义材料常数 MP,EX,1,2.07e11 MP,PRXY,1,0.3 MP,DENS,1,78 !定义材料属性 real,1 * do,i,1,3 * do,j,1,12 n,j+12*(i-1),( r-3*(i-1))*cos(pi/6*(j-1)),( r-3*(i-1))*sin(pi/6*(j-1)) *enddo *enddo * do,i,1,3 * do,j,1,12 n,108+j+12*(i-1),( r-1.5-3*(i-1))*cos(pi/6*(j-1)+pi/12),( r-1.5-3*(i-1))*sin(pi/6*(j-1)+ pi/12),3 *enddo *enddo *do,i,1,3 !上层纬向杆连接 *do,j,1,e,12*(i-1)+j,12*(i-1)+j+1 *enddo *enddo e,1,12 e,13,24 e,25,36 *do,i,1,3 !下层纬向杆连接 *do,j,1,e,108+12*(i-1)+j,108+12*(i-1)+j+1 *enddo *enddo e,109,120 e,121,132 e,133,144 * do,i,1,2 !上层径向杆连接 *do,j,1,e,12*(i-1)+j,12*(i-1)+j+12 *enddo *enddo *do,i,2,11 !上层斜向杆连接 e,,12+i-1 e,I,12+i+1 *enddo E,1,24, E,1,14 E,12,23 E,12,13 *do,i,14,23 !上层斜向杆连接 e,,12+i-1 e,I,12+i+1 *enddo E,13,26, E,13,36 E,24,25 E,24,35 *do,i,25,30 !上层斜向杆连接 e,,i+6 *enddo * do,i,109,132 !下层径向杆连接 e, ,i+12 *enddo * do,i,2,12 !腹杆杆连接 e, ,109+i-1 e,I,109+i-2 *enddo e, 1,109 e,1,120 * do,i,110,119 !腹杆杆连接 e, ,14+i-110 e,I,15+i-110 *enddo e, 13,109 e, 14,109 e, 13,120 e,24,120 * do,i,14,24 !腹杆杆连接 e, ,109+i-1 e,I,109+i-2 *enddo e, 13,121 e,13,132 * do,i,122,131 !腹杆杆连接 e, ,14+i-110 e,I,15+i-110 *enddo e, 25,121 e, 26,121 e, 25,132 e,36,132 * do,i,26,36 !腹杆杆连接 e, ,109+i-1 e,I,109+i-2 *enddo e, 133,25 e,144,25 d,1,all * do,i,1,12,2 * do,j,109,120,3 d,i,all d,j,all *enddo *enddo F,all,FZ,10000 NSEL,ALL FINISH /SOLU ANTYPE,0 SOLVE FINISH /Post1 ! 进入一般后处理模块 ETABLE,VOLU,VOLU ! 体积为内容定义单元表 Ssum *GET,TVOL,SSUM,,ITEM, VOLU ETABLE, smax, LS, 1 !轴向应力建立单元表 Esort,etab, smax !轴向应力排序 *GET, smax1, sort,, max ！取最大应力 *GET, smin1, sort,, min ！取最小应力 smin1=ABS(smin1) PRNSOL,U,z ! 列表显示节点位移UY值 NSORT,U,z,,1 ! 对节点计算结果UY排序 PRNSOL,U,z ! 列表显示节点位移UY值 *GET,DEFL,SORT,,MAX ! 提取排序后的节点位移结果,保存到变量DEFL FINISH *END FINISH /CLEAR, NOSTART /FILENAME, EX11.4 *USE,