Another one for double porosity and advanced solid model
Pre-release
This release is actually a zip file that contains the numerical scheme for a coupled double porosity hydromechanical model. The anisotropic modified Cam-Clay model is adopted to simulate the behavior of solid (shale-like material). It differs from the original MCC model even in the isotropic case. Both the solid anisotropy and permeability anisotropy are affected by the orientation theta (counterclockwise is positive), whose default value is 45 [deg]. If you want to change theta, you only need to change it at the beginning.
There are two main simulation files. The second file approximately simulates the example from Zhao and Borja (2021), since we DO NOT handle a finite Ks (solid grain bulk modulus) case. In that paper, the A matrix (constant strain storage coefficient matrix) only considers the fluid compressibility effect. The first file considers the internal structure evolution through the matrix A, but the fluid compressibility is ignored, see Zhang et al. (2022).
If you want to try anisotropic elasticity, you only need to change the following sentence in the AMCC_UMAT.m file:
if YIELD_VALUE < 1E-12 % Elastic
to
if YIELD_VALUE < 1E12 % Elastic
Post_processing.m has not been modified properly. The original purpose was to save the MATLAB results as a bunch of vtk files and import them into Paraview. paraview_script.py is also for Paraview.
It also contains the data CoolWarmFloat257.csv for cool-warm colormap.
Following figures (p1, p2, p1, p2) and (0.6 d, 0.6 d, 56 d, 56 d). The data unit is MPa. In double porosity theory,
