/solution
ANTYPE,1
BUCOPT,LANB,1,0
SOLVE  
finish
/solution
EXPASS,on
OUTRES,all,all
finish

/post1
CSYS,1
SET,FIRST
*GET,PBUCKL,active,0,set,freq	! første bulingslast
NSORT,u,x,0,1
*GET,maxudboj,sort,,max		! maksimal udbøjning fra første mode
*GET,mnodemax,sort,,imax	! knude, hvor maksimal udbøjning er
NSEL,s,loc,z,HOJDE-.1,HOJDE+.1
NSORT,u,z,0,1
*GET,tnodemax,sort,,imax	! knude, hvor maksimal flytning i toppen er

ALLSEL

/prep7

*IF,1.5*Lm/100,gt,1.5*TYKKELSE,then
w0 = 1.5*Lm/100
*ELSE
w0 = 1.5*TYKKELSE
*ENDIF
udboj = w0/maxudboj

*STATUS,PARM


UPGEOM,w0/maxudboj,LAST,LAST,shell,rst,
FDELE,all,all

Q = PBUCKL

LSEL,s,loc,z,HOJDE-1,HOJDE+1
NSLL,s,1
NSLE,r,mid
*GET,MIDNODES,NODE,0,COUNT
F,all,FZ,Q*(-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,Q*(-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,Q*(-6.283185307*RADIUS/ENDNODES*1/6)/4

ALLSEL

CSYS,1
/solution
ANTYPE,0
NROPT,full
NLGEOM,ON
NSUBST,50,0,0
OUTRES,ALL,ALL
ARCLEN,1,.1,0.001
ARCTRM,L
AUTOTS,-1.0
solve

save