*HEADING
.INP: TESTING CPE4IH ELEMENTS
** Dynamic Features Tested:
**      *STEADY STATE DYNAMICS,DIRECT
**      *FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
**      *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,3
*ELEMENT,TYPE=CPE4IH,ELSET=E1
1,1,2,5,4
2,2,3,6,5
3,4,5,8,7
4,5,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
2,P2,30000.
4,P2,30000.
*EL PRINT, ELSET=EMID
S, 
*NODE PRINT, NSET=NMID
U, 
*NODE FILE, NSET=NMID
U, 
*EL FILE, ELSET=EMID, POSITION=AVERAGED AT NODES, FREQ=44
S, 
*ENERGY PRINT
*ENERGY FILE
*OUTPUT,FIELD,OP=ADD,FREQ=44
*NODE OUTPUT,NSET=NMID
U,V,A,RF,CF
*OUTPUT,HISTORY,OP=ADD
*NODE OUTPUT,NSET=NMID
U, 
*OUTPUT,HISTORY,OP=ADD
*ELEMENT OUTPUT,ELSET=EMID
S11,E11
*OUTPUT,HISTORY, OP=ADD
*ENERGY OUTPUT
ALLIE, 
*OUTPUT,FIELD, OP=ADD,FREQ=44
*ENERGY OUTPUT
ALLIE, 
*END STEP
****
*STEP
*FREQUENCY, eigensolver=lanczos, SIM=NO
12, , 
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*OUTPUT,HISTORY,FREQ=0
*EL FILE, ELSET=EMID,FREQ=0
*END STEP
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION=EIGENFREQUENCY,FREQUENCY SCALE = LINEAR,
 INTERVAL=RANGE
10,15.,100,1
*DLOAD
2,P2,30000.
4,P2,30000.
*EL PRINT
S, 
E, 
*NODE PRINT, NSET=NMID
U, 
*NODE FILE, NSET=NMID
U, 
*EL FILE, ELSET=EMID, POSITION=AVERAGED AT NODES, FREQ=44
S, 
*OUTPUT,FIELD,VAR=PRESELECT,FREQ=50
*END STEP
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION=EIGENFREQUENCY,FREQUENCY SCALE = LINEAR,
 INTERVAL=RANGE
10,15.,100,1
*DLOAD
2,TRVEC2,30000.,-1,0,0
4,TRVEC2,30000.,-1,0,0
*END STEP