CohesiveElement::CohesiveInterfaceElement::OpCalculateGapGlobal Struct Reference

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

#include <users_modules/basic_finite_elements/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 342 of file CohesiveInterfaceElement.hpp.

Constructor & Destructor Documentation

OpCalculateGapGlobal()

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

inline

Definition at line 345 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 349 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,1,ss.str().c_str());
        }
      MoFEMFunctionReturnHot(0);
    }

Member Data Documentation

commonData

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

Definition at line 344 of file CohesiveInterfaceElement.hpp.

---

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

• users_modules/basic_finite_elements/cohesive_interface/src/CohesiveInterfaceElement.hpp