Télécharger dycglr.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCGLR SOURCE CHAT 05/01/12 23:00:14 5004
  2. C DYCGLR SOURCE BEAU 90/07/26 21:16:02
  3. SUBROUTINE DYCGLR(IDIM,EPSI,XPALB,NLIAB,I,I0,I1,XGLIS,XFN,XFT,
  4. & XPUS)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. * *
  9. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  10. * ________________________________________________ *
  11. * *
  12. * Calcul de la force de glissement *
  13. * *
  14. * Param}tres: *
  15. * e EPSI Vitesse tangentielle de regularisation d' Oden * *
  16. * es XPALB Tableau contenant les param}tres de la liaison *
  17. * e NLIAB Nombre de liaisons sur la base B *
  18. * e I Num{ro de la liaison *
  19. * e I0 Pointeur de la vitesse tangentielle dans XPALB *
  20. * e I1 Pointeur de la force tangentielle dans XPALB *
  21. * e XFN Valeur de la force de choc suivant la normale. *
  22. * e XGLIS Valeur du coefficient de glissement *
  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. * I. Politopoulos ,Juin 1993 *
  30. * *
  31. *--------------------------------------------------------------------*
  32. DIMENSION XPALB(NLIAB,*)
  33. XFTT = - XGLIS * ABS(XFN)
  34. *
  35. PSVITT = 0.D0
  36. DO 10 ID = 1,IDIM
  37. PSVITT = PSVITT + XPALB(I,I0+ID) * XPALB(I,I0+ID)
  38. 10 CONTINUE
  39. * end do
  40. VITT = SQRT(PSVITT)
  41. PSXFT = 0.D0
  42. DO 11 ID = 1,IDIM
  43. XPALB(I,I1+ID) = XFTT *
  44. & (2- (VITT/EPSI)) * XPALB(I,I0+ID)/EPSI
  45. PSXFT = PSXFT + XPALB(I,I1+ID) * XPALB(I,I1+ID)
  46. 11 CONTINUE
  47. * end do
  48. XFT = SQRT(PSXFT)
  49. XPUS = ABS(XFN * VITT) * (VITT/EPSI)
  50.  
  51.  
  52. RETURN
  53. END
  54.  
  55.  

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