!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=2
adiv=3.2
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

!---Spiegeln
asel,all
arsym,y,all		!Spiegeln -> 4 Elemente
allsel

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

!---L27/L24
lsel,s,,,27
nsll,s,1
type,2
mat,5
real,3
esurf
alls

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

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

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

!---L23/L17
lsel,s,,,23
nsll,s,1
type,2
mat,5
real,3
esurf
alls

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

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

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

!---L16/L1
lsel,s,,,16
nsll,s,1
type,2
mat,5
real,3
esurf
alls

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

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

lsel,s,,,1
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,-2*E        ! 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

!-------------Postprocesing------------------
/post26

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

abs,2,2,,,,,,-1
add,3,3

xvar,2
plvar,3

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

plvar,3