Télécharger transfer.procedur

Retour à la liste

Numérotation des lignes :

  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.  
  205.  

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