1 : C BGFDKT SOURCE CHAT 05/01/12 21:39:27 5004 2 : SUBROUTINE BGFDKT(IGAU,ITEL,LRE,IFOU,XEL,BGF) 3 : C|=================================================================| 4 : C| | 5 : C| CALCULE LA MATRICE BGF DE GRADIENS DE FLEXION ( GRAFLEXI ) | 6 : C| ROUTINE FORTRAN PUR | 7 : C| CODE SUO X.Z. MARS 87 | 8 : C| INPUT | 9 : C| IGAU=NUMERO DU POINT DE GAUSS | 10 : C| ITEL=NUMERO DE L ELEMENT DANS NOMTP | 11 : C| LRE =NOMBRE DE COLONNES DE LA MATRICE B | 12 : C| IFOU=IFOUR DE CCOPTIO | 13 : C| XEL=COORDONNEES DE L ELEMENT | 14 : C| OUTPUT | 15 : C| BGF(9,LRE)=MATRICE BGF | 16 : C| DJAC=JACOBIEN DE L ELEMENT | 17 : C|=================================================================| 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : DIMENSION XEL(3,*),BGF(9,*),CC(2,9) 21 : DIMENSION BB(4,9),GEOM(20),XX(3),YY(3) 22 : DATA XX/.5D0,.0D0,.5D0/ 23 : DATA YY/.0D0,.5D0,.5D0/ 24 : C 24 : CALL ZERO(BGF,9,LRE) 25 : C 26 : IFR=IFOU+4 27 : IF (ITEL.NE.28) RETURN 28 : C 29 : C 28 IEME ELEMENT : DKT 30 : C 31 : 28 CONTINUE 31 : CALL GEOCST(XEL,GEOM) 31 : CALL BBGFDK(XX(IGAU),YY(IGAU),GEOM,BB) 32 : DJAC=GEOM(17) 33 : K=2 34 : KK=0 35 : DO 327 NPOI=1,3 36 : DO 527 IY=1,3 37 : BGF(1,K+IY)=BB(1,IY+KK) 38 : BGF(2,K+IY)=BB(2,IY+KK) 39 : BGF(4,K+IY)=BB(3,IY+KK) 40 : BGF(5,K+IY)=BB(4,IY+KK) 41 : 527 CONTINUE 42 : KK=KK+3 43 : K=K+6 44 : 327 CONTINUE 44 : CALL GERADK(XEL,GEOM) 44 : CALL BBGRDK(XX(IGAU),YY(IGAU),GEOM,CC) 45 : K=2 46 : KK=0 47 : DO 310 NPOI=1,3 48 : DO 510 IY=1,3 49 : BGF(3,K+IY)=+CC(1,IY+KK) 50 : BGF(6,K+IY)=+CC(2,IY+KK) 51 : 510 CONTINUE 52 : KK=KK+3 53 : K=K+6 54 : 310 CONTINUE 55 : RETURN 56 : END 57 :
© Cast3M 2003 - All rights reserved.
Disclaimer