Télécharger joiglv.eso

Retour à la liste

Numérotation des lignes :

  1. C JOIGLV SOURCE AM 11/05/03 21:17:16 6955
  2. SUBROUTINE JOIGLV(REL,BPSS,RELO,P,LRE,IDIM)
  3. C-----------------------------------------------------------------------
  4. C CALCUL DU VECTEUR CONTENANT LES CONTRAINTES DANS LE REPERE GLOBAL
  5. C
  6. C ENTREE
  7. C REL VECTEUR EXPRIME DANS LE REPERE GLOBAL
  8. C BPSS MATRICE DE PASSAGE (IDIM x IDIM)
  9. C SORTIE
  10. C REL VECTEUR EXPRIME DANS LE REPERE LOCAL
  11. C-----------------------------------------------------------------------
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8(A-H,O-Z)
  14. DIMENSION REL(LRE),BPSS(3,3),RELO(LRE),P(LRE,LRE)
  15. *
  16. CALL ZERO(P,LRE,LRE)
  17. CALL ZERO(RELO,LRE,1)
  18. *
  19. * ASSEMBLAGE DE PASSAGE P
  20. *
  21. IF(IDIM.EQ.3) THEN
  22. DO K=0,3
  23. DO I=1,3
  24. DO J=1,3
  25. P(I+3*K,J+3*K)=BPSS(J,I)
  26. ENDDO
  27. ENDDO
  28. ENDDO
  29. ELSE IF(IDIM.EQ.2) THEN
  30. DO K=0,1
  31. DO I=1,3
  32. DO J=1,3
  33. P(I+3*K,J+3*K)=BPSS(J,I)
  34. ENDDO
  35. ENDDO
  36. ENDDO
  37. ENDIF
  38. C
  39. C REG = P*REL
  40. C
  41. DO I=1,LRE
  42. DO J=1,LRE
  43. RELO(I) = RELO(I) + P(I,J)*REL(J)
  44. ENDDO
  45. ENDDO
  46. *
  47. CALL ZERO(REL,LRE,1)
  48. C
  49. C REL = RELO
  50. C
  51. DO I=1,LRE
  52. REL(I) = RELO(I)
  53. ENDDO
  54. *
  55. RETURN
  56. END
  57.  
  58.  
  59.  

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