Télécharger la plaquette Cast3M


fr en

Télécharger adve_bcn.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  adve_bcn.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 =  DALLER LI1 LI2 LI3 LI4;
  32. MOD1= MODE SU THERMIQUE ISOTROPE advection;
  33.  
  34. * La densite de puissance thermique :
  35. vec1 = 'SOUR' MOD1 SU 4. ;
  36.  
  37. * Le vecteur vitesse :
  38. jxx  = MANU 'CHPO' SU 1 'UX' 1. ;
  39. JJxx = CHAN 'CHAM' Jxx MOD1 'RIGIDITE' 'CARACTERISTIQUES' ;
  40. jyy  = MANU 'CHPO' SU 1 'UY' 1. ;
  41. JJyy = CHAN 'CHAM' Jyy MOD1 'RIGIDITE' 'CARACTERISTIQUES' ; 
  42.  
  43. jj = mate MOD1 'K' 1. 'RHO' 1. 'C' 1. 'UX' jjxx 'UY' jjyy;
  44. matno1 = ADVE mod1 JJ;
  45. list matno1;
  46.  
  47. *Matrice élémentaire non-symétrique
  48. *-0.33333E+00  0.83333E-01  0.16667E+00  0.83333E-01 ;
  49. *-0.25000E+00  0.00000E+00  0.25000E+00  0.00000E+00 ;
  50. *-0.16667E+00 -0.83333E-01  0.33333E+00 -0.83333E-01 ;
  51. *-0.25000E+00  0.00000E+00  0.25000E+00  0.00000E+00 ;
  52.  
  53. *A linear system is solved with a non-symmetric matrix
  54.  
  55. * Les temperatures imposees :                                          *
  56. BL2 = BLOQ 'T' P2;
  57. BL4 = BLOQ 'T' P4;
  58. SOL7 = RESOU (matno1 et (bl2 et bl4)) vec1;
  59.  
  60. * Verification :
  61. TP1  = -2. ;
  62. TP3  = +2. ;
  63. ERR1 = ABS((EXTRAIRE SOL7 'T' P1) - TP1) ;
  64. ERR2 = ABS((EXTRAIRE SOL7 'T' P3) - TP3) ;
  65. SI ((ERR1+ERR2) > 1.E-10);
  66.   ERRE 5 ;
  67. * MESS '###### TEST NOT PASSED ######';
  68. SINON;
  69.   ERRE 0 ;
  70. * MESS '###### TEST PASSED ######';
  71. FINSI;
  72.  
  73. *****************************************
  74. matno1 = ADVE mod1 JJ 'SYMM';
  75. list matno1;
  76.  
  77. *Matrice élémentaire symétrique
  78. *-0.33333E+00 -0.83333E-01 0.00000E+00 -0.83333E-01
  79. *-0.83333E-01  0.00000E+00 0.83333E-01  0.00000E+00 
  80. * 0.00000E+00  0.83333E-01 0.33333E+00  0.83333E-01
  81. *-0.83333E-01  0.00000E+00 0.83333E-01  0.00000E+00
  82.  
  83. FIN;
  84.  
  85.  
  86.  

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