Télécharger dychfr.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHFR SOURCE BP208322 16/06/06 21:15:01 8944
  2. SUBROUTINE DYCHFR(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,I4,
  3. & XGLIS,XADHE,XRAIT,XAMOT,XFN,XFT,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. * Gestion des phases glissement et adh{rence du frottement de *
  12. * COULOMB *
  13. * *
  14. * Param}tres: *
  15. * *
  16. * e IGP Indice de glissement au pas pr{c{dent *
  17. * es XPALB Tableau contenant les param}tres de la liaison *
  18. * e NLIAB Nombre de liaisons sur la base B *
  19. * e I Num{ro de la liaison *
  20. * e I0 Pointeur de la vitesse tangentielle dans XPALB *
  21. * e I1 Pointeur de la variation de d{placement par rapport a *
  22. * la position origine d'adh{rence, projet{e sur le plan *
  23. * tangent, dans XPALB *
  24. * e I2 Pointeur de la force tangentielle dans XPALB *
  25. * e I3 Pointeur de la vitesse tangentielle au pas pr{c{dent *
  26. * dans XPALB *
  27. * e I4 Pointeur de la force tangentielle @ la fin du glis- *
  28. * -sement dans XPALB *
  29. * e XFN Valeur de la force de choc suivant la normale. *
  30. * e XGLIS Valeur du coefficient de glissement *
  31. * e XADHE Valeur du coefficient d'adh{rence *
  32. * e XRAIT Valeur de la raideur tangentielle *
  33. * e XAMOT Valeur de l'amortissement tangentiel *
  34. * s XFT Module de la force de choc projet{e sur le plan tangent*
  35. * s XPUS Valeur de la puissance d'usure instantan{e *
  36. * *
  37. * *
  38. * Auteur, date de cr{ation: *
  39. * *
  40. * Bertrand BEAUFILS, le 31 mai 1990. *
  41. * I. POLITOPOULOS (bricolage pour le model d'Oden) ,juin 1993 *
  42. *--------------------------------------------------------------------*
  43. DIMENSION XPALB(NLIAB,*)
  44. XFT = 0.D0
  45. XPUS = 0.D0
  46.  
  47. *---- modele d' Antunes ---------------------------------------------*
  48.  
  49. cbp,2016 IF (XRAIT.GT.0.D0) THEN
  50. IF (XRAIT.GE.0.D0) THEN
  51.  
  52. IF (IGP.EQ.1) THEN
  53. * Il y avait glissement au pas pr{cedent ; calcul du produit scalaire
  54. * de la vitesse et de la vitesse au pas pr{c{dent ; calcul du module
  55. * de la vitesse
  56. *
  57. PS = 0.D0
  58. DO 10 ID = 1,IDIM
  59. PS = PS + XPALB(I,I0+ID) * XPALB(I,I3+ID)
  60. 10 CONTINUE
  61. * end do
  62. IF (PS.GT.0.D0) THEN
  63. * Il y a glissement
  64. CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
  65. ELSE
  66. * Il y a possibilit{ d'adh{rence
  67. CALL DYCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I4,XADHE,XRAIT,
  68. & XAMOT,XFN,XFT,XPUS)
  69. IF (IGP.EQ.1) THEN
  70. * En fait, il y a glissement.
  71. CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,
  72. & XPUS)
  73. ELSE
  74. * Ca adh}re
  75. CONTINUE
  76. ENDIF
  77. ENDIF
  78. ELSE
  79. *
  80. * Il y avait adh{rence au pas pr{c{dent ; \a peut continuer
  81. *
  82. CALL DYCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I4,XADHE,XRAIT,
  83. & XAMOT,XFN,XFT,XPUS)
  84. IF (IGP.EQ.1) THEN
  85. * En fait, il y a glissement.
  86. CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
  87. ELSE
  88. * Ca adh}re
  89. CONTINUE
  90. ENDIF
  91. ENDIF
  92.  
  93.  
  94. *---- modele d'Oden -------------------------------------------------*
  95.  
  96. ELSE
  97. *
  98. IGP = 1
  99. EPSI = XAMOT
  100. * calcul du module de la vitesse
  101. *
  102. PSV = 0.D0
  103. DO 20 ID = 1,IDIM
  104. PSV = PSV + XPALB(I,I0+ID) * XPALB(I,I0+ID)
  105. 20 CONTINUE
  106.  
  107. * end do
  108. VITT = SQRT(PSV)
  109. IF (VITT.GT.EPSI) THEN
  110. * Il y a glissement classique
  111. CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
  112. ELSE
  113. * on est dans le domaine regularise
  114. CALL DYCGLR(IDIM,EPSI,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
  115. ENDIF
  116. ENDIF
  117.  
  118. RETURN
  119. END
  120.  
  121.  
  122.  

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