一个接触分析的命令流，带动画

nieliping0910

今天在一本书上见到一个典型的接触问题，做完了，发到网上和大家分享一下，命令流如下：
/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
!*

pfsr

感谢楼主分享

shin1209

收下看看,谢谢楼主分享

neumining

收下看看,谢谢楼主分享