Télécharger la plaquette Cast3M


fr en

Télécharger dyche3.eso

Retour à la liste

Numérotation des lignes : 

dyche3

  1. C DYCHE3    SOURCE    CHAT      05/01/12    23:04:09     5004
  2.       SUBROUTINE DYCHE3(XDEP,IDIM,IGP,XPALB,NLIAB,I,INTER
  3.      &                 ,XFN,XFT,XPUS,iannul)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. *                                                                    *
  8. *     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *
  9. *     ________________________________________________               *
  10. *                                                                    *
  11. *     Calcul de la force de choc pour un choc {l{mentaire.           *
  12. *                                                                    *
  13. *     Param}tres:                                                    *
  14. *                                                                    *
  15. * e   XDEP    Valeur du d{placement sur la normale                   *
  16. * e   IGP     Indice de glissement au pas pr{c{dent                  *
  17. * es  XPALB   Tableau contenant les param}tres de la liaison         *
  18. * e   NLIAB   Nombre de liaison sur la base B                        *
  19. * e   I       Num{ro de la liaison                                   *
  20. * e   INTER   Indice de liaison interieure / exterieure              *
  21. *  s  XFN     Valeur de la force de choc suivant la normale.         *
  22. *  s  XFT     Module de la force de choc projetee sur le plan tangent*
  23. *  s  XPUS    Valeur de la puissance d'usure instantanee             *
  24. *                                                                    *
  25. *                                                                    *
  26. *     Auteur, date de cr{ation:                                      *
  27. *                                                                    *
  28. *     Bertrand BEAUFILS, le 30 mai 1990.                             *
  29. *                                                                    *
  30. *--------------------------------------------------------------------*
  31.       DIMENSION XPALB(NLIAB,*)
  32.       XFN = 0.D0
  33.       XFT = 0.D0
  34.       XPUS = 0.D0
  35.       XJEU  = XPALB(I,2)
  36.       XAD = ABS(XDEP)
  37.       XAJ = ABS(XJEU)
  38. *
  39.       IF ( (((XAD.LE.XAJ).AND.(INTER.EQ.0)).OR.
  40.      &      ((XAD.GE.XAJ).AND.(INTER.EQ.1))).AND.
  41.      &      (iannul.EQ.0)) THEN
  42.          ID4 = 6 + 3*IDIM
  43.          ID5 = 6 + 4*IDIM
  44.          ID6 = 6 + 5*IDIM
  45.          ID8 = 6 + 7*IDIM
  46.          ID9 = 6 + 8*IDIM
  47.          XRAIN = XPALB(I,1)
  48.          XGLIS = XPALB(I,3)
  49.          XADHE = XPALB(I,4)
  50.          XRAIT = XPALB(I,5)
  51.          XAMOT = XPALB(I,6)
  52.          XVAL = XDEP - XJEU
  53.          XFN = -XRAIN * XVAL
  54.          CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID4,ID5,ID6,ID8,ID9,XGLIS,
  55.      &               XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  56.       ELSE
  57.          ID6 = 6 + 5*IDIM
  58.          IGP = -1
  59.          DO 10 ID=1,IDIM
  60.             XPALB(I,ID6+ID) = 0.D0
  61.  10         CONTINUE
  62. *        end do
  63.       ENDIF
  64. *
  65.       RETURN
  66.       END
  67.  
  68.  
  69.  
  70.  

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