1 : C BDBST SOURCE CHAT 05/01/12 21:37:58 5004 2 : 3 : C======================================================================= 4 : C= B D B S T = 5 : C= --------- = 6 : C= = 7 : C= Fonction : = 8 : C= ---------- = 9 : C= Calcul de la contribution d'un point de Gauss a la matrice de = 10 : C= RIGIDITE elastique REL d'un element fini. = 11 : C= Contribution : Jacobien.transpose(BGENE)*DHOOK*BGENE = 12 : C= = 13 : C= Parametres : (E)=Entree (S)=Sortie = 14 : C= ------------ = 15 : C= BGENE (E) Matrice B = 16 : C= DJac (E) "Volume" associe au point de Gauss = 17 : C= DHOOK (E) Matrice de Hooke au point de Gauss = 18 : C= LRE (E) Nombre de DDL de l'element fini = 19 : C= NSTRS (E) Nombre de composantes de contraintes/deformations = 20 : C= REL (E/S) Matrice de rigidite de l'element fini = 21 : C= = 22 : C= Remarque : Actuellement, NSTRS doit etre inferieur ou egal a 10 ! = 23 : C======================================================================= 24 : 25 : SUBROUTINE BDBST(BGENE,DJAC,DHOOK,LRE,NSTRS,REL) 26 : 27 : IMPLICIT INTEGER(I-N) 28 : IMPLICIT REAL*8 (A-H,O-Z) 29 : 30 : -INC CCREEL 31 : 32 : DIMENSION DHOOK(NSTRS,*),REL(LRE,*),BGENE(NSTRS,*) 33 : DIMENSION WORK1(10),WORK2(10) 34 : 35 : C= Boucle sur les DDL 36 : DO IX=1,LRE 37 : C== Calcul de JACobien*Poids*B pour le DDL IX 38 : DO i=1,NSTRS 39 : WORK1(i)=DJAC*BGENE(i,IX) 40 : ENDDO 41 : C== Calcul de Vol.*B*D pour le DDL IX 42 : DO i=1,NSTRS 43 : CC=XZero 44 : DO j=1,NSTRS 45 : CC=CC+WORK1(j)*DHOOK(i,j) 46 : ENDDO 47 : WORK2(i)=CC 48 : ENDDO 49 : C== Calcul de Vol.Bt.D.B pour les DDL IX et IY 50 : DO IY=1,IX 51 : CC=XZero 52 : DO i=1,NSTRS 53 : CC=CC+BGENE(i,IY)*WORK2(i) 54 : ENDDO 55 : REL(IX,IY)=REL(IX,IY)+CC 56 : REL(IY,IX)=REL(IX,IY) 57 : ENDDO 58 : ENDDO 59 : 60 : RETURN 61 : END 62 : 63 : 64 :
© Cast3M 2003 - All rights reserved.
Disclaimer