Exemples/ex19/titi.comm

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