Проблема контактной задачи

Добрый день!

Проблема в следующем:

Между двумя контактными поверхностями имеется зазор, при задании давления на одну из частей конструкции, для уменьшения зазора, и осуществления контакта, контакт не находит себя и пролетает мимо.

Задачу реализую в Mechanical APDL, прикрепляю код программы, задача решается в 3 этапа (вычисление расширения кольца и вала под определенной температурой; вычисление расширения шарика под определенной температурой; относительно вычисленных значений построение модели, прогрев и последним шагом необходимо приложить давление, что не выходит):

finish
/clear
/filname,Kolco_Sharic_Val
/prep7            !переход в препроцессор
/pnum,area,1    !включение нумерации областей
/pnum,line,1    !включение нумерации линий
/pnum,kp,1    !включение нумерации ключевых точек
/psf,pres,,1
/eshape,1

DVal = 95.065e-3
DKolcoVnutr = 94.980e-3
DKolcoVnesh = 112.648e-3
Shirina = 32e-3

RadiusVdavl = 12.3e-3
DVdavl = 108.742e-3
HVdavl = (DKolcoVnesh-DVdavl)/2

DShara = 23.820e-3

P=28000

CYL4,0,0,DVal/2-1e-3,90,
CYL4,0,0,DVal/2-1e-3,,DVal/2,90
CYL4,0,0,DKolcoVnutr/2,,DKolcoVnutr/2+1e-3,90
CYL4,0,0,DKolcoVnutr/2+1e-3,,DKolcoVnesh/2,90
vext,all,,,,,-Shirina/2
k,1000,,DVdavl/2+RadiusVdavl
kwpave,1000
wprota,,,90
CYL4,0,0,RadiusVdavl
vrota,24,,,,,,3,18,90
vsbv,4,5

k,100,,DVdavl/2+DShara/2
kwpave,100
sphere,DShara/2
vsbw,4
alls
wprota,,,90
vsbw,7
alls
vdele,5,8,3
alls
k,101,,DVdavl/2+DShara/5+1e-4
wpcsys,0
kwpave,101
wprota,,90
vsbw,4
alls
k,102,kx(37),
wpcsys,0
kwpave,102
wprota,,,90
vsbw,6
alls

vglue,1,2
alls
vsel,s,,,3,4
vsel,a,,,8
vglue,all
alls
vsel,s,,,5,7,2
vglue,all
alls

!vsel,s,,,5,9,4
!vglue,all
!============================================
!    Задание типов конечных элементов
!============================================
ET,1,SOLID227
KEYOPT,1,1,11
!KEYOPT,1,6,1

ET,2,conta173

ET,3,TARGE170

ET,4,conta174
KEYOPT,4,3,1
KEYOPT,4,5,1 
KEYOPT,4,12,5

ESharKolco = 2.11e11
KLTR2 = 10.08e-6
EVal = 2.11e11
KLTR1 = 11.5e-6
KXX1 = 40
nu = 0.3

MP,EX,1,EVal
MP,PRXY,1,nu
MP,ALPX,1,KLTR1
MP,KXX,1,KXX1

MP,EX,2,ESharKolco
MP,PRXY,2,nu
MP,ALPX,2,KLTR2
MP,KXX,2,KXX1

MP,EX,3,ESharKolco
MP,PRXY,3,nu
MP,ALPX,3,KLTR2
MP,KXX,3,KXX1

VSEL,S,,,1,6,5
VATT,1,,1
ALLS
VSEL,s,,,2
vsel,a,,,9,10
VATT,2,,1
alls
VSEL,u,mat,,1,2
VATT,3,,1
alls
ASEL,s,,,10
AATT,1,,3
alls
ASEL,S,,,17
AATT,2,,2
alls
ASEL,S,,,41,42
AATT,2,,3
alls
ASEL,S,,,21
AATT,3,,2
alls

esize,(DShara/5)
MSHAPE,1,2D 
MSHKEY,0
asel,s,type,,2,3
!vsel,s,,,6
!aslv,u,1
amesh,all
alls
esize,(DShara/5)
vsel,s,,,2,6,4
vsel,a,,,5,9,4
MSHAPE,1,3D 
MSHKEY,0
vmesh,all
alls
esize,(DShara/2)
vsel,u,,,2,6,4
vsel,u,,,5,9,4
vmesh,all
alls
!esize,(DShara/50)

asel,s,,,10,17,7
asel,a,,,41
areverse,all
alls
wpcsys,0

nsel,s,loc,x,0
d,all,ux,0
nsel,s,loc,y,0
d,all,uy,0
nsel,s,loc,z,0
d,all,uz,0
esel,s,type,,2
esel,r,mat,,2
ekill,all
alls
esel,s,mat,,3
ekill,all
/solu
ANTYPE,0
AUTOTS,ON
TIME,1E-5
OUTRES,ALL,ALL
NLGEOM,ON
KBC,0
alls
esel,s,mat,,2
nsle
d,all,temp,400
esel,s,mat,,1
nsle
d,all,temp,20
alls
solve
time,10
KBC,0
esel,s,type,,2
esel,r,mat,,2
ealive,all
alls
d,all,temp,0
alls
solve
time,110
!nsubst,10,10,10
alls
esel,s,mat,,2
nsle
d,all,temp,70
esel,s,mat,,1
nsle
d,all,temp,50
alls
solve
/post1
set,last
esel,s,type,,3
esel,r,mat,,2
nsle,s
nsel,r,loc,z,0
nsel,r,loc,x,0
*GET,NN,NODE,0,NUM,MAX
*get,ymax1,node,NN,u,y
!ymax1=0.707e-4

fini
/prep7
alls
vclear,all
aclear,all
vdele,all
adele,all
ldele,all
kdele,all
/prep7
sphere,DShara/2
vsbw,all
wprota,,,90
vsbw,all
wprota,,90
vsbw,all
vsel,s,loc,x,-10,0
vsel,a,loc,y,-10,0
vsel,a,loc,z,-10,0
vdele,all
alls
type,1
mat,3
esize,(DShara/5)
vmesh,all

nsel,s,loc,x,0
d,all,ux
nsel,s,loc,y,0
d,all,uy
nsel,s,loc,z,0
d,all,uz
/solu
!esel,s,mat,,1,2
!ekill,all
esel,s,mat,,3
nsle
d,all,temp,70
alls
solve
/post1
set,last
plnsol,u,y
*get,ymax2,plnsol,0,max
!ymax2=0.0031
vplo

fini
/prep7
alls
vclear,all
aclear,all
vdele,all
adele,all
ldele,all
kdele,all
csys,0
wpcsys,0

CYL4,0,0,DVal/2-1e-3,
CYL4,0,0,DVal/2-1e-3,,DVal/2
CYL4,0,0,DKolcoVnutr/2,,DKolcoVnutr/2+1e-3
CYL4,0,0,DKolcoVnutr/2+1e-3,,DKolcoVnesh/2
wprota,,90
asbw,all
alls
asel,s,loc,y,0.00001,10
adele,all
lsel,s,loc,y,0.00001,100
ldele,all
ksel,s,loc,y,0.001,100
kdele,all
alls
vext,all,,,,,-Shirina/2
wpcsys,0
k,1000,,-(DVdavl/2+RadiusVdavl)
kwpave,1000
wprota,,,90
CYL4,0,0,RadiusVdavl
k,1001,
k,1002,,,10
vrota,30,,,,,,1001,1002
vsbv,4,5
vsel,s,loc,y,0.0001,100
vdele,all
alls
vsbv,9,8
wpcsys,0
k,1003,,-(DVdavl/2+ymax1+DShara/2+ymax2)
kwpave,1003
sphere,DShara/2
vsbw,5
alls
vsel,s,loc,z,0,10
vdele,all
alls
wpcsys,0
kwpave,1003
wprota,,,90
vsbw,6
wprota,,90
vsel,s,,,5,7,2
vsbw,all
alls
wpcsys,0
kwpave,1003
wprota,,,90
wprota,,-60
vsbw,8
wprota,,120
vsbw,10
vadd,5,6
alls
vadd,7,8
alls
vadd,9,11

cswpla,11,0
wpcsys,11
csys,11
k,1004,,-(DVdavl/2+ymax1+DShara/2+ymax2)
kwpave,1004
wprota,,,90
CYL4,0,0,DVal/2-1e-3,
CYL4,0,0,DVal/2-1e-3,,DVal/2
CYL4,0,0,DKolcoVnutr/2,,DKolcoVnutr/2+1e-3
CYL4,0,0,DKolcoVnutr/2+1e-3,,DKolcoVnesh/2
asel,s,loc,y,-(DVdavl/2+ymax1+DShara/2+ymax2),-(DVdavl/2+ymax1+ymax2)
vext,all,,,Shirina/2

wpcsys,0
kwpave,1004
cswpla,12,0
wpcsys,12
csys,12
k,1000,,-(DVdavl/2+RadiusVdavl)
kwpave,1000
!wprota,,,90
CYL4,0,0,RadiusVdavl
k,1001,
k,1002,10
vrota,95,,,,,,1001,1002
alls
vsbv,11,12
vsbv,16,13
vsbv,11,14
vsbv,12,15
wpcsys,12
wprota,,90
vsbw,all
vsel,s,loc,y,0,-10
vdele,all

csys,11
wpcsys,11
wprota,,-120
cswpla,13,0
wpcsys,13
csys,13

k,1004,,-(DVdavl/2+ymax1+DShara/2+ymax2)
kwpave,1004
wprota,,,90
CYL4,0,0,DVal/2-1e-3,
CYL4,0,0,DVal/2-1e-3,,DVal/2
CYL4,0,0,DKolcoVnutr/2,,DKolcoVnutr/2+1e-3
CYL4,0,0,DKolcoVnutr/2+1e-3,,DKolcoVnesh/2
asel,s,loc,y,-(DVdavl/2+ymax1+DShara/2+ymax2),-(DVdavl/2+ymax1+ymax2)
vext,all,,,Shirina/2

wpcsys,0
kwpave,1004
cswpla,14,0
wpcsys,14
csys,14
k,1000,,-(DVdavl/2+RadiusVdavl)
kwpave,1000
!wprota,,,90
CYL4,0,0,RadiusVdavl
k,1001,
k,1002,10
vrota,150,,,,,,1001,1002
alls
vsbv,11,13
vsbv,20,15
vsbv,11,17
vsbv,13,18
wpcsys,14
wprota,,90
vsbw,all
vsel,s,loc,y,0,-10
vdele,all
alls
vglue,1,2
vglue,3,4
vsel,s,,,5,6
vsel,a,,,10
vglue,all
alls
vglue,12,14
vglue,16,19
vglue,13,17
vglue,20,23
vsel,s,,,1,7,6
vsel,a,,,4,12,8
vsel,a,,,9,13,4
VATT,1,,1
vsel,s,,,5,6
vsel,a,,,10
VATT,3,,1
alls
vsel,u,mat,,1,3,2
VATT,2,,1
alls
ASEL,s,,,7,9,2
asel,a,,,69,70
asel,a,,,124,125
AATT,1,,3

ASEL,s,,,20,21
asel,a,,,84,85
asel,a,,,139,140
AATT,2,,2

asel,s,,,12,52,40
asel,a,,,89,98,9
asel,a,,,144,153,9
AATT,2,,3

asel,s,,,48,49
asel,a,,,46,59,13
asel,a,,,30,62,32
AATT,3,,4
alls

esize,(DShara/5)
MSHAPE,1,2D 
MSHKEY,0
asel,s,type,,2,4
!vsel,s,,,6
!aslv,u,1
amesh,all
alls
esize,(DShara/5)
vsel,s,,,4,7,3
vsel,a,,,9
vsel,a,mat,,2
vsel,a,mat,,3
MSHAPE,1,3D 
MSHKEY,0
vmesh,all
alls
esize,(DShara/2)
vsel,u,,,4,7,3
vsel,u,,,9
vsel,u,mat,,2
vsel,u,mat,,3
vmesh,all

asel,s,,,20,21
asel,a,,,7,9,2
asel,a,,,89,98,9
asel,a,,,144,153
areverse,all

csys,0
wpcsys,0
esel,s,mat,,3
ekill,all

csys,12
wpcsys,12
wprota,,,-90
cswpla,12,0
wpcsys,12
csys,12

csys,14
wpcsys,14
wprota,,,-90
cswpla,14,0
wpcsys,14
csys,14

csys,0
wpcsys,0
nsel,s,loc,z,0
d,all,uz

nsel,s,loc,y,0
d,all,uy

vsel,s,,,1,7,6
vsel,a,,,2,3
nslv,s,1
nsel,r,loc,x,0
d,all,ux

csys,12
wpcsys,12
nsel,s,loc,y,0
d,all,uy
vsel,s,,,4,12,8
vsel,a,,,8,16,8
nslv,s,1
nsel,r,loc,x,0
d,all,ux
vsel,s,,,4,12,8
vsel,a,,,8,16,8
nslv,s,1
nrotat,all

csys,14
wpcsys,14
nsel,s,loc,y,0
d,all,uy
vsel,s,,,9,13,4
vsel,a,,,11,20,9
nslv,s,1
nsel,r,loc,x,0
d,all,ux
vsel,s,,,9,13,4
vsel,a,,,11,20,9
nslv,s,1
nrotat,all

nsel,s,loc,y,(DVdavl/2+ymax1+DShara/2+ymax2)
nsel,r,loc,z,0
nsel,r,loc,x,0
d,all,all

!Сделать закрепление!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/solu
ANTYPE,0
AUTOTS,ON
TIME,1E-5
NLGEOM,ON
LNSRCH,ON
OUTRES,ALL,ALL
!KBC,0
esel,s,type,,3
ekill,all
esel,s,mat,,3
ekill,all
esel,s,mat,,2
nsle
d,all,temp,400
alls
solve

time,10
esel,s,type,,3
ealive,all
d,all,temp,0
alls
solve

time,55
esel,s,mat,,3
ealive,all
esel,s,mat,,1
nsle
d,all,temp,50
esel,s,mat,,2,3
nsle
d,all,temp,50
vsel,s,mat,,2
nslv,s,1
ddele,all,ux

alls
solve

time,110
esel,s,mat,,1
nsle
d,all,temp,50
esel,s,mat,,2,3
nsle
d,all,temp,70
alls
solve

time,111
!CNCHECK, AUTO
csys,0
wpcsys,0
nsel,s,loc,y,0
ddele,all,uy
vsel,s,mat,,1
aslv,s,1
asel,r,loc,y,0
nsla,s,1
sfa,all,1,pres,19

vsel,s,mat,,2
aslv,s,1
asel,r,loc,y,0
nsla,s,1
sfa,all,1,pres,19

csys,12
wpcsys,12
nsel,s,loc,y,0
ddele,all,uy
vsel,s,mat,,1
aslv,s,1
asel,r,loc,y,0
nsla,s,1
sfa,all,1,pres,19

vsel,s,mat,,2
aslv,s,1
asel,r,loc,y,0
nsla,s,1
sfa,all,1,pres,19

csys,14
wpcsys,14
nsel,s,loc,y,0
ddele,all,uy
vsel,s,mat,,1
aslv,s,1
asel,r,loc,y,0
nsla,s,1
sfa,all,1,pres,19

vsel,s,mat,,2
aslv,s,1
asel,r,loc,y,0
nsla,s,1
sfa,all,1,pres,19

csys,0
wpcsys,0
alls
solve

/post1
set,last
plnsol,u,sum