*HEADING
ANTIRESONANT FREQUENCIES OF OPEN-CIRCUITED PIEZOELECTRIC CUBE    
QUARTER OF CUBE IS MODELED USING SYMMETRY CONDITIONS 
REFINED MODEL WITH C3D8E 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=C3D8E
1,1,101,102,2,11,111,112,12
*ELGEN,ELSET=ALL
1,8,1,1,8,10,10,16,100,100
**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
21,25
31,35
41,45
** COMMENT OUT THE NEXT 9 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
6,9
16,19
26,29
36,39
46,49
51,59
61,69
71,79
81,89
*NSET,NSET=TOPELECT,GENERATE
1601,1605
1611,1615
1621,1625
1631,1635
1641,1645
** COMMENT OUT THE NEXT 9 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
1606,1609                
1616,1619
1626,1629
1636,1639
1646,1649
1651,1659
1661,1669
1671,1679
1681,1689
*NSET,NSET=MID,GENERATE
801,809
811,819
821,829
831,839
841,849
851,859
861,869
871,879
881,889
*NSET,NSET=XSIDE,GENERATE
1,9
101,109
201,209
301,309
401,409
501,509
601,609
701,709
801,809
901,909
1001,1009
1101,1109
1201,1209
1301,1309
1401,1409
1501,1509
1601,1609
*NSET,NSET=YSIDE,GENERATE
1,81,10
101,181,10
201,281,10
301,381,10
401,481,10
501,581,10
601,681,10
701,781,10
801,881,10
901,9081,10
1001,1081,10
1101,1181,10
1201,1281,10
1301,1381,10
1401,1481,10
1501,1581,10
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 OPEN-CIRCUITED CASE, SET BOTTOM ELECTRODE ONLY TO ZERO. 
BOTELECT,9
** GENERATE NSETS USED IN RESTRAINTS FOR BREATHING-TYPE MODES ONLY
*NSET,NSET=DIAG,GENERATE
12,1612,100
23,1623,100
34,1634,100
45,1645,100
56,1656,100
67,1667,100
78,1678,100
89,1689,100
*NSET,NSET=XCORNERS,GENERATE
5,1605,100
*NSET,NSET=YCORNERS,GENERATE
41,1641,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,801,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
*END STEP