用ansys作abaqus前处理

/COM, ==========================================================
/COM, Output ABAQUS/Standard input file
/COM, ==========================================================

*CFOPEN,punch,inp !打开ABAQUS的输入INP文件

*VWRITE
('*HEADING')
*VWRITE
('INDENTATION OF A VISCOELASTIC FOAM SPECIMEN')
*VWRITE
('*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO,CONTACT=NO')

*VWRITE
('*NODE,NSET=FOAM')
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM,NX(NNUM),NY(NNUM),NZ(NNUM)
(F7.0,TL1,3(',',E16.9))
*ENDDO

*VWRITE
('*ELEMENT,TYPE=CAX4,ELSET=FOAM')
*GET,NELE,ELEM,,COUNT
ENUM = 0
*DO,I,1,NELE,1
ENUM = ELNEXT(ENUM)
NI = NELEM(ENUM,1)
NJ = NELEM(ENUM,2)
NK = NELEM(ENUM,3)
NL = NELEM(ENUM,4)
*VWRITE,ENUM,NI,NJ,NK,NL
(F7.0,TL1,4(',',F7.0,TL1),' ')
*ENDDO

*GET,NMAX,NODE,,NUM,MAX
NREF = NMAX+1
*VWRITE
('*NODE,NSET=NREF')
*VWRITE,NREF,0.0,FTHK+PRAD,0.0
(F7.0,TL1,4(',',F7.0,TL1),' ')

*VWRITE
('*NSET,NSET=NBOT')
LSEL,S,LOC,Y,0.0
NSLL,S,1
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,',')
*ENDDO
ALLSEL

*VWRITE
('*NSET,NSET=NTOP')
LSEL,S,LOC,Y,FTHK
NSLL,S,1
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,',')
*ENDDO
ALLSEL

*VWRITE
('*NSET,NSET=NCEN')
LSEL,S,LOC,X,0.0
NSLL,S,1
*GET,NNOD,NODE,,COUNT
NNUM = 0
*DO,I,1,NNOD,1
NNUM = NDNEXT(NNUM)
*VWRITE,NNUM
(F7.0,TL1,',')
*ENDDO
ALLSEL

*VWRITE
('*ELSET,ELSET=ETOP')
LSEL,S,LOC,Y,FTHK
NSLL,S,1
ESLN
*GET,NELE,ELEM,,COUNT
ENUM = 0
*DO,I,1,NELE,1
ENUM = ELNEXT(ENUM)
*VWRITE,ENUM
(F7.0,TL1,',')
*ENDDO
ALLSEL

*VWRITE
('*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM,ORIENTATION=RECT')
*VWRITE
('*ORIENTATION,NAME=RECT,SYSTEM=RECTANGULAR')
*VWRITE
('1.0,0.0,0.0,0.0,1.0,0.0')
*VWRITE
('3,0.0')

*VWRITE
('*MATERIAL,NAME=FOAM')
*VWRITE
('*HYPERFOAM,N=2,TEST DATA INPUT,POISSON=0.0,MODULI=INSTANTANEOUS')
*VWRITE
('*UNIAXIAL TEST DATA')
*VWRITE,SUNI(1),EUNI(1)
(E16.9,',',E16.9)
*VWRITE
('*SIMPLE SHEAR TEST DATA')
*VWRITE,SSHR(1),ESHR(1),TSHR(1)
(E16.9,2(',',E16.9))
*VWRITE
('*VISCOELASTIC,TIME=RELAXATION TEST DATA')
*VWRITE
('*SHEAR TEST DATA,SHRINF=0.50')
*VWRITE,SRLX(1),TRLX(1)
(E16.9,',',E16.9)

*VWRITE
('*SURFACE,NAME=FSRF')
*VWRITE
('ETOP,S3')

*VWRITE,NREF
('*RIGID BODY,ANALYTICAL SURFACE=RSRF,REF NODE=',F7.0,TL1,' ')

*VWRITE
('*SURFACE,TYPE=SEGMENTS,NAME=RSRF')
*AFUN,DEG
*VWRITE,(PRAD*COS(45)),FTHK+PRAD+(PRAD*SIN(45))
('START',2(',',E12.4))
*VWRITE,-1.0,FTHK,0.0,FTHK+PRAD
('CIRCL',4(',',E12.4))

*VWRITE
('*CONTACT PAIR,INTERACTION=ROUGH')
*VWRITE
('FSRF,RSRF')

*VWRITE
('*SURFACE INTERACTION,NAME=ROUGH')
*VWRITE
('*FRICTION')
*VWRITE,CFRC
(F7.4,',')

*VWRITE
('*BOUNDARY')
*VWRITE
('NBOT,1,2,0.0')
*VWRITE
('NCEN,1, ,0.0')
*VWRITE

('NREF,1, ,0.0')
*VWRITE
('NREF,6, ,0.0')

*VWRITE
('*STEP,NLGEOM=YES,INC=200,AMPLITUDE=RAMP,UNSYMM=YES')
*VWRITE
('*VISCO,CETOL=0.01')
*VWRITE,TINI,TTOT,TMAX
(E12.4,',',E12.4,',,',E12.4)

*VWRITE
('*PRINT,CONTACT=YES')
*VWRITE
('*CONTACT CONTROLS,FRICTION ONSET=DELAY')
*VWRITE
('*CONTACT FILE,SLAVE=FSRF,FREQUENCY=10')

*VWRITE
('*BOUNDARY')
*VWRITE,PDSP
('NREF,2,,',E12.4)

*VWRITE
('*RESTART,WRITE,FREQUENCY=10')

*VWRITE
('*OUTPUT,FIELD,FREQUENCY=5,VARIABLE=PRESELECT')
*VWRITE
('*OUTPUT,HISTORY,FREQUENCY=1,VARIABLE=PRESELECT')
*VWRITE
('*CONTACT OUTPUT,NSET=NTOP,VARIABLE=PRESELECT')
*VWRITE
('*NODE OUTPUT,NSET=NREF')
*VWRITE
('U,RF')

*VWRITE
('*END STEP')

*VWRITE
('*STEP,NLGEOM=YES,INC=200,AMPLITUDE=RAMP,UNSYMM=YES')
*VWRITE
('*VISCO,CETOL=0.01')
*VWRITE,TINI,TTOT,TMAX
(E12.4,',',E12.4,',,',E12.4)

*VWRITE
('*BOUNDARY,OP=MOD')
*VWRITE
('NREF,2,,0.0')

*VWRITE
('*END STEP')

*CFCLOS !Close command file

FINI !Exit preprocessor

/COM, ==========================================================
/COM, 运行ABAQUS求解模型
/COM, ==========================================================

/SYS,abaqus job=punch

/EXIT !退出ANSYS