src/TestsEltsVisu/Nicolas/ssnl132a.comm

   1 # MODIF  DATE 07/10/2005   AUTEUR CIBHHPD L.SALMONA
   2 # RESPONSABLE GODARD V.GODARD
   3 # TITRE SIMULATION NUMERIQUE DU COMPORTEMENT D'UN ASSEMBLAGE COMBUSTIBLE
   4 #            CONFIGURATION MANAGEMENT OF EDF VERSION
   5 # ======================================================================
   6 # COPYRIGHT (C) 1991 - 2005  EDF R&D                  WWW.CODE-ASTER.ORG
   7 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
   8 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
   9 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
  10 # (AT YOUR OPTION) ANY LATER VERSION.
  11 #
  12 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  13 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  14 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  15 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  16 #
  17 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
  18 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
  19 #    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
  20 # ======================================================================
  21 #       EN COMPRESSION ET EN FLEXION
  22 #            CONFIGURATION MANAGEMENT OF EDF VERSION
  23 # ======================================================================
  24 # UTILISATION : COMPARAISON DE LA SIMULATION NUMERIQUE D'UN ASSEMBLAGE
  25 #               COMBUSTIBLE EN COMPRESSION AVEC LES RESULTATS EXPERIMENTAUX
  26
  27 DEBUT(CODE=_F( NOM = 'SSNL132A',NIV_PUB_WEB='INTRANET'),
  28       PAR_LOT='NON');
  29
  30
  31 #####################################################
  32 #####################################################
  33 ##     DONNEES GENERALES DU PROBLEME               ##
  34 #####################################################
  35 #####################################################
  36
  37 ###########################################
  38 #        Caracteristiques materiaux       #
  39 ###########################################
  40
  41 # embout superieur
  42 EES0= 1.99e+11;
  43 EES1= 1.75e+11;
  44 NUES0= 0.3;
  45 NUES1= 0.3;
  46 Kesup= 1.5E+8;
  47 RHOES= 898162.6 ;
  48 # embout inferieur
  49 EEI0= 1.99e+11 ;
  50 EEI1= 1.75e+11 ;
  51 NUEI0= 0.3 ;
  52 NUEI1= 0.3 ;
  53 Keinf= 2.6E+8;
  54 RHOEI= 351651.1 ;
  55 # Tube guide
  56 ETU0= 9.84e+10 ;
  57 ETU1= 7.80e+10 ;
  58 NUTU0= 0.3 ;
  59 RHOTU= 6310. ;
  60 # Crayon
  61 ECR0= 9.84e+10 ;
  62 ECR1= 7.80e+10 ;
  63 NUCR0= 0.3 ;
  64 RHOCR= 38050. ;
  65 # systeme de maintien
  66 EMA0= 2.0E11 ;
  67 EMA1= 1.84E11 ;
  68 NUMA0= 0.3 ;
  69 NUMA1= 0.3 ;
  70 MPZMA0= 349. ;
  71 SYMA0= 1.03E9 ;
  72 DSMA0= 0.0 ;
  73 CAYMA0= 0.84 ;
  74 CAZMA0= 0.84 ;
  75 CBYMA0= 0.105 ;
  76 CBZMA0= 0.105 ;
  77 RHOMA = 0. ;
  78 UELA  = 0.105 ;
  79 FELA  = 1480. ;
  80 #####################################
  81 #   Caracteristiques geometriques      #
  82 #####################################
  83 # Caracteristiques geometriques des tubes-guides (longueur,rayon,epaisseur)
  84 LONTU = 4.5445 ;
  85 XINFT = 0.140 ;
  86 #  Partie courante
  87 RAY1GU = 0.006225 ;
  88 EP1GU = 0.0005 ;
  89 #  Partie retreinte
  90 RAY2GU = 0.006225 ;
  91 EP2GU = 0.00118 ;
  92 # Caracteristiques geometriques des crayons (longueur,rayon, epaisseur)
  93 LONCR = 4.4973;
  94 XINFC = 0.1660;
  95 RAYCRA = 0.00475;
  96 EPCRA = 0.00057;
  97 # longueur de la colonne fissile d'un crayon pour profil axial flux
  98 L_COLON    = 4.2672;
  99 L_BOUCHI = 0.0145;
 100 # Caracteristiques geometriques du systeme de maintien
 101 # composante ressort dans le plan horiz lgres (cf fichier fichier maillage gibi
 102 LGRES= 0.1165;
 103 ## hauteur de ressort de maintien comme dans gibi
 104 HRESS= 0.073;
 105 #longueur des embouts
 106 Leinf = 0.14;
 107 Lesup = 0.0988;
 108
 109 ##############################################################
 110 # CARACTERISTIQUES GEOMETRIQUES EMBOUTS                      #
 111 ##############################################################
 112 Heinf = ((Leinf*Keinf)/EEI0)**0.5;
 113 Hesup = ((Lesup*Kesup)/EES0)**0.5;
 114 ##############################################################
 115 # AFFECTATION DES CARACTERISTIQUES GEOMETRIQUES ELEMENTAIRES #
 116 ##############################################################
 117 #
 118 ##############################################################
 119 # TUBES-GUIDES
 120 # PARTIE COURANTE
 121 S_TG_C=3.1415*(RAY1GU**2-(RAY1GU-EP1GU)**2);
 122 I_TG_C=3.1415/4*(RAY1GU**4-(RAY1GU-EP1GU)**4);
 123 # PARTIE RETREINTE
 124 S_TG_R=3.1415*(RAY2GU**2-(RAY2GU-EP2GU)**2);
 125 I_TG_R=3.1415/4*(RAY2GU**4-(RAY2GU-EP2GU)**4);
 126 # LAISON PARTIE COURANTE-RETREINT : BIAIS
 127 # on fait une aproximation en prenant le d moyen sur la hauteur du cone
 128 EPMOY=(EP1GU+EP2GU)/2;
 129 S_TG_B=3.1415*(RAY2GU**2-(RAY2GU-EPMOY)**2);
 130 I_TG_B=3.1415/4*(RAY2GU**4-(RAY2GU-EPMOY)**4);
 131 #
 132 # CRAYONS
 133 S_CR=3.1415*(RAYCRA**2-(RAYCRA-EPCRA)**2);
 134 I_CR=3.1415/4*(RAYCRA**4-(RAYCRA-EPCRA)**4);
 135
 136 #####################################
 137 #   Rigidites des composants         #
 138 #####################################
 139 ########    Grilles de melanges  ########
 140 # Rigidite en rotation des liaisons grille-tube guide
 141 KR_GM = 990. ;
 142 # NOMBRE DE TUBES-GUIDES
 143 NBTG=25;
 144 # NOMBRE DE CRAYONS
 145 NBCR=264.;
 146 # TOTAL
 147 NBTU =289.;
 148 #
 149 ########################################
 150 # CARACTERISTIQUES DU SYST DE MAINTIEN #
 151 ########################################
 152 LSYMA=((LGRES*LGRES)+(HRESS*HRESS))**(1/2) ;
 153 HYMA= (4.*SYMA0*LGRES*LSYMA)/(6.*EMA0*UELA) ;
 154 HZMA= (6.*LGRES*FELA)/(SYMA0*HYMA*HYMA) ;
 155 MEZMA0= SYMA0*HZMA*(HYMA*HYMA)/6. ;
 156 NPMA0= HYMA*HZMA*SYMA0 ;
 157 MEYMA0= SYMA0*HYMA*HZMA*HZMA/6.;
 158 MPYMA0= 1.5*MEYMA0 ;
 159 MPXMA0= 100.*NPMA0 ;
 160 SYMA1=  (SYMA0*EMA1)/EMA0  ;
 161 NPMA1=  (NPMA0*EMA1)/EMA0  ;
 162 MEYMA1= (MEYMA0*EMA1)/EMA0 ;
 163 MPYMA1= (MPYMA0*EMA1)/EMA0 ;
 164 MEZMA1= (MEZMA0*EMA1)/EMA0 ;
 165 MPZMA1= (MPZMA0*EMA1)/EMA0 ;
 166 MPXMA1= (MPXMA0*EMA1)/EMA0 ;
 167 #
 168
 169 #####################################################
 170 #####################################################
 171 ##            COURBES DE RECALAGE                  ##
 172 #####################################################
 173 #####################################################
 174
 175 # Effort axial max sur embout Sup
 176 EFF_COMP = 40000. ;
 177 # Definition de la courbe de reference essai de compression
 178 # charge en compression
 179 COMP_CH=DEFI_FONCTION(
 180    NOM_PARA='DX',
 181    VALE=(
 182          0.        ,0.,
 183          2.5E-04   ,1.0331E+04,
 184          6.63E-04  ,2.5E+04,
 185          1.E-03    ,3.4702E+04,
 186          1.23E-03  ,4.E+04,
 187          ),
 188    PROL_DROITE='CONSTANT',
 189    PROL_GAUCHE='CONSTANT',);
 190 # decharge en compression
 191 COMP_DCH=DEFI_FONCTION(
 192    NOM_PARA='DX',
 193    VALE=(
 194          2.35E-04,   0.,
 195          3.87E-04,   5000.,
 196          5.18E-04,   1.0331E+04,
 197          8.49E-04,   2.5E+04,
 198          1.09E-03,   3.5E+04,
 199          1.23E-03,   4.E+04,
 200    ),
 201    PROL_DROITE='CONSTANT',
 202    PROL_GAUCHE='CONSTANT',);
 203
 204 # Definition des courbes de reference essai de flexion
 205 # Effort axial constant sur embout Sup
 206 EFF_FLEX = 9000. ;
 207 # Charge
 208
 209 #####################################################
 210 #####################################################
 211 ##    MAILLAGE,MODELE,CHARGEMENT,MATERIAU            ##
 212 #####################################################
 213 #####################################################
 214
 215
 216 ###########################################
 217 #        LECTURE DU MAILLAGE              #
 218 ###########################################
 219 #
 220 MA1=LIRE_MAILLAGE();
 221 #
 222 #
 223 MA1=DEFI_GROUP(
 224    reuse =MA1,
 225    MAILLAGE=MA1,
 226    CREA_GROUP_NO=(_F(TOUT_GROUP_MA='OUI',),),
 227               );
 228
 229 ###########################################
 230 #        AFFECTATION DU MODELE            #
 231 ###########################################
 232 #
 233 MO1=AFFE_MODELE(MAILLAGE=MA1,
 234  AFFE=(_F(GROUP_MA=('PQSUP'),
 235           PHENOMENE='MECANIQUE',
 236           MODELISATION='DKT',),
 237        _F(GROUP_MA='CRAYON',
 238           PHENOMENE='MECANIQUE',
 239           MODELISATION='POU_D_TGM',),
 240        _F(GROUP_MA='T_GUIDE',
 241           PHENOMENE='MECANIQUE',
 242           MODELISATION='POU_D_TGM',),
 243        _F(GROUP_MA=('EBOSUP','EBOINF'),
 244           PHENOMENE='MECANIQUE',
 245           MODELISATION='POU_D_T',),
 246        _F(GROUP_MA='MAINTIEN',
 247           PHENOMENE='MECANIQUE',
 248           MODELISATION='POU_D_T',),
 249        _F(GROUP_MA=('ELA','RIG'),
 250           PHENOMENE='MECANIQUE',
 251           MODELISATION='DIS_TR',),
 252           ),);
 253
 254 ##########################################################
 255 #        AFFECTATION DES CARACTERISTIQUES DES ELEMENTS   #
 256 ##########################################################
 257 #
 258 carel=[0.]*78;
 259 #
 260
 261 CARA1=AFFE_CARA_ELEM(
 262    MODELE=MO1,
 263    POUTRE=(
 264      _F(GROUP_MA='CRAYON',
 265         SECTION='GENERALE',
 266         CARA=('A','IZ','IY','JX','AY','AZ','EY','EZ'),
 267         VALE=(S_CR,I_CR,I_CR,I_CR*2,1,1,0,0),),
 268      _F(GROUP_MA='LGUIDE',
 269         SECTION='GENERALE',
 270         CARA=('A','IZ','IY','JX','AY','AZ'),
 271         VALE=(S_TG_C,I_TG_C,I_TG_C,I_TG_C*2,1,1),),
 272      _F(GROUP_MA='RETRE',
 273         SECTION='GENERALE',
 274         CARA=('A','IZ','IY','JX','AY','AZ'),
 275    VALE=(S_TG_R,I_TG_R,I_TG_R,I_TG_R*2,1,1),),
 276 _F(GROUP_MA='BIAIS',
 277         SECTION='GENERALE',
 278         CARA=('A','IZ','IY','JX','AY','AZ'),
 279         VALE=(S_TG_B,I_TG_B,I_TG_B,I_TG_B*2,1,1),),
 280      _F(GROUP_MA=('MAINTI_Y','MAINTI_Z'),
 281         SECTION='RECTANGLE',
 282         CARA=('HY','HZ'),
 283         VALE=(HYMA,HZMA),),
 284      _F(GROUP_MA=('EBOINF'),
 285         SECTION='RECTANGLE',
 286         CARA='H',
 287         VALE= Heinf,),
 288      _F(GROUP_MA=('EBOSUP'),
 289         SECTION='RECTANGLE',
 290         CARA='H',
 291         VALE= Hesup,),
 292            ),
 293    COQUE=(
 294       _F(GROUP_MA=('PQSUP'),
 295          COQUE_NCOU=1,
 296          EPAIS=1.,
 297          ANGL_REP=(90.,0.,),
 298           ),
 299            ),
 300    AFFE_FIBRE =(
 301           _F(GROUP_MA='CRAYON',
 302           COOR_AXE_POUTRE = (0.,0.,),
 303           CARA='SURFACE',
 304           VALE=(
 305              0.,(RAYCRA-EPCRA/2),S_CR/4,
 306              0.,-(RAYCRA-EPCRA/2),S_CR/4,
 307              (RAYCRA-EPCRA/2),0.,S_CR/4,
 308             -(RAYCRA-EPCRA/2),0.,S_CR/4,
 309                 )
 310                ),
 311           _F(GROUP_MA='LGUIDE',
 312                COOR_AXE_POUTRE = (0.,0.,),
 313                CARA='SURFACE',
 314                VALE=(
 315                 0.,(RAY1GU-EP1GU/2),S_TG_C/4,
 316                 0.,-(RAY1GU-EP1GU/2),S_TG_C/4,
 317                 (RAY1GU-EP1GU/2),0.,S_TG_C/4,
 318                -(RAY1GU-EP1GU/2),0.,S_TG_C/4,
 319                     )
 320                ),
 321           _F(GROUP_MA='BIAIS',
 322                COOR_AXE_POUTRE = (0.,0.,),
 323                 CARA='SURFACE',
 324                 VALE=(
 325                 0.,(RAY2GU-EPMOY/2),S_TG_B/4,
 326                 0.,-(RAY2GU-EPMOY/2),S_TG_B/4,
 327                 (RAY2GU-EPMOY/2),0.,S_TG_B/4,
 328                -(RAY2GU-EPMOY/2),0.,S_TG_B/4,
 329                     )
 330                ),
 331           _F(GROUP_MA='RETRE',
 332                COOR_AXE_POUTRE = (0.,0.,),
 333                 CARA='SURFACE',
 334                 VALE=(
 335                 0.,(RAY2GU-EP2GU/2),S_TG_R/4,
 336                 0.,-(RAY2GU-EP2GU/2),S_TG_R/4,
 337                 (RAY2GU-EP2GU/2),0.,S_TG_R/4,
 338                -(RAY2GU-EP2GU/2),0.,S_TG_R/4,
 339                     )
 340                ),
 341                     ),
 342      DISCRET=(
 343                _F(GROUP_MA='ELA_ME',
 344                     REPERE='LOCAL',
 345                     CARA = 'K_TR_L',
 346                     VALE = carel,),
 347                _F(GROUP_MA='ELA_EX',
 348                     REPERE='LOCAL',
 349                     CARA = 'K_TR_L',
 350                     VALE = carel,),
 351                _F(GROUP_MA='RIG_ME',
 352                     REPERE='LOCAL',
 353                     CARA='K_TR_D_L',
 354                     VALE=((1.E9*NBTG/4),(1.E9*NBTG/4),0.,(KR_GM*NBTG/8.),(1.E9*NBTG/4),(KR_GM*NBTG/8.),),),
 355                _F(GROUP_MA='RIG_EX',
 356                     REPERE='LOCAL',
 357                     CARA='K_TR_D_L',
 358                     VALE=((1.E9*NBTG/4),(1.E9*NBTG/4),0.,(KR_GM*NBTG/8.),(1.E9*NBTG/4),(KR_GM*NBTG/8.),),),),
 359
 360      ORIENTATION=(
 361                _F(GROUP_MA=('ELA_EX','ELA_ME','RIG_EX','RIG_ME'),
 362                   CARA='VECT_Y',
 363                   VALE=(1.,0.,0.),),
 364                _F(GROUP_MA=('MAINTI_Y'),
 365                   CARA='VECT_Y',
 366                   VALE=(0. ,1. ,0.),),
 367                _F(GROUP_MA=('MAINTI_Z'),
 368                   CARA='VECT_Y',
 369                   VALE=(0. ,0. ,1. ),),
 370                     ),
 371                );
 372
 373 MA1=MODI_MAILLAGE(
 374    reuse =MA1,
 375    MAILLAGE=MA1,
 376    MODELE=MO1,
 377    ORIE_NORM_COQUE=
 378       _F(   GROUP_MA='PQSUP',
 379             VECT_NORM=(-1.,0.,0.),
 380             GROUP_NO='PQ1',),
 381   INFO=1,);
 382
 383 #####################################################
 384 #              DEFINITION DES MATERIAUX             #
 385 #####################################################
 386 #
 387 #     EMBOUT SUPERIEUR
 388 MT_ES=DEFI_MATERIAU(ELAS=_F(E=EES0,
 389                                NU=NUES0,
 390                                RHO=RHOES,
 391                     ),);
 392
 393 #   EMBOUT INFERIEUR
 394 MT_EI=DEFI_MATERIAU(ELAS=_F(E=EEI0,
 395                                NU=NUEI0,
 396                                RHO=RHOEI,
 397                     ),);
 398 MT_RIG=DEFI_MATERIAU(ELAS=_F(E=1.0,
 399                      NU=0.0,
 400                      RHO=0.0,
 401                      ALPHA=0.0,),);
 402
 403 #    LIAISON GRILLE/CRAYON
 404 kn  =37000.*NBCR;
 405 kt  =140000.*NBCR  ;
 406 rn  = 3700.*NBCR;
 407 rt  =70000.*NBCR  ;
 408 lb  =0.04;
 409 fno  =26.*NBCR/1.25   ;
 410 coul_bos  =0.51  ;
 411 coul_res  =coul_bos*rt/2./kt  ;
 412 ecrob  =kt/1e9 ;
 413 ecror  =rt/1e9  ;
 414
 415 M_GM=DEFI_MATERIAU(DIS_GRICRA=_F(KN_BOS=kn,
 416                                  KT_BOS=kt,
 417                                  KN_RES=rn,
 418                                  KT_RES= rt,
 419                                  DIST_BOS= lb,
 420                                  FORC_SER= fno,
 421                                  COUL_BOS=coul_bos,
 422                                  COUL_RES=coul_res,
 423                                  ECRO_BOS =ecrob,
 424                                  ECRO_RES =ecror,
 425                                    ),
 426                    ELAS=_F( E = 1.E5, NU = 0.3,),
 427                                     );
 428
 429 #    SYSTEME DE MAINTIEN
 430 MT_MA=DEFI_MATERIAU(ELAS=_F(E=EMA0,
 431                                NU=NUMA0,
 432                                RHO=RHOMA,),
 433                     ECRO_LINE=_F(D_SIGM_EPSI=DSMA0,
 434                                     SY=SYMA0,),
 435                     VMIS_POUTRE=_F(NP=NPMA0,
 436                                       MEY=MEYMA0,
 437                                       MPY=MPYMA0,
 438                                       CAY=CAYMA0,
 439                                       CBY=CBYMA0,
 440                                       MEZ=MEZMA0,
 441                                       MPZ=MPZMA0,
 442                                       CAZ=CAZMA0,
 443                                       CBZ=CBZMA0,
 444                                       MPX=MPXMA0,),);
 445
 446 #     TUBE GUIDE
 447 MT_TU1=DEFI_MATERIAU(
 448    ELAS=_F(
 449       E     =  (ETU0*NBTG),
 450       NU    =  NUTU0,
 451       RHO   =  (RHOTU*NBTG),),
 452       );
 453
 454 #     CRAYON
 455 MT_CR=DEFI_MATERIAU(
 456    ELAS=_F(
 457       E     =  (ECR0*NBCR),
 458       NU    =  NUCR0,
 459       RHO   =  (RHOCR*NBCR),
 460        ),
 461       );
 462
 463 #     PQSUP
 464 MT_PLAQ=DEFI_MATERIAU(
 465    ELAS=_F(E=1.E14,
 466            NU=0.3,
 467            RHO=0.,
 468            ALPHA=0.,),);
 469
 470 ##############################################################
 471 #           AFFECTATION DES MATERIAUX                        #
 472 ##############################################################
 473 A_MAT=AFFE_MATERIAU(MAILLAGE=MA1,
 474                     AFFE=(_F(GROUP_MA=('PQSUP'),
 475                              MATER=MT_PLAQ,),
 476                           _F(GROUP_MA=('MAINTIEN'),
 477                              MATER=MT_MA,),
 478                           _F(GROUP_MA='CRAYON',
 479                              MATER=MT_CR,),
 480                           _F(GROUP_MA='T_GUIDE',
 481                              MATER=MT_TU1,),
 482                           _F(GROUP_MA='EBOINF',
 483                              MATER=MT_EI,),
 484                           _F(GROUP_MA=('EBOSUP'),
 485                              MATER=MT_ES,),
 486                           _F(GROUP_MA='RIG',
 487                              MATER=MT_RIG,),
 488                           _F(GROUP_MA='ELA_EX_B',
 489                              MATER=M_GM,),
 490                           _F(GROUP_MA='ELA_EX_T',
 491                              MATER=M_GM,),
 492                           _F(GROUP_MA='ELA_ME',
 493                              MATER=M_GM,),
 494                             ),
 495                             );
 496
 497 ####################################################
 498 #   DEFINITION DES CHARGEMENTS ET DES COND. LIM.   #
 499 ####################################################
 500
 501 # CONDITIONS LIMITES
 502 # ---------------------
 503 CL_CTC=AFFE_CHAR_MECA(
 504    MODELE=MO1,
 505    VERI_NORM='NON',
 506    CONTACT=(_F(APPARIEMENT='NODAL',
 507               VECT_NORM_ESCL=(1.,0.,0.),
 508               GROUP_MA_ESCL='PINF',
 509               GROUP_MA_MAIT='CR_BAS',)
 510               ,
 511            _F(APPARIEMENT='NODAL',
 512               VECT_NORM_ESCL=(-1.,0.,0.),
 513               GROUP_MA_ESCL='PSUP',
 514               GROUP_MA_MAIT='CR_HAU',),),
 515    );
 516
 517
 518 CL_PERM=AFFE_CHAR_MECA(MODELE=MO1,
 519                        DDL_IMPO=(_F(GROUP_MA='PQSUP',
 520                                     DX=0.,   DY=0.,   DZ=0.,),
 521                                  _F(GROUP_MA='ASSEM2',
 522                                     DY=0.,),
 523                                  _F(GROUP_NO='FIX',
 524                                     DX=0.,
 525                                     DZ=0.,
 526                                     DRX=0.,
 527                                       ),
 528                                  _F(GROUP_NO='CRAYON',
 529                                     DRX=0.,
 530                                       ),
 531                                  _F(GROUP_MA=('EBOSUP','EBOINF'),
 532                                     DZ=0.0,
 533                                     DRX=0.0,
 534                                     ),
 535                                  _F(GROUP_NO=('PGRIL_1','PGRIL_2','PGRIL_3','PGRIL_4',
 536                                            'PGRIL_5','PGRIL_6','PGRIL_7','PGRIL_8',
 537                                            'PGRIL_9','PGRIL_10',),
 538                                     DY=0.,DRX=0.,DRY=0.,DRZ=0.),
 539                                     ),
 540                       LIAISON_UNIF=(
 541                                     _F(GROUP_MA='GRIL_1',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 542                                     _F(GROUP_MA='GRIL_2',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 543                                     _F(GROUP_MA='GRIL_3',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 544                                     _F(GROUP_MA='GRIL_4',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 545                                     _F(GROUP_MA='GRIL_5',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 546                                     _F(GROUP_MA='GRIL_6',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 547                                     _F(GROUP_MA='GRIL_7',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 548                                     _F(GROUP_MA='GRIL_8',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 549                                     _F(GROUP_MA='GRIL_9',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 550                                     _F(GROUP_MA='GRIL_10',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 551                                     ),
 552                                     );
 553
 554
 555 BLOQ=AFFE_CHAR_MECA(MODELE=MO1,
 556                     DDL_IMPO=_F(GROUP_NO=('PGRIL_1','PGRIL_2','PGRIL_3','PGRIL_4',
 557                                           'PGRIL_5','PGRIL_6','PGRIL_7','PGRIL_8',
 558                                           'PGRIL_9','PGRIL_10',),
 559                                 DX=0.,),
 560                     );
 561
 562
 563 #
 564 #      CHARGEMENT
 565 # ---------------------
 566 PESANT=AFFE_CHAR_MECA(MODELE=MO1,
 567                       PESANTEUR=(9.81,-1.,0.,0.),);
 568
 569
 570 # Compression initiale
 571 TIR_CR=AFFE_CHAR_MECA(
 572    MODELE=MO1,
 573    DDL_IMPO=
 574    _F(GROUP_NO='CR_BAS',DX=-1.E-3,),);
 575 DEPL_CR=DEFI_FONCTION(
 576    NOM_PARA='INST',
 577    VALE=(
 578      -2.,  0. ,
 579      -1.,  1. ,),);
 580
 581
 582 CL_DI_FR=AFFE_CHAR_MECA(
 583       MODELE=MO1,
 584       FORCE_NODALE=
 585          _F(GROUP_NO='PEBO_S',
 586             FX=(-1.),),);
 587
 588 CL_GRI_6=AFFE_CHAR_MECA(
 589       MODELE=MO1,
 590       DDL_IMPO=
 591          _F(GROUP_NO='PGRIL_6',
 592             DZ=1.E-3,),);
 593
 594 FOR_MAIN=DEFI_FONCTION(NOM_PARA='INST',
 595                        VALE=(
 596                   0.,   0. ,
 597                   1.,   1.e-10 ,
 598                   2.,   1.0331E+04,
 599                   3.,   2.5E+04,
 600                   4.,   3.4702E+04,
 601                   5.,   4.E+04,
 602                   6.,   3.5E+04,
 603                   7.,   2.5E+04,
 604                   8.,   1.0331E+04,
 605                   9.,   5000.,
 606                   10.,   0.
 607                   ),
 608                   PROL_DROITE='CONSTANT',
 609                   PROL_GAUCHE='CONSTANT',);
 610
 611
 612
 613 ###############################################
 614 #  DEFINITION DES LISTES D INSTANTS           #
 615 ###############################################
 616
 617 L_INST1=DEFI_LIST_REEL(
 618    DEBUT=-2.,
 619    INTERVALLE=(
 620       _F(JUSQU_A=-1.,   NOMBRE=10,),
 621                ),);
 622
 623 L_INST11=DEFI_LIST_REEL(
 624    DEBUT=-1.,
 625    INTERVALLE=(
 626       _F(JUSQU_A=0.,   NOMBRE=10,),
 627                ),);
 628
 629 L_INST2=DEFI_LIST_REEL(DEBUT=0.0,
 630                       INTERVALLE=(_F(JUSQU_A=1.0,
 631                                      NOMBRE=2,),
 632                                   _F(JUSQU_A=10.0,
 633                                      NOMBRE=45,),),);
 634
 635 L_INST=DEFI_LIST_REEL(DEBUT=-2.0,
 636                       INTERVALLE=(_F(JUSQU_A=-1.,
 637                                        NOMBRE=10,),
 638                                   _F(JUSQU_A=0.,
 639                                        NOMBRE=10,),
 640                                   _F(JUSQU_A=1.0,
 641                                        NOMBRE=2,),
 642                                   _F(JUSQU_A=10.0,
 643                                        NOMBRE=45,),
 644                                   ),);
 645
 646 #####################################################################
 647 #####################################################################
 648 #     CALCUL INITIAL COMMUN A L ESSAI DE COMPRESSION ET DE FLEXION  #
 649 #####################################################################
 650 #####################################################################
 651
 652 RESU=STAT_NON_LINE(
 653                    MODELE=MO1,
 654                    CHAM_MATER=A_MAT,
 655                    CARA_ELEM=CARA1,
 656                    EXCIT=(_F(CHARGE=TIR_CR, FONC_MULT=DEPL_CR,),
 657                           _F(CHARGE=BLOQ,),
 658                           _F(CHARGE=CL_PERM,),
 659                           _F(CHARGE=CL_CTC,),
 660                           _F(CHARGE=PESANT,),
 661                           ),
 662                    COMP_INCR=(
 663                           _F(RELATION='VMIS_POU_LINE',
 664                              GROUP_MA='MAINTIEN',
 665                              ITER_INTE_MAXI=500),
 666                           _F(RELATION='ELAS',
 667                              GROUP_MA=('CRAYON','T_GUIDE'),),
 668                           _F(RELATION='DIS_GRICRA',
 669                              ITER_INTE_MAXI=500,
 670                              GROUP_MA='ELA'),
 671                           _F(RELATION='ELAS',
 672                              GROUP_MA=('EBOINF','EBOSUP',
 673                                        'RIG','PQSUP'),),),
 674                    INCREMENT=_F(
 675                                  LIST_INST=L_INST1,
 676                                  SUBD_PAS=10,
 677                                  SUBD_PAS_MINI=1.E-6,
 678                                  COEF_SUBD_PAS_1=0.0001,
 679                                  ),
 680                    NEWTON=_F(MATRICE='TANGENTE',
 681                              REAC_ITER=1,),
 682                    CONVERGENCE=_F(RESI_GLOB_RELA=1.e-06,
 683                                   ITER_GLOB_MAXI=50,),
 684                    ARCHIVAGE=_F(LIST_INST=L_INST1,),);
 685
 686 RESU=STAT_NON_LINE(
 687                    reuse = RESU,
 688                    MODELE=MO1,
 689                    CHAM_MATER=A_MAT,
 690                    CARA_ELEM=CARA1,
 691                    EXCIT=(_F(CHARGE=CL_PERM,),
 692                           _F(CHARGE=CL_CTC,),
 693                           _F(CHARGE=PESANT,),
 694                           ),
 695                    COMP_INCR=(
 696                           _F(RELATION='VMIS_POU_LINE',
 697                              GROUP_MA='MAINTIEN',
 698                              ITER_INTE_MAXI=500,),
 699                           _F(RELATION='ELAS',
 700                              GROUP_MA=('CRAYON','T_GUIDE'),),
 701                           _F(RELATION='DIS_GRICRA',
 702                              GROUP_MA='ELA',
 703                              ITER_INTE_MAXI=500,),
 704                           _F(RELATION='ELAS',
 705                              GROUP_MA=('EBOINF','EBOSUP',
 706                                        'RIG','PQSUP'),),),
 707                    INCREMENT=_F(
 708                                 LIST_INST=L_INST11,
 709                                 SUBD_PAS=10,
 710                                 SUBD_PAS_MINI=1.E-6,
 711                                 COEF_SUBD_PAS_1=0.0001,),
 712                    ETAT_INIT=_F(
 713                                 EVOL_NOLI=RESU,
 714                                 INST=-1.,),
 715                    NEWTON=_F(MATRICE='TANGENTE',
 716                                 REAC_ITER=1,),
 717                    CONVERGENCE=_F(RESI_GLOB_RELA=1.e-06,
 718                                   ITER_GLOB_MAXI=50,),
 719                    ARCHIVAGE=_F(LIST_INST=L_INST11,),);
 720
 721
 722 #####################################################################
 723 #####################################################################
 724 #     POURSUITE ESSAI DE COMPRESSION                                #
 725 #####################################################################
 726 #####################################################################
 727
 728 # CALCUL
 729
 730 COMPR=STAT_NON_LINE(
 731                     MODELE=MO1,
 732                     CHAM_MATER=A_MAT,
 733                     CARA_ELEM=CARA1,
 734                     EXCIT=(_F(CHARGE=CL_DI_FR,FONC_MULT=FOR_MAIN,),
 735                            _F(CHARGE=CL_PERM,),
 736                            _F(CHARGE=CL_CTC,),
 737                            _F(CHARGE=PESANT,),
 738                            ),
 739                     COMP_INCR=(
 740                            _F(RELATION='VMIS_POU_LINE',
 741                               GROUP_MA='MAINTIEN',
 742                               ITER_INTE_MAXI=500),
 743                            _F(RELATION='ELAS',
 744                               GROUP_MA=('CRAYON','T_GUIDE'),),
 745                            _F(RELATION='DIS_GRICRA',
 746                               ITER_INTE_MAXI=500,
 747                               GROUP_MA='ELA',),
 748                            _F(RELATION='ELAS',
 749                               GROUP_MA=('EBOINF','EBOSUP',
 750                                         'RIG','PQSUP'),),),
 751                     INCREMENT=_F(
 752                            LIST_INST=L_INST2,
 753                            SUBD_PAS=10,
 754                            SUBD_PAS_MINI=1.E-6,
 755                            COEF_SUBD_PAS_1=0.0001,
 756                            ),
 757                     ETAT_INIT=_F(EVOL_NOLI=RESU,
 758                                  INST=0.,),
 759                     NEWTON=_F(MATRICE='TANGENTE',
 760                               REAC_ITER=1,),
 761                     CONVERGENCE=_F(RESI_GLOB_RELA=1.e-06,
 762                                   ITER_GLOB_MAXI=50,),
 763                     ARCHIVAGE=_F(LIST_INST=L_INST2,),);
 764
 765
 766 # POST_TRAITEMENT
 767
 768 COMPR=CALC_NO(reuse =COMPR,
 769   RESULTAT=COMPR,
 770   MODELE=MO1,
 771   TOUT_ORDRE='OUI',
 772   OPTION='FORC_NODA',
 773   CHAM_MATER=A_MAT,
 774   CARA_ELEM=CARA1,
 775          );
 776
 777 # Deplacement de l'embout
 778 DX1=RECU_FONCTION(
 779    RESULTAT=COMPR,
 780    NOM_CHAM='DEPL',
 781    NOM_CMP='DX',
 782    GROUP_NO='PEBO_S',
 783    PROL_DROITE='CONSTANT',
 784    PROL_GAUCHE='CONSTANT',);
 785 DXTOTA=CALC_FONCTION(
 786                      COMB=(_F(FONCTION=DX1,COEF=-1.,),),
 787                      );
 788
 789 # Force exercee sur l'embout
 790 POST3=POST_RELEVE_T(
 791    ACTION=_F(INTITULE='FORCES',
 792              GROUP_NO=('PEBO_S'),
 793              RESULTAT=COMPR,
 794              NOM_CHAM='FORC_NODA',
 795              RESULTANTE=('DX','DY','DZ'),
 796              REPERE='GLOBAL',
 797              OPERATION='EXTRACTION',),);
 798 FORC3=RECU_FONCTION(
 799              TABLE=POST3,
 800              PARA_X='INST',
 801              PARA_Y='DX',
 802              PROL_DROITE='CONSTANT',
 803              PROL_GAUCHE='CONSTANT',);
 804
 805 FORXA=CALC_FONCTION(
 806    COMB=_F(FONCTION=FORC3, COEF=-1.,),
 807    );
 808
 809
 810 titi=DXTOTA.Absc();
 811 toto=DXTOTA.Ordo();
 812 tata=FORXA.Ordo();
 813 LINS=titi[0:];
 814 DEPLA=toto[0:];
 815 FORCE=tata[0:];
 816
 817 ini=DEPLA[0];
 818 DEPL=[];
 819
 820 for i in range(len(DEPLA)):
 821   DEPL.append(DEPLA[i]-ini)
 822
 823
 824 # IMPR_FONCTION(
 825 #    FORMAT='XMGRACE',PILOTE='INTERACTIF',
 826 #    COURBE=(
 827 #           _F(ABSCISSE=DEPL,
 828 #              ORDONNEE=FORCE,
 829 #              COULEUR=2,
 830 #              ),
 831 #            _F(FONCTION=COMP_CH,COULEUR=3,),
 832 #            _F(FONCTION=COMP_DCH,COULEUR=3,),
 833 #    ),
 834 #    LEGENDE_X='Deplacement',
 835 #    LEGENDE_Y='Force',
 836 # )
 837
 838 FODEP=CREA_TABLE(LISTE=(_F(LISTE_R=tuple(DEPL),PARA='X'),
 839                         _F(LISTE_R=tuple(FORCE),PARA='Y'),
 840                         _F(LISTE_R=tuple(LINS),PARA='K'),
 841                         ),
 842                 )
 843
 844 IMPR_TABLE(TABLE=FODEP)
 845
 846 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 2.00000E+00),),
 847            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 848            PRECISION=0.05,VALE=1.0331E+04)
 849 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 2.00000E+00),),
 850            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 851            PRECISION=0.06,VALE=2.5E-04)
 852 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 3.00000E+00),),
 853            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 854            PRECISION=0.05,VALE=2.5E+04)
 855 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 3.00000E+00),),
 856            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 857            PRECISION=0.05,VALE=6.63E-04)
 858 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 4.00000E+00),),
 859            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 860            PRECISION=0.05,VALE=3.4702E+04)
 861 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 4.00000E+00),),
 862            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 863            PRECISION=0.05,VALE=1.E-03)
 864 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 5.00000E+00),),
 865            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 866            PRECISION=0.05,VALE=4.E+04)
 867 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 5.00000E+00),),
 868            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 869            PRECISION=0.05,VALE=1.23E-03)
 870 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 6.00000E+00),),
 871            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 872            PRECISION=0.05,VALE=3.5E+04)
 873 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 6.00000E+00),),
 874            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 875            PRECISION=0.05,VALE=1.09E-03)
 876 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 7.00000E+00),),
 877            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 878            PRECISION=0.05,VALE=2.5E+04)
 879 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 7.00000E+00),),
 880            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 881            PRECISION=0.05,VALE=8.49E-04)
 882 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 8.00000E+00),),
 883            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 884            PRECISION=0.05,VALE=1.0331E+04)
 885 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 8.00000E+00),),
 886            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 887            PRECISION=0.12,VALE=5.18E-04)
 888 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 9.00000E+00),),
 889            NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 890            PRECISION=0.05,VALE=5000.)
 891 TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 9.00000E+00),),
 892            NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 893            PRECISION=0.18,VALE=3.87E-04)
 894
 895
 896 FIN();