thermal_elastic_isotropic_3d.h

Go to the documentation of this file.

// KRATOS ___                _   _ _         _   _             __                       _
//       / __\___  _ __  ___| |_(_) |_ _   _| |_(_)_   _____  / /  __ ___      _____   /_\  _ __  _ __
//      / /  / _ \| '_ \/ __| __| | __| | | | __| \ \ / / _ \/ /  / _` \ \ /\ / / __| //_\\| '_ \| '_  |
//     / /__| (_) | | | \__ \ |_| | |_| |_| | |_| |\ V /  __/ /__| (_| |\ V  V /\__ \/  _  \ |_) | |_) |
//     \____/\___/|_| |_|___/\__|_|\__|\__,_|\__|_| \_/ \___\____/\__,_| \_/\_/ |___/\_/ \_/ .__/| .__/
//                                                                                         |_|   |_|
//
//  License:         BSD License
//                   license: structural_mechanics_application/license.txt
//
//  Main authors:    Alejandro Cornejo
//
 
# pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "custom_constitutive/elastic_isotropic_3d.h"
 
namespace Kratos
{
 
 
 
 
 
class KRATOS_API(CONSTITUTIVE_LAWS_APPLICATION) ThermalElasticIsotropic3D
    : public ElasticIsotropic3D
{
public:
 
 
    using BaseType = ElasticIsotropic3D;
 
    KRATOS_CLASS_POINTER_DEFINITION(ThermalElasticIsotropic3D);
 
 
    ThermalElasticIsotropic3D() 
    {}
 
    ConstitutiveLaw::Pointer Clone() const override
    {
        return Kratos::make_shared<ThermalElasticIsotropic3D>(*this);
    }
 
 
    ~ThermalElasticIsotropic3D() override {}
 
    ThermalElasticIsotropic3D(const ThermalElasticIsotropic3D &rOther)
        : BaseType(rOther),
        mReferenceTemperature(rOther.mReferenceTemperature)
    {
    }
 
 
 
    void CalculateMaterialResponsePK2(ConstitutiveLaw::Parameters &rValues) override;
 
    int Check(
        const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
    double& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<double>& rThisVariable,
        double& rValue
        ) override;
 
    void InitializeMaterial(
        const Properties &rMaterialProperties,
        const GeometryType &rElementGeometry,
        const Vector &rShapeFunctionsValues) override;
 
 
    void CalculatePK2Stress(
        const ConstitutiveLaw::StrainVectorType &rStrainVector,
        ConstitutiveLaw::StressVectorType &rStressVector,
        ConstitutiveLaw::Parameters &rValues) override;
 
    void CalculateElasticMatrix(
        ConstitutiveLaw::VoigtSizeMatrixType& rConstitutiveMatrix,
        ConstitutiveLaw::Parameters& rValues
        ) override;
 
    virtual void SubstractThermalStrain(
        ConstitutiveLaw::StrainVectorType &rStrainVector,
        const double ReferenceTemperature,
        ConstitutiveLaw::Parameters &rParameters,
        const bool IsPlaneStrain = false);
 
 
 
 
 
protected:
 
 
 
 
 
 
private:
 
 
 
    double mReferenceTemperature = 0.0;
 
 
 
 
    double& GetReferenceTemperature()
    {
        return mReferenceTemperature;
    }
 
    void SetReferenceTemperature(const double ToRefTemperature)
    {
        mReferenceTemperature = ToRefTemperature;
    }
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, ElasticIsotropic3D)
        rSerializer.save("ReferenceTemperature", mReferenceTemperature);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, ElasticIsotropic3D)
        rSerializer.load("ReferenceTemperature", mReferenceTemperature);
    }
 
 
}; // Class LinearPlaneStrain
}  // namespace Kratos.