small_strain_plane_stress_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_PLANE_STRESS_2D_LAW_H_INCLUDED )
#define  KRATOS_SMALL_STRAIN_PLANE_STRESS_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) SmallStrainPlaneStress2DLaw : public SmallStrain3DLaw
    {
    public:
 
      KRATOS_CLASS_POINTER_DEFINITION(SmallStrainPlaneStress2DLaw);
 
 
      SmallStrainPlaneStress2DLaw() : SmallStrain3DLaw() {}
 
      SmallStrainPlaneStress2DLaw(ModelTypePointer pModel) : SmallStrain3DLaw(pModel) {};
 
      SmallStrainPlaneStress2DLaw(const SmallStrainPlaneStress2DLaw& rOther) : SmallStrain3DLaw(rOther) {}
 
      SmallStrainPlaneStress2DLaw& operator=(SmallStrainPlaneStress2DLaw const& rOther)
      {
    SmallStrain3DLaw::operator=(rOther);
    return *this;
      }
 
      ConstitutiveLaw::Pointer Clone() const override
      {
        return Kratos::make_shared<SmallStrainPlaneStress2DLaw>(*this);
      }
 
      ~SmallStrainPlaneStress2DLaw() override{}
 
 
 
      SizeType WorkingSpaceDimension() override { return 2; }
 
      SizeType GetStrainSize() const override { return 3; }
 
      void GetLawFeatures(Features& rFeatures) override
      {
    KRATOS_TRY
 
        //Set the type of law
    rFeatures.mOptions.Set( PLANE_STRESS_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(" ")
      }
 
 
      void CalculateMaterialResponseKirchhoff(Parameters& rValues) override
      {
    KRATOS_TRY
 
    //0.- Check if the constitutive parameters are passed correctly to the law calculation
    //CheckParameters(rValues);
 
        const Flags& rOptions = rValues.GetOptions();
 
    const Properties& rProperties  = rValues.GetMaterialProperties();
 
    Vector& rStrainVector                  = rValues.GetStrainVector();
    Vector& rStressVector                  = rValues.GetStressVector();
 
    // Calculate total Kirchhoff stress
 
    if( rOptions.Is( ConstitutiveLaw::COMPUTE_STRESS ) ){
 
      Matrix& rConstitutiveMatrix = rValues.GetConstitutiveMatrix();
 
      this->CalculateConstitutiveMatrix( rConstitutiveMatrix, rProperties);
 
      this->CalculateStress( rStrainVector, rConstitutiveMatrix, rStressVector );
 
    }
    else if( rOptions.Is( ConstitutiveLaw::COMPUTE_CONSTITUTIVE_TENSOR ) ){
 
      Matrix& rConstitutiveMatrix  = rValues.GetConstitutiveMatrix();
      this->CalculateConstitutiveMatrix(rConstitutiveMatrix, rProperties);
 
    }
 
    // std::cout<<" StrainVector "<<rValues.GetStrainVector()<<std::endl;
    // std::cout<<" StressVector "<<rValues.GetStressVector()<<std::endl;
    // std::cout<<" ConstitutiveMatrix "<<rValues.GetConstitutiveMatrix()<<std::endl;
 
    KRATOS_CATCH(" ")
 
      }
 
 
 
 
 
 
      std::string Info() const override
      {
    std::stringstream buffer;
        buffer << "SmallStrainPlaneStress2DLaw" ;
        return buffer.str();
      }
 
      void PrintInfo(std::ostream& rOStream) const override {rOStream << "SmallStrainPlaneStress2DLaw";}
 
      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];
 
    // Plane stress constitutive matrix
    rConstitutiveMatrix ( 0 , 0 ) = (rYoungModulus)/(1.0-rPoissonCoefficient*rPoissonCoefficient);
    rConstitutiveMatrix ( 1 , 1 ) = rConstitutiveMatrix ( 0 , 0 );
 
    rConstitutiveMatrix ( 2 , 2 ) = rConstitutiveMatrix ( 0 , 0 )*(1.0-rPoissonCoefficient)*0.5;
 
    rConstitutiveMatrix ( 0 , 1 ) = rConstitutiveMatrix ( 0 , 0 )*rPoissonCoefficient;
    rConstitutiveMatrix ( 1 , 0 ) = rConstitutiveMatrix ( 0 , 1 );
 
    rConstitutiveMatrix ( 0 , 2 ) = 0.0;
    rConstitutiveMatrix ( 2 , 0 ) = 0.0;
    rConstitutiveMatrix ( 1 , 2 ) = 0.0;
    rConstitutiveMatrix ( 2 , 1 ) = 0.0;
 
    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 SmallStrainPlaneStress2DLaw
 
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_SMALL_STRAIN_PLANE_STRESS_2D_LAW_H_INCLUDED  defined