Télécharger palier.eso

Retour à la liste

Numérotation des lignes :

  1. C PALIER SOURCE BP208322 15/07/22 21:15:35 8586
  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)
  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 P0C(NBMAIL),P0S(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.  
  52. DO 1 I=1,NBMAIL
  53. *
  54. * ----- Calcul de l'angle et epaisseur du film
  55. TETA = DBLE(I-1)*DTETA
  56. COSTE = COS(TETA)
  57. SINTE = SIN(TETA)
  58. H = XJEU - X*COSTE - Y*SINTE
  59.  
  60. * ----- Calcul de variables auxiliaires
  61. H3 = H**3
  62. DHDTET = X*SINTE - Y*COSTE
  63. DHDT = VX*COSTE + VY*SINTE
  64. *
  65. * ----- Calcul de la pression en Z=0
  66. P0 = ( 2.D0*DHDT - VITROT*DHDTET) / H3
  67. IF(PALM.NE.0.D0) P0 = PALMAD + P0
  68. * si condition de Gumbel on prends partie positive de p
  69. IF(ICAVIT.EQ.2) P0 = MAX(0.D0,P0)
  70. P0C(I)= P0 * COSTE
  71. P0S(I)= P0 * SINTE
  72.  
  73. 1 CONTINUE
  74. *
  75. * --- On integre avec les trapezes
  76. *
  77. FXAD = 0.D0
  78. FYAD = 0.D0
  79. c simplification car fonction periodique --> p0c(1)=p0c(NBMAIL+1)
  80. DO 2 I=1,NBMAIL
  81. FXAD = FXAD + P0C(I)
  82. FYAD = FYAD + P0S(I)
  83. 2 CONTINUE
  84. *
  85. * --- Calcul des efforts (valeur dimensionnee)
  86. *
  87. * on multiplie par l integrale selon Z du profil parabolique
  88. * et le R du R dtheta
  89. ZCOEF = -2.D0*XLONG/3.D0 * DTETA * RARBRE
  90.  
  91. FX = ZCOEF * PCOEF * FXAD
  92. FY = ZCOEF * PCOEF * FYAD
  93. *
  94. RETURN
  95. END
  96.  
  97.  
  98.  

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