FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs Struct Reference

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

Inheritance diagram for FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs:

Collaboration diagram for FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs:

Public Member Functions
	OpHeterogeneousGcLhs (int tag, BlockData &block_data, CommonData &common_data, const double beta_gc)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_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
MatrixDouble	mat
const double	betaGC
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

Definition at line 718 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

OpHeterogeneousGcLhs()

FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs::OpHeterogeneousGcLhs ( int  tag, BlockData &  block_data, CommonData &  common_data, const double  beta_gc  )

inline

Definition at line 722 of file GriffithForceElement.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              "MESH_NODE_POSITIONS", "MESH_NODE_POSITIONS",
              UserDataOperator::OPROWCOL),
          betaGC(beta_gc), AuxOp(tag, block_data, common_data) {}

Member Function Documentation

doWork()

MoFEMErrorCode FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

inline

Definition at line 732 of file GriffithForceElement.hpp.

                                                                     {
      MoFEMFunctionBegin;
 
      if (row_type != MBVERTEX) {
        MoFEMFunctionReturnHot(0);
      }
 
      FTensor::Index<'i', 3> i;
      FTensor::Index<'j', 3> j;
 
      AuxFunctions<double> auxFun;
      auxFun.currentCoords.resize(9, false);
 
      int nb_dofs = row_data.getFieldData().size();
      for (int dd = 0; dd != nb_dofs; dd++)
        auxFun.currentCoords[dd] = row_data.getFieldData()[dd];
 
      CHKERR setIndices(row_data);
      CHKERR setVariables(this, row_data);
      int nb_gauss_pts = row_data.getN().size1();
      int nb_base_row = row_data.getFieldData().size() / 3;
      int nb_base_col = col_data.getFieldData().size() / 3;
 
      int row_nb_dofs = row_data.getIndices().size();
      mat.resize(9, 9, false);
      mat.clear();
 
      const double &gc = blockData.gc;
      auto rho = getFTensor0FromVec(commonData.densityRho);
      auto get_ftensor_from_vec = [](auto &v) {
        return FTensor::Tensor1<double *, 3>(&v(0), &v(1), &v(2), 3);
      };
 
      auto get_tensor2 = [](MatrixDouble &m, const int r, const int c) {
        return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 3, 3>(
            &m(3 * r + 0, 3 * c + 0), &m(3 * r + 0, 3 * c + 1),
            &m(3 * r + 0, 3 * c + 2), &m(3 * r + 1, 3 * c + 0),
            &m(3 * r + 1, 3 * c + 1), &m(3 * r + 1, 3 * c + 2),
            &m(3 * r + 2, 3 * c + 0), &m(3 * r + 2, 3 * c + 1),
            &m(3 * r + 2, 3 * c + 2));
      };
 
      for (int gg = 0; gg != nb_gauss_pts; gg++) {
 
        double val = getGaussPts()(2, gg) * 0.5;
 
        auxFun.calculateGriffithForce(1, 1, row_data.getDiffN(gg));
 
        auto t_griffith = get_ftensor_from_vec(auxFun.griffithForce);
        auto t_diff_rho = getFTensor1FromMat<3>(commonData.diffRho);
        for (int rrr = 0; rrr != nb_base_row; rrr++) {
 
          const double a = val * gc * pow(rho, betaGC - 1.) * betaGC;
          FTensor::Tensor2<double, 3, 3> k;
          k(i, j) = t_diff_rho(j) * t_griffith(i) * a;
          // there is no dependence on the other nodes so row == col
          auto tt_mat = get_tensor2(mat, rrr, rrr);
          tt_mat(i, j) += k(i, j);
 
          ++t_diff_rho;
          ++t_griffith;
          ++rho;
        }
      }
 
      int col_nb_dofs = col_data.getIndices().size();
 
      CHKERR MatSetValues(getFEMethod()->snes_B, row_nb_dofs,
                          &*rowIndices.data().begin(), col_nb_dofs,
                          &*col_data.getIndices().data().begin(),
                          &*mat.data().begin(), ADD_VALUES);
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

betaGC

const double FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs::betaGC

Definition at line 730 of file GriffithForceElement.hpp.

mat

MatrixDouble FractureMechanics::GriffithForceElement::OpHeterogeneousGcLhs::mat

Definition at line 729 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp