rve_periodicity_utility.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Riccardo Rossi
//
 
#if !defined (KRATOS_RVE_PERIODICITY_UTILITY_H_INCLUDED)
#define KRATOS_RVE_PERIODICITY_UTILITY_H_INCLUDED
 
// System includes
 
// External includes
#include <tuple>
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "utilities/binbased_fast_point_locator_conditions.h"
 
namespace Kratos
{
 
class KRATOS_API(KRATOS_CORE) RVEPeriodicityUtility
{
  public:
 
    KRATOS_CLASS_POINTER_DEFINITION(RVEPeriodicityUtility);
 
    typedef std::size_t IndexType;
 
    typedef std::tuple<std::vector<IndexType>, std::vector<double>, Vector> DataTupletype;
 
    typedef Dof<double> DofType;
 
    typedef std::vector< DofType::Pointer > DofPointerVectorType;
 
    typedef Node NodeType;
 
    typedef Variable<double> DoubleVariableType;
 
 
    RVEPeriodicityUtility(
        ModelPart& rDestinationModelPart,
        std::size_t EchoLevel = 0
        ) : mrModelPart(rDestinationModelPart)
          , mEchoLevel(EchoLevel)
    {
    }
 
    RVEPeriodicityUtility(RVEPeriodicityUtility const& rOther) = delete;
 
    virtual ~RVEPeriodicityUtility() = default;
 
 
    RVEPeriodicityUtility &operator=(RVEPeriodicityUtility const& rOther) = delete;
 
 
    void AssignPeriodicity(ModelPart& rMasterModelPart,
                           ModelPart& rSlaveModelPart,
                           const Matrix& rStrainTensor,
                           const Vector& rDirection,
                           const double SearchTolerance = 1.0e-6
                           );
 
    void Finalize(const Variable<array_1d<double, 3>>& rVariable);
 
 
    virtual std::string Info() const
    {
        std::stringstream buffer;
        buffer << "RVEPeriodicityUtility";
        return buffer.str();
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const
    {
        rOStream << "RVEPeriodicityUtility";
    }
 
    virtual void PrintData(std::ostream& rOStream) const {}
 
 
  protected:
 
    void AppendIdsAndWeights(
        std::map<IndexType, DataTupletype>& rAux,
        const IndexType MasterId,
        const double MasterWeight,
        std::vector<IndexType>& rFinalMastersIds,
        std::vector<double>& rFinalMastersWeights,
        Vector& rFinalT
        );
 
 
private:
 
    ModelPart& mrModelPart; 
 
    std::map<IndexType, DataTupletype> mAuxPairings; 
 
    std::size_t mEchoLevel = 0; 
 
 
    MasterSlaveConstraint::Pointer  GenerateConstraint(
        IndexType& rConstraintId,
        const DoubleVariableType& rVar,
        NodeType::Pointer pSlaveNode,
        const std::vector<IndexType>& rMasterIds,
        const Matrix& rRelationMatrix,
        const Vector& rTranslationVector
        );
 
    template<std::size_t TDim>
    void AuxiliaryAssignPeriodicity(
        ModelPart& rMasterModelPart,
        ModelPart& rSlaveModelPart,
        const Matrix& rStrainTensor,
        const Vector& rDirection,
        const double SearchTolerance = 1.0e-6);
 
 
}; // Class RVEPeriodicityUtility
 
 
inline std::istream &operator>>(std::istream& rIStream,
                                RVEPeriodicityUtility& rThis)
{
    return rIStream;
}
 
inline std::ostream &operator<<(std::ostream& rOStream,
                                const RVEPeriodicityUtility& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
} // namespace Kratos.
 
#endif // KRATOS_RVE_PERIODICITY_UTILITY_H_INCLUDED defined