Télécharger la plaquette Cast3M


fr en

Télécharger dycha31.eso

Retour à la liste

Numérotation des lignes : 

dycha31

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

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