This data is a part of ongoing research on evolution on hydramechanical response in high porous sandstone during cataclastic flow.

1. Localization_seepage2_300MPa.obd: output of FE simulation software ABAQUS
It simulates drained triaxial compression on high porosity sandsone (Bentheim) under 300 MPa. A damage breakage constitutive model is used for this simulation. Along with this, a steady state flows simulation is also carried out in order to determine flow variation during localized deformation. 
The unit system is Pressure: MPa, Length: m, Time: S, Force: MN. 

	The state variables are following:

	P=STATEV(1)		; mean stress
	Q=STATEV(2)		; deviatoric stress
	EPSV=STATEV(3)		; volumetric strain
	EPSS=STATEV(4)		; deviatoric strain
	B=STATEV(5)		; grain breakage index
	EPSEV=STATEV(6)		; elastic volumetric strain
	EPSPV=STATEV(7)		; plastic volumetric strain
	EPSES=STATEV(8)		; elastic deviatoric strain
	EPSPS=STATEV(9)		; plastic deviatoric strain
	D=STATEV(10)		; bond damage index
      	DD=STATEV(11)		; incremental breakage index
	DB=STATEV(12)		; incremental damage index
	EPSIJP(1:NTENS)=STATEV(13:18)	; plastic strain tensor
      	ECRAND = STATEV(19)	; randemized initial breakage energy
        POR  = STATEV (20)	; porosity

	Parameters for Bentheim sandstone are following:

      	THETAG =PROPS(1)            !GRADING INDEX FOR GRAINS
      	KG = PROPS(2)               !NON-DIMENSIONAL BULK MODULUS FOR GRAINS
      	GG = PROPS(3)               !SHEAR MODULUS GRAIN
      	THETAC = 1.0                !GRADING INDEX CEMENT, BY DEFAULT 1.0
	KC = PROPS(4)               !NON-DIMENSIONAL BULK MODULUS FOR CEMENT
      	GC = PROPS(5)               !SHEAR MODULUS CEMENT
      	OMEGAG = PROPS(6)*PI/180.0  !COUPLING ANGLE GRAIN IN DEGREE 
	OMEGAC = 0.0                !COUPLING ANGLE CEMENT IN DEGREE, BY DEFAULT 0.0
      	M = PROPS(7)                !CRITICAL STATE PARAMETER
      	FC = PROPS(8)               !CEMENT VOLUME FRACTION
      	FG = 1.0-FC                 !GRAIN VOLUME FRACTION = 1- FC
      	PR = PROPS(9)               !REFERENCE PRESSURE, DEFAULT 1KPA, SO CHANGE ACCORDING TO UNIT
      	MM = 0.5                    !DEGREE OF PRESSURE DEPENDNET NON-LINEARITY, DEFAULT 0.5
      	EPSVR = 0.0                 !REFERENCE VOLUMETRIC PRESSURE, DEFAULT 0.0
      	C = PROPS(10)               !COHESION
      	EC = PROPS(11)              !GRAIN CRITICAL ENERGY
      	ET = PROPS(12)              !CEMENT CRITICAL ENERGY
      	ETA = PROPS(13)             !VISCOSITY PARAMETER FOR RATE DEPENDENT MODEL, S/KPA
      	DT = DTIME                  !STEP TIME INCREMEN
      	CON1 = PROPS(14)            !CONSTANT FOR MODIFIED OMEGAG RELATION, USED IN CALCARENITE, OTHERWISE USE ZERO.
      	EC = ECRAND                 !USE IT ONLY WHEN NEED THE RANDDOM DISTRIBUTION

 	PROPS (1:14) are following: 0.9, 42000.,  7588., 42000.,  7588.,    84.,    1.7,   0.05,  0.001,    50.,   3.85,     3., 0.0004,     1.0

2. *.txt files: Are output of duscrete element simulation on high porosity sandsone (Bentheim) using PFC3D code.
The unit system is Pressure: Pa, Length: m, Time: S, Force: N.
The state variables are following:
	deax	:Axial strain
	wszz	:Axial or vertical stress	
	conf	:Mean stress	
	devol	:volumetric strain
	devi	:deviatoric stress (read only low confining compression)

	PFC3D Properties are following
	prop pb_kn=0.7e15 pb_ks=2.8e15 pb_rad=0.6 range z = 0.0,2e-3 
	prop pb_nstren=4.5e8 pb_sstren=14.0e8 range z = 0.0,2e-3 
  	n_stiff = 0.5e7
  	s_stiff = 2.0e7

3. *.dat files: PFC3D ver. 5 input files that consists of ball location and diameter during 1D compression at different axial strain levels (eps). These files directly generates .stl 3D particle assembly which can be read using any image processing tool like (AVIZO, imageJ etc) for further pore structure analysis.