Tests/testcomm9/titi.comm

   1 DEBUT()
   2 YOUNG = DEFI_PARA_SENSI(VALE=1.)
   3 FORC = DEFI_PARA_SENSI(VALE=1.)
   4 MA=LIRE_MAILLAGE()
   5 MA=DEFI_GROUP( reuse=MA,   MAILLAGE=MA,
   6                            CREA_GROUP_NO=_F(
   7                          GROUP_MA = ( 'LP9P10',  'LP9P10B',
   8                                       'LP8P9',   'LP8P9B',
   9                                          ))
  10                        )
  11
  12 MO=AFFE_MODELE( MAILLAGE=MA,
  13                              AFFE=_F(  TOUT = 'OUI',
  14                                     PHENOMENE = 'MECANIQUE',
  15                                     MODELISATION = 'C_PLAN')
  16                        )
  17 ZERO =    DEFI_CONSTANTE(VALE=0.)
  18 NU =    DEFI_CONSTANTE(VALE=0.3)
  19 ALPHA =    DEFI_CONSTANTE(VALE=0.)
  20
  21 MAT=DEFI_MATERIAU(ELAS_FO=_F(  E = YOUNG, NU = NU,
  22                                 TEMP_DEF_ALPHA =20.,
  23                                ALPHA = ALPHA, RHO=7800.,))
  24
  25 CHMAT=AFFE_MATERIAU(    MAILLAGE=MA,
  26                                   AFFE=_F(  TOUT = 'OUI',
  27                                          MATER = MAT)
  28                           )
  29
  30 CH=AFFE_CHAR_MECA_F(           MODELE=MO,
  31                        FORCE_INTERNE=_F(TOUT='OUI',FY=FORC,),
  32                                DDL_IMPO=(
  33                                     _F(  GROUP_NO = 'LP8P9',
  34                                                DX = ZERO),
  35                                              _F(  GROUP_NO = 'LP8P9B',
  36                                                DX = ZERO),
  37                                              _F(  GROUP_NO = 'P11',
  38                                                DY = ZERO))
  39                             )
  40
  41 DEP=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT,
  42                       SENSIBILITE = (FORC),
  43                               EXCIT=_F( CHARGE = CH),
  44                      )
  45
  46 DEP=CALC_ELEM(reuse=DEP,RESULTAT=DEP,
  47                           SENSIBILITE = (FORC),
  48               OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)
  49
  50 FOND=DEFI_FOND_FISS(    MAILLAGE=MA,
  51                               FOND_FISS=_F( GROUP_NO = ('P0',)),
  52                            NORMALE=(0., 1., 0.,)
  53                       )
  54
  55 THETA0=CALC_THETA(       MODELE=MO,
  56                          THETA_2D=_F(  GROUP_NO = ('P0',),
  57                                     MODULE = 1.,
  58                                     R_INF = 3.75E-5,
  59                                     R_SUP = 7.50E-5),
  60                         DIRECTION=(1., 0., 0.,)
  61                      )
  62
  63
  64 FIN()
  65 MA1=LIRE_MAILLAGE()
  66 MO1=AFFE_MODELE( MAILLAGE=MA1, AFFE=_F(  TOUT = 'OUI',
  67                  PHENOMENE = 'MECANIQUE', MODELISATION = 'C_PLAN'))
  68
  69 FIN()