96c0acb825
* FEM: clean code clean code * [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci --------- Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
134 lines
5.1 KiB
C++
134 lines
5.1 KiB
C++
/***************************************************************************
|
|
* Copyright (c) 2013 Jürgen Riegel <FreeCAD@juergen-riegel.net> *
|
|
* *
|
|
* This file is part of the FreeCAD CAx development system. *
|
|
* *
|
|
* This library is free software; you can redistribute it and/or *
|
|
* modify it under the terms of the GNU Library General Public *
|
|
* License as published by the Free Software Foundation; either *
|
|
* version 2 of the License, or (at your option) any later version. *
|
|
* *
|
|
* This library is distributed in the hope that it will be useful, *
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
|
|
* GNU Library General Public License for more details. *
|
|
* *
|
|
* You should have received a copy of the GNU Library General Public *
|
|
* License along with this library; see the file COPYING.LIB. If not, *
|
|
* write to the Free Software Foundation, Inc., 59 Temple Place, *
|
|
* Suite 330, Boston, MA 02111-1307, USA *
|
|
* *
|
|
***************************************************************************/
|
|
|
|
#include "PreCompiled.h"
|
|
#include <SMESH_Version.h>
|
|
|
|
#ifndef _PreComp_
|
|
#include <Python.h>
|
|
#include <SMESHDS_Mesh.hxx>
|
|
#include <SMESH_Mesh.hxx>
|
|
|
|
#ifdef FCWithNetgen
|
|
#include <NETGENPlugin_Hypothesis.hxx>
|
|
#include <NETGENPlugin_Mesher.hxx>
|
|
#endif
|
|
#endif
|
|
|
|
#include <App/DocumentObjectPy.h>
|
|
#include <Base/Console.h>
|
|
#include <Mod/Part/App/PartFeature.h>
|
|
|
|
#include "FemMesh.h"
|
|
#include "FemMeshShapeNetgenObject.h"
|
|
|
|
|
|
using namespace Fem;
|
|
using namespace App;
|
|
|
|
PROPERTY_SOURCE(Fem::FemMeshShapeNetgenObject, Fem::FemMeshShapeBaseObject)
|
|
|
|
const char* FinenessEnums[] =
|
|
{"VeryCoarse", "Coarse", "Moderate", "Fine", "VeryFine", "UserDefined", nullptr};
|
|
|
|
FemMeshShapeNetgenObject::FemMeshShapeNetgenObject()
|
|
{
|
|
ADD_PROPERTY_TYPE(MaxSize, (1000), "MeshParams", Prop_None, "Maximum element size");
|
|
ADD_PROPERTY_TYPE(MinSize, (0), "MeshParams", Prop_None, "Minimum element size");
|
|
ADD_PROPERTY_TYPE(SecondOrder, (true), "MeshParams", Prop_None, "Create quadric elements");
|
|
ADD_PROPERTY_TYPE(Fineness, (2), "MeshParams", Prop_None, "Fineness level of the mesh");
|
|
Fineness.setEnums(FinenessEnums);
|
|
ADD_PROPERTY_TYPE(
|
|
GrowthRate,
|
|
(0.3),
|
|
"MeshParams",
|
|
Prop_None,
|
|
" allows defining how much the linear dimensions of two adjacent cells can differ");
|
|
ADD_PROPERTY_TYPE(
|
|
NbSegsPerEdge,
|
|
(1),
|
|
"MeshParams",
|
|
Prop_None,
|
|
"allows defining the minimum number of mesh segments in which edges will be split");
|
|
ADD_PROPERTY_TYPE(
|
|
NbSegsPerRadius,
|
|
(2),
|
|
"MeshParams",
|
|
Prop_None,
|
|
"allows defining the minimum number of mesh segments in which radii will be split");
|
|
ADD_PROPERTY_TYPE(Optimize, (true), "MeshParams", Prop_None, "Optimize the resulting mesh");
|
|
}
|
|
|
|
FemMeshShapeNetgenObject::~FemMeshShapeNetgenObject() = default;
|
|
|
|
App::DocumentObjectExecReturn* FemMeshShapeNetgenObject::execute()
|
|
{
|
|
#ifdef FCWithNetgen
|
|
|
|
Fem::FemMesh newMesh;
|
|
|
|
const Part::Feature* feat = Shape.getValue<Part::Feature*>();
|
|
if (!feat) {
|
|
return App::DocumentObject::StdReturn;
|
|
}
|
|
|
|
TopoDS_Shape shape = feat->Shape.getValue();
|
|
|
|
NETGENPlugin_Mesher myNetGenMesher(newMesh.getSMesh(), shape, true);
|
|
#if SMESH_VERSION_MAJOR >= 9
|
|
NETGENPlugin_Hypothesis* tet = new NETGENPlugin_Hypothesis(0, newMesh.getGenerator());
|
|
#else
|
|
NETGENPlugin_Hypothesis* tet = new NETGENPlugin_Hypothesis(0, 0, newMesh.getGenerator());
|
|
#endif
|
|
tet->SetMaxSize(MaxSize.getValue());
|
|
tet->SetMinSize(MinSize.getValue());
|
|
tet->SetSecondOrder(SecondOrder.getValue());
|
|
tet->SetOptimize(Optimize.getValue());
|
|
int iFineness = Fineness.getValue();
|
|
tet->SetFineness((NETGENPlugin_Hypothesis::Fineness)iFineness);
|
|
if (iFineness == 5) {
|
|
tet->SetGrowthRate(GrowthRate.getValue());
|
|
tet->SetNbSegPerEdge(NbSegsPerEdge.getValue());
|
|
tet->SetNbSegPerRadius(NbSegsPerRadius.getValue());
|
|
}
|
|
myNetGenMesher.SetParameters(tet);
|
|
newMesh.getSMesh()->ShapeToMesh(shape);
|
|
|
|
myNetGenMesher.Compute();
|
|
|
|
SMESHDS_Mesh* data = const_cast<SMESH_Mesh*>(newMesh.getSMesh())->GetMeshDS();
|
|
const SMDS_MeshInfo& info = data->GetMeshInfo();
|
|
int numFaces = data->NbFaces();
|
|
int numNode = info.NbNodes();
|
|
int numVolu = info.NbVolumes();
|
|
|
|
Base::Console().log("NetgenMesh: %i Nodes, %i Volumes, %i Faces\n", numNode, numVolu, numFaces);
|
|
|
|
FemMesh.setValue(newMesh);
|
|
return App::DocumentObject::StdReturn;
|
|
#else
|
|
return new App::DocumentObjectExecReturn(
|
|
"The FEM module is built without NETGEN support. Meshing will not work!!!",
|
|
this);
|
|
#endif
|
|
}
|