Télécharger dgchfr.eso

Retour à la liste

Numérotation des lignes :

  1. C DGCHFR SOURCE CHAT 05/01/12 22:48:33 5004
  2. SUBROUTINE DGCHFR(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,I4,
  3. + XGLIS,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS,XPAS)
  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. * LE PAS DE TEMPS DOIT ETRE CHOISI DE TELLE FACON QUE LE *
  12. * RESSORT TANGENTIEL NE PUISSE PAS PASSER DE LA COMPRESSION *
  13. * ET A L EXTENSION DURANT UN PAS DE TEMPS CE QUI EXCLUT LE *
  14. * PASSAGE BRUTAL D UN GLISSEMENT DANS UN SENS A UN GLISSEMENT *
  15. * DANS L AUTRE SENS *
  16. * *
  17. * *
  18. * Gestion des phases glissement et adh{rence du frottement de *
  19. * COULOMB *
  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. DIMENSION XPALB(NLIAB,*)
  46. DIMENSION XFORT (3)
  47. DIMENSION XVITG (3)
  48. XFT = 0.D0
  49. XPUS = 0.D0
  50.  
  51. * modele d' Antunez
  52. IF (IGP.EQ.1) THEN
  53. *
  54. *CAS1 IL Y AVAIT GLISSEMENT AU PAS PRECEDENT ON SUPPOSE QUE LE *
  55. *GLISSEMENT CONTINUE DANS LE MEME SENS *
  56. CALL DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XGLIS,
  57. + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
  58. *ON FAIT LE TEST SUR LES VITESSES DU GLISSEUR *
  59. *PRODUIT SCALAIRE DE LA VITESSE DU GLISSEUR AU PAS n ET AU *
  60. *PAS n-1 *
  61. *
  62. PSVITG = 0.D0
  63. DO 10 ID = 1,IDIM
  64. PSVITG = PSVITG + XVITG(ID) * XPALB(I,I3+ID)
  65. *
  66. 10 CONTINUE
  67. * end do
  68. IF (PSVITG.GT.0.D0) THEN
  69. *IL Y A GLISSEMENT
  70. IGP = 1
  71. *
  72. DO 11 ID = 1,IDIM
  73. XPALB(I,I2+ID) = XFORT(ID)
  74. XPALB(I,I3+ID) = XVITG(ID)
  75.  
  76. 11 CONTINUE
  77. * end do
  78. ELSE
  79. CALL DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XADHE,
  80. + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
  81. IGP = 0
  82. DO 30 ID = 1,IDIM
  83. XPALB(I,I2+ID) = XFORT(ID)
  84. XPALB(I,I3+ID) = XVITG(ID)
  85. 30 CONTINUE
  86.  
  87. ENDIF
  88. ELSE
  89. *(ELSE DE IGP=1 DU DEBUT)
  90. *CAS2 IL Y AVAIT ADHERENCE AU PAS PRECEDENT ,ON SUPPOSE QUE *
  91. *L ADHERENCE CONTINUE
  92.  
  93. CALL DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XADHE,
  94. + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
  95. *
  96.  
  97. *ON FAIT LE TEST
  98. IF (XFT.GT.(XADHE * ABS(XFN))) THEN
  99. *L ADHERENCE NE TIENT PAS
  100. IGP = 1
  101. CALL DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XGLIS,
  102. + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
  103.  
  104. DO 40 ID = 1,IDIM
  105. XPALB(I,I2+ID) = XFORT(ID)
  106. XPALB(I,I3+ID) = XVITG(ID)
  107. ***ICI ON RAJOUTE BOUCLE 40
  108.  
  109. 40 CONTINUE
  110. *End do
  111. ELSE
  112. IGP = 0
  113.  
  114. DO 50 ID = 1,IDIM
  115. XPALB(I,I2+ID) = XFORT(ID)
  116. XPALB(I,I3+ID) = XVITG(ID)
  117. ***ICI ON RAJOUTE BOUCLE 50
  118.  
  119. 50 CONTINUE
  120.  
  121. endif
  122. endif
  123. end
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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