FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs Struct Reference

Calculate density at the crack fron nodes and multiplity gc. More...

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

Inheritance diagram for FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs:

Collaboration diagram for FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs:

Public Member Functions
	OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs (string rho_tag_name, const double beta_gc, boost::shared_ptr< MWLSApprox > mwls_approx, MoFEM::Interface &m_field, int tag, BlockData &block_data, CommonData &common_data)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &row_data)
Public Member Functions inherited from FractureMechanics::GriffithForceElement::AuxOp
	AuxOp (int tag, BlockData &block_data, CommonData &common_data)
MoFEMErrorCode	setIndices (DataForcesAndSourcesCore::EntData &data)
MoFEMErrorCode	setVariables (FaceElementForcesAndSourcesCore::UserDataOperator *fe_ptr, DataForcesAndSourcesCore::EntData &data)
MoFEMErrorCode	setLambdaNodes (FaceElementForcesAndSourcesCore::UserDataOperator *fe_ptr, const std::string &lambda_field_name)
MoFEMErrorCode	setLambdaIndices (FaceElementForcesAndSourcesCore::UserDataOperator *fe_ptr, const std::string &lambda_field_name)

Public Attributes
const double	betaGC
std::string	rhoTagName
boost::shared_ptr< MWLSApprox >	mwlsApprox
MoFEM::Interface &	mField
Public Attributes inherited from FractureMechanics::GriffithForceElement::AuxOp
int	tAg
BlockData &	blockData
CommonData &	commonData
ublas::vector< int >	rowIndices
ublas::vector< int >	rowIndicesLocal
VectorInt3	rowLambdaIndices
VectorInt3	rowLambdaIndicesLocal
VectorDouble3	lambdaAtNodes
VectorDouble	activeVariables

Detailed Description

Calculate density at the crack fron nodes and multiplity gc.

For heterogeneous case, gc is multiplied by a coefficient depending on spatial distribution of density as follows:

\[ g_c(\mathbf X) = g_c \cdot \rho(\mathbf X)^{\beta^{gc}} \]

Definition at line 432 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs()

FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs ( string  rho_tag_name, const double  beta_gc, boost::shared_ptr< MWLSApprox >  mwls_approx, MoFEM::Interface &  m_field, int  tag, BlockData &  block_data, CommonData &  common_data  )

inline

Definition at line 440 of file GriffithForceElement.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              "MESH_NODE_POSITIONS", UserDataOperator::OPROW),
          rhoTagName(rho_tag_name), betaGC(beta_gc), mwlsApprox(mwls_approx),
          mField(m_field), AuxOp(tag, block_data, common_data) {}

Member Function Documentation

doWork()

MoFEMErrorCode FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs::doWork ( int  side, EntityType  type, DataForcesAndSourcesCore::EntData &  row_data  )

inline

Definition at line 449 of file GriffithForceElement.hpp.

                                                                     {
      MoFEMFunctionBegin;
      FTensor::Index<'i', 3> i;
      FTensor::Index<'j', 3> j;
 
      if (type != MBVERTEX)
        MoFEMFunctionReturnHot(0);
      const int row_nb_dofs = row_data.getIndices().size();
      if (!row_nb_dofs)
        MoFEMFunctionReturnHot(0);
 
      CHKERR setIndices(row_data);
      CHKERR setVariables(this, row_data);
 
      auto &inv_AB_map = mwlsApprox->invABMap.at(getFEEntityHandle());
      auto &influence_nodes_map =
          mwlsApprox->influenceNodesMap.at(getFEEntityHandle());
      auto &dm_nodes_map = mwlsApprox->dmNodesMap[this->getFEEntityHandle()];
 
      auto get_ftensor_from_vec = [](auto &v) {
        return FTensor::Tensor1<double *, 3>(&v(0), &v(1), &v(2), 3);
      };
 
      Tag th;
      CHKERR mwlsApprox->mwlsMoab.tag_get_handle(rhoTagName.c_str(), th);
 
      int nb_dofs = row_data.getFieldData().size();
 
      commonData.diffRho.resize(3, nb_dofs / 3, false);
      std::fill(commonData.diffRho.data().begin(),
                commonData.diffRho.data().end(), 0);
 
      auto t_coords = get_ftensor_from_vec(row_data.getFieldData());
 
      commonData.densityRho.resize(nb_dofs / 3, false);
      std::fill(commonData.densityRho.begin(), commonData.densityRho.end(), 0);
      auto rho = getFTensor0FromVec(commonData.densityRho);
      VectorDouble gc_coeff(nb_dofs, false);
      auto t_coeff = get_ftensor_from_vec(gc_coeff);
      auto t_approx_base_point =
          getFTensor1FromMat<3>(mwlsApprox->approxBasePoint);
      const int nb_gauss_pts = row_data.getN().size1();
 
      for (int gg = 0; gg != nb_gauss_pts; ++gg) {
        for (int dd = 0; dd != nb_dofs / 3; ++dd) {
 
          // go only through nodes at the crack front
          if (rowIndices[dd * 3] != -1) {
            FTensor::Tensor1<double, 3> t_pos;
            t_pos(i) = t_coords(i);
 
            for (int dd : {0, 1, 2})
              mwlsApprox->approxPointCoords[dd] = t_approx_base_point(dd);
 
            mwlsApprox->getInfluenceNodes() = influence_nodes_map[gg];
            mwlsApprox->getInvAB() = inv_AB_map[gg];
            mwlsApprox->getShortenDm() = dm_nodes_map[gg];
 
            // approximate at MESH_NODE_POSITIONS
            CHKERR mwlsApprox->evaluateApproxFun(&t_pos(0));
            CHKERR mwlsApprox->evaluateFirstDiffApproxFun(&t_pos(0), false);
            CHKERR mwlsApprox->getTagData(th);
            rho = mwlsApprox->getDataApprox()[0];
            const auto &diff_vals = mwlsApprox->getDiffDataApprox();
            for (int ii = 0; ii != 3; ++ii) {
              commonData.diffRho(ii, dd) = diff_vals(ii, 0);
            }
            t_coeff(i) = pow(rho, betaGC);
          }
 
          ++rho;
          ++t_coords;
          ++t_coeff;
        }
      }
      // multiply griffith force with coefficients
      commonData.griffithForce =
          element_prod(commonData.griffithForce, gc_coeff);
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

betaGC

const double FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs::betaGC

Definition at line 435 of file GriffithForceElement.hpp.

mField

MoFEM::Interface& FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs::mField

Definition at line 438 of file GriffithForceElement.hpp.

mwlsApprox

boost::shared_ptr<MWLSApprox> FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs::mwlsApprox

Definition at line 437 of file GriffithForceElement.hpp.

rhoTagName

std::string FractureMechanics::GriffithForceElement::OpCalculateRhoAtCrackFrontUsingPrecalulatedCoeffs::rhoTagName

Definition at line 436 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp