Télécharger la plaquette Cast3M


fr en

Télécharger dgchgl.eso

Retour à la liste

Numérotation des lignes : 

dgchgl

  1. C DGCHGL    SOURCE    CHAT      05/01/12    22:48:36     5004
  2.  
  3.         SUBROUTINE DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,
  4.      +  XGLIS,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                        *
  46. *ON AVAIT GLISSEMENT AU PAS (n-1) AVEC IGP=1                         *
  47. *HYPOTHESE: ON CONTINUE A GLISSER DANS LE MEME SENS                  *
  48. *
  49. *
  50.         DIMENSION XPALB(NLIAB,*)
  51.         DIMENSION XFORT (3)
  52.         DIMENSION XVITG (3)
  53.         XFT = - XGLIS * ABS(XFN)
  54.         IGP = 1
  55. *
  56.  
  57. *
  58. *PSFT = PRODUIT SCALAIRE DE LA FORCE PAR ELLE MEME                   *
  59.         PSFT = 0.D0
  60.         DO 30 ID = 1,IDIM
  61.         PSFT = PSFT + XPALB(I,I2+ID) * XPALB(I,I2+ID)
  62.   30    CONTINUE
  63. *END DO
  64. *
  65.         PSFT = SQRT(PSFT)
  66. *
  67. *CALCUL DE LA FORCE DE FROTTEMENT XFORT(n)
  68.         DO 40 ID = 1,IDIM
  69.         XFORT(ID) = XGLIS * ABS(XFN) * XPALB(I,I2+ID) / PSFT
  70.  
  71.   40    CONTINUE
  72. *END DO
  73. *
  74. *CALCUL DE LA VITESSE DU GLISSEUR VG(n)
  75.         DO 50 ID = 1,IDIM
  76.         XVITG(ID) = XPALB(I,I0+ID) + ((XFORT(ID)
  77.      +     - XPALB(I,I2+ID)) / (XRAIT * XPAS + XAMOT))
  78.      +     - (XAMOT * (XPALB(I,I1+ID)
  79.      +     - XPALB(I,I3+ID)) / (XRAIT * XPAS + XAMOT))
  80. *
  81.  
  82.   50    CONTINUE
  83. *END DO
  84.  
  85.  
  86.        RETURN
  87.        END
  88.  
  89.  
  90.  
  91.  

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