EshelbyInclusionBase.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MoDELib is distributed without any warranty under the
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef model_EshelbyInclusionBase_H_
#define model_EshelbyInclusionBase_H_
 
#ifndef _MODEL_NON_SINGULAR_DD_
#define _MODEL_NON_SINGULAR_DD_ 1
#endif
 
#include <Eigen/Dense>
//#include <Material.h>
#include <StaticID.h>
 
// http://solidmechanics.org/text/Chapter5_4/Chapter5_4.htm
#include <SlipSystem.h>
#include <SecondPhase.h>
 
namespace model
{
    
    
    template <int dim>
    class EshelbyInclusionBase : public StaticID<EshelbyInclusionBase<dim>>
    {
        typedef Eigen::Matrix<double,dim,dim> MatrixDim;
        typedef Eigen::Matrix<double,dim,1>   VectorDim;
        
    public:
        const double nu;
        const double mu;
        const double lambda;
        const double mobilityReduction;
        const int phaseID;
        const std::shared_ptr<SecondPhase<dim>> secondPhase;
//        const double L;
//        const double M;
//        const double K;
 
//        VectorDim C;
//        double a;
//        double a2;
        MatrixDim eT;
        double eTNorm;
        MatrixDim pT;
        
        
    public:
        
//        static std::map<const GammaSurface*,GammaSurface> gammaSurfaceMap;
//
//        static void addSlipSystems(const std::vector<std::shared_ptr<SlipSystem>>& slipSystems);
        
        
//        double misfitEnergy(const Eigen::Matrix<double,3,1>& b, const GammaSurface* const matrixGammaSurface) const;
 
        
        /**********************************************************************/
        EshelbyInclusionBase(
                         const MatrixDim& _eT,
                         const double& _nu,
                         const double& _mu,
                         const double& _mobilityReduction,
                         const int& _phaseID,
                         const std::shared_ptr<SecondPhase<dim>>& sph);
  
        virtual ~EshelbyInclusionBase();
        virtual bool contains(const VectorDim& x) const =0 ;
        virtual MatrixDim stress(const VectorDim& x) const =0;
        
    };
    
    
//    template <int dim>
//    std::map<const GammaSurface*,GammaSurface> EshelbyInclusionBase<dim>::gammaSurfaceMap;
}
#endif