两圆球的Hertz接触
求解结果错误，请大家帮帮忙，谢谢！/BATCH
/COM,ANSYS RELEASE 10.0 UP20050718 22:00:12 05/29/2008
/input,menust,tmp,'',,,,,,,,,,,,,,,,1
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1
/REPLOT,RESIZE
WPSTYLE,,,,,,,,0
/REPLOT,RESIZE
/PREP7
!*
ET,1,SOLID185
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,30E3
MPDATA,PRXY,1,,0.3
SPH4, , ,10
SPH4, ,20-1e-5,10
MSHAPE,1,3D
MSHKEY,0
!*
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,5,1,5,ORDE,1
FITEM,5,1
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
CMDELE,_Y
!*
!*
AREFINE,_Y1, , ,2,0,1,1
CMDELE,_Y1
!*
FLST,5,1,5,ORDE,1
FITEM,5,4
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
CMDELE,_Y
!*
!*
AREFINE,_Y1, , ,2,0,1,1
CMDELE,_Y1
!*
!*
!*
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,
MAT,1
R,3
REAL,3
ET,2,170
ET,3,174
KEYOPT,3,9,0
KEYOPT,3,10,2
R,3,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
ASEL,S,,,1
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,4
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
/MREP,EPLOT
EPLOT
ALLSEL,ALL
/UI,MESH,OFF
K, ,0,30-1e-5,,
klist,all,,,coord
NSEL,S,LOC,X,
FINISH
/SOL
!*
ANTYPE,0
FLST,2,4,1,ORDE,4
FITEM,2,1
FITEM,2,8
FITEM,2,348
FITEM,2,355
!*
/GO
D,P51X, , , , , ,UX, , , , ,
NSEL,S,LOC,Y,
FLST,2,72,1,ORDE,50
FITEM,2,1
FITEM,2,-24
FITEM,2,744
FITEM,2,-745
FITEM,2,1000
FITEM,2,-1001
FITEM,2,1004
FITEM,2,-1005
FITEM,2,1104
FITEM,2,-1105
FITEM,2,1121
FITEM,2,-1122
FITEM,2,1197
FITEM,2,-1198
FITEM,2,1208
FITEM,2,-1209
FITEM,2,1213
FITEM,2,-1214
FITEM,2,1219
FITEM,2,-1220
FITEM,2,1225
FITEM,2,-1226
FITEM,2,1231
FITEM,2,-1232
FITEM,2,1278
FITEM,2,-1279
FITEM,2,1289
FITEM,2,-1290
FITEM,2,1297
FITEM,2,-1298
FITEM,2,1310
FITEM,2,-1311
FITEM,2,1319
FITEM,2,-1320
FITEM,2,1325
FITEM,2,-1326
FITEM,2,1335
FITEM,2,-1336
FITEM,2,1354
FITEM,2,-1355
FITEM,2,1365
FITEM,2,-1366
FITEM,2,1497
FITEM,2,-1498
FITEM,2,1523
FITEM,2,-1524
FITEM,2,1558
FITEM,2,-1559
FITEM,2,1670
FITEM,2,-1671
!*
/GO
D,P51X, , , , , ,UY, , , , ,
NSEL,S,LOC,Z,
FLST,2,4,1,ORDE,4
FITEM,2,2
FITEM,2,14
FITEM,2,349
FITEM,2,361
!*
/GO
D,P51X, , , , , ,UZ, , , , ,
FLST,2,1,3,ORDE,1
FITEM,2,9
!*
/GO
FK,P51X,FY,-3000
/STATUS,SOLU
! SOLVE
ALLSEL,ALL
/STATUS,SOLU
SOLVE
FINISH
/POST1
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0 Verification manual中有Hertz接触的例子 其中的VM191就是你所说的两圆球接触的问题，你自己可以看看！ 谢谢！VM191是两圆柱体的接触，圆柱体的各个截面是相同的，但圆球的截面不太一样。我个人觉得是不是加约束的时候存在问题，还在考虑！
回复 1# 的帖子
楼主能加我好友么？我也在做这个问题396127582 注明 赫兹接触 :handshake 同仁 以后 多 请教赫兹接触
这个问题放了很长一段时间了。之前计算了最大接触应力和接触半径，最大接触应力1122Mpa和理论值1848MPa相差较大，接触半径0.7902与理论值0.8804相差不大。请各位多多指导啊！页:
[1]