MetaEdgeForces Struct Reference

#include <users_modules/basic_finite_elements/src/EdgeForce.hpp>

Static Public Member Functions
static MoFEMErrorCode	addElement (MoFEM::Interface &m_field, const std::string field_name, Range *intersect_ptr=NULL)
	Add element taking information from NODESET. More...
static MoFEMErrorCode	setOperators (MoFEM::Interface &m_field, boost::ptr_map< std::string, EdgeForce > &edge_forces, Vec F, const std::string field_name, std::string mesh_node_positions="MESH_NODE_POSITIONS")
	Set integration point operators. More...

Detailed Description

Definition at line 59 of file EdgeForce.hpp.

Member Function Documentation

addElement()

static MoFEMErrorCode MetaEdgeForces::addElement ( MoFEM::Interface &  m_field, const std::string  field_name, Range *  intersect_ptr = NULL  )

inlinestatic

Add element taking information from NODESET.

Examples

Remodeling.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, and nonlinear_dynamics.cpp.

Definition at line 62 of file EdgeForce.hpp.

                                                                {
    MoFEMFunctionBegin;
    CHKERR m_field.add_finite_element("FORCE_FE", MF_ZERO);
    CHKERR m_field.modify_finite_element_add_field_row("FORCE_FE", field_name);
    CHKERR m_field.modify_finite_element_add_field_col("FORCE_FE", field_name);
    CHKERR m_field.modify_finite_element_add_field_data("FORCE_FE", field_name);
    if (m_field.check_field("MESH_NODE_POSITIONS")) {
      CHKERR m_field.modify_finite_element_add_field_data(
          "FORCE_FE", "MESH_NODE_POSITIONS");
    }
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field, NODESET | FORCESET,
                                                    it)) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_type(it->meshset, MBTRI, tris,
                                                     true);
      Range edges;
      CHKERR m_field.get_moab().get_entities_by_type(it->meshset, MBEDGE, edges,
                                                     true);
      Range tris_edges;
      CHKERR m_field.get_moab().get_adjacencies(tris, 1, false, tris_edges,
                                                moab::Interface::UNION);
      edges = subtract(edges, tris_edges);
      if (intersect_ptr) {
        edges = intersect(edges, *intersect_ptr);
      }
      CHKERR m_field.add_ents_to_finite_element_by_type(edges, MBEDGE,
                                                        "FORCE_FE");
    }
    MoFEMFunctionReturn(0);
  }

setOperators()

static MoFEMErrorCode MetaEdgeForces::setOperators ( MoFEM::Interface &  m_field, boost::ptr_map< std::string, EdgeForce > &  edge_forces, Vec  F, const std::string  field_name, std::string  mesh_node_positions = "MESH_NODE_POSITIONS"  )

inlinestatic

Set integration point operators.

Examples

Remodeling.cpp, elasticity.cpp, and elasticity_mixed_formulation.cpp.

Definition at line 97 of file EdgeForce.hpp.

                                                                      {
    MoFEMFunctionBegin;
    string fe_name = "FORCE_FE";
    edge_forces.insert(fe_name, new EdgeForce(m_field));
    if (m_field.check_field(mesh_node_positions)) {
      auto &fe = edge_forces.at(fe_name).getLoopFe();
      fe.getOpPtrVector().push_back(
          new OpGetHOTangentsOnEdge(mesh_node_positions));
    }
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field, NODESET | FORCESET,
                                                    it)) {
      CHKERR edge_forces.at(fe_name).addForce(field_name, F,
                                              it->getMeshsetId());
    }
    MoFEMFunctionReturn(0);
  }

---

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

• users_modules/basic_finite_elements/src/EdgeForce.hpp