Tests/testelem/azv6.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 P6 = (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(REEL='(a+z) = (a,z)',);
  21
  22 MAIL=LIRE_MAILLAGE(UNITE=P3,);
  23 # 'LIRE_MAILLAGE', 'UNITE'            --> uniquebasepanel
  24
  25 az = FORMULE(REEL='(aaaaa) = (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,P6,),),);
  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                        P4[1],P3,
  61                        ),);
  62
  63 F3=DEFI_FONCTION(NOM_PARA='DRX',VALE_C=(5.0,7.0,9.0,
  64                          9.0,8.0,7.0,
  65                          ),);
  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();