做的一个接触分析的问题,虽然计算收敛了,但是感觉不对,但又找不到到底那块不妥?还请高手们帮忙看看,多谢了。
/PREP7
ET,1,182,1,,2
ET,2,169
ET,3,171
R,2,0,0,1,0.1
RMORE,0.1,0.01
!*
!* Material Properties
!*
MP,EX,1,110E9
MP,NUXY,1,0.3
MP,EX,2,439.43E9
MP,NUXY,2,0.3
MP,MU,2,0.3
!*
!* keypoints
!*
K,1
K,2,100
K,3,100,90
K,4,0,90
K,5,0,120
K,6,100,120
!*
!* areas
!*
A,1,2,3,4
TYPE,1
MAT,1
A,3,4,5,6
TYPE,1
MAT,2
!*
!* glue areas
!*
AGLUE,1,2
!*
!* mesh
!*
LESIZE,4,,,9
LESIZE,2,,,9
LESIZE,3,,,20
LESIZE,1,,,20
AMESH,1
LESIZE,5,,,6
LESIZE,7,,,6
LESIZE,6,,,20
AMESH,2
!*
!* indenter
!*
K,20,50,170
KWPAVE,20
CSWPLA,11,1
K,21,50,180
K,22,50,-90
K,23,50,0
K,24,0,0
L,21,22
L,22,23
TYPE,2
REAL,2
LMESH,8,9
KMESH,24
KSEL,S,KP,,24
NSLK,S
*GET,N_LOAD,NODE,,NUM,MAX
ALLS
CSYS,0
NSEL,S,LOC,Y,120
TYPE,3
ESURF
!*
!* Boundary Conditions
!*
NSEL,S,LOC,Y,0
D,ALL,UY
NSEL,S,LOC,X,50
NSEL,R,LOC,Y,0
D,ALL,UX
NSEL,ALL
FINISH
/SOLU
ANTYPE,STATIC
NLGEOM,ON
SOLC,ON
NSUBST,20,500,10
OUTRES,ALL,ALL
F,N_LOAD,FY,-20000 |