Exemples/ex18/totalmod.comm

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