运行出错，请教（附命令流和图）
明明我在前面已经定义了材料3，为什么计算前检查说没有材料3呢，真奇怪命令流如下
是个振动计算的
fini
/clear
/title,zhengti
r1=6.85
r2=7.5
/prep7
et,1,plane42
mp,ex,1,3e10
mp,prxy,1,0.25
mp,dens,1,2500
et,2,combin14,,,2!用于土弹簧
mp,ex,2,5000
mp,prxy,2,0.1667
mp,dens,2,1800
et,3,solid45
mp,ex,3,3e10
mp,prxy,3,0.25
mp,dens,3,2500
rectng,-6.55,6.55,-2.71,-2.11
rectng,-r2,-5.32,-r2,-2.7
rectng,-2.06,-1.76,-r2,-2.7
rectng,1.76,2.06,-r2,-2.7
rectng,5.32,r2,-r2,-2.7
rectng,2.0,5.0,-r2,-5.8
rectng,-5.0,-2.0,-r2,-5.8
rectng,-1.8,1.8,-r2,-6.2
rectng,-0.77,-0.67,-6.2,-6.05
rectng,0.67,0.77,-6.2,-6.05
cyl4,,,6.6,288,r1,327
cyl4,,,6.6,216,r1,252
cyl4,,,6.6,252,r1,288
aadd,all
cm,in_area,area
asel,none
cyl4,,,r1,0,15,360
cm,out_c,area
asel,all
asba,in_area,out_c
cm,in_area,area
allsel,all
!asel,none
!rectng,-5.6,5.6,4.62,4.89
!cm,liang,area
!asel,none
!cyl4,,,6.85,0,10,360
!cm,out_c,area
!allsel,all
!asba,liang,out_c
nummrg,all
numcmp,all
asel,all
aplot
*do,i,1,4
cyl4,,,r1,(i-1)*90,r2,i*90
*enddo
aplot
!aadd,all
/pnum,area,1
nummrg,all
numcmp,all
asel,all
aglue,all
nummrg,all
numcmp,all
vext,all,,,,,20
asel,s,loc,z,0
aplo
asel,s,,,1
aplot
lsla,s,1
lplot
mat,1
lesize,all,0.15
nummrg,all
numcmp,all
amesh,all
allsel
asel,s,,,2
aplo
lsla,s,1
/pnum,line,1
lplo
lesize,30,,,90
lesize,32,,,90
lesize,31,,,7
lesize,33,,,7
mshkey,1
amesh,2
allsel
asel,s,,,3
aplo
lsla,s,1
/pnum,line,1
lplo
lesize,34,,,90
lesize,36,,,90
lesize,35,,,7
mshkey,1
amesh,3
allsel
asel,s,,,4
aplo
lsla,s,1
/pnum,line,1
lplo
lesize,39,,,90
lesize,all,0.15
mshkey,0
amesh,4
allsel
asel,s,,,5
aplo
lsla,s,1
/pnum,line,1
lplo
lesize,37,,,90
lesize,all,0.15
mshkey,0
amesh,5
alls
csys,1
asel,s,loc,x,7,8
asel,r,loc,y,-5,5
aplo
lsla,s,1
lplo
lsel,s,loc,x,7.4,7.6
lsel,r,loc,y,-1,1
lplo
lesize,all,,,20
alls
type,3
mat,3
extopt,aclear,20
vsweep,all!体扫掠
!建立土弹簧
csys,1
nsel,s,loc,x,7.45,7.6
nplo
*get,node_min,node,,num,min
*get,node_max,node,,num,max
ngen,2,500000,all,,,0.5
!加弹簧单元
type,2
mat,2
numreal=0
*do,i,node_min,node_max
*if,nsel(i),eq,1,then
numeral=numreal+1
r,numreal,5000000
real,numreal
e,i,i+500000
*endif
*enddo
allsel
eplo
!计算分析
finish
/solu
antype,trans
time,0.1
deltim,0.01
autots,off
timint,on
kbc,1
outres,all,all
allsel
csys,1
nsel,s,loc,x,7.7,8.1
nplo
d,all,ux
d,all,uy
alls
*DIM,train_t,ARRAY,1401,1,1,time,
*DIM,train_t,TABLE,1401,1,1,time,
*TREAD,TRAIN_T,'train_t','txt',' ',2,
csys,0
allsel,all
nsel,s,loc,x,-0.77,-0.67
nsel,r,loc,y,-6.06,-6.04
f,all,fy,%train_t%
alls
nsel,s,loc,x,0.67,0.77
nsel,r,loc,y,-6.06,-6.04
f,all,fy,%train_t%
alls
acel,,-9.8
solve
请高人指点，谢谢
[ 本帖最后由 fisher-net 于 2008-10-28 16:26 编辑 ]
回复 楼主 fisher-net 的帖子
这里有错误type,3
mat,3
extopt,aclear,20
vsweep,all!体扫掠
!建立土弹簧
csys,1
nsel,s,loc,x,7.45,7.6
nplo
*get,node_min,node,,num,min
*get,node_max,node,,num,max
ngen,2,500000,all,,,0.5
!加弹簧单元
type,2
mat,2
没有给ngen生成的节点附加任何材料属性。 非常感谢楼上的回复
原来是用numcmp和nummrg所造成两种数据完全相同的材料变成一种了
页:
[1]