Télécharger la plaquette Cast3M


fr en

Télécharger bgfdkt.eso

Retour à la liste

Numérotation des lignes : 

bgfdkt

  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
  25.       CALL ZERO(BGF,9,LRE)
  26. C
  27.       IFR=IFOU+4
  28.       IF (ITEL.NE.28) RETURN
  29. C
  30. C  28 IEME   ELEMENT :   DKT
  31. C
  32.   28  CONTINUE
  33.       CALL GEOCST(XEL,GEOM)
  34.       CALL BBGFDK(XX(IGAU),YY(IGAU),GEOM,BB)
  35.       DJAC=GEOM(17)
  36.       K=2
  37.       KK=0
  38.       DO 327 NPOI=1,3
  39.       DO 527 IY=1,3
  40.       BGF(1,K+IY)=BB(1,IY+KK)
  41.       BGF(2,K+IY)=BB(2,IY+KK)
  42.       BGF(4,K+IY)=BB(3,IY+KK)
  43.       BGF(5,K+IY)=BB(4,IY+KK)
  44.   527 CONTINUE
  45.       KK=KK+3
  46.       K=K+6
  47.   327 CONTINUE
  48.       CALL GERADK(XEL,GEOM)
  49.       CALL BBGRDK(XX(IGAU),YY(IGAU),GEOM,CC)
  50.       K=2
  51.       KK=0
  52.       DO 310 NPOI=1,3
  53.       DO 510 IY=1,3
  54.       BGF(3,K+IY)=+CC(1,IY+KK)
  55.       BGF(6,K+IY)=+CC(2,IY+KK)
  56.   510 CONTINUE
  57.       KK=KK+3
  58.       K=K+6
  59.   310 CONTINUE
  60.       RETURN
  61.       END
  62.  
  63.  

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