!PREPROCESSING

/filnam,Ringfeder

!---Geometriedaten---
!Werte für Fmax=100
Di=75.5
De=90
se=2
E=11
b=18
pi=3.14159
phi=2*pi*12/360
x1=E/2*tan(phi)
x2=(E-2*se)/2*tan(phi)


!------------------------------PREP7------------------
/prep7

!---Materialvereinbarung---
mp,ex,1,210000		!E-Modul für Material1
mp,nuxy,1,0.3

mp,ex,2,210000		!E-Modul für Material2
mp,nuxy,2,0.3

mp,ex,3,210000		!E-Modul für Material3
mp,nuxy,3,0.3

mp,ex,4,210000		!E-Modul für Material4
mp,nuxy,4,0.3

!-----Contact-Reibungszahl
mp,mu,5,0.1									

!---Elementvereinbarung---
et,1,plane42
et,2,targe169,,1
et,3,conta171
keyopt,1,3,1

r,1
r,2
r,3,,,0.1

!---Geometrievereinbarung---

!keypoints
k,1,Di/2
k,3,(De-Di)/4+x1+Di/2
k,4,(De-Di)/4-x2+Di/2,b/2
k,5,Di/2,b/2
k,6,(De-Di)/4+x2+Di/2,se
k,7,De/2,se
k,8,De/2,E
k,9,De/2,se+b
k,10,(De-Di)/4+x2+Di/2,se+b
k,11,(De-Di)/4-x1+Di/2,E
k,13,Di/2,E+se
k,14,(De-Di)/4-x2+Di/2,E+se
k,15,(De-Di)/4+x1+Di/2,2*E
k,17,Di/2,2*E

!lines
l,1,3
l,3,4
l,4,5
l,5,1

l,6,7
l,7,8
l,8,11
l,11,6

l,8,9
l,9,10
l,10,11
l,11,8

l,14,15
l,15,17
l,17,13
l,13,14

al,1,2,3,4
al,5,6,7,8
al,7,9,10,11
al,12,13,14,15

!---Vernetzung---
type,1
mat,1
real,1
mshape,0
esize,,9
lesize,1,,,5
lesize,3,,,5
amesh,1
eplo

type,1
mat,2
real,1
mshape,0
esize,,9
lesize,7,,,5
lesize,5,,,5
amesh,2
eplo

type,1
mat,3
real,1
mshape,0
esize,,9
lesize,10,,,5
amesh,3
eplo

type,1
mat,4
real,1
mshape,0
esize,,9
lesize,13,,,5
lesize,15,,,5
amesh,4
eplo



!---Führungen
!n,10000,93/2
!n,10001,93/2,2*e
!n,10002,72.5/2
!n,10003,72.5/2,2*e
!fill,10000,10001,21,30000,1
!fill,10002,10003,21,20000,1




!---Kontakte erstellen---
!---L2/L8
lsel,s,,,2
nsll,s,1
type,2
mat,5
real,2
esurf
alls

lsel,s,,,8
nsll,s,1
type,3
esurf
alls

lsel,s,,,2
nsll,s,1
type,3
mat,5
real,3
esurf
alls

lsel,s,,,8
nsll,s,1
type,2
esurf
alls

!---L11/L12
lsel,s,,,11
nsll,s,1
type,2
mat,5
real,2
esurf
alls

lsel,s,,,12
nsll,s,1
type,3
esurf
alls

lsel,s,,,11
nsll,s,1
type,3
mat,5
real,3
esurf
alls

lsel,s,,,12
nsll,s,1
type,2
esurf
alls


!---------------SOLU----------------
/solu
!---nichtlineare Berechnung---
antype, static, new	!statische Analyse
eqslv, spar		!Specifies the type of equation solver
nlgeom, on		!nichtlineare Geometrie   
pred, -1		!Praediktor Method
			!chosen by Program for all DOFs

nsub, 40, 100, 10
outres, all, all

!---Randbedingungen---
!lsel, s, , ,1
!NSLL,S,1
!d, all, uy
!alls

NSEL,S,LOC,Y,0        ! Randbedingung unten
CP,1,UY,ALL            ! 
*GET,NC,NODE,,NUM,MIN 
d,NC,UY 

!nsel,s,loc,x,93/2
!d,all,all
!alls

!nsel,s,loc,x,72.5/2
!d,all,all
!alls

!------1.Lastfall - Ini-Lastfall
time, 1
!lsel, s, , , 13		!der kontakt wird initialisiert durch eine kleine Verschiebung
!NSLL,S,1
!d, all, uy, -0.05
!allsel
!lswrite,1

NSEL,S,LOC,Y,22        ! !!!!!
CP,2,UY,ALL            ! 
*GET,NL,NODE,,NUM,MIN 
d,NL,UY,-0.1
F,NL,FY,-300
allsel
lswrite,1 
!------2.Lastfall

time, 2	
!lsel, s, , , 13		!nach der ini wird die verschiebung gelöscht
!NSLL,S,1
!nsel,s,loc,y,2*E
!d,all,uy,-1
DDELE,NL,UY
!allsel
lswrite,2