Télécharger hbmzfx.eso

Retour à la liste

Numérotation des lignes :

hbmzfx
  1. C HBMZFX SOURCE BP208322 20/09/18 21:16:57 10718
  2.  
  3. SUBROUTINE HBMZFX(NT,NDDL,OMEG,Q1,AiDi,Di,LAMBD,VCTCS,ZX)
  4.  
  5. *=======================================================================
  6. * Rajoute les forces de couplage Fluide-Elastique : ZX=ZX+Ff = Zf*X
  7. *=======================================================================
  8.  
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8(A-H,O-Z)
  11. INTEGER NHBM,NT,NDDL,H
  12. REAL*8 ZX(NT),Q1(NT),VCTCS(7),LAMBD(NDDL),AiDi(2),Di(2),bi(2)
  13. REAL*8 OMEG,V,cf0,mf0,sf0,dmfej,cfei,dniR,dniI
  14.  
  15. *-----------------------------------------------------------------------
  16. * 0. Recuperation des coefficients
  17. V = VCTCS(4)
  18. * Coefficient de couplage: mf0 = -(rho_f/2)*(dCLdy)*V^2
  19. mf0 = -VCTCS(1)*VCTCS(7)*(V**2)
  20. * Coefficient d'amortissement: cf0 = (rho_f/2)*(DIAM)*(CD)*V
  21. cf0 = VCTCS(1)*VCTCS(5)*VCTCS(6)*V
  22. * Coefficient de raideur: sf0 = mf0*alpha_0
  23. sf0 = mf0*VCTCS(3)
  24. DO I = 1,2
  25. bi(I) = Di(I)*V/VCTCS(5)
  26. ENDDO
  27. NHBM = (NT/NDDL - 1)/2
  28.  
  29. *-----------------------------------------------------------------------
  30. * 1. Termes de raideur/amortissement ajoutes
  31. * 1.1 Composante statique
  32. DO I=1,NDDL
  33. ZX(I) = ZX(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. ZX(NDDL*(1+(J-2))+I) = ZX(NDDL*(1+(J-2))+I)
  41. & + sf0 *LAMBD(I)*Q1(NDDL*(1+(J-2))+I)
  42. & + dmfej*LAMBD(I)*Q1(NDDL*(1+(J-1))+I)
  43. * Sin
  44. ZX(NDDL*(1+(J-1))+I) = ZX(NDDL*(1+(J-1))+I)
  45. & + sf0 *LAMBD(I)*Q1(NDDL*(1+(J-1))+I)
  46. & - dmfej*LAMBD(I)*Q1(NDDL*(1+(J-2))+I)
  47. ENDDO
  48. ENDDO
  49.  
  50. *-----------------------------------------------------------------------
  51. ** 2. Effet memoire
  52. ** Boucle sur le nombre de termes
  53. DO JJ = 1,2
  54. cfei = mf0*AiDi(JJ)*V/VCTCS(5)
  55. * 2.1 Composante statique
  56. DO I=1,NDDL
  57. ZX(I) = ZX(I) + cfei*LAMBD(I)*Q1(I)/bi(JJ)
  58. ENDDO
  59. * 2.2 Boucle sur les harmoniques >0
  60. DO H = 2,2*NHBM,2
  61. dniR = bi(JJ)/((bi(JJ))**2+((H/2)*OMEG)**2)
  62. dniI = OMEG*(H/2)/((bi(JJ))**2+((H/2)*OMEG)**2)
  63. ** Boucle sur les DDL
  64. DO K = 1,NDDL
  65. * Cos
  66. ZX(NDDL*(1+(H-2))+K) = ZX(NDDL*(1+(H-2))+K)
  67. & + LAMBD(K)*cfei* ( dniR*Q1(NDDL*(1+(H-2))+K)
  68. & -dniI*Q1(NDDL*(1+(H-1))+K) )
  69. * Sin
  70. ZX(NDDL*(1+(H-1))+K) = ZX(NDDL*(1+(H-1))+K)
  71. & + LAMBD(K)*cfei*( dniR*Q1(NDDL*(1+(H-1))+K)
  72. & +dniI*Q1(NDDL*(1+(H-2))+K) )
  73. ENDDO
  74. ENDDO
  75. ENDDO
  76.  
  77. END
  78.  
  79.  
  80.  
  81.  

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