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< VectorDouble(double delta_t, double t, int nb_gauss_pts, MatrixDouble &spatial_coords, MatrixDouble &normals_at_pts, int block_id)>
	[Common data] More...
using	GradSurfaceDistanceFunction = boost::function< MatrixDouble(double delta_t, double t, int nb_gauss_pts, MatrixDouble &spatial_coords, MatrixDouble &normals_at_pts, int block_id)>
using	HessSurfaceDistanceFunction = boost::function< MatrixDouble(double delta_t, double t, int nb_gauss_pts, MatrixDouble &spatial_coords, MatrixDouble &normals_at_pts, int block_id)>
using	PostProcEleDomain = PostProcEleByDim< SPACE_DIM >::PostProcEleDomain
using	SideEle = PostProcEleByDim< SPACE_DIM >::SideEle
using	PostProcEleBdy = PostProcEleByDim< SPACE_DIM >::PostProcEleBdy

Functions
VectorDouble	surface_distance_function (double delta_t, double t, int nb_gauss_pts, MatrixDouble &m_spatial_coords, MatrixDouble &m_normals_at_pts, int block_id)
MatrixDouble	grad_surface_distance_function (double delta_t, double t, int nb_gauss_pts, MatrixDouble &m_spatial_coords, MatrixDouble &m_normals_at_pts, int block_id)
MatrixDouble	hess_surface_distance_function (double delta_t, double t, int nb_gauss_pts, MatrixDouble &m_spatial_coords, MatrixDouble &m_normals_at_pts, int block_id)
template<typename T1 , typename T2 , int DIM1, int DIM2>
auto	get_spatial_coords (FTensor::Tensor1< T1, DIM1 > &&t_coords, FTensor::Tensor1< T2, DIM2 > &&t_disp, size_t nb_gauss_pts)
template<typename T1 , int DIM1>
auto	get_normalize_normals (FTensor::Tensor1< T1, DIM1 > &&t_normal_at_pts, size_t nb_gauss_pts)
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, bool is_axisymmetric=false)
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, bool is_axisymmetric=false)
template<int DIM, AssemblyType A, IntegrationType I, typename BoundaryEleOp >
MoFEMErrorCode	opFactoryBoundaryLhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma, std::string u, bool is_axisymmetric=false)
template<int DIM, AssemblyType A, IntegrationType I, typename BoundaryEleOp >
MoFEMErrorCode	opFactoryBoundaryRhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma, std::string u, bool is_axisymmetric=false)
template<int DIM, IntegrationType I, typename BoundaryEleOp >
MoFEMErrorCode	opFactoryCalculateTraction (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string sigma, bool is_axisymmetric=false)

Variables
double	cn_contact = 1
double	alpha_contact_const = 1
double	alpha_contact_quadratic = 1
double	scale = 1

Typedef Documentation

AssemblyBoundaryEleOp

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

Examples

ContactOps.hpp.

Definition at line 17 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

using ContactOps::GradSurfaceDistanceFunction = typedef boost::function<MatrixDouble( double delta_t, double t, int nb_gauss_pts, MatrixDouble &spatial_coords, MatrixDouble &normals_at_pts, int block_id)>

Examples

ContactOps.hpp.

Definition at line 210 of file ContactOps.hpp.

HessSurfaceDistanceFunction

using ContactOps::HessSurfaceDistanceFunction = typedef boost::function<MatrixDouble( double delta_t, double t, int nb_gauss_pts, MatrixDouble &spatial_coords, MatrixDouble &normals_at_pts, int block_id)>

Examples

ContactOps.hpp.

Definition at line 214 of file ContactOps.hpp.

PostProcEleBdy

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

Definition at line 26 of file PostProcContact.hpp.

PostProcEleDomain

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

Definition at line 24 of file PostProcContact.hpp.

SideEle

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

Definition at line 25 of file PostProcContact.hpp.

SurfaceDistanceFunction

using ContactOps::SurfaceDistanceFunction = typedef boost::function<VectorDouble( double delta_t, double t, int nb_gauss_pts, MatrixDouble &spatial_coords, MatrixDouble &normals_at_pts, int block_id)>

[Common data]

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

Examples

ContactOps.hpp.

Definition at line 206 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

Examples

ContactOps.hpp.

Definition at line 572 of file ContactOps.hpp.

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

get_normalize_normals()

template<typename T1 , int DIM1>

auto ContactOps::get_normalize_normals ( FTensor::Tensor1< T1, DIM1 > &&  t_normal_at_pts, size_t  nb_gauss_pts  )

inline

Examples

ContactOps.hpp.

Definition at line 422 of file ContactOps.hpp.

                                                       {
  MatrixDouble m_normals_at_pts(3, nb_gauss_pts);
  m_normals_at_pts.clear();
  FTensor::Index<'i', DIM1> i;
  auto t_set_normal = getFTensor1FromMat<3>(m_normals_at_pts);
  for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
    t_set_normal(i) = t_normal_at_pts(i) / t_normal_at_pts.l2();
    ++t_set_normal;
    ++t_normal_at_pts;
  }
  return m_normals_at_pts;
}

get_spatial_coords()

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

auto ContactOps::get_spatial_coords ( FTensor::Tensor1< T1, DIM1 > &&  t_coords, FTensor::Tensor1< T2, DIM2 > &&  t_disp, size_t  nb_gauss_pts  )

inline

Examples

ContactOps.hpp.

Definition at line 405 of file ContactOps.hpp.

                                                    {
  MatrixDouble m_spatial_coords(nb_gauss_pts, 3);
  m_spatial_coords.clear();
  auto t_spatial_coords = getFTensor1FromPtr<3>(&m_spatial_coords(0, 0));
  FTensor::Index<'i', DIM2> i;
  for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
    t_spatial_coords(i) = t_coords(i) + t_disp(i);
    ++t_spatial_coords;
    ++t_coords;
    ++t_disp;
  }
  return m_spatial_coords;
}

grad_surface_distance_function()

MatrixDouble ContactOps::grad_surface_distance_function ( double  delta_t, double  t, int  nb_gauss_pts, MatrixDouble &  m_spatial_coords, MatrixDouble &  m_normals_at_pts, int  block_id  )

inline

Examples

ContactOps.hpp.

Definition at line 273 of file ContactOps.hpp.

                                                                             {
#ifdef PYTHON_SDF
  if (auto sdf_ptr = sdfPythonWeakPtr.lock()) {
 
    VectorDouble v_spatial_coords = m_spatial_coords.data();
    VectorDouble v_normal_at_pts = m_normals_at_pts.data();
 
    bp::list python_coords;
    bp::list python_normals;
 
    for (int idx = 0; idx < 3; ++idx) {
      python_coords.append(
          convert_to_numpy(v_spatial_coords, nb_gauss_pts, idx));
      python_normals.append(
          convert_to_numpy(v_normal_at_pts, nb_gauss_pts, idx));
    }
 
    np::ndarray np_grad_sdf = np::empty(bp::make_tuple(nb_gauss_pts, 3),
                                        np::dtype::get_builtin<double>());
    CHK_MOAB_THROW(sdf_ptr->evalGradSdf(
                       delta_t, t, bp::extract<np::ndarray>(python_coords[0]),
                       bp::extract<np::ndarray>(python_coords[1]),
                       bp::extract<np::ndarray>(python_coords[2]),
                       bp::extract<np::ndarray>(python_normals[0]),
                       bp::extract<np::ndarray>(python_normals[1]),
                       bp::extract<np::ndarray>(python_normals[2]), block_id,
                       np_grad_sdf),
                   "Failed python call");
 
    double *grad_ptr = reinterpret_cast<double *>(np_grad_sdf.get_data());
 
    MatrixDouble m_grad_sdf;
    m_grad_sdf.resize(3, nb_gauss_pts, false);
    for (size_t gg = 0; gg < nb_gauss_pts; ++gg) {
      for (int idx = 0; idx < 3; ++idx)
        m_grad_sdf(idx, gg) = *(grad_ptr + (3 * gg + idx));
    }
    return m_grad_sdf;
  }
#endif
  MatrixDouble m_grad_sdf;
  m_grad_sdf.resize(3, nb_gauss_pts, false);
  FTensor::Index<'i', 3> i;
  FTensor::Tensor1<double, 3> t_grad_sdf_set{0.0, 0.0, 1.0};
  auto t_grad_sdf = getFTensor1FromMat<3>(m_grad_sdf);
 
  for (size_t gg = 0; gg < nb_gauss_pts; ++gg) {
    t_grad_sdf(i) = t_grad_sdf_set(i);
    ++t_grad_sdf;
  }
 
  return m_grad_sdf;
}

hess_surface_distance_function()

MatrixDouble ContactOps::hess_surface_distance_function ( double  delta_t, double  t, int  nb_gauss_pts, MatrixDouble &  m_spatial_coords, MatrixDouble &  m_normals_at_pts, int  block_id  )

inline

Examples

ContactOps.hpp.

Definition at line 330 of file ContactOps.hpp.

                                                                             {
#ifdef PYTHON_SDF
  if (auto sdf_ptr = sdfPythonWeakPtr.lock()) {
 
    VectorDouble v_spatial_coords = m_spatial_coords.data();
    VectorDouble v_normal_at_pts = m_normals_at_pts.data();
 
    bp::list python_coords;
    bp::list python_normals;
 
    for (int idx = 0; idx < 3; ++idx) {
      python_coords.append(
          convert_to_numpy(v_spatial_coords, nb_gauss_pts, idx));
      python_normals.append(
          convert_to_numpy(v_normal_at_pts, nb_gauss_pts, idx));
    };
 
    np::ndarray np_hess_sdf = np::empty(bp::make_tuple(nb_gauss_pts, 6),
                                        np::dtype::get_builtin<double>());
    CHK_MOAB_THROW(sdf_ptr->evalHessSdf(
                       delta_t, t, bp::extract<np::ndarray>(python_coords[0]),
                       bp::extract<np::ndarray>(python_coords[1]),
                       bp::extract<np::ndarray>(python_coords[2]),
                       bp::extract<np::ndarray>(python_normals[0]),
                       bp::extract<np::ndarray>(python_normals[1]),
                       bp::extract<np::ndarray>(python_normals[2]), block_id,
                       np_hess_sdf),
                   "Failed python call");
 
    double *hess_ptr = reinterpret_cast<double *>(np_hess_sdf.get_data());
 
    MatrixDouble m_hess_sdf;
    m_hess_sdf.resize(6, nb_gauss_pts, false);
    for (size_t gg = 0; gg < nb_gauss_pts; ++gg) {
      for (int idx = 0; idx < 6; ++idx)
        m_hess_sdf(idx, gg) =
            *(hess_ptr + (6 * gg + idx));
    }
    return m_hess_sdf;
  }
#endif
  MatrixDouble m_hess_sdf;
  m_hess_sdf.resize(6, nb_gauss_pts, false);
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Tensor2_symmetric<double, 3> t_hess_sdf_set{0., 0., 0., 0., 0., 0.};
  auto t_hess_sdf = getFTensor2SymmetricFromMat<3>(m_hess_sdf);
 
  for (size_t gg = 0; gg < nb_gauss_pts; ++gg) {
    t_hess_sdf(i, j) = t_hess_sdf_set(i, j);
    ++t_hess_sdf;
  }
  return m_hess_sdf;
}

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,
		bool	is_axisymmetric = false
	)

Examples

ContactOps.hpp.

Definition at line 1196 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, is_axisymmetric));
  pip.push_back(new typename C::template Assembly<A>::
                    template OpConstrainBoundaryLhs_dTraction<DIM, GAUSS>(
                        sigma, sigma, common_data_ptr, is_axisymmetric));
 
  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,
		bool	is_axisymmetric = false
	)

Examples

ContactOps.hpp.

Definition at line 1220 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, is_axisymmetric));
 
  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,
		bool	is_axisymmetric = false
	)

Examples

ContactOps.hpp.

Definition at line 1142 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 OpMixDivUCylLhs =
      typename B::template OpMixDivTimesVec<DIM, CYLINDRICAL>;
  using OpLambdaGraULhs = typename B::template OpMixTensorTimesGrad<DIM>;
 
  using OpMixDivULhsSide = OpMixLhsSide<OpMixDivULhs>;
  using OpMixDivUCylLhsSide = OpMixLhsSide<OpMixDivUCylLhs>;
  using OpLambdaGraULhsSide = OpMixLhsSide<OpLambdaGraULhs>;
 
  auto unity = []() { return 1; };
  auto jacobian = [is_axisymmetric](const double r, const double,
                                    const double) {
    if (is_axisymmetric)
      return 2. * M_PI * r;
    else
      return 1.;
  };
 
  if (!is_axisymmetric) {
    op_loop_side->getOpPtrVector().push_back(
        new OpMixDivULhsSide(sigma, u, unity, jacobian, true));
  } else {
    op_loop_side->getOpPtrVector().push_back(
        new OpMixDivUCylLhsSide(sigma, u, unity, jacobian, true));
  }
  op_loop_side->getOpPtrVector().push_back(
      new OpLambdaGraULhsSide(sigma, u, unity, jacobian, 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,
		bool	is_axisymmetric = false
	)

Examples

ContactOps.hpp.

Definition at line 1241 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, is_axisymmetric));
 
  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,
		bool	is_axisymmetric = false
	)

Examples

ContactOps.hpp.

Definition at line 1063 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 OpMixDivUCylRhs =
      typename B::template OpMixDivTimesU<3, DIM, DIM, CYLINDRICAL>;
 
  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>();
 
  auto jacobian = [is_axisymmetric](const double r, const double,
                                    const double) {
    if (is_axisymmetric)
      return 2. * M_PI * r;
    else
      return 1.;
  };
 
  pip.push_back(new OpCalculateVectorFieldValues<DIM>(
      u, common_data_ptr->contactDispPtr()));
  pip.push_back(
      new OpCalculateHVecTensorField<DIM, DIM>(sigma, contact_stress_ptr));
 
  if (!is_axisymmetric) {
    pip.push_back(
        new OpCalculateHVecTensorDivergence<DIM, DIM>(sigma, div_stress_ptr));
  } else {
    pip.push_back(new OpCalculateHVecTensorDivergence<DIM, DIM, CYLINDRICAL>(
        sigma, div_stress_ptr));
  }
 
  pip.push_back(new OpCalculateVectorFieldGradient<DIM, DIM>(u, mat_grad_ptr));
 
  if (!is_axisymmetric) {
    pip.push_back(
        new OpMixDivURhs(sigma, common_data_ptr->contactDispPtr(), jacobian));
  } else {
    pip.push_back(new OpMixDivUCylRhs(sigma, common_data_ptr->contactDispPtr(),
                                      jacobian));
  }
 
  pip.push_back(new OpMixLambdaGradURhs(sigma, mat_grad_ptr, jacobian));
  pip.push_back(new OpMixUTimesDivLambdaRhs(u, div_stress_ptr, jacobian));
  pip.push_back(new OpMixUTimesLambdaRhs(u, contact_stress_ptr, jacobian));
 
  MoFEMFunctionReturn(0);
}

sign()

double ContactOps::sign ( double  x )

inline

Examples

ContactOps.hpp, and PlasticOpsGeneric.hpp.

Definition at line 549 of file ContactOps.hpp.

                             {
  constexpr auto eps = std::numeric_limits<float>::epsilon();
  if (std::abs(x) < eps)
    return 0;
  else if (x > eps)
    return 1;
  else
    return -1;
};

surface_distance_function()

VectorDouble ContactOps::surface_distance_function ( double  delta_t, double  t, int  nb_gauss_pts, MatrixDouble &  m_spatial_coords, MatrixDouble &  m_normals_at_pts, int  block_id  )

inline

Examples

ContactOps.hpp.

Definition at line 216 of file ContactOps.hpp.

                                                            {
 
#ifdef PYTHON_SDF
  if (auto sdf_ptr = sdfPythonWeakPtr.lock()) {
 
    VectorDouble v_spatial_coords = m_spatial_coords.data();
    VectorDouble v_normal_at_pts = m_normals_at_pts.data();
 
    bp::list python_coords;
    bp::list python_normals;
 
    for (int idx = 0; idx < 3; ++idx) {
      python_coords.append(
          convert_to_numpy(v_spatial_coords, nb_gauss_pts, idx));
      python_normals.append(
          convert_to_numpy(v_normal_at_pts, nb_gauss_pts, idx));
    }
 
    np::ndarray np_sdf = np::empty(bp::make_tuple(nb_gauss_pts),
                                   np::dtype::get_builtin<double>());
    CHK_MOAB_THROW(sdf_ptr->evalSdf(delta_t, t,
                                    bp::extract<np::ndarray>(python_coords[0]),
                                    bp::extract<np::ndarray>(python_coords[1]),
                                    bp::extract<np::ndarray>(python_coords[2]),
                                    bp::extract<np::ndarray>(python_normals[0]),
                                    bp::extract<np::ndarray>(python_normals[1]),
                                    bp::extract<np::ndarray>(python_normals[2]),
                                    block_id, np_sdf),
                   "Failed python call");
 
    double *sdf_val_ptr = reinterpret_cast<double *>(np_sdf.get_data());
 
    VectorDouble v_sdf;
    v_sdf.resize(nb_gauss_pts, false);
 
    for (size_t gg = 0; gg < nb_gauss_pts; ++gg)
      v_sdf[gg] = *(sdf_val_ptr + gg);
 
    return v_sdf;
  }
#endif
  VectorDouble v_sdf;
  v_sdf.resize(nb_gauss_pts, false);
  m_spatial_coords.resize(3, nb_gauss_pts);
  auto t_coords = getFTensor1FromMat<3>(m_spatial_coords);
 
  for (size_t gg = 0; gg < nb_gauss_pts; ++gg)
    v_sdf[gg] = t_coords(1) + 0.5;
 
  return v_sdf;
}

w()

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

inline

Examples

ContactOps.hpp.

Definition at line 559 of file ContactOps.hpp.

                                                   {
  return sdf - cn_contact * tn;
}

Variable Documentation

alpha_contact_const

double ContactOps::alpha_contact_const = 1

Definition at line 20 of file EshelbianContact.hpp.

alpha_contact_quadratic

double ContactOps::alpha_contact_quadratic = 1

Definition at line 21 of file EshelbianContact.hpp.

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

adolc_plasticity.cpp, cell_forces.cpp, ContactOps.hpp, EshelbianOperators.cpp, phase.cpp, PlasticOps.hpp, test_jacobian_of_simple_contact_element.cpp, and UnsaturatedFlow.hpp.

Definition at line 22 of file EshelbianContact.hpp.