neo_hookean_J_1_squared_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_NEO_HOOKEAN_J_1_SQUARED_MODEL_H_INCLUDED)
#define  KRATOS_NEO_HOOKEAN_J_1_SQUARED_MODEL_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_models/elasticity_models/neo_hookean_model.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
  class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) NeoHookeanJ_1SquaredModel : public NeoHookeanModel
  {
  public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(NeoHookeanJ_1SquaredModel);
 
 
    NeoHookeanJ_1SquaredModel() : NeoHookeanModel() {}
 
    NeoHookeanJ_1SquaredModel(NeoHookeanJ_1SquaredModel const& rOther) : NeoHookeanModel(rOther) {}
 
    NeoHookeanJ_1SquaredModel& operator=(NeoHookeanJ_1SquaredModel const& rOther)
    {
      NeoHookeanModel::operator=(rOther);
      return *this;
    }
 
    ConstitutiveModel::Pointer Clone() const override
    {
      return Kratos::make_shared<NeoHookeanJ_1SquaredModel>(*this);
    }
 
    ~NeoHookeanJ_1SquaredModel() override {}
 
 
 
 
 
 
 
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "NeoHookeanJ_1SquaredModel";
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "NeoHookeanJ_1SquaredModel";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
      rOStream << "NeoHookeanJ_1SquaredModel Data";
    }
 
 
 
 
  protected:
 
 
 
 
 
 
    // specialized methods:
 
    void CalculateVolumetricFactor(HyperElasticDataType& rVariables, double& rFactor) override
    {
      KRATOS_TRY
 
      rFactor = 0.5 * rVariables.Strain.Invariants.J * (rVariables.Strain.Invariants.J - 1.0);
 
      KRATOS_CATCH(" ")
    }
 
    void CalculateConstitutiveMatrixFactor(HyperElasticDataType& rVariables, double& rFactor) override
    {
      KRATOS_TRY
 
      rFactor = rVariables.Strain.Invariants.I3;
 
      KRATOS_CATCH(" ")
    }
 
    //************// W
 
    //option: g(J) = (lambda/2)*(J-1)² - (mu)*lnJ (neo_hookean_J_1_model)
 
    void CalculateAndAddVolumetricStrainEnergy(HyperElasticDataType& rVariables, double& rVolumetricDensityFunction) override
    {
      KRATOS_TRY
 
      const MaterialDataType& rMaterial = rVariables.GetMaterialParameters();
 
      //g(J) = (lambda/2)*(J-1)² - (mu)*lnJ
      rVolumetricDensityFunction = rMaterial.GetLameLambda() * 0.5 * (rVariables.Strain.Invariants.J - 1.0) * (rVariables.Strain.Invariants.J - 1.0);
      rVolumetricDensityFunction -= rMaterial.GetLameMu() * std::log( rVariables.Strain.Invariants.J );
 
      KRATOS_CATCH(" ")
    }
 
 
    //************// dW
 
 
    double& GetFunction1stI3Derivative(HyperElasticDataType& rVariables, double& rDerivative) override //dW/dI3
    {
      KRATOS_TRY
 
      const MaterialDataType&  rMaterial = rVariables.GetMaterialParameters();
 
      //derivative of "g(J) = (lambda/2)*(J-1)² - (mu)*lnJ"
      //dg(J)/dI3 = (lambda/4)*(J-1)/J - mu*(1/J²)/2
      rDerivative  = 0.25 * rMaterial.GetLameLambda() * ( rVariables.Strain.Invariants.I3 - rVariables.Strain.Invariants.J);
      rDerivative -= 0.5 * rMaterial.GetLameMu();
      rDerivative /= rVariables.Strain.Invariants.I3;
 
      return rDerivative;
 
      KRATOS_CATCH(" ")
    }
 
 
    double& GetFunction2ndI3Derivative(HyperElasticDataType& rVariables, double& rDerivative) override //ddW/dI3dI3
    {
      KRATOS_TRY
 
      const MaterialDataType&   rMaterial = rVariables.GetMaterialParameters();
      //ddg(J)/dI3dI3 = (lambda/4)*(1/J³) + (mu/2)*(1/J⁴)
      rDerivative  = 0.25 * rMaterial.GetLameLambda() * (rVariables.Strain.Invariants.J);
      rDerivative += 0.5 * rMaterial.GetLameMu();
      rDerivative /= (rVariables.Strain.Invariants.I3 * rVariables.Strain.Invariants.I3);
 
      return rDerivative;
 
      KRATOS_CATCH(" ")
    }
 
 
 
 
 
 
 
 
 
  private:
 
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
 
    void save(Serializer& rSerializer) const override
    {
      KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, NeoHookeanModel )
    }
 
    void load(Serializer& rSerializer) override
    {
      KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, NeoHookeanModel )
    }
 
 
 
 
 
  }; // Class NeoHookeanJ_1SquaredModel
 
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_NEO_HOOKEAN_J_1_SQUARED_MODEL_H_INCLUDED  defined