Télécharger la plaquette Cast3M


fr en

Télécharger ddfour.procedur

Retour à la liste

Numérotation des lignes : 

  1. * DDFOUR    PROCEDUR  TTMF3     96/06/04    21:15:12     2157           
  2.  'DEBPROC'  DDFOUR GEO*'MAILLAGE' CHARM*'MAILLAGE'  NHARM*'ENTIER'              
  3.   SOL*'CHPOINT'  RHARM*'FLOTTANT'  ORIG/'POINT' LIS*'LOGIQUE' 
  4. ITYPE/'ENTIER' ;
  5. *********************************************************************   
  6. *        ANALYSE HARMONIQUE  DU POTENTIEL VECTEUR                       
  7. *        GEO    MAILLAGE SUPPORT SOLUTION GENERALE                  *   
  8. *        CHARM  ARC DE CERCLE SUR LEQUEL ON A LE POTENTIEL          *   
  9. *        NHARM  NOMBRE D HARMONIQUES A CALCULER                     *   
  10. *        SOL  SOLUTION  EN POTENTIEL                                *   
  11. *        RHARM  RAYON DE NORMALISATION                              *   
  12. *        ORIG   ORIGINE CERCLE D ANALYSE                            *    
  13. *        LIS LOGIQUE  FAUX SI PAS DE LISSAGE POLYNOMIAL             * 
  14. *        ITYPE  ORDRE DE SYMETRIE  POUR L EXPANSION POLY SI LIS=VRAI*  
  15. *********************************************************************   
  16.  1P = CHARM 'POINT' 'INITIAL' ;2P= CHARM 'POINT' 'FINAL' ;                      
  17.         'SI' ( 'EXISTE' ORIG ) ;                                            
  18.         RCIRC = 'NORM' ( 1P 'MOINS' ORIG ) ;                                
  19.         'SINON' ;                                                         
  20.         RCIRC = 'NORM' 1P ;ORIG = 0. 0. ;                                 
  21.         'FINS' ;                                                         
  22.  X1 Y1 = 'COOR' ( 1P 'MOINS' ORIG );X2 Y2 = 'COOR' 2P ;                       
  23.  PT = 2P 'PROJETER' (( Y1 * -1. )  X1 )  'DROITE' ORIG 1P ;                 
  24.  H = 'NORM' ( 2P 'MOINS' PT ) ;                                             
  25.      'SI' ( (X1 + X2 ) < 1.E-3 ) ;                                        
  26.  NUM = 2 ;                                                              
  27.  'MESS'   'SOLUTION DONNEE SUR  180  DEGRES ' ;                           
  28.      'SINON' ;                                                            
  29.  ANG = 'ATG' H ( 'NORM' ( PT 'MOINS' ORIG)) ;                                 
  30.  NUM = ( 'ENTI' ( 360.1 / ANG )) ;                                        
  31.  'MESS'   'SOLUTION DONNEE SUR ' ANG ' DEGRES '   ;                       
  32.      'FINSI' ;                                                            
  33.    'SI' ( 'NON' ( LIS )) ;                                                  
  34.  SCIRC = INT_COMP  GEO   SOL  CHARM    ;                                
  35.  SCIRC = SCIRC 'NOMC' 'T' ;                                               
  36.    'SINON' ;                                                              
  37. * lissage polynomial                                                    
  38. * CHLI  = PROI POLY CHARM GEO   SOL    ITYPE  'PLAN' ; 
  39.  CHLI  = 'PROI' POLY GEO CHARM   SOL    ITYPE  'PLAN' ; 
  40.  SCIRC = ('EXCO' CHLI  'A' ) 'NOMC' 'T' ;                                   
  41.    'FINS' ;                                                              
  42.  RRAP = RHARM / RCIRC ;                                                 
  43.  CQTT RTOT = F_S2PI CHARM  SCIRC  NUM ORIG  ;                           
  44.   EQTT = 'EVOL' 'ROUG' 'CHPO' CQTT  'T' RTOT ;                                
  45.   AAAA=  'EXTR' EQTT 'ABSC' ;                                             
  46. * QTOT=  'EXTR' EQTT 'ORDO' ;                                             
  47. * LIST QTOT ;                                                           
  48. * NI1= (DIME AAAA ) - 1  ;                                              
  49. * IVAL = VALE ELEM ;                                                    
  50. * OPTION ELEM SEG2 ;I1 = 0. 0. ; I2 = 360. 0. ;                         
  51. * RTOT = DROI  NI1  I1 I2  ;  
  52.  OBA_RMA = 'MODEL'  RTOT  'MECANIQUE'  'COQ2';                                  
  53. * OBA_RMA= 'AFFECTE' RTOT ('MODELE' 'STANDARD') COQ2;                          
  54.  PERIOD= 'MAXI' AAAA ;                                                    
  55.  MULT= 360. / PERIOD ;                                                  
  56.  ANGVEC= AAAA * MULT  ;                                                 
  57. * MESS ' PERIODE '  PERIOD ;                                            
  58.  N=0;                                                                   
  59.  'MESS'   '  ANALYSE CONFORME A POISSON  HARMON ';                        
  60.  'MESS'   ' HARMONIQUE     POTENTIEL    CHAMP ';                          
  61.  'MESS' ' NOMBRE DE SEGMENTS DE DR ' NI1 ;                                
  62.  'MESS' ' DIME DE ANGVEC ' ( DIME ANGVEC ) ;                              
  63.  CHPA= 'MANU' 'CHPO' RTOT 1 'SCAL' ANGVEC ;                                   
  64. * CQTT = MANU CHPO RTOT 1 'POT' QTOT ;                                  
  65.  'REPETER' BLOCIT NHARM ;                                                 
  66.  N= N + 1 ;                                                             
  67.  RAN= RRAP ** N ;                                                       
  68.  ANGLEN= CHPA  *  N ;                                                   
  69.  CNX=  'COS'  ANGLEN ;                                                    
  70.  SNX=  'SIN'  ANGLEN ;                                                    
  71.  EVREL= CQTT * CNX ;EVIMA= CQTT * SNX ;     
  72.   CEVREL =EVREL 'CHAN' 'CHAM' OBA_RMA 'RIGIDITE' ;                              
  73.    CEVIMA = EVIMA 'CHAN' 'CHAM' OBA_RMA 'RIGIDITE' ;                            
  74.  AK= (( 'INTG' OBA_RMA CEVREL ) / PERIOD ) * 2. * RAN ;
  75.  BK= (( 'INTG' OBA_RMA CEVIMA ) / PERIOD ) * 2. * RAN ;
  76.  CK= (( AK * AK ) + ( BK * BK ) ) ** .5 ;                               
  77.  FF= ( 'FLOT' N) / RHARM ;                                                
  78.  AKK= AK * FF ;                                                         
  79.  BKK= BK * FF ;                                                         
  80.  CKK= CK * FF ;                                                         
  81.  'SI' (N 'EGA' 1 );AKN= 'PROG' AKK ; 'SINON' ;
  82. AKN = AKN 'ET' ('PROG' AKK ) ; 'FINSI' ;
  83.  'SAUTER' 1 'LIGNE' ;                                                       
  84. * MESS   N AK BK CK AKK BKK CKK ;                                       
  85.  'MESS'   N   AK  BK  CK  ;                                               
  86.  'MESS'   N  AKK BKK CKK  ;                                               
  87.  'FIN' BLOCIT ;                                                           
  88. * OPTION ELEM IVAL ;                                                    
  89.  'FINPROC' AKN  ANGVEC ; 
  90.  
  91.  

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