接触分析中，速度和位移的问题，求助高手

spiderliu

我做接触分析
模型是球在平面上滑动，在球上加载，考虑摩擦生热，计算接触面上的应力和温度，用热应力耦合单元计算，摩擦生热部分计算出热流密度，将热流密度加载在接触面上，滑动速度是18.85m/s，计算时间是0.03,（计算出来的热达到稳态需要的时间），加载时，给球一个18.85*0.03的位移，这样就相当于速度是18.85
但是最近计算大载荷时将网格密度画得很细，所以将模型网格太多，所以想减小滑动距离，但有稳态时间的限制，如果滑动距离减小了，相应的速度就小于18.85了
想到一个办法就是分两步加载，第一个载荷步计算时间为0.025，此时不滑动，第二个载荷步计算时间为0.005，滑动距离为0.005*18.85
不知道这样能否实现
求高手指点
附程序段
/SOLU
ANTYPE,TRANS
NLGEOM,ON   
ASEL,S,,,5,6,1
ASEL,A,,,10,11,1
DA,ALL,UX,0
ASEL,S,,,1,8,7
DA,ALL,UZ,0
ALLSEL
DA,3,UY,0
ASEL,S,,,2,9,7
ASEL,A,,,12,16,4
DA,ALL,SYMM
ALLSEL
KBC,1
time,0.045
ASEL,S,,,15
NSLA,S,1
D,ALL,UX,0
allsel
ASEL,S,,,7
NSLA,S,1
NSEL,R,LOC,Z,-20*L,0
SF,ALL,HFLUX,FLUX
ALLSEL  
NSUBST,80,100,80
lswrite,1
time,0.005
ASEL,S,,,15
NSLA,S,1
D,ALL,UX,s
ALLSEL
NSEL,S,LOC,X,RAD
NSEL,R,LOC,Y,RAD+H
NSEL,R,LOC,Z,0
F,ALL,FY,-FORCE
ASEL,S,,,7
NSLA,S,1
NSEL,R,LOC,Z,-20*L,0
SF,ALL,HFLUX,FLUX
ALLSEL
NSUBST,80,120,80
lswrite,2
AUTO,ON
EQSLV,PCG
!NROP,UNSYM  
lsSOLVE,1,2,1   
SAVE