Télécharger hbmrv.eso

Retour à la liste

Numérotation des lignes :

hbmrv
  1. C HBMRV SOURCE BP208322 20/09/18 21:16:49 10718
  2.  
  3. SUBROUTINE HBMRV(NT,NDDL,OMEG,Q1,AiDi,Di,LAMBD,VCTCS,Rv)
  4.  
  5. *=======================================================================
  6. * Derivee du residu par rapport a la vitesse moyenne d'ecoulement Vmoy:
  7. * Rv = Zv*X
  8. * (pour a continuation en fonction de Vmoy)
  9. *=======================================================================
  10.  
  11. IMPLICIT INTEGER(H-N)
  12. IMPLICIT REAL*8(A-G,O-Z)
  13. REAL*8 OMEG,mf0,cf0,sf0
  14. REAL*8 Rv(NT),Q1(NT),VCTCS(7),LAMBD(NDDL),AiDi(2),Di(2),bi(2)
  15.  
  16. *-----------------------------------------------------------------------
  17. * 0. Recuperation des coefficients
  18. V = VCTCS(4)
  19. * Coefficient de couplage: mf0dV = -(rho_f/2)*(dCLdy)*V^2
  20. mf0 = -VCTCS(1)*VCTCS(7)*(V**2)
  21. * Coefficient d'amortissement: cf0 = (rho_f/2)*(DIAM)*(CD)
  22. cf0 = VCTCS(1)*VCTCS(5)*VCTCS(6)
  23. DO I = 1,2
  24. bi(I) = Di(I)*V/VCTCS(5)
  25. ENDDO
  26. sf0 = 2.D0*mf0*VCTCS(3)/V
  27. NHBM = (NT/NDDL - 1)/2
  28.  
  29. *-----------------------------------------------------------------------
  30. * 1. Terme de raideur/amortissement ajoutes
  31. * 1.1 Composante statique:
  32. DO I=1,NDDL
  33. Rv(I) = sf0*LAMBD(I)*Q1(I)
  34. ENDDO
  35. * 1.2 Boucle sur les harmoniques
  36. DO J=2,2*NHBM,2
  37. dmfej = (J/2)*OMEG*cf0
  38. DO I=1,NDDL
  39. * Cos
  40. Rv(NDDL*(1+(J-2))+I) = sf0*LAMBD(I)*Q1(NDDL*(1+(J-2))+I) +
  41. & dmfej*LAMBD(I)*Q1(NDDL*(1+(J-1))+I)
  42. * Sin
  43. Rv(NDDL*(1+(J-1))+I) = sf0*LAMBD(I)*Q1(NDDL*(1+(J-1))+I) -
  44. & dmfej*LAMBD(I)*Q1(NDDL*(1+(J-2))+I)
  45. ENDDO
  46. ENDDO
  47.  
  48. *-----------------------------------------------------------------------
  49. ** 2. Effet memoire
  50. ** Boucle sur le nombre de termes
  51. DO JJ = 1,2
  52. cfei = mf0*AiDi(JJ)*V/VCTCS(5)
  53. * 2.1 Composante statique:
  54. DO I = 1,NDDL
  55. Rv(I) = Rv(I)+2.*cfei*LAMBD(I)*Q1(I)/(bi(JJ)*V)
  56. ENDDO
  57. * 2.2 Boucle sur les harmoniques
  58. DO H = 2,2*NHBM,2
  59. cii = 1.D0/((bi(JJ))**2+((H/2)*OMEG)**2)
  60. dniR = bi(JJ)*(-2.D0*(bi(JJ)*cii)**2 + 4.D0*cii)/V
  61. dniI = (H/2)*OMEG*(2.D0*(cii*bi(JJ))**2 - 3.D0*cii)/V
  62. ** Boucle sur les DDL
  63. DO K = 1,NDDL
  64. * Cos
  65. Rv(NDDL*(1+(H-2))+K)=LAMBD(K)*cfei*(dniR*Q1(NDDL*(1+(H-2))+K)
  66. & +dniI*Q1(NDDL*(1+(H-1))+K))+Rv(NDDL*(1+(H-2))+K)
  67. * Sin
  68. Rv(NDDL*(1+(H-1))+K)=LAMBD(K)*cfei*(dniR*Q1(NDDL*(1+(H-1))+K)
  69. & -dniI*Q1(NDDL*(1+(H-2))+K))+Rv(NDDL*(1+(H-1))+K)
  70. ENDDO
  71. ENDDO
  72. ENDDO
  73. END
  74.  
  75.  
  76.  

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