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 from .fissError import fissError
14 def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe):
16 pipe de fond de fissure, prolongé, partition face fissure par pipe
17 identification des edges communes pipe et face fissure
20 shapeDefaut = shapesFissure[0] # faces de fissure, débordant
21 fondFiss = shapesFissure[4] # groupe d'edges de fond de fissure
22 centreFondFiss = elementsDefaut[15]
23 tgtCentre = elementsDefaut[16]
25 if geompy.NumberOfFaces(shapeDefaut) == 1:
26 plan = geompy.MakePlane(centreFondFiss, tgtCentre, 10000)
27 shapeDefaut = geompy.MakePartition([shapeDefaut], [plan], [], [], geompy.ShapeType["FACE"], 0, [], 0)
28 #fondFissCoupe = geompy.GetInPlaceByHistory(shapeDefaut, fondFiss) #= inutile
29 geomPublish(initLog.debug, shapeDefaut, 'shapeDefaut_coupe')
30 #geomPublishInFather(initLog.debug,shapeDefaut, fondFissCoupe, 'fondFiss_coupe')
32 extrem, norms = findWireEndVertices(fondFiss, True)
33 logging.debug("extrem: %s, norm: %s",extrem, norms)
34 cercle = geompy.MakeCircle(extrem[0], norms[0], rayonPipe)
35 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
36 geomPublish(initLog.debug, cercle, 'cercle')
37 fondFissProlonge = prolongeWire(fondFiss, extrem, norms, 2*rayonPipe)
39 pipeFiss = geompy.MakePipe(cercle, fondFissProlonge)
41 texte = "génération du pipe le long de la ligne de fond de fissure prolongée impossible. "
42 texte += "Cause possible : la ligne s'autointersecte lorsqu'on la prolonge."
43 raise fissError(traceback.extract_stack(),texte)
44 geomPublish(initLog.debug, pipeFiss, 'pipeFiss')
45 partFissPipe = geompy.MakePartition([shapeDefaut, pipeFiss], [], [], [], geompy.ShapeType["FACE"], 0, [], 1)
46 geomPublish(initLog.debug, partFissPipe, 'partFissPipe')
47 fissPipe = geompy.GetInPlaceByHistory(partFissPipe, shapeDefaut)
48 geomPublish(initLog.debug, fissPipe, 'fissPipe')
49 partPipe = geompy.GetInPlaceByHistory(partFissPipe, pipeFiss)
50 geomPublish(initLog.debug, partPipe, 'partPipe')
52 edgesPipeFiss = geompy.GetSharedShapesMulti([fissPipe, partPipe], geompy.ShapeType["EDGE"])
53 for i, edge in enumerate(edgesPipeFiss):
55 geomPublishInFather(initLog.debug,fissPipe, edge, name)
57 wirePipeFiss = geompy.MakeWire(edgesPipeFiss)
59 wirePipeFiss = geompy.MakeCompound(edgesPipeFiss)
60 logging.debug("wirePipeFiss construit sous forme de compound")
61 geomPublish(initLog.debug, wirePipeFiss, "wirePipeFiss")
63 wireFondFiss = geompy.GetInPlace(partFissPipe,fondFiss)
64 edgesFondFiss = geompy.GetSharedShapesMulti([fissPipe, wireFondFiss], geompy.ShapeType["EDGE"])
65 for i, edge in enumerate(edgesFondFiss):
66 name = "edgeFondFiss%d"%i
67 geomPublishInFather(initLog.debug,fissPipe, edge, name)
68 wireFondFiss = geompy.MakeWire(edgesFondFiss)
69 geomPublish(initLog.debug, wireFondFiss,"wireFondFiss")
71 return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss)