modified_cam_clay_yield_criterion.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:        BSD License
//                  Kratos default license: kratos/license.txt
//
//  Main authors:    Bodhinanda Chandra
//
 
 
#if !defined(KRATOS_MODIFIED_CAM_CLAY_YIELD_CRITERION_H_INCLUDED )
#define      KRATOS_MODIFIED_CAM_CLAY_YIELD_CRITERION_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_constitutive/yield_criteria/particle_yield_criterion.hpp"
#include "custom_constitutive/hardening_laws/cam_clay_hardening_law.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(PARTICLE_MECHANICS_APPLICATION) ModifiedCamClayYieldCriterion
    : public ParticleYieldCriterion
{
    public:
 
        KRATOS_CLASS_POINTER_DEFINITION( ModifiedCamClayYieldCriterion );
 
 
        ModifiedCamClayYieldCriterion();
 
        ModifiedCamClayYieldCriterion(ModifiedCamClayYieldCriterion const& rOther);
 
        ModifiedCamClayYieldCriterion(HardeningLawPointer pHardeningLaw);
 
        ModifiedCamClayYieldCriterion& operator=(ModifiedCamClayYieldCriterion const& rOther);
 
 
        virtual ~ModifiedCamClayYieldCriterion();
 
 
 
 
 
 
        /*
        * @brief This function return the Modified cam clay yield criterion at the given principal stress condition
        * @param[in/out] rStateFunction Mohr coulomb yield criterion
        * @param[in] rStressVector Principal stresses
        * @param[in] rAlpha Plastic volumetric strain
        * @param[in] rOldPreconsolidationPressure The value of Preconsolidation Stress at the previous time step
        * @return Modified cam clay yield criterion
        */
        double& CalculateYieldCondition(double & rStateFunction, const Vector& rStressVector, const double& rAlpha, const double& rOldPreconsolidationPressure, const Properties& rProp) override;
 
 
        /*
        * @brief This function return the first derivative of modified cam clay yield criterion at the given principal stress condition
        * @param[in] rStressVector Principal stresses
        * @param[in/out] rFirstDerivative First stress derivative value of yield function
        * @param[in] rAlpha Plastic volumetric strain
        * @param[in] rOldPreconsolidationPressure The value of Preconsolidation Stress at the previous time step
        */
        void CalculateYieldFunctionDerivative(const Vector& rStressVector, Vector& rFirstDerivative, const double& rAlpha, const double& rOldPreconsolidationPressure, const Properties& rProp) override;
 
 
        /*
        * @brief This function return the second derivative of modified cam clay yield criterion at the given principal stress condition
        * @param[in] rStressVector Principal stresses
        * @param[in/out] rSecondDerivative Second stress derivative value of yield function
        */
        void CalculateYieldFunctionSecondDerivative(const Vector& rStressVector, Vector& rSecondDerivative, const Properties& rProp) 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:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
    private:
 
 
 
 
 
 
 
 
 
    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_MODIFIED_CAM_CLAY_YIELD_CRITERION_H_INCLUDED  defined