Exemples/ex22/toto.comm

   1
   2
   3 DEBUT();
   4 # LECTURE DES PARAMETRES
   5 #
   6 # LECTURE DES TEMPS
   7 #CYCLE 1
   8 # LECTURE DU MAILLAGE ET DEFINITIONS DES PROPRIETES ELEMENTS FINIS
   9
  10
  11 PRE_GIBI();
  12 #
  13
  14
  15 MA1=LIRE_MAILLAGE();
  16 #
  17
  18
  19 MA1=DEFI_GROUP(reuse =MA1,
  20                CREA_GROUP_NO=_F(GROUP_MA=('FIXATI','APPUI','LIAS_1','LIAS_2','EBGUIDE','MAINTIEN','MAIN_SUP','EMBO_SUP','TUB_GUID',
  21                                            'GRILLES','EMBO_INF','PIEDS','CRAYONS',),),
  22                MAILLAGE=MA1,);
  23 #
  24 # AFFECTATION DU MODELE
  25 #
  26
  27
  28 MO1=AFFE_MODELE(AFFE=(_F(GROUP_MA='CRAYONS',
  29                          PHENOMENE='MECANIQUE',
  30                          MODELISATION='POU_D_E',),
  31                       _F(GROUP_MA='TUB_GUID',
  32                          PHENOMENE='MECANIQUE',
  33                          MODELISATION='POU_D_E',),
  34                       _F(GROUP_MA=('EMBO_SUP','EMBO_INF',),
  35                          PHENOMENE='MECANIQUE',
  36                          MODELISATION='DKT',),
  37                       _F(GROUP_MA='PIEDS',
  38                          PHENOMENE='MECANIQUE',
  39                          MODELISATION='POU_D_E',),
  40                       _F(GROUP_MA='GRILLES',
  41                          PHENOMENE='MECANIQUE',
  42                          MODELISATION='POU_D_E',),
  43                       _F(GROUP_MA='LI_GR_CR',
  44                          PHENOMENE='MECANIQUE',
  45                          MODELISATION='DIS_TR',),
  46                       _F(GROUP_MA='LI_GR_GU',
  47                          PHENOMENE='MECANIQUE',
  48                          MODELISATION='DIS_TR',),
  49                       _F(GROUP_MA='MAINTIEN',
  50                          PHENOMENE='MECANIQUE',
  51                          MODELISATION='POU_D_E',),),
  52                 MAILLAGE=MA1,);
  53 #
  54 # AFFECTATION DES CARACTERISTIQUES GEOMETRIQUES ELEMENTAIRES
  55 #
  56
  57
  58 CARA1=AFFE_CARA_ELEM(MODELE=MO1,
  59                      COQUE=(_F(ANGL_REP=90.0,
  60                                GROUP_MA='EBOINF',
  61                                EPAIS=2.0,),
  62                             _F(ANGL_REP=90.0,
  63                                GROUP_MA='EBOSUP',
  64                                EPAIS=3.0,),),
  65                      DISCRET=(_F(VALE=(1.0,2.0,3.0,4.0,5.0,6.0,),
  66                                  REPERE='LOCAL',
  67                                  CARA='K_TR_D_L',
  68                                  GROUP_MA='ELASTI',),
  69                               _F(VALE=(1.0,2.0,3.0,4.0,5.0,6.0,),
  70                                  REPERE='LOCAL',
  71                                  CARA='K_TR_D_L',
  72                                  GROUP_MA='RIGID',),),
  73                      ORIENTATION=(_F(VALE=(1.0,0.0,0.0,),
  74                                      CARA='VECT_Y',
  75                                      GROUP_MA=('GRIL_I','GRIL_E',),),
  76                                   _F(VALE=(0.0,1.0,0.0,),
  77                                      CARA='VECT_Y',
  78                                      GROUP_MA='MAINTI_Y',),
  79                                   _F(VALE=(0.0,0.0,1.0,),
  80                                      CARA='VECT_Y',
  81                                      GROUP_MA='MAINTI_Z',),
  82                                   _F(VALE=(0.0,1.0,0.0,),
  83                                      CARA='VECT_Y',
  84                                      GROUP_MA=('PI_PLE','PI_TRO',),),
  85                                   _F(VALE=(1.0,0.0,0.0,),
  86                                      CARA='VECT_Y',
  87                                      GROUP_MA=('ELASTI','RIGID',),),),
  88                      POUTRE=(_F(SECTION='CERCLE',
  89                                 GROUP_MA='CRAY',
  90                                 VALE=(3.0,3.0,),
  91                                 CARA=('R','EP',),),
  92                              _F(SECTION='CERCLE',
  93                                 GROUP_MA='GUID',
  94                                 VALE=(3.0,3.0,),
  95                                 CARA=('R','EP',),),
  96                              _F(SECTION='CERCLE',
  97                                 GROUP_MA='RETRE',
  98                                 VALE=(3.0,3.0,),
  99                                 CARA=('R','EP',),),
 100                              _F(SECTION='RECTANGLE',
 101                                 GROUP_MA='GRIL_I',
 102                                 VALE=(3.0,3.0,),
 103                                 CARA=('HY','HZ',),),
 104                              _F(SECTION='RECTANGLE',
 105                                 GROUP_MA='GRIL_E',
 106                                 VALE=(3.0,3.0,),
 107                                 CARA=('HY','HZ',),),
 108                              _F(SECTION='RECTANGLE',
 109                                 GROUP_MA=('MAINTI_Y','MAINTI_Z',),
 110                                 VALE=(3.0,3.0,),
 111                                 CARA=('HY','HZ',),),
 112                              _F(SECTION='GENERALE',
 113                                 GROUP_MA='PI_PLE',
 114                                 VALE=(1.0,2.0,3.0,4.0,),
 115                                 CARA=('A','IZ','IY','JX',),),
 116                              _F(SECTION='GENERALE',
 117                                 GROUP_MA='PI_TRO',
 118                                 VALE=(1.0,2.0,3.0,4.0,),
 119                                 CARA=('A','IZ','IY','JX',),),),);
 120 #
 121 # DEFINITION DES CHARGES ET CONDITIONS LIMITES
 122 # DEFINITION DES TEMPERATURES NODALES EVOLUTIVES
 123 #
 124
 125
 126 F_TP1_1=DEFI_FONCTION(NOM_RESU='TEMP',
 127                       PROL_GAUCHE='CONSTANT',
 128                       NOM_PARA='X',
 129                       VALE=(0.0,1.0,5.0,4.0,),
 130                       PROL_DROIT='CONSTANT',);
 131
 132 F_TP2_1=DEFI_FONCTION(NOM_RESU='TEMP',
 133                       PROL_GAUCHE='CONSTANT',
 134                       NOM_PARA='X',
 135                       VALE=(0.0,1.0,5.0,4.0,),
 136                       PROL_DROIT='CONSTANT',);
 137
 138 F_TP3_1=DEFI_FONCTION(NOM_RESU='TEMP',
 139                       PROL_GAUCHE='CONSTANT',
 140                       NOM_PARA='X',
 141                       VALE=(0.0,1.0,5.0,4.0,),
 142                       PROL_DROIT='CONSTANT',);
 143
 144 CHTEM1_1=AFFE_CHAM_NO(GRANDEUR='TEMP_F',
 145                       AFFE=(_F(FONCTION=F_TP1_1,
 146                                GROUP_NO=('MAINTIEN','PIEDS','TUB_GUID','GRILLES','CRAYONS',),
 147                                NOM_CMP='TEMP',),
 148                             _F(FONCTION=(F_TP1_1,F_TP1_1,F_TP1_1,),
 149                                GROUP_NO=('EMBO_SUP','EMBO_INF',),
 150                                NOM_CMP=('TEMP_INF','TEMP','TEMP_SUP',),),),
 151                       MAILLAGE=MA1,);
 152
 153 CHTEM2_1=AFFE_CHAM_NO(GRANDEUR='TEMP_F',
 154                       AFFE=(_F(FONCTION=F_TP2_1,
 155                                GROUP_NO=('MAINTIEN','PIEDS','TUB_GUID','GRILLES','CRAYONS',),
 156                                NOM_CMP='TEMP',),
 157                             _F(FONCTION=(F_TP2_1,F_TP2_1,F_TP2_1,),
 158                                GROUP_NO=('EMBO_SUP','EMBO_INF',),
 159                                NOM_CMP=('TEMP_INF','TEMP','TEMP_SUP',),),),
 160                       MAILLAGE=MA1,);
 161
 162 CHTEM3_1=AFFE_CHAM_NO(GRANDEUR='TEMP_F',
 163                       AFFE=(_F(FONCTION=F_TP3_1,
 164                                GROUP_NO=('MAINTIEN','PIEDS','TUB_GUID','GRILLES',),
 165                                NOM_CMP='TEMP',),
 166                             _F(FONCTION=(F_TP3_1,F_TP3_1,F_TP3_1,),
 167                                GROUP_NO=('EMBO_SUP','EMBO_INF',),
 168                                NOM_CMP=('TEMP_INF','TEMP','TEMP_SUP',),),
 169                             _F(FONCTION=F_TP3_1,
 170                                GROUP_NO='CRAYONS',
 171                                NOM_CMP='TEMP',),),
 172                       MAILLAGE=MA1,);
 173
 174 CHTH_1=CREA_RESU(NOM_CHAM='TEMP',
 175                  CHAM_GD=(_F(INST=0.0,
 176                              CHAM_NO=CHTEM1_1,),
 177                           _F(INST=1000.0,
 178                              CHAM_NO=CHTEM1_1,),
 179                           _F(INST=2000.0,
 180                              CHAM_NO=CHTEM2_1,),
 181                           _F(INST=3000.0,
 182                              CHAM_NO=CHTEM3_1,),
 183                           _F(INST=4000.0,
 184                              CHAM_NO=CHTEM3_1,),
 185                           _F(INST=5000.0,
 186                              CHAM_NO=CHTEM2_1,),
 187                           _F(INST=6000.0,
 188                              CHAM_NO=CHTEM1_1,),
 189                           _F(INST=7000.0,
 190                              CHAM_NO=CHTEM1_1,),),
 191                  TYPE_RESU='EVOL_THER',);
 192 # CHARGEMENT: CHARGES THERMO-MECANIQUES
 193 # ----------------------------------------
 194
 195
 196 THERM_1=AFFE_CHAR_MECA(MODELE=MO1,
 197                        TEMP_CALCULEE=CHTH_1,);
 198 # DEFINITION DES CONDITIONS LIMITES
 199 #            DES CONDITIONS DE CONTACT
 200 #
 201 # CHARGEMENT: PESANTEUR
 202 #
 203
 204
 205 PESANT_1=AFFE_CHAR_MECA(MODELE=MO1,
 206                         PESANTEUR=(9.8100000000000005,-1.0,0.0,0.0,),);
 207 #
 208 # CONDITIONS AUX LIMITES
 209 #
 210
 211
 212 CL_PER_1=AFFE_CHAR_MECA(MODELE=MO1,
 213                         DDL_IMPO=(_F(DZ=0.0,
 214                                      DX=0.0,
 215                                      DY=0.0,
 216                                      GROUP_NO='FIXATI',
 217                                      DRZ=0.0,
 218                                      DRX=0.0,
 219                                      DRY=0.0,),
 220                                   _F(DX=0.0,
 221                                      GROUP_NO='APPUI',),
 222                                   _F(DZ=0.0,
 223                                      DY=0.0,
 224                                      GROUP_NO='EBGUIDE',
 225                                      DRZ=0.0,
 226                                      DRX=0.0,
 227                                      DRY=0.0,),
 228                                   _F(GROUP_NO=('LIAS_1','LIAS_2',),
 229                                      DRX=0.0,),),
 230                         LIAISON_GROUP=(_F(COEF_MULT_1=1.0,
 231                                           DDL_1='DX',
 232                                           COEF_MULT_2=-1.0,
 233                                           COEF_IMPO=0.0,
 234                                           GROUP_NO_1='LIAS_1',
 235                                           GROUP_NO_2='LIAS_2',
 236                                           DDL_2='DX',),
 237                                        _F(COEF_MULT_1=1.0,
 238                                           DDL_1='DY',
 239                                           COEF_MULT_2=-1.0,
 240                                           COEF_IMPO=0.0,
 241                                           GROUP_NO_1='LIAS_1',
 242                                           GROUP_NO_2='LIAS_2',
 243                                           DDL_2='DY',),
 244                                        _F(COEF_MULT_1=1.0,
 245                                           DDL_1='DZ',
 246                                           COEF_MULT_2=-1.0,
 247                                           COEF_IMPO=0.0,
 248                                           GROUP_NO_1='LIAS_1',
 249                                           GROUP_NO_2='LIAS_2',
 250                                           DDL_2='DZ',),
 251                                        _F(COEF_MULT_1=1.0,
 252                                           DDL_1='DRY',
 253                                           COEF_MULT_2=-1.0,
 254                                           COEF_IMPO=0.0,
 255                                           GROUP_NO_1='LIAS_1',
 256                                           GROUP_NO_2='LIAS_2',
 257                                           DDL_2='DRY',),
 258                                        _F(COEF_MULT_1=1.0,
 259                                           DDL_1='DRZ',
 260                                           COEF_MULT_2=-1.0,
 261                                           COEF_IMPO=0.0,
 262                                           GROUP_NO_1='LIAS_1',
 263                                           GROUP_NO_2='LIAS_2',
 264                                           DDL_2='DRZ',),),);
 265 #
 266
 267
 268 CL_CTC=AFFE_CHAR_MECA(INFO=1,
 269                       MODELE=MO1,
 270                       CONTACT=_F(GROUP_MA_1='EBOINF',
 271                                  VECT_NORM_2=(-1.0,0.0,0.0,),
 272                                  APPARIEMENT='NODAL',
 273                                  GROUP_MA_2='MCRBAS',),);
 274 #
 275 # DEPLACEMENT DU SYSTEME DE MAINTIEN
 276 #
 277
 278
 279 CLDI_F1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 280                       NOM_PARA='INST',
 281                       VALE=
 282                       (0.0,0.0,1000.0,0.10000000000000001,2000.0,0.20000000000000001,3000.0,
 283                              0.20000000000000001,4000.0,0.20000000000000001,5000.0,0.20000000000000001,6000.0,
 284                              0.10000000000000001,7000.0,0.10000000000000001,),
 285                       PROL_DROIT='CONSTANT',);
 286
 287 CL_DI_1=AFFE_CHAR_MECA(MODELE=MO1,
 288                        DDL_IMPO=_F(DX=-1.0,
 289                                    GROUP_NO='MAIN_SUP',),);
 290 # DEFINITION DES FORCES EXTERIEURES
 291 #
 292 # FONCTION MULTIPLICATIVE DE LA FORCE HYDRODYNAMIQUE NORMEE
 293 #
 294
 295
 296 HYDR_F1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 297                       NOM_PARA='INST',
 298                       VALE=
 299                       (0.0,0.0,1000.0,0.10000000000000001,2000.0,0.20000000000000001,3000.0,
 300                              0.20000000000000001,4000.0,0.20000000000000001,5000.0,0.20000000000000001,6000.0,
 301                              0.10000000000000001,7000.0,0.10000000000000001,),
 302                       PROL_DROIT='CONSTANT',);
 303 #
 304 # PROFIL DES FORCES HYDRODYNAMIQUES TUBES GUIDES A FROID
 305 #
 306
 307
 308 HYFRTU_1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 309                        NOM_PARA='X',
 310                        VALE=(0.0,1.0,5.0,2.0,),
 311                        PROL_DROIT='CONSTANT',);
 312 #
 313 # PROFIL DES FORCES HYDRODYNAMIQUES CRAYON A FROID
 314 #
 315
 316
 317 HYFRCR_1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 318                        NOM_PARA='X',
 319                        VALE=(0.0,1.0,5.0,2.0,),
 320                        PROL_DROIT='CONSTANT',);
 321 #
 322 # CHARGEMENT : FORCES HYDRODYNAMIQUES
 323 #
 324
 325
 326 HYDR_1=AFFE_CHAR_MECA(FORCE_POUTRE=(_F(GROUP_MA='GRIL10',
 327                                        FX=10.0,),
 328                                     _F(GROUP_MA='GRIL1',
 329                                        FX=10.0,),
 330                                     _F(GROUP_MA='GRIL2',
 331                                        FX=10.0,),
 332                                     _F(GROUP_MA='GRIL3',
 333                                        FX=10.0,),
 334                                     _F(GROUP_MA='GRIL4',
 335                                        FX=10.0,),
 336                                     _F(GROUP_MA='GRIL5',
 337                                        FX=10.0,),
 338                                     _F(GROUP_MA='GRIL6',
 339                                        FX=10.0,),
 340                                     _F(GROUP_MA='GRIL7',
 341                                        FX=10.0,),
 342                                     _F(GROUP_MA='GRIL8',
 343                                        FX=10.0,),
 344                                     _F(GROUP_MA='GRIL9',
 345                                        FX=10.0,),),
 346                       MODELE=MO1,
 347                       FORCE_COQUE=(_F(FX=10.0,
 348                                       GROUP_MA='EBOSUP',),
 349                                    _F(FX=10.0,
 350                                       GROUP_MA='EBOINF',),),);
 351
 352 TOTO_1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 353                      NOM_PARA='X',
 354                      VALE=(0.0,10.0,5.0,10.0,),
 355                      PROL_DROIT='CONSTANT',);
 356 #
 357
 358
 359 FOHYDR_1=AFFE_CHAR_MECA_F(FORCE_POUTRE=(_F(GROUP_MA='GUID',
 360                                            FX=TOTO_1,),
 361                                         _F(GROUP_MA='CRAY',
 362                                            FX=TOTO_1,),),
 363                           MODELE=MO1,);
 364 #
 365 # PROFIL DES FORCES D'ARCHIMEDE TUBE GUIDE
 366 #
 367
 368
 369 ARFRTU_1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 370                        NOM_PARA='X',
 371                        VALE=(0.0,10.0,5.0,10.0,),
 372                        PROL_DROIT='CONSTANT',);
 373 #
 374 # PROFIL DES FORCES D'ARCHIMEDE CRAYON
 375 #
 376
 377
 378 ARFRCR_1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 379                        NOM_PARA='X',
 380                        VALE=(0.0,10.0,5.0,10.0,),
 381                        PROL_DROIT='CONSTANT',);
 382 #
 383 # FONCTION MULTIPLICATIVE DE LA FORCE D'ARCHIMEDE
 384 #
 385
 386
 387 ARCH_F1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 388                       NOM_PARA='INST',
 389                       VALE=
 390                       (0.0,0.0,1000.0,0.10000000000000001,2000.0,0.20000000000000001,3000.0,
 391                              0.20000000000000001,4000.0,0.20000000000000001,5000.0,0.20000000000000001,6000.0,
 392                              0.10000000000000001,7000.0,0.10000000000000001,),
 393                       PROL_DROIT='CONSTANT',);
 394 #
 395 #
 396 # CHARGEMENT : FORCES D'ARCHIMEDE
 397 # ----------------------------------------
 398 #
 399
 400
 401 ARCH_1=AFFE_CHAR_MECA(FORCE_POUTRE=(_F(GROUP_MA='MAINTIEN',
 402                                        FX=10.0,),
 403                                     _F(GROUP_MA='GRIL10',
 404                                        FX=10.0,),
 405                                     _F(GROUP_MA='GRIL1',
 406                                        FX=10.0,),
 407                                     _F(GROUP_MA='GRIL2',
 408                                        FX=10.0,),
 409                                     _F(GROUP_MA='GRIL3',
 410                                        FX=10.0,),
 411                                     _F(GROUP_MA='GRIL4',
 412                                        FX=10.0,),
 413                                     _F(GROUP_MA='GRIL5',
 414                                        FX=10.0,),
 415                                     _F(GROUP_MA='GRIL6',
 416                                        FX=10.0,),
 417                                     _F(GROUP_MA='GRIL7',
 418                                        FX=10.0,),
 419                                     _F(GROUP_MA='GRIL8',
 420                                        FX=10.0,),
 421                                     _F(GROUP_MA='GRIL9',
 422                                        FX=10.0,),),
 423                       MODELE=MO1,
 424                       FORCE_COQUE=(_F(FX=10.0,
 425                                       GROUP_MA='EBOSUP',),
 426                                    _F(FX=10.0,
 427                                       GROUP_MA='EBOINF',),),);
 428
 429 FOARCH_1=AFFE_CHAR_MECA_F(FORCE_POUTRE=(_F(GROUP_MA='GUID',
 430                                            FX=ARCH_F1,),
 431                                         _F(GROUP_MA='CRAY',
 432                                            FX=ARCH_F1,),),
 433                           MODELE=MO1,);
 434 #
 435 #
 436 # FORCE DE SERRAGE DE LA LIAISON GRILLE - CRAYON
 437 #
 438 #
 439 #
 440 #
 441 #
 442 # DEFINITION DE LA LISTE D'INSTANTS
 443 #
 444
 445
 446 L_FRIG=DEFI_LIST_REEL(INTERVALLE=_F(JUSQU_A=1000.0,
 447                                     NOMBRE=100,),
 448                       DEBUT=0.0,);
 449 #
 450 # INRTERPRETATION DES FONCTIONS
 451 #
 452
 453
 454 F1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 455                  NOM_PARA='INST',
 456                  VALE=(0.0,0.0,43200.0,0.0,),
 457                  PROL_DROIT='CONSTANT',);
 458
 459 F2=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 460                  NOM_PARA='INST',
 461                  VALE=(0.0,0.0,43200.0,0.0,),
 462                  PROL_DROIT='CONSTANT',);
 463
 464 F3=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 465                  NOM_PARA='INST',
 466                  VALE=(0.0,0.0,43200.0,0.0,),
 467                  PROL_DROIT='CONSTANT',);
 468
 469 F4=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 470                  NOM_PARA='INST',
 471                  VALE=(0.0,0.0,43200.0,0.0,),
 472                  PROL_DROIT='CONSTANT',);
 473
 474 F5=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 475                  NOM_PARA='INST',
 476                  VALE=(0.0,0.0,43200.0,0.0,),
 477                  PROL_DROIT='CONSTANT',);
 478
 479 F6=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 480                  NOM_PARA='INST',
 481                  VALE=(0.0,0.0,43200.0,0.0,),
 482                  PROL_DROIT='CONSTANT',);
 483
 484 F7=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 485                  NOM_PARA='INST',
 486                  VALE=(0.0,0.0,43200.0,0.0,),
 487                  PROL_DROIT='CONSTANT',);
 488
 489 F8=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 490                  NOM_PARA='INST',
 491                  VALE=(0.0,0.0,43200.0,0.0,),
 492                  PROL_DROIT='CONSTANT',);
 493
 494 F9=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 495                  NOM_PARA='INST',
 496                  VALE=(0.0,0.0,43200.0,0.0,),
 497                  PROL_DROIT='CONSTANT',);
 498
 499 F10=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 500                   NOM_PARA='INST',
 501                   VALE=(0.0,0.0,43200.0,0.0,),
 502                   PROL_DROIT='CONSTANT',);
 503
 504 FRIG1=CALC_FONC_INTERP(FONCTION=F1,
 505                        PROL_DROIT='CONSTANT',
 506                        INTERPOL='LIN',
 507                        PROL_GAUCHE='CONSTANT',
 508                        LIST_PARA=L_FRIG,);
 509
 510 FRIG2=CALC_FONC_INTERP(FONCTION=F2,
 511                        PROL_DROIT='CONSTANT',
 512                        INTERPOL='LIN',
 513                        PROL_GAUCHE='CONSTANT',
 514                        LIST_PARA=L_FRIG,);
 515
 516 FRIG3=CALC_FONC_INTERP(FONCTION=F3,
 517                        PROL_DROIT='CONSTANT',
 518                        INTERPOL='LIN',
 519                        PROL_GAUCHE='CONSTANT',
 520                        LIST_PARA=L_FRIG,);
 521
 522 FRIG4=CALC_FONC_INTERP(FONCTION=F4,
 523                        PROL_DROIT='CONSTANT',
 524                        INTERPOL='LIN',
 525                        PROL_GAUCHE='CONSTANT',
 526                        LIST_PARA=L_FRIG,);
 527
 528 FRIG5=CALC_FONC_INTERP(FONCTION=F5,
 529                        PROL_DROIT='CONSTANT',
 530                        INTERPOL='LIN',
 531                        PROL_GAUCHE='CONSTANT',
 532                        LIST_PARA=L_FRIG,);
 533
 534 FRIG6=CALC_FONC_INTERP(FONCTION=F6,
 535                        PROL_DROIT='CONSTANT',
 536                        INTERPOL='LIN',
 537                        PROL_GAUCHE='CONSTANT',
 538                        LIST_PARA=L_FRIG,);
 539
 540 FRIG7=CALC_FONC_INTERP(FONCTION=F7,
 541                        PROL_DROIT='CONSTANT',
 542                        INTERPOL='LIN',
 543                        PROL_GAUCHE='CONSTANT',
 544                        LIST_PARA=L_FRIG,);
 545
 546 FRIG8=CALC_FONC_INTERP(FONCTION=F8,
 547                        PROL_DROIT='CONSTANT',
 548                        INTERPOL='LIN',
 549                        PROL_GAUCHE='CONSTANT',
 550                        LIST_PARA=L_FRIG,);
 551
 552 FRIG9=CALC_FONC_INTERP(FONCTION=F9,
 553                        PROL_DROIT='CONSTANT',
 554                        INTERPOL='LIN',
 555                        PROL_GAUCHE='CONSTANT',
 556                        LIST_PARA=L_FRIG,);
 557
 558 FRIG10=CALC_FONC_INTERP(FONCTION=F10,
 559                         PROL_DROIT='CONSTANT',
 560                         INTERPOL='LIN',
 561                         PROL_GAUCHE='CONSTANT',
 562                         LIST_PARA=L_FRIG,);
 563 #
 564 # DEFINITION DU FLUX NEUTRONIQUE RADIAL
 565
 566
 567 FLUX_F1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 568                       NOM_PARA='INST',
 569                       VALE=
 570                       (0.0,0.0,43200.0,0.0,974520.0,453280.0,2638800.0,2077850.0,4087800.0,3512790.0,
 571                              5428440.0,4856770.0,7245000.0,6676270.0,8853480.0,8279770.0,10460900.0,9872610.0,12202900.0,
 572                              11586900.0,13944600.0,13295200.0,14143700.0,13490200.0,15563200.0,14875300.0,17236100.0,
 573                              16506900.0,17948200.0,17199200.0,18190100.0,17433400.0,18193700.0,17435200.0,18197300.0,
 574                              17435200.0,18200900.0,17435200.0,),
 575                       PROL_DROIT='CONSTANT',);
 576 # DEFINITION DU FLUX AXIAL POUR LE CYCLE 1
 577
 578
 579 FLUXAX1=LIRE_FONCTION(PROL_GAUCHE='CONSTANT',
 580                       UNITE=60,
 581                       PROL_DROIT='CONSTANT',);
 582
 583 FLY_1_1=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 584                       NOM_PARA='Y',
 585                       VALE=
 586                       (-0.1008,3.8032834757835002e-07,-0.088200000000000001,3.8083689458688999e-07,
 587                              -0.075600000000000001,3.8296082621083001e-07,-0.063,3.858400997151e-07,-0.0504,
 588                              3.8946100427349999e-07,-0.0378,3.9259455128205e-07,-0.0252,3.9258457977207998e-07,-0.0126,
 589                              3.9294106125356001e-07,0.0,3.9399554843304999e-07,0.0126,3.9293233618234e-07,0.0252,
 590                              3.9256712962963002e-07,0.0378,3.9256712962963002e-07,0.0504,3.8942610398860001e-07,0.063,
 591                              3.8579522792023001e-07,0.075600000000000001,3.8290972222222002e-07,0.088200000000000001,
 592                              3.8077706552707e-07,0.1008,3.8026103988603998e-07,),
 593                       PROL_DROIT='CONSTANT',);
 594
 595 FLY_1_2=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 596                       NOM_PARA='Y',
 597                       VALE=
 598                       (-0.1008,4.3595566239316001e-07,-0.088200000000000001,4.3866542022792e-07,
 599                              -0.075600000000000001,4.4331463675214e-07,-0.063,4.4911556267806002e-07,-0.0504,
 600                              4.5750658831909e-07,-0.0378,4.6647222222221997e-07,-0.0252,4.6002938034188002e-07,-0.0126,
 601                              4.6004309116809001e-07,0.0,4.6686235754985997e-07,0.0126,4.6003311965811999e-07,0.0252,
 602                              4.6000943732193999e-07,0.0378,4.6643981481481001e-07,0.0504,4.5746420940171002e-07,0.063,
 603                              4.4906445868945999e-07,0.075600000000000001,4.4325356125356e-07,0.088200000000000001,
 604                              4.3859686609687001e-07,0.1008,4.3587838319088001e-07,),
 605                       PROL_DROIT='CONSTANT',);
 606
 607 FLY_1_3=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 608                       NOM_PARA='Y',
 609                       VALE=
 610                       (-0.1008,4.9431891025640996e-07,-0.088200000000000001,5.0005128205127998e-07,
 611                              -0.075600000000000001,5.1337197293447004e-07,-0.063,5.2915313390313001e-07,-0.0504,
 612                              5.3788069800570004e-07,-0.0378,5.3504629629630004e-07,-0.0252,5.3221189458688998e-07,-0.0126,
 613                              5.3136431623932e-07,0.0,5.3135870726495999e-07,0.0126,5.3135309829059998e-07,0.0252,
 614                              5.3218696581197002e-07,0.0378,5.3500827991452996e-07,0.0504,5.3782959401709e-07,0.063,
 615                              5.2909205840456002e-07,0.075600000000000001,5.1330217236466999e-07,0.088200000000000001,
 616                              4.9997150997151001e-07,0.1008,4.9423165954416005e-07,),
 617                       PROL_DROIT='CONSTANT',);
 618
 619 FLY_1_4=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 620                       NOM_PARA='Y',
 621                       VALE=
 622                       (-0.1008,5.5489957264956999e-07,-0.088200000000000001,5.6477884615384998e-07,
 623                              -0.075600000000000001,5.9007781339030998e-07,-0.063,5.9726415598291003e-07,-0.0504,
 624                              6.0445049857550001e-07,-0.0378,5.9840776353275996e-07,-0.0252,5.8440651709402002e-07,-0.0126,
 625                              5.8263906695157003e-07,0.0,5.9153240740741004e-07,0.0126,5.8262660256409999e-07,0.0252,
 626                              5.8438034188033999e-07,0.0378,5.9836663105413e-07,0.0504,6.0439316239316004e-07,0.063,
 627                              5.9719497863248e-07,0.075600000000000001,5.8999679487179001e-07,0.088200000000000001,
 628                              5.6469159544159996e-07,0.1008,5.5480110398860005e-07,),
 629                       PROL_DROIT='CONSTANT',);
 630
 631 FLY_1_5=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 632                       NOM_PARA='Y',
 633                       VALE=
 634                       (-0.1008,6.1267200854701004e-07,-0.088200000000000001,6.2962357549858005e-07,
 635                              -0.075600000000000001,6.5655413105413e-07,-0.063,6.6164583333332999e-07,-0.0504,
 636                              6.5190242165241998e-07,-0.0378,6.5212678062678005e-07,-0.0252,6.3746616809116999e-07,-0.0126,
 637                              6.3547435897436004e-07,0.0,6.4509686609687001e-07,0.0126,6.3546064814814995e-07,0.0252,
 638                              6.3743625356125e-07,0.0378,6.5207941595442004e-07,0.0504,6.5183885327634997e-07,0.063,
 639                              6.6156980056979999e-07,0.075600000000000001,6.5646563390312997e-07,0.088200000000000001,
 640                              6.2952386039885999e-07,0.1008,6.1256356837606998e-07,),
 641                       PROL_DROIT='CONSTANT',);
 642
 643 FLY_1_6=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 644                       NOM_PARA='Y',
 645                       VALE=
 646                       (-0.1008,6.5755128205127996e-07,-0.088200000000000001,6.8336253561253996e-07,
 647                              -0.075600000000000001,6.9012072649573003e-07,-0.063,6.9687891737892e-07,-0.0504,
 648                              6.9362571225071001e-07,-0.0378,6.9074768518518997e-07,-0.0252,6.8786965811965998e-07,-0.0126,
 649                              6.8586413817664004e-07,0.0,6.8585603632478996e-07,0.0126,6.8584793447293003e-07,0.0252,
 650                              6.8783600427350001e-07,0.0378,6.9069720441594996e-07,0.0504,6.9355840455840002e-07,0.063,
 651                              6.9679789886040002e-07,0.075600000000000001,6.9002412749288002e-07,0.088200000000000001,
 652                              6.8325035612536003e-07,0.1008,6.5743162393161998e-07,),
 653                       PROL_DROIT='CONSTANT',);
 654
 655 FLY_1_7=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 656                       NOM_PARA='Y',
 657                       VALE=
 658                       (-0.1008,6.9684526353276003e-07,-0.088200000000000001,7.1406606125356005e-07,
 659                              -0.075600000000000001,7.3236502849003e-07,-0.063,7.2076442307691997e-07,-0.0504,
 660                              7.1793500712250995e-07,-0.0378,7.2835648148147998e-07,-0.0252,7.1444747150996998e-07,-0.0126,
 661                              7.1130893874644005e-07,0.0,7.2107727920227996e-07,0.0126,7.1129398148148005e-07,0.0252,
 662                              7.1441257122506995e-07,0.0378,7.2830288461537998e-07,0.0504,7.1786396011395999e-07,0.063,
 663                              7.2067841880341996e-07,0.075600000000000001,7.3226157407406997e-07,0.088200000000000001,
 664                              7.1395138888889005e-07,0.1008,6.9671812678063005e-07,),
 665                       PROL_DROIT='CONSTANT',);
 666
 667 FLY_1_8=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 668                       NOM_PARA='Y',
 669                       VALE=
 670                       (-0.1008,7.3679362535612996e-07,-0.088200000000000001,7.5423130341880005e-07,
 671                              -0.075600000000000001,7.7215633903134002e-07,-0.063,7.5865740740741005e-07,-0.0504,
 672                              7.5547774216523998e-07,-0.0378,7.6662464387464e-07,-0.0252,7.5088461538461995e-07,-0.0126,
 673                              7.4313924501423998e-07,0.0,7.486323005698e-07,0.0126,7.4312054843304996e-07,0.0252,
 674                              7.5084722222222e-07,0.0378,7.6656730769230998e-07,0.0504,7.5540170940170999e-07,0.063,
 675                              7.5856267806267997e-07,0.075600000000000001,7.7204415954415998e-07,0.088200000000000001,
 676                              7.5410790598291005e-07,0.1008,7.3665776353275997e-07,),
 677                       PROL_DROIT='CONSTANT',);
 678
 679 FLY_1_9=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 680                       NOM_PARA='Y',
 681                       VALE=
 682                       (-0.1008,7.7370940170939995e-07,-0.088200000000000001,8.0144764957264998e-07,
 683                              -0.075600000000000001,8.0374358974358997e-07,-0.063,8.0603952991452995e-07,-0.0504,
 684                              8.0238372507123005e-07,-0.0378,7.9936672008546997e-07,-0.0252,7.9634971509971995e-07,-0.0126,
 685                              7.8317361111111004e-07,0.0,7.8316426282051005e-07,0.0126,7.8315491452990996e-07,0.0252,
 686                              7.9631107549857995e-07,0.0378,7.9930689102564003e-07,0.0504,8.0230270655270996e-07,0.063,
 687                              8.0593856837607005e-07,0.075600000000000001,8.0362580128205003e-07,0.088200000000000001,
 688                              8.0131303418803e-07,0.1008,7.7356356837607e-07,),
 689                       PROL_DROIT='CONSTANT',);
 690
 691 FLY_1_10=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 692                        NOM_PARA='Y',
 693                        VALE=
 694                        (-0.1008,7.9424198717948999e-07,-0.088200000000000001,8.1225925925926001e-07,
 695                               -0.075600000000000001,8.3075516381766e-07,-0.063,8.1540651709402001e-07,-0.0504,
 696                               8.1131071937321997e-07,-0.0378,8.2297489316238996e-07,-0.0252,8.0577403846154003e-07,-0.0126,
 697                               7.9717236467236002e-07,0.0,8.0284615384615001e-07,0.0126,7.9715117521368004e-07,0.0252,
 698                               8.0573165954416005e-07,0.0378,8.2291257122506995e-07,0.0504,8.1122845441594998e-07,0.063,
 699                               8.1530306267805998e-07,0.075600000000000001,8.3063051994301999e-07,0.088200000000000001,
 700                               8.1211841168090999e-07,0.1008,7.9408742877492998e-07,),
 701                        PROL_DROIT='CONSTANT',);
 702
 703 FLY_1_11=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 704                        NOM_PARA='Y',
 705                        VALE=
 706                        (-0.1008,8.1647845441594995e-07,-0.088200000000000001,8.3510648148147996e-07,
 707                               -0.075600000000000001,8.5487873931623998e-07,-0.063,8.3969088319087998e-07,-0.0504,
 708                               8.3502920227919995e-07,-0.0378,8.4651014957265e-07,-0.0252,8.2973931623932002e-07,-0.0126,
 709                               8.2549269943020005e-07,0.0,8.3637909544160002e-07,0.0126,8.2547275641026002e-07,0.0252,
 710                               8.2969319800569997e-07,0.0378,8.4644408831909002e-07,0.0504,8.3493945868945997e-07,0.063,
 711                               8.3957995014244996e-07,0.075600000000000001,8.5474786324785998e-07,0.088200000000000001,
 712                               8.3495690883191003e-07,0.1008,8.1631641737892005e-07,),
 713                        PROL_DROIT='CONSTANT',);
 714
 715 FLY_1_12=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 716                        NOM_PARA='Y',
 717                        VALE=
 718                        (-0.1008,8.3948272792023003e-07,-0.088200000000000001,8.7002546296295996e-07,
 719                               -0.075600000000000001,8.7608564814814997e-07,-0.063,8.8214583333333004e-07,-0.0504,
 720                               8.7589245014244996e-07,-0.0378,8.7128872863248004e-07,-0.0252,8.6668500712251001e-07,-0.0126,
 721                               8.6321616809117001e-07,0.0,8.6320432692307996e-07,0.0126,8.6319248575499001e-07,0.0252,
 722                               8.6663764245014005e-07,0.0378,8.7121705840456005e-07,0.0504,8.7579647435896998e-07,0.063,
 723                               8.8202617521368001e-07,0.075600000000000001,8.7594667022791999e-07,0.088200000000000001,
 724                               8.6986716524217003e-07,0.1008,8.3931196581196996e-07,),
 725                        PROL_DROIT='CONSTANT',);
 726
 727 FLY_1_13=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 728                        NOM_PARA='Y',
 729                        VALE=
 730                        (-0.1008,8.5508315527066002e-07,-0.088200000000000001,8.7546616809116998e-07,
 731                               -0.075600000000000001,9.0944533475783005e-07,-0.063,9.1292663817663998e-07,-0.0504,
 732                               8.9653721509972004e-07,-0.0378,8.9547275641025999e-07,-0.0252,8.7405644586895003e-07,-0.0126,
 733                               8.7004789886040001e-07,0.0,8.8226549145299002e-07,0.0126,8.7002546296295996e-07,0.0252,
 734                               8.7400534188033999e-07,0.0378,8.9539672364672004e-07,0.0504,8.9643376068376001e-07,0.063,
 735                               9.1279825498575999e-07,0.075600000000000001,9.0929576210826002e-07,0.088200000000000001,
 736                               8.7530288461538001e-07,0.1008,8.5490491452990995e-07,),
 737                        PROL_DROIT='CONSTANT',);
 738
 739 FLY_1_14=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 740                        NOM_PARA='Y',
 741                        VALE=
 742                        (-0.1008,8.5798237179486999e-07,-0.088200000000000001,8.6880893874644002e-07,
 743                               -0.075600000000000001,9.0299750712251001e-07,-0.063,9.0948646723646997e-07,-0.0504,
 744                               9.1597542735043003e-07,-0.0378,9.0500302706552995e-07,-0.0252,8.8214084757835001e-07,-0.0126,
 745                               8.7779825498574995e-07,0.0,8.8993732193732e-07,0.0126,8.7777207977207998e-07,0.0252,
 746                               8.8208974358973997e-07,0.0378,9.0492325498575998e-07,0.0504,9.1586574074073996e-07,0.063,
 747                               9.0935434472933995e-07,0.075600000000000001,9.0284294871795e-07,0.088200000000000001,
 748                               8.6864066951567002e-07,0.1008,8.5779665242165003e-07,),
 749                        PROL_DROIT='CONSTANT',);
 750
 751 FLY_1_15=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 752                        NOM_PARA='Y',
 753                        VALE=
 754                        (-0.1008,8.5773557692308004e-07,-0.088200000000000001,8.6180893874644003e-07,
 755                               -0.075600000000000001,8.7862838319088002e-07,-0.063,8.9924074074074004e-07,-0.0504,
 756                               9.0872489316239005e-07,-0.0378,9.0174545940170996e-07,-0.0252,8.9476602564102997e-07,-0.0126,
 757                               8.9122364672365005e-07,0.0,8.9120993589743996e-07,0.0126,8.9119622507122997e-07,0.0252,
 758                               8.9471118233617996e-07,0.0378,9.0166132478632003e-07,0.0504,9.0861146723646996e-07,0.063,
 759                               8.9910363247862999e-07,0.075600000000000001,8.7847257834757996e-07,0.088200000000000001,
 760                               8.6163693019942995e-07,0.1008,8.5754362535613004e-07,),
 761                        PROL_DROIT='CONSTANT',);
 762
 763 FLY_1_16=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 764                        NOM_PARA='Y',
 765                        VALE=
 766                        (-0.1008,8.5903935185184999e-07,-0.088200000000000001,8.5690046296295998e-07,
 767                               -0.075600000000000001,8.5831018518518999e-07,-0.063,8.6167307692307999e-07,-0.0504,
 768                               8.7125320512820998e-07,-0.0378,8.8563586182335997e-07,-0.0252,8.7089298433048003e-07,-0.0126,
 769                               8.6842877492877004e-07,0.0,8.7942111823362e-07,0.0126,8.6840135327634996e-07,0.0252,
 770                               8.7083814102563997e-07,0.0378,8.8554985754985996e-07,0.0504,8.7113977920228004e-07,0.063,
 771                               8.6153846153846001e-07,0.075600000000000001,8.5815313390312996e-07,0.088200000000000001,
 772                               8.5672471509971999e-07,0.1008,8.5884116809116995e-07,),
 773                        PROL_DROIT='CONSTANT',);
 774
 775 FLY_1_17=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 776                        NOM_PARA='Y',
 777                        VALE=
 778                        (-0.1008,8.6388675213674999e-07,-0.088200000000000001,8.5563782051281999e-07,
 779                               -0.075600000000000001,8.5087019230768996e-07,-0.063,8.4755965099714995e-07,-0.0504,
 780                               8.4761698717949003e-07,-0.0378,8.5133636039885996e-07,-0.0252,8.4844088319088001e-07,-0.0126,
 781                               8.4633689458689002e-07,0.0,8.4644658119657999e-07,0.0126,8.4630947293447004e-07,0.0252,
 782                               8.4838479344729004e-07,0.0378,8.512516025641e-07,0.0504,8.4750480769231e-07,0.063,
 783                               8.4742378917379002e-07,0.075600000000000001,8.5071064814815003e-07,0.088200000000000001,
 784                               8.5545584045584e-07,0.1008,8.6368108974359e-07,),
 785                        PROL_DROIT='CONSTANT',);
 786
 787 FLUXRD1=DEFI_NAPPE(FONCTION=(FLY_1_1,FLY_1_2,FLY_1_3,FLY_1_4,FLY_1_5,FLY_1_6,FLY_1_7,FLY_1_8,FLY_1_9,FLY_1_10,FLY_1_11,FLY_1_12,
 788                               FLY_1_13,FLY_1_14,FLY_1_15,FLY_1_16,FLY_1_17,),
 789                    PARA=
 790                    (-0.1008,-0.088200000000000001,-0.075600000000000001,-0.063,-0.0504,-0.0378,-0.0252,
 791                           -0.0126,0.0,0.0126,0.0252,0.0378,0.0504,0.063,0.075600000000000001,0.088200000000000001,0.1008,),
 792                    PROL_GAUCHE='CONSTANT',
 793                    NOM_PARA='Z',
 794                    PROL_DROIT='CONSTANT',);
 795 #
 796 # DEFINITION DES MATERIAUX ET AFFECTATIONS
 797 #
 798 #           **** DEFINITION DES MATERIAUX SANS FLUX NEUTRONIQUE ****
 799 #
 800 #
 801 # FAMILLE ASSEMBLAGE : NUMERO_FAM
 802 #
 803 # EMBOUT SUPERIEUR
 804
 805
 806 E_ES=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 807                    NOM_PARA='TEMP',
 808                    VALE=(100.0,2.0,300.0,3.0,),
 809                    PROL_DROIT='CONSTANT',);
 810
 811 NU_ES=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 812                     NOM_PARA='TEMP',
 813                     VALE=(100.0,2.0,300.0,3.0,),
 814                     PROL_DROIT='CONSTANT',);
 815
 816 AL_ES=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 817                     NOM_PARA='TEMP',
 818                     VALE=(100.0,2.0,300.0,3.0,),
 819                     PROL_DROIT='CONSTANT',);
 820
 821 MT_ES=DEFI_MATERIAU(ELAS_FO=_F(NU=NU_ES,
 822                                ALPHA=AL_ES,
 823                                RHO=7800.0,
 824                                E=E_ES,),);
 825 # EMBOUT INFERIEUR
 826
 827
 828 E_EI=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 829                    NOM_PARA='TEMP',
 830                    VALE=(100.0,2.0,300.0,3.0,),
 831                    PROL_DROIT='CONSTANT',);
 832
 833 NU_EI=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 834                     NOM_PARA='TEMP',
 835                     VALE=(100.0,2.0,300.0,3.0,),
 836                     PROL_DROIT='CONSTANT',);
 837
 838 AL_EI=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 839                     NOM_PARA='TEMP',
 840                     VALE=(100.0,2.0,300.0,3.0,),
 841                     PROL_DROIT='CONSTANT',);
 842
 843 MT_EI=DEFI_MATERIAU(ELAS_FO=_F(NU=NU_EI,
 844                                ALPHA=AL_EI,
 845                                RHO=12.0,
 846                                E=E_EI,),);
 847
 848 MT_RIG=DEFI_MATERIAU(ELAS=_F(NU=0.0,
 849                              ALPHA=0.0,
 850                              RHO=0.0,
 851                              E=1.0,),);
 852 # FONCTION COUPLE EN FONCTION DE DRZ POUR LA LIAISON GRILLE/CRAYON
 853
 854
 855 FOTRAC=DEFI_FONCTION(PROL_GAUCHE='LINEAIRE',
 856                      NOM_PARA='DRZ',
 857                      VALE=(1.0,1.0,0.0,0.0,0.0,0.0,0.0,0.0,),
 858                      PROL_DROIT='LINEAIRE',);
 859
 860 MT_RE1=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 861                                     RELA_MZ=FOTRAC,
 862                                     RIGI_N_FO=FRIG1,
 863                                     COULOMB=1.0,),
 864                      ELAS=_F(NU=0.0,
 865                              ALPHA=0.0,
 866                              RHO=0.0,
 867                              E=1.0,),);
 868
 869 MT_RE2=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 870                                     RELA_MZ=FOTRAC,
 871                                     RIGI_N_FO=FRIG2,
 872                                     COULOMB=1.0,),
 873                      ELAS=_F(NU=0.0,
 874                              ALPHA=0.0,
 875                              RHO=0.0,
 876                              E=1.0,),);
 877
 878 MT_RE3=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 879                                     RELA_MZ=FOTRAC,
 880                                     RIGI_N_FO=FRIG3,
 881                                     COULOMB=1.0,),
 882                      ELAS=_F(NU=0.0,
 883                              ALPHA=0.0,
 884                              RHO=0.0,
 885                              E=1.0,),);
 886
 887 MT_RE4=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 888                                     RELA_MZ=FOTRAC,
 889                                     RIGI_N_FO=FRIG4,
 890                                     COULOMB=1.0,),
 891                      ELAS=_F(NU=0.0,
 892                              ALPHA=0.0,
 893                              RHO=0.0,
 894                              E=1.0,),);
 895
 896 MT_RE5=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 897                                     RELA_MZ=FOTRAC,
 898                                     RIGI_N_FO=FRIG5,
 899                                     COULOMB=1.0,),
 900                      ELAS=_F(NU=0.0,
 901                              ALPHA=0.0,
 902                              RHO=0.0,
 903                              E=1.0,),);
 904
 905 MT_RE6=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 906                                     RELA_MZ=FOTRAC,
 907                                     RIGI_N_FO=FRIG6,
 908                                     COULOMB=1.0,),
 909                      ELAS=_F(NU=0.0,
 910                              ALPHA=0.0,
 911                              RHO=0.0,
 912                              E=1.0,),);
 913
 914 MT_RE7=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 915                                     RELA_MZ=FOTRAC,
 916                                     RIGI_N_FO=FRIG7,
 917                                     COULOMB=1.0,),
 918                      ELAS=_F(NU=0.0,
 919                              ALPHA=0.0,
 920                              RHO=0.0,
 921                              E=1.0,),);
 922
 923 MT_RE8=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 924                                     RELA_MZ=FOTRAC,
 925                                     RIGI_N_FO=FRIG8,
 926                                     COULOMB=1.0,),
 927                      ELAS=_F(NU=0.0,
 928                              ALPHA=0.0,
 929                              RHO=0.0,
 930                              E=1.0,),);
 931
 932 MT_RE9=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 933                                     RELA_MZ=FOTRAC,
 934                                     RIGI_N_FO=FRIG9,
 935                                     COULOMB=1.0,),
 936                      ELAS=_F(NU=0.0,
 937                              ALPHA=0.0,
 938                              RHO=0.0,
 939                              E=1.0,),);
 940
 941 MT_RE10=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 942                                      RELA_MZ=FOTRAC,
 943                                      RIGI_N_FO=FRIG10,
 944                                      COULOMB=1.0,),
 945                       ELAS=_F(NU=0.0,
 946                               ALPHA=0.0,
 947                               RHO=0.0,
 948                               E=1.0,),);
 949 # SYSTEME DE MAINTIEN
 950
 951
 952 E_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 953                    NOM_PARA='TEMP',
 954                    VALE=(0.0,1.0,5.0,2.0,),
 955                    PROL_DROIT='CONSTANT',);
 956
 957 NU_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 958                     NOM_PARA='TEMP',
 959                     VALE=(0.0,1.0,5.0,2.0,),
 960                     PROL_DROIT='CONSTANT',);
 961
 962 AL_MA=DEFI_CONSTANTE(VALE=0.29999999999999999,);
 963
 964 SY_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 965                     NOM_PARA='TEMP',
 966                     VALE=(0.0,1.0,5.0,2.0,),
 967                     PROL_DROIT='CONSTANT',);
 968
 969 DS_MA=DEFI_CONSTANTE(VALE=0.29999999999999999,);
 970
 971 NP_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 972                     NOM_PARA='TEMP',
 973                     VALE=(0.0,1.0,5.0,2.0,),
 974                     PROL_DROIT='CONSTANT',);
 975
 976 MEY_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 977                      NOM_PARA='TEMP',
 978                      VALE=(0.0,1.0,5.0,2.0,),
 979                      PROL_DROIT='CONSTANT',);
 980
 981 MPY_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 982                      NOM_PARA='TEMP',
 983                      VALE=(0.0,1.0,5.0,2.0,),
 984                      PROL_DROIT='CONSTANT',);
 985
 986 MEZ_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 987                      NOM_PARA='TEMP',
 988                      VALE=(0.0,1.0,5.0,2.0,),
 989                      PROL_DROIT='CONSTANT',);
 990
 991 MPZ_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
 992                      NOM_PARA='TEMP',
 993                      VALE=(0.0,1.0,5.0,2.0,),
 994                      PROL_DROIT='CONSTANT',);
 995
 996 CAY_MA=DEFI_CONSTANTE(VALE=0.29999999999999999,);
 997
 998 CBY_MA=DEFI_CONSTANTE(VALE=0.29999999999999999,);
 999
1000 CAZ_MA=DEFI_CONSTANTE(VALE=0.29999999999999999,);
1001
1002 CBZ_MA=DEFI_CONSTANTE(VALE=0.29999999999999999,);
1003
1004 MPX_MA=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1005                      NOM_PARA='TEMP',
1006                      VALE=(0.0,1.0,5.0,2.0,),
1007                      PROL_DROIT='CONSTANT',);
1008
1009 MT_MA=DEFI_MATERIAU(VMIS_POUTRE_FO=_F(MEZ=MEZ_MA,
1010                                       MPX=MPX_MA,
1011                                       MPZ=MPZ_MA,
1012                                       CBY=CBY_MA,
1013                                       NP=NP_MA,
1014                                       CBZ=CBZ_MA,
1015                                       MEY=MEY_MA,
1016                                       CAY=CAY_MA,
1017                                       CAZ=CAZ_MA,
1018                                       MPY=MPY_MA,),
1019                     ECRO_LINE_FO=_F(D_SIGM_EPSI=DS_MA,
1020                                     SY=SY_MA,),
1021                     ELAS_FO=_F(NU=NU_MA,
1022                                ALPHA=AL_MA,
1023                                RHO=7800.0,
1024                                E=E_MA,),);
1025 #
1026 #           **** DEFINITION DES MATERIAUX AVEC FLUX NEUTRONIQUE ****
1027 #
1028 #
1029 # GRILLE
1030
1031
1032 E_GRE=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1033                     NOM_PARA='TEMP',
1034                     VALE=(0.0,1.0,5.0,2.0,),
1035                     PROL_DROIT='CONSTANT',);
1036
1037 E_GRM=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1038                     NOM_PARA='TEMP',
1039                     VALE=(0.0,1.0,5.0,2.0,),
1040                     PROL_DROIT='CONSTANT',);
1041
1042 NU_GRE=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1043                      NOM_PARA='TEMP',
1044                      VALE=(0.0,1.0,5.0,2.0,),
1045                      PROL_DROIT='CONSTANT',);
1046
1047 NU_GRM=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1048                      NOM_PARA='TEMP',
1049                      VALE=(0.0,1.0,5.0,2.0,),
1050                      PROL_DROIT='CONSTANT',);
1051
1052 AL_GRE=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1053                      NOM_PARA='TEMP',
1054                      VALE=(0.0,1.0,5.0,2.0,),
1055                      PROL_DROIT='CONSTANT',);
1056
1057 AL_GRM=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1058                      NOM_PARA='TEMP',
1059                      VALE=(0.0,1.0,5.0,2.0,),
1060                      PROL_DROIT='CONSTANT',);
1061 # TUBE GUIDE
1062
1063
1064 E_TU=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1065                    NOM_PARA='TEMP',
1066                    VALE=(0.0,1.0,5.0,2.0,),
1067                    PROL_DROIT='CONSTANT',);
1068
1069 NU_TU=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1070                     NOM_PARA='TEMP',
1071                     VALE=(0.0,1.0,5.0,2.0,),
1072                     PROL_DROIT='CONSTANT',);
1073
1074 AL_TU=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1075                     NOM_PARA='TEMP',
1076                     VALE=(0.0,1.0,5.0,2.0,),
1077                     PROL_DROIT='CONSTANT',);
1078 # CRAYON
1079
1080
1081 E_CR=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1082                    NOM_PARA='TEMP',
1083                    VALE=(0.0,1.0,5.0,2.0,),
1084                    PROL_DROIT='CONSTANT',);
1085
1086 NU_CR=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1087                     NOM_PARA='TEMP',
1088                     VALE=(0.0,1.0,5.0,2.0,),
1089                     PROL_DROIT='CONSTANT',);
1090
1091 AL_CR=DEFI_FONCTION(PROL_GAUCHE='CONSTANT',
1092                     NOM_PARA='TEMP',
1093                     VALE=(0.0,1.0,5.0,2.0,),
1094                     PROL_DROIT='CONSTANT',);
1095 # GRANDISSEMENT DES GRILLES
1096 # GRILLES EXTREMES 1 ET 10
1097
1098
1099 MT_GRE1=DEFI_MATERIAU(GRAN_IRRA=_F(FONC_MULT=FLUX_F1,
1100                                    FLUX_L=FLUXAX1,
1101                                    FLUX_TN=FLUXRD1,
1102                                    S=1.0,
1103                                    A=0.0,),
1104                       LEMAITRE=_F(N=0.0,
1105                                   UN_SUR_M=0.0,
1106                                   QSR_K=1.0,
1107                                   UN_SUR_K=0.0,),
1108                       ELAS_FO=_F(NU=NU_GRE,
1109                                  ALPHA=AL_GRE,
1110                                  RHO=7800.0,
1111                                  E=E_GRE,),);
1112 # GRILLES DE MELANGE
1113
1114
1115 MT_GRM1=DEFI_MATERIAU(GRAN_IRRA=_F(FONC_MULT=FLUX_F1,
1116                                    FLUX_L=FLUXAX1,
1117                                    FLUX_TN=FLUXRD1,
1118                                    S=1.0,
1119                                    A=0.0,),
1120                       LEMAITRE=_F(N=0.0,
1121                                   UN_SUR_M=0.0,
1122                                   QSR_K=2.0,
1123                                   UN_SUR_K=0.0,),
1124                       ELAS_FO=_F(NU=NU_GRM,
1125                                  ALPHA=AL_GRM,
1126                                  RHO=7800.0,
1127                                  E=E_GRM,),);
1128 # DEFINITION DES MATERIAUX VISCOPLASTIQUES ( TUBES ET CRAYONS )
1129
1130
1131 MT_TU1=DEFI_MATERIAU(GRAN_IRRA=_F(FONC_MULT=FLUX_F1,
1132                                   FLUX_L=FLUXAX1,
1133                                   FLUX_TN=FLUXRD1,
1134                                   S=1.0,
1135                                   A=0.0,),
1136                      LEMAITRE=_F(N=0.0,
1137                                  UN_SUR_M=0.0,
1138                                  QSR_K=2.0,
1139                                  UN_SUR_K=0.0,),
1140                      ELAS_FO=_F(NU=NU_TU,
1141                                 ALPHA=AL_TU,
1142                                 RHO=7800.0,
1143                                 E=E_TU,),);
1144
1145 MT_CR1=DEFI_MATERIAU(GRAN_IRRA=_F(FONC_MULT=FLUX_F1,
1146                                   FLUX_L=FLUXAX1,
1147                                   FLUX_TN=FLUXRD1,
1148                                   S=1.0,
1149                                   A=0.0,),
1150                      LEMAITRE=_F(N=0.0,
1151                                  UN_SUR_M=0.0,
1152                                  QSR_K=2.0,
1153                                  UN_SUR_K=0.0,),
1154                      ELAS_FO=_F(NU=NU_CR,
1155                                 ALPHA=AL_CR,
1156                                 RHO=7800.0,
1157                                 E=E_CR,),);
1158 #
1159 #           AFFECTATION DES MATERIAUX A CHAUD - CYCLE
1160 #
1161 #
1162
1163
1164 A_MAT_1=AFFE_MATERIAU(AFFE=(_F(TEMP_REF=100.0,
1165                                MATER=MT_MA,
1166                                GROUP_MA=('MAINTI_Y','MAINTI_Z',),),
1167                             _F(TEMP_REF=100.0,
1168                                MATER=MT_GRM1,
1169                                GROUP_MA='GRIL_I',),
1170                             _F(TEMP_REF=100.0,
1171                                MATER=MT_GRE1,
1172                                GROUP_MA='GRIL_E',),
1173                             _F(TEMP_REF=100.0,
1174                                MATER=MT_CR1,
1175                                GROUP_MA='CRAY',),
1176                             _F(TEMP_REF=100.0,
1177                                MATER=MT_TU1,
1178                                GROUP_MA='GUID',),
1179                             _F(TEMP_REF=100.0,
1180                                MATER=MT_EI,
1181                                GROUP_MA='EBOINF',),
1182                             _F(TEMP_REF=100.0,
1183                                MATER=MT_EI,
1184                                GROUP_MA=('PI_PLE','PI_TRO',),),
1185                             _F(TEMP_REF=100.0,
1186                                MATER=MT_ES,
1187                                GROUP_MA='EBOSUP',),
1188                             _F(TEMP_REF=100.0,
1189                                MATER=MT_RIG,
1190                                GROUP_MA='RIGID',),
1191                             _F(TEMP_REF=100.0,
1192                                MATER=MT_RE1,
1193                                GROUP_MA='ELA1',),
1194                             _F(TEMP_REF=100.0,
1195                                MATER=MT_RE2,
1196                                GROUP_MA='ELA2',),
1197                             _F(TEMP_REF=100.0,
1198                                MATER=MT_RE3,
1199                                GROUP_MA='ELA3',),
1200                             _F(TEMP_REF=100.0,
1201                                MATER=MT_RE4,
1202                                GROUP_MA='ELA4',),
1203                             _F(TEMP_REF=100.0,
1204                                MATER=MT_RE5,
1205                                GROUP_MA='ELA5',),
1206                             _F(TEMP_REF=100.0,
1207                                MATER=MT_RE6,
1208                                GROUP_MA='ELA6',),
1209                             _F(TEMP_REF=100.0,
1210                                MATER=MT_RE7,
1211                                GROUP_MA='ELA7',),
1212                             _F(TEMP_REF=100.0,
1213                                MATER=MT_RE8,
1214                                GROUP_MA='ELA8',),
1215                             _F(TEMP_REF=100.0,
1216                                MATER=MT_RE9,
1217                                GROUP_MA='ELA9',),
1218                             _F(TEMP_REF=100.0,
1219                                MATER=MT_RE10,
1220                                GROUP_MA='ELA10',),),
1221                       MAILLAGE=MA1,);
1222 #
1223
1224
1225 L_INST_1=DEFI_LIST_REEL(INTERVALLE=(_F(JUSQU_A=0.5,
1226                                        NOMBRE=1,),
1227                                     _F(JUSQU_A=1000.0,
1228                                        NOMBRE=1,),
1229                                     _F(JUSQU_A=2000.0,
1230                                        NOMBRE=1,),
1231                                     _F(JUSQU_A=3000.0,
1232                                        NOMBRE=5,),
1233                                     _F(JUSQU_A=4000.0,
1234                                        NOMBRE=1,),
1235                                     _F(JUSQU_A=5000.0,
1236                                        NOMBRE=1,),
1237                                     _F(JUSQU_A=6000.0,
1238                                        NOMBRE=1,),
1239                                     _F(JUSQU_A=7000.0,
1240                                        NOMBRE=1,),
1241                                     _F(JUSQU_A=8000.0,
1242                                        NOMBRE=1,),),
1243                         DEBUT=0.0,);
1244
1245 L_ARCH_1=DEFI_LIST_REEL(INTERVALLE=(_F(JUSQU_A=0.5,
1246                                        NOMBRE=1,),
1247                                     _F(JUSQU_A=1000.0,
1248                                        NOMBRE=1,),
1249                                     _F(JUSQU_A=2000.0,
1250                                        NOMBRE=1,),
1251                                     _F(JUSQU_A=3000.0,
1252                                        NOMBRE=5,),
1253                                     _F(JUSQU_A=4000.0,
1254                                        NOMBRE=1,),
1255                                     _F(JUSQU_A=5000.0,
1256                                        NOMBRE=1,),
1257                                     _F(JUSQU_A=6000.0,
1258                                        NOMBRE=1,),
1259                                     _F(JUSQU_A=7000.0,
1260                                        NOMBRE=1,),
1261                                     _F(JUSQU_A=8000.0,
1262                                        NOMBRE=1,),),
1263                         DEBUT=0.0,);
1264 #
1265 # CALCUL PHASE I A V DU CYCLE 1
1266 #
1267
1268
1269 RESU=STAT_NON_LINE(CHAM_MATER=A_MAT_1,
1270                    MODELE=MO1,
1271                    ARCHIVAGE=_F(LIST_INST=L_ARCH_1,),
1272                    CONVERGENCE=_F(ITER_GLOB_MAXI=50,
1273                                   ITER_INTE_MAXI=500,
1274                                   RESI_GLOB_RELA=1.E-4,
1275                                   ARRET='OUI',),
1276                    COMP_INCR=(_F(RELATION='VMIS_POU_LINE',
1277                                  GROUP_MA='MAINTIEN',),
1278                               _F(RELATION='ASSE_COMBU',
1279                                  GROUP_MA=('CRAYONS','TUB_GUID',),),
1280                               _F(RELATION='DIS_CONTACT',
1281                                  GROUP_MA='LI_GR_CR',),
1282                               _F(COQUE_NCOU=1,
1283                                  RELATION='ELAS',
1284                                  GROUP_MA=('EMBO_INF','EMBO_SUP','PIEDS','GRILLES','LI_GR_GU',),),),
1285                    INCREMENT=_F(NUME_INST_FIN=20,
1286                                 LIST_INST=L_INST_1,),
1287                    CARA_ELEM=CARA1,
1288                    RECH_LINEAIRE=_F(ITER_LINE_MAXI=3,),
1289                    NEWTON=_F(REAC_ITER=1,
1290                              MATRICE='TANGENTE',),
1291                    SOLVEUR=_F(METHODE='MULT_FRONT',
1292                               RENUM='MD',),
1293                    EXCIT=(_F(FONC_MULT=ARCH_F1,
1294                              CHARGE=ARCH_1,),
1295                           _F(FONC_MULT=ARCH_F1,
1296                              CHARGE=FOARCH_1,),
1297                           _F(FONC_MULT=HYDR_F1,
1298                              CHARGE=HYDR_1,),
1299                           _F(FONC_MULT=HYDR_F1,
1300                              CHARGE=FOHYDR_1,),
1301                           _F(FONC_MULT=CLDI_F1,
1302                              CHARGE=CL_DI_1,),
1303                           _F(CHARGE=CL_PER_1,),
1304                           _F(CHARGE=CL_CTC,),
1305                           _F(CHARGE=PESANT_1,),
1306                           _F(CHARGE=THERM_1,),),
1307                    PARM_THETA=0.5,);
1308
1309 FIN();