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_U_P_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_U_P_ELEMENT_H_INCLUDED)
11 #define KRATOS_LARGE_DISPLACEMENT_U_P_ELEMENT_H_INCLUDED
12 
13 // System includes
14 
15 // External includes
16 
17 // Project includes
19 
20 
21 namespace Kratos
22 {
37 
39 
45 class KRATOS_API(SOLID_MECHANICS_APPLICATION) LargeDisplacementUPElement
47 {
48 public:
49 
55  typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
64 
68 
71 
74 
76  LargeDisplacementUPElement(IndexType NewId, GeometryType::Pointer pGeometry);
77 
78  LargeDisplacementUPElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
79 
82 
83 
85  ~LargeDisplacementUPElement() override;
86 
90 
93 
97 
108  Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
109 
117  Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
118 
119  //************* GETTING METHODS
120 
124  void GetDofList(DofsVectorType& rElementalDofList, const ProcessInfo& rCurrentProcessInfo) const override;
125 
129  void EquationIdVector(EquationIdVectorType& rResult, const ProcessInfo& rCurrentProcessInfo) const override;
130 
134  void GetValuesVector(Vector& rValues, int Step = 0) const override;
135 
139  void GetFirstDerivativesVector(Vector& rValues, int Step = 0) const override;
140 
144  void GetSecondDerivativesVector(Vector& rValues, int Step = 0) const override;
145 
146 
147  //************* COMPUTING METHODS
148 
155  void CalculateMassMatrix(MatrixType& rMassMatrix, const ProcessInfo& rCurrentProcessInfo) override;
156 
163  void CalculateDampingMatrix(MatrixType& rDampingMatrix, const ProcessInfo& rCurrentProcessInfo) override;
164 
165 
166 
167  //************************************************************************************
168  //************************************************************************************
176  int Check(const ProcessInfo& rCurrentProcessInfo) const override;
177 
181 
192 
193 protected:
199 
203 
207 
208 
213  void CalculateAndAddLHS(LocalSystemComponents& rLocalSystem,
214  ElementDataType& rVariables,
215  double& rIntegrationWeight) override;
216 
221  void CalculateAndAddRHS(LocalSystemComponents& rLocalSystem,
222  ElementDataType& rVariables,
223  Vector& rVolumeForce,
224  double& rIntegrationWeight) override;
225 
230  void CalculateAndAddDynamicLHS(MatrixType& rLeftHandSideMatrix,
231  ElementDataType& rVariables,
232  const ProcessInfo& rCurrentProcessInfo,
233  double& rIntegrationWeight) override;
234 
239  void CalculateAndAddDynamicRHS(VectorType& rRightHandSideVector,
240  ElementDataType& rVariables,
241  const ProcessInfo& rCurrentProcessInfo,
242  double& rIntegrationWeight) override;
243 
247  void CalculateAndAddKuum(MatrixType& rK,
248  ElementDataType & rVariables,
249  double& rIntegrationWeight
250  ) override;
251 
255  void CalculateAndAddKuug(MatrixType& rK,
256  ElementDataType & rVariables,
257  double& rIntegrationWeight
258  ) override;
259 
263  virtual void CalculateAndAddKup (MatrixType& rK,
264  ElementDataType & rVariables,
265  double& rIntegrationWeight
266  );
267 
271  virtual void CalculateAndAddKpu(MatrixType& rK,
272  ElementDataType & rVariables,
273  double& rIntegrationWeight
274  );
275 
276 
280  virtual void CalculateAndAddKpp(MatrixType& rK,
281  ElementDataType & rVariables,
282  double& rIntegrationWeight
283  );
284 
285 
289  virtual void CalculateAndAddKppStab(MatrixType& rK,
290  ElementDataType & rVariables,
291  double& rIntegrationWeight
292  );
293 
294 
295 
299  void CalculateAndAddExternalForces(VectorType& rRightHandSideVector,
300  ElementDataType& rVariables,
301  Vector& rVolumeForce,
302  double& rIntegrationWeight
303  ) override;
304 
305 
309  void CalculateAndAddInternalForces(VectorType& rRightHandSideVector,
310  ElementDataType & rVariables,
311  double& rIntegrationWeight
312  ) override;
313 
314 
318  virtual void CalculateAndAddPressureForces(VectorType& rRightHandSideVector,
319  ElementDataType & rVariables,
320  double& rIntegrationWeight
321  );
322 
323 
327  virtual void CalculateAndAddStabilizedPressure(VectorType& rRightHandSideVector,
328  ElementDataType & rVariables,
329  double& rIntegrationWeight);
330 
334  SizeType GetDofsSize() const override;
335 
336 
340  virtual double& CalculatePUCoefficient(double& rCoefficient, ElementDataType & rVariables);
341 
345  virtual double& CalculatePUDeltaCoefficient(double& rCoefficient, ElementDataType & rVariables);
346 
357 
358 private:
359 
365 
366 
370 
371 
375 
376 
381 
385  friend class Serializer;
386 
387  // A private default constructor necessary for serialization
388 
389  void save(Serializer& rSerializer) const override;
390 
391  void load(Serializer& rSerializer) override;
392 
393 
400 
401 }; // Class LargeDisplacementUPElement
402 
410 
411 } // namespace Kratos.
412 #endif // KRATOS_LARGE_DISPLACEMENT_U_P_ELEMENT_H_INCLUDED defined
PeriodicInterfaceProcess & operator=(const PeriodicInterfaceProcess &)=delete
Definition: constitutive_law.h:47
StressMeasure
Definition: constitutive_law.h:69
std::size_t SizeType
Definition: element.h:94
std::vector< DofType::Pointer > DofsVectorType
Definition: element.h:100
std::vector< std::size_t > EquationIdVectorType
Definition: element.h:98
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
Large Displacement Lagrangian U-P Element for 3D and 2D geometries. Linear Triangles and Tetrahedra (...
Definition: large_displacement_U_P_element.hpp:47
GeometryData::SizeType SizeType
Type for size.
Definition: large_displacement_U_P_element.hpp:61
GeometryData::IntegrationMethod IntegrationMethod
Type definition for integration methods.
Definition: large_displacement_U_P_element.hpp:59
LargeDisplacementElement::ElementDataType ElementDataType
Type for element variables.
Definition: large_displacement_U_P_element.hpp:63
ConstitutiveLaw ConstitutiveLawType
Definition: large_displacement_U_P_element.hpp:53
ConstitutiveLawType::StressMeasure StressMeasureType
StressMeasure from constitutive laws.
Definition: large_displacement_U_P_element.hpp:57
KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(LargeDisplacementUPElement)
Counted pointer of LargeDisplacementUPElement.
ConstitutiveLawType::Pointer ConstitutiveLawPointerType
Pointer type for constitutive laws.
Definition: large_displacement_U_P_element.hpp:55
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
Modeler::Pointer Create(const std::string &ModelerName, Model &rModel, const Parameters ModelParameters)
Checks if the modeler is registered.
Definition: modeler_factory.cpp:30
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21
def load(f)
Definition: ode_solve.py:307
Definition: solid_element.hpp:83
Definition: solid_element.hpp:233