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

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

rad=50 !Radius der Kugel

!Real constants
R,2,rad			!bei Kugel rad

!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
Pilot1=Nmax+1		!Variable Pilot 1 definieren
n,pilot1,,,rad		!Knoten definieren, z-Wert Radius der Kugel

!Target festelgen
TYPE,3
REAL,2
TSHAP,SPHE		!defines simple 2D and 3D geometric surfaces for target segment elements
E,Pilot1,		!Element an Pilotknoten definieren
nsel,s,loc,z,0		!Alle Knoten mit z=0 auswählen
TYPE,2			!Kontaktelemente definieren? z=0 wegen obigem Koordinatensystem
ESURF			!vernetzen

ALLS

/SOLU
NLGEOm,ON
SOLCON,ON
NEQIT,50
EQSLV,FRONT
NSUB,50,200
OUTRES,ALL,-5
CSYS

NSEL,S,LOC,Z
D,ALL,all
ALLS
D,Pilot1,ALL
D,Pilot1,UZ,-10		!hier aufbringen der Verschiebung auf den Pilotknoten
ALLS
SOLVE

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