1 : C TRANSG SOURCE CHAT 05/01/13 03:44:57 5004 2 : SUBROUTINE TRANSG(VKE,B,LRE,NBNO,ICOM) 3 : C===================================================== 4 : C 5 : C TRANSFORME VKE(LRE,LRE) DU REPERE LOCAL VKE ANTISYMETRIQUE 6 : C VKE = CORIOLIS 7 : C AU REPERE GENERAL A L AIDE DE LA MATRICE DE PASSAGE 8 : C B(3,3) 9 : C NBNO NOMBRE DE NOEUDS DE L ELEMENT 3 NOEUDS DKT 10 : C 4 NOEUDS DKQ 11 : C LRE NOMBRE DE DDL DE L ELEMENT 12 : C 6 DDL PAR NOEUD 13 : C ENTREES 14 : C VKE(LRE,LRE) = MATRICE DE RAIDEUR APRES CHGMNT DE REPERE 15 : C B(3,3) = MATRICE DE PASSAGE 16 : C LRE = NOMBRE DE D.D.L. DANS LA MATRICE DE RIGIDITE 17 : C NBNO = NOMBRE DE NOEUDS ( LRE = NBNO*6 18 : C ICOM = 0 SI LES D.D.L. DE FLEXION SONT DECOUPLES DES 19 : C D.D.L. DE MEMBRANE 1 SINON 20 : C SORTIES 21 : C VKE(LRE,LRE) = MATRICE DE CORIOLIS APRES CHGMNT DE REPERE 22 : C 23 : C D'APRES TRANSK COMBESCURE 2004 24 : C===================================================== 25 : IMPLICIT INTEGER(I-N) 26 : IMPLICIT REAL*8(A-H,O-Z) 27 : DIMENSION VKE(LRE,*),B(3,*) 28 : DIMENSION AA(3,3),A(3,3) 29 : DIMENSION IMLT(36) 30 : DATA IMLT/1,2,1,4,3,1,2,4,2,1,4,3,4,3,1,2,4,2,4,2,1, 31 : 1 4,3,4,3,4,3,1,2,4,2,4,2,4,2,1/ 32 : NB=NBNO+NBNO 33 : IJ=1 34 : JTROIS=0 35 : DO 100 J=1,NB 36 : ITROIS=0 37 : DO 200 I=1,NB 38 : ITROIS =3*I - 3 39 : JTROIS =3*J - 3 40 : C TRANSVASE UN BLOC 3 3 DE VKE DANS LA MATRICE DE TRAVAIL A 41 : DO 300 IA=1,3 42 : DO 300 IB=1,3 43 : A(IA,IB)=VKE(ITROIS+IA,JTROIS+IB) 44 : 300 CONTINUE 45 : C MATRICE NON SYMETRIQUE 45 : CALL VMULT4(A,AA,B) 46 : 350 CONTINUE 47 : C TRANSVASE LA MATRICE AA DANS LE BLOC 3 3 DE VKE 48 : DO 400 IA=1,3 49 : DO 400 IB=1,3 50 : VKE(ITROIS+IA,JTROIS+IB)=AA(IA,IB) 51 : C VKE(JTROIS+IB,ITROIS+IA)=(-1.D0)*AA(IA,IB) 52 : 400 CONTINUE 53 : IJ=IJ+1 54 : 200 CONTINUE 55 : 100 CONTINUE 56 : C 200 ITROIS=ITROIS+3 57 : C 100 JTROIS=JTROIS+3 58 : RETURN 59 : END 60 : 61 :
© Cast3M 2003 - All rights reserved.
Disclaimer