!/CLEAR

/UNITS,SI

/PREP7

!! Elementtyp auswählen
 
ET,1,PIPE16

!! Materialeigenschaften zuordnen

MP,EX,1,2.1E11		!Elastizitätsmodul in x-Richtung [N/m^2]
MP,NUXY,1,0.3		!Querkontraktionszahl xy-Richtung
MP,DENS,1,7.800E-3	!Dichte [kg/m^3]


!! Erstellung der Geometrie

! Definition der geometrischen Parameter [m]

dw=10			!Durchmesser der Welle
lw=1000			!Laenge der Welle

R,1,dw,dw/2

! Definiere Knoten

N,1,0
N,9,lw

FILL,1,9


!! Vernetzung

ESIZE,5
TYPE,1
REAL,1
E,1,2
E,2,3
E,3,4
E,4,5
E,5,6
E,6,7
E,7,8
E,8,9


!! Loesungsabschnitt

/SOLU
ANTY,TRANS

TRNOPT,FULL

NSEL,S,LOC,X,0
D,ALL,UX
D,ALL,UY
D,ALL,UZ			!Knoten am linken Rand (0; Knoten 1) auswaehlen für feste Lagerung


ALLS


NSEL,S,LOC,X,lw			!Knoten am rechten Rand (lw; Knoten 9) fuer Loslagerung
D,ALL,UY
D,ALL,UZ


ALLS

NSEL,S,LOC,X,lw/2

TIME,1

DELTIM,1

TIMINT,OFF

F,ALL,FY,-100			!Kraft auf Wellenmitte (lw/2; Knoten 5) bis 1s

ALLS


SOLVE

TIMINT,ON

TIME,1.5

AUTOTS,ON

DELTIM,0.005,0.001,0.01

OUTRES,ALL,ALL

FDEL,ALL,FY			!Entfernen der Kraft bis 1.5s

ALLSEL

SOLVE

/POST26

FILE,,RDSP

NSOL,2,4,U,Y

/AXLAB,Y,Verschiebung

PLVAR,2				!Verschiebung-Zeit von Knoten 4