Made perm parameter of "perm-holmes-mow" a FEParamDouble.

febiosoftware · Feb 12, 2025 · 384f363 · 384f363
1 parent 54a9e5e
commit 384f363
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 5 deletions.
diff --git a/FEBioMix/FEPermHolmesMow.cpp b/FEBioMix/FEPermHolmesMow.cpp
@@ -84,8 +84,8 @@ mat3ds FEPermHolmesMow::Permeability(FEMaterialPoint& mp)
 	}
 
     // --- strain-dependent isotropic permeability ---
-
-	return mat3dd(m_perm*pow((J-phi0)/(1.0-phi0),m_alpha)*exp(m_M*(J*J-1.0)/2.0));
+	double perm = m_perm(mp);
+	return mat3dd(perm*pow((J-phi0)/(1.0-phi0),m_alpha)*exp(m_M*(J*J-1.0)/2.0));
 }
 
 //-----------------------------------------------------------------------------
@@ -113,8 +113,9 @@ tens4dmm FEPermHolmesMow::Tangent_Permeability_Strain(FEMaterialPoint &mp)
 	}
 
 	mat3dd I(1);	// Identity
-
-	double k0 = m_perm*pow((J-phi0)/(1.0-phi0),m_alpha)*exp(m_M*(J*J-1.0)/2.0);
+	double perm = m_perm(mp);
+
+	double k0 = perm*pow((J-phi0)/(1.0-phi0),m_alpha)*exp(m_M*(J*J-1.0)/2.0);
 	double K0prime = (J*J*m_M+(J*(m_alpha+1)-phi0)/(J-phi0))*k0;
 	mat3ds k0hat = I*K0prime;
 

diff --git a/FEBioMix/FEPermHolmesMow.h b/FEBioMix/FEPermHolmesMow.h
@@ -49,7 +49,7 @@ class FEBIOMIX_API FEPermHolmesMow :	public FEHydraulicPermeability
 	tens4dmm Tangent_Permeability_Strain(FEMaterialPoint& mp) override;
 
 public:
-	double	m_perm;			//!< permeability
+	FEParamDouble	m_perm;			//!< permeability
 	double	m_M;			//!< nonlinear exponential coefficient
 	double	m_alpha;		//!< nonlinear power exponent