!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
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---
ldiv=8
adiv=0.4
type,1
mat,1
real,1
mshkey,2
esize,adiv
lesize,1,,,ldiv
lesize,3,,,ldiv
amesh,1
eplo

type,1
mat,2
real,3
mshkey,2
esize,adiv
lesize,7,,,ldiv
lesize,5,,,ldiv
amesh,2
eplo

type,1
mat,3
real,3
mshkey,2
esize,adiv
lesize,10,,,ldiv
amesh,3
eplo

type,1
mat,4
real,3
mshkey,2
esize,adiv
lesize,13,,,ldiv
lesize,15,,,ldiv
amesh,4
eplo


!---Kontakte erstellen---
!---L2/L8
lsel,s,,,2
nsll,s,1
type,2
mat,5
real,3
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,3
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


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


!---Lastfall---
NSEL,S,LOC,Y,22        ! !!!!!
CP,2,UY,ALL            ! 
*GET,NL,NODE,,NUM,MIN

*do,k,1,10
time,k
d,NL,UY,-k*0.1
LSWRITE,k
*enddo

*do,k,11,20
time,k
d,NL,UY,-1+(k-10)*0.1
LSWRITE,k
*enddo


LSSOLVE,1,20


!POSTPROCESSING

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

!nn=node(Di/2,2*E,0)

!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	

!xvar,2
!plvar,3

!/grid,1				!Diagramm beschriften
!/axla,x,Weg s (mm)
!/axla,y,Kraft F (N)
!/xran,0,12
!/yran,0,21000

!plvar,3