Aster/az.comm

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