Tutorial for exchanging meshes with Kratos
Table of Contents
This tutorial shows how to exchange meshes with Kratos. It is required to do this tutorial first.
Mesh Exchange with Kratos
Here we try to send the mesh (in the form of CoSimIO::ModelPart
) to Kratos and get it back from it. Then we can check if both meshes are the same. Again the python file for Kratos side is very similar to the one described in the python tutorial. Importing the Kratos library is the first step:
import KratosMultiphysics as KM
from KratosMultiphysics.CoSimulationApplication import CoSimIO
Next is the creation of the Kratos::Model
(which is a container of Kratos::ModelPart
) and Kratos::ModelPart
. The ModelPart
will be used to import the mesh through CoSimIO
model = KM.Model()
model_part = model.CreateModelPart("mp_test")
Connecting is required before importing the mesh:
connection_settings = CoSimIO.Info()
connection_settings.SetString("my_name", "Kratos")
connection_settings.SetString("connect_to", "my_code")
connection_settings.SetInt("echo_level", 0)
info = CoSimIO.Connect(connection_settings)
connection_name = info.GetString("connection_name")
if info.GetInt("connection_status") != CoSimIO.ConnectionStatus.Connected:
raise Exception("Connecting failed")
Now everything is prepared to import the mesh, using the previously created Kratos::ModelPart
. After importing it is exported back.
import_info = CoSimIO.Info()
import_info.SetString("connection_name", connection_name)
import_info.SetString("identifier", "mesh_exchange_1")
CoSimIO.ImportMesh(import_info, model_part)
# print(model_part)
export_info = CoSimIO.Info()
export_info.SetString("connection_name", connection_name)
export_info.SetString("identifier", "mesh_exchange_2")
CoSimIO.ExportMesh(export_info, model_part)
Now as before the disconnection is done:
disconnect_settings = CoSimIO.Info()
disconnect_settings.SetString("connection_name", connection_name)
info = CoSimIO.Disconnect(disconnect_settings)
if info.GetInt("connection_status") != CoSimIO.ConnectionStatus.Disconnected:
raise Exception("Disconnecting failed")
The python file with the full script can be found here
On the other side we use first export mesh and then import it back, following what was done in this tutorial:
// System includes
#include <cmath> // std::abs
// CoSimulation includes
#include "co_sim_io.hpp"
#define COSIMIO_CHECK_EQUAL(a, b) \
if (a != b) { \
std::cout << "in line " << __LINE__ << " : " << a \
<< " is not equal to " << b << std::endl; \
return 1; \
}
int main()
{
CoSimIO::Info settings;
settings.Set("my_name", "my_code");
settings.Set("connect_to", "Kratos");
settings.Set("echo_level", 1);
auto info = CoSimIO::Connect(settings);
COSIMIO_CHECK_EQUAL(info.Get<int>("connection_status"), CoSimIO::ConnectionStatus::Connected);
const std::string connection_name = info.Get<std::string>("connection_name");
// send mesh to Kratos
CoSimIO::ModelPart model_part("my_mesh");
// fill the ModelPart with nodes & elements ...
for (int i=0; i<10; ++i) {
model_part.CreateNewNode(i+1, 1.1*i,0,0);
}
for (int i=0; i<5; ++i) {
model_part.CreateNewElement(i+1, CoSimIO::ElementType::Line2D2, {i+1, i+2});
}
info.Clear();
info.Set("identifier", "mesh_exchange_1");
info.Set("connection_name", connection_name);
info = CoSimIO::ExportMesh(info, model_part);
// receive the same data from Kratos
CoSimIO::ModelPart model_part_received("my_mesh_2");
info.Clear();
info.Set("identifier", "mesh_exchange_1");
info.Set("connection_name", connection_name);
info = CoSimIO::ImportMesh(info, model_part_received);
// check that the meshes are the same
// e.g. similar to "CoSimIO/tests/co_sim_io/cpp/co_sim_io_testing.hpp"
if (model_part.NumberOfNodes() != model_part_received.NumberOfNodes()) {
return 1;
}
if (model_part.NumberOfElements() != model_part_received.NumberOfElements()) {
return 1;
}
// do more detailed checks if necessary ...
// disconnecting afterwards
CoSimIO::Info disconnect_settings;
disconnect_settings.Set("connection_name", connection_name);
info = CoSimIO::Disconnect(disconnect_settings); // disconnect afterwards
COSIMIO_CHECK_EQUAL(info.Get<int>("connection_status"), CoSimIO::ConnectionStatus::Disconnected);
return 0;
}
This file can be found in integration_tutorials/cpp/mesh_exchange.cpp
Then you may run your executable with python script of Kratos from your working directory:
path/to/bin/mesh_exchange_cpp_test & python3 path/to/import_export_mesh.py
Next steps
In the next tutorial, meshes are exchanged with Kratos and mapping is used.