Télécharger la plaquette Cast3M


fr en

Télécharger hbmrwf.eso

Retour à la liste

Numérotation des lignes : 

hbmrwf

  1. C HBMRWF    SOURCE    CB215821  26/06/25    21:15:12     12581          
  2.  
  3. *=======================================================================
  4. * Rajoute la derivee des efforts de couplage fluidelastique:
  5. * R,w = R,w - dFf/dw
  6. *     = R,w - Zwf*X
  7. *=======================================================================
  8.  
  9.       SUBROUTINE HBMRWF(NT,NDDL,OMEG,Q1,AiDi,Di,LAMBD,VCTCS,Rw)
  10.  
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8(A-H,O-Z)
  13.  
  14.       INTEGER NT,NDDL
  15.       REAL*8 OMEG
  16. C      DIMENSION Rw(NT),Q1(NT),VCTCS(7),LAMBD(NDDL),AiDi(2),Di(2)
  17.       REAL*8 Rw(*),Q1(*),VCTCS(*),LAMBD(*),AiDi(*),Di(*)
  18.  
  19.       INTEGER NHBM,H
  20.       REAL*8 V,cf0,mf0,dmfej,cfei,cii,dniR,dniI
  21.       REAL*8 bi(2)
  22.  
  23. *-----------------------------------------------------------------------
  24. * 0. Recuperation des coefficients
  25.       V = VCTCS(4)
  26. *     Coefficient de couplage: mf0dV = -(rho_f/2)*(dCLdy)*V^2
  27.       mf0  = -VCTCS(1)*VCTCS(7)*(V**2)
  28. *     Coefficient d'amortissement: cf0 = (rho_f/2)*(DIAM)*(CD)*V
  29.       cf0  = VCTCS(1)*VCTCS(5)*VCTCS(6)*V
  30.       DO I = 1,2
  31.         bi(I) = Di(I)*V/VCTCS(5)
  32.       ENDDO
  33.       NHBM = (NT/NDDL - 1)/2
  34.  
  35. *-----------------------------------------------------------------------
  36. * 1. Terme d'amortissement ajoute
  37. *  1.1 Composante statique: nulle
  38. *  1.2 Boucle sur les harmoniques
  39.       DO J=2,2*NHBM,2
  40.         dmfej =REAL(J/2)*cf0
  41.         DO I=1,NDDL
  42. *         Cos
  43.           Rw(NDDL*(1+(J-2))+I) = Rw(NDDL*(1+(J-2))+I)
  44.      &                         + dmfej*LAMBD(I)*Q1(NDDL*(1+(J-1))+I)
  45. *         Sin
  46.           Rw(NDDL*(1+(J-1))+I) = Rw(NDDL*(1+(J-1))+I)
  47.      &                         - dmfej*LAMBD(I)*Q1(NDDL*(1+(J-2))+I)
  48.         ENDDO
  49.       ENDDO
  50.  
  51. *-----------------------------------------------------------------------
  52. ** 2. Effet memoire
  53. ** Boucle sur le nombre de termes
  54.       DO JJ = 1,2
  55.         cfei = mf0*AiDi(JJ)*V/VCTCS(5)
  56. *   2.1 Composante statique: nulle
  57. *   2.2 Boucle sur les harmoniques
  58.         DO H = 2,2*NHBM,2
  59.           xHs2=REAL(H/2)
  60.           cii = 1./((bi(JJ))**2+(xHs2*OMEG)**2)
  61.           dniR = -2.D0*bi(JJ)*OMEG*(xHs2*cii)**2
  62.           dniI = 2.D0*(xHs2**3)*(OMEG*cii)**2 -cii*xHs2
  63. **        Boucle sur les DDL
  64.           DO K = 1,NDDL
  65. *          Cos
  66.             Rw(NDDL*(1+(H-2))+K) = Rw(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.             Rw(NDDL*(1+(H-1))+K) = Rw(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.       RETURN
  78.       END
  79.  
  80.  
  81.  
  82.  

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