Télécharger la plaquette Cast3M


fr en

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