1 // Copyright (c) 1997-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
3 //
4 // This file is part of Open CASCADE Technology software library.
5 //
6 // This library is free software; you can redistribute it and/or modify it under
7 // the terms of the GNU Lesser General Public License version 2.1 as published
8 // by the Free Software Foundation, with special exception defined in the file
9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
10 // distribution for complete text of the license and disclaimer of any warranty.
11 //
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
14
15 #include <FilletSurf_InternalBuilder.hxx>
16
17 #include <Adaptor3d_TopolTool.hxx>
18 #include <BRep_Tool.hxx>
19 #include <BRepAdaptor_Curve.hxx>
20 #include <BRepAdaptor_Curve2d.hxx>
21 #include <BRepAdaptor_Surface.hxx>
22 #include <BRepBlend_ConstRad.hxx>
23 #include <BRepBlend_ConstRadInv.hxx>
24 #include <BRepBlend_Line.hxx>
25 #include <BRepExtrema_DistShapeShape.hxx>
26 #include <ChFi3d_Builder_0.hxx>
27 #include <ChFi3d_FilBuilder.hxx>
28 #include <ChFiDS_CircSection.hxx>
29 #include <ChFiDS_CommonPoint.hxx>
30 #include <ChFiDS_FaceInterference.hxx>
31 #include <ChFiDS_FilSpine.hxx>
32 #include <ChFiDS_HData.hxx>
33 #include <ChFiDS_ElSpine.hxx>
34 #include <ChFiDS_ListIteratorOfListOfStripe.hxx>
35 #include <ChFiDS_ListOfStripe.hxx>
36 #include <ChFiDS_SecHArray1.hxx>
37 #include <ChFiDS_SequenceOfSurfData.hxx>
38 #include <ChFiDS_Spine.hxx>
39 #include <ChFiDS_Stripe.hxx>
40 #include <ChFiDS_SurfData.hxx>
41 #include <ElSLib.hxx>
42 #include <FilletSurf_ErrorTypeStatus.hxx>
43 #include <FilletSurf_StatusDone.hxx>
44 #include <FilletSurf_StatusType.hxx>
45 #include <Geom2d_Curve.hxx>
46 #include <Geom_Circle.hxx>
47 #include <Geom_Curve.hxx>
48 #include <Geom_Plane.hxx>
49 #include <Geom_Surface.hxx>
50 #include <Geom_TrimmedCurve.hxx>
51 #include <GeomAbs_Shape.hxx>
52 #include <GeomAdaptor_Curve.hxx>
53 #include <GeomAdaptor_Surface.hxx>
54 #include <gp_Circ.hxx>
55 #include <gp_Pln.hxx>
56 #include <IntCurveSurface_HInter.hxx>
57 #include <IntCurveSurface_IntersectionPoint.hxx>
58 #include <Precision.hxx>
59 #include <Standard_ConstructionError.hxx>
60 #include <TopoDS.hxx>
61 #include <TopoDS_Face.hxx>
62 #include <TopoDS_Shape.hxx>
63 #include <TopOpeBRepDS_HDataStructure.hxx>
64 #include <TopTools_ListIteratorOfListOfShape.hxx>
65
isinlist(const TopoDS_Shape & E,const TopTools_ListOfShape & L)66 static Standard_Boolean isinlist(const TopoDS_Shape& E,
67 const TopTools_ListOfShape& L){
68 TopTools_ListIteratorOfListOfShape It;
69 for (It.Initialize(L); It.More(); It.Next()){
70 if(E.IsSame(It.Value())) return 1;
71 }
72 return 0;
73 }
74
IntPlanEdge(Handle (BRepAdaptor_Curve)& Ed,const gp_Pln & P,Standard_Real & w,const Standard_Real tol3d)75 static Standard_Boolean IntPlanEdge(Handle(BRepAdaptor_Curve)& Ed,
76 const gp_Pln& P,
77 Standard_Real& w,
78 const Standard_Real tol3d)
79 {
80 Standard_Boolean done = 0;
81 Standard_Real f = Ed->FirstParameter();
82 Standard_Real l = Ed->LastParameter();
83 gp_Pnt Or = P.Location();
84 Handle(Geom_Plane) Pln = new Geom_Plane(P);
85 Handle(GeomAdaptor_Surface)
86 Plan = new GeomAdaptor_Surface(GeomAdaptor_Surface(Pln));
87
88 IntCurveSurface_HInter Intersection;
89 Standard_Integer nbp ,iip;
90 IntCurveSurface_IntersectionPoint IP;
91 Standard_Real dist = RealLast();
92
93 Intersection.Perform(Ed,Plan);
94
95 if(Intersection.IsDone()){
96 nbp = Intersection.NbPoints();
97 for (iip = 1; iip <= nbp; iip++){
98 IP = Intersection.Point(iip);
99 gp_Pnt pint = IP.Pnt();
100 Standard_Real d = pint.Distance(Or);
101 if(d<dist){
102 done = 1;
103 w = IP.W();
104 dist = d;
105 }
106 }
107 }
108 gp_Pnt pdeb = Ed->Value(f);
109 gp_Pnt pfin = Ed->Value(l);
110 Standard_Real u,v;
111 //check if the extremities are not solution
112 ElSLib::Parameters(P,pdeb,u,v);
113 gp_Pnt projdeb = ElSLib::Value(u,v,P);
114 Standard_Real dprojdeb = pdeb.Distance(projdeb);
115 if(dprojdeb<tol3d){
116 Standard_Real d = pdeb.Distance(Or);
117 if(d<dist){
118 done = 1;
119 w = f;
120 dist = d;
121 }
122 }
123 ElSLib::Parameters(P,pfin,u,v);
124 gp_Pnt projfin = ElSLib::Value(u,v,P);
125 Standard_Real dprojfin = pfin.Distance(projfin);
126 if(dprojfin<tol3d){
127 Standard_Real d = pfin.Distance(Or);
128 if(d<dist){
129 done = 1;
130 w = l;
131 dist = d;
132 }
133 }
134 return done;
135 }
136
ComputeEdgeParameter(const Handle (ChFiDS_Spine)& Spine,const Standard_Integer ind,const Standard_Real pelsp,Standard_Real & ped,const Standard_Real tol3d)137 static Standard_Boolean ComputeEdgeParameter(const Handle(ChFiDS_Spine)& Spine,
138 const Standard_Integer ind,
139 const Standard_Real pelsp,
140 Standard_Real& ped,
141 const Standard_Real tol3d)
142 {
143 Handle(ChFiDS_ElSpine) Guide = Spine->ElSpine(ind);
144 gp_Pnt P; gp_Vec V;
145 Guide->D1(pelsp,P,V);
146 gp_Pln pln(P,V);
147 Handle(BRepAdaptor_Curve) ed = new BRepAdaptor_Curve (Spine->CurrentElementarySpine(ind));
148 return IntPlanEdge(ed,pln,ped,tol3d);
149 }
150
151 //=======================================================================
152 //function : FilletSurf_InternalBuilder
153 //purpose :
154 //=======================================================================
155
FilletSurf_InternalBuilder(const TopoDS_Shape & S,const ChFi3d_FilletShape FShape,const Standard_Real Ta,const Standard_Real Tapp3d,const Standard_Real Tapp2d)156 FilletSurf_InternalBuilder::FilletSurf_InternalBuilder
157 (const TopoDS_Shape& S,
158 const ChFi3d_FilletShape FShape,
159 const Standard_Real Ta,
160 const Standard_Real Tapp3d,
161 const Standard_Real Tapp2d):ChFi3d_FilBuilder(S,FShape,Ta)
162 {
163 SetParams(Ta,Tapp3d,Tapp2d,
164 Tapp3d,Tapp2d,1.e-3);
165 SetContinuity(GeomAbs_C2,Ta);
166 }
167
168 //=======================================================================
169 //function : Add
170 //purpose : creation of spine on a set of edges
171 //
172 // 0 : no problem
173 // 1 : empty list
174 // 2 : non g1 edges
175 // 3 : non G1 adjacent faces
176 // 4 : edge is not on the shape
177 // 5 : edge is not alive
178 //=======================================================================
179
Add(const TopTools_ListOfShape & E,const Standard_Real R)180 Standard_Integer FilletSurf_InternalBuilder::Add(const TopTools_ListOfShape& E,
181 const Standard_Real R)
182 {
183 if(E.IsEmpty()) return 1;
184 TopTools_ListIteratorOfListOfShape It;
185 for(It.Initialize(E); It.More(); It.Next()){
186 TopoDS_Edge cured = TopoDS::Edge(It.Value());
187 if(cured.IsNull()) return 4;
188 if(!myEFMap.Contains(cured)) return 4;
189 //check if the edge is a fracture edge
190 TopoDS_Face ff1,ff2;
191 for(It.Initialize(myEFMap(cured));It.More();It.Next()){
192 if (ff1.IsNull()) {
193 ff1 = TopoDS::Face(It.Value());
194 }
195 else {
196 ff2 = TopoDS::Face(It.Value());
197 if(!ff2.IsSame(ff1)){
198 break;
199 }
200 }
201 }
202 if(ff1.IsNull() || ff2.IsNull()) return 5;
203 if(ff1.IsSame(ff2)) return 5;
204 if(BRep_Tool::Continuity(cured,ff1,ff2) != GeomAbs_C0) return 5;
205 }
206 TopoDS_Edge ed = TopoDS::Edge(E.First());
207 ed.Orientation(TopAbs_FORWARD);
208 ChFi3d_FilBuilder::Add(R,ed);
209 Handle(ChFiDS_Stripe) st = myListStripe.First();
210 Handle(ChFiDS_Spine)& sp = st->ChangeSpine();
211 Standard_Boolean periodic = sp->IsPeriodic();
212
213 //It is checked if edges of list E are in the contour,
214 //the edges that arenot in the list are removed from the contour,
215 //it is checked that the remainder is monoblock.
216
217 for(It.Initialize(E); It.More(); It.Next()){
218 TopoDS_Edge cured = TopoDS::Edge(It.Value());
219 if(!Contains(cured)) return 2;
220 }
221
222 Handle(ChFiDS_FilSpine) newsp = new ChFiDS_FilSpine();
223 Standard_Boolean debut = 0;
224 Standard_Integer premierquinyestpas = 0;
225 Standard_Integer yatrou = 0;
226 for(Standard_Integer i = 1; i <= sp->NbEdges(); i++){
227 TopoDS_Edge cured = sp->Edges(i);
228 if(isinlist(cured,E)){
229 debut = 1;
230 if(premierquinyestpas) {
231 yatrou = 1;
232 break;
233 }
234 newsp->SetEdges(cured);
235 }
236 else if(debut && !premierquinyestpas) premierquinyestpas = i;
237 }
238 if(!periodic && yatrou) return 2;
239 if(periodic && yatrou){
240 Standard_Boolean vraitrou = 0, aLocalDebut = 0;
241 for(Standard_Integer i = sp->NbEdges(); i > yatrou; i--){
242 TopoDS_Edge cured = sp->Edges(i);
243 if(isinlist(cured,E)){
244 if(vraitrou) return 2;
245 newsp->PutInFirst(cured);
246 }
247 else if(aLocalDebut) vraitrou = 1;
248 aLocalDebut = 1;
249 }
250 }
251
252 if(newsp->NbEdges() != sp->NbEdges()) {
253 newsp->Load();
254 newsp->SetRadius(R);
255 sp = newsp;
256 }
257
258 //ElSpine is immediately constructed
259 Handle(ChFiDS_ElSpine) hels = new ChFiDS_ElSpine();
260 gp_Vec TFirst,TLast;
261 gp_Pnt PFirst,PLast;
262 sp->D1(sp->FirstParameter(),PFirst,TFirst);
263 sp->D1(sp->LastParameter(),PLast,TLast);
264 hels->FirstParameter(sp->FirstParameter());
265 hels->SetFirstPointAndTgt(PFirst,TFirst);
266 hels->LastParameter(sp->LastParameter());
267 hels->SetLastPointAndTgt(PLast,TLast);
268 ChFi3d_PerformElSpine(hels,sp,myConti,tolesp);
269 sp->AppendElSpine(hels);
270 sp->SplitDone(Standard_True);
271 return 0;
272 }
273
274 //=======================================================================
275 //function : Perform
276 //purpose :
277 //=======================================================================
278
Perform()279 void FilletSurf_InternalBuilder::Perform()
280 {
281 //PerformSetOfSurfOnElSpine is enough.
282
283 Handle(ChFiDS_Stripe) Stripe = myListStripe.First();
284 Handle(ChFiDS_HData)& HData = Stripe->ChangeSetOfSurfData();
285 HData = new ChFiDS_HData();
286 Handle(ChFiDS_Spine)& Spine = Stripe->ChangeSpine();
287 TopAbs_Orientation RefOr1,RefOr2;
288 Standard_Integer RefChoix;
289 StripeOrientations(Spine,RefOr1,RefOr2,RefChoix);
290 Stripe->OrientationOnFace1(RefOr1);
291 Stripe->OrientationOnFace2(RefOr2);
292 Stripe->Choix(RefChoix);
293 PerformSetOfKGen(Stripe,0);
294 }
295
296 //=======================================================================
297 //function : PerformSurf
298 //purpose :
299 //=======================================================================
300
301 Standard_Boolean
PerformSurf(ChFiDS_SequenceOfSurfData & SeqData,const Handle (ChFiDS_ElSpine)& Guide,const Handle (ChFiDS_Spine)& Spine,const Standard_Integer Choix,const Handle (BRepAdaptor_Surface)& S1,const Handle (Adaptor3d_TopolTool)& I1,const Handle (BRepAdaptor_Surface)& S2,const Handle (Adaptor3d_TopolTool)& I2,const Standard_Real MaxStep,const Standard_Real Fleche,const Standard_Real TolGuide,Standard_Real & First,Standard_Real & Last,const Standard_Boolean Inside,const Standard_Boolean Appro,const Standard_Boolean Forward,const Standard_Boolean RecOnS1,const Standard_Boolean RecOnS2,const math_Vector & Soldep,Standard_Integer & Intf,Standard_Integer & Intl)302 FilletSurf_InternalBuilder::PerformSurf(ChFiDS_SequenceOfSurfData& SeqData,
303 const Handle(ChFiDS_ElSpine)& Guide,
304 const Handle(ChFiDS_Spine)& Spine,
305 const Standard_Integer Choix,
306 const Handle(BRepAdaptor_Surface)& S1,
307 const Handle(Adaptor3d_TopolTool)& I1,
308 const Handle(BRepAdaptor_Surface)& S2,
309 const Handle(Adaptor3d_TopolTool)& I2,
310 const Standard_Real MaxStep,
311 const Standard_Real Fleche,
312 const Standard_Real TolGuide,
313 Standard_Real& First,
314 Standard_Real& Last,
315 const Standard_Boolean Inside,
316 const Standard_Boolean Appro,
317 const Standard_Boolean Forward,
318 const Standard_Boolean RecOnS1,
319 const Standard_Boolean RecOnS2,
320 const math_Vector& Soldep,
321 Standard_Integer& Intf,
322 Standard_Integer& Intl)
323 {
324 Handle(ChFiDS_SurfData) Data = SeqData(1);
325 Handle(ChFiDS_FilSpine) fsp = Handle(ChFiDS_FilSpine)::DownCast(Spine);
326 if(fsp.IsNull()) throw Standard_ConstructionError("PerformSurf : this is not the spine of a fillet");
327 Handle(BRepBlend_Line) lin;
328 TopAbs_Orientation Or = S1->Face().Orientation();
329 if(!fsp->IsConstant()) throw Standard_ConstructionError("PerformSurf : no variable radiuses");
330 // Standard_Boolean maybesingular; //pour scinder les Surfdata singulieres
331
332 Handle(ChFiDS_ElSpine) EmptyGuide;
333
334 BRepBlend_ConstRad Func(S1,S2,Guide);
335 BRepBlend_ConstRadInv FInv(S1,S2,Guide);
336 Func.Set(fsp->Radius(),Choix);
337 FInv.Set(fsp->Radius(),Choix);
338 switch (GetFilletShape()) {
339 case ChFi3d_Rational:
340 Func.Set(BlendFunc_Rational);
341 break;
342 case ChFi3d_QuasiAngular:
343 Func.Set(BlendFunc_QuasiAngular);
344 break;
345 case ChFi3d_Polynomial:
346 Func.Set(BlendFunc_Polynomial);
347 }
348 Standard_Real PFirst = First;
349 done = SimulData(Data,Guide,EmptyGuide,lin,S1,I1,
350 S2,I2,Func,FInv,PFirst,MaxStep,Fleche,
351 TolGuide,First,Last,Inside,Appro,Forward,Soldep,
352 20,RecOnS1,RecOnS2);
353 if(!done) return Standard_False;
354 if(lin->StartPointOnFirst().NbPointOnRst() !=0){
355 ChFi3d_FilCommonPoint(lin->StartPointOnFirst(),lin->TransitionOnS1(),
356 Standard_True, Data->ChangeVertexFirstOnS1(), tolesp);
357 }
358 if(lin->EndPointOnFirst().NbPointOnRst() !=0){
359 ChFi3d_FilCommonPoint(lin->EndPointOnFirst(),lin->TransitionOnS1(),
360 Standard_False,Data->ChangeVertexLastOnS1(), tolesp);
361 }
362 if(lin->StartPointOnSecond().NbPointOnRst() !=0){
363 ChFi3d_FilCommonPoint(lin->StartPointOnSecond(),lin->TransitionOnS2(),
364 Standard_True, Data->ChangeVertexFirstOnS2(), tolesp);
365 }
366 if(lin->EndPointOnSecond().NbPointOnRst() !=0){
367 ChFi3d_FilCommonPoint(lin->EndPointOnSecond(),lin->TransitionOnS2(),
368 Standard_False, Data->ChangeVertexLastOnS2(), tolesp);
369 }
370 done = CompleteData(Data,Func,lin,S1,S2,Or,0,0,0,0);
371 if(!done) throw Standard_Failure("PerformSurf : Failed approximation!");
372 // maybesingular = (Func.GetMinimalDistance()<=100*tolapp3d);
373 Standard_Boolean ok = Standard_False;
374 if(!Forward){
375 Intf = 0;
376 const ChFiDS_CommonPoint& cpf1 = Data->VertexFirstOnS1();
377 if(cpf1.IsOnArc()){
378 TopoDS_Face F1 = S1->Face();
379 TopoDS_Face bid;
380 Intf = !SearchFace(Spine,cpf1,F1,bid);
381 ok = Intf != 0;
382 }
383 const ChFiDS_CommonPoint& cpf2 = Data->VertexFirstOnS2();
384 if(cpf2.IsOnArc() && !ok){
385 TopoDS_Face F2 = S2->Face();
386 TopoDS_Face bid;
387 Intf = !SearchFace(Spine,cpf2,F2,bid);
388 }
389 }
390 Intl = 0;
391 ok = Standard_False;
392 const ChFiDS_CommonPoint& cpl1 = Data->VertexLastOnS1();
393 if(cpl1.IsOnArc()){
394 TopoDS_Face F1 = S1->Face();
395 TopoDS_Face bid;
396 Intl = !SearchFace(Spine,cpl1,F1,bid);
397 ok = Intl != 0;
398 }
399 const ChFiDS_CommonPoint& cpl2 = Data->VertexLastOnS2();
400 if(cpl2.IsOnArc() && !ok){
401 TopoDS_Face F2 = S2->Face();
402 TopoDS_Face bid;
403 Intl = !SearchFace(Spine,cpl2,F2,bid);
404 }
405 Data->FirstSpineParam(First);
406 Data->LastSpineParam (Last);
407
408 // if (maybesingular) SplitSurf(SeqData, lin);
409 // Necessite de trimer resultats : A faire
410 return Standard_True;
411 }
412
PerformSurf(ChFiDS_SequenceOfSurfData &,const Handle (ChFiDS_ElSpine)&,const Handle (ChFiDS_Spine)&,const Standard_Integer,const Handle (BRepAdaptor_Surface)&,const Handle (Adaptor3d_TopolTool)&,const Handle (BRepAdaptor_Curve2d)&,const Handle (BRepAdaptor_Surface)&,const Handle (BRepAdaptor_Curve2d)&,Standard_Boolean &,const Handle (BRepAdaptor_Surface)&,const Handle (Adaptor3d_TopolTool)&,const TopAbs_Orientation,const Standard_Real,const Standard_Real,const Standard_Real,Standard_Real &,Standard_Real &,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const math_Vector &)413 void FilletSurf_InternalBuilder::PerformSurf (ChFiDS_SequenceOfSurfData& , const Handle(ChFiDS_ElSpine)& , const Handle(ChFiDS_Spine)& , const Standard_Integer , const Handle(BRepAdaptor_Surface)& , const Handle(Adaptor3d_TopolTool)& , const Handle(BRepAdaptor_Curve2d)& , const Handle(BRepAdaptor_Surface)& , const Handle(BRepAdaptor_Curve2d)& , Standard_Boolean& , const Handle(BRepAdaptor_Surface)& , const Handle(Adaptor3d_TopolTool)& , const TopAbs_Orientation , const Standard_Real , const Standard_Real , const Standard_Real , Standard_Real& , Standard_Real& , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const math_Vector& )
414 {
415 throw Standard_DomainError("BlendFunc_CSConstRad::Section : Not implemented");
416 }
417
PerformSurf(ChFiDS_SequenceOfSurfData &,const Handle (ChFiDS_ElSpine)&,const Handle (ChFiDS_Spine)&,const Standard_Integer,const Handle (BRepAdaptor_Surface)&,const Handle (Adaptor3d_TopolTool)&,const TopAbs_Orientation,const Handle (BRepAdaptor_Surface)&,const Handle (Adaptor3d_TopolTool)&,const Handle (BRepAdaptor_Curve2d)&,const Handle (BRepAdaptor_Surface)&,const Handle (BRepAdaptor_Curve2d)&,Standard_Boolean &,const Standard_Real,const Standard_Real,const Standard_Real,Standard_Real &,Standard_Real &,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const math_Vector &)418 void FilletSurf_InternalBuilder::PerformSurf (ChFiDS_SequenceOfSurfData& , const Handle(ChFiDS_ElSpine)& , const Handle(ChFiDS_Spine)& , const Standard_Integer , const Handle(BRepAdaptor_Surface)& , const Handle(Adaptor3d_TopolTool)& , const TopAbs_Orientation , const Handle(BRepAdaptor_Surface)& , const Handle(Adaptor3d_TopolTool)& , const Handle(BRepAdaptor_Curve2d)& , const Handle(BRepAdaptor_Surface)& , const Handle(BRepAdaptor_Curve2d)& , Standard_Boolean& , const Standard_Real , const Standard_Real , const Standard_Real , Standard_Real& , Standard_Real& , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const math_Vector& )
419 {
420 throw Standard_DomainError("BlendFunc_CSConstRad::Section : Not implemented");
421 }
422
PerformSurf(ChFiDS_SequenceOfSurfData &,const Handle (ChFiDS_ElSpine)&,const Handle (ChFiDS_Spine)&,const Standard_Integer,const Handle (BRepAdaptor_Surface)&,const Handle (Adaptor3d_TopolTool)&,const Handle (BRepAdaptor_Curve2d)&,const Handle (BRepAdaptor_Surface)&,const Handle (BRepAdaptor_Curve2d)&,Standard_Boolean &,const TopAbs_Orientation,const Handle (BRepAdaptor_Surface)&,const Handle (Adaptor3d_TopolTool)&,const Handle (BRepAdaptor_Curve2d)&,const Handle (BRepAdaptor_Surface)&,const Handle (BRepAdaptor_Curve2d)&,Standard_Boolean &,const TopAbs_Orientation,const Standard_Real,const Standard_Real,const Standard_Real,Standard_Real &,Standard_Real &,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const Standard_Boolean,const math_Vector &)423 void FilletSurf_InternalBuilder::PerformSurf (ChFiDS_SequenceOfSurfData& , const Handle(ChFiDS_ElSpine)& , const Handle(ChFiDS_Spine)& , const Standard_Integer , const Handle(BRepAdaptor_Surface)& , const Handle(Adaptor3d_TopolTool)& , const Handle(BRepAdaptor_Curve2d)& , const Handle(BRepAdaptor_Surface)& , const Handle(BRepAdaptor_Curve2d)& , Standard_Boolean& , const TopAbs_Orientation , const Handle(BRepAdaptor_Surface)& , const Handle(Adaptor3d_TopolTool)& , const Handle(BRepAdaptor_Curve2d)& , const Handle(BRepAdaptor_Surface)& , const Handle(BRepAdaptor_Curve2d)& , Standard_Boolean& , const TopAbs_Orientation , const Standard_Real , const Standard_Real , const Standard_Real , Standard_Real& , Standard_Real& , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const Standard_Boolean , const math_Vector& )
424 {
425 throw Standard_DomainError("BlendFunc_CSConstRad::Section : Not implemented");
426 }
427
Done() const428 Standard_Boolean FilletSurf_InternalBuilder::Done() const
429 {
430 return done;
431 }
432 //=======================================================================
433 //function : NbSurface
434 //purpose : gives the number of NUBS surfaces of the Fillet
435 //=======================================================================
436
NbSurface() const437 Standard_Integer FilletSurf_InternalBuilder::NbSurface() const
438 {
439 return myListStripe.First()->SetOfSurfData()->Length();
440 }
441
442 //=======================================================================
443 //function : SurfaceFillet
444 //purpose : gives the NUBS surface of index Index
445 //=======================================================================
446
Handle(Geom_Surface)447 const Handle(Geom_Surface)& FilletSurf_InternalBuilder::SurfaceFillet(const Standard_Integer Index) const
448 {
449 Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->Surf();
450
451 return myDS->Surface(isurf).Surface();
452 }
453
454 //=======================================================================
455 //function : TolApp3d
456 //purpose : gives the 3d tolerance reached during approximation
457 //=======================================================================
458
TolApp3d(const Standard_Integer Index) const459 Standard_Real FilletSurf_InternalBuilder::TolApp3d(const Standard_Integer Index) const
460 {
461 Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->Surf();
462
463 return myDS->Surface(isurf).Tolerance();
464 }
465
466 //=======================================================================
467 //function : SupportFace1
468 //purpose : gives the first support face relative to SurfaceFillet(Index)
469 //=======================================================================
SupportFace1(const Standard_Integer Index) const470 const TopoDS_Face& FilletSurf_InternalBuilder::SupportFace1(const Standard_Integer Index) const
471 {
472 Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->IndexOfS1();
473
474 return TopoDS::Face(myDS->Shape(isurf));
475 }
476 //=======================================================================
477 //function : SupportFace2
478 //purpose : gives the second support face relative to SurfaceFillet(Index)
479 //=======================================================================
SupportFace2(const Standard_Integer Index) const480 const TopoDS_Face& FilletSurf_InternalBuilder::SupportFace2(const Standard_Integer Index) const
481 {
482 Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->IndexOfS2();
483
484 return TopoDS::Face(myDS->Shape(isurf));
485 }
486 //===============================================================================
487 //function : CurveOnSup1
488 //purpose : gives the 3d curve of SurfaceFillet(Index) on SupportFace1(Index)
489 //===============================================================================
490
Handle(Geom_Curve)491 const Handle(Geom_Curve)& FilletSurf_InternalBuilder::CurveOnFace1(const Standard_Integer Index) const
492 {
493 Standard_Integer icurv = myListStripe.First()->SetOfSurfData()->Value(Index)->
494 InterferenceOnS1().LineIndex();
495 return myDS->Curve(icurv).Curve();
496 }
497
498 //=======================================================================
499 //function : CurveOnSup2
500 //purpose : gives the 3d curve of SurfaceFillet(Index) on SupportFace2(Index
501 //=======================================================================
502
Handle(Geom_Curve)503 const Handle(Geom_Curve)& FilletSurf_InternalBuilder::CurveOnFace2(const Standard_Integer Index) const
504 {
505 Standard_Integer icurv=myListStripe.First()->SetOfSurfData()->Value(Index)->
506 InterferenceOnS2().LineIndex();
507 return myDS->Curve(icurv).Curve();
508 }
509 //=======================================================================
510 //function : PCurveOnFace1
511 //purpose : gives the PCurve associated to CurvOnSup1(Index) on the support face
512 //=======================================================================
Handle(Geom2d_Curve)513 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurveOnFace1(const Standard_Integer Index) const
514 {
515 return myListStripe.First()->SetOfSurfData()->Value(Index)->
516 InterferenceOnS1().PCurveOnFace();
517
518 }
519
520 //=======================================================================
521 //function : PCurveOnFillet1
522 //purpose : gives the PCurve associated to CurveOnFace1(Index) on the Fillet
523 //=======================================================================
524
Handle(Geom2d_Curve)525 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurve1OnFillet(const Standard_Integer Index) const
526 {return myListStripe.First()->SetOfSurfData()->Value(Index)->
527 InterferenceOnS1().PCurveOnSurf();
528
529 }
530
531 //=======================================================================
532 //function : PCurveOnFace2
533 //purpose : gives the PCurve associated to CurvOnSup2(Index) on the support face
534 //=======================================================================
Handle(Geom2d_Curve)535 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurveOnFace2(const Standard_Integer Index) const
536 {
537 return myListStripe.First()->SetOfSurfData()->Value(Index)->
538 InterferenceOnS2().PCurveOnFace();
539 }
540
541 //=======================================================================
542 //function : PCurveOnFillet2
543 //purpose : gives the PCurve associated to CurveOnFace2(Index) on the Fillet
544 //=======================================================================
Handle(Geom2d_Curve)545 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurve2OnFillet(const Standard_Integer Index) const
546 {return myListStripe.First()->SetOfSurfData()->Value(Index)->
547 InterferenceOnS2().PCurveOnSurf();
548
549 }
550
551 //=======================================================================
552 //function : FirstParameter
553 //purpose : gives the parameter of the fillet on the first edge
554 //=======================================================================
555
FirstParameter() const556 Standard_Real FilletSurf_InternalBuilder::FirstParameter() const
557 {
558 const Handle(ChFiDS_Stripe)& st = myListStripe.First();
559 const Handle(ChFiDS_Spine)& sp = st->Spine();
560 const Handle(ChFiDS_SurfData)& sd = st->SetOfSurfData()->Value(1);
561 Standard_Real p = sd->FirstSpineParam();
562 Standard_Integer ind = 1;
563 if(sp->IsPeriodic()) ind = sp->Index(p);
564 Standard_Real ep;
565 if(ComputeEdgeParameter(sp,ind,p,ep,tolesp)) return ep;
566 return 0.0;
567 }
568 //=======================================================================
569 //function : LastParameter
570 //purpose : gives the parameter of the fillet on the last edge
571 //=======================================================================
LastParameter() const572 Standard_Real FilletSurf_InternalBuilder::LastParameter() const
573 {
574 const Handle(ChFiDS_Stripe)& st = myListStripe.First();
575 const Handle(ChFiDS_Spine)& sp = st->Spine();
576 const Handle(ChFiDS_SurfData)& sd = st->SetOfSurfData()->Value(NbSurface());
577 Standard_Real p = sd->LastSpineParam();
578 Standard_Integer ind = sp->NbEdges();
579 if(sp->IsPeriodic()) ind = sp->Index(p);
580 Standard_Real ep;
581 if(ComputeEdgeParameter(sp,ind,p,ep,tolesp)) return ep;
582 return 0.0;
583 }
584
585 //=======================================================================
586 //function : StatusStartSection
587 //purpose : returns:
588 // twoExtremityonEdge: each extremity of start section of the Fillet is
589 // on the edge of the corresponding support face.
590 // OneExtremityOnEdge: only one of the extremities of start section of the Fillet
591 // is on the edge of the corresponding support face.
592 // NoExtremityOnEdge: any extremity of start section of the fillet is on
593 // the edge of the corresponding support face.
594 //=======================================================================
595
596
StartSectionStatus() const597 FilletSurf_StatusType FilletSurf_InternalBuilder::StartSectionStatus() const
598 {
599
600 Standard_Boolean isonedge1 = myListStripe.First()->SetOfSurfData()->Value(1)->
601 VertexFirstOnS1().IsOnArc();
602 Standard_Boolean isonedge2= myListStripe.First()->SetOfSurfData()->Value(1)->
603 VertexFirstOnS2().IsOnArc();
604
605 if (isonedge1 && isonedge2)
606 {return FilletSurf_TwoExtremityOnEdge;}
607 else if ((!isonedge1)&& (!isonedge2))
608 {return FilletSurf_NoExtremityOnEdge ;}
609 else
610 {return FilletSurf_OneExtremityOnEdge;}
611 }
612 //=======================================================================
613 //function : StatusEndSection()
614 //purpose : returns:
615 // twoExtremityonEdge: each extremity of end section of the Fillet is
616 // on the edge of the corresponding support face.
617 // OneExtremityOnEdge: only one of the extremities of end section of the Fillet
618 // is on the edge of the corresponding support face.
619 // NoExtremityOnEdge: any extremity of end section of the fillet is on
620 // the edge of the corresponding support face.
621 //=======================================================================
EndSectionStatus() const622 FilletSurf_StatusType FilletSurf_InternalBuilder::EndSectionStatus() const
623 {
624 Standard_Boolean isonedge1 = myListStripe.First()->SetOfSurfData()->Value(NbSurface())->
625 VertexLastOnS1().IsOnArc();
626 Standard_Boolean isonedge2= myListStripe.First()->SetOfSurfData()->Value(NbSurface())->
627 VertexLastOnS2().IsOnArc();
628
629 if (isonedge1 && isonedge2)
630 { return FilletSurf_TwoExtremityOnEdge;}
631 else if ((!isonedge1)&& (!isonedge2))
632 { return FilletSurf_NoExtremityOnEdge;}
633 else
634 { return FilletSurf_OneExtremityOnEdge;}
635 }
636
637
638 //=======================================================================
639 //function : Simulate
640 //purpose : computes only the sections used in the computation of the fillet
641 //=======================================================================
Simulate()642 void FilletSurf_InternalBuilder::Simulate()
643 {
644 //ChFi3d_FilBuilder::Simulate(1);
645 Handle(ChFiDS_Stripe) Stripe = myListStripe.First();
646 Handle(ChFiDS_HData)& HData = Stripe->ChangeSetOfSurfData();
647 HData = new ChFiDS_HData();
648 Handle(ChFiDS_Spine)& Spine = Stripe->ChangeSpine();
649 TopAbs_Orientation RefOr1,RefOr2;
650 Standard_Integer RefChoix;
651 StripeOrientations(Spine,RefOr1,RefOr2,RefChoix);
652 Stripe->OrientationOnFace1(RefOr1);
653 Stripe->OrientationOnFace2(RefOr2);
654 Stripe->Choix(RefChoix);
655 PerformSetOfKGen(Stripe,1);
656 }
657
658
659 //=======================================================================
660 //function : NbSection
661 //purpose : gives the number of sections relative to SurfaceFillet(IndexSurf)
662 //=======================================================================
NbSection(const Standard_Integer IndexSurf) const663 Standard_Integer FilletSurf_InternalBuilder::NbSection(const Standard_Integer IndexSurf) const
664 {
665 return Sect(1,IndexSurf)->Length();
666 }
667
668 //=======================================================================
669 //function : Section
670 //purpose : gives the arc of circle corresponding to section number
671 // IndexSec of SurfaceFillet(IndexSurf) (The basis curve of the
672 // trimmed curve is a Geom_Circle)
673 //=======================================================================
Section(const Standard_Integer IndexSurf,const Standard_Integer IndexSec,Handle (Geom_TrimmedCurve)& Circ) const674 void FilletSurf_InternalBuilder::Section(const Standard_Integer IndexSurf, const Standard_Integer IndexSec, Handle(Geom_TrimmedCurve)& Circ) const
675 {
676 gp_Circ c;
677 Standard_Real deb,fin;
678 Sect(1,IndexSurf)->Value(IndexSec).Get(c,deb,fin);
679 Handle(Geom_Circle) Gc = new Geom_Circle(c);
680 Circ= new Geom_TrimmedCurve(Gc,deb,fin);
681 }
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697