*HEADING
 COMPLEX EIGENSOLVER - SUBSTRUCTURE TEST
SUBSTRUCTURE GENERATION FILE pcfreq_sup_gen.inp SHOULD BE RUN FIRST.
*NODE, NSET=NROTOR
*ELEMENT,TYPE=Z1, ELSET=EROTOR, FILE=pcfreq_sup_gen
*SUBSTRUCTURE PROPERTY,ELSET=EROTOR
*NSET ,NSET= NSET_SIDE1
*NSET ,NSET= NSET_SIDE2
*SURFACE, TYPE=NODE, NAME= SURF_SIDE1
  NSET_SIDE1,
*SURFACE, TYPE=NODE, NAME= SURF_SIDE2
  NSET_SIDE2,
*MATERIAL  ,NAME=A1
*ELASTIC
 200, .3
*DENSITY
 1.,
*ORIENTATION      ,SYS=CYLI,NAME= ORAXI
*NODE, NSET=NPAD1
*ELEMENT, TYPE=M3D4, ELSET=ELPAD1
*NODE, NSET=NPAD2
*ELEMENT, TYPE=M3D4, ELSET=ELPAD2
*MATERIAL, NAME=A2
*ELASTIC
 1000, .3
*DENSITY
 1.,
*MEMBRANESECTION, ELSET=ELPAD1, MATERIAL=A2
*MEMBRANESECTION, ELSET=ELPAD2, MATERIAL=A2
*SURFACE, NAME=SURF_PAD1
 ELPAD1, SPOS
*SURFACE, NAME=SURF_PAD2
 ELPAD2, SNEG
*CONTACT PAIR, INTERACTION=ROTOR_PAD, SMALL SLIDING, NOTHICK
  SURF_SIDE1, SURF_PAD1
  SURF_SIDE2, SURF_PAD2
*SURFACE INTERACTION, NAME=ROTOR_PAD
*FRICTION
 0.0,
*BOUNDARY
 NPAD1, 1,2
 NPAD2, 1,2
*****************************
*CONSTRAINT CONTROL, no checks
*STEP, INC=100, NLGEOM
 ESTABLISH CONTACT
*STATIC
 1., 1.0
*BOUNDARY
 NPAD1, 3,3,  0.05
 NPAD2, 3,3, -0.05
*OUTPUT, FIELD, FREQ=10
*NODE OUTPUT
 U,
*END STEP
****
*STEP, INC=100, NLGEOM
*STATIC
 .5, 1.0
*CHANGE FRICTION, INTERACTION=ROTOR_PAD
*FRICTION
 0.3,
*MOTION, ROTATION
 NROTOR, 2.5, .0,.0,.0, .0,.0,1.
*OUTPUT, FIELD, FREQ=10
*NODE OUTPUT
 U,
*END STEP
**
*STEP
 REAL EIGENSOLVER
*FREQUENCY
 12,
*SUBSTRUCTURE PATH,ENTER=1
*NODE PRINT
 U,
*OUTPUT,FIELD
*NODE OUTPUT
 U,
*Substructure path, leave
*END STEP
**
*STEP
*COMPLEX FREQUENCY
*SELECT EIGENMODES, GENERATE
*SUBSTRUCTURE PATH,ENTER=1
*NODE PRINT
 U,
*OUTPUT,FIELD
*NODE OUTPUT
 U,
*Substructure path, leave
*END STEP