Merge pull request #14070 from KratosMultiphysics/TFM_Jaume

[CL App] Exposing MINIMUM_FATIGUE_REDUCTION_FACTOR for High cycle fatigue CL

Author: SergioJimenezReyes

applications/ConstitutiveLawsApplication/constitutive_laws_application.cpp

@@ -444,6 +444,7 @@ void KratosConstitutiveLawsApplication::Register()
     KRATOS_REGISTER_VARIABLE(DAMAGE_ACTIVATION)
     KRATOS_REGISTER_VARIABLE(PREVIOUS_CYCLE);
     KRATOS_REGISTER_VARIABLE(CYCLE_PERIOD)
+    KRATOS_REGISTER_VARIABLE(MINIMUM_FATIGUE_REDUCTION_FACTOR)
     KRATOS_REGISTER_VARIABLE(ADVANCE_STRATEGY_APPLIED)
 
 

applications/ConstitutiveLawsApplication/constitutive_laws_application_variables.cpp

@@ -36,6 +36,7 @@ KRATOS_CREATE_VARIABLE(double, TIME_INCREMENT)
 KRATOS_CREATE_VARIABLE(bool, DAMAGE_ACTIVATION)
 KRATOS_CREATE_VARIABLE(double, PREVIOUS_CYCLE)
 KRATOS_CREATE_VARIABLE(double, CYCLE_PERIOD)
+KRATOS_CREATE_VARIABLE(double, MINIMUM_FATIGUE_REDUCTION_FACTOR)
 KRATOS_CREATE_VARIABLE(bool, ADVANCE_STRATEGY_APPLIED)
 
 

applications/ConstitutiveLawsApplication/constitutive_laws_application_variables.h

@@ -46,6 +46,7 @@ namespace Kratos
     KRATOS_DEFINE_APPLICATION_VARIABLE(CONSTITUTIVE_LAWS_APPLICATION, bool, DAMAGE_ACTIVATION)
     KRATOS_DEFINE_APPLICATION_VARIABLE(CONSTITUTIVE_LAWS_APPLICATION, double, PREVIOUS_CYCLE)
     KRATOS_DEFINE_APPLICATION_VARIABLE(CONSTITUTIVE_LAWS_APPLICATION, double, CYCLE_PERIOD)
+    KRATOS_DEFINE_APPLICATION_VARIABLE(CONSTITUTIVE_LAWS_APPLICATION, double, MINIMUM_FATIGUE_REDUCTION_FACTOR)
     KRATOS_DEFINE_APPLICATION_VARIABLE(CONSTITUTIVE_LAWS_APPLICATION, bool, ADVANCE_STRATEGY_APPLIED)
 
 

applications/ConstitutiveLawsApplication/custom_constitutive/auxiliary_files/cl_integrators/high_cycle_fatigue_law_integrator.h

@@ -16,6 +16,8 @@
 // System includes
 
 // Project includes
+#include "includes/define.h"
+#include "includes/checks.h"
 #include "constitutive_laws_application_variables.h"
 #include "custom_utilities/advanced_constitutive_law_utilities.h"
 
@@ -278,7 +280,7 @@ class HighCycleFatigueLawIntegrator
         }
         if (MaxStress > Sth) { //In those cases with no fatigue in course (MaxStress < Sth), tbe fatigue reduction factor does not evolve.
             rFatigueReductionFactor = std::min(rFatigueReductionFactor, std::exp(-B0 * std::pow(std::log10(static_cast<double>(LocalNumberOfCycles)), FatigueReductionFactorSmoothness * (BETAF * BETAF))));
-            const double min_fatigue_reduction_factor = rMaterialParameters[HIGH_CYCLE_FATIGUE_COEFFICIENTS][0];
+            const double min_fatigue_reduction_factor = rMaterialParameters.Has(MINIMUM_FATIGUE_REDUCTION_FACTOR) ? rMaterialParameters[MINIMUM_FATIGUE_REDUCTION_FACTOR] : rMaterialParameters[HIGH_CYCLE_FATIGUE_COEFFICIENTS][0];
             rFatigueReductionFactor = (rFatigueReductionFactor < min_fatigue_reduction_factor) ? min_fatigue_reduction_factor : rFatigueReductionFactor;
         }
     }
@@ -305,6 +307,18 @@ class HighCycleFatigueLawIntegrator
         return ultimate_stress;
     }
 
+    /**
+     * @brief This method defines in the CL integrator
+     * @return 0 if OK, 1 otherwise
+     */
+    static int Check(const Properties& rMaterialProperties)
+    {
+        const double min_fatigue_reduction_factor = rMaterialProperties.Has(MINIMUM_FATIGUE_REDUCTION_FACTOR) ? rMaterialProperties[MINIMUM_FATIGUE_REDUCTION_FACTOR] : rMaterialProperties[HIGH_CYCLE_FATIGUE_COEFFICIENTS][0];
+        KRATOS_ERROR_IF(min_fatigue_reduction_factor <= 0.0 || min_fatigue_reduction_factor >= 1.0) << "Minimum fatigue reduction factor must be in (0,1). Provided: " << min_fatigue_reduction_factor << std::endl;
+
+        return 0;
+    }
+
     ///@}
     ///@name Access
     ///@{

applications/ConstitutiveLawsApplication/custom_constitutive/small_strains/fatigue/generic_small_strain_high_cycle_fatigue_law.cpp

@@ -653,6 +653,22 @@ Matrix& GenericSmallStrainHighCycleFatigueLaw<TConstLawIntegratorType>::Calculat
 /***********************************************************************************/
 /***********************************************************************************/
 
+template <class TConstLawIntegratorType>
+int GenericSmallStrainHighCycleFatigueLaw<TConstLawIntegratorType>::Check(
+    const Properties& rMaterialProperties,
+    const GeometryType& rElementGeometry,
+    const ProcessInfo& rCurrentProcessInfo
+    ) const
+{
+    const int check_base = BaseType::Check(rMaterialProperties, rElementGeometry, rCurrentProcessInfo);
+    const int check_integrator = HighCycleFatigueLawIntegrator<6>::Check(rMaterialProperties);
+    if ((check_base + check_integrator) > 0) return 1;
+    return 0;
+}
+
+/***********************************************************************************/
+/***********************************************************************************/
+
 template class GenericSmallStrainHighCycleFatigueLaw <GenericConstitutiveLawIntegratorDamage<VonMisesYieldSurface<VonMisesPlasticPotential<6>>>>;
 template class GenericSmallStrainHighCycleFatigueLaw <GenericConstitutiveLawIntegratorDamage<ModifiedMohrCoulombYieldSurface<ModifiedMohrCoulombPlasticPotential<6>>>>;
 template class GenericSmallStrainHighCycleFatigueLaw <GenericConstitutiveLawIntegratorDamage<TrescaYieldSurface<TrescaPlasticPotential<6>>>>;

applications/ConstitutiveLawsApplication/custom_constitutive/small_strains/fatigue/generic_small_strain_high_cycle_fatigue_law.h

@@ -317,6 +317,21 @@ class KRATOS_API(CONSTITUTIVE_LAWS_APPLICATION) GenericSmallStrainHighCycleFatig
         Matrix& rValue
         ) override;
 
+    /**
+     * @brief This function provides the place to perform checks on the completeness of the input.
+     * @details It is designed to be called only once (or anyway, not often) typically at the beginning
+     * of the calculations, so to verify that nothing is missing from the input or that no common error is found.
+     * @param rMaterialProperties The properties of the material
+     * @param rElementGeometry The geometry of the element
+     * @param rCurrentProcessInfo The current process info instance
+     * @return 0 if OK, 1 otherwise
+     */
+    int Check(
+        const Properties& rMaterialProperties,
+        const GeometryType& rElementGeometry,
+        const ProcessInfo& rCurrentProcessInfo
+        ) const override;
+
     /**
      * @brief If the CL requires to initialize the material response, called by the element in InitializeSolutionStep.
      */
@@ -507,7 +522,7 @@ class KRATOS_API(CONSTITUTIVE_LAWS_APPLICATION) GenericSmallStrainHighCycleFatig
         rSerializer.load("PreviousCycleTime", mPreviousCycleTime);
         rSerializer.load("Period", mPeriod);
         rSerializer.load("ReferenceDamage", mReferenceDamage);
-        rSerializer.save("PreviousCycleDamage", mPreviousCycleDamage);
+        rSerializer.load("PreviousCycleDamage", mPreviousCycleDamage);
         rSerializer.load("FirstCycleOfANewLoad", mFirstCycleOfANewLoad);
         rSerializer.load("CFactor", mCFactor);
     }

applications/ConstitutiveLawsApplication/custom_python/constitutive_laws_python_application.cpp

@@ -62,6 +62,7 @@ PYBIND11_MODULE(KratosConstitutiveLawsApplication,m)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, DAMAGE_ACTIVATION)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, PREVIOUS_CYCLE);
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, CYCLE_PERIOD)
+    KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, MINIMUM_FATIGUE_REDUCTION_FACTOR)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, ADVANCE_STRATEGY_APPLIED)