ConstitutiveStaticData: Some design choices

Setup

• each constitutive law (NuTo::ConstitutiveBase) is coupled with its corresponding NuTo::IPConstitutiveLawBase
• there is one instance of NuTo::IPConstitutiveLawBase at each integration point
• it is created via the constitutive law with the method

virtual std::unique_ptr<NuTo::Constitutive::IPConstitutiveLawBase> 
   NuTo::ConstitutiveBase::CreateIPLaw() = 0;

• CASE A: constitutive laws with static data:
  • they return a NuTo::Constitutive::IPConstitutiveLaw\<SpecificLaw> in CreateIPLaw()
  • they define the public typedef StaticDataType

typedef <specific static data type> StaticDataType
// e.g. NuTo::MoistureTransport:
typedef NuTo::Constitutive::StaticData::DataMoistureTransport StaticDatatType
// e.g. NuTo::GradientDamageEngineeringStress:
typedef double StaticDataType

• CASE B: constitutive laws without:
  • they return a NuTo::Constitutive::IPConstitutiveLawWithoutData\<SpecificLaw> in CreateIPLaw()
  • they do not need the typedef

Benefits:

• the IPConstitutiveLawBase automatically allocates, deallocates and provides the correct static data, if the constitutive law needs them.

Evaluate()

The element calls IPConstitutiveLawBase::Evaluate(inout, output). This calls is forwarded to the constitutive law providing an additional vector of StaticDataType, if needed.

Access

via

template <typename TLaw>
DataContainer<TLaw> NuTo::Contitutive::IPConstitutiveLawBase::GetData();

The class DataContainer<TLaw> contains accessors for the data.

E.g:

ElementBase* e;
unsigned ipIndex = 2;
e->GetIPData().GetIPConstitutiveLaw(ipIndex).GetData<NuTo::PhaseField>.SetData(42.);

Serialization

see NuToSerializeStream: How to use it.

Make sure that your StaticDataType is

• serializable by default, that is
  • Eigen::Matrix/Vector
  • NuTo::FullMatrix/FullVector
  • double, bool, int
• provides a NuToSerializeSave(...) and NuToSerializeLoad(...)