FractureMechanics::OpAnalyticalSpatialTraction Struct Reference

#include <users_modules/fracture_mechanics/src/CrackFrontElement.hpp>

Inheritance diagram for FractureMechanics::OpAnalyticalSpatialTraction:

Collaboration diagram for FractureMechanics::OpAnalyticalSpatialTraction:

Public Member Functions
	OpAnalyticalSpatialTraction (MoFEM::Interface &m_field, const bool &set_singular_coordinates, const Range &crack_front_nodes, const Range &crack_front_nodes_edges, const Range &crack_front_elements, PetscBool add_singularity, const std::string field_name, const Range tris, boost::ptr_vector< MethodForForceScaling > &methods_op, boost::ptr_vector< NeumannForcesSurface::MethodForAnalyticalForce > &analytical_force_op)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
const Range	tRis
boost::ptr_vector< MethodForForceScaling > &	methodsOp
boost::ptr_vector< NeumannForcesSurface::MethodForAnalyticalForce > &	analyticalForceOp
CrackFrontSingularBase< VolumeElementForcesAndSourcesCoreOnSide, VolumeElementForcesAndSourcesCore >	volSideFe
VectorDouble	sNrm
	Length of the normal vector. More...
bool	setSingularCoordinates
VectorDouble	Nf

Detailed Description

Definition at line 404 of file CrackFrontElement.hpp.

Constructor & Destructor Documentation

OpAnalyticalSpatialTraction()

FractureMechanics::OpAnalyticalSpatialTraction::OpAnalyticalSpatialTraction	(	MoFEM::Interface &	m_field,
		const bool &	set_singular_coordinates,
		const Range &	crack_front_nodes,
		const Range &	crack_front_nodes_edges,
		const Range &	crack_front_elements,
		PetscBool	add_singularity,
		const std::string	field_name,
		const Range	tris,
		boost::ptr_vector< MethodForForceScaling > &	methods_op,
		boost::ptr_vector< NeumannForcesSurface::MethodForAnalyticalForce > &	analytical_force_op
	)

Definition at line 881 of file CrackFrontElement.cpp.

    : FaceElementForcesAndSourcesCore::UserDataOperator(
          field_name, UserDataOperator::OPROW),
      tRis(tris), methodsOp(methods_op), analyticalForceOp(analytical_force_op),
      volSideFe(m_field, set_singular_coordinates, crack_front_nodes,
                crack_front_nodes_edges, crack_front_elements, add_singularity),
      setSingularCoordinates(set_singular_coordinates) {}

Member Function Documentation

doWork()

MoFEMErrorCode FractureMechanics::OpAnalyticalSpatialTraction::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

Definition at line 897 of file CrackFrontElement.cpp.

                                                                           {
  MoFEMFunctionBeginHot;
 
  if (data.getIndices().size() == 0)
    MoFEMFunctionReturnHot(0);
  EntityHandle ent = getNumeredEntFiniteElementPtr()->getEnt();
  if (tRis.find(ent) == tRis.end())
    MoFEMFunctionReturnHot(0);
 
  if (type == MBVERTEX) {
    CHKERR loopSideVolumes("ELASTIC", volSideFe);
  }
 
  int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
  int nb_row_dofs = data.getIndices().size() / rank;
 
  Nf.resize(data.getIndices().size(), false);
  Nf.clear();
 
  VectorDouble3 coords(3);
  VectorDouble3 normal(3);
  VectorDouble3 force(3);
  VectorDouble3 FTforce(3);
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
 
  //  If singular element apply transformation to h tensor
  if (type == MBVERTEX && volSideFe.singularElement) {
    double *inv_jac_ptr = &*volSideFe.invSJac.data().begin();
    FTensor::Tensor2<double *, 3, 3> t_inv_s_jac(
        inv_jac_ptr, &inv_jac_ptr[1], &inv_jac_ptr[2], &inv_jac_ptr[3],
        &inv_jac_ptr[4], &inv_jac_ptr[5], &inv_jac_ptr[6], &inv_jac_ptr[7],
        &inv_jac_ptr[8], 9);
    FTensor::Tensor1<double *, 3> t_normal = getFTensor1Normal();
    double nrm = 2 * getArea();
    FTensor::Tensor1<double, 3> t_s_normal;
    sNrm.resize(data.getN().size1(), false);
    for (unsigned int gg = 0; gg != data.getN().size1(); gg++) {
      // Push normal to tranformed singular element
      double det = volSideFe.detS[gg];
      t_s_normal(i) = det * t_inv_s_jac(j, i) * t_normal(j) / nrm;
      sNrm[gg] = sqrt(t_s_normal(i) * t_s_normal(i));
      ++t_inv_s_jac;
    }
  }
 
  FTensor::Tensor1<double *, 3> t_force(&force[0], &force[1], &force[2]);
 
  const auto &coords_at_gauss_pts = getCoordsAtGaussPts();
  const auto &normals_at_gauss_pts = getNormalsAtGaussPts();
 
  for (unsigned int gg = 0; gg < data.getN().size1(); gg++) {
 
    // get integration weight and Jacobian of integration point (area of face)
    double val = getGaussPts()(2, gg);
    val *= 0.5 * cblas_dnrm2(3, &getNormalsAtGaussPts()(gg, 0), 1);
    for (int dd = 0; dd != 3; dd++) {
      coords[dd] = coords_at_gauss_pts(gg, dd);
      normal[dd] = normals_at_gauss_pts(gg, dd);
    }
    if (volSideFe.singularElement) {
      for (int dd = 0; dd != 3; dd++) {
        coords[dd] += volSideFe.singularDisp(gg, dd);
      }
      val *= sNrm[gg];
    }
 
    for (boost::ptr_vector<
             NeumannForcesSurface::MethodForAnalyticalForce>::iterator vit =
             analyticalForceOp.begin();
         vit != analyticalForceOp.end(); vit++) {
 
      // Calculate force
      CHKERR vit->getForce(ent, coords, normal, force);
 
      for (int rr = 0; rr != 3; rr++) {
        cblas_daxpy(nb_row_dofs, val * force[rr], &data.getN()(gg, 0), 1,
                    &Nf[rr], rank);
      }
    }
  }
 
  // Scale force using user defined scaling operator
  CHKERR MethodForForceScaling::applyScale(getFEMethod(), methodsOp, Nf);
 
  // Assemble force into vector
  CHKERR VecSetValues(getFEMethod()->snes_f, data.getIndices().size(),
                      &data.getIndices()[0], &Nf[0], ADD_VALUES);
 
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

analyticalForceOp

boost::ptr_vector<NeumannForcesSurface::MethodForAnalyticalForce>& FractureMechanics::OpAnalyticalSpatialTraction::analyticalForceOp

Definition at line 410 of file CrackFrontElement.hpp.

methodsOp

boost::ptr_vector<MethodForForceScaling>& FractureMechanics::OpAnalyticalSpatialTraction::methodsOp

Definition at line 408 of file CrackFrontElement.hpp.

Nf

VectorDouble FractureMechanics::OpAnalyticalSpatialTraction::Nf

Definition at line 429 of file CrackFrontElement.hpp.

setSingularCoordinates

bool FractureMechanics::OpAnalyticalSpatialTraction::setSingularCoordinates

Definition at line 418 of file CrackFrontElement.hpp.

sNrm

VectorDouble FractureMechanics::OpAnalyticalSpatialTraction::sNrm

Length of the normal vector.

Definition at line 417 of file CrackFrontElement.hpp.

tRis

const Range FractureMechanics::OpAnalyticalSpatialTraction::tRis

Definition at line 407 of file CrackFrontElement.hpp.

volSideFe

CrackFrontSingularBase<VolumeElementForcesAndSourcesCoreOnSide, VolumeElementForcesAndSourcesCore> FractureMechanics::OpAnalyticalSpatialTraction::volSideFe

Definition at line 414 of file CrackFrontElement.hpp.

---

The documentation for this struct was generated from the following files:

• users_modules/fracture_mechanics/src/CrackFrontElement.hpp
• users_modules/fracture_mechanics/src/impl/CrackFrontElement.cpp