1 : C POURO2 SOURCE CHAT 05/01/13 02:19:46 5004 2 : C POUROT SOURCE CASTEM 85/06/13 12:34:25 3 : SUBROUTINE POURO2(REL,LRE,P,RP) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES 6 : C GLOBAUX EN 2D 7 : C D APRES TEDEL 8 : C 9 : C ENTREE 10 : C LRE NOMBRE DE LIGNES DU TABLEAU REL 11 : C P MATRICE DE PASSAGE 12 : C RP MATRICE ELEMENTAIRE 12X12 EN AXES LOCAUX 13 : C 14 : C SORTIE 15 : C REL MATRICE ELEMENTAIRE EN AXES GLOBAUX 16 : C 17 : C VIENT DE TEDEL AOUT 84 18 : C 19 : C D'APRES POUROT DC 98 20 : C----------------------------------------------------------------------- 21 : IMPLICIT INTEGER(I-N) 22 : IMPLICIT REAL*8(A-H,O-Z) 23 : DIMENSION REL(LRE,1),P(2,2),RP(6,6) 24 : DIMENSION NI(6),NJ(6) 25 : DATA ZER,UN/0.D0,1.0D0/ 26 : C 27 : DATA NI/1,1,1,4,4,4/ 28 : DATA NJ/1,2,3,1,2,3/ 29 : DO 10 K=1,6 30 : I3=NI(K) 31 : I4=I3+2 32 : LK=NJ(K) 33 : DO 11 L=K,6 34 : REL(K,L)= 0.D00 35 : I1=NI(L) 36 : I2=I1+2 37 : LL=NJ(L) 38 : LM=0 39 : DO 12 M=I3,I4 40 : LM=LM+1 41 : S=0. 42 : II=0 43 : DO 13 N=I1,I2 44 : II=II+1 45 : C 46 : IF (II.EQ.3) THEN 47 : IF (LL.EQ.3) THEN 48 : XFACT1 =UN 49 : ELSE 50 : XFACT1 = ZER 51 : ENDIF 52 : ELSEIF (LL.EQ.3) THEN 53 : IF (II.EQ.3) THEN 54 : XFACT1 =UN 55 : ELSE 56 : XFACT1 = ZER 57 : ENDIF 58 : ELSE 59 : XFACT1 = P(II,LL) 60 : ENDIF 61 : C 62 : IF (LM.EQ.3) THEN 63 : IF (LK.EQ.3) THEN 64 : XFACT2 =UN 65 : ELSE 66 : XFACT2 = ZER 67 : ENDIF 68 : ELSEIF (LK.EQ.3) THEN 69 : IF (LM.EQ.3) THEN 70 : XFACT2 =UN 71 : ELSE 72 : XFACT2 = ZER 73 : ENDIF 74 : ELSE 75 : XFACT2 = P(LM,LK) 76 : ENDIF 77 : C 78 : 13 S=S+ RP(M,N)*XFACT1 79 : 12 REL(K,L)= REL(K,L)+ S*XFACT2 80 : 11 REL(L,K)= REL(K,L) 81 : 10 CONTINUE 82 : RETURN 83 : END 84 : 85 : 86 :
© Cast3M 2003 - All rights reserved.
Disclaimer