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.
axisym_updated_lagrangian_U_J_W_wP_element.hpp
Go to the documentation of this file.
1 //
2 // Project Name: KratosPfemSolidMechanicsApplication $
3 // Created by: $Author: PNavas $
4 // Last modified by: $Co-Author: LMonforte $
5 // Date: $Date: October 2017 $
6 // Revision: $Revision: -0.1 $
7 //
8 //
9 
10 #if !defined(KRATOS_AXISYM_UPDATED_LAGRANGIAN_U_J_W_WP_ELEMENT_H_INCLUDED )
11 #define KRATOS_AXISYM_UPDATED_LAGRANGIAN_U_J_W_WP_ELEMENT_H_INCLUDED
12 
13 // System includes
14 
15 // External includes
16 
17 // Project includes
19 
20 namespace Kratos
21 {
36 
38 
39 
40 class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) AxisymUpdatedLagrangianUJWwPElement
42 {
43 public:
44 
50  typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
55 
59 
62 
65 
67  AxisymUpdatedLagrangianUJWwPElement(IndexType NewId, GeometryType::Pointer pGeometry);
68 
69  AxisymUpdatedLagrangianUJWwPElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
70 
73 
74 
77 
81 
84 
85 
89 
100  Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
101 
109  Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
110 
111  //************* STARTING - ENDING METHODS
112 
113 
117  void GetDofList(DofsVectorType& rElementalDofList, ProcessInfo& rCurrentProcessInfo) override;
118 
122  void EquationIdVector(EquationIdVectorType& rResult, ProcessInfo& rCurrentProcessInfo) override;
123 
127  void GetValuesVector(Vector& rValues, int Step = 0) const override;
128 
132  void GetFirstDerivativesVector(Vector& rValues, int Step = 0) const override;
133 
137  void GetSecondDerivativesVector(Vector& rValues, int Step = 0) const override;
138 
139 
140  //************************************************************************************
141  //************************************************************************************
149  int Check(const ProcessInfo& rCurrentProcessInfo) override;
150 
151 
158  void CalculateMassMatrix(MatrixType& rMassMatrix,
159  ProcessInfo& rCurrentProcessInfo) override;
160 
167  void CalculateDampingMatrix(MatrixType& rDampingMatrix,
168  ProcessInfo& rCurrentProcessInfo) override;
169 
173 
184 protected:
190 
191 
192  /****
193  the time step (requiered). It shall be somewhere else.
194  ****/
195  double mTimeStep;
196 
200 
204 
205 
210  virtual void CalculateAndAddLHS(LocalSystemComponents& rLocalSystem,
211  ElementDataType& rVariables,
212  double& rIntegrationWeight) override;
213 
217  virtual void CalculateAndAddKWwP(MatrixType& rK,
218  ElementDataType & rVariables,
219  double& rIntegrationWeight
220  );
221 
225  virtual void CalculateAndAddKUwP(MatrixType& rK,
226  ElementDataType & rVariables,
227  double& rIntegrationWeight
228  );
232  virtual void CalculateAndAddKPPStab(MatrixType& rK,
233  ElementDataType & rVariables,
234  double& rIntegrationWeight
235  );
240  virtual void CalculateAndAddRHS(LocalSystemComponents& rLocalSystem,
241  ElementDataType& rVariables,
242  Vector& rVolumeForce,
243  double& rIntegrationWeight) override;
244 
248  virtual void CalculateAndAddInternalWaterForces(VectorType& rRightHandSideVector,
249  ElementDataType & rVariables,
250  double& rIntegrationWeight
251  );
252 
256  void CalculateAndAddExternalForcesUJWwP(VectorType& rRightHandSideVector,
257  ElementDataType & rVariables,
258  Vector & rVolumeForces,
259  double& rIntegrationWeight
260  );
264  void CalculateAndAddFluidLinearMomentum(VectorType& rRightHandSideVector,
265  ElementDataType & rVariables,
266  double& rIntegrationWeight
267  );
271  void CalculateAndAddMassBalanceEquation(VectorType& rRightHandSideVector,
272  ElementDataType & rVariables,
273  double& rIntegrationWeight
274  );
275 
279  virtual void CalculateAndAddStabilizationRHS(VectorType& rRightHandSideVector,
280  ElementDataType & rVariables,
281  double& rIntegrationWeight
282  );
286  virtual void CalculateAndAddMassStabilizationMatrix(MatrixType& rMassMatrix,
287  ElementDataType & rVariables,
288  double& rIntegrationWeight
289  );
290 
294  virtual void CalculateAndAddDampingStabilizationMatrix(MatrixType& rDampingMatrix,
295  ElementDataType & rVariables,
296  double& rIntegrationWeight
297  );
298 
302  virtual void CalculateAndAddHighOrderDampingMatrix(MatrixType& rDampingMatrix,
303  ElementDataType & rVariables,
304  double& rIntegrationWeight
305  );
306 
310  virtual void CalculateAndAddHighOrderKPP(MatrixType& rK,
311  ElementDataType & rVariables,
312  double& rIntegrationWeight
313  );
314 
318  virtual void CalculateAndAddHighOrderRHS(VectorType& rRightHandSideVector,
319  ElementDataType & rVariables,
320  double& rIntegrationWeight
321  );
322 
323 
327  virtual void InitializeElementData(ElementDataType & rVariables, const ProcessInfo& rCurrentProcessInfo) override;
328 
329 
330 
334  void InitializeSystemMatrices(MatrixType& rLeftHandSideMatrix,
335  VectorType& rRightHandSideVector,
336  Flags& rCalculationFlags) override;
337 
338 
342  double & CalculateStabilizationFactor( ElementDataType & rVariables, double & rStabilizationFactor);
343 
354 
355 private:
356 
362 
363 
367 
368 
372 
373 
378 
382  friend class Serializer;
383 
384  // A private default constructor necessary for serialization
385 
386  virtual void save(Serializer& rSerializer) const override;
387 
388  virtual void load(Serializer& rSerializer) override;
389 
390 
397 
398 
399 }; // Class UpdatedLagrangianUJWwPElement
400 
401 
402 
403 } // namespace Kratos
404 #endif // KRATOS_AXISYM_UPDATED_LAGRANGIAN_U_J_W_wP_ELEMENT_H_INCLUDED
405 
PeriodicInterfaceProcess & operator=(const PeriodicInterfaceProcess &)=delete
Definition: axisym_updated_lagrangian_U_J_element.hpp:40
Updated Lagrangian Large Displacement Lagrangian U-W Element for 3D and 2D geometries....
Definition: axisym_updated_lagrangian_U_J_W_wP_element.hpp:42
KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(AxisymUpdatedLagrangianUJWwPElement)
Counted pointer of LargeDisplacementUPElement.
ConstitutiveLawType::Pointer ConstitutiveLawPointerType
Pointer type for constitutive laws.
Definition: axisym_updated_lagrangian_U_J_W_wP_element.hpp:50
ConstitutiveLaw ConstitutiveLawType
Definition: axisym_updated_lagrangian_U_J_W_wP_element.hpp:48
double mTimeStep
Definition: axisym_updated_lagrangian_U_J_W_wP_element.hpp:195
GeometryData::IntegrationMethod IntegrationMethod
Type definition for integration methods.
Definition: axisym_updated_lagrangian_U_J_W_wP_element.hpp:54
ConstitutiveLawType::StressMeasure StressMeasureType
StressMeasure from constitutive laws.
Definition: axisym_updated_lagrangian_U_J_W_wP_element.hpp:52
Definition: constitutive_law.h:47
StressMeasure
Definition: constitutive_law.h:69
std::vector< DofType::Pointer > DofsVectorType
Definition: element.h:100
std::vector< std::size_t > EquationIdVectorType
Definition: element.h:98
Definition: flags.h:58
std::size_t IndexType
Definition: flags.h:74
IntegrationMethod
Definition: geometry_data.h:76
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