Commit | Line | Data |
---|---|---|
b311480e | 1 | // Created on: 1996-02-13 |
2 | // Created by: Jacques GOUSSARD | |
3 | // Copyright (c) 1996-1999 Matra Datavision | |
973c2be1 | 4 | // Copyright (c) 1999-2014 OPEN CASCADE SAS |
b311480e | 5 | // |
973c2be1 | 6 | // This file is part of Open CASCADE Technology software library. |
b311480e | 7 | // |
d5f74e42 | 8 | // This library is free software; you can redistribute it and/or modify it under |
9 | // the terms of the GNU Lesser General Public License version 2.1 as published | |
973c2be1 | 10 | // by the Free Software Foundation, with special exception defined in the file |
11 | // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT | |
12 | // distribution for complete text of the license and disclaimer of any warranty. | |
b311480e | 13 | // |
973c2be1 | 14 | // Alternatively, this file may be used under the terms of Open CASCADE |
15 | // commercial license or contractual agreement. | |
7fd59977 | 16 | |
7fd59977 | 17 | |
7fd59977 | 18 | #include <Bnd_Box.hxx> |
7fd59977 | 19 | #include <BRep_Builder.hxx> |
20 | #include <BRep_Tool.hxx> | |
7fd59977 | 21 | #include <BRepAlgoAPI_Cut.hxx> |
22 | #include <BRepAlgoAPI_Fuse.hxx> | |
42cf5bc1 | 23 | #include <BRepBndLib.hxx> |
24 | #include <BRepBuilderAPI.hxx> | |
25 | #include <BRepFeat.hxx> | |
26 | #include <BRepFeat_MakeRevol.hxx> | |
7fd59977 | 27 | #include <BRepLib_MakeFace.hxx> |
42cf5bc1 | 28 | #include <BRepSweep_Revol.hxx> |
7fd59977 | 29 | #include <BRepTools_Modifier.hxx> |
30 | #include <BRepTools_TrsfModification.hxx> | |
42cf5bc1 | 31 | #include <ElCLib.hxx> |
32 | #include <Geom2d_Curve.hxx> | |
33 | #include <Geom_Circle.hxx> | |
34 | #include <Geom_Curve.hxx> | |
35 | #include <Geom_Plane.hxx> | |
36 | #include <Geom_RectangularTrimmedSurface.hxx> | |
37 | #include <Geom_Surface.hxx> | |
38 | #include <gp_Ax1.hxx> | |
39 | #include <gp_Pln.hxx> | |
40 | #include <gp_Pnt.hxx> | |
41 | #include <gp_Pnt2d.hxx> | |
42 | #include <gp_Vec.hxx> | |
43 | #include <LocOpe.hxx> | |
44 | #include <LocOpe_BuildShape.hxx> | |
45 | #include <LocOpe_CSIntersector.hxx> | |
46 | #include <LocOpe_FindEdges.hxx> | |
47 | #include <LocOpe_Gluer.hxx> | |
48 | #include <LocOpe_PntFace.hxx> | |
49 | #include <LocOpe_Revol.hxx> | |
50 | #include <LocOpe_SequenceOfCirc.hxx> | |
51 | #include <Precision.hxx> | |
7fd59977 | 52 | #include <Standard_ConstructionError.hxx> |
42cf5bc1 | 53 | #include <TColgp_SequenceOfPnt.hxx> |
54 | #include <TopExp.hxx> | |
55 | #include <TopExp_Explorer.hxx> | |
7fd59977 | 56 | #include <TopoDS_Compound.hxx> |
42cf5bc1 | 57 | #include <TopoDS_Edge.hxx> |
58 | #include <TopoDS_Face.hxx> | |
59 | #include <TopoDS_Shape.hxx> | |
7fd59977 | 60 | #include <TopoDS_Shell.hxx> |
42cf5bc1 | 61 | #include <TopoDS_Solid.hxx> |
62 | #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx> | |
63 | #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx> | |
64 | #include <TopTools_ListIteratorOfListOfShape.hxx> | |
65 | #include <TopTools_MapIteratorOfMapOfShape.hxx> | |
66 | #include <TopTools_MapOfShape.hxx> | |
7fd59977 | 67 | |
0797d9d3 | 68 | #ifdef OCCT_DEBUG |
1d0a9d4d | 69 | extern Standard_Boolean BRepFeat_GettraceFEAT(); |
7fd59977 | 70 | #endif |
71 | ||
72 | static void MajMap(const TopoDS_Shape&, // base | |
d7988ee1 | 73 | const LocOpe_Revol&, |
74 | TopTools_DataMapOfShapeListOfShape&, // myMap | |
75 | TopoDS_Shape&, // myFShape | |
76 | TopoDS_Shape&); // myLShape | |
7fd59977 | 77 | |
78 | ||
7fd59977 | 79 | static void VerifGluedFaces(const TopoDS_Face& theSkface, |
d7988ee1 | 80 | const TopoDS_Shape& thePbase, |
81 | Handle(Geom_Curve)& theBCurve, | |
82 | TColGeom_SequenceOfCurve& theCurves, | |
83 | LocOpe_Revol& theRevol, | |
84 | TopTools_DataMapOfShapeShape& theMap); | |
7fd59977 | 85 | |
86 | ||
87 | static Standard_Boolean ToFuse(const TopoDS_Face& , | |
d7988ee1 | 88 | const TopoDS_Face&); |
7fd59977 | 89 | |
90 | ||
91 | ||
92 | ||
93 | //======================================================================= | |
94 | //function : Init | |
95 | //purpose : | |
96 | //======================================================================= | |
97 | ||
98 | void BRepFeat_MakeRevol::Init(const TopoDS_Shape& Sbase, | |
d7988ee1 | 99 | const TopoDS_Shape& Pbase, |
100 | const TopoDS_Face& Skface, | |
101 | const gp_Ax1& Axis, | |
102 | const Standard_Integer Mode, | |
103 | const Standard_Boolean Modify) | |
7fd59977 | 104 | { |
0797d9d3 | 105 | #ifdef OCCT_DEBUG |
7fd59977 | 106 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 107 | if (trc) std::cout << "BRepFeat_MakeRevol::Init" << std::endl; |
7fd59977 | 108 | #endif |
109 | myAxis = Axis; | |
110 | myPbase = Pbase; | |
111 | mySbase = Sbase; | |
112 | BasisShapeValid(); | |
113 | mySkface = Skface; | |
114 | SketchFaceValid(); | |
115 | myPbase = Pbase; | |
116 | mySlface.Clear(); | |
117 | if(Mode == 0) { | |
118 | myFuse = Standard_False; | |
119 | myJustFeat = Standard_False; | |
120 | } | |
121 | else if(Mode == 1) { | |
122 | myFuse = Standard_True; | |
123 | myJustFeat = Standard_False; | |
124 | } | |
125 | else if(Mode == 2) { | |
126 | myFuse = Standard_True; | |
127 | myJustFeat = Standard_True; | |
128 | } | |
129 | else { | |
130 | } | |
131 | myModify = Modify; | |
132 | myJustGluer = Standard_False; | |
133 | ||
134 | //-------------- ifv | |
d7988ee1 | 135 | // mySkface.Nullify(); |
7fd59977 | 136 | //-------------- ifv |
137 | ||
138 | ||
139 | myShape.Nullify(); | |
140 | myMap.Clear(); | |
141 | myFShape.Nullify(); | |
142 | myLShape.Nullify(); | |
143 | TopExp_Explorer exp; | |
144 | for (exp.Init(mySbase,TopAbs_FACE);exp.More();exp.Next()) { | |
145 | TopTools_ListOfShape thelist; | |
146 | myMap.Bind(exp.Current(), thelist); | |
147 | myMap(exp.Current()).Append(exp.Current()); | |
148 | } | |
0797d9d3 | 149 | #ifdef OCCT_DEBUG |
7fd59977 | 150 | if (trc) { |
04232180 | 151 | if (myJustFeat) std::cout << " Just Feature" << std::endl; |
152 | if (myFuse) std::cout << " Fuse" << std::endl; | |
153 | if (!myFuse) std::cout << " Cut" << std::endl; | |
154 | if (!myModify) std::cout << " Modify = 0" << std::endl; | |
7fd59977 | 155 | } |
156 | #endif | |
157 | } | |
158 | ||
159 | ||
160 | //======================================================================= | |
161 | //function : Add | |
0d969553 | 162 | //purpose : add faces add edges of sliding |
7fd59977 | 163 | //======================================================================= |
164 | ||
165 | void BRepFeat_MakeRevol::Add(const TopoDS_Edge& E, | |
d7988ee1 | 166 | const TopoDS_Face& F) |
7fd59977 | 167 | { |
0797d9d3 | 168 | #ifdef OCCT_DEBUG |
7fd59977 | 169 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 170 | if (trc) std::cout << "BRepFeat_MakeRevol::Add(Edge,face)" << std::endl; |
7fd59977 | 171 | #endif |
172 | TopExp_Explorer exp; | |
173 | for (exp.Init(mySbase,TopAbs_FACE);exp.More();exp.Next()) { | |
174 | if (exp.Current().IsSame(F)) { | |
175 | break; | |
176 | } | |
177 | } | |
178 | if (!exp.More()) { | |
9775fa61 | 179 | throw Standard_ConstructionError(); |
7fd59977 | 180 | } |
181 | ||
182 | for (exp.Init(myPbase,TopAbs_EDGE);exp.More();exp.Next()) { | |
183 | if (exp.Current().IsSame(E)) { | |
184 | break; | |
185 | } | |
186 | } | |
187 | if (!exp.More()) { | |
9775fa61 | 188 | throw Standard_ConstructionError(); |
7fd59977 | 189 | } |
190 | ||
191 | if (!mySlface.IsBound(F)) { | |
192 | TopTools_ListOfShape thelist; | |
193 | mySlface.Bind(F, thelist); | |
194 | } | |
195 | TopTools_ListIteratorOfListOfShape itl(mySlface(F)); | |
196 | for (; itl.More();itl.Next()) { | |
197 | if (itl.Value().IsSame(E)) { | |
198 | break; | |
199 | } | |
200 | } | |
201 | if (!itl.More()) { | |
202 | mySlface(F).Append(E); | |
203 | } | |
204 | } | |
205 | ||
206 | ||
207 | //======================================================================= | |
208 | //function : Perform | |
209 | //purpose : | |
210 | //======================================================================= | |
211 | ||
212 | void BRepFeat_MakeRevol::Perform(const Standard_Real Angle) | |
213 | { | |
0797d9d3 | 214 | #ifdef OCCT_DEBUG |
7fd59977 | 215 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 216 | if (trc) std::cout << "BRepFeat_MakeRevol::Perform(Angle)" << std::endl; |
7fd59977 | 217 | #endif |
218 | mySFrom.Nullify(); | |
219 | ShapeFromValid(); | |
220 | mySUntil.Nullify(); | |
221 | ShapeUntilValid(); | |
222 | myGluedF.Clear(); | |
223 | myPerfSelection = BRepFeat_NoSelection; | |
224 | PerfSelectionValid(); | |
c6541a0c | 225 | Standard_Boolean RevolComp = (2*M_PI-Abs(Angle) <= Precision::Angular()); |
7fd59977 | 226 | LocOpe_Revol theRevol; |
227 | Standard_Real angledec = 0.; | |
228 | TopExp_Explorer exp; | |
229 | if(RevolComp) { | |
230 | /* | |
231 | if (!mySkface.IsNull() || !mySlface.IsEmpty()) { | |
d7988ee1 | 232 | for (exp.Init(mySbase,TopAbs_FACE); exp.More(); exp.Next()) { |
233 | if (exp.Current().IsSame(mySkface)) { | |
234 | angledec = M_PI/5; // pourquoi pas | |
235 | if (myFuse) angledec = -angledec; | |
236 | break; | |
237 | } | |
238 | } | |
7fd59977 | 239 | } |
240 | */ | |
241 | mySkface.Nullify(); | |
242 | } | |
243 | if(angledec == 0.) theRevol.Perform(myPbase, myAxis, Angle); | |
244 | else theRevol.Perform(myPbase, myAxis, Angle, angledec); | |
245 | ||
246 | TopoDS_Shape VraiRevol = theRevol.Shape(); | |
d7988ee1 | 247 | |
7fd59977 | 248 | MajMap(myPbase,theRevol,myMap,myFShape,myLShape); |
249 | ||
250 | myGShape = VraiRevol; | |
251 | GeneratedShapeValid(); | |
252 | TopoDS_Shape Base = theRevol.FirstShape(); | |
253 | exp.Init(Base, TopAbs_FACE); | |
254 | TopoDS_Face theBase = TopoDS::Face(exp.Current()); | |
255 | exp.Next(); | |
256 | if(exp.More()) { | |
257 | NotDone(); | |
258 | myStatusError = BRepFeat_InvFirstShape; | |
259 | return; | |
260 | } | |
261 | ||
262 | TopoDS_Face FFace; | |
d7988ee1 | 263 | |
7fd59977 | 264 | Standard_Boolean found = Standard_False; |
265 | ||
266 | if(!mySkface.IsNull() || !mySlface.IsEmpty()) { | |
267 | if(myLShape.ShapeType() == TopAbs_WIRE) { | |
268 | TopExp_Explorer ex1(VraiRevol, TopAbs_FACE); | |
269 | for(; ex1.More(); ex1.Next()) { | |
d7988ee1 | 270 | TopExp_Explorer ex2(ex1.Current(), TopAbs_WIRE); |
271 | for(; ex2.More(); ex2.Next()) { | |
272 | if(ex2.Current().IsSame(myLShape)) { | |
273 | FFace = TopoDS::Face(ex1.Current()); | |
274 | found = Standard_True; | |
275 | break; | |
276 | } | |
277 | } | |
278 | if(found) break; | |
7fd59977 | 279 | } |
280 | } | |
d7988ee1 | 281 | |
51740958 | 282 | TopExp_Explorer anExp(mySbase, TopAbs_FACE); |
283 | for(; anExp.More(); anExp.Next()) { | |
284 | const TopoDS_Face& ff = TopoDS::Face(anExp.Current()); | |
7fd59977 | 285 | if(ToFuse(ff, FFace)) { |
d7988ee1 | 286 | TopTools_DataMapOfShapeListOfShape sl; |
287 | if(!FFace.IsSame(myPbase) && BRepFeat::IsInside(ff, FFace)) | |
288 | break; | |
7fd59977 | 289 | } |
290 | } | |
291 | } | |
7fd59977 | 292 | GluedFacesValid(); |
293 | if (!mySkface.IsNull()) { | |
294 | VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF); | |
295 | } | |
296 | ||
297 | if(myGluedF.IsEmpty()) { | |
298 | if(myFuse == 1) { | |
7fd59977 | 299 | BRepAlgoAPI_Fuse f(mySbase, myGShape); |
300 | myShape = f.Shape(); | |
301 | UpdateDescendants(f, myShape, Standard_False); | |
7fd59977 | 302 | Done(); |
303 | } | |
304 | else if(myFuse == 0) { | |
7fd59977 | 305 | BRepAlgoAPI_Cut c(mySbase, myGShape); |
306 | myShape = c.Shape(); | |
307 | UpdateDescendants(c, myShape, Standard_False); | |
7fd59977 | 308 | Done(); |
309 | } | |
310 | else { | |
311 | myShape = myGShape; | |
312 | Done(); | |
313 | } | |
314 | } | |
315 | else { | |
316 | theRevol.Curves(myCurves); | |
317 | myBCurve = theRevol.BarycCurve(); | |
318 | GlobalPerform(); | |
319 | } | |
320 | } | |
321 | ||
322 | ||
323 | //======================================================================= | |
324 | //function : Perform | |
0d969553 | 325 | //purpose : feature till shape Until |
7fd59977 | 326 | //======================================================================= |
327 | ||
328 | void BRepFeat_MakeRevol::Perform(const TopoDS_Shape& Until) | |
329 | { | |
0797d9d3 | 330 | #ifdef OCCT_DEBUG |
7fd59977 | 331 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 332 | if (trc) std::cout << "BRepFeat_MakeRevol::Perform(Until)" << std::endl; |
7fd59977 | 333 | #endif |
334 | Standard_Real Angle = 0.; | |
335 | Standard_Boolean TourComplet = Standard_False; | |
336 | ||
337 | if (Until.IsNull()) { | |
9775fa61 | 338 | throw Standard_ConstructionError(); |
7fd59977 | 339 | } |
340 | TopExp_Explorer exp(Until, TopAbs_FACE); | |
341 | if (!exp.More()) { | |
9775fa61 | 342 | throw Standard_ConstructionError(); |
7fd59977 | 343 | } |
344 | if (!mySkface.IsNull() && Until.IsSame(mySkface)) { | |
c6541a0c | 345 | Angle = 2*M_PI; |
7fd59977 | 346 | TourComplet = Standard_True; |
347 | } | |
348 | myGluedF.Clear(); | |
349 | myPerfSelection = BRepFeat_SelectionU; | |
350 | PerfSelectionValid(); | |
351 | mySFrom.Nullify(); | |
352 | ShapeFromValid(); | |
353 | mySUntil = Until; | |
354 | Standard_Boolean Trf = TransformShapeFU(1); | |
355 | ShapeUntilValid(); | |
356 | ||
d7988ee1 | 357 | // Do systematically almost complete revolution |
358 | // BRepSweep_Revol theRevol(myPbase,myAxis,2.*M_PI-10.*Precision::Angular()); | |
7fd59977 | 359 | LocOpe_Revol theRevol; |
360 | if(!TourComplet) { | |
c6541a0c | 361 | Angle = 2.*M_PI- 3*M_PI/180.; |
0797d9d3 | 362 | #ifdef OCCT_DEBUG |
04232180 | 363 | if (trc) std::cout << " No complete Revolution" << std::endl; |
7fd59977 | 364 | #endif |
365 | } | |
366 | theRevol.Perform(myPbase, myAxis, Angle); | |
367 | TopoDS_Shape VraiRevol = theRevol.Shape(); | |
368 | MajMap(myPbase,theRevol,myMap,myFShape,myLShape); | |
369 | ||
370 | ||
371 | if(!Trf) { | |
d7988ee1 | 372 | |
7fd59977 | 373 | myGShape = VraiRevol; |
374 | GeneratedShapeValid(); | |
d7988ee1 | 375 | |
7fd59977 | 376 | TopoDS_Shape Base = theRevol.FirstShape(); |
377 | exp.Init(Base, TopAbs_FACE); | |
378 | TopoDS_Face theBase = TopoDS::Face(exp.Current()); | |
379 | exp.Next(); | |
380 | if(exp.More()) { | |
381 | NotDone(); | |
382 | myStatusError = BRepFeat_InvFirstShape; | |
383 | return; | |
384 | } | |
7fd59977 | 385 | GluedFacesValid(); |
386 | //VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF); | |
387 | ||
388 | theRevol.Curves(myCurves); | |
389 | myBCurve = theRevol.BarycCurve(); | |
390 | GlobalPerform(); | |
391 | } | |
392 | else { | |
393 | TColGeom_SequenceOfCurve scur; | |
394 | theRevol.Curves(myCurves); | |
395 | myBCurve = theRevol.BarycCurve(); | |
396 | scur.Clear(); | |
397 | scur.Append(myBCurve); | |
398 | LocOpe_CSIntersector ASI(mySUntil); | |
399 | ASI.Perform(scur); | |
400 | if (ASI.IsDone() && ASI.NbPoints(1) >=1) { | |
401 | TopAbs_Orientation Or = ASI.Point(1,1).Orientation(); | |
402 | TopoDS_Face FUntil = ASI.Point(1,1).Face(); | |
403 | TopoDS_Shape Comp; | |
404 | BRep_Builder B; | |
405 | B.MakeCompound(TopoDS::Compound(Comp)); | |
406 | TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, Or); | |
407 | if (!S.IsNull()) B.Add(Comp,S); | |
7fd59977 | 408 | BRepAlgoAPI_Cut trP(VraiRevol,Comp); |
7fd59977 | 409 | TopoDS_Shape Cutsh = trP.Shape(); |
410 | TopExp_Explorer ex(Cutsh, TopAbs_SOLID); | |
411 | for(; ex.More(); ex.Next()) { | |
d7988ee1 | 412 | TopExp_Explorer ex1(ex.Current(), TopAbs_FACE); |
413 | for(; ex1.More(); ex1.Next()) { | |
414 | const TopoDS_Face& fac = TopoDS::Face(ex1.Current()); | |
415 | if(fac.IsSame(myPbase)) { | |
416 | VraiRevol = ex.Current(); | |
417 | break; | |
418 | } | |
419 | } | |
7fd59977 | 420 | } |
421 | if(myFuse == 1) { | |
d7988ee1 | 422 | BRepAlgoAPI_Fuse f(mySbase, VraiRevol); |
423 | myShape = f.Shape(); | |
424 | UpdateDescendants(f, myShape, Standard_False); | |
d7988ee1 | 425 | Done(); |
7fd59977 | 426 | } |
427 | else if(myFuse == 0) { | |
d7988ee1 | 428 | BRepAlgoAPI_Cut c(mySbase, VraiRevol); |
429 | myShape = c.Shape(); | |
430 | UpdateDescendants(c, myShape, Standard_False); | |
d7988ee1 | 431 | Done(); |
7fd59977 | 432 | } |
433 | else { | |
d7988ee1 | 434 | myShape = VraiRevol; |
435 | Done(); | |
7fd59977 | 436 | } |
437 | } | |
438 | } | |
d7988ee1 | 439 | } |
7fd59977 | 440 | |
441 | ||
442 | //======================================================================= | |
443 | //function : Perform | |
0d969553 | 444 | //purpose : feature limited by two shapes |
7fd59977 | 445 | //======================================================================= |
446 | ||
447 | void BRepFeat_MakeRevol::Perform(const TopoDS_Shape& From, | |
d7988ee1 | 448 | const TopoDS_Shape& Until) |
7fd59977 | 449 | { |
0797d9d3 | 450 | #ifdef OCCT_DEBUG |
7fd59977 | 451 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 452 | if (trc) std::cout << "BRepFeat_MakeRevol::Perform(From,Until)" << std::endl; |
7fd59977 | 453 | #endif |
454 | if (From.IsNull() || Until.IsNull()) { | |
9775fa61 | 455 | throw Standard_ConstructionError(); |
7fd59977 | 456 | } |
457 | if (!mySkface.IsNull()) { | |
458 | if (From.IsSame(mySkface)) { | |
459 | myJustGluer = Standard_True; | |
460 | Perform(Until); | |
461 | if (myJustGluer) return; | |
462 | } | |
463 | else if (Until.IsSame(mySkface)) { | |
464 | myJustGluer = Standard_True; | |
465 | myAxis.Reverse(); | |
466 | Perform(From); | |
467 | if (myJustGluer) return; | |
468 | } | |
469 | } | |
470 | ||
471 | myGluedF.Clear(); | |
472 | myPerfSelection = BRepFeat_SelectionFU; | |
473 | PerfSelectionValid(); | |
474 | ||
475 | TopExp_Explorer exp(From, TopAbs_FACE); | |
476 | if (!exp.More()) { | |
9775fa61 | 477 | throw Standard_ConstructionError(); |
7fd59977 | 478 | } |
479 | exp.Init(Until, TopAbs_FACE); | |
480 | if (!exp.More()) { | |
9775fa61 | 481 | throw Standard_ConstructionError(); |
7fd59977 | 482 | } |
483 | ||
484 | mySFrom = From; | |
485 | Standard_Boolean Trff = TransformShapeFU(0); | |
486 | ShapeFromValid(); | |
487 | mySUntil = Until; | |
488 | Standard_Boolean Trfu = TransformShapeFU(1); | |
489 | ShapeUntilValid(); | |
d7988ee1 | 490 | |
7fd59977 | 491 | if(Trfu != Trff) { |
492 | NotDone(); | |
493 | myStatusError = BRepFeat_IncTypes; | |
494 | return; | |
495 | } | |
496 | ||
497 | LocOpe_Revol theRevol; | |
c6541a0c | 498 | theRevol.Perform(myPbase, myAxis, 2*M_PI); |
7fd59977 | 499 | TopoDS_Shape VraiRevol = theRevol.Shape(); |
500 | ||
501 | MajMap(myPbase,theRevol,myMap,myFShape,myLShape); | |
502 | ||
503 | if(!Trff) { | |
504 | myGShape = VraiRevol; | |
505 | GeneratedShapeValid(); | |
7fd59977 | 506 | GluedFacesValid(); |
d7988ee1 | 507 | // VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF); |
7fd59977 | 508 | |
509 | theRevol.Curves(myCurves); | |
510 | myBCurve = theRevol.BarycCurve(); | |
511 | GlobalPerform(); | |
512 | } | |
513 | else { | |
514 | theRevol.Curves(myCurves); | |
515 | myBCurve = theRevol.BarycCurve(); | |
516 | TColGeom_SequenceOfCurve scur; | |
517 | scur.Clear(); | |
518 | scur.Append(myBCurve); | |
519 | LocOpe_CSIntersector ASI1(mySUntil); | |
520 | LocOpe_CSIntersector ASI2(mySFrom); | |
521 | ASI1.Perform(scur); | |
522 | ASI2.Perform(scur); | |
523 | TopAbs_Orientation OrU, OrF; | |
524 | TopoDS_Face FFrom, FUntil; | |
525 | Standard_Real PrF, PrU; | |
526 | if (ASI1.IsDone() && ASI1.NbPoints(1) >=1) { | |
527 | OrU = ASI1.Point(1,1).Orientation(); | |
528 | FUntil = ASI1.Point(1,1).Face(); | |
529 | PrU = ASI1.Point(1,1).Parameter(); | |
530 | } | |
531 | else { | |
532 | NotDone(); | |
533 | myStatusError = BRepFeat_NoIntersectU; | |
534 | return; | |
535 | } | |
536 | if (ASI2.IsDone() && ASI2.NbPoints(1) >=1) { | |
537 | Standard_Real pr1 = ASI2.Point(1,1).Parameter(); | |
c6541a0c | 538 | pr1 = ElCLib::InPeriod(pr1,PrU-2*M_PI,PrU); |
7fd59977 | 539 | Standard_Real pr2 = ASI2.Point(1,ASI2.NbPoints(1)).Parameter(); |
c6541a0c | 540 | pr2 = ElCLib::InPeriod(pr2,PrU-2*M_PI,PrU); |
d7988ee1 | 541 | //OrF = OrU; |
542 | OrF = TopAbs::Reverse(OrU); | |
7fd59977 | 543 | FFrom = ASI2.Point(1,1).Face(); |
544 | PrF = Max(pr1, pr2); | |
545 | } | |
546 | else { | |
547 | NotDone(); | |
548 | myStatusError = BRepFeat_NoIntersectF; | |
549 | return; | |
550 | } | |
551 | if(!(PrU > PrF)) { | |
552 | NotDone(); | |
553 | myStatusError = BRepFeat_IncParameter; | |
554 | return; | |
555 | } | |
556 | TopoDS_Shape Comp; | |
557 | BRep_Builder B; | |
558 | B.MakeCompound(TopoDS::Compound(Comp)); | |
559 | TopoDS_Solid SF = BRepFeat::Tool(mySFrom, FFrom, OrF); | |
560 | if (!SF.IsNull()) B.Add(Comp,SF); | |
561 | TopoDS_Solid SU = BRepFeat::Tool(mySUntil, FUntil, OrU); | |
562 | if (!SU.IsNull()) B.Add(Comp,SU); | |
7fd59977 | 563 | BRepAlgoAPI_Cut trP(VraiRevol,Comp); |
7fd59977 | 564 | TopoDS_Shape Cutsh = trP.Shape(); |
565 | TopExp_Explorer ex(Cutsh, TopAbs_SOLID); | |
7fd59977 | 566 | VraiRevol = ex.Current(); |
567 | for(; ex.More(); ex.Next()) { | |
568 | Standard_Real PrCur = BRepFeat:: | |
d7988ee1 | 569 | ParametricBarycenter(ex.Current(), myBCurve); |
7fd59977 | 570 | if(PrF <= PrCur && PrU >= PrCur) { |
d7988ee1 | 571 | VraiRevol = ex.Current(); |
572 | break; | |
7fd59977 | 573 | } |
574 | } | |
d7988ee1 | 575 | if(myFuse == 1 && !myJustFeat) { |
7fd59977 | 576 | BRepAlgoAPI_Fuse f(mySbase, VraiRevol); |
577 | myShape = f.Shape(); | |
578 | UpdateDescendants(f, myShape, Standard_False); | |
7fd59977 | 579 | Done(); |
580 | } | |
d7988ee1 | 581 | else if(myFuse == 0 && !myJustFeat) { |
7fd59977 | 582 | BRepAlgoAPI_Cut c(mySbase, VraiRevol); |
583 | myShape = c.Shape(); | |
584 | UpdateDescendants(c, myShape, Standard_False); | |
7fd59977 | 585 | Done(); |
586 | } | |
587 | else { | |
588 | myShape = VraiRevol; | |
589 | Done(); | |
590 | } | |
591 | } | |
592 | } | |
593 | ||
594 | ||
595 | //======================================================================= | |
596 | //function : PerformThruAll | |
0d969553 | 597 | //purpose : feature throughout the initial shape |
7fd59977 | 598 | //======================================================================= |
599 | ||
600 | void BRepFeat_MakeRevol::PerformThruAll() | |
601 | { | |
0797d9d3 | 602 | #ifdef OCCT_DEBUG |
7fd59977 | 603 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 604 | if (trc) std::cout << "BRepFeat_MakeRevol::PerformThruAll()" << std::endl; |
7fd59977 | 605 | #endif |
c6541a0c | 606 | Perform(2.*M_PI); |
7fd59977 | 607 | } |
608 | ||
609 | //======================================================================= | |
610 | //function : PerformUntilAngle | |
0d969553 | 611 | //purpose : feature till shape Until defined with the angle |
7fd59977 | 612 | //======================================================================= |
613 | ||
614 | void BRepFeat_MakeRevol::PerformUntilAngle(const TopoDS_Shape& Until, | |
d7988ee1 | 615 | const Standard_Real Angle) |
7fd59977 | 616 | { |
0797d9d3 | 617 | #ifdef OCCT_DEBUG |
7fd59977 | 618 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 619 | if (trc) std::cout << "BRepFeat_MakeRevol::PerformUntilAngle(Until,Angle)" << std::endl; |
7fd59977 | 620 | #endif |
621 | if (Until.IsNull()) { | |
622 | Perform(Angle); | |
623 | } | |
624 | if(Angle == 0) { | |
625 | Perform(Until); | |
626 | } | |
627 | TopExp_Explorer exp(Until, TopAbs_FACE); | |
628 | if (!exp.More()) { | |
9775fa61 | 629 | throw Standard_ConstructionError(); |
7fd59977 | 630 | } |
631 | if (!mySkface.IsNull() && Until.IsSame(mySkface)) { | |
632 | Perform(Angle); | |
633 | return; | |
634 | } | |
635 | myGluedF.Clear(); | |
636 | myPerfSelection = BRepFeat_NoSelection; | |
637 | PerfSelectionValid(); | |
638 | mySFrom.Nullify(); | |
639 | ShapeFromValid(); | |
640 | mySUntil = Until; | |
641 | Standard_Boolean Trf = TransformShapeFU(1); | |
642 | ShapeUntilValid(); | |
643 | ||
d7988ee1 | 644 | // Produce systematicallt an almost complete revolution |
645 | // BRepSweep_Revol theRevol(myPbase,myAxis,2.*M_PI-10.*Precision::Angular()); | |
7fd59977 | 646 | LocOpe_Revol theRevol; |
647 | theRevol.Perform(myPbase, myAxis, Angle); | |
648 | TopoDS_Shape VraiRevol = theRevol.Shape(); | |
649 | ||
650 | MajMap(myPbase,theRevol,myMap,myFShape,myLShape); | |
651 | ||
652 | if(Trf) { | |
653 | myGShape = VraiRevol; | |
654 | GeneratedShapeValid(); | |
d7988ee1 | 655 | |
7fd59977 | 656 | TopoDS_Shape Base = theRevol.FirstShape(); |
657 | exp.Init(Base, TopAbs_FACE); | |
658 | TopoDS_Face theBase = TopoDS::Face(exp.Current()); | |
659 | exp.Next(); | |
660 | if(exp.More()) { | |
661 | NotDone(); | |
662 | myStatusError = BRepFeat_InvFirstShape; | |
663 | return; | |
664 | } | |
7fd59977 | 665 | GluedFacesValid(); |
666 | //VerifGluedFaces(mySkface, theBase, myBCurve, myCurves, theRevol, myGluedF); | |
667 | ||
668 | ||
669 | theRevol.Curves(myCurves); | |
670 | myBCurve = theRevol.BarycCurve(); | |
671 | GlobalPerform(); | |
672 | } | |
673 | else { | |
674 | TColGeom_SequenceOfCurve scur; | |
675 | theRevol.Curves(myCurves); | |
676 | myBCurve = theRevol.BarycCurve(); | |
677 | scur.Clear(); | |
678 | scur.Append(myBCurve); | |
679 | LocOpe_CSIntersector ASI(mySUntil); | |
680 | ASI.Perform(scur); | |
681 | if (ASI.IsDone() && ASI.NbPoints(1) >=1) { | |
682 | TopAbs_Orientation Or = ASI.Point(1,1).Orientation(); | |
683 | TopoDS_Face FUntil = ASI.Point(1,1).Face(); | |
684 | TopoDS_Shape Comp; | |
685 | BRep_Builder B; | |
686 | B.MakeCompound(TopoDS::Compound(Comp)); | |
687 | TopoDS_Solid S = BRepFeat::Tool(mySUntil, FUntil, Or); | |
688 | if (!S.IsNull()) B.Add(Comp,S); | |
7fd59977 | 689 | BRepAlgoAPI_Cut trP(VraiRevol,Comp); |
7fd59977 | 690 | TopoDS_Shape Cutsh = trP.Shape(); |
691 | TopExp_Explorer ex(Cutsh, TopAbs_SOLID); | |
692 | for(; ex.More(); ex.Next()) { | |
d7988ee1 | 693 | TopExp_Explorer ex1(ex.Current(), TopAbs_FACE); |
694 | for(; ex1.More(); ex1.Next()) { | |
695 | const TopoDS_Face& fac = TopoDS::Face(ex1.Current()); | |
696 | if(fac.IsSame(myPbase)) { | |
697 | VraiRevol = ex.Current(); | |
698 | break; | |
699 | } | |
700 | } | |
7fd59977 | 701 | } |
702 | if(myFuse == 1) { | |
d7988ee1 | 703 | BRepAlgoAPI_Fuse f(mySbase, VraiRevol); |
704 | myShape = f.Shape(); | |
705 | UpdateDescendants(f, myShape, Standard_False); | |
d7988ee1 | 706 | Done(); |
7fd59977 | 707 | } |
708 | else if(myFuse == 0) { | |
d7988ee1 | 709 | BRepAlgoAPI_Cut c(mySbase, VraiRevol); |
710 | myShape = c.Shape(); | |
711 | UpdateDescendants(c, myShape, Standard_False); | |
d7988ee1 | 712 | Done(); |
7fd59977 | 713 | } |
714 | else { | |
d7988ee1 | 715 | myShape = VraiRevol; |
716 | Done(); | |
7fd59977 | 717 | } |
718 | } | |
719 | } | |
720 | } | |
721 | //======================================================================= | |
722 | //function : Curves | |
0d969553 | 723 | //purpose : circles parallel to the generating edge of revolution |
7fd59977 | 724 | //======================================================================= |
725 | ||
726 | void BRepFeat_MakeRevol::Curves(TColGeom_SequenceOfCurve& scur) | |
727 | { | |
728 | scur = myCurves; | |
729 | } | |
730 | ||
731 | //======================================================================= | |
732 | //function : BarycCurve | |
0d969553 | 733 | //purpose : pass through the center of mass of the primitive |
7fd59977 | 734 | //======================================================================= |
735 | ||
736 | Handle(Geom_Curve) BRepFeat_MakeRevol::BarycCurve() | |
737 | { | |
738 | return myBCurve; | |
739 | } | |
740 | ||
7fd59977 | 741 | //======================================================================= |
742 | //function : VerifGluedFaces | |
0d969553 Y |
743 | //purpose : Check intersection Tool/theSkface = thePbase |
744 | // if yes -> OK otherwise -> case without gluing | |
7fd59977 | 745 | //======================================================================= |
746 | ||
747 | static void VerifGluedFaces(const TopoDS_Face& theSkface, | |
d7988ee1 | 748 | const TopoDS_Shape& thePbase, |
749 | Handle(Geom_Curve)& theBCurve, | |
750 | TColGeom_SequenceOfCurve& theCurves, | |
751 | LocOpe_Revol& theRevol, | |
752 | TopTools_DataMapOfShapeShape& theMap) | |
7fd59977 | 753 | { |
754 | Standard_Boolean GluedFaces = Standard_True; | |
755 | TopoDS_Shape VraiRevol = theRevol.Shape(); | |
d7988ee1 | 756 | |
7fd59977 | 757 | TColGeom_SequenceOfCurve scur; |
758 | theRevol.Curves(theCurves); | |
759 | theBCurve = theRevol.BarycCurve(); | |
760 | scur.Clear(); | |
761 | scur.Append(theBCurve); | |
762 | LocOpe_CSIntersector ASI(theSkface); | |
763 | ASI.Perform(scur); | |
764 | if (ASI.IsDone() && ASI.NbPoints(1) >=1) { | |
765 | TopAbs_Orientation Or = ASI.Point(1,1).Orientation(); | |
766 | TopoDS_Face FSk = ASI.Point(1,1).Face(); | |
767 | TopoDS_Shape Comp; | |
768 | BRep_Builder B; | |
769 | B.MakeCompound(TopoDS::Compound(Comp)); | |
770 | TopoDS_Solid S = BRepFeat::Tool(theSkface, FSk, Or); | |
771 | if (!S.IsNull()) B.Add(Comp,S); | |
7fd59977 | 772 | BRepAlgoAPI_Cut trP(VraiRevol,Comp); |
7fd59977 | 773 | TopoDS_Shape Cutsh = trP.Shape(); |
774 | TopExp_Explorer ex(Cutsh, TopAbs_SOLID); | |
775 | for(; ex.More(); ex.Next()) { | |
776 | TopExp_Explorer ex1(ex.Current(), TopAbs_FACE); | |
777 | for(; ex1.More(); ex1.Next()) { | |
d7988ee1 | 778 | const TopoDS_Face& fac1 = TopoDS::Face(ex1.Current()); |
779 | TopExp_Explorer ex2(thePbase, TopAbs_FACE); | |
780 | for(; ex2.More(); ex2.Next()) { | |
781 | const TopoDS_Face& fac2 = TopoDS::Face(ex2.Current()); | |
782 | if(fac1.IsSame(fac2)) break; | |
783 | } | |
784 | if (ex2.More()) break; | |
7fd59977 | 785 | } |
786 | if (ex1.More()) continue; | |
787 | GluedFaces = Standard_False; | |
788 | break; | |
789 | } | |
790 | if (!GluedFaces) { | |
0797d9d3 | 791 | #ifdef OCCT_DEBUG |
7fd59977 | 792 | Standard_Boolean trc = BRepFeat_GettraceFEAT(); |
04232180 | 793 | if (trc) std::cout << " Intersection Revol/skface : no gluing" << std::endl; |
7fd59977 | 794 | #endif |
795 | theMap.Clear(); | |
796 | } | |
797 | } | |
798 | } | |
799 | ||
800 | //======================================================================= | |
801 | //function : MajMap | |
0d969553 | 802 | //purpose : management of descendants |
7fd59977 | 803 | //======================================================================= |
804 | ||
805 | static void MajMap(const TopoDS_Shape& theB, | |
d7988ee1 | 806 | const LocOpe_Revol& theP, |
807 | TopTools_DataMapOfShapeListOfShape& theMap, // myMap | |
808 | TopoDS_Shape& theFShape, // myFShape | |
809 | TopoDS_Shape& theLShape) // myLShape | |
7fd59977 | 810 | { |
811 | TopExp_Explorer exp(theP.FirstShape(),TopAbs_WIRE); | |
812 | if (exp.More()) { | |
813 | theFShape = exp.Current(); | |
814 | TopTools_ListOfShape thelist; | |
815 | theMap.Bind(theFShape, thelist); | |
816 | for (exp.Init(theP.FirstShape(),TopAbs_FACE);exp.More();exp.Next()) { | |
817 | theMap(theFShape).Append(exp.Current()); | |
818 | } | |
819 | } | |
d7988ee1 | 820 | |
7fd59977 | 821 | exp.Init(theP.LastShape(),TopAbs_WIRE); |
822 | if (exp.More()) { | |
823 | theLShape = exp.Current(); | |
824 | TopTools_ListOfShape thelist1; | |
825 | theMap.Bind(theLShape, thelist1); | |
826 | for (exp.Init(theP.LastShape(),TopAbs_FACE);exp.More();exp.Next()) { | |
827 | theMap(theLShape).Append(exp.Current()); | |
828 | } | |
829 | } | |
830 | ||
831 | for (exp.Init(theB,TopAbs_EDGE); exp.More(); exp.Next()) { | |
832 | if (!theMap.IsBound(exp.Current())) { | |
833 | TopTools_ListOfShape thelist2; | |
834 | theMap.Bind(exp.Current(), thelist2); | |
835 | theMap(exp.Current()) = theP.Shapes(exp.Current()); | |
836 | } | |
837 | } | |
838 | } | |
839 | ||
840 | ||
841 | ||
842 | //======================================================================= | |
843 | //function : ToFuse | |
0d969553 | 844 | //purpose : two faces samedomaine or not |
7fd59977 | 845 | //======================================================================= |
846 | ||
847 | Standard_Boolean ToFuse(const TopoDS_Face& F1, | |
d7988ee1 | 848 | const TopoDS_Face& F2) |
7fd59977 | 849 | { |
850 | if (F1.IsNull() || F2.IsNull()) { | |
851 | return Standard_False; | |
852 | } | |
853 | ||
854 | Handle(Geom_Surface) S1,S2; | |
855 | TopLoc_Location loc1, loc2; | |
856 | Handle(Standard_Type) typS1,typS2; | |
857 | const Standard_Real tollin = Precision::Confusion(); | |
858 | const Standard_Real tolang = Precision::Angular(); | |
859 | ||
860 | S1 = BRep_Tool::Surface(F1,loc1); | |
861 | S2 = BRep_Tool::Surface(F2,loc2); | |
862 | ||
863 | typS1 = S1->DynamicType(); | |
864 | typS2 = S2->DynamicType(); | |
865 | ||
866 | if (typS1 == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) { | |
c5f3a425 | 867 | S1 = Handle(Geom_RectangularTrimmedSurface)::DownCast (S1)->BasisSurface(); |
7fd59977 | 868 | typS1 = S1->DynamicType(); |
869 | } | |
870 | ||
871 | if (typS2 == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) { | |
c5f3a425 | 872 | S2 = Handle(Geom_RectangularTrimmedSurface)::DownCast (S2)->BasisSurface(); |
7fd59977 | 873 | typS2 = S2->DynamicType(); |
874 | } | |
875 | ||
876 | if (typS1 != typS2) { | |
877 | return Standard_False; | |
878 | } | |
879 | ||
880 | ||
881 | Standard_Boolean ValRet = Standard_False; | |
882 | if (typS1 == STANDARD_TYPE(Geom_Plane)) { | |
0d969553 | 883 | S1 = BRep_Tool::Surface(F1); // to apply the location. |
7fd59977 | 884 | S2 = BRep_Tool::Surface(F2); |
c5f3a425 | 885 | gp_Pln pl1( Handle(Geom_Plane)::DownCast (S1)->Pln()); |
886 | gp_Pln pl2( Handle(Geom_Plane)::DownCast (S2)->Pln()); | |
7fd59977 | 887 | |
888 | if (pl1.Position().IsCoplanar(pl2.Position(),tollin,tolang)) { | |
889 | ValRet = Standard_True; | |
890 | } | |
891 | } | |
892 | ||
893 | return ValRet; | |
894 | } | |
895 | ||
896 | ||
897 | ||
898 | ||
899 | ||
900 | ||
901 | ||
902 | ||
903 |