Télécharger dycha5.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHA5 SOURCE CHAT 05/01/12 23:02:11 5004
  2. C DYCHA4 SOURCE LAVARENN 96/08/05 21:17:45 2237
  3. SUBROUTINE DYCHA5(XNET,XTE,XXXN,XVITN,IDIM,IGP
  4. &,XPALB,NLIAB,I,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. * Calcul de la force de choc pour un choc elementaire *
  13. * avec amortissement pour la liaison ligne_cercle *
  14. * Param}tres: *
  15. * *
  16. * e XDEP Valeur du d{placement sur la normale *
  17. * e XVITN Valeur de la vitesse normale au plan tangent du contact*
  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 liaisons 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.  
  35. REAL*8 XNET,XTE(3),XXXN(3)
  36.  
  37.  
  38.  
  39. DIMENSION XPALB(NLIAB,*)
  40. INTEGER IPALB(NLIAB,*)
  41. XFN = 0.D0
  42. XFT = 0.D0
  43. XPUS = 0.D0
  44. XRAY = XPALB(I,2)
  45. XRAIT = XPALB(I,5)
  46.  
  47.  
  48. ITYP =IPALB(I,1)
  49.  
  50. *
  51. IF ((XNET.GE.XRAY).AND.(iannul.EQ.0)) THEN
  52. ID2 = 7 + IDIM
  53. ID3 = 7 + 2*IDIM
  54. ID4 = 7 + 3*IDIM
  55. ID6 = 7 + 5*IDIM
  56. ID7 = 7 + 6*IDIM
  57. XRAIN = XPALB(I,1)
  58. XAMON = XPALB(I,7)
  59. XGLIS = XPALB(I,3)
  60. XADHE = XPALB(I,4)
  61. XRAIT = XPALB(I,5)
  62. XAMOT = XPALB(I,6)
  63.  
  64.  
  65. *----dans le cas ou on suppose la normale de contact
  66. *----dans le plan du cercle on ne change pas xval
  67.  
  68. XVAL = XNET - XRAY
  69.  
  70.  
  71. *----sinon il faut projeter sur la normale reelle de contact
  72. IF (ITYP.EQ.40) THEN
  73.  
  74.  
  75. XVAL=0.D0
  76. DO 4 ID=1,IDIM
  77. XVAL=XVAL+XTE(ID)*XXXN(ID)
  78. 4 CONTINUE
  79. XVAL=XVAL*(XNET-XRAY)/XNET
  80. ENDIF
  81.  
  82.  
  83. *----CALCUL DE LA FORCE DE CHOC NORMALE
  84.  
  85. XFN = - XRAIN * XVAL - XAMON * XVITN
  86.  
  87.  
  88. CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6
  89. &,ID7,XGLIS,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  90.  
  91. IF ( XFN.GT.0.D0) XFN = 0.D0
  92. ELSE
  93. IGP = -1
  94. ID4 = 7 + 3*IDIM
  95. DO 10 ID=1,IDIM
  96. XPALB(I,ID4+ID) = 0.D0
  97. 10 CONTINUE
  98. * end do
  99. ENDIF
  100. *
  101.  
  102. RETURN
  103.  
  104.  
  105. END
  106.  
  107.  
  108.  
  109.  
  110.  

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