src/TestsEltsVisu/Nicolas/forma01b.comm

   1 # MODIF  DATE 06/07/2004   AUTEUR CIBHHLV L.VIVAN
   2 # RESPONSABLE JMBHH01 J.M.PROIX
   3 #
   4 #
   5 DEBUT(CODE=_F(NOM='FORMA01B',
   6               NIV_PUB_WEB='INTERNET',),);
   7 # TITRE TRAVAUX PRATIQUES FORMATION DE BASE
   8 #            CONFIGURATION MANAGEMENT OF EDF VERSION
   9 # ======================================================================
  10 # COPYRIGHT (C) 1991 - 2001  EDF R&D                  WWW.CODE-ASTER.ORG
  11 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
  12 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
  13 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
  14 # (AT YOUR OPTION) ANY LATER VERSION.
  15 #
  16 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  17 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  18 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  19 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  20 #
  21 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
  22 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
  23 #    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
  24 # ======================================================================
  25 #        TUYAUTERIE COMPORTANT UN COUDE
  26 #     MODELISATION PAR DES ELEMENTS COQUES DKT
  27 #     ----------------------------------------
  28 #
  29 #
  30 #
  31 #
  32
  33
  34 PRE_GIBI();
  35
  36 MAIL=LIRE_MAILLAGE();
  37 #
  38 # DEFINITION DES ELEMENTS FINIS UTILISES
  39 #
  40 #
  41 #
  42 #
  43
  44
  45 MODMECA=AFFE_MODELE(MAILLAGE=MAIL,
  46                     AFFE=_F(GROUP_MA=('TUYAU','EFOND',),
  47                             PHENOMENE='MECANIQUE',
  48                             MODELISATION='DKT',),);
  49 # ORIENTATION DES NORMALES AUX COQUES RENTRANTES DANS LE TUYAU
  50 #
  51 #
  52 #
  53
  54
  55 MAIL=MODI_MAILLAGE(reuse =MAIL,
  56                    MAILLAGE=MAIL,
  57                    ORIE_NORM_COQUE=_F(GROUP_MA='TUYAU',
  58                                       VECT_NORM=(1.0,0.0,0.0,),
  59                                       GROUP_NO='A1',),
  60                    MODELE=MODMECA,);
  61 # DEFINITION DU MATERIAU UTILISE
  62 #
  63 #
  64 #
  65 #
  66
  67
  68 ACIER=DEFI_MATERIAU(ELAS=_F(E=204000000000.0,
  69                             NU=0.3,),);
  70
  71 CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
  72                     AFFE=_F(TOUT='OUI',
  73                             MATER=ACIER,),);
  74 #
  75 # CARACTERISTIQUES DES COQUES
  76 #
  77 #
  78 #
  79 #
  80
  81
  82 CARA_COQ=AFFE_CARA_ELEM(MODELE=MODMECA,
  83                         COQUE=_F(GROUP_MA=('TUYAU','EFOND',),
  84                                  EPAIS=0.02,
  85                                  ANGL_REP=(0.0,-90.0,),),);
  86 #
  87 # DEFINITION DES CONDITIONS AUX LIMITES
  88 #    ENCASTREMENT SUR LA SECTION ENCAST
  89 #
  90 #
  91 #
  92 #
  93
  94
  95 BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODMECA,
  96                         DDL_IMPO=(_F(GROUP_MA='ENCAST',
  97                                      DX=0.0,
  98                                      DY=0.0,
  99                                      DZ=0.0,
 100                                      DRX=0.0,
 101                                      DRY=0.0,
 102                                      DRZ=0.0,),
 103                                   _F(GROUP_MA='SYMETRIE',
 104                                      DZ=0.0,
 105                                      DRX=0.0,
 106                                      DRY=0.0,),),);
 107 #
 108 # DEFINITION DES CHARGEMENTS
 109 #    CHARGEMENT 1 : FY=100KN (FX,FZ,MX,MY,MZ NULS)
 110 #    CHARGEMENT 2 : PRESSION INTERNE DE 15 MPA
 111 #
 112 #
 113 #
 114 #
 115
 116 FYTOT = 100000.0;
 117
 118 EPTUB = 0.02;
 119
 120 REXT = 0.2;
 121
 122
 123 RMOY=REXT-(EPTUB/2)
 124
 125 FYREP=FYTOT/2./pi/RMOY
 126
 127 RINT=RMOY-(EPTUB/2)
 128
 129 SINT=pi*(RINT*RINT)
 130 #
 131 # 2*PI*R*FY = 100 KN
 132 #
 133 #
 134 #
 135 #
 136
 137
 138 CHARG1=AFFE_CHAR_MECA(MODELE=MODMECA,
 139                       FORCE_ARETE=_F(GROUP_MA='EFOND',
 140                                      FY=FYREP,),);
 141 #
 142 # 2*PI*R*FX = PINT * SINT ;
 143 #
 144 #
 145 #
 146
 147 PINT = 15000000.0;
 148
 149
 150 PRES=PINT*RINT/RMOY
 151
 152 CHARG2=AFFE_CHAR_MECA(MODELE=MODMECA,
 153                       FORCE_COQUE=_F(GROUP_MA='TUYAU',
 154                                      PRES=PRES,),);
 155
 156 RESU1=MECA_STATIQUE(MODELE=MODMECA,
 157                     CHAM_MATER=CHMAT,
 158                     CARA_ELEM=CARA_COQ,
 159                     EXCIT=(_F(CHARGE=BLOQUAGE,),
 160                            _F(CHARGE=CHARG1,),),);
 161
 162 RESU1=CALC_ELEM(reuse =RESU1,
 163                 OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
 164                 RESULTAT=RESU1,
 165                 NIVE_COUCHE='SUP',);
 166
 167 RESU2=MECA_STATIQUE(MODELE=MODMECA,
 168                     CHAM_MATER=CHMAT,
 169                     CARA_ELEM=CARA_COQ,
 170                     EXCIT=(_F(CHARGE=BLOQUAGE,),
 171                            _F(CHARGE=CHARG2,),),);
 172
 173 RESU2=CALC_ELEM(reuse =RESU2,
 174                 OPTION='SIGM_ELNO_DEPL',
 175                 RESULTAT=RESU2,
 176                 NIVE_COUCHE='SUP',);
 177
 178 RESU2=CALC_ELEM(reuse =RESU2,
 179                 OPTION='EQUI_ELNO_SIGM',
 180                 RESULTAT=RESU2,);
 181
 182 IMPR_RESU(MODELE=MODMECA,
 183           RESU=(_F(RESULTAT=RESU1,
 184                    NOM_CHAM='DEPL',
 185                    GROUP_MA='EFOND',),
 186                 _F(RESULTAT=RESU2,
 187                    NOM_CHAM='DEPL',
 188                    GROUP_MA='EFOND',),),);
 189
 190 IMPR_RESU(MODELE=MODMECA,FORMAT='CASTEM',
 191           RESU=(_F(MAILLAGE=MAIL,),
 192                 _F(RESULTAT=RESU1,
 193                    NOM_CHAM=('DEPL','SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),),
 194                 _F(RESULTAT=RESU2,
 195                    NOM_CHAM=('DEPL','SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),),),);
 196
 197 MAIL=DEFI_GROUP(reuse =MAIL,
 198                 MAILLAGE=MAIL,
 199                 CREA_GROUP_NO=(_F(GROUP_MA='EFOND',),
 200                                _F(OPTION='PLAN',
 201                                   NOM='PEXT',
 202                                   POINT=(3.6,3.6,0.0,),
 203                                   VECT_NORMALE=(0.0,1.0,0.0,),
 204                                   PRECISION=0.001,),
 205                                _F(INTERSEC=('EFOND','PEXT',),
 206                                   NOM='PB',),),);
 207
 208 TABDEP1=POST_RELEVE_T(ACTION=_F(INTITULE='DEPB1',
 209                                 GROUP_NO='PB',
 210                                 RESULTAT=RESU1,
 211                                 NOM_CHAM='DEPL',
 212                                 TOUT_CMP='OUI',
 213                                 OPERATION='EXTRACTION',),);
 214
 215 TABDEP2=POST_RELEVE_T(ACTION=_F(INTITULE='DEPB2',
 216                                 GROUP_NO='PB',
 217                                 RESULTAT=RESU2,
 218                                 NOM_CHAM='DEPL',
 219                                 TOUT_CMP='OUI',
 220                                 OPERATION='EXTRACTION',),);
 221
 222 IMPR_TABLE(TABLE=TABDEP1,);
 223
 224 IMPR_TABLE(TABLE=TABDEP2,);
 225
 226 TEST_TABLE(TABLE=TABDEP1,
 227            NOM_PARA='DX',
 228            VALE=-0.02901,
 229            PRECISION=0.0001,
 230            REFERENCE='NON_REGRESSION',
 231            VERSION='6.3.20',);
 232
 233 TEST_TABLE(TABLE=TABDEP1,
 234            NOM_PARA='DY',
 235            VALE=0.1060178,
 236            PRECISION=0.0001,
 237            REFERENCE='NON_REGRESSION',
 238            VERSION='6.3.20',);
 239
 240 TEST_TABLE(TABLE=TABDEP1,
 241            NOM_PARA='DRZ',
 242            VALE=0.032629,
 243            PRECISION=0.0001,
 244            REFERENCE='NON_REGRESSION',
 245            VERSION='6.3.20',);
 246
 247 TEST_TABLE(TABLE=TABDEP2,
 248            NOM_PARA='DX',
 249            VALE=-0.290284,
 250            PRECISION=0.0001,
 251            REFERENCE='NON_REGRESSION',
 252            VERSION='6.3.20',);
 253
 254 TEST_TABLE(TABLE=TABDEP2,
 255            NOM_PARA='DY',
 256            VALE=0.4687,
 257            PRECISION=0.0001,
 258            REFERENCE='NON_REGRESSION',
 259            VERSION='6.3.20',);
 260
 261 TEST_TABLE(TABLE=TABDEP2,
 262            NOM_PARA='DRZ',
 263            VALE=0.13034,
 264            PRECISION=0.0001,
 265            REFERENCE='NON_REGRESSION',
 266            VERSION='6.3.20',);
 267
 268 RESU1=CALC_NO(reuse =RESU1,
 269               RESULTAT=RESU1,
 270               OPTION='FORC_NODA');
 271
 272 MAIL=DEFI_GROUP(reuse =MAIL,
 273                 MAILLAGE=MAIL,
 274                 CREA_GROUP_NO=_F(GROUP_MA='ENCAST',),);
 275
 276 REAC=POST_RELEVE_T(ACTION=_F(INTITULE='REAC',
 277                              GROUP_NO='ENCAST',
 278                              RESULTAT=RESU1,
 279                              NOM_CHAM='FORC_NODA',
 280                              RESULTANTE=('DX','DY','DZ',),
 281                              MOMENT=('DRX','DRY','DRZ',),
 282                              POINT=(0.0,0.0,0.0,),
 283                              OPERATION='EXTRACTION',),);
 284
 285 FAPP=POST_RELEVE_T(ACTION=_F(INTITULE='REAC',
 286                              GROUP_NO='EFOND',
 287                              RESULTAT=RESU1,
 288                              NOM_CHAM='FORC_NODA',
 289                              RESULTANTE=('DX','DY','DZ',),
 290                              MOMENT=('DRX','DRY','DRZ',),
 291                              POINT=(3.6,3.6,0.0,),
 292                              OPERATION='EXTRACTION',),);
 293
 294 IMPR_TABLE(TABLE=REAC,);
 295
 296 IMPR_TABLE(TABLE=FAPP,);
 297
 298 FIN();
 299 #
 300 #
 301 #
 302 #
 303 #
 304 #
 305