finish
/clear

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

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

!Geometrie
D2=1		!Kugeldurchmesser
D1=6		!Drahtdurchmesser

!Material
E1=206000	!Emodul
V1=0.3		!Querkontraktionszahl
E2=300000
V2=0.3


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

/prep7

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

et,1,plane182,,,2	!Elemente Grundkörper
et,2,plane182,,,2

et,3,conta171,,,2	!Kontaktelemente Draht/Kugel
et,4,targe169		

et,5,conta171,,2,,1!KOntakt zur Steuerung der Verdrehung
KEYOPT,5,12,5      !Bonded always
et,6,targe169,,0,,111		
!-------------Material------------------------------

mp,ex,1,E1
mp,nuxy,1,V1

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

mp,ex,2,E2
mp,nuxy,2,V2

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

cyl4,0,0,,0,D1/2,180			!Kreis,x,y,Radius innen,Startwinkel,Radius außen, Endwinkel
cyl4,0,D1/2+D2/2,D2/2,180,,360
!----------Vernetzung-------------------------------

type,1
mat,1
esize,0.3
amesh,1			!Vernetzung von Draht

type,2
mat,2
esize,.1
amesh,2    !Vernetzung Kugel

!----Erstellung von Komponeneten--------------------

local,11,1              !Definition von KOS
local,12,1,,D1/2+D2/2

csys,11
esel,s,mat,,1
nsle
nsel,r,loc,x,D1/2
nsel,r,loc,y,60,120
cm,TARGET,node

csys,12
esel,s,mat,,2
nsle
nsel,r,loc,x,D2/2
nsel,r,loc,y,240,300
cm,CONTACT,node

csys,0
nsel,s,loc,y,0
cm,SLAVE,node

nsel,s,loc,y,D1/2+D2/2
cm,ANGRIFF,node
allsel



!------KONTAKTDEFINITION KUGEL/DRAHT Kontakt-----------------

type,3
real,3
cmsel,s,CONTACT
esurf

type,4
cmsel,s,TARGET
esurf

allsel

!-----Kontakt Pilot und untere Linie--------------------

type,6
real,5
r,5

*get,nmax,node,,num,max
PILOT_N=nmax+1
n,PILOT_N
tshap,pilot
e,PILOT_N

NSEL,S,,,SLAVE  
TYPE,5  
ESLN,S,0
ESURF   
ALLSEL

finish

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

/solu

ddele,all,all
d,Pilot_N,all,0

nlgeom,on		!Berücksichtigung geometrischer Nichtlinearitäten 
nsubst,100
autots,on		!Automatische Zeitschrittsteuerung aktivieren
neqit,500		!maximale 500 Iterationsschritte zum erreichen der Lösung
outres,all,all

time,1
d,ANGRIFF,uy,-0.00001
solve

time,2
d,ANGRIFF,uy,0.1
solve

time,3
d,PILOT_N,rotz,-0.05
solve


!********************Postprozessor

/post1
set,last
plnsol,s,eqv
	
ANTIME,50,0.1, ,1,1,0,0