Exemples/ex22/tyty.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 PRE_GIBI();
  11 #
  12
  13 MA1=LIRE_MAILLAGE();
  14 #
  15
  16 MA1=DEFI_GROUP(reuse =MA1,
  17                MAILLAGE=MA1,
  18                CREA_GROUP_NO=_F(GROUP_MA=('FIXATI','APPUI','LIAS_1','LIAS_2','EBGUIDE','MAINTIEN','MAIN_SUP','EMBO_SUP','TUB_GUID',
  19                                            'GRILLES','EMBO_INF','PIEDS','CRAYONS',),),);
  20 #
  21 # AFFECTATION DU MODELE
  22
  23 MO1=AFFE_MODELE(MAILLAGE=MA1,
  24                 AFFE=(_F(GROUP_MA='CRAYONS',
  25                          PHENOMENE='MECANIQUE',
  26                          MODELISATION='POU_D_E',),
  27                       _F(GROUP_MA='TUB_GUID',
  28                          PHENOMENE='MECANIQUE',
  29                          MODELISATION='POU_D_E',),
  30                       _F(GROUP_MA=('EMBO_SUP','EMBO_INF',),
  31                          PHENOMENE='MECANIQUE',
  32                          MODELISATION='DKT',),
  33                       _F(GROUP_MA='PIEDS',
  34                          PHENOMENE='MECANIQUE',
  35                          MODELISATION='POU_D_E',),
  36                       _F(GROUP_MA='GRILLES',
  37                          PHENOMENE='MECANIQUE',
  38                          MODELISATION='POU_D_E',),
  39                       _F(GROUP_MA='LI_GR_CR',
  40                          PHENOMENE='MECANIQUE',
  41                          MODELISATION='DIS_TR',),
  42                       _F(GROUP_MA='LI_GR_GU',
  43                          PHENOMENE='MECANIQUE',
  44                          MODELISATION='DIS_TR',),
  45                       _F(GROUP_MA='MAINTIEN',
  46                          PHENOMENE='MECANIQUE',
  47                          MODELISATION='POU_D_E',),),);
  48 #
  49 # AFFECTATION DES CARACTERISTIQUES GEOMETRIQUES ELEMENTAIRES
  50
  51 CARA1=AFFE_CARA_ELEM(MODELE=MO1,
  52                      POUTRE=(_F(GROUP_MA='CRAY',
  53                                 SECTION='CERCLE',
  54                                 CARA=('R','EP',),
  55                                 VALE=(3.0,3.0,),),
  56                              _F(GROUP_MA='GUID',
  57                                 SECTION='CERCLE',
  58                                 CARA=('R','EP',),
  59                                 VALE=(3.0,3.0,),),
  60                              _F(GROUP_MA='RETRE',
  61                                 SECTION='CERCLE',
  62                                 CARA=('R','EP',),
  63                                 VALE=(3.0,3.0,),),
  64                              _F(GROUP_MA='GRIL_I',
  65                                 SECTION='RECTANGLE',
  66                                 CARA=('HY','HZ',),
  67                                 VALE=(3.0,3.0,),),
  68                              _F(GROUP_MA='GRIL_E',
  69                                 SECTION='RECTANGLE',
  70                                 CARA=('HY','HZ',),
  71                                 VALE=(3.0,3.0,),),
  72                              _F(GROUP_MA=('MAINTI_Y','MAINTI_Z',),
  73                                 SECTION='RECTANGLE',
  74                                 CARA=('HY','HZ',),
  75                                 VALE=(3.0,3.0,),),
  76                              _F(GROUP_MA='PI_PLE',
  77                                 SECTION='GENERALE',
  78                                 CARA=('A','IZ','IY','JX',),
  79                                 VALE=(1.0,2.0,3.0,4.0,),),
  80                              _F(GROUP_MA='PI_TRO',
  81                                 SECTION='GENERALE',
  82                                 CARA=('A','IZ','IY','JX',),
  83                                 VALE=(1.0,2.0,3.0,4.0,),),),
  84                      ORIENTATION=(_F(GROUP_MA=('GRIL_I','GRIL_E',),
  85                                      VALE=(1.0,0.0,0.0,),
  86                                      CARA='VECT_Y',),
  87                                   _F(GROUP_MA='MAINTI_Y',
  88                                      VALE=(0.0,1.0,0.0,),
  89                                      CARA='VECT_Y',),
  90                                   _F(GROUP_MA='MAINTI_Z',
  91                                      VALE=(0.0,0.0,1.0,),
  92                                      CARA='VECT_Y',),
  93                                   _F(GROUP_MA=('PI_PLE','PI_TRO',),
  94                                      VALE=(0.0,1.0,0.0,),
  95                                      CARA='VECT_Y',),
  96                                   _F(GROUP_MA=('ELASTI','RIGID',),
  97                                      VALE=(1.0,0.0,0.0,),
  98                                      CARA='VECT_Y',),),
  99                      DISCRET=(_F(GROUP_MA='ELASTI',
 100                                  CARA='K_TR_D_L',
 101                                  VALE=(1.0,2.0,3.0,4.0,5.0,6.0,),
 102                                  REPERE='LOCAL',),
 103                               _F(GROUP_MA='RIGID',
 104                                  CARA='K_TR_D_L',
 105                                  VALE=(1.0,2.0,3.0,4.0,5.0,6.0,),
 106                                  REPERE='LOCAL',),),
 107                      COQUE=(_F(EPAIS=2.0,
 108                                GROUP_MA='EBOINF',
 109                                ANGL_REP=90.0,),
 110                             _F(EPAIS=3.0,
 111                                GROUP_MA='EBOSUP',
 112                                ANGL_REP=90.0,),),);
 113 #
 114 # DEFINITION DES CHARGES ET CONDITIONS LIMITES
 115 # DEFINITION DES TEMPERATURES NODALES EVOLUTIVES
 116
 117 F_TP1_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
 118                             5.0,4.0,
 119                             ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',NOM_RESU='TEMP',);
 120
 121 F_TP2_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
 122                             5.0,4.0,
 123                             ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',NOM_RESU='TEMP',);
 124
 125 F_TP3_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
 126                             5.0,4.0,
 127                             ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',NOM_RESU='TEMP',);
 128
 129 CHTEM1_1=AFFE_CHAM_NO(MAILLAGE=MA1,
 130                       GRANDEUR='TEMP_F',
 131                       AFFE=(_F(NOM_CMP='TEMP',
 132                                GROUP_NO=('MAINTIEN','PIEDS','TUB_GUID','GRILLES','CRAYONS',),
 133                                FONCTION=F_TP1_1,),
 134                             _F(NOM_CMP=('TEMP_INF','TEMP','TEMP_SUP',),
 135                                GROUP_NO=('EMBO_SUP','EMBO_INF',),
 136                                FONCTION=(F_TP1_1,F_TP1_1,F_TP1_1,),),),);
 137
 138 CHTEM2_1=AFFE_CHAM_NO(MAILLAGE=MA1,
 139                       GRANDEUR='TEMP_F',
 140                       AFFE=(_F(NOM_CMP='TEMP',
 141                                GROUP_NO=('MAINTIEN','PIEDS','TUB_GUID','GRILLES','CRAYONS',),
 142                                FONCTION=F_TP2_1,),
 143                             _F(NOM_CMP=('TEMP_INF','TEMP','TEMP_SUP',),
 144                                GROUP_NO=('EMBO_SUP','EMBO_INF',),
 145                                FONCTION=(F_TP2_1,F_TP2_1,F_TP2_1,),),),);
 146
 147 CHTEM3_1=AFFE_CHAM_NO(MAILLAGE=MA1,
 148                       GRANDEUR='TEMP_F',
 149                       AFFE=(_F(NOM_CMP='TEMP',
 150                                GROUP_NO=('MAINTIEN','PIEDS','TUB_GUID','GRILLES',),
 151                                FONCTION=F_TP3_1,),
 152                             _F(NOM_CMP=('TEMP_INF','TEMP','TEMP_SUP',),
 153                                GROUP_NO=('EMBO_SUP','EMBO_INF',),
 154                                FONCTION=(F_TP3_1,F_TP3_1,F_TP3_1,),),
 155                             _F(NOM_CMP='TEMP',
 156                                GROUP_NO='CRAYONS',
 157                                FONCTION=F_TP3_1,),),);
 158
 159 CHTH_1=CREA_RESU(TYPE_RESU='EVOL_THER',
 160                  NOM_CHAM='TEMP',
 161                  CHAM_GD=(_F(CHAM_NO=CHTEM1_1,
 162                              INST=0.0,),
 163                           _F(CHAM_NO=CHTEM1_1,
 164                              INST=1000.0,),
 165                           _F(CHAM_NO=CHTEM2_1,
 166                              INST=2000.0,),
 167                           _F(CHAM_NO=CHTEM3_1,
 168                              INST=3000.0,),
 169                           _F(CHAM_NO=CHTEM3_1,
 170                              INST=4000.0,),
 171                           _F(CHAM_NO=CHTEM2_1,
 172                              INST=5000.0,),
 173                           _F(CHAM_NO=CHTEM1_1,
 174                              INST=6000.0,),
 175                           _F(CHAM_NO=CHTEM1_1,
 176                              INST=7000.0,),),);
 177 # CHARGEMENT: CHARGES THERMO-MECANIQUES
 178 # ----------------------------------------
 179
 180 THERM_1=AFFE_CHAR_MECA(TEMP_CALCULEE=CHTH_1,
 181                        MODELE=MO1,);
 182 # DEFINITION DES CONDITIONS LIMITES
 183 #            DES CONDITIONS DE CONTACT
 184 #
 185 # CHARGEMENT: PESANTEUR
 186
 187 PESANT_1=AFFE_CHAR_MECA(PESANTEUR=(9.81,-1.0,0.0,0.0,),
 188                         MODELE=MO1,);
 189 #
 190 # CONDITIONS AUX LIMITES
 191
 192 CL_PER_1=AFFE_CHAR_MECA(MODELE=MO1,
 193                         DDL_IMPO=(_F(DRY=0.0,
 194                                      DRX=0.0,
 195                                      DRZ=0.0,
 196                                      GROUP_NO='FIXATI',
 197                                      DY=0.0,
 198                                      DZ=0.0,
 199                                      DX=0.0,),
 200                                   _F(GROUP_NO='APPUI',
 201                                      DX=0.0,),
 202                                   _F(DRY=0.0,
 203                                      DRX=0.0,
 204                                      DRZ=0.0,
 205                                      GROUP_NO='EBGUIDE',
 206                                      DY=0.0,
 207                                      DZ=0.0,),
 208                                   _F(DRX=0.0,
 209                                      GROUP_NO=('LIAS_1','LIAS_2',),),),
 210                         LIAISON_GROUP=(_F(GROUP_NO_1='LIAS_1',
 211                                           GROUP_NO_2='LIAS_2',
 212                                           DDL_2='DX',
 213                                           COEF_IMPO=0.0,
 214                                           COEF_MULT_2=-1.0,
 215                                           DDL_1='DX',
 216                                           COEF_MULT_1=1.0,),
 217                                        _F(GROUP_NO_1='LIAS_1',
 218                                           GROUP_NO_2='LIAS_2',
 219                                           DDL_2='DY',
 220                                           COEF_IMPO=0.0,
 221                                           COEF_MULT_2=-1.0,
 222                                           DDL_1='DY',
 223                                           COEF_MULT_1=1.0,),
 224                                        _F(GROUP_NO_1='LIAS_1',
 225                                           GROUP_NO_2='LIAS_2',
 226                                           DDL_2='DZ',
 227                                           COEF_IMPO=0.0,
 228                                           COEF_MULT_2=-1.0,
 229                                           DDL_1='DZ',
 230                                           COEF_MULT_1=1.0,),
 231                                        _F(GROUP_NO_1='LIAS_1',
 232                                           GROUP_NO_2='LIAS_2',
 233                                           DDL_2='DRY',
 234                                           COEF_IMPO=0.0,
 235                                           COEF_MULT_2=-1.0,
 236                                           DDL_1='DRY',
 237                                           COEF_MULT_1=1.0,),
 238                                        _F(GROUP_NO_1='LIAS_1',
 239                                           GROUP_NO_2='LIAS_2',
 240                                           DDL_2='DRZ',
 241                                           COEF_IMPO=0.0,
 242                                           COEF_MULT_2=-1.0,
 243                                           DDL_1='DRZ',
 244                                           COEF_MULT_1=1.0,),),);
 245 #
 246
 247 CL_CTC=AFFE_CHAR_MECA(MODELE=MO1,
 248                       INFO=1,
 249                       CONTACT=_F(APPARIEMENT='NODAL',
 250                                  VECT_NORM_2=(-1.0,0.0,0.0,),
 251                                  GROUP_MA_1='EBOINF',
 252                                  GROUP_MA_2='MCRBAS',),);
 253 #
 254 # DEPLACEMENT DU SYSTEME DE MAINTIEN
 255
 256 CLDI_F1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 257                             1000.0,0.1,
 258                             2000.0,0.2,
 259                             3000.0,0.2,
 260                             4000.0,0.2,
 261                             5000.0,0.2,
 262                             6000.0,0.1,
 263                             7000.0,0.1,
 264                             ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 265
 266 CL_DI_1=AFFE_CHAR_MECA(MODELE=MO1,
 267                        DDL_IMPO=_F(GROUP_NO='MAIN_SUP',
 268                                    DX=-1.0,),);
 269 # DEFINITION DES FORCES EXTERIEURES
 270 #
 271 # FONCTION MULTIPLICATIVE DE LA FORCE HYDRODYNAMIQUE NORMEE
 272
 273 HYDR_F1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 274                             1000.0,0.1,
 275                             2000.0,0.2,
 276                             3000.0,0.2,
 277                             4000.0,0.2,
 278                             5000.0,0.2,
 279                             6000.0,0.1,
 280                             7000.0,0.1,
 281                             ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 282 #
 283 # PROFIL DES FORCES HYDRODYNAMIQUES TUBES GUIDES A FROID
 284
 285 HYFRTU_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
 286                              5.0,2.0,
 287                              ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',);
 288 #
 289 # PROFIL DES FORCES HYDRODYNAMIQUES CRAYON A FROID
 290
 291 HYFRCR_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
 292                              5.0,2.0,
 293                              ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',);
 294 #
 295 # CHARGEMENT : FORCES HYDRODYNAMIQUES
 296
 297 HYDR_1=AFFE_CHAR_MECA(MODELE=MO1,
 298                       FORCE_COQUE=(_F(FX=10.0,
 299                                       GROUP_MA='EBOSUP',),
 300                                    _F(FX=10.0,
 301                                       GROUP_MA='EBOINF',),),
 302                       FORCE_POUTRE=(_F(GROUP_MA='GRIL10',
 303                                        FX=10.0,),
 304                                     _F(GROUP_MA='GRIL1',
 305                                        FX=10.0,),
 306                                     _F(GROUP_MA='GRIL2',
 307                                        FX=10.0,),
 308                                     _F(GROUP_MA='GRIL3',
 309                                        FX=10.0,),
 310                                     _F(GROUP_MA='GRIL4',
 311                                        FX=10.0,),
 312                                     _F(GROUP_MA='GRIL5',
 313                                        FX=10.0,),
 314                                     _F(GROUP_MA='GRIL6',
 315                                        FX=10.0,),
 316                                     _F(GROUP_MA='GRIL7',
 317                                        FX=10.0,),
 318                                     _F(GROUP_MA='GRIL8',
 319                                        FX=10.0,),
 320                                     _F(GROUP_MA='GRIL9',
 321                                        FX=10.0,),),);
 322
 323 TOTO_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,10.0,
 324                            5.0,10.0,
 325                            ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',);
 326 #
 327
 328 FOHYDR_1=AFFE_CHAR_MECA_F(MODELE=MO1,
 329                           FORCE_POUTRE=(_F(GROUP_MA='GUID',
 330                                            FX=TOTO_1,),
 331                                         _F(GROUP_MA='CRAY',
 332                                            FX=TOTO_1,),),);
 333 #
 334 # PROFIL DES FORCES D'ARCHIMEDE TUBE GUIDE
 335
 336 ARFRTU_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,10.0,
 337                              5.0,10.0,
 338                              ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',);
 339 #
 340 # PROFIL DES FORCES D'ARCHIMEDE CRAYON
 341
 342 ARFRCR_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,10.0,
 343                              5.0,10.0,
 344                              ),NOM_PARA='X',PROL_GAUCHE='CONSTANT',);
 345 #
 346 # FONCTION MULTIPLICATIVE DE LA FORCE D'ARCHIMEDE
 347
 348 ARCH_F1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 349                             1000.0,0.1,
 350                             2000.0,0.2,
 351                             3000.0,0.2,
 352                             4000.0,0.2,
 353                             5000.0,0.2,
 354                             6000.0,0.1,
 355                             7000.0,0.1,
 356                             ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 357 #
 358 #
 359 # CHARGEMENT : FORCES D'ARCHIMEDE
 360 # ----------------------------------------
 361
 362 ARCH_1=AFFE_CHAR_MECA(MODELE=MO1,
 363                       FORCE_COQUE=(_F(FX=10.0,
 364                                       GROUP_MA='EBOSUP',),
 365                                    _F(FX=10.0,
 366                                       GROUP_MA='EBOINF',),),
 367                       FORCE_POUTRE=(_F(GROUP_MA='MAINTIEN',
 368                                        FX=10.0,),
 369                                     _F(GROUP_MA='GRIL10',
 370                                        FX=10.0,),
 371                                     _F(GROUP_MA='GRIL1',
 372                                        FX=10.0,),
 373                                     _F(GROUP_MA='GRIL2',
 374                                        FX=10.0,),
 375                                     _F(GROUP_MA='GRIL3',
 376                                        FX=10.0,),
 377                                     _F(GROUP_MA='GRIL4',
 378                                        FX=10.0,),
 379                                     _F(GROUP_MA='GRIL5',
 380                                        FX=10.0,),
 381                                     _F(GROUP_MA='GRIL6',
 382                                        FX=10.0,),
 383                                     _F(GROUP_MA='GRIL7',
 384                                        FX=10.0,),
 385                                     _F(GROUP_MA='GRIL8',
 386                                        FX=10.0,),
 387                                     _F(GROUP_MA='GRIL9',
 388                                        FX=10.0,),),);
 389
 390 FOARCH_1=AFFE_CHAR_MECA_F(MODELE=MO1,
 391                           FORCE_POUTRE=(_F(GROUP_MA='GUID',
 392                                            FX=ARCH_F1,),
 393                                         _F(GROUP_MA='CRAY',
 394                                            FX=ARCH_F1,),),);
 395 #
 396 #
 397 # FORCE DE SERRAGE DE LA LIAISON GRILLE - CRAYON
 398 #
 399 #
 400 #
 401 #
 402 #
 403 # DEFINITION DE LA LISTE D'INSTANTS
 404
 405 L_FRIG=DEFI_LIST_REEL(DEBUT=0.0,
 406                       INTERVALLE=_F(NOMBRE=100,
 407                                     JUSQU_A=1000.0,),);
 408 #
 409 # INRTERPRETATION DES FONCTIONS
 410
 411 F1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 412                        43200.0,0.0,
 413                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 414
 415 F2=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 416                        43200.0,0.0,
 417                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 418
 419 F3=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 420                        43200.0,0.0,
 421                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 422
 423 F4=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 424                        43200.0,0.0,
 425                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 426
 427 F5=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 428                        43200.0,0.0,
 429                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 430
 431 F6=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 432                        43200.0,0.0,
 433                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 434
 435 F7=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 436                        43200.0,0.0,
 437                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 438
 439 F8=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 440                        43200.0,0.0,
 441                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 442
 443 F9=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 444                        43200.0,0.0,
 445                        ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 446
 447 F10=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 448                         43200.0,0.0,
 449                         ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 450
 451 FRIG1=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 452                        PROL_DROIT='CONSTANT',
 453                        PROL_GAUCHE='CONSTANT',
 454                        INTERPOL='LIN',
 455                        FONCTION=F1,);
 456
 457 FRIG2=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 458                        PROL_DROIT='CONSTANT',
 459                        PROL_GAUCHE='CONSTANT',
 460                        INTERPOL='LIN',
 461                        FONCTION=F2,);
 462
 463 FRIG3=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 464                        PROL_DROIT='CONSTANT',
 465                        PROL_GAUCHE='CONSTANT',
 466                        INTERPOL='LIN',
 467                        FONCTION=F3,);
 468
 469 FRIG4=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 470                        PROL_DROIT='CONSTANT',
 471                        PROL_GAUCHE='CONSTANT',
 472                        INTERPOL='LIN',
 473                        FONCTION=F4,);
 474
 475 FRIG5=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 476                        PROL_DROIT='CONSTANT',
 477                        PROL_GAUCHE='CONSTANT',
 478                        INTERPOL='LIN',
 479                        FONCTION=F5,);
 480
 481 FRIG6=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 482                        PROL_DROIT='CONSTANT',
 483                        PROL_GAUCHE='CONSTANT',
 484                        INTERPOL='LIN',
 485                        FONCTION=F6,);
 486
 487 FRIG7=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 488                        PROL_DROIT='CONSTANT',
 489                        PROL_GAUCHE='CONSTANT',
 490                        INTERPOL='LIN',
 491                        FONCTION=F7,);
 492
 493 FRIG8=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 494                        PROL_DROIT='CONSTANT',
 495                        PROL_GAUCHE='CONSTANT',
 496                        INTERPOL='LIN',
 497                        FONCTION=F8,);
 498
 499 FRIG9=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 500                        PROL_DROIT='CONSTANT',
 501                        PROL_GAUCHE='CONSTANT',
 502                        INTERPOL='LIN',
 503                        FONCTION=F9,);
 504
 505 FRIG10=CALC_FONC_INTERP(LIST_PARA=L_FRIG,
 506                         PROL_DROIT='CONSTANT',
 507                         PROL_GAUCHE='CONSTANT',
 508                         INTERPOL='LIN',
 509                         FONCTION=F10,);
 510 #
 511 # DEFINITION DU FLUX NEUTRONIQUE RADIAL
 512
 513 FLUX_F1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,0.0,
 514                             43200.0,0.0,
 515                             974520.0,453280.0,
 516                             2638800.0,2077850.0,
 517                             4087800.0,3512790.0,
 518                             5428440.0,4856770.0,
 519                             7245000.0,6676270.0,
 520                             8853480.0,8279770.0,
 521                             10460900.0,9872610.0,
 522                             12202900.0,11586900.0,
 523                             13944600.0,13295200.0,
 524                             14143700.0,13490200.0,
 525                             15563200.0,14875300.0,
 526                             17236100.0,16506900.0,
 527                             17948200.0,17199200.0,
 528                             18190100.0,17433400.0,
 529                             18193700.0,17435200.0,
 530                             18197300.0,17435200.0,
 531                             18200900.0,17435200.0,
 532                             ),NOM_PARA='INST',PROL_GAUCHE='CONSTANT',);
 533 # DEFINITION DU FLUX AXIAL POUR LE CYCLE 1
 534
 535 FLUXAX1=LIRE_FONCTION(PROL_DROIT='CONSTANT',
 536                       UNITE=60,
 537                       PROL_GAUCHE='CONSTANT',);
 538
 539 FLY_1_1=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,3.80328347578e-07,
 540                             -0.0882,3.80836894587e-07,
 541                             -0.0756,3.82960826211e-07,
 542                             -0.063,3.85840099715e-07,
 543                             -0.0504,3.89461004273e-07,
 544                             -0.0378,3.92594551282e-07,
 545                             -0.0252,3.92584579772e-07,
 546                             -0.0126,3.92941061254e-07,
 547                             0.0,3.93995548433e-07,
 548                             0.0126,3.92932336182e-07,
 549                             0.0252,3.9256712963e-07,
 550                             0.0378,3.9256712963e-07,
 551                             0.0504,3.89426103989e-07,
 552                             0.063,3.8579522792e-07,
 553                             0.0756,3.82909722222e-07,
 554                             0.0882,3.80777065527e-07,
 555                             0.1008,3.80261039886e-07,
 556                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 557
 558 FLY_1_2=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,4.35955662393e-07,
 559                             -0.0882,4.38665420228e-07,
 560                             -0.0756,4.43314636752e-07,
 561                             -0.063,4.49115562678e-07,
 562                             -0.0504,4.57506588319e-07,
 563                             -0.0378,4.66472222222e-07,
 564                             -0.0252,4.60029380342e-07,
 565                             -0.0126,4.60043091168e-07,
 566                             0.0,4.6686235755e-07,
 567                             0.0126,4.60033119658e-07,
 568                             0.0252,4.60009437322e-07,
 569                             0.0378,4.66439814815e-07,
 570                             0.0504,4.57464209402e-07,
 571                             0.063,4.49064458689e-07,
 572                             0.0756,4.43253561254e-07,
 573                             0.0882,4.38596866097e-07,
 574                             0.1008,4.35878383191e-07,
 575                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 576
 577 FLY_1_3=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,4.94318910256e-07,
 578                             -0.0882,5.00051282051e-07,
 579                             -0.0756,5.13371972934e-07,
 580                             -0.063,5.29153133903e-07,
 581                             -0.0504,5.37880698006e-07,
 582                             -0.0378,5.35046296296e-07,
 583                             -0.0252,5.32211894587e-07,
 584                             -0.0126,5.31364316239e-07,
 585                             0.0,5.31358707265e-07,
 586                             0.0126,5.31353098291e-07,
 587                             0.0252,5.32186965812e-07,
 588                             0.0378,5.35008279915e-07,
 589                             0.0504,5.37829594017e-07,
 590                             0.063,5.29092058405e-07,
 591                             0.0756,5.13302172365e-07,
 592                             0.0882,4.99971509972e-07,
 593                             0.1008,4.94231659544e-07,
 594                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 595
 596 FLY_1_4=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,5.5489957265e-07,
 597                             -0.0882,5.64778846154e-07,
 598                             -0.0756,5.9007781339e-07,
 599                             -0.063,5.97264155983e-07,
 600                             -0.0504,6.04450498576e-07,
 601                             -0.0378,5.98407763533e-07,
 602                             -0.0252,5.84406517094e-07,
 603                             -0.0126,5.82639066952e-07,
 604                             0.0,5.91532407407e-07,
 605                             0.0126,5.82626602564e-07,
 606                             0.0252,5.8438034188e-07,
 607                             0.0378,5.98366631054e-07,
 608                             0.0504,6.04393162393e-07,
 609                             0.063,5.97194978632e-07,
 610                             0.0756,5.89996794872e-07,
 611                             0.0882,5.64691595442e-07,
 612                             0.1008,5.54801103989e-07,
 613                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 614
 615 FLY_1_5=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,6.12672008547e-07,
 616                             -0.0882,6.29623575499e-07,
 617                             -0.0756,6.56554131054e-07,
 618                             -0.063,6.61645833333e-07,
 619                             -0.0504,6.51902421652e-07,
 620                             -0.0378,6.52126780627e-07,
 621                             -0.0252,6.37466168091e-07,
 622                             -0.0126,6.35474358974e-07,
 623                             0.0,6.45096866097e-07,
 624                             0.0126,6.35460648148e-07,
 625                             0.0252,6.37436253561e-07,
 626                             0.0378,6.52079415954e-07,
 627                             0.0504,6.51838853276e-07,
 628                             0.063,6.6156980057e-07,
 629                             0.0756,6.56465633903e-07,
 630                             0.0882,6.29523860399e-07,
 631                             0.1008,6.12563568376e-07,
 632                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 633
 634 FLY_1_6=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,6.57551282051e-07,
 635                             -0.0882,6.83362535613e-07,
 636                             -0.0756,6.90120726496e-07,
 637                             -0.063,6.96878917379e-07,
 638                             -0.0504,6.93625712251e-07,
 639                             -0.0378,6.90747685185e-07,
 640                             -0.0252,6.8786965812e-07,
 641                             -0.0126,6.85864138177e-07,
 642                             0.0,6.85856036325e-07,
 643                             0.0126,6.85847934473e-07,
 644                             0.0252,6.87836004274e-07,
 645                             0.0378,6.90697204416e-07,
 646                             0.0504,6.93558404558e-07,
 647                             0.063,6.9679789886e-07,
 648                             0.0756,6.90024127493e-07,
 649                             0.0882,6.83250356125e-07,
 650                             0.1008,6.57431623932e-07,
 651                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 652
 653 FLY_1_7=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,6.96845263533e-07,
 654                             -0.0882,7.14066061254e-07,
 655                             -0.0756,7.3236502849e-07,
 656                             -0.063,7.20764423077e-07,
 657                             -0.0504,7.17935007123e-07,
 658                             -0.0378,7.28356481481e-07,
 659                             -0.0252,7.1444747151e-07,
 660                             -0.0126,7.11308938746e-07,
 661                             0.0,7.21077279202e-07,
 662                             0.0126,7.11293981481e-07,
 663                             0.0252,7.14412571225e-07,
 664                             0.0378,7.28302884615e-07,
 665                             0.0504,7.17863960114e-07,
 666                             0.063,7.20678418803e-07,
 667                             0.0756,7.32261574074e-07,
 668                             0.0882,7.13951388889e-07,
 669                             0.1008,6.96718126781e-07,
 670                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 671
 672 FLY_1_8=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,7.36793625356e-07,
 673                             -0.0882,7.54231303419e-07,
 674                             -0.0756,7.72156339031e-07,
 675                             -0.063,7.58657407407e-07,
 676                             -0.0504,7.55477742165e-07,
 677                             -0.0378,7.66624643875e-07,
 678                             -0.0252,7.50884615385e-07,
 679                             -0.0126,7.43139245014e-07,
 680                             0.0,7.4863230057e-07,
 681                             0.0126,7.43120548433e-07,
 682                             0.0252,7.50847222222e-07,
 683                             0.0378,7.66567307692e-07,
 684                             0.0504,7.55401709402e-07,
 685                             0.063,7.58562678063e-07,
 686                             0.0756,7.72044159544e-07,
 687                             0.0882,7.54107905983e-07,
 688                             0.1008,7.36657763533e-07,
 689                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 690
 691 FLY_1_9=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,7.73709401709e-07,
 692                             -0.0882,8.01447649573e-07,
 693                             -0.0756,8.03743589744e-07,
 694                             -0.063,8.06039529915e-07,
 695                             -0.0504,8.02383725071e-07,
 696                             -0.0378,7.99366720085e-07,
 697                             -0.0252,7.963497151e-07,
 698                             -0.0126,7.83173611111e-07,
 699                             0.0,7.83164262821e-07,
 700                             0.0126,7.8315491453e-07,
 701                             0.0252,7.96311075499e-07,
 702                             0.0378,7.99306891026e-07,
 703                             0.0504,8.02302706553e-07,
 704                             0.063,8.05938568376e-07,
 705                             0.0756,8.03625801282e-07,
 706                             0.0882,8.01313034188e-07,
 707                             0.1008,7.73563568376e-07,
 708                             ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 709
 710 FLY_1_10=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,7.94241987179e-07,
 711                              -0.0882,8.12259259259e-07,
 712                              -0.0756,8.30755163818e-07,
 713                              -0.063,8.15406517094e-07,
 714                              -0.0504,8.11310719373e-07,
 715                              -0.0378,8.22974893162e-07,
 716                              -0.0252,8.05774038462e-07,
 717                              -0.0126,7.97172364672e-07,
 718                              0.0,8.02846153846e-07,
 719                              0.0126,7.97151175214e-07,
 720                              0.0252,8.05731659544e-07,
 721                              0.0378,8.22912571225e-07,
 722                              0.0504,8.11228454416e-07,
 723                              0.063,8.15303062678e-07,
 724                              0.0756,8.30630519943e-07,
 725                              0.0882,8.12118411681e-07,
 726                              0.1008,7.94087428775e-07,
 727                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 728
 729 FLY_1_11=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.16478454416e-07,
 730                              -0.0882,8.35106481481e-07,
 731                              -0.0756,8.54878739316e-07,
 732                              -0.063,8.39690883191e-07,
 733                              -0.0504,8.35029202279e-07,
 734                              -0.0378,8.46510149573e-07,
 735                              -0.0252,8.29739316239e-07,
 736                              -0.0126,8.2549269943e-07,
 737                              0.0,8.36379095442e-07,
 738                              0.0126,8.2547275641e-07,
 739                              0.0252,8.29693198006e-07,
 740                              0.0378,8.46444088319e-07,
 741                              0.0504,8.34939458689e-07,
 742                              0.063,8.39579950142e-07,
 743                              0.0756,8.54747863248e-07,
 744                              0.0882,8.34956908832e-07,
 745                              0.1008,8.16316417379e-07,
 746                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 747
 748 FLY_1_12=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.3948272792e-07,
 749                              -0.0882,8.70025462963e-07,
 750                              -0.0756,8.76085648148e-07,
 751                              -0.063,8.82145833333e-07,
 752                              -0.0504,8.75892450142e-07,
 753                              -0.0378,8.71288728632e-07,
 754                              -0.0252,8.66685007123e-07,
 755                              -0.0126,8.63216168091e-07,
 756                              0.0,8.63204326923e-07,
 757                              0.0126,8.63192485755e-07,
 758                              0.0252,8.6663764245e-07,
 759                              0.0378,8.71217058405e-07,
 760                              0.0504,8.75796474359e-07,
 761                              0.063,8.82026175214e-07,
 762                              0.0756,8.75946670228e-07,
 763                              0.0882,8.69867165242e-07,
 764                              0.1008,8.39311965812e-07,
 765                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 766
 767 FLY_1_13=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.55083155271e-07,
 768                              -0.0882,8.75466168091e-07,
 769                              -0.0756,9.09445334758e-07,
 770                              -0.063,9.12926638177e-07,
 771                              -0.0504,8.965372151e-07,
 772                              -0.0378,8.9547275641e-07,
 773                              -0.0252,8.74056445869e-07,
 774                              -0.0126,8.7004789886e-07,
 775                              0.0,8.82265491453e-07,
 776                              0.0126,8.70025462963e-07,
 777                              0.0252,8.7400534188e-07,
 778                              0.0378,8.95396723647e-07,
 779                              0.0504,8.96433760684e-07,
 780                              0.063,9.12798254986e-07,
 781                              0.0756,9.09295762108e-07,
 782                              0.0882,8.75302884615e-07,
 783                              0.1008,8.5490491453e-07,
 784                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 785
 786 FLY_1_14=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.57982371795e-07,
 787                              -0.0882,8.68808938746e-07,
 788                              -0.0756,9.02997507123e-07,
 789                              -0.063,9.09486467236e-07,
 790                              -0.0504,9.1597542735e-07,
 791                              -0.0378,9.05003027066e-07,
 792                              -0.0252,8.82140847578e-07,
 793                              -0.0126,8.77798254986e-07,
 794                              0.0,8.89937321937e-07,
 795                              0.0126,8.77772079772e-07,
 796                              0.0252,8.8208974359e-07,
 797                              0.0378,9.04923254986e-07,
 798                              0.0504,9.15865740741e-07,
 799                              0.063,9.09354344729e-07,
 800                              0.0756,9.02842948718e-07,
 801                              0.0882,8.68640669516e-07,
 802                              0.1008,8.57796652422e-07,
 803                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 804
 805 FLY_1_15=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.57735576923e-07,
 806                              -0.0882,8.61808938746e-07,
 807                              -0.0756,8.78628383191e-07,
 808                              -0.063,8.99240740741e-07,
 809                              -0.0504,9.08724893162e-07,
 810                              -0.0378,9.01745459402e-07,
 811                              -0.0252,8.94766025641e-07,
 812                              -0.0126,8.91223646724e-07,
 813                              0.0,8.91209935897e-07,
 814                              0.0126,8.91196225071e-07,
 815                              0.0252,8.94711182336e-07,
 816                              0.0378,9.01661324786e-07,
 817                              0.0504,9.08611467236e-07,
 818                              0.063,8.99103632479e-07,
 819                              0.0756,8.78472578348e-07,
 820                              0.0882,8.61636930199e-07,
 821                              0.1008,8.57543625356e-07,
 822                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 823
 824 FLY_1_16=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.59039351852e-07,
 825                              -0.0882,8.56900462963e-07,
 826                              -0.0756,8.58310185185e-07,
 827                              -0.063,8.61673076923e-07,
 828                              -0.0504,8.71253205128e-07,
 829                              -0.0378,8.85635861823e-07,
 830                              -0.0252,8.7089298433e-07,
 831                              -0.0126,8.68428774929e-07,
 832                              0.0,8.79421118234e-07,
 833                              0.0126,8.68401353276e-07,
 834                              0.0252,8.70838141026e-07,
 835                              0.0378,8.8554985755e-07,
 836                              0.0504,8.71139779202e-07,
 837                              0.063,8.61538461538e-07,
 838                              0.0756,8.58153133903e-07,
 839                              0.0882,8.567247151e-07,
 840                              0.1008,8.58841168091e-07,
 841                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 842
 843 FLY_1_17=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(-0.1008,8.63886752137e-07,
 844                              -0.0882,8.55637820513e-07,
 845                              -0.0756,8.50870192308e-07,
 846                              -0.063,8.47559650997e-07,
 847                              -0.0504,8.47616987179e-07,
 848                              -0.0378,8.51336360399e-07,
 849                              -0.0252,8.48440883191e-07,
 850                              -0.0126,8.46336894587e-07,
 851                              0.0,8.46446581197e-07,
 852                              0.0126,8.46309472934e-07,
 853                              0.0252,8.48384793447e-07,
 854                              0.0378,8.51251602564e-07,
 855                              0.0504,8.47504807692e-07,
 856                              0.063,8.47423789174e-07,
 857                              0.0756,8.50710648148e-07,
 858                              0.0882,8.55455840456e-07,
 859                              0.1008,8.63681089744e-07,
 860                              ),NOM_PARA='Y',PROL_GAUCHE='CONSTANT',);
 861
 862 FLUXRD1=DEFI_NAPPE(PROL_DROIT='CONSTANT',
 863                    PARA=
 864                    (-0.1008,-0.0882,-0.0756,-0.063,-0.0504,-0.0378,-0.0252,-0.0126,0.0,0.0126,0.0252,
 865                           0.0378,0.0504,0.063,0.0756,0.0882,0.1008,),
 866                    NOM_PARA='Z',
 867                    PROL_GAUCHE='CONSTANT',
 868                    FONCTION=
 869                    (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,
 870                               FLY_1_11,FLY_1_12,FLY_1_13,FLY_1_14,FLY_1_15,FLY_1_16,FLY_1_17,),);
 871 #
 872 # DEFINITION DES MATERIAUX ET AFFECTATIONS
 873 #
 874 #           **** DEFINITION DES MATERIAUX SANS FLUX NEUTRONIQUE ****
 875 #
 876 #
 877 # FAMILLE ASSEMBLAGE : NUMERO_FAM
 878 #
 879 # EMBOUT SUPERIEUR
 880
 881 E_ES=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(100.0,2.0,
 882                          300.0,3.0,
 883                          ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
 884
 885 NU_ES=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(100.0,2.0,
 886                           300.0,3.0,
 887                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
 888
 889 AL_ES=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(100.0,2.0,
 890                           300.0,3.0,
 891                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
 892
 893 MT_ES=DEFI_MATERIAU(ELAS_FO=_F(E=E_ES,
 894                                RHO=7800.0,
 895                                ALPHA=AL_ES,
 896                                NU=NU_ES,),);
 897 # EMBOUT INFERIEUR
 898
 899 E_EI=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(100.0,2.0,
 900                          300.0,3.0,
 901                          ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
 902
 903 NU_EI=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(100.0,2.0,
 904                           300.0,3.0,
 905                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
 906
 907 AL_EI=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(100.0,2.0,
 908                           300.0,3.0,
 909                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
 910
 911 MT_EI=DEFI_MATERIAU(ELAS_FO=_F(E=E_EI,
 912                                RHO=12.0,
 913                                ALPHA=AL_EI,
 914                                NU=NU_EI,),);
 915
 916 MT_RIG=DEFI_MATERIAU(ELAS=_F(E=1.0,
 917                              RHO=0.0,
 918                              ALPHA=0.0,
 919                              NU=0.0,),);
 920 # FONCTION COUPLE EN FONCTION DE DRZ POUR LA LIAISON GRILLE/CRAYON
 921
 922 FOTRAC=DEFI_FONCTION(PROL_DROIT='LINEAIRE',VALE=(1.0,1.0,
 923                            0.0,0.0,
 924                            0.0,0.0,
 925                            0.0,0.0,
 926                            ),NOM_PARA='DRZ',PROL_GAUCHE='LINEAIRE',);
 927
 928 MT_RE1=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 929                                     RIGI_N_FO=FRIG1,
 930                                     COULOMB=1.0,
 931                                     RELA_MZ=FOTRAC,),
 932                      ELAS=_F(E=1.0,
 933                              RHO=0.0,
 934                              ALPHA=0.0,
 935                              NU=0.0,),);
 936
 937 MT_RE2=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 938                                     RIGI_N_FO=FRIG2,
 939                                     COULOMB=1.0,
 940                                     RELA_MZ=FOTRAC,),
 941                      ELAS=_F(E=1.0,
 942                              RHO=0.0,
 943                              ALPHA=0.0,
 944                              NU=0.0,),);
 945
 946 MT_RE3=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 947                                     RIGI_N_FO=FRIG3,
 948                                     COULOMB=1.0,
 949                                     RELA_MZ=FOTRAC,),
 950                      ELAS=_F(E=1.0,
 951                              RHO=0.0,
 952                              ALPHA=0.0,
 953                              NU=0.0,),);
 954
 955 MT_RE4=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 956                                     RIGI_N_FO=FRIG4,
 957                                     COULOMB=1.0,
 958                                     RELA_MZ=FOTRAC,),
 959                      ELAS=_F(E=1.0,
 960                              RHO=0.0,
 961                              ALPHA=0.0,
 962                              NU=0.0,),);
 963
 964 MT_RE5=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 965                                     RIGI_N_FO=FRIG5,
 966                                     COULOMB=1.0,
 967                                     RELA_MZ=FOTRAC,),
 968                      ELAS=_F(E=1.0,
 969                              RHO=0.0,
 970                              ALPHA=0.0,
 971                              NU=0.0,),);
 972
 973 MT_RE6=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 974                                     RIGI_N_FO=FRIG6,
 975                                     COULOMB=1.0,
 976                                     RELA_MZ=FOTRAC,),
 977                      ELAS=_F(E=1.0,
 978                              RHO=0.0,
 979                              ALPHA=0.0,
 980                              NU=0.0,),);
 981
 982 MT_RE7=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 983                                     RIGI_N_FO=FRIG7,
 984                                     COULOMB=1.0,
 985                                     RELA_MZ=FOTRAC,),
 986                      ELAS=_F(E=1.0,
 987                              RHO=0.0,
 988                              ALPHA=0.0,
 989                              NU=0.0,),);
 990
 991 MT_RE8=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
 992                                     RIGI_N_FO=FRIG8,
 993                                     COULOMB=1.0,
 994                                     RELA_MZ=FOTRAC,),
 995                      ELAS=_F(E=1.0,
 996                              RHO=0.0,
 997                              ALPHA=0.0,
 998                              NU=0.0,),);
 999
1000 MT_RE9=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
1001                                     RIGI_N_FO=FRIG9,
1002                                     COULOMB=1.0,
1003                                     RELA_MZ=FOTRAC,),
1004                      ELAS=_F(E=1.0,
1005                              RHO=0.0,
1006                              ALPHA=0.0,
1007                              NU=0.0,),);
1008
1009 MT_RE10=DEFI_MATERIAU(DIS_CONTACT=_F(EFFO_N_INIT=-35.0,
1010                                      RIGI_N_FO=FRIG10,
1011                                      COULOMB=1.0,
1012                                      RELA_MZ=FOTRAC,),
1013                       ELAS=_F(E=1.0,
1014                               RHO=0.0,
1015                               ALPHA=0.0,
1016                               NU=0.0,),);
1017 # SYSTEME DE MAINTIEN
1018
1019 E_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1020                          5.0,2.0,
1021                          ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1022
1023 NU_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1024                           5.0,2.0,
1025                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1026
1027 AL_MA=DEFI_CONSTANTE(VALE=0.3,);
1028
1029 SY_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1030                           5.0,2.0,
1031                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1032
1033 DS_MA=DEFI_CONSTANTE(VALE=0.3,);
1034
1035 NP_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1036                           5.0,2.0,
1037                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1038
1039 MEY_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1040                            5.0,2.0,
1041                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1042
1043 MPY_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1044                            5.0,2.0,
1045                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1046
1047 MEZ_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1048                            5.0,2.0,
1049                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1050
1051 MPZ_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1052                            5.0,2.0,
1053                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1054
1055 CAY_MA=DEFI_CONSTANTE(VALE=0.3,);
1056
1057 CBY_MA=DEFI_CONSTANTE(VALE=0.3,);
1058
1059 CAZ_MA=DEFI_CONSTANTE(VALE=0.3,);
1060
1061 CBZ_MA=DEFI_CONSTANTE(VALE=0.3,);
1062
1063 MPX_MA=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1064                            5.0,2.0,
1065                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1066
1067 MT_MA=DEFI_MATERIAU(ELAS_FO=_F(E=E_MA,
1068                                RHO=7800.0,
1069                                ALPHA=AL_MA,
1070                                NU=NU_MA,),
1071                     ECRO_LINE_FO=_F(SY=SY_MA,
1072                                     D_SIGM_EPSI=DS_MA,),
1073                     VMIS_POUTRE_FO=_F(CAZ=CAZ_MA,
1074                                       CAY=CAY_MA,
1075                                       MEY=MEY_MA,
1076                                       NP=NP_MA,
1077                                       CBZ=CBZ_MA,
1078                                       CBY=CBY_MA,
1079                                       MPZ=MPZ_MA,
1080                                       MPY=MPY_MA,
1081                                       MPX=MPX_MA,
1082                                       MEZ=MEZ_MA,),);
1083 #
1084 #           **** DEFINITION DES MATERIAUX AVEC FLUX NEUTRONIQUE ****
1085 #
1086 #
1087 # GRILLE
1088
1089 E_GRE=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1090                           5.0,2.0,
1091                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1092
1093 E_GRM=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1094                           5.0,2.0,
1095                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1096
1097 NU_GRE=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1098                            5.0,2.0,
1099                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1100
1101 NU_GRM=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1102                            5.0,2.0,
1103                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1104
1105 AL_GRE=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1106                            5.0,2.0,
1107                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1108
1109 AL_GRM=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1110                            5.0,2.0,
1111                            ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1112 # TUBE GUIDE
1113
1114 E_TU=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1115                          5.0,2.0,
1116                          ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1117
1118 NU_TU=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1119                           5.0,2.0,
1120                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1121
1122 AL_TU=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1123                           5.0,2.0,
1124                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1125 # CRAYON
1126
1127 E_CR=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1128                          5.0,2.0,
1129                          ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1130
1131 NU_CR=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1132                           5.0,2.0,
1133                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1134
1135 AL_CR=DEFI_FONCTION(PROL_DROIT='CONSTANT',VALE=(0.0,1.0,
1136                           5.0,2.0,
1137                           ),NOM_PARA='TEMP',PROL_GAUCHE='CONSTANT',);
1138 # GRANDISSEMENT DES GRILLES
1139 # GRILLES EXTREMES 1 ET 10
1140
1141 MT_GRE1=DEFI_MATERIAU(LEMAITRE=_F(UN_SUR_K=0.0,
1142                                   QSR_K=1.0,
1143                                   UN_SUR_M=0.0,
1144                                   N=0.0,),
1145                       ELAS_FO=_F(E=E_GRE,
1146                                  RHO=7800.0,
1147                                  ALPHA=AL_GRE,
1148                                  NU=NU_GRE,),
1149                       GRAN_IRRA=_F(A=0.0,
1150                                    FLUX_L=FLUXAX1,
1151                                    FONC_MULT=FLUX_F1,
1152                                    S=1.0,
1153                                    FLUX_TN=FLUXRD1,),);
1154 # GRILLES DE MELANGE
1155
1156 MT_GRM1=DEFI_MATERIAU(LEMAITRE=_F(UN_SUR_K=0.0,
1157                                   QSR_K=2.0,
1158                                   UN_SUR_M=0.0,
1159                                   N=0.0,),
1160                       ELAS_FO=_F(E=E_GRM,
1161                                  RHO=7800.0,
1162                                  ALPHA=AL_GRM,
1163                                  NU=NU_GRM,),
1164                       GRAN_IRRA=_F(A=0.0,
1165                                    FLUX_L=FLUXAX1,
1166                                    FONC_MULT=FLUX_F1,
1167                                    S=1.0,
1168                                    FLUX_TN=FLUXRD1,),);
1169 # DEFINITION DES MATERIAUX VISCOPLASTIQUES ( TUBES ET CRAYONS )
1170
1171 MT_TU1=DEFI_MATERIAU(LEMAITRE=_F(UN_SUR_K=0.0,
1172                                  QSR_K=2.0,
1173                                  UN_SUR_M=0.0,
1174                                  N=0.0,),
1175                      ELAS_FO=_F(E=E_TU,
1176                                 RHO=7800.0,
1177                                 ALPHA=AL_TU,
1178                                 NU=NU_TU,),
1179                      GRAN_IRRA=_F(A=0.0,
1180                                   FLUX_L=FLUXAX1,
1181                                   FONC_MULT=FLUX_F1,
1182                                   S=1.0,
1183                                   FLUX_TN=FLUXRD1,),);
1184
1185 MT_CR1=DEFI_MATERIAU(LEMAITRE=_F(UN_SUR_K=0.0,
1186                                  QSR_K=2.0,
1187                                  UN_SUR_M=0.0,
1188                                  N=0.0,),
1189                      ELAS_FO=_F(E=E_CR,
1190                                 RHO=7800.0,
1191                                 ALPHA=AL_CR,
1192                                 NU=NU_CR,),
1193                      GRAN_IRRA=_F(A=0.0,
1194                                   FLUX_L=FLUXAX1,
1195                                   FONC_MULT=FLUX_F1,
1196                                   S=1.0,
1197                                   FLUX_TN=FLUXRD1,),);
1198 #
1199 #           AFFECTATION DES MATERIAUX A CHAUD - CYCLE
1200
1201 A_MAT_1=AFFE_MATERIAU(MAILLAGE=MA1,
1202                       AFFE=(_F(GROUP_MA=('MAINTI_Y','MAINTI_Z',),
1203                                TEMP_REF=100.0,
1204                                MATER=MT_MA,),
1205                             _F(GROUP_MA='GRIL_I',
1206                                TEMP_REF=100.0,
1207                                MATER=MT_GRM1,),
1208                             _F(GROUP_MA='GRIL_E',
1209                                TEMP_REF=100.0,
1210                                MATER=MT_GRE1,),
1211                             _F(GROUP_MA='CRAY',
1212                                TEMP_REF=100.0,
1213                                MATER=MT_CR1,),
1214                             _F(GROUP_MA='GUID',
1215                                TEMP_REF=100.0,
1216                                MATER=MT_TU1,),
1217                             _F(GROUP_MA='EBOINF',
1218                                TEMP_REF=100.0,
1219                                MATER=MT_EI,),
1220                             _F(GROUP_MA=('PI_PLE','PI_TRO',),
1221                                TEMP_REF=100.0,
1222                                MATER=MT_EI,),
1223                             _F(GROUP_MA='EBOSUP',
1224                                TEMP_REF=100.0,
1225                                MATER=MT_ES,),
1226                             _F(GROUP_MA='RIGID',
1227                                TEMP_REF=100.0,
1228                                MATER=MT_RIG,),
1229                             _F(GROUP_MA='ELA1',
1230                                TEMP_REF=100.0,
1231                                MATER=MT_RE1,),
1232                             _F(GROUP_MA='ELA2',
1233                                TEMP_REF=100.0,
1234                                MATER=MT_RE2,),
1235                             _F(GROUP_MA='ELA3',
1236                                TEMP_REF=100.0,
1237                                MATER=MT_RE3,),
1238                             _F(GROUP_MA='ELA4',
1239                                TEMP_REF=100.0,
1240                                MATER=MT_RE4,),
1241                             _F(GROUP_MA='ELA5',
1242                                TEMP_REF=100.0,
1243                                MATER=MT_RE5,),
1244                             _F(GROUP_MA='ELA6',
1245                                TEMP_REF=100.0,
1246                                MATER=MT_RE6,),
1247                             _F(GROUP_MA='ELA7',
1248                                TEMP_REF=100.0,
1249                                MATER=MT_RE7,),
1250                             _F(GROUP_MA='ELA8',
1251                                TEMP_REF=100.0,
1252                                MATER=MT_RE8,),
1253                             _F(GROUP_MA='ELA9',
1254                                TEMP_REF=100.0,
1255                                MATER=MT_RE9,),
1256                             _F(GROUP_MA='ELA10',
1257                                TEMP_REF=100.0,
1258                                MATER=MT_RE10,),),);
1259 #
1260
1261 L_INST_1=DEFI_LIST_REEL(DEBUT=0.0,
1262                         INTERVALLE=(_F(NOMBRE=1,
1263                                        JUSQU_A=0.5,),
1264                                     _F(NOMBRE=1,
1265                                        JUSQU_A=1000.0,),
1266                                     _F(NOMBRE=1,
1267                                        JUSQU_A=2000.0,),
1268                                     _F(NOMBRE=5,
1269                                        JUSQU_A=3000.0,),
1270                                     _F(NOMBRE=1,
1271                                        JUSQU_A=4000.0,),
1272                                     _F(NOMBRE=1,
1273                                        JUSQU_A=5000.0,),
1274                                     _F(NOMBRE=1,
1275                                        JUSQU_A=6000.0,),
1276                                     _F(NOMBRE=1,
1277                                        JUSQU_A=7000.0,),
1278                                     _F(NOMBRE=1,
1279                                        JUSQU_A=8000.0,),),);
1280
1281 L_ARCH_1=DEFI_LIST_REEL(DEBUT=0.0,
1282                         INTERVALLE=(_F(NOMBRE=1,
1283                                        JUSQU_A=0.5,),
1284                                     _F(NOMBRE=1,
1285                                        JUSQU_A=1000.0,),
1286                                     _F(NOMBRE=1,
1287                                        JUSQU_A=2000.0,),
1288                                     _F(NOMBRE=5,
1289                                        JUSQU_A=3000.0,),
1290                                     _F(NOMBRE=1,
1291                                        JUSQU_A=4000.0,),
1292                                     _F(NOMBRE=1,
1293                                        JUSQU_A=5000.0,),
1294                                     _F(NOMBRE=1,
1295                                        JUSQU_A=6000.0,),
1296                                     _F(NOMBRE=1,
1297                                        JUSQU_A=7000.0,),
1298                                     _F(NOMBRE=1,
1299                                        JUSQU_A=8000.0,),),);
1300 #
1301 # CALCUL PHASE I A V DU CYCLE 1
1302
1303 RESU=STAT_NON_LINE(PARM_THETA=0.5,
1304                    CARA_ELEM=CARA1,
1305                    MODELE=MO1,
1306                    CHAM_MATER=A_MAT_1,
1307                    EXCIT=(_F(CHARGE=ARCH_1,
1308                              FONC_MULT=ARCH_F1,),
1309                           _F(CHARGE=FOARCH_1,
1310                              FONC_MULT=ARCH_F1,),
1311                           _F(CHARGE=HYDR_1,
1312                              FONC_MULT=HYDR_F1,),
1313                           _F(CHARGE=FOHYDR_1,
1314                              FONC_MULT=HYDR_F1,),
1315                           _F(CHARGE=CL_DI_1,
1316                              FONC_MULT=CLDI_F1,),
1317                           _F(CHARGE=CL_PER_1,),
1318                           _F(CHARGE=CL_CTC,),
1319                           _F(CHARGE=PESANT_1,),
1320                           _F(CHARGE=THERM_1,),),
1321                    SOLVEUR=_F(METHODE='MULT_FRONT',
1322                               RENUM='MD',),
1323                    NEWTON=_F(MATRICE='TANGENTE',
1324                              REAC_ITER=1,),
1325                    RECH_LINEAIRE=_F(ITER_LINE_MAXI=3,),
1326                    INCREMENT=_F(LIST_INST=L_INST_1,
1327                                 NUME_INST_FIN=20,),
1328                    COMP_INCR=(_F(GROUP_MA='MAINTIEN',
1329                                  RELATION='VMIS_POU_LINE',),
1330                               _F(GROUP_MA=('CRAYONS','TUB_GUID',),
1331                                  RELATION='ASSE_COMBU',),
1332                               _F(GROUP_MA='LI_GR_CR',
1333                                  RELATION='DIS_CONTACT',),
1334                               _F(GROUP_MA=('EMBO_INF','EMBO_SUP','PIEDS','GRILLES','LI_GR_GU',),
1335                                  RELATION='ELAS',
1336                                  COQUE_NCOU=1,),),
1337                    CONVERGENCE=_F(ARRET='OUI',
1338                                   RESI_GLOB_RELA=0.0001,
1339                                   ITER_INTE_MAXI=500,
1340                                   ITER_GLOB_MAXI=50,),
1341                    ARCHIVAGE=_F(LIST_INST=L_ARCH_1,),);
1342
1343 FIN();