moving_load_condition.h

Go to the documentation of this file.

// KRATOS  ___|  |                   |                   |
//       \___ \  __|  __| |   |  __| __| |   |  __| _` | |
//             | |   |    |   | (    |   |   | |   (   | |
//       _____/ \__|_|   \__,_|\___|\__|\__,_|_|  \__,_|_| MECHANICS
//
//  License:         BSD License
//                   license: StructuralMechanicsApplication/license.txt
//
//  Main authors:    Aron Noordam
//
 
#pragma once 
 
// System includes
 
// System includes
 
// External includes
 
// Project includes
#include "custom_conditions/base_load_condition.h"
 
namespace Kratos
{
 
 
 
 
 
 
template<std::size_t TDim, std::size_t TNumNodes>
class KRATOS_API(STRUCTURAL_MECHANICS_APPLICATION)  MovingLoadCondition
    : public BaseLoadCondition
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(MovingLoadCondition );
 
 
    MovingLoadCondition(
        IndexType NewId,
        GeometryType::Pointer pGeometry
        );
 
    MovingLoadCondition(
        IndexType NewId,
        GeometryType::Pointer pGeometry,
        PropertiesType::Pointer pProperties
        );
 
    ~MovingLoadCondition() override;
 
 
 
 
    Condition::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        PropertiesType::Pointer pProperties
        ) const override;
 
    Condition::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        PropertiesType::Pointer pProperties
        ) const override;
 
    Condition::Pointer Clone (
        IndexType NewId,
        NodesArrayType const& ThisNodes
        ) const override;
 
 
    bool HasRotDof() const override;
 
 
    void InitializeSolutionStep(const ProcessInfo & rCurrentProcessInfo) override;
 
    void InitializeNonLinearIteration(const ProcessInfo & rCurrentProcessInfo) override;
 
 
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "Point load Condition #" << Id();
        return buffer.str();
    }
 
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "Point load Condition #" << Id();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        pGetGeometry()->PrintData(rOStream);
    }
 
 
 
 
protected:
 
 
 
 
 
 
 
    void CalculateAll(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo,
        const bool CalculateStiffnessMatrixFlag,
        const bool CalculateResidualVectorFlag
        ) override;
 
    void CalculateExactNormalShapeFunctions(VectorType& rShapeFunctionsVector, const double LocalXCoord) const;
 
    void CalculateExactShearShapeFunctions(VectorType& rShapeFunctionsVector, const double LocalXCoord) const;
 
    void CalculateExactRotationalShapeFunctions(VectorType& rShapeFunctionsVector, const double LocalXCoord) const;
 
    void CalculateExactShearShapeFunctionsDerivatives(VectorType& rShapeFunctionsVector, const double LocalXCoord) const;
 
    void CalculateExactRotationalShapeFunctionsDerivatives(VectorType& rShapeFunctionsVector, const double LocalXCoord) const;
 
    void CalculateRotationMatrix(BoundedMatrix<double, TDim, TDim>& rRotationMatrix, const GeometryType& rGeom);
 
 
    Matrix CalculateGlobalMomentMatrix(const VectorType& rRotationalShapeFunctionVector, const array_1d<double, TDim>& rLocalMovingLoad) const;
 
 
 
 
 
 
    // A protected default constructor necessary for serialization
    MovingLoadCondition() {};
 
 
private:
 
 
    bool mIsMovingLoad = false;
 
 
 
    void GetRotationsVector(Vector& rRotationsVector, const int Step) const;
 
 
 
    Vector CalculateLoadPointDisplacementVector();
 
    Vector CalculateLoadPointRotationVector();
 
 
 
 
 
    friend class Serializer;
 
    void save( Serializer& rSerializer ) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, BaseLoadCondition );
        rSerializer.save("mIsMovingLoad", mIsMovingLoad);
    }
 
    void load( Serializer& rSerializer ) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, BaseLoadCondition );
        rSerializer.load("mIsMovingLoad", mIsMovingLoad);
    }
 
    
 
}; // Class MovingLoadCondition
 
 
 
 
 
template<std::size_t TDim, std::size_t TNumNodes>
inline std::istream& operator >> (std::istream& rIStream,
        MovingLoadCondition<TDim, TNumNodes>& rThis);
 
template<std::size_t TDim, std::size_t TNumNodes>
inline std::ostream& operator << (std::ostream& rOStream,
        const MovingLoadCondition<TDim, TNumNodes>& rThis)
{
  rThis.PrintInfo(rOStream);
  rOStream << std::endl;
  rThis.PrintData(rOStream);
 
  return rOStream;
}
 
 
}  // namespace Kratos.