Aster/az.comm

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