CohesiveElement::ArcLengthElement Struct Reference

Inheritance diagram for CohesiveElement::ArcLengthElement:

Collaboration diagram for CohesiveElement::ArcLengthElement:

Public Member Functions
	ArcLengthElement (MoFEM::Interface &m_field, boost::shared_ptr< ArcLengthCtx > &arc_ptr)
MoFEMErrorCode	postProcessLoadPath ()
Public Member Functions inherited from CohesiveElement::ArcLengthIntElemFEMethod
	ArcLengthIntElemFEMethod (moab::Interface &moab, boost::shared_ptr< ArcLengthCtx > arcptr)
	~ArcLengthIntElemFEMethod ()
MoFEMErrorCode	remove_damaged_prisms_nodes ()
	remove nodes of prims which are fully damaged More...
MoFEMErrorCode	preProcess ()
MoFEMErrorCode	calculate_lambda_int (double &_lambda_int_)
virtual MoFEMErrorCode	calculate_db ()
MoFEMErrorCode	operator() ()
MoFEMErrorCode	postProcess ()
MoFEMErrorCode	calculate_dx_and_dlambda (Vec &x)
MoFEMErrorCode	calculate_init_dlambda (double *dlambda)
MoFEMErrorCode	set_dlambda_to_x (Vec &x, double dlambda)

Public Attributes
MoFEM::Interface &	mField
Range	postProcNodes
Public Attributes inherited from CohesiveElement::ArcLengthIntElemFEMethod
moab::Interface &	mOab
boost::shared_ptr< ArcLengthCtx >	arcPtr
Vec	GhostLambdaInt
Range	Faces3
Range	Faces4
Range	Edges3
Range	Edges4
Range	Nodes3
Range	Nodes4
Tag	thDamagedPrism
double	lambda_int

Detailed Description

Definition at line 31 of file arc_length_interface.cpp.

Constructor & Destructor Documentation

ArcLengthElement()

CohesiveElement::ArcLengthElement::ArcLengthElement ( MoFEM::Interface &  m_field, boost::shared_ptr< ArcLengthCtx > &  arc_ptr  )

inline

Definition at line 34 of file arc_length_interface.cpp.

      : ArcLengthIntElemFEMethod(m_field.get_moab(), arc_ptr), mField(m_field) {
 
    for (_IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_(mField, "LoadPath", cit)) {
      EntityHandle meshset = cit->getMeshset();
      Range nodes;
      rval = mOab.get_entities_by_type(meshset, MBVERTEX, nodes, true);
      MOAB_THROW(rval);
      postProcNodes.merge(nodes);
    }
 
    PetscPrintf(PETSC_COMM_WORLD, "Nb. PostProcNodes %lu\n",
                postProcNodes.size());
  };

Member Function Documentation

postProcessLoadPath()

MoFEMErrorCode CohesiveElement::ArcLengthElement::postProcessLoadPath ( )

inline

Definition at line 50 of file arc_length_interface.cpp.

                                       {
    MoFEMFunctionBegin;
    FILE *datafile;
    PetscFOpen(PETSC_COMM_SELF, DATAFILENAME, "a+", &datafile);
    const auto bit_number_lambda = mField.get_field_bit_number("LAMBDA");
    
    boost::shared_ptr<NumeredDofEntity_multiIndex> numered_dofs_rows =
        problemPtr->getNumeredRowDofsPtr();
    auto lit = numered_dofs_rows->lower_bound(
        FieldEntity::getLoBitNumberUId(bit_number_lambda));
    auto hi_lit = numered_dofs_rows->upper_bound(
        FieldEntity::getHiBitNumberUId(bit_number_lambda));
 
    if (lit == numered_dofs_rows->end()) {
      fclose(datafile);
      MoFEMFunctionReturnHot(0);
    } else if (std::distance(lit, hi_lit) != 1) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
              "Only one DOF is expected");
    }
 
    Range::iterator nit = postProcNodes.begin();
    for (; nit != postProcNodes.end(); nit++) {
      NumeredDofEntityByEnt::iterator dit, hi_dit;
      dit = numered_dofs_rows->get<Ent_mi_tag>().lower_bound(*nit);
      hi_dit = numered_dofs_rows->get<Ent_mi_tag>().upper_bound(*nit);
      double coords[3];
      CHKERR mOab.get_coords(&*nit, 1, coords);
      for (; dit != hi_dit; dit++) {
        PetscPrintf(PETSC_COMM_WORLD, "%s [ %d ] %6.4e -> ",
                    lit->get()->getName().c_str(), lit->get()->getDofCoeffIdx(),
                    lit->get()->getFieldData());
        PetscPrintf(PETSC_COMM_WORLD, "%s [ %d ] %6.4e ",
                    dit->get()->getName().c_str(), dit->get()->getDofCoeffIdx(),
                    dit->get()->getFieldData());
        PetscPrintf(PETSC_COMM_WORLD, "-> %3.4f %3.4f %3.4f\n", coords[0],
                    coords[1], coords[2]);
        PetscFPrintf(PETSC_COMM_WORLD, datafile, "%6.4e %6.4e ",
                     dit->get()->getFieldData(), lit->get()->getFieldData());
      }
    }
    PetscFPrintf(PETSC_COMM_WORLD, datafile, "\n");
    fclose(datafile);
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

mField

MoFEM::Interface& CohesiveElement::ArcLengthElement::mField

Definition at line 32 of file arc_length_interface.cpp.

postProcNodes

Range CohesiveElement::ArcLengthElement::postProcNodes

Definition at line 33 of file arc_length_interface.cpp.

---

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

• users_modules/basic_finite_elements/cohesive_interface/arc_length_interface.cpp