*HEADING
RESONANT FREQUENCIES OF CLOSED-CIRCUITED PIEZOELECTRIC CUBE    
QUARTER OF CUBE IS MODELED USING SYMMETRY CONDITIONS 
REFINED MODEL WITH C3D20E ELEMENTS
*RESTART,WRITE
*NODE
 1, 0.0,0.0,0.0
 9, .01,0.0,0.0
81, 0.0,.01,0.0
89, .01,.01,0.0
1601,0.0,0.0,.02
1609,.01,0.0,.02
1681,0.0,.01,.02
1689,.01,.01,.02
*NGEN,NSET=A
1,9
*NGEN,NSET=B
81,89
*NGEN,NSET=C
1601,1609
*NGEN,NSET=D
1681,1689
*NFILL,NSET=SIDEA
A,B,8,10
*NFILL,NSET=SIDEB
C,D,8,10
*NFILL
SIDEA,SIDEB,16,100
*ELEMENT,TYPE=C3D20E
1,1,3,23,21,201,203,223,221,2,13,22,11,202,213,222,
211,101,103,123,121
*ELGEN,ELSET=ALL
1,4,2,1,4,20,4,8,200,16
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,MATERIAL=PZT,ELSET=ALL,ORIENT=RECT
*MATERIAL,NAME=PZT
*ELASTIC,TYPE=ORTHOTROPIC
13.2E10,7.1E10,13.2E10,7.3E10,7.3E10,11.5E10,3.0E10,2.6E10
2.6E10, 
*PIEZOELECTRIC,TYPE=S
0.,0.,0.,0.,10.5,0.,0.,0.
0.,0.,0.,10.5,-4.1,-4.1,14.1,0.,
0.,0.
*DIELECTRIC,TYPE=ORTHO
7.124E-9,7.124E-9,5.841E-9
*DENSITY
7.5E3, 
** GENERATE NODE SETS FOR THE TOP AND BOTTOM ELECTRODES.
** CURRENTLY IT IS SETUP FOR FULLY COVERING ELECTRODES. COMMENT OUT THE 
** NODE SPECIFICATIONS IN THE FOLLOWING WHERE NOTED TO OBTAIN THE 
** PARTIALLY COVERING ELECTRODES.
*NSET,NSET=BOTELECT,GENERATE
1,5
11,15,2
21,25
31,35,2
41,45
** COMMENT OUT THE NEXT 9 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
6,9
17,19,2
26,29
37,39,2
46,49
51,59,2
61,69
71,79,2
81,89
*NSET,NSET=TOPELECT,GENERATE
1601,1605
1611,1615,2
1621,1625
1631,1635,2
1641,1645
** COMMENT OUT THE NEXT 9 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
1606,1609                
1617,1619,2
1626,1629
1637,1639,2
1646,1649
1651,1659,2
1661,1669
1671,1679,2
1681,1689
*NSET,NSET=MID,GENERATE
801,809
811,819,2
821,829
831,839,2
841,849
851,859,2
861,869
871,879,2
881,889
*NSET,NSET=XSIDE,GENERATE
1,9
101,109,2
201,209
301,309,2
401,409
501,509,2
601,609
701,709,2
801,809
901,909,2
1001,1009
1101,1109,2
1201,1209
1301,1309,2
1401,1409
1501,1509,2
1601,1609
*NSET,NSET=YSIDE,GENERATE
1,81,10
101,181,20
201,281,10
301,381,20
401,481,10
501,581,20
601,681,10
701,781,20
801,881,10
901,981,20
1001,1081,10
1101,1181,20
1201,1281,10
1301,1381,20
1401,1481,10
1501,1581,20
1601,1681,10
*BOUNDARY
** SYMMETRY ON PLANES OF SYMMETRY
XSIDE,YSYMM
YSIDE,XSYMM
** FOR BREATHING-TYPE MODES, DISPLACEMENTS HAVE SYMMETRY ABOUT XY-PLANE
MID,3
** FOR CLOSED-CIRCUITED CASE, SET BOTTOM AND TOP ELECTRODES TO ZERO. 
TOPELECT,9
BOTELECT,9
** GENERATE NSETS USED IN RESTRAINTS FOR BREATHING-TYPE MODES ONLY
*NSET,NSET=DIAG,GENERATE
23,1623,100
45,1645,100
67,1667,100
89,1689,100
*NSET,NSET=XCORNERS,GENERATE
9,1609,100
*NSET,NSET=YCORNERS,GENERATE
81,1681,100
** APPLY EQUATIONS TO ENSURE BREATHING-TYPE MODES ONLY ARE EXTRACTED
*EQUATION
2, 
DIAG,1,1.,DIAG,2,-1.
2, 
XCORNERS,1,1.,YCORNERS,2,-1.
2, 
XCORNERS,3,1.,YCORNERS,3,-1.
*STEP,PERTURBATION
EXTRACT BREATHING MODES UP TO 100000 HZ
*FREQUENCY, eigensolver=lanczos, SIM=NO
5, ,100000.,10000.
*ELSET,ELSET=ELOUT
1, 
*NSET,NSET=NOUT
1,12,1023,1601,1612
*EL PRINT,ELSET=ELOUT
S, 
E, 
EFLX,EPG       
*NODE PRINT,NSET=NOUT
U,EPOT
*NODE FILE, LAST MODE=1, NSET=NOUT
U,EPOT,RF,RCHG
*EL FILE, LAST MODE=1, ELSET=ELOUT
S,E,EFLX,EPG
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NOUT
U,EPOT,RF,RCHG
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NOUT
U,EPOT,RF,RCHG
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=ELOUT
S,E,EFLX,EPG
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=ELOUT
S,E,EFLX,EPG
*END STEP