1 : * TRANSFER PROCEDUR DC 06/02/06 21:16:07 5304 2 : * TRANSFER PROCEDUR JEANV 97/10/27 21:16:57 2887 3 : 'DEBPROC' TRANSFER MOD/SOLUTION TA/TABLE LAMOR*LISTREEL MOT1*MOT; 4 : * 5 : SEIS=FAUX; 6 : 'SI' ('EGA' MOT1 'SEISME'); 7 : 'ARGU' DIRECT*MOT; 8 : SEIS=VRAI; 9 : 'SI' ('EGA' DIRECT 'UX');QDIR='QX';EDIR = 1;'FINSI'; 10 : 'SI' ('EGA' DIRECT 'UY');QDIR='QY';EDIR = 2;'FINSI'; 11 : 'SI' ('EGA' DIRECT 'UZ');QDIR='QZ';EDIR = 3;'FINSI'; 12 : 'SINON'; 13 : 'ARGU' POINT1*POINT; 14 : 'FINSI'; 15 : * 16 : 'ARGU' MOT2*MOT POINT2*POINT MOSOR*MOT LFREQ*LISTREEL MOCHO*MOT ; 17 : *********************************************************************** 18 : * CALCUL D UNE FONCTION DE TRANSFERT * 19 : * * 20 : * MOD : LES MODES DE LA STRUCTURE * 21 : * LAMOR : LES AMORTISSEMENTS MODAUX (ENTRE 0 ET 1) * 22 : * MOT1 : DIRECTION DE LA SOLLICITATION (UX, UY, UZ ...) * 23 : * OU 'SEISME' POUR ACCELERATION D'ENSEMBLE * 24 : * POINT1 : POINT D APPLICATION DE LA SOLLICITATION * 25 : * OU DIRECTION DU SEISME * 26 : * MOT2 : DIRECTION DE LA REPONSE * 27 : * POINT2 : POINT OU ON CHERCHE LA REPONSE * 28 : * MOSOR : NATURE DES SORTIES (DEPL, VITE OU ACCE) * 29 : * LFREQ : LES FREQUENCES DE CALCUL * 30 : * MOCHO : TYPE DES SORTIES (MOPH OU REIM) * 31 : * * 32 : * EVENTUELLEMENT * 33 : * MOCOU : COULEUR DES COURBES * 34 : * * 35 : * EN SORTIE * 36 : * TRA : OBJET EVOLUTION COMPLEXE * 37 : *********************************************************************** 38 : 'REPETER' BEX 1; 39 : * 40 : 'SI' ('EXIS' TA) ; 41 : TAB_LOG = VRAI; 42 : 'SINON'; 43 : TAB_LOG = FAUX; 44 : 'SI' ('NON' ('EXIS' MOD)) ; 45 : 'MESS' '***** ERREUR *****'; 46 : 'MESS' 'IL MANQUE LA DONNEE DES MODES PROPRES'; 47 : 'QUITTER' BEX; 48 : 'FINSI'; 49 : 'FINSI'; 50 : * 51 : ISOR=0;ISO2=0; 52 : 'SI' SEIS; 53 : 'SI' ('EGA' MOSOR 'ACCE');ISO2=1;'FINSI'; 54 : 'SI' ('EGA' MOSOR 'VITE');ISO2=2;'FINSI'; 55 : 'SI' ('EGA' MOSOR 'DEPL');ISO2=3;'FINSI'; 56 : 'SINON'; 57 : 'SI' ('EGA' MOSOR 'ACCE');ISOR=1;'FINSI'; 58 : 'SI' ('EGA' MOSOR 'VITE');ISOR=2;'FINSI'; 59 : 'SI' ('EGA' MOSOR 'DEPL');ISOR=3;'FINSI'; 60 : 'FINSI'; 61 : 'SI' ((ISOR 'EGA' 0) 'ET' (ISO2 'EGA' 0)) ; 62 : 'MESS' '***** ERREUR *****'; 63 : 'MESS' 'IL FAUT DONNER DEPL, VITE OU ACCE'; 64 : 'QUITTER' BEX; 65 : 'FINSI'; 66 : * 67 : ICHO=0; 68 : 'SI' ('EGA' MOCHO 'MOPH');ICHO=1;'FINSI'; 69 : 'SI' ('EGA' MOCHO 'REIM');ICHO=2;'FINSI'; 70 : 'SI' (ICHO 'EGA' 0); 71 : 'MESS' '***** ERREUR *****'; 72 : 'MESS' 'IL FAUT DONNER MOPH OU REIM'; 73 : 'QUITTER' BEX; 74 : 'FINSI'; 75 : * 76 : 'SI' TAB_LOG; 77 : NMOD = ('DIME' (TA . 'MODES')) - 2; 78 : 'SINON'; 79 : NMOD='DIME' MOD; 80 : 'FINSI'; 81 : NAMO='DIME' LAMOR; 82 : 'SI' (NMOD '>' NAMO); 83 : 'MESS' '***** ERREUR *****'; 84 : 'MESS' 'IL MANQUE UNE OU PLUSIEURS VALEURS D AMORTISSEMENT MODAL'; 85 : 'QUITTER' BEX; 86 : 'FINSI'; 87 : * 88 : DEUPI=6.28318530718; 89 : Q4PI2=39.4784176044; 90 : FRE2=LFREQ**2; 91 : I=0; 92 : * 93 : 'REPETER' BMOD NMOD; 94 : I=I+1; 95 : 'SI' TAB_LOG; 96 : XI = TA . 'MODES' . I . 'DEFORMEE_MODALE'; 97 : FI = TA . 'MODES' . I . 'FREQUENCE'; 98 : MI = TA . 'MODES' . I . 'MASSE_GENERALISEE'; 99 : 'SINON'; 100 : XI='TIRE' MOD 'DEPL' 'RANG' I; 101 : FI='TIRE' MOD 'FREQ' 'RANG' I; 102 : MI='TIRE' MOD 'MGEN' 'RANG' I; 103 : 'FINSI'; 104 : EI='EXTR' LAMOR I; 105 : FI2=FI*FI; 106 : Q4E2F2=4*(FI2*(EI**2)); 107 : M2EIFI=-2*(EI*FI); 108 : 'SI' SEIS; 109 : 'SI' TAB_LOG; 110 : X1=-1. * (TA.'MODES'. I . 'DEPLACEMENTS_GENERALISES' . EDIR); 111 : 'SINON'; 112 : X1=-1. * ('TIRE' MOD QDIR 'RANG' I); 113 : 'FINSI'; 114 : 'SINON'; 115 : X1='EXTR' XI MOT1 POINT1; 116 : 'FINSI'; 117 : X2='EXTR' XI MOT2 POINT2; 118 : C2=X1*X2; 119 : COEF=C2/MI; 120 : NUMRE='PROG' 'LINE' 'A' -1 'B' FI2 FRE2; 121 : NUMIM='PROG' 'LINE' 'A' M2EIFI LFREQ; 122 : DEN1='PROG' 'LINE' 'A' Q4E2F2 FRE2; 123 : DEN2=NUMRE**2; 124 : DENOM=DEN1+DEN2;'DETR' DEN1;'DETR' DEN2; 125 : UNSDE=DENOM**-1;'DETR' DENOM; 126 : EV1=NUMRE*UNSDE;EV2=NUMIM*UNSDE; 127 : 'DETR' NUMRE;'DETR' NUMIM;'DETR' UNSDE; 128 : RES1=COEF*EV1;'DETR' EV1; 129 : IMS1=COEF*EV2;'DETR' EV2; 130 : * 131 : 'SI' (I 'EGA' 1); 132 : RES='COPI' RES1;'DETR' RES1; 133 : IMS='COPI' IMS1;'DETR' IMS1; 134 : 'SINON'; 135 : RES2=RES+RES1;'DETR' RES;'DETR' RES1; 136 : RES='COPI' RES2;'DETR' RES2; 137 : IMS2=IMS+IMS1;'DETR' IMS;'DETR' IMS1; 138 : IMS='COPI' IMS2;'DETR' IMS2; 139 : 'FINSI'; 140 : * 141 : 'FIN' BMOD; 142 : * 143 : 'SI' (ISOR 'EGA' 1); 144 : PRES=RES*FRE2;'DETR' RES; 145 : PIMS=IMS*FRE2;'DETR' IMS; 146 : RES=-1*PRES;'DETR' PRES; 147 : IMS=-1*PIMS;'DETR' PIMS; 148 : 'FINSI'; 149 : 'SI' (ISOR 'EGA' 2); 150 : PRES=IMS*LFREQ;'DETR' IMS; 151 : PIMS=RES*LFREQ;'DETR' RES; 152 : RES=PRES/(-1*DEUPI);'DETR' PRES; 153 : IMS=PIMS/DEUPI;'DETR' PIMS; 154 : 'FINSI'; 155 : 'SI' (ISOR 'EGA' 3); 156 : RES=RES/Q4PI2; 157 : IMS=IMS/Q4PI2; 158 : 'FINSI'; 159 : * 160 : 'SI' SEIS; 161 : 'SI' (ISO2 'EGA' 1); 162 : PRES=RES*FRE2;'DETR' RES; 163 : PIMS=IMS*FRE2;'DETR' IMS; 164 : RES=-1*PRES;'DETR' PRES; 165 : IMS=-1*PIMS;'DETR' PIMS; 166 : SI (EGA DIRECT MOT2); 167 : RES=RES+('PROG' ('DIME' RES) * 1.); 168 : FINSI; 169 : 'FINSI'; 170 : 'SI' (ISO2 'EGA' 2); 171 : PRES=(IMS*LFREQ)/DEUPI;'DETR' IMS; 172 : PIMS=(RES*LFREQ)/DEUPI;'DETR' RES; 173 : RES=-1*PRES;'DETR' PRES; 174 : IMS= 1*PIMS;'DETR' PIMS; 175 : 'FINSI'; 176 : 'SI' (ISO2 'EGA' 3); 177 : PRES=RES/Q4PI2;'DETR' RES; 178 : PIMS=IMS/Q4PI2;'DETR' IMS; 179 : RES=1*PRES;'DETR' PRES; 180 : IMS=1*PIMS;'DETR' PIMS; 181 : 'FINSI'; 182 : 'FINSI'; 183 : * 184 : 'DETR' FRE2; 185 : * 186 : 'ARGUMENT' MOCOU/MOT; 187 : 'SI' ('NON' ('EXISTE' MOCOU)); 188 : MOCOU= VALEUR 'COUL'; 189 : 'FINSI'; 190 : * 191 : 'SI' (ICHO 'EGA' 1); 192 : TRA='EVOL' MOCOU 'COMP' 'FREQUENCE' LFREQ 'MODULE' RES 'PHASE' IMS; 193 : 'RIMP' TRA; 194 : 'SINON'; 195 : TRA='EVOL' MOCOU 'COMP' 'FREQUENCE' LFREQ 'P REELLE' RES 196 : 'P IMAGINAIRE' IMS; 197 : 'FINSI'; 198 : * 199 : 'FIN' BEX; 200 : * 201 : 'FINPROC' TRA; 202 : 203 : 204 :
© Cast3M 2003 - All rights reserved.
Disclaimer