Télécharger la plaquette Cast3M


fr en

Télécharger flambage.procedur

Retour à la liste

Numérotation des lignes : 

  1. * FLAMBAGE  PROCEDUR  BP208322  17/10/03    21:15:27     9580           
  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. XX = 2. * PI ;  
  123. *       
  124. *--------------------------------------------   
  125. *     CALCUL DES FREQUENCES A PARTIR DES COEFFICIENTS MULTIPLICATEURS   
  126. *--------------------------------------------   
  127. XX0 = X0 ** .5 ; XX1 = X1 ** .5;
  128. XX0= XX0 / XX ; XX1 = XX1 / XX ;
  129. *       
  130. *---------------------------------------------  
  131. *       BOUCLE SUR LES MODES DE FOURIER 
  132. *---------------------------------------------  
  133. *       
  134. TAMON = 'TABLE';
  135. TAMOM = 'TABLE';
  136. TACOE = 'TABLE';
  137. TAPERM = 'TABLE';       
  138. STAB11 = 'TABLE' ;      
  139.       II = 0 ;  
  140.       IK = 0 ;  
  141.       MODMIN= -10;      
  142.    COEMIN=1.E30;
  143. *       
  144. *  STAB11 EST UNE TABLE DE TABLES INDICEE PAR LE NUMERO DU MODE SORTI   
  145. *       
  146. IFOU=0  ;       
  147. 'REPETER' BOUCFOU;      
  148. IFOU = IFOU + 1 ;       
  149. 'SI' ( IFOU >   NBMOD) ;
  150. 'QUITTER' BOUCFOU;      
  151. 'FINSI' ;       
  152.       N = 0;    
  153. 'SI' (IFOURI 'EGA' 1) ; 
  154. N = 'EXTR' LMOD IFOU;   
  155. 'OPTION' 'MODE' 'FOUR' N ;      
  156. 'FINSI';
  157. 'SI' (IRIGN 'EGA'  0) ; 'SI' (ICARA 'EGA' 1) ;MATN1=MATN 'ET' CAR1;     
  158. 'SINON'  ;MATN1=MATN  ;'FINSI' ; RIGN = 'RIGI' OBJMOD MATN1;'FINSI';    
  159. TST = 'EXISTE' ETAB11 'SIG0';   
  160. 'SI' ( TST ) ;  
  161. SIG10 = ETAB11.'SIG0'; 
  162. KSIG0 = 'KSIGMA' OBJMOD SIG10 CAR1 FLAM ;       
  163. RIGN = RIGN 'ET' KSIG0 ;
  164. 'FINSI' ;       
  165. 'SI' (IKSIG 'EGA' 0 ) ; KSIG1 = 'KSIG' OBJMOD SIG11  CAR1 FLAM ;
  166.   'FINSI';      
  167. 'SI' ( 'EXISTE' ETAB11 'CLIM' ) ; RIGP = RIGN 'ET' KLIM;
  168.       'SINON';RIGP = RIGN ;     
  169.  'FINSI';       
  170. *------------------------------------------------------------   
  171. *       CALCUL DES FREQUENCES   
  172. *-----------------------------------------------------------    
  173. TMOD = 'VIBR' 'INTERVALLE' XX0 XX1 'BASSE' NMO RIGP KSIG1 ; 
  174. XMOD = TMOD . 'MODES';
  175. *-----------------------------------------------------------    
  176. *       COMBIEN A T ON DE MODES 
  177. *-------------------------------------------------------
  178.    NMOCAL  = 'DIMENSION' XMOD - 2;  
  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 =  XMOD . JJ . FREQUENCE ;       
  200. DD =  XMOD . JJ . DEFORMEE_MODALE ;       
  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.  
  278.  
  279.  
  280.  

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