Aster/az.comm

   1
   2
   3 DEBUT();
   4 P1 = 9.8;
   5
   6 P2 = 8.8;
   7
   8 P3 = 7;
   9
  10 P5 = P3*P1;
  11
  12 P4 = [2,3,4];
  13
  14 #Pas trouve             shellpanel
  15
  16 aaa = FORMULE(VALE='a+z',
  17               NOM_PARA=('a','z',),);
  18
  19 MAIL=LIRE_MAILLAGE(UNITE=P3,);
  20
  21 MAILLA2=LIRE_MAILLAGE(UNITE=P4[1],);
  22 # 'LIRE_MAILLAGE', 'UNITE'            --> uniquebasepanel
  23
  24 AFFE1=AFFE_MODELE(MAILLAGE=MAILLA2,
  25                   AFFE=(_F(GROUP_MA=('RESSORT','eee',),
  26                            PHENOMENE='MECANIQUE',
  27                            MODELISATION='DIS_T',),
  28                         _F(GROUP_MA='MASSES',
  29                            PHENOMENE='MECANIQUE',
  30                            MODELISATION='DIS_T',),
  31                         _F(GROUP_NO=('GNP3','GNP5','GNP6','GNP7','GNP8','GNP9','GNP10','GNP11','GNP12',),
  32                            PHENOMENE='ACOUSTIQUE',
  33                            MODELISATION='PLAN',),),);
  34
  35 MOD=AFFE_MODELE(MAILLAGE=MAIL,
  36                 AFFE=(_F(GROUP_MA='RESSORT',
  37                          PHENOMENE='MECANIQUE',
  38                          MODELISATION='DIS_T',),
  39                       _F(GROUP_MA='MASSES',
  40                          PHENOMENE='MECANIQUE',
  41                          MODELISATION='DIS_T',),
  42                       _F(TOUT='OUI',
  43                          PHENOMENE='THERMIQUE',
  44                          MODELISATION='COQUE',),),);
  45
  46 CARA=AFFE_CARA_ELEM(MODELE=MOD,
  47                     POUTRE=_F(GROUP_MA='MA',
  48                               SECTION='CERCLE',
  49                               CARA='R',
  50                               VALE=1.0,),);
  51 # 'AFFE_MODELE', 'MAILLAGE'           --> uniqueassdpanel
  52 #  AFFE_MODELE', 'AFFE', 'GROUP_MA'   --> plusieursbasepanel
  53 # 'AFFE_MODELE', 'AFFE', 'PHENOMENE'  --> uniqueintopanel
  54 # 'AFFE_MODELE', 'AFFE', 'b_mecanique'--> plusieursintopanel
  55
  56 F1=DEFI_FONCTION(NOM_PARA='DX',
  57                  VALE=(5.0,3.0,P4[1],P3,),);
  58
  59 F3=DEFI_FONCTION(NOM_PARA='DRX',
  60                  VALE_C=(5.0,7.0,9.0,9.0,8.0,7.0,),);
  61 # 'DEFI_FONCTION', 'VALE'             --> fonctionpanel
  62
  63 MATER2=DEFI_MATERIAU(ELAS=_F(E=100000000000.0,
  64                              NU=0.0,),
  65                      ECRO_ASYM_LINE=_F(DC_SIGM_EPSI=0.0,
  66                                        SY_C=200000000.0,
  67                                        DT_SIGM_EPSI=0.0,
  68                                        SY_T=50000000.0,),);
  69
  70 PS1=DEFI_PARA_SENSI(VALE=1.0,);
  71
  72 PS2=DEFI_PARA_SENSI(VALE=1.0,);
  73
  74 PS3=DEFI_PARA_SENSI(VALE=1.0,);
  75
  76 CHMAT2=AFFE_MATERIAU(MAILLAGE=MAIL,
  77                      AFFE=_F(TOUT='OUI',
  78                              MATER=MATER2,),);
  79
  80 AAAZ=AFFE_CHAR_THER(MODELE=AFFE1,
  81                     TEMP_IMPO=_F(TOUT='OUI',
  82                                  TEMP=0.0,),);
  83
  84 TH1=THER_LINEAIRE(MODELE=AFFE1,
  85                   CHAM_MATER=CHMAT2,
  86                   EXCIT=_F(CHARGE=AAAZ,),
  87                   SENSIBILITE=(PS1,PS2,),);
  88 # 'THER_LINEAIRE', 'SENSIBILITE'       --> plusieursassdpanel
  89
  90 ACA1=AFFE_CHAR_ACOU(MODELE=AFFE1,
  91                     PRES_IMPO=_F(TOUT='OUI',
  92                                  PRES=('RI',3.0,3.0,),),);
  93 # 'AFFE_CHAR_ACOU', 'PRES_IMPO', 'PRES' --> uniquecomppanel
  94
  95 MACRO_MATR_ASSE(MODELE=AFFE1,
  96                 NUME_DDL=CO('DDL1'),
  97                 MATR_ASSE=_F(MATRICE=CO('MAT1'),
  98                              OPTION='RIGI_THER',),);
  99 # 'MACRO_MATR_ASSE', 'MATR_ASSE', 'MATRICE'  --> uniquesdcopanel
 100
 101 FIN();