Télécharger palier.eso

Retour à la liste

Numérotation des lignes :

  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