Télécharger dgchgl.eso

Retour à la liste

Numérotation des lignes :

  1. C DGCHGL SOURCE CHAT 05/01/12 22:48:36 5004
  2.  
  3. SUBROUTINE DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,
  4. + XGLIS,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 *
  46. *ON AVAIT GLISSEMENT AU PAS (n-1) AVEC IGP=1 *
  47. *HYPOTHESE: ON CONTINUE A GLISSER DANS LE MEME SENS *
  48. *
  49. *
  50. DIMENSION XPALB(NLIAB,*)
  51. DIMENSION XFORT (3)
  52. DIMENSION XVITG (3)
  53. XFT = - XGLIS * ABS(XFN)
  54. IGP = 1
  55. *
  56.  
  57. *
  58. *PSFT = PRODUIT SCALAIRE DE LA FORCE PAR ELLE MEME *
  59. PSFT = 0.D0
  60. DO 30 ID = 1,IDIM
  61. PSFT = PSFT + XPALB(I,I2+ID) * XPALB(I,I2+ID)
  62. 30 CONTINUE
  63. *END DO
  64. *
  65. PSFT = SQRT(PSFT)
  66. *
  67. *CALCUL DE LA FORCE DE FROTTEMENT XFORT(n)
  68. DO 40 ID = 1,IDIM
  69. XFORT(ID) = XGLIS * ABS(XFN) * XPALB(I,I2+ID) / PSFT
  70.  
  71. 40 CONTINUE
  72. *END DO
  73. *
  74. *CALCUL DE LA VITESSE DU GLISSEUR VG(n)
  75. DO 50 ID = 1,IDIM
  76. XVITG(ID) = XPALB(I,I0+ID) + ((XFORT(ID)
  77. + - XPALB(I,I2+ID)) / (XRAIT * XPAS + XAMOT))
  78. + - (XAMOT * (XPALB(I,I1+ID)
  79. + - XPALB(I,I3+ID)) / (XRAIT * XPAS + XAMOT))
  80. *
  81.  
  82. 50 CONTINUE
  83. *END DO
  84.  
  85.  
  86. RETURN
  87. END
  88.  
  89.  
  90.  
  91.  

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