Editeur/Patrons/elasticite.salome.comm

   1 DEBUT();
   2
   3 MA=DEFI_MATERIAU(ELAS=_F(E=210000000000,
   4                          NU=0.3,),);
   5
   6
   7
   8 MAIL=LIRE_MAILLAGE(UNITE=21,
   9                    FORMAT='MED',
  10                    INFO_MED=2,);
  11
  12 MODE=AFFE_MODELE(MAILLAGE=MAIL,
  13                 AFFE=_F(TOUT='OUI',
  14                                 PHENOMENE='MECANIQUE',
  15                                 MODELISATION= '3D',),);
  16
  17
  18
  19 MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
  20                    AFFE=_F(TOUT='OUI',
  21                            MATER=MA,),);
  22
  23 CHAR=AFFE_CHAR_MECA(MODELE=MODE,
  24                         FACE_IMPO=(
  25                                 _F(GROUP_MA='ABB1A1',
  26                                                 DX=0,
  27                                                 DY=0,
  28                                                 DZ=0,),),
  29                         PRES_REP=(
  30                                 _F(GROUP_MA='CDD1C1',
  31                                                 PRES=10000,),),);
  32
  33 RESU=MECA_STATIQUE(MODELE=MODE,
  34                    CHAM_MATER=MATE,
  35                    EXCIT=_F(CHARGE=CHAR,),);
  36
  37 RESU=CALC_ELEM(reuse =RESU,
  38                MODELE=MODE,
  39                CHAM_MATER=MATE,
  40                RESULTAT=RESU,
  41                OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
  42                EXCIT=_F(
  43                CHARGE=CHAR,),);
  44
  45 RESU=CALC_NO(reuse =RESU,
  46              RESULTAT=RESU,
  47              OPTION=('SIGM_NOEU_DEPL', 'EQUI_NOEU_SIGM', ),);
  48
  49 IMPR_RESU(FORMAT='MED',
  50           UNITE=80,
  51           RESU=_F(MAILLAGE=MAIL,
  52                   RESULTAT=RESU,
  53                   NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','DEPL',),),);
  54
  55 FIN();