src/TestsEltsVisu/Nicolas/forma12c.comm

   1 # MODIF  DATE 12/05/2005   AUTEUR DURAND C.DURAND
   2 # TITRE TP ANALYSE SISMIQUE
   3 #            CONFIGURATION MANAGEMENT OF EDF VERSION
   4 # ======================================================================
   5 # COPYRIGHT (C) 1991 - 2004  EDF R&D                  WWW.CODE-ASTER.ORG
   6 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
   7 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
   8 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
   9 # (AT YOUR OPTION) ANY LATER VERSION.
  10 #
  11 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  12 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  13 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  14 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  15 #
  16 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
  17 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
  18 #    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
  19 # ======================================================================
  20 # RESPONSABLE BOYERE E.BOYERE
  21 #-----------------------------------------------------------------------
  22 #
  23 #  ANALYSE RIS12 TR4 REGLEMENTAIRE
  24 #  SCENARIO 1
  25 #  GL1(X,Y) + GL2(Z)
  26 #
  27 #-----------------------------------------------------------------------
  28
  29 DEBUT(CODE=_F(NOM ='FORMA12C',NIV_PUB_WEB='INTERNET'))
  30
  31 MAILLAGE=LIRE_MAILLAGE(INFO=2)
  32
  33 MODELE=AFFE_MODELE(
  34                   MAILLAGE=MAILLAGE,
  35                    AFFE=(
  36                          _F(GROUP_MA='GPOU_D_T',
  37                             PHENOMENE='MECANIQUE',
  38                             MODELISATION='POU_D_T',),
  39                          _F(GROUP_MA='GPOU_C_T',
  40                             PHENOMENE='MECANIQUE',
  41                             MODELISATION='POU_C_T',),
  42                          _F(GROUP_MA='GDIS_TR',
  43                             PHENOMENE='MECANIQUE',
  44                             MODELISATION='DIS_TR',),
  45                             ),
  46                               )
  47
  48 IMPR_RESU(FORMAT='CASTEM',UNITE=37,
  49           MODELE=MODELE,
  50           RESU=_F(
  51                   MAILLAGE=MAILLAGE))
  52
  53
  54
  55 #-----------------------------------------------------------------------
  56 # CARACTERISTIQUES DES ELEMENTS
  57 #-----------------------------------------------------------------------
  58
  59
  60 CARA=AFFE_CARA_ELEM(
  61                     MODELE=MODELE,
  62                     POUTRE=(
  63                            _F(GROUP_MA='GMEL0301',
  64                                SECTION='CERCLE',
  65                                CARA=('R','EP'),
  66                                VALE=(0.08415,0.00711),),
  67                             _F(GROUP_MA='GMEL0501',
  68                                SECTION='CERCLE',
  69                                CARA=('R','EP'),
  70                                VALE=(0.08415,0.0034),),
  71                             _F(GROUP_MA='GMEL0601',
  72                                SECTION='CERCLE',
  73                                CARA=('R','EP'),
  74                                VALE=(0.08415,0.0034),),
  75                             _F(GROUP_MA='GMEL0701',
  76                                SECTION='CERCLE',
  77                                CARA=('R','EP'),
  78                                VALE=(0.08415,0.0034),),
  79                             _F(GROUP_MA='GMEL1001',
  80                                SECTION='CERCLE',
  81                                CARA=('R','EP'),
  82                                VALE=(0.08415,0.0034),),
  83                             _F(GROUP_MA='GMEL1101',
  84                                SECTION='CERCLE',
  85                                CARA=('R','EP'),
  86                                VALE=(0.08415,0.0034),),
  87                             _F(GROUP_MA='GMEL1201',
  88                                SECTION='CERCLE',
  89                                CARA=('R','EP'),
  90                                VALE=(0.08415,0.0034),),
  91                             _F(GROUP_MA='GMEL1501',
  92                                SECTION='CERCLE',
  93                                CARA=('R','EP'),
  94                                VALE=(0.08415,0.0034),),
  95                             _F(GROUP_MA='GMEL1701',
  96                                SECTION='CERCLE',
  97                                CARA=('R','EP'),
  98                                VALE=(0.08415,0.0034),),
  99                             _F(GROUP_MA='GMEL1801',
 100                                SECTION='CERCLE',
 101                                CARA=('R','EP'),
 102                                VALE=(0.08415,0.0034),),
 103                             _F(GROUP_MA='GMEL1901',
 104                                SECTION='CERCLE',
 105                                CARA=('R','EP'),
 106                                VALE=(0.08415,0.0034),),
 107                             _F(GROUP_MA='GMEL2001',
 108                                SECTION='CERCLE',
 109                                CARA=('R','EP'),
 110                                VALE=(0.08415,0.0034),),
 111                             _F(GROUP_MA='GMEL2101',
 112                                SECTION='CERCLE',
 113                                CARA=('R','EP'),
 114                                VALE=(0.08415,0.0034),),
 115                             _F(GROUP_MA='GMEL2201',
 116                                SECTION='CERCLE',
 117                                CARA=('R','EP'),
 118                                VALE=(0.08415,0.0034),),
 119                             _F(GROUP_MA='GMEL2401',
 120                                SECTION='CERCLE',
 121                                CARA=('R','EP'),
 122                                VALE=(0.08415,0.00711),),),
 123                     DISCRET=(_F(MAILLE='M0801I02',
 124                                 CARA='M_TR_D_N',
 125                                 VALE=(36.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0),),
 126                              _F(MAILLE='M1301I02',
 127                                 CARA='M_TR_D_N',
 128                                 VALE=(36.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0),),),
 129                     DEFI_ARC=(_F(GROUP_MA='GMEL0601',
 130                                  CENTRE=(466.39260860000002,107.5500031,0.27747499939999998),
 131                                  RAYON=0.22859999540000001,
 132                                  COEF_FLEX_XY=14.43145561,
 133                                  COEF_FLEX_XZ=14.43145561,),
 134                               _F(GROUP_MA='GMEL1101',
 135                                  CENTRE=(466.16400149999998,107.3214035,3.7724525930000001),
 136                                  RAYON=0.22859999540000001,
 137                                  COEF_FLEX_XY=14.43145561,
 138                                  COEF_FLEX_XZ=14.43145561,),
 139                               _F(GROUP_MA='GMEL1701',
 140                                  CENTRE=(466.16400149999998,106.4785995,4.2335472110000003),
 141                                  RAYON=0.22859999540000001,
 142                                  COEF_FLEX_XY=14.43145561,
 143                                  COEF_FLEX_XZ=14.43145561,),
 144                               _F(GROUP_MA='GMEL1901',
 145                                  CENTRE=(466.16400149999998,106.4785995,8.7655010220000005),
 146                                  RAYON=0.22859999540000001,
 147                                  COEF_FLEX_XY=14.43145561,
 148                                  COEF_FLEX_XZ=14.43145561,),
 149                               _F(GROUP_MA='GMEL2101',
 150                                  CENTRE=(466.16400149999998,107.4713974,9.2274990080000006),
 151                                  RAYON=0.22859999540000001,
 152                                  COEF_FLEX_XY=14.43145561,
 153                                  COEF_FLEX_XZ=14.43145561,),),)
 154
 155 #-----------------------------------------------------------------------
 156 # DEFINITION DES MATERIAUX
 157 #-----------------------------------------------------------------------
 158
 159 MATP01=DEFI_MATERIAU(ELAS=_F(E=1.97E11,
 160                              NU=0.3,
 161                              RHO=1.3108749019999999E4,
 162                              ALPHA=1.639999937E-05,),)
 163 MATP03=DEFI_MATERIAU(ELAS=_F(E=1.97E11,
 164                              NU=0.3,
 165                              RHO=1.956015625E4,
 166                              ALPHA=1.639999937E-05,),)
 167
 168 #-----------------------------------------------------------------------
 169 # DEFINITION DES SUPPORTS
 170 #-----------------------------------------------------------------------
 171
 172 SUPPORT=AFFE_CHAR_MECA(
 173       MODELE=MODELE,
 174       DDL_IMPO=(
 175 # PF1 ET PF2
 176      _F( NOEUD = ('N02__I00', 'N01__I00',),
 177                DX = 0., DY = 0., DZ = 0., DRX = 0., DRY = 0., DRZ = 0.),
 178 # GL1
 179      _F( NOEUD = 'N0701I00',
 180                DX = 0., DY = 0.),
 181 # GL2
 182      _F( NOEUD = 'N1201I00',
 183                DZ = 0.),
 184               ),
 185                   )
 186
 187 #-----------------------------------------------------------------------
 188 # DEFINITION DES CHAMPS DE MATERIAUX
 189 #-----------------------------------------------------------------------
 190
 191 CHMATP=AFFE_MATERIAU(MAILLAGE=MAILLAGE,
 192                      AFFE=(_F(GROUP_MA='GMAT01',
 193                               MATER=MATP01,
 194                               TEMP_REF=20.0,),
 195                            _F(GROUP_MA='GMAT03',
 196                               MATER=MATP03,
 197                               TEMP_REF=20.0,),),)
 198
 199 #-----------------------------------------------------------------------
 200 # ASSEMBLAGE MATRICE MASSE ET RIGIDITE
 201 #-----------------------------------------------------------------------
 202
 203 MACRO_MATR_ASSE(MODELE=MODELE,
 204                 CHAM_MATER=CHMATP,
 205                 CARA_ELEM=CARA,
 206                 CHARGE=SUPPORT,
 207                 NUME_DDL=CO('NUMDDL'),
 208                 MATR_ASSE=(_F(MATRICE=CO("MATRRIGI"),
 209                               OPTION='RIGI_MECA',),
 210                            _F(MATRICE=CO("MATRMASS"),
 211                               OPTION='MASS_MECA',),),);
 212
 213
 214 #-----------------------------------------------------------------------
 215 # ASSEMBLAGE DU SECOND MEMBRE
 216 #-----------------------------------------------------------------------
 217
 218 INCLUDE(UNITE=22)
 219
 220 CHSEIX=CALC_CHAR_SEISME(  MATR_MASS=MATRMASS,
 221      DIRECTION=(1., 0., 0.,), MONO_APPUI='OUI'
 222       )
 223 CHSEIY=CALC_CHAR_SEISME(  MATR_MASS=MATRMASS,
 224      DIRECTION=(0., 1., 0.,), MONO_APPUI='OUI'
 225       )
 226 CHSEIZ=CALC_CHAR_SEISME(  MATR_MASS=MATRMASS,
 227      DIRECTION=(0., 0., 1.,), MONO_APPUI='OUI'
 228       )
 229
 230 #-----------------------------------------------------------------------
 231 # CALCUL DES MODES
 232 #-----------------------------------------------------------------------
 233 MASSINER=POST_ELEM(
 234                MODELE=MODELE,
 235                CHARGE=SUPPORT,
 236                CARA_ELEM=CARA,
 237                CHAM_MATER=CHMATP,
 238                MASS_INER=_F( TOUT = 'OUI',)
 239                         )
 240 MODES=MACRO_MODE_MECA(
 241                MATR_A=MATRRIGI,
 242                MATR_B=MATRMASS,
 243                CALC_FREQ=_F(
 244                FREQ_MIN = 0.0,
 245                FREQ_MAX = 33.0,
 246                NB_BLOC_FREQ = 1),
 247                VERI_MODE=_F( STOP_ERREUR = 'NON'),
 248                NORM_MODE=_F(
 249                      NORME = 'MASS_GENE',
 250                      MASS_INER =MASSINER),
 251                  IMPRESSION=_F( TOUT_PARA = 'OUI')
 252                   )
 253
 254 # TEST_RESU UNIQUEMENT POUR FAIRE CAS TEST
 255 TEST_RESU(RESU=(_F(RESULTAT=MODES,
 256                    NUME_ORDRE=1,
 257                    PARA='FREQ',
 258                    VALE= 4.5135794612152,
 259                    PRECISION=1e-05,
 260                    REFERENCE='NON_REGRESSION',
 261                    VERSION='7.3.19',)))
 262
 263 MODES=CALC_ELEM(
 264                 reuse=MODES,
 265                 MODELE=MODELE,
 266                 CHAM_MATER=CHMATP,
 267                 CARA_ELEM=CARA,
 268                 OPTION=(
 269                 'EFGE_ELNO_DEPL',
 270                   ),
 271                 RESULTAT=MODES,)
 272
 273 ###########################################################
 274 #----  REPONSE TRANSITOIRE PAR SUPERPOSITION MODALE ------
 275 #----  MONO APPUI GRANDEURS RELATIVES ------
 276 ###########################################################
 277
 278 # DEFINITION DES INSTANTS DE RECUPERATION
 279
 280 L_RECU=DEFI_LIST_REEL(  DEBUT=0.,
 281                       INTERVALLE=_F(
 282                         JUSQU_A = 20.47,
 283                         PAS = 0.01)
 284                               )
 285
 286 #----------------------------------------------------
 287 #------------- PROJECTION SUR LA BASE MODALE --------
 288 #----------------------------------------------------
 289 MACRO_PROJ_BASE(
 290                BASE=MODES,
 291                MATR_ASSE_GENE=(
 292                         _F(
 293                         MATRICE=CO("MASSEGEN"),
 294                         MATR_ASSE=MATRMASS,
 295                         ),
 296                         _F(
 297                         MATRICE=CO("RIGIDGEN"),
 298                         MATR_ASSE=MATRRIGI,
 299                         ),
 300                ),
 301                VECT_ASSE_GENE=(
 302                         _F(
 303                         VECTEUR=CO("EFGENX"),
 304                         VECT_ASSE=CHSEIX,
 305                         ),
 306                         _F(
 307                         VECTEUR=CO("EFGENY"),
 308                         VECT_ASSE=CHSEIY,
 309                         ),
 310                         _F(
 311                         VECTEUR=CO("EFGENZ"),
 312                         VECT_ASSE=CHSEIZ,
 313                         ),
 314                       ),
 315                      )
 316
 317 #------------------------------------------------------------------
 318 #-------------- RESOLUTION DU PB GENERALISE -----------------------
 319 #-------------- SANS CORRECTION STATIQUE --------------------------
 320 #------------------------------------------------------------------
 321 GSCR=DYNA_TRAN_MODAL(
 322                   MASS_GENE=MASSEGEN,
 323                   RIGI_GENE=RIGIDGEN,
 324                   AMOR_REDUIT=0.02,
 325                   EXCIT=(
 326                   _F(
 327                      VECT_GENE=EFGENX,
 328                      FONC_MULT=ACCELH1,
 329                         ),
 330                   _F(
 331                      VECT_GENE=EFGENY,
 332                      FONC_MULT=ACCELH1,
 333                         ),
 334                   _F(
 335                      VECT_GENE=EFGENZ,
 336                      FONC_MULT=ACCELV1,
 337                         ),
 338                            ),
 339                   INCREMENT=_F(
 340                      INST_INIT = 0.,
 341                      INST_FIN =20.47,
 342                      PAS = 0.001
 343                            ),
 344                         )
 345
 346 #--------------------------------------------------------------
 347 #---- RESTITUTION GRANDEUR PHYSIQUE DANS LE REPERE RELATIF
 348 #--- SANS CORRECTION STATIQUE
 349 #--------------------------------------------------------------
 350 # RESTITUTION GLOBALE DU CHAMP DE DEPLACEMENT RELATIF A CHAQUE INSTANT
 351 # DE LA LISTE L_RECU
 352 RGSCR=REST_BASE_PHYS(
 353                RESU_GENE=GSCR,
 354                INTERPOL='LIN',
 355                LIST_INST=L_RECU,
 356                #TOUT_INST='OUI',
 357                TOUT_CHAM='OUI',
 358                       )
 359 DPSCR=CREA_CHAMP(
 360               TYPE_CHAM='NOEU_DEPL_R',
 361               OPERATION='EXTR',
 362               RESULTAT=RGSCR,
 363               NOM_CHAM='DEPL',
 364               TYPE_MAXI='MAXI',
 365               TYPE_RESU='VALE',
 366               TOUT_ORDRE='OUI'
 367                      )
 368 TDPSCR=POST_RELEVE_T(
 369               ACTION=
 370               _F(
 371             INTITULE = 'DEPL_MAX_R',
 372             GROUP_NO='GNSTR01',
 373             CHAM_GD=DPSCR,
 374             TOUT_CMP='OUI',
 375             OPERATION ='EXTRACTION',
 376                      ),
 377             )
 378
 379
 380 # RESTITUTION GLOBALE DU CHAMP EFGE_ELNO_DEPL A CHAQUE INSTANT
 381 # DE LA LISTE L_RECU
 382 RGSCR=CALC_ELEM(
 383                 reuse=RGSCR,
 384                 MODELE=MODELE,
 385                 CHAM_MATER=CHMATP,
 386                 CARA_ELEM=CARA,
 387                 OPTION=(
 388                 'EFGE_ELNO_DEPL',
 389                   ),
 390                 RESULTAT=RGSCR,)
 391 EFSCR=CREA_CHAMP(
 392               TYPE_CHAM='ELNO_SIEF_R',
 393               OPERATION='EXTR',
 394               RESULTAT=RGSCR,
 395               NOM_CHAM='EFGE_ELNO_DEPL',
 396               TYPE_MAXI='MAXI',
 397               TYPE_RESU='VALE',
 398               TOUT_ORDRE='OUI'
 399                      )
 400 TEFSCR=POST_RELEVE_T(
 401               ACTION=
 402               _F(
 403             INTITULE = 'EFGE_MAX_R',
 404             GROUP_NO='GNSTR01',
 405             CHAM_GD=EFSCR,
 406             TOUT_CMP='OUI',
 407             OPERATION ='EXTRACTION',
 408                      ),
 409             )
 410 IMPR_TABLE(
 411           TABLE=TEFSCR,
 412          FORMAT='TABLEAU',
 413                )
 414
 415 ###########################################################
 416 #----  REPONSE TRANSITOIRE PAR SUPERPOSITION MODALE ------
 417 #----  MULTI APPUIS GRANDEURS ABSOLUES ------
 418 ###########################################################
 419
 420 # DEFINITION EXCITATION MULTIPLE
 421
 422 MSTA=MODE_STATIQUE(
 423          MATR_RIGI=MATRRIGI,
 424          MATR_MASS=MATRMASS,
 425          MODE_STAT=_F(TOUT='OUI',AVEC_CMP=('DX','DY','DZ'),),)
 426
 427
 428 CHSEIXM=CALC_CHAR_SEISME(
 429      MATR_MASS=MATRMASS,
 430      DIRECTION=(1., 0., 0.,),
 431      MODE_STAT=MSTA,
 432      NOEUD=('N02__I00','N01__I00','N0701I00',),
 433       )
 434
 435 CHSEIYM=CALC_CHAR_SEISME(
 436      MATR_MASS=MATRMASS,
 437      DIRECTION=(0., 1., 0.,),
 438      MODE_STAT=MSTA,
 439      NOEUD=('N02__I00','N01__I00','N0701I00',),
 440       )
 441
 442 CHSEIZM=CALC_CHAR_SEISME(
 443      MATR_MASS=MATRMASS,
 444      DIRECTION=(0., 0., 1.,),
 445      MODE_STAT=MSTA,
 446      NOEUD=('N02__I00','N01__I00','N1201I00',),
 447       )
 448
 449 MACRO_PROJ_BASE(
 450                BASE=MODES,
 451                MATR_ASSE_GENE=(
 452                         _F(
 453                         MATRICE=CO("MASGEN2"),
 454                         MATR_ASSE=MATRMASS,
 455                         ),
 456                         _F(
 457                         MATRICE=CO("RIGGEN2"),
 458                         MATR_ASSE=MATRRIGI,
 459                         ),
 460                ),
 461                VECT_ASSE_GENE=(
 462                         _F(
 463                         VECTEUR=CO("EFX2"),
 464                         VECT_ASSE=CHSEIXM,
 465                         ),
 466                         _F(
 467                         VECTEUR=CO("EFY2"),
 468                         VECT_ASSE=CHSEIYM,
 469                         ),
 470                         _F(
 471                         VECTEUR=CO("EFZ2"),
 472                         VECT_ASSE=CHSEIZM,
 473                         ),
 474                       ),
 475                      )
 476
 477
 478 VITESH1=CALC_FONCTION(
 479                   INTEGRE=_F(FONCTION=ACCELH1,METHODE='TRAPEZE',),
 480                   PROL_DROITE='CONSTANT'         )
 481 DEPLAH1=CALC_FONCTION(
 482                   INTEGRE=_F(FONCTION=VITESH1,METHODE='TRAPEZE',),
 483                   PROL_DROITE='CONSTANT'         )
 484 VITESV1=CALC_FONCTION(
 485                   INTEGRE=_F(FONCTION=ACCELV1,METHODE='TRAPEZE',),
 486                   PROL_DROITE='CONSTANT'         )
 487 DEPLAV1=CALC_FONCTION(
 488                   INTEGRE=_F(FONCTION=VITESV1,METHODE='TRAPEZE',),
 489                   PROL_DROITE='CONSTANT'         )
 490
 491 GSCA=DYNA_TRAN_MODAL(
 492                   MASS_GENE=MASGEN2,
 493                   RIGI_GENE=RIGGEN2,
 494                   AMOR_REDUIT=0.02,
 495                   MODE_STAT=MSTA,
 496                   EXCIT=(
 497                   _F(
 498                    VECT_GENE=EFX2,
 499                    ACCE=ACCELH1,
 500                    VITE=VITESH1,
 501                    DEPL=DEPLAH1,
 502                    MULT_APPUI='OUI',
 503                    DIRECTION=(1.,0.,0.,),
 504                    NOEUD=('N02__I00','N01__I00','N0701I00',),
 505                        ),
 506                   _F(
 507                    VECT_GENE=EFY2,
 508                    ACCE=ACCELH1,
 509                    VITE=VITESH1,
 510                    DEPL=DEPLAH1,
 511                    MULT_APPUI='OUI',
 512                    DIRECTION=(0.,1.,0.,),
 513                    NOEUD=('N02__I00','N01__I00','N0701I00',),
 514                        ),
 515                   _F(
 516                    VECT_GENE=EFZ2,
 517                    ACCE=ACCELV1,
 518                    VITE=VITESV1,
 519                    DEPL=DEPLAV1,
 520                    MULT_APPUI='OUI',
 521                    DIRECTION=(0.,0.,1.,),
 522                    NOEUD=('N02__I00','N01__I00','N1201I00',),
 523                        ),
 524                            ),
 525                   INCREMENT=_F(
 526                      INST_INIT = 0.,
 527                      INST_FIN =20.47,
 528                      PAS = 0.001
 529                            ),
 530                         )
 531
 532
 533 # RESTITUTION GRANDEURS ABSOLUES
 534 #-------------------------------
 535
 536 RGSCA=REST_BASE_PHYS(
 537                RESU_GENE=GSCA,
 538                INTERPOL='LIN',
 539                LIST_INST=L_RECU,
 540                TOUT_CHAM='OUI',
 541                MULT_APPUI='OUI',
 542                       )
 543
 544 # DEPLACEMENTS
 545
 546 DPSCA=CREA_CHAMP(
 547               TYPE_CHAM='NOEU_DEPL_R',
 548               OPERATION='EXTR',
 549               RESULTAT=RGSCA,
 550               NOM_CHAM='DEPL',
 551               TYPE_MAXI='MAXI',
 552               TYPE_RESU='VALE',
 553               TOUT_ORDRE='OUI'
 554                      )
 555
 556 TDPSCA=POST_RELEVE_T(
 557               ACTION=
 558               _F(
 559             INTITULE = 'DEPL_MAX_A',
 560             GROUP_NO='GNSTR01',
 561             CHAM_GD=DPSCA,
 562             TOUT_CMP='OUI',
 563             OPERATION ='EXTRACTION',
 564                      ),
 565             )
 566
 567 IMPR_TABLE(
 568           TABLE=TDPSCA,
 569          FORMAT='TABLEAU',
 570                )
 571
 572 # EFFORTS
 573
 574 RGSCA=CALC_ELEM(
 575                 reuse=RGSCA,
 576                 MODELE=MODELE,
 577                 CHAM_MATER=CHMATP,
 578                 CARA_ELEM=CARA,
 579                 OPTION=(
 580                 'EFGE_ELNO_DEPL',
 581                   ),
 582                 RESULTAT=RGSCA,)
 583
 584 EFSCA=CREA_CHAMP(
 585               TYPE_CHAM='ELNO_SIEF_R',
 586               OPERATION='EXTR',
 587               RESULTAT=RGSCA,
 588               NOM_CHAM='EFGE_ELNO_DEPL',
 589               TYPE_MAXI='MAXI',
 590               TYPE_RESU='VALE',
 591               TOUT_ORDRE='OUI'
 592                      )
 593
 594 TEFSCA=POST_RELEVE_T(
 595               ACTION=
 596               _F(
 597             INTITULE = 'EFGE_MAX_A',
 598             GROUP_NO='GNSTR01',
 599             CHAM_GD=EFSCA,
 600             TOUT_CMP='OUI',
 601             OPERATION ='EXTRACTION',
 602                      ),
 603             )
 604 IMPR_TABLE(
 605           TABLE=TEFSCA,
 606          FORMAT='TABLEAU',
 607                )
 608
 609
 610 # RESTITUTION GRANDEURS RELATIVES
 611 #-------------------------------
 612
 613 # DEPLACEMENTS
 614
 615 RGSCR2=REST_BASE_PHYS(
 616                RESU_GENE=GSCA,
 617                INTERPOL='LIN',
 618                LIST_INST=L_RECU,
 619                TOUT_CHAM='OUI',
 620                       )
 621
 622 DPSCR2=CREA_CHAMP(
 623               TYPE_CHAM='NOEU_DEPL_R',
 624               OPERATION='EXTR',
 625               RESULTAT=RGSCR2,
 626               NOM_CHAM='DEPL',
 627               TYPE_MAXI='MAXI',
 628               TYPE_RESU='VALE',
 629               TOUT_ORDRE='OUI'
 630                      )
 631
 632 TDPSCR2=POST_RELEVE_T(
 633               ACTION=
 634               _F(
 635             INTITULE = 'DEPL_MAX_R2',
 636             GROUP_NO='GNSTR01',
 637             CHAM_GD=DPSCR2,
 638             TOUT_CMP='OUI',
 639             OPERATION ='EXTRACTION',
 640                      ),
 641             )
 642
 643 IMPR_TABLE(
 644           TABLE=TDPSCR2,
 645          FORMAT='TABLEAU',
 646                )
 647
 648 # EFFORTS
 649
 650 RGSCR2=CALC_ELEM(
 651                 reuse=RGSCR2,
 652                 MODELE=MODELE,
 653                 CHAM_MATER=CHMATP,
 654                 CARA_ELEM=CARA,
 655                 OPTION=(
 656                 'EFGE_ELNO_DEPL',
 657                   ),
 658                 RESULTAT=RGSCR2,)
 659
 660 EFSCR2=CREA_CHAMP(
 661               TYPE_CHAM='ELNO_SIEF_R',
 662               OPERATION='EXTR',
 663               RESULTAT=RGSCR2,
 664               NOM_CHAM='EFGE_ELNO_DEPL',
 665               TYPE_MAXI='MAXI',
 666               TYPE_RESU='VALE',
 667               TOUT_ORDRE='OUI'
 668                      )
 669
 670 TEFSCR2=POST_RELEVE_T(
 671               ACTION=
 672               _F(
 673             INTITULE = 'EFGE_MAX_R2',
 674             GROUP_NO='GNSTR01',
 675             CHAM_GD=EFSCR2,
 676             TOUT_CMP='OUI',
 677             OPERATION ='EXTRACTION',
 678                      ),
 679             )
 680
 681 IMPR_TABLE(
 682           TABLE=TEFSCR2,
 683          FORMAT='TABLEAU',
 684                )
 685
 686 ###########################################################
 687 #------REPONSE SPECTRALE PAR SUPERPOSITION MODALE---------
 688 ###########################################################
 689
 690 # CALCUL DU SPECTRE DE REPONSE OSCILLATEUR EN PSEUDO ACCELERATION ABSOLUE
 691 # NAPPE SRO FONCTION DE AMORTISSEMENT
 692 #------------------------------------------------------------------------
 693 SROXY=CALC_FONCTION(
 694                   SPEC_OSCI=_F(
 695                      FONCTION=ACCELH1,
 696                      AMOR_REDUIT=(0.02,0.05,0.10),
 697                      NORME=9.81,
 698                      ),
 699                   )
 700 SROZ=CALC_FONCTION(
 701                   SPEC_OSCI=_F(
 702                      FONCTION=ACCELV1,
 703                      AMOR_REDUIT=(0.02,0.05,0.10),
 704                      NORME=9.81,
 705                      ),
 706                   )
 707
 708
 709 #-----------------------------------------------
 710 # REPONSE SPECTRALE SANS CORRECTION STATIQUE
 711 #-----------------------------------------------
 712 SPEC=COMB_SISM_MODAL(
 713                        MODE_MECA=MODES,
 714                        AMOR_REDUIT=0.02,
 715                        MASS_INER=MASSINER,
 716                        EXCIT=_F(
 717                            MONO_APPUI='OUI',
 718                            TRI_SPEC='OUI',
 719                            SPEC_OSCI=(SROXY,SROXY,SROZ),
 720                            ECHELLE=(9.81,9.81,9.81),),
 721                        COMB_MODE=_F(TYPE='CQC',),
 722                        COMB_DIRECTION=_F(  TYPE = 'QUAD',),
 723                        OPTION=(
 724                        'DEPL',
 725                        'EFGE_ELNO_DEPL',
 726                        ),
 727                         )
 728
 729 EFSPEC=POST_RELEVE_T(
 730               ACTION=
 731               _F(
 732             INTITULE = 'EFGE_SPEC',
 733             GROUP_NO='GNSTR01',
 734             RESULTAT=SPEC,
 735             NOM_CHAM='EFGE_ELNO_DEPL',
 736             NUME_ORDRE=4,
 737             TOUT_CMP='OUI',
 738             OPERATION ='EXTRACTION',
 739                      ),
 740             )
 741
 742 IMPR_TABLE(
 743           TABLE=EFSPEC,
 744          FORMAT='TABLEAU',
 745                )
 746
 747 DPSPEC=POST_RELEVE_T(
 748               ACTION=
 749               _F(
 750             INTITULE = 'DEPL_SPEC',
 751             GROUP_NO='GNSTR01',
 752             RESULTAT=SPEC,
 753             NOM_CHAM='DEPL',
 754             NUME_ORDRE=4,
 755             TOUT_CMP='OUI',
 756             OPERATION ='EXTRACTION',
 757                      ),
 758             )
 759
 760 IMPR_TABLE(
 761           TABLE=DPSPEC,
 762          FORMAT='TABLEAU',
 763                )
 764
 765
 766
 767 FIN()
 768
 769
 770
 771
 772
 773