FractureMechanics::OpGetDensityFieldForTesting Struct Reference

Op to generate artificial density field. More...

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

Inheritance diagram for FractureMechanics::OpGetDensityFieldForTesting:

Collaboration diagram for FractureMechanics::OpGetDensityFieldForTesting:

Public Member Functions
	OpGetDensityFieldForTesting (const std::string row_field, boost::shared_ptr< MatrixDouble > mat_coords_ptr, boost::shared_ptr< VectorDouble > density_at_pts, boost::shared_ptr< MatrixDouble > rho_grad_at_gauss_pts_ptr, boost::shared_ptr< MatrixDouble > rho_grad_grad_at_gauss_pts_ptr, boost::shared_ptr< CrackFrontElement > fe_singular_ptr, MatrixDouble &singular_disp, boost::shared_ptr< MatrixDouble > mat_grad_pos_at_pts_ptr)
MoFEMErrorCode	doWork (int row_side, EntityType row_type, HookeElement::EntData &row_data)

Public Attributes
boost::shared_ptr< MatrixDouble >	matCoordsPtr
boost::shared_ptr< VectorDouble >	rhoAtGaussPtsPtr
boost::shared_ptr< MatrixDouble >	rhoGradAtGaussPtsPtr
boost::shared_ptr< MatrixDouble >	rhoGradGradAtGaussPtsPtr
boost::shared_ptr< CrackFrontElement >	feSingularPtr
MatrixDouble &	singularDisp
boost::shared_ptr< MatrixDouble >	matGradPosAtPtsPtr

Detailed Description

Op to generate artificial density field.

Definition at line 760 of file CrackFrontElement.hpp.

Constructor & Destructor Documentation

OpGetDensityFieldForTesting()

FractureMechanics::OpGetDensityFieldForTesting::OpGetDensityFieldForTesting ( const std::string  row_field, boost::shared_ptr< MatrixDouble >  mat_coords_ptr, boost::shared_ptr< VectorDouble >  density_at_pts, boost::shared_ptr< MatrixDouble >  rho_grad_at_gauss_pts_ptr, boost::shared_ptr< MatrixDouble >  rho_grad_grad_at_gauss_pts_ptr, boost::shared_ptr< CrackFrontElement >  fe_singular_ptr, MatrixDouble &  singular_disp, boost::shared_ptr< MatrixDouble >  mat_grad_pos_at_pts_ptr  )

inline

Definition at line 769 of file CrackFrontElement.hpp.

      : HookeElement::VolUserDataOperator(row_field, OPROW),
        matCoordsPtr(mat_coords_ptr), rhoAtGaussPtsPtr(density_at_pts),
        rhoGradAtGaussPtsPtr(rho_grad_at_gauss_pts_ptr),
        rhoGradGradAtGaussPtsPtr(rho_grad_grad_at_gauss_pts_ptr),
        feSingularPtr(fe_singular_ptr), singularDisp(singular_disp),
        matGradPosAtPtsPtr(mat_grad_pos_at_pts_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode FractureMechanics::OpGetDensityFieldForTesting::doWork	(	int	row_side,
		EntityType	row_type,
		HookeElement::EntData &	row_data
	)

Definition at line 1485 of file CrackFrontElement.cpp.

                                                                   {
  MoFEMFunctionBegin;
  if (row_type != MBVERTEX)
    MoFEMFunctionReturnHot(0);
  // get number of integration points
  const int nb_integration_pts = getGaussPts().size2();
  rhoAtGaussPtsPtr->resize(nb_integration_pts, false);
  rhoAtGaussPtsPtr->clear();
  rhoGradAtGaussPtsPtr->resize(3, nb_integration_pts, false);
  rhoGradAtGaussPtsPtr->clear();
  rhoGradGradAtGaussPtsPtr->resize(9, nb_integration_pts, false);
  rhoGradGradAtGaussPtsPtr->clear();
 
  auto t_rho = getFTensor0FromVec(*rhoAtGaussPtsPtr);
  auto t_grad_rho = getFTensor1FromMat<3>(*rhoGradAtGaussPtsPtr);
  auto t_grad_grad_rho = getFTensor2FromMat<3, 3>(*rhoGradGradAtGaussPtsPtr);
 
  MatrixDouble &coords = *matCoordsPtr;
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
 
  singularDisp.resize(3, nb_integration_pts, false);
  singularDisp.clear();
 
  auto t_material_positions = getFTensor1FromMat<3>(coords);
  auto t_mwls_material_positions = getFTensor1FromMat<3>(coords);
  if (feSingularPtr->singularElement) {
    MatrixDouble singular_current_displacement;
 
    singular_current_displacement.resize(3, nb_integration_pts, false);
    singular_current_displacement.clear();
 
    FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3> t_singular_displacement(
        &feSingularPtr->singularDisp(0, 0), &feSingularPtr->singularDisp(0, 1),
        &feSingularPtr->singularDisp(0, 2));
    auto t_inital_singular_displacement = getFTensor1FromMat<3>(singularDisp);
    auto t_current_singular_displacement =
        getFTensor1FromMat<3>(singular_current_displacement);
    if (!matGradPosAtPtsPtr)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Matrix for gradient of positions not allocated");
    auto t_H = getFTensor2FromMat<3, 3>(*matGradPosAtPtsPtr);
    for (int gg = 0; gg != nb_integration_pts; gg++) {
 
      t_inital_singular_displacement(i) = t_singular_displacement(i);
      t_current_singular_displacement(i) =
          t_H(i, j) * t_singular_displacement(j);
 
      t_material_positions(i) += t_current_singular_displacement(i);
 
      ++t_material_positions;
      ++t_H;
      ++t_singular_displacement;
      ++t_inital_singular_displacement;
      ++t_current_singular_displacement;
    }
  }
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
    t_rho =
        1 + t_mwls_material_positions(i) *
                t_mwls_material_positions(i); // density is equal to
                                              // 1+x^2+y^2+z^2 (x,y,z coordines)
    t_grad_rho(i) = 2 * t_mwls_material_positions(i);
    t_grad_grad_rho(0, 0) = 2.;
    t_grad_grad_rho(1, 1) = 2.;
    t_grad_grad_rho(2, 2) = 2.;
 
    ++t_rho;
    ++t_mwls_material_positions;
    ++t_grad_rho;
    ++t_grad_grad_rho;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

feSingularPtr

boost::shared_ptr<CrackFrontElement> FractureMechanics::OpGetDensityFieldForTesting::feSingularPtr

Definition at line 766 of file CrackFrontElement.hpp.

matCoordsPtr

boost::shared_ptr<MatrixDouble> FractureMechanics::OpGetDensityFieldForTesting::matCoordsPtr

Definition at line 762 of file CrackFrontElement.hpp.

matGradPosAtPtsPtr

boost::shared_ptr<MatrixDouble> FractureMechanics::OpGetDensityFieldForTesting::matGradPosAtPtsPtr

Definition at line 768 of file CrackFrontElement.hpp.

rhoAtGaussPtsPtr

boost::shared_ptr<VectorDouble> FractureMechanics::OpGetDensityFieldForTesting::rhoAtGaussPtsPtr

Definition at line 763 of file CrackFrontElement.hpp.

rhoGradAtGaussPtsPtr

boost::shared_ptr<MatrixDouble> FractureMechanics::OpGetDensityFieldForTesting::rhoGradAtGaussPtsPtr

Definition at line 764 of file CrackFrontElement.hpp.

rhoGradGradAtGaussPtsPtr

boost::shared_ptr<MatrixDouble> FractureMechanics::OpGetDensityFieldForTesting::rhoGradGradAtGaussPtsPtr

Definition at line 765 of file CrackFrontElement.hpp.

singularDisp

MatrixDouble& FractureMechanics::OpGetDensityFieldForTesting::singularDisp

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