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.
contact_domain_LM_2D_condition.hpp
Go to the documentation of this file.
1 //
2 // Project Name: KratosContactMechanicsApplication $
3 // Created by: $Author: JMCarbonell $
4 // Last modified by: $Co-Author: $
5 // Date: $Date: July 2016 $
6 // Revision: $Revision: 0.0 $
7 //
8 //
9 
10 #if !defined(KRATOS_CONTACT_DOMAIN_LM_2D_CONDITION_H_INCLUDED )
11 #define KRATOS_CONTACT_DOMAIN_LM_2D_CONDITION_H_INCLUDED
12 
13 // System includes
14 
15 // External includes
16 
17 // Project includes
19 
20 namespace Kratos
21 {
30 
34 
38 
39 
40 class KRATOS_API(CONTACT_MECHANICS_APPLICATION) ContactDomainLM2DCondition
41  : public ContactDomainCondition
42 {
43 public:
44 
45 
51  typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
54 
56  typedef Node NodeType;
61 
70 
71 
74 
78 
79 
82 
83  ContactDomainLM2DCondition(IndexType NewId, GeometryType::Pointer pGeometry);
84 
85  ContactDomainLM2DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
86 
89 
90 
92  virtual ~ContactDomainLM2DCondition();
93 
97 
100 
101 
105 
113  Condition::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
114 
122  Condition::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
123 
124 
125  //************************************************************************************
126  //************************************************************************************
134  //int Check(const ProcessInfo& rCurrentProcessInfo);
135 
136  //std::string Info() const;
137 
141 
148 
150  // String Info() const override;
151 
153  // void PrintInfo(std::ostream& rOStream) const override;
154 
156  // void PrintData(std::ostream& rOStream) const override;
161 
162 protected:
171 
175  void SetMasterGeometry() override;
176 
177 
181  void CalculateContactFactor(const ProcessInfo& rCurrentProcessInfo) override;
182 
183 
187  void CalculatePreviousGap() override;
188 
189 
193  void CalculateExplicitFactors(ConditionVariables& rVariables,
194  const ProcessInfo& rCurrentProcessInfo) override;
198  void CalculateDomainShapeN(ConditionVariables& rVariables) override;
199 
200 
204  double& CalculateIntegrationWeight(double& rIntegrationWeight) override;
205 
209  void CalculateContactStiffness (double &Kcont,ConditionVariables& rVariables,
210  unsigned int& ndi,unsigned int& ndj,
211  unsigned int& idir,unsigned int& jdir) override;
212 
213 
217  void CalculateNormalForce (double &F,ConditionVariables& rVariables,
218  unsigned int& ndi,unsigned int& idir) override;
219 
223  void CalculateTangentStickForce (double &F,ConditionVariables& rVariables,
224  unsigned int& ndi,unsigned int& idir) override;
228  void CalculateTangentSlipForce (double &F,ConditionVariables& rVariables,
229  unsigned int& ndi,unsigned int& idir) override;
230 
234 
235  inline bool CheckFictiousContacts(ConditionVariables& rVariables);
236 
237  PointType& CalculateCurrentTangent(PointType &rTangent) override;
238 
239  void FSigmaP(ConditionVariables& rVariables, std::vector<Vector > &SigmaP, PointType& AuxVector,unsigned int &ndi,unsigned int &ndj,unsigned int &ndk,unsigned int &ndr);
240 
241  void FSigmaPnd(ConditionVariables& rVariables, std::vector<Vector > &SigmaP, PointType& AuxVector,unsigned int &ndi,unsigned int &ndj);
242 
243 
244 
251  friend class Serializer;
252 
253  void save(Serializer& rSerializer) const override;
254 
255  void load(Serializer& rSerializer) override;
256 
264 
265 private:
266 
291 
292 }; // Class ContactDomainLM2DCondition
293 
301 /* inline std::istream& operator >> (std::istream& rIStream,
302  ContactDomainLM2DCondition& rThis);
303 */
305 /* inline std::ostream& operator << (std::ostream& rOStream,
306  const ContactDomainLM2DCondition& rThis)
307  {
308  rThis.PrintInfo(rOStream);
309  rOStream << std::endl;
310  rThis.PrintData(rOStream);
311 
312  return rOStream;
313  }*/
315 
316 } // namespace Kratos.
317 #endif // KRATOS_CONTACT_DOMAIN_LM_2D_CONDITION_H_INCLUDED defined
PeriodicInterfaceProcess & operator=(const PeriodicInterfaceProcess &)=delete
Definition: constitutive_law.h:47
Definition: contact_domain_condition.hpp:49
Definition: contact_domain_LM_2D_condition.hpp:42
Node NodeType
NodeType.
Definition: contact_domain_LM_2D_condition.hpp:56
GeometryData::IntegrationMethod IntegrationMethod
Type definition for integration methods.
Definition: contact_domain_LM_2D_condition.hpp:53
ContactDomainUtilities::PointType PointType
Tensor order 1 definition.
Definition: contact_domain_LM_2D_condition.hpp:63
ContactDomainLM2DCondition()
Default constructors.
Definition: contact_domain_LM_2D_condition.hpp:81
ContactDomainUtilities::BaseLengths BaseLengths
BaseLengths.
Definition: contact_domain_LM_2D_condition.hpp:69
Element::ElementType ElementType
Element Type.
Definition: contact_domain_LM_2D_condition.hpp:60
ConstitutiveLaw ConstitutiveLawType
Definition: contact_domain_LM_2D_condition.hpp:49
Geometry< NodeType > GeometryType
Geometry Type.
Definition: contact_domain_LM_2D_condition.hpp:58
ContactDomainUtilities::SurfaceVector SurfaceVector
SurfaceVector.
Definition: contact_domain_LM_2D_condition.hpp:65
ConstitutiveLawType::Pointer ConstitutiveLawPointerType
Pointer type for constitutive laws.
Definition: contact_domain_LM_2D_condition.hpp:51
KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(ContactDomainLM2DCondition)
Counted pointer of ContactDomainLM2DCondition.
ContactDomainUtilities::SurfaceScalar SurfaceScalar
SurfaceScalar.
Definition: contact_domain_LM_2D_condition.hpp:67
Base class for all Elements.
Definition: element.h:60
IntegrationMethod
Definition: geometry_data.h:76
Geometry base class.
Definition: geometry.h:71
This class defines the node.
Definition: node.h:65
Point class.
Definition: point.h:59
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
std::size_t IndexType
The definition of the index type.
Definition: key_hash.h:35
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
ModelPart::NodesContainerType NodesArrayType
Definition: gid_gauss_point_container.h:42
F
Definition: hinsberg_optimization.py:144
def load(f)
Definition: ode_solve.py:307
Definition: contact_domain_utilities.hpp:74
Definition: contact_domain_utilities.hpp:91
Definition: contact_domain_utilities.hpp:83