X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FAdvancedEngine%2FGEOMImpl_IAdvancedOperations.cxx;fp=src%2FAdvancedEngine%2FGEOMImpl_IAdvancedOperations.cxx;h=0000000000000000000000000000000000000000;hb=aeddae77192d933be1aee8855d8f7c525878cad9;hp=59c4d95f8c659e6eda538167cd598a8098ec6b2e;hpb=aaea242262db5525ffcb4b90b65515048c0c00b5;p=modules%2Fgeom.git diff --git a/src/AdvancedEngine/GEOMImpl_IAdvancedOperations.cxx b/src/AdvancedEngine/GEOMImpl_IAdvancedOperations.cxx deleted file mode 100644 index 59c4d95f8..000000000 --- a/src/AdvancedEngine/GEOMImpl_IAdvancedOperations.cxx +++ /dev/null @@ -1,3508 +0,0 @@ -// Copyright (C) 2007-2014 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, or (at your option) any later version. -// -// This library is distributed in the hope that it will be useful, -// but WITHOUT ANY WARRANTY; without even the implied warranty of -// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU -// 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 - -#include -#include -#include - -#include "GEOM_Function.hxx" -#include "GEOM_PythonDump.hxx" -#include "GEOMUtils.hxx" -#include "GEOMAlgo_Splitter.hxx" -#include "GEOMAlgo_FinderShapeOn1.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_IGroupOperations.hxx" - -#include "GEOMImpl_GlueDriver.hxx" -#include "GEOMImpl_PipeTShapeDriver.hxx" -#include "GEOMImpl_IPipeTShape.hxx" -#include "GEOMImpl_DividedDiskDriver.hxx" -#include "GEOMImpl_IDividedDisk.hxx" -#include -#include -/*@@ insert new functions before this line @@ do not remove this line @@ do not remove this line @@*/ - -#include -#include -#include -#include -#include - -#include -#include -#include -#include -#include -#include -#include - -#include -#include - -#include -#include -#include -#include -#include -#include -#include -#include - -#include -#include -#include -#include -#include -#include - -#include - -#include - -#include "AdvancedEngine_Types.hxx" - -#include -#include -#include -#include // 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()); - myGroupOperations = new GEOMImpl_IGroupOperations(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; - delete myGroupOperations; -} - -//============================================================================= -/*! - * 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, Standard_False); // 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 aCompoundOfFacesList; - - for (int i=0 ; i<= aCommonFaces->Length()-4 ; i+=4) { - std::list 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; -} - -//============================================================================= -/*! - * Return faces that are laying on surface. - */ -//============================================================================= -bool GEOMImpl_IAdvancedOperations::GetFacesOnSurf - (const TopoDS_Shape &theShape, - const Handle_Geom_Surface& theSurface, - const Standard_Real theTolerance, - TopTools_ListOfShape &theFaces) -{ - GEOMAlgo_FinderShapeOn1 aFinder; - - aFinder.SetShape(theShape); - aFinder.SetTolerance(theTolerance); - aFinder.SetSurface(theSurface); - aFinder.SetShapeType(TopAbs_FACE); - aFinder.SetState(GEOMAlgo_ST_ON); - - // Sets the minimal number of inner points for the faces that do not have own - // inner points at all (for e.g. rectangular planar faces have just 2 triangles). - // Default value=3 - aFinder.SetNbPntsMin(3); - // Sets the maximal number of inner points for edges or faces. - // It is usefull for the cases when this number is very big (e.g =2000) to improve - // the performance. If this value =0, all inner points will be taken into account. - // Default value=0 - aFinder.SetNbPntsMax(100); - aFinder.Perform(); - - // Interprete results - Standard_Integer iErr = aFinder.ErrorStatus(); - // the detailed description of error codes is in GEOMAlgo_FinderShapeOn1.cxx - if (iErr) { - MESSAGE(" iErr : " << iErr); - TCollection_AsciiString aMsg (" iErr : "); - aMsg += TCollection_AsciiString(iErr); - SetErrorCode(aMsg); - return false; - } - Standard_Integer iWrn = aFinder.WarningStatus(); - // the detailed description of warning codes is in GEOMAlgo_FinderShapeOn1.cxx - if (iWrn) { - MESSAGE(" *** iWrn : " << iWrn); - } - - const TopTools_ListOfShape &aListRes = aFinder.Shapes(); // the result - TopTools_ListIteratorOfListOfShape anIter (aListRes); - - for (; anIter.More(); anIter.Next()) { - theFaces.Append(anIter.Value()); - } - - return true; -} - -//============================================================================= -/*! - * Creates and returns conical face. - */ -//============================================================================= -TopoDS_Shape GEOMImpl_IAdvancedOperations::MakeConicalFace - (const gp_Ax2 &theAxis, - const double theRadius, - const double theRadiusThin, - const double theHeight, - const gp_Trsf &theTrsf) -{ - BRepPrimAPI_MakeCone aMkCone (theAxis, theRadius, theRadiusThin, theHeight); - TopoDS_Shape aResult; - - aMkCone.Build(); - if (aMkCone.IsDone()) { - TopExp_Explorer anExp(aMkCone.Shape(), TopAbs_FACE); - - for (; anExp.More(); anExp.Next()) { - TopoDS_Face aFace = TopoDS::Face(anExp.Current()); - - if (aFace.IsNull() == Standard_False) { - BRepAdaptor_Surface anAdaptor(aFace, Standard_False); - - if (anAdaptor.GetType() == GeomAbs_Cone) { - // This is a conical face. Transform and return it. - BRepBuilderAPI_Transform aTransf(aFace, theTrsf, Standard_False); - - aResult = aTransf.Shape(); - break; - } - } - } - } - - return aResult; -} - -//============================================================================= -/*! - * Generate the internal group of a Pipe T-Shape - */ -//============================================================================= -bool GEOMImpl_IAdvancedOperations::MakeInternalGroup - (const Handle(GEOM_Object) &theShape, - const double theR1, const double theLen1, - const double theR2, const double theLen2, - const double theRL, double theTransLenL, - const double theRR, double theTransLenR, - const double theRI, double theTransLenI, - const Handle(TColStd_HSequenceOfTransient) &theSeq, - const gp_Trsf &theTrsf) -{ - SetErrorCode(KO); - - if (theShape.IsNull()) { - return false; - } - - TopoDS_Shape aShape = theShape->GetValue(); - - if (aShape.IsNull()) { - SetErrorCode("Shape is not defined"); - return false; - } - - // Compute tolerance - Standard_Real aMaxTol = -RealLast(); - TopExp_Explorer anExp(aShape, TopAbs_VERTEX); - - for (; anExp.More(); anExp.Next()) { - TopoDS_Vertex aVertex = TopoDS::Vertex(anExp.Current()); - - if (aVertex.IsNull() == Standard_False) { - const Standard_Real aTol = BRep_Tool::Tolerance(aVertex); - - if (aTol > aMaxTol) { - aMaxTol = aTol; - } - } - } - - // Construct internal surfaces. - Standard_Integer i = 0; - const Standard_Integer aMaxNbSurf = 5; - Handle(Geom_Surface) aSurface[aMaxNbSurf]; - TopTools_ListOfShape aConicalFaces; - Standard_Real aTolConf = Precision::Confusion(); - - // 1. Construct the internal surface of main pipe. - gp_Ax2 anAxis1 (gp::Origin(), gp::DX(), gp::DZ()); - gp_Ax2 anAxis2 (gp::Origin(), gp::DZ(), gp::DX()); - - aSurface[i++] = new Geom_CylindricalSurface(anAxis1, theR1); - - // 2. Construct the internal surface of incident pipe. - aSurface[i++] = new Geom_CylindricalSurface(anAxis2, theR2); - - // 3. Construct the internal surface of left reduction pipe. - if (theRL > aTolConf) { - aSurface[i++] = new Geom_CylindricalSurface(anAxis1, theRL); - - if (theTransLenL > aTolConf) { - // 3.1. Construct the internal surface of left transition pipe. - gp_Pnt aPLeft (-theLen1, 0., 0.); - gp_Ax2 anAxisLeft (aPLeft, -gp::DX(), gp::DZ()); - TopoDS_Shape aConeLeft = - MakeConicalFace(anAxisLeft, theR1, theRL, theTransLenL, theTrsf); - - if (aConeLeft.IsNull() == Standard_False) { - aConicalFaces.Append(aConeLeft); - } - } - } - - // 4. Construct the internal surface of right reduction pipe. - if (theRR > aTolConf) { - // There is no need to construct another cylinder of the same radius. Skip it. - if (Abs(theRR - theRL) > aTolConf) { - aSurface[i++] = new Geom_CylindricalSurface(anAxis1, theRR); - } - - if (theTransLenL > aTolConf) { - // 4.1. Construct the internal surface of right transition pipe. - gp_Pnt aPRight (theLen1, 0., 0.); - gp_Ax2 anAxisRight (aPRight, gp::DX(), gp::DZ()); - TopoDS_Shape aConeRight = - MakeConicalFace(anAxisRight, theR1, theRR, theTransLenR, theTrsf); - - if (aConeRight.IsNull() == Standard_False) { - aConicalFaces.Append(aConeRight); - } - } - } - - // 5. Construct the internal surface of incident reduction pipe. - if (theRI > aTolConf) { - aSurface[i++] = new Geom_CylindricalSurface(anAxis2, theRI); - - if (theTransLenI > aTolConf) { - // 5.1. Construct the internal surface of incident transition pipe. - gp_Pnt aPInci (0., 0., theLen2); - gp_Ax2 anAxisInci (aPInci, gp::DZ(), gp::DX()); - TopoDS_Shape aConeInci = - MakeConicalFace(anAxisInci, theR2, theRI, theTransLenI, theTrsf); - - if (aConeInci.IsNull() == Standard_False) { - aConicalFaces.Append(aConeInci); - } - } - } - - // Get faces that are laying on cylindrical surfaces. - TopTools_ListOfShape aFaces; - gp_Trsf anInvTrsf = theTrsf.Inverted(); - - for (i = 0; i < aMaxNbSurf; i++) { - if (aSurface[i].IsNull()) { - break; - } - - aSurface[i]->Transform(theTrsf); - - TopTools_ListOfShape aLocalFaces; - - if (!GetFacesOnSurf(aShape, aSurface[i], aMaxTol, aLocalFaces)) { - return false; - } - - if (i < 2) { - // Check if the result contains outer cylinders. - // It is required for main and incident pipes. - TopTools_ListIteratorOfListOfShape anIter(aLocalFaces); - - while (anIter.More()) { - TopExp_Explorer anExp(anIter.Value(), TopAbs_VERTEX); - Standard_Boolean isInside = Standard_False; - - // Get a vertex from this shape - if (anExp.More()) { - TopoDS_Vertex aVtx = TopoDS::Vertex(anExp.Current()); - - if (aVtx.IsNull() == Standard_False) { - gp_Pnt aPnt = BRep_Tool::Pnt(aVtx); - - aPnt.Transform(anInvTrsf); - - if (i == 0) { - // Check if the point is inside the main pipe. - isInside = (Abs(aPnt.X()) <= theLen1); - } else { // i == 1 - // Check if the point is inside the incident pipe. - isInside = (aPnt.Z() <= theLen2); - } - } - } - - if (isInside) { - // Keep this face. - anIter.Next(); - } else { - // Remove this face. - aLocalFaces.Remove(anIter); - } - } - } - - aFaces.Append(aLocalFaces); - } - - // Get faces that are laying on conical faces. - if (aConicalFaces.IsEmpty() == Standard_False) { - Handle(GEOM_Object) aCone = - GetEngine()->AddObject(GetDocID(), GEOM_TSHAPE); - Handle(GEOM_Function) aFunction = - aCone->AddFunction(GEOMImpl_PipeTShapeDriver::GetID(), TSHAPE_BASIC); - TopTools_ListIteratorOfListOfShape aFIter(aConicalFaces); - Handle(GEOM_Object) aConeFromShape; - - for (; aFIter.More(); aFIter.Next()) { - aFunction->SetValue(aFIter.Value()); - aConeFromShape = myShapesOperations->GetInPlace(theShape, aCone); - - if (aConeFromShape.IsNull() == Standard_False) { - aConeFromShape->GetLastFunction()->SetDescription(""); - TopoDS_Shape aConeFaces = aConeFromShape->GetValue(); - TopExp_Explorer anExp(aConeFaces, TopAbs_FACE); - - for (; anExp.More(); anExp.Next()) { - TopoDS_Face aConeFace = TopoDS::Face(anExp.Current()); - - if (aConeFace.IsNull() == Standard_False) { - aFaces.Append(aConeFace); - } - } - } - } - } - - // Create a group of internal faces. - if (aFaces.IsEmpty() == Standard_False) { - Handle(GEOM_Object) aGroup = myGroupOperations->CreateGroup(theShape, TopAbs_FACE); - - if (aGroup.IsNull() == Standard_False) { - aGroup->GetLastFunction()->SetDescription(""); - aGroup->SetName("INTERNAL_FACES"); - - TopTools_IndexedMapOfShape anIndices; - Handle(TColStd_HSequenceOfInteger) aSeqIDs = new TColStd_HSequenceOfInteger; - - TopExp::MapShapes(aShape, anIndices); - - TopTools_ListIteratorOfListOfShape anIter(aFaces); - - for (; anIter.More(); anIter.Next()) { - const TopoDS_Shape &aFace = anIter.Value(); - const Standard_Integer anIndex = anIndices.FindIndex(aFace); - - if (anIndex > 0) { - aSeqIDs->Append(anIndex); - } - } - - myGroupOperations->UnionIDs(aGroup, aSeqIDs); - aGroup->GetLastFunction()->SetDescription(""); - theSeq->Append(aGroup); - } - } - - 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, arete_intersect_ext; - 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; - } - - double d1min = theR2+theW2, d2min=theR2+theW2; - 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); - if (Abs(aP.X()) <= Precision::Confusion()) { - if (Abs(aP.Y()) < d1min) { - vi1 = v; - d1min = Abs(aP.Y()); - } - } else if (Abs(aP.Y()) <= Precision::Confusion()) { - if (Abs(aP.X()) < d2min) { - vi2 = v; - d2min = Abs(aP.X()); - } - } - } - if (vi1.IsNull() || vi2.IsNull()) { - SetErrorCode("Cannot find internal intersection vertices"); - return false; - } - - std::list theShapes; - - if (isNormal) { - Handle(GEOM_Object) ve1, ve2; - TopoDS_Vertex vertex1, vertex2; - - 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(""); - - // 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; - } - - double d1max = 0, d2max = 0; - 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); - if (Abs(aP.X()) <= Precision::Confusion()) { - if (Abs(aP.Y()) > d1max) { - ve1 = v; - vertex1 = aVertex; - d1max = Abs(aP.Y()); - } - } else if (Abs(aP.Y()) <= Precision::Confusion()) { - if (Abs(aP.X()) > d2max) { - ve2 = v; - vertex2 = aVertex; - d2max = Abs(aP.X()); - } - } - } - if (ve1.IsNull() || ve2.IsNull()) { - SetErrorCode("Cannot find external intersection vertices"); - return false; - } - Handle(GEOM_Object) edge_e1, edge_e2; - - // 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; - } - ShapeAnalysis_Edge sae; - for (int i=1; i<=edges_e->Length();i++) { - Handle(GEOM_Object) anObj = Handle(GEOM_Object)::DownCast(edges_e->Value(i)); - anObj->GetLastFunction()->SetDescription(""); - TopoDS_Edge anEdge = TopoDS::Edge(anObj->GetValue()); - if ( !anEdge.IsNull() && - (sae.FirstVertex(anEdge).IsSame(vertex1) || sae.LastVertex(anEdge).IsSame(vertex1)) && - (sae.FirstVertex(anEdge).IsSame(vertex2) || sae.LastVertex(anEdge).IsSame(vertex2))) { - arete_intersect_ext = anObj; - } - } - - 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 edge_e_elist; - edge_e_elist.push_back(arete_intersect_int); - edge_e_elist.push_back(edge_e1); - edge_e_elist.push_back(arete_intersect_ext); - 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 LX; - std::vector 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 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 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, Standard_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 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 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.AddArgument(aReduction); - if (isThinPart) - PS.AddArgument(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); - - try { - if (theHexMesh) { - // Get the groups - if (!MakeGroups(aShape, TSHAPE_BASIC, theR1, theW1, theL1, theR2, theW2, theL2, - 0., 0., 0., aSeq, gp_Trsf())) - return NULL; - } - - // Get internal group. - if (!MakeInternalGroup(aShape, theR1, theL1, theR2, theL2, theRL, theLtransL, - theRR, theLtransR, theRI, theLtransI, - 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); - - try { - if (theHexMesh) { - // Get the groups - if (!MakeGroups(aShape,TSHAPE_BASIC, theR1, theW1, theL1, theR2, theW2, theL2, - 0., 0., 0., aSeq, aTrsf)) { - return NULL; - } - } - - // Get internal group. - if (!MakeInternalGroup(aShape, theR1, theL1, theR2, theL2, theRL, theLtransL, - theRR, theLtransR, theRI, theLtransI, - 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 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); - - try { - if (theHexMesh) { - // Get the groups - if (!MakeGroups(aShape, TSHAPE_CHAMFER, theR1, theW1, theL1, theR2, theW2, theL2, - theH, theW, 0., aSeq, gp_Trsf())) - return NULL; - } - - // Get internal group. - if (!MakeInternalGroup(aShape, theR1, theL1, theR2, theL2, theRL, theLtransL, - theRR, theLtransR, theRI, theLtransI, - 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 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); - - try { - if (theHexMesh) { - // Get the groups - if (!MakeGroups(aShape, TSHAPE_CHAMFER, theR1, theW1, theL1, theR2, theW2, theL2, - theH, theW, 0., aSeq, aTrsf)) - return NULL; - } - - // Get internal group. - if (!MakeInternalGroup(aShape, theR1, theL1, theR2, theL2, theRL, theLtransL, - theRR, theLtransR, theRI, theLtransI, - 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 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); - - try { - if (theHexMesh) { - // Get the groups - if (!MakeGroups(aShape, TSHAPE_FILLET, theR1, theW1, theL1, theR2, theW2, theL2, - 0., 0., theRF, aSeq, gp_Trsf())) - return NULL; - } - - // Get internal group. - if (!MakeInternalGroup(aShape, theR1, theL1, theR2, theL2, theRL, theLtransL, - theRR, theLtransR, theRI, theLtransI, - 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 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); - - try { - if (theHexMesh) { - // Get the groups - if (!MakeGroups(aShape, TSHAPE_FILLET, theR1, theW1, theL1, theR2, theW2, theL2, - 0., 0., theRF, aSeq, aTrsf)) - return NULL; - } - - // Get internal group. - if (!MakeInternalGroup(aShape, theR1, theL1, theR2, theL2, theRL, theLtransL, - theRR, theLtransR, theRI, theLtransI, - 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; -} -//============================================================================= -/*! - * Create a smoothing surface from a set of points - * \param thelPoints list of points or compounds of points - * \param theNbMax maximum number of Bezier pieces in the resulting surface. - * \param theDegMax maximum degree of the resulting BSpline surface - * \param theDMax specifies maximum value of the GeomPlate_PlateG0Criterion criterion. - * \return New GEOM_Object, containing the created shape. - */ -//============================================================================= -Handle(GEOM_Object) GEOMImpl_IAdvancedOperations::MakeSmoothingSurface (std::list thelPoints, - int theNbMax, - int theDegMax, - double theDMax) -{ - SetErrorCode(KO); - - //Add a new object - Handle(GEOM_Object) aShape = GetEngine()->AddObject(GetDocID(), GEOM_SMOOTHINGSURFACE); - - //Add a new shape function with parameters - Handle(GEOM_Function) aFunction = aShape->AddFunction(GEOMImpl_SmoothingSurfaceDriver::GetID(), SMOOTHINGSURFACE_LPOINTS); - if (aFunction.IsNull()) return NULL; - - //Check if the function is set correctly - if (aFunction->GetDriverGUID() != GEOMImpl_SmoothingSurfaceDriver::GetID()) return NULL; - - GEOMImpl_ISmoothingSurface aData (aFunction); - - int aLen = thelPoints.size(); - aData.SetLength(aLen); - int ind = 1; - std::list::iterator it = thelPoints.begin(); - for (; it != thelPoints.end(); it++, ind++) { - Handle(GEOM_Function) aRefObj = (*it)->GetLastFunction(); - if (aRefObj.IsNull()) { - SetErrorCode("NULL point or compound for bSplineFaceShape"); - return NULL; - } - aData.SetPntOrComp(ind, aRefObj); - } - - aData.SetNbMax(theNbMax); - aData.SetDegMax(theDegMax); - aData.SetDMax(theDMax); - - //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("SmoothingSurface driver failed"); - return NULL; - } - } - catch (Standard_Failure) { - Handle(Standard_Failure) aFail = Standard_Failure::Caught(); - SetErrorCode(aFail->GetMessageString()); - return NULL; - } - - //Make a Python command - GEOM::TPythonDump pd (aFunction); - pd << aShape << " = geompy.MakeSmoothingSurface(["; - it = thelPoints.begin(); - pd << (*it++); - while (it != thelPoints.end()) { - pd << ", " << (*it++); - } - pd << "], " - << theNbMax << ", " - << theDegMax << ", " - << theDMax <<")"; - - SetErrorCode(OK); - - return aShape; -} -/*@@ insert new functions before this line @@ do not remove this line @@ do not remove this line @@*/