Exemples/ex23/toto.comm

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