fini
/cle

/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=1		!Abstand Kugeln
	T=0.1		!Eindrucktiefe


!*****************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.3
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
n,1000,DD/2+DK/2,90	!Knoten erzeugen mit Nummer 100000,x,y
e,1000			!Element durch den Knoten 100000
TSHAP,PILOT
csys,0
n,1001,0,0,0
e,1001



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

csys,0
nsel,s,loc,y,0		!alle Knoten bei y=0 neuselektieren
d,all,all		!alle Freiheitsgrade sperren
alls

! csys,1
! 
! nrotat,100000

fini

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

/solu			!Solutionprozessor starten

!---------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

*do,i,1,Z

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

	time,(i-1)*3+i+1		!Zurückfahren
	d,1001,uy,4*T
	solve

	time,(i-1)*3+i+2		!Verdrehen
	d,1001,rotz,-0.5*i
	solve

*enddo

fini

/post1
set,2,last
/dscale,,1
/pbc,rot,,2
plnsol,s,x
set,3,last
/rep
ANTIME,50,0.1, ,1,1,0,0