Télécharger la plaquette Cast3M


fr en

Télécharger dycglr.eso

Retour à la liste

Numérotation des lignes : 

dycglr

  1. C DYCGLR    SOURCE    CHAT      05/01/12    23:00:14     5004
  2. C DYCGLR    SOURCE    BEAU      90/07/26    21:16:02
  3.       SUBROUTINE DYCGLR(IDIM,EPSI,XPALB,NLIAB,I,I0,I1,XGLIS,XFN,XFT,
  4.      &                  XPUS)
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. *                                                                    *
  9. *     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *
  10. *     ________________________________________________               *
  11. *                                                                    *
  12. *     Calcul de la force de glissement                               *
  13. *                                                                    *
  14. *     Param}tres:                                                    *
  15. * e   EPSI    Vitesse tangentielle de regularisation d' Oden         *                                                                   *
  16. * es  XPALB   Tableau contenant les param}tres de la liaison         *
  17. * e   NLIAB   Nombre de liaisons sur la base B                       *
  18. * e   I       Num{ro de la liaison                                   *
  19. * e   I0      Pointeur de la vitesse tangentielle dans XPALB         *
  20. * e   I1      Pointeur de la force tangentielle dans XPALB           *
  21. * e   XFN     Valeur de la force de choc suivant la normale.         *
  22. * e   XGLIS   Valeur du coefficient de glissement                    *
  23. *  s  XFT     Module de la force de choc projet{e sur le plan tangent*
  24. *  s  XPUS    Valeur de la puissance d'usure instantan{e             *
  25. *                                                                    *
  26. *                                                                    *
  27. *     Auteur, date de cr{ation:                                      *
  28. *                                                                    *
  29. *    I. Politopoulos ,Juin 1993                                      *
  30. *                                                                    *
  31. *--------------------------------------------------------------------*
  32.       DIMENSION XPALB(NLIAB,*)
  33.       XFTT = - XGLIS * ABS(XFN)
  34. *
  35.       PSVITT = 0.D0
  36.       DO 10 ID = 1,IDIM
  37.           PSVITT = PSVITT + XPALB(I,I0+ID) * XPALB(I,I0+ID)
  38.  10       CONTINUE
  39. *     end do
  40.       VITT = SQRT(PSVITT)
  41.       PSXFT = 0.D0
  42.       DO 11 ID = 1,IDIM
  43.           XPALB(I,I1+ID) = XFTT *
  44.      &                   (2- (VITT/EPSI)) * XPALB(I,I0+ID)/EPSI
  45.           PSXFT = PSXFT + XPALB(I,I1+ID) * XPALB(I,I1+ID)
  46.  11       CONTINUE
  47. *     end do
  48.       XFT = SQRT(PSXFT)
  49.       XPUS = ABS(XFN * VITT) * (VITT/EPSI)
  50.  
  51.  
  52.       RETURN
  53.       END
  54.  
  55.  

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