Télécharger la plaquette Cast3M


fr en

Télécharger palier.eso

Retour à la liste

Numérotation des lignes : 

palier

  1. C PALIER    SOURCE    BP208322  18/07/11    21:15:18     9879           
  2.  
  3. **********************************************************************
  4. *                                                                    *
  5. * Sous-programme de calcul des efforts engendrés par un palier       *
  6. * avec le modèle de palier court                                     *
  7. * Benoit PRABEL, 29/04/2015                                          *
  8. *                                                                    *
  9. **********************************************************************
  10.  
  11.       SUBROUTINE PALIER(X,Y,VX,VY,XAD,YAD,VXAD,VYAD,FX,FY,PDT,ICAVIT,
  12.      &     XLONG,VISCDY,RARBRE,XJEU,VITROT,NBMAIL,PALM,XPROG)
  13. *
  14. **********************************************************************
  15. *
  16.       IMPLICIT INTEGER(I-N)
  17.       IMPLICIT REAL*8 (A-H,O-Z)
  18. *
  19.  
  20. -INC PPARAM
  21. -INC CCOPTIO
  22. -INC CCREEL
  23. *
  24.       REAL*8 X,Y,VX,VY,XAD,YAD,VXAD,VYAD,FX,FY,PDT,AMPLIT,XLONG,VISCDY,
  25.      &     RARBRE,XJEU,VITROT,PALM
  26.       REAL*8 EMAX,EMIN
  27.       REAL*8 XPROG(2,NBMAIL)
  28.  
  29.       PARAMETER (AMPLIT=2.D0*XPI)
  30. *     max d'excentrement admissible  et min pour detecter 0
  31.       PARAMETER (EMAX=0.9999D0)
  32.       PARAMETER (EMIN=1.0D-16)
  33. *
  34. **********************************************************************
  35. *
  36. *     coef de P0
  37.       PCOEF = ((0.5D0*XLONG)**2) * 3.D0 * VISCDY
  38.       PALMAD = PALM / PCOEF
  39.  
  40. * --- calcul de E, phi(COSPHI et SINPHI), VE et VPHI
  41. *
  42.       E     = SQRT(XAD**2 + YAD**2)
  43. *    -jeu entierement consomme
  44.       IF(E.GE.EMAX) THEN
  45.         write(ioimp,*) 'Excentrement relatif=',E,' trop proche de 1 :',
  46.      &                  'on le limite a ',EMAX
  47.         E=EMAX
  48.       ENDIF
  49. *
  50. * --- Boucle sur les angles et calcul de la pression
  51. *
  52.       DTETA=AMPLIT/DBLE(NBMAIL)
  53.       FXAD = 0.D0
  54.       FYAD = 0.D0
  55.  
  56.       DO 1 I=1,NBMAIL
  57. *
  58. * -----  Calcul de l'angle et epaisseur du film
  59. c          TETA  = DBLE(I-1)*DTETA
  60. c          COSTE = COS(TETA)
  61. c          SINTE = SIN(TETA)
  62.          COSTE = XPROG(1,I)
  63.          SINTE = XPROG(2,I)
  64.          H     = XJEU - X*COSTE - Y*SINTE
  65.  
  66. * -----  Calcul de variables auxiliaires
  67.          H3     = H**3
  68.          DHDTET =  X*SINTE -  Y*COSTE
  69.          DHDT   = VX*COSTE + VY*SINTE
  70. *
  71. * -----  Calcul de la pression en Z=0
  72.          P0  = ( 2.D0*DHDT - VITROT*DHDTET) / H3
  73.          IF(PALM.NE.0.D0) P0 = PALMAD + P0
  74. *        si condition de Gumbel on prends partie positive de p
  75.          IF(ICAVIT.EQ.2) P0 = MAX(0.D0,P0)
  76.  
  77. * ----- On integre la force avec les trapezes
  78. *       simplification car fonction periodique :
  79. *       p0*cos(1)=p0*cos(NBMAIL+1)  et  idem pour p0*sin
  80.         FXAD = FXAD + (P0 * COSTE)
  81.         FYAD = FYAD + (P0 * SINTE)
  82.  
  83.  1    CONTINUE
  84. *
  85. *
  86. * --- Calcul des efforts (valeur dimensionnee)
  87. *
  88. *     on multiplie par l integrale selon Z du profil parabolique
  89. *     et le R du R dtheta
  90.       ZCOEF = -2.D0*XLONG/3.D0 * DTETA * RARBRE
  91.  
  92.       FX = ZCOEF * PCOEF * FXAD
  93.       FY = ZCOEF * PCOEF * FYAD
  94. *
  95.       RETURN
  96.       END
  97.  
  98.  
  99.  
  100.  

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