Aster/totalmod.comm

   1 DEBUT();
   2 % LECTURE DES PARAMETRES
   3 %
   4
   5
   6 % LECTURE DES TEMPS
   7
   8 %cycle 1
   9
  10
  11 % LECTURE DU MAILLAGE ET DEFINITIONS DES PROPRIETES ELEMENTS FINIS
  12
  13 PRE_GIBI ();
  14 %
  15  ma1 = LIRE_MAILLAGE ();
  16 %
  17 &ma1 = DEFI_GROUP ( MAILLAGE : ma1
  18                    CREA_GROUP_NO : (GROUP_MA : ( fixati
  19                                                  appui
  20                                                  lias_1
  21                                                  lias_2
  22                                                  ebguide
  23                                                  maintien
  24                                                  main_sup
  25                                                  embo_sup
  26                                                  tub_guid
  27                                                  grilles
  28                                                  embo_inf
  29                                                  pieds
  30                                                  crayons ))
  31                   ) ;
  32 %
  33 % AFFECTATION DU MODELE
  34 %
  35  mo1 = AFFE_MODELE ( Maillage: ma1
  36            % crayons combustibles
  37                     AFFE : ( GROUP_MA :     crayons
  38                              PHENOMENE :    'MECANIQUE'
  39                              MODELISATION : 'POU_D_E'
  40                             )
  41            % tubes guides
  42                     AFFE : ( GROUP_MA :     tub_guid
  43                              PHENOMENE :    'MECANIQUE'
  44                              MODELISATION : 'POU_D_E'
  45                             )
  46            % embouts
  47                     AFFE : ( GROUP_MA :     (embo_sup, embo_inf)
  48                              PHENOMENE :    'MECANIQUE'
  49                              MODELISATION : 'DKT'
  50                             )
  51           % pieds
  52
  53                     AFFE : ( GROUP_MA :    pieds
  54                              PHENOMENE :    'MECANIQUE'
  55                              MODELISATION : 'POU_D_E'
  56                             )
  57
  58            % grilles
  59                     AFFE : ( GROUP_MA :     grilles
  60                              PHENOMENE :    'MECANIQUE'
  61                              MODELISATION : 'POU_D_E'
  62                             )
  63            % liaison grille crayon
  64                     AFFE : ( GROUP_MA :     li_gr_cr
  65                              PHENOMENE :    'MECANIQUE'
  66                              MODELISATION : 'DIS_TR'
  67                             )
  68            % liaison grille tube
  69                     AFFE : ( GROUP_MA :     li_gr_gu
  70                              PHENOMENE :    'MECANIQUE'
  71                              MODELISATION : 'DIS_TR'
  72                             )
  73            % systeme de maintien
  74                     AFFE : ( GROUP_MA :     maintien
  75                              PHENOMENE:'MECANIQUE',
  76                              MODELISATION:'POU_D_E'
  77                            % MODELISATION:'POU_D_T'
  78                             )
  79                    );
  80 %
  81 % AFFECTATION DES CARACTERISTIQUES GEOMETRIQUES ELEMENTAIRES
  82 %
  83  cara1 = AFFE_CARA_ELEM ( MODELE :  mo1
  84
  85 % Famille assemblage : numero_fam
  86 %
  87            % crayons
  88                         POUTRE : ( GROUP_MA : cray
  89                                    SECTION : 'CERCLE'
  90                                    CARA    : ( 'R','EP' )
  91                                    VALE    : ( 3.   3. )
  92                                   )
  93            % tube guide partie courante
  94                         POUTRE : ( GROUP_MA : guid
  95                                    SECTION : 'CERCLE'
  96                                    CARA    : ( 'R','EP' )
  97                                    VALE    : ( 3.   3. )
  98                                   )
  99            % tube guide partie retrecie
 100                         POUTRE : ( GROUP_MA : retre
 101                                    SECTION : 'CERCLE'
 102                                    CARA    : ( 'R','EP' )
 103                                    VALE    : ( 3.   3. )
 104                                   )
 105            % grilles
 106                         POUTRE :  ( GROUP_MA : ( gril_i        )
 107                                     SECTION : 'RECTANGLE'
 108                                     CARA    : ( 'HY','HZ' )
 109                                     VALE    : ( 3.   3.   )
 110                                   )
 111                         POUTRE :  ( GROUP_MA : ( gril_e  )
 112                                     SECTION : 'RECTANGLE'
 113                                     CARA    : ( 'HY','HZ' )
 114                                     VALE    : ( 3.   3.   )
 115                                   )
 116                         ORIENTATION : ( GROUP_MA : (gril_i  , gril_e )
 117                                         CARA : 'VECT_Y'
 118                                         VALE : (1. 0. 0.)
 119                                       )
 120            % systeme de maintien
 121                         POUTRE : ( GROUP_MA : (mainti_y , mainti_z )
 122                                    SECTION : 'RECTANGLE'
 123                                    CARA    : ( 'HY','HZ' )
 124                                    VALE    : ( 3.   3. )
 125                                   )
 126                         ORIENTATION : ( GROUP_MA : (mainti_y )
 127                                         CARA : 'VECT_Y'
 128                                         VALE : (0. 1. 0.)
 129                                       )
 130                         ORIENTATION : ( GROUP_MA : (mainti_z )
 131                                         CARA : 'VECT_Y'
 132                                         VALE : (0. 0. 1.)
 133                                       )
 134            % embout inferieur
 135                         COQUE : ( GROUP_MA : eboinf
 136                                   EPAIS : 2.
 137                                   ANGL_REP : 90.
 138                                   )
 139           % pieds
 140                      POUTRE :  ( GROUP_MA : pi_ple
 141                                           SECTION : 'GENERALE'
 142                      CARA : ( 'A' , 'IZ' , 'IY' , 'JX' )
 143                     VALE : (1. 2. 3. 4.)
 144                                       )
 145                      POUTRE :  ( GROUP_MA : pi_tro
 146                                           SECTION : 'GENERALE'
 147                      CARA : ( 'A' , 'IZ' , 'IY' , 'JX' )
 148                     VALE : (1. 2. 3. 4. )
 149                                 )
 150                         ORIENTATION : ( GROUP_MA : (pi_ple  pi_tro )
 151                                         CARA : 'VECT_Y'
 152                                         VALE : (0. 1. 0.)
 153                                       )
 154
 155            % embout superieur
 156                         COQUE : ( GROUP_MA  : ebosup
 157                                   EPAIS : 3.
 158                                   ANGL_REP : 90.
 159                                   )
 160
 161            % liaison grille/crayon
 162                         DISCRET : ( GROUP_MA : elasti
 163                                     CARA : 'K_TR_D_L'
 164                                     VALE : ( 1.
 165                                              2.
 166                                              3.
 167                                              4.
 168                                              5.
 169                                              6.
 170                                            )
 171                                     REPERE : 'LOCAL'
 172                                   )
 173            % liaison grille/tube guide
 174                         DISCRET : ( GROUP_MA : rigid
 175                                     CARA : 'K_TR_D_L'
 176                                     VALE : ( 1.
 177                                              2.
 178                                              3.
 179                                              4.
 180                                              5.
 181                                              6.
 182                                            )
 183                                     REPERE : 'LOCAL'
 184                                   )
 185
 186                         ORIENTATION : ( GROUP_MA : (elasti
 187                                                     rigid        )
 188                                         CARA : 'VECT_Y'
 189                                         VALE : (1. 0. 0.) )
 190
 191     );
 192 %
 193 % DEFINITION DES CHARGES ET CONDITIONS LIMITES
 194 % DEFINITION DES TEMPERATURES NODALES EVOLUTIVES
 195 %
 196  f_tp1_1 = DEFI_FONCTION ( NOM_PARA: 'X'
 197                             NOM_RESU : 'TEMP'
 198                             PROL_DROIT: 'CONSTANT'
 199                             PROL_GAUCHE: 'CONSTANT'
 200                             VALE : ( 0. 1.
 201                                      5. 4.  )
 202                          );
 203  f_tp2_1 = DEFI_FONCTION ( NOM_PARA: 'X'
 204                           NOM_RESU : 'TEMP'
 205                           PROL_DROIT: 'CONSTANT'
 206                           PROL_GAUCHE: 'CONSTANT'
 207                           VALE : ( 0. 1.
 208                                    5. 4. )
 209                          );
 210  f_tp3_1 = DEFI_FONCTION ( NOM_PARA: 'X'
 211                           NOM_RESU : 'TEMP'
 212                           PROL_DROIT: 'CONSTANT'
 213                           PROL_GAUCHE: 'CONSTANT'
 214                           VALE : ( 0. 1.
 215                                    5. 4. )
 216                          );
 217
 218  chtem1_1 = AFFE_CHAM_NO ( MAILLAGE : ma1
 219                          GRANDEUR : 'TEMP_F'
 220                          AFFE : ( GROUP_NO: (maintien, pieds
 221                                              tub_guid, grilles
 222                                              crayons)
 223                                   NOM_CMP : 'TEMP'
 224                                   FONCTION : f_tp1_1 )
 225                          AFFE : ( GROUP_NO: (embo_sup, embo_inf)
 226                                   NOM_CMP : ('TEMP_INF' 'TEMP'
 227                                              'TEMP_SUP')
 228                                   FONCTION :(f_tp1_1,f_tp1_1,f_tp1_1))
 229                         );
 230
 231  chtem2_1 = AFFE_CHAM_NO ( MAILLAGE : ma1
 232                          GRANDEUR : 'TEMP_F'
 233                          AFFE : ( GROUP_NO: (maintien, pieds
 234                                              tub_guid, grilles
 235                                              crayons)
 236                                   NOM_CMP : 'TEMP'
 237                                   FONCTION : f_tp2_1 )
 238                          AFFE : ( GROUP_NO: (embo_sup, embo_inf)
 239                                   NOM_CMP : ('TEMP_INF' 'TEMP'
 240                                              'TEMP_SUP')
 241                                   FONCTION :(f_tp2_1,f_tp2_1,f_tp2_1))
 242                         );
 243
 244  chtem3_1 = AFFE_CHAM_NO ( MAILLAGE : ma1
 245                          GRANDEUR : 'TEMP_F'
 246                          AFFE : ( GROUP_NO: (maintien, pieds
 247                                              tub_guid, grilles)
 248                                   NOM_CMP : 'TEMP'
 249                                   FONCTION : f_tp3_1 )
 250                          AFFE : ( GROUP_NO: (embo_sup, embo_inf)
 251                                   NOM_CMP : ('TEMP_INF' 'TEMP'
 252                                              'TEMP_SUP')
 253                                   FONCTION : (f_tp3_1,f_tp3_1,f_tp3_1))
 254                          AFFE : ( GROUP_NO: crayons
 255                                   NOM_CMP : 'TEMP'
 256                                   FONCTION :f_tp3_1 )
 257                         );
 258
 259  chth_1 = CREA_RESU ( TYPE_RESU : 'EVOL_THER'
 260                     NOM_CHAM : 'TEMP'
 261                     CHAM_GD : ( CHAM_NO : chtem1_1
 262                                 INST : 0. )
 263                     CHAM_GD : ( CHAM_NO : chtem1_1
 264                                 INST : 1000. )
 265                     CHAM_GD : ( CHAM_NO : chtem2_1
 266                                 INST :2000.  )
 267                     CHAM_GD : ( CHAM_NO : chtem3_1
 268                                 INST : 3000.  )
 269                     CHAM_GD : ( CHAM_NO : chtem3_1
 270                                 INST : 4000.    )
 271                     CHAM_GD : ( CHAM_NO : chtem2_1
 272                                 INST : 5000.   )
 273                     CHAM_GD : ( CHAM_NO : chtem1_1
 274                                 INST : 6000.   )
 275                     CHAM_GD : ( CHAM_NO : chtem1_1
 276                                 INST : 7000.    )
 277                        );
 278
 279 % Chargement: charges thermo-mecaniques
 280 % ----------------------------------------
 281
 282  therm_1 = AFFE_CHAR_MECA ( MODELE: mo1
 283                              TEMP_CALCULEE : chth_1
 284                              );
 285
 286
 287
 288 % DEFINITION DES CONDITIONS LIMITES
 289 %            DES CONDITIONS DE CONTACT
 290 %
 291 % Chargement: pesanteur
 292 %
 293  pesant_1 = AFFE_CHAR_MECA ( MODELE : mo1
 294                              PESANTEUR : (9.81 -1. 0. 0. )
 295                             );
 296
 297 %
 298 % Conditions aux limites
 299 %
 300  cl_per_1 = AFFE_CHAR_MECA ( MODELE : mo1
 301                 % embout inferieur encastrement
 302                          DDL_IMPO : ( GROUP_NO : fixati
 303                                       DX : 0.  DY:0.    DZ : 0.
 304                                       DRX : 0. DRY : 0. DRZ : 0.
 305                                     )
 306                 % embout inferieur appui simple
 307                          DDL_IMPO : ( GROUP_NO : appui
 308                                       DX : 0.
 309                                     )
 310                 % embout superieur liaison guidee
 311                          DDL_IMPO : ( GROUP_NO : ebguide
 312                                       DY:0. DZ : 0.
 313                                       DRX : 0. DRY : 0. DRZ : 0.
 314                                     )
 315                 % liaison systeme de maintien embout superieur
 316                          DDL_IMPO : ( GROUP_NO : (lias_1, lias_2)
 317                                       DRX : 0.
 318                                     )
 319                 % liaison systeme de maintien embout superieur
 320                          LIAISON_GROUP : ( GROUP_NO_1 : lias_1
 321                                            GROUP_NO_2 : lias_2
 322                                            DDL_1 : 'DX' DDL_2 : 'DX'
 323                                            COEF_MULT_1 : 1.
 324                                            COEF_MULT_2 : -1.
 325                                            COEF_IMPO : 0. )
 326                          LIAISON_GROUP : ( GROUP_NO_1 : lias_1
 327                                            GROUP_NO_2 : lias_2
 328                                            DDL_1 : 'DY' DDL_2 : 'DY'
 329                                            COEF_MULT_1 : 1.
 330                                            COEF_MULT_2 : -1.
 331                                            COEF_IMPO : 0. )
 332                          LIAISON_GROUP : ( GROUP_NO_1 : lias_1
 333                                            GROUP_NO_2 : lias_2
 334                                            DDL_1 : 'DZ' DDL_2 : 'DZ'
 335                                            COEF_MULT_1 : 1.
 336                                            COEF_MULT_2 : -1.
 337                                            COEF_IMPO : 0. )
 338                          LIAISON_GROUP : ( GROUP_NO_1 : lias_1
 339                                            GROUP_NO_2 : lias_2
 340                                            DDL_1 : 'DRY' DDL_2 : 'DRY'
 341                                            COEF_MULT_1 : 1.
 342                                            COEF_MULT_2 : -1.
 343                                            COEF_IMPO : 0. )
 344                          LIAISON_GROUP : ( GROUP_NO_1 : lias_1
 345                                            GROUP_NO_2 : lias_2
 346                                            DDL_1 : 'DRZ' DDL_2 : 'DRZ'
 347                                            COEF_MULT_1 : 1.
 348                                            COEF_MULT_2 : -1.
 349                                            COEF_IMPO : 0. )
 350                             );
 351 %
 352  cl_ctc = AFFE_CHAR_MECA ( MODELE : mo1
 353
 354
 355                          CONTACT : ( APPARIEMENT : 'NODAL'
 356                                      GROUP_MA_1 : eboinf
 357                                      GROUP_MA_2 : mcrbas
 358                                      VECT_NORM_2 : (-1. 0. 0.)
 359                                      )
 360                                      INFO : 1
 361                          );
 362 %
 363 % Deplacement du systeme de maintien
 364 %
 365 cldi_f1 = DEFI_FONCTION ( NOM_PARA   : 'INST'
 366                            PROL_DROIT : 'CONSTANT'
 367                            PROL_GAUCHE: 'CONSTANT'
 368                            VALE : ( 0.     0.
 369                                     1000.  0.1
 370                                     2000.  0.2
 371                                     3000.  0.2
 372                                     4000.  0.2
 373                                     5000.  0.2
 374                                     6000.  0.1
 375                                     7000.  0.1      )
 376                          );
 377
 378  cl_di_1  = AFFE_CHAR_MECA ( MODELE : mo1
 379                 % deplacement du systeme de maintien
 380                          DDL_IMPO : ( GROUP_NO : main_sup
 381                                       DX : -1.
 382                                     )
 383                             );
 384
 385 % DEFINITION DES FORCES EXTERIEURES
 386 %
 387 % Fonction multiplicative de la force hydrodynamique normee
 388 %
 389 hydr_f1  = DEFI_FONCTION ( NOM_PARA    : 'INST'
 390                             PROL_DROIT  : 'CONSTANT'
 391                             PROL_GAUCHE : 'CONSTANT'
 392                             VALE        : ( 0.     0.
 393                                             1000.  0.1
 394                                             2000.  0.2
 395                                             3000.  0.2
 396                                             4000.  0.2
 397                                             5000.  0.2
 398                                             6000.  0.1
 399                                             7000.  0.1  )
 400                           );
 401 %
 402 % Profil des forces hydrodynamiques tubes guides a froid
 403 %
 404
 405
 406  hyfrtu_1   = DEFI_FONCTION (
 407                NOM_PARA    : 'X'
 408                VALE        : ( 0. 1.
 409                                5. 2. )
 410                PROL_DROIT  : 'CONSTANT'
 411                PROL_GAUCHE : 'CONSTANT'
 412                            );
 413 %
 414 % Profil des forces hydrodynamiques crayon a froid
 415 %
 416
 417
 418  hyfrcr_1   = DEFI_FONCTION (
 419                NOM_PARA    : 'X'
 420                VALE        : (0. 1.
 421                               5. 2. )
 422                PROL_DROIT  : 'CONSTANT'
 423                PROL_GAUCHE : 'CONSTANT'
 424                            );
 425 %
 426
 427 % Chargement : Forces hydrodynamiques
 428 %
 429  hydr_1   = AFFE_CHAR_MECA ( MODELE : mo1
 430
 431              % embout superieur
 432                           FORCE_COQUE  : ( GROUP_MA : ebosup
 433                                          FX : 10. )
 434              % embout inferieur
 435                           FORCE_COQUE  : ( GROUP_MA : eboinf
 436                                          FX : 10.  )
 437              % grille superieure
 438                           FORCE_POUTRE : ( GROUP_MA : (gril10 )
 439                                          FX : 10.   )
 440              % grille inferieure
 441                           FORCE_POUTRE : ( GROUP_MA : (gril1 )
 442                                          FX : 10.    )
 443              % grille de melange
 444                           FORCE_POUTRE : ( GROUP_MA : (gril2 )
 445                                          FX : 10.    )
 446                           FORCE_POUTRE : ( GROUP_MA : (gril3 )
 447                                          FX : 10.    )
 448                           FORCE_POUTRE : ( GROUP_MA : (gril4 )
 449                                          FX : 10.    )
 450                           FORCE_POUTRE : ( GROUP_MA : (gril5 )
 451                                          FX : 10.   )
 452                           FORCE_POUTRE : ( GROUP_MA : (gril6 )
 453                                          FX : 10.   )
 454                           FORCE_POUTRE : ( GROUP_MA : (gril7 )
 455                                          FX : 10.    )
 456                           FORCE_POUTRE : ( GROUP_MA : (gril8 )
 457                                          FX : 10.    )
 458                           FORCE_POUTRE : ( GROUP_MA : (gril9 )
 459                                          FX : 10.    )
 460                           );
 461  toto_1   = DEFI_FONCTION (
 462                NOM_PARA    : 'X'
 463                VALE        : ( 0. 10.
 464                                5. 10.  )
 465                PROL_DROIT  : 'CONSTANT'
 466                PROL_GAUCHE : 'CONSTANT'
 467                            );
 468 %
 469  fohydr_1 = AFFE_CHAR_MECA_F ( MODELE : mo1
 470
 471              % tubes guides et tube d'instrumentation
 472                           FORCE_POUTRE : ( GROUP_MA : (guid        )
 473                                            FX : toto_1  )
 474              % crayons
 475                           FORCE_POUTRE : ( GROUP_MA : cray
 476                                            FX :  toto_1  )
 477                             );
 478 %
 479 % Profil des forces d'Archimede tube guide
 480 %
 481
 482
 483  arfrtu_1   = DEFI_FONCTION (
 484                NOM_PARA    : 'X'
 485                VALE        : ( 0. 10.
 486                                5. 10.  )
 487                PROL_DROIT  : 'CONSTANT'
 488                PROL_GAUCHE : 'CONSTANT'
 489                            );
 490 %
 491 % Profil des forces d'Archimede crayon
 492 %
 493
 494
 495  arfrcr_1   = DEFI_FONCTION (
 496                NOM_PARA    : 'X'
 497                VALE        : ( 0. 10.
 498                                5. 10. )
 499                PROL_DROIT  : 'CONSTANT'
 500                PROL_GAUCHE : 'CONSTANT'
 501                            );
 502 %
 503 % Fonction multiplicative de la force d'Archimede
 504 %
 505 arch_f1  = DEFI_FONCTION ( NOM_PARA    : 'INST'
 506                             PROL_DROIT  : 'CONSTANT'
 507                             PROL_GAUCHE : 'CONSTANT'
 508                             VALE        : ( 0.     0.
 509                                             1000.  0.1
 510                                             2000.  0.2
 511                                             3000.  0.2
 512                                             4000.  0.2
 513                                             5000.  0.2
 514                                             6000.  0.1
 515                                             7000.  0.1 )
 516                           );
 517 %
 518
 519 %
 520 % Chargement : Forces d'Archimede
 521 % ----------------------------------------
 522 %
 523    arch_1 = AFFE_CHAR_MECA ( MODELE : mo1
 524              % systeme de maintien
 525                      FORCE_POUTRE : ( GROUP_MA : maintien
 526                                FX : 10. )
 527              % embout superieur
 528                           FORCE_COQUE  : ( GROUP_MA : ebosup
 529                                     FX : 10. )
 530              % embout inferieur
 531                           FORCE_COQUE  : ( GROUP_MA : eboinf
 532                                     FX : 10. )
 533              % grille superieure
 534                           FORCE_POUTRE : ( GROUP_MA : (gril10  )
 535                                     FX : 10.  )
 536              % grille inferieure
 537                           FORCE_POUTRE : ( GROUP_MA : (gril1 )
 538                                            FX : 10.  )
 539              % grille de melange
 540                           FORCE_POUTRE : ( GROUP_MA : (gril2 )
 541                                            FX : 10.  )
 542                           FORCE_POUTRE : ( GROUP_MA : (gril3 )
 543                                            FX : 10.   )
 544                           FORCE_POUTRE : ( GROUP_MA : (gril4 )
 545                                            FX : 10.   )
 546                           FORCE_POUTRE : ( GROUP_MA : (gril5 )
 547                                            FX : 10.   )
 548                           FORCE_POUTRE : ( GROUP_MA : (gril6 )
 549                                            FX : 10.   )
 550                           FORCE_POUTRE : ( GROUP_MA : (gril7 )
 551                                            FX : 10.   )
 552                           FORCE_POUTRE : ( GROUP_MA : (gril8 )
 553                                            FX : 10.   )
 554                           FORCE_POUTRE : ( GROUP_MA : (gril9 )
 555                                            FX : 10.  )
 556                           );
 557  foarch_1 = AFFE_CHAR_MECA_F ( MODELE : mo1
 558
 559              % tubes guides et tube d'instrumentation
 560                           FORCE_POUTRE : ( GROUP_MA : guid
 561                                            FX : arch_f1  )
 562              % crayons
 563                           FORCE_POUTRE : ( GROUP_MA : cray
 564                                            FX :arch_f1  )
 565                             );
 566 %
 567 %
 568 % FORCE DE SERRAGE DE LA LIAISON GRILLE - CRAYON
 569 %
 570 %
 571
 572 %
 573 %
 574 %
 575 % definition de la liste d'instants
 576 %
 577   l_frig   = DEFI_LIST_REEL( Debut: 0.0
 578                              Intervalle : (Jusqu_a: 1000. Nombre: 100)
 579                            );
 580 %
 581 % Inrterpretation des fonctions
 582 %
 583 f1=DEFI_FONCTION ( NOM_PARA    : 'INST'
 584                          PROL_DROIT  : 'CONSTANT'
 585                          PROL_GAUCHE : 'CONSTANT'
 586                          VALE : (
 587                          0.00000E+00   0.00000E+00
 588                          4.32000E+04   0.00000E+00 ) );
 589 f2=DEFI_FONCTION ( NOM_PARA    : 'INST'
 590                          PROL_DROIT  : 'CONSTANT'
 591                          PROL_GAUCHE : 'CONSTANT'
 592                          VALE : (
 593                          0.00000E+00   0.00000E+00
 594                          4.32000E+04   0.00000E+00 ) );
 595 f3=DEFI_FONCTION ( NOM_PARA    : 'INST'
 596                          PROL_DROIT  : 'CONSTANT'
 597                          PROL_GAUCHE : 'CONSTANT'
 598                          VALE : (
 599                          0.00000E+00   0.00000E+00
 600                          4.32000E+04   0.00000E+00 ) );
 601 f4=DEFI_FONCTION ( NOM_PARA    : 'INST'
 602                          PROL_DROIT  : 'CONSTANT'
 603                          PROL_GAUCHE : 'CONSTANT'
 604                          VALE : (
 605                          0.00000E+00   0.00000E+00
 606                          4.32000E+04   0.00000E+00 ) );
 607 f5=DEFI_FONCTION ( NOM_PARA    : 'INST'
 608                          PROL_DROIT  : 'CONSTANT'
 609                          PROL_GAUCHE : 'CONSTANT'
 610                          VALE : (
 611                          0.00000E+00   0.00000E+00
 612                          4.32000E+04   0.00000E+00 ) );
 613 f6=DEFI_FONCTION ( NOM_PARA    : 'INST'
 614                          PROL_DROIT  : 'CONSTANT'
 615                          PROL_GAUCHE : 'CONSTANT'
 616                          VALE : (
 617                          0.00000E+00   0.00000E+00
 618                          4.32000E+04   0.00000E+00 ) );
 619 f7=DEFI_FONCTION ( NOM_PARA    : 'INST'
 620                          PROL_DROIT  : 'CONSTANT'
 621                          PROL_GAUCHE : 'CONSTANT'
 622                          VALE : (
 623                          0.00000E+00   0.00000E+00
 624                          4.32000E+04   0.00000E+00 ) );
 625 f8=DEFI_FONCTION ( NOM_PARA    : 'INST'
 626                          PROL_DROIT  : 'CONSTANT'
 627                          PROL_GAUCHE : 'CONSTANT'
 628                          VALE : (
 629                          0.00000E+00   0.00000E+00
 630                          4.32000E+04   0.00000E+00 ) );
 631 f9=DEFI_FONCTION ( NOM_PARA    : 'INST'
 632                          PROL_DROIT  : 'CONSTANT'
 633                          PROL_GAUCHE : 'CONSTANT'
 634                          VALE : (
 635                          0.00000E+00   0.00000E+00
 636                          4.32000E+04   0.00000E+00 ) );
 637 f10=DEFI_FONCTION ( NOM_PARA    : 'INST'
 638                          PROL_DROIT  : 'CONSTANT'
 639                          PROL_GAUCHE : 'CONSTANT'
 640                          VALE : (
 641                          0.00000E+00   0.00000E+00
 642                          4.32000E+04   0.00000E+00 ) );
 643
 644  frig1 = CALC_FONC_INTERP ( FONCTION    : f1
 645                             LIST_PARA   : l_frig
 646                             PROL_DROIT  : 'CONSTANT'
 647                             PROL_GAUCHE : 'CONSTANT'
 648                             INTERPOL    : 'LIN'
 649                            );
 650  frig2 = CALC_FONC_INTERP ( FONCTION    : f2
 651                             LIST_PARA   : l_frig
 652                             PROL_DROIT  : 'CONSTANT'
 653                             PROL_GAUCHE : 'CONSTANT'
 654                             INTERPOL    : 'LIN'
 655                            );
 656  frig3 = CALC_FONC_INTERP ( FONCTION    : f3
 657                             LIST_PARA   : l_frig
 658                             PROL_DROIT  : 'CONSTANT'
 659                             PROL_GAUCHE : 'CONSTANT'
 660                             INTERPOL    : 'LIN'
 661                            );
 662  frig4 = CALC_FONC_INTERP ( FONCTION    : f4
 663                             LIST_PARA   : l_frig
 664                             PROL_DROIT  : 'CONSTANT'
 665                             PROL_GAUCHE : 'CONSTANT'
 666                             INTERPOL    : 'LIN'
 667                            );
 668  frig5 = CALC_FONC_INTERP ( FONCTION    : f5
 669                             LIST_PARA   : l_frig
 670                             PROL_DROIT  : 'CONSTANT'
 671                             PROL_GAUCHE : 'CONSTANT'
 672                              INTERPOL    : 'LIN'
 673                            );
 674  frig6 = CALC_FONC_INTERP ( FONCTION    : f6
 675                             LIST_PARA   : l_frig
 676                             PROL_DROIT  : 'CONSTANT'
 677                             PROL_GAUCHE : 'CONSTANT'
 678                             INTERPOL    : 'LIN'
 679                            );
 680  frig7 = CALC_FONC_INTERP ( FONCTION    : f7
 681                             LIST_PARA   : l_frig
 682                             PROL_DROIT  : 'CONSTANT'
 683                             PROL_GAUCHE : 'CONSTANT'
 684                             INTERPOL    : 'LIN'
 685                            );
 686  frig8 = CALC_FONC_INTERP ( FONCTION    : f8
 687                             LIST_PARA   : l_frig
 688                             PROL_DROIT  : 'CONSTANT'
 689                             PROL_GAUCHE : 'CONSTANT'
 690                             INTERPOL    : 'LIN'
 691                            );
 692  frig9 = CALC_FONC_INTERP ( FONCTION    : f9
 693                             LIST_PARA   : l_frig
 694                             PROL_DROIT  : 'CONSTANT'
 695                             PROL_GAUCHE : 'CONSTANT'
 696                             INTERPOL    : 'LIN'
 697                            );
 698  frig10= CALC_FONC_INTERP ( FONCTION    : f10
 699                             LIST_PARA   : l_frig
 700                             PROL_DROIT  : 'CONSTANT'
 701                             PROL_GAUCHE : 'CONSTANT'
 702                             INTERPOL    : 'LIN'
 703                            );
 704 %
 705 % DEFINITION DU FLUX NEUTRONIQUE RADIAL
 706
 707 flux_f1=DEFI_FONCTION ( NOM_PARA    : 'INST'
 708                          PROL_DROIT  : 'CONSTANT'
 709                          PROL_GAUCHE : 'CONSTANT'
 710                          VALE : (
 711                          0.00000E+00   0.00000E+00
 712                          4.32000E+04   0.00000E+00
 713                          9.74520E+05   4.53280E+05
 714                          2.63880E+06   2.07785E+06
 715                          4.08780E+06   3.51279E+06
 716                          5.42844E+06   4.85677E+06
 717                          7.24500E+06   6.67627E+06
 718                          8.85348E+06   8.27977E+06
 719                          1.04609E+07   9.87261E+06
 720                          1.22029E+07   1.15869E+07
 721                          1.39446E+07   1.32952E+07
 722                          1.41437E+07   1.34902E+07
 723                          1.55632E+07   1.48753E+07
 724                          1.72361E+07   1.65069E+07
 725                          1.79482E+07   1.71992E+07
 726                          1.81901E+07   1.74334E+07
 727                          1.81937E+07   1.74352E+07
 728                          1.81973E+07   1.74352E+07
 729                          1.82009E+07   1.74352E+07
 730                                 )
 731                        ) ;
 732 % Definition du flux axial pour le cycle 1
 733
 734 fluxax1 = LIRE_FONCTION ( UNITE       : 60
 735                            PROL_GAUCHE : 'CONSTANT'
 736                            PROL_DROIT  : 'CONSTANT'
 737                          );
 738
 739  fly_1_1  = DEFI_FONCTION (
 740                      NOM_PARA : 'Y'
 741                      PROL_DROIT  :'CONSTANT'
 742                      PROL_GAUCHE :'CONSTANT'
 743                      VALE : (
 744                 -1.0080000000000D-01    3.8032834757835D-07
 745                 -8.8200000000000D-02    3.8083689458689D-07
 746                 -7.5600000000000D-02    3.8296082621083D-07
 747                 -6.3000000000000D-02    3.8584009971510D-07
 748                 -5.0400000000000D-02    3.8946100427350D-07
 749                 -3.7800000000000D-02    3.9259455128205D-07
 750                 -2.5200000000000D-02    3.9258457977208D-07
 751                 -1.2600000000000D-02    3.9294106125356D-07
 752                0.    3.9399554843305D-07
 753                  1.2600000000000D-02    3.9293233618234D-07
 754                  2.5200000000000D-02    3.9256712962963D-07
 755                  3.7800000000000D-02    3.9256712962963D-07
 756                  5.0400000000000D-02    3.8942610398860D-07
 757                  6.3000000000000D-02    3.8579522792023D-07
 758                  7.5600000000000D-02    3.8290972222222D-07
 759                  8.8200000000000D-02    3.8077706552707D-07
 760                  1.0080000000000D-01    3.8026103988604D-07
 761                      ) ) ;
 762
 763  fly_1_2  = DEFI_FONCTION (
 764                      NOM_PARA : 'Y'
 765                      PROL_DROIT  :'CONSTANT'
 766                      PROL_GAUCHE :'CONSTANT'
 767                      VALE : (
 768                 -1.0080000000000D-01    4.3595566239316D-07
 769                 -8.8200000000000D-02    4.3866542022792D-07
 770                 -7.5600000000000D-02    4.4331463675214D-07
 771                 -6.3000000000000D-02    4.4911556267806D-07
 772                 -5.0400000000000D-02    4.5750658831909D-07
 773                 -3.7800000000000D-02    4.6647222222222D-07
 774                 -2.5200000000000D-02    4.6002938034188D-07
 775                 -1.2600000000000D-02    4.6004309116809D-07
 776                0.    4.6686235754986D-07
 777                  1.2600000000000D-02    4.6003311965812D-07
 778                  2.5200000000000D-02    4.6000943732194D-07
 779                  3.7800000000000D-02    4.6643981481481D-07
 780                  5.0400000000000D-02    4.5746420940171D-07
 781                  6.3000000000000D-02    4.4906445868946D-07
 782                  7.5600000000000D-02    4.4325356125356D-07
 783                  8.8200000000000D-02    4.3859686609687D-07
 784                  1.0080000000000D-01    4.3587838319088D-07
 785                      ) ) ;
 786
 787  fly_1_3  = DEFI_FONCTION (
 788                      NOM_PARA : 'Y'
 789                      PROL_DROIT  :'CONSTANT'
 790                      PROL_GAUCHE :'CONSTANT'
 791                      VALE : (
 792                 -1.0080000000000D-01    4.9431891025641D-07
 793                 -8.8200000000000D-02    5.0005128205128D-07
 794                 -7.5600000000000D-02    5.1337197293447D-07
 795                 -6.3000000000000D-02    5.2915313390313D-07
 796                 -5.0400000000000D-02    5.3788069800570D-07
 797                 -3.7800000000000D-02    5.3504629629630D-07
 798                 -2.5200000000000D-02    5.3221189458689D-07
 799                 -1.2600000000000D-02    5.3136431623932D-07
 800                0.    5.3135870726496D-07
 801                  1.2600000000000D-02    5.3135309829060D-07
 802                  2.5200000000000D-02    5.3218696581197D-07
 803                  3.7800000000000D-02    5.3500827991453D-07
 804                  5.0400000000000D-02    5.3782959401709D-07
 805                  6.3000000000000D-02    5.2909205840456D-07
 806                  7.5600000000000D-02    5.1330217236467D-07
 807                  8.8200000000000D-02    4.9997150997151D-07
 808                  1.0080000000000D-01    4.9423165954416D-07
 809                      ) ) ;
 810
 811  fly_1_4  = DEFI_FONCTION (
 812                      NOM_PARA : 'Y'
 813                      PROL_DROIT  :'CONSTANT'
 814                      PROL_GAUCHE :'CONSTANT'
 815                      VALE : (
 816                 -1.0080000000000D-01    5.5489957264957D-07
 817                 -8.8200000000000D-02    5.6477884615385D-07
 818                 -7.5600000000000D-02    5.9007781339031D-07
 819                 -6.3000000000000D-02    5.9726415598291D-07
 820                 -5.0400000000000D-02    6.0445049857550D-07
 821                 -3.7800000000000D-02    5.9840776353276D-07
 822                 -2.5200000000000D-02    5.8440651709402D-07
 823                 -1.2600000000000D-02    5.8263906695157D-07
 824                0.    5.9153240740741D-07
 825                  1.2600000000000D-02    5.8262660256410D-07
 826                  2.5200000000000D-02    5.8438034188034D-07
 827                  3.7800000000000D-02    5.9836663105413D-07
 828                  5.0400000000000D-02    6.0439316239316D-07
 829                  6.3000000000000D-02    5.9719497863248D-07
 830                  7.5600000000000D-02    5.8999679487179D-07
 831                  8.8200000000000D-02    5.6469159544160D-07
 832                  1.0080000000000D-01    5.5480110398860D-07
 833                      ) ) ;
 834
 835  fly_1_5  = DEFI_FONCTION (
 836                      NOM_PARA : 'Y'
 837                      PROL_DROIT  :'CONSTANT'
 838                      PROL_GAUCHE :'CONSTANT'
 839                      VALE : (
 840                 -1.0080000000000D-01    6.1267200854701D-07
 841                 -8.8200000000000D-02    6.2962357549858D-07
 842                 -7.5600000000000D-02    6.5655413105413D-07
 843                 -6.3000000000000D-02    6.6164583333333D-07
 844                 -5.0400000000000D-02    6.5190242165242D-07
 845                 -3.7800000000000D-02    6.5212678062678D-07
 846                 -2.5200000000000D-02    6.3746616809117D-07
 847                 -1.2600000000000D-02    6.3547435897436D-07
 848                0.    6.4509686609687D-07
 849                  1.2600000000000D-02    6.3546064814815D-07
 850                  2.5200000000000D-02    6.3743625356125D-07
 851                  3.7800000000000D-02    6.5207941595442D-07
 852                  5.0400000000000D-02    6.5183885327635D-07
 853                  6.3000000000000D-02    6.6156980056980D-07
 854                  7.5600000000000D-02    6.5646563390313D-07
 855                  8.8200000000000D-02    6.2952386039886D-07
 856                  1.0080000000000D-01    6.1256356837607D-07
 857                      ) ) ;
 858
 859  fly_1_6  = DEFI_FONCTION (
 860                      NOM_PARA : 'Y'
 861                      PROL_DROIT  :'CONSTANT'
 862                      PROL_GAUCHE :'CONSTANT'
 863                      VALE : (
 864                 -1.0080000000000D-01    6.5755128205128D-07
 865                 -8.8200000000000D-02    6.8336253561254D-07
 866                 -7.5600000000000D-02    6.9012072649573D-07
 867                 -6.3000000000000D-02    6.9687891737892D-07
 868                 -5.0400000000000D-02    6.9362571225071D-07
 869                 -3.7800000000000D-02    6.9074768518519D-07
 870                 -2.5200000000000D-02    6.8786965811966D-07
 871                 -1.2600000000000D-02    6.8586413817664D-07
 872                0.    6.8585603632479D-07
 873                  1.2600000000000D-02    6.8584793447293D-07
 874                  2.5200000000000D-02    6.8783600427350D-07
 875                  3.7800000000000D-02    6.9069720441595D-07
 876                  5.0400000000000D-02    6.9355840455840D-07
 877                  6.3000000000000D-02    6.9679789886040D-07
 878                  7.5600000000000D-02    6.9002412749288D-07
 879                  8.8200000000000D-02    6.8325035612536D-07
 880                  1.0080000000000D-01    6.5743162393162D-07
 881                      ) ) ;
 882
 883  fly_1_7  = DEFI_FONCTION (
 884                      NOM_PARA : 'Y'
 885                      PROL_DROIT  :'CONSTANT'
 886                      PROL_GAUCHE :'CONSTANT'
 887                      VALE : (
 888                 -1.0080000000000D-01    6.9684526353276D-07
 889                 -8.8200000000000D-02    7.1406606125356D-07
 890                 -7.5600000000000D-02    7.3236502849003D-07
 891                 -6.3000000000000D-02    7.2076442307692D-07
 892                 -5.0400000000000D-02    7.1793500712251D-07
 893                 -3.7800000000000D-02    7.2835648148148D-07
 894                 -2.5200000000000D-02    7.1444747150997D-07
 895                 -1.2600000000000D-02    7.1130893874644D-07
 896                0.    7.2107727920228D-07
 897                  1.2600000000000D-02    7.1129398148148D-07
 898                  2.5200000000000D-02    7.1441257122507D-07
 899                  3.7800000000000D-02    7.2830288461538D-07
 900                  5.0400000000000D-02    7.1786396011396D-07
 901                  6.3000000000000D-02    7.2067841880342D-07
 902                  7.5600000000000D-02    7.3226157407407D-07
 903                  8.8200000000000D-02    7.1395138888889D-07
 904                  1.0080000000000D-01    6.9671812678063D-07
 905                      ) ) ;
 906
 907  fly_1_8  = DEFI_FONCTION (
 908                      NOM_PARA : 'Y'
 909                      PROL_DROIT  :'CONSTANT'
 910                      PROL_GAUCHE :'CONSTANT'
 911                      VALE : (
 912                 -1.0080000000000D-01    7.3679362535613D-07
 913                 -8.8200000000000D-02    7.5423130341880D-07
 914                 -7.5600000000000D-02    7.7215633903134D-07
 915                 -6.3000000000000D-02    7.5865740740741D-07
 916                 -5.0400000000000D-02    7.5547774216524D-07
 917                 -3.7800000000000D-02    7.6662464387464D-07
 918                 -2.5200000000000D-02    7.5088461538462D-07
 919                 -1.2600000000000D-02    7.4313924501424D-07
 920                0.    7.4863230056980D-07
 921                  1.2600000000000D-02    7.4312054843305D-07
 922                  2.5200000000000D-02    7.5084722222222D-07
 923                  3.7800000000000D-02    7.6656730769231D-07
 924                  5.0400000000000D-02    7.5540170940171D-07
 925                  6.3000000000000D-02    7.5856267806268D-07
 926                  7.5600000000000D-02    7.7204415954416D-07
 927                  8.8200000000000D-02    7.5410790598291D-07
 928                  1.0080000000000D-01    7.3665776353276D-07
 929                      ) ) ;
 930
 931  fly_1_9  = DEFI_FONCTION (
 932                      NOM_PARA : 'Y'
 933                      PROL_DROIT  :'CONSTANT'
 934                      PROL_GAUCHE :'CONSTANT'
 935                      VALE : (
 936                 -1.0080000000000D-01    7.7370940170940D-07
 937                 -8.8200000000000D-02    8.0144764957265D-07
 938                 -7.5600000000000D-02    8.0374358974359D-07
 939                 -6.3000000000000D-02    8.0603952991453D-07
 940                 -5.0400000000000D-02    8.0238372507123D-07
 941                 -3.7800000000000D-02    7.9936672008547D-07
 942                 -2.5200000000000D-02    7.9634971509972D-07
 943                 -1.2600000000000D-02    7.8317361111111D-07
 944                0.    7.8316426282051D-07
 945                  1.2600000000000D-02    7.8315491452991D-07
 946                  2.5200000000000D-02    7.9631107549858D-07
 947                  3.7800000000000D-02    7.9930689102564D-07
 948                  5.0400000000000D-02    8.0230270655271D-07
 949                  6.3000000000000D-02    8.0593856837607D-07
 950                  7.5600000000000D-02    8.0362580128205D-07
 951                  8.8200000000000D-02    8.0131303418803D-07
 952                  1.0080000000000D-01    7.7356356837607D-07
 953                      ) ) ;
 954
 955  fly_1_10 = DEFI_FONCTION (
 956                      NOM_PARA : 'Y'
 957                      PROL_DROIT  :'CONSTANT'
 958                      PROL_GAUCHE :'CONSTANT'
 959                      VALE : (
 960                 -1.0080000000000D-01    7.9424198717949D-07
 961                 -8.8200000000000D-02    8.1225925925926D-07
 962                 -7.5600000000000D-02    8.3075516381766D-07
 963                 -6.3000000000000D-02    8.1540651709402D-07
 964                 -5.0400000000000D-02    8.1131071937322D-07
 965                 -3.7800000000000D-02    8.2297489316239D-07
 966                 -2.5200000000000D-02    8.0577403846154D-07
 967                 -1.2600000000000D-02    7.9717236467236D-07
 968                0.    8.0284615384615D-07
 969                  1.2600000000000D-02    7.9715117521368D-07
 970                  2.5200000000000D-02    8.0573165954416D-07
 971                  3.7800000000000D-02    8.2291257122507D-07
 972                  5.0400000000000D-02    8.1122845441595D-07
 973                  6.3000000000000D-02    8.1530306267806D-07
 974                  7.5600000000000D-02    8.3063051994302D-07
 975                  8.8200000000000D-02    8.1211841168091D-07
 976                  1.0080000000000D-01    7.9408742877493D-07
 977                      ) ) ;
 978
 979  fly_1_11 = DEFI_FONCTION (
 980                      NOM_PARA : 'Y'
 981                      PROL_DROIT  :'CONSTANT'
 982                      PROL_GAUCHE :'CONSTANT'
 983                      VALE : (
 984                 -1.0080000000000D-01    8.1647845441595D-07
 985                 -8.8200000000000D-02    8.3510648148148D-07
 986                 -7.5600000000000D-02    8.5487873931624D-07
 987                 -6.3000000000000D-02    8.3969088319088D-07
 988                 -5.0400000000000D-02    8.3502920227920D-07
 989                 -3.7800000000000D-02    8.4651014957265D-07
 990                 -2.5200000000000D-02    8.2973931623932D-07
 991                 -1.2600000000000D-02    8.2549269943020D-07
 992                0.    8.3637909544160D-07
 993                  1.2600000000000D-02    8.2547275641026D-07
 994                  2.5200000000000D-02    8.2969319800570D-07
 995                  3.7800000000000D-02    8.4644408831909D-07
 996                  5.0400000000000D-02    8.3493945868946D-07
 997                  6.3000000000000D-02    8.3957995014245D-07
 998                  7.5600000000000D-02    8.5474786324786D-07
 999                  8.8200000000000D-02    8.3495690883191D-07
1000                  1.0080000000000D-01    8.1631641737892D-07
1001                      ) ) ;
1002
1003  fly_1_12 = DEFI_FONCTION (
1004                      NOM_PARA : 'Y'
1005                      PROL_DROIT  :'CONSTANT'
1006                      PROL_GAUCHE :'CONSTANT'
1007                      VALE : (
1008                 -1.0080000000000D-01    8.3948272792023D-07
1009                 -8.8200000000000D-02    8.7002546296296D-07
1010                 -7.5600000000000D-02    8.7608564814815D-07
1011                 -6.3000000000000D-02    8.8214583333333D-07
1012                 -5.0400000000000D-02    8.7589245014245D-07
1013                 -3.7800000000000D-02    8.7128872863248D-07
1014                 -2.5200000000000D-02    8.6668500712251D-07
1015                 -1.2600000000000D-02    8.6321616809117D-07
1016                0.    8.6320432692308D-07
1017                  1.2600000000000D-02    8.6319248575499D-07
1018                  2.5200000000000D-02    8.6663764245014D-07
1019                  3.7800000000000D-02    8.7121705840456D-07
1020                  5.0400000000000D-02    8.7579647435897D-07
1021                  6.3000000000000D-02    8.8202617521368D-07
1022                  7.5600000000000D-02    8.7594667022792D-07
1023                  8.8200000000000D-02    8.6986716524217D-07
1024                  1.0080000000000D-01    8.3931196581197D-07
1025                      ) ) ;
1026
1027  fly_1_13 = DEFI_FONCTION (
1028                      NOM_PARA : 'Y'
1029                      PROL_DROIT  :'CONSTANT'
1030                      PROL_GAUCHE :'CONSTANT'
1031                      VALE : (
1032                 -1.0080000000000D-01    8.5508315527066D-07
1033                 -8.8200000000000D-02    8.7546616809117D-07
1034                 -7.5600000000000D-02    9.0944533475783D-07
1035                 -6.3000000000000D-02    9.1292663817664D-07
1036                 -5.0400000000000D-02    8.9653721509972D-07
1037                 -3.7800000000000D-02    8.9547275641026D-07
1038                 -2.5200000000000D-02    8.7405644586895D-07
1039                 -1.2600000000000D-02    8.7004789886040D-07
1040                0.    8.8226549145299D-07
1041                  1.2600000000000D-02    8.7002546296296D-07
1042                  2.5200000000000D-02    8.7400534188034D-07
1043                  3.7800000000000D-02    8.9539672364672D-07
1044                  5.0400000000000D-02    8.9643376068376D-07
1045                  6.3000000000000D-02    9.1279825498576D-07
1046                  7.5600000000000D-02    9.0929576210826D-07
1047                  8.8200000000000D-02    8.7530288461538D-07
1048                  1.0080000000000D-01    8.5490491452991D-07
1049                      ) ) ;
1050
1051  fly_1_14 = DEFI_FONCTION (
1052                      NOM_PARA : 'Y'
1053                      PROL_DROIT  :'CONSTANT'
1054                      PROL_GAUCHE :'CONSTANT'
1055                      VALE : (
1056                 -1.0080000000000D-01    8.5798237179487D-07
1057                 -8.8200000000000D-02    8.6880893874644D-07
1058                 -7.5600000000000D-02    9.0299750712251D-07
1059                 -6.3000000000000D-02    9.0948646723647D-07
1060                 -5.0400000000000D-02    9.1597542735043D-07
1061                 -3.7800000000000D-02    9.0500302706553D-07
1062                 -2.5200000000000D-02    8.8214084757835D-07
1063                 -1.2600000000000D-02    8.7779825498575D-07
1064                0.    8.8993732193732D-07
1065                  1.2600000000000D-02    8.7777207977208D-07
1066                  2.5200000000000D-02    8.8208974358974D-07
1067                  3.7800000000000D-02    9.0492325498576D-07
1068                  5.0400000000000D-02    9.1586574074074D-07
1069                  6.3000000000000D-02    9.0935434472934D-07
1070                  7.5600000000000D-02    9.0284294871795D-07
1071                  8.8200000000000D-02    8.6864066951567D-07
1072                  1.0080000000000D-01    8.5779665242165D-07
1073                      ) ) ;
1074
1075  fly_1_15 = DEFI_FONCTION (
1076                      NOM_PARA : 'Y'
1077                      PROL_DROIT  :'CONSTANT'
1078                      PROL_GAUCHE :'CONSTANT'
1079                      VALE : (
1080                 -1.0080000000000D-01    8.5773557692308D-07
1081                 -8.8200000000000D-02    8.6180893874644D-07
1082                 -7.5600000000000D-02    8.7862838319088D-07
1083                 -6.3000000000000D-02    8.9924074074074D-07
1084                 -5.0400000000000D-02    9.0872489316239D-07
1085                 -3.7800000000000D-02    9.0174545940171D-07
1086                 -2.5200000000000D-02    8.9476602564103D-07
1087                 -1.2600000000000D-02    8.9122364672365D-07
1088                0.    8.9120993589744D-07
1089                  1.2600000000000D-02    8.9119622507123D-07
1090                  2.5200000000000D-02    8.9471118233618D-07
1091                  3.7800000000000D-02    9.0166132478632D-07
1092                  5.0400000000000D-02    9.0861146723647D-07
1093                  6.3000000000000D-02    8.9910363247863D-07
1094                  7.5600000000000D-02    8.7847257834758D-07
1095                  8.8200000000000D-02    8.6163693019943D-07
1096                  1.0080000000000D-01    8.5754362535613D-07
1097                      ) ) ;
1098
1099  fly_1_16 = DEFI_FONCTION (
1100                      NOM_PARA : 'Y'
1101                      PROL_DROIT  :'CONSTANT'
1102                      PROL_GAUCHE :'CONSTANT'
1103                      VALE : (
1104                 -1.0080000000000D-01    8.5903935185185D-07
1105                 -8.8200000000000D-02    8.5690046296296D-07
1106                 -7.5600000000000D-02    8.5831018518519D-07
1107                 -6.3000000000000D-02    8.6167307692308D-07
1108                 -5.0400000000000D-02    8.7125320512821D-07
1109                 -3.7800000000000D-02    8.8563586182336D-07
1110                 -2.5200000000000D-02    8.7089298433048D-07
1111                 -1.2600000000000D-02    8.6842877492877D-07
1112                0.    8.7942111823362D-07
1113                  1.2600000000000D-02    8.6840135327635D-07
1114                  2.5200000000000D-02    8.7083814102564D-07
1115                  3.7800000000000D-02    8.8554985754986D-07
1116                  5.0400000000000D-02    8.7113977920228D-07
1117                  6.3000000000000D-02    8.6153846153846D-07
1118                  7.5600000000000D-02    8.5815313390313D-07
1119                  8.8200000000000D-02    8.5672471509972D-07
1120                  1.0080000000000D-01    8.5884116809117D-07
1121                      ) ) ;
1122
1123  fly_1_17 = DEFI_FONCTION (
1124                      NOM_PARA : 'Y'
1125                      PROL_DROIT  :'CONSTANT'
1126                      PROL_GAUCHE :'CONSTANT'
1127                      VALE : (
1128                 -1.0080000000000D-01    8.6388675213675D-07
1129                 -8.8200000000000D-02    8.5563782051282D-07
1130                 -7.5600000000000D-02    8.5087019230769D-07
1131                 -6.3000000000000D-02    8.4755965099715D-07
1132                 -5.0400000000000D-02    8.4761698717949D-07
1133                 -3.7800000000000D-02    8.5133636039886D-07
1134                 -2.5200000000000D-02    8.4844088319088D-07
1135                 -1.2600000000000D-02    8.4633689458689D-07
1136                0.    8.4644658119658D-07
1137                  1.2600000000000D-02    8.4630947293447D-07
1138                  2.5200000000000D-02    8.4838479344729D-07
1139                  3.7800000000000D-02    8.5125160256410D-07
1140                  5.0400000000000D-02    8.4750480769231D-07
1141                  6.3000000000000D-02    8.4742378917379D-07
1142                  7.5600000000000D-02    8.5071064814815D-07
1143                  8.8200000000000D-02    8.5545584045584D-07
1144                  1.0080000000000D-01    8.6368108974359D-07
1145                      ) ) ;
1146
1147  fluxrd1 = DEFI_NAPPE (
1148                      NOM_PARA : 'Z'
1149                      PROL_DROIT  :'CONSTANT'
1150                      PROL_GAUCHE :'CONSTANT'
1151                 PARA : (
1152                          -1.0080000000000D-01
1153                          -8.8200000000000D-02
1154                          -7.5600000000000D-02
1155                          -6.3000000000000D-02
1156                          -5.0400000000000D-02
1157                          -3.7800000000000D-02
1158                          -2.5200000000000D-02
1159                          -1.2600000000000D-02
1160                         0.
1161                           1.2600000000000D-02
1162                           2.5200000000000D-02
1163                           3.7800000000000D-02
1164                           5.0400000000000D-02
1165                           6.3000000000000D-02
1166                           7.5600000000000D-02
1167                           8.8200000000000D-02
1168                           1.0080000000000D-01
1169                 )
1170                 FONCTION : (
1171                       fly_1_1
1172                       fly_1_2
1173                       fly_1_3
1174                       fly_1_4
1175                       fly_1_5
1176                       fly_1_6
1177                       fly_1_7
1178                       fly_1_8
1179                       fly_1_9
1180                       fly_1_10
1181                       fly_1_11
1182                       fly_1_12
1183                       fly_1_13
1184                       fly_1_14
1185                       fly_1_15
1186                       fly_1_16
1187                       fly_1_17
1188                 )
1189                 );
1190 %
1191 % DEFINITION DES MATERIAUX ET AFFECTATIONS
1192 %
1193 %           **** DEFINITION DES MATERIAUX SANS FLUX NEUTRONIQUE ****
1194 %
1195 %
1196
1197
1198 % Famille assemblage : numero_fam
1199 %
1200 % Embout superieur
1201
1202 e_es  = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1203                             PROL_DROIT  : 'CONSTANT'
1204                             PROL_GAUCHE : 'CONSTANT'
1205                             VALE        : ( 100. 2.
1206                                             300. 3.   )
1207                           );
1208 nu_es  = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1209                             PROL_DROIT  : 'CONSTANT'
1210                             PROL_GAUCHE : 'CONSTANT'
1211                             VALE        : ( 100. 2.
1212                                             300. 3.     )
1213                           );
1214 al_es  = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1215                             PROL_DROIT  : 'CONSTANT'
1216                             PROL_GAUCHE : 'CONSTANT'
1217                             VALE        : ( 100. 2.
1218                                             300. 3.  )
1219                           );
1220  mt_es  = DEFI_MATERIAU ( ELAS_FO : (    E  : e_es
1221                                          NU : nu_es
1222                                       ALPHA : al_es
1223                                         RHO : 7800.
1224                                    )
1225                           );
1226 % Embout inferieur
1227
1228 e_ei  = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1229                             PROL_DROIT  : 'CONSTANT'
1230                             PROL_GAUCHE : 'CONSTANT'
1231                             VALE        : ( 100. 2.
1232                                             300. 3.  )
1233                           );
1234 nu_ei = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1235                             PROL_DROIT  : 'CONSTANT'
1236                             PROL_GAUCHE : 'CONSTANT'
1237                             VALE        : ( 100. 2.
1238                                             300. 3.    )
1239                           );
1240 al_ei = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1241                             PROL_DROIT  : 'CONSTANT'
1242                             PROL_GAUCHE : 'CONSTANT'
1243                             VALE        : ( 100. 2.
1244                                             300. 3. )
1245                           );
1246  mt_ei  = DEFI_MATERIAU ( ELAS_FO : (    E     : e_ei
1247                                          NU    : nu_ei
1248                                          ALPHA : al_ei
1249                                          RHO   : 12.
1250                                    )
1251                           );
1252  mt_rig  = DEFI_MATERIAU ( ELAS : (   E     : 1.
1253                                       NU    : 0.
1254                                       ALPHA : 0.
1255                                       RHO   : 0.)
1256                          );
1257 % Fonction couple en fonction de drz pour la liaison grille/crayon
1258
1259  fotrac= DEFI_FONCTION ( NOM_PARA   : 'DRZ'
1260                            VALE       : ( 1. 1.
1261                                           0. 0.
1262                                           0. 0.
1263                                           0. 0.
1264                                          )
1265                          PROL_DROIT : 'LINEAIRE'
1266                          PROL_GAUCHE: 'LINEAIRE' );
1267
1268
1269  mt_re1 = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1270                                             RELA_MZ  : fotrac
1271                                             RIGI_N_FO:frig1
1272                                      EFFO_N_INIT : -35.)
1273                            ELAS : ( E     : 1.
1274                                     NU    : 0.
1275                                     ALPHA : 0.
1276                                     RHO   : 0. )
1277                          );
1278  mt_re2 = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB      : 1.
1279                                           RELA_MZ     : fotrac
1280                                           RIGI_N_FO   : frig2
1281                                      EFFO_N_INIT : -35.)
1282                            ELAS : ( E     : 1.
1283                                     NU    : 0.
1284                                     ALPHA : 0.
1285                                     RHO   : 0. )
1286                          );
1287  mt_re3  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1288                                           RELA_MZ     : fotrac
1289                                           RIGI_N_FO   : frig3
1290                                      EFFO_N_INIT : -35.)
1291                            ELAS : ( E     : 1.
1292                                     NU    : 0.
1293                                     ALPHA : 0.
1294                                     RHO   : 0. )
1295                          );
1296  mt_re4  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1297                                           RELA_MZ     : fotrac
1298                                           RIGI_N_FO   : frig4
1299                                      EFFO_N_INIT : -35.)
1300                            ELAS : ( E     : 1.
1301                                     NU    : 0.
1302                                     ALPHA : 0.
1303                                     RHO   : 0. )
1304                          );
1305  mt_re5  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1306                                           RELA_MZ     : fotrac
1307                                           RIGI_N_FO   : frig5
1308                                      EFFO_N_INIT : -35.)
1309                            ELAS : ( E     : 1.
1310                                     NU    : 0.
1311                                     ALPHA : 0.
1312                                     RHO   : 0. )
1313                          );
1314  mt_re6  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1315                                           RELA_MZ     : fotrac
1316                                           RIGI_N_FO   : frig6
1317                                      EFFO_N_INIT :-35.)
1318                            ELAS : ( E     : 1.
1319                                     NU    : 0.
1320                                     ALPHA : 0.
1321                                     RHO   : 0. )
1322                          );
1323  mt_re7  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1324                                           RELA_MZ     : fotrac
1325                                           RIGI_N_FO   : frig7
1326                                      EFFO_N_INIT : -35.)
1327                            ELAS : ( E     : 1.
1328                                     NU    : 0.
1329                                     ALPHA : 0.
1330                                     RHO   : 0. )
1331                          );
1332  mt_re8  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1333                                           RELA_MZ     : fotrac
1334                                           RIGI_N_FO   : frig8
1335                                      EFFO_N_INIT : -35.)
1336                            ELAS : ( E     : 1.
1337                                     NU    : 0.
1338                                     ALPHA : 0.
1339                                     RHO   : 0. )
1340                          );
1341  mt_re9  = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1342                                           RELA_MZ     : fotrac
1343                                           RIGI_N_FO   :frig9
1344                                      EFFO_N_INIT : -35.)
1345                            ELAS : ( E     : 1.
1346                                     NU    : 0.
1347                                     ALPHA : 0.
1348                                     RHO   : 0. )
1349                          );
1350  mt_re10 = DEFI_MATERIAU ( DIS_CONTACT : (COULOMB     : 1.
1351                                           RELA_MZ     : fotrac
1352                                           RIGI_N_FO   : frig10
1353                                      EFFO_N_INIT : -35.)
1354                            ELAS : ( E     : 1.
1355                                     NU    : 0.
1356                                     ALPHA : 0.
1357                                     RHO   : 0. )
1358                          );
1359
1360
1361 % Systeme de maintien
1362
1363  e_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1364                               PROL_DROIT  : 'CONSTANT'
1365                               PROL_GAUCHE : 'CONSTANT'
1366                               VALE        : ( 0. 1.
1367                                               5. 2. )
1368                             );
1369
1370  nu_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1371                                PROL_DROIT  : 'CONSTANT'
1372                                PROL_GAUCHE : 'CONSTANT'
1373                                VALE        : ( 0. 1.
1374                                                5. 2.   )
1375                              );
1376
1377  al_ma    = DEFI_CONSTANTE (  VALE        :  0.3  );
1378
1379  sy_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1380                               PROL_DROIT  : 'CONSTANT'
1381                                PROL_GAUCHE : 'CONSTANT'
1382                                VALE        : ( 0. 1.
1383                                                5. 2. )
1384                               );
1385
1386  ds_ma    = DEFI_CONSTANTE (  VALE        :  0.3 );
1387
1388  np_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1389                                PROL_DROIT  : 'CONSTANT'
1390                                PROL_GAUCHE : 'CONSTANT'
1391                                VALE        : ( 0. 1.
1392                                                5. 2. )
1393                               );
1394
1395  mey_ma  = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1396                                PROL_DROIT  : 'CONSTANT'
1397                                 PROL_GAUCHE : 'CONSTANT'
1398                                 VALE        : ( 0. 1.
1399                                                 5. 2.  )
1400                               );
1401
1402  mpy_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1403                                 PROL_DROIT  : 'CONSTANT'
1404                                 PROL_GAUCHE : 'CONSTANT'
1405                                 VALE        : ( 0. 1.
1406                                                 5. 2.  )
1407                               );
1408
1409  mez_ma     = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1410                                 PROL_DROIT  : 'CONSTANT'
1411                                 PROL_GAUCHE : 'CONSTANT'
1412                                 VALE        : ( 0. 1.
1413                                                 5. 2.  )
1414                                );
1415
1416  mpz_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1417                                PROL_DROIT  : 'CONSTANT'
1418                                PROL_GAUCHE : 'CONSTANT'
1419                                VALE        : (0. 1.
1420                                               5. 2.  )
1421                              );
1422
1423  cay_ma    = DEFI_CONSTANTE (  VALE        :  0.3 );
1424
1425  cby_ma    = DEFI_CONSTANTE (  VALE        :   0.3  );
1426
1427  caz_ma    = DEFI_CONSTANTE (  VALE        :   0.3 );
1428
1429  cbz_ma    = DEFI_CONSTANTE (  VALE        :   0.3 );
1430
1431  mpx_ma    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1432                              PROL_DROIT  : 'CONSTANT'
1433                              PROL_GAUCHE : 'CONSTANT'
1434                              VALE        : ( 0. 1.
1435                                              5. 2.  )
1436                           );
1437
1438 mt_ma     = DEFI_MATERIAU ( ELAS_FO : (  E       : e_ma
1439                                             NU    : nu_ma
1440                                             ALPHA : al_ma
1441                                             RHO   : 7800. )
1442                             ECRO_LINE_FO:   ( D_SIGM_EPSI:ds_ma
1443                                               SY : sy_ma   )
1444                             VMIS_POUTRE_FO: ( NP : np_ma
1445                                               MEY: mey_ma
1446                                               MPY: mpy_ma
1447                                               CAY: cay_ma
1448                                               CBY: cby_ma
1449                                               MEZ: mez_ma
1450                                               MPZ: mpz_ma
1451                                               CAZ: caz_ma
1452                                               CBZ: cbz_ma
1453                                               MPX: mpx_ma )
1454                           );
1455
1456 %
1457 %           **** DEFINITION DES MATERIAUX AVEC FLUX NEUTRONIQUE ****
1458 %
1459 %
1460
1461 % Grille
1462
1463
1464
1465 e_gre     = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1466                              PROL_DROIT  : 'CONSTANT'
1467                              PROL_GAUCHE : 'CONSTANT'
1468                              VALE        : ( 0. 1.
1469                                              5. 2.   )
1470                           );
1471
1472 e_grm     = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1473                              PROL_DROIT  : 'CONSTANT'
1474                              PROL_GAUCHE : 'CONSTANT'
1475                              VALE        : ( 0. 1.
1476                                              5. 2.   )
1477                           );
1478
1479 nu_gre    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1480                              PROL_DROIT  : 'CONSTANT'
1481                              PROL_GAUCHE : 'CONSTANT'
1482                              VALE        : ( 0. 1.
1483                                              5. 2.       )
1484                           );
1485
1486 nu_grm    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1487                              PROL_DROIT  : 'CONSTANT'
1488                              PROL_GAUCHE : 'CONSTANT'
1489                              VALE        : ( 0. 1.
1490                                              5. 2.      )
1491                           );
1492
1493 al_gre    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1494                              PROL_DROIT  : 'CONSTANT'
1495                              PROL_GAUCHE : 'CONSTANT'
1496                              VALE        : ( 0. 1.
1497                                              5. 2.  )
1498                           );
1499
1500 al_grm    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1501                              PROL_DROIT  : 'CONSTANT'
1502                              PROL_GAUCHE : 'CONSTANT'
1503                              VALE        : ( 0. 1.
1504                                              5. 2.  )
1505                           );
1506
1507 % Tube guide
1508
1509 e_tu     = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1510                             PROL_DROIT  : 'CONSTANT'
1511                             PROL_GAUCHE : 'CONSTANT'
1512                             VALE        : ( 0. 1.
1513                                             5. 2.
1514                                           )
1515                           );
1516
1517 nu_tu    = DEFI_FONCTION ( NOM_PARA    : 'TEMP'
1518                             PROL_DROIT  : 'CONSTANT'
1519                             PROL_GAUCHE : 'CONSTANT'
1520                             VALE        : ( 0. 1.
1521                                             5. 2.     )
1522                           );
1523
1524 al_tu    = DEFI_FONCTION (  NOM_PARA    : 'TEMP'
1525                             PROL_DROIT  : 'CONSTANT'
1526                             PROL_GAUCHE : 'CONSTANT'
1527                             VALE        : ( 0. 1.
1528                                             5. 2.  )
1529                           );
1530
1531
1532
1533
1534 % Crayon
1535
1536 e_cr     = DEFI_FONCTION (  NOM_PARA    : 'TEMP'
1537                             PROL_DROIT  : 'CONSTANT'
1538                             PROL_GAUCHE : 'CONSTANT'
1539                             VALE        : ( 0. 1.
1540                                             5. 2.
1541                                           )
1542                           );
1543
1544 nu_cr    = DEFI_FONCTION (  NOM_PARA    : 'TEMP'
1545                             PROL_DROIT  : 'CONSTANT'
1546                             PROL_GAUCHE : 'CONSTANT'
1547                             VALE        : ( 0. 1.
1548                                             5. 2.    )
1549                           );
1550
1551 al_cr    = DEFI_FONCTION (  NOM_PARA    : 'TEMP'
1552                             PROL_DROIT  : 'CONSTANT'
1553                             PROL_GAUCHE : 'CONSTANT'
1554                             VALE        : ( 0. 1.
1555                                             5. 2.  )
1556                           );
1557
1558
1559 % GRANDISSEMENT DES GRILLES
1560
1561
1562
1563 % Grilles extremes 1 et 10
1564
1565 mt_gre1  = DEFI_MATERIAU ( ELAS_FO : ( E     : e_gre
1566                                         NU    : nu_gre
1567                                         ALPHA : al_gre
1568                                         RHO   : 7800. )
1569                                            GRAN_IRRA:( FLUX_L   :fluxax1
1570                                         FLUX_TN  :fluxrd1
1571                                         FONC_MULT:flux_f1
1572                                         A        : 0.
1573                                         S        : 1. )
1574                             LEMAITRE :( N        : 0.
1575                                         UN_SUR_K : 0.
1576                                         UN_SUR_M : 0.
1577                                         QSR_K    : 1.  )
1578                           );
1579
1580 % Grilles de melange
1581
1582 mt_grm1= DEFI_MATERIAU ( ELAS_FO : ( E        : e_grm
1583                                       NU       : nu_grm
1584                                       ALPHA    : al_grm
1585                                       RHO      : 7800.  )
1586                            GRAN_IRRA: ( FLUX_L  :fluxax1
1587                                        FLUX_TN  :fluxrd1
1588                                        FONC_MULT:flux_f1
1589                                        A        : 0.
1590                                        S        : 1. )
1591                            LEMAITRE : (N        : 0.
1592                                        UN_SUR_K : 0.
1593                                        UN_SUR_M : 0.
1594                                        QSR_K    : 2. )
1595                           );
1596
1597 % Definition des materiaux viscoplastiques ( tubes et crayons )
1598
1599
1600 mt_tu1  = DEFI_MATERIAU ( ELAS_FO : ( E         : e_tu
1601                                         NU       : nu_tu
1602                                         ALPHA    : al_tu
1603                                         RHO      : 7800. )
1604                            GRAN_IRRA: ( FLUX_L   : fluxax1
1605                                         FLUX_TN  : fluxrd1
1606                                         FONC_MULT: flux_f1
1607                                         A        : 0.
1608                                         S        : 1.)
1609                            LEMAITRE : ( N        : 0.
1610                                         UN_SUR_K : 0.
1611                                         UN_SUR_M : 0.
1612                                         QSR_K    : 2.)
1613                           );
1614
1615
1616  mt_cr1  = DEFI_MATERIAU (ELAS_FO : ( E        : e_cr
1617                                         NU       : nu_cr
1618                                         ALPHA    : al_cr
1619                                         RHO      :7800.  )
1620                            GRAN_IRRA: ( FLUX_L   : fluxax1
1621                                         FLUX_TN  : fluxrd1
1622                                         FONC_MULT: flux_f1
1623                                         A        : 0.
1624                                         S        : 1.)
1625                                           LEMAITRE : ( N        : 0.
1626                                         UN_SUR_K : 0.
1627                                         UN_SUR_M : 0.
1628                                         QSR_K    : 2.)
1629
1630                                      );
1631
1632 %
1633 %           AFFECTATION DES MATERIAUX A CHAUD - CYCLE
1634 %
1635 %
1636  a_mat_1  = AFFE_MATERIAU ( MAILLAGE : ma1
1637
1638   % systeme de maintien
1639                          AFFE : ( GROUP_MA : (mainti_y, mainti_z)
1640                                   MATER : mt_ma
1641                                   TEMP_REF: 100.
1642                                 )
1643            % grilles de melange
1644                          AFFE : ( GROUP_MA : gril_i
1645                                   MATER : mt_grm1
1646                                   TEMP_REF: 100.
1647                                 )
1648            % grilles de extremes
1649                          AFFE : ( GROUP_MA : gril_e
1650                                   MATER : mt_gre1
1651                                   TEMP_REF: 100.
1652                                 )
1653            % crayons combustibles
1654                          AFFE : ( GROUP_MA : cray
1655                                   MATER : mt_cr1
1656                                   TEMP_REF: 100.
1657                                 )
1658            % tubes guides
1659                          AFFE : ( GROUP_MA : guid
1660                                   MATER : mt_tu1
1661                                   TEMP_REF: 100.
1662                                 )
1663            % embout inferieur
1664                          AFFE : ( GROUP_MA : eboinf
1665                                   MATER : mt_ei
1666                                   TEMP_REF: 100.
1667                                 )
1668            %pieds
1669                          AFFE : ( GROUP_MA :(pi_ple, pi_tro)
1670                                   MATER : mt_ei
1671                                  TEMP_REF: 100.
1672                                 )
1673
1674            % embout superieur
1675                          AFFE : ( GROUP_MA : ebosup
1676                                   MATER : mt_es
1677                                   TEMP_REF: 100.
1678                                 )
1679            % liaison grille tubes guides
1680                          AFFE : ( GROUP_MA : rigid
1681                                   MATER : mt_rig
1682                                   TEMP_REF: 100.
1683                                 )
1684
1685            % liaison grille crayon
1686                          AFFE : ( GROUP_MA : ela1
1687                                   MATER : mt_re1
1688                                   TEMP_REF: 100.
1689                                 )
1690                          AFFE : ( GROUP_MA : ela2
1691                                   MATER : mt_re2
1692                                   TEMP_REF: 100.
1693                                 )
1694                          AFFE : ( GROUP_MA : ela3
1695                                   MATER : mt_re3
1696                                   TEMP_REF: 100.
1697                                 )
1698                          AFFE : ( GROUP_MA : ela4
1699                                   MATER : mt_re4
1700                                   TEMP_REF: 100.
1701                                 )
1702                          AFFE : ( GROUP_MA : ela5
1703                                   MATER : mt_re5
1704                                   TEMP_REF: 100.
1705                                 )
1706                          AFFE : ( GROUP_MA : ela6
1707                                   MATER : mt_re6
1708                                   TEMP_REF: 100.
1709                                 )
1710                          AFFE : ( GROUP_MA : ela7
1711                                   MATER : mt_re7
1712                                   TEMP_REF: 100.
1713                                 )
1714                          AFFE : ( GROUP_MA : ela8
1715                                   MATER : mt_re8
1716                                   TEMP_REF: 100.
1717                                 )
1718                          AFFE : ( GROUP_MA : ela9
1719                                   MATER : mt_re9
1720                                   TEMP_REF: 100.
1721                                 )
1722                          AFFE : ( GROUP_MA : ela10
1723                                   MATER : mt_re10
1724                                   TEMP_REF: 100.
1725                                 )
1726                      );
1727 %
1728
1729
1730  L_INST_1 = DEFI_LIST_REEL(
1731                          DEBUT: 0.0
1732                          INTERVALLE : (JUSQU_A: 0.5    NOMBRE: 1  )
1733                          INTERVALLE : (JUSQU_A: 1000.  NOMBRE: 1   )
1734                          INTERVALLE : (JUSQU_A: 2000.  NOMBRE: 1   )
1735                          INTERVALLE : (JUSQU_A: 3000.  NOMBRE: 5   )
1736                          INTERVALLE : (JUSQU_A: 4000.  NOMBRE: 1   )
1737                          INTERVALLE : (JUSQU_A: 5000.  NOMBRE: 1   )
1738                          INTERVALLE : (JUSQU_A: 6000.  NOMBRE: 1   )
1739                          INTERVALLE : (JUSQU_A: 7000.  NOMBRE: 1   )
1740                          INTERVALLE : (JUSQU_A: 8000.  NOMBRE: 1   ));
1741
1742  L_ARCH_1 = DEFI_LIST_REEL(
1743                          DEBUT: 0.0
1744                          INTERVALLE : (JUSQU_A: 0.5     NOMBRE: 1  )
1745                          INTERVALLE : (JUSQU_A: 1000.   NOMBRE: 1  )
1746                          INTERVALLE : (JUSQU_A: 2000.   NOMBRE: 1  )
1747                          INTERVALLE : (JUSQU_A: 3000.   NOMBRE: 5  )
1748                          INTERVALLE : (JUSQU_A: 4000.   NOMBRE: 1  )
1749                          INTERVALLE : (JUSQU_A: 5000.   NOMBRE: 1  )
1750                          INTERVALLE : (JUSQU_A: 6000.   NOMBRE: 1  )
1751                          INTERVALLE : (JUSQU_A: 7000.   NOMBRE: 1  )
1752                          INTERVALLE : (JUSQU_A: 8000.   NOMBRE: 1  ));
1753
1754 %
1755 % CALCUL PHASE I A V DU CYCLE 1
1756 %
1757
1758   RESU = STAT_NON_LINE( MODELE     : mo1
1759                         CHAM_MATER : a_mat_1
1760                         CARA_ELEM  : cara1
1761                         EXCIT :( CHARGE : arch_1   FONC_MULT: arch_f1 )
1762                         EXCIT :( CHARGE : foarch_1 FONC_MULT: arch_f1 )
1763                         EXCIT :( CHARGE : hydr_1   FONC_MULT: hydr_f1 )
1764                         EXCIT :( CHARGE : fohydr_1 FONC_MULT: hydr_f1 )
1765                         EXCIT :( CHARGE : cl_di_1  FONC_MULT: cldi_f1 )
1766                         EXCIT :( CHARGE : cl_per_1 )
1767                         EXCIT :( CHARGE : cl_ctc   )
1768                         EXCIT :( CHARGE : pesant_1 )
1769                         EXCIT :( CHARGE : therm_1  )
1770                         COMP_INCR :(RELATION :'VMIS_POU_LINE'
1771                                     GROUP_MA : maintien)
1772                         COMP_INCR :(RELATION :'ASSE_COMBU'
1773                                     GROUP_MA : (crayons tub_guid     ))
1774                         COMP_INCR :( RELATION :'DIS_CONTACT'
1775                                      GROUP_MA : li_gr_cr)
1776                         COMP_INCR :(RELATION :'ELAS'
1777                                     COQUE_NCOU : 1
1778                                     GROUP_MA :(embo_inf embo_sup pieds
1779                                                grilles li_gr_gu ))
1780
1781                         INCREMENT :( LIST_INST : l_inst_1
1782                                      NUME_INST_FIN  :  20
1783                                    )
1784
1785                         NEWTON : ( MATRICE       : 'TANGENTE'
1786                                    REAC_ITER    : 1          )
1787                         RECH_LINEAIRE : (ITER_LINE_MAXI: 3)
1788                         CONVERGENCE : (RESI_GLOB_RELA: 1.E-4
1789                                        ITER_GLOB_MAXI: 50
1790                                        ITER_INTE_MAXI: 500
1791                                        ARRET         : 'OUI' )
1792                         SOLVEUR : ( METHODE : 'MULT_FRONT'
1793                                     RENUM   : 'MD' )
1794                         PARM_THETA : 0.5
1795                         ARCHIVAGE : ( LIST_INST : l_arch_1)
1796                       );
1797
1798 fin();