Télécharger la plaquette Cast3M


fr en

Télécharger vhynoe.eso

Retour à la liste

Numérotation des lignes : 

vhynoe

  1. C VHYNOE    SOURCE    CHAT      05/01/13    04:07:09     5004
  2.       SUBROUTINE VHYNOE(UELEM,IZNOEU,IZUN,NDIM,XY,ITP,IEL)
  3. C
  4. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  5. C
  6. C     CALCUL DE LA VITESSE REELE POUR UN NOEUD DU MAILLAGE EN
  7. C     FORMULATION HYBRIDE
  8. C
  9. C     ENTREES
  10. C          IZNOEU POINTEUR
  11. C              XELE(NDIM,NOEL) COORDONNES DES NOEUDS DE L ELEMENT CONSIDERE
  12. C          IZUN POINTEUR FLUX AUX CENTRES DES FACES DE L ELEMENT
  13. C          ITP             TYPE DE L'ELEMENT CONSIDERE
  14. C          XY          coordonnees DU NOEUD POUR LEQUEL ON CALCULERA
  15. C                         LA VITESSE
  16. C     SORTIE
  17. C          UELEM VITESSE REELLE AU NOEUD INOELO
  18. C
  19. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  20. C
  21. C
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8 (A-H,O-Z)
  24. C
  25. C
  26. C
  27.       SEGMENT  IZNOEU
  28.             REAL*8 XELE(IDIM,NOEL)
  29.             INTEGER NOEGLO(NOEL)
  30.       ENDSEGMENT
  31.       SEGMENT IZUN
  32.            REAL*8 UN(I1,I2,I3)
  33.       ENDSEGMENT
  34.       SEGMENT HYTRAV
  35.             REAL*8 XEL(3,NOEL),SHP(6,NOEL),SHY(12,NOEL)
  36.             REAL*8 ZJAC(NDIM,NDIM)
  37.       ENDSEGMENT
  38. C
  39. C
  40.       DIMENSION XY(3),UELEM(*),VECT(3)
  41. C
  42. C
  43. C                     FORMULATION HYBRIDE ON VA CALCULER LA VITESSE
  44. C                   REELLE AU NOEUD DU MAILLAGE
  45.            CALL INITD(UELEM,NDIM,0.D0)
  46.            NOEL=XELE(/2)
  47.            SEGINI HYTRAV
  48.            CALL SHAPE(XY(1),XY(2),XY(3),ITP,SHP,IRET)
  49. C          write(6,*)'shp ',(shp(1,j),j=1,noel)
  50.            DO 40 J=1,NOEL
  51.            DO 35 I=1,NDIM
  52.            XEL(I,J)=XELE(I,J)
  53.    35      CONTINUE
  54.    40      CONTINUE
  55.            CALL MATJAC(XEL,SHP,NDIM,NOEL,ZJAC,DJAC)
  56. C          write(6,*)' zjac',((zjac(i,j),i=1,ndim),j=1,ndim)
  57.            NTYP=NUMHYB(ITP)
  58.            CALL HYSHP(XY(1),XY(2),XY(3),NTYP,NOEL,12,SHY,IRET)
  59. C           write(6,*)' shy ',((shy(i,j),i=1,ndim),j=1,noel)
  60.            I2=UN(/2)
  61.            DO 50 IDDL=1,I2
  62.                 CALL INITD(VECT,NDIM,0.D0)
  63.                 DO 60 J=1,NDIM
  64.                      DO 55 I=1,NDIM
  65.                         VECT(I)=VECT(I)+ZJAC(I,J)*SHY(J,IDDL)
  66.    55                CONTINUE
  67.    60           CONTINUE
  68.                 COEF=UN(1,IDDL,IEL)/DJAC
  69. C       write(6,*)' flux ',iddl,iel,UN(1,IDDL,IEL),coef,djac
  70.                      DO 65 I=1,NDIM
  71.                           UELEM(I)=UELEM(I)+COEF*VECT(I)
  72.    65                CONTINUE
  73.    50           CONTINUE
  74.            SEGSUP HYTRAV
  75.       RETURN
  76.       END
  77.  
  78.  

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