Material.h

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
// Copyright (c) 2021-2025 MPM-Geomechanics Development Team
 
#ifndef MATERIAL_H_
#define MATERIAL_H_
 
#include "Eigen/Core"
using Eigen::Vector3d;
 
class Particle;
 
class Material {
 
public:
 
    enum MaterialType{ NONE, ELASTIC, MOHRCOULOMB, ELASTOPLASTIC, NEWTONIAN };
 
    Material(int id=-1, double density=0, double mu = 0, MaterialType type = MaterialType::NONE);
    
    virtual ~Material();
    
    inline int getId() const { return this->id; }
    
    inline double getDensity() const { return this->density; }
 
    inline double getDensityFluid() const { return this->densityFluid; }
 
    inline double getPorosity() const { return this->porosity; }
 
    inline Vector3d getHydraulicConductivity() const { return this->hydraulicConductivity; }
 
    inline double getBulkModulusFluid() const { return this->bulkModulusFluid; }
 
    inline double getFrictionCoefficient() const { return this->mu; }
 
    inline void setId(int material_id) { this->id=material_id; }
    
    inline void setDensity(double material_density) { this->density=material_density; }
    
    inline void setDensityFluid(double fluid_density) { this->densityFluid=fluid_density; }
 
    inline void setPorosity(double material_porosity) { this->porosity=material_porosity; }
 
    inline void setBulkModulusFluid(double bulk_fluid) { this->bulkModulusFluid=bulk_fluid; }
 
    inline void setFrictionCoefficient(double _mu) { this->mu = _mu; }
 
    inline void setHydraulicConductivity(Vector3d hydraulic_conductivity) { this->hydraulicConductivity=hydraulic_conductivity; }
 
    inline void setType(MaterialType material_type) { this->type=material_type; }
    
    inline virtual MaterialType getType() const { return this->type; }
 
    virtual void updateStress(Particle* particle) const = 0;
 
    virtual double getSoundSpeed( ) const = 0;
 
protected:
    
    int id; 
 
    double density; 
 
    double densityFluid; 
 
    double porosity; 
 
    double bulkModulusFluid; 
 
    double mu;  
 
    Vector3d hydraulicConductivity; 
    
    MaterialType type; 
};
 
inline Material::Material(int id, double density, double mu, MaterialType type) {
    
    this->id=id;
    this->density=density;
    this->densityFluid=0.0;
    this->porosity=0.0;
    this->bulkModulusFluid=0.0;
    this->mu = mu;
    this->hydraulicConductivity.setZero();
    this->type=type;
}
 
inline Material::~Material() { }
 
#endif /* MATERIAL_H_ */