pfem_particle.h

Go to the documentation of this file.

/*
==============================================================================
Kratos
A General Purpose Software for Multi-Physics Finite Element Analysis
Version 1.0 (Released on march 05, 2007).
 
Copyright 2007
Pooyan Dadvand, Riccardo Rossi
pooyan@cimne.upc.edu
rrossi@cimne.upc.edu
CIMNE (International Center for Numerical Methods in Engineering),
Gran Capita' s/n, 08034 Barcelona, Spain
 
Permission is hereby granted, free  of charge, to any person obtaining
a  copy  of this  software  and  associated  documentation files  (the
"Software"), to  deal in  the Software without  restriction, including
without limitation  the rights to  use, copy, modify,  merge, publish,
distribute,  sublicense and/or  sell copies  of the  Software,  and to
permit persons to whom the Software  is furnished to do so, subject to
the following condition:
 
Distribution of this code for  any  commercial purpose  is permissible
ONLY BY DIRECT ARRANGEMENT WITH THE COPYRIGHT OWNER.
 
The  above  copyright  notice  and  this permission  notice  shall  be
included in all copies or substantial portions of the Software.
 
THE  SOFTWARE IS  PROVIDED  "AS  IS", WITHOUT  WARRANTY  OF ANY  KIND,
EXPRESS OR  IMPLIED, INCLUDING  BUT NOT LIMITED  TO THE  WARRANTIES OF
MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
IN NO EVENT  SHALL THE AUTHORS OR COPYRIGHT HOLDERS  BE LIABLE FOR ANY
CLAIM, DAMAGES OR  OTHER LIABILITY, WHETHER IN AN  ACTION OF CONTRACT,
TORT  OR OTHERWISE, ARISING  FROM, OUT  OF OR  IN CONNECTION  WITH THE
SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
 
==============================================================================
*/
 
//
//   Project Name:        Kratos
//   Last Modified by:    $Author: pablo $
//   Date:                $Date: 2013-03-05 11:34:12 $
//   Revision:            $Revision: 1.5 $
//
//
 
 
#if !defined(KRATOS_PFEM_PARTICLE_H_INCLUDED )
#define  KRATOS_PFEM_PARTICLE_H_INCLUDED
 
 
 
// System includes
#include <string>
#include <iostream>
#include <sstream>
#include <cstddef>
 
 
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "containers/array_1d.h"
#include "includes/serializer.h"
#include "geometries/point.h"
#include "includes/model_part.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
 
/*
@see Geometry
@see Node
@see IntegrationPoint
*/
//template<std::size_t TDimension, class TDataType = double> //always size 3!
class PFEM_Particle : public Point
{
public: 
 
    KRATOS_CLASS_POINTER_DEFINITION(PFEM_Particle);
    
    typedef double TDataType; 
 
    typedef array_1d<double,3> BaseType;
 
    typedef Point Type;
 
    typedef BaseType CoordinatesArrayType;
 
    typedef std::size_t SizeType;
 
    typedef std::size_t IndexType;
 
 
 
    PFEM_Particle(TDataType const& NewX, TDataType const& NewY, TDataType const& NewZ) : Point(NewX, NewY, NewZ)
    {
        this->ERASE_FLAG=true; //initializing as useless particle
        this->VELOCITY=ZeroVector(3);
        this->ACCELERATION=ZeroVector(3);
        this->DISTANCE=0.0;
        this->PRESSURE=0.0;
        this->OLD_PRESSURE=0.0;
        this->PLASTIC_PRESSURE=0.0;
        this->TEMPERATURE=0.0;
        this->SHEAR_MODULUS=0.0;
        this->BULK_MODULUS=0.0;
        this->COHESION=0.0;
        this->THETA=0.0;
        this->DENSITY=0.0;
        this->TOTAL_PLASTIC_DEFORMATION=ZeroVector(6);
        this->OLD_TOTAL_PLASTIC_DEFORMATION=ZeroVector(6);
        this->OLD_SIGMA=ZeroVector(6);
        this->SIGMA=ZeroVector(6);
        this->HAS_UPDATED_STRESSES=false;
        this->PLASTICIZED=false;
    }
    
    PFEM_Particle() : Point(0.0, 0.0, 0.0)
    {
        this->ERASE_FLAG=true;
        this->VELOCITY=ZeroVector(3);
        this->ACCELERATION=ZeroVector(3);
        this->DISTANCE=0.0;
        this->PRESSURE=0.0;
        this->OLD_PRESSURE=0.0;
        this->PLASTIC_PRESSURE=0.0;
        this->TEMPERATURE=0.0;
        this->SHEAR_MODULUS=0.0;
        this->BULK_MODULUS=0.0;
        this->COHESION=0.0;
        this->THETA=0.0;
        this->DENSITY=0.0;
        this->TOTAL_PLASTIC_DEFORMATION=ZeroVector(6);
        this->OLD_SIGMA=ZeroVector(6);
        this->SIGMA=ZeroVector(6);
        this->HAS_UPDATED_STRESSES=false;
        this->PLASTICIZED=false;
    }
    
    ~PFEM_Particle()
        {}
    //returning references
    array_1d<double,3>& GetVelocity()
    {
        return this->VELOCITY;
    }
        
    double& GetVelocity(const unsigned int i)
    {
        return this->VELOCITY[i];
    }
    
    //returning references
    array_1d<double,3>& GetAcceleration()
    {
        return this->ACCELERATION;
    }
    
    double& GetAcceleration(const unsigned int i)
    {
        return this->ACCELERATION[i];
    }
    
    double& GetDistance()
    {
        return this->DISTANCE;
    }
    
    double& GetShearModulus()
    {
        return this->SHEAR_MODULUS;
    }
    
    double& GetBulkModulus()
    {
        return this->BULK_MODULUS;
    }
    double& GetSecantBulkModulus()
    {
        return this->SECANT_BULK_MODULUS;
    }
    
    //returning references
    array_1d<double,6>& GetSigma() //should be 6 in 3d
    {
        return this->SIGMA;
    }
    
    double& GetSigma(const unsigned int i)
    {
        return this->SIGMA[i];
    }
 
    array_1d<double,6>& GetOldSigma()
    {
        return this->OLD_SIGMA;
    }
    double& GetOldSigma(const unsigned int i)
    {
        return this->OLD_SIGMA[i];
    }
    
    array_1d<double,6>& GetTotalPlasticDeformation()
    {
        return this->TOTAL_PLASTIC_DEFORMATION;
    }
    double& GetTotalPlasticDeformation(const unsigned int i)
    {
        return this->TOTAL_PLASTIC_DEFORMATION[i];
    }
    
    array_1d<double,6>& GetOldTotalPlasticDeformation()
    {
        return this->OLD_TOTAL_PLASTIC_DEFORMATION;
    }
    double& GetOldTotalPlasticDeformation(const unsigned int i)
    {
        return this->OLD_TOTAL_PLASTIC_DEFORMATION[i];
    }
    //double& GetTemperature()
    //{
    //  return this->TEMPERATURE;
    //}
    
    //double& GetOxygen()
    //{
    //  return this->OXYGEN;
    //}
    
    double& GetPressure()
    {
        return this->PRESSURE;
    }
    
    double& GetOldPressure()
    {
        return this->OLD_PRESSURE;
    }
    
    double& GetPlasticPressure()
    {
        return this->PLASTIC_PRESSURE;
    }
    
    double& GetTemperature()
    {
        return this->TEMPERATURE;
    }
    
    double& GetDensity()
    {
        return this->DENSITY;
    }
    double& GetTheta()
    {
        return this->THETA;
    }
    double& GetCohesion()
    {
        return this->COHESION;
    }
    
    /*
    double& GetGradientDiscontinuity()
    {
        return this->GRADIENT_DISCONTINUITY;
    }
    */
    //Element::Pointer& GetElement()
    //{
    //  return this->ELEMENT_WEAKPOINTER;
    //}
    /*
    unsigned int& GetElementId()
    {
        return this->ELEMENT_ID;
    }
    */
    bool& GetEraseFlag()
    {
        return this->ERASE_FLAG;
    }
    
    bool& HasUpdatedStresses()
    {
        return this->HAS_UPDATED_STRESSES;
    }
    
    bool& IsPlasticized()
    {
        return this->PLASTICIZED;
    }
    
    
    //Copy constructor. Initialize this point with the coordinates
    //of given particle.
    
    //PFEM_Particle(PFEM_Particle const& rOtherParticle)
     //   : Point(rOtherParticle) {}
      
       
private: 
    array_1d<double,3> VELOCITY;
    array_1d<double,3> ACCELERATION;
    double DISTANCE;
    double PRESSURE;
    double TEMPERATURE;
    double SHEAR_MODULUS;
    double BULK_MODULUS;
    double SECANT_BULK_MODULUS;
    double COHESION;
    double THETA;
    double DENSITY;
    double OLD_PRESSURE;
    double PLASTIC_PRESSURE;
    array_1d<double,6> OLD_TOTAL_PLASTIC_DEFORMATION;
    array_1d<double,6> TOTAL_PLASTIC_DEFORMATION;
    array_1d<double,6> OLD_SIGMA;
    array_1d<double,6> SIGMA; //should be of size 6 in 3d!
    //Element::Pointer ELEMENT_WEAKPOINTER;
    
    //unsigned int ELEMENT_ID;
    //double GRADIENT_DISCONTINUITY;
    bool ERASE_FLAG;
    bool HAS_UPDATED_STRESSES;
    bool PLASTICIZED;
    
 
 
};
 
} //namespace Kratos
 
#endif //KRATOS_PFEM_PARTICLE_H_INCLUDED defined