Télécharger pouve2.eso

Retour à la liste

Numérotation des lignes :

  1. C POUVE2 SOURCE CHAT 05/01/13 02:20:45 5004
  2. C POUVE2 SOURCE CASTEM 85/06/13 12:34:35
  3. SUBROUTINE POUVE2(U,V,P,ICOD)
  4. C-----------------------------------------------------------------------
  5. C ROUTINE DE PASSAGE DU VECTEUR DEPLACEMENT GLOBAL AU LOCAL
  6. C OU DU VECTEUR FORCE LOCAL AU GLOBAL
  7. C
  8. C CAS 2D - D.COMBESCURE 1998
  9. C
  10. C ENTREE
  11. C ------ U DEPLACEMENTS GLOBAUX AUX 2 NOEUDS
  12. C OU FORCES LOCALES AUX 2 NOEUDS
  13. C P MATRICE DE PASSAGE
  14. C ICOD 1 POUR Q-LOC =P*Q-GLOB
  15. C 2 POUR F-GLOB=PT*F-LOC
  16. C
  17. C SORTIE
  18. C ------ V DEPLACEMENTS LOCAUX AUX 2 NOEUDS
  19. C OU FORCES GLOBALES AUX 2 NOEUDS
  20. C-----------------------------------------------------------------------
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. DIMENSION U(6),V(6),P(2,2)
  24. C
  25. GO TO (100,200),ICOD
  26. 100 CONTINUE
  27. DO 6 J=1,2
  28. JJ=J + 3
  29. V(J)=0.D0
  30. V(JJ)=0.D0
  31. DO 6 I=1,2
  32. II=I+3
  33. V(J)=V(J) + P(J,I)*U(I)
  34. V(JJ)=V(JJ) + P(J,I)*U(II)
  35. 6 CONTINUE
  36. V(3) = U(3)
  37. V(6) = U(6)
  38. RETURN
  39. 200 CONTINUE
  40. DO 7 J=1,2
  41. JJ=J + 3
  42. V(J)=0.D0
  43. V(JJ)=0.D0
  44. DO 7 I=1,2
  45. II=I+3
  46. V(J)=V(J) + P(I,J)*U(I)
  47. V(JJ)=V(JJ) + P(I,J)*U(II)
  48. 7 CONTINUE
  49. V(3) = U(3)
  50. V(6) = U(6)
  51. RETURN
  52. END
  53.  
  54.  
  55.  

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