Numérotation des lignes :

C BBGRDK    SOURCE    CHAT      05/01/12    21:36:32     5004      SUBROUTINE BBGRDK(QSI,ETA,GEOM,BB)C |====================================================================|C |                                                                    |C |   EVALUATION DE LA MATRICE BB(2,9) AU POINT QSI ETA                |C |   POUR LES TERMS DE GRADIAN DE L'ELEMENT DKT                       |C |     DUXDZ  DUYDZ   DUZDX  DUZDY        MAIS                        |C |     DUXDZ=-DUZDX=+BETA(X)=BB(1,I)   (I=1,2...9)                    |C |     DUYDZ=-DUZDY=+BETA(Y)=BB(2,I)   (I=1,2...9)                    |C |                                                                    |C |   CODE SUO X.Z.  23/01/1987                                        |C |====================================================================|      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)      DIMENSION GEOM(*),BB(2,*)C      DATA TRDE/1.5D0/,DEUX/2.D0/,QUATRE/4.D0/,UN/1.D0/C      A4=GEOM(1)      A5=GEOM(2)      A6=GEOM(3)      B4=GEOM(4)      B5=GEOM(5)      B6=GEOM(6)      C4= GEOM(7)      C5= GEOM(8)      C6= GEOM(9)      D4=GEOM(10)      D5=GEOM(11)      D6=GEOM(12)      E4=GEOM(13)      E5=GEOM(14)      E6=GEOM(15)      AUX=UN-QSI-ETACC       FONCTIONS DE FORME DE TRIANGLE A 6 NOEUDSC        SN1=-AUX*(1-DEUX*AUX)        SN2=-QSI*(1-DEUX*QSI)        SN3=-ETA*(1-DEUX*ETA)        SN4=QUATRE*QSI*ETA        SN5=QUATRE*ETA*AUX        SN6=QUATRE*QSI*AUXCC       FONCTIONS DE FORME DE TRIANGLE DKTC        BB(1,1)= TRDE*(A6*SN6-A5*SN5)        BB(1,2)= B5*SN5+B6*SN6        BB(1,3)= SN1-C5*SN5-C6*SN6        BB(1,4)= TRDE*(A4*SN4-A6*SN6)        BB(1,5)= B6*SN6+B4*SN4        BB(1,6)= SN2-C6*SN6-C4*SN4        BB(1,7)= TRDE*(A5*SN5-A4*SN4)        BB(1,8)= B4*SN4+B5*SN5        BB(1,9)= SN3-C4*SN4-C5*SN5        BB(2,1)= TRDE*(D6*SN6-D5*SN5)        BB(2,2)= -SN1+E5*SN5+E6*SN6        BB(2,3)= -B5*SN5-B6*SN6        BB(2,4)= TRDE*(D4*SN4-D6*SN6)        BB(2,5)= -SN2+E6*SN6+E4*SN4        BB(2,6)= -B6*SN6-B4*SN4        BB(2,7)= TRDE*(D5*SN5-D4*SN4)        BB(2,8)= -SN3+E4*SN4+E5*SN5        BB(2,9)= -B4*SN4-B5*SN5CC      RETURN      END

© Cast3M 2003 - Tous droits réservés.
Mentions légales