constitutive_model.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosConstitutiveModelsApplication $
//   Created by:          $Author:                JMCarbonell $
//   Last modified by:    $Co-Author:                         $
//   Date:                $Date:                   April 2017 $
//   Revision:            $Revision:                      0.0 $
//
//
 
#if !defined(KRATOS_CONSTITUTIVE_MODEL_H_INCLUDED )
#define  KRATOS_CONSTITUTIVE_MODEL_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
 
// External includes
 
// Project includes
#include "custom_utilities/constitutive_model_utilities.hpp"
 
#include "custom_models/constitutive_model_data.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
  class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) ConstitutiveModel
  {
  protected:
 
    using VoigtIndexType = const unsigned int(*)[2];
 
  public:
 
    //state flags
    KRATOS_DEFINE_LOCAL_FLAG( ADD_HISTORY_VECTOR );
    KRATOS_DEFINE_LOCAL_FLAG( HISTORY_STRAIN_MEASURE );
    KRATOS_DEFINE_LOCAL_FLAG( HISTORY_STRESS_MEASURE );
 
    typedef ConstitutiveModelData::SizeType                    SizeType;
    typedef ConstitutiveModelData::VectorType                VectorType;
    typedef ConstitutiveModelData::MatrixType                MatrixType;
    typedef ConstitutiveModelData::ModelData              ModelDataType;
    typedef ConstitutiveModelData::MaterialData        MaterialDataType;
 
    typedef ConstitutiveModelData::StrainMeasureType  StrainMeasureType;
    typedef ConstitutiveModelData::StressMeasureType  StressMeasureType;
 
    KRATOS_CLASS_POINTER_DEFINITION( ConstitutiveModel );
 
 
    ConstitutiveModel();
 
    ConstitutiveModel(ConstitutiveModel const& rOther);
 
    virtual ConstitutiveModel::Pointer Clone() const;
 
    ConstitutiveModel& operator=(ConstitutiveModel const& rOther);
 
    virtual ~ConstitutiveModel();
 
 
 
 
 
    virtual void InitializeMaterial(const Properties& rProperties);
 
 
    virtual void InitializeModel(ModelDataType& rValues);
 
    virtual void FinalizeModel(ModelDataType& rValues);
 
 
    virtual void CalculateStrainEnergy(ModelDataType& rValues, double& rDensityFunction);
 
 
    virtual void CalculateStressTensor(ModelDataType& rValues, MatrixType& rStressMatrix);
 
    virtual void CalculateIsochoricStressTensor(ModelDataType& rValues, MatrixType& rStressMatrix);
 
    virtual void CalculateVolumetricStressTensor(ModelDataType& rValues, MatrixType& rStressMatrix);
 
 
    virtual void CalculateConstitutiveTensor(ModelDataType& rValues, Matrix& rConstitutive);
 
    virtual void CalculateIsochoricConstitutiveTensor(ModelDataType& rValues, Matrix& rConstitutive);
 
    virtual void CalculateVolumetricConstitutiveTensor(ModelDataType& rValues, Matrix& rConstitutive);
 
 
    virtual void CalculateStressAndConstitutiveTensors(ModelDataType& rValues, MatrixType& rStressMatrix, Matrix& rConstitutive);
 
    virtual void CalculateIsochoricStressAndConstitutiveTensors(ModelDataType& rValues, MatrixType& rStressMatrix, Matrix& rConstitutive);
 
    virtual void CalculateVolumetricStressAndConstitutiveTensors(ModelDataType& rValues, MatrixType& rStressMatrix, Matrix& rConstitutive);
 
 
    virtual int Check(const Properties& rProperties, const ProcessInfo& rCurrentProcessInfo);
 
 
    virtual bool Has(const Variable<double>& rThisVariable);
 
    virtual void SetValue(const Variable<double>& rVariable, const double& rValue,
              const ProcessInfo& rCurrentProcessInfo);
 
    virtual void SetValue(const Variable<Vector>& rThisVariable, const Vector& rValue,
              const ProcessInfo& rCurrentProcessInfo);
 
    virtual void SetValue(const Variable<Matrix>& rThisVariable, const Matrix& rValue,
              const ProcessInfo& rCurrentProcessInfo);
 
    virtual double& GetValue(const Variable<double>& rThisVariable, double& rValue);
 
    virtual void GetDomainVariablesList(std::vector<Variable<double> >& rScalarVariables,
                    std::vector<Variable<array_1d<double,3> > >& rComponentVariables);
 
 
 
 
    virtual std::string Info() const
    {
        std::stringstream buffer;
        buffer << "ConstitutiveModel";
        return buffer.str();
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const
    {
        rOStream << "ConstitutiveModel";
    }
 
    virtual void PrintData(std::ostream& rOStream) const
    {
      rOStream << "ConstitutiveModel Data";
    }
 
 
 
 
 
  protected:
 
 
 
    Flags mOptions;
 
    //initial or historical strains/stresses
    VectorType mHistoryVector;
 
 
 
 
 
 
 
 
 
 
 
 
  private:
 
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
 
    virtual void save(Serializer& rSerializer) const
    {
      rSerializer.save("mOptions",mOptions);
      rSerializer.save("mHistoryVector",mHistoryVector);
    }
 
    virtual void load(Serializer& rSerializer)
    {
      rSerializer.load("mOptions",mOptions);
      rSerializer.load("mHistoryVector",mHistoryVector);
    }
 
 
 
 
 
  }; // Class ConstitutiveModel
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_CONSTITUTIVE_MODEL_H_INCLUDED  defined