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.
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 # RESPONSABLE BOYERE E.BOYERE
21 #-----------------------------------------------------------------------
23 # ANALYSE RIS12 TR4 REGLEMENTAIRE
27 #-----------------------------------------------------------------------
29 DEBUT(CODE=_F(NOM ='FORMA12C',NIV_PUB_WEB='INTERNET'))
31 MAILLAGE=LIRE_MAILLAGE(INFO=2)
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',),
48 IMPR_RESU(FORMAT='CASTEM',UNITE=37,
55 #-----------------------------------------------------------------------
56 # CARACTERISTIQUES DES ELEMENTS
57 #-----------------------------------------------------------------------
63 _F(GROUP_MA='GMEL0301',
66 VALE=(0.08415,0.00711),),
67 _F(GROUP_MA='GMEL0501',
70 VALE=(0.08415,0.0034),),
71 _F(GROUP_MA='GMEL0601',
74 VALE=(0.08415,0.0034),),
75 _F(GROUP_MA='GMEL0701',
78 VALE=(0.08415,0.0034),),
79 _F(GROUP_MA='GMEL1001',
82 VALE=(0.08415,0.0034),),
83 _F(GROUP_MA='GMEL1101',
86 VALE=(0.08415,0.0034),),
87 _F(GROUP_MA='GMEL1201',
90 VALE=(0.08415,0.0034),),
91 _F(GROUP_MA='GMEL1501',
94 VALE=(0.08415,0.0034),),
95 _F(GROUP_MA='GMEL1701',
98 VALE=(0.08415,0.0034),),
99 _F(GROUP_MA='GMEL1801',
102 VALE=(0.08415,0.0034),),
103 _F(GROUP_MA='GMEL1901',
106 VALE=(0.08415,0.0034),),
107 _F(GROUP_MA='GMEL2001',
110 VALE=(0.08415,0.0034),),
111 _F(GROUP_MA='GMEL2101',
114 VALE=(0.08415,0.0034),),
115 _F(GROUP_MA='GMEL2201',
118 VALE=(0.08415,0.0034),),
119 _F(GROUP_MA='GMEL2401',
122 VALE=(0.08415,0.00711),),),
123 DISCRET=(_F(MAILLE='M0801I02',
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',
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,),),)
155 #-----------------------------------------------------------------------
156 # DEFINITION DES MATERIAUX
157 #-----------------------------------------------------------------------
159 MATP01=DEFI_MATERIAU(ELAS=_F(E=1.97E11,
161 RHO=1.3108749019999999E4,
162 ALPHA=1.639999937E-05,),)
163 MATP03=DEFI_MATERIAU(ELAS=_F(E=1.97E11,
166 ALPHA=1.639999937E-05,),)
168 #-----------------------------------------------------------------------
169 # DEFINITION DES SUPPORTS
170 #-----------------------------------------------------------------------
172 SUPPORT=AFFE_CHAR_MECA(
176 _F( NOEUD = ('N02__I00', 'N01__I00',),
177 DX = 0., DY = 0., DZ = 0., DRX = 0., DRY = 0., DRZ = 0.),
179 _F( NOEUD = 'N0701I00',
182 _F( NOEUD = 'N1201I00',
187 #-----------------------------------------------------------------------
188 # DEFINITION DES CHAMPS DE MATERIAUX
189 #-----------------------------------------------------------------------
191 CHMATP=AFFE_MATERIAU(MAILLAGE=MAILLAGE,
192 AFFE=(_F(GROUP_MA='GMAT01',
195 _F(GROUP_MA='GMAT03',
199 #-----------------------------------------------------------------------
200 # ASSEMBLAGE MATRICE MASSE ET RIGIDITE
201 #-----------------------------------------------------------------------
203 MACRO_MATR_ASSE(MODELE=MODELE,
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',),),);
214 #-----------------------------------------------------------------------
215 # ASSEMBLAGE DU SECOND MEMBRE
216 #-----------------------------------------------------------------------
220 CHSEIX=CALC_CHAR_SEISME( MATR_MASS=MATRMASS,
221 DIRECTION=(1., 0., 0.,), MONO_APPUI='OUI'
223 CHSEIY=CALC_CHAR_SEISME( MATR_MASS=MATRMASS,
224 DIRECTION=(0., 1., 0.,), MONO_APPUI='OUI'
226 CHSEIZ=CALC_CHAR_SEISME( MATR_MASS=MATRMASS,
227 DIRECTION=(0., 0., 1.,), MONO_APPUI='OUI'
230 #-----------------------------------------------------------------------
232 #-----------------------------------------------------------------------
238 MASS_INER=_F( TOUT = 'OUI',)
240 MODES=MACRO_MODE_MECA(
247 VERI_MODE=_F( STOP_ERREUR = 'NON'),
251 IMPRESSION=_F( TOUT_PARA = 'OUI')
254 # TEST_RESU UNIQUEMENT POUR FAIRE CAS TEST
255 TEST_RESU(RESU=(_F(RESULTAT=MODES,
258 VALE= 4.5135794612152,
260 REFERENCE='NON_REGRESSION',
273 ###########################################################
274 #---- REPONSE TRANSITOIRE PAR SUPERPOSITION MODALE ------
275 #---- MONO APPUI GRANDEURS RELATIVES ------
276 ###########################################################
278 # DEFINITION DES INSTANTS DE RECUPERATION
280 L_RECU=DEFI_LIST_REEL( DEBUT=0.,
286 #----------------------------------------------------
287 #------------- PROJECTION SUR LA BASE MODALE --------
288 #----------------------------------------------------
293 MATRICE=CO("MASSEGEN"),
297 MATRICE=CO("RIGIDGEN"),
303 VECTEUR=CO("EFGENX"),
307 VECTEUR=CO("EFGENY"),
311 VECTEUR=CO("EFGENZ"),
317 #------------------------------------------------------------------
318 #-------------- RESOLUTION DU PB GENERALISE -----------------------
319 #-------------- SANS CORRECTION STATIQUE --------------------------
320 #------------------------------------------------------------------
321 GSCR=DYNA_TRAN_MODAL(
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
352 RGSCR=REST_GENE_PHYS(
360 TYPE_CHAM='NOEU_DEPL_R',
368 TDPSCR=POST_RELEVE_T(
371 INTITULE = 'DEPL_MAX_R',
375 OPERATION ='EXTRACTION',
380 # RESTITUTION GLOBALE DU CHAMP EFGE_ELNO_DEPL A CHAQUE INSTANT
392 TYPE_CHAM='ELNO_SIEF_R',
395 NOM_CHAM='EFGE_ELNO_DEPL',
400 TEFSCR=POST_RELEVE_T(
403 INTITULE = 'EFGE_MAX_R',
407 OPERATION ='EXTRACTION',
415 ###########################################################
416 #---- REPONSE TRANSITOIRE PAR SUPERPOSITION MODALE ------
417 #---- MULTI APPUIS GRANDEURS ABSOLUES ------
418 ###########################################################
420 # DEFINITION EXCITATION MULTIPLE
425 MODE_STAT=_F(TOUT='OUI',AVEC_CMP=('DX','DY','DZ'),),)
428 CHSEIXM=CALC_CHAR_SEISME(
430 DIRECTION=(1., 0., 0.,),
432 NOEUD=('N02__I00','N01__I00','N0701I00',),
435 CHSEIYM=CALC_CHAR_SEISME(
437 DIRECTION=(0., 1., 0.,),
439 NOEUD=('N02__I00','N01__I00','N0701I00',),
442 CHSEIZM=CALC_CHAR_SEISME(
444 DIRECTION=(0., 0., 1.,),
446 NOEUD=('N02__I00','N01__I00','N1201I00',),
453 MATRICE=CO("MASGEN2"),
457 MATRICE=CO("RIGGEN2"),
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' )
491 GSCA=DYNA_TRAN_MODAL(
503 DIRECTION=(1.,0.,0.,),
504 NOEUD=('N02__I00','N01__I00','N0701I00',),
512 DIRECTION=(0.,1.,0.,),
513 NOEUD=('N02__I00','N01__I00','N0701I00',),
521 DIRECTION=(0.,0.,1.,),
522 NOEUD=('N02__I00','N01__I00','N1201I00',),
533 # RESTITUTION GRANDEURS ABSOLUES
534 #-------------------------------
536 RGSCA=REST_GENE_PHYS(
547 TYPE_CHAM='NOEU_DEPL_R',
556 TDPSCA=POST_RELEVE_T(
559 INTITULE = 'DEPL_MAX_A',
563 OPERATION ='EXTRACTION',
585 TYPE_CHAM='ELNO_SIEF_R',
588 NOM_CHAM='EFGE_ELNO_DEPL',
594 TEFSCA=POST_RELEVE_T(
597 INTITULE = 'EFGE_MAX_A',
601 OPERATION ='EXTRACTION',
610 # RESTITUTION GRANDEURS RELATIVES
611 #-------------------------------
615 RGSCR2=REST_GENE_PHYS(
623 TYPE_CHAM='NOEU_DEPL_R',
632 TDPSCR2=POST_RELEVE_T(
635 INTITULE = 'DEPL_MAX_R2',
639 OPERATION ='EXTRACTION',
661 TYPE_CHAM='ELNO_SIEF_R',
664 NOM_CHAM='EFGE_ELNO_DEPL',
670 TEFSCR2=POST_RELEVE_T(
673 INTITULE = 'EFGE_MAX_R2',
677 OPERATION ='EXTRACTION',
686 ###########################################################
687 #------REPONSE SPECTRALE PAR SUPERPOSITION MODALE---------
688 ###########################################################
690 # CALCUL DU SPECTRE DE REPONSE OSCILLATEUR EN PSEUDO ACCELERATION ABSOLUE
691 # NAPPE SRO FONCTION DE AMORTISSEMENT
692 #------------------------------------------------------------------------
696 AMOR_REDUIT=(0.02,0.05,0.10),
703 AMOR_REDUIT=(0.02,0.05,0.10),
709 #-----------------------------------------------
710 # REPONSE SPECTRALE SANS CORRECTION STATIQUE
711 #-----------------------------------------------
712 SPEC=COMB_SISM_MODAL(
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',),
729 EFSPEC=POST_RELEVE_T(
732 INTITULE = 'EFGE_SPEC',
735 NOM_CHAM='EFGE_ELNO_DEPL',
738 OPERATION ='EXTRACTION',
747 DPSPEC=POST_RELEVE_T(
750 INTITULE = 'DEPL_SPEC',
756 OPERATION ='EXTRACTION',