AuxiliaryLogarithmicStressOperators.hpp

Go to the documentation of this file.

/**
 * @file AuxiliaryLogarithmicStressOperators.hpp
 * @brief Material and stress-work blocks for the independent D/theta/Td fields.
 */
 
#ifndef AUXILIARY_LOGARITHMIC_STRESS_OPERATORS_HPP
#define AUXILIARY_LOGARITHMIC_STRESS_OPERATORS_HPP
 
#include <AuxiliaryLogarithmicStress.hpp>
#include <functional>
 
namespace EshelbianPlasticity {
 
struct DataAtIntegrationPts;
 
/** Derived values for one trial state, all stored as GaussByCoeffs.
 * Field interpolation remains in AuxiliaryLogarithmicStressData. These values
 * are overwritten on every element evaluation, including rejected trials.
 */
struct AuxiliaryLogarithmicStressMaterialData {
  using MatrixPtr = boost::shared_ptr<MatrixDouble>;
  enum EnergyComponent { DEVIATORIC, VOLUMETRIC, CONJUGATE };
  enum VolumeComponent { JACOBIAN, FIRST_DERIVATIVE, SECOND_DERIVATIVE };
 
  MatrixPtr materialDeviator = boost::make_shared<MatrixDouble>(); ///< 5
  MatrixPtr elasticStress =
      boost::make_shared<MatrixDouble>(); ///< 5, Td-alphaU*dotD
  MatrixPtr compliance = boost::make_shared<MatrixDouble>(); ///< 5 x 5
  MatrixPtr energy = boost::make_shared<MatrixDouble>();     ///< 3
  MatrixPtr volume = boost::make_shared<MatrixDouble>();     ///< 3
 
  MatrixPtr residualDeviator = boost::make_shared<MatrixDouble>(); ///< 5
  MatrixPtr residualJacobian = boost::make_shared<MatrixDouble>(); ///< 1
  MatrixPtr residualStress = boost::make_shared<MatrixDouble>();   ///< 5
 
  MatrixPtr tangentDeviator = boost::make_shared<MatrixDouble>(); ///< 5 x 5
  MatrixPtr tangentDeviatorJacobian =
      boost::make_shared<MatrixDouble>(); ///< 5 x 1
  MatrixPtr tangentDeviatorStress =
      boost::make_shared<MatrixDouble>(); ///< 5 x 5
  MatrixPtr tangentStressDeviator =
      boost::make_shared<MatrixDouble>();                          ///< 5 x 5
  MatrixPtr tangentJacobian = boost::make_shared<MatrixDouble>();  ///< 1 x 1
  MatrixPtr tangentStress = boost::make_shared<MatrixDouble>();    ///< 5 x 5
  MatrixPtr deviatorRotation = boost::make_shared<MatrixDouble>(); ///< 5 x 3
  MatrixPtr jacobianRotation = boost::make_shared<MatrixDouble>(); ///< 1 x 3
  MatrixPtr deformationDeviator =
      boost::make_shared<MatrixDouble>(); ///< 3 x 3 x 5, dF_ij/dD_A
  MatrixPtr deformationJacobian =
      boost::make_shared<MatrixDouble>(); ///< 3 x 3 x 1, dF_ij/dtheta
};
 
/** Evaluate Te:D-f*(Te)+g(J), Te=Td-alphaU*dotD, without another inversion.
 * This is the mixed stored-energy contribution, not the full stationary
 * mechanical value, which also contains constraint, load and boundary work.
 */
MoFEMErrorCode evaluateAuxiliaryMixedStoredEnergy(
    const AuxiliaryLogarithmicStressData &fields,
    const AuxiliaryLogarithmicStressMaterialData &material,
    VectorDouble &mixed_energy);
 
/// Return parameters already checked by the material's setup validation.
using AuxiliaryLogarithmicMaterialParameters =
    std::function<NeoHookeanLogarithmicMaterial::Parameters(EntityHandle)>;
 
/** Evaluate material diagnostics without assembling any equation blocks.
 * Current field reconstruction and geometry must precede this producer.
 */
MoFEMErrorCode pushAuxiliaryLogarithmicMaterialEvaluation(
    const EshelbianCore &ep,
    boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pipeline,
    boost::shared_ptr<DataAtIntegrationPts> data,
    boost::shared_ptr<AuxiliaryLogarithmicStressMaterialData> material_data,
    AuxiliaryLogarithmicMaterialParameters parameters, bool lhs = false);
 
/** Assemble material data prepared by
 * pushAuxiliaryLogarithmicMaterialEvaluation. The lhs group supplies both
 * partners of each auxiliary/mechanical off-diagonal block; mechanical-only
 * operators remain in the caller's existing pipeline.
 */
MoFEMErrorCode pushAuxiliaryLogarithmicMaterialOps(
    const EshelbianCore &ep,
    boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pipeline,
    boost::shared_ptr<AuxiliaryLogarithmicStressMaterialData> material_data,
    bool lhs);
 
} // namespace EshelbianPlasticity
 
#endif // AUXILIARY_LOGARITHMIC_STRESS_OPERATORS_HPP