This workflow outlines the modeling steps used to simulate the electromechanical properties freeze-cast 2-2 structures.
The FE models demonstrate the effects of 1) the thickness (i.e., grain clamping), 2) wall defects (i.e., pores within 
ceramic walls) and 3) pore channel defects (i.e., grains within ceramic channels) on the bulk electromechanical properties.

Geometry and material property set naming convention:
t - thickness
b - percentage volume of pore channel defects
w - percentage volume of wall defects
v - variant of model

For example: materials_t_2_b_0_w_1_v_3 means the material property set of a 2-2 model with thickness of 2, 0 vol.% of
pore channel defects and 10 vol.% of wall defects, model number 3.

------------------------
1. Structure generation (MATLAB)
   - `Cube_set_up.m` creates a network model with randomly distributed pores within a ceramic matrix.
   - The output structure is saved as `assigned.txt`.

2. Polarization initialization
   - `Cube_set_up.m` calls `transform.m` to assign a randomized initial polarization state.
   - This polarization data is saved in `materials.txt`, which corresponds to `assigned.txt`.

3. Simulation setup (ANSYS APDL)
   - `master.txt` orchestrates the APDL simulation workflow.
   - `start_script.txt` initializes a coupled-field piezoelectric simulation.

4. Material registration
   - Material properties are imported into APDL from `materials.txt`.

5. Geometry definition
   - A cubic geometry is defined using `geom.txt`.

6. Structure import
   - The pore structure from `assigned.txt` is read into APDL to define the model's microstructure.

7. Boundary conditions
   - `boundary_conditions.txt` applies appropriate mechanical and electrical boundary conditions for the piezoelectric model.

8. Electromechanical response
   - `determine_res.txt` applies switching conditions and solves for the electromechanical properties of the model.echanical properties of the model.