Seismo-VLAB  1.3
An Open-Source Finite Element Software for Meso-Scale Simulations
lin2DFrame2.hpp
Go to the documentation of this file.
1 //==============================================================================
2 //
3 // Seismo Virtual Laboratory
4 // Module for Serial and Parallel Analysis of seismic
5 // wave propagation and soil-structure interaction simulation
6 // Copyright (C) 2018-2021, The California Institute of Technology
7 // All Rights Reserved.
8 //
9 // Commercial use of this program without express permission of the California
10 // Institute of Technology, is strictly prohibited. See file "COPYRIGHT" in
11 // main directory for information on usage and redistribution, and for a
12 // DISCLAIMER OF ALL WARRANTIES.
13 //
14 //==============================================================================
15 //
16 // Written by:
17 // Danilo S. Kusanovic (dkusanov@caltech.edu)
18 // Elnaz E. Seylabi (elnaze@unr.edu)
19 //
20 // Supervised by:
21 // Domniki M. Asimaki (domniki@caltech.edu)
22 //
23 // References :
24 // [1] Finite Element Procedures, Bathe, K.J., Chapter 5: pages 400-408.
25 // Prentice-Hall, 1996.
26 //
27 // Description:
30 //------------------------------------------------------------------------------
31 
32 #ifndef _LIN2DFRAME2_HPP_
33 #define _LIN2DFRAME2_HPP_
34 
35 #include <map>
36 #include <memory>
37 #include <string>
38 #include <Eigen/Dense>
39 
40 #include "Node.hpp"
41 #include "Load.hpp"
42 #include "Section.hpp"
43 #include "Element.hpp"
44 #include "Damping.hpp"
45 #include "QuadratureRule.hpp"
46 
54 class lin2DFrame2 : public Element{
55 
56  public:
65  lin2DFrame2(const std::vector<unsigned int> nodes, std::unique_ptr<Section> &section, bool formulation=false, const std::string quadrature="GAUSS", unsigned int nGauss=3);
66 
68  ~lin2DFrame2();
69 
72  void CommitState();
73 
76  void ReverseState();
77 
80  void InitialState();
81 
84  void UpdateState();
85 
90  void SetDomain(std::map<unsigned int, std::shared_ptr<Node> > &nodes);
91 
95  void SetDamping(const std::shared_ptr<Damping> &damping);
96 
99  std::vector<unsigned int> GetTotalDegreeOfFreedom() const;
100 
104  Eigen::MatrixXd GetStrain() const;
105 
109  Eigen::MatrixXd GetStress() const;
110 
114  Eigen::MatrixXd GetStrainRate() const;
115 
121  Eigen::MatrixXd GetStrainAt(double x3, double x2) const;
122 
128  Eigen::MatrixXd GetStressAt(double x3, double x2) const;
129 
134  Eigen::VectorXd GetVTKResponse(std::string response) const;
135 
138  double ComputeEnergy();
139 
144  Eigen::MatrixXd ComputeMassMatrix();
145 
150  Eigen::MatrixXd ComputeStiffnessMatrix();
151 
156  Eigen::MatrixXd ComputeDampingMatrix();
157 
162  Eigen::MatrixXd ComputePMLMatrix();
163 
168  Eigen::VectorXd ComputeInternalForces();
169 
174  Eigen::VectorXd ComputeInternalDynamicForces();
175 
182  Eigen::VectorXd ComputeSurfaceForces(const std::shared_ptr<Load> &surface, unsigned int face);
183 
190  Eigen::VectorXd ComputeBodyForces(const std::shared_ptr<Load> &body, unsigned int k=0);
191 
198  Eigen::VectorXd ComputeDomainReductionForces(const std::shared_ptr<Load> &drm, unsigned int k);
199 
200  private:
202  double L;
203 
206 
208  double Phi;
209 
211  std::shared_ptr<Damping> theDamping;
212 
214  std::vector<std::shared_ptr<Node> > theNodes;
215 
217  std::vector<std::unique_ptr<Section> > theSection;
218 
220  std::unique_ptr<QuadratureRule> QuadraturePoints;
221 
223  double ComputeLength() const;
224 
228  Eigen::MatrixXd ComputeLocalAxes() const;
229 
233  Eigen::VectorXd ComputeStrain(const Eigen::MatrixXd &Bij) const;
234 
238  Eigen::VectorXd ComputeStrainRate(const Eigen::MatrixXd &Bij) const;
239 
243  Eigen::MatrixXd ComputeShapeFunctionMatrix(const double ri) const;
244 
248  Eigen::MatrixXd ComputeStrainDisplacementMatrix(const double ri) const;
249 
253  Eigen::MatrixXd ComputeInitialStiffnessMatrix() const;
254 };
255 
256 #endif
void CommitState()
Save the material states in the element.
~lin2DFrame2()
Destroys this lin2DFrame2 object.
double ComputeLength() const
Compute the current length of the element.
Eigen::MatrixXd ComputeStiffnessMatrix()
Compute the stiffness matrix of the element using gauss-integration.
Eigen::MatrixXd ComputeStrainDisplacementMatrix(const double ri) const
Evaluates the strain-displacement matrix at a given Gauss point.
This file contains the "Load" class declarations, which defines a load that can act in a node or an e...
Eigen::MatrixXd GetStrainRate() const
Gets the material/section (generalised) strain-rate.
std::vector< std::shared_ptr< Node > > theNodes
The Element&#39;s Nodes.
Definition: lin2DFrame2.hpp:214
Eigen::MatrixXd GetStress() const
Gets the material/section (generalised) stress.
void SetDomain(std::map< unsigned int, std::shared_ptr< Node > > &nodes)
Sets the finite element dependance among objects.
Eigen::VectorXd ComputeStrainRate(const Eigen::MatrixXd &Bij) const
Update strain rate in the element.
Eigen::MatrixXd ComputeLocalAxes() const
Compute/update the local axis-1 of the element.
Eigen::VectorXd ComputeStrain(const Eigen::MatrixXd &Bij) const
Update strain in the element.
double L
Length of the element.
Definition: lin2DFrame2.hpp:202
Eigen::VectorXd ComputeInternalDynamicForces()
Compute the elastic, inertial, and viscous forces acting on the element.
This file contains the abstract "Section object" declarations, which computes the strain...
void ReverseState()
Reverse the material/section states to previous converged state in this element.
std::vector< unsigned int > GetTotalDegreeOfFreedom() const
Gets the list of total-degree of freedom of this Element.
Eigen::MatrixXd ComputeShapeFunctionMatrix(const double ri) const
Evaluates the shape function matrix at a given Gauss point.
std::unique_ptr< QuadratureRule > QuadraturePoints
Coordinate of Gauss points.
Definition: lin2DFrame2.hpp:220
std::vector< std::unique_ptr< Section > > theSection
The Element&#39;s Section.
Definition: lin2DFrame2.hpp:217
void SetDamping(const std::shared_ptr< Damping > &damping)
Sets the damping model.
Eigen::MatrixXd GetStrainAt(double x3, double x2) const
Gets the material strain in section at coordinate (x3,x2).
Class for creating a 2D linearized two-node frame element in a mesh.
Definition: lin2DFrame2.hpp:54
std::shared_ptr< Damping > theDamping
The Damping model.
Definition: lin2DFrame2.hpp:211
Virtual class for creating an element in a mesh.
Definition: Element.hpp:51
lin2DFrame2(const std::vector< unsigned int > nodes, std::unique_ptr< Section > &section, bool formulation=false, const std::string quadrature="GAUSS", unsigned int nGauss=3)
Creates a lin2DFrame2 in a finite element Mesh.
Eigen::VectorXd ComputeInternalForces()
Compute the internal (elastic) forces acting on the element.
This file contains the "Node object" declarations, which is stores the coordinates, state variables, degrees-of-freedom, and total-degree of freedom lists of a node in a finite element mesh.
This file contains the "Damping object" declarations to update damping matrix to account for Caughey-...
Eigen::MatrixXd ComputePMLMatrix()
Compute the PML history matrix using gauss-integration.
Eigen::MatrixXd ComputeMassMatrix()
Compute the lumped/consistent mass matrix of the element.
Eigen::MatrixXd GetStrain() const
Gets the material/section (generalised) strain.
Eigen::MatrixXd ComputeInitialStiffnessMatrix() const
Compute the initial stiffness matrix of the element.
Eigen::VectorXd GetVTKResponse(std::string response) const
Gets the element internal response in VTK format for Paraview display.
void InitialState()
Brings the material/section state to its initial state in this element.
Eigen::MatrixXd ComputeDampingMatrix()
Compute the damping matrix of the element using gauss-integration.
bool Formulation
The element formulation.
Definition: lin2DFrame2.hpp:205
void UpdateState()
Update the material states in the element.
double ComputeEnergy()
Computes the element energy for a given deformation.
Eigen::VectorXd ComputeSurfaceForces(const std::shared_ptr< Load > &surface, unsigned int face)
Compute the surface forces acting on the element.
Eigen::MatrixXd GetStressAt(double x3, double x2) const
Gets the material stress in section at coordinate (x3,x2).
Eigen::VectorXd ComputeDomainReductionForces(const std::shared_ptr< Load > &drm, unsigned int k)
Compute the domain reduction forces acting on the element.
double Phi
Section Shear Factor.
Definition: lin2DFrame2.hpp:208
Eigen::VectorXd ComputeBodyForces(const std::shared_ptr< Load > &body, unsigned int k=0)
Compute the body forces acting on the element.
This file contains the abstract "Quadrature Rule" declarations to integrate quantities in the iso-par...
This file contains the "Element" object declarations, which defines an element in a finite element me...