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.
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.
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
22 DEBUT(CODE=_F( NOM = 'SDLL130B',NIV_PUB_WEB='INTERNET'))
23 # REPONSE SISMIQUE D'UNE POUTRE A COMPORTEMENT LINEAIRE (POU_D_EM)
28 # LECTURE MAILLAGE LINEIQUE
32 # LECTURE MAILLAGE SECTION1
33 PRE_GIBI(UNITE_GIBI=18, UNITE_MAILLAGE=21,);
34 MASEC=LIRE_MAILLAGE(UNITE=21,);
36 # ON ENRICHIT LE MAILLAGE LINEIQUE EN MULTIPLIANT LES LIGNES
37 MAPOU=CREA_MAILLAGE(MAILLAGE=POUT,
38 CREA_GROUP_MA=_F(NOM='POU_AC',
43 MAPOU=DEFI_GROUP(reuse =MAPOU,
45 CREA_GROUP_NO=_F(GROUP_MA='POUTRE',
49 MOPOU=AFFE_MODELE(MAILLAGE=MAPOU,
51 PHENOMENE='MECANIQUE',
52 MODELISATION='POU_D_EM',),);
57 BETON=DEFI_MATERIAU(ELAS=_F(E=3.7272000000E10,
60 AMOR_ALPHA=8.5000000000000006E-05,
61 AMOR_BETA=18.984999999999999,),);
63 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E11,
66 AMOR_ALPHA=8.5000000000000006E-05,
67 AMOR_BETA=18.984999999999999,),);
69 # LES CHAMPS DE MATERIAUX
73 CHMAT=AFFE_MATERIAU(MAILLAGE=MAPOU,
74 AFFE=(_F(GROUP_MA='POUTRE',
79 # ORIENTATIONS ET RIGIDITE DE TORSION
83 POUCA=AFFE_CARA_ELEM(MODELE=MOPOU,
84 POUTRE=_F(GROUP_MA=('POUTRE','POU_AC'),
87 PREC_INERTIE=2.948e-1,
89 VALE=(0.20000000000000001,0.5),),
90 ORIENTATION=_F(GROUP_MA=('POUTRE','POU_AC'),
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,
101 -0.066, 0.218, 8.E-3,),),
108 BLOCAGE=AFFE_CHAR_MECA(MODELE=MOPOU,
109 DDL_IMPO=(_F(GROUP_NO='A',
114 _F(GROUP_NO='TOUS_NO',
122 RIGI_ELE=CALC_MATR_ELEM(OPTION='RIGI_MECA',
128 MASS_ELE=CALC_MATR_ELEM(OPTION='MASS_MECA',
134 AMOR_ELE=CALC_MATR_ELEM(OPTION='AMOR_MECA',
145 NUME_DDL=NUME_DDL(MATR_RIGI=RIGI_ELE,
152 RIGIDITE=ASSE_MATRICE(MATR_ELEM=RIGI_ELE,
155 MASSE=ASSE_MATRICE(MATR_ELEM=MASS_ELE,
158 AMORT=ASSE_MATRICE(MATR_ELEM=AMOR_ELE,
161 MODE_MEC=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
163 CALC_FREQ=_F(OPTION='PLUS_PETITE',
165 # ON DEFINIT L'ACCELEROGRAMME DU SEISME
169 ACCELERO=CALC_FONCTION(COMB=_F(FONCTION=ACS2C11,
175 DIRSEISM=CALC_CHAR_SEISME(MATR_MASS=MASSE,
176 DIRECTION=(0.0,-1.0,0.0),
178 # CALCUL DYNAMIQUE LINEAIRE TRANSITOIRE
182 LINST=DEFI_LIST_REEL(DEBUT=0.0,
183 INTERVALLE=_F(JUSQU_A=5.0,
186 U1=DYNA_LINE_TRAN(MATR_MASS=MASSE,
190 EXCIT=_F(VECT_ASSE=DIRSEISM,
191 FONC_MULT=ACCELERO,),
192 INCREMENT=_F(LIST_INST=LINST,),);
194 U1=CALC_ELEM(reuse =U1,
200 OPTION='SIEF_ELGA_DEPL',
203 U1=CALC_NO(reuse =U1,
212 # QUELQUES TEST DE RESULTATS, REFERENCE SDLL130A
213 TEST_RESU(RESU=(_F(RESULTAT=MODE_MEC,
216 VALE=37.602600000000002,
217 PRECISION=6.0000000000000001E-3,
218 REFERENCE='AUTRE_ASTER',),
219 _F(RESULTAT=MODE_MEC,
223 PRECISION=1.0000000000000001E-2,
224 REFERENCE='AUTRE_ASTER',),
225 _F(RESULTAT=MODE_MEC,
228 VALE=197.09100000000001,
229 PRECISION=0.017000000000000001,
230 REFERENCE='AUTRE_ASTER',),
233 NOM_CHAM='REAC_NODA',
237 PRECISION=0.021999999999999999,
238 REFERENCE='AUTRE_ASTER',),
244 VALE=-6.0693599999999998E-4,
246 REFERENCE='AUTRE_ASTER',),
249 NOM_CHAM='REAC_NODA',
254 REFERENCE='AUTRE_ASTER',),
260 VALE=-2.3507300000000001E-3,
262 REFERENCE='AUTRE_ASTER',),
265 NOM_CHAM='REAC_NODA',
269 PRECISION=0.035000000000000003,
270 REFERENCE='AUTRE_ASTER',),
276 VALE=8.5790399999999996E-3,
277 PRECISION=0.021999999999999999,
278 REFERENCE='AUTRE_ASTER',),
281 NOM_CHAM='REAC_NODA',
285 PRECISION=0.029999999999999999,
286 REFERENCE='AUTRE_ASTER',),
292 VALE=-9.1083999999999991E-3,
294 REFERENCE='AUTRE_ASTER',),),);