Télécharger la plaquette Cast3M


fr en

Télécharger hbmrv.eso

Retour à la liste

Numérotation des lignes : 

hbmrv

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

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