Procedures | Cast3M

   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 : *------------------------------------------------------------