Syrthes/cata_syrthes.py

   1 # debut entete
   2 import Accas
   3 from Accas import *
   4 #
   5 JdC = JDC_CATA(code='SYRTHES',
   6                execmodul=None,
   7                niveaux=(NIVEAU(nom='generalites',label='Généralités',valide_vide=0),
   8                         NIVEAU(nom='gest_calc',label='Gestion du calcul',valide_vide=0,actif="generalites == 1"),
   9                         NIVEAU(nom='gest_sort',label='Gestion des sorties'),
  10                         NIVEAU(nom='choix_num',label='Choix numériques',valide_vide=0),\r
  11                         NIVEAU(nom='prop_phys',label='Propriétés physiques',valide_vide=0,actif='choix_num==1'),\r
  12                         NIVEAU(nom='cond_init',label='Conditions initiales',valide_vide=0,actif='prop_phys==1'),\r
  13                         NIVEAU(nom='cond_limi',label='Conditions aux limites',valide_vide=0,actif='cond_init==1'),\r
  14                         NIVEAU(nom='flux_volu',label='Flux volumiques',valide_vide=0,actif='cond_limi==1'),\r
  15                         NIVEAU(nom='soli_mouv',label='Solides en mouvement'),\r
  16                         NIVEAU(nom='couplages',label='Couplages'),\r
  17                         NIVEAU(nom='rayonnement',label='Rayonnement',valide_vide=0,actif='flux_volu==1'),\r
  18                         NIVEAU(nom='cond_lim_ray',label='Conditions limites rayonnement',valide_vide=0,actif='rayonnement==1'),\r
  19                         NIVEAU(nom='in_out_ray',label='Entrées/Sorties rayonnement',valide_vide=0,actif='rayonnement==1'),\r
  20                         NIVEAU(nom='fichiers',label='Fichiers pour le calcul'))
  21                )
  22
  23
  24 GENERALITES = PROC(nom='GENERALITES',op=10,\r
  25                    fr="Définition des principales options du calcul",\r
  26                    niveau='generalites',
  27                    TYPE_DE_CALCUL=SIMP(typ='TXM',\r
  28                                        statut='o',\r
  29                                        into=('Syrthes thermique solide seule','Syrthes couplé'),\r
  30                                        defaut='Syrthes thermique solide seule',\r
  31                                        position='global_jdc'),
  32                    TITRE = SIMP(statut='f',typ='TXM'),
  33                    DIMENSION = SIMP(statut='o',\r
  34                                     typ='I',\r
  35                                     into=(2,3),\r
  36                                     defaut=2,\r
  37                                     position='global_jdc'),
  38                    b_axys = BLOC(condition="TYPE_DE_CALCUL=='Syrthes thermique solide seule'",\r
  39                                  AXE_SYMETRIE = SIMP(statut='o',\r
  40                                                      typ='TXM',\r
  41                                                      into=('AUCUN','OX','OY'),\r
  42                                                      defaut='AUCUN')));\r
  43 \r
  44 GESTION_DU_CALCUL = PROC(nom='GESTION_DU_CALCUL',op=10,\r
  45                          fr="Gestion des pas de temps et éventuellement du couplage fluide/solide",\r
  46                          niveau='gest_calc',\r
  47                          CALCUL_SYRTHES_SEUL = BLOC(condition="TYPE_DE_CALCUL == 'Syrthes thermique solide seule'",\r
  48                                                     PAS_DE_TEMPS_SOLIDE = SIMP(typ='I',statut='o'),\r
  49                                                     NB_PAS_DE_TEMPS = SIMP(typ='I',statut='o'),\r
  50                                                     SUITE_DE_CALCUL = SIMP(typ='TXM',statut='o',into=('OUI','NON'),defaut='NON')),\r
  51                          CALCUL_SYRTHES_COUPLE = BLOC(condition="TYPE_DE_CALCUL =='Syrthes couplé'",\r
  52                                                       CORRESPONDANCE_MAILLAGES_F_S = SIMP(statut='o',\r
  53                                                                                           typ='TXM',\r
  54                                                                                           into=('calcul','lecture sur fichier','calcul et écriture sur fichier'),\r
  55                                                                                           defaut='calcul et écriture sur fichier'),\r
  56                                                       PAS_DE_TEMPS_SOLIDE = SIMP(typ='TXM',\r
  57                                                                                  statut='o',\r
  58                                                                                  into=('fixé indépendamment','égal à celui du fluide'),\r
  59                                                                                  defaut='égal à celui du fluide'),\r
  60                                                       CHOIX_FREQUENCE_CALCUL_SOLIDE = SIMP(typ='TXM',\r
  61                                                                                      statut='o',\r
  62                                                                                      into=('calcul solide à chaque pas de temps',\r
  63                                                                                            'calcul solide tous les n pas de temps fluide',\r
  64                                                                                            'calcul fluide tous les n pas de temps solide'),\r
  65                                                                                      defaut='calcul solide à chaque pas de temps'),\r
  66                                                       frequence_calcul = BLOC(condition="CHOIX_FREQUENCE_CALCUL_SOLIDE != 'calcul solide à chaque pas de temps'",\r
  67                                                                               FREQUENCE = SIMP(statut='o',\r
  68                                                                                                typ='I')))\r
  69                                                       \r
  70                        );\r
  71 \r
  72 GESTION_DES_SORTIES = PROC(nom="GESTION_DES_SORTIES",op=10,\r
  73                            fr="Gestion des résultats du code",\r
  74                            niveau='gest_sort',\r
  75                            IMPRESSIONS_SOLIDE=SIMP(statut='o',typ='TXM',\r
  76                                                    into=('aucune','minimales','normales','détaillées'),\r
  77                                                    defaut='normales'),\r
  78                            RESULTATS_DOMAINE_SOLIDE=FACT(statut='o',\r
  79                                                          FICHIER_CHRONOLOGIQUE=SIMP(typ='TXM',statut='o',\r
  80                                                                                     into=('aucune sortie','tous les n pas de temps'),\r
  81                                                                                     defaut='aucune sortie'),\r
  82                                                          pas_de_temps=BLOC(condition="FICHIER_CHRONOLOGIQUE=='tous les n pas de temps'",\r
  83                                                                            PAS_DE_TEMPS=SIMP(statut='o',typ='I')),\r
  84                                                          FICHIER_LISTING=SIMP(statut='o',typ='TXM',\r
  85                                                                               into=('aucune sortie',),\r
  86                                                                               defaut='aucune sortie'),\r
  87                                                          HISTORIQUE_SOLIDE=SIMP(statut='o',typ='TXM',\r
  88                                                                                 into=('OUI','NON'),\r
  89                                                                                 defaut='NON'),\r
  90                                                          historique=BLOC(condition="HISTORIQUE_SOLIDE=='OUI'",\r
  91                                                                          NOMBRE_CAPTEURS=SIMP(statut='o',typ='I'),\r
  92                                                                          FREQUENCE_ECRITURE=SIMP(statut='o',typ='I',defaut=1),\r
  93                                                                          LISTE_NOEUDS = FACT(statut='o',max='**',\r
  94                                                                                              NOEUDS=SIMP(statut='o',max='**',typ='I'),\r
  95                                                                                              COMMENTAIRE=SIMP(statut='o',typ='TXM',defaut=""))),\r
  96                                                          ),\r
  97                            RESULTATS_DOMAINE_FLUIDE=BLOC(condition="TYPE_DE_CALCUL=='Syrthes couplé'",\r
  98                                                          ECRITURE_MAILLAGE_PEAU_FLUIDE=SIMP(statut='o',into=('OUI','NON'),\r
  99                                                                                             defaut='NON',typ='TXM'),\r
 100                                                          ECRITURE_RESULTATS_PEAU_FLUIDE=SIMP(statut='o',into=('OUI','NON'),\r
 101                                                                                              defaut='NON',typ='TXM'),\r
 102                                                          ECRITURE_CHRONO_PEAU_FLUIDE=SIMP(statut='o',into=('OUI','NON'),\r
 103                                                                                           defaut='NON',typ='TXM'))\r
 104                            );\r
 105 \r
 106 CHOIX_NUMERIQUES=PROC(nom="CHOIX_NUMERIQUES",op=10,\r
 107                       fr="Modification des options du solveur",\r
 108                       niveau='choix_num',\r
 109                       SOLVEUR_SOLIDE = FACT(statut='o',\r
 110                                             NOMBRE_ITERATIONS=SIMP(statut='o',typ='I'),\r
 111                                             PRECISION=SIMP(statut='o',typ='R'))\r
 112                       );\r
 113 \r
 114 PROPRIETES_PHYSIQUES=PROC(nom="PROPRIETES_PHYSIQUES",op=10,\r
 115                           fr="Définition des propriétés physiques des matériaux",\r
 116                           niveau='prop_phys',\r
 117                           DEFINITION_PAR = SIMP(statut='o',typ='TXM',into=('NOEUD','ELEMENT','NOEUD PAR ELEMENT'),\r
 118                                                 defaut='ELEMENT'),\r
 119                           ISOTROPIE = SIMP(statut='o', typ='TXM', into=('Matériaux isotropes',\r
 120                                                                         'Matériaux anisotropes',\r
 121                                                                         'Matériaux orthotropes'),\r
 122                                            defaut='Matériaux isotropes'),\r
 123                           proprietes_des_materiaux_isotropes = BLOC(condition="ISOTROPIE=='Mat\351riaux isotropes'",                          \r
 124                               DEFINITION_MASSE_VOLUMIQUE=FACT(statut='o',max='**',\r
 125                                                           MASSE_VOLUMIQUE = SIMP(statut='o',typ='R'),\r
 126                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 127                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 128                               DEFINITION_CHALEUR_SPECIFIQUE=FACT(statut='o',max='**',\r
 129                                                           CHALEUR_SPECIFIQUE = SIMP(statut='o',typ='R'),\r
 130                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 131                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 132                               DEFINITION_CONDUCTIVITE_THERMIQUE=FACT(statut='o',max='**',\r
 133                                                           CONDUCTIVITE_THERMIQUE = SIMP(statut='o',typ='R'),\r
 134                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 135                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM'))\r
 136                                                                     ),\r
 137                           proprietes_des_materiaux_orthotropes_2D = BLOC(condition="ISOTROPIE=='Mat\351riaux orthotropes' and DIMENSION==2",\r
 138                               DEFINITION_MASSE_VOLUMIQUE=FACT(statut='o',max='**',\r
 139                                                           MASSE_VOLUMIQUE = SIMP(statut='o',typ='R'),\r
 140                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 141                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 142                               DEFINITION_CHALEUR_SPECIFIQUE=FACT(statut='o',max='**',\r
 143                                                           CHALEUR_SPECIFIQUE = SIMP(statut='o',typ='R'),\r
 144                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 145                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 146                               DEFINITION_CONDUCTIVITE_THERMIQUE=FACT(statut='o',max='**',\r
 147                                                           CONDUCTIVITE_EN_X = SIMP(statut='o',typ='R'),\r
 148                                                           CONDUCTIVITE_EN_Y = SIMP(statut='o',typ='R'),\r
 149                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 150                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM'))\r
 151                                                                          ),\r
 152                           proprietes_des_materiaux_orthotropes_3D = BLOC(condition="ISOTROPIE=='Mat\351riaux orthotropes' and DIMENSION==3",\r
 153                               DEFINITION_MASSE_VOLUMIQUE=FACT(statut='o',max='**',\r
 154                                                           MASSE_VOLUMIQUE = SIMP(statut='o',typ='R'),\r
 155                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 156                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 157                               DEFINITION_CHALEUR_SPECIFIQUE=FACT(statut='o',max='**',\r
 158                                                           CHALEUR_SPECIFIQUE = SIMP(statut='o',typ='R'),\r
 159                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 160                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 161                               DEFINITION_CONDUCTIVITE_THERMIQUE=FACT(statut='o',max='**',\r
 162                                                           CONDUCTIVITE_EN_X = SIMP(statut='o',typ='R'),\r
 163                                                           CONDUCTIVITE_EN_Y = SIMP(statut='o',typ='R'),\r
 164                                                           CONDUCTIVITE_EN_Z = SIMP(statut='o',typ='R'),\r
 165                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 166                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM'))\r
 167                                                                          ),\r
 168                           proprietes_des_materiaux_anisotropes_2D = BLOC(condition="ISOTROPIE=='Mat\351riaux anisotropes' and DIMENSION==2",\r
 169                               DEFINITION_MASSE_VOLUMIQUE=FACT(statut='o',max='**',\r
 170                                                           MASSE_VOLUMIQUE = SIMP(statut='o',typ='R'),\r
 171                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 172                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 173                               DEFINITION_CHALEUR_SPECIFIQUE=FACT(statut='o',max='**',\r
 174                                                           CHALEUR_SPECIFIQUE = SIMP(statut='o',typ='R'),\r
 175                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 176                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 177                               DEFINITION_CONDUCTIVITE_THERMIQUE=FACT(statut='o',max='**',\r
 178                                                           CONDUCTIVITE_EN_X = SIMP(statut='o',typ='R'),\r
 179                                                           CONDUCTIVITE_EN_Y = SIMP(statut='o',typ='R'),\r
 180                                                           ANGLE_SUIVANT_Z = SIMP(statut='o',typ='R'),\r
 181                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 182                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM'))\r
 183                                                                          ),\r
 184                           proprietes_des_materiaux_anisotropes_3D = BLOC(condition="ISOTROPIE=='Mat\351riaux anisotropes' and DIMENSION==3",\r
 185                               DEFINITION_MASSE_VOLUMIQUE=FACT(statut='o',max='**',\r
 186                                                           MASSE_VOLUMIQUE = SIMP(statut='o',typ='R'),\r
 187                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 188                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 189                               DEFINITION_CHALEUR_SPECIFIQUE=FACT(statut='o',max='**',\r
 190                                                           CHALEUR_SPECIFIQUE = SIMP(statut='o',typ='R'),\r
 191                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 192                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM')),\r
 193                               DEFINITION_CONDUCTIVITE_THERMIQUE=FACT(statut='o',max='**',\r
 194                                                           CONDUCTIVITE_EN_X = SIMP(statut='o',typ='R'),\r
 195                                                           CONDUCTIVITE_EN_Y = SIMP(statut='o',typ='R'),\r
 196                                                           CONDUCTIVITE_EN_Z = SIMP(statut='o',typ='R'),\r
 197                                                           ANGLE_SUIVANT_X = SIMP(statut='o',typ='R'),\r
 198                                                           ANGLE_SUIVANT_Y = SIMP(statut='o',typ='R'),\r
 199                                                           ANGLE_SUIVANT_Z = SIMP(statut='o',typ='R'),\r
 200                                                           REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 201                                                           COMMENTAIRES=SIMP(statut='f',typ='TXM'))\r
 202                                                                          ),\r
 203                           );\r
 204 \r
 205 CONDITIONS_INITIALES=PROC(nom="CONDITIONS_INITIALES",op=10,\r
 206                           fr="Définition des conditions initiales",\r
 207                           niveau='cond_init',\r
 208                           Temperature_initiale = FACT(statut='o',max='**',\r
 209                                                       TEMPERATURE = SIMP(typ='R',statut='o',defaut=20),\r
 210                                                       REFERENCES = SIMP(statut='o',typ='I',max='**',defaut=-1),\r
 211                                                       COMMENTAIRES = SIMP(statut='f',typ='TXM')\r
 212                                                       )\r
 213                           );\r
 214 \r
 215 \r
 216 CONDITIONS_LIMITES = PROC(nom='CONDITIONS_LIMITES',op=10,\r
 217                           niveau='cond_limi',\r
 218                           fr="Définition des conditions aux limites",\r
 219                           NOEUD_OU_FACE = SIMP(statut='o',typ='TXM',into=('NOEUD','FACE'),defaut='FACE'),\r
 220                           FLUX_SURFACIQUE=FACT(max='**',\r
 221                                                FLUX=SIMP(statut='o',typ='R',defaut=0),\r
 222                                                REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 223                                                COMMENTAIRES = SIMP(statut='f',typ='TXM')\r
 224                                                ),\r
 225                           DIRICHLET=FACT(max='**',\r
 226                                          TEMPERATURE=SIMP(statut='o',typ='R',defaut=20),\r
 227                                          REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 228                                          COMMENTAIRES = SIMP(statut='f',typ='TXM')\r
 229                                          ),\r
 230 \r
 231                           COEFFICIENT_ECHANGE=FACT(max='**',\r
 232                                                    TEMPERATURE=SIMP(statut='o',typ='R',defaut=20),\r
 233                                                    COEFFICIENT = SIMP(statut='o',typ='R',defaut=0),\r
 234                                                    REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 235                                                    COMMENTAIRES = SIMP(statut='f',typ='TXM')\r
 236                                                   ),\r
 237                           RESISTANCE_CONTACT=FACT(max='**',\r
 238                                                   RESISTANCE = SIMP(statut='o',typ='R',defaut=0),\r
 239                                                   REFERENCES = SIMP(statut='o',typ='I',max='**'),\r
 240                                                   COMMENTAIRES = SIMP(statut='f',typ='TXM')\r
 241                                                  ),\r
 242                          );
 243