truss_element.h

Go to the documentation of this file.

//  KRATOS  _____________
//         /  _/ ____/   |
//         / // / __/ /| |
//       _/ // /_/ / ___ |
//      /___/\____/_/  |_| Application
 
 
#if !defined(KRATOS_TRUSS_ELEMENT_H_INCLUDED)
#define KRATOS_TRUSS_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/element.h"
 
// Application includes
#include "iga_application_variables.h"
 
namespace Kratos
{
 
class TrussElement
    : public Element
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( TrussElement );
 
    typedef std::size_t SizeType;
    typedef std::size_t IndexType;
 
 
    TrussElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry)
        : Element(NewId, pGeometry)
    {};
 
    TrussElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry,
        PropertiesType::Pointer pProperties)
        : Element(NewId, pGeometry, pProperties)
    {};
 
    TrussElement()
        : Element()
    {};
 
    ~TrussElement() override = default;
 
 
    Element::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        PropertiesType::Pointer pProperties
        ) const override
    {
        return Kratos::make_intrusive<TrussElement>(
            NewId, pGeom, pProperties);
    };
 
    Element::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        PropertiesType::Pointer pProperties
        ) const override
    {
        return Kratos::make_intrusive<TrussElement>(
            NewId, GetGeometry().Create(ThisNodes), pProperties);
    };
 
 
    void GetDofList(
        DofsVectorType& rElementalDofList,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
 
    void Initialize(const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override
    {
        const SizeType nb_nodes = GetGeometry().size();
        const SizeType nb_dofs = nb_nodes * 3;
 
        if (rRightHandSideVector.size() != nb_dofs)
            rRightHandSideVector.resize(nb_dofs);
        noalias(rRightHandSideVector) = ZeroVector(nb_dofs);
 
        MatrixType left_hand_side_matrix;
 
        CalculateAll(left_hand_side_matrix, rRightHandSideVector,
            rCurrentProcessInfo, false, true);
    }
 
    void CalculateLeftHandSide(
        MatrixType& rLeftHandSideMatrix,
        const ProcessInfo& rCurrentProcessInfo) override
    {
        const SizeType nb_nodes = GetGeometry().size();
        const SizeType nb_dofs = nb_nodes * 3;
 
        VectorType right_hand_side_vector;
 
        if (rLeftHandSideMatrix.size1() != nb_dofs && rLeftHandSideMatrix.size2() != nb_dofs)
            rLeftHandSideMatrix.resize(nb_dofs, nb_dofs);
        noalias(rLeftHandSideMatrix) = ZeroMatrix(nb_dofs, nb_dofs);
 
        CalculateAll(rLeftHandSideMatrix, right_hand_side_vector,
            rCurrentProcessInfo, true, false);
    }
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override
    {
        const SizeType nb_nodes = GetGeometry().size();
        const SizeType nb_dofs = nb_nodes * 3;
 
        if (rRightHandSideVector.size() != nb_dofs)
            rRightHandSideVector.resize(nb_dofs);
        noalias(rRightHandSideVector) = ZeroVector(nb_dofs);
 
        if (rLeftHandSideMatrix.size1() != nb_dofs && rLeftHandSideMatrix.size2() != nb_dofs)
            rLeftHandSideMatrix.resize(nb_dofs, nb_dofs);
        noalias(rLeftHandSideMatrix) = ZeroMatrix(nb_dofs, nb_dofs);
 
        CalculateAll(rLeftHandSideMatrix, rRightHandSideVector,
            rCurrentProcessInfo, true, true);
    }
 
    void CalculateAll(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo,
        const bool ComputeLeftHandSide,
        const bool ComputeRightHandSide);
 
    void FinalizeSolutionStep(
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateBodyForces(Vector& rBodyForces);
 
    bool HasSelfWeight() const;
 
 
    void AddExplicitContribution(
        const VectorType& rRHSVector,
        const Variable<VectorType>& rRHSVariable,
        const Variable<double >& rDestinationVariable,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void AddExplicitContribution(
        const VectorType& rRHSVector, const Variable<VectorType>& rRHSVariable,
        const Variable<array_1d<double, 3>>& rDestinationVariable,
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateDampingMatrix(
        MatrixType& rDampingMatrix,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateMassMatrix(
        MatrixType& rMassMatrix,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLumpedMassVector(
        VectorType &rLumpedMassVector,
        const ProcessInfo &rCurrentProcessInfo) const override;
 
    void GetValuesVector(
        Vector& rValues,
        int Step = 0) const override;
 
    void GetFirstDerivativesVector(
        Vector& rValues,
        int Step = 0) const override;
 
    void GetSecondDerivativesVector(
        Vector& rValues,
        int Step = 0) const override;
 
 
    void CalculateOnIntegrationPoints(
        const Variable<double>& rVariable,
        std::vector<double>& rOutput,
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
    void PrintInfo(std::ostream& rOStream) const override {
        rOStream << "\"TrussElement\" #" << Id()
            << " with geometry #" << this->GetGeometry().Id() << " with center in: "
            << this->GetGeometry().Center() << std::endl;
    }
 
private:
 
    std::vector<array_1d<double, 3>> mReferenceBaseVector;
 
    std::vector<ConstitutiveLaw::Pointer> mConstitutiveLawVector;
 
 
    void InitializeMaterial();
 
    array_1d<double, 3> CalculateActualBaseVector(IndexType IntegrationPointIndex) const;
 
    void CalculateGreenLagrangeStrain(
        std::vector<double>& rGreenLagrangeVector) const;
 
    void CalculateTangentModulus(
        std::vector<double>& rTangentModulusVector,
        const ProcessInfo& rCurrentProcessInfo);
 
    double CalculatePrestressPK2(double reference_a, double actual_a) const;
 
    void CalculateStressPK2(
        std::vector<double>& rStressVector,
        const ProcessInfo& rCurrentProcessInfo) const;
 
    void CalculateStressCauchy(
        std::vector<double>& rStressVector,
        const ProcessInfo& rCurrentProcessInfo) const;
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, TrussElement);
        rSerializer.save("ReferenceBaseVector", mReferenceBaseVector);
        rSerializer.save("ConstitutiveLawVector", mConstitutiveLawVector);
    }
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, TrussElement);
        rSerializer.load("ReferenceBaseVector", mReferenceBaseVector);
        rSerializer.load("ConstitutiveLawVector", mConstitutiveLawVector);
    }
 
};
 
} // namespace Kratos
 
#endif // !defined(KRATOS_TRUSS_ELEMENT_H_INCLUDED)