Node.hpp

Go to the documentation of this file.

//==============================================================================
//
//                       Seismo Virtual Laboratory
//             Module for Serial and Parallel Analysis of seismic 
//         wave propagation and soil-structure interaction simulation
//         Copyright (C) 2018-2021, The California Institute of Technology
//                         All Rights Reserved.
//
// Commercial use of this program without express permission of the California
// Institute of Technology, is strictly  prohibited. See  file "COPYRIGHT"  in
// main  directory  for  information on  usage  and  redistribution, and for a
// DISCLAIMER OF ALL WARRANTIES.
//
//==============================================================================
//
// Written by:
//   Danilo S. Kusanovic (dkusanov@caltech.edu)
//   Elnaz E. Seylabi    (elnaze@unr.edu)
//
// Supervised by:
//   Domniki M. Asimaki  (domniki@caltech.edu)
//
// References : 
//   [1]
//
// Description:
//------------------------------------------------------------------------------

#ifndef _NODE_HPP_
#define _NODE_HPP_

#include <map>
#include <vector>
#include <Eigen/Dense>

class Node{

    public:
        Node(unsigned int numberDofs, Eigen::VectorXd coordinates, bool isFixed);

        ~Node();

        bool IsFixed() const;

        void InitialState();

        void SetMass(Eigen::VectorXd &mass);

        void SetAsFixed(bool isFixed);

        void SetReaction(Eigen::VectorXd &reaction);

        void SetProgressiveForces(Eigen::VectorXd &force);

        void SetCoordinates(Eigen::VectorXd &coordinates);

        void SetFreeDegreeOfFreedom(std::vector<int> &freeDOF);

        void SetTotalDegreeOfFreedom(std::vector<int> &totalDOF);

        void SetDisplacements(Eigen::VectorXd &Uo);

        void SetVelocities(Eigen::VectorXd &Vo);

        void SetAccelerations(Eigen::VectorXd &Ao);

        void SetIncrementalDisplacements(Eigen::VectorXd &dU);

        void SetDomainReductionMotion(Eigen::MatrixXd &Uo);

        void SetPMLVector(Eigen::VectorXd &Ub);

        void SetSupportMotion(unsigned int k, std::vector<double> &Uo);

        void DelSupportMotion();

        Eigen::VectorXd GetMass() const;

        Eigen::VectorXd GetReaction() const;

        const Eigen::VectorXd& GetCoordinates() const;

        const Eigen::VectorXd& GetDisplacements() const;

        const Eigen::VectorXd& GetVelocities() const;

        const Eigen::VectorXd& GetAccelerations() const;

        Eigen::VectorXd GetInertialForces() const;

        const Eigen::VectorXd& GetProgressiveForces() const;

        const Eigen::VectorXd& GetPMLVector() const;

        const Eigen::VectorXd& GetIncrementalDisplacements() const;

        Eigen::VectorXd GetDomainReductionMotion(unsigned int k) const;

        Eigen::VectorXd GetSupportMotion(unsigned int k);

        unsigned int GetNumberOfSupportMotion();

        unsigned int GetNumberOfDegreeOfFreedom() const;

        const std::vector<int>& GetFreeDegreeOfFreedom() const;

        const std::vector<int>& GetTotalDegreeOfFreedom() const;

    private:
        //The node condition.
        bool Fixed;

        unsigned int NumberOfDegreeOfFreedom;

        Eigen::VectorXd Coordinates;

        Eigen::VectorXd Displacements;

        Eigen::VectorXd Velocities;

        Eigen::VectorXd Accelerations;

        Eigen::VectorXd IncrementalDisplacements;

        Eigen::MatrixXd DomainReductionMotion;

        Eigen::VectorXd Mass;

        Eigen::VectorXd Reaction;

        Eigen::VectorXd CumulatedForce;

        Eigen::VectorXd PMLIntegratedVector;

        std::vector<int> FreeDegreeOfFreedom;

        std::vector<int> TotalDegreeOfFreedom;

        std::map<unsigned int, std::vector<double> > SupportMotion;
};

#endif