Homard/homard_cata_STA7.py

   1 # -*- coding: utf-8 -*-
   2 # --------------------------------------------------
   3 # debut entete
   4 # --------------------------------------------------
   5
   6 import Accas
   7 from Accas import *
   8
   9 #CONTEXT.debug=1
  10
  11 JdC = JDC_CATA(code='HOMARD',
  12                execmodul=None,
  13                regles = (AU_MOINS_UN('DONNEES_HOMARD'),),
  14                         )
  15
  16 # Type le plus general
  17 class entier  (ASSD):pass
  18 class reel    (ASSD):pass
  19 class complexe(ASSD):pass
  20 class liste   (ASSD):pass
  21 class chaine  (ASSD):pass
  22
  23 # Types geometriques
  24 class no  (GEOM):pass
  25 class grno(GEOM):pass
  26 class ma  (GEOM):pass
  27 class grma(GEOM):pass
  28
  29 # --------------------------------------------------
  30 # fin entete
  31 # --------------------------------------------------
  32
  33
  34 DONNEES_HOMARD=PROC(nom="DONNEES_HOMARD",op= 189, docu="U7.04.01-b",
  35             UIinfo={"groupes":("Fonction",)},
  36                     fr="Imprime le fichier de configuration de HOMARD.",
  37                     ang="Writes the configuration file for HOMARD.",
  38 #
  39 # 1. Langue des messages issus de HOMARD
  40 #
  41          LANGUE = SIMP(statut='f',typ='TXM',defaut="FRANCAIS",
  42                                into=("FRANCAIS","FRENCH","ANGLAIS","ENGLISH",),
  43                            fr="Langue des messages issus de HOMARD.",
  44                            ang="Language for HOMARD messages." ),
  45 #
  46 # 2. Les fichiers en entree/sortie
  47 #
  48          FICHIERS = FACT(statut='o',
  49                 NOM_MED_MAILLAGE_N = SIMP(statut='o' ,typ='TXM',
  50                          fr="Nom du fichier contenant le maillage à adapter",
  51                          ang="File name including the in-mesh"),
  52                 NOM_MED_MAILLAGE_NP1 = SIMP(statut='o' ,typ='TXM',
  53                          fr="Nom du fichier qui contiendra le maillage adapté",
  54                          ang="File name for resulting out-mesh"),
  55                 INFORMATION =  SIMP(statut='o' ,typ='TXM',
  56                          fr="Nom du fichier contenant les informations de sortie",
  57                          ang="Log File"),
  58                      ),
  59 #
  60 # 3. Le type de traitement :
  61 #
  62          TRAITEMENT      =FACT(statut='o',
  63 #
  64 # 3.1. DEUX CHOIX EXCLUSIFS :
  65 #
  66 #      A. ADAPTATION AVEC DES VARIANTES SUR LE MODE DE RAFFINEMENT/DERAFFINEMENT
  67 #         . RAFFINEMENT ET DERAFFINEMENT
  68 #         . RAFFINEMENT SEUL
  69 #         . DERAFFINEMENT SEUL
  70 #      B. INFORMATION SUR UN MAILLAGE
  71 #
  72            regles=( AU_MOINS_UN('RAFFINEMENT','DERAFFINEMENT','INFORMATION'),
  73                     EXCLUS('RAFFINEMENT','INFORMATION'),
  74                     EXCLUS('DERAFFINEMENT','INFORMATION'),),
  75            RAFFINEMENT      =SIMP(statut='f',typ='TXM',
  76                                  fr="Choix du mode de raffinement.",
  77                                  ang="Choice of refinement mode.",
  78                                  defaut="LIBRE",
  79                                  into=("LIBRE","UNIFORME","NON","NON-CONFORME","NON-CONFORME-INDICATEUR") ),
  80            DERAFFINEMENT   =SIMP(statut='f',typ='TXM',
  81                                  fr="Choix du mode de deraffinement.",
  82                                  ang="Choice of unrefinement mode.",
  83                                  defaut="LIBRE",
  84                                  into=("LIBRE","UNIFORME","NON") ),
  85            INFORMATION     =SIMP(statut='f',typ='TXM',
  86                                  fr="Information sur un maillage",
  87                                  ang="Information on a mesh",
  88                                  into=("OUI",) ),
  89 #
  90 # 3.2. LES CONTRAINTES :
  91 #
  92 # 3.2.1. POUR DE L'ADAPTATION LIBRE, IL FAUT :
  93 #      A. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
  94 #      B. LE NOM MED DU MAILLAGE D'ENTREE
  95 #      C. LE NOM MED DE L'INDICATEUR D'ERREUR
  96 #      D. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
  97 #      E. LA MISE A JOUR DE SOLUTION
  98 #      F. LE NOM MED DU MAILLAGE DE SORTIE
  99 #      REMARQUE : IL FAUT DES CRITERES, MAIS ON NE SAIT PAS LESQUELS
 100 #
 101 # 3.2.2. POUR DE L'ADAPTATION UNIFORME
 102 #          IL FAUT :
 103 #      A. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
 104 #      B. LE NOM MED DU MAILLAGE DE SORTIE
 105 #          IL NE FAUT PAS :
 106 #      A. LE NOM MED DE L'INDICATEUR D'ERREUR
 107 #      B. LE NOM DE LA COMPOSANTE DE L'INDICATEUR D'ERREUR
 108 #      C. LES CRITERES
 109 #      REMARQUE : A L'ITERATION 0, OU AUX ITERATIONS SUIVANTES SI MAJ DE SOLUTION,
 110 #                 IL FAUT LE NOM MED DU MAILLAGE D'ENTREE
 111 #
 112 # 3.2.3. POUR DE L'INFORMATION :
 113 #          IL FAUT :
 114 #      A. LE NOM MED DU MAILLAGE D'ENTREE
 115 #          IL NE FAUT PAS :
 116 #      A. LE NOM MED DE L'INDICATEUR D'ERREUR
 117 #      B. LE NOM DE LA COMPOSANTE DE L'INDICATEUR D'ERREUR
 118 #      C. LES CRITERES
 119 #      D. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
 120 #      E. LA MISE A JOUR DE SOLUTION
 121 #
 122            NOM_MED_MAILLAGE_N    =SIMP(statut='o',typ='TXM',
 123                                  fr="Nom MED du maillage en entrée",
 124                                  ang="MED name of the in-mesh",),
 125 #
 126            b_iteration_maj_champ =BLOC(condition = "( RAFFINEMENT != None ) or ( DERAFFINEMENT != None ) ",
 127                            fr="Nom MED du maillage en sortie, numero d'iteration et mise à jour de champs",
 128                            ang="MED name of the out-mesh, iteration rank and field updating",
 129                            NITER                =SIMP(statut='o',typ='I',
 130                            fr="Numéro d'itération avant l'adaptation.",
 131                            ang="Iteration number before adaptation." ),
 132                            NOM_MED_MAILLAGE_NP1 =SIMP(statut='o',typ='TXM',
 133                            fr="Nom MED du maillage en sortie",
 134                            ang="MED name of the out-mesh" ),
 135                            ) ,
 136 #
 137            b_indicateur_d_erreur  =BLOC(condition = "( RAFFINEMENT == 'LIBRE' ) or ( RAFFINEMENT == 'NON-CONFORME' ) or \
 138                                                      ( RAFFINEMENT == 'NON-CONFORME-INDICATEUR' ) or ( DERAFFINEMENT == 'LIBRE' ) ",
 139                            fr="Indicateur d'erreur",
 140                            ang="Error indicator",
 141                            NOM_MED  =SIMP(statut='o',typ='TXM',
 142                            fr="Nom MED de l'indicateur d'erreur.",
 143                            ang="MED name of error indicator.",),
 144                            COMPOSANTE  =SIMP(statut='o',typ='TXM',
 145                            fr="Nom de la composante de l'indicateur d'erreur retenue.",
 146                            ang="Name of the selected component of the error indicator.",),
 147                            NUME_ORDRE  =SIMP(statut='f',typ='I',
 148                            fr="Numero d'ordre de l'indicateur.",
 149                            ang="Rank number of the error indicator.",),
 150                            ) ,
 151 #
 152            b_critere_de_raffinement =BLOC( condition = "( RAFFINEMENT == 'LIBRE' ) or ( RAFFINEMENT == 'NON-CONFORME' ) or \
 153                                                         ( RAFFINEMENT == 'NON-CONFORME-INDICATEUR' ) ",
 154                            fr="Critère de raffinement.",
 155                            ang="Refinement threshold.",
 156                            regles=(UN_PARMI ( 'CRIT_RAFF_ABS', 'CRIT_RAFF_REL', 'CRIT_RAFF_PE' ),),
 157                            CRIT_RAFF_ABS   =SIMP(statut='f',typ='R',
 158                                                  fr="Critère absolu",
 159                                                  ang="Absolute threshold"  ),
 160                            CRIT_RAFF_REL   =SIMP(statut='f',typ='R',
 161                                                  fr="Critère relatif",
 162                                                  ang="Relative threshold" ),
 163                            CRIT_RAFF_PE    =SIMP(statut='f',typ='R',
 164                                                  fr="Pourcentage d'éléments",
 165                                                  ang="Percentage of elements" ),
 166                            ) ,
 167 #
 168            b_critere_de_deraffinement =BLOC ( condition = "( DERAFFINEMENT == 'LIBRE' ) ",
 169                            fr="Critère de déraffinement.",
 170                            ang="Unrefinement threshold.",
 171                            regles=(UN_PARMI ( 'CRIT_DERA_ABS', 'CRIT_DERA_REL', 'CRIT_DERA_PE' ),),
 172                            CRIT_DERA_ABS   =SIMP(statut='f',typ='R' ,
 173                                                  fr="Critère absolu",
 174                                                  ang="Absolute threshold" ),
 175                            CRIT_DERA_REL   =SIMP(statut='f',typ='R',
 176                                                  fr="Critère relatif",
 177                                                  ang="Relative threshold" ),
 178                            CRIT_DERA_PE    =SIMP(statut='f',typ='R',
 179                                                  fr="Pourcentage d'éléments",
 180                                                  ang="Percentage of elements" ),
 181                            ) ,
 182 #
 183            b_niveau_maximum =BLOC ( condition = " ( RAFFINEMENT == 'LIBRE' ) or ( RAFFINEMENT == 'NON-CONFORME' ) or \
 184                                                   ( RAFFINEMENT == 'NON-CONFORME-INDICATEUR' ) or ( RAFFINEMENT == 'UNIFORME' ) ",
 185                            NIVE_MAX        =SIMP(statut='f',typ='I',
 186                                                  fr="Niveau maximum de profondeur de raffinement",
 187                                                  ang="Maximum level for refinement"),
 188                            ) ,
 189 #
 190            b_niveau_minimum =BLOC ( condition = " ( DERAFFINEMENT == 'LIBRE' ) or ( DERAFFINEMENT == 'UNIFORME' ) ",
 191                            NIVE_MIN        =SIMP(statut='f',typ='I',
 192                                                  fr="Niveau minimum de déraffinement",
 193                                                  ang="Minimum level for unrefinement" ),
 194                            ) ,
 195 #
 196 # 3.3. Le suivi de frontiere eventuel :
 197 #
 198          NOM_MED_MAILLAGE_FRONTIERE =SIMP(statut='f',typ='TXM',
 199                            fr="Nom MED du maillage de la frontiere à suivre",
 200                            ang="MED name of the boundary mesh" ),
 201 #
 202          b_frontiere_1 =BLOC ( condition = "NOM_MED_MAILLAGE_FRONTIERE != None" ,
 203                            GROUP_MA        =SIMP(statut='f',typ=grma,validators=NoRepeat(),max='**',
 204                                                  fr="Groupes définissant la frontière",
 205                                                  ang="Groups which define the boundary" ),
 206                                ) ,
 207          fichier_frontiere=BLOC ( condition = "NOM_MED_MAILLAGE_FRONTIERE != None" ,
 208                           FIC_FRON = SIMP(statut='f',typ='TXM',
 209                                          fr="Nom du fichier MED contenant le maillage frontiere",
 210                                          ang="MED File including the boundary mesh" ),
 211                                ),
 212 #
 213          ),
 214 #
 215 # 4. L'ANALYSE DU MAILLAGE
 216 #
 217          ANALYSE         =FACT(statut='f',
 218                                fr="Analyse du maillage.",
 219                                ang="Mesh analysis.",
 220 #
 221 #    5 CHOIX NON EXCLUSIFS, AVEC DEUX VARIANTES (OUI/NON) :
 222 #    A. NOMBRE DES ELEMENTS
 223 #    B. QUALITE DES ELEMENTS
 224 #    C. INTERPENETRATION DES ELEMENTS
 225 #    D. CONNEXITE DU MAILLAGE
 226 #    E. TAILLE DES DIFFERENTS SOUS-DOMAINES
 227 #
 228            regles=(AU_MOINS_UN('NOMBRE','QUALITE','INTERPENETRATION','CONNEXITE','TAILLE'),),
 229 #
 230          NOMBRE          =SIMP(statut='f',typ='TXM',defaut="OUI",into=("OUI","NON"),
 231                           fr="Nombre de noeuds et éléments du maillage",
 232                           ang="Number of nodes and elements in the mesh" ),
 233 #
 234          QUALITE         =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 235                           fr="Qualité du maillage",
 236                           ang="Mesh quality" ),
 237 #
 238          INTERPENETRATION=SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 239                           fr="Controle de la non interpénétration des éléments.",
 240                           ang="Overlapping checking." ),
 241 #
 242          CONNEXITE       =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 243                           fr="Connexité du maillage.",
 244                           ang="Mesh connexity." ),
 245 #
 246          TAILLE          =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 247                           fr="Tailles des sous-domaines du maillage.",
 248                           ang="Sizes of mesh sub-domains." ),
 249 #
 250          ),
 251 #
 252 )  ;