Télécharger la plaquette Cast3M


fr en

Télécharger dgchad.eso

Retour à la liste

Numérotation des lignes : 

dgchad

  1. C DGCHAD    SOURCE    CHAT      05/01/12    22:48:29     5004
  2.  
  3.         SUBROUTINE DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,
  4.      +  XADHE,XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
  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. *     LE PAS DE TEMPS DOIT ETRE CHOISI DE TELLE FACON QUE LE         *
  13. *     RESSORT TANGENTIEL NE PUISSE PAS PASSER DE LA COMPRESSION      *
  14. *     ET A L EXTENSION DURANT UN PAS DE TEMPS CE QUI EXCLUT LE       *
  15. *     PASSAGE BRUTAL D UN GLISSEMENT DANS UN SENS A UN GLISSEMENT    *
  16. *     DANS L AUTRE SENS
  17. *                                                                    *
  18. *                                                                    *
  19. *                                                                    *
  20. *                                                                    *
  21. *     Param}tres:                                                    *
  22. *                                                                    *
  23. * e   IGP     Indice de glissement au pas pr{c{dent                  *
  24. * es  XPALB   Tableau contenant les param}tres de la liaison         *
  25. * e   NLIAB   Nombre de liaisons sur la base B
  26. * e   I       Num{ro de la liaison                                   *
  27. * e   I0      Pointeur de la vitesse tangentielle dans XPALB         *
  28. * e   I1      Pointeur de la vitesse tangentielle au pas pr{c{dent   *
  29. *             dans XPALB                                             *
  30. * e   I2      Pointeur de la force tangentielle au pas               *
  31. *             pr{c{dent dans XPALB                                   *
  32. * e   I3      Pointeur de la vitesse du glisseur au pas
  33. *             pr{c{dent dans XPALB                                   *
  34. * e   XFN     Valeur de la force de choc suivant la normale.         *
  35. * e   XGLIS   Valeur du coefficient de glissement                    *
  36. * e   XADHE   Valeur du coefficient d'adh{rence                      *
  37. * e   XRAIT   Valeur de la raideur tangentielle                      *
  38. * e   XAMOT   Valeur de l'amortissement tangentiel                   *
  39. *  s  XFT     Module de la force de choc projet{e sur le plan tangent*
  40. *  s  XPUS    Valeur de la puissance d'usure instantan{e             *
  41. *  e  XPAS    Pas de temps                                           *                                                                *
  42. *                                                                    *
  43. * Auteurs Heidi NEDJAI  Rene GIBERT                                               *
  44. *--------------------------------------------------------------------*
  45. *CALCUL DES FORCES FT(n)=XFORT ET VG(n)=XVITG = 0                    *
  46. *ON AVAIT ADHERENCE AU PAS (n-1) AVEC IGP=0                          *
  47. *HYPOTHESE: ON CONTINUE A ADHERER AU PAS (n)                         *
  48. *
  49. *
  50.         DIMENSION XPALB(NLIAB,*)
  51.         DIMENSION XFORT (3)
  52.         DIMENSION XVITG (3)
  53.         DO 10 ID = 1,IDIM
  54.         XVITG(ID) = 0.D0
  55.   10    CONTINUE
  56. *END DO
  57. * RECHERCHE DE XRAIT XPAS XAMOT
  58.  
  59.         PSXFT = 0.D0
  60.         DO 20 ID = 1,IDIM
  61.         XFORT(ID) = XPALB(I,I2+ID)
  62.      +        - (XRAIT * XPAS + XAMOT) * XPALB(I,I0+ID)
  63.      +        - (XAMOT * (XPALB(I,I1+ID) - XPALB(I,I3+ID)))
  64. *
  65.          PSXFT = PSXFT + XFORT(ID) * XFORT(ID)
  66.  
  67.   20     CONTINUE
  68.         XFT = SQRT(PSXFT)
  69.  
  70.         RETURN
  71.         END
  72.  
  73.  

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