1 : C BSIG SOURCE CHAT 05/01/12 21:42:58 5004 2 : 3 : C======================================================================= 4 : C= B S I G = 5 : C= ------- = 6 : C= = 7 : C= Fonction : = 8 : C= ---------- = 9 : C= Calcul de la contribution d'un point de Gauss aux forces internes = 10 : C= B*Sigma d'un element fini = 11 : C= = 12 : C= Parametres : (E)=Entree (S)=Sortie = 13 : C= ------------ = 14 : C= BGENE (E) Matrice de gradient (B) calculee au point de Gauss = 15 : C= XSTRS (E) Valeurs des contraintes au point de Gauss = 16 : C= NSTRS (E) Nombre de composantes de contraintes/deformations = 17 : C= LRE (E) Nombre de DDL de l'element fini = 18 : C= DJAC (S) Jacobien*poids au point de Gauss etudie = 19 : C= FORCE (E/S) Forces internes aux noeuds de l'element = 20 : C======================================================================= 21 : 22 : SUBROUTINE BSIG (BGENE,XSTRS,NSTRS,LRE,DJAC,XFORC) 23 : 24 : IMPLICIT INTEGER(I-N) 25 : IMPLICIT REAL*8 (A-H,O-Z) 26 : 27 : -INC CCREEL 28 : 29 : DIMENSION BGENE(NSTRS,*),XSTRS(*),XFORC(*) 30 : 31 : DO i=1,LRE 32 : cc=XZero 33 : DO j=1,NSTRS 34 : cc=cc+BGENE(j,i)*XSTRS(j) 35 : ENDDO 36 : XFORC(i)=cc*DJAC+XFORC(i) 37 : ENDDO 38 : 39 : RETURN 40 : END 41 : 42 : 43 :
© Cast3M 2003 - All rights reserved.
Disclaimer