Télécharger gflex2.eso

Retour à la liste

Numérotation des lignes :

  1. C GFLEX2 SOURCE PV 05/04/25 21:15:21 5077
  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. -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=====================================C
  65. C BOUCLE SUR LES PAS DE TEMPS C
  66. C=====================================C
  67. C
  68. DO 40 L=1,LANBN
  69. T=L*DELTAT*CTC/RF
  70. C
  71. C 1- CALCUL DES FONCTIONS DE GREEN POUR X=0
  72. C ======================================
  73. C
  74. C------ CALCUL DU PAS D INTEGRATION DK ET DU NOMBRE DE PAS NP0
  75. DK= MIN(AA/T,(AK2-AK1)/100)
  76. NP0=IDINT((AK2-AK1)/DK)
  77. C
  78. C------ INITIALISATION DES VARIABLES GREEI
  79. GREE0 = .5D0*DK*SIN(U1*T)/U1
  80. GREE2 = U1*GREE0
  81. GREE4 = U1*GREE2
  82. C
  83. C------ INTEGRATION SUR K
  84. C
  85. DO 10 P=1,NP0-1
  86. AK=AK1+P*DK
  87. U=AK*AK
  88. UT=U*T
  89. RESU=DK*SIN(UT)
  90. GREE0 = GREE0 + RESU/U
  91. GREE2 = GREE2 + RESU
  92. GREE4 = GREE4 + RESU*U
  93. 10 CONTINUE
  94. C
  95. RESU2=.5*DK*SIN(U2*T)
  96. GREE0 = GREE0 + RESU2/U2
  97. GREE2 = GREE2 + RESU2
  98. GREE4 = GREE4 + RESU2*U2
  99. C
  100. C
  101. C
  102. C 2- CALCUL DES FONCTIONS DE GREEN POUR X=L
  103. C ======================================
  104. C
  105. C------ CALCUL DU PAS DK0 ET DU NOMBRE DE PAS NPL
  106. DK0=MIN(AA/T,BB)
  107. DK0= MIN (DK0,(AK2-AK1)/100)
  108. NPL=IDINT((AK2-AK1)/DK0)
  109. C
  110. C------ INITIALISATION DES GREELI
  111. DKL=DK0/2.
  112. AL1=DKL*SIN(U1*T)
  113. BL1=AL1*COS(AK1*XI)
  114. BL2=AL1*SIN(AK1*XI)
  115. GREEL0 = BL1/U1
  116. GREEL1 = BL2/AK1
  117. GREEL2 = BL1
  118. GREEL3 = BL2*AK1
  119. GREEL4 = BL1*U1
  120. C
  121. C------ INTEGRATION SUR K
  122. C
  123. DO 20 P=1,NPL-1
  124. AK=AK1+P*DK0
  125. U=AK**2
  126. AV=U*T
  127. AW=AK*XI
  128. AY=AV+AW
  129. AZ=AV-AW
  130. SIY=SIN(AY)
  131. SIZ=SIN(AZ)
  132. CIY=COS(AY)
  133. CIZ=COS(AZ)
  134. C
  135. GREEL0 = GREEL0 + DKL*(SIY+SIZ)/U
  136. GREEL1 = GREEL1 + DKL*(CIZ-CIY)/AK
  137. GREEL2 = GREEL2 + DKL*(SIY+SIZ)
  138. GREEL3 = GREEL3 + DKL*(CIZ-CIY)*AK
  139. GREEL4 = GREEL4 + DKL*(SIY+SIZ)*U
  140. 20 CONTINUE
  141. C
  142. CL1=DKL*SIN(U2*T)
  143. DL1=CL1*COS(AK2*XI)
  144. DL2=CL1*SIN(AK2*XI)
  145. C
  146. GREEL0 = GREEL0 + DL1/U2
  147. GREEL1 = GREEL1 + DL2/AK2
  148. GREEL2 = GREEL2 + DL1
  149. GREEL3 = GREEL3 + DL2*AK2
  150. GREEL4 = GREEL4 + DL1*U2
  151. C
  152. C
  153. C 3- CALCUL DES AB
  154. C =============
  155. C
  156. AB(1,L) = DTPI*GREE0
  157. AB(2,L) = 0.D0
  158. AB(3,L) = -DTPI*GREE2
  159. AB(4,L) = 0.D0
  160. AB(5,L) = DTPI*GREE4
  161. AB(6,L) = DTPI*GREEL0
  162. AB(7,L) = -DTPI*GREEL1
  163. AB(8,L) = -DTPI*GREEL2
  164. AB(9,L) = DTPI*GREEL3
  165. AB(10,L) = DTPI*GREEL4
  166. C
  167. 40 CONTINUE
  168. C
  169. RETURN
  170. C
  171. END
  172.  
  173.  
  174.  
  175.  
  176.  

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