|
| v0.14.0
|
#include <users_modules/tutorials/adv-1/src/ContactOps.hpp>
template<int DIM, typename AssemblyBoundaryEleOp>
struct ContactOps::OpConstrainBoundaryRhsImpl< DIM, GAUSS, AssemblyBoundaryEleOp >
Definition at line 466 of file ContactOps.hpp.
◆ OpConstrainBoundaryRhsImpl()
template<int DIM, typename AssemblyBoundaryEleOp >
◆ iNtegrate()
template<int DIM, typename AssemblyBoundaryEleOp >
- Examples
- ContactOps.hpp.
Definition at line 726 of file ContactOps.hpp.
735 const size_t nb_gauss_pts = AssemblyBoundaryEleOp::getGaussPts().size2();
739 auto t_normal_at_pts = AssemblyBoundaryEleOp::getFTensor1NormalsAtGaussPts();
741 auto t_w = AssemblyBoundaryEleOp::getFTensor0IntegrationWeight();
742 auto t_disp = getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp);
743 auto t_traction = getFTensor1FromMat<DIM>(
commonDataPtr->contactTraction);
744 auto t_coords = AssemblyBoundaryEleOp::getFTensor1CoordsAtGaussPts();
746 size_t nb_base_functions = data.getN().size2() / 3;
747 auto t_base = data.getFTensor1N<3>();
751 getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp), nb_gauss_pts);
755 auto t_normal = getFTensor1FromMat<3>(m_normals_at_pts);
757 auto ts_time = AssemblyBoundaryEleOp::getTStime();
758 auto ts_time_step = AssemblyBoundaryEleOp::getTStimeStep();
765 m_spatial_coords, m_normals_at_pts, block_id);
769 m_spatial_coords, m_normals_at_pts, block_id);
772 auto t_grad_sdf = getFTensor1FromMat<3>(m_grad_sdf);
774 for (
size_t gg = 0; gg != nb_gauss_pts; ++gg) {
776 auto t_nf = getFTensor1FromPtr<DIM>(&nf[0]);
778 double jacobian = 1.;
780 jacobian = 2. * M_PI * t_coords(0);
782 const double alpha = t_w * jacobian * AssemblyBoundaryEleOp::getMeasure();
784 auto tn = -t_traction(
i) * t_grad_sdf(
i);
788 t_cP(
i,
j) = (
c * t_grad_sdf(
i)) * t_grad_sdf(
j);
799 t_cP(
i,
j) * t_disp(
j) +
800 c * (t_sdf * t_grad_sdf(
i));
804 const double beta = alpha * (t_base(
i) * t_normal(
i));
805 t_nf(
i) -= beta * t_rhs(
i);
810 for (; bb < nb_base_functions; ++bb)
◆ commonDataPtr
template<int DIM, typename AssemblyBoundaryEleOp >
◆ gradSurfaceDistanceFunction
template<int DIM, typename AssemblyBoundaryEleOp >
◆ isAxisymmetric
template<int DIM, typename AssemblyBoundaryEleOp >
◆ surfaceDistanceFunction
template<int DIM, typename AssemblyBoundaryEleOp >
The documentation for this struct was generated from the following file:
VectorDouble locF
local entity vector
const double c
speed of light (cm/ns)
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
FTensor::Index< 'i', SPACE_DIM > i
constexpr auto field_name
int nbRows
number of dofs on rows
Tensor2_Expr< Kronecker_Delta< T >, T, Dim0, Dim1, i, j > kronecker_delta(const Index< i, Dim0 > &, const Index< j, Dim1 > &)
Rank 2.
FTensor::Index< 'j', 3 > j
auto getFTensor1NormalsAtGaussPts()
get normal at integration points
FTensor::Index< 'k', 3 > k
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
FTensor::Index< 'l', 3 > l
auto getFTensor1CoordsAtGaussPts()
Get coordinates at integration points assuming linear geometry.