model_volume_calculation_process.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosDelaunayMeshingApplication $
//   Created by:          $Author:             JMCarbonell $
//   Last modified by:    $Co-Author:                      $
//   Date:                $Date:                April 2018 $
//   Revision:            $Revision:                   0.0 $
//
//
 
#if !defined(KRATOS_VOLUME_CALCULATION_PROCESS_H_INCLUDED )
#define  KRATOS_VOLUME_CALCULATION_PROCESS_H_INCLUDED
 
 
// External includes
 
// System includes
 
// Project includes
#include "includes/model_part.h"
 
//Data:
//StepData:
//Flags:    (checked)
//          (set)
//          (modified)
//          (reset)
 
 
namespace Kratos
{
 
 
class ModelVolumeCalculationProcess
  : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION( ModelVolumeCalculationProcess );
 
    typedef ModelPart::ConditionType         ConditionType;
    typedef ModelPart::PropertiesType       PropertiesType;
    typedef ConditionType::GeometryType       GeometryType;
 
    ModelVolumeCalculationProcess(ModelPart& rModelPart, bool Axisymmetric, int EchoLevel)
      : mrModelPart(rModelPart)
    {
      mAxisymmetric = Axisymmetric;
      mEchoLevel = EchoLevel;
    }
 
    ModelVolumeCalculationProcess(ModelPart& rModelPart)
      : mrModelPart(rModelPart)
    {
      mAxisymmetric = false;
    }
 
    virtual ~ModelVolumeCalculationProcess() {}
 
 
 
 
 
 
    void ExecuteInitializeSolutionStep() override
    {
 
      KRATOS_TRY
 
      mModelVolume = ComputeModelVolume(mMeshVolume);
 
      if( mEchoLevel > 0 )
    std::cout<<"  [ Model Volume : "<<mModelVolume<<" ]"<<std::endl;
 
      KRATOS_CATCH( "" )
 
    }
 
    void ExecuteFinalizeSolutionStep() override
    {
      KRATOS_TRY
 
      std::vector<double> MeshVolume;
      double ModelVolume = ComputeModelVolume(MeshVolume);
 
      if( mEchoLevel > 0 )
    std::cout<<"  [ Model Volume : "<<ModelVolume<<" ] [ Step increment : "<<ModelVolume-mModelVolume<<" ] "<<std::endl;
 
 
      KRATOS_CATCH( "" )
    }
 
 
 
 
 
 
    std::string Info() const override
    {
        return "ModelVolumeCalculationProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "ModelVolumeCalculationProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
 
private:
 
    ModelPart&  mrModelPart;
 
    bool mAxisymmetric;
 
    double mModelVolume;
 
    std::vector<double> mMeshVolume;
 
    int mEchoLevel;
 
 
    ModelVolumeCalculationProcess& operator=(ModelVolumeCalculationProcess const& rOther);
 
 
 
    double ComputeModelVolume(std::vector<double> & rMeshVolume)
    {
 
      KRATOS_TRY
 
      unsigned int start=0;
      unsigned int NumberOfMeshes=mrModelPart.NumberOfMeshes();
      if(NumberOfMeshes>1)
    start=1;
 
 
      rMeshVolume.resize(NumberOfMeshes+start);
      std::fill( rMeshVolume.begin(), rMeshVolume.end(), 0.0 );
 
      double ModelVolume = 0;
 
      if( mAxisymmetric ){
 
    //std::cout<<"  AXISYMMETRIC MODEL "<<std::endl;
 
    double radius = 0;
    double two_pi = 6.28318530717958647693;
 
    //By the way: set meshes options from bools
    for(unsigned int MeshId=start; MeshId<NumberOfMeshes; ++MeshId)
      {
        for(ModelPart::ElementsContainerType::const_iterator ie = mrModelPart.ElementsBegin(MeshId); ie != mrModelPart.ElementsEnd(MeshId); ++ie)
          {
 
        const unsigned int dimension = ie->GetGeometry().WorkingSpaceDimension();
 
        if( dimension > 2 )
          std::cout<<" Axisymmetric problem with dimension: "<<dimension<<std::endl;
 
        radius = 0;
        for( unsigned int i=0; i<ie->GetGeometry().size(); i++ )
          radius += ie->GetGeometry()[i].X();
 
        radius/=double(ie->GetGeometry().size());
 
        rMeshVolume[MeshId] += ie->GetGeometry().Area() * two_pi * radius ;
 
 
          }
 
        ModelVolume += rMeshVolume[MeshId];
      }
 
      }
      else{
 
    //By the way: set meshes options from bools
    for(unsigned int MeshId=start; MeshId<NumberOfMeshes; ++MeshId)
      {
        for(ModelPart::ElementsContainerType::const_iterator ie = mrModelPart.ElementsBegin(MeshId); ie != mrModelPart.ElementsEnd(MeshId); ++ie)
          {
        const unsigned int dimension = ie->GetGeometry().WorkingSpaceDimension();
        if( dimension == 2){
          rMeshVolume[MeshId] += ie->GetGeometry().Area() *  ie->GetProperties()[THICKNESS];
        }
        else if( dimension == 3){
          rMeshVolume[MeshId] += ie->GetGeometry().Volume();
        }
        else{
          //do nothing.
        }
 
          }
 
        ModelVolume += rMeshVolume[MeshId];
 
      }
 
      }
 
      return ModelVolume;
 
      KRATOS_CATCH( "" )
 
    }
 
 
    //Process(Process const& rOther);
 
 
 
}; // Class Process
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  ModelVolumeCalculationProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const ModelVolumeCalculationProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_MODEL_VOLUME_CALCULATION_PROCESS_H_INCLUDED  defined