1 : C MTOGYF SOURCE DC 05/02/15 21:15:51 5030 2 : SUBROUTINE MTOGYF(LRE,NDDL,REL,VROT,RELB) 3 : C======================================================================= 4 : C 5 : C Calcul de la matrice du couplage de Coriolis à partir de la matrice de masse 6 : C pour les éléments massifs en 3D 7 : C 8 : C Didier COMBESCURE Janvier 2005 9 : C 10 : C======================================================================= 11 : IMPLICIT INTEGER(I-N) 12 : IMPLICIT REAL*8 (A-H,O-Z) 13 : DIMENSION REL(LRE,1),RELB(2*LRE,1),VROT(*),ROME(6,6) 14 : C 15 : NNOEU = LRE/NDDL 16 : C 17 : C 18 : DO 10 I=1,6 19 : DO 10 J=1,6 20 : ROME(I,J)=0.D0 21 : 10 CONTINUE 22 : * 23 : ROME(1,6) = (-1.)*VROT(2) 24 : ROME(2,5) = (-1.)*VROT(2) 25 : ROME(6,1) = VROT(2) 26 : ROME(5,2) = VROT(2) 27 : C 28 : C 29 : DO 100 IN=1,NNOEU 30 : C 31 : DO 100 IM=1,NNOEU 32 : C 33 : DO 100 IDD1=1,NDDL 34 : C 35 : DO 100 IDD2=1,NDDL 36 : C 37 : RELB((6*IN)-7+2*IDD1,(6*IM)-7+2*IDD2)= 38 : . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2-1) 39 : RELB((6*IN)-6+2*IDD1,(6*IM)-6+2*IDD2)= 40 : . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2) 41 : RELB((6*IN)-6+2*IDD1,(6*IM)-7+2*IDD2)= 42 : . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2-1) 43 : RELB((6*IN)-7+2*IDD1,(6*IM)-6+2*IDD2)= 44 : . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2) 45 : C 46 : 100 CONTINUE 47 : C 48 : RETURN 49 : END 50 : 51 : 52 :
© Cast3M 2003 - All rights reserved.
Disclaimer