Editeur/Patrons/statlin_pou.comm

   1
   2 # PATRON
   3 # ANALYSE STATIQUE
   4 # POUTRE
   5
   6 DEBUT();
   7 # MISE EN DONNEES
   8
   9 MAIL=LIRE_MAILLAGE();
  10
  11 MAIL=DEFI_GROUP(reuse =MAIL,
  12                 MAILLAGE=MAIL,
  13                 CREA_GROUP_MA=_F(NOM='TOUT',
  14                                  TOUT='OUI',),);
  15
  16 MODMECA=AFFE_MODELE(MAILLAGE=MAIL,
  17                     AFFE=_F(GROUP_MA='TOUT',
  18                             PHENOMENE='MECANIQUE',
  19                             MODELISATION='POU_D_T',),);
  20
  21 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E10,
  22                             NU=0.3,),);
  23
  24 CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
  25                     AFFE=_F(TOUT='OUI',
  26                             MATER=ACIER,),);
  27 R1 = 1.0;
  28
  29 EP1 = 0.001;
  30
  31
  32 CARA_POU=AFFE_CARA_ELEM(MODELE=MODMECA,
  33                         POUTRE=_F(GROUP_MA='TOUT',
  34                                   SECTION='CERCLE',
  35                                   CARA=('R','EP',),
  36                                   VALE=(R1,EP1,),),);
  37
  38 # CONDITIONS AUX LIMITES ET CHARGEMENT
  39
  40 #POINT A : ENCASTREMENT
  41
  42 BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODMECA,
  43                         DDL_IMPO=_F(GROUP_NO='PA',
  44                                     LIAISON='ENCASTRE',),);
  45 #POINT B : FORCE IMPOSEE
  46
  47 CHARGT=AFFE_CHAR_MECA(MODELE=MODMECA,
  48                       FORCE_NODALE=_F(GROUP_NO='PB',
  49                                       FY=1000.0,),);
  50 # RESOLUTION DU PROBLEME ELASTIQUE
  51
  52 RESU1=MECA_STATIQUE(MODELE=MODMECA,
  53                     CHAM_MATER=CHMAT,
  54                     CARA_ELEM=CARA_POU,
  55                     EXCIT=(_F(CHARGE=BLOQUAGE,),
  56                            _F(CHARGE=CHARGT,),),);
  57 # CALCUL DES EFFORTS
  58
  59 RESU1=CALC_ELEM(reuse =RESU1,
  60                 RESULTAT=RESU1,
  61                 OPTION='EFGE_ELNO_DEPL',);
  62 # RESULTATS FORMAT GMSH
  63
  64 IMPR_RESU(MODELE=MODMECA,
  65           FORMAT='GMSH',
  66           RESU=_F(MAILLAGE=MAIL,
  67                   RESULTAT=RESU1,
  68                   NOM_CHAM=('DEPL','EFGE_ELNO_DEPL',),),);
  69
  70 FIN();