doc/files/tutorial_4.py

   1 #!/usr/bin/env python
   2 # -*- coding: iso-8859-1 -*-
   3
   4 # Copyright (C) 2011-2014  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-R&D 1996, 2011, 2014
  26 """
  27 __revision__ = "V2.4"
  28 #
  29 import os
  30 import sys
  31 #
  32 # ==================================
  33 # Repertoire a personnaliser
  34 # Ce repertoire contiendra les fichiers de resultats : maill.01.med, maill.02.med, maill.03.med
  35 if os.environ.has_key("LOGNAME") :
  36   user = os.environ ["LOGNAME"]
  37 else :
  38   user = "anonymous"
  39 dircase = os.path.join( os.sep, "tmp", "HOMARD_"+user)
  40 if not os.path.isdir(dircase) :
  41     os.mkdir (dircase)
  42 dircase = os.path.join( dircase, "tutorial_4" )
  43 if not os.path.isdir(dircase) :
  44   os.mkdir (dircase)
  45 # ==================================
  46 # Ce repertoire contient les fichiers de donnees : tutorial_4.00.med, tutorial_4.fr.med
  47 pathHomard = os.getenv('HOMARD_ROOT_DIR')
  48 data_dir = os.path.join(pathHomard, "share", "doc", "salome", "gui", "HOMARD", "fr", "_downloads")
  49 sys.path.append(data_dir)
  50 from tutorial_util import gzip_gunzip
  51 gzip_gunzip(data_dir, 4, -1)
  52 # ==================================
  53 #
  54 import salome
  55 salome.salome_init()
  56 import HOMARD
  57 #
  58 homard = salome.lcc.FindOrLoadComponent("FactoryServer", "HOMARD")
  59 study_main = salome.myStudyManager.NewStudy("HOMARD")
  60 homard.SetCurrentStudy(salome.myStudy)
  61 #
  62 # Frontieres
  63 # ==========
  64 Boun_4_1 = homard.CreateBoundaryDi('intersection', 'PIQUAGE', data_dir+'/tutorial_4.fr.med')
  65 #
  66 Boun_4_2 = homard.CreateBoundaryCylinder('cyl_1_ext', 0.0, 25., -25., 25., 50., 75., 100.)
  67 #
  68 Boun_4_3 = homard.CreateBoundaryCylinder('cyl_2_ext', 17.5, -2.5, -12.5, -100., -75., -25., 50.)
  69 #
  70 Boun_4_4 = homard.CreateBoundaryCylinder('cyl_1_int', 0.0, 25., -25., 25., 50., 75., 75.)
  71 #
  72 Boun_4_5 = homard.CreateBoundaryCylinder('cyl_2_int', 17.5, -2.5, -12.5, -100., -75., -25., 25.)
  73 #
  74 # Hypotheses
  75 # ==========
  76 # Creation of the hypothesis Hypo_4
  77 Hypo_4 = homard.CreateHypothesis('Hypo_4')
  78 Hypo_4.SetUnifRefinUnRef(1)
  79 Hypo_4.AddGroup('T1_INT_I')
  80 Hypo_4.AddGroup('T1_INT_O')
  81 Hypo_4.AddGroup('T2_INT')
  82 # Creation of the hypothesis Hypo_4_bis
  83 Hypo_4_bis = homard.CreateHypothesis('Hypo_4_bis')
  84 Hypo_4_bis.SetUnifRefinUnRef(1)
  85 Hypo_4_bis.AddGroup('T1_EXT_I')
  86 Hypo_4_bis.AddGroup('T1_EXT_O')
  87 Hypo_4_bis.AddGroup('T2_EXT')
  88 #
  89 # Cas
  90 # ===
  91 Case_4 = homard.CreateCase('Case_4', 'PIQUAGE', data_dir+'/tutorial_4.00.med')
  92 Case_4.SetDirName(dircase)
  93 Case_4.AddBoundaryGroup( 'intersection', '' )
  94 Case_4.AddBoundaryGroup( 'cyl_1_int', 'T1_INT_I' )
  95 Case_4.AddBoundaryGroup( 'cyl_1_ext', 'T1_EXT_I' )
  96 Case_4.AddBoundaryGroup( 'cyl_1_int', 'T1_INT_O' )
  97 Case_4.AddBoundaryGroup( 'cyl_1_ext', 'T1_EXT_O' )
  98 Case_4.AddBoundaryGroup( 'cyl_2_int', 'T2_INT' )
  99 Case_4.AddBoundaryGroup( 'cyl_2_ext', 'T2_EXT' )
 100 #
 101 # Iterations
 102 # ==========
 103 # Iteration Iter_4_1 : raffinement selon les faces internes
 104 Iter_4_1 = Case_4.NextIteration('Iter_4_1')
 105 Iter_4_1.SetMeshName('PIQUAGE_1')
 106 Iter_4_1.SetMeshFile(dircase+'/maill.01.med')
 107 Iter_4_1.AssociateHypo('Hypo_4')
 108 codret = Iter_4_1.Compute(1, 2)
 109 # Iteration Iter_4_2 : raffinement selon les faces externes
 110 Iter_4_2 = Iter_4_1.NextIteration('Iter_4_2')
 111 Iter_4_2.SetMeshName('PIQUAGE_2')
 112 Iter_4_2.SetMeshFile(dircase+'/maill.02.med')
 113 Iter_4_2.AssociateHypo('Hypo_4_bis')
 114 codret = Iter_4_2.Compute(1, 2)
 115 # Iteration Iter_4_3 : second raffinement selon les faces externes
 116 Iter_4_3 = Iter_4_2.NextIteration('Iter_4_3')
 117 Iter_4_3.SetMeshName('PIQUAGE_3')
 118 Iter_4_3.SetMeshFile(dircase+'/maill.03.med')
 119 Iter_4_3.AssociateHypo('Hypo_4_bis')
 120 codret = Iter_4_3.Compute(1, 2)
 121
 122 # ==================================
 123 gzip_gunzip(data_dir, 4, 1)
 124 # ==================================
 125
 126 if salome.sg.hasDesktop():
 127   salome.sg.updateObjBrowser(1)