Tests/testcomm9/hplp101b.comm

   1 # MODIF  DATE 09/05/2006   AUTEUR REZETTE C.REZETTE
   2 # TITRE FISSURE AU CENTRE D'UNE PLAQUE MINCE RECTANGULAIRE FAISANT
   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 #       OBSTACLE A UN FLUC DE CHALEUR UNIFORME EN MILIEU ISOTROPE
  21 #
  22 # TESTE LES DERIVEES DE G DANS LE CAS D'UN CHARGEMENT EN FORCE DE VOLUME
  23
  24 DEBUT(  CODE=_F( NOM = 'HPLP101B',NIV_PUB_WEB='INTERNET'))
  25
  26 YOUNG = DEFI_PARA_SENSI(VALE=1.)
  27 FORC = DEFI_PARA_SENSI(VALE=1.)
  28
  29 MA=LIRE_MAILLAGE( )
  30
  31 MA=DEFI_GROUP( reuse=MA,   MAILLAGE=MA,
  32                            CREA_GROUP_NO=_F(
  33                          GROUP_MA = ( 'LP9P10',  'LP9P10B',
  34                                       'LP8P9',   'LP8P9B',
  35                                          ))
  36                        )
  37
  38 MO=AFFE_MODELE( MAILLAGE=MA,
  39                              AFFE=_F(  TOUT = 'OUI',
  40                                     PHENOMENE = 'MECANIQUE',
  41                                     MODELISATION = 'C_PLAN')
  42                        )
  43 ZERO =    DEFI_CONSTANTE(VALE=0.)
  44 NU =    DEFI_CONSTANTE(VALE=0.3)
  45 ALPHA =    DEFI_CONSTANTE(VALE=0.)
  46
  47 MAT=DEFI_MATERIAU(ELAS_FO=_F(  E = YOUNG, NU = NU,
  48                                 TEMP_DEF_ALPHA =20.,
  49                                ALPHA = ALPHA, RHO=7800.,))
  50
  51 CHMAT=AFFE_MATERIAU(    MAILLAGE=MA,
  52                                   AFFE=_F(  TOUT = 'OUI',
  53                                          MATER = MAT)
  54                           )
  55
  56 CH=AFFE_CHAR_MECA_F(           MODELE=MO,
  57                        FORCE_INTERNE=_F(TOUT='OUI',FY=FORC,),
  58                                DDL_IMPO=(
  59                                     _F(  GROUP_NO = 'LP8P9',
  60                                                DX = ZERO),
  61                                              _F(  GROUP_NO = 'LP8P9B',
  62                                                DX = ZERO),
  63                                              _F(  GROUP_NO = 'P11',
  64                                                DY = ZERO))
  65                             )
  66
  67 DEP=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT,
  68                       SENSIBILITE = (FORC),
  69                               EXCIT=_F( CHARGE = CH),
  70                      )
  71
  72 DEP=CALC_ELEM(reuse=DEP,RESULTAT=DEP,
  73                           SENSIBILITE = (FORC),
  74               OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)
  75
  76 FOND=DEFI_FOND_FISS(    MAILLAGE=MA,
  77                               FOND_FISS=_F( GROUP_NO = ('P0',)),
  78                            NORMALE=(0., 1., 0.,)
  79                       )
  80
  81 THETA0=CALC_THETA(       MODELE=MO,
  82                          THETA_2D=_F(  GROUP_NO = ('P0',),
  83                                     MODULE = 1.,
  84                                     R_INF = 3.75E-5,
  85                                     R_SUP = 7.50E-5),
  86                         DIRECTION=(1., 0., 0.,)
  87                      )
  88
  89 G0=CALC_G(      RESULTAT=DEP,
  90                         THETA=_F(THETA=THETA0),
  91                         SYME_CHAR='SANS',
  92                         COMP_ELAS=_F( RELATION = 'ELAS',
  93                                       DEFORMATION = 'PETIT')
  94                       )
  95
  96 DG0=CALC_G(     RESULTAT=DEP,
  97                         SENSIBILITE = (FORC),
  98                         THETA=_F(THETA=THETA0),
  99                         SYME_CHAR='SANS',
 100                         COMP_ELAS=_F( RELATION = 'ELAS',
 101                                       DEFORMATION = 'PETIT')
 102                       )
 103
 104 TEST_TABLE(       TABLE=DG0,
 105                        NOM_PARA='DG/DF',
 106                           SENSIBILITE = (FORC),
 107                            VALE=1.06E-2,
 108                       REFERENCE='ANALYTIQUE',
 109                       PRECISION=0.0005    )
 110 FIN()
 111
 112 THETA01=CALC_THETA(       MODELE=MO,
 113                           THETA_2D=_F(  GROUP_NO = ('P0',),
 114                                      MODULE = 1.,
 115                                      R_INF = 7.50E-5,
 116                                      R_SUP = 1.125E-4),
 117                          DIRECTION=(1., 0., 0.,)
 118                      )
 119
 120 G1=CALC_G(        RESULTAT=DEP,
 121                           THETA=_F(THETA=THETA01),
 122                           COMP_ELAS=_F( RELATION = 'ELAS',
 123                                      DEFORMATION = 'PETIT'),
 124                           SYME_CHAR='SANS'
 125                      )
 126
 127 DG1=CALC_G(      RESULTAT=DEP,
 128                          SENSIBILITE = (FORC),
 129                          THETA=_F(THETA=THETA01),
 130                          SYME_CHAR='SANS',
 131                           COMP_ELAS=_F( RELATION = 'ELAS',
 132                                       DEFORMATION = 'PETIT')
 133                       )
 134
 135 TEST_TABLE(       TABLE=DG1,
 136                        NOM_PARA='G',
 137                           SENSIBILITE = (FORC),
 138                       REFERENCE='ANALYTIQUE',
 139                            VALE=1.06E-2,
 140                       PRECISION=0.0005    )
 141
 142 #
 143 #  DERIVEE PAR RAPPORT A E
 144 #
 145
 146 DEP2=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT,
 147                       SENSIBILITE = (YOUNG),
 148                               EXCIT=_F( CHARGE = CH),
 149                    )
 150
 151 DEP2=CALC_ELEM(reuse=DEP2,RESULTAT=DEP2,
 152                           SENSIBILITE = (YOUNG),
 153               OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)
 154
 155
 156 DG2=CALC_G(     RESULTAT=DEP2,
 157                         SENSIBILITE = (YOUNG),
 158                         THETA=_F(THETA=THETA0),
 159                         SYME_CHAR='SANS',
 160                          COMP_ELAS=_F( RELATION = 'ELAS',
 161                                       DEFORMATION = 'PETIT')
 162                       )
 163
 164 TEST_TABLE(       TABLE=DG2,
 165                        NOM_PARA='G',
 166                           SENSIBILITE = (YOUNG),
 167                       REFERENCE='ANALYTIQUE',
 168                            VALE=-5.3E-3,
 169                       PRECISION=0.0005    )
 170
 171
 172 DG3=CALC_G(     RESULTAT=DEP2,
 173                         SENSIBILITE = (YOUNG),
 174                         THETA=_F(THETA=THETA01),
 175                         SYME_CHAR='SANS',
 176                         COMP_ELAS=_F( RELATION = 'ELAS',
 177                                       DEFORMATION = 'PETIT')
 178                       )
 179
 180 TEST_TABLE(       TABLE=DG3,
 181                        NOM_PARA='G',
 182                            VALE=-5.3E-3,
 183                           SENSIBILITE = (YOUNG),
 184                       REFERENCE='ANALYTIQUE',
 185                       PRECISION=0.0005    )
 186
 187
 188 FIN()
 189