Télécharger la plaquette Cast3M


fr en

Télécharger joiglv.eso

Retour à la liste

Numérotation des lignes : 

joiglv

  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