Télécharger la plaquette Cast3M


fr en

Télécharger annoimp.procedur

Retour à la liste

Numérotation des lignes : 

  1. * ANNOIMP   PROCEDUR  NOBODY    92/01/01    00:00:00     001                    
  2. ************************************************************                    
  3. *                                                                               
  4. *--------------------------------------------------------------                 
  5. *         PROCEDURE ANNOIMP                                                     
  6. *                                                                               
  7. *    GENERATION D UN MAILLAGE D ANNEAU IMPARFAIT D ANGLE ANGM                   
  8. *                                                                               
  9. *     EQUATION R = RAY + DO COS(MOD *TETA)                                      
  10. *--------------------------------------------------------------                 
  11. *    ENTREE  DANS L'ORDRE                                                       
  12. *                                                                               
  13. *      NELEM    NOMBRE D ELEMNTS                                                
  14. *      MOD      MODE DU DEFAUT                                                  
  15. *      RAY      RAYON                                                           
  16. *      H        HAUTEUR                                                         
  17. *      DO       DEFAUT                                                          
  18. *      ANGM     ANGLE DU MAILLAGE                                               
  19. *                                                                               
  20. *    SORTIE                                                                     
  21. *                                                                               
  22. *      CYL     MAILLAGE DU CYLINDRE                                             
  23. *            CAS TRIDIM COTE2 FACE SUPERIEURE                                   
  24. *                       LIBAS FACE INFERIEURE                                   
  25. *                       COTE4 FACE A ANG =0                                     
  26. *                       COTE3 FACE A ANG =ANGM                                  
  27. *                                                                               
  28. *--------------------------------------------------------------                 
  29.                DEBPROC ANNOIMP                                                  
  30. NELEM*ENTIER   MOD*ENTIER   RAY*FLOTTANT   H*FLOTTANT                           
  31. DO*FLOTTANT  ANGM*FLOTTANT ;                                                    
  32. *--------------------------------------------------------------                 
  33. *                                                                               
  34. ************* MAILLAGE   QUART DE CYLINDRE AVEC DEFAUT D0 RAYON R               
  35. *             NOMBRE D ELEMTS NELEM                                             
  36. *                                                                               
  37.       IDIM=  VALE DIME;                                                         
  38.    XELEM=NELEM * 1.;                                                            
  39.     XMODE = MOD ;                                                               
  40.  DANG=ANGM ;DANG= DANG *( (XELEM ) ** (-1.));                                   
  41.   DANGM=XMODE * DANG;                                                           
  42. SI ( IDIM EGA 3 );                                                              
  43. PA= (RAY+DO) 0 0;PB=(RAY+DO) 0 H;                                               
  44.       FINSI;                                                                    
  45. SI ( IDIM EGA 2 );                                                              
  46. PA= (RAY+DO) 0 ;PB=(RAY+DO) 0 ;                                                 
  47.       FINSI;                                                                    
  48. P0=PA;  I=1;                                                                    
  49.  CSA = COS ( I * DANGM) ;                                                       
  50.       XI=RAY +(DO * CSA ); ANGI = DANG * I;  CSB= COS(ANGI);                    
  51.       SSB = SIN(ANGI);                                                          
  52. P1 = (XI * CSB ) (XI * SSB ) 0.;                                                
  53. LIBAS= DROIT 1 P0 P1;                                                           
  54.       NELM1=NELEM - 1;                                                          
  55.       REPETER  BOUCL1 NELM1 ;                                                   
  56. I = I + 1;                                                                      
  57. P0=P1;                                                                          
  58.  CSA = COS ( I * DANGM) ;                                                       
  59.       XI=RAY +(DO * CSA) ; ANGI = DANG * I;  CSB= COS(ANGI);                    
  60.       SSB = SIN(ANGI);                                                          
  61. SI ( IDIM EGA 3 );                                                              
  62. P1 = (XI * CSB ) (XI * SSB ) 0.;                                                
  63.       FINSI;                                                                    
  64. SI ( IDIM EGA 2 );                                                              
  65. P1 = (XI * CSB ) (XI * SSB ) ;                                                  
  66.       FINSI;                                                                    
  67. LIBAS= P0 D 1 P1 ET LIBAS ;                                                     
  68. FIN BOUCL1;                                                                     
  69. SI ( IDIM EGA 3)    ;                                                           
  70. CYL = LIBAS TRANS 1 (0 0 H)  ;                                                  
  71.       O1 = 1000*RAY ;                                                           
  72.       OEIL = O1 O1 O1;                                                          
  73.        TRAC OEIL CYL;                                                           
  74. FINSI;                                                                          
  75.  
  76. SI ( IDIM EGA 2)    ;                                                           
  77. CYL = LIBAS;                                                                    
  78.   TRAC      CYL;                                                                
  79.  
  80. FINSI;                                                                          
  81.  
  82. *------------------------------------------------------------                   
  83.                FINPROC                                                          
  84.   CYL   ;                                                                       
  85. *------------------------------------------------------------                   
  86.  

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