Télécharger la plaquette Cast3M


fr en

Télécharger trjcou.eso

Retour à la liste

Numérotation des lignes : 

trjcou

  1. C TRJCOU    SOURCE    CHAT      05/01/13    03:49:24     5004
  2.       SUBROUTINE TRJCOU(UE,IZCOU,IZFACE,IEL,JEL,NUCENT)
  3. C
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8 (A-H,O-Z)
  6. C
  7. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  8. C
  9. C     CALCUL DE DTCO(JEL) POUR AVOIR UN NOMBRE DE COURANT DE 1 DANS IEL
  10. C     UE  VITESSES MAX POUR L'ELEMENT JEL
  11. C     IZCOU POINTEUR DU SEGMENT CONTENANT LE NOMBRE DE COURANT
  12. C     IZFACE POINTEUR DU MAILLAGE CONTENANT LES CENTRES DES FACES
  13. C     IEL NUMERO DE L ELEMENT CONSIDERE DANS LE SOUS MAILLAGE
  14. C     JEL NUMERO DE L ELEMENT CONSIDERE DANS LE MAILLAGE COMPLET
  15. C     NUCENT NUMERO DU NOEUD CENTRE DE L ELEMENT
  16. C
  17. C
  18. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  19. C
  20. C
  21.  
  22. -INC PPARAM
  23. -INC CCOPTIO
  24. -INC SMCOORD
  25. -INC SMELEME
  26.       POINTEUR IZFACE.MELEME
  27.       SEGMENT IZUMAX
  28.            REAL*8 UMAX(NBREL)
  29.       ENDSEGMENT
  30. C
  31.       SEGMENT  IZCOU
  32.            REAL*8 DTCO(NEL),COU
  33.       ENDSEGMENT
  34. C
  35. C
  36.       NF=IZFACE.NUM(/1)
  37. C        ON VA DETERMINER DIAM : LA PLUS PETITE DISTANCE ENTRE LE
  38. C        POINT CENTRE DE L ELEMENT ET LE MILIEU DES FACES
  39.       IPCENT=(NUCENT-1)*(IDIM+1)
  40.       NF1=(IZFACE.NUM(1,IEL)-1)*(IDIM+1)
  41.       VINT=0.D0
  42.       DO 5 I=1,IDIM
  43.         VINT=VINT+(XCOOR(IPCENT+I)-XCOOR(NF1+I))**2
  44.     5 CONTINUE
  45.       DIAM=VINT
  46.       DO 20 J=2,NF
  47.            NF1=(IZFACE.NUM(J,IEL)-1)*(IDIM+1)
  48.            VINT=0.D0
  49.            DO 15 I=1,IDIM
  50.              VINT=VINT+(XCOOR(IPCENT+I)-XCOOR(NF1+I))**2
  51.    15      CONTINUE
  52.            IF(VINT.LT.DIAM)DIAM=VINT
  53.    20 CONTINUE
  54.       DIAM=SQRT(DIAM)*2.D0
  55. C
  56. C
  57. C     UE=UMAX(JEL)
  58. C     write(6,*)' DIAM UE ',DIAM,UE,JEL
  59.       DTCO(JEL)=DIAM/UE
  60.       RETURN
  61.       END
  62.  
  63.  
  64.  

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