save_structure_conditions_process.h

Go to the documentation of this file.

//
//   Project Name:        Kratos
//   Last Modified by:    $Author: anonymous $
//   Date:                $Date: 2008-10-02 15:38:16 $
//   Revision:            $Revision: 1.1 $
//
//  this process save structural elements in a separate list
 
#if !defined(KRATOS_SAVE_STRUCTURE_CONDITIONS_PROCESS_INCLUDED )
#define  KRATOS_SAVE_STRUCTURE_CONDITIONS_PROCESS_INCLUDED
 
 
 
// System includes
#include <string>
#include <iostream>
#include <algorithm>
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "processes/process.h"
#include "includes/node.h"
#include "includes/condition.h"
#include "includes/model_part.h"
//#include "custom_elements/updated_lagrangian_fluid.h"
//#include "custom_elements/updated_lagrangian_fluid3D.h"
//#include "custom_elements/updated_lagrangian_fluid_inc.h"
//#include "custom_elements/updated_lagrangian_fluid3D_inc.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
 
 
class SaveStructureConditionsProcess
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(SaveStructureConditionsProcess);
 
 
    SaveStructureConditionsProcess()
    {
    }
 
    ~SaveStructureConditionsProcess() override
    {
    }
 
 
 
//      void operator()()
//      {
//          SaveStructure();
//      }
 
 
 
    void  SaveStructureConditions(ModelPart& fluid_model_part, ModelPart& structure_model_part)
    {
        KRATOS_TRY
 
        //number of structure nodes
        KRATOS_WATCH("SAVING STRUCTURE CONDITION")
        for(ModelPart::ConditionsContainerType::iterator ic = fluid_model_part.ConditionsBegin() ;
                ic != fluid_model_part.ConditionsEnd() ; ++ic)
        {
            //PointerVector<Element> struct_elements_list;
            //check number of structure nodes
            unsigned int n_struct=0;
            //unsigned int n_fluid = 0;
            for (unsigned int i=0; i<ic->GetGeometry().size(); i++)
            {
                n_struct += int( ic->GetGeometry()[i].FastGetSolutionStepValue(IS_STRUCTURE) );
                //n_fluid += int( ic->GetGeometry()[i].FastGetSolutionStepValue(IS_FLUID) );
 
            }
 
            if(n_struct==ic->GetGeometry().size())// && n_fluid != im->GetGeometry().size())
            {
                structure_model_part.Conditions().push_back(*(ic.base()));
                /*
                for (unsigned int i=0;i<ic->GetGeometry().size();i++)
                {
                ic->GetGeometry()[i].FastGetSolutionStepValue(IS_FLUID)=0.0;
                }
                */
            }
 
 
        }
        KRATOS_CATCH("")
    }
 
 
 
 
 
 
 
    std::string Info() const override
    {
        return "SaveStructureConditionsProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "SaveStructureConditionsProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    //ModelPart& mr_fluid_model_part;
    //ModelPart& mr_structure_model_part;
 
 
 
 
 
 
 
 
 
 
//      SaveStructureConditionsProcess& operator=(SaveStructureConditionsProcess const& rOther);
 
//      SaveStructureConditionsProcess(SaveStructureConditionsProcess const& rOther);
 
 
 
}; // Class SaveStructureConditionsProcess
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  SaveStructureConditionsProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const SaveStructureConditionsProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_SAVE_STRUCTURE_CONDITIONS_PROCESS_INCLUDED  defined