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);
     }