KratosMultiphysics
KRATOS Multiphysics (Kratos) is a framework for building parallel, multi-disciplinary simulation software, aiming at modularity, extensibility, and high performance. Kratos is written in C++, and counts with an extensive Python interface.
large_displacement_element.hpp
Go to the documentation of this file.
1 //
2 // Project Name: KratosSolidMechanicsApplication $
3 // Created by: $Author: JMCarbonell $
4 // Last modified by: $Co-Author: $
5 // Date: $Date: July 2013 $
6 // Revision: $Revision: 0.0 $
7 //
8 //
9 
10 #if !defined(KRATOS_LARGE_DISPLACEMENT_ELEMENT_H_INCLUDED)
11 #define KRATOS_LARGE_DISPLACEMENT_ELEMENT_H_INCLUDED
12 
13 // System includes
14 
15 // External includes
16 
17 // Project includes
19 
20 namespace Kratos
21 {
36 
38 
44 class KRATOS_API(SOLID_MECHANICS_APPLICATION) LargeDisplacementElement
45  : public SolidElement
46 {
47 public:
48 
51 
55  typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
64 
67 
71 
74 
76  LargeDisplacementElement(IndexType NewId, GeometryType::Pointer pGeometry);
77 
78  LargeDisplacementElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
79 
82 
84  ~LargeDisplacementElement() override;
85 
89 
92 
96 
104  Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
105 
113  Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
114 
115 
119  void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
120 
121 
125  void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
126 
127 
128 
129  //on integration points:
130 
134  void CalculateOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
135 
139  void CalculateOnIntegrationPoints(const Variable<Vector>& rVariable, std::vector<Vector>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
140 
144  void CalculateOnIntegrationPoints(const Variable<Matrix >& rVariable, std::vector< Matrix >& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
145 
146 
147 
148  //************************************************************************************
149  //************************************************************************************
157  int Check(const ProcessInfo& rCurrentProcessInfo) const override;
158 
159 
163 
171  std::string Info() const override
172  {
173  std::stringstream buffer;
174  buffer << "Large Displacement Element #" << Id();
175  return buffer.str();
176  }
177 
179  void PrintInfo(std::ostream& rOStream) const override
180  {
181  rOStream << "Large Displacement Element #" << Id();
182  }
183 
185  void PrintData(std::ostream& rOStream) const override
186  {
187  GetGeometry().PrintData(rOStream);
188  }
193 
194 protected:
206 
210  void CalculateAndAddKuug(MatrixType& rLeftHandSideMatrix,
211  ElementDataType& rVariables,
212  double& rIntegrationWeight) override;
213 
217  void SetElementData(ElementDataType& rVariables,
219  const int & rPointNumber) override;
220 
224  virtual void GetHistoricalVariables(ElementDataType& rVariables,
225  const double& rPointNumber );
226 
230  virtual void CalculateGreenLagrangeStrain(const Matrix& rF,
231  Vector& rStrainVector);
232 
236  virtual void CalculateAlmansiStrain(const Matrix& rF,
237  Vector& rStrainVector);
238 
239 
250 
251 private:
252 
258 
259 
263 
264 
268 
269 
274 
278  friend class Serializer;
279 
280  // A private default constructor necessary for serialization
281 
282  void save(Serializer& rSerializer) const override;
283 
284  void load(Serializer& rSerializer) override;
285 
286 
293 
294 }; // Class LargeDisplacementElement
295 
303 
304 } // namespace Kratos.
305 #endif // KRATOS_LARGE_DISPLACEMENT_ELEMENT_H_INCLUDED defined
PeriodicInterfaceProcess & operator=(const PeriodicInterfaceProcess &)=delete
Definition: constitutive_law.h:47
StressMeasure
Definition: constitutive_law.h:69
std::size_t IndexType
Definition: flags.h:74
IntegrationMethod
Definition: geometry_data.h:76
std::size_t SizeType
Definition: geometry_data.h:173
Large Displacement Lagrangian Element for 3D and 2D geometries. (base class)
Definition: large_displacement_element.hpp:46
SolidElement::ElementDataType ElementDataType
Type for element variables.
Definition: large_displacement_element.hpp:63
void PrintData(std::ostream &rOStream) const override
Print object's data.
Definition: large_displacement_element.hpp:185
void PrintInfo(std::ostream &rOStream) const override
Print information about this object.
Definition: large_displacement_element.hpp:179
std::string Info() const override
Turn back information as a string.
Definition: large_displacement_element.hpp:171
KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(LargeDisplacementElement)
Counted pointer of LargeDisplacementElement.
GeometryData::IntegrationMethod IntegrationMethod
Type definition for integration methods.
Definition: large_displacement_element.hpp:59
ConstitutiveLawType::Pointer ConstitutiveLawPointerType
Pointer type for constitutive laws.
Definition: large_displacement_element.hpp:55
ConstitutiveLawType::StressMeasure StressMeasureType
StressMeasure from constitutive laws.
Definition: large_displacement_element.hpp:57
ConstitutiveLaw ConstitutiveLawType
Reference type definition for constitutive laws.
Definition: large_displacement_element.hpp:53
GeometryData::SizeType SizeType
Type for size.
Definition: large_displacement_element.hpp:61
PointerVector is a container like stl vector but using a vector to store pointers to its data.
Definition: pointer_vector.h:72
ProcessInfo holds the current value of different solution parameters.
Definition: process_info.h:59
The serialization consists in storing the state of an object into a storage format like data file or ...
Definition: serializer.h:123
Large Displacement Lagrangian Element for 3D and 2D geometries. (base class)
Definition: solid_element.hpp:49
Matrix MatrixType
Definition: geometrical_transformation_utilities.h:55
Modeler::Pointer Create(const std::string &ModelerName, Model &rModel, const Parameters ModelParameters)
Checks if the modeler is registered.
Definition: modeler_factory.cpp:30
pybind11::list CalculateOnIntegrationPoints(TObject &dummy, const Variable< TDataType > &rVariable, const ProcessInfo &rProcessInfo)
Definition: add_mesh_to_python.cpp:142
void InitializeSolutionStep(ConstructionUtility &rThisUtil, std::string ThermalSubModelPartName, std::string MechanicalSubModelPartName, std::string HeatFluxSubModelPartName, std::string HydraulicPressureSubModelPartName, bool thermal_conditions, bool mechanical_conditions, int phase)
Definition: add_custom_utilities_to_python.cpp:45
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21
def load(f)
Definition: ode_solve.py:307
Definition: constitutive_law.h:189
Definition: solid_element.hpp:83