From bec769943a9fae9103eb28f52d413d691784849b Mon Sep 17 00:00:00 2001
From: Vahid Ziaei Rad <vahidzrad@gmail.com>
Date: Wed, 26 Oct 2022 18:23:33 +0200
Subject: [PATCH] new outputting variables for splitted PhF models
---
.../PhaseField/LocalAssemblerInterface.h | 18 ++++++++++++++++
ProcessLib/PhaseField/PhaseFieldProcess.cpp | 21 +++++++++++++++++++
2 files changed, 39 insertions(+)
diff --git a/ProcessLib/PhaseField/LocalAssemblerInterface.h b/ProcessLib/PhaseField/LocalAssemblerInterface.h
index dde11de5ead..28915ce8aa6 100644
--- a/ProcessLib/PhaseField/LocalAssemblerInterface.h
+++ b/ProcessLib/PhaseField/LocalAssemblerInterface.h
@@ -29,12 +29,30 @@ struct PhaseFieldLocalAssemblerInterface
std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
std::vector<double>& cache) const = 0;
+ virtual std::vector<double> const& getIntPtSigmaTensile(
+ const double t,
+ std::vector<GlobalVector*> const& x,
+ std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
+ std::vector<double>& cache) const = 0;
+
+ virtual std::vector<double> const& getIntPtSigmaCompressive(
+ const double t,
+ std::vector<GlobalVector*> const& x,
+ std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
+ std::vector<double>& cache) const = 0;
+
virtual std::vector<double> const& getIntPtEpsilon(
const double t,
std::vector<GlobalVector*> const& x,
std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
std::vector<double>& cache) const = 0;
+ virtual std::vector<double> const& getIntPtEpsilonTensile(
+ const double t,
+ std::vector<GlobalVector*> const& x,
+ std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
+ std::vector<double>& cache) const = 0;
+
virtual void computeCrackIntegral(
std::size_t mesh_item_id,
std::vector<
diff --git a/ProcessLib/PhaseField/PhaseFieldProcess.cpp b/ProcessLib/PhaseField/PhaseFieldProcess.cpp
index 1408afc5f14..000f99472b6 100644
--- a/ProcessLib/PhaseField/PhaseFieldProcess.cpp
+++ b/ProcessLib/PhaseField/PhaseFieldProcess.cpp
@@ -137,6 +137,27 @@ void PhaseFieldProcess<DisplacementDim>::initializeConcreteProcess(
getExtrapolator(), _local_assemblers,
&LocalAssemblerInterface::getIntPtEpsilon));
+ _secondary_variables.addSecondaryVariable(
+ "sigma_tensile",
+ makeExtrapolator(MathLib::KelvinVector::KelvinVectorType<
+ DisplacementDim>::RowsAtCompileTime,
+ getExtrapolator(), _local_assemblers,
+ &LocalAssemblerInterface::getIntPtSigmaTensile));
+
+ _secondary_variables.addSecondaryVariable(
+ "sigma_compressive",
+ makeExtrapolator(MathLib::KelvinVector::KelvinVectorType<
+ DisplacementDim>::RowsAtCompileTime,
+ getExtrapolator(), _local_assemblers,
+ &LocalAssemblerInterface::getIntPtSigmaCompressive));
+
+ _secondary_variables.addSecondaryVariable(
+ "eps_tensile",
+ makeExtrapolator(MathLib::KelvinVector::KelvinVectorType<
+ DisplacementDim>::RowsAtCompileTime,
+ getExtrapolator(), _local_assemblers,
+ &LocalAssemblerInterface::getIntPtEpsilonTensile));
+
// Initialize local assemblers after all variables have been set.
GlobalExecutor::executeMemberOnDereferenced(
&LocalAssemblerInterface::initialize, _local_assemblers,
--
GitLab