FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db Struct Reference

assemble internal residual derivative More...

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

Inheritance diagram for FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db:

Collaboration diagram for FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db:

Public Member Functions
	OpAssemble_db (int tag, 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::FrontArcLengthControl::OpIntLambda
	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
int	tAg
AuxFunctions< adouble >	a_auxFun
VectorDouble	gradDelta
VectorDouble	activeVariables
Public Attributes inherited from FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpIntLambda
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 derivative

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 1396 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

OpAssemble_db()

FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db::OpAssemble_db ( int  tag, 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 1398 of file GriffithForceElement.hpp.

          : OpIntLambda(m_field, block_data, common_data, lambda_field_name,
                        arc_front),
            tAg(tag) {}

Member Function Documentation

doWork()

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

inline

Definition at line 1411 of file GriffithForceElement.hpp.

                                                                   {
        MoFEMFunctionBeginHot;
        if (type != MBVERTEX) {
          MoFEMFunctionReturnHot(0);
        }
        CHKERR setIndices(data);
        // get indices of Lagrange multiplier DOFs
        CHKERR setLambdaIndices(this, lambdaFieldName);
        double d_delta;
        adouble a_delta;
        // 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);
        a_auxFun.referenceCoords.resize(9, false);
        a_auxFun.currentCoords.resize(9, false);
 
        activeVariables.resize(18, false);
        trace_on(tAg);
        for (int dd = 0; dd != 9; dd++) {
          int loc_idx = data.getLocalIndices()[dd];
          if (loc_idx != -1) {
            a_auxFun.referenceCoords[dd] <<= x0[loc_idx];
            activeVariables[dd] = x0[loc_idx];
          } else {
            SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Note expected");
          }
        }
        for (int dd = 0; dd != 9; dd++) {
          int loc_idx = data.getLocalIndices()[dd];
          if (loc_idx != -1) {
            a_auxFun.currentCoords[dd] <<= x0[dd] + dx[loc_idx];
            activeVariables[9 + dd] = x0[loc_idx] + dx[loc_idx];
          } else {
            SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Note expected");
          }
        }
        // calculate crack area increment
 
        CHKERR a_auxFun.calculateGriffithForce(1, getGaussPts()(2, 0) * 0.5,
                                               data.getDiffN(0));
        // calculate area change
        a_delta = 0;
        for (int nn = 0; nn != 3; nn++) {
          if (rowLambdaIndices[nn] == -1)
            continue;
          auto griffith_force =
              getVectorAdaptor(&(a_auxFun.griffithForce[3 * nn]), 3);
          // griffith_force /= sqrt(griffith_force[0] * griffith_force[0] +
          //                         griffith_force[1] * griffith_force[1] +
          //                         griffith_force[2] * griffith_force[2]);
          for (int dd = 0; dd != 3; dd++) {
            const int idx = 3 * nn + dd;
            a_delta += griffith_force[dd] * (a_auxFun.currentCoords[idx] -
                                             a_auxFun.referenceCoords[idx]);
          }
        }
        a_delta >>= d_delta;
        trace_off();
        CHKERR VecRestoreArrayRead(arcFront.arcPtr->dx, &dx);
        CHKERR VecRestoreArrayRead(arcFront.arcPtr->x0, &x0);
 
        gradDelta.resize(18, false);
        {
          int r;
          // play recorder for jacobian
          r = ::gradient(tAg, 18, &activeVariables[0], &gradDelta[0]);
          if (r < 3) { // function is locally analytic
            SETERRQ1(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
                     "ADOL-C function evaluation with error r = %d", r);
          }
        }
 
        for (int dd = 0; dd != 9; dd++) {
          double val = arcFront.arcPtr->alpha * gradDelta[9 + dd];
          CHKERR VecSetValue(arcFront.arcPtr->db, data.getIndices()[dd], val,
                             ADD_VALUES);
        }
 
        MoFEMFunctionReturnHot(0);
      }

Member Data Documentation

a_auxFun

AuxFunctions<adouble> FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db::a_auxFun

Definition at line 1407 of file GriffithForceElement.hpp.

activeVariables

VectorDouble FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db::activeVariables

Definition at line 1409 of file GriffithForceElement.hpp.

gradDelta

VectorDouble FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db::gradDelta

Definition at line 1408 of file GriffithForceElement.hpp.

tAg

int FractureMechanics::GriffithForceElement::FrontArcLengthControl::OpAssemble_db::tAg

Definition at line 1397 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp