/clear,nostart 									! löscht die Anzeige

! allgemeine Eingaben
/filename,aufgabe-b 							! Dateiname
/titel,Drehfeder oben eingespannt (siehe Aufgabe Abbildung b)             		  

/eshape,0.25                      !überhöhte darstellung

!##############################################
!Variablen Block
N         = 4         !windungszahl
di        = 5         !Durchmesser mittellinie
steig     = 2         !steigung
verschieb = 30        !des neuen punktes y-achse
teil      = 360/verschieb
H         = 5         ! länge des steges
punkt1    = 4         ! angriffspunkt 1 der kraft
punkt2    = 1         ! angriffspunkt 2 der kraft
!##############################################

/prep7
csys,1          !zylindrisches KOS

!keypoints erstellen
!keypoint,nr,x,y,z
k,1,di/2,0,0
!anzahl der widerholungen,keypoint1,x-verschiebung,y-verschiebung(abstand zwischen zwei keypoints),
!z-verschiebung(anstieg)
*repeat,(teil*N)+1,1,0,verschieb,steig/teil,
max_kp = (teil*N)+1



et,1,BEAM4
r,1,2,1,1,1,1                                   !real constants
keyopt,1,2,0
keyopt,1,9,0
mp,ex,1,300000
!mp,GXY,1,0.3						
!mp,DAMP,1,2

PSTRES,ON !???

!linien für die windung erstellen erstellen
l,1,2
*repeat,(teil*N),1,1

!alle linien slektieren
lsel,all
!alle zu einer linie zusammenfassen
lcomb,all

!****************************************
!  unterer steg
!****************************************
csys,0      !kart.KOS
clocal,11,0,di/2,0,0 !kos in den ersten punkt legen
k,,H,0,0
l,1,2


!linien am 1 punkt teilen und den ersten punkt löschen um einen radiuseinzufügen
!ldiv verwenden um linien zu teilen
!mit ratio verschieben
!spline anstiege ausrechnen
ldiv,1,0.005
ldiv,2,0.08
ldele,1,2,1   !löscht die linien
csys,0

!anstiegsberechnung
ldiv,3,0.0001
x_buf=0
y_buf=0
z_buf=0

x_dif=0
y_dif=0
z_dif=0

*get,x_buf,kp,5,loc,x
*get,y_buf,kp,5,loc,y
*get,z_buf,kp,5,loc,z

*get,x_dif,kp,3,loc,x
*get,y_dif,kp,3,loc,y
*get,z_dif,kp,3,loc,z

x_dif = x_buf - x_dif
y_dif = y_buf - y_dif
z_dif = z_buf - z_dif


spline,3,,,,,4,x_dif,y_dif,z_dif,1,0,0


!****************************************
!  oberer steg
!****************************************
csys,0      !kart.KOS
!koordinaten von punkt aufnehmen
x_kp=0
y_kp=0
z_kp=0

*get,x_kp,kp,max_kp,loc,x
*get,y_kp,kp,max_kp,loc,y
*get,z_kp,kp,max_kp,loc,z

!mittelpunkt erstellen
k,3000,0,0,z_kp-0.00001
k,2999,0,0,z_kp
cskp,11,0,max_kp,3000,2999

k,,-H,0,0
l,max_kp,6

!linien am 1 punkt teilen und den ersten punkt löschen um einen radiuseinzufügen
!ldiv verwenden um linien zu teilen
!mit ratio verschieben
!spline anstiege ausrechnen
ldiv,1,0.995
ldiv,5,0.08
ldele,5,6,1   !löscht die linien..... linie1,linie2,inkrement

csys,0


!anstiegsberechnung der windung
ldiv,1,0.9999
x_buf=0
y_buf=0
z_buf=0

x_dif=0
y_dif=0
z_dif=0

*get,x_buf,kp,7,loc,x
*get,y_buf,kp,7,loc,y
*get,z_buf,kp,7,loc,z

*get,x_dif,kp,9,loc,x
*get,y_dif,kp,9,loc,y
*get,z_dif,kp,9,loc,z

x_dif = x_buf - x_dif
y_dif = y_buf - y_dif
z_dif = z_buf - z_dif

!anstiegsberechnung der abzweigung
ldiv,7,0.00001
x2_buf=0
y2_buf=0
z2_buf=0

x2_dif=0
y2_dif=0
z2_dif=0

*get,x2_buf,kp,8,loc,x
*get,y2_buf,kp,8,loc,y
*get,z2_buf,kp,8,loc,z

*get,x2_dif,kp,10,loc,x
*get,y2_dif,kp,10,loc,y
*get,z2_dif,kp,10,loc,z

x2_dif = x2_buf - x2_dif
y2_dif = y2_buf - y2_dif
z2_dif = z2_buf - z2_dif

spline,8,,,,,7,-x2_dif,-y2_dif,-z2_dif,-x_dif,-y_dif,-z_dif

!****************************************

!angriffslinien teilen
ldiv,4,punkt1/H,
H_stern = H - (H-punkt1)
ldiv,4,punkt2/H_stern,


*get,x_buf,kp,12,loc,x
*get,y_buf,kp,12,loc,y
*get,z_buf,kp,12,loc,z

*get,x2_buf,kp,11,loc,x
*get,y2_buf,kp,11,loc,y
*get,z2_buf,kp,11,loc,z

LCOMB,1,3
LCOMB,1,5
LCOMB,1,8
LCOMB,1,2
LCOMB,1,7

nkpt,,2
nkpt,,11
nkpt,,12
nkpt,,4

e,1,2
e,2,3
e,3,4

ESIZE,0.10,0,   
LMESH,1

!Oberensteg meshen
lmesh,6

!oberesteg-knoten selektieren
nsel,s,loc,x,2.9,N*steig,0.001
nsel,r,loc,z,N*steig

!knoten DOF festsetzen
d,all,all

!selektion aufheben
allsel,all


!zwei knoten müssen verbunden werden
nsel,s,node,,4,5,1
nummrg,node

allsel,all

! unten kraft anbringen
f,2,fy,-500
f,3,fy,-500


!Start des Solutionprozessors
/solu
antype,Static
solve													! beginne mit berechnung
finish


! Start des Postprozessors
/post1
pldisp,1

!Nodal Solution X-Stress
prnsol,s,comp