Télécharger la plaquette Cast3M


fr en

Télécharger dycha41.eso

Retour à la liste

Numérotation des lignes : 

dycha41

  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