CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal Struct Reference

Operator calculate gap, normal vector and rotation matrix. More...

#include "tutorials/cor-12_cohesive_interface/src/CohesiveInterfaceElement.hpp"

Inheritance diagram for CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal:

Collaboration diagram for CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal:

Public Member Functions
	OpCalculateGapGlobal (const std::string field_name, CommonData &common_data)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)

Public Attributes
CommonData &	commonData

Detailed Description

Operator calculate gap, normal vector and rotation matrix.

Definition at line 351 of file CohesiveInterfaceElement.hpp.

Constructor & Destructor Documentation

OpCalculateGapGlobal()

CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal::OpCalculateGapGlobal ( const std::string  field_name, CommonData &  common_data  )

inline

Definition at line 355 of file CohesiveInterfaceElement.hpp.

        : FlatPrismElementForcesAndSourcesCore::UserDataOperator(
              field_name, ForcesAndSourcesCore::UserDataOperator::OPROW),
          commonData(common_data) {}

Member Function Documentation

doWork()

MoFEMErrorCode CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inline

Definition at line 360 of file CohesiveInterfaceElement.hpp.

                                                          {
      MoFEMFunctionBeginHot;
      try {
        int nb_dofs = data.getIndices().size();
        if (nb_dofs == 0)
          MoFEMFunctionReturnHot(0);
        int nb_gauss_pts = data.getN().size1();
        if (type == MBVERTEX) {
          commonData.R.resize(nb_gauss_pts);
          for (int gg = 0; gg < nb_gauss_pts; gg++) {
            commonData.R[gg].resize(3, 3);
            double nrm2_normal = 0;
            double nrm2_tangent1 = 0;
            double nrm2_tangent2 = 0;
            for (int dd = 0; dd < 3; dd++) {
              nrm2_normal += pow(getNormalsAtGaussPtsF3()(gg, dd), 2);
              nrm2_tangent1 += pow(getTangent1AtGaussPtF3()(gg, dd), 2);
              nrm2_tangent2 += pow(getTangent2AtGaussPtF3()(gg, dd), 2);
            }
            nrm2_normal = sqrt(nrm2_normal);
            nrm2_tangent1 = sqrt(nrm2_tangent1);
            nrm2_tangent2 = sqrt(nrm2_tangent2);
            for (int dd = 0; dd < 3; dd++) {
              commonData.R[gg](0, dd) =
                  getNormalsAtGaussPtsF3()(gg, dd) / nrm2_normal;
              commonData.R[gg](1, dd) =
                  getTangent1AtGaussPtF3()(gg, dd) / nrm2_tangent1;
              commonData.R[gg](2, dd) =
                  getTangent2AtGaussPtF3()(gg, dd) / nrm2_tangent2;
            }
          }
        }
        if (type == MBVERTEX) {
          commonData.gapGlob.resize(nb_gauss_pts, 3);
          commonData.gapGlob.clear();
        }
        for (int gg = 0; gg < nb_gauss_pts; gg++) {
          for (int dd = 0; dd < 3; dd++) {
            commonData.gapGlob(gg, dd) += cblas_ddot(
                nb_dofs / 3, &data.getN(gg)[0], 1, &data.getFieldData()[dd], 3);
          }
        }
      } catch (const std::exception &ex) {
        std::ostringstream ss;
        ss << "throw in method: " << ex.what() << std::endl;
        SETERRQ(PETSC_COMM_SELF, MOFEM_STD_EXCEPTION_THROW, "%s",
                ss.str().c_str());
      }
      MoFEMFunctionReturnHot(0);
    }

Member Data Documentation

commonData

CommonData& CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal::commonData

Definition at line 354 of file CohesiveInterfaceElement.hpp.

---

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

• tutorials/cor-12_cohesive_interface/src/CohesiveInterfaceElement.hpp