ContactOps Namespace Reference

Classes
struct	CommonData
	[Common data] More...
struct	ContactIntegrators
struct	Monitor
struct	OpAssembleTotalContactTractionImpl
struct	OpAssembleTotalContactTractionImpl< DIM, GAUSS, BoundaryEleOp >
struct	OpConstrainBoundaryLhs_dTractionImpl
struct	OpConstrainBoundaryLhs_dTractionImpl< DIM, GAUSS, AssemblyBoundaryEleOp >
struct	OpConstrainBoundaryLhs_dUImpl
struct	OpConstrainBoundaryLhs_dUImpl< DIM, GAUSS, AssemblyBoundaryEleOp >
struct	OpConstrainBoundaryRhsImpl
struct	OpConstrainBoundaryRhsImpl< DIM, GAUSS, AssemblyBoundaryEleOp >
struct	OpEvaluateSDFImpl
struct	OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >
struct	OpMixLhsSide
struct	PostProcEleByDim
struct	PostProcEleByDim< 2 >
struct	PostProcEleByDim< 3 >

Typedefs
using	EntData = EntitiesFieldData::EntData
using	BoundaryEle = FaceElementForcesAndSourcesCore
using	BoundaryEleOp = BoundaryEle::UserDataOperator
using	AssemblyBoundaryEleOp = FormsIntegrators< BoundaryEleOp >::Assembly< A >::OpBase
using	SurfaceDistanceFunction = boost::function< double(double t, double x, double y, double z, double tx, double ty, double tz)>
	[Common data] More...
using	GradSurfaceDistanceFunction = boost::function< FTensor::Tensor1< double, 3 >(double t, double x, double y, double z, double tx, double ty, double tz)>
using	HessSurfaceDistanceFunction = boost::function< FTensor::Tensor2_symmetric< double, 3 >(double t, double x, double y, double z, double tx, double ty, double tz)>
using	PostProcEleDomain = PostProcEleByDim< SPACE_DIM >::PostProcEleDomain
using	SideEle = PostProcEleByDim< SPACE_DIM >::SideEle
using	PostProcEleBdy = PostProcEleByDim< SPACE_DIM >::PostProcEleBdy

Functions
double	surface_distance_function (double t, double x, double y, double z, double tx, double ty, double tz)
FTensor::Tensor1< double, 3 >	grad_surface_distance_function (double t, double x, double y, double z, double tx, double ty, double tz)
FTensor::Tensor2_symmetric< double, 3 >	hess_surface_distance_function (double t, double x, double y, double z, double tx, double ty, double tz)
double	sign (double x)
double	w (const double sdf, const double tn)
double	constrain (double sdf, double tn)
	constrain function More...
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEleOp >
MoFEMErrorCode	opFactoryDomainRhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma, std::string u)
template<int DIM, AssemblyType A, IntegrationType I, typename DomainEle >
MoFEMErrorCode	opFactoryBoundaryToDomainLhs (MoFEM::Interface &m_field, boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string fe_domain_name, std::string sigma, std::string u, std::string geom, ForcesAndSourcesCore::RuleHookFun rule)
template<int DIM, AssemblyType A, IntegrationType I, typename BoundaryEleOp >
MoFEMErrorCode	opFactoryBoundaryLhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma, std::string u)
template<int DIM, AssemblyType A, IntegrationType I, typename BoundaryEleOp >
MoFEMErrorCode	opFactoryBoundaryRhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma, std::string u)
template<int DIM, IntegrationType I, typename BoundaryEleOp >
MoFEMErrorCode	opFactoryCalculateTraction (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma)

Variables
double	cn_contact = 1
double	scale = 1

Typedef Documentation

AssemblyBoundaryEleOp

using ContactOps::AssemblyBoundaryEleOp = typedef FormsIntegrators<BoundaryEleOp>::Assembly<A>::OpBase

Examples

ContactOps.hpp.

Definition at line 16 of file EshelbianContact.hpp.

BoundaryEle

using ContactOps::BoundaryEle = typedef FaceElementForcesAndSourcesCore

Definition at line 13 of file EshelbianContact.hpp.

BoundaryEleOp

using ContactOps::BoundaryEleOp = typedef BoundaryEle::UserDataOperator

Examples

ContactOps.hpp.

Definition at line 14 of file EshelbianContact.hpp.

EntData

using ContactOps::EntData = typedef EntitiesFieldData::EntData

Definition at line 12 of file EshelbianContact.hpp.

GradSurfaceDistanceFunction

typedef boost::function< FTensor::Tensor1< double, 3 >(double t, double x, double y, double z, double tx, double ty, double tz)> ContactOps::GradSurfaceDistanceFunction

Examples

ContactOps.hpp.

Definition at line 206 of file ContactOps.hpp.

HessSurfaceDistanceFunction

typedef boost::function< FTensor::Tensor2_symmetric< double, 3 >(double t, double x, double y, double z, double tx, double ty, double tz)> ContactOps::HessSurfaceDistanceFunction

Examples

ContactOps.hpp.

Definition at line 209 of file ContactOps.hpp.

PostProcEleBdy

typedef PostProcEleByDim< SPACE_DIM >::PostProcEleBdy ContactOps::PostProcEleBdy

Definition at line 26 of file PostProcContact.hpp.

PostProcEleDomain

typedef PostProcEleByDim< SPACE_DIM >::PostProcEleDomain ContactOps::PostProcEleDomain

Definition at line 24 of file PostProcContact.hpp.

SideEle

typedef PostProcEleByDim< SPACE_DIM >::SideEle ContactOps::SideEle

Definition at line 25 of file PostProcContact.hpp.

SurfaceDistanceFunction

typedef boost::function< double(double t, double x, double y, double z, double tx, double ty, double tz)> ContactOps::SurfaceDistanceFunction

[Common data]

[Surface distance function from python] [Surface distance function from python]

Examples

ContactOps.hpp.

Definition at line 203 of file ContactOps.hpp.

Function Documentation

constrain()

double ContactOps::constrain ( double  sdf, double  tn  )

inline

constrain function

return 1 if negative sdf or positive tn

Parameters

sdf	signed distance
tn	traction

Returns

double

return 1 if negative sdn or positive tn

Parameters

sdf	signed distance
tn	traction

Returns

double

Examples

ContactOps.hpp.

Definition at line 400 of file ContactOps.hpp.

                                               {
  const auto s = sign(w(sdf, tn));
  return (1 - s) / 2;
}

grad_surface_distance_function()

FTensor::Tensor1< double, 3 > ContactOps::grad_surface_distance_function ( double  t, double  x, double  y, double  z, double  tx, double  ty, double  tz  )

inline

Examples

ContactOps.hpp.

Definition at line 229 of file ContactOps.hpp.

                                                                {
#ifdef PYTHON_SFD
  if (auto sdf_ptr = sdfPythonWeakPtr.lock()) {
    std::vector<double> grad_sdf;
    CHK_MOAB_THROW(sdf_ptr->evalGradSdf(t, x, y, z, tx, ty, tz, grad_sdf),
                   "Failed python call");
    return FTensor::Tensor1<double, 3>{grad_sdf[0], grad_sdf[1], grad_sdf[2]};
  }
#endif
  return FTensor::Tensor1<double, 3>{0., 1., 0.};
}

hess_surface_distance_function()

FTensor::Tensor2_symmetric< double, 3 > ContactOps::hess_surface_distance_function ( double  t, double  x, double  y, double  z, double  tx, double  ty, double  tz  )

inline

Examples

ContactOps.hpp.

Definition at line 243 of file ContactOps.hpp.

                                                                {
#ifdef PYTHON_SFD
  if (auto sdf_ptr = sdfPythonWeakPtr.lock()) {
    std::vector<double> hess_sdf;
    CHK_MOAB_THROW(sdf_ptr->evalHessSdf(t, x, y, z, tx, ty, tz, hess_sdf),
                   "Failed python call");
    return FTensor::Tensor2_symmetric<double, 3>{hess_sdf[0], hess_sdf[1],
                                                 hess_sdf[2], hess_sdf[3],
                                                 hess_sdf[4], hess_sdf[5]};
  }
#endif
  return FTensor::Tensor2_symmetric<double, 3>{0., 0., 0., 0., 0., 0.};
}

opFactoryBoundaryLhs()

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

MoFEMErrorCode ContactOps::opFactoryBoundaryLhs	(	boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	sigma,
		std::string	u
	)

Examples

ContactOps.hpp.

Definition at line 918 of file ContactOps.hpp.

                                  {
  MoFEMFunctionBegin;
 
  using C = ContactIntegrators<BoundaryEleOp>;
 
  auto common_data_ptr = boost::make_shared<ContactOps::CommonData>();
 
  pip.push_back(new OpCalculateVectorFieldValues<DIM>(
      u, common_data_ptr->contactDispPtr()));
  pip.push_back(new OpCalculateHVecTensorTrace<DIM, BoundaryEleOp>(
      sigma, common_data_ptr->contactTractionPtr()));
  pip.push_back(
      new typename C::template Assembly<A>::template OpConstrainBoundaryLhs_dU<
          DIM, GAUSS>(sigma, u, common_data_ptr));
  pip.push_back(new typename C::template Assembly<A>::
                    template OpConstrainBoundaryLhs_dTraction<DIM, GAUSS>(
                        sigma, sigma, common_data_ptr));
 
  MoFEMFunctionReturn(0);
}

opFactoryBoundaryRhs()

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

MoFEMErrorCode ContactOps::opFactoryBoundaryRhs	(	boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	sigma,
		std::string	u
	)

Examples

ContactOps.hpp.

Definition at line 942 of file ContactOps.hpp.

                                  {
  MoFEMFunctionBegin;
 
  using C = ContactIntegrators<BoundaryEleOp>;
 
  auto common_data_ptr = boost::make_shared<ContactOps::CommonData>();
 
  pip.push_back(new OpCalculateVectorFieldValues<DIM>(
      u, common_data_ptr->contactDispPtr()));
  pip.push_back(new OpCalculateHVecTensorTrace<DIM, BoundaryEleOp>(
      sigma, common_data_ptr->contactTractionPtr()));
  pip.push_back(
      new typename C::template Assembly<A>::template OpConstrainBoundaryRhs<
          DIM, GAUSS>(sigma, common_data_ptr));
 
  MoFEMFunctionReturn(0);
}

opFactoryBoundaryToDomainLhs()

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

MoFEMErrorCode ContactOps::opFactoryBoundaryToDomainLhs	(	MoFEM::Interface &	m_field,
		boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	fe_domain_name,
		std::string	sigma,
		std::string	u,
		std::string	geom,
		ForcesAndSourcesCore::RuleHookFun	rule
	)

Examples

ContactOps.hpp.

Definition at line 882 of file ContactOps.hpp.

                                                          {
  MoFEMFunctionBegin;
 
  using DomainEleOp = typename DomainEle::UserDataOperator;
 
  auto op_loop_side = new OpLoopSide<DomainEle>(
      m_field, fe_domain_name, DIM, Sev::noisy,
      boost::make_shared<ForcesAndSourcesCore::UserDataOperator::AdjCache>());
  pip.push_back(op_loop_side);
 
  CHKERR AddHOOps<DIM, DIM, DIM>::add(op_loop_side->getOpPtrVector(),
                                      {H1, HDIV}, geom);
 
  using B = typename FormsIntegrators<DomainEleOp>::template Assembly<
      A>::template BiLinearForm<I>;
 
  using OpMixDivULhs = typename B::template OpMixDivTimesVec<DIM>;
  using OpLambdaGraULhs = typename B::template OpMixTensorTimesGrad<DIM>;
  using OpMixDivULhsSide = OpMixLhsSide<OpMixDivULhs>;
  using OpLambdaGraULhsSide = OpMixLhsSide<OpLambdaGraULhs>;
 
  auto unity = []() { return 1; };
  op_loop_side->getOpPtrVector().push_back(
      new OpMixDivULhsSide(sigma, u, unity, true));
  op_loop_side->getOpPtrVector().push_back(
      new OpLambdaGraULhsSide(sigma, u, unity, true));
 
  op_loop_side->getSideFEPtr()->getRuleHook = rule;
  MoFEMFunctionReturn(0);
}

opFactoryCalculateTraction()

template<int DIM, IntegrationType I, typename BoundaryEleOp >

MoFEMErrorCode ContactOps::opFactoryCalculateTraction	(	boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	sigma
	)

Examples

ContactOps.hpp.

Definition at line 963 of file ContactOps.hpp.

                     {
  MoFEMFunctionBegin;
 
  using C = ContactIntegrators<BoundaryEleOp>;
 
  auto common_data_ptr = boost::make_shared<ContactOps::CommonData>();
  pip.push_back(new OpCalculateHVecTensorTrace<DIM, BoundaryEleOp>(
      sigma, common_data_ptr->contactTractionPtr()));
  pip.push_back(new typename C::template OpAssembleTotalContactTraction<DIM, I>(
      common_data_ptr, 1. / scale));
 
  MoFEMFunctionReturn(0);
}

opFactoryDomainRhs()

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

MoFEMErrorCode ContactOps::opFactoryDomainRhs	(	boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	sigma,
		std::string	u
	)

Examples

ContactOps.hpp.

Definition at line 826 of file ContactOps.hpp.

                                  {
  MoFEMFunctionBegin;
  using B = typename FormsIntegrators<DomainEleOp>::template Assembly<
      A>::template LinearForm<I>;
  using OpMixDivURhs = typename B::template OpMixDivTimesU<3, DIM, DIM>;
  using OpMixLambdaGradURhs = typename B::template OpMixTensorTimesGradU<DIM>;
  using OpMixUTimesDivLambdaRhs =
      typename B::template OpMixVecTimesDivLambda<SPACE_DIM>;
  using OpMixUTimesLambdaRhs =
      typename B::template OpGradTimesTensor<1, DIM, DIM>;
 
  auto common_data_ptr = boost::make_shared<ContactOps::CommonData>();
  auto mat_grad_ptr = boost::make_shared<MatrixDouble>();
  auto div_stress_ptr = boost::make_shared<MatrixDouble>();
  auto contact_stress_ptr = boost::make_shared<MatrixDouble>();
 
  pip.push_back(new OpCalculateVectorFieldValues<DIM>(
      u, common_data_ptr->contactDispPtr()));
  pip.push_back(
      new OpCalculateHVecTensorField<DIM, DIM>(sigma, contact_stress_ptr));
  pip.push_back(
      new OpCalculateHVecTensorDivergence<DIM, DIM>(sigma, div_stress_ptr));
 
  pip.push_back(new OpCalculateVectorFieldGradient<DIM, DIM>(u, mat_grad_ptr));
 
  pip.push_back(
      new OpMixDivURhs(sigma, common_data_ptr->contactDispPtr(),
                       [](double, double, double) constexpr { return 1; }));
  pip.push_back(new OpMixLambdaGradURhs(sigma, mat_grad_ptr));
  pip.push_back(new OpMixUTimesDivLambdaRhs(u, div_stress_ptr));
  pip.push_back(new OpMixUTimesLambdaRhs(u, contact_stress_ptr));
 
  MoFEMFunctionReturn(0);
}

sign()

double ContactOps::sign ( double  x )

inline

Examples

ContactOps.hpp.

Definition at line 378 of file ContactOps.hpp.

                             {
  if (x == 0)
    return 0;
  else if (x > 0)
    return 1;
  else
    return -1;
};

surface_distance_function()

double ContactOps::surface_distance_function ( double  t, double  x, double  y, double  z, double  tx, double  ty, double  tz  )

inline

Examples

ContactOps.hpp.

Definition at line 214 of file ContactOps.hpp.

                                                                         {
 
#ifdef PYTHON_SFD
  if (auto sdf_ptr = sdfPythonWeakPtr.lock()) {
    double sdf;
    CHK_MOAB_THROW(sdf_ptr->evalSdf(t, x, y, z, tx, ty, tz, sdf),
                   "Failed python call");
    return sdf;
  }
#endif
  return y + 0.5;
}

w()

double ContactOps::w ( const double  sdf, const double  tn  )

inline

Examples

ContactOps.hpp.

Definition at line 387 of file ContactOps.hpp.

                                                   {
  return sdf - cn_contact * tn;
}

Variable Documentation

cn_contact

double ContactOps::cn_contact = 1

Examples

ContactOps.hpp, and plastic.cpp.

Definition at line 19 of file EshelbianContact.hpp.

scale

double ContactOps::scale = 1

Examples

ContactOps.hpp.

Definition at line 20 of file EshelbianContact.hpp.