与Hertz接触的计算结果不符合？

nickler

按道理最大应力不应该出现在表面，而应该是表层下面，为什么总是在表面最大？请高手指教。
*SET,R1,0.2  !Define constant of sphere,R  
*SET,R2,0.1*R1 !Define constant of sphere,0.1R
*SET,R3,0.025*R1 !Define constant of sphere,0.025R
*SET,R4,0.01*R1 !Define constant of sphere,0.01R
*SET,Wc,1.0666e-4
*SET,Wn,1
*SET,Ac,4.6e-3
/PREP7  
ET,1,PLANE183   
MP,DENS,1,7930e6
MP,EX,1,1.813e8
MP,NUXY,1,0.3  
TB,BISO,1   
TBDATA,1,1.813e6 ! Yield stress (10^3Pa)  
TBDATA,2,3.626e6 ! Tangent modulus (10^3Pa)  
CSYS,1  
K,6
K,1,R1  
K,2,R1,90  
L,1,2
  
LOCAL,11,1,,R1,,-90
WPAVE,0,0,0
K,3,R2
K,4,R3
K,5,R4
L,6,3
L,3,4
L,4,5
L,5,2
K,7,R2,115
K,8,R3,115
K,9,R4,115
L,3,7
L,4,8  
L,5,9  
BTOL,0.1e-006,
FLST,2,4,4,ORDE,3   
FITEM,2,1   
FITEM,2,6   
FITEM,2,-8  
LOVLAP,P51X
LDELE,14,,,1
LDELE,15,,,1
LDELE,16,,,1
CSYS,0  
L,1,6
Al,1,9,11,2  
Al,11,17,12,3   
Al,12,18,13,4  
Al,13,10,5  
CSYS,1
MSHKEY,1
MSHAPE,0,2d
esize,0.02*Ac
MAT,1   
Amesh,4
esize,0.08*Ac
Amesh,3
MSHKEY,0
esize,R1/200
Amesh,2
MSHKEY,0
esize,R1/20
Amesh,1
CSYS,0
K,101,-0.1*R1,R1
K,102,0.8*R1,R1
LSTR,101,102
ET,2,169              ! Target169 2D Target Segment
ET,3,172
MP,MU,1,0
KEYOPT,3,5,3
TYPE,2
MAT,1
LESIZE,6,,,1
LMESH,6
KMESH,101
TYPE,3
R,1,,,10,0.1
REAL,1
LSEL,S,LINE,,10
LSEL,A,LINE,,18
NSLL,S,1
ESURF
LSEL,ALL
NSEL,ALL
!LSEL,,LOC,X,0
!DL,ALL,,UX,0
NSEL,S,LOC,X
D,ALL,UX,0
NSEL,ALL
NSEL,S,LOC,Y
D,ALL,UY,0
D,ALL,UX,0
NSEL,ALL
n_load=node(kx(101),ky(101),0)
SBCTRAN
FINISH
/SOLU
ANTYPE,STATIC   
nlgeom,on
solc,on                    ! N.L求解自动控制打开（缺省）
NROPT,FULL,,OFF  
time,1
!F,n_load,fy,-130.8701
d,n_load,uy,-Wn*Wc
nsubst,25,500,10
AUTOTS,1
outres,all,all
SOLVE

nickler

请大家帮忙啊

nickler

没人帮忙啊

fengerchai

帮你顶
顺便看结果

nickler

完了，没人看了，呜呜

aoshiqingyi

三年前的帖子了。
按照接触力学理论，剪应力是造成接触破坏的主要因素，对于两圆柱体接触情况，最大剪应力出现在接触中心下方，深度为接触面半宽的0.78倍，所以，要查看剪应力云图。个人简介，仅供参考。