Assembler.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 _ASSEMBLER_HPP_
#define _ASSEMBLER_HPP_

#include <vector>
#include <memory>
#include <Eigen/Dense>
#include <Eigen/SparseCore>

#include "Mesh.hpp"
#include "LoadCombo.hpp"

class Assembler{

    public:
        Assembler();

        ~Assembler();

        void SetMassTolerance(double tol);

        void SetForceTolerance(double tol);

        void SetStiffnessTolerance(double tol);

        void SetLoadCombination(std::shared_ptr<LoadCombo> &combo);

        Eigen::SparseMatrix<double> ComputeMassMatrix(std::shared_ptr<Mesh> &mesh); 

        Eigen::SparseMatrix<double> ComputeStiffnessMatrix(std::shared_ptr<Mesh> &mesh);

        Eigen::SparseMatrix<double> ComputeDampingMatrix(std::shared_ptr<Mesh> &mesh);

        Eigen::SparseMatrix<double> ComputePMLHistoryMatrix(std::shared_ptr<Mesh> &mesh);

        Eigen::VectorXd ComputeProgressiveForceVector(std::shared_ptr<Mesh> &mesh);

        Eigen::VectorXd ComputeInternalForceVector(std::shared_ptr<Mesh> &mesh);

        Eigen::VectorXd ComputeDynamicInternalForceVector(std::shared_ptr<Mesh> &mesh);

        Eigen::VectorXd ComputeExternalForceVector(std::shared_ptr<Mesh> &mesh, unsigned int k);

        Eigen::VectorXd ComputeSupportMotionIncrement(std::shared_ptr<Mesh> &mesh, unsigned int k);

        Eigen::VectorXd ComputePMLHistoryVector(std::shared_ptr<Mesh> &mesh);

    protected:
        void AddNodeInertiaForces(std::shared_ptr<Mesh> &mesh, Eigen::VectorXd &DynamicForces);

        void AddElementDynamicForces(std::shared_ptr<Mesh> &mesh, Eigen::VectorXd &DynamicForces);

        void AssembleNodalMass(std::shared_ptr<Mesh> &mesh, Eigen::SparseMatrix<double>& nodemass);

        void AssembleElementMass(std::shared_ptr<Mesh> &mesh, Eigen::SparseMatrix<double>& elemmass);

    private:
        double MassTolerance;

        double StiffnessTolerance;

        double ForceTolerance;

        std::shared_ptr<LoadCombo> LoadCombination;
};

#endif