!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,1,0.3

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

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

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

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

r,1
r,2

!---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, 10, 100, 10
outres, all, all

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

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

!------2.Lastfall

time, 2	
lsel, s, , , 13		!nach der ini wird die verschiebung gelöscht
NSLL,S,1
ddele, all, all
allsel

nsel,s,loc,y,2*E
d,all,uy,-1