/FILNAME,shell,1

HOJDE = 1000
RADIUS = 250
TYKKELSE = 1
MESHSIZE1 = 20
MESHSIZE2 = 20

Lm = 4*sqrt(RADIUS*TYKKELSE)

/NOPR
KEYW,PR_STRUC,1

/NERR,0,0

/PREP7  
MPTEMP,,,,,,,,  		!
MPTEMP,1,0  			! materialedata
MPDATA,EX,1,,210000 		!
MPDATA,PRXY,1,,.3		!

TB,bkin,1,1,2
TBTEMP,0
TBDATA,,100,210

ET,1,SHELL93			! elementtype

R,1,TYKKELSE			! real constants

k,1,				!
k,2,,,HOJDE			!
k,3,RADIUS			!
k,4,RADIUS,,HOJDE		!
k,5,RADIUS-5,,HOJDE		!
K,6,RADIUS+5,,HOJDE		!

L,3,4

AROTAT,all,,,,,,1,2,90

LSEL,s,loc,z,1,HOJDE-1
LESIZE,all,MESHSIZE1
LSEL,inve
LESIZE,all,MESHSIZE2
ALLSEL
AMESH,all

ALLSEL

LSEL,s,loc,z,HOJDE-1,HOJDE+1
NSLL,s,1
NSLE,r,mid
*GET,MIDNODES,NODE,0,COUNT
F,all,FZ,(-6.283185307*RADIUS/MIDNODES*2/3)/4

CSYS,0
LSEL,s,loc,z,HOJDE-1,HOJDE+1
NSLL,s,1
NSLE,r,corner
*GET,ENDNODES,NODE,0,COUNT
NSEL,u,loc,y,-.1,.1
NSEL,u,loc,x,-.1,.1
F,all,FZ,(-6.283185307*RADIUS/ENDNODES*1/3)/4

NSEL,s,loc,y,-.1,.1
NSEL,a,loc,x,-.1,.1
NSEL,r,loc,z,HOJDE-1,HOJDE+1
F,all,FZ,(-6.283185307*RADIUS/ENDNODES*1/6)/4

ALLSEL

LSEL,s,loc,z,-1,1
DL,all,,ux
DL,all,,uy
DL,all,,uz
DL,all,,rotz

LSEL,s,loc,z,HOJDE-1,HOJDE+1
DL,all,,ux
DL,all,,uy
DL,all,,rotz

LSEL,s,loc,y,-.1,.1
LSEL,a,loc,x,-.1,.1
DL,all,,symm

ALLSEL

/solution
ANTYPE,0
PSTRES,ON   
SOLVE   
finish