Télécharger la plaquette Cast3M


fr en

Télécharger jo3loc.eso

Retour à la liste

Numérotation des lignes : 

jo3loc

  1. C JO3LOC    SOURCE    CHAT      05/01/13    00:49:44     5004
  2.       SUBROUTINE JO3LOC(XE,SHPTOT,IGAU,NBNO,BPSS)
  3. C=======================================================================
  4. C
  5. C    -CALCUL DE LA MATRICE DE PASSAGE BPSS
  6. C          ROUTINE FORTRAN PUR
  7. C    Pierre Pegon (ISPRA) OCTOBER 93
  8. C=======================================================================
  9. C  INPUT
  10. C     XE     = COORDONNEES  DE L ELEMENT
  11. C     SHPTOT = FONCTIONS DE FORME
  12. C            = SHPTOT(1,...) = FONCTIONS DE FORME
  13. C            = SHPTOT(2,...) = DERIVEE PAR RAPPORT A QSI
  14. C            = SHPTOT(3,...) = DERIVEE PAR RAPPORT A ETA
  15. C     IGAU   = NUMERO DU POINT DE GAUSS
  16. C     NBNO   = NOMBRE DE NOEUDS DE L'ELEMENT
  17. C  OUTPUT
  18. C     BPSS   = MATRICE DE PASAGE REPERE GLOBAL/REPERE LOCAL
  19. C
  20. C=======================================================================
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8(A-H,O-Z)
  23.       DIMENSION XE(3,*),BPSS(3,3),SHPTOT(6,NBNO,*)
  24.       DIMENSION XD(2,4)
  25.       DIMENSION S1(2),SN(2)
  26. C
  27.       NOQUAL = 0
  28. C
  29. C----------     CALCUL DE LA MATRICE DE  PASSAGE
  30. C
  31.       DO 6 I=1,2
  32.          S1(I)=0.0D0
  33.          SN(I)=0.0D0
  34.     6 CONTINUE
  35. C
  36.       DO 7 I=1,NBNO
  37. C
  38. C-------------------TANGENTE AU POINT DE GAUSS IGAU SELON QSI
  39. C
  40.          S1(1) = S1(1) + ( SHPTOT(2,I,IGAU)*XE(1,I) )
  41.          S1(2) = S1(2) + ( SHPTOT(2,I,IGAU)*XE(2,I) )
  42.     7  CONTINUE
  43. C
  44. C-------------------NORMALISATION DE S1
  45. C
  46.        XNORME = SQRT((S1(1)**2) + (S1(2)**2))
  47.        S1(1) = S1(1) / XNORME
  48.        S1(2) = S1(2) / XNORME
  49. C
  50. C-------------------NORMALE AU JOINT (PAR ROTATION DE 90 DEGRES)
  51. C
  52.        SN(1) =-S1(2)
  53.        SN(2) = S1(1)
  54. C
  55. C-------------------STOCKAGE DE LA MATRICE DE PASSAGE
  56. C
  57.        DO 10 I=1,2
  58.           BPSS(1,I) = S1(I)
  59.           BPSS(2,I) = SN(I)
  60.    10  CONTINUE
  61. C
  62.       RETURN
  63.       END
  64.  
  65.  
  66.  

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