Operator calculate element stiffens matrix. More...

#include "users_modules/basic_finite_elements/cohesive_interface/src/CohesiveInterfaceElement.hpp"

Inheritance diagram for CohesiveElement::CohesiveInterfaceElement::OpLhs:

Collaboration diagram for CohesiveElement::CohesiveInterfaceElement::OpLhs:

Public Member Functions
	OpLhs (const std::string field_name, CommonData &common_data, PhysicalEquation &physical_eqations)

MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Public Attributes
CommonData &	commonData

PhysicalEquation &	physicalEqations

MatrixDouble	K

MatrixDouble	D

MatrixDouble	ND

Detailed Description

Operator calculate element stiffens matrix.

Definition at line 475 of file CohesiveInterfaceElement.hpp.

Constructor & Destructor Documentation

◆ OpLhs()

CohesiveElement::CohesiveInterfaceElement::OpLhs::OpLhs	(	const std::string	field_name,
		CommonData &	common_data,
		PhysicalEquation &	physical_eqations )

inline

Definition at line 479 of file CohesiveInterfaceElement.hpp.

                                                                                                 :
    FlatPrismElementForcesAndSourcesCore::UserDataOperator(field_name,ForcesAndSourcesCore::UserDataOperator::OPROWCOL),
    commonData(common_data),physicalEqations(physical_eqations) { sYmm = false; }

Member Function Documentation

◆ doWork()

MoFEMErrorCode CohesiveElement::CohesiveInterfaceElement::OpLhs::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data )

inline

Definition at line 484 of file CohesiveInterfaceElement.hpp.

      {
      MoFEMFunctionBeginHot;
 
      try {
        int nb_row = row_data.getIndices().size();
        if(nb_row == 0) MoFEMFunctionReturnHot(0);
        int nb_col = col_data.getIndices().size();
        if(nb_col == 0) MoFEMFunctionReturnHot(0);
        if(physicalEqations.pRisms.find(getNumeredEntFiniteElementPtr()->getEnt())
        == physicalEqations.pRisms.end()) {
          MoFEMFunctionReturnHot(0);
        }
        //std::cerr << row_side << " " << row_type << " " << row_data.getN() << std::endl;
        //std::cerr << col_side << " " << col_type << " " << col_data.getN() << std::endl;
        ND.resize(nb_row,3);
        K.resize(nb_row,nb_col);
        K.clear();
        int nb_gauss_pts = row_data.getN().size1();
        for(int gg = 0;gg<nb_gauss_pts;gg++) {
          ierr = physicalEqations.calculateTangentStiffeness(D,gg,commonData,getFEMethod()); CHKERRG(ierr);
          double w = getGaussPts()(2,gg)*cblas_dnrm2(3,&getNormalsAtGaussPtsF3()(gg,0),1)*0.5;
          ND.clear();
          for(int nn = 0; nn<nb_row/3;nn++) {
            for(int dd = 0;dd<3;dd++) {
              for(int DD = 0;DD<3;DD++) {
                ND(3*nn+dd,DD) += row_data.getN(gg)[nn]*D(dd,DD);
              }
            }
          }
          for(int nn = 0; nn<nb_row/3; nn++) {
            for(int dd = 0;dd<3;dd++) {
              for(int NN = 0; NN<nb_col/3; NN++) {
                for(int DD = 0; DD<3;DD++) {
                  K(3*nn+dd,3*NN+DD) += w*ND(3*nn+dd,DD)*col_data.getN(gg)[NN];
                }
              }
            }
          }
        }
        ierr = MatSetValues(getFEMethod()->snes_B,
          nb_row,&row_data.getIndices()[0],
          nb_col,&col_data.getIndices()[0],
          &K(0,0),ADD_VALUES
        ); CHKERRG(ierr);
      } 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);
    }