帮我看看这个命令流哪里出错了？该如何改？

tiexie1986

二维轴对称模型，模拟管子胀接到管板上，管子管板初始间隙0.2mm，加载分2步：一是180Mpa胀接压力下，管子胀至管板；二是卸掉压力，管子管板产生残余接触压力。



finish
/clear

/filname,2d  
/title,axisymmetric
/prep7  
et,1,plane42   
et,2,targe169   
et,3,conta172   
keyopt,1,3,1        !plane设置为轴对称单元
keyopt,3,5,3        !contact设置icont adjustment为gap/penetration
keyopt,3,9,0        !包含初始化穿透
keyppt,3,10,1
r,3
real,3
R,3,,,0.3,0.15,-0.2,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0

mp,ex,1,1.9e5   
mp,prxy,1,0.3  
tb,miso,1           !设置multilinear，输入应力应变曲线、摩擦系数
tbpt,,0.001,190
tbpt,,0.008,290
tbpt,,0.027,380
tbpt,,0.054,430
tbpt,,0.080,480
tbpt,,0.107,510
tbpt,,0.134,520
tbpt,,0.157,550
mp,mu,1,0.1         
mp,ex,2,2.1e5  
mp,prxy,2,0.3  
tb,miso,2           
tbpt,,0.0023,483
tbpt,,0.013,485
tbpt,,0.033,570
tbpt,,0.053,640
tbpt,,0.076,680
tbpt,,0.116,720
tbpt,,0.155,730
mp,mu,2,0.1

rectng,10.5,12.5,0,85
rectng,12.7,21.7,2,52

lsel,s,,,2,4,2      
lesize,all,,,100
lsel,s,,,1,3,2      
lesize,all,,,4
lsel,s,,,5,7,2      
lesize,all,,,6
lsel,s,,,6,8,2      
lesize,all,,,50
mat,1               
amesh,1              
mat,2
amesh,2
allsel
lsel,s,,,8         
cm,_target,line
type,2              
nsll,s,1            
esurf,all            
allsel
lsel,s,,,2
cm,_contact,line
type,3
nsll,s,1
esurf,all
allsel
finish

/solu
antype,static
NROPT,full
nlgeom,on
autots,on
pred,on
lnsrch,on
outres,all,all
nsub,25,100,1
lsel,s,,,1,5,4
nsll,s,1
d,all,uy

lsel,s,,,4
nsll,s,1
nsel,r,loc,y,4,50
sf,all,pres,180
time,10
allsel
lswrite,1

lsel,s,,,4  
nsll,s,1
nsel,r,loc,y,4,50
sf,all,pres,0
time,20
allsel
lswrite,2
lssolve,1,2,1

[ 本帖最后由 tiexie1986 于 2009-9-18 09:24 编辑 ]

tiexie1986

没人解答啊 。。。。