按道理最大应力不应该出现在表面,而应该是表层下面,为什么总是在表面最大?请高手指教。
*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 |