ogden_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_OGDEN_MODEL_H_INCLUDED )
#define  KRATOS_OGDEN_MODEL_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
 
// External includes
 
// Project includes
#include "custom_models/elasticity_models/hyper_elastic_model.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
  class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) OgdenModel : public HyperElasticModel
  {
  public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION( OgdenModel );
 
 
    OgdenModel();
 
    OgdenModel(OgdenModel const& rOther);
 
    OgdenModel& operator=(OgdenModel const& rOther);
 
    ConstitutiveModel::Pointer Clone() const override;
 
 
    ~OgdenModel() override;
 
 
 
 
 
    void CalculateStrainEnergy(ModelDataType& rValues, double& rDensityFunction) override;
 
    void CalculateConstitutiveTensor(ModelDataType& rValues, Matrix& rConstitutiveMatrix) override;
 
 
    int Check(const Properties& rProperties, const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "OgdenModel";
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "OgdenModel";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
      rOStream << "OgdenModel Data";
    }
 
 
 
 
 
  protected:
 
 
 
 
 
    void CalculateAndAddStressTensor(HyperElasticDataType& rVariables, MatrixType& rStressMatrix) override;
 
    void CalculateAndAddConstitutiveTensor(HyperElasticDataType& rVariables, Matrix& rConstitutiveMatrix) override;
 
    virtual void CalculateAndAddConstitutiveTensorB(HyperElasticDataType& rVariables, Matrix& rConstitutiveMatrix);
 
    double& AddConstitutiveComponent(HyperElasticDataType& rVariables, double &rCabcd,
                     const array_1d<double,3>& rVectorDerivative,
                     const unsigned int& a, const unsigned int& b,
                     const unsigned int& c, const unsigned int& d);
 
    void GetEigenCoincidence(const array_1d<double,3>& rStrainEigenValues,
                 array_1d<unsigned int,3>& Order,
                 unsigned int& option);
 
    double& CalculateIsotropicTensorDerivative(const MatrixType& rStrainMatrix,
                           const MatrixType& rStrainEigenVectors,
                           const array_1d<double,3>& rStrainEigenValues,
                           const MatrixType& rStressDerivatives,
                           const array_1d<double,3>& rStressEigenValues,
                           const array_1d<double,6>& rOptionFactors,
                           const unsigned int& rOption,
                           double &rCabcd,
                           const unsigned int& a, const unsigned int& b,
                           const unsigned int& c, const unsigned int& d);
 
 
    void CalculateDerivativeFactors(array_1d<double,6>& rDerivativeFactors, const MatrixType& rStressDerivatives, const array_1d<double,3>& rStressEigenValues, const array_1d<double,3>& rStrainEigenValues, const array_1d<unsigned int,3>& rOrder);
 
 
    //************// Strain Data
 
 
    void CalculateStrainData(ModelDataType& rValues, HyperElasticDataType& rVariables) override;
 
    virtual void CalculateMainStresses(HyperElasticDataType& rVariables, array_1d<double,3>& rMainStresses);
 
    virtual void CalculateMainStressDerivatives(HyperElasticDataType& rVariables, MatrixType& rStressDerivatives);
 
    virtual double& CalculateStressDerivativesI(HyperElasticDataType& rVariables, double& rValue, const unsigned int& i, const unsigned int& j);
    virtual double& CalculateStressDerivativesII(HyperElasticDataType& rVariables, double& rValue, const unsigned int& i, const unsigned int& j);
 
    //************//W
 
    void CalculateAndAddVolumetricStrainEnergy(HyperElasticDataType& rVariables, double& rVolumetricDensityFunction) override;
 
 
 
 
 
 
 
 
  private:
 
 
 
 
 
 
 
 
    using HyperElasticModel::AddConstitutiveComponent;
 
 
 
    friend class Serializer;
 
 
    void save(Serializer& rSerializer) const override
    {
      KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, HyperElasticModel )
    }
 
    void load(Serializer& rSerializer) override
    {
      KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, HyperElasticModel )
    }
 
 
 
 
 
  }; // Class OgdenModel
 
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_OGDEN_MODEL_H_INCLUDED  defined