请大家帮我看看，怎么划分才最好？

装配式空心斜板桥，空心板和铰缝用45号单元模拟，采用扫掠网格划分，感觉网格划分得太细，内存不足，而我将网格划分粗略一点，却又不行，很迷茫！请各位不吝指教！命令流如下：
finish
/clear
/prep7
/title External Transverse Prestress for Sressgthen Girder Bridges
/uis,msgpop,3
*SET,qmpz,2*0.8625*3.5*23000
*SET,rxdlg,8000
*SET,g1,(2*rxdlg+qmpz)/9
*SET,jfarea,0.01*0.6
*SET,g2,jfarea*26000
*SET,kxbarea,0.3123
*SET,g3,kxbarea*26000
*SET,g,g1+g2+g3
*SET,B,1
*SET,H,1
*SET,A1,B*H
*SET,I1,B*H**3/12
*SET,I2,H*B**3/12
*SET,L,17.0
*SET,Lo,15.6 !梁的计算跨径
et,1,solid45
et,2,solid45
mp,dens,1,2600
mp,ex,1,3.45e10
mp,prxy,1,0.1667
mp,dens,2,0.1*g/(kxbarea+jfarea)
mp,ex,2,3.45e10
mp,prxy,2,0.1667
k,1,0.99,0,0
k,2,1,0,0
k,3,1,0.6,0
k,4,0.99,0.6,0
a,1,2,3,4
rectng,0,0.99,0,0.6
cyl4,0.27,0.26,0.19,-180
cyl4,0.27,0.34,0.19,180
rectng,0.08,0.46,0.26,0.34
agen,2,3,5,1,0.45,0,0,100
asel,s,,,3,8,1
cm,ar,area
allsel
asba,2,ar
vext,1,9,8,0,0,L
vsel,s,loc,x,0,0.99
vatt,2,2,2
allsel
vsel,s,loc,x,0.99,1
vatt,1,1,1
allsel
wprota,0,0,-30
cswpla,100,0
vsbw,all
asbw,all
lsbw,all
csys,100
vsel,s,loc,y,0,0.6
vsel,r,loc,x,0,1/cos(30*acos(-1)/180)
vsel,r,loc,z,-sin(30*acos(-1)/180),0-0.1
cm,v1,volu
vdele,v1
asel,s,loc,y,0,0.6
asel,r,loc,x,0,1/cos(30*acos(-1)/180)
asel,r,loc,z,-sin(30*acos(-1)/180),0-0.1
cm,a1,area
adele,a1
asel,s,loc,y,0,0.6
asel,r,loc,x,0,1/cos(30*acos(-1)/180)
asel,r,loc,z,-sin(30*acos(-1)/180),0-0.1
cm,a1,area
adele,a1
lsel,s,loc,y,0,0.6
lsel,r,loc,x,0,1/cos(30*acos(-1)/180)
lsel,r,loc,z,-sin(30*acos(-1)/180),0-0.1
cm,l1,line
ldele,l1
ksel,s,loc,y,0,0.6
ksel,r,loc,x,0,1/cos(30*acos(-1)/180)
ksel,r,loc,z,-sin(30*acos(-1)/180),0-0.1
cm,k1,kp
kdele,k1
allsel
wpcsys,-1,0
csys,0
wpoff,0,0,Lo
wprota,0,0,-30
cswpla,101,0
vsbw,all
asbw,all
lsbw,all
wpcsys,-1,0
allsel
csys,101
vsel,s,loc,y,0,0.6
vsel,r,loc,x,0,(L-Lo)/cos(60*acos(-1)/180)
vsel,r,loc,z,0.1,(L-Lo)*sin(60*acos(-1)/180)
cm,v1,volu
vdele,v1
asel,s,loc,y,0,0.6
asel,r,loc,x,0,(L-Lo)/cos(60*acos(-1)/180)
asel,r,loc,z,0.1,(L-Lo)*sin(60*acos(-1)/180)
cm,a1,area
adele,a1
asel,s,loc,y,0,0.6
asel,r,loc,x,0,(L-Lo)/cos(60*acos(-1)/180)
asel,r,loc,z,0.1,(L-Lo)*sin(60*acos(-1)/180)
cm,a1,area
adele,a1
lsel,s,loc,y,0,0.6
lsel,r,loc,x,0,(L-Lo)/cos(60*acos(-1)/180)
lsel,r,loc,z,0.1,(L-Lo)*sin(60*acos(-1)/180)
cm,l1,line
ldele,l1
ksel,s,loc,y,0,0.6
ksel,r,loc,x,0,(L-Lo)/cos(60*acos(-1)/180)
ksel,r,loc,z,0.1,(L-Lo)*sin(60*acos(-1)/180)
cm,k1,kp
kdele,k1
wpcsys,-1,0
csys,0
vsel,all
nummrg,all
numcmp,all
wpoffs,0.27,0,0
wprota,0,0,90
vsbw,all
wpoffs,0,0,0.45
vsbw,all
wpoffs,0,0.26,0
wprota,0,90,0
vsbw,all
wpoffs,0,0,-0.08
vsbw,all
wpcsys,-1,0
wpoffs,0.495,0,0
wprota,0,0,90
vsbw,all
vsel,all
vglue,all
wpcsys,-1,0
et,5,mesh200,6
esize,0.2
amap,70,27,56,52,74
amap,10,79,28,27,74
amap,95,28,58,54,79
amap,105,56,30,84,52
amap,51,29,83,84,30
amap,108,29,58,54,83
amap,81,31,66,59,84
amap,72,83,32,31,84
amap,101,32,64,61,83
amap,96,66,34,67,59
amap,19,1,2,68,67
amap,53,33,34,67,76
amap,86,76,75,68,67
amap,78,33,64,61,76
amap,1,76,75,3,4
lsel,s,length,,15.6
cm,l2,line
lesize,l2,,,26
vsel,all
vsweep,all
vsel,all
nummrg,all
numcmp,all
csys,0
vsel,s,loc,x,0,0.99
vsel,r,loc,y,0,0.6
vsel,r,loc,z,0,20.8
cm,vvv,volu
csys,100
vgen,9,vvv,,,1/cos(30*acos(-1)/180),,,1000
csys,0
allsel
vsel,s,loc,x,0.99,1
vsel,r,loc,z,0,20.8
vsel,r,loc,y,0,0.6
cm,vvvv,volu
csys,100
vgen,8,vvvv,,,1/cos(30*acos(-1)/180),,,10000
*dim,zzx,array,18
*dim,zx,array,18
*SET,zzx(1),0.27
*SET,zzx(2),0.72
*SET,zzx(3),1.27
*SET,zzx(4),1.72
*SET,zzx(5),2.27
*SET,zzx(6),2.72
*SET,zzx(7),3.27
*SET,zzx(8),3.72
*SET,zzx(9),4.27
*SET,zzx(10),4.72
*SET,zzx(11),5.27
*SET,zzx(12),5.72
*SET,zzx(13),6.27
*SET,zzx(14),6.72
*SET,zzx(15),7.27
*SET,zzx(16),7.72
*SET,zzx(17),8.27
*SET,zzx(18),8.72
*do,i,1,18
*SET,zx(i),zzx(i)/cos(30*acos(-1)/180)
*enddo
wprota,0,0,-30
cswpla,100,0
wpcsys,-1,0
csys,100
d,node(zx(1),0,0),ux,,,,,,uy,uz
*do,i,2,18
d,node(zx(i),0,0),,,,,,,uy,uz
*enddo
csys,0
wpoffs,0,0,Lo
wprota,0,0,-30
cswpla,101,0
csys,101
d,node(zx(1),0,0),ux,,,,,,uy
*do,i,2,18
d,node(zx(i),0,0),,,,,,,uy
*enddo
/SOL
allsel
antype,static
acel,,9.8
solve
save
/POST1
prrsol,fy

[ 本帖最后由 twg2012 于 2010-11-17 19:32 编辑 ]