Télécharger la plaquette Cast3M


fr en

Télécharger jacob.eso

Retour à la liste

Numérotation des lignes : 

jacob

  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