src/TestsEltsVisu/Nicolas/forma11d.comm

   1 # MODIF  DATE 13/02/2006   AUTEUR DURAND C.DURAND
   2 # TITRE TP ANALYSE MODALE TPAM4
   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.
  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 # RESPONSABLE BOYERE E.BOYERE
  21 #-------------------------------------------------------------------
  22 #                 CORRIGE TP ANALYSE MODAL N4
  23 #                            MODE_ITER_INV
  24 #
  25 #-------------------------------------------------------------------
  26 #directory = "/home/xxx/ASTER/TP_MODAL/TPAM4/"
  27
  28 #-------------------------------------------------------------------
  29 #                       CALCUL ASTER
  30 #-------------------------------------------------------------------
  31 DEBUT(CODE=_F(NOM ='FORMA11D',NIV_PUB_WEB='INTERNET'))
  32
  33 # LECTURE MAILLAGE ET AFFECTATION MATERIAU/MODELE
  34 MAYA=LIRE_MAILLAGE();
  35 MODL=AFFE_MODELE(MAILLAGE=MAYA,
  36                  AFFE=(_F(GROUP_MA='MANCHONS',
  37                           PHENOMENE='MECANIQUE',
  38                           MODELISATION='POU_D_T',),
  39                        _F(GROUP_MA='CUVE_TUB',
  40                           PHENOMENE='MECANIQUE',
  41                           MODELISATION='DKT',),
  42                        _F(GROUP_MA='SUPGUID',
  43                           PHENOMENE='MECANIQUE',
  44                           MODELISATION='DKT',),
  45                        _F(GROUP_MA='RENF_SUP',
  46                           PHENOMENE='MECANIQUE',
  47                           MODELISATION='DKT',),),);
  48 Z = 0.0;
  49 IMANCH=(2.6876E-04/64);
  50 JMANCH=(2.6876E-04/32);
  51 SMANCH=(2.8098E-01/64);
  52
  53 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.10000000000E11,
  54                             NU=0.29999999999999999,
  55                             RHO=7850.0,),);
  56 MANCHON=DEFI_MATERIAU(ELAS=_F(E=3.15000000000E11,
  57                               NU=0.29999999999999999,
  58                               RHO=7850.0,),);
  59 CHMAT=AFFE_MATERIAU(MAILLAGE=MAYA,
  60                     AFFE=(_F(GROUP_MA='CUVE_TUB',
  61                              MATER=ACIER,),
  62                           _F(GROUP_MA='SUPGUID',
  63                              MATER=ACIER,),
  64                           _F(GROUP_MA='RENF_SUP',
  65                              MATER=ACIER,),
  66                           _F(GROUP_MA='MANCHONS',
  67                              MATER=MANCHON,),),);
  68 PARAM=AFFE_CARA_ELEM(MODELE=MODL,
  69                      POUTRE=_F(GROUP_MA='MANCHONS',
  70                                SECTION='GENERALE',
  71                                CARA=('A','IZ','IY','JX','AY','AZ',),
  72                                VALE=(SMANCH,IMANCH,IMANCH,JMANCH,Z,Z,),),
  73                      COQUE=(_F(GROUP_MA='CUVE_128',
  74                                EPAIS=0.128,),
  75                             _F(GROUP_MA='CUVE_060',
  76                                EPAIS=0.059999999999999998,),
  77                             _F(GROUP_MA='CUVE_200',
  78                                EPAIS=0.20000000000000001,),
  79                             _F(GROUP_MA='CUVE_230',
  80                                EPAIS=0.23000000000000001,),
  81                             _F(GROUP_MA='CUVE_450',
  82                                EPAIS=0.45000000000000001,),
  83                             _F(GROUP_MA='CUVE_160',
  84                                EPAIS=0.16,),
  85                             _F(GROUP_MA='CUVE_252',
  86                                EPAIS=0.252,),
  87                             _F(GROUP_MA='CUVE_264',
  88                                EPAIS=0.26400000000000001,),
  89                             _F(GROUP_MA='CUVE_166',
  90                                EPAIS=0.16600000000000001,),
  91                             _F(GROUP_MA='CUVE_170',
  92                                EPAIS=0.17000000000000001,),
  93                             _F(GROUP_MA='CUVE_81',
  94                                EPAIS=0.081000000000000003,),
  95                             _F(GROUP_MA='CUVE_77',
  96                                EPAIS=0.076999999999999999,),
  97                             _F(GROUP_MA='SUPGUID',
  98                                EPAIS=0.1143,),
  99                             _F(GROUP_MA='RENF_SUP',
 100                                EPAIS=0.042999999999999997,),),);
 101
 102 # AFFECTATION CL
 103 CONDLIM=AFFE_CHAR_MECA(MODELE=MODL,
 104                        DDL_IMPO=(_F(GROUP_NO='SUPP_CUV',
 105                                     DX=0.0,
 106                                     DY=0.0,
 107                                     DZ=0.0,
 108                                     DRX=0.0,
 109                                     DRY=0.0,
 110                                     DRZ=0.0,),
 111                                  _F(GROUP_NO='ENCASS1',
 112                                     DX=0.0,
 113                                     DY=0.0,
 114                                     DZ=0.0,
 115                                     DRX=0.0,
 116                                     DRY=0.0,
 117                                     DRZ=0.0,),),);
 118
 119 # CALCULS DES MATRICES DE MASSE ET DE RIGIDITE
 120 K_ELEM1=CALC_MATR_ELEM(OPTION='RIGI_MECA',
 121                        MODELE=MODL,
 122                        CHAM_MATER=CHMAT,
 123                        CARA_ELEM=PARAM,
 124                        CHARGE=CONDLIM,);
 125 M_ELEM1=CALC_MATR_ELEM(OPTION='MASS_MECA',
 126                        MODELE=MODL,
 127                        CHAM_MATER=CHMAT,
 128                        CARA_ELEM=PARAM,
 129                        CHARGE=CONDLIM,);
 130 NUM1=NUME_DDL(MATR_RIGI=K_ELEM1,);
 131 K_ASSE1=ASSE_MATRICE(MATR_ELEM=K_ELEM1,
 132                      NUME_DDL=NUM1,);
 133 M_ASSE1=ASSE_MATRICE(MATR_ELEM=M_ELEM1,
 134                      NUME_DDL=NUM1,);
 135
 136 #--------------------------------------------------------------------
 137 # ANALYSE MODALE
 138 # QUESTION 1 : CALCUL BANDE PAR SORENSEN VIA MODE_ITER_SIMULT
 139
 140
 141 MODES=MODE_ITER_SIMULT(MATR_A=K_ASSE1,
 142                        MATR_B=M_ASSE1,
 143                        CALC_FREQ=_F(OPTION='BANDE',
 144                                     FREQ=(1.0,117.0,),),);
 145
 146 # TEST_RESU UNIQUEMENT POUR FAIRE CAS TEST
 147 TEST_RESU(RESU=(_F(RESULTAT=MODES,
 148                    NUME_ORDRE=1,
 149                    PARA='FREQ',
 150                    VALE=46.7756,
 151                    PRECISION=1e-05,
 152                    REFERENCE='NON_REGRESSION',
 153                    VERSION='7.3.11',)))
 154
 155 IMPR_RESU(RESU=_F(RESULTAT=MODES,
 156                   TOUT_CHAM='NON',
 157                   TOUT_PARA='OUI',),);
 158
 159 # IMPRESSION FORMAT CASTEM DES MODES PROPRES
 160 IMPR_RESU(FORMAT='CASTEM',UNITE=55,
 161           RESU=_F(
 162                   MAILLAGE=MAYA,
 163                   RESULTAT=MODES,),);
 164
 165 # AFFICHAGE GIBI INTERACTIF VIA PYTHON
 166 #os.system('cp %spost.dcas .' %directory)
 167 #os.system('/logiciels/aster/outils/gibi2000.x post.dcas')
 168
 169 # FIN DE LA QUESTION 1.
 170 #FIN();
 171
 172 # QUESTION 2: CALCUL OPTION PROCHE VIA MODE_ITER_INV
 173
 174 MODEINV=MODE_ITER_INV(MATR_A=K_ASSE1,
 175                       MATR_B=M_ASSE1,
 176                       CALC_FREQ=_F(OPTION='PROCHE',
 177                                    FREQ=106.0,),);
 178 IMPR_RESU(RESU=_F(RESULTAT=MODEINV,
 179                   TOUT_CHAM='NON',
 180                   TOUT_PARA='OUI',),);
 181
 182 FIN();