Télécharger la plaquette Cast3M


fr en

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