FractureMechanics::GriffithForceElement::MyTriangleFEConstrains Struct Reference

#include <users_modules/fracture_mechanics/src/GriffithForceElement.hpp>

Inheritance diagram for FractureMechanics::GriffithForceElement::MyTriangleFEConstrains:

Collaboration diagram for FractureMechanics::GriffithForceElement::MyTriangleFEConstrains:

Public Member Functions
	MyTriangleFEConstrains (MoFEM::Interface &m_field, const std::string &lambda_field_name, BlockData &block_data, SmartPetscObj< Vec > &delta_vec)
MoFEMErrorCode	preProcess ()
MoFEMErrorCode	postProcess ()
Public Member Functions inherited from FractureMechanics::GriffithForceElement::MyTriangleFE
	MyTriangleFE (MoFEM::Interface &m_field)
int	getRule (int order)

Public Attributes
const std::string	lambdaFieldName
BlockData &	blockData
SmartPetscObj< Vec >	deltaVec
Public Attributes inherited from FractureMechanics::GriffithForceElement::MyTriangleFE
Mat	B
Vec	F

Detailed Description

Definition at line 915 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

MyTriangleFEConstrains()

FractureMechanics::GriffithForceElement::MyTriangleFEConstrains::MyTriangleFEConstrains ( MoFEM::Interface &  m_field, const std::string &  lambda_field_name, BlockData &  block_data, SmartPetscObj< Vec > &  delta_vec  )

inline

Definition at line 921 of file GriffithForceElement.hpp.

        : MyTriangleFE(m_field), lambdaFieldName(lambda_field_name),
          blockData(block_data), deltaVec(delta_vec) {}

Member Function Documentation

postProcess()

MoFEMErrorCode FractureMechanics::GriffithForceElement::MyTriangleFEConstrains::postProcess ( )

inline

Definition at line 987 of file GriffithForceElement.hpp.

                                 {
      MoFEMFunctionBeginHot;
      MoFEMFunctionReturnHot(0);
    }

preProcess()

MoFEMErrorCode FractureMechanics::GriffithForceElement::MyTriangleFEConstrains::preProcess ( )

inline

Definition at line 927 of file GriffithForceElement.hpp.

                                {
      MoFEMFunctionBegin;
 
      MOFEM_LOG_CHANNEL("WORLD");
      MOFEM_LOG_FUNCTION();
      MOFEM_LOG_TAG("WORLD", "GriffithConstrain");
 
      const double *delta;
      CHKERR VecGetArrayRead(deltaVec, &delta);
      double sum_of_delta = 0;
      double sum_of_lambda = 0;
      const auto bit_field_number =
          mField.get_field_bit_number(lambdaFieldName);
      switch (snes_ctx) {
      case CTX_SNESSETFUNCTION: {
        for (_IT_NUMEREDDOF_ROW_BY_BITNUMBER_FOR_LOOP_(problemPtr,
                                                       bit_field_number, dit)) {
          if (static_cast<int>(dit->get()->getPart()) ==
              mField.get_comm_rank()) {
            int local_idx = dit->get()->getPetscLocalDofIdx();
            double lambda = dit->get()->getFieldData();
            double W = lambda - blockData.E * delta[local_idx];
            double val = lambda - calMax(W, 0, blockData.r);
            int global_idx = dit->get()->getPetscGlobalDofIdx();
            CHKERR VecSetValue(snes_f, global_idx, val, ADD_VALUES);
          }
        }
      } break;
      case CTX_SNESSETJACOBIAN: {
        for (_IT_NUMEREDDOF_ROW_BY_BITNUMBER_FOR_LOOP_(problemPtr,
                                                       bit_field_number, dit)) {
          if (static_cast<int>(dit->get()->getPart()) ==
              mField.get_comm_rank()) {
            int local_idx = dit->get()->getPetscLocalDofIdx();
            double lambda = dit->get()->getFieldData();
            double W = lambda - blockData.E * delta[local_idx];
            double diffW = 1;
            double val = 1 - diffCalMax_a(W, 0, blockData.r) * diffW;
            int global_idx = dit->get()->getPetscGlobalDofIdx();
            MOFEM_LOG_C("WORLD", Sev::verbose,
                        "Constrains on node %lu diag = %+3.5e "
                        "delta = %+3.5e lambda = %+3.5e",
                        dit->get()->getEnt(), val, delta[local_idx], lambda);
            CHKERR MatSetValue(snes_B, global_idx, global_idx, val, ADD_VALUES);
            sum_of_delta += delta[local_idx];
            sum_of_lambda += lambda;
          }
        }
        MOFEM_LOG_C("WORLD", Sev::noisy,
                    "Sum of delta = %+6.4e Sum of lambda = %+6.4e",
                    sum_of_delta, sum_of_lambda);
 
      } break;
      default:
        break;
      }
      CHKERR VecRestoreArrayRead(deltaVec, &delta);
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

blockData

BlockData& FractureMechanics::GriffithForceElement::MyTriangleFEConstrains::blockData

Definition at line 918 of file GriffithForceElement.hpp.

deltaVec

SmartPetscObj<Vec> FractureMechanics::GriffithForceElement::MyTriangleFEConstrains::deltaVec

Definition at line 919 of file GriffithForceElement.hpp.

lambdaFieldName

const std::string FractureMechanics::GriffithForceElement::MyTriangleFEConstrains::lambdaFieldName

Definition at line 917 of file GriffithForceElement.hpp.

---

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

• users_modules/fracture_mechanics/src/GriffithForceElement.hpp