Cuve2dg/Cuve2dg_Cata_V1.py

   1 # -*- coding: utf-8 -*-
   2
   3 # --------------------------------------------------
   4 # debut entete
   5 # --------------------------------------------------
   6
   7 import Accas
   8 from Accas import *
   9
  10 class loi ( ASSD ) : pass
  11 class variable ( ASSD ) : pass
  12
  13
  14 #CONTEXT.debug = 1
  15 JdC = JDC_CATA ( code = 'CUVE1D-DEFAILLGLOB',
  16                  execmodul = None,
  17                  regles = ( AU_MOINS_UN ('OPTIONS','DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'),
  18                             AU_MOINS_UN ( 'FIN' ),
  19                             A_CLASSER ( ('OPTIONS', 'DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'),'FIN')
  20                           )
  21                  ) # Fin JDC_CATA
  22
  23 # --------------------------------------------------
  24 # fin entete
  25 # --------------------------------------------------
  26
  27
  28
  29
  30
  31
  32 #================================
  33 # 1. Definition des OPTIONS
  34 #================================
  35
  36 # Nota : les variables de type OPER doivent etre en majuscules !
  37 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
  38
  39 OPTIONS = OPER ( nom = "OPTIONS",
  40                  sd_prod = loi,
  41                  op = 68,
  42                  fr = "Definitions des options",
  43
  44 #===
  45 # Liste des paramètres
  46 #===
  47
  48   INCRTPS = SIMP ( statut = "o",
  49                    typ = "I",
  50                    defaut = "1",
  51                    max = 1,
  52                    val_max = 100,
  53                    fr = "Increment temporel (=1 pour calcul deterministe)",
  54                   ),
  55
  56   DTPREC = SIMP ( statut = "o",
  57                   typ = "R",
  58                   defaut = "0.1",
  59                   max = 1,
  60                   val_max = 1.,
  61                   fr = "Increment maximum d'evolution de la temperature par noeud et par instant (°C)",
  62                 ),
  63
  64   DTARCH = SIMP ( statut = "o",
  65                   typ = "R",
  66                   defaut = "1000.",
  67                   max = 1,
  68                   val_max = 1000.,
  69                   fr = "Increment maximum de temps pour l'affichage (s)",
  70                 ),
  71
  72   NBO = SIMP ( statut = "o",
  73                typ = "R",
  74                max=1,
  75                val_max = 1000.,
  76                fr = "Nombre de noeuds a considerer dans le maillage interne",
  77               ),
  78
  79   Liste_instants = SIMP ( statut = "o",
  80                           typ = "R",
  81                           max = "**",
  82                           fr = "Liste des instants ",
  83                         ),
  84
  85 ) # Fin OPER OPTIONS
  86
  87 #================================
  88 # 2. Caracteristiques du DEFAUT
  89 #================================
  90
  91 # Nota : les variables de type OPER doivent etre en majuscules !
  92 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
  93 DEFAUT = OPER ( nom = "DEFAUT",
  94                 sd_prod = loi,
  95                 op = 68,
  96                 fr = "Caracteristiques du defaut",
  97
  98 #===
  99 # Liste des paramètres
 100 #===
 101
 102   TYPEDEF = SIMP ( statut = "o", typ = "TXM",
 103                    into = ( "DSR",
 104                             "DD",
 105                           ),
 106                    #defaut = "DSR",
 107                    fr = "Type de defaut : sous revetement ou debouchant",
 108                   ),
 109
 110 #====
 111 # Definition des parametres selon le type du defaut
 112 #====
 113
 114   Parametres_DSR = BLOC ( condition = " TYPEDEF in ( 'DSR', ) ",
 115
 116                   ORIEDEF = SIMP ( statut = "o",
 117                                    typ = "TXM",
 118                                    into = ( "LONGITUD", "CIRCONF" ),
 119                                    #defaut = "LONGITUD",
 120                                    fr = "Orientation du defaut : longitudinale ou circonferentielle",
 121                                  ),
 122
 123                   PROFDEF = SIMP ( statut = "o",
 124                                    typ = "R",
 125                                    #defaut = "0.006",
 126                                    max = 1,
 127                                    val_max = 1.,
 128                                    fr = "Profondeur radiale du defaut (m)",
 129                                  ),
 130
 131                   OPTLONG = SIMP ( statut = "o",
 132                                    typ = "TXM",
 133                                    into = ( "VALEUR", "RAPPORT" ),
 134                                    #defaut = "VALEUR",
 135                                    fr = "Option pour caracteriser la longueur du defaut : soit par valeur, soit par un rapport LONG/PROF",
 136                                  ),
 137
 138                     Option_Valeur = BLOC ( condition = "OPTLONG in ( 'VALEUR', ) ",
 139
 140                             LONGDEF = SIMP ( statut = "o",
 141                                              typ = "R",
 142                                              #defaut = "0.060",
 143                                              max = 1,
 144                                              val_max = 1.,
 145                                              fr = "Longueur du defaut sous revetement (m)",
 146                                             ),
 147
 148                                           ), # Fin BLOC Option_Valeur
 149
 150                     Option_Rapport = BLOC ( condition = "OPTLONG in ( 'RAPPORT', ) ",
 151
 152                             LONGSURPROF = SIMP ( statut = "o",
 153                                                  typ = "R",
 154                                                  #defaut = "6.",
 155                                                  max = 1,
 156                                                  val_max = 100.,
 157                                                  fr = "Rapport longueur/profondeur du defaut sous revetement",
 158                                                 ),
 159
 160                                            ), # Fin BLOC Option_Rapport
 161
 162                   DECADEF = SIMP ( statut = "o",
 163                                    typ = "R",
 164                                    #defaut = "-0.00001",
 165                                    fr = "Decalage radial du defaut sous revetement (m)",
 166                                   ),
 167
 168                   ANGLDEF = SIMP ( statut = "o",
 169                                    typ = "R",
 170                                    defaut = "0.",
 171                                    fr = "Coordonnee angulaire du defaut (degres)",
 172                                   ),
 173
 174                   ALTIDEF = SIMP ( statut = "o",
 175                                    typ = "R",
 176                                    defaut = "2.",
 177                                    fr = "Altitude du defaut (m)",
 178                                   ),
 179
 180                   POINDEF = SIMP ( statut = "o",
 181                                    typ = "TXM",
 182                                    into = ( "A", "B" ),
 183                                    defaut = "A",
 184                                    fr = "Choix du point considere du defaut sous revetement",
 185                                   ),
 186
 187                   ARRETFISSURE = SIMP ( statut = "o",
 188                                         typ = "TXM",
 189                                         into = ( "OUI", "NON" ),
 190                                         defaut = "NON",
 191                                         fr = "Prise en compte de l arret de fissure",
 192                                        ),
 193
 194                   INCRDEF = SIMP ( statut = "o",
 195                                    typ = "R",
 196                                    defaut = "0.005",
 197                                    fr = "Increment de la taille de fissure (m)",
 198                                   ),
 199
 200                   CORRECPLASTIC = SIMP ( statut = "o",
 201                                          typ = "TXM",
 202                                          into = ( "OUI", "NON" ),
 203                                          defaut = "NON",
 204                                          fr = "Prise en compte de la correction plastique BETA ",
 205                                         ),
 206
 207                          ), # Fin BLOC Parametres_DSR
 208
 209   Parametres_DD = BLOC ( condition = " TYPEDEF in ( 'DD', ) ",
 210
 211                   ORIEDEF = SIMP ( statut = "o",
 212                                    typ = "TXM",
 213                                    into = ( "LONGITUD", "CIRCONF" ),
 214                                    #defaut = "LONGITUD",
 215                                    fr = "Orientation du defaut : longitudinale ou circonferentielle",
 216                                   ),
 217
 218                   PROFDEF = SIMP ( statut = "o",
 219                                    typ = "R",
 220                                    #defaut = "0.006",
 221                                    max = 1,
 222                                    val_max = 1.,
 223                                    fr = "Profondeur radiale du defaut (m)",
 224                                   ),
 225
 226                   ANGLDEF = SIMP ( statut = "o",
 227                                    typ = "R",
 228                                    defaut = "0.",
 229                                    fr = "Coordonnee angulaire du defaut (degres)",
 230                                   ),
 231
 232                   ALTIDEF = SIMP ( statut = "o",
 233                                    typ = "R",
 234                                    defaut = "2.",
 235                                    fr = "Altitude du defaut (m)",
 236                                   ),
 237
 238                   ARRETFISSURE = SIMP ( statut = "o",
 239                                         typ = "TXM",
 240                                         into = ( "OUI", "NON" ),
 241                                         defaut = "NON",
 242                                         fr = "Prise en compte de l arret de fissure",
 243                                        ),
 244
 245                   INCRDEF = SIMP ( statut = "o",
 246                                    typ = "R",
 247                                    defaut = "0.005",
 248                                    fr = "Increment de la taille de fissure (m)",
 249                                   ),
 250
 251                   IRWIN = SIMP ( statut = "o",
 252                                  typ = "TXM",
 253                                  into = ( "OUI", "NON" ),
 254                                  defaut = "NON",
 255                                  fr = "Prise en compte de la correction plastique d'Irwin ",
 256                                 ),
 257
 258                   CORRECPLASTIC = SIMP ( statut = "o",
 259                                          typ = "TXM",
 260                                          into = ( "OUI", "NON" ),
 261                                          defaut = "NON",
 262                                          fr = "Prise en compte de la correction plastique BETA ",
 263                                         ),
 264
 265                         ), # Fin BLOC Parametres_DD
 266
 267 ) # Fin OPER DEFAUT
 268
 269
 270 #================================
 271 # 3. Caracteristiques de la CUVE
 272 #================================
 273
 274 # Nota : les variables de type OPER doivent etre en majuscules !
 275 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
 276 CUVE = OPER (nom = "CUVE",
 277              sd_prod = loi,
 278              op = 68,
 279              fr = "Caracteristiques de la cuve",
 280
 281 #===
 282 # Liste des paramètres
 283 #===
 284
 285   TYPEGEOM = SIMP ( statut = "o",
 286                     typ = "TXM",
 287                     into = ( "GEOMETRIE", "MAILLAGE"),
 288                     #defaut = "GEOMETRIE",
 289                     fr = "Traitement de la geometrie d'une cuve",
 290                    ),
 291
 292
 293 #====
 294 # Definition des parametres selon le type de traitement de la geometrie
 295 #====
 296
 297   Geometrie = BLOC ( condition = " TYPEGEOM in ( 'GEOMETRIE', ) ",
 298
 299                   RINT = SIMP ( statut = "o",
 300                                 typ = "R",
 301                                 defaut = "1.994",
 302                                 fr = "Rayon interne de la cuve (m)",
 303                                ),
 304
 305                   REXT = SIMP ( statut = "o",
 306                                 typ = "R",
 307                                 defaut = "2,2015",
 308                                 fr = "Rayon externe de la cuve (m)",
 309                                ),
 310
 311                   LREV = SIMP ( statut = "o",
 312                                 typ = "R",
 313                                 defaut = "0.0075",
 314                                 fr = "Epaisseur du revetement (m)",
 315                                ),
 316
 317                   LIGMIN = SIMP ( statut = "o",
 318                                   typ = "R",
 319                                   defaut = "0.75",
 320                                   fr = "Ligament externe minimal avant rupture (% de l'epaisseur de cuve)",
 321                                  ),
 322
 323                      ), # Fin BLOC  Geometrie
 324
 325   Maillage = BLOC ( condition = " TYPEGEOM in ( 'MAILLAGE', ) ",
 326
 327                   Liste_abscisses = SIMP ( statut = "o",
 328                                            typ = "R",
 329                                            max = "**",
 330                                            fr = "Liste des abscisses (m) A FAIRE",
 331                                           ),
 332                    ), # Fin BLOC Maillage
 333
 334 ) # Fin OPER CUVE
 335
 336 #====================================================
 337 # 4. Modeles de fluence, d'irradiation et de tenacite
 338 #====================================================
 339
 340 #=======================
 341 # 4.1 Modeles de fluence
 342 #=======================
 343
 344 # Nota : les variables de type OPER doivent etre en majuscules !
 345 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
 346 MODELES = OPER ( nom = "MODELES",
 347                  sd_prod = loi,
 348                  op = 68,
 349                  fr = "Modeles de fluence, d'irradiation et de tenacite",
 350
 351
 352 #===
 353 # Liste des paramètres
 354 #===
 355
 356   MODELFLUENCE = SIMP ( statut = "o",
 357                         typ = "TXM",
 358                         into = ( "Reglementaire", "France", "ValeurImposee", "SDM", "USNRC", "REV_2", "SDM_Lissage", "GrandeDev", "GD_Cuve"),
 359                         #defaut = "Reglementaire",
 360                         fr = "Modele d'attenuation de la fluence dans l'epaisseur de la cuve",
 361                        ),
 362
 363
 364 #====
 365 # Definition des parametres selon le modele de fluence
 366 #====
 367
 368   Reglementaire = BLOC ( condition = " MODELFLUENCE in ( 'Reglementaire', ) ",
 369
 370                   fmax = SIMP ( statut = "o",
 371                                 typ = "R",
 372                                 defaut = "6.5",
 373                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 374                                ),
 375
 376                         ), # Fin BLOC Reglementaire
 377
 378   France = BLOC ( condition = " MODELFLUENCE in ( 'France', ) ",
 379
 380                   fmax = SIMP ( statut = "o",
 381                                 typ = "R",
 382                                 defaut = "6.5",
 383                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 384                                ),
 385
 386                   KPFRANCE = SIMP ( statut = "o",
 387                                     typ = "R",
 388                                     defaut = "12.7",
 389                                     fr = "Parametre exponentiel du modele France",
 390                                    ),
 391
 392                    ), # Fin BLOC France
 393
 394   ValeurImposee = BLOC ( condition = " MODELFLUENCE in ( 'ValeurImposee', ) ",
 395
 396                   fmax = SIMP ( statut = "o",
 397                                 typ = "R",
 398                                 defaut = "6.5",
 399                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 400                                ),
 401
 402                          ), # Fin BLOC ValeurImposee
 403
 404   SDM = BLOC ( condition = " MODELFLUENCE in ( 'SDM', ) ",
 405
 406                   fmax = SIMP ( statut = "o",
 407                                 typ = "R",
 408                                 defaut = "6.5",
 409                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 410                                ),
 411
 412               ), # Fin BLOC SDM
 413
 414   USNRC = BLOC ( condition = " MODELFLUENCE in ( 'USNRC', ) ",
 415
 416                   fmax = SIMP ( statut = "o",
 417                                 typ = "R",
 418                                 defaut = "6.5",
 419                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 420                                ),
 421
 422                   KPUS = SIMP ( statut = "o",
 423                                 typ = "R",
 424                                 defaut = "9.4488",
 425                                 fr = "Parametre exponentiel du modele US",
 426                                ),
 427
 428                 ), # Fin BLOC USNRC
 429
 430   REV_2 = BLOC ( condition = " MODELFLUENCE in ( 'REV_2', ) ",
 431
 432                   fmax = SIMP ( statut = "o",
 433                                 typ = "R",
 434                                 defaut = "6.5",
 435                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 436                                ),
 437
 438                 ), # Fin BLOC REV_2
 439
 440   SDM_Lissage = BLOC ( condition = " MODELFLUENCE in ( 'SDM_Lissage', ) ",
 441
 442                   fmax = SIMP ( statut = "o",
 443                                 typ = "R",
 444                                 defaut = "6.5",
 445                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 446                                ),
 447
 448                       ), # Fin BLOC SDM_Lissage
 449
 450   GrandeDev = BLOC ( condition = " MODELFLUENCE in ( 'GrandeDev', ) ",
 451
 452                   fmax = SIMP ( statut = "o",
 453                                 typ = "R",
 454                                 defaut = "6.5",
 455                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 456                                ),
 457
 458                      ), # Fin BLOC GrandeDev
 459
 460   GD_Cuve = BLOC ( condition = " MODELFLUENCE in ( 'GD_Cuve', ) ",
 461
 462                   fmax = SIMP ( statut = "o",
 463                                 typ = "R",
 464                                 defaut = "6.5",
 465                                 fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
 466                                ),
 467
 468                   COEFFLUENCE1 = SIMP ( statut = "o",
 469                                         typ = "R",
 470                                         defaut = "5.8",
 471                                         fr = "Fluence a l'azimut 0 (10^19 n/cm)",
 472                                        ),
 473
 474                   COEFFLUENCE2 = SIMP ( statut = "o",
 475                                         typ = "R",
 476                                         defaut = "5.48",
 477                                         fr = "Fluence a l'azimut 5 (10^19 n/cm)",
 478                                        ),
 479
 480                   COEFFLUENCE3 = SIMP ( statut = "o",
 481                                         typ = "R",
 482                                         defaut = "4.46",
 483                                         fr = "Fluence a l'azimut 10 (10^19 n/cm)",
 484                                        ),
 485
 486                   COEFFLUENCE4 = SIMP ( statut = "o",
 487                                         typ = "R",
 488                                         defaut = "3.41",
 489                                         fr = "Fluence a l'azimut 15 (10^19 n/cm)",
 490                                        ),
 491
 492                   COEFFLUENCE5 = SIMP ( statut = "o",
 493                                         typ = "R",
 494                                         defaut = "3.37",
 495                                         fr = "Fluence a l'azimut 20 (10^19 n/cm)",
 496                                        ),
 497
 498                   COEFFLUENCE6 = SIMP ( statut = "o",
 499                                         typ = "R",
 500                                         defaut = "3.16",
 501                                         fr = "Fluence a l'azimut 25 (10^19 n/cm)",
 502                                        ),
 503
 504                   COEFFLUENCE7 = SIMP ( statut = "o",
 505                                         typ = "R",
 506                                         defaut = "2.74",
 507                                         fr = "Fluence a l'azimut 30 (10^19 n/cm)",
 508                                        ),
 509
 510                   COEFFLUENCE8 = SIMP ( statut = "o",
 511                                         typ = "R",
 512                                         defaut = "2.25",
 513                                         fr = "Fluence a l'azimut 35 (10^19 n/cm)",
 514                                        ),
 515
 516                   COEFFLUENCE9 = SIMP ( statut = "o",
 517                                         typ = "R",
 518                                         defaut = "1.89",
 519                                         fr = "Fluence a l'azimut 40 (10^19 n/cm)",
 520                                        ),
 521
 522                   COEFFLUENCE10 = SIMP ( statut = "o",
 523                                         typ = "R",
 524                                         defaut = "1.78",
 525                                         fr = "Fluence a l'azimut 45 (10^19 n/cm)",
 526                                        ),
 527
 528                   ), # Fin BLOC GD_Cuve
 529
 530 #==========================
 531 # 4.2 Modeles d'irradiation
 532 #==========================
 533
 534   TYPIRR = SIMP ( statut = "o",
 535                   typ = "TXM",
 536                   into = ( "RTNDT", "FLUENCE" ),
 537                   #defaut = "RTNDT",
 538                   fr = "Type d'irradiation ",
 539                  ),
 540
 541 #====
 542 # Definition des parametres selon le type d'irradiation
 543 #====
 544
 545   Parametres_RTNDT = BLOC ( condition = " TYPIRR in ( 'RTNDT', ) ",
 546
 547                   RTNDT = SIMP ( statut = "o",
 548                                  typ = "R",
 549                                  defaut = "73.",
 550                                  fr = "RTNDT finale (°C)",
 551                                 ),
 552
 553                             ), # Fin BLOC Parametres_RTNDT
 554
 555   Parametres_FLUENCE = BLOC ( condition = " TYPIRR in ( 'FLUENCE', ) ",
 556
 557                   MODELIRR = SIMP ( statut = "o",
 558                                     typ = "TXM",
 559                                     into = ( "HOUSSIN", "PERSOZ", "LEFEBVRE", "USNRCmdb", "BRILLAUD", "USNRCsoud" ),
 560                                     #defaut = "HOUSSIN",
 561                                     fr = "Modele d'irradiation pour virole ou joint soude",
 562                                    ),
 563
 564                   CU = SIMP ( statut = "o",
 565                               typ = "R",
 566                               defaut = "0.",
 567                               fr = "Teneur en cuivre (%)",
 568                              ),
 569
 570                   Ni = SIMP ( statut = "o",
 571                               typ = "R",
 572                               defaut = "0.",
 573                               fr = "Teneur en nickel (%)",
 574                              ),
 575
 576                   P = SIMP ( statut = "o",
 577                              typ = "R",
 578                              defaut = "0.",
 579                              fr = "Teneur en phosphore (%)",
 580                             ),
 581
 582                   RTimoy = SIMP ( statut = "o",
 583                                   typ = "R",
 584                                   defaut = "0.",
 585                                   fr = "Moyenne de RTNDT : virole C1 de cuve Chinon : mdb=>-17.°C et js=>42.°C (HT-56/05/038 : p.52)",
 586                                  ),
 587
 588                   RTicov = SIMP ( statut = "o",
 589                                   typ = "R",
 590                                   defaut = "0.",
 591                                   fr = "Coefficient de variation de la RTNDT initiale",
 592                                  ),
 593
 594   Parametres_USNRC = BLOC ( condition = " MODELIRR in ( 'USNRCsoud', 'USNRCmdb' , ) ",
 595
 596                   USectDRT = SIMP ( statut = "o",
 597                                     typ = "R",
 598                                     defaut = "28.",
 599                                     fr = "pour modeles USNRCsoud ou USNRCmdb, ecart-type du decalage de RTNDT (°F) (28. pour js et 17. pour mdb)",
 600                                    ),
 601
 602                            ), # Fin BLOC Parametres_USNRC
 603
 604                   nbectDRTNDT = SIMP ( statut = "o",
 605                                        typ = "R",
 606                                        defaut = "2.",
 607                                        fr = "Nombre d ecart-type par rapport a la moyenne de DRTNDT",
 608                                       ),
 609
 610                            ), # Fin BLOC Parametres_FLUENCE
 611
 612 #========================
 613 # 4.3 Modeles de tenacite
 614 #========================
 615
 616   MODELKIC = SIMP ( statut = "o",
 617                     typ = "TXM",
 618                     into = ( "RCC-M", "RCC-M_pal", "RCC-M_exp", "Houssin_RC", "Wallin", "REME", "ORNL", "Frama", "WEIB3", "WEIB2", "LOGWOLF", "WEIB_GEN" ),
 619                     #defaut = "RCC-M",
 620                     fr = "Modele de tenacite ",
 621                    ),
 622
 623 #====
 624 # Definition des parametres selon le modele de tenacité
 625 #====
 626
 627   Parametres_RCCM = BLOC ( condition = " MODELKIC in ( 'RCC-M', ) ",
 628
 629                   nbectKIc = SIMP ( statut = "o",
 630                                     typ = "R",
 631                                     defaut = "-2.",
 632                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 633                                    ),
 634
 635                   fractKIc = SIMP ( statut = "o",
 636                                     typ = "R",
 637                                     defaut = "5.",
 638                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 639                                    ),
 640
 641                   KICPAL = SIMP ( statut = "o",
 642                                   typ = "R",
 643                                   defaut = "195.",
 644                                   fr = "palier deterministe de K1c quand modele RCC-M  (MPa(m^0.5)) ",
 645                                  ),
 646
 647                   KICCDV = SIMP ( statut = "o",
 648                                   typ = "R",
 649                                   defaut = "0.15",
 650                                   fr = "coef de variation de la loi normale de K1c quand modele RCC-M ",
 651                                  ),
 652
 653                   nbectKIa = SIMP ( statut = "o",
 654                                     typ = "R",
 655                                     defaut = "-2.",
 656                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 657                                    ),
 658
 659
 660                   KIAPAL = SIMP ( statut = "o",
 661                                   typ = "R",
 662                                   defaut = "195.",
 663                                   fr = "palier deterministe de K1a -tenacite a l'arret- quand modele RCC-M  (MPa(m^0.5)) ",
 664                                  ),
 665
 666                   KIACDV = SIMP ( statut = "o",
 667                                   typ = "R",
 668                                   defaut = "0.10",
 669                                   fr = "coef de variation de la loi normale de K1a -tenacite a l'arret- quand modele RCC-M ",
 670                                  ),
 671
 672                           ), # Fin BLOC Parametres_RCC-M
 673
 674   Parametres_RCCMpal = BLOC ( condition = " MODELKIC in ( 'RCC-M_pal', ) ",
 675
 676                   nbectKIc = SIMP ( statut = "o",
 677                                     typ = "R",
 678                                     defaut = "-2.",
 679                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 680                                    ),
 681
 682                   fractKIc = SIMP ( statut = "o",
 683                                     typ = "R",
 684                                     defaut = "5.",
 685                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 686                                    ),
 687
 688                   nbectKIa = SIMP ( statut = "o",
 689                                     typ = "R",
 690                                     defaut = "-2.",
 691                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 692                                    ),
 693
 694                              ), # Fin BLOC Parametres_RCCMpal
 695
 696   Parametres_RCCMexp = BLOC ( condition = " MODELKIC in ( 'RCC-M_exp', ) ",
 697
 698                   nbectKIc = SIMP ( statut = "o",
 699                                     typ = "R",
 700                                     defaut = "-2.",
 701                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 702                                    ),
 703
 704                   fractKIc = SIMP ( statut = "o",
 705                                     typ = "R",
 706                                     defaut = "5.",
 707                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 708                                    ),
 709
 710                   nbectKIa = SIMP ( statut = "o",
 711                                     typ = "R",
 712                                     defaut = "-2.",
 713                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 714                                    ),
 715
 716                              ), # Fin BLOC Parametres_RCCMexp
 717
 718   Parametres_Houssin_RC = BLOC ( condition = " MODELKIC in ( 'Houssin_RC', ) ",
 719
 720                   nbectKIc = SIMP ( statut = "o",
 721                                     typ = "R",
 722                                     defaut = "-2.",
 723                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 724                                    ),
 725
 726                   fractKIc = SIMP ( statut = "o",
 727                                     typ = "R",
 728                                     defaut = "5.",
 729                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 730                                    ),
 731
 732                   nbectKIa = SIMP ( statut = "o",
 733                                     typ = "R",
 734                                     defaut = "-2.",
 735                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 736                                    ),
 737
 738                                 ), # Fin BLOC Parametres_Houssin_RC
 739
 740   Parametres_Wallin = BLOC ( condition = " MODELKIC in ( 'Wallin', ) ",
 741
 742                   T0WALLIN = SIMP ( statut = "o",
 743                                     typ = "I",
 744                                     defaut = "-27",
 745                                     fr = "parametre T0 du modele Wallin (°C)",
 746                                    ),
 747
 748                   nbectKIc = SIMP ( statut = "o",
 749                                     typ = "R",
 750                                     defaut = "-2.",
 751                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 752                                    ),
 753
 754                   fractKIc = SIMP ( statut = "o",
 755                                     typ = "R",
 756                                     defaut = "5.",
 757                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 758                                    ),
 759
 760                   nbectKIa = SIMP ( statut = "o",
 761                                     typ = "R",
 762                                     defaut = "-2.",
 763                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 764                                    ),
 765
 766                             ), # Fin BLOC Parametres_Wallin
 767
 768   Parametres_REME = BLOC ( condition = " MODELKIC in ( 'REME', ) ",
 769
 770                   nbectKIc = SIMP ( statut = "o",
 771                                     typ = "R",
 772                                     defaut = "-2.",
 773                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 774                                    ),
 775
 776                   fractKIc = SIMP ( statut = "o",
 777                                     typ = "R",
 778                                     defaut = "5.",
 779                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 780                                    ),
 781
 782                   nbectKIa = SIMP ( statut = "o",
 783                                     typ = "R",
 784                                     defaut = "-2.",
 785                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 786                                    ),
 787
 788                           ), # Fin BLOC Parametres_REME
 789
 790   Parametres_ORNL = BLOC ( condition = " MODELKIC in ( 'ORNL', ) ",
 791
 792                   nbectKIc = SIMP ( statut = "o",
 793                                     typ = "R",
 794                                     defaut = "-2.",
 795                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 796                                    ),
 797
 798                   fractKIc = SIMP ( statut = "o",
 799                                     typ = "R",
 800                                     defaut = "5.",
 801                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 802                                    ),
 803
 804                   nbectKIa = SIMP ( statut = "o",
 805                                     typ = "R",
 806                                     defaut = "-2.",
 807                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 808                                    ),
 809
 810                           ), # Fin BLOC Parametres_ORNL
 811
 812   Parametres_Frama = BLOC ( condition = " MODELKIC in ( 'Frama', ) ",
 813
 814                   nbectKIc = SIMP ( statut = "o",
 815                                     typ = "R",
 816                                     defaut = "-2.",
 817                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 818                                    ),
 819
 820                   fractKIc = SIMP ( statut = "o",
 821                                     typ = "R",
 822                                     defaut = "5.",
 823                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 824                                    ),
 825
 826                   nbectKIa = SIMP ( statut = "o",
 827                                     typ = "R",
 828                                     defaut = "-2.",
 829                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 830                                    ),
 831
 832                            ), # Fin BLOC Parametres_Frama
 833
 834   Parametres_WEIB3 = BLOC ( condition = " MODELKIC in ( 'WEIB3', ) ",
 835
 836                   nbectKIc = SIMP ( statut = "o",
 837                                     typ = "R",
 838                                     defaut = "-2.",
 839                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 840                                    ),
 841
 842                   fractKIc = SIMP ( statut = "o",
 843                                     typ = "R",
 844                                     defaut = "5.",
 845                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 846                                    ),
 847
 848                   nbectKIa = SIMP ( statut = "o",
 849                                     typ = "R",
 850                                     defaut = "-2.",
 851                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 852                                    ),
 853
 854                            ), # Fin BLOC Parametres_WEIB3
 855
 856   Parametres_WEIB2 = BLOC ( condition = " MODELKIC in ( 'WEIB2', ) ",
 857
 858                   nbectKIc = SIMP ( statut = "o",
 859                                     typ = "R",
 860                                     defaut = "-2.",
 861                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 862                                    ),
 863
 864                   fractKIc = SIMP ( statut = "o",
 865                                     typ = "R",
 866                                     defaut = "5.",
 867                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 868                                    ),
 869
 870                   nbectKIa = SIMP ( statut = "o",
 871                                     typ = "R",
 872                                     defaut = "-2.",
 873                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 874                                    ),
 875
 876                            ), # Fin BLOC Parametres_WEIB2
 877
 878   Parametres_LOGWOLF = BLOC ( condition = " MODELKIC in ( 'LOGWOLF', ) ",
 879
 880                   nbectKIc = SIMP ( statut = "o",
 881                                     typ = "R",
 882                                     defaut = "-2.",
 883                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 884                                    ),
 885
 886                   fractKIc = SIMP ( statut = "o",
 887                                     typ = "R",
 888                                     defaut = "5.",
 889                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 890                                    ),
 891
 892                   nbectKIa = SIMP ( statut = "o",
 893                                     typ = "R",
 894                                     defaut = "-2.",
 895                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 896                                    ),
 897
 898                              ), # Fin BLOC Parametres_LOGWOLF
 899
 900   Parametres_WEIB_GEN = BLOC ( condition = " MODELKIC in ( 'WEIB_GEN',) ",
 901
 902                   A1 = SIMP ( statut = "o",
 903                               typ = "R",
 904                               defaut = "21.263",
 905                               fr = "coef du coef a(T) d'une Weibull generale",
 906                              ),
 907
 908                   A2 = SIMP ( statut = "o",
 909                               typ = "R",
 910                               defaut = "9.159",
 911                               fr = "coef du coef a(T) d'une Weibull generale",
 912                              ),
 913
 914                   A3 = SIMP ( statut = "o",
 915                               typ = "R",
 916                               defaut = "0.04057",
 917                               fr = "coef du coef a(T) d'une Weibull generale",
 918                              ),
 919
 920                   B1 = SIMP ( statut = "o",
 921                               typ = "R",
 922                               defaut = "17.153",
 923                               fr = "coef du coef b(T) d'une Weibull generale",
 924                              ),
 925
 926                   B2 = SIMP ( statut = "o",
 927                               typ = "R",
 928                               defaut = "55.089",
 929                               fr = "coef du coef b(T) d'une Weibull generale",
 930                              ),
 931
 932                   B3 = SIMP ( statut = "o",
 933                               typ = "R",
 934                               defaut = "0.0144",
 935                               fr = "coef du coef b(T) d'une Weibull generale",
 936                              ),
 937
 938                   C1 = SIMP ( statut = "o",
 939                               typ = "R",
 940                               defaut = "4.",
 941                               fr = "coef du coef c(T) d'une Weibull generale",
 942                              ),
 943
 944                   C2 = SIMP ( statut = "o",
 945                               typ = "R",
 946                               defaut = "0.",
 947                               fr = "coef du coef c(T) d'une Weibull generale",
 948                              ),
 949
 950                   C3 = SIMP ( statut = "o",
 951                               typ = "R",
 952                               defaut = "0.",
 953                               fr = "coef du coef c(T) d'une Weibull generale",
 954                              ),
 955
 956                   nbectKIc = SIMP ( statut = "o",
 957                                     typ = "R",
 958                                     defaut = "-2.",
 959                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
 960                                    ),
 961
 962                   fractKIc = SIMP ( statut = "o",
 963                                     typ = "R",
 964                                     defaut = "5.",
 965                                     fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
 966                                    ),
 967
 968                   nbectKIa = SIMP ( statut = "o",
 969                                     typ = "R",
 970                                     defaut = "-2.",
 971                                     fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
 972                                    ),
 973
 974                               ), # Fin BLOC Parametres_WEIB_GEN
 975
 976 ) # Fin OPER MODELES
 977
 978
 979 #==================
 980 # 5. Initialisation
 981 #==================
 982
 983 # Nota : les variables de type OPER doivent etre en majuscules !
 984 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
 985 INITIALISATION = OPER ( nom = "INITIALISATION",
 986                         sd_prod = loi,
 987                         op = 68,
 988                         fr = "Initialisation",
 989
 990   Liste_tempinit = SIMP ( statut = "o",
 991                           typ = "R",
 992                           max = "**",
 993                           fr = "Profil radial de la temperature initiale dans la cuve (m) (°C) ",
 994                          ),
 995
 996   Prolongation_tempsinit = SIMP ( statut = "o",
 997                                   typ = "TXM",
 998                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
 999                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1000                                  ),
1001
1002   Liste_sigmainit = SIMP ( statut = "o",
1003                            typ = "R",
1004                            max = "**",
1005                            fr = "Profil radial des contraintes residuelles dans la cuve (m) (xx) (xx) (xx) ",
1006                           ),
1007
1008   Prolongation_sigmainit = SIMP ( statut = "o",
1009                                   typ = "TXM",
1010                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1011                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1012                                  ),
1013
1014   INSTINIT = SIMP ( statut = "o",
1015                     typ = "R",
1016                     defaut = "-1.",
1017                     fr = "Instant initial (s) ",
1018                    ),
1019
1020 ) # Fin OPER INITIALISATION
1021
1022
1023 #==================================
1024 # 6. CARACTERISTIQUES DU REVETEMENT
1025 #==================================
1026
1027 # Nota : les variables de type OPER doivent etre en majuscules !
1028 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
1029 REVETEMENT = OPER ( nom = "REVETEMENT",
1030                     sd_prod = loi,
1031                     op = 68,
1032                     fr = "Caracteristiques du revetement",
1033
1034   KTHREV = SIMP ( statut = "o",
1035                   typ = "TXM",
1036                   into = ( "ENTHALPIE", "CHALEUR",),
1037                   #defaut = "CHALEUR",
1038                   fr = "Options pour definir les caracteristiques du revetement ",
1039                  ),
1040
1041   Parametres_ENTHALPIErev = BLOC ( condition = " KTHREV in ( 'ENTHALPIE', ) ",
1042
1043                   Liste_enthalpie_rev = SIMP ( statut = "o",
1044                                                typ = "R",
1045                                                max = "**",
1046                                                fr = "Temperature (°C) / enthalpie (J/kg) ",
1047                                               ),
1048
1049                   Prolongation_enthalpie_rev = SIMP ( statut = "o",
1050                                                       typ = "TXM",
1051                                                       into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1052                                                       fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1053                                                      ),
1054
1055                                   ), # Fin BLOC Parametres_ENTHALPIErev
1056
1057
1058   Parametres_CHALEURrev = BLOC ( condition = " KTHREV in ( 'CHALEUR', ) ",
1059
1060                   Liste_chaleur_rev = SIMP ( statut = "o",
1061                                              typ = "R",
1062                                              max = "**",
1063                                              fr = "Temperature (°C) / chaleur volumique (J/kg/K) ",
1064                                             ),
1065
1066                   Prolongation_chaleur_rev = SIMP ( statut = "o",
1067                                                     typ = "TXM",
1068                                                     into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1069                                                     fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1070                                                    ),
1071
1072                                 ), # Fin BLOC Parametres_CHALEURrev
1073
1074    Liste_conductivite_rev = SIMP ( statut = "o",
1075                                    typ = "R",
1076                                    max = "**",
1077                                    fr = "Temperature (°C) / conductivite thermique (W/m/°C) ",
1078                                   ),
1079
1080    Prolongation_conductivite_rev = SIMP ( statut = "o",
1081                                           typ = "TXM",
1082                                           into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1083                                           fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1084                                          ),
1085
1086    Liste_young_rev = SIMP ( statut = "o",
1087                             typ = "R",
1088                             max = "**",
1089                             fr = "Temperature (°C) / module d'Young (MPa) ",
1090                            ),
1091
1092    Prolongation_young_rev = SIMP ( statut = "o",
1093                                    typ = "TXM",
1094                                    into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1095                                    fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1096                                   ),
1097
1098    Liste_dilatation_rev = SIMP ( statut = "o",
1099                                  typ = "R",
1100                                  max = "**",
1101                                  fr = "Temperature (°C) / coefficient de dilatation thermique (°C-1) ",
1102                                 ),
1103
1104    Prolongation_dilatation_rev = SIMP ( statut = "o",
1105                                         typ = "TXM",
1106                                         into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1107                                         fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1108                                        ),
1109
1110    Liste_elasticite_rev = SIMP ( statut = "o",
1111                                  typ = "R",
1112                                  max = "**",
1113                                  fr = "Temperature (°C) / limite d'elasticite (MPa) ",
1114                                 ),
1115
1116    Prolongation_elasticite_rev = SIMP ( statut = "o",
1117                                         typ = "TXM",
1118                                         into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1119                                         fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1120                                        ),
1121
1122    TREFREV = SIMP ( statut = "o",
1123                     typ = "R",
1124                     defaut = "20.",
1125                     fr = "Temperature de deformation nulle (°C) ",
1126                    ),
1127
1128    TDETREV = SIMP ( statut = "o",
1129                     typ = "R",
1130                     defaut = "287.",
1131                     fr = "Temperature de definition du coefficient de dilatation thermique (°C) ",
1132                    ),
1133
1134    NUREV = SIMP ( statut = "o",
1135                   typ = "R",
1136                   defaut = "0.3",
1137                   fr = "Coefficient de Poisson ",
1138                  ),
1139
1140 ) # Fin OPER REVETEMENT
1141
1142
1143 #=====================================
1144 # 7. CARACTERISTIQUES DU METAL DE BASE
1145 #=====================================
1146
1147 # Nota : les variables de type OPER doivent etre en majuscules !
1148 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
1149 METAL_BASE = OPER ( nom = "METAL_BASE",
1150                     sd_prod = loi,
1151                     op = 68,
1152                     fr = "Caracteristiques du metal de base",
1153
1154   KTHMDB = SIMP ( statut = "o",
1155                   typ = "TXM",
1156                   into = ( "ENTHALPIE", "CHALEUR",),
1157                   #defaut = "CHALEUR",
1158                   fr = "Options pour definir les caracteristiques du revetement ",
1159                  ),
1160
1161   Parametres_ENTHALPIEmdb = BLOC ( condition = " KTHMDB in ( 'ENTHALPIE', ) ",
1162
1163                   Liste_enthalpie_mdb = SIMP ( statut = "o",
1164                                                typ = "R",
1165                                                max = "**",
1166                                                fr = "Temperature (°C) / enthalpie (J/kg) ",
1167                                               ),
1168
1169                   Prolongation_enthalpie_mdb = SIMP ( statut = "o",
1170                                                       typ = "TXM",
1171                                                       into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1172                                                       fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1173                                                      ),
1174
1175                                   ), # Fin BLOC Parametres_ENTHALPIEmdb
1176
1177
1178   Parametres_CHALEURmdb = BLOC ( condition = " KTHMDB in ( 'CHALEUR', ) ",
1179
1180                   Liste_chaleur_mdb = SIMP ( statut = "o",
1181                                              typ = "R",
1182                                              max = "**",
1183                                              fr = "Temperature (°C) / chaleur volumique (J/kg/K) ",
1184                                             ),
1185
1186                   Prolongation_chaleur_mdb = SIMP ( statut = "o",
1187                                                     typ = "TXM",
1188                                                     into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1189                                                     fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1190                                                    ),
1191
1192                                 ), # Fin BLOC Parametres_CHALEURmdb
1193
1194    Liste_conductivite_mdb = SIMP ( statut = "o",
1195                                    typ = "R",
1196                                    max = "**",
1197                                    fr = "Temperature (°C) / conductivite thermique (W/m/°C) ",
1198                                   ),
1199
1200    Prolongation_conductivite_mdb = SIMP ( statut = "o",
1201                                           typ = "TXM",
1202                                           into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1203                                           fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1204                                          ),
1205
1206    Liste_young_mdb = SIMP ( statut = "o",
1207                             typ = "R",
1208                             max = "**",
1209                             fr = "Temperature (°C) / module d'Young (MPa) ",
1210                            ),
1211
1212    Prolongation_young_mdb = SIMP ( statut = "o",
1213                                    typ = "TXM",
1214                                    into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1215                                    fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1216                                   ),
1217
1218    Liste_dilatation_mdb = SIMP ( statut = "o",
1219                                  typ = "R",
1220                                  max = "**",
1221                                  fr = "Temperature (°C) / coefficient de dilatation thermique (°C-1) ",
1222                                 ),
1223
1224    Prolongation_dilatation_mdb = SIMP ( statut = "o",
1225                                         typ = "TXM",
1226                                         into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1227                                         fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1228                                        ),
1229
1230    TREFMDB = SIMP ( statut = "o",
1231                     typ = "R",
1232                     defaut = "20.",
1233                     fr = "Temperature de deformation nulle (°C) ",
1234                    ),
1235
1236    TDETMDB = SIMP ( statut = "o",
1237                     typ = "R",
1238                     defaut = "287.",
1239                     fr = "Temperature de definition du coefficient de dilatation thermique (°C) ",
1240                    ),
1241
1242    NUMDB = SIMP ( statut = "o",
1243                   typ = "R",
1244                   defaut = "0.3",
1245                   fr = "Coefficient de Poisson ",
1246                  ),
1247
1248 ) # Fin OPER METAL_BASE
1249
1250
1251 #================
1252 # 8. TRANSITOIRES
1253 #================
1254
1255 # Nota : les variables de type OPER doivent etre en majuscules !
1256 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
1257 TRANSITOIRE = OPER ( nom = "TRANSITOIRE",
1258                      sd_prod = loi,
1259                      op = 68,
1260                      fr = "Description du transitoire",
1261
1262   Liste_pression = SIMP ( statut = "o",
1263                           typ = "R",
1264                           max = "**",
1265                           fr = "instant (s) / pression (MPa) ",
1266                          ),
1267
1268   Prolongation_pression = SIMP ( statut = "o",
1269                                  typ = "TXM",
1270                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1271                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1272                                 ),
1273
1274   TYPCLTH = SIMP ( statut = "o",
1275                    typ = "TXM",
1276                    into = ( "TEMP_IMPO", "FLUX_REP", "ECHANGE", "DEBIT", "TEMP_FLU"),
1277                    #defaut = "ECHANGE",
1278                    fr = "Type de condition thermique en paroi interne ",
1279                   ),
1280
1281   Parametres_TEMP_IMPO = BLOC ( condition = " TYPCLTH in ( 'TEMP_IMPO', ) ",
1282
1283                   Liste_temp_impo = SIMP ( statut = "o",
1284                                            typ = "R",
1285                                            max = "**",
1286                                            fr = "Instant (s) / Temperature imposee (°C) ",
1287                                           ),
1288
1289                   Prolongation_temp_impo = SIMP ( statut = "o",
1290                                                   typ = "TXM",
1291                                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1292                                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1293                                                  ),
1294
1295                                ), # Fin BLOC Parametres_TEMP_IMPO
1296
1297   Parametres_FLUX_REP = BLOC ( condition = " TYPCLTH in ( 'FLUX_REP', ) ",
1298
1299                   Liste_flux_rep = SIMP ( statut = "o",
1300                                           typ = "R",
1301                                           max = "**",
1302                                           fr = "Instant (s) / Flux de chaleur impose (W/m2) ",
1303                                          ),
1304
1305                   Prolongation_flux_rep = SIMP ( statut = "o",
1306                                                  typ = "TXM",
1307                                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1308                                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1309                                                 ),
1310
1311                               ), # Fin BLOC Parametres_FLUX_REP
1312
1313   Parametres_ECHANGE = BLOC ( condition = " TYPCLTH in ( 'ECHANGE', ) ",
1314
1315                   Liste_echange_temp = SIMP ( statut = "o",
1316                                               typ = "R",
1317                                               max = "**",
1318                                               fr = "Instant (s) / Temperature impose (°C) ",
1319                                              ),
1320
1321                   Prolongation_echange_temp = SIMP ( statut = "o",
1322                                                      typ = "TXM",
1323                                                      into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1324                                                      fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1325                                                     ),
1326
1327                   Liste_echange_coef = SIMP ( statut = "o",
1328                                               typ = "R",
1329                                               max = "**",
1330                                               fr = "Instant (s) / Coefficient d echange (W/m2/K) ",
1331                                              ),
1332
1333                   Prolongation_echange_coef = SIMP ( statut = "o",
1334                                                      typ = "TXM",
1335                                                      into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1336                                                      fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1337                                                     ),
1338
1339                              ), # Fin BLOC Parametres_ECHANGE
1340
1341   Parametres_DEBIT = BLOC ( condition = " TYPCLTH in ( 'DEBIT', ) ",
1342
1343                   Liste_debit = SIMP ( statut = "o",
1344                                        typ = "R",
1345                                        max = "**",
1346                                        fr = "Instant (s) / Debit massique (kg/s) ",
1347                                       ),
1348
1349                   Prolongation_debit = SIMP ( statut = "o",
1350                                               typ = "TXM",
1351                                               into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1352                                               fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1353                                              ),
1354
1355                   Liste_temp_injection = SIMP ( statut = "o",
1356                                                 typ = "R",
1357                                                 max = "**",
1358                                                 fr = "Instant (s) / Temperature d injection de securite  (°C) ",
1359                                                ),
1360
1361                   Prolongation_temp_injection = SIMP ( statut = "o",
1362                                                        typ = "TXM",
1363                                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1364                                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1365                                                       ),
1366
1367                   DH = SIMP ( statut = "o",
1368                               typ = "R",
1369                               defaut = "-2.",
1370                               fr = "Diametre hydraulique (m) ",
1371                              ),
1372
1373                   SECTION = SIMP ( statut = "o",
1374                                    typ = "R",
1375                                    defaut = "-2.",
1376                                    fr = "Section espace annulaire (m2) ",
1377                                   ),
1378
1379                   DELTA = SIMP ( statut = "o",
1380                                  typ = "R",
1381                                  defaut = "-2.",
1382                                  fr = "Hauteur caracteristique convection naturelle (m) ",
1383                                 ),
1384
1385                   ALPHA_CF = SIMP ( statut = "o",
1386                                     typ = "R",
1387                                     defaut = "1.",
1388                                     fr = "Coefficient Vestale convection forcee (-) ",
1389                                    ),
1390
1391                   ALPHA_CN = SIMP ( statut = "o",
1392                                     typ = "R",
1393                                     defaut = "1.",
1394                                     fr = "Coefficient Vestale convection naturelle (-) ",
1395                                    ),
1396
1397                   EPS = SIMP ( statut = "o",
1398                                typ = "R",
1399                                defaut = "0.00001",
1400                                fr = "Critere convergence relative (-) ",
1401                               ),
1402
1403                   VM = SIMP ( statut = "o",
1404                               typ = "R",
1405                               defaut = "-2.",
1406                               fr = "Volume de melange CREARE (m3) ",
1407                              ),
1408
1409                   T0 = SIMP ( statut = "o",
1410                               typ = "R",
1411                               defaut = "-2.",
1412                               fr = "Temperature initiale CREARE (degC) ",
1413                              ),
1414
1415                   SE = SIMP ( statut = "o",
1416                               typ = "R",
1417                               defaut = "-2.",
1418                               fr = "Surface d'echange fluide/structure (m2) ",
1419                              ),
1420
1421                            ), # Fin BLOC Parametres_DEBIT
1422
1423   Parametres_TEMP_FLU = BLOC ( condition = " TYPCLTH in ( 'TEMP_FLU', ) ",
1424
1425                   Liste_temp_flu = SIMP ( statut = "o",
1426                                           typ = "R",
1427                                           max = "**",
1428                                           fr = "Instant (s) / Debit massique (kg/s) ",
1429                                          ),
1430
1431                   Prolongation_temp_flu = SIMP ( statut = "o",
1432                                                  typ = "TXM",
1433                                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1434                                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1435                                                 ),
1436
1437                   Liste_debit_injection = SIMP ( statut = "o",
1438                                                  typ = "R",
1439                                                  max = "**",
1440                                                  fr = "Instant (s) / Debit d injection de securite (kg/s) ",
1441                                                 ),
1442
1443                   Prolongation_debit_injection = SIMP ( statut = "o",
1444                                                         typ = "TXM",
1445                                                         into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
1446                                                         fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
1447                                                        ),
1448
1449                   DH = SIMP ( statut = "o",
1450                               typ = "R",
1451                               defaut = "-2.",
1452                               fr = "Diametre hydraulique (m) ",
1453                              ),
1454
1455                   SECTION = SIMP ( statut = "o",
1456                                    typ = "R",
1457                                    defaut = "-2.",
1458                                    fr = "Section espace annulaire (m2) ",
1459                                   ),
1460
1461                   DELTA = SIMP ( statut = "o",
1462                                  typ = "R",
1463                                  defaut = "-2.",
1464                                  fr = "Hauteur caracteristique convection naturelle (m) ",
1465                                 ),
1466
1467                   ALPHA_CF = SIMP ( statut = "o",
1468                                     typ = "R",
1469                                     defaut = "1.",
1470                                     fr = "Coefficient Vestale convection forcee (-) ",
1471                                    ),
1472
1473                   ALPHA_CN = SIMP ( statut = "o",
1474                                     typ = "R",
1475                                     defaut = "1.",
1476                                     fr = "Coefficient Vestale convection naturelle (-) ",
1477                                    ),
1478
1479                   EPS = SIMP ( statut = "o",
1480                                typ = "R",
1481                                defaut = "0.00001",
1482                                fr = "Critere convergence relative (-) ",
1483                               ),
1484
1485                               ), # Fin BLOC Parametres_TEMP_FLU
1486
1487 ) # Fin OPER TRANSITOIRE
1488
1489
1490 #================================
1491 # 3. Definition de l'etude
1492 #================================
1493
1494 # Nota : les variables de type PROC doivent etre en majuscules !
1495 CRITERIA = PROC ( nom = "CRITERIA",
1496                   op = None,
1497                   docu = "",
1498                   fr = "Mise en donnee pour le fichier de configuration de OPENTURNS.",
1499                   ang = "Writes the configuration file for OPENTURNS.",
1500
1501
1502
1503   Type = SIMP ( statut = "o",
1504                 typ = "TXM",
1505                 into = ( "Min/Max", "Central Uncertainty", "Threshold Exceedence" ),
1506                 fr = "Type d'Analyse",
1507                 ang = "Analysis",
1508                 ),
1509
1510
1511
1512
1513   MinMax = BLOC ( condition = " Type in ( 'Min/Max', ) ",
1514
1515                   Method = SIMP ( statut = "o",
1516                                   typ = "TXM",
1517                                   into = ( "Experiment Plane", "Random Sampling" ),
1518                                   fr = "Methode",
1519                                   ang = "Method",
1520                                   ),
1521                   # UC 3.1.1
1522                   ExperimentPlaneSettings = BLOC ( condition = " Method in ( 'Experiment Plane', ) ",
1523
1524                           ExperimentPlane = SIMP ( statut = "o",
1525                                                    typ = "TXM",
1526                                                    into = ( "Axial", "Factorial", "Composite", ),
1527                                                    fr = "Methode",
1528                                                    ang = "Method",
1529                                                    ),
1530
1531                           Levels = SIMP ( statut = "o",
1532                                           typ = "R",
1533                                           val_min = 0.0,
1534                                           max = '**',
1535                                           fr = "Nombre de niveaux dans chaque direction",
1536                                           ang = "Levels in each direction",
1537                                           ),
1538
1539                           # Scaled Vector
1540                           UnitsPerDimension = SIMP ( statut = "o",
1541                                           typ = "R",
1542                                           max = '**',
1543                                           fr = "Unite par dimension (autant que de variables declarees)",
1544                                           ang = "Units per dimension (as much as declared variables)",
1545                                           ),
1546
1547                           # Translation Vector
1548                           Center = SIMP ( statut = "o",
1549                                           typ = "R",
1550                                           max = '**',
1551                                           fr = "Unite par dimension",
1552                                           ang = "Units per dimension",
1553                                           ),
1554
1555                     ), # Fin BLOC ExperimentPlaneSettings
1556
1557
1558
1559                   RandomSamplingSettings = BLOC ( condition = " Method in ( 'Random Sampling', ) ",
1560
1561                           PointsNumber = SIMP ( statut = "o",
1562                                                 typ = "I",
1563                                                 val_min = 1,
1564                                                 fr = "Nombre de points",
1565                                                 ang = "Points number",
1566                                                 ),
1567
1568                     ), # Fin BLOC RandomSamplingSettings
1569
1570                   Result = SIMP (  statut = "o",
1571                                    typ = "TXM",
1572                                    into = ( "Min/Max", ),
1573                                    defaut = "Min/Max",
1574                                    fr = "Le minimum et le maximum",
1575                                    ang = "The min and max values",
1576                                    ),
1577
1578
1579   ), # Fin BLOC MinMax
1580
1581
1582
1583
1584   CentralUncertainty = BLOC ( condition = " Type in ( 'Central Uncertainty', ) ",
1585
1586                   Method = SIMP ( statut = "o",
1587                                   typ = "TXM",
1588                                   into = ( "Taylor Variance Decomposition", "Random Sampling" ),
1589                                   fr = "Methode",
1590                                   ang = "Method",
1591                                   ),
1592
1593                   # UC 3.2.
1594                   TaylorVarianceDecompositionSettings = BLOC ( condition = " Method in ( 'Taylor Variance Decomposition', ) ",
1595
1596                       Result = FACT ( statut = "o",
1597                                       min = 1,
1598                                       max = "**",
1599
1600                               MeanFirstOrder = SIMP ( statut = "o",
1601                                                 typ = 'TXM',
1602                                                 into = ( 'yes', 'no' ),
1603                                                 defaut = 'yes',
1604                                                 max = 1,
1605                                                 fr = "Moyenne au premier ordre",
1606                                                 ang = "MeanFirstOrder",
1607                                                 ),
1608
1609                               StandardDeviationFirstOrder = SIMP ( statut = "o",
1610                                                                    typ = 'TXM',
1611                                                                    into = ( 'yes', 'no' ),
1612                                                                    defaut = 'yes',
1613                                                                    max = 1,
1614                                                                    fr = "Ecart-type au premier ordre",
1615                                                                    ang = "StandardDeviationFirstOrder",
1616                                                                    ),
1617
1618                               MeanSecondOrder = SIMP ( statut = "o",
1619                                                        typ = 'TXM',
1620                                                        into = ( 'yes', 'no' ),
1621                                                        defaut = 'yes',
1622                                                        max = 1,
1623                                                        fr = "Moyenne au second ordre",
1624                                                        ang = "MeanSecondOrder",
1625                                                        ),
1626
1627                               ImportanceFactor = SIMP ( statut = "o",
1628                                                         typ = 'TXM',
1629                                                         into = ( 'yes', 'no' ),
1630                                                         defaut = 'no',
1631                                                         max = 1,
1632                                                         fr = "Facteur d'importance pour variable de sortie scalaire",
1633                                                         ang = "ImportanceFactor",
1634                                                         ),
1635
1636                              ImportanceFactorSettings = BLOC ( condition = " ImportanceFactor in ( 'yes', ) ",
1637
1638                                     NumericalResults  = SIMP ( statut = "o",
1639                                                                typ = 'TXM',
1640                                                                into = ( 'yes', 'no' ),
1641                                                                defaut = 'yes',
1642                                                                max = 1,
1643                                                                fr = "Resultats numeriques",
1644                                                                ang = "NumericalResults",
1645                                                                ),
1646
1647                                      GraphicalResults  = SIMP ( statut = "o",
1648                                                                 typ = 'TXM',
1649                                                                 into = ( 'yes', 'no' ),
1650                                                                 defaut = 'no',
1651                                                                 max = 1,
1652                                                                 fr = "Resultats graphiques",
1653                                                                 ang = "GraphicalResults",
1654                                                                 ),
1655
1656                             ), # Fin BLOC ImportanceFactorSettings
1657
1658                       ), # Fin FACT Result
1659
1660                   ), # Fin BLOC TaylorVarianceDecompositionSettings
1661
1662
1663
1664                   RandomSamplingSettings = BLOC ( condition = " Method in ( 'Random Sampling', ) ",
1665
1666                           PointsNumber = SIMP ( statut = "o",
1667                                                 typ = "I",
1668                                                 val_min = 1,
1669                                                 fr = "Nombre de points",
1670                                                 ang = "Points number",
1671                                                 ),
1672
1673                        Result = FACT ( statut = "o",
1674                                       min = 1,
1675                                       max = "**",
1676
1677                               EmpiricalMean = SIMP ( statut = "o",
1678                                                      typ = 'TXM',
1679                                                      into = ( 'yes', 'no' ),
1680                                                      defaut = 'yes',
1681                                                      max = 1,
1682                                                      fr = "Moyenne empirique",
1683                                                      ang = "Empirical mean",
1684                                                      ),
1685
1686                               EmpiricalStandardDeviation = SIMP ( statut = "o",
1687                                                                   typ = 'TXM',
1688                                                                   into = ( 'yes', 'no' ),
1689                                                                   defaut = 'yes',
1690                                                                   max = 1,
1691                                                                   fr = "Ecart-type empirique",
1692                                                                   ang = "Empirical standard deviation",
1693                                                                   ),
1694
1695                               EmpiricalQuantile = SIMP ( statut = "o",
1696                                                          typ = 'R',
1697                                                          defaut = 0.0,
1698                                                          max = 1,
1699                                                          val_min = 0.0,
1700                                                          val_max = 1.0,
1701                                                          fr = "Quantile empirique",
1702                                                          ang = "Empirical quantile",
1703                                                          ),
1704
1705                               AnalysedCorrelations = SIMP ( statut = "o",
1706                                                             typ = 'TXM',
1707                                                             into = ( 'yes', 'no' ),
1708                                                             defaut = 'no',
1709                                                             max = 1,
1710                                                             fr = "Correlations analysees",
1711                                                             ang = "Analysed correlations",
1712                                                             ),
1713
1714                               KernelSmoothing = SIMP ( statut = "o",
1715                                                        typ = 'TXM',
1716                                                        into = ( 'yes', 'no' ),
1717                                                        defaut = 'no',
1718                                                        max = 1,
1719                                                        fr = "Kernel smoothing de l'echantillon",
1720                                                        ang = "Kernel smoothing of the sample",
1721                                                        ),
1722
1723                       ), # Fin FACT Result
1724
1725                   ), # Fin BLOC RandomSamplingSettings
1726
1727   ), # Fin BLOC CentralUncertainty
1728
1729
1730
1731
1732   ThresholdExceedence = BLOC ( condition = " Type in ( 'Threshold Exceedence', ) ",
1733
1734          Event =  FACT ( statut = "o",
1735                          min = 1,
1736                          max = 1,
1737
1738                          Threshold = SIMP ( statut = "o",
1739                                             typ = "R",
1740                                             max = 1,
1741                                             fr = "Le seuil de defaillance",
1742                                             ang = "Failure threshold",
1743                                             ),
1744
1745                          ComparisonOperator = SIMP ( statut = "o",
1746                                                      typ = "TXM",
1747                                                      max = 1,
1748                                                      into = ( "Less", "LessOrEqual", "Equal", "GreaterOrEqual", "Greater" ),
1749                                                      fr = "Que faut-il ne pas depasser : un maximum ou un minimum",
1750                                                      ang = "What is the failure threshold : maximum or minimum",
1751                                                      ),
1752          ), # Fin FACT Event
1753
1754
1755          Method = SIMP ( statut = "o",
1756                          typ = "TXM",
1757                          into = ( "Simulation", "Analytical" ),
1758                          fr = "Methode",
1759                          ang = "Method",
1760                          ),
1761
1762          SimulationSettings = BLOC ( condition = " Method in ( 'Simulation', ) ",
1763
1764                Algorithm = SIMP ( statut = "o",
1765                                   typ = "TXM",
1766                                   into = ( "MonteCarlo", "LHS", "ImportanceSampling" ),
1767                                   fr = "Algorithme de simulation",
1768                                   ang = "Simulation algorithm",
1769                                   ),
1770
1771
1772                RandomGenerator = FACT ( statut = "o",
1773                                         min = 1,
1774                                         max = 1,
1775
1776                            SeedToBeSet = SIMP ( statut = "o",
1777                                                 typ = 'TXM',
1778                                                 into = ( 'yes', 'no' ),
1779                                                 defaut = 'no',
1780                                                 max = 1,
1781                                                 fr = "La racine du generateur aleatoire doit-elle etre positionnee ?",
1782                                                 ang = "Does the random generator seed need to be set ?",
1783                                                 ),
1784
1785                            SeedSettings = BLOC ( condition = " SeedToBeSet in ( 'yes', ) ",
1786
1787                                                  RandomGeneratorSeed = SIMP ( statut = "o",
1788                                                                               typ = "I",
1789                                                                               max = 1,
1790                                                                               fr = "Racine du generateur aleatoire",
1791                                                                               ang = "Random generator seed",
1792                                                                               ),
1793
1794                                                ), # Fin BLOC SeedSettings
1795
1796                ), # Fin FACT RandomGenerator
1797
1798
1799                BlockSize = SIMP ( statut = "f",
1800                                   typ = "I",
1801                                   max = 1,
1802                                   val_min = 1,
1803                                   defaut = 1,
1804                                   fr = "Nombre de calculs realises en bloc",
1805                                   ang = "Number of computations as a block",
1806                                   ),
1807
1808                MaximumOuterSampling = SIMP ( statut = "o",
1809                                              typ = "I",
1810                                              max = 1,
1811                                              val_min = 1,
1812                                              fr = "Maximum d'iterations externes",
1813                                              ang = "Maximum outer Sampling value",
1814                                              ),
1815
1816                MaximumCoefficientOfVariation = SIMP ( statut = "f",
1817                                                       typ = "R",
1818                                                       max = 1,
1819                                                       defaut = 0.1,
1820                                                       val_min = 0.0,
1821                                                       fr = " maximum ...",
1822                                                       ang = "Absolute maximum ...."
1823                                                       ),
1824
1825                ImportanceSamplingSettings = BLOC ( condition = " Algorithm in ( 'ImportanceSampling', ) ",
1826
1827                             MeanVector = SIMP ( statut = "o",
1828                                                 typ = "R",
1829                                                 max = "**",
1830                                                 fr = "Moyenne",
1831                                                 ang = "Mean vector",
1832                                                 ),
1833
1834                             Correlation = SIMP ( statut = "o",
1835                                                  typ = 'TXM',
1836                                                  into = ( 'Independent', 'Linear' ),
1837                                                  defaut = 'Linear',
1838                                                  max = 1,
1839                                                  fr = "Le type de correlation entre les variables",
1840                                                  ang = "Correlation between variables",
1841                                                  ),
1842
1843                ), # Fin BLOC ImportanceSamplingSettings
1844
1845                Result = FACT ( statut = "o",
1846                                       min = 1,
1847                                       max = "**",
1848
1849                     Probability = SIMP ( statut = "o",
1850                                          typ = 'TXM',
1851                                          into = ( 'yes', ),
1852                                          defaut = 'yes',
1853                                          max = 1,
1854                                          fr = "Probabiblite",
1855                                          ang = "Probability",
1856                                          ),
1857
1858                     ConfidenceInterval = SIMP ( statut = "o",
1859                                                 typ = 'TXM',
1860                                                 into = ( 'yes', 'no' ),
1861                                                 defaut = 'yes',
1862                                                 max = 1,
1863                                                 fr = "Ecart-type empirique",
1864                                                 ang = "Empirical standard deviation",
1865                                                 ),
1866
1867                     ConfidenceIntervalSettings = BLOC ( condition = " ConfidenceInterval in ( 'yes', ) ",
1868
1869                           Level = SIMP ( statut = "o",
1870                                          typ = 'R',
1871                                          defaut = 0.0,
1872                                          max = 1,
1873                                          val_min = 0.0,
1874                                          val_max = 1.0,
1875                                          fr = "Niveau de confiance",
1876                                          ang = "Confidence level",
1877                                          ),
1878
1879                     ), # Fin BLOC ConfidenceIntervalSettings
1880
1881                     VariationCoefficient = SIMP ( statut = "o",
1882                                                   typ = 'TXM',
1883                                                   into = ( 'yes', 'no' ),
1884                                                   defaut = 'yes',
1885                                                   max = 1,
1886                                                   fr = "Coefficient de variation",
1887                                                   ang = "VariationCoefficient",
1888                                                   ),
1889
1890                     IterationNumber = SIMP ( statut = "o",
1891                                              typ = 'TXM',
1892                                              into = ( 'yes', 'no' ),
1893                                              defaut = 'yes',
1894                                              max = 1,
1895                                              fr = "Nombre d'iterations",
1896                                              ang = "Iteration number",
1897                                              ),
1898
1899                     ConvergenceGraph = SIMP ( statut = "o",
1900                                              typ = 'TXM',
1901                                              into = ( 'yes', 'no' ),
1902                                              defaut = 'yes',
1903                                              max = 1,
1904                                              fr = "Graphe de convergence",
1905                                              ang = "Convergence graph",
1906                                              ),
1907
1908                ), # Fin FACT Result
1909
1910
1911
1912          ), # Fin BLOC SimulationSettings
1913
1914
1915
1916          AnalyticalSettings = BLOC ( condition = " Method in ( 'Analytical', ) ",
1917
1918                 Approximation = SIMP ( statut = "o",
1919                                        typ = "TXM",
1920                                        into = ( "FORM", "SORM" ),
1921                                        fr = "Approximation",
1922                                        ang = "Approximation",
1923                                        ),
1924
1925                 OptimizationAlgorithm = SIMP ( statut = "o",
1926                                                typ = "TXM",
1927                                                into = ( "Cobyla", "AbdoRackwitz" ),
1928                                                fr = "Methode d'optimisation",
1929                                                ang = "Optimisation method",
1930                                                ),
1931
1932
1933                 PhysicalStartingPoint = SIMP ( statut = "f",
1934                                                typ = "R",
1935                                                max = "**",
1936                                                fr = "Point de demarrage de l'algorithme iteratif",
1937                                                ang = "Initial point for iterative process",
1938                                                ),
1939
1940                 MaximumIterationsNumber = SIMP ( statut = "f",
1941                                                  typ = "I",
1942                                                  max = 1,
1943                                                  val_min = 1,
1944                                                  fr = "Nombre maximum d iterations",
1945                                                  ang = "Maximum number of iterations",
1946                                                  ),
1947
1948                 regles = ( EXCLUS ( "MaximumAbsoluteError", "RelativeAbsoluteError" ),  ),
1949
1950                 MaximumAbsoluteError = SIMP ( statut = "f",
1951                                               typ = "R",
1952                                               max = 1,
1953                                               defaut = 1E-6,
1954                                               val_min = 0.0,
1955                                               fr = "Distance maximum absolue entre 2 iterations successifs",
1956                                               ang = "Absolute maximum distance between 2 successive iterates",
1957                                               ),
1958
1959                 RelativeAbsoluteError = SIMP ( statut = "f",
1960                                                typ = "R",
1961                                                max = 1,
1962                                                defaut = 1E-6,
1963                                                val_min = 0.0,
1964                                                fr = "Distance maximum relative entre 2 iterations successives",
1965                                                ang = "Relative maximum distance between 2 successive iterates",
1966                                                ),
1967
1968                 MaximumConstraintError = SIMP ( statut = "f",
1969                                                 typ = "R",
1970                                                 max = 1,
1971                                                 val_min = 0.0,
1972                                                 fr = "Valeur maximum absolue de la fonction moins la valeur du niveau",
1973                                                 ang = "Maximum absolute value of the constraint function minus the level value",
1974                                                 ),
1975
1976                 ImportanceSampling = SIMP ( statut = "o",
1977                                             typ = 'TXM',
1978                                             into = ( 'yes', 'no' ),
1979                                             defaut = 'no',
1980                                             max = 1,
1981                                             fr = "Tirage d'importance au point de conception",
1982                                             ang = "Importance sampling at design point",
1983                                             ),
1984
1985                 FORM = BLOC ( condition = " Approximation in ( 'FORM', ) ",
1986
1987                     Probability = SIMP ( statut = "o",
1988                                          typ = 'TXM',
1989                                          into = ( 'yes', ),
1990                                          defaut = 'yes',
1991                                          max = 1,
1992                                          fr = "Probabiblite",
1993                                          ang = "Probability",
1994                                          ),
1995
1996                     DesignPoint = SIMP ( statut = "o",
1997                                          typ = 'TXM',
1998                                          into = ( 'yes', 'no' ),
1999                                          defaut = 'yes',
2000                                          max = 1,
2001                                          fr = "Point de conception",
2002                                          ang = "Design point",
2003                                          ),
2004
2005                     HasReliabilityIndex = SIMP ( statut = "o",
2006                                                  typ = 'TXM',
2007                                                  into = ( 'yes', 'no' ),
2008                                                  defaut = 'yes',
2009                                                  max = 1,
2010                                                  fr = "Indice de fiabilite",
2011                                                  ang = "Reliability index",
2012                                                  ),
2013
2014                     ImportanceFactor = SIMP ( statut = "o",
2015                                               typ = 'TXM',
2016                                               into = ( 'yes', 'no' ),
2017                                               defaut = 'no',
2018                                               max = 1,
2019                                               fr = "Facteur d'importance pour variable de sortie scalaire",
2020                                               ang = "ImportanceFactor",
2021                                               ),
2022
2023                     ImportanceFactorSettings = BLOC ( condition = " ImportanceFactor in ( 'yes', ) ",
2024
2025                             NumericalResults  = SIMP ( statut = "o",
2026                                                        typ = 'TXM',
2027                                                        into = ( 'yes', 'no' ),
2028                                                        defaut = 'yes',
2029                                                        max = 1,
2030                                                        fr = "Resultats numeriques",
2031                                                        ang = "NumericalResults",
2032                                                        ),
2033
2034                              GraphicalResults  = SIMP ( statut = "o",
2035                                                         typ = 'TXM',
2036                                                         into = ( 'yes', 'no' ),
2037                                                         defaut = 'no',
2038                                                         max = 1,
2039                                                         fr = "Resultats graphiques",
2040                                                         ang = "GraphicalResults",
2041                                                         ),
2042
2043                     ), # Fin BLOC ImportanceFactorSettings
2044
2045
2046                     SensitivityAnalysis = SIMP ( statut = "o",
2047                                                  typ = 'TXM',
2048                                                  into = ( 'yes', 'no' ),
2049                                                  defaut = 'no',
2050                                                  max = 1,
2051                                                  fr = "Analyse de sensibilite",
2052                                                  ang = "Sensitivity analysis",
2053                                                  ),
2054
2055                     SensitivityAnalysisSettings = BLOC ( condition = " SensitivityAnalysis in ( 'yes', ) ",
2056
2057                             HasoferReliabilityIndex = SIMP ( statut = "o",
2058                                                              typ = 'TXM',
2059                                                              into = ( 'yes', 'no' ),
2060                                                              defaut = 'no',
2061                                                              max = 1,
2062                                                              fr = "Indice de fiabilite de Hasofer",
2063                                                              ang = "Hasofer reliability index",
2064                                                              ),
2065
2066                             HasoferReliabilityIndexSettings = BLOC ( condition = " HasoferReliabilityIndex in ( 'yes', ) ",
2067
2068                                     NumericalResults  = SIMP ( statut = "o",
2069                                                                typ = 'TXM',
2070                                                                into = ( 'yes', 'no' ),
2071                                                                defaut = 'yes',
2072                                                                max = 1,
2073                                                                fr = "Resultats numeriques",
2074                                                                ang = "NumericalResults",
2075                                                                ),
2076
2077                                      GraphicalResults  = SIMP ( statut = "o",
2078                                                                 typ = 'TXM',
2079                                                                 into = ( 'yes', 'no' ),
2080                                                                 defaut = 'no',
2081                                                                 max = 1,
2082                                                                 fr = "Resultats graphiques",
2083                                                                 ang = "GraphicalResults",
2084                                                                 ),
2085
2086                             ), # Fin BLOC HasoferReliabilityIndexSettings
2087
2088                     ), # Fin BLOC SensitivityAnalysisSettings
2089
2090                     FunctionCallsNumber = SIMP ( statut = "o",
2091                                                  typ = 'TXM',
2092                                                  into = ( 'yes', 'no' ),
2093                                                  defaut = 'no',
2094                                                  max = 1,
2095                                                  fr = "Nombre d'appels a la fonction",
2096                                                  ang = "Function calls number",
2097                                                  ),
2098
2099
2100                 ), # Fin BLOC FORM
2101
2102
2103                 SORM = BLOC ( condition = " Approximation in ( 'SORM', ) ",
2104
2105
2106                     TvedtApproximation = SIMP ( statut = "o",
2107                                                 typ = 'TXM',
2108                                                 into = ( 'yes', 'no' ),
2109                                                 defaut = 'no',
2110                                                 max = 1,
2111                                                 fr = "Approximation de Tvedt",
2112                                                 ang = "Tvedt approximation",
2113                                                 ),
2114
2115                     HohenBichlerApproximation = SIMP ( statut = "o",
2116                                                        typ = 'TXM',
2117                                                        into = ( 'yes', 'no' ),
2118                                                        defaut = 'no',
2119                                                        max = 1,
2120                                                        fr = "Approximation de HohenBichler",
2121                                                        ang = "HohenBichler approximation",
2122                                                        ),
2123
2124                     BreitungApproximation = SIMP ( statut = "o",
2125                                                    typ = 'TXM',
2126                                                    into = ( 'yes', 'no' ),
2127                                                    defaut = 'no',
2128                                                    max = 1,
2129                                                    fr = "Approximation de Breitung",
2130                                                    ang = "Breitung approximation",
2131                                                    ),
2132
2133                     DesignPoint = SIMP ( statut = "o",
2134                                          typ = 'TXM',
2135                                          into = ( 'yes', 'no' ),
2136                                          defaut = 'yes',
2137                                          max = 1,
2138                                          fr = "Point de conception",
2139                                          ang = "Design point",
2140                                          ),
2141
2142                     ImportanceFactor = SIMP ( statut = "o",
2143                                               typ = 'TXM',
2144                                               into = ( 'yes', 'no' ),
2145                                               defaut = 'no',
2146                                               max = 1,
2147                                               fr = "Facteur d'importance pour variable de sortie scalaire",
2148                                               ang = "ImportanceFactor",
2149                                               ),
2150
2151                     ImportanceFactorSettings = BLOC ( condition = " ImportanceFactor in ( 'yes', ) ",
2152
2153                             NumericalResults  = SIMP ( statut = "o",
2154                                                        typ = 'TXM',
2155                                                        into = ( 'yes', 'no' ),
2156                                                        defaut = 'yes',
2157                                                        max = 1,
2158                                                        fr = "Resultats numeriques",
2159                                                        ang = "NumericalResults",
2160                                                        ),
2161
2162                              GraphicalResults  = SIMP ( statut = "o",
2163                                                         typ = 'TXM',
2164                                                         into = ( 'yes', 'no' ),
2165                                                         defaut = 'no',
2166                                                         max = 1,
2167                                                         fr = "Resultats graphiques",
2168                                                         ang = "GraphicalResults",
2169                                                         ),
2170
2171                     ), # Fin BLOC ImportanceFactorSettings
2172
2173
2174                     SensitivityAnalysis = SIMP ( statut = "o",
2175                                                  typ = 'TXM',
2176                                                  into = ( 'yes', 'no' ),
2177                                                  defaut = 'no',
2178                                                  max = 1,
2179                                                  fr = "Analyse de sensibilite",
2180                                                  ang = "Sensitivity analysis",
2181                                                  ),
2182
2183                     SensitivityAnalysisSettings = BLOC ( condition = " SensitivityAnalysis in ( 'yes', ) ",
2184
2185                             HasoferReliabilityIndex = SIMP ( statut = "o",
2186                                                              typ = 'TXM',
2187                                                              into = ( 'yes', 'no' ),
2188                                                              defaut = 'no',
2189                                                              max = 1,
2190                                                              fr = "Indice de fiabilite de Hasofer",
2191                                                              ang = "Hasofer reliability index",
2192                                                              ),
2193
2194                             HasoferReliabilityIndexSettings = BLOC ( condition = " HasoferReliabilityIndex in ( 'yes', ) ",
2195
2196                                     NumericalResults  = SIMP ( statut = "o",
2197                                                                typ = 'TXM',
2198                                                                into = ( 'yes', 'no' ),
2199                                                                defaut = 'yes',
2200                                                                max = 1,
2201                                                                fr = "Resultats numeriques",
2202                                                                ang = "NumericalResults",
2203                                                                ),
2204
2205                                      GraphicalResults  = SIMP ( statut = "o",
2206                                                                 typ = 'TXM',
2207                                                                 into = ( 'yes', 'no' ),
2208                                                                 defaut = 'no',
2209                                                                 max = 1,
2210                                                                 fr = "Resultats graphiques",
2211                                                                 ang = "GraphicalResults",
2212                                                                 ),
2213
2214                             ), # Fin BLOC HasoferReliabilityIndexSettings
2215
2216                     ), # Fin BLOC SensitivityAnalysisSettings
2217
2218                     FunctionCallsNumber = SIMP ( statut = "o",
2219                                                  typ = 'TXM',
2220                                                  into = ( 'yes', 'no' ),
2221                                                  defaut = 'no',
2222                                                  max = 1,
2223                                                  fr = "Nombre d'appels a la fonction",
2224                                                  ang = "Function calls number",
2225                                                  ),
2226
2227
2228                 ), # Fin BLOC SORM
2229
2230
2231
2232         ), # Fin BLOC AnalyticalSettings
2233
2234
2235
2236   ), # Fin BLOC ThresholdExceedence
2237
2238
2239
2240 ) # Fin PROC CRITERIA
2241
2242
2243 #===============================
2244 # 5. Definition des parametres
2245 #===============================
2246 VARI = OPER ( nom = "VARI",
2247                       sd_prod = variable,
2248                       op = None,
2249                       fr = "Definitions des lois marginales utilisees par les variables d'entree",
2250                       type=SIMP(statut='f',defaut="IN",into=("IN","OUT"), typ = "TXM",)
2251               )
2252
2253 ESSAI=PROC(nom="ESSAI",
2254        op=None,
2255        fr="Essai",
2256        ang = "Test",
2257
2258        MALOI       = SIMP(statut='o',typ=(loi,),),
2259        MAVARIABLE  = SIMP(statut='o',typ=(variable,),),
2260 ) ;
2261
2262
2263
2264