Télécharger dyche5.eso

Retour à la liste

Numérotation des lignes :

dyche5
  1. C DYCHE5 SOURCE CHAT 05/01/12 23:04:59 5004
  2. C DYCHE4 SOURCE LAVARENN 96/08/05 21:18:27 2237
  3. SUBROUTINE DYCHE5(XNET,XTE,XXXN,IDIM,IGP,XPALB,NLIAB,I
  4. &,IPALB,XFN,XFT,XPUS,iannul)
  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. * Appel du module de frottement pour la liaison *
  13. * ligne_cercle sans amortissement *
  14. * *
  15. * Param}tres: *
  16. * *
  17. * e XDEP Valeur du d{placement sur la normale *
  18. * e IGP Indice de glissement au pas pr{c{dent *
  19. * es XPALB Tableau contenant les param}tres de la liaison *
  20. * e NLIAB Nombre de liaison sur la base B *
  21. * e I Num{ro de la liaison *
  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. * IBRAHIM PINTO, 05/97 *
  30. * *
  31. *--------------------------------------------------------------------*
  32.  
  33.  
  34. REAL*8 XNET,XTE(3),XXXN(3)
  35.  
  36. DIMENSION XPALB(NLIAB,*)
  37. INTEGER IPALB(NLIAB,*)
  38. XFN = 0.D0
  39. XFT = 0.D0
  40. XPUS = 0.D0
  41. XVAL=0.D0
  42. XRAY = XPALB(I,2)
  43.  
  44. ITYP =IPALB(I,1)
  45.  
  46.  
  47. *
  48. IF ((XNET.GE.XRAY).AND.(iannul.EQ.0)) THEN
  49. ID2 = 6 + IDIM
  50. ID3 = 6 + 2*IDIM
  51. ID4 = 6 + 3*IDIM
  52. ID6 = 6 + 5*IDIM
  53. ID7 = 6 + 6*IDIM
  54. XRAIN = XPALB(I,1)
  55. XGLIS = XPALB(I,3)
  56. XADHE = XPALB(I,4)
  57. XRAIT = XPALB(I,5)
  58. XAMOT = XPALB(I,6)
  59.  
  60. *----dans le cas ou on suppose la normale de contact
  61. *----dans le plan du cercle on ne change pas xval
  62.  
  63. XVAL = XNET - XRAY
  64.  
  65. *----sinon il faut projeter sur la normale reelle de contact
  66. IF ( ITYP.EQ.39 ) THEN
  67. xval = 0.d0
  68. DO 4 ID=1,IDIM
  69. XVAL=XVAL+XTE(ID)*XXXN(ID)
  70. 4 CONTINUE
  71. XVAL=XVAL*(XNET-XRAY)/XNET
  72. ENDIF
  73. *----CALCUL DE LA FORCE DE CHOC NORMALE
  74.  
  75.  
  76. XFN = - XRAIN * XVAL
  77.  
  78. CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6,ID7,XGLIS
  79. &,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  80. ELSE
  81. ID4 = 6 + 3*IDIM
  82. IGP = -1
  83. DO 10 ID=1,IDIM
  84. XPALB(I,ID4+ID) = 0.D0
  85. 10 CONTINUE
  86. * end do
  87. ENDIF
  88. *
  89. RETURN
  90.  
  91.  
  92. END
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  

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