CommandList-HW6img.txt

!*********************
! PROBLEM HOMEWORK 6 *
!*********************

FINISH
/CLEAR, START, NEW
/TITLE, HOMEWORK 6
/FILNAME,HOMEWORK6
/EFACET, 4
/GFILE, 2400
/PLOTS,DATE,0
/PNUM,LINE,1
/PNUM,AREA,1
! >>> PARAMETERS MODEL <<<
*SET,l, 50
*SET,g, 2
*SET,t, 1.5
*SET,Wp, 25
*SET,Rp, 4
*SET,Rd, 3
*SET,delta, 4*t
*SET,EPS, 1E-4
*SET,n_div, 12
*SET,e_lenght, 1
*SET,step,100
*SET,x, delta/step

! >>> PROPERTIES MATERIAL <<<
*SET,E_Young, 205E3
*SET,ni, 0.3
*SET,SigmaYield, 355
*SET,E_p, 4E3
*SET,friction, 0.1

! >>> LOAD <<<
*SET,F, 100E3

! >>> PRE PROCESSING <<<
/PREP7
ET,1,PLANE183,,,1

MP,EX,1,E_Young
MP,PRXY,1,ni
MP,MU,1,friction
TB,BKIN,1,1
TBTEMP,0.0
TBDATA,1,SigmaYield,E_p

!***PUNCH MODEL
K,1,0,0
K,2,Wp,0
K,3,Wp,Wp
K,4,0,Wp
L,1,2
L,2,3
L,3,4
L,4,1
LFILLT,1,2,Rp
AL,1,2,3,4,5

!***BLANCK HOLDER MODEL
K,7,Wp+g,0
K,8,Wp+g+(l/4),0
K,9,l+g,0
K,10,l+g,Wp/2
K,11,Wp+g,Wp/2
L,7,8
L,8,9
L,9,10
L,10,11
L,11,7
AL,6,7,8,9,10

!***DIE MODEL
K,12,Wp+g,-t
K,13,Wp+g+(l/4),-t
K,14,l+g,-t
K,15,l+g,-t-Wp
K,16,Wp+g,-t-Wp
L,12,13
L,13,14
L,14,15
L,15,16
L,16,12
LFILLT,15,11,Rd
AL,11,12,13,14,15,16
K,19,0,0
K,20,l,0
K,21,l,-t
K,22,0,-t
L,19,20
L,20,21
L,21,22
L,22,19
AL,17,18,19,20
SAVE
APLOT
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15   
/REPLOT   
/IMAGE, SAVE, HW6-Geometry, bmp
/RGB,INDEX, 0, 0, 0, 0  
/RGB,INDEX, 60, 60, 60,13   
/RGB,INDEX, 80, 80, 80,14   
/RGB,INDEX,100,100,100,15   
/REPLOT 

! >>> MESHING <<<
LESIZE,18,,,n_div
LESIZE,20,,,n_div

ESIZE,e_lenght
MSHAPE,0
MSHKEY,1
AMESH,4
SAVE

! >>> DEFINITION CONTACT <<<
ET,2,TARGE169  
ET,3,CONTA172
KEYOPT, 3, 2, 0
KEYOPT, 3, 7, 1  
KEYOPT, 3, 9, 2	  
KEYOPT, 3,10, 2

!****PUNCH - BLANK
!***DEFORMABLE BLANK
TYPE,3
REAL,1
LSEL,s,LINE,,17
NSLL,s,1
NSEL,R,LOC,x,0,Wp+g
ESURF
ALLSEL,all
SAVE

!***RIGIDBODY PUNCH
TYPE,2
REAL,1
TSHAP,LINE
LMESH,1
TSHAP,ARC
LMESH,5
SAVE
TSHAP,PILO 
KMESH,1

!***CONTACT BLANKHOLDER-BLANK
TYPE,3
REAL,2
LSEL,s,LINE,,17
NSLL,s,1
NSEL,R,LOC,x,Wp+g,l+g
ESURF
ALLSEL,all
SAVE

TYPE,2
REAL,2
TSHAP,LINE
LMESH,6,7
TSHAP,PILO
KMESH,8
SAVE

!***CONTACT DIE-BLANK
TYPE,3
REAL,3
LSEL,S,LINE,,19
NSLL,S,1
NSEL,R,LOC,x,Wp,l
ESURF
ALLSEL,all
SAVE

TYPE,2
REAL,3
TSHAP,LINE
LMESH,11,12
LSEL,S,LINE,,11,12
ESLL,S
ESURF,,REVERSE
ALLSEL,ALL
SAVE

TSHAP,ARC
LMESH,16
ESEL,S,,,684
ESURF,,REVERSE  
ALLSEL,ALL
TSHAP,PILO
KMESH,13
SAVE
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15   
/REPLOT   
/IMAGE, SAVE, HW6-Mesh, bmp
/RGB,INDEX, 0, 0, 0, 0  
/RGB,INDEX, 60, 60, 60,13   
/RGB,INDEX, 80, 80, 80,14   
/RGB,INDEX,100,100,100,15   
/REPLOT 

! >>>CONSTRAIN <<<
LSEL,S,,,20
NSLL,S
DSYM,SYMM,X
ALLSEL,ALL,ALL
SAVE

! >>> LOAD CONDITION <<<
!***LOAD STEP 1: PRELOAD
TIME,1 
F,1931,fy,-F 
OUTRES,ALL,ALL 
LSWRITE,1
!***LOAD STEP 2: DESCENT OF PUNCH
*DO,i,1,step
 TIME,1+i
 D,1926,uy,-x*i
 OUTRES,ALL,ALL 
 LSWRITE,1+i 
*ENDDO
!***LOAD STEP 3: ASCENT OF PUNCH
*DO,i,step+2,step*2
 TIME,i 
 D,1926,uy,(-x*step)+(i-step)*x
 OUTRES,ALL,ALL
 LSWRITE,i
*ENDDO
SAVE

! >>> SOLUTION <<<
/SOLU
SOLCONTROL,ON
NSUBST,10,100,2	 
NROPT,FULL,,OFF	 
AUTOTS,ON
EQSLV,PCG
NLGEOM,ON
LSSOLVE,1,2*step
FINISH

! >>> POST PROCESSING <<<
/POST1
*DO,N,1,2*step
SET,N
PLNSOL,S,EQV
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15   
/REPLOT   
/IMAGE, SAVE, HW6-Solu_%N%, bmp
*ENDDO
/RGB,INDEX, 0, 0, 0, 0  
/RGB,INDEX, 60, 60, 60,13   
/RGB,INDEX, 80, 80, 80,14   
/RGB,INDEX,100,100,100,15   
/REPLOT 

!***VON MISES EQUIVALENT STRESS AT THE END OF THE TRAVEL OF THE PUNCH
SET,step+1
PLNSOL,S,EQV 
PRRFOR,FY
*DO,N,1,2*step
 SET,N
 *GET,F,NODE,1926,RF,FY
 *CFOPEN,ForcePunch,txt,,APPEND 
 *VWRITE,N,F
 (F20.10,F20.10)
 *CFCLOS
*ENDDO

!***PUNCH STROKE THAT MAINTAINS THE MAXIMUM ABSOLUTE HOOP STRAIN BELOW 5%.
*DO,N,1,2*step
 SET,N
 PLNSOL,EPTO,Z
 *GET,MAXSTR,PLNSOL,0,MAX
 *GET,MINSTR,PLNSOL,0,MIN
 *CFOPEN,HoopStrain,txt,,APPEND
 *VWRITE,N,MAXSTR,MINSTR
 (F20.10,F20.10,F20.10)
 *CFCLOS
*ENDDO
SAVE
!***RESIDUAL STRESS DISTRIBUTION SURFACE
SET,103
PATH,TOP,2,,1000
PPATH,1,,0,0
PPATH,2,,L,0 
PDEF,TOP_SEQV,S,EQV
PLPATH,TOP_SEQV
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15   
/REPLOT   
/IMAGE, SAVE, HW6-TOP_SEQV, bmp
/RGB,INDEX, 0, 0, 0, 0  
/RGB,INDEX, 60, 60, 60,13   
/RGB,INDEX, 80, 80, 80,14   
/RGB,INDEX,100,100,100,15   
/REPLOT 
!***RESIDUAL STRESS DISTRIBUTION BOTTOM SURFACE
PATH,BOTTOM,2,,1000
PPATH,1,,0,-T
PPATH,2,,L,-T
PDEF,BOT_SEQV,S,EQV
PLPATH,BOT_SEQV
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15   
/REPLOT   
/IMAGE, SAVE, HW6-BOT_SEQV, bmp
/RGB,INDEX, 0, 0, 0, 0  
/RGB,INDEX, 60, 60, 60,13   
/RGB,INDEX, 80, 80, 80,14   
/RGB,INDEX,100,100,100,15   
/REPLOT 

!***AXIAL DISPLACEMENT AFTER THE PUNCH REMOVE
SET,step+1
NSEL,S,LOC,y,-(1.5/2)-0.05,-(1.5/2)+0.05
*GET,nnodi,NODE,,count
*DIM,dy,ARRAY,nnodi
*DIM,xn,ARRAY,nnodi 
CM,nodi,NODE
*do,i,1,nnodi
 *GET,xmin,NODE,,mnloc,x
 nsel,r,loc,x,xmin-1e-3,xmin+1e-3
 *GET,nmin,NODE,,NUM,MAX 
 xn(i)=xmin 
 *GET,dy(i),NODE,nmin,u,y
 CMSEL,S,nodi
 NSEL,U,NODE,,nmin
 CM,nodi,NODE
 *CFOPEN,set%step+1%,txt
 *VWRITE,xn(1),dy(1)
 (F20.10,F20.10)
 *CFCLOS
*ENDDO 
*DEL,ALL

SET,103
NSEL,S,LOC,y,-(1.5/2)-0.05,-(1.5/2)+0.05
*GET,nnodi,NODE,,COUNT
*DIM,dy,ARRAY,nnodi
*DIM,xn,ARRAY,nnodi 
CM,NODI,NODE
*DO,i,1,nnodi
 *GET,XMIN,NODE,,MNLOC,X
 NSEL,R,LOC,X,XMIN-1E-3,XMIN+1E-3
 *GET,NMIN,NODE,,NUM,MAX 
 xn(i)=XMIN
 *GET,DY(i),NODE,NMIN,U,Y
 CMSEL,S,NODI
 NSEL,U,NODE,,NMIN
 CM,NODI,NODE
 *CFOPEN,set103,txt
 *VWRITE,xn(1),dy(1)
 (F20.10,F20.10) 
 *CFCLOS
*ENDDO
FINISH
/EXIT,ALL