*HEADING
: TESTING CPS3 ELEMENTS
** Dynamic Features Tested:
**      *STEADY STATE DYNAMICS,DIRECT
**      *FREQUENCY,EIGENSOLVER=LANCZOS
**      *STEADY STATE DYNAMICS,SUBSPACE PROJECTION
**
*NODE, Nset=ODB_NSET
1,0,0
2,0.5,0
3,1.,0
4,0,0.5
5,0.5,0.5
6,1.,0.5
7,0.,1.
8,0.5,1.
9,1.,1.
*NSET, NSET=NMID
6, 
*ELSET, ELSET=EMID
1,2,5,6
*ELEMENT,TYPE=CPS3,ELSET=E1
1,1,2,4
2,2,5,4
3,2,6,5
4,2,3,6
5,4,8,7
6,4,5,8
7,5,6,8
8,6,9,8
*SOLID SECTION, ELSET=E1, MATERIAL=MAT
*MATERIAL, NAME=MAT
*ELASTIC,TYPE=ISO
2.E7,0.
*DENSITY
8000., 
*DAMPING,ALPHA=5.36,BETA=7.46E-5
*BOUNDARY
1,1,2
4,1,2
7,1,2
**
** Sets added by addodb
**
*Elset,Elset=ODB_ESET
E1,
**
**
*Nset, Nset=ODB_NSET
ODB_NSET, NMID
**
**
*STEP,PERT
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE = LINEAR
10,15.,100,1
*DLOAD
4,P2,30000.
8,P1,30000.
*EL PRINT
S, 
E, 
*NODE PRINT, NSET=NMID
U, 
*NODE FILE, NSET=NMID
U, 
*EL FILE, ELSET=EMID, POSITION=AVERAGED AT NODES, FREQ=42
S, 
*ENERGY PRINT
*ENERGY FILE
*OUTPUT,FIELD,VARIABLE=PRESELECT,OP=ADD,FREQ=42
*OUTPUT,HISTORY,OP=ADD
*NODE OUTPUT,NSET=NMID
U, 
*OUTPUT,HISTORY,OP=ADD
*ELEMENT OUTPUT,ELSET=EMID
S, 
*OUTPUT,HISTORY, OP=ADD
*ENERGY OUTPUT
ALLIE, 
*OUTPUT,FIELD, OP=ADD,FREQ=42
*ENERGY OUTPUT
ALLIE, 
*END STEP
**
*STEP
*FREQUENCY, eigensolver=lanczos
15, , 
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*EL FILE,FREQ=0
*END STEP
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION=EIGENFREQUENCY,FREQUENCY SCALE = LINEAR,
 INTERVAL=RANGE
10,15.,100,1
*DLOAD
4,P2,30000.
8,P1,30000.
*EL PRINT
S, 
E, 
*NODE PRINT, NSET=NMID
U, 
*NODE FILE, NSET=NMID
U, 
*EL FILE, ELSET=EMID, POSITION=AVERAGED AT NODES, FREQ=42
S, 
*END STEP