1 : C BGFCQ4 SOURCE CHAT 06/03/29 21:15:48 5360 2 : SUBROUTINE BGFCQ4(IGAU,XEL,SHPTOT,SHP,BGF,DJAC,IARR) 3 : C |====================================================================| 4 : C | | 5 : C | CALCULE LA MATRICE DE GRADIAN EN FLEXION BGF POUR COQ4 6 : C | CODE SUO X.Z. MARS 87 7 : C |INPUT | 8 : C | IGAU=NUMERO DU POINT DE GAUSS | 9 : C | NBNO=NOMBRE DE NOEUDS | 10 : C | LRE =NOMBRE DE COLONNES DE LA MATRICE BGR | 11 : C | XEL =COORDONNEES LOCALES DE L ELEMENT | 12 : C | SHPTOT(6,NBNO,NBGAU)=FONCTIONS DE FORMES ET DERIVEES | 13 : C |ZONE DE TRAVAIL | 14 : C | SHP(6,NBNO)=TABLEAU DE TRAVAIL | 15 : C |OUTPUT | 16 : C | DJAC=JACOBIEN | 17 : C | BGF(9,LRE)=MATRICE BGF | 18 : C |====================================================================| 19 : IMPLICIT INTEGER(I-N) 20 : IMPLICIT REAL*8(A-H,O-Z) 21 : PARAMETER (LRE=24,NBNO=4) 22 : DIMENSION XEL(3,*),BGF(9,*),SHP(6,*),SHPTOT(6,NBNO,*) 23 : DATA ZER/0.D0/ 24 : IARR=0 25 : C 25 : CALL ZERO(BGF,9,LRE) 26 : C 27 : C 28 : DO 101 NP=1,NBNO 29 : SHP(1,NP)=SHPTOT(1,NP,IGAU) 30 : SHP(2,NP)=SHPTOT(2,NP,IGAU) 31 : SHP(3,NP)=SHPTOT(3,NP,IGAU) 32 : 101 CONTINUE 32 : CALL JACOBI(XEL,SHP,2,NBNO,DJAC) 33 : IF(DJAC.LE.ZER) IARR=1 34 : IF(DJAC.LE.ZER) RETURN 35 : C 36 : C ====================== CALCUL 37 : C 38 : K=1 39 : DO 102 NP=1,NBNO 40 : BGF(1,K+4)=SHP(2,NP) 41 : BGF(2,K+4)=SHP(3,NP) 42 : BGF(3,K+4)=SHP(1,NP) 43 : BGF(4,K+3)=-SHP(2,NP) 44 : BGF(5,K+3)=-SHP(3,NP) 45 : BGF(6,K+3)=-SHP(1,NP) 46 : K=K+6 47 : 102 CONTINUE 48 : RETURN 49 : END 50 : 51 :
© Cast3M 2003 - All rights reserved.
Disclaimer