ADOLCPlasticity Namespace Reference

Classes
struct	ADOLCPlasticityIntegrators
struct	ClosestPointProjection
	Closest point projection algorithm. More...
struct	CommonData
	common data used by volume elements More...
struct	J2Plasticity
	J2 plasticity (Kinematic Isotropic (Linear) Hardening) More...
struct	J2Plasticity< 2 >
	[J2 2D] More...
struct	J2Plasticity< 3 >
	J2 (Von Misses) plasticity. More...
struct	MakeB
	Calculate tensorial base functions. Apply bBar method when needed. More...
struct	OpCalculateStress
struct	OpLhsImpl
	Assemble left hand side. More...
struct	OpLhsImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpRhsImpl
	Assemble right hand side. More...
struct	OpRhsImpl< DIM, GAUSS, AssemblyDomainEleOp >
struct	OpUpdate
struct	ParaboloidalPlasticity
	[J2 2D] More...
struct	TSUpdate
	Update internal fluxes (update history variables) More...
struct	TSUpdateImpl
	[TSUpdateImpl] More...

Typedefs
using	Pip = boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator >

Functions
FTensor::Dg< double, 3, 6 >	voight_to_strain_op ()
	Op convert Vight strain vector to strain tensor. More...
FTensor::Dg< double, 3, 6 >	strain_to_voight_op ()
	Op convert strain tensor to Vight strain vector. More...
FTensor::Dg< double, 3, 6 >	voight_to_stress_op ()
	Op convert Vight stress vector to stress tensor. More...
MoFEMErrorCode	ADOLCPlasticityRes (SNES snes, Vec chi, Vec r, void *ctx)
	Internal SNES function used at integration points to calulate stress. More...
MoFEMErrorCode	ADOLCPlasticityJac (SNES, Vec, Mat, Mat, void *ctx)
	Internal SNES function used at integration points to calulate tangent matrix. More...
template<int DIM>
ForcesAndSourcesCore::UserDataOperator *	getRawPtrOpCalculateStress (MoFEM::Interface &m_field, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr, bool calc_lhs)
	Get opreator to calulate stress. More...
template<>
ForcesAndSourcesCore::UserDataOperator *	getRawPtrOpCalculateStress< 3 > (MoFEM::Interface &m_field, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr, bool calc_lhs)
template<>
ForcesAndSourcesCore::UserDataOperator *	getRawPtrOpCalculateStress< 2 > (MoFEM::Interface &m_field, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr, bool calc_lhs)
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >
MoFEMErrorCode	opFactoryDomainRhs (MoFEM::Interface &m_field, string field_name, Pip &pip, std::string block_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr, Sev sev=Sev::inform)
	Assemble the left hand side, i.e. tangent matrix. More...
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >
MoFEMErrorCode	opFactoryDomainLhs (MoFEM::Interface &m_field, string field_name, Pip &pip, std::string block_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr)
	Assemble the left hand side, i.e. tangent matrix. More...
template<int DIM>
MoFEMErrorCode	opFactoryDomainUpdate (MoFEM::Interface &m_field, Pip &pip, std::string block_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr)
	Push operators to update history variables. More...
template<>
MoFEMErrorCode	opFactoryDomainUpdate< 3 > (MoFEM::Interface &m_field, Pip &pip, std::string block_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr)
	Push operators to update history variables. More...
template<>
MoFEMErrorCode	opFactoryDomainUpdate< 2 > (MoFEM::Interface &m_field, Pip &pip, std::string block_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr)
	Push operators to update history variables. More...
boost::shared_ptr< TSUpdate >	createTSUpdate (std::string fe_name, boost::shared_ptr< FEMethod > fe_ptr)
template<typename MODEL >
auto	createMaterial (std::array< int, ClosestPointProjection::LAST_TAPE > tapes_tags={0, 1, 2})
template<typename T1 , typename T2 , int DIM1, int DIM2>
double	calcAveStrain (bool b_bar, const int nb_gauss_pts, FTensor::Tensor2< T1, DIM1, DIM2 > &&t_grad, FTensor::Tensor0< T2 > &&t_w)
template<typename T1 , typename T2 , int DIM1, int DIM2, int DIM3>
FTensor::Tensor1< T2, DIM3 >	calcStrain (bool b_bar, double ave_tr_strain, FTensor::Tensor2< T1, DIM1, DIM2 > &t_grad, FTensor::Tensor1< T2, DIM3 > &&t_voight_strain)
template<int DIM>
FTensor::Tensor2_symmetric< FTensor::PackPtr< double *, 6 >, 3 >	getFTensor2SymmetricDiffBaseImpl (DataForcesAndSourcesCore::EntData &data, MatrixDouble &storage, const bool b_bar, const int nb_integration_pts, double *w_ptr, FTensor::Number< DIM >)
	[BBar method] More...
template<int DIM>
MoFEMErrorCode	opFactoryDomainUpdateImpl (MoFEM::Interface &m_field, Pip &pip, std::string block_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< ClosestPointProjection > cp_ptr)
boost::shared_ptr< TSUpdate >	createTSUpdate (std::string fe_name, boost::shared_ptr< FEMethod > fe_ptr)
	Create test update object. More...

Variables
static FTensor::Dg< double, 3, 6 >	tStrainToVoightOp = strain_to_voight_op()

Typedef Documentation

Pip

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

Examples

ADOLCPlasticity.hpp.

Definition at line 399 of file ADOLCPlasticity.hpp.

Function Documentation

ADOLCPlasticityJac()

MoFEMErrorCode ADOLCPlasticity::ADOLCPlasticityJac	(	SNES	snes,
		Vec	chi,
		Mat	A,
		Mat	,
		void *	ctx
	)

Internal SNES function used at integration points to calulate tangent matrix.

Function executed by nested SENES at evaluationg tangent matrix.

Examples

ADOLCPlasticity.hpp.

Definition at line 584 of file ClosetPointProjection.cpp.

                                                    {
  MoFEMFunctionBegin;
  ClosestPointProjection *cp;
  cp = (ClosestPointProjection *)ctx;
 
  auto plastic_strain = cp->getPlasticStrain();
  auto internal_variables = cp->getInternalVariables();
 
  {
    const double *array;
    CHKERR VecGetArrayRead(chi, &array);
    for (auto ii = 0; ii < 6; ii++) {
      plastic_strain[ii] = array[ii];
    }
    for (auto ii = 0; ii < cp->nbInternalVariables; ii++) {
      internal_variables[ii] = array[6 + ii];
    }
    cp->deltaGamma = array[6 + cp->nbInternalVariables];
    CHKERR VecRestoreArrayRead(chi, &array);
  }
  CHKERR cp->playPotentials();
  CHKERR cp->calculateA();
  MoFEMFunctionReturn(0);
}

ADOLCPlasticityRes()

MoFEMErrorCode ADOLCPlasticity::ADOLCPlasticityRes	(	SNES	snes,
		Vec	chi,
		Vec	r,
		void *	ctx
	)

Internal SNES function used at integration points to calulate stress.

Function executed by nested SENES at evaluationg residual.

Examples

ADOLCPlasticity.hpp.

Definition at line 558 of file ClosetPointProjection.cpp.

                                                                        {
  MoFEMFunctionBegin;
  ClosestPointProjection *cp;
  cp = (ClosestPointProjection *)ctx;
 
  auto plastic_strain = cp->getPlasticStrain();
  auto internal_variables = cp->getInternalVariables();
 
  {
    const double *array;
    CHKERR VecGetArrayRead(chi, &array);
    for (auto ii = 0; ii < 6; ii++) {
      plastic_strain[ii] = array[ii];
    }
    for (auto ii = 0; ii < cp->nbInternalVariables; ii++) {
      internal_variables[ii] = array[6 + ii];
    }
    cp->deltaGamma = array[6 + cp->nbInternalVariables];
    CHKERR VecRestoreArrayRead(chi, &array);
  }
  CHKERR cp->playPotentials_NoHessian();
  CHKERR cp->calculateR(r);
 
  MoFEMFunctionReturn(0);
}

calcAveStrain()

template<typename T1 , typename T2 , int DIM1, int DIM2>

double ADOLCPlasticity::calcAveStrain	(	bool	b_bar,
		const int	nb_gauss_pts,
		FTensor::Tensor2< T1, DIM1, DIM2 > &&	t_grad,
		FTensor::Tensor0< T2 > &&	t_w
	)

Definition at line 18 of file ADOLCPlasticity.cpp.

                                                      {
  if constexpr (DIM1 != DIM2)
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "Case of mixed dimension by gradient not implemented");
 
 
  FTensor::Index<'i', DIM1> i;
  double ave_tr_strain = 0;
  if (b_bar) {
    double v = 0;
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
      v += t_w;
      ave_tr_strain += t_w * t_grad(i, i);
      ++t_grad;
      ++t_w;
    }
    ave_tr_strain /= (DIM1 * v);
  }
  return ave_tr_strain;
}

calcStrain()

template<typename T1 , typename T2 , int DIM1, int DIM2, int DIM3>

FTensor::Tensor1<T2, DIM3> ADOLCPlasticity::calcStrain	(	bool	b_bar,
		double	ave_tr_strain,
		FTensor::Tensor2< T1, DIM1, DIM2 > &	t_grad,
		FTensor::Tensor1< T2, DIM3 > &&	t_voight_strain
	)

Definition at line 43 of file ADOLCPlasticity.cpp.

                                                       {
 
  if constexpr (DIM1 != DIM2)
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "Case of mixed dimessnion by gradient not implemented");
 
  FTensor::Index<'i', DIM1> i;
  FTensor::Index<'j', DIM2> j;
  FTensor::Index<'Z', DIM3> Z;
  FTensor::Tensor2_symmetric<double, 3> t_strain{0., 0., 0., 0., 0., 0.};
  t_strain(i, j) = (t_grad(i, j) || t_grad(j, i)) / 2;
  double tr_strain = t_grad(i, i) / DIM1;
  if (b_bar) {
    FTensor::Index<'i', DIM1> i;
    FTensor::Index<'j', DIM2> j;
    auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
    t_strain(i, j) += (ave_tr_strain - tr_strain) * t_kd(i, j);
  } 
  {
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    t_voight_strain(Z) = tStrainToVoightOp(i, j, Z) * t_strain(i, j);
  }
  return t_voight_strain;
};

createMaterial()

template<typename MODEL >

auto ADOLCPlasticity::createMaterial ( std::array< int, ClosestPointProjection::LAST_TAPE >  tapes_tags = {0, 1, 2} )

inline

Examples

ADOLCPlasticityMaterialModels.hpp.

Definition at line 721 of file ADOLCPlasticityMaterialModels.hpp.

                                                                {0, 1, 2}) {
  auto cp_ptr = boost::make_shared<MODEL>();
  cp_ptr->getDefaultMaterialParameters();
  cp_ptr->createMatAVecR();
  cp_ptr->snesCreate();
  cp_ptr->recordTapes();
  cp_ptr->tapesTags = tapes_tags;
  return cp_ptr;
}

createTSUpdate() [1/2]

boost::shared_ptr<TSUpdate> ADOLCPlasticity::createTSUpdate	(	std::string	fe_name,
		boost::shared_ptr< FEMethod >	fe_ptr
	)

Examples

adolc_plasticity.cpp, and ADOLCPlasticity.hpp.

createTSUpdate() [2/2]

boost::shared_ptr<TSUpdate> ADOLCPlasticity::createTSUpdate	(	std::string	fe_name,
		boost::shared_ptr< FEMethod >	fe_ptr
	)

Create test update object.

Parameters

fe_name
fe_ptr

Returns

boost::shared_ptr<TSUpdate>

Definition at line 630 of file ADOLCPlasticity.cpp.

                                                                             {
  return boost::make_shared<TSUpdateImpl>(fe_name, fe_ptr);
}

getFTensor2SymmetricDiffBaseImpl()

template<int DIM>

FTensor::Tensor2_symmetric<FTensor::PackPtr<double *, 6>, 3> ADOLCPlasticity::getFTensor2SymmetricDiffBaseImpl	(	DataForcesAndSourcesCore::EntData &	data,
		MatrixDouble &	storage,
		const bool	b_bar,
		const int	nb_integration_pts,
		double *	w_ptr,
		FTensor::Number< DIM >
	)

[BBar method]

Calculate base for B-bar method

Template Parameters

DIM	Dimension of problem

Parameters

data	access to base functions
storage	storage for base functions at integration points
b_bar	flag to sith on B-bar method
nb_integration_pts	number of integration points
w_ptr	integration weights

Returns

FTensor::Tensor2_symmetric<FTensor::PackPtr<double *, 6>, 3>

Definition at line 85 of file ADOLCPlasticity.cpp.

                                                     {
  const auto nb_dofs = data.getFieldData().size();
  const auto nb_bases = data.getN().size2();
 
  storage.resize(nb_integration_pts, 6 * nb_dofs, false);
 
  FTensor::Index<'I', DIM> I;
  FTensor::Index<'J', DIM> J;
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Number<0> N0;
  FTensor::Number<1> N1;
  FTensor::Number<2> N2;
 
  // get tensorial base
  auto get_ftensor2_symmetric = [&](const int gg, const int rr) {
    return getFTensor2SymmetricFromPtr<3>(&storage(gg, 6 * rr));
  };
 
  // calculate base 
  auto calc_base = [&]() {
    auto t_t2_diff = get_ftensor2_symmetric(0, 0);
    for (auto gg = 0; gg != nb_integration_pts; ++gg) {
      auto t_diff = data.getFTensor1DiffN<DIM>(gg, 0);
      for (auto b = 0; b != nb_dofs / DIM; ++b) {
        FTensor::Tensor2<double, 3, 3> t_grad;
 
        t_grad(i, j) = 0;
 
        t_grad(N0, J) = t_diff(J);
        t_t2_diff(i, j) = (t_grad(i, j) || t_grad(j, i)) / 2;
        ++t_t2_diff;
 
        t_grad(i, j) = 0;
        t_grad(N1, J) = t_diff(J);
        t_t2_diff(i, j) = (t_grad(i, j) || t_grad(j, i)) / 2;
        ++t_t2_diff;
 
        t_grad(i, j) = 0;
        if constexpr (DIM == 3) {
          t_grad(N2, J) = t_diff(J);
          t_t2_diff(i, j) = (t_grad(i, j) || t_grad(j, i)) / 2;
          ++t_t2_diff;
        }
 
        ++t_diff;
      }
    }
    return get_ftensor2_symmetric(0, 0);
  };
 
  // calculate volume average of trace
  auto calc_vol = [&](auto &&t_t2_dff) {
    std::vector<double> vol(nb_dofs, 0);
    for (auto gg = 0; gg != nb_integration_pts; ++gg) {
      for (auto b = 0; b != nb_dofs; ++b) {
        vol[b] += w_ptr[gg] * t_t2_dff(i, i) / DIM;
        ++t_t2_dff;
      }
    }
    double sum = 0;
    for (auto gg = 0; gg != nb_integration_pts; ++gg) {
      sum += w_ptr[gg];
    }
    for (auto &v : vol) {
      v /= sum;
    }
    return vol;
  };
 
  // modify base for B-bar method
  auto make_b_bar = [&](auto &&vol) {
    auto t_t2_diff = get_ftensor2_symmetric(0, 0);
    constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
    for (auto gg = 0; gg != nb_integration_pts; ++gg) {
      for (auto b = 0; b != nb_dofs; ++b) {
        const auto trace = t_t2_diff(J, J) / DIM;
        t_t2_diff(I, J) += (vol[b] - trace) * t_kd(I, J);
        ++t_t2_diff;
      }
    }
    return get_ftensor2_symmetric(0, 0);
  };
 
  return b_bar ? make_b_bar(calc_vol(calc_base())) : calc_base();
};

getRawPtrOpCalculateStress()

template<int DIM>

ForcesAndSourcesCore::UserDataOperator* ADOLCPlasticity::getRawPtrOpCalculateStress	(	MoFEM::Interface &	m_field,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< ClosestPointProjection >	cp_ptr,
		bool	calc_lhs
	)

Get opreator to calulate stress.

Examples

ADOLCPlasticity.hpp.

getRawPtrOpCalculateStress< 2 >()

template<>

ForcesAndSourcesCore::UserDataOperator * ADOLCPlasticity::getRawPtrOpCalculateStress< 2 >	(	MoFEM::Interface &	m_field,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< ClosestPointProjection >	cp_ptr,
		bool	calc_lhs
	)

Examples

ADOLCPlasticity.hpp.

Definition at line 554 of file ADOLCPlasticity.cpp.

                                                                   {
  return new OpCalculateStress<2>(m_field, common_data_ptr, cp_ptr, calc_lhs);
}

getRawPtrOpCalculateStress< 3 >()

template<>

ForcesAndSourcesCore::UserDataOperator * ADOLCPlasticity::getRawPtrOpCalculateStress< 3 >	(	MoFEM::Interface &	m_field,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< ClosestPointProjection >	cp_ptr,
		bool	calc_lhs
	)

Examples

ADOLCPlasticity.hpp.

Definition at line 547 of file ADOLCPlasticity.cpp.

                                                                   {
  return new OpCalculateStress<3>(m_field, common_data_ptr, cp_ptr, calc_lhs);
}

opFactoryDomainUpdate< 2 >()

template<>

MoFEMErrorCode ADOLCPlasticity::opFactoryDomainUpdate< 2 >	(	MoFEM::Interface &	m_field,
		Pip &	pip,
		std::string	block_name,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< ClosestPointProjection >	cp_ptr
	)

Push operators to update history variables.

Template Parameters

DIM	dimension of the problem

Parameters

m_field	core interface
pip
block_name	mesh block name caring material parameters
common_data_ptr
cp_ptr

Returns

MoFEMErrorCode

Examples

ADOLCPlasticity.hpp.

Definition at line 586 of file ADOLCPlasticity.cpp.

                                                                         {
  return opFactoryDomainUpdateImpl<2>(m_field, pip, block_name, common_data_ptr,
                                      cp_ptr);
};

opFactoryDomainUpdate< 3 >()

template<>

MoFEMErrorCode ADOLCPlasticity::opFactoryDomainUpdate< 3 >	(	MoFEM::Interface &	m_field,
		Pip &	pip,
		std::string	block_name,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< ClosestPointProjection >	cp_ptr
	)

Push operators to update history variables.

Template Parameters

DIM	dimension of the problem

Parameters

m_field	core interface
pip
block_name	mesh block name caring material parameters
common_data_ptr
cp_ptr

Returns

MoFEMErrorCode

Examples

ADOLCPlasticity.hpp.

Definition at line 576 of file ADOLCPlasticity.cpp.

                                                                         {
  return opFactoryDomainUpdateImpl<3>(m_field, pip, block_name, common_data_ptr,
                                      cp_ptr);
};

opFactoryDomainUpdateImpl()

template<int DIM>

MoFEMErrorCode ADOLCPlasticity::opFactoryDomainUpdateImpl	(	MoFEM::Interface &	m_field,
		Pip &	pip,
		std::string	block_name,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< ClosestPointProjection >	cp_ptr
	)

Definition at line 562 of file ADOLCPlasticity.cpp.

                                                                          {
  MoFEMFunctionBegin;
  CHKERR cp_ptr->addMatBlockOps(m_field, pip, block_name, Sev::noisy);
  pip.push_back(
      getRawPtrOpCalculateStress<DIM>(m_field, common_data_ptr, cp_ptr, false));
  pip.push_back(new OpUpdate(common_data_ptr));
  MoFEMFunctionReturn(0);
};

strain_to_voight_op()

FTensor::Dg<double, 3, 6> ADOLCPlasticity::strain_to_voight_op ( )

inline

Op convert strain tensor to Vight strain vector.

Examples

ADOLCPlasticity.hpp.

Definition at line 39 of file ADOLCPlasticity.hpp.

                                                     {
  return FTensor::Dg<double, 3, 6>{1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                                   1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0,
                                   0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                                   0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0};
};

voight_to_strain_op()

FTensor::Dg<double, 3, 6> ADOLCPlasticity::voight_to_strain_op ( )

inline

Op convert Vight strain vector to strain tensor.

Examples

ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 29 of file ADOLCPlasticity.hpp.

                                                     {
  return FTensor::Dg<double, 3, 6>{1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                                   0.5, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.5, 0.0,
                                   0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                                   0.0, 0.0, 0.5, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0};
};

voight_to_stress_op()

FTensor::Dg<double, 3, 6> ADOLCPlasticity::voight_to_stress_op ( )

inline

Op convert Vight stress vector to stress tensor.

Examples

ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 49 of file ADOLCPlasticity.hpp.

                                                     {
  return FTensor::Dg<double, 3, 6>{1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                                   1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0,
                                   0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                                   0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0};
};

Variable Documentation

tStrainToVoightOp

FTensor::Dg<double, 3, 6> ADOLCPlasticity::tStrainToVoightOp = strain_to_voight_op()

static

Definition at line 15 of file ADOLCPlasticity.cpp.