From: Pascale Noyret Date: Fri, 17 Nov 2006 15:33:17 +0000 (+0000) Subject: PN: version mergée entre AMA et eltStructure X-Git-Tag: mergefrom_V3_2_0_maintenance_18nov06 X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=3277807704ca065ae8b46cf81295e1d5bb832748;p=modules%2Feficas.git PN: version mergée entre AMA et eltStructure --- diff --git a/src/EFICASGUI/cabri.py b/src/EFICASGUI/cabri.py new file mode 100644 index 00000000..71406b38 --- /dev/null +++ b/src/EFICASGUI/cabri.py @@ -0,0 +1,737 @@ +# -*- coding: iso-8859-1 -*- +""" +outil métier Cabri pour Salome + +""" + + +def tetra( name, **param ): + """ + construction géométrie Cabri en Tetra + """ + #************************************************************* + # Paramètres physiques + #************************************************************* + # Nombre déléments circonférentiels (NBR_CIR) + ncir = param['ncir'] + # Temps danalyse + temps = param['temps'] + # Nombre déléments de lalésage (NBR_ALE) + nsect = param['nsect'] + # Nombre déléments radiaux (NBR_RAD) + nrad = param['nrad'] + # Nombre déléments verticaux (NBR_VER) + nver = param['nver'] + + + # Epaisseur de la rondelle (GOUJ_E_RONDEL) + he = param['he'] + # Hauteur de la bride (BRID_H) + hc1 = param['hc1'] + # Epaisseur de lécrou (GOUJ_E_ECROU) + e = param['e'] + # Nombre de goujons de la jonction boulonnée (GOUJ_N_GOUJON) + nbgouj = param['nbgouj'] + # Diamètre des goujons de la jonction boulonnée (GOUJ_D_GOUJON) + dg = param['dg'] + # Hauteur des alésages de la bride permettant lintroduction des goujons (BRID_H_ALESAG) + htrou = param['htrou'] + # Epaisseur du joint au niveau de linterface détanchéité (ETAN_E_JOINT) + j = param['j'] + # Hauteur de lépaulement de la bride au niveau de linterface détanchéité (BRID_H_EPAUL) + hb = param['hb'] + # Rayon du congé de la bride (BRID_R_CONGE) + rcong = param['rcong'] + # Position des alésages de la bride permettant lintroduction des goujons (BRID_P_ALESAG) + dtrou = param['dtrou'] + # Profondeur (épaisseur) des filets des goujons (GOUJ_E_FILET) + pf = param['pf'] + # Hauteur de conduite (TUBU_H) + hcg1 = param['hcg1'] + # Diamètre intérieur de la bride (BRID_D_INT) + dint = param['dint'] + # Diamètre extérieur de la bride (BRID_D_EXT) + dext = param['dext'] + # Diamètre des alésages de la bride permettant lintroduction des goujons (BRID_D_ALESAG) + dt = param['dt'] + # Diamètre de lécrou (GOUJ_D_ECROU) + dec = param['dec'] + # Diamétre extérieur de la conduite (TUBU_D_EXT) + dex1 = param['dex1'] + # Diamètre de la rondelle (GOUJ_D_RONDEL) + drd = param['drd'] + # Diamètre de lépaulement de la bride au niveau de linterface détanchéité (BRID_D_EPAUL) + dex3 = param['dex3'] + # Position (diamètre) du congé de la bride (BRID_D_CONGE) + dex2 = param['dex2'] + + + + import math + import time + + # imports spécifiques à salome: + + import geompy + import smesh + import salome + #Import de l'interface graphique de GEOM (nécessaire pour affecter de la couleur) + gg = salome.ImportComponentGUI("GEOM") + + time_init = time.time() + + #************************** + # test sur les parametres * + #************************** + + #opti trac psc + + nrad = abs(nrad) + ncir = abs(ncir) + nver = abs(nver) + temps = abs(temps) + nsect = abs(nsect) + nbgouj = abs(nbgouj) + dint = abs(dint) + dex1 = abs(dex1) + dex2 = abs(dex2) + dex3 = abs(dex3) + dtrou = abs(dtrou) + dext = abs(dext) + dt = abs(dt) + drd = abs(drd) + dg = abs(dg) + dec = abs(dec) + rcong = abs(rcong) + he = abs(he) + e = abs(e) + j = abs(j) + hc1 = abs(hc1) + hcg1 = abs(hcg1) + hb = abs(hb) + htrou = abs(htrou) + pf = abs(pf) + j = abs(j) + + if (nbgouj <= 2): + nbgouj = 2 + + if ((nbgouj == 2) and (ncir <= 4)): + ncir = 4 + + if (dex1 <= dint): + dex1 = dint+10. + + if (dex2 <= dex1): + dex2 = dex1+20 + + if (dg >= dt): + dg = dt-2. + + if (dec <= dt): + if (drd > dt): + dec = (dt+drd)/2 + else: + dec = dt+2 + drd = dec+2 + + if ((dtrou-dt) <= dex2): + if (dtrou <= dex2): + dtrou = dex2+(2*dt) + else: + dt = (dtrou-dex2)/2 + + if (dex3 < dint): + dex3 = dtrou-drd + if (dex3 > (dtrou-dt)): + dex3 = dtrou-drd + + if ((dtrou+drd) <= (2*dtrou-dex3)): + a = 2*dtrou-dex3 + else: + a = dtrou+drd + if (dext < a): + dext = a + + if (hb >= hc1): + hb = 1 + if (htrou >= hc1): + htrou = hc1/2 + if ((hb+htrou) >= hc1): + hb = (hc1-htrou)/10 + + if (rcong >= 15.): + rcong = 15. + + rg = dg / 2. + if (pf >= rg): + pf = rg/2. + + if (nrad <= 0): + nrad = 1 + + if (ncir <= 1): + ncir = 2 + + if (nver <= 0): + nver = 1 + + if (nsect <= 3): + nsect = 4 + + if (temps <= 0): + nsect = 3 + + #************************** + # parametres intrinsèques * + #************************** + #opti dime 3 elem cub8 + #dens 1 + + # critere pour elim + crit = 0.0001 + + # rayons + rint = dint / 2. + rex1 = dex1 / 2. + rex2 = dex2 / 2. + rex3 = dex3 / 2. + rtrou = dtrou / 2. + rext = dext / 2. + rt = dt / 2. + rrd = drd / 2. + rg = dg / 2. + rec = dec / 2. + + # angle de coupe + stet = rrd/rtrou + ctet = stet**2 + ctet = 1-ctet + ctet = ctet**0.5 + tet = stet/ctet + tet = math.atan(tet) + beta = math.pi / nbgouj + # if (tet >= beta): + # tet = 1.1*tet + # beta = tet + # else: + # if ((beta-tet) < (0.1*tet)): + # tet = beta + # else: + # tet = (tet+beta)/2 + + # hauteurs + ht = hc1 + hcg1; + h = (hb + htrou); + mj = 0-j; + + p0 = geompy.MakeVertex(0., 0., 0.) + + Vx = geompy.MakeVectorDXDYDZ(1., 0., 0.) + Vy = geompy.MakeVectorDXDYDZ(0., 1., 0.) + Vz = geompy.MakeVectorDXDYDZ(0., 0., 1.) + + + p1 = geompy.MakeVertex(rint, 0., 0.) + p2 = geompy.MakeVertex(rint, 0., -j/2.) + p3 = geompy.MakeVertex(rex3, 0., -j/2.) + p4 = geompy.MakeVertex(rex3, 0., 0.) + + edge1 = geompy.MakeEdge(p1,p2) + edge2 = geompy.MakeEdge(p2,p3) + edge3 = geompy.MakeEdge(p3,p4) + edge4 = geompy.MakeEdge(p4,p1) + + wire_joint = geompy.MakeWire([edge1, edge2, edge3, edge4]) + face_joint = geompy.MakeFace(wire_joint,1) + + p5 = geompy.MakeVertex(rex3, 0., hb) + + edge5 = geompy.MakeEdge(p4, p5) + + p7 = geompy.MakeVertex(rext, 0., hb) + p8 = geompy.MakeVertex(rext, 0., hb + htrou) + p9 = geompy.MakeVertex(rex2, 0., hb + htrou) + + edge6 = geompy.MakeEdge(p5, p7) + edge9 = geompy.MakeEdge(p7, p8) + # construction du congé + # Rq: voir ce que l'on fait si la rondelle s'appuye sur le congé + + an = math.atan ((hc1-h)/(rex2-rex1)) + alpha = (math.pi - an)/2. + + conge_p1 = geompy.MakeVertex(rex2 + rcong/math.tan(alpha), 0., h) + conge_centre = geompy.MakeVertex(rex2 + rcong/math.tan(alpha), 0., h + rcong) + + gamma = math.pi - 2*alpha + + conge_axe = geompy.MakePrismVecH(conge_centre, Vy, 1) + conge = geompy.MakeRevolution(conge_p1, conge_axe, gamma) + # geompy.addToStudy(conge, "conge") + + conge_p2 = geompy.MakeRotation(conge_p1, conge_axe, gamma) + + ## fin construction du congé + + #avant le congé + edge10 = geompy.MakeEdge(p8, conge_p1) + # après le congé + p12 = geompy.MakeVertex(rex1, 0., hc1) + edge11 = geompy.MakeEdge(conge_p2, p12) + + p11 = geompy.MakeVertex(rint, 0., hc1) + p13 = geompy.MakeVertex(rex1, 0., hc1 + hcg1) + p14 = geompy.MakeVertex(rint, 0., hc1 + hcg1) + + edge16 = geompy.MakeEdge(p11, p14) + edge17 = geompy.MakeEdge(p14, p13) + edge18 = geompy.MakeEdge(p13, p12) + + edge19 = geompy.MakeEdge(p11, p1) + + wire_front = geompy.MakeWire([edge4, edge5, edge6, edge9, edge10, conge, + edge11, edge18, edge17, edge16, edge19]) + + face_front = geompy.MakeFace(wire_front,1) + # geompy.addToStudy(face_front, "face_front") + # on ajoutera le joint après la révolution + + + # Création du bloc solide + + bride_revol = geompy.MakeRevolution(face_front, Vz, beta) + # geompy.addToStudy(bride_revol, "bride_revol") + + # outil de coupe + p_goujon = geompy.MakeVertex(rtrou, 0., -j/2.) + cut_tool = geompy.MakeCylinder(p_goujon, Vz, rt, hc1 + j) + # geompy.addToStudy(cut_tool, "cut_tool") + + # coupe + bride_cut = geompy.MakeCut(bride_revol, cut_tool) + # geompy.addToStudy(bride_cut, "bride_cut") + + + + # goujon + axe_goujon = geompy.MakePrismVecH(p_goujon, Vz, 1.) + p_rayon_goujon = geompy.MakeTranslation(p_goujon, rg, 0., 0.) + rayon_goujon = geompy.MakeEdge(p_goujon, p_rayon_goujon) + # base_goujon = geompy.MakeRevolution(rayon_goujon, axe_goujon, math.pi) + h_goujon = j/2. + h + e + 1.5 * he + p1_filet_goujon = geompy.MakeVertex(rtrou + rg, 0., h_goujon - pf) + p2_filet_goujon = geompy.MakeVertex(rtrou + rg - pf, 0., h_goujon) + p_goujon_h = geompy.MakeVertex(rtrou, 0., h_goujon) + edge1_goujon = geompy.MakeEdge(p_rayon_goujon, p1_filet_goujon) + edge2_goujon = geompy.MakeEdge(p1_filet_goujon, p2_filet_goujon) + edge3_goujon = geompy.MakeEdge(p2_filet_goujon, p_goujon_h) + wire_goujon = geompy.MakeWire([rayon_goujon, edge1_goujon, edge2_goujon, edge3_goujon]) + face_ext_goujon = geompy.MakeRevolution(wire_goujon, axe_goujon, math.pi) + # geompy.addToStudy(face_ext_goujon, "face_ext_goujon") + wire_gauche_goujon = geompy.MakeRotation(wire_goujon, axe_goujon, math.pi) + + p2_filet_goujon_g = geompy.MakeRotation(p2_filet_goujon, axe_goujon, math.pi) + p_rayon_goujon_g = geompy.MakeTranslation(p_goujon, - rg, 0., 0.) + diametre_haut = geompy.MakeEdge(p2_filet_goujon_g, p2_filet_goujon) + diametre_bas = geompy.MakeEdge(p_rayon_goujon_g, p_rayon_goujon) + edge1_goujon_g = geompy.MakeRotation(edge1_goujon, axe_goujon, math.pi) + edge2_goujon_g = geompy.MakeRotation(edge2_goujon, axe_goujon, math.pi) + wire_face_int_goujon2 = geompy.MakeWire([diametre_bas, edge1_goujon, edge2_goujon, + diametre_haut, edge1_goujon_g, edge2_goujon_g]) + # geompy.addToStudy(wire_face_int_goujon2, "wire_face_int_goujon2") + face_int_goujon = geompy.MakeFace(wire_face_int_goujon2, 1) + + shell_goujon = geompy.MakeShell([face_ext_goujon, face_int_goujon]) + goujon = geompy.MakeSolid([shell_goujon]) + # geompy.addToStudy(goujon, "goujon") + + # ecrou + p1_ecrou = geompy.MakeVertex(rtrou + rg, 0., hb + htrou + e) + p2_ecrou = geompy.MakeVertex(rtrou + rec, 0., hb + htrou + e) + rayon_ecrou = geompy.MakeEdge(p1_ecrou, p2_ecrou) + base_ecrou = geompy.MakeRevolution(rayon_ecrou, axe_goujon, math.pi) + ecrou = geompy.MakePrismVecH(base_ecrou, Vz, he) + # geompy.addToStudy(ecrou, "ecrou") + + # assemblage goujon, ecrou: + goujon_ecrou = geompy.MakeCompound([goujon, ecrou]) + # geompy.addToStudy(goujon_ecrou, "goujon_ecrou") + + + # rondelle + p1_rondelle = geompy.MakeVertex(rtrou + rt, 0., hb + htrou) + p2_rondelle = geompy.MakeVertex(rtrou + rrd, 0., hb + htrou) + rayon_rondelle = geompy.MakeEdge(p1_rondelle, p2_rondelle) + # base_rondelle = geompy.MakeRevolution(rayon_rondelle, axe_goujon, math.pi) + face_rondelle = geompy.MakePrismVecH(rayon_rondelle, Vz, e) + rondelle = geompy.MakeRevolution(face_rondelle, axe_goujon, math.pi) + # geompy.addToStudy(rondelle, "rondelle") + + + # assemblage goujon-ecrou avec rondelle + goujon_ecrou_rondelle = geompy.MakeFuse(goujon_ecrou, rondelle) + # geompy.addToStudy(goujon_ecrou_rondelle, "goujon_ecrou_rondelle") + + + # assemblage sur la piece principale + bride_but_joint_tmp = geompy.MakeFuse(bride_cut, goujon_ecrou_rondelle) + # geompy.addToStudy(bride_but_joint_tmp, "bride_but_joint_tmp") + + # On partitionne avec la rondelle pour récupérer les solides physiques + bride_but_joint = geompy.MakePartition([bride_but_joint_tmp], [rondelle]) + # geompy.addToStudy(bride_but_joint, "bride_but_joint") + + # on ajoute le joint (compound pour que les faces communes soient en double) + joint = geompy.MakeRevolution(face_joint, Vz, beta) + # geompy.addToStudy(joint, "joint") + + bride_tmp = geompy.MakeCompound([bride_but_joint, joint]) + # geompy.addToStudy(bride_tmp, "bride_tmp") + + time_0 = time.time() + print "Temps Geometrie = ", (time_0-time_init) + + + # on partitionne pour obtenir certains points (P_BRI, P_GOU) + p6 = geompy.MakeVertex(rint, 0., hb) + edge_part = geompy.MakeEdge(p6, p7) + plan_part = geompy.MakeRevolution(edge_part, Vz, beta) + bride_part = geompy.MakePartition([bride_tmp], [plan_part]) + # geompy.addToStudy(bride_part, "bride_part") + + edge_p5_p6 = geompy.GetEdge(bride_part, p5, p6) + edge_p5_p6_ind = geompy.GetSubShapeID(bride_part, edge_p5_p6) + bride_vertex1 = geompy.DivideEdge(bride_part, edge_p5_p6_ind, 1./3., 1) + # geompy.addToStudy(bride_vertex1, "bride_vertex1") + + p_gouj_part1 = geompy.MakeVertex(rtrou - rg, 0., hb) + p_gouj_part2 = geompy.MakeVertex(rtrou + rg, 0., hb) + p_gouj_mid = geompy.MakeVertex(rtrou, 0., hb) + edge_gouj_part = geompy.GetEdgeNearPoint(bride_vertex1, p_gouj_mid) + edge_gouj_part_ind = geompy.GetSubShapeID(bride_vertex1, edge_gouj_part) + bride = geompy.DivideEdge(bride_vertex1, edge_gouj_part_ind, 0.5, 1) +# geompy.addToStudy(bride, "bride") + geompy.addToStudy(bride, name ) + + + time_1 = time.time() + print "Temps Partition = ", (time_1-time_0) + + # Détermination des différents solides et affichage en couleur + + idToDisplay=[] + + GOUJON1 = geompy.GetBlockNearPoint(bride, p_goujon) + GOUJON2 = geompy.GetBlockNearPoint(bride, p_goujon_h) + list_GOUJON = [GOUJON1, GOUJON2] + + GOUJON = geompy.CreateGroup(bride, geompy.ShapeType["SOLID"]) + id_GOUJON = geompy.addToStudyInFather(bride, GOUJON, "GOUJON") + idToDisplay.append(id_GOUJON) + + for solid in list_GOUJON: + f_ind_tmp = geompy.GetSubShapeID(bride, solid) + geompy.AddObject(GOUJON, f_ind_tmp) + + + JOINT = geompy.GetBlockNearPoint(bride, p2) + id_JOINT = geompy.addToStudyInFather(bride, JOINT, "JOINT") + idToDisplay.append(id_JOINT) + + p_rondelle = geompy.MakeTranslation(p2_rondelle, 0., 0., e) + ROND = geompy.GetBlockNearPoint(bride, p_rondelle) + id_ROND = geompy.addToStudyInFather(bride, ROND, "ROND") + idToDisplay.append(id_ROND) + + p_ecrou = geompy.MakeTranslation(p2_ecrou, 0., 0., he) + ECROU = geompy.GetBlockNearPoint(bride, p_ecrou) + id_ECROU = geompy.addToStudyInFather(bride, ECROU, "ECROU") + idToDisplay.append(id_ECROU) + + list_BRIDE = [] + + BRIDE1 = geompy.GetBlockNearPoint(bride, p14) + list_BRIDE.append(BRIDE1) + + p_BRIDE2 = geompy.MakeVertex((rint+rex3)/2., 0., hb/2.) + p_BRIDE2_rota = geompy.MakeRotation(p_BRIDE2, Vz, beta/2.) + BRIDE2 = geompy.GetBlockNearPoint(bride, p_BRIDE2_rota) + list_BRIDE.append(BRIDE2) + + BRIDE = geompy.CreateGroup(bride, geompy.ShapeType["SOLID"]) + id_BRIDE = geompy.addToStudyInFather(bride, BRIDE, "BRIDE") + idToDisplay.append(id_BRIDE) + + for solid in list_BRIDE: + f_ind_tmp = geompy.GetSubShapeID(bride, solid) + geompy.AddObject(BRIDE, f_ind_tmp) + + RedGreenBlue = [[189,97,0],[255,215,0],[255,0,0],[0,176,0],[0,0,255]] + + for i in range(len(idToDisplay)): + gg.createAndDisplayGO(idToDisplay[i]) + gg.setDisplayMode(idToDisplay[i],1) + gg.setColor(idToDisplay[i],RedGreenBlue[i][0],RedGreenBlue[i][1],RedGreenBlue[i][2]) + + list_VTOT = geompy.SubShapeAllSorted(bride, geompy.ShapeType["SOLID"]) + VTOT = geompy.CreateGroup(bride, geompy.ShapeType["SOLID"]) + for solid in list_VTOT: + f_ind_tmp = geompy.GetSubShapeID(bride, solid) + geompy.AddObject(VTOT, f_ind_tmp) + geompy.addToStudyInFather(bride, VTOT, "VTOT") + + time_2 = time.time() + print "Temps affichage solides = ", (time_2-time_1) + + # Détermination des différentes faces + + p_SCEG = geompy.MakeVertex(rtrou, rg, h + e + he/2.) + SCEG = geompy.GetFaceNearPoint(ECROU, p_SCEG) + geompy.addToStudyInFather(bride, SCEG, "SCEG") + + SCGE = geompy.GetFaceNearPoint(GOUJON, p_SCEG) + geompy.addToStudyInFather(bride, SCGE, "SCGE") + + p_2_3 = geompy.MakeVertexOnCurve(edge2, 0.5) + p_2_3_rota = geompy.MakeRotation(p_2_3, Vz, beta/2.) + M_JOI = geompy.GetFaceNearPoint(JOINT, p_2_3_rota) + geompy.addToStudyInFather(bride, M_JOI, "M_JOI") + + p_1_4 = geompy.MakeVertexOnCurve(edge4, 0.5) + p_1_4_rota = geompy.MakeRotation(p_1_4, Vz, beta/2.) + SCJB = geompy.GetFaceNearPoint(JOINT, p_1_4_rota) + geompy.addToStudyInFather(bride, SCJB, "SCJB") + + SCBJ = geompy.GetFaceNearPoint(BRIDE, p_1_4_rota) + geompy.addToStudyInFather(bride, SCBJ, "SCBJ") + + p1_rota = geompy.MakeRotation(p1, Vz, beta/2.) + Vint = geompy.MakeVector(p1_rota, p0) + list_M_INT = geompy.GetShapesOnPlane(bride, geompy.ShapeType["FACE"], + Vint, geompy.GEOM.ST_ONOUT) + + M_INT = geompy.CreateGroup(bride, geompy.ShapeType["FACE"]) + geompy.addToStudyInFather(bride, M_INT, "M_INT") + + for face in list_M_INT: + f_ind_tmp = geompy.GetSubShapeID(bride, face) + geompy.AddObject(M_INT, f_ind_tmp) + + + p_13_14 = geompy.MakeVertexOnCurve(edge17, 0.5) + p_13_14_rota = geompy.MakeRotation(p_13_14, Vz, beta/2.) + M_TUB = geompy.GetFaceNearPoint(BRIDE, p_13_14_rota) + geompy.addToStudyInFather(bride, M_TUB, "M_TUB") + + p_M_GOU = geompy.MakeVertex(rtrou, rg/2., -j/2.) + M_GOU = geompy.GetFaceNearPoint(GOUJON, p_M_GOU) + geompy.addToStudyInFather(bride, M_GOU, "M_GOU") + + Vy_rota = geompy.MakeRotation(Vy, Vz, beta) + list_M_L_SA = geompy.GetShapesOnPlane(bride, geompy.ShapeType["FACE"], + Vy_rota, geompy.GEOM.ST_ONOUT) + + M_L_SA = geompy.CreateGroup(bride, geompy.ShapeType["FACE"]) + geompy.addToStudyInFather(bride, M_L_SA, "M_L_SA") + + for face in list_M_L_SA: + f_ind_tmp = geompy.GetSubShapeID(bride, face) + geompy.AddObject(M_L_SA, f_ind_tmp) + + + moins_Vy = geompy.ChangeOrientation(Vy) + list_M_L_AA = geompy.GetShapesOnPlane(bride, geompy.ShapeType["FACE"], + moins_Vy, geompy.GEOM.ST_ONOUT) + + M_L_AA = geompy.CreateGroup(bride, geompy.ShapeType["FACE"]) + geompy.addToStudyInFather(bride, M_L_AA, "M_L_AA") + + for face in list_M_L_AA: + f_ind_tmp = geompy.GetSubShapeID(bride, face) + geompy.AddObject(M_L_AA, f_ind_tmp) + + + list_M_EXT = [] + + p_4_5 = geompy.MakeVertexOnCurve(edge5, 0.5) + p_4_5_rota = geompy.MakeRotation(p_4_5, Vz, beta/2.) + M_EXT1 = geompy.GetFaceNearPoint(bride, p_4_5_rota) + list_M_EXT.append(M_EXT1) + + p_M_EXT2 = geompy.MakeVertex((rex3 + rtrou-rt)/2., 0., hb) + p_M_EXT2_rota = geompy.MakeRotation(p_M_EXT2, Vz, beta/2.) + M_EXT2 = geompy.GetFaceNearPoint(bride, p_M_EXT2_rota) + list_M_EXT.append(M_EXT2) + + p_7_8 = geompy.MakeVertexOnCurve(edge9, 0.5) + p_7_8_rota = geompy.MakeRotation(p_7_8, Vz, beta/2.) + M_EXT3 = geompy.GetFaceNearPoint(bride, p_7_8_rota) + list_M_EXT.append(M_EXT3) + + p_M_EXT4 = geompy.MakeVertex((rext + rtrou+rrd)/2., 0., h) + p_M_EXT4_rota = geompy.MakeRotation(p_M_EXT4, Vz, beta/2.) + M_EXT4 = geompy.GetFaceNearPoint(bride, p_M_EXT4_rota) + list_M_EXT.append(M_EXT4) + + p_mid_conge = geompy.MakeVertexOnCurve(conge, 0.5) + p_mid_conge_rota = geompy.MakeRotation(p_mid_conge, Vz, beta/2.) + M_EXT5 = geompy.GetFaceNearPoint(bride, p_mid_conge_rota) + list_M_EXT.append(M_EXT5) + + p_9_12 = geompy.MakeVertexOnCurve(edge11, 0.5) + p_9_12_rota = geompy.MakeRotation(p_9_12, Vz, beta/2.) + M_EXT6 = geompy.GetFaceNearPoint(bride, p_9_12_rota) + list_M_EXT.append(M_EXT6) + + p_12_13 = geompy.MakeVertexOnCurve(edge18, 0.5) + p_12_13_rota = geompy.MakeRotation(p_12_13, Vz, beta/2.) + M_EXT7 = geompy.GetFaceNearPoint(bride, p_12_13_rota) + list_M_EXT.append(M_EXT7) + + p_haut_gouj = geompy.MakeTranslation(p_goujon_h, 0., (rg-pf)/2., 0.) + M_EXT8 = geompy.GetFaceNearPoint(bride, p_haut_gouj) + list_M_EXT.append(M_EXT8) + + p_filet_gouj = geompy.MakeVertexOnCurve(edge2_goujon, 0.5) + p_filet_gouj_rota = geompy.MakeRotation(p_filet_gouj, axe_goujon, math.pi/2.) + M_EXT9 = geompy.GetFaceNearPoint(bride, p_filet_gouj_rota) + list_M_EXT.append(M_EXT9) + + p_cote_gouj = geompy.MakeVertex(rtrou, rg, (h+e+he + h_goujon) / 2.) + M_EXT10 = geompy.GetFaceNearPoint(bride, p_cote_gouj) + list_M_EXT.append(M_EXT10) + + p_haut_ecrou = geompy.MakeVertex(rtrou, (rg+rec)/2., h + e + he) + M_EXT11 = geompy.GetFaceNearPoint(bride, p_haut_ecrou) + list_M_EXT.append(M_EXT11) + + p_cote_ecrou = geompy.MakeVertex(rtrou, rec, h + e + he/2.) + M_EXT12 = geompy.GetFaceNearPoint(bride, p_cote_ecrou) + list_M_EXT.append(M_EXT12) + + p_haut_rondelle = geompy.MakeVertex(rtrou, (rec+rrd)/2., h + e) + M_EXT13 = geompy.GetFaceNearPoint(bride, p_haut_rondelle) + list_M_EXT.append(M_EXT13) + + p_cote_rondelle = geompy.MakeVertex(rtrou, rrd, h + e/2.) + M_EXT14 = geompy.GetFaceNearPoint(bride, p_cote_rondelle) + list_M_EXT.append(M_EXT14) + + M_EXT = geompy.CreateGroup(bride, geompy.ShapeType["FACE"]) + geompy.addToStudyInFather(bride, M_EXT, "M_EXT") + + for face in list_M_EXT: + f_ind_tmp = geompy.GetSubShapeID(bride, face) + geompy.AddObject(M_EXT, f_ind_tmp) + + time_3 = time.time() + print "Temps explosion des faces = ", (time_3-time_2) + + + # Détermination des différents points + + epsilon = 1e-7 + + #P4 + PJE_OUV = geompy.GetPoint(JOINT, rex3, 0., 0., epsilon) + geompy.addToStudyInFather(bride, PJE_OUV, "PJE_OUV") + + #P1 + PJI_OUV = geompy.GetPoint(JOINT, rint, 0., 0., epsilon) + geompy.addToStudyInFather(bride, PJI_OUV, "PJI_OUV") + + #P4 + PBE_OUV = geompy.GetPoint(BRIDE, rex3, 0., 0., epsilon) + geompy.addToStudyInFather(bride, PBE_OUV, "PBE_OUV") + + #P1 + PBI_OUV = geompy.GetPoint(BRIDE, rint, 0., 0., epsilon) + geompy.addToStudyInFather(bride, PBI_OUV, "PBI_OUV") + + P_ECR = geompy.GetPoint(bride, rtrou + rt, 0., h + e, epsilon) + geompy.addToStudyInFather(bride, P_ECR, "P_ECR") + + P_GOU = geompy.GetPoint(bride, rtrou, 0., hb, epsilon) + geompy.addToStudyInFather(bride, P_GOU, "P_GOU") + + P_BRI = geompy.GetPoint(bride, rint + (rex3-rint)/3., 0., hb, epsilon) + geompy.addToStudyInFather(bride, P_BRI, "P_BRI") + + time_4 = time.time() + print "Temps explosion des vertices = ", (time_4-time_3) + + +# #=============== MAILLAGE ====================== +# +# AverageLength = (rex1-rint)/nrad +# +# # Creation du maillage +# # -------------------- +# +# maillageBride = smesh.Mesh(bride, "MeshBride") +# +# # Algorithmes et hypotheses globales +# # ---------------------------------- +# +# # 1D +# +# algo = maillageBride.Segment() +# algo.LocalLength(AverageLength) +# # On veut un maillage quadratique +# algo.QuadraticMesh() +# +# # 2D +# +# algo = maillageBride.Triangle() +# algo.LengthFromEdges() +# +# # 3D +# +# maillageBride.Tetrahedron(smesh.NETGEN) +# +# +# # Calcul +# # ------ +# +# maillageBride.Compute() +# time_5 = time.time() +# print "Temps Maillage = ", (time_5-time_4) +# +# # Création des groupes +# # -------------------- +# +# maillageBride.Group(P_GOU, "P_GOU") +# maillageBride.Group(PBI_OUV, "PBI_OUV") +# maillageBride.Group(P_BRI, "P_BRI") +# maillageBride.Group(PBE_OUV, "PBE_OUV") +# maillageBride.Group(P_ECR, "P_ECR") +# maillageBride.Group(PJI_OUV, "PJI_OUV") +# maillageBride.Group(SCEG, "SCEG") +# maillageBride.Group(SCGE, "SCGE") +# maillageBride.Group(M_JOI, "M_JOI") +# maillageBride.Group(SCJB, "SCJB") +# maillageBride.Group(SCBJ, "SCBJ") +# maillageBride.Group(M_EXT, "M_EXT") +# maillageBride.Group(M_INT, "M_INT") +# maillageBride.Group(M_TUB, "M_TUB") +# maillageBride.Group(M_GOU, "M_GOU") +# maillageBride.Group(M_L_SA, "M_L_SA") +# maillageBride.Group(M_L_AA, "M_L_AA") +# maillageBride.Group(GOUJON, "GOUJON") +# maillageBride.Group(ROND, "ROND") +# maillageBride.Group(ECROU, "ECROU") +# maillageBride.Group(BRIDE, "BRIDE") +# maillageBride.Group(JOINT, "JOINT") +# maillageBride.Group(VTOT, "VTOT") +# +# time_6 = time.time() +# +# print "Temps Groupes Maillage = ", (time_6-time_5) +# # Mise à jour de l'arbre d'étude + + salome.sg.updateObjBrowser(1) + + + + \ No newline at end of file