MohrCoulomb.h

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
// Copyright (c) 2021-2025 MPM-Geomechanics Development Team
 
#ifndef INC_MATERIALS_MOHRCOULOMB_H_
#define INC_MATERIALS_MOHRCOULOMB_H_
 
#include <Materials/ElasticJaumann.h>
 
class MohrCoulomb: public ElasticJaumann {
public:
 
    // softening struct
    class Softening {
 
    public :
 
        enum SofteningType { NONE, EXPONENTIAL };
 
        // softening type
        SofteningType softening_type;
 
        // exponential softening
        double exponential_shape_factor;
 
        // residual values
        double friction_residual;
        double cohesion_residual;
        double tensile_residual;
        double dilation_residual;
 
        // parameters to be affected by the softening
        bool friction_softening_active;
        bool cohesion_softening_active;
        bool tensile_softening_active;
        bool dilation_softening_active;
 
        inline double exponentialSoftening (double x, double eta, double y_i, double y_f) const {return y_f+(y_i-y_f)*std::exp(-eta*x);};
 
        Softening() {
 
            // softening type
            this->softening_type = MohrCoulomb::Softening::SofteningType::NONE;
 
            // exponential softening
            this->exponential_shape_factor = 0.0;
 
            // residual values
            this->friction_residual = 0.0;
            this->cohesion_residual = 0.0;
            this->tensile_residual = 0.0;
            this->dilation_residual = 0.0;
 
            // parameters to be affected by the softening
            this->friction_softening_active = false;
            this->cohesion_softening_active = false;
            this->tensile_softening_active = false;
            this->dilation_softening_active = false;
        }
    };
 
    MohrCoulomb(int id, double density, double young, double poisson, double friction, double cohesion, double dilation, double tensile, MohrCoulomb::Softening softening);
 
    virtual void updateStress(Particle* particle) const;
 
    void setCohesion(double c){this->cohesion=c;}
    
    double getCohesion() { return this->cohesion; }
    
    void setFriction(double f){this->friction=f;}
 
    void setSuVerticalStressFactor(double su){this->su_factor=su;}
 
    double getSuVerticalStressFactor(){ return this->su_factor;}
 
    virtual ~MohrCoulomb();
 
protected:
 
    double friction; 
    double cohesion; 
    double dilation; 
    double tensile; 
    double su_factor; 
    MohrCoulomb::Softening softening; 
};
 
#endif /* INC_MATERIALS_MOHRCOULOMB_H_ */