PlasticOps Namespace Reference

Classes
struct	CommonData
	[Common data] More...
struct	Monitor
struct	OpCalculateConstraintsLhs_dEPImpl
struct	OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculateConstraintsLhs_dTAUImpl
struct	OpCalculateConstraintsLhs_dTAUImpl< GAUSS, AssemblyDomainEleOp >
struct	OpCalculateConstraintsLhs_dUImpl
struct	OpCalculateConstraintsLhs_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculateConstraintsLhs_LogStrain_dUImpl
struct	OpCalculateConstraintsLhs_LogStrain_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculateConstraintsRhsImpl
struct	OpCalculateConstraintsRhsImpl< GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticFlowLhs_dEPImpl
struct	OpCalculatePlasticFlowLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticFlowLhs_dTAUImpl
struct	OpCalculatePlasticFlowLhs_dTAUImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticFlowLhs_dUImpl
struct	OpCalculatePlasticFlowLhs_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticFlowLhs_LogStrain_dUImpl
struct	OpCalculatePlasticFlowLhs_LogStrain_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticFlowRhsImpl
struct	OpCalculatePlasticFlowRhsImpl< DIM, GAUSS, AssemblyDomainEleOp >
	[Calculate stress] More...
struct	OpCalculatePlasticInternalForceLhs_dEPImpl
struct	OpCalculatePlasticInternalForceLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl
struct	OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpCalculatePlasticityImpl
struct	OpCalculatePlasticityImpl< DIM, GAUSS, DomainEleOp >
struct	OpCalculatePlasticSurfaceImpl
struct	OpCalculatePlasticSurfaceImpl< DIM, GAUSS, DomainEleOp >
	[Auxiliary functions functions More...
struct	OpPlasticStressImpl
struct	OpPlasticStressImpl< DIM, GAUSS, DomainEleOp >
struct	PlasticityIntegrators

Typedefs
using	Pip = boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator >
using	CommonPlasticPtr = boost::shared_ptr< PlasticOps::CommonData >
using	CommonHenkyPtr = boost::shared_ptr< HenckyOps::CommonData >

Functions
template<int DIM>
MoFEMErrorCode	addMatBlockOps (MoFEM::Interface &m_field, std::string block_name, Pip &pip, boost::shared_ptr< MatrixDouble > mat_D_Ptr, boost::shared_ptr< CommonData::BlockParams > mat_params_ptr, double scale_value, Sev sev)
template<int DIM, IntegrationType I, typename DomainEleOp >
auto	createCommonPlasticOps (MoFEM::Interface &m_field, std::string block_name, boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string u, std::string ep, std::string tau, double scale, Sev sev)
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >
MoFEMErrorCode	opFactoryDomainRhs (MoFEM::Interface &m_field, std::string block_name, Pip &pip, std::string u, std::string ep, std::string tau)
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >
MoFEMErrorCode	opFactoryDomainLhs (MoFEM::Interface &m_field, std::string block_name, Pip &pip, std::string u, std::string ep, std::string tau)
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >
MoFEMErrorCode	opFactoryDomainReactions (MoFEM::Interface &m_field, std::string block_name, Pip &pip, std::string u, std::string ep, std::string tau)
template<int DIM>
auto	diff_tensor (FTensor::Number< DIM >)
	[Lambda functions] More...
template<int DIM>
auto	symm_L_tensor (FTensor::Number< DIM >)
template<int DIM>
auto	diff_symmetrize (FTensor::Number< DIM >)
template<typename T >
double	trace (FTensor::Tensor2_symmetric< T, 2 > &t_stress)
template<typename T >
double	trace (FTensor::Tensor2_symmetric< T, 3 > &t_stress)
template<typename T , int DIM>
auto	deviator (FTensor::Tensor2_symmetric< T, DIM > &t_stress, double trace, FTensor::Tensor2_symmetric< double, DIM > &t_alpha, FTensor::Number< DIM >)
template<typename T >
auto	deviator (FTensor::Tensor2_symmetric< T, 2 > &t_stress, double trace, FTensor::Tensor2_symmetric< double, 2 > &&t_alpha)
template<typename T >
auto	deviator (FTensor::Tensor2_symmetric< T, 3 > &t_stress, double trace, FTensor::Tensor2_symmetric< double, 3 > &&t_alpha)
template<int DIM>
auto	diff_deviator (FTensor::Ddg< double, DIM, DIM > &&t_diff_stress, FTensor::Number< DIM >)
auto	diff_deviator (FTensor::Ddg< double, 2, 2 > &&t_diff_stress)
auto	diff_deviator (FTensor::Ddg< double, 3, 3 > &&t_diff_stress)
double	platsic_surface (FTensor::Tensor2_symmetric< double, 3 > &&t_stress_deviator)
template<int DIM>
auto	plastic_flow (long double f, FTensor::Tensor2_symmetric< double, 3 > &&t_dev_stress, FTensor::Ddg< double, 3, DIM > &&t_diff_deviator)
template<typename T , int DIM>
auto	diff_plastic_flow_dstress (long double f, FTensor::Tensor2_symmetric< T, DIM > &t_flow, FTensor::Ddg< double, 3, DIM > &&t_diff_deviator)
template<typename T , int DIM>
auto	diff_plastic_flow_dstrain (FTensor::Ddg< T, DIM, DIM > &t_D, FTensor::Ddg< double, DIM, DIM > &&t_diff_plastic_flow_dstress)
double	constrian_sign (double x)
double	constrain_abs (double x)
double	w (double eqiv, double dot_tau, double f, double sigma_y, double sigma_Y)
double	constraint (double eqiv, double dot_tau, double f, double sigma_y, double abs_w, double vis_H, double sigma_Y)
double	diff_constrain_ddot_tau (double sign, double eqiv, double dot_tau, double vis_H, double sigma_Y)
double	diff_constrain_deqiv (double sign, double eqiv, double dot_tau, double sigma_Y)
auto	diff_constrain_df (double sign)
auto	diff_constrain_dsigma_y (double sign)
template<typename T , int DIM>
auto	diff_constrain_dstress (double diff_constrain_df, FTensor::Tensor2_symmetric< T, DIM > &t_plastic_flow)
template<typename T1 , typename T2 , int DIM>
auto	diff_constrain_dstrain (T1 &t_D, T2 &&t_diff_constrain_dstress)
template<typename T , int DIM>
auto	equivalent_strain_dot (FTensor::Tensor2_symmetric< T, DIM > &t_plastic_strain_dot)
template<typename T1 , typename T2 , typename T3 , int DIM>
auto	diff_equivalent_strain_dot (const T1 eqiv, T2 &t_plastic_strain_dot, T3 &t_diff_plastic_strain, FTensor::Number< DIM >)
static auto	get_mat_tensor_sym_dvector (size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
	[Lambda functions] More...
static auto	get_mat_tensor_sym_dvector (size_t rr, MatrixDouble &mat, FTensor::Number< 3 >)
static auto	get_mat_tensor_sym_dtensor_sym (size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
static auto	get_mat_tensor_sym_dtensor_sym (size_t rr, MatrixDouble &mat, FTensor::Number< 3 >)
static auto	get_mat_tensor_sym_dscalar (size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
static auto	get_mat_tensor_sym_dscalar (size_t rr, MatrixDouble &mat, FTensor::Number< 3 >)
auto	get_mat_scalar_dtensor_sym (MatrixDouble &mat, FTensor::Number< 2 >)
auto	get_mat_scalar_dtensor_sym (MatrixDouble &mat, FTensor::Number< 3 >)
static FTensor::Tensor2< FTensor::PackPtr< double *, 3 >, 2, 3 >	get_mat_vector_dtensor_sym (size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
static FTensor::Tensor2< FTensor::PackPtr< double *, 6 >, 3, 6 >	get_mat_vector_dtensor_sym (size_t rr, MatrixDouble &mat, FTensor::Number< 3 >)

Variables
FTensor::Index< 'I', 3 >	I
	[Common data] More...
FTensor::Index< 'J', 3 >	J
FTensor::Index< 'M', 3 >	M
FTensor::Index< 'N', 3 >	N

Typedef Documentation

CommonHenkyPtr

using PlasticOps::CommonHenkyPtr = typedef boost::shared_ptr<HenckyOps::CommonData>

Examples

PlasticOps.hpp.

Definition at line 244 of file PlasticOps.hpp.

CommonPlasticPtr

using PlasticOps::CommonPlasticPtr = typedef boost::shared_ptr<PlasticOps::CommonData>

Examples

PlasticOps.hpp.

Definition at line 243 of file PlasticOps.hpp.

Pip

using PlasticOps::Pip = typedef boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator>

Examples

PlasticOps.hpp.

Definition at line 242 of file PlasticOps.hpp.

Function Documentation

addMatBlockOps()

template<int DIM>

MoFEMErrorCode PlasticOps::addMatBlockOps	(	MoFEM::Interface &	m_field,
		std::string	block_name,
		Pip &	pip,
		boost::shared_ptr< MatrixDouble >	mat_D_Ptr,
		boost::shared_ptr< CommonData::BlockParams >	mat_params_ptr,
		double	scale_value,
		Sev	sev
	)

Extract block data from meshsets

Parameters

m_field
meshset_vec_ptr
sev

Returns

MoFEMErrorCode

Get elasticity tensor

Calculate elasticity tensor for given material parameters

Parameters

mat_D_ptr
bulk_modulus_K
shear_modulus_G

Returns

MoFEMErrorCode

[Calculate elasticity tensor]

[Calculate elasticity tensor]

Examples

PlasticOps.hpp.

Definition at line 248 of file PlasticOps.hpp.

                                            {
  MoFEMFunctionBegin;
 
  struct OpMatBlocks : public DomainEleOp {
    OpMatBlocks(boost::shared_ptr<MatrixDouble> m_D_ptr,
                boost::shared_ptr<CommonData::BlockParams> mat_params_ptr,
                double scale_value, MoFEM::Interface &m_field, Sev sev,
                std::vector<const CubitMeshSets *> meshset_vec_ptr)
        : DomainEleOp(NOSPACE, DomainEleOp::OPSPACE), matDPtr(m_D_ptr),
          matParamsPtr(mat_params_ptr), scaleVal(scale_value) {
      CHK_THROW_MESSAGE(extractBlockData(m_field, meshset_vec_ptr, sev),
                        "Can not get data from block");
    }
 
    MoFEMErrorCode doWork(int side, EntityType type,
                          EntitiesFieldData::EntData &data) {
      MoFEMFunctionBegin;
 
      auto getK = [](auto &p) {
        auto young_modulus = p[CommonData::YOUNG_MODULUS];
        auto poisson_ratio = p[CommonData::POISSON_RATIO];
        return young_modulus / (3 * (1 - 2 * poisson_ratio));
      };
 
      auto getG = [](auto &p) {
        auto young_modulus = p[CommonData::YOUNG_MODULUS];
        auto poisson_ratio = p[CommonData::POISSON_RATIO];
        return young_modulus / (2 * (1 + poisson_ratio));
      };
 
      auto scale_fun = [this](auto &p) {
        p[CommonData::YOUNG_MODULUS] *= scaleVal;
        p[CommonData::SIGMA_Y] *= scaleVal;
        p[CommonData::H] *= scaleVal;
        p[CommonData::VIS_H] *= scaleVal;
        p[CommonData::QINF] *= scaleVal;
        p[CommonData::C1_k] *= scaleVal;
      };
 
      for (auto &b : blockData) {
        if (b.blockEnts.find(getFEEntityHandle()) != b.blockEnts.end()) {
          *matParamsPtr = b.bParams;
          scale_fun(*matParamsPtr);
          CHKERR getMatDPtr(matDPtr, getK(*matParamsPtr), getG(*matParamsPtr));
          MoFEMFunctionReturnHot(0);
        }
      }
 
      (*matParamsPtr) = {young_modulus, poisson_ratio, sigmaY, H,
                         visH,          Qinf,          b_iso,  C1_k};
      scale_fun(*matParamsPtr);
      CHKERR getMatDPtr(matDPtr, getK(*matParamsPtr), getG(*matParamsPtr));
 
      MoFEMFunctionReturn(0);
    }
 
  private:
    boost::shared_ptr<MatrixDouble> matDPtr;
    boost::shared_ptr<CommonData::BlockParams> matParamsPtr;
    const double scaleVal;
 
    struct BlockData {
      std::array<double, CommonData::LAST_PARAM> bParams;
      Range blockEnts;
    };
    std::vector<BlockData> blockData;
 
    /**
     * @brief Extract block data from meshsets
     *
     * @param m_field
     * @param meshset_vec_ptr
     * @param sev
     * @return MoFEMErrorCode
     */
    MoFEMErrorCode
    extractBlockData(MoFEM::Interface &m_field,
                     std::vector<const CubitMeshSets *> meshset_vec_ptr,
                     Sev sev) {
      MoFEMFunctionBegin;
 
      for (auto m : meshset_vec_ptr) {
        MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock") << *m;
        std::vector<double> block_data;
        CHKERR m->getAttributes(block_data);
        if (block_data.size() != CommonData::LAST_PARAM) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "Wrong number of block attribute");
        }
        auto get_block_ents = [&]() {
          Range ents;
          CHKERR m_field.get_moab().get_entities_by_handle(m->meshset, ents,
                                                           true);
          return ents;
        };
 
        CommonData::BlockParams block_params;
        for (auto i = 0; i != CommonData::LAST_PARAM; ++i) {
          block_params[i] = block_data[i];
        }
 
        MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock")
            << "E = " << block_params[CommonData::YOUNG_MODULUS]
            << " nu = " << block_params[CommonData::POISSON_RATIO];
        MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock")
            << std::endl
            << "sigma_y = " << block_params[CommonData::SIGMA_Y] << std::endl
            << "h = " << block_params[CommonData::H] << std::endl
            << "vis_h = " << block_params[CommonData::VIS_H] << std::endl
            << "qinf = " << block_params[CommonData::QINF] << std::endl
            << "biso = " << block_params[CommonData::BISO] << std::endl
            << "C1_k = " << block_params[CommonData::C1_k] << std::endl;
 
        blockData.push_back({block_params, get_block_ents()});
      }
      MOFEM_LOG_CHANNEL("WORLD");
      MoFEMFunctionReturn(0);
    }
 
    /**
     * @brief Get elasticity tensor
     *
     * Calculate elasticity tensor for given material parameters
     *
     * @param mat_D_ptr
     * @param bulk_modulus_K
     * @param shear_modulus_G
     * @return MoFEMErrorCode
     *
     */
    MoFEMErrorCode getMatDPtr(boost::shared_ptr<MatrixDouble> mat_D_ptr,
                              double bulk_modulus_K, double shear_modulus_G) {
      MoFEMFunctionBegin;
      //! [Calculate elasticity tensor]
      auto set_material_stiffness = [&]() {
        FTensor::Index<'i', DIM> i;
        FTensor::Index<'j', DIM> j;
        FTensor::Index<'k', DIM> k;
        FTensor::Index<'l', DIM> l;
        constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
        double A = (DIM == 2)
                       ? 2 * shear_modulus_G /
                             (bulk_modulus_K + (4. / 3.) * shear_modulus_G)
                       : 1;
        auto t_D = getFTensor4DdgFromMat<DIM, DIM, 0>(*mat_D_ptr);
        t_D(i, j, k, l) =
            2 * shear_modulus_G * ((t_kd(i, k) ^ t_kd(j, l)) / 4.) +
            A * (bulk_modulus_K - (2. / 3.) * shear_modulus_G) * t_kd(i, j) *
                t_kd(k, l);
      };
      //! [Calculate elasticity tensor]
      constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
      mat_D_ptr->resize(size_symm * size_symm, 1);
      set_material_stiffness();
      MoFEMFunctionReturn(0);
    }
  };
 
  // push operator to calculate material stiffness matrix for each block
  pip.push_back(new OpMatBlocks(
      mat_D_Ptr, mat_params_ptr, scale_value, m_field, sev,
 
      // Get blockset using regular expression
      m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
          (boost::format("%s(.*)") % block_name).str()
 
              ))
 
          ));
 
  MoFEMFunctionReturn(0);
}

constrain_abs()

double PlasticOps::constrain_abs ( double  x )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 265 of file PlasticOpsGeneric.hpp.

                                      {
  const auto y = -x / zeta;
  if (y > std::numeric_limits<float>::max_exponent10 ||
      y < std::numeric_limits<float>::min_exponent10) {
    return std::abs(x);
  } else {
    const double e = std::exp(y);
    return x + 2 * zeta * std::log1p(e);
  }
};

constraint()

double PlasticOps::constraint ( double  eqiv, double  dot_tau, double  f, double  sigma_y, double  abs_w, double  vis_H, double  sigma_Y  )

inline

\[ \dot{\tau} - \frac{1}{2}\left\{\dot{\tau} + (f(\pmb\sigma) - \sigma_y) + \| \dot{\tau} + (f(\pmb\sigma) - \sigma_y) \|\right\} = 0 \\ c_n \sigma_y \dot{\tau} - \frac{1}{2}\left\{c_n\sigma_y \dot{\tau} + (f(\pmb\sigma) - \sigma_y) + \| c_n \sigma_y \dot{\tau} + (f(\pmb\sigma) - \sigma_y) \|\right\} = 0 \]

Examples

PlasticOpsGeneric.hpp.

Definition at line 292 of file PlasticOpsGeneric.hpp.

                                                                     {
  return vis_H * dot_tau +
         (sigma_Y / 2) * ((cn0 * (dot_tau - eqiv) + cn1 * (eqiv) -
                           (f - sigma_y) / sigma_Y) -
                          abs_w);
};

constrian_sign()

double PlasticOps::constrian_sign ( double  x )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 251 of file PlasticOpsGeneric.hpp.

                                       {
  const auto y = x / zeta;
  if (y > std::numeric_limits<float>::max_exponent10 ||
      y < std::numeric_limits<float>::min_exponent10) {
    if (x > 0)
      return 1.;
    else
      return -1.;
  } else {
    const auto e = std::exp(y);
    return (e - 1) / (1 + e);
  }
};

createCommonPlasticOps()

template<int DIM, IntegrationType I, typename DomainEleOp >

auto PlasticOps::createCommonPlasticOps	(	MoFEM::Interface &	m_field,
		std::string	block_name,
		boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	u,
		std::string	ep,
		std::string	tau,
		double	scale,
		Sev	sev
	)

Examples

PlasticOps.hpp.

Definition at line 426 of file PlasticOps.hpp.

                                                                     {
 
  using P = PlasticityIntegrators<DomainEleOp>;
 
  auto common_plastic_ptr = boost::make_shared<PlasticOps::CommonData>();
  common_plastic_ptr = boost::make_shared<PlasticOps::CommonData>();
 
  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
  auto make_d_mat = []() {
    return boost::make_shared<MatrixDouble>(size_symm * size_symm, 1);
  };
 
  common_plastic_ptr->mDPtr = make_d_mat();
  common_plastic_ptr->mGradPtr = boost::make_shared<MatrixDouble>();
  common_plastic_ptr->mStrainPtr = boost::make_shared<MatrixDouble>();
  common_plastic_ptr->mStressPtr = boost::make_shared<MatrixDouble>();
 
  auto m_D_ptr = common_plastic_ptr->mDPtr;
 
  CHK_THROW_MESSAGE(addMatBlockOps<DIM>(m_field, block_name, pip, m_D_ptr,
                                        common_plastic_ptr->getParamsPtr(),
                                        scale, sev),
                    "add mat block ops");
 
  pip.push_back(new OpCalculateScalarFieldValues(
      tau, common_plastic_ptr->getPlasticTauPtr()));
  pip.push_back(new OpCalculateTensor2SymmetricFieldValues<DIM>(
      ep, common_plastic_ptr->getPlasticStrainPtr()));
  pip.push_back(new OpCalculateVectorFieldGradient<DIM, DIM>(
      u, common_plastic_ptr->mGradPtr));
 
  CommonHenkyPtr common_henky_ptr;
 
  if (is_large_strains) {
    common_henky_ptr = boost::make_shared<HenckyOps::CommonData>();
    common_henky_ptr->matGradPtr = common_plastic_ptr->mGradPtr;
    common_henky_ptr->matDPtr = common_plastic_ptr->mDPtr;
    common_henky_ptr->matLogCPlastic =
        common_plastic_ptr->getPlasticStrainPtr();
    common_plastic_ptr->mStrainPtr = common_henky_ptr->getMatLogC();
    common_plastic_ptr->mStressPtr = common_henky_ptr->getMatHenckyStress();
 
    using H = HenckyOps::HenkyIntegrators<DomainEleOp>;
 
    pip.push_back(new typename H::template OpCalculateEigenVals<DIM, I>(
        u, common_henky_ptr));
    pip.push_back(
        new typename H::template OpCalculateLogC<DIM, I>(u, common_henky_ptr));
    pip.push_back(new typename H::template OpCalculateLogC_dC<DIM, I>(
        u, common_henky_ptr));
    pip.push_back(new
                  typename H::template OpCalculateHenckyPlasticStress<DIM, I>(
                      u, common_henky_ptr, m_D_ptr));
    pip.push_back(new typename H::template OpCalculatePiolaStress<DIM, I>(
        u, common_henky_ptr));
  } else {
 
    pip.push_back(new OpSymmetrizeTensor<SPACE_DIM>(
        u, common_plastic_ptr->mGradPtr, common_plastic_ptr->mStrainPtr));
    pip.push_back(new typename P::template OpPlasticStress<DIM, I>(
        u, common_plastic_ptr, m_D_ptr));
  }
 
  pip.push_back(new typename P::template OpCalculatePlasticSurface<DIM, I>(
      u, common_plastic_ptr));
 
  return std::make_tuple(common_plastic_ptr, common_henky_ptr);
}

deviator() [1/3]

template<typename T >

auto PlasticOps::deviator ( FTensor::Tensor2_symmetric< T, 2 > &  t_stress, double  trace, FTensor::Tensor2_symmetric< double, 2 > &&  t_alpha  )

inline

Definition at line 113 of file PlasticOpsGeneric.hpp.

                                                                    {
  return deviator(t_stress, trace, t_alpha, FTensor::Number<2>());
};

deviator() [2/3]

template<typename T >

auto PlasticOps::deviator ( FTensor::Tensor2_symmetric< T, 3 > &  t_stress, double  trace, FTensor::Tensor2_symmetric< double, 3 > &&  t_alpha  )

inline

Definition at line 119 of file PlasticOpsGeneric.hpp.

                                                                    {
  return deviator(t_stress, trace, t_alpha, FTensor::Number<3>());
};

deviator() [3/3]

template<typename T , int DIM>

auto PlasticOps::deviator ( FTensor::Tensor2_symmetric< T, DIM > &  t_stress, double  trace, FTensor::Tensor2_symmetric< double, DIM > &  t_alpha, FTensor::Number< DIM >    )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 92 of file PlasticOpsGeneric.hpp.

  {
  FTensor::Tensor2_symmetric<double, 3> t_dev;
  t_dev(I, J) = 0;
  for (int ii = 0; ii != DIM; ++ii)
    for (int jj = ii; jj != DIM; ++jj)
      t_dev(ii, jj) = t_stress(ii, jj);
  t_dev(0, 0) -= trace;
  t_dev(1, 1) -= trace;
  t_dev(2, 2) -= trace;
  for (int ii = 0; ii != DIM; ++ii)
    for (int jj = ii; jj != DIM; ++jj)
      t_dev(ii, jj) -= t_alpha(ii, jj);
  return t_dev;
};

diff_constrain_ddot_tau()

double PlasticOps::diff_constrain_ddot_tau ( double  sign, double  eqiv, double  dot_tau, double  vis_H, double  sigma_Y  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 300 of file PlasticOpsGeneric.hpp.

                                                                    {
  return vis_H + (sigma_Y / 2) * (cn0);
};

diff_constrain_deqiv()

double PlasticOps::diff_constrain_deqiv ( double  sign, double  eqiv, double  dot_tau, double  sigma_Y  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 305 of file PlasticOpsGeneric.hpp.

                                                   {
  return (sigma_Y / 2) * (-cn0 + cn1 * (1 - sign));
};

diff_constrain_df()

auto PlasticOps::diff_constrain_df ( double  sign )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 310 of file PlasticOpsGeneric.hpp.

{ return (-1 - sign) / 2; };

diff_constrain_dsigma_y()

auto PlasticOps::diff_constrain_dsigma_y ( double  sign )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 312 of file PlasticOpsGeneric.hpp.

{ return (1 + sign) / 2; }

diff_constrain_dstrain()

template<typename T1 , typename T2 , int DIM>

auto PlasticOps::diff_constrain_dstrain ( T1 &  t_D, T2 &&  t_diff_constrain_dstress  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 326 of file PlasticOpsGeneric.hpp.

                                                                           {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Tensor2_symmetric<double, DIM> t_diff_constrain_dstrain;
  t_diff_constrain_dstrain(k, l) =
      t_diff_constrain_dstress(i, j) * t_D(i, j, k, l);
  return t_diff_constrain_dstrain;
};

diff_constrain_dstress()

template<typename T , int DIM>

auto PlasticOps::diff_constrain_dstress ( double  diff_constrain_df, FTensor::Tensor2_symmetric< T, DIM > &  t_plastic_flow  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 316 of file PlasticOpsGeneric.hpp.

                                                                         {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Tensor2_symmetric<double, DIM> t_diff_constrain_dstress;
  t_diff_constrain_dstress(i, j) = diff_constrain_df * t_plastic_flow(i, j);
  return t_diff_constrain_dstress;
};

diff_deviator() [1/3]

auto PlasticOps::diff_deviator ( FTensor::Ddg< double, 2, 2 > &&  t_diff_stress )

inline

Definition at line 157 of file PlasticOpsGeneric.hpp.

                                                                    {
  return diff_deviator(std::move(t_diff_stress), FTensor::Number<2>());
}

diff_deviator() [2/3]

auto PlasticOps::diff_deviator ( FTensor::Ddg< double, 3, 3 > &&  t_diff_stress )

inline

Definition at line 161 of file PlasticOpsGeneric.hpp.

                                                                    {
  return diff_deviator(std::move(t_diff_stress), FTensor::Number<3>());
}

diff_deviator() [3/3]

template<int DIM>

auto PlasticOps::diff_deviator ( FTensor::Ddg< double, DIM, DIM > &&  t_diff_stress, FTensor::Number< DIM >    )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 125 of file PlasticOpsGeneric.hpp.

                                              {
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Ddg<double, 3, DIM> t_diff_deviator;
  t_diff_deviator(I, J, k, l) = 0;
  for (int ii = 0; ii != DIM; ++ii)
    for (int jj = ii; jj != DIM; ++jj)
      for (int kk = 0; kk != DIM; ++kk)
        for (int ll = kk; ll != DIM; ++ll)
          t_diff_deviator(ii, jj, kk, ll) = t_diff_stress(ii, jj, kk, ll);
 
  constexpr double third = boost::math::constants::third<double>();
 
  t_diff_deviator(0, 0, 0, 0) -= third;
  t_diff_deviator(0, 0, 1, 1) -= third;
 
  t_diff_deviator(1, 1, 0, 0) -= third;
  t_diff_deviator(1, 1, 1, 1) -= third;
 
  t_diff_deviator(2, 2, 0, 0) -= third;
  t_diff_deviator(2, 2, 1, 1) -= third;
 
  if constexpr (DIM == 3) {
    t_diff_deviator(0, 0, 2, 2) -= third;
    t_diff_deviator(1, 1, 2, 2) -= third;
    t_diff_deviator(2, 2, 2, 2) -= third;
  }
 
  return t_diff_deviator;
};

diff_equivalent_strain_dot()

template<typename T1 , typename T2 , typename T3 , int DIM>

auto PlasticOps::diff_equivalent_strain_dot ( const T1  eqiv, T2 &  t_plastic_strain_dot, T3 &  t_diff_plastic_strain, FTensor::Number< DIM >    )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 349 of file PlasticOpsGeneric.hpp.

                                                           {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  constexpr double A = 2. / 3;
  FTensor::Tensor2_symmetric<double, DIM> t_diff_eqiv;
  t_diff_eqiv(i, j) = A * (t_plastic_strain_dot(k, l) / eqiv) *
                      t_diff_plastic_strain(k, l, i, j);
  return t_diff_eqiv;
};

diff_plastic_flow_dstrain()

template<typename T , int DIM>

auto PlasticOps::diff_plastic_flow_dstrain ( FTensor::Ddg< T, DIM, DIM > &  t_D, FTensor::Ddg< double, DIM, DIM > &&  t_diff_plastic_flow_dstress  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 224 of file PlasticOpsGeneric.hpp.

                                                                {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Index<'m', DIM> m;
  FTensor::Index<'n', DIM> n;
  FTensor::Ddg<double, DIM, DIM> t_diff_flow;
  t_diff_flow(i, j, k, l) =
      t_diff_plastic_flow_dstress(i, j, m, n) * t_D(m, n, k, l);
  return t_diff_flow;
};

diff_plastic_flow_dstress()

template<typename T , int DIM>

auto PlasticOps::diff_plastic_flow_dstress ( long double  f, FTensor::Tensor2_symmetric< T, DIM > &  t_flow, FTensor::Ddg< double, 3, DIM > &&  t_diff_deviator  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 208 of file PlasticOpsGeneric.hpp.

                                                                        {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Ddg<double, DIM, DIM> t_diff_flow;
  t_diff_flow(i, j, k, l) =
      (1.5 * (t_diff_deviator(M, N, i, j) * t_diff_deviator(M, N, k, l) -
              (2. / 3.) * t_flow(i, j) * t_flow(k, l))) /
      f;
  return t_diff_flow;
};

diff_symmetrize()

template<int DIM>

auto PlasticOps::diff_symmetrize ( FTensor::Number< DIM >  )

inline

Examples

PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, and PlasticOpsSmallStrains.hpp.

Definition at line 43 of file PlasticOpsGeneric.hpp.

                                                                   {
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
 
  FTensor::Tensor4<double, DIM, DIM, DIM, DIM> t_diff;
 
  t_diff(i, j, k, l) = 0;
  t_diff(0, 0, 0, 0) = 1;
  t_diff(1, 1, 1, 1) = 1;
 
  t_diff(1, 0, 1, 0) = 0.5;
  t_diff(1, 0, 0, 1) = 0.5;
 
  t_diff(0, 1, 0, 1) = 0.5;
  t_diff(0, 1, 1, 0) = 0.5;
 
  if constexpr (DIM == 3) {
    t_diff(2, 2, 2, 2) = 1;
 
    t_diff(2, 0, 2, 0) = 0.5;
    t_diff(2, 0, 0, 2) = 0.5;
    t_diff(0, 2, 0, 2) = 0.5;
    t_diff(0, 2, 2, 0) = 0.5;
 
    t_diff(2, 1, 2, 1) = 0.5;
    t_diff(2, 1, 1, 2) = 0.5;
    t_diff(1, 2, 1, 2) = 0.5;
    t_diff(1, 2, 2, 1) = 0.5;
  }
 
  return t_diff;
};

diff_tensor()

template<int DIM>

auto PlasticOps::diff_tensor ( FTensor::Number< DIM >  )

inline

[Lambda functions]

Examples

PlasticOpsGeneric.hpp.

Definition at line 10 of file PlasticOpsGeneric.hpp.

                                                               {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Ddg<double, DIM, DIM> t_diff;
  constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
  t_diff(i, j, k, l) = (t_kd(i, k) ^ t_kd(j, l)) / 4.;
  return t_diff;
};

equivalent_strain_dot()

template<typename T , int DIM>

auto PlasticOps::equivalent_strain_dot ( FTensor::Tensor2_symmetric< T, DIM > &  t_plastic_strain_dot )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 338 of file PlasticOpsGeneric.hpp.

                                                            {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  constexpr double A = 2. / 3;
  return std::sqrt(A * t_plastic_strain_dot(i, j) *
                   t_plastic_strain_dot(i, j)) +
         std::numeric_limits<double>::epsilon();
};

get_mat_scalar_dtensor_sym() [1/2]

auto PlasticOps::get_mat_scalar_dtensor_sym	(	MatrixDouble &	mat,
		FTensor::Number< 2 >
	)

Examples

PlasticOpsGeneric.hpp.

Definition at line 1170 of file PlasticOpsGeneric.hpp.

                                                                     {
  return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>{
      &mat(0, 0), &mat(0, 1), &mat(0, 2)};
}

get_mat_scalar_dtensor_sym() [2/2]

auto PlasticOps::get_mat_scalar_dtensor_sym	(	MatrixDouble &	mat,
		FTensor::Number< 3 >
	)

Definition at line 1175 of file PlasticOpsGeneric.hpp.

                                                                     {
  return FTensor::Tensor1<FTensor::PackPtr<double *, 6>, 6>{
      &mat(0, 0), &mat(0, 1), &mat(0, 2), &mat(0, 3), &mat(0, 4), &mat(0, 5)};
}

get_mat_tensor_sym_dscalar() [1/2]

static auto PlasticOps::get_mat_tensor_sym_dscalar ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 2 >    )

inlinestatic

Examples

PlasticOpsGeneric.hpp.

Definition at line 1076 of file PlasticOpsGeneric.hpp.

                                                                {
  return FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 3>{
      &mat(3 * rr + 0, 0), &mat(3 * rr + 1, 0), &mat(3 * rr + 2, 0)};
}

get_mat_tensor_sym_dscalar() [2/2]

static auto PlasticOps::get_mat_tensor_sym_dscalar ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 3 >    )

inlinestatic

Definition at line 1082 of file PlasticOpsGeneric.hpp.

                                                                {
  return FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 6>{
      &mat(6 * rr + 0, 0), &mat(6 * rr + 1, 0), &mat(6 * rr + 2, 0),
      &mat(6 * rr + 3, 0), &mat(6 * rr + 4, 0), &mat(6 * rr + 5, 0)};
}

get_mat_tensor_sym_dtensor_sym() [1/2]

static auto PlasticOps::get_mat_tensor_sym_dtensor_sym ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 2 >    )

inlinestatic

Examples

PlasticOpsGeneric.hpp.

Definition at line 944 of file PlasticOpsGeneric.hpp.

                                                                    {
  return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 3, 3>{
      &mat(3 * rr + 0, 0), &mat(3 * rr + 0, 1), &mat(3 * rr + 0, 2),
      &mat(3 * rr + 1, 0), &mat(3 * rr + 1, 1), &mat(3 * rr + 1, 2),
      &mat(3 * rr + 2, 0), &mat(3 * rr + 2, 1), &mat(3 * rr + 2, 2)};
}

get_mat_tensor_sym_dtensor_sym() [2/2]

static auto PlasticOps::get_mat_tensor_sym_dtensor_sym ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 3 >    )

inlinestatic

Definition at line 952 of file PlasticOpsGeneric.hpp.

                                                                    {
  return FTensor::Tensor2<FTensor::PackPtr<double *, 6>, 6, 6>{
      &mat(6 * rr + 0, 0), &mat(6 * rr + 0, 1), &mat(6 * rr + 0, 2),
      &mat(6 * rr + 0, 3), &mat(6 * rr + 0, 4), &mat(6 * rr + 0, 5),
      &mat(6 * rr + 1, 0), &mat(6 * rr + 1, 1), &mat(6 * rr + 1, 2),
      &mat(6 * rr + 1, 3), &mat(6 * rr + 1, 4), &mat(6 * rr + 1, 5),
      &mat(6 * rr + 2, 0), &mat(6 * rr + 2, 1), &mat(6 * rr + 2, 2),
      &mat(6 * rr + 2, 3), &mat(6 * rr + 2, 4), &mat(6 * rr + 2, 5),
      &mat(6 * rr + 3, 0), &mat(6 * rr + 3, 1), &mat(6 * rr + 3, 2),
      &mat(6 * rr + 3, 3), &mat(6 * rr + 3, 4), &mat(6 * rr + 3, 5),
      &mat(6 * rr + 4, 0), &mat(6 * rr + 4, 1), &mat(6 * rr + 4, 2),
      &mat(6 * rr + 4, 3), &mat(6 * rr + 4, 4), &mat(6 * rr + 4, 5),
      &mat(6 * rr + 5, 0), &mat(6 * rr + 5, 1), &mat(6 * rr + 5, 2),
      &mat(6 * rr + 5, 3), &mat(6 * rr + 5, 4), &mat(6 * rr + 5, 5)};
}

get_mat_tensor_sym_dvector() [1/2]

static auto PlasticOps::get_mat_tensor_sym_dvector ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 2 >    )

inlinestatic

[Lambda functions]

[Auxiliary functions functions

Examples

PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, and PlasticOpsSmallStrains.hpp.

Definition at line 366 of file PlasticOpsGeneric.hpp.

                                                                {
  return FTensor::Tensor2<FTensor::PackPtr<double *, 2>, 3, 2>{
      &mat(3 * rr + 0, 0), &mat(3 * rr + 0, 1), &mat(3 * rr + 1, 0),
      &mat(3 * rr + 1, 1), &mat(3 * rr + 2, 0), &mat(3 * rr + 2, 1)};
}

get_mat_tensor_sym_dvector() [2/2]

static auto PlasticOps::get_mat_tensor_sym_dvector ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 3 >    )

inlinestatic

Definition at line 373 of file PlasticOpsGeneric.hpp.

                                                                {
  return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 6, 3>{
      &mat(6 * rr + 0, 0), &mat(6 * rr + 0, 1), &mat(6 * rr + 0, 2),
      &mat(6 * rr + 1, 0), &mat(6 * rr + 1, 1), &mat(6 * rr + 1, 2),
      &mat(6 * rr + 2, 0), &mat(6 * rr + 2, 1), &mat(6 * rr + 2, 2),
      &mat(6 * rr + 3, 0), &mat(6 * rr + 3, 1), &mat(6 * rr + 3, 2),
      &mat(6 * rr + 4, 0), &mat(6 * rr + 4, 1), &mat(6 * rr + 4, 2),
      &mat(6 * rr + 5, 0), &mat(6 * rr + 5, 1), &mat(6 * rr + 5, 2)};
}

get_mat_vector_dtensor_sym() [1/2]

static FTensor::Tensor2< FTensor::PackPtr< double *, 3 >, 2, 3 > PlasticOps::get_mat_vector_dtensor_sym ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 2 >    )

inlinestatic

Examples

PlasticOpsLargeStrains.hpp, and PlasticOpsSmallStrains.hpp.

Definition at line 38 of file PlasticOpsSmallStrains.hpp.

                                                                           {
  return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 2, 3>{
      &mat(2 * rr + 0, 0), &mat(2 * rr + 0, 1), &mat(2 * rr + 0, 2),
      &mat(2 * rr + 1, 0), &mat(2 * rr + 1, 1), &mat(2 * rr + 1, 2)};
}

get_mat_vector_dtensor_sym() [2/2]

static FTensor::Tensor2< FTensor::PackPtr< double *, 6 >, 3, 6 > PlasticOps::get_mat_vector_dtensor_sym ( size_t  rr, MatrixDouble &  mat, FTensor::Number< 3 >    )

inlinestatic

Definition at line 45 of file PlasticOpsSmallStrains.hpp.

                                                                           {
  return FTensor::Tensor2<FTensor::PackPtr<double *, 6>, 3, 6>{
      &mat(3 * rr + 0, 0), &mat(3 * rr + 0, 1), &mat(3 * rr + 0, 2),
      &mat(3 * rr + 0, 3), &mat(3 * rr + 0, 4), &mat(3 * rr + 0, 5),
 
      &mat(3 * rr + 1, 0), &mat(3 * rr + 1, 1), &mat(3 * rr + 1, 2),
      &mat(3 * rr + 1, 3), &mat(3 * rr + 1, 4), &mat(3 * rr + 1, 5),
 
      &mat(3 * rr + 2, 0), &mat(3 * rr + 2, 1), &mat(3 * rr + 2, 2),
      &mat(3 * rr + 2, 3), &mat(3 * rr + 2, 4), &mat(3 * rr + 2, 5)};
}

opFactoryDomainLhs()

template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >

MoFEMErrorCode PlasticOps::opFactoryDomainLhs	(	MoFEM::Interface &	m_field,
		std::string	block_name,
		Pip &	pip,
		std::string	u,
		std::string	ep,
		std::string	tau
	)

Examples

PlasticOps.hpp.

Definition at line 547 of file PlasticOps.hpp.

                                                             {
  MoFEMFunctionBegin;
 
  using namespace HenckyOps;
 
  using B = typename FormsIntegrators<DomainEleOp>::template Assembly<
      A>::template BiLinearForm<I>;
  using OpKPiola = typename B::template OpGradTensorGrad<1, DIM, DIM, 1>;
  using OpKCauchy = typename B::template OpGradSymTensorGrad<1, DIM, DIM, 0>;
 
  using P = PlasticityIntegrators<DomainEleOp>;
 
  auto [common_plastic_ptr, common_henky_ptr] =
      createCommonPlasticOps<DIM, I, DomainEleOp>(m_field, block_name, pip, u,
                                                  ep, tau, scale, Sev::verbose);
 
  auto m_D_ptr = common_plastic_ptr->mDPtr;
 
  pip.push_back(new OpCalculateTensor2SymmetricFieldValuesDot<DIM>(
      ep, common_plastic_ptr->getPlasticStrainDotPtr()));
  pip.push_back(new OpCalculateScalarFieldValuesDot(
      tau, common_plastic_ptr->getPlasticTauDotPtr()));
  pip.push_back(new typename P::template OpCalculatePlasticity<DIM, I>(
      u, common_plastic_ptr, m_D_ptr));
 
  if (common_henky_ptr) {
    using H = HenckyOps::HenkyIntegrators<DomainEleOp>;
    pip.push_back(new typename H::template OpHenckyTangent<DIM, I>(
        u, common_henky_ptr, m_D_ptr));
    pip.push_back(new OpKPiola(u, u, common_henky_ptr->getMatTangent()));
    pip.push_back(
        new typename P::template Assembly<A>::
            template OpCalculatePlasticInternalForceLhs_LogStrain_dEP<DIM, I>(
                u, ep, common_plastic_ptr, common_henky_ptr, m_D_ptr));
  } else {
    pip.push_back(new OpKCauchy(u, u, m_D_ptr));
    pip.push_back(new typename P::template Assembly<A>::
                      template OpCalculatePlasticInternalForceLhs_dEP<DIM, I>(
                          u, ep, common_plastic_ptr, m_D_ptr));
  }
 
  if (common_henky_ptr) {
    pip.push_back(
        new typename P::template Assembly<A>::
            template OpCalculateConstraintsLhs_LogStrain_dU<DIM, I>(
                tau, u, common_plastic_ptr, common_henky_ptr, m_D_ptr));
    pip.push_back(
        new typename P::template Assembly<A>::
            template OpCalculatePlasticFlowLhs_LogStrain_dU<DIM, I>(
                ep, u, common_plastic_ptr, common_henky_ptr, m_D_ptr));
  } else {
    pip.push_back(
        new
        typename P::template Assembly<A>::template OpCalculateConstraintsLhs_dU<
            DIM, I>(tau, u, common_plastic_ptr, m_D_ptr));
    pip.push_back(
        new
        typename P::template Assembly<A>::template OpCalculatePlasticFlowLhs_dU<
            DIM, I>(ep, u, common_plastic_ptr, m_D_ptr));
  }
 
  pip.push_back(
      new
      typename P::template Assembly<A>::template OpCalculatePlasticFlowLhs_dEP<
          DIM, I>(ep, ep, common_plastic_ptr, m_D_ptr));
  pip.push_back(
      new
      typename P::template Assembly<A>::template OpCalculatePlasticFlowLhs_dTAU<
          DIM, I>(ep, tau, common_plastic_ptr, m_D_ptr));
  pip.push_back(
      new
      typename P::template Assembly<A>::template OpCalculateConstraintsLhs_dEP<
          DIM, I>(tau, ep, common_plastic_ptr, m_D_ptr));
  pip.push_back(
      new
      typename P::template Assembly<A>::template OpCalculateConstraintsLhs_dTAU<
          I>(tau, tau, common_plastic_ptr));
 
  MoFEMFunctionReturn(0);
}

opFactoryDomainReactions()

template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >

MoFEMErrorCode PlasticOps::opFactoryDomainReactions	(	MoFEM::Interface &	m_field,
		std::string	block_name,
		Pip &	pip,
		std::string	u,
		std::string	ep,
		std::string	tau
	)

Examples

PlasticOps.hpp.

Definition at line 630 of file PlasticOps.hpp.

                                                       {
  MoFEMFunctionBegin;
 
  using B = typename FormsIntegrators<DomainEleOp>::template Assembly<
      A>::template LinearForm<I>;
  using OpInternalForceCauchy =
      typename B::template OpGradTimesSymTensor<1, DIM, DIM>;
  using OpInternalForcePiola =
      typename B::template OpGradTimesTensor<1, DIM, DIM>;
 
  auto [common_plastic_ptr, common_henky_ptr] =
      createCommonPlasticOps<DIM, I, DomainEleOp>(m_field, block_name, pip, u,
                                                  ep, tau, 1., Sev::inform);
 
  // Calculate internal forces
  if (common_henky_ptr) {
    pip.push_back(new OpInternalForcePiola(
        u, common_henky_ptr->getMatFirstPiolaStress()));
  } else {
    pip.push_back(new OpInternalForceCauchy(u, common_plastic_ptr->mStressPtr));
  }
 
  MoFEMFunctionReturn(0);
}

opFactoryDomainRhs()

template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >

MoFEMErrorCode PlasticOps::opFactoryDomainRhs	(	MoFEM::Interface &	m_field,
		std::string	block_name,
		Pip &	pip,
		std::string	u,
		std::string	ep,
		std::string	tau
	)

Examples

PlasticOps.hpp.

Definition at line 500 of file PlasticOps.hpp.

                                                             {
  MoFEMFunctionBegin;
 
  using B = typename FormsIntegrators<DomainEleOp>::template Assembly<
      A>::template LinearForm<I>;
  using OpInternalForceCauchy =
      typename B::template OpGradTimesSymTensor<1, DIM, DIM>;
  using OpInternalForcePiola =
      typename B::template OpGradTimesTensor<1, DIM, DIM>;
 
  using P = PlasticityIntegrators<DomainEleOp>;
 
  auto [common_plastic_ptr, common_henky_ptr] =
      createCommonPlasticOps<DIM, I, DomainEleOp>(m_field, block_name, pip, u,
                                                  ep, tau, scale, Sev::inform);
 
  auto m_D_ptr = common_plastic_ptr->mDPtr;
 
  pip.push_back(new OpCalculateTensor2SymmetricFieldValuesDot<DIM>(
      ep, common_plastic_ptr->getPlasticStrainDotPtr()));
  pip.push_back(new OpCalculateScalarFieldValuesDot(
      tau, common_plastic_ptr->getPlasticTauDotPtr()));
  pip.push_back(new typename P::template OpCalculatePlasticity<DIM, I>(
      u, common_plastic_ptr, m_D_ptr));
 
  // Calculate internal forces
  if (common_henky_ptr) {
    pip.push_back(new OpInternalForcePiola(
        u, common_henky_ptr->getMatFirstPiolaStress()));
  } else {
    pip.push_back(new OpInternalForceCauchy(u, common_plastic_ptr->mStressPtr));
  }
 
  pip.push_back(
      new
      typename P::template Assembly<A>::template OpCalculateConstraintsRhs<I>(
          tau, common_plastic_ptr, m_D_ptr));
  pip.push_back(
      new typename P::template Assembly<A>::template OpCalculatePlasticFlowRhs<
          DIM, I>(ep, common_plastic_ptr, m_D_ptr));
 
  MoFEMFunctionReturn(0);
}

plastic_flow()

template<int DIM>

auto PlasticOps::plastic_flow ( long double  f, FTensor::Tensor2_symmetric< double, 3 > &&  t_dev_stress, FTensor::Ddg< double, 3, DIM > &&  t_diff_deviator  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 195 of file PlasticOpsGeneric.hpp.

                                                                       {
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Tensor2_symmetric<double, DIM> t_diff_f;
  t_diff_f(k, l) =
      (1.5 * (t_dev_stress(I, J) * t_diff_deviator(I, J, k, l))) / f;
  return t_diff_f;
};

platsic_surface()

double PlasticOps::platsic_surface ( FTensor::Tensor2_symmetric< double, 3 > &&  t_stress_deviator )

inline

\[ \begin{split} f&=\sqrt{s_{ij}s_{ij}}\\ A_{ij}&=\frac{\partial f}{\partial \sigma_{ij}}= \frac{1}{f} s_{kl} \frac{\partial s_{kl}}{\partial \sigma_{ij}}\\ \frac{\partial A_{ij}}{\partial \sigma_{kl}}&= \frac{\partial^2 f}{\partial \sigma_{ij}\partial\sigma_{mn}}= \frac{1}{f} \left( \frac{\partial s_{kl}}{\partial \sigma_{mn}}\frac{\partial s_{kl}}{\partial \sigma_{ij}} -A_{mn}A_{ij} \right)\\ \frac{\partial f}{\partial \varepsilon_{ij}}&=A_{mn}D_{mnij} \\ \frac{\partial A_{ij}}{\partial \varepsilon_{kl}}&= \frac{\partial A_{ij}}{\partial \sigma_{mn}} \frac{\partial \sigma_{mn}}{\partial \varepsilon_{kl}}= \frac{\partial A_{ij}}{\partial \sigma_{mn}} D_{mnkl} \end{split} \]

Examples

PlasticOpsGeneric.hpp.

Definition at line 189 of file PlasticOpsGeneric.hpp.

                                                                         {
  return std::sqrt(1.5 * t_stress_deviator(I, J) * t_stress_deviator(I, J)) +
         std::numeric_limits<double>::epsilon();
};

symm_L_tensor()

template<int DIM>

auto PlasticOps::symm_L_tensor ( FTensor::Number< DIM >  )

inline

Examples

PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, and PlasticOpsSmallStrains.hpp.

Definition at line 21 of file PlasticOpsGeneric.hpp.

                                                                 {
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
  FTensor::Index<'L', size_symm> L;
  FTensor::Dg<double, DIM, size_symm> t_L;
  t_L(i, j, L) = 0;
  if constexpr (DIM == 2) {
    t_L(0, 0, 0) = 1;
    t_L(1, 0, 1) = 1;
    t_L(1, 1, 2) = 1;
  } else {
    t_L(0, 0, 0) = 1;
    t_L(1, 0, 1) = 1;
    t_L(2, 0, 2) = 1;
    t_L(1, 1, 3) = 1;
    t_L(2, 1, 4) = 1;
    t_L(2, 2, 5) = 1;
  }
  return t_L;
}

trace() [1/2]

template<typename T >

double PlasticOps::trace ( FTensor::Tensor2_symmetric< T, 2 > &  t_stress )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 80 of file PlasticOpsGeneric.hpp.

                                                              {
  constexpr double third = boost::math::constants::third<double>();
  return (t_stress(0, 0) + t_stress(1, 1)) * third;
};

trace() [2/2]

template<typename T >

double PlasticOps::trace ( FTensor::Tensor2_symmetric< T, 3 > &  t_stress )

inline

Definition at line 86 of file PlasticOpsGeneric.hpp.

                                                              {
  constexpr double third = boost::math::constants::third<double>();
  return (t_stress(0, 0) + t_stress(1, 1) + t_stress(2, 2)) * third;
};

w()

double PlasticOps::w ( double  eqiv, double  dot_tau, double  f, double  sigma_y, double  sigma_Y  )

inline

Examples

PlasticOpsGeneric.hpp.

Definition at line 276 of file PlasticOpsGeneric.hpp.

                                {
  return (f - sigma_y) / sigma_Y + cn1 * (eqiv);
};

Variable Documentation

I

FTensor::Index<'I', 3> PlasticOps::I

[Common data]

Examples

PlasticOps.hpp, and PlasticOpsGeneric.hpp.

Definition at line 115 of file PlasticOps.hpp.

J

FTensor::Index<'J', 3> PlasticOps::J

Examples

PlasticOps.hpp, and PlasticOpsGeneric.hpp.

Definition at line 116 of file PlasticOps.hpp.

M

FTensor::Index<'M', 3> PlasticOps::M

Examples

PlasticOps.hpp, and PlasticOpsGeneric.hpp.

Definition at line 117 of file PlasticOps.hpp.

N

FTensor::Index<'N', 3> PlasticOps::N

Examples

PlasticOps.hpp, and PlasticOpsGeneric.hpp.

Definition at line 118 of file PlasticOps.hpp.