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.
non_linear_hencky_plastic_axisym_2D_law.hpp
Go to the documentation of this file.
1 //
2 // Project Name: KratosPfemSolidMechanicsApplication $
3 // Created by: $Author: LMonforte $
4 // Last modified by: $Co-Author: $
5 // Date: $Date: July 2015 $
6 // Revision: $Revision: 0.0 $
7 //
8 //
9 
10 #if !defined (KRATOS_NON_LINEAR_HENCKY_PLASTIC_AXISYM_2D_LAW_H_INCLUDED)
11 #define KRATOS_NON_LINEAR_HENCKY_PLASTIC_AXISYM_2D_LAW_H_INCLUDED
12 
13 
14 // System includes
15 
16 // External includes
17 
18 // Project includes
20 
21 namespace Kratos
22 {
29 class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) NonLinearHenckyElasticPlasticAxisym2DLaw: public NonLinearHenckyElasticPlastic3DLaw
30 {
31 
32 public:
38  typedef std::size_t SizeType;
39 
40  typedef FlowRule::Pointer FlowRulePointer;
41  typedef YieldCriterion::Pointer YieldCriterionPointer;
42  typedef HardeningLaw::Pointer HardeningLawPointer;
43  typedef Properties::Pointer PropertiesPointer;
44 
50 
59 
60 
62 
67 
68 
73  //HyperElasticPlastic3DLaw& operator=(const HyperElasticPlastic3DLaw& rOther);
74 
79  ConstitutiveLaw::Pointer Clone() const override;
80 
85 
98  {
99  return 2;
100  };
101 
105  SizeType GetStrainSize() const override
106  {
107  return 4;
108  };
109 
110  void GetLawFeatures(Features& rFeatures) override;
111 
112 /* bool Has( const Variable<double>& rThisVariable );
113  bool Has( const Variable<Vector>& rThisVariable );
114  bool Has( const Variable<Matrix>& rThisVariable );
115 
116  double& GetValue( const Variable<double>& rThisVariable, double& rValue );
117  Vector& GetValue( const Variable<Vector>& rThisVariable, Vector& rValue );
118  Matrix& GetValue( const Variable<Matrix>& rThisVariable, Matrix& rValue );
119 
120 
121  void SetValue( const Variable<double>& rVariable,
122  const double& Value,
123  const ProcessInfo& rCurrentProcessInfo );
124  void SetValue( const Variable<Vector>& rThisVariable,
125  const Vector& rValue,
126  const ProcessInfo& rCurrentProcessInfo );
127  void SetValue( const Variable<Matrix>& rThisVariable,
128  const Matrix& rValue,
129  const ProcessInfo& rCurrentProcessInfo );
130 */
139 // void CalculateMaterialResponsePK1 (Parameters & rValues);
140 // Esta función sirve como está definida en la clase madre
141 
148 // void CalculateMaterialResponsePK2 (Parameters & rValues);
149 // Esta función sirve como está definida en la clase madre
150 
157 // void CalculateMaterialResponseKirchhoff (Parameters & rValues);
158 
159 
166 // void CalculateMaterialResponseCauchy (Parameters & rValues);
167 // Esta función sirve como está definida
168 
175 // void FinalizeMaterialResponsePK1 (Parameters & rValues);
176 
183 // void FinalizeMaterialResponsePK2 (Parameters & rValues);
184 
198 // void FinalizeMaterialResponseCauchy (Parameters & rValues);
199 
209 // int Check(const Properties& rProperties, const GeometryType& rGeometry, const ProcessInfo& rCurrentProcessInfo);
210 
211 
212 
219  //String Info() const override;
223  //void PrintInfo(std::ostream& rOStream) const override;
227  //void PrintData(std::ostream& rOStream) const override;
228 
229 protected:
230 
236 
237 // Matrix mElasticLeftCauchyGreen;
238 
239 // FlowRulePointer mpFlowRule;
240 
241 // YieldCriterionPointer mpYieldCriterion;
242 
243 // HardeningLawPointer mpHardeningLaw;
244 
251 
252  Matrix SetConstitutiveMatrixToAppropiateDimension(const Matrix& rElastoPlasticTangentMatrix) override;
253 
254 
258 // Matrix& DeformationGradient3D (Matrix & Matrix2D);
259 
265  void CalculateGreenLagrangeStrain( const Matrix & rRightCauchyGreen,
266  Vector& rStrainVector ) override;
267 
268 
274  void CalculateAlmansiStrain( const Matrix & rLeftCauchyGreen,
275  Vector& rStrainVector ) override;
302 /* virtual void CalculateIsochoricConstitutiveMatrix (const MaterialResponseVariables& rElasticVariables,
303  const Matrix & rInverseDeformationGradientF,
304  const Vector & rIsoStressVector,
305  Matrix& rConstitutiveMatrix);
306 */
307 
312 /* double& IsochoricConstitutiveComponent( double & rCabcd,
313  const MaterialResponseVariables& rElasticVariables,
314  const Matrix & rIsoStressMatrix,
315  const unsigned int& a, const unsigned int& b,
316  const unsigned int& c, const unsigned int& d);
317 */
322 /* double& IsochoricConstitutiveComponent( double & rCabcd,
323  const MaterialResponseVariables& rElasticVariables,
324  const Matrix & rInverseDeformationGradientF,
325  const Matrix & rIsoStressMatrix,
326  const unsigned int& a, const unsigned int& b,
327  const unsigned int& c, const unsigned int& d);
328 
329 */
330 
337 /* virtual void CalculateVolumetricConstitutiveMatrix (const MaterialResponseVariables& rElasticVariables,
338  Matrix& rConstitutiveMatrix);
339 
340 */
350 /* virtual void CalculateVolumetricConstitutiveMatrix (const MaterialResponseVariables& rElasticVariables,
351  const Matrix & rInverseDeformationGradientF,
352  Matrix& rConstitutiveMatrix);
353 */
354 
359 /* double& VolumetricConstitutiveComponent( double & rCabcd,
360  const MaterialResponseVariables& rElasticVariables,
361  const Vector& rFactors,
362  const unsigned int& a, const unsigned int& b,
363  const unsigned int& c, const unsigned int& d);
364 */
369 /* double& VolumetricConstitutiveComponent( double & rCabcd,
370  const MaterialResponseVariables& rElasticVariables,
371  const Matrix & rInverseDeformationGradientF,
372  const Vector & rFactors,
373  const unsigned int& a, const unsigned int& b,
374  const unsigned int& c, const unsigned int& d);
375 
376 */
384 /* virtual void CalculatePlasticConstitutiveMatrix (const MaterialResponseVariables& rElasticVariables,
385  FlowRule::RadialReturnVariables & rReturnMappingVariables,
386  Matrix& rConstitutiveMatrix);
387 */
388 
397 /* virtual void CalculatePlasticConstitutiveMatrix (const MaterialResponseVariables& rElasticVariables,
398  const Matrix & rInverseDeformationGradientF,
399  FlowRule::RadialReturnVariables & rReturnMappingVariables,
400  Matrix& rConstitutiveMatrix);
401 */
402 
407 /* double& PlasticConstitutiveComponent( double & rCabcd,
408  const MaterialResponseVariables& rElasticVariables,
409  const Matrix & rIsoStressMatrix,
410  const FlowRule::PlasticFactors & rScalingFactors,
411  const unsigned int& a, const unsigned int& b,
412  const unsigned int& c, const unsigned int& d);
413 */
418 /* double& PlasticConstitutiveComponent( double & rCabcd,
419  const MaterialResponseVariables& rElasticVariables,
420  const Matrix & rInverseDeformationGradientF,
421  const Matrix & rIsoStressMatrix,
422  const FlowRule::PlasticFactors & rScalingFactors,
423  const unsigned int& a, const unsigned int& b,
424  const unsigned int& c, const unsigned int& d);
425 */
426 
434 /* virtual void CalculatePlasticIsochoricStress( MaterialResponseVariables & rElasticVariables,
435  FlowRule::RadialReturnVariables & rReturnMappingVariables,
436  StressMeasure rStressMeasure,
437  Matrix& rIsoStressMatrix,
438  Vector& rIsoStressVector);
439 */
447 /* virtual void CalculateVolumetricStress( const MaterialResponseVariables & rElasticVariables,
448  Vector& rVolStressVector );
449 */
450 
451 
458 // virtual double& CalculateDomainPressure (const MaterialResponseVariables & rElasticVariables,
459 // double & rPressure);
460 
466 // virtual Vector& CalculateDomainPressureFactors (const MaterialResponseVariables & rElasticVariables,
467 // Vector & rFactors);
468 
469 
476 // virtual bool CheckParameters(Parameters& rValues);
477 //Ll: a pensar què s'ha de fer amb els constructors, els CheckParameters i alguna altra funció virtual
478 private:
479 
482 
483 
487 
488 
492 
493 
497 
498 
503 
504 
508  friend class Serializer;
509 
510  void save(Serializer& rSerializer) const override
511  {
513  }
514 
515  void load(Serializer& rSerializer) override
516  {
518  }
519 
520 
521 
522 
523 }; // Class HenckyElasticPlastic3DLaw
524 
525 } //namespace Kratos
526 
527 #endif //KRATOS_HENCKY_PLASTIC_3D_LAW_H_INCLUDED
Definition: constitutive_law.h:47
std::size_t SizeType
Definition: constitutive_law.h:82
Definition: non_linear_hencky_plastic_3D_law.hpp:32
HardeningLaw::Pointer HardeningLawPointer
Definition: non_linear_hencky_plastic_3D_law.hpp:44
FlowRule::Pointer FlowRulePointer
Definition: non_linear_hencky_plastic_3D_law.hpp:42
YieldCriterion::Pointer YieldCriterionPointer
Definition: non_linear_hencky_plastic_3D_law.hpp:43
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:30
FlowRule::Pointer FlowRulePointer
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:40
YieldCriterion::Pointer YieldCriterionPointer
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:41
Properties::Pointer PropertiesPointer
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:43
ConstitutiveLaw BaseType
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:37
ProcessInfo ProcessInfoType
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:36
SizeType GetStrainSize() const override
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:105
HardeningLaw::Pointer HardeningLawPointer
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:42
KRATOS_CLASS_POINTER_DEFINITION(NonLinearHenckyElasticPlasticAxisym2DLaw)
std::size_t SizeType
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:38
SizeType WorkingSpaceDimension() override
Definition: non_linear_hencky_plastic_axisym_2D_law.hpp:97
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
#define KRATOS_SERIALIZE_SAVE_BASE_CLASS(Serializer, BaseType)
Definition: define.h:812
#define KRATOS_SERIALIZE_LOAD_BASE_CLASS(Serializer, BaseType)
Definition: define.h:815
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21
def load(f)
Definition: ode_solve.py:307