*HEADING
 Usage level: FAN=HUB+4xBLADES
 Substructure analysis of a fan
 Static analysis : Two level substructure.
**
   Z100: a substructure for a single blade;
   Z200: a substructure which contains the 1/4 hub 
         and a substructure blade (Z100);
   Z300: a substructure which contains the 1/4 hub 
         and a full blade.
Requires substructure generation file fansubstr_gen2.inp
**
*RESTART,WRITE,FRE=1
**
*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
**  
   31,   6.00000,   3.00000,   4.00000  
   33,   5.79830,   4.08144,   4.00000
   37,   4.08144,   5.79830,   4.00000   
   39,   3.00000,   6.00000,   4.00000
**
   91,   6.00000,   3.00000,   0.00000   
   93,   5.79830,   4.08144,   0.00000
   97,   4.08144,   5.79830,   0.00000 
   99,   3.00000,   6.00000,   0.00000
**
  100,   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
 31,33,1,9998
 33,37,1,9998
 37,39,1,9998
*NGEN,LINE=C,NSET=RINGB
 91,93,1,9999
 93,97,1,9999
 97,99,1,9999
*NFILL,NSET=HUB
 NHUB,RINGF,3,10
 RINGF,RINGB,6,10
**
*NSET,NSET=PART1,GENERATE
 1,91,10
 9,99,10
 2,8,1
 32,38,1
 100,100,1
**
*NCOPY,CHANGE NUMBER=100,OLD SET=PART1,SHIFT
 0.,0.,0.
 3.,3.,-1.,3.,3.,1.,90.
*NCOPY,CHANGE NUMBER=200,OLD SET=PART1,SHIFT
 0.,0.,0.
 3.,3.,-1.,3.,3.,1.,180.
*NCOPY,CHANGE NUMBER=300,OLD SET=PART1,SHIFT
 0.,0.,0.
 3.,3.,-1.,3.,3.,1.,270.
**
*NSET,NSET=PART1X,GENERATE
 1,91,10
*NSET,NSET=PART1Y,GENERATE
 9,99,10
*NSET,NSET=PART2X,GENERATE
 109,199,10
*NSET,NSET=PART2Y,GENERATE
 101,191,10
*NSET,NSET=PART3X,GENERATE
 201,291,10
*NSET,NSET=PART3Y,GENERATE
 209,299,10
*NSET,NSET=PART4X,GENERATE
 309,399,10
*NSET,NSET=PART4Y,GENERATE
 301,391,10
*MPC
 TIE,PART4X,PART1X
 TIE,PART2Y,PART1Y
 TIE,PART2X,PART3X
 TIE,PART4Y,PART3Y
**
*NSET,NSET=NODEHUB,GENERATE
   1,   9, 1
 102, 108, 1
 201, 209, 1
 302, 308, 1
**
*ELEMENT,TYPE=Z200,ELSET=P1,FILE=FAN
 901,1,2,3,4,5,6,7,8,9,11,19,21,29,31,32,
 33,34,35,36,37,38,39,41,49,51,59,61,69,71,79,81,
 89,91,99,100
*ELCOPY,ELEMENT SHIFT=1,OLD SET=P1,SHIFT NODES=100,NEW SET=P2
*ELCOPY,ELEMENT SHIFT=2,OLD SET=P1,SHIFT NODES=200,NEW SET=P3
*ELCOPY,ELEMENT SHIFT=3,OLD SET=P1,SHIFT NODES=300,NEW SET=P4
**
*SUBSTRUCTURE PROPERTY,ELSET=P1
 0.,0.,0.
*SUBSTRUCTURE PROPERTY,ELSET=P2
 0.,0.,0.
 3.,3.,0.,3.,3.,3.,90.
*SUBSTRUCTURE PROPERTY,ELSET=P3
 0.,0.,0.
 3.,3.,0.,3.,3.,3.,180.
*SUBSTRUCTURE PROPERTY,ELSET=P4
 0.,0.,0.
 3.,3.,0.,3.,3.,3.,270.
**
*NSET,NSET=NPRINT,GENERATE
 31,39,1
 61,69,1
**
*STEP
 Pressure loading on the blade
*STATIC
 .1,1.
*BOUNDARY
 NODEHUB,ENCASTRE
*SLOAD
 P1, SPRESS, 1.
 P2, SPRESS, 1.
 P3, SPRESS, 1.
 P4, SPRESS, 1.
**
*EL PRINT,FREQ=0
*EL FILE,FREQ=0 
*NODE PRINT,NSET=NPRINT,FREQ=100
 U,
 RF, 
*NODE FILE,FREQ=0 
**
*OUTPUT,VAR=PRESELECT,FIELD
*SUBSTRUCTURE PATH, ENTER ELEMENT=901
*OUTPUT,FIELD,VAR=PRESELECT
*EL PRINT,FREQ=0
*NODE PRINT,NSET=NPRINT,FRE=100
 U,
 RF,
*NODE FILE,NSET=NPRINT,FRE=100
 U,
 RF,
*SUBSTRUCTURE PATH, ENTER ELEMENT=800
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT
 U,
*EL PRINT,FREQ=0
**NODE PRINT,FRE=0
*NODE PRINT,NSET=NLOOK,FRE=100
 U,
 RF,
*SUBSTRUCTURE PATH, LEAVE
*SUBSTRUCTURE PATH, LEAVE
**
*END STEP