SOLUTION
  solver=gdsw
  transient
    # Transient analysis with set time step size
    # and number of steps
    time_step = 1.000000e-06
    nsteps =  1001
END

FILE
  geometry_file 'single_patch.exo'
END

LOADS
END

GDSW
  default_solver = nopivot
END

BOUNDARY
  sideset 5       // symmetry boundary condition
    x = 0   
  sideset 4       // symmetry boundary condition
    y = 0
  sideset 6       // voltage input
    transV = 1
    function voltage_input
  sideset 7       // grounded voltage
    V = 0
END

# all nodes in sideset 8 are set to solve to an
# identical voltage.
RIGIDSET set1
  voltage
  sideset 8
END


FUNCTION voltage_input // voltage input in scaled units (Vin * 1e-9) 
  type linear
  name "voltage_in"
  {include(create_input_deck/voltage_input.inp)}#
END

ECHO
END

OUTPUTS
 disp
 # This is the node-by-node solved for voltage
 # output
 voltage
END

BLOCK 1
  material Aluminum
  # Non-default selective deviatoric hex element formulation
  hex8u
END

BLOCK 2
  material Piezoelectric
  hex8u
END

BLOCK 3
  material Piezoelectric 
  hex8u
END

MATERIAL ALUMINUM
  # Elastic structural material
  density = 2700
  E = {70 * 10^9}
  nu = 0.33
END

// {C11 = 1.38999e+11}
// {C12 = .778366e+11}
// {C13 = .742836e+11}
// {C33 = 1.15412e+11}
// {C44 = 2.5641e+10}
// {C66 = 3.0581e+10}

// {scale = 1e9}

// {ep = 8.85418782e-12 * scale * scale}
// {D11 = ep * 762.5}
// {D33 = ep * 663.2}
 
// {E11 = -5.20279 * scale}
// {E33 = 15.0804 * scale}
// {E15 = 12.7179 * scale}

MATERIAL PIEZOELECTRIC
  ORTHOTROPIC_PIEZOELECTRIC
    # Piezo electric material properties.  Defines
    # the stiffness, conductivity, and coupling between
    # deformation and voltage generation for the material
    Cij = {C11} {C12} {C13}  
                {C11} {C13}
                      {C33}
                      {C44}  
                      {C44}
                      {C66}

    permittivity_ij {D11} 0    0
                    0   {D11}  0  
                    0     0  {D33}
 
     e_ij =  0    0    {E11} 
             0    0    {E11} 
             0    0    {E33} 
             0   {E15}  0
           {E15}   0    0 
            0      0    0	
    density = {7500}
END