*HEADING
BUCKLING OF A CLAMPED PZT-5H BEAM, Z - Length Polarized (0.005 x 0.006 x 0.400)
**
*NODE, NSET=BEAM 
  1, 0.000, 0.000, 0.000
  2, 0.006, 0.000, 0.000
  3, 0.006, 0.005, 0.000
  4, 0.000, 0.005, 0.000
  5, 0.003, 0.000, 0.000
  6, 0.006, 0.0025, 0.000
  7, 0.003, 0.005, 0.000
  8, 0.000, 0.0025, 0.000
4001, 0.000, 0.000, 0.400
4002, 0.006, 0.000, 0.400
4003, 0.006, 0.005, 0.400
4004, 0.000, 0.005, 0.400
4005, 0.003, 0.000, 0.400
4006, 0.006, 0.0025, 0.400
4007, 0.003, 0.005, 0.400
4008, 0.000, 0.0025, 0.400
**
*NGEN, NSET=BEAM
1, 4001, 100
2, 4002, 100
3, 4003, 100
4, 4004, 100
5, 4005, 100
6, 4006, 100
7, 4007, 100
8, 4008, 100
**
*ELEMENT, TYPE=C3D20E, ELSET=BEAM
1, 1, 2, 3, 4, 201, 202, 203, 204, 5, 6, 7, 8, 205, 206, 207,
208, 101, 102, 103, 104
**
*ELGEN, ELSET=BEAM
1, 20, 200, 1
**
**
**
*SOLID SECTION, ELSET=BEAM, MATERIAL=PZT-5H-Z, ORIENTATION=Z-1
1.
*Orientation, name=Z-1
1.,  0.,   0.,   0.,  1., 0.
1, 0.
**
**
*NSET, NSET=END-1M
1,
*NSET, NSET=END-2M
4001,
*NSET, NSET=END-1S, GENERATE
2,8,1
*NSET, NSET=END-2S, GENERATE
4002,4008,1
*NSET, NSET=CENTER
2003, 2004
**
*NSET, NSET=N-OUT1
2003
*NSET, NSET=N-OUT2
1, 4001
**
**
**
**
*EQUATION
2
END-1S, 9, -1.0, END-1M, 9, 1.0
2
END-2S, 9, -1.0, END-2M, 9, 1.0
3
2007, 1, -1.0, 2003, 1, 0.5, 2004, 1, 0.5
3
2007, 2, -1.0, 2003, 2, 0.5, 2004, 2, 0.5
3
2007, 3, -1.0, 2003, 3, 0.5, 2004, 3, 0.5
**
*MPC
BEAM, END-1S, END-1M
BEAM, END-2S, END-2M
**
** ###################################################################
**
*MATERIAL, NAME=PZT-5H-Z
**
**        For this example the material is Z [3] axis polarized
**    
*ELASTIC, TYPE=ENGINEERING CONSTANTS
6.061e10, 6.061e10, 4.831e10, 0.289, 0.512, 0.512, 2.35E10, 2.30E10
2.30e10,
**  
*DIELECTRIC, TYPE=ORTHO
1.505E-8, 1.505E-8, 1.301E-8
**
*PIEZOELECTRIC, TYPE=E
0.0,   0.0,   0.0,       0.0,  741.0E-12,        0.0,         0.0,   0.0
0.0,   0.0,   0.0, 741.0E-12, -274.0E-12, -274.0E-12,   593.0E-12,   0.0
0.0,   0.0
**
*DENSITY
7500.0
**
**
*BOUNDARY
END-2M, 1, 6
END-1M, 1, 6
**
** -----------------------------------------------------------------
**
*STEP, NAME=STEP-1, NLGEOM=YES, INC=100
Apply a 0.175N Center Load to Induce a Small Geometric Imperfection
*STATIC
.0001,.0001,0.00001,.0001
**
**  Short Circuit Voltage Condition (zero effective Gradient)
**  
*BOUNDARY, OP=MOD
END-1M, 9, 9, 0.0
END-2M, 9, 9, 0.0
**
*CLOAD, OP=NEW
CENTER, 2, -.0875
**
**
*OUTPUT,FIELD,FREQ=0
*EL PRINT, FREQUENCY=0
*NODE PRINT, FREQUENCY=0
**
*END STEP
**
** -----------------------------------------------------------------
**
*STEP, NAME=STEP-2, NLGEOM=YES, INC=1000
IMPOSED VOLTAGE GRADIENT  (Ramp up Voltage Differential to 600000) 
*STATIC
.2,3.5,0.00001,.2
**
*BOUNDARY, OP=MOD
END-2M, 9, 9, -500000.0
**
**
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
EFLX,EPG
**
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,nset=N-OUT2
RF,
*NODE OUTPUT,nset=N-OUT1
U,
**
*END STEP