d0/d60/mooney__rivlin__model_8hpp_source.html

 //

 //   Project Name:        KratosConstitutiveModelsApplication $

 //   Created by:          $Author:                JMCarbonell $

 //   Last modified by:    $Co-Author:                         $

 //   Date:                $Date:                   April 2017 $

 //   Revision:            $Revision:                      0.0 $

 //

 //


 #if !defined(KRATOS_MOONEY_RIVLIN_MODEL_H_INCLUDED )

 #define  KRATOS_MOONEY_RIVLIN_MODEL_H_INCLUDED


 // System includes

 #include <string>

 #include <iostream>


 // External includes


 // Project includes

 #include "custom_models/elasticity_models/hyper_elastic_model.hpp"


 namespace Kratos

 {


   class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) MooneyRivlinModel : public HyperElasticModel

   {

   public:


     KRATOS_CLASS_POINTER_DEFINITION( MooneyRivlinModel );


     MooneyRivlinModel();


     MooneyRivlinModel(MooneyRivlinModel const& rOther);


     MooneyRivlinModel& operator=(MooneyRivlinModel const& rOther);


     ConstitutiveModel::Pointer Clone() const override;


     ~MooneyRivlinModel() override;


     int Check(const Properties& rProperties, const ProcessInfo& rCurrentProcessInfo) override;


     std::string Info() const override

     {

         std::stringstream buffer;

         buffer << "MooneyRivlinModel";

         return buffer.str();

     }


     void PrintInfo(std::ostream& rOStream) const override

     {

         rOStream << "MooneyRivlinModel";

     }


     void PrintData(std::ostream& rOStream) const override

     {

       rOStream << "MooneyRivlinModel Data";

     }


   protected:


     void CalculateAndAddStressTensor(HyperElasticDataType& rVariables, MatrixType& rStressMatrix) override;


     double& AddConstitutiveComponent(HyperElasticDataType& rVariables, double &rCabcd,

                          const unsigned int& a, const unsigned int& b,

                          const unsigned int& c, const unsigned int& d) override;


     //************// Strain Data


     void CalculateStrainData(ModelDataType& rValues, HyperElasticDataType& rVariables) override;


     void CalculateScalingFactors(HyperElasticDataType& rVariables) override;


     //************// dW


     double& GetVolumetricFunction1stJDerivative(HyperElasticDataType& rVariables, double& rDerivative) override; //dU/dJ


     double& GetVolumetricFunction2ndJDerivative(HyperElasticDataType& rVariables, double& rDerivative) override; //ddU/dJdJ


     virtual double& GetFunction1stI1Derivative(HyperElasticDataType& rVariables, double& rDerivative); //dW/dI1


     virtual double& GetFunction1stI2Derivative(HyperElasticDataType& rVariables, double& rDerivative); //dW/dI2


     virtual double& GetFunction1stI3Derivative(HyperElasticDataType& rVariables, double& rDerivative); //dW/dI3


     virtual double& GetFunction2ndI1Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI1dI1


     virtual double& GetFunction2ndI2Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI2dI2


     virtual double& GetFunction2ndI3Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI3dI3


     // the implementation of the crossed derivatives have to be added for a more general form (usually they are zero)

     // virtual double& GetFunction2ndI2I1Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI2dI1

     // virtual double& GetFunction2ndI3I1Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI3dI1

     // virtual double& GetFunction2ndI1I2Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI1dI2

     // virtual double& GetFunction2ndI3I2Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI3dI2

     // virtual double& GetFunction2ndI1I3Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI1dI3

     // virtual double& GetFunction2ndI2I3Derivative(HyperElasticDataType& rVariables, double& rDerivative); //ddW/dI2dI3


     //isochoric volumetric slit


     MatrixType& GetIsochoricRightCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative); //dC'/dC


     double& GetIsochoricRightCauchyGreenDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b,

                            const double& c,

                            const double& d); //dC'/dC


     MatrixType& GetIsochoricLeftCauchyGreenDerivative(const StrainData& rStrain,

                               MatrixType& rDerivative); //db'/db


     double& GetIsochoricLeftCauchyGreenDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b,

                           const double& c,

                           const double& d); //db'/db


     //************// right cauchy green: C


     MatrixType& GetI1RightCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative); //dI1/dC


     MatrixType& GetI2RightCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative); //dI2/dC


     MatrixType& GetI3RightCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative); //dI3/dC


     double& GetFourthOrderUnitTensor(double& rValue,

                      const double& a,

                      const double& b,

                      const double& c,

                      const double& d); //ddC/dCdC or ddb/dbdb


     double& GetInverseRightCauchyGreenDerivative(const StrainData& rStrain,

                          double& rDerivative,

                          const double& a,

                          const double& b,

                          const double& c,

                          const double& d); //dC^-1/dC


     //Invariants 1st derivatives by components

     double& GetI1RightCauchyGreen1stDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b); //dI1/dC


     double& GetI2RightCauchyGreen1stDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b); //dI2/dC


     double& GetI3RightCauchyGreen1stDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b); //dI3/dC


     //Invariants Square of the 1st derivatives by components

     double& GetI1RightCauchyGreenSquare1stDerivative(const StrainData& rStrain,

                              double& rDerivative,

                              const double& a,

                              const double& b,

                              const double& c,

                              const double& d); //dI1/dC * dI2/dC


     double& GetI2RightCauchyGreenSquare1stDerivative(const StrainData& rStrain,

                              double& rDerivative,

                              const double& a,

                              const double& b,

                              const double& c,

                              const double& d); //dI2/dC * dI3/dC


     double& GetI3RightCauchyGreenSquare1stDerivative(const StrainData& rStrain,

                              double& rDerivative,

                              const double& a,

                              const double& b,

                              const double& c,

                              const double& d); //dI3/dC * dI3/dC


     //Invariants 2nd derivatives by components

     double& GetI1RightCauchyGreen2ndDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b,

                            const double& c,

                            const double& d); //ddI1/dCdC


     double& GetI2RightCauchyGreen2ndDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b,

                            const double& c,

                            const double& d); //ddI2/dCdC


     double& GetI3RightCauchyGreen2ndDerivative(const StrainData& rStrain,

                            double& rDerivative,

                            const double& a,

                            const double& b,

                            const double& c,

                            const double& d); //ddI3/dCdC


     //************// left cauchy green : b


     MatrixType& GetI1LeftCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative); //dI1/db


     MatrixType& GetI2LeftCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative);  //dI2/db


     MatrixType& GetI3LeftCauchyGreenDerivative(const StrainData& rStrain, MatrixType& rDerivative); //dI3/db


     //Invariants 1st derivatives by components

     double& GetI1LeftCauchyGreen1stDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b); //dI1/db


     double& GetI2LeftCauchyGreen1stDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b); //dI2/db


     double& GetI3LeftCauchyGreen1stDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b); //dI3/db


     //Invariants Square of the 1st derivatives by components

     double& GetI1LeftCauchyGreenSquare1stDerivative(const StrainData& rStrain,

                             double& rDerivative,

                             const double& a,

                             const double& b,

                             const double& c,

                             const double& d); //dI1/db * dI1/db


     double& GetI2LeftCauchyGreenSquare1stDerivative(const StrainData& rStrain,

                             double& rDerivative,

                             const double& a,

                             const double& b,

                             const double& c,

                             const double& d); //dI2/db * dI2/db


     double& GetI3LeftCauchyGreenSquare1stDerivative(const StrainData& rStrain,

                             double& rDerivative,

                             const double& a,

                             const double& b,

                             const double& c,

                             const double& d); //dI3/db * dI3/db


     //Invariants 2nd derivatives by components

     double& GetI1LeftCauchyGreen2ndDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b,

                           const double& c,

                           const double& d); //ddI1/dbdb


     double& GetI2LeftCauchyGreen2ndDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b,

                           const double& c,

                           const double& d); //ddI2/dbdb


     double& GetI3LeftCauchyGreen2ndDerivative(const StrainData& rStrain,

                           double& rDerivative,

                           const double& a,

                           const double& b,

                           const double& c,

                           const double& d); //ddI3/dbdb


   private:


     friend class Serializer;


     void save(Serializer& rSerializer) const override

     {

       KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, HyperElasticModel )

     }


     void load(Serializer& rSerializer) override

     {

       KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, HyperElasticModel )

     }


   }; // Class MooneyRivlinModel


 }  // namespace Kratos.


 #endif // KRATOS_MOONEY_RIVLIN_MODEL_H_INCLUDED  defined

operator=
PeriodicInterfaceProcess & operator=(const PeriodicInterfaceProcess &)=delete

Kratos::HyperElasticModel
Short class definition.
Definition: hyper_elastic_model.hpp:50

Kratos::MooneyRivlinModel
Short class definition.
Definition: mooney_rivlin_model.hpp:50

Kratos::MooneyRivlinModel::PrintData
void PrintData(std::ostream &rOStream) const override
Print object's data.
Definition: mooney_rivlin_model.hpp:124

Kratos::MooneyRivlinModel::Info
std::string Info() const override
Turn back information as a string.
Definition: mooney_rivlin_model.hpp:110

Kratos::MooneyRivlinModel::KRATOS_CLASS_POINTER_DEFINITION
KRATOS_CLASS_POINTER_DEFINITION(MooneyRivlinModel)
Pointer definition of MooneyRivlinModel.

Kratos::MooneyRivlinModel::PrintInfo
void PrintInfo(std::ostream &rOStream) const override
Print information about this object.
Definition: mooney_rivlin_model.hpp:118

Kratos::ProcessInfo
ProcessInfo holds the current value of different solution parameters.
Definition: process_info.h:59

Kratos::Properties
Properties encapsulates data shared by different Elements or Conditions. It can store any type of dat...
Definition: properties.h:69

Kratos::Serializer
The serialization consists in storing the state of an object into a storage format like data file or ...
Definition: serializer.h:123

KRATOS_SERIALIZE_SAVE_BASE_CLASS
#define KRATOS_SERIALIZE_SAVE_BASE_CLASS(Serializer, BaseType)
Definition: define.h:812

KRATOS_SERIALIZE_LOAD_BASE_CLASS
#define KRATOS_SERIALIZE_LOAD_BASE_CLASS(Serializer, BaseType)
Definition: define.h:815

hyper_elastic_model.hpp

Kratos::GeometricalTransformationUtilities::MatrixType
Matrix MatrixType
Definition: geometrical_transformation_utilities.h:55

Kratos
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21

generate_stokes_twofluid_element.a
a
Definition: generate_stokes_twofluid_element.py:77

generate_total_lagrangian_mixed_volumetric_strain_element.b
b
Definition: generate_total_lagrangian_mixed_volumetric_strain_element.py:31

generate_weakly_compressible_navier_stokes_element.c
c
Definition: generate_weakly_compressible_navier_stokes_element.py:108

ode_solve.load
def load(f)
Definition: ode_solve.py:307

ode_solve.d
int d
Definition: ode_solve.py:397