1 : * FLAMBAGE PROCEDUR AM 02/01/14 21:15:01 4286 2 : DEBPROC FLAMBAGE ETAB11*TABLE ; 3 : * CETTE PROCEDURE EFFECTUE UN CALCUL DE FLAMBAGE 4 : * 5 : * EN ENTREE : 6 : * --------- 7 : * ETAB TABLE CONTENANT : 8 : * INDICE 'OBJM' : OBJET MODELE 9 : * INDICE 'SIG0' : CHAMP CONSTANT DE CONTRAINTES 10 : * INDICE 'SIG1' : CONTRAINTE DONT ON CHERCHE LE 11 : * COEFFICIENT MULTIPLICATEUR 12 : * INDICE 'RIGI' : MATRICE DE RAIDEUR 13 : * INDICE 'KSIG' : MATRICE DE RIGIDITE GEOMETRIQUE 14 : * INDICE 'LAM1' : COEFFICIENT MINI 15 : * INDICE 'LAM2' : COEFFICIENT MAXI 16 : * INDICE 'NMOD' : NOMBRE DE MODES AZIMUTAL 17 : * INDICE 'MODE' : LISTE DES MODES DE FOURIER 18 : * INDICE 'MATE' : MATERIAU 19 : * INDICE 'CARA' : CARACTERISTIQUES 20 : * INDICE 'CLIM' : CONDITIONS LIMITES 21 : * INDICE 'PLAS' : SI FLAMBAGE PLASTIQUE 22 : * INDICE 'MATP' : MATERIAU PLASTIQUE 23 : * 24 : * EN SORTIE : 25 : * --------- 26 : * STAB11 TABLE CONTENANT AUTANT DE TABLES QUE DE MODES 27 : * DE SORTIE 28 : * INDICE N : NUMERO DU MODE 29 : * 30 : * LA TABLE STAB11 N CONTIENT : 31 : * INDICE 'LAMB' : COEFFICIENT MULTIPLICATEUR 32 : * INDICE 'DEPL' : LE DEPLACEMENT MODAL ( MAXI 1) 33 : * INDICE 'MODN' : LE MODE DE FOURIER 34 : * INDICE 'MODM' : LE MODE AZIMUTHAL 35 : * 36 : *--------------------------------------------------------------------- 37 : * 38 : * 39 : *---------------------------------------------------------------------- 40 : * ON RECUPERE L'INFORMATION 41 : *------------------------------ 42 : * 43 : IERR = 0; 44 : * VERIF DE DONNEES 45 : NBMOD = 1; IFOURI= 0; 46 : 'SI' ('EXISTE' ETAB11 'MODE' ) ; 47 : IFOURI =1; 48 : LMOD = ETAB11.'MODE' ; 49 : NBMOD = 'DIME' LMOD ; 50 : 'FINSI'; 51 : * 52 : 'SI' ('EXISTE' ETAB11 'OBJM'); 53 : OBJMOD = ETAB11.'OBJM' ; 54 : 'SINON'; 'MESSAGE' ' ****ERREUR IL FAUT DONNER L OBJET MODELE '; 55 : 'FINSI'; 56 : * 57 : IKSIG=0; 58 : 'SI' ('EXISTE' ETAB11 'SIG1'); 59 : SIG10= 'COPIER' ETAB11.'SIG1' ; 60 : SIG11 = -1. * SIG10; 'DETR' SIG10; 61 : 'SINON' ;'SI' ( IFOURI 'EGA' 1) ; 62 : IERR = 1 ; 'MESSAGE' ' ***** ERREUR IL FAUT DONNER SIG1 '; 'SINON'; 63 : 'SI' ('EXISTE' ETAB11 'KSIG'); IKSIG=1; KSII0 = ETAB11.'KSIG'; 64 : KSIG1 = -1. * KSII0 ; 'SINON'; 65 : 'MESSAGE' ' ***** ERREUR IL FAUT DONNER KSIG OU SIG1'; 66 : IERR = 1; 67 : 'FINSI'; 68 : 'FINSI'; 69 : 'FINSI'; 70 : 'SI' ('EXISTE' ETAB11 'CLIM'); 71 : KLIM = ETAB11.'CLIM' ; 72 : 'SINON'; 73 : 'MESSAGE' ' **** VOUS N AVEZ PAS DONNE DE CONDITIONS LIMITES '; 74 : 'FINSI'; 75 : IRIGN = 0; 76 : 'SI' ( 'EXISTE' ETAB11 'MATE'); 77 : MAT1 = ETAB11.'MATE' ; 78 : 'SINON'; 'SI' ( IFOURI 'EGA' 1); IERR = 1; 79 : 'MESSAGE' ' ****** ERREUR IL FAUT DONNER MATE'; 'SINON'; 80 : 'SI' ( 'EXISTE' ETAB11 'RIGI' ) ; IRIGN =1; 81 : RIGN = ETAB11.'RIGI' ; 82 : 'SINON'; IERR = 1; 83 : 'MESSAGE' ' ***** ERREUR IL FAUT DONNER MATE OU RIGI'; 84 : 'FINSI'; 85 : 'FINSI'; 86 : 'FINSI'; 87 : ICARA = 0; 88 : CAR1 = 'TEXTE' ' ' ; 89 : 'SI' ( 'EXISTE' ETAB11 'CARA'); 90 : ICARA = 1; 91 : CAR1 = ETAB11.'CARA'; 92 : 'FINSI'; 93 : X1 = 1.E30; 94 : 'SI' ( 'EXISTE' ETAB11 'LAM2'); 95 : X1 = ETAB11.'LAM2' ; 96 : 'FINSI' ; 97 : X0 = 0.; 98 : 'SI' ( 'EXISTE' ETAB11 'LAM1'); 99 : X0 = ETAB11.'LAM1' ; 100 : 'FINSI' ; 101 : NMO = 1;'SI' ('EXISTE' ETAB11 'NMOD') ; 102 : NMO = ETAB11.'NMOD'; 'FINSI'; 103 : MATN = MAT1; 104 : * 105 : * LISTE DES MODES DE FOURIER 106 : * 107 : IFLPLA= 0; 108 : 'SI' ('EXISTE' ETAB11 'PLAS' ) ; 109 : IFLPLA = 1; 110 : 'SI' ('EXISTE' ETAB11 'MATP' ); 111 : MATN = ETAB11.'MATP'; 112 : 'SINON'; 113 : IERR = 1; 114 : 'MESSAGE' ' ***** IL FAUT DONNER LE MATERIAU PLASTIQUE ' ; 115 : 'FINSI'; 116 : 'FINSI'; 117 : 'SI'(IERR 'EGA' 0 ) ; 118 : * 119 : *-------------------------------------------- 120 : * INITIALISATIONS 121 : *-------------------------------------------- 122 : PI = 3.14159 ; 123 : XX = 2. * PI ; 124 : * 125 : *-------------------------------------------- 126 : * CALCUL DES FREQUENCES A PARTIR DES COEFFICIENTS MULTIPLICATEURS 127 : *-------------------------------------------- 128 : XX0 = X0 ** .5 ; XX1 = X1 ** .5; 129 : XX0= XX0 / XX ; XX1 = XX1 / XX ; 130 : * 131 : *--------------------------------------------- 132 : * BOUCLE SUR LES MODES DE FOURIER 133 : *--------------------------------------------- 134 : * 135 : TAMON = 'TABLE'; 136 : TAMOM = 'TABLE'; 137 : TACOE = 'TABLE'; 138 : TAPERM = 'TABLE'; 139 : STAB11 = 'TABLE' ; 140 : II = 0 ; 141 : IK = 0 ; 142 : MODMIN= -10; 143 : COEMIN=1.E30; 144 : * 145 : * STAB11 EST UNE TABLE DE TABLES INDICEE PAR LE NUMERO DU MODE SORTI 146 : * 147 : IFOU=0 ; 148 : 'REPETER' BOUCFOU; 149 : IFOU = IFOU + 1 ; 150 : 'SI' ( IFOU > NBMOD) ; 151 : 'QUITTER' BOUCFOU; 152 : 'FINSI' ; 153 : N = 0; 154 : 'SI' (IFOURI 'EGA' 1) ; 155 : N = 'EXTR' LMOD IFOU; 156 : 'OPTION' 'MODE' 'FOUR' N ; 157 : 'FINSI'; 158 : 'SI' (IRIGN 'EGA' 0) ; 'SI' (ICARA 'EGA' 1) ;MATN1=MATN 'ET' CAR1; 159 : 'SINON' ;MATN1=MATN ;'FINSI' ; RIGN = 'RIGI' OBJMOD MATN1;'FINSI'; 160 : TST = 'EXISTE' ETAB11 'SIG0'; 161 : 'SI' ( TST ) ; 162 : SIG10 = ETAB11.'SIG0'; 163 : KSIG0 = 'KSIGMA' OBJMOD SIG10 CAR1 FLAM ; 164 : RIGN = RIGN 'ET' KSIG0 ; 165 : 'FINSI' ; 166 : 'SI' (IKSIG 'EGA' 0 ) ; KSIG1 = 'KSIG' OBJMOD SIG11 CAR1 FLAM ; 167 : 'FINSI'; 168 : 'SI' ( 'EXISTE' ETAB11 'CLIM' ) ; RIGP = RIGN 'ET' KLIM; 169 : 'SINON';RIGP = RIGN ; 170 : 'FINSI'; 171 : *------------------------------------------------------------ 172 : * CALCUL DES FREQUENCES 173 : *----------------------------------------------------------- 174 : XMOD = 'VIBR' 'INTERVALLE' XX0 XX1 'BASSE' NMO RIGP KSIG1 ; 175 : *----------------------------------------------------------- 176 : * COMBIEN A T ON DE MODES 177 : *------------------------------------------------------- 178 : NMOCAL = 'DIMENSION' XMOD; 179 : 'SI' (NMOCAL 'EGA' 0); 180 : 'MESS' ' IL N''Y A PAS DE FREQUENCE DANS L''INTERVALLE CHOISI'; 181 : 'MESS' ' POUR LE MODE DE FOURIER ' N ; 182 : 'SINON'; 183 : *------------------------------ 184 : * EXTRACTION DES MODES 185 : *------------------------------ 186 : NMOO = NMO; 187 : 'SI' ( NMOCAL < NMOO ) ; 188 : 'MESS' ' ***** IL N Y A QUE ' NMOCAL ' MODES DE FLAMBAGE ENTRE ' X0 189 : ' ET ' X1 ; 190 : NMOO= NMOCAL; 191 : 'FINSI'; 192 : JJ = 0; 193 : 'REPETER' BOUCMO; 194 : JJ = JJ + 1; 195 : 'SI' ( JJ > NMOO) ; 196 : 'QUITTER' BOUCMO; 197 : 'FINSI'; 198 : II = II + 1; 199 : FF = 'TIRE' XMOD 'FREQ' RANG JJ ; 200 : DD = 'TIRE' XMOD 'DEPL' RANG JJ ; 201 : XMUL = XX * FF; XMUL = XMUL * XMUL ; 202 : TAMON.II = N ; 203 : TAMOM.II = JJ ; 204 : TACOE.II = XMUL; 205 : * 206 : ****** RECHERCHE DU MODE A COEFFICIENT MINI 207 : * 208 : 'SI'(IFOURI 'EGA' 1 ) ; 209 : 'SI'( XMUL < COEMIN ) ; 210 : MODMIN = N ; 211 : COEMIN = XMUL; 212 : 'FINS'; 213 : 'FINSI'; 214 : * 215 : *-------------------------------------- 216 : * ARCHIVAGE DES RESULTATS SI BESOIN 217 : *-------------------------------------- 218 : STN = 'TABLE'; 219 : STN.'LAMB' = XMUL; 220 : STN.'DEPL' = DD; 221 : STN.'MODN' = N ; 222 : STN.'MODM' = JJ ; 223 : IK = IK + 1; 224 : STAB11.IK = 'TABLE' STN ; 225 : 'FIN' BOUCMO; 226 : 'SI' ( IKSIG 'EGA' 0) ;'DETR' KSIG1 'ELEM';'FINSI'; 227 : 'SI'(IRIGN 'EGA' 0 ) ; 'DETR' RIGN 'ELEM'; 'FINSI'; 228 : 'FINSI'; 229 : 'FIN' BOUCFOU; 230 : * 231 : ********** RANGEMENT PAR ORDRE CROISSANT 232 : * 233 : * TAPERM.1 = 1; 234 : * SI ( II > 1); 235 : *KK = 0; 236 : *REPETER BOUCPER; 237 : * KK = KK + 1; 238 : *SI(KK > II ) ; QUITTER BOUCPER; FINSI; 239 : *XMU = TACOE.KK ; 240 : *III = 0; KKK = 0; 241 : *REPETER BOUCIN; 242 : *III = III + 1; SI (III > II) ; QUITTER BOUCIN; FINSI; 243 : *XMU1 = TACOE.III; 244 : *SI( XMU >EG XMU1); KKK = KKK + 1; FINSI; 245 : *FIN BOUCIN; 246 : *TAPERM.KK = KKK; 247 : *FIN BOUCPER; 248 : * FINSI; 249 : * 250 : * IMPRESSSION DES RESULTATS 251 : 'MESS' ' ***********************************************************'; 252 : 'MESS' ' * MODE FOURIER * MODE M * COEFFICIENT MULTIPLICATEUR *'; 253 : 'MESS' ' ***********************************************************'; 254 : KK = 0; 255 : 'REPETER' BOUC; 256 : KK = KK + 1; 257 : 'SI' (KK > II ) ; 'QUITTER' BOUC ; 'FINSI' ; 258 : * KKK = TAPERM.KK; 259 : KKK = KK ; 260 : N = TAMON.KKK; 261 : XMUL = TACOE.KKK; 262 : JJ = TAMOM.KKK; 263 : 'MESSAGE' ' * ' N ' * ' JJ '* ' XMUL ' *'; 264 : 'FIN' BOUC; 265 : 'MESS' ' ***********************************************************'; 266 : 'SI' ( IFOURI 'EGA' 1) ; 267 : MESSAGE ' '; 268 : MESSAGE ' '; 269 : 'MESS' ' ***** COEFFICIENT MULTIPLICATEUR MINI ' COEMIN ; 270 : 'MESS' ' SUR LE MODE DE FOURIER ' MODMIN ; 271 : 'FINSI'; 272 : 'FINSI'; 273 : *---------------------------------------------------------------------- 274 : * SORTIE DE LA PROCEDURE 275 : *---------------------------------------------------------------------- 276 : 'FINPROC' STAB11 ; 277 :
© Cast3M 2003 - All rights reserved.
Disclaimer