/PREP7

SHPP, ON

ET,1,SOLID226,11
*  
KEYOPT,1,1,1001


!ET, 1, SOLID5, 3
EMUNIT, EPZRO, 8.854E-12

con_no = 20                         ! Convergence number
area = lx*ly                                            ! Area of piezo block
f = 1e9                                                  ! Force to apply (100N)
p = f/area                                              ! Pressure applied (Force / area)
*DIM,res,ARRAY,iterations,10
! Initialize variables
d33 = 0
d31 = 0
d32 = 0
d33v = 0
d31v = 0
d32v = 0
per33 = 0
a = 0  

pzcomp = 100 - airpore                          !Calculations to find volume fractions of tri-phasic system
phaseone = (pzcomp/100)*density
phasetwo = pzcomp - phaseone

*DIM,con_res,ARRAY,con_no,4    



elem_total = mesh_size*mesh_size*mesh_size
Lx = 1
Ly = Lx
Lz = Lx
elem_air = 0
elem_poled = 0




