MetaNeumannForces Struct Reference

Set of high-level function declaring elements and setting operators to apply forces/fluxes. More...

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

Static Public Member Functions
static MoFEMErrorCode	addNeumannBCElements (MoFEM::Interface &m_field, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS", Range *intersect_ptr=NULL)
	Declare finite element. More...
static MoFEMErrorCode	setMomentumFluxOperators (MoFEM::Interface &m_field, boost::ptr_map< std::string, NeumannForcesSurface > &neumann_forces, Vec F, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")
	Set operators to finite elements calculating right hand side vector. More...
static MoFEMErrorCode	addNeumannFluxBCElements (MoFEM::Interface &m_field, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")
static MoFEMErrorCode	setMassFluxOperators (MoFEM::Interface &m_field, boost::ptr_map< std::string, NeumannForcesSurface > &neumann_forces, Vec F, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")

Detailed Description

Set of high-level function declaring elements and setting operators to apply forces/fluxes.

Definition at line 1261 of file SurfacePressure.hpp.

Member Function Documentation

addNeumannBCElements()

MoFEMErrorCode MetaNeumannForces::addNeumannBCElements ( MoFEM::Interface &  m_field, const std::string  field_name, const std::string  mesh_nodals_positions = "MESH_NODE_POSITIONS", Range *  intersect_ptr = NULL  )

static

Declare finite element.

Search cubit sidesets and blocksets with pressure bc and declare surface elemen

Block set has to have name “PRESSURE”. Can have name “PRESSURE_01” or any other name with prefix. The first attribute of block set is pressure value.

Parameters

m_field	Interface insurance
field_name	Field name (f.e. DISPLACEMENT)
mesh_nodals_positions	Name of field on which ho-geometry is defined
intersect_ptr	Pointer to range to interect meshset entities

Returns

Error code

Examples

elasticity.cpp, elasticity_mixed_formulation.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, nonlinear_dynamics.cpp, simple_contact.cpp, and simple_contact_thermal.cpp.

Definition at line 1974 of file SurfacePressure.cpp.

                                                                 {
  MoFEMFunctionBegin;
 
  // Define boundary element that operates on rows, columns and data of a
  // given field
  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_nodals_positions)) {
    CHKERR m_field.modify_finite_element_add_field_data("FORCE_FE",
                                                        mesh_nodals_positions);
  }
  // Add entities to that element, here we add all triangles with FORCESET
  // from cubit
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field, NODESET | FORCESET,
                                                  it)) {
    Range tris;
    CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                        true);
    if (intersect_ptr)
      tris = intersect(tris, *intersect_ptr);
    CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "FORCE_FE");
  }
 
  CHKERR m_field.add_finite_element("PRESSURE_FE", MF_ZERO);
  CHKERR m_field.modify_finite_element_add_field_row("PRESSURE_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_col("PRESSURE_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_data("PRESSURE_FE",
                                                      field_name);
  if (m_field.check_field(mesh_nodals_positions)) {
    CHKERR m_field.modify_finite_element_add_field_data("PRESSURE_FE",
                                                        mesh_nodals_positions);
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field,
                                                  SIDESET | PRESSURESET, it)) {
    Range tris;
    CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                        true);
    if (intersect_ptr)
      tris = intersect(tris, *intersect_ptr);
    CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "PRESSURE_FE");
  }
 
  // Reading forces from BLOCKSET
 
  const string block_set_force_name("FORCE");
  // search for block named FORCE and add its attributes to FORCE_FE element
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    if (it->getName().compare(0, block_set_force_name.length(),
                              block_set_force_name) == 0) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                          true);
      if (intersect_ptr)
        tris = intersect(tris, *intersect_ptr);
      CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "FORCE_FE");
    }
  }
  // search for block named PRESSURE and add its attributes to PRESSURE_FE
  // element
  const string block_set_pressure_name("PRESSURE");
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    if (it->getName().compare(0, block_set_pressure_name.length(),
                              block_set_pressure_name) == 0) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                          true);
      if (intersect_ptr)
        tris = intersect(tris, *intersect_ptr);
      CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "PRESSURE_FE");
    }
  }
 
  // search for block named LINEAR_PRESSURE and add its attributes to
  // PRESSURE_FE element
  const string block_set_linear_pressure_name("LINEAR_PRESSURE");
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    if (it->getName().compare(0, block_set_linear_pressure_name.length(),
                              block_set_linear_pressure_name) == 0) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                          true);
      if (intersect_ptr)
        tris = intersect(tris, *intersect_ptr);
      CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "PRESSURE_FE");
    }
  }
 
  MoFEMFunctionReturn(0);
}

addNeumannFluxBCElements()

MoFEMErrorCode MetaNeumannForces::addNeumannFluxBCElements ( MoFEM::Interface &  m_field, const std::string  field_name, const std::string  mesh_nodals_positions = "MESH_NODE_POSITIONS"  )

static

Definition at line 2141 of file SurfacePressure.cpp.

                                           {
  MoFEMFunctionBegin;
 
  CHKERR m_field.add_finite_element("FLUX_FE", MF_ZERO);
  CHKERR m_field.modify_finite_element_add_field_row("FLUX_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_col("FLUX_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_data("FLUX_FE", field_name);
  if (m_field.check_field(mesh_nodals_positions)) {
    CHKERR m_field.modify_finite_element_add_field_data("FLUX_FE",
                                                        mesh_nodals_positions);
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field,
                                                  SIDESET | PRESSURESET, it)) {
    Range tris;
    CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                        true);
    CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "FLUX_FE");
  }
 
  MoFEMFunctionReturn(0);
}

setMassFluxOperators()

MoFEMErrorCode MetaNeumannForces::setMassFluxOperators ( MoFEM::Interface &  m_field, boost::ptr_map< std::string, NeumannForcesSurface > &  neumann_forces, Vec  F, const std::string  field_name, const std::string  mesh_nodals_positions = "MESH_NODE_POSITIONS"  )

static

Definition at line 2166 of file SurfacePressure.cpp.

                                                                       {
  MoFEMFunctionBegin;
  bool ho_geometry = m_field.check_field(mesh_nodals_positions);
 
  string fe_name;
  fe_name = "FLUX_FE";
  neumann_forces.insert(fe_name, new NeumannForcesSurface(m_field));
 
  auto &nf = neumann_forces.at(fe_name);
  auto &fe = nf.getLoopFe();
 
  if (ho_geometry)
    CHKERR addHOOpsFace3D(mesh_nodals_positions, fe, false, false);
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field,
                                                  SIDESET | PRESSURESET, it)) {
    CHKERR nf.addFlux(field_name, F, it->getMeshsetId(), ho_geometry);
  }
  MoFEMFunctionReturn(0);
}

setMomentumFluxOperators()

MoFEMErrorCode MetaNeumannForces::setMomentumFluxOperators ( MoFEM::Interface &  m_field, boost::ptr_map< std::string, NeumannForcesSurface > &  neumann_forces, Vec  F, const std::string  field_name, const std::string  mesh_nodals_positions = "MESH_NODE_POSITIONS"  )

static

Set operators to finite elements calculating right hand side vector.

Parameters

m_field	Interface
neumann_forces	Map of pointers to force/pressure elements
F	Right hand side vector
field_name	Field name (f.e. DISPLACEMENT)
mesh_nodals_positions	Name of field on which ho-geometry is defined

Returns

Error code

Examples

elasticity.cpp, elasticity_mixed_formulation.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, simple_contact.cpp, and simple_contact_thermal.cpp.

Definition at line 2069 of file SurfacePressure.cpp.

                                                                       {
  MoFEMFunctionBegin;
  bool ho_geometry = m_field.check_field(mesh_nodals_positions);
  // Add forces
  {
    std::string fe_name = "FORCE_FE";
    neumann_forces.insert(fe_name, new NeumannForcesSurface(m_field));
    auto &nf = neumann_forces.at(fe_name);
    auto &fe = nf.getLoopFe();
 
    if (ho_geometry)
      CHKERR addHOOpsFace3D(mesh_nodals_positions, fe, false, false);
 
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field, NODESET | FORCESET,
                                                    it)) {
      CHKERR nf.addForce(field_name, F, it->getMeshsetId(), ho_geometry, false);
    }
    // Reading forces from BLOCKSET
    const string block_set_force_name("FORCE");
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, block_set_force_name.length(),
                                block_set_force_name) == 0) {
        CHKERR nf.addForce(field_name, F, it->getMeshsetId(), ho_geometry,
                           true);
      }
    }
  }
 
  // Add pressures
  {
    std::string fe_name = "PRESSURE_FE";
    neumann_forces.insert(fe_name, new NeumannForcesSurface(m_field));
 
    auto &nf = neumann_forces.at(fe_name);
    auto &fe = nf.getLoopFe();
 
    if (ho_geometry)
      CHKERR addHOOpsFace3D(mesh_nodals_positions, fe, false, false);
 
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
             m_field, SIDESET | PRESSURESET, it)) {
      CHKERR nf.addPressure(field_name, F, it->getMeshsetId(), ho_geometry,
                            false);
    }
 
    // Reading pressures from BLOCKSET
    const string block_set_pressure_name("PRESSURE");
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, block_set_pressure_name.length(),
                                block_set_pressure_name) == 0) {
        CHKERR nf.addPressure(field_name, F, it->getMeshsetId(), ho_geometry,
                              true);
      }
    }
 
    // Reading pressures from BLOCKSET
    const string block_set_linear_pressure_name("LINEAR_PRESSURE");
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, block_set_linear_pressure_name.length(),
                                block_set_linear_pressure_name) == 0) {
        CHKERR nf.addLinearPressure(field_name, F, it->getMeshsetId(),
                                    ho_geometry);
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

---

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

• users_modules/basic_finite_elements/src/SurfacePressure.hpp
• users_modules/basic_finite_elements/src/impl/SurfacePressure.cpp