*HEADING
 Quarter model: FAN=HUB/4+BLADE    fan_cyclicsymmodel
 Cyclic symmetry model option and steady state dynamics
**
*RESTART,WRITE,FREQUENCY=1
**
** NODE DEFINITION: HUB/4+BLADE
** Nodes for the hub: 1--90 
**
*NODE
    1,   3.50000,   3.00000,   4.00000   
    3,   3.46638,   3.18024,   4.00000
    7,   3.18024,   3.46638,   4.00000
    9,   3.00000,   3.50000,   4.00000
**  
   28,   6.00000,   3.00000,   4.00000
   30,   5.79830,   4.08144,   4.00000
   34,   4.08144,   5.79830,   4.00000   
   36,   3.00000,   6.00000,   4.00000
**
   82,   6.00000,   3.00000,   0.00000
   84,   5.79830,   4.08144,   0.00000
   88,   4.08144,   5.79830,   0.00000         
   90,   3.00000,   6.00000,   0.00000
**
 9998,   3.00000,   3.00000,   4.00000
 9999,   3.00000,   3.00000,   0.00000
*NGEN,LINE=C,NSET=NHUB
 1,3,1,9998
 3,7,1,9998
 7,9,1,9998
*NGEN,LINE=C,NSET=RINGF
 28,30,1,9998
 30,34,1,9998
 34,36,1,9998
*NGEN,LINE=C,NSET=RINGB
 82,84,1,9999
 84,88,1,9999
 88,90,1,9999
*NFILL,NSET=HUB
 NHUB,RINGF,3,9
 RINGF,RINGB,6,9
**
** Total nodes for blade: 91--120, 
**     plus common nodes  57--61
**
*NODE,NSET=BLADE
   91,   6.34825,   4.61208,   2.00000
   92,   6.04542,   5.17363,   2.04275 
   93,   5.64799,   5.65929,   2.11107
   94,   5.16740,   6.05420,   2.19011
   95,   4.61648,   6.34954,   2.25973
   96,   6.97098,   5.04709,   2.00000
   97,   6.60952,   5.64419,   2.08975
   98,   6.14537,   6.18965,   2.22829
   99,   5.61424,   6.64337,   2.39622
  100,   5.04534,   6.97776,   2.55559
  101,   7.71735,   5.22982,   2.00000
  102,   7.23607,   5.92235,   2.13475
  103,   6.57818,   6.65456,   2.33376
  104,   5.87275,   7.29336,   2.58456
  105,   5.23136,   7.72660,   2.83877
  106,   8.44058,   5.48119,   2.00000
  107,   7.86154,   6.25188,   2.17913
  108,   7.04773,   7.13818,   2.43316
  109,   6.20063,   7.93303,   2.76013
  110,   5.49725,   8.45477,   3.10513
  111,   8.94199,   6.04637,   2.00000
  112,   8.35346,   6.80121,   2.24584
  113,   7.55742,   7.64823,   2.59979
  114,   6.73147,   8.41155,   3.04013
  115,   6.03050,   8.93642,   3.47975
  116,   8.88910,   6.78398,   2.00000
  117,   8.41737,   7.43225,   2.26100
  118,   7.89643,   8.00249,   2.59810
  119,   7.35361,   8.48209,   3.02957
  120,   6.78399,   8.88910,   3.50000
**
** NODE SETS
**
*NSET,NSET=XAXIS,GENERATE
 10,82,9
*NSET,NSET=YAXIS,GENERATE
 18,90,9
*surface,type=node,name=secondary1
YAXIS,1.5
*surface,type=node,name=main1
XAXIS,1.5
*tie,cyclicsym,position=0.5,name=tie
 secondary1,main1
*cyclic symmetry model,n=4
3., 3., 0.,  3., 3., 10.
*ELEMENT,TYPE=S4R
 1,1,10,11,2
*ELGEN,ELSET=HUB
 1,8,1,1,9,9,8
*ELEMENT,TYPE=S4R
 73,57,91,92,58
 77,91,96,97,92
*ELGEN,ELSET=CYLBLADE
 73,4,1,1
 77,4,1,1,5,5,4
**
*ELSET,ELSET=HUBF,GENERATE
   1,  24, 1
*ELSET,ELSET=HUBT,GENERATE
  25,  72, 1
**
*SHELL SECTION,ELSET=HUBF,MATERIAL=MAT
 0.250, 3
*SHELL SECTION,ELSET=HUBT,MATERIAL=MAT
 0.125, 3
*SHELL SECTION,ELSET=CYLBLADE,MATERIAL=MAT
 0.125, 3
*MATERIAL,NAME=MAT
*DENSITY
 7850.,
*ELASTIC,TYPE=ISOTROPIC
 2.E+11, 0.29
*NSET,NSET=NPRINT,GENERATE
 28,36,1
 55,63,1
**
*NSET,NSET=PRINT,GEN
116,120,1
*BOUNDARY
 NHUB,ENCASTRE
**
*STEP,nlgeom,INC=1000
 Step 1: centrifugal loading
*STATIC
 .1,1.
**
*DLOAD
 CYLBLADE, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
 HUBF, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
 HUBT, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
**
*EL PRINT,FREQ=0
*NODE FILE,FREQ=2,NSET=PRINT
U
*NODE PRINT,FREQ=0 
*NODE FILE,FREQ=0 
**
*END STEP
**
*STEP
** *FREQUENCY,eigensolver=lanczos
*FREQUENCY,eigensolver=lanczos,normalization=mass
 26,,,,8
*SELECT CYCLIC SYMMETRY MODES, NMIN=0, NMAX=0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*END STEP
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
6.0,137.0,10,4.
*DLOAD, LOAD CASE=1, OP=NEW, CYCLIC MODE=0
 CYLBLADE, P, 6.3e6
 HUBF, P, 6.3e6
 HUBT, P, 6.3e6
*DLOAD, LOAD CASE=2, OP=NEW, CYCLIC MODE=0
 CYLBLADE, P, 6.3e6
 HUBF, P, 6.3e6
 HUBT, P, 6.3e6
*MODAL DAMPING,MODAL=DIRECT
 1,26,0.03
*SELECT EIGENMODES, GENERATE
1,26,1
*EL PRINT,ELSET=CYLBLADE,FREQUENCY=10
S11,E11
*NODE PRINT,NSET=PRINT,FREQUENCY=10
U, 
*NODE FILE,NSET=PRINT,FREQUENCY=10
U,V,A,RF
*MODAL PRINT,FREQUENCY=10
GU, 
GA, 
GPU, 
*EL FILE,ELSET=CYLBLADE,FREQUENCY=10
S, 
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S, 
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S, 
*END STEP
*STEP
** *FREQUENCY,eigensolver=lanczos
*FREQUENCY,eigensolver=lanczos,normalization=mass
 56,,,,8
*SELECT CYCLIC SYMMETRY MODES, NMIN=1, NMAX=1
*OUTPUT, FIELD, FREQUENCY=0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*NODE FILE,FREQ=0
*EL FILE,FREQ=0
*MODAL FILE
**
*END STEP
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
6.0,137.0,10,4.
*DLOAD, LOAD CASE=1, OP=NEW, CYCLIC MODE=1
 CYLBLADE, P, 6.3e6
 HUBF, P, 6.3e6
 HUBT, P, 6.3e6
*DLOAD, LOAD CASE=2, OP=NEW, CYCLIC MODE=1
 CYLBLADE, P, 6.3e6
 HUBF, P, 6.3e6
 HUBT, P, 6.3e6
*MODAL DAMPING,MODAL=DIRECT
 1,26,0.03
*SELECT EIGENMODES, GENERATE
1,26,1
*EL PRINT,ELSET=CYLBLADE,FREQUENCY=10
S11,E11
*NODE PRINT,NSET=PRINT,FREQUENCY=10
U, 
*NODE FILE,NSET=PRINT,FREQUENCY=10
U,V,A,RF
*MODAL PRINT,FREQUENCY=10
GU, 
GA, 
GPU, 
*EL FILE,ELSET=CYLBLADE,FREQUENCY=10
S, 
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S, 
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S, 
*END STEP
*STEP
** *FREQUENCY,eigensolver=lanczos
*FREQUENCY,eigensolver=lanczos,normalization=mass
 56,,,,8
*SELECT CYCLIC SYMMETRY MODES, NMIN=2, NMAX=2
*OUTPUT, FIELD, FREQUENCY=0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*EL FILE,FREQ=0
*NODE FILE,FREQ=0
*MODAL FILE
**
*END STEP
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
6.0,137.0,10,4.
*DLOAD, LOAD CASE=1, OP=NEW, CYCLIC MODE=2
 CYLBLADE, P, 6.3e6
 HUBF, P, 6.3e6
 HUBT, P, 6.3e6
*DLOAD, LOAD CASE=2, OP=NEW, CYCLIC MODE=2
 CYLBLADE, P, 6.3e6
 HUBF, P, 6.3e6
 HUBT, P, 6.3e6
*MODAL DAMPING,MODAL=DIRECT
 1,26,0.03
*SELECT EIGENMODES, GENERATE
 1,26,1
*EL PRINT,ELSET=CYLBLADE,FREQUENCY=10
S11,E11
*NODE PRINT,NSET=PRINT,FREQUENCY=10
U, 
*NODE FILE,NSET=PRINT,FREQUENCY=10
U,V,A,RF
*MODAL PRINT,FREQUENCY=10
GU, 
GA, 
GPU, 
*EL FILE,ELSET=CYLBLADE,FREQUENCY=10
S, 
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S, 
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S, 
*END STEP