Télécharger dycha11.eso

Retour à la liste

Numérotation des lignes :

dycha11
  1. C DYCHA11 SOURCE BP208322 20/03/26 21:15:42 10562
  2. c
  3. C DYCHA41 SOURCE CHAT 05/01/12 23:01:19 5004
  4. SUBROUTINE DYCHA11(XDEP,XVITN,IDIM,IGP,XPALB,NLIAB,I,XFN,XFT,XPUS,
  5. & XABSCI,XORDON,NIP,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 de type : *
  12. * - POINT_PLAN_FROTTEMENT avec LOI_DE_COMPORTEMENT *
  13. * *
  14. *--------------------------------------------------------------------*
  15. * *
  16. * Parametres: *
  17. * *
  18. * e XDEP Valeur du deplacement sur la normale *
  19. * e XVITN Valeur de la vitesse normale au plan tangent du contact*
  20. * e IGP Indice de glissement au pas precedent *
  21. * e/s XPALB Tableau contenant les parametres de la liaison *
  22. * e NLIAB Nombre de liaisons sur la base B *
  23. * e I Numero de la liaison *
  24. * e XABSCI Abscisses loi elastique *
  25. * e XORDON Ordonnees loi elastique *
  26. * s XFN Valeur de la force de choc suivant la normale. *
  27. * s XFT Module de la force de choc projetee sur le plan tangent*
  28. * s XPUS Valeur de la puissance d'usure instantanee *
  29. * *
  30. *--------------------------------------------------------------------*
  31. DIMENSION XPALB(NLIAB,*)
  32. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  33. XFN = 0.D0
  34. XFT = 0.D0
  35. XPUS = 0.D0
  36. XJEU = XPALB(I,2)
  37. *
  38. IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
  39. * on recupere les grandeurs scalaires
  40. XAMON = XPALB(I,7)
  41. XGLIS = XPALB(I,3)
  42. XADHE = XPALB(I,4)
  43. XRAIT = XPALB(I,5)
  44. XAMOT = XPALB(I,6)
  45. c ID2 = 7 + IDIM
  46. c ID3 = 7 + 2*IDIM
  47. c ID4 = 7 + 3*IDIM
  48. c ID6 = 7 + 5*IDIM
  49. c ID7 = 7 + 6*IDIM
  50. ID2 = 9 + 2*IDIM
  51. ID3 = 9 + 3*IDIM
  52. ID4 = 9 + 4*IDIM
  53. ID6 = 9 + 6*IDIM
  54. ID7 = 9 + 7*IDIM
  55. c XVAL = XDEP - XJEU
  56. CALL LIRANG(XDEP,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  57. XFN = - (XORDON(I,NRG) + XPENTE*(XDEP-XABSCI(I,NRG)))
  58. & - XAMON * XVITN
  59. CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6,ID7,XGLIS,
  60. & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  61. IF (XFN.GT.0.D0) XFN = 0.D0
  62. ELSE
  63. IGP = -1
  64. c ID4 = 7 + 3*IDIM
  65. ID4 = 9 + 4*IDIM
  66. DO 10 ID=1,IDIM
  67. XPALB(I,ID4+ID) = 0.D0
  68. 10 CONTINUE
  69. ENDIF
  70. *
  71. RETURN
  72. END
  73.  
  74.  
  75.  
  76.  
  77.  

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