Télécharger biovol.procedur

Retour à la liste

Numérotation des lignes :

  1. * BIOVOL PROCEDUR FANDEUR 14/10/10 21:15:02 8178
  2. 'DEBPROC' BIOVOL BOBI*'MAILLAGE' PCAL*'MAILLAGE' J*'CHPOINT'
  3. MU0* 'FLOTTANT' ;
  4. ***********************************************************************
  5. * CALCUL DU CHAMP BS CREE SUR PCAL PAR L INDUCTEUR VOLUMIQUE BOBI *
  6. * PARCOURU PAR DES COURANTS DES DENSITES J (COMPOSANTES JX JY JZ) *
  7. * MU0 PERMEABILITE DU VIDE *
  8. ***********************************************************************
  9. fact = 1. / (4. * pi) ;
  10. nni = bobi 'NBNO' ;
  11. JIX = 'EXCO' J 'JX' ;
  12. JIY = 'EXCO' J 'JY';
  13. JIZ = 'EXCO' J 'JZ' ;
  14. *
  15. X Y Z = 'COOR' BOBI ;
  16. NGEO = PCAL 'NBNO' ;IGEO = 0;
  17. 'REPETER' BO_GEO NGEO ;
  18. IGEO = IGEO + 1 ;
  19. PCAL1 = PCAL 'POINT' IGEO ;
  20. XP YP ZP = 'COOR' PCAL1 ;
  21. RX = XP - X ;RY = YP - Y ; RZ = ZP - Z ;
  22. RM2 = (RX * RX) + ( RY * RY ) + ( RZ * RZ ) ;
  23. R3 = RM2 ** 1.5 ;
  24. HX = (((JIY * RZ ) - ( JIZ * RY) )/ R3) ;
  25. HY = (((JIZ * RX ) - ( JIX * RZ) )/ R3) ;
  26. HZ = (((JIX * RY ) - ( JIY * RX ) )/ R3) ;
  27. HPT =(HX 'NOMC' 'HX')+( HY 'NOMC' 'HY' )+( HZ 'NOMC' 'HZ' ) ;
  28. OBJM = 'MODE' bobi 'THERMIQUE' 'ISOTROPE' ;
  29. HEPT = 'CHAN' CHAM HPT OBJM 'NOEUDS' ;
  30. hheX = 'EXCO' hept 'HX' ;
  31. hhey = 'EXCO'hept 'HY' ;
  32. hhez = 'EXCO' hept 'HZ' ;
  33. HTX = 'INTG' HHEX OBJM ;
  34. HTY = 'INTG' HHEY OBJM ;
  35. HTZ = 'INTG' HHEZ OBJM ;
  36. *
  37. HT ='MANU' 'CHPO' pcal1 3 'HX' HTX 'HY' HTY 'HZ' HTZ ;
  38. 'SI' (IGEO 'EGA' 1) ;
  39. HGEO = HT ;
  40. 'SINON' ;
  41. HGEO =HGEO + HT ;
  42. 'FINSI' ;
  43. 'FIN' BO_GEO ;
  44. ******************************************************
  45. HGEO = HGEO * fact ;
  46. 'FINPROC' HGEO ;
  47.  
  48.  

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