*HEADING
 TEST OF ELEMENT MATRIX INPUT 
 ELEMENT: USER ELEMENT
*NODE
 1, .0110,.0
 2, .01175, 0.0
 3, .0125, 0.0
13, .0110, .0125
14, .01175,.0125
15, .0125,.0125
16, .0110, .014
17, .01175,.014
18, .0125,.014
19, .0110, .0155
20, .01175,.0155
21, .0125,.0155
22, 0.0,.0125
25,.00825,.0125
26, 0.0,.014
27, .0055,.014
28, 0.0,.0155
31,.00825,.0155
*NGEN
1,13,3
2,14,6
3,15,3
22,25
28,31                                   
*USER ELEMENT,TYPE=U1,FILE=eigenpiezotrans_elmatrixout,OLD ELEMENT=1,STEP=1,INCREMENT=1
*USER ELEMENT,TYPE=U2,FILE=eigenpiezotrans_elmatrixout,OLD ELEMENT=2,STEP=1,INCREMENT=1
*USER ELEMENT,TYPE=U3,FILE=eigenpiezotrans_elmatrixout,OLD ELEMENT=3,STEP=1,INCREMENT=1
*USER ELEMENT,TYPE=U4,FILE=eigenpiezotrans_elmatrixout,OLD ELEMENT=4,STEP=1,INCREMENT=1
*USER ELEMENT,TYPE=U5,FILE=eigenpiezotrans_elmatrixout,OLD ELEMENT=5,STEP=1,INCREMENT=1
*ELEMENT,TYPE=U1
1, 1,3,9,7,2,6,8,4
*ELEMENT,TYPE=U2
2, 7,9,15,13,8,12,14,10
*ELEMENT,TYPE=U3
3, 13,15,21,19,14,18,20,16
*ELEMENT,TYPE=U4
4, 22,24,30,28,23,27,29,26
*ELEMENT,TYPE=U5
5, 24,13,19,30,25,16,31,27
*ELSET,ELSET=PZT4
1,2
*ELSET,ELSET=BRASS
3,4,5
*ELSET,ELSET=EALL,GENERATE
1,5,1
*UEL PROPERTY,ELSET=PZT4
*UEL PROPERTY,ELSET=BRASS
*NSET,NSET=IN,GENERATE
1,13,3
*NSET,NSET=SYM
1,2,3
*NSET,NSET=ASYM
22,26,28
*BOUNDARY
SYM,2
ASYM,1
IN,9
*RESTART,WRITE
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos, SIM=NO
5, ,100000.,10. 
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*MODAL FILE
*ELEMENT MATRIX OUTPUT,ELSET=EALL,STIFFNESS=YES,MASS=YES,
OUTPUT FILE=USER DEFINED
*END STEP