!Ebene spannung (ET,1, plane183,,, 2)

hoz=40 			!Kantenlänge Zylinder
reib1=0.4 			!Reibungskoeffizient Material1
abst=0			!Abstand Indenter zu  Oberfläche
bweg=1.7			!Bewegung des Indenters
weg=bweg-abst 		!Eindringweg 
fink=1				!Feinheit des Netzes um Eindringpunkt0.05
fdge=8			!Feinheit des Netzes 
zst=5				!Feinheit des Netzes an den Rändern des Rechtecks2
/PREP7		!Preprocessor

!---Elementtypen---
ET,1, plane183,,, 2 			!plane strain
ET,2,CONTA172 ,,,1
KEYOPT, 2, 9, 1		! keypt9/1 Exclude both initial geometrical penetration !or gap and offset

KEYOPT, 2, 12, 2	!keypt12/2 means no separation but sliding permitted



ET,3,TARGE169,,1,,,


!---Materialkennwert Stahl---
MP,EX,1,2.1E5			!E-Modul
MP,NUXY,1,0.3			!Querkontraktionszahl
MP, MU, 1, reib1, 			!Reibungskoeffizient 



!Spannungs-DehnungsDiagramm
TB,MISO,1,1,7                		! Activate a data table
TBTEMP,20                		! Temperature = 20
TBPT,DEFI,1E-3, 210    		! Strain, stress at temperature = 20
TBPT,DEFI,2E-3, 235    
TBPT,DEFI,2.6E-2, 236
TBPT,DEFI,5E-2, 260
TBPT,DEFI,10E-2, 310
TBPT,DEFI,15E-2, 350 
TBPT,DEFI,17.5E-2, 360
TBPLOT,MISO,1



!---Geometrie----
!---Keypoints---
K, 1, 4.195, 4.85+abst, 
K, 2, -4.195, 4.85+abst, 
K, 3, -0.20977, 0.1+abst, 
K, 4, 0.20977, 0.1+abst, 

!---Lines
LSTR, 1, 2
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 1

AL, 1, 2, 3, 4, 
K, 5, 0, abst, 
BSPLIN, 3, 5, 4,
AL, 3, 5,  
AADD, ALL, , , , , , , , 

!----------Teilen der Indenterfläche---
K, 100, , 4.85+abst, 
LSTR, 5, 100
ASBL, 3, 3, --, KEEP, KEEP

!--------Subtraktion der geteilten Fläche
ADELE, 1, , , 1
ADELE, 3, , , 1

!---Rechteck---
rectng,0,hoz,0,-hoz 

!---Zugstab---

rectng,-20 ,-15,50,25 


!----Vernetzung------Probe---
esize, fdge
type,1
mat,1
amesh,1

EREFINE, 21, , , 2, 2, off, off

EREFINE, 112, , , 2, 5, off, off


!----Vernetzung------Zugstab---
esize, zst
type,1
mat,1
amesh,3


!---Kontakte----
type,3			! Type-Befehl bezieht sich auf ET; TARGE169 wird aktiviert
TSHAP,ARC		!Bogen
LREVERSE, 9, 0	!Umkehrung der Normale
lmesh,9
TSHAP, LINE	!Gerade
lmesh, 4

*GET,NMAX,NODE,,NUM,max
Pilot1= NMAX+1
n, Pilot1,, 0.25
TSHAP,PILOT
e,Pilot1 

!-------------------sel line4, sel nodes, overlais elements over sel line
lsel,s,,,5
nsll,s,1
type,2
esurf

!----Randbedingung----- 

!----Randbedingung-----Zugstab
nsel,s,loc,y,50 
d,all,uy,0
d,all,ux,0

!----Randbedingung-----Rest

nsel,s,loc,y,-hoz 
d,all,uy,0
d,all,ux,0
nsel,s,loc,x,0
d,all,ux,0


!---Lasten----- 
!d=dof s for pilot1 nur y in -weg
d,pilot1,uy,-weg
ALLSEL, , 

!---Lasten----- Zugstab
nsel,s,loc,y,25
F, all, FY, -800, 


/SOLU		!Solution

CNCHECK, detail,1 , , 


! analysis, EQSLV=type of equation solver(frontaldirect)
! outres=solution data written to db

SOLCON,ON
!NLGEOm,ON	! nlgeom includes large deflections(wtf), NLGEOM,ON for large ! ! displacement/strain 

CUTCONTROL, PLSLIMIT  , 0.15,

!NEQIT, 500		! Maximum number of equilibrium iterations allowed each 
! substep. OFF
EQSLV,FRONT	! Frontal direct equation solver
NSUBST, 1000,	! Specifies the number of substeps to be taken this load step. 

OUTRES,ALL,all
CSYS 
ALLSEL, ,
SOLVE 

/POST1		!Postprocessor
/EDGE,,1
/DSCA,,1
PLNSOL,S,EQV	!plnsol,s equivalent stress (wtf),nodal sol v mise