Télécharger jo3loc.eso

Retour à la liste

Numérotation des lignes :

  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