Tests/testcomm7/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 FIN()
  90
  91 THETA01=CALC_THETA(       MODELE=MO,
  92                           THETA_2D=_F(  GROUP_NO = ('P0',),
  93                                      MODULE = 1.,
  94                                      R_INF = 7.50E-5,
  95                                      R_SUP = 1.125E-4),
  96                          DIRECTION=(1., 0., 0.,)
  97                      )
  98
  99 #
 100 #  DERIVEE PAR RAPPORT A E
 101 #
 102
 103 DEP2=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT,
 104                       SENSIBILITE = (YOUNG),
 105                               EXCIT=_F( CHARGE = CH),
 106                    )
 107
 108 DEP2=CALC_ELEM(reuse=DEP2,RESULTAT=DEP2,
 109                           SENSIBILITE = (YOUNG),
 110               OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)
 111
 112
 113
 114 FIN()
 115