*HEADING           
RESONANT FREQUENCIES OF CLOSED-CIRCUITED PIEZOELECTRIC CUBE    
QUARTER OF CUBE IS MODELED USING SYMMETRY CONDITIONS 
COARSE MODEL WITH C3D20E 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=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,2,2,1,2,20,2,4,200,4
**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,2
21,23
** COMMENT OUT THE NEXT 5 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
4,5
15,15
24,25
31,35,2
41,45
*NSET,NSET=TOPELECT,GENERATE
801,803
811,813,2
821,823
** COMMENT OUT THE NEXT 5 LINES TO OBTAIN THE PARTIALLY COVERING ELECTRODE
804,805
815,815
824,825
831,835,2
841,845
*NSET,NSET=MID,GENERATE
401,405
411,415,2
421,425
431,435,2
441,445
*NSET,NSET=XSIDE,GENERATE
1,5
101,105,2
201,205
301,305,2
401,405
501,505,2
601,605
701,705,2
801,805
*NSET,NSET=YSIDE,GENERATE
1,41,10
101,141,20
201,241,10
301,341,20
401,441,10
501,541,20
601,641,10
701,741,20
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 CLOSED-CIRCUITED CASE, SET BOTH ELECTRODES TO ZERO. 
BOTELECT,9
TOPELECT,9
** GENERATE NSETS USED IN RESTRAINTS FOR BREATHING-TYPE MODES ONLY
*NSET,NSET=DIAG,GENERATE
23,823,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, SIM=NO
5, ,100000.,10000. 
*ELSET,ELSET=ELOUT
1, 
*NSET,NSET=NOUT
1,12,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
*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