This commit is contained in:
Paul RASCLE 2014-09-21 21:49:53 +02:00
parent 106415d543
commit 249a5808d8
9 changed files with 21 additions and 7 deletions

View File

@ -131,7 +131,7 @@ runall = True
if runall: if runall:
torun = [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,] torun = [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,]
else: #prob 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27 else: #prob 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27
torun = [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,] torun = [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,]
for i in range(len(problemes)): for i in range(len(problemes)):
if torun[i]: if torun[i]:

View File

@ -3,6 +3,7 @@
import logging import logging
from geomsmesh import geompy from geomsmesh import geompy
import GEOM
from sortEdges import sortEdges from sortEdges import sortEdges
def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks, def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks,
@ -48,6 +49,9 @@ def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks,
for i, nappes in enumerate(listNappes): for i, nappes in enumerate(listNappes):
if facesDebouchantes[i]: if facesDebouchantes[i]:
for k, face in enumerate(facesPipePeau): for k, face in enumerate(facesPipePeau):
#logging.debug('i, k, face, nappes[0] %s %s %s %s', i, k, face, nappes[0])
#geompy.addToStudy(nappes[0], 'lanappe')
#geompy.addToStudy(face, 'laface')
edge = geompy.MakeSection(face, nappes[0]) edge = geompy.MakeSection(face, nappes[0])
if geompy.NbShapes(edge, geompy.ShapeType["EDGE"]) > 0: if geompy.NbShapes(edge, geompy.ShapeType["EDGE"]) > 0:
idFacesDebouchantes[i] = k idFacesDebouchantes[i] = k

View File

@ -154,8 +154,8 @@ def construitFissureGenerale(maillagesSains,
nbFacesFilling = len(partitionsPeauFissFond) nbFacesFilling = len(partitionsPeauFissFond)
ptEdgeFond = [ [] for i in range(nbFacesFilling)] # pour chaque face [points edge fond de fissure aux débouchés du pipe] ptEdgeFond = [ [] for i in range(nbFacesFilling)] # pour chaque face [points edge fond de fissure aux débouchés du pipe]
facesPipePeau = [ [] for i in range(nbFacesFilling)] # pour chaque face [faces du pipe débouchantes] fsPipePeau = [ [] for i in range(nbFacesFilling)] # pour chaque face [faces du pipe débouchantes]
edgeRadFacePipePeau = [ [] for i in range(nbFacesFilling)] # pour chaque face [edge radiale des faces du pipe débouchantes ] edRadFPiPo = [ [] for i in range(nbFacesFilling)] # pour chaque face [edge radiale des faces du pipe débouchantes ]
fsFissuExt = [ [] for i in range(nbFacesFilling)] # pour chaque face [faces de fissure externes au pipe] fsFissuExt = [ [] for i in range(nbFacesFilling)] # pour chaque face [faces de fissure externes au pipe]
edFisExtPe = [ [] for i in range(nbFacesFilling)] # pour chaque face [edge en peau des faces de fissure externes (pas subshape facePeau)] edFisExtPe = [ [] for i in range(nbFacesFilling)] # pour chaque face [edge en peau des faces de fissure externes (pas subshape facePeau)]
edFisExtPi = [ [] for i in range(nbFacesFilling)] # pour chaque face [edge commun au pipe des faces de fissure externes] edFisExtPi = [ [] for i in range(nbFacesFilling)] # pour chaque face [edge commun au pipe des faces de fissure externes]
@ -174,8 +174,8 @@ def construitFissureGenerale(maillagesSains,
facesDefaut, centreFondFiss, rayonPipe, facesDefaut, centreFondFiss, rayonPipe,
aretesVivesCoupees) aretesVivesCoupees)
ptEdgeFond[ifil] = dataPPFF['endsEdgeFond'] ptEdgeFond[ifil] = dataPPFF['endsEdgeFond']
facesPipePeau[ifil] = dataPPFF['facesPipePeau'] fsPipePeau[ifil] = dataPPFF['facesPipePeau']
edgeRadFacePipePeau[ifil] = dataPPFF['edgeRadFacePipePeau'] edRadFPiPo[ifil] = dataPPFF['edgeRadFacePipePeau']
fsFissuExt[ifil] = dataPPFF['facesFissExt'] fsFissuExt[ifil] = dataPPFF['facesFissExt']
edFisExtPe[ifil] = dataPPFF['edgesFissExtPeau'] edFisExtPe[ifil] = dataPPFF['edgesFissExtPeau']
edFisExtPi[ifil] = dataPPFF['edgesFissExtPipe'] edFisExtPi[ifil] = dataPPFF['edgesFissExtPipe']
@ -187,6 +187,11 @@ def construitFissureGenerale(maillagesSains,
edFissPeau[ifil] = dataPPFF['edgesFissurePeau'] edFissPeau[ifil] = dataPPFF['edgesFissurePeau']
ptFisExtPi[ifil] = dataPPFF['verticesPipePeau'] ptFisExtPi[ifil] = dataPPFF['verticesPipePeau']
facesPipePeau = []
edgeRadFacePipePeau = []
for ifil in range(nbFacesFilling):
facesPipePeau += fsPipePeau[ifil]
edgeRadFacePipePeau += edRadFPiPo[ifil]
for i, avc in enumerate(aretesVivesCoupees): for i, avc in enumerate(aretesVivesCoupees):
name = "areteViveCoupee%d"%i name = "areteViveCoupee%d"%i

View File

@ -4,6 +4,7 @@ import logging
from geomsmesh import geompy from geomsmesh import geompy
from geomsmesh import smesh from geomsmesh import smesh
import SMESH
def construitMaillagePipe(gptsdisks, idisklim, nbsegCercle, nbsegRad): def construitMaillagePipe(gptsdisks, idisklim, nbsegCercle, nbsegRad):
""" """

View File

@ -54,7 +54,7 @@ def identifieElementsGeometriquesPeau(ifil, partitionPeauFissFond, edgesPipeFiss
edgesPeauFondIn, fillingFaceExterne, aretesVivesC, aretesVivesCoupees) edgesPeauFondIn, fillingFaceExterne, aretesVivesC, aretesVivesCoupees)
dataPPFF = dict(endsEdgeFond = endsEdgeFond, # pour chaque face [points edge fond de fissure aux débouchés du pipe] dataPPFF = dict(endsEdgeFond = endsEdgeFond, # pour chaque face [points edge fond de fissure aux débouchés du pipe]
facesPipePeau = facesPipePeau, # pour chaque face [face facesPipePeau = facesPipePeau, # pour chaque face [faces du pipe débouchantes]
edgeRadFacePipePeau = edgeRadFacePipePeau, # pour chaque face [edge radiale des faces du pipe débouchantes ] edgeRadFacePipePeau = edgeRadFacePipePeau, # pour chaque face [edge radiale des faces du pipe débouchantes ]
facesFissExt = facesFissExt, # pour chaque face [faces de fissure externes au pipe] facesFissExt = facesFissExt, # pour chaque face [faces de fissure externes au pipe]
edgesFissExtPeau = edgesFissExtPeau, # pour chaque face [edge en peau des faces de fissure externes (pas subshape facePeau)] edgesFissExtPeau = edgesFissExtPeau, # pour chaque face [edge en peau des faces de fissure externes (pas subshape facePeau)]

View File

@ -13,7 +13,7 @@ def identifieFacesEdgesFissureExterne(fsFissuExt, edFisExtPe, edFisExtPi, edgesP
facesFissExt = [] facesFissExt = []
edgesFissExtPeau = [] edgesFissExtPeau = []
edgesFissExtPipe = [] edgesFissExtPipe = []
for ifil in len(fsFissuExt): # TODO: éliminer les doublons (comparer tous les vertices triés, avec mesure de distance ?) for ifil in range(len(fsFissuExt)): # TODO: éliminer les doublons (comparer tous les vertices triés, avec mesure de distance ?)
facesFissExt += fsFissuExt[ifil] facesFissExt += fsFissuExt[ifil]
edgesFissExtPeau += edFisExtPe[ifil] edgesFissExtPeau += edFisExtPe[ifil]
edgesFissExtPipe += edFisExtPi[ifil] edgesFissExtPipe += edFisExtPi[ifil]

View File

@ -4,6 +4,7 @@ import logging
from geomsmesh import geompy from geomsmesh import geompy
from geomsmesh import smesh from geomsmesh import smesh
import SMESH
from putName import putName from putName import putName
@ -30,6 +31,7 @@ def mailleAretesEtJonction(internalBoundary, aretesVivesCoupees, dmoyen):
# --- maillage des éventuelles arêtes vives entre faces reconstruites # --- maillage des éventuelles arêtes vives entre faces reconstruites
grpAretesVives = None
if len(aretesVivesCoupees) > 0: if len(aretesVivesCoupees) > 0:
aretesVivesC = geompy.MakeCompound(aretesVivesCoupees) aretesVivesC = geompy.MakeCompound(aretesVivesCoupees)
meshAretesVives = smesh.Mesh(aretesVivesC) meshAretesVives = smesh.Mesh(aretesVivesC)

View File

@ -5,6 +5,7 @@ import logging
from geomsmesh import geompy from geomsmesh import geompy
from geomsmesh import smesh from geomsmesh import smesh
from salome.smesh import smeshBuilder from salome.smesh import smeshBuilder
import SMESH
from putName import putName from putName import putName

View File

@ -5,6 +5,7 @@ import logging
from geomsmesh import geompy from geomsmesh import geompy
from geomsmesh import smesh from geomsmesh import smesh
from salome.smesh import smeshBuilder from salome.smesh import smeshBuilder
import SMESH
from putName import putName from putName import putName