Télécharger h_b.procedur

Retour à la liste

Numérotation des lignes :

  1. * H_B PROCEDUR BAZE 97/01/30 21:15:07 2468
  2. 'DEBPROC' H_B MU0*'FLOTTANT' EVO1/'EVOLUTION' F/'FLOTTANT' ;
  3. *************************************************************************
  4. * definition de la courbe de permeabilite magnetostatique 2d 3d *
  5. * *
  6. * ENTREES :
  7. * MU0 permeabilité du vide ( S.I. 4 pi 10-7 ) *
  8. * EVO1 évolution abscisse H ordonne B *
  9. * F Packing factor 1 par defaut
  10. *************************************************************************
  11. * SORTIES :
  12. * *
  13. * MU(H) 3d methode a 2 potentiels scalaires *
  14. * 1./MU(B) 2d plan ou axi potentiel vecteur *
  15. * *
  16. * *
  17. *************************************************************************
  18.  
  19. 'SI' ('NON' ('EXISTE' EVO1 )) ;
  20. TTTT = 'B(H)__STANDARD' ;
  21. MESS '***********************************************************';
  22. MESS ' ON UTILISE UNE COURBE B(H) STANDARD TESLA-AMPERE/METRE ' ;
  23. MESS '***********************************************************';
  24. B= 'PROG' 0. .8 1.2 1.4 1.5 'PAS' .05 2.15 2.175 2.2 2.25 2.28 2.3
  25. 2.3443
  26. 2.3996 2.4905 2.5627 2.6706 2.8498 3.2074 3.5644 4.2782 4.8134 5.7052
  27. 6.4186 7.4887 17.48 27.48 ;
  28. H= 'PROG' 0. 159.2 294.4 501.3 795.8 1154. 1795. 2862. 4383. 6044.
  29. 8122.
  30. 10590. 13610. 17220. 21170. 26750. 33760. 43800. 52440. 66000. 99470.
  31. 120960. 141210. 169600. 212170. 283130. 339890. 425040. 566950.
  32. 850760.
  33. 1134600. 1702300. 2128000. 2837700. 3405100. 4256700. 12215578.
  34. 20174457. ;
  35. 'SINON' ;
  36. TTTT = 'B(H)__UTILISATEUR' ;
  37. H = 'EXTR' EVO1 'ABSC' ; B = 'EXTR' EVO1 'ORDO' ;
  38. 'FINSI' ;
  39. 'SI' ('NON' ('EXISTE' F)) ; F= 1. ;'FINSI' ;
  40.  
  41. TUTU = CHAIN TTTT '_______PACKING ' F ;
  42. 'MESS' '********************************************' ;
  43. 'MESS' TUTU ;
  44. 'MESS' '********************************************' ;
  45. NN = DIME H ;
  46. H_MAX = MAXI H ;
  47. B_MAX = MAXI B ;
  48. B_EXT = 5.e4 ; H_EXT = H_MAX + (( B_EXT - B_MAX)/MU0 ) ;
  49. H = H ET ( 'PROG' H_EXT) ;
  50. B = B ET ( 'PROG' B_EXT) ;
  51. 'SI' (f < 1.) ;
  52. ip = 1 ;
  53. BRED = PROG 0. ;
  54. 'REPETER' BUB nn ;
  55. ip = ip + 1 ;
  56. Bi = EXTR B IP ; Hi = EXTR H IP ; MUI = Bi/Hi ;
  57. MURED =( MU0 * ( 1. - f )) + (MUI * f ) ;
  58. BRED = BRED et (PROG ( Hi * MURED)) ;
  59. 'FIN' BUB ;
  60. B = BRED ;
  61. 'FINSI' ;
  62.  
  63. BHEVO=EVOL MANU 'H' H 'B' B ;
  64. * DESS BHEVO;
  65. H_1= 'ENLEVER' H 1 ;B_1= 'ENLEVER' B 1 ;
  66. MUV= H_1 / B_1;TU= 'EXTRA' MUV 1;MUV= ( 'PROG' TU ) 'ET' MUV ;
  67. IFOR = 'VALE' DIME ;
  68. 'SI' ('EGA' IFOR 3 );
  69. BOBO= 'TEXTE' ' MU(H)______POTENTIEL_SCALAIRE' ;
  70. 'TITRE' BOBO ;
  71. REVOL = 'EVOL' MANU 'TEMPERATURE ' H 'CONDUCTIVITE ' ( MUV ** -1.) ;
  72. 'SINON' ;
  73. BOBO= 'TEXTE' ' MU(B)______POTENTIEL_VECTEUR ';
  74. 'TITRE' BOBO ;
  75. REVOL='EVOL' 'MANU' 'TEMPERATURE' B 'CONDUCTIVITE' MUV ;
  76. 'FINSI' ;
  77. *
  78. 'SAUTER' 1 'LIGNE' ;
  79. 'MESS' '********************************************' ;
  80. 'MESS' BOBO ;
  81. 'MESS' '********************************************' ;
  82. 'SAUTER' 1 'LIGNE' ;
  83. 'FINPROC' REVOL ;
  84.  
  85.  
  86.  

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