Editeur/Patrons/statlin_3D.comm

   1
   2 # PATRON
   3 # ANALYSE STATIQUE
   4 # 3D
   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='3D',),);
  20
  21 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E11,
  22                             NU=0.3,),);
  23
  24 CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
  25                     AFFE=_F(TOUT='OUI',
  26                             MATER=ACIER,),);
  27 # CONDITIONS AUX LIMITES ET CHARGEMENT
  28 #POINT A : ENCASTREMENT
  29
  30 BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODMECA,
  31                         DDL_IMPO=_F(GROUP_MA='ENCAST',
  32                                     LIAISON='ENCASTRE',),);
  33 #POINT B : FORCE IMPOSEE
  34
  35 CHARGT=AFFE_CHAR_MECA(MODELE=MODMECA,
  36                       FORCE_NODALE=_F(GROUP_NO='PB',
  37                                       FY=1000.0,),);
  38 # RESOLUTION DU PROBLEME ELASTIQUE
  39
  40 RESU1=MECA_STATIQUE(MODELE=MODMECA,
  41                     CHAM_MATER=CHMAT,
  42                     EXCIT=(_F(CHARGE=BLOQUAGE,),
  43                            _F(CHARGE=CHARGT,),),);
  44 # CALCUL DES EFFORTS
  45
  46 RESU1=CALC_ELEM(reuse =RESU1,
  47                 RESULTAT=RESU1,
  48                 OPTION='EFGE_ELNO_DEPL',);
  49 # RESULTATS FORMAT GMSH
  50
  51 IMPR_RESU(MODELE=MODMECA,
  52           FORMAT='GMSH',
  53           RESU=_F(MAILLAGE=MAIL,
  54                   RESULTAT=RESU1,
  55                   NOM_CHAM=('DEPL','EFGE_ELNO_DEPL',),),);
  56
  57 FIN();