Télécharger dychad.eso

Retour à la liste

Numérotation des lignes :

  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