1 : C GFLEX3 SOURCE PV 05/04/25 21:15:23 5077 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=IDINT((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=IDINT((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 :
© Cast3M 2003 - All rights reserved.
Disclaimer