1 # -*- coding: utf-8 -*-
5 from geomsmesh import geompy
6 from findWireEndVertices import findWireEndVertices
7 from prolongeWire import prolongeWire
9 def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe):
11 pipe de fond de fissure, prolongé, partition face fissure par pipe
12 identification des edges communes pipe et face fissure
15 shapeDefaut = shapesFissure[0] # faces de fissure, débordant
16 fondFiss = shapesFissure[4] # groupe d'edges de fond de fissure
17 centreFondFiss = elementsDefaut[15]
18 tgtCentre = elementsDefaut[16]
20 if geompy.NumberOfFaces(shapeDefaut) == 1:
21 plan = geompy.MakePlane(centreFondFiss, tgtCentre, 10000)
22 shapeDefaut = geompy.MakePartition([shapeDefaut], [plan], [], [], geompy.ShapeType["FACE"], 0, [], 0)
23 #fondFissCoupe = geompy.GetInPlaceByHistory(shapeDefaut, fondFiss) #= inutile
24 geompy.addToStudy(shapeDefaut, 'shapeDefaut_coupe')
25 #geompy.addToStudyInFather(shapeDefaut, fondFissCoupe, 'fondFiss_coupe')
27 extrem, norms = findWireEndVertices(fondFiss, True)
28 logging.debug("extrem: %s, norm: %s",extrem, norms)
29 cercle = geompy.MakeCircle(extrem[0], norms[0], rayonPipe)
30 cercle = geompy.MakeRotation(cercle, norms[0], math.pi/3.0 ) # éviter d'avoir l'arête de couture du pipe presque confondue avec la face fissure
31 geompy.addToStudy(cercle, 'cercle')
32 fondFissProlonge = prolongeWire(fondFiss, extrem, norms, 2*rayonPipe)
33 pipeFiss = geompy.MakePipe(cercle, fondFissProlonge)
34 geompy.addToStudy(pipeFiss, 'pipeFiss')
35 partFissPipe = geompy.MakePartition([shapeDefaut, pipeFiss], [], [], [], geompy.ShapeType["FACE"], 0, [], 1)
36 geompy.addToStudy(partFissPipe, 'partFissPipe')
37 fissPipe = geompy.GetInPlaceByHistory(partFissPipe, shapeDefaut)
38 geompy.addToStudy(fissPipe, 'fissPipe')
39 partPipe = geompy.GetInPlaceByHistory(partFissPipe, pipeFiss)
40 geompy.addToStudy(partPipe, 'partPipe')
42 edgesPipeFiss = geompy.GetSharedShapesMulti([fissPipe, partPipe], geompy.ShapeType["EDGE"])
43 for i, edge in enumerate(edgesPipeFiss):
45 geompy.addToStudyInFather(fissPipe, edge, name)
47 wirePipeFiss = geompy.MakeWire(edgesPipeFiss)
49 wirePipeFiss = geompy.MakeCompound(edgesPipeFiss)
50 logging.debug("wirePipeFiss construit sous forme de compound")
51 geompy.addToStudy(wirePipeFiss, "wirePipeFiss")
53 wireFondFiss = geompy.GetInPlace(partFissPipe,fondFiss)
54 edgesFondFiss = geompy.GetSharedShapesMulti([fissPipe, wireFondFiss], geompy.ShapeType["EDGE"])
55 for i, edge in enumerate(edgesFondFiss):
56 name = "edgeFondFiss%d"%i
57 geompy.addToStudyInFather(fissPipe, edge, name)
58 wireFondFiss = geompy.MakeWire(edgesFondFiss)
59 geompy.addToStudy(wireFondFiss,"wireFondFiss")
61 return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss)