Télécharger la plaquette Cast3M


fr en

Télécharger dychad.eso

Retour à la liste

Numérotation des lignes : 

dychad

  1. C DYCHAD    SOURCE    CHAT      05/01/12    23:02:39     5004
  2.       SUBROUTINE DYCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XADHE,XRAIT,
  3.      &                  XAMOT,XFN,XFT,XPUS)
  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 d'adh{rence et test de sa validit{          *
  12. *                                                                    *
  13. *     Param}tres:                                                    *
  14. *                                                                    *
  15. * es  XPALB   Tableau contenant les param}tres de la liaison         *
  16. * e   NLIAB   Nombre de liaisons sur la base B                       *
  17. * e   I       Num{ro de la liaison                                   *
  18. * e   I0      Pointeur de la vitesse tangentielle dans XPALB         *
  19. * e   I1      Pointeur de la variation de d{placement par rapport a  *
  20. *             la position origine d'adh{rence, projet{e sur le plan  *
  21. *             tangent, dans XPALB                                    *
  22. * e   I2      Pointeur de la force tangentielle dans XPALB           *
  23. * e   I3      Pointeur de la force tangentielle @ la fin du glis-    *
  24. *             -sement dans XPALB                                     *
  25. * e   XFN     Valeur de la force de choc suivant la normale.         *
  26. * e   XADHE   Valeur du coefficient d'adh{rence                      *
  27. * e   XRAIT   Valeur de la raideur tangentielle                      *
  28. * e   XAMOT   Valeur de l'amortissement tangentiel                   *
  29. * es  IGP     Indice de glissement au pas pr{c{dent                  *
  30. *  s  XFT     Module de la force de choc projet{e sur le plan tangent*
  31. *  s  XPUS    Valeur de la puissance d'usure instantan{e             *
  32. *                                                                    *
  33. *                                                                    *
  34. *     Auteur, date de cr{ation:                                      *
  35. *                                                                    *
  36. *     Bertrand BEAUFILS, le 31 mai 1990.                             *
  37. *                                                                    *
  38. *--------------------------------------------------------------------*
  39.       DIMENSION XPALB(NLIAB,*)
  40. *  Calcul de la force tangentielle introduite par l'oscillateur
  41. *  {quivalent
  42.       PSXFT = 0.D0
  43.       DO 10 ID = 1,IDIM
  44.           XPALB(I,I2+ID) = XPALB(I,I3+ID) - XRAIT * XPALB(I,I1+ID)
  45.      &                   - XAMOT * XPALB(I,I0+ID)
  46.           PSXFT = PSXFT + XPALB(I,I2+ID) * XPALB(I,I2+ID)
  47.  10       CONTINUE
  48. *     end do
  49.       XFT = SQRT(PSXFT)
  50. *
  51.       IF (XFT.GT.( XADHE * ABS(XFN) ) ) THEN
  52. *  L'adh{rence ne tient pas
  53.          IGP = 1
  54.       ELSE
  55. *  Ca colle
  56.          XPUS = 0.D0
  57.          IGP = 0
  58.       ENDIF
  59. *
  60.       RETURN
  61.       END
  62.  
  63.  

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