FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda Struct Reference

assemble internal residual More...

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

Inheritance diagram for FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda:

Collaboration diagram for FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda:

Public Member Functions
	OpIntLambda (MoFEM::Interface &m_field, GriffithForceElement::BlockData &block_data, GriffithForceElement::CommonData &common_data, const std::string &lambda_field_name, FrontArcLengthControl &arc_front)
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
MoFEM::Interface &	mField
std::string	lambdaFieldName
FrontArcLengthControl &	arcFront
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

assemble internal residual

Parameters

block_data	griffith block data
common_data	griffith common data
lambda_field_name	lagrange multipliers controling crack surface area
arc_front	reference to FrontArcLengthControl object

Definition at line 1312 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

OpIntLambda()

FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda::OpIntLambda ( MoFEM::Interface &  m_field, GriffithForceElement::BlockData &  block_data, GriffithForceElement::CommonData &  common_data, const std::string &  lambda_field_name, FrontArcLengthControl &  arc_front  )

inline

Definition at line 1319 of file GriffithForceElement.hpp.

          : FaceElementForcesAndSourcesCore::UserDataOperator(
                "MESH_NODE_POSITIONS", lambda_field_name,
                UserDataOperator::OPROW,
                false // not symmetric operator
                ),
            AuxOp(0, block_data, common_data), mField(m_field),
            lambdaFieldName(lambda_field_name), arcFront(arc_front) {}

Member Function Documentation

doWork()

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

inline

Definition at line 1334 of file GriffithForceElement.hpp.

                                                                   {
        MoFEMFunctionBegin;
        if (type != MBVERTEX) {
          MoFEMFunctionReturnHot(0);
        }
        // get indices of Lagrange multiplier DOFs
        CHKERR setLambdaIndices(this, lambdaFieldName);
        // get dofs increment
        const double *dx;
        CHKERR VecGetArrayRead(arcFront.arcPtr->dx, &dx);
        // get dofs at begining of load step
        const double *x0;
        CHKERR VecGetArrayRead(arcFront.arcPtr->x0, &x0);
        auxFun.referenceCoords.resize(9, false);
        auxFun.currentCoords.resize(9, false);
        for (int dd = 0; dd != 9; dd++) {
          int loc_idx = data.getLocalIndices()[dd];
          if (loc_idx != -1) {
            auxFun.referenceCoords[dd] = x0[loc_idx];
            auxFun.currentCoords[dd] = auxFun.referenceCoords[dd] + dx[loc_idx];
          } else {
            SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Note expected");
          }
        }
        CHKERR VecRestoreArrayRead(arcFront.arcPtr->dx, &dx);
        CHKERR VecRestoreArrayRead(arcFront.arcPtr->x0, &x0);
 
        // calculate crack area increment
        CHKERR auxFun.calculateGriffithForce(1, getGaussPts()(2, 0) * 0.5,
                                             data.getDiffN(0));
 
        // calculate area change
        double delta = 0;
        for (int nn = 0; nn != 3; nn++) {
          if (rowLambdaIndices[nn] == -1)
            continue;
          auto griffith_force =
              getVectorAdaptor(&(auxFun.griffithForce[3 * nn]), 3);
          // griffith_force /= norm_2(griffith_force);
          for (int dd = 0; dd != 3; dd++) {
            const int idx = 3 * nn + dd;
            delta += griffith_force[dd] *
                     (auxFun.currentCoords[idx] - auxFun.referenceCoords[idx]);
          }
        }
        if (delta != delta) {
          SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY, "nan value");
        }
        arcFront.intLambdaDelta += delta;
        MoFEMFunctionReturn(0);
      };

Member Data Documentation

arcFront

FrontArcLengthControl& FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda::arcFront

Definition at line 1317 of file GriffithForceElement.hpp.

auxFun

AuxFunctions<double> FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda::auxFun

Definition at line 1332 of file GriffithForceElement.hpp.

lambdaFieldName

std::string FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda::lambdaFieldName

Definition at line 1316 of file GriffithForceElement.hpp.

mField

MoFEM::Interface& FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda::mField

Definition at line 1315 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp