Télécharger la plaquette Cast3M


fr en

Télécharger dychgl.eso

Retour à la liste

Numérotation des lignes : 

dychgl

  1. C DYCHGL    SOURCE    CHAT      05/01/12    23:06:07     5004
  2.       SUBROUTINE DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,XGLIS,XFN,XFT,
  3.      &                  XPUS)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. *                                                                    *
  8. *     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *
  9. *     ________________________________________________               *
  10. *                                                                    *
  11. *     Calcul de la force de glissement                               *
  12. *                                                                    *
  13. *     Param}tres:                                                    *
  14. *                                                                    *
  15. * es  XPALB   Tableau contenant les param}tres de la liaison         *
  16. * e   NLIAB   Nombre de liaisons sur la base B                       *
  17. * e   I       Num{ro de la liaison                                   *
  18. * e   I0      Pointeur de la vitesse tangentielle dans XPALB         *
  19. * e   I1      Pointeur de la force tangentielle dans XPALB           *
  20. * e   XFN     Valeur de la force de choc suivant la normale.         *
  21. * e   XGLIS   Valeur du coefficient d'adh{rence                      *
  22. *  s  IGP     Indice de glissement au pas pr{c{dent                  *
  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. *     Bertrand BEAUFILS, le 31 mai 1990.                             *
  30. *                                                                    *
  31. *--------------------------------------------------------------------*
  32.       DIMENSION XPALB(NLIAB,*)
  33.       XFT = - XGLIS * ABS(XFN)
  34.       IGP = 1
  35. *
  36.       PSVITT = 0.D0
  37.       DO 10 ID = 1,IDIM
  38.           PSVITT = PSVITT + XPALB(I,I0+ID) * XPALB(I,I0+ID)
  39.  10       CONTINUE
  40. *     end do
  41.       VITT = SQRT(PSVITT)
  42.       XPUS = ABS(XFN * VITT)
  43. *
  44.       IF (VITT.GT.1.D-20) THEN
  45.          DO 11 ID = 1,IDIM
  46.             XPALB(I,I1+ID) = XFT * XPALB(I,I0+ID) / VITT
  47.  11         CONTINUE
  48. *        end do
  49.       ELSE
  50.          DO 12 ID = 1,IDIM
  51.             XPALB(I,I1+ID) = 0.D0
  52.  12         CONTINUE
  53. *        end do
  54.       ENDIF
  55. *
  56.       RETURN
  57.       END
  58.  
  59.  

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