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.