急请教接触方面的问题

masliu

问题和参数描述如下

就是一个块体在一段梁上的移动
长度单位均为米
条形体
长50，截面 长2.51188，宽1.09788
杨氏模量 2.1e11
泊松比 0.3
密度2000
摩擦因子0.25
移动块体
长1，截面长0.627925，宽0.27447
杨氏模量2.1e9
泊松比0.3
密度7800  

梁为两车道
块体在其中一个车道的中间运行
2秒钟运行完

在调试程序的过程中出现如下错误
*** WARNING ***                         CP =      10.094   TIME= 16:38:01
Overconstraint may occur for Lagrange multiplier or MPC based contact   
algorithm.                                                              
*** WARNING ***                         CP =      59.062   TIME= 16:39:00
There are 4 small equation solver pivot terms.                          
*** WARNING ***                         CP =      87.625   TIME= 16:39:36
There are 4 small equation solver pivot terms.                          
*** WARNING ***                         CP =     150.625   TIME= 16:41:03
Contact element 13624 (real ID 1) has too much penetration/gap         
2.386212893E-04 contact (with target element 13611).                    
*** WARNING ***                         CP =     211.250   TIME= 16:42:36
Contact element 13624 (real ID 1) has too much penetration/gap         
6.603508813E-04 contact (with target element 13611).                    
*** WARNING ***   SUPPRESSED MESSAGE    CP =     212.328   TIME= 16:42:54
Contact element 13624 (real ID 1) has too much penetration/gap         
2.590607259E-04 contact (with target element 13611).                    
*** WARNING ***   SUPPRESSED MESSAGE    CP =     487.938   TIME= 16:49:04
Contact element 13624 (real ID 1) has too much penetration/gap         
1.885876052E-04 contact (with target element 13611).                    
*** ERROR ***     SUPPRESSED MESSAGE    CP =    1343.734   TIME= 17:05:17
Solution not converged at time 1.6E-02 (load step 1 substep 3).         
  Run terminated.                                                        
*** WARNING ***   SUPPRESSED MESSAGE    CP =    1343.734   TIME= 17:05:17
The unconverged solution (identified as time 2 substep 999999) is      
output for analysis debug purposes.  Results should not be used for     
any other purpose.                                                      

拜读了AaronSpark一篇文章后
计算了该模型的模态，居然都为零
可是不知道如何修改
还希望各位高手倾囊相助
程序如下
/title,Lai Wang   !指定标题
*set,g,9.8  
*set,M,50000
*set,m,2000
/prep7                              !进入前处理器
et,1,solid45                        !定义实体单元类型
et,2,targe170                       !定义单元用于三维刚性目标面
et,3,conta174                       !定义接触面
keyopt,3,4,2                              
keyopt,3,2,3
keyopt,3,5,3
keyopt,3,7,1
keyopt,3,9,4
keyopt,3,10,
keyopt,3,12,5   !设置关键字
R,1,,,0.1,0.001,,   
RMORE,,,, ,,
RMORE,,,, ,,
RMORE,,,, ,,
RMORE,        !设置实常数FKT,FKN,FTOLN
mptemp,,,,,,,,                      !设置材料属性
mptemp,1,0
mp,ex,1,2.1e11                      !弹性模量   
mp,prxy,1,0.3                       !泊松比
mp,dens,1,m                         !密度
mp,mu,1,0.25                        !摩擦因数
mp,ex,2,2.1e9                       !弹性模量
mp,prxy,2,0.3                       !泊松比
mp,dens,2,7800                     !密度   
block,,50,,1.09788,,2.51188,        !生成梁
block,,1,1.09788,1.37235,1.569925,2.197895,         !生成方块
CM,_Y,VOLU  
VSEL, , , ,       2
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,14  
FITEM,5,16  
FITEM,5,-17
FITEM,5,19  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,2, , , , ,1  
!*  
FLST,5,4,4,ORDE,2   
FITEM,5,21  
FITEM,5,-24
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,4, , , , ,1  
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,13  
FITEM,5,15  
FITEM,5,18  
FITEM,5,20  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,2, , , , ,1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,2   
FITEM,5,4   
FITEM,5,-5  
FITEM,5,7   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,100, , , , ,1
!*  
FLST,5,4,4,ORDE,2   
FITEM,5,9   
FITEM,5,-12
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,16, , , , ,1
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,1   
FITEM,5,3   
FITEM,5,6   
FITEM,5,8   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,8, , , , ,1  
!*  
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*
type,2                              !改变单元类型
vsel,s,,,1                          !选择编号为1的体
mat,1                               !   
real,1
nslv,s,1                            !选择体上的所有节点
nsel,r,loc,y,1.09788                !在所选节点中再选择y=1.09788
nsel,r,loc,z,1.25594,2.51188        !在所选节点中再选择1.25594<z<2.51188
esurf,,top,                         !根据所选的节点生成刚性接触单元面
allsel,all                          !选择所有实体
vsel,s,,,2                          !选择编号为2的体
nslv,s,1                            !选择体2中所有的节点
nsel,r,loc,y,1.09788                !在所选节点中再选择满足y=1.09788的节点
type,3                              !定义第三种单元类型
mat,2                               !指定材料编号
real,1                              !指定实常数编号
tshap,quad                          !对于接触单元指定几何形状
esurf,,top,                         !在柔性体上生成接触单元
allsel,all                          !选择所有实体
finish
/sol
solc,on
pred,on                             !打开预测，可以增进收敛
CNVTOL,F,,0.05,2,0.01                    !定义收敛条件,力收敛，精度5％
vsel,s,,,1
acel,,9.8,
asel,r,loc,x,0                      !选择x=0的面
asel,a,loc,x,50                     !选择x=50的面
asel,r,loc,y,0,1.09788
nsla,s,1  
d,all,,,,,,all,,,,,                 !在所选节点上施加约束
allsel,all
vsel,s,,,2
acel,,9.8,
asel,r,loc,y,1.09788,1.37235        !选择介于1.09788~1.37235之间的面
asel,r,loc,x,0                      !在所选面中再选择满足x=0的面
nsla,s,1                            !选择所选面上的所有节点
d,all,,50,,,,ux,,,,,                 !在所选面节点上施加ux=50的
allsel,all

nlgeom,1  
nropt,full,,off                     !将牛顿-拉普辛森选项设置full，关闭自下降因子
neqit,25                            !设置合理的平衡迭代次数
Pred,on,,on
ncnv,2   
time,2                              
autots,1                           !   
nsubst,100,1000,10,               !   
outres,all,all,                     !   
solve

masliu

各位大侠，命令流调试以后，出现了两个错误
*** WARNING ***   SUPPRESSED MESSAGE    CP =    2402.281   TIME= 19:43:09
There are 1 small equation solver pivot terms.                          

*** WARNING ***   SUPPRESSED MESSAGE    CP =    2410.391   TIME= 19:43:17
The calculated reference FORCE CONVERGENCE VALUE = 1.654992899E-08 is   
less than 1E-2 a threshold zero, A value of 1E-2 or specified MINREF   
is used.  Check results carefully.   

另外不知道怎么回事
我定义的是d,all,,5,,,,ux,,,,,  
但是算完以后发现质量块先后退同样的距离再前进5m
还有就是什么力都不加载，包括摩擦因子也不用
但是质量块不是沿着直线前进
慢慢开始偏移

希望高手能够指点一二