]> SALOME platform Git repositories - modules/shaper.git/blob - src/GeomAlgoAPI/GeomAlgoAPI_ConeSegment.cpp
Salome HOME
0535a36e3d4f62c65c755f6bff4c8e3adfb7ce3f
[modules/shaper.git] / src / GeomAlgoAPI / GeomAlgoAPI_ConeSegment.cpp
1 // Copyright (C) 2014-2024  CEA, EDF
2 //
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License, or (at your option) any later version.
7 //
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
11 // Lesser General Public License for more details.
12 //
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
16 //
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
18 //
19
20 #include <GeomAlgoAPI_ConeSegment.h>
21
22 #include <gp_Dir.hxx>
23 #include <gp_Ax1.hxx>
24 #include <gp_Pnt.hxx>
25 #include <TopoDS_Edge.hxx>
26
27 #include <BRepPrimAPI_MakeRevol.hxx>
28 #include <BRepBuilderAPI_MakeEdge.hxx>
29 #include <BRepBuilderAPI_MakeWire.hxx>
30 #include <BRepBuilderAPI_MakeFace.hxx>
31 #include <Precision.hxx>
32
33 #include <math.h>
34
35 //=================================================================================================
36
37 GeomAlgoAPI_ConeSegment::GeomAlgoAPI_ConeSegment(const double theRMin1,
38                                                  const double theRMax1,
39                                                  const double theRMin2,
40                                                  const double theRMax2,
41                                                  const double theZ,
42                                                  const double theStartPhi,
43                                                  const double theDeltaPhi)
44 {
45   myRMin1 = theRMin1;
46   myRMax1 = theRMax1;
47   myRMin2 = theRMin2;
48   myRMax2 = theRMax2;
49   myZ = theZ;
50   myStartPhi = theStartPhi;
51   myDeltaPhi = theDeltaPhi;
52 }
53
54 //=================================================================================================
55 bool GeomAlgoAPI_ConeSegment::check()
56 {
57   if (myRMin1 < 0.)
58   {
59     myError = "Cone Segment builder :: rmin1 is negative.";
60     return false;
61   }
62   if (myRMin2 < 0.)
63   {
64     myError = "Cone Segment builder :: rmin2 is negative.";
65     return false;
66   }
67   if ((myRMax1-myRMin1) < Precision::Confusion())
68   {
69     myError = "Cone Segment builder :: rmin1 is larger than rmax1.";
70     return false;
71   }
72   if ((myRMax2-myRMin2) < Precision::Confusion())
73   {
74     myError = "Cone Segment builder :: rmin2 is larger than rmax2.";
75     return false;
76   }
77   if (myZ < Precision::Confusion())
78   {
79     myError = "Cone Segment builder :: z is negative or null.";
80     return false;
81   }
82   if (myDeltaPhi < Precision::Angular() * 180./M_PI)
83   {
84     myError = "Cone Segment builder :: deltaphi is negative or null.";
85     return false;
86   }
87   if (myDeltaPhi > 360)
88   {
89     myError = "Cone Segment builder :: deltaphi is larger than 360 degree.";
90     return false;
91   }
92   return true;
93 }
94
95 //=================================================================================================
96 void GeomAlgoAPI_ConeSegment::build()
97 {
98   myCreatedFaces.clear();
99
100   const double aStartPhiRad = myStartPhi * M_PI/180.;
101   BRepBuilderAPI_MakeWire aWireBuilder;
102
103   // Define the section from the 4 vertices
104   gp_Pnt aPointOuterBase(myRMax1 * cos(aStartPhiRad), myRMax1 * sin(aStartPhiRad), -myZ/2.);
105   gp_Pnt aPointInnerBase(myRMin1 * cos(aStartPhiRad), myRMin1 * sin(aStartPhiRad), -myZ/2.);
106   gp_Pnt aPointOuterTop(myRMax2 * cos(aStartPhiRad), myRMax2 * sin(aStartPhiRad), myZ/2.);
107   gp_Pnt aPointInnerTop(myRMin2 * cos(aStartPhiRad), myRMin2 * sin(aStartPhiRad), myZ/2.);
108
109   if ((myRMax1 - myRMin1) >= Precision::Confusion()){
110     BRepBuilderAPI_MakeEdge anEdgeBuilderBase(aPointOuterBase, aPointInnerBase);
111     anEdgeBuilderBase.Build();
112     aWireBuilder.Add(anEdgeBuilderBase.Edge());
113   }
114
115   BRepBuilderAPI_MakeEdge anEdgeBuilderOuter(aPointOuterBase, aPointOuterTop);
116   anEdgeBuilderOuter.Build();
117   aWireBuilder.Add(anEdgeBuilderOuter.Edge());
118   if ((myRMax2 - myRMin2) >= Precision::Confusion()){
119     BRepBuilderAPI_MakeEdge anEdgeBuilderTop(aPointOuterTop, aPointInnerTop);
120     anEdgeBuilderTop.Build();
121     aWireBuilder.Add(anEdgeBuilderTop.Edge());
122   }
123
124   BRepBuilderAPI_MakeEdge anEdgeBuilderInner(aPointInnerTop, aPointInnerBase);
125   anEdgeBuilderInner.Build();
126   aWireBuilder.Add(anEdgeBuilderInner.Edge());
127
128   aWireBuilder.Build();
129   BRepBuilderAPI_MakeFace aFaceBuilder(aWireBuilder.Wire());
130   aFaceBuilder.Build();
131
132   if (!aFaceBuilder.IsDone()){
133     myError = "Cone Segment builder :: section is not valid";
134     return;
135   }
136
137   // Perform a revolution based on the section to build the solid
138   gp_Dir aZDir(0., 0., 1.);
139   gp_Pnt anOrigin(0., 0., 0.);
140   gp_Ax1 aZAxis(anOrigin, aZDir);
141   BRepPrimAPI_MakeRevol* aRevolBuilder =
142     new BRepPrimAPI_MakeRevol(aFaceBuilder.Face(), aZAxis, myDeltaPhi * M_PI/180., Standard_True);
143   if(!aRevolBuilder) {
144     myError = "Cone Segment builder :: section revolution did not succeed";
145     return;
146   }
147   if(!aRevolBuilder->IsDone()) {
148     myError = "Cone Segment builder :: section revolution did not succeed";
149     return;
150   }
151
152   setImpl(aRevolBuilder);
153   setBuilderType(OCCT_BRepBuilderAPI_MakeShape);
154
155   std::shared_ptr<GeomAPI_Shape> aResultShape =
156     std::shared_ptr<GeomAPI_Shape>(new GeomAPI_Shape());
157   aResultShape->setImpl(new TopoDS_Shape(aRevolBuilder->Shape()));
158
159   // Test on the shapes
160   if (!(aResultShape).get() || aResultShape->isNull()) {
161     myError = "Cone Segment builder  :: resulting shape is null.";
162     return;
163   }
164
165   setShape(aResultShape);
166   setDone(true);
167 }