能不能指点一下我这为题出在哪?也是同样的问题
finish
/clear
/title,the optimization of the grid structure
/replot
*set,t4,8
*set,t5,6
*set,t6,10
*set,md,7.8e-9
*set,pi,acos(-1)
*set,d1,114 !杆1-6直径和壁厚
*set,t1,4
*set,d2,88.5
*set,d3,60
*set,t3,3.5
*set,d4,140
*set,t4,8
*set,d5,114
*set,t5,6
*set,d6,159
*set,t6,10
*set,area1,pi*(d1*d1/4-(d1/2-t1)*(d1/2-t1)) !杆1-6截面面积
*set,area2,pi*(d2*d2/4-(d2/2-t1)*(d2/2-t1))
*set,area3,pi*(d3*d3/4-(d3/2-t3)*(d3/2-t3))
*set,area4,pi*(d4*d4/4-(d4/2-t4)*(d4/2-t4))
*set,area5,pi*(d5*d5/4-(d5/2-t5)*(d5/2-t5))
*set,area6,pi*(d6*d6/4-(d6/2-t6)*(d6/2-t6))
*create,opgrid
/PREP7 !进入前处理模块: 定义模型
ET, 1, LINK8 !定义第一类单元为三维杆单元LINK8
MP, EX, 1, 207E3 !定义第一类材料弹性模量EX
MP, PRXY,1, 0.3 !材料密度
MP, DENS, 1, 7.8e-9 !材料密度
R, 1, area1 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米(114*4=1381)斜杆
R, 2, area2 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米(88.5*4=1061)上弦杆
R, 3, area3 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米(60*3.5=621)下弦杆
R, 4, area4 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米(140*8=3315)斜杆
R, 5, area5 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米(114*6=2034)上弦杆
R, 6, area6 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米(159*10=4678)斜杆
K, 1 ,-5000 ,-5000 , 0 !定义网架的所有结点,单位毫米
K, 2 ,-3000 ,-5000 , 0
K, 3 ,-1000 ,-5000 , 0
K, 4 , 1000 ,-5000 , 0
K, 5 , 3000 ,-5000 , 0
K, 6 , 5000 ,-5000 , 0
K, 7 ,-4000 ,-4000 ,-800
K, 8 ,-2000 ,-4000 ,-800
K, 9 , 0 ,-4000 ,-800
K, 10 , 2000 ,-4000 ,-800
K, 11 , 4000 ,-4000 ,-800
K, 12 ,-5000 ,-3000 , 0
K, 13 ,-3000 ,-3000 , 0
K, 14 ,-1000 ,-3000 , 0
K, 15 , 1000 ,-3000 , 0
K, 16 , 3000 ,-3000 , 0
K, 17 , 5000 ,-3000 , 0
K, 18 ,-4000 ,-2000 ,-800
K, 19 ,-2000 ,-2000 ,-800
K, 20 , 0 ,-2000 ,-800
K, 21 , 2000 ,-2000 ,-800
K, 22 , 4000 ,-2000 ,-800
K, 23 ,-5000 , -1000 , 0
K, 24 ,-3000 , -1000 , 0
K, 25 ,-1000 , -1000 , 0
K, 26 , 1000 , -1000 , 0
K, 27 , 3000 , -1000 , 0
K, 28 , 5000 , -1000, 0
K, 29 ,-4000 , 0,-800
K, 30 ,-2000 , 0 ,-800
K, 31 , 0 , 0 ,-800
K, 32 , 2000 , 0 ,-800
K, 33 , 4000 , 0 ,-800
K, 34 ,-5000 , 1000 , 0
K, 35 ,-3000 , 1000 , 0
K, 36 ,-1000 , 1000 , 0
K, 37 , 1000 , 1000 , 0
K, 38 , 3000 , 1000 , 0
K, 39 , 5000 , 1000 , 0
K, 40 ,-4000 , 2000 ,-800
K, 41 ,-2000 , 2000 ,-800
K, 42 , 0 , 2000 ,-800
K, 43 , 2000 , 2000 ,-800
K, 44 , 4000 , 2000 ,-800
K, 45 ,-5000 , 3000 , 0
K, 46 ,-3000 , 3000 , 0
K, 47 ,-1000 , 3000 , 0
K, 48 , 1000 , 3000 , 0
K, 49 , 3000 , 3000 , 0
K, 50 , 5000 , 3000 , 0
K, 51 ,-4000 , 4000 ,-800
K, 52 ,-2000 , 4000 ,-800
K, 53 , 0 , 4000 ,-800
K, 54 , 2000 , 4000 ,-800
K, 55 , 4000 , 4000 ,-800
K, 56 ,-5000 , 5000 , 0
K, 57 ,-3000 , 5000 , 0
K, 58 ,-1000 , 5000 , 0
K, 59 , 1000 , 5000 , 0
K, 60 , 3000 , 5000 , 0
K, 61, 5000 , 5000 , 0
*DO,I,1,5 !利用循环创建多段下弦杆直线
L, I , I+1
*ENDDO
*DO,I,12,16
L, I , I+1
*ENDDO
*DO,I,23,27
L, I , I+1
*ENDDO
*DO,I,34,38
L, I , I+1
*ENDDO
*DO,I,45,49
L, I , I+1
*ENDDO
*DO,I,56,60
L, I , I+1
*ENDDO
*DO,I,1,6
L, I , I+11
*ENDDO
*DO,I,12,17
L, I , I+11
*ENDDO
*DO,I,23,28
L, I , I+11
*ENDDO
*DO,I,34,39
L, I , I+11
*ENDDO
*DO,I,45,50
L, I , I+11
*ENDDO !利用循环创建多段下弦杆直线
*DO,I,7,11 !利用循环创建多段斜杆直线
L, I , I-6
*ENDDO
*DO,I,18,22
L, I , I-6
*ENDDO
*DO,I,29,33
L, I , I-6
*ENDDO
*DO,I,40,44
L, I , I-6
*ENDDO
*DO,I,51,55
L, I , I-6
*ENDDO
*DO,I,7,11
L, I , I-5
*ENDDO
*DO,I,18,22
L, I , I-5
*ENDDO
*DO,I,29,33
L, I , I-5
*ENDDO
*DO,I,40,44
L, I , I-5
*ENDDO
*DO,I,51,55
L, I , I-5
*ENDDO
*DO,I,7,11
L, I , I+5
*ENDDO
*DO,I,18,22
L, I , I+5
*ENDDO
*DO,I,29,33
L, I , I+5
*ENDDO
*DO,I,40,44
L, I , I+5
*ENDDO
*DO,I,51,55
L, I , I+5
*ENDDO
*DO,I,7,11
L, I , I+6
*ENDDO
*DO,I,18,22
L, I , I+6
*ENDDO
*DO,I,29,33
L, I , I+6
*ENDDO
*DO,I,40,44
L, I , I+6
*ENDDO
*DO,I,51,55
L, I , I+6
*ENDDO !利用循环创建多段斜杆直线
*DO,I,7,10 !利用循环创建多段上弦杆直线
L, I , I+1
*ENDDO
*DO,I,18,21
L, I , I+1
*ENDDO
*DO,I,29,32
L, I , I+1
*ENDDO
*DO,I,40,43
L, I , I+1
*ENDDO
*DO,I,51,54
L, I , I+1
*ENDDO
*DO,I,7,11
L, I , I+11
*ENDDO
*DO,I,18,22
L, I , I+11
*ENDDO
*DO,I,29,33
L, I , I+11
*ENDDO
*DO,I,40,44
L, I , I+11
*ENDDO !利用循环创建多段上弦杆直线
LSEL,S,,,1,60,1, !选择上弦杆指定单元属性
LATT,1,2,1
LSEL,S,,,61,160,1, !选择斜杆指定单元属性
LATT,1,1,1
LSEL,S,,,161,200,1, !选择下弦杆指定单元属性
LATT,1,3,1
LSEL,S,,,61,136,25, !选择斜杆指定单元属性
LATT,1,6,1
LSEL,S,,,71,146,25, !选择斜杆指定单元属性
LATT,1,6,1
LSEL,S,,,81,156,25, !选择斜杆指定单元属性
LATT,1,6,1
LSEL,S,,,63,138,25, !选择斜杆指定单元属性
LATT,1,4,1
LSEL,S,,,73,148,25, !选择斜杆指定单元属性
LATT,1,4,1
LSEL,S,,,83,158,25, !选择斜杆指定单元属性
LATT,1,4,1
LSEL,S,,,1,26,5, !选择上弦杆指定单元属性
LATT,1,5,1
LSEL,S,,,31,55,6, !选择上弦杆指定单元属性
LATT,1,5,1
LSEL,S,,,32,56,6, !选择上弦杆指定单元属性
LATT,1,5,1
LSEL,S,,,2,27,5, !选择上弦杆指定单元属性
LATT,1,5,1
ALLSEL,all !网格划分
LESIZE,all,,,1,,,,,1
LMESH,all
!接下来可以利用命令定义边界条件
D,37,all,0 !完全固定结点37,39,41,47,49,51,57,59,61
D,39,all,0
D,41,all,0
D,47,all,0
D,49,all,0
D,51,all,0
D,57,all,0
D,59,all,0
D,61,all,0
ALLSEL$FINISH !退出前处理模块
/SOLU$ACEL,9800,0,0$SOLVE !自重为第一荷载步
F,1,FZ,-140000
F,2,FZ,-242500
F,3,FZ,-150000
F,4,FZ,-75000
F,5,FZ,-19000
F,6,FZ,0
F,7,FZ,-280000
F,8,FZ,-485000
F,9,FZ,-300000
F,10,FZ,-150000
F,11,FZ,-38000
F,12,FZ,0
F,13,FZ,-280000
F,14,FZ,-485000
F,15,FZ,-300000
F,16,FZ,-150000
F,17,FZ,-38000
F,18,FZ,0
F,19,FZ,-280000
F,20,FZ,-485000
F,21,FZ,-300000
F,22,FZ,-150000
F,23,FZ,-38000
F,24,FZ,0
F,25,FZ,-280000
F,26,FZ,-485000
F,27,FZ,-300000
F,28,FZ,-150000
F,29,FZ,-38000
F,30,FZ,0
F,31,FZ,-140000
F,32,FZ,-242500
F,33,FZ,-150000
F,34,FZ,-75000
F,35,FZ,-19000
F,36,FZ,0
SOLVE !节点荷载为第二荷载步
FINISH !结束求解状态
/post1
etable,evolume,volu
ssum
*get,vtot,ssum,,item,evolume
etable,smax_e,ls,1
esort,etab,smax_e,0,1,,
*get,smaxe,sort,,max
*get,dzmax,node,53,u,z
dzmax=abs(dzmax)
finish
*end
*use,opgrid
/opt
opanl,opgrid,
opvar,t1,dv,1,5
opvar,t3,dv,3,5
opvar,t4,dv,5,10
opvar,t5,dv,4,8
opvar,t6,dv,7,12
opvar,dzmax,sv,5,30,0.1
opvar,smaxe,sv,180,210,1
opvar,vtot,obj,,,
optype,first
opfrst,25
opsave,,,
opexe
*status
finish
|