cam_clay_yield_criterion.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosPfemSolidMechanicsApplication $
//   Created by:          $Author:                  LMonforte $
//   Last modified by:    $Co-Author:                         $
//   Date:                $Date:                    July 2015 $
//   Revision:            $Revision:                      0.0 $
//
//
 
#if !defined(KRATOS_CAM_CLAY_YIELD_CRITERION_H_INCLUDED )
#define      KRATOS_CAM_CLAY_YIELD_CRITERION_H_INCLUDED
 
 
 
// System includes
 
// External includes
 
// Project includes
#include "custom_constitutive/custom_yield_criteria/yield_criterion.hpp"
#include "custom_constitutive/custom_hardening_laws/cam_clay_hardening_law.hpp"
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) CamClayYieldCriterion
    : public YieldCriterion 
{
    public:
 
        KRATOS_CLASS_POINTER_DEFINITION( CamClayYieldCriterion );
 
 
        CamClayYieldCriterion();
 
        CamClayYieldCriterion(CamClayYieldCriterion const& rOther);
 
        CamClayYieldCriterion(HardeningLawPointer pHardeningLaw);
 
        CamClayYieldCriterion& operator=(CamClayYieldCriterion const& rOther);
 
 
        virtual ~CamClayYieldCriterion();
 
 
 
 
 
        double& CalculateYieldCondition(double & rStateFunction, const Vector& rStressVector, const double& rAlpha) override;
 
    void CalculateYieldFunctionDerivative(const Vector& rStressVector, Vector& rFirstDerivative, const double& rAlpha) override;
        
 
 
 
 
        // /// Turn back information as a string.
        // std::string Info() const override;
 
        // /// Print information about this object.
        // void PrintInfo(std::ostream& rOStream) const override;
 
        // /// Print object's data.
        // void PrintData(std::ostream& rOStream) const override;
 
 
 
 
 
    protected:
 
 
    
    
        //double EvaluateThirdInvariantEffect(const Vector& rStressVector);
        double EvaluateThirdInvariantEffect( const double & rLodeAngle );
 
        void CalculateAndAddThirdInvDerivative(const Vector& rStressVector, Vector& rYieldFunctionD);
 
 
 
 
 
 
 
 
 
    private:
 
 
 
 
 
 
 
        double GetPI();
 
        double GetSmoothingLodeAngle();
 
        void GetSmoothingConstants(double& rA, double& rB, const double& rLode);
 
 
    
    friend class Serializer;
 
    // A private default constructor necessary for serialization
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
 
 
 
    }; // Class MisesHuberYieldCriterion
 
 
 
 
 
 
    // /// input stream function
    // inline std::istream& operator >> (std::istream& rIStream,
    //                                   MisesHuberYieldCriterion& rThis);
 
    // /// output stream function
    // inline std::ostream& operator << (std::ostream& rOStream,
    //                                   const MisesHuberYieldCriterion& rThis)
    // {
    //     rThis.PrintInfo(rOStream);
    //     rOStream << std::endl;
    //     rThis.PrintData(rOStream);
 
    //     return rOStream;
    // }
 
 
}  // namespace Kratos.
 
#endif // KRATOS_CAM_CLAY_YIELD_CRITERION_H_INCLUDED  defined