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