Télécharger la plaquette Cast3M


fr en

Télécharger bjo3.eso

Retour à la liste

Numérotation des lignes : 

bjo3

  1. C BJO3      SOURCE    CHAT      05/01/12    21:40:29     5004
  2.       SUBROUTINE BJO3(IGAU,MFR,IFOU,NIFOU,XEL,BPSS,SHPTOT,SHP,
  3.      .                                  BGENE,DJAC,IRRT)
  4. C=======================================================================
  5. C
  6. C     CALCUL DE LA MATRICE B = TETA * ( N , -N )
  7. C     ET DU JACOBIEN EN IGAU POUR UN JOINT 2D
  8. C          ROUTINE FORTRAN PUR
  9. C     Pierre Pegon (ISPRA) 28/10/93
  10. C=======================================================================
  11. C  INPUT
  12. C     IGAU  = NUMERO DU POINT DE GAUSS
  13. C     MFR   = NUMERO DE LA FORMULATION
  14. C     XEL   = COORDONNEES DES NOEUDS DE L'ELEMENT
  15. C     BPSS  = MATRICE DE PASSAGE
  16. C             BPSS(,) = AXE S1
  17. C             BPSS(,) = AXE SN
  18. C     SHPTOT= FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE DE REFERENCE
  19. C             SHPTOT(1, ) = FONCTION DE FORME
  20. C             SHPTOT(2, ) = DERIVEES PAR RAPPORT A QSI
  21. C             SHPTOT(3, ) = DERIVEES PAR RAPPORT A ETA
  22. C  OUTPUT
  23. C     SHP   = FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE GEOMETRIQUE
  24. C             SHP(1, ) = FONCTION DE FORME
  25. C             SHP(2, ) = DERIVEES PAR RAPPORT A X LOCAL
  26. C             SHP(3, ) = DERIVEES PAR RAPPORT A Y LOCAL
  27. C     DJAC  = JACOBIEN AU POINT D'INTEGRATION IGAU
  28. C     BGENE = MATRICE B AU POINT D'INTEGRATION IGAU
  29. C     IRRT  = DIFFERENT DE ZERO SI ERREUR
  30. C  CONSTANTES
  31. C     NBNO  = NOMBRE DE NOEUDS
  32. C     LRE   = NOMBRE DE COLONNES DE LA MATRICE B
  33. C     NST   = NOMBRE DE COMPOSANTES DE CONTRAINTES
  34. C
  35. C  REMARQUE : CE SOUS-PROGRAMME EST VALABLE EN CONTRAINTES PLANES,
  36. C             DEFORMATIONS PLANES ET AXISYMETRIQUES. LA MATRICE
  37. C                     BGENE = TETA ( Ntop , - Nbot )
  38. C             EST UNE MATRICE (2x12). EN EFFET, TETA EST UNE MATRICE
  39. C             (2x2), ET ( Ntop , - Nbot ) EST UNE MATRICE (2x12).
  40. C             LA MATRICE TETA EST ICI APPELEE BPSS.
  41. C
  42. C=======================================================================
  43.       IMPLICIT INTEGER(I-N)
  44.       IMPLICIT REAL*8(A-H,O-Z)
  45.       PARAMETER (LRE=12,NST=2,NBNO=6)
  46.       DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)
  47.       DIMENSION BPSS(3,3)
  48. C
  49.       IRRT = 0
  50. C
  51. C     JACOBIEN
  52. C
  53.       DO 1 I=1,NBNO
  54.        SHP(1,I) = SHPTOT(1,I,IGAU)
  55.        SHP(2,I) = SHPTOT(2,I,IGAU)
  56.     1 CONTINUE
  57. C
  58.       NBNONN=NBNO/2
  59.       CALL DEVOLU(XEL,SHP,MFR,NBNONN,IFOU,NIFOU,86,1.D0,RR,DJAC)
  60. C
  61.       IF (DJAC.EQ.0.0D0) THEN
  62.          IRRT = 2
  63.       ENDIF
  64. C
  65. C     MATRICE B ATTENTION AU RESPECT DE LA NUMEROTATION LOCALE
  66. C     DES NOEUDS !!!
  67. C
  68.       CALL ZERO(BGENE,NST,LRE)
  69.       DO 2 I=1,2
  70.        DO 3 J=1,3
  71.         DO 4 K=1,2
  72.          L=2*(J-1)+K
  73. C CE N'EST PAS M=L+6 ....
  74.          M=L+2*(2*(3-J)+1)
  75.          BGENE(I,L)=BPSS(I,K)*SHP(1,J)
  76.          BGENE(I,M)=-BGENE(I,L)
  77.     4   CONTINUE
  78.     3  CONTINUE
  79.     2 CONTINUE
  80. C
  81.       RETURN
  82.       END
  83.  
  84.  
  85.  

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