cubit_bc_test.cpp

Go to the documentation of this file.

/** \file cubit_bc_test.cpp
 * \example cubit_bc_test.cpp
 * \brief Atom test for getting boundary conditions from blocksets, sidesets and
 * nodesets.
 *
*/
 
#include <MoFEM.hpp>
 
using namespace MoFEM;
 
static char help[] = "Read file and print boundary conditions (ex. "
                     "./cubit_bc_test -my_file disp01.h5m) \n\n";
 
int main(int argc, char *argv[]) {
 
  MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
 
  try {
 
    moab::Core mb_instance;
    moab::Interface &moab = mb_instance;
    int rank;
    MPI_Comm_rank(PETSC_COMM_WORLD, &rank);
 
    // Read parameters from line command
    PetscBool flg = PETSC_TRUE;
    char mesh_file_name[255];
#if PETSC_VERSION_GE(3, 6, 4)
    CHKERR PetscOptionsGetString(PETSC_NULL, "", "-my_file", mesh_file_name,
                                 255, &flg);
#else
    CHKERR PetscOptionsGetString(PETSC_NULL, PETSC_NULL, "-my_file",
                                 mesh_file_name, 255, &flg);
#endif
    if (flg != PETSC_TRUE) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
              "*** ERROR -my_file (MESH FILE NEEDED)");
    }
 
    // Read mesh to MOAB
    const char *option;
    option = ""; //"PARALLEL=BCAST;";//;DEBUG_IO";
    CHKERR moab.load_file(mesh_file_name, 0, option);
 
    // Create MoFEM (Joseph) database
    MoFEM::Core core(moab);
    MoFEM::Interface &m_field = core;
 
    // Open mesh_file_name.txt for writing
    std::ofstream myfile;
    myfile.open((std::string(mesh_file_name) + ".txt").c_str());
 
    std::cout << "<<<< NODESETs >>>>>" << std::endl;
    // NODESETs
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, NODESET, it)) {
      std::cout << *it << std::endl;
      CHKERR it->printBcData(std::cout);
      std::vector<char> bc_data;
      CHKERR it->getBcData(bc_data);
      if (bc_data.empty())
        continue;
 
      // Displacement
      if (strcmp(&bc_data[0], "Displacement") == 0) {
        DisplacementCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      // Force
      else if (strcmp(&bc_data[0], "Force") == 0) {
        ForceCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      // Velocity
      else if (strcmp(&bc_data[0], "Velocity") == 0) {
        VelocityCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      // Acceleration
      else if (strcmp(&bc_data[0], "Acceleration") == 0) {
        AccelerationCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      // Temperature
      else if (strcmp(&bc_data[0], "Temperature") == 0) {
        TemperatureCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      else
        SETERRQ(PETSC_COMM_SELF, 1, "Error: Unrecognizable BC type");
    }
 
    std::cout << "<<<< SIDESETs >>>>>" << std::endl;
    // SIDESETs
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, SIDESET, it)) {
      std::cout << *it << std::endl;
      CHKERR it->printBcData(std::cout);
      std::vector<char> bc_data;
      CHKERR it->getBcData(bc_data);
      if (bc_data.empty())
        continue;
 
      // Pressure
      if (strcmp(&bc_data[0], "Pressure") == 0) {
        PressureCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      // Heat Flux
      else if (strcmp(&bc_data[0], "HeatFlux") == 0) {
        HeatFluxCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      }
 
      // cfd_bc
      else if (strcmp(&bc_data[0], "cfd_bc") == 0) {
        CfgCubitBcData mydata;
        CHKERR it->getBcDataStructure(mydata);
 
        // Interface bc (Hex:6 Dec:6)
        if (mydata.data.type == 6) { // 6 is the decimal value of the
                                     // corresponding value (hex) in bc_data
          // Print data
          std::cout << std::endl << "Interface" << std::endl;
          myfile << std::endl << "Interface" << std::endl;
          std::cout << mydata;
          myfile << mydata;
        }
        // Pressure inlet (Hex:f Dec:15)
        else if (mydata.data.type == 15) { // 15 is the decimal value of the
                                           // corresponding value (hex) in
                                           // bc_data
          // Print data
          std::cout << std::endl << "Pressure Inlet" << std::endl;
          myfile << std::endl << "Pressure Inlet" << std::endl;
          std::cout << mydata;
          myfile << mydata;
        }
        // Pressure outlet (Hex:10 Dec:16)
        else if (mydata.data.type == 16) { // 16 is the decimal value of the
                                           // corresponding value (hex) in
                                           // bc_data
          // Print data
          std::cout << std::endl << "Pressure Outlet" << std::endl;
          myfile << std::endl << "Pressure Outlet" << std::endl;
          std::cout << mydata;
          myfile << mydata;
        }
      }
 
      else
        SETERRQ(PETSC_COMM_SELF, 1, "Error: Unrecognizable BC type");
    }
 
    MeshsetsManager *meshsets_manager_ptr;
    CHKERR m_field.getInterface(meshsets_manager_ptr);
 
    std::cout << "<<<< BLOCKSETs >>>>>" << std::endl;
    // BLOCKSETs
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      std::cout << std::endl << *it << std::endl;
 
      // Get and print block name
      CHKERR it->printName(std::cout);
      CHKERR it->printName(myfile);
 
      // Get and print block attributes
      std::vector<double> attributes;
      CHKERR it->getAttributes(attributes);
      CHKERR it->printAttributes(std::cout);
      CHKERR it->printAttributes(myfile);
    }
 
    // Get block attributes and assign them as material properties/solution
    // parameters based on the name of each block
 
    // Conventions:
    //----------------------------------------------------------------------------------------
    // Materials are defined with block names starting with MAT_ e.g.
    // MAT_ELASTIC_abcd,
    // MAT_FRACTcdef etc.
    // Solution procedures are defined with block names starting with SOL_ e.g.
    // SOL_ELASTIC_xx, SOL_NLELASTICxx, SOL_FRACTabcd etc.
    //----------------------------------------------------------------------------------------
 
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
 
      std::cout << std::endl << *it << std::endl;
 
      // Get block name
      std::string name = it->getName();
 
      // Elastic material
      if (name.compare(0, 20, "MAT_ELASTIC_TRANSISO") == 0) {
        Mat_Elastic_TransIso mydata;
        CHKERR it->getAttributeDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      } else if (name.compare(0, 11, "MAT_ELASTIC") == 0) {
        Mat_Elastic mydata;
        CHKERR it->getAttributeDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      } else if (name.compare(0, 10, "MAT_INTERF") == 0) {
        Mat_Interf mydata;
        CHKERR it->getAttributeDataStructure(mydata);
        // Print data
        std::cout << mydata;
        myfile << mydata;
      } else
        SETERRQ(PETSC_COMM_SELF, 1, "Error: Unrecognizable Material type");
    }
 
    // Close mesh_file_name.txt
    myfile.close();
 
  } CATCH_ERRORS;
 
  return MoFEM::Core::Finalize();
}