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 : * 15 : CALL ZERO(P,LRE,LRE) 15 : CALL ZERO(RELO,LRE,1) 16 : * 17 : * ASSEMBLAGE DE PASSAGE P 18 : * 19 : IF(IDIM.EQ.3) THEN 20 : DO K=0,3 21 : DO I=1,3 22 : DO J=1,3 23 : P(I+3*K,J+3*K)=BPSS(J,I) 24 : ENDDO 25 : ENDDO 26 : ENDDO 27 : ELSE IF(IDIM.EQ.2) THEN 28 : DO K=0,1 29 : DO I=1,3 30 : DO J=1,3 31 : P(I+3*K,J+3*K)=BPSS(J,I) 32 : ENDDO 33 : ENDDO 34 : ENDDO 35 : ENDIF 36 : C 37 : C REG = P*REL 38 : C 39 : DO I=1,LRE 40 : DO J=1,LRE 41 : RELO(I) = RELO(I) + P(I,J)*REL(J) 42 : ENDDO 43 : ENDDO 44 : * 44 : CALL ZERO(REL,LRE,1) 45 : C 46 : C REL = RELO 47 : C 48 : DO I=1,LRE 49 : REL(I) = RELO(I) 50 : ENDDO 51 : * 52 : RETURN 53 : END 54 : 55 :
© Cast3M 2003 - All rights reserved.
Disclaimer