linear_associative_plasticity_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_LINEAR_ASSOCIATIVE_PLASTICITY_MODEL_H_INCLUDED )
#define  KRATOS_LINEAR_ASSOCIATIVE_PLASTICITY_MODEL_H_INCLUDED
 
 
// System includes
 
// External includes
 
// Project includes
#include "custom_models/plasticity_models/non_linear_associative_plasticity_model.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
  template<class TElasticityModel, class TYieldSurface>
  class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) LinearAssociativePlasticityModel : public NonLinearAssociativePlasticityModel<TElasticityModel,TYieldSurface>
  {
  public:
 
 
    //elasticity model
    typedef TElasticityModel                               ElasticityModelType;
    typedef typename  ElasticityModelType::Pointer      ElasticityModelPointer;
 
    //yield surface
    typedef TYieldSurface                                     YieldSurfaceType;
    typedef typename YieldSurfaceType::Pointer             YieldSurfacePointer;
 
    //derived type
    typedef NonLinearAssociativePlasticityModel<ElasticityModelType,YieldSurfaceType>   DerivedType;
 
    //base type
    typedef PlasticityModel<ElasticityModelType,YieldSurfaceType>     BaseType;
 
    //common types
    typedef typename BaseType::Pointer                         BaseTypePointer;
    typedef typename BaseType::SizeType                               SizeType;
    typedef typename BaseType::VoigtIndexType                   VoigtIndexType;
    typedef typename BaseType::MatrixType                           MatrixType;
    typedef typename BaseType::ModelDataType                     ModelDataType;
    typedef typename BaseType::MaterialDataType               MaterialDataType;
    typedef typename BaseType::PlasticDataType                 PlasticDataType;
    typedef typename BaseType::InternalVariablesType     InternalVariablesType;
 
    KRATOS_CLASS_POINTER_DEFINITION( LinearAssociativePlasticityModel );
 
 
    LinearAssociativePlasticityModel() : DerivedType() {}
 
    LinearAssociativePlasticityModel(ElasticityModelPointer pElasticityModel, YieldSurfacePointer pYieldSurface) : DerivedType(pElasticityModel, pYieldSurface) {}
 
    LinearAssociativePlasticityModel(LinearAssociativePlasticityModel const& rOther) : DerivedType(rOther) {}
 
    LinearAssociativePlasticityModel& operator=(LinearAssociativePlasticityModel const& rOther)
    {
      DerivedType::operator=(rOther);
      return *this;
    }
 
    virtual ConstitutiveModel::Pointer Clone() const override
    {
      return Kratos::make_shared<LinearAssociativePlasticityModel>(*this);
    }
 
    virtual ~LinearAssociativePlasticityModel() {}
 
 
 
 
 
 
 
 
 
 
 
    virtual std::string Info() const override
    {
      std::stringstream buffer;
      buffer << "LinearAssociativePlasticityModel" ;
      return buffer.str();
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const override
    {
      rOStream << "LinearAssociativePlasticityModel";
    }
 
    virtual void PrintData(std::ostream& rOStream) const override
    {
      rOStream << "LinearAssociativePlasticityModel Data";
    }
 
 
 
 
  protected:
 
 
 
 
 
 
    // calculate return mapping
    virtual bool CalculateReturnMapping(PlasticDataType& rVariables, MatrixType& rStressMatrix)
    {
      KRATOS_TRY
 
      //start
      double DeltaStateFunction = 0;
 
      double& rEquivalentPlasticStrain     = rVariables.Internal.Variables[0];
      double& rDeltaGamma                  = rVariables.DeltaInternal.Variables[0];
 
 
      double StateFunction                = rVariables.TrialStateFunction;
 
      //Calculate Delta State Function:
      DeltaStateFunction = this->mYieldSurface.CalculateDeltaStateFunction( rVariables, DeltaStateFunction );
 
      //Calculate DeltaGamma:
      rDeltaGamma  = StateFunction/DeltaStateFunction;
 
      //Update Equivalent Plastic Strain:
      DeltaPlasticStrain        = sqrt(2.0/3.0) * rDeltaGamma;
      rEquivalentPlasticStrain += DeltaPlasticStrain;
 
      KRATOS_CATCH(" ")
    }
 
 
 
 
 
 
 
 
 
  private:
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const override
    {
      KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, DerivedType )
    }
 
    virtual void load(Serializer& rSerializer) override
    {
      KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, DerivedType )
    }
 
 
 
 
 
  }; // Class LinearAssociativePlasticityModel
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_LINEAR_ASSOCIATIVE_PLASTICITY_MODEL_H_INCLUDED  defined