/filnam, Drehfeder_Modell_1_6_BEAM4_SOLID45_CONTAC49
/title, Drehfeder-Modell mit BEAM4=Drehfeder  SOLID45=Bolzen	CONTAC49 

/prep7

/com Material 1, Material der Drehfeder
et,1,beam4,,				! Elementtyp BEAM4
mp,ex,1,206000				! E-modul 206000 N/mm²	
mp,nuxy,1,0.3				! Querkontraktionszahl 0.3
mp,gxy,1,81500				! G-Modul 81500 N/mm²

/com Material 2, Material des Bolzens
et,2,solid45
mp,ex,2,210000
mp,nuxy,2,0.3

/com Material 3, Kontaktmaterial
et,3,contac49
mp,ex,3,210000	
mp,nuxy,3,0.3
mp,mu,2,0.25
r,3,210000

/com Realkonstanten festlegen! r,nummer,Fläche,FlächenträgheitsmomentIZZ,Flächenträgheitsmomenr IYY,tky,tkz=Dicken nur 
d1=0.7					! Durchmesser 0.7 mm
pi=3.14159265359

r,1,d1**2*pi/4,d1**4*pi/64,d1**4*pi/64,d1,d1

/com Keypoints festlegen
k,1,0,0,0
k,2,0,0,-15
k,3,4,-0.2625,-15
k,4,4.06563,0.7375,-15
k,5,4,-0.7875,-15
k,6,3.94375,0.2125,-15
k,7,4,-1.3125,-15
k,8,4.065625,-0.3125,-15
k,9,4,-1.8375,-15
k,10,3.94375,-0.8375,-15
k,11,4,-2.3625,-15
k,12,4.065625,-1.3625,-15
k,13,4,-2.8875,-15
k,14,3.94375,-1.8875,-15

! letzter Achtelbogen
k,15,4,-3.4125,-15
k,16,4.065625,-2.4125,-15

/com Linien spannen
l,1,2 
lrotat,2,,,,,,4,3,180
lrotat,18,,,,,,6,5,180
lrotat,20,,,,,,8,7,180
lrotat,22,,,,,,10,9,180
lrotat,24,,,,,,12,11,180
lrotat,26,,,,,,14,13,180
lrotat,28,,,,,,16,15,45

! Jetzt existiert die Spur der Drehfeder
csys,0

local,12,1,1.129570,-3.018007,-17.89010,-90,45,90	!zylindrisches Koordinatensystem zur Generierung von Knoten30

csys,12

k,30,10

csys,0

l,29,30

/com Standardelementgrösse
esize,,6

lesize,1,,,15
lesize,15,,,10

type,1	$mat,1	$real,1 $lmesh,all

nsel,s,node,,13,141	$nsel,a,node,,2	$cm,Drehwindungen_surface,node 	!Knoten auf den Drehwindungen selektier

/com Volumenkörper erstellen

!lokales Koo-system erstellen (kartesisch,90 grad um die X-Achse gedreht)

!local,11,1,4,-4,-15,,-90
local,11,1,3.29289321881,-4,-15.707106781187,,-90
!erstelltes Koordinatensystem aktivieren

WPCSYS,,11

csys,11

k,31,0,0
k,32,3,0
k,33,3,90
k,34,3,180
k,35,3,270
k,36,0,0,6

a,32,33,34,35

l,31,36

vdrag,1,,,,,,20

lesize,22,,,12
lesize,23,,,12
lesize,25,,,12
lesize,27,,,12
lesize,20,,,12

type,2	$mat,2	$esize,,6	$mshape,0	$vmesh,1

/com Jetzt Kontaktelemente

asel,s,,,2,5	$nsla,,1	$cm,Bolzen_surface,node		! Knoten auf der bolzenaußenhaut selektieren und zusammenfassen

allsel


allsel

type,3	$mat,3	$real,3		$gcgen,Drehwindungen_surface,Bolzen_surface

allsel

wpCSYS,,0			!Koordinatensystem zurückschalten
csys,0



! Jetzt Koordinatensystem am Knoten 30 drehen, um Kraft richtig aufzubringen

!local,12,0,8.242640687,-3.28125,-15-9.899494937,,,135
wpcsys,,12
csys,12

nrotat,95
nrotat,104

d,95,uz,0
d,104,uz,0

!fk,30,fy,1.5
sfbeam,103,y,pres,-0.75

wpcsys,,0
csys,0

csys,11

! Bolzen lagern

nsel,s,loc,z,0

d,all,all,0

allsel

csys,0

d,1,ux,0
d,1,uy,0
d,1,uz,0
nsel,s,loc,y,-4,-2.5
nsel,r,node,,85,94
nsel,a,node,,96,103
d,all,uy,0



allsel

fini

/solu

! kopierter Teil für Nichtlineare Geometrie

$nlgeom,on $pred,on $autots,on	
time,1 $deltim,0.02,0.01,0.05