850                                    {
  852 
  857 
  858  const size_t nb_gauss_pts = AssemblyBoundaryEleOp::getGaussPts().size2();
  859 
  860  auto &nf = AssemblyBoundaryEleOp::locF;
  861 
  862  auto t_normal_at_pts = AssemblyBoundaryEleOp::getFTensor1NormalsAtGaussPts();
  863 
  864  auto t_w = AssemblyBoundaryEleOp::getFTensor0IntegrationWeight();
  865  auto t_disp = getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp);
 
  866  auto t_traction = getFTensor1FromMat<DIM>(
commonDataPtr->contactTraction);
 
  867  auto t_coords = AssemblyBoundaryEleOp::getFTensor1CoordsAtGaussPts();
  868 
  869  size_t nb_base_functions = data.getN().size2() / 3;
  870  auto t_base = data.getFTensor1N<3>();
  871 
  873      BoundaryEleOp::getFTensor1CoordsAtGaussPts(),
  874      getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp), nb_gauss_pts);
 
  876      BoundaryEleOp::getFTensor1NormalsAtGaussPts(), nb_gauss_pts);
  877 
  878  auto t_normal = getFTensor1FromMat<3>(m_normals_at_pts);
  879 
  880  auto ts_time = AssemblyBoundaryEleOp::getTStime();
  881  auto ts_time_step = AssemblyBoundaryEleOp::getTStimeStep();
  882 
  883  
  884  int block_id = 0;
  885 
  886  auto v_sdf =
  888                              m_spatial_coords, m_normals_at_pts, block_id);
  889 
  890  auto m_grad_sdf =
  892                                  m_spatial_coords, m_normals_at_pts, block_id);
  893 
  894  auto t_sdf = getFTensor0FromVec(v_sdf);
  895  auto t_grad_sdf = getFTensor1FromMat<3>(m_grad_sdf);
  896 
  897  for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
  898 
  899    auto t_nf = getFTensor1FromPtr<DIM>(&nf[0]);
  900 
  901    double jacobian = 1.;
  903      jacobian = 2. * M_PI * t_coords(0);
  904    }
  905    const double alpha = t_w * jacobian * AssemblyBoundaryEleOp::getMeasure();
  906 
  907    auto tn = -t_traction(
i) * t_grad_sdf(
i);
 
  909 
  911    t_cP(
i, 
j) = (
c * t_grad_sdf(
i)) * t_grad_sdf(
j);
 
  914 
  917 
  919 
  920        +
  921 
  922        t_cP(
i, 
j) * t_disp(
j) +
 
  923        c * (t_sdf * t_grad_sdf(
i)); 
 
  924 
  925    size_t bb = 0;
  926    for (; bb != AssemblyBoundaryEleOp::nbRows / DIM; ++bb) {
  927      const double beta = alpha * (t_base(
i) * t_normal(
i));
 
  928      t_nf(
i) -= beta * t_rhs(
i);
 
  929 
  930      ++t_nf;
  931      ++t_base;
  932    }
  933    for (; bb < nb_base_functions; ++bb)
  934      ++t_base;
  935 
  936    ++t_disp;
  937    ++t_traction;
  938    ++t_coords;
  939    ++t_w;
  940    ++t_normal;
  941    ++t_sdf;
  942    ++t_grad_sdf;
  943  }
  944 
  946}
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
const double c
speed of light (cm/ns)
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
Tensor2_Expr< Kronecker_Delta< T >, T, Dim0, Dim1, i, j > kronecker_delta(const Index< i, Dim0 > &, const Index< j, Dim1 > &)
Rank 2.