FractureMechanics::GriffithForceElement::OpCalculateGriffithForce Struct Reference

calculate griffith force vector More...

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

Inheritance diagram for FractureMechanics::GriffithForceElement::OpCalculateGriffithForce:

Collaboration diagram for FractureMechanics::GriffithForceElement::OpCalculateGriffithForce:

Public Member Functions
	OpCalculateGriffithForce (int tag, BlockData &block_data, CommonData &common_data)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &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
AuxFunctions< double >	auxFun
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 griffith force vector

\[ \mathbf{f}_{\mathrm{m}}^{\mathrm{h}}=\frac{1}{2}\left(\tilde{\mathbf{A}}_{\Gamma}^{\mathrm{h}}\right)^{\mathrm{T}} \mathbf{g}_{c}(\rho(\mathbf{X})) ]\f The matrix $\mathbf{A}_{\Gamma}^h$ comprising of direction vectors along the crack front that are normal to the crack front and tangent to the crack surface is defined as follows: \f[ A^h_{\Gamma} = \| \mathbf{N}(\tilde{\mathbf{X}}) \| = \left\| \epsilon_{ijk} \frac{\partial \Phi^\alpha_p}{\partial \xi_i} \frac{\partial \Phi^\beta_r}{\partial \xi_j} \tilde{X}^\alpha_p \tilde{X}^\beta_r \right\| \]

Definition at line 378 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

OpCalculateGriffithForce()

FractureMechanics::GriffithForceElement::OpCalculateGriffithForce::OpCalculateGriffithForce ( int  tag, BlockData &  block_data, CommonData &  common_data  )

inline

Definition at line 381 of file GriffithForceElement.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              "MESH_NODE_POSITIONS", UserDataOperator::OPROW),
          AuxOp(tag, block_data, common_data) {}

Member Function Documentation

doWork()

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

inline

Definition at line 387 of file GriffithForceElement.hpp.

                                                                 {
      MoFEMFunctionBegin;
 
      if (type != MBVERTEX)
        MoFEMFunctionReturnHot(0);
      CHKERR setIndices(data);
      CHKERR setVariables(this, data);
      int nb_dofs = data.getIndices().size();
      int nb_gauss_pts = data.getN().size1();
 
      auxFun.referenceCoords.resize(nb_dofs, false);
      auxFun.currentCoords.resize(nb_dofs, false);
 
      for (int dd = 0; dd != nb_dofs; dd++) {
        auxFun.referenceCoords[dd] = getCoords()[dd];
        auxFun.currentCoords[dd] = data.getFieldData()[dd];
      }
 
      for (int gg = 0; gg != nb_gauss_pts; gg++) {
        double val = getGaussPts()(2, gg) * 0.5;
        CHKERR auxFun.calculateGriffithForce(blockData.gc, val,
                                             data.getDiffN(gg));
      }
      commonData.griffithForce.resize(nb_dofs, false);
      std::fill(commonData.griffithForce.begin(),
                commonData.griffithForce.end(), 0);
      std::copy(auxFun.griffithForce.begin(), auxFun.griffithForce.end(),
                commonData.griffithForce.begin());
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

auxFun

AuxFunctions<double> FractureMechanics::GriffithForceElement::OpCalculateGriffithForce::auxFun

Definition at line 386 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp