/title,mehrere Kugeln (starr) 2D
/filename,Mehrfacheindruck

!******************Variablen**********************

!Geometrie
	DK=1		!Kugeldurchmesser
	DD=6		!Drahtdurchmesser

!Material
	E=206000	!Emodul
	V=0.3		!Querkontraktionszahl

!Parameter

	Z=3		!Anzahl Kugeln
	A=0.5		!Abstand Kugeln
	T=0.1		!Eindrucktiefe
	Abstand=0.02

!*****************Prep7******************************

/prep7

!--------------Elemente-----------------------------

et,1,plane182,,,2	!EbenerSpannungszustand
et,2,conta171,,,2	!Ebener Spannungszustand
et,3,targe169		!2D Zielkontaktfläche

!-------------Material------------------------------

mp,ex,1,E
mp,nuxy,1,V

tb,biso,1
tbdata,1,1200,E1/100

!----------Geometrievereinbarung--------------------

r,2,DK/2		!Realkonstante 1 von Element 3, Zielkontaktfläche mit Kreisradius D/2

!----------Modellgenerierung

cyl4,0,0,,0,DD/2,180			!Kreis,x,y,Radius innen,Startwinkel,Radius außen, Endwinkel

!----------Vernetzung-------------------------------

type,1
mat,1
esize,0.1
amesh,1			!Vernetzung von Rechteck

!--------Kontaktdefinition--------------------------

type,2			!Elementyp 2 aktivieren
real,2			!Festlegen der Realkonstante für nachfolgende Elemente
lsel,s,line,,1		!Linie 3 neuselektieren
nsll,s,1		!alle Knoten die vollständig zu den ausgewählten Linien gehören
esurf			!Erzeugen von Kontaktelementen auf den Flächenelementen von Linie 1

!--------Zielkontaktdefinition------------------------

csys,1
type,3			!Elemttyp 3 aktivieren
tshap,circle		!Zuweisen eines Kreises als 2D Geometrie der Zielkontaktfläche
alls
*GET,NMAX,NODE,,NUM,MAX

Circ_p=Nmax+1
Pilot1=Nmax+2
Pilot2=Nmax+3

n,Circ_p,DD/2+DK/2+Abstand,90	!Knoten erzeugen mit Nummer 100000,x,y
e,Circ_p		!Element durch den Knoten 100000
TSHAP,PILOT
csys,0
n,Pilot1,0,0,0
e,Pilot1



!-------Randbedingungen und Lasten-------------------------

csys,0
nsel,s,loc,y,0		!alle Knoten bei y=0 neuselektieren
cm,tar_2,node
n,Pilot2,

MAT,1   
R,3 
REAL,3  
ET,4,169
ET,5,172
KEYOPT,5,12,5   
KEYOPT,5,4,1
KEYOPT,5,2,2
KEYOPT,4,2,0
KEYOPT,4,4,111  
TYPE,4  
! Create a pilot node   
TSHAP,PILO  
E,Pilot2  
! Generate the contact surface  
NSEL,S,,,TAR_2  
TYPE,5  
ESLN,S,0
ESURF   
ALLSEL  
! d,all,all		!alle Freiheitsgrade sperren
! alls
! /eof
! csys,1
! 
! nrotat,100000

fini

!**********************Solution*************************

/solu			!Solutionprozessor starten
ddele,all,all
d,Pilot1,all
d,Pilot2,all
!---------Lösungseinstellungen für nichtlineare Berechnung und Lastschritte----------

nlgeom,on		!Berücksichtigung geometrischer Nichtlinearitäten 

nsubst,100
! deltim,0.01,0.01,0.5	!Substepeinstellungen
autots,on		!Automatische Zeitschrittsteuerung aktivieren
neqit,500		!maximale 500 Iterationsschritte zum erreichen der Lösung
outres,all,all		!Abspeichern alles Zwischenergebnisse all 10 Schritte

I=1

wink_=10			!winkel in Grad
wink_rad=wink_*PI/180
wink_akt=0
*afun,rad

*do,i,1,Z

	time,(i-1)*3+i			!Eintauchen
	d,Pilot1,uy,-T
	solve

	time,(i-1)*3+i+1		!Zurückfahren
	d,Pilot1,uy,0

	solve
	
	time,(i-1)*3+i+2		!Verdrehen
	d,Pilot2,rotz,wink_rad*i
	solve
*enddo
/post1
set,last
plnsol,s,eqv