Télécharger gflex2.eso

Retour à la liste

Numérotation des lignes :

  1. C GFLEX2 SOURCE CB215821 17/11/30 21:16:18 9639
  2. C
  3. C GFLEX2 SOURCE VACE 25/11/88
  4. C
  5. SUBROUTINE GFLEX2(AB,DLL,RF,CTC,DELTAT,LANBN,FREQ1,FREQ2)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. integer p
  9. C
  10. C
  11. C =====================================================================
  12. C MODELE DE BERNOUILLI-EULER
  13. C --------------------------
  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)/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)/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
  40. C
  41. C CREATION : 25/11/88
  42. C PROGRAMMEUR : VACELET
  43. C =====================================================================
  44. C
  45. C 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.  
  50. -INC PPARAM
  51. -INC CCOPTIO
  52. -INC CCREEL
  53. DIMENSION AB(10,*)
  54. C
  55. C---- DEFINITION DE QUELQUES CONSTANTES
  56. AK1=SQRT(2.*XPI*RF*FREQ1/CTC)
  57. AK2=SQRT(2.*XPI*RF*FREQ2/CTC)
  58. XI=DLL/RF
  59. U1=AK1**2
  60. U2=AK2**2
  61. AA=XPI/(16.*AK2)
  62. BB=XPI/(8.*XI)
  63. DT=CTC*DELTAT/RF
  64. DTPI=DT/XPI
  65. C
  66. C=====================================C
  67. C BOUCLE SUR LES PAS DE TEMPS C
  68. C=====================================C
  69. C
  70. DO 40 L=1,LANBN
  71. T=L*DELTAT*CTC/RF
  72. C
  73. C 1- CALCUL DES FONCTIONS DE GREEN POUR X=0
  74. C ======================================
  75. C
  76. C------ CALCUL DU PAS D INTEGRATION DK ET DU NOMBRE DE PAS NP0
  77. DK= MIN(AA/T,(AK2-AK1)/100)
  78. NP0=INT((AK2-AK1)/DK)
  79. C
  80. C------ INITIALISATION DES VARIABLES GREEI
  81. GREE0 = .5D0*DK*SIN(U1*T)/U1
  82. GREE2 = U1*GREE0
  83. GREE4 = U1*GREE2
  84. C
  85. C------ INTEGRATION SUR K
  86. C
  87. DO 10 P=1,NP0-1
  88. AK=AK1+P*DK
  89. U=AK*AK
  90. UT=U*T
  91. RESU=DK*SIN(UT)
  92. GREE0 = GREE0 + RESU/U
  93. GREE2 = GREE2 + RESU
  94. GREE4 = GREE4 + RESU*U
  95. 10 CONTINUE
  96. C
  97. RESU2=.5*DK*SIN(U2*T)
  98. GREE0 = GREE0 + RESU2/U2
  99. GREE2 = GREE2 + RESU2
  100. GREE4 = GREE4 + RESU2*U2
  101. C
  102. C
  103. C
  104. C 2- CALCUL DES FONCTIONS DE GREEN POUR X=L
  105. C ======================================
  106. C
  107. C------ CALCUL DU PAS DK0 ET DU NOMBRE DE PAS NPL
  108. DK0=MIN(AA/T,BB)
  109. DK0= MIN (DK0,(AK2-AK1)/100)
  110. NPL=INT((AK2-AK1)/DK0)
  111. C
  112. C------ INITIALISATION DES GREELI
  113. DKL=DK0/2.
  114. AL1=DKL*SIN(U1*T)
  115. BL1=AL1*COS(AK1*XI)
  116. BL2=AL1*SIN(AK1*XI)
  117. GREEL0 = BL1/U1
  118. GREEL1 = BL2/AK1
  119. GREEL2 = BL1
  120. GREEL3 = BL2*AK1
  121. GREEL4 = BL1*U1
  122. C
  123. C------ INTEGRATION SUR K
  124. C
  125. DO 20 P=1,NPL-1
  126. AK=AK1+P*DK0
  127. U=AK**2
  128. AV=U*T
  129. AW=AK*XI
  130. AY=AV+AW
  131. AZ=AV-AW
  132. SIY=SIN(AY)
  133. SIZ=SIN(AZ)
  134. CIY=COS(AY)
  135. CIZ=COS(AZ)
  136. C
  137. GREEL0 = GREEL0 + DKL*(SIY+SIZ)/U
  138. GREEL1 = GREEL1 + DKL*(CIZ-CIY)/AK
  139. GREEL2 = GREEL2 + DKL*(SIY+SIZ)
  140. GREEL3 = GREEL3 + DKL*(CIZ-CIY)*AK
  141. GREEL4 = GREEL4 + DKL*(SIY+SIZ)*U
  142. 20 CONTINUE
  143. C
  144. CL1=DKL*SIN(U2*T)
  145. DL1=CL1*COS(AK2*XI)
  146. DL2=CL1*SIN(AK2*XI)
  147. C
  148. GREEL0 = GREEL0 + DL1/U2
  149. GREEL1 = GREEL1 + DL2/AK2
  150. GREEL2 = GREEL2 + DL1
  151. GREEL3 = GREEL3 + DL2*AK2
  152. GREEL4 = GREEL4 + DL1*U2
  153. C
  154. C
  155. C 3- CALCUL DES AB
  156. C =============
  157. C
  158. AB(1,L) = DTPI*GREE0
  159. AB(2,L) = 0.D0
  160. AB(3,L) = -DTPI*GREE2
  161. AB(4,L) = 0.D0
  162. AB(5,L) = DTPI*GREE4
  163. AB(6,L) = DTPI*GREEL0
  164. AB(7,L) = -DTPI*GREEL1
  165. AB(8,L) = -DTPI*GREEL2
  166. AB(9,L) = DTPI*GREEL3
  167. AB(10,L) = DTPI*GREEL4
  168. C
  169. 40 CONTINUE
  170. C
  171. RETURN
  172. C
  173. END
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  

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