src/TestsEltsVisu/Nicolas/sdll130b.comm

   1 # MODIF  DATE 08/11/2004   AUTEUR LEBOUVIE F.LEBOUVIER
   2 # TITRE  REPONSE SISMIQUE D'UNE POUTRE EN BETON ARME RECTANGULAIRE
   3 #            CONFIGURATION MANAGEMENT OF EDF VERSION
   4 # ======================================================================
   5 # COPYRIGHT (C) 1991 - 2002  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 #        ELEMENT POU_D_EM COMPORTEMENT LINEAIRE
  21 #
  22 DEBUT(CODE=_F( NOM = 'SDLL130B',NIV_PUB_WEB='INTERNET'))
  23 # REPONSE SISMIQUE D'UNE POUTRE A COMPORTEMENT LINEAIRE (POU_D_EM)
  24
  25 # LECTURE DU SEISME
  26 INCLUDE(UNITE=92,);
  27
  28 # LECTURE MAILLAGE LINEIQUE
  29 PRE_GIBI();
  30 POUT=LIRE_MAILLAGE();
  31
  32 # LECTURE MAILLAGE SECTION1
  33 PRE_GIBI(UNITE_GIBI=18, UNITE_MAILLAGE=21,);
  34 MASEC=LIRE_MAILLAGE(UNITE=21,);
  35
  36 # ON ENRICHIT LE MAILLAGE LINEIQUE EN MULTIPLIANT LES LIGNES
  37 MAPOU=CREA_MAILLAGE(MAILLAGE=POUT,
  38                     CREA_GROUP_MA=_F(NOM='POU_AC',
  39                                      TOUT='OUI',
  40                                      PREF_MAILLE='A',
  41                                      PREF_NUME=100,),);
  42
  43 MAPOU=DEFI_GROUP(reuse =MAPOU,
  44                  MAILLAGE=MAPOU,
  45                  CREA_GROUP_NO=_F(GROUP_MA='POUTRE',
  46                                   NOM='TOUS_NO',),
  47                  INFO=1,);
  48
  49 MOPOU=AFFE_MODELE(MAILLAGE=MAPOU,
  50                   AFFE=_F(TOUT='OUI',
  51                           PHENOMENE='MECANIQUE',
  52                           MODELISATION='POU_D_EM',),);
  53 # LES MATERIAUX
  54 #
  55
  56
  57 BETON=DEFI_MATERIAU(ELAS=_F(E=3.7272000000E10,
  58                             NU=0.0,
  59                             RHO=2400.0,
  60                             AMOR_ALPHA=8.5000000000000006E-05,
  61                             AMOR_BETA=18.984999999999999,),);
  62
  63 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E11,
  64                             NU=0.0,
  65                             RHO=7800.0,
  66                             AMOR_ALPHA=8.5000000000000006E-05,
  67                             AMOR_BETA=18.984999999999999,),);
  68 #
  69 # LES CHAMPS DE MATERIAUX
  70 #
  71
  72
  73 CHMAT=AFFE_MATERIAU(MAILLAGE=MAPOU,
  74                     AFFE=(_F(GROUP_MA='POUTRE',
  75                              MATER=BETON,),
  76                           _F(GROUP_MA='POU_AC',
  77                              MATER=ACIER,),),);
  78 #
  79 # ORIENTATIONS ET RIGIDITE DE TORSION
  80 #
  81
  82
  83 POUCA=AFFE_CARA_ELEM(MODELE=MOPOU,
  84                      POUTRE=_F(GROUP_MA=('POUTRE','POU_AC'),
  85                                SECTION='RECTANGLE',
  86                                PREC_AIRE=1.71e-2,
  87                                PREC_INERTIE=2.948e-1,
  88                                CARA=('HY','HZ'),
  89                                VALE=(0.20000000000000001,0.5),),
  90                      ORIENTATION=_F(GROUP_MA=('POUTRE','POU_AC'),
  91                                     CARA='ANGL_VRIL',
  92                                     VALE=-90.0,),
  93            AFFE_SECT = _F( GROUP_MA = 'POUTRE', NOM='BETON',
  94                            MAILLAGE_SECT = MASEC, TOUT_SECT = 'OUI',
  95                            COOR_AXE_POUTRE = (0., 0.,),),
  96            AFFE_FIBRE = _F( GROUP_MA = 'POU_AC', CARA = 'DIAMETRE',
  97                             NOM='ARMATURE' ,COOR_AXE_POUTRE = (0.,0.,),
  98                             VALE =( 0.066, -0.218, 32.E-3,
  99                                    -0.066, -0.218, 32.E-3,
 100                                     0.066, 0.218, 8.E-3,
 101                                    -0.066, 0.218, 8.E-3,),),
 102
 103                                     );
 104
 105
 106
 107
 108 BLOCAGE=AFFE_CHAR_MECA(MODELE=MOPOU,
 109                        DDL_IMPO=(_F(GROUP_NO='A',
 110                                     DX=0.0,
 111                                     DY=0.0,),
 112                                  _F(GROUP_NO='B',
 113                                     DY=0.0,),
 114                                  _F(GROUP_NO='TOUS_NO',
 115                                     DZ=0.0,
 116                                     DRX=0.0,
 117                                     DRY=0.0,),),);
 118
 119
 120
 121
 122 RIGI_ELE=CALC_MATR_ELEM(OPTION='RIGI_MECA',
 123                         MODELE=MOPOU,
 124                         CHAM_MATER=CHMAT,
 125                         CARA_ELEM=POUCA,
 126                         CHARGE=BLOCAGE,);
 127
 128 MASS_ELE=CALC_MATR_ELEM(OPTION='MASS_MECA',
 129                         MODELE=MOPOU,
 130                         CHAM_MATER=CHMAT,
 131                         CARA_ELEM=POUCA,
 132                         CHARGE=BLOCAGE,);
 133
 134 AMOR_ELE=CALC_MATR_ELEM(OPTION='AMOR_MECA',
 135                         MODELE=MOPOU,
 136                         CARA_ELEM=POUCA,
 137                         CHAM_MATER=CHMAT,
 138                         RIGI_MECA=RIGI_ELE,
 139                         MASS_MECA=MASS_ELE,
 140                         CHARGE=BLOCAGE,);
 141 #
 142 #
 143
 144
 145 NUME_DDL=NUME_DDL(MATR_RIGI=RIGI_ELE,
 146                   METHODE='LDLT',
 147                   RENUM='SANS',);
 148 #
 149 #
 150
 151
 152 RIGIDITE=ASSE_MATRICE(MATR_ELEM=RIGI_ELE,
 153                       NUME_DDL=NUME_DDL,);
 154
 155 MASSE=ASSE_MATRICE(MATR_ELEM=MASS_ELE,
 156                    NUME_DDL=NUME_DDL,);
 157
 158 AMORT=ASSE_MATRICE(MATR_ELEM=AMOR_ELE,
 159                    NUME_DDL=NUME_DDL,);
 160
 161 MODE_MEC=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
 162                           MATR_B=MASSE,
 163                           CALC_FREQ=_F(OPTION='PLUS_PETITE',
 164                                        NMAX_FREQ=3,),);
 165 # ON DEFINIT L'ACCELEROGRAMME DU SEISME
 166 #
 167
 168
 169 ACCELERO=CALC_FONCTION(COMB=_F(FONCTION=ACS2C11,
 170                                COEF=137.0,),);
 171 #
 172 #
 173
 174
 175 DIRSEISM=CALC_CHAR_SEISME(MATR_MASS=MASSE,
 176                           DIRECTION=(0.0,-1.0,0.0),
 177                           MONO_APPUI='OUI',);
 178 # CALCUL DYNAMIQUE LINEAIRE TRANSITOIRE
 179 #
 180
 181
 182 LINST=DEFI_LIST_REEL(DEBUT=0.0,
 183                      INTERVALLE=_F(JUSQU_A=5.0,
 184                                    NOMBRE=500,),);
 185
 186 U1=DYNA_LINE_TRAN(MATR_MASS=MASSE,
 187                   MATR_RIGI=RIGIDITE,
 188                   MATR_AMOR=AMORT,
 189                   NEWMARK=_F(),
 190                   EXCIT=_F(VECT_ASSE=DIRSEISM,
 191                            FONC_MULT=ACCELERO,),
 192                   INCREMENT=_F(LIST_INST=LINST,),);
 193
 194 U1=CALC_ELEM(reuse =U1,
 195              MODELE=MOPOU,
 196              CHAM_MATER=CHMAT,
 197              CARA_ELEM=POUCA,
 198              TOUT='OUI',
 199              TOUT_ORDRE='OUI',
 200              OPTION='SIEF_ELGA_DEPL',
 201              RESULTAT=U1,);
 202
 203 U1=CALC_NO(reuse =U1,
 204            MODELE=MOPOU,
 205            RESULTAT=U1,
 206            TOUT_ORDRE='OUI',
 207            OPTION='REAC_NODA',
 208            CHAM_MATER=CHMAT,
 209            CARA_ELEM=POUCA,
 210            TOUT='OUI',);
 211
 212 # QUELQUES TEST DE RESULTATS, REFERENCE SDLL130A
 213 TEST_RESU(RESU=(_F(RESULTAT=MODE_MEC,
 214                    NUME_MODE=1,
 215                    PARA='FREQ',
 216                    VALE=37.602600000000002,
 217                    PRECISION=6.0000000000000001E-3,
 218                    REFERENCE='AUTRE_ASTER',),
 219                 _F(RESULTAT=MODE_MEC,
 220                    NUME_MODE=2,
 221                    PARA='FREQ',
 222                    VALE=150.411,
 223                    PRECISION=1.0000000000000001E-2,
 224                    REFERENCE='AUTRE_ASTER',),
 225                 _F(RESULTAT=MODE_MEC,
 226                    NUME_MODE=3,
 227                    PARA='FREQ',
 228                    VALE=197.09100000000001,
 229                    PRECISION=0.017000000000000001,
 230                    REFERENCE='AUTRE_ASTER',),
 231                 _F(RESULTAT=U1,
 232                    NUME_ORDRE=100,
 233                    NOM_CHAM='REAC_NODA',
 234                    NOM_CMP='DY',
 235                    GROUP_NO='A',
 236                    VALE=1.8878E4,
 237                    PRECISION=0.021999999999999999,
 238                    REFERENCE='AUTRE_ASTER',),
 239                 _F(RESULTAT=U1,
 240                    NUME_ORDRE=100,
 241                    NOM_CHAM='DEPL',
 242                    NOM_CMP='DY',
 243                    GROUP_NO='C',
 244                    VALE=-6.0693599999999998E-4,
 245                    PRECISION=0.014,
 246                    REFERENCE='AUTRE_ASTER',),
 247                 _F(RESULTAT=U1,
 248                    NUME_ORDRE=200,
 249                    NOM_CHAM='REAC_NODA',
 250                    NOM_CMP='DY',
 251                    GROUP_NO='A',
 252                    VALE=6.3393E4,
 253                    PRECISION=0.02,
 254                    REFERENCE='AUTRE_ASTER',),
 255                 _F(RESULTAT=U1,
 256                    NUME_ORDRE=200,
 257                    NOM_CHAM='DEPL',
 258                    NOM_CMP='DY',
 259                    GROUP_NO='C',
 260                    VALE=-2.3507300000000001E-3,
 261                    PRECISION=1.E-2,
 262                    REFERENCE='AUTRE_ASTER',),
 263                 _F(RESULTAT=U1,
 264                    NUME_ORDRE=268,
 265                    NOM_CHAM='REAC_NODA',
 266                    NOM_CMP='DY',
 267                    GROUP_NO='A',
 268                    VALE=-2.32223E5,
 269                    PRECISION=0.035000000000000003,
 270                    REFERENCE='AUTRE_ASTER',),
 271                 _F(RESULTAT=U1,
 272                    NUME_ORDRE=268,
 273                    NOM_CHAM='DEPL',
 274                    NOM_CMP='DY',
 275                    GROUP_NO='C',
 276                    VALE=8.5790399999999996E-3,
 277                    PRECISION=0.021999999999999999,
 278                    REFERENCE='AUTRE_ASTER',),
 279                 _F(RESULTAT=U1,
 280                    NUME_ORDRE=468,
 281                    NOM_CHAM='REAC_NODA',
 282                    NOM_CMP='DY',
 283                    GROUP_NO='A',
 284                    VALE=2.46923E5,
 285                    PRECISION=0.029999999999999999,
 286                    REFERENCE='AUTRE_ASTER',),
 287                 _F(RESULTAT=U1,
 288                    NUME_ORDRE=468,
 289                    NOM_CHAM='DEPL',
 290                    NOM_CMP='DY',
 291                    GROUP_NO='C',
 292                    VALE=-9.1083999999999991E-3,
 293                    PRECISION=0.02,
 294                    REFERENCE='AUTRE_ASTER',),),);
 295
 296 FIN();
 297 #
 298 #