FractureMechanics::ConstantArea::OpAreaJacobian Struct Reference

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

Inheritance diagram for FractureMechanics::ConstantArea::OpAreaJacobian:

Collaboration diagram for FractureMechanics::ConstantArea::OpAreaJacobian:

Classes
struct	AuxFunctions

Public Member Functions
	OpAreaJacobian (int tag, CommonData &common_data)
PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
int	tAg
CommonData &	commonData
bool	isTapeRecorded
AuxFunctions< adouble >	a_auxFun

Detailed Description

Definition at line 240 of file ConstantArea.hpp.

Constructor & Destructor Documentation

OpAreaJacobian()

FractureMechanics::ConstantArea::OpAreaJacobian::OpAreaJacobian ( int  tag, CommonData &  common_data  )

inline

Definition at line 247 of file ConstantArea.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              "MESH_NODE_POSITIONS", UserDataOperator::OPROW),
          tAg(tag), commonData(common_data), isTapeRecorded(false) {
      sYmm = false;
    }

Member Function Documentation

doWork()

PetscErrorCode FractureMechanics::ConstantArea::OpAreaJacobian::doWork ( int  side, EntityType  type, DataForcesAndSourcesCore::EntData &  data  )

inline

Definition at line 423 of file ConstantArea.hpp.

                                                                 {
      MoFEMFunctionBegin;
 
      if (isTapeRecorded) {
        MoFEMFunctionReturnHot(0);
      }
      if (type != MBVERTEX) {
        MoFEMFunctionReturnHot(0);
      }
 
      int nb_dofs = data.getFieldData().size();
      int nb_gauss_pts = data.getN().size1();
 
      a_auxFun.referenceCoords.resize(9, false);
      a_auxFun.currentCoords.resize(9, false);
 
      trace_on(tAg);
 
      for (int dd = 0; dd != nb_dofs; dd++) {
        a_auxFun.referenceCoords[dd] <<= getCoords()[dd];
      }
      for (int dd = 0; dd != nb_dofs; dd++) {
        a_auxFun.currentCoords[dd] <<= data.getFieldData()[dd];
      }
 
      a_auxFun.crackAreaIncrement.resize(9, false);
      std::fill(a_auxFun.crackAreaIncrement.begin(),
                a_auxFun.crackAreaIncrement.end(), 0);
 
      for (int gg = 0; gg != nb_gauss_pts; ++gg) {
 
        double val = getGaussPts()(2, gg) * 0.5;
 
        CHKERR a_auxFun.matrixB(gg, data);
        CHKERR a_auxFun.dIffX();
        CHKERR a_auxFun.nOrmal();
        CHKERR a_auxFun.matrixA();
        CHKERR a_auxFun.calculateAreaIncrementDirection(val);
      }
 
      // cerr << "Griffith Force " << auxFun.griffithForce << endl;
 
      commonData.crackAreaIncrement.resize(nb_dofs, false);
      for (int dd = 0; dd != nb_dofs; dd++) {
        a_auxFun.crackAreaIncrement[dd] >>= commonData.crackAreaIncrement[dd];
      }
 
      trace_off();
 
      isTapeRecorded = true;
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

a_auxFun

AuxFunctions<adouble> FractureMechanics::ConstantArea::OpAreaJacobian::a_auxFun

Definition at line 421 of file ConstantArea.hpp.

commonData

CommonData& FractureMechanics::ConstantArea::OpAreaJacobian::commonData

Definition at line 244 of file ConstantArea.hpp.

isTapeRecorded

bool FractureMechanics::ConstantArea::OpAreaJacobian::isTapeRecorded

Definition at line 245 of file ConstantArea.hpp.

tAg

int FractureMechanics::ConstantArea::OpAreaJacobian::tAg

Definition at line 243 of file ConstantArea.hpp.

---

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

• users_modules/fracture_mechanics/src/ConstantArea.hpp