]> SALOME platform Git repositories - modules/eficas.git/commitdiff
Salome HOME
PN : mise a jour pour include_materiau
authoreficas <>
Fri, 30 Apr 2004 14:47:21 +0000 (14:47 +0000)
committereficas <>
Fri, 30 Apr 2004 14:47:21 +0000 (14:47 +0000)
Aster/Tests/Recette/efica01a.11
Aster/Tests/Recette/efica01a.com0
Aster/Tests/Recette/efica01a.comm

index eef038b5932ae2ab08bee35696aaa3f6b38a1902..e323531ee4eac9ca8c6404eadea540433ddae37e 100755 (executable)
@@ -2,9 +2,9 @@ DEPL2 = FORMULE(REEL="""(REEL:INST) =
 SIN(OMEGAA*INST)/(OMEGAA**2) """)
 DEPLACE2=CALC_FONC_INTERP( FONCTION=DEPL2,
                            LIST_PARA=L_INST,
+                          NOM_PARA='INST',
                            PROL_DROITE='LINEAIRE',
                            PROL_GAUCHE='LINEAIRE',
-                           NOM_PARA='INST',
                            NOM_RESU='DEPL'       )
 
 MONO_X=CALC_CHAR_SEISME(  MATR_MASS=MASSE,
@@ -18,4 +18,3 @@ MACRO_PROJ_BASE(BASE=MODE_MEC,MATR_ASSE_GENE=(
                 VECT_ASSE_GENE=_F( VECTEUR = CO("VECT_X"), VECT_ASSE = MONO_X)
                )
 
-
index bad17ff2395af15730073bd938ad0357d4dbb742..214e3b8920289f4b0c056980ea573e79154caab4 100755 (executable)
@@ -1,5 +1,13 @@
 POURSUITE(CODE=_F(  NOM = 'EFICA01A'))
 
+INCLUDE_MATERIAU(   NOM_AFNOR='18MND5',      TYPE_MODELE='REF',
+                      VARIANTE='A',          TYPE_VALE='NOMI',
+                     NOM_MATER='MAT3',       INFO=2               )
+
+CHMAT3=AFFE_MATERIAU(  MAILLAGE=MAILLA3,
+                       AFFE=_F(  TOUT='OUI', MATER = MAT3, TEMP_REF = 20.))
+
+riginor = 2.88E7 
 
 TRAN_GE2=DYNA_TRAN_MODAL(  MASS_GENE=MGEN_BIC,   RIGI_GENE=RGEN_BIC,
                  METHODE='EULER',
@@ -20,7 +28,7 @@ TRAN_GE2=DYNA_TRAN_MODAL(  MASS_GENE=MGEN_BIC,   RIGI_GENE=RGEN_BIC,
                        NORM_OBST = (0., 0., 1.,),
                        DIST_1 = 0.4495,
                        DIST_2 = 0.4495,
-                       RIGI_NOR = 2.88E7,
+                       RIGI_NOR = riginor,
                        AMOR_NOR = 0.,
                        RIGI_TAN = 0.,
                        COULOMB = 0.),
@@ -38,40 +46,4 @@ RESUA1=REST_BASE_PHYS(  RESU_GENE=TRAN_GE1,
                         MULT_APPUI='OUI'
                        )
 
-RESUR1=REST_BASE_PHYS(  RESU_GENE=TRAN_GE1,
-                        LIST_INST=LISTIMP,
-                        INTERPOL='LIN',
-                        NOM_CHAM='DEPL',
-#                        MULT_APPUI='NON'
-                       )
-
-
-NO2_DXA=RECU_FONCTION(  NOEUD='NO2',  NOM_CMP='DX',  NOM_CHAM='DEPL',
-                        RESULTAT=RESUA1,  
-                        LIST_INST=LISTIMP,   INTERPOL='LIN' )
-
-NO2_DX_A=RECU_FONCTION(  NOEUD='NO2',  NOM_CMP='DX',  NOM_CHAM='DEPL',
-                         RESU_GENE=TRAN_GE2,  
-                         MULT_APPUI='OUI',        INTERP_NUME='LIN',
-                         LIST_INST=LISTIMP,   INTERPOL='LIN' )
-
-DIFF=CALC_FONCTION(COMB=(
-               _F(  FONCTION = NO2_DX_A,  COEF = 1.),
-               _F(  FONCTION = NO2_DXA,   COEF = -1.))
-                    )
-
-TEST_FONCTION(VALEUR=(
-    _F(  FONCTION = DIFF, VALE_PARA = 0.10, VALE_REFE = 0.,
-             CRITERE = 'ABSOLU', PRECISION = 3.E-5),
-    _F(  FONCTION = DIFF, VALE_PARA = 0.30, VALE_REFE = 0.,
-             CRITERE = 'ABSOLU', PRECISION = 3.E-5),
-    _F(  FONCTION = DIFF, VALE_PARA = 0.50, VALE_REFE = 0.,
-             CRITERE = 'ABSOLU', PRECISION = 3.E-5),
-    _F(  FONCTION = DIFF, VALE_PARA = 0.70, VALE_REFE = 0.,
-             CRITERE = 'ABSOLU', PRECISION = 3.E-5),
-    _F(  FONCTION = DIFF, VALE_PARA = 1.00, VALE_REFE = 0.,
-             CRITERE = 'ABSOLU', PRECISION = 3.E-5))
-             )
-
-
 FIN()
index d3555c27768c6e1a0afbff44ba1d497c7f54ce91..46507be0482aeaac7d2dc4b924c3c673893dbcb3 100755 (executable)
@@ -1,10 +1,21 @@
-# MODIF  DATE 16/01/2002   AUTEUR DURAND C.DURAND 
+# MODIF  DATE 13/01/2004   AUTEUR DURAND C.DURAND 
 # TITRE TEST DE NON REGRESSION DE L IHM EFICAS - DERIVE DE SDND102A
 #            CONFIGURATION MANAGEMENT OF EDF VERSION
 # ======================================================================
 # COPYRIGHT (C) 1991 - 2001  EDF R&D                  WWW.CODE-ASTER.ORG
-#              SEE THE FILE "LICENSE.TERMS" FOR INFORMATION ON USAGE AND
-#              REDISTRIBUTION OF THIS FILE.
+# THIS PROGRAM IS FREE SOFTWARE YOU CAN REDISTRIBUTE IT AND/OR MODIFY
+# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
+# THE FREE SOFTWARE FOUNDATION EITHER VERSION 2 OF THE LICENSE, OR   
+# (AT YOUR OPTION) ANY LATER VERSION.                                 
+#
+# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT 
+# WITHOUT ANY WARRANTY WITHOUT EVEN THE IMPLIED WARRANTY OF          
+# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU    
+# GENERAL PUBLIC LICENSE FOR MORE DETAILS.                            
+#
+# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE   
+# ALONG WITH THIS PROGRAM IF NOT, WRITE TO EDF R&D CODE_ASTER,       
+#    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.      
 # ======================================================================
 # Ce cas test est gere en configuration dans la base ASTER, il sert de
 # modele pour prononcer la recette de l IHM d EFICAS : l objectif est de
 # - verifier qu on peut bien supprimer une commande, un mot cle simple et facteur
 # - verifier les acces a la doc
 #
+#
 
-DEBUT()
-
-MAILLAGE=LIRE_MAILLAGE( )
 
-MAILLA2=LIRE_MAILLAGE(  UNITE=19 )
+DEBUT(CODE=_F(NOM='EFICA01A',NIV_PUB_WEB='INTERNET'),)
 
-MODELE=AFFE_MODELE(  MAILLAGE=MAILLAGE,
-                     AFFE=(
-              _F(  PHENOMENE = 'MECANIQUE', MODELISATION = 'DIS_T',
-                   GROUP_MA = 'RESSORT'),
-              _F(  PHENOMENE = 'MECANIQUE', MODELISATION = 'DIS_T',
-                   GROUP_NO = 'MASSES')    )
-                    )
+MAILLAGE=LIRE_MAILLAGE()
 
-BICHOC=AFFE_MODELE(  MAILLAGE=MAILLA2,
-                     AFFE=(
-              _F(  PHENOMENE = 'MECANIQUE', MODELISATION = 'DIS_T',
-                   GROUP_MA = 'RESSORTS'),
-              _F(  PHENOMENE = 'MECANIQUE', MODELISATION = 'DIS_T',
-                   GROUP_NO = ( 'MASSES1', 'MASSES2',)))
-                    )
+MAILLAGE=DEFI_GROUP(reuse=MAILLAGE,
+                    MAILLAGE=MAILLAGE,
+                    CREA_GROUP_NO=(_F(NOEUD='NO2',
+                                      NOM='MASSES',),
+                                   _F(NOEUD='NO1',
+                                      NOM='ENCASTRE',),),)
 
-VAL = 98696.0
+MAILLA2=LIRE_MAILLAGE(UNITE=21,)
 
-CARA_ELE=AFFE_CARA_ELEM(  MODELE=MODELE,
-                          DISCRET=(
-           _F(  CARA = 'K_T_D_L',  REPERE = 'GLOBAL', GROUP_MA = 'RESSORT',
-                VALE = (VAL, 0., 0., )),
-           _F(  CARA = 'M_T_D_N',  GROUP_NO = 'MASSES', VALE = 25.0))
-                         )
-
-CARA_BIC=AFFE_CARA_ELEM(  MODELE=BICHOC,
-                          DISCRET=(
-           _F(  CARA = 'K_T_D_L',  REPERE = 'GLOBAL', GROUP_MA = 'RESSORTS',
-                VALE = (VAL, 0., 0., )),
-           _F(  CARA = 'M_T_D_N',  GROUP_NO = 'MASSES1', VALE = 25.0),
-           _F(  CARA = 'M_T_D_N',  GROUP_NO = 'MASSES2', VALE = 25.0))
-                         )
-
-
-CON_LIM=AFFE_CHAR_MECA(  MODELE=MODELE,DDL_IMPO=(
-              _F(  GROUP_NO = 'ENCASTRE',  DX = 0.,  DY = 0.,  DZ = 0.),
-              _F(  GROUP_NO = 'MASSES',              DY = 0.,  DZ = 0.))
-                       )
-
-CL_BICHO=AFFE_CHAR_MECA(  MODELE=BICHOC,DDL_IMPO=(
-              _F(  GROUP_NO = 'ENCBICHO',     DX = 0.,  DY = 0.,  DZ = 0.),
-              _F(  GROUP_NO = ( 'MASSES1', 'MASSES2',), DY = 0.,  DZ = 0.))
-                       )
-
-MACRO_MATR_ASSE(  MODELE=MODELE,
-                  CHARGE=CON_LIM,
-                  CARA_ELEM=CARA_ELE,
-                  NUME_DDL=CO("NUMEDDL"),
-                  MATR_ASSE=(
-                 _F(  MATRICE = CO("RIGIDITE"),  OPTION = 'RIGI_MECA'),
-                 _F(  MATRICE = CO("MASSE"),     OPTION = 'MASS_MECA'))
-               )
-
-MACRO_MATR_ASSE(  MODELE=BICHOC,
-                  CHARGE=CL_BICHO,
-                  CARA_ELEM=CARA_BIC,
-                  NUME_DDL=CO("NUMDDLC"),
-                  MATR_ASSE=(
-                 _F(  MATRICE = CO("RIGI_BIC"),  OPTION = 'RIGI_MECA'),
-                 _F(  MATRICE = CO("MASS_BIC"),  OPTION = 'MASS_MECA'))
-               )
-
-MODE_MEC=MODE_ITER_SIMULT(  MATR_A=RIGIDITE,   MATR_B=MASSE )
-
-MODE_MEC=NORM_MODE(reuse=MODE_MEC,  MODE=MODE_MEC,   NORME='MASS_GENE' )
-
-MODE_BIC=MODE_ITER_SIMULT(  MATR_A=RIGI_BIC,   MATR_B=MASS_BIC,
-                  METHODE='JACOBI',
-                  OPTION='SANS',
-                  CALC_FREQ=_F(  OPTION = 'BANDE',  FREQ = (1., 10., ))
-                            )
-
-MODE_BIC=NORM_MODE(reuse=MODE_BIC,  MODE=MODE_BIC,   NORME='MASS_GENE' )
-
-MODE_STA=MODE_STATIQUE(  MATR_RIGI=RIGIDITE,   MATR_MASS=MASSE,
-                          MODE_STAT=_F(  TOUT = 'OUI', AVEC_CMP = 'DX') )
-
-MSTA_BIC=MODE_STATIQUE(  MATR_RIGI=RIGI_BIC,   MATR_MASS=MASS_BIC,
-                          MODE_STAT=_F(  TOUT = 'OUI', AVEC_CMP = 'DX') )
-
-
-L_INST=DEFI_LIST_REEL(  DEBUT=0.,
-                           INTERVALLE=_F(  JUSQU_A = 1., PAS = 0.0001) )
-
-OMEGAA=DEFI_VALEUR( R8=EVAL("""2.*PI*10.""") )
-
-ACCE1 = FORMULE(REEL="""(REEL:INST) = SIN(OMEGAA*INST) """)
-ACCELER1=CALC_FONC_INTERP( FONCTION=ACCE1, LIST_PARA=L_INST,
-                               PROL_DROITE='LINEAIRE',
-                              PROL_GAUCHE='LINEAIRE',
-                              NOM_PARA='INST',
-                               NOM_RESU='ACCE'       )
-
-ACCE2 = FORMULE(REEL="""(REEL:INST) =- SIN(OMEGAA*INST) """)
-ACCELER2=CALC_FONC_INTERP( FONCTION=ACCE2, LIST_PARA=L_INST,
-                               PROL_DROITE='LINEAIRE',
-                              PROL_GAUCHE='LINEAIRE',
-                              NOM_PARA='INST',
-                               NOM_RESU='ACCE'       )
-
-VITE1 = FORMULE(REEL="""(REEL:INST) =-COS(OMEGAA*INST)/OMEGAA """)
-VITESSE1=CALC_FONC_INTERP( FONCTION=VITE1, LIST_PARA=L_INST,
-                       PROL_DROITE='LINEAIRE',
-                      PROL_GAUCHE='LINEAIRE',
-                              NOM_PARA='INST',
-                      NOM_RESU='VITE'       )
-
-DEPL1 = FORMULE(REEL="""(REEL:INST) =-SIN(OMEGAA*INST)/(OMEGAA**2) """)
-DEPLACE1=CALC_FONC_INTERP( FONCTION=DEPL1, LIST_PARA=L_INST,
-                       PROL_DROITE='LINEAIRE',
-                      PROL_GAUCHE='LINEAIRE',
-                              NOM_PARA='INST',
-                      NOM_RESU='DEPL'       )
-
-VITE2 = FORMULE(REEL="""(REEL:INST) =COS(OMEGAA*INST)/OMEGAA """)
-VITESSE2=CALC_FONC_INTERP( FONCTION=VITE2, LIST_PARA=L_INST,
-                       PROL_DROITE='LINEAIRE',
-                      PROL_GAUCHE='LINEAIRE',
-                              NOM_PARA='INST',
-                      NOM_RESU='VITE'       )
-
-INCLUDE(   UNITE=11,   INFO=1 )
-
-MUR=DEFI_OBSTACLE(   TYPE='PLAN_Z' )
-
-TRAN_GE1=DYNA_TRAN_MODAL(  MASS_GENE=MASS_GEN,   RIGI_GENE=RIGI_GEN,
-                METHODE='EULER',
-                AMOR_REDUIT=0.07, 
-                MODE_STAT=MODE_STA,
-                EXCIT=_F(
-                       VECT_GENE = VECT_X,
-                       ACCE = ACCELER1,
-                       VITE = VITESSE1,
-                       DEPL = DEPLACE1,
-                       MULT_APPUI = 'OUI',
-                       DIRECTION = ( 1., 0., 0.,),
-                       NOEUD = 'NO1'),
-                CHOC=_F(  GROUP_NO_1 = 'MASSES',
-                      OBSTACLE = MUR,
-                      INTITULE = 'NO2/MUR',
-                      ORIG_OBST = ( -1., 0., 0., ),
-                      NORM_OBST = (  0., 0., 1., ),
-                      JEU = 1.1005,
-                      RIGI_NOR = 5.76E7,
-                      AMOR_NOR = 0.,
-                      RIGI_TAN = 0.,
-                      COULOMB = 0.0),
-                INCREMENT=_F( INST_INIT = 0.,  INST_FIN = 1.,  PAS = 0.0002),
-                ARCHIVAGE=_F(  PAS_ARCH = 8)
-                          )
-
-
-MULT_X1=CALC_CHAR_SEISME(  MATR_MASS=MASS_BIC,  DIRECTION=( 1., 0., 0.,),
-                           MODE_STAT=MSTA_BIC,  NOEUD='NO1' )
-
-MULT_X2=CALC_CHAR_SEISME(  MATR_MASS=MASS_BIC,  DIRECTION=( 1., 0., 0.,),
-                           MODE_STAT=MSTA_BIC,  NOEUD='NO11' )
+MODELE=AFFE_MODELE(MAILLAGE=MAILLAGE,
+                   AFFE=(_F(GROUP_MA='RESSORT',
+                            PHENOMENE='MECANIQUE',
+                            MODELISATION='DIS_T',),
+                         _F(GROUP_NO='MASSES',
+                            PHENOMENE='MECANIQUE',
+                            MODELISATION='DIS_T',),),)
 
-MACRO_PROJ_BASE(BASE=MODE_BIC,
-                MATR_ASSE_GENE=(
-                _F( MATRICE = CO("MGEN_BIC"), MATR_ASSE = MASS_BIC),
-                _F( MATRICE = CO("RGEN_BIC"), MATR_ASSE = RIGI_BIC)),
-                VECT_ASSE_GENE=(
-                _F( VECTEUR = CO("VECT_X1"),  VECT_ASSE = MULT_X1),
-                _F( VECTEUR = CO("VECT_X2"),  VECT_ASSE = MULT_X2))
-               )
+BICHOC=AFFE_MODELE(MAILLAGE=MAILLA2,
+                   AFFE=(_F(GROUP_MA='RESSORTS',
+                            PHENOMENE='MECANIQUE',
+                            MODELISATION='DIS_T',),
+                         _F(GROUP_NO=('MASSES1','MASSES2',),
+                            PHENOMENE='MECANIQUE',
+                            MODELISATION='DIS_T',),),)
+VAL = 98696.0
 
 
-GRILLE=DEFI_OBSTACLE(   TYPE='BI_PLAN_Z' )
+CARA_ELE=AFFE_CARA_ELEM(MODELE=MODELE,
+                        DISCRET=(_F(GROUP_MA='RESSORT',
+                                    REPERE='GLOBAL',
+                                    CARA='K_T_D_L',
+                                    VALE=(VAL,0.0,0.0,),),
+                                 _F(GROUP_NO='MASSES',
+                                    CARA='M_T_D_N',
+                                    VALE=25.0,),),)
+
+CARA_BIC=AFFE_CARA_ELEM(MODELE=BICHOC,
+                        DISCRET=(_F(GROUP_MA='RESSORTS',
+                                    REPERE='GLOBAL',
+                                    CARA='K_T_D_L',
+                                    VALE=(VAL,0.0,0.0,),),
+                                 _F(GROUP_NO='MASSES1',
+                                    CARA='M_T_D_N',
+                                    VALE=25.0,),
+                                 _F(GROUP_NO='MASSES2',
+                                    CARA='M_T_D_N',
+                                    VALE=25.0,),),)
+
+CON_LIM=AFFE_CHAR_MECA(MODELE=MODELE,
+                       DDL_IMPO=(_F(GROUP_NO='ENCASTRE',
+                                    DX=0.0,
+                                    DY=0.0,
+                                    DZ=0.0,),
+                                 _F(GROUP_NO='MASSES',
+                                    DY=0.0,
+                                    DZ=0.0,),),)
+
+CL_BICHO=AFFE_CHAR_MECA(MODELE=BICHOC,
+                        DDL_IMPO=(_F(GROUP_NO='ENCBICHO',
+                                     DX=0.0,
+                                     DY=0.0,
+                                     DZ=0.0,),
+                                  _F(GROUP_NO=('MASSES1','MASSES2',),
+                                     DY=0.0,
+                                     DZ=0.0,),),)
+
+MACRO_MATR_ASSE(MODELE=MODELE,
+                CARA_ELEM=CARA_ELE,
+                CHARGE=CON_LIM,
+                NUME_DDL=CO('NUMEDDL'),
+                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
+                              OPTION='RIGI_MECA',),
+                           _F(MATRICE=CO('MASSE'),
+                              OPTION='MASS_MECA',),),)
+
+MACRO_MATR_ASSE(MODELE=BICHOC,
+                CARA_ELEM=CARA_BIC,
+                CHARGE=CL_BICHO,
+                NUME_DDL=CO('NUMDDLC'),
+                MATR_ASSE=(_F(MATRICE=CO('RIGI_BIC'),
+                              OPTION='RIGI_MECA',),
+                           _F(MATRICE=CO('MASS_BIC'),
+                              OPTION='MASS_MECA',),),)
+
+MODE_MEC=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
+                          MATR_B=MASSE,)
+
+MODE_MEC=NORM_MODE(reuse =MODE_MEC,
+                   MODE=MODE_MEC,
+                   NORME='MASS_GENE',)
+
+MODE_BIC=MODE_ITER_SIMULT(MATR_A=RIGI_BIC,
+                          MATR_B=MASS_BIC,
+                          METHODE='JACOBI',
+                          OPTION='SANS',
+                          CALC_FREQ=_F(OPTION='BANDE',
+                                       FREQ=(1.0,10.0,),),)
+
+MODE_BIC=NORM_MODE(reuse =MODE_BIC,
+                   MODE=MODE_BIC,
+                   NORME='MASS_GENE',)
+
+MODE_STA=MODE_STATIQUE(MATR_RIGI=RIGIDITE,
+                       MATR_MASS=MASSE,
+                       MODE_STAT=_F(TOUT='OUI',
+                                    AVEC_CMP='DX',),)
+
+MSTA_BIC=MODE_STATIQUE(MATR_RIGI=RIGI_BIC,
+                       MATR_MASS=MASS_BIC,
+                       MODE_STAT=_F(TOUT='OUI',
+                                    AVEC_CMP='DX',),)
+
+L_INST=DEFI_LIST_REEL(DEBUT=0.0,
+                      INTERVALLE=_F(JUSQU_A=1.0,
+                                    PAS=1.E-4,),)
+
+OMEGAA=DEFI_VALEUR(R8=EVAL("""2.*PI*10."""),)
+
+ACCE1 = FORMULE(REEL = '''(REEL:INST) = SIN(OMEGAA*INST)''')
+
+ACCELER1=CALC_FONC_INTERP(FONCTION=ACCE1,
+                          NOM_RESU='ACCE',
+                          LIST_PARA=L_INST,
+                          NOM_PARA = 'INST',
+                          PROL_DROITE='LINEAIRE',
+                          PROL_GAUCHE='LINEAIRE',)
+
+ACCE2 = FORMULE(REEL = '''(REEL:INST) = - SIN(OMEGAA*INST)''')
+
+ACCELER2=CALC_FONC_INTERP(FONCTION=ACCE2,
+                          NOM_RESU='ACCE',
+                          LIST_PARA=L_INST,
+                          NOM_PARA = 'INST',
+                          PROL_DROITE='LINEAIRE',
+                          PROL_GAUCHE='LINEAIRE',)
+
+VITE1 = FORMULE(REEL = '''(REEL:INST) = -COS(OMEGAA*INST)/OMEGAA''')
+
+VITESSE1=CALC_FONC_INTERP(FONCTION=VITE1,
+                          NOM_RESU='VITE',
+                          LIST_PARA=L_INST,
+                          NOM_PARA = 'INST',
+                          PROL_DROITE='LINEAIRE',
+                          PROL_GAUCHE='LINEAIRE',)
+
+DEPL1 = FORMULE(REEL = '''(REEL:INST) = -SIN(OMEGAA*INST)/(OMEGAA**2)''')
+
+DEPLACE1=CALC_FONC_INTERP(FONCTION=DEPL1,
+                          NOM_RESU='DEPL',
+                          LIST_PARA=L_INST,
+                          NOM_PARA = 'INST',
+                          PROL_DROITE='LINEAIRE',
+                          PROL_GAUCHE='LINEAIRE',)
+
+VITE2 = FORMULE(REEL = '''(REEL:INST) = COS(OMEGAA*INST)/OMEGAA''')
+
+VITESSE2=CALC_FONC_INTERP(FONCTION=VITE2,
+                          NOM_RESU='VITE',
+                          LIST_PARA=L_INST,
+                          NOM_PARA = 'INST',
+                          PROL_DROITE='LINEAIRE',
+                          PROL_GAUCHE='LINEAIRE',)
+
+INCLUDE(UNITE=11,
+        INFO=1,)
+
+MUR=DEFI_OBSTACLE(TYPE='PLAN_Z',)
+
+TRAN_GE1=DYNA_TRAN_MODAL(METHODE='EULER',
+                         MASS_GENE=MASS_GEN,
+                         RIGI_GENE=RIGI_GEN,
+                         AMOR_REDUIT=0.07,
+                         MODE_STAT=MODE_STA,
+                         INCREMENT=_F(INST_INIT=0.0,
+                                      INST_FIN=1.0,
+                                      PAS=2.0E-4,),
+                         ARCHIVAGE=_F(PAS_ARCH=8,),
+                         EXCIT=_F(VECT_GENE=VECT_X,
+                                  ACCE=ACCELER1,
+                                  VITE=VITESSE1,
+                                  DEPL=DEPLACE1,
+                                  MULT_APPUI='OUI',
+                                  DIRECTION=(1.0,0.0,0.0,),
+                                  NOEUD='NO1',),
+                         CHOC=_F(INTITULE='NO2/MUR',
+                                 GROUP_NO_1='MASSES',
+                                 OBSTACLE=MUR,
+                                 ORIG_OBST=(-1.0,0.0,0.0,),
+                                 NORM_OBST=(0.0,0.0,1.0,),
+                                 JEU=1.1005,
+                                 RIGI_NOR=5.76E7,
+                                 AMOR_NOR=0.0,
+                                 RIGI_TAN=0.0,
+                                 COULOMB=0.0,),)
+
+MULT_X1=CALC_CHAR_SEISME(MATR_MASS=MASS_BIC,
+                         DIRECTION=(1.0,0.0,0.0,),
+                         MODE_STAT=MSTA_BIC,
+                         NOEUD='NO1',)
+
+MULT_X2=CALC_CHAR_SEISME(MATR_MASS=MASS_BIC,
+                         DIRECTION=(1.0,0.0,0.0,),
+                         MODE_STAT=MSTA_BIC,
+                         NOEUD='NO11',)
 
+MACRO_PROJ_BASE(BASE=MODE_BIC,
+                MATR_ASSE_GENE=(_F(MATRICE=CO('MGEN_BIC'),
+                                   MATR_ASSE=MASS_BIC,),
+                                _F(MATRICE=CO('RGEN_BIC'),
+                                   MATR_ASSE=RIGI_BIC,),),
+                VECT_ASSE_GENE=(_F(VECTEUR=CO('VECT_X1'),
+                                   VECT_ASSE=MULT_X1,),
+                                _F(VECTEUR=CO('VECT_X2'),
+                                   VECT_ASSE=MULT_X2,),),)
+
+GRILLE=DEFI_OBSTACLE(TYPE='BI_PLAN_Z',)
+
+MAILLA3=LIRE_MAILLAGE(UNITE=22,)
 
 FIN()