!PREPROCESSING

/filnam,Ringfeder

!---Geometriedaten---
!Werte für Fmax=100kN
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


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

NSEL,S,LOC,Y,0        		! Randbedingung unten
d,all,UY,0 
allsel

NSEL,S,LOC,Y,2*E        
d,all,UY,-1
allsel


nlgeom,on			!nichtlineare geometrie

autots,on			!Last- oder Zeitschrittautomatik
nsub,100,1000,30

outrse,all,all			!Ergebnis jeden Zwischenschritts abspeichern

pred,on				!Folge der Iterationen schneller auf Konvergenz hinführen

save

solve


!---------------------------POSTPROCESSING----------------------------
/post26				!Verlauf einzelner Ergebnisgrößen über Folge der Substeps auswerten

nn=node(Di/2,2*E,0)		!Knoten oben innen wählen

nsol,2,nn,u,y			!Tabelle mit Y-Verschiebung des gewünschten Knotens
rfor,3,nn,f,y			!Reaktionskraft in Y-Richtung in Tabelle einbinden

abs,2,2
add,3,3,,,,,,-1*12

xvar,2
plvar,3

/grid,1				!Diagramm beschriften
/axla,x,Federweg s (mm)
/axla,y,Kraft F (N)
plvar,3