diff --git a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp index aec517c09e33b7bc6768b2eff0964eddb99322e8..ec6f085e20f484b134a43fca21bb66453e7892e2 100644 --- a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp +++ b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp @@ -19,6 +19,7 @@ #include "MaterialLib/MPL/Medium.h" #include "MaterialLib/SolidModels/CreateConstitutiveRelation.h" #include "MaterialLib/SolidModels/MechanicsBase.h" +#include "MeshLib/Utils/Is2DMeshOnRotatedVerticalPlane.h" #include "ParameterLib/Utils.h" #include "ProcessLib/Output/CreateSecondaryVariables.h" #include "ProcessLib/Utils/ProcessUtils.h" @@ -42,6 +43,18 @@ std::unique_ptr<Process> createHydroMechanicsProcess( config.checkConfigParameter("type", "HYDRO_MECHANICS"); DBUG("Create HydroMechanicsProcess."); + if (DisplacementDim == 2) + { + if (mesh.isAxiallySymmetric() && + MeshLib::is2DMeshOnRotatedVerticalPlane(mesh)) + { + OGS_FATAL( + "Mesh {:s} is on a plane rotated around the vertical axis. The " + "axisymmetric problem can not use such mesh.", + mesh.getName()); + } + } + auto const staggered_scheme = //! \ogs_file_param{prj__processes__process__HYDRO_MECHANICS__coupling_scheme} config.getConfigParameterOptional<std::string>("coupling_scheme");