ConstitutiveSetting.cpp

Go to the documentation of this file.

 
#include "ConstitutiveSetting.h"
 
#include "Invoke.h"
 
namespace ProcessLib::LargeDeformation
{
namespace ConstitutiveRelations
{
template <int DisplacementDim>
void ConstitutiveSetting<DisplacementDim>::eval(
    ConstitutiveModels<DisplacementDim>& models, double const t,
    double const dt, ParameterLib::SpatialPosition const& x_position,
    MaterialPropertyLib::Medium const& medium, double const T_ref,
    DeformationGradientData<DisplacementDim> const& deformation_gradient_data,
    GradientVectorType const& deformation_gradient_prev,
    StatefulData<DisplacementDim>& state,
    StatefulDataPrev<DisplacementDim> const& prev_state,
    MaterialStateData<DisplacementDim>& mat_state,
    ConstitutiveTempData<DisplacementDim>& tmp,
    ConstitutiveData<DisplacementDim>& cd) const
{
    namespace MPL = MaterialPropertyLib;
 
    auto& deformation_gradient_data_prev = tmp.deformation_gradient_data_prev;
    deformation_gradient_data_prev->deformation_gradient =
        deformation_gradient_prev;
    auto& rho_SR = tmp.rho_SR;
 
    auto& s_mech_data = cd.s_mech_data;
    auto& volumetric_body_force = cd.volumetric_body_force;
 
    Temperature const T{T_ref};
    SpaceTimeData const x_t{x_position, t, dt};
    MediaData const media_data{medium};
 
    assertEvalArgsUnique(models.s_mech_model);
    models.s_mech_model.eval(
        x_t, T, deformation_gradient_data, deformation_gradient_data_prev,
        mat_state, prev_state.stress_data, state.stress_data, s_mech_data);
 
    assertEvalArgsUnique(models.rho_S_model);
    models.rho_S_model.eval(x_t, media_data, T, rho_SR);
 
    assertEvalArgsUnique(models.gravity_model);
    models.gravity_model.eval(rho_SR, volumetric_body_force);
}
 
template struct ConstitutiveSetting<2>;
template struct ConstitutiveSetting<3>;
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::LargeDeformation