/PREP7
! Elementtype
ET,1,SOLID185
ET,2,CONTA173
!ET,2,CONTA173,,1 !für penalty func
ET,3,TARGE170

!Materialkennwert Stahl
MP,EX,1,2E5
MP,NUXY,1,0.3
!
mp,mu,3,0.1		!Kontaktreibung
!
TB,BISO,1		!data table for material properties, bilinear isotropic hardening using von mises or hill plastic
TBDATA,1,200,1000


!Real constants
R,1
R,2			
R,3

!Geometrie
blc4,0,0,100,100,100	!Würfel erstellen mit 100x100x100
esize,10		!Elementgroesse 10
VMES,ALL		
!local,11,0,50,50,100	!lokales Koordinatensystem definieren, 11=Nummer, 0=Cartesian, x, y, z
*GET,NMAX,NODE,,NUM,max !größte verwendete Knotenzahl auswählen
I=Nmax+1		!Variable Pilot 1 definieren
J=Nmax+2
K=Nmax+3
L=Nmax+4
n,I,120,-20,100.3		!I
n,J,-20,-20,100.3
n,K,-20,120,100.3
n,L,120,120,100.3

!Target festelgen
nsel,s,loc,z,100	!Alle Knoten mit z=100 auswählen
TYPE,2
Real,2			
ESURF
allsel
TYPE,3
TSHAP,Quad		!defines simple 2D and 3D geometric surfaces for target segment elements
E,I,J,K,L		!Element an Pilotknoten definieren

asel,s,,,4
nsla,s,0
d,all,UX,0,0,,,UY
allsel
asel,s,,,5
nsla,s,0
d,all,UX,0,0,,,UY
allsel

CNCHECK,all

/SOLU
NLGEOm,ON
SOLCON,ON
NEQIT,50
EQSLV,FRONT
Autots,on
OUTRES,ALL,-5
CSYS

time,1
asel,s,,,1
nsla,s,1
d,all,uz,2
allsel
solve

time,2
asel,s,,,1
nsla,s,1
ddele,all,all
allsel
asel,s,,,1
nsla,s,1
f,all,fz,100
nsubst,1000
!KBC,1
allsel
solve

/POST1
/EDGE,,1
/VIEW,,3,1,1
ESEL,U,TYPE,,2,3
EPLO
SET,LAST
/DSCA,,1
/GRA,FULL
PLNS,S,EQV