FINI
/CLEAR
*set,rw,78     !Radius Welle
*set,irh,77.73 !Innenradius Hülse
*set,arh,85    !Außenradius Hülse
*set,aw,60     !axiale Länge Welle
*set,ah,60     !axiale Länge Hülse
!
/PREP7
CSYS,5
ET,1,PLANE42
ET,2,SOLID45
ET,3,CONTA178
KBC,0         !Rampenweise Anstieg
KEYOPT,3,2,1  !Penalty-Verfahren
KEYOPT,3,5,1  !Richtung Kontaktnormale 
KEYOPT,3,9,1  !Rampenweises Schliessen des Gaps.
!
MP,EX,1,2.1e5  !Material
MP,NUXY,1,.3
MP,MU,1,.2   
!Globales KS im Mittelpunkt bei halber axialer Länge
RECT,0,rw,-aw/2,aw/2 !Geo Welle 
RECT,irh,arh,-ah/2,ah/2 !Geo Hülse
!
!Netz
!Anzahl der Elemente
TYPE,1
MAT,1
LESIZE,1,,,16    !radiale Elementeinteilung Welle 
LESIZE,3,,,16
LESIZE,5,,,2     !radiale Elementeinteilung Hülse
LESIZE,7,,,2

LESIZE,2,,,12    !axiale  Elementeinteilung Welle 
LESIZE,4,,,12
LESIZE,6,,,12    !axiale  Elementeinteilung Hülse
LESIZE,8,,,12
!
!Netz Welle
AMESH,1
TYPE,2
ESIZE,,10 !10 Unterteilung pro Erstellten VOlumen
VROTAT,1,,,,,,1,4,360,4
!
!Netz Hülse
!*GET,NNN,NODE,,NUM,MAX !Max. Knotennummer
NUMSTR,NODE,100001
TYPE,1
AMESH,2
!Rotieren der Hülse
TYPE,2
VROTAT,2,,,,,,1,4,360,4
!
ACLEAR,1,2
EPLOT
!
CSYS,5
!NSEL,S,LOC,X,irh,rw
NROTAT,ALL
TYPE,3
EINTF,,,Low,5,0.27
!
ANTYPE,STATIC
/NERR,100,100
/SOLU
LSCLEAR,ALL
NROPT,UNSYM
csys,5
nsel,s,loc,x,0
csys,0
d,all,ux
d,all,uy
d,all,uz
alls
!
time,100
autots,off
nsubst,1
lswrite
solve