1 : C BBGRDK SOURCE CHAT 05/01/12 21:36:32 5004 2 : SUBROUTINE BBGRDK(QSI,ETA,GEOM,BB) 3 : C |====================================================================| 4 : C | | 5 : C | EVALUATION DE LA MATRICE BB(2,9) AU POINT QSI ETA | 6 : C | POUR LES TERMS DE GRADIAN DE L'ELEMENT DKT | 7 : C | DUXDZ DUYDZ DUZDX DUZDY MAIS | 8 : C | DUXDZ=-DUZDX=+BETA(X)=BB(1,I) (I=1,2...9) | 9 : C | DUYDZ=-DUZDY=+BETA(Y)=BB(2,I) (I=1,2...9) | 10 : C | | 11 : C | CODE SUO X.Z. 23/01/1987 | 12 : C |====================================================================| 13 : IMPLICIT INTEGER(I-N) 14 : IMPLICIT REAL*8 (A-H,O-Z) 15 : DIMENSION GEOM(*),BB(2,*) 16 : C 17 : DATA TRDE/1.5D0/,DEUX/2.D0/,QUATRE/4.D0/,UN/1.D0/ 18 : C 19 : A4=GEOM(1) 20 : A5=GEOM(2) 21 : A6=GEOM(3) 22 : B4=GEOM(4) 23 : B5=GEOM(5) 24 : B6=GEOM(6) 25 : C4= GEOM(7) 26 : C5= GEOM(8) 27 : C6= GEOM(9) 28 : D4=GEOM(10) 29 : D5=GEOM(11) 30 : D6=GEOM(12) 31 : E4=GEOM(13) 32 : E5=GEOM(14) 33 : E6=GEOM(15) 34 : AUX=UN-QSI-ETA 35 : C 36 : C FONCTIONS DE FORME DE TRIANGLE A 6 NOEUDS 37 : C 38 : SN1=-AUX*(1-DEUX*AUX) 39 : SN2=-QSI*(1-DEUX*QSI) 40 : SN3=-ETA*(1-DEUX*ETA) 41 : SN4=QUATRE*QSI*ETA 42 : SN5=QUATRE*ETA*AUX 43 : SN6=QUATRE*QSI*AUX 44 : C 45 : C FONCTIONS DE FORME DE TRIANGLE DKT 46 : C 47 : BB(1,1)= TRDE*(A6*SN6-A5*SN5) 48 : BB(1,2)= B5*SN5+B6*SN6 49 : BB(1,3)= SN1-C5*SN5-C6*SN6 50 : BB(1,4)= TRDE*(A4*SN4-A6*SN6) 51 : BB(1,5)= B6*SN6+B4*SN4 52 : BB(1,6)= SN2-C6*SN6-C4*SN4 53 : BB(1,7)= TRDE*(A5*SN5-A4*SN4) 54 : BB(1,8)= B4*SN4+B5*SN5 55 : BB(1,9)= SN3-C4*SN4-C5*SN5 56 : BB(2,1)= TRDE*(D6*SN6-D5*SN5) 57 : BB(2,2)= -SN1+E5*SN5+E6*SN6 58 : BB(2,3)= -B5*SN5-B6*SN6 59 : BB(2,4)= TRDE*(D4*SN4-D6*SN6) 60 : BB(2,5)= -SN2+E6*SN6+E4*SN4 61 : BB(2,6)= -B6*SN6-B4*SN4 62 : BB(2,7)= TRDE*(D5*SN5-D4*SN4) 63 : BB(2,8)= -SN3+E4*SN4+E5*SN5 64 : BB(2,9)= -B4*SN4-B5*SN5 65 : C 66 : C 67 : RETURN 68 : END 69 :
© Cast3M 2003 - All rights reserved.
Disclaimer