*HEADING
MORGAN MATROC PZT-5H TEST BLOCK FOR MATERIAL VERIFICATION, STATIC NLGEOM=NO
**
*RESTART,WRITE
** ####################################################################
**
**    MODEL DATA SECTION
**
*NODE, NSET=BLOCK
1, 0., 0., 0.
2, 1., 0., 0.
3, 1., 1., 0.
4, 0., 1., 0.
5, 0., 0., 1.
6, 1., 0., 1.
7, 1., 1., 1.
8, 0., 1., 1.
**
*ELEMENT, TYPE=C3D8E, ELSET=BLOCK
1, 1, 2, 3, 4, 5, 6, 7, 8
**
**
*SOLID SECTION, ELSET=BLOCK, MATERIAL=PZT-5H, ORIENTATION=global-123
**
*Orientation, Name=global-123, System=Rectangular
1.0, 0.0, 0.0, 0.0, 1.0, 0.0
1, 0.0
**
**
** ###################################################################
**
**  Define Node Sets to be used for BC and output control. 
**
*NSET, NSET=TOP
5,6,7,8
*NSET, NSET=BOT
1,2,3,4
*NSET, NSET=Left
1,4,5,8
*NSET, NSET=Right
2,3,6,7
**
** ###################################################################
**
*MATERIAL, NAME=PZT-5H
**
**        For this example the material is Z [3] axis polarized
**    
*ELASTIC, TYPE=ENGINEERING CONSTANTS
5.061e10, 5.061e10, 3.831e10, 0.279, 0.502, 0.502, 1.35E10, 1.30E10
1.30e10,0.
7.061e10, 7.061e10, 5.831e10, 0.299, 0.522, 0.522, 3.35E10, 3.30E10
3.30e10,200.
**  
*DIELECTRIC, TYPE=ORTHO
0.505E-8, 0.505E-8, 0.301E-8,0.
2.505E-8, 2.505E-8, 2.301E-8,200.
**
*PIEZOELECTRIC, TYPE=E
0.0,   0.0,   0.0,       0.0,  641.0E-12,        0.0,         0.0,   0.0
0.0,   0.0,   0.0, 641.0E-12, -174.0E-12, -174.0E-12,   493.0E-12,   0.0
0.0,   0.0, 0.
0.0,   0.0,   0.0,       0.0,  841.0E-12,        0.0,         0.0,   0.0
0.0,   0.0,   0.0, 841.0E-12, -374.0E-12, -374.0E-12,   693.0E-12,   0.0
0.0,   0.0, 200.
**
*DENSITY
7500.0
**
** ###################################################################
**
**  RESTRAIN 6 RIGID BODY MODES
**
*BOUNDARY
1, 1, 3, 0.0
2, 3, 3, 0.0
4, 3, 3, 0.0
2, 2, 2, 0.0
**
** ###################################################################
**
*INITIAL CONDITIONS,TYPE=TEMPERATURE
BLOCK,100.
**  Note, for these tests the 3 axis is the polarized axis.
**
*STEP, NAME=STEP-1, NLGEOM=YES
IMPOSED VOLTAGE GRADIENT (Z AXIS): Top = 0, Bot = 1 
*STATIC
1., 1.
**
**   Verify Material d33, d31, g33, and g31 values
**
**   d333 = Strain in 3 direction divided by the Potential Gradient.  
**
**   d311 = Strain in 1 & 2 directions divided by the Potential Gradient.  
**   
**   g333 = Strain in 3 direction divided by the Reaction Charge per
**         unit area (Electrical Flux).
** 
**   g311 = Strain in 1 or 2 direction divided by the Reaction Charge
**         per unit area (Electrical Flux)
**
**   Also, Verify the value for the effective relative dielectric constant 
**         in the polarized direction for a system that is free to expand
**         (no external forces either applied or via boundary conditions).
**         This is equal to the Electrical Flux (Charge per Unit Area) 
**         associated with the applied voltage divided by the magnitude of
**         Voltage Gradient divided by the dielectric constant of free space. 
**
**             
**           Manufacturer's Data Sheet   
**					 
**             d333 =  593.0e-12          
**             d311 = -274.0e-12          
**             g333 =  19.7e-3            
**             g311 = -9.11e-3            
**             DC-33(free)= 3400         
**
*BOUNDARY, OP=MOD
BOT, 9, 9, 1.0
TOP, 9, 9, 0.0
**
*OUTPUT, FIELD, FREQUENCY=1
*NODE OUTPUT
U, RF, EPOT, RCHG, CECHG
*ELEMENT OUTPUT
S, E, EPG, EFLX
**
***NODE FILE, FREQUENCY=1
**U,  EPOT
**RF, RCHG, CECHG
**
***EL FILE, FREQUENCY=1
**S,
**E,
**EPG,
**EFLX,
**CHRGS,
**
*END STEP
**
**   #################################################################
**
*STEP, NAME=STEP-2, NLGEOM=YES
IMPOSED VOLTAGE GRADIENT (Z AXIS ): Top = 1, Bot = 0 
**
*STATIC
1.,1.
**
**   Validate Material Data (see step 1)
**
*BOUNDARY, OP=MOD
BOT, 9, 9, 0.0
TOP, 9, 9, 1.0
**
*END STEP
**
**   #################################################################
**
**  Verify that nodal charge loading is working correctly
**
*STEP, NAME=STEP-3, NLGEOM=YES
Voltage Bottom = 0.0, Top Charge = 3.042164e-8 Coulombs
**
*STATIC
1.,1.
**
*BOUNDARY, OP=NEW
BOT, 9, 9,   0.0
1, 1, 3, 0.0
2, 3, 3, 0.0
4, 3, 3, 0.0
2, 2, 2, 0.0
**
**  Apply Nodal Charges
**
*CECHARGE, OP=NEW
TOP, , 7.60541E-9
**
*END STEP
**
**
** ###################################################################
**
**  Verify that distributed charge loading is working correctly
**
*STEP, NAME=STEP-4, NLGEOM=YES
Voltage Bottom = 0.0, Top Charge = 3.042164e-8 Coulombs
**
*STATIC
1.,1.
**
**  Apply charge as a distributed load 
**
*CECHARGE, OP=NEW
*DECHARGE, OP=NEW
1, ES2, 3.042164E-8
**
*END STEP
**
** ###################################################################
**
*STEP, NAME=STEP-5, NLGEOM=YES
Validate Material d15 & g15 Parameter: 1V Gradient in 1 Direction
**
**  d113 = Shear strain about 2 axis 1-3 plane divided by the Potential
**        Gradient along the 1 axis (3-axis is the axis of polarization).
**
**  g113 = Shear strain about 2 axis (E13) divided by the Reaction Charge
**        per Unit Area (Electrical Flux).
**
**   Also, Verify the value for the effective relative dialectric constant 
**         in the cross-polarized direction for a system free to expand
**         (no external forces either applied or via boundary conditions).
**         This is equal to the Electrical Flux (Charge per Unit Area) 
**         associated with the applied voltage divided by the magnitude of
**         Voltage Gradient divided by the dielectric constant of free space. 
**
**             
**           Manufacturer's Data Sheet   
**					 
**             d113 =  741.0e-12          
**             g113 =   26.8e-3           
**             DC-11(free)= 3130         
**
**
*STATIC
1.,1.
**
*BOUNDARY, OP=NEW
LEFT,  9, 9, 1.0
RIGHT, 9, 9, 0.0
1, 1, 3, 0.0
2, 3, 3, 0.0
4, 3, 3, 0.0
2, 2, 2, 0.0
*DECHARGE, OP=NEW
*CECHARGE, OP=NEW
**
*END STEP
**
** ###################################################################
**
*STEP, NAME=STEP-6, NLGEOM=YES
Validate The Open Circuit Compliance: 1 Pa Tension along Polarized (3) axis
**
**  Note:  can also be used to verify g material values.
**
**      Manufacturer's Data Sheet   
**         			    
**          SD-33 =  8.99e-12       
**          SD-31 = -3.05e-12       
**
*STATIC
1.,1.
**
*BOUNDARY, OP=NEW
BOT, 9, 9, 0.0
BOT, 3, 3, 0.0
1,2,2,0.0
2,2,2,0.0
5,2,2,0.0
6,2,2,0.0
1,1,1,0.0
4,1,1,0.0
5,1,1,0.0
8,1,1,0.0
*DECHARGE, OP=NEW
*CECHARGE, OP=NEW
*DLOAD, OP=NEW
1,P2,-1.
**
*END STEP
**
** ###################################################################
**
*STEP, NAME=STEP-7, NLGEOM=YES
Validate The Open Circuit Compliance: 1 Pa Tension along (1) axis
**
**       Manufacturer's Data Sheet        
**         				  
**         SD-11 =  14.05e-12             
**         SD-12 =  -7.27e-12             
**         SD-13 =  -3.05e-12           
**
*STATIC
1.,1.
**
*BOUNDARY, OP=NEW
BOT, 9, 9, 0.0
BOT, 3, 3, 0.0
1,2,2,0.0
2,2,2,0.0
5,2,2,0.0
6,2,2,0.0
1,1,1,0.0
4,1,1,0.0
5,1,1,0.0
8,1,1,0.0
*DECHARGE, OP=NEW
*CECHARGE, OP=NEW
*DLOAD, OP=NEW
1,P4,-1.
**
*END STEP
**
**   #################################################################
**
*STEP, NAME=STEP-8, NLGEOM=YES
Short Circuit Condition -- Hydrostatic Pressure Load of 1 PA 
**
**   Validate Material dh value 
**   dh = 45.0e-12 = Reaction Charge per Unit Area (Electrical Flux)
**                   divided by the value of the Hydrostatic Pressure.
**                   (0.0% Diff with manufacturer's Data Sheet) 
**
*STATIC
1.,1.
**
*BOUNDARY, OP=NEW
TOP, 9, 9,   0.0
BOT, 9, 9,   0.0
BOT, 3, 3,   0.0
1,2,2,0.0
2,2,2,0.0
5,2,2,0.0
6,2,2,0.0
1,1,1,0.0
4,1,1,0.0
5,1,1,0.0
8,1,1,0.0
*DECHARGE, OP=NEW
*CECHARGE, OP=NEW
*DLOAD, OP=NEW
1,P2,1.
1,P4,1.
1,P5,1.
**
*END STEP