Télécharger dycha41.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHA41 SOURCE CHAT 05/01/12 23:01:19 5004
  2. C DYCHA41 SOURCE LAVARENN 96/08/05 21:17:45 2237
  3. SUBROUTINE DYCHA41(XDEP,XVITN,IDIM,IGP,XPALB,NLIAB,I,XFN,XFT
  4. & ,XPUS,XABSCI,XORDON,NIP,iannul)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. * *
  9. * Op{rateur DYNE : *
  10. * ________________________________________________ *
  11. * *
  12. * Calcul de la force de choc pour un choc {l{mentaire. *
  13. * Cas d'une loi elastique *
  14. * *
  15. * Param}tres: *
  16. * *
  17. * e XDEP Valeur du d{placement sur la normale *
  18. * e XVITN Valeur de la vitesse normale au plan tangent du contact*
  19. * e IGP Indice de glissement au pas pr{c{dent *
  20. * es XPALB Tableau contenant les param}tres de la liaison *
  21. * e NLIAB Nombre de liaisons sur la base B *
  22. * e I Num{ro de la liaison *
  23. * e XABSCI Abscisses loi elastique *
  24. * e XORDON Ordonnees loi elastique *
  25. * s XFN Valeur de la force de choc suivant la normale. *
  26. * s XFT Module de la force de choc projetee sur le plan tangent*
  27. * s XPUS Valeur de la puissance d'usure instantanee *
  28. * *
  29. * *
  30. * Auteur, date de cr{ation: *
  31. * *
  32. * Bruno FONTAINE, le 26/02/98 *
  33. * *
  34. *--------------------------------------------------------------------*
  35. DIMENSION XPALB(NLIAB,*)
  36. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  37. XFN = 0.D0
  38. XFT = 0.D0
  39. XPUS = 0.D0
  40. XJEU = XPALB(I,2)
  41. *
  42. IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
  43. ID2 = 7 + IDIM
  44. ID3 = 7 + 2*IDIM
  45. ID4 = 7 + 3*IDIM
  46. ID6 = 7 + 5*IDIM
  47. ID7 = 7 + 6*IDIM
  48. XAMON = XPALB(I,7)
  49. XGLIS = XPALB(I,3)
  50. XADHE = XPALB(I,4)
  51. XRAIT = XPALB(I,5)
  52. XAMOT = XPALB(I,6)
  53. XVAL = XDEP - XJEU
  54. CALL LIRANG(XDEP,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  55. XFN = - (XORDON(I,NRG) + XPENTE*(XDEP-XABSCI(I,NRG)))
  56. & - XAMON * XVITN
  57. CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6,ID7,XGLIS,
  58. & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  59. IF ( XFN.GT.0.D0) XFN = 0.D0
  60. ELSE
  61. IGP = -1
  62. ID4 = 7 + 3*IDIM
  63. DO 10 ID=1,IDIM
  64. XPALB(I,ID4+ID) = 0.D0
  65. 10 CONTINUE
  66. * end do
  67. ENDIF
  68. *
  69. RETURN
  70. END
  71.  
  72.  
  73.  
  74.  

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