1 # MODIF DATE 31/08/2004 AUTEUR JMBHH01 J.M.PROIX
2 #-----------------------------------------------------------------------
3 # TITRE VALIDATION DU MODELE GRILLE EN ELASTICITE LINEAIRE
4 # CONFIGURATION MANAGEMENT OF EDF VERSION
5 # ======================================================================
6 # COPYRIGHT (C) 1991 - 2004 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.
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.
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 # PLAQUE D BETON AVEC DEUX NAPPES D'ARMATURE
22 #-----------------------------------------------------------------------
23 DEBUT( CODE=_F( NOM = 'SSLS109C',NIV_PUB_WEB='INTERNET') )
26 MAILL1=LIRE_MAILLAGE( )
40 MAILTOT=CREA_MAILLAGE(
53 # -----------------------------------------------------
54 # DEFINITION DES MATERIAUX
55 # -----------------------------------------------------
59 ACIER=DEFI_MATERIAU( ELAS=_F( E = 1.E11, NU = 0., ) ,
62 BETON=DEFI_MATERIAU( ELAS=_F( E = 1.E9, NU = 0., ) ,
66 # -------------------------------------------------
67 # AFFECTATION DES MATERIAUX
68 # -------------------------------------------------
70 MATAF=AFFE_MATERIAU( MAILLAGE=MAILTOT,AFFE=(
71 _F( GROUP_MA = ('GEOB'),
74 _F( GROUP_MA = ('GEOB1','GEOB2'),
81 # -------------------------------------------------
82 # AFFECTATION DES MODELES
83 # -------------------------------------------------
85 MODGRIL=AFFE_MODELE( MAILLAGE=MAILTOT, VERIF='MAILLE',
87 _F( GROUP_MA = ('GEOB1','GEOB2'),
88 MODELISATION = 'GRILLE',
89 PHENOMENE = 'MECANIQUE',
91 _F( GROUP_MA = ('GEOB','BORD1'),
93 PHENOMENE = 'MECANIQUE',
98 MODDKT=AFFE_MODELE( MAILLAGE=MAILTOT, VERIF='MAILLE',AFFE=(
99 _F( GROUP_MA = ('GEOB','GEOB1','GEOB2'),
100 MODELISATION = 'DKT',
101 PHENOMENE = 'MECANIQUE',
103 _F( GROUP_MA = ('BORD1'),
104 MODELISATION = 'DKT',
105 PHENOMENE = 'MECANIQUE',
113 # -------------------------------------------------
114 # AFFECTATION DES CARACTERISTIQUES
115 # -------------------------------------------------
117 CARAGRIL=AFFE_CARA_ELEM( MODELE=MODGRIL,
119 _F( GROUP_MA = ('GEOB', 'BORD1',),
121 ANGL_REP = (0.0, 0.0,),
125 _F( GROUP_MA = 'GEOB1',
127 ANGL_REP = (0.0, 0.0,),
130 _F( GROUP_MA = 'GEOB2',
132 ANGL_REP = (0.0, 0.0,),
137 CARADKT=AFFE_CARA_ELEM( MODELE=MODDKT,
139 _F( GROUP_MA = ('GEOB', 'BORD1',),
141 ANGL_REP = (0.0, 0.0,),
143 _F( GROUP_MA = 'GEOB1',
145 ANGL_REP = (0.0, 0.0,),
149 _F( GROUP_MA = 'GEOB2',
151 ANGL_REP = (0.0, 0.0,),
158 # -------------------------------------------------
159 # AFFECTATION DU CHARGEMENT
160 # -------------------------------------------------
162 CHARGRIL=AFFE_CHAR_MECA( MODELE=MODGRIL,DDL_IMPO=(
163 _F( NOEUD = ('NO1','NO4'),
164 DX = 0., DY = 0., DZ = 0.,
165 DRX = 0., DRY = 0., DRZ = 0.,
168 _F( NOEUD = ('NO3','NO5','NO2'), DRZ = 0., DRX=0., DY=0.,
171 FORCE_ARETE=_F(GROUP_MA='BORD1',FZ=1)
175 CHARDKT=AFFE_CHAR_MECA( MODELE=MODDKT,DDL_IMPO=(
176 _F( NOEUD = ('NO1','NO4'),
177 DX = 0., DY = 0., DZ = 0.,
178 DRX = 0., DRY = 0., DRZ = 0.,
181 _F( NOEUD = ('NO3','NO5','NO2'), DRZ = 0., DRX=0., DY=0.,
184 FORCE_ARETE=_F(GROUP_MA='BORD1',FZ=1)
188 # -------------------------------------------------
189 # CALCUL MECANIQUE STATIQUE LINEAIRE GRILLES
190 # -------------------------------------------------
193 RESGRIL1=MECA_STATIQUE( MODELE=MODGRIL, CHAM_MATER=MATAF,
195 EXCIT=_F( CHARGE = CHARGRIL, ) )
197 RESGRIL1=CALC_ELEM(reuse=RESGRIL1, RESULTAT=RESGRIL1,
198 OPTION='EFGE_ELNO_DEPL' )
200 RESGRIL1=CALC_ELEM(reuse=RESGRIL1, RESULTAT=RESGRIL1,
201 OPTION='SIGM_ELNO_DEPL' )
204 RESGRIL1=CALC_NO( reuse=RESGRIL1, RESULTAT=RESGRIL1,
209 RESGRIL1=CALC_ELEM(reuse=RESGRIL1, RESULTAT=RESGRIL1,
210 OPTION='EPSI_ELNO_DEPL' )
212 RESGRIL1=CALC_ELEM(reuse=RESGRIL1, RESULTAT=RESGRIL1,
213 OPTION='DEGE_ELNO_DEPL' )
215 RESGRIL1=CALC_ELEM(reuse=RESGRIL1, RESULTAT=RESGRIL1,
216 OPTION='SIEF_ELNO_ELGA' )
220 # --------------------------------------------------------
221 # CALCUL MECANIQUE STATIQUE LINEAIRE STAT_NON_LINE GRILLES
222 # --------------------------------------------------------
225 LINST=DEFI_LIST_REEL( DEBUT=0.,INTERVALLE=(
226 _F( JUSQU_A = 1.00, NOMBRE = 1, ),))
228 RESGRIL2=STAT_NON_LINE( MODELE=MODGRIL, CHAM_MATER=MATAF,
230 EXCIT=_F( CHARGE = CHARGRIL ) ,
231 NEWTON=_F(REAC_ITER=1,
232 # PREDICTION='ELASTIQUE'
234 INCREMENT=_F(LIST_INST = LINST),
235 COMP_INCR=_F(RELATION = 'ELAS' ,
238 RESGRIL2=CALC_NO( reuse=RESGRIL2, RESULTAT=RESGRIL2,
240 OPTION='FORC_NODA', )
242 RESGRIL2=CALC_ELEM(reuse=RESGRIL2, RESULTAT=RESGRIL2,
243 OPTION='SIEF_ELNO_ELGA' )
245 RESGRIL2=CALC_ELEM(reuse=RESGRIL2, RESULTAT=RESGRIL2,
246 OPTION='SIGM_ELNO_COQU' )
249 RESGRIL2=CALC_ELEM(reuse=RESGRIL2, RESULTAT=RESGRIL2,
250 OPTION='EPSI_ELNO_DEPL' )
252 RESGRIL2=CALC_ELEM(reuse=RESGRIL2, RESULTAT=RESGRIL2,
253 OPTION='DEGE_ELNO_DEPL' )
256 # -------------------------------------------------
257 # CALCUL MECANIQUE STATIQUE LINEAIRE DKT
258 # -------------------------------------------------
261 RESUDKT=MECA_STATIQUE( MODELE=MODDKT, CHAM_MATER=MATAF,
263 EXCIT=_F( CHARGE = CHARDKT, ) )
265 RESUDKT=CALC_ELEM(reuse=RESUDKT, RESULTAT=RESUDKT,
266 OPTION='EFGE_ELNO_DEPL' )
268 RESUDKT=CALC_ELEM(reuse=RESUDKT, RESULTAT=RESUDKT,
269 OPTION='SIGM_ELNO_DEPL' )
272 #RESUDKT=CALC_NO( reuse=RESUDKT, RESULTAT=RESUDKT,
274 # OPTION='FORC_NODA',)
277 RESUDKT=CALC_ELEM(reuse=RESUDKT, RESULTAT=RESUDKT,
278 OPTION='DEGE_ELNO_DEPL' )
280 #RESUDKT=CALC_ELEM(reuse=RESUDKT, RESULTAT=RESUDKT,
281 # OPTION='SIEF_ELNO_ELGA' )
286 IMPR_RESU(RESU=_F(RESULTAT=RESGRIL1))
288 IMPR_RESU(RESU=_F(RESULTAT=RESGRIL2))
290 IMPR_RESU(RESU=_F(RESULTAT=RESUDKT))
295 # -- TEST SUR LES DEPLACEMENTS ET LES ROTATIONS
303 REFERENCE='NON_REGRESSION',VALE= 1.56354E-07,
311 REFERENCE='NON_REGRESSION',VALE= 1.56354E-07,
319 REFERENCE='NON_REGRESSION',VALE= -2.18657E-07,
327 REFERENCE='NON_REGRESSION',VALE= -2.18657E-07,
331 # -- TEST SUR LES FORCES NODALES
336 NOM_CHAM = 'FORC_NODA',
339 REFERENCE='ANALYTIQUE',VALE= -0.5,
344 NOM_CHAM = 'FORC_NODA',
347 REFERENCE='ANALYTIQUE',VALE=0.5,
352 NOM_CHAM = 'FORC_NODA',
355 REFERENCE='ANALYTIQUE',VALE=0.5,
360 NOM_CHAM = 'FORC_NODA',
363 REFERENCE='ANALYTIQUE',VALE=-0.5,
368 NOM_CHAM = 'FORC_NODA',
371 REFERENCE='ANALYTIQUE',VALE=0.5,
376 NOM_CHAM = 'FORC_NODA',
379 REFERENCE='ANALYTIQUE',VALE=0.5,
382 # -- TEST SUR LES EFFORTS
387 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA4',
390 REFERENCE='NON_REGRESSION',VALE= -0.297,
395 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA4',
398 REFERENCE='NON_REGRESSION',VALE=-0.297,
404 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA1',
407 REFERENCE='NON_REGRESSION',VALE= -3.33876E-01,
412 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA3',
415 REFERENCE='NON_REGRESSION',VALE=-3.33876E-01,
422 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA4',
425 REFERENCE='NON_REGRESSION',VALE=4.01305E+00,
430 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA4',
433 REFERENCE='NON_REGRESSION',VALE=4.01305E+00,
438 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA4',
441 REFERENCE='NON_REGRESSION',VALE=-4.01305E+00,
446 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA4',
449 REFERENCE='NON_REGRESSION',VALE=-4.01305E+00,
455 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA1',
458 REFERENCE='NON_REGRESSION',VALE=4.50733E+00,
463 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA3',
466 REFERENCE='NON_REGRESSION',VALE=4.50733E+00,
471 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA1',
474 REFERENCE='NON_REGRESSION',VALE=-4.50733E+00,
479 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA3',
482 REFERENCE='NON_REGRESSION',VALE=-4.50733E+00,
498 # -- TEST SUR LES DEPLACEMENTS ET LES ROTATIONS
506 REFERENCE='NON_REGRESSION',VALE= 1.56354E-07,
514 REFERENCE='NON_REGRESSION',VALE= 1.56354E-07,
522 REFERENCE='NON_REGRESSION',VALE= -2.18657E-07,
530 REFERENCE='NON_REGRESSION',VALE= -2.18657E-07,
534 # -- TEST SUR LES FORCES NODALES
539 NOM_CHAM = 'FORC_NODA',
542 REFERENCE='ANALYTIQUE',VALE= -0.5,
547 NOM_CHAM = 'FORC_NODA',
550 REFERENCE='ANALYTIQUE',VALE=0.5,
555 NOM_CHAM = 'FORC_NODA',
558 REFERENCE='ANALYTIQUE',VALE=0.5,
563 NOM_CHAM = 'FORC_NODA',
566 REFERENCE='ANALYTIQUE',VALE=-0.5,
571 NOM_CHAM = 'FORC_NODA',
574 REFERENCE='ANALYTIQUE',VALE=0.5,
579 NOM_CHAM = 'FORC_NODA',
582 REFERENCE='ANALYTIQUE',VALE=0.5,
585 # -- TEST SUR LES EFFORTS
590 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA4',
593 REFERENCE='NON_REGRESSION',VALE= -0.297,
598 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA4',
601 REFERENCE='NON_REGRESSION',VALE=-0.297,
607 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA1',
610 REFERENCE='NON_REGRESSION',VALE= -3.33876E-01,
615 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='MA3',
618 REFERENCE='NON_REGRESSION',VALE=-3.33876E-01,
625 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA4',
628 REFERENCE='NON_REGRESSION',VALE=4.01305E+00,
633 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA4',
636 REFERENCE='NON_REGRESSION',VALE=4.01305E+00,
641 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA4',
644 REFERENCE='NON_REGRESSION',VALE=-4.01305E+00,
649 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA4',
652 REFERENCE='NON_REGRESSION',VALE=-4.01305E+00,
658 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA1',
661 REFERENCE='NON_REGRESSION',VALE=4.50733E+00,
666 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='BMA3',
669 REFERENCE='NON_REGRESSION',VALE=4.50733E+00,
674 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA1',
677 REFERENCE='NON_REGRESSION',VALE=-4.50733E+00,
682 NOM_CHAM = 'SIEF_ELNO_ELGA',MAILLE='CMA3',
685 REFERENCE='NON_REGRESSION',VALE=-4.50733E+00,
694 # TEST DKT + EXCENTREMENT
699 # -- TEST SUR LES DEPLACEMENTS ET LES ROTATIONS
700 # -- LA DIFFERENCE SUR LES DEPLACEMENTS VIENT DE LA ROTATION
701 # -- PROPRE, PRISE EN COMPTE POUR DKT MAIS PAS POUR GRILLE
708 NOM_CMP = 'DZ', PRECISION=0.03,
709 REFERENCE='AUTRE_ASTER',VALE= 1.56354E-07,
716 NOM_CMP = 'DZ', PRECISION=0.03,
717 REFERENCE='AUTRE_ASTER',VALE= 1.56354E-07,
725 REFERENCE='AUTRE_ASTER',VALE= -2.18657E-07,
733 REFERENCE='AUTRE_ASTER',VALE= -2.18657E-07,
736 # -- TEST SUR LES EFFORTS
741 NOM_CHAM = 'EFGE_ELNO_DEPL',MAILLE='MA4',
744 REFERENCE='NON_REGRESSION',VALE= -0.297, PRECISION=0.08,
749 NOM_CHAM = 'EFGE_ELNO_DEPL',MAILLE='MA4',
752 REFERENCE='NON_REGRESSION',VALE=-0.297, PRECISION=0.08,
758 NOM_CHAM = 'EFGE_ELNO_DEPL',MAILLE='MA1',
761 REFERENCE='NON_REGRESSION',VALE= -3.33876E-01, PRECISION=0.05,
766 NOM_CHAM = 'EFGE_ELNO_DEPL',MAILLE='MA3',
769 REFERENCE='NON_REGRESSION',VALE=-3.33876E-01,PRECISION=0.05,
782 # -- TEST SUR LES FORCES NODALES : FAUSSES EN DKT + EXCENTREMENT
783 # 70% D'ECART, EQUILIBRE NON VERIFIE
790 NOM_CHAM = 'FORC_NODA',
793 REFERENCE='ANALYTIQUE',VALE= -0.5,
798 NOM_CHAM = 'FORC_NODA',
801 REFERENCE='ANALYTIQUE',VALE=0.5,
806 NOM_CHAM = 'FORC_NODA',
809 REFERENCE='ANALYTIQUE',VALE=0.5,
814 NOM_CHAM = 'FORC_NODA',
817 REFERENCE='ANALYTIQUE',VALE=-0.5,
822 NOM_CHAM = 'FORC_NODA',
825 REFERENCE='ANALYTIQUE',VALE=0.5,
830 NOM_CHAM = 'FORC_NODA',
833 REFERENCE='ANALYTIQUE',VALE=0.5,