Télécharger jacob.eso

Retour à la liste

Numérotation des lignes :

  1. C JACOB SOURCE CHAT 05/01/13 00:48:20 5004
  2. C***********************************************************************
  3. C* *
  4. C* PROJET : Opérateur LIMI *
  5. C* NOM : JACOB, sous-programme *
  6. C* DESCRIPTION : Ce sous-programme calcule la matrice jacobienne et *
  7. C* son inverse, pour le système de 3 équations à trois *
  8. C* inconnues utilisé dans la méthode de résolution des *
  9. C* relations de Cousteix *
  10. C* LANGAGE : Esope *
  11. C* AUTEUR : Guillaume VENCO - DRN/DMT/SEMT/LTMF *
  12. C* *
  13. C***********************************************************************
  14. C* *
  15. C* APPELES : Sous-programme INVMAT3.eso *
  16. C* *
  17. C***********************************************************************
  18. C* *
  19. C* ENTREES : G : -1ère inconnue du système, paramètre de Clauser *
  20. C* H : -2ème inconnue, facteur de forme *
  21. C* ga : -3ème inconnue, facteur de frottement gamma *
  22. C* *
  23. C* SORTIES : W : -Matrice jacobienne du système *
  24. C* Q : -Matrice inverse de la matrice jacobienne *
  25. C* *
  26. C***********************************************************************
  27. C* *
  28. C* VERSION : 20/04/2000 *
  29. C* CREATION : *
  30. C* *
  31. C***********************************************************************
  32.  
  33. SUBROUTINE JACOB(H,G,GA,W,Q)
  34.  
  35. C* *********************************
  36. C* *** Déclaration des variables ***
  37. C* *********************************
  38. IMPLICIT INTEGER(I-N)
  39.  
  40. REAL*8 G,H,ga,W(3,3),Q(3,3)
  41. REAL*8 F1,DF1,D,DD,beta,dbeta
  42.  
  43. C* *********************************
  44. C* ******** Sous-programme *********
  45. C* *********************************
  46.  
  47. C---- Calcul des fonctions de similitude F et D et leurs dérivées
  48. F1=0.613*G-(3.6+76.86*(((1./G)-0.154)**2))/G
  49. c write(6,*) 'F1(G)=',F1
  50. DF1=0.613+(1./(G**2))*
  51. & ((2*76.86*((1./G)-0.154)/G)+(3.6+76.86*(((1./G)-0.154)**2)))
  52. c write(6,*) 'F1''(G)=',DF1
  53. D=2*G-4.25*sqrt(G)+2.12
  54. c write(6,*) 'D*(G)=',D
  55. DD=2-4.25/(2*sqrt(G))
  56. c write(6,*) 'D*''(G)=',DD
  57.  
  58. C---- Calcul du cofficient intervenant dans la première équation
  59. beta=(G/F1)-1
  60. dbeta=(F1-G*DF1)/(F1**2)
  61.  
  62. C---- Matrice jacobienne, dérivées partielles des équations
  63. W(1,1)=-(beta*(H**2)-2*beta*H-1)/((H-1)**2)
  64. W(1,2)=-dbeta*(H**2)/(H-1)
  65. W(1,3)=0.
  66. W(2,1)=0.
  67. W(2,2)=DD
  68. W(2,3)=1./(ga**2)
  69. W(3,1)=-1./(ga*(H**2))
  70. W(3,2)=1.
  71. W(3,3)=(H-1)/(H*(ga**2))
  72.  
  73. C---- Inverse de la matrice jacobienne
  74. CALL INV(W,Q)
  75.  
  76. c write(6,1004) Q(1,1),Q(1,2),Q(1,3)
  77. c write(6,1004) Q(1,1),Q(1,2),Q(1,3)
  78. c write(6,1004) Q(1,1),Q(1,2),Q(1,3)
  79.  
  80. RETURN
  81.  
  82. c 1004 FORMAT(3(1X,1PE11.4))
  83.  
  84. END
  85.  
  86. C* *********************************
  87. C* **** Fin Sous-Programme *******
  88. C* *********************************
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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