Télécharger dychgl.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHGL SOURCE CHAT 05/01/12 23:06:07 5004
  2. SUBROUTINE DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,XGLIS,XFN,XFT,
  3. & 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 de glissement *
  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 force tangentielle dans XPALB *
  20. * e XFN Valeur de la force de choc suivant la normale. *
  21. * e XGLIS Valeur du coefficient d'adh{rence *
  22. * s IGP Indice de glissement au pas pr{c{dent *
  23. * s XFT Module de la force de choc projet{e sur le plan tangent*
  24. * s XPUS Valeur de la puissance d'usure instantan{e *
  25. * *
  26. * *
  27. * Auteur, date de cr{ation: *
  28. * *
  29. * Bertrand BEAUFILS, le 31 mai 1990. *
  30. * *
  31. *--------------------------------------------------------------------*
  32. DIMENSION XPALB(NLIAB,*)
  33. XFT = - XGLIS * ABS(XFN)
  34. IGP = 1
  35. *
  36. PSVITT = 0.D0
  37. DO 10 ID = 1,IDIM
  38. PSVITT = PSVITT + XPALB(I,I0+ID) * XPALB(I,I0+ID)
  39. 10 CONTINUE
  40. * end do
  41. VITT = SQRT(PSVITT)
  42. XPUS = ABS(XFN * VITT)
  43. *
  44. IF (VITT.GT.1.D-20) THEN
  45. DO 11 ID = 1,IDIM
  46. XPALB(I,I1+ID) = XFT * XPALB(I,I0+ID) / VITT
  47. 11 CONTINUE
  48. * end do
  49. ELSE
  50. DO 12 ID = 1,IDIM
  51. XPALB(I,I1+ID) = 0.D0
  52. 12 CONTINUE
  53. * end do
  54. ENDIF
  55. *
  56. RETURN
  57. END
  58.  
  59.  

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