Télécharger gflex3.eso

Retour à la liste

Numérotation des lignes :

  1. C GFLEX3 SOURCE CB215821 17/11/30 21:16:19 9639
  2. C
  3. C GFLEX3 SOURCE GAUT 25/11/88
  4. C
  5. SUBROUTINE GFLEX3(AB,DLL,RF,CTC,DELTAT,LANBN,FREQ1,FREQ2,EPS)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. integer p
  9. C
  10. C
  11. C =====================================================================
  12. C
  13. C MODELE DE BERNOUILLI-EULER AVEC AMORTISSEMENT PHYSIQUE
  14. C ------------------------------------------------------
  15. C
  16. C SOUS-PROGRAMME FORTRAN APPELE PAR GREEN1 POUR CALCULER LES FCTS DE
  17. C GREEN FILTREES EN FLEXION PAR INTEGRATION PAR LA METHODE DES TRAPEZES
  18. C
  19. C LA FONCTION DE GREEN EST DEFINI PAR :
  20. C -------------------------------------
  21. C G(O,T) = INTEGRALE SUR K ENTRE AK1 ET AK2 DE :
  22. C SIN(K*K*T)*EXP(-EPS*K**2*T)/K*K*XPI
  23. C G(L,T) = INTEGRALE SUR K ENTRE AK1 ET AK2 DE :
  24. C SIN(K*K*T)*COS(K*XI)*EXP(-EPS*K**2*T)/K*K*XPI
  25. C
  26. C AVEC AK=SQRT(2*XPI*RF*FREQ/CTC) OU FREQ EST LA FREQUENCE
  27. C
  28. C DEFINITION DES ARGUMENTS :
  29. C --------------------------
  30. C AB : TABLEAU 10*LANBN DONNANT LES INTEGRALES DES FONCTIONS
  31. C DE GREEN ET DE LEURS DERIVEES A CHAQUE PAS DE TEMPS
  32. C DLL : LONGUEUR DE LA POUTRE
  33. C RF : DIMENSION TRANSVERSE DE LA POUTRE
  34. C CTC : TEMPS DE CALCUL
  35. C DELTAT : PAS DE TEMPS
  36. C LANBN : NOMBRE DE PAS DE TEMPS
  37. C FREQ1 : FREQUENCE BASSE DU FILTRE
  38. C FREQ2 : FREQUENCE HAUTE DU FILTRE
  39. C EPS : COEFFICIENT D AMORTISSEMENT
  40. C
  41. C CREATION : 10/12/88
  42. C PROGRAMMEUR : VACELET
  43. C =====================================================================
  44. C
  45. C DEFINITION DES VARIABLES UTILISEES :
  46. C GREEI : INTEGRALE DONNANT LA DERIVEE IEME DE GREEN EN 0
  47. C GREELI : INTEGRALE DONNANT LA DERIVEE IEME DE GREEN EN L
  48. C
  49. -INC CCOPTIO
  50. -INC CCREEL
  51. DIMENSION AB(10,*)
  52. C
  53. C---- DEFINITION DE QUELQUES CONSTANTES
  54. AK1=SQRT(2.*XPI*RF*FREQ1/CTC)
  55. AK2=SQRT(2.*XPI*RF*FREQ2/CTC)
  56. XI=DLL/RF
  57. U1=AK1**2
  58. U2=AK2**2
  59. AA=XPI/(16.*AK2)
  60. BB=XPI/(8.*XI)
  61. DT=CTC*DELTAT/RF
  62. DTPI=DT/XPI
  63. C
  64. C
  65. C=====================================C
  66. C BOUCLE SUR LES PAS DE TEMPS C
  67. C=====================================C
  68. C
  69. DO 40 L=1,LANBN
  70. T=L*DELTAT*CTC/RF
  71. C
  72. C 1- CALCUL DES FONCTIONS DE GREEN POUR X=0
  73. C ======================================
  74. C
  75. C------ CALCUL DU PAS D INTEGRATION DK ET DU NOMBRE DE PAS NP0
  76. DK=AA/T
  77. NP0=INT((AK2-AK1)/DK)
  78. C
  79. C------ INITIALISATION DE GREEI
  80. GREE0 = .5D0*DK*SIN(U1*T)*EXP(-EPS*U1*T)/U1
  81. GREE2 = U1*GREE0
  82. GREE4 = U1*GREE2
  83. C
  84. C------ INTEGRATION SUR K
  85. C
  86. DO 10 P=1,NP0-1
  87. AK=AK1+P*DK
  88. U=AK*AK
  89. UT=U*T
  90. RESU=DK*SIN(UT)*EXP(-EPS*U*T)
  91. GREE0 = GREE0 + RESU/U
  92. GREE2 = GREE2 + RESU
  93. GREE4 = GREE4 + RESU*U
  94. 10 CONTINUE
  95. C
  96. RESU2=.5*DK*SIN(U2*T)*EXP(-EPS*U2*T)
  97. GREE0 = GREE0 + RESU2/U2
  98. GREE2 = GREE2 + RESU2
  99. GREE4 = GREE4 + RESU2*U2
  100. C
  101. C
  102. C 2- CALCUL DES FONCTIONS DE GREEN POUR X=L
  103. C ======================================
  104. C
  105. C------ CALCUL DU PAS D INTEGRATION DK0 ET DU NOMBRE DE PAS NPL
  106. DK0=MIN(AA/T,BB)
  107. NPL=INT((AK2-AK1)/DK0)
  108. C
  109. C------ INITIALISATION DE GREELI
  110. DKL=DK0/2.
  111. AL1=DKL*SIN(U1*T)*EXP(-EPS*U1*T)
  112. BL1=AL1*COS(AK1*XI)
  113. BL2=AL1*SIN(AK1*XI)
  114. GREEL0 = BL1/U1
  115. GREEL1 = BL2/AK1
  116. GREEL2 = BL1
  117. GREEL3 = BL2*AK1
  118. GREEL4 = BL1*U1
  119. C
  120. C------ INTEGRATION SUR K
  121. C
  122. DO 20 P=1,NPL-1
  123. AK=AK1+P*DK0
  124. U=AK**2
  125. AV=U*T
  126. AW=AK*XI
  127. AY=AV+AW
  128. AZ=AV-AW
  129. SIY=SIN(AY)
  130. SIZ=SIN(AZ)
  131. CIY=COS(AY)
  132. CIZ=COS(AZ)
  133. DKEXP=DKL*EXP(-EPS*U*T)
  134. C
  135. GREEL0 = GREEL0 + DKEXP*(SIY+SIZ)/U
  136. GREEL1 = GREEL1 + DKEXP*(CIZ-CIY)/AK
  137. GREEL2 = GREEL2 + DKEXP*(SIY+SIZ)
  138. GREEL3 = GREEL3 + DKEXP*(CIZ-CIY)*AK
  139. GREEL4 = GREEL4 + DKEXP*(SIY+SIZ)*U
  140. 20 CONTINUE
  141. C
  142. CL1=DKL*SIN(U2*T)*EXP(-EPS*U2*T)
  143. DL1=CL1*COS(AK2*XI)
  144. DL2=CL1*SIN(AK2*XI)
  145. GREEL0 = GREEL0 + DL1/U2
  146. GREEL1 = GREEL1 + DL2/AK2
  147. GREEL2 = GREEL2 + DL1
  148. GREEL3 = GREEL3 + DL2*AK2
  149. GREEL4 = GREEL4 + DL1*U2
  150. C
  151. C
  152. C 3- CALCUL DES AB
  153. C =============
  154. C
  155. AB(1,L) = DTPI*GREE0
  156. AB(2,L) = 0.D0
  157. AB(3,L) = -DTPI*GREE2
  158. AB(4,L) = 0.D0
  159. AB(5,L) = DTPI*GREE4
  160. AB(6,L) = DTPI*GREEL0
  161. AB(7,L) = -DTPI*GREEL1
  162. AB(8,L) = -DTPI*GREEL2
  163. AB(9,L) = DTPI*GREEL3
  164. AB(10,L) = DTPI*GREEL4
  165. C
  166. 40 CONTINUE
  167. C
  168. RETURN
  169. C
  170. END
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  

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