small_strain_axisymmetric_2D_law.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_SMALL_STRAIN_AXISYMMETRIC_2D_LAW_H_INCLUDED)
#define  KRATOS_SMALL_STRAIN_AXISYMMETRIC_2D_LAW_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_laws/small_strain_laws/small_strain_3D_law.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
  class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) SmallStrainAxisymmetric2DLaw : public SmallStrain3DLaw
    {
    public:
 
      KRATOS_CLASS_POINTER_DEFINITION(SmallStrainAxisymmetric2DLaw);
 
 
      SmallStrainAxisymmetric2DLaw() : SmallStrain3DLaw() {}
 
      SmallStrainAxisymmetric2DLaw(ModelTypePointer pModel) : SmallStrain3DLaw(pModel) {};
 
      SmallStrainAxisymmetric2DLaw(const SmallStrainAxisymmetric2DLaw& rOther) : SmallStrain3DLaw(rOther) {}
 
      SmallStrainAxisymmetric2DLaw& operator=(SmallStrainAxisymmetric2DLaw const& rOther)
      {
    SmallStrain3DLaw::operator=(rOther);
    return *this;
      }
 
      ConstitutiveLaw::Pointer Clone() const override
      {
        return Kratos::make_shared<SmallStrainAxisymmetric2DLaw>(*this);
      }
 
      ~SmallStrainAxisymmetric2DLaw() override {}
 
 
 
      SizeType WorkingSpaceDimension() override { return 2; }
 
      SizeType GetStrainSize() const override { return 4; }
 
      void GetLawFeatures(Features& rFeatures) override
      {
    KRATOS_TRY
 
        //Set the type of law
    rFeatures.mOptions.Set( AXISYMMETRIC_LAW );
    rFeatures.mOptions.Set( INFINITESIMAL_STRAINS );
    rFeatures.mOptions.Set( ISOTROPIC );
 
    //Get model features
    GetModelFeatures(rFeatures);
 
    //Set strain measure required by the consitutive law
    rFeatures.mStrainMeasures.push_back(StrainMeasure_Infinitesimal);
    rFeatures.mStrainMeasures.push_back(StrainMeasure_Deformation_Gradient);
 
    //Set the strain size
    rFeatures.mStrainSize = GetStrainSize();
 
    //Set the spacedimension
    rFeatures.mSpaceDimension = WorkingSpaceDimension();
 
    KRATOS_CATCH(" ")
      }
 
 
 
 
 
 
 
 
      std::string Info() const override
      {
    std::stringstream buffer;
        buffer << "SmallStrainAxisymmetric2DLaw" ;
        return buffer.str();
      }
 
      void PrintInfo(std::ostream& rOStream) const override {rOStream << "SmallStrainAxisymmetric2DLaw";}
 
      void PrintData(std::ostream& rOStream) const override {}
 
 
 
 
 
    protected:
 
 
 
 
 
 
      void CalculateConstitutiveMatrix(Matrix& rConstitutiveMatrix,
                    const Properties& rProperties) override
      {
    KRATOS_TRY
 
    // Lame constants
    const double& rYoungModulus          = rProperties[YOUNG_MODULUS];
    const double& rPoissonCoefficient    = rProperties[POISSON_RATIO];
 
    rConstitutiveMatrix ( 0 , 0 ) = (rYoungModulus*(1.0-rPoissonCoefficient)/((1.0+rPoissonCoefficient)*(1.0-2*rPoissonCoefficient)));
    rConstitutiveMatrix ( 1 , 1 ) = rConstitutiveMatrix ( 0 , 0 );
    rConstitutiveMatrix ( 2 , 2 ) = rConstitutiveMatrix ( 0 , 0 );
 
    rConstitutiveMatrix ( 3 , 3 ) = rConstitutiveMatrix ( 0 , 0 )*(1-2*rPoissonCoefficient)/(2*(1.0-rPoissonCoefficient));
 
    rConstitutiveMatrix ( 0 , 1 ) = rConstitutiveMatrix ( 0 , 0 )*rPoissonCoefficient/(1.0-rPoissonCoefficient);
    rConstitutiveMatrix ( 1 , 0 ) = rConstitutiveMatrix ( 0 , 1 );
 
    rConstitutiveMatrix ( 0 , 2 ) = rConstitutiveMatrix ( 0 , 1 );
 
    rConstitutiveMatrix ( 2 , 0 ) = rConstitutiveMatrix ( 0 , 1 );
 
    rConstitutiveMatrix ( 1 , 2 ) = rConstitutiveMatrix ( 0 , 1 );
 
    rConstitutiveMatrix ( 2 , 1 ) = rConstitutiveMatrix ( 0 , 1 );
 
    KRATOS_CATCH(" ")
      }
 
 
 
 
 
 
 
 
 
 
 
    private:
 
 
 
 
 
 
 
 
 
 
 
 
      friend class Serializer;
 
      void save(Serializer& rSerializer) const override
      {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, SmallStrain3DLaw )
      }
 
      void load(Serializer& rSerializer) override
      {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, SmallStrain3DLaw )
      }
 
 
 
 
 
    }; // Class SmallStrainAxisymmetric2DLaw
 
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_SMALL_STRAIN_AXISYMMETRIC_2D_LAW_H_INCLUDED  defined