FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags Struct Reference

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

Collaboration diagram for FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags:

Public Member Functions
	SaveGriffithForceOnTags (MoFEM::Interface &m_field, Tag th, Range &tris, Tag *th_coords=nullptr)
MoFEMErrorCode	operator() ()

Public Attributes
MoFEM::Interface &	mField
Tag	tH
Range &	tRis
Tag *	thCoordsPtr
AuxFunctions< double >	auxFun

Detailed Description

Definition at line 532 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

SaveGriffithForceOnTags()

FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags::SaveGriffithForceOnTags ( MoFEM::Interface &  m_field, Tag  th, Range &  tris, Tag *  th_coords = nullptr  )

inline

Definition at line 540 of file GriffithForceElement.hpp.

        : mField(m_field), tH(th), tRis(tris), thCoordsPtr(th_coords) {}

Member Function Documentation

operator()()

MoFEMErrorCode FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags::operator() ( )

inline

Definition at line 546 of file GriffithForceElement.hpp.

                                {
      MoFEMFunctionBegin;
 
      Range verts;
      CHKERR mField.get_moab().get_connectivity(tRis, verts, true);
      {
        std::vector<double> data(verts.size() * 3, 0);
        CHKERR mField.get_moab().tag_set_data(tH, verts, &*data.begin());
      }
      MatrixDouble diff_n(3, 2);
      CHKERR ShapeDiffMBTRI(&*diff_n.data().begin());
 
      for (Range::iterator tit = tRis.begin(); tit != tRis.end(); ++tit) {
        int num_nodes;
        const EntityHandle *conn;
        CHKERR mField.get_moab().get_connectivity(*tit, conn, num_nodes, true);
        VectorDouble9 coords(9);
 
        int nb_dofs = num_nodes * 3;
        auxFun.referenceCoords.resize(9, false);
        auxFun.currentCoords.resize(9, false);
 
        if (thCoordsPtr) {
          CHKERR mField.get_moab().tag_get_data(*thCoordsPtr, conn, num_nodes,
                                                &*coords.begin());
        } else {
          CHKERR mField.get_moab().get_coords(conn, num_nodes,
                                              &*coords.begin());
        }
        for (int dd = 0; dd != nb_dofs; dd++) {
          auxFun.currentCoords[dd] = coords[dd];
        }
        if (thCoordsPtr)
          CHKERR mField.get_moab().get_coords(conn, num_nodes,
                                              &*coords.begin());
        for (int dd = 0; dd != nb_dofs; dd++) {
          auxFun.referenceCoords[dd] = coords[dd];
        }
 
        CHKERR auxFun.calculateGriffithForce(1, 0.5, diff_n);
 
        double *tag_vals[3];
        CHKERR mField.get_moab().tag_get_by_ptr(tH, conn, 3,
                                                (const void **)tag_vals);
        for (int nn = 0; nn != 3; ++nn) {
          for (int dd = 0; dd != 3; ++dd) {
            tag_vals[nn][dd] += auxFun.griffithForce[3 * nn + dd];
          }
        }
      }
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

auxFun

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

Definition at line 544 of file GriffithForceElement.hpp.

mField

MoFEM::Interface& FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags::mField

Definition at line 534 of file GriffithForceElement.hpp.

tH

Tag FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags::tH

Definition at line 536 of file GriffithForceElement.hpp.

thCoordsPtr

Tag* FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags::thCoordsPtr

Definition at line 538 of file GriffithForceElement.hpp.

tRis

Range& FractureMechanics::GriffithForceElement::SaveGriffithForceOnTags::tRis

Definition at line 537 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp