finish
/clear


!!======Preprocessing=======

/prep7
et,1,mesh200,7                   !Elementtyp1: 8 Knoten Flächenelement
et,2,solid186                    !Elementtyp2: 20 Knoten Volumenelement
mp,ex,1,210000
mp,nuxy,1,0.3
shpp,off
mshape,0
mshkey,2
type,1
*afun,rad
Pi=4*atan(1)

!*******Fläche*******
pcirc,11.5,12.75,0,360           !Innenring
pcirc,12.75,14,0,360             !Außenring
lesize,all,,,8
type,1
aglue,all
amesh,all

!*****Volumen*****
type,2
esize,,20
vext,all,,,,,100
aclear,all

!*****Pilotknoten*****
*get,nmax,node,,num,max
pilot1=nmax+1
pilot2=nmax+2


!-----1. Pilotknoten-----
nsel,s,loc,z,0

et,3,TARGE170                    !Zielelement
et,4,CONTA175                    !Kontaktelement

keyopt,3,2,1                     !boundary conditions for rigid target node specified by user
keyopt,3,4,111111                !all DOFs active
keyopt,4,12,5                    !bounded (always)
keyopt,4,4,1                     !normal from Contact node
keyopt,4,2,2                     !MPC (Multipoint constraint)

n,pilot1,0,0,0

r,3
real,3
type,3
tshap,pilo
e,pilot1

type,4
esln,s,0
esurf

!-----2. Pilotknoten-----
nsel,s,loc,z,100

et,5,TARGE170
et,6,CONTA175

keyopt,5,2,1
keyopt,5,4,111111
keyopt,6,12,5
keyopt,6,4,1
keyopt,6,2,2

n,pilot2,0,0,100

r,5
real,5
type,5
tshap,pilo
e,pilot2

type,6
esln,s,0
esurf
allsel


!*****Balkenanbindung*****
!Knotennummern
b1=nmax+3
b2=nmax+4
b3=nmax+5
b4=nmax+6
b5=nmax+7
b6=nmax+8
b7=nmax+9
b8=nmax+10


!Elementnummern
*get,emax,elem,,num,max
e1=emax+1
e2=emax+2
e3=emax+3
e4=emax+4
e5=emax+5
e6=emax+6
e7=emax+7
e8=emax+8

n,b1,0,0,-10
n,b2,0,0,-10
n,b3,5,5,-11
n,b4,10,10,-12
n,b5,10,10,-12
n,b6,10,10,-25
n,b7,0,-50,-25
n,b8,0,-50,-35

!*****CTube*****
et,7,beam188                     !Balkenelement
sectype,7,beam,ctube             !Typ: Rohr
secdata,30,50                    !Innen- und Außendurchmesser
secnum,7
type,7

en,e1,pilot1,b1
en,e2,b5,b6
en,e3,b6,b7
en,e4,b7,b8



!*****CSolid*****
et,8,beam188                     !Balkenelement
sectype,8,beam,csolid            !Typ: Stab
secdata,5,8,2                    !Radius und Teilungen
secnum,8
type,8

en,e5,b2,b3
en,e6,b3,b4

!-----Hilfkoordinatensystem-----
*get,kmax,kp,,num,max
*get,lmax,line,,num,max

k,kmax+1,0,0,-10
k,kmax+2,10,10,-12
l,kmax+1,kmax+2

lwplan,11,lmax+1,0.5             !WP in Linienmitte legen
csys,4
clocal,13,0,0,0,0,0,0,0

ldele,lmax+1
kdele,kmax+1,kmax+2

nrotat,b3                        !Knoten b3 in KOOS 13 rotieren

!*****Kugellager*****
et,9,MPC184                      !MPC184
keyopt,9,1,15                    !Spherical Joint Element
sectype,9,joint,sphe,sphball
secnum,9
type,9

!-----Rohrseite-----
local,11,0,0,0,-10,0,0,0
secjoint,lsys,11                 !Define local coordinate systems at joint element
en,e7,b1,b2

!-----andere Seite-----
local,12,0,10,10,-12,0,0,0
secjoint,lsys,12                 !Define local coordinate systems at joint element
en,e8,b4,b5





!!=====Solution=====
/solu
d,pilot2,ux,0,,,,uy,uz,rotx,roty,rotz        !feste Einspannung

d,b3,rotz,0                                  !Rotation des Links sperren

d,b7,ux,0,,,,uy,uz,rotx,roty                 !Lagerung der Welle
d,b8,ux,0,,,,uy,uz,rotx,roty

d,b8,rotz,10*(pi/180)                        !Auslenkung um 10°

allsel

save
solve
save