doc/files/yacs_aster_01.comm

   1 # -*- coding: iso-8859-1 -*-
   2
   3 DEBUT();
   4
   5 BOX=LIRE_MAILLAGE(FORMAT='MED',);
   6
   7 BETON=DEFI_MATERIAU(THER=_F(LAMBDA=8000.,
   8                             RHO_CP=2.400E+6,),);
   9
  10 MOTH=AFFE_MODELE(MAILLAGE=BOX,
  11                  AFFE=_F(TOUT='OUI',
  12                          PHENOMENE='THERMIQUE',
  13                          MODELISATION='3D',),);
  14
  15 #BOX=MODI_MAILLAGE( reuse=BOX,
  16 #                  MAILLAGE=BOX,
  17 #                  ORIE_PEAU_3D=_F(GROUP_MA=('CL_1','CL_2','XMIN','XMAX','YMIN','YMAX','ZMIN','ZMAX',),)
  18 #                )
  19
  20 CHMAT=AFFE_MATERIAU(MAILLAGE=BOX,
  21                     AFFE=_F(TOUT='OUI',
  22                             MATER=BETON,),);
  23
  24 CHTH=AFFE_CHAR_THER(MODELE=MOTH,
  25                     TEMP_IMPO=(_F(GROUP_MA='CL_1',
  26                                   TEMP=50.0,),
  27                                _F(GROUP_MA='CL_2',
  28                                   TEMP=5.0,),),
  29                     FLUX_REP=(_F(GROUP_MA='XMIN',
  30                                  FLUN=0.0,),
  31                               _F(GROUP_MA='XMAX',
  32                                  FLUN=0.0,),
  33                               _F(GROUP_MA='YMIN',
  34                                  FLUN=0.0,),
  35                               _F(GROUP_MA='YMAX',
  36                                  FLUN=0.0,),
  37                               _F(GROUP_MA='ZMIN',
  38                                  FLUN=0.0,),
  39                               _F(GROUP_MA='ZMAX',
  40                                  FLUN=0.0,),),);
  41
  42 RESU=THER_LINEAIRE(MODELE=MOTH,
  43                    CHAM_MATER=CHMAT,
  44                    EXCIT=_F(CHARGE=CHTH,),);
  45
  46 IMPR_RESU(FORMAT='MED',
  47           RESU=_F(RESULTAT=RESU,
  48                   NOM_CHAM='TEMP',
  49                   NOM_CHAM_MED='TEMPERATURE',),);
  50
  51 # Extraction de la temperature sur un noeud predefini
  52
  53 TA=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
  54                            INTITULE='V_TEST',
  55                            RESULTAT=RESU,
  56                            NOM_CHAM='TEMP',
  57                            NUME_ORDRE=0,
  58                            GROUP_NO='A',
  59                            NOM_CMP='TEMP',),);
  60
  61 #
  62
  63 IMPR_TABLE(TABLE=TA,
  64            NOM_PARA=('INTITULE','TEMP',),);
  65
  66 # Calcul de l'indicateur d'erreur
  67
  68 RESU=CALC_ERREUR(reuse =RESU,
  69                  RESULTAT=RESU,
  70                  OPTION='ERTH_ELEM',);
  71
  72
  73 IMPR_RESU(FORMAT='MED',
  74           RESU=_F(RESULTAT=RESU,
  75                   NUME_ORDRE=0,
  76                   NOM_CHAM='ERTH_ELEM',
  77                   NOM_CMP='ERTABS',
  78                   NOM_CHAM_MED='ERREUR',),);
  79
  80
  81 FIN();
  82
  83 #
  84
  85 #CHECKSUM:5021673eed6e81bf9074892cd639b2ad  -:FIN CHECKSUM