1 # -*- coding: utf-8 -*-
5 from geomsmesh import geompy
6 from geomsmesh import geomPublish
7 from geomsmesh import geomPublishInFather
9 from findWireEndVertices import findWireEndVertices
10 from prolongeWire import prolongeWire
12 def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe):
14 pipe de fond de fissure, prolongé, partition face fissure par pipe
15 identification des edges communes pipe et face fissure
18 shapeDefaut = shapesFissure[0] # faces de fissure, débordant
19 fondFiss = shapesFissure[4] # groupe d'edges de fond de fissure
20 centreFondFiss = elementsDefaut[15]
21 tgtCentre = elementsDefaut[16]
23 if geompy.NumberOfFaces(shapeDefaut) == 1:
24 plan = geompy.MakePlane(centreFondFiss, tgtCentre, 10000)
25 shapeDefaut = geompy.MakePartition([shapeDefaut], [plan], [], [], geompy.ShapeType["FACE"], 0, [], 0)
26 #fondFissCoupe = geompy.GetInPlaceByHistory(shapeDefaut, fondFiss) #= inutile
27 geomPublish(initLog.debug, shapeDefaut, 'shapeDefaut_coupe')
28 #geomPublishInFather(initLog.debug,shapeDefaut, fondFissCoupe, 'fondFiss_coupe')
30 extrem, norms = findWireEndVertices(fondFiss, True)
31 logging.debug("extrem: %s, norm: %s",extrem, norms)
32 cercle = geompy.MakeCircle(extrem[0], norms[0], rayonPipe)
33 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
34 geomPublish(initLog.debug, cercle, 'cercle')
35 fondFissProlonge = prolongeWire(fondFiss, extrem, norms, 2*rayonPipe)
36 pipeFiss = geompy.MakePipe(cercle, fondFissProlonge)
37 geomPublish(initLog.debug, pipeFiss, 'pipeFiss')
38 partFissPipe = geompy.MakePartition([shapeDefaut, pipeFiss], [], [], [], geompy.ShapeType["FACE"], 0, [], 1)
39 geomPublish(initLog.debug, partFissPipe, 'partFissPipe')
40 fissPipe = geompy.GetInPlaceByHistory(partFissPipe, shapeDefaut)
41 geomPublish(initLog.debug, fissPipe, 'fissPipe')
42 partPipe = geompy.GetInPlaceByHistory(partFissPipe, pipeFiss)
43 geomPublish(initLog.debug, partPipe, 'partPipe')
45 edgesPipeFiss = geompy.GetSharedShapesMulti([fissPipe, partPipe], geompy.ShapeType["EDGE"])
46 for i, edge in enumerate(edgesPipeFiss):
48 geomPublishInFather(initLog.debug,fissPipe, edge, name)
50 wirePipeFiss = geompy.MakeWire(edgesPipeFiss)
52 wirePipeFiss = geompy.MakeCompound(edgesPipeFiss)
53 logging.debug("wirePipeFiss construit sous forme de compound")
54 geomPublish(initLog.debug, wirePipeFiss, "wirePipeFiss")
56 wireFondFiss = geompy.GetInPlace(partFissPipe,fondFiss)
57 edgesFondFiss = geompy.GetSharedShapesMulti([fissPipe, wireFondFiss], geompy.ShapeType["EDGE"])
58 for i, edge in enumerate(edgesFondFiss):
59 name = "edgeFondFiss%d"%i
60 geomPublishInFather(initLog.debug,fissPipe, edge, name)
61 wireFondFiss = geompy.MakeWire(edgesFondFiss)
62 geomPublish(initLog.debug, wireFondFiss,"wireFondFiss")
64 return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss)