Operator calculate element stiffens matrix. More...

#include "tutorials/cor-12_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 501 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 505 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 514 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);
        }
        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++) {
          CHKERR physicalEqations.calculateTangentStiffeness(D, gg, commonData,
                                                             getFEMethod());
          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];
                }
              }
            }
          }
        }
        CHKERR MatSetValues(getFEMethod()->snes_B, nb_row,
                            &row_data.getIndices()[0], nb_col,
                            &col_data.getIndices()[0], &K(0, 0), ADD_VALUES);
      } 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);
    }