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 :
© Cast3M 2003 - All rights reserved.
Disclaimer