diff --git a/MaterialLib/MPL/Properties/PermeabilityMohrCoulombFailureIndexModel.cpp b/MaterialLib/MPL/Properties/PermeabilityMohrCoulombFailureIndexModel.cpp index 662a1ba677de18119db7c56a2c0f6db443deaa9e..fd0a1066ecfcf827337dd0025248d8d0e78e7289 100644 --- a/MaterialLib/MPL/Properties/PermeabilityMohrCoulombFailureIndexModel.cpp +++ b/MaterialLib/MPL/Properties/PermeabilityMohrCoulombFailureIndexModel.cpp @@ -150,10 +150,15 @@ PermeabilityMohrCoulombFailureIndexModel<DisplacementDim>::value( template <int DisplacementDim> PropertyDataType PermeabilityMohrCoulombFailureIndexModel<DisplacementDim>::dValue( - VariableArray const& /*variable_array*/, Variable const /*variable*/, + VariableArray const& /*variable_array*/, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { + if (variable == Variable::mechanical_strain) + { + return 0.; + } + OGS_FATAL( "The derivative of the intrinsic permeability k(sigma, ...) with " "respect to stress tensor (sigma) is not implemented because that " diff --git a/MaterialLib/MPL/Properties/StrainDependentPermeability.cpp b/MaterialLib/MPL/Properties/StrainDependentPermeability.cpp index b141bbc1e47446606e83ff2246108efa4737cd48..5769a7cf60c87d5ee1ec72a37d4b637c30bf3919 100644 --- a/MaterialLib/MPL/Properties/StrainDependentPermeability.cpp +++ b/MaterialLib/MPL/Properties/StrainDependentPermeability.cpp @@ -97,10 +97,15 @@ PropertyDataType StrainDependentPermeability<DisplacementDim>::value( template <int DisplacementDim> PropertyDataType StrainDependentPermeability<DisplacementDim>::dValue( - VariableArray const& /*variable_array*/, Variable const /*variable*/, + VariableArray const& /*variable_array*/, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { + if (variable == Variable::mechanical_strain) + { + return 0.; + } + OGS_FATAL( "The derivative of the intrinsic permeability of " "StrainDependentPermeability"