Télécharger la plaquette Cast3M


fr en

Télécharger test_thermique_1D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : test_thermique_1d.dgibi
  2. ************************************************************************
  3. * Cas test THERMIQUE 1D : test des operateurs en 1D - a completer !
  4. ************************************************************************
  5. TITRE 'THER1D - C.L. TEMPERATURE';     
  6.  
  7. OPTION 'DIME' 1 'ELEM' SEG2;
  8.  
  9. XMAX = 1.;
  10. P0 = 'POIN' 0.;
  11. P1 = 'POIN' XMAX;
  12. NX = 5;
  13. D1 = 'DROIT' NX P0 P1;
  14. opti sort 'testavs1D' ; sort avs D1 ;
  15. ** opti sort 'toto' ;
  16. *'TRACER' D1;
  17.  
  18. * DONNEES DU PB
  19. UI = 0;
  20. FI = 25;
  21. COND = 130.;
  22.  
  23. * MODELE
  24. MOD1 = MODE D1 THERMIQUE ISOTROPE ;
  25. * MATERIAU                                                                      
  26. MAT1 = MATE MOD1 K COND ;
  27.  
  28. * MATRICE DE CONDUCTIVITE
  29. CND1 = CONDUCTIVITE MOD1 MAT1 ;
  30.  
  31. * CL DIRICHLET
  32. B0 = BLOQUE 'T   ' P0 ;    
  33. B1 = BLOQUE 'T   ' P1 ;
  34.  
  35. * TERME SOURCE
  36. S  = 'MANUEL' CHPO D1 'F' FI;
  37. ES = 'SOURCE' MOD1 S ;
  38.  
  39. * ASSEMBLAGE DES PREMIER ET SECOND MEMBRES:
  40. RIG1 = CND1 'ET'B0 'ET' B1;
  41. FF1 = ES;
  42.  
  43. CHTER1 = RESOUDRE RIG1 FF1 ;
  44. *'TRACER' CHTER1 D1;
  45.  
  46. *calcul de l'erreur
  47. *-- sol analytique
  48. *     T(X) = -FI/2K * X  * (X '-' 1)
  49. X= 'COORDONNEE' 1 D1;
  50. CHUN = 'MANUEL' CHPO D1 'SCAL' 1;
  51. CHP = X*(X-CHUN);
  52. SOL = CHP*(FI/(2*COND))*(-1);
  53. *'TRACER'SOL D1;
  54. ECART = ('NOMC' 'T' SOL) '-'CHTER1;
  55. NOR = ('XTX'ECART '/' NX )'**'(0.5);
  56. mess ' Calcul de l erreur NOR=' NOR ;
  57. * => super-convergence
  58. * -- non super convergence en un noeud hors discretisation
  59.  
  60. * projection valeur en 0.5
  61. pm = 'POIN' 0.5;
  62. MPM = 'MANUEL' 'POI1' pm;
  63. chter2 = 'CHANGER' 'CHAM' CHTER1 MOD1 'NOEUD';
  64. val = proi mpm chter2;
  65. list val;
  66. * => 2.30769E-02
  67. * valeur analytique = 2.4038 E-02
  68. *****
  69.  
  70. * Conditions de Neumann homogene partout
  71. MP0= 'MANUEL' 'POI1' P0;
  72. MP1= 'MANUEL' 'POI1' P1;
  73. F0 = FLUX MOD1 0. MP0 ;
  74. F1 = FLUX MOD1 0. MP1 ;
  75. * ASSEMBLAGE DES PREMIER ET SECOND MEMBRES:
  76. RIG1 = CND1 ;
  77. FF1 = F0 'ET' F1;
  78. CHTER1 = RESOUDRE RIG1 FF1 ;
  79. * Pas de blocages donc un mode rigide en T -> bloquer a 0 par RESOU -> solution constante nulle !
  80. * Si terme source nul : FF1 = F0 'ET' F1 'ET' ES ; -> pas de solution (2 points a bloquer en T)
  81. *'TRACER' CHTER1 D1;
  82.  
  83. opti lect 'testavs1D.inp' ; mail1D = lire avs ;
  84. list mail1D.'LEMAILLA' ;
  85.  
  86. 'FIN';
  87.  
  88.  
  89.  
  90.  
  91.  

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