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

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