Télécharger dgchad.eso

Retour à la liste

Numérotation des lignes :

  1. C DGCHAD SOURCE CHAT 05/01/12 22:48:29 5004
  2.  
  3. SUBROUTINE DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,
  4. + XADHE,XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
  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. * LE PAS DE TEMPS DOIT ETRE CHOISI DE TELLE FACON QUE LE *
  13. * RESSORT TANGENTIEL NE PUISSE PAS PASSER DE LA COMPRESSION *
  14. * ET A L EXTENSION DURANT UN PAS DE TEMPS CE QUI EXCLUT LE *
  15. * PASSAGE BRUTAL D UN GLISSEMENT DANS UN SENS A UN GLISSEMENT *
  16. * DANS L AUTRE SENS
  17. * *
  18. * *
  19. * *
  20. * *
  21. * Param}tres: *
  22. * *
  23. * e IGP Indice de glissement au pas pr{c{dent *
  24. * es XPALB Tableau contenant les param}tres de la liaison *
  25. * e NLIAB Nombre de liaisons sur la base B
  26. * e I Num{ro de la liaison *
  27. * e I0 Pointeur de la vitesse tangentielle dans XPALB *
  28. * e I1 Pointeur de la vitesse tangentielle au pas pr{c{dent *
  29. * dans XPALB *
  30. * e I2 Pointeur de la force tangentielle au pas *
  31. * pr{c{dent dans XPALB *
  32. * e I3 Pointeur de la vitesse du glisseur au pas
  33. * pr{c{dent dans XPALB *
  34. * e XFN Valeur de la force de choc suivant la normale. *
  35. * e XGLIS Valeur du coefficient de glissement *
  36. * e XADHE Valeur du coefficient d'adh{rence *
  37. * e XRAIT Valeur de la raideur tangentielle *
  38. * e XAMOT Valeur de l'amortissement tangentiel *
  39. * s XFT Module de la force de choc projet{e sur le plan tangent*
  40. * s XPUS Valeur de la puissance d'usure instantan{e *
  41. * e XPAS Pas de temps * *
  42. * *
  43. * Auteurs Heidi NEDJAI Rene GIBERT *
  44. *--------------------------------------------------------------------*
  45. *CALCUL DES FORCES FT(n)=XFORT ET VG(n)=XVITG = 0 *
  46. *ON AVAIT ADHERENCE AU PAS (n-1) AVEC IGP=0 *
  47. *HYPOTHESE: ON CONTINUE A ADHERER AU PAS (n) *
  48. *
  49. *
  50. DIMENSION XPALB(NLIAB,*)
  51. DIMENSION XFORT (3)
  52. DIMENSION XVITG (3)
  53. DO 10 ID = 1,IDIM
  54. XVITG(ID) = 0.D0
  55. 10 CONTINUE
  56. *END DO
  57. * RECHERCHE DE XRAIT XPAS XAMOT
  58.  
  59. PSXFT = 0.D0
  60. DO 20 ID = 1,IDIM
  61. XFORT(ID) = XPALB(I,I2+ID)
  62. + - (XRAIT * XPAS + XAMOT) * XPALB(I,I0+ID)
  63. + - (XAMOT * (XPALB(I,I1+ID) - XPALB(I,I3+ID)))
  64. *
  65. PSXFT = PSXFT + XFORT(ID) * XFORT(ID)
  66.  
  67. 20 CONTINUE
  68. XFT = SQRT(PSXFT)
  69.  
  70. RETURN
  71. END
  72.  
  73.  

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