Télécharger adve_03.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : adve_03.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * *
  5. * On resout : v.gradT = 4, *
  6. * *
  7. * Avec : vx=1, vy=1 *
  8. * on a : dT/dx + dT/dy = 4 *
  9. * *
  10. * De + : T(1,0)=0, T(0,1)=0 *
  11. * Sur 1 element QUA4 de geometrie carree et de cote 1. *
  12. * *
  13. * La solution analytique est : T(x,y) = 2(x+y-1) *
  14. * Soit : T(0,0)=T(P1)=-2 et T(1,1)=T(P3)=2 *
  15. * *
  16. ***************************************************
  17. ** TESTING FILE FOR THE OPERATOR ADVE BUILDING A *
  18. ** NON-SYMMETRIC STIFFNESS MATRIX *
  19. ***************************************************
  20.  
  21. OPTI DIME 2 ELEM qua4 MODE PLAN CONT;
  22. p1 = 0. 0.;
  23. p2 = 1. 0.;
  24. p3 = 1. 1.;
  25. p4 = 0. 1.;
  26. LI1 = droite 1 p1 p2;
  27. LI2 = droite 1 p2 p3;
  28. LI3 = droite 1 p3 p4;
  29. LI4 = droite 1 p1 p4;
  30.  
  31. SU = DALL LI1 LI2 LI3 LI4;
  32. MOD1 = MODE SU THERMIQUE ISOTROPE advection;
  33.  
  34. jj = MATE MOD1 'VITX' 1. 'VITY' 1. ;
  35. matno1 = ADVE mod1 JJ;
  36. list matno1;
  37.  
  38. *Matrice élémentaire non-symétrique
  39. *-0.33333E+00 0.83333E-01 0.16667E+00 0.83333E-01 ;
  40. *-0.25000E+00 0.00000E+00 0.25000E+00 0.00000E+00 ;
  41. *-0.16667E+00 -0.83333E-01 0.33333E+00 -0.83333E-01 ;
  42. *-0.25000E+00 0.00000E+00 0.25000E+00 0.00000E+00 ;
  43.  
  44. * La densite de puissance thermique :
  45. vec1 ='SOUR' MOD1 SU 4. ;
  46.  
  47. * Les temperatures imposees : *
  48. BL2 = BLOQ 'T' P2;
  49. BL4 = BLOQ 'T' P4;
  50.  
  51. *A linear system is solved with a non-symmetric matrix
  52. SOL7 = RESOU (matno1 et (bl2 et bl4)) vec1;
  53.  
  54. * Verification :
  55. TP1 = -2. ;
  56. TP3 = +2. ;
  57. ERR1 = ABS((EXTRAIRE SOL7 'T' P1) - TP1) ;
  58. ERR2 = ABS((EXTRAIRE SOL7 'T' P3) - TP3) ;
  59. SI ((ERR1+ERR2) > 1.E-10);
  60. ERRE 5 ;
  61. * MESS '###### TEST NOT PASSED ######';
  62. SINON;
  63. ERRE 0 ;
  64. * MESS '###### TEST PASSED ######';
  65. FINSI;
  66.  
  67. *****************************************
  68. matno1 = ADVE mod1 JJ 'SYMM';
  69. list matno1;
  70.  
  71. *Matrice élémentaire symétrique
  72. *-0.33333E+00 -0.83333E-01 0.00000E+00 -0.83333E-01
  73. *-0.83333E-01 0.00000E+00 0.83333E-01 0.00000E+00
  74. * 0.00000E+00 0.83333E-01 0.33333E+00 0.83333E-01
  75. *-0.83333E-01 0.00000E+00 0.83333E-01 0.00000E+00
  76.  
  77. FIN;
  78.  
  79.  
  80.  

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