Télécharger dycha31.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHA31 SOURCE BP208322 20/03/26 21:15:43 10562
  2.  
  3. C DYCHA3 SOURCE CHAT 05/01/12 23:00:43 5004
  4. SUBROUTINE DYCHA31(XDEP,XVITN,IDIM,IGP,XPALB,NLIAB,I,INTER
  5. & ,XFN,XFT,XPUS,iannul)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. *--------------------------------------------------------------------*
  9. * *
  10. * Operateur DYNE : *
  11. * Calcul de la force de choc pour une liaison B de type : *
  12. * - POINT_CERCLE_FROTTEMENT *
  13. * *
  14. *--------------------------------------------------------------------*
  15. * Parametres: *
  16. * *
  17. * e XDEP Valeur du deplacement sur la normale *
  18. * e XVITN Valeur de la vitesse normale au plan tangent du contact*
  19. * e IGP Indice de glissement au pas precedent *
  20. * es XPALB Tableau contenant les parametres de la liaison *
  21. * e NLIAB Nombre de liaisons sur la base B *
  22. * e I Numero de la liaison *
  23. * e INTER Indice de liaison interieure / exterieure *
  24. * s XFN Valeur de la force de choc suivant la normale. *
  25. * s XFT Module de la force de choc projetee sur le plan tangent*
  26. * s XPUS Valeur de la puissance d'usure instantanee *
  27. * *
  28. *--------------------------------------------------------------------*
  29. DIMENSION XPALB(NLIAB,*)
  30. XFN = 0.D0
  31. XFT = 0.D0
  32. XPUS = 0.D0
  33. XJEU = XPALB(I,2)
  34. *
  35. XAD = ABS(XDEP)
  36. XAJ = ABS(XJEU)
  37. *
  38. IF ( (((XAD.LE.XAJ).AND.(INTER.EQ.0)).OR.
  39. & ((XAD.GE.XAJ).AND.(INTER.EQ.1))).AND.
  40. & (iannul.EQ.0)) THEN
  41. c ID4 = 7 + 3*IDIM
  42. c ID5 = 7 + 4*IDIM
  43. c ID6 = 7 + 5*IDIM
  44. c ID8 = 7 + 7*IDIM
  45. c ID9 = 7 + 8*IDIM
  46. ID4 = 10 + 3*IDIM
  47. ID5 = 10 + 4*IDIM
  48. ID6 = 10 + 5*IDIM
  49. ID8 = 10 + 7*IDIM
  50. ID9 = 10 + 8*IDIM
  51. XRAIN = XPALB(I,1)
  52. XGLIS = XPALB(I,3)
  53. XADHE = XPALB(I,4)
  54. XRAIT = XPALB(I,5)
  55. XAMOT = XPALB(I,6)
  56. XAMON = XPALB(I,7)
  57. XVAL = XDEP - XJEU
  58. XFN = - XRAIN * XVAL - XAMON * XVITN
  59. CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID4,ID5,ID6,ID8,ID9,XGLIS,
  60. & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  61. ELSE
  62. IGP = -1
  63. c ID6 = 7 + 5*IDIM
  64. ID6 = 10 + 5*IDIM
  65. DO 10 ID=1,IDIM
  66. XPALB(I,ID6+ID) = 0.D0
  67. 10 CONTINUE
  68. ENDIF
  69. *
  70. RETURN
  71. END
  72.  
  73.  
  74.  
  75.  
  76.  

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