src/TestsEltsVisu/Nicolas/ssnl132ab.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 #
 259
 260 M1= I_CR*2 ;
 261 M2= I_TG_C*2 ;
 262 M3= I_TG_R*2 ;
 263 M4= I_TG_B*2;
 264 CARA1=AFFE_CARA_ELEM(
 265    MODELE=MO1,
 266    POUTRE=(
 267      _F(GROUP_MA='CRAYON',
 268         SECTION='GENERALE',
 269         CARA=('A','IZ','IY','JX','AY','AZ','EY','EZ'),
 270         VALE=(S_CR,I_CR,I_CR,M1,1,1,0,0),),
 271      _F(GROUP_MA='LGUIDE',
 272         SECTION='GENERALE',
 273         CARA=('A','IZ','IY','JX','AY','AZ'),
 274         VALE=(S_TG_C,I_TG_C,I_TG_C,M2,1,1),),
 275      _F(GROUP_MA='RETRE',
 276         SECTION='GENERALE',
 277         CARA=('A','IZ','IY','JX','AY','AZ'),
 278    VALE=(S_TG_R,I_TG_R,I_TG_R,M3,1,1),),
 279 _F(GROUP_MA='BIAIS',
 280         SECTION='GENERALE',
 281         CARA=('A','IZ','IY','JX','AY','AZ'),
 282         VALE=(S_TG_B,I_TG_B,I_TG_B,M4,1,1),),
 283      _F(GROUP_MA=('MAINTI_Y','MAINTI_Z'),
 284         SECTION='RECTANGLE',
 285         CARA=('HY','HZ'),
 286         VALE=(HYMA,HZMA),),
 287      _F(GROUP_MA=('EBOINF'),
 288         SECTION='RECTANGLE',
 289         CARA='H',
 290         VALE= Heinf,),
 291      _F(GROUP_MA=('EBOSUP'),
 292         SECTION='RECTANGLE',
 293         CARA='H',
 294         VALE= Hesup,),
 295            ),
 296    COQUE=(
 297       _F(GROUP_MA=('PQSUP'),
 298          COQUE_NCOU=1,
 299          EPAIS=1.,
 300          ANGL_REP=(90.,0.,),
 301           ),
 302            ),
 303    AFFE_FIBRE =(
 304           _F(GROUP_MA='CRAYON',
 305           COOR_AXE_POUTRE = (0.,0.,),
 306           CARA='SURFACE',
 307           VALE=(
 308              0.,(RAYCRA-EPCRA/2),S_CR/4,
 309              0.,-(RAYCRA-EPCRA/2),S_CR/4,
 310              (RAYCRA-EPCRA/2),0.,S_CR/4,
 311             -(RAYCRA-EPCRA/2),0.,S_CR/4,
 312                 )
 313                ),
 314           _F(GROUP_MA='LGUIDE',
 315                COOR_AXE_POUTRE = (0.,0.,),
 316                CARA='SURFACE',
 317                VALE=(
 318                 0.,(RAY1GU-EP1GU/2),S_TG_C/4,
 319                 0.,-(RAY1GU-EP1GU/2),S_TG_C/4,
 320                 (RAY1GU-EP1GU/2),0.,S_TG_C/4,
 321                -(RAY1GU-EP1GU/2),0.,S_TG_C/4,
 322                     )
 323                ),
 324           _F(GROUP_MA='BIAIS',
 325                COOR_AXE_POUTRE = (0.,0.,),
 326                 CARA='SURFACE',
 327                 VALE=(
 328                 0.,(RAY2GU-EPMOY/2),S_TG_B/4,
 329                 0.,-(RAY2GU-EPMOY/2),S_TG_B/4,
 330                 (RAY2GU-EPMOY/2),0.,S_TG_B/4,
 331                -(RAY2GU-EPMOY/2),0.,S_TG_B/4,
 332                     )
 333                ),
 334           _F(GROUP_MA='RETRE',
 335                COOR_AXE_POUTRE = (0.,0.,),
 336                 CARA='SURFACE',
 337                 VALE=(
 338                 0.,(RAY2GU-EP2GU/2),S_TG_R/4,
 339                 0.,-(RAY2GU-EP2GU/2),S_TG_R/4,
 340                 (RAY2GU-EP2GU/2),0.,S_TG_R/4,
 341                -(RAY2GU-EP2GU/2),0.,S_TG_R/4,
 342                     )
 343                ),
 344                     ),
 345      DISCRET=(
 346                _F(GROUP_MA='RIG_ME',
 347                     REPERE='LOCAL',
 348                     CARA='K_TR_D_L',
 349                     VALE=((1.E9*NBTG/4),(1.E9*NBTG/4),0.,(KR_GM*NBTG/8.),(1.E9*NBTG/4),(KR_GM*NBTG/8.),),),
 350                _F(GROUP_MA='RIG_EX',
 351                     REPERE='LOCAL',
 352                     CARA='K_TR_D_L',
 353                     VALE=((1.E9*NBTG/4),(1.E9*NBTG/4),0.,(KR_GM*NBTG/8.),(1.E9*NBTG/4),(KR_GM*NBTG/8.),),),),
 354
 355      ORIENTATION=(
 356                _F(GROUP_MA=('ELA_EX','ELA_ME','RIG_EX','RIG_ME'),
 357                   CARA='VECT_Y',
 358                   VALE=(1.,0.,0.),),
 359                _F(GROUP_MA=('MAINTI_Y'),
 360                   CARA='VECT_Y',
 361                   VALE=(0. ,1. ,0.),),
 362                _F(GROUP_MA=('MAINTI_Z'),
 363                   CARA='VECT_Y',
 364                   VALE=(0. ,0. ,1. ),),
 365                     ),
 366                );
 367
 368 MA1=MODI_MAILLAGE(
 369    reuse =MA1,
 370    MAILLAGE=MA1,
 371    MODELE=MO1,
 372    ORIE_NORM_COQUE=
 373       _F(   GROUP_MA='PQSUP',
 374             VECT_NORM=(-1.,0.,0.),
 375             GROUP_NO='PQ1',),
 376   INFO=1,);
 377
 378 #####################################################
 379 #              DEFINITION DES MATERIAUX             #
 380 #####################################################
 381 #
 382 #     EMBOUT SUPERIEUR
 383 MT_ES=DEFI_MATERIAU(ELAS=_F(E=EES0,
 384                                NU=NUES0,
 385                                RHO=RHOES,
 386                     ),);
 387
 388 #   EMBOUT INFERIEUR
 389 MT_EI=DEFI_MATERIAU(ELAS=_F(E=EEI0,
 390                                NU=NUEI0,
 391                                RHO=RHOEI,
 392                     ),);
 393 MT_RIG=DEFI_MATERIAU(ELAS=_F(E=1.0,
 394                      NU=0.0,
 395                      RHO=0.0,
 396                      ALPHA=0.0,),);
 397
 398 #    LIAISON GRILLE/CRAYON
 399 #kn  =37000.*NBCR;
 400 #kt  =140000.*NBCR  ;
 401 #rn  = 3700.*NBCR;
 402 #rt  =70000.*NBCR  ;
 403 #lb  =0.04;
 404 #fno  =26.*NBCR/1.25   ;
 405 #coul_bos  =0.51  ;
 406 #coul_res  =coul_bos*rt/2./kt  ;
 407 #ecrob  =kt/1e9 ;
 408 #ecror  =rt/1e9  ;
 409 #
 410 #M_GM=DEFI_MATERIAU(DIS_GRICRA=_F(KN_BOS=kn,
 411 #                                 KT_BOS=kt,
 412 #                                 KN_RES=rn,
 413 #                                 KT_RES= rt,
 414 #                                 DIST_BOS= lb,
 415 #                                 FORC_SER= fno,
 416 #                                 COUL_BOS=coul_bos,
 417 #                                 COUL_RES=coul_res,
 418 #                                 ECRO_BOS =ecrob,
 419 #                                 ECRO_RES =ecror,
 420 #                                   ),
 421 #                   ELAS=_F( E = 1.E5, NU = 0.3,),
 422 #                                    );
 423 #
 424 ##    SYSTEME DE MAINTIEN
 425 #MT_MA=DEFI_MATERIAU(ELAS=_F(E=EMA0,
 426 #                               NU=NUMA0,
 427 #                               RHO=RHOMA,),
 428 #                    ECRO_LINE=_F(D_SIGM_EPSI=DSMA0,
 429 #                                    SY=SYMA0,),
 430 #                    VMIS_POUTRE=_F(NP=NPMA0,
 431 #                                      MEY=MEYMA0,
 432 #                                      MPY=MPYMA0,
 433 #                                      CAY=CAYMA0,
 434 #                                      CBY=CBYMA0,
 435 #                                      MEZ=MEZMA0,
 436 #                                      MPZ=MPZMA0,
 437 #                                      CAZ=CAZMA0,
 438 #                                      CBZ=CBZMA0,
 439 #                                      MPX=MPXMA0,),);
 440 #
 441 ##     TUBE GUIDE
 442 #MT_TU1=DEFI_MATERIAU(
 443 #   ELAS=_F(
 444 #      E     =  (ETU0*NBTG),
 445 #      NU    =  NUTU0,
 446 #      RHO   =  (RHOTU*NBTG),),
 447 #      );
 448 #
 449 ##     CRAYON
 450 #MT_CR=DEFI_MATERIAU(
 451 #   ELAS=_F(
 452 #      E     =  (ECR0*NBCR),
 453 #      NU    =  NUCR0,
 454 #      RHO   =  (RHOCR*NBCR),
 455 #       ),
 456 #      );
 457 #
 458 ##     PQSUP
 459 #MT_PLAQ=DEFI_MATERIAU(
 460 #   ELAS=_F(E=1.E14,
 461 #           NU=0.3,
 462 #           RHO=0.,
 463 #           ALPHA=0.,),);
 464 #
 465 ###############################################################
 466 ##           AFFECTATION DES MATERIAUX                        #
 467 ###############################################################
 468 #A_MAT=AFFE_MATERIAU(MAILLAGE=MA1,
 469 #                    AFFE=(_F(GROUP_MA=('PQSUP'),
 470 #                             MATER=MT_PLAQ,),
 471 #                          _F(GROUP_MA=('MAINTIEN'),
 472 #                             MATER=MT_MA,),
 473 #                          _F(GROUP_MA='CRAYON',
 474 #                             MATER=MT_CR,),
 475 #                          _F(GROUP_MA='T_GUIDE',
 476 #                             MATER=MT_TU1,),
 477 #                          _F(GROUP_MA='EBOINF',
 478 #                             MATER=MT_EI,),
 479 #                          _F(GROUP_MA=('EBOSUP'),
 480 #                             MATER=MT_ES,),
 481 #                          _F(GROUP_MA='RIG',
 482 #                             MATER=MT_RIG,),
 483 #                          _F(GROUP_MA='ELA_EX_B',
 484 #                             MATER=M_GM,),
 485 #                          _F(GROUP_MA='ELA_EX_T',
 486 #                             MATER=M_GM,),
 487 #                          _F(GROUP_MA='ELA_ME',
 488 #                             MATER=M_GM,),
 489 #                            ),
 490 #                            );
 491 #
 492 #####################################################
 493 ##   DEFINITION DES CHARGEMENTS ET DES COND. LIM.   #
 494 #####################################################
 495 #
 496 ## CONDITIONS LIMITES
 497 ## ---------------------
 498 #CL_CTC=AFFE_CHAR_MECA(
 499 #   MODELE=MO1,
 500 #   VERI_NORM='NON',
 501 #   CONTACT=(_F(APPARIEMENT='NODAL',
 502 #              VECT_NORM_ESCL=(1.,0.,0.),
 503 #              GROUP_MA_ESCL='PINF',
 504 #              GROUP_MA_MAIT='CR_BAS',)
 505 #              ,
 506 #           _F(APPARIEMENT='NODAL',
 507 #              VECT_NORM_ESCL=(-1.,0.,0.),
 508 #              GROUP_MA_ESCL='PSUP',
 509 #              GROUP_MA_MAIT='CR_HAU',),),
 510 #   );
 511 #
 512 #
 513 #CL_PERM=AFFE_CHAR_MECA(MODELE=MO1,
 514 #                       DDL_IMPO=(_F(GROUP_MA='PQSUP',
 515 #                                    DX=0.,   DY=0.,   DZ=0.,),
 516 #                                 _F(GROUP_MA='ASSEM2',
 517 #                                    DY=0.,),
 518 #                                 _F(GROUP_NO='FIX',
 519 #                                    DX=0.,
 520 #                                    DZ=0.,
 521 #                                    DRX=0.,
 522 #                                      ),
 523 #                                 _F(GROUP_NO='CRAYON',
 524 #                                    DRX=0.,
 525 #                                      ),
 526 #                                 _F(GROUP_MA=('EBOSUP','EBOINF'),
 527 #                                    DZ=0.0,
 528 #                                    DRX=0.0,
 529 #                                    ),
 530 #                                 _F(GROUP_NO=('PGRIL_1','PGRIL_2','PGRIL_3','PGRIL_4',
 531 #                                           'PGRIL_5','PGRIL_6','PGRIL_7','PGRIL_8',
 532 #                                           'PGRIL_9','PGRIL_10',),
 533 #                                    DY=0.,DRX=0.,DRY=0.,DRZ=0.),
 534 #                                    ),
 535 #                      LIAISON_UNIF=(
 536 #                                    _F(GROUP_MA='GRIL_1',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 537 #                                    _F(GROUP_MA='GRIL_2',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 538 #                                    _F(GROUP_MA='GRIL_3',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 539 #                                    _F(GROUP_MA='GRIL_4',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 540 #                                    _F(GROUP_MA='GRIL_5',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 541 #                                    _F(GROUP_MA='GRIL_6',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 542 #                                    _F(GROUP_MA='GRIL_7',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 543 #                                    _F(GROUP_MA='GRIL_8',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 544 #                                    _F(GROUP_MA='GRIL_9',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 545 #                                    _F(GROUP_MA='GRIL_10',DDL=('DX','DY','DZ','DRX','DRY','DRZ',),),
 546 #                                    ),
 547 #                                    );
 548 #
 549 #
 550 #BLOQ=AFFE_CHAR_MECA(MODELE=MO1,
 551 #                    DDL_IMPO=_F(GROUP_NO=('PGRIL_1','PGRIL_2','PGRIL_3','PGRIL_4',
 552 #                                          'PGRIL_5','PGRIL_6','PGRIL_7','PGRIL_8',
 553 #                                          'PGRIL_9','PGRIL_10',),
 554 #                                DX=0.,),
 555 #                    );
 556 #
 557 #
 558 ##
 559 ##      CHARGEMENT
 560 ## ---------------------
 561 ##PESANT=AFFE_CHAR_MECA(MODELE=MO1,
 562 ##                      PESANTEUR=(9.81,-1.,0.,0.),);
 563 ##
 564 ##
 565 ### Compression initiale
 566 ##TIR_CR=AFFE_CHAR_MECA(
 567 ##   MODELE=MO1,
 568 ##   DDL_IMPO=
 569 ##   _F(GROUP_NO='CR_BAS',DX=-1.E-3,),);
 570 ##DEPL_CR=DEFI_FONCTION(
 571 ##   NOM_PARA='INST',
 572 ##   VALE=(
 573 ##     -2.,  0. ,
 574 ##     -1.,  1. ,),);
 575 ##
 576 ##
 577 ##CL_DI_FR=AFFE_CHAR_MECA(
 578 ##      MODELE=MO1,
 579 ##      FORCE_NODALE=
 580 ##         _F(GROUP_NO='PEBO_S',
 581 ##            FX=(-1.),),);
 582 ##
 583 ##CL_GRI_6=AFFE_CHAR_MECA(
 584 ##      MODELE=MO1,
 585 ##      DDL_IMPO=
 586 ##         _F(GROUP_NO='PGRIL_6',
 587 ##            DZ=1.E-3,),);
 588 ##
 589 ##FOR_MAIN=DEFI_FONCTION(NOM_PARA='INST',
 590 ##                       VALE=(
 591 ##                  0.,   0. ,
 592 ##                  1.,   1.e-10 ,
 593 ##                  2.,   1.0331E+04,
 594 ##                  3.,   2.5E+04,
 595 ##                  4.,   3.4702E+04,
 596 ##                  5.,   4.E+04,
 597 ##                  6.,   3.5E+04,
 598 ##                  7.,   2.5E+04,
 599 ##                  8.,   1.0331E+04,
 600 ##                  9.,   5000.,
 601 ##                  10.,   0.
 602 ##                  ),
 603 ##                  PROL_DROITE='CONSTANT',
 604 ##                  PROL_GAUCHE='CONSTANT',);
 605 ##
 606 ##
 607 ##
 608 #################################################
 609 ###  DEFINITION DES LISTES D INSTANTS           #
 610 #################################################
 611 ##
 612 ##L_INST1=DEFI_LIST_REEL(
 613 ##   DEBUT=-2.,
 614 ##   INTERVALLE=(
 615 ##      _F(JUSQU_A=-1.,   NOMBRE=10,),
 616 ##               ),);
 617 ##
 618 ##L_INST11=DEFI_LIST_REEL(
 619 ##   DEBUT=-1.,
 620 ##   INTERVALLE=(
 621 ##      _F(JUSQU_A=0.,   NOMBRE=10,),
 622 ##               ),);
 623 ##
 624 ##L_INST2=DEFI_LIST_REEL(DEBUT=0.0,
 625 ##                      INTERVALLE=(_F(JUSQU_A=1.0,
 626 ##                                     NOMBRE=2,),
 627 ##                                  _F(JUSQU_A=10.0,
 628 ##                                     NOMBRE=45,),),);
 629 ##
 630 ##L_INST=DEFI_LIST_REEL(DEBUT=-2.0,
 631 ##                      INTERVALLE=(_F(JUSQU_A=-1.,
 632 ##                                       NOMBRE=10,),
 633 ##                                  _F(JUSQU_A=0.,
 634 ##                                       NOMBRE=10,),
 635 ##                                  _F(JUSQU_A=1.0,
 636 ##                                       NOMBRE=2,),
 637 ##                                  _F(JUSQU_A=10.0,
 638 ##                                       NOMBRE=45,),
 639 ##                                  ),);
 640 ##
 641 #######################################################################
 642 #######################################################################
 643 ###     CALCUL INITIAL COMMUN A L ESSAI DE COMPRESSION ET DE FLEXION  #
 644 #######################################################################
 645 #######################################################################
 646 ##
 647 ##RESU=STAT_NON_LINE(
 648 ##                   MODELE=MO1,
 649 ##                   CHAM_MATER=A_MAT,
 650 ##                   CARA_ELEM=CARA1,
 651 ##                   EXCIT=(_F(CHARGE=TIR_CR, FONC_MULT=DEPL_CR,),
 652 ##                          _F(CHARGE=BLOQ,),
 653 ##                          _F(CHARGE=CL_PERM,),
 654 ##                          _F(CHARGE=CL_CTC,),
 655 ##                          _F(CHARGE=PESANT,),
 656 ##                          ),
 657 ##                   COMP_INCR=(
 658 ##                          _F(RELATION='VMIS_POU_LINE',
 659 ##                             GROUP_MA='MAINTIEN',
 660 ##                             ITER_INTE_MAXI=500),
 661 ##                          _F(RELATION='ELAS',
 662 ##                             GROUP_MA=('CRAYON','T_GUIDE'),),
 663 ##                          _F(RELATION='DIS_GRICRA',
 664 ##                             ITER_INTE_MAXI=500,
 665 ##                             GROUP_MA='ELA'),
 666 ##                          _F(RELATION='ELAS',
 667 ##                             GROUP_MA=('EBOINF','EBOSUP',
 668 ##                                       'RIG','PQSUP'),),),
 669 ##                   INCREMENT=_F(
 670 ##                                 LIST_INST=L_INST1,
 671 ##                                 SUBD_PAS=10,
 672 ##                                 SUBD_PAS_MINI=1.E-6,
 673 ##                                 COEF_SUBD_PAS_1=0.0001,
 674 ##                                 ),
 675 ##                   NEWTON=_F(MATRICE='TANGENTE',
 676 ##                             REAC_ITER=1,),
 677 ##                   CONVERGENCE=_F(RESI_GLOB_RELA=1.e-06,
 678 ##                                  ITER_GLOB_MAXI=50,),
 679 ##                   ARCHIVAGE=_F(LIST_INST=L_INST1,),);
 680 ##
 681 ##RESU=STAT_NON_LINE(
 682 ##                   reuse = RESU,
 683 ##                   MODELE=MO1,
 684 ##                   CHAM_MATER=A_MAT,
 685 ##                   CARA_ELEM=CARA1,
 686 ##                   EXCIT=(_F(CHARGE=CL_PERM,),
 687 ##                          _F(CHARGE=CL_CTC,),
 688 ##                          _F(CHARGE=PESANT,),
 689 ##                          ),
 690 ##                   COMP_INCR=(
 691 ##                          _F(RELATION='VMIS_POU_LINE',
 692 ##                             GROUP_MA='MAINTIEN',
 693 ##                             ITER_INTE_MAXI=500,),
 694 ##                          _F(RELATION='ELAS',
 695 ##                             GROUP_MA=('CRAYON','T_GUIDE'),),
 696 ##                          _F(RELATION='DIS_GRICRA',
 697 ##                             GROUP_MA='ELA',
 698 ##                             ITER_INTE_MAXI=500,),
 699 ##                          _F(RELATION='ELAS',
 700 ##                             GROUP_MA=('EBOINF','EBOSUP',
 701 ##                                       'RIG','PQSUP'),),),
 702 ##                   INCREMENT=_F(
 703 ##                                LIST_INST=L_INST11,
 704 ##                                SUBD_PAS=10,
 705 ##                                SUBD_PAS_MINI=1.E-6,
 706 ##                                COEF_SUBD_PAS_1=0.0001,),
 707 ##                   ETAT_INIT=_F(
 708 ##                                EVOL_NOLI=RESU,
 709 ##                                INST=-1.,),
 710 ##                   NEWTON=_F(MATRICE='TANGENTE',
 711 ##                                REAC_ITER=1,),
 712 ##                   CONVERGENCE=_F(RESI_GLOB_RELA=1.e-06,
 713 ##                                  ITER_GLOB_MAXI=50,),
 714 ##                   ARCHIVAGE=_F(LIST_INST=L_INST11,),);
 715 ##
 716 ##
 717 #######################################################################
 718 #######################################################################
 719 ###     POURSUITE ESSAI DE COMPRESSION                                #
 720 #######################################################################
 721 #######################################################################
 722 ##
 723 ### CALCUL
 724 ##
 725 ##COMPR=STAT_NON_LINE(
 726 ##                    MODELE=MO1,
 727 ##                    CHAM_MATER=A_MAT,
 728 ##                    CARA_ELEM=CARA1,
 729 ##                    EXCIT=(_F(CHARGE=CL_DI_FR,FONC_MULT=FOR_MAIN,),
 730 ##                           _F(CHARGE=CL_PERM,),
 731 ##                           _F(CHARGE=CL_CTC,),
 732 ##                           _F(CHARGE=PESANT,),
 733 ##                           ),
 734 ##                    COMP_INCR=(
 735 ##                           _F(RELATION='VMIS_POU_LINE',
 736 ##                              GROUP_MA='MAINTIEN',
 737 ##                              ITER_INTE_MAXI=500),
 738 ##                           _F(RELATION='ELAS',
 739 ##                              GROUP_MA=('CRAYON','T_GUIDE'),),
 740 ##                           _F(RELATION='DIS_GRICRA',
 741 ##                              ITER_INTE_MAXI=500,
 742 ##                              GROUP_MA='ELA',),
 743 ##                           _F(RELATION='ELAS',
 744 ##                              GROUP_MA=('EBOINF','EBOSUP',
 745 ##                                        'RIG','PQSUP'),),),
 746 ##                    INCREMENT=_F(
 747 ##                           LIST_INST=L_INST2,
 748 ##                           SUBD_PAS=10,
 749 ##                           SUBD_PAS_MINI=1.E-6,
 750 ##                           COEF_SUBD_PAS_1=0.0001,
 751 ##                           ),
 752 ##                    ETAT_INIT=_F(EVOL_NOLI=RESU,
 753 ##                                 INST=0.,),
 754 ##                    NEWTON=_F(MATRICE='TANGENTE',
 755 ##                              REAC_ITER=1,),
 756 ##                    CONVERGENCE=_F(RESI_GLOB_RELA=1.e-06,
 757 ##                                  ITER_GLOB_MAXI=50,),
 758 ##                    ARCHIVAGE=_F(LIST_INST=L_INST2,),);
 759 ##
 760 ##
 761 ### POST_TRAITEMENT
 762 ##
 763 ##COMPR=CALC_NO(reuse =COMPR,
 764 ##  RESULTAT=COMPR,
 765 ##  MODELE=MO1,
 766 ##  TOUT_ORDRE='OUI',
 767 ##  OPTION='FORC_NODA',
 768 ##  CHAM_MATER=A_MAT,
 769 ##  CARA_ELEM=CARA1,
 770 ##         );
 771 ##
 772 ### Deplacement de l'embout
 773 ##DX1=RECU_FONCTION(
 774 ##   RESULTAT=COMPR,
 775 ##   NOM_CHAM='DEPL',
 776 ##   NOM_CMP='DX',
 777 ##   GROUP_NO='PEBO_S',
 778 ##   PROL_DROITE='CONSTANT',
 779 ##   PROL_GAUCHE='CONSTANT',);
 780 ##DXTOTA=CALC_FONCTION(
 781 ##                     COMB=(_F(FONCTION=DX1,COEF=-1.,),),
 782 ##                     );
 783 ##
 784 ### Force exercee sur l'embout
 785 ##POST3=POST_RELEVE_T(
 786 ##   ACTION=_F(INTITULE='FORCES',
 787 ##             GROUP_NO=('PEBO_S'),
 788 ##             RESULTAT=COMPR,
 789 ##             NOM_CHAM='FORC_NODA',
 790 ##             RESULTANTE=('DX','DY','DZ'),
 791 ##             REPERE='GLOBAL',
 792 ##             OPERATION='EXTRACTION',),);
 793 ##FORC3=RECU_FONCTION(
 794 ##             TABLE=POST3,
 795 ##             PARA_X='INST',
 796 ##             PARA_Y='DX',
 797 ##             PROL_DROITE='CONSTANT',
 798 ##             PROL_GAUCHE='CONSTANT',);
 799 ##
 800 ##FORXA=CALC_FONCTION(
 801 ##   COMB=_F(FONCTION=FORC3, COEF=-1.,),
 802 ##   );
 803 ##
 804 ##
 805 ##titi=DXTOTA.Absc();
 806 ##toto=DXTOTA.Ordo();
 807 ##tata=FORXA.Ordo();
 808 ##LINS=titi[0:];
 809 ##DEPLA=toto[0:];
 810 ##FORCE=tata[0:];
 811 ##
 812 ##ini=DEPLA[0];
 813 ##DEPL=[];
 814 ##
 815 ##for i in range(len(DEPLA)):
 816 ##  DEPL.append(DEPLA[i]-ini)
 817 ##
 818 ##
 819 ### IMPR_FONCTION(
 820 ###    FORMAT='XMGRACE',PILOTE='INTERACTIF',
 821 ###    COURBE=(
 822 ###           _F(ABSCISSE=DEPL,
 823 ###              ORDONNEE=FORCE,
 824 ###              COULEUR=2,
 825 ###              ),
 826 ###            _F(FONCTION=COMP_CH,COULEUR=3,),
 827 ###            _F(FONCTION=COMP_DCH,COULEUR=3,),
 828 ###    ),
 829 ###    LEGENDE_X='Deplacement',
 830 ###    LEGENDE_Y='Force',
 831 ### )
 832 ##
 833 ##FODEP=CREA_TABLE(LISTE=(_F(LISTE_R=tuple(DEPL),PARA='X'),
 834 ##                        _F(LISTE_R=tuple(FORCE),PARA='Y'),
 835 ##                        _F(LISTE_R=tuple(LINS),PARA='K'),
 836 ##                        ),
 837 ##                )
 838 ##
 839 ##IMPR_TABLE(TABLE=FODEP)
 840 ##
 841 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 2.00000E+00),),
 842 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 843 ##           PRECISION=0.05,VALE=1.0331E+04)
 844 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 2.00000E+00),),
 845 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 846 ##           PRECISION=0.06,VALE=2.5E-04)
 847 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 3.00000E+00),),
 848 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 849 ##           PRECISION=0.05,VALE=2.5E+04)
 850 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 3.00000E+00),),
 851 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 852 ##           PRECISION=0.05,VALE=6.63E-04)
 853 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 4.00000E+00),),
 854 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 855 ##           PRECISION=0.05,VALE=3.4702E+04)
 856 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 4.00000E+00),),
 857 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 858 ##           PRECISION=0.05,VALE=1.E-03)
 859 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 5.00000E+00),),
 860 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 861 ##           PRECISION=0.05,VALE=4.E+04)
 862 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 5.00000E+00),),
 863 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 864 ##           PRECISION=0.05,VALE=1.23E-03)
 865 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 6.00000E+00),),
 866 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 867 ##           PRECISION=0.05,VALE=3.5E+04)
 868 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 6.00000E+00),),
 869 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 870 ##           PRECISION=0.05,VALE=1.09E-03)
 871 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 7.00000E+00),),
 872 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 873 ##           PRECISION=0.05,VALE=2.5E+04)
 874 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 7.00000E+00),),
 875 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 876 ##           PRECISION=0.05,VALE=8.49E-04)
 877 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 8.00000E+00),),
 878 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 879 ##           PRECISION=0.05,VALE=1.0331E+04)
 880 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 8.00000E+00),),
 881 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 882 ##           PRECISION=0.12,VALE=5.18E-04)
 883 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 9.00000E+00),),
 884 ##           NOM_PARA='Y',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 885 ##           PRECISION=0.05,VALE=5000.)
 886 ##TEST_TABLE(TABLE=FODEP,FILTRE=(_F(NOM_PARA = 'K',VALE = 9.00000E+00),),
 887 ##           NOM_PARA='X',REFERENCE = 'SOURCE_EXTERNE',CRITERE='RELATIF',
 888 ##           PRECISION=0.18,VALE=3.87E-04)
 889 ##
 890 ##
 891 FIN();