mirror of
https://git.salome-platform.org/gitpub/modules/smesh.git
synced 2025-04-01 14:34:30 +05:00
97 lines
3.5 KiB
Python
97 lines
3.5 KiB
Python
# -*- coding: utf-8 -*-
|
|
|
|
import logging
|
|
from geomsmesh import smesh
|
|
import SMESH
|
|
import SALOMEDS
|
|
|
|
from creeZoneDefautMaillage import creeZoneDefautMaillage
|
|
from peauInterne import peauInterne
|
|
from quadranglesToShapeNoCorner import quadranglesToShapeNoCorner
|
|
from creeZoneDefautFilling import creeZoneDefautFilling
|
|
from creeZoneDefautGeom import creeZoneDefautGeom
|
|
from getCentreFondFiss import getCentreFondFiss
|
|
|
|
# -----------------------------------------------------------------------------
|
|
# ---
|
|
|
|
def creeZoneDefautDansObjetSain(geometriesSaines, maillagesSains, shapesFissure,
|
|
shapeFissureParams, maillageFissureParams):
|
|
"""
|
|
#TODO: a compléter
|
|
"""
|
|
logging.info('start')
|
|
|
|
#smesh.SetCurrentStudy(salome.myStudy)
|
|
|
|
geometrieSaine = geometriesSaines[0]
|
|
maillageSain = maillagesSains[0]
|
|
isHexa = maillagesSains[1]
|
|
shapeDefaut = shapesFissure[0]
|
|
tailleDefaut = shapesFissure[2]
|
|
coordsNoeudsFissure = shapesFissure[3]
|
|
|
|
isElliptique = False
|
|
if shapeFissureParams.has_key('elliptique'):
|
|
isElliptique = shapeFissureParams['elliptique']
|
|
if isElliptique:
|
|
if shapeFissureParams.has_key('demiGrandAxe'):
|
|
demiGrandAxe = shapeFissureParams['demiGrandAxe']
|
|
else:
|
|
demiGrandAxe = shapeFissureParams['longueur']
|
|
lgExtrusion = 2.0*demiGrandAxe
|
|
else:
|
|
lgExtrusion = 50.
|
|
|
|
nomRep = maillageFissureParams['nomRep']
|
|
nomFicSain = maillageFissureParams['nomFicSain']
|
|
|
|
fichierMaillageSain = nomRep + '/' + nomFicSain + '.med'
|
|
|
|
# --- centre de fond de fissure et tangente
|
|
|
|
edgeFondExt, centreFondFiss, tgtCentre = getCentreFondFiss(shapesFissure)
|
|
|
|
|
|
# --- zone de défaut
|
|
nomZones = "zoneDefaut"
|
|
|
|
[origShapes, verticesShapes, dmoyen] = \
|
|
creeZoneDefautMaillage(maillagesSains, shapeDefaut, tailleDefaut, nomZones, coordsNoeudsFissure)
|
|
|
|
maillageSain.ExportMED( fichierMaillageSain, 0, SMESH.MED_V2_2, 1 )
|
|
logging.debug("fichier maillage sain %s", fichierMaillageSain)
|
|
[maillageSain, internalBoundary, zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges] = \
|
|
peauInterne(fichierMaillageSain, shapeDefaut, nomZones)
|
|
|
|
facesDefaut = []
|
|
centresDefaut = []
|
|
normalsDefaut =[]
|
|
extrusionsDefaut = []
|
|
isPlane = False
|
|
if isHexa and not isPlane:
|
|
meshQuad = smesh.CopyMesh( zoneDefaut_skin, 'meshQuad', 0, 0)
|
|
|
|
fillings, noeuds_bords, bordsPartages, fillconts, idFilToCont = quadranglesToShapeNoCorner(meshQuad, shapeFissureParams, centreFondFiss)
|
|
|
|
for filling in fillings:
|
|
[faceDefaut, centreDefaut, normalDefaut, extrusionDefaut] = \
|
|
creeZoneDefautFilling(filling, shapeDefaut, lgExtrusion)
|
|
facesDefaut.append(faceDefaut)
|
|
centresDefaut.append(centreDefaut)
|
|
normalsDefaut.append(normalDefaut)
|
|
extrusionsDefaut.append(extrusionDefaut)
|
|
else:
|
|
[facesDefaut, centreDefaut, normalDefaut, extrusionDefaut] = \
|
|
creeZoneDefautGeom( geometrieSaine, shapeDefaut, origShapes, verticesShapes, dmoyen, lgExtrusion)
|
|
bordsPartages = []
|
|
for face in facesDefaut:
|
|
bordsPartages.append([None,None]) # TODO : traitement des arêtes vives ?
|
|
fillconts = facesDefaut
|
|
idFilToCont = range(len(facesDefaut))
|
|
|
|
return [facesDefaut, centresDefaut, normalsDefaut, extrusionsDefaut, dmoyen, bordsPartages, fillconts, idFilToCont,
|
|
maillageSain, internalBoundary, zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges,
|
|
edgeFondExt, centreFondFiss, tgtCentre]
|
|
|