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 :
© Cast3M 2003 - All rights reserved.
Disclaimer