ex2-3-1_3d.txt

《ANSYS工程分析进阶实例》-王呼佳-源代码 看文件吧.不多说了.

fini
/cle

/FILNAME,rebound_3D
/PREP7  
*SET,exx_p,7e4
*SET,exx_t,7e10  
*SET,rp,65 
*SET,rt,5  
*SET,h,6.5  
*SET,sts1,55
*SET,stn1,sts1/exx_p  
*SET,sts2,112   
*SET,stn2,0.00575   
*SET,sts3,172   
*SET,stn3,0.02925   
*SET,sts4,241     
*SET,stn4,0.1   

ET,1,solid185
et,2,targe170   
et,3,conta174
ET,4,plane42
MP,EX,1,exx_p
MP,PRXY,1,0.325
MP,EX,2,exx_t
MP,PRXY,2,0.3
TB,KINH,1,1,4,0 
TBPT,,stn1,sts1 
TBPT,,stn2,sts2 
TBPT,,stn3,sts3 
TBPT,,stn4,sts4

CYL4,0,0,rp
RPR4,4,0,0,rp/2
/PNUM,KP,1  
/PNUM,LINE,1
/PNUM,AREA,1
AOVLAP,all
L,2,6  
L,1,5  
L,4,8  
L,3,7
csys,1  
LSEL,S,LOC,Y,0  
LSEL,A,LOC,Y,90 
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
ASBL,3,all    
LESIZE,all,,,12
LSEL,S,LOC,X,rp
LESIZE,all,,,20
LSEL,S,LOC,X,,rp/2
LESIZE,all,,,12
csys
type,4
MSHAPE,0,2D 
MSHKEY,0
AMESH,all
TYPE,1  
EXTOPT,ESIZE,2
EXTOPT,ACLEAR,1 
VEXT,all,,,0,0,h
allsel,all

wpoff,0,0,-0.1  
CYL4,0,0,rt
RPR4,4,0,0,rt/2
asel,s,loc,z,-100,-0.1
AOVLAP,all
L,18,22  
L,17,21  
L,20,24  
L,19,23

LSEL,s,LOC,z,-100,-0.1  
csys,1
lsel,u,loc,x,rt
lsel,u,loc,x,0,rt/2
asel,s,loc,z,-100,-0.1
ASBL,all,all    
LESIZE,all,,,3
LSEL,S,LOC,X,rt
LESIZE,all,,,5
LSEL,S,LOC,X,,rt/2
LESIZE,all,,,3
csys
type,4
MSHAPE,0,2D 
MSHKEY,0
AMESH,all
TYPE,1 
mat,2 
EXTOPT,ESIZE,2
EXTOPT,ACLEAR,1 
VEXT,all,,,0,0,-h
wpoff,0,0,0.1
allsel,all

MP,MU,1,0.35
MAT,1   
R,3 
REAL,3  
R,3,,,1.0,0.1
ASEL,S,loc,z,-0.1,-0.1 
TYPE,2  
NSLA,S,1
ESLN,S,0
ESURF,all
allsel,all
ASEL,S,,,2  
TYPE,3  
NSLA,S,1
ESLN,S,0
ESURF,all
ALLSEL,all

nsel,s,loc,z,-100,-0.1
D,all,ux
D,all,uy
CSWPLA,11,1,1,1
NROTAT,all
D,all,ux
D,all,uy
nsel,s,loc,z,h/2,h/2
nsel,r,loc,x,rp
NROTAT,all
d,all,uz
d,all,ux
ALLSEL,all
fini

/sol
csys
antype,0
OUTRES,ALL,all
nlgeom,on
autots,1
time,1       
nsubst,100
nsel,s,loc,z,-100,-0.1
d,all,uz,4
allsel,all
solve
*get,def1,node,413,u,z
time,2
nsubst,100
nsel,s,loc,z,-100,-0.1
d,all,uz,0
allsel,all
solve
*get,def2,node,413,u,z
disp_3d=def1-def2
fini

/post1
! 冲压完毕
set,1,last
PLNSOL, S,EQV
!切片
wprota,,90
/type,1,1
/cplane,1
PLNSOL, S,EQV
PLNSOL, NL,EPEQ		!等效塑性应变

! 卸载完毕
set,2,last
PLNSOL, S,EQV
PLNSOL, NL,EPEQ		!等效塑性应变

/type,1,6
/cplane,0
!找到最大应力的节点编号，75号节点
NSORT,S,EQV 
*GET,n_max_3d,SORT,0,IMAX

/post26
ANSOL,2,n_max_3d,S,EQV,eqv_3d
plvar,2