跪请高手帮助!!!

我是一个即将毕业的研究生,但是由于自己的水平和语言问题(我是俄语生),所以学习ANSYS非常吃力,但是还是坚持学习了命令流的方法,觉得受益非浅,但是,目前在做到接触的时候总是不成功,希望能得到大家的帮助,小弟在这里先谢过了 !!!!
问题主要有:
1 求解的过程中程序自动退出,没有任何提示信息
2 警告稀疏矩阵求解失败
3 我用的是solid 45,但是划分的时候提示要我用 10 节点的单元,不过看书上的例子是用的这样的单元!!!!
4 有的时候能计算到最后,但是只能看到位移的状况,看不了应力,好象没有!!!!
5 接触对的设置中,关于接触刚度的大小,最初渗透的大小的把握上不知道!!!!!!!!

真心的希望能得到大家的帮助,谢谢!!!!!!!
命令流见下面:
/filname,zhuangtu
/title，xuebao
/units,si !国际单位制
fe=0.2 !摩擦系数
fk=4 !法向接触刚度
fd=3 !切向接触刚度
ftoln=0.1 !初始渗透因子

/PREP7
et,1,plane42
ET,2,SOLID45
ET,3,TARGE170
ET,4,CONTA173

KEYOPT,4,12,2 !!!设置摩擦类型，法向不分开，但是允许滑动

MP,EX,1,2.5E+010 !桩的参数
MP,PRXY,1,0.2
MP,DENS,1,2500

MP,EX,2,3.8E+006 ！第一层土的参数
MP,PRXY,2,0.3
TB,DP,2,,,
TBMODIF,1,1,15000
TBMODIF,1,2,21.5
TBMODIF,1,3,
MP,DENS,2,1857

MP,EX,3,2.5E+006 ！第二层土的参数
MP,PRXY,3,0.4
TB,DP,3,,,
TBMODIF,1,1,16800
TBMODIF,1,2,19.2
TBMODIF,1,3,
MP,DENS,3,1690

MP,DENS,4,1890
MP,EX,4,1.0E+010
MP,PRXY,4,0.23
TB,DP,4,,,
TBMODIF,1,1,0
TBMODIF,1,2,30
TBMODIF,1,3,

MP,EX,5,8.8E+006
MP,PRXY,5,0.3
TB,DP,5,,,
TBMODIF,1,1,25000
TBMODIF,1,2,23
TBMODIF,1,3,
MP,DENS,5,1800

MP,EX,6,1.2E+007 ！第五层土的参数
MP,PRXY,6,0.21
TB,DP,6,,,
TBMODIF,1,1,0
TBMODIF,1,2,35
TBMODIF,1,3,
MP,DENS,6,1900

!!!!!!!!!!!!!!!!!!!!!!建模

block,0,2,0,-0.4,0,-2
block,0,2,-0.4,-1.6,0,-2
block,0,2,-1.6,-5.6,0,-2
block,0,2,-5.6,-10.6,0,-2
block,0,2,-10.6,-15,0,-2
block,0,2,-15,-30,0,-2
wpro,,90.000000,
CYLIND, ,0.35,0,30,270,360,
allsel
FLST,2,6,6,ORDE,2
FITEM,2,1
FITEM,2,-6
VSBV,P51X,7
! vplot
CYLIND, 0.2,0.35,0.4,1.6,270,360,
CYLIND, 0.2,0.35,1.6,5.6,270,360,
CYLIND, 0.2,0.35,5.6,10.6,270,360,
CYLIND, 0.2,0.35,10.6,15,270,360,
CYLIND, 0.2,0.35,15,30,270,360,
allsel
vglue,all
nummrg,all
numcmp,all
CYLIND, ,0.35,0,0.4,270,360,
CYLIND, ,0.2,0.4,1.6,270,360,
CYLIND, ,0.2,1.6,5.6,270,360,
CYLIND, ,0.2,5.6,10.6,270,360,
CYLIND, ,0.2,10.6,15,270,360,
CYLIND, ,0.2,15,30,270,360,
wpro,,-90.000000,

WPSTYLE,,,,,,,,0

vadd,12,13,14,15,16
vglue,17,5
NUMCMP,LINE
NUMCMP,VOLU
！NUMCMP,KP

!!!!!!!!!!!!!!!!!!!!!!!!!土体网格的划分

lsel,s,,,40,67,27 ！！！！vsel,s,,, aslv,all, lsla,all lplot
lsel,a,,,39,65,26
lsel,s,,,66,68,2
lsel,a,,,1
lsel,a,,,3,6,3
lsel,s,,,9,38,29
allsel
accat,1,6
lccat,2,11
lccat,10,4
lsel,s,,,10,11,1
lsel,a,,,2,4,2
lesize,all,,,5
lsel,s,,,40,67,27
lsel,a,,,39,65,26
lesize,all,,,10,0.3
lsel,s,,,66,68,2
lsel,a,,,1
lsel,a,,,3,6,3
lesize,all,,,4
lsel,s,,,9,38,29
lesize,all,,,10
type,2
mat,2
vmesh,6

allsel
accat,34,35
lccat,10,4
lccat,13,18
lsel,s,,,10,18,8
lsel,a,,,4,13,9
! lplot
lesize,all,,,5
lsel,s,,,67,72,5
lsel,a,,,65,73,8
! lplot
lesize,all,,,10,0.3
lsel,s,,,17,19,2
lsel,a,,,69
lsel,a,,,70,71,1
lesize,all,,,12
lsel,s,,,9,12,3
lesize,all,,,10
type,2
mat,2
vmesh,7

allsel
accat,40,39
lccat,20,23
lsel,s,,,20,23,3
lesize,all,,,5
lsel,s,,,77,78,1
lesize,all,,,10,0.3
lsel,s,,,49,50,1
lsel,a,,,74,76,1
lesize,all,,,40
lsel,s,,,25
lesize,all,,,10
type,2
mat,3
vmesh,8

allsel
accat,45,44
lccat,28,24
lsel,s,,,24,28,4
lesize,all,,,5
lsel,s,,,82,83,1
lesize,all,,,10,0.3
lsel,s,,,53,54,1
lsel,a,,,80,81,1
lsel,a,,,79
lesize,all,,,50
lsel,s,,,33
lesize,all,,,10
type,2
mat,4
vmesh,9

allsel
accat,49,50
lccat,30,31
lsel,s,,,30,31,1
lesize,all,,,5
lsel,a,,,87,88,1
lesize,all,,,10,0.3
lsel,s,,,57,58,1
lsel,a,,,84,86,1,
lesize,all,,,44
lsel,s,,,41
lesize,all,,,10
type,2
mat,5
vmesh,10

allsel
accat,54,55
lccat,32,34
lsel,s,,,32,34,2
lesize,all,,,5
lsel,a,,,92,93,1
lesize,all,,,10,0.3
lsel,s,,,89,91,1
lsel,a,,,61,62,1,
lesize,all,,,150
lsel,s,,,45
lesize,all,,,10
type,2
mat,6
vmesh,11

allsel
lsel,s,,,7,9,2
lsel,a,,,12,15,3
lesize,all,,,10
lsel,s,,,21,22,1
lsel,a,,,17,19,2
lesize,all,,,12
lsel,s,,,14,16,2
lsel,a,,,5,8,3
lesize,all,,,3
type,2
mat,2
vmesh,1

allsel
lsel,s,,,25,27,2
lesize,all,,,10
lsel,s,,,49,52,1
lesize,all,,,40
lsel,s,,,26,29,3
lesize,all,,,3
type,2
mat,3
vmesh,2
allsel
lsel,s,,,33,36,3
lesize,all,,,10
lsel,s,,,55,56,1
lesize,all,,,50
lsel,s,,,35,37,2
lesize,all,,,3
type,2
mat,4
vmesh,3
allsel
lsel,s,,,41,43,2
lesize,all,,,10
lsel,s,,,59,60,1
lesize,all,,,44
lsel,s,,,42,44,2
lesize,all,,,3
type,2
mat,5
vmesh,4

allsel
lsel,s,,,45,47,2
lesize,all,,,10
lsel,s,,,46,48,2
lesize,all,,,3
type,2
mat,6
vmesh,5

TYPE,2
MAT, 6
ESIZE,0.08
VSEL,s,,,12
VMESH,12

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下是桩的网格划分

TYPE,2
MAT,1
ESIZE,0.05, ,
MSHAPE,1,3D
MSHKEY,0 !!!! 用的是自由划分
VSEL,s,,,13
vplot
VMESH,13

!!!!!!!!!!下面是用鼠标操作的去除网格划分时连接的面和线

/PREP7
CM,_Y,AREA
ASEL,R,ACCA
*SET,_Z1,ARINQR(0,13)
*IF,_z1,ne,0,then
ADELE,ALL
*SET,_Z2,ARINQR(0,13)
*SET,_Z3,_Z1-_Z2
*IF,_Z3,NE,0,THEN
CMSEL,S,_Y
CMDELE,_Y
*ENDIF
*ELSE
CMSEL,S,_Y
CMDELE,_Y
*ENDIF
!*
CM,_Y,LINE
*SET,_Z1,LSINQR(0,13)
*IF,_z1,ne,0,then
LSEL,R,LCCA
*SET,_Z2,LSINQR(0,13)
LDELE,ALL
*SET,_Z3,_Z1-_Z2
*IF,_Z3,NE,0,THEN
CMSEL,S,_Y
CMDELE,_Y
*ENDIF
*ELSE
CMSEL,S,_Y
CMDELE,_Y
*ENDIF


allsel
keyopt,4,9,0
keyopt,4,12,2

mp,mu,1,0.2
mp,mu,2,0.2
mp,mu,3,0.2
mp,mu,4,0.2
mp,mu,5,0.2
mp,mu,6,0.3

!!!!!!做六个接触对，每层土一个，其中第一层土作了两个，桩的底面单独作了一个接触对

R,10, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,10
vsel,s,,,12
vplot
asel,s,,,61,89,28
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
！eplot
allsel
vsel,s,,,1,6,5
vplot
asel,s,,,2,32,30
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,11, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,11
asel,s,,,66
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
allsel
asel,s,,,5
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
eplot
allsel

r,12, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,12
asel,s,,,90
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
allsel
asel,s,,,16
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,13, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,13
asel,s,,,95
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
allsel
asel,s,,,20
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel
r,14, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,14
asel,s,,,98
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
！eplot
allsel
asel,s,,,24
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,15, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,15
asel,s,,,80
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
！eplot
allsel
asel,s,,,60
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel



/solu
asel,s,loc,x,0
aplot
da,all,symm
asel,s,loc,x,2
aplot
da,all,all
asel,s,loc,z,0
aplot
da,all,symm
asel,s,loc,z,-2
da,all,all
asel,s,loc,y,-30
da,all,all
allsel

ACEL,0,9.8,0
csys,5
nsel,s,loc,z,0
nsel,r,loc,x,0.01,0.34
esln,s
eplot
csys,0
sfe,all,3,pres,,1000

nropt,full,,off ！！ 全牛顿-拉普拉斯解选项
antype,static ！！ 分析类型
!nlgeom,on ！！ 打开大变形
time,1
autot,on ! ！ 自动载荷步打开
!nsub,400
pred,on ！！ 预测求解器
outres,all,all ！！ 记录所有的结果文件
allsel
solv
finish