python_outfitted_constitutive_law.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosConstitutiveModelsApplication $
//   Created by:          $Author:                JMCarbonell $
//   Last modified by:    $Co-Author:                         $
//   Date:                $Date:                   April 2017 $
//   Revision:            $Revision:                      0.0 $
//
//
 
#if !defined (KRATOS_PYTHON_OUTFITTED_CONSTITUTIVE_LAW_H_INCLUDED)
#define  KRATOS_PYTHON_OUTFITTED_CONSTITUTIVE_LAW_H_INCLUDED
 
// System includes
#include <iostream>
 
// External includes
#include <cmath>
#include <boost/python.h>
 
// Project includes
#include "includes/constitutive_law.h"
 
#include "constitutive_models_application_variables.h"
 
namespace Kratos
{
class KRATOS_API(CONSTITUTIVE_MODELS_APPLICATION) PythonOutfittedConstitutiveLaw : public ConstitutiveLaw
{
public:
    typedef ProcessInfo      ProcessInfoType;
    typedef ConstitutiveLaw         BaseType;
    typedef std::size_t             SizeType;
 
    KRATOS_CLASS_POINTER_DEFINITION( PythonOutfittedConstitutiveLaw );
 
    PythonOutfittedConstitutiveLaw();
 
    PythonOutfittedConstitutiveLaw(PyObject* pPyConstitutiveLaw);
 
 
    ConstitutiveLaw::Pointer Clone() const;
 
    PythonOutfittedConstitutiveLaw (const PythonOutfittedConstitutiveLaw& rOther);
 
 
    //PythonOutfittedConstitutiveLaw& operator=(const PythonOutfittedConstitutiveLaw& rOther);
 
 
    virtual ~PythonOutfittedConstitutiveLaw();
 
    SizeType WorkingSpaceDimension()
    {
      return boost::python::call_method<int>(mpPyConstitutiveLaw->ptr(),"WorkingSpaceDimension");
    };
 
    SizeType GetStrainSize()
    {
      return boost::python::call_method<int>(mpPyConstitutiveLaw->ptr(),"GetStrainSize");
    };
 
 
    bool Has( const Variable<double>& rThisVariable );
    bool Has( const Variable<Vector>& rThisVariable );
    bool Has( const Variable<Matrix>& rThisVariable );
 
    double& GetValue( const Variable<double>& rThisVariable, double& rValue );
    Vector& GetValue( const Variable<Vector>& rThisVariable, Vector& rValue );
    Matrix& GetValue( const Variable<Matrix>& rThisVariable, Matrix& rValue );
 
 
    void SetValue( const Variable<double>& rVariable,
                   const double& rValue,
                   const ProcessInfo& rCurrentProcessInfo );
    void SetValue( const Variable<Vector>& rThisVariable,
                   const Vector& rValue,
                   const ProcessInfo& rCurrentProcessInfo );
    void SetValue( const Variable<Matrix>& rThisVariable,
                   const Matrix& rValue,
                   const ProcessInfo& rCurrentProcessInfo );
    void InitializeMaterial( const Properties& rProperties,
                             const GeometryType& rElementGeometry,
                             const Vector& rShapeFunctionsValues );
 
    virtual void CalculateMaterialResponsePK1 (Parameters & rValues);
 
    virtual void CalculateMaterialResponsePK2 (Parameters & rValues);
 
    virtual void CalculateMaterialResponseKirchhoff (Parameters & rValues);
 
 
    virtual void CalculateMaterialResponseCauchy (Parameters & rValues);
 
 
    virtual void FinalizeMaterialResponsePK1 (Parameters & rValues);
 
    virtual void FinalizeMaterialResponsePK2 (Parameters & rValues);
 
    virtual void FinalizeMaterialResponseKirchhoff (Parameters & rValues);
 
    virtual void FinalizeMaterialResponseCauchy (Parameters & rValues);
 
 
    void GetLawFeatures(Features& rFeatures);
 
    int Check(const Properties& rProperties, const GeometryType& rElementGeometry, const ProcessInfo& rCurrentProcessInfo);
 
    //virtual String Info() const;
    //virtual void PrintInfo(std::ostream& rOStream) const;
    //virtual void PrintData(std::ostream& rOStream) const;
 
protected:
 
 
    Matrix mInverseTotalDeformationMatrix;  //F0
 
    double mTotalDeformationDet;  //detF0
 
    double mStrainEnergy;
 
 
    virtual void UpdateInternalVariables (Parameters & rValues);
 
 
private:
 
 
    //PyObject* mpPyConstitutiveLaw;
    boost::shared_ptr<boost::python::object> mpPyConstitutiveLaw;
 
 
 
 
    Matrix& Transform2DTo3D (Matrix& rMatrix);
 
 
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, ConstitutiveLaw )
    rSerializer.save("mInverseTotalDeformationMatrix",mInverseTotalDeformationMatrix);
    rSerializer.save("mTotalDeformationDet",mTotalDeformationDet);
    rSerializer.save("mStrainEnergy",mStrainEnergy);
    }
 
    virtual void load(Serializer& rSerializer)
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, ConstitutiveLaw )
    rSerializer.load("mInverseTotalDeformationMatrix",mInverseTotalDeformationMatrix);
    rSerializer.load("mTotalDeformationDet",mTotalDeformationDet);
    rSerializer.load("mStrainEnergy",mStrainEnergy);
    }
 
 
 
}; // Class PythonOutfittedConstitutiveLaw
}  // namespace Kratos.
#endif // KRATOS_PYTHON_OUTFITTED_CONSTITUTIVE_LAW_H_INCLUDED  defined