今天在一本书上见到一个典型的接触问题,做完了,发到网上和大家分享一下,命令流如下:
/PREP7
BLOCK,-2,2,-2,2,2.5,3.5,
CYLIND,0.49, ,2.5,3.5,0,360,
VSBV, 1, 2
CYLIND,0.5, ,2,4.5,0,360,
/PNUM,KP,0
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0
!*
/PNUM,ELEM,0
/REPLOT
!*
wpstyle,0.05,0.1,-1,1,0.003,0,0,,5
WPSTYLE,,,,,,,,1
wpro,,,90.000000
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,3
VSBW,P51X
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,6
VDELE,P51X, , ,1
WPCSYS,-1,0
wpro,,90.000000,
FLST,2,2,6,ORDE,2
FITEM,2,4
FITEM,2,-5
VSBW,P51X
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,3
VDELE,P51X, , ,1
SAVE
!*
ET,1,SOLID185
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,36E6
MPDATA,PRXY,1,,0.3
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
FLST,5,2,4,ORDE,2
FITEM,5,4
FITEM,5,10
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,3, , , , ,0
!*
FLST,5,1,4,ORDE,1
FITEM,5,18
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,4, , , , ,0
!*
FLST,5,2,6,ORDE,2
FITEM,5,2
FITEM,5,6
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
/UI,MESH,OFF
!*
/SHRINK,0
/ESHAPE,0.0
/EFACET,2
/RATIO,1,1,1
/CFORMAT,32,0
/REPLOT
!*
!*
!*
/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,0.2
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,,,25
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,,,21
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
APLOT
FINISH
/SOL
FLST,2,4,5,ORDE,4
FITEM,2,3
FITEM,2,8
FITEM,2,10
FITEM,2,27
DA,P51X,SYMM
FLST,2,1,5,ORDE,1
FITEM,2,13
!*
/GO
DA,P51X,ALL,
ANTYPE,0
NLGEOM,1
NSUBST,100,1000,10
OUTRES,ERASE
OUTRES,ALL,ALL
AUTOTS,1
PSTRES,1
LSWRITE,1,
ALLSEL,ALL
APLOT
NSEL,S,LOC,Z,4.5
FLST,2,18,1,ORDE,4
FITEM,2,37
FITEM,2,-42
FITEM,2,73
FITEM,2,-84
!*
/GO
D,P51X, ,1.7, , , ,UZ, , , , ,
ALLSEL,ALL
NSUBST,1000,10000,10
LSWRITE,2,
LSSOLVE,1,2,1
FINISH
/POST1
SET,FIRST
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0
SET,1,LAST,1,
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0
SET,LAST
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0
SET,1,LAST,1,
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0
!*
/EXPAND,4,POLAR,HALF,,90
/REPLOT
!*
SET,LAST
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0
!*
PLNS,S,EQV
ANDATA,0.5, ,1,1,2,1,0,1
!* |