*HEADING
 Substructure generations: FAN=HUB+4xBLADES  
 Five retained modes.
**
   Z200: a substructure which contains the 1/4 hub 
	 and a substructure blade (Z100);
**
*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
**
  118,   7.89643,   8.00249,   2.59810
 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
**
*NSET,NSET=NBLADE,GENERATE
 57,61,1
*NSET,NSET=NXSYMM,GENERATE
 9,90,9
*NSET,NSET=NYSYMM,GENERATE
 1,82,9
*NSET,NSET=RETAINED,UNSORTED
   1,   2,   3,   4,   5,   6,   7,   8,   9,  10,
  18,  19,  27,  28,  29,  30,  31,  32,  33,  34,
  35,  36,  37,  45,  46,  54,  55,  63,  64,  72,
  73,  81,  82,  90, 118
**
*ELEMENT,TYPE=Z100,ELSET=SBLADE,FILE=FAN
 800,57,58,59,60,61,118
*SUBSTRUCTURE PROPERTY,ELSET=SBLADE
**
*ELEMENT,TYPE=S4R
 1,1,10,11,2
*ELGEN,ELSET=HUB
 1,8,1,1,9,9,8
**
** The top hub and front hub have different thickness.
**
*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
*MATERIAL,NAME=MAT
*DENSITY
 7850.,
*ELASTIC,TYPE=ISOTROPIC
 2.E+11, 0.29
**
*NSET,NSET=NPRINT,GENERATE
 28,36,1
 55,63,1
**
** PRELOAD
**
*STEP
*STATIC
 .1,1.
*BOUNDARY
 NXSYMM,XSYMM
 NYSYMM,YSYMM
 NHUB,ENCASTRE
*SLOAD
 800,SCENTFGL,1.
*DLOAD
 HUB, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
*EL PRINT,FREQ=0 
*NODE PRINT,NSET=NPRINT,FREQ=100
 U,
 RF, 
*SUBSTRUCTURE PATH, ENTER ELEMENT=800
*NODE PRINT,NSET=NLOOK,FREQ=99
 U,
 RF,
*SUBSTRUCTURE PATH, LEAVE
*END STEP
**
*STEP
*FREQUENCY
 20,
*BOUNDARY
RETAINED, 1,6
*EL PRINT,FREQ=0 
*NODE PRINT,FRE=0
*END STEP
**
**
*STEP
*SUBSTRUCTURE GENERATE,TYPE=Z200,LIBRARY=FAN,RECOVERY MATRIX=YES,MASS MATRIX=YES,OVERWRITE
*RETAINED NODAL DOFS,SORTED=NO
RETAINED,
** *RETAINED EIGENMODES, GENERATE
** 1,20,1
**
*BOUNDARY,OP=NEW
*SUBSTRUCTURE LOAD CASE, NAME=SPRESS
*SLOAD
 800,SDLOAD,1.
**
*SUBSTRUCTURE LOAD CASE, NAME=SHUB
*DLOAD
 HUB, CENT, 6.3E6, 3., 3., 0., 0., 0., -1.
**
*END STEP