doc/files/tutorial_3.py

   1 #!/usr/bin/env python
   2 # -*- coding: iso-8859-1 -*-
   3
   4 # Copyright (C) 2011-2012  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.
  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, 2010
  26 """
  27 __revision__ = "V1.3"
  28 #
  29 # ==================================
  30 # Repertoire a personnaliser
  31 # Ce repertoire contient les fichiers de donnees : tutorial_3.00.med, tutorial_3.01.med
  32 # Ce repertoire contiendra les fichiers de resultats : maill.01.med, maill.02.med
  33 dircase = "/tmp"
  34 # ==================================
  35 #
  36 import salome
  37 salome.salome_init()
  38 import HOMARD
  39 #
  40 homard = salome.lcc.FindOrLoadComponent("FactoryServer", "HOMARD")
  41 study_main = salome.myStudyManager.NewStudy("HOMARD")
  42 homard.SetCurrentStudy(salome.myStudy)
  43 #
  44 # Hypothesis "Hypo_0vers1"
  45 # ========================
  46 Hypo_0vers1 = homard.CreateHypothesis('Hypo_0vers1')
  47 Hypo_0vers1.SetAdapRefinUnRef(1, 1, 0)
  48 # Characterization of the field
  49 Hypo_0vers1.SetField('SOLU_0__QIRE_ELEM_SIGM__________')
  50 Hypo_0vers1.SetUseComp(0)
  51 Hypo_0vers1.AddComp('ERREST          ')
  52 Hypo_0vers1.SetRefinThr(3, 1.0)
  53 Hypo_0vers1.SetTypeFieldInterp(2)
  54 Hypo_0vers1.AddFieldInterp('SOLU_0__DEPL____________________')
  55 Hypo_0vers1.AddFieldInterp('SOLU_0__ERRE_ELEM_SIGM__________')
  56 #
  57 # Hypothesis "Hypo_1vers2"
  58 # ========================
  59 Hypo_1vers2 = homard.CreateHypothesis('Hypo_1vers2')
  60 Hypo_1vers2.SetAdapRefinUnRef(1, 1, 1)
  61 # Characterization of the field
  62 Hypo_1vers2.SetField('SOLU_1__QIRE_ELEM_SIGM__________')
  63 Hypo_1vers2.SetUseComp(0)
  64 Hypo_1vers2.AddComp('ERREST          ')
  65 Hypo_1vers2.SetRefinThr(3, 1.5)
  66 Hypo_1vers2.SetUnRefThr(3, 6.)
  67 Hypo_1vers2.SetTypeFieldInterp(2)
  68 Hypo_1vers2.AddFieldInterp('SOLU_1__DEPL____________________')
  69 Hypo_1vers2.AddFieldInterp('SOLU_1__QIRE_ELEM_SIGM__________')
  70 #
  71 # Hypothesis "Hypo_1vers2_bis"
  72 # ============================
  73 Hypo_1vers2_bis = homard.CreateHypothesis('Hypo_1vers2_bis')
  74 Hypo_1vers2_bis.SetAdapRefinUnRef(1, 1, 1)
  75 # Characterization of the field
  76 Hypo_1vers2_bis.SetField('SOLU_1__DEPL____________________')
  77 Hypo_1vers2_bis.SetUseComp(1)
  78 Hypo_1vers2_bis.AddComp('DX')
  79 Hypo_1vers2_bis.AddComp('DY')
  80 Hypo_1vers2_bis.AddComp('DZ')
  81 Hypo_1vers2_bis.SetRefinThr(1, 0.0001)
  82 Hypo_1vers2_bis.SetUnRefThr(1, 0.000001)
  83 Hypo_1vers2_bis.SetTypeFieldInterp(0)
  84 #
  85 # Case "Case_1"
  86 # =============
  87 Case_1 = homard.CreateCase('Case_1', 'G_0', dircase+'/tutorial_3.00.med')
  88 Case_1.SetDirName(dircase)
  89 #
  90 # Iteration "Iter_1"
  91 # ==================
  92 Iter_1 = homard.CreateIteration('Iter_1', Case_1.GetIter0Name())
  93 Iter_1.SetMeshName('H_1')
  94 Iter_1.SetMeshFile(dircase+'/maill.01.med')
  95 Iter_1.SetFieldFile(dircase+'/tutorial_3.00.med')
  96 Iter_1.SetTimeStepRank( 1, 1)
  97 homard.AssociateIterHypo('Iter_1', 'Hypo_0vers1')
  98 codret = Iter_1.Compute(1)
  99 #
 100 # Iteration "Iter_2"
 101 # ==================
 102 Iter_2 = homard.CreateIteration('Iter_2', 'Iter_1')
 103 Iter_2.SetMeshName('H_2')
 104 Iter_2.SetMeshFile(dircase+'/maill.02.med')
 105 Iter_2.SetFieldFile(dircase+'/tutorial_3.01.med')
 106 Iter_2.SetTimeStepRank(1, 1)
 107 homard.AssociateIterHypo('Iter_2', 'Hypo_1vers2')
 108 codret = Iter_2.Compute(1)
 109 #
 110 # Iteration "Iter_2_bis"
 111 # ======================
 112 Iter_2_bis = homard.CreateIteration('Iter_2_bis', 'Iter_1')
 113 Iter_2_bis.SetMeshName('H_2_bis')
 114 Iter_2_bis.SetMeshFile(dircase+'/maill.02.bis.med')
 115 Iter_2_bis.SetFieldFile(dircase+'/tutorial_3.01.med')
 116 Iter_2_bis.SetTimeStepRank(1, 1)
 117 homard.AssociateIterHypo('Iter_2_bis', 'Hypo_1vers2_bis')
 118 codret = Iter_2_bis.Compute(1)
 119 #
 120 if salome.sg.hasDesktop():
 121   salome.sg.updateObjBrowser(1)