Télécharger fronabs.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fronabs.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************************
  5. *
  6. * test des frontieres absorbantes
  7. *
  8. * on teste la reultante pour un champ de vitesses donne
  9. *
  10. *
  11. ***************************************************************
  12.  
  13. option dime 3 elem cub8;
  14.  
  15. E = 1.;
  16. POIS = 0. ;
  17. RO = 1.;
  18.  
  19. P0 = 0. 0. 0. ;
  20. P1 = 1. 0. 0. ;
  21. P2 = 0. 0. -1. ;
  22.  
  23. D1 = P0 DROI 5 P1 ;
  24. D2 = P0 DROI 5 P2 ;
  25. C1 = CERC 5 P1 P0 P2 ;
  26.  
  27. SUR1 = SURF (D1 ET D2 ET C1) PLANE ;
  28.  
  29. SOL1 = SUR1 VOLU 5 ROTA 90. ( 0.0 0.0 0.0 ) ( 0.0 0.0 1.0 ) ;
  30.  
  31. *TRAC QUAL SOL1 cach;
  32.  
  33. BOR1 = C1 ROTA 5 90. ( 0.0 0.0 0.0 ) ( 0.0 0.0 1.0 ) ;
  34.  
  35. ELIM 0.01 (BOR1 et SOL1) ;
  36.  
  37. MOD_S = MODE SOL1 MECANIQUE ;
  38. MAT_S1 = MATE MOD_S YOUN E NU POIS RHO RO ;
  39.  
  40. AMOT1 = AMOR MOD_S BOR1 MAT_S1 ;
  41.  
  42. E = 2.;
  43. POIS = 0. ;
  44. RO = 1.;
  45.  
  46. MAT_S2 = MATE MOD_S YOUN E NU POIS RHO RO ;
  47.  
  48. AMOT2 = AMOR MOD_S BOR1 MAT_S2 ;
  49.  
  50. X1 = COOR 1 BOR1 ;
  51. Y1 = COOR 2 BOR1 ;
  52. Z1 = COOR 3 BOR1 ;
  53. X1 = NOMC 'UX' X1;
  54. Y1 = NOMC 'UY' Y1;
  55. Z1 = NOMC 'UZ' Z1;
  56.  
  57. DEP1 = X1 + Y1 ;
  58. DEP1 = DEP1 + Z1 ;
  59.  
  60. RIG1 = RELA CORI DEPL SOL1 ;
  61. RIG2 = BLOQ (P0 et P2) DEPL;
  62. RIG3 = BLOQ P1 UY ;
  63. F1 = DEPI RIG3 1. ;
  64. DEP2 = RESO (RIG1 ET RIG2 et rig3 ) F1 ;
  65.  
  66. vec1 = vect dep2 ux uy uz 1. jaun ;
  67. *trac vec1 sol1 cach ;
  68.  
  69. FOR1=RESUL (AMOT1 * DEP1);
  70. FOR2=RESUL (AMOT2 * DEP2);
  71. for1 = (xtx for1)**.5 ;
  72. for2 = (xtx for2)**.5 ;
  73. ** resultat analytiques
  74. fan1 = 3**.5 * pi / 4 ; fan2 = 2**.5 * pi / 4 ;
  75. err1 = for1 - fan1 / fan1 * 100 ;
  76. err2 = for2 - fan2 / fan2 * 100 ;
  77. errmax = maxi abs (prog err1 err2) ;
  78.  
  79. list errmax ;
  80.  
  81. MESS 'La résultante des forces visqueuses dues au déplacement
  82. radial est' FOR1;
  83.  
  84. MESS 'La résultante des forces visqueuses dues au déplacement
  85. tangentiel est' FOR2;
  86.  
  87. SI(errmax <EG 4.);
  88. ERRE 0;
  89. SINO;
  90. ERRE 5;
  91. FINSI;
  92.  
  93. FIN;
  94.  
  95.  
  96.  
  97.  
  98.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales