AuxiliaryLogarithmicStress.hpp

Go to the documentation of this file.

/**
 * @file AuxiliaryLogarithmicStress.hpp
 * @brief Evaluation of independent logarithmic material fields.
 */
 
#ifndef AUXILIARY_LOGARITHMIC_STRESS_HPP
#define AUXILIARY_LOGARITHMIC_STRESS_HPP
 
#include <NeoHookeanLogarithmicMaterial.hpp>
 
namespace EshelbianPlasticity {
 
struct EshelbianCore;
 
/** One trial state's 5/1/5 coordinates and evaluated symmetric tensors,
 * all in GaussByCoeffs storage. */
struct AuxiliaryLogarithmicStressData {
  using MatrixPtr = boost::shared_ptr<MatrixDouble>;
  MatrixPtr logDeviator = boost::make_shared<MatrixDouble>();
  MatrixPtr logDeviatorDot = boost::make_shared<MatrixDouble>();
  MatrixPtr logJacobian = boost::make_shared<MatrixDouble>();
  MatrixPtr stress = boost::make_shared<MatrixDouble>();
  MatrixPtr logDeviatorTensor = boost::make_shared<MatrixDouble>();
  MatrixPtr stressTensor = boost::make_shared<MatrixDouble>();
};
 
/** Evaluate the selected state and reconstruct H=D+theta*I/3.
 * A null vector reads mesh data; explicit vectors follow the existing field
 * evaluators' ghost-vector convention. Each state owns separate output data.
 */
MoFEMErrorCode pushAuxiliaryLogarithmicFields(
    const EshelbianCore &ep,
    boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pipeline,
    boost::shared_ptr<AuxiliaryLogarithmicStressData> values,
    boost::shared_ptr<MatrixDouble> reconstructed_log_stretch,
    SmartPetscObj<Vec> state = nullptr);
 
} // namespace EshelbianPlasticity
 
#endif // AUXILIARY_LOGARITHMIC_STRESS_HPP