Télécharger murloc.eso

Retour à la liste

Numérotation des lignes :

  1. C MURLOC SOURCE CHAT 05/01/13 01:55:49 5004
  2. C MURLOC SOURCE PP 27/04/04 21:25:46 4501
  3. SUBROUTINE MURLOC(XE,NBNO,LHOOK,LRE,BPSS,XH,BGENE)
  4. C=======================================================================
  5. C
  6. C -TEST DE PLANEITE DES FACES DE L'ELEMENT JOI2
  7. C -CALCUL DE LA MATRICE DE PASSAGE BPSS
  8. C -CALCUL DES COORDONNEES LOCALES XEL
  9. C ROUTINE FORTRAN PUR
  10. C CODE S. FELIX SEPTEMBRE 92
  11. C=======================================================================
  12. C INPUT
  13. C XE = COORDONNEES DE L ELEMENT
  14. C NBNO = NOMBRE DE NOEUDS DE L'ELEMENT
  15. C OUTPUT
  16. C BPSS = MATRICE DE PASAGE REPERE GLOBAL/REPERE LOCAL
  17. C XH = HAUTEUR DE L'ELEMENT
  18. C BGENE = MATRICE B
  19. C
  20. C=======================================================================
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. DIMENSION XE(3,*),BPSS(3,3),BGENE(LHOOK,LRE)
  24. DIMENSION S1(2),SN(2)
  25. PARAMETER(ZERO=0.D0, UN=1.D0)
  26. C
  27. C---------- CALCUL DE LA MATRICE DE PASSAGE
  28. C
  29. S1(1) = XE(1,3) - XE(1,1)
  30. S1(2) = XE(2,3) - XE(2,1)
  31. C
  32. XH = SQRT((S1(1)**2) + (S1(2)**2))
  33. S1(1) = S1(1) / XH
  34. S1(2) = S1(2) / XH
  35. C
  36. SN(1) =-S1(2)
  37. SN(2) = S1(1)
  38. C
  39. DO I=1,2
  40. BPSS(1,I) = S1(I)
  41. BPSS(2,I) = SN(I)
  42. ENDDO
  43. C
  44. C---------- CALCUL DE BGENE
  45. C
  46. DO I=1,LHOOK
  47. DO J=1,LRE
  48. BGENE(I,J)=ZERO
  49. ENDDO
  50. ENDDO
  51. C
  52. BGENE(1,1)=-BPSS(1,1)
  53. BGENE(1,2)=-BPSS(1,2)
  54. BGENE(1,7)= UN
  55. C
  56. BGENE(2,3)=-UN
  57. BGENE(2,8)= UN
  58. C
  59. BGENE(3,1)=-BPSS(2,1)
  60. BGENE(3,2)=-BPSS(2,2)
  61. BGENE(3,4)= BPSS(2,1)
  62. BGENE(3,5)= BPSS(2,2)
  63. BGENE(3,8)= XH
  64. C
  65. BGENE(4,4)= BPSS(1,1)
  66. BGENE(4,5)= BPSS(1,2)
  67. BGENE(4,7)=-UN
  68. C
  69. BGENE(5,6)= UN
  70. BGENE(5,8)=-UN
  71. C
  72. RETURN
  73. END
  74.  
  75.  
  76.  
  77.  

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