ForkMaster/geom
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/geom.git synced 2024-12-28 18:30:36 +05:00
Code Issues Packages Projects Releases Wiki Activity
4ab1782428
geom/src/GEOMImpl/GEOMImpl_IAdvancedOperations.cxx
vsr 4ab1782428 Merge from V6_main 15/03/2013
2013-03-15 14:28:11 +00:00

3024 lines
110 KiB
C++
Raw Blame History

// Copyright (C) 2007-2012 CEA/DEN, EDF R&D, OPEN CASCADE
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License.
//
// 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
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
//
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//
// File : GEOMImpl_IAdvancedOperations.cxx
// Author : Vadim SANDLER, Open CASCADE S.A.S. (vadim.sandler@opencascade.com)
#include "GEOMImpl_IAdvancedOperations.hxx"
#include <Basics_OCCTVersion.hxx>
#include <utilities.h>
#include <OpUtil.hxx>
#include <Utils_ExceptHandlers.hxx>
#include "GEOM_Function.hxx"
#include "GEOM_PythonDump.hxx"
#include "GEOMUtils.hxx"
#include "GEOMAlgo_Splitter.hxx"
#include "GEOMImpl_Gen.hxx"
#include "GEOMImpl_Types.hxx"
#include "GEOMImpl_IBasicOperations.hxx"
#include "GEOMImpl_IBooleanOperations.hxx"
#include "GEOMImpl_IShapesOperations.hxx"
#include "GEOMImpl_ITransformOperations.hxx"
#include "GEOMImpl_IBlocksOperations.hxx"
#include "GEOMImpl_I3DPrimOperations.hxx"
#include "GEOMImpl_ILocalOperations.hxx"
#include "GEOMImpl_IHealingOperations.hxx"
#include "GEOMImpl_GlueDriver.hxx"
#include "GEOMImpl_PipeTShapeDriver.hxx"
#include "GEOMImpl_IPipeTShape.hxx"
#include "GEOMImpl_DividedDiskDriver.hxx"
#include "GEOMImpl_IDividedDisk.hxx"
// #include "GEOMImpl_DividedCylinderDriver.hxx"
// #include "GEOMImpl_IDividedCylinder.hxx"
/*@@ insert new functions before this line @@ do not remove this line @@ do not remove this line @@*/
#include <TDF_Tool.hxx>
#include <TFunction_DriverTable.hxx>
#include <TFunction_Driver.hxx>
#include <TFunction_Logbook.hxx>
#include <TNaming_CopyShape.hxx>
#include <TopExp.hxx>
#include <TopExp_Explorer.hxx>
#include <TopoDS.hxx>
#include <TopoDS_Vertex.hxx>
#include <TopTools_IndexedMapOfShape.hxx>
#include <TopTools_ListIteratorOfListOfShape.hxx>
#include <TColStd_IndexedDataMapOfTransientTransient.hxx>
#include <BRep_Builder.hxx>
#include <BRep_Tool.hxx>
#include <BRepAlgoAPI_Cut.hxx>
#include <BRepAlgoAPI_Fuse.hxx>
#include <BRepBuilderAPI_MakeFace.hxx>
#include <BRepBuilderAPI_MakeVertex.hxx>
#include <BRepBuilderAPI_Transform.hxx>
#include <BRepPrimAPI_MakeCone.hxx>
#include <BRepPrimAPI_MakeCylinder.hxx>
#include <gp_Ax3.hxx>
#include <gp_Pln.hxx>
#include <gp_Pnt.hxx>
#include <gp_Vec.hxx>
#include <cmath>
#include <Standard_Stream.hxx>
#include <Standard_Failure.hxx>
#include <StdFail_NotDone.hxx>
#include <Standard_ErrorHandler.hxx> // CAREFUL ! position of this file is critic : see Lucien PIGNOLONI / OCC
#define HALF_LENGTH_MAIN_PIPE "Main pipe half length" //"Tuyau principal - demi longueur"
#define HALF_LENGTH_INCIDENT_PIPE "Incident pipe half length" //"Tuyau incident - demi longueur"
#define CIRCULAR_QUARTER_PIPE "Circular quarter of pipe" //"Circulaire - quart de tuyau"
#define THICKNESS "Thickness" //"Epaisseur"
#define FLANGE "Flange" // "Collerette"
#define CHAMFER_OR_FILLET "Chamfer or fillet" //"Chanfrein ou Raccord"
#define JUNCTION_FACE_1 "Junction 1" //"Face de jonction 1"
#define JUNCTION_FACE_2 "Junction 2" //"Face de jonction 2"
#define JUNCTION_FACE_3 "Junction 3" //"Face de jonction 3"
#define FIND_GROUPS_BY_POINTS 1
//=============================================================================
/*!
* Constructor
*/
//=============================================================================
GEOMImpl_IAdvancedOperations::GEOMImpl_IAdvancedOperations(GEOM_Engine* theEngine, int theDocID) :
GEOM_IOperations(theEngine, theDocID)
{
MESSAGE("GEOMImpl_IAdvancedOperations::GEOMImpl_IAdvancedOperations");
myBasicOperations = new GEOMImpl_IBasicOperations(GetEngine(), GetDocID());
myBooleanOperations = new GEOMImpl_IBooleanOperations(GetEngine(), GetDocID());
myShapesOperations = new GEOMImpl_IShapesOperations(GetEngine(), GetDocID());
myTransformOperations = new GEOMImpl_ITransformOperations(GetEngine(), GetDocID());
myBlocksOperations = new GEOMImpl_IBlocksOperations(GetEngine(), GetDocID());
my3DPrimOperations = new GEOMImpl_I3DPrimOperations(GetEngine(), GetDocID());
myLocalOperations = new GEOMImpl_ILocalOperations(GetEngine(), GetDocID());
myHealingOperations = new GEOMImpl_IHealingOperations(GetEngine(), GetDocID());
}
//=============================================================================
/*!
* Destructor
*/
//=============================================================================
GEOMImpl_IAdvancedOperations::~GEOMImpl_IAdvancedOperations()
{
MESSAGE("GEOMImpl_IAdvancedOperations::~GEOMImpl_IAdvancedOperations");
delete myBasicOperations;
delete myBooleanOperations;
delete myShapesOperations;
delete myTransformOperations;
delete myBlocksOperations;
delete my3DPrimOperations;
delete myLocalOperations;
delete myHealingOperations;
}
//=============================================================================
/*!
* SetPosition
*/
//=============================================================================
gp_Trsf GEOMImpl_IAdvancedOperations::GetPositionTrsf(double theL1, double theL2,
Handle(GEOM_Object) theP1,
Handle(GEOM_Object) theP2,
Handle(GEOM_Object) theP3)
{
// Old Local Coordinates System oldLCS
gp_Pnt P0(0, 0, 0);
gp_Pnt P1(-theL1, 0, 0);
gp_Pnt P2(theL1, 0, 0);
gp_Pnt P3(0, 0, theL2);
gp_Dir oldX(gp_Vec(P1, P2));
gp_Dir oldZ(gp_Vec(P0, P3));
gp_Ax3 oldLCS(P0, oldZ, oldX);
// New Local Coordinates System newLCS
double LocX, LocY, LocZ;
gp_Pnt newP1 = BRep_Tool::Pnt(TopoDS::Vertex(theP1->GetValue()));
gp_Pnt newP2 = BRep_Tool::Pnt(TopoDS::Vertex(theP2->GetValue()));
gp_Pnt newP3 = BRep_Tool::Pnt(TopoDS::Vertex(theP3->GetValue()));
LocX = (newP1.X() + newP2.X()) / 2.;
LocY = (newP1.Y() + newP2.Y()) / 2.;
LocZ = (newP1.Z() + newP2.Z()) / 2.;
gp_Pnt newO(LocX, LocY, LocZ);
gp_Dir newX(gp_Vec(newP1, newP2)); // P1P2 Vector
gp_Dir newZ(gp_Vec(newO, newP3)); // OP3 Vector
gp_Ax3 newLCS = gp_Ax3(newO, newZ, newX);
gp_Trsf aTrsf;
aTrsf.SetDisplacement(oldLCS, newLCS);
return aTrsf;
}
//=============================================================================
/*!
* CheckCompatiblePosition
*
*/
//=============================================================================
bool GEOMImpl_IAdvancedOperations::CheckCompatiblePosition(double& theL1, double& theL2,
Handle(GEOM_Object) theP1,
Handle(GEOM_Object) theP2,
Handle(GEOM_Object) theP3,
double theTolerance)
{
SetErrorCode(KO);
gp_Pnt P1 = BRep_Tool::Pnt(TopoDS::Vertex(theP1->GetValue()));
gp_Pnt P2 = BRep_Tool::Pnt(TopoDS::Vertex(theP2->GetValue()));
gp_Pnt P3 = BRep_Tool::Pnt(TopoDS::Vertex(theP3->GetValue()));
double d12 = P1.Distance(P2);
double d13 = P1.Distance(P3);
double d23 = P2.Distance(P3);
// double d2 = newO.Distance(P3);
if (Abs(d12) <= Precision::Confusion()) {
SetErrorCode("Junctions points P1 and P2 are identical");
return false;
}
if (Abs(d13) <= Precision::Confusion()) {
SetErrorCode("Junctions points P1 and P3 are identical");
return false;
}
if (Abs(d23) <= Precision::Confusion()) {
SetErrorCode("Junctions points P2 and P3 are identical");
return false;
}
double newL1 = 0.5 * d12;
double newL2 = sqrt(pow(d13,2)-pow(newL1,2));
//
// theL1*(1-theTolerance) <= newL1 <= theL1*(1+theTolerance)
//
if (fabs(newL1 - theL1) > Precision::Approximation()) {
if ( (newL1 * (1 - theTolerance) -theL1 <= Precision::Approximation()) &&
(newL1 * (1 + theTolerance) -theL1 >= Precision::Approximation()) ) {
// std::cerr << "theL1 = newL1" << std::endl;
theL1 = newL1;
} else {
theL1 = -1;
SetErrorCode("Dimension for main pipe (L1) is incompatible with new position");
return false;
}
}
//
// theL2*(1-theTolerance) <= newL2 <= theL2*(1+theTolerance)
//
if (fabs(newL2 - theL2) > Precision::Approximation()) {
if ( (newL2 * (1 - theTolerance) -theL2 <= Precision::Approximation()) &&
(newL2 * (1 + theTolerance) -theL2 >= Precision::Approximation()) ) {
theL2 = newL2;
} else {
theL2 = -1;
SetErrorCode("Dimension for incident pipe (L2) is incompatible with new position");
return false;
}
}
SetErrorCode(OK);
return true;
}
//=============================================================================
/*!
* Generate the propagation groups of a Pipe T-Shape used for hexa mesh
*/
//=============================================================================
bool GEOMImpl_IAdvancedOperations::MakeGroups(Handle(GEOM_Object) theShape, int shapeType,
double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theH, double theW, double theRF,
Handle(TColStd_HSequenceOfTransient) theSeq,
gp_Trsf aTrsf)
{
SetErrorCode(KO);
if (theShape.IsNull()) return false;
TopoDS_Shape aShape = theShape->GetValue();
if (aShape.IsNull()) {
SetErrorCode("Shape is not defined");
return false;
}
gp_Trsf aTrsfInv = aTrsf.Inverted();
// int expectedGroups = 0;
// if (shapeType == TSHAPE_BASIC)
// if (Abs(theR2+theW2-theR1-theW1) <= Precision::Approximation())
// expectedGroups = 10;
// else
// expectedGroups = 11;
// else if (shapeType == TSHAPE_CHAMFER || shapeType == TSHAPE_FILLET)
// expectedGroups = 12;
double aR1Ext = theR1 + theW1;
double aR2Ext = theR2 + theW2;
/////////////////////////
//// Groups of Faces ////
/////////////////////////
//
// Comment the following lines when GetInPlace bug is solved
// == BEGIN
// Workaround of GetInPlace bug
// Create a bounding box that fits the shape
Handle(GEOM_Object) aBox = my3DPrimOperations->MakeBoxDXDYDZ(2*theL1, 2*aR1Ext, aR1Ext+theL2);
aBox->GetLastFunction()->SetDescription("");
myTransformOperations->TranslateDXDYDZ(aBox, -theL1, -aR1Ext, -aR1Ext);
aBox->GetLastFunction()->SetDescription("");
// Apply transformation to box
BRepBuilderAPI_Transform aTransformationBox(aBox->GetValue(), aTrsf, Standard_False);
TopoDS_Shape aBoxShapeTrsf = aTransformationBox.Shape();
aBox->GetLastFunction()->SetValue(aBoxShapeTrsf);
// Get the shell of the box
Handle(GEOM_Object) aShell = Handle(GEOM_Object)::DownCast
(myShapesOperations->MakeExplode(aBox, TopAbs_SHELL, true)->Value(1));
aBox->GetLastFunction()->SetDescription("");
aShell->GetLastFunction()->SetDescription("");
// Get the common shapes between shell and shape
Handle(GEOM_Object) aCommonCompound = myBooleanOperations->MakeBoolean (theShape, aShell, 1); // MakeCommon
if (aCommonCompound.IsNull()) {
SetErrorCode(myBooleanOperations->GetErrorCode());
return false;
}
aCommonCompound->GetLastFunction()->SetDescription("");
// Explode the faces of common shapes => 3 faces
Handle(TColStd_HSequenceOfTransient) aCommonFaces =
myShapesOperations->MakeExplode(aCommonCompound, TopAbs_FACE, true);
aCommonCompound->GetLastFunction()->SetDescription("");
std::list<Handle(GEOM_Object)> aCompoundOfFacesList;
for (int i=0 ; i<= aCommonFaces->Length()-4 ; i+=4) {
std::list<Handle(GEOM_Object)> aFacesList;
for (int j = 1 ; j <= 4 ; j++) {
Handle(GEOM_Object) aFace = Handle(GEOM_Object)::DownCast(aCommonFaces->Value(i+j)); // Junction faces
if (!aFace.IsNull()) {
aFace->GetLastFunction()->SetDescription("");
aFacesList.push_back(aFace);
}
}
Handle(GEOM_Object) aCompoundOfFaces = myShapesOperations->MakeCompound(aFacesList);
if (!aCompoundOfFaces.IsNull()) {
aCompoundOfFaces->GetLastFunction()->SetDescription("");
aCompoundOfFacesList.push_back(aCompoundOfFaces);
}
}
if (aCompoundOfFacesList.size() == 3) {
Handle(GEOM_Object) aPln1 = aCompoundOfFacesList.front();
aCompoundOfFacesList.pop_front();
Handle(GEOM_Object) aPln2 = aCompoundOfFacesList.front();
aCompoundOfFacesList.pop_front();
Handle(GEOM_Object) aPln3 = aCompoundOfFacesList.front();
aCompoundOfFacesList.pop_front();
// == END
//
// Uncomment the following lines when GetInPlace bug is solved
// == BEGIN
// Handle(GEOM_Object) aP1 = myBasicOperations->MakePointXYZ(-theL1, 0, 0);
// Handle(GEOM_Object) aP2 = myBasicOperations->MakePointXYZ(-0, 0, theL2);
// Handle(GEOM_Object) aP3 = myBasicOperations->MakePointXYZ(theL1, 0, 0);
// aP1->GetLastFunction()->SetDescription("");
// aP2->GetLastFunction()->SetDescription("");
// aP3->GetLastFunction()->SetDescription("");
// Handle(GEOM_Object) aV1 = myBasicOperations->MakeVectorDXDYDZ(-1, 0, 0);
// Handle(GEOM_Object) aV2 = myBasicOperations->MakeVectorDXDYDZ(0, 0, 1);
// Handle(GEOM_Object) aV3 = myBasicOperations->MakeVectorDXDYDZ(1, 0, 0);
// aV1->GetLastFunction()->SetDescription("");
// aV2->GetLastFunction()->SetDescription("");
// aV3->GetLastFunction()->SetDescription("");
// Handle(GEOM_Object) aPln1 = myBasicOperations->MakePlanePntVec(aP1, aV1, 2*(aR1Ext+theL2));
// Handle(GEOM_Object) aPln2 = myBasicOperations->MakePlanePntVec(aP2, aV2, 2*(aR2Ext));
// Handle(GEOM_Object) aPln3 = myBasicOperations->MakePlanePntVec(aP3, aV3, 2*(aR1Ext+theL2));
// aPln1->GetLastFunction()->SetDescription("");
// aPln2->GetLastFunction()->SetDescription("");
// aPln3->GetLastFunction()->SetDescription("");
//
// BRepBuilderAPI_Transform aTransformation1(aPln1->GetValue(), aTrsf, Standard_False);
// TopoDS_Shape aTrsf_Shape1 = aTransformation1.Shape();
// aPln1->GetLastFunction()->SetValue(aTrsf_Shape1);
// BRepBuilderAPI_Transform aTransformation2(aPln2->GetValue(), aTrsf, Standard_False);
// TopoDS_Shape aTrsf_Shape2 = aTransformation2.Shape();
// aPln2->GetLastFunction()->SetValue(aTrsf_Shape2);
// BRepBuilderAPI_Transform aTransformation3(aPln3->GetValue(), aTrsf, Standard_False);
// TopoDS_Shape aTrsf_Shape3 = aTransformation3.Shape();
// aPln3->GetLastFunction()->SetValue(aTrsf_Shape3);
// == END
//
Handle(GEOM_Object) junctionFaces1 = myShapesOperations->GetInPlace(theShape, aPln1);
if (junctionFaces1.IsNull())
junctionFaces1 = myShapesOperations->GetShapesOnShapeAsCompound
(aPln1, theShape, TopAbs_FACE, GEOMAlgo_ST_ONIN);
if (!junctionFaces1.IsNull()) {
junctionFaces1->GetLastFunction()->SetDescription("");
junctionFaces1->SetName("JUNCTION_FACE_1");
theSeq->Append(junctionFaces1);
}
else {
SetErrorCode("Junction face 1 not found");
// theSeq->Append(aPln1);
// return false;
}
Handle(GEOM_Object) junctionFaces2 = myShapesOperations->GetInPlace(theShape, aPln2);
if (junctionFaces2.IsNull())
junctionFaces2 = myShapesOperations->GetShapesOnShapeAsCompound
(aPln2, theShape, TopAbs_FACE, GEOMAlgo_ST_ONIN);
if (!junctionFaces2.IsNull()) {
junctionFaces2->GetLastFunction()->SetDescription("");
junctionFaces2->SetName("JUNCTION_FACE_2");
theSeq->Append(junctionFaces2);
}
else {
SetErrorCode("Junction face 2 not found");
// theSeq->Append(aPln2);
// return false;
}
Handle(GEOM_Object) junctionFaces3 = myShapesOperations->GetInPlace(theShape, aPln3);
if (junctionFaces3.IsNull())
junctionFaces3 = myShapesOperations->GetShapesOnShapeAsCompound
(aPln3, theShape, TopAbs_FACE, GEOMAlgo_ST_ONIN);
if (!junctionFaces3.IsNull()) {
junctionFaces3->GetLastFunction()->SetDescription("");
junctionFaces3->SetName("JUNCTION_FACE_3");
theSeq->Append(junctionFaces3);
}
else {
SetErrorCode("Junction face 3 not found");
// theSeq->Append(aPln3);
// return false;
}
// Comment the following lines when GetInPlace bug is solved
// == BEGIN
}
// == END
/////////////////////////
//// Groups of Edges ////
/////////////////////////
// Result of propagate
Handle(GEOM_Function) aFunction = theShape->GetLastFunction();
TCollection_AsciiString theDesc = aFunction->GetDescription();
Handle(TColStd_HSequenceOfTransient) aSeqPropagate = myBlocksOperations->Propagate(theShape);
if (aSeqPropagate.IsNull() || aSeqPropagate->Length() == 0) {
SetErrorCode("Propagation groups not found");
return false;
}
Standard_Integer aNbGroups = aSeqPropagate->Length();
// Recover previous description to get rid of Propagate dump
aFunction->SetDescription(theDesc);
#ifdef FIND_GROUPS_BY_POINTS
// BEGIN: new groups search
// W2 R2
// .----.-----.----.
// e| | | | |
// | | | | |
// . | | | .
// g / ''..| | |..'' \
// f / ''''''' \
// .---.--'.. | | | ..'--.---.
// |a \ '''...........''' / |
// |-------\------' | '------/-------.
// | \ | / |
// c| \ | / |
// | R1 \ | / |
// | \ | / |
// ._________________|_________________.
// | L1 | |
// | | |
// | | |
// b| | |
// | | |
// |-----------------|-----------------|
// | W1 | |
// '-----------------'-----------------'
// d
// "Thickness" group (a)
gp_Pnt aPntA (-theL1, 0, theR1 + theW1/2.);
aPntA.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexA (aPntA);
TopoDS_Vertex aVertA = TopoDS::Vertex(mkVertexA.Shape());
TopoDS_Shape anEdgeA = GEOMUtils::GetEdgeNearPoint(aShape, aVertA);
// "Circular quarter of pipe" group (b)
gp_Pnt aPntB (-theL1, -aR1Ext * sin(M_PI/4.), -aR1Ext * sin(M_PI/4.));
aPntB.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexB (aPntB);
TopoDS_Vertex aVertB = TopoDS::Vertex(mkVertexB.Shape());
TopoDS_Shape anEdgeB = GEOMUtils::GetEdgeNearPoint(aShape, aVertB);
// "Circular quarter of pipe" group (c)
gp_Pnt aPntC (-theL1, -aR1Ext * sin(M_PI/4.), aR1Ext * sin(M_PI/4.));
aPntC.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexC (aPntC);
TopoDS_Vertex aVertC = TopoDS::Vertex(mkVertexC.Shape());
TopoDS_Shape anEdgeC = GEOMUtils::GetEdgeNearPoint(aShape, aVertC);
// "Main pipe half length" group (d)
gp_Pnt aPntD (-theL1/2., 0, -aR1Ext);
aPntD.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexD (aPntD);
TopoDS_Vertex aVertD = TopoDS::Vertex(mkVertexD.Shape());
TopoDS_Shape anEdgeD = GEOMUtils::GetEdgeNearPoint(aShape, aVertD);
// "Incident pipe half length" group (e)
double aTol10 = Precision::Confusion() * 10.;
gp_Pnt aPntE (-aR2Ext, 0, theL2 - aTol10);
aPntE.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexE (aPntE);
TopoDS_Vertex aVertE = TopoDS::Vertex(mkVertexE.Shape());
TopoDS_Shape anEdgeE = GEOMUtils::GetEdgeNearPoint(aShape, aVertE);
// "Flange" group (f)
double aFx = - aR2Ext - aTol10;
if (shapeType == TSHAPE_CHAMFER)
aFx -= theW;
else if (shapeType == TSHAPE_FILLET)
aFx -= theRF;
gp_Pnt aPntF (aFx, 0, aR1Ext);
aPntF.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexF (aPntF);
TopoDS_Vertex aVertF = TopoDS::Vertex(mkVertexF.Shape());
TopoDS_Shape anEdgeF = GEOMUtils::GetEdgeNearPoint(aShape, aVertF);
// "Chamfer or Fillet" group (g)
TopoDS_Shape anEdgeG;
if (shapeType == TSHAPE_CHAMFER) {
gp_Pnt aPntG (-aR2Ext - theW/2., 0, aR1Ext + theH/2.);
aPntG.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexG (aPntG);
TopoDS_Vertex aVertG = TopoDS::Vertex(mkVertexG.Shape());
anEdgeG = GEOMUtils::GetEdgeNearPoint(aShape, aVertG);
}
else if (shapeType == TSHAPE_FILLET) {
gp_Pnt aPntG (-aR2Ext - theRF/2., 0, aR1Ext + theRF/2.);
aPntG.Transform(aTrsf);
BRepBuilderAPI_MakeVertex mkVertexG (aPntG);
TopoDS_Vertex aVertG = TopoDS::Vertex(mkVertexG.Shape());
anEdgeG = GEOMUtils::GetEdgeNearPoint(aShape, aVertG);
}
for (int i = 1 ; i <= aNbGroups; i++) {
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(aSeqPropagate->Value(i));
if (aGroup.IsNull())
continue;
TopoDS_Shape aGroupShape = aGroup->GetValue();
TopTools_IndexedMapOfShape anEdgesMap;
TopExp::MapShapes(aGroupShape, TopAbs_EDGE, anEdgesMap);
if (anEdgesMap.Contains(anEdgeA)) { // a
aGroup->SetName("THICKNESS");
theSeq->Append(aGroup);
}
else if (anEdgesMap.Contains(anEdgeB)) { // b
aGroup->SetName("CIRCULAR_QUARTER_PIPE");
theSeq->Append(aGroup);
}
else if (anEdgesMap.Contains(anEdgeC)) { // c
aGroup->SetName("CIRCULAR_QUARTER_PIPE");
theSeq->Append(aGroup);
}
else if (anEdgesMap.Contains(anEdgeD)) { // d
aGroup->SetName("HALF_LENGTH_MAIN_PIPE");
theSeq->Append(aGroup);
}
else if (anEdgesMap.Contains(anEdgeE)) { // e
aGroup->SetName("HALF_LENGTH_INCIDENT_PIPE");
theSeq->Append(aGroup);
}
else if (anEdgesMap.Contains(anEdgeF)) { // f
aGroup->SetName("FLANGE");
theSeq->Append(aGroup);
}
else if (shapeType == TSHAPE_CHAMFER) { // g
if (anEdgesMap.Contains(anEdgeG)) {
aGroup->SetName("CHAMFER");
theSeq->Append(aGroup);
}
}
else if (shapeType == TSHAPE_FILLET) { // g
if (anEdgesMap.Contains(anEdgeG)) {
aGroup->SetName("FILLET");
theSeq->Append(aGroup);
}
}
else {
}
}
// END: new groups search
#else
bool addGroup;
bool circularFoundAndAdded = false;
bool circularFound10 = false;
bool incidentPipeFound = false;
bool mainPipeFound = false;
bool mainPipeFoundAndAdded = false;
bool radialFound =false;
bool flangeFound = false;
bool flangeFoundAndAdded = false;
bool chamferOrFilletFound = false;
for (int i = 1 ; i <= aNbGroups; i++) {
addGroup = false;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(aSeqPropagate->Value(i));
if (aGroup.IsNull())
continue;
TopoDS_Shape aGroupShape = aGroup->GetValue();
BRepBuilderAPI_Transform aTransformationShapeInv (aGroupShape, aTrsfInv, Standard_False);
TopoDS_Shape aGroupShapeTrsfInv = aTransformationShapeInv.Shape();
TopTools_IndexedMapOfShape anEdgesMap;
TopExp::MapShapes(aGroupShapeTrsfInv,TopAbs_EDGE, anEdgesMap);
Standard_Integer nbEdges = anEdgesMap.Extent();
if (shapeType == TSHAPE_BASIC) {
if ((nbEdges >= 21) || /*R1Ext = R2Ext*/(nbEdges == 17)) { // 17, 17+8*{1,2,3}, 21, 21+8*{1,2,3}
addGroup = true;
aGroup->SetName("THICKNESS");
}
else if (nbEdges == 6) {
if (!circularFoundAndAdded) {
circularFoundAndAdded = true;
addGroup = true;
aGroup->SetName("CIRCULAR_QUARTER_PIPE");
}
}
else if (nbEdges == 8) {
incidentPipeFound = true;
mainPipeFound = false;
radialFound = false;
flangeFound = false;
TopExp_Explorer Ex(aGroupShapeTrsfInv,TopAbs_VERTEX);
while (Ex.More()) {
gp_Pnt aP = BRep_Tool::Pnt(TopoDS::Vertex(Ex.Current()));
double x=aP.X(), y=aP.Y(), z=aP.Z();
if ((Abs(x) > aR2Ext + Precision::Confusion()) ||
(Abs(y) > aR2Ext + Precision::Confusion())) {
incidentPipeFound = false;
}
if ( z < -Precision::Confusion()) {
// length of main pipe
mainPipeFound = true;
if (!mainPipeFoundAndAdded) {
mainPipeFoundAndAdded = true;
addGroup = true;
aGroup->SetName("HALF_LENGTH_MAIN_PIPE");
}
}
else if (Abs(x) > (theL1-Precision::Confusion())) {
// discretisation circulaire
radialFound = true;
if (!circularFoundAndAdded) {
circularFoundAndAdded = true;
addGroup = true;
aGroup->SetName("CIRCULAR_QUARTER_PIPE");
}
}
Ex.Next();
}
if (incidentPipeFound) {
addGroup = true;
aGroup->SetName("HALF_LENGTH_INCIDENT_PIPE");
}
if (!addGroup && (!incidentPipeFound &&
!radialFound &&
!mainPipeFound &&
!flangeFound)) {
// Flange (collerette)
flangeFound = true;
addGroup = true;
aGroup->SetName("FLANGE");
}
}
else
continue;
}
else if (shapeType == TSHAPE_CHAMFER || shapeType == TSHAPE_FILLET) {
if (nbEdges >= 25) { // 25, 25+8, 25+16, 25+24
addGroup = true;
aGroup->SetName("THICKNESS");
}
else if ((nbEdges == 10) || (nbEdges == 6)) {
if (!circularFoundAndAdded) {
addGroup = true;
circularFoundAndAdded = true;
aGroup->SetName("CIRCULAR_QUARTER_PIPE");
if (nbEdges == 10) {
circularFound10 = true;
}
}
else if (!circularFound10 && nbEdges == 10) {
circularFound10 = true;
addGroup = true;
aGroup->SetName("CIRCULAR_QUARTER_PIPE");
}
}
else if (nbEdges == 8) {
incidentPipeFound = true;
mainPipeFound = true;
flangeFound = false;
bool isNearZ0 = false;
bool isBelowZ0 = false;
TopExp_Explorer Ex (aGroupShapeTrsfInv,TopAbs_VERTEX);
while (Ex.More()) {
gp_Pnt aP = BRep_Tool::Pnt(TopoDS::Vertex(Ex.Current()));
double x=aP.X(), y=aP.Y(), z=aP.Z();
// tuy_princ_long_avant & tuy_princ_long_apres
//bool isMain = (((z < Precision::Confusion()) || (x < Precision::Confusion())) &&
// ((y <= aR1Ext + Precision::Confusion()) ||
// (y <= -(aR1Ext + Precision::Confusion())) ||
// (y <= theR1 + Precision::Confusion()) ||
// (y == -(theR1 + Precision::Confusion()))));
bool isMain = ((z < Precision::Confusion() || x < Precision::Confusion()) &&
(fabs(y) > theR1 - Precision::Confusion() ||
fabs(y) < Precision::Confusion()));
if (!isMain) {
mainPipeFound = false;
}
// collerette
//if (z < Precision::Confusion() && !isMain) {
// flangeFound = true;
// if (!flangeFoundAndAdded) {
// flangeFoundAndAdded = true;
// addGroup = true;
// aGroup->SetName("FLANGE");
// }
//}
if (fabs(z) < Precision::Confusion()) isNearZ0 = true;
if (z < - Precision::Confusion()) isBelowZ0 = true;
// tuyau incident
if ((Abs(x) > aR2Ext + Precision::Confusion()) ||
(Abs(y) > aR2Ext + Precision::Confusion())) {
incidentPipeFound = false;
}
Ex.Next();
}
if (mainPipeFound) {
addGroup = true;
aGroup->SetName("HALF_LENGTH_MAIN_PIPE");
}
if (incidentPipeFound) {
addGroup = true;
aGroup->SetName("HALF_LENGTH_INCIDENT_PIPE");
}
if (isNearZ0 && !isBelowZ0) {
flangeFound = true;
if (!flangeFoundAndAdded) {
flangeFoundAndAdded = true;
addGroup = true;
aGroup->SetName("FLANGE");
}
}
if (!addGroup && (!incidentPipeFound &&
!mainPipeFound &&
!flangeFound &&
!chamferOrFilletFound)) {
addGroup = true;
chamferOrFilletFound = true;
if (shapeType == TSHAPE_CHAMFER)
aGroup->SetName("CHAMFER");
else
aGroup->SetName("FILLET");
}
}
else
continue;
}
// Add group to the list
if (addGroup)
theSeq->Append(aGroup);
}
#endif
SetErrorCode(OK);
return true;
}
bool GEOMImpl_IAdvancedOperations::MakePipeTShapePartition(Handle(GEOM_Object) theShape,
double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theH, double theW,
double theRF, bool isNormal)
{
SetErrorCode(KO);
// Build tools for partition operation:
// 1 face and 2 planes
// Face
Handle(GEOM_Object) arete_intersect_int;
Handle(GEOM_Object) wire_t, wire_t2, face_t, face_t2;
Handle(GEOM_Object) chan_racc;
Handle(GEOM_Object) vi1, vi2;
Handle(GEOM_Object) Te3;
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
Handle(GEOM_Object) Vector_Z = myBasicOperations->MakeVectorDXDYDZ(0, 0, 1);
Vector_Z->GetLastFunction()->SetDescription("");
// Useful values
double aSize = 2*(theL1 + theL2);
double aR1Ext = theR1 + theW1;
double aR2Ext = theR2 + theW2;
double theVertCylinderRadius = aR2Ext + theW + theRF;
double theHoriCylinderRadius = aR1Ext + theH + theRF;
// Common edges on internal cylinder
Handle(GEOM_Object) box_i = my3DPrimOperations->MakeBoxDXDYDZ(theR2, theR2, theR1);
box_i->GetLastFunction()->SetDescription("");
box_i = myTransformOperations->TranslateDXDYDZ(box_i, -theR2, -theR2, 0);
box_i->GetLastFunction()->SetDescription("");
Handle(GEOM_Function) aFunction = theShape->GetLastFunction();
TCollection_AsciiString theDesc = aFunction->GetDescription();
Handle(TColStd_HSequenceOfTransient) edges_i =
myShapesOperations->GetShapesOnBox(box_i, theShape, TopAbs_EDGE, GEOMAlgo_ST_IN);
// Recover previous description to get rid of Propagate dump
aFunction->SetDescription(theDesc);
if (edges_i.IsNull() || edges_i->Length() == 0) {
SetErrorCode("Internal edges not found");
return false;
}
for (int i=1; i<=edges_i->Length();i++) {
Handle(GEOM_Object) anObj = Handle(GEOM_Object)::DownCast(edges_i->Value(i));
anObj->GetLastFunction()->SetDescription("");
}
arete_intersect_int = Handle(GEOM_Object)::DownCast(edges_i->Value(1));
// search for vertices located on both internal pipes
aFunction = theShape->GetLastFunction();
theDesc = aFunction->GetDescription();
Handle(TColStd_HSequenceOfTransient) vertices_i =
myShapesOperations->GetShapesOnBox(box_i, theShape, TopAbs_VERTEX, GEOMAlgo_ST_ONIN);
// Recover previous description to get rid of Propagate dump
aFunction->SetDescription(theDesc);
if (vertices_i.IsNull() || vertices_i->Length() == 0) {
SetErrorCode("Internal vertices not found");
return false;
}
for (int i = 1; i <= vertices_i->Length(); i++) {
Handle(GEOM_Object) v = Handle(GEOM_Object)::DownCast(vertices_i->Value(i));
v->GetLastFunction()->SetDescription("");
TopoDS_Vertex aVertex = TopoDS::Vertex(v->GetValue());
gp_Pnt aP = BRep_Tool::Pnt(aVertex);
// std::cout << "Coords: " << aP.X() << ", " << aP.Y() << ", " << aP.Z() << std::endl;
if (Abs(aP.X()) <= Precision::Confusion()) {
if (Abs(aP.Y()) - theR1 <= Precision::Confusion()) {
vi1 = v;
}
} else if (Abs(aP.Y()) <= Precision::Confusion()) {
if (Abs(aP.X()) - theR1 <= Precision::Confusion()) {
vi2 = v;
}
}
}
std::list<Handle(GEOM_Object)> theShapes;
if (isNormal) {
Handle(GEOM_Object) ve1, ve2;
Handle(GEOM_Object) box_e = my3DPrimOperations->MakeBoxDXDYDZ(aR2Ext, aR2Ext, aR1Ext);
box_e->GetLastFunction()->SetDescription("");
box_e = myTransformOperations->TranslateDXDYDZ(box_e, -aR2Ext, -aR2Ext, 0);
box_e->GetLastFunction()->SetDescription("");
// Common edges on external cylinder
aFunction = theShape->GetLastFunction();
theDesc = aFunction->GetDescription();
Handle(TColStd_HSequenceOfTransient) edges_e =
myShapesOperations->GetShapesOnBox(box_e, theShape, TopAbs_EDGE, GEOMAlgo_ST_IN);
// Recover previous description to get rid of Propagate dump
aFunction->SetDescription(theDesc);
if (edges_e.IsNull() || edges_e->Length() == 0) {
SetErrorCode("External edges not found");
return false;
}
for (int i=1; i<=edges_e->Length();i++) {
Handle(GEOM_Object) anObj = Handle(GEOM_Object)::DownCast(edges_e->Value(i));
anObj->GetLastFunction()->SetDescription("");
}
// search for vertices located on both external pipes
aFunction = theShape->GetLastFunction();
theDesc = aFunction->GetDescription();
Handle(TColStd_HSequenceOfTransient) vertices_e =
myShapesOperations->GetShapesOnBox(box_e, theShape, TopAbs_VERTEX, GEOMAlgo_ST_ONIN);
// Recover previous description to get rid of Propagate dump
aFunction->SetDescription(theDesc);
if (vertices_e.IsNull() || vertices_e->Length() == 0) {
SetErrorCode("External vertices not found");
return false;
}
for (int i = 1; i <= vertices_e->Length(); i++) {
Handle(GEOM_Object) v = Handle(GEOM_Object)::DownCast(vertices_e->Value(i));
v->GetLastFunction()->SetDescription("");
TopoDS_Vertex aVertex = TopoDS::Vertex(v->GetValue());
gp_Pnt aP = BRep_Tool::Pnt(aVertex);
// std::cout << "Coords: " << aP.X() << ", " << aP.Y() << ", " << aP.Z() << std::endl;
if (Abs(aP.X()) <= Precision::Confusion()) {
if (Abs(aP.Y()) - theR2 > Precision::Confusion()) {
ve1 = v;
}
} else if (Abs(aP.Y()) <= Precision::Confusion()) {
if (Abs(aP.X()) - theR2 > Precision::Confusion()) {
ve2 = v;
}
}
if ( !ve1.IsNull() && !ve2.IsNull())
break;
}
Handle(GEOM_Object) edge_e1, edge_e2;
edge_e1 = myBasicOperations->MakeLineTwoPnt(ve1, vi1);
if (edge_e1.IsNull()) {
SetErrorCode("Edge 1 could not be built");
return false;
}
edge_e2 = myBasicOperations->MakeLineTwoPnt(ve2, vi2);
if (edge_e2.IsNull()) {
SetErrorCode("Edge 2 could not be built");
return false;
}
edge_e1->GetLastFunction()->SetDescription("");
edge_e2->GetLastFunction()->SetDescription("");
std::list<Handle(GEOM_Object)> edge_e_elist;
edge_e_elist.push_back(arete_intersect_int);
edge_e_elist.push_back(edge_e1);
edge_e_elist.push_back(Handle(GEOM_Object)::DownCast(edges_e->Value(1)));
edge_e_elist.push_back(edge_e2);
wire_t = myShapesOperations->MakeWire(edge_e_elist, 1e-7);
if (wire_t.IsNull()) {
SetErrorCode("Impossible to build wire");
return false;
}
wire_t->GetLastFunction()->SetDescription("");
face_t = myShapesOperations->MakeFace(wire_t, false);
if (face_t.IsNull()) {
SetErrorCode("Impossible to build face");
return false;
}
face_t->GetLastFunction()->SetDescription("");
theShapes.push_back(theShape);
theShapes.push_back(vi1);
theShapes.push_back(vi2);
theShapes.push_back(ve1);
theShapes.push_back(ve2);
theShapes.push_back(edge_e1);
theShapes.push_back(edge_e2);
theShapes.push_back(wire_t);
theShapes.push_back(face_t);
}
else {
Handle(GEOM_Object) P1, P2, P3, P4, P5, P6;
int idP1, idP2, idP3, idP4;
int PZX, PZY;
double ZX=0, ZY=0;
std::vector<int> LX;
std::vector<int> LY;
Handle(GEOM_Object) box_e = my3DPrimOperations->MakeBoxDXDYDZ
(theVertCylinderRadius, theVertCylinderRadius, theHoriCylinderRadius);
box_e->GetLastFunction()->SetDescription("");
box_e = myTransformOperations->TranslateDXDYDZ
(box_e, -theVertCylinderRadius, -theVertCylinderRadius, 0);
box_e->GetLastFunction()->SetDescription("");
aFunction = theShape->GetLastFunction();
theDesc = aFunction->GetDescription();
Handle(TColStd_HSequenceOfTransient) extremVertices =
myShapesOperations->GetShapesOnBox(box_e, theShape, TopAbs_VERTEX, GEOMAlgo_ST_ONIN);
// Recover previous description to get rid of Propagate dump
aFunction->SetDescription(theDesc);
if (extremVertices.IsNull() || extremVertices->Length() == 0) {
if (theRF == 0)
SetErrorCode("Vertices on chamfer not found");
else
SetErrorCode("Vertices on fillet not found");
return false;
}
theShapes.push_back(theShape);
theShapes.push_back(box_e);
if (extremVertices->Length() != 6) {
// for (int i=1; i<=extremVertices->Length(); i++){
// theShapes.push_back(Handle(GEOM_Object)::DownCast(extremVertices->Value(i)));
// }
// Handle(GEOM_Object) aCompound = myShapesOperations->MakeCompound(theShapes);
// TopoDS_Shape aCompoundShape = aCompound->GetValue();
// theShape->GetLastFunction()->SetValue(aCompoundShape);
SetErrorCode("Bad number of vertices on chamfer found");
return false;
}
for (int i=1; i<=extremVertices->Length(); i++){
Handle(GEOM_Object) aV = Handle(GEOM_Object)::DownCast(extremVertices->Value(i));
aV->GetLastFunction()->SetDescription("");
gp_Pnt aP = BRep_Tool::Pnt(TopoDS::Vertex(aV->GetValue()));
if (Abs(aP.X()) <= Precision::Confusion()) {
if (Abs(aP.Y()) - theR2 > Precision::Confusion()) {
LX.push_back(i);
if (aP.Z()-ZX > Precision::Confusion()) {
ZX = aP.Z();
PZX = i;
}
}
}
else {
if (Abs(aP.X()) - theR2 > Precision::Confusion()) {
LY.push_back(i);
if (aP.Z() - ZY > Precision::Confusion()) {
ZY = aP.Z();
PZY = i;
}
}
}
}
idP2 = PZX;
idP4 = PZY;
idP1 = LX.at(0);
if (LX.at(0) == PZX)
idP1 = LX.at(1);
idP3 = LY.at(0);
if (LY.at(0) == PZY)
idP3 = LY.at(1);
P1 = Handle(GEOM_Object)::DownCast(extremVertices->Value(idP1));
P2 = Handle(GEOM_Object)::DownCast(extremVertices->Value(idP2));
P3 = Handle(GEOM_Object)::DownCast(extremVertices->Value(idP3));
P4 = Handle(GEOM_Object)::DownCast(extremVertices->Value(idP4));
Handle(GEOM_Object) Cote_1 = myBasicOperations->MakeLineTwoPnt(P1, vi1);
if (Cote_1.IsNull()) {
SetErrorCode("Impossible to build edge in thickness");
return false;
}
Cote_1->GetLastFunction()->SetDescription("");
Handle(GEOM_Object) Cote_2 = myBasicOperations->MakeLineTwoPnt(vi2, P3);
if (Cote_2.IsNull()) {
SetErrorCode("Impossible to build edge in thickness");
return false;
}
Cote_2->GetLastFunction()->SetDescription("");
// edge_chan_princ = arete du chanfrein (ou raccord) sur le tuyau principal
// edge_chan_inc = arete du chanfrein (ou raccord) sur le tuyau incident
// std::cerr << "Getting chamfer edge on main pipe" << std::endl;
Handle(GEOM_Object) edge_chan_princ = myBlocksOperations->GetEdge(theShape, P1, P3);
if (edge_chan_princ.IsNull()) {
SetErrorCode("Impossible to find edge on main pipe");
return false;
}
edge_chan_princ->GetLastFunction()->SetDescription("");
Handle(GEOM_Object) edge_chan_inc = myBlocksOperations->GetEdge(theShape, P2, P4);
if (edge_chan_inc.IsNull()) {
SetErrorCode("Impossible to find edge on incident pipe");
return false;
}
edge_chan_inc->GetLastFunction()->SetDescription("");
std::list<Handle(GEOM_Object)> edgeList1;
edgeList1.push_back(edge_chan_princ);
edgeList1.push_back(Cote_1);
edgeList1.push_back(arete_intersect_int);
edgeList1.push_back(Cote_2);
// std::cerr << "Creating wire 1" << std::endl;
wire_t = myShapesOperations->MakeWire(edgeList1, 1e-7);
if (wire_t.IsNull()) {
SetErrorCode("Impossible to build wire");
return false;
}
wire_t->GetLastFunction()->SetDescription("");
// std::cerr << "Creating face 1" << std::endl;
face_t = myShapesOperations->MakeFace(wire_t, false);
if (face_t.IsNull()) {
SetErrorCode("Impossible to build face");
return false;
}
face_t->GetLastFunction()->SetDescription("");
theShapes.push_back(face_t);
gp_Pnt aP2 = BRep_Tool::Pnt(TopoDS::Vertex(P2->GetValue()));
gp_Pnt aP5 = BRep_Tool::Pnt(TopoDS::Vertex(vi1->GetValue()));
double deltaZ = aP2.Z() - aP5.Z();
// std::cerr << "Creating new point from vi1 with deltaZ = " << deltaZ << std::endl;
Handle(GEOM_Object) P5bis = myTransformOperations->TranslateDXDYDZCopy(vi1, 0, 0, deltaZ);
if (P5bis.IsNull()) {
SetErrorCode("Impossible to translate vertex");
return false;
}
P5bis->GetLastFunction()->SetDescription("");
gp_Pnt aP4 = BRep_Tool::Pnt(TopoDS::Vertex(P4->GetValue()));
gp_Pnt aP6 = BRep_Tool::Pnt(TopoDS::Vertex(vi2->GetValue()));
deltaZ = aP4.Z() - aP6.Z();
// std::cerr << "Creating new point from vi2 with deltaZ = " << deltaZ << std::endl;
Handle(GEOM_Object) P6bis = myTransformOperations->TranslateDXDYDZCopy(vi2, 0, 0, deltaZ);
if (P6bis.IsNull()) {
SetErrorCode("Impossible to translate vertex");
return false;
}
P6bis->GetLastFunction()->SetDescription("");
// std::cerr << "Creating new line 1 from 2 previous points" << std::endl;
Handle(GEOM_Object) Cote_3 = myBasicOperations->MakeLineTwoPnt(P5bis, P2);
if (Cote_3.IsNull()) {
SetErrorCode("Impossible to build edge in thickness");
return false;
}
Cote_3->GetLastFunction()->SetDescription("");
// std::cerr << "Creating new line 2 from 2 previous points" << std::endl;
Handle(GEOM_Object) Cote_4 = myBasicOperations->MakeLineTwoPnt(P6bis, P4);
if (Cote_4.IsNull()) {
SetErrorCode("Impossible to build edge in thickness");
return false;
}
Cote_4->GetLastFunction()->SetDescription("");
// std::cerr << "Creating new line 3 from 2 previous points" << std::endl;
Handle(GEOM_Object) Cote_5 = myBasicOperations->MakeLineTwoPnt(P5bis, P6bis);
if (Cote_4.IsNull()) {
SetErrorCode("Impossible to build edge in thickness");
return false;
}
Cote_5->GetLastFunction()->SetDescription("");
//std::list<Handle(GEOM_Object)> edgeList2;
//edgeList2.push_back(edge_chan_inc);
//edgeList2.push_back(Cote_3);
//edgeList2.push_back(Cote_5);
//edgeList2.push_back(Cote_4);
// std::cerr << "Creating wire 2" << std::endl;
//wire_t2 = myShapesOperations->MakeWire(edgeList2, 1e-7);
//if (wire_t2.IsNull()) {
// SetErrorCode("Impossible to build wire");
// return false;
//}
//wire_t2->GetLastFunction()->SetDescription("");
// std::cerr << "Creating face 2" << std::endl;
//face_t2 = myShapesOperations->MakeFace(wire_t2, false);
// Mantis issue 0021682
face_t2 = my3DPrimOperations->MakePrismVecH(edge_chan_inc, Cote_4, - (theR2 + theW2));
//face_t2 = my3DPrimOperations->MakePrismVecH(edge_chan_inc, Cote_4, - 2.0*theR2);
if (face_t2.IsNull()) {
SetErrorCode("Impossible to build face");
return false;
}
face_t2->GetLastFunction()->SetDescription("");
theShapes.push_back(face_t2);
}
// Planes
Handle(GEOM_Object) aP0 = myBasicOperations->MakePointXYZ(0, 0, 0);
Handle(GEOM_Object) aVZ = myBasicOperations->MakeVectorDXDYDZ(0, 0, 1);
Handle(GEOM_Object) aVXZ = myBasicOperations->MakeVectorDXDYDZ(aR1Ext, 0, 0.5*(theL1+theVertCylinderRadius));
Handle(GEOM_Object) aPlnOZ = myBasicOperations->MakePlanePntVec(aP0, aVZ, aSize);
Handle(GEOM_Object) aPlnOXZ = myBasicOperations->MakePlanePntVec(aP0, aVXZ, aSize);
aP0->GetLastFunction()->SetDescription("");
aVZ->GetLastFunction()->SetDescription("");
aVXZ->GetLastFunction()->SetDescription("");
aPlnOZ->GetLastFunction()->SetDescription("");
aPlnOXZ->GetLastFunction()->SetDescription("");
theShapes.push_back(aPlnOZ);
theShapes.push_back(aPlnOXZ);
// Partition
Handle(TColStd_HSequenceOfTransient) partitionShapes = new TColStd_HSequenceOfTransient;
Handle(TColStd_HSequenceOfTransient) theTools = new TColStd_HSequenceOfTransient;
Handle(TColStd_HSequenceOfTransient) theKeepInside = new TColStd_HSequenceOfTransient;
Handle(TColStd_HSequenceOfTransient) theRemoveInside = new TColStd_HSequenceOfTransient;
Handle(TColStd_HArray1OfInteger) theMaterials;
partitionShapes->Append(theShape);
theTools->Append(aPlnOZ);
if (Abs(aR1Ext - aR2Ext) > Precision::Confusion())
theTools->Append(aPlnOXZ);
theTools->Append(face_t);
if (!isNormal)
theTools->Append(face_t2);
Te3 = myBooleanOperations->MakePartition
(partitionShapes, theTools, theKeepInside, theRemoveInside,
TopAbs_SOLID, false, theMaterials, 0, false);
if (Te3.IsNull()) {
SetErrorCode("Impossible to build partition of TShape");
return false;
}
Te3->GetLastFunction()->SetDescription("");
// Last verification: result should be a block
std::list<GEOMImpl_IBlocksOperations::BCError> errList;
if (!myBlocksOperations->CheckCompoundOfBlocks(Te3,errList)) {
SetErrorCode("TShape is not a compound of block");
return false;
}
// // BEGIN Compound of created shapes - Only for debug purpose
// theShapes.clear();
// theShapes.push_back(theShape);
// theShapes.push_back(aPlnOZ);
// if (Abs(aR1Ext - aR2Ext) > Precision::Confusion() )
// theShapes.push_back(aPlnOXZ);
// theShapes.push_back(face_t);
// if (!isNormal)
// theShapes.push_back(face_t2);
//
// Handle(GEOM_Object) aCompound = myShapesOperations->MakeCompound(theShapes);
// TopoDS_Shape aCompoundShape = aCompound->GetValue();
// theShape->GetLastFunction()->SetValue(aCompoundShape);
// // END Compound of created shapes - Only for debug purpose
TopoDS_Shape aShape = Te3->GetValue();
theShape->GetLastFunction()->SetValue(aShape);
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return false;
}
SetErrorCode(OK);
return true;
}
// Mirror and glue faces
bool GEOMImpl_IAdvancedOperations::MakePipeTShapeMirrorAndGlue(Handle(GEOM_Object) theShape,
double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2)
{
SetErrorCode(KO);
// Useful values
double aSize = 2*(theL1 + theL2);
double aR1Ext = theR1 + theW1;
// Planes
Handle(GEOM_Object) aP0 = myBasicOperations->MakePointXYZ(0, 0, 0);
aP0->GetLastFunction()->SetDescription("");
Handle(GEOM_Object) aVX = myBasicOperations->MakeVectorDXDYDZ(1, 0, 0);
Handle(GEOM_Object) aVY = myBasicOperations->MakeVectorDXDYDZ(0, 1, 0);
aVX->GetLastFunction()->SetDescription("");
aVY->GetLastFunction()->SetDescription("");
Handle(GEOM_Object) aPlane_OX = myBasicOperations->MakePlanePntVec(aP0, aVX, 2*(aR1Ext + theL2));
Handle(GEOM_Object) aPlane_OY = myBasicOperations->MakePlanePntVec(aP0, aVY, aSize);
aPlane_OX->GetLastFunction()->SetDescription("");
aPlane_OY->GetLastFunction()->SetDescription("");
Handle(GEOM_Object) Te4 = myTransformOperations->MirrorPlaneCopy(theShape, aPlane_OX);
if (Te4.IsNull()) {
SetErrorCode("Impossible to build mirror of quarter TShape");
return false;
}
Handle(GEOM_Object) Te5 = myTransformOperations->MirrorPlaneCopy(theShape, aPlane_OY);
if (Te5.IsNull()) {
SetErrorCode("Impossible to build mirror of half TShape");
return false;
}
Handle(GEOM_Object) Te6 = myTransformOperations->MirrorPlaneCopy(Te4, aPlane_OY);
if (Te6.IsNull()) {
SetErrorCode("Impossible to build mirror of half TShape");
return false;
}
std::list<Handle(GEOM_Object)> aShapesList;
aShapesList.push_back(theShape);
aShapesList.push_back(Te4);
aShapesList.push_back(Te5);
aShapesList.push_back(Te6);
Handle(GEOM_Object) Te7 = myShapesOperations->MakeCompound(aShapesList);
if (Te7.IsNull()) {
SetErrorCode("Impossible to build compound");
return false;
}
// Copy source shape
TopoDS_Shape aShapeCopy;
TColStd_IndexedDataMapOfTransientTransient aMapTShapes;
TNaming_CopyShape::CopyTool(Te7->GetValue(), aMapTShapes, aShapeCopy);
Handle(GEOM_Object) Te8 = myShapesOperations->MakeGlueFaces(Te7, 1e-7, true);
if (Te8.IsNull()) {
SetErrorCode("Impossible to glue faces of TShape");
return false;
}
TopoDS_Shape aShape = Te8->GetValue();
BRepCheck_Analyzer anAna (aShape, Standard_True);
if (!anAna.IsValid()) {
// Try to do gluing with the tolerance equal to maximal
// tolerance of vertices of the source shape.
Standard_Real aTolMax = -RealLast();
for (TopExp_Explorer ExV (aShapeCopy, TopAbs_VERTEX); ExV.More(); ExV.Next()) {
TopoDS_Vertex aVertex = TopoDS::Vertex(ExV.Current());
Standard_Real aTol = BRep_Tool::Tolerance(aVertex);
if (aTol > aTolMax) {
aTolMax = aTol;
}
}
// Perform gluing
Te7->GetLastFunction()->SetValue(aShapeCopy);
Te8 = myShapesOperations->MakeGlueFaces(Te7, aTolMax, true);
if (Te8.IsNull()) {
SetErrorCode("Impossible to glue faces of TShape");
return false;
}
aShape = Te8->GetValue();
}
theShape->GetLastFunction()->SetValue(aShape);
Te4->GetLastFunction()->SetDescription("");
Te5->GetLastFunction()->SetDescription("");
Te6->GetLastFunction()->SetDescription("");
Te7->GetLastFunction()->SetDescription("");
Te8->GetLastFunction()->SetDescription("");
SetErrorCode(OK);
return true;
}
//=======================================================================
//function : MakePipeTShapeThicknessReduction
//purpose : Static method. Add thiskness reduction elements at the three
// open ends of the T-Shape.
//=======================================================================
TopoDS_Shape GEOMImpl_IAdvancedOperations::MakePipeTShapeThicknessReduction
(TopoDS_Shape theShape,
double r1, double w1, double l1,
double r2, double w2, double l2,
double rL, double wL, double ltransL, double lthinL,
double rR, double wR, double ltransR, double lthinR,
double rI, double wI, double ltransI, double lthinI,
bool fuseReductions)
{
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
//
// ---------------------.
// W \
// ---------------------. \
// ^ \ '-----------------.
// |R \ Wthin |
// | '-----------------'
// v Rthin
// --.--.--.--.--.--.--.--.--.--.--.--.--.--.--
// Ltrans Lthin
TopoDS_Shape aResult = theShape;
double aTol = Precision::Confusion();
gp_Vec aVX = gp::DX(), aVZ = gp::DZ();
// Left reduction (rL, wL, ltransL, lthinL)
if (rL > aTol && wL > aTol && ltransL > aTol) {
gp_Pnt aPLeft (-l1, 0, 0);
gp_Ax2 anAxesLeft (aPLeft, -aVX, aVZ);
TopoDS_Shape aReductionLeft = GEOMImpl_IAdvancedOperations::MakeThicknessReduction
(anAxesLeft, r1, w1, rL, wL, ltransL, lthinL, fuseReductions);
if (fuseReductions) {
BRepAlgoAPI_Fuse fuseL (aResult, aReductionLeft);
if (!fuseL.IsDone())
StdFail_NotDone::Raise("Cannot fuse Te with left reduction");
aResult = fuseL.Shape();
}
else {
BRep_Builder B;
TopoDS_Compound C;
B.MakeCompound(C);
B.Add(C, aResult);
B.Add(C, aReductionLeft);
aResult = C;
}
}
// Right reduction
if (rR > aTol && wR > aTol && ltransR > aTol) {
gp_Pnt aPRight (l1, 0, 0);
gp_Ax2 anAxesRight (aPRight, aVX, aVZ);
TopoDS_Shape aReductionRight = GEOMImpl_IAdvancedOperations::MakeThicknessReduction
(anAxesRight, r1, w1, rR, wR, ltransR, lthinR, fuseReductions);
if (fuseReductions) {
BRepAlgoAPI_Fuse fuseR (aResult, aReductionRight);
if (!fuseR.IsDone())
StdFail_NotDone::Raise("Cannot fuse Te with right reduction");
aResult = fuseR.Shape();
}
else {
BRep_Builder B;
TopoDS_Compound C;
B.MakeCompound(C);
B.Add(C, aResult);
B.Add(C, aReductionRight);
aResult = C;
}
}
// Incident reduction
if (rI > aTol && wI > aTol && ltransI > aTol) {
gp_Pnt aPInci (0, 0, l2);
gp_Ax2 anAxesInci (aPInci, aVZ, aVX);
TopoDS_Shape aReductionInci = GEOMImpl_IAdvancedOperations::MakeThicknessReduction
(anAxesInci, r2, w2, rI, wI, ltransI, lthinI, fuseReductions);
if (fuseReductions) {
BRepAlgoAPI_Fuse fuseInci (aResult, aReductionInci);
if (!fuseInci.IsDone())
StdFail_NotDone::Raise("Cannot fuse Te with incident reduction");
aResult = fuseInci.Shape();
}
else {
BRep_Builder B;
TopoDS_Compound C;
B.MakeCompound(C);
B.Add(C, aResult);
B.Add(C, aReductionInci);
aResult = C;
}
}
// Get rid of extra compounds
TopTools_ListOfShape listShapeRes;
GEOMUtils::AddSimpleShapes(aResult, listShapeRes);
aResult = listShapeRes.First(); // useful for the case "fuseReductions == true"
if (!fuseReductions && listShapeRes.Extent() > 1) {
// Simplify T-Shape compound (get rid of sub-compounds) and glue duplicated faces
BRep_Builder B;
TopoDS_Compound C;
B.MakeCompound(C);
TopTools_ListIteratorOfListOfShape itSub (listShapeRes);
for (; itSub.More(); itSub.Next())
B.Add(C, itSub.Value());
// GlueFaces
aResult = GEOMImpl_GlueDriver::GlueFaces(C, Precision::Confusion(), Standard_True);
}
return aResult;
}
//=======================================================================
//function : MakeThicknessReduction
//purpose : Static method. Create one thickness reduction element.
//=======================================================================
TopoDS_Shape GEOMImpl_IAdvancedOperations::MakeThicknessReduction (gp_Ax2 theAxes,
const double R, const double W,
const double Rthin, const double Wthin,
const double Ltrans, const double Lthin,
bool fuse)
{
double aTol = Precision::Confusion();
if (Rthin < aTol || Wthin < aTol || Ltrans < aTol) {
StdFail_NotDone::Raise("Cannot build thickness reduction: too small values");
}
bool isThinPart = (Lthin > aTol);
// .
// W |\
// . \
// ^ \ '-----------------.
// |R \| | Wthin
// | '-----------------'
// v Rthin
// --.--.--.--.--.--.--.--.--.--.--.--.--> theAxes.Direction()
// Ltrans Lthin
double RExt = R + W;
double RthinExt = Rthin + Wthin;
gp_Dir aNormal = theAxes.Direction();
gp_Dir anXDir = theAxes.XDirection();
gp_Pnt aPntCyl (theAxes.Location().XYZ() + aNormal.XYZ()*Ltrans);
gp_Ax2 anAxesCyl (aPntCyl, aNormal, anXDir);
// Build the transition part
BRepPrimAPI_MakeCone ConeExt (theAxes, RExt, RthinExt, Ltrans);
BRepPrimAPI_MakeCone ConeInt (theAxes, R, Rthin, Ltrans);
ConeExt.Build();
ConeInt.Build();
if (!ConeExt.IsDone() || !ConeInt.IsDone())
StdFail_NotDone::Raise("Cannot build cones of thickness reduction");
BRepAlgoAPI_Cut cut1 (ConeExt.Shape(), ConeInt.Shape());
if (!cut1.IsDone())
StdFail_NotDone::Raise("Coudn't build transition part of thickness reduction");
TopoDS_Shape aReduction = cut1.Shape();
// Build the thin part, if required
TopoDS_Shape aThinPart;
if (isThinPart) {
BRepPrimAPI_MakeCylinder CExt (anAxesCyl, RthinExt, Lthin);
BRepPrimAPI_MakeCylinder CInt (anAxesCyl, Rthin, Lthin);
CExt.Build();
CInt.Build();
if (!CExt.IsDone() || !CInt.IsDone())
StdFail_NotDone::Raise("Cannot build cylinders of thickness reduction");
BRepAlgoAPI_Cut cut2 (CExt.Shape(), CInt.Shape());
if (!cut2.IsDone())
StdFail_NotDone::Raise("Coudn't build thin part of thickness reduction");
aThinPart = cut2.Shape();
}
// Join parts
if (fuse) {
if (isThinPart) {
BRepAlgoAPI_Fuse fuse1 (aReduction, aThinPart);
if (!fuse1.IsDone())
StdFail_NotDone::Raise("Cannot fuse parts of thickness reduction");
aReduction = fuse1.Shape();
}
}
else {
// Partition the reduction on blocks
gp_Ax3 anAxesPln1 (aPntCyl, theAxes.XDirection(), aNormal);
gp_Ax3 anAxesPln2 (aPntCyl, theAxes.YDirection(), aNormal);
gp_Pln aPln1 (anAxesPln1);
gp_Pln aPln2 (anAxesPln2);
double aSize = Ltrans + Lthin + R + Rthin + Wthin; // to guarantee enough size in all directions
TopoDS_Shape aTool1 = BRepBuilderAPI_MakeFace(aPln1, -aSize, +aSize, -aSize, +aSize).Shape();
TopoDS_Shape aTool2 = BRepBuilderAPI_MakeFace(aPln2, -aSize, +aSize, -aSize, +aSize).Shape();
GEOMAlgo_Splitter PS;
PS.AddShape(aReduction);
if (isThinPart)
PS.AddShape(aThinPart);
PS.AddTool(aTool1);
PS.AddTool(aTool2);
PS.SetLimit(TopAbs_SOLID);
PS.Perform();
aReduction = PS.Shape();
}
return aReduction;
}
//=============================================================================
/*!
* MakePipeTShape
* \brief Create a T-shape object with specified caracteristics for the main and
* the incident pipes (radius, width, half-length).
* Center of the shape is (0,0,0). The main plane of the T-shape is XOY.
* \param theR1 Internal radius of main pipe
* \param theW1 Width of main pipe
* \param theL1 Half-length of main pipe
* \param theR2 Internal radius of incident pipe (R2 < R1)
* \param theW2 Width of incident pipe (R2+W2 < R1+W1)
* \param theL2 Half-length of incident pipe
* \param theHexMesh Boolean indicating if shape is prepared for hex mesh
* \return List of GEOM_Objects, containing the created shape and propagation groups.
*/
//=============================================================================
Handle(TColStd_HSequenceOfTransient)
GEOMImpl_IAdvancedOperations::MakePipeTShape(double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theRL, double theWL, double theLtransL, double theLthinL,
double theRR, double theWR, double theLtransR, double theLthinR,
double theRI, double theWI, double theLtransI, double theLthinI,
bool theHexMesh)
{
MESSAGE("GEOMImpl_IAdvancedOperations::MakePipeTShape");
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_BASIC);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_PipeTShapeDriver::GetID()) return NULL;
GEOMImpl_IPipeTShape aData (aFunction);
aData.SetR1(theR1);
aData.SetW1(theW1);
aData.SetL1(theL1);
aData.SetR2(theR2);
aData.SetW2(theW2);
aData.SetL2(theL2);
aData.SetHexMesh(theHexMesh);
bool isTRL = (theRL + theWL + theLtransL + theLthinL) > Precision::Confusion();
bool isTRR = (theRR + theWR + theLtransR + theLthinR) > Precision::Confusion();
bool isTRI = (theRI + theWI + theLtransI + theLthinI) > Precision::Confusion();
//Compute the resulting value
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("TShape driver failed");
return NULL;
}
if (theHexMesh) {
if (!MakePipeTShapePartition(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
if (!MakePipeTShapeMirrorAndGlue(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
}
if (isTRL || isTRR || isTRI) {
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
TopoDS_Shape aResShape =
MakePipeTShapeThicknessReduction(aShape->GetValue(), theR1, theW1, theL1, theR2, theW2, theL2,
theRL, theWL, theLtransL, theLthinL,
theRR, theWR, theLtransR, theLthinR,
theRI, theWI, theLtransI, theLthinI,
!theHexMesh);
aFunction->SetValue(aResShape);
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
Handle(TColStd_HSequenceOfTransient) aSeq = new TColStd_HSequenceOfTransient;
aSeq->Append(aShape);
if (theHexMesh) {
// Get the groups
try {
if (!MakeGroups(aShape, TSHAPE_BASIC, theR1, theW1, theL1, theR2, theW2, theL2,
0., 0., 0., aSeq, gp_Trsf()))
return NULL;
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
}
//Make a Python command
TCollection_AsciiString anEntry, aListRes("[");
// Iterate over the sequence aSeq
Standard_Integer aNbGroups = aSeq->Length();
Standard_Integer i = 1;
for (; i <= aNbGroups; i++) {
Handle(Standard_Transient) anItem = aSeq->Value(i);
if (anItem.IsNull()) continue;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(anItem);
if (aGroup.IsNull()) continue;
//Make a Python command
TDF_Tool::Entry(aGroup->GetEntry(), anEntry);
aListRes += anEntry + ", ";
}
aListRes.Trunc(aListRes.Length() - 2);
GEOM::TPythonDump pd (aFunction);
pd << aListRes.ToCString() << "] = geompy.MakePipeTShape("
<< theR1 << ", " << theW1 << ", " << theL1 << ", "
<< theR2 << ", " << theW2 << ", " << theL2 << ", "
<< theHexMesh;
// thickness reduction
if (isTRL)
pd << ", theRL=" << theRL << ", theWL=" << theWL
<< ", theLtransL=" << theLtransL << ", theLthinL=" << theLthinL;
if (isTRR)
pd << ", theRR=" << theRR << ", theWR=" << theWR
<< ", theLtransR=" << theLtransR << ", theLthinR=" << theLthinR;
if (isTRI)
pd << ", theRI=" << theRI << ", theWI=" << theWI
<< ", theLtransI=" << theLtransI << ", theLthinI=" << theLthinI;
pd << ")";
SetErrorCode(OK);
return aSeq;
}
//=============================================================================
/*!
* MakePipeTShapeWithPosition
* Create a T-shape object with specified caracteristics for the main and
* the incident pipes (radius, width, half-length).
* The extremities of the main pipe are located on junctions points P1 and P2.
* The extremity of the incident pipe is located on junction point P3.
* \param theR1 Internal radius of main pipe
* \param theW1 Width of main pipe
* \param theL1 Half-length of main pipe
* \param theR2 Internal radius of incident pipe (R2 < R1)
* \param theW2 Width of incident pipe (R2+W2 < R1+W1)
* \param theL2 Half-length of incident pipe
* \param theHexMesh Boolean indicating if shape is prepared for hex mesh
* \param theP1 1st junction point of main pipe
* \param theP2 2nd junction point of main pipe
* \param theP3 Junction point of incident pipe
* \return List of GEOM_Objects, containing the created shape and propagation groups..
*/
//=============================================================================
Handle(TColStd_HSequenceOfTransient)
GEOMImpl_IAdvancedOperations::MakePipeTShapeWithPosition
(double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theRL, double theWL, double theLtransL, double theLthinL,
double theRR, double theWR, double theLtransR, double theLthinR,
double theRI, double theWI, double theLtransI, double theLthinI,
bool theHexMesh,
Handle(GEOM_Object) theP1,
Handle(GEOM_Object) theP2,
Handle(GEOM_Object) theP3)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE);
/////////////////
// TSHAPE CODE
/////////////////
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_BASIC);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_PipeTShapeDriver::GetID()) return NULL;
// Check new position
if (!CheckCompatiblePosition(theL1, theL2, theP1, theP2, theP3, 0.01)) {
return NULL;
}
GEOMImpl_IPipeTShape aData(aFunction);
aData.SetR1(theR1);
aData.SetW1(theW1);
aData.SetL1(theL1);
aData.SetR2(theR2);
aData.SetW2(theW2);
aData.SetL2(theL2);
aData.SetHexMesh(theHexMesh);
bool isTRL = (theRL + theWL + theLtransL + theLthinL) > Precision::Confusion();
bool isTRR = (theRR + theWR + theLtransR + theLthinR) > Precision::Confusion();
bool isTRI = (theRI + theWI + theLtransI + theLthinI) > Precision::Confusion();
//Compute the resulting value
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("TShape driver failed");
return NULL;
}
if (theHexMesh) {
if (!MakePipeTShapePartition(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
if (!MakePipeTShapeMirrorAndGlue(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
}
if (isTRL || isTRR || isTRI) {
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
TopoDS_Shape aResShape =
MakePipeTShapeThicknessReduction(aShape->GetValue(), theR1, theW1, theL1, theR2, theW2, theL2,
theRL, theWL, theLtransL, theLthinL,
theRR, theWR, theLtransR, theLthinR,
theRI, theWI, theLtransI, theLthinI,
!theHexMesh);
aFunction->SetValue(aResShape);
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
TopoDS_Shape Te = aShape->GetValue();
// Set Position
gp_Trsf aTrsf = GetPositionTrsf(theL1, theL2, theP1, theP2, theP3);
BRepBuilderAPI_Transform aTransformation(Te, aTrsf, Standard_False);
TopoDS_Shape aTrsf_Shape = aTransformation.Shape();
aFunction->SetValue(aTrsf_Shape);
Handle(TColStd_HSequenceOfTransient) aSeq = new TColStd_HSequenceOfTransient;
aSeq->Append(aShape);
if (theHexMesh) {
// Get the groups
try {
if (!MakeGroups(aShape,TSHAPE_BASIC, theR1, theW1, theL1, theR2, theW2, theL2,
0., 0., 0., aSeq, aTrsf)) {
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
}
//Make a Python command
TCollection_AsciiString anEntry, aListRes("[");
// Iterate over the sequence aSeq
Standard_Integer aNbGroups = aSeq->Length();
Standard_Integer i = 1;
for (; i <= aNbGroups; i++) {
Handle(Standard_Transient) anItem = aSeq->Value(i);
if (anItem.IsNull()) continue;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(anItem);
if (aGroup.IsNull()) continue;
//Make a Python command
TDF_Tool::Entry(aGroup->GetEntry(), anEntry);
aListRes += anEntry + ", ";
}
aListRes.Trunc(aListRes.Length() - 2);
GEOM::TPythonDump pd (aFunction);
pd << aListRes.ToCString() << "] = geompy.MakePipeTShape("
<< theR1 << ", " << theW1 << ", " << theL1 << ", "
<< theR2 << ", " << theW2 << ", " << theL2 << ", "
<< theHexMesh << ", " << theP1 << ", " << theP2 << ", " << theP3;
// thickness reduction
if (isTRL)
pd << ", theRL=" << theRL << ", theWL=" << theWL
<< ", theLtransL=" << theLtransL << ", theLthinL=" << theLthinL;
if (isTRR)
pd << ", theRR=" << theRR << ", theWR=" << theWR
<< ", theLtransR=" << theLtransR << ", theLthinR=" << theLthinR;
if (isTRI)
pd << ", theRI=" << theRI << ", theWI=" << theWI
<< ", theLtransI=" << theLtransI << ", theLthinI=" << theLthinI;
pd << ")";
SetErrorCode(OK);
return aSeq;
}
//=============================================================================
/*!
* MakePipeTShapeChamfer
* Create a T-shape object with specified caracteristics for the main and
* the incident pipes (radius, width, half-length). A chamfer is created
* on the junction of the pipes.
* Center of the shape is (0,0,0). The main plane of the T-shape is XOY.
* \param theR1 Internal radius of main pipe
* \param theW1 Width of main pipe
* \param theL1 Half-length of main pipe
* \param theR2 Internal radius of incident pipe (R2 < R1)
* \param theW2 Width of incident pipe (R2+W2 < R1+W1)
* \param theL2 Half-length of incident pipe
* \param theH Height of chamfer.
* \param theW Width of chamfer.
* \param theHexMesh Boolean indicating if shape is prepared for hex mesh
* \return List of GEOM_Objects, containing the created shape and propagation groups.
*/
//=============================================================================
Handle(TColStd_HSequenceOfTransient)
GEOMImpl_IAdvancedOperations::MakePipeTShapeChamfer
(double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theRL, double theWL, double theLtransL, double theLthinL,
double theRR, double theWR, double theLtransR, double theLthinR,
double theRI, double theWI, double theLtransI, double theLthinI,
double theH, double theW,
bool theHexMesh)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_CHAMFER);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_PipeTShapeDriver::GetID()) return NULL;
GEOMImpl_IPipeTShape aData(aFunction);
aData.SetR1(theR1);
aData.SetW1(theW1);
aData.SetL1(theL1);
aData.SetR2(theR2);
aData.SetW2(theW2);
aData.SetL2(theL2);
aData.SetH(theH);
aData.SetW(theW);
aData.SetHexMesh(theHexMesh);
bool isTRL = (theRL + theWL + theLtransL + theLthinL) > Precision::Confusion();
bool isTRR = (theRR + theWR + theLtransR + theLthinR) > Precision::Confusion();
bool isTRI = (theRI + theWI + theLtransI + theLthinI) > Precision::Confusion();
//Compute the resulting value
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("TShape driver failed");
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
// BEGIN of chamfer
TopoDS_Shape aShapeShape = aShape->GetValue();
TopTools_IndexedMapOfShape anEdgesIndices;
TopExp::MapShapes(aShapeShape, anEdgesIndices);
// Common edges on external cylinders
Handle(GEOM_Object) box_e;
if (theHexMesh) {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(theR2+theW2, theR2+theW2, theR1+theW1);
}
else {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(2*(theR2+theW2), 2*(theR2+theW2), theR1+theW1);
}
box_e->GetLastFunction()->SetDescription("");
box_e = myTransformOperations->TranslateDXDYDZ(box_e, -(theR2+theW2), -(theR2+theW2), 0);
box_e->GetLastFunction()->SetDescription("");
Handle(TColStd_HSequenceOfInteger) edges_e =
myShapesOperations->GetShapesOnBoxIDs(box_e, aShape, TopAbs_EDGE, GEOMAlgo_ST_IN);
box_e->GetLastFunction()->SetDescription("");
if (edges_e.IsNull() || edges_e->Length() == 0) {
SetErrorCode("External edges not found");
return NULL;
}
int nbEdgesInChamfer = 0;
std::list<int> theEdges;
for (int i=1; i<=edges_e->Length();i++) {
int edgeID = edges_e->Value(i);
TopoDS_Shape theEdge = anEdgesIndices.FindKey(edgeID);
TopExp_Explorer Ex(theEdge,TopAbs_VERTEX);
int iv=0;
while (Ex.More()) {
iv ++;
gp_Pnt aPt = BRep_Tool::Pnt(TopoDS::Vertex(Ex.Current()));
if (Abs(aPt.Z() - (theR1+theW1)) <= Precision::Confusion()) {
nbEdgesInChamfer ++;
theEdges.push_back(edgeID);
}
Ex.Next();
}
if (theHexMesh && nbEdgesInChamfer == 1)
break;
}
Handle(GEOM_Object) aChamfer;
try {
aChamfer = myLocalOperations->MakeChamferEdges(aShape, theW, theH, theEdges);
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
if (aChamfer.IsNull()) {
SetErrorCode("Chamfer can not be computed on the given shape with the given parameters");
return NULL;
}
aChamfer->GetLastFunction()->SetDescription("");
TopoDS_Shape aChamferShape = aChamfer->GetValue();
aFunction->SetValue(aChamferShape);
// END of chamfer
if (theHexMesh) {
if (!MakePipeTShapePartition(aShape, theR1, theW1, theL1, theR2, theW2, theL2, theH, theW, 0, false))
return NULL;
if (!MakePipeTShapeMirrorAndGlue(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
}
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (isTRL || isTRR || isTRI) {
TopoDS_Shape aResShape =
MakePipeTShapeThicknessReduction(aShape->GetValue(), theR1, theW1, theL1, theR2, theW2, theL2,
theRL, theWL, theLtransL, theLthinL,
theRR, theWR, theLtransR, theLthinR,
theRI, theWI, theLtransI, theLthinI,
!theHexMesh);
aFunction->SetValue(aResShape);
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
Handle(TColStd_HSequenceOfTransient) aSeq = new TColStd_HSequenceOfTransient;
aSeq->Append(aShape);
if (theHexMesh) {
// Get the groups
try {
if (!MakeGroups(aShape, TSHAPE_CHAMFER, theR1, theW1, theL1, theR2, theW2, theL2,
theH, theW, 0., aSeq, gp_Trsf()))
return NULL;
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
}
//Make a Python command
TCollection_AsciiString anEntry, aListRes("[");
// Iterate over the sequence aSeq
Standard_Integer aNbGroups = aSeq->Length();
Standard_Integer i = 1;
for (; i <= aNbGroups; i++) {
Handle(Standard_Transient) anItem = aSeq->Value(i);
if (anItem.IsNull()) continue;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(anItem);
if (aGroup.IsNull()) continue;
//Make a Python command
TDF_Tool::Entry(aGroup->GetEntry(), anEntry);
aListRes += anEntry + ", ";
}
aListRes.Trunc(aListRes.Length() - 2);
GEOM::TPythonDump pd (aFunction);
pd << aListRes.ToCString() << "] = geompy.MakePipeTShapeChamfer("
<< theR1 << ", " << theW1 << ", " << theL1 << ", "
<< theR2 << ", " << theW2 << ", " << theL2 << ", "
<< theH << ", " << theW << ", " << theHexMesh;
// thickness reduction
if (isTRL)
pd << ", theRL=" << theRL << ", theWL=" << theWL
<< ", theLtransL=" << theLtransL << ", theLthinL=" << theLthinL;
if (isTRR)
pd << ", theRR=" << theRR << ", theWR=" << theWR
<< ", theLtransR=" << theLtransR << ", theLthinR=" << theLthinR;
if (isTRI)
pd << ", theRI=" << theRI << ", theWI=" << theWI
<< ", theLtransI=" << theLtransI << ", theLthinI=" << theLthinI;
pd << ")";
SetErrorCode(OK);
return aSeq;
}
//=============================================================================
/*!
* MakePipeTShapeChamferWithPosition
* Create a T-shape object with specified caracteristics for the main and
* the incident pipes (radius, width, half-length). A chamfer is created
* on the junction of the pipes.
* The extremities of the main pipe are located on junctions points P1 and P2.
* The extremity of the incident pipe is located on junction point P3.
* \param theR1 Internal radius of main pipe
* \param theW1 Width of main pipe
* \param theL1 Half-length of main pipe
* \param theR2 Internal radius of incident pipe (R2 < R1)
* \param theW2 Width of incident pipe (R2+W2 < R1+W1)
* \param theL2 Half-length of incident pipe
* \param theH Height of chamfer.
* \param theW Width of chamfer.
* \param theHexMesh Boolean indicating if shape is prepared for hex mesh
* \param theP1 1st junction point of main pipe
* \param theP2 2nd junction point of main pipe
* \param theP3 Junction point of incident pipe
* \return List of GEOM_Objects, containing the created shape and propagation groups.
*/
//=============================================================================
Handle(TColStd_HSequenceOfTransient)
GEOMImpl_IAdvancedOperations::MakePipeTShapeChamferWithPosition
(double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theRL, double theWL, double theLtransL, double theLthinL,
double theRR, double theWR, double theLtransR, double theLthinR,
double theRI, double theWI, double theLtransI, double theLthinI,
double theH, double theW,
bool theHexMesh,
Handle(GEOM_Object) theP1,
Handle(GEOM_Object) theP2,
Handle(GEOM_Object) theP3)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_CHAMFER);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_PipeTShapeDriver::GetID()) return NULL;
// Check new position
if (!CheckCompatiblePosition(theL1, theL2, theP1, theP2, theP3, 0.01)) {
return NULL;
}
GEOMImpl_IPipeTShape aData(aFunction);
aData.SetR1(theR1);
aData.SetW1(theW1);
aData.SetL1(theL1);
aData.SetR2(theR2);
aData.SetW2(theW2);
aData.SetL2(theL2);
aData.SetH(theH);
aData.SetW(theW);
aData.SetHexMesh(theHexMesh);
bool isTRL = (theRL + theWL + theLtransL + theLthinL) > Precision::Confusion();
bool isTRR = (theRR + theWR + theLtransR + theLthinR) > Precision::Confusion();
bool isTRI = (theRI + theWI + theLtransI + theLthinI) > Precision::Confusion();
//Compute the resulting value
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("TShape driver failed");
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
// BEGIN of chamfer
TopoDS_Shape aShapeShape = aShape->GetValue();
TopTools_IndexedMapOfShape anEdgesIndices;
TopExp::MapShapes(aShapeShape, anEdgesIndices);
// Common edges on external cylinders
Handle(GEOM_Object) box_e;
if (theHexMesh) {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(theR2+theW2, theR2+theW2, theR1+theW1);
}
else {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(2*(theR2+theW2), 2*(theR2+theW2), theR1+theW1);
}
box_e->GetLastFunction()->SetDescription("");
box_e = myTransformOperations->TranslateDXDYDZ(box_e, -(theR2+theW2), -(theR2+theW2), 0);
box_e->GetLastFunction()->SetDescription("");
Handle(TColStd_HSequenceOfInteger) edges_e =
myShapesOperations->GetShapesOnBoxIDs(box_e, aShape, TopAbs_EDGE, GEOMAlgo_ST_IN);
box_e->GetLastFunction()->SetDescription("");
if (edges_e.IsNull() || edges_e->Length() == 0) {
SetErrorCode("External edges not found");
return NULL;
}
int nbEdgesInChamfer = 0;
std::list<int> theEdges;
for (int i=1; i<=edges_e->Length();i++) {
int edgeID = edges_e->Value(i);
TopoDS_Shape theEdge = anEdgesIndices.FindKey(edgeID);
TopExp_Explorer Ex(theEdge,TopAbs_VERTEX);
while (Ex.More()) {
gp_Pnt aPt = BRep_Tool::Pnt(TopoDS::Vertex(Ex.Current()));
if (Abs(aPt.Z() - (theR1+theW1)) <= Precision::Confusion()) {
nbEdgesInChamfer ++;
theEdges.push_back(edgeID);
}
Ex.Next();
}
if (theHexMesh && nbEdgesInChamfer == 1)
break;
}
Handle(GEOM_Object) aChamfer;
try {
aChamfer = myLocalOperations->MakeChamferEdges(aShape, theW, theH, theEdges);
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
if (aChamfer.IsNull()) {
SetErrorCode("Chamfer can not be computed on the given shape with the given parameters");
return NULL;
}
aChamfer->GetLastFunction()->SetDescription("");
TopoDS_Shape aChamferShape = aChamfer->GetValue();
aFunction->SetValue(aChamferShape);
// END of chamfer
if (theHexMesh) {
if (!MakePipeTShapePartition(aShape, theR1, theW1, theL1, theR2, theW2, theL2, theH, theW, 0, false))
return NULL;
if (!MakePipeTShapeMirrorAndGlue(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
}
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (isTRL || isTRR || isTRI) {
TopoDS_Shape aResShape =
MakePipeTShapeThicknessReduction(aShape->GetValue(), theR1, theW1, theL1, theR2, theW2, theL2,
theRL, theWL, theLtransL, theLthinL,
theRR, theWR, theLtransR, theLthinR,
theRI, theWI, theLtransI, theLthinI,
!theHexMesh);
aFunction->SetValue(aResShape);
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
// Set Position
gp_Trsf aTrsf = GetPositionTrsf(theL1, theL2, theP1, theP2, theP3);
BRepBuilderAPI_Transform aTransformation (aShape->GetValue(), aTrsf, Standard_False);
TopoDS_Shape aTrsf_Shape = aTransformation.Shape();
aFunction->SetValue(aTrsf_Shape);
Handle(TColStd_HSequenceOfTransient) aSeq = new TColStd_HSequenceOfTransient;
aSeq->Append(aShape);
if (theHexMesh) {
// Get the groups
try {
if (!MakeGroups(aShape, TSHAPE_CHAMFER, theR1, theW1, theL1, theR2, theW2, theL2,
theH, theW, 0., aSeq, aTrsf))
return NULL;
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
}
//Make a Python command
TCollection_AsciiString anEntry, aListRes("[");
// Iterate over the sequence aSeq
Standard_Integer aNbGroups = aSeq->Length();
Standard_Integer i = 1;
for (; i <= aNbGroups; i++) {
Handle(Standard_Transient) anItem = aSeq->Value(i);
if (anItem.IsNull()) continue;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(anItem);
if (aGroup.IsNull()) continue;
//Make a Python command
TDF_Tool::Entry(aGroup->GetEntry(), anEntry);
aListRes += anEntry + ", ";
}
aListRes.Trunc(aListRes.Length() - 2);
GEOM::TPythonDump pd (aFunction);
pd << aListRes.ToCString() << "] = geompy.MakePipeTShapeChamfer("
<< theR1 << ", " << theW1 << ", " << theL1 << ", "
<< theR2 << ", " << theW2 << ", " << theL2 << ", "
<< theH << ", " << theW << ", " << theHexMesh << ", "
<< theP1 << ", " << theP2 << ", " << theP3;
// thickness reduction
if (isTRL)
pd << ", theRL=" << theRL << ", theWL=" << theWL
<< ", theLtransL=" << theLtransL << ", theLthinL=" << theLthinL;
if (isTRR)
pd << ", theRR=" << theRR << ", theWR=" << theWR
<< ", theLtransR=" << theLtransR << ", theLthinR=" << theLthinR;
if (isTRI)
pd << ", theRI=" << theRI << ", theWI=" << theWI
<< ", theLtransI=" << theLtransI << ", theLthinI=" << theLthinI;
pd << ")";
SetErrorCode(OK);
return aSeq;
}
//=============================================================================
/*!
* MakePipeTShapeFillet
* Create a T-shape object with specified caracteristics for the main and
* the incident pipes (radius, width, half-length). A fillet is created
* on the junction of the pipes.
* Center of the shape is (0,0,0). The main plane of the T-shape is XOY.
* \param theR1 Internal radius of main pipe
* \param theW1 Width of main pipe
* \param theL1 Half-length of main pipe
* \param theR2 Internal radius of incident pipe (R2 < R1)
* \param theW2 Width of incident pipe (R2+W2 < R1+W1)
* \param theL2 Half-length of incident pipe
* \param theRF Radius of curvature of fillet.
* \param theHexMesh Boolean indicating if shape is prepared for hex mesh
* \return List of GEOM_Objects, containing the created shape and propagation groups.
*/
//=============================================================================
Handle(TColStd_HSequenceOfTransient)
GEOMImpl_IAdvancedOperations::MakePipeTShapeFillet
(double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theRL, double theWL, double theLtransL, double theLthinL,
double theRR, double theWR, double theLtransR, double theLthinR,
double theRI, double theWI, double theLtransI, double theLthinI,
double theRF, bool theHexMesh)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_FILLET);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_PipeTShapeDriver::GetID()) return NULL;
GEOMImpl_IPipeTShape aData(aFunction);
aData.SetR1(theR1);
aData.SetW1(theW1);
aData.SetL1(theL1);
aData.SetR2(theR2);
aData.SetW2(theW2);
aData.SetL2(theL2);
aData.SetRF(theRF);
aData.SetHexMesh(theHexMesh);
bool isTRL = (theRL + theWL + theLtransL + theLthinL) > Precision::Confusion();
bool isTRR = (theRR + theWR + theLtransR + theLthinR) > Precision::Confusion();
bool isTRI = (theRI + theWI + theLtransI + theLthinI) > Precision::Confusion();
//Compute the resulting value
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("TShape driver failed");
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
// BEGIN of fillet
TopoDS_Shape aShapeShape = aShape->GetValue();
TopTools_IndexedMapOfShape anEdgesIndices;
TopExp::MapShapes(aShapeShape, anEdgesIndices);
// Common edges on external cylinders
Handle(GEOM_Object) box_e;
if (theHexMesh) {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(theR2+theW2, theR2+theW2, theR1+theW1);
}
else {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(2*(theR2+theW2), 2*(theR2+theW2), theR1+theW1);
}
box_e->GetLastFunction()->SetDescription("");
box_e = myTransformOperations->TranslateDXDYDZ(box_e, -(theR2+theW2), -(theR2+theW2), 0);
box_e->GetLastFunction()->SetDescription("");
Handle(TColStd_HSequenceOfInteger) edges_e =
myShapesOperations->GetShapesOnBoxIDs(box_e, aShape, TopAbs_EDGE, GEOMAlgo_ST_IN);
box_e->GetLastFunction()->SetDescription("");
if (edges_e.IsNull() || edges_e->Length() == 0) {
SetErrorCode("External edges not found");
return NULL;
}
int nbEdgesInFillet = 0;
std::list<int> theEdges;
for (int i=1; i<=edges_e->Length();i++) {
int edgeID = edges_e->Value(i);
TopoDS_Shape theEdge = anEdgesIndices.FindKey(edgeID);
TopExp_Explorer Ex(theEdge,TopAbs_VERTEX);
while (Ex.More()) {
gp_Pnt aPt = BRep_Tool::Pnt(TopoDS::Vertex(Ex.Current()));
if (Abs(aPt.Z() - (theR1+theW1)) <= Precision::Confusion()) {
nbEdgesInFillet ++;
theEdges.push_back(edgeID);
}
Ex.Next();
}
if (theHexMesh && nbEdgesInFillet == 1)
break;
}
Handle(GEOM_Object) aFillet;
try {
aFillet = myLocalOperations->MakeFilletEdges(aShape, theRF, theEdges);
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
if (aFillet.IsNull()) {
//SetErrorCode("Fillet can not be computed on the given shape with the given parameters");
SetErrorCode(myLocalOperations->GetErrorCode());
return NULL;
}
aFillet->GetLastFunction()->SetDescription("");
TopoDS_Shape aFilletShape = aFillet->GetValue();
aFunction->SetValue(aFilletShape);
// END of fillet
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - BEGIN (1)
// the following block, when enabled, leads to partitioning problems
#if 0
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - END (1)
// BEGIN: Limit tolerances (debug)
Handle(GEOM_Object) aCorr1 = myHealingOperations->LimitTolerance(aShape, 1e-07);
TopoDS_Shape aCorr1Shape = aCorr1->GetValue();
aShape->GetLastFunction()->SetValue(aCorr1Shape);
aCorr1->GetLastFunction()->SetDescription("");
// END: Limit tolerances (debug)
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - BEGIN (2)
#endif
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - END (2)
if (theHexMesh) {
if (!MakePipeTShapePartition(aShape, theR1, theW1, theL1, theR2, theW2, theL2, 0, 0, theRF, false))
return NULL;
if (!MakePipeTShapeMirrorAndGlue(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
}
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (isTRL || isTRR || isTRI) {
TopoDS_Shape aResShape =
MakePipeTShapeThicknessReduction(aShape->GetValue(), theR1, theW1, theL1, theR2, theW2, theL2,
theRL, theWL, theLtransL, theLthinL,
theRR, theWR, theLtransR, theLthinR,
theRI, theWI, theLtransI, theLthinI,
!theHexMesh);
aFunction->SetValue(aResShape);
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
Handle(TColStd_HSequenceOfTransient) aSeq = new TColStd_HSequenceOfTransient;
aSeq->Append(aShape);
if (theHexMesh) {
// Get the groups
try {
if (!MakeGroups(aShape, TSHAPE_FILLET, theR1, theW1, theL1, theR2, theW2, theL2,
0., 0., theRF, aSeq, gp_Trsf()))
return NULL;
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
}
//Make a Python command
TCollection_AsciiString anEntry, aListRes("[");
// Iterate over the sequence aSeq
Standard_Integer aNbGroups = aSeq->Length();
Standard_Integer i = 1;
for (; i <= aNbGroups; i++) {
Handle(Standard_Transient) anItem = aSeq->Value(i);
if (anItem.IsNull()) continue;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(anItem);
if (aGroup.IsNull()) continue;
//Make a Python command
TDF_Tool::Entry(aGroup->GetEntry(), anEntry);
aListRes += anEntry + ", ";
}
aListRes.Trunc(aListRes.Length() - 2);
GEOM::TPythonDump pd (aFunction);
pd << aListRes.ToCString() << "] = geompy.MakePipeTShapeFillet("
<< theR1 << ", " << theW1 << ", " << theL1 << ", "
<< theR2 << ", " << theW2 << ", " << theL2 << ", "
<< theRF << ", " << theHexMesh;
// thickness reduction
if (isTRL)
pd << ", theRL=" << theRL << ", theWL=" << theWL
<< ", theLtransL=" << theLtransL << ", theLthinL=" << theLthinL;
if (isTRR)
pd << ", theRR=" << theRR << ", theWR=" << theWR
<< ", theLtransR=" << theLtransR << ", theLthinR=" << theLthinR;
if (isTRI)
pd << ", theRI=" << theRI << ", theWI=" << theWI
<< ", theLtransI=" << theLtransI << ", theLthinI=" << theLthinI;
pd << ")";
SetErrorCode(OK);
return aSeq;
}
//=============================================================================
/*!
* MakePipeTShapeFilletWithPosition
* \brief Create a T-shape object with specified caracteristics for the main and
* the incident pipes (radius, width, half-length). A fillet is created
* on the junction of the pipes.
* The extremities of the main pipe are located on junctions points P1 and P2.
* The extremity of the incident pipe is located on junction point P3.
* \param theR1 Internal radius of main pipe
* \param theW1 Width of main pipe
* \param theL1 Half-length of main pipe
* \param theR2 Internal radius of incident pipe (R2 < R1)
* \param theW2 Width of incident pipe (R2+W2 < R1+W1)
* \param theL2 Half-length of incident pipe
* \param theRF Radius of curvature of fillet
* \param theHexMesh Boolean indicating if shape is prepared for hex mesh
* \param theP1 1st junction point of main pipe
* \param theP2 2nd junction point of main pipe
* \param theP3 Junction point of incident pipe
* \return List of GEOM_Objects, containing the created shape and propagation groups.
*/
//=============================================================================
Handle(TColStd_HSequenceOfTransient)
GEOMImpl_IAdvancedOperations::MakePipeTShapeFilletWithPosition
(double theR1, double theW1, double theL1,
double theR2, double theW2, double theL2,
double theRL, double theWL, double theLtransL, double theLthinL,
double theRR, double theWR, double theLtransR, double theLthinR,
double theRI, double theWI, double theLtransI, double theLthinI,
double theRF, bool theHexMesh,
Handle(GEOM_Object) theP1,
Handle(GEOM_Object) theP2,
Handle(GEOM_Object) theP3)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_FILLET);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_PipeTShapeDriver::GetID()) return NULL;
// Check new position
if (!CheckCompatiblePosition(theL1, theL2, theP1, theP2, theP3, 0.01)) {
return NULL;
}
GEOMImpl_IPipeTShape aData(aFunction);
aData.SetR1(theR1);
aData.SetW1(theW1);
aData.SetL1(theL1);
aData.SetR2(theR2);
aData.SetW2(theW2);
aData.SetL2(theL2);
aData.SetRF(theRF);
aData.SetHexMesh(theHexMesh);
bool isTRL = (theRL + theWL + theLtransL + theLthinL) > Precision::Confusion();
bool isTRR = (theRR + theWR + theLtransR + theLthinR) > Precision::Confusion();
bool isTRI = (theRI + theWI + theLtransI + theLthinI) > Precision::Confusion();
//Compute the resulting value
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("TShape driver failed");
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
// BEGIN of fillet
TopoDS_Shape aShapeShape = aShape->GetValue();
TopTools_IndexedMapOfShape anEdgesIndices;
TopExp::MapShapes(aShapeShape, anEdgesIndices);
// Common edges on external cylinders
Handle(GEOM_Object) box_e;
if (theHexMesh) {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(theR2+theW2, theR2+theW2, theR1+theW1);
}
else {
box_e = my3DPrimOperations->MakeBoxDXDYDZ(2*(theR2+theW2), 2*(theR2+theW2), theR1+theW1);
}
box_e->GetLastFunction()->SetDescription("");
box_e = myTransformOperations->TranslateDXDYDZ(box_e, -(theR2+theW2), -(theR2+theW2), 0);
box_e->GetLastFunction()->SetDescription("");
Handle(TColStd_HSequenceOfInteger) edges_e =
myShapesOperations->GetShapesOnBoxIDs(box_e, aShape, TopAbs_EDGE, GEOMAlgo_ST_IN);
box_e->GetLastFunction()->SetDescription("");
if (edges_e.IsNull() || edges_e->Length() == 0) {
SetErrorCode("External edges not found");
return NULL;
}
int nbEdgesInFillet = 0;
std::list<int> theEdges;
for (int i=1; i<=edges_e->Length();i++) {
int edgeID = edges_e->Value(i);
TopoDS_Shape theEdge = anEdgesIndices.FindKey(edgeID);
TopExp_Explorer Ex(theEdge,TopAbs_VERTEX);
while (Ex.More()) {
gp_Pnt aPt = BRep_Tool::Pnt(TopoDS::Vertex(Ex.Current()));
if (Abs(aPt.Z() - (theR1+theW1)) <= Precision::Confusion()) {
nbEdgesInFillet ++;
theEdges.push_back(edgeID);
}
Ex.Next();
}
if (theHexMesh && nbEdgesInFillet == 1)
break;
}
Handle(GEOM_Object) aFillet;
try {
aFillet = myLocalOperations->MakeFilletEdges(aShape, theRF, theEdges);
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
if (aFillet.IsNull()) {
SetErrorCode("Fillet can not be computed on the given shape with the given parameters");
return NULL;
}
aFillet->GetLastFunction()->SetDescription("");
TopoDS_Shape aFilletShape = aFillet->GetValue();
aFunction->SetValue(aFilletShape);
// END of fillet
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - BEGIN (3)
// the following block, when enabled, leads to partitioning problems
#if 0
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - END (3)
// BEGIN: Limit tolerances (debug)
Handle(GEOM_Object) aCorr1 = myHealingOperations->LimitTolerance(aShape, 1e-07);
TopoDS_Shape aCorr1Shape = aCorr1->GetValue();
aShape->GetLastFunction()->SetValue(aCorr1Shape);
aCorr1->GetLastFunction()->SetDescription("");
// END: Limit tolerances (debug)
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - BEGIN (4)
#endif
// VSR: debug issues 0021568 and 0021550 (15/05/2012) - END (4)
if (theHexMesh) {
if (!MakePipeTShapePartition(aShape, theR1, theW1, theL1, theR2, theW2, theL2, 0, 0, theRF, false))
return NULL;
if (!MakePipeTShapeMirrorAndGlue(aShape, theR1, theW1, theL1, theR2, theW2, theL2))
return NULL;
}
// Add thickness reduction elements
// at the three extremities: Left, Right and Incident
try {
#if OCC_VERSION_LARGE > 0x06010000
OCC_CATCH_SIGNALS;
#endif
if (isTRL || isTRR || isTRI) {
TopoDS_Shape aResShape =
MakePipeTShapeThicknessReduction(aShape->GetValue(), theR1, theW1, theL1, theR2, theW2, theL2,
theRL, theWL, theLtransL, theLthinL,
theRR, theWR, theLtransR, theLthinR,
theRI, theWI, theLtransI, theLthinI,
!theHexMesh);
aFunction->SetValue(aResShape);
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
// Set Position
gp_Trsf aTrsf = GetPositionTrsf(theL1, theL2, theP1, theP2, theP3);
BRepBuilderAPI_Transform aTransformation (aShape->GetValue(), aTrsf, Standard_False);
TopoDS_Shape aTrsf_Shape = aTransformation.Shape();
aFunction->SetValue(aTrsf_Shape);
Handle(TColStd_HSequenceOfTransient) aSeq = new TColStd_HSequenceOfTransient;
aSeq->Append(aShape);
if (theHexMesh) {
// Get the groups
try {
if (!MakeGroups(aShape, TSHAPE_FILLET, theR1, theW1, theL1, theR2, theW2, theL2,
0., 0., theRF, aSeq, aTrsf))
return NULL;
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
}
//Make a Python command
TCollection_AsciiString anEntry, aListRes("[");
// Iterate over the sequence aSeq
Standard_Integer aNbGroups = aSeq->Length();
Standard_Integer i = 1;
for (; i <= aNbGroups; i++) {
Handle(Standard_Transient) anItem = aSeq->Value(i);
if (anItem.IsNull()) continue;
Handle(GEOM_Object) aGroup = Handle(GEOM_Object)::DownCast(anItem);
if (aGroup.IsNull()) continue;
//Make a Python command
TDF_Tool::Entry(aGroup->GetEntry(), anEntry);
aListRes += anEntry + ", ";
}
aListRes.Trunc(aListRes.Length() - 2);
GEOM::TPythonDump pd (aFunction);
pd << aListRes.ToCString() << "] = geompy.MakePipeTShapeFillet("
<< theR1 << ", " << theW1 << ", " << theL1 << ", "
<< theR2 << ", " << theW2 << ", " << theL2 << ", "
<< theRF << ", " << theHexMesh << ", "
<< theP1 << ", " << theP2 << ", " << theP3;
// thickness reduction
if (isTRL)
pd << ", theRL=" << theRL << ", theWL=" << theWL
<< ", theLtransL=" << theLtransL << ", theLthinL=" << theLthinL;
if (isTRR)
pd << ", theRR=" << theRR << ", theWR=" << theWR
<< ", theLtransR=" << theLtransR << ", theLthinR=" << theLthinR;
if (isTRI)
pd << ", theRI=" << theRI << ", theWI=" << theWI
<< ", theLtransI=" << theLtransI << ", theLthinI=" << theLthinI;
pd << ")";
SetErrorCode(OK);
return aSeq;
}
//=============================================================================
/*!
* This function allows to create a disk already divided into blocks. It can be
* used to create divided pipes for later meshing in hexaedra.
* \param theR Radius of the disk
* \param theRatio Relative size of the central square diagonal against the disk diameter
* \param theOrientation Plane on which the disk will be built
* \param thePattern The division pattern of the disk (hexagon or square in the center)
* \return New GEOM_Object, containing the created shape.
*/
//=============================================================================
Handle(GEOM_Object) GEOMImpl_IAdvancedOperations::MakeDividedDisk (double theR, double theRatio,
int theOrientation, int thePattern)
{
SetErrorCode(KO);
if (theOrientation != 1 &&
theOrientation != 2 &&
theOrientation != 3)
{
SetErrorCode("theOrientation must be 1(=OXY), 2(=OYZ) or 3(=OZX)");
return NULL;
}
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_DIVIDEDDISK);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_DividedDiskDriver::GetID(), DIVIDEDDISK_R_RATIO);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_DividedDiskDriver::GetID()) return NULL;
GEOMImpl_IDividedDisk aData (aFunction);
aData.SetR(theR);
aData.SetRatio(theRatio);
aData.SetOrientation(theOrientation);
aData.SetType(thePattern);
//Compute the resulting value
try {
#if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("DividedDisk driver failed");
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
std::string aPatternStr;
switch(thePattern)
{
case 0:
aPatternStr = "GEOM.SQUARE";
break;
case 1:
aPatternStr = "GEOM.HEXAGON";
break;
}
//Make a Python command
GEOM::TPythonDump(aFunction) << aShape << " = geompy.MakeDividedDisk(" << theR << ", " << theOrientation << ", " << aPatternStr.c_str() << ")";
SetErrorCode(OK);
return aShape;
}
//=============================================================================
/*!
* This function allows to create a disk already divided into blocks. It can be
* used to create divided pipes for later meshing in hexaedra.
* \param theR Radius of the disk
* \param theRatio Relative size of the central square diagonal against the disk diameter
* \return New GEOM_Object, containing the created shape.
*/
//=============================================================================
Handle(GEOM_Object) GEOMImpl_IAdvancedOperations::MakeDividedDiskPntVecR (Handle(GEOM_Object) thePnt,
Handle(GEOM_Object) theVec,
double theR,
double theRatio,
int thePattern)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_DIVIDEDDISK);
//Add a new shape function with parameters
Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_DividedDiskDriver::GetID(), DIVIDEDDISK_R_VECTOR_PNT);
if (aFunction.IsNull()) return NULL;
//Check if the function is set correctly
if (aFunction->GetDriverGUID() != GEOMImpl_DividedDiskDriver::GetID()) return NULL;
GEOMImpl_IDividedDisk aData (aFunction);
Handle(GEOM_Function) aRefPnt = thePnt->GetLastFunction();
Handle(GEOM_Function) aRefVec = theVec->GetLastFunction();
if (aRefPnt.IsNull() || aRefVec.IsNull()) return NULL;
aData.SetCenter(aRefPnt);
aData.SetVector(aRefVec);
aData.SetR(theR);
aData.SetRatio(theRatio);
aData.SetType(thePattern);
//Compute the resulting value
try {
#if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
OCC_CATCH_SIGNALS;
#endif
if (!GetSolver()->ComputeFunction(aFunction)) {
SetErrorCode("DividedDisk driver failed");
return NULL;
}
}
catch (Standard_Failure) {
Handle(Standard_Failure) aFail = Standard_Failure::Caught();
SetErrorCode(aFail->GetMessageString());
return NULL;
}
std::string aPatternStr;
switch(thePattern)
{
case 0:
aPatternStr = "GEOM.SQUARE";
break;
case 1:
aPatternStr = "GEOM.HEXAGON";
break;
}
//Make a Python command
GEOM::TPythonDump(aFunction) << aShape << " = geompy.MakeDividedDiskPntVecR(" << thePnt << ", " << theVec << ", " << theR << ", " << aPatternStr.c_str() << ")";
SetErrorCode(OK);
return aShape;
}
//=============================================================================
/*!
* Builds a cylinder prepared for hexa meshes
* \param theR Radius of the cylinder
* \param theH Height of the cylinder
* \return New GEOM_Object, containing the created shape.
*/
//=============================================================================
Handle(GEOM_Object) GEOMImpl_IAdvancedOperations::MakeDividedCylinder (double theR,
double theH,
int thePattern)
{
SetErrorCode(KO);
//Add a new object
Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_DIVIDEDCYLINDER);
Handle(GEOM_Object) aBaseShape = MakeDividedDisk(theR, 67.0, 1, thePattern);
aBaseShape->GetLastFunction()->SetDescription(""); // Erase dump of MakeDividedDisk
aShape = my3DPrimOperations->MakePrismDXDYDZ(aBaseShape,0.0,0.0,theH, -1.0);
Handle(GEOM_Function) aFunction = aShape->GetLastFunction();
aFunction->SetDescription(""); // Erase dump of MakePrismDXDYDZ
aShape->SetType(GEOM_DIVIDEDCYLINDER);
std::string aPatternStr;
switch(thePattern)
{
case 0:
aPatternStr = "GEOM.SQUARE";
break;
case 1:
aPatternStr = "GEOM.HEXAGON";
break;
}
//Make a Python command
GEOM::TPythonDump(aFunction) << aShape << " = geompy.MakeDividedCylinder(" << theR << ", " << theH << ", " << aPatternStr.c_str() << ")";
SetErrorCode(OK);
return aShape;
}
/*@@ insert new functions before this line @@ do not remove this line @@ do not remove this line @@*/
Reference in New Issue View Git Blame Copy Permalink