Télécharger bgfdkt.eso

Retour à la liste

Numérotation des lignes :

  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