*HEADING
ANTIRESONANT FREQUENCIES OF OPEN-CIRCUITED PIEZOELECTRIC CUBE    
QUARTER OF CUBE IS MODELED USING SYMMETRY CONDITIONS 
COARSE MODEL WITH C3D8E ELEMENTS
*RESTART,WRITE
*NODE
 1, 0.0, 0.0,0.
 5, .01, 0.0,0.0
41, 0.0, .01,0.
45,.01, .01,0.0
801,0.,0.,.02
805,.01,0.,.02
841,0.,.01,.02
845,.01,.01,.02
*NGEN,NSET=A
1,5
*NGEN,NSET=B
41,45
*NGEN,NSET=C
801,805
*NGEN,NSET=D
841,845
*NFILL,NSET=SIDEA
A,B,4,10
*NFILL,NSET=SIDEB
C,D,4,10
*NFILL
SIDEA,SIDEB,8,100
*ELEMENT,TYPE=C3D8E
1,1,101,102,2,11,111,112,12
*ELGEN,ELSET=ALL
1,4,1,1,4,10,10,8,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,3
11,13
21,23
** COMMENT OUT THE NEXT 5 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
4,5
14,15
24,25
31,35
41,45
*NSET,NSET=TOPELECT,GENERATE
801,803
811,813
821,823
** COMMENT OUT THE NEXT 5 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
804,805
814,815
824,825
831,835
841,845
*NSET,NSET=MID,GENERATE
401,405
411,415
421,425
431,435
441,445
*NSET,NSET=XSIDE,GENERATE
1,5
101,105
201,205
301,305
401,405
501,505
601,605
701,705
801,805
*NSET,NSET=YSIDE,GENERATE
1,41,10
101,141,10
201,241,10
301,341,10
401,441,10
501,541,10
601,641,10
701,741,10
801,841,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,812,100
23,823,100
34,834,100
45,845,100
*NSET,NSET=XCORNERS,GENERATE
5,805,100
*NSET,NSET=YCORNERS,GENERATE
41,841,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
5, ,100000.,10000. 
*ELSET,ELSET=ELOUT
1, 
*NSET,NSET=NOUT
1,12,401,623,801,812
*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