Exemples/ex23/titi.comm

   1
   2
   3 DEBUT(CODE=_F(NOM='AHLV100A',
   4               NIV_PUB_WEB=None,),);
   5
   6 MAIL=LIRE_MAILLAGE();
   7
   8 AIR=DEFI_MATERIAU(FLUIDE=_F(CELE_C=('RI',343.0,0.0,),
   9                             RHO=1.3,),);
  10
  11 CHAMPMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
  12                        AFFE=_F(MATER=AIR,
  13                                TOUT='OUI',),);
  14
  15 GUIDE=AFFE_MODELE(MAILLAGE=MAIL,
  16                   VERIF='MAILLE',
  17                   AFFE=_F(PHENOMENE='ACOUSTIQUE',
  18                           TOUT='OUI',
  19                           MODELISATION='3D',),);
  20
  21 CHARACOU=AFFE_CHAR_ACOU(MODELE=GUIDE,
  22                         VITE_FACE=_F(GROUP_MA='ENTREE',
  23                                      VNOR=('RI',0.014,0.0,),),);
  24
  25 IMPEACOU=AFFE_CHAR_ACOU(MODELE=GUIDE,
  26                         IMPE_FACE=_F(GROUP_MA='SORTIE',
  27                                      IMPE=('RI',445.9,0.0,),),);
  28
  29 MACRO_MATR_ASSE(CHARGE=IMPEACOU,
  30                 MODELE=GUIDE,
  31                 NUME_DDL=CO('NUM'),
  32                 CHAM_MATER=CHAMPMAT,
  33                 MATR_ASSE=(_F(OPTION='RIGI_ACOU',
  34                               MATRICE=CO('MATASK'),),
  35                            _F(OPTION='MASS_ACOU',
  36                               MATRICE=CO('MATASM'),),
  37                            _F(OPTION='AMOR_ACOU',
  38                               MATRICE=CO('MATASI'),),),);
  39
  40 VECTELEM=CALC_VECT_ELEM(OPTION='CHAR_ACOU',
  41                         CHARGE=CHARACOU,
  42                         CHAM_MATER=CHAMPMAT,);
  43
  44 IMPR_MATRICE(MATR_ELEM=_F(FICHIER='RESULTAT',
  45                           FORMAT='RESULTAT',
  46                           MATRICE=VECTELEM,
  47                           GRAIN='MAILLE',),);
  48
  49 VECTASS=ASSE_VECTEUR(NUME_DDL=NUM,
  50                      VECT_ELEM=VECTELEM,);
  51
  52 MATASKR=COMB_MATR_ASSE(COMB_R=_F(MATR_ASSE=MATASK,
  53                                  COEF_R=1.0,
  54                                  PARTIE='REEL',),);
  55
  56 MATASMR=COMB_MATR_ASSE(COMB_R=_F(MATR_ASSE=MATASM,
  57                                  COEF_R=1.0,
  58                                  PARTIE='REEL',),);
  59
  60 MODES=MODE_ITER_SIMULT(MATR_B=MATASMR,
  61                        MATR_A=MATASKR,
  62                        CALC_FREQ=_F(OPTION='BANDE',
  63                                     FREQ=(1.0,1000.0,),),);
  64
  65 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
  66                    PARA='FREQ',
  67                    RESULTAT=MODES,
  68                    VALE=171.5,
  69                    CRITERE='RELATIF',
  70                    PRECISION=0.0001,),
  71                 _F(NUME_ORDRE=2,
  72                    PARA='FREQ',
  73                    RESULTAT=MODES,
  74                    VALE=343.0,
  75                    CRITERE='RELATIF',
  76                    PRECISION=0.0001,),
  77                 _F(NUME_ORDRE=3,
  78                    PARA='FREQ',
  79                    RESULTAT=MODES,
  80                    VALE=514.5,
  81                    CRITERE='RELATIF',
  82                    PRECISION=0.001,),
  83                 _F(NUME_ORDRE=4,
  84                    PARA='FREQ',
  85                    RESULTAT=MODES,
  86                    VALE=686.0,
  87                    CRITERE='RELATIF',
  88                    PRECISION=0.01,),
  89                 _F(NUME_ORDRE=5,
  90                    PARA='FREQ',
  91                    RESULTAT=MODES,
  92                    VALE=857.5,
  93                    CRITERE='RELATIF',
  94                    PRECISION=0.001,),
  95                 _F(NUME_ORDRE=6,
  96                    PARA='FREQ',
  97                    RESULTAT=MODES,
  98                    VALE=857.5,
  99                    CRITERE='RELATIF',
 100                    PRECISION=0.005,),
 101                 _F(NUME_ORDRE=7,
 102                    PARA='FREQ',
 103                    RESULTAT=MODES,
 104                    VALE=874.482,
 105                    CRITERE='RELATIF',
 106                    PRECISION=0.005,),
 107                 _F(NUME_ORDRE=8,
 108                    PARA='FREQ',
 109                    RESULTAT=MODES,
 110                    VALE=923.556,
 111                    CRITERE='RELATIF',
 112                    PRECISION=0.005,),),);
 113
 114 MATRES=COMB_MATR_ASSE(COMB_C=(_F(MATR_ASSE=MATASK,
 115                                  COEF_R=1.0,),
 116                               _F(MATR_ASSE=MATASM,
 117                                  COEF_R=-9869604.4,),
 118                               _F(MATR_ASSE=MATASI,
 119                                  COEF_C=('RI',0.0,3141.5926,),),),);
 120
 121 IMPR_MATRICE(MATR_ASSE=_F(FICHIER='RESULTAT',
 122                           OPTION='LIGNE',
 123                           FORMAT='RESULTAT',
 124                           MATRICE=MATRES,
 125                           GRAIN='NOEUD',),);
 126
 127 IMPR_MATRICE(MATR_ASSE=_F(FICHIER='RESULTAT',
 128                           OPTION='LIGNE',
 129                           FORMAT='RESULTAT',
 130                           MATRICE=MATRES,
 131                           GRAIN='VALEUR',),);
 132
 133 VECRES=COMB_CHAM_NO(COMB_C=_F(CHAM_NO=VECTASS,
 134                               COEF_C=('RI',0.0,-3141.5926,),),);
 135
 136 MATRES=FACT_LDLT(reuse =MATRES,
 137                  MATR_ASSE=MATRES,);
 138
 139 VECRES=RESO_LDLT(reuse =VECRES,
 140                  CHAM_NO=VECRES,
 141                  MATR_FACT=MATRES,);
 142
 143 TEST_RESU(CHAM_NO=(_F(NOEUD='NO1',
 144                       VALE_C=('RI',-6.2426,0.0,),
 145                       NOM_CMP='PRES',
 146                       CHAM_GD=VECRES,),
 147                    _F(NOEUD='NO780',
 148                       VALE_C=('RI',-6.2426,0.0,),
 149                       NOM_CMP='PRES',
 150                       CHAM_GD=VECRES,),
 151                    _F(NOEUD='NO751',
 152                       VALE_C=('RI',6.02368,1.6387,),
 153                       NOM_CMP='PRES',
 154                       CHAM_GD=VECRES,),
 155                    _F(NOEUD='NO763',
 156                       VALE_C=('RI',6.02368,1.6387,),
 157                       NOM_CMP='PRES',
 158                       CHAM_GD=VECRES,),),);
 159
 160 CHNODBEL=CALC_CHAM_ELEM(OPTION='PRES_ELNO_DBEL',
 161                         CHAM_MATER=CHAMPMAT,
 162                         MODELE=GUIDE,
 163                         PRES=VECRES,);
 164
 165 TEST_RESU(CHAM_ELEM=(_F(NOEUD='NO1',
 166                         NOM_CMP='DB',
 167                         VALE=109.9,
 168                         MAILLE='MA1',
 169                         CHAM_GD=CHNODBEL,),
 170                      _F(NOEUD='NO780',
 171                         NOM_CMP='DB',
 172                         VALE=109.9,
 173                         MAILLE='MA1',
 174                         CHAM_GD=CHNODBEL,),
 175                      _F(NOEUD='NO751',
 176                         NOM_CMP='DB',
 177                         VALE=109.9,
 178                         MAILLE='MA57',
 179                         CHAM_GD=CHNODBEL,),
 180                      _F(NOEUD='NO763',
 181                         NOM_CMP='DB',
 182                         VALE=109.9,
 183                         MAILLE='MA57',
 184                         CHAM_GD=CHNODBEL,),),);
 185
 186 PREEL=CALC_CHAM_ELEM(OPTION='PRES_ELNO_REEL',
 187                      CHAM_MATER=CHAMPMAT,
 188                      MODELE=GUIDE,
 189                      PRES=VECRES,);
 190
 191 TEST_RESU(CHAM_ELEM=(_F(NOEUD='NO1',
 192                         NOM_CMP='PRES',
 193                         VALE=-6.2426,
 194                         MAILLE='MA1',
 195                         CHAM_GD=PREEL,),
 196                      _F(NOEUD='NO780',
 197                         NOM_CMP='PRES',
 198                         VALE=-6.2426,
 199                         MAILLE='MA1',
 200                         CHAM_GD=PREEL,),
 201                      _F(NOEUD='NO751',
 202                         NOM_CMP='PRES',
 203                         VALE=6.02368,
 204                         MAILLE='MA57',
 205                         CHAM_GD=PREEL,),
 206                      _F(NOEUD='NO763',
 207                         NOM_CMP='PRES',
 208                         VALE=6.02368,
 209                         MAILLE='MA57',
 210                         CHAM_GD=PREEL,),),);
 211
 212 PIMAG=CALC_CHAM_ELEM(OPTION='PRES_ELNO_IMAG',
 213                      CHAM_MATER=CHAMPMAT,
 214                      MODELE=GUIDE,
 215                      PRES=VECRES,);
 216
 217 TEST_RESU(CHAM_ELEM=(_F(NOEUD='NO1',
 218                         CRITERE='ABSOLU',
 219                         NOM_CMP='PRES',
 220                         VALE=0.0,
 221                         MAILLE='MA1',
 222                         CHAM_GD=PIMAG,),
 223                      _F(NOEUD='NO780',
 224                         CRITERE='ABSOLU',
 225                         NOM_CMP='PRES',
 226                         VALE=0.0,
 227                         MAILLE='MA1',
 228                         CHAM_GD=PIMAG,),
 229                      _F(NOEUD='NO751',
 230                         NOM_CMP='PRES',
 231                         PRECISION=0.0033,
 232                         VALE=1.6387,
 233                         MAILLE='MA57',
 234                         CHAM_GD=PIMAG,),
 235                      _F(NOEUD='NO763',
 236                         NOM_CMP='PRES',
 237                         PRECISION=0.0033,
 238                         VALE=1.6387,
 239                         MAILLE='MA57',
 240                         CHAM_GD=PIMAG,),),);
 241
 242 PRESHARM=DYNA_LINE_HARM(MATR_RIGI=MATASK,
 243                         MATR_MASS=MATASM,
 244                         MATR_AMOR=MATASI,
 245                         FREQ=500.0,
 246                         EXCIT=_F(VECT_ASSE=VECTASS,
 247                                  PHAS_DEG=90.0,
 248                                  COEF_MULT=-1.0,
 249                                  PUIS_PULS=1,),);
 250
 251 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
 252                    VALE_C=('RI',-6.2426,0.0,),
 253                    RESULTAT=PRESHARM,
 254                    NOM_CHAM='PRES',
 255                    NOEUD='NO1',
 256                    NOM_CMP='PRES',),
 257                 _F(NUME_ORDRE=1,
 258                    VALE_C=('RI',-6.2426,0.0,),
 259                    RESULTAT=PRESHARM,
 260                    NOM_CHAM='PRES',
 261                    NOEUD='NO780',
 262                    NOM_CMP='PRES',),
 263                 _F(NUME_ORDRE=1,
 264                    VALE_C=('RI',6.02368,1.6387,),
 265                    RESULTAT=PRESHARM,
 266                    NOM_CHAM='PRES',
 267                    NOEUD='NO751',
 268                    NOM_CMP='PRES',),
 269                 _F(NUME_ORDRE=1,
 270                    VALE_C=('RI',6.02368,1.6387,),
 271                    RESULTAT=PRESHARM,
 272                    NOM_CHAM='PRES',
 273                    NOEUD='NO763',
 274                    NOM_CMP='PRES',),),);
 275
 276 PRESHARM=CALC_ELEM(reuse =PRESHARM,
 277                    OPTION='INTE_ELNO_ACTI',
 278                    CHAM_MATER=CHAMPMAT,
 279                    RESULTAT=PRESHARM,
 280                    TOUT_ORDRE='OUI',
 281                    MODELE=GUIDE,);
 282
 283 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
 284                    RESULTAT=PRESHARM,
 285                    NOM_CHAM='INTE_ELNO_ACTI',
 286                    NOEUD='NO1',
 287                    NOM_CMP='INTX',
 288                    VALE=0.0437,
 289                    PRECISION=0.03,
 290                    MAILLE='MA1',),
 291                 _F(NUME_ORDRE=1,
 292                    RESULTAT=PRESHARM,
 293                    NOM_CHAM='INTE_ELNO_ACTI',
 294                    NOEUD='NO780',
 295                    NOM_CMP='INTX',
 296                    VALE=0.0437,
 297                    PRECISION=0.03,
 298                    MAILLE='MA1',),
 299                 _F(NUME_ORDRE=1,
 300                    RESULTAT=PRESHARM,
 301                    NOM_CHAM='INTE_ELNO_ACTI',
 302                    NOEUD='NO751',
 303                    NOM_CMP='INTX',
 304                    VALE=0.0437,
 305                    PRECISION=0.03,
 306                    MAILLE='MA57',),
 307                 _F(NUME_ORDRE=1,
 308                    RESULTAT=PRESHARM,
 309                    NOM_CHAM='INTE_ELNO_ACTI',
 310                    NOEUD='NO763',
 311                    NOM_CMP='INTX',
 312                    VALE=0.0437,
 313                    PRECISION=0.03,
 314                    MAILLE='MA57',),),);
 315
 316 PRESHARM=CALC_ELEM(reuse =PRESHARM,
 317                    OPTION='INTE_ELNO_REAC',
 318                    CHAM_MATER=CHAMPMAT,
 319                    RESULTAT=PRESHARM,
 320                    TOUT_ORDRE='OUI',
 321                    MODELE=GUIDE,);
 322
 323 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
 324                    RESULTAT=PRESHARM,
 325                    NOM_CHAM='INTE_ELNO_REAC',
 326                    NOEUD='NO1',
 327                    NOM_CMP='INTX',
 328                    VALE=0.0,
 329                    CRITERE='ABSOLU',
 330                    PRECISION=0.00035,
 331                    MAILLE='MA1',),
 332                 _F(NUME_ORDRE=1,
 333                    RESULTAT=PRESHARM,
 334                    NOM_CHAM='INTE_ELNO_REAC',
 335                    NOEUD='NO780',
 336                    NOM_CMP='INTX',
 337                    VALE=0.0,
 338                    CRITERE='ABSOLU',
 339                    PRECISION=0.00035,
 340                    MAILLE='MA1',),
 341                 _F(NUME_ORDRE=1,
 342                    RESULTAT=PRESHARM,
 343                    NOM_CHAM='INTE_ELNO_REAC',
 344                    NOEUD='NO751',
 345                    NOM_CMP='INTX',
 346                    VALE=0.0,
 347                    CRITERE='ABSOLU',
 348                    PRECISION=0.00035,
 349                    MAILLE='MA57',),
 350                 _F(NUME_ORDRE=1,
 351                    RESULTAT=PRESHARM,
 352                    NOM_CHAM='INTE_ELNO_REAC',
 353                    NOEUD='NO763',
 354                    NOM_CMP='INTX',
 355                    VALE=0.0,
 356                    CRITERE='ABSOLU',
 357                    PRECISION=0.00035,
 358                    MAILLE='MA57',),),);
 359
 360 PRESHARM=CALC_ELEM(reuse =PRESHARM,
 361                    OPTION='PRES_ELNO_DBEL',
 362                    CHAM_MATER=CHAMPMAT,
 363                    RESULTAT=PRESHARM,
 364                    TOUT_ORDRE='OUI',
 365                    MODELE=GUIDE,);
 366
 367 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
 368                    RESULTAT=PRESHARM,
 369                    NOM_CHAM='PRES_ELNO_DBEL',
 370                    NOEUD='NO1',
 371                    NOM_CMP='DB',
 372                    VALE=109.9,
 373                    MAILLE='MA1',),
 374                 _F(NUME_ORDRE=1,
 375                    RESULTAT=PRESHARM,
 376                    NOM_CHAM='PRES_ELNO_DBEL',
 377                    NOEUD='NO780',
 378                    NOM_CMP='DB',
 379                    VALE=109.9,
 380                    MAILLE='MA1',),
 381                 _F(NUME_ORDRE=1,
 382                    RESULTAT=PRESHARM,
 383                    NOM_CHAM='PRES_ELNO_DBEL',
 384                    NOEUD='NO751',
 385                    NOM_CMP='DB',
 386                    VALE=109.9,
 387                    MAILLE='MA57',),
 388                 _F(NUME_ORDRE=1,
 389                    RESULTAT=PRESHARM,
 390                    NOM_CHAM='PRES_ELNO_DBEL',
 391                    NOEUD='NO763',
 392                    NOM_CMP='DB',
 393                    VALE=109.9,
 394                    MAILLE='MA57',),),);
 395
 396 PRESHARM=CALC_ELEM(reuse =PRESHARM,
 397                    OPTION='PRES_ELNO_REEL',
 398                    CHAM_MATER=CHAMPMAT,
 399                    RESULTAT=PRESHARM,
 400                    TOUT_ORDRE='OUI',
 401                    MODELE=GUIDE,);
 402
 403 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
 404                    RESULTAT=PRESHARM,
 405                    NOM_CHAM='PRES_ELNO_REEL',
 406                    NOEUD='NO1',
 407                    NOM_CMP='PRES',
 408                    VALE=-6.2426,
 409                    MAILLE='MA1',),
 410                 _F(NUME_ORDRE=1,
 411                    RESULTAT=PRESHARM,
 412                    NOM_CHAM='PRES_ELNO_REEL',
 413                    NOEUD='NO780',
 414                    NOM_CMP='PRES',
 415                    VALE=-6.2426,
 416                    MAILLE='MA1',),
 417                 _F(NUME_ORDRE=1,
 418                    RESULTAT=PRESHARM,
 419                    NOM_CHAM='PRES_ELNO_REEL',
 420                    NOEUD='NO751',
 421                    NOM_CMP='PRES',
 422                    VALE=6.02368,
 423                    MAILLE='MA57',),
 424                 _F(NUME_ORDRE=1,
 425                    RESULTAT=PRESHARM,
 426                    NOM_CHAM='PRES_ELNO_REEL',
 427                    NOEUD='NO763',
 428                    NOM_CMP='PRES',
 429                    VALE=6.02368,
 430                    MAILLE='MA57',),),);
 431
 432 PRESHARM=CALC_ELEM(reuse =PRESHARM,
 433                    OPTION='PRES_ELNO_IMAG',
 434                    CHAM_MATER=CHAMPMAT,
 435                    RESULTAT=PRESHARM,
 436                    TOUT_ORDRE='OUI',
 437                    MODELE=GUIDE,);
 438
 439 TEST_RESU(RESU=(_F(NUME_ORDRE=1,
 440                    RESULTAT=PRESHARM,
 441                    NOM_CHAM='PRES_ELNO_IMAG',
 442                    NOEUD='NO1',
 443                    NOM_CMP='PRES',
 444                    VALE=0.0,
 445                    CRITERE='ABSOLU',
 446                    MAILLE='MA1',),
 447                 _F(NUME_ORDRE=1,
 448                    RESULTAT=PRESHARM,
 449                    NOM_CHAM='PRES_ELNO_IMAG',
 450                    NOEUD='NO780',
 451                    NOM_CMP='PRES',
 452                    VALE=0.0,
 453                    CRITERE='ABSOLU',
 454                    MAILLE='MA1',),
 455                 _F(NUME_ORDRE=1,
 456                    RESULTAT=PRESHARM,
 457                    NOM_CHAM='PRES_ELNO_IMAG',
 458                    NOEUD='NO751',
 459                    NOM_CMP='PRES',
 460                    VALE=1.6387,
 461                    CRITERE='RELATIF',
 462                    PRECISION=0.0033,
 463                    MAILLE='MA57',),
 464                 _F(NUME_ORDRE=1,
 465                    RESULTAT=PRESHARM,
 466                    NOM_CHAM='PRES_ELNO_IMAG',
 467                    NOEUD='NO763',
 468                    NOM_CMP='PRES',
 469                    VALE=1.6387,
 470                    PRECISION=0.0033,
 471                    MAILLE='MA57',),),);
 472
 473 FIN();