1 : C BJT3 SOURCE CHAT 05/01/12 21:40:46 5004 2 : SUBROUTINE BJT3(IGAU,MFR,IFOU,NIFOU,XEL,BPSS,SHPTOT,SHP, 3 : . BGENE,DJAC,IRRT) 4 : C======================================================================= 5 : C 6 : C CALCUL DE LA MATRICE B = TETA * ( N , -N ) 7 : C ET DU JACOBIEN EN IGAU 8 : C ROUTINE FORTRAN PUR 9 : C DERIVEE DE LA ROUTINE BJO4 PAR S. FELIX 10 : C======================================================================= 11 : C INPUT 12 : C IGAU = NUMERO DU POINT DE GAUSS 13 : C XEL = COORDONNEES DES NOEUDS DE L'ELEMENT 14 : C BPSS = MATRICE DE PASSAGE 15 : C BPSS(,) = AXE S1 16 : C BPSS(,) = AXE S2 17 : C BPSS(,) = AXE SN 18 : C SHPTOT= FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE DE REFERENCE 19 : C SHPTOT(1, ) = FONCTION DE FORME 20 : C SHPTOT(2, ) = DERIVEES PAR RAPPORT A QSI 21 : C SHPTOT(3, ) = DERIVEES PAR RAPPORT A ETA 22 : C OUTPUT 23 : C SHP = FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE GEOMETRIQUE 24 : C SHP(1, ) = FONCTION DE FORME 25 : C SHP(2, ) = DERIVEES PAR RAPPORT A X LOCAL 26 : C SHP(3, ) = DERIVEES PAR RAPPORT A Y LOCAL 27 : C DJAC = JACOBIEN AU POINT D'INTEGRATION IGAU 28 : C BGENE = MATRICE B AU POINT D'INTEGRATION IGAU 29 : C IRRT = DIFFERENT DE ZERO SI ERREUR 30 : C CONSTANTES 31 : C NBNO = NOMBRE DE NOEUDS 32 : C LRE = NOMBRE DE COLONNES DE LA MATRICE B 33 : C NST = NOMBRE DE COMPOSANTES DE CONTRAINTES 34 : C 35 : C REMARQUE : ATTENTION : DANS LES CAS CONTRAINTES PLANES ET 36 : C AXISYMETRIQUE, LA MATRICE BPSS ( RESPECTIVEMENT BGENE ) 37 : C N'A PLUS LA DIMENSION (3X3) ( RESPECTIVEMENT (3X18) ). 38 : C TROUVER LEURS DIMENSIONS CORRECTES ET MODIFIER LES 39 : C PARAMETRES LRE, NST, NBNO. CEUX CI CORRESPONDRONT A 40 : C CEUX DU CAS BIDIMENSIONNELS. 41 : C 42 : C======================================================================= 43 : IMPLICIT INTEGER(I-N) 44 : IMPLICIT REAL*8(A-H,O-Z) 45 : PARAMETER (LRE=18,NST=3,NBNO=6) 46 : DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*) 47 : DIMENSION BPSS(3,3) 48 : C 49 : IRRT = 0 50 : C 51 : C MATRICE JACOBIENNE 52 : C 53 : DO 1 I=1,NBNO 54 : SHP(1,I) = SHPTOT(1,I,IGAU) 55 : SHP(2,I) = SHPTOT(2,I,IGAU) 56 : SHP(3,I) = SHPTOT(3,I,IGAU) 57 : 1 CONTINUE 58 : C 59 : C !!!!!! ATTENTION : IL FAUT CALCULER LE JACOBIEN AVEC NBNO=3 !!!!!! 60 : C IL NE FAUT SURTOUT PAS METTRE NBNO=6 CAR CA FAUSSE ALORS 61 : C LES RESULTATS ... NBNO DOIT ETRE EGAL A 6 CAR IL YA 6 NOEUDS 62 : C 63 : NBNONN=NBNO/2 63 : CALL DEVOLU(XEL,SHP,MFR,NBNONN,IFOU,NIFOU,2,1.D0,RR,DJAC) 64 : IF (DJAC.LT.0.0D0) THEN 65 : IRRT = 1 66 : ELSE IF (DJAC.EQ.0.0D0) THEN 67 : IRRT = 2 68 : ENDIF 69 : C 70 : C MATRICE B 71 : C 71 : CALL ZERO(BGENE,NST,LRE) 72 : DO 2 I=1,3 73 : DO 3 J=1,3 74 : DO 4 K=1,3 75 : L=3*(J-1)+K 76 : M=L+9 77 : BGENE(I,L)=BPSS(I,K)*SHP(1,J) 78 : BGENE(I,M)=-BGENE(I,L) 79 : 4 CONTINUE 80 : 3 CONTINUE 81 : 2 CONTINUE 82 : C 83 : RETURN 84 : END 85 : 86 : 87 : 88 :
© Cast3M 2003 - All rights reserved.
Disclaimer