这段命令的时间步长减少问题

byb

这段命令中，随着计算的进行，time step不断减少，导致了估算时间不断增加。那位大侠可以帮忙提出好的办法解决这个问题呢？希望与大家一起讨论进步，谢谢。

C************************************************************
C*** PARAMETERS
C************************************************************
width=971.3
length=1000
vel_p=-18800
h=10.5
degree=48
x1=710.33/2
r=150.5
upr=90
downr=82
vtime0=0
vtime1=2
vuprotate0=-vel_p/upr
vuprotate1=vuprotate0
vdownrotate0=vel_p/downr
vdownrotate1=vdownrotate0
C************************************************************
C*** ATTRIBUTES
C************************************************************
/PREP7  
ET,1,SHELL163   
ET,2,SHELL163   
ET,3,SHELL163
R, 1, 5/6, 4, h
R,2,6/6, 4,0.01
R,3,6/6, 4,0.01
keyopt,1,1,10
keyopt,2,1,10
keyopt,3,1,10
MP,DENS,1,7.8e-6
MP,EX,1,2.06e8  
MP,NUXY,1,0.3   
TB,BKIN,1,,,,   
TBDAT,1,2.45e5  
TBDAT,2,5.8e5   
EDMP,RIGI,2,7,5
MP,DENS,2,7.8e-6
MP,EX,2,2.06e8  
MP,NUXY,2,0.3   
EDMP,RIGI,2,7,5
MPCOPY, ,2,3
TBCOPY,ALL,2,3
finish
C************************************************************
C*** MODEL WORKPEICE (PLATE)
C************************************************************
/prep7
*afun,deg
asel,none
wprota,,-90,0
rectng,0,width/2,-r*(1-cos(degree))-downr-20,-r*(1-cos(degree))-downr-20-length
wprota,,90,0
*get,aplate,area,,num,min
ESIZE,8
TYPE,   1   
MAT,       1
REAL,       1   
MSHAPE,0,2D
MSHKEY,1
AMESH,1
finish
C************************************************************
C*** MODEL  ROLLER
C************************************************************
/prep7
*afun,deg
k,,0,h/2
k,,x1,h/2
lstr,kp(0,h/2,0),kp(x1,h/2,0)
lnumber1=_return
k,,0,h/2+upr
k,,x1,h/2+upr
k,,x1,h/2+r
circle,kp(x1,h/2+r,0),r,,kp(x1,h/2,0),degree,1
lnumber2=_return
k,,0,-h/2
k,,x1,-h/2
k,,0,-h/2-downr
k,,x1,-h/2-downr
lstr,kp(0,-h/2,0),kp(x1,-h/2,0)
lnumber3=_return
circle,kp(x1,h/2+r,0),r+h,,kp(x1,-h/2,0),degree,1
lnumber4=_return
nummrg,kp
asel,none
arotat,lnumber1,lnumber2,,,,,kp(0,h/2+upr,0),kp(x1,h/2+upr,0)
*get,anumber1,area,,num,min
*get,anumber2,area,,num,max
asel,none
arotat,lnumber3,lnumber4,,,,,kp(0,-h/2-downr,0),kp(x1,-h/2-downr,0)
*get,anumber3,area,,num,min
*get,anumber4,area,,num,max
allsel,all

TYPE,   2   
MAT,       2
REAL,       2   
MSHAPE,0,2D
MSHKEY,1
ASEL,S,,,anumber1,anumber2,1
AMESH,ALL
allsel,all
TYPE,3
MAT,3
REAL,3
MSHAPE,0,2D
MSHKEY,1
ASEL,S,,,anumber3,anumber4,1
AMESH,ALL
allsel,all
finish

C************************************************************
C*** CONTACT
C************************************************************
/prep7
*afun,deg
EDPART,CREATE  
EDCGEN,ASTS,       1,       2,0.2,0.1,0,0,10, , , , ,0,10000000,0,0
EDCGEN,ASTS,       1,       3,0.2,0.1,0,0,10, , , , ,0,10000000,0,0
FINISH
C************************************************************
C*** BOUNDARY CONDITIONS
C************************************************************
/solu
ASEL,S,,,aplate
NSLA,S,1
CM,slab,NODE
EDVE,VELO,SLAB,0,0,vel_p,0,0,0, , , , , ,
NSEL,R,LOC,X,0
D,ALL, , , , , ,roty,rotz
D,ALL, , , , , ,Ux
ALLSEL,ALL
FINISH
C************************************************************
C*** RILLER ROLLING
C************************************************************
/PREP7
*afun,rad
*dim,vtime,array,2
*dim, vuprotate, array, 2
*dim, vdownrotate, array, 2
vtime(1)=vtime0,vtime1
vuprotate(1)=vuprotate0,vuprotate1
vdownrotate(1)=vdownrotate0,vdownrotate1
EDLOAD,ADD,RBRX,0,       2,VTIME,VUPROTATE, 0, , , , ,
EDLOAD,ADD,RBRX,0,       3,VTIME,VDOWNROTATE, 0, , , , ,   
FINISH
C************************************************************
C*** SOLVE
C************************************************************
/solu
time,0.1
EDINT,4,4
EDENERGY,1,1,1,1
EDOPT,ADD,blank,ansys
EDRST,200
EDHTIME,1000
EDDUMP,1,   
EDOUT,RCFORC
solve
finish

byb

高手贵人赶快出手哦。