*HEADING
 Full model: FAN=HUB+4xBLADES  
 Frequency analysis
**
*RESTART,WRITE,FREQUENCY=1
**
** NODE DEFINITION: HUB+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
**
*NSET,NSET=PART1,GENERATE
 1,120,1
*NCOPY,CHANGE NUMBER=1000,OLD SET=PART1,SHIFT
 0.,0.,0.
 3.,3.,-1.,3.,3.,1.,90.
*NCOPY,CHANGE NUMBER=2000,OLD SET=PART1,SHIFT
 0.,0.,0.
 3.,3.,-1.,3.,3.,1.,180.
*NCOPY,CHANGE NUMBER=3000,OLD SET=PART1,SHIFT
 0.,0.,0.
 3.,3.,-1.,3.,3.,1.,270.
**
** NODE SETS
**
*NSET,NSET=NHUBN,GENERATE
    1,   9, 1
 1002,1008, 1 
 2001,2009, 1 
 3002,3008, 1 
*NSET,NSET=NBLADE,GENERATE
   57,  61, 1
 1057,1061, 1
 2057,2061, 1
 3057,3061, 1
*NSET,NSET=PART1X,GENERATE
 1,82,9
*NSET,NSET=PART1Y,GENERATE
 9,90,9
*NSET,NSET=PART2X,GENERATE
 1009,1090,9
*NSET,NSET=PART2Y,GENERATE
 1001,1082,9
*NSET,NSET=PART3X,GENERATE
 2001,2082,9
*NSET,NSET=PART3Y,GENERATE
 2009,2090,9
*NSET,NSET=PART4X,GENERATE
 3009,3090,9
*NSET,NSET=PART4Y,GENERATE
 3001,3082,9
*MPC
 TIE,PART4X,PART1X
 TIE,PART2Y,PART1Y
 TIE,PART2X,PART3X
 TIE,PART4Y,PART3Y
*NSET,NSET=NXSYMM
 PART1Y,PART3Y
*NSET,NSET=NYSYMM
 PART1X,PART3X
**
** ELEMENT DEFINITION: HUB+4xBLADE
**
*ELEMENT,TYPE=S4R
 1,1,10,11,2
*ELGEN,ELSET=HUB1
 1,8,1,1,9,9,8
**
*ELEMENT,TYPE=S4R
 73,57,91,92,58
 77,91,96,97,92
*ELGEN,ELSET=BLADE1
 73,4,1,1
 77,4,1,1,5,5,4
**
*ELSET,ELSET=P1
 HUB1,BLADE1
*ELCOPY,ELEMENT SHIFT=100,OLD SET=P1,SHIFT NODES=1000,NEW SET=P2
*ELCOPY,ELEMENT SHIFT=200,OLD SET=P1,SHIFT NODES=2000,NEW SET=P3
*ELCOPY,ELEMENT SHIFT=300,OLD SET=P1,SHIFT NODES=3000,NEW SET=P4
*ELSET,ELSET=HUBF,GENERATE
   1,  24, 1
 101, 124, 1
 201, 224, 1
 301, 324, 1
*ELSET,ELSET=HUBT,GENERATE
  25,  72, 1
 125, 172, 1
 225, 272, 1
 325, 372, 1
*ELSET,ELSET=CYLBLADE,GENERATE
  73,  96, 1
 173, 196, 1
 273, 296, 1
 373, 396, 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
**
** STEP 1
**
*STEP,nlgeom,INC=1000
 Step 1: centrifugal loading
*STATIC
 .1,1.
**
*BOUNDARY
 NHUBN,ENCASTRE
*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 
*EL FILE,FREQ=0 
*NODE PRINT,FREQ=0 
*NODE FILE,FREQ=0 
**
*END STEP
**
** STEP 2
**
*STEP
 Step 2: Eigenfrequency extraction
*FREQUENCY, SIM=NO
 20,
**
*EL PRINT,FREQ=0 
*EL FILE,FREQ=0 
*NODE PRINT,FREQ=0 
*NODE FILE,FREQ=0 
*MODAL FILE,FREQ=99
*END STEP