Tests/testcomm7/e.comm

   1 DEBUT()
   2 MAYA=LIRE_MAILLAGE()
   3
   4 MAYA=DEFI_GROUP( reuse=MAYA,   MAILLAGE=MAYA,
   5   CREA_GROUP_MA=_F(  NOM = 'TOUT', TOUT = 'OUI'))
   6
   7 BARRE1=AFFE_MODELE(  MAILLAGE=MAYA,
   8                             AFFE=_F( GROUP_MA='SS1',
   9                                   MODELISATION = 'POU_D_E',
  10                                   PHENOMENE = 'MECANIQUE'))
  11 #
  12 MATERIO1=DEFI_MATERIAU(  ELAS=_F( RHO = 1.E08,  NU = 0.3,  E = 1.E10))
  13
  14 MATERIO2=DEFI_MATERIAU(  ELAS=_F( RHO = 0.,  NU = 0.,   E = 1.E15))
  15 #
  16 CHMAT=AFFE_MATERIAU(  MAILLAGE=MAYA,AFFE=(
  17                               #_F( TOUT = 'OUI',      MATER = MATERIO1),
  18                               _F( GROUP_MA = 'SS1',   MATER = MATERIO1),
  19                               _F( GROUP_MA = 'SS2',   MATER = MATERIO1),
  20                               _F( GROUP_MA = 'SS3',   MATER = MATERIO1),
  21                               _F( GROUP_MA = 'RIGIDE', MATER = MATERIO2),
  22 )) # commentaire test
  23
  24 CARA1=AFFE_CARA_ELEM(  MODELE=BARRE1,
  25                            POUTRE=_F(
  26                                     GROUP_MA = 'TOUT',
  27                                     SECTION = 'CERCLE',
  28                                     CARA = ('R', 'EP',),
  29                                     VALE = (0.1, 0.01,)))
  30
  31 GUIDAGE1=AFFE_CHAR_MECA(    MODELE=BARRE1,DDL_IMPO=(
  32                            _F( GROUP_MA='SS1',
  33                                      DX = 0., DZ = 0., DRX = 0., DRY = 0.),
  34                            _F( GROUP_NO = 'L1',  DY = 0., DRZ = 0.),
  35                            ))
  36
  37 K_ELEM1=CALC_MATR_ELEM(      MODELE=BARRE1,
  38                              CARA_ELEM=CARA1,
  39                             CHAM_MATER=CHMAT,
  40                                OPTION='RIGI_MECA',
  41                                CHARGE=GUIDAGE1)
  42
  43 NUM1=NUME_DDL(  MATR_RIGI=K_ELEM1)
  44
  45 #
  46 K_ASSE1=ASSE_MATRICE(  MATR_ELEM=K_ELEM1,
  47                            NUME_DDL=NUM1)
  48
  49 MODESTA1=MODE_STATIQUE(  MATR_RIGI=K_ASSE1,
  50                      MODE_STAT=_F(
  51                      GROUP_NO=('L12'),
  52                      AVEC_CMP = ('DY','DRZ'),
  53                                  ))
  54
  55 # parse: -affectation
  56
  57 DS1=[None]*5
  58 DS2=[None]*5
  59 DS3=[None]*5
  60 DS4=[None]*5
  61 CHS1=[None]*5
  62 CHS2=[None]*5
  63
  64 # parse: +affectation
  65
  66 #for k in range(1,5):
  67 #  DS1[k] = CREA_CHAMP( OPERATION='EXTR', TYPE_CHAM='NOEU_DEPL_R',
  68 #                  RESULTAT= MODESTA1, NUME_ORDRE=k, NOM_CHAM = 'DEPL');
  69
  70 FIN()