Homard/homard_cata_STA7.py

   1
   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                                  into=("LIBRE","UNIFORME","NON","NON-CONFORME","NON-CONFORME-INDICATEUR") ),
  79            DERAFFINEMENT   =SIMP(statut='f',typ='TXM',
  80                                  fr="Choix du mode de deraffinement.",
  81                                  ang="Choice of unrefinement mode.",
  82                                  into=("LIBRE","UNIFORME","NON") ),
  83            INFORMATION     =SIMP(statut='f',typ='TXM',
  84                                  fr="Information sur un maillage",
  85                                  ang="Information on a mesh",
  86                                  into=("OUI",) ),
  87 #
  88 # 3.2. LES CONTRAINTES :
  89 #
  90 # 3.2.1. POUR DE L'ADAPTATION LIBRE, IL FAUT :
  91 #      A. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
  92 #      B. LE NOM MED DU MAILLAGE D'ENTREE
  93 #      C. LE NOM MED DE L'INDICATEUR D'ERREUR
  94 #      D. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
  95 #      E. LA MISE A JOUR DE SOLUTION
  96 #      F. LE NOM MED DU MAILLAGE DE SORTIE
  97 #      REMARQUE : IL FAUT DES CRITERES, MAIS ON NE SAIT PAS LESQUELS
  98 #
  99 # 3.2.2. POUR DE L'ADAPTATION UNIFORME
 100 #          IL FAUT :
 101 #      A. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
 102 #      B. LE NOM MED DU MAILLAGE DE SORTIE
 103 #          IL NE FAUT PAS :
 104 #      A. LE NOM MED DE L'INDICATEUR D'ERREUR
 105 #      B. LE NOM DE LA COMPOSANTE DE L'INDICATEUR D'ERREUR
 106 #      C. LES CRITERES
 107 #      REMARQUE : A L'ITERATION 0, OU AUX ITERATIONS SUIVANTES SI MAJ DE SOLUTION,
 108 #                 IL FAUT LE NOM MED DU MAILLAGE D'ENTREE
 109 #
 110 # 3.2.3. POUR DE L'INFORMATION :
 111 #          IL FAUT :
 112 #      A. LE NOM MED DU MAILLAGE D'ENTREE
 113 #          IL NE FAUT PAS :
 114 #      A. LE NOM MED DE L'INDICATEUR D'ERREUR
 115 #      B. LE NOM DE LA COMPOSANTE DE L'INDICATEUR D'ERREUR
 116 #      C. LES CRITERES
 117 #      D. LE NUMERO D'ITERATION DU MAILLAGE DE DEPART
 118 #      E. LA MISE A JOUR DE SOLUTION
 119 #
 120            NOM_MED_MAILLAGE_N    =SIMP(statut='o',typ='TXM',
 121                                  fr="Nom MED du maillage en entrée",
 122                                  ang="MED name of the in-mesh",),
 123 #
 124            b_iteration_maj_champ =BLOC(condition = "( RAFFINEMENT != None ) or ( DERAFFINEMENT != None ) ",
 125                            fr="Nom MED du maillage en sortie, numero d'iteration et mise à jour de champs",
 126                            ang="MED name of the out-mesh, iteration rank and field updating",
 127                            NITER                =SIMP(statut='o',typ='I',
 128                            fr="Numéro d'itération avant l'adaptation.",
 129                            ang="Iteration number before adaptation." ),
 130                            NOM_MED_MAILLAGE_NP1 =SIMP(statut='o',typ='TXM',
 131                            fr="Nom MED du maillage en sortie",
 132                            ang="MED name of the out-mesh" ),
 133                            ) ,
 134 #
 135            b_indicateur_d_erreur  =BLOC(condition = "( RAFFINEMENT == 'LIBRE' ) or ( RAFFINEMENT == 'NON-CONFORME' ) or \
 136                                                      ( RAFFINEMENT == 'NON-CONFORME-INDICATEUR' ) or ( DERAFFINEMENT == 'LIBRE' ) ",
 137                            fr="Indicateur d'erreur",
 138                            ang="Error indicator",
 139                            NOM_MED  =SIMP(statut='o',typ='TXM',
 140                            fr="Nom MED de l'indicateur d'erreur.",
 141                            ang="MED name of error indicator.",),
 142                            COMPOSANTE  =SIMP(statut='o',typ='TXM',
 143                            fr="Nom de la composante de l'indicateur d'erreur retenue.",
 144                            ang="Name of the selected component of the error indicator.",),
 145                            NUME_ORDRE  =SIMP(statut='f',typ='I',
 146                            fr="Numero d'ordre de l'indicateur.",
 147                            ang="Rank number of the error indicator.",),
 148                            ) ,
 149 #
 150            b_critere_de_raffinement =BLOC( condition = "( RAFFINEMENT == 'LIBRE' ) or ( RAFFINEMENT == 'NON-CONFORME' ) or \
 151                                                         ( RAFFINEMENT == 'NON-CONFORME-INDICATEUR' ) ",
 152                            fr="Critère de raffinement.",
 153                            ang="Refinement threshold.",
 154                            regles=(UN_PARMI ( 'CRIT_RAFF_ABS', 'CRIT_RAFF_REL', 'CRIT_RAFF_PE' ),),
 155                            CRIT_RAFF_ABS   =SIMP(statut='f',typ='R',
 156                                                  fr="Critère absolu",
 157                                                  ang="Absolute threshold"  ),
 158                            CRIT_RAFF_REL   =SIMP(statut='f',typ='R',
 159                                                  fr="Critère relatif",
 160                                                  ang="Relative threshold" ),
 161                            CRIT_RAFF_PE    =SIMP(statut='f',typ='R',
 162                                                  fr="Pourcentage d'éléments",
 163                                                  ang="Percentage of elements" ),
 164                            ) ,
 165 #
 166            b_critere_de_deraffinement =BLOC ( condition = "( DERAFFINEMENT == 'LIBRE' ) ",
 167                            fr="Critère de déraffinement.",
 168                            ang="Unrefinement threshold.",
 169                            regles=(UN_PARMI ( 'CRIT_DERA_ABS', 'CRIT_DERA_REL', 'CRIT_DERA_PE' ),),
 170                            CRIT_DERA_ABS   =SIMP(statut='f',typ='R' ,
 171                                                  fr="Critère absolu",
 172                                                  ang="Absolute threshold" ),
 173                            CRIT_DERA_REL   =SIMP(statut='f',typ='R',
 174                                                  fr="Critère relatif",
 175                                                  ang="Relative threshold" ),
 176                            CRIT_DERA_PE    =SIMP(statut='f',typ='R',
 177                                                  fr="Pourcentage d'éléments",
 178                                                  ang="Percentage of elements" ),
 179                            ) ,
 180 #
 181            b_niveau_maximum =BLOC ( condition = " ( RAFFINEMENT == 'LIBRE' ) or ( RAFFINEMENT == 'NON-CONFORME' ) or \
 182                                                   ( RAFFINEMENT == 'NON-CONFORME-INDICATEUR' ) or ( RAFFINEMENT == 'UNIFORME' ) ",
 183                            NIVE_MAX        =SIMP(statut='f',typ='I',
 184                                                  fr="Niveau maximum de profondeur de raffinement",
 185                                                  ang="Maximum level for refinement"),
 186                            ) ,
 187 #
 188            b_niveau_minimum =BLOC ( condition = " ( DERAFFINEMENT == 'LIBRE' ) or ( DERAFFINEMENT == 'UNIFORME' ) ",
 189                            NIVE_MIN        =SIMP(statut='f',typ='I',
 190                                                  fr="Niveau minimum de déraffinement",
 191                                                  ang="Minimum level for unrefinement" ),
 192                            ) ,
 193 #
 194 # 3.3. Le suivi de frontiere eventuel :
 195 #
 196          NOM_MED_MAILLAGE_FRONTIERE =SIMP(statut='f',typ='TXM',
 197                            fr="Nom MED du maillage de la frontiere à suivre",
 198                            ang="MED name of the boundary mesh" ),
 199 #
 200          b_frontiere_1 =BLOC ( condition = "NOM_MED_MAILLAGE_FRONTIERE != None" ,
 201                            GROUP_MA        =SIMP(statut='f',typ=grma,validators=NoRepeat(),max='**',
 202                                                  fr="Groupes définissant la frontière",
 203                                                  ang="Groups which define the boundary" ),
 204                                ) ,
 205          fichier_frontiere=BLOC ( condition = "NOM_MED_MAILLAGE_FRONTIERE != None" ,
 206                           FIC_FRON = SIMP(statut='f',typ='TXM',
 207                                          fr="Nom du fichier MED contenant le maillage frontiere",
 208                                          ang="MED File including the boundary mesh" ),
 209                                ),
 210 #
 211          ),
 212 #
 213 # 4. L'ANALYSE DU MAILLAGE
 214 #
 215          ANALYSE         =FACT(statut='f',
 216                                fr="Analyse du maillage.",
 217                                ang="Mesh analysis.",
 218 #
 219 #    5 CHOIX NON EXCLUSIFS, AVEC DEUX VARIANTES (OUI/NON) :
 220 #    A. NOMBRE DES ELEMENTS
 221 #    B. QUALITE DES ELEMENTS
 222 #    C. INTERPENETRATION DES ELEMENTS
 223 #    D. CONNEXITE DU MAILLAGE
 224 #    E. TAILLE DES DIFFERENTS SOUS-DOMAINES
 225 #
 226            regles=(AU_MOINS_UN('NOMBRE','QUALITE','INTERPENETRATION','CONNEXITE','TAILLE'),),
 227 #
 228          NOMBRE          =SIMP(statut='f',typ='TXM',defaut="OUI",into=("OUI","NON"),
 229                           fr="Nombre de noeuds et éléments du maillage",
 230                           ang="Number of nodes and elements in the mesh" ),
 231 #
 232          QUALITE         =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 233                           fr="Qualité du maillage",
 234                           ang="Mesh quality" ),
 235 #
 236          INTERPENETRATION=SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 237                           fr="Controle de la non interpénétration des éléments.",
 238                           ang="Overlapping checking." ),
 239 #
 240          CONNEXITE       =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 241                           fr="Connexité du maillage.",
 242                           ang="Mesh connexity." ),
 243 #
 244          TAILLE          =SIMP(statut='f',typ='TXM',defaut="NON",into=("OUI","NON"),
 245                           fr="Tailles des sous-domaines du maillage.",
 246                           ang="Sizes of mesh sub-domains." ),
 247 #
 248          ),
 249 #
 250 )  ;