Télécharger dycha3.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHA3 SOURCE CHAT 05/01/12 23:00:43 5004
  2. SUBROUTINE DYCHA3(XDEP,XVITN,IDIM,IGP,XPALB,NLIAB,I,INTER
  3. & ,XFN,XFT,XPUS,iannul)
  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. * Calcul de la force de choc pour un choc {l{mentaire. *
  12. * *
  13. * Param}tres: *
  14. * *
  15. * e XDEP Valeur du d{placement sur la normale *
  16. * e XVITN Valeur de la vitesse normale au plan tangent du contact*
  17. * e IGP Indice de glissement au pas pr{c{dent *
  18. * es XPALB Tableau contenant les param}tres de la liaison *
  19. * e NLIAB Nombre de liaisons sur la base B *
  20. * e I Num{ro de la liaison *
  21. * e INTER Indice de liaison interieure / exterieure *
  22. * s XFN Valeur de la force de choc suivant la normale. *
  23. * s XFT Module de la force de choc projetee sur le plan tangent*
  24. * s XPUS Valeur de la puissance d'usure instantanee *
  25. * *
  26. * *
  27. * Auteur, date de cr{ation: *
  28. * *
  29. * Bertrand BEAUFILS, le 31 mai 1990. *
  30. * *
  31. *--------------------------------------------------------------------*
  32. DIMENSION XPALB(NLIAB,*)
  33. XFN = 0.D0
  34. XFT = 0.D0
  35. XPUS = 0.D0
  36. XJEU = XPALB(I,2)
  37. *
  38. XAD = ABS(XDEP)
  39. XAJ = ABS(XJEU)
  40. *
  41. IF ( (((XAD.LE.XAJ).AND.(INTER.EQ.0)).OR.
  42. & ((XAD.GE.XAJ).AND.(INTER.EQ.1))).AND.
  43. & (iannul.EQ.0)) THEN
  44. ID4 = 7 + 3*IDIM
  45. ID5 = 7 + 4*IDIM
  46. ID6 = 7 + 5*IDIM
  47. ID8 = 7 + 7*IDIM
  48. ID9 = 7 + 8*IDIM
  49. XRAIN = XPALB(I,1)
  50. XAMON = XPALB(I,7)
  51. XGLIS = XPALB(I,3)
  52. XADHE = XPALB(I,4)
  53. XRAIT = XPALB(I,5)
  54. XAMOT = XPALB(I,6)
  55. XVAL = XDEP - XJEU
  56. XFN = - XRAIN * XVAL - XAMON * XVITN
  57. CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID4,ID5,ID6,ID8,ID9,XGLIS,
  58. & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  59. ELSE
  60. IGP = -1
  61. ID6 = 7 + 5*IDIM
  62. DO 10 ID=1,IDIM
  63. XPALB(I,ID6+ID) = 0.D0
  64. 10 CONTINUE
  65. * end do
  66. ENDIF
  67. *
  68. RETURN
  69. END
  70.  
  71.  
  72.  
  73.  

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