Télécharger la plaquette Cast3M


fr en

Télécharger murloc.eso

Retour à la liste

Numérotation des lignes : 

murloc

  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