*HEADING
COMPLEX EIGENSOLVER - SUBSTRUCTURE TEST
SUBSTRUCTURE GENERATION FILE pcfreq_sup_gen.inp SHOULD BE RUN FIRST.
*NODE, NSET=NROTOR
352  ,  1.0607    ,  1.0607    ,  0.0000    ,
602  , 2.42331E-15,  1.5000    ,  0.0000    ,
852  , -1.0607    ,  1.0607    ,  0.0000    ,
353  ,  1.4142    ,  1.4142    ,  0.0000    ,
603  , 3.23109E-15,  2.0000    ,  0.0000    ,
853  , -1.4142    ,  1.4142    ,  0.0000    ,
452  ,  1.0607    ,  1.0607    ,  1.0000    ,
702  , 2.42331E-15,  1.5000    ,  1.0000    ,
952  , -1.0607    ,  1.0607    ,  1.0000    ,
453  ,  1.4142    ,  1.4142    ,  1.0000    ,
703  , 3.23109E-15,  2.0000    ,  1.0000    ,
953  , -1.4142    ,  1.4142    ,  1.0000    ,
*ELEMENT,TYPE=Z1, ELSET=EROTOR, FILE=pcfreq_sup_gen
1, 352,353,602,603,852,853,452,453,702,703,952,953
*SUBSTRUCTURE PROPERTY,ELSET=EROTOR
*NSET ,NSET= NSET_SIDE1
352,353,602,603,852,853
*NSET ,NSET= NSET_SIDE2
452,453,702,703,952,953
*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
.0,.0,.0,  .0,.0,1.
1,  90.0,
*NODE, NSET=NPAD1
5, -1.8,  0.2, 0.
6,  1.8,  0.2, 0.
7,  1.8,  2.1, 0.
8, -1.8,  2.1, 0.
*ELEMENT, TYPE=M3D4, ELSET=ELPAD1
11, 5,6,7,8
*NODE, NSET=NPAD2
15, -1.8,  0.2, 1.
16,  1.8,  0.2, 1.
17,  1.8,  2.1, 1.
18, -1.8,  2.1, 1.
*ELEMENT, TYPE=M3D4, ELSET=ELPAD2
12, 15,16,17,18
*MATERIAL, NAME=A2
*ELASTIC
1000, .3
*DENSITY
1.,
*MEMBRANESECTION, ELSET=ELPAD1, MATERIAL=A2
1.,
*MEMBRANESECTION, ELSET=ELPAD2, MATERIAL=A2
1.,
*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,Convert Sdi=YES
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,Convert Sdi=YES
*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,Convert Sdi=YES
REAL EIGENSOLVER
*FREQUENCY, SIM
12,
*SUBSTRUCTURE PATH,ENTER=1
*NODE PRINT
U,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*Substructure path, leave
*END STEP
**
*STEP,Convert Sdi=YES
*COMPLEX FREQUENCY
*SELECT EIGENMODES, GENERATE
1,10
*SUBSTRUCTURE PATH,ENTER=1
*NODE PRINT
U,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*Substructure path, leave
*END STEP