1 : C CUBHR2 SOURCE CHAT 05/01/12 22:32:40 5004 2 : C CUBHR2 SOURCE LOFT 88/06/23 21:07:05 3 : SUBROUTINE CUBHR2(XGENE,DJAC,VKL1,VKL2,LRE,NST,NBNO,IFOU,REL) 4 : C======================================================================= 5 : C 6 : C INTEGRE LE PRODUIT NT KL N 7 : C POUR LES ELEMENTS HOMOGENEISES TRIH 8 : C 9 : C ROUTINE FORTRAN PUR 10 : C 11 : C======================================================================= 12 : C NBDL = NOMBRE DE DDL PAR NOEUD 13 : C INPUT 14 : C XGENE(NBNO,NST)=MATRICE N 15 : C DJAC =JACOBIEN 16 : C VKL1=(COEFPR**2)/(RHO*C**2)*SFLU/SCEL 17 : C VKL2= EI/((EPS**2)*SCEL) 18 : C LRE=TAILLE DE LA MATRICE DE RIGIDITE 19 : C NST=NOMBRE DE COLONNES DE LA MATRICE XGENE 20 : C NBNO=NOMBRE DE NOEUDS DE L ELEMENT CONSIDERE 21 : C IFOU=IFOUR DE CCOPTIO 22 : C OUTPUT 23 : C REL=MATRICE DE RIGIDITE 24 : C======================================================================= 25 : IMPLICIT INTEGER(I-N) 26 : IMPLICIT REAL*8 (A-H,O-Z) 27 : DIMENSION REL(LRE,*),XGENE((NBNO+NBNO),*) 28 : NBDL=LRE/NBNO 29 : C 30 : C 31 : C TERMES EN P*P 32 : C 33 : IX1=0 34 : IY1=0 35 : DO 300 IX=1,LRE ,NBDL 36 : IX1=IX1 + 1 37 : DO 400 IY=1,IX ,NBDL 38 : IY1=IY1 + 1 39 : REL(IX,IY) = REL(IX,IY) + DJAC*VKL1*XGENE(IX1,1)*XGENE(IY1,1) 40 : REL(IY,IX) = REL(IX,IY) 41 : 400 CONTINUE 42 : IY1=0 43 : 300 CONTINUE 44 : C 45 : C 46 : C TERMES EN (UX,RY )*(UX,RY ) 47 : C 48 : IX1=0 49 : IY1=0 50 : DO 302 IX=3,LRE ,NBDL 51 : IX1=IX1 + 1 52 : DO 402 IY=3,IX ,NBDL 53 : IY1=IY1 + 1 54 : REL(IX,IY) = REL(IX,IY) + DJAC*VKL2*XGENE(IX1,2)*XGENE(IY1,2) 55 : REL(IX+1,IY) = REL(IX+1,IY) + DJAC*VKL2*XGENE(IX1+8,2) 56 : #*XGENE(IY1,2) 57 : REL(IX,IY+1) = REL(IX,IY+1) + DJAC*VKL2*XGENE(IX1,2) 58 : #*XGENE(IY1+8,2) 59 : REL(IX+1,IY+1) = REL(IX+1,IY+1) + DJAC*VKL2*XGENE(IX1+8,2)* 60 : #XGENE(IY1+8,2) 61 : REL(IY,IX) = REL(IX,IY) 62 : REL(IY,IX+1) = REL(IX+1,IY) 63 : REL(IY+1,IX) = REL(IX,IY+1) 64 : REL(IY+1,IX+1) = REL(IX+1,IY+1) 65 : 402 CONTINUE 66 : IY1=0 67 : 302 CONTINUE 68 : C 69 : C TERMES EN (UY,RX )*(UY,RX ) 70 : C 71 : IX1=0 72 : IY1=0 73 : DO 303 IX=5,LRE ,NBDL 74 : IX1=IX1 + 1 75 : DO 403 IY=5,IX ,NBDL 76 : IY1=IY1 + 1 77 : REL(IX,IY) = REL(IX,IY) + DJAC*VKL2*XGENE(IX1,2)*XGENE(IY1,2) 78 : REL(IX+1,IY) = REL(IX+1,IY) + DJAC*VKL2*XGENE(IX1+8,2) 79 : #*XGENE(IY1,2)*(-1.D0) 80 : REL(IX,IY+1) = REL(IX,IY+1) + DJAC*VKL2*XGENE(IX1,2) 81 : #*XGENE(IY1+8,2)*(-1.D0) 82 : REL(IX+1,IY+1) = REL(IX+1,IY+1) + DJAC*VKL2*XGENE(IX1+8,2)* 83 : #XGENE(IY1+8,2) 84 : REL(IY,IX) = REL(IX,IY) 85 : REL(IY,IX+1) = REL(IX+1,IY) 86 : REL(IY+1,IX) = REL(IX,IY+1) 87 : REL(IY+1,IX+1) = REL(IX+1,IY+1) 88 : 403 CONTINUE 89 : IY1=0 90 : 303 CONTINUE 91 : C 92 : C ENDIF 93 : RETURN 94 : END 95 : 96 :
© Cast3M 2003 - All rights reserved.
Disclaimer