doc/files/tutorial_6.py

   1 #!/usr/bin/env python3
   2 # -*- coding: utf-8 -*-
   3
   4 # Copyright (C) 2011-2023  CEA/DEN, EDF R&D
   5 #
   6 # This library is free software; you can redistribute it and/or
   7 # modify it under the terms of the GNU Lesser General Public
   8 # License as published by the Free Software Foundation; either
   9 # version 2.1 of the License, or (at your option) any later version.
  10 #
  11 # This library is distributed in the hope that it will be useful,
  12 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 # Lesser General Public License for more details.
  15 #
  16 # You should have received a copy of the GNU Lesser General Public
  17 # License along with this library; if not, write to the Free Software
  18 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 #
  20 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 #
  22
  23 """
  24 Exemple de couplage HOMARD-Salome
  25 Copyright EDF 1996, 2011, 2019
  26 """
  27 __revision__ = "V4.06"
  28 #
  29 import os
  30 import sys
  31 #
  32 # ==================================
  33 PATH_HOMARD = os.getenv("HOMARD_ROOT_DIR")
  34 # Repertoire des donnees du tutorial
  35 DATA_TUTORIAL = os.path.join(PATH_HOMARD, "share", "doc", "salome", "gui", "HOMARD", "fr", "_downloads")
  36 DATA_TUTORIAL = os.path.normpath(DATA_TUTORIAL)
  37 sys.path.append(DATA_TUTORIAL)
  38 from tutorial_util import creation_dircase
  39 # ==================================
  40 DIRCASE = creation_dircase(6)
  41 # ==================================
  42 #
  43 import salome
  44 salome.salome_init()
  45 import HOMARD
  46 #
  47 homard = salome.lcc.FindOrLoadComponent("FactoryServer", "HOMARD")
  48 homard.UpdateStudy()
  49 #
  50 #============================= Début des commandes =============================
  51 #
  52 # Frontières
  53 # ==========
  54 boun_6_1 = homard.CreateBoundaryDi('intersection', 'COURBES', os.path.join(DATA_TUTORIAL, "tutorial_6.fr.med"))
  55 #
  56 boun_6_2 = homard.CreateBoundaryCylinder('cyl_1_ext', 0.0, 25., -25., 25., 50., 75., 100.)
  57 #
  58 boun_6_3 = homard.CreateBoundaryCylinder('cyl_2_ext', 17.5, -2.5, -12.5, -100., -75., -25., 50.)
  59 #
  60 boun_6_4 = homard.CreateBoundaryCylinder('cyl_1_int', 0.0, 25., -25., 25., 50., 75., 75.)
  61 #
  62 boun_6_5 = homard.CreateBoundaryCylinder('cyl_2_int', 17.5, -2.5, -12.5, -100., -75., -25., 25.)
  63 #
  64 # Cas
  65 # ===
  66 le_cas = homard.CreateCase('Case_6', 'PIQUAGE', os.path.join(DATA_TUTORIAL, "tutorial_4.00.med"))
  67 le_cas.SetDirName(DIRCASE)
  68 le_cas.AddBoundary( 'intersection' )
  69 le_cas.AddBoundaryGroup( 'cyl_1_int', 'T1_INT_I' )
  70 le_cas.AddBoundaryGroup( 'cyl_1_ext', 'T1_EXT_I' )
  71 le_cas.AddBoundaryGroup( 'cyl_1_int', 'T1_INT_O' )
  72 le_cas.AddBoundaryGroup( 'cyl_1_ext', 'T1_EXT_O' )
  73 le_cas.AddBoundaryGroup( 'cyl_2_int', 'T2_INT' )
  74 le_cas.AddBoundaryGroup( 'cyl_2_ext', 'T2_EXT' )
  75 #
  76 #============================== Fin des commandes ==============================
  77 #
  78 # Hypotheses
  79 # ==========
  80 # Creation of the hypothesis hypo_6
  81 l_hypothese = homard.CreateHypothesis('hypo_6')
  82 l_hypothese.SetUnifRefinUnRef(1)
  83 l_hypothese.AddGroup('IN1')
  84 l_hypothese.AddGroup('IN2')
  85 l_hypothese.AddGroup('T1_INT_I')
  86 l_hypothese.AddGroup('T1_INT_O')
  87 l_hypothese.AddGroup('T2_INT')
  88 # Creation of the hypothesis hypo_6_bis
  89 l_hypothese_bis = homard.CreateHypothesis('hypo_6_bis')
  90 l_hypothese_bis.SetUnifRefinUnRef(1)
  91 l_hypothese_bis.AddGroup('T1_EXT_I')
  92 l_hypothese_bis.AddGroup('T1_EXT_O')
  93 l_hypothese_bis.AddGroup('T2_EXT')
  94 #
  95 # Iterations
  96 # ==========
  97 # Iteration iter_6_1 : raffinement selon les faces internes
  98 iter_6_1 = le_cas.NextIteration('iter_6_1')
  99 iter_6_1.SetMeshName('PIQUAGE_6_1')
 100 iter_6_1.SetMeshFile(os.path.join(DIRCASE, "maill.01.med"))
 101 iter_6_1.AssociateHypo('hypo_6')
 102 erreur = iter_6_1.Compute(1, 2)
 103 # Iteration iter_6_2 : raffinement selon les faces externes
 104 iter_6_2 = iter_6_1.NextIteration('iter_6_2')
 105 iter_6_2.SetMeshName('PIQUAGE_6_2')
 106 iter_6_2.SetMeshFile(os.path.join(DIRCASE, "maill.02.med"))
 107 iter_6_2.AssociateHypo('hypo_6_bis')
 108 erreur = iter_6_2.Compute(1, 2)
 109 # Iteration iter_6_3 : second raffinement selon les faces externes
 110 iter_6_3 = iter_6_2.NextIteration('iter_6_3')
 111 iter_6_3.SetMeshName('PIQUAGE_6_3')
 112 iter_6_3.SetMeshFile(os.path.join(DIRCASE, "maill.03.med"))
 113 iter_6_3.AssociateHypo('hypo_6_bis')
 114 erreur = iter_6_3.Compute(1, 2)
 115 #
 116 if salome.sg.hasDesktop():
 117   salome.sg.updateObjBrowser()