This workflow outlines the modeling steps used to simulate the electric field distribution in a perfect 2-2 ceramic 
model with wall defects (pores within ceramic walls)
------------------------

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_5.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_5.txt`, which corresponds to `assigned_5.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_5.txt`.

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

6. Structure import
   - The pore structure from `assigned_5.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.