kindred/create
1
1
Fork 0
Code Issues 53 Pull Requests Actions Packages Projects 9 Releases 2 Wiki Activity

Add ViewProvider and Dialogs for Analysis object and the Netgen mesher object

Browse Source
This commit is contained in:
jriegel
2013-04-06 23:08:23 +02:00
parent 720b4a2807
commit 2818bc4ddf
26 changed files with 1817 additions and 188 deletions
Download Patch File Download Diff File Expand all files Collapse all files

320
src/Mod/Fem/App/FemMeshShapeObject.cpp
View File

@@ -32,50 +32,6 @@
#include <Base/Placement.h>
#include <Mod/Part/App/PartFeature.h>
#include <SMESH_Gen.hxx>
#include <SMESH_Mesh.hxx>
#include <SMDS_PolyhedralVolumeOfNodes.hxx>
#include <SMDS_VolumeTool.hxx>
#include <StdMeshers_Arithmetic1D.hxx>
#include <StdMeshers_AutomaticLength.hxx>
#include <StdMeshers_MaxLength.hxx>
#include <StdMeshers_LocalLength.hxx>
#include <StdMeshers_MaxElementArea.hxx>
#include <StdMeshers_NotConformAllowed.hxx>
#include <StdMeshers_QuadranglePreference.hxx>
#include <StdMeshers_Quadrangle_2D.hxx>
#include <StdMeshers_Regular_1D.hxx>
#include <StdMeshers_UseExisting_1D2D.hxx>
#include <StdMeshers_CompositeSegment_1D.hxx>
#include <StdMeshers_Deflection1D.hxx>
#include <StdMeshers_Hexa_3D.hxx>
#include <StdMeshers_LayerDistribution.hxx>
#include <StdMeshers_LengthFromEdges.hxx>
#include <StdMeshers_MaxElementVolume.hxx>
#include <StdMeshers_MEFISTO_2D.hxx>
#include <StdMeshers_NumberOfLayers.hxx>
#include <StdMeshers_NumberOfSegments.hxx>
#include <StdMeshers_Prism_3D.hxx>
#include <StdMeshers_Projection_1D.hxx>
#include <StdMeshers_Projection_2D.hxx>
#include <StdMeshers_Projection_3D.hxx>
#include <StdMeshers_QuadraticMesh.hxx>
#include <StdMeshers_RadialPrism_3D.hxx>
#include <StdMeshers_SegmentAroundVertex_0D.hxx>
#include <StdMeshers_TrianglePreference.hxx>
#include <StdMeshers_ProjectionSource1D.hxx>
#include <StdMeshers_ProjectionSource2D.hxx>
#include <StdMeshers_ProjectionSource3D.hxx>
#include <StdMeshers_SegmentLengthAroundVertex.hxx>
#include <StdMeshers_StartEndLength.hxx>
//#include <StdMeshers_Propagation.hxx>
#include <StdMeshers_CompositeHexa_3D.hxx>
#include <NETGENPlugin_SimpleHypothesis_3D.hxx>
#include <NETGENPlugin_Mesher.hxx>
#include <BRepBuilderAPI_Copy.hxx>
#include <BRepTools.hxx>
using namespace Fem;
using namespace App;
@@ -91,140 +47,140 @@ FemMeshShapeObject::~FemMeshShapeObject()
{
}
App::DocumentObjectExecReturn *FemMeshShapeObject::execute(void)
{
Fem::FemMesh newMesh;
Part::Feature *feat = Shape.getValue<Part::Feature*>();
#if 0
TopoDS_Shape oshape = feat->Shape.getValue();
BRepBuilderAPI_Copy copy(oshape);
const TopoDS_Shape& shape = copy.Shape();
BRepTools::Clean(shape); // remove triangulation
#else
TopoDS_Shape shape = feat->Shape.getValue();
#endif
newMesh.getSMesh()->ShapeToMesh(shape);
SMESH_Gen *myGen = newMesh.getGenerator();
int hyp=0;
#if 0
SMESH_HypothesisPtr len(new StdMeshers_MaxLength(hyp++, 1, myGen));
static_cast<StdMeshers_MaxLength*>(len.get())->SetLength(1.0);
newMesh.addHypothesis(shape, len);
SMESH_HypothesisPtr loc(new StdMeshers_LocalLength(hyp++, 1, myGen));
static_cast<StdMeshers_LocalLength*>(loc.get())->SetLength(1.0);
newMesh.addHypothesis(shape, loc);
SMESH_HypothesisPtr area(new StdMeshers_MaxElementArea(hyp++, 1, myGen));
static_cast<StdMeshers_MaxElementArea*>(area.get())->SetMaxArea(1.0);
newMesh.addHypothesis(shape, area);
SMESH_HypothesisPtr segm(new StdMeshers_NumberOfSegments(hyp++, 1, myGen));
static_cast<StdMeshers_NumberOfSegments*>(segm.get())->SetNumberOfSegments(1);
newMesh.addHypothesis(shape, segm);
SMESH_HypothesisPtr defl(new StdMeshers_Deflection1D(hyp++, 1, myGen));
static_cast<StdMeshers_Deflection1D*>(defl.get())->SetDeflection(0.01);
newMesh.addHypothesis(shape, defl);
SMESH_HypothesisPtr reg(new StdMeshers_Regular_1D(hyp++, 1, myGen));
newMesh.addHypothesis(shape, reg);
//SMESH_HypothesisPtr sel(new StdMeshers_StartEndLength(hyp++, 1, myGen));
//static_cast<StdMeshers_StartEndLength*>(sel.get())->SetLength(1.0, true);
//newMesh.addHypothesis(shape, sel;
SMESH_HypothesisPtr qdp(new StdMeshers_QuadranglePreference(hyp++,1,myGen));
newMesh.addHypothesis(shape, qdp);
//SMESH_HypothesisPtr q2d(new StdMeshers_Quadrangle_2D(hyp++,1,myGen));
//newMesh.addHypothesis(shape, q2d);
SMESH_HypothesisPtr h3d(new StdMeshers_Hexa_3D(hyp++,1,myGen));
newMesh.addHypothesis(shape, h3d);
// create mesh
newMesh.compute();
#endif
#if 0 // Surface quad mesh
SMESH_HypothesisPtr len(new StdMeshers_MaxLength(hyp++, 1, myGen));
static_cast<StdMeshers_MaxLength*>(len.get())->SetLength(1.0);
newMesh.addHypothesis(shape, len);
SMESH_HypothesisPtr loc(new StdMeshers_LocalLength(hyp++, 1, myGen));
static_cast<StdMeshers_LocalLength*>(loc.get())->SetLength(1.0);
newMesh.addHypothesis(shape, loc);
SMESH_HypothesisPtr area(new StdMeshers_MaxElementArea(hyp++, 1, myGen));
static_cast<StdMeshers_MaxElementArea*>(area.get())->SetMaxArea(1.0);
newMesh.addHypothesis(shape, area);
SMESH_HypothesisPtr segm(new StdMeshers_NumberOfSegments(hyp++, 1, myGen));
static_cast<StdMeshers_NumberOfSegments*>(segm.get())->SetNumberOfSegments(1);
newMesh.addHypothesis(shape, segm);
SMESH_HypothesisPtr defl(new StdMeshers_Deflection1D(hyp++, 1, myGen));
static_cast<StdMeshers_Deflection1D*>(defl.get())->SetDeflection(0.01);
newMesh.addHypothesis(shape, defl);
SMESH_HypothesisPtr reg(new StdMeshers_Regular_1D(hyp++, 1, myGen));
newMesh.addHypothesis(shape, reg);
//SMESH_HypothesisPtr sel(new StdMeshers_StartEndLength(hyp++, 1, myGen));
//static_cast<StdMeshers_StartEndLength*>(sel.get())->SetLength(1.0, true);
//newMesh.addHypothesis(shape, sel;
SMESH_HypothesisPtr qdp(new StdMeshers_QuadranglePreference(hyp++,1,myGen));
newMesh.addHypothesis(shape, qdp);
SMESH_HypothesisPtr q2d(new StdMeshers_Quadrangle_2D(hyp++,1,myGen));
newMesh.addHypothesis(shape, q2d);
// create mesh
newMesh.compute();
#endif
#if 1 // NETGEN test
NETGENPlugin_Mesher myNetGenMesher(newMesh.getSMesh(),shape,true);
//NETGENPlugin_SimpleHypothesis_2D * tet2 = new NETGENPlugin_SimpleHypothesis_2D(hyp++,1,myGen);
//static_cast<NETGENPlugin_SimpleHypothesis_2D*>(tet2.get())->SetNumberOfSegments(5);
//static_cast<NETGENPlugin_SimpleHypothesis_2D*>(tet2.get())->SetLocalLength(0.1);
//static_cast<NETGENPlugin_SimpleHypothesis_2D*>(tet2.get())->LengthFromEdges();
//myNetGenMesher.SetParameters(tet2);
//NETGENPlugin_SimpleHypothesis_3D* tet= new NETGENPlugin_SimpleHypothesis_3D(hyp++,1,myGen);
//static_cast<NETGENPlugin_SimpleHypothesis_3D*>(tet.get())->LengthFromFaces();
//static_cast<NETGENPlugin_SimpleHypothesis_3D*>(tet.get())->SetMaxElementVolume(0.1);
//myNetGenMesher.SetParameters( tet);
myNetGenMesher.Compute();
#endif
SMESHDS_Mesh* data = const_cast<SMESH_Mesh*>(newMesh.getSMesh())->GetMeshDS();
const SMDS_MeshInfo& info = data->GetMeshInfo();
int numNode = info.NbNodes();
int numTria = info.NbTriangles();
int numQuad = info.NbQuadrangles();
int numPoly = info.NbPolygons();
int numVolu = info.NbVolumes();
int numTetr = info.NbTetras();
int numHexa = info.NbHexas();
int numPyrd = info.NbPyramids();
int numPris = info.NbPrisms();
int numHedr = info.NbPolyhedrons();
// set the value to the object
FemMesh.setValue(newMesh);
return App::DocumentObject::StdReturn;
}
//App::DocumentObjectExecReturn *FemMeshShapeObject::execute(void)
//{
// Fem::FemMesh newMesh;
//
// Part::Feature *feat = Shape.getValue<Part::Feature*>();
//
//#if 0
// TopoDS_Shape oshape = feat->Shape.getValue();
// BRepBuilderAPI_Copy copy(oshape);
// const TopoDS_Shape& shape = copy.Shape();
// BRepTools::Clean(shape); // remove triangulation
//#else
// TopoDS_Shape shape = feat->Shape.getValue();
//#endif
//
// newMesh.getSMesh()->ShapeToMesh(shape);
// SMESH_Gen *myGen = newMesh.getGenerator();
//
// int hyp=0;
//#if 0
// SMESH_HypothesisPtr len(new StdMeshers_MaxLength(hyp++, 1, myGen));
// static_cast<StdMeshers_MaxLength*>(len.get())->SetLength(1.0);
// newMesh.addHypothesis(shape, len);
//
// SMESH_HypothesisPtr loc(new StdMeshers_LocalLength(hyp++, 1, myGen));
// static_cast<StdMeshers_LocalLength*>(loc.get())->SetLength(1.0);
// newMesh.addHypothesis(shape, loc);
//
// SMESH_HypothesisPtr area(new StdMeshers_MaxElementArea(hyp++, 1, myGen));
// static_cast<StdMeshers_MaxElementArea*>(area.get())->SetMaxArea(1.0);
// newMesh.addHypothesis(shape, area);
//
// SMESH_HypothesisPtr segm(new StdMeshers_NumberOfSegments(hyp++, 1, myGen));
// static_cast<StdMeshers_NumberOfSegments*>(segm.get())->SetNumberOfSegments(1);
// newMesh.addHypothesis(shape, segm);
//
// SMESH_HypothesisPtr defl(new StdMeshers_Deflection1D(hyp++, 1, myGen));
// static_cast<StdMeshers_Deflection1D*>(defl.get())->SetDeflection(0.01);
// newMesh.addHypothesis(shape, defl);
//
// SMESH_HypothesisPtr reg(new StdMeshers_Regular_1D(hyp++, 1, myGen));
// newMesh.addHypothesis(shape, reg);
//
// //SMESH_HypothesisPtr sel(new StdMeshers_StartEndLength(hyp++, 1, myGen));
// //static_cast<StdMeshers_StartEndLength*>(sel.get())->SetLength(1.0, true);
// //newMesh.addHypothesis(shape, sel;
//
// SMESH_HypothesisPtr qdp(new StdMeshers_QuadranglePreference(hyp++,1,myGen));
// newMesh.addHypothesis(shape, qdp);
//
// //SMESH_HypothesisPtr q2d(new StdMeshers_Quadrangle_2D(hyp++,1,myGen));
// //newMesh.addHypothesis(shape, q2d);
//
// SMESH_HypothesisPtr h3d(new StdMeshers_Hexa_3D(hyp++,1,myGen));
// newMesh.addHypothesis(shape, h3d);
//
// // create mesh
// newMesh.compute();
//#endif
//#if 0 // Surface quad mesh
// SMESH_HypothesisPtr len(new StdMeshers_MaxLength(hyp++, 1, myGen));
// static_cast<StdMeshers_MaxLength*>(len.get())->SetLength(1.0);
// newMesh.addHypothesis(shape, len);
//
// SMESH_HypothesisPtr loc(new StdMeshers_LocalLength(hyp++, 1, myGen));
// static_cast<StdMeshers_LocalLength*>(loc.get())->SetLength(1.0);
// newMesh.addHypothesis(shape, loc);
//
// SMESH_HypothesisPtr area(new StdMeshers_MaxElementArea(hyp++, 1, myGen));
// static_cast<StdMeshers_MaxElementArea*>(area.get())->SetMaxArea(1.0);
// newMesh.addHypothesis(shape, area);
//
// SMESH_HypothesisPtr segm(new StdMeshers_NumberOfSegments(hyp++, 1, myGen));
// static_cast<StdMeshers_NumberOfSegments*>(segm.get())->SetNumberOfSegments(1);
// newMesh.addHypothesis(shape, segm);
//
// SMESH_HypothesisPtr defl(new StdMeshers_Deflection1D(hyp++, 1, myGen));
// static_cast<StdMeshers_Deflection1D*>(defl.get())->SetDeflection(0.01);
// newMesh.addHypothesis(shape, defl);
//
// SMESH_HypothesisPtr reg(new StdMeshers_Regular_1D(hyp++, 1, myGen));
// newMesh.addHypothesis(shape, reg);
//
// //SMESH_HypothesisPtr sel(new StdMeshers_StartEndLength(hyp++, 1, myGen));
// //static_cast<StdMeshers_StartEndLength*>(sel.get())->SetLength(1.0, true);
// //newMesh.addHypothesis(shape, sel;
//
// SMESH_HypothesisPtr qdp(new StdMeshers_QuadranglePreference(hyp++,1,myGen));
// newMesh.addHypothesis(shape, qdp);
//
// SMESH_HypothesisPtr q2d(new StdMeshers_Quadrangle_2D(hyp++,1,myGen));
// newMesh.addHypothesis(shape, q2d);
//
// // create mesh
// newMesh.compute();
//#endif
//#if 1 // NETGEN test
// NETGENPlugin_Mesher myNetGenMesher(newMesh.getSMesh(),shape,true);
//
// //NETGENPlugin_SimpleHypothesis_2D * tet2 = new NETGENPlugin_SimpleHypothesis_2D(hyp++,1,myGen);
// //static_cast<NETGENPlugin_SimpleHypothesis_2D*>(tet2.get())->SetNumberOfSegments(5);
// //static_cast<NETGENPlugin_SimpleHypothesis_2D*>(tet2.get())->SetLocalLength(0.1);
// //static_cast<NETGENPlugin_SimpleHypothesis_2D*>(tet2.get())->LengthFromEdges();
// //myNetGenMesher.SetParameters(tet2);
//
// //NETGENPlugin_SimpleHypothesis_3D* tet= new NETGENPlugin_SimpleHypothesis_3D(hyp++,1,myGen);
// //static_cast<NETGENPlugin_SimpleHypothesis_3D*>(tet.get())->LengthFromFaces();
// //static_cast<NETGENPlugin_SimpleHypothesis_3D*>(tet.get())->SetMaxElementVolume(0.1);
// //myNetGenMesher.SetParameters( tet);
//
// myNetGenMesher.Compute();
//#endif
//
//
//
// SMESHDS_Mesh* data = const_cast<SMESH_Mesh*>(newMesh.getSMesh())->GetMeshDS();
// const SMDS_MeshInfo& info = data->GetMeshInfo();
// int numNode = info.NbNodes();
// int numTria = info.NbTriangles();
// int numQuad = info.NbQuadrangles();
// int numPoly = info.NbPolygons();
// int numVolu = info.NbVolumes();
// int numTetr = info.NbTetras();
// int numHexa = info.NbHexas();
// int numPyrd = info.NbPyramids();
// int numPris = info.NbPrisms();
// int numHedr = info.NbPolyhedrons();
//
// // set the value to the object
// FemMesh.setValue(newMesh);
//
//
// return App::DocumentObject::StdReturn;
//}
//short FemMeshShapeObject::mustExecute(void) const
//{
@@ -240,7 +196,7 @@ App::DocumentObjectExecReturn *FemMeshShapeObject::execute(void)
// return Py::new_reference_to(PythonObject);
//}
void FemMeshShapeObject::onChanged(const Property* prop)
{
App::GeoFeature::onChanged(prop);
}
//void FemMeshShapeObject::onChanged(const Property* prop)
//{
// App::GeoFeature::onChanged(prop);
//}